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 F1 mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F1 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 F1 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.

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.

Environmental and genealogical effects on DNA methylation in a widespread apomictic dandelion lineage.

DNA methylation in plant genomes occurs in different sequences and genomic contexts that have very different properties. DNA methylation that occurs in CG (mCG) sequence context shows transgenerational stability and high epimutation rate, and can thus provide genealogical information at short time scales. However, due to meta-stability and because mCG variants may arise due to other factors than epimutation, such as environmental stress exposure, it is not clear how well mCG captures genealogical information at micro-evolutionary time scales. Here, we analysed DNA methylation variation between accessions from a geographically widespread, apomictic common dandelion (Taraxacum officinale) lineage when grown experimentally under different light conditions. Using a reduced-representation bisulphite sequencing approach, we show that the light treatment induced differentially methylated cytosines (DMCs) in all sequence contexts, with a bias towards transposable elements. Accession differences were associated mainly with DMCs in CG context. Hierarchical clustering of samples based on total mCG profiles revealed a perfect clustering of samples by accession identity, irrespective of light conditions. Using microsatellite information as a benchmark of genetic divergence within the clonal lineage, we show that genetic divergence between accessions correlates strongly with overall mCG profiles. However, our results suggest that environmental effects that do occur in CG context may produce a heritable signal that partly dilutes the genealogical signal. Our study shows that methylation information in plants can be used to reconstruct micro-evolutionary genealogy, providing a useful tool in systems that lack genetic variation such as clonal and vegetatively propagated plants.

Characterisation of the complete mitochondrial genome of Taraxacum mongolicum revealed five repeat-mediated recombinations.

KEY MESSAGE: We reported the mitochondrial genome of Taraxacum mongolicum for the first time. Five pairs of repeats that can mediate recombination were validated, leading to multiple conformations of genome. Taraxacum mongolicum belongs to the Asteraceae family and has important pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled the complete mitochondrial genome (mitogenome) of T. mongolicum using Illumina and Oxford Nanopore sequencing data. This genome corresponded to a circular molecule 304,467 bp long. It encodes 52 unique genes including 31 protein-coding, 3 ribosomal RNA (rRNA) and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, the existence of alternative conformations mediated by five repetitive sequences in the mitogenome was identified and validated. Recombination mediated by the inverted repeats resulted in two conformations. Conversely, recombination mediated by the two direct repeats broke one large circular molecule into two subgenomic circular molecules. Furthermore, we identified 12 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. Lastly, we identified a total of 278 RNA-editing sites in protein-coding sequences based on RNA-seq data. Among them, cox1 and nad5 gene has the most sites (21), followed by the nad2 gene with 19 sites. We successfully validated 213 predicted RNA-editing sites using PCR amplification and Sanger sequencing. This project reported the first mitogenome of T. mongolicum and demonstrated its multiple conformations generated by repeat-mediated recombination. This genome could provide critical information for the molecular breeding of T. mongolicum, and also be used as a reference genome for other species of the genus Taraxacum.

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.

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.

Metabolic and Transcriptomic Analyses Reveal the Effects of Ethephon on Taraxacum kok-saghyz Rodin.

The roots of Taraxacum kok-saghyz Rodin (TKS) are well-known and valued for their rubber-producing ability. Therefore, research on the analysis and detection of metabolites from the roots of TKS have been reported in previous studies. However, all of these studies have the shortcoming of focusing on only the rubber of TKS, without profiling the other metabolites in a systematic and comprehensive way. Here, the primary and secondary metabolites from the leaves of TKS were investigated using UPLC-ESI-MS/MS, and a total of 229 metabolites were characterized. Carboxylic acid derivatives, fatty acyls, phenols, and organooxygen compounds were found to be the major metabolites of TKS. The transcriptome data indicated that ribosomal, glycolysis/gluconeogenesis, phenylpropanoid biosynthesis, and linoleic acid metabolism genes were significantly differentially expressed. This study is the first to report the differences in the metabolic and transcriptome profiles of TKS leaves under exogenous ethephon spray, which improves our understanding of the main metabolites and their molecular mechanisms in TKS leaves.

Niche divergence contributes to geographical parthenogenesis in two dandelion taxa.

