UPLC/Q-TOF MS-based metabolomics and qRT-PCR in enzyme gene screening with key role in triterpenoid saponin biosynthesis of Polygala tenuifolia

Molecular cloning, functional characterization and expression of the β-amyrin synthase gene involved in saikosaponin biosynthesis in Bupleurum chinense DC

Bupleurum chinense DC. is a commonly used plant in traditional Chinese medicine, and saikosaponins(SSs) are the main active oleanane-typetriterpene saponins in B. chinense. β-Amyrin synthase (β-AS) is an important enzyme in oleanane-type triterpenoid saponin synthesis, but its role in saikosaponin synthesis has rarely been studied. Here, the putative β-AS gene BcBAS1(Accession No.ON890382) selected according to metabolomic and transcriptomic analyses was cloned and functionally characterized by heterologous expression in Escherichia coli and Pichia pastoris, and its subcellular localization and expression patterns were examined. The molecular weight of the BcBAS1 recombinant protein was approximately 87 kDa, and this protein could catalyse the production of β-amyrin, the precursor of SSs. Furthermore, BcBAS1 was located in the cytosol, and relative expression in four tissues of the four genotypes was positively correlated with SSa and SSd contents. Our results indicate that BcBAS1 is a β-AS gene and may play an important role in saikosaponin biosynthesis and regulation. This study sheds light on the role of β-AS genes in the synthesis of SSs and provides insights for the metabolic engineering of SSs.

Recent Advances in Separation and Analysis of Saponins in Natural Products

To better control the quality of saponins, ensure their biological activity and clinical therapeutic effect, and expand the development and application of saponins, this paper systematically and comprehensively reviews the separation and analytical methods of saponins in the past decade. Since 2010, the electronic databases of PubMed, Google Scholar, ISI Web of Science, Science Direct, Wiley, Springer, CNKI (National Knowledge Infrastructure, CNKI), Wanfang Med online, and other databases have been searched systematically. As a result, it is found that ionic liquids and high-performance countercurrent chromatography are the most popular extraction and separation techniques for saponins, and the combined chromatography technique is the most widely used method for the analysis of saponins. Liquid chromatography can be used in combination with different detectors to achieve qualitative or quantitative analysis and quality control of saponin compounds in medicinal materials and their preparations. This paper provides the latest valuable insights and references for the analytical methods and continued development and application of saponins.

Cloning, Yeast Expression, and Characterization of a β-Amyrin C-28 Oxidase (CYP716A249) Involved in Triterpenoid Biosynthesis in Polygala tenuifolia

Polygala tenuifolia Willd. is a traditional Chinese herbal medicine that is widely used in treating nervous system disorders. Triterpene saponins in P. tenuifolia (polygala saponins) have excellent biological activity. As a precursor for the synthesis of presenegin, oleanolic acid (OA) plays an important role in the biosynthesis of polygala saponins. However, the mechanism behind the biosynthesis of polygala saponins remains to be elucidated. In this study, we found that CYP716A249 (GenBank: ASB17946) oxidized the C-28 position of β-amyrin to produce OA. Using quantitative real-time PCR, we observed that CYP716A249 had the highest expression in the roots of 2-year-old P. tenuifolia, which provided a basis for the selection of samples for gene cloning. To identify the function of CYP716A249, the strain R-BE-20 was constructed by expressing β-amyrin synthase in yeast. Then, CYP716A249 was co-expressed with β-amyrin synthase to construct the strain R-BPE-20 by using the lithium acetate method. Finally, we detected β-amyrin and OA by ultra-HPLC-Q Exactive hybrid quadrupole-Orbitrap high-resolution accurate mass spectrometry and GC-MS. The results of this study provide insights into the biosynthesis pathway of polygala saponins.

Global transcriptome analysis uncovers the gene co-expression regulation network and key genes involved in grain development of wheat (Triticum aestivum L.)

Wheat grain development is a robust biological process that largely determines grain quality and yield. In this study, we investigated the grain transcriptome of winter wheat cv. Xiaoyan-6 at four developmental stages (5, 10, 15, and 20 days post-anthesis), using high-throughput RNA sequencing (RNA-Seq). We identified 427 grain-specific transcription factors (TFs) and 1653 differentially expressed TFs during grain development as well as a grain co-expression regulation network (GrainNet) of the TFs and their predicted co-expressed genes. Our study identified ten putative key TFs and the predicted regulatory genes of these TFs in wheat grain development of Xiaoyan-6. The analysis was given a firm basis through the study of additional wheat tissues, including root, stem, leaf, flag leaf, grain, spikes (from wheat plants at booting or heading stages) to provide a dataset of 92,478 high-confidence protein-coding genes that were mostly evenly distributed among subgenomes, but unevenly distributed across each of the chromosomes or each of the seven homeologous groups. Within this larger framework of the transcriptomes, we identified 4659 grain-specific genes (SEGs) and 26,500 differentially expressed genes (DEGs) throughout grain development stages tested. The SEGs identified mainly associate with regulation and signaling-related biological processes, while the DEGs mainly involve in cellular component organization or biogenesis and nutrient reservoir activity during grain development of Xiaoyan-6. This study establishes new targets for modifying genes related to grain development and yield, to fine-tune expression in different varieties and environments.

