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

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Comparative transcriptional analysis of metabolic pathways and mechanisms regulating essential oil biosynthesis in four elite Cymbopogon spp

Cymbopogon is an important aromatic and medicinal grass with several species of ethnopharmaceutical importance. The genus is extremely rich in secondary metabolites, monoterpenes like geraniol and citral being principal constituents, also used as biomarker for classification and identification of Cymbopogon chemotypes. In the light of this, present study involved RNA sequencing and comparison of expression profiles of four contrasting Cymbopogon species namely C. flexuosus var. Chirharit (citral rich and frost resistant), C. martinii var. PRC-1 (geraniol rich), C. pendulus var. Praman (the most stable and citral-rich genotype), and Jamrosa (a hybrid of C. nardus var. confertiflorus × C. jwarancusa (rich in geraniol and geranyl acetate). The transcriptome profiles revealed marked differences in gene expression patterns of 28 differentially expressed genes (DEGs) of terpenoid metabolic pathways between the four Cymbopogon sp. The major DEGs were Carotenoid Cleavage Dioxygenases (CCD), Aspartate aminotransferase (ASP amino), Mevalonate E-4 hydroxy, AKR, GGPS, FDPS, and AAT. In addition, few TFs related to different regulatory pathways were also identified. The gene expression profiles of DEGs were correlated to the EO yield and their monoterpene compositions. Overall, the PRC-1 (C. martinii) shows distinguished gene expression profiles from all other genotypes. Thus, the transcriptome sequence database expanded our understanding of terpenoid metabolism and its molecular regulation in Cymbopogon species. Additionally, this data also serves as an important source of knowledge for enhancing oil yield and quality in Cymbopogon and closely related taxa. KEY MESSAGE: Unfolding the new secretes surrounding EO biosynthesis and regulation in four contrasting Cymbopogon species.

Full-length transcriptome profiling reveals insight into the cold response of two kiwifruit genotypes (A. arguta) with contrasting freezing tolerances

BACKGROUND

Kiwifruit (Actinidia Lindl.) is considered an important fruit species worldwide. Due to its temperate origin, this species is highly vulnerable to freezing injury while under low-temperature stress. To obtain further knowledge of the mechanism underlying freezing tolerance, we carried out a hybrid transcriptome analysis of two A. arguta (Actinidi arguta) genotypes, KL and RB, whose freezing tolerance is high and low, respectively. Both genotypes were subjected to - 25 °C for 0 h, 1 h, and 4 h.

RESULTS

SMRT (single-molecule real-time) RNA-seq data were assembled using the de novo method, producing 24,306 unigenes with an N50 value of 1834 bp. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs showed that they were involved in the 'starch and sucrose metabolism', the 'mitogen-activated protein kinase (MAPK) signaling pathway', the 'phosphatidylinositol signaling system', the 'inositol phosphate metabolism', and the 'plant hormone signal transduction'. In particular, for 'starch and sucrose metabolism', we identified 3 key genes involved in cellulose degradation, trehalose synthesis, and starch degradation processes. Moreover, the activities of beta-GC (beta-glucosidase), TPS (trehalose-6-phosphate synthase), and BAM (beta-amylase), encoded by the abovementioned 3 key genes, were enhanced by cold stress. Three transcription factors (TFs) belonging to the AP2/ERF, bHLH (basic helix-loop-helix), and MYB families were involved in the low-temperature response. Furthermore, weighted gene coexpression network analysis (WGCNA) indicated that beta-GC, TPS5, and BAM3.1 were the key genes involved in the cold response and were highly coexpressed together with the CBF3, MYC2, and MYB44 genes.

CONCLUSIONS

Cold stress led various changes in kiwifruit, the 'phosphatidylinositol signaling system', 'inositol phosphate metabolism', 'MAPK signaling pathway', 'plant hormone signal transduction', and 'starch and sucrose metabolism' processes were significantly affected by low temperature. Moreover, starch and sucrose metabolism may be the key pathway for tolerant kiwifruit to resist low temperature damages. These results increase our understanding of the complex mechanisms involved in the freezing tolerance of kiwifruit under cold stress and reveal a series of candidate genes for use in breeding new cultivars with enhanced freezing tolerance.

Transcriptomics and Metabolomics Reveal Purine and Phenylpropanoid Metabolism Response to Drought Stress in Dendrobium sinense, an Endemic Orchid Species in Hainan Island

Drought stress is a bottleneck factor for plant growth and development, especially in epiphytic orchids that absorb moisture mainly from the air. Recent studies have suggested that there are complex transcriptional regulatory networks related to drought stress in Dendrobium sinense. In this study, the transcription and metabolite alterations involved in drought stress response in D. sinense were investigated through RNA-seq and metabolomics. A total of 856 metabolites were identified from stressed and control samples, with 391 metabolites showing significant differences. With PacBio and Illumina RNA sequencing, 72,969 genes were obtained with a mean length of 2,486 bp, and 622 differentially expressed genes (DEGs) were identified. Correlation analysis showed 7 differential genes, and 39 differential metabolites were involved in interaction networks. The network analysis of differential genes and metabolites suggested that the pathways of purine metabolism and phenylpropanoid biosynthesis may play an important role in drought response in D. sinense. These results provide new insights and reference data for culturally important medicinal plants and the protection of endangered orchids.

Isolation, structural properties, bioactivities of polysaccharides from Dendrobium officinale Kimura et. Migo: A review

Dendrobium officinale Kimura et Migo (D. officinale) is used as herbal medicine and new food resource in China, which is nontoxic and harmless, and can be used as common food. Polysaccharide as one of the main bioactive components in D. officinale, mainly composed of glucose and mannose (Manp: Glcp = 2.01:1.00-8.82:1.00), along with galactose, xylose, arabinose, and rhamnose in different molar ratios and types of glycosidic bonds. Polysaccharides of D. officinale exhibit a variety of biological effects, including immunomodulatory, anti-tumor, gastro-protective, hypoglycemic, anti-inflammatory, hepatoprotective, and vasodilating effects. This paper presents the extraction, purification, structural characteristics, bioactivities, structure-activity relationships and analyzes gaps in the current research on D. officinale polysaccharides. In addition, based on in vitro and in vivo experiments, the possible mechanisms of bioactivities of D. officinale polysaccharides were summarized. We hope that this work may provide helpful references and promising directions for further study and development of D. officinale polysaccharides.

Comparative physiological and full-length transcriptome analyses reveal the molecular mechanism of melatonin-mediated salt tolerance in okra (Abelmoschus esculentus L.)

BACKGROUND

Melatonin, a multifunctional signal molecule, has been reported to play crucial roles in growth and development and stress responses in various plant species. Okra (Abelmoschus esculentus L.) is a food crop with extremely high values of nutrition and healthcare. Recent reports have revealed the protective role of melatonin in alleviating salt stress. However, little is known about its regulatory mechanisms in response to salt stress in okra.

RESULTS

In this study, we explored whether exogenous melatonin pretreatment could alleviate salt stress (300 mM NaCl) of okra plants. Results showed that exogenous application of melatonin (50 μM) significantly enhanced plant tolerance to salt stress, as demonstrated by the plant resistant phenotype, as well as by the higher levels of the net photosynthetic rate, chlorophyll fluorescence and chlorophyll content in comparison with nontreated salt-stressed plants. Additionally, melatonin pretreatment remarkably decreased the levels of lipid peroxidation and H₂O₂ content and scavenged O₂^•- in melatonin-pretreated plants, which may be attributed to the higher levels of enzyme activities including POD and GR. Moreover, a combination of third- (PacBio) and second-generation (Illumina) sequencing technologies was applied to sequence full-length transcriptomes of okra. A total of 121,360 unigenes was obtained, and the size of transcript lengths ranged from 500 to 6000 bp. Illumina RNA-seq analysis showed that: Comparing with control, 1776, 1063 and 1074 differential expression genes (DEGs) were identified from the three treatments (NaCl, MT50 and MT + NaCl, respectively). These genes were enriched in more than 10 GO terms and 34 KEGG pathways. Nitrogen metabolism, sulfur metabolism, and alanine, aspartate and glutamate metabolism were significantly enriched in all three treatments. Many transcription factors including MYB, WRKY, NAC etc., were also identified as DEGs.

CONCLUSIONS

Our preliminary results suggested that melatonin pretreatment enhanced salt tolerance of okra plants for the first time. These data provide the first set of full-length isoforms in okra and more comprehensive insights into the molecular mechanism of melatonin responses to salt stress.

Research Advances in Multi-Omics on the Traditional Chinese Herb Dendrobium officinale

Dendrobium officinale Kimura et Migo is an important epiphytic plant, belonging to the Orchidaceae family. There are various bioactive components in D. officinale plants, mainly including polysaccharides, alkaloids, and phenolic compounds. These compounds have been demonstrated to possess multiple functions, such as anti-oxidation, immune regulation, and anti-cancer. Due to serious shortages of wild resources, deterioration of cultivated germplasm and the unstable quality of D. officinale, the study has been focused on the biosynthetic pathway and regulation mechanisms of bioactive compounds. In recent years, with rapid developments in detection technologies and analysis tools, omics research including genomics, transcriptomics, proteomics and metabolomics have all been widely applied in various medicinal plants, including D. officinale. Many important advances have been achieved in D. officinale research, such as chromosome-level reference genome assembly and the identification of key genes involved in the biosynthesis of active components. In this review, we summarize the latest research advances in D. officinale based on multiple omics studies. At the same time, we discuss limitations of the current research. Finally, we put forward prospective topics in need of further study on D. officinale.

Full-length transcriptome sequencing analysis and development of EST-SSR markers for the endangered species Populus wulianensis

Populus wulianensis is an endangered species endemic to Shandong Province, China. Despite the economic and ornamental value of this species, few genomics and genetic studies have been performed. In this study, we performed a relevant analysis of the full-length transcriptome sequencing data of P. wulianensis and obtained expressed sequence tag (EST)-simple sequence repeat (SSR) markers with polymorphisms that can be used for further genetic research. In total, 8.18 Gb (3,521,665) clean reads with an average GC content of 42.12% were obtained. From the corrected 64,737 high-quality isoforms, 42,323 transcript sequences were obtained after redundancy analysis with CD-HIT. Among these transcript sequences, 41,876 sequences were annotated successfully. A total of 23,539 potential EST-SSRs were identified from 16,057 sequences. Excluding mononucleotides, the most abundant motifs were trinucleotide SSRs (47.80%), followed by di- (46.80%), tetra- (2.98%), hexa- (1.58%) and pentanucleotide SSRs (0.84%). Among the 100 designed EST-SSRs, 18 were polymorphic with high PIC values (0.721 and 0.683) and could be used for analyses of the genetic diversity and population structure of P. wulianensis. These full-length transcriptome sequencing data will facilitate gene discovery and functional genomics research in P. wulianensis, and the novel EST-SSRs developed in our study will promote molecular-assisted breeding, genetic diversity and conservation biology research in this species.

Multi-omics sequencing provides insight into floral transition in Catalpa bungei. C.A. Mey

BACKGROUND

Floral transition plays an important role in development, and proper time is necessary to improve the value of valuable ornamental trees. The molecular mechanisms of floral transition remain unknown in perennial woody plants. "Bairihua" is a type of C. bungei that can undergo floral transition in the first planting year.

RESULTS

Here, we combined short-read next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to provide a more complete view of transcriptome regulation during floral transition in C. bungei. The circadian rhythm-plant pathway may be the critical pathway during floral transition in early flowering (EF) C. bungei, according to horizontal and vertical analysis in EF and normal flowering (NF) C. bungei. SBP and MIKC-MADS-box were seemingly involved in EF during floral transition. A total of 61 hub genes were associated with floral transition in the MEturquoise model with Weighted Gene Co-expression Network Analysis (WGCNA). The results reveal that ten hub genes had a close connection with the GASA homologue gene (Cbu.gene.18280), and the ten co-expressed genes belong to five flowering-related pathways. Furthermore, our study provides new insights into the complexity and regulation of alternative splicing (AS). The ratio or number of isoforms of some floral transition-related genes is different in different periods or in different sub-genomes.

CONCLUSIONS

Our results will be a useful reference for the study of floral transition in other perennial woody plants. Further molecular investigations are needed to verify our sequencing data.

Transcriptomic Landscape of Medicinal Dendrobium Reveals Genes Associated With the Biosynthesis of Bioactive Components

Many plants of Dendrobium genus are precious traditional herbs with high commercial value and excellent medicinal effects. They are perennial aerophytes or epiphytes of terrestrial orchids growing on cliffs and tree trunks covered with mosses in forests throughout the tropical and subtropical Asia and eastern Australia. The stem contains a variety of bioactive components, including polysaccharides and alkaloids, with strong antioxidant, neuroprotective, and immunomodulatory effects. Great attention has been drawn to the Dendrobium genus regarding its medicinal effectiveness, and the related researches have been accumulating rapidly in recent years. The bioactive components are mainly the intermediates or final products produced in specialized metabolite biosynthesis. Thus far, the activity, molecular structure, and composition of major medicinal ingredients have been partially elucidated, and the sequencing of several transcriptomes has been starting to shed new light on the biosynthesis regulation mechanism. This paper reviewed the advances of researches concerning the biosynthetic pathways of medicinal specialized metabolites from Dendrobium, especially the large number of related genes, with the hope of further promoting the development and utilization of those components and correspondingly protecting the Dendrobium resources in more effective ways.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Hybrid Sequencing of Full-Length cDNA Transcripts of the Medicinal Plant Scutellaria baicalensis

Scutellaria baicalensis is a well-known medicinal plant that produces biologically active flavonoids, such as baicalin, baicalein, and wogonin. Pharmacological studies have shown that these compounds have anti-inflammatory, anti-bacterial, and anti-cancer activities. Therefore, it is of great significance to investigate the genetic information of S. baicalensis, particularly the genes related to the biosynthetic pathways of these compounds. Here, we constructed the full-length transcriptome of S. baicalensis using a hybrid sequencing strategy and acquired 338,136 full-length sequences, accounting for 93.3% of the total reads. After the removal of redundancy and correction with Illumina short reads, 75,785 nonredundant transcripts were generated, among which approximately 98% were annotated with significant hits in the protein databases, and 11,135 sequences were classified as lncRNAs. Differentially expressed gene (DEG) analysis showed that most of the genes related to flavonoid biosynthesis were highly expressed in the roots, consistent with previous reports that the flavonoids were mainly synthesized and accumulated in the roots of S. baicalensis. By constructing unique transcription models, a total of 44,071 alternative splicing (AS) events were identified, with intron retention (IR) accounting for the highest proportion (44.5%). A total of 94 AS events were present in five key genes related to flavonoid biosynthesis, suggesting that AS may play important roles in the regulation of flavonoid biosynthesis in S. baicalensis. This study provided a large number of highly accurate full-length transcripts, which represents a valuable genetic resource for further research of the molecular biology of S. baicalensis, such as the development, breeding, and biosynthesis of active ingredients.

Comprehensive transcriptome analysis reveals genes potentially involved in isoflavone biosynthesis in Pueraria thomsonii Benth

Pueraria thomsonii Benth is an important medicinal plant. Transcriptome sequencing, unigene assembly, the annotation of transcripts and the study of gene expression profiles play vital roles in gene function research. However, the full-length transcriptome of P. thomsonii remains unknown. Here, we obtained 44,339 nonredundant transcripts of P. thomsonii by using the PacBio RS II Isoform and Illumina sequencing platforms, of which 43,195 were annotated genes. Compared with the expression levels in the plant roots, those of transcripts with a |fold change| ≥ 4 and FDR < 0.01 in the leaves or stems were assigned as differentially expressed transcripts (DETs). In total, we found 9,225 DETs, 32 of which came from structural genes that were potentially involved in isoflavone biosynthesis. The expression profiles of 8 structural genes from the RNA-Seq data were validated by qRT-PCR. We identified 437 transcription factors (TFs) that were positively or negatively correlated with at least 1 of the structural genes involved in isoflavone biosynthesis using Pearson correlation coefficients (r) (r > 0.8 or r < -0.8). We also identified a total of 32 microRNAs (miRNAs), which targeted 805 transcripts. These miRNAs caused enriched function in ‘ATP binding’, ‘defense response’, ‘ADP binding’, and ‘signal transduction’. Interestingly, MIR156a potentially promoted isoflavone biosynthesis by repressing SBP, and MIR319 promoted isoflavone biosynthesis by repressing TCP and HB-HD-ZIP. Finally, we identified 2,690 alternative splicing events, including that of the structural genes of trans-cinnamate 4-monooxygenase and pullulanase, which are potentially involved in the biosynthesis of isoflavone and starch, respectively, and of three TFs potentially involved in isoflavone biosynthesis. Together, these results provide us with comprehensive insight into the gene expression and regulation of P. thomsonii.

Full-Length Transcriptome Analysis of the Genes Involved in Tocopherol Biosynthesis in Torreya grandis

The seeds of Torreya grandis (Cephalotaxaceae) are rich in tocopherols, which are essential components of the human diet as a result of their function in scavenging reactive oxygen and free radicals. Different T. grandis cultivars (10 cultivars selected in this study were researched, and their information is shown in Table S1 of the Supporting Information) vary enormously in their tocopherol contents (0.28-11.98 mg/100 g). However, little is known about the molecular basis and regulatory mechanisms of tocopherol biosynthesis in T. grandis kernels. Here, we applied single-molecule real-time (SMRT) sequencing to T. grandis (X08 cultivar) for the first time and obtained a total of 97 211 full-length transcripts. We proposed the biosynthetic pathway of tocopherol and identified eight full-length transcripts encoding enzymes potentially involved in tocopherol biosynthesis in T. grandis. The results of the correlation analysis between the tocopherol content and gene expression level in the 10 selected cultivars and different kernel developmental stages of the X08 cultivar suggested that homogentisate phytyltransferase coding gene ( TgVTE2b) and γ-tocopherol methyltransferase coding gene ( TgVTE4) may be key players in tocopherol accumulation in the kernels of T. grandis. Subcellular localization assays showed that both TgVTE2b and TgVTE4 were localized to the chloroplast. We also identified candidate regulatory genes similar to WRI1 and DGAT1 in Arabidopsis that may be involved in the regulation of tocopherol biosynthesis. Our findings provide valuable genetic information for T. grandis using full-length transcriptomic analysis, elucidating the candidate genes and key regulatory genes involved in tocopherol biosynthesis. This information will be critical for further molecular-assisted screening and breeding of T. grandis genotypes with high tocopherol contents.

Genome-wide researches and applications on Dendrobium

This review summarizes current knowledge of chromosome characterization, genetic mapping, genomic sequencing, quality formation, floral transition, propagation, and identification in Dendrobium. The widely distributed Dendrobium has been studied for a long history, due to its important economic values in both medicine and ornamental. In recent years, some species of Dendrobium and other orchids had been reported on genomic sequences, using the next-generation sequencing technology. And the chloroplast genomes of many Dendrobium species were also revealed. The chromosomes of most Dendrobium species belong to mini-chromosomes, and showed 2n = 38. Only a few of genetic studies were reported in Dendrobium. After revealing of genomic sequences, the techniques of transcriptomics, proteomics and metabolomics could be employed on Dendrobium easily. Some other molecular biological techniques, such as gene cloning, gene editing, genetic transformation and molecular marker developing, had also been applied on the basic research of Dendrobium, successively. As medicinal plants, insights into the biosynthesis of some medicinal components were the most important. As ornamental plants, regulation of flower related characteristics was the most important. More, knowledge of growth and development, environmental interaction, evolutionary analysis, breeding of new cultivars, propagation, and identification of species and herbs were also required for commercial usage. All of these studies were improved using genomic sequences and related technologies. To answer some key scientific issues in Dendrobium, quality formation, flowering, self-incompatibility and seed germination would be the focus of future research. And genome related technologies and studies would be helpful.

Full-Length Transcriptome Survey and Expression Analysis of Cassia obtusifolia to Discover Putative Genes Related to Aurantio-Obtusin Biosynthesis, Seed Formation and Development, and Stress Response

The seed is the pharmaceutical and breeding organ of Cassia obtusifolia, a well-known medical herb containing aurantio-obtusin (a kind of anthraquinone), food, and landscape. In order to understand the molecular mechanism of the biosynthesis of aurantio-obtusin, seed formation and development, and stress response of C. obtusifolia, it is necessary to understand the genomics information. Although previous seed transcriptome of C. obtusifolia has been carried out by short-read next-generation sequencing (NGS) technology, the vast majority of the resulting unigenes did not represent full-length cDNA sequences and supply enough gene expression profile information of the various organs or tissues. In this study, fifteen cDNA libraries, which were constructed from the seed, root, stem, leaf, and flower (three repetitions with each organ) of C. obtusifolia, were sequenced using hybrid approach combining single-molecule real-time (SMRT) and NGS platform. More than 4,315,774 long reads with 9.66 Gb sequencing data and 361,427,021 short reads with 108.13 Gb sequencing data were generated by SMRT and NGS platform, respectively. 67,222 consensus isoforms were clustered from the reads and 81.73% (61,016) of which were longer than 1000 bp. Furthermore, the 67,222 consensus isoforms represented 58,106 nonredundant transcripts, 98.25% (57,092) of which were annotated and 25,573 of which were assigned to specific metabolic pathways by KEGG. CoDXS and CoDXR genes were directly used for functional characterization to validate the accuracy of sequences obtained from transcriptome. A total of 658 seed-specific transcripts indicated their special roles in physiological processes in seed. Analysis of transcripts which were involved in the early stage of anthraquinone biosynthesis suggested that the aurantio-obtusin in C. obtusifolia was mainly generated from isochorismate and Mevalonate/methylerythritol phosphate (MVA/MEP) pathway, and three reactions catalyzed by Menaquinone-specific isochorismate synthase (ICS), 1-deoxy-d-xylulose-5-phosphate synthase (DXS) and isopentenyl diphosphate (IPPS) might be the limited steps. Several seed-specific CYPs, SAM-dependent methyltransferase, and UDP-glycosyltransferase (UDPG) supplied promising candidate genes in the late stage of anthraquinone biosynthesis. In addition, four seed-specific transcriptional factors including three MYB Transcription Factor (MYB) and one MADS-box Transcription Factor (MADS) transcriptional factors) and alternative splicing might be involved with seed formation and development. Meanwhile, most members of Hsp20 genes showed high expression level in seed and flower; seven of which might have chaperon activities under various abiotic stresses. Finally, the expressional patterns of genes with particular interests showed similar trends in both transcriptome assay and qRT-PCR. In conclusion, this is the first full-length transcriptome sequencing reported in Caesalpiniaceae family, and thus providing a more complete insight into aurantio-obtusin biosynthesis, seed formation and development, and stress response as well in C. obtusifolia.