Ethylene Promotes Expression of the Appressorium- and Pathogenicity-Related Genes via GPCR- and MAPK-Dependent Manners in Colletotrichum gloeosporioides

Ethylene (ET) represents a signal that can be sensed by plant pathogenic fungi to accelerate their spore germination and subsequent infection. However, the molecular mechanisms of responses to ET in fungi remain largely unclear. In this study, Colletotrichum gloeosporioides was investigated via transcriptomic analysis to reveal the genes that account for the ET-regulated fungal development and virulence. The results showed that ET promoted genes encoding for fungal melanin biosynthesis enzymes, extracellular hydrolases, and appressorium-associated structure proteins at 4 h after treatment. When the germination lasted until 24 h, ET induced multiple appressoria from every single spore, but downregulated most of the genes. Loss of selected ET responsive genes encoding for scytalone dehydratase (CgSCD1) and cerato-platanin virulence protein (CgCP1) were unable to alter ET sensitivity of C. gloeosporioides in vitro but attenuated the influence of ET on pathogenicity. Knockout of the G-protein-coupled receptors CgGPCR3-1/2 and the MAPK signaling pathway components CgMK1 and CgSte11 resulted in reduced ET sensitivity. Taken together, this study in C. gloeosporioides reports that ET can cause transcription changes in a large set of genes, which are mainly responsible for appressorium development and virulence expression, and these processes are dependent on the GPCR and MAPK pathways.


