New α-Glucosidase Inhibitors from the Whole Plant of Hypericum beanii Based on Ligand Fishing and Molecular Networking Analysis.

In this work, a new rapid and targeted method for screening α-glucosidase inhibitors from Hypericum beanii was developed and verified. Ten new polycyclic polyprenylated acylphloroglucinols (PPAPs), hyperlagarol A-J (1-10), and nine known PPAPs (11-19) were obtained from H. beanii. Their structures were identified by using comprehensive analyses involving mass spectrometry, ultraviolet spectroscopy, infrared spectroscopy, nuclear magnetic resonance spectroscopy, and electron capture dissociation calculations. 1 and 2 are two new rare 2,3-seco-spirocyclic PPAPs, 3 and 4 are two novel 12,13-seco-spirocyclic PPAPs, 5 and 6 are two novel spirocyclic PPAPs, 7 and 8 are two new unusual spirocyclic PPAPs with complex bridged ring systems, and 9 and 10 are two novel nonspirocyclic PPAPs. α-GC inhibitory activities of all isolated compounds were tested. Most of them displayed inhibitory activities against α-glucosidase, with the IC50 values ranging from 6.85 ± 0.65 to 112.5 ± 9.03 µM. Moreover, the inhibitory type and mechanism of the active compounds were further analyzed using kinetic studies and molecular docking.

Discovering a uniform functional trade-off of the CBC-type 2,3-oxidosqualene cyclases and deciphering its chemical logic.

Many functionally promiscuous plant 2,3-oxidosqualene cyclases (OSCs) have been found, but complete functional reshaping is rarely reported. In this study, we have identified two new plant OSCs: a unique protostadienol synthase (AoPDS) and a common cycloartenol synthase (AoCAS) from Alisma orientale (Sam.) Juzep. Multiscale simulations and mutagenesis experiments revealed that threonine-727 is an essential residue responsible for protosta-13 (17),24-dienol biosynthesis in AoPDS and that the F726T mutant completely reshapes the native function of AoCAS into a PDS function to yield almost exclusively protosta-13 (17),24-dienol. Unexpectedly, various native functions were uniformly reshaped into a PDS function by introducing the phenylalanine → threonine substitution at this conserved position in other plant and non-plant chair-boat-chair-type OSCs. Further computational modeling elaborated the trade-off mechanisms of the phenylalanine → threonine substitution that leads to the PDS activity. This study demonstrates a general strategy for functional reshaping by using a plastic residue based on the decipherment of the catalytic mechanism.

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.

Oxidosqualene Cyclases Involved in the Biosynthesis of Diverse Triterpenes in Camellia sasanqua.

Camellia sasanqua is an important economic plant that is rich in lipophilic triterpenols with pharmacological activities including antiallergic, anti-inflammatory, and anticancer activities. However, the key enzymes related to triterpene biosynthesis have seldom been studied in C. sasanqua. Oxidosqualene cyclases (OSCs) are the rate-limiting enzymes related to triterpene biosynthesis. In this study, seven putative OSC genes (CsOSC1-7) were mined from the C. sasanqua transcriptome. Six CsOSCs were characterized for the biosynthesis of diverse triterpene skeletons, including α-amyrin, ß-amyrin, δ-amyrin, dammarenediol-II, ψ-taraxasterol, taraxasterol, and cycloartenol by the heterologous expression system. CsOSC3 was a multiple functional α-amyrin synthase. Three key residues, Trp260, Tyr262, and Phe415, are critical to the catalytic performance of CsOSC3 judging from the results of molecular docking and site-directed mutagenesis. These findings provide important insights into the biosynthesis pathway of triterpenes in C. sasanqua.

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.

Six Unusual Meroterpenoids from the Leaves of Psidium guajava L. and Their PTP1B Inhibitory Activities.

Six unusual meroterpenoids, psidiguajadiol A-J (1-6), and three known meroterpenoids (7-9) were isolated from the leaves of Psidium guajava L. Compounds 2-6 represent the first examples of 6/8-formyl-5,7-dihydroxy-4-phenylchromane-coupled sesquiterpenoids. The structures of the undescribed compounds, including their absolute configurations, were elucidated by spectroscopic analyses, X-ray diffraction, and computational calculations. Compounds 3, 4, and 6 exhibited inhibitory activities against PTP1B with IC50 values of 9.83, 18.52, and 16.87 µM, respectively. In light of these findings, we performed molecular docking studies to predict their inhibition mechanisms at the atomic level.

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.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants.

BACKGROUND: Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits. RESULTS: Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene. CONCLUSIONS: This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

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.

An Unexpected Oxidosqualene Cyclase Active Site Architecture in the Iris tectorum Multifunctional α-Amyrin Synthase.

