Genome-Wide Identification of the AP2/ERF Gene Family Involved in Active Constituent Biosynthesis in

Transcription factors in tanshinones: Emerging mechanisms of transcriptional regulation

Transcription factors play a crucial role in the biosynthesis of tanshinones, which are significant secondary metabolites derived from Salvia miltiorrhiza, commonly known as Danshen. These compounds have extensive pharmacological properties, including anti-inflammatory and cardioprotective effects. This review delves into the roles of various transcription factor families, such as APETALA2/ethylene response factor, basic helix-loop-helix, myeloblastosis, basic leucine zipper, and WRKY domain-binding protein, in regulating the biosynthetic pathways of tanshinones. We discuss the emerging mechanisms by which these transcription factors influence the synthesis of tanshinones, both positively and negatively, by directly regulating gene expression or forming complex regulatory networks. Additionally, the review highlights the potential applications of these insights in enhancing tanshinone production through genetic and metabolic engineering, setting the stage for future advancements in medicinal plant research.

Genome-wide identification of CRF gene family members in four rice subspecies and expression analysis of OsCRF members in response to cold stress at seedling stage

Cytokinin Response Factors (CRFs) play a crucial role in plant growth and development, hormone signaling, and responses to biotic and abiotic stresses. However, there have been no reports on CRF genes in rice until now. We analyzed the CRF families in four rice subspecies: cultivated rice Oryza sativa Japonica Group, Oryza sativa Indica Group, and Oryza sativa (circum-Aus1 var. N22), as well as wild rice Oryza rufipogon. We identified 7, 6, 6, and 7 CRF in their genomes, respectively, distributed across different chromosomes. The protein motifs and gene structures of CRF in these four types of rice show high conservation. Cis-regulatory element analysis revealed that the promoter regions of the CRF contain numerous hormone and stress-related elements. The number of CRF in these four types of rice is not influenced by gene duplication. The expression pattern showed that OsCRF exhibit significant tissue-specific expression. The qRT-PCR results showed that OsCRF strongly responded to low-temperature stress and can be induced by melatonin and cytokinin to increase expression levels. In addition, the nuclear localisation of OsCRF4/5 was confirmed to be as predicted. The results above will provide a foundation for further and deeper investigation of CRFs.

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.

DoAP2/ERF89 activated the terpene synthase gene DoPAES in Dendrobium officinale and participated in the synthesis of β-patchoulene

Dendrobium officinale Kimura et Migo is a tonic plant that has both ornamental and medicinal properties. Terpenoids are significant and diverse secondary metabolites in plants, and are one of the important natural active ingredients in D. officinale. The AP2/ERF gene family plays a major role in primary and secondary metabolism. However, the AP2/ERF transcription factor family has not been identified in D. officinale, and it is unclear if it is involved in the regulation of terpenoid biosynthesis. This study identified a sesquiterpene synthetase-β-patchoulene synthase (DoPAES) using transcriptome and terpenic metabolic profile analyses. A total of 111 members of the AP2/ERF family were identified through the whole genome of D. officinale. The tissue-specific expression and gene co-expression pattern of the DoAP2/ERF family members were analyzed. The results showed that the expression of DoPAES was highly correlated with the expression of DoAP2/ERF89 and DoAP2/ERF47. The yeast one-hybrid (Y1H) assays and dual-luciferase experiments demonstrated that DoAP2/ERF89 and DoAP2/ERF47 could regulate the expression of DoPAES. The transcriptional regulatory effects were examined using homologous transient expression of DoAP2/ERF89 in protocorms of D. officinale. DoAP2/ERF89 positively regulated the biosynthesis of β-patchoulene. This study showed that DoAP2/ERF89 can bind to the promoter region of DoPAES to control its expression and further regulate the biosynthesis of β-patchoulene in D. officinale. These results provide new insights on the regulation of terpenoid biosynthesis.

The SmERF1b-like regulates tanshinone biosynthesis in Salvia miltiorrhiza hairy root

The ethylene response factor family genes are involved in the regulation of secondary metabolism in Salvia miltiorrhiza, but the mechanism underlying this regulation remains elusive. In the present study, based on the cDNA library of S. miltiorrhiza, an AP2/ERF gene was cloned and named SmERF1b-like. This gene exhibited a significant response to exogenous ethylene supply, such that ethylene remarkably upregulated SmERF1b-like expression levels in the leaves of S. miltiorrhiza. Subcellular localization showed that SmERF1b-like is located in the nucleus. Furthermore, SmERF1b-like showed a binding affinity with a GCC-box motif in the promoter region of genes associated with tanshinone biosynthesis in S. miltiorrhiza. Overexpression of SmERF1b-like in hairy roots of S. miltiorrhiza substantially upregulated SmCPS1 and SmKSL1 expression levels, resulting in increased biosynthesis of tanshinone I and cryptotanshinone contents. This finding provides valuable theoretical support for the utilization of a plant genetic engineering strategy to enhance S. miltiorrhiza resources.

Genome-wide identification of Apetala2 gene family in Hypericum perforatum L and expression profiles in response to different abiotic and hormonal treatments

The Apetala2 (AP2) gene family of transcription factors (TFs) play important functions in plant development, hormonal response, and abiotic stress. To reveal the biological functions and the expression profiles of AP2 genes in Hypericum perforatum, genome-wide identification of HpAP2 family members was conducted.

Methods

We identified 21 AP2 TFs in H. perforatum using bioinformatic methods; their physical and chemical properties, gene structures, conserved motifs, evolutionary relationships, cis-acting elements, and expression patterns were investigated.

Results

We found that based on the structural characteristics and evolutionary relationships, the HpAP2 gene family can be divided into three subclasses: euANT, baselANT, and euAP2. A canonical HpAP2 TF shared a conserved protein structure, while a unique motif 6 was found in HpAP2_1, HpAP2_4, and HpAP2_5 from the euANT subgroup, indicating potential biological and regulatory functions of these genes. Furthermore, a total of 59 cis-acting elements were identified, most of which were associated with growth, development, and resistance to stress in plants. Transcriptomics data showed that 57.14% of the genes in the AP2 family were differentially expressed in four organs. For example, HpAP2_18 was specifically expressed in roots and stems, whereas HpAP2_17 and HpAP2_11 were specifically expressed in leaves and flowers, respectively. HpAP2_5, HpAP2_11, and HpAP2_18 showed tissue-specific expression patterns and responded positively to hormones and abiotic stresses.

Conclusion

These results demonstrated that the HpAP2 family genes are involved in diverse developmental processes and generate responses to abiotic stress conditions in H. perforatum. This article, for the first time, reports the identification and expression profiles of the AP2 family genes in H. perforatum, laying the foundation for future functional studies with these genes.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

Jasmonic acid regulates the biosynthesis of medicinal metabolites via the JAZ9-MYB76 complex in Salvia miltiorrhiza

Jasmonic acid (JA) signaling pathway plays an important role in tanshinone and phenolic acid biosynthesis in Salvia miltiorrhiza. However, the specific regulatory mechanism remains largely unclear. Previous work showed that a JASMONATE ZIM-domain (JAZ) protein, SmJAZ9, acted as a repressor of tanshinone production in S. miltiorrhiza. In this study, we revealed that SmJAZ9 reduced both phenolic acid accumulation and related biosynthetic gene expression, confirming that SmJAZ9 also negatively affected phenolic acid biosynthesis. Then, we identified a novel MYB transcription factor, SmMYB76, which interacted with SmJAZ9. SmMYB76 repressed phenolic acid biosynthesis by directly downregulating SmPAL1, Sm4CL2, and SmRAS1. Further investigation demonstrated that JA mediated phenolic acids biosynthesis via SmJAZ9-SmMYB76 complex. Taken together, these findings state the molecular mechanism that SmJAZ9-SmMYB76 regulated phenolic acid biosynthesis at the transcriptional and protein levels, which provided new insights into JA signaling pathway regulating plant metabolism.

Genome-Wide Identification of the ERF Transcription Factor Family for Structure Analysis, Expression Pattern, and Response to Drought Stress in Populus alba × Populus glandulosa

The Ethylene Responsive Factor (ERF) transcription factor family is important for regulating plant growth and stress responses. Although the expression patterns of ERF family members have been reported in many plant species, their role in Populus alba × Populus glandulosa, an important model plant for forest research, remains unclear. In this study, we identified 209 PagERF transcription factors by analyzing the P. alba × P. glandulosa genome. We analyzed their amino acid sequences, molecular weight, theoretical pI (Isoelectric point), instability index, aliphatic index, grand average of hydropathicity, and subcellular localization. Most PagERFs were predicted to localize in the nucleus, with only a few PagERFs localized in the cytoplasm and nucleus. Phylogenetic analysis divided the PagERF proteins into ten groups, Class I to X, with those belonging to the same group containing similar motifs. Cis-acting elements associated with plant hormones, abiotic stress responses, and MYB binding sites were analyzed in the promoters of PagERF genes. We used transcriptome data to analyze the expression patterns of PagERF genes in different tissues of P. alba × P. glandulosa, including axillary buds, young leaves, functional leaves, cambium, xylem, and roots, and the results indicated that PagERF genes are expressed in all tissues of P. alba × P. glandulosa, especially in roots. Quantitative verification results were consistent with transcriptome data. When P. alba × P. glandulosa seedlings were treated with 6% polyethylene glycol 6000 (PEG6000), the results of RT-qRCR showed that nine PagERF genes responded to drought stress in various tissues. This study provides a new perspective on the roles of PagERF family members in regulating plant growth and development, and responses to stress in P. alba × P. glandulosa. Our study provides a theoretical basis for ERF family research in the future.

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%).

Methyl Jasmonate Activates the 2C Methyl-D-erithrytol 2,4-cyclodiphosphate Synthase Gene and Stimulates Tanshinone Accumulation in Salvia miltiorrhiza Solid Callus Cultures

The present study characterizes the 5′ regulatory region of the SmMEC gene. The isolated fragment is 1559 bp long and consists of a promoter, 5′UTR and 31 nucleotide 5′ fragments of the CDS region. In silico bioinformatic analysis found that the promoter region contains repetitions of many potential cis-active elements. Cis-active elements associated with the response to methyl jasmonate (MeJa) were identified in the SmMEC gene promoter. Co-expression studies combined with earlier transcriptomic research suggest the significant role of MeJa in SmMEC gene regulation. These findings were in line with the results of the RT-PCR test showing SmMEC gene expression induction after 72 h of MeJa treatment. Biphasic total tanshinone accumulation was observed following treatment of S. miltiorrhiza solid callus cultures with 50–500 μM methyl jasmonate, with peaks observed after 10–20 and 50–60 days. An early peak of total tanshinone concentration (0.08%) occurred after 20 days of 100 μM MeJa induction, and a second, much lower one, was observed after 50 days of 50 μM MeJa stimulation (0.04%). The dominant tanshinones were cryptotanshinone (CT) and dihydrotanshinone (DHT). To better understand the inducing effect of MeJa treatment on tanshinone biosynthesis, a search was performed for methyl jasmonate-responsive cis-active motifs in the available sequences of gene proximal promoters associated with terpenoid precursor biosynthesis. The results indicate that MeJa has the potential to induce a significant proportion of the presented genes, which is in line with available transcriptomic and RT-PCR data.

Function Analysis of the ERF and DREB Subfamilies in Tomato Fruit Development and Ripening

APETALA2/ethylene responsive factors (AP2/ERF) are unique regulators in the plant kingdom and are involved in the whole life activity processes such as development, ripening, and biotic and abiotic stresses. In tomato (Solanum lycopersicum), there are 140 AP2/ERF genes; however, their functionality remains poorly understood. In this work, the 14th and 19th amino acid differences in the AP2 domain were used to distinguish DREB and ERF subfamily members. Even when the AP2 domain of 68 ERF proteins from 20 plant species and motifs in tomato DREB and ERF proteins were compared, the binding ability of DREB and ERF proteins with DRE/CRT and/or GCC boxes remained unknown. During fruit development and ripening, the expressions of 13 DREB and 19 ERF subfamily genes showed some regular changes, and the promoters of most genes had ARF, DRE/CRT, and/or GCC boxes. This suggests that these genes directly or indirectly respond to IAA and/or ethylene (ET) signals during fruit development and ripening. Moreover, some of these may feedback regulate IAA or ET biosynthesis. In addition, 16 EAR motif-containing ERF genes in tomato were expressed in many organs and their total transcripts per million (TPM) values exceeded those of other ERF genes in most organs. To determine whether the EAR motif in EAR motif-containing ERF proteins has repression function, their EAR motifs were retained or deleted in a yeast one-hybrid (YIH) assay. The results indicate that most of EAR motif-containing ERF proteins lost repression activity after deleting the EAR motif. Moreover, some of these were expressed during ripening. Thus, these EAR motif-containing ERF proteins play vital roles in balancing the regulatory functions of other ERF proteins by completing the DRE/CRT and/or GCC box sites of target genes to ensure normal growth and development in tomato.

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

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.

Transcription Factor: A Powerful Tool to Regulate Biosynthesis of Active Ingredients in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a common Chinese herbal medicine, and its major active ingredients are phenolic acids and tanshinones, which are widely used to treat vascular diseases. However, the wild form of S. miltiorrhiza possess low levels of these important pharmaceutical agents; thus, improving their levels is an active area of research. Transcription factors, which promote or inhibit the expressions of multiple genes involved in one or more biosynthetic pathways, are powerful tools for controlling gene expression in biosynthesis. Several families of transcription factors have been reported to participate in regulating phenolic acid and tanshinone biosynthesis and influence their accumulation. This review summarizes the current status in this field, with focus on the transcription factors which have been identified in recent years and their functions in the biosynthetic regulation of phenolic acids and tanshinones. Otherwise, the new insight for further research is provided. Finally, the application of the biosynthetic regulation of active ingredients by the transcription factors in S. miltiorrhiza are discussed, and new insights for future research are explored.

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.

Whole-genome characterization of Rosa chinensis AP2/ERF transcription factors and analysis of negative regulator RcDREB2B in Arabidopsis

BACKGROUND

Rose (Rosa chinensis) is a traditional famous flower with valuable ornamental characteristics. However, drought stress restricts its growth and development, leading to an abnormal phenotype. One of the main transcription factor (TF) protein groups in the plant kingdom are the APETALA2/ethylene-responsive factor (AP2/ERF) proteins and are potentially involved in the growth and stress responses of various plants.

RESULTS

Our investigation mainly focused on exploring the genome of rose and thereby we discovered 135 apparent AP2/ERF TFs. Phylogenic analyses revealed that RcAP2/ERF genes are categorized into DREB, Soloist, AP2, and ERF subfamilies, and are further classified these into 17 groups, with the same as Malus domestica and Arabidopsis thaliana. The analysis of the gene structure revealed that the introns ranged from 0 to 9 in number. Pattern examination demonstrated that the RcAP2/ERF predominantly consists of typical AP2 domains, of which the 2nd motif is the most ubiquitous. Distributions of cis-acting elements indicated that members of the AP2/ERF family are frequently involved in growth and development, phytohormone and stress response in rose species. Also, the distribution mapping of the rose chromosomes indicated that AP2/ERF class genes are dispersed among all seven chromosomes. Additionally, we isolated a novel DREB A2 subgroup gene and named it RcDREB2B. Subsequently, the RcDREB2B transcript accumulation was repressed under the mild and severe drought stress in the root samples of rose. RcDREB2B was targeted to the nucleus and exhibited transactivation in yeast cells. The overexpression of RcDREB2B was found to promote sensitivity to a higher salt concentration, ABA, and PEG at the germination and post-germination stages. Twelve putative osmotic and ABA-related genes were impaired in RcDREB2B-overexpressing plants.

CONCLUSIONS

The results provide comprehensive information regarding the gene structure, phylogenic, and distribution of the rose AP2/ERF family and bring insight into the complex transcriptional gene regulation of RcAP2/ERF. Findings in this study would also contribute to further understanding of the RcDREB2B gene in rose.

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 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 identification and expression analysis of the superoxide dismutase (SOD) gene family in Salvia miltiorrhiza

Adverse environmental conditions, such as salinity, cold, drought, heavy metals, and pathogens affect the yield and quality of Salvia miltiorrhiza, a well-known medicinal plant used for the treatment of cardiovascular and cerebrovascular diseases. Superoxide dismutase (SOD), a key enzyme of antioxidant system in plants, plays a vital role in protecting plants against various biotic and abiotic stresses via scavenging the reactive oxygen species produced by organisms. However, little is known about the SOD gene family in S. miltiorrhiza. In this study, eight SOD genes, including three Cu/Zn-SODs, two Fe-SODs and three Mn-SODs, were identified in the S. miltiorrhiza genome. Their gene structures, promoters, protein features, phylogenetic relationships, and expression profiles were comprehensively investigated. Gene structure analysis implied that most SmSODs have different introns/exons distrbution patterns. Many cis-elements related to different stress responses or plant hormones were found in the promoter of each SmSOD. Expression profile analysis indicated that SmSODs exhibited diverse responses to cold, salt, drought, heavy metal, and plant hormones. Additionally, 31 types of TFs regulating SmSODs were predicted and analyzed. These findings provided valuable information for further researches on the functions and applications of SmSODs in S. miltiorrhiza growth and adaptation to stress.

The expanding roles of APETALA2/Ethylene Responsive Factors and their potential applications in crop improvement

Understanding the molecular basis of the gene-regulatory networks underlying agronomic traits or plant responses to abiotic/biotic stresses is very important for crop improvement. In this context, transcription factors, which either singularly or in conjugation directly control the expression of many target genes, are suitable candidates for improving agronomic traits via genetic engineering. In this regard, members of one of the largest class of plant-specific APETALA2/Ethylene Response Factor (AP2/ERF) superfamily, which is implicated in various aspects of development and plant stress adaptation responses, are considered high-value targets for crop improvement. Besides their long-known regulatory roles in mediating plant responses to abiotic stresses such as drought and submergence, the novel roles of AP2/ERFs during fruit ripening or secondary metabolites production have also recently emerged. The astounding functional plasticity of AP2/ERF members is considered to be achieved by their interplay with other regulatory networks and signalling pathways. In this review, we have integrated the recently accumulated evidence from functional genomics studies and described their newly emerged functions in plants. The key structural features of AP2/ERF proteins and the modes of their action are briefly summarized. The importance of AP2/ERFs in plant development and stress responses and a summary of the event of their successful applications in crop improvement programs are also provided. Altogether, we envisage that the synthesized information presented in this review will be useful to design effective strategies for improving agronomic traits in crop plants.

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

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.

The biosynthesis of phenolic acids is positively regulated by the JA-responsive transcription factor ERF115 in Salvia miltiorrhiza

Phenolic acids are important secondary metabolites produced in the Chinese medicinal plant Salvia miltiorrhiza, but little is known about the transcription factors involved in the regulation of tanshinone and phenolic acid biosynthesis. Here, a novel AP2/ERF transcription factor SmERF115 was isolated and functionally characterized. SmERF115 was most responsive to methyl jasmonate (MeJA) treatment and was localized in the nucleus. The phenolic acid production was increased in SmERF115-overexpressing hairy roots, but with a decrease in tanshinone content. In contrast, silencing of SmERF115 reduced the phenolic acid level, but increased tanshinone content. The expression of the key biosynthetic gene SmRAS1 was up-regulated in SmERF115 overexpression lines but was down-regulated in SmERF115-RNAi lines. Yeast one-hybrid (Y1H) assay and EMSA showed that SmERF115 directly binds to the promoter of SmRAS1, while dual-luciferase assays showed that SmERF115 could activate expression of SmRAS1 in vivo. Furthermore, global transcriptomic analysis by RNA sequencing revealed that expression of other genes such as PAL3, 4CL5, TAT3, and RAS4 was also increased in the overexpression line, implying that they were potentially involved in the SmERF115-mediated pathway. Our data show that SmERF115 is a positive regulator of phenolic acid biosynthesis, and may be a potential target for further metabolic engineering of phenolic acid biosynthesis in S. miltiorrhiza.

Systematic analysis of DEMETER-like DNA glycosylase genes shows lineage-specific Smi-miR7972 involved in SmDML1 regulation in Salvia miltiorrhiza

DEMETER-like DNA glycosylases (DMLs) initiate the base excision repair-dependent DNA demethylation to regulate a wide range of biological processes in plants. Six putative SmDML genes, termed SmDML1-SmDML6, were identified from the genome of S. miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies. Integrated analysis of gene structures, sequence features, conserved domains and motifs, phylogenetic analysis and differential expression showed the conservation and divergence of SmDMLs. SmDML1, SmDML2 and SmDML4 were significantly down-regulated by the treatment of 5Aza-dC, a general DNA methylation inhibitor, suggesting involvement of SmDMLs in genome DNA methylation change. SmDML1 was predicted and experimentally validated to be target of Smi-miR7972. Computational analysis of forty whole genome sequences and almost all of RNA-seq data from Lamiids revealed that MIR7972s were only distributed in some plants of the three orders, including Lamiales, Solanales and Boraginales, and the number of MIR7972 genes varied among species. It suggests that MIR7972 genes underwent expansion and loss during the evolution of some Lamiids species. Phylogenetic analysis of MIR7972s showed closer evolutionary relationships between MIR7972s in Boraginales and Solanales in comparison with Lamiales. These results provide a valuable resource for elucidating DNA demethylation mechanism 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.

Transcriptomic analyses reveal biosynthetic genes related to rosmarinic acid in Dracocephalum tanguticum

Dracocephalum tanguticum Maxim, a Lamiaceae species endemic to the Qinghai-Tibetan Plateau and adjacent regions, is an important ornamental, medicinal and aromatic herb. In this study, a comprehensive transcriptome of 18 libraries from six organs namely, roots, stems, leaves, sepals, flowers and seeds of D. tanguticum were generated. More than 100 Gb of sequence data were obtained and assembled de novo into 187,447 transcripts, including 151,463 unigenes, among which the six organs shared 17.7% (26,841). In addition, all unigenes were assigned to 362 pathways, in which 'biosynthesis of secondary metabolites' is the second enriched pathway. Furthermore, rosmarinic acid (RA) is one of the multifunctional phenolic bioactive compounds produced in some Lamiaceae species. The six organs of D. tanguticum were confirmed to produce RA. A total of 22 predicted biosynthetic genes related to RA from the transcriptome were further isolated. Two of these genes were identified as candidates by evaluating the correlation coefficient between the RA contents and the expression of the predicted biosynthetic genes in the six organs. The new sequence information will improve the knowledge of D. tanguticum, as well as provide a reference tool for future studies of biosynthetic genes related to RA in this species.

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.

Summary: Genome-wide analysis identified 25 ARF gene members (seven transcriptional activators and 18 repressors) in S. miltiorrhiza. The gene structures, functional domains, miRNA targets and expression patterns were analyzed in detail.