Overexpression of CmWRKY8-1–VP64 Fusion Protein Reduces Resistance in Response to Fusarium oxysporum by Modulating the Salicylic Acid Signaling Pathway in Chrysanthemum morifolium

Chrysanthemum Fusarium wilt, caused by the pathogenic fungus Fusarium oxysporum, severely reduces ornamental quality and yields. WRKY transcription factors are extensively involved in regulating disease resistance pathways in a variety of plants; however, it is unclear how members of this family regulate the defense against Fusarium wilt in chrysanthemums. In this study, we characterized the WRKY family gene CmWRKY8-1 from the chrysanthemum cultivar ‘Jinba’, which is localized to the nucleus and has no transcriptional activity. We obtained CmWRKY8-1 transgenic chrysanthemum lines overexpressing the CmWRKY8-1-VP64 fusion protein that showed less resistance to F. oxysporum. Compared to Wild Type (WT) lines, CmWRKY8-1 transgenic lines had lower endogenous salicylic acid (SA) content and expressed levels of SA-related genes. RNA-Seq analysis of the WT and CmWRKY8-1-VP64 transgenic lines revealed some differentially expressed genes (DEGs) involved in the SA signaling pathway, such as PAL, AIM1, NPR1, and EDS1. Based on Gene Ontology (GO) enrichment analysis, the SA-associated pathways were enriched. Our results showed that CmWRKY8-1-VP64 transgenic lines reduced the resistance to F. oxysporum by regulating the expression of genes related to the SA signaling pathway. This study demonstrated the role of CmWRKY8-1 in response to F. oxysporum, which provides a basis for revealing the molecular regulatory mechanism of the WRKY response to F. oxysporum infestation in chrysanthemum.


Introduction
Fusarium wilt is a severe soil-borne disease that causes yellowing and wilting of plant leaves by damaging their vascular bundles. Fusarium oxysporum is the culprit of wilt disease, which usually invades the roots and multiplies the vascular bundles. F. oxysporum can block the vascular bundle and water cannot be transported, eventually leading to plant death [1][2][3]. The plant defense against pathogens relies on the natural immune system of plants. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is triggered when PAMPs are sensed by pattern recognition receptors (PRRs) located on the surface of plant cells [4]. Intracellular nucleotide-binding leucine-rich repeat (NB-LRR) sensing specific effectors of pathogens can, in turn, trigger effector-triggered immunity (ETI) [5]. Among the transcription factors studied, WRKY transcription factors play an indispensable role in PTI and ETI immune pathways.
WRKY transcription factors are a family of transcription factors that are unique to plants and are widely distributed. WRKY transcription factors are widely involved in the regulation of plant growth, development, senescence, organ synthesis, and various hormone-mediated signaling pathways [6]. More importantly, WRKY transcription factors isolated and characterized CmWRKY8-1 and studied its function in response to F. oxysporum using transgenic technology. Finally, we observed that CmWRKY8-1 responds to F. oxysporum infestation by regulating the SA signaling pathway in chrysanthemum.

Isolation and Sequence Analysis of CmWRKY8-1
Previously, CmWRKY8 (KC615362) was found to respond positively to F. oxysporum inoculation in chrysanthemums 'Jinba', suggesting that CmWRKY8 may be involved in resistance to F. oxysporum in chrysanthemums [28]. We isolated a 744 bp full-length open reading frame (ORF) encoding a polypeptide of 247 amino acid residues from 'Jinba' (Table S1). Compared to the protein sequence of CmWRKY8, only two amino acids were found to be different. Therefore, we named the gene CmWRKY8-1. CmWRKY8-1 contains 27 negatively charged residues and 33 positively charged residues. The molecular weight is 27,704.47 Da, and the theoretical pI is 8.80. The instability index is computed as 48.14, which classifies the protein as unstable. The grand average hydropathicity is −1.079, which makes it a hydrophilic protein. CmWRKY8-1 contains a WRKY domain of approximately 60 amino acids that contains a WRKYGQK motif and a zinc finger motif (C-X7-C-X23-H-X1-C), belonging to class III WRKY transcription factors ( Figure 1A). Phylogenetic tree analysis showed that CmWRKY8-1 had the highest sequence similarity to Artemisia annua AaWRKY8 ( Figure 1B).

