The Chrysanthemum nankingense Genome Provides Insights into the Evolution and Diversification of Chrysanthemum Flowers and Medicinal Traits

Chrysanthemum morifolium as a traditional herb: A review of historical development, classification, phytochemistry, pharmacology and application

ETHNOPHARMACOLOGICAL RELEVANCE

In recent years, Chinese herbal medicine has gained more and more recognition in disease prevention and control due to its low toxicity and comprehensive treatment. C. morifolium (Chrysanthemum morifolium Ramat.), as the medicine food homology plant with the bioactivity of anti-oxidation, anti-inflammatory, neuroprotection and cardiovascular protection, has important therapeutic effects and health benefits for colds, inflammation, cardiovascular diseases and various chronic diseases.

AIM OF THE STUDY

By reviewing the historical development, classification and distribution of germplasm resources, phytochemistry, pharmacology, and modern application of C. morifolium, the paper provides a reliable basis for the further research and application of chrysanthemum as therapeutic agents and functional additives.

MATERIALS AND METHODS

The literature and information about C. morifolium published in the last ten years were collected from various platforms, including Google Scholar, PubMed, ScienceDirect, Web of Science and China Knowledge Network.

RESULTS

A comprehensive analysis confirmed that C. morifolium originated in China, and it went through the development process from food and tea to medicine for more than 3000 years. During this period, different cultivars emerged through several breeding techniques and were distributed throughout the world. Moreover, A variety of chemical components such as flavonoids, phenolic acids, volatile oils, and terpenes in chrysanthemum have been proven they possess various pharmacology of anti-inflammatory, anti-oxidant, and prevention of chronic diseases by regulating inflammatory cytokines, oxidative stress responses and signaling pathways, which are the essential conditions to play a role in TCM, nutraceuticals and diet.

CONCLUSION

This paper provides a comprehensive review of historical development, classification, phytochemistry, pharmacology, and modern application of C. morifolium. However, future studies should continue to focus on the bioactive compounds and the synergistic mechanism of the "multi-component, multi-target, and multi-pathway" of chrysanthemum, and it is necessary to develop more innovative products with therapeutic effects.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family

BACKGROUND

Terpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated.

RESULTS

We identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality.

CONCLUSIONS

We systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Development and evolution of the Asteraceae capitulum

Asteraceae represent one of the largest and most diverse families of plants. The evolutionary success of this family has largely been contributed to their unique inflorescences, capitula that mimic solitary flowers but are typically aggregates of multiple florets. Here, we summarize the recent molecular and genetic level studies that have promoted our understanding of the development and evolution of capitula. We focus on new results on patterning of the enlarged meristem resulting in the iconic phyllotactic arrangement of florets in Fibonacci numbers of spirals. We also summarize the current understanding of the genetic networks regulating the characteristic reproductive traits in the family such as floral dimorphism and differentiation of highly specialized floral organs. So far, developmental studies in Asteraceae are still limited to a very narrow selection of model species. Along with the recent advancements in genomics and phylogenomics, Asteraceae and its relatives provide an outstanding model clade for extended evo-devo studies to exploit the morphological diversity and the underlying molecular networks and to translate this knowledge to the breeding of the key crops in the family.

Comparative metabolomic analysis reveals nutritional properties and pigmentation mechanism of tea-scented rosehips

BACKGROUND

The fruits of the genus Rosa, commonly known as rosehips, have attracted significant attention owing to their rich content of various bioactive compounds. However, their utility is generally secondary to the ornamental appeal of their flowers. This study aimed to explore the quality differences among tea-scented rosehips found in Yunnan, China, including those of Rosa odorata var. odorata (RO), Rosa odorata var. gigantea (RG), and Rosa yangii (RY). Morphological characteristics, chemical composition, and antioxidant activity of their fruits were evaluated.

RESULTS

The study revealed significant variability in composition and biological activities based on fruit color. RO exhibited the highest levels of polyphenols, flavonoids, anthocyanins, carotenoids, and vitamin C, with the strongest antioxidant activity (10.99 μmol Trolox·g-1 ), followed by RG (7.91 μmol Trolox·g-1 ) and RY (6.52 μmol Trolox·g-1 ). This supports RO's potential as a functional food source. Untargeted metabolomics identified and quantified 502 metabolites, with flavonoids (171) and phenolic acids (147) as the main metabolites. The differential metabolites among the fruits are primarily enriched for flavonoid biosynthesis and phenylpropanoid biosynthesis pathways. Insights into color formation supported the role of anthocyanins, flavones, and flavonols in fruit color variation.

CONCLUSION

Tea-scented rosehips offer vibrant colors and high nutritional value with potent biological activities. Rosa odorata var. odorata stands out as a functional food source owing to its rich bioactive compounds. These findings lay the groundwork for utilizing rosehips in functional foods, health supplements, and food additives, emphasizing the practical and beneficial applications of Rosa spp. independent of their ornamental value. © 2023 Society of Chemical Industry.

cla-miR164-NO APICAL MERISTEM (ClNAM) regulates the inflorescence architecture development of Chrysanthemum lavandulifolium

Chrysanthemum × morifolium has great ornamental and economic value on account of its exquisite capitulum. However, previous studies have mainly focused on the corolla morphology of the capitulum. Such an approach cannot explain the variable inflorescence architecture of the chrysanthemum. Previous research from our group has shown that NO APICAL MERISTEM (ClNAM) is likely to function as a hub gene in capitulum architecture in the early development stage. In the present study, ClNAM was used to investigate the function of these boundary genes in the capitulum architecture of Chrysanthemum lavandulifolium, a closely related species of C. × morifolium in the genus. Modification of ClNAM in C. lavandulifolium resulted in an advanced initiation of the floral primordium at the capitulum. As a result, the receptacle morphology was altered and the number of florets decreased. The ray floret corolla was shortened, but the disc floret was elongated. The number of capitula increased significantly, arranged in more densely compounded corymbose synflorescences. The yeast and luciferase reporter system revealed that ClAP1, ClRCD2, and ClLBD18 target and activate ClNAM. Subsequently, ClNAM targets and activates ClCUC2a/c, which regulates the initiation of floral and inflorescence in C. lavandulifolium. ClNAM was also targeted and cleaved by cla-miR164 in this process. In conclusion, this study established a boundary gene regulatory network with cla-miR164-ClNAM as the hub. This network not only influences the architecture of capitulum, but also affects compound corymbose synflorescences of the C. lavandulifolium. These results provide new insights into the mechanisms regulating inflorescence architecture in chrysanthemum.

