Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues of Salvia miltiorrhiza and tanshinone biosynthesis

The methyl jasmonate-responsive transcription factor SmERF106 promotes tanshinone accumulation in Salvia miltiorrhiza

Understanding the regulatory biosynthesis mechanisms of active compounds in herbs is vital for the preservation and sustainable use of natural medicine resources. Diterpenoids, which play a key role in plant growth and resistance, also serve as practical products for humans. Tanshinone, a class of abietane-type diterpenes unique to the Salvia genus, such as Salvia miltiorrhiza, is an excellent model for studying diterpenoids. In this study, we discovered that a transcription factor, SmERF106, responds to MeJA induction and is located in the nucleus. It exhibits a positive correlation with the expression of SmKSL1 and SmIDI1, which are associated with tanshinone biosynthesis. We performed DNA affinity purification sequencing (DAP-seq) to predict genes that may be transcriptionally regulated by SmERF106. Our cis-elements analysis suggested that SmERF106 might bind to GCC-boxes in the promoters of SmKSL1 and SmIDI1. This indicates that SmKSL1 and SmIDI1 could be potential target genes regulated by SmERF106 in the tanshinone biosynthesis pathway. Their interaction was then demonstrated through a series of in vitro and in vivo binding experiments, including Y1H, EMSA, and Dual-LUC. Overexpression of SmERF106 in the hairy root of S. miltiorrhiza led to a significant increase in tanshinone content and the transcriptional levels of SmKSL1 and SmIDI1. In summary, we found that SmERF106 can activate the transcription of SmKSL1 and SmIDI1 in response to MeJA induction, thereby promoting tanshinone biosynthesis. This discovery provides new insights into the regulatory mechanisms of tanshinones in response to JA and offers a potential gene tool for tanshinone metabolic engineering strategy.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.)

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

Biosynthesis and signal transduction of plant growth regulators and their effects on bioactive compound production in Salvia miltiorrhiza (Danshen)

Plant growth regulators (PGRs) are involved in multiple aspects of plant life, including plant growth, development, and response to environmental stimuli. They are also vital for the formation of secondary metabolites in various plants. Salvia miltiorrhiza is a famous herbal medicine and has been used commonly for > 2000 years in China, as well as widely used in many other countries. S. miltiorrhiza is extensively used to treat cardiovascular and cerebrovascular diseases in clinical practices and has specific merit against various diseases. Owing to its outstanding medicinal and commercial potential, S. miltiorrhiza has been extensively investigated as an ideal model system for medicinal plant biology. Tanshinones and phenolic acids are primary pharmacological constituents of S. miltiorrhiza. As the growing market for S. miltiorrhiza, the enhancement of its bioactive compounds has become a research hotspot. S. miltiorrhiza exhibits a significant response to various PGRs in the production of phenolic acids and tanshinones. Here, we briefly review the biosynthesis and signal transduction of PGRs in plants. The effects and mechanisms of PGRs on bioactive compound production in S. miltiorrhiza are systematically summarized and future research is discussed. This article provides a scientific basis for further research, cultivation, and metabolic engineering in S. miltiorrhiza.

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome

BACKGROUND

The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms.

RESULTS

A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing.

CONCLUSION

In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

Genome-Wide Identification and Characterization of the PHT1 Gene Family and Its Response to Mycorrhizal Symbiosis in Salvia miltiorrhiza under Phosphate Stress

Phosphorus (P) is a vital nutrient element that is essential for plant growth and development, and arbuscular mycorrhizal fungi (AMF) can significantly enhance P absorption. The phosphate transporter protein 1 (PHT1) family mediates the uptake of P in plants. However, the PHT1 gene has not yet been characterized in Salvia miltiorrhiza. In this study, to gain insight into the functional divergence of PHT1 genes, nine SmPHT1 genes were identified in the S. miltiorrhiza genome database via bioinformatics tools. Phylogenetic analysis revealed that the PHT1 proteins of S. miltiorrhiza, Arabidopsis thaliana, and Oryza sativa could be divided into three groups. PHT1 in the same clade has a similar gene structure and motif, suggesting that the features of each clade are relatively conserved. Further tissue expression analysis revealed that SmPHT1 was expressed mainly in the roots and stems. In addition, phenotypic changes, P content, and PHT1 gene expression were analyzed in S. miltiorrhiza plants inoculated with AMF under different P conditions (0 mM, 0.1 mM, and 10 mM). P stress and AMF significantly affected the growth and P accumulation of S. miltiorrhiza. SmPHT1;6 was strongly expressed in the roots colonized by AMF, implying that SmPHT1;6 was a specific AMF-inducible PHT1. Taken together, these results provide new insights into the functional divergence and genetic redundancy of the PHT1 genes in response to P stress and AMF symbiosis in S. miltiorrhiza.

Advances in pharmacology, biosynthesis, and metabolic engineering of Scutellaria-specialized metabolites

Scutellaria Linn., which belongs to the family Lamiaceae, is a commonly used medicinal plant for heat clearing and detoxification. In particular, the roots of S. baicalensis and the entire herb of S. barbata have been widely used in traditional medicine for thousands of years. The main active components of Scutellaria, including: baicalein, wogonin, norwogonin, scutellarein, and their glycosides have potential or existing drug usage. However, the wild resources of Scutellaria plants have been overexploited, and degenerated germplasm resources cannot fulfill the requirements of chemical extraction and clinical usage. Metabolic engineering and green production via microorganisms provide alternative strategies for greater efficiency in the production of natural products. Here, we review the progress of: pharmacological investigations, multi-omics, biosynthetic pathways, and metabolic engineering of various Scutellaria species and their active compounds. In addition, based on multi-omics data, we systematically analyze the phylogenetic relationships of Scutellaria and predict candidate transcription factors related to the regulation of active flavonoids. Finally, we propose the prospects of directed evolution of core enzymes and genome-assisted breeding to alleviate the shortage of plant resources of Scutellaria. This review provides important insights into the sustainable utilization and development of Scutellaria resources.

Full-Length Transcriptome Sequencing and RNA-Seq Analysis Offer Insights into Terpenoid Biosynthesis in Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., an important economic and medicinal herb, has a long history of being used as a traditional Chinese medicine. Its leaves have always been used as a raw material for the extraction of essential oils, comprising large amounts of terpenoids, which have good therapeutic effects on many diseases, such as eczema, bacterial infection, and hypertension. However, the genetic basis of terpenoid biosynthesis in this plant is virtually unknown on account of the lack of genomic data. Here, a combination of next-generation sequencing (NGS) and full-length transcriptome sequencing was applied to identify genes involved in terpenoid biosynthesis at five developmental stages. Then, the main components of essential oils in B. balsamifera were identified using GC-MS. Overall, 16 monoterpenoids and 20 sesquiterpenoids were identified and 333,860 CCS reads were generated, yielding 65,045 non-redundant transcripts. Among these highly accurate transcripts, 59,958 (92.18%) transcripts were successfully annotated using NR, eggNOG, Swissprot, KEGG, KOG, COG, Pfam, and GO databases. Finally, a total of 56 differently expressed genes (DEGs) involved in terpenoid biosynthesis were identified, including 38 terpenoid backbone genes and 18 TPSs, which provide a significant amount of genetic information for B. balsamifera. These results build a basis for resource protection, molecular breeding, and the metabolic engineering of this plant.