Introduction
Fresh fruits are valuable dietary sources for human health; however, they are also attractive targets for phytopathogenic fungi, leading to fruit decay and economic losses. As the fruits ripen and senescence, the physicochemical properties of the host tissues become increasingly susceptible to microorganism infections. On the other side, the pathogens should be capable of recognizing fruit ripening status to orchestrate their infection activities, although the mechanisms by which they do this remain unclear.
The gaseous phytohormone ethylene (ET) is involved in many aspects of plant growth and development, including senescence, organ abscission, and fruit ripening [1][2][3]. In plants, ET is perceived by transmembrane-receptor proteins localized in the endoplasmic reticulum (ER) belonging to the ET receptor family that harbors homology to bacterial two-component histidine kinases [4,5]. The ET signal transduction pathway in plants ultimately leads to regulating the expression of ET-dependent genes [6,7]. In flesh fruits, ET regulates ripening by coordinating the expression of genes that participate in multiple biological processes, such as the increase in respiration, autocatalytic ET production, pigment accumulation, change in texture, and the buildup of overall fruit quality traits [8]. As such, ET is known as the fruit-ripening hormone. Because of its gaseous nature, ET can diffuse rapidly throughout plant materials, cross cell membranes, and escape the plant as an environmental signal, which could play important roles in interactions between host plants and pathogenic microbes [9].
In contrast to healthy fruits, the diseased ones usually release notably increased ET, which can thereby stimulate the ripening and rotting of nearby fruits, fulfilling the probable capacity of ET to accelerate fruit spoilage in holistic scales [10,11]. ET is also involved in a variety of plant defense responses to microbial attacks [12,13]. A large number of studies have shown that, depending on the nature of the plant-microbe interaction, ET may reduce, enhance, or have no effect on disease development [14]. Most studies have focused on the roles of ET signaling in the immunity of host plants, leading to the commonly accepted view that ET and jasmonic acid signaling components are involved in resistance against necrotrophic pathogens [15,16]. However, the effects of ET on the pathogen side are relatively less understood.
It is noteworthy that plant pathogenic fungi are capable of sensing ET. For example, ET can promote the spore germination and appressorium formation of Colletotrichum spp. [17]; In contrast, 100 ppm and 1000 ppm exogenous ET released by ethephon mildly inhibited the spore germination of Botrytis spp., while low concentrations of ET showed a negligible effect [18]. Treatment of B. cinerea with 2,5-norbornadiene, a kind of ET receptor competition inhibitor, can inhibit the germination of spores, while ET can eliminate the inhibitory effect [19]. Similar phenomena have been found in Alternaria alternata [20]. Additionally, the expression of the pathogenic gene bcspl1 of B. cinerea could be induced by ET [21]. In Aspergillus flavus, ET can suppress aflatoxin production, and inhibit the expression of aflatoxin-synthetic genes (aflR and aflD) [20,22]. However, the effect of ET on the overall gene expression profiles of plant pathogenic fungi has not been reported yet, leaving a knowledge gap in our understanding of the mechanisms underlying the effects of ET on fungal growth, development, and pathogenicity.
Colletotrichum gloeosporioides is one of the most common pathogens causing postharvest diseases in varied fruits and vegetables [23][24][25]. C. gloeosporioides is a characteristically hemibiotrophic fungus that can survive in the plant cuticle or epidermal cells until the fruits ripen [26]. In other words, C. gloeosporioides prefers to remain under the latent infection stage in immature fruits and vegetables, with the anthracnose symptoms occurring as the host matures. The conidia of C. gloeosporioides can form a typical dome-shaped infection structure called appressorium [27]. The fungal spores need to spend their substance storage inside to provide energy for the formation of appressoria, and then rely on the newly developed cell structures to infect plants and obtain nutrients to maintain the growth and reproduction [28,29]. C. gloeosporioides represents a typical plant pathogenic fungus showing obvious ET-sensitive responses, which include spore germination and appressorium formation [17]. Therefore, ET may be recognized as a signal of fruit ripening by this pathogenic fungus to initiate the occurrence and development of anthracnose disease.
The G-protein-coupled receptors (GPCRs) are the largest class of cell surface receptors in eukaryotes, sensing sense environmental cues to initiate intracellular G-protein signaling that coordinates biological responses [30,31]. In Magnaporthe grisea, the G protein subunits, MagB and Mgb1, are required in a cAMP-dependent manner for recognition of host surface hydrophobicity during appressorium formation [32,33]. However, cAMP does not rescue defects in appressorial penetration and infectious growth in Mgb1 and Mgg1 mutation strains [33,34]. In addition, the putative GPCR Pth11 also regulates appressorium formation, presumably through sensing surface hydrophobicity and plant cutin monomers [35]. Thus, Pth11/MagB/Mgb1/Mgg1 could comprise a GPCR/G protein system that controls appressorium differentiation in response to molecular cues from the host plant. Interestingly, a Gα null mutant (∆bcg1) in B. cinerea was ET insensitive [21], and one identified GPCR (BcGPR3) demonstrated the ability to bind ET [36]. Cumulatively, these studies imply that certain GPCR/G protein systems could be involved in ET sensing in fungal pathogens. Upon perceiving exogenous stimuli by the GPCR/G protein system, the signals could be further relayed to intracellular cyclic adenosine monophosphate (cAMP) or mitogenactivated protein kinase (MAPK) signaling pathways, vital regulators for the establishment of appressorial structure and function in plant pathogenic fungi [37][38][39]. The studies with C. gloeosporioides revealed that deletion of the gene encoding the catalytic subunit of cAMPdependent protein kinase A retarded the formation of appressorium [39]. Meanwhile, C. gloeosporioides containing mutated mitogen-activated protein kinase (CgMK1) could not form appressorium. Exogenous addition of cAMP failed to restore the mutant's defect, indicating that CgMK1 may function downstream or independent of cAMP-mediated signaling for appressorium formation [40]. Thus, whether the intracellular secondary signaling pathways are involved in ET responses to regulate appressorium development is worth investigating.
The present study aimed to uncover the molecular events regulated by ET to affect the development and pathogenicity of C. gloeosporioides. Through transcriptome analysis, the ET-regulated genes related to appressorium formation and pathogenicity of C. gloeosporioides were systematically explored. In addition, reverse genetic study determined the functions of specific ET signaling and responsive genes, leading to the speculation that ET may promote appressorium formation and pathogenicity of C. gloeosporioides through the GPCRs and MAPK signaling pathways.

Fungal Pathogen and Host Plant Materials
Single-spore strain EX2016-02 of C. gloeosporioides has been isolated and characterized in a previous study [41], and this strain was used as a recipient for genetic modification. All the fungal strains used in this study are listed in Table 1. To generate fungal conidial suspension, approximately two-week-old, sporulating cultures on potato dextrose agar (PDA) plates were flooded with sterile water. The resultant conidial suspension was measured for concentration and adjusted to 3 × 10 5 conidia mL −1 using a hemocytometer. Arabidopsis thaliana (Col-0) was grown in vermiculite soil at 23 • C under short-day conditions (8 h light/16 h dark) for 4 weeks to obtain detached leaves for the inoculation assay.