Characteristics of CmWRKY8-1
To determine the subcellular localization of CmWRKY8-1, the pORE R4-GFP-CmWRKY8-1 vector or pORE R4-GFP empty vector was infiltrated into the epidermal cells

Characteristics of CmWRKY8-1
To determine the subcellular localization of CmWRKY8-1, the pORE R4-GFP-CmWRKY8-1 vector or pORE R4-GFP empty vector was infiltrated into the epidermal cells of N. benthamiana leaves. The GFP fluorescence signal of tobacco cells transformed with the pORE R4-GFP-CmWRKY8-1 vector appeared only in the nucleus and overlapped with the nuclear marker ( Figure 2A). In conclusion, CmWRKY8-1 was localized in the nucleus. In order to explore the expression pattern of CmWRKY8-1 in various tissues during the vegetative period of chrysanthemums, the terminal bud, stem, leaf and root of 'Jinba' were sampled and analyzed by qRT-PCR. The results showed that the relative expression level was highest in the roots, followed by the leaves ( Figure 2D). The pCL1 and pGBKT7 plasmids were used as positive and negative controls, respectively. (C) GLOMAX chemiluminescence determination of relative luciferase activity. The different letters mean significant differences according to Duncan's multiple range test at p < 0.05; the same scheme applies below. (D) The relative expression of CmWRK8-1 in bud, stem, leaf and root during the nutritional growth period. The 2 −ΔΔCt method was used to calculate relative transcript abundances. for RNA extraction, cDNA reverse transcription, and qRT-PCR data analysis. The results showed a significant decrease in CmWRKY8-1 expression in the experimental group compared to that in the control ( Figure 3A), suggesting that CmWRKY8-1 may play an important role in the response to F. oxysporum.

Expression Pattern of
Recent studies have found that when F. oxysporum infests chrysanthemums, the salicylic acid O-β-glucoside (SAG) content in chrysanthemums increases, triggering systemic defense [29]. To investigate whether SA affected the resistance of chrysanthemums to F. oxysporum, we sprayed exogenous SA. The experimental group was inoculated 24 h after spraying with 200 μm SA, and the control group was inoculated 24 h after spraying with sterile water. By observation, we found that chrysanthemums sprayed with SA showed The pCL1 and pGBKT7 plasmids were used as positive and negative controls, respectively. (C) GLO-MAX chemiluminescence determination of relative luciferase activity. The different letters mean significant differences according to Duncan's multiple range test at p < 0.05; the same scheme applies below. (D) The relative expression of CmWRK8-1 in bud, stem, leaf and root during the nutritional growth period. The 2 −∆∆Ct method was used to calculate relative transcript abundances.
The pDEST-GBKT7-CmWRKY8-1 plasmid was transformed into the yeast strain Y2H to verify its transcriptional activation activity. Yeast transformed with the positive control pCL1 plasmid grew normally, whereas yeast transformed with the negative control pDEST-GBKT7 plasmid and pDEST-GBKT7-CmWRKY8-1 plasmid did not grow normally on SD/-Ade/-His deficient medium ( Figure 2B). These results indicate that CmWRKY8-1 has no transcriptional activation activity. To further determine the transcriptional activity of CmWRKY8-1, we performed a transcriptional activity analysis of chrysanthemum protoplasts. The 35S::GAL4DB-CmWRKY8-1, 35S::GAL4DB-AtARF5, and 35S::GAL4DB plasmids were cotransformed with equal amounts of 5×GAL4-LUC in chrysanthemum protoplasts. CmWRKY8-1 luciferase activity values were significantly lower than those of the positive control, as measured using the GLOMAX chemiluminescence meter ( Figure 2C). Therefore, CmWRKY8-1 showed no transcriptional activation activity.
In order to explore the expression pattern of CmWRKY8-1 in various tissues during the vegetative period of chrysanthemums, the terminal bud, stem, leaf and root of 'Jinba' were sampled and analyzed by qRT-PCR. The results showed that the relative expression level was highest in the roots, followed by the leaves ( Figure 2D).

