Dynamic transcriptional profiling provides insights into tuberous root development in Rehmannia glutinosa

Identification of MeC3HDZ1/MeCNA as a potential regulator of cassava storage root development

The storage root (SR) of cassava is the main staple food in sub-Saharan Africa, where it feeds over 500 million people. However, little is known about the genetic and molecular regulation underlying its development. Unraveling such regulation would pave the way for biotechnology approaches aimed at enhancing cassava productivity. Anatomical studies indicate that SR development relies on the massive accumulation of xylem parenchyma, a cell-type derived from the vascular cambium. The C3HDZ family of transcription factors regulate cambial cells proliferation and xylem differentiation in Arabidopsis and other species. We thus aimed at identifying C3HDZ proteins in cassava and determining whether any of them shows preferential activity in the SR cambium and/or xylem. Using phylogeny and synteny studies, we identified eight C3HDZ proteins in cassava, namely MeCH3DZ1-8. We observed that MeC3HDZ1 is the MeC3HDZ gene displaying the highest expression in SR and that, within that organ, the gene also shows high expression in cambium and xylem. In-silico analyses revealed the existence of a number of potential C3HDZ targets displaying significant preferential expression in the SR. Subsequent Y1H analyses proved that MeC3HDZ1 can bind canonical C3HDZ binding sites, present in the promoters of these targets. Transactivation assays demonstrated that MeC3HDZ1 can regulate the expression of genes downstream of promoters harboring such binding sites, thereby demonstrating that MeC3HDZ1 has C3HDZ transcription factor activity. We conclude that MeC3HDZ1 may be a key factor for the regulation of storage root development in cassava, holding thus great promise for future biotechnology applications.

Transcriptome Analysis Provides Insights into Catalpol Biosynthesis in the Medicinal Plant Rehmannia glutinosa and the Functional Characterization of RgGES Genes

Rehmannia glutinosa, a member of the Scrophulariaceae family, has been widely used in traditional Chinese medicine since ancient times. The main bioactive component of R. glutinosa is catalpol. However, the biogenesis of catalpol, especially its downstream pathway, remains unclear. To identify candidate genes involved in the biosynthesis of catalpol, transcriptomes were constructed from R. glutinosa using the young leaves of three cultivars, Beijing No. 3, Huaifeng, and Jin No. 9, as well as the tuberous roots and adventitious roots of the Jin No. 9 cultivar. As a result, 71,142 unigenes with functional annotations were generated. A comparative analysis of the R. glutinosa transcriptomes identified over 200 unigenes of 13 enzymes potentially involved in the downstream steps of catalpol formation, including 9 genes encoding UGTs, 13 for aldehyde dehydrogenases, 70 for oxidoreductases, 44 for CYP450s, 22 for dehydratases, 30 for decarboxylases, 19 for hydroxylases, and 10 for epoxidases. Moreover, two novel genes encoding geraniol synthase (RgGES), which is the first committed enzyme in catalpol production, were cloned from R. glutinosa. The purified recombinant proteins of RgGESs effectively converted GPP to geraniol. This study is the first to discover putative genes coding the tailoring enzymes mentioned above in catalpol biosynthesis, and functionally characterize the enzyme-coding gene in this pathway in R. glutinosa. The results enrich genetic resources for engineering the biosynthetic pathway of catalpol and iridoids.

Transcriptomic analysis of Chinese yam (Dioscorea polystachya Turcz.) variants indicates brassinosteroid involvement in tuber development

