Keywords:

Global food security is endangered by fungal phytopathogens causing devastating crop production losses. Many of these pathogens use specialized appressoria cells to puncture plant cuticles. Here, we unveil a pair of alcohol oxidase-peroxidase enzymes to be essential for pathogenicity. Using Colletotrichum orbiculare, we show that the enzyme pair is co-secreted by the fungus early during plant penetration, and that single and double mutants have impaired penetration ability. Molecular modeling, biochemical and biophysical approaches revealed a fine-tuned interplay between these metalloenzymes, which oxidize plant cuticular long-chain alcohols into aldehydes. We show that the enzyme pair is involved in transcriptional regulation of genes necessary for host penetration. The identification of these infection-specific metalloenzymes opens new avenues on the role of wax-derived compounds and the design of oxidase-specific inhibitors for crop protection.

Tandem metalloenzymes gate plant cell entry by pathogenic fungi

Bastien Bissaro 1 † , Sayo Kodama 2 † , Takumi Nishiuchi 3 , Anna Maria Díaz-Rovira 4 , Hayat Hage 1 , David Ribeaucourt 1,6,7 , Mireille Haon 1 , Sacha Grisel 1 , A. Jalila Simaan 6 , Fred Beisson 8 , Stephanie M. Forget 9 , Harry Brumer 9 , Marie-Noëlle Rosso 1 , Victor Guallar 4,5 , Richard O´Connell 10 , Mickaël Lafond 6 , Yasuyuki Kubo 2* and Jean-Guy Berrin 1*

Introduction

Fungal phytopathogens represent a serious threat to plant health [START_REF] Kamoun | Plant health emergencies demand open science: Tackling a cereal killer on the run[END_REF] and global food security [START_REF] Fisher | Worldwide emergence of resistance to antifungal drugs challenges human health and food security[END_REF]. Colletotrichum and Magnaporthe species rank among the top 10 most devastating fungal phytopathogens in the world and reduce crop yield by up to 30% [START_REF] Dean | The Top 10 fungal pathogens in molecular plant pathology[END_REF]. Despite being separated by ca. 300 million years of evolution [START_REF] Beimforde | Estimating the Phanerozoic history of the Ascomycota lineages: Combining fossil and molecular data[END_REF][START_REF] Liang | Pathogenic adaptations of Colletotrichum fungi revealed by genome wide gene family evolutionary analyses[END_REF], these fungi share remarkable similarities in their infection strategy, notably the formation of a specialized cell dedicated to host penetration called an appressorium (6-8) (Fig. 1A). This dome-shaped, darkly melanized cell generates a high internal turgor and directs this mechanical pressure onto a needle-like penetration peg [START_REF] Chen | Appressorium Turgor Pressure of Colletotrichum kahawae Might Have a Role in Coffee Cuticle Penetration[END_REF][START_REF] Ryder | A sensor kinase controls turgor-driven plant infection by the rice blast fungus[END_REF][START_REF] He | Discovery of broad-spectrum fungicides that block septin-dependent infection processes of pathogenic fungi[END_REF], which emerges from a 2-500 nm pore at the appressorial base to puncture the plant's outer defensive barriers, namely the cuticle and epidermal cell wall.

Despite major advances in our understanding of the cellular processes preceding [START_REF] Ryder | Regulation of appressorium development in pathogenic fungi[END_REF][START_REF] De Jong | Glycerol generates turgor in rice blast[END_REF][START_REF] He | Discovery of broad-spectrum fungicides that block septin-dependent infection processes of pathogenic fungi[END_REF][START_REF] Dagdas | Septin-mediated plant cell invasion by the rice blast fungus, Magnaporthe oryzae[END_REF][START_REF] Kodama | The morphogenesisrelated NDR kinase pathway of Colletotrichum orbiculare is required for translating plant surface signals into infection-related morphogenesis and pathogenesis[END_REF][START_REF] Rocha | Spermine-mediated tight sealing of the Magnaporthe oryzae appressorial pore-rice leaf surface interface[END_REF] and following [START_REF] Giraldo | Filamentous plant pathogen effectors in action[END_REF][START_REF] Djamei | Metabolic priming by a secreted fungal effector[END_REF] plant cell entry, the (bio)chemical reactions occurring at this nanoscale plant-fungus interface and their role in host penetration are not fully elucidated.

Copper radical oxidases (CROs) are enzymes with diverse substrate specificities, and have been extensively studied since the late 1960s [START_REF] Cooper | Galactose Oxidase from Polyporus circinatus[END_REF]. Today, CROs include galactose 6oxidases (GalOx) [START_REF] Cooper | Galactose Oxidase from Polyporus circinatus[END_REF], glyoxal oxidases [START_REF] Kersten | Involvement of a new enzyme, glyoxal oxidase, in extracellular H2O2 production by Phanerochaete chrysosporium[END_REF], and broad specificity primary alcohol oxidases (AlcOx) [START_REF] Tyler) Yin | Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family[END_REF][START_REF] Mathieu | Discovery of a Fungal Copper Radical Oxidase with High Catalytic Efficiency toward 5-Hydroxymethylfurfural and Benzyl Alcohols for Bioprocessing[END_REF]. Despite the detailed knowledge available on the enzymology and structure of CROs, their biological function remains largely unknown. It has only been proposed that glyoxal oxidases play a role in lignin degradation by fungal saprotrophs [START_REF] Kersten | Involvement of a new enzyme, glyoxal oxidase, in extracellular H2O2 production by Phanerochaete chrysosporium[END_REF]. During the course of our previous work [START_REF] Tyler) Yin | Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family[END_REF], we noted that genes encoding secreted AlcOx orthologs are particularly widespread among phytopathogenic ascomycete fungi and are absent in plants.

Given that long-chain primary alcohols are components of the waxy cuticle of aerial plant surfaces [START_REF] Lee | Advances in the understanding of cuticular waxes in Arabidopsis thaliana and crop species[END_REF], we hypothesized that AlcOx could play a role in fungal pathogenesis.

In this study, we used a combination of 'omics analyses to unveil the pairing of AlcOx with a redox partner, namely heme-peroxidase, in some fungal plant pathogens of high agricultural importance. We used wet enzymology and molecular modeling to demonstrate and characterize the interplay between these metalloenzymes. Reverse genetics, live-cell imaging and fungal transcriptomics allowed us to probe the in vivo function of the enzyme pair, providing new molecular insights into the host penetration cascade.

Results

Discovery of the tandem Perox-AlcOx

While studying the enzymology of Colletotrichum AlcOx enzymes for biotechnological applications [START_REF] Tyler) Yin | Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family[END_REF][START_REF] Ribeaucourt | Controlled Oxidation of Long-Chain Primary Aliphatic Alcohols by a Fungal Copper-Radical Alcohol Oxidase[END_REF], we noticed the presence of a gene encoding a putative peroxidase located adjacent to an AlcOx-encoding gene (Fig. 1B). To strengthen this initial observation, we searched for alcox orthologs in 30 sequenced Colletotrichum genomes, which revealed the near-ubiquitous presence of a putative peroxidase (hereafter called "Tandem Peroxidase"). Interestingly, the perox and alcox genes were found in a head-to-head arrangement (fig. S1), suggesting the presence of a bidirectional promoter for tight coexpression of the genes. Interestingly, these Tandem Peroxidases are never found in combination with other types of CROs (fig. S1C) and both proteins encoded by the perox-alcox pair are predicted to be secreted (Table S1). These observations aroused our interest because it is known that CROs require activation by horseradish peroxidase (HRP) for maximum activity in vitro [START_REF] Parikka | Oxidation of methyl α-D-galactopyranoside by galactose oxidase: products formed and optimization of reaction conditions for production of aldehyde[END_REF][START_REF] Forget | Determination of biocatalytic parameters of a copper radical oxidase using real-Time reaction progress monitoring[END_REF].

Our phylogenetic analysis of the peroxidase-catalase superfamily showed that the Tandem Peroxidases cluster together in a sister clade within the under-explored ascomycete Class II peroxidases [START_REF] Zerva | Characterization and application of a novel class II thermophilic peroxidase from Myceliophthora thermophila in biosynthesis of polycatechol[END_REF] (fig. S2A). Furthermore, Tandem Peroxidases form a distinct clade amongst the 333 Class II peroxidases found in Colletotrichum species (fig. S2B), suggestive of neofunctionalization. A broad search for the co-occurrence of perox and alcox orthologs across fungal genomes revealed that the pair is also present in Magnaporthe species, including the infamous causal agent of rice blast, Magnaporthe oryzae (syn. Pyricularia oryzae) (Table S2). Mapping the occurrence of Perox-AlcOx protein pairs and their corresponding genomic neighborhoods onto a phylogeny of representative pathogenic ascomycetes (Fig. 1B, fig. S3) allowed us to conclude that the pair is present in most Colletotrichum species complexes for which genome sequences are available and in Magnaporthe spp., and that the head-to-head organization of the pair is conserved in all these fungi, suggesting there is selection pressure to retain the pairing and that it is has a critical role in the biology of these pathogens (Fig. 1B).