Ordered polycyclization catalyzed by oxidosqualene synthases (OSCs) morph a common linear precursor into structurally complex and diverse triterpene scaffolds with varied bioactivities. We identified three OSCs from Iris tectorum. ItOSC2 is a rare multifunctional α-amyrin synthase. Sequence comparisons, site-directed mutagenesis and multiscale simulations revealed that three spatially clustered residues, Y531/L256/L258 form an unusual Y-LL triad at the active site, replacing the highly conserved W-xY triad occurring in other amyrin synthases. The discovery of this unprecedented active site architecture in ItOSC2 underscores the plasticity of terpene cyclase catalytic mechanisms and opens new avenues for protein engineering towards custom designed OSCs.

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 Characterization of WRKY Transcription Factors Related to Andrographolide Biosynthesis in Andrographis paniculata.

Andrographolide, which is enriched in the leaves of Andrographis paniculata, has been known as "natural antibiotic" due to its pharmacological activities such as anti-inflammatory, antimicrobial and antioxidant effects. Several key enzymes in andrographolide biosynthetic pathway have been studied since the genome sequences were released, but its regulatory mechanism remains unknown. WRKY transcription factors proteins have been reported to regulate plant secondary metabolism, development as well as biotic and abiotic stresses. Here, WRKY transcription factors related to andrographolide biosynthesis were systematically identified, including sequences alignment, phylogenetic analysis, chromosomal distribution, gene structure, conserved motifs, synteny, alternative splicing event and Gene ontology (GO) annotation. A total of 58 WRKYs were identified in Chuanxinlian genome and phylogenetically classified into three groups. Moreover, nine WRKY genes underwent alternative splicing events. Furthermore, the combination of binding site prediction, gene-specific expression patterns, and phylogenetic analysis suggested that 7 WRKYs (ApWRKY01, ApWRKY08, ApWRKY12, ApWRKY14, ApWRKY19, ApWRKY20, and ApWRKY50) might regulate andrographolide biosynthesis. This study laid a foundation for understanding the regulatory mechanism of andrographolide biosynthesis and the improvement and breeding of Andrographis paniculata 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.

The AP2/ERF transcription factor SmERF128 positively regulates diterpenoid biosynthesis in Salvia miltiorrhiza.

KEY MESSAGE: The novel AP2/ERF transcription factor SmERF128 positively regulates diterpenoid tanshinone biosynthesis by activating the expression of SmCPS1, SmKSL1, and SmCYP76AH1 in Salvia miltiorrhiza. Certain members of the APETALA2/ethylene-responsive factor (AP2/ERF) family regulate plant secondary metabolism. Although it is clearly documented that AP2/ERF transcription factors (TFs) are involved in sesquiterpenoid biosynthesis, the regulation of diterpenoid biosynthesis by AP2/ERF TFs remains elusive. Here, we report that the novel AP2/ERF TF SmERF128 positively regulates diterpenoid tanshinone biosynthesis in Salvia miltiorrhiza. Overexpression of SmERF128 increased the expression levels of copalyl diphosphate synthase 1 (SmCPS1), kaurene synthase-like 1 (SmKSL1) and cytochrome P450 monooxygenase 76AH1 (SmCYP76AH1), whereas their expression levels were decreased when SmERF128 was silenced. Accordingly, the content of tanshinone was reduced in SmERF128 RNA interference (RNAi) hairy roots and dramatically increased in SmERF128 overexpression hairy roots, as demonstrated through Ultra Performance Liquid Chromatography (UPLC) and Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) analysis. Furthermore, SmERF128 activated the expression of SmCPS1, SmKSL1, and SmCYP76AH1 by binding to the GCC box, and to the CRTDREHVCBF2 (CBF2) and RAV1AAT (RAA) motifs within their promoters during in vivo and in vitro assays. Our findings not only reveal the molecular basis of how the AP2/ERF transcription factor SmERF128 regulates diterpenoid biosynthesis, but also provide useful information for improving tanshinone production through genetic engineering.

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 analysis of auxin response factor gene family members in medicinal model plant Salvia miltiorrhiza.

Auxin response factors (ARFs) can function as transcriptional activators or repressors to regulate the expression of auxin response genes by specifically binding to auxin response elements (AuxREs) during plant development. Based on a genome-wide strategy using the medicinal model plant Salvia miltiorrhiza, 25 S. miltiorrhiza ARF (SmARF) gene family members in four classes (class Ia, IIa, IIb and III) were comprehensively analyzed to identify characteristics including gene structures, conserved domains, phylogenetic relationships and expression patterns. In a hybrid analysis of the phylogenetic tree, microRNA targets, and expression patterns of SmARFs in different organs, root tissues, and methyl jasmonate or indole-3-acetic acid treatment conditions, we screened for candidate SmARFs involved in various developmental processes of S. miltiorrhiza Based on this analysis, we predicted that SmARF25, SmARF7, SmARF16 and SmARF20 are involved in flower, leaf, stem and root development, respectively. With the further insight into the targets of miR160 and miR167, specific SmARF genes in S. miltiorrhiza might encode products that participate in biological processes as described for ARF genes in Arabidopsis Our results provide a foundation for understanding the molecular basis and regulatory mechanisms of SmARFs in S. miltiorrhiza.