Dioscorea is an important but underutilized genus of flowering plants that grows predominantly in tropical and subtropical regions. Several species, known as yam, develop large underground tubers and aerial bulbils that are used as food. The Chinese yam (D. polystachya Turcz.) is one of the few Dioscorea species that grows well in temperate regions and has been proposed as a climate-resilient crop to enhance food security in Europe. However, the fragile, club-like tubers are unsuitable for mechanical harvesting, which is facilitated by shorter and thicker storage organs. Brassinosteroids (BRs) play a key role in plant cell division, cell elongation and proliferation, as well as in the gravitropic response. We collected RNA-Seq data from the head, middle and tip of two tuber shape variants: F60 (long, thin) and F2000 (short, thick). Comparative transcriptome analysis of F60 vs. F2000 revealed 30,229 differentially expressed genes (DEGs), 1,393 of which were differentially expressed in the growing tip. Several DEGs are involved in steroid/BR biosynthesis or signaling, or may be regulated by BRs. The quantification of endogenous BRs revealed higher levels of castasterone (CS), 28-norCS, 28-homoCS and brassinolide in F2000 compared to F60 tubers. The highest BR levels were detected in the growing tip, and CS was the most abundant (439.6 ± 196.41 pmol/g in F2000 and 365.6 ± 112.78 pmol/g in F60). Exogenous 24-epi-brassinolide (epi-BL) treatment (20 nM) in an aeroponic system significantly increased the width-to-length ratio (0.045 ± 0.002) compared to the mock-treated plants (0.03 ± 0.002) after 7 weeks, indicating that exogenous epi-BL produces shorter and thicker tubers. In this study we demonstrate the role of BRs in D. polystachya tuber shape, providing insight into the role of plant hormones in yam storage organ development. We found that BRs can influence tuber shape in Chinese yam by regulating the expression of genes involved cell expansion. Our data can help to improve the efficiency of Chinese yam cultivation, which could provide an alternative food source and thus contribute to future food security in Europe.

Integrative Analysis of the Transcriptome and Metabolome Reveals the Developmental Mechanisms and Metabolite Biosynthesis of the Tuberous Roots of Tetrastigma hemsleyanum

The tuberous root of Tetrastigma hemsleyanum Diels et Gilg (T. hemsleyanum) is a traditional Chinese medicine with a wide range of clinical applications. However, the scarcity of its wild resources, its low yield, and the variable quality that results from its artificial cultivation leads to expensive market prices that are not conducive to the further industrial development of T. hemsleyanum. In this study, transcriptomic and non-targeted metabolomic analyses were integrated to explore the underlying molecular mechanisms and metabolite biosynthesis that occur during its root development. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes (DEGs) were predominantly enriched for processes associated with flavonoid and phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant hormone signal transduction. Genes related to lignin were downregulated in tuberous roots (TRs), resulting in a decrease in lignification and the downregulation of metabolites related to flavonoids and phenylpropanoid biosynthesis. In addition, the expression levels of starch- and sucrose-related genes were upregulated in TRs. The root development of SYQ is also related to IAA, GA, ABA, and JA signaling pathways. Collectively, this study lays the foundation for analyzing the root development and quality-modulating mechanisms employed by T. hemsleyanum; this will be beneficial in conducting molecular-assisted breeding and controlling its secondary metabolite production.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages)

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

The Effects of NAA on the Tuberous Root Yield and Quality of Rehmannia glutinosa and Its Regulatory Mechanism by Transcriptome and Metabolome Profiling

Naphthylacetic acid (NAA) was used to increase the tuberous root yield of Rehmannia glutinosa, but the differences between its NAA-treated and control tuberous roots (NT and CG) and the regulatory mechanism of NAA effect remain unclear. In order to investigate them, NTs and CGs were used as materials, and both yield-related indices were measured; the metabolomics and transcriptomics were used to capture differentially accumulated metabolites (DAM) and to validate them via mining differentially expressed genes (DEGs), respectively. The effects of NAA treatment: increased NT mass per plant by 21.14%, through increasing the number of roots and increasing the mean root diameter; increased catalpol content by 1.2234% (p < 0.05); up-regulated 11DAMs and 596DEGs; and down-regulated 18 DAMs and 517DEGs. In particular, we discovered that NAA regulated its DAMs and biomass via 10 common metabolic pathways, and that the number of NAA-down-regulated DAMs was more than that of NAA-up-regulated DAMs in its tuberous root. Furthermore, HPLC validated the changes of several DAMs and 15 DEGs (4CL, ARF, CCoAOMT, ARGOS, etc.) associated with the yield increase and DAMs were verified by RT-qPCR. This study provided some valuable resources, such as tuberous root indices, key genes, and DAMs of Rehmannia glutinosa in response to NAA for distinguishing the CGs from NTs, and novel insights into the regulatory mechanism of NAA effects on both at the transcriptomic and metabolomic levels, so it will lay a theoretical foundation for NAA-regulated plant yield and quality, and provide references for prohibiting the uses of NAA as a swelling agent in medicinal tuber plants in China.

Set up from the beginning: The origin and early development of cassava storage roots