Expression
Pattern of CmWRKY8-1 after F. oxysporum Infection in C. moriflium 'Jinba' and Spraying Exogenous SA to Improve the Resistance of C. moriflium 'Jinba' to F. oxysporum To verify whether CmWRKY8-1 responds to F. oxysporum infection, we collected root samples at 0 h, 3 h, 6 h, 12 h, 24 h, 72 h, and 120 h after F. oxysporum inoculation in 'Jinba' for RNA extraction, cDNA reverse transcription, and qRT-PCR data analysis. The results showed a significant decrease in CmWRKY8-1 expression in the experimental group compared to that in the control ( Figure 3A), suggesting that CmWRKY8-1 may play an important role in the response to F. oxysporum. greater resistance ( Figure 3B). These results indicate that SA plays a positive defensive role in chrysanthemums in response to F. oxysporum infestation.

Overexpression of CmWRKY8-1-VP64 Fusion Protein Increases the Susceptibility of Chrysanthemums to F. oxysporum
To clarify the function of CmWRKY8-1, we obtained CmWRKY8-1 transgenic lines that overexpressed the CmWRKY8-1-VP64 fusion protein ( Figure 4A). We selected FVuv-1 and FVuv-2 as subjects for the CmWRKY8-1 transgenic lines. Through PCR amplification with the 35S forward primer and CmWRKY8-1 reverse primer, transgenic lines of CmWRKY8-1 were identified ( Figure 4B). The transgenic lines of CmWRKY8-1 were further analyzed by qRT-PCR ( Figure 4C). Recent studies have found that when F. oxysporum infests chrysanthemums, the salicylic acid O-β-glucoside (SAG) content in chrysanthemums increases, triggering systemic defense [29]. To investigate whether SA affected the resistance of chrysanthemums to F. oxysporum, we sprayed exogenous SA. The experimental group was inoculated 24 h after spraying with 200 µm SA, and the control group was inoculated 24 h after spraying with sterile water. By observation, we found that chrysanthemums sprayed with SA showed greater resistance ( Figure 3B). These results indicate that SA plays a positive defensive role in chrysanthemums in response to F. oxysporum infestation.
To clarify the function of CmWRKY8-1, we obtained CmWRKY8-1 transgenic lines that overexpressed the CmWRKY8-1-VP64 fusion protein ( Figure 4A). We selected FVuv-1 and FVuv-2 as subjects for the CmWRKY8-1 transgenic lines. Through PCR amplification with the 35S forward primer and CmWRKY8-1 reverse primer, transgenic lines of CmWRKY8-1 were identified ( Figure 4B). The transgenic lines of CmWRKY8-1 were further analyzed by qRT-PCR ( Figure 4C). After inoculation with F. oxysporum, we observed that FVuv-1 and FVuv-2 already showed a clear wilting state on the 8th day, whereas WT lines had only a slight wilting at After inoculation with F. oxysporum, we observed that FVuv-1 and FVuv-2 already showed a clear wilting state on the 8th day, whereas WT lines had only a slight wilting at the base ( Figure 5A). On the 12th day, the transgenic seedling disease of the FVuv-1 and FVuv-2 lines had spread from the middle to the top of the stem, and the entire plant had died ( Figure 5B). The WT lines still had only a few leaves withered at the base. We evaluated and graded chrysanthemum seedlings manually for the degree of root browning and the degree of yellowing and browning of leaves in diseased chrysanthemum seedlings. The DSI of WT lines was significantly lower than that of FVuv-1 and FVuv-2 ( Figure 5C). More importantly, we measured the activities of peroxidase (POD), catalase (CAT), phenylalaninammo-nialyase (PAL), and polyphenol oxidase (PPO) enzymes in plants after infection on the 8th day. POD, CAT, PAL, and PPO are critical enzymes in plant defense systems, and can be used as standards to measure plant resistance [30,31]. The results showed that enzyme contents in FVuv-1 and FVuv-2 were significantly lower than those in WT lines ( Figure 5D).
These results suggest that the overexpression of CmWRKY8-1-VP64 fusion protein can improve the susceptibility of chrysanthemums to F. oxysporum.