Many sexual-asexual complexes show a distinct pattern where the asexuals have larger and more northerly ranges than closely related sexuals. A prime candidate to explain this so-called "geographical parthenogenesis" is ecological niche divergence between the sexuals and asexuals. Modern niche modelling techniques allow testing niche divergence by directly comparing the niches of sexuals and asexuals. In this study, I use such techniques to perform range-wide tests of whether nine bioclimatic variables, including annual mean temperature and annual precipitation, contribute to geographical parthenogenesis in two dandelion taxa: Taraxacum section Ruderalia and Taraxacum section Erythrosperma, which are both comprised of sexual diploids and asexual triploids. For both sections, I found evidence of niche divergence, though the exact nature of this divergence was different for the two sections. In section Ruderalia, the sexuals preferred warmer and wetter conditions, whereas in section Erythrosperma, the sexuals preferred dryer conditions. Using Species Distribution Modelling, consistent differences between the sexuals and asexuals were found when looking at the niche determinants: the variables that are most important for modelling the distribution. Furthermore, and in contrast with theoretical expectations that predict that the sexuals should have a wider niche, in section Erythrosperma the asexuals were found to have a wider niche than the sexuals. In conclusion, differences in niche optima, niche determinants, and niche width all contribute to the pattern of geographical parthenogenesis of these two dandelion taxa. However, the results also indicate that the exact causation of geographical parthenogenesis is not uniform across taxa.

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.

Small rubber particle proteins from Taraxacum brevicorniculatum promote stress tolerance and influence the size and distribution of lipid droplets and artificial poly(cis-1,4-isoprene) bodies.

Natural rubber biosynthesis occurs on rubber particles, i.e. organelles resembling small lipid droplets localized in the laticifers of latex-containing plant species, such as Hevea brasiliensis and Taraxacum brevicorniculatum. The latter expresses five small rubber particle protein (SRPP) isoforms named TbSRPP1-5, the most abundant proteins in rubber particles. These proteins maintain particle stability and are therefore necessary for rubber biosynthesis. TbSRPP1-5 were transiently expressed in Nicotiana benthamiana protoplasts and the proteins were found to be localized on lipid droplets and in the endoplasmic reticulum, with TbSRPP1 and TbSRPP3 also present in the cytosol. Bimolecular fluorescence complementation confirmed pairwise interactions between all proteins except TbSRPP2. The corresponding genes showed diverse expression profiles in young T. brevicorniculatum plants exposed to abiotic stress, and all except TbSRPP4 and TbSRPP5 were upregulated. Young Arabidopsis thaliana plants that overexpressed TbSRPP2 and TbSRPP3 tolerated drought stress better than wild-type plants. Furthermore, we used rubber particle extracts and standards to investigate the affinity of the TbSRPPs for different phospholipids, revealing a preference for negatively charged head groups and 18:2/16:0 fatty acid chains. This finding may explain the effect of TbSRPP3-5 on the dispersity of artificial poly(cis-1,4-isoprene) bodies and on the lipid droplet distribution we observed in N. benthamiana leaves. Our data provide insight into the assembly of TbSRPPs on rubber particles, their role in rubber particle structure, and the link between rubber biosynthesis and lipid droplet-associated stress responses, suggesting that SRPPs form the basis of evolutionarily conserved intracellular complexes in plants.

Cloning and Characterization of Oxidosqualene Cyclases Involved in Taraxasterol, Taraxerol and Bauerenol Triterpene Biosynthesis in Taraxacum coreanum.

Triterpenes, consisting of six isoprene units, are one of the largest classes of natural compounds in plants. The genus Taraxacum is in the family Asteraceae and is widely distributed in the Northern Hemisphere. Various triterpenes, especially taraxerol and taraxasterol, are present in Taraxacum plants. Triterpene biosynthesis occurs through the action of oxidosqualene cyclase (OSC), which generates various types of triterpenes from 2,3-oxidosqualene after the rearrangement of the triterpene skeleton. However, no functional characterization of the OSC genes involved in triterpene biosynthesis, except for a lupeol synthase in Taraxacum officinale, has been performed. Taraxacum coreanum, or Korean dandelion, grows in Korea and China. Putative OSC genes in T. coreanum plants were isolated by transcriptome analysis, and four of these (TcOSC1, TcOSC2, TcOSC3 and TcOSC4) were functionally characterized by heterologous expression in yeast. Both TcOSC1 and TcOSC2 were closely related to dammarenediol-II synthases. TcOSC3 and TcOSC4 were strongly grouped with ß-amyrin synthases. Functional analysis revealed that TcOSC1 produced several triterpenes, including taraxasterol; Ψ-taraxasterol; α-, ß- and δ-amyrin; and dammarenediol-II. TcOSC2 catalyzed the production of bauerenol and another unknown triterpene, TcOSC3 catalyzed the production of ß-amyrin. TcOSC4 catalyzed the production of taraxerol. Moreover, we identified taraxasterol, ψ-taraxasterol, taraxerol, lupeol, δ-amyrin, α-amyrin, ß-amyrin and bauerenol in the roots and leaves of T. coreanum. Our results suggest that TcOSC1, TcOSC2, TcOSC3 and TcOSC4 are key triterpene biosynthetic enzymes in T. coreanum. These enzymes are novel triterpene synthases involved in the production of taraxasterol, bauerenol and taraxerol.

