Genomic characterization of bZIP transcription factors related to andrographolide biosynthesis in Andrographis paniculata.

The basic leucine zipper (bZIP) transcription factor family plays an important role in various biological processes in plants. Andrographis paniculata (Burm.f) Nees, belonging to the family Acanthaceae, has been widely used as an important traditional herb with a wide range of pharmacological activities, such as antivenom, antiretroviral, anticancer and so on. However, there was no comprehensive analysis of bZIP gene family in the Andrographis paniculata been reported. In this study, we identified 62 bZIPs in Andrographis paniculata and grouped them into 12 subfamilies through the phylogenetic tree analysis. The bZIPs in the same groups have similar motif composition, exon-intron structure and domain distribution. In addition, the RNA-seq data gave a reference for selecting candidate bZIPs to make further function verification. Lastly, qRT-PCR analyses revealed seven ApbZIPs (ApbZIP4, ApbZIP19, ApbZIP30, ApbZIP42, ApbZIP50, ApbZIP52, ApbZIP62) were the most highly expressed in leaf and significantly up-regulated with MeJA and ABA treatment which may be involved in biosynthesis regulation of andrographolide. These data pave the way for further revealing the function of the bZIPs in Andrographis paniculata.

Transcriptome Analysis and Identification of the Cholesterol Side Chain Cleavage Enzyme BbgCYP11A1 From Bufo bufo gargarizans.

Bufo bufo gargarizans Cantor are precious medicinal animals in traditional Chinese medicine (TCM). Bufadienolides as the major pharmacological components are generated from the venomous glands of B. bufo gargarizans. Bufadienolides are one type of cardiac aglycone with a six-member lactone ring and have properties of antitumor, cardiotonic, tonsillitis, and anti-inflammatory. The biosynthesis of bufadienolides is complex and unclear. This study explored the transcriptome of three different tissues (skin glands, venom glands, and muscles) of B. bufo gargarizans by high-throughput sequencing. According to the gene tissue-specific expression profile, 389 candidate genes were predicted possibly participating in the bufadienolides biosynthesis pathway. Then, BbgCYP11A1 was identified as a cholesterol side chain cleaving the enzyme in engineering yeast producing cholesterol. Furthermore, the catalytic activity of BbgCYP11A1 was studied with various redox partners. Interestingly, a plant NADPH-cytochrome P450 reductase (CPR) from Anemarrhena asphodeloides showed notably higher production than BbgAdx-2A-BbgAdR from B. bufo gargarizans. These results will provide certainly molecular research to reveal the bufadienolides biosynthesis pathway in B. bufo gargarizans.

Two cardenolide glycosides from the seed fairs of Asclepias curassavica and their cytotoxic activities.

Two cardenolide glycosides, corotoxigenin 3-O-[ß-D-glucopyranosyl-(1→4)-6-deoxy-ß-D-glucopyranoside] (1) and coroglaucigenin 3-O-[ß-D-glucopyranosyl-(1→4)-6-deoxy-ß-D-glucopyranoside] (2), were isolated from the seed fairs of Asclepias curassavica. The structures of 1-2 were determined based on the combination of the analysis of their MS, NMR spectroscopic data and acid hydrolysis. The inhibitory effects of compounds 1 and 2 on human colorectal carcinoma cells (HCT116), non-small cell lung carcinoma cells (A549) and hepatic cancer cells (SMMC-7721) were evaluated. The results showed that both compounds 1 and 2 significantly inhibited the viability, proliferation, and migration of A549, HCT116 and SMMC-7721 cells, suggesting that compounds 1 and 2 can be applied in the treatment of lung, colon and liver cancers in clinical practice. This study may not only provide a scientific basis for clarifying the active ingredients in A. curassavica, but also help to understand its antitumor activity, which can promote the application of A. curassavica in clinical treatment of various cancers.

Transcriptomic investigation of the biochemical function of 7-dehydrocholesterol reductase 1 from the traditional Chinese medicinal plant Anemarrhena asphodeloides Bunge.

