Next Article in Journal
Alternative Splicing in Breast Cancer and the Potential Development of Therapeutic Tools
Previous Article in Journal
Exonization of an Intronic LINE-1 Element Causing Becker Muscular Dystrophy as a Novel Mutational Mechanism in Dystrophin Gene
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 8
Issue 10
10.3390/genes8100257
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Sequencing of Full-Length cDNA Transcripts of Stems and Leaves in Dendrobium officinale

by
Liu He
1,
Shuhua Fu
2,
Zhichao Xu
1,
Jun Yan
3,
Jiang Xu
4,
Hong Zhou
1,
Jianguo Zhou
1,
Xinlian Chen
1,
Ying Li
1,
Kin Fai Au
2,5,* and
Hui Yao
1,*
1
Key Laboratory of Bioactive Substances and Resources Utilization of Chinese Herbal Medicine, Ministry of Education, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100193, China
2
Department of Internal Medicine, University of Iowa, Iowa City, IA 52242, USA
3
Department of Plant Genetic Breeding, College of Agriculture and Biotechnology, China Agricultural University, Beijing 100193, China
4
Institute of Chinese Materia Medica, Academy of Chinese Medical Sciences, Beijing 100700, China
5
Department of Biostatistics, University of Iowa, Iowa City, IA 52242, USA
*
Authors to whom correspondence should be addressed.
Genes 2017, 8(10), 257; https://doi.org/10.3390/genes8100257
Submission received: 18 August 2017 / Revised: 16 September 2017 / Accepted: 2 October 2017 / Published: 5 October 2017
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Dendrobium officinale is an extremely valuable orchid used in traditional Chinese medicine, so sought after that it has a higher market value than gold. Although the expression profiles of some genes involved in the polysaccharide synthesis have previously been investigated, little research has been carried out on their alternatively spliced isoforms in D. officinale. In addition, information regarding the translocation of sugars from leaves to stems in D. officinale also remains limited. We analyzed the polysaccharide content of D. officinale leaves and stems, and completed in-depth transcriptome sequencing of these two diverse tissue types using second-generation sequencing (SGS) and single-molecule real-time (SMRT) sequencing technology. The results of this study yielded a digital inventory of gene and mRNA isoform expressions. A comparative analysis of both transcriptomes uncovered a total of 1414 differentially expressed genes, including 844 that were up-regulated and 570 that were down-regulated in stems. Of these genes, one sugars will eventually be exported transporter (SWEET) and one sucrose transporter (SUT) are expressed to a greater extent in D. officinale stems than in leaves. Two glycosyltransferase (GT) and four cellulose synthase (Ces) genes undergo a distinct degree of alternative splicing. In the stems, the content of polysaccharides is twice as much as that in the leaves. The differentially expressed GT and transcription factor (TF) genes will be the focus of further study. The genes DoSWEET4 and DoSUT1 are significantly expressed in the stem, and are likely to be involved in sugar loading in the phloem.

1. Introduction

