Comparative transcriptome analysis provides molecular insights into heterosis of waterlogging tolerance in Chrysanthemum indicum

BACKGROUND

Heterosis breeding is one of the most important breeding methods for chrysanthemum. To date, the genetic mechanisms of heterosis for waterlogging tolerance in chrysanthemum are still unclear. This study aims to analyze the expression profiles and potential heterosis-related genes of two hybrid lines and their parents with extreme differences in waterlogging tolerance under control and waterlogging stress conditions by RNA-seq.

RESULTS

A population of 140 F1 progeny derived from Chrysanthemum indicum (Nanchang) (waterlogging-tolerant) and Chrysanthemum indicum (Nanjing) (waterlogging-sensitive) was used to characterize the extent of genetic variation in terms of seven waterlogging tolerance-related traits across two years. Lines 98 and 95, respectively displaying positive and negative overdominance heterosis for the waterlogging tolerance traits together with their parents under control and waterlogging stress conditions, were used for RNA-seq. In consequence, the maximal number of differentially expressed genes (DEGs) occurred in line 98. Gene ontology (GO) enrichment analysis revealed multiple stress-related biological processes for the common up-regulated genes. Line 98 had a significant increase in non-additive genes under waterlogging stress, with transgressive up-regulation and paternal-expression dominant patterns being the major gene expression profiles. Further, GO analysis identified 55 and 95 transgressive up-regulation genes that overlapped with the up-regulated genes shared by two parents in terms of responses to stress and stimulus, respectively. 6,640 genes in total displaying maternal-expression dominance patterns were observed in line 95. In addition, 16 key candidate genes, including SAP12, DOX1, and ERF017 which might be of significant importance for the formation of waterlogging tolerance heterosis in line 98, were highlighted.

CONCLUSION

The current study provides a comprehensive overview of the root transcriptomes among F1 hybrids and their parents under waterlogging stress. These findings lay the foundation for further studies on molecular mechanisms underlying chrysanthemum heterosis on waterlogging tolerance.

CmDOF18 positively regulates salinity tolerance in Chrysanthemum morifolium by activating the oxidoreductase system

BACKGROUND

Chrysanthemum, one of the four major cut flowers all over the world, is very sensitive to salinity during cultivation. DNA binding with one finger (DOF) transcription factors play important roles in biological processes in plants. The response mechanism of CmDOF18 from chrysanthemum to salt stress remains unclear.

RESULTS

In this study, CmDOF18 was cloned from Chrysanthemum morifolium, and its expression was induced by salinity stress. The gene encodes a 291-amino acid protein with a typical DOF domain. CmDOF18 was localized to the nucleus in onion epidermal cells and showed transcriptional activation in yeast. CmDOF18 transgenic plants were generated to identify the role of this gene in resistance to salinity treatment. Chrysanthemum plants overexpressing CmDOF18 were more resistant to salinity stress than wild-type plants. Under salinity stress, the malondialdehyde content and leaf electrolyte conductivity in CmDOF18-overexpressing transgenic plants were lower than those in wild-type plants, while the proline content, chlorophyll content, superoxide dismutase activity and peroxidase activity were higher than those in wild-type plants. The opposite findings were observed in gene-silenced plants compared with wild-type plants. The gene expression levels of oxidoreductase increased in CmDOF18-overexpressing transgenic plants but decreased in CmDOF18-SRDX gene-silenced transgenic plants.

CONCLUSION

In summary, we analyzed the function of CmDOF18 from chrysanthemum, which may regulate salinity stress in plants, possibly due to its role in the regulation of oxidoreductase.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

A group VIIIa ethylene-responsive factor, CmERF4, negatively regulates waterlogging tolerance in chrysanthemum

Ethylene-responsive factors (ERF) play an important role in plant responses to waterlogging stress. However, the function and mechanism of action of ERFVIII in response to waterlogging stress remain poorly understood. In this study, we found that expression of the ERF VIIIa gene CmERF4 in chrysanthemum was induced by waterlogging stress. CmERF4 localized to the nucleus when expressed in tobacco leaves. Yeast two-hybrid and luciferase assays showed that CmERF4 is a transcriptional inhibitor. CmERF4 overexpression in chrysanthemum reduced plant waterlogging tolerance, whereas overexpression of the chimeric activator CmERF4-VP64 reversed its transcriptional activity, promoting higher waterlogging tolerance than that observed in wild-type plants, indicating that CmERF4 negatively regulates waterlogging tolerance. Transcriptome profiling showed that energy metabolism and reactive oxygen species (ROS) pathway-associated genes were differentially expressed between CmERF4-VP64 and wild-type plants. RT-qPCR analysis of selected energy metabolism and reactive oxygen species-related genes showed that the gene expression patterns were consistent with the expression levels obtained from RNA-seq analysis. Overall, we identified new functions of CmERF4 in negatively regulating chrysanthemum waterlogging tolerance by modulating energy metabolism and ROS pathway genes.

Genomic insights into the evolution of flavonoid biosynthesis and O-methyltransferase and glucosyltransferase in Chrysanthemum indicum