The genomes of Dahlia pinnata, Cosmos bipinnatus, and Bidens alba in tribe Coreopsideae provide insights into polyploid evolution and inulin biosynthesis

BACKGROUND

The Coreopsideae tribe, a subset of the Asteraceae family, encompasses economically vital genera like Dahlia, Cosmos, and Bidens, which are widely employed in medicine, horticulture, ecology, and food applications. Nevertheless, the lack of reference genomes hinders evolutionary and biological investigations in this tribe.

RESULTS

Here, we present 3 haplotype-resolved chromosome-level reference genomes of the tribe Coreopsideae, including 2 popular flowering plants (Dahlia pinnata and Cosmos bipinnatus) and 1 invasive weed plant (Bidens alba), with assembled genome sizes 3.93 G, 1.02 G, and 1.87 G, respectively. We found that Gypsy transposable elements contribute mostly to the larger genome size of D. pinnata, and multiple chromosome rearrangements have occurred in tribe Coreopsideae. Besides the shared whole-genome duplication (WGD-2) in the Heliantheae alliance, our analyses showed that D. pinnata and B. alba each underwent an independent recent WGD-3 event: in D. pinnata, it is more likely to be a self-WGD, while in B. alba, it is from the hybridization of 2 ancestor species. Further, we identified key genes in the inulin metabolic pathway and found that the pseudogenization of 1-FEH1 and 1-FEH2 genes in D. pinnata and the deletion of 3 key residues of 1-FFT proteins in C. bipinnatus and B. alba may probably explain why D. pinnata produces much more inulin than the other 2 plants.

CONCLUSIONS

Collectively, the genomic resources for the Coreopsideae tribe will promote phylogenomics in Asteraceae plants, facilitate ornamental molecular breeding improvements and inulin production, and help prevent invasive weeds.

Identification and characterization of CYC2-like genes related to floral symmetric development in Tagetes erecta (Asteraceae)

Marigold (Tagetes erecta) is an annual herbaceous flower belonging to Asteraceae, whose capitulum is composed of bilateral symmetry ray florets on the outer periphery and radial symmetry disk florets on the inside. The flower symmetry evolution from radial symmetry to bilateral symmetry has changed the morphology, inflorescence architecture and function of florets among several lineages in Asteraceae. Several studies have identified that CYC2 genes in TCP transcription factor family are the key genes regulating the flower morphogenesis, such as corolla symmetry and stamen development. Here, seven TeCYC2 genes were cloned and phylogenetically grouped into the CYC2 branch of TCP transcription family. TeCYC2c and TeCYC2d were found to be expressed specifically in ray florets, TeCYC2b was strongly expressed in both ray and disk florets, TeCYC2g was significantly higher expressed in ray florets than in disk florets, while TeCYC2a, TeCYC2e1 and TeCYC2e2 were significantly expressed in disk florets, according to an examination of the expression profile. Among the ectopic expression lines of seven TeCYC2 genes in Arabidopsis thaliana, the flower symmetry of all transgenic lines was changed from radial symmetry to bilateral symmetry, and only the reproductive growth of TeCYC2c lines was affected. In TeCYC2c transgenic Arabidopsis, the pollen sac was difficult to crack, and the filaments were shorter than the pistils, resulting in a significant decrease in the seed setting rate. All TeCYC2 proteins were localized in the nucleus. Eight pairs of interactions between TeCYC2 proteins were validated by Y2H and BiFC assays, indicating the possibility of TeCYC2 proteins forming homodimers or heterodimers to improve functional specificity. Our findings verified the main regulatory role of TeCYC2c on the development of corollas and stamen in marigold, and analyzed the interaction network of the formation mechanism of floral symmetry in two florets, which provided more insights into the expansion of CYC2 genes in the evolution of Asteraceae inflorescence and contributed to elucidate the complex regulatory network, as well as the molecular breeding concerning flower form diversity in marigold.

Designing of future ornamental crops: a biotechnological driven perspective

With a basis in human appreciation of beauty and aesthetic values, the new era of ornamental crops is based on implementing innovative technologies and transforming symbols into tangible assets. Recent advances in plant biotechnology have attracted considerable scientific and industrial interest, particularly in terms of modifying desired plant traits and developing future ornamental crops. By utilizing omics approaches, genomic data, genetic engineering, and gene editing tools, scientists have successively explored the underlying molecular mechanism and potential gene(s) behind trait regulation such as floral induction, plant architecture, stress resistance, plasticity, adaptation, and phytoremediation in ornamental crop species. These signs of progress lay a theoretical and practical foundation for designing and enhancing the efficiency of ornamental plants for a wide range of applications. In this review, we briefly summarized the existing literature and advances in biotechnological approaches for the improvement of vital traits in ornamental plants. The future ornamental plants, such as light-emitting plants, biotic/abiotic stress detectors, and pollution abatement, and the introduction of new ornamental varieties via domestication of wild species are also discussed.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Chromosome-level reference genome assembly provides insights into the evolution of Pennisetum alopecuroides