To further test the hypothesis of a functional linkage of the perox and alcox gene products, we parsed transcriptomic data available for those fungal species harboring the pair.

Remarkably, this analysis revealed that both genes are always tightly co-transcribed at the appressorium stage in various pathosystems involving Colletotrichum species attacking maize, fruits, and model plants [START_REF] O'connell | Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses[END_REF][START_REF] Gan | Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi[END_REF], and in M. oryzae attacking rice and barley (28, 29) (Fig. 1C andD). In each case, the transcript levels are relatively low and detected within a narrow time window, which may explain why these genes were overlooked in previous studies. S1) during plant infection for different pathosystems (Myc., Mycelium) [START_REF] O'connell | Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses[END_REF][START_REF] Gan | Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi[END_REF][START_REF] Dong | Global Genome and Transcriptome Analyses of Magnaporthe oryzae Epidemic Isolate 98-06 Uncover Novel Effectors and Pathogenicity-Related Genes, Revealing Gene Gain and Lose Dynamics in Genome Evolution[END_REF][START_REF] Shimizu | RNA-Seq of in planta-expressed Magnaporthe oryzae genes identifies MoSVP as a highly expressed gene required for pathogenicity at the initial stage of infection[END_REF].

Actual time points associated with each infection stages are provided in the Material and Methods section. (D) Illustration of targeted plant hosts.

Tandem Perox-AlcOx oxidize plant long chain alcohols

As a prelude to analyzing their biological function in vivo, we studied the substrate specificity and enzyme interplay of the Perox-AlcOx pair in vitro (Fig. 2A). C. orbiculare was selected as a model because it not only causes the economically important anthracnose disease of cucurbits (e.g., melons, cucumber) but also has been used for decades as a model system for studying fungal pathogenesis [START_REF] Kubo | Dynamics of infection-related morphogenesis and pathogenesis in Colletotrichum orbiculare[END_REF]. Despite the notorious difficulties associated with heterologous expression of such metalloenzymes, we successfully produced in the yeast Pichia pastoris recombinant copper-radical AlcOx and heme-iron tandem peroxidase from C. orbiculare (hereafter CorAlcOx and CorPerox, respectively). Similar to the previously studied AlcOx orthologs from C. graminicola and C. gloeosporioides [START_REF] Tyler) Yin | Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family[END_REF] (fig. S1A), CorAlcOx oxidized both aromatic and long-chain aliphatic primary alcohols (fig. S4A). This finding raises the possibility that fatty primary alcohols present in the cuticle of many plant species [START_REF] Lee | Advances in the understanding of cuticular waxes in Arabidopsis thaliana and crop species[END_REF], including cucumber, could be the native substrates of these enzymes.

On the other hand, CorPerox was confirmed to be a peroxidase, albeit with moderate catalytic efficiency (kcat = 1.52 ± 0.02 s -1 , KM H2O2 = 80 ± 3 µM, kcat/KM = 1.9 × 10 4 ± 0.1 s -1 M - 1 ) compared to the only previously characterized ascomycete Class II peroxidase [START_REF] Zerva | Characterization and application of a novel class II thermophilic peroxidase from Myceliophthora thermophila in biosynthesis of polycatechol[END_REF], commercial horseradish peroxidase (HRP) [START_REF] Rodriguez-Lopez | Role of Arginine 38 in Horseradish Peroxidase. A critical residue for substrate binding and catalysis[END_REF], and well-studied basidiomycete lignin-active peroxidases (32) (kcat/KM = 10 4 -10 7 s -1 M -1 ). CorPerox was only active on low redox-potential substrates and not on any of the substrates of canonical lignin-active peroxidases (fig. S4B), and required the presence of calcium ions for stability (fig. S4C). These observations are in agreement with structural predictions (fig. S5), which indicate the presence of two conserved calcium ion binding sites, but absence of the manganese binding site and the surface-exposed tryptophan involved in long-range electron transfer, which are two key features of lignin-active peroxidases [START_REF] Ayuso-Fernández | Evolutionary convergence in lignin-degrading enzymes[END_REF].

Despite its comparatively low peroxidase activity, CorPerox activates CorAlcOx for oxidation of primary alcohols in a dose-dependent manner, and to a much greater extent than the plant peroxidase HRP does (Fig. 2B). The pH optima of CorAlcOx and CorPerox were markedly different (ca. 8 and 4, respectively, fig. S4D). Importantly, the pH of the environment measured on cucumber cotyledons at the time of triggering appressorium penetration by C. orbiculare was between 7.5 and 8.0. Taken together, these results indicate that AlcOx activation is not dependent on highly efficient peroxidase activity. We also heterologously produced the Perox-AlcOx pair from the rice blast pathogen M. oryzae, of which MorAlcOx was recently confirmed to be a primary alcohol oxidase [START_REF] Oide | Carbohydrate-binding property of a cell wall integrity and stress response component (WSC) domain of an alcohol oxidase from the rice blast pathogen Pyricularia oryzae[END_REF]. Here, we obtained an activation profile for the MorPerox-MorAlcOx pair resembling that observed for the C. orbiculare pair (fig. S4E cf. Fig. 2B).

Having determined optimal enzyme activation conditions, we then probed further the activity of the CorPerox-CorAlcOx pair on biologically relevant aliphatic alcohols. Despite challenges associated with substrate solubility in aqueous buffer, we clearly detected activity on hexadecan-1,16-diol and octadecan-1-ol, as well as on a crude preparation of waxes extracted from cucumber cotyledons (Fig. 2C). Importantly, such activity was detected only when both enzymes were present. Product analysis by gas chromatography unambiguously indicated that octadecan-1-ol was oxidized to the corresponding aldehyde (fig. S4F).

To obtain a deeper understanding of the activation of AlcOx by the Tandem Peroxidase, we analyzed electron transfer by electron paramagnetic resonance (EPR) spectroscopy (Fig. 2D) and state-of-the art molecular modeling (Fig. 2E and fig. S6). Reduction in the EPR signal of the inactive Cu(II)-non-radical form of CorAlcOx upon addition of CorPerox was supportive of one-electron oxidation leading to the EPR silent, active Cu(II)-radical form (Fig. 2A andD). This change in electronic structure, observed in the absence of any substrate, confirms the activity-independent activating role of CorPerox. It also indicates a close contact between the enzymes during the activation process, which is concordant with our modeling studies (see below) and with the sigmoidal, titration-like curves observed during activity assays carried out in the presence of substrate (Fig. 2B). Further, CorAlcOx-CorPerox top 5 models predicted by two independent computational techniques, PIPER [START_REF] Kozakov | PIPER: An FFT-based protein docking program with pairwise potentials[END_REF] and AlphaFold2-Multimer (36) consistently placed CorPerox structures in front of the CorAlcOx active site (fig. S6A). Refinement of these models with the all atom Monte Carlo software PELE [START_REF] Borrelli | PELE: Protein energy landscape exploration. A novel Monte Carlo based technique[END_REF], which includes protein small rotations and translations followed by an exhaustive side-chain prediction at the interface, resulted in a clear minimum (fig. S6B) where the heme group of CorPerox is oriented towards the CorAlcOx copper ion (Fig. 2E and fig. S6C). Notably, the binding surfaces of AlcOx orthologs are considerably more hydrophobic than that of GalOx (fig. S7A andB), which is consistent with the idea that AlcOx may have evolved to interact with the hydrophobic plant cuticle. Moreover, in the top model (Fig. 2E), there is a significant decrease of solvent exposure of the CorAlcOx active site, defining a cavity between both enzymes that could facilitate diffusion and binding of a long-chain alcohol substrate. C18 docking and PELE induced fit simulations confirm this point, revealing a pronounced local minimum in which the alcoholic group of the substrate is well-positioned for catalysis (fig. S6D). Lowest energy CorAlcOx-CorPerox complex obtained by protein-protein modeling simulation with PELE (see fig. S6 for more details). The copper atom is shown as a blue sphere, the heme group as magenta sticks and calcium ions as purple spheres.

The Perox-AlcOx pair gates plant penetration