Anemarrhena asphodeloides Bunge (Liliaceae) is an important Traditional Chinese Medicine herb, which contains up to 6 % total steroidal saponins (timosaponins) and has multiple pharmacological properties. However, the timosaponin biosynthetic pathway has not been extensively investigated. Here we conducted de novo transcriptome sequencing and analysis of A. asphodeloides Bunge and screened for candidate genes involved in the timosaponin biosynthetic pathway. Targeted metabolite analysis showed that timosaponins primarily accumulated in rhizomes, while phytosterols (including cholesterol) were distributed throughout various organs. Most of the identified candidate genes of the timosaponin biosynthetic pathway were also highly expressed in the rhizome, consistent with the results of metabolic analysis. Based on the transcriptome results, two candidate 7-dehydrocholesterol reductase genes were cloned and heterologously expressed in the yeast Saccharomyces cerevisiae. The purified and identified products supported that Aa7DR1 possessed Δ7-reduction activity in yeast and therefore may be involved in the timosaponins biosynthetic pathway in A. asphodeloides Bunge. Phylogenetic analysis showed Aa7DR1 belongs to monocotyledonous Δ7 reductase of phytosterol biosynthesis. These data expand our understanding of timosaponin biosynthesis.

Functional Characterization of a Novel Glycosyltransferase (UGT73CD1) from Iris tectorum Maxim. for the Substrate promiscuity.

Glycosylflavonoids are a class of natural products with multiple pharmacological activities and a lot of glycosyltransferases from various plant species have been reported that they were involved in the biosynthesis of these phytochemicals. However, no corresponding glycosyltransferase has been identified from the famous horticultural and medicinal plant Iris tectorum Maxim. Here, UGT73CD1, a novel glycosyltransferase, was identified from I. tectorum. based on transcriptome analysis and functional identification. Phylogenetic analysis revealed that UGT73CD1 grouped into the clade of flavonoid 7-OH OGTs. Biochemical analysis showed that UGT73CD1 was able to glycosylate tectorigenin at 7-OH to produce tectoridin, and thus assigned as a 7-O-glycosyltransferase. In addition, it also possessed robust catalytic promiscuity toward 12 structurally diverse flavonoid scaffolds and 3, 4-dichloroaniline, resulting in forming O- and N-glycosides. This work will provide insights into efficient biosynthesis of structurally diverse flavonoid glycosides for drug discovery.

Genome-wide identification and analysis of AP2/ERF transcription factors related to camptothecin biosynthesis in Camptotheca acuminata.

Camptotheca acuminata produces camptothecin (CPT), a monoterpene indole alkaloid (MIA) that is widely used in the treatment of lung, colorectal, cervical, and ovarian cancers. Its biosynthesis pathway has attracted significant attention, but the regulation of CPT biosynthesis by the APETALA2/ethylene-responsive factor (AP2/ERF) transcription factors (TFs) remains unclear. In this study, a systematic analysis of the AP2/ERF TFs family in C. acuminata was performed, including phylogeny, gene structure, conserved motifs, and gene expression profiles in different tissues and organs (immature bark, cotyledons, young flower, immature fruit, mature fruit, mature leaf, roots, upper stem, and lower stem) of C. acuminata. A total of 198 AP2/ERF genes were identified and divided into five relatively conserved subfamilies, including AP2 (26 genes), DREB (61 genes), ERF (92 genes), RAV (18 genes), and Soloist (one gene). The combination of gene expression patterns in different C. acuminata tissues and organs, the phylogenetic tree, the co-expression analysis with biosynthetic genes, and the analysis of promoter sequences of key enzymes genes involved in CPT biosynthesis pathways revealed that eight AP2/ERF TFs in C. acuminata might be involved in CPT synthesis regulation, which exhibit relatively high expression levels in the upper stem or immature bark. Among these, four genes (CacAP2/ERF123, CacAP2/ERF125, CacAP2/ERF126, and CacAP2/ERF127) belong to the ERF-B2 subgroup; two genes (CacAP2/ERF149 and CacAP2/ERF152) belong to the ERF-B3 subgroup; and two more genes (CacAP2/ERF095 and CacAP2/ERF096) belong to the DREB-A6 subgroup. These results provide a foundation for future functional characterization of the AP2/ERF genes to enhance the biosynthesis of CPT compounds of C. acuminata.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Crocin Biosynthesis in Gardenia jasminoides Ellis (Rubiaceae).