Changes in Genes Involved in the SA Signaling Pathway and Alterations in Endogenous SA
Plants can synthesize SA via ICS and PAL pathways. EDS5, PBS3, EDS1, and EPS1 are involved in the SA pathway [32]. To investigate whether CmWRKY8-1 is involved in regulating the response to F. oxysporum through SA, we determined the genes related to the SA pathway using qRT-PCR after inoculation with F. oxysporum. The results showed that the relative expression levels of ICS1, PAL, EDS5, PBS3, EPS1, and EDS1 in CmWRKY8-1 transgenic lines were significantly lower than those in WT lines ( Figure 6A). We also determined the endogenous SA content of plants after inoculation with F. oxysporum. We observed that the content of endogenous SA in CmWRKY8-1 transgenic lines was significantly lower than that in the WT lines ( Figure 6B). Therefore, CmWRKY8-1 may participate in the SA pathway in response to F. osporium infection. To determine whether CmWRKY8-1 could directly respond to F. oxysporum infection through the SA pathway, we measured the expression of disease-resistant defense genes downstream of SA. Here, we measured the relative expression of PR1, PR2, and PR5 after inoculation with F. oxysporum by qRT-PCR. The results showed that the transcription levels of PR1, PR2, and PR5 in CmWRKY8-1 transgenic lines were generally lower than those in the WT lines ( Figure 6C).
More importantly, we measured the activities of peroxidase (POD), catalase (CAT), phenylalaninammo-nialyase (PAL), and polyphenol oxidase (PPO) enzymes in plants after infection on the 8th day. POD, CAT, PAL, and PPO are critical enzymes in plant defense systems, and can be used as standards to measure plant resistance [30,31]. The results showed that enzyme contents in FVuv-1 and FVuv-2 were significantly lower than those in WT lines ( Figure 5D). These results suggest that the overexpression of CmWRKY8-1-VP64 fusion protein can improve the susceptibility of chrysanthemums to F. oxysporum.

Changes in Genes Involved in the SA Signaling Pathway and Alterations in Endogenous SA
Plants can synthesize SA via ICS and PAL pathways. EDS5, PBS3, EDS1, and EPS1 are involved in the SA pathway [32]. To investigate whether CmWRKY8-1 is involved in regulating the response to F. oxysporum through SA, we determined the genes related to the SA pathway using qRT-PCR after inoculation with F. oxysporum. The results showed that the relative expression levels of ICS1, PAL, EDS5, PBS3, EPS1, and EDS1 in CmWRKY8-1 transgenic lines were significantly lower than those in WT lines ( Figure 6A). We also determined the endogenous SA content of plants after inoculation with F. oxysporum. We observed that the content of endogenous SA in CmWRKY8-1 transgenic lines was significantly lower than that in the WT lines ( Figure 6B). Therefore, CmWRKY8-1 may participate in the SA pathway in response to F. osporium infection. To determine whether CmWRKY8-1 could directly respond to F. oxysporum infection through the SA pathway, we measured the expression of disease-resistant defense genes downstream of SA. Here, we measured the relative expression of PR1, PR2, and PR5 after inoculation with F. oxysporum by qRT-PCR. The results showed that the transcription levels of PR1, PR2, and PR5 in CmWRKY8-1 transgenic lines were generally lower than those in the WT lines ( Figure 6C). Therefore, we concluded that CmWRKY8-1 transgenic lines could potentially reduce endogenous SA by inhibiting the expression of SA-related genes, which decreases the expression of downstream disease-related proteins, leading to susceptibility to disease.

Transcriptome Sequencing Analysis and Functional Enrichment of DEGs in CmWRKY8-1 Transgenic Lines
We performed transcriptome sequencing to better elucidate the expression profiles of WT lines and CmWRKY8-1 transgenic chrysanthemum lines. The roots of WT lines and CmWRKY8-1 transgenic lines at 0 h, 3 h, and 72 h after infection were sampled with three Therefore, we concluded that CmWRKY8-1 transgenic lines could potentially reduce endogenous SA by inhibiting the expression of SA-related genes, which decreases the expression of downstream disease-related proteins, leading to susceptibility to disease.

