Cloning and functional verification of the CmHSP17.9 gene from chrysanthemum.

Small molecular heat shock proteins (sHSPs) belong to the HSP family of molecular chaperones. Under high-temperature stress, they can prevent the aggregation of irreversible proteins and maintain the folding of denatured proteins to enhance heat resistance. In this study, the CmHSP17.9-1 and CmHSP17.9-2 genes, which were cloned from chrysanthemum (Chrysanthemum×morifolium 'Jinba') by homologous cloning, had a complete open reading frame of 480 bp each, encoding 159 amino acids. The protein subcellular localization analysis showed that CmHSP17.9-1 and CmHSP17.9-2 were located in the cytoplasm and mostly aggregated in granules, especially around the nucleus. Real-time quantitative PCR (qRT-PCR) analysis showed that the relative expression level of the CmHSP17.9-1 and CmHSP17.9-2 genes was highest in the terminal buds of the chrysanthemum, followed by the leaves. CmHSP17.9-1 and CmHSP17.9-2 overex-pression vectors were constructed and used to transform the chrysanthemum; overexpression of these genes led to the chrysanthemum phenotypes being less affected by high-temperature, and the antioxidant capacity was enhanced. The results showed that chrysanthemum with overex-pression of the CmHSP17.9-1 and CmHSP17.9-2 genes had stronger tolerance than the wild type chrysanthemum after high-temperature treatment or some degree of heat exercise, and overex-pression of the CmHSP17.9-1 gene led to stronger heat resistance than that of the CmHSP17.9-2 gene, providing an important theoretical basis for the subsequent molecular breeding and pro-duction applications of chrysanthemum.

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.

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.

The A subunit of vacuolar H+-ATPase gene (CmVHA-A) plays opposite roles in plant growth and drought tolerance of chrysanthemum under different growing conditions.

As the most prominent proton pumps in plants, vacuolar H+-ATPases (VHAs) comprise multiple subunits that are important for physiological processes and stress tolerance in plants. However, few studies on the roles of subunit genes of VHAs in chrysanthemum have been reported to date. In this study, the gene of A subunit of V-ATPase in chrysanthemum (CmVHA-A) was cloned and identified. CmVHA-A was conserved with VHA-A proteins from other plants. Expression analysis showed that CmVHA-A was highly expressed in most tissues of chrysanthemum except for the flower bud, and was readily induced by polyethylene glycol (PEG) treatment. Functional analysis demonstrated that CmVHA-A exerted a negative influence on the growth and development of shoot and root of chrysanthemum under normal conditions. RNA-sequencing (RNA-seq) analysis revealed the possible explanations for phenotypic differences between transgenic and wild-type (WT) plants. Under drought conditions, CmVHA-A positively affected the drought tolerance of chrysanthemum by enhancing antioxidase activity and alleviating photosynthetic disruption. Overall, CmVHA-A plays opposite roles in plant growth and drought tolerance of chrysanthemums under different growing conditions.

[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.

A GASA Protein Family Gene, CmGEG, Inhibits Petal Growth in Chrysanthemum.

The diversity in the petal morphology of chrysanthemums makes this species an excellent model for investigating the regulation mechanisms of petal size. However, our understanding of the molecular regulation of petal growth in chrysanthemums remains limited. The GASA (gibberellic acid [GA]-stimulated Arabidopsis) protein plays a significant role in various aspects of plant growth and development. Previous studies have indicated that GEG (a gerbera homolog of the gibberellin-stimulated transcript 1 [GAST1] from tomato) is involved in regulating ray petal growth by inhibiting cell expansion in gerberas. In this study, we successfully cloned the GASA family gene from chrysanthemums, naming it CmGEG, which shares 81.4% homology with GEG. Our spatiotemporal expression analysis revealed that CmGEG is expressed in all tissues, with the highest expression levels observed in the ray florets, particularly during the later stages of development. Through transformation experiments, we demonstrated that CmGEG inhibits petal elongation in chrysanthemums. Further observations indicated that CmGEG restricts cell elongation in the top, middle, and basal regions of the petals. To investigate the relationship between CmGEG and GA in petal growth, we conducted a hormone treatment assay using detached chrysanthemum petals. Our results showed that GA promotes petal elongation while downregulating CmGEG expression. In conclusion, the constrained growth of chrysanthemum petals may be attributed to the inhibition of cell elongation by CmGEG, a process regulated by GA.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.