Ethylene Treatment and Sample Preparation
Gaseous ethylene (ET) was purchased from Pujiang Special Gas Co., Ltd., Shanghai, China. To examine the effects of ET, the samples were placed in humidified plastic boxes with tight sealing caps. Pure ET was injected into the boxes to achieve specified concentrations, which were confirmed via gas chromatographic analysis according to the reported method [43]. Meanwhile, the samples placed in the boxes without ET input were used as controls.
To harvest fungal samples for transcriptome study, conidia suspension was poured into 90 mm diameter glass Petri dishes, with 20 mL per dish. The Petri dishes were placed in sealed plastic boxes. Control and ET (200 ppm) treatment were performed as shown above. After incubation for 4 or 24 h, the supernatant liquid was gently poured out, and the remaining conidia and germlings were scratched from the bottom of the dishes by using a cell scraper. Three independent samples were collected for each treatment, and the experiment was repeated three times. The resulted samples were immediately frozen in liquid nitrogen, followed by thorough lyophilization before RNA extraction.

Total RNA Extraction
Total RNA was extracted using QIAGEN Reagent (Hilden, Germany) and then treated with DNase I according to the manufacturer's instructions. The quality of extracted total RNA was verified using a Nanodrop 2000 (Thermo Fisher Scientific Inc., Waltham, MA, USA) and Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). The total amount of RNA required for a single sample of transcriptome test was not less than 5 µg, and the optical density value of 260 over 280 was between 1.8 and 2.2.

Library Preparation and RNA Sequencing
The RNA library was constructed using the Illumina Truseq TM RNA sample prep kit method (Origingene biotechnology company, Shanghai, China). Briefly, three independent samples were collected for each treatment. Total RNA was extracted and enriched by oligo-dT selection. The mRNA samples were fragmented and reverse-transcribed to doublestrand cDNA using the random primers and reverse transcriptase. The sticky ends of synthesized cDNA were repaired to blunt ends using the End Repair Mix, and then to the 3' end, an A base was added to connect to the Y-shaped adaptor. Paired-end reads of 150 bases were obtained using Illumina HiSeq 2000 and converted into sequence data through Base Calling.

Transcriptomic Data Analysis
The clean reads were filtered using cutadapt (v1.16) from the low-quality reads, which were stored in FASTQ format [44]. Genome assemblies of C. gloeosporioides Cg-14 (NCBI accession number SUB133583) were used as a reference sequence, and the reads were aligned to the reference sequence with hisat2 software (v2.1.0) [45]. Fragments per kilobase per million (FPKM) for differential expression using StringTie software (v1.3.3b) were calculated [46]. Differentially expressed genes (DEGs) (fold change ≥ 2, adjusted p-value ≤ 0.05) were identified with the DEseq2 algorithm [47]. Heat map visualization was performed using ggplot2 and pheatmap packages, and the overall number of differential genes were counted using the UpSetR package [48,49]. Principal component analysis (PCA) was calculated from the expression data using ggpubr, ggthemes, and gmodels packages. Blast2GO software was used for functional annotation of the transcripts [50]. GO analysis and visualization were performed using the R language GOcircle package.

Quantitative Real-Time Reverse Transcription PCR (qRT-PCR) Analysis
For qRT-PCR analysis, 1 µg of each RNA sample was used for reverse transcription with the Prime Script™ RT reagent kit (Perfect Real Time) (TakaRa Biotechnology, Co., Dalian, China). The real-time PCR amplifications were conducted in a CFX96TM Real-Time System (BIO-BAD, Inc., California, USA) using TakaRa SYBR Premix Ex Taq (TakaRa Biotechnology, Co., Dalian, China). For each sample, the expression of the Actin gene was used as an internal reference. The primers for the qRT-PCR assays are listed in Table S1. There were three replicates for each sample, and the experiment was repeated three times.
Relative expression levels of the tested genes were calculated using the 2 −∆∆Ct method [51].

Fungal Transformation
Molecular manipulation and PEG-mediated protoplast transformation as described in a previous study were used for the purposes of deleting target genes in C. gloeosporioides [42]. For the purposes of expressing exogenous genes, Agrobacterium tumefaciens-mediated transformation (ATMT) was used. The recombinant plasmids based on the backbone pAg1-H3 vector were constructed and transformed into A. tumefaciens strain Agl1 for ATMT according to the reported method [52].