A Latex Metabolite Benefits Plant Fitness under Root Herbivore Attack.

Plants produce large amounts of secondary metabolites in their shoots and roots and store them in specialized secretory structures. Although secondary metabolites and their secretory structures are commonly assumed to have a defensive function, evidence that they benefit plant fitness under herbivore attack is scarce, especially below ground. Here, we tested whether latex secondary metabolites produced by the common dandelion (Taraxacum officinale agg.) decrease the performance of its major native insect root herbivore, the larvae of the common cockchafer (Melolontha melolontha), and benefit plant vegetative and reproductive fitness under M. melolontha attack. Across 17 T. officinale genotypes screened by gas and liquid chromatography, latex concentrations of the sesquiterpene lactone taraxinic acid ß-D-glucopyranosyl ester (TA-G) were negatively associated with M. melolontha larval growth. Adding purified TA-G to artificial diet at ecologically relevant concentrations reduced larval feeding. Silencing the germacrene A synthase ToGAS1, an enzyme that was identified to catalyze the first committed step of TA-G biosynthesis, resulted in a 90% reduction of TA-G levels and a pronounced increase in M. melolontha feeding. Transgenic, TA-G-deficient lines were preferred by M. melolontha and suffered three times more root biomass reduction than control lines. In a common garden experiment involving over 2,000 T. officinale individuals belonging to 17 different genotypes, high TA-G concentrations were associated with the maintenance of high vegetative and reproductive fitness under M. melolontha attack. Taken together, our study demonstrates that a latex secondary metabolite benefits plants under herbivore attack, a result that provides a mechanistic framework for root herbivore driven natural selection and evolution of plant defenses below ground.

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.

Oxidosqualene Cyclase Knock-Down in Latex of Taraxacum koksaghyz Reduces Triterpenes in Roots and Separated Natural Rubber.

In addition to natural rubber (NR), several triterpenes are synthesized in laticifers of the Russian dandelion (Taraxacum koksaghyz). Detailed analysis of NR and resin contents revealed different concentrations of various pentacyclic triterpenes such as α-, ß-amyrin and taraxasterol, which strongly affect the mechanical properties of the resulting rubber material. Therefore, the reduction of triterpene content would certainly improve the industrial applications of dandelion NR. We developed T. koksaghyz plants with reduced triterpene contents by tissue-specific downregulation of major laticifer-specific oxidosqualene cyclases (OSCs) by RNA interference, resulting in an almost 67% reduction in the triterpene content of NR. Plants of the T1 generation showed no obvious phenotypic changes and the rubber yield also remained unaffected. Hence, this study will provide a solid basis for subsequent modern breeding programs to develop Russian dandelion plants with low and stable triterpene levels.

Comparative transcriptomics between high and low rubber producing Taraxacum kok-saghyz R. plants.