Comparative transcriptomics of two Salvia subg. Perovskia species contribute towards molecular background of abietane-type diterpenoid biosynthesis

Tanshinones, are a group of diterpenoid red pigments present in Danshen - an important herbal drug of Traditional Chinese Medicine which is a dried root of Salvia miltiorrhiza Bunge. Some of the tanshinones are sought after as pharmacologically active natural products. To date, the biosynthetic pathway of tanshinones has been only partially elucidated. These compounds are also present in some of the other Salvia species, i.a. from subgenus Perovskia, such as S. abrotanoides (Kar.) Sytsma and S. yangii B.T. Drew. Despite of the close genetic relationship between these species, significant qualitative differences in their diterpenoid profile have been discovered. In this work, we have used the Liquid Chromatography-Mass Spectrometry analysis to follow the content of diterpenoids during the vegetation season, which confirmed our previous observations of a diverse diterpenoid profile. As metabolic differences are reflected in different transcript profile of a species or tissues, we used metabolomics-guided transcriptomic approach to select candidate genes, which expression possibly led to observed chemical differences. Using an RNA-sequencing technology we have sequenced and de novo assembled transcriptomes of leaves and roots of S. abrotanoides and S. yangii. As a result, 134,443 transcripts were annotated by UniProt and 56,693 of them were assigned as Viridiplantae. In order to seek for differences, the differential expression analysis was performed, which revealed that 463, 362, 922 and 835 genes indicated changes in expression in four comparisons. GO enrichment analysis and KEGG functional analysis of selected DEGs were performed. The homology and expression of two gene families, associated with downstream steps of tanshinone and carnosic acid biosynthesis were studied, namely: cytochromes P-450 and 2-oxoglutarate-dependend dioxygenases. Additionally, BLAST analysis revealed existence of 39 different transcripts related to abietane diterpenoid biosynthesis in transcriptomes of S. abrotanoides and S. yangii. We have used quantitative real-time RT-PCR analysis of selected candidate genes, to follow their expression levels over the vegetative season. A hypothesis of an existence of a multifunctional CYP76AH89 in transcriptomes of S. abrotanoides and S. yangii is discussed and potential roles of other CYP450 homologs are speculated. By using the comparative transcriptomic approach, we have generated a dataset of candidate genes which provides a valuable resource for further elucidation of tanshinone biosynthesis. In a long run, our investigation may lead to optimization of diterpenoid profile in S. abrotanoides and S. yangii, which may become an alternative source of tanshinones for further research on their bioactivity and pharmacological therapy.

Comparative analysis of three SmHPPR genes encoding hydroxyphenylpyruvate reductases in Salvia miltiorrhiza

Hydroxyphenylpyruvate reductase (HPPR) is an enzyme that is involved in the biosynthesis of hydrophilic phenolic acids in Salvia miltiorrhiza, which is a model medicinal plant. Three SmHPPR genes have been identified in the S. miltiorrhiza genome; however, only one has been functionally analyzed. Here, we cloned three SmHPPR genes (SmHPPR1, SmHPPR2, and SmHPPR3) from the cDNA of S. miltiorrhiza, and their expression profiles were studied. The expression levels of SmHPPR1 were significantly higher than those of SmHPPR2 and SmHPPR3, where SmHPPR1 revealed the highest level in stems, while SmHPPR2 and SmHPPR3 exhibited the highest level in flowers. SmHPPR1, SmHPPR2, and SmHPPR3 are localized in the cytoplasm. All three recombinant enzymes had HPPR activities and catalyzed the reduction of 4-hydroxyphenylpyruvic acid (pHPP) to 4-hydroxyphenyllactic acid (pHPL), with SmHPPR1 showing the highest activity. The transient over-expression of SmHPPR1, SmHPPR2, and SmHPPR3 in the leaves of Nicotiana benthamiana promoted the production of pHPL, which indicated that all three SmHPPRs had in vivo activities. Overall, between the three homologs, SmHPPR1 plays a dominant role in catalyzing pHPP to pHPL, which provides new insights into the biosynthesis of phenolic acids in S. miltiorrhiza.

Multi-omics landscape to decrypt the distinct flavonoid biosynthesis of Scutellaria baicalensis across multiple tissues

Scutellaria baicalensis Georgi, also known as huang-qin in traditional Chinese medicine, is a widely used herbal remedy due to its anticancer, antivirus, and hepatoprotective properties. The S. baicalensis genome was sequenced many years ago; by contrast, the proteome as the executer of most biological processes of S. baicalensis in the aerial parts, as well as the secondary structure of the roots (xylem, phloem, and periderm), is far less comprehensively characterized. Here we attempt to depict the molecular landscape of the non-model plant S. baicalensis through a multi-omics approach, with the goal of constructing a highly informative and valuable reference dataset. Furthermore, we provide an in-depth characterization dissection to explain the two distinct flavonoid biosynthesis pathways that exist in the aerial parts and root, at the protein and phosphorylated protein levels. Our study provides detailed spatial proteomic and phosphoproteomic information in the context of secondary structures, with implications for the molecular profiling of secondary metabolite biosynthesis in non-model medicinal plants.

A comparative full-length transcriptomic resource provides insight into the perennial monocarpic mass flowering

