Peroxin FgPEX22-Like Is Involved in FgPEX4 Tethering and Fusarium graminearum Pathogenicity

Peroxisomes are essential organelles that play important roles in a variety of biological processes in eukaryotic cells. To understand the synthesis of peroxisomes comprehensively, we identified the gene FgPEX22-like, encoding FgPEX22-like, a peroxin, in Fusarium graminearum. Our results showed that although FgPEX22-like was notably different from other peroxins (PEX) in Saccharomyces cerevisiae, it contained a predicted PEX4-binding site and interacted with FgPEX4 as a rivet protein of FgPEX4. To functionally characterize the roles of FgPEX22-like in F. graminearum, we performed homologous recombination to construct a deletion mutant (ΔPEX22-like). Analysis of the mutant showed that FgPEX22-like was essential for sexual and asexual reproduction, fatty acid utilization, pathogenicity, and production of the mycotoxin deoxynivalenol. Deletion of FgPEX22-like also led to increased production of lipid droplets and decreased elimination of reactive oxygen species. In addition, FgPEX22-like was required for the biogenesis of Woronin bodies. Taken together, our data demonstrate that FgPEX22-like is a peroxin in F. graminearum that interacts with PEX4 by anchoring PEX4 at the peroxisomal membrane and contributes to the peroxisome function in F. graminearum.


INTRODUCTION
Fusarium head blight (FHB), caused by F. graminearum (teleomorph Gibberella zeae), is an acute disease of wheat and barley worldwide. In addition to causing severe crop yield losses, F. graminearum produces deoxynivalenol (DON), a mycotoxin that also acts as a virulence factor to facilitate wheat infection, and is a menace to the health of humans and animals (Desjardins et al., 1993;Pestka and Smolinski, 2005;Dean et al., 2012).
Peroxisomes are single-membrane-bound organelles present in most eukaryotic organisms and are relevant to multifarious metabolic conversions. For instance, in eukaryotic cells, peroxisomes are involved in methanol oxidation, disposal of ROS, and utilization of carbon sources (Lazarow and Fujiki, 1985;Wanders, 2004;Gould et al., 2010). In yeasts, peroxisomes are the unique sites of fatty acid β-oxidation (Hiltunen et al., 2003), whereas in plants, they are essential for host resistance, embryo development, synthesis of phytohormones, and the glyoxylate cycle (Hu et al., 2012;Shabab, 2013). In addition, peroxisomes are involved in several physiological processes in mammals, such as the synthesis of cholesterol, plasmalogens, and bile acids (Wanders and Waterham, 2010). Zellweger syndrome (ZS), a prototypic peroxisome biogenesis disorder (PBD) with the most severe phenotype in humans, may result in the absence of functional peroxisomes (Wanders, 2004;Faust et al., 2005).
Proteins related to peroxisomal biogenesis are termed peroxins and are encoded by PEX genes. So far, more than 30 PEXs have been identified in various organisms (Distel et al., 1996;Pieuchot and Jedd, 2012). As peroxisomes do not contain any genetic material, their peroxisomal membrane proteins (PMPs) and matrix proteins are encoded in the nucleus, synthesized in the cytoplasm, and then imported into the peroxisomes. For example, PEX5 is a cycling receptor for the import of PMPs containing the peroxisomal targeting signal type 1 (PTS1). The PEX5 import cycle involves the following steps: protein containing PTS1 is recognized by PEX5 in the cytosol, the PEX5-cargo complex docks at the peroxisomal membrane, the cargo is translocated into the peroxisomal lumen, and after dissociation PEX5 is recycled back to the cytosol for a new import cycle (Albertini et al., 1997;Wang et al., 2003;Stanley et al., 2006;Platta et al., 2007Platta et al., , 2008Meinecke et al., 2010). During the last step, PEX5 must be monoubiquitinated by the ubiquitin-conjugating enzyme PEX4 and its membrane-anchor PEX22 (Collins et al., 2000;Zolman et al., 2005). We previously showed that PEX4 is indispensable for peroxisome function in F. graminearum and its pathogenesis (Zhang et al., 2019b). More attention has been paid to the function of peroxins in fungi that are pathogenic to plants. There are several peroxins in F. graminearum involved in the peroxisome life cycle that have been characterized. For instance, FgPEX1,FgPEX2,FgPEX4,FgPEX5,FgPEX6,FgPEX7,FgPEX10,FgPEX12,FgPEX13,FgPEX14, and FgPEX33 are involved in mycotoxin biosynthesis, pathogenicity, and pexophagy (Min et al., 2012;Chen et al., 2018;Zhang et al., 2019a,b;Wang et al., 2020). In Magnaporthe oryzae, MoPEX5, MoPEX6, MoPEX7, MoPEX14, MoPEX19, and MoPEX11 family peroxins are involved in matrix protein import and peroxisomal fission processes (Deng et al., 2013;Wang et al., 2013Wang et al., , 2015Li et al., 2014). Peroxisome studies in other filamentous fungi, such as Neurospora crassa, Colletotrichum orbiculare, and Aspergillus nidulans, are becoming more common (Hynes et al., 2008;Fujihara et al., 2010;Managadze et al., 2010).
In S. cerevisiae, PEX22 plays a key role in tethering the PEX4, a ubiquitin-conjugating enzyme, to the peroxisome and is associated with PEX5 receptor recycling (Williams et al., 2012). Previous studies have shown that most of the PEX genes was contained in filamentous fungi, but except for PEX15, PEX17, PEX18, PEX21, and PEX22 (Heinemann et al., 2004). Until the functions of FAM1 in C. orbiculare were described, functions similar to those of PEX22 was discovered (Kubo et al., 2015). However, whether a similar protein exists in F. graminearum and whether it would possess PEX22 function is unclear.
Herein, we identified FgPEX22-like, which encodes a peroxin FgPEX22-like, which is the functional ortholog of the PEX22 of Saccharomyces cerevisiae. FgPEX22-like was able to interact with FgPEX4 and was essential for the subcellular localization of FgPEX4. Functional analysis provided evidence that FgPEX22like plays important roles in sexual and asexual reproduction, carbon source utilization, pathogenicity, and cell wall integrity. Importantly, FgPEX22-like was essential for the biosynthesis of Woronin bodies.