BACKGROUND: Taraxacum kok-saghyz R. (Tks) is a promising alternative species to Hevea brasiliensis for production of high quality natural rubber (NR). A comparative transcriptome analysis of plants with differential production of NR will contribute to elucidate which genes are involved in the synthesis, regulation and accumulation of this natural polymer and could help to develop Tks into a rubber crop. RESULTS: We measured rubber content in the latex of 90 individual Tks plants from 9 accessions, observing a high degree of variability. We carried out de novo root transcriptome sequencing, assembly, annotation and comparison of gene expression of plants with the lower (LR plants) and the higher rubber content (HR plants). The transcriptome analysis also included one plant that did not expel latex, in principle depleted of latex transcripts. Moreover, the transcription of some genes well known to play a major role in rubber biosynthesis, was probed by qRT-PCR. Our analysis showed a high modulation of genes involved in the synthesis of NR between LR and HR plants, and evidenced that genes involved in sesquiterpenoids, monoterpenoids and phenylpropanoid biosynthesis are upregulated in LR plants. CONCLUSIONS: Our results show that a higher amount of rubber in the latex in HR plants is positively correlated with high expression levels of a number of genes directly involved in rubber synthesis showing that NR production is highly controlled at transcriptional level. On the other hand, lower amounts of rubber in LR plants is related with higher expression of genes involved in the synthesis of other secondary metabolites that, we hypothesize, may compete towards NR biosynthesis. This dataset represents a fundamental genomic resource for the study of Tks and the comprehension of the synthesis of NR and other biochemically and pharmacologically relevant compounds in the Taraxacum genus.

A specific amino acid residue in the catalytic site of dandelion polyphenol oxidases acts as 'selector' for substrate specificity.

KEY MESSAGE: Successful site-directed mutagenesis combined with in silico modeling and docking studies for the first time offers experimental proof of the role of the 'substrate selector' residue in plant polyphenol oxidases. The plant and fungi enzymes responsible for tissue browning are called polyphenol oxidases (PPOs). In plants, PPOs often occur as families of isoenzymes which are differentially expressed, but little is known about their physiological roles or natural substrates. In a recent study that explored these structure-function relationships, the eleven known dandelion (Taraxacum officinale) PPOs were shown to separate into two different phylogenetic groups differing in catalytic cavity architecture, kinetic parameters, and substrate range. The same study proposed that the PPOs' substrate specificity is controlled by one specific amino acid residue positioned at the entrance to the catalytic site: whereas group 1 dandelion PPOs possess a hydrophobic isoleucine (I) at position HB2+1, group 2 PPOs exhibit a larger, positively charged arginine (R). However, this suggestion was only based on bioinformatic analyses, not experiments. To experimentally investigate this hypothesis, we converted group 1 ToPPO-2 and group 2 ToPPO-6 into PPO-2-I244R and PPO-6-R254I, respectively, and expressed them in E. coli. By performing detailed kinetic characterization and in silico docking studies, we found that replacing this single amino acid significantly changed the PPO's substrate specificity. Our findings therefore proof the role of the 'substrate selector' in plant PPOs.

Small RNAs Reflect Grandparental Environments in Apomictic Dandelion.

Plants can show long-term effects of environmental stresses and in some cases a stress "memory" has been reported to persist across generations, potentially mediated by epigenetic mechanisms. However, few documented cases exist of transgenerational effects that persist for multiple generations and it remains unclear if or how epigenetic mechanisms are involved. Here, we show that the composition of small regulatory RNAs in apomictic dandelion lineages reveals a footprint of drought stress and salicylic acid treatment experienced two generations ago. Overall proportions of 21 and 24 nt RNA pools were shifted due to grandparental treatments. While individual genes did not show strong up- or downregulation of associated sRNAs, the subset of genes that showed the strongest shifts in sRNA abundance was significantly enriched for several GO terms including stress-specific functions. This suggests that a stress-induced signal was transmitted across multiple unexposed generations leading to persistent changes in epigenetic gene regulation.

Intergenerational environmental effects: functional signals in offspring transcriptomes and metabolomes after parental jasmonic acid treatment in apomictic dandelion.

Parental environments can influence offspring traits. However, the magnitude of the impact of parental environments on offspring molecular phenotypes is poorly understood. Here, we test the direct effects and intergenerational effects of jasmonic acid (JA) treatment, which is involved in herbivory-induced defense signaling, on transcriptomes and metabolomes in apomictic common dandelion (Taraxacum officinale). In a full factorial crossed design with parental and offspring JA and control treatments, we performed leaf RNA-seq gene expression analysis, LC-MS metabolomics and total phenolics assays in offspring plants. Expression analysis, leveraged by a de novo assembled transcriptome, revealed an induced response to JA exposure that is consistent with known JA effects. The intergenerational effect of treatment was considerable: 307 of 858 detected JA-responsive transcripts were affected by parental JA treatment. In terms of the numbers of metabolites affected, the magnitude of the chemical response to parental JA exposure was c. 10% of the direct JA treatment response. Transcriptome and metabolome analyses both identified the phosphatidylinositol signaling pathway as a target of intergenerational JA effects. Our results highlight that parental environments can have substantial effects in offspring generations. Transcriptome and metabolome assays provide a basis for zooming in on the potential mechanisms of inherited JA effects.