Perennial monocarpic mass flowering represents as a key developmental innovation in flowering time diversity in several biological and economical essential families, such as the woody bamboos and the shrubby Strobilanthes. However, molecular and genetic mechanisms underlying this important biodiversity remain poorly investigated. Here, we generated a full-length transcriptome resource incorporated into the BlueOmics database (http://blueomics.iflora.cn) for two Strobilanthes species, which feature contrasting flowering time behaviors. Using about 112 and 104 Gb Iso-seq reads together with ~185 and ~75 Gb strand-specific RNA seq data, we annotated 80 971 and 79 985 non-redundant full-length transcripts for the perennial polycarpic Strobilanthes tetrasperma and the perennial monocarpic Strobilanthes biocullata, respectively. In S. tetrasperma, we identified 8794 transcripts showing spatiotemporal expression in nine tissues. In leaves and shoot apical meristems at two developmental stages, 977 and 1121 transcripts were differentially accumulated in S. tetrasperma and S. biocullata, respectively. Interestingly, among the 33 transcription factors showing differential expression in S. tetrasperma but without differential expression in S. biocullata, three were involved potentially in the photoperiod and circadian-clock pathway of flowering time regulation (FAR1 RELATED SEQUENCE 12, FRS12; NUCLEAR FACTOR Y A1, NFYA1; PSEUDO-RESPONSE REGULATOR 5, PRR5), hence provides an important clue in deciphering the flowering diversity mechanisms. Our data serve as a key resource for further dissection of molecular and genetic mechanisms underpinning key biological innovations, here, the perennial monocarpic mass flowering.

Functional characterization of CYP81C16 involved in the tanshinone biosynthetic pathway in Salvia miltiorrhiza

Danshen, the dried roots and rhizomes of Salvia miltiorrhiza Bunge (S. miltiorrhiza), is widely used in the treatment of cardiovascular and cerebrovascular diseases. Tanshinones, the bioactive compounds from Danshen, exhibit a wide spectrum of pharmacological properties, suggesting their potential for future therapeutic applications. Tanshinone biosynthesis is a complex process involving at least six P450 enzymes that have been identified and characterized, most of which belong to the CYP76 and CYP71 families. In this study, CYP81C16, a member of the CYP71 clan, was identified in S. miltiorrhiza. An in vitro assay revealed that it could catalyze the hydroxylation of four para-quinone-type tanshinones, namely neocryptotanshinone, deoxyneocryptotanshinone, and danshenxinkuns A and B. SmCYP81C16 emerged as a potential broad-spectrum oxidase targeting the C-18 position of para-quinone-type tanshinones with an impressive relative conversion rate exceeding 90%. Kinetic evaluations andin vivo assays underscored its highest affinity towards neocryptotanshinone among the tested substrates. The overexpression of SmCYP81C16 promoted the accumulation of (iso)tanshinone in hairy root lines. The characterization of SmCYP81C16 in this study accentuates its potential as a pivotal tool in the biotechnological production of tanshinones, either through microbial or plant metabolic engineering.

Functions of Representative Terpenoids and Their Biosynthesis Mechanisms in Medicinal Plants

Terpenoids are the broadest and richest group of chemicals obtained from plants. These plant-derived terpenoids have been extensively utilized in various industries, including food and pharmaceuticals. Several specific terpenoids have been identified and isolated from medicinal plants, emphasizing the diversity of biosynthesis and specific functionality of terpenoids. With advances in the technology of sequencing, the genomes of certain important medicinal plants have been assembled. This has improved our knowledge of the biosynthesis and regulatory molecular functions of terpenoids with medicinal functions. In this review, we introduce several notable medicinal plants that produce distinct terpenoids (e.g., Cannabis sativa, Artemisia annua, Salvia miltiorrhiza, Ginkgo biloba, and Taxus media). We summarize the specialized roles of these terpenoids in plant-environment interactions as well as their significance in the pharmaceutical and food industries. Additionally, we highlight recent findings in the fields of molecular regulation mechanisms involved in these distinct terpenoids biosynthesis, and propose future opportunities in terpenoid research, including biology seeding, and genetic engineering in medicinal plants.

Integrating Full-Length and Second-Generation Transcriptomics to Reveal Differentially Expressed Genes Associated with the Development of Corydalis yanhusuo Tuber

Corydalis yanhusuo is a medicinal herb in China that has been widely used to treat various kinds of pain. The tuber is the main organ of C. yanhusuo used for medicinal purposes, but changes in related genes during the development of the tuber have rarely been reported. To identify the differentially expressed genes during tuber development, C. yanhusuo full-length transcriptomic sequencing was performed using single-molecule real-time technology, and tubers at three development stages were selected for comparative transcriptome analysis. A total of 90,496 full-length non-chimeric transcripts were obtained, and 19,341 transcripts were annotated in at least one public database. A total of 9221 differentially expressed genes were identified during the swelling process of C. yanhusuo tuber. A Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes associated with a "starch and sucrose metabolism pathway", "phenylpropanoid biosynthesis pathway", "isoquinoline alkaloid biosynthesis pathway", "zeatin biosynthesis pathway", and "brassinosteroid biosynthesis pathway" were predominantly enriched. In addition, the genes involved in cell wall metabolism were potentially associated with tuber swelling. These processes regulated and were involved in C. yanhusuo tuber development. The results provide a foundation for further research on tuber formation in medicinal plants.

Construction of heat stress regulation networks based on Illumina and SMRT sequencing data in potato

Potato (Solanum tuberosum L.) is one of the most important tuber food crops in the world; however, the cultivated potatoes are susceptible to high temperature, by which potato production is adversely affected. Understanding the coping mechanism of potato to heat stress is essential to secure yield and expand adaptability under environmental conditions with rising temperature. However, the lack of heat-related information has significantly limited the identification and application of core genes. To gain deeper insights into heat tolerance genes, next-generation sequencing and single-molecule real-time sequencing were used to learn the transcriptional response of potato to heat stress and 13,159 differentially expressed genes (DEGs) were identified in this study. All DEGs were grouped into 12 clusters using the K-means clustering algorithm. Gene Ontology enrichment analysis revealed that they were involved in temperature signaling, phytohormone, and protein modification. Among them, there were 950 differentially expressed transcription factors (DETFs). According to the network analysis of DETFs at the sixth hour under heat stress, we found some genes that were previously reported to be associated with photoperiodic tuberization, StCO (CONSTANS), tuber formation, StBEL11 (BEL1-LIKE 11), and earliness in potato, StCDF1 (CYCLING DOF FACTOR 1) responding to temperature. Furthermore, we verified the relative expression levels using quantitative real-time polymerase chain reaction, and the results were consistent with the inferences from transcriptomes. In addition, there were 22,125 alternative splicing events and 2,048 long non-coding RNAs. The database and network established in this study will extend our understanding of potato response to heat stress. It ultimately provided valuable resources for molecular analysis of heat stress response in potato and cultivation of potato varieties with heat tolerance.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing

Drought, as a widespread environmental factor in nature, has become one of the most critical factors restricting the yield of forage grass. Sudangrass (Sorghum sudanense (Piper) Stapf.), as a tall and large grass, has a large biomass and is widely used as forage and biofuel. However, its growth and development are limited by drought stress. To obtain novel insight into the molecular mechanisms underlying the drought response and excavate drought tolerance genes in sudangrass, the first full-length transcriptome database of sudangrass under drought stress at different time points was constructed by combining single-molecule real-time sequencing (SMRT) and next-generation transcriptome sequencing (NGS). A total of 32.3 Gb of raw data was obtained, including 20,199 full-length transcripts with an average length of 1628 bp after assembly and correction. In total, 11,921 and 8559 up- and down-regulated differentially expressed genes were identified between the control group and plants subjected to drought stress. Additionally, 951 transcription factors belonging to 50 families and 358 alternative splicing events were found. A KEGG analysis of 158 core genes exhibiting continuous changes over time revealed that 'galactose metabolism' is a hub pathway and raffinose synthase 2 and β-fructofuranosidase are key genes in the response to drought stress. This study revealed the molecular mechanism underlying drought tolerance in sudangrass. Furthermore, the genes identified in this study provide valuable resources for further research into the response to drought stress.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

Chromosome-level genome assembly of Salvia miltiorrhiza with orange roots uncovers the role of Sm2OGD3 in catalyzing 15,16-dehydrogenation of tanshinones

Salvia miltiorrhiza is well known for its clinical practice in treating heart and cardiovascular diseases. Its roots, used for traditional Chinese medicine materials, are usually brick-red due to accumulation of red pigments, such as tanshinone IIA and tanshinone I. Here we report a S. miltiorrhiza line (shh) with orange roots. Compared with the red roots of normal S. miltiorrhiza plants, the contents of tanshinones with a single bond at C-15,16 were increased, whereas those with a double bond at C-15,16 were significantly decreased in shh. We assembled a high-quality chromosome-level genome of shh. Phylogenomic analysis showed that the relationship between two S. miltiorrhiza lines with red roots was closer than the relationship with shh. It indicates that shh could not be the mutant of an extant S. miltiorrhiza line with red roots. Comparative genomic and transcriptomic analyses showed that a 1.0 kb DNA fragment was deleted in shh Sm2OGD3m. Complementation assay showed that overexpression of intact Sm2OGD3 in shh hairy roots recovered furan D-ring tanshinone accumulation. Consistently, in vitro protein assay showed that Sm2OGD3 catalyzed the conversion of cyptotanshinone, 15,16-dihydrotanshinone I and 1,2,15,16-tetrahydrotanshinone I into tanshinone IIA, tanshinone I and 1,2-dihydrotanshinone I, respectively. Thus, Sm2OGD3 functions as tanshinone 15,16-dehydrogenase and is a key enzyme in tanshinone biosynthesis. The results provide novel insights into the metabolic network of medicinally important tanshinone compounds.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

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 analysis of Zanthoxylum nitidum (Roxb.) DC

Zanthoxylum nitidum (Roxb.) DC. (Z. nitidum) is a type of Chinese Dao-di herb, also called Liangmianzhen, which is widely used to treat arthralgia, rheumatic arthralgia, and stomach pain. However, genomic resources for Z. nitidum are still scarce. This study provides transcriptomic resources for Z. nitidum by applying single-molecule real-time (SMRT) sequencing technology. In total, 456,109 circular consensus sequencing (CCS) reads were generated with a mean length of 2,216 bp from Z. nitidum roots, old stems, young branches, leaves, flowers, and fruits. Of these total reads, 353,932 were full-length nonchimeric (FLNC) reads with an average length of 1,996 bp. A total of 16,163 transcripts with a mean length of 1,171 bp were acquired. Of these transcripts, 14,231 (88%) were successfully annotated using public databases. Across all the 16,163 transcripts, we identified 6,255 long non-coding RNAs (lncRNAs) and 22,780 simple sequence repeats (SSRs). Furthermore, 3,482 transcription factors were identified. Among the SSR loci, 1-3 nucleotide repeats were dominant, occupying 99.36% of the total SSR loci, with mono-, di-, and tri-nucleotide repeats accounting for 61.80%, 19.89%, and 5.02% of the total SSR loci, respectively. A total of 36 out of 100 randomly selected primer pairs were verified to be positive, 20 of which showed polymorphism. These findings enrich the genetic resources available for facilitating future studies and research on relevant topics such as population genetics in Z. nitidum.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Isolation of Salvia miltiorrhiza Kaurene Synthase-like (KSL) Gene Promoter and Its Regulation by Ethephon and Yeast Extract

The presented study describes the regulation of the promoter region of the Salvia miltiorrhiza kaurene synthase-like gene (SmKSL) by ethylene and yeast extract. The isolated fragment is 897 bp and is composed of a promoter (763 bp), 5'UTR (109 bp), and a short CDS (25 bp). The initial in silico analysis revealed the presence of numerous putative cis-active sites for trans-factors responding to different stress conditions. However, this study examines the influence of ethylene and yeast extract on SmKSL gene expression and tanshinone biosynthesis regulation. The results of 72h RT-PCR indicate an antagonistic interaction between ethylene, provided as ethephon (0.05, 0.10, 0.25, and 0.50 mM), and yeast extract (0.5%) on SmKSL gene expression in callus cultures of S. miltiorrhiza. A similar antagonistic effect was observed on total tanshinone concentration for up to 60 days. Ethylene provided as ethephon (0.05, 0.10, 0.25, and 0.50 mM) is a weak inducer of total tanshinone biosynthesis, increasing them only up to the maximum value of 0.67 ± 0.04 mg g^-1 DW (60-day induction with 0.50 mM ethephon). Among the tanshinones elicited by ethephon, cryptotanshinone (52.21%) dominates, followed by dihydrotanshinone (45.00%) and tanshinone IIA (3.79%). In contrast, the 0.5% yeast extract strongly increases the total tanshinone concentration up to a maximum value of 13.30 ± 1.09 mg g^-1 DW, observed after 50 days of induction. Yeast extract and ethylene appear to activate different fragments of the tanshinone biosynthesis route; hence the primary tanshinones induced by yeast extract were cryptotanshinone (81.42%), followed by dihydrotanshinone (17.06%) and tanshinone IIA (1.52%).

Metabolome and transcriptome associated analysis of sesquiterpenoid metabolism in Nardostachys jatamansi

BACKGROUND

Nardostachys jatamansi, an extremely endangered valuable plant of the alpine Himalayas, can synthesize specific sesquiterpenoids with multiple effective therapies and is widely exploited for the preparation of drugs, cosmetics and even religious functions (e.g., well-known spikenard). However, how accumulation trend of the sesquiterpenoids in tissues and the molecular mechanisms underlying the production of the active ingredients are not well understood.

METHODS

The single-molecule real-time (SMRT) and RNA-seq transcriptome sequencing were combined to analyse the roots, rhizomes, leaves, flowers and anthocaulus of N. jatamansi. The phytochemical analysis was performed by gas chromatography‒mass spectrometry (GC‒MS) and ultrahigh-performance liquid chromatography (UPLC).

RESULTS

A high-quality full-length reference transcriptome with 26,503 unigenes was generated for the first time. For volatile components, a total of sixty-five compounds were successfully identified, including fifty sesquiterpenoids. Their accumulation levels in five tissues were significantly varied, and most of the sesquiterpenoids were mainly enriched in roots and rhizomes. In addition, five aromatic compounds were only detected in flowers, which may help the plant attract insects for pollination. For nonvolatile ingredients, nardosinone-type sesquiterpenoids (nardosinone, kanshone C, and isonardosinone) were detected almost exclusively in roots and rhizomes. The candidate genes associated with sesquiterpenoid biosynthesis were identified by transcriptome analysis. Consistently, it was found that most biosynthesis genes were abundantly expressed in the roots and rhizomes according to the functional enrichment and expression patterns results. There was a positive correlation between the expression profile of genes related to the biosynthesis and the accumulation level of sesquiterpenoids in tissues. Gene family function analysis identified 28 NjTPSs and 43 NjCYPs that may be involved in the biosynthesis of the corresponding sesquiterpenoids. Furthermore, gene family functional analysis and gene coexpression network analysis revealed 28 NjTPSs and 43 NjCYPs associated with nardosinone-type sesquiterpenoid biosynthesis.

CONCLUSION

Our research results reveal the framework of sesquiterpenoids accumulation and biosynthesis in plant tissues and provide valuable support for further studies to elucidate the molecular mechanisms of sesquiterpenoid regulation and accumulation in N. jatamansi and will also contribute to the comprehensive utilization of this alpine plant.

From single- to multi-omics: future research trends in medicinal plants

Medicinal plants are the main source of natural metabolites with specialised pharmacological activities and have been widely examined by plant researchers. Numerous omics studies of medicinal plants have been performed to identify molecular markers of species and functional genes controlling key biological traits, as well as to understand biosynthetic pathways of bioactive metabolites and the regulatory mechanisms of environmental responses. Omics technologies have been widely applied to medicinal plants, including as taxonomics, transcriptomics, metabolomics, proteomics, genomics, pangenomics, epigenomics and mutagenomics. However, because of the complex biological regulation network, single omics usually fail to explain the specific biological phenomena. In recent years, reports of integrated multi-omics studies of medicinal plants have increased. Until now, there have few assessments of recent developments and upcoming trends in omics studies of medicinal plants. We highlight recent developments in omics research of medicinal plants, summarise the typical bioinformatics resources available for analysing omics datasets, and discuss related future directions and challenges. This information facilitates further studies of medicinal plants, refinement of current approaches and leads to new ideas.

Spatiotemporal and Transcriptional Characterization on Tanshinone Initial Synthesis in Salvia miltiorrhiza Roots

Tanshinones are the bioactive constituents of Danshen (Salvia miltiorrhiza Bunge), which is used in Traditional Chinese Medicine to treat cardiovascular and other diseases, and they synthesize and accumulate in the root periderm of S. miltiorrhiza. However, there is no relevant report on the initial stage of tanshinone synthesis, as well as the root structure and gene expression characteristics. The present study aims to provide new insights into how these bioactive principles begin to synthesize by characterizing possible differences in their biosynthesis and accumulation during early root development from both spatial and temporal aspects. The morphological characteristics and the content of tanshinones in roots of S. miltiorrhiza were investigated in detail by monitoring the seedlings within 65 days after germination (DAGs). The ONT transcriptome sequencing was applied to investigate gene expression patterns. The periderm of the S. miltiorrhiza storage taproot initially synthesized tanshinone on about 30 DAGs. Three critical stages of tanshinone synthesis were preliminarily determined: preparation, the initial synthesis, and the continuous rapid synthesis. The difference of taproots in the first two stages was the smallest, and the differentially expressed genes (DEGs) were mainly enriched in terpene synthesis. Most genes involved in tanshinone synthesis were up regulated during the gradual formation of the red taproot. Plant hormone signal transduction and ABC transport pathways were widely involved in S. miltiorrhiza taproot development. Five candidate genes that may participate in or regulate tanshinone synthesis were screened according to the co-expression pattern. Moreover, photosynthetic ferredoxin (FD), cytochrome P450 reductase (CPR), and CCAAT binding transcription factor (CBF) were predicted to interact with the known downstream essential enzyme genes directly. The above results provide a necessary basis for analyzing the initial synthesis and regulation mechanism of Tanshinones.

De novo full length transcriptome analysis and gene expression profiling to identify genes involved in phenylethanol glycosides biosynthesis in Cistanche tubulosa

Background

The dried stem of Cistanche, is a famous Chinese traditional medicine. The main active pharmacodynamic components are phenylethanol glycosides (PhGs). Cistanche tubulosa produces higher level of PhGs in its stems than that of Cistanche deserticola. However, the key genes in the PhGs biosynthesis pathway is not clear in C. tubulosa.

Results

In this study, we performed the full-length transcriptome sequencing and gene expression profiling of C. tubulosa using PacBio combined with BGISEQ-500 RNA-seq technology. Totally, 237,772 unique transcripts were obtained, ranging from 199 bp to 31,857 bp. Among the unique transcripts, 188,135 (79.12%) transcripts were annotated. Interestingly, 1080 transcripts were annotated as 22 enzymes related to PhGs biosynthesis. We measured the content of echinacoside, acteoside and total PhGs at two development stages, and found that the content of PhGs was 46.74% of dry matter in young fleshy stem (YS1) and then decreased to 31.22% at the harvest stage (HS2). To compare with YS1, 13,631 genes were up-regulated, and 15,521 genes were down regulated in HS2. Many differentially expressed genes (DEGs) were identified to be involved in phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, and tyrosine metabolism pathway.

Conclusions

This is the first report of transcriptome study of C. tubulosa which provided the foundation for understanding of PhGs biosynthesis. Based on these results, we proposed a potential model for PhGs biosynthesis in C. tubulosa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08921-x.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Full-length transcriptome-referenced analysis reveals crucial roles of hormone and wounding during induction of aerial bulbils in lily

Aerial bulbils are important vegetative reproductive organs in Lilium. They are often perpetually dormant in most Lilium species, and little is known about the induction of these vegetative structures. The world-famous Oriental hybrid lily cultivar 'Sorbonne', which blooms naturally devoid of aerial bulbils, is known for its lovely appearance and sweet fragrance. We found that decapitation stimulated the outgrowth of aerial bulbils at lower stems (LSs) and then application of low and high concentrations of IAA promoted aerial bulbils emergence around the wound at upper stems (USs) of 'Sorbonne'. However, the genetic basis of aerial bulbil induction is still unclear. Herein, 'Sorbonne' transcriptome has been sequenced for the first time using the combination of third-generation long-read and next-generation short-read technology. A total of 46,557 high-quality non-redundant full-length transcripts were generated. Transcriptomic profiling was performed on seven tissues and stems with treatments of decapitation and application of low and high concentrations of IAA, respectively. Functional annotation of 1918 DEGs within stem samples of different treatments showed that hormone signaling, sugar metabolism and wound-induced genes were crucial to bulbils outgrowth. The expression pattern of auxin-, shoot branching hormone-, plant defense hormone- and wound-inducing-related genes indicated their crucial roles in bulbil induction. Then we established five hormone- and wounding-regulated co-expression modules and identified some candidate transcriptional factors, such as MYB, bZIP, and bHLH, that may function in inducing bulbils. High connectivity was observed among hormone signaling genes, wound-induced genes, and some transcriptional factors, suggesting wound- and hormone-invoked signals exhibit extensive cross-talk and regulate bulbil initiation-associated genes via multilayered regulatory cascades. We propose that the induction of aerial bulbils at LSs after decapitation can be explained as the release of apical dominance. In contrast, the induction of aerial bulbils at the cut surface of USs after IAA application occurs via a process similar to callus formation. This study provides abundant candidate genes that will deepen our understanding of the regulation of bulbil outgrowth, paving the way for further molecular breeding of lily.

Transcriptomic Insight into Viviparous Growth in Water Lily

Water lily is an important ornamental flower plant which is capable of viviparous plantlet development. But no study has been reported on the molecular basis of viviparity in water lily. Hence, we performed a comparative transcriptome study between viviparous water lily Nymphaea micrantha and a nonviviparous species Nymphaea colorata at four developmental stages. The higher expression of highly conserved AUX/IAA, ARF, GH3, and SAUR gene families in N. micrantha compared to N. colorata is predicted to have a major impact on the development and evolution of viviparity in water lily. Likewise, differential regulation of hormone signaling, brassinosteroid, photosynthesis, and energy-related pathways in the two species provide clues of their involvement in viviparity phenomenon. This study revealed the complex mechanism of viviparity trait in water lily. The transcriptomic signatures identified are important basis for future breeding and research of viviparity in water lily and other plant species.

The Expression Profiles of the Salvia miltiorrhiza 3-Hydroxy-3-methylglutaryl-coenzyme A Reductase 4 Gene and Its Influence on the Biosynthesis of Tanshinones

Salvia miltiorrhiza is a medicinal plant that synthesises biologically-active tanshinones with numerous therapeutic properties. An important rate-limiting enzyme in the biosynthesis of their precursors is 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR). This study presents the organ-specific expression profile of the S. miltiorrhiza HMGR4 gene and its sensitivity to potential regulators, viz. gibberellic acid (GA₃), indole-3-acetic acid (IAA) and salicylic acid (SA). In addition, it demonstrates the importance of the HMGR4 gene, the hormone used, the plant organ, and the culture environment for the biosynthesis of tanshinones. HMGR4 overexpression was found to significantly boost the accumulation of dihydrotanshinone I (DHTI), cryptotanshinone (CT), tanshinone I (TI) and tanshinone IIA (TIIA) in roots by 0.44 to 5.39 mg/g dry weight (DW), as well as TIIA in stems and leaves. S. miltiorrhiza roots cultivated in soil demonstrated higher concentrations of the examined metabolites than those grown in vitro. GA₃ caused a considerable increase in the quantity of CT (by 794.2 µg/g DW) and TIIA (by 88.1 µg/g DW) in roots. In turn, IAA significantly inhibited the biosynthesis of the studied tanshinones in root material.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

A Global Analysis of Alternative Splicing of Dichocarpum Medicinal Plants, Ranunculales

Background: The multiple isoforms are often generated from a single gene via Alternative Splicing (AS) in plants, and the functional diversity of the plant genome is significantly increased. Despite well-studied gene functions, the specific functions of isoforms are little known, therefore, the accurate prediction of isoform functions is exceedingly wanted. Methods: Here we perform the first global analysis of AS of Dichocarpum, a medicinal genus of Ranunculales, by utilizing full-length transcriptome datasets of five Chinese endemic Dichocarpum taxa. Multiple software were used to identify AS events, the gene function was annotated based on seven databases, and the protein-coding sequence of each AS isoform was translated into an amino acid sequence. The self-developed software DIFFUSE was used to predict the functions of AS isoforms. Results: Among 8,485 genes with AS events, the genes with two isoforms were the most (6,038), followed by those with three isoforms and four isoforms. Retained intron (RI, 551) was predominant among 1,037 AS events, and alternative 3' splice sites and alternative 5' splice sites were second. The software DIFFUSE was effective in predicting functions of Dichocarpum isoforms, which have not been unearthed. When compared with the sequence alignment-based database annotations, DIFFUSE performed better in differentiating isoform functions. The DIFFUSE predictions on the terms GO:0003677 (DNA binding) and GO: 0010333 (terpene synthase activity) agreed with the biological features of transcript isoforms. Conclusion: Numerous AS events were for the first time identified from full-length transcriptome datasets of five Dichocarpum taxa, and functions of AS isoforms were successfully predicted by the self-developed software DIFFUSE. The global analysis of Dichocarpum AS events and predicting isoform functions can help understand the metabolic regulations of medicinal taxa and their pharmaceutical explorations.

Integrated SMRT and Illumina Sequencing Provide New Insights into Crocin Biosynthesis of Gardenia jasminoides

Crocins are valuable bioactive components of gardenia fruit, and their biosynthesis and accumulation have attracted widespread interest. Studies have investigated the biosynthesis and accumulation of crocin based on Illumina sequencing, but there is a lack of reports based on full-length transcriptome sequencing. Utilising SMRT sequencing and high-performance liquid chromatography (HPLC), we explored crocin biosynthesis and accumulation in the fruit of Gardenia jasminoides. HPLC analysis showed that crocins specifically exist in fruit and that the content of crocins increases gradually during fruit development. SMRT sequencing generated 46,715 high-quality full-length isoforms, including 5230 novel isoforms that are not present in the G. jasminoides genome. Furthermore, a total of 46 genes and 91 lncRNAs were involved in the biosynthesis and accumulation of crocin. The qRT-PCR indicated that genes involved in crocin biosynthesis reached a peak in the NOV stage. These findings contributed to our understanding of crocin biosynthesis and accumulation.

SmGDB: genome database of Salvia miltiorrhiza, an important TCM Plant

BACKGROUND

Salvia miltiorrhiza is an important traditional Chinese medicinal (TCM) plant and a model plant in the genetic study of TCM. A series of omics related to Danshen have been published. Integrating, managing, storing, and sharing data has become an urgent problem to be solved in S. miltiorrhiza genetic studies.

OBJECTIVES

The genome database is the link for the exchange, acquisition, and use of different omics data between data producers and users, maximizing value and utilization of data.

METHODS

The genome database included DSS3 genome and five RNA-Seq data. The back-end performs data search and retrieval through the LAMP (Linux, Apache, MySQL, PHP) framework.

RESULTS

Here, we present SmGDB (S. miltiorrhiza genome database; http://8.140.162.85/ ), which houses the latest version of genome sequence and annotation data for S. miltiorrhiza, combining three unpublished RNA-Seq data from our group and two released RNA-Seq data. We also identified a novel gene cluster including seven CYP71D genes involved in the tanshinone synthesis pathway based on genome sequences and expression data. Besides, SmGDB provides user-friendly web interfaces for querying and browsing gene annotation, structure, location, and expression profiles for concerned genes. Popular bioinformatics tools such as 'BLAST', 'Search', 'Heatmap', 'JBrowse', etc., were also provided in SmGDB.

CONCLUSIONS

SmGDB will provide utility for characterizing the structure of the S. miltiorrhiza genome and better understanding gene functions and biological processes underlying complex secondary metabolism in Danshen.

Integration of full-length transcriptomes and anthocyanin metabolite analysis for understanding fruit coloration mechanism in Schisandra chinensis

Coloration directly affects the commercial value of Schisandra chinensis fruits. The composition and content of anthocyanin determine the S. chinensis fruit coloration. However, the molecular mechanism of anthocyanin biosynthesis and regulation in this fruit remains unknown. In this study, we performed integrative full-length transcriptomics and targeted metabolomics analyses in S. chinensis fruits at four different developmental stages to elucidate the coloration mechanism. Cyanidin3-O-xyl-rutinoside is the key anthocyanin, which is responsible for the reddening of S. chinensis fruits, and its accumulation gradually accelerated from the 80th day after fluorescence. Overall, 122,289 unigenes with an average length of 2592 bp and an N50 of 4232 bp were obtained through single-molecule real-time sequencing; a total of 16,456 differentially expressed genes were identified. Moreover, 10 full-length structural genes related to anthocyanin biosynthesis were found to be significantly differentially expressed with fruit ripening. Moreover, 10 glycosyltransferases (GTs) that may possess the activities of anthocyanidin 3-O-glucosyltransferase, anthocyanidin 3-O-glucoside rhamnosyltransferase, and xylosyltransferases, which are involved in the final three steps for cyanidin3-O-xyl-rutinoside synthesis, were identified through phylogenetic analysis. Based on these findings, we constructed the complete anthocyanin biosynthetic pathway in S. chinensis fruits; five ScMYBs, three ScbHLHs, and two ScWD40s potentially involved in regulating anthocyanin biosynthesis in S. chinensis fruits were also selected. Our study provides the foundation for further research on the molecular mechanism of anthocyanin biosynthesis and regulation for improving the quality of S. chinensis fruits. The results of full-length transcriptomes would provide researchers with novel insights into the molecular cloning of enzymes and their activity.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01179-3.

Tanshinones: Leading the way into Lamiaceae labdane-related diterpenoid biosynthesis

Tanshinones are the bioactive diterpenoid constituents of the traditional Chinese medicinal herb Danshen (Salvia miltiorrhiza), and are examples of the phenolic abietanes widely found within the Lamiaceae plant family. Due to the significant interest in these labdane-related diterpenoid natural products, their biosynthesis has been intensively investigated. In addition to providing the basis for metabolic engineering efforts, this work further yielded pioneering insights into labdane-related diterpenoid biosynthesis in the Lamiaceae more broadly. This includes stereochemical foreshadowing of aromatization, with novel protein domain loss in the relevant diterpene synthase, as well as broader phylogenetic conservation of the relevant enzymes. Beyond such summary of more widespread metabolism, formation of the furan ring that characterizes the tanshinones also has been recently elucidated. Nevertheless, the biocatalysts for the pair of demethylations remain unknown, and the intriguing potential connection of these reactions to the further aromatization observed in the tanshinones are speculated upon here.

Transcriptome Analysis of Stephania tetrandra and Characterization of Norcoclaurine-6-O-Methyltransferase Involved in Benzylisoquinoline Alkaloid Biosynthesis

Stephania tetrandra (S. Moore) is a source of traditional Chinese medicine that is widely used to treat rheumatism, rheumatoid arthritis, edema, and hypertension. Benzylisoquinoline alkaloids (BIAs) are the main bioactive compounds. However, the current understanding of the biosynthesis of BIAs in S. tetrandra is poor. Metabolite and transcriptomic analyses of the stem, leaf, xylem, and epidermis of S. tetrandra were performed to identify candidate genes associated with BIAs biosynthesis. According to the metabolite analysis, the majority of the BIAs accumulated in the root, especially in the epidermis. Transcriptome sequencing revealed a total of 113,338 unigenes that were generated by de novo assembly. Among them, 79,638 unigenes were successfully annotated, and 42 candidate structural genes associated with 15 steps of BIA biosynthesis identified. Additionally, a new (S)-norcoclaurine-6-O-methyltransferase (6OMT) gene was identified in S. tetrandra, named St6OMT2. Recombinant St6OMT2 catalyzed (S)-norcoclaurine methylation to form (S)-coclaurine in vitro. Maximum activity of St6OMT2 was determined at 30°C and pH 6.0 in NaAc-HAc buffer. Its half-life at 50°C was 22 min with the K_m and k_cat of 28.2 μM and 1.5 s^-1, respectively. Our results provide crucial transcriptome information for S. tetrandra, shedding light on the understanding of BIAs biosynthesis and further gene functional characterization.

Phytotoxic effect and molecular mechanism induced by graphene towards alfalfa (Medicago sativa L.) by integrating transcriptomic and metabolomics analysis

Although the widespread use of nanoparticles has been reported in various fields, the toxic mechanisms of molecular regulation involved in the alfalfa treated by nanomaterials is still in the preliminary research stage. In this study, Bara 310 SC (Bara, tolerant genotype) and Gold Empress (Gold, susceptible genotype) were used to investigate how the leaves of alfalfa interpret the physiological responses to graphene stress based on metabolome and transcriptome characterizations. Herein, graphene at different concentrations (0, 1% and 2%, w/w) were selected as the analytes. Physiological results showed antioxidant defence system and photosynthesis was significantly disturbed under high environmental concentration of graphene. With Ultra high performance liquid chromatography electrospray tandem mass spectrometry (UPLC-ESI-MS/MS), 406 metabolites were detected and 62/13 and 110/58 metabolites significantly changed in the leaves of Gold/Bara under the 1% and 2%-graphene treatments (w/w), respectively. The most important metabolites which were accumulated under graphene stress includes amino acids, flavonoids, organic acids and sugars. Transcriptomic analysis reveals 1125 of core graphene-responsive genes in alfalfa that was robustly differently expressed in both genotypes. And differential expression genes (DEGs) potentially related to photosynthetic enzymes, antioxidant enzymes, amino acids metabolism, and sucrose and starch metabolic which finding was supported by the metabolome study. Gold was more disturbed by graphene stress at both transcriptional and metabolic levels, since more stress-responsive genes/metabolites were identified in Gold. A comprehensive analysis of transcriptomic and metabolomic data highlights the important role of amino acid metabolism and nicotinate and nicotinamide metabolism pathways for graphene tolerance in alfalfa. Our study provide necessary information for better understanding the phytotoxicity molecular mechanism underlying nanomaterials tolerance of plant.

Endophytic fungus Cladosporium tenuissimum DF11, an efficient inducer of tanshinone biosynthesis in Salvia miltiorrhiza roots

Salvia miltiorrhiza is a traditional medicinal plant mainly used for cardiovascular and cerebrovascular disease treatment. Tanshinones are the main bioactive constituents of S. miltiorrhiza, which mainly accumulate around its root periderm tissue. Endophytic fungi are important bioelicitors or probiotics that can promote the accumulation of secondary metabolites and sustainable cultivation of medicinal plants. Among them, endophytic Cladosporium spp., possessing a variety of biotransformation and metabolic abilities, is an ideal elicitor source. Here, we used a gnotobiotic system to investigate the effects of the endophytic fungus Cladosporium tenuissimum DF11 on tanshinone biosynthesis in S. miltiorrhiza roots. The results showed that C. tenuissimum DF11 mainly colonizes the intercellular space of the root tissues and promotes tanshinone biosynthesis and accumulation in S. miltiorrhiza roots by upregulating the expression of the genes encoding for key enzymes HMGR, DXS, DXR, GGPPS, CPS, KSL and CYP76AH1 of the tanshinone biosynthesis pathway. The expression levels of almost all genes encoding for key enzymes reached the response peak in the first or second week after DF11 colonization. Taken together, the endophytic fungus C. tenuissimum DF11 could promote secondary metabolite accumulation in S. miltiorrhiza roots. These results indicate that DF11 will be a potential biofertilizer fungus to regulate and stabilize the quality of cultivated S. miltiorrhiza medicinal materials.

Application of High-Throughput Sequencing on the Chinese Herbal Medicine for the Data-Mining of the Bioactive Compounds

The Chinese Herbal Medicine (CHM) has been used worldwide in clinic to treat the vast majority of human diseases, and the healing effect is remarkable. However, the functional components and the corresponding pharmacological mechanism of the herbs are unclear. As one of the main means, the high-throughput sequencing (HTS) technologies have been employed to discover and parse the active ingredients of CHM. Moreover, a tremendous amount of effort is made to uncover the pharmacodynamic genes associated with the synthesis of active substances. Here, based on the genome-assembly and the downstream bioinformatics analysis, we present a comprehensive summary of the application of HTS on CHM for the synthesis pathways of active ingredients from two aspects: active ingredient properties and disease classification, which are important for pharmacological, herb molecular breeding, and synthetic biology studies.

Transcriptome analysis of Kentucky bluegrass subject to drought and ethephon treatment

Kentucky bluegrass (Poa pratensis L.) is an excellent cool-season turfgrass utilized widely in Northern China. However, turf quality of Kentucky bluegrass declines significantly due to drought. Ethephon seeds-soaking treatment has been proved to effectively improve the drought tolerance of Kentucky bluegrass seedlings. In order to investigate the effect of ethephon leaf-spraying method on drought tolerance of Kentucky bluegrass and understand the underlying mechanism, Kentucky bluegrass plants sprayed with and without ethephon are subjected to either drought or well watered treatments. The relative water content and malondialdehyde conent were measured. Meanwhile, samples were sequenced through Illumina. Results showed that ethephon could improve the drought tolerance of Kentucky bluegrass by elevating relative water content and decreasing malondialdehyde content under drought. Transcriptome analysis showed that 58.43% transcripts (254,331 out of 435,250) were detected as unigenes. A total of 9.69% (24,643 out of 254,331) unigenes were identified as differentially expressed genes in one or more of the pairwise comparisons. Differentially expressed genes due to drought stress with or without ethephon pre-treatment showed that ethephon application affected genes associated with plant hormone, signal transduction pathway and plant defense, protein degradation and stabilization, transportation and osmosis, antioxidant system and the glyoxalase pathway, cell wall and cuticular wax, fatty acid unsaturation and photosynthesis. This study provides a theoretical basis for revealing the mechanism for how ethephon regulates drought response and improves drought tolerance of Kentucky bluegrass.

Dissection of full-length transcriptome and metabolome of Dichocarpum (Ranunculaceae): implications in evolution of specialized metabolism of Ranunculales medicinal plants

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

Characterization of NAC family genes in Salvia miltiorrhiza and NAC2 potentially involved in the biosynthesis of tanshinones

The NAC (NAM, ATAF, and CUC) family members are specific transcription factors in plants. The large family is involved in many plant growth and developmental processes, as well as in abiotic/biotic stress responses. It has been well studied in the genomes of various plants, including Arabidopsis thaliana, tomato, and quinoa. However, identification and functional studies of NAC family members in medicinal Salvia miltiorrhiza are limited. Here, we systematically identified 84 NAC genes and named them according to their gene IDs in the recently sequenced genome. The phylogeny of NAC family protein sequences was analyzed using bioinformatics methods, which divided them into nine subfamilies. Then, their chromosomal locations, gene structures and conserved domains were analyzed comprehensively. To further investigate the regulatory functions of NACs in S. miltiorrhiza, we analyzed the response of 10 selected NAC genes to methyl jasmonate and used NAC2 for transgenic experiments. The overexpression of Sm-NAC2 decreased the tanshinone I and IIA contents by 56% and 62%, respectively. However, Sm-NAC2-RNAi promoted the accumulation of four tanshinones, tanshinone I, tanshinone IIA, cryptotanshinone, and dihydrotanshinone I, which increased 3.68-, 4.1-, 3.13- and 5.9- fold, respectively, compared with wild type. In the tanshinone biosynthetic pathways, the overexpression of Sm-NAC2 down-regulated CYP76AH1, and the silencing of Sm-NAC2 up-regulated the expression levels of HMGR1, DXS2, KSL2, and CYP76AH1. This study provides information on the evolution of Sm-NAC genes and their possible functions, and it lays a foundation for further research into the NAC family-associated regulation of tanshinone biosynthesis.

Functional identification of the terpene synthase family involved in diterpenoid alkaloids biosynthesis in Aconitum carmichaelii

Aconitum carmichaelii is a high-value medicinal herb widely used across China, Japan, and other Asian countries. Aconitine-type diterpene alkaloids (DAs) are the characteristic compounds in Aconitum. Although six transcriptomes, based on short-read next generation sequencing technology, have been reported from the Aconitum species, the terpene synthase (TPS) corresponding to DAs biosynthesis remains unidentified. We apply a combination of Pacbio isoform sequencing and RNA sequencing to provide a comprehensive view of the A. carmichaelii transcriptome. Nineteen TPSs and five alternative splicing isoforms belonging to TPS-b, TPS-c, and TPS-e/f subfamilies were identified. In vitro enzyme reaction analysis functional identified two sesqui-TPSs and twelve diTPSs. Seven of the TPS-c subfamily genes reacted with GGPP to produce the intermediate ent-copalyl diphosphate. Five AcKSLs separately reacted with ent-CPP to produce ent-kaurene, ent-atiserene, and ent-13-epi-sandaracopimaradie: a new diterpene found in Aconitum. AcTPSs gene expression in conjunction DAs content analysis in different tissues validated that ent-CPP is the sole precursor to all DAs biosynthesis, with AcKSL1, AcKSL2s and AcKSL3-1 responsible for C20 atisine and napelline type DAs biosynthesis, respectively. These data clarified the molecular basis for the C20-DAs biosynthetic pathway in A. carmichaelii and pave the way for further exploration of C19-DAs biosynthesis in the Aconitum species.

Single-Molecule Real-Time and Illumina Sequencing to Analyze Transcriptional Regulation of Flavonoid Synthesis in Blueberry

Blueberries (Vaccinium corymbosum) contain large amounts of flavonoids, which play important roles in the plant's ability to resist stress and can also have beneficial effects on human health when the fruits are eaten. However, the molecular mechanisms that regulate flavonoid synthesis in blueberries are still unclear. In this study, we combined two different transcriptome sequencing platforms, single-molecule real-time (SMRT) and Illumina sequencing, to elucidate the flavonoid synthetic pathways in blueberries. We analyzed transcript quantity, length, and the number of annotated genes. We mined genes associated with flavonoid synthesis (such as anthocyanins, flavonols, and proanthocyanidins) and employed fluorescence quantitative PCR to analyze the expression of these genes and their correlation with flavonoid synthesis. We discovered one R2R3 MYB transcription factor from the sequencing library, VcMYB1, that can positively regulate anthocyanin synthesis in blueberries. VcMYB1 is mainly expressed in colored (mature) fruits. Experiments showed that overexpression and transient expression of VcMYB1 promoted anthocyanin synthesis in Arabidopsis, tobacco (Nicotiana benthamiana) plants and green blueberry fruits. Yeast one-hybrid (Y1H) assay, electrophoretic mobility shift assay, and transient expression experiments showed that VcMYB1 binds to the MYB binding site on the promoter of the structural gene for anthocyanin synthesis, VcMYB1 to positively regulate the transcription of VcDFR, thereby promoting anthocyanin synthesis. We also performed an in-depth investigation of transcriptional regulation of anthocyanin synthesis. This study provides background information and data for studying the synthetic pathways of flavonoids and other secondary metabolites in blueberries.

Genome-Wide Analysis of Light-Regulated Alternative Splicing in Artemisia annua L

Artemisinin is currently the most effective ingredient in the treatment of malaria, which is thus of great significance to study the genetic regulation of Artemisia annua. Alternative splicing (AS) is a regulatory process that increases the complexity of transcriptome and proteome. The most common mechanism of alternative splicing (AS) in plant is intron retention (IR). However, little is known about whether the IR isoforms produced by light play roles in regulating biosynthetic pathways. In this work we would explore how the level of AS in A. annua responds to light regulation. We obtained a new dataset of AS by analyzing full-length transcripts using both Illumina- and single molecule real-time (SMRT)-based RNA-seq as well as analyzing AS on various tissues. A total of 5,854 IR isoforms were identified, with IR accounting for the highest proportion (48.48%), affirming that IR is the most common mechanism of AS. We found that the number of up-regulated IR isoforms (1534/1378, blue and red light, respectively) was more than twice that of down-regulated (636/682) after treatment of blue or red light. In the artemisinin biosynthetic pathway, 10 genes produced 16 differentially expressed IR isoforms. This work demonstrated that the differential expression of IR isoforms induced by light has the potential to regulate sesquiterpenoid biosynthesis. This study also provides high accuracy full-length transcripts, which can be a valuable genetic resource for further research of A. annua, including areas of development, breeding, and biosynthesis of active compounds.