Fluorescence-Based Reporter Assay
The fluorescence-based reporter assay was conducted to compare the degree of ET responses in C. gloeosporioides. Briefly, the promoters of relevant ET-responsive genes were fused with GFP and assembled into the pAg1-H3 vector at the SacI digestion site via a one-step cloning kit (Yeasen, Shanghai, China). The resulted recombinant plasmids were transformed into the fungus by the ATMT method. The conidial suspension (3 × 10 5 conidia mL −1 ) of the transformants was dropped on sterilized glass slides and incubated in humidified plastic boxes for 24 h. Control and ET treatments were conducted as described above. Observation of the fluorescent expression of transformants with a fluorescence microscope (Axio Imager 2, ZEISS, Jena, Germany). For each treatment, three replicates of 200 conidia were randomly selected for observation recording.

Appressorial Formation and Pathogenicity Assay
To examine the effects of ET on conidial germination and subsequent hyphal development, a conidial suspension (3 × 10 5 conidia mL −1 ) was dropped on sterilized glass slides and incubated in humidified plastic boxes. Control and ET treatments were conducted as described above. After incubation for 4 or 24 h, conidial germination and appressorium formation rates were recorded. For each treatment, three replicates of 300 conidia were randomly selected for statistical analysis.
In the pathogenicity assay, A. thaliana leaves were detached and inoculated with 6 µL conidia suspension. The inoculated materials were kept in humidified boxes containing either 0 or 200 ppm ET at the beginning of the assay. The boxes were placed under light and incubated at 25 • C. After 24 h, the samples were subjected to DAB staining analysis.

Staining Analysis for Reactive Oxygen Species in Plant Leaves
Reactive oxygen species (ROS) accumulation upon pathogen infection in plant leaves were visualized by DAB (3,3 -diaminobenzidine) staining. Briefly, 220 mg DAB was added to 200 mL sterile water. The pH of the solution was adjusted to 3.7 to completely dissolve the contents. The solution was mixed with 110 mL Tween 20 (0.05% v/v) and 0.8 g Na 2 HPO 4 to obtain the DAB dye solution. The samples for analysis were stained in the freshly prepared solution by shaking at 100 rpm in the dark for 4 h. After staining, the samples were boiled in 95% ethanol to decolorize. Photographs of the samples were analyzed with ImageJ software to calculate the depth of staining according to the reported method [53].