Pennisetum alopecuroides is an important forage grass resource, which plays a vital role in ecological environment improvement. Therefore, the acquisition of P. alopecuroides genome resources is conducive to the study of the adaptability of Pennisetum species in ecological remediation and forage breeding development. Here we assembled a P. alopecuroides cv. 'Liqiu' genome at the chromosome level with a size of approximately 845.71 Mb, contig N50 of 84.83Mb, and genome integrity of 99.13% as assessed by CEGMA. A total of 833.41-Mb sequences were mounted on nine chromosomes by Hi-C technology. In total, 60.66% of the repetitive sequences and 34,312 genes were predicted. The genomic evolution analysis showed that P. alopecuroides cv. 'Liqiu' was isolated from Setaria 7.53-13.80 million years ago and from Cenchrus 5.33-8.99 million years ago, respectively. The whole-genome event analysis showed that P. alopecuroides cv. 'Liqiu' underwent two whole-genome duplication (WGD) events in the evolution process, and the duplication events occurred at a similar time to that of Oryza sativa and Setaria viridis. The completion of the genome sequencing of P. alopecuroides cv. 'Liqiu' provides data support for mining high-quality genetic resources of P. alopecuroides and provides a theoretical basis for the origin and evolutionary characteristics of Pennisetum.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Genome-wide identification of the MIKCc-type MADS-box gene family in Chrysanthemum lavandulifolium reveals their roles in the capitulum development

Chrysanthemum ×morifolium is well known throughout the world for its diverse and exquisite flower types. However, due to the complicated genetic background of C. ×morifolium, it is difficult to understand the molecular mechanism of its flower development. And it limits the molecular breeding of improving chrysanthemum flower types. C. ×morifolium has the typical radial capitulum, and many researches showed that the members of the MIKCc-type MADS box gene family play a key role in the formation and development of the capitulum. However, it has been difficult to isolate the important MIKCc and investigate their roles in this process due to the lack of genomic information in chrysanthemum. Here, we identified MIKCc-type MADS box genes at whole genome-wide level in C. lavandulifolium, a diploid species closely related to C. ×morifolium, and investigated their roles in capitulum development by gene expression pattern analysis and protein interaction analysis. A total of 40 ClMIKCc were identified and were phylogenetically grouped into 12 clades. Members of all clades showed different enriched expression patterns during capitulum formation. We speculate that the E-class genes in C. lavandulifolium underwent subfunctionalization because they have a significantly expanded, more diverse expression patterns, and specifically tissue expression than AtSEPs. Meanwhile, we detected the C-class expressed in disc floret corolla, which could be the clue to explore the morphological differences between disc and ray floret corolla. In addition, the potential roles of some MIKCcs in complex inflorescence formation were explored by comparing the number and phylogenetic relationship of MIKCc subfamily members in Asteraceae with different capitulum types. Members of the FLC branch in Asteraceae were found to be possibly related to the differentiation and development of the ray floret.

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.

IGTCM: An integrative genome database of traditional Chinese medicine plants

Fully understanding traditional Chinese medicines (TCMs) is still challenging because of the extreme complexity of their chemical components and mechanisms of action. The TCM Plant Genome Project aimed to obtain genetic information, determine gene functions, discover regulatory networks of herbal species, and elucidate the molecular mechanisms involved in the disease prevention and treatment, thereby accelerating the modernization of TCMs. A comprehensive database that contains TCM-related information will provide a vital resource. Here, we present an integrative genome database of TCM plants (IGTCM) that contains 14,711,220 records of 83 annotated TCM-related herb genomes, including 3,610,350 genes, 3,534,314 proteins and corresponding coding sequences, and 4,032,242 RNAs, as well as 1033 non-redundant component records for 68 herbs, downloaded and integrated from the GenBank and RefSeq databases. For minimal interconnectivity, each gene, protein, and component was annotated using the eggNOG-mapper tool and Kyoto Encyclopedia of Genes and Genomes database to acquire pathway information and enzyme classifications. These features can be linked across several species and different components. The IGTCM database also provides visualization and sequence similarity search tools for data analyses. These annotated herb genome sequences in IGTCM database are a necessary resource for systematically exploring genes related to the biosynthesis of compounds that have significant medicinal activities and excellent agronomic traits that can be used to improve TCM-related varieties through molecular breeding. It also provides valuable data and tools for future research on drug discovery and the protection and rational use of TCM plant resources. The IGTCM database is freely available at http://yeyn.group:96/.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

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.

Intraspecific variation in genome size in Artemisia argyi determined using flow cytometry and a genome survey

Different collections and accessions of Artemisia argyi (Chinese mugwort) harbour considerable diversity in morphology and bioactive compounds, but no mechanisms have been reported that explain these variations. We studied genome size in A. argyi accessions from different regions of China by flow cytometry. Genome size was significantly distinct among origins of these 42 Chinese mugwort accessions, ranging from 8.428 to 11.717 pg. There were no significant intraspecific differences among the 42 accessions from the five regions of China. The clustering analysis showed that these 42 A. argyi accessions could be divided into three groups, which had no significant relationship with geographical location. In a genome survey, the total genome size of A. argyi (A15) was estimated to be 7.852 Gb (or 8.029 pg) by K-mer analysis. This indicated that the results from the two independent methods are consistent, and that the genome survey can be used as an adjunct to flow cytometry to compensate for its deficiencies. In addition, genome survey can provide the information about heterozygosity, repeat sequences, GC content and ploidy of A. argyi genome. The nuclear DNA contents determined here provide a new reference for intraspecific variation in genome size in A. argyi, and may also be a potential resource for the study of genetic diversity and for breeding new cultivar.

Biosynthesis of α-Bisabolol by Farnesyl Diphosphate Synthase and α-Bisabolol Synthase and Their Related Transcription Factors in Matricaria recutita L