Nitric oxide alleviates cell death through protein S-nitrosylation and transcriptional regulation during the ageing of elm seeds

Seed ageing is a major problem in the conservation of germplasm resources. The involvement of possible signalling molecules during seed deterioration needs to be identified. In this study, we confirmed that nitric oxide (NO), a key signalling molecule in plants, plays a positive role in the resistance of elm seeds to deterioration. To explore which metabolic pathways were affected by NO, an untargeted metabolomic analysis was conducted, and 163 metabolites could respond to both NO and the ageing treatment. The primary altered pathways include glutathione, methionine, and carbohydrate metabolism. The genes involved in glutathione and methionine metabolism were up-regulated by NO at the transcriptional level. Using a biotin switch method, proteins with an NO-dependent post-translational modification were screened during seed deterioration, and 82 putative S-nitrosylated proteins were identified. Eleven of these proteins were involved in carbohydrate metabolism, and the activities of the three enzymes were regulated by NO. In combination, the results of the metabolomic and S-nitrosoproteomic studies demonstrated that NO could activate glycolysis and inhibit the pentose phosphate pathway. In summary, the combination of these results demonstrated that NO could modulate carbohydrate metabolism at the post-translational level and regulate glutathione and methionine metabolism at the transcriptional level. It provides initial insights into the regulatory mechanisms of NO in seed deterioration.

UPLC Quantitative Analysis of Multi-Components by Single Marker and Quality Evaluation of Polygala tenuifolia Wild. Extracts

The quality control of Polygala tenuifolia Wild. is a major challenge in its clinical application. In this paper, a new strategy for the quality evaluation of P. tenuifolia extracts was verified through reverse-phase ultra-performance liquid chromatography (UPLC). The quantitative analysis of multi-components by a single marker (QAMS) was conducted with 3,6′-disinapoyl sucrose as an internal reference substance. Eight components (i.e., sibiricose A5, sibiricose A6, glomeratose A, tenuifoliside A, tenuifoliside B, tenuifoliside C, sibiricaxanthone B, and polygalaxanthone III) were determined based on the relative correction factors. The concentrations of these components were also determined by applying a conventional external standard method. The cosine value confirmed the consistency of the two methods (cosine ratio value >0.999920). Hierarchical cluster analysis, radar plots, and discriminant analysis were performed to classify 23 batches of P. tenuifolia extracts from Shanxi, Hebei, and Shaanxi in China. Results revealed that QAMS combined with radar plots and multivariate data analysis could accurately measure and clearly distinguish the different quality samples of P. tenuifolia. Hence, QAMS is a feasible and promising method for the quality control of P. tenuifolia.

Utilization of UPLC/Q-TOF-MS-Based Metabolomics and AFLP-Based Marker-Assisted Selection to Facilitate/Assist Conventional Breeding of Polygala tenuifolia

As one of the most important traditional Chinese medicine, the quality of Polygala tenuifolia is difficult to control and a new method must be established to facilitate/assist the breeding of P. tenuifolia. In this study, UPLC/Q-TOF-MS-based metabolomics analysis was performed to determine the chemical composition and screen metabolite biomarkers according to agronomic traits. A total of 29 compounds and 18 metabolite biomarkers were found. AFLP-based marker-assisted selection (MAS) was used to identify molecular marker bands and screen characteristic bands associated with specific agronomic traits. 184 bands and 76 characteristic AFLP bands were found. The correlation network between compounds and characteristic AFLP bands was built, so we may directly breed certain P. tenuifolia herbs with special agronomic traits (or characteristic AFLP bands), which exhibit specific pharmacological functions depending on the content of the active compounds. The proposed method of metabolomics coupled with MAS could facilitate/assist the breeding of P. tenuifolia.

Predicting the Potential Distribution of Polygala tenuifolia Willd. under Climate Change in China

Global warming has created opportunities and challenges for the survival and development of species. Determining how climate change may impact multiple ecosystem levels and lead to various species adaptations is necessary for both biodiversity conservation and sustainable biological resource utilization. In this study, we employed Maxent to predict changes in the habitat range and altitude of Polygala tenuifolia Willd. under current and future climate scenarios in China. Four representative concentration pathways (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) were modeled for two time periods (2050 and 2070). The model inputs included 732 presence points and nine sets of environmental variables under the current conditions and the four RCPs in 2050 and 2070. The area under the receiver-operating characteristic (ROC) curve (AUC) was used to evaluate model performance. All of the AUCs were greater than 0.80, thereby placing these models in the "very good" category. Using a jackknife analysis, the precipitation in the warmest quarter, annual mean temperature, and altitude were found to be the top three variables that affect the range of P. tenuifolia. Additionally, we found that the predicted highly suitable habitat was in reasonable agreement with its actual distribution. Furthermore, the highly suitable habitat area was slowly reduced over time.

Transcriptomic Analysis of Multipurpose Timber Yielding Tree Neolamarckia cadamba during Xylogenesis Using RNA-Seq

Neolamarckia cadamba is a fast-growing tropical hardwood tree that is used extensively for plywood and pulp production, light furniture fabrication, building materials, and as a raw material for the preparation of certain indigenous medicines. Lack of genomic resources hampers progress in the molecular breeding and genetic improvement of this multipurpose tree species. In this study, transcriptome profiling of differentiating stems was performed to understand N. cadamba xylogenesis. The N. cadamba transcriptome was sequenced using Illumina paired-end sequencing technology. This generated 42.49 G of raw data that was then de novo assembled into 55,432 UniGenes with a mean length of 803.2bp. Approximately 47.8% of the UniGenes (26,487) were annotated against publically available protein databases, among which 21,699 and 7,754 UniGenes were assigned to Gene Ontology categories (GO) and Clusters of Orthologous Groups (COG), respectively. 5,589 UniGenes could be mapped onto 116 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Among 6,202 UniGenes exhibiting differential expression during xylogenesis, 1,634 showed significantly higher levels of expression in the basal and middle stem segments compared to the apical stem segment. These genes included NAC and MYB transcription factors related to secondary cell wall biosynthesis, genes related to most metabolic steps of lignin biosynthesis, and CesA genes involved in cellulose biosynthesis. This study lays the foundation for further screening of key genes associated with xylogenesis in N. cadamba as well as enhancing our understanding of the mechanism of xylogenesis in fast-growing trees.

Analysis of Polygala tenuifolia Transcriptome and Description of Secondary Metabolite Biosynthetic Pathways by Illumina Sequencing

Radix polygalae, the dried roots of Polygala tenuifolia and P. sibirica, is one of the most well-known traditional Chinese medicinal plants. Radix polygalae contains various saponins, xanthones, and oligosaccharide esters and these compounds are responsible for several pharmacological properties. To provide basic breeding information, enhance molecular biological analysis, and determine secondary metabolite biosynthetic pathways of P. tenuifolia, we applied Illumina sequencing technology and de novo assembly. We also applied this technique to gain an overview of P. tenuifolia transcriptome from samples with different years. Using Illumina sequencing, approximately 67.2% of unique sequences were annotated by basic local alignment search tool similarity searches against public sequence databases. We classified the annotated unigenes by using Nr, Nt, GO, COG, and KEGG databases compared with NCBI. We also obtained many candidates CYP450s and UGTs by the analysis of genes in the secondary metabolite biosynthetic pathways, including putative terpenoid backbone and phenylpropanoid biosynthesis pathway. With this transcriptome sequencing, future genetic and genomics studies related to the molecular mechanisms associated with the chemical composition of P. tenuifolia may be improved. Genes involved in the enrichment of secondary metabolite biosynthesis-related pathways could enhance the potential applications of P. tenuifolia in pharmaceutical industries.

Integration of a Decrescent Transcriptome and Metabolomics Dataset of Peucedanum praeruptorum to Investigate the CYP450 and MDR Genes Involved in Coumarins Biosynthesis and Transport

Peucedanum praeruptorum Dunn is well-known traditional Chinese medicine. However, little is known in the biosynthesis and the transport mechanisms of its coumarin compounds at the molecular level. Although transcriptomic sequence is playing an increasingly significant role in gene discovery, it is not sufficient in predicting the specific function of target gene. Furthermore, there is also a huge database to be analyzed. In this study, RNA sequencing assisted transcriptome dataset and high-performance liquid chromatography (HPLC) coupled with electrospray-ionization quadrupole time-of-flight mass spectrometry (Q-TOF MS)-based metabolomics dataset of P. praeruptorum were firstly constructed for gene discovery and compound identification. Subsequently, methyl jasmonate (MeJA)-induced gene expression analysis and metabolomics analysis were conducted to narrow-down the dataset for selecting the candidate genes and the potential marker metabolites. Finally, the genes involved in coumarins biosynthesis and transport were predicted with parallel analysis of transcript and metabolic profiles. As a result, a total of 40,952 unigenes and 19 coumarin compounds were obtained. Based on the results of gene expression and metabolomics analysis, 7 cytochrome-P450 and 8 multidrug resistance transporter unigenes were selected as candidate genes and 8 marker compounds were selected as biomarkers, respectively. The parallel analysis of gene expression and metabolites accumulation indicated that the gene labeled as 23,746, 228, and 30,922 were related to the formation of the coumarin core compounds whereas 36,276 and 9533 participated in the prenylation, hydroxylation, cyclization or structural modification. Similarly, 1462, 20,815, and 15,318 participated in the transport of coumarin core compounds while 124,029 and 324,293 participated in the transport of the modified compounds. This finding suggested that integration of a decrescent transcriptome and metabolomics dataset could largely narrow down the number of gene to be investigated and significantly improve the efficiency of functional gene predication. In addition, the large amount of transcriptomic data produced from P. praeruptorum and the genes discovered in this study would provide useful information in investigating the biosynthesis and transport mechanism of coumarins.