FDA-PI Dual Fluorescence Staining
FDA (fluorescein diacetate) was dissolved in acetone to obtain a final concentration of 30 mg mL −1 , and PI (propidium iodide) was dissolved in DPBS buffer (Dulbecco's phosphate-buffered saline) to obtain a final concentration of 10 mg mL −1 . Five microliters of the above solutions were diluted to 1 mL with sterile water to obtain a fresh staining solution. The conidial suspension was placed on the glass slide to germinate for 4 h or 24 h, and the staining solution was added to the samples. After staining at room temperature for 30 min, the dyes were drained, and the samples were washed 3 times with sterile water. Images were acquired with the fluorescence microscopic analysis as follows: live cells were stained with FDA to show a bright green fluorescent signal (530 nm) upon excitation light (488 nm), while in contrast, dead cells were stained with PI to emit a red fluorescent signal (617 nm) under green excitation light (535 nm).

Calcofluor White-EosinY Double Fluorescence Staining
The double staining assay was conducted according to the previous method [54], in which chitin and its deacetylated form chitosan were stained by Calcofluor white (18909-100ML-F, Sigma-Aldrich, St. Louis, Missouri, USA) and EosinY (A600190-0025, Sangon Biotech, Shanghai, China), respectively. The conidial suspension was placed on the glass slide for germination. After staining in the dark at room temperature for 10 min, the staining solution was removed, and the samples were analyzed with a fluorescent microscope. Fluorescent signals of calcofluor white and Eosin Y were captured using 405 and 561 nm excitation wavelengths, respectively.

Statistical Analysis
Statistical analyses were carried out with the GraphPad Prism 8.0 software. Bar charts represent mean values with standard deviations; differences were considered significant at p < 0.05 (*).

Ethylene Accelerated Conidial Germination, Appressorium Formation, and Pathogenicity of C. gloeosporioides
The conidia of C. gloeosporioides treated with ET demonstrated an appressorium formation rate of 71.96% on glass slides at 4 h post-inoculation (hpi), in contrast to that of the control (CK) value of 10.64%. At 24 hpi, the gap between the CK (86.08%) and ET treatment (91.39%) groups in appressorium formation rates was not comparable as it had been at 4 hpi; however, the appressorium yield per single conidium was higher in the ET group than that of CK group at 24 hpi ( Figure 1A,B). FDA and PI double staining showed that germinated spores in the CK-and ET-treated groups survived at 4 hpi, as these cells were filled with FDA staining signal (bright green). However, at 24 hpi, the proportion of PI-stained cell areas was higher in the ET-treated samples, in other words, the prolonged ET treatment resulted in stimulating appressorium development, which was accompanied by apoptotic processes ( Figure S1). Fluorescent microscopic analysis with the Cg-gfp strain demonstrated that ET could also promote the appressorium formation on the Arabidopsis leaf surface. Additionally, the accumulation of reactive oxygen species (ROS) at the inoculation site as visualized by DAB staining showed that there were stronger ROS in the ET group than in the CK ( Figure 1C). These results indicate that ET can significantly promote the appressorium differentiation and host infection of C. gloeosporioides.

ET Treatment Altered Transcript Levels of a Large Set of Genes in Germinating Conidia of C. gloeosporioides
In order to reveal the molecular events related to the effects of ET on C. gloeosporioides, RNA-seq assays were conducted with the germinating conidia at 4 and 24 hpi. Compared with the CK group, ET treatment resulted in 833 upregulated and 997 downregulated genes at 4 hpi. However, ET treatment only upregulated 78 genes but downregulated 2649 genes at 24 hpi (Figure 2A). The up-and downregulated genes by ET treatment at both 4 hpi and 24 hpi were 9 and 518, respectively ( Figure 2B). PCA revealed that there were significant differences in gene expression between the CK-4hpi VS CK-24hpi, ET-4hpi VS CK-4hpi, and ET-24hpi VS CK-24hpi groups. However, the ET-4hpi and CK-24hpi groups were overlapped in PCA (Table S2, Figure 2C). The gene expression values obtained from RNA-seq were validated using the qRT-PCR assay. The expression levels of 27 randomly selected C. gloeosporioides genes demonstrated a strong correlation (R 2 = 0.9089) between the results obtained with the two techniques ( Figures 2D, S2 and S3, Table S3).

Functional Enrichment Analysis of Differentially Expressed Genes
The differentially expressed genes (DEGs) as influenced by ET treatment were sub-  Table S4).
As the "extracellular region" GO enrichment included the largest number of upregulated DEGs at 4 hpi, it was further found that these DEGs mainly encode two types of products: extracellular hydrolases and extracellular structural components. As shown in Figure S4, ET treatment can significantly promote the expression of the genes encoding for polysaccharide deacetylase, cutinase, multicopper oxidase, and cerato-platanin protein. At the same time, ET could also promote transcript levels of the genes encoding for cell wall proteins, appressoria structure proteins, glycolipid-anchored surface protein, and other structural proteins.

Hydrophobic Surface Binding Protein A and Cutinase-Related Genes Were Induced by ET Treatment at the Early Germinating Stage
Twelve most strongly induced DEGs by ET at 4 hpi are described in Table 2. Via Protein Blast analysis to characterize the conserved domains encoded by these DEGs, three genes (CGLO_03844, CGLO_00547, CGLO_13252) encoding for hydrophobic surface binding protein A domains (HsbA) were identified. Via the GFP-based reporter analysis, it was shown that ET treatment promoted the GFP-expressing signal when this reporter gene was driven by the promoters of HsbA genes (CGLO_03844 and CGLO_13252), suggesting that these HsbA genes are indeed highly sensitive to ET stimulus ( Figure 4A,C, Table S5).
Filamentous fungal cell walls are generally composed of chitin polysaccharides. However, this cell wall component can be modified by the fungi themselves for certain purposes, such as evading the host immune recognition and thus enhancing the pathogenicity of the fungi. Chitin deacetylases (CDA) are enzymes that catalyze the transformation of chitin into chitosan. The genome of C. gloeosporioides contains seven putative CDA-encoding genes, among which five CDA genes (CGLO_07057, CGLO_16772, CGLO_00212, CGLO_16754, and CGLO_02140) were upregulated by ET. Double-staining assay with Calcofluor and Eosin Y revealed that as appressorium formation proceeded, chitin was deacetylated to form chitosan at the appressorial zones ( Figure 5B, Table S5). These results indicate that the appressorial cell wall of C. gloeosporioides is characterized by accumulating melanin and chitosan, which can be highly promoted by the presence of ET.