Fungal Strains and Growth Conditions
For performing mycelial growth assays, the WT F. graminearum strain PH-1, FgPEX4 deletion mutants, and other transformants generated in this study were grown on potato dextrose agar (PDA) medium in a 25 • C incubator. For the aerial hyphal growth assay, all strains were inoculated into test tubes (1.5 cm diameter) containing 5 mL PDA medium and grown at 25 • C for 5 days. To study the integrity of cell membranes and cell walls, all strains were grown as previously reported on CM medium supplemented with 0.01% SDS as a cell membrane-damaging agent and 0.2% Congo red as a cell wall-damaging agent (Gavric et al., 2007;Chayakulkeeree et al., 2008;Qin et al., 2015).

Yeast Two-Hybrid Assay and Co-immunoprecipitation
The Y2H assay was performed using the Matchmaker GAL4 Two-Hybrid System 3 (Clontech) according to the manufacturer's instructions. Full-length complementary DNA (cDNA) of FgPEX22-like and FgPEX4 was PCR amplified using the primer pairs P22-AD-F/22R and P4-AD-F/R, respectively. The resulting PCR products were cloned into pGADT7 and pGBKT7 which were digested with XholI to create FgPEX22like-AD and FgPEX4-BD as the prey vector and bait vector, respectively. A similar method was used to generate FgPEX4-AD and FgPEX22-like-BD. The resulting bait and prey vectors were confirmed by sequencing and co-transfected as pairs. To explore the key action regions of FgPEX4 and FgPEX22-like, truncated cDNA of FgPEX4 and FgPEX22-like were amplified with their respective primers (Supplementary Table 2). These shorter PCR products were cloned into pGADT7 to construct different FgPEX22-like (a-b)-ADs as the prey vector. The same method was used to create FgPEX4 (c-d)-BDs as bait vectors. The truncated bait and prey vectors were also confirmed by sequencing and were co-transfected into yeast in various combinations (Supplementary Table 2). The interaction between pGBKT7-53 (BD-3) and pGADT7-T (AD-1) was used as a positive control.
For co-immunoprecipitation assays, the intracellular region of FgPEX4 and FgPEX22-like were PCR-amplified with the primer pairs 4-GFP-F/R and 22-Flag-F/R (Supplementary Table 2), and then inserted into the vector pFL2 and pFL7 which were digested with Xhol I respectively. The resulting fusion constructs PEX4-GFP and FgPEX22-like-3 × FLAG were verified by DNA sequencing. Constructs PEX22-like-3 × FLAG and PEX22-like-3 × FLAG with PEX4-GFP were transfected into strain PH-1. For transformant selection, G418 (Geneticin) was added at a final concentration of 200 mg/mL. Total proteins were extracted and incubated with GFP beads as previously reported (Li et al., 2017). Total proteins and proteins eluted from the GFP beads (elution) (Kangti Life Technology Co., Ltd., KTSM1334) were analyzed by western blotting using monoclonal anti-FLAG antibody (Abimate medical technology (Shanghai) Co., Ltd., PA9020S) and anti-GFP antibody (Abimate medical technology (Shanghai) Co., Ltd., PA9056S), accordingly. The results were visualized using the enhanced chemiluminescent (ECL) detection system.

