Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

Aristolochic acids (AAs) and their derivatives are present in multiple Aristolochiaceae species that have been or are being used as medicinal materials. During the past decades, AAs have received increasing attention because of their nephrotoxicity and carcinogenicity. Elimination of AAs from medicinal materials using biotechnological approaches is important for improving medication safety. However, it has not been achieved because of the limited information available on AA biosynthesis. Here, we report a high-quality, reference-grade genome assembly of the AA-containing vine Aristolochia contorta. The total size of the assembly is 209.27 Mb, and it is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution, and 4DTv suggest an absence of whole-genome duplication (WGD) events in Aristolochia contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic, and metabolic data, pathways and candidate genes were proposed for benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta. Five O-methyltransferase genes, including AcOMT1–3, AcOMT5, and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. They lay a solid foundation for further elucidation of AA biosynthesis and regulation and for the molecular breeding of Aristolochiaceae medicinal materials.

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs in plants. They play critical functions in various biological processes during plant growth and development. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant with significant medicinal, economic, and academic values. In order to elucidate the role of miRNAs in S. miltiorrhiza, six small RNA libraries from mature roots, young roots, stems, mature leaves, young leaves and flowers of S. miltiorrhiza and one degradome library from mixed tissues were constructed. A total of 184 miRNA precursors, generating 137 known and 49 novel miRNAs, were genome-widely identified. The identified miRNAs were predicted to play diversified regulatory roles in plants through regulating 891 genes. qRT-PCR and 5′ RLM-RACE assays validated the negative regulatory role of smi-miR159a in SmMYB62, SmMYB78, and SmMYB80. To elucidate the function of smi-miR159a in bioactive compound biosynthesis, smi-miR159a transgenic hairy roots were generated and analyzed. The results showed that overexpression of smi-miR159a caused a significant decrease in rosmarinic acid and salvianolic acid B contents. qRT-PCR analysis showed that the targets of smi-miR159a, including SmMYB62, SmMYB78, and SmMYB80, were significantly down-regulated, accompanied by the down-regulation of SmPAL1, SmC4H1, Sm4CL1, SmTAT1, SmTAT3, SmHPPR1, SmRAS, and SmCYP98A14 genes involved in phenolic acid biosynthesis. It suggests that smi-miR159a is a significant negative regulator of phenolic acid biosynthesis in S. miltiorrhiza.

Alkaline/neutral invertases (NINVs), which irreversibly catalyze the hydrolysis of sucrose into fructose and glucose, play crucial roles in carbohydrate metabolism and plant development. Comprehensive insights into NINV genes are lacking in Salvia miltiorrhiza, a well-known traditional Chinese medicinal (TCM) plant with significant medicinal and economic value. Through genome-wide prediction, nine putative SmNINV genes, termed SmNINV1-SmNINV9, were identified. Integrated analysis of gene structures, sequence features, conserved domains, conserved motifs and phylogenetic trees revealed the conservation and divergence of SmNINVs. The identified SmNINVs were differentially expressed in roots, stems, leaves, flowers, and different root tissues. They also responded to drought, salicylic acid, yeast extract, and methyl jasmonate treatments. More importantly, computational prediction and experimental validation showed that SmNINV3 and SmNINV4 were targets of Smi-miR399, a conserved miRNA previously shown to affect Pi uptake and translocation through the cleavage of PHOSPHATE2 (PHO2). Consistently, analysis of 43 NINV genes and 26 miR399 sequences from Arabidopsis thaliana, Populus trichocarpa, Manihot esculenta, and Solanum lycopersicum showed that various AtNINV, PtNINV, MeNINV, and SlNINV genes were regulated by miR399. It indicates that the miR399-NINV module exists widely in plants. Furthermore, Smi-miR399 also cleaved SmPHO2 transcripts in S. miltiorrhiza, suggesting the complexity of NINVs, PHO2, and miR399 networks.

Aristolochia contorta Bunge has been widely used as traditional Chinese medicine materials. However, its utility faces a great challenge due to the presence of aristolochic acids (AAs), a class of benzylisoquinoline alkaloid (BIA) derivatives. The first step in BIA skeleton formation is catalysis by norcoclaurine synthase (NCS). To gain knowledge of BIA and AA biosynthesis in A. contorta, genome-wide characterizations of NCS genes were carried out. This resulted in the identification of 15 A. contorta NCSs, namely, AcNCS1–AcNCS15. The AcNCS1–AcNCS8 proteins contained one catalytic domain, whereas the AcNCS9–AcNCS15 proteins had two. Phylogenetic analysis shows that AcNCS proteins can be classified into two clades. Gene expression analysis shows that five AcNCSs, including AcNCS2, AcNCS4, AcNCS5, AcNCS14, and AcNCS15, exhibited relatively high expression in roots and flowers, where norcoclaurine accumulated. An enzyme catalytic activity assay shows that all five of the AcNCSs can catalyze norcoclaurine formation with AcNCS14 and AcNCS15, exhibiting higher catalytic efficiency. Precolumn derivatization analysis shows that the formed norcoclaurine included (S)- and (R)-norcoclaurine, with more (S)-configuration. The results provide useful information for further understanding BIA and AA biosynthesis in A. contorta and for AA elimination and bioactive compound improvement in AA-containing medicinal materials.