Pathogenicity Enhancement of C. gloeosporioides in Response to ET Was Partially Dependent on the Upregulation of Melanin Synthase and Effector Genes
Hitherto, the study confirms that ET can promote appressorium formation and infection of C. gloeosporioides. Accompanying ET-induced genes were further screened for their contribution to the ET-regulated pathogenicity. Three selected ET-responsive genes were targeted for knockout: CgCAP22 (CGLO_01483), predicted to encode the homolog of appressorium-specific protein CAP22 [58]; CgCP1 (CGLO_12973), encoding the virulent effector cerato-platanin [59]; and CgSCD1 (CGLO_09685), encoding the melanin biosynthesis enzyme [42]. The resultant mutant strains of ∆Cgcap22, ∆Cgscd1, and ∆Cgcp1 showed no significant difference in the growth rate compared to the WT ( Figure 6A).
Furthermore, in vitro assays demonstrated that ET induced appressorium development in all the tested mutants by comparable level in the WT, suggesting that the mutants retain sensitivity to ET ( Figure 6B). In contrast, in the infection assay on Arabidopsis leaves, ET treatment enhanced ROS accumulation at the inoculation sites of the WT and ∆Cgcap22 mutant spores, while the ∆Cgscd1 and ∆Cgcp1 inoculation sites resulted in negligible ROS accumulation, irrespective of ET treatment or not ( Figure 6C). Consequently, upregulating the melanin synthesis and the effector genes is likely of vital importance for ET to promote the pathogenicity of C. gloeosporioides.
Concerning the ET-responsive phenotypic characters, the ∆Cgmk1 and ∆Cgste11 mutants showed defects in appressorium formation, but spore germination was slightly promoted by ET ( Figure 7D). Moreover, the ∆Cggpcr3-1, ∆Cggpcr3-2, and ∆Cggpcr3-1/2 showed elevated spore germination rates in response to ET; however, induction of appressorium formation by ET was attenuated in these GPCR mutants, among which the double mutant ∆Cggpcr3-1/2 was the most extremely depressed ( Figure 7B,D,E).
DAB staining further demonstrated that upon ET treatment, the inoculation sites with WT and ∆Cggpcr3-1/2 showed stronger ROS accumulation than their corresponding control samples. In contrast, the ∆Cgste11 inoculation sites resulted in negligible ROS accumulation in both control and ET treatment conditions ( Figure 7F). Collectively, these results imply that the MAPK-and GPCR-mediated signaling pathways are likely involved in the ETaccelerated appressorium development; however, it is the MAPK pathway component that is required for mediating ET to promote the pathogenicity of C. gloeosporioides in the plant. Moreover, gene expression analysis indicated that several ET-induced genes, including CGLO_09685 (CgSCD1), CGLO_12793 (CgCP1), CGLO_16754 (encoding for chitin deacetylase), CGLO_03844 and CGLO_13252 (encoding for hydrophobic surface binding protein A), and CGLO_15181 (hypothetical protein), were extensively impaired in ET responsiveness in the ∆Cgmk1 and ∆Cgste11 mutants. As for the ∆Cggpcr3-1/2 double mutant, expression levels of the above genes except for CGLO_12793 (CgCP1) were also not sensitive to ET stimulus (Figure 8).

