DgMYB2 improves cold resistance in chrysanthemum by directly targeting DgGPX1

Abstract MYB transcription factors play important roles in plant responses to cold stress, but the associated underlying mechanisms remain unclear. In this study, a cold-induced MYB transcription factor, DgMYB2, was isolated from chrysanthemum (Chrysanthemum morifolium Ramat). DgMYB2 was localized to the nucleus and exhibited transactivational activity. Overexpression of DgMYB2 improved cold tolerance in chrysanthemum, while cold tolerance in the antisense suppression lines decreased compared to that of the wild type. Additionally, electrophoretic mobility shift assays, chromatin immunoprecipitation, luciferase complementary imaging analysis, and dual-luciferase reporter gene detection experiments confirmed that DgMYB2 directly targets DgGPX1 and increases the activity of glutathione peroxidase to reduce the accumulation of reactive oxygen species, thereby improving cold resistance in chrysanthemum.


Introduction
Cold stress is a common abiotic stress factor that negatively affects plant growth and is one of the leading factors restricting the distribution of plants. Cold stress can be classified as either chilling (0-15 • C) stress or freezing (< 0 • C) stress [1]. Plants have evolved defense mechanisms to tolerate cold stress. When plants are subjected to cold stress, transcription factors are activated through a series of signal transduction pathways, ultimately improving plant cold tolerance. Activated transcription factors specifically bind to their corresponding cis-acting elements and activate the expression of a series of downstream cold resistance-related genes, which then trigger a series of physiological and biochemical reactions, thereby improving plant cold resistance. Members of transcription factor families have been identified in plants, including members of NAC, WRKY, and AP2/ER families, which are some of the largest transcription factor families in plants. Specifically, members of the MYB transcription factor family are largely involved in plant cell morphology, plant physiological metabolism, and stress response processes [2][3][4][5].
Members of the MYB transcription factor family are named after their N-terminal conserved DNA-binding region, the MYB domain, which is composed of 1-4 incomplete repetitive (R) sequences consisting of 50-52 amino groups in a series [6]. Based on the different number of R repetitions in a sequence, MYB transcription factors can be divided into four subfamilies [7].
Members of the R1-MYB subfamily contain only one R domain and are referred to as R1-MYB/MYB-related transcription factors. The R2R3-MYB subfamily has the most members among the MYB transcription factor family, the members of which contain two R structures at the N-terminus. R1R2R3-MYB transcription factors contain three R domains and compose a relatively small MYB transcription factor subfamily, and research on the functions of R1R2R3-MYB transcription factors is scarce. The fourth subfamily, the 4R-MYB subfamily, has the fewest members, and their functions remain poorly understood [8,9].
Of these subfamilies, the most widely studied is the R2R3-type MYB subfamily [10]. Research on the regulatory mechanisms of R2R3-MYBs in the plant cold tolerance process has mainly been conducted in apple and pear. Overexpression of the cold-inducible genes MdMYB88/124 significantly increased cold tolerance in apple, while a reduction in the transcript levels of these two genes weakened cold tolerance [11]. The expression of MdMYB23 was upregulated under cold stress, and its overexpression enhanced cold tolerance in transgenic apple calli and Arabidopsis plants [7]. Moreover, the overexpression of MdMYB308L significantly increased cold tolerance in apple plants [12]. In pear (Pyrus betulifolia Bunge), cold stress induced the transcription of PbrMYB5, and its overexpression improved the cold resistance of transgenic tobacco [13]. However, the mechanisms of R1-MYBs in response to cold stress in plants remain unclear, necessitating further research.
Reactive oxygen species (ROS) are the products of plant metabolism and have extremely high activity and toxicity [14]. Under normal growth conditions, ROS can be eliminated by the antioxidant defense mechanisms of plants. Under environmental stress, the balance between ROS production and elimination in plants is disrupted and leads to the rapid accumulation of ROS and oxidative damage [15]. Glutathione peroxidase (GPX) is an important enzyme in the ROS-scavenging system and plays a role in reducing the accumulation of ROS caused by abiotic stress and maintaining the integrity of the plasma membrane [16]. However, the regulatory mechanisms underlying the effects of MYB transcription factors on GPX remain unclear.
Chrysanthemums are often used as fresh cut f lowers and have important ornamental and economic value [17]. As a common stress factor, cold stress hinders the productivity and economic value of chrysanthemum [18]. Studying the responses of chrysanthemum to cold stress will improve chrysanthemum production. In this study, a cold-induced MYB transcription factor, DgMYB2, was isolated from chrysanthemum. We found that DgMYB2 was localized in the nucleus and exhibited transcriptional activity. Moreover, overexpression of DgMYB2 enhanced cold tolerance in chrysanthemum, while reducing its expression reduced the survival rate of transgenic plants under cold stress. Based on the collective findings, we provide evidence that DgMYB2 plays a positive role in plants under cold stress by regulating the expression of DgGPX1 directly to increase GPX activity.

Isolation and expression patterns of DgMYB2
Cold-treated (4 • C for 24 h followed by −4 • C for 4 h) and untreated (25 • C) chrysanthemum seedlings were sampled for transcriptome analyses (NCBI accession No. GSE117262). The results revealed that 43 MYB genes were significantly (P < 0.05) induced at the transcriptional level in response to cold treatment (Table S1). Among these genes, DgMYB2 (GenBank accession No. MW853963) was notably induced (log 2 [fold change] = 2.7918) and was thus selected for further investigation. The full-length cDNA of DgMYB2 contained an open reading frame (ORF) of 885 bp and encoded a predicted protein of 294 amino acids with a molecular weight of 32.23 kDa. The amino acid sequence of DgMYB2 contained a SANT/MYB conserved domain, indicating that DgMYB2 is a typical R1-MYB transcription factor (Fig. S1a). Phylogenetic analysis revealed that the DgMYB2 sequence is highly homologous with other R1-MYB transcription factor sequences of various plant species and is closely related to AaMYB-like from Artemisia annua (Fig. S1b).
The expression patterns of DgMYB2 in different tissues showed that DgMYB2 was most highly expressed in the leaves, followed by the stems, and it was expressed the least in the roots (Fig. 1a). We also found that the expres-sion of DgMYB2 increased under cold treatment, reaching the highest value after 12 h and then decreasing slightly ( Fig. 1 b).
A dual-LUC reporter gene assay (DLA) and luciferase complementation imaging (LCI) were performed to verify the function of the natural promoter of DgMYB2 (1500 bp). The results demonstrated that DgMYB2 responded to cold stress (Fig. 1c, d).

DgMYB2 was localized to the nucleus and had transactivational activity
To determine the subcellular localization of DgMYB2, pSuper1300-DgMYB2-GFP and the nuclear marker pSuper::NF-YA4-mCherry were transiently expressed in tobacco leaves, while pSuper1300-GFP and the nuclear marker pSuper::NF-YA4-mCherry were coexpressed in the epidermal cells of tobacco leaves as controls. The results indicated that DgMYB2 was localized within the nucleus (Fig. 2a).
A dual-luciferase reporter gene analysis system was used to evaluate the transcriptional activity of DgMYB2. The results revealed that the LUC/REN value of the experimental group whose members were coinjected with the reporter gene and pBD-DgMYB2 was significantly (P < 0.05) higher than that of the negative control (pBD-Empty) (Fig. 2b), indicating that DgMYB2 had transcriptional activation activity. To further analyze the location of the transcriptional activation region of DgMYB2, a recombinant plasmid containing the corresponding sequence fragment was transferred into yeast onehybrid (Y1H) competent cells and then cultured on corresponding media. Yeasts transformed with the complete ORF sequence (1-294 aa) of DgMYB2 and yeasts transformed with the C-terminal sequence (145-294 aa) of DgMYB2 grew normally and appeared blue on SD/-H/−A plates including X-α-gal (Fig. 2c). Taken together, these results indicated that DgMYB2 had transcriptional activation activity and that the responsible region which was located within the C-terminal region (145-294 aa) of DgMYB2.