The perennial epiphytic herb Dendrobium officinale Kimura et Migo is a rare lithophytic orchid that lives in dark and damp rocky environments in Chinese subtropical mountainous regions. Dendrobii Officinalis Caulis, the dried stem of D. officinale, has been used as a tonic in traditional Chinese medicine (TCM) for thousands of years to diminish inflammation, counteract senescence, and improve immunity [1]. Water-soluble polysaccharides are the dominant active ingredients in D. officinale [2,3,4,5], comprising glucose, mannose, galactose, xylose, and arabinose in different proportions [6,7]. As carbohydrates are produced in green plant leaves as a result of photosynthetic activity, assimilated carbon is moved immediately out of the leaves and is mainly converted to polysaccharides in the stem and roots [8].
The D. officinale genes associated with polysaccharide synthesis, including glucose transporters (GT) and glycoside hydrolases, have been evaluated in previous studies [9,10,11]. In addition to polysaccharide synthesis, sugar translocation from the leaves to the stems is particularly noteworthy, because the sugar content of plants generally depends on metabolism and transport [12].
Sucrose is actively transported from the leaves to other tissues to sustain growth and metabolism. SWEET (sugars will eventually be exported transporters) refers to a new class of transporters that are essential for the movement of sugars between cells in both plants and animals [13]. Investigations have shown that SWEETs are widespread in plants; 17 have been identified in Arabidopsis thaliana, 21 in Oryza sativa, 23 in Sorghum bicolor, and 29 in Solanum lycopersicum [13,14,15].
In A. thaliana, the SWEETs of sucrose, glucose, galactose, fructose, and 2-deoxyglucose have been shown to act as bidirectional transporters, mediating both sugar uptake and efflux across the plasma membrane in leaves, roots, flowers, pollens, and seeds [13,16,17]. In rice, for example, the transporters OsSWEET11 and OsSWEET14 participated in the efflux of photosynthesized sucrose from the leaves [13,18], with the former required for pollen development. In addition to SWEETs, sucrose transporters (SUTs) are also responsible for sucrose importation from the apoplast into phloem cells; ZmSUT1, for example, is known to be responsible for sucrose loading and unloading in the phloem [19,20]. Despite this background, however, information regarding sugar accumulation in the stem remains limited as the model plants investigated to date accumulate just small amounts of these carbohydrates in this region. Nevertheless, a larger polysaccharide volume has been identified in stems of D. officinale [1]; but the molecular mechanisms that underlie this accumulation necessitate further investigation.
Because of interest in the medicinal properties of D. officinale, extensive DNA and RNA sequencing studies have been carried out recently on this extremely valuable plant [9,10,21,22]. However, molecular information on the differing gene expression in relation to sugar translocation between the leaves and stems of D. officinale remains limited. The advent of PacBio sequencing has meant that full length gene transcripts are now available while the use of second generation sequencing (SGS) has led to improvements in sequence quality. Hybrid sequencing strategies are now generally both more affordable and of higher quality than the use of single sequencing alone; indeed, the combined use of single-molecule real-time (SMRT) long reads and SGS short reads has enabled the detection of sensitive isoforms and predictions based on genome sequences [23,24]. The production of gene isoforms has also led to an increase in the structural and functional variety of gene products.
In this study, we utilized the combination of SGS and SMRT sequencing technology to generate two organ-specific full-length transcriptomes for D. officinale. We applied this strategy to dissected leaf and stem samples, thereby enabling more correlation of specific expression information from resultant transcriptional data to stems, where polysaccharides are accumulated. The results of this study are significant for future research on the molecular mechanisms of sugar translocation and polysaccharide accumulation in D. officinale.

2. Material and Methods

2.1. Plant Materials and RNA Sample Preparation

All the D. officinale seedlings (three-year-old) were harvested from a plastic greenhouse (Lin’an City, Zhejiang Province, China). A total of forty independent seedlings were collected. The leaves and stem of each seedling were separated and wrapped individually in pieces of aluminum paper, frozen in liquid nitrogen, and stored at −80 °C prior to RNA isolation and polysaccharide extraction. Each leaf (stem) RNA sample came from five seedlings. Dendrobium fimbriatum Hook. (three-year-old) leaves and stems were cut from a local greenhouse (Nanning City, Guangxi Province, China). The collection method used is the same as above.

2.2. Polysaccharide Determination

We followed the method outlined in the Pharmacopeia of the People’s Republic of China (2015) to determine the polysaccharide content of Dendrobium nobile, performing three parallel experiments under the same conditions. The anhydrous d-glucose (batch number 110833–201205) used for all experiments was purchased from the Chinese Food and Drug Inspection Institute, and an initial reference solution was prepared at a concentration of 93 µg/mL.

2.3. RNA Preparation and cDNA Synthesis