To investigate the role of the Perox-AlcOx pair in plant infection, we isolated single and double gene deletion mutants of C. orbiculare (fig. S8). Inoculation of spore suspensions onto intact cucumber cotyledons showed that fewer and smaller lesions were formed by all the mutants compared to the wild-type strains (Fig. 3A andB). For instance, the proportion of lesions with a diameter > 4 mm fell from 95 % to < 20%. Furthermore, similar phenotypes were obtained for single and double mutants, suggesting that both oxidases are crucial for fungal pathogenicity.

Further investigations indicated that neither single nor double perox/alcox deletion mutants were affected in mycelial growth (fig. S8A). Microscopy revealed that the significant loss of pathogenicity of the mutants was due to a large decrease in the frequency of host penetration (Fig. 3C). However, morphogenesis, cell wall melanization, and turgor build-up within appressoria cells were virtually indistinguishable from those of the wild-type strain (fig.

S8B-D).

The next step in the infection process is the emergence of a needle-like penetration peg through a pore in the basal cell wall of the appressorium (Fig. 1A), during which actin assembly at the pore provides rigidity [START_REF] Fukada | Colletotrichum orbiculare Regulates Cell Cycle G1/SProgression via a Two-Component GAP and a GTPase to Establish Plant Infection[END_REF]. Using a red fluorescent protein/actin-binding protein fusion (Lifeact-RFP), we found normal actin assembly at the appressorium pore for both the single and double mutants inoculated onto cucumber cotyledons (fig. S8E andF).

Remarkably, the mutants could penetrate and form hyphae inside inert cellophane membranes (fig. S9A-C) and caused wild-type-like lesions when inoculated on mechanically wounded cucumber cotyledons (Fig. S9D andE), in contrast to the crippled invasive capacity observed on intact cotyledons (Fig. 3C). This points to a mechanism involving plant surface compounds, concordant with the catalytic activity of the Perox-AlcOx pair (Fig. 2C). These results collectively indicate that CorAlcOx and CorPerox play a crucial role during the early penetration stage, but are not involved in either appressorium or peg formation.

To further examine the function of the Perox-AlcOx pair during plant infection, we attempted to localize the proteins by live-cell imaging of CorAlcOx-mCherry and CorPerox-GFP driven by their native promoters. Although the CorAlcOx-mCherry and CorPerox-GFP complemented the defect in pathogenicity of the deletion mutants, fluorescence of CorAlcOx-mCherry and CorPerox-GFP was not detectable during appressorium formation on cucumber cotyledons (fig. S10A), suggesting that gene expression was too low, consistent with transcriptomic data (Fig. 1C), or that the gene products were secreted and diffused away from the penetration site. However, the constitutive overexpression of CorAlcOx-mCherry and CorPerox-GFP, driven by the translation elongation factor (TEF) promoter, revealed that both proteins accumulated specifically at the appressorial penetration pore and that signal intensity increased during penetration peg formation (Fig. 3D, fig. S10B andC, and Supplementary Movie 1). This protein co-localization observed in vivo is consistent with gene co-expression data (Fig. 1C) and rationalizes the co-operative activity demonstrated by biochemical assays (Fig. 2). Remarkably, CorAlcOx-mCherry and CorPerox-GFP were detected at the plant surface beneath detached appressoria, at the penetration site (Fig. 3E), suggesting that the tandem metalloenzymes are secreted from appressoria into the plant epidermis. Strikingly, CorAlcOx-mCherry and CorPerox-GFP were not detected at the penetration site on cellophane membranes (fig. S10D), suggesting that interaction of the fungus with the plant surface triggers local and specific recruitment of the tandem metalloenzymes to the pore.

To probe the presence of potential natural substrates of the fungal AlcOx in the neighborhood of the penetration site, we carried out a compositional analysis of waxes present at the surface of uninoculated cucumber cotyledons (fig. S11). This analysis showed that the extracted waxes are mainly composed of odd-numbered alkanes (C27-C33) and evennumbered long-chain primary alcohols (C24-C32). Fatty aldehydes were found only as traces.

Thus, this experiment demonstrates that potential substrates of AlcOx represent a major part of the plant cuticular compounds while AlcOx reaction products are very minor components.

To further explore the role of Perox-AlcOx in the fungus-plant dialogue, we exposed the C. orbiculare single and double mutants to a product of the AlcOx, viz. the aliphatic long chain aldehyde n-octadecanal. For all gene deletion mutants, the addition of n-octadecanal partially restored appressorium penetration ability and lesion formation on cucumber leaves (Fig. 4A andB), suggesting that the role of the fungal Perox-AlcOx pair is to generate long chain aldehydes to prime the fungus for efficient plant infection.

Taken together, our results suggest that the role of the fungal Perox-AlcOx pair is to increase locally the concentration of long-chain aldehydes. These aliphatic compounds, members of the volatile organic compounds, are well-known to function as signal molecules [START_REF] Pennerman | Eight-carbon volatiles: prominent fungal and plant interaction compounds[END_REF], raising the possibility that the Perox-AlcOx pair generates signals to prime the fungus for efficient plant infection.

To gain insights into the possible steps controlled by the Perox/AlcOx activities during the early infection process, we performed comparative transcriptomic analyses of C. orbiculare wild-type and double mutant strains at the appressorial stage on cucumber leaves and cellophane (Fig. 5). Analysis of the differentially-expressed genes suggested that the Perox-AlcOx pair contributes to the regulation of a subset of 32 plant-inducible genes predicted to encode small secreted proteins (SSPs), carbohydrate-active enzymes (CAZymes) and membrane transporters (Table S6). SSPs and CAZymes are well-known fungal effectors playing a key role in the molecular dialog with host plants. A phylogenetic analysis of CAZymes (Table S6) present in this subset of genes revealed that the pair is required for the up-regulation of genes encoding proteins directed towards the plant cell wall (PCW) and the fungal cell wall (FCW) (Fig. 5). Among the former, we detected cellulose-active enzymes known to display an enzymatic interplayviz. the cello-oligosaccharide dehydrogenase AA7

and the lytic polysaccharide monooxygenase (LPMO) AA9 ( 40)as well as pectin-active enzymes from the PL1, PL3 and GH93 CAZy families. Regarding the FCW-targeting proteins, we detected three proteins with carbohydrate binding modules (CBM50s; also called LysM domains) that bind to chitin and function to evade recognition by host immune receptors during infection [START_REF] Sanz-Martín | A highly conserved metalloprotease effector enhances virulence in the maize anthracnose fungus Colletotrichum graminicola[END_REF]. predicted to localize in the plant cell (in the apoplast (A) and/or the cytoplasm (C)). See Table S6 for more details.

Discussion

Research on plant invasion by appressorium-forming fungi teaches us that the development of these specialized infection structures is an extremely complex, finely regulated process. In this study, we have shown that fungal AlcOx, which are only encountered in phytopathogens, have a very different biological function than the few distantly-related CROs for which roles in morphogenesis [START_REF] Leuthner | A H2O2-producing glyoxal oxidase is required for filamentous growth and pathogenicity in Ustilago maydis[END_REF][START_REF] Chaplin | GlxA is a new structural member of the radical copper oxidase family and is required for glycan deposition at hyphal tips and morphogenesis of Streptomyces lividans[END_REF] or lignin degradation [START_REF] Kersten | Involvement of a new enzyme, glyoxal oxidase, in extracellular H2O2 production by Phanerochaete chrysosporium[END_REF] were proposed. We showed that the tandem Perox-AlcOx metalloenzymes are specifically deployed at the initiation of penetration peg formation and their localization at the pore formed at the fungusplant cell interface. Our results suggest that the oxidative action of Perox-AlcOx on long-chain alcohols triggers a yet-to-be-elucidated biochemical cascade leading to penetration. Indeed, functional complementation of the deletion mutants, which were defective in their ability to puncture intact plant cuticles, by the product of the reaction suggests that cuticular alcohol oxidation provides a chemical cue required for plant cell entry (Fig. 6). We speculate that acquisition of the Perox-AlcOx pair provided the ecologically widespread Colletotrichum and Magnaporthe species with an advantage, namely an in-house "locksmith" (i.e. the tandem metalloenzymes) ensuring localized production of the aldehyde. This entry key is needed to move to the next pathogenesis stage with the expression of PCW-degrading enzymes to facilitate penetration, and of FCW-binding proteins and carbohydrate oxidases to evade host immunity. For example, it has been suggested that, upon oxidation by fungal LPMOs and oligosaccharide oxidases, plant oligosaccharides can no longer play their role as inducers of host immune response [START_REF] Haddad Momeni | Discovery of fungal oligosaccharide-oxidising flavo-enzymes with previously unknown substrates, redox-activity profiles and interplay with LPMOs[END_REF][START_REF] Vandhana | On the expansion of biological functions of Lytic Polysaccharides Monooxygenases[END_REF][START_REF] Sabbadin | Secreted pectin monooxygenases drive plant infection by pathogenic oomycetes[END_REF].