DgMYB2 played a positive regulatory role in the chrysanthemum response to cold stress
Transgenic chrysanthemum lines were obtained by Agrobacterium-mediated transformation of wild-type (WT) chrysanthemum. Two overexpression lines (DgMYB2 OE-6 and DgMYB2 OE-31) and two RNA interference (RNAi) lines (DgMYB2 Ri-11 and DgMYB2 Ri-45) were selected for further experimentation. Quantitative real-time PCR (qRT-PCR) was used to determine the expression level of DgMYB2 in WT and transgenic lines (Fig. 3a). Following cold treatment, the RNAi lines were the most severely withered, followed by the WT, while the overexpression lines maintained relatively good growth (Fig. 3b). Following cold treatment, plants were subjected to a 2-week recovery period at room temperature (25 • C) to analyze the survival rate of the WT, overexpression, and RNAi lines. The survival rate of the overexpression  lines was the highest, followed by the WT lines, while the RNAi lines had the lowest survival rate (Fig. 3c).
The histochemical staining results showed that, compared with those of the WT, the spots on the leaves of the overexpression lines were fewer in number and lighter in color. In contrast, compared with those of the WT, the spots on the leaves of the RNAi lines were more numerous and darker (Fig. 3d,   Under cold stress, the malondialdehyde (MDA) content and the relative electrolyte leakage in each line gradually tended to increase. Moreover, the MDA content and relative electrolyte leakage rate in the overexpression lines were the lowest, followed by the WT and RNAi lines (Fig. 3h, i). We measured the activities of antioxidant enzymes, including superoxide dismutase (SOD), peroxidase (POD), and glutathione peroxidase (GPX), which are components of the ROS-scavenging system. The results revealed that the activity of these enzymes increased to a certain degree when compared to values prior to cold treatment. No statistically significant differences in the SOD or POD activity between the transgenic, WT, and RNAi lines were detected (Fig. S2), but there were significant differences in GPX activity. The GPX activity of the overexpression lines was the highest, followed by that of the WT and RNAi lines (Fig. 3j). Furthermore, the qRT-PCR results were consistent with the GPX activity results. The overexpression lines had the highest expression level of DgGPX1, followed by the WT and RNAi lines (Fig. 3k).
Taken together, the results indicated that DgMYB2 positively affected the expression of DgGPX1 to regulate the activity of GPX, thereby reducing the accumulation of ROS to reduce plasma membrane damage, thus playing a positive regulatory role in the chrysanthemum response to cold stress.

DgMYB2 promoted DgGPX1 transcription by directly binding to the MYB-binding motif in the DgGPX1 promoter
The natural promoter sequence of DgGPX1 (1125 bp) contained a predicted MYB-binding site (CAACCA). The EMSA results confirmed that DgMYB2 was able to bind to the MYB-binding site (CAACCA) in the DgGPX1 promoter (Fig. 4a). The green fluorescent protein (GFP) expression levels in the transgenic lines were identified by western blots (Fig. S3). The transgenic lines were subjected to chromatin immunoprecipitation assay (ChIP)-qPCR experiments. ChIP-qPCR, LCI, and DLA experiments demonstrated that there was a protein-DNA interaction between DgMYB2 and the MYB-binding site (CAACCA) in the natural promoter sequence of DgGPX1 (Fig. 4b-d).

Discussion
To date, many MYB transcription factors have been isolated from different species [19][20][21]. The genes encoding these proteins are widely involved in the abiotic and biotic stress responses of plants. In this study, we isolated a cold-induced MYB transcription factor, DgMYB2, which was upregulated at the transcriptional level under cold stress (Fig. 1b). Phylogenetic analysis revealed that DgMYB2 was highly homologous with other plant R1-MYB proteins and had the closest relationship with that of A. annua (Fig. S1a, b). The expression levels of many novel MYB transcription factors have been shown to be induced by low temperature. In rice, the expression of OsMYB3R-2 is induced by low temperature [22]. In sheepgrass, LcMYB4 is upregulated in response to cold induction [23]. In tomato, the R2R3-MYB transcription factor LeAN2 is induced by low temperature and oxidative stress [24]. However, the functions and mechanisms of MYB transcription factors in response to cold stress in plants remain poorly understood.
Previous research on the function of MYB transcription factors at low temperatures revealed that these proteins affect the cold resistance responses of plants by regulating the ROS-scavenging system. LcMYB4 overexpression in sheepgrass resulted in higher POD activity, suggesting that the antioxidant ability of transgenic plants improved [23]. A MYB-related transcription factor, AQUILO, was isolated from grapevine (Vitis vinifera L.), and overexpression of AQUILO in Arabidopsis increased the POD and SOD activities and upregulated the expression levels of their encoding genes [25]. To verify the function of DgMYB2 in chrysanthemum under cold stress, we obtained overexpression lines (DgMYB2 OE-6 and DgMYB2 OE-31) and RNAi lines (DgMYB2 Ri-11 and DgMYB2 Ri-45). The expression level of DgMYB2 in the overexpression lines was significantly (P < 0.05) higher than that in the WT and RNAi lines, while the expression level of DgMYB2 in the WT lines was significantly (P < 0.05) higher than that in the RNAi lines (Fig. 3a). Compared with the WT lines, the overexpression lines experienced less damage from cold stress, and their survival rate was significantly (P < 0.05) higher, while the RNAi lines had the lowest survival (Fig. 3c). The accumulation of ROS in the overexpression lines was significantly (P < 0.05) lower than that in the WT lines, while the RNAi lines accumulated more ROS than the WT lines ( Fig. 3d-g). Under cold conditions, the relative electrolyte leakage and MDA content of the RNAi lines were the highest, followed by the WT and overexpression lines (Fig. 3h, i), indicating that, compared with the RNAi lines, the overexpression lines had a lower degree of membrane lipid peroxidation, experienced less peroxidation damage, and exhibited higher cold resistance. Additionally, the GPX activity of the overexpression lines was the highest, followed by the WT and RNAi lines (Fig. 3j). Moreover, the expression of DgGPX1 tended to be consistent with GPX activity (Fig. 3k). We found that DgMYB2 affected GPX activity by regulating the expression of DgGPX1, thereby inf luencing the ROS-scavenging process and degree of oxidative damage caused by ROS accumulation, thus positively regulating plant responses to cold stress.
MYB transcription factors bind to specific motifs in gene promoters to regulate the expression of cold stress- related genes to affect the cold stress resistance of plants [5]. In apple, the downstream target genes MdCCA1 and MdCSP3 of MdMYB88/124 were identified, and the overexpression of MdMYB88/124 contributed to increased cold tolerance [11]. MdMYB23 directly binds to the promoter regions of MdCBF1 and MdCBF2 and activates their expression to enhance plant cold resistance [7]. In pear, PbrMYB5 directly binds to the promoter region of PbrDHAR2 and activates its expression, which promotes the synthesis of ascorbic acid (AsA), in turn increasing plant cold resistance [13]. In this study, the EMSA, DLA, LCI, and ChIP-qPCR results demonstrated that DgMYB2 specifically bound to the MYB-binding site in the promoter region of DgGPX1. DgMYB2 directly targeted DgGPX1 to increase GPX activity, thus enhancing chrysanthemum cold resistance.
In conclusion, a cold-induced R1-MYB transcription factor, DgMYB2, was isolated from chrysanthemum, and its overexpression enhanced chrysanthemum cold resistance. DgMYB2 directly targeted DgGPX1 to increase GPX enzyme activity and thereby reduce ROS accumulation, thus improving chrysanthemum cold resistance (Fig. 5).

Plant materials and cold treatment
WT chrysanthemum was cultured on Murashige and Skoog (MS) media for 30 d under normal conditions (25 • C/22 • C day/night temperature conditions with a 16/8-h photoperiod). Afterward, these plant materials were used to obtain transgenic plants using Agrobacterium-mediated methods. After 30 d of cultivation, the chrysanthemum plantlets were transplanted into pots with a 1:1 peat:perlite mixture and cultivated under the aforementioned conditions for 7 d. Then, these plant materials were used for subsequent physiological experiments.
Cold-treated (4 • C for 24 h followed by −4 • C for 4 h) and untreated (25 • C) chrysanthemum seedlings were sampled for transcriptome analyses. The plant materials were subjected to 4 • C for 24 h and sampled after 0, 3, 6, 12, and 24 h of cold treatment for target gene expression level determination; relative electrolyte leakage assays; MDA, H 2 O 2 , and O 2 − content determination; and GPX activity assays. The samples were f lash-frozen in liquid nitrogen and used directly for subsequent experiments or otherwise stored at −80 • C. To verify the survival rates of transgenic and WT lines and record phenotypic changes, plant materials were exposed to chilling (4 • C for 24 h) and freezing (−6 • C for 6 h) treatments, followed by a recovery period of 2 weeks at room temperature (25 • C) to simulate the temperature conditions of the chrysanthemum production process.

RNA preparation and qRT-PCR
Chrysanthemum plants were sampled at the corresponding time points after cold stress treatment. An RNAprep Pure Plant Kit (Tiangen Biotech, Beijing, China) was used to obtain the total RNA from the chrysanthemum samples, and the obtained total RNA was reverse transcribed into cDNA using a One-Step gDNA Removal kit (TransGen Biotech, Beijing, China). The obtained cDNA was used as a template for qRT-PCR using a Two-Step qRT-PCR Kit (TransGen Biotech, Beijing, China); the elongation factor 1α gene (EF1α) was used as a stable reference gene. The results were detected using a CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The 2 -CT method was used for analysis, and the primers used for qRT-PCR are listed in Table S1.

Construction of expression vectors and generation of transgenic chrysanthemum
The National Center for Biotechnology Information (NCBI) website was used for the conserved domain analysis of DgMYB2. Homology analysis of DgMYB2 was performed by DNAMAN (Lynnon Biosoft, San Ramon, CA, USA), and phylogenetic analysis was conducted using MEGA (Molecular Evolutionary Genetics Analysis) v7. The ORF sequence of DgMYB2 was amplified using the primer pair M2-F and M2-R (Table S2). Then, the DgMYB2 ORF sequence was inserted into a modified pSuper1300 plasmid together with GFP to obtain a 35S::DgMYB2-GFP recombinant plasmid. The 35S::DgMYB2-GFP plasmid was used in subsequent subcellular localization experiments and generation of DgMYB2-overexpressing lines. To obtain RNAi lines, a 300 bp fragment of the ORF region containing no structural domain was selected and inserted into the pCAMBIA2301-GW-RNAi vector to construct pCAMBIA2301-DgMYB2-RNAi vectors. Construction of the recombinant plasmids was performed by Shanghai Sangon Biotech Co., Ltd. (Shanghai, China). To obtain transgenic chrysanthemum plants, the overexpression vector 35S::DgMYB2-GFP and RNAi vector pCAMBIA2301-DgMYB2-RNAi were transformed into WT chrysanthemum via Agrobacterium-mediated transformation. A DNAsecure Plant Kit (Tiangen, Beijing, China) was used to extract DNA from the obtained chrysanthemum plants to identify positive transgenic plants. Then, qRT-PCR-based detection of positive transgenic plants was performed. The operational steps were conducted following the manufacturer's instructions.

Dual-LUC reporter gene assays
To verify the function of the natural promoter of DgMYB2, the sequence of the DgMYB2 promoter (1500 bp) was inserted into a pSuper1300 vector together with the LUC reporter gene to form the recombinant plasmid ProDgMYB2-LUC. Renilla luciferase (REN) was used as the internal reference gene and inserted into a pSuper1300 vector to form the recombinant plasmid pSuper1300-REN. Construction of the recombinant plasmids was performed by Shanghai Sangon Biotech Co., Ltd. (Shanghai, China). The recombinant plasmids ProDgMYB2-LUC and pSuper1300-REN were then transferred into Agrobacterium, which were subsequently coinjected into tobacco for transient expression. Samples were collected from the leaves after transient expression, and the LUC activity was detected. A Dual Luciferase Reporter Gene Assay Kit (Beyotime Biotechnology, Shanghai, China) was used for analysis.
The natural DgGPX1 promoter sequence (1250 bp) was cloned and inserted into a pSuper1300 vector together with the LUC reporter gene to form the recombinant plasmid ProDgGPX1-LUC. To verify the direct regulatory activity of DgMYB2 on DgGPX1, ProDgGPX1-LUC and 35S::DgMYB2 were coinjected into tobacco leaves as a test group; tobacco leaves injected with ProDgGPX1-LUC were used as a control group. As an internal control, pSuper1300-REN was injected into tobacco leaves for transient expression. Samples were collected 48 h after transient expression, and a Dual Luciferase Reporter Gene Assay Kit (Beyotime Biotechnology, Shanghai, China) was used for analysis.

Transient expression assays in tobacco and luciferase complementation imaging
To detect the response of DgMYB2 to cold stress, ProDgMYB2-LUC was transformed into Agrobacterium GV3101, which was subsequently coinjected into tobacco leaves. After 36 h, D-luciferin potassium salt was injected, and fluorescence detection was performed using a plant luciferase in vivo imaging system.
To detect the regulatory activity of DgMYB2 on the DgGPX1 promoter, ProDgGPX1-LUC and 35S::DgMYB2-GFP were first transformed into Agrobacterium GV3101, which was then coinjected into tobacco leaves as a test group; tobacco leaves injected with the ProDgGPX1-LUC vector served as the control group. After 36 h, D-luciferin potassium salt was injected, and f luorescence detection was performed using a plant luciferase in vivo imaging system.

Transient expression assays and subcellular localization of DgMYB2
Bacterial liquid containing 35S::DgMYB2-GFP and bacterial liquid containing the nuclear marker pSuper::NF-YA4-mCherry were mixed together at a 1:1 ratio, and then the combined bacterial liquid was injected into tobacco leaves. The injected leaves were used as the experimental group. Concurrently, tobacco leaves coinjected with a pSuper1300-GFP empty vector and the nuclear marker pSuper::NF-YA4-mCherry were used as the control group. The method used for transient expression in tobacco leaves was described in a previous study [26]. After the injected tobacco plants were cultured for 36 h (23 • C ± 2 • C, 75% relative humidity), their leaves were cut at the injection site in the form of a square fragment (2 mm × 2 mm). A laser scanning confocal microscope (LSCM; Nikon Corporation, Tokyo, Japan) was used to observe the experimental results.

Transactivation analysis and Y1H assays
A pGreenII0800-LUC vector was used to insert 5 GAL4 sequences and 1 TATA sequence in front of LUC to construct a reporter gene vector. The ORF region of DgMYB2 was inserted into a modified pSuper1300 vector containing the GAL4 BD to construct the recombinant plasmid pBD-DgMYB2. VP16 was inserted into a modified pSuper1300 vector containing GAL4 BD as a recombinant plasmid pBD-VP16, which was used as a positive control. The modified pSuper1300 empty vector pBD-Empty containing GAL4 BD was used as a negative control. The obtained recombinant plasmids were used to transform Agrobacterium GV3101 competent cells, and the transformed bacterial liquid containing the reporter gene and the transformed bacterial liquid containing pBD-DgMYB2, pBD-VP16, and pBD-Empty were mixed at a 1:1 ratio and injected into tobacco leaves. After culturing for 3 d, the LUC/REN value was determined using a dual-luciferase reporter gene detection kit (Beyotime Biotechnology, Shanghai, China) following the manufacturer's instructions.

Electrophoretic mobility shift assays
The full-length ORF region of DgMYB2 was inserted into the pET28a expression vector to form a pET28a-DgMYB2 fusion construct. The recombinant plasmid was transformed into DE3 Escherichia coli competent cells, and 0.5 mM isopropyl bD-1 thiogalactopyranoside (IPTG) was added to induce the protein expression. Protein was extracted from the bacterial solution using a Bacterial Protein Extraction Kit (Sangon Biotech, Shanghai, China) following the manufacturer's instructions. Afterward, a Ni-NTA Sefinose Resin Kit (Sangon Biotech, Shanghai, China) was used for protein purification.
To synthesize general probes, a 100 bp DNA fragment containing a MYB element and an A base at the 5 end of the DgGPX1 promoter region was selected, and a 6-FAM tag was added at the 5 end. The sequence of the competitive probe was the same as that of the general probe, but there was no 6-FAM tag at the 5 end. When synthesizing the mutation probes, all sequences of the MYB site (CAACCA) were mutated to A bases, and the remaining sequences were the same as those of the general probe. Additionally, a His protein served as a negative control. The DNA probes used in the EMSA were all synthesized by Shanghai Sangon Biotech Co., Ltd. (Shanghai, China). The probes and samples were loaded onto a 6% polyacrylamide gel for the binding reaction and electrophoresis, and TBE buffer was used as the electrophoresis solution. The results were recorded using a charge-coupled device (CCD) camera.

Chromatin immunoprecipitation assays
The leaves of DgMYB2-overexpressing chrysanthemum plants were sampled and subjected to 1% formaldehyde cross-linking. Immunoprecipitation was performed using the SimpleChIP chromatin experimental procedure (magnetic beads) (Cell Signaling Technology, Danvers, Massachusetts, USA) following the manufacturer's instructions. ChIP-qPCR analysis was performed after purifying the resulting chromatin preparation, and a fragment containing the MYB-binding site (CAACCA) in the DgGPX1 promoter region was considered the functional area. GFP tag monoclonal antibodies were produced by Proteintech Group, Inc. (Chicago, IL, USA). The primers used in the ChIP-qPCR analysis are listed in Table S1.
The process and operating steps were described in a previous publication [27].

Determination of the physiological data of transgenic chrysanthemum
A quantitative measurement kit (Keming Biological Co., Ltd., Suzhou, China) was used for the quantitative determination of the H 2 O 2 and O 2 − contents. A GSH-PX assay kit (Jiancheng Bioengineering Institute, Nanjing, China) was used to determine the GPX activity. MDA accumulation was determined using an MDA assay kit (Jiancheng Bioengineering Institute, Nanjing, China), and relative electrolyte leakage was measured following previously described methods [28].