Transcriptome Assembly and Gene Expression Analysis of Leaf and Rhizome tissues of Dioscorea opposita

Transcriptome Assembly and Gene Expression Analysis of Leaf and Rhizome tissues of Dioscorea opposita | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Transcriptome Assembly and Gene Expression Analysis of Leaf and Rhizome tissues of Dioscorea opposita Dianyun Hou, Ruiyang Cheng, Zhihong Li, Jingya Xu, Zhanqiang Ma, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-28829/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: Dioscorea opposita is a monocotyledonous plant in the Dioscoreaceae family. The rhizome of D. opposita , known as “dioscoreae rhiaoma” or “Chinese yam”, is not only a well-known traditional Chinese medicine, but also a popular health food. Methods: In this study, the leaf and rhizome tissues of D. opposita were collected. Their transcriptome sequenceand were obtained using an Illumina Hiseq 4000 platform. The data was analyzed by bioinformatics software. Likewise, the consistent expression pattern of DEGs was validated with qRT-PCR analysis. Results: A total of 46,982,031 and 53,689,620 clean reads from leaf and rhizome tissues were acquired, respectively (SRP: PRJNA628588). 50,765 unigenes were obtained with an average length of 1050 base pairs (bp). A total of 13,259 significant differentially expressed genes (DEGs) were obtained from leaf versus rhizome tissues, including 5536 up-regulated and 7723 down-regulated unigenes. In addition, 20,607 simple sequence repeats (SSRs) were identified in this study. Conclusions: This study not only enriched the transcriptomic data of D. opposita , but also laid a foundation for the analysis of terpenoid synthesis pathway. Integrative & Complementary Medicine Dioscorea opposita leaf and rhizome transcriptome assembly expression terpene biosynthesis Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Background Dioscorea opposita Thunb (Family Dioscoreaceae), is an important staple in tropical regions and edible crop mainly distributed in parts of Africa, South America, and Asia [ 1 – 3 ]. The rhizome of D. opposita , known as “dioscoreae rhiaoma” or “Chinese yam”, has been listed as both a nutritious food source and an invigorant in traditional Chinese medicine for thousands of years [ 4 – 6 ]. Various types of natural constituents, including dopamine, starch, dioscin, polysaccharides, mucin, allantoin, choline, phytosterols, and essential amino acids, have been reported to be in D. opposita [ 7 – 9 ]. There were reports that 2-C-Methyl-D-Erythritol-4-Phosphate (MEP) and mevalonate acid (MVA) pathway are the early biosynthetic pathways of dioscin in many species [ 10 – 13 ].It is well known that 1-deoxy-D-xylulose-5-phosphate synthase (DXS), 1-deoxy-D-xylulose 5-phosphate reductoisomerase (DXR), geranylgeranyl pyrophosphate synthase (GGPS) and hydroxymethylglutaryl-CoA reductase (HMGR) are the key enzyme genes, related to the biosynthesis of dioscin[ 13 ]. However, there are no reports on key enzyme genes, molecular biology studies or comprehensive genomic and transcriptomic data of D. opposita . RNA-Seq is an efficient method for transcriptome analysis, which can also provide a highly sensitive and accurate platform at a low cost and is efficient for gene discovery and biosynthesis pathway analysis in non-model plants without knowing their genome sequences [ 14 – 17 ]. RNA-Seq is currently a very popular method that is applied in the exploration of plants at transcriptomic levels. The ultimate aim is often to reveal genes that are involved in the biosynthetic pathways in secondary metabolism of many species [ 18 – 21 ]. Until now, the genomic resources of D . alata [ 22 ], D. rotundata [ 23 ], and D. dumetorum [ 24 ] have been reported. Nevertheless, as a non-model species, there is lack of enomic resources, in particular, information on active substance metabolism of D. opposita . In the study, the sample of leaf and rhizome tissues of D. opposita were collected from Wenxian County, called “Tie Gui Shan Yao”. It is a local cultivar of D. opposita , belonging to one of the “four huaiqing Chinese medicines”. There have reported that the content of terpenoids in“Tie Gun Shan Yao” was higher than other D. opposita [ 25 ]. In this study, transcriptome analysis was performed from leaf and rhizome tissues of D. opposita . The results will enrich the known genetic information for further gene function research and analyses, such as differential transcript abundance analysis. We also try to understand the transcript abundance difference and explore the molecular mechanisms of biosynthetic pathways in secondary metabolites in D. opposita . Materials And Methods Plant materials and RNA extraction Leaf and rhizome tissues of D. opposita were collected from Wenxian County, Jiaozuo City, Henan Province of China, on October 2017（Wenxian County Weather Station： N34°56′49.54″, E113°10′18.54″; T:15℃）. Three biological replicates from three leaf tissues (HSYY_1, HSYY_2, and HSYY_3) and three rhizome tissues (HSYK_1, HSYK_2, and HSYK_3) were used for the analysis. All samples were flash frozen in liquid nitrogen after cleaning and stored at –80 °C until RNA extraction. Total RNAs were extracted from the leaf and rhizome tissues using the TRIzol® Reagent (Invitrogen, USA) following the product description with slight improvement [26, 27]. The extracted RNAs were checked through NanoDrop™2000 spectrophotometer (NanoDrop Technologies, USA) and RNA integrity was assessed using electrophoresis in 1.0% agarose gel. Furthermore, the RNA concentration was determined by the RNA Agilent Bioanalyzer 2100 system (Agilent, CA, USA) to guarantee the RNA Integrity Number (RIN) values were not less than 7.0. Transcriptome sequencing and de novo assembly The cDNA libraries were constructed using the TruSeq-TM RNA sample preparation Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. The cDNA libraries were enriched through a polymerase chain reaction (PCR) with 15 cycles. Then, the target band was extracted and recycled by a 2% low range ultra-agarose gel, and quantified by TBS380 Mini-Fluorometer. The suitable cDNA libraries were sent to Shanghai Majorbio Biopharm Biotechnology Co., Ltd. (Shanghai, China) where the RNA-Seq was performed using the Illumina HiSeq4000 SBS Kit (300 cycles, Illumina, San Diego, CA, USA) on an Illumina HiSeq4000 platform. The raw reads were cleaned and assessed following the steps. Firstly, the adaptor sequence of reads was removed, and discarded the reads that have no inserted fragments due to the reason of adaptor self-connection; Secondly, low quality bases (quality value less than 20) were trimed off at the end of the sequence (3 '). If there are still bases with quality value less than 10 in the remaining reads, the whole reads are eliminated, otherwise it is retained.Thirdly, reads with “N” bases was also removed. Lastly, the reads less than 100bp in length were discarded after removing adapter and quality pruning. All the above mentioned processes were performed carried out using the SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (Version 1.33) software [28]. The high quality clean reads of leaves and rhizomes of D. opposita were de novo assembled by Trinity (Version v2.8.5) software with operation parameters of k-mer size of 25, 27, 30, 31and 32(max normalize read 50，Min contig length 200). The assembly results were assessed based on the assembly size, number of contigs, N50 and complete BUSCOs. Then, 2x150bp reads were generated and filtered, and more than 200 bp contigs were produced. At last, assembly completeness and quality also were assessed with Benchmarking Universal Single-Copy Orthologs (BUSCO, version 3.0.2) and TransRate package (version 1.0.3) [29, 30]. Functional annotation and classification of unigenes The unigenes were compared with the NCBI BlastX program (E-value <1 E-5) against the non-redundant (NR), KEGG (Kyoto Encyclopedia of Genes and Genomes), KOG (clusters of orthologous groups for eukaryotic complete genomes), GO (gene ontology) and Pfam [31-33]. The functional classification of these unigenes was performed using Blast2GO and KEGG. The GO analysis was categorized into biological processes, molecular functions, and cellular component [34]. The pathway distributions were determined based on the KEGG database [35]. Unigenes were also aligned to the Evolutionary Genealogy of Genes: Non-Supervised Orthologous Groups (eggNOG) database for further functional prediction and classification [36]. Differential transcript abundance analysis The clean reads were used to determine alignment counts and to quantify transcript abundance by the RSEM (RNA-seq by Expectation Maximizattion, http://deweylab.github.io/RSEM/) package [37]. The unigene and transcript abundance levels were analyzed by TPM (transcripts per million reads) [38]. The software DESeq2 (Version 1.24.0, http://bioconductor.org/packages/ stats/bioc/DESeq2/) was used to analyze the differentially expressed genes (DEGs) between leaf and rhizome tissues of D. opposita . The false discovery rate (FDR) can be used to assess p -value thresholds in multiple testing [39, 40]. The significance of DEGs was defined with a threshold of FDR < 0.05 and an absolute value of log2-fold change of ≥2.0 (|log2fc|≥1.0) [41]. The DEGs enrichment analyses were carried out for GO and KEGG terms by the GOatools and KOBAS software [42, 43]. SSR Mining and analysis The unigenes were searched for the prediction and detection of simple sequence repeats (SSR) through the MISA 1.0 software with default parameters [44]. The motifs of mono-, di-, tri-, tetra-, penta-, and hexa- nucleotide repeats were identified. A minimum of five repetitions were considered as the search criteria in MISA. The parameters of definement of microsatellites (unit size/minimum number of repeats) were used to predict SSRs. That is Mono- (1/10), Di- (2/6), Tri-(3/5), Tetra- (4/5), Penta- (5/5) and Hexa-nucleotide (6/5). And the maximal number of bases interrupting 2 SSRs in a compound microsatellite was 100. Expression analysis of terpene biosynthetic related DEGs We selected 2 DEGs, DN39197_c0_g1 and DN37057_c0_g1, which were the key enzyme genes in the terpene synthesis from D. opposita , to validate their expression using quantitative real-time PCR (qRT-PCR). Total RNAs from leaf and rhizome tissues of D. opposita were extracted through the same procedures as used for RNA-seq in the study. The RNA from each sample was reverse transcribed using a reverse transcription kit (PC17-TRUE script, Aidlab Biotechnologies Co., Ltd, China). The DoActin gene was used as an internal control. Primers for the 2 DEGs and DoActin were designed in the Primer Premier 5.0 software. All the primers sequences are listed in Table 1 . The reaction system (10 µL) included 5 µL of SYBR Premix Ex TaqTM, 4 µL of cDNA template, and 0.5 µL of each primer (10 µM). The qRT-PCR reactions were conducted on a Roche Light Cycle 96 real-time PCR system (Roche, USA) with SYBR green fluorescent dye (TaKaRa, China). PCR amplification was initiated at 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s, 60 °C for 30 s, and 72 °C for 30 s. Three replicates were set for the experiment. Relative transcript abundance levels were calculated by the 2-ΔΔCt method [ 45 ]. Table 1 Primers for qRT-PCR of 2 DEGs and DoActin Unigene ID Annotation Primer sequence(5’ − 3’) DN39197_c0_g1 DXR DXR-F: GGTGCCGATTACGACCATATC DXR-R: GGATTGGCAGACGCATATCA DN37057_c0_g1 GGPS GGPS-F: ATCGGTGGTGGTTCTGATTC GGPS-R: CTGCAGTCTTCCCAAGTTCT Actin Actin-F: CTAAGGCCAACAGAGAGAAGATG Actin-R: GGTACGGCCACTTGCATAAA Results Transcriptome sequencing and de novo assembly We obtained averages of 48,786,977 and 55,202,249 raw reads for the rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues and leaf (HSYY_1, HSYY_2, and HSYY_3) tissues, respectively, which consisted of 7,366,833,477 bp and 8,335,539,549 bp. We obtained 46,982,031 and 53,689,620 clean reads, including 7,053,884,682 bp and 8,067,994,222 bp clean bases, respectively. The sequences can be found in NCBI SRA with the accession number PRJNA628588. Moreover, the average GC content of rhizome and leaf tissues was 46.8% and 46.7%, respectively. The k-mer size of 25, 27, 30, 31 and 32 were selected to evaluate the assembly results based on the assembly size, number of contigs, N50 and complete BUSCOs. The results showed that the k-mer size of 30 for the final assembly results in the best final assembly. Trinity assembly resulted in 50,765 unigenes with fragment size ranging from 201 bp to 16,991 bp, and the average length of each unigene was 1049.6 bp. In addition, all of the unigenes consisted of a total of 53,283,125 bp. The unigenes were retained with an N50 length of 1825 bp and N90 length of 2797 bp. The TransRate assembly and Benchmarking UniversalSingle-Copy Orthologs (BUSCO) completeness assessment score were 0.26405 and 68.9%, respectively. The detailed sequencing and assembly results for the leaf and rhizome tissues of D. opposita are shown in Table 2 . Table 2 Optimized assembly results of unigenes from D. opposita rhizome and leaf tissues Type Unigene Total sequence num 50,765 Total sequence base (bp) 53,283,125 Largest (bp) 16,991 Smallest (bp) 201 Average length(bp) 1,050 N50 (bp) 1,825 N90 (bp) 2,797 TransRate score 0.26405 Complete single copy BUSCOs 68.9% Functional annotation of D. opposita unigenes The functional annotation results indicated that a total of 24,840 unigenes were annotated (Table 3 ). Among them, different unigenes were matched in different databases, for example, 23,848 (46.98%) unigenes were matched in the NR database.In addition, 18,531 (36.50%) unigenes were matched in the Swiss-Prot database. Table 3 The annotation statistics of D. opposita Type Unigene number (percent%) NR 23,848(46.98) Swiss-Prot 18,531(36.50) Pfam 18,805(37.04) COG 21,973(43.28) GO 20,014(39.42) KEGG 10,617(20.91) Total_annot 24,840(48.93) Total 50,765(1) The length distributions of annotated and unannotated unigenes results indicated that most annotated and unannotated unigenes were distrbuted in the 200–499 bp range (annotated unigenes: 5970/24,840; unannotated unigenes: 16,040/25,925), followed by the 500–799 bp ranges (annotated unigenes: 3446/24,840; unannotated unigenes:5566/25,925), and 800–1099 bp range (annotated unigenes: 2394/24,840; unannotated unigenes: 2194/25,925). Long sequences of unigens were clearly more likely to be annotated than those with short sequences. A total of 7151 unigenes measuring over 2000 bp in length were annotated, and only 371 unigenes were unannotated. Functional classification results of unigenes To further evaluate the annotated unigene data, the unigene sequences were searched against the cluster of orthologous groups of proteins (COG) databases. In total, 23,139 unigenes assigned to 23 COG categories under the following types: information storage and processing, cellular processes and signaling, metabolism, and poorly characterized (Fig. 1). Among them, alphabetized“Posttranslational modification, protein turnover, chaperones”(O) (1373; 5.93 %) was the categories with the most members in the COG databases. Using the Blast2GO database, the 20,014 unigenes in the gene ontology (GO) analysis could be summarized into three main GO categories biological process (BP), molecular function (MF), and cellular component (CC) on two levels. In the BP category, the cellular process (7345;36.70%), metabolic process (6480; 32.389%), and biological regulation (2125; 10.62%) were represented the most. Likewise, cell part (7756; 38.75%), membrane part (7018; 35.07%),organelle (4292; 21.44%) and organelle part (3154; 15.76%)were the most abundant in the CC category. Furthermore, the MF category mainly consisted of binding (10,949; 54.71%), catalytic activity (10,794; 53.93%), transporter activity (1074; 5.37%), and structural molecule activity (730; 3.65%) (Fig. S1). In order to perform an in-depth analysis of unigene function based on the GO database,all of the three categories (BP, MF, and CC) were subdivided into level 3 and level 4 terms in order to determine their function. For example, cellular process is the main level 2 term from the BP category. In the level 3 terms, cellular process was assigned to the following specific functions: cellular metabolic process (GO: 0044237; 5615 unigenes; 28.06%), cellular component organization (GO:0016043; 1147 unigenes; 5.73%), cellular response to stimulus (GO:0051716; 555 unigenes; 2.78%), and signal transduction (GO:0007165; 526 unigenes; 2.63%), etc. (Fig. S2). Through the KEGG database BLAST, 10,617 unigene sequences were assigned to 145 pathways (Table S1). Among them, the three pathways with the greatest number of unigene sequences were ribosome (map03010; 605 unigenes; 5.70%), plant-pathogen interaction (map04626; 281 unigenes; 2.65%),protein processing in endoplasmic reticulum (map04141; 276 unigenes; 2.60%), and thermogenesis (map04714; 248 unigenes; 2.34%) (Table S1). In addition, all the pathways were divided into five main types of KEGG metabolic pathways (Fig. 2). Expression analysis of unigenes from leaf and rhizome tissues In order to explore the differential expression in leaf and rhizome tissues of D. opposit a , the clean reads from the leaf and rhizome library were compared with the unigene database by using the RNA-seq by RSEM package. The results indicated that 41,692,444 reads, respectively 77.62% of the original clean reads (53,689,620), were aligned in leaf tissues (HSYY: average reads of HSYY_1, HSYY_2, and HSYY_3). Likewise, 37,978,278 reads, representing 80.83% of the original clean reads (46,982.031), were aligned in rhizome tissues (HSYK: average reads of HSYK_1, HSYK_2, and HSYK_3). In this study, the sample correlation analysis was performed using RSEM software (Tables S1). A high correlation coefficient was obtained among the three biological replicates of each tissue. The correlation coefficients of rhizome_1 (HSYK_1) with rhizome_2 (HSYK_2) and rhizome_3 (HSYK_3) were 0.935 and 0.913, respectively. The correlation coefficients of leaf_1 (HSYY_1) with leaf_2 (HSYY_2) and leaf_3 (HSYY_3) were 0.919 and 0.935, respectively. The correlation coefficient between two different of leaf and rhizome tissues was significantly lower than that among the three biological replicates of each leaf and rhizome tissue. The results indicated that all three biological replicates of each tissue are all clustered together, thus suggesting that the sequencing data are suitable for further biological analysis (Fig. S3). Venn analysis of samples with a unigene expression greater than 2 was performed using RSEM software (Fig. S4). A total of 14,627 unigenes were shared among the six samples. Moreover, a total of 1552, 2295, 1271, 1195, 1588, and 1409 unigenes were specifically expressed in leaf_1 (HSYY_1), leaf_2 (HSYY_2), leaf_3 (HSYY_3), rhizome_1 (HSYK_1), rhizome_2 (HSYK_2), and rhizome_3 (HSYK_3), respectively. Overall, more unigenes were assembled in leaf tissues than in the rhizome tissues. Principal component analysis (PCA) of the expression levels of total transcriptome (Fig. S5), revealed that the leaf (HSYY_1, HSYY_2, HSYY_3) and rhizome (HSYK_1, HSYK_2, HSYK_3) treatment groups could be obviously distinguished. In other words, the three biological replicates of each sampling group clustered well together. Hence, the sequencing data can be used for further analysis. Identification and functional enrichment analysis of DEGs To identify the differentially expressed genes (DEGs), we compared the expression unigenes in leaf and rhizome tissues based on the false discovery rate (FDR) through the DESeq2 software. The results showed that 13,259 significant DEGs (FDR< 0.05 and |log2FC|≥1) were identified in the leaf versus rhizome tissues, which included 5536 up-regulated unigenes and 7723 down-regulated unigenes (Fig. 3, Table S2). A total of 13259 DEGs were enriched to 179 GO terms based on thresholds of FDR < 0.05, |log2FC| ≥ 1, p -value < 0.05 and p - value_corrected < 0.05 (Table S2). Fig. 4 indicates the top 20 enriched GO enrichment terms of DEGs. Among them, the membrane part (2964/8239), intrinsic component of membrane (2887/8239), integral component of membrane (2882/8239), and oxidoreductase activity (778/8239) were the most enriched GO terms of the DEGs (Fig. 4). By contrast, oxidoreductase activity, acting on iron-sulfur proteins as donors (13/8239), ferredoxin-NADP+reductase activity (10/8239) and omega peptidase activity (10/8239) were the least-enriched GO terms ( p -value_corrected < 0.05). To thoroughly understand the functions and interactions of DEGs between leaf and rhizome tissues, all DEGs were subjected to KEGG pathway enrichment analysis by mapping against the KEGG database using KEGG orthology based annotation system (KOBAS) software. The results indicated that 13,259 DEGs were involved in 144 metabolic pathways. There 23 pathways were enriched based on the p -value _uncorrected< 0.05 (Table S3). We found that these pathways were mainly related to the biosynthesis of other secondary metabolites, genetic information processing, and signal transduction, such as tropane, piperidine and pyridine alkaloid biosynthesis (map00960;28/4059), DNA replication (map03030; 12/4059), and phosphatidylinositol signaling system (map04070; 22/4059). The top 20 enriched pathways with a p -value of < 0.05 are presented in Fig. 5. We found that phenylpropanoid biosynthesis (map00940; 76/4059; 1.87%), photosynthesis (map04626; 54/4059; 1.33%), glutathione metabolism (map00480; 51/4059; 1.26%), glycine, serine and threonine metabolism (map00260; 49/4059; 1.21%), and flavonoid biosynthesis (map00941; 40/4059; 0.99%) were the highest enriched pathways. Detection and Characterization of SSR marker Using MISA 1.0 software, a total of 20,607 potential simple sequence repeats (SSRs) were identified across 13,770 unigenes (27.12%), of which 3464 unigenes contained more than one SSR. Among the 20,607 SSRs, the mono- and di-nucleotide repeats were the most abundant types (10,031, 48.68%; 5312, 25.78%, respectively), followed by tri- (4603, 22.34%), tetra- (440, 2.14%), hexa - (128, 0.62%), penta - (93, 0.45%), and nucleotide repeats. In addition, the repetitions of different SSRs types were calculated for four levels of repeat numbers 1–5, 6–10, 11–15, and >15 and are presented in Fig. S6. The results indicated that most mono-nucleotide repetitions were distributed in repeat number 6–10 (4385), followed by mono-nucleotide in 6–10 (4385), di-nucleotide in 6–10 (4095), tri-nucleotide 6–10 (2344), tri-nucleotide 1–5 (1991), di-nucleotide in 11–15 (874), and tetra-nucleotide 1–5 (290). Analysis of the SSR motif length distribution showed that most of the SSR motifs are less than 100 bp in length. The results indicated that 17,412，713，and 342 SSRs had motif length in the ranges of 10-50 , 51-100 and 101-150 bp, respectively. Identification of unigenes involved in terpene biosynthesis Comparison to the unigenes with the KEGG pathway database, indicated, a total of 174 unigenes related to terpene biosynthesis (Table S1), including terpenoid backbone biosynthesis (61), monoterpenoid biosynthesis (37), diterpenoid biosynthesis (19), sesquiterpenoid and triterpenoid biosynthesis (8), biosynthesis of terpenoids and steroids (3), and ubiquinone and other terpenoid-quinone biosynthesis (43). Many unigenes were involved in the pathway of terpenoid backbone biosynthesis, which belongs to the universal metabolic pathway map00900. Among these unigenes, DN2456_c0_g1 (DXS, K01662) is the first enzyme in the MEP pathway. DXS catalyze D-glyceraldehyde 3-phosphate (GA-3P) and pyruvate to form 1-deoxy-D-xylulose-5-phosphate (DOXP). The DXS gene is significantly less expressed in the rhizome than in the leaf. DN26885_c0_g1 (HMGR, K00021) is the main enzyme in the MVA pathway. It catalyzes 3-hydroxy-3-ethylglutaryl CoA (HMG-CoA) to form MVA. The expression of HMGR gene shows no difference in the leaf and rhizome tissues of D. opposita . qRT-PCR validation s of DEGs related to terpene biosynthesis To further analyze the consistency of RNA-seq data in this study, the DN39197_c0_g1 and DN37057_c0_g12 genes involved in terpene biosynthesis were selected for qRT-PCR validation. The results showed that a similar expression pattern in qRT-PCR analysis was observed in RNA-seq data from the leaf and rhizome tissues of D. opposita (Fig. 6). Further analysis showed that DXR (DN4169_c0_g1, K00099) and GGPS (DN17971_c0_g1, K13789) were the rate-limiting enzymes in the MVP and MEP pathways, which are the basic processes of a terpene synthesis pathway. Discussion In this study, six transcriptome sequencing analyses of RNA from leaf (HSYY_1, HSYY_2, HSYY_3) and rhizome (HSYK_1, HSYK_2, HSYK_3) tissues of D. opposita were conducted using the Illumina platform. More than five million raw reads and four million clean reads were obtained from life and rhizome library, respectively, which indicated that there more reads were assembled in leaf than in rhizome. Based on the large number of clean reads, 50,765 unigenes were assembled after evaluating the assembly results, of which, 46.98% unigenes were annotated by BLASTX analysis using the publicly available protein database. Otherwise, there are 49.5% unigenes were annotated from the D. alata [ 22 ]. This phenomenon shows that genomic data of Dioscorea species is relatively scarce in publicly database. There only D. rotundata [ 23 ], and D. dumetorum [ 24 ] genomic sequences were reported. The N50 is one of the most important indices used to evaluate the assembly quality. In the study, the N50 length of unigenes was 1825 bp, indicating that the sequence assembly was high in quality and suitable for intensive research. Moreover, the N50 could be compared with those of other species [ 46 – 48 ]. The average length of the unigenes was larger (1050 bp) (Table 2 ), than those in previous studies, such as those of Myrica rubra (531 bp) [ 49 ] and Sorbus pohuashanensis (770 bp) [ 50 ]. In this study, 24,840 and 25,925 unigenes were annotated and unannotated, respectively, in the NR, KEGG, COG, GO and Pfam database. However, the proportion of annotated long sequences was significantly higher than that of annotated short sequences. Similar results have been reported for other species [ 17 , 51 ]. The identified unigenes were compared among different Dioscorea species. A total of 5575 unigenes were matched with 315 Dioscorea species. Unigenes with the highese similarity were found in D. polystachya (16.41%), followed by D. bulbifera (5.70%), D. alatar (4.77%), and D. rotundata (3.28%), this result may indicate the evolutionary relationship among these species. The BLAST results (e-value 0.00001) of these Dioscorea species confirmed the correctness of transcriptome assembly and analysis in the present study. However, 296 Dioscorea species were matched to less than 100 unigenes of D. opposita. Although a small number of genes from Dioscorea species are available in several public database, complete transcriptome or genomic data are scarce. Indeed, to date, only the genomic resources for D. alata [ 22 ], D. rotundata [ 23 ], and D. dumetorum [ 24 ] have been reported. Therefore, enhancing research efforts on the genome of Dioscorea species is important and the main focus of the presenr study. Most D. alata tubers have white flesh, and purple-flesh tubers with high anthocyanidin content are produced by spontaneous variation. In the present study, elite purple-flesh cultivar (DP) and conventional white-flesh cultivar (DW) of D. alata were performed using Illumina sequencing [ 22 ]. The transcriptomes were compared with the two species of D. alata [ 22 ], and we compared with two tissues of D. opposita . A total of 49.5% (60,02/125,123) [ 22 ] and 48.93% (24,840/ 50,765) unigenes were annotated in D. alata and D. opposita , respectively, which indicates that few genome sequences of Dioscorea species are available in public database; thus, the ratio of annotation is fairly low. For example, the alignment lengths of gene NC_039707.1 (AYI69196.1) in D. alata and gene DN1708_c1_g2 in D. opposita are identical at 2250 bp (750 aa), and the similarity of these genes is 100%. Thus, these genes are believed to be the psbA gene, which encodes Photosystem I P700 chlorophyll a apoprotein A1 protein. 97 unigenes that were compared with D. dumetorum. Among these unigenes, the unigene DN2638_c1_g5 in D. opposita was identified to be gene NC_039691.1 (YP_009528139.1) in D. dumetorum ; these genes are 465 bp (155 aa), in length and encode ribosomal protein S7 in the chloroplast genome. A previous report matched a total of 11,115, 3542, 3304, 2943, and 2859 unigenes of D. alata with sequences from Vitis vinifera , Oryza sativa , Populus trichocarpa , Zea mays , and Ricinus communis , respectively, with high similarity [ 22 ]. However, in the present study, only 510, 170, 38, 280, and 80 unigenes were matched with sequences from V. vinifera , O. sativa , P. trichocarpa , Z. mays , and R. communis , respectively. Through the Blast2GO database, a large number of unigenes GO categories were assigned (Table 3 ), suggesting that the assembled unigenes represented a wide diversity of transcripts in the D. opposita genome, similar to in many previous reports [ 52 , 53 ]. In the three GO categories, reproduction and metabolic processes were the most abundant classes in BP cell and cell part activity were the most abundant classes in CC, which was consistent with some previous research [ 17 , 51 , 52 ]. Likewise, compared with the transcriptome of differentially expressed genes of leaf versus rhizome tissues by the DESeq2 software, there were 13,259 unigenes existing in significant DEGs, with 5536 up-regulated and 7723 down-regulated unigenes. These results showed that more unigenes abundance down-regulated in rhizome than in leaf tissue, which may be related to secondary metabolites of rhizome tissue. Moreover, most of these unigenes were related to transferase activity, membrane part, photosynthesis, plant–pathogen interaction, and plant hormone signal transduction. Trinity can be used with different k-mer size and the standed flag for de novo transcriptome assembly [ 54 ]. The developers describe comprehensive workflow for assembly validation on the trinity software. In this study, analyses by using five k-mer sizes of 25, 27, 30, 31 and 32 were performed with Trinity (Version v2.8.5) website. The de novo transcriptome assembly was generated using the stranded flag and 2 × 150nt paired-end reads. The assembly results were evaluated on the basis of assembly size, number of contigs, and BUSCO (Table S4). The number of complete BUSCOs is highest for a k-mer size of 30. Meanwhile, the N50 of 1825 bp is in the expected range for a transcriptome assembly. So, 30 were selected as the best k-mer size for the assembly and analysis of the transcriptome data of leaf and enzymes tissues from D. opposita. Based on previous reports [ 55 , 56 ], we speculate the expression levels of dioscin in different tissues of D. opposita may be different and are related to key enzymes of the biosynthetic pathway. Therefore, to explore the dioscin accumulation and translocation in leaf and rhizome tissues, we compared the differentially expressed genes (DEGs) and elementary identified terpene biosynthetic genes through the RNA-seq, which would lay a foundation to study the terpenoids biosynthesis pathway. We analyzed SSR sequences with length of 1–30 bp, and detected 16,417 unigenes. Among these unigenes, 2338 and 5782 were protein coding and non-coding sequences, respectively. The results also indicated that SSRs with lengths of 3 or multiples of 3 were identified more often within coding sequences than SSRs of other lengths, likely because a complete codon would not lead to a frameshift, but to the addition/deletion of an amino acid. Terpenoid biosynthesis requires common precursor substances, including isopentenyl diphosphate (IPP) and its allylic isomer dimethylallyl diphosphate (DMAPP), which are generated through the MVA and MEP pathways [ 55 , 56 ]. In the MEP pathway, DXR is a rata-limiting enzyme, that catalyze 1-deoxy-D-xylulose-5-phosphate (DXP) to generate MEP [ 57 ]. GGPS is the main enzyme in the terpene biosynthesis pathway, and can catalyze IPP to generate geranylgeranyl diphosphate (GGPP) [ 57 ]. In this study, many DEGs involved in terpenoids synthesi, such as diterpenoid biosynthesis, monoterpenoid biosynthesis and tropane, piperidine and pyridine alkaloid biosynthesis were detected. The expression patterns of the DXR (DN4169_c0_g1) and GGPS (DN17971_c0_g1) genes were studied through qRT-PCR, the results of which were consistent with the RNA-Seq tata. Our results also indicated that the abundance of DXR and GGPS genes is higher in the leaf tissue than in the rhizome tissue of D. opposita . This phenomenon may be related to the synthesis and accumulation of secondary metabolites in these tissues. Other enzymes, such as HMGR and hydroxymethylglutaryl-CoA synthase (HMGS) may also participate in the biosynthesis of terpenoids. These genes will be studied in great depth in future work. Conclusions In this study, the RNA-seq of leaf and rhizome tissues of D. opposita were obtained using an Illumina Hiseq 4000 platform. The results enriched the genomic data of D. opposita, and explore the different expression genes between leaf and rhizome tissues of D. opposita . Abbreviations MEP: 2-C-Methyl-D-Erythritol-4-Phosphate; GGPS: Geranylgeranyl pyrophosphate synthase; HMGR: Hydroxymethylglutaryl-CoA reductase; DXR:1-deoxy-D-xylulose 5-phosphate reductoisomerase; DXS: 1-deoxy-D-xylulose-5-phosphate synthase; GO:gene ontology; KEGG:Kyoto Encyclopedia of Genes and Genomes; PCR:polymerase chain reaction; MVA:Mevalonate acid; TPM:Transcripts per million reads; NR: Non-redundant; KOG:clusters of orthologous groups for eukaryotic complete genomes; RSEM: RNA-seq by Expectation Maximizattion; DEGs: differentially expressed genes; FDR: False discovery rate; SSRs: Simple sequence repeats. Declarations Ethics approval and consent to participate Not applicable Consent for publication Not applicable Competing interests Not applicable Acknowledgments We thank Dr. Shuang Zhou for her help with the data analysis. We would also thank manager Jinjing Chen（Shanghai Majorbio Bio-pham Technology Co., LTD, Shanghai,China） for her guidance and assistance in data analysis. Authors’ contributions All authors have participated extensively in the study and had proofread the final manuscript. DYH, JX and XGH conceived and designed the research. ZHL, JYX and RYC conducted the experiments. ZQM and HWX analyzed the data. DYH, RYC and ZHL wrote the manuscript. HWX and XLZ approved and reviewed the final manuscript. All authors read and approved the final manuscript. Funding This research was supported by the National Natural Science Foundation of China (U1404829)；Key project at central government level: The ability establishment of sustainable use for valuable Chinese medicine resources (2060302)；Science and technology project of Henan province(202102110156)；This research was also supported by a special fund for the construction of technical systems of the traditional Chinese medicine industry in Henan Province. Availability of data and materials The sequences can be found in the National Center for Biotechnology Information database (NCBI) SRA with the accession number PRJNA628588. Other datasets used and/or analyzed during the current study would be available from the corresponding author on reasonable request. Author details 1 Agricultural College, Henan University of Science and Technology, Luoyang 471023, China. 2 The Luoyang Engineering Research Center of Breeding and Utilization of Dao-di Herbs, Luoyang 471023, China. 3 Institute of Chinese Materia Medical, China Academy of Chinese Medical Sciences, Beijing, China. References Xie YF, Chang PR, Wang SJ, Yu JG, Ma XF. Preparation and properties of halloysite nanotubes/plasticized Dioscorea opposita Thunb. starch composites. Carbohydr Polym. 2011;83(1):186–91. Zhang LM, Liu P, Wang YG, Gao WY. Study on physico-chemical properties of dialdehyde yam starch with different aldehyde group contents. Thermochim Acta. 2011;512(1–2):196–201. Zhao GH, Li ZX, Chen ZD. Structural analysis and antitumor activity of RDPS-I polysaccharide from Chinese yam. Acta Pharm Sin. 2003;38(1):37–41. Wang SJ, Yu JL, Gao WY, Liu HY, Xiao PG. New starches from traditional Chinese medicine (TCM)-Chinese yam (Dioscorea opposita Thunb.) cultivars. Carbohydr Res. 2006;341(2):289–93. Chang SJ, Lee YC, Liu SY, Chang TW. Chinese yam (Dioscorea alata cv.Tainung No. 2) feeding exhibited antioxidative effects in hyperhomocysteinemia rats. J Agric Food Chem. 2004;52(6):1720–5. Ma FY, Zhang Y, Wen YR, Yao YN, Zhu JH, Liu XH, et al. Emulsification properties of polysaccharides from Dioscorea opposita Thunb. Food Chem. 2017;221(15):919–25. Chen YF, Zhu Q, Wu SJ. Preparation of oligosaccharides from Chinese yam and their antioxidant activity. Food Chem. 2015;173:1107–10. Herlina H. Deproteinase effect of hydrocolloid flour made of ‘‘Gembili Tuber” (Dioscorea esculenta L.) on chemical and technical functional properties. International Journal on Advanced Science Engineering Information Technology. 2015;5(4):298–302. Chen JT, Huang WJ, Mu Q. Active components and mineral elements analysis of three kinds of Chinese yam. Sci Tech food Industry. 2017;38(7):352–5. Guitton Y, Nicolè F, Moja S, Valot N, Legrand S, Jullien F, et al. Differential accumulation of volatile terpene and terpene synthase mRNAs during lavender (Lavandula angustifolia and L. x intermedia) inflorescence development. Physiol plant. 2010;138(2):150–63. Sun C, Li Y, Wu Q, Luo HM, Sun YZ, Song JY, et al. De novo sequencing and analysis of American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genom. 2010;11:262. Zhang HL, Yang B, Liu J, Guo DL, Hou J. Analysis of structural genes and key transcription factors related to anthocyanin biosynthesis in potato tubers. Sci Hortic. 2017;225(1):310–6. Li J, Liang Q, Li CF, Liu MD, Zhang YS. Comparative Transcriptome Analysis Identifies Putative Genes Involved in Dioscin Biosynthesis in Dioscorea zingiberensis . Molecules. 2018;23(2):454. Piskol R, Ramaswami G, Li JB. Reliable Identification of Genomic Variants from RNA-Seq Data. Am J Hum Genet. 2013;93(4):641–51. Li CF, Zhu YJ, Guo X, Sun C, Luo HM, Song JY, et al. Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng. BMC Genom. 2013;14(1):245. Xiao M, Zhang Y, Chen X, Lee EJ, Barber CJ, Chakrabarty R, et al. Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest. J Biotechnol. 2013;166(3):122–34. Wang DJ, Yang CL, Dong L, Zhu JC, Wang JP, Zhang SF. Comparative transcriptome analyses of drought-resistant and susceptible Brassica napus L. and development of EST-SSR markers by RNA-SEq. J Plant Biol. 2015;58:259–69. Pandey A, Swarnkar V, Pandey T, Srivastava P, Kanojiya S, Mishra DK, et al. Transcriptome and metabolite analysis reveal candidate genes of the cardiac Glycoside biosynthetic pathway from Calotropis procera. Sci Rep. 2016;6:34464. Akhtar MQ, Qamar N, Yadav P, Kulkarni P, Kumar A, Shasany AK. Comparative glandular trichome transcriptome based gene characterization reveals reasons for differential (–)-menthol biosynthesis in Mentha species. Physiol Plant. 2017;160(2):128–41. Wang XH, Hui F, Yang YC, Yang SH. Deep sequencing and transcriptome analysis to identify genes related to biosynthesis of aristolochic acid in Asarum heterotropoides . Sci Rep. 2018;8:17850. Li WG, Xu RR, Yan XG, Liang DM, Zhang L, Qin XY, et al. De novo leaf and root transcriptome analysis to explore biosynthetic pathway of Celangulin V in Celastrus angulatus maxim. BMC Genom. 2019;20:7. Wu ZG, Jiang W, Mantri N, Bao XQ, Chen SL, Tao ZM. Transciptome analysis reveals flavonoid biosynthesis regulation and simple sequence repeats in yam (Dioscorea alata L.) tubers. BMC Genom. 2015;16(1):1–12. Tamiru M, Natsume S, Takagi H, White B, Yaegashi H, Shimizu MY, et al. Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination. BMC Biol. 2017;15(1):86. Siadjeu C, Pucker B, Viehöver P, Albach DC, Weisshaar B. High Contiguity De Novo Genome Sequence Assembly of Trifoliate Yam (Dioscorea dumetorum) Using Long Read Sequencing. Genes. 2020;11(3):274. Chen YZ, Chen JZ. Comparative study on the pharmacological active components of four kinds of Chinese yam. Lishizhen Medicine Materia Medica Research. 2014;25(10):2389-91. Xu HW, Mo YW, Wang W, Wang H. OsPIN1a gene participates in regulating negative phototropism of rice roots. Rice Sci. 2014;21(2):83–9. Jia XP, Dai LF, Quan JZ, Li ZY, Fan BY. Clong and mutation sites analysis of a putative HD3A_like gene in eleven foxtail millet cultivars. Pak J Agr sci. 2018;55(1):29–32. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference Genome. Nat Biotechnol. 2011;29:644–52. Waterhouse RM, Seppey M, Simão FA, Manni M, Ioannidis P, Klioutchnikov G, et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol Biol Evol. 2018;35(3):543–8. Smith-Unna R, Boursnell C, Patro R, Hibberd JM, Kelly S. TransRate:reference-free quality assessment of de novo transcriptome assemblies. Genome Res. 2016;26(8):1134–44. El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47(D1):D427-32. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(suppl_1):D61-5. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(suppl_1):D480-4. Du F, Nu Er zi ya Ya Li mai mai ti, Hu QX, Ye ZYJ, Zhang D HJ. A Comparative Transcriptome Analysis Reveals Physiological Maturation Properties of Mycelia in Pleurotus tuoliensis . Genes. 2019; 10(9):703. Li B, Dewey CN. RSEM:accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. Damasco G, Shivakumar VS, Misciewicz TM, Daly DC, Fine PVA. Leaf Transcriptome Assembly of Protium copal (Burseraceae) and Annotation of Terpene Biosynthetic Genes. Genes. 2019;10(5):392. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84. Zhu YY, Li YQ, Xin DD, Chen WR, Shao X, Wang Y, et al. RNA-Seq-based transcriptome analysis of dormant flower buds of Chinese cherry (Prunus pseudocerasus). Gene. 2015;555(2):362–76. Niu J, Hou XY, Fang CL, An JY, Ha DL, Qiu L, et al. Transcriptome analysis of distinct Lindera glauca tissues revealed the differences in the unigenes related to terpenoid biosynthesis. Gene. 2015;559(1):22–30. Tang HB, Wang XY, Bowers JE, Ming R, Alam M, Paterson AH. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008;18(12):1944–54. Xie C, Mao XZ, Huang JJ, Ding Y, Wu JM, Dong S,et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(suppl_2):W316-22. Thiel T, Michalek W, Varshney R, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106:411–22. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCt method. Methods. 2001;25(4):402–8. Hu LS, Hao CY, Fan R, Wu BD, Tan LH, Wu HS. De Novo Assembly and Characterization of Fruit Transcriptome in Black Pepper (Piper nigrum). PLoS One. 2015;10(6):e0129822. Niño MC, Kang KK, Cho YG. Genome wide transcriptional response of papain like cysteine protease mediated resistance against xanthomonas oryzae pv. oryzae in rice. Plant Cell Rep. 2020;39:457–72. Liu SH, Yue TT, Ahmad MJ, Hu XW, Zhang XX, Deng TX, et al. Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle. Genes. 2020;11(1):68. Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, et al. Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-SEq. BMC Genom. 2012;13:19. Liu CC, Dou Y, Guan XL, Fu Q, Zhang Z, Hu ZH, et al. De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl. PLoS One. 2017;12(6):e0179219. Hou DY, Shi LC, Yang MM, Li J, Zhou S, Zhang HX, et al. De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using illumina platform. PLoS One. 2018;13(2):e0192610. Tang LD, Liu JM, Liu LH, Yu YH, Zhao HY, Lu W. De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead). Genes. 2020;11(2):144. Kim SH, Kim SW, Lim GH, Lyu JI, Choi HI, Jo YD, et al. Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation. PLoS One. 2020;15(1):e0228078. Haak M, Vinke S, Keller W, Droste JL, Rückert C, Kalinowski J, et al. High quality de Novo transcriptome assembly of Croton tiglium. Front Mol Biosci. 2018;5:62. Xu C, Wei H, Movahedi A, Sun WB, Ma XX, Li DW, et al. Evaluation, characterization, expression profiling, and functional analysis of DXS and DXR genes of Populus trichocarpa. Plant Physiol Biochem. 2019;142:94–105. Tong YR, Su P, Zhao YJ, Zhang M, Wang XJ, Liu YJ, et al. Molecular Cloning and Characterization of DXS and DXR Genes in the Terpenoid Biosynthetic Pathway of Tripterygium wilfordii. Int J Mol Sci. 2015;16(10):25516–35. Ma J, Xu ZS, Tan GF, Wang F, Xiong AS. Distinct transcription profile of genes involved in carotenoid biosynthesis among six different color carrot (Daucus carota L.) cultivars. Acta Biochim Biophys Sin. 2017;49(9):817–26. Additional Files Additional file 1: Table S1. The 145 pathways statistics Additional file 2: Table S2. Summary of significant DEGs with leaf versus rhizometissues Additional file 3: Table S3. Summary of KEGG pathway enrichment of significant DEGs with leaf versus rhizome tissues Additional file4: Table S4. Results of different k-mer sizes were assessed for the de novo transcriptome of leaf and rhizome tissues. Additional file 5: Fig. S1 Gene ontology (GO) level 2 term category classification of unigenes. The upper x-axis indicates the percentage of unigenes, and the lower x-axis indicates the number of unigenes. The y-axis indicates the specific categories of genes in the main categories. Additional file 6: Fig. S2 Gene ontology (GO) terms level 3 categories classification of unigenes Additional file 7: Fig. S3 The correlation coefficients heat map of the comparative transcriptome analysis of leaf (HSYY) and rhizome (HSYK) tissues. HSYY_1, HSYY_2, HSYY_3 and HSYK_1, HSYK_2, HSYK_3 indicate three biological replications of leaf and rhizome tissues, respectively, Additional file 8: Fig. S4 The venn diagram of six samples, including three biological replicates of leaf (HSYY_1, HSYY_2, and HSYY_3) and rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues, Additional file 9: Fig. S5 The PCA analysis of leaf (HSYY_1, HSYY_2, and HSYY_3) and rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues (A score plot of the first two PCs for samples grouped, the ellipses represent group membership assuming the 95% confidence). The x-axis and y-axis indicate the contribution degree of a principal component N (PCN) and principal component M (PCM) to the differentiated samples in the two-dimensional figure, respectively. Additional file 10: Fig. S6 The distribution statistics of SSR types. The x-axis indicates the different SSR types. The y-axis indicates the number of SSRs of different types corresponding to different repetition times, Supplementary Files Fig.S6.pdf Fig.S5.pdf Fig.S4.pdf Fig.S3.pdf Fig.S2.pdf TableS1S4.xlsx Fig.S1.pdf Fig.S5.pdf Fig.S6.pdf Fig.S4.pdf Fig.S3.pdf Fig.S2.pdf TableS1S4.xlsx Fig.S1.pdf Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-28829","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research","associatedPublications":[],"authors":[{"id":583416,"identity":"bab1676f-d55c-4d8f-b3b0-84241461806b","order_by":1,"name":"Dianyun Hou","email":"","orcid":"","institution":"Henan University of science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Dianyun","middleName":"","lastName":"Hou","suffix":""},{"id":583417,"identity":"508299d7-4d41-4fec-a726-9adc0d007c02","order_by":2,"name":"Ruiyang Cheng","email":"","orcid":"","institution":"China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Ruiyang","middleName":"","lastName":"Cheng","suffix":""},{"id":583418,"identity":"b910c0c7-2546-4d61-9e18-be4f769fa6a3","order_by":3,"name":"Zhihong Li","email":"","orcid":"","institution":"Henan university of science and technology","correspondingAuthor":false,"prefix":"","firstName":"Zhihong","middleName":"","lastName":"Li","suffix":""},{"id":583419,"identity":"c9fa4e55-955a-4d23-ba68-a3eb554e8c74","order_by":4,"name":"Jingya Xu","email":"","orcid":"","institution":"Henan university of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Jingya","middleName":"","lastName":"Xu","suffix":""},{"id":583420,"identity":"6bdca1a8-7151-4a12-b8d2-0d6cbf29bc95","order_by":5,"name":"Zhanqiang Ma","email":"","orcid":"","institution":"Henan university of science and technology","correspondingAuthor":false,"prefix":"","firstName":"Zhanqiang","middleName":"","lastName":"Ma","suffix":""},{"id":583421,"identity":"d586a094-99f1-43c8-bb4b-5755976e176f","order_by":6,"name":"Xingli Zhao","email":"","orcid":"","institution":"Henan university of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Xingli","middleName":"","lastName":"Zhao","suffix":""},{"id":583422,"identity":"4e544855-a564-4a7b-9686-3d73df548443","order_by":7,"name":"Huawei Xu","email":"","orcid":"","institution":"Henan University of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Huawei","middleName":"","lastName":"Xu","suffix":""},{"id":583423,"identity":"4e7f8204-7684-4224-924e-6f1a1149a401","order_by":8,"name":"Jiang Xu","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAo0lEQVRIiWNgGAWjYBACxgYGxgcJpGphNiBNCxCwSZCmnnlGjlnFg5pt0Qzshx8w/NxBjMN6zpjdSDh2O7eBJ82AsfcMMVrae8xuJDYAtTDkMDAzthGjpZnHrACshf8NsVqAtjCAtUgQbUvPsWIJkF/aJJ4ZHOwlRovhjOSNH3/U3M7t509++OAnUVoaoAw2ID5AhAYGBnmiVI2CUTAKRsHIBgC+njUefddBewAAAABJRU5ErkJggg==","orcid":"","institution":"China Academy of Chinese Medical Sciences","correspondingAuthor":true,"prefix":"","firstName":"Jiang","middleName":"","lastName":"Xu","suffix":""},{"id":583424,"identity":"f959f58e-f5e5-4ac4-80fd-1ad5dd2a6ef9","order_by":9,"name":"Xiaogai Hou","email":"","orcid":"","institution":"Henan University of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Xiaogai","middleName":"","lastName":"Hou","suffix":""}],"badges":[],"createdAt":"2020-05-14 09:07:53","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-28829/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-28829/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":1131527,"identity":"cc880886-f450-42f7-af34-c8d6dc7a3390","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":131358,"visible":true,"origin":"","legend":"Clusters of COG functional classification of D. opposita unigenes (E-value≤1.0 E-5) . The capital letters on the x-axis indicate the COG categories as listed to the right of the histogram. The y-axis indicates the number of unigenes.","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/1.png"},{"id":1131529,"identity":"79e05ea8-e005-463c-a446-65ff6ff7cd9b","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":152762,"visible":true,"origin":"","legend":"Classification map of KEGG metabolic pathways. The different colors indicate the different metabolic pathways.","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/2.jpg"},{"id":1131531,"identity":"76cf8305-5390-4712-bbe9-08a744691c12","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":56729,"visible":true,"origin":"","legend":"Volcano plots showing the expression level of each unigene of leaf (HSYY) and rhizome (HSYK) tissues (p-value≤ 0.05, |log2FC|≥1).","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/3.jpg"},{"id":1131534,"identity":"dfc149ff-092d-4ab6-93a3-7cabc60c815a","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":138088,"visible":true,"origin":"","legend":"The top 20 enriched GO enrichment terms of differentially expressed genes (DEGs) of leaf (HSYY) and rhizome (HSYK) tissues. The upper x-axis indicates the number of unigenes, and the lower x-axis indicates the false discovery rate (FDR). The larger the FDR value, the more significant the enrichment. The y-axis indicates the different metabolic pathways. Each point on the fold represents the number of each metabolic pathway.","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/4.jpg"},{"id":1131535,"identity":"cd339bf2-235c-4dce-b823-8cddd8b5973a","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":40061,"visible":true,"origin":"","legend":"The top 20 enriched KEGG pathways of DEGs between leaf (HSYY) and rhizome (HSYK) tissues. ","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/5.png"},{"id":1131537,"identity":"5162201f-bf4c-42c3-aeaa-811297ab76e2","added_by":"auto","created_at":"2020-05-19 17:33:44","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":39411,"visible":true,"origin":"","legend":"Expression of two selected unigenes as determined through qRT-PCR in comparison to the RNA-seq data.","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/6.png"},{"id":15667081,"identity":"e999d614-343c-4314-9316-27b5d5bb6c30","added_by":"auto","created_at":"2021-11-18 13:40:13","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":884634,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/ec018940-2930-4509-8ef8-a8798bfeadd3.pdf"},{"id":1131538,"identity":"93eddece-da18-4e5e-abb0-a2f04e5df2be","added_by":"auto","created_at":"2020-05-19 17:33:44","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":118952,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S6.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S6.pdf"},{"id":1131536,"identity":"52df34c6-22e7-4d76-b9ea-fc171a228f61","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":23294,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S5.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S5.pdf"},{"id":1131533,"identity":"812c1eb7-4390-4497-9177-03cabf42c801","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":29155,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S4.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S4.pdf"},{"id":1131532,"identity":"97d6233a-96a3-405f-9636-049d6b120396","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":18144,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S3.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S3.pdf"},{"id":1131530,"identity":"e7384f54-1760-40cf-a313-0e35110aa9d9","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":16026,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S2.pdf"},{"id":1131528,"identity":"e250c101-66fb-4a0e-9b87-2d39ea74c52f","added_by":"auto","created_at":"2020-05-19 17:33:43","extension":"xlsx","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":3871420,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1S4.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/TableS1S4.xlsx"},{"id":1131526,"identity":"e8ab20de-fd5c-4e07-afff-1487f3faa5b9","added_by":"auto","created_at":"2020-05-19 17:33:42","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":109902,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S1.pdf"},{"id":1131522,"identity":"338182e4-dcd1-4569-af43-b4228add09a7","added_by":"auto","created_at":"2020-05-19 17:33:40","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":23294,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S5.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S5.pdf"},{"id":1131524,"identity":"c7a42532-4570-4747-925f-53dee34d423e","added_by":"auto","created_at":"2020-05-19 17:33:40","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":118952,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S6.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S6.pdf"},{"id":1131519,"identity":"58bb832f-4066-4195-b0b8-e944764d8954","added_by":"auto","created_at":"2020-05-19 17:33:39","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":29155,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S4.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S4.pdf"},{"id":1131518,"identity":"51296b09-c118-442e-a9bd-1cd74db05333","added_by":"auto","created_at":"2020-05-19 17:33:39","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":18144,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S3.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S3.pdf"},{"id":1131516,"identity":"fac4a9a5-2c76-434b-be03-7ec03e3e000b","added_by":"auto","created_at":"2020-05-19 17:33:39","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":16026,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S2.pdf"},{"id":1131514,"identity":"2c62985a-61dd-46d7-b6ac-dda4d82e1ee7","added_by":"auto","created_at":"2020-05-19 17:33:39","extension":"xlsx","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":3871420,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1S4.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/TableS1S4.xlsx"},{"id":1131512,"identity":"5e00685d-25e0-460a-99e1-56dbed61cbd2","added_by":"auto","created_at":"2020-05-19 17:33:38","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":109902,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-28829/v1/Fig.S1.pdf"}],"financialInterests":"","formattedTitle":"\u003cp\u003eTranscriptome Assembly and Gene Expression Analysis of Leaf and Rhizome tissues of \u003cem\u003eDioscorea opposita\u003c/em\u003e\u003c/p\u003e","fulltext":[{"header":"Background","content":" \u003cp\u003e \u003cem\u003eDioscorea opposita\u003c/em\u003e Thunb (Family Dioscoreaceae), is an important staple in tropical regions and edible crop mainly distributed in parts of Africa, South America, and Asia [\u003cspan additionalcitationids=\"CR2\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e–\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. The rhizome of \u003cem\u003eD. opposita\u003c/em\u003e, known as “dioscoreae rhiaoma” or “Chinese yam”, has been listed as both a nutritious food source and an invigorant in traditional Chinese medicine for thousands of years [\u003cspan additionalcitationids=\"CR5\" citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e–\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Various types of natural constituents, including dopamine, starch, dioscin, polysaccharides, mucin, allantoin, choline, phytosterols, and essential amino acids, have been reported to be in \u003cem\u003eD. opposita\u003c/em\u003e [\u003cspan additionalcitationids=\"CR8\" citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e–\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThere were reports that 2-C-Methyl-D-Erythritol-4-Phosphate (MEP) and mevalonate acid (MVA) pathway are the early biosynthetic pathways of dioscin in many species [\u003cspan additionalcitationids=\"CR11 CR12\" citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e–\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e].It is well known that 1-deoxy-D-xylulose-5-phosphate synthase (DXS), 1-deoxy-D-xylulose 5-phosphate reductoisomerase (DXR), geranylgeranyl pyrophosphate synthase (GGPS) and hydroxymethylglutaryl-CoA reductase (HMGR) are the key enzyme genes, related to the biosynthesis of dioscin[\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. However, there are no reports on key enzyme genes, molecular biology studies or comprehensive genomic and transcriptomic data of \u003cem\u003eD. opposita\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eRNA-Seq is an efficient method for transcriptome analysis, which can also provide a highly sensitive and accurate platform at a low cost and is efficient for gene discovery and biosynthesis pathway analysis in non-model plants without knowing their genome sequences [\u003cspan additionalcitationids=\"CR15 CR16\" citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e–\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. RNA-Seq is currently a very popular method that is applied in the exploration of plants at transcriptomic levels. The ultimate aim is often to reveal genes that are involved in the biosynthetic pathways in secondary metabolism of many species [\u003cspan additionalcitationids=\"CR19 CR20\" citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e–\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Until now, the genomic resources of \u003cem\u003eD\u003c/em\u003e. alata [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], \u003cem\u003eD. rotundata\u003c/em\u003e [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], and \u003cem\u003eD. dumetorum\u003c/em\u003e [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e] have been reported. Nevertheless, as a non-model species, there is lack of enomic resources, in particular, information on active substance metabolism of \u003cem\u003eD. opposita\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eIn the study, the sample of leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e were collected from Wenxian County, called “Tie Gui Shan Yao”. It is a local cultivar of \u003cem\u003eD. opposita\u003c/em\u003e, belonging to one of the “four huaiqing Chinese medicines”. There have reported that the content of terpenoids in“Tie Gun Shan Yao” was higher than other \u003cem\u003eD. opposita\u003c/em\u003e [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn this study, transcriptome analysis was performed from leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e. The results will enrich the known genetic information for further gene function research and analyses, such as differential transcript abundance analysis. We also try to understand the transcript abundance difference and explore the molecular mechanisms of biosynthetic pathways in secondary metabolites in \u003cem\u003eD. opposita\u003c/em\u003e.\u003c/p\u003e "},{"header":"Materials And Methods","content":"\u003cp\u003e\u003cstrong\u003ePlant materials and RNA extraction\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLeaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e were collected from Wenxian County, Jiaozuo City, Henan Province of China, on October 2017（Wenxian County Weather Station： N34°56′49.54″, E113°10′18.54″; T:15℃）. Three biological replicates from three leaf tissues (HSYY_1, HSYY_2, and HSYY_3) and three rhizome tissues (HSYK_1, HSYK_2, and HSYK_3) were used for the analysis. All samples were flash frozen in liquid nitrogen after cleaning and stored at –80 °C until RNA extraction.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Total RNAs were extracted from the leaf and rhizome tissues using the TRIzol® Reagent (Invitrogen, USA) following the product description with slight improvement [26, 27]. The extracted RNAs were checked through NanoDrop™2000 spectrophotometer (NanoDrop Technologies, USA) and RNA integrity was assessed using electrophoresis in 1.0% agarose gel. Furthermore, the RNA concentration was determined by the RNA Agilent Bioanalyzer 2100 system (Agilent, CA, USA) to guarantee the RNA Integrity Number (RIN) values were not less than 7.0.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTranscriptome sequencing and de novo assembly \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe cDNA libraries were constructed using the TruSeq-TM RNA sample preparation Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. The cDNA libraries were enriched through a polymerase chain reaction (PCR) with 15 cycles. Then, the target band was extracted and recycled by a 2% low range ultra-agarose gel, and quantified by TBS380 Mini-Fluorometer. The suitable cDNA libraries were sent to Shanghai Majorbio Biopharm Biotechnology Co., Ltd. (Shanghai, China) where the RNA-Seq was performed using the Illumina HiSeq4000 SBS Kit (300 cycles, Illumina, San Diego, CA, USA) on an Illumina HiSeq4000 platform.\u003c/p\u003e\n\u003cp\u003eThe raw reads were cleaned and assessed following the steps. Firstly, the adaptor sequence of reads was removed, and discarded the reads that have no inserted fragments due to the reason of adaptor self-connection; Secondly, low quality bases (quality value less than 20) were trimed off at the end of the sequence (3 '). If there are still bases with quality value less than 10 in the remaining reads, the whole reads are eliminated, otherwise it is retained.Thirdly, reads with “N” bases was also removed. Lastly, the reads less than 100bp in length were discarded after removing adapter and quality pruning. All the above mentioned processes were performed carried out using the SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (Version 1.33) software [28]. The high quality clean reads of leaves and rhizomes of \u003cem\u003eD. opposita\u003c/em\u003e were de novo assembled by Trinity (Version v2.8.5) software with operation parameters of k-mer size of 25, 27, 30, 31and 32(max normalize read 50，Min\u0026nbsp;contig\u0026nbsp;length 200). The assembly results were assessed based on the assembly size, number of contigs, N50 and complete BUSCOs. Then, 2x150bp reads were generated and filtered, and more than 200 bp contigs were produced. At last, assembly completeness and quality also were assessed with Benchmarking Universal Single-Copy Orthologs (BUSCO, version 3.0.2) and TransRate package (version 1.0.3) [29, 30].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunctional annotation and classification of unigenes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe unigenes were compared with the NCBI BlastX program (E-value \u0026lt;1 E-5) against the non-redundant (NR), KEGG (Kyoto Encyclopedia of Genes and Genomes), KOG (clusters of orthologous groups for eukaryotic complete genomes), GO (gene ontology) and Pfam [31-33].\u003c/p\u003e\n\u003cp\u003eThe functional classification of these unigenes was performed using Blast2GO and KEGG. The GO analysis was categorized into biological processes, molecular functions, and cellular component [34]. The pathway distributions were determined based on the KEGG database [35]. Unigenes were also aligned to the Evolutionary Genealogy of Genes: Non-Supervised Orthologous Groups (eggNOG) database for further functional prediction and classification [36].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDifferential transcript abundance analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe clean reads were used to determine alignment counts and to quantify transcript abundance by the RSEM (RNA-seq by Expectation Maximizattion, http://deweylab.github.io/RSEM/) package [37]. The unigene and transcript abundance levels were analyzed by TPM (transcripts per million reads) [38].\u003c/p\u003e\n\u003cp\u003eThe software DESeq2 (Version 1.24.0, http://bioconductor.org/packages/ stats/bioc/DESeq2/) was used to analyze the differentially expressed genes (DEGs) between leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e. The false discovery rate (FDR) can be used to assess \u003cem\u003ep\u003c/em\u003e-value thresholds in multiple testing [39, 40]. The significance of DEGs was defined with a threshold of FDR \u0026lt; 0.05 and an absolute value of log2-fold change of ≥2.0 (|log2fc|≥1.0) [41]. The DEGs enrichment analyses were carried out for GO and KEGG terms by the GOatools and KOBAS software [42, 43].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSSR Mining and\u003c/strong\u003e\u003cstrong\u003e analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe unigenes were searched for the prediction and detection of simple sequence repeats (SSR) through the MISA 1.0 software with default parameters [44]. The motifs of mono-, di-, tri-, tetra-, penta-, and hexa- nucleotide repeats were identified. A minimum of five repetitions were considered as the search criteria in MISA. The parameters of definement of microsatellites (unit size/minimum number of repeats) were used to predict SSRs. That is Mono- (1/10), Di- (2/6), Tri-(3/5), Tetra- (4/5), Penta- (5/5) and Hexa-nucleotide (6/5). And the maximal number of bases interrupting 2 SSRs in a compound microsatellite was 100.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eExpression analysis of terpene biosynthetic related DEGs\u003c/strong\u003e\u003c/p\u003e \u003cp\u003eWe selected 2 DEGs, DN39197_c0_g1 and DN37057_c0_g1, which were the key enzyme genes in the terpene synthesis from \u003cem\u003eD. opposita\u003c/em\u003e, to validate their expression using quantitative real-time PCR (qRT-PCR). Total RNAs from leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e were extracted through the same procedures as used for RNA-seq in the study. The RNA from each sample was reverse transcribed using a reverse transcription kit (PC17-TRUE script, Aidlab Biotechnologies Co., Ltd, China). The \u003cem\u003eDoActin\u003c/em\u003e gene was used as an internal control. Primers for the 2 DEGs and \u003cem\u003eDoActin\u003c/em\u003e were designed in the Primer Premier 5.0 software. All the primers sequences are listed in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. The reaction system (10\u0026nbsp;µL) included 5\u0026nbsp;µL of SYBR Premix Ex TaqTM, 4\u0026nbsp;µL of cDNA template, and 0.5\u0026nbsp;µL of each primer (10\u0026nbsp;µM). The qRT-PCR reactions were conducted on a Roche Light Cycle 96 real-time PCR system (Roche, USA) with SYBR green fluorescent dye (TaKaRa, China). PCR amplification was initiated at 95\u0026nbsp;°C for 30\u0026nbsp;s, followed by 40 cycles of 95\u0026nbsp;°C for 5\u0026nbsp;s, 60\u0026nbsp;°C for 30\u0026nbsp;s, and 72\u0026nbsp;°C for 30\u0026nbsp;s. Three replicates were set for the experiment. Relative transcript abundance levels were calculated by the 2-ΔΔCt method [\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003e\u003cb\u003ePrimers for qRT-PCR of 2 DEGs and\u003c/b\u003e \u003cspan type=\"BoldItalic\" class=\"BoldItalic\" name=\"Emphasis\"\u003eDoActin\u003c/span\u003e\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003c/colgroup\u003e\u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eUnigene ID\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnnotation\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePrimer sequence(5’ − 3’)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eDN39197_c0_g1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eDXR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eDXR-F: GGTGCCGATTACGACCATATC\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eDXR-R: GGATTGGCAGACGCATATCA\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eDN37057_c0_g1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eGGPS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eGGPS-F: ATCGGTGGTGGTTCTGATTC\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eGGPS-R: CTGCAGTCTTCCCAAGTTCT\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003eActin\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eActin-F: CTAAGGCCAACAGAGAGAAGATG\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eActin-R: GGTACGGCCACTTGCATAAA\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e "},{"header":"Results","content":" \u003cp\u003e \u003cstrong\u003eTranscriptome sequencing and de novo assembly\u003c/strong\u003e\u003c/p\u003e \u003cp\u003eWe obtained averages of 48,786,977 and 55,202,249 raw reads for the rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues and leaf (HSYY_1, HSYY_2, and HSYY_3) tissues, respectively, which consisted of 7,366,833,477\u0026nbsp;bp and 8,335,539,549\u0026nbsp;bp. We obtained 46,982,031 and 53,689,620 clean reads, including 7,053,884,682\u0026nbsp;bp and 8,067,994,222\u0026nbsp;bp clean bases, respectively. The sequences can be found in NCBI SRA with the accession number PRJNA628588. Moreover, the average GC content of rhizome and leaf tissues was 46.8% and 46.7%, respectively. The k-mer size of 25, 27, 30, 31 and 32 were selected to evaluate the assembly results based on the assembly size, number of contigs, N50 and complete BUSCOs. The results showed that the k-mer size of 30 for the final assembly results in the best final assembly. Trinity assembly resulted in 50,765 unigenes with fragment size ranging from 201\u0026nbsp;bp to 16,991\u0026nbsp;bp, and the average length of each unigene was 1049.6\u0026nbsp;bp. In addition, all of the unigenes consisted of a total of 53,283,125\u0026nbsp;bp. The unigenes were retained with an N50 length of 1825\u0026nbsp;bp and N90 length of 2797\u0026nbsp;bp. The TransRate assembly and Benchmarking UniversalSingle-Copy Orthologs (BUSCO) completeness assessment score were 0.26405 and 68.9%, respectively. The detailed sequencing and assembly results for the leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e are shown in Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eOptimized assembly results of unigenes from \u003cem\u003eD. opposita\u003c/em\u003e rhizome and leaf tissues\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"2\"\u003e \u003c/colgroup\u003e\u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eType\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnigene\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal sequence num\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e50,765\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal sequence base (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e53,283,125\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLargest (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e16,991\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSmallest (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e201\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAverage length(bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1,050\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eN50 (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1,825\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eN90 (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,797\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTransRate score\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e0.26405\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eComplete single copy BUSCOs\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e68.9%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003e \u003cb\u003eFunctional annotation of\u003c/b\u003e \u003cspan type=\"BoldItalic\" class=\"BoldItalic\" name=\"Emphasis\"\u003eD. opposita\u003c/span\u003e \u003cb\u003eunigenes\u003c/b\u003e\u003c/p\u003e \u003cp\u003eThe functional annotation results indicated that a total of 24,840 unigenes were annotated (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Among them, different unigenes were matched in different databases, for example, 23,848 (46.98%) unigenes were matched in the NR database.In addition, 18,531 (36.50%) unigenes were matched in the Swiss-Prot database.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe annotation statistics of \u003cem\u003eD. opposita\u003c/em\u003e\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"2\"\u003e \u003c/colgroup\u003e\u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eType\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eUnigene number (percent%)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e23,848(46.98)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSwiss-Prot\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e18,531(36.50)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePfam\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e18,805(37.04)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eCOG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e21,973(43.28)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGO\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e20,014(39.42)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKEGG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e10,617(20.91)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal_annot\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e24,840(48.93)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTotal\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e50,765(1)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003eThe length distributions of annotated and unannotated unigenes results indicated that most annotated and unannotated unigenes were distrbuted in the 200–499\u0026nbsp;bp range (annotated unigenes: 5970/24,840; unannotated unigenes: 16,040/25,925), followed by the 500–799\u0026nbsp;bp ranges (annotated unigenes: 3446/24,840; unannotated unigenes:5566/25,925), and 800–1099\u0026nbsp;bp range (annotated unigenes: 2394/24,840; unannotated unigenes: 2194/25,925). Long sequences of unigens were clearly more likely to be annotated than those with short sequences. A total of 7151 unigenes measuring over 2000\u0026nbsp;bp in length were annotated, and only 371 unigenes were unannotated.\u003c/p\u003e \u003cp\u003e\u003cstrong\u003eFunctional classification results of\u003c/strong\u003e \u003cstrong\u003eunigenes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo further evaluate the annotated unigene data, the unigene sequences were searched against the cluster of orthologous groups of proteins (COG) databases. In total, 23,139 unigenes assigned to 23 COG categories under the following types: information storage and processing, cellular processes and signaling, metabolism, and poorly characterized (Fig. 1). Among them, alphabetized“Posttranslational\u0026nbsp;modification, protein turnover, chaperones”(O) (1373; 5.93 %) was the categories with the most members in the COG databases.\u003c/p\u003e\n\u003cp\u003eUsing the Blast2GO database, the 20,014 unigenes in the gene ontology (GO) analysis could be summarized into three main GO categories biological process (BP), molecular function (MF), and cellular component (CC) on two levels. In the BP category, the cellular process (7345;36.70%), metabolic process (6480; 32.389%), and biological regulation (2125; 10.62%) were represented the most. Likewise, cell part (7756; 38.75%), membrane part (7018; 35.07%),organelle (4292; 21.44%) and organelle part (3154; 15.76%)were the most abundant in the CC category. Furthermore, the MF category mainly consisted of binding (10,949; 54.71%), catalytic activity (10,794; 53.93%), transporter activity (1074; 5.37%), and structural molecule activity (730; 3.65%) (Fig. S1).\u003c/p\u003e\n\u003cp\u003eIn order to perform an in-depth analysis of unigene function based on the GO database,all of the three categories (BP, MF, and CC) were subdivided into level 3 and level 4 terms in order to determine their function. For example, cellular process is the main level 2 term from the BP category. In the level 3 terms, cellular process was assigned to the following specific functions: cellular metabolic process (GO: 0044237; 5615 unigenes; 28.06%), cellular component organization (GO:0016043; 1147 unigenes; 5.73%), cellular response to stimulus (GO:0051716; 555 unigenes; 2.78%), and signal transduction (GO:0007165; 526 unigenes; 2.63%), etc. (Fig. S2).\u003c/p\u003e\n\u003cp\u003eThrough the KEGG database BLAST, 10,617 unigene sequences were assigned to 145 pathways (Table S1). Among them, the three pathways with the greatest number of unigene sequences were ribosome (map03010; 605 unigenes; 5.70%), plant-pathogen interaction (map04626; 281 unigenes; 2.65%),protein processing in endoplasmic reticulum (map04141; 276 unigenes; 2.60%), and thermogenesis (map04714; 248 unigenes; 2.34%) (Table S1). In addition, all the pathways were divided into five main types of KEGG metabolic pathways (Fig. 2).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eExpression analysis of unigenes from \u003c/strong\u003e\u003cstrong\u003eleaf\u003c/strong\u003e\u003cstrong\u003e and rhizome tissues\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to explore the differential expression in leaf and rhizome tissues of \u003cem\u003eD. opposit\u003c/em\u003e\u003cem\u003ea\u003c/em\u003e, the clean reads from the leaf and rhizome library were compared with the unigene database by using the RNA-seq by RSEM package. The results indicated that 41,692,444 reads, respectively 77.62% of the original clean reads (53,689,620), were aligned in leaf tissues (HSYY: average reads of HSYY_1, HSYY_2, and HSYY_3). Likewise, 37,978,278 reads, representing 80.83% of the original clean reads (46,982.031), were aligned in rhizome tissues (HSYK: average reads of HSYK_1, HSYK_2, and HSYK_3).\u003c/p\u003e\n\u003cp\u003eIn this study, the sample correlation analysis was performed using RSEM software (Tables S1). A high correlation coefficient was obtained among the three biological replicates of each tissue. The correlation coefficients of rhizome_1 (HSYK_1) with rhizome_2 (HSYK_2) and rhizome_3 (HSYK_3) were 0.935 and 0.913, respectively. The correlation coefficients of leaf_1 (HSYY_1) with leaf_2 (HSYY_2) and leaf_3 (HSYY_3) were 0.919 and 0.935, respectively. The correlation coefficient between two different of leaf and rhizome tissues was significantly lower than that among the three biological replicates of each leaf and rhizome tissue. The results indicated that all three biological replicates of each tissue are all clustered together, thus suggesting that the sequencing data are suitable for further biological analysis (Fig. S3).\u003c/p\u003e\n\u003cp\u003eVenn analysis of samples with a unigene expression greater than 2 was performed using RSEM software (Fig. S4). A total of 14,627 unigenes were shared among the six samples. Moreover, a total of 1552, 2295, 1271, 1195, 1588, and 1409 unigenes were specifically expressed in leaf_1 (HSYY_1), leaf_2 (HSYY_2), leaf_3 (HSYY_3), rhizome_1 (HSYK_1), rhizome_2 (HSYK_2), and rhizome_3 (HSYK_3), respectively. Overall, more unigenes were assembled in leaf tissues than in the rhizome tissues.\u003c/p\u003e\n\u003cp\u003ePrincipal component analysis (PCA) of the expression levels of total transcriptome (Fig. S5), revealed that the leaf (HSYY_1, HSYY_2, HSYY_3) and rhizome (HSYK_1, HSYK_2, HSYK_3) treatment groups could be obviously distinguished. In other words, the three biological replicates of each sampling group clustered well together. Hence, the sequencing data can be used for further analysis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification and functional enrichment analysis of DEGs \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo identify the differentially expressed genes (DEGs), we compared the expression unigenes in leaf and rhizome tissues based on the false discovery rate (FDR) through the DESeq2 software. The results showed that 13,259 significant DEGs (FDR\u0026lt; 0.05 and |log2FC|≥1) were identified in the leaf versus rhizome tissues, which included 5536 up-regulated unigenes and 7723 down-regulated unigenes (Fig. 3, Table S2).\u003c/p\u003e\n\u003cp\u003eA total of 13259 DEGs were enriched to 179 GO terms based on thresholds of FDR \u0026lt; 0.05, |log2FC| ≥ 1, \u003cem\u003ep\u003c/em\u003e-value \u0026lt; 0.05 and \u003cem\u003ep\u003c/em\u003e\u003cem\u003e-\u003c/em\u003evalue_corrected \u0026lt; 0.05 (Table S2). Fig. 4 indicates the top 20 enriched GO enrichment terms of DEGs. Among them, the membrane part (2964/8239), intrinsic component of membrane (2887/8239), integral component of membrane (2882/8239), and oxidoreductase activity (778/8239) were the most enriched GO terms of the DEGs (Fig. 4). By contrast, oxidoreductase activity, acting on iron-sulfur proteins as donors (13/8239), ferredoxin-NADP+reductase activity (10/8239) and omega peptidase activity (10/8239) were the least-enriched GO terms (\u003cem\u003ep\u003c/em\u003e-value_corrected \u0026lt; 0.05).\u003c/p\u003e\n\u003cp\u003eTo thoroughly understand the functions and interactions of DEGs between leaf and rhizome tissues, all DEGs were subjected to KEGG pathway enrichment analysis by mapping against the KEGG database using KEGG orthology based annotation system (KOBAS) software. The results indicated that 13,259 DEGs were involved in 144 metabolic pathways. There 23 pathways were enriched based on the \u003cem\u003ep\u003c/em\u003e-value _uncorrected\u0026lt; 0.05 (Table S3). \u0026nbsp;We found that these pathways were mainly related to the biosynthesis of other secondary metabolites, genetic information processing, and signal transduction, such as tropane, piperidine and pyridine alkaloid biosynthesis (map00960;28/4059), DNA replication (map03030; 12/4059), and phosphatidylinositol signaling system (map04070; 22/4059).\u003c/p\u003e\n\u003cp\u003eThe top 20 enriched pathways with a \u003cem\u003ep\u003c/em\u003e-value of \u0026lt; 0.05 are presented in Fig. 5. We found that phenylpropanoid biosynthesis (map00940; 76/4059; 1.87%), photosynthesis (map04626; 54/4059; 1.33%), glutathione metabolism (map00480; 51/4059; 1.26%), glycine, serine and threonine metabolism (map00260; 49/4059; 1.21%), and flavonoid biosynthesis (map00941; 40/4059; 0.99%) were the highest enriched pathways.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDetection and Characterization of SSR marker\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eUsing MISA 1.0 software, a total of 20,607 potential simple sequence repeats (SSRs) were identified across 13,770 unigenes (27.12%), of which 3464 unigenes contained more than one SSR. Among the 20,607 SSRs, the mono- and di-nucleotide repeats were the most abundant types (10,031, 48.68%; 5312, 25.78%, respectively), followed by tri- (4603, 22.34%), tetra- (440, 2.14%), \u0026nbsp;hexa - (128, 0.62%), penta - (93, 0.45%), and nucleotide repeats. In addition, the repetitions of different SSRs types were calculated for four levels of repeat numbers 1–5, 6–10, 11–15, and \u0026gt;15 and are presented in Fig. S6. The results indicated that most mono-nucleotide repetitions were distributed in repeat number 6–10 (4385), followed by mono-nucleotide in 6–10 (4385), di-nucleotide in 6–10 (4095), tri-nucleotide 6–10 (2344), tri-nucleotide 1–5 (1991), di-nucleotide in 11–15 (874), and tetra-nucleotide 1–5 (290). Analysis of the SSR motif length distribution showed that most of the SSR motifs are less than 100 bp in length. The results indicated that 17,412，713，and 342 SSRs had motif length in the ranges of 10-50 , 51-100 and 101-150 bp, respectively.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of unigenes involved in terpene biosynthesis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u0026nbsp; \u0026nbsp;\u0026nbsp;Comparison to the unigenes with the KEGG pathway database, indicated, a total of 174 unigenes related to terpene biosynthesis (Table S1), including terpenoid backbone biosynthesis (61), monoterpenoid biosynthesis (37), diterpenoid biosynthesis (19), sesquiterpenoid and triterpenoid biosynthesis (8), biosynthesis of terpenoids and steroids (3), and ubiquinone and other terpenoid-quinone biosynthesis (43). Many unigenes were involved in the pathway of terpenoid backbone biosynthesis, which belongs to the universal metabolic pathway map00900. Among these unigenes, DN2456_c0_g1 (DXS, K01662) is the first enzyme in the MEP pathway. DXS catalyze D-glyceraldehyde 3-phosphate (GA-3P) and pyruvate to form 1-deoxy-D-xylulose-5-phosphate (DOXP). The \u003cem\u003eDXS\u003c/em\u003e gene is significantly less expressed in the rhizome than in the leaf. DN26885_c0_g1 (HMGR, K00021) is the main enzyme in the MVA pathway. It catalyzes 3-hydroxy-3-ethylglutaryl CoA (HMG-CoA) to form MVA. The expression of \u003cem\u003eHMGR gene shows no difference \u003c/em\u003ein the leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eqRT-PCR validation\u003c/strong\u003e\u003cstrong\u003es\u003c/strong\u003e\u003cstrong\u003e of DEGs related to terpene biosynthesis \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo further analyze the consistency of RNA-seq data in this study, the DN39197_c0_g1 and DN37057_c0_g12 genes involved in terpene biosynthesis were selected for qRT-PCR validation. The results showed that a similar expression pattern in qRT-PCR analysis was observed in RNA-seq data from the leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e (Fig. 6). Further analysis showed that DXR (DN4169_c0_g1, K00099) and GGPS (DN17971_c0_g1, K13789) were the rate-limiting enzymes in the MVP and MEP pathways, which are the basic processes of a terpene synthesis pathway.\u003c/p\u003e"},{"header":"Discussion","content":" \u003cp\u003eIn this study, six transcriptome sequencing analyses of RNA from leaf (HSYY_1, HSYY_2, HSYY_3) and rhizome (HSYK_1, HSYK_2, HSYK_3) tissues of \u003cem\u003eD. opposita\u003c/em\u003e were conducted using the Illumina platform. More than five million raw reads and four million clean reads were obtained from life and rhizome library, respectively, which indicated that there more reads were assembled in leaf than in rhizome. Based on the large number of clean reads, 50,765 unigenes were assembled after evaluating the assembly results, of which, 46.98% unigenes were annotated by BLASTX analysis using the publicly available protein database. Otherwise, there are 49.5% unigenes were annotated from the \u003cem\u003eD. alata\u003c/em\u003e [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. This phenomenon shows that genomic data of \u003cem\u003eDioscorea\u003c/em\u003e species is relatively scarce in publicly database. There only \u003cem\u003eD. rotundata\u003c/em\u003e [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], and \u003cem\u003eD. dumetorum\u003c/em\u003e [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e] genomic sequences were reported. The N50 is one of the most important indices used to evaluate the assembly quality. In the study, the N50 length of unigenes was 1825\u0026nbsp;bp, indicating that the sequence assembly was high in quality and suitable for intensive research. Moreover, the N50 could be compared with those of other species [\u003cspan additionalcitationids=\"CR47\" citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e–\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e]. The average length of the unigenes was larger (1050\u0026nbsp;bp) (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), than those in previous studies, such as those of \u003cem\u003eMyrica rubra\u003c/em\u003e (531\u0026nbsp;bp) [\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e] and \u003cem\u003eSorbus pohuashanensis\u003c/em\u003e (770\u0026nbsp;bp) [\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn this study, 24,840 and 25,925 unigenes were annotated and unannotated, respectively, in the NR, KEGG, COG, GO and Pfam database. However, the proportion of annotated long sequences was significantly higher than that of annotated short sequences. Similar results have been reported for other species [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe identified unigenes were compared among different \u003cem\u003eDioscorea\u003c/em\u003e species. A total of 5575 unigenes were matched with 315 \u003cem\u003eDioscorea\u003c/em\u003e species. Unigenes with the highese similarity were found in \u003cem\u003eD. polystachya\u003c/em\u003e (16.41%), followed by \u003cem\u003eD. bulbifera\u003c/em\u003e (5.70%), \u003cem\u003eD. alatar\u003c/em\u003e (4.77%), and \u003cem\u003eD. rotundata\u003c/em\u003e (3.28%), this result may indicate the evolutionary relationship among these species. The BLAST results (e-value 0.00001) of these \u003cem\u003eDioscorea\u003c/em\u003e species confirmed the correctness of transcriptome assembly and analysis in the present study. However, 296 \u003cem\u003eDioscorea\u003c/em\u003e species were matched to less than 100 unigenes of \u003cem\u003eD. opposita.\u003c/em\u003e Although a small number of genes from \u003cem\u003eDioscorea\u003c/em\u003e species are available in several public database, complete transcriptome or genomic data are scarce. Indeed, to date, only the genomic resources for \u003cem\u003eD. alata\u003c/em\u003e [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], \u003cem\u003eD. rotundata\u003c/em\u003e [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], and \u003cem\u003eD. dumetorum\u003c/em\u003e [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e] have been reported. Therefore, enhancing research efforts on the genome of \u003cem\u003eDioscorea\u003c/em\u003e species is important and the main focus of the presenr study.\u003c/p\u003e \u003cp\u003eMost \u003cem\u003eD. alata\u003c/em\u003e tubers have white flesh, and purple-flesh tubers with high anthocyanidin content are produced by spontaneous variation. In the present study, elite purple-flesh cultivar (DP) and conventional white-flesh cultivar (DW) of \u003cem\u003eD. alata\u003c/em\u003e were performed using Illumina sequencing [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. The transcriptomes were compared with the two species of \u003cem\u003eD. alata\u003c/em\u003e [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], and we compared with two tissues of \u003cem\u003eD. opposita\u003c/em\u003e. A total of 49.5% (60,02/125,123) [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e] and 48.93% (24,840/ 50,765) unigenes were annotated in \u003cem\u003eD. alata\u003c/em\u003e and \u003cem\u003eD. opposita\u003c/em\u003e, respectively, which indicates that few genome sequences of \u003cem\u003eDioscorea\u003c/em\u003e species are available in public database; thus, the ratio of annotation is fairly low. For example, the alignment lengths of gene NC_039707.1 (AYI69196.1) in \u003cem\u003eD. alata\u003c/em\u003e and gene DN1708_c1_g2 in \u003cem\u003eD. opposita\u003c/em\u003e are identical at 2250\u0026nbsp;bp (750 aa), and the similarity of these genes is 100%. Thus, these genes are believed to be the \u003cem\u003epsbA\u003c/em\u003e gene, which encodes Photosystem I P700 chlorophyll a apoprotein A1 protein. 97 unigenes that were compared with \u003cem\u003eD. dumetorum.\u003c/em\u003e Among these unigenes, the unigene DN2638_c1_g5 in \u003cem\u003eD. opposita\u003c/em\u003e was identified to be gene NC_039691.1 (YP_009528139.1) in \u003cem\u003eD. dumetorum\u003c/em\u003e; these genes are 465\u0026nbsp;bp (155 aa), in length and encode ribosomal protein S7 in the chloroplast genome.\u003c/p\u003e \u003cp\u003eA previous report matched a total of 11,115, 3542, 3304, 2943, and 2859 unigenes of \u003cem\u003eD. alata\u003c/em\u003e with sequences from \u003cem\u003eVitis vinifera\u003c/em\u003e, \u003cem\u003eOryza sativa\u003c/em\u003e, \u003cem\u003ePopulus trichocarpa\u003c/em\u003e, \u003cem\u003eZea mays\u003c/em\u003e, and \u003cem\u003eRicinus communis\u003c/em\u003e, respectively, with high similarity [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. However, in the present study, only 510, 170, 38, 280, and 80 unigenes were matched with sequences from \u003cem\u003eV. vinifera\u003c/em\u003e, \u003cem\u003eO. sativa\u003c/em\u003e, \u003cem\u003eP. trichocarpa\u003c/em\u003e, \u003cem\u003eZ. mays\u003c/em\u003e, and \u003cem\u003eR. communis\u003c/em\u003e, respectively.\u003c/p\u003e \u003cp\u003eThrough the Blast2GO database, a large number of unigenes GO categories were assigned (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e), suggesting that the assembled unigenes represented a wide diversity of transcripts in the \u003cem\u003eD. opposita\u003c/em\u003e genome, similar to in many previous reports [\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. In the three GO categories, reproduction and metabolic processes were the most abundant classes in BP cell and cell part activity were the most abundant classes in CC, which was consistent with some previous research [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e, \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eLikewise, compared with the transcriptome of differentially expressed genes of leaf versus rhizome tissues by the DESeq2 software, there were 13,259 unigenes existing in significant DEGs, with 5536 up-regulated and 7723 down-regulated unigenes. These results showed that more unigenes abundance down-regulated in rhizome than in leaf tissue, which may be related to secondary metabolites of rhizome tissue. Moreover, most of these unigenes were related to transferase activity, membrane part, photosynthesis, plant–pathogen interaction, and plant hormone signal transduction.\u003c/p\u003e \u003cp\u003eTrinity can be used with different k-mer size and the standed flag for de novo transcriptome assembly [\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. The developers describe comprehensive workflow for assembly validation on the trinity software. In this study, analyses by using five k-mer sizes of 25, 27, 30, 31 and 32 were performed with Trinity (Version v2.8.5) website. The de novo transcriptome assembly was generated using the stranded flag and 2 × 150nt paired-end reads. The assembly results were evaluated on the basis of assembly size, number of contigs, and BUSCO (Table S4). The number of complete BUSCOs is highest for a k-mer size of 30. Meanwhile, the N50 of 1825\u0026nbsp;bp is in the expected range for a transcriptome assembly. So, 30 were selected as the best k-mer size for the assembly and analysis of the transcriptome data of leaf and enzymes tissues from \u003cem\u003eD. opposita.\u003c/em\u003e\u003c/p\u003e \u003cp\u003eBased on previous reports [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e, \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e], we speculate the expression levels of dioscin in different tissues of \u003cem\u003eD. opposita\u003c/em\u003e may be different and are related to key enzymes of the biosynthetic pathway. Therefore, to explore the dioscin accumulation and translocation in leaf and rhizome tissues, we compared the differentially expressed genes (DEGs) and elementary identified terpene biosynthetic genes through the RNA-seq, which would lay a foundation to study the terpenoids biosynthesis pathway.\u003c/p\u003e \u003cp\u003eWe analyzed SSR sequences with length of 1–30\u0026nbsp;bp, and detected 16,417 unigenes. Among these unigenes, 2338 and 5782 were protein coding and non-coding sequences, respectively. The results also indicated that SSRs with lengths of 3 or multiples of 3 were identified more often within coding sequences than SSRs of other lengths, likely because a complete codon would not lead to a frameshift, but to the addition/deletion of an amino acid.\u003c/p\u003e \u003cp\u003eTerpenoid biosynthesis requires common precursor substances, including isopentenyl diphosphate (IPP) and its allylic isomer dimethylallyl diphosphate (DMAPP), which are generated through the MVA and MEP pathways [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e, \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. In the MEP pathway, DXR is a rata-limiting enzyme, that catalyze 1-deoxy-D-xylulose-5-phosphate (DXP) to generate MEP [\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e]. GGPS is the main enzyme in the terpene biosynthesis pathway, and can catalyze IPP to generate geranylgeranyl diphosphate (GGPP) [\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e]. In this study, many DEGs involved in terpenoids synthesi, such as diterpenoid biosynthesis, monoterpenoid biosynthesis and tropane, piperidine and pyridine alkaloid biosynthesis were detected. The expression patterns of the DXR (DN4169_c0_g1) and GGPS (DN17971_c0_g1) genes were studied through qRT-PCR, the results of which were consistent with the RNA-Seq tata. Our results also indicated that the abundance of DXR and GGPS genes is higher in the leaf tissue than in the rhizome tissue of \u003cem\u003eD. opposita\u003c/em\u003e. This phenomenon may be related to the synthesis and accumulation of secondary metabolites in these tissues. Other enzymes, such as HMGR and hydroxymethylglutaryl-CoA synthase (HMGS) may also participate in the biosynthesis of terpenoids. These genes will be studied in great depth in future work.\u003c/p\u003e "},{"header":"Conclusions","content":" \u003cp\u003eIn this study, the RNA-seq of leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e were obtained using an Illumina Hiseq 4000 platform. The results enriched the genomic data of D. opposita, and explore the different expression genes between leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e.\u003c/p\u003e "},{"header":"Abbreviations","content":"\u003cp\u003eMEP: 2-C-Methyl-D-Erythritol-4-Phosphate; GGPS: Geranylgeranyl pyrophosphate synthase; HMGR: Hydroxymethylglutaryl-CoA reductase; DXR:1-deoxy-D-xylulose 5-phosphate reductoisomerase; DXS:\n1-deoxy-D-xylulose-5-phosphate synthase; GO:gene ontology; KEGG:Kyoto Encyclopedia of Genes and Genomes; PCR:polymerase chain reaction; MVA:Mevalonate acid; TPM:Transcripts per million reads; NR: Non-redundant; KOG:clusters of orthologous groups for eukaryotic complete genomes; RSEM: RNA-seq by Expectation Maximizattion; DEGs: differentially expressed genes; FDR: False discovery rate; SSRs: Simple sequence repeats.\n\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Not applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;Consent for publication \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Not applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe thank Dr. Shuang Zhou for her help with the data analysis. We would also thank manager Jinjing Chen（Shanghai Majorbio Bio-pham Technology Co., LTD, Shanghai,China） for her guidance and assistance in data analysis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors’ contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll authors have participated extensively in the study and had proofread the final manuscript. DYH, JX and XGH conceived and designed the research. ZHL, JYX and RYC conducted the experiments. ZQM and HWX analyzed the data. DYH, RYC and ZHL wrote the manuscript. HWX and XLZ approved and reviewed the final manuscript. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research was supported by the National Natural Science Foundation of China (U1404829)；Key project at central government level: The ability establishment of sustainable use for valuable Chinese medicine resources (2060302)；Science and technology project of Henan province(202102110156)；This research was also supported by a special fund for the construction of technical systems of the traditional Chinese medicine industry in Henan Province.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe sequences can be found in the National Center for Biotechnology Information database (NCBI) SRA with the accession number PRJNA628588. Other datasets used and/or analyzed during the current study would be available from the corresponding author on reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor details\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e1\u003c/sup\u003eAgricultural College, Henan University of Science and Technology, Luoyang 471023, China.\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e2\u003c/sup\u003eThe Luoyang Engineering Research Center of Breeding and Utilization of Dao-di Herbs, Luoyang 471023, China. \u003csup\u003e3\u003c/sup\u003e Institute of Chinese Materia Medical, China Academy of Chinese Medical Sciences, Beijing, China.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e \u003cspan\u003eXie YF, Chang PR, Wang SJ, Yu JG, Ma XF. Preparation and properties of halloysite nanotubes/plasticized \u003cem\u003eDioscorea opposita\u003c/em\u003e Thunb. starch composites. Carbohydr Polym. 2011;83(1):186–91.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eZhang LM, Liu P, Wang YG, Gao WY. Study on physico-chemical properties of dialdehyde yam starch with different aldehyde group contents. Thermochim Acta. 2011;512(1–2):196–201.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eZhao GH, Li ZX, Chen ZD. Structural analysis and antitumor activity of RDPS-I polysaccharide from Chinese yam. Acta Pharm Sin. 2003;38(1):37–41.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eWang SJ, Yu JL, Gao WY, Liu HY, Xiao PG. New starches from traditional Chinese medicine (TCM)-Chinese yam (Dioscorea opposita Thunb.) cultivars. Carbohydr Res. 2006;341(2):289–93.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eChang SJ, Lee YC, Liu SY, Chang TW. Chinese yam (Dioscorea alata cv.Tainung No. 2) feeding exhibited antioxidative effects in hyperhomocysteinemia rats. J Agric Food Chem. 2004;52(6):1720–5.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eMa FY, Zhang Y, Wen YR, Yao YN, Zhu JH, Liu XH, et al. Emulsification properties of polysaccharides from Dioscorea opposita Thunb. Food Chem. 2017;221(15):919–25.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eChen YF, Zhu Q, Wu SJ. Preparation of oligosaccharides from Chinese yam and their antioxidant activity. Food Chem. 2015;173:1107–10.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eHerlina H. Deproteinase effect of hydrocolloid flour made of ‘‘Gembili Tuber” (Dioscorea esculenta L.) on chemical and technical functional properties. International Journal on Advanced Science Engineering Information Technology. 2015;5(4):298–302.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eChen JT, Huang WJ, Mu Q. Active components and mineral elements analysis of three kinds of Chinese yam. Sci Tech food Industry. 2017;38(7):352–5.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eGuitton Y, Nicolè F, Moja S, Valot N, Legrand S, Jullien F, et al. Differential accumulation of volatile terpene and terpene synthase mRNAs during lavender (Lavandula angustifolia and L. x intermedia) inflorescence development. Physiol plant. 2010;138(2):150–63.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eSun C, Li Y, Wu Q, Luo HM, Sun YZ, Song JY, et al. De novo sequencing and analysis of American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genom. 2010;11:262.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eZhang HL, Yang B, Liu J, Guo DL, Hou J. Analysis of structural genes and key transcription factors related to anthocyanin biosynthesis in potato tubers. Sci Hortic. 2017;225(1):310–6.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLi J, Liang Q, Li CF, Liu MD, Zhang YS. Comparative Transcriptome Analysis Identifies Putative Genes Involved in Dioscin Biosynthesis in \u003cem\u003eDioscorea zingiberensis\u003c/em\u003e. Molecules. 2018;23(2):454.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003ePiskol R, Ramaswami G, Li JB. Reliable Identification of Genomic Variants from RNA-Seq Data. Am J Hum Genet. 2013;93(4):641–51.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLi CF, Zhu YJ, Guo X, Sun C, Luo HM, Song JY, et al. Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng. BMC Genom. 2013;14(1):245.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eXiao M, Zhang Y, Chen X, Lee EJ, Barber CJ, Chakrabarty R, et al. Transcriptome analysis based on next-generation sequencing of non-model plants producing specialized metabolites of biotechnological interest. J Biotechnol. 2013;166(3):122–34.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eWang DJ, Yang CL, Dong L, Zhu JC, Wang JP, Zhang SF. Comparative transcriptome analyses of drought-resistant and susceptible Brassica napus L. and development of EST-SSR markers by RNA-SEq. J Plant Biol. 2015;58:259–69.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003ePandey A, Swarnkar V, Pandey T, Srivastava P, Kanojiya S, Mishra DK, et al. Transcriptome and metabolite analysis reveal candidate genes of the cardiac Glycoside biosynthetic pathway from Calotropis procera. Sci Rep. 2016;6:34464.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eAkhtar MQ, Qamar N, Yadav P, Kulkarni P, Kumar A, Shasany AK. Comparative glandular trichome transcriptome based gene characterization reveals reasons for differential (–)-menthol biosynthesis in Mentha species. Physiol Plant. 2017;160(2):128–41.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eWang XH, Hui F, Yang YC, Yang SH. Deep sequencing and transcriptome analysis to identify genes related to biosynthesis of aristolochic acid in \u003cem\u003eAsarum heterotropoides\u003c/em\u003e. Sci Rep. 2018;8:17850.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLi WG, Xu RR, Yan XG, Liang DM, Zhang L, Qin XY, et al. De novo leaf and root transcriptome analysis to explore biosynthetic pathway of Celangulin V in \u003cem\u003eCelastrus angulatus\u003c/em\u003e maxim. BMC Genom. 2019;20:7.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eWu ZG, Jiang W, Mantri N, Bao XQ, Chen SL, Tao ZM. Transciptome analysis reveals flavonoid biosynthesis regulation and simple sequence repeats in yam (Dioscorea alata L.) tubers. BMC Genom. 2015;16(1):1–12.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eTamiru M, Natsume S, Takagi H, White B, Yaegashi H, Shimizu MY, et al. Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination. BMC Biol. 2017;15(1):86.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eSiadjeu C, Pucker B, Viehöver P, Albach DC, Weisshaar B. High Contiguity De Novo Genome Sequence Assembly of Trifoliate Yam (Dioscorea dumetorum) Using Long Read Sequencing. Genes. 2020;11(3):274.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eChen YZ, Chen JZ. Comparative study on the pharmacological active components of four kinds of Chinese yam. Lishizhen Medicine Materia Medica Research. 2014;25(10):2389-91.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eXu HW, Mo YW, Wang W, Wang H. OsPIN1a gene participates in regulating negative phototropism of rice roots. Rice Sci. 2014;21(2):83–9.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eJia XP, Dai LF, Quan JZ, Li ZY, Fan BY. Clong and mutation sites analysis of a putative HD3A_like gene in eleven foxtail millet cultivars. Pak J Agr sci. 2018;55(1):29–32.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eGrabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference Genome. Nat Biotechnol. 2011;29:644–52.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eWaterhouse RM, Seppey M, Simão FA, Manni M, Ioannidis P, Klioutchnikov G, et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol Biol Evol. 2018;35(3):543–8.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eSmith-Unna R, Boursnell C, Patro R, Hibberd JM, Kelly S. TransRate:reference-free quality assessment of de novo transcriptome assemblies. Genome Res. 2016;26(8):1134–44.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eEl-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47(D1):D427-32.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003ePruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(suppl_1):D61-5.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eCamacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eConesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eKanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(suppl_1):D480-4.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eDu F, Nu Er zi ya Ya Li mai mai ti, Hu QX, Ye ZYJ, Zhang D HJ. A Comparative Transcriptome Analysis Reveals Physiological Maturation Properties of Mycelia in \u003cem\u003ePleurotus tuoliensis\u003c/em\u003e. Genes. 2019; 10(9):703.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLi B, Dewey CN. RSEM:accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eDamasco G, Shivakumar VS, Misciewicz TM, Daly DC, Fine PVA. Leaf Transcriptome Assembly of Protium copal (Burseraceae) and Annotation of Terpene Biosynthetic Genes. Genes. 2019;10(5):392.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eBenjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eZhu YY, Li YQ, Xin DD, Chen WR, Shao X, Wang Y, et al. RNA-Seq-based transcriptome analysis of dormant flower buds of Chinese cherry (Prunus pseudocerasus). Gene. 2015;555(2):362–76.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eNiu J, Hou XY, Fang CL, An JY, Ha DL, Qiu L, et al. Transcriptome analysis of distinct Lindera glauca tissues revealed the differences in the unigenes related to terpenoid biosynthesis. Gene. 2015;559(1):22–30.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eTang HB, Wang XY, Bowers JE, Ming R, Alam M, Paterson AH. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008;18(12):1944–54.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eXie C, Mao XZ, Huang JJ, Ding Y, Wu JM, Dong S,et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(suppl_2):W316-22.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eThiel T, Michalek W, Varshney R, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106:411–22.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLivak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCt method. Methods. 2001;25(4):402–8.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eHu LS, Hao CY, Fan R, Wu BD, Tan LH, Wu HS. De Novo Assembly and Characterization of Fruit Transcriptome in Black Pepper (Piper nigrum). PLoS One. 2015;10(6):e0129822.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eNiño MC, Kang KK, Cho YG. Genome wide transcriptional response of papain like cysteine protease mediated resistance against xanthomonas oryzae pv. oryzae in rice. Plant Cell Rep. 2020;39:457–72.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLiu SH, Yue TT, Ahmad MJ, Hu XW, Zhang XX, Deng TX, et al. Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle. Genes. 2020;11(1):68.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eFeng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, et al. Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-SEq. BMC Genom. 2012;13:19.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eLiu CC, Dou Y, Guan XL, Fu Q, Zhang Z, Hu ZH, et al. De novo transcriptomic analysis and development of EST-SSRs for Sorbus pohuashanensis (Hance) Hedl. PLoS One. 2017;12(6):e0179219.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eHou DY, Shi LC, Yang MM, Li J, Zhou S, Zhang HX, et al. De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using illumina platform. PLoS One. 2018;13(2):e0192610.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eTang LD, Liu JM, Liu LH, Yu YH, Zhao HY, Lu W. De Novo Transcriptome Identifies Olfactory Genes in Diachasmimorpha longicaudata (Ashmead). Genes. 2020;11(2):144.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eKim SH, Kim SW, Lim GH, Lyu JI, Choi HI, Jo YD, et al. Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation. PLoS One. 2020;15(1):e0228078.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eHaak M, Vinke S, Keller W, Droste JL, Rückert C, Kalinowski J, et al. High quality de Novo transcriptome assembly of Croton tiglium. Front Mol Biosci. 2018;5:62.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eXu C, Wei H, Movahedi A, Sun WB, Ma XX, Li DW, et al. Evaluation, characterization, expression profiling, and functional analysis of DXS and DXR genes of Populus trichocarpa. Plant Physiol Biochem. 2019;142:94–105.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eTong YR, Su P, Zhao YJ, Zhang M, Wang XJ, Liu YJ, et al. Molecular Cloning and Characterization of DXS and DXR Genes in the Terpenoid Biosynthetic Pathway of Tripterygium wilfordii. Int J Mol Sci. 2015;16(10):25516–35.\u003c/span\u003e \u003c/li\u003e \u003cli\u003e \u003cspan\u003eMa J, Xu ZS, Tan GF, Wang F, Xiong AS. Distinct transcription profile of genes involved in carotenoid biosynthesis among six different color carrot (Daucus carota L.) cultivars. Acta Biochim Biophys Sin. 2017;49(9):817–26.\u003c/span\u003e \u003c/li\u003e\u003c/ol\u003e"},{"header":"Additional Files","content":"\u003cp\u003eAdditional file 1: Table S1. The 145 pathways statistics\u003c/p\u003e\n\u003cp\u003eAdditional file 2: Table S2. Summary of significant DEGs with leaf versus rhizometissues\u003c/p\u003e\n\u003cp\u003eAdditional file 3: Table S3. Summary of KEGG pathway enrichment of significant DEGs with leaf versus rhizome tissues\u003c/p\u003e\n\u003cp\u003eAdditional file4: Table S4. Results of different k-mer sizes were assessed for the de novo transcriptome of leaf and rhizome tissues.\u003c/p\u003e\n\u003cp\u003eAdditional file 5: Fig. S1 Gene ontology (GO) level 2 term category classification of unigenes. The upper x-axis indicates the percentage of unigenes, and the lower x-axis indicates the number of unigenes. The y-axis indicates the specific categories of genes in the main categories.\u003c/p\u003e\n\u003cp\u003eAdditional file 6: Fig. S2 Gene ontology (GO) terms level 3 categories classification of unigenes\u003c/p\u003e\n\u003cp\u003eAdditional file 7: Fig. S3 The correlation coefficients heat map of the comparative transcriptome analysis of leaf (HSYY) and rhizome (HSYK) tissues. HSYY_1, HSYY_2, HSYY_3 and HSYK_1, HSYK_2, HSYK_3 indicate three biological replications of leaf and rhizome tissues, respectively,\u003c/p\u003e\n\u003cp\u003eAdditional file 8: Fig. S4 The venn diagram of six samples, including three biological replicates of leaf (HSYY_1, HSYY_2, and HSYY_3) and rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues,\u003c/p\u003e\n\u003cp\u003eAdditional file 9: Fig. S5 The PCA analysis of leaf (HSYY_1, HSYY_2, and HSYY_3) and rhizome (HSYK_1, HSYK_2, and HSYK_3) tissues (A score plot of the first two PCs for samples grouped, the ellipses represent group membership assuming the 95% confidence). The x-axis and y-axis indicate the contribution degree of a principal component N (PCN) and principal component M (PCM) to the differentiated samples in the two-dimensional figure, respectively.\u003c/p\u003e\n\u003cp\u003eAdditional file 10: Fig. S6 The distribution statistics of SSR types. The x-axis indicates the different SSR types. The y-axis indicates the number of SSRs of different types corresponding to different repetition times,\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Dioscorea opposita, leaf and rhizome, transcriptome assembly, expression, terpene biosynthesis","lastPublishedDoi":"10.21203/rs.3.rs-28829/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-28829/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e\u003cem\u003e Dioscorea opposita\u003c/em\u003e is a monocotyledonous plant in the Dioscoreaceae family. The rhizome of \u003cem\u003eD. opposita\u003c/em\u003e, known as “dioscoreae rhiaoma” or “Chinese yam”, is not only a well-known traditional Chinese medicine, but also a popular health food.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eIn this study, the leaf and rhizome tissues of \u003cem\u003eD. opposita\u003c/em\u003e were collected. Their transcriptome sequenceand were obtained using an Illumina Hiseq 4000 platform. The data was analyzed by bioinformatics software. Likewise, the consistent expression pattern of DEGs was validated with qRT-PCR analysis.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003e\u0026nbsp;A total of 46,982,031 and 53,689,620 clean reads from leaf and rhizome tissues were acquired, respectively (SRP: PRJNA628588). 50,765 unigenes were obtained with an average length of 1050 base pairs (bp). A total of 13,259 significant differentially expressed genes (DEGs) were obtained from leaf versus rhizome tissues, including 5536 up-regulated and 7723 down-regulated unigenes. In addition, 20,607 simple sequence repeats (SSRs) were identified in this study. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eThis study not only enriched the transcriptomic data of \u003cem\u003eD. opposita\u003c/em\u003e, but also laid a foundation for the analysis of terpenoid synthesis pathway.\u003c/p\u003e","manuscriptTitle":"Transcriptome Assembly and Gene Expression Analysis of Leaf and Rhizome tissues of Dioscorea opposita","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2020-05-19 17:33:38","doi":"10.21203/rs.3.rs-28829/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"2b4de7b8-28a6-4b40-a53e-51058ca5396b","owner":[],"postedDate":"May 19th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":102957,"name":"Integrative \u0026 Complementary Medicine"}],"tags":[],"updatedAt":"2020-05-29T02:41:39+00:00","versionOfRecord":[],"versionCreatedAt":"2020-05-19 17:33:38","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-28829","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-28829","identity":"rs-28829","version":["v1"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}