It is interesting to speculate why Nature has evolved such a complex enzymatic system relying on heme-iron peroxidase and elaborate copper radical chemistry when simpler, monoenzyme systems could be employed for the same purpose (e.g. FAD-dependent oxidases active on long-chain alcohols [START_REF] Pickl | The substrate tolerance of alcohol oxidases[END_REF]). We show that the functional integrity of the secreted tandem metalloenzymes system endows pathogenic fungi with fine control over oxidation reactions in the extracellular space. We propose that the copper-radical center equips CROs with a "redox switch" to turn on activity with high spatio-temporal resolution, when paired with a cognate peroxidase (Fig. 6). More generally, our study suggests that fungi have evolved ways of controlling oxidative reactions that are seemingly out of reach, i.e. beyond the bounds of the fungal plasma membrane and cell wall, through tight genetic regulation and protein interplay between secreted oxidoreductases.

In conclusion, using a combination of in silico, in vitro and in vivo approaches, we have unveiled the existence of a natural redox partner, i.e. a heme-peroxidase, for AlcOx-type CROs, and showed that the pair acts as a secreted virulence factor during early infection (Fig. 6). It is noteworthy that the enzymology-driven approach pursued here was essential to bring to light this unique mechanism, because low expression levels, fine temporal tuning, and highly

Materials and Methods

Chemicals and commercial enzymes

Most chemicals were purchased from Sigma (Darmstadt, Germany) or VWR (Fontenay-sous-Bois, France) unless stated otherwise. Molar concentrations of type II HRP (Sigma-Aldrich; MW 33.89 kDa) were estimated by Bradford assay. n-octadecanal was purchased from TCI-Europe (Zwijndrecht, Belgium). All alcohol substrates stock solutions were prepared either in H2O or in acetone and stored at -20°C. The concentration of H2O2 stock solution was verified at 240 nm ( 240 = 43.6 M -1 .cm -1 ).

Bioinformatics

Species tree of the 18 ascomycete genomes (shown in Fig. 1b) was constructed as previously described [START_REF] Morin | Comparative genomics of Rhizophagus irregularis, R. cerebriforme, R. diaphanus and Gigaspora rosea highlights specific genetic features in Glomeromycotina[END_REF]. Core clusters containing only one protein-coding gene per species were identified using FastOrtho [START_REF] Wattam | PATRIC, the bacterial bioinformatics database and analysis resource[END_REF] with the following parameters: 50% identity and 50% coverage.

Each S3 for more details). The phylogenetic analysis of the peroxidase-catalase superfamily relied on 150 sequences encompassing Class I (intracellular peroxidases), II (fungal secreted peroxidases) and III (plant secreted peroxidases) peroxidases and from so-called Hybrid B peroxidases, as previously reported [START_REF] Zámocký | Fungal Hybrid B heme peroxidasesunique fusions of a heme peroxidase domain with a carbohydrate-binding domain[END_REF]. Subsequently to manual curation (removing signal peptides), the sequences were first aligned using MAFFT-DASH (L-INS-i method) (53) and the resulting multiple sequence alignment (MSA) was used to infer a phylogenetic tree using RAxML (1000 bootstraps). The trees were then visualized in iTOL [START_REF] Letunic | Interactive Tree of Life (iTOL) v4: Recent updates and new developments[END_REF] and edited in Illustrator®.

For the gene neighborhood survey, we retrieved 4 genes located upstream and downstream of each Colletotrichum's AA5_2 query (72 sequences). The resulting 569 genes were assigned to 78 different PFAM domains. The frequency of occurrence of a given type of domain in the neighborhood of each AA5_2 phylogenetic clade was then computed and visualized in Excel.

For interrogating the co-occurrence of both Perox and AlcOx-coding genes beyond Colletotrichum species two independent BLAST searches were run against the NCBI nonredundant database, using CorAlcOx and CorPerox as query sequences. 1,000 AlcOx-like (down to 37% sequence identity) and 1,000 Perox-like (25% sequence identity) sequences along with their corresponding source microorganism were retrieved. A cross-comparison of both lists of microorganisms, applying different sequence identity-based thresholds (60% for AlcOx and 30% for Perox) returned the list of species harboring both type of enzymes.

Transcriptomics data were retrieved from publicly available datasets [START_REF] O'connell | Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses[END_REF][START_REF] Gan | Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi[END_REF][START_REF] Dong | Global Genome and Transcriptome Analyses of Magnaporthe oryzae Epidemic Isolate 98-06 Uncover Novel Effectors and Pathogenicity-Related Genes, Revealing Gene Gain and Lose Dynamics in Genome Evolution[END_REF][START_REF] Shimizu | RNA-Seq of in planta-expressed Magnaporthe oryzae genes identifies MoSVP as a highly expressed gene required for pathogenicity at the initial stage of infection[END_REF]. To normalize different dataset longitudinally, we expressed reported sampling time points as infection stages in Fig. 1

Structure Prediction and Preparation

Surface hydrophobicity

Structural homology models were generated with AlphaFold (55) and surface hydrophobicity of selected CROs was computed with the "protein-sol patches" online software [START_REF] Hebditch | Web-based display of protein surface and pH-dependent properties for assessing the developability of biotherapeutics[END_REF]. Average hydrophobicity of the binding surface was determined as follows: using PyMOL 2.4, we selected the residues constituting the entire binding surface of FgrGalOx, CgrAlcOx and equivalent residues (based on MSA) in orthologous enzymes (5 GalOx and 11 AlcOx in total)

as well as those of 4 characterized hydrophobins (PDB 2N4O, 2LSH, 1R2M and 2FZ6) as hydrophobic protein reference. Average hydrophobicity of the selected residues was computed as a GRAVY index score (Kyte-Doolittle method) [START_REF] Kyte | A simple method for displaying the hydropathic character of a protein[END_REF].

Preparation of models for docking experiments

We used Alphafold2 ( 55) to obtain a model for CorAlcOx and CorPerOx. To add the metals and cofactors to AlphaFold2 models, we performed a BLAST search in the Protein Data Bank proteins database and selected the PDB 2EIC (Sequence Identity 47%) to add the copper to CorAlcOx, and the PDB 1MN2 (Sequence Identity 28%) to add the calcium ions and the heme group to CorPerox. Subsequently, we prepared the systems with Schrödinger Protein Preparation Wizard [START_REF] Sastry | Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments[END_REF] to determine the protonation states at pH 7 using PROPKA [START_REF] Olsson | PROPKA3: Consistent treatment of internal and surface residues in empirical p K a predictions[END_REF], and finally relax the systems performing a restrained minimization with convergence criteria for heavy atoms to 0.30 Å using the OPLS_2005 force field.

Protein-Protein Docking

For the generation of protein-protein poses (PPPs), we used PIPER [START_REF] Kozakov | PIPER: An FFT-based protein docking program with pairwise potentials[END_REF] to generate 70,000

PPPs, which were clustered with the in-built clustering protocol using the default RMSD threshold of 9.0 Å and a minimum population of 10 poses per cluster. In parallel, we also used the recently developed Alphafold2 Multimer [START_REF] Evans | Protein complex prediction with AlphaFold-Multimer[END_REF], which only uses sequence information, together with ion and cofactor placements, as described above.

PPPs refinement with PELE

We used the all atom Monte Carlo (MC) software PELE (Protein Energy Landscape exploration), to map intermolecular interactions [START_REF] Borrelli | PELE: Protein energy landscape exploration. A novel Monte Carlo based technique[END_REF][START_REF] Gilabert | Monte Carlo Techniques for Drug Design : The success case of PELE in Methods and Principles in Medicinal Chemistry[END_REF]. PELE follows a heuristic MC approach, generating new conformational proposals using vibration modes of the proteins with translations and rotations of the ligand (CorPerox in this case), and relaxing the system with structure prediction methods, so that the probability of acceptance in the Metropolis criterion (61) remains high [START_REF] Gilabert | Monte Carlo Techniques for Drug Design : The success case of PELE in Methods and Principles in Medicinal Chemistry[END_REF]. Here, we refined the top five PIPER and five AlphaFold2-Multimer models, using a simulation of 250 PELE steps with 256 computing cores (about 25 independent trajectories per model). Each MC PELE step consists of a perturbation and a relaxation stage.