We isolated total RNA using a RNeasy plus mini kit (#74134, Qiagen, Duesseldorf, Germany). RNA quality was then assessed using EtBr-stained 0.8% agarose gels, while concentration and integrity were evaluated using an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA). Transcriptome sequencing was carried out using Illumina Hiseq 1500 (Illumina, San Diego, CA, USA) and PacBio RS II systems (Pacific Biosciences, Menlo Park, CA, USA), following the manufacturer’s instructions. SMRT sequence and Illumina HiSeq 1500 data have been submitted to the Sequence Read Archive (SRA) of the NCBI database under study accession number SRP094520 of bioproject PRJNA356162.

2.4. Transcriptome Sequencing Data Analysis

Transcriptome assembly was carried out using all clean data in the software Trinity2 [25]. The parameter min kmer cov was set to 2 and all other parameters were set to default. Open reading frame (ORF) prediction of assembled transcripts was performed using TransDecoder (Find Coding Regions within Transcripts). Gene function was annotated using Trinotate2.0 based on Pfam [26], egg-NOG [27], KEGG [28] and GO [29].

2.5. Differential Expression Analysis

Firstly, the read counts were corrected by edgeR [30] program package through one scaling normalized factor. After that, the differential expression analysis of the two samples was performed using the DEGSeq R package (1.12.0) [31]. The Benjamini–Hochberg method was used to adjust the p-value. Corrected p-value of 0.005 and log2 (Fold change) of 1 were set as the threshold for determining whether expression was significantly differential.

2.6. Isoform Detection and Prediction

The quality of short-read libraries was evaluated using the software FastQC 0.11.4 [32], with 13 nucleotides from the 3′ end on each read trimmed. Adapters among short-read libraries were then removed using the software Cutadapt [33], and errors in PacBio long reads were corrected with short-read alignment using the software LSC 1.beta [34]. Following error correction, long and short reads were aligned to the genome using the software HISAT [35], before being input into the software IDP 0.1.9 [23] for the detection and prediction of isoforms.

2.7. Real-Time PCR

Two RNA samples were isolated from D. officinale stems and leaf tissues and reverse transcription was performed using PrimeScriptTM Reverse Transcriptase (TaKaRa, Tokyo, Japan). At the same time, qPCR primers were designed using PRIMER PREMIER 6 (Premier Biosoft International, Palo Alto, CA, USA), and their specificity was verified using PCR. qPCR analysis was then conducted in triplicate using SYBR® Premix Ex TaqTM (TaKaRa, Tokyo, Japan) and a 7500 Real-time PCR system (ABI), with DoActin used as the reference gene [36].

2.8. Phylogenetic Analysis

We downloaded amino acid sequences for 17 AtSWEET and 21 OsSWEET genes from the National Center for Biotechnology Information (NCBI) database, pooled them with sequences for eight DoSWEET genes, and performed an alignment using the software MEGA 5.0 [37]. An unrooted phylogenetic tree for full-length amino acid sequences was constructed using neighbor-joining and the bootstrap method with 1000 replicates.

3. Results

3.1. Comparing the Polysaccharide Content of Stems and Leaves

As a result of their high polysaccharide content, D. officinale plants are used in TCM when they reach heights between 20 cm and 30 cm in the greenhouse (Figure 1a). Indeed, ‘Feng Dou’, the dried stems of D. officinale (Figure 1b), have long been used as crude drugs in TCM as they have good efficacy as a tonic. The leaves and stems of seedlings are separated, and polysaccharide contents is detected respectively (Figure 1c,d); the results have shown that the stems of D. officinale contain twice as much polysaccharide as the leaves, and that concentrations in D. officinale are much higher than in D. fimbriatum at the same growth stage (Figure 1e). These differences mean that D. officinale is well suited for identifying the key genes involved in sugar translocation between tissues via the phloem. Thus, in order to obtain more information on the molecular mechanisms of sugar accumulation in D. officinale, we performed RNA sequencing (RNA-seq) on leaves and stems using SGS and SMRT sequencing platforms.

3.2. Hybrid-Seq and Sequence Assembly

We subjected two mRNA samples from leaves and stems, respectively, to 2 × 150 paired-end sequencing using the HiSeq 1500 platform, generating 68,977,848 reads. Full-length complementary DNA (cDNA) sequences from the same two mRNA samples were then normalized and subjected to SMRT sequencing using the PacBio RS II platform, generating a total of 19,941,155 pre-filter reads. We then carried out post-filtering using the RS_Subreads.1 protocol in the software PacBio RS, generating 18,772,392 subreads. These were then combined with SGS sequences to improve the number of correct subreads. This process generated a total of 197,029 Illumina screened reads which were then used as input data to generate 321,433 corrected reads. The transcript length distribution generated using these two platforms, as well as the quality of PacBio reads, can be seen in Figure 2. Combining high quality SMRT long reads, 1,083,400 corrected reads were then used as input for further analysis; the use of long SMRT sequencing reads (N50 = 1863 base pairs, bp) generated the full-length transcripts relative to assemblies using just SGS reads (SGS-assembled unigene N50 = 658 bp).

3.3. The Reads per Kilobase per Million Mapped Reads Method to Determine Differential Gene Expression between Leaves and Stems

The reads per kilobase per million mapped reads (RPKM) method was used to calculate read densities and to compare differences in gene expression between leaves and stems. SGS data show that the use of a RPKM cut-off value greater than 10 reveals that 1414 genes are differentially expressed between leaves and stems (Table S1), and that of these, 844 and 570 are up-regulated in the two tissue types, respectively. Next, we classified the unigenes involved in biosynthetic pathways using KEGG [38,39,40] metabolic categories to enable a more detailed analysis. The results of this classification showed that 161 differentially expressed genes (DEGs) are involved in metabolic pathways; of these, 49, 6, and 12 DEGs belong to the carbohydrate metabolism, secondary metabolism, and glycan biosynthesis and metabolism pathways, respectively. These include three glucomannan 4-beta-mannosyltransferase (CSLA9) genes (TR50039, TR79639, and TR80798), one GT gene (TR86109), and one Ces A catalytic subunit 6 (CesA6) gene (TR63142), all exhibiting higher levels of expression in stems than in leaves (Figure 3a). The high expression level of these transcripts may be linked to the biosynthesis of polysaccharides in D. officinale. We also identified three GT genes (TR54509, TR80949, and TR82092) that are up-regulated in leaves, and confirmed all expression levels using real-time PCR (RT-PCR) (Figure 3b). In addition, we classified 37 transcription factor (TF) unigenes as well as 13 TF DEGs that are up-regulated in stems. The expression profiles of these TFs that correlated with genes involved in polysaccharide synthesis or sugar transportation will be a focus of future research.

3.4. Candidate Genes for Sugar Translocation in the Stem

Sugar transport proteins play crucial roles in both the cell-to-cell and long-distance movement of sugars throughout plants. In particular, both SWEETs and SUTs are intimately involved in sucrose transport; the transcriptome data presented in this study show that three SWEET genes and one SUT gene are expressed to a greater extent in the stems of D. officinale than in leaves. We confirmed these levels of differential expression using quantitative RT-PCR. The results of this study show that the genes DoSWEET1 and DoSWEET14b are expressed in stems to about twice the level seen in leaves, while stem expression of DoSWEET4 is more than 20-fold higher than in leaves (Figure 4).
Analysis with the software TMHMM 2.0 [41] reveals that the secondary SWEET protein structures of this orchid comprise seven transmembrane helices, while that of SUT contain 12 transmembrane domains (Figure S1). We therefore screened our unigene dataset and, based on annotation results, discovered the presence of five additional D. officinale SWEET genes (Table S2). Amino acid sequence alignments were compared phylogenetically with their homologs in Arabidopsis and rice using the software MEGA 5.0 [37] (Figure 5), and the new genes identified in this study were named DoSWEET1, DoSWEET2, DoSWEET4, DoSWEET14a, DoSWEET14b, DoSWEET15, DoSWEET16, DoSWEET17, and DoSUT1. Results also reveal similar patterns of expression in DoSWEET4 compared to SbSWEET4–3; both these genes are more significantly expressed in stems than they are in leaves during sucrose accumulation [15]. Indeed, our phylogenetic study shows that DoSWEET4 and SbSWEET4–3 cluster together in a single clade with 65.2% identity.

3.5. Alternatively Spliced Isoforms

Alternative splicing has played an extremely important role in understanding protein diversity [42,43]; full-length transcripts can be used to identify alternative splicing events via the hybrid-seq isoform detection and prediction (IDP) pipeline [23]. IDP analysis combined with the D. officinale genome sequence (http://202.203.187.112/genome/dendrobe/) results reveal the presence of 12,910 isoforms in leaves and stem transcriptomes; in total, out of 12,910 multi-exon genes, 2280 (31%) exhibit alternatively spliced isoforms, while two isoforms were found in 19% of cases, three in 7% of cases, and more than four isoforms in 5% of cases (Figure 6). Based on the genome annotation, SMRT sequencing compensates for the defects of SGS assembly strategy in the alternative splicing predicting.
The genes that encode sugar transporters do not seem to undergo alternative splicing to any significant extent, although two GT and four Ces genes, probably involved in polysaccharide metabolism pathways, were identified in the IDP pipeline using a genome-wide approach. Results show that both of these GT genes have two isoforms, while the four Ces genes identified in this study possess five, four, three, and two isoforms, respectively (Figure 7). The alternatively splicing variants possess more than one function. The effect of alternative splicing is not limited to protein isoforms, but extends to gene expression patterns [43].

4. Discussion

Natural plant polysaccharides have an extensive history of application in medicine and pharmaceutics. Indeed, polysaccharide bioactivity is increasingly valued in the context of human health. The orchid D. officinale has long been used as a highly-prized tonic in TCM, and thus has considerable economic value. Polysaccharide is the main active component in D. officinale; the molecular mechanisms underlying bioactive polysaccharide synthesis in this plant have attracted considerable research attention. In earlier work, Zhang and colleagues reported a comparison of transcriptomes from two development stages of D. officinale, generated using an Illumina HiSeq™ 2000 platform [10]. Comparative analysis of the transcriptomes of juvenile seedlings and adult plants revealed that numerous candidate synthases are involved in polysaccharide biosynthesis. Thus, to further address the question of polysaccharide synthesis and sugar translocation from leaves to stems, the short-read SGS and long-read SMRT sequencing of the two distinct tissue transcriptomes available for D. officinale were combined in this study. Our results show that the expression levels of the genes DoSWEET4 and DoSUT1 are markedly increased in stems compared to leaves. Previous work has shown that the gene ZmSUT1 is critical for the efficient loading of sucrose in the phloem in maize leaves [19]. Thus, SWEET4 genes that encode hexose transporters in maize and rice might be also responsible for sugar transport along the phloem unloading path [44]. DoSWEET4 and DoSUT1 potentially act synergistically in the process of sugar translocation from the leaves to the stems when monosaccharides are combined to form polysaccharides. As discussed, while SWEETs comprise a group of newly identified sugar transporters in Arabidopsis and rice [13], this study is the first to report the presence of these transporters in a medicinal plant that is rich in water-soluble polysaccharides.
Sucrose is exported from leaves to stems to sustain plant growth and metabolism; redundant sugars are stored as soluble polysaccharides in the stems. In this respect, D. officinale is a typical mycorrhizal plant, dependent on fungi for seed germination and growth during its seedling and early growth stages. The abundant sugars present in the stem of this plant are therefore likely to sustain the growth of mycorrhizal fungi associated with D. officinale; these fungi obtain energy from plants in exchange for inorganic substances and other nutrient elements. This form of symbiosis is a key environmental adaptation, and SWEETs are likely to play a significant role in providing mycorrhizae with sugar.
In summary, this study reports the functional genomics of a combination of SGS short-read and SMRT long-read transcriptomes of D. officinale. Based on a comparison of leaf and stem transcriptomes, our results reveal the presence of several GT, TF and sugar transporter genes. In addition, the genes encoding GT and Ces undergo a distinct degree of alternative splicing. Both transcriptional and post-transcriptional regulation may be involved in the synthesis of polysaccharides. We also report the differential expression of genes DoSWEET4 and DoSUT1 in the leaves and stems of D. officinale, and note that these sugar transporting genes should be targets for identifying quality polysaccharide production in medicinal plants. The results of this study reveal a new direction for investigations into polysaccharide accumulation in medicinal plants, via the process of sugar translocation in D. officinale.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4425/8/10/257/s1. Figure S1 Prediction of transmembrane helices in SWEETs and SUT from D. officinale. Table S1 Summary of the differentially expressed genes between leaves and stems. Table S2 Summary of putative SWEETs in D. officinale.

Acknowledgments

This work was supported by National Natural Science Foundation of China (NO. 81373916), CAMS Innovation Fund for Medical Sciences (CIFMS) (NO. 2016-I2M-3-016), the Key Projects in the National Science and Technology Pillar Program (NO. 2011BAI07B08), an institutional fund of Department of Internal Medicine, University of Iowa and the National Human Genome Research Institute (R01HG008759). We thank Qingsong Shao from School of Forestry and Bio-technology, Zhejiang A & F University and Yaosheng Wu from Department of Biochemistry and Molecular Biology, Guangxi Medical University for sample collection.

Author Contributions

H.Y. and L.H. designed the experiments. L.H. and H.Z. collected samples. L.H., H.Z. and Z.X. contributed to the RNA extraction, sequencing libraries preparation. J.X. contributed to produce the SGS data. J.Y. contributed to the SGS data analysis. S.F., K.F.A., Z.X. and Y.L. contributed to the SMRT data analysis. Z.X. and X.C. contributed to RT-PCR experiments. Y.H. and J.Z. contributed to detect the content of polysaccharides in Dendrobium. L.H. contributed to the manuscript writing.

Conflicts of Interest

The authors declare no competing financial interests.

Availability of Materials and Data

The datasets generated during the current study are available in the NCBI database (the accession number SRP094520 of bioproject PRJNA356162).

References

  1. Chinese Pharmacopoeia Commission. Pharmacopeia of the People’s Republic of China; China Medical Science Press: Beijing, China, 2015; pp. 92–93. [Google Scholar]
  2. Liu, X.F.; Zhu, J.; Ge, S.Y.; Xia, L.J.; Yang, H.Y.; Qian, Y.T.; Ren, F.Z. Orally administered Dendrobium officinale and its polysaccharides enhance immune functions in BALB/c mice. Nat. Prod. Commun. 2011, 6, 867–870. [Google Scholar] [PubMed]
  3. Ng, T.B.; Liu, J.; Wong, J.H.; Ye, X.; Wing Sze, S.C.; Tong, Y.; Zhang, K.Y. Review of research on Dendrobium, a prized folk medicine. Appl. Microbiol. Biotechnol. 2012, 93, 1795–1803. [Google Scholar] [CrossRef] [PubMed]
  4. Xie, S.Z.; Liu, B.; Zhang, D.D.; Zha, X.Q.; Pan, L.H.; Luo, J.P. Intestinal immunomodulating activity and structural characterization of a new polysaccharide from stems of Dendrobium officinale. Food Funct. 2016, 7, 2789–2799. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, X.J.; Nie, S.P.; Cai, H.L.; Zhang, G.Y.; Cui, S.W.; Xie, M.Y.; Phillips, G.O. Study on Dendrobium officinale O-acetyl-glucomannan (Dendronan): Part, IV. Immunomodulatory activity in vivo. J. Funct. Foods 2015, 15, 525–532. [Google Scholar] [CrossRef]
  6. Hua, Y.F.; Zhang, M.; Fu, C.X.; Chen, Z.H.; Chan, G.Y.S. Structural characterization of a 2-O-acetyl glucomannan from Dendrobium officinale stem. Carbohydr. Res. 2004, 339, 2219–2224. [Google Scholar] [CrossRef] [PubMed]
  7. Luo, A.X.; He, X.J.; Zhou, S.D.; Fan, Y.J.; Luo, A.S.; Chun, Z. Purification, composition analysis and antioxidant activity of the polysaccharides from Dendrobium nobile Lindl. Carbohydr. Polym. 2010, 79, 1014–1019. [Google Scholar] [CrossRef]
  8. Porter, H.K. Synthesis of polysaccharides of higher plants. Ann. Rev. Plant Physiol. 2003, 13, 303–326. [Google Scholar] [CrossRef]
  9. Zhang, G.Q.; Xu, Q.; Bian, C.; Tsai, W.C.; Yeh, C.M.; Liu, K.W.; Yoshida, K.; Zhang, L.S.; Chang, S.B.; Chen, F.; et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Zhang, J.X.; He, C.M.; Wu, K.L.; Teixeira da Silva, J.A.; Zeng, S.J.; Zhang, X.H.; Yu, Z.; Xia, H.; Duan, J. Transcriptome analysis of Dendrobium officinale and its application to the identification of genes associated with polysaccharide synthesis. Front. Plant Sci. 2016, 7, 5. [Google Scholar] [CrossRef] [PubMed]
  11. He, C.M.; Wu, K.L.; Zhang, J.X.; Liu, X.C.; Zeng, S.J.; Yu, Z.M.; Zhang, X.H.; Teixeira da Silva, J.A.; Deng, R.F.; Tan, J.W.; et al. Cytochemical localization of polysaccharides in Dendrobium officinale and the involvement of DoCSLA6 in the synthesis of mannan polysaccharides. Front. Plant Sci. 2017, 8, 173. [Google Scholar] [CrossRef] [PubMed]
  12. Ludewig, F.; Flügge, U.I. Role of metabolite transporters in source-sink carbon allocation. Front. Plant Sci. 2013, 4, 231. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, L.Q.; Hou, B.H.; Lalonde, S.; Takanaga, H.; Hartung, M.L.; Qu, X.Q.; Guo, W.J.; Kim, J.G.; Underwood, W.; Chaudhuri, B.; et al. Sugar transporters for intercellular exchange and nutrition of pathogens. Nature 2010, 468, 527–532. [Google Scholar] [CrossRef] [PubMed]
  14. Feng, C.Y.; Han, J.X.; Han, X.X.; Jiang, J. Genome-wide identification, phylogeny, and expression analysis of the SWEET gene family in tomato. Gene 2015, 573, 261–272. [Google Scholar] [CrossRef] [PubMed]
  15. Mizuno, H.; Kasuga, S.; Kawahigashi, H. The sorghum SWEET gene family: Stem sucrose accumulation as revealed through transcriptome profiling. Biotechnol. Biofuels 2016, 9, 127. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, L.Q.; Qu, X.Q.; Hou, B.H.; Sosso, D.; Osorio, S.; Fernie, A.R.; Frommer, W.B. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science 2012, 335, 207–211. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, L.Q.; Cheung, L.S.; Feng, L.; Tanner, W.; Frommer, W.B. Transport of sugars. Ann. Rev. Biochem. 2015, 84, 865–894. [Google Scholar] [CrossRef] [PubMed]
  18. Antony, G.; Zhou, J.H.; Huang, S.; Li, T.; Liu, B.; White, F.; Yang, B. Rice xa13 recessive resistance to bacterial blight is defeated by induction of the disease susceptibility gene Os-11N3. Plant Cell 2010, 22, 3864–3876. [Google Scholar] [CrossRef] [PubMed]
  19. Slewinski, T.L.; Meeley, R.; Braun, D.M. Sucrose transporter 1 functions in phloem loading in maize leaves. J. Exp. Bot. 2009, 60, 881–892. [Google Scholar] [CrossRef] [PubMed]
  20. Carpaneto, A.; Koepsell, H.; Bamberg, E.; Hedrich, R.; Geiger, D. Sucrose- and H-dependent charge movements associated with the gating of sucrose transporter ZmSUT1. PLoS ONE 2010, 5, e12605. [Google Scholar] [CrossRef] [PubMed]
  21. Yan, L.; Wang, X.; Liu, H.; Tian, Y.; Lian, J.; Yang, R.; Hao, S.; Wang, X.; Yang, S.; Li, Q.; et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese Orchid herb. Mol. Plant 2015, 8, 922–934. [Google Scholar] [CrossRef] [PubMed]
  22. Meng, Y.J.; Yu, D.L.; Xue, J.; Lu, J.J.; Feng, S.G.; Shen, C.J.; Wang, H.Z. A transcriptome-wide, organ-specific regulatory map of Dendrobium officinale, an important traditional Chinese orchid herb. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [PubMed]
  23. Au, K.F.; Sebastiano, V.; Afshar, P.T.; Durruthy, D.J.; Lee, L.; Williams, B.A.; Bakel, H.; Schadt, E.E.; Rejio-Pera, R.A.; Underwood, J.G.; et al. Characterization of the human ESC transcriptome by hybrid sequencing. Proc. Natl. Acad. Sci. USA 2013, 110, E4821–E4830. [Google Scholar] [CrossRef] [PubMed]
  24. Sharon, D.; Tilgner, H.; Grubert, F.; Snyder, M. A single-molecule long-read survey of the human transcriptome. Nat. Biotechnol. 2013, 31, 1009–1014. [Google Scholar] [CrossRef] [PubMed]
  25. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  26. Punta, M.; Coggill, P.C.; Eberhardt, R.Y.; Mistry, J.; Tate, J.; Boursnell, C.; Pang, N.; Forslund, K.; Ceric, G.; Clements, J.; et al. The Pfam protein families database. Nucleic Acids Res. 2012, 40, D290–D301. [Google Scholar] [CrossRef] [PubMed]
  27. Powell, S.; Szklarczyk, D.; Trachana, K.; Roth, A.; Kuhn, M.; Muller, J.; Arnold, R.; Rattei, T.; Letunic, I.; Doerks, T.; et al. eggNOG v3.0: Orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res. 2012, 40, D284–D289. [Google Scholar] [CrossRef] [PubMed]
  28. Kanehisa, M.; Goto, S.; Sato, Y.; Furumichi, M.; Tanabe, M. KEGG for integration and interpretation of large-scale molecular datasets. Nucleic Acids Res. 2012, 40, D109–D114. [Google Scholar] [CrossRef] [PubMed]
  29. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Gene Ontol. Consort. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  30. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  31. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  32. Patel, R.K.; Jain, M. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data. PLoS ONE 2012, 7, e30619. [Google Scholar] [CrossRef] [PubMed]
  33. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  34. Au, K.F.; Underwood, J.G.; Lee, L.; Wong, W.H. Improving PacBio long read accuracy by short read alignment. PLoS ONE 2012, 7, e46679. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed]
  36. Kozera, B.; Rapacz, M. Reference genes in real-time PCR. J. Appl. Genet. 2013, 54, 391–406. [Google Scholar] [CrossRef] [PubMed]
  37. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  38. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed]
  39. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, D457–D462. [Google Scholar] [CrossRef] [PubMed]
  40. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  41. Krogh, A.; Larsson, B.; Von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed]
  42. Graveley, B.R. Alternative splicing: Increasing diversity in the proteomic world. Trends Genet. 2001, 17, 100–107. [Google Scholar] [CrossRef]
  43. Kelemen, O.; Convertini, P.; Zhang, Z.Y.; Wen, Y.; Shen, M.L.; Falaleeva, M.; Stefan, S. Function of alternative splicing. Gene 2013, 514, 1–30. [Google Scholar] [CrossRef] [PubMed]
  44. Sosso, D.; Luo, D.; Li, Q.B.; Sasse, J.; Yang, J.; Gendrot, G.; Suzuki, M.; Koch, K.E.; McCarty, D.R.; Chourey, P.S.; et al. Seed filling in domesticated maize and rice depends on SWEET-mediated hexose transport. Nat. Genet. 2015, 47, 1489–1493. [Google Scholar] [CrossRef] [PubMed]