Generation of FgPEX22-Like Gene Deletion and the PEX22-Like Complementation Strain
The mutants were generated using the split-marker method (Wang et al., 2011). For the FgPEX22-like gene (FGSG_11970), the sequence was obtained from the F. graminearum database. The upstream (931 bp) and downstream (855 bp) regions flanking the gene were PCR-amplified using primer pairs AF/AR and BF/BR, respectively (Supplementary Table 2). The plasmid pCB1003 harbored hygromycin B resistance gene (HPH) and the primer pair HPH-F/HPH-R was used to amplify HPH (Supplementary Table 2). A fusion cassette containing the FgPEX22-like flanking sequences and the HPH gene was transfected into protoplasts of PH-1 to generate the PEX22-like mutant (Catlett et al., 2003). Southern blot assay for FgPEX22-like deletion mutants was performed using the digoxigenin (DIG)labeled probe and the High Prime DNA Labeling and Detection Starter Kit I (Roche Diagnostics, Mannheim, Germany), as instructed by the manufacturer. Hybridization was performed using a DIG-labeled specific probe (Supplementary Figure 1).
To prepare the complementation strain, a 2,674 bp fragment containing the full-length FgPEX22-like gene sequence and its promoter sequence was PCR-amplified using the primer pair 22CF/22CR (Supplementary Table 2) and then amplimer inserted into pYF11 using the yeast in vivo recombination approach (Bruno et al., 2004;Zhou et al., 2011). The recombinant plasmid pYF11-PEX22-like was transfected into protoplasts of the PEX22-like mutant to produce PEX22like complementation ( PEX22-like -C) strains, which were identified by PCR analysis (Supplementary Figure 1).
To generate the FgPEX22-like and FgPEX4 double-knockout mutants, a FgPEX22-like gene replacement construct was generated with the G418 amplified from pFL2 (Zhou et al., 2011) and transfected into the FgPEX4 mutant PEX4 (Table 1). Transformants resistant to both hygromycin and G418 were screened by PCR analysis .

Conidiation, Germination, and Sexual Reproduction Assays
For conidiation assays, strains were inoculated into CMC medium as previously described (Hou et al., 2002). The chemical compound 4 ,6-diamidino-2-phenylindole (DAPI, 10 µg/mL) and calcofluor white (CFW, 1 µg/mL) were used to stain nuclei and septa of conidia, respectively. After staining, images were captured using a fluorescence microscope (Eclipse 90i, Nikon). To determine the germination rate of conidia, freshly harvested conidia were transferred to the sterile YEPD medium (yeast extract, 10.0 g; peptone, 20.0 g; glucose, 20.0 g; and distilled water to make up the volume to 1,000 mL) for 6 h and observed under a fluorescent microscope (Eclipse 90i, Nikon). For sexual reproduction assays, 7-day-old aerial hyphae growing on specific sporulation medium were compressed in 1 mL sterile 2.5% Tween 60 solution as previously described (Bowden and Leslie, 1999;Jenczmionka et al., 2003). Perithecium formation was examined after 2-3 weeks of incubation at 25 • C.

Plant Infection and Deoxynivalenol Production Assays
Plant infections were assayed on wheat heads and corn silks. Susceptible wheat cultivar Jimai 22 was used in wheat infection assays and was sprayed with a conidial suspension (2 × 10 5 spores/mL) collected from 5-day-old liquid CMC medium (Qin et al., 2015). Wheat heads were photographed and assayed 14 dpi. Fresh corn silks were infected with hyphal plugs, incubated at 25 • C, and examined at 5 dpi (Seong et al., 2005;Gale et al., 2007). To determine DON production, we inoculated three mycelial plugs from each strain into 5 g healthy and aseptic rice grains. After incubating at 25 • C, DON was extracted at 20 dpi as described previously and quantified using a liquid chromatography-mass spectrometer/mass spectrometer (HPLC-MS/MS) system (AB Sciex 5500) as previously described (Mirocha et al., 1998).

Analysis of Fatty Acid Utilization
The carbon source utilization was evaluated using minimal medium containing various carbon sources in lieu of sucrose. The following concentrations were used: 2.5 mM myristic acid (C14), 2.5 mM palmitic acid (C16), 2.5 mM oleic acid (C18), and 2.5 mM erucic acid (C22) as previously described (Leslie and Summerell, 2007). Emulsifier NP40 was added to the minimal medium containing palmitic acid, oleic acid, and erucic acid. Colony diameters were measured after 3.5 days incubation at 25 • C.

Light Microscopy and Transmission Electron Microscopy Observations
Lipid droplets (LD) in hyphae were stained using Nile red (50 µg/mL) as reported previously (Lu et al., 2009). Hyphae grown on PDA plates at 25 • C for 3 days were collected and subjected to ultrastructural analysis. The collected fungal mass was treated and examined by transmission electron microscopy (TEM) (JEM-1400 Plus, JEOL, Tokyo, Japan).
For subcellular localization of FgPEX4, peroxisome membrane protein 70 (PMP70), and Woronin body protein hex1 (HEX1), coding sequences of FgPEX4, FgPMP70, and FgHEX1 were PCR-amplified using the primer pair 4-GFP-F/R, PMP70-GFP-F/R, and HEX1-GFP-F/R, respectively (Supplementary Table 1). The amplimers were then inserted into pYF11 using the yeast in vivo recombination approach (Bruno et al., 2004;Zhou et al., 2011). The recombinant plasmids were transfected into protoplasts of the WT and mutant strains. The transformants were verified by PCR analysis using the appropriate primers. Transformants were observed using a fluorescence microscope.