The essential oil of German chamomile (Matricaria recutita L.) is widely used in food, cosmetics, and the pharmaceutical industry. α-Bisabolol is the main active substance in German chamomile. Farnesyl diphosphate synthase (FPS) and α-bisabolol synthase (BBS) are key enzymes related to the α-bisabolol biosynthesis pathway. However, little is known about the α-bisabolol biosynthesis pathway in German chamomile, especially the transcription factors (TFs) related to the regulation of α-bisabolol synthesis. In this study, we identified MrFPS and MrBBS and investigated their functions by prokaryotic expression and expression in hairy root cells of German chamomile. The results suggest that MrFPS is the key enzyme in the production of sesquiterpenoids, and MrBBS catalyzes the reaction that produces α-bisabolol. Subcellular localization analysis showed that both MrFPS and MrBBS proteins were located in the cytosol. The expression levels of both MrFPS and MrBBS were highest in the extension period of ray florets. Furthermore, we cloned and analyzed the promoters of MrFPS and MrBBS. A large number of cis-acting elements related to light responsiveness, hormone response elements, and cis-regulatory elements that serve as putative binding sites for specific TFs in response to various biotic and abiotic stresses were identified. We identified and studied TFs related to MrFPS and MrBBS, including WRKY, AP2, and MYB. Our findings reveal the biosynthesis and regulation of α-bisabolol in German chamomile and provide novel insights for the production of α-bisabolol using synthetic biology methods.

Creating novel ornamentals via new strategies in the era of genome editing

Ornamental breeding has traditionally focused on improving novelty, yield, quality, and resistance to biotic or abiotic stress. However, achieving these goals has often required laborious crossbreeding, while precise breeding techniques have been underutilized. Fortunately, recent advancements in plant genome sequencing and editing technology have opened up exciting new frontiers for revolutionizing ornamental breeding. In this review, we provide an overview of the current state of ornamental transgenic breeding and propose four promising breeding strategies that have already proven successful in crop breeding and could be adapted for ornamental breeding with the help of genome editing. These strategies include recombination manipulation, haploid inducer creation, clonal seed production, and reverse breeding. We also discuss in detail the research progress, application status, and feasibility of each of these tactics.

PlantNLRatlas: a comprehensive dataset of full- and partial-length NLR resistance genes across 100 chromosome-level plant genomes

Plants have evolved two layers of protection against biotic stress: PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). The primary mechanism of ETI involves nucleotide-binding leucine-rich repeat immune receptors (NLRs). Although NLR genes have been studied in several plant species, a comprehensive database of NLRs across a diverse array of species is still lacking. Here, we present a thorough analysis of NLR genes across 100 high-quality plant genomes (PlantNLRatlas). The PlantNLRatlas includes a total of 68,452 NLRs, of which 3,689 are full-length and 64,763 are partial-length NLRs. The majority of NLR groups were phyletically clustered. In addition, the domain sequences were found to be highly conserved within each NLR group. Our PlantNLRatlas dataset is complementary to RefPlantNLR, a collection of NLR genes which have been experimentally confirmed. The PlantNLRatlas should prove helpful for comparative investigations of NLRs across a range of plant groups, including understudied taxa. Finally, the PlantNLRatlas resource is intended to help the field move past a monolithic understanding of NLR structure and function.

Chromosome-level genome of a leaf vegetable Glebionis coronaria provides insights into the biosynthesis of monoterpenoids contributing to its special aroma

Glebionis coronaria is a popular vegetable with special aroma and a medical plant in East Asia and Mediterranean, but its biological studies and breeding have been hindered by the lack of reference genome. Here, we present a chromosome-level reference genome of G. coronaria, with assembled genome size of 6.8 Gb, which is the largest among all the published genomes of diploid Asteraceae species. The large genome size of G. coronaria is mainly caused by the recent widespread explosions of long-terminal-repeat retrotransposons. Analyses of macro-synteny and synonymous mutation rate distribution indicate that the G. coronaria genome experienced a whole-genome triplication at 40-45 million years ago, shared with all Asteraceae species. In subtribe Artemisiinae, Glebionis arose before the divergence of Chrysanthemum from Artemisia, and Glebionis species evolved much faster than Chrysanthemum and Artemisia species. In G. coronaria, the synthesis genes of monoterpenoids 8-oxocitronellyl enol and isopiperitenone were expanded, and the higher expressions of these expanded genes in leaves and stems may contribute to its special aroma. The G. coronaria genomic resources will promote the evolution studies of Asteraceae, the metabolism mechanism studies of bioactive compounds, and the breeding improvement of agronomic traits in G. coronaria.

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.

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.

The genomes of chicory, endive, great burdock and yacon provide insights into Asteraceae palaeo-polyploidization history and plant inulin production

Inulin is an important reserve polysaccharide in Asteraceae plants, and is also widely used as a sweetener, a source of dietary fibre and prebiotic. Nevertheless, a lack of genomic resources for inulin-producing plants has hindered extensive studies on inulin metabolism and regulation. Here, we present chromosome-level reference genomes for four inulin-producing plants: chicory (Cichorium intybus), endive (Cichorium endivia), great burdock (Arctium lappa) and yacon (Smallanthus sonchifolius), with assembled genome sizes of 1.28, 0.89, 1.73 and 2.72 Gb, respectively. We found that the chicory, endive and great burdock genomes were shaped by whole genome triplication (WGT-1), and the yacon genome was shaped by WGT-1 and two subsequent whole genome duplications (WGD-2 and WGD-3). A yacon unique whole genome duplication (WGD-3) occurred 5.6-5.8 million years ago. Our results also showed the genome size difference between chicory and endive is largely due to LTR retrotransposons, and rejected a previous hypothesis that chicory is an ancestor of endive. Furthermore, we identified fructan-active-enzyme and transcription-factor genes, and found there is one copy in chicory, endive and great burdock but two copies in yacon for most of these genes, except for the 1-FEH II gene which is significantly expanded in chicory. Interestingly, inulin synthesis genes 1-SST and 1-FFT are located close to each other, as are the degradation genes 1-FEH I and 1-FEH II. Finally, we predicted protein structures for 1-FFT genes to explore the mechanism determining inulin chain length.

Genome-wide identification and expression analysis of SBP-box gene family reveal their involvement in hormone response and abiotic stresses in Chrysanthemum nankingense