Crocins, enriched in Gardenia jasminoides fruits, have a pharmacological activity against central nervous system diseases, cardiovascular diseases, and cancer cell growth. The biosynthesis of crocins has been widely explored, but its regulatory mechanism remains unknown. Here, the basic helix-loop-helix (bHLH) transcription factors related to crocin biosynthesis were systematically identified on the basis of the genome of G. jasminoides. A total of 95 GjbHLH transcription factor genes were identified, and their phylogenetic analysis indicated that they could be classified into 23 subfamilies. The combination of gene-specific bHLH expression patterns, the coexpression analysis of biosynthesis genes, and the analysis of promoter sequences in crocin biosynthesis pathways suggested that nine bHLHs in G. jasminoides might negatively regulate crocin biosynthesis. This study laid a foundation for understanding the regulatory mechanism of crocin biosynthesis and the improvement and breeding of G. jasminoides varieties.

Genome-wide analysis of LBD(lateral organ boundaries domain) gene family in Cannabis sativa of traditional Chinese medicine hemp seed / 中国中药杂志

LBD(lateral organ boundaries)transcription factors play an important role in the regulation of plant growth, development and secondary metabolism. In order to explore the function of LBD genes in cannabis, the Cannabis sativa genome and transcriptome were used to identify the C. sativa LBD gene family, and analyzed their expression patterns. Our results showed that the cannabis LBD contains 32 members, which were divided into two major categories, seven sub-families. Class Ⅰ was divided into 5 sub-families, named Class Ⅰ_a to Class Ⅰ_e, while Class Ⅱ was divided into 2 sub-families, including Class Ⅱ_a and Class Ⅱ_b. Analysis showed that the number of amino acids encoded LBDs was between 172 and 356, and the isoelectric point was between 4.92 and 9.43. The mole-cular weight of LBD was between 18 862.92 Da and 40 081.33 Da, and most members are located in the nucleus. Chromosome positioning of LBD showed that 32 members were unevenly distributed on 10 chromosomes of C. sativa LBD transcription factor domain, gene structure and motifs are relatively conservative, and the characteristics of different class members are similar. The upstream promoter region of the gene contains a variety of cis-acting elements related to plant hormones and environmental factors, C. sativa LBD genes have different expression patterns in the stems, leaves, and flowers of ZYS varieties(low tetrahydrocannabinol, high cannabidiol). The members of the LBD gene family are mainly expressed in the flowers and stems of ZYS varieties, while members expressed in the leaves very few; Class Ⅱ members CsLBD21 and CsLBD23 are expressed in flowers and stems, and CsLBD8 and CsLBD18 are expressed in flowers, stems and leaves. These genes may participate in the growth and development of cannabis and affect the biosynthesis of cannabinoids. This study laid the foundation for the subsequently functional research of the cannabis LBD gene family.

bHLH transcription factor SmbHLH92 negatively regulates biosynthesis of phenolic acids and tanshinones in Salvia miltiorrhiza.