Flavonoids are a class of secondary metabolites widely distributed in plants. Regiospecific modification by methylation and glycosylation determines flavonoid diversity. A rare flavone glycoside, diosmin (luteolin-4'-methoxyl-7-O-glucosyl-rhamnoside), occurs in Chrysanthemum indicum. How Chrysanthemum plants evolve new biosynthetic capacities remains elusive. Here, we assemble a 3.11-Gb high-quality C. indicum genome with a contig N50 value of 4.39 Mb and annotate 50,606 protein-coding genes. One (CiCOMT10) of the tandemly repeated O-methyltransferase genes undergoes neofunctionalization, preferentially transferring the methyl group to the 4'-hydroxyl group of luteolin with ortho-substituents to form diosmetin. In addition, CiUGT11 (UGT88B3) specifically glucosylates 7-OH group of diosmetin. Next, we construct a one-pot cascade biocatalyst system by combining CiCOMT10, CiUGT11, and our previously identified rhamnosyltransferase, effectively producing diosmin with over 80% conversion from luteolin. This study clarifies the role of transferases in flavonoid diversity and provides important gene elements essential for producing rare flavone.

[Identification and expression of genome of uridine diphosphate glycosyltransferase (UGT) gene family from Chrysanthemum indicum]

Uridine diphosphate glycosyltransferase(UGT) is involved in the glycosylation of a variety of secondary metabolites in plants and plays an important role in plant growth and development and regulation of secondary metabolism. Based on the genome of a diploid Chrysanthemum indicum, the UGT gene family from Ch. indicum was identified by bioinformatics methods, and the physical and chemical properties, subcellular localization prediction, conserved motif, phylogeny, chromosome location, gene structure, and gene replication events of UGT protein were analyzed. Transcriptome and real-time fluorescence quantitative polymerase chain reaction(PCR) were used to analyze the expression pattern of the UGT gene in flowers and leaves of Ch. indicum. Quasi-targeted metabolomics was used to analyze the differential metabolites in flowers and leaves. The results showed that a total of 279 UGT genes were identified in the Ch. indicum genome. Phylogenetic analysis showed that these UGT genes were divided into 8 subfamilies. Members of the same subfamily were distributed in clusters on the chromosomes. Tandem duplications were the main driver of the expansion of the UGT gene family from Ch. indicum. Structural domain analysis showed that 262 UGT genes had complete plant secondary metabolism signal sequences(PSPG box). The analysis of cis-acting elements indicated that light-responsive elements were the most ubiquitous elements in the promoter regions of UGT gene family members. Quasi-targeted metabolome analysis of floral and leaf tissue revealed that most of the flavonoid metabolites, including luteolin-7-O-glucoside and kaempferol-7-O-glucoside, had higher accumulation in flowers. Comparative transcriptome analysis of flower and leaf tissue showed that there were 72 differentially expressed UGT genes, of which 29 genes were up-regulated in flowers, and 43 genes were up-regulated in leaves. Correlation network and phylogenetic analysis showed that CindChr9G00614970.1, CindChr2G00092510.1, and CindChr2G00092490.1 may be involved in the synthesis of 7-O-flavonoid glycosides in Ch. indicum, and real-time fluorescence quantitative PCR analysis further confirmed the reliability of transcriptome data. The results of this study are helpful to understand the function of the UGT gene family from Ch. indicum and provide data reference and theoretical basis for further study on the molecular regulation mechanism of flavonoid glycosides synthesis in Ch. indicum.

The histone deacetylase SRT2 enhances the tolerance of chrysanthemum to low temperatures through the ROS scavenging system

Low temperatures can severely affect plant growth and reduce their ornamental value. A family of plant histone deacetylases allows plants to cope with both biotic and abiotic stresses. In this study, we screened and cloned the cDNA of DgSRT2 obtained from transcriptome sequencing of chrysanthemum leaves under low-temperature stress. Sequence analysis showed that DgSRT2 belongs to the sirtuin family of histone deacetylases. We obtained the stable transgenic chrysanthemum lines OE-2 and OE-12. DgSRT2 showed tissue specificity in wild-type chrysanthemum and was most highly expressed in leaves. Under low-temperature stress, the OE lines showed higher survival rates, proline content, solute content, and antioxidant enzyme activities, and lower relative electrolyte leakage, malondialdehyde, hydrogen peroxide, and superoxide ion accumulation than the wild-type lines. This work suggests that DgSRT2 can serve as an essential gene for enhancing cold resistance in plants. In addition, a series of cold-responsive genes in the OE line were compared with WT. The results showed that DgSRT2 exerted a positive regulatory effect by up-regulating the transcript levels of cold-responsive genes. The above genes help to increase antioxidant activity, maintain membrane stability and improve osmoregulation, thereby enhancing survival under cold stress. It can be concluded from the above work that DgSRT2 enhances chrysanthemum tolerance to low temperatures by scavenging the ROS system.

The BTB/TAZ domain-containing protein CmBT1-mediated CmANR1 ubiquitination negatively regulates root development in chrysanthemum

Roots are fundamental for plants to adapt to variable environmental conditions. The development of a robust root system is orchestrated by numerous genetic determinants and, among them, the MADS-box gene ANR1 has garnered substantial attention. Prior research has demonstrated that, in chrysanthemum, CmANR1 positively regulates root system development. Nevertheless, the upstream regulators involved in the CmANR1-mediated regulation of root development remain unidentified. In this study, we successfully identified bric-a-brac, tramtrack and broad (BTB) and transcription adapter putative zinc finger (TAZ) domain protein CmBT1 as the interacting partner of CmANR1 through a yeast-two-hybrid (Y2H) screening library. Furthermore, we validated this physical interaction through bimolecular fluorescence complementation and pull-down assays. Functional assays revealed that CmBT1 exerted a negative influence on root development in chrysanthemum. In both in vitro and in vivo assays, it was evident that CmBT1 mediated the ubiquitination of CmANR1 through the ubiquitin/26S proteasome pathway. This ubiquitination subsequently led to the degradation of the CmANR1 protein and a reduction in the transcription of CmANR1-targeted gene CmPIN2, which was crucial for root development in chrysanthemum. Genetic analysis suggested that CmBT1 modulated root development, at least in part, by regulating the level of CmANR1 protein. Collectively, these findings shed new light on the regulatory role of CmBT1 in degrading CmANR1 through ubiquitination, thereby repressing the expression of its targeted gene and inhibiting root development in chrysanthemum.

Transcription factor DgMYB recruits H3K4me3 methylase to DgPEROXIDASE to enhance chrysanthemum cold tolerance

Cold affects the growth and development of plants. MYB transcription factors and histone H3K4me3 transferase ARABIDOPSIS TRITHORAXs (ATXs) play important regulatory functions in the process of plant resistance to low-temperature stress. In this study, DgMYB expression was responsive to low temperature, and overexpression of DgMYB led to increased tolerance, whereas the dgmyb mutant resulted in decreased tolerance of Chrysanthemum morifolium (Dendranthema grandiflorum var. Jinba) to cold stresses. Interestingly, we found that only peroxidase (POD) activity differed substantially between wild type (WT), overexpression lines, and the mutant line. A DgATX H3K4me3 methylase that interacts with DgMYB was isolated by further experiments. DgATX expression was also responsive to low temperature. Overexpression of DgATX led to increased tolerance, whereas the dgatx mutant resulted in decreased tolerance of chrysanthemum to cold stresses. Moreover, the dgmyb, dgatx, and dgmyb dgatx double mutants all led to reduced H3K4me3 levels at DgPOD, thus reducing DgPOD expression. Together, our results show that DgMYB interacts with DgATX, allowing DgATX to specifically target DgPOD, altering H3K4me3 levels, increasing DgPOD expression, and thereby reducing the accumulation of reactive oxygen species (ROS) in chrysanthemum.

Chrysanthemum (Chrysanthemum morifolium) CmHRE2-like negatively regulates the resistance of chrysanthemum to the aphid (Macrosiphoniella sanborni)

BACKGROUND

The growth and ornamental value of chrysanthemums are frequently hindered by aphid attacks. The ethylene-responsive factor (ERF) gene family is pivotal in responding to biotic stress, including insect stress. However, to date, little is known regarding the involvement of ERF transcription factors (TFs) in the response of chrysanthemum to aphids.

RESULTS

In the present study, CmHRE2-like from chrysanthemum (Chrysanthemum morifolium), a transcription activator that localizes mainly to the nucleus, was cloned. Expression is induced by aphid infestation. Overexpression of CmHRE2-like in chrysanthemum mediated its susceptibility to aphids, whereas CmHRE2-like-SRDX dominant repressor transgenic plants enhanced the resistance of chrysanthemum to aphids, suggesting that CmHRE2-like contributes to the susceptibility of chrysanthemum to aphids. The flavonoids in CmHRE2-like-overexpression plants were decreased by 29% and 28% in two different lines, whereas they were increased by 42% and 29% in CmHRE2-like-SRDX dominant repressor transgenic plants. The expression of Chrysanthemum-chalcone-synthase gene(CmCHS), chalcone isomerase gene (CmCHI), and flavonoid 3'-hydroxylase gene(CmF3'H) was downregulated in CmHRE2-like overexpression plants and upregulated in CmHRE2-like-SRDX dominant repressor transgenic plants, suggesting that CmHRE2-like regulates the resistance of chrysanthemum to aphids partially through the regulation of flavonoid biosynthesis.

CONCLUSION

CmHRE2-like was a key gene regulating the vulnerability of chrysanthemum to aphids. This study offers fresh perspectives on the molecular mechanisms of chrysanthemum-aphid interactions and may bear practical significance for developing new strategies to manage aphid infestation in chrysanthemums.

Complete nucleotide sequence of chrysanthemum virus D, a polero-like virus

A putative new polerovirus, named "chrysanthemum virus D" (ChVD), was detected in a Chrysanthemum morifolium plant in South Korea. The virus was identified by high-throughput sequencing and confirmed by reverse transcription polymerase chain reaction. The entire ChVD genome is composed of 5,963 nucleotides and contains seven open reading frames (ORF0-5 and ORF3a), which are arranged similarly to those of other poleroviruses. These ORFs encode the putative proteins P0-5 and P3a, respectively. Pairwise amino acid sequence comparisons showed that the ChVD P0-5 and P3a proteins have 30.45-75% sequence identity to the corresponding proteins of other members of the genus Polerovirus. Since one of the species demarcation criteria for the genus Polerovirus is > 10% difference in the amino acid sequence of any gene product, the sequence comparisons indicate that ChVD represents a new species in this genus. Phylogenetic analysis of the P1-P2 and P3 amino acid sequences further indicate that ChVD is a novel polerovirus.

DgHDA6 enhances the cold tolerance in chrysanthemum by improving ROS scavenging capacity

Histone deacetylases have been demonstrated to play an important role in responding to low-temperature stress, but the related response mechanism in chrysanthemum remains unclear. In this study, we isolated a cold-induced gene, DgHDA6, from chrysanthemum (Chrysanthemum morifolium Ramat). DgHDA6 contains 474 amino acids and shares a typical deacetylation domain with RPD3/HDA1 family members. The overexpression of DgHDA6 enhanced cold resistance in chrysanthemums. After low-temperature stress, the overexpression lines showed a higher survival rate. The contents of proline, soluble proteins and sugars, and the activities of antioxidant enzymes were significantly increased while the contents of H2O2, O2- and MDA were lower. Moreover, cold-stress-responding genes such as DgCuZnSOD, DgCAT, DgP5CS, and DgFAD were upregulated after cold stress. These results suggest that the overexpression of DgHDA6 can improve cold tolerance in chrysanthemum by enhancing ROS scavenging capacity.

Transcriptome analysis during axillary bud growth in chrysanthemum (chrysanthemum×morifolium)

The chrysanthemum DgLsL gene, homologous with tomato Ls, is one of the earliest expressed genes controlling axillary meristem initiation. In this study, the wild-type chrysanthemum (CW) and DgLsL-overexpressed line 15 (C15) were used to investigate the regulatory mechanism of axillary bud development in chrysanthemum. Transcriptome sequencing was carried out to detect the differentially expressed genes of the axillary buds 0 h, 24 h and 48 h after decapitation. The phenotypic results showed that the number of axillary buds of C15 was significantly higher than CW. A total of 9,224 DEGs were identified in C15-0 vs. CW-0, 10,622 DEGs in C15-24 vs. CW-24, and 8,929 DEGs in C15-48 vs. CW-48.GO and KEGG pathway enrichment analyses showed that the genes of the flavonoid, phenylpropanoids and plant hormone pathways appeared to be differentially expressed, indicating their important roles in axillary bud germination. DgLsL reduces GA content in axillary buds by promoting GA2ox expression.These results confirmed previous studies on axillary bud germination and growth, and revealed the important roles of genes involved in plant hormone biosynthesis and signal transduction, aiding in the study of the gene patterns involved in axillary bud germination and growth.