In the first one, a perturbation of the ligand, CorAlcOx in this case, is first performed, including random rotations of 0.01-0.04 rad and translations of 0.25-0.50 Å, and is followed by a backbone perturbation of both proteins following the normal modes directions predicted by an Anisotropic Network Model (ANM). In the second stage, the system is first relaxed by a high resolution side chain prediction including all protein-protein interphase side chains, defined by the region within 3 Å of any heavy atom of the other protein. Afterwards, a global minimization is performed to relax the entire system, providing a final conformation and energy; this energy is then used in a Metropolis importance sampling to accept/reject the MC step.

Exploring the substrate interactions

The C18 substrate was created with Maestro 3D Builder, followed by a Glide docking (62) on the lowest energy PPPs, as predicted with PELE. The docking grid was defined with a cubic box of 30 Å centered in between the copper ion and Tyr120(OH) of the CorAlcOx. Following the Glide rigid receptor docking, we performed a C18 induced fit rescoring simulation involving 40 PELE steps (using 48 computing cores). Each PELE step consisted of a random rotation of 0.01-0.04 rad and translation of 0.05-0.15 Å, using a spherically restrained search space of radius 12 Å centered on the copper ion. The perturbation step was then followed by a relaxation phase including side chain prediction (all side chains inside the spherical space) and a full system minimization.

DNA cloning and strain production

DNA cloning and strain production of the AA5_2 alcohol oxidases (AlcOx) from

Colletotrichum graminearum (CgrAlcOx, Genbank ID XM_008096275.1, Uniprot ID E3QHV8) was already carried out in previous study [START_REF] Tyler) Yin | Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family[END_REF]. The intron-free sequences of the 63). The best-producing transformants were conserved as glycerol stock at -80°C.

Heterologous protein production in flasks

All proteins were first produced in 2 L Erlenmeyer flasks. To this end, single colonies of P.

pastoris X33 expressing each gene of interest were individually streaked on a YPD agar plate containing Zeocin (100 μg.mL -1 ) and incubated 3 days at 30°C. A single colony was then used to inoculate 5 mL of YPD, in a 50 mL sterile Falcon tube and incubated during 5 h (30°C, 160 rpm). This pre-culture was used to inoculate at 0.2% (vol/vol) 500 mL of BMGY medium, in a 2 L Erlenmeyer flask, incubated during approximately 16 h (30°C, 200 rpm) until the OD600 nm reached 4-6. The produced cellular biomass was then harvested by centrifugation (5 min, 16°C, 3,000 x g). For the AlcOx, the cell pellet was then resuspended in 100 mL BMMY medium supplemented with methanol (1%, vol/vol) and CuSO4 (500 µM). The culture was incubated for 3 days (16°C, 200 rpm), with daily additions of methanol (1% added, vol/vol).

The Tandem Peroxidases production conditions were optimized and varied from the standard protocol as follows: the BMMY was supplemented with methanol (3% vol/vol), hemin (25 µM) and CaCl2 (2 mM). The culture was incubated for 3 days (20°C, 200 rpm), with daily additions of methanol (3%, vol/vol) and hemin (25 µM). Then, the extracellular medium was recovered by centrifugation (10 min, 4°C, 3,000 x g) and the supernatant filtrated on 0.45 µm membrane (Millipore, Massachusetts, USA) and stored at 4°C prior to purification.

Heterologous protein production in bioreactors

The upscaled production of CorPerox was carried out in 1.3 L and 7.5 L bioreactors (New Brunswick BioFlo 115 fermentor, Eppendorf, Germany) as per the P. pastoris fermentation process guidelines (Invitrogen) with the following optimizations: the glycerol fed-batch phase was replaced by a sorbitol and methanol transition phase, besides 200 µM (1.3 L bioreactor)

and 150 µM (7.5 L bioreactor) of hemin were added to the methanol solution. CaCl2 (10 mM final) was added to the crude protein solution before being either directly purified or flashfrozen in liquid nitrogen and stored at -80°C. We verified that flash-freezing did not cause any activity loss, for both AlcOx and Perox enzymes.

Protein purification

The

Enzyme assays

For screening the substrate specificity of CRO-AlcOx enzymes, the alcohol substrates were prepared in sodium phosphate buffer ( = 20,000 M -1 .cm -1 ) converted into non-chromogenic product RB5 ox ; veratryl alcohol (500 M) converted into veratraldehyde ( 310 = 9,300 M -1 .cm -1 ). For testing the manganese peroxidase activity, CorPerox was mixed with Mn(II)SO4 (1 mM final) in tartrate buffer (50 mM, pH 2.0 to 5.0) and the formation of Mn 3+ -tartrate complex upon addition of H2O2 (100 µM) was followed at 238 nm ( 238 = 6,500 M -1 .cm -1 ), as previously described [START_REF] Martínez | Purification and catalytic properties of two manganese peroxidase isoenzymes from Pleurotus eryngii[END_REF].

All activities were expressed as Vi/E (s -1 ), i.e. the initial rate (Vi, moles of H2O2 consumed per second) divided by the amount of enzyme (in moles).