ROS Detection
The ROS generated in the hyphae of PH-1 and PEX22-like strains was assessed using NBT (nitroblue tetrazolium chloride) following the growth of the strains on CM at 25 • C for 3 days.
*Number of conidia in 100 mL of carboxymethylcellulose (CMC) cultures were examined after incubation for 5 days. **Conidia were incubated in YEPD medium at 25 • C for 6 h. † HPLC-MS/MS analysis of DON produced in the PH-1 and mutants. ‡ Expression levels of several genes at the level of transcription. The relative expression level of GAPDH gene was used as an internal control. The gene expression in PH-1 was set to 1.0 (P < 0.05). The different letter on the bars for each treatment indicates significant difference at P < 0.05 by Duncan's multiple range test.
Each plate was then stained with 20 mL of 0.2% NBT solution and incubated in the dark at 28 • C for 45 min. The liquid stain was drained from the plates, which were then washed with ethanol. The plates were incubated again for 45 min in the dark at 28 • C prior to imaging.

Quantitative Real-Time PCR
Total RNA was isolated from hyphae of the WT, mutant, and complementation strains using TransZol Up (TransGen Biotech, Beijing, China). The quantitative real-time PCR (qRT-PCR) experiments were performed as the manufacturer's instructions (Vazyme Biotech Co., Jiangsu, China). The glyceraldehyde 3phosphate dehydrogenase (GAPDH) gene of F. graminearum was used as the internal control. Relative expression levels for each gene were calculated using the 2 − CT method (Livak and Schmittgen, 2001).

Statistical Analysis
Each experiment was performed in three individual replicates. Data are presented as mean ± standard error values and the differences among variables were analyzed using Duncan's multiple range test. Results with p < 0.05 were considered statistically significant.

Identification of FgPEX22-Like in Fusarium graminearum
The homologs of PEX22 (FGSG_11970) were identified using the Fusarium genome database 1 . The predicted gene was 1,631 bp, encoding 363 amino acids. FgPEX22-like shared only 14.34% identity with C. orbiculare and 17.93% identity with S. cerevisiae PEX22. Therefore, we named the gene FgPEX22-like. SMART-PFAM analysis revealed the protein encoded by PEX22-like contained a transmembrane domain and a possible PEX4 interacting region similar to that in C. orbiculare and S. cerevisiae (Supplementary Figure 1A). The deduced FgPEX22-like protein was distributed in one branch in the phylogenetic tree along with other fungal PEX22 proteins (Supplementary Figure 1B). To investigate whether FgPEX22-like interacted with FgPEX4 in F. graminearum, a yeast two-hybrid assay was used. Yeast transformants expressing binding domain (BD)-FgPEX4 as the bait and the activation domain (AD) pGADT7 alone as the prey, could not grow on the Sabouraud dextrose (SD)-Leu-Trp-His-Ade plates, excluding the possibility of self-activation. This was also true for the reverse situation. The sets of interacting pairs, BD-FgPEX4 and AD-FgPEX22-like or AD-FgPEX4 and BD-FgPEX22-like, were found to bind with each other (Figure 1A).
To further verify these results, the intracellular region of FgPEX4 was fused with pFL2 [PEX4-green fluorescent protein (GFP)] and co-transfected into F. graminearum strain PH-1 with the FgPEX22-like-3 × FLAG fusion construct. In western blot analysis of total proteins from the strain containing PEX22-3 × FLAG, a 40-kDa band was detected using the anti-FLAG antibody, while the elution proteins were not detected. This showed the beads had been washed prior to protein elution and that there were no other contaminants, confirming the subsequent results were reliable. In the transformants coexpressing FgPEX22-like-3 × FLAG and FgPEX4-GFP, 19-kDa and 40-kDa bands were detected using the anti-FLAG and anti-GFP antibodies, respectively ( Figure 1B). Therefore, FgPEX22like interacted with FgPEX4 in F. graminearum. Taken together, these results demonstrated the two proteins could interact directly with each other.
Based on the predicted secondary structure and domain structure, we explored the key functional regions of FgPEX4 and FgPEX22-like. FgPEX22-like contained an N-terminal transmembrane (TM) fragment and a large region spanning an unknown fold that was exposed in cytosol. Several truncated variants of FgPEX22-like were established and their interactions with FgPEX22-like and FgPEX4 were evaluated using the yeast two-hybrid system. Our results showed that the region spanning residues 102-282 was essential for FgPEX4 binding; moreover, residues were part of the soluble region of FgPEX22like. Meanwhile, we also discovered that the first 11 amino acids of FgPEX4 were not necessary for FgPEX22-like binding (Supplementary Figure 2). FIGURE 1 | Interaction of FgPEX4 and FgPEX22-like. (A) Y2H assays were performed to detect interactions between FgPEX4 and FgPEX22-like. All constructed yeast transformants were tested for growth on SD-Leu-Trp-His-Ade and evaluated for β-galactosidase (LacZ) activity. The interaction of AD-1 and BD-3 was used as the positive control. All transformants were diluted to different concentrations (cells/ml) and plated onto SD-Leu-Trp-His-Ade media. 22 and 4 referring to FgPEX22-like and FgPEX4 respectively. (B) Co-immunoprecipitation assays. Total proteins were extracted from transformants co-expressing PEX4-GFP and PEX22-3 × FLAG constructs. The proteins were then eluted from the anti-GFP M2 beads (elution). The immunoblots were incubated with monoclonal anti-FLAG or anti-GFP antibody, as indicated.