SQUAMOSA promoter-binding-protein (SBP)-box family proteins are a class of plant-specific transcription factors, and widely regulate the development of floral and leaf morphology in plant growth and involve in environment and hormone signal response. In this study, we isolated and identified 21 non-redundant SBP-box genes in Chrysanthemum nankingense with bioinformatics analysis. Sequence alignments of 21 CnSBP proteins discovered a highly conserved SBP domain including two zinc finger-like structures and a nuclear localization signal region. According to the amino acid sequence alignments, 67 SBP-box genes from Arabidopsis thaliana, rice, Artemisia annua and C. nankingense were clustered into eight groups, and the motif and gene structure analysis also sustained this classification. The gene evolution analysis indicated the CnSBP genes experienced a duplication event about 10 million years ago (Mya), and the CnSBP and AtSPL genes occurred a divergence at 24 Mya. Transcriptome data provided valuable information for tissue-specific expression profiles of the CnSBPs, which highly expressed in floral tissues and differentially expressed in leaf, root and stem organs. Quantitative Real-time Polymerase Chain Reaction data showed expression patterns of the CnSBPs under exogenous hormone and abiotic stress treatments, separately abscisic acid, salicylic acid, gibberellin A3, methyl jasmonate and ethylene spraying as well as salt and drought stresses, indicating that the candidate CnSBP genes showed differentiated spatiotemporal expression patterns in response to hormone and abiotic stresses. Our study provides a systematic genome-wide analysis of the SBP-box gene family in C. nankingense. In general, it provides a fundamental theoretical basis that SBP-box genes may regulate the resistance of stress physiology in chrysanthemum via exogenous hormone pathways.

Phylogeographic structure of Heteroplexis (Asteraceae), an endangered endemic genus in the limestone karst regions of southern China

Though the karst regions in south and southwest China are plant diversity hotspots, our understanding of the phylogeography and evolutionary history of the plants there remains limited. The genus Heteroplexis (Asteraceae) is one of the typical representative plants isolated by karst habitat islands, and is also an endangered and endemic plant to China. In this study, species-level phylogeographic analysis of the genus Heteroplexis was conducted using restriction site-associated DNA sequencing (RADseq). The genetic structure showed a clear phylogeographic structure consistent with the current species boundaries in the H. microcephala, H. incana, H. vernonioides, H. sericophylla, and H. impressinervia. The significant global (R = 0.37, P < 0.01) and regional (R = 0.650.95, P < 0.05) isolation by distance (IBD) signals among species indicate strong geographic isolation in the karst mountains, which may result in chronically restricted gene flow and increased genetic drift and differentiation. Furthermore, the phylogeographic structure of Heteroplexis suggested a southward migration since the last glacial period. Demographic analysis revealed the karst mountains as a refuge for Heteroplexis species. Finally, both Treemix and ABBA-BABA statistic detected significant historical gene flow between species. Significant historical gene flow and long-term stability of effective population size (Ne) together explain the high genome-wide genetic diversity among species (π = 0.05370.0838). However, the recent collapse of Ne, widespread inbreeding within populations, and restricted contemporary gene flow suggest that Heteroplexis species are probably facing a high risk of genetic diversity loss. Our results help to understand the evolutionary history of karst plants and guide conservation.

Combination of long-read and short-read sequencing provides comprehensive transcriptome and new insight for Chrysanthemum morifolium ray-floret colorization

Chrysanthemum morifolium is one of the most popular ornamental plants globally. Owing to its large and complex genome (around 10 Gb, segmental hexaploid), it has been difficult to obtain comprehensive transcriptome, which will promote to perform new breeding technique, such as genome editing, in C. morifolium. In this study, we used single-molecule real-time (SMRT) sequencing and RNA-seq technologies, combined them with an error-correcting process, and obtained high-coverage ray-floret transcriptome. The SMRT-seq data increased the ratio of long mRNAs containing complete open-reading frames, and the combined dataset provided a more complete transcriptomic data than those produced from either SMRT-seq or RNA-seq-derived transcripts. We finally obtained 'Sei Arabella' transcripts containing 928,645 non-redundant mRNA, which showed 96.6% Benchmarking Universal Single-Copy Orthologs (BUSCO) score. We also validated the reliability of the dataset by analyzing a mapping rate, annotation and transcript expression. Using the dataset, we searched anthocyanin biosynthesis gene orthologs and performed a qRT-PCR experiment to assess the usability of the dataset. The assessment of the dataset and the following analysis indicated that our dataset is reliable and useful for molecular biology. The combination of sequencing methods provided genetic information and a way to analyze the complicated C. morifolium transcriptome.

The genome of single-petal jasmine (Jasminum sambac) provides insights into heat stress tolerance and aroma compound biosynthesis

Jasmine [Jasminum sambac (L.) Aiton] is a commercially important cultivated plant species known for its fragrant flowers used in the perfume industry, medicine and cosmetics. In the present study, we obtained a draft genome for the J. sambac cultivar 'Danbanmoli' (JSDB, a single-petal phenotype). We showed that the final genome of J. sambac was 520.80 Mb in size (contig N50 = 145.43 kb; scaffold N50 = 145.53 kb) and comprised 35,363 genes. Our analyses revealed that the J. sambac genome has undergone only an ancient whole-genome duplication (WGD) event. We estimated that the lineage that has given rise to J. sambac diverged from the lineage leading to Osmanthus fragrans and Olea europaea approximately 31.1 million years ago (Mya). On the basis of a combination of genomic and transcriptomic analyses, we identified 92 transcription factors (TFs) and 206 genes related to heat stress response. Base on a combination of genomic, transcriptomic and metabolomic analyses, a range of aroma compounds and genes involved in the benzenoid/phenylpropanoid and terpenoid biosynthesis pathways were identified. In the newly assembled J. sambac genome, we identified a total of 122 MYB, 122 bHLH and 69 WRKY genes. Our assembled J. sambac JSDB genome provides fundamental knowledge to study the molecular mechanism of heat stress tolerance, and improve jasmine flowers and dissect its fragrance.

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.

Genome sequencing reveals chromosome fusion and extensive expansion of genes related to secondary metabolism in Artemisia argyi