Objective: Salvia miltiorrhiza is a valuable herbal medicine with tanshinone and phenolic acid as the main biological active ingredients. The biosynthetic regulation of these bioactive compounds is controlled by a set of transcription factors (TFs). The basic helix-loop-helix (bHLH) transcription factor plays an important role in various physiological and biochemical processes in plants. However, research on bHLH TFs regulating phenolic acid or tanshinone biosynthesis in S. miltiorrhiza is limited. Methods: qRT-PCR was used for gene expression analysis. The subcellular localization of SmbHLH92 was detected by SmbHLH92-GFP transient transformation into tobacco leaves, and its fluorescence was observed using a confocal laser scanning microscope. The transcriptional activity of SmbHLH92 was confirmed in the AH109 yeast strain. RNA interference hairy roots of SmbHLH92-RNAi transgenic lines were obtained through Agrobacterium-mediated genetic transformation. Ultra performance liquid chromatography (UPLC) was used to detect the changes of phenolic acids and tanshinones. Results: SmbHLH92 is a bHLH transcription factor that is highly expressed in the root and phloem of S. miltiorrhiza. The subcellular localization and transcriptional activity of SmbHLH92 indicated that SmbHLH92 was located in the nucleus and may be a transcription factor. RNA interference (RNAi) of SmbHLH92 in hairy roots of S. miltiorrhiza significantly increased the accumulation of phenolic acid and tanshinone. Quantitative RT-PCR (RT-qPCR) analysis showed the transcription level of genes encoding the key enzymes involved in the phenolic acid and tanshinone biosynthetic pathways was increased in the hairy roots of the SmbHLH92-RNAi transgenic line, comparing with the control line. Conclusion: These data indicate that SmbHLH92 is a negative regulator involved in the regulation of phenolic acid and tanshinone biosynthesis in S. miltiorrhiza.

Genomic survey of bZIP transcription factor genes related to tanshinone biosynthesis in Salvia miltiorrhiza.

Tanshinones are a class of bioactive components in the traditional Chinese medicine Salvia miltiorrhiza, and their biosynthesis and regulation have been widely studied. Current studies show that basic leucine zipper (bZIP) proteins regulate plant secondary metabolism, growth and developmental processes. However, the bZIP transcription factors involved in tanshinone biosynthesis are unknown. Here, we conducted the first genome-wide survey of the bZIP gene family and analyzed the phylogeny, gene structure, additional conserved motifs and alternative splicing events in S. miltiorrhiza. A total of 70 SmbZIP transcription factors were identified and categorized into 11 subgroups based on their phylogenetic relationships with those in Arabidopsis. Moreover, seventeen SmbZIP genes underwent alternative splicing events. According to the transcriptomic data, the SmbZIP genes that were highly expressed in the Danshen root and periderm were selected. Based on the prediction of bZIP binding sites in the promoters and the co-expression analysis and co-induction patterns in response to Ag+ treatment via quantitative real-time polymerase chain reaction (qRT-PCR), we concluded that SmbZIP7 and SmbZIP20 potentially participate in the regulation of tanshinone biosynthesis. These results provide a foundation for further functional characterization of the candidate SmbZIP genes, which have the potential to increase tanshinone production.

Transcriptome-Guided Mining of Genes Involved in Crocin Biosynthesis.

Gardenia jasminoides is used in traditional Chinese medicine and has drawn attention as a rich source of crocin, a compound with reported activity against various cancers, depression and cardiovascular disease. However, genetic information on the crocin biosynthetic pathway of G. jasminoides is scarce. In this study, we performed a transcriptome analysis of the leaves, green fruits, and red fruits of G. jasminoides to identify and predict the genes that encode key enzymes responsible for crocin production, compared with Crocus sativus. Twenty-seven putative pathway genes were specifically expressed in the fruits, consistent with the distribution of crocin in G. jasminoides. Twenty-four of these genes were reported for the first time, and a novel CCD4a gene was predicted that encodes carotenoid cleavage dioxygenase leading to crocin synthesis, in contrast to CCD2 of C. sativus. In addition, 6 other candidate genes (ALDH12, ALDH14, UGT94U1, UGT86D1, UGT71H4, and UGT85K18) were predicted to be involved in crocin biosynthesis following phylogenetic analysis and different gene expression profiles. Identifying the genes that encode key enzymes should help elucidate the crocin biosynthesis pathway.

Genome-wide characterisation and analysis of bHLH transcription factors related to tanshinone biosynthesis in Salvia miltiorrhiza.