Chloroplast genome resources and molecular markers differentiate rubber dandelion species from weedy relatives.

BACKGROUND: Rubber dandelion (Taraxacum kok-saghyz, TK) is being developed as a domestic source of natural rubber to meet increasing global demand. However, the domestication of TK is complicated by its colocation with two weedy dandelion species, Taraxacum brevicorniculatum (TB) and the common dandelion (Taraxacum officinale, TO). TB is often present as a seed contaminant within TK accessions, while TO is a pandemic weed, which may have the potential to hybridize with TK. To discriminate these species at the molecular level, and facilitate gene flow studies between the potential rubber crop, TK, and its weedy relatives, we generated genomic and marker resources for these three dandelion species. RESULTS: Complete chloroplast genome sequences of TK (151,338 bp), TO (151,299 bp), and TB (151,282 bp) were obtained using the Illumina GAII and MiSeq platforms. Chloroplast sequences were analyzed and annotated for all the three species. Phylogenetic analysis within Asteraceae showed that TK has a closer genetic distance to TB than to TO and Taraxacum species were most closely related to lettuce (Lactuca sativa). By sequencing multiple genotypes for each species and testing variants using gel-based methods, four chloroplast Single Nucleotide Polymorphism (SNP) variants were found to be fixed between TK and TO in large populations, and between TB and TO. Additionally, Expressed Sequence Tag (EST) resources developed for TO and TK permitted the identification of five nuclear species-specific SNP markers. CONCLUSIONS: The availability of chloroplast genomes of these three dandelion species, as well as chloroplast and nuclear molecular markers, will provide a powerful genetic resource for germplasm differentiation and purification, and the study of potential gene flow among Taraxacum species.

Isoprenoid biosynthesis in dandelion latex is enhanced by the overexpression of three key enzymes involved in the mevalonate pathway.

BACKGROUND: Latex from the dandelion species Taraxacum brevicorniculatum contains many high-value isoprenoid end products, e.g. triterpenes and polyisoprenes such as natural rubber. The isopentenyl pyrophosphate units required as precursors for these isoprenoids are provided by the mevalonate (MVA) pathway. The key enzyme in this pathway is 3-hydroxy-methyl-glutaryl-CoA reductase (HMGR) and its activity has been thoroughly characterized in many plant species including dandelion. However, two enzymes acting upstream of HMGR have not been characterized in dandelion latex: ATP citrate lyase (ACL), which provides the acetyl-CoA utilized in the MVA pathway, and acetoacetyl-CoA thiolase (AACT), which catalyzes the first step in the pathway to produce acetoacetyl-CoA. Here we isolated ACL and AACT genes from T. brevicorniculatum latex and characterized their expression profiles. We also overexpressed the well-characterized HMGR, ACL and AACT genes from Arabidopsis thaliana in T. brevicorniculatum to determine their impact on isoprenoid end products in the latex. RESULTS: The spatial and temporal expression profiles of T. brevicorniculatum ACL and AACT revealed their pivotal role in the synthesis of precursors necessary for isoprenoid biosynthesis in latex. The overexpression of A. thaliana ACL and AACT and HMGR in T. brevicorniculatum latex resulted in the accumulation of all three enzymes, increased the corresponding enzymatic activities and ultimately increased sterol levels by ~5-fold and pentacyclic triterpene and cis-1,4-isoprene levels by ~2-fold. Remarkably high levels of the triterpene precursor squalene were also detected in the triple-transgenic lines (up to 32 mg/g root dry weight) leading to the formation of numerous lipid droplets which were observed in root cross-sections. CONCLUSIONS: We could show the effective expression of up to three transgenes in T. brevicorniculatum latex which led to increased enzymatic activity and resulted in high level squalene accumulation in the dandelion roots up to an industrially relevant amount. Our data provide insight into the regulation of the MVA pathway in dandelion latex and can be used as a basis for metabolic engineering to enhance the production of isoprenoid end products in this specialized tissue.