Artemisia argyi, as famous as Artemisia annua, is a medicinal plant with huge economic value in the genus of Artemisia and has been widely used in the world for about 3000 years. However, a lack of the reference genome severely hinders the understanding of genetic basis for the active ingredient synthesis of A. argyi. Here, we firstly report a complex chromosome-level genome assembly of A. argyi with a large size of 8.03 Gb, with features of high heterozygosity (2.36%), high repetitive sequences (73.59%) and a huge number of protein-coding genes (279 294 in total). The assembly reveals at least three rounds of whole-genome duplication (WGD) events, including a recent WGD event in the A. argyi genome, and a recent burst of transposable element, which may contribute to its large genome size. The genomic data and karyotype analyses confirmed that A. argyi is an allotetraploid with 34 chromosomes. Intragenome synteny analysis revealed that chromosomes fusion event occurred in the A. argyi genome, which elucidates the changes in basic chromosome numbers in Artemisia genus. Significant expansion of genes related to photosynthesis, DNA replication, stress responses and secondary metabolism were identified in A. argyi, explaining the extensive environmental adaptability and rapid growth characteristics. In addition, we analysed genes involved in the biosynthesis pathways of flavonoids and terpenoids, and found that extensive gene amplification and tandem duplication contributed to the high contents of metabolites in A. argyi. Overall, the reference genome assembly provides scientific support for evolutionary biology, functional genomics and breeding in A. argyi and other Artemisia species.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Characterization and expression analysis of WRKY genes during leaf and corolla senescence of Petunia hybrida plants

Several families of transcription factors (TFs) control the progression of senescence. Many key TFs belonging to the WRKY family have been described to play crucial roles in the regulation of leaf senescence, mainly in Arabidopsis thaliana. However, little is known about senescence-associated WRKY members in floricultural species. Delay of senescence in leaves and petals of Petunia hybrida, a worldwide ornamental crop are highly appreciated traits. In this work, starting from 28 differentially expressed WRKY genes of A. thaliana during the progression of leaf senescence, we identified the orthologous in P. hybrida and explored the expression profiles of 20 PhWRKY genes during the progression of natural (age-related) leaf and corolla senescence as well as in the corollas of flowers undergoing pollination-induced senescence. Simultaneous visualization showed consistent and similar expression profiles of PhWRKYs during natural leaf and corolla senescence, although weak expression changes were observed during pollination-induced senescence. Comparable expression trends between PhWRKYs and the corresponding genes of A. thaliana were observed during leaf senescence, although more divergence was found in petals of pollinated petunia flowers. Integration of expression data with phylogenetics, conserved motif and cis-regulatory element analyses were used to establish a list of candidates that could regulate more than one senescence process. Our results suggest that several members of the WRKY family of TFs are tightly linked to the regulation of senescence in P. hybrida.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01243-y.

The genus Chrysanthemum: Phylogeny, biodiversity, phytometabolites, and chemodiversity

The ecologically and economically important genus Chrysanthemum contains around 40 species and many hybrids and cultivars. The dried capitulum of Chrysanthemum morifolium (CM) Ramat. Tzvel, i.e., Flos Chrysanthemi, is frequently used in traditional Chinese medicine (TCM) and folk medicine for at least 2,200 years. It has also been a popular tea beverage for about 2,000 years since Han Dynasty in China. However, the origin of different cultivars of CM and the phylogenetic relationship between Chrysanthemum and related Asteraceae genera are still elusive, and there is a lack of comprehensive review about the association between biodiversity and chemodiversity of Chrysanthemum. This article aims to provide a synthetic summary of the phylogeny, biodiversity, phytometabolites and chemodiversity of Chrysanthemum and related taxonomic groups, focusing on CM and its wild relatives. Based on extensive literature review and in light of the medicinal value of chrysanthemum, we give some suggestions for its relationship with some genera/species and future applications. Mining chemodiversity from biodiversity of Chrysanthemum containing subtribe Artemisiinae, as well as mining therapeutic efficacy and other utilities from chemodiversity/biodiversity, is closely related with sustainable conservation and utilization of Artemisiinae resources. There were eight main cultivars of Flos Chrysanthemi, i.e., Hangju, Boju, Gongju, Chuju, Huaiju, Jiju, Chuanju and Qiju, which differ in geographical origins and processing methods. Different CM cultivars originated from various hybridizations between multiple wild species. They mainly contained volatile oils, triterpenes, flavonoids, phenolic acids, polysaccharides, amino acids and other phytometabolites, which have the activities of antimicrobial, anti-viral, antioxidant, anti-aging, anticancer, anti-inflammatory, and closely related taxonomic groups could also be useful as food, medicine and tea. Despite some progresses, the genetic/chemical relationships among varieties, species and relevant genera have yet to be clarified; therefore, the roles of pharmacophylogeny and omics technology are highlighted.

A genome-wide association and fine-mapping study of white rust resistance in hexaploid chrysanthemum cultivars with a wild diploid reference genome

White rust caused by Puccinia horiana is one of the most serious diseases of chrysanthemum (Chrysanthemum × morifolium). In this study, we report the DNA markers associated with resistance against P. horiana via a simple approach using the genome of a wild diploid relative, Chrysanthemum seticuspe. First, we identified the important region of the genome in the resistant cultivar "Ariesu" via a genome-wide association study. Simplex single nucleotide polymorphism (SNP) markers mined from ddRAD-Seq were used in a biparental population originating from crosses between resistant "Ariesu" and susceptible "Yellow Queen". The C. seticuspe genome was used as a reference. For the fine mapping of P. horiana resistance locus 2 (Phr2), a comparative whole genome sequencing study was conducted. Although the genome sequences of chrysanthemum cultivars assembled via the short-read approach were fragmented, reliable genome alignments were reconstructed by mapping onto the chromosome level of the C. seticuspe pseudomolecule. Base variants were then identified by comparing the assembled genome sequences of resistant "Ariesu" and susceptible "Yellow Queen". Consequently, SNP markers that were closer to Phr2 compared with ddRAD-Seq markers were obtained. These SNP markers co-segregated with resistance in F₁ progenies originating from resistant "Ariesu" and showed robust transferability for detecting Phr2-conferring resistance among chrysanthemum genetic resources. The wild C. seticuspe pseudomolecule, a de facto monoploid genome used for ddRAD-Seq analysis and assembled genome sequence comparison, demonstrated this method's utility as a model for developing DNA markers in hexaploid chrysanthemum cultivars.