Transcriptome Sequencing Analysis and Functional Enrichment of DEGs in CmWRKY8-1 Transgenic Lines
We performed transcriptome sequencing to better elucidate the expression profiles of WT lines and CmWRKY8-1 transgenic chrysanthemum lines. The roots of WT lines and CmWRKY8-1 transgenic lines at 0 h, 3 h, and 72 h after infection were sampled with three biological replicates, and total RNA was extracted for sequencing. In total, 861,033,134 reads were generated, resulting in 856,387,816 clean reads after filtration.
Biological replicates in different groups showed a high degree of consistency, while both WT lines and CmWRKY8-1 transgenic lines in 72 h showed different expression profiles compared to the other groups ( Figure S1). Differentially expressed genes (DEGs) were identified using the following threshold criteria: FDR ≤ 0.05 and |log2| ≥ 1. There were 5745, 5840, and 1527 DEGs between the WT and CmWRKY8-1 transgenic lines at 0 h, 3 h, and 72 h, respectively ( Figure 7A). The Venn diagram showed that there were 191 DEGs in the three comparisons ( Figure S2). To investigate the function of differential genes at each time point, the DEGs of three time-point comparisons between WT lines and transgenic lines were annotated using GO analysis ( Figure 7B). In the GO enrichment analysis, SA-associated pathways were enriched, including "Response to salicylic acid," "Salicylic acid mediated signaling pathway," "Salicylic acid metabolic process," "Cellular response to salicylic acid stimulus," "Systemic acquired resistance", and "salicylic acid mediated signaling pathway". To further verify the function of SA, we quantified the expression of genes in the SA pathway ( Figure S3). Compared with WT plants, the expression levels of two transcripts of PAL (evm.TU.scaffold_665.57, evm.TU.scaffold_1462.189), four transcripts of EDS1 (evm.TU.scaffold_6081.14, evm.TU.scaffold_248.42, evm.TU.scaffold_248.45, evm.TU.scaffold_248.46), one transcript of NPR1 (evm.TU.scaffold_11524.9), and two transcripts of AIM1 (MSTRG.77477, MSTRG.142194) were lower in CmWRKY8-1 transgenic lines. To verify the authenticity of the transcriptome data, we validated the DEGs using qRT-PCR against the SA pathway described above. The results showed that the trends of transcript expression changes obtained by qRT-PCR were consistent with the transcriptomic data ( Figure 8). This indicated that CmWRKY8-1 could influence F. oxysporum infection by mediating the SA pathway genes. transgenic lines were annotated using GO analysis ( Figure 7B). In the GO enrichment analysis, SA-associated pathways were enriched, including "Response to salicylic acid," "Salicylic acid mediated signaling pathway," "Salicylic acid metabolic process," "Cellular response to salicylic acid stimulus," "Systemic acquired resistance", and "salicylic acid mediated signaling pathway". To further verify the function of SA, we quantified the expression of genes in the SA pathway ( Figure S3). To verify the authenticity of the transcriptome data, we validated the DEGs using qRT-PCR against the SA pathway described above. The results showed that the trends of transcript expression changes obtained by qRT-PCR were consistent with the transcriptomic data ( Figure 8). This indicated that CmWRKY8-1 could influence F. oxysporum infection by mediating the SA pathway genes.

Discussion
In this experiment, we demonstrated that the CmWRKY8-1 transgenic strain increased susceptibility to F. oxysproum by overexpressing CmWRKY8-1-VP64. Furthermore, we found that CmWRKY8-1 could influence the resistance of 'Jinba' to F. oxysproum by regulating the SA pathway through transcriptomic data analysis. WRKY transcription factors play a complex regulatory role in plant defence signaling pathways. Currently, an increasing number of studies on WRKY responses to F. oxysporum infestation are also being conducted. LrWRKY3 heterologously specializes in tobacco, which exhibits enhanced resistance to F. oxysporum infestation, and transient expression of the LrWRKY3 RNAi vector in Lilium regale scales increased the susceptibility to F. oxysporum [33]. In cotton, group IIc WRKY TFs enhance resistance to F. oxysporum by promoting flavonoid synthesis via the WRKY-MAPK pathway [34]. In chickpeas, CaWRKY70 inhibits multiple immune signaling pathways, including CaMPK9-CaWRKY40 signaling, thereby negatively regulating resistance to F. oxysporum [17,35]. However, there are few reports on WRKY-responsive F. oxysporum in chrysanthemums. To investigate the function of WRKY transcription factors in response to F. oxysporum infestation, we selected and cloned CmWRKY8-1 based on previous results [28]. CmWRKY8-1 is a transcription factor that is localized in the nucleus of cells (Figure 2A,B). Based on quantitative tissue analysis, CmWRKY8-1 was found to be highly expressed in roots ( Figure 2D). We speculate that the high expression of CmWRKY8-1 in the roots may facilitate its positive response to F. oxysporum infestation.
In the study of transcription factors, the fusion of transcription factors with VP64 proteins can have a more significant regulatory effect if VP64 can convert transcriptional repressors into transcriptional activators [36]. We took advantage of the characteristics of

Discussion
In this experiment, we demonstrated that the CmWRKY8-1 transgenic strain increased susceptibility to F. oxysproum by overexpressing CmWRKY8-1-VP64. Furthermore, we found that CmWRKY8-1 could influence the resistance of 'Jinba' to F. oxysproum by regulating the SA pathway through transcriptomic data analysis. WRKY transcription factors play a complex regulatory role in plant defence signaling pathways. Currently, an increasing number of studies on WRKY responses to F. oxysporum infestation are also being conducted. LrWRKY3 heterologously specializes in tobacco, which exhibits enhanced resistance to F. oxysporum infestation, and transient expression of the LrWRKY3 RNAi vector in Lilium regale scales increased the susceptibility to F. oxysporum [33]. In cotton, group IIc WRKY TFs enhance resistance to F. oxysporum by promoting flavonoid synthesis via the WRKY-MAPK pathway [34]. In chickpeas, CaWRKY70 inhibits multiple immune signaling pathways, including CaMPK9-CaWRKY40 signaling, thereby negatively regulating resistance to F. oxysporum [17,35]. However, there are few reports on WRKY-responsive F. oxysporum in chrysanthemums. To investigate the function of WRKY transcription factors in response to F. oxysporum infestation, we selected and cloned CmWRKY8-1 based on previous results [28]. CmWRKY8-1 is a transcription factor that is localized in the nucleus of cells (Figure 2A,B). Based on quantitative tissue analysis, CmWRKY8-1 was found to be highly expressed in roots ( Figure 2D). We speculate that the high expression of CmWRKY8-1 in the roots may facilitate its positive response to F. oxysporum infestation.
In the study of transcription factors, the fusion of transcription factors with VP64 proteins can have a more significant regulatory effect if VP64 can convert transcriptional repressors into transcriptional activators [36]. We took advantage of the characteristics of VP64 and fused CmWKRY8-1, which has no transcriptional activation activity, to VP64 to make it a transcriptional activator. Transgenic chrysanthemum lines overexpressing CmWKRY8-1-VP64 were generated using transgenic techniques. CmWKRY8-1 transgenic lines exhibited reduced resistance to F. oxysporum ( Figure 5).
Plants have intricate defense systems and plant hormones play a vital role. SA plays an important role in plant immune processes by inducing systemic acquired resistance (SAR) in plants [37,38]. The synthesis of SA in plants is divided into PAL and ICS pathways, and it has been shown that the ICS pathway produces most of the SA produced by disease resistance induction [39,40]. SA biosynthesis and signal transduction are enhanced in plants upon pathogen infestation, and SA induces the expression of disease resistancerelated genes, thereby improving plant resistance to disease. When F. oxysporum infested chrysanthemums, the SAG content in chrysanthemums increased [29]. Our experiments also demonstrated that chrysanthemum resistance to F. oxysporum was enhanced when exogenous SA was sprayed ( Figure 3B). In soybean, GmWRKY31 mediates resistance to Peronospora manshurica by regulating the expression of the GmSAGT1 gene and participating in the SA pathway [41]. In a study on cotton, the GhMKK4-GhMPK20-GhWRKY40 cascade was found to reduce resistance to F. oxysporum by blocking the SA-mediated defense pathway [42]. Silencing GhWRKY70 in cotton can reduce the sensitivity of cotton to Verticillium dahliae by downregulating SA gene expression [43]. In woody plants, the SA-mediated defense genes PR1, PR2, and PAD4 were significantly upregulated in poplar after the overexpression of PtrWRKY73 in vivo [44]. It has also been shown that CmWRKY15-1 can enhance resistance to Puccinia horiana Henn through the SA pathway by interacting with NPR1 [45]. In the present study, the key genes of the SA synthesis pathway were detected by qRT-PCR. We observed that the relative expression of all these genes was lower in CmWRKY8-1 transgenic lines than in WT lines after inoculation ( Figure 6A). In the transgenic CmWRKY8-1 lines, the endogenous SA content was lower than that of the WT lines ( Figure 6B). Therefore, we tentatively hypothesized that CmWRKY8-1 might be involved in the SA signaling pathway and may affect the synthesis and degradation of SA. We also focused on PR genes downstream of SA. The relative expression levels of PR1, PR2, and PR5 were lower in the CmWRKY8-1 transgenic lines than in the WT lines ( Figure 6C). RNA-seq analysis revealed that the relative expression of DEGs associated with the SA pathway was lower in CmWRKY8-1 transgenic lines (Figure 8). A GO analysis revealed that GOs were associated with the SA signaling pathway ( Figure 7B). In summary, we suggest that CmWRKY8-1 affects resistance to F. oxysporum by regulating the SA pathway. However, the mechanism of how CmWRKY8-1 directly or indirectly regulates the SA pathway needs to be further explored.
Under the conditions of high temperature, high humidity and continuous cultivation, large amounts of F. oxysporum will accumulate in the soil, resulting in the frequent occurrence of chrysanthemum Fusarium wilt. We expect to apply transgenic chrysanthemum to field cultivation, screen the resistant germplasm of chrysanthemum by transgenic means, and cultivate resistant and high-quality varieties. We also hope to link the environment, molecule and F. oxysporum to study the disease resistance mechanism of chrysanthemum, improving the quality of chrysanthemum from many angles.

Plant Materials and Growth Conditions
The chrysanthemum variety 'Jinba' used in this experiment was provided by the Chrysanthemum Germplasm Resource Preserving Center, Nanjing Agricultural University (Nanjing, China). Chrysanthemum cuttings were planted in a 1:2 (v/v) mixture of soil and vermiculite. Chrysanthemums were grown in a greenhouse with a photoperiod of 16 h/8 h (light/dark), a temperature of 25 • C, and a humidity of 70%.

Isolation, Identification, Culture and Inoculation of Pathogenic Fungus
The diseased parts of plants were observed in the field, and the diseased plants were collected. The fungus was isolated by the tissue block method and inoculated on a PDA plate after purification. The colony morphology was observed, and the genomic DNA of fungus was extracted and sequenced. Finally, the fungal species was observed by re-inoculation and phenotyping.
F. oxysporum used in the experiment was isolated from the root and soil of chrysanthemum variety 'Jinba' at the Chrysanthemum Germplasm Resource Preserving Center (Nanjing, China) and stored at −80 • C. Before inoculation, the preserved fungal cakes were inoculated on PDA plates and cultured at 28 • C for six days. Then, 10 fungal cakes with a diameter of 0.7 cm were cut and inoculated into 500 mL PDB medium at 28 • C and 180 rpm, and the culture was shaken for 5 days. Spore concentration was determined using a blood counting plate. For inoculation, chrysanthemum roots were cut lightly and soaked in a spore suspension at a concentration of 10 7 CFU ml −1 for 30 min. When each plant finally colonized, the seedlings were inoculated with 1 × 10 7 spores per gram of substrate. The chrysanthemums were cultured in an environment of 16 h/8 h (light/dark), 28 • C and 80% humidity in a photo-culture chamber.

Isolation and Sequence Analysis of CmWRKY8-1
Total RNA was extracted from snap-frozen roots of the cultivar 'Jinba' using an RNA extraction kit (HuaYueyang), and reverse transcription was performed using a rapid reverse transcription kit (TaKaRa) to obtain cDNA as a template for gene cloning. Primers (Table S2) were designed for PCR amplification according to the gene sequences logged in the NCBI database (KC615362). PCR products were purified, constructed into the pMD19-T (TaKaRa) vector, and sequenced. The CmWRKY8-1 homologue peptide sequence was retrieved from TAIR (https://www.arabidopsis.org/, accessed on 10 October 2022) and the NCBI web site (https://www.ncbi.nlm.nih.gov, accessed on 10 October 2022). CmWRKY8-1 was used to perform multiple sequence alignments of homologous sequences using DNAMAN6, and MEGA X was used to construct phylogenetic trees using the adjacency method. The online tool ExPASy (http://expasy.org, accessed on 10 October 2022) was used to predict the physicochemical properties of CmWRKY8-1.

Transactivation Assays of CmWRKY8-1
A yeast assay was performed to examine the transcriptional activity of CmWRKY8-1 [47]. The pGBKT7-CmWRKY8-1 vector was constructed. pCL1 (positive control), pG-BKT7 (negative control), and pGBKT7-CmWRKY8-1 were transformed into the yeast cells Y2HGold. pCL1 was coated in SD/Leu-media, and the pGBKT7 and pGBKT7-CmWRKY8-1 vectors were coated in SD/Trp-media and incubated for 3-5 days at 28 • C in an inverted position. Finally, colonies were picked and spotted according to different concentration gradients plates onto SD/Ade-His-media with and without X-α-gal, and incubated overnight at 28 • C in an inverted position. Colonies were observed and photographed.

Quantitative Real-Time PCR (qRT-PCR)
Quantitative primers (Table S2) were designed using Primer5, and EF1α (Table S2) was used as the internal reference gene [49]. qRT-PCR was performed using the SYBR Green PCR Master Mix (TaKaRa). The 96-well plate was placed in a Mastercycler ep realplex 2S (Eppendorf, Hamburg, Germany) qRT-PCR instrument, and the fluorescence acquisition channel and fluorescence reading step were set. PCR reactions were performed according to the following reaction conditions: 95 • C for 2 min, 95 • C 15 s, 55 • C 15 s, 72 • C 20 s, 40 cycles. Finally, the dissolution curve program was developed. Calculations were performed using the 2 −∆∆Ct method [50].

Analysis of CmWRKY8-1 Expression under Stress Treatments
40 days old seedlings were used to determine the expression patterns of CmWRKY8-1 under different stress treatments. For the inoculation treatment, the experimental group was treated with F. oxysporum, whereas the control group was treated with sterile water. Root samples were collected at 0 h, 3 h, 6 h, 12 h, 24 h, 72 h, and 120 h after treatment with F. oxysporum. For the SA treatment, the experimental group was sprayed with 200 µm SA, while the control group was sprayed with sterile water.

Analysis of F. oxysporum Resistance in CmWRKY8-1 Transgenic Chrysanthemums
Twenty seedlings of 40 days old wild-type and transgenic chrysanthemums were used for inoculation treatment. The disease severity index (DSI) was calculated on days 8 and 12 after inoculation, and the phenotypes were photographed and recorded. The DSI was calculated according to the following formula: DSI = ∑ (number of diseased plants per level × number of corresponding levels) ×100 / (total number of plants surveyed × highest disease level) (Table S3).

Enzyme Activity Assay and Endogenous SA Determination
Root samples were taken from transgenic lines and WT lines on the 8th day after inoculation. The enzyme activity in chrysanthemums was measured using POD, CAT, PAL, and PPO kits (Comin, Suzhou, China). In the meantime, the endogenous SA content was determined by a Plant Hormone Salicylic Acid ELISA Kit (Lengton Bioscience Co., Shanghai, China).

RNA-Seq Analysis
Seedlings of the WT lines and CmWKRY8-1 transgenic lines were planted in a greenhouse at Nanjing Agricultural University for 40 days. After inoculation at 0 h, 3 h, and 72 h, roots of WT lines and CmWKRY8-1 transgenic lines were sampled with three biological replicates. Total RNA was extracted using a TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA) and assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Transcriptome sequencing was performed using an Illumina Novaseq6000 made by the Gene Denovo Biotechnology Co. (Guangzhou, China). The reads were filtered and aligned using Bowtie2 (version 2.2.8) [52]. Paired-end clean reads were mapped to the reference genome using the HISAT2 software version 2.4 [53]. Differentially expressed genes (DEGs) between two different groups were identified using DESeq2 software [54]. The transcripts were annotated using the BLASTX algorithm [55] and Gene Ontology database [56]. The quantification of transcripts was transformed into FPKM values with RSEM [57].

Data Analysis
Statistical analyses were performed using SPSS version 25.0. All data were analyzed using analysis of variance (ANOVA) and t-tests to determine significant differences.

Conclusions
In conclusion, we identified a WRKY gene in chrysanthemums, CmWRKY8-1, and found that it was responsive to F. oxysporum infestation. By overexpressing CmWRKY8-1-VP64, we found that the resistance of chrysanthemums to F. oxysporum was reduced. Moreover, the expression of genes related to the SA signaling pathway and endogenous SA content was decreased in CmWRKY8-1 transgenic lines. RNA-seq analysis showed that the expression of DEGs of the SA signaling pathway was decreased in CmWRKY8-1 transgenic lines. Therefore, we demonstrated that CmWRKY8-1 responds to F. oxysporum infection by regulating the SA signaling pathway.