Sucrose-delaying flower color fading associated with delaying anthocyanin accumulation decrease in cut chrysanthemum

As fresh ornamental crops, vase life and post-harvested quality of cut flowers have attracted much attention. Flower color fading is the prominent defect in red and purple cut flowers, especially in cut chrysanthemum which have a relative long vase life. Here, the effect of sucrose on change in anthocyanin contents during the vase life of 'Dante Purple' cut chrysanthemum was studied. Results showed that 500 mM sucrose as holding solution could significantly delay the decrease in anthocyanin content and maintain the ornamental value for as long as 38 vase days. Moreover, the sucrose also increased the flower diameter, soluble sugar contents and total antioxidant capacity, while decreasing the malondialdehyde contents. Further studies suggested that the transcript levels of anthocyanin biosynthetic genes and transcription factors, CmMYB6 and CmMYB#7, had continuously decreased during the vase life. The changes in these genes expression patterns was retarded by the sucrose treatment, except for CmMYB#7 which is a repressor of anthocyanin biosynthesis gene expression. The decline in relative expression of CmMYB#7 was accelerated by sucrose. These results have supplied clues to study the mechanism whereby sucrose serves as a signal molecule to regulate anthocyanin biosynthesis.

Multi-locus genome-wide association study and genomic prediction for flowering time in chrysanthemum

MAIN CONCLUSION

Multi-locus GWAS detected several known and candidate genes responsible for flowering time in chrysanthemum. The associations could greatly increase the predictive ability of genome selection that accelerates the possible application of GS in chrysanthemum breeding. Timely flowering is critical for successful reproduction and determines the economic value for ornamental plants. To investigate the genetic architecture of flowering time in chrysanthemum, a multi-locus genome-wide association study (GWAS) was performed using a collection of 200 accessions and 330,710 single-nucleotide polymorphisms (SNPs) via 3VmrMLM method. Five flowering time traits including budding (FBD), visible colouring (VC), early opening (EO), full-bloom (OF) and senescing (SF) stages, plus five derived conditional traits were recorded in two environments. Extensive phenotypic variations were observed for these flowering time traits with coefficients of variation ranging from 6.42 to 38.27%, and their broad-sense heritability ranged from 71.47 to 96.78%. GWAS revealed 88 stable quantitative trait nucleotides (QTNs) and 93 QTN-by-environment interactions (QEIs) associated with flowering time traits, accounting for 0.50-8.01% and 0.30-10.42% of the phenotypic variation, respectively. Amongst the genes around these stable QTNs and QEIs, 21 and 10 were homologous to known flowering genes in Arabidopsis; 20 and 11 candidate genes were mined by combining the functional annotation and transcriptomics data, respectively, such as MYB55, FRIGIDA-like, WRKY75 and ANT. Furthermore, genomic selection (GS) was assessed using three models and seven unique marker datasets. We found the prediction accuracy (PA) using significant SNPs identified by GWAS under SVM model exhibited the best performance with PA ranging from 0.90 to 0.95. Our findings provide new insights into the dynamic genetic architecture of flowering time and the identified significant SNPs and candidate genes will accelerate the future molecular improvement of chrysanthemum.

Comparative Physiological and Transcriptome Analysis of Crossostephium chinense Reveals Its Molecular Mechanisms of Salt Tolerance

Crossostephium chinense is a wild species with strong salt tolerance that has great potential to improve the salt tolerance of cultivated chrysanthemums. Conversely, the unique salt-tolerant molecular mechanisms of Cr. chinense are still unclear. This study performed a comparative physiological and transcriptome analysis of Cr. chinense, Chrysanthemum lavandulifolium, and three hybrids to investigate the salt-tolerant molecular mechanisms of Cr. chinense. The physiological results showed that Cr. chinense maintained higher superoxide dismutase (SOD) activity, alleviating oxidative damage to the membrane. KEGG enrichment analysis showed that plant hormone signaling transduction and the MAPK signaling pathway were mostly enriched in Cr. chinense and hybrids under salt stress. Further weighted gene co-expression network analysis (WGCNA) of DEGs suggested that abscisic acid (ABA) signaling transduction may play a significant role in the salt-tolerant mechanisms of Cr. chinense and hybrids. The tissue-specific expression patterns of the candidate genes related to ABA signaling transduction and the MAPK signaling pathway indicate that genes related to ABA signaling transduction demonstrated significant expression levels under salt stress. This study offers important insights into exploring the underlying salt-tolerant mechanisms of Cr. chinense mediated by ABA signaling transduction and broadens our understanding of the breeding strategies for developing salt-tolerant cultivars utilizing salt-tolerant chrysanthemum germplasms.

PHOTOLYASE/BLUE LIGHT RECEPTOR2 regulates chrysanthemum flowering by compensating for gibberellin perception

The gibberellins (GAs) receptor GA INSENSITIVE DWARF1 (GID1) plays a central role in GA signal perception and transduction. The typical photoperiodic plant chrysanthemum (Chrysanthemum morifolium) only flowers when grown in short-day photoperiods. In addition, chrysanthemum flowering is also controlled by the aging pathway, but whether and how GAs participate in photoperiod- and age-dependent regulation of flowering remain unknown. Here, we demonstrate that photoperiod affects CmGID1B expression in response to GAs and developmental age. Moreover, we identified PHOTOLYASE/BLUE LIGHT RECEPTOR2, an atypical photocleavage synthase, as a CRYPTOCHROME-INTERACTING bHLH1 interactor with which it forms a complex in response to short days to activate CmGID1B transcription. Knocking down CmGID1B raised endogenous bioactive GA contents and GA signal perception, in turn modulating the expression of the aging-related genes MicroRNA156 and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE3. We propose that exposure to short days accelerates the juvenile-to-adult transition by increasing endogenous GA contents and response to GAs, leading to entry into floral transformation.

Flowering repressor CmSVP recruits the TOPLESS corepressor to control flowering in chrysanthemum

Plant flowering time is induced by environmental and endogenous signals perceived by the plant. The MCM1-AGAMOUSDEFICIENS-Serum Response Factor-box (MADS-box) protein SHORT VEGETATIVE PHASE (SVP) is a pivotal repressor that negatively regulates the floral transition during the vegetative phase; however, the transcriptional regulatory mechanism remains poorly understood. Here, we report that CmSVP, a chrysanthemum (Chrysanthemum morifolium Ramat.) homolog of SVP, can repress the expression of a key flowering gene, a chrysanthemum FLOWERING LOCUS T-like gene (CmFTL3), by binding its promoter CArG element to delay flowering in the ambient temperature pathway in chrysanthemum. Protein-protein interaction assays identified an interaction between CmSVP and CmTPL1-2, a chrysanthemum homologue of TOPLESS (TPL) that plays critical roles as transcriptional corepressor in many aspects of plant life. Genetic analyses revealed the CmSVP-CmTPL1-2 transcriptional complex is a prerequisite for CmSVP to act as a floral repressor. Furthermore, overexpression of CmSVP rescued the phenotype of the svp-31 mutant in Arabidopsis (Arabidopsis thaliana), overexpression of AtSVP or CmSVP in the Arabidopsis dominant-negative mutation tpl-1 led to ineffective late flowering, and AtSVP interacted with AtTPL, confirming the conserved function of SVP in chrysanthemum and Arabidopsis. We have validated a conserved machinery wherein SVP partially relies on TPL to inhibit flowering via a thermosensory pathway.

CmNAC25 targets CmMYB6 to positively regulate anthocyanin biosynthesis during the post-flowering stage in chrysanthemum

BACKGROUND

Anthocyanin is a class of important secondary metabolites that determines colorful petals in chrysanthemum, a famous cut flower. 'Arctic Queen' is a white chrysanthemum cultivar that does not accumulate anthocyanin during the flowering stage. During the post-flowering stage, the petals of 'Arctic Queen' accumulate anthocyanin and turn red. However, the molecular mechanism underlying this flower color change remains unclear.

RESULTS

In this study, by using transcriptome analysis, we identified CmNAC25 as a candidate gene promoting anthocyanin accumulation in the post-flowering stage of 'Arctic Queen'. CmNAC25 is directly bound to the promoter of CmMYB6, a core member of the MBW protein complex that promotes anthocyanin biosynthesis in chrysanthemum, to activate its expression. CmNAC25 also directly activates the promoter of CmDFR, which encodes the key enzyme in anthocyanin biosynthesis. CmNAC25 was highly expressed during the post-flowering stage, while the expression level of CmMYB#7, a known R3 MYB transcription factor interfering with the formation of the CmMYB6-CmbHLH2 complex, significantly decreased. Genetic transformation of both chrysanthemum and Nicotiana tabacum verified that CmNAC25 was a positive regulator of anthocyanin biosynthesis. Another two cultivars that turned red during the post-flowering stages also demonstrated a similar mechanism.

CONCLUSIONS