FgPEX22-Like Acted as a Rivet Protein of FgPEX4
To further understand the relationship between the two proteins, we compared the subcellular localization of wild-type (WT) FgPEX4 to that of mutant strains. Interestingly, the results showed that a deficiency of FgPEX22-like could result in abnormal subcellular localization of FgPEX4. FgPEX4-GFP fusions was distributed in punctate patterns in transformed PH-1, but were dispersed in the cytoplasm of the PEX22-like deletion mutant PEX22-like (Figure 2). These results indicated FgPEX22-like was necessary for the location of FgPEX4 and that it functioned as a rivet of FgPEX4.

FgPEX22-Like Was Involved in Sexual and Asexual Reproduction
To elucidate the biological functions of FgPEX22-like in F. graminearum, we generated the single-gene knockout mutant PEX22-like and the double-gene knockout mutant PEX4/22-like by split-marker polymerase chain reaction (PCR) and protoplast transformation. A schematic diagram depicting the strategy used to create the FgPEX22-like gene deletion mutant and molecular analysis of PEX22-like is shown in Supplementary Figure 3. The pYF11-PEX22-like complementation construct was generated by PCR amplification using the primer pair 22hfF/22hfR (Supplementary Table 1), followed by transformation of PEX22-like mutant protoplasts to produce the PEX22-like complementation strain PEX22like-C. No obvious vegetative growth or colony morphology defects were observed for PEX22-like compared to that of the WT strain. The FgPEX4 and FgPEX22-like double mutants exhibited a significant reduction in growth rate compared to that of the FgPEX22-like single mutant and produced compact colonies with limited aerial hyphae (Figures 3A-D). These results suggested FgPEX22-like was not involved alone in the regulation of hyphal growth of F. graminearum.
When assayed for conidiation in carboxymethyl cellulose (CMC) medium cultures, PEX22-like and PEX4/22-like were reduced by 25.0 and 82.3%, respectively, compared with that FIGURE 2 | Subcellular localization of FgPEX4 in F. graminearum. Mycelium of all the strains were collected from 5-day-old complete medium (CM) plates and observed using confocal fluorescence microscopy. observe, but FgPEX4-GFP in PEX22-like exhibited a completely cytoplasmic distribution. Bar = 5 mm.  of the WT strain (Table 1). In addition, conidia of the mutants exhibited normal conidium morphology, displaying conidia that lacked intracellular content and devoid of nuclei ( Figure 4A). Phialides of the PEX22-like and PEX4/22-like mutants were rarely clustered together, unlike those of PH-1 (Figure 4B), which may have been directly responsible for reduced conidiation in the PEX22-like mutant. When incubated in yeast extract peptone dextrose (YEPD) medium, the conidium germination of PEX22-like and PEX4/22-like at 6-h post-incubation was reduced by 26.2 and 31.2%, respectively, compared with that of PH-1 and the germ tubes were shorter than those of PH-1 ( Figure 4C and Table 1). The ability of the FgPEX22-like mutants to undergo sexual reproduction was also investigated. We found the number of perithecia for PEX22-like and PEX4/22-like were reduced by 65.0 and 56.6%, respectively, compared with that of the WT strain at 2-wk post-fertilization, and the ascocarp is no difference between them ( Figure 4D). Together, these results demonstrated that FgPEX22-like played vital roles in conidiation, conidial germination, and sexual reproduction.

