Extensive sequence divergence between the reference genomes of Taraxacum kok-saghyz and Taraxacum mongolicum

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.

Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the Cretaceous Eupteleaceae

Living fossils are evidence of long-term sustained ecological success. However, whether living fossils have little molecular changes remains poorly known, particularly in plants. Here, we have introduced a novel method that integrates phylogenomic, comparative genomic, and ecological niche modeling analyses to investigate the rate of molecular evolution of Eupteleaceae, a Cretaceous relict angiosperm family endemic to East Asia. We assembled a high-quality chromosome-level nuclear genome, and the chloroplast and mitochondrial genomes of a member of Eupteleaceae (Euptelea pleiosperma). Our results show that Eupteleaceae is most basal in Ranunculales, the earliest-diverging order in eudicots, and shares an ancient whole-genome duplication event with the other Ranunculales. We document that Eupteleaceae has the slowest rate of molecular changes in the observed angiosperms. The unusually low rate of molecular evolution of Eupteleaceae across all three independent inherited genomes and genes within each of the three genomes is in association with its conserved genome architecture, ancestral woody habit, and conserved niche requirements. Our findings reveal the evolution and adaptation of living fossil plants through large-scale environmental change and also provide new insights into early eudicot diversification.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

Study of Dandelion (Taraxacum mongolicum Hand.-Mazz.) Salt Response and Caffeic Acid Metabolism under Saline Stress by Transcriptome Analysis

Utilizing salt-tolerant plants is a cost-effective strategy for agricultural production on salinized land. However, little is known about the mechanism of dandelion (Taraxacum mongolicum Hand.-Mazz.) in response to saline stress and caffeic acid biosynthesis. We investigated the morphological and physiological variations of two dandelions, namely, "BINPU2" (dandelion A) and "TANGHAI" (dandelion B) under gradient NaCl concentrations (0, 0.3%, 0.5%, 0.7%, and 0.9%), and analyzed potential mechanisms through a comparison analysis of transcriptomes in the two dandelions. Dandelion A had a high leaf weight; high ρ-coumaric acid, caffeic acid, ferulic acid, and caffeoyl shikimic acid contents; and high activities of POD and Pro. The maximum content of four kinds of phenolic acids mostly occurred in the 0.7% NaCl treatment. In this saline treatment, 2468 and 3238 differentially expressed genes (DEGs) in dandelion A and B were found, of which 1456 and 1369 DEGs in the two dandelions, respectively, showed up-regulation, indicating that more up-regulated DEGs in dandelion A may cause its high salt tolerance. Further, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that dandelion salt response and caffeic acid metabolism were mainly enriched in the phenylpropanoid biosynthesis pathway (ko00940) and response to ethylene (GO: 0009723). The caffeic acid biosynthesis pathway was reconstructed based on DEGs which were annotated to PAL, C4H, 4CL, HCT, C3'H, and CSE. Most of these genes showed a down-regulated mode, except for parts of DEGs of 4CL (TbA05G077650 and TbA07G073600), HCT (TbA03G009110, TbA03G009080, and novel.16880), and COMT (novel.13839). In addition, more up-regulated transcription factors (TFs) of ethylene TFs in dandelion A were found, but the TFs of ERF104, CEJ1, and ERF3 in the two dandelions under saline stress showed an opposite expression pattern. These up-regulated genes could enhance dandelion salt tolerance, and down-regulated DEGs in the caffeic acid biosynthesis pathway, especially CSE (TbA08G014310) and COMT (TbA04G07330), could be important candidate genes in the synthesis of caffeic acid under saline stress. The above findings revealed the potential mechanisms of salt response and caffeic acid metabolism in dandelion under saline stress, and provide references for salt-tolerant plant breeding and cultivation on saline-alkali land in the future.

Transcriptome analysis reveals the molecular mechanisms of rubber biosynthesis and laticifer differentiation during rubber seed germination

The molecular mechanisms underlying the initiation of natural rubber synthesis and laticifer differentiation have not been fully elucidated. In this study, we conducted a time-series transcriptome analysis of five rubber tree tissues at four stages of seed germination. A total of 161,199 DEGs were identified between the two groups, including most 16,673 DEGs (A3 vs B3 and A3 vs C3) and lest 1,210 DEGs (C2 vs D2). We found that the maturation of the seed is accompanied by the formation of laticifer cells in cotyledon. Meanwhile, the analysis of hormones related genes expression may provide effective clues for us to promote the differentiation of laticifer cells in seeds by hormones in the future. In this study, hormone-related gene enrichment analyses revealed that IAA, GA, and CTK were activated in laticifer containing tissues. Similarly, GO and GEGG analysis showed that hormone pathways, especially the auxin pathway, are enriched. Gene expression clustering was analyzed using the short time-series expression miner (STEM), and the analysis revealed four distinct trends in the gene expression profiles. Moreover, we enriched transcription factor (TF) enrichment in cotyledon and embryonic axis tissues, and the MYB type exhibited the most significant difference. Furthermore, our findings revealed that genes related to rubber synthesis exhibited tissue-specific expression patterns during seed germination. Notably, key genes associated with rubber biosynthesis, specifically small rubber particle protein (SRPP) and cis-prenyltransferase (CPT), exhibited significant changes in expression in cotyledon and embryonic axis tissues, suggesting synchronous rubber synthesis with seed germination. Our staining results reveled that laticifer cells were exits in the cotyledon before seed imbibition stage. In conclusion, these results lay the foundation for exploring the molecular mechanisms underlying laticifer differentiation and rubber synthesis during seed germination, deepening our understanding of the initiation stages of rubber biosynthesis and laticifer differentiation.

Advances in Genome Sequencing and Natural Rubber Biosynthesis in Rubber-Producing Plants

Natural rubber (cis-1,4-polyisoprene, NR) is an important raw material utilized widely in the manufacturing of medical, agricultural, and industrial products. Rubber tree (Hevea brasiliensis) and several alternative rubber-producing plants (Taraxacum kok-saghyz, Lactuca sativa, and Parthenium argentatum) have the capability to produce high-quality NR. With the progress of genome sequencing, similar rubber biosynthesis pathways have been discovered among different rubber-producing plant species. NR is synthesized and stored in rubber particles, which are specialized organelles comprising a hydrophobic NR core surrounded by a lipid monolayer and membrane-bound proteins. The rubber transferase complex is considered to be the pivotal enzyme involved in catalyzing NR biosynthesis. However, the exact compositions of the RT complex in rubber-producing plants remain elusive and poorly understood. Here, we review the progress of genome sequencing, natural rubber biosynthesis, and the components of the RT complex in rubber-producing plants. We emphasize that identifying the detailed components of the RT complex holds great significance for exploring the mechanism of NR biosynthesis and accelerating molecular breeding in rubber-producing plants.

Identification of miRNAs and their targets in two Taraxacum species with contrasting rubber-producing ability

MicroRNAs (miRNAs) are widely involved in various aspects of plant growth and development. However, how miRNAs and their targets regulate natural rubber metabolism remains unclear in the rubber-producing dandelions, which are being developed as alternative commercial sources of natural rubber. Here, we combined small RNA sequencing, degradome sequencing, target gene prediction, and mRNA sequencing to identify miRNAs and their targets in two dandelion species, the high rubber-yielding Taraxacum kok-saghyz (Tk) and the low rubber-yielding T. spadiceum (Ts). A total of 142 miRNAs, including 108 known and 34 novel ones, were discovered, with 53 identified as differentially expressed (DE) between the latex of Tk and Ts. Degradome sequencing identified 145 targets corresponding to 74 miRNAs. TAPIR and psRNATarget, respectively, predicted 165 and 164 non-redundant targets for the 53 aforementioned DE miRNAs. Gene ontology (GO) enrichment analysis indicated the DE miRNAs and their targets might affect natural rubber production via regulating macromolecular biosynthesis and metabolism in latex. Four critical types of regulatory modules, including miR172-AP2/ERF, miR164-NAC, miR160-ARF, and miRN19-protein kinase, were identified and their interaction networks were constructed, indicating a potential involvement in natural rubber production. The findings and the large miRNA dataset presented here are beneficial to further deciphering the roles of miRNAs in the biosynthesis of natural rubber and medicinal metabolites in dandelion.

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Large-scale long terminal repeat insertions produced a significant set of novel transcripts in cotton

Genomic analysis has revealed that the 1,637-Mb Gossypium arboreum genome contains approximately 81% transposable elements (TEs), while only 57% of the 735-Mb G. raimondii genome is occupied by TEs. In this study, we investigated whether there were unknown transcripts associated with TE or TE fragments and, if so, how these new transcripts were evolved and regulated. As sequence depths increased from 4 to 100 G, a total of 10,284 novel intergenic transcripts (intergenic genes) were discovered. On average, approximately 84% of these intergenic transcripts possibly overlapped with the long terminal repeat (LTR) insertions in the otherwise untranscribed intergenic regions and were expressed at relatively low levels. Most of these intergenic transcripts possessed no transcription activation markers, while the majority of the regular genic genes possessed at least one such marker. Genes without transcription activation markers formed their+1 and -1 nucleosomes more closely (only (117±1.4)bp apart), while twice as big spaces (approximately (403.5±46.0) bp apart) were detected for genes with the activation markers. The analysis of 183 previously assembled genomes across three different kingdoms demonstrated systematically that intergenic transcript numbers in a given genome correlated positively with its LTR content. Evolutionary analysis revealed that genic genes originated during one of the whole-genome duplication events around 137.7 million years ago (MYA) for all eudicot genomes or 13.7 MYA for the Gossypium family, respectively, while the intergenic transcripts evolved around 1.6 MYA, resultant of the last LTR insertion. The characterization of these low-transcribed intergenic transcripts can facilitate our understanding of the potential biological roles played by LTRs during speciation and diversifications.

Genome-Wide Identification of MADS-Box Genes in Taraxacum kok-saghyz and Taraxacum mongolicum: Evolutionary Mechanisms, Conserved Functions and New Functions Related to Natural Rubber Yield Formation

MADS-box transcription regulators play important roles in plant growth and development. However, very few MADS-box genes have been isolated in the genus Taraxacum, which consists of more than 3000 species. To explore their functions in the promising natural rubber (NR)-producing plant Taraxacum kok-saghyz (TKS), MADS-box genes were identified in the genome of TKS and the related species Taraxacum mongolicum (TM; non-NR-producing) via genome-wide screening. In total, 66 TkMADSs and 59 TmMADSs were identified in the TKS and TM genomes, respectively. From diploid TKS to triploid TM, the total number of MADS-box genes did not increase, but expansion occurred in specific subfamilies. Between the two genomes, a total of 11 duplications, which promoted the expansion of MADS-box genes, were identified in the two species. TkMADS and TmMADS were highly conserved, and showed good collinearity. Furthermore, most TkMADS genes exhibiting tissue-specific expression patterns, especially genes associated with the ABCDE model, were preferentially expressed in the flowers, suggesting their conserved and dominant functions in flower development in TKS. Moreover, by comparing the transcriptomes of different TKS lines, we identified 25 TkMADSs related to biomass formation and 4 TkMADSs related to NR content, which represented new targets for improving the NR yield of TKS.

Phylogenomic analysis provides insights into MADS-box and TCP gene diversification and floral development of the Asteraceae, supported by de novo genome and transcriptome sequences from dandelion (Taraxacum officinale)

The Asteraceae is the largest angiosperm family with more than 25,000 species. Individual studies have shown that MADS-box and TCP transcription factors are regulators of the development and symmetry of flowers, contributing to their iconic flower-head (capitulum) and floret. However, a systematic study of MADS-box and TCP genes across the Asteraceae is lacking. We performed a comparative analysis of genome sequences of 33 angiosperm species including our de novo assembly of diploid sexual dandelion (Taraxacum officinale) and 11 other Asteraceae to investigate the lineage-specific evolution of MADS-box and TCP genes in the Asteraceae. We compared the phylogenomic results of MADS-box and TCP genes with their expression in T. officinale floral tissues at different developmental stages to demonstrate the regulation of genes with Asteraceae-specific attributes. Here, we show that MADS-box MIKC^c and TCP-CYCLOIDEA (CYC) genes have expanded in the Asteraceae. The phylogenomic analysis identified AGAMOUS-like (AG-like: SEEDSTICK [STK]-like), SEPALATA-like (SEP3-like), and TCP-PROLIFERATING CELL FACTOR (PCF)-like copies with lineage-specific genomic contexts in the Asteraceae, Cichorioideae, or dandelion. Different expression patterns of some of these gene copies suggest functional divergence. We also confirm the presence and revisit the evolutionary history of previously named “Asteraceae-Specific MADS-box genes (AS-MADS).” Specifically, we identify non-Asteraceae homologs, indicating a more ancient origin of this gene clade. Syntenic relationships support that AS-MADS is paralogous to FLOWERING LOCUS C (FLC) as demonstrated by the shared ancient duplication of FLC and SEP3.

Molecular Genetic Research and Genetic Engineering of Taraxacum kok-saghyz L.E. Rodin

Natural rubber (NR) remains an indispensable raw material with unique properties that is used in the manufacture of a large number of products and the global demand for it is growing every year. The only industrially important source of NR is the tropical tree Hevea brasiliensis (Willd. ex A.Juss.) Müll.Arg., thus alternative sources of rubber are required. For the temperate zone, the most suitable source of high quality rubber is the Russian (Kazakh) dandelion Taraxacum kok-saghyz L.E. Rodin (TKS). An obstacle to the widespread industrial cultivation of TKS is its high heterozygosity, poor growth energy, and low competitiveness in the field, as well as inbreeding depression. Rapid cultivation of TKS requires the use of modern technologies of marker-assisted and genomic selection, as well as approaches of genetic engineering and genome editing. This review is devoted to describing the progress in the field of molecular genetics, genomics, and genetic engineering of TKS. Sequencing and annotation of the entire TKS genome made it possible to identify a large number of SNPs, which were subsequently used in genotyping. To date, a total of 90 functional genes have been identified that control the rubber synthesis pathway in TKS. The most important of these proteins are part of the rubber transferase complex and are encoded by eight genes for cis-prenyltransferases (TkCPT), two genes for cis-prenyltransferase-like proteins (TkCPTL), one gene for rubber elongation factor (TkREF), and nine genes for small rubber particle proteins (TkSRPP). In TKS, genes for enzymes of inulin metabolism have also been identified and genome-wide studies of other gene families are also underway. Comparative transcriptomic and proteomic studies of TKS lines with different accumulations of NR are also being carried out, which help to identify genes and proteins involved in the synthesis, regulation, and accumulation of this natural polymer. A number of authors already use the knowledge gained in the genetic engineering of TKS and the main goal of these works is the rapid transformation of the TKS into an economically viable rubber crop. There are no great successes in this area so far, therefore work on genetic transformation and genome editing of TKS should be continued, considering the recent results of genome-wide studies.

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

Genome-wide analysis of the SWEET genes in Taraxacum kok-saghyz Rodin: An insight into two latex-abundant isoforms

Taraxacum kok-saghyz Rodin (Tk) is a promising alternative rubber-producing grass. However, low biomass and rubber-producing capability limit its commercial application. As a carbon source transporter in plants, sugar will eventually be exported transporters (SWEETs) have been reported to play pivotal roles in diverse physiological events in the context of carbon assimilate transport and utilization. Theoretically, SWEETs would participate in Tk growth, development and response to environmental cues with relation to the accumulation of rubber and biomass, both of which rely on the input of carbon assimilates. Here, we identified 22 TkSWEETs through homology searching of the Tk genomes and bioinformatics analyses. RNA-seq and qRT-PCR analysis revealed these TkSWEETs to have overlapping yet distinct tissue expression patterns. Two TkSWEET isofroms, TkSWEET1 and TkSWEET12 expressed substantially in the latex, the cytoplasm of rubber-producing laticifers as well as the rubber source. As revealed by the transient expression analysis using Tk mesophyll protoplasts, both TkSWEET1 and TkSWEET12 were located in the plasma membrane. Heterologous expressions of the two TkSWEETs in a yeast mutant revealed that only TkSWEET1 exhibited apparent sugar transport activities, with a preference for monosaccharides. Interestingly, TkSWEET12, the latex-predominant TkSWEET isoform, seemed to have evolved from a tandem duplication event that results in a cluster of six TkSWEET genes with the TkSWEET12 therein, suggesting its specialized roles in the laticifers.

The integrated genomics of crop domestication and breeding

As a major event in human civilization, wild plants were successfully domesticated to be crops, largely owing to continuing artificial selection. Here, we summarize new discoveries made during the past decade in crop domestication and breeding. The construction of crop genome maps and the functional characterization of numerous trait genes provide foundational information. Approaches to read, interpret, and write complex genetic information are being leveraged in many plants for highly efficient de novo or re-domestication. Understanding the underlying mechanisms of crop microevolution and applying the knowledge to agricultural productions will give possible solutions for future challenges in food security.

Phenols and terpenoids: natural products as inhibitors of NLRP3 inflammasome in cardiovascular diseases

Inflammatory infiltration has been implicated in the pathogenesis of cardiovascular diseases (CVDs). The NLRP3 inflammasome is involved in the development of several types of CVDs, including myocardial infarction, myocardial ischemia-reperfusion damage, heart failure, atrial fibrillation, and hypertension. Inhibiting the activity of NLRP3 inflammasome can inhibit the progress of CVDs. However, there is no NLRP3 inflammasome inhibitor in clinic, and it is very important to find a safe and effective NLRP3 inhibitor. Phenols and terpenoids are naturally natural products that have many anti-inflammatory effects in CVDs by modulating the NLRP3 inflammatory pathway. Thus, 20 natural products from phenols and terpenoids for the treatment of cardiovascular disease based on the inhibition of NLRP3 inflammasome were summarized and screened. Docking results showed salvianolic acid B and ellagic acid in phenols, and oridonin and triptolide in terpenoids had a better binding activity with NLRP3, which can provide theoretical support for finding novel NLRP3 inflammasome inhibitors or lead compounds in the future.