Salvia miltiorrhiza Bunge (Labiatae) is an emerging model plant for traditional medicine, and tanshinones are among the pharmacologically active constituents of this plant. Although extensive chemical and pharmaceutical studies of these compounds have been performed, studies on the basic helix-loop-helix (bHLH) transcription factors that regulate tanshinone biosynthesis are limited. In our study, 127 bHLH transcription factor genes were identified in the genome of S. miltiorrhiza, and phylogenetic analysis indicated that these SmbHLHs could be classified into 25 subfamilies. A total of 19 sequencing libraries were constructed for expression pattern analyses using RNA-Seq. Based on gene-specific expression patterns and up-regulated expression patterns in response to MeJA treatment, 7 bHLH genes were revealed as potentially involved in the regulation of tanshinone biosynthesis. Among them, the gene expression of SmbHLH37, SmbHLH74 and SmbHLH92 perfectly matches the accumulation pattern of tanshinone biosynthesis in S. miltiorrhiza. Our results provide a foundation for understanding the molecular basis and regulatory mechanisms of bHLH transcription factors in S. miltiorrhiza.

Chromosome-level genome map provides insights into diverse defense mechanisms in the medicinal fungus Ganoderma sinense.

Fungi have evolved powerful genomic and chemical defense systems to protect themselves against genetic destabilization and other organisms. However, the precise molecular basis involved in fungal defense remain largely unknown in Basidiomycetes. Here the complete genome sequence, as well as DNA methylation patterns and small RNA transcriptomes, was analyzed to provide a holistic overview of secondary metabolism and defense processes in the model medicinal fungus, Ganoderma sinense. We reported the 48.96 Mb genome sequence of G. sinense, consisting of 12 chromosomes and encoding 15,688 genes. More than thirty gene clusters involved in the biosynthesis of secondary metabolites, as well as a large array of genes responsible for their transport and regulation were highlighted. In addition, components of genome defense mechanisms, namely repeat-induced point mutation (RIP), DNA methylation and small RNA-mediated gene silencing, were revealed in G. sinense. Systematic bioinformatic investigation of the genome and methylome suggested that RIP and DNA methylation combinatorially maintain G. sinense genome stability by inactivating invasive genetic material and transposable elements. The elucidation of the G. sinense genome and epigenome provides an unparalleled opportunity to advance our understanding of secondary metabolism and fungal defense mechanisms.

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.

Danshen, Salvia miltiorrhiza Bunge, is one of the most widely used herbs in traditional Chinese medicine, wherein its rhizome/roots are particularly valued. The corresponding bioactive components include the tanshinone diterpenoids, the biosynthesis of which is a subject of considerable interest. Previous investigations of the S. miltiorrhiza transcriptome have relied on short-read next-generation sequencing (NGS) technology, and the vast majority of the resulting isotigs do not represent full-length cDNA sequences. Moreover, these efforts have been targeted at either whole plants or hairy root cultures. Here, we demonstrate that the tanshinone pigments are produced and accumulate in the root periderm, and apply a combination of NGS and single-molecule real-time (SMRT) sequencing to various root tissues, particularly including the periderm, to provide a more complete view of the S. miltiorrhiza transcriptome, with further insight into tanshinone biosynthesis as well. In addition, the use of SMRT long-read sequencing offered the ability to examine alternative splicing, which was found to occur in approximately 40% of the detected gene loci, including several involved in isoprenoid/terpenoid metabolism.

[Progress in the study of Velvet and LaeA proteins and their relation to the development and bioactive compounds in medicinal fungi].

The medicinal fungi, which are of great importance in traditional medicine, are facing the problems of wild resources scarcity and low concentration of bioactive compounds. Velvet family and LaeA global regulator play a vital role in secondary metabolism and developmental programs, which are found in a wide variety of fungi ranging from Chytridiomycota to Basidiomycota. This review elaborates the structures and functions between Velvet family and LaeA protein. The Velvet family which shares the Velvet protein domain, including VeA (Velvet), VelB (Velvet like B), VosA (viability of spores A) and VelC (Velvet like C), acts on the regulation function is secondary metabolism and developmental programs such as asexual and sexual development. Furthermore, the function is affected by environmental factors such as light and temperature. LaeA protein which owns S-adenosylmethionine-dependent methyltransferase domain, coordinately regulates development and secondary metabolism by regulating and modifying the Velvet proteins. The regulation of LaeA is mediated by light receptor proteins. Therefore, clarifying the mechanism of Velvet and LaeA proteins in medicinal fungi will pave the way for nurturing medicinal fungi and improving production of bioactive compounds.