Genome-wide identification, phylogenetic analysis, and expression profiles of trihelix transcription factor family genes in quinoa (Chenopodium quinoa Willd.) under abiotic stress conditions

Background

The trihelix family of transcription factors plays essential roles in the growth, development, and abiotic stress response of plants. Although several studies have been performed on the trihelix gene family in several dicots and monocots, this gene family is yet to be studied in Chenopodium quinoa (quinoa).

Results

In this study, 47 C. quinoa trihelix (CqTH) genes were in the quinoa genome. Phylogenetic analysis of the CqTH and trihelix genes from Arabidopsis thaliana and Beta vulgaris revealed that the genes were clustered into five subfamilies: SIP1, GTγ, GT1, GT2, and SH4. Additionally, synteny analysis revealed that the CqTH genes were located on 17 chromosomes, with the exception of chromosomes 8 and 11, and 23 pairs of segmental duplication genes were detected. Furthermore, expression patterns of 10 CqTH genes in different plant tissues and at different developmental stages under abiotic stress and phytohormone treatment were examined. Among the 10 genes, CqTH02, CqTH25, CqTH18, CqTH19, CqTH25, CqTH31, and CqTH36, were highly expressed in unripe achenes 21 d after flowering and in mature achenes compared with other plant tissues. Notably, the 10 CqTH genes were upregulated in UV-treated leaves, whereas CqTH36 was consistently upregulated in the leaves under all abiotic stress conditions.

Conclusions

The findings of this study suggest that gene duplication could be a major driver of trihelix gene evolution in quinoa. These findings could serve as a basis for future studies on the roles of CqTH transcription factors and present potential genetic markers for breeding stress-resistant and high-yielding quinoa varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08726-y.

The chromosome-scale assembly of endive (Cichorium endivia) genome provides insights into the sesquiterpenoid biosynthesis

Endive (Cichorium endivia L.) is a leafy vegetable in the Asteraceae family. Sesquiterpene lactones (STLs) in endive leaves bring a bitter taste that varies between varieties. Despite their importance in breeding varieties with unique flavours, sesquiterpenoid biosynthesis pathways in endive are poorly understood. We assembled a chromosome-scale endive genome of 641 Mb with a contig N50 of 5.16 Mb and annotated 46,711 protein-coding genes. Several gene families, especially terpene synthases (TPS) genes, expanded significantly in the C. endivia genome. STLs biosynthesis-related genes and TPS genes in more bitter varieties have shown a higher level of expression, which could be attributed to genomic variations. Our results penetrate the origin and diversity of bitter taste and facilitate the molecular breeding of endive varieties with unique bitter tastes. The high-quality endive assembly would provide a reference genome for studying the evolution and diversity of Asteraceae.

The genome sequencing and comparative analysis of a wild kiwifruit Actinidia eriantha

The current kiwifruit industry is mainly based on the cultivars derived from the species Actinidia chinensis (Ac) which may bring risks such as canker disease. Introgression of desired traits from wild relatives is an important method for improving kiwifruit cultivars. Actinidia eriantha (Ae) is a particularly important taxon used for hybridization or introgressive breeding of new kiwifruit cultivars because of its valued species-specific traits. Here, we assembled a chromosome-scale high-quality genome of a Ae sample which was directly collected from its wild populations. Our analysis revealed that 41.3% of the genome consists of repetitive elements, comparable to the percentage in Ac and Ae cultivar "White" genomes. The genomic structural variation, including the presence/absence-variation (PAV) of genes, is distinct between Ae and Ac, despite both sharing the same two kiwifruit-specific whole genome duplication (WGD) events. This suggests that a post-WGD divergence mechanism occurred during their evolution. We further investigated genes involved in ascorbic acid biosynthesis and disease-resistance of Ae, and we found introgressive genome could contribute to the complex relationship between Ae and other representative kiwifruit taxa. Collectively, the Ae genome offers valuable genetic resource to accelerate kiwifruit breeding applications.

Comparative metabolomics analysis revealed biomarkers and distinct flavonoid biosynthesis regulation in Chrysanthemum mongolicum and C. rhombifolium

INTRODUCTION

Chrysanthemums are traditional flowers that originated in China and have high ornamental, economic and medicinal value. They are widely used as herbal remedies and consumed as food or beverages in folk medicine. However, little is known about their metabolic composition.

OBJECTIVES

The aims of this work were to determine the metabolic composition of and natural variation among different species of Chrysanthemum and to explore new potential resources for drug discovery and sustainable utilisation of wild Chrysanthemum.

METHODS

The metabolomes of Chrysanthemum mongolicum (Ling) Tzvel. and Chrysanthemum rhombifolium H. Ohashi & Yonek. were compared using a widely targeted metabolomics approach based on liquid chromatography-tandem mass spectrometry (LC-MS/MS).

RESULTS

In total, 477 metabolites were identified, of which 72 showed significant differences in expression between C. mongolicum and C. rhombifolium, mainly in flavonoids, organic acids and nucleotides. The flavone and flavonol biosynthesis pathway showed significant enrichment among the differentially expressed metabolites. The contents of genkwanin, trigonelline, diosmin, narcissoside, 3,4-dihydroxyphenylacetic acid, linarin, N',N'-p-coumarin, C-hexosyl-tricetin O-pentoside, chrysoeriol, acacetin and kaempferol-3-O-gentiobioside were significantly different between the two species and represent potential biomarkers.

CONCLUSION

The types of flavonoid-related metabolites in the flavonoid biosynthesis pathway differed between C. mongolicum and C. rhombifolium. The mechanisms underlying the unique adaptations of these two species to their environments may involve variations in the composition and abundance of flavonoids, organic acids, and nucleotides. These methods are promising to identify functional compounds in Chrysanthemum species and can provide potential resources for drug discovery and the sustainable utilisation of Chrysanthemum plants.

