Original Research ARTICLE
Dynamic transcriptional profiling provides insights into tuberous root development in Rehmannia glutinosa
- 1Center for Medicinal Plant Cultivation, Institute of Medicinal Plant Development, Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China
- 2State Key Laboratory of Plant Physiology and Biochemistry, College of Biological Sciences, China Agricultural University, Beijing, China
- 3Department of Agronomy, College of Agriculture and Biotechnology, China Agricultural University, Beijing, China
Rehmannia glutinosa, an herb of the Scrophulariaceae family, is widely cultivated in the Northern part of China. The tuberous root has well-known medicinal properties; however, yield and quality are threatened by abiotic and biotic stresses. Understanding the molecular process of tuberous root development may help identify novel targets for its control. In the present study, we used Illumina sequencing and de novo assembly strategies to obtain a reference transcriptome that is relevant to tuberous root development. We then conducted RNA-seq quantification analysis to determine gene expression profiles of the adventitious root (AR), thickening adventitious root (TAR), and the developing tuberous root (DTR). Expression profiling identified a total of 6794 differentially expressed unigenes during root development. Bioinformatics analysis and gene expression profiling revealed changes in phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant hormone biosynthesis during root development. Moreover, we identified and allocated putative functions to the genes involved in tuberous root development, including genes related to major carbohydrate metabolism, hormone metabolism, and transcription regulation. The present study provides the initial description of gene expression profiles of AR, TAR, and DTR, which facilitates identification of genes of interest. Moreover, our work provides insights into the molecular mechanisms underlying tuberous root development and may assist in the design and development of improved breeding schemes for different R. glutinosa varieties through genetic manipulation.
Rehmannia glutinosa, commonly known as the Chinese foxglove (Scrophulariaceae), is native to China, Japan, and Korea. The Chinese medical classic “Shen Nong Ben Cao Jing” (the earliest monograph on Chinese materia medica) lists R. glutinosa as a “top grade” herb (Yang, 2007). Traditional Chinese medicine uses three types of tuberous roots, namely, fresh, dried, and treated, to nourish the Yin (the negative, dark, and feminine, according to Chinese philosophy), abate pathogenic heat, strengthen the body, and replenish vitality (Li, 2014). Modern pharmacological studies indicate that Rehmannia extracts exert broad-spectrum pharmacological effects on various human systems such as the circulatory, endocrine, cardiovascular, nervous, and immune systems (Zhang et al., 2008). There has been a surge in the use of Rehmannia in treating ailments such as angiocardiopathy, diabetes, hypertension, vascular dementia, cerebral infarction, urinary tract problems, and menstrual disorders (Li, 2014).
R. glutinosa has been commercially cultivated in China for over a thousand years, and several high yielding cultivars such as “Beijing No. 1” and “85-5” were developed by Chinese breeders through cross breeding and selection. The cultivar “85-5,” which was developed in the 1980s, is characterized by big tuberous roots and high yield. It has an average fresh weight of 24.5 tons per hectare (Li et al., 2007), which is comparable with that of sweet potato, a globally important food crop. Under high-density planting and plastic film sheeting conditions, the yield of “85-5” can even reach 95.0 tons per hectare (Li et al., 2009a). Rehmannia has a high yield potential; it suffers a significant reduction in tuberous root yield due to viral infections and replanting diseases (Wen et al., 2001). Replanting occurs during consecutive plantings and causes a severe decline in the productivity and quality of the tuberous products (Liu et al., 2006; Zhang et al., 2011). Researchers have studied the anatomical and hormonal changes associated with induced tuberous root formation in vitro to combat diseases and improve its productivity (Xue et al., 2004a,b). The molecular basis of replanting diseases had piqued the interest of researchers (Yang et al., 2011, 2014; Li et al., 2013); however, research studies on this particular area are limited. Understanding of the molecular process of tuberous root development will greatly facilitate the genetic improvement of Rehmannia in terms of yield, quality, and medicinal values.
Despite the extensive research efforts to identify genes involved in tuberous root development, the molecular mechanism underlying tuberous root development in Rehmannia remains unclear (Sun et al., 2010). The limited genomic information on Rehmannia consists of 1301 EST sequences in the GenBank database (as of December 2014). The lack of extensive transcriptomic and functional genomic resources is one of the major limiting factors of molecular research using non-model species. RNA-seq is a useful tool for studying differentially expressed transcripts in various tissues as well as developmental stages. Large-scale transcriptome sequencing of the tuberous roots of the sweet potato has facilitated gene discovery and has provided sequence resources for further molecular research (Schafleitner et al., 2010; Wang et al., 2010; Tao et al., 2012; Firon et al., 2013). Next-generation sequencing technology has also helped understand the molecular basis of replanting disease (Yang et al., 2011, 2014; Li et al., 2013). Moreover, it enables researchers to identify putative genes that might be involved in iridoid biosynthesis in R. glutinosa (Sun et al., 2012). However, the dynamic root transcriptome profile and the functional genes related to tuberous root development have not yet been explored. Therefore, the present study was conducted to improve our understanding of the molecular basis of tuberous root development in Rehmannia. We performed the first global analysis of the Rehmannia transcriptome during tuberous root development using the Illumina RNA-Seq method to identify genes related to sugar metabolism, hormone metabolism, and transcriptional regulation.
Materials and Methods
Plant Materials and RNA Extraction
R. glutinosa plants (cultivar 85-5) were planted in pots containing a ratio of 1:1 vermiculite and humus soil, and were grown in a natural light-greenhouse. The root samples were collected at three stages: adventitious root (AR), thickening adventitious root (TAR), and the developing tuberous root (DTR), sampled at 30, 45, and 60 days after sprouting (DAS), respectively. For consistency, the roots were classified according to diameter: AR ≤ 2 mm, 3 mm ≤ TAR ≤ 5 mm, and 1 cm ≤ DTR ≤ 5 cm (Figure S1). Leaf and stem samples were collected at 60 DAS. Flower samples were collected during the floral season from 2-year-old plants. Five individual plants at various root developmental stages were collected and pooled together. Total RNAs were extracted using the Trizol™ reagent (Invitrogen, Carlsbad, CA, USA) and treated with DNase I (TaKaRa, Dalian, China).
Illumina Sequencing, de novo Assembly, and Functional Annotation
To obtain a global view of the transcriptome that is relevant to tuberous root development in R. glutinosa, a reference transcriptome was built as described below. First, RNAs extracted from the ARs, TARs, and DTRs were equally pooled and sequenced on an Illumina HiSeq2000™ platform according to the manufacture's protocol. First-strand cDNA synthesis was performed by using random hexamer primers, whereas second-strand cDNA synthesis was conducted with buffer, dNTPs, RNaseH, and DNA polymerase I. These products were purified by using agarose gel electrophoresis and enriched by PCR to create the final cDNA library. The Beijing Genomics Institute commercially sequenced the cDNA library on an Illumina HiSeq2000™ platform, using paired-end 90 base pair (bp) reads to generate the raw sequences. The generated data are available at the NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra/SRX357355[accn]). Sequence assembly was conducted using the short reads assembling program SOAPdenovo, following the provided software instructions (Li et al., 2010). Second, the unique sequences produced by the SOAPdenovo pipeline were combined with the Rehmannia transcripts generated on the 454 GS FLX Titanium platform in a previous study (Sun et al., 2012) by using the CAP3 program (Huang and Madan, 1999). Assembly criteria included a 100-bp minimum match and a 95% minimum identity at the overlapping region. Third, all unigenes were then annotated by BLAST against various protein databases [i.e., Nr, SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Cluster of Orthologous Groups of proteins (COG)], with an e-value cutoff of 1e-5. All unigenes generated in the present study were subjected to Gene Ontology (GO) functional analysis. We used the Blast2GO program (Conesa and Götz, 2008) to identify GO terms for all assembled unigenes. Then, the WEGO software was used to perform GO functional classification for all unigenes to determine the distribution of gene functions of the species at the macromolecular level.
RNA-seq Quantification and KEGG Pathway Enrichment Analysis
Three independent cDNA libraries prepared from three developmental staged roots (AR, TAR, and DTR) were subjected to RNA-Seq analysis. We generated 7.25 million (M), 7.13 M, and 7.50 M clean sequence reads for each cDNA library, with each read approximately 49-bp in length. The raw reads were deposited to the NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/Traces/sra/) under Accession numbers SRX700627, SRX700629, and SRX700631. The clean reads were then aligned to the reference transcriptome using SOAPaligner/soap2, allowing two-base mismatches. Reads that mapped to multiple locations were excluded from the analysis. Only unique matched reads were collected for statistical analysis of gene expression. The gene expression level of every unigene was calculated using the reads per kb per million reads (RPKM) method. The differentially expressed genes (DEG) in the ARs, TARs, and DTRs were analyzed as described elsewhere (Audic and Claverie, 1997), with some modifications. False discovery rate (FDR) was used to control the p-values. A threshold value of FDR ≤ 0.001 and an absolute value of log2Ratio ≥ 1 were used to judge the significance of the differences in gene expression. The DEGs were subjected to pathway enrichment analysis by comparing the DGE number (in a specific pathway) with that of the whole genome. The hypergeometric test was used to find significantly enriched pathways. The calculated p-value was submitted for Bonferroni correction, taking the corrected p-value ≤ 0.05 as threshold. Pathways fulfilling these criteria were defined as significantly enriched pathways.
Identification of Transcription Factors (TFs) in R. glutinosa
Rehmannia TFs were searched and organized based on their sequence homology in PlantTFDB3.0 (http://planttfdb.cbi.pku.edu.cn/, Center for Bioinformatics, Peking University, China), which provides comprehensive information on various plant TF families (Jin et al., 2014).
Quantitative Real-Time PCR
First-strand cDNA synthesis was performed with 500 ng of total RNA using the PrimeScript RT® reagent kit with gDNA Eraser (TaKaRa, Dalian, China), according to the manufacturer's instruction. The transcriptional profiles of unigenes were analyzed by qRT-PCR by using the SYBR® Premix Ex Taq II (TaKaRa, Dalian, China) and a real-time PCR system (Bio-Rad, USA). Gene-specific primers were designed by using the software Primer premier 5. The PCR conditions were as follows: a 95°C hold for 3 min; followed 95°C for 15 s, 58°C for 30 s, and 72°C for 20 s for a total of 40 cycles. Gene expression levels were normalized against the internal reference gene, TIP41. The relative expression ratio was calculated for each root sample relative to that of the ARs by using the 2−ΔΔCt method.
Quantification of Endogenous Plant Hormones
Based on the significant changes in the levels of expression of hormone biosynthetic genes during root development, we examined the contents of plant hormones to determine their roles in tuberous root development. The levels of endogenous indole-3-acetic acid (IAA), gibberellic acid (GA3), zeatin riboside (ZR), and abscisic acid (ABA) were determined by using the method of Teng et al. (2006). The root samples were cleaned, weighed, flash frozen in liquid nitrogen, and stored at −80°C until analysis. The samples were ground in liquid nitrogen by using a mortar and pestle, extracted with ice-cold 80% methanol (v/v) containing 1 mol/L butylated hydroxytoluene to prevent oxidation, and then stored overnight at 4°C. The extracts were then centrifuged at 12,000 g for 15 min at 4°C. The residues were resuspended in the same ice-cold extraction solution and stored at 4°C for 1 h, and then centrifuged again at 12,000 g for 15 min at 4°C. The two resulting supernatants were combined and passed through a C-18 Sep-Pak cartridge (Waters, Milford, MA, USA). The efflux was collected and dried under nitrogen. The residues were then dissolved in 0.01 mol/L phosphate buffer solution (pH 7.4) and the concentrations of IAA, GA3, ZR, and ABA were determined by using enzyme-linked immunosorbent assay (ELISA), as desribed elsewhere (Teng et al., 2006). All the measurements were performed with three biological replicates for each development stage, and at least three samples were pooled together for each replicate.
Establishment of a Root Reference Transcriptome for Rehmannia
We generated a total of 81,003 unique sequences from the cDNA library constructed from the ARs, TARs, and DTRs of Rehmannia by using the Illumina sequencing platform. We combined the 81,003 unique sequences with the transcripts generated by 454 GS FLX Titanium platform in a previous study by using the Cap3 program to increase the representativeness of the reference sequences. Our results identified a total of 93,172 unigenes with an average length of 383-bp that served as reference for subsequent analyses.
Alignment of unigenes with the nine full-length cDNA sequences resulted in eight significant hits in our assembly. These hits were well-covered by the unigenes (Supplementary Table S1), indicating that the assembly was well-represented at the genomic level and could be a useful source for gene discovery and global gene expression analysis.
Functional Annotation Analysis
Gene annotation showed that of the 93,172 unigenes, 54,105 (58.07%) had significant matches in the public databases. The annotation ratio was similar to those of previous reports (Yang et al., 2011, 2014; Sun et al., 2012; Li et al., 2013). This transcriptome annotation of unigenes provides valuable information for the functional assignment of genes as well as for novel gene discovery.
GO term assignment helped assign at least one GO term to a total of 24,428 unigenes. All GO terms were classified into 44 functional groups based on biological processes, cellular components, and molecular function categories. The GO analysis provided comprehensive information on unigene functions of R. glutinosa. Among the biological processes, unigenes assigned to metabolic (11,712, 47.94%) and cellular (10,172, 41.64%) processes were the most abundant. A high percentage of genes were grouped into the “response to stimuli” (3785, 15.49%), “biological regulation” (2646, 10.83%), and “developmental process” (2119, 8.67%) categories. These categories feature signal induction, perception, transduction, and modulation of functional genes, which are the main steps of plant root growth in response to environmental and nutrient stimuli, and eventually determine the root fate. These genes are potentially important topics that could be examined in research studies on root development (Figure 1, Supplementary Table S2).
Metabolic pathway analysis assigned 23,171 significantly matched unigenes to the KEGG pathways. Among these, a total of 5644 unigenes were assigned to the metabolic pathways of carbohydrate metabolism, amino acid metabolism, biosynthesis of secondary metabolites, lipid metabolism, and metabolism of terpenoids and polyketides. The tuberous root stores energy (in the form of carbohydrates), as well as nutrients that are essential for growth and reproduction in the subsequent season. Carbohydrate metabolism was the most represented pathway during tuberous root development (Figure 2A). The unigenes involved in carbohydrate metabolism were further grouped into 12 sub-pathways. Among these, starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, amino sugar and nucleotide sugar metabolism, and galactose metabolism were the top five represented groups (Figure 2B).
Figure 2. Kyoto Encyclopedia of Genes and Genomes classification of assembled genes. (A) Genes mapped to metabolic pathways; (B) genes involved in carbohydrate metabolism.
Analysis of Differentially Expressed Genes
The cDNA libraries from ARs, TARs, and DTRs, respectively generated 7.25, 7.13, and 7.50 M clean reads that were then aligned to the reference sequences to map 84,746 unigenes (Figure 3A). Most of these expressed genes (66,618 out of 84,746) were present in all three libraries. We used the RPKM method to calculate the gene expression values of each unigene in the three libraries. A pairwise comparison of the gene expression levels in the ARs, TARs, and DTRs revealed that there were a small number of DEGs between ARs and TARs. On the other hand, the comparison of AR/DTR revealed a high number of DEGs (Figure 3B). This result was in agreement with the root development sequence, i.e., TAR is developmentally close to AR. A total of 6794 unigenes were filtered as DEGs in the three pairwise comparisons (i.e., AR vs. TAR, AR vs. DTR, and TAR vs. DTR). The DEGs identified in the present study (Supplementary Table S3) were further used for gene enrichment analysis and identification of genes involved in tuberous root development.
Figure 3. Venn diagrams of unigenes of three libraries and statistical analysis of the differentially expressed genes (DEGs). (A) Distribution of the unigenes of the three libraries; (B) The red columns indicate the upregulated DEGs and the green columns represent the downregulated DEGs in three pair-wise comparisons (FDR ≤ 0.001 and an absolute value of log2Ratio ≥ 1 was used as the significant threshold for DEGs).
To gain an overall view of the expression profiles of all DEGs across various developmental stages, cluster analysis was performed using the k-means method. Genes were divided into six groups based on their expression patterns (Figure S2). Clusters 1 and 2 comprised genes negatively or positively regulated during root development. Cluster 3 contained genes that were downregulated in the ARs. Cluster 4 was the most abundant group, with 2285 genes that were downregulated in DTRs. Cluster 5 consisted of genes that were downregulated in TARs and then upregulated in DTRs, and Cluster 6 consisted of genes that were upregulated in TARs and then downregulated in DTRs.
Validation of RNA-seq-Based Gene Expression Analysis
We randomly selected 15 differentially expressed unigenes for qRT-PCR validation of our RNA-seq data. The qRT-PCR profiles of 11 genes were completely in agreement with those obtained by using RNA-Seq (Figure 4A), although four differed at one or two stages (Figure 4B). Most of the congruent genes, including unigene5525, unigene66,196, and unigene75,484, were downregulated during tuberous root development (Figure 4A). Unigene5525 encodes a nodulin family protein. Recent advances highlight the importance of nodulin proteins in the transport of nutrients, solutes, amino acids, or hormones, as well as in major aspects of plant development in non-leguminous species (Denancé et al., 2014).
Figure 4. Validation of DEGs by qRT-PCR. (A) Genes displaying congruent results as revealed by qRT-PCR and RNA-seq analyses; (B) Genes showing inconsistent results revealed by qRT-PCR and RNA-seq analyses. The values indicate the means of three biological replicates ± SD. RPKM, reads per kb per million reads.
Pathway Enrichment Analysis for Differentially Expressed Genes
The results of our pathway enrichment analysis suggested that the tuberous root development involves the integration of several pathways. Phenylpropanoid biosynthesis, starch and sucrose metabolism, biosynthesis of plant hormones, and biosynthesis of terpenoids and steroids pathways were significantly affected by the functional transformation of the roots (Supplementary Table S4). Among these, it is notable that phenylpropanoid biosynthesis was downregulated throughout the developmental process. Genes related to lignin biosynthesis in the phenylpropanoid biosynthesis pathway such as cinnamate 4-hydroxylase (C4H), 4-coumarate CoA ligase (4CL), and cinnamyl alcohol dehydrogenase (CAD) were significantly downregulated in DTRs. Similar results were also observed in sweet potato (Firon et al., 2013), which suggested that lignification inhibits tuberous root formaion (Togari, 1950). However, changes in the biosynthesis of plant hormones, as well as that of terpenoids and steroids, were observed during the tuberous root enlargement stage (i.e., these pathways were enriched by DEGs, as indicated by the results of comparative analyses of DTR vs. TAR and DTR vs. AR).
Candidate Genes Related to Tuberous Root Development
Genes Related to Major Carbohydrate Metabolism
We investigated the expression profiles of genes involved in starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, the citrate cycle (TCA cycle), and galactose metabolism to determine the molecular and biochemical changes in the levels of carbohydrates during root development. In higher plants, the entry of carbon from sucrose into cellular metabolism is catalyzed by either sucrose synthase (SUS) or invertase (INV). Interestingly, most of the transcripts encoding SUS were upregulated in DTRs when compared to ARs and TARs. On the other hand, the expression levels of the INV enzymes declined during root development. We did not find any transcript encoding sucrose-phosphatase, which converts sucrose-6-phosphate into sucrose. The transcripts encoding for sucrose-phosphate synthase showed either an upregulated or downregulated expression pattern during tuberous root development. Although two unigenes (unigene80,774 and 1365) that putatively encoded for starch synthase were upregulated in TARs and DTRs relative to that of the ARs, the enzyme ADP-glucose pyrophosphorylase (AGPase, a key enzyme in the starch biosynthetic pathway) was downregulated in DTRs. Moreover, the starch branching enzymes were downregulated during root development, suggesting that starch biosynthesis is not predominant in DTRs (Figure 5A). Most of the transcripts involved in glycolysis/gluconeogenesis, pyruvate metabolism, and TCA cycle were upregulated in ARs and TARs (data not shown), which suggests that added energy, reductants, and carbon building blocks for anabolism are required for root growth at AR and TAR stages.
Figure 5. Heat map of genes in ARs, TARs, and DTRs of R. glutinosa. The expression values were calculated by using the RPKM method and then were log2 transformed prior to the generation of heat maps. (A) The genes related to sucrose and starch metabolism; (B) key genes involved in stachyose metabolism.
According to the KEGG, a total of 249 unigenes were mapped to the galactose metabolism pathway, which features raffinose family oligosaccharide (RFO) metabolism in R. glutinosa. Several unigenes in this category were involved in RFO metabolism, including transcripts encoding stachyose synthase (STS), raffinose synthase, and alkaline alpha-galactosidase, which are the key enzymes responsible for RFO degradation. We found five STS-encoding unigenes and one raffinose synthase-encoding unique in this category, which were enriched in DTRs (Figure 5B). A gene cloning experiment demonstrated that the five unigenes related to stachyose synthase corresponded to one full-length cDNA named RgSTS that is primarily expressed in tuberous roots (Figure 6).
Figure 6. The expression profile of RgSTS as revealed by qRT-PCR analysis. AR, adventitious root; TAR, thickening adventitious root; 60-DAS DTR, developing tuberous root at 60 DAS; 90-DAS DTR, developing tuberous root at 90 DAS; F, flower; St, stem; L, leaf.
Genes Involved in Hormone Metabolism
The endogenous levels of phytohormones (IAA, ABA, ZR, and GA3) were measured in four root stages, which included 30-DAS ARs, 45-DAS TARs, 60-DAS DTRs, and 90-DAS DTRs. The ABA concentration increased during tuberous root development, reaching a maximum at 90 DAS, whereas GA content decreased during root development. The levels of IAA and ZR fluctuated during root development. In general, the highest IAA level was observed in ARs, whereas the ZR level increased at initial stages of root thickening (45 and 60 DAS) and then declined at 90 DAS (Figure 7). Based on the Arabidopsis Hormone Database (http://ahd.cbi.pku.edu.cn) (Jiang et al., 2011) and the KEGG annotation, we identified 82 putative genes related to ABA, GA, IAA, ZR, and ethylene metabolism. Most of these genes were homologs of Arabidopsis genes associated with various hormone actions (Jiang and Guo, 2010). The most significant were the genes crucial for biosynthesis of ABA, GAs, and cytokinins such as 9-cis-epoxycarotenoid dioxygenase, gibberellin 20-oxidase, adenylate isopentenyltransferase, and cytokinin hydroxylase (CYP735A) (Table 1). Besides genes involved in hormone biosynthesis, we also identified several genes responsible for hormone catabolism and conjugation, which function in hormone balance. We did not detect the hormone ethylene in our samples. However, we identified ethylene biosynthetic genes including S-adenosylmethionine synthase, 1-aminocyclopropane-1-carboxylate oxidase, and 1-aminocyclopropane-1-carboxylate synthase (Table 1). The expression patterns of the ethylene biosynthetic genes showed that ethylene content increased during tuberous root development (Table 1).
We identified 2059 TFs through a BLAST search against the Plant Transcription Factor Database 3.0 (http://planttfdb.cbi.pku.edu.cn) (Jin et al., 2014). These TFs accounted for 2.21% (2059/93,172) of the total sequences. We classified the TFs in R. glutinosa into 55 families according to the Plant TFDB category system. Of these, bHLH (8.35%), C3H (6.31%), MYB (5.93%), WRKY (5.73%), and C2H2 (5.39%) were the most abundant types. The distribution of TFs across the genome of R. glutinosa was similar to those of Arabidopsis and rice, except that a higher number of C3H and GRAS family TFs were present in R. glutinosa as compared with that of Arabidopsis and rice (Supplementary Table S5). Among these, a total of 267 were identified as DEGs. The bHLH (34, 12.73%), WRKY (23, 8.61%), GRAS (23, 8.61%), ERF (21, 7.87%), MYB (16, 5.99%), C3H (15, 5.62%), G2-like (15, 5.62%), and HD-ZIP (11, 4.12%) were the top differentially expressed categories; the remaining categories contained <10 differentially expressed TFs (Supplementary Table S5).
We also determined the expression profiles of the differentially expressed TFs. Most transcripts in the C3H, GRAS, and ERF families were upregulated in DTRs compared to that of ARs and TARs, whereas the expression of most G2-like and ARF transcripts was downregulated in DTRs (Supplementary Table S6). Interestingly, the majority of the transcripts encoding ARF proteins were upregulated in ARs and TARs. The expression levels of these transcripts were positively correlated with IAA content, suggesting that ARFs played a role in AR formation (Supplementary Table S6).
A First Insight into Tuberous Root Development in R. glutinosa
Tuberous roots of species such as sweet potato, cassava, R. glutinosa, and Polygonum multiflorum have significant medicinal value, in addition to being an important food resource to humans. Despite their economic importance, research on tuberous roots has not received the same attention as the root systems of model plant species. Lack of extensive genomic sequences and information on functional genes, especially gene expression profiles and functional assignment of genes, has hindered our understanding of the molecular processes of tuberous root development in related species. In the present study, the application of the RNA-Seq technology generated a total of 93,172 transcriptome sequences. These sequences provide a valuable resource for root specific gene discovery and for the assembly of a Rehmannia reference at the genomic level. Our data, which was generated from developing roots, may facilitate in the analysis of transcriptomic changes associated with root development. Moreover, it may help in the genetic manipulation of varieties to control the development of tuberous roots. The global gene expression profiles of ARs, TARs, and DTRs (Supplementary Table S3) revealed the dynamic characteristics of the developing root transcriptome. This RNA-seq based transcriptome analysis provides an initial view of the tuberous root development process, assisting in the identification of candidate genes related to hormone biosynthesis, sugar metabolism, and TFs, thus paving the way for further elucidation of the molecular mechanism underlying tuberous root development.
The Relationship of Sugars and Tuberous Root Development
Glucose and sucrose control various developmental events such as seed germination, flowering, and senescence, and modulate organ formation. Sucrose affects the formation of sink organs and is a notable inducer of tuberization in potato (Perl et al., 1991) and yam (Ovono et al., 2009). Similarly, a high sucrose concentration is required for the in vitro induction of tuberous roots in R. glutinosa (Xue et al., 2002). An exogenous supply of sucrose enhances the production of sweet potato tuberous roots (Tsubone et al., 2000). Sucrose cleavage functions in the allocation of crucial carbon resources, as well as in the initiation of hexose-based sugar signals in importing structures, which have profound developmental effects (Koch, 2004). The channeling of sucrose into a different pathway requires its cleavage. The isoforms of INV and SUS cleave sucrose to channel it into different pathways (Sturm and Tang, 1999). INVs are crucial for sink initiation and expansion, whereas SUSs are mainly responsible for reserve storage and organ maturation (Koch, 2004). The present study has shown that most of the unigenes encoding for SUS were upregulated in DTRs as compared to that of the ARs and TARs; a few were downregulated in DTRs, which may represent the different isoforms of SUS. In contrast, INVs were downregulated during root development. Similar results were observed in sweet potato (Li and Zhang, 2003; McGregor, 2006). The expression profile of INVs and SUSs indicate that these sequentially contribute to the development of tuberous roots. The mechanism underlying how sugars regulate the gene activity of the enzymes related to tuberous root development, including the sucrose cleaving enzymes, remains unclear. Moreover, the crosstalk between sugar and hormone signaling pathways in modulating tuberous root architecture has yet to be determined. Sugars are active effectors of multiple processes, including the biosynthesis and perception of hormones such as ABA (Rolland et al., 2002; Koch, 2004). Recent studies have shown that sugar and hormone signals are tightly coupled to each other during plant growth and development (Gazzarrini and McCourt, 2001; Gibson, 2005; Eveland and Jackson, 2012). The present study identified one transcript that encoded for a sucrose transporter (unigene5098) and another transcript encoding for an INV inhibitor (unigene34,077).
Previous studies have shown that the expression of AGPase and starch-branching enzymes are upregulated in tuberous roots of sweet potato (McGregor, 2006; Wang et al., 2006). However, both enzymes were downregulated during root development in R. glutinosa (Figure 5A). These differences can be attributed to the fact that sweet potato accumulates starch, whereas R. glutinosa stores stachyose in tuberous roots. The metabolism and physiological function of stachyose during tuberous root development are not well-understood. RFOs are transportable carbohydrates found in seeds, tubers, and leaves (Handley et al., 1983; Keller and Matile, 1985; Horbowicz and Obendorf, 1994; Bachmann and Keller, 1995). These are regarded as reserve carbohydrates for seedling germination; these also influence seed viability, as well as contribute to tolerance to various stresses (Bachmann et al., 1994; Horbowicz and Obendorf, 1994; Obendorf, 1997). In the present study, the tuberous root-specific expression of RgSTS suggested that stachyose functions as a reserve carbohydrate that provides energy and the carbon skeleton for plant propagation in favorable conditions, rather than serve as a transportable and osmotic sugar. Stachys sieboldii stores a high amount of stachyose in their tubers, which significantly decreases during the sprouting period (Keller and Matile, 1985), thus suggesting its role in energy storage. Transgenic experiments have demonstrated that reserve starch is not necessary for tuber formation in potato (Müller-Röber et al., 1992); therefore, stachyose might also not be a prerequisite for tuber root formation in Rehmannia. Transgenic Rehmannia lines will help elucidate the physiological role of stachyose in tuberous root formation.
Hormones, Hormone-Related Genes, and Tuberous Root Development
Our experiments showed that ABA levels increased during tuberous root development. Furthermore, the expression patterns of most unigenes encoding 9-cis-epoxycarotenoid dioxygenase (NCED) and all the sequences encoding abscisic-aldehyde oxidase (AAO) were positively correlated with changes in the level of ABA during root development. The endogenous levels of ABA show a positive correlation with the thickening potential of sweet potato tuberous roots (Nakatani et al., 1988). Therefore, ABA might play a crucial role in tuberous root enlargement and/or the sink potential establishment. Recent studies (Arenas-Huertero et al., 2000; Cheng et al., 2002; Travaglia et al., 2007) have implicated ABA in plant growth and carbohydrate accumulation, in addition to its function as a stress hormone. Endogenous ABA levels promote plant growth and development in the absence of severe stress (Arenas-Huertero et al., 2000; Cheng et al., 2002). Exogenous application of ABA to wheat leaves increases yield by accumulating carbohydrates and redistributing grains (Travaglia et al., 2007). Therefore, ABA plays a complex physiological role in plant growth and development.
In sweet potato, high cytokinin and auxin levels promote the onset and subsequent primary thickening and growth of tuberous roots (Matsuo et al., 1983, 1988; Nakatani and Komeichi, 1991, 1992). The present study showed that the level of endogenous ZR increased in 45-DAS TARs and 60-DAS DTRs, indicating that ZR promotes tuberous root formation and its initial enlargement. The expression levels of adenylate isopentenyltransferases (IPTs) and cytokinin hydroxylases (CYP735As) increased in 45-DAS TARs, but dramatically declined in 60-DAS DTRs, which was not in agreement with the ZR contents. Auxin has been shown to control AR formation (Geiss et al., 2009; Li et al., 2009b). Hormone content analysis conducted in the present study has confirmed the role of auxin in AR formation, as AR has the highest IAA content in all the tested roots. Previous studies have implicated multiple pathways in the de novo biosynthesis of auxins; however, the details have yet to be elucidated (Mano and Nemoto, 2012). We identified transcripts encoding the enzymes tryptophan N-monooxygenase and indole-3-pyruvate monooxygenase, which are involved in auxin biosynthesis. However, their expression patterns did not coincide with the IAA profiles. Both genes showed a high expression level in ARs compared to that in DTRs. In the present study, ZR and IAA contents were not consistent with the levels of their biosynthetic genes, implicating the complexity of their regulation. GA3 levels declined during root development and were corroborated by the expression studies of gibberellin 20-oxidase and gibberellin 3-beta-dioxygenase, which also declined during the process. The effects of GA on tuberous root have not been fully determined. GA promotes stolon elongation and inhibits tuber formation in potato (Smith and Rappaport, 1960; Kumar and Wareing, 1972; Railton and Wareing, 1973; Krauss and Marschner, 1982). Tubers and tuberous roots are of different origins; however, these may recruit common physiological factors such as hormones to trigger downstream actions when initiating sink development. Several GRAS-domain containing TFs were upregulated during root development (Supplementary Table S6). GRAS proteins such as DELLA and SCARECROW are essential components of the GA signaling pathway and function in multiple plant process, including root development (Bolle, 2004; Hirsch and Oldroyd, 2009). Ethylene biosynthetic genes and ERFs were also upregulated during tuberous root development, suggesting a role for ethylene in the process. Similar results were reported in sweet potato tuberous roots (Firon et al., 2013). ERFs control cell dedifferentiation, metabolite biosynthesis, and trait development (Song et al., 2013). The findings of the present study may facilitate in better understanding the role of ethylene in tuberous root development.
We also identified various genes involved in hormone homeostasis, which include abscisic acid 8′-hydroxylase, gibberellin 2-beta-dioxygenase, GH3, cytokinin hydroxylase, cytokinin dehydrogenase, and UDP-glycosyltransferase 76C2. Among these, GH3 catalyzes the synthesis of IAA-amino acid conjugates, thus providing a mechanism for the plant to cope with the presence of excess auxin (Staswick et al., 2005). A recent report demonstrated that auxin-inducible proteins GH3.3, GH3.5, and GH3.6, which can conjugate JA, do not act by regulating auxin levels, but control JA homeostasis in Arabidopsis. Thus, auxin controls AR initiation through a complex network that requires GH3, which regulates JA homeostasis that ultimately controls AR formation (Gutierrez et al., 2012). The Arabidopsis study (Gutierrez et al., 2012) can partially explain the abundance of GH3 genes in ARs, which have high IAA content. Auxin promotes root growth by modulating the gibberellin response (Fu and Harberd, 2003). Moreover, ethylene (perhaps via an auxin-dependent pathway) and auxin both alter the abundance/stability of DELLA (Achard et al., 2003). Hormones also interact with other signals to regulate plant development. Sugars and ABA tend to act synergistically during embryo growth, in transitioning from a phase of rapid cell division to cell enlargement, and in the accumulation of storage reserves (Eveland and Jackson, 2012). Therefore, a complex, dynamic, and possibly interconnected signaling network is involved in tuberous root development. Our results may help determine the role of hormones in tuberous root development at the molecular level.
Overview of the Biological Processes Related to Tuberous Root Development
Tuberous root formation and enlargement are of two vital important processes of Rehmannia tuberous root development. However, the root development programs in R. glutinosa are highly complex and need to be finely controlled and coordinated by the intervention of several crosstalks. Based on the data presented in the present study, we propose a preliminary overview of the important biological processes occurring during Rehmannia root development, which are in part schematically represented in Figure S3. In favorable conditions, including environment factors such as light, temperature, humidity, and nutrient status, the plants detect these signals and initiate various developmental processes such as root development, which involves hormone regulation (including hormone biosynthesis, conjugation, transport, and degradation), carbohydrate metabolism, and upregulation of a series of root development-related genes. In general, at the AR formation stage, the level of the plant hormone auxin significantly increases in ARs, which has been shown to be intimately involved in the process of adventitious rooting (Geiss et al., 2009; Li et al., 2009b). Crosstalks between auxin and the plant hormones GA and JA finely tune the adventitious rooting process (Fu and Harberd, 2003; Gutierrez et al., 2012). Once the carbohydrate supply in the whole plant is sufficient, ARs transform into tuberous roots under high cytokinin concentrations, which has been implicated in the onset of tuberous root in sweet potato (Matsuo et al., 1983, 1988; Nakatani and Komeichi, 1991). Then, the tuberous roots become the largest sinks; the plants trigger a series of biological processes such as ABA and ethylene biosynthesis, conversion of sugars, which act as developmental signals to regulate the transport and metabolism of carbohydrates, and the transcription of developmental genes in DTRs. As earlier described, certain crosstalks between the hormone and sugar signaling pathways might be involved in the regulation of tuberous root development. Tuberous root development is accompanied by the downregulation of genes in the lignin biosynthesis pathway. The degree of lignification in roots might be a limiting factor for tuberous root formation and enlargement, as these were also observed in sweet potato (Togari, 1950; Firon et al., 2013).
In the present study, we generated the root transcriptome of R. glutinosa. We compared the expression profiles of ARs, TARs, and DTRs to discover a series of genes involved in tuberous root development. Our study provides valuable insights into the molecular basis of tuberous root development and presents critical information for breeding and engineering that may yield large amount of Rehmannia varieties. We also showed that tuberous root development is a complicated process that involves the interaction of several genes and various physiological factors. Our large transcriptome dataset also provides information that could be employed in comparative analysis studies of closely related species.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the National Natural Science Foundation of China (Grant No. 30901969 and 81373915) and the Ph.D. Programs Founded by the Ministry of Education of China (20091106120047).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.00396/abstract
Achard, P., Vriezen, W. H., van der Straeten, D., and Harberd, N. P. (2003). Ethylene regulates Arabidopsis development via the modulation of DELLA protein growth repressor function. Plant Cell 15, 2816–2825. doi: 10.1105/tpc.015685
Arenas-Huertero, F., Arroyo, A., Zhou, L., Sheen, J., and León, P. (2000). Analysis of Arabidopsis glucose insensitive mutants, gin5 and gin6, reveals a central role of the plant hormone ABA in the regulation of plant vegetative development by sugar. Genes Dev. 14, 2085–2096. doi: 10.1101/gad.14.16.2085
Bachmann, M., Matile, P., and Keller, F. (1994). Metabolism of the raffinose family oligosaccharides in leaves of Ajuga reptans L. cold acclimation, translocation, and sink to source transition: discovery of a chain elongation enzyme. Plant Physiol. 105, 1335–1345.
Cheng, W. H., Endo, A., Zhou, L., Penney, J., Chen, H. C., Arroyo, A., et al. (2002). A unique short-chain dehydrogenase/reductase in Arabidopsis glucose signaling and abscisic acid biosynthesis and functions. Plant Cell 14, 2723–2743. doi: 10.1105/tpc.006494
Firon, N., LaBonte, D., Villordon, A., Kfir, Y., Solis, J., and Lapis, E., et al. (2013). Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation. BMC Genomics 14:460. doi: 10.1186/1471-2164-14-460
Gutierrez, L., Mongelard, G., Floková, K., Pacurar, D. I., Novák, O., Staswick, P., et al. (2012). Auxin controls Arabidopsis adventitious root initiation by regulating jasmonic acid homeostasis. Plant Cell 24, 2515–2527. doi: 10.1105/tpc.112.099119
Handley, L. W., Pharr, D. M., and McFeeters, R. F. (1983). Carbohydrate changes during maturation of cucumber fruit: implications for sugar metabolism and transport. Plant Physiol. 72, 498–502. doi: 10.1104/pp.72.2.498
Horbowicz, M., and Obendorf, R. L. (1994). Seed desiccation tolerance and storability: dependence on flatulence-producing oligosaccharides and cyclitols—review and survey. Seed Sci. Res. 4, 385–405. doi: 10.1017/S0960258500002440
Jiang, Z., Liu, X., Peng, Z., Wan, Y., Ji, Y., and He, W. (2011). HD2.0: an update version of Arabidopsis Hormone Database for plant systematic studies. Nucleic Acids Res. 39, D1123-D1129. doi: 10.1093/nar/gkq1066
Jin, J., Zhang, H., Kong, L., Gao, G., and Luo, J. (2014). PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 42, D1182-D1187. doi: 10.1093/nar/gkt1016
Krauss, A., and Marschner, H. (1982). Influence of nitrogen nutrition, daylength, and temperature on contents of gibberellic and abscisic acid and on tuberization in potato plants. Potato Res. 25, 13–21. doi: 10.1007/BF02357269
Li, J. J., Fan, H. J., Yao, H. L., Sun, H., Zhou, Y. Q., and Zhang, Z. Y. (2009a). Comparative studies on the yields and qualities of differentcultivars of Rehmannia glutinosa f. hueichingensis under plastic film mulching cultivation. J. Henan Agri. Univ. 43, 252–255.
Li, J. J., Fan, X. M., Sun, H., Gu, F. P., Zhang, C. Y., Gao, Z. M., et al. (2007). The correlation analyses of yield related six popularly cultivated cultivars of Rehmannia glutinosa f. hueichingensis. J. Henan Agric. Sci. 9, 86–96.
Li, M. J., Yang, Y. H., Chen, X. J., Wang, F. Q., Lin, W. X., and Yi, Y. J., et al. (2013). Transcriptome/degradome-wide identification of R. glutinosa miRNAs and their targets: the role of miRNA activity in the replanting disease. PLoS ONE 8:e68531. doi: 10.1371/journal.pone.0068531
Li, R., Zhu, H., Ruan, J., Qian, W., Fang, X., Shi, Z., et al. (2010). De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 20, 265–272. doi: 10.1101/gr.097261.109
Li, X. Q., and Zhang, D. (2003). Gene expression activity and pathway selection for sucrose metabolism in developing storage root of sweet potato. Plant Cell Physiol. 44, 630–636. doi: 10.1093/pcp/pcg080
Matsuo, T., Mitsuzono, H., Okada, R., and Itoo, S. (1988). Variations in the levels of major free cytokinins and free abscisic acid during tuber development of sweet potato. J. Plant Growth Regul. 7, 249–258. doi: 10.1007/BF02025267
Matsuo, T., Yonedo, T., and Itoo, S. (1983). Identification of free cytokinins as the changes in endogenous levels during tuber development of sweet potato (Ipomoea batatas, Lam.). Plant Cell Physiol. 24, 1305–1312.
McGregor, C. (2006). Differential Expression and Detection of Transcripts in Sweetpotato (Ipomoea batatas (L.) LAM.) using cDNA Microarrays. PhD thesis, Louisiana State University and Agricultural and Mechanical College.
Müller-Röber, B., Sonnewald, U., and Willmitzer, L. (1992). Inhibition of the ADP-glucose pyrophosphorylase in transgenic potatoes leads to sugar-storing tubers and influences tuber formation and expression of tuber storage protein genes. EMBO J. 11, 1229–1238.
Nakatani, M., and Komeichi, M. (1991). Changes in the endogenous level of zeatin riboside, abscisic acid and indole acetic acid during formation and thickening of tuberous roots in sweetpotato. Jpn. J. Crop Sci. 60, 91–100. doi: 10.1626/jcs.60.91
Nakatani, M., Komeichi, M., and Watanabe, Y. (1988). Tuber sink potential in sweetpotato (Ipomoea batatas Lam.). II. Estimation of tuber sink potential of cultivars using single leaf grafts. Jpn. J. Crop Sci. 57, 544–552. doi: 10.1626/jcs.57.544
Ovono, P. O., Kevers, C., and Dommes, J. (2009). Effects of reducing sugar concentration on in vitro tuber formation and sprouting in yam (Dioscorea cayenensis-D. rotundata complex). Plant Cell Tissue Organ Cult. 99, 55–59. doi: 10.1007/s11240-009-9575-1
Perl, A., Aviv, D., Willmitzer, L., and Galun, E. (1991). In vitro tuberization in transgenic potatoes harboring β-glucuronidase linked to a patatin promoter: effects of sucrose levels and photoperiods. Plant Sci. 73, 87–95. doi: 10.1016/0168-9452(91)90129-V
Schafleitner, R., Tincopa, L. R., Palomino, O., Rossel, G., Robles, R. F., and Alagon, R., et al. (2010). A sweetpotato gene index established by de novo assembly of pyrosequencing and Sanger sequences and mining for gene-based microsatellite markers. BMC Genomics 11:604. doi: 10.1186/1471-2164-11-604
Smith, O. E., and Rappaport, L. (1960). “Endogenous gibberellins in resting and sprouting potato tubers,” in Advances in Chemistry, ed R. F. Gould (Washington, DC: American Chemical Society press), 42–48.
Song, X., Li, Y., and Hou, X. (2013). Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genomics 14:573. doi: 10.1186/1471-2164-14-573
Staswick, P. E., Serban, B., Rowe, M., Tiryaki, I., and Maldonado, M. T., Maldonado, M. C., et al. (2005). Characterization of an Arabidopsis enzyme family that conjugates amino acids to indole-3-acetic acid. Plant Cell 17, 616–627. doi: 10.1105/tpc.104.026690
Sun, P., Guo, Y. H., Qi, J. J., Zhou, L. L., and Li, X. E. (2010). Isolation and expression analysis of tuberous root development related genes in Rehmannia glutinosa. Mol. Biol. Rep. 37, 1069–1079. doi: 10.1007/s11033-009-9834-6
Sun, P., Song, S. H., Zhou, L. L., Zhang, B., Qi, J. J., and Li, X. E. (2012). Transcriptome analysis reveals putative genes involved in iridoid biosynthesis in Rehmannia glutinosa. Int. J. Mol. Sci. 13, 13748–13763. doi: 10.3390/ijms131013748
Tao, X., Gu, Y. H., Wang, H. Y., Zheng, W., Li, X., and Zhao, C W. et al. (2012). Digital gene expression analysis based on integrated de novo transcriptome assembly of sweetpotato [Ipomoea batatas (L.) Lam.]. PLoS ONE 7:e36234. doi: 10.1371/journal.pone.0036234
Teng, N., Wang, J., Chen, T., Wu, X., Wang, Y., and Lin, J. (2006). Elevated CO2 induces physiological, biochemical and structural changes in leaves of Arabidopsis thaliana. New Phytol. 172, 92–103. doi: 10.1111/j.1469-8137.2006.01818.x
Travaglia, C., Cohen, A. C., Reinoso, H., Castillo, C., and Bottini, R. (2007). Exogenous abscisic acid increases carbohydrate accumulation and redistribution to the grains in wheat grown under field conditions of soil water restriction. J. Plant Growth Regul. 26, 285–289. doi: 10.1007/s00344-007-9018-3
Tsubone, M., Kubota, F., Saitou, K., and Kadowaki, M. (2000). Enhancement of tuberous root production and adenosine-5w-diphosphate pyrophosphorylase (AGPase) activity in sweet potato (Ipomoea batatas Lam.) by exogenous injection of sucrose solution. J. Agron Crop Sci. 84, 181–186. doi: 10.1046/j.1439-037x.2000.00396.x
Wang, S. J., Chen, M. H., Yeh, K. W., and Tsai, C. Y. (2006). Changes in carbohydrate content and gene expression during tuberous toot development of sweet potato. J. plant Biochem. Biotechnol. 15, 21–25. doi: 10.1007/BF03321896
Wang, Z., Fang, B., Chen, J., Zhang, X., Luo, Z., and Huang, L., et al. (2010). De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato (Ipomoea batatas). BMC Genomics 11:726. doi: 10.1186/1471-2164-11-726
Wen, X. S., Li, X. E., and Yang, S. L. (2001). Viral diseases of Rehmannia glutinosa and problems demanding prompt solution. Chin. Trad. Herb. Drug. 7, 23–24. doi: 10.3321/j.issn:0253-2670.2001.07.046
Yang, Y., Chen, X., Chen, J., Xu, H., Li, J., and Zhang, Z. (2011). Differential miRNA expression in Rehmannia glutinosa plants subjected to continuous cropping. BMC Plant Biol. 11:53. doi: 10.1186/1471-2229-11-53
Yang, Y. H., Li, M. J., Li, X. Y., Chen, X. J., Lin, W. X., and Zhang, Z. Y. (2014). Transcriptome-wide identification of the genes responding to replanting disease in Rehmannia glutinosa L. roots. Mol. Biol. Rep. 42, 881–892. doi: 10.1007/s11033-014-3825-y
Keywords: Rehmannia glutinosa, tuberous root, development, transcriptome, RNA-seq
Citation: Sun P, Xiao X, Duan L, Guo Y, Qi J, Liao D, Zhao C, Liu Y, Zhou L and Li X (2015) Dynamic transcriptional profiling provides insights into tuberous root development in Rehmannia glutinosa. Front. Plant Sci. 6:396. doi: 10.3389/fpls.2015.00396
Received: 15 February 2015; Accepted: 18 May 2015;
Published: 10 June 2015.
Edited by:Tiegang Lu, Chinese Academy of Agricultural Sciences, China
Reviewed by:Jingfa Xiao, Chinese Academy of Sciences, China
Rita Sharma, Jawaharlal Nehru University, India
Copyright © 2015 Sun, Xiao, Duan, Guo, Qi, Liao, Zhao, Liu, Zhou and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xingguo Xiao, State Key Laboratory of Plant Physiology and Biochemistry, College of Biological Sciences, China Agricultural University, No. 2 Yuanmingyuan West Road, Beijing 10093, China, firstname.lastname@example.org;
Xianen Li, Center for Medicinal Plant Cultivation, Institute of Medicinal Plant Development, Peking Union Medical College, Chinese Academy of Medical Sciences, No. 151 Malianwa North Road, Beijing 100193, China, email@example.com