Despite the importance of storage root (SR) organs for cassava and the other root crops yield, their developmental origin is poorly understood. Here we use multiple approaches to shed light on the initial stages of root development demonstrating that SR and fibrous roots (FR) follow different rhizogenic processes. Transcriptome analysis carried out on roots collected before, during and after root bulking highlighted early and specific activation of a number of functions essential for root swelling and identified root-specific genes able to effectively discriminate emerging FR and SR. Starch and sugars start to accumulate at a higher rate in SR before they swell but only after parenchyma tissue has been produced. Finally, using non-destructive computed tomography measurements, we show that SR (but not FR) contain, since their emergence from the stem, an inner channel structure in continuity with the stem secondary xylem, indicating that SR derive from a distinct rhizogenic process compared with FR.

Roles of Abscisic Acid and Gibberellins in Stem/Root Tuber Development

Root and tuber crops are of great importance. They not only contribute to feeding the population but also provide raw material for medicine and small-scale industries. The yield of the root and tuber crops is subject to the development of stem/root tubers, which involves the initiation, expansion, and maturation of storage organs. The formation of the storage organ is a highly intricate process, regulated by multiple phytohormones. Gibberellins (GAs) and abscisic acid (ABA), as antagonists, are essential regulators during stem/root tuber development. This review summarizes the current knowledge of the roles of GA and ABA during stem/root tuber development in various tuber crops.

Study of Terpenoid Synthesis and Prenyltransferase in Roots of Rehmannia glutinosa Based on iTRAQ Quantitative Proteomics

Rehmannia glutinosa has important medicinal value; terpenoid is one of the main active components in R. glutinosa. In this study, iTRAQ technique was used to analyze the relative abundance of proteins in roots of R. glutinosa, and 6,752 reliable proteins were quantified. GO enrichment results indicated that most proteins were involved in metabolic process or cellular process, 57.63% proteins had catalytic activity, and 65.80% proteins were enriched in membrane-bounded organelle. In roots of R. glutinosa, there were 38 KEGG enrichments with significance, more DEPs were found in some pathways, especially the proteasome pathway and TCA cycle with 15.0% DEPs between elongation stage and expansion stage of roots. Furthermore, five KEGG pathways of terpenoid synthesis were found. Most prenyltransferases belong to FPP/GGPP synthase family, involved in terpenoid backbone biosynthesis, and all interacted with biotin carboxylase CAC2. Compared with that at the elongation stage, many prenyltransferases exhibited higher expression at the expansion stage or maturation stage of roots. In addition, eight FPP/GGPP synthase encoding genes were cloned from R. glutinosa, namely FPPS, FPPS1, GGPS, GGPS3, GGPS4, GGPS5, GPPS and GPPS2, introns were also found in FPPS, FPPS1, GGPS5 and GGPS2, and FPP/GPP synthases were more conservative in organisms, especially in viridiplantae, in which the co-occurrence of GPPS or GPPS2 was significantly higher in plants. Further analysis found that FPP/GGPP synthases of R. glutinosa were divided into three kinds, GGPS, GPPS and FPPS, and their gene expression was significantly diverse in different varieties, growth periods, or tissues of R. glutinosa. Compared with that of GGPS, the expression of GPPS and FPPS was much higher in R. glutinosa, especially at the expansion stage and maturation stage. Thus, the synthesis of terpenoids in roots of R. glutinosa is intricately regulated and needs to be further studied.

Over-immunity mediated abnormal deposition of lignin arrests the normal enlargement of the root tubers of Rehmannia glutinosa under consecutive monoculture stress

The rapid accumulation of lignin in the cell wall is one of important immune defense mechanism in response to adversity stress in plants. In this study, we found that the enlargement of the root tubers of Rehmannia glutinosa (R. glutinosa) is arrested under consecutive monoculture stress, and this process is accompanied by abnormal accumulation of lignin. Meanwhile, the function of key catalytic enzyme genes in lignin biosynthetic pathway under consecutive monoculture stress was systematically analyzed, of which roles of core genes were validated using reverse genetics. We elucidated that an abnormal deposition of lignin in R. glutinosa roots, induced by consecutive monoculture stress, and arrested the enlargement of root tubers. Additionally, by manipulating the key catalytic enzyme gene RgCCR6, we were able to alter lignin content of roots of R. glutinosa, thereby affecting tuber enlargement. We speculate that cell lignification is an important defense strategy in resistance against consecutive monoculture stress, but the overreacted defense hindered the normal enlargement of root tubers. The findings provide new insights for effectively improving yield reductions of root crops subjected to environmental stress.

Revealing Microbiome Structure and Assembly Process in Three Rhizocompartments of Achyranthes bidentata Under Continuous Monoculture Regimes

The complex composition and interaction of root-associated microbes are critical to plant health and performance. In this study, we presented a detailed characterization of three rhizocompartment (rhizosphere, rhizoplane, and root) microbiomes of Achyranthes bidentata under different years of consecutive monoculture by deep sequencing in order to determine keystone microorganisms via co-occurrence network analysis. The network analysis showed that multiple consecutive monoculture (MCM, represented 5Y and 10Y) soils generated some distinct beneficial bacterial taxa such as Bacillus, Fictibacillus, Bradyrhizobium, Shinella, and Herbaspirillum. For fungi, Mortierella substituted for Fusarium in occupying an important position in different rhizocompartments under A. bidentate monoculture. Quantitative PCR analysis confirmed a significant increase in Bacillus, Pseudomonas, and Burkholderia spp. The results of the inoculation assay showed that addition of beneficial bacteria Bacillus subtilis 74 and Bacillus halodurans 75 significantly increased the root length and fresh weight of A. bidentata. Furthermore, three types of phytosterones, as the main allochemicals, were identified both in the rhizosphere soil and in culture medium under sterile conditions by LC-MS/MS. When looking at in vitro interactions, it was found that phytosterones displayed a positive interaction with dominant beneficial species (Bacillus amyloliquefaciens 4 and B. halodurans 75) and had a negative effect on the presence of the pathogenic fungi Fusarium solani and Fusarium oxysporum. Overall, this study demonstrated that consecutive monoculture of A. bidentata can alter the bacterial and fungal community by secreting root exudates, leading to recruitment of beneficial microbes and replacement of plant-specific pathogenic fungi with plant beneficial fungi.

Hormonal and transcriptional analyses provides new insights into the molecular mechanisms underlying root thickening and isoflavonoid biosynthesis in Callerya speciosa (Champ. ex Benth.) Schot

Callerya speciosa (Champ. ex Benth.) Schot is a traditional Chinese medicine characterized by tuberous roots as the main organ of isoflavonoid accumulation. Root thickening and isoflavonoid accumulation are two major factors for yield and quality of C. speciosa. However, the underlying mechanisms of root thickening and isoflavonoid biosynthesis have not yet been elucidated. Here, integrated morphological, hormonal and transcriptomic analyses of C. speciosa tuberous roots at four different ages (6, 12, 18, 30 months after germination) were performed. The growth cycle of C. speciosa could be divided into three stages: initiation, rapid-thickening and stable-thickening stage, which cued by the activity of vascular cambia. Endogenous changes in phytohormones were associated with developmental changes during root thickening. Jasmonic acid might be linked to the initial development of tuberous roots. Abscisic acid seemed to be essential for tuber maturation, whereas IAA, cis-zeatin and gibberellin 3 were considered essential for rapid thickening of tuberous roots. A total of 4337 differentially expressed genes (DEGs) were identified during root thickening, including 15 DEGs participated in isoflavonoid biosynthesis, and 153 DEGs involved in starch/sucrose metabolism, hormonal signaling, transcriptional regulation and cell wall metabolism. A hypothetical model of genetic regulation associated with root thickening and isoflavonoid biosynthesis in C. speciosa is proposed, which will help in understanding the underlying mechanisms of tuberous root formation and isoflavonoid biosynthesis.

Seasonal Variation in Transcriptomic Profiling of Tetrastigma hemsleyanum Fully Developed Tuberous Roots Enriches Candidate Genes in Essential Metabolic Pathways and Phytohormone Signaling

Tetrastigma hemsleyanum Diels et Gilg (Sanyeqing, SYQ) is a perennial climbing liana and an endemic plant to southern China. Its tuberous roots (TRs) are used in traditional Chinese medicine for treating some diseases such as high fever, pneumonia, asthma, hepatitis, and cancers. However, the mechanisms underlying the development of TR and the content of flavonoids and phenylpropanoids (FPs) are not well-understood. In this study, we performed a transcriptomic analysis of 12 fully developed TR (FD-TR) samples harvested in four seasons [spring (Sp), summer (Su), autumn (Au), and winter (Wi)] using the RNA-Sequencing (RNA-Seq). We obtained a total of 78.54 Gb raw data and 65,578 unigenes. Then, the unigenes were annotated by using six databases such as non-redundant protein database (NR), Pfam, eggNOG, SWISSProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene ontology (GO). The transcriptomic profiling showed closer relationships between the samples obtained in Su and Au than those obtained in Sp and Wi based on the results of both total unigenes and differentially expressed genes (DEGs). Three pathways, including the biosynthesis of FPs, metabolism of starch and sucrose, and signaling of phytohormones, were highly enriched, suggesting a gene-level seasonal variation. Based on the numbers of DEGs, brassinosteroid (BR) signal transduction factors appeared to play a key role in modulating the development of TRs while most of the auxin signaling genes were mainly activated in Wi and Sp FD-TRs. Most genes in the biosynthesis and biodegradation of starch and biodegradation of cellulose were activated in Wi FD-TRs. As determined by the high performance liquid chromatography (HPLC) and aluminum nitrate colorimetric method, the contents of total flavonoids and most detected FP components increased from Sp to Au but decreased in Wi. Enhanced expression levels of some genes in the biosynthetic pathways of FPs were detected in Su and Au samples, which corroborated well with metabolite content. Our findings provide the first transcriptomic and biochemical data on a seasonal variation in the composition of medically important metabolites in SYQ FD-TRs.

Root diversity in sesame (Sesamum indicum L.): insights into the morphological, anatomical and gene expression profiles

Sesame harbors a large diversity in root morphological and anatomical traits and a high root biomass improves the plant aboveground biomass as well as the seed yield. Sesame provides one of the most nutritious and healthy vegetable oils, sparking an increasing demand of its seeds. However, with the low yield and productivity of sesame, there is still a huge gap between the seed demand and supply. Improving the root system has a high potential to increase crop productivity, but information on the diversity of the sesame root systems is still lacking. In this study, 40 diverse sesame varieties were grown in soil and hydroponics systems and the diversity of the root system was investigated. The results showed that sesame holds a large root morphological and anatomical diversity, which can be harnessed in breeding programmes. Based on the clustering of the genotypes in hydroponics and soil culture systems, we found that similar genotypes were commonly clustered either in the small-root or in the big-root group, indicating that the hydroponics system can be employed for a large-scale root phenotyping. Our results further revealed that the root biomass positively contributes to increased seed yield in sesame, based on multi-environmental trials. By comparing the root transcriptome of two contrasting genotypes, 2897 differentially expressed genes were detected and they were enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, stilbenoid, diarylheptanoid and gingerol biosynthesis, flavonoid biosynthesis, suggesting that these pathways are crucial for sesame root growth and development. Overall, this study sheds light on the diversity of sesame root system and offers the basis for improving root traits and increasing sesame seed yield.

Molecular Regulation of Catalpol and Acteoside Accumulation in Radial Striation and non-Radial Striation of Rehmannia glutinosa Tuberous Root