CorPerox stability over time was carried out by monitoring the peroxidase activity of CorPerox samples (50 M) stored in sodium acetate buffer (50 mM, pH 5.2), at 4°C, in the presence of varying concentrations of CaCl2 (0-500 mM). The peroxidase activity of these samples was measured as described above (final concentration of 0.

Electron Paramagnetic Resonance (EPR)

EPR spectra were recorded on frozen solutions (120K) using a Bruker Elexsys E500 spectrometer operating at X-band equipped with a BVT 3000 digital temperature controller.

The following acquisition parameters were used: modulation frequency 100 kHz; modulation amplitude 5 G; gain 87 dB; and microwave power, 20 mW. EPR spectra were simulated using the EasySpin toolbox developed for Matlab [START_REF] Stoll | EasySpin, a comprehensive software package for spectral simulation and analysis in EPR[END_REF]. CorAlcOx (100 µM final), prepared in sodium phosphate buffer (50 mM, pH 7.0), in the absence or presence of CorPerox (100 µM final), was flash-frozen in liquid nitrogen and continuous-wave EPR spectra were recorded.

CorAlcOx and CorPerox were placed in contact for various amount of time (2.5 min, 15 min) before flash-freezing the solution. Controls containing buffer only or the CorPerox were also carried out.

Analysis of cuticular waxes from cucumber cotyledons

Cuticular waxes were extracted from 2-weeks old cotyledons by immersing 6 intact cotyledons for 30 s in chloroform in a glass beaker. Chloroform was evaporated under a stream of nitrogen gas and wax extracts were derivatized using N,O-Bis (trimethylsilyl)trifluoroacetamide (BSTFA) and analyzed by gas chromatography coupled to mass spectrometry (GC-MS) as previously described [START_REF] Ménard | Histone H2B monoubiquitination is involved in the regulation of cutin and wax composition in Arabidopsis thaliana[END_REF].

Strains and media

Strain 104-T (MAFF240422) of C. orbiculare was used as the wild-type strain. All strains used in this study are listed in Table S4 [START_REF] Kodama | The morphogenesisrelated NDR kinase pathway of Colletotrichum orbiculare is required for translating plant surface signals into infection-related morphogenesis and pathogenesis[END_REF][START_REF] Fukada | Colletotrichum orbiculare Regulates Cell Cycle G1/SProgression via a Two-Component GAP and a GTPase to Establish Plant Infection[END_REF][START_REF] Tanaka | Saccharomyces cerevisiae SSD1 orthologues are essential for host infection by the ascomycete plant pathogens Colletotrichum lagenarium and Magnaporthe grisea[END_REF][START_REF] Kimura | A Novel Transformation System for Pyricularia Oryzae: Adhesion of Regenerating Fungal Protoplasts to Collagen-Coated Dishes[END_REF][START_REF] Ishida | Relation of Temperature to Germination of Conidia and Appressorium Formation in Colletotrichum lagenarium[END_REF]. C. orbiculare strains were cultured on 3.9% PDA (Nissui) at 24°C in darkness. For genetic manipulation, Escherichia coli DH5αcompetent cells were maintained on Luria-Bertani (LB) agar at 37°C. For fungal transformation, Agrobacterium tumefaciens C58C1 was maintained on LB agar at 28°C.

Transformations of C. orbiculare [START_REF] Kodama | The morphogenesisrelated NDR kinase pathway of Colletotrichum orbiculare is required for translating plant surface signals into infection-related morphogenesis and pathogenesis[END_REF][START_REF] Fukada | Colletotrichum orbiculare Regulates Cell Cycle G1/SProgression via a Two-Component GAP and a GTPase to Establish Plant Infection[END_REF] were carried out as previously described.

Strain construction

Primers and plasmids used in this study are listed in Table S5. For construction of CorAlcOx deletion strains, 1.1-kb upstream and 1.0-kb downstream flanking sequences and a 1.0-kb fragment of the neomycin-resistance cassette was amplified with the respective primer pairs.

For construction of CorPerox deletion strains, 1.1-kb upstream and downstream flanking sequences and a 1.4-kb fragment of the hygromycin-resistance cassette was amplified. These three fragments were inserted into linearized pPZP-PvuII using the In-Fusion HD cloning kit (Clontech). The same procedures were used for construction of M. oryzae gene deletion strains.

For construction of CorAlcOx-mCherry and CorPerox-GFP gene fusion, a 5.9-kb CorAlcOx-CorPerox fragment containing 1.1-kb downstream flanking sequences was inserted into linearized pPZP-PvuII-SUR, and the mCherry and GFP fragments were inserted.

For construction of CorAlcOx-mCherry overexpression strains, a 4.0-kb CorAlcOx-mCherry fragment containing its 1.1-kb downstream flanking sequence was amplified from pPZP-AlcOx-mCherry-Perox-GFP-S, and fused to linearized pCAMSUR-TEF (71) containing the translation elongation factor promoter of Aureobasidium pullulans [START_REF] Vanden Wymelenberg | Expression of green fluorescent protein in Aureobasidium pullulans and quantification of the fungus on leaf surfaces[END_REF]. The same procedures were used for construction of CorPerox-GFP overexpression strains.

Plant infection

Infection assays on detached cucumber leaves (Cucumis sativus L. 'Suyo') with conidial suspension (1×10 5 conidia/ml in distilled water) of C. orbiculare were performed as previously described (73). The inoculated leaves were incubated in a humid environment for 5 days at 24°C. For testing the effect of n-octadecanal supplementation, conidial suspensions (1×10 5 conidia/mL in 10 μM n-octadecanal dissolved in 1% ethanol or 1% ethanol as a control) were spotted onto detached cucumber leaves, and incubated in a humid environment for 5 days at 24°C.

Microscopy

For observation of appressorium formation of C. orbiculare, a conidial suspension (5×10 5 conidia/ml) was placed on a multiwell glass slide or cover glass (Matsunami Glass), respectively. Cells were incubated in a humid box for 24 h at 24°C in the dark. Observation of penetration hyphae on cucumber cotyledons and cellophane membranes were performed as previously described [START_REF] Kodama | The morphogenesisrelated NDR kinase pathway of Colletotrichum orbiculare is required for translating plant surface signals into infection-related morphogenesis and pathogenesis[END_REF]. Appressorial cytorrhysis assay was conducted based on a previous procedure [START_REF] Tanaka | Saccharomyces cerevisiae SSD1 orthologues are essential for host infection by the ascomycete plant pathogens Colletotrichum lagenarium and Magnaporthe grisea[END_REF].

A confocal laser scanning microscope LSM900 with Airyscan 2 (Carl Zeiss) equipped with a Plan Apochromat 63×/1.4 Oil DIC objective (Carl Zeiss) was used to acquire confocal microscopic images. Excitation/emission wavelengths were 488 nm/490-556 nm for GFP and 561 nm/565-630 nm for mCherry. Images were acquired and processed using ZEN Software (Version 3.1; Carl Zeiss) and Imaris (Version 9.3.1; Bitplane). For detection of appressorial actin assembly, cells were observed using a Zeiss Axio Imager M2 Upright microscope (Carl Zeiss) equipped with a Plan Apochromat 100× oil immersion lens, an Axio Cam MRm digital camera and excitation/barrier filter set of 595 nm/620 nm for RFP. Images were acquired using Axiovision 4.8. Bright-field microscopy was performed using a Nikon ECLIPSE E600 microscope equipped with a 40× water immersion lens (Nikon) and an OLYMPUS DP74 digital camera system.

Microarray analysis

For sampling appressoria, the abaxial surface of cucumber cotyledons or cellophane membranes (Wako Chemicals) were inoculated with 10 μL droplets or 10 mL, respectively, of a conidial suspension (1 × 10 6 conidia/ml), and then incubated at 24°C in a humid box. After 16 h, the cellophane was frozen in liquid nitrogen. After 24 h, the lower epidermis of the cotyledons was peeled off. All samples were ground in liquid nitrogen, and total RNA was prepared using the Maxwell RSC Plant RNA Kit (Promega) and the Agilent Plant RNA Isolation MiniKit (Agilent Technologies). Microarray analyses were performed as described previously (13) using the C. orbiculare (8 × 60,000, 13,352 independent probes, Design ID:

Fig. 1 .

 1 Fig. 1. Genomic and transcriptomic analysis of the Perox-AlcOx pair. (A) The multistage plant infection process of appressorium-forming fungi: appressorium (AP), penetration (PP), biotrophic (BP) and necrotrophic (NP) phases. (B) Phylogenomic occurrence of the perox-alcox pair and consensus genomic environment (within each clade) amongst pathogenic ascomycetes (black bold lettering indicates species for which transcriptomics data are shown in panel C, and grey lettering indicates absence of the gene pair. Note that analogous tandem oxidases systems may occur in other fungi but were possibly not detected due to the stringent sequence identity thresholds used to define AlcOx and Perox-encoding genes. Selected Colletotrichum species and associated consensus sequences are representative of their respective species complexes; see fig. S3). (C) Time-course transcriptomic analysis of the Tandem Peroxidases (Perox, green) and AlcOx (orange) encoding genes (gene accession numbers in TableS1) during plant infection for different pathosystems (Myc., Mycelium)[START_REF] O'connell | Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses[END_REF][START_REF] Gan | Comparative genomic and transcriptomic analyses reveal the hemibiotrophic stage shift of Colletotrichum fungi[END_REF][START_REF] Dong | Global Genome and Transcriptome Analyses of Magnaporthe oryzae Epidemic Isolate 98-06 Uncover Novel Effectors and Pathogenicity-Related Genes, Revealing Gene Gain and Lose Dynamics in Genome Evolution[END_REF][START_REF] Shimizu | RNA-Seq of in planta-expressed Magnaporthe oryzae genes identifies MoSVP as a highly expressed gene required for pathogenicity at the initial stage of infection[END_REF].

Fig. 2 .

 2 Fig. 2. Biochemical and biophysical evidence for interplay between CorPerox and CorAlcOx. (A) Reaction mechanism of CROs showing activation of the resting, inactive form of the enzyme via formation of a tyrosine radical, yielding the Cu(II)-radical active form. The latter will oxidize an alcohol into the corresponding aldehyde followed by regeneration of the active form via the two-electron reduction of O2 into H2O2. (B) CorAlcOx oxidation rate of benzyl alcohol in presence of varying amounts of HRP or CorPerox. (C) Activity of CorPerox-CorAlcOx on crude extract of cucumber waxes and derived long-chain aliphatic alcohols, monitored via the production of chromogenic resorufin (RF), product of the Perox-catalyzed oxidation of Amplex-Red by H2O2, the latter being the co-product of AlcOx-catalyzed oxidation of primary alcohols into aldehydes. (D) EPR spectra of inactive CorAlcOx before (orange curve) and after mixture with CorPerox (black curve). EPR parameters of the Cu(II) inactive form: gz = 2.270, 𝐴 𝑧 𝐶𝑢 = 171 10 -4 cm -1 , gx = 2.047, 𝐴 𝑥 𝐶𝑢 < 50 10 -4 cm -1 , gy = 2.054, 𝐴 𝑦 𝐶𝑢 < 50 10 - 4 cm -1 and super-hyperfine coupling constant corresponding to two N-ligands A N = 43 10 -4 cm -1 . (E)

Fig. 3 .

 3 Fig. 3. In vivo characterization of the role of the Perox-AlcOx pair during plant infection. A-B, Infection phenotypes (A) and violin plot of the necrotic lesions size (B), at 5 dpi, for the wild-type (WT) and perox/alcox deletion mutants of C. orbiculare on intact cucumber (Cucumis sativus) cotyledons. For each strain/plant combination, at least 60 inoculations were carried out. (C) Relative frequency (WT set to 100%) of development of normal appressorium (AP) and penetrating hyphae (PH) by C. orbiculare strains on cucumber cotyledons. Data are presented as average values (>300 appressoria for each replicate, n = 3 independent biological replicates) and error bars show s.d.. (D-E) Localization of the C. orbiculare tandem oxidases on cotyledons at 48 hpi, in the presence of appressoria cells (D) or after detachment of the appressoria from the surface of the leaf (E). See movie S1 for a 3D view. In panel D, white arrowheads indicate enzyme accumulation at the appressorium pore. In panel E, enlarged insets show a ring-like localization of the enzymes (scale bar = 1 µm). In panels B-C, a one-tailed independent t-test for each mutant vs WT was applied (***P<0.001, ****P<0.0001).

Fig. 4 .

 4 Fig. 4. Effect of n-octadecanal on the pathogenicity of perox/alcox deletion mutants. The figure shows the effect of the addition of 10 μM n-octadecanal (denoted as +C18) on necrotic lesions size formed by C. orbiculare WT and mutant strains on cucumber cotyledons (24°C, 5 days) (A) and on the formation of penetration hyphae (B). Control experiments, i.e. 1% ethanol without addition of noctadecanal, are denoted as Ref. In panel A, for each strain, we show one representative cotyledon image of the infection phenotype and, below, a violin plot of lesions sizes based on 66 inoculation sites per condition (carried out over 4 independent biological replicates). In panel B, data are presented as mean percentages (>300 appressoria per replicate, n = 4 independent biological replicates) and error bars show s.d. In panels A and B, for each strain, a one-tailed independent t-test of +C18 vs Ref was applied.

Figure 5 .

 5 Figure 5. Loss of the Perox-AlcOx pair affects fungal gene expression during appressoriummediated penetration of C. orbiculare. (A) Venn diagrams illustrating the differentially regulated (upor down-regulated) genes between cucumber vs cellophane that are under control of the AlcOx-Perox pair. Top: green circle indicates the number of downregulated genes (46 genes) in Δpair compared with wild-type strain (WT), both incubated on cucumber; magenta circle denotes the number of genes (1608 genes) upregulated on cucumber vs cellophane in the WT strain. Bottom: Magenta circle indicates the number of upregulated genes (86 genes) in Δpair compared with WT, both incubated on cucumber; green circle shows the number of downregulated genes (2887 genes) in the WT strain on cucumber vs cellophane. In the top and bottom Venn diagrams, the overlap indicates that 32 and 26 genes are respectively upregulated or downregulated in the presence of the plant and also under control of the AlcOx-Perox pair. (B) Class-wise distribution of the number of genes that were differentially expressed in appressoria of the Δpair double mutant compared with WT on cucumber cotyledons. (C) Hierarchical clustering heatmap of differentially expressed genes encoding CAZymes and SSPs in the Δpair mutant compared to WT on cucumber cotyledon. The expression values per gene were median normalized. The columns were clustered by Euclidean distance. The expression levels of up-(magenta) and downregulated (green) genes are shown as log2-transformed values. On the right-hand side of the figure, we show the presence of a predicted signal peptide (Y, Yes; N, No), phylogeny-based substrate specificity predictions for CAZymes (PCW, plant cell wall; FCW, fungal cell wall; GAL, galactose-containing compounds; Nd, not determined) and SSPs/CBM50s that are putative effectors (pale green circles)

  localized co-secretion of the Perox-AlcOx would have evaded classical 'omic approaches. The specific occurrence of the Perox-AlcOx pair in most Colletotrichum and Magnaporthe species raises the possibility that functionally equivalent, coupled oxidative enzymatic mechanisms may operate in other appressorium-forming fungal pathogens. We anticipate that the present discovery will open new research avenues, notably on the role of wax-derived compounds in the cascade of events leading to successful infection as well as on the development of oxidasespecific inhibitors as surface-acting, anti-penetrant drugs for crop protection.

Fig. 6 .

 6 Fig. 6. Schematic summary of the recruitment of the fungal Perox-AlcOx pair during early plant penetration and its proposed role in the induction of a biochemical cascade. The bottom panel is a zoom-in view illustrating the Perox-AlcOx reaction occurring at the fungus-plant interface and the triggered downstream cascade, highlighting the induction of genes encoding proteins targeting the plant cell wall (pcw) and fungal cell wall (fcw).

  genes coding for the AlcOx from Colletotrichum orbiculare MAFF 240422 (CorAlcOx, Genbank ID TDZ17043.1, Uniprot ID N4UTF2), the AlcOx from Magnaporthe oryzae (MorAlcOx, Genbank ID XM_003719321.1, Uniprot ID G4NG45), the Tandem Peroxidase (Perox) from Colletotrichum orbiculare (CorPerox, Genbank ID TDZ17044.1, Uniprot ID N4UUY4) and the Tandem Peroxidase from Magnaporthe oryzae (MorPerox, Genbank ID XM_003719322.1, Uniprot ID G4NG46) were synthesized after codon optimization for expression in P. pastoris and inserted into a modified pPICZαC vector using XhoI* and NotI restriction sites in frame with the α secretion factor at N-terminus (i.e. without native signal peptide) and with a (His)6-tag at the C-terminus (without c-myc epitope) (Genewiz, Leipzig, Germany). Transformation of competent P. pastoris X33, selection of zeocin-resistant P. pastoris transformants screened for protein production was carried out as described by Haon et al. (

  filtered CorAlcOx and MorAlcOx crude supernatants were adjusted to pH 8.5, filtered on 0.22 µm filters (Millipore, Molsheim, France), and purified by anion exchange chromatography (DEAE) on a HiPrep FF 16/10 column (GE Healthcare, USA). Elution was performed by applying a linear gradient from 0 to 500 mM NaCl (in Tris-HCl buffer 50 mM, pH 8.5) over 20 column volumes, with a flow rate set to 5 mL.min -1 .The filtered CorPerox and MorPerox culture supernatant was adjusted to pH 7.8 just before purification and filtered on 0.22 µm filters (Millipore, Molsheim, France). Depending on the volume to purify, the crude protein sample was either loaded on a His-Trap HP 5-mL column (GE Healthcare, Buc, France) on a HisPrep FF 16/10 column (GE Healthcare) connected to an ÄKTAxpress system (GE Healthcare) equilibrated with HEPES (10 mM, pH 8.0), NaCl (100 mM), CaCl2 (2 mM) and imidazole (10 mM) buffer. Each (His)6-tagged recombinant enzyme was eluted with HEPES (10 mM, pH 8.0), NaCl (100 mM), CaCl2 (2 mM) and imidazole (500 mM) buffer. The Tandem Peroxidases were further purified by size exclusion chromatography, using a HiLoad 26/600 Superdex 200 pg column (GE Healthcare) operated at 2.5 mL/min and with a running buffer containing HEPES (10 mM, pH 8.0), NaCl (100 mM) and CaCl2 (2 mM).After SDS-PAGE analysis, fractions containing the recombinant enzyme were pooled, concentrated and buffer exchanged in sodium phosphate (50 mM, pH 7.0) for the AlcOx or in HEPES (10 mM, pH 8.0), NaCl (100 mM) and CaCl2 (2 mM) buffer for the Tandem Peroxidases.Protein concentrations of CgrAlcOx (52,337 Da,  280 = 101,215 M -1 .cm -1 ), CorAlcOx[START_REF] Zámocký | Fungal Hybrid B heme peroxidasesunique fusions of a heme peroxidase domain with a carbohydrate-binding domain[END_REF]317 Da,  280 = 92,735 M -1 .cm -1 ), MorAlcOx (62,894 Da,  280 = 90,020 M -1 .cm -1 ), CorPerox[START_REF] O'connell | Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses[END_REF]137 Da,  280 = 21,345 M -1 .cm -1 ), and MorPerox (26,290 Da,  280 = 24,450 M -1 .cm -1 ) were determined by the Bradford assay (64) using BSA as reference protein as well as by UV absorption at 280 nm using a Nanodrop ND-200 spectrophotometer (Thermo Fisher Scientific, Massachusetts, USA).

  50 mM, pH 7.0) in 96-well microplates and reactions were initiated by the addition of a pre-mix of CRO-AlcOx (1 nM final concentration), HRP (0.1 mg.mL -1 ) and 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid (ABTS, 500 M) in sodium phosphate buffer (50 mM, pH 7.0). The tested substrates included D-glucose(50 mM final concentration, D-Glc), D-galactose (50 mM, D-Gal), D-raffinose (50 mM, D-Raf), xyloglucan (0.1% mM, XG), butan-1-ol (3 mM), butan-2-ol (3 mM), octan-1-ol (3 mM), decan-1-ol (3 mM), 2,4-hexadiene-1-ol (3 mM, HD-OH), glycol aldehyde dimer (3 mM, GAD), benzyl alcohol (3 mM, BnOH), 4-hydroxybenzyl alcohol (3 mM, p-OH BnOH), vanillic alcohol (3 mM, Van-OH), syringic alcohol (3 mM, Syr-OH) and cinnamyl alcohol (3 mM, Cin-OH). The absorbance of the final reaction (100 L total volume) was monitored at 414 nm using a microplate spectrophotometer (TECAN), thermostated at 23°C. The 414 nm absorbance allows to determine the concentration of ABTS cation radical over time (ABTS •+ ,  414 = 31,100 M -1 .cm -1 ), and in turn the rate of alcohol oxidation, considering a peroxidase reaction stoichiometry for (H2O2:ABTS •+ ) of 1:2 and a CRO-AlcOx reaction stoichiometry for (alcohol:H2O2) of 1:1. For screening the substrate specificity of CorPerox, unless stated otherwise the enzyme (0.125 µM final) was prepared in citrate-phosphate buffer (50 mM, pH 4.0 to 7.0) in 96-well microplates (for wavelength in the visible range) or in 1 mL Quartz cuvettes (for UV range), in the presence of various substrates (vide infra). Reactions were initiated by the addition of H2O2 (100 µM final), incubated at 23°C, and monitored spectrophotometrically at the wavelengths indicated below. The tested peroxidase substrates included: ABTS (500 M) converted into ABTS •+ ( 414 = 31,100 M -1 .cm -1 ); 2,6-dimethoxyphenol (2,6-DMP, 500 M) converted into hydrocoerulignone ( 469 = 53,200 M -1 .cm -1 ); guaiacol (500 M) converted into the final product tetraguaiacol ( 470 = 26,600 M -1 .cm -1 ); Reactive Black 5 (RB5, 100 M,  600

a

  previously published protocol[START_REF] Ribeaucourt | Controlled Oxidation of Long-Chain Primary Aliphatic Alcohols by a Fungal Copper-Radical Alcohol Oxidase[END_REF], the reaction mixture was then acidified by addition of 10 µL HCl (12 M). Products and possible remaining substrate were extracted by adding 500 µL of hexane (containing 1 mM of internal standard dodecane), followed by shaking and centrifugation for 5 min at 3,000 x g. The organic layer was transferred into a new vial and analyzed with a GC-2010 Plus apparatus (Shimadzu, Japan) equipped with a flame ionization detector (FID) and a DB-5 capillary column (30 m x 0.25 mm x 0.25 µm; Agilent). Nitrogen (200 kPa) was used as carrier gas. The injector and detector temperatures were set at 250°C. After injection (2 µL sample), the analytes were separated by applying the following temperature program: step 1) from 65°C to 250°C over 9.25 min (i.e. 20°C/min); step 2) plateau at 250°C for 6 min. For quantitation, standard curves of octadecanol, n-octadecanal and octadecanoic acid were prepared by following the same procedure.

  060762) oligo microarray, according to the Agilent 60-mer Oligo Microarray ProcessingProtocol (Agilent Technologies). The normalization condition: i) intensity-dependent Lowess normalization; ii) data transformation, measurements less than 0.01 were set to 0.01; iii) perchip 75 th -percentile normalization of each array; and iv) the expression values per gene were median normalized. The normalized data were subjected to a t test, with statistically significant gene sets defined as those giving P values less than 0.05. The differentially-regulated genes (fold change > 2 and P< 0.05) were selected and were used for further analysis. Functional classification was based on the Gene Ontology (GO), protein families (Pfam) and C. orbiculare genome information of CAZymes and SSPs. Note that to avoid redundancy in the count of upand downregulated genes, all sequences considered as CAZymes were not counted in the broader classes of Hydrolases/Transferases/Oxidoreductases.

  M. oryzae/Rice[START_REF] Dong | Global Genome and Transcriptome Analyses of Magnaporthe oryzae Epidemic Isolate 98-06 Uncover Novel Effectors and Pathogenicity-Related Genes, Revealing Gene Gain and Lose Dynamics in Genome Evolution[END_REF], AP (8 hpi), BP (24 hpi) and NP (48 hpi); M. oryzae/Barley (29), AP (12 hpi), PP (24 hpi), BP (36 hpi) and NP (48 hpi).