FgPEX22-Like Was Important for Virulence and Deoxynivalenol Production
The pathogenicity assays were performed to ascertain the effect of PEX22-like on flowering wheat heads and corn silks. Wheat heads inoculated with either PH-1 or the PEX22-like -C strain presented typical scabs, spreading from the inoculated spikelets to almost the entire head by 14 days post inoculation (dpi). In contrast, wheat heads inoculated with the PEX22like mutant presented scabs mainly on or near the inoculated spikelets ( Figure 5A). The disease index of PEX22-like and PEX4/22-like exhibited approximately a 71.9 and 74.0% reduction, respectively, compared with that of PH-1 (Figure 5B). In corn silk infection, the extended length of brown lesions caused by PH-1 and the PEX22-like-C strain were longer than those of PEX22-like and PEX4/22-like. The ability of PEX22-like or PEX4/22-like to infect corn silk was only 1/5 that of the WT strain at 5 dpi (Figures 5C,D).
We also assayed PEX4 and PEX4/22-like for DON production. We found that infected rice seeds inoculated with the PEX22-like mutant produced significantly less DON than that of PH-1 or the PEX22-like-C strain ( Table 1). These results indicated FgPEX22-like had an important role in pathogenicity and DON production.

FgPEX22-Like Deficiency Caused Abnormal Organelle Development
To determine the effect of FgPEX22-like on F. graminearum organelles, the ultrastructure of PH-1 and PEX22-like was evaluated using TEM. Spherical peroxisomes were observed in the periphery of mycelial cells of PH-1, but not PEX22-like. In addition, Woronin bodies, which were present around the mycelial septum in the WT strain, were absent in the PEX22like mutant ( Figure 6A). Besides, Compared with WT strain, PEX22-like accumulated more lipid droplets. And Nile Red staining revealed that the WT strain degraded most of the lipid droplets, while mutant strains still contained numerous bright lipid droplets (Figure 6B).
To confirm the existence of peroxisome and Woronin body, GFP fusion constructs with the Woronin body protein HEX1 or PMP70 were transfected into the WT and mutant strains and evaluated using laser-scanning confocal microscopy. For PMP70-GFP, PEX22-like displayed a green punctate distribution in hyphal cells similar to that in the WT strain (Figure 7). However, the fusion of GFP with the Woronin body protein HEX1 resulted in a punctate distribution in the WT strain while being dispersed in the cytoplasm of PEX22-like (Figure 7). These results indicated FgPEX22-like was indispensable for maintaining Woronin bodies, but not peroxisomes.
In addition, based on the TEM observations, the mutant accumulated more lipid drops, indicating the utilization rate of cellular lipid drops in PEX22-like was slower than that in PH-1. Further verification was carried out by performing Nile red staining. Brighter and larger lipid droplets were observed in PEX22-like compared to that of the WT strain. This was consistent with the TEM findings.

Deletion of FgPEX22-Like Altered the Cell Wall Integrity of Fusarium graminearum
To investigate whether FgPEX22-like participated in environmental stress responses, we examined the sensitivity of mutants to the cell wall-damaging agent Congo red and the cell membrane-damaging agent sodium dodecyl sulfate (SDS). We found that PEX22-like exhibited increased sensitivity to both SDS and Congo red compared to that of PH-1 (Supplementary  Figures 4A,B). To verify these results, hyphae were treated with lysozyme and driselase. After incubation at 30 • C for 30 min, hyphae of the PEX22-like mutant were almost completely digested and had released abundant numbers of protoplasts, whereas few protoplasts were observed among hyphae of the PH-1 and PEX22-like -C strains (Supplementary Figure 4C). Together, these results illustrated that deletion of FgPEX22-like resulted in reduced cell wall integrity in F. graminearum.

FgPEX22-Like Mutant Was More Sensitive to ROS and Involved in Lipid Metabolism
To investigate whether FgPEX22-like was involved in the response to oxidative stress, the tolerance of PEX22-like to ROS was measured. The radial growth of PEX22-like and PEX4/22-like were inhibited by 60.0 and 56.8%, respectively, when exposed to 20 mM H 2 O 2 , which was greater than that of PH-1 at 42.7% (Supplementary Figures 5A,B). Cellular ROS production in the hyphae was qualitatively analyzed by staining with nitroblue tetrazolium (NBT). The assay showed that the staining of PEX22-like and PEX4/22-like hyphae were darker than those of PH-1 (Supplementary Figure 5C). Taken together, these results demonstrated the capacity of PEX22-like and PEX4/22-like to eliminate ROS was decreased. To determine the effect of deleting FgPEX4 on fatty acid utilization, vegetative growth was assessed using minimal medium with different fatty acids as sole carbon sources. The fatty acids tested included long-chain fatty acids myristic acid (C14), palmitic acid (C16), and oleic acid (C18), and very longchain fatty acid erucic acid (C22). After incubation with the particular fatty acids for 3 days, the radial growth of the mutants was significantly reduced on the media containing the long-chain fatty acids and very long-chain fatty acids as the sole carbon sources (Figures 8A,B). These results indicated that deletion of FgPEX22-like in F. graminearum resulted in a defect in fatty acid metabolism.