Discussion
Ethylene (ET) is the ripening hormone for fruit crops to achieve full maturation and attractive qualities. Microbial contamination and rotting are favored by fruit ripening, and usually lead to the release of more ET, thereby stimulating the ripening and rotting of nearby fruits. In this context, the influence of ET on both fruits and their rotting agent fungi should be precisely understood to achieve better management of postharvest losses. The ET signaling mechanisms in plants have been extensively elucidated since the 1990s, mostly due to a series of milestone findings based on elegant genetic analysis with the model plant Arabidopsis [60,61], while the ET action mode and mechanisms on plant pathogens are still largely unknown. The current study used transcriptomic and reverse genetic analysis to reveal the molecular response mechanisms to ET in the anthracnose fungus C. gloeosporioides.
C. gloeosporioides is a typical hemibiotrophic fungus that infects plants by forming appressoria. Spore germination and subsequent appressorium development of this fungus can be accelerated by ET [17]. Using the transgenic strain of Cg-gfp expressing GFP, this study confirms that ET can significantly promote the formation of appressoria in vitro and on Arabidopsis leaves. Fungal ET responses such as accelerated spore germination appear similar to the postulated ancestral phenomenon of ET motivating cell elongation in plants [62,63]. Moreover, appressoria are known as specialized cells or adhesion structures from which a penetration peg emerges to pierce the host tissue [64]. We found that ET enhanced appressoria development and ROS accumulation at the inoculation site on leaves, suggesting that ET could influence the outcome of the pathogen-host interaction by affecting the fungal side. These ET-sensing responses in fungi, exemplified by C. gloeosporioides, are thought to provide the pathogen an adaptive advantage by allowing them to accurately time their infections to the host's ripening stages [17].
In order to reveal the molecular mechanisms underlying the ET responsive phenotypes of C. gloeosporioides, the transcriptome data of the germinating spores were analyzed. ET was found to cause a large number of upregulated genes at 4 hpi, with the genes coding for "extracellular region"-related genes, such as the hydrophobic-surface-binding protein A (HsbA) and cutinase-related proteins, cell wall structural proteins, melanin synthesis enzymes, chitin deacetylases, and virulence proteins being extremely enriched. These DEGs are closely related to the phenotypic responses to ET in this phytopathogenic fungus.
In filamentous fungi, cell walls are commonly darkened by accumulating DHNmelanins. These pigments are secondary metabolites produced via the polyketide synthase pathway [56]. In this study, microscopic observation of the melanized appressoria is in accordance with the gene expression data, indicating that the melanin biosynthesis activity is accelerated by ET during the spore germinating and appressorium developing stages of C. gloeosporioides. The roles of appressorial melanization in fungal pathogenicity has been well documented [65,66]. Melanin may contribute to the cross-linking of cell wall components in order to reinforce the rigidity of appressoria, allowing them to withstand hyper-osmosis and generate mechanical infection force [67]. Meanwhile, hyperaccumulation of melanin pigments at the appressorial surface may also contribute to fungal survival during the oxidative burst of host plants [68,69]. This study found that ET can promote the melanin-defective mutant ∆Cgscd1 to germinate and develop appressoria as it does in the wild-type strain; however, ET had little effect on the pathogenicity of the ∆Cgscd1 mutant. These findings suggest that melanin could be an indispensable cellular constituent being responsive to ET stimulus to determine the pathogenicity of C. gloeosporioides.
Another representative gene induced by ET is CgCP1 coding for the cerato-platanin effector, reminiscent of the scenario wherein the homologous gene Bcspl1 in B. cinerea was also sensitive to ET [21]. BcSpl1 is one of the most abundant proteins in the secretory group of B. cinerea, and the loss of this gene led to a decrease in the pathogenicity of the fungus [70]. Moreover, the cerato-platanin homologs have also been found to be required for the pathogenicity of M. grisea [71] and C. gloeosporioides [59]. This study shows that the ∆Cgcp1 mutant is still ET-sensitive in vitro; however, exogenous ET could not promote the pathogenicity of ∆Cgcp1, suggesting that cerato-platanin is a vital virulence factor being stimulated by ET to accelerate host infection.
Secreted enzymes were included in the DEGs. Genes coding for chitin deacetylases (CDAs) were upregulated by ET at 4 hpi, which was phenotypically verified by differential staining of chitin and chitosan, indicating the presence of chitosan in appressoria but not in vegetative hyphae. Chitin undergoing deacetylation into chitosan in the cell wall has been associated with hyphal growth and differentiation, such as polar growth in Aspergillus fumigatus [72], appressorium development in Magnaporthe oryzae [73] and Uromyces fabae [74], and formation of the infection structure in Botrytis cinerea [54]. The transformation of chitin into chitosan in C. gloeosporioides thus might also be responsive to ET to accelerate appressorium development and pathogenicity. The production of chitosan plays a role in cell wall integrity and anchoring of melanins [75][76][77]. Moreover, chitin deacetylation is thought to prevent the detection of chitin as a pathogen-associated molecular pattern (PAMP) by the host immune system [78]. In the powdery mildew fungus, RNAi silencing or chemical inhibition of CDA resulted in a dramatic reduction in fungal growth, which was linked to a rapid elicitation of chitin-triggered immunity; thus, CDA represents a promising target for the control of phytopathogenic fungi [79]. However, due to the redundancy of CDA genes, the exact roles of chitin deacetylation in C. gloeosporioides remain open to be characterized.
Moreover, fungal extracellular cutinases can help the pathogen penetrate through the outermost cuticular barriers of host plants. Cutin monomers, produced by low levels of constitutively expressed cutinase, are prone to inducing high levels of cutinase, which can aid pathogenic fungi in penetrating the host through the cuticle, with cutin as the major structural polymer [80]. Cutinase transcription factor 1α (CTF1α) and CTF1β are involved in strongly inducing and constitutively maintaining cutinase expression, respectively [80,81]. Similarly, the HsbA is known to recruit cutinases to the hydrophobic surface [55]. The transcriptomic analysis in this study revealed that ET can stimulate the expression levels of both CTF1α, CTF1β, and HsbAs and their related cutinases during the early germinating stage of C. gloeosporioides. These findings suggest that ET signaling is capable of triggering the invasive weapons of the fungal pathogen to pierce the host.
Programmed cell death and autophagy are also known to be involved in conidial germination and appressoria differentiation in some fungal plant pathogens [28,82]. The ET-motivated appressorium formation is a series of morphogenesis involving cell cycling processes. ET treatment leads to the development of multiple appressoria from each individual spore, accompanied by cellular apoptosis as evidenced by FAD and PI staining assay. These findings imply that ET may be recognized by the fungal pathogens as a triggering signal for ripening and susceptibility in host plants, and upon the ET stimulus, the fungal spore germlings should exert all their energies until sacrificed to establish infections.
Besides ET, many other host-plant-related materials may stimulate spore germination and appressoria differentiation of pathogenic fungi, such as adhesives, cutin monomers, surface hydrophobicity, surface hardness, and topographic signals [83,84]. Several signal pathways participate in regulating the appressorium formation, such as Gα and Gβ proteins [33], the adenylyl cyclase and PKA-mediated cAMP pathway [85], calmodulin signaling [86], and the MAPK signaling pathway [37]. It is proposed that the G-proteincoupled receptor BcGPCR3 of B. cinerea possesses the ability to bind ET and regulate responses to ET [36], while the MAPK cascade and their protein phosphorylation are involved in regulating the ET responses of C. gloeosporioides [87]. This study reveals that the MAPK pathway components (Cgmk1 and Cgste11) play core roles in mediating morphological and molecular responses to ET. However, the identified GPCR3 homologs (Cggpcr3-1/2) are responsible for partial ET responsiveness, suggesting that there are still remaining undiscovered components that could participate in ET sensing in this fungus.
Collectively, this paper reports that ET could alter transcript levels of the genes that are closely relevant to accelerating spore germination, appressoria formation, and virulence expression of C. gloeosporioides. Moreover, ET promoting the gene expression, morphogenesis, and pathogenicity of C. gloeosporioides seems to be dependent on the GPCRs and the MAPK pathway, although the upstream perception mechanisms for ET signal in the fungus still remain obscure. These findings correspond well with the hypothesis that phytopathogenic fungi could response to ET via the pathway being different from in plant [88]. Thus, it is possible to selectively block the ET perception in fungal pathogens for the purpose for disease management during fruit ripening, although achieving this goal should be based full understanding the exact ET sensing and signaling mechanisms in fungi in the future.
Supplementary Materials: The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/jof8060570/s1, Figure S1: Statistics of red fluorescence area in FDA-PI dual fluorescence staining. Figure S2: qRT-PCR verification of upregulated differential genes. Figure S3: qRT-PCR verification of upregulated differential genes. Figure S4: Annotation information of genes enriched in extracellular region (GO: 0005576). Table S1: Oligonucleotide primers used for target amplification in this study. Table S2: Principal component analysis. Table S3: Correlation analysis of QPCR and RNAseq. Table S4: GO terms of fungal transcripts that are significantly differentially expressed by ethylene treatment. Table S5: Colletotrichum gloeosporioides gene expression profile during ethylene treatment 4 hpi.