as follows (hpi = hours post-infection): for C. fructicola nara gc5/strawberry (27), AP (24 hpi), BP (72 hpi) and NP (144 hpi); for C. orbiculare/N. benthamiana (27), AP (24 hpi), BP (72 hpi) and NP (168 hpi); for C. higginsianum/A. thaliana (26), AP (22 hpi), BP (40 hpi) and NP (60 hpi); C. graminicola/Maize (26), AP (22 hpi), BP (40 hpi) and NP (60 hpi);

  5 M CorPerox), using ABTS (500 M) and H2O2 (100 M) as substrates, in citrate-phosphate buffer (50 mM, pH 4.0), at 23°C. Michaelis-Menten kinetic parameters of CorPerox were determined by measuring the peroxidase initial rate, as described above, in the presence of ABTS (500 M) and varying concentrations of H2O2 (0-1200 M), in citrate-phosphate buffer (50 mM, pH 4.0), at 23°C. Enzymatic reactions were carried out in 4 mL-clear borosilicate glass vials closed by screw caps with PTFE septum (500 µL final reaction volume). CorAlcOx (2 µM final) was mixed with CorPerox (2 µM) in sodium phosphate buffer (50 mM, pH 7.0). The reaction was initiated by the addition of octadecanol (0.3 mg.mL -1 , eq. 1.1 mM final) and the mixture was incubated at 23°C at 190 rpm in an Innova 42R incubator (New Brunswick, USA), during 1 h. Following

	Experimental data could be fit to the standard Michaelis-Menten equation (residual standard
	error = 0.019).
	AlcOx activation by the peroxidases was assayed by monitoring changes in absorbance at 254
	nm upon oxidation of benzyl alcohol (1.5 mM) into benzaldehyde by the CRO (10 nM final
	concentration), in the presence of varying concentrations of peroxidase (0-1000 nM). Reactions
	were carried out in sodium-phosphate buffer (50 mM, pH 7.0), at 23°C, in UV-transparent
	cuvettes (1 mL reaction volume). The reactions were initiated by addition of the CRO and
	vigorously mixed by pipetting up and down. The absorbance was measured using an Evolution
	201 UV-Vis spectrophotometer (Thermo-Fisher). The concentration of benzaldehyde was
	calculated as [Benzaldehyde]t = (Abs 254 nm	t -Abs 254 nm t0)/( 254	benzaldehyde- 254 BnOH), where
	 254 benzaldehyde = 8,500 M -1 .cm -1 and  254 BnOH = 150 M -1 .cm -1 .
	Gas chromatography (GC) analysis