Genes 08 00257 g001 550
Figure 1. Morphology and chemical profiling of the leaves and stems of Dendrobium officinale. (a) Cultivar morphology in the greenhouse; (b) Morphology of the dry stem used in traditional Chinese medicine (TCM); (c) Stems that were sampled for RNA-seq and chemical profiling; (d) Leaves that were sampled for RNA-seq and chemical profiling; (e) Analyses of the polysaccharide contents of leaves and stems.
Figure 1. Morphology and chemical profiling of the leaves and stems of Dendrobium officinale. (a) Cultivar morphology in the greenhouse; (b) Morphology of the dry stem used in traditional Chinese medicine (TCM); (c) Stems that were sampled for RNA-seq and chemical profiling; (d) Leaves that were sampled for RNA-seq and chemical profiling; (e) Analyses of the polysaccharide contents of leaves and stems.
Genes 08 00257 g001
Genes 08 00257 g002 550
Figure 2. Distribution of D. officinale second-generation sequencing (SGS) and single-molecule real-time (SMRT) reads. (a) Transcript length distributions from the two different sequencing platforms; (b) PacBio quality of subreads and corrected reads.
Figure 2. Distribution of D. officinale second-generation sequencing (SGS) and single-molecule real-time (SMRT) reads. (a) Transcript length distributions from the two different sequencing platforms; (b) PacBio quality of subreads and corrected reads.
Genes 08 00257 g002
Genes 08 00257 g003 550
Figure 3. Comparison of polysaccharide biosynthesis candidate gene expression in leaves and stems of D. officinale. (a) Heat map of candidate gene expression profiles from SGS data; (b) RT-qPCR analysis of candidate genes in leaves and stems. Data are shown as means ± standard error (n = 3).
Figure 3. Comparison of polysaccharide biosynthesis candidate gene expression in leaves and stems of D. officinale. (a) Heat map of candidate gene expression profiles from SGS data; (b) RT-qPCR analysis of candidate genes in leaves and stems. Data are shown as means ± standard error (n = 3).
Genes 08 00257 g003
Genes 08 00257 g004 550
Figure 4. RT-qPCR analysis of DoSWEETs in leaves and stems. Data are shown as means ± standard error (n = 3).
Figure 4. RT-qPCR analysis of DoSWEETs in leaves and stems. Data are shown as means ± standard error (n = 3).
Genes 08 00257 g004
Genes 08 00257 g005 550
Figure 5. Neighbor-joining phylogenetic tree of sugars will eventually be exported transporters (SWEETs) from D. officinale, Arabidopsis, and rice. This analysis comprised amino acid sequences for 17 AtSWEETs, 21 OsSWEETs, and 8 DoSWEETs (in red).
Figure 5. Neighbor-joining phylogenetic tree of sugars will eventually be exported transporters (SWEETs) from D. officinale, Arabidopsis, and rice. This analysis comprised amino acid sequences for 17 AtSWEETs, 21 OsSWEETs, and 8 DoSWEETs (in red).
Genes 08 00257 g005
Genes 08 00257 g006 550
Figure 6. Detection and prediction of D. officinale gene isoforms using isoform detection and prediction (IDP). (a) Venn diagram to show the detection and prediction of 2680 and 10,230 isoforms, respectively; (b) Pie chart to show different alternatively spliced types. ES: exon skipping; IR: intron retention; A3’S: alternative 3′ splice site; A5’S: alternative 5′ splice site; (c) The distribution of alternatively spliced isoforms from each gene locus.
Figure 6. Detection and prediction of D. officinale gene isoforms using isoform detection and prediction (IDP). (a) Venn diagram to show the detection and prediction of 2680 and 10,230 isoforms, respectively; (b) Pie chart to show different alternatively spliced types. ES: exon skipping; IR: intron retention; A3’S: alternative 3′ splice site; A5’S: alternative 5′ splice site; (c) The distribution of alternatively spliced isoforms from each gene locus.
Genes 08 00257 g006
Genes 08 00257 g007 550
Figure 7. Alternative splicing isoforms of genes involved in polysaccharide biosynthesis encompassed by PacBio long reads. GT: glycosyltransferases gene; Ces: cellulose synthase gene.
Figure 7. Alternative splicing isoforms of genes involved in polysaccharide biosynthesis encompassed by PacBio long reads. GT: glycosyltransferases gene; Ces: cellulose synthase gene.
Genes 08 00257 g007