DISCUSSION
Previous studies have analyzed the whole-genome sequences of 17 fungal species and found that PEX22 is less conserved than that of most peroxins (Kiel et al., 2006). In the current study, we identified a new peroxin in F. graminearum, FgPEX22-like, which has a PEX4-binding site and can directly interact with FgPEX4. Deletion of FgPEX22-like resulted in mislocation of FgPEX4. In addition, we found that FgPEX22-like was involved in the regulation of development, carbon source utilization, cell wall integrity, and pathogenicity of F. graminearum. To our knowledge, this is the first report of a PEX22 homolog playing an important role in F. graminearum.
A previous study found that yeast pex22 mutant can be fully complemented by FAM1, which is a functional ortholog of PEX22 in C. orbiculare (Kubo et al., 2015). The peroxin described in this study shared 14.34% identity with C. orbiculare and named FgPEX22-like. Although there is low sequence similarity between FgPEX22-like, FAM1 and ScPEX22, FgPEX22like is structurally similar to other PEX22 proteins with a single predicted TM domain near the N terminus and a possible PEX4 binding interface near the C terminus. In S. cerevisiae and Arabidopsis, PEX4 and PEX22 can interact with each other directly (Zolman et al., 2005;Williams et al., 2012). To determine whether FgPEX22-like functioned similar to that of PEX22 in S. cerevisiae and Arabidopsis, a yeast two-hybrid system and immunoprecipitation were employed. The results showed that FgPEX22-like could interact directly with FgPEX4 in F. graminearum.
To explore the relationship between FgPEX22-like and FgPEX4, the subcellular localization of FgPEX4 in the WT and FgPEX22-like gene deletion strain was evaluated. We FIGURE 7 | Subcellular localization of FgPMP70 and FgHEX1 in F. graminearum. Subcellular localization of FgPMP70 and HEX1. Mycelium of all the strains were harvested from 5-day-old complete medium (CM) plates and detected using confocal fluorescence microscopy. GFP-PMP70 was predominantly present in a punctate pattern (peroxisomal distribution) in the wild-type and PEX22-like. Bar = 5 mm. GFP-HEX1 was bimodally distributed in a punctate pattern (peroxisomal distribution) in the wild-type, but GFP-HEX1 in PEX22-like was completely cytoplasmically distributed. Bar = 10 mm.
found FgPEX4 in the FgPEX22-like deletion strain mislocated, which suggested that FgPEX22-like anchored FgPEX4 on the peroxisomal membrane and acted as a rivet protein for PEX4. These results were consistent with PEX22 in Pichia pastoris in which PEX4 is unstable in a PEX22 strain (Koller et al., 1999). The results illustrated that FgPEX22-like is the homologous protein of PEX22 in F. graminearum and is an important peroxin at the peroxisomal membrane that recruits and holds PEX4 at this location.
In S. cerevisiae, PEX22 plays an important role in many developmental processes (Negoro et al., 2018). In this study, we evaluated FgPEX22-like, the FgPEX22-like mutant PEX22-like, and the FgPEX4 double-knockout mutant PEX4/22-like. The results showed that FgPEX22-like on its own lacked any notable phenotype regarding hyphal growth, but it enhanced PEX4 mutant defects. This suggested that FgPEX4 and FgPEX22-like interact with each other at the molecular level and FgPEX22-like contributes to peroxisome function.
Taken together, we can conclude that FgPEX22-like contributes to the development of F. graminearum.
Peroxisome functions are important factors in plant infections. In C. orbiculare, FAM1, which encodes a woroninbody-associated PEX22, plays important roles in appressorium development and pathogenicity (Kubo et al., 2015). In this study, the PEX22-like mutant showed a significant reduction in disease severity. In addition, similar to FgPEX1, FgPEX2, FgPEX3, FgPEX4, FgPEX10, FgPEX12, FgPEX13, FgPEX14, and FgPEX33, we found that toxin production, which is known to be essential for the virulence of F. graminearum (Desjardins et al., 1993;Chen et al., 2018;Kong et al., 2019;Zhang et al., 2019a,b;Wang et al., 2020), was obviously reduced in the FgPEX22 deletion strain. Similarly, most peroxins in M. oryzae, except for PEX5, have critical roles in virulence (Wang et al., 2007Goh et al., 2011;Chen et al., 2016;Li et al., 2017). In Alternaria alternata and Colletotrichum gloeosporioides, peroxisomes are also important for the biosynthesis of AK-toxin and plant infection (Imazaki et al., 2010;Zhao et al., 2020). These results indicate that different peroxins have different regulatory mechanisms in regulating pathogenicity, and FgPEX22-like is an important pathogenic factor of F. graminearum.
Interestingly, the results of electron microscopy showed the microstructure of PEX22-like displayed an abnormal morphology in the following four aspects. First, the peroxisome structures were absent in the FgPEX22-like deletion mutant. In humans, mutations in the AAA-complex PEX26, PEX1, and PEX6 result in a decrease in peroxisome number and function (Law et al., 2017). To verify the presence or absence of peroxisome structures in the FgPEX22-like deletion mutant, we determined the localization in the mycelia of a known peroxisomal membrane protein, PMP70. The results showed that PMP70-GFP in the PEX22-like mutant displayed a punctate distribution in hyphal cells, similar to that in PH-1, indicating that the deletion of FgPEX22-like had no effect on the existence of peroxisomes. Second, Woronin bodies were also absent in the FgPEX22like deletion mutant. To verify this result, we determined the localization of the Woronin body protein HEX1 in the mycelia. The results showed that HEX1-GFP was dispersed uniformly in the cytoplasm of PEX22-like. This phenotype was in agreement with the functions of FgPEX1, FgPEX2, and FgPEX10 in the biogenesis of Woronin bodies. In addition, MoPEX11, MoPEX19, and MoPEX14/17 are also essential in M. oryzae for the biogenesis of Woronin bodies (Li et al., 2014(Li et al., , 2017. Woronin body, which is derived from peroxisomes, is a characteristic organelle specifically present in filamentous ascomycetes. Woronin bodies can plug septa into intact growing hyphae to maintain hyphal heterogeneity in a fungal mycelium by impeding cytoplasmic continuity (Liu et al., 2008).
Third, deletion of FgPEX22-like resulted in reduced cell wall integrity, displaying considerable leakage of hyphae and conidia. This result showed that the absence of FgPEX22-like influenced F. graminearum cell wall integrity. This finding was further confirmed by the sensitivity measurement of PH-1 and PEX22 to cell wall-damaging agents and degrading enzymes. Similar results were obtained for FgPEX2, FgPEX4, MoEX5, MoPEX6, MoPEX14, MoPEX17, and MoEX19, indicating they too are essential for cell wall integrity in F. graminearum and M. oryzae, accordingly (Li et al., 2014(Li et al., , 2017Zhang et al., 2019b;Wang et al., 2020). Studies have shown that the fungal cell wall glucan and chitin are derived from acetyl-CoA, which is a product of β-oxidation in peroxisomes (Ramospamplona and Naqvi, 2006). Thus, we conjectured that the damaged cell wall was related to decreased acetyl-CoA. In addition, damaged cell walls of spores and mycelia may have confirmed the functional disorders in Woronin bodies, which serve as a plug to impede cytoplasmic continuity when the mycelium is damaged to maintain hyphal heterogeneity (Tenney et al., 2000). Taken together, FgPEX22-like plays a crucial role in regulating F. graminearum cell wall integrity.
Finally, based on TEM observations, lipid droplets accumulated and increased in both quantity and size in the PEX22-like mutant. Nile red staining further verified this result. Previous studies have emphasized that peroxisomes maintain a close association with lipid bodies, which constitute the intracellular storage sites of triacylglycerol and cholesterol ester (Binns et al., 2006;Farese and Walther, 2009;Beller et al., 2010). In M. oryzae, F. graminearum, and Aspergillus flavus, the translocation and degradation of lipid droplets are also damaged by the deletion of MoPEX1, FgPEX1, FgPEX2, and AflPex5 (Deng et al., 2016;Zhang et al., 2018Zhang et al., , 2019aWang et al., 2020). The deletion of PEX6 in Colletotrichum lagenarium abrogates its ability to use long-chain fatty acids (Matsuzono et al., 1999). When cultured on medium using long-chain fatty acids or very long-chain fatty acids as sole carbon sources, the PEX22-like mutant exhibited slower growth rates, indicating FgPEX22-like was involved in utilizing long-chain and very long-chain fatty acids. Taken together, these results indicate that deletion of FgPEX22-like results in a deficiency in the utilization of lipids and long-chain fatty acids. Combining these results, we suggest that FgPEX22-like plays essential roles in maintaining normal organelle development in F. graminearum, which probably accounts for the loss of pathogenicity in the FgPEX22-like mutant.
Herein, we show that FgPEX22-like encodes a peroxin protein in F. graminearum that interacts directly with FgPEX4 and acts as a rivet protein of FgPEX4, and the loss of FgPEX22like leads to the abnormal subcellular localization of FgPEX4 protein. Moreover, FgPEX22-like is involved in the regulation of asexual and sexual reproduction, pathogenicity, cell wall integrity, oxidative stress, and organelle integrity. Our study has established, for the first time, the comprehensive biological functions of a homologous protein of PEX22 in F. graminearum.

DATA AVAILABILITY STATEMENT
The data that supports the findings of this study are available in the Supplementary Material of this article.

AUTHOR CONTRIBUTIONS
JY: writing-review and editing. LZ and CL: investigation and writing-original draft preparation. MW: data curation. YT: software. YL: resources. All authors contributed to manuscript revision, read, and approved the submitted version.