Altogether, our data revealed that CmNAC25 positively regulates anthocyanin biosynthesis in chrysanthemum petals during the post-flowering stages by directly activating CmMYB6 and CmDFR. Our results thus revealed a crucial role of CmNAC25 in regulating flower color change during petal senescence and provided a target gene for molecular design breeding of flower color in chrysanthemum.

Jasmonate signaling drives defense responses against Alternaria alternata in chrysanthemum

BACKGROUND

Black spot disease caused by the necrotrophic fungus Alternaria spp. is one of the most devastating diseases affecting Chrysanthemum morifolium. There is currently no effective way to prevent chrysanthemum black spot.

RESULTS

We revealed that pre-treatment of chrysanthemum leaves with the methy jasmonate (MeJA) significantly reduces their susceptibility to Alternaria alternata. To understand how MeJA treatment induces resistance, we monitored the dynamics of metabolites and the transcriptome in leaves after MeJA treatment following A. alternata infection. JA signaling affected the resistance of plants to pathogens through cell wall modification, Ca2+ regulation, reactive oxygen species (ROS) regulation, mitogen-activated protein kinase cascade and hormonal signaling processes, and the accumulation of anti-fungal and anti-oxidant metabolites. Furthermore, the expression of genes associated with these functions was verified by reverse transcription quantitative PCR and transgenic assays.

CONCLUSION

Our findings indicate that MeJA pre-treatment could be a potential orchestrator of a broad-spectrum defense response that may help establish an ecologically friendly pest control strategy and offer a promising way of priming plants to induce defense responses against A. alternata.

[Analysis of the Complete Tomato Aspermy Virus Genomes Suggests Reassortment in Russian Isolates from Chrysanthemum]

Tomato aspermy virus (TAV, genus Cucumovirus from the family Bromoviridae) is one of the most common and harmful chrysanthemum viruses, causing severe flower distortion, size reduction, and color breaking. Metatranscriptome sequencing of chrysanthemum plants of the Ribonette and Golden Standard cultivars from the collection of the Nikita Botanical Garden (Yalta, Republic of Crimea) generated TAV-related RNA reads. The complete genomes of two Russian isolates of the virus were assembled from the reads. This is the first report of full-length TAV genomes from Russia. Typically of cucumoviruses, the segmented TAV genome is represented by three single-stranded positive-sense linear RNA molecules of 3412 (RNA1), 3097 (RNA2) and 2219 (RNA3) nucleotides. Five open reading frames (ORF) have been identified that encode replicase (ORF1), RNA-dependent RNA polymerase (ORF2a), silencing suppressor protein (OFR2b), movement protein (OFR3a) and the coat protein (ORF3b). The identity of TAV genomes from the two chrysanthemum cultivars was 99.8% for all three viral RNAs; with other TAV isolates from GenBank it was 97.5-99.7% (RNA1), 93.8-99.8% (RNA2), and 89.3-99.3% (RNA3). Phylogenetic analysis showed that RNA1 and RNA3 of the Russian isolates were assigned to heterogeneous groups of TAV isolates found on various plant species in different regions of the world. At the same time, RNA2 clearly clustered with tomato isolates SKO20ST2 from Slovenia and PV-0220 from Bulgaria and, to a lesser extent, with the Iranian isolate Ker.Mah.P from petunia and the Chinese isolate Henan from chrysanthemum. The incongruence of phylogenetic trees reconstructed from different genome segments suggests pseudo-recombination (reassortment) in the Russian TAV isolates.

Assembly and co-occurrence patterns of rhizosphere bacterial communities are closely linked to soil fertility during continuous cropping of cut chrysanthemum (Chrysanthemum morifolium Ramat)

AIMS

Continuous cropping is known to have profound effects on the soil microbial community in different planting systems. However, we lack an understanding of how different years of continuous cropping affects rhizosphere soil bacterial community co-occurrence pattern and assembly processes in the cut chrysanthemum (Chrysanthemum morifolium Ramat.) field.

METHODS AND RESULTS

We collected the soils from cut chrysanthemum rhizospheres with planting for 1 year (PY1) and continuous cropping for 6 years (CY6) and 12 years (CY12). Real-time quantitative PCR and flow cytometry (FCM) techniques were used to test the 16S rRNA gene copy number and bacterial cell count, respectively. The bacterial community structure was analysed by using high-throughput sequencing technology. The CY12 had a significantly decreased soil fertility index and rhizosphere bacterial living cell counts and gene copy numbers compared to CY6 and PY1 (P < 0.05). The rhizosphere bacterial community dissimilarity increased as the continuous cropping years increased. Three main ecological clusters (modules #1, #2, and #3) were observed in the bacterial co-occurrence network across all samples, and only the relative abundance of module #1 (enriched in the CY12) was significantly correlated with soil fertility (P < 0.05). Moreover, the rhizosphere bacterial community assembly was primarily governed by the deterministic process under 12 years of continuous cropping.

CONCLUSIONS

Soil fertility decline correlates with ecological network modularization and the deterministic assembly process of the rhizosphere bacterial community of cut chrysanthemum during continuous cropping.

[Genetic diversity analysis and fingerprints of Chrysanthemum×morifolium based on SSR molecular markers]

The present study aims to explore the genetic diversity of germplasm resources of Chrysanthemum×morifolium (hereinafter, C.×morifolium) at the molecular level and to establish a fingerprint database of C.×morifolium varieties. We employed 12 pairs of primers with high levels of polymorphism, clear bands, and high degrees of reproducibility to analyze the SSR molecular markers and genetic diversity of 91 C.×morifolium materials and 14 chrysanthemum- related materials. With regard to constructing the fingerprints of the tested materials, we chose 9 pairs of core primers. The findings revealed that 12 primer pairs detected 104 alleles in 105 samples, ranging from 2 to 26. The average number of observed alleles (Na) per site was 9.25. The average number of effective alleles (Ne) per site was 2.745 6, with its range being 1.276 0 to 4.742 5. Shannon genetic diversity index (I) values ranged between 0.513 3 and 2.239 9 (M=1.209 0). Nei's gene diversity index (H) ranged between 0.216 3 and 0.789 1 (M=0.578 0). The observed heterozygosity (Ho) ranged between 0.223 3 and 0.895 2 (M=0.557 5). The expected heterozygosity (He) ranged between 0.217 4 and 0.793 3 (M=0.580 8). The polymorphism information content (PIC) ranged between 0.211 5 and 0.774 0 (M=0.532 9). The genetic similarity (GS) ranged between 0.228 5 and 1.000 0 (M=0.608 3). Cluster analysis revealed that when the genetic distance (GD) equals to 0.30, the tested materials can be classified into 2 groups. When the GD equals to 0.27, the first group can be divided into 6 subgroups; accordingly, 105 tested materials can be divided into 7 subgroups. The cophenetic correlation test was carried out based on the cluster analysis, and the corresponding results showed that the cluster map correlated with the genetic similarity coefficient (r=0.952 73). According to the results of Structure population analysis, we obtained the optimal population number, with the true number of populations (K) being 3 and the population being divided concerning Q≥0.5. Three subgroups, i.e., Q1, Q2 and Q3, included 34, 33 and 28 germplasms, respectively, and the remaining 10 germplasms were identified as the mixed population. During the experiment, 9 pairs of core primers were screened among the total of 12 for a complete differentiation regarding 105 tested materials, and the fingerprints of 91 C.×morifolium materials and 14 chrysanthemum-related materials were further constructed. Overall, there were significant genetic differences and rich genetic diversity among C.×morifolium materials, which would shed light on the garden application and variety selection fields of C.×morifolium. The fingerprint database of 105 C.×morifolium varieties and chrysanthemum-related species may provide technical support for future research regarding the identification and screening system of C.×morifolium varieties.

The CmMYB3 transcription factors isolated from the Chrysanthemum morifolium regulate flavonol biosynthesis in Arabidopsis thaliana

Chrysanthemum morifolium MYB3 factors are transcriptional activators for the regulation of flavonol biosynthesis. Flavonol was not only the critical secondary metabolite participating in the growth and development of plants but also the main active ingredient in medicinal chrysanthemum. However, few pieces of research revealed the transcriptional regulation of flavonol biosynthesis in Chrysanthemum morifolium. Here, we isolated two CmMYB3 transcription factors (CmMYB3a and CmMYB3b) from the capitulum of Chrysanthemum morifolium cv 'Hangju'. According to the sequence characteristics, the CmMYB3a and CmMYB3b belonged to the R2R3-MYB subgroup 7, whose members were often reported to regulate flavonol biosynthesis positively. CmMYB3a and CmMYB3b factors were identified to localize in the nucleus by subcellular localization assay. Besides, both of them have obvious transcriptional self-activation activity in their C-terminal. After the overexpression of CmMYB3 genes in Nicotiana benthamiana and Arabidopsis thaliana, the flavonol contents in plants were increased, and the expression of AtCHS, AtCHI, AtF3H, and AtFLS genes in A. thaliana was also improved. Interestingly, the CmMYB3a factor had stronger functions in improving flavonol contents and related gene expression levels than CmMYB3b. The interaction analysis between transcription factors and promoters suggested that CmMYB3 could bind and activate the promoters of CmCHI and CmFLS genes in C. morifolium, and CmMYB3a also functioned more powerfully. Overall, these results indicated that CmMYB3a and CmMYB3b work as transcriptional activators in controlling flavonol biosynthesis.

Genome-wide characterization and expression analysis of MYB transcription factors in Chrysanthemum nankingense

BACKGROUND

Chrysanthemum is a popular ornamental plant worldwide. MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factors play an important role in everything from stress resistance to plant growth and development. However, the MYB family of chrysanthemums has not been the subject of a detailed bioinformatics and expression investigation.

RESULTS

In this study, we examined 324 CnMYB transcription factors from Chrysanthemum nankingense genome data, which contained 122 Cn1R-MYB, 183 CnR2R3-MYB, 12 Cn3R-MYB, 2 Cn4R-MYB, and 5 atypical CnMYB. The protein motifs and classification of CnMYB transcription factors were analyzed. Among them, motifs 1, 2, 3, and 4 were found to encode the MYB DNA-binding domain in R2R3-MYB proteins, while in other-MYB proteins, the motifs 1, 2, 3, 4, 5, 6, 7, and 8 encode the MYB DNA-binding domain. Among all CnMYBs, 44 genes were selected due to the presence of CpG islands, while methylation is detected in three genes, including CnMYB9, CnMYB152, and CnMYB219. We analyzed the expression levels of each CnMYB gene in ray floret, disc floret, flower bud, leaf, stem, and root tissues. Based on phylogenetic analysis and gene expression analysis, three genes appeared likely to control cellulose and lignin synthesis in stem tissue, and 16 genes appeared likely to regulate flowering time, anther, pollen development, and flower color. Fifty-one candidate genes that may be involved in stress response were identified through phylogenetic, stress-responseve motif of promoter, and qRT-PCR analyses. According to genes expression levels under stress conditions, six CnMYB genes (CnMYB9, CnMYB172, CnMYB186, CnMYB199, CnMYB219, and CnMYB152) were identified as key stress-responsive genes.

CONCLUSIONS

This research provides useful information for further functional analysis of the CnMYB gene family in chrysanthemums, as well as offers candidate genes for further study of cellulose and lignin synthesis, flowering traits, salt and drought stress mechanism.

Molecular module of CmMYB15-like-Cm4CL2 regulating lignin biosynthesis of chrysanthemum (Chrysanthemum morifolium) in response to aphid (Macrosiphoniella sanborni) feeding

Lignin is a major component of plant cell walls and a conserved basic defense mechanism in higher plants deposited in response to aphid infection. However, the molecular mechanisms of lignin biosynthesis in response to aphid infection and the effect of lignin on aphid feeding behavior remain unclear. We report that 4-Coumarate:coenzyme A ligase 2 (Cm4CL2), a gene encoding a key enzyme in the lignin biosynthesis pathway, is induced by aphid feeding, resulting in lignin deposition and reduced aphid attack. Upstream regulator analysis showed that the expression of Cm4CL2 in response to aphid feeding was directly upregulated by CmMYB15-like, an SG2-type R2R3-MYB transcription factor. CmMYB15-like binds directly to the AC cis-element in the promoter region of Cm4CL2. Genetic validation demonstrated that CmMYB15-like was induced by aphid infection and contributed to lignin deposition and cell wall thickening, which consequently enhanced aphid resistance in a Cm4CL2-dependent manner. This study is the first to show that the CmMYB15-like-Cm4CL2 module regulates lignin biosynthesis in response to aphid feeding.

CmWRKY41 activates CmHMGR2 and CmFPPS2 to positively regulate sesquiterpenes synthesis in Chrysanthemum morifolium

Chrysanthemum morifolium is one of the most significant multipurpose crops with ornamental, medicinal, and edible value. Terpenoids, an essentials component of volatile oils, are abundant in chrysanthemum. However, the transcriptional regulation of terpenoid biosynthesis in chrysanthemums remains unclear. In the present investigation, we identified CmWRKY41, whose expression pattern is similar to that of terpenoid content in chrysanthemum floral scent, as a candidate gene that may promote terpenoid biosynthesis in chrysanthemum. Two structural genes 3-hydroxy-3-methylglutaryl-CoA reductase 2 (CmHMGR2) and farnesyl pyrophosphate synthase 2 (CmFPPS2), play key role in terpene biosynthesis in chrysanthemum. CmWRKY41 can directly bind to the promoters of CmHMGR2 or CmFPPS2 through GTGACA or CTGACG elements and activate its expression to promote sesquiterpene biosynthesis. In summary, these results indicate that CmWRKY41 targets CmHMGR2 and CmFPPS2 to positively regulate sesquiterpene biosynthesis in chrysanthemums. This study preliminarily revealed the molecular mechanism of terpenoid biosynthesis in chrysanthemum while enriching the secondary metabolism regulatory network.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes

BACKGROUND AND AIMS

Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution.

METHODS

Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs.

KEY RESULTS

Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements.

CONCLUSIONS

The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

CmERF5-CmRAP2.3 transcriptional cascade positively regulates waterlogging tolerance in Chrysanthemum morifolium

Waterlogging stress affects plant growth by limiting root respiration and reducing yield and economic value. Therefore, identifying genes involved in regulating waterlogging stress is vital. This study reports the ethylene-responsive VII transcription factor (CmRAP2.3) in the chrysanthemum. Subcellular localization and transactivation assay analyses revealed that CmRAP2.3 was localized in the nucleus and possessed transactivation activity. Overexpression of CmRAP2.3 in chrysanthemum was found to enhance waterlogging tolerance by decreasing reactive oxygen species (ROS) levels. Furthermore, we found that the transcription factor CmERF5 binds to GCC-like motifs in the CmRAP2.3 promoter region and activates CmRAP2.3 expression. Additionally, CmERF5 overexpression maintained a low ROS level and improved chrysanthemum waterlogging tolerance. Taken together, this study shows a molecular mechanism by which CmERF5 transcriptionally activates CmRAP2.3 to reduce waterlogging stress via the ROS pathway in the chrysanthemum.

An ethylene-responsive transcription factor and a B-box protein coordinate vegetative growth and photoperiodic flowering in chrysanthemum

Plants employ several endogenous and exogenous signals to guarantee timely floral transitions with floral integrators. To avoid premature flowering, flowering plants must control the balance between vegetative and floral development. As a Group II member of BBX family, CmBBX8 promotes flowering by directly activating CmFTL1 in summer-flowering chrysanthemum. However, the mechanisms underlying this floral transition is yet to be elucidated. Here, we report that the chrysanthemum ERF3 homologue, CmERF3, physically interacts with CmBBX8 through yeast two-hybrid (Y2H), bimolecular fluorescence complementation (BiFC), pull-down, and luciferase complementation (LCI) assays. We found that CmERF3 was highly expressed at the vegetative stage and rarely expressed in the reproductive phase, indicating that CmERF3 may play a critical role in maintaining vegetative growth to prevent premature flowering. Rhythm analysis revealed that CmERF3 had a different response to rhythm compared to CmBBX8. Knockdown of CmERF3 facilitated floral initiation, whereas overexpression of CmERF3 delayed floral transition. We further found that CmERF3 repressed the transactivation activity of CmBBX8 on the downstream CmFTL1 gene. Collectively, our results indicate that the CmERF3-CmBBX8 transcriptional complex is a crucial module that balances the vegetative growth and reproductive development of chrysanthemum.

The TIFY family protein CmJAZ1-like negatively regulates petal size via interaction with the bHLH transcription factor CmBPE2 in Chrysanthemum morifolium

Petals are the second floral whorl of angiosperms, exhibiting astonishing diversity in their size between and within species. This variation is essential for protecting their inner reproductive organs and attracting pollinators for fertilization. However, currently, the genetic and developmental control of petal size remains unexplored. Chrysanthemum (Chrysanthemum morifolium) belongs to the Asteraceae family, the largest group of angiosperms, and the extraordinary diversity of petal size in chrysanthemums makes it an ideal model for exploring the regulation mechanism of petal size. Here, we reveal that overexpression of a JAZ repressor CmJAZ1-like exhibits decreased petal size compared to that of the wild-type as a result of repressed cell expansion. Through further in-depth exploration, we confirm an interaction pair between CmJAZ1-like and the bHLH transcription factor CmBPE2. The inhibition of CmBPE2 expression negatively regulates petal size by downregulating the expression of genes involved in cell expansion. Furthermore, CmJAZ1-like significantly reduced the activation ability of CmBPE2 on its target gene CmEXPA7 by directly interacting with it, thus participating in the regulation of petal size development in chrysanthemum. Our results will provide insights into the molecular mechanisms of petal size regulation in flowering plants.

DWARF AND ROBUST PLANT regulates plant height via modulating gibberellin biosynthesis in chrysanthemum

YABBY (YAB) genes are specifically expressed in abaxial cells of lateral organs and determine abaxial cell fate. However, most studies have focused on few model plants, and the molecular mechanisms of YAB genes are not well understood. Here, we identified a YAB transcription factor in chrysanthemum (Chrysanthemum morifolium), Dwarf and Robust Plant (CmDRP), that belongs to a distinct FILAMENTOUS FLOWER (FlL)/YAB3 sub-clade lost in Brassicaceae. CmDRP was expressed in various tissues but did not show any polar distribution in chrysanthemum. Overexpression of CmDRP resulted in a semi-dwarf phenotype with a significantly decreased active GA3 content, while reduced expression generated the opposite phenotype. Furthermore, plant height of transgenic plants was partially rescued through the exogenous application of GA3 and Paclobutrazol, and expression of the GA biosynthesis gene CmGA3ox1 was significantly altered in transgenic plants. Yeast one-hybrid, luciferase, and chromatin immunoprecipitation-qPCR analyses showed that CmDRP could directly bind to the CmGA3ox1 promoter and suppress its expression. Our research reveals a nonpolar expression pattern of a YAB family gene in dicots and demonstrates it regulates plant height through the GA pathway, which will deepen the understanding of the genetic and molecular mechanisms of YAB genes.

Rhamnosyltransferases involved in the biosynthesis of flavone rutinosides in Chrysanthemum species

Linarin (acacetin-7-O-rutinoside), isorhoifolin (apigenin-7-O-rutinoside), and diosmin (diosmetin-7-O-rutinoside) are chemically and structurally similar flavone rutinoside (FR) compounds found in Chrysanthemum L. (Anthemideae, Asteraceae) plants. However, their biosynthetic pathways remain largely unknown. In this study, we cloned and compared FRs and genes encoding rhamnosyltransferases (RhaTs) among eight accessions of Chrysanthemum polyploids. We also biochemically characterized RhaTs of Chrysanthemum plants and Citrus (Citrus sinensis and Citrus maxima). RhaTs from these two genera are substrate-promiscuous enzymes catalyzing the rhamnosylation of flavones, flavanones, and flavonols. Substrate specificity analysis revealed that Chrysanthemum 1,6RhaTs preferred flavone glucosides (e.g. acacetin-7-O-glucoside), whereas Cs1,6RhaT preferred flavanone glucosides. The nonsynonymous substitutions of RhaTs found in some cytotypes of diploids resulted in the loss of catalytic function. Phylogenetic analysis and specialized pathways responsible for the biosynthesis of major flavonoids in Chrysanthemum and Citrus revealed that rhamnosylation activity might share a common evolutionary origin. Overexpression of RhaT in hairy roots resulted in 13-, 2-, and 5-fold increases in linarin, isorhoifolin, and diosmin contents, respectively, indicating that RhaT is mainly involved in the biosynthesis of linarin. Our findings not only suggest that the substrate promiscuity of RhaTs contributes to the diversity of FRs in Chrysanthemum species but also shed light on the evolution of flavone and flavanone rutinosides in distant taxa.

Genome-Wide Identification, Characterization, and Expression Analysis of CHS Gene Family Members in Chrysanthemum nankingense

The chalcone synthase (CHS) gene family catalyzes the first committed step in the biosynthesis of flavonoids and plays key roles in various biological processes in plants. However, systematic studies of the CHS gene family in chrysanthemum remain unknown to date. In this study, 16 CnCHS genes were identified by searching the complete genome sequence of Chrysanthemum nankingense. Most contained two exons and one intron with Chal-sti-synt_N and Chal-sti-synt_C domains. A phylogenetic tree of CnCHSs indicated divergence into three major groups, including I, II, and III. Analyses of the genes and promoters of these genes indicated that there are many cis-acting elements that respond to light, phytohormones, stress, and developmental stages. The CnCHS genes have extensive patterns of expression in various tissues and stages of flower development. Tandemly repeated and segmental repeat genes were expressed at higher levels in different tissues than most of the CnCHS genes that have been identified. CnCHS10 is expressed at higher levels in various flower organs than in vegetative tissues, particularly in disc floret petals and pistils. Our study provides valuable information for the systematic analysis of the CnCHS gene family, which also contributes to further research on flavonoid synthesis and petal colors of chrysanthemum.

Phenotypic differences among and within extant populations of Chrysanthemum arcticum L. and C. a. subsp. arcticum

BACKGROUND

Chrysanthemum arcticum, arctic daisy and its two subspecies (Chrysanthemum arcticum subsp. arcticum, Chrysanthemum arcticum subsp. polaré) are the only chrysanthemum species native to North America. A study on species' variation in morphological and diagnostic traits is important to link morphological traits with previously described single nucleotide polymorphism (SNP) markers, particularly when the genomes are sequenced. The purpose of this study was to establish phenotypic differences and soil conditions among wild C. arcticum and C. a. subsp. arcticum populations, when grown in a uniform environment for two years, for potential linkages with our SNP library. Sixteen quantitative morphological traits and five qualitative morphological traits were investigated for 255 individuals from nine C. arcticum populations and 326 individuals from 21 C. a. subsp. arcticum populations.

RESULTS

In long-day controlled environment, C. arcticum flowering rate was 0% in Year 1, increased to 2.7% in Year 2, while C. a. subsp. arcticum flowering rate was 98.5% in Year 2. Two distinct clusters, distributed by taxonomic classification, were detected by Principal component analysis (PCoA) for 551 individuals from C. arcticum and C. a. subsp. arcticum. Pearson's correlation coefficient analysis indicated a positive and significant correlation between plant height, flower fresh and dry weights. Flower fresh weights were correlated with Δflower weight, while inflorescence length had showed a negative correlation with leaf number. Soil samples had high Na levels along with heavy metals. Thus, the species are salt-tolerant.

CONCLUSION

A high level of salt tolerance (Na) is tolerated by these maritime species which is a unique trait in Chrysanthemum. A new diagnostic trait of inflorescence length was discovered to distinguish among C. arcticum and C. a. subsp. arcticum. Significant flowering differences occurred among the species C. arcticum and C. a. subsp. arcticum under same photoperiodic environment, including flowering rates and visible bud date. This study on the species' variation in morphological and diagnostic traits is of importance to link morphological traits with single nucleotide polymorphism (SNP) markers.

Whole-transcriptome profiles of Chrysanthemum seticuspe improve genome annotation and shed new light on mRNA-miRNA-lncRNA networks in ray florets and disc florets

BACKGROUND

Chrysanthemum seticuspe has emerged as a model plant species of cultivated chrysanthemums, especially for studies involving diploid and self-compatible pure lines (Gojo-0). Its genome was sequenced and assembled into chromosomes. However, the genome annotation of C. seticuspe still needs to be improved to elucidate the complex regulatory networks in this species.

RESULTS

In addition to the 74,259 mRNAs annotated in the C. seticuspe genome, we identified 18,265 novel mRNAs, 51,425 novel lncRNAs, 501 novel miRNAs and 22,065 novel siRNAs. Two C-class genes and YABBY family genes were highly expressed in disc florets, while B-class genes were highly expressed in ray florets. A WGCNA was performed to identify the hub lncRNAs and mRNAs in ray floret- and disc floret-specific modules, and CDM19, BBX22, HTH, HSP70 and several lncRNAs were identified. ceRNA and lncNAT networks related to flower development were also constructed, and we found a latent functional lncNAT-mRNA combination, LXLOC_026470 and MIF2.

CONCLUSIONS

The annotations of mRNAs, lncRNAs and small RNAs in the C. seticuspe genome have been improved. The expression profiles of flower development-related genes, ceRNA networks and lncNAT networks were identified, laying a foundation for elucidating the regulatory mechanisms underlying disc floret and ray floret formation.

Co-option of a carotenoid cleavage dioxygenase gene (CCD4a) into the floral symmetry gene regulatory network contributes to the polymorphic floral shape-color combinations in Chrysanthemum sensu lato

Morphological novelties, including formation of trait combinations, may result from de novo gene origination and/or co-option of existing genes into other developmental contexts. A variety of shape-color combinations of capitular florets occur in Chrysanthemum and its allies. We hypothesized that co-option of a carotenoid cleavage dioxygenase gene into the floral symmetry gene network would generate a white zygomorphic ray floret. We tested this hypothesis in an evolutionary context using species in Chrysanthemum sensu lato, a monophyletic group with diverse floral shape-color combinations, based on morphological investigation, interspecific crossing, molecular interaction and transgenic experiments. Our results showed that white color was significantly associated with floret zygomorphy. Specific expression of the carotenoid cleavage dioxygenase gene CCD4a in marginal florets resulted in white color. Crossing experiments between Chrysanthemum lavandulifolium and Ajania pacifica indicated that expression of CCD4a is trans-regulated. The floral symmetry regulator CYC2g can activate expression of CCD4a with a dependence on TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING (TCP) binding element 8 on the CCD4a promoter. Based on all experimental findings, we propose that gene co-option of carotenoid degradation into floral symmetry regulation, and the subsequent dysfunction or loss of either CYC2g or CCD4a, may have led to evolution of capitular shape-color patterning in Chrysanthemum sensu lato.

Mitotically heritable epigenetic modifications of CmMYB6 control anthocyanin biosynthesis in chrysanthemum

Flower color, which is determined by various chemical pigments, is a vital trait for ornamental plants, in which anthocyanin is a major component. However, the epigenetic regulation of anthocyanin biosynthesis remains poorly understood. During chrysanthemum cultivation, we found a heterochromatic chrysanthemum accession (YP) whose progeny generated by asexual reproduction contained both yellow-flowered (YP-Y) and pink-flowered (YP-P) plants. In this study, we aimed to elucidate the epigenetic mechanisms of different flower colors in the YP plant progeny. Metabolome and transcriptome analyses revealed that the difference in flower color between YP-Y and YP-P was caused by expression variation of the anthocyanin biosynthesis gene CmMYB6. Bisulfite sequencing revealed that methylation at the CmMYB6 promoter, especially in the CHH context, was higher in YP-Y than YP-P. After demethylation of the CmMYB6 promoter using the dCas9-TET1cd system, the flower color returned from yellow to pink. Furthermore, the methylation status of the CmMYB6 promoter was higher in YP-Y over three consecutive generations, indicating that this methylation status was heritable mitotically. Finally, investigation of other chrysanthemum cultivars showed that the methylation of CmMYB6 decreased gradually with the increase in anthocyanin content. These results lay an epigenetic foundation for the improvement of flower color in horticultural plants.

Towards the Improvement of Ornamental Attributes in Chrysanthemum: Recent Progress in Biotechnological Advances

Incessant development and introduction of novel cultivars with improved floral attributes are vital in the dynamic ornamental industry. Chrysanthemum (Chrysanthemum morifolium) is a highly favored ornamental plant, ranking second globally in the cut flower trade, after rose. Development of new chrysanthemum cultivars with improved and innovative modifications in ornamental attributes, including floral color, shape, plant architecture, flowering time, enhanced shelf life, and biotic and abiotic stress tolerance, is a major goal in chrysanthemum breeding. Despite being an economically important ornamental plant, the application of conventional and molecular breeding approaches to various key traits of chrysanthemum is hindered owing to its genomic complexity, heterozygosity, and limited gene pool availability. Although classical breeding of chrysanthemum has resulted in the development of several hundreds of cultivars with various morphological variations, the genetic and transcriptional control of various important ornamental traits remains unclear. The coveted blue colored flowers of chrysanthemums cannot be achieved through conventional breeding and mutation breeding due to technical limitations. However, blue-hued flower has been developed by genetic engineering, and transgenic molecular breeding has been successfully employed, leading to substantial progress in improving various traits. The recent availability of whole-genome sequences of chrysanthemum offers a platform to extensively employ MAS to identify a large number of markers for QTL mapping, and GWAS to dissect the genetic control of complex traits. The combination of NGS, multi-omic platforms, and genome editing technologies has provided a tremendous scope to decipher the molecular and regulatory mechanisms. However, the application and integration of these technologies remain inadequate for chrysanthemum. This review, therefore, details the significance of floral attributes, describes the efforts of recent advancements, and highlights the possibilities for future application towards the improvement of crucial ornamental traits in the globally popular chrysanthemum plant.

Development of a Transformation System and Locus Identification Pipeline for T-DNA in Chrysanthemum seticuspe, A Model Species for Hexaploid Cultivated Chrysanthemum

Chrysanthemum is one of the most popular flowers worldwide and has high aesthetic and commercial value. However, the cultivated varieties of chrysanthemum are hexaploid and highly heterozygous, which makes gene editing and gene function research difficult. Gojo-0 is a diploid homozygous line bred from a self-compatible mutant of Chrysanthemum seticuspe and is expected to become a model plant of the genus Chrysanthemum. After assessment of different growth regulator combinations, the optimal concentrations of α-naphthaleneacetic acid (NAA) and 6-benzyladenine (6-BA) in the regeneration system were 1.0 mg·L−1 and 0.2 mg·L−1, respectively. In the genetic transformation system, the selected concentrations of kanamycin, hygromycin and glufosinate-ammonium were 10 mg·L−1, 2.5 mg·L−1 and 0.6 mg·L−1 for bud generation and 12 mg L−1, 1.5 mg·L−1 and 0.5 mg·L−1 for rooting. The transgenic plants were verified by not only PCR detection and GUS staining, but also identification of the T-DNA insertion locus using high-throughput sequencing. Our results lay the foundation for gene functional research on chrysanthemum and will help with the identification of transgenic plants.

Transcription factor CmbHLH16 regulates petal anthocyanin homeostasis under different lights in Chrysanthemum

Light is essential to plant survival and elicits a wide range of plant developmental and physiological responses under different light conditions. A low red-to-far red (R/FR) light ratio induces shade-avoidance responses, including decreased anthocyanin accumulation, whereas a high R/FR light ratio promotes anthocyanin biosynthesis. However, the detailed molecular mechanism underpinning how different R/FR light ratios regulate anthocyanin homeostasis remains elusive, especially in non-model species. Here, we demonstrate that a low R/FR light ratio induced the expression of CmMYB4, which suppressed the anthocyanin activator complex CmMYB6-CmbHLH2, leading to the reduction of anthocyanin accumulation in Chrysanthemum (Chrysanthemum morifolium) petals. Specifically, CmMYB4 recruited the corepressor CmTPL (TOPLESS) to directly bind the CmbHLH2 promoter and suppressed its transcription by impairing histone H3 acetylation. Moreover, the low R/FR light ratio inhibited the PHYTOCHROME INTERACTING FACTOR family transcription factor CmbHLH16, which can competitively bind to CmMYB4 and destabilize the CmMYB4-CmTPL protein complex. Under the high R/FR light ratio, CmbHLH16 was upregulated, which impeded the formation of the CmMYB4-CmTPL complex and released the suppression of CmbHLH2, thus promoting anthocyanin accumulation in Chrysanthemum petals. Our findings reveal a mechanism by which different R/FR light ratios fine-tune anthocyanin homeostasis in flower petals.

The chrysanthemum DEAD-box RNA helicase CmRH56 regulates rhizome outgrowth in response to drought stress

Plants have evolved complex mechanisms to reprogram growth in response to drought stress. In herbaceous perennial plant species, the rhizome, which is normally an organ for propagation and food storage, can also support plant growth in stressful environments, and allows the plant to perennate and survive stress damage. However, the mechanisms that regulate rhizome growth in perennial herbs during abiotic stresses are unknown. Here, we identified a chrysanthemum (Chrysanthemum morifolium) DEAD-box RNA helicase gene, CmRH56, that is specifically expressed in the rhizome shoot apex. Knock down of CmRH56 transcript levels decreased the number of rhizomes and enhanced drought stress tolerance. We determined that CmRH56 represses the expression of a putative gibberellin (GA) catabolic gene, GA2 oxidase6 (CmGA2ox6). Exogenous GA treatment and silencing of CmGA2ox6 resulted in more rhizomes. These results demonstrate that CmRH56 suppresses rhizome outgrowth under drought stress conditions by blocking GA biosynthesis.

Flower color mutation, pink to orange, through CmGATA4 - CCD4a-5 module regulates carotenoids degradation in chrysanthemum

The carotenoids biosynthesis pathway in plants has been studied extensively, yet little is known about the regulatory mechanisms underlying this process, especially for ornamental horticulture plants. In this study, a natural variation of chrysanthemum with orange coloration was identified and compared with the wild type with pink coloration; the content and component of carotenoids were largely enriched in the mutant with orange coloration. CmCCD4a-5, the DNA sequence in both 'Pink yan' and the mutant, was identified and shown to function as a carotenoid degradation enzyme. Compared with 'Pink yan', the mutant shows lower expression level of CmCCD4a-5. Furthermore, CmGATA4 was found to have an opposite expression trend to CmCCD4a-5, and it could directly bind with the CmCCD4a-5 promoter. Taken together, this study demonstrates that CmGATA4 acts as a negative regulator of CmCCD4a-5 and, furthermore, low expression of CmCCD4a-5 resulted in carotenoid accumulation in the mutant.

A natural antisense RNA improves chrysanthemum cold tolerance by regulating the transcription factor DgTCP1

Long noncoding RNAs (lncRNAs) are widely involved in the regulation of plant growth and development, but their mechanism of action in response to cold stress in plants remains unclear. Here, we found an lncRNA transcribed from the antisense strand of DgTCP1 (class I Teosinte branched1/Cycloidea/Proliferating [TCP] transcription factor) of chrysanthemum (Chrysanthemum morifolium Ramat.), named DglncTCP1. During the response of chrysanthemum to cold stress, overexpression of DgTCP1 improved the cold tolerance of chrysanthemum, while the DgTCP1 editing line (dgtcp1) showed decreased tolerance to cold stress. Overexpression of DglncTCP1 also increased the cold tolerance of chrysanthemum, while the DglncTCP1 amiRNA lines (DglncTCP1 amiR-18/38) also showed decreased tolerance to cold stress. Additionally, the overexpression of DglncTCP1 upregulated the expression of DgTCP1. This indicated that DglncTCP1 may play a cis-regulatory role in the regulatory process of DgTCP1 in cold tolerance. DglncTCP1 acts as a scaffold to recruit the histone modification protein DgATX (ARABIDOPSIS TRITHORAX from chrysanthemum) to DgTCP1 to enhance H3K4me3 levels, thereby activating DgTCP1 expression. Moreover, DgTCP1 can directly target DgPOD (peroxidase gene from chrysanthemum) to promote its expression and reduce reactive oxygen species accumulation, thereby improving the cold tolerance of chrysanthemum. In conclusion, these results suggest that natural antisense lncRNA plays a key role in improving the cold tolerance of chrysanthemum.

Increased susceptibility to Chrysanthemum Yellows phytoplasma infection in Atcals7ko plants is accompanied by enhanced expression of carbohydrate transporters

Loss of CALS7 appears to confer increased susceptibility to phytoplasma infection in Arabidopsis, altering expression of genes involved in sugar metabolism and membrane transport. Callose deposition around sieve pores, under control of callose synthase 7 (CALS7), has been interpreted as a mechanical response to limit pathogen spread in phytoplasma-infected plants. Wild-type and Atcals7ko mutants were, therefore, employed to unveil the mode of involvement of CALS7 in the plant's response to phytoplasma infection. The fresh weights of healthy and CY-(Chrysanthemum Yellows) phytoplasma-infected Arabidopsis wild type and mutant plants indicated two superimposed effects of the absence of CALS7: a partial impairment of photo-assimilate transport and a stimulated phytoplasma proliferation as illustrated by a significantly increased phytoplasma titre in Atcal7ko mutants. Further studies solely dealt with the effects of CALS7 absence on phytoplasma growth. Phytoplasma infection affected sieve-element substructure to a larger extent in mutants than in wild-type plants, which was also true for the levels of some free carbohydrates. Moreover, infection induced a similar upregulation of gene expression of enzymes involved in sucrose cleavage (AtSUS5, AtSUS6) and transmembrane transport (AtSWEET11) in mutants and wild-type plants, but an increased gene expression of carbohydrate transmembrane transporters (AtSWEET12, AtSTP13, AtSUC3) in infected mutants only. It remains still unclear how the absence of AtCALS7 leads to gene upregulation and how an increased intercellular mobility of carbohydrates and possibly effectors contributes to a higher susceptibility. It is also unclear if modified sieve-pore structures in mutants allow a better spread of phytoplasmas giving rise to higher titre.

Uneven Levels of 5S and 45S rDNA Site Number and Loci Variations across Wild Chrysanthemum Accessions

Ribosomal DNA (rDNA) is an excellent cytogenetic marker owing to its tandem arrangement and high copy numbers. However, comparative studies have focused more on the number of rDNA site variations within the Chrysanthemum genus, and studies on the types of rDNA sites with the same experimental procedures at the species levels are lacking. To further explore the number and types of rDNA site variations, we combined related data to draw ideograms of the rDNA sites of Chrysanthemum accessions using oligonucleotide fluorescence in situ hybridization (Oligo-FISH). Latent variations (such as polymorphisms of 45S rDNA sites and co-localized 5S-45S rDNA) also occurred among the investigated accessions. Meanwhile, a significant correlation was observed between the number of 5S rDNA sites and chromosome number. Additionally, the clumped and concentrated geographical distribution of different ploidy Chrysanthemum accessions may significantly promote the karyotype evolution. Based on the results above, we identified the formation mechanism of rDNA variations. Furthermore, these findings may provide a reliable method to examine the sites and number of rDNA variations among Chrysanthemum and its related accessions and allow researchers to further understand the evolutionary and phylogenetic relationships of the Chrysanthemum genus.

Chrysanthemum × grandiflora leaf and root transcript profiling in response to salinity stress

As high soil salinity threatens the growth and development of plants, understanding the mechanism of plants' salt tolerance is critical. The Chrysanthemum × grandiflora is a newly developed species with a strong salt resistance that possesses multiple genes controlling its quantitative salt resistance. Because of this multigene control, we chose to investigate the plant stress genes overall responses at the transcriptome level. C. grandiflora were treated with a 200 mM NaCl solution for 12 h to study its effect on the roots and leaves via Illumina RNA sequencing. PAL, CYP73A, and 4CL in the phenylpropanoid biosynthesis pathway were upregulated in roots and leaves. In the salicylic acid signal transduction pathway, TGA7 was upregulated in the roots and leaves, while in the jasmonic acid signal transduction pathway, TIFY9 was upregulated in the roots and leaves. In the ion transporter gene, we identified HKT1 that showed identical expression patterns in the roots and leaves. The impact of NaCl imposition for 12 h was largely due to osmotic effect of salinity on C. grandiflora, and most likely the transcript abundance changes in this study were due to the osmotic effect. In order to verify the accuracy of the Illumina sequencing data, we selected 16 DEGs for transcription polymerase chain reaction (qRT-PCR) analysis. qRT-PCR and transcriptome sequencing analysis revealed that the transcriptome sequencing results were reliable.

An ethylene-responsive transcription factor and a flowering locus KH domain homologue jointly modulate photoperiodic flowering in chrysanthemum

The timely transition from vegetative to reproductive development is coordinated through the quantitative regulation of floral pathway genes in response to physiological and environmental cues. The function of ethylene-responsive element-binding protein (ERF) transcription factors in the regulation of flowering in chrysanthemum (Chrysanthemum morifolium Ramat.) is not well understood. Here, chrysanthemum overexpressing CmERF110 flowered earlier than the wild-type plants, while those in which CmERF110 was suppressed flowered later. RNA-seq results revealed that several genes involved in the circadian rhythm were transcribed differently in CmERF110 transgenic plants from that of the wild-type plants. The rhythm peak of the circadian clock genes in transgenic plants was delayed. Yeast two-hybrid screening of CmERF110 interactors identified a chrysanthemum FLOWERING LOCUS KH DOMAIN (FLK) homologue CmFLK, which was further confirmed with both in vitro and in vivo assays. KEGG pathway enrichment also revealed that CmFLK is involved in the regulation of circadian rhythm-related genes. CmFLK transgenic plants showed a change in flowering time and delayed rhythm peak of the circadian rhythm genes. Taken together, the present data not only suggest that CmERF110 interacts with CmFLK to promote floral transition by tuning the circadian clock, but also provides evidence for the evolutionary conservation of the components in the autonomous pathway in chrysanthemum.

Regulation of lignin biosynthesis by an atypical bHLH protein CmHLB in Chrysanthemum

Stem mechanical strength is one of the most important agronomic traits that affects the resistance of plants against insects and lodging, and plays an essential role in the quality and yield of plants. Several transcription factors regulate mechanical strength in crops. However, mechanisms of stem strength formation and regulation remain largely unexplored, especially in ornamental plants. In this study, we identified an atypical bHLH transcription factor CmHLB (HLH PROTEIN INVOLVED IN LIGNIN BIOSYNTHESIS) in chrysanthemum, belonging to a small bHLH sub-family - the PACLOBUTRAZOL RESISTANCE (PRE) family. Overexpression of CmHLB in chrysanthemum significantly increased mechanical strength of the stem, cell wall thickness, and lignin content, compared with the wild type. In contrast, CmHLB RNA interference lines exhibited the opposite phenotypes. RNA-seq analysis indicated that CmHLB promoted the expression of genes involved in lignin biosynthesis. Furthermore, we demonstrated that CmHLB interacted with Chrysanthemum KNOTTED ARABIDOPSIS THALIANA7 (CmKNAT7) through the KNOX2 domain, which has a conserved function, i.e. it negatively regulates secondary cell wall formation of fibres and lignin biosynthesis. Collectively, our results reveal a novel role for CmHLB in regulating lignin biosynthesis by interacting with CmKNAT7 and affecting stem mechanical strength in Chrysanthemum.

Lysine decrotonylation of glutathione peroxidase at lysine 220 site increases glutathione peroxidase activity to resist cold stress in chrysanthemum

Lysine crotonylation is a protein post-translational modification that has been newly discovered in recent years. There are few studies on the lysine crotonylation of proteins in plants, and their functions in response to cold stress are still unclear. In this study, the chrysanthemum (Chrysanthemum morifolium Ramat.) glutathione peroxidase (GPX) gene was selected and named DgGPX1, and was found to be responsive to low temperature. Overexpression of DgGPX1 improved the cold resistance of transgenic chrysanthemum by increasing GPX activity to reduce the accumulation of reactive oxygen species (ROS) under low-temperature conditions. Furthermore, the level of DgGPX1 lysine crotonylation at lysine (K) 220 decreased under low temperature in chrysanthemum. Lysine decrotonylation of DgGPX1 at K220 further increased GPX activity to reduce ROS accumulation under cold stress, and thereby enhanced the cold resistance of chrysanthemum. The above results show that lysine decrotonylation of DgGPX1 at K220 increases GPX activity to resist cold stress in chrysanthemum.