Share and Cite

MDPI and ACS Style

He, L.; Fu, S.; Xu, Z.; Yan, J.; Xu, J.; Zhou, H.; Zhou, J.; Chen, X.; Li, Y.; Au, K.F.; et al. Hybrid Sequencing of Full-Length cDNA Transcripts of Stems and Leaves in Dendrobium officinale. Genes 2017, 8, 257. https://doi.org/10.3390/genes8100257

AMA Style

He L, Fu S, Xu Z, Yan J, Xu J, Zhou H, Zhou J, Chen X, Li Y, Au KF, et al. Hybrid Sequencing of Full-Length cDNA Transcripts of Stems and Leaves in Dendrobium officinale. Genes. 2017; 8(10):257. https://doi.org/10.3390/genes8100257

Chicago/Turabian Style

He, Liu, Shuhua Fu, Zhichao Xu, Jun Yan, Jiang Xu, Hong Zhou, Jianguo Zhou, Xinlian Chen, Ying Li, Kin Fai Au, and et al. 2017. "Hybrid Sequencing of Full-Length cDNA Transcripts of Stems and Leaves in Dendrobium officinale" Genes 8, no. 10: 257. https://doi.org/10.3390/genes8100257

APA Style

He, L., Fu, S., Xu, Z., Yan, J., Xu, J., Zhou, H., Zhou, J., Chen, X., Li, Y., Au, K. F., & Yao, H. (2017). Hybrid Sequencing of Full-Length cDNA Transcripts of Stems and Leaves in Dendrobium officinale. Genes, 8(10), 257. https://doi.org/10.3390/genes8100257

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop