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.

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.

Characterization of natural rubber concerning its components and molecular weight in Taraxacum kok-saghyz Rodin using pyrolysis gas chromatography-mass spectrometry

Taraxacum kok-saghyz Rodin (TKS) has abundant natural rubber in its root and the molecular weight of its natural rubber is higher than that in Hevea brasiliensis. Thus, TKS is an excellent alternative for the commercial production of natural rubber. The content and molecular weight of natural rubber are two qualitative indicators. Efficient determination for both indicators is still a challenge. In this study, we developed a method to simultaneously determine the content and molecular weight of natural rubber in TKS with pyrolysis－gas chromatography-mass spectrometry. The content of natural rubber was quantified by internal standard method. We optimized the pyrolysis temperature and chromatographic method during content determination. The limits of detection and quantification were 0.47 and 1.56 μg, respectively. In addition, the arachidonic acid methyl ester, an unsaturated fatty acid proposed from the α-end group of natural rubber, was quantified to obtain the number of natural rubber polymers. Based on the content and the polymer number, we also quantified the molecular weight of natural rubber. Thus, the content and molecular weight of natural rubber were simultaneously determined in TKS. Our study provides a new perspective for the high throughput analysis of natural rubber.

Overexpression of a pseudo-etiolated-in-light-like protein in Taraxacum koksaghyz leads to a pale green phenotype and enables transcriptome-based network analysis of photomorphogenesis and isoprenoid biosynthesis

Introduction

Plant growth and greening in response to light require the synthesis of photosynthetic pigments such as chlorophylls and carotenoids, which are derived from isoprenoid precursors. In Arabidopsis, the pseudo-etiolated-in-light phenotype is caused by the overexpression of repressor of photosynthetic genes 2 (RPGE2), which regulates chlorophyll synthesis and photosynthetic genes.

Methods

We investigated a homologous protein in the Russian dandelion (Taraxacum koksaghyz) to determine its influence on the rich isoprenoid network in this species, using a combination of in silico analysis, gene overexpression, transcriptomics and metabolic profiling.

Results

Homology-based screening revealed a gene designated pseudo-etiolated-in-light-like (TkPEL-like), and in silico analysis identified a light-responsive G-box element in its promoter. TkPEL-like overexpression in dandelion plants and other systems reduced the levels of chlorophylls and carotenoids, but this was ameliorated by the mutation of one or both conserved cysteine residues. Comparative transcriptomics in dandelions overexpressing TkPEL-like showed that genes responsible for the synthesis of isoprenoid precursors and chlorophyll were downregulated, probably explaining the observed pale green leaf phenotype. In contrast, genes responsible for carotenoid synthesis were upregulated, possibly in response to feedback signaling. The evaluation of additional differentially expressed genes revealed interactions between pathways.

Discussion

We propose that TkPEL-like negatively regulates chlorophyll- and photosynthesis-related genes in a light-dependent manner, which appears to be conserved across species. Our data will inform future studies addressing the regulation of leaf isoprenoid biosynthesis and photomorphogenesis and could be used in future breeding strategies to optimize selected plant isoprenoid profiles and generate suitable plant-based production platforms.

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.

Physiological and Transcriptomic Analyses Reveal the Effects of Carbon-Ion Beam on Taraxacum kok-saghyz Rodin Adventitious Buds

Taraxacum kok-saghyz Rodin (TKS) has great potential as an alternative natural-rubber (NR)-producing crop. The germplasm innovation of TKS still faces great challenges due to its self-incompatibility. Carbon-ion beam (CIB) irradiation is a powerful and non-species-specific physical method for mutation creation. Thus far, the CIB has not been utilized in TKS. To better inform future mutation breeding for TKS by the CIB and provide a basis for dose-selection, adventitious buds, which not only can avoid high levels of heterozygosity, but also further improve breeding efficiency, were irradiated here, and the dynamic changes of the growth and physiologic parameters, as well as gene expression pattern were profiled, comprehensively. The results showed that the CIB (5-40 Gy) caused significant biological effects on TKS, exhibiting inhibitory effects on the fresh weight and the number of regenerated buds and roots. Then,15 Gy was chosen for further study after comprehensive consideration. CIB-15 Gy resulted in significant oxidative damages (hydroxyl radical (OH^•) generation activity, 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical-scavenging activity and malondialdehyde (MDA) content) and activated the antioxidant system (superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), and ascorbate peroxidase (APX)) of TKS. Based on RNA-seq analysis, the number of differentially expressed genes (DEGs) peaked at 2 h after CIB irradiation. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that DNA-replication-/repair- (mainly up-regulated), cell-death- (mainly up-regulated), plant-hormone- (auxin and cytokinin, which are related to plant morphogenesis, were mainly down-regulated), and photosynthesis- (mainly down-regulated) related pathways were involved in the response to the CIB. Furthermore, CIB irradiation can also up-regulate the genes involved in NR metabolism, which provides an alternative strategy to elevate the NR production in TKS in the future. These findings are helpful to understand the radiation response mechanism and further guide the future mutation breeding for TKS by the CIB.

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.

Transcriptome analysis of Taraxacum kok-saghyz reveals the role of exogenous methyl jasmonate in regulating rubber biosynthesis and drought tolerance

Taraxacum kok-saghyz has been identified as one of the most promising alternative rubber crops, with laticifer cells that produce high-quality rubber. To uncover the underlying molecular mechanisms regulating natural rubber biosynthesis under MeJA induction, a reference transcriptome was constructed from nine samples of T. kok-saghyz. MeJA treatment was applied for 0 h (control), 6 h, and 24 h. A total of 7452 differentially expressed genes (DEGs) were identified in response to MeJA stress, relative to the control. Functional enrichment showed that these DEGs were primarily related to hormone signaling, defensive responses, and secondary metabolism. Combined analysis of the DEGs induced by MeJA and high-expression genes in laticifer cells further identified seven DEGs related to natural rubber biosynthesis that were upregulated in latex tissue, suggesting that these candidate genes could prove valuable in studying the mechanism of MeJA-mediated natural rubber biosynthesis. In addition, 415 MeJA-responsive DEGs were from several transcription factor families associated with drought resistance. This study helps to elucidate the mechanism of natural rubber biosynthesis in T. kok-saghyz in response to MeJA stress and identifies key candidate MeJA-induced DEGs in laticifer tissue, as well as a candidate drought-response target gene, whose knowledge will promote the breeding of T. kok-saghyz in the aspect of rubber yields and quality, and drought tolerance.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

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.

Functional Characterization of TkSRPP Promoter in Response to Hormones and Wounding Stress in Transgenic Tobacco

Taraxacum kok-saghyz is a model species for studying natural rubber biosynthesis because its root can produce high-quality rubber. Small rubber particle protein (SRPP), a stress-related gene to multiple stress responses, involves in natural rubber biosynthesis. To investigate the transcriptional regulation of the TkSRPP promoter, the full-length promoter PR0 (2188 bp) and its four deletion derivatives, PR1 (1592 bp), PR2 (1274 bp), PR3 (934 bp), and PR4 (450 bp), were fused to β-glucuronidase (GUS) reporter gene and transformed into tobacco. The GUS tissue staining showed that the five promoters distinctly regulated GUS expression utilizing transient transformation of tobacco. The GUS activity driven by a PR0 promoter was detected in transgenic tobacco leaves, stem and roots, suggesting that the TkSRPP promoter was not tissue-specific. Deletion analyses in transgenic tobacco have demonstrated that the PR3 from -934 bp to -450 bp core region responded strongly to the hormones, methyl jasmonate (MeJA), abscisic acid (ABA), and salicylic acid (SA), and also to injury induction. The TkSRPP gene was highly expressed under hormones and wound-induced conditions. This study reveals the regulation pattern of the SRPP promoter, and provides valuable information for studying natural rubber biosynthesis under hormones and wounding stress.

Natural rubber reduces herbivory and alters the microbiome below ground

Laticifers are hypothesized to mediate both plant-herbivore and plant-microbe interactions. However, there is little evidence for this dual function. We investigated whether the major constituent of natural rubber, cis-1,4-polyisoprene, a phylogenetically widespread and economically important latex polymer, alters plant resistance and the root microbiome of the Russian dandelion (Taraxacum koksaghyz) under attack of a root herbivore, the larva of the May cockchafer (Melolontha melolontha). Rubber-depleted transgenic plants lost more shoot and root biomass upon herbivory than normal rubber content near-isogenic lines. Melolontha melolontha preferred to feed on artificial diet supplemented with rubber-depleted rather than normal rubber content latex. Likewise, adding purified cis-1,4-polyisoprene in ecologically relevant concentrations to diet deterred larval feeding and reduced larval weight gain. Metagenomics and metabarcoding revealed that abolishing biosynthesis of natural rubber alters the structure but not the diversity of the rhizosphere and root microbiota (ecto- and endophytes) and that these changes depended on M. melolontha damage. However, the assumption that rubber reduces microbial colonization or pathogen load is contradicted by four lines of evidence. Taken together, our data demonstrate that natural rubber biosynthesis reduces herbivory and alters the plant microbiota, which highlights the role of plant-specialized metabolites and secretory structures in shaping multitrophic interactions.

A jasmonate-responsive bHLH transcription factor TaMYC2 positively regulates triterpenes biosynthesis in Taraxacum antungense Kitag

Dandelion is a well-known traditional medical herb, also used as functional food. Dandelion possesses many medical properties, such as anti-bacterial and antioxidant activity and contains a variety of triterpenes, such as α-amyrin, β-amyrin, taraxerol and taraxasterol. In this study, we found that triterpenes biosynthesis was promoted by methyl jasmonate (MeJA), while the transcriptional mechanism underlying triterpenes biosynthesis was rarely investigated. Here, a MeJA-induced bHLH transcription factor TaMYC2 was identified. The content of taraxasterol and taraxerol in dandelion was obviously enhanced in overexpression TaMYC2 transgenic lines and expression level of the squalene synthase gene (TaSS) was elevated to about 3-5 folds compared with the control lines. Dual-LUC, Y1H and EMSA experiments revealed that TaMYC2 bound to the E-box motif in the promoter of TaSS and activated its transcription. Taken together, this study suggested that TaMYC2 acted as a positive regulator for bioengineering approaches to produce high content triterpenes-producing dandelions.

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.

Cut-dip-budding delivery system enables genetic modifications in plants without tissue culture

Of the more than 370 000 species of higher plants in nature, fewer than 0.1% can be genetically modified due to limitations of the current gene delivery systems. Even for those that can be genetically modified, the modification involves a tedious and costly tissue culture process. Here, we describe an extremely simple cut-dip-budding (CDB) delivery system, which uses Agrobacterium rhizogene to inoculate explants, generating transformed roots that produce transformed buds due to root suckering. We have successfully used CDB to achieve the heritable transformation of plant species in multiple plant families, including two herbaceous plants (Taraxacum kok-saghyz and Coronilla varia), a tuberous root plant (sweet potato), and three woody plant species (Ailanthus altissima, Aralia elata, and Clerodendrum chinense). These plants have previously been difficult or impossible to transform, but the CDB method enabled efficient transformation or gene editing in them using a very simple explant dipping protocol, under non-sterile conditions and without the need for tissue culture. Our work suggests that large numbers of plants could be amenable to genetic modifications using the CDB method.

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.

Identification of the WRKY Gene Family and Characterization of Stress-Responsive Genes in Taraxacum kok-saghyz Rodin

WRKY transcription factors present unusual research value because of their critical roles in plant physiological processes and stress responses. Taraxacum kok-saghyz Rodin (TKS) is a perennial herb of dandelion in the Asteraceae family. However, the research on TKS WRKY TFs is limited. In this study, 72 TKS WRKY TFs were identified and named. Further comparison of the core motifs and the structure of the WRKY motif was analyzed. These TFs were divided into three groups through phylogenetic analysis. Genes in the same group of TkWRKY usually exhibit a similar exon-intron structure and motif composition. In addition, virtually all the TKS WRKY genes contained several cis-elements related to stress response. Expression profiling of the TkWRKY genes was assessed using transcriptome data sets and Real-Time RT-PCR data in tissues during physiological development, under abiotic stress and hormonal treatments. For instance, the TkWRKY18, TkWRKY23, and TkWRKY38 genes were significantly upregulated during cold stress, whereas the TkWRKY21 gene was upregulated under heat-stress conditions. These results could provide a basis for further studies on the function of the TKS WRKY gene family and genetic amelioration of TKS germplasm.

Genome-Wide Identification and Spatial Expression Analysis of Histone Modification Gene Families in the Rubber Dandelion Taraxacum kok-saghyz

Taraxacum kok-saghyz (Tks), also known as the Russian dandelion, is a recognized alternative source of natural rubber quite comparable, for quality and use, to the one obtained from the so-called rubber tree, Hevea brasiliensis. In addition to that, Tks roots produce several other compounds, including inulin, whose use in pharmaceutical and dietary products is quite extensive. Histone-modifying genes (HMGs) catalyze a series of post-translational modifications that affect chromatin organization and conformation, which, in turn, regulate many downstream processes, including gene expression. In this study, we present the first analysis of HMGs in Tks. Altogether, we identified 154 putative Tks homologs: 60 HMTs, 34 HDMs, 42 HATs, and 18 HDACs. Interestingly, whilst most of the classes showed similar numbers in other plant species, including M. truncatula and A. thaliana, HATs and HMT-PRMTs were indeed more abundant in Tks. Composition and structure analysis of Tks HMG proteins showed, for some classes, the presence of novel domains, suggesting a divergence from the canonical HMG model. The analysis of publicly available transcriptome datasets, combined with spatial expression of different developmental tissues, allowed us to identify several HMGs with a putative role in metabolite biosynthesis. Overall, our work describes HMG genomic organization and sets the premises for the functional characterization of epigenetic modifications in rubber-producing plants.

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.

Genome-wide identification of CBF genes and their responses to cold acclimation in Taraxacum kok-saghyz

C-repeat binding factors (CBFs) are transcription factors that are known to play important roles in plant cold acclimation. They are highly conserved in most higher plants. Taraxacum kok-saghyz (TKS) is an herb native to China and Kazakhstan and is well-known for its production of rubber silk with industrial and economic value. To understand cold acclimation mechanisms, we conducted a genome-wide discovery of the CBF family genes in TKS and revealed ten CBF genes. A bioinformatic analysis of the CBF genes was carried out to analyze the phylogenetic relationship, protein conservative motifs, protein physicochemical properties, gene structure, promoter cis-acting elements, and the gene expression patterns under cold acclimation and control conditions. It was found that most of these genes were highly responsive at the late stage of cold acclimation, indicating that they play important roles in the cold acclimation processes of TKS. This study provides a theoretical basis for the study of the molecular functions of the CBF gene family in TKS, and a useful guidance for the genetic improvement of the cold tolerance traits of TKS and other plants, including crops.

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

Plants belonging to the genus Taraxacum are widespread all over the world, which contain rubber-producing and non-rubber-producing species. However, the genomic basis underlying natural rubber (NR) biosynthesis still needs more investigation. Here, we presented high-quality genome assemblies of rubber-producing T. kok-saghyz TK1151 and non-rubber-producing T. mongolicum TM5. Comparative analyses uncovered a large number of genetic variations, including inversions, translocations, presence/absence variations, as well as considerable protein divergences between the two species. Two gene duplication events were found in these two Taraxacum species, including one common ancestral whole-genome triplication and one subsequent round of gene amplification. In genomes of both TK1151 and TM5, we identified the genes encoding for each step in the NR biosynthesis pathway and found that the SRPP and CPT gene families have experienced a more obvious expansion in TK1151 compared to TM5. This study will have large-ranging implications for the mechanism of NR biosynthesis and genetic improvement of NR-producing crops.

Comparative Transcriptome Analysis in Taraxacum koksaghyz to Identify Genes that Determine Root Volume and Root Length

The Russian dandelion (Taraxacum koksaghyz, family Asteraceae) produces large amounts of natural rubber in the laticifers of its roots. This species has been proposed as an alternative source of natural rubber to augment or partly replace the rubber tree (Hevea brasiliensis) but domestication would require genetic improvement to increase rubber yields and agronomic optimization to facilitate harvesting and processing. Optimization has focused thus far on the size and shape of the roots, the primary storage organ for natural rubber and inulin. However, the corresponding genetic factors are poorly understood. Here we describe the comparative transcriptomic analysis of root tissues from T. koksaghyz plant sets featuring different root sizes and shapes, aiming to identify differentially expressed genes correlating with root length or root diameter in the upper root and root tip. The resulting datasets revealed multiple candidate genes for each trait and root part, including a glucan endo-1,3-β-d-glucosidase, an allene oxide synthase 3, and a TIFY10A/JAZ1 homolog. These three genes were tested by qRT-PCR in outdoor-grown plants with diverse root morphology, and the expression of two genes correlated with the appropriate root morphotype, confirming the effectiveness of our method. We evaluated the candidate genes to gain insight into their potential functions in root development. Such candidate genes could be suitable for marker-assisted breeding programs in the future.

The basic helix-loop-helix transcription factor TabHLH1 increases chlorogenic acid and luteolin biosynthesis in Taraxacum antungense Kitag

Polyphenols are the main active components of the anti-inflammatory compounds in dandelion, and chlorogenic acid (CGA) is one of the primary polyphenols. However, the molecular mechanism underlying the transcriptional regulation of CGA biosynthesis remains unclear. Hydroxycinnamoyl-CoA:quinate hydroxycinnamoyl transferase (HQT2) is the last rate-limiting enzyme in chlorogenic acid biosynthesis in Taraxacum antungense. Therefore, using the TaHQT2 gene promoter as a probe, a yeast one-hybrid library was performed, and a basic helix-loop-helix (bHLH) transcription factor, TabHLH1, was identified that shared substantial homology with Gynura bicolor DC bHLH1. The TabHLH1 transcript was highly induced by salt stress, and the TabHLH1 protein was localized in the nucleus. CGA and luteolin concentrations in TabHLH1-overexpression transgenic lines were significantly higher than those in the wild type, while CGA and luteolin concentrations in TabHLH1-RNA interference (RNAi) transgenic lines were significantly lower. Quantitative real-time polymerase chain reaction demonstrated that overexpression and RNAi of TabHLH1 in T. antungense significantly affected CGA and luteolin concentrations by upregulating or downregulating CGA and luteolin biosynthesis pathway genes, especially TaHQT2, 4-coumarate-CoA ligase (Ta4CL), chalcone isomerase (TaCHI), and flavonoid-3'-hydroxylase (TaF3'H). Dual-luciferase, yeast one-hybrid, and electrophoretic mobility shift assays indicated that TabHLH1 directly bound to the bHLH-binding motifs of proTaHQT2 and proTa4CL. This study suggests that TabHLH1 participates in the regulatory network of CGA and luteolin biosynthesis in T. antungense and might be useful for metabolic engineering to promote plant polyphenol biosynthesis.

Sepal Identity of the Pappus and Floral Organ Development in the Common Dandelion (Taraxacum officinale; Asteraceae)

The dry one-seeded fruits (cypselae) of the Asteraceae are often crowned with a pappus, an appendage of hairs or scales that assists in dispersal. It is generally assumed, but little investigated, that the pappus represents the outer floral whorl where the sepals are usually located. We analysed pappus–sepal homology in dandelions using micromorphological and floral gene expression analyses. We show that the pappus initiates from a ring primordium at the base of the corolla, heterochronic to the petals. Pappus parts form from this ring, with those in the alternipetalaous position usually being ahead in growth, referring to sepal identity. Tof-APETALLA1 expression increased during floret development and was highest in mature pappus. Tof-PISTILLATA expression was high and confined to the floral tissues containing the petals and stamens, consistent with expectations for sepals. Apart from the pappus, the dispersal structure of dandelion consists of the upper part of the fruit, the beak, which originates from the inner floral whorl. Thus, our results support the homology of the pappus with the sepals, but show that it is highly derived. Together with our floral stage definitions and verified qPCR reference genes, they provide a basis for evolution and development studies in dandelions and possibly other Asteraceae.

Genome Warehouse: A Public Repository Housing Genome-scale Data

The Genome Warehouse (GWH) is a public repository housing genome assembly data for a wide range of species and delivering a series of web services for genome data submission, storage, release, and sharing. As one of the core resources in the National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB; https://ngdc.cncb.ac.cn), GWH accepts both full and partial (chloroplast, mitochondrion, and plasmid) genome sequences with different assembly levels, as well as an update of existing genome assemblies. For each assembly, GWH collects detailed genome-related metadata of biological project, biological sample, and genome assembly, in addition to genome sequence and annotation. To archive high-quality genome sequences and annotations, GWH is equipped with a uniform and standardized procedure for quality control. Besides basic browse and search functionalities, all released genome sequences and annotations can be visualized with JBrowse. By May 21, 2021, GWH has received 19,124 direct submissions covering a diversity of 1108 species and has released 8772 of them. Collectively, GWH serves as an important resource for genome-scale data management and provides free and publicly accessible data to support research activities throughout the world. GWH is publicly accessible at https://ngdc.cncb.ac.cn/gwh.

Genetic transformation technologies for the common dandelion, Taraxacum officinale

BACKGROUND

Taraxacum officinale, or the common dandelion, is a widespread perennial species recognized worldwide as a common lawn and garden weed. Common dandelion is also cultivated for use in teas, as edible greens, and for use in traditional medicine. It produces latex and is closely related to the Russian dandelion, T. kok-saghyz, which is being developed as a rubber crop. Additionally, the vast majority of extant common dandelions reproduce asexually through apomictically derived seeds- an important goal for many major crops in modern agriculture. As such, there is increasing interest in the molecular control of important pathways as well as basic molecular biology and reproduction of common dandelion.

RESULTS

Here we present an improved Agrobacterium-based genetic transformation and regeneration protocol, a protocol for generation and transformation of protoplasts using free DNA, and a protocol for leaf Agrobacterium infiltration for transient gene expression. These protocols use easily obtainable leaf explants from soil-grown plants and reagents common to most molecular plant laboratories. We show that common markers used in many plant transformation systems function as expected in common dandelion including fluorescent proteins, GUS, and anthocyanin regulation, as well as resistance to kanamycin, Basta, and hygromycin.

CONCLUSION

Reproducible, stable and transient transformation methods are presented that will allow for needed molecular structure and function studies of genes and proteins in T. officinale.

Genetic diversity and evolutionary patterns of Taraxacum kok-saghyz Rodin

Taraxacum kok-saghyz Rodin (TKS) is an important potential alternative source of natural inulin and rubber production, which has great significance for the production of industrial products. In this study, we sequenced 58 wild TKS individuals collected from four different geography regions worldwide to elucidate the population structure, genetic diversity, and the patterns of evolution. Also, the first flowering time, crown diameter, morphological characteristics of leaf, and scape of all TKS individuals were measured and evaluated statistically. Phylogenetic analysis based on SNPs and cluster analysis based on agronomic traits showed that all 58 TKS individuals could be roughly divided into three distinct groups: (a) Zhaosu County in Xinjiang (population AB, including a few individuals from population C and D); (b) Tekes County in Xinjiang (population C); and (c) Tuzkol lake in Kazakhstan (population D). Population D exhibited a closer genetic relationship with population C compared with population AB. Genetic diversity analysis further revealed that population expansion from C and D to AB occurred, as well as gene flow between them. Additionally, some natural selection regions were identified in AB population. Function annotation of candidate genes identified in these regions revealed that they mainly participated in biological regulation processes, such as transporter activity, structural molecule activity, and molecular function regulator. We speculated that the genes identified in selective sweep regions may contribute to TKS adaptation to the Yili River Valley of Xinjiang. In general, this study provides new insights in clarifying population structure and genetic diversity analysis of TKS using SNP molecular markers and agronomic traits.

Natural Polymer Eucommia Ulmoides Rubber: A Novel Material

As the second natural rubber resource, Eucommia ulmoides rubber (EUR) from Eucommia ulmoides Oliver is mainly composed of trans-1,4-polyisoprene, which is the isomer of natural rubber cis-1,4-polyisoprene from Hevea brasiliensis. In the past few years, the great potential application of EUR has received increasing attention, and there is a growing awareness that the natural polymer EUR could become an emerging research topic in field of the novel materials due to its unique and excellent duality of both rubber and plastic. To gain insight into its further development, in this review, the extraction, structure, physicochemical properties, and modification of EUR are discussed in detail. More emphasis on the potential applications in the fields of the environment, agriculture, engineering, and biomedical engineering is summarized. Finally, some insights into the challenges and perspectives of EUR are also suggested.

Microscopic and Transcriptomic Analysis of Pollination Processes in Self-Incompatible Taraxacum koksaghyz

The transition of the Russian dandelion Taraxacum koksaghyz (Asteraceae) to a profitable, alternative crop producing natural rubber and inulin requires the optimization of several agronomic traits, cultivation conditions and harvesting procedures to improve the yield. However, efficient breeding is hindered by the obligatory sexual outcrossing of this species. Several other asters have been investigated to determine the mechanism of self-incompatibility, but the underlying molecular basis remains unclear. We therefore investigated the self-pollination and cross-pollination of two compatible T. koksaghyz varieties (TkMS2 and TkMS3) by microscopy and transcriptomic analysis to shed light on the pollination process. Self-pollination showed typical sporophytic self-incompatibility characteristics, with the rare pollen swelling at the pollen tube apex. In contrast, cross-pollination was characterized by pollen germination and penetration of the stigma by the growing pollen tubes. RNA-Seq was used to profile gene expression in the floret tissue during self-pollination and cross-pollination, and the differentially expressed genes were identified. This revealed three candidates for the early regulation of pollination in T. koksaghyz, which can be used to examine self-incompatibility mechanisms in more detail and to facilitate breeding programs.

SMARTdenovo: a de novo assembler using long noisy reads

Long-read single-molecule sequencing has revolutionized de novo genome assembly and enabled the automated reconstruction of reference-quality genomes. It has also been widely used to study structural variants, phase haplotypes and more. Here, we introduce the assembler SMARTdenovo, a single-molecule sequencing (SMS) assembler that follows the overlap-layout-consensus (OLC) paradigm. SMARTdenovo (RRID: SCR_017622) was designed to be a rapid assembler, which, unlike contemporaneous SMS assemblers, does not require highly accurate raw reads for error correction. It has performed well in the evaluation of congeneric assemblers and has been successfully users for various assembly projects. It is compatible with Canu for assembling high-quality genomes, and several of the assembly strategies in this program have been incorporated into subsequent popular assemblers. The assembler has been in use since 2015; here we provide information on the development of SMARTdenovo and how to implement its algorithms into current projects.

Designing future crops: challenges and strategies for sustainable agriculture

Crop production is facing unprecedented challenges. Despite the fact that the food supply has significantly increased over the past half-century, ~8.9 and 14.3% people are still suffering from hunger and malnutrition, respectively. Agricultural environments are continuously threatened by a booming world population, a shortage of arable land, and rapid changes in climate. To ensure food and ecosystem security, there is a need to design future crops for sustainable agriculture development by maximizing net production and minimalizing undesirable effects on the environment. The future crops design projects, recently launched by the National Natural Science Foundation of China and Chinese Academy of Sciences (CAS), aim to develop a roadmap for rapid design of customized future crops using cutting-edge technologies in the Breeding 4.0 era. In this perspective, we first introduce the background and missions of these projects. We then outline strategies to design future crops, such as improvement of current well-cultivated crops, de novo domestication of wild species and redomestication of current cultivated crops. We further discuss how these ambitious goals can be achieved by the recent development of new integrative omics tools, advanced genome-editing tools and synthetic biology approaches. Finally, we summarize related opportunities and challenges in these projects.

Genome Survey Sequencing of In Vivo Mother Plant and In Vitro Plantlets of Mikania cordata

Mikania cordata, the only native congener of the invasive weed Mikania micrantha in China, is an ideal species for comparative study to reveal the invasion mechanism. However, its genome resources are lagging far behind its congener, which limits the comparative genomic analysis. Our goal is to characterize the genome of M. cordata by next-generation sequencing and propose a scheme for long-read genome sequencing. Previous studies have shown that the genomic resources of the host plant would be affected by the endophytic microbial DNA. An aseptic sample of M. cordata will ensure the proper genome in downstream analysis. Because endophytes are ubiquitous in the greenhouse-grown M. cordata, the in vitro culture with cefotaxime or timentin treatment was undertaken to obtain the aseptic plantlets. The in vivo mother plant and in vitro plantlets were used to survey the genome. The microbial contamination in M. cordata was recognized by blast search and eliminated from the raw reads. The decontaminated sequencing reads were used to predict the genome size, heterozygosity, and repetitive rate. The in vivo plant was so contaminated that microbes occupied substantial sequencing resources and misled the scaffold assembly. Compared with cefotaxime, treatment with timentin performed better in cultivating robust in vitro plantlets. The survey result from the in vitro plantlets was more accurate due to low levels of contamination. The genome size was estimated to be 1.80 Gb with 0.50% heterozygosity and 78.35% repetitive rate. Additionally, 289,831 SSRs were identified in the genome. The genome is heavily contaminated and repetitive; therefore, the in vitro culture technique and long-read sequencing technology are recommended to generate a high-quality and highly contiguous genome.

Comparison of Morphological Characteristics and Determination of Different Patterns for Rubber Particles in Dandelion and Different Rubber Grass Varieties

Russian dandelion Taraxacum kok-saghyz (TKS) is one promising alternative crop for natural rubber production. However, it is easily confused with other dandelions. In this study, we performed a systematical comparison of the morphological characteristics for different TKS varieties and common dandelion Taraxacum officinale (TO). Our results demonstrated that several obvious differences in morphology can be found between TKS and TO. TO leaf is a pinnate shape, its margin is heavily jagged and its base is cuneate, but TKS leaf is more cuneate and its leaf margin is nearly smooth and round. There are obvious differences for the outer bracts of TO and TKS flower buds. TKS bracts are oblanceolate, apex obtuse, margin smooth and sinuate, and its outer layer of flower buds and faceplate involucre sepal is buckled inward to form a certain angle. TKS is self-incompatible, and its seeds are spindle-shaped achene and show upright plumpness. A large amount of laticifer cells and rubber particles can be detected from many TKS tissues, and dry roots of TKS contain high contents of natural rubber. Laticifer cells and rubber particles can only be examined in the vein, stem, and roots of TKS. Our statical results also revealed that the numbers of laticifer cells and rubber particles have a positive relationship with the rubber content in TKS roots. These morphological features can help us to easily distinguish TKS from common dandelion and approximately estimate the rubber content in the roots of different TKS varieties for TKS breeding in future.

High-quality de novo assembly of the Eucommia ulmoides haploid genome provides new insights into evolution and rubber biosynthesis

We report the acquisition of a high-quality haploid chromosome-scale genome assembly for the first time in a tree species, Eucommia ulmoides, which is known for its rubber biosynthesis and medicinal applications. The assembly was obtained by applying PacBio and Hi-C technologies to a haploid that we specifically generated. Compared to the initial genome release, this one has significantly improved assembly quality. The scaffold N50 (53.15 MB) increased 28-fold, and the repetitive sequence content (520 Mb) increased by 158.24 Mb, whereas the number of gaps decreased from 104,772 to 128. A total of 92.87% of the 26,001 predicted protein-coding genes identified with multiple strategies were anchored to the 17 chromosomes. A new whole-genome duplication event was superimposed on the earlier γ paleohexaploidization event, and the expansion of long terminal repeats contributed greatly to the evolution of the genome. The more primitive rubber biosynthesis of this species, as opposed to that in Hevea brasiliensis, relies on the methylerythritol-phosphate pathway rather than the mevalonate pathway to synthesize isoprenyl diphosphate, as the MEP pathway operates predominantly in trans-polyisoprene-containing leaves and central peels. Chlorogenic acid biosynthesis pathway enzymes were preferentially expressed in leaves rather than in bark. This assembly with higher sequence contiguity can foster not only studies on genome structure and evolution, gene mapping, epigenetic analysis and functional genomics but also efforts to improve E. ulmoides for industrial and medical uses through genetic engineering.

Combinatorial Metabolic Engineering in Saccharomyces cerevisiae for the Enhanced Production of the FPP-Derived Sesquiterpene Germacrene

Farnesyl diphosphate (FPP)-derived isoprenoids represent a diverse group of plant secondary metabolites with great economic potential. To enable their efficient production in the heterologous host Saccharomyces cerevisiae, we refined a metabolic engineering strategy using the CRISPR/Cas9 system with the aim of increasing the availability of FPP for downstream reactions. The strategy included the overexpression of mevalonate pathway (MVA) genes, the redirection of metabolic flux towards desired product formation and the knockout of genes responsible for competitive reactions. Following the optimisation of culture conditions, the availability of the improved FPP biosynthesis for downstream reactions was demonstrated by the expression of a germacrene synthase from dandelion. Subsequently, biosynthesis of significant amounts of germacrene-A was observed in the most productive strain compared to the wild type. Thus, the presented strategy is an excellent tool to increase FPP-derived isoprenoid biosynthesis in yeast.

The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation

Tea plant is an important economic crop, which is used to produce the world's oldest and most widely consumed tea beverages. Here, we present a high-quality reference genome assembly of the tea plant (Camellia sinensis var. sinensis) consisting of 15 pseudo-chromosomes. LTR retrotransposons (LTR-RTs) account for 70.38% of the genome, and we present evidence that LTR-RTs play critical roles in genome size expansion and the transcriptional diversification of tea plant genes through preferential insertion in promoter regions and introns. Genes, particularly those coding for terpene biosynthesis proteins, associated with tea aroma and stress resistance were significantly amplified through recent tandem duplications and exist as gene clusters in tea plant genome. Phylogenetic analysis of the sequences of 81 tea plant accessions with diverse origins revealed three well-differentiated tea plant populations, supporting the proposition for the southwest origin of the Chinese cultivated tea plant and its later spread to western Asia through introduction. Domestication and modern breeding left significant signatures on hundreds of genes in the tea plant genome, particularly those associated with tea quality and stress resistance. The genomic sequences of the reported reference and resequenced tea plant accessions provide valuable resources for future functional genomics study and molecular breeding of improved cultivars of tea plants.

Molecular adaptation to salinity fluctuation in tropical intertidal environments of a mangrove tree Sonneratia alba

BACKGROUND

Mangroves have adapted to intertidal zones - the interface between terrestrial and marine ecosystems. Various studies have shown adaptive evolution in mangroves at physiological, ecological, and genomic levels. However, these studies paid little attention to gene regulation of salt adaptation by transcriptome profiles.

RESULTS

We sequenced the transcriptomes of Sonneratia alba under low (fresh water), medium (half the seawater salinity), and high salt (seawater salinity) conditions and investigated the underlying transcriptional regulation of salt adaptation. In leaf tissue, 64% potential salinity-related genes were not differentially expressed when salinity increased from freshwater to medium levels, but became up- or down-regulated when salt concentrations further increased to levels found in sea water, indicating that these genes are well adapted to the medium saline condition. We inferred that both maintenance and regulation of cellular environmental homeostasis are important adaptive processes in S. alba. i) The sulfur metabolism as well as flavone and flavonol biosynthesis KEGG pathways were significantly enriched among up-regulated genes in leaves. They are both involved in scavenging ROS or synthesis and accumulation of osmosis-related metabolites in plants. ii) There was a significantly increased percentage of transcription factor-encoding genes among up-regulated transcripts. High expressions of salt tolerance-related TF families were found under high salt conditions. iii) Some genes up-regulated in response to salt treatment showed signs of adaptive evolution at the amino acid level and might contribute to adaptation to fluctuating intertidal environments.

CONCLUSIONS

This study first elucidates the mechanism of high-salt adaptation in mangroves at the whole-transcriptome level by salt gradient experimental treatments. It reveals that several candidate genes (including salt-related genes, TF-encoding genes, and PSGs) and major pathways are involved in adaptation to high-salt environments. Our study also provides a valuable resource for future investigation of adaptive evolution in extreme environments.

Molecular Evolution in Small Steps under Prevailing Negative Selection: A Nearly Universal Rule of Codon Substitution

The widely accepted view that evolution proceeds in small steps is based on two premises: 1) negative selection acts strongly against large differences and 2) positive selection favors small-step changes. The two premises are not biologically connected and should be evaluated separately. We now extend a previous approach to studying codon evolution in the entire genome. Codon substitution rate is a function of the physicochemical distance between amino acids (AAs), equated with the step size of evolution. Between nine pairs of closely related species of plants, invertebrates, and vertebrates, the evolutionary rate is strongly and negatively correlated with a set of AA distances (Δ_U, scaled to [0, 1]). Δ_U, a composite measure of evolutionary rates across diverse taxa, is influenced by almost all of the 48 physicochemical properties used here. The new analyses reveal a crucial trend hidden from previous studies: Δ_U is strongly correlated with the evolutionary rate (R² > 0.8) only when the genes are predominantly under negative selection. Because most genes in most taxa are strongly constrained by negative selection, Δ_U has indeed appeared to be a nearly universal measure of codon evolution. In conclusion, molecular evolution at the codon level generally takes small steps due to the prevailing negative selection. Whether positive selection may, or may not, follow the small-step rule is addressed in a companion study.

Comparative proteome and metabolome analyses of latex-exuding and non-exuding Taraxacum koksaghyz roots provide insights into laticifer biology

Taraxacum koksaghyz has been identified as one of the most promising alternative rubber crops. Its high-quality rubber is produced in the latex of laticifers, a specialized cell type that is organized in a network of elongated tubules throughout the entire plant body. In order to gain insights into the physiological role(s) of latex and hence laticifer biology, we examine the effects of barnase-induced latex RNA degradation on the metabolite and protein compositions in the roots. We established high-quality datasets that enabled precise discrimination between cellular and physiological processes in laticifers and non-laticifer cell types of roots at different vegetative stages. We identified numerous latex-specific proteins, including a perilipin-like protein that has not been studied in plants yet. The barnase-expressing plants revealed a phenotype that did not exude latex, which may provide a valuable genetic basis for future studies of plant-environment interactions concerning latex and also help to clarify the evolution and arbitrary distribution of latex throughout the plant kingdom. The overview of temporal changes in composition and protein abundance provided by our data opens the way for a deeper understanding of the molecular interactions, reactions, and network relationships that underlie the different metabolic pathways in the roots of this potential rubber crop.

The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

Evolutionary Analysis of Calcium-Dependent Protein Kinase in Five Asteraceae Species

Calcium-dependent protein kinase (CPK) is crucial in Ca²⁺ signal transduction, and is a large gene family in plants. In our previous work, we reported Hevea brasiliensis CPKs were important for natural rubber biosynthesis. However, this CPK gene family in other rubber producing plants has not been investigated. Here, we report the CPKs in five representative Asteraceae species, including three rubber-producing and two non-rubber species. A total of 34, 34, 40, 34 and 30 CPKs were identified from Taraxacum koksaghyz, Lactuca sativa, Helianthus annuus, Chrysanthemum nankingense and Cynara cardunculus, respectively. All CPKs were classified into four individual groups (group I to IV). In addition, 10 TkCPKs, 11 LsCPKs, 20 HaCPKs, 13 CnCPKs and 7 CcCPKs duplicated paralogs were identified. Further evolutionary analysis showed that, compared to other subfamilies, the group III had been expanded in the Asteraceae species, especially in the rubber-producing species. Meanwhile, the CPKs in group III from Asteraceae species tend to expand with low calcium binding capacity. This study provides a systematical evolutionary investigation of the CPKs in five representative Asteraceae species, suggesting that the sub-family specific expansion of CPKs might be related to natural rubber producing.

Uncovering mechanisms of rubber biosynthesis in Taraxacum koksaghyz - role of cis-prenyltransferase-like 1 protein

The Russian dandelion Taraxacum koksaghyz synthesizes considerable amounts of high-molecular-weight rubber in its roots. The characterization of factors that participate in natural rubber biosynthesis is fundamental for the establishment of T. koksaghyz as a rubber crop. The cis-1,4-isoprene polymers are stored in rubber particles. Located at the particle surface, the rubber transferase complex, member of the cis-prenyltransferase (cisPT) enzyme family, catalyzes the elongation of the rubber chains. An active rubber transferase heteromer requires a cisPT subunit (CPT) as well as a CPT-like subunit (CPTL), of which T. koksaghyz has two homologous forms: TkCPTL1 and TkCPTL2, which potentially associate with the rubber transferase complex. Knockdown of TkCPTL1, which is predominantly expressed in latex, led to abolished poly(cis-1,4-isoprene) synthesis but unaffected dolichol content, whereas levels of triterpenes and inulin were elevated in roots. Analyses of latex from these TkCPTL1-RNAi plants revealed particles that were similar to native rubber particles regarding their particle size, phospholipid composition, and presence of small rubber particle proteins (SRPPs). We found that the particles encapsulated triterpenes in a phospholipid shell stabilized by SRPPs. Conversely, downregulating the low-expressed TkCPTL2 showed no altered phenotype, suggesting its protein function is redundant in T. koksaghyz. MS-based comparison of latex proteomes from TkCPTL1-RNAi plants and T. koksaghyz wild-types discovered putative factors that convert metabolites in biosynthetic pathways connected to isoprenoids or that synthesize components of the rubber particle shell.

Natural rubber biosynthesis in plants, the rubber transferase complex, and metabolic engineering progress and prospects

Natural rubber (NR) is a nonfungible and valuable biopolymer, used to manufacture ~50 000 rubber products, including tires and medical gloves. Current production of NR is derived entirely from the para rubber tree (Hevea brasiliensis). The increasing demand for NR, coupled with limitations and vulnerability of H. brasiliensis production systems, has induced increasing interest among scientists and companies in potential alternative NR crops. Genetic/metabolic pathway engineering approaches, to generate NR-enriched genotypes of alternative NR plants, are of great importance. However, although our knowledge of rubber biochemistry has significantly advanced, our current understanding of NR biosynthesis, the biosynthetic machinery and the molecular mechanisms involved remains incomplete. Two spatially separated metabolic pathways provide precursors for NR biosynthesis in plants and their genes and enzymes/complexes are quite well understood. In contrast, understanding of the proteins and genes involved in the final step(s)-the synthesis of the high molecular weight rubber polymer itself-is only now beginning to emerge. In this review, we provide a critical evaluation of recent research developments in NR biosynthesis, in vitro reconstitution, and the genetic and metabolic pathway engineering advances intended to improve NR content in plants, including H. brasiliensis, two other prospective alternative rubber crops, namely the rubber dandelion and guayule, and model species, such as lettuce. We describe a new model of the rubber transferase complex, which integrates these developments. In addition, we highlight the current challenges in NR biosynthesis research and future perspectives on metabolic pathway engineering of NR to speed alternative rubber crop commercial development.

Integration of Transcriptomes, Small RNAs, and Degradome Sequencing to Identify Putative miRNAs and Their Targets Related to Eu-Rubber Biosynthesis in Eucommia ulmoides

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

Proteomic Landscape of the Mature Roots in a Rubber-Producing Grass Taraxacum Kok-saghyz

The rubber grass Taraxacum kok-saghyz (TKS) contains large amounts of natural rubber (cis-1,4-polyisoprene) in its enlarged roots and it is an alternative crop source of natural rubber. Natural rubber biosynthesis (NRB) and storage in the mature roots of TKS is a cascade process involving many genes, proteins and their cofactors. The TKS genome has just been annotated and many NRB-related genes have been determined. However, there is limited knowledge about the protein regulation mechanism for NRB in TKS roots. We identified 371 protein species from the mature roots of TKS by combining two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). Meanwhile, a large-scale shotgun analysis of proteins in TKS roots at the enlargement stage was performed, and 3545 individual proteins were determined. Subsequently, all identified proteins from 2-DE gel and shotgun MS in TKS roots were subject to gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and most proteins were involved in carbon metabolic process with catalytic activity in membrane-bounded organelles, followed by proteins with binding ability, transportation and phenylpropanoid biosynthesis activities. Fifty-eight NRB-related proteins, including eight small rubber particle protein (SRPP) and two rubber elongation factor(REF) members, were identified from the TKS roots, and these proteins were involved in both mevalonate acid (MVA) and methylerythritol phosphate (MEP) pathways. To our best knowledge, it is the first high-resolution draft proteome map of the mature TKS roots. Our proteomics of TKS roots revealed both MVA and MEP pathways are important for NRB, and SRPP might be more important than REF for NRB in TKS roots. These findings would not only deepen our understanding of the TKS root proteome, but also provide new evidence on the roles of these NRB-related proteins in the mature TKS roots.

Loss of function mutation of the Rapid Alkalinization Factor (RALF1)-like peptide in the dandelion Taraxacum koksaghyz entails a high-biomass taproot phenotype

The Russian dandelion (Taraxacum koksaghyz) is a promising source of inulin and natural rubber because large amounts of both feedstocks can be extracted from its roots. However, the domestication of T. koksaghyz requires the development of stable agronomic traits such as higher yields of inulin and natural rubber, a higher root biomass, and an agronomically preferable root morphology which is more suitable for cultivation and harvesting. Arabidopsis thaliana Rapid Alkalinisation Factor 1 (RALF1) has been shown to suppress root growth. We identified the T. koksaghyz orthologue TkRALF-like 1 and knocked out the corresponding gene (TkRALFL1) using the CRISPR/Cas9 system to determine its impact on root morphology, biomass, and inulin and natural rubber yields. The TkRALFL1 knockout lines more frequently developed a taproot phenotype which is easier to cultivate and harvest, as well as a higher root biomass and greater yields of both inulin and natural rubber. The TkRALFL1 gene could therefore be suitable as a genetic marker to support the breeding of profitable new dandelion varieties with improved agronomic traits. To our knowledge, this is the first study addressing the root system of T. koksaghyz to enhance the agronomic performance.

The enzymes OSC1 and CYP716A263 produce a high variety of triterpenoids in the latex of Taraxacum koksaghyz

Only very little is known about the resin composition of natural rubber from the dandelion species Taraxacum koksaghyz, thus its full characterization could provide new insights into how the isoprenoid end-products influence the physical properties of natural rubber, and this resin might be a good source of highly diverse triterpenoids. Here, we present a comprehensive analysis of the triterpenoid composition in an acetone extract and identified 13 triterpenes and triterpenoids also including the so far unknown pentacyclic compounds lup-19(21)-en-3-ol (1) and its ketone lup-19(21)-en-3-one (2). We purified single triterpenes from the acetone extract by developing a two-step HPLC system that is adapted to the structural differences of the described triterpenoids. Furthermore, we isolated six different oxidosqualene cyclases (OSCs) and two P450 enzymes, and we functionally characterized TkOSC1 and CYP716A263 in Nicotiana benthamiana and Saccharomyces cerevisiae in detail. TkOSC1 is a multifunctional OSC that was capable of synthesizing at least four of the latex-predominant pentacyclic triterpenes (taraxasterol, α-, β-amyrin and lup-19(21)-en-3-ol) while CYP716A263 oxidized pentacyclic triterpenes at the C-3 position. The identified enzymes responsible for biosynthesis and modification of pentacyclic triterpenes in T. koksaghyz latex may represent excellent tools for bioengineering approaches to produce pentacyclic triterpenes heterologously.

Taraxacum kok-saghyz (rubber dandelion) genomic microsatellite loci reveal modest genetic diversity and cross-amplify broadly to related species

Taraxacum kok-saghyz (TKS) carries great potential as alternative natural rubber source. To better inform future breeding efforts with TKS and gain a deeper understanding of its genetic diversity, we utilized de novo sequencing to generate novel genomic simple sequence repeats markers (gSSRs). We utilized 25 gSSRs on a collection of genomic DNA (gDNA) samples from germplasm bank, and two gDNA samples from historical herbarium specimens. PCR coupled with capillary electrophoresis and an array of population genetics tools were employed to analyze the dataset of our study as well as a dataset of the recently published genic SSRs (eSSRs) generated on the same germplasm. Our results using both gSSRs and eSSRs revealed that TKS has low- to- moderate genetic diversity with most of it partitioned to the individuals and individuals within populations, whereas the species lacked population structure. Nineteen of the 25 gSSR markers cross-amplified to other Taraxacum spp. collected from Southeastern United States and identified as T. officinale by ITS sequencing. We used a subset of 14 gSSRs to estimate the genetic diversity of the T. officinale gDNA collection. In contrast to the obligatory outcrossing TKS, T. officinale presented evidence for population structure and clonal reproduction, which agreed with the species biology. We mapped the molecular markers sequences from this study and several others to the well-annotated sunflower genome. Our gSSRs present a functional tool for the biodiversity analyses in Taraxacum, but also in the related genera, as well as in the closely related tribes of the Asteraceae.

Structural variation during dog domestication: insights from gray wolf and dhole genomes

Several processes like phenotypic evolution, disease susceptibility and environmental adaptations, which fashion the domestication of animals, are largely attributable to structural variations (SVs) in the genome. Here, we present high-quality draft genomes of the gray wolf (Canis lupus) and dhole (Cuon alpinus) with scaffold N50 of 6.04 Mb and 3.96 Mb, respectively. Sequence alignment comprising genomes of three canid species reveals SVs specific to the dog, particularly 16 315 insertions, 2565 deletions, 443 repeats, 16 inversions and 15 translocations. Functional annotation of the dog SVs associated with genes indicates their enrichments in energy metabolisms, neurological processes and immune systems. Interestingly, we identify and verify at population level an insertion fully covering a copy of the AKR1B1 (Aldo-Keto Reductase Family 1 Member B) transcript. Transcriptome analysis reveals a high level of expression of the new AKR1B1 copy in the small intestine and liver, implying an increase in de novo fatty acid synthesis and antioxidant ability in dog compared to gray wolf, likely in response to dietary shifts during the agricultural revolution. For the first time, we report a comprehensive analysis of the evolutionary dynamics of SVs during the domestication step of dogs. Our findings demonstrate that retroposition can birth new genes to facilitate domestication, and affirm the importance of large-scale genomic variants in domestication studies.