Characterization of the TCP Gene Family in Chrysanthemum nankingense and the Role of CnTCP4 in Cold Tolerance

Plant-specific TCP transcription factors play a key role in plant development and stress responses. Chrysanthemum nankingense shows higher cold tolerance than its ornamental polyploid counterpart. However, whether the TCP gene family plays a role in conferring cold tolerance upon C. nankingense remains unknown. Here, we identified 23 CnTCP genes in C. nankingense, systematically analyzed their phylogenetic relationships and synteny with TCPs from other species, and evaluated their expression profiles at low temperature. Phylogenetic analysis of the protein sequences suggested that CnTCP proteins fall into two classes and three clades, with a typical bHLH domain. However, differences between C. nankingense and Arabidopsis in predicted protein structure and binding sites suggested a unique function of CnTCPs in C. nankingense. Furthermore, expression profiles showed that expression of most CnTCPs were downregulated under cold conditions, suggesting their importance in plant responses to cold stress. Notably, expression of miR319 and of its predicted target genes, CnTCP2/4/14, led to fast responses to cold. Overexpression of Arabidopsis CnTCP4 led to hypersensitivity to cold, suggesting that CnTCP4 might play a negative role in C. nankingense responses to cold stress. Our results provide a foundation for future functional genomic studies on this gene family in chrysanthemum.

Chromosome-Scale Genome Assembly Provides Insights into Speciation of Allotetraploid and Massive Biomass Accumulation of Elephant Grass (Pennisetum purpureum Schum.)

Elephant grass (Pennisetum purpureum Schum) is an important forage, biofuels and industrial plant widely distributed in tropical and subtropical areas globally. It is characterized with robust growth and high biomass. We sequenced its allopolyploid genome and assembled 2.07 Gb into A' and B sub-genomes of 14 chromosomes with scaffold N50 of 8.47 Mb, yielding a total of 77,139 genes. The allotetraploid speciation occurred approximately 15 MYA after the divergence between Setaria italica and Pennisetum glaucum, according to a phylogenetic analysis of Pennisetum species. Double whole-genome duplication (WGD) and polyploidization events resulted in large scale gene expansion, especially in the key steps of growth and biomass accumulation. Integrated transcriptome profiling revealed the functional divergence between sub-genomes A' and B. A' sub-genome mainly contributed to plant growth, development and photosynthesis, whereas the B sub-genome was primarily responsible for effective transportation and resistance to stimulation. Some key gene families related to cellulose biosynthesis were expanded and highly expressed in stems, which could explain the high cellulose content in elephant grass. Our findings provide deep insights into genetic evolution of elephant grass and will aid future biological research and breeding, even for other grasses in the family Poaceae.

SSR Loci Analysis in Transcriptome and Molecular Marker Development in Polygonatum sibiricum

To study the SSR loci information and develop molecular markers, a total of 435,858 unigenes in transcriptome of Polygonatum sibiricum were used to explore SSR. The distribution frequency of SSR and the basic characteristics of repeat motifs were analyzed using MISA software, and SSR primers were designed by Primer 3.0 software and then validated by PCR. Moreover, the gene function analysis of SSR Unigene was obtained by Blast. The results showed that 112,728 SSR loci were found in the transcriptome of Polygonatum sibiricum, which distributed in 435,858 unigenes with a distribution frequency of 25.86%. Mo-nucleotide and Di-nucleotide repeat were the main types, accounted for 83.83% of all SSRs. The repeat motifs of A/T and AC/GT were the predominant repeat types of Mo-nucleotide and Di-nucleotide, respectively. A total of 113,305 pairs of SSR primers with the potential to produce polymorphism were designed for maker development. One hundred and fifty-four of the 500 randomly selected primers not only produced fragments with expected molecular size but also had high polymorphism, which could accurately separate the tested varieties. The gene function of unigenes containing SSR was mostly related to the molecular function of Polygonatum sibiricum. The SSR markers in transcriptome of Polygonatum sibiricum show rich type, strong specificity, and high potential of polymorphism, which will benefit the candidate gene mining and marker-assisted breeding. The developed markers can also provide technical methods for molecular identification of intraspecific species of Polygonatum Mill. and maker-assisted breeding of superior varieties of Polygonatum Mill.

Comparative genomics reveal the convergent evolution of CYP82D and CYP706X members related to flavone biosynthesis in Lamiaceae and Asteraceae

Distant species producing the same secondary metabolites is an interesting and common phenomenon in nature. A classic example of this is scutellarein whose derivatives have been used clinically for more than 30 years. Scutellarein occurs in significant amounts in species of two different orders, Scutellaria baicalensis and Erigeron breviscapus, which diverged more than 100 million years ago. Here, according to the genome-wide selection and functional identification of 39 CYP450 genes from various angiosperms, we confirmed that only seven Scutellaria-specific CYP82D genes and one Erigeron CYP706X gene could perform the catalytic activity of flavone 6-hydroxylase (F6H), suggesting that the convergent evolution of scutellarein production in these two distant species was caused by two independently evolved CYP450 families. We also identified seven Scutellaria-specific CYP82D genes encoding flavone 8-hydroxylase (F8H). The evolutionary patterns of CYP82 and CYP706 families via kingdom-wide comparative genomics highlighted the evolutionary diversity of CYP82D and the specificity of CYP706X in angiosperms. Multi-collinearity and phylogenetic analysis of CYP82D in Scutellaria confirmed that the function of F6H evolved from F8H. Furthermore, the SbaiCYP82D1^A319D , EbreCYP706X^R130A , EbreCYP706X^F312D and EbreCYP706X^A318D mutants can significantly decrease the catalytic activity of F6H, revealing the contribution of crucial F6H amino acids to the scutellarein biosynthesis of distant species. This study provides important insights into the multi-origin evolution of the same secondary metabolite biosynthesis in the plant kingdom.