Acknowledgments: We thank Drs. Morin and Lebreton (INRAE Nancy, France) for allowing us to use their in-house tool "Genocomp" for species tree construction. We thank Pamela Gan (RIKEN-CSRS, Japan) for sharing transcriptome data. Funding: This study was supported by the "Agence Nationale de la Recherche" and by the Natural Sciences and Engineering Research Council of Canada through the ANR-NSERC project "FUNTASTIC" (ANR-17-CE07-0047, STPGP 493781-16). We are grateful to MANE & Fils and the "Association Nationale Recherche Technologie" (ANRT) for funding the Ph.D. fellowship of D.R. (grant no. 2017/1169). Work in Japan was supported by the Japan Society for the Promotion of Science Grants-in-Aid for Scientific Research -KAKENHI, Grant Numbers 15H05780 and 20H02989 to Y.K., and 20K15529 to S.K. Author contributions: B.B., S.K., R.O'C., M.L., Y.K., and J.G.B. conceived the work. J.G.B and Y.K. coordinated the study. B.B., S.K., H.H., D.R., A.J.S., V.G., R.O'C., M.L., Y.K., and J.G.B. designed the experiments. B.B., S.K., H.H., D.R., M.H., S.G., F.B., T.N., A.M.D.R., and A.J.S. performed the experiments. B.B., S.K., H.H., D.R., A.J.S., S.M.F., H.B., V.G., M.N.R., R.O'C., M.L., Y.K., and J.G.B. analyzed the data. B.B., S.K., Y.K., and J.G.B. wrote the original draft with review and editing from M.L., H.B., and R.O'C. All authors have approved the final version of the paper. Competing interests: Authors declare no competing interests; and Data and materials availability: All data is available in the main text or the supplementary materials.

Supplementary Materials:

Movies S1