Rehmannia glutinosa L., a perennial plant of Scrophulariaceae, is one of the most commonly used herbs in traditional Chinese medicine (TCM) that have been widely cultivated in China. However, to date, the biosynthetic pathway of its two quality-control components, catalpol and acteoside, are only partially elucidated and the mechanism for their tissue-specific accumulation remains unknown. To facilitate the basic understanding of the key genes and transcriptional regulators involved in the biosynthesis of catalpol and acteoside, transcriptome sequencing of radial striation (RS) and non-radial striation (nRS) from four R. glutinosa cultivars was performed. A total of 715,158,202 (~107.27 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 150,405 transcripts. Functional annotation with multiple public databases identified 155 and 223 unigenes involved in catalpol and acteoside biosynthesis, together with 325 UGTs, and important transcription factor (TF) families. Comparative analysis of the transcriptomes identified 362 unigenes, found to be differentially expressed in all RS vs. nRS comparisons, with 143 upregulated unigenes, including those encoding enzymes of the catalpol and acteoside biosynthetic pathway, such as geranyl diphosphate synthase (RgGPPS), geraniol 8-hydroxylase (RgG10H), and phenylalanine ammonia-lyase (RgPAL). Other differentially expressed unigenes predicted to be related to catalpol and acteoside biosynthesis fall into UDP-dependent glycosyltransferases (UGTs), as well as transcription factors. In addition, 16 differentially expressed genes were selectively confirmed by real-time PCR. In conclusion, a large unigene dataset of R. glutinosa generated in the current study will serve as a resource for the identification of potential candidate genes for investigation of the tuberous root development and biosynthesis of active components.

Metabolite accumulation and metabolic network in developing roots of Rehmannia glutinosa reveals its root developmental mechanism and quality

Rehmannia glutinosa root contains many compounds with important medicinal properties and nutritional benefits, but only more than 140 compounds have been reported so far. Many other compounds and their accumulation and metabolic networks during its development remain unclear. In order to clarify them, its metabolic profiles at three different developmental stages were analyzed using untargeted LC-MS analysis. Multivariate analysis revealed that 434 metabolites differently accumulated in its different stages, suggesting different change trends. The metabolites having the same trend share common metabolic pathways, the metabolites showing increasing contents during its development have medical and nutritional values, and some mature root-specific metabolites may be better candidates for its quality control; 434 metabolites were mapped onto 111 KEGG pathways including 62 enzymes, whose increasing and decreasing patterns were shown during its development. Some metabolites complicatedly interacted with some enzymes and the top-10 pathways enriched from 111 KEGG pathways in network analysis. These findings extended the dataset of its identified compounds, and revealed that its development and quality were associated with the accumulation of different metabolites. Our work will lay the foundation for the better understanding of its chemical constituents, quality and developmental mechanism.

Differential proteomic analysis of replanted Rehmannia glutinosa roots by iTRAQ reveals molecular mechanisms for formation of replant disease

BACKGROUND

The normal growth of Rehmannia glutinosa, a widely used medicinal plant in China, is severely disturbed by replant disease. The formation of replant disease commonly involves interactions among plants, allelochemicals and microbes; however, these relationships remain largely unclear. As a result, no effective measures are currently available to treat replant disease.

RESULTS

In this study, an integrated R. glutinosa transcriptome was constructed, from which an R. glutinosa protein library was obtained. iTRAQ technology was then used to investigate changes in the proteins in replanted R. glutinosa roots, and the proteins that were expressed in response to replant disease were identified. An integrated R. glutinosa transcriptome from different developmental stages of replanted and normal-growth R. glutinosa produced 65,659 transcripts, which were accurately translated into 47,818 proteins. Using this resource, a set of 189 proteins was found to be significantly differentially expressed between normal-growth and replanted R. glutinosa. Of the proteins that were significantly upregulated in replanted R. glutinosa, most were related to metabolism, immune responses, ROS generation, programmed cell death, ER stress, and lignin synthesis.

CONCLUSIONS

By integrating these key events and the results of previous studies on replant disease formation, a new picture of the damaging mechanisms that cause replant disease stress emerged. Replant disease altered the metabolic balance of R. glutinosa, activated immune defence systems, increased levels of ROS and antioxidant enzymes, and initiated the processes of cell death and senescence in replanted R. glutinosa. Additionally, lignin deposition in R. glutinosa roots that was caused by replanting significantly inhibited tuberous root formation. These key processes provide important insights into the underlying mechanisms leading to the formation of replant disease and also for the subsequent development of new control measures to improve production and quality of replanted plants.

Transcriptome Analysis of Salicylic Acid Treatment in Rehmannia glutinosa Hairy Roots Using RNA-seq Technique for Identification of Genes Involved in Acteoside Biosynthesis

Rehmannia glutinosa is a common bulk medicinal material that has been widely used in China due to its active ingredients. Acteoside, one of the ingredients, has antioxidant, antinephritic, anti-inflammatory, hepatoprotective, immunomodulatory, and neuroprotective effects, is usually selected as a quality-control component for R. glutinosa herb in the Chinese Pharmacopeia. The acteoside biosynthesis pathway in R. glutinosa has not yet been clearly established. Herein, we describe the establishment of a genetic transformation system for R. glutinosa mediated by Agrobacterium rhizogenes. We screened the optimal elicitors that markedly increased acteoside accumulation in R. glutinosa hairy roots. We found that acteoside accumulation dramatically increased with the addition of salicylic acid (SA); the optimal SA dose was 25 μmol/L for hairy roots. RNA-seq was applied to analyze the transcriptomic changes in hairy roots treated with SA for 24 h in comparison with an untreated control. A total of 3,716, 4,018, and 2,715 differentially expressed transcripts (DETs) were identified in 0 h-vs.-12 h, 0 h-vs.-24 h, and 12 h-vs.-24 h libraries, respectively. KEGG pathway-based analysis revealed that 127 DETs were enriched in "phenylpropanoid biosynthesis." Of 219 putative unigenes involved in acteoside biosynthesis, 54 were found to be up-regulated at at least one of the time points after SA treatment. Selected candidate genes were analyzed by quantitative real-time PCR (qRT-PCR) in hairy roots with SA, methyl jasmonate (MeJA), AgNO₃ (Ag⁺), and putrescine (Put) treatment. All genes investigated were up-regulated by SA treatment, and most candidate genes were weakly increased by MeJA to some degree. Furthermore, transcription abundance of eight candidate genes in tuberous roots of the high-acteoside-content (HA) cultivar QH were higher than those of the low-acteoside-content (LA) cultivar Wen 85-5. These results will pave the way for understanding the molecular basis of acteoside biosynthesis in R. glutinosa, and can serve as a basis for future validation studies.

De Novo Sequencing and Assembly Analysis of the Pseudostellaria heterophylla Transcriptome

Pseudostellaria heterophylla (Miq.) Pax is a mild tonic herb widely cultivated in the Southern part of China. The tuberous roots of P. heterophylla accumulate high levels of secondary metabolism products of medicinal value such as saponins, flavonoids, and isoquinoline alkaloids. Despite numerous studies on the pharmacological importance and purification of these compounds in P. heterophylla, their biosynthesis is not well understood. In the present study, we used Illumina HiSeq 4000 sequencing platform to sequence the RNA from flowers, leaves, stem, root cortex and xylem tissues of P. heterophylla. We obtained 616,413,316 clean reads that we assembled into 127, 334 unique sequences with an N50 length of 951 bp. Among these unigenes, 53,184 unigenes (41.76%) were annotated in a public database and 39, 795 unigenes were assigned to 356 KEGG pathways; 23,714 unigenes (8.82%) had high homology with the genes from Beta vulgaris. We discovered 32, 095 DEGs in different tissues and performed GO and KEGG enrichment analysis. The most enriched KEGG pathway of secondary metabolism showed up-regulated expression in tuberous roots as compared with the ground parts of P. heterophylla. Moreover, we identified 72 candidate genes involved in triterpenoids saponins biosynthesis in P. heterophylla. The expression profiles of 11 candidate unigenes were analyzed by quantitative real-time PCR (RT-qPCR). Our study established a global transcriptome database of P. heterophylla for gene identification and regulation. We also identified the candidate unigenes involved in triterpenoids saponins biosynthesis. Our results provide an invaluable resource for the secondary metabolites and physiological processes in different tissues of P. heterophylla.

Identification of candidate domestication regions in the radish genome based on high-depth resequencing analysis of 17 genotypes

This study provides high-quality variation data of diverse radish genotypes. Genome-wide SNP comparison along with RNA-seq analysis identified candidate genes related to domestication that have potential as trait-related markers for genetics and breeding of radish. Radish (Raphanus sativus L.) is an annual root vegetable crop that also encompasses diverse wild species. Radish has a long history of domestication, but the origins and selective sweep of cultivated radishes remain controversial. Here, we present comprehensive whole-genome resequencing analysis of radish to explore genomic variation between the radish genotypes and to identify genetic bottlenecks due to domestication in Asian cultivars. High-depth resequencing and multi-sample genotyping analysis of ten cultivated and seven wild accessions obtained 4.0 million high-quality homozygous single-nucleotide polymorphisms (SNPs)/insertions or deletions. Variation analysis revealed that Asian cultivated radish types are closely related to wild Asian accessions, but are distinct from European/American cultivated radishes, supporting the notion that Asian cultivars were domesticated from wild Asian genotypes. SNP comparison between Asian genotypes identified 153 candidate domestication regions (CDRs) containing 512 genes. Network analysis of the genes in CDRs functioning in plant signaling pathways and biochemical processes identified group of genes related to root architecture, cell wall, sugar metabolism, and glucosinolate biosynthesis. Expression profiling of the genes during root development suggested that domestication-related selective advantages included a main taproot with few branched lateral roots, reduced cell wall rigidity and favorable taste. Overall, this study provides evolutionary insights into domestication-related genetic selection in radish as well as identification of gene candidates with the potential to act as trait-related markers for background selection of elite lines in molecular breeding.

Transcriptome Analysis of Storage Roots and Fibrous Roots of the Traditional Medicinal Herb Callerya speciosa (Champ.) ScHot

Callerya speciosa (Champ.) ScHot is a woody perennial plant in Fabaceae, the roots of which are used medicinally. The storage roots of C. speciosa are derived from fibrous roots, but not all fibrous roots can develop into storage roots. To detect key genes involved in storage roots formation, we performed Illumina sequencing of the C. speciosa storage roots and fibrous roots. De novo assembly resulted in 161,926 unigenes, which were subsequently annotated by BLAST, GO and KEGG analyses. After expression profiling, 4538 differentially expressed genes were identified. The KEGG pathway enrichment analysis revealed changes in the biosynthesis of cytokinin, phenylpropanoid, starch, sucrose, flavone and other secondary metabolites. Transcription factor-related differentially expressed genes (DEGs) were also identified, including such gene families as GRAS, COL, MIKC, ERF, LBD, and NAC. The DEGs related to light signaling, starch, sugar, photohormones and cell wall-loosening might be involved in the formation of storage roots. This study provides the first transcriptome profiling of C. speciosa roots, data that will facilitate future research of root development and metabolites with medicinal value as well as the breeding of C. speciosa.

Comparative Transcriptome Analysis of Primary Roots of Brassica napus Seedlings with Extremely Different Primary Root Lengths Using RNA Sequencing

Primary root (PR) development is a crucial developmental process that is essential for plant survival. The elucidation of the PR transcriptome provides insight into the genetic mechanism controlling PR development in crops. In this study, we performed a comparative transcriptome analysis to investigate the genome-wide gene expression profiles of the seedling PRs of four Brassica napus genotypes that were divided into two groups, short group (D43 and D61), and long group (D69 and D72), according to their extremely different primary root lengths (PRLs). The results generated 55,341,366-64,631,336 clean reads aligned to 62,562 genes (61.9% of the current annotated genes) in the B. napus genome. We provide evidence that at least 44,986 genes are actively expressed in the B. napus PR. The majority of the genes that were expressed during seedling PR development were associated with metabolism, cellular processes, response to stimulus, biological regulation, and signaling. Using a pairwise comparison approach, 509 differentially expressed genes (DEGs; absolute value of log2 fold-change ≥1 and p ≤ 0.05) between the long and short groups were revealed, including phytohormone-related genes, protein kinases and phosphatases, oxygenase, cytochrome P450 proteins, etc. Combining GO functional category, KEGG, and MapMan pathway analyses indicated that the DEGs involved in cell wall metabolism, carbohydrate metabolism, lipid metabolism, secondary metabolism, protein modification and degradation, hormone pathways and signaling pathways were the main causes of the observed PRL differences. We also identified 16 differentially expressed transcription factors (TFs) involved in PR development. Taken together, these transcriptomic datasets may serve as a foundation for the identification of candidate genes and may provide valuable information for understanding the molecular and cellular events related to PR development.

Comparative Transcriptomic Analyses of Vegetable and Grain Pea (Pisum sativum L.) Seed Development

Understanding the molecular mechanisms regulating pea seed developmental process is extremely important for pea breeding. In this study, we used high-throughput RNA-Seq and bioinformatics analyses to examine the changes in gene expression during seed development in vegetable pea and grain pea, and compare the gene expression profiles of these two pea types. RNA-Seq generated 18.7 G of raw data, which were then de novo assembled into 77,273 unigenes with a mean length of 930 bp. Our results illustrate that transcriptional control during pea seed development is a highly coordinated process. There were 459 and 801 genes differentially expressed at early and late seed maturation stages between vegetable pea and grain pea, respectively. Soluble sugar and starch metabolism related genes were significantly activated during the development of pea seeds coinciding with the onset of accumulation of sugar and starch in the seeds. A comparative analysis of genes involved in sugar and starch biosynthesis in vegetable pea (high seed soluble sugar and low starch) and grain pea (high seed starch and low soluble sugar) revealed that differential expression of related genes at late development stages results in a negative correlation between soluble sugar and starch biosynthetic flux in vegetable and grain pea seeds. RNA-Seq data was validated by using real-time quantitative RT-PCR analysis for 30 randomly selected genes. To our knowledge, this work represents the first report of seed development transcriptomics in pea. The obtained results provide a foundation to support future efforts to unravel the underlying mechanisms that control the developmental biology of pea seeds, and serve as a valuable resource for improving pea breeding.