Molecular cloning and expression analysis of WRKY transcription factor genes in Salvia miltiorrhiza

Genome-wide profiling of WRKY genes involved in flavonoid biosynthesis in Erigeron breviscapus

The transcription factors of WRKY genes play essential roles in plant growth, stress responses, and metabolite biosynthesis. Erigeron breviscapus, a traditional Chinese herb, is abundant in flavonoids and has been used for centuries to treat cardiovascular and cerebrovascular diseases. However, the WRKY transcription factors that regulate flavonoid biosynthesis in E. breviscapus remain unknown. In this study, a total of 75 EbWRKY transcription factors were predicted through comprehensive genome-wide characterization of E. breviscapus and the chromosomal localization of each EbWRKY gene was investigated. RNA sequencing revealed transient responses of 74 predicted EbWRKY genes to exogenous abscisic acid (ABA), salicylic acid (SA), and gibberellin 3 (GA3) after 4 h of treatment. In contrast, the expression of key structural genes involved in flavonoid biosynthesis increased after 4 h in GA3 treatment. However, the content of flavonoid metabolites in leaves significantly increased at 12 h. The qRT-PCR results showed that the expression patterns of EbWRKY11, EbWRKY30, EbWRKY31, EbWRKY36, and EbWRKY44 transcription factors exhibited a high degree of similarity to the 11 structural genes involved in flavonoid biosynthesis. Protein-DNA interactions were performed between the key genes involved in scutellarin biosynthesis and candidate WRKYs. The result showed that F7GAT interacts with EbWRKY11, EbWRKY36, and EbWRKY44, while EbF6H has a self-activation function. This study provides comprehensive information on the regulatory control network of flavonoid accumulation mechanisms, offering valuable insights for breeding E. breviscapus varieties with enhanced scutellarin content.

Genome-Wide Identification and Expression Analysis of Sucrose Nonfermenting 1-Related Protein Kinase (SnRK) Genes in Salvia miltiorrhiza in Response to Hormone

The SnRK gene family is the chief component of plant stress resistance and metabolism through activating the phosphorylation of downstream proteins. S. miltiorrhiza is widely used for the treatment of cardiovascular diseases in Asian countries. However, information about the SnRK gene family of S. miltiorrhiza is not clear. The aim of this study is to comprehensively analyze the SnRK gene family of S. miltiorrhiza and its response to phytohormone. Here, 33 SmSnRK genes were identified and divided into three subfamilies (SmSnRK1, SmSnRK2 and SmSnRK3) according to phylogenetic analysis and domain. SmSnRK genes within same subgroup shared similar protein motif composition and were unevenly distributed on eight chromosomes of S. miltiorrhiza. Cis-acting element analysis showed that the promoter of SmSnRK genes was enriched with ABRE motifs. Expression pattern analysis revealed that SmSnRK genes were preferentially expressed in leaves and roots. Most SmSnRK genes were induced by ABA and MeJA treatment. Correlation analysis showed that SmSnRK3.15 and SmSnRK3.18 might positively regulate tanshinone biosynthesis; SmSnRK3.10 and SmSnRK3.12 might positively regulate salvianolic acid biosynthesis. RNAi-based silencing of SmSnRK2.6 down-regulated the biosynthesis of tanshinones and biosynthetic genes expression. An in vitro phosphorylation assay verified that SmSnRK2.2 interacted with and phosphorylated SmAREB1. These findings will provide a valuable basis for the functional characterization of SmSnRK genes and quality improvement of S. miltiorrhiza.

Phenolic Acids and Flavonoids Play Important Roles in Flower Bud Differentiation in Mikania micrantha: Transcriptomics and Metabolomics

Mikania micrantha is a highly invasive vine, and its ability to sexually reproduce is a major obstacle to its eradication. The long-distance dissemination of M. micrantha depends on the distribution of seeds; therefore, inhibiting M. micrantha flowering and seed production is an effective control strategy. The number of blooms of M. micrantha differs at different altitudes (200, 900, and 1300 m). In this study, we used a combination of metabolomics and transcriptomics methods to study the patterns of metabolite accumulation in the flower buds of M. micrantha. Using LC-MS/MS, 658 metabolites were found in the flower buds of M. micrantha at three different altitudes (200, 900, and 1300 m). Flavonoids and phenolic acids were found to be the main differential metabolites, and their concentrations were lower at 900 m than at 200 m and 1300 m, with the concentrations of benzoic acid, ferulic acid, and caffeic acid being the lowest. The biosynthesis pathways for flavonoids and phenolic compounds were significantly enriched for differentially expressed genes (DEGs), according to the results of transcriptome analysis. The production of flavonoid and phenolic acids was strongly linked with the expressions of phenylalanine ammonia-lyase (PAL), caffeoyl-CoA O-methyltransferase (COMT), and 4-coumarate-CoA ligase (4CL), according to the results of the combined transcriptome and metabolome analysis. These genes' roles in the regulation of distinct phenolic acids and flavonoids during M. micrantha bud differentiation are still unknown. This study adds to our understanding of how phenolic acids and flavonoids are regulated in M. micrantha flower buds at various altitudes and identifies regulatory networks that may be involved in this phenomenon, offering a new approach for the prevention and management of M. micrantha.

Genome-Wide Identification and Characterization of the WRKY Gene Family in Cucurbita maxima

In higher plants, WRKY transcription factors are broadly involved in a variety of life activities and play an important role in both biotic and abiotic stress responses. However, little is known about the functions of WRKY genes in the popular species, such as Cucurbita maxima (pumpkin), which is planted worldwide. In the present study, 102 CmWRKY genes were identified in the C. maxima genome. Chromosome location, multiple sequence alignment, phylogenetic analysis, and synteny analysis of the CmWRKYs were performed. Notably, we found that silencing CmWRKY22 promoted cucumber mosaic virus (CMV) infection, whereas overexpression of CmWRKY22 inhibited the CMV infection. Subsequently, an electrophoretic mobility shift assay (EMSA) confirmed that CmWRKY22 was able to bind to the W-box at the promoter of CmPR1b, which is a responsive gene of the salicylic acid (SA) signaling pathway. In summary, this study has provided a foundation for the antiviral functions of WRKY transcription factors in C. maxima.

Overexpression of a WRKY transcription factor McWRKY57-like from Mentha canadensis L. enhances drought tolerance in transgenic Arabidopsis

BACKGROUND

Drought has become a major environmental problem affecting crop production. Members of the WRKY family play important roles in plant development and stress responses. However, their roles in mint have been barely explored.

RESULTS

In this study, we isolated a drought-inducible gene McWRKY57-like from mint and investigated its function. The gene encodes a group IIc WRKY transcription factor, McWRKY57-like, which is a nuclear protein with a highly conserved WRKY domain and a C2H2 zinc-finger structure, and has transcription factor activity. Its expression levels were examined in different tissues of mint and under the treatment of mannitol, NaCl, abscisic acid, and methyl jasmonate. We found that McWRKY57-like overexpression in Arabidopsis significantly increased drought tolerance. Further studies showed that under drought stress, McWRKY57-like-overexpressing plants had higher chlorophyll, soluble sugar, soluble protein, and proline contents but lower water loss rate and malondialdehyde content than wild-type plants. Moreover, the activities of antioxidant enzymes catalase, superoxide dismutase, and peroxidase were enhanced in McWRKY57-like transgenic plants. Furthermore, qRT-PCR analysis revealed that the drought-related genes AtRD29A, AtRD29B, AtRD20, AtRAB18, AtCOR15A, AtCOR15B, AtKIN2, and AtDREB1A were upregulated in McWRKY57-like transgenic plants than in wild-type Arabidopsis under simulated drought conditions.

CONCLUSION

These data demonstrated that McWRKY57-like conferred drought tolerance in transgenic Arabidopsis by regulating plant growth, osmolyte accumulation and antioxidant enzyme activities, and the expression of stress-related genes. The study indicates that McWRKY57-like plays a positive role in drought response in plants.

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.

WRKY transcription factors: evolution, regulation, and functional diversity in plants

The recent advancements in sequencing technologies and informatic tools promoted a paradigm shift to decipher the hidden biological mysteries and transformed the biological issues into digital data to express both qualitative and quantitative forms. The transcriptomic approach, in particular, has added new dimensions to the versatile essence of plant genomics through the large and deep transcripts generated in the process. This has enabled the mining of super families from the sequenced plants, both model and non-model, understanding their ancestry, diversity, and evolution. The elucidation of the crystal structure of the WRKY proteins and recent advancement in computational prediction through homology modeling and molecular dynamic simulation has provided an insight into the DNA-protein complex formation, stability, and interaction, thereby giving a new dimension in understanding the WRKY regulation. The present review summarizes the functional aspects of the high volume of sequence data of WRKY transcription factors studied from different species, till date. The review focuses on the dynamics of structural classification and lineage in light of the recent information. Additionally, a comparative analysis approach was incorporated to understand the functions of the identified WRKY transcription factors subjected to abiotic (heat, cold, salinity, senescence, dark, wounding, UV, and carbon starvation) stresses as revealed through various sets of studies on different plant species. The review will be instrumental in understanding the events of evolution and the importance of WRKY TFs under the threat of climate change, considering the new scientific evidences to propose a fresh perspective.

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.

Genome-wide identification and expression analysis response to GA3 stresses of WRKY gene family in seed hemp (Cannabis sativa L)

WRKY transcription factor is one of the largest transcription factor families in higher plants. However, the investigations of the WRKY gene family have not yet been reported in seed hemp. In the present study, we identified 39 CasWRKYs at the genome-wide level and analyzed phylogenetic relationship, chromosome location, cis-acting elements, gene structure, conserved motif, and expression pattern. Based on the gene structure and phylogenetic analyses, CasWRKY proteins were divided into 3 groups and 7 subgroups. The gene duplication investigation revealed that 6 and 5 pairs of CasWRKY genes underwent tandem and segmental duplication events, respectively. These events may contribute to the diversity and expansion of the CasWRKY gene family. The regulatory elements in the promoter regions of CasWRKYs contained diverse cis-regulatory elements, among which P-box cis-regulatory elements showed high frequency, indicating that CasWRKYs can respond to the regulation of gibberellin. The expression profiles derived from RNA-seq and qRT-PCR showed that 13 CasWRKY genes could respond to GA₃ stress and affect fiber development, as well as play significant roles in stem growth and development. This study will serve as molecular basis and practical reference for further exploring the genetic evolution and biological function of CasWRKY genes in seed hemp.

Genome-Wide Identification and Characterization of the WRKY Gene Family in Scutellaria baicalensis Georgi under Diverse Abiotic Stress

The WRKY gene family is an important inducible regulatory factor in plants, which has been extensively studied in many model plants. It has progressively become the focus of investigation for the secondary metabolites of medicinal plants. Currently, there is no systematic analysis of the WRKY gene family in Scutellaria baicalensis Georgi. For this study, a systematic and comprehensive bioinformatics analysis of the WRKY gene family was conducted based on the genomic data of S. baicalensis. A total of 77 WRKY members were identified and 75 were mapped onto nine chromosomes, respectively. Their encoded WRKY proteins could be classified into three subfamilies: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III, based on the characteristics of the amino acid sequences of the WRKY domain and genetic structure. Syntenic analysis revealed that there were 35 pairs of repetitive fragments. Furthermore, the transcriptome data of roots, stems, leaves, and flowers showed that the spatial expression profiles of WRKYs were different. qRT-PCR analysis revealed that 11 stress-related WRKYs exhibited specific expression patterns under diverse treatments. In addition, sub cellular localization analysis indicated that SbWRKY26 and SbWRKY41 were localized in nucleus. This study is the first to report the identification and characterization of the WRKY gene family in S. baicalensis, which is valuable for the further exploration of the biological function of SbWRKYs. It also provides valuable bioinformatics data for S. baicalensis and provides a reference for assessing the medicinal properties of the genus.

Comparative transcriptome analyses of different Salvia miltiorrhiza varieties during the accumulation of tanshinones

Salvia miltiorrhiza (Labiatae) is an important medicinal plant in traditional Chinese medicine. Tanshinones are one of the main active components of S. miltiorrhiza. It has been found that the intraspecific variation of S. miltiorrhiza is relatively large and the content of tanshinones in its roots of different varieties is also relatively different. To investigate the molecular mechanisms that responsible for the differences among these varieties, the tanshinones content was determined and comparative transcriptomics analysis was carried out during the tanshinones accumulation stage. A total of 52,216 unigenes were obtained from the transcriptome by RNA sequencing among which 23,369 genes were differentially expressed among different varieties, and 2,016 genes including 18 diterpenoid biosynthesis-related genes were differentially expressed during the tanshinones accumulation stage. Functional categorization of the differentially expressed genes (DEGs) among these varieties revealed that the pathway related to photosynthesis, oxidative phosphorylation, secondary metabolite biosynthesis, diterpenoid biosynthesis, terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis are the most differentially regulated processes in these varieties. The six tanshinone components in these varieties showed different dynamic changes in tanshinone accumulation stage. In addition, combined with the analysis of the dynamic changes, 277 DEGs (including one dehydrogenase, three CYP450 and 24 transcription factors belonging to 12 transcription factor families) related to the accumulation of tanshinones components were obtained. Furthermore, the KEGG pathway enrichment analysis of these 277 DEGs suggested that there might be an interconnection between the primary metabolic processes, signaling processes and the accumulation of tanshinones components. This study expands the vision of intraspecific variation and gene regulation mechanism of secondary metabolite biosynthesis pathways in medicinal plants from the “omics” perspective.

Identification and Application of BhAPRR2 Controlling Peel Colour in Wax Gourd (Benincasa hispida)

Peel color is an important factor affecting commodity quality in vegetables; however, the genes controlling this trait remain unclear in wax gourd. Here, we used two F₂ genetic segregation populations to explore the inheritance patterns and to clone the genes associated with green and white skin in wax gourd. The F₂ and BC₁ trait segregation ratios were 3:1 and 1:1, respectively, and the trait was controlled by nuclear genes. Bulked segregant analysis of both F₂ plants revealed peaks on Chr5 exceeding the confidence interval. Additionally, 6,244 F₂ plants were used to compress the candidate interval into a region of 179 Kb; one candidate gene, Bch05G003950 (BhAPRR2), encoding two-component response regulator-like protein Arabidopsis pseudo-response regulator2 (APRR2), which is involved in the regulation of peel color, was present in this interval. Two bases (GA) present in the coding sequence of BhAPRR2 in green-skinned wax gourd were absent from white-skinned wax gourd. The latter contained a frameshift mutation, a premature stop codon, and lacked 335 residues required for the protein functional region. The chlorophyll content and BhAPRR2 expression were significantly higher in green-skinned than in white-skinned wax gourd. Thus, BhAPRR2 may regulate the peel color of wax gourd. This study provides a theoretical foundation for further studies of the mechanism of gene regulation for the fruit peel color of wax gourd.

WRKY transcription factors and plant defense responses: latest discoveries and future prospects

WRKY transcription factors are among the largest families of transcriptional regulators. In this review, their pivotal role in modulating various signal transduction pathways during biotic and abiotic stresses is discussed. Transcription factors (TFs) are important constituents of plant signaling pathways that define plant responses against biotic and abiotic stimuli besides playing a role in response to internal signals which coordinate different interacting partners during developmental processes. WRKY TFs, deriving their nomenclature from their signature DNA-binding sequence, represent one of the largest families of transcriptional regulators found exclusively in plants. By modulating different signal transduction pathways, these TFs contribute to various plant processes including nutrient deprivation, embryogenesis, seed and trichome development, senescence as well as other developmental and hormone-regulated processes. A growing body of research suggests transcriptional regulation of WRKY TFs in adapting plant to a variety of stressed environments. WRKY TFs can regulate diverse biological functions from receptors for pathogen triggered immunity, modulator of chromatin for specific interaction and signal transfer through a complicated network of genes. Latest discoveries illustrate the interaction of WRKY proteins with other TFs to form an integral part of signaling webs that regulate several seemingly disparate processes and defense-related genes, thus establishing their significant contributions to plant immune response. The present review starts with a brief description on the structural characteristics of WRKY TFs followed by the sections that present recent evidence on their roles in diverse biological processes in plants. We provide a comprehensive overview on regulatory crosstalks involving WRKY TFs during multiple stress responses in plants and future prospects of WRKY TFs as promising molecular diagnostics for enhancing crop improvement.

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.

Anaplasma phagocytophilum AptA enhances the UPS, autophagy, and anti-apoptosis of host cells by PSMG3

Anaplasma phagocytophilum is an obligate intracellular bacterium and a common tick-borne infectious pathogen that can cause human granulocytic anaplasmosis (HGA). Effector proteins play an important role in the pathogenic mechanism of A. phagocytophilum, but the specifics of the disease mechanism are unclear. We studied the effector protein AptA (A. phagocytophilum toxin A) using yeast two hybrid assays to screen its interacting protein proteasome assembly chaperone 3 (PSMG3, PAC3), and identified new mechanisms for the pathogenicity of A. phagocytophilum in HEK293T cells. After AptA enters the host cell, it interacts with PSMG3 to enhance the activity of the proteasome, causing ubiquitination and autophagy in the host cell and thereby increasing cross-talk between the ubiquitination-proteasome system (UPS) and autophagy. AptA also reduces the apoptotic efficiency of the host cells. These results offer new clues as to the pathogenic mechanism of A. phagocytophilum and support the hypothesis that AptA interacts with host PSMG3.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Identification and expression analysis of the WRKY gene family in Isatis indigotica

The WRKY proteins, which represent one of the largest families of transcriptional regulators in plants, play pivotal roles in regulating multiple processes of growth and development, particularly in diverse stress responses. Isatis indigotica is widely used in Traditional Chinese Medicine and is famous for its use as a dye for the color indigo. However, reports of the WRKY gene family in I. indigotica are limited. In this study, 64 IiWRKY genes encoding proteins with the complete WRKY domain were identified from genome of I. indigotica. Based on their structure and phylogenetic relationships of this gene family in I. indigotica, the IiWRKY genes were classified into three groups: Group I (n = 13), Group II (n = 35) and Group III (n = 16). Sequence alignment revealed that IiWRKY proteins harbored two variants, WRKYRQK and WRKYGKK, of the highly conserved WRKYGQK motif. The number of exons in IiWRKY genes varied from two to 14, with most of IiWRKY genes containing three exons. Investigation of gene duplication demonstrated that 10 and 14 IiWRKY genes were incorporated in tandem and segmental duplication events, respectively. Finally, the expression profiles derived from transcriptome data and quantitative real-time PCR analysis showed distinct expression patterns of these IiWRKY gene in five different organs or in response to four abiotic stresses. Taken together, our results will contribute to functional analysis of IiWRKY genes, and also provide a basis for further clarification of the molecular mechanism of stress responses in this important herb.

Functional analysis of CgWRKY57 from Cymbidium goeringii in ABA response

Background

The orchid is one of the top ten Chinese flowers and has high ornamental value and elegant color. However, orchids are vulnerable to abiotic stresses during their growth and development, and the molecular mechanism of the abiotic stress response in orchids is unclear. WRKY proteins belong to a transcription factor family that plays important roles in biotic stress, abiotic stress, growth and development in plants, but little is known about the WRKY family in Cymbidium goeringii.

Methods

The specific fragment of the CgWRKY57 gene of C. goeringii was analyzed by bioinformatics. The expression of the CgWRKY57 gene of C. goeringii under 4 °C, 42 °C water and ABA stress as well as different tissues was detected by real-time fluorescence quantitative PCR. CgWRKY57 gene was overexpressed in wild type Arabidopsis thaliana by inflorescence infection method, and the function of transgenic lines under ABA stress was analyzed.

Results

CgWRKY57 was cloned from C. goeringii and found to encode 303 amino acids. The CgWRKY57 protein is an acidic, nonsecreted hydrophilic protein without a signal peptide or transmembrane domain. The CgWRKY57 protein is located to the nucleus and may function intracellularly according to its predicted subcellular localization. A domain analysis and homology comparison showed that the CgWRKY57 protein has a "WRKYGQK" domain and belongs to Group III of the WRKY family, and a phylogenetic analysis demonstrated that CgWRKY57 is closely related to OsWRKY47. CgWRKY57 was expressed in the roots, stems, leaves and floral organs of C. goeringii, and its expression level was highest in the roots according to real-time qPCR analysis. There were significant differences in CgWRKY57 expression under 4 °C, 42 °C ABA and water stress treatments, and its expression changed greatly under ABA stress. The expression of CgWRKY57 in transgenic plants was significantly higher than that in wild type plants under ABA stress, and the root length and germination rate were reduced in transgenic plants compared to wild type plants.

Conclusions

These results indicate that CgWRKY57 overexpression is responsive to ABA stress, and they provide a foundation for future analyses of the biological functions of the WRKY family in C. goeringii.

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.

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 identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

In-silico analysis of cucumber (Cucumis sativus L.) Genome for WRKY transcription factors and cis-acting elements

WRKY genes, comprises one among a large clan of transcription factor (TFs) genes in the plant kingdom, playing a fundamental role in the vegetative and reproductive growth, development and stress responses of a plant. In spite of several studies on cucumber (Cucumis sativus L.), WRKY genes and their interaction with stress response is limited. The present study, on the whole genome of cucumber was analyzed for WRKY genes which recognized 62 CsWRKY genes associated with the proteins obtained from lineages of supplementary plants. The physicochemical properties reveal the CsWRKY gene is ser-rich TF (6.70-18.40 %). The chromosomal distribution showed that all putative CsWRKY genes were distributed in seven chromosomes, enriched on chromosome 3 and 6 and least on chromosome 5. Based on phylogenetic analysis, along with motif determination and gene structure analysis, CsWRKYs are categorized as a Group I, II and III. The Group II further subdivided as Groups IIa-e. In the present study, it was observed that Group II WRKY-TFs was the largest group containing 43 WRKY genes containing a single WD (WRKY domain - WRKYGQK/WRKYGKK) and C₂H₂ type zinc finger structure (C-X_4-5-C-X₂₃-H-X₁-H). The data also revealed that chromosome 3 and 5 contained all the three major groups and chromosome 6 contained I and II WRKY genes with uneven distribution. STRING analysis of selected CsWRKY proteins expressed in response to abiotic stress interacts with the CsMAPK proteins. Analysis of cis-acting elements and results suggest that CSWRKY genes play important role in response to biotic and abiotic stress. Response also predicted the candidate gene expression in cucumber during its development under different cellular condition.

Transcriptomic analysis reveals potential genes involved in tanshinone biosynthesis in Salvia miltiorrhiza

Tanshinones are important bioactive components in Salvia miltiorrhiza and mainly accumulate in the periderms of mature roots. Tanshinone biosynthesis is a complicated process, and little is known about the third stage of the pathway. To investigate potential genes that are responsible for tanshinone biosynthesis, we conducted transcriptome profiling analysis of two S. miltiorrhiza cultivars. Differential expression analysis provided 2,149 differentially expressed genes (DEGs) for further analysis. GO and KEGG analysis showed that the DEGs were mainly associated with the biosynthesis of secondary metabolites. Weighted gene coexpression network analysis (WGCNA) was further performed to identify a “cyan” module associated with tanshinone biosynthesis. In this module, 25 cytochromes P450 (CYPs), three 2-oxoglutarate-dependent dioxygenases (2OGDs), one short-chain alcohol dehydrogenases (SDRs) and eight transcription factors were found to be likely involved in tanshinone biosynthesis. Among these CYPs, 14 CYPs have been reported previously, and 11 CYPs were identified in this study. Expression analysis showed that four newly identified CYPs were upregulated upon application of MeJA, suggesting their possible roles in tanshinone biosynthesis. Overall, this study not only identified candidate genes involved in tanshinone biosynthesis but also provided a basis for characterization of genes involved in important active ingredients of other traditional Chinese medicinal plants.

Analysis of the laccase gene family and miR397-/miR408-mediated posttranscriptional regulation in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most commonly used traditional Chinese medicine materials. It contains important bioactive phenolic compounds, such as salvianolic acids, flavonoids and anthocyanins. Elucidation of phenolic compound biosynthesis and its regulatory mechanism is of great significance for S. miltiorrhiza quality improvement. Laccases (LACs) are multicopper-containing enzymes potentially involved in the polymerization of phenolic compounds. So far, little has been known about LAC genes in S. miltiorrhiza. Through systematic investigation of the whole genome sequence and transcriptomes of S. miltiorrhiza, we identified 65 full-length SmLAC genes (SmLAC1–SmLAC65). Phylogenetic analysis showed that 62 of the identified SmLACs clustered with LACs from Arabidopsis and Populus trichocarpa in seven clades (C1–C7), whereas the other three fell into one S. miltiorrhiza-specific clade (C8). All of the deduced SmLAC proteins contain four conserved signature sequences and three typical Cu-oxidase domains, and gene structures of most LACs from S. miltiorrhiza, Arabidopsis and P. trichocarpa were highly conserved, however SmLACs encoding C8 proteins showed distinct intron-exon structures. It suggests the conservation and diversity of plant LACs in gene structures. The majority of SmLACs exhibited tissue-specific expression patterns, indicates manifold functions of SmLACs played in S. miltiorrhiza. Analysis of high-throughput small RNA sequences and degradome data and experimental validation using the 5′ RACE method showed that 23 SmLACs were targets of Smi-miR397. Among them, three were also targeted by Smi-miR408. It suggests the significance of miR397 and miR408 in posttranscriptional regulation of SmLAC genes. Our results provide a foundation for further demonstrating the functions of SmLACs in the production of bioactive phenolic compounds in S. miltiorrhiza.

Genome-Wide Comprehensive Analysis of the SABATH Gene Family in Arabidopsis and Rice

Low molecular weight metabolites are important plant hormones and signaling molecules, and play an important part among the processes of plant development. Their activities may also be affected by the chemical modifications of methylation performed by SABATH. In this study, a total of 24 and 21 SABATH genes in Arabidopsis and rice, respectively, were identified and taken a comprehensive study. Phylogenetic analysis showed that AtSABATH and OsSABATH genes could be classified into 4 major groups and 6 subgroups. Gene expansion analysis showed that the main expansion mechanism of SABATH gene family in Arabidopsis and rice was tandem duplication and segmental duplication. The ratios of nonsynonymous (Ka) and synonymous (Ks) substitution rates of 12 pairs paralogous of AtSABATH and OsSABATH genes indicated that the SABATH gene family in Arabidopsis and rice had gone through purifying selection. Positive selection analysis with site models and branch-site models revealed that AtSABATH and OsSABATH genes had undergone selective pressure for adaptive evolution. Motif analysis showed that certain motifs only existed in specific subgroups or species, which indicated that the SABATH proteins of Arabidopsis and rice appear divergence in different species and subgroups. Functional divergence analysis also suggested that the AtSABATH and OsSABATH subgroup genes had functional differences, and the positive selection sites which contributed to functional divergence among subgroups were detected. These results provide insights into functional conservation and diversification of SABATH gene family, and are useful information for further elucidating SABATH gene family functions.

Signalling molecules responsive to ozone-induced oxidative stress in Salvia officinalis

Tropospheric ozone (O₃) is the most important gaseous pollutant and induces a mass of negative impacts on vegetation at functional and genic levels. The aim of the present study was to investigate the role of reactive oxygen species and signalling molecules in sage plants exposed to O₃ (200 ppb, 5 h). Ozone exposure induced only a transient oxidative burst, as confirmed by the rapid peak of anion superoxide during the first hours of exposure (+16% compared to controls). The spontaneous reaction of O₃ with membrane fatty acids stimulates peroxidative processes, as demonstrated by the rise of thiobarbituric acid reactive substances concentration starting after 1 h of exposure (+25%). The formation of lipid-based signalling molecules (e.g. jasmonic acid) may be regarded as a sort of O₃-perception. The concomitant accumulation of salicylic acid suggests that sage responds early to O₃ by inducing cellular antioxidants mechanisms in order to minimize O₃-oxidative burst. The transient increase of abscisic acid (+25% at the end of the treatment) twinned with the maximal ethylene emission (about two-fold higher than controls) could be interpreted as a first attempt by plants to regulate the signalling responses induced by O₃. In order to investigate the involvement of transcription factors in managing oxidative protection, BLASTX analysis against the Salvia miltiorrhiza sequence genome was carried out using Arabidopsis thaliana WRKY sequences as queries. Six gene sequences were identified for sage WRKYs and their relative gene expression analyses were characterized. WRKY4, WRKY5, WRKY11 and WRKY46 were up-regulated by O₃ at 2 and 5 h of exposure and they showed similarity with AtWRKY48, AtWRKY22 and AtWRKY53 in A. thaliana. These results suggest that WRKYs could play a pivotal role in the signalling mechanisms during the responses of plants to O₃.

Coactivation of MEP-biosynthetic genes and accumulation of abietane diterpenes in Salvia sclarea by heterologous expression of WRKY and MYC2 transcription factors

Plant abietane diterpenoids (e.g. aethiopinone, 1- oxoaethiopinone, salvipisone and ferruginol), synthesized in the roots of several Salvia spp, have antibacterial, antifungal, sedative and anti-proliferative properties. Recently we have reported that content of these compounds in S. sclarea hairy roots is strongly depending on transcriptional regulation of genes belonging to the plastidial MEP-dependent terpenoid pathway, from which they mostly derive. To boost the synthesis of this interesting class of compounds, heterologous AtWRKY18, AtWRKY40, and AtMYC2 TFs were overexpressed in S. sclarea hairy roots and proved to regulate in a coordinated manner the expression of several genes encoding enzymes of the MEP-dependent pathway, especially DXS, DXR, GGPPS and CPPS. The content of total abietane diterpenes was enhanced in all overexpressing lines, although in a variable manner due to a negative pleiotropic effect on HR growth. Interestingly, in the best performing HR lines overexpressing the AtWRKY40 TF induced a significant 4-fold increase in the final yield of aethiopinone, for which we have reported an interesting anti-proliferative activity against resistant melanoma cells. The present results are also informative and instrumental to enhance the synthesis of abietane diterpenes derived from the plastidial MEP-derived terpenoid pathway in other Salvia species.

Transcriptional Profiles of SmWRKY Family Genes and Their Putative Roles in the Biosynthesis of Tanshinone and Phenolic Acids in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a Chinese traditional herb for treating cardiovascular and cerebrovascular diseases, and tanshinones and phenolic acids are the dominated medicinal and secondary metabolism constituents of this plant. WRKY transcription factors (TFs) can function as regulators of secondary metabolites biosynthesis in many plants. However, studies on the WRKY that regulate tanshinones and phenolics biosynthesis are limited. In this study, 69 SmWRKYs were identified in the transcriptome database of S. miltiorrhiza, and phylogenetic analysis indicated that some SmWRKYs had closer genetic relationships with other plant WRKYs, which were involved in secondary metabolism. Hairy roots of S. miltiorrhiza were treated by methyl jasmonate (MeJA) to detect the dynamic change trend of SmWRKY, biosynthetic genes, and medicinal ingredients accumulation. Base on those date, a correlation analysis using Pearson’s correlation coefficient was performed to construct gene-to-metabolite network and identify 9 SmWRKYs (SmWRKY1, 7, 19, 29, 45, 52, 56, 58, and 68), which were most likely to be involved in tanshinones and phenolic acids biosynthesis. Taken together, this study has provided a significant resource that could be used for further research on SmWRKY in S. miltiorrhiza and especially could be used as a cue for further investigating SmWRKY functions in secondary metabolite accumulation.

Identification and characterization of the cytosine-5 DNA methyltransferase gene family in Salvia miltiorrhiza

Cytosine DNA methylation is highly conserved epigenetic modification involved in a wide range of biological processes in eukaryotes. It was established and maintained by cytosine-5 DNA methyltransferases (C5-MTases) in plants. Through genome-wide identification, eight putative SmC5-MTase genes were identified from the genome of Salvia miltiorrhiza, a well-known traditional Chinese medicine material and an emerging model medicinal plant. Based on conserved domains and phylogenetic analysis, eight SmC5-MTase genes were divided into four subfamilies, including MET, CMT, DRM and DNMT2. Genome-wide comparative analysis of the C5-MTase gene family in S. miltiorrhiza and Arabidopsis thaliana, including gene structure, sequence features, sequence alignment and conserved motifs, was carried out. The results showed conservation and divergence of the members of each subfamily in plants. The length of SmC5-MTase open reading frames ranges widely from 1,152 (SmDNMT2) to 5,034 bp (SmMET1). The intron number of SmC5-MTases varies between 7 (SmDRM1) and 20 (SmCMT1 and SmCMT2b). These features were similar to their counterparts from Arabidopsis. Sequence alignment and conserved motif analysis showed the existence of highly conserved and subfamily-specific motifs in the C5-MTases analyzed. Differential transcript abundance was detected for SmC5-MTases, implying genome-wide variance of DNA methylation in different organs and tissues. Transcriptome-wide analysis showed that the transcript levels of all SmC5-MTase genes was slightly changed under yeast extract and methyl jasmonate treatments. Six SmC5-MTases, including SmMET1, SmCMT1, SmCMT2a, SmCMT2b, SmCMT3 and SmDRM1, were salicylic acid-responsive, suggesting the involvement of SmC5-MTases in salicylic acid-dependent immunity. These results provide useful information for demonstrating the role of DNA methylation in bioactive compound biosynthesis and Dao-di herb formation in medicinal plants.

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.

Identification of the WRKY gene family and functional analysis of two genes in Caragana intermedia

BACKGROUND

WRKY transcription factors, one of the largest families of transcriptional regulators in plants, play important roles in plant development and various stress responses. The WRKYs of Caragana intermedia are still not well characterized, although many WRKYs have been identified in various plant species.

RESULTS

We identified 53 CiWRKY genes from C. intermedia transcriptome data, 28 of which exhibited complete open reading frames (ORFs). These CiWRKYs were divided into three groups via phylogenetic analysis according to their WRKY domains and zinc finger motifs. Conserved domain analysis showed that the CiWRKY proteins contain a highly conserved WRKYGQK motif and two variant motifs (WRKYGKK and WKKYEEK). The subcellular localization of CiWRKY26 and CiWRKY28-1 indicated that these two proteins localized exclusively to nuclei, supporting their role as transcription factors. The expression patterns of the 28 CiWRKYs with complete ORFs were examined through quantitative real-time PCR (qRT-PCR) in various tissues and under different abiotic stresses (drought, cold, salt, high-pH and abscisic acid (ABA)). The results showed that each CiWRKY responded to at least one stress treatment. Furthermore, overexpression of CiWRKY75-1 and CiWRKY40-4 in Arabidopsis thaliana suppressed the drought stress tolerance of the plants and delayed leaf senescence, respectively.

CONCLUSIONS

Fifty-three CiWRKY genes from the C. intermedia transcriptome were identified and divided into three groups via phylogenetic analysis. The expression patterns of the 28 CiWRKYs under different abiotic stresses suggested that each CiWRKY responded to at least one stress treatment. Overexpression of CiWRKY75-1 and CiWRKY40-4 suppressed the drought stress tolerance of Arabidopsis and delayed leaf senescence, respectively. These results provide a basis for the molecular mechanism through which CiWRKYs mediate stress tolerance.

Genome-Wide Comprehensive Analysis the Molecular Phylogenetic Evaluation and Tissue-Specific Expression of SABATH Gene Family in Salvia miltiorrhiza

The plant SABATH gene family is a group of O-methyltransferases (O-MTs), which belongs to the S-adenosyl-l-methionine-dependent methyltransferases (SAM-MTs). The resulting reaction products of SABATH genes play an important role in various processes of plant development. In this study, a total of 30 SABATH genes were detected in Salvia miltiorrhiza, which is an important medicinal plant, widely used to treat cardiovascular disease. Multiple sequence alignment and phylogenetic analyses showed that SmSABATH genes could be classified into three groups. The ratios of non-synonymous (Ka) and synonymous (Ks) substitution rates of 11 pairs paralogous of SmSABATH genes revealed that the SmSABATH genes had gone through purifying selection. Positive selection analyses using site models and branch-site models indicated that SmSABATH genes had undergone selective pressure for adaptive evolution. Functional divergence analyses suggested that the SmSABATH subgroup genes were divergent in terms of functions and positive selection sites that contributed to a functional divergence among the subgroups that were detected. Tissue-specific expression showed that the SABATH gene family in S. miltiorrhiza was primarily expressed in stems and leaves.

Genome-wide analysis of WRKY transcription factors in wheat (Triticum aestivum L.) and differential expression under water deficit condition

BACKGROUND

WRKY proteins, which comprise one of the largest transcription factor (TF) families in the plant kingdom, play crucial roles in plant development and stress responses. Despite several studies on WRKYs in wheat (Triticum aestivum L.), functional annotation information about wheat WRKYs is limited.

RESULTS

Here, 171 TaWRKY TFs were identified from the whole wheat genome and compared with proteins from 19 other species representing nine major plant lineages. A phylogenetic analysis, coupled with gene structure analysis and motif determination, divided these TaWRKYs into seven subgroups (Group I, IIa-e, and III). Chromosomal location showed that most TaWRKY genes were enriched on four chromosomes, especially on chromosome 3B. In addition, 85 (49.7%) genes were either tandem (5) or segmental duplication (80), which suggested that though tandem duplication has contributed to the expansion of TaWRKY family, segmental duplication probably played a more pivotal role. Analysis of cis-acting elements revealed putative functions of WRKYs in wheat during development as well as under numerous biotic and abiotic stresses. Finally, the expression of TaWRKY genes in flag leaves, glumes, and lemmas under water-deficit condition were analyzed. Results showed that different TaWRKY genes preferentially express in specific tissue during the grain-filling stage.

CONCLUSION

Our results provide a more extensive insight on WRKY gene family in wheat, and also contribute to the screening of more candidate genes for further investigation on function characterization of WRKYs under various stresses.

Characterization of the polyphenol oxidase gene family reveals a novel microRNA involved in posttranscriptional regulation of PPOs in Salvia miltiorrhiza

Salvia miltiorrhiza is a well-known material of traditional Chinese medicine. Understanding the regulatory mechanisms of phenolic acid biosynthesis and metabolism are important for S. miltiorrhiza quality improvement. We report here that S. miltiorrhiza contains 19 polyphenol oxidases (PPOs), forming the largest PPO gene family in plant species to our knowledge. Analysis of gene structures and sequence features revealed the conservation and divergence of SmPPOs. SmPPOs were differentially expressed in plant tissues and eight of them were predominantly expressed in phloem and xylem, indicating that some SmPPOs are functionally redundant, whereas the others are associated with different physiological processes. Expression patterns of eighteen SmPPOs were significantly altered under MeJA treatment, and twelve were yeast extract and Ag⁺-responsive, suggesting the majority of SmPPOs are stress-responsive. Analysis of high-throughput small RNA sequences and degradome data showed that miR1444-mediated regulation of PPOs existing in P. trichocarpa is absent from S. miltiorrhiza. Instead, a subset of SmPPOs was posttranscriptionally regulated by a novel miRNA, termed Smi-miR12112. It indicates the specificity and significance of miRNA-mediated regulation of PPOs. The results shed light on the regulation of SmPPO expression and suggest the complexity of SmPPO-associated phenolic acid biosynthesis and metabolism.

OCF can repress tumor metastasis by inhibiting epithelial-mesenchymal transition involved in PTEN/PI3K/AKT pathway in lung cancer cells

A component formula with definite compositions provides a new approach to treat various diseases. Salvia miltiorrhiza and Panax ginseng are widely used in China because of their antitumor properties. In the previous study, the optimizing component formula (OCF), prepared with salvianolic acids, ginsenosides, and ginseng polysaccharides (5, 10, and 5 mg·L⁻¹, respectively) extracted from S. miltiorrhiza and P. ginseng on the basis of IC₅₀ in lung cancer A549 cells and damage minimization on human bronchial epithelial cells in vitro. Currently, we also have demonstrated the inhibitory effect of OCF on A549 cell migration and invasion in vitro. According to Lewis lung cancer cells (LLC) allograft in C57BL/6 mice and A549 xenograft in nude mice experiment, we found that the anti-tumor and anti-metastasis effects of OCF treatment were related to the inhibition of epithelial–mesenchymal transition (EMT). Further studies showed that the inhibitory effect of OCF on EMT was associated with the PTEN/PI3K/AKT pathway. Therefore, all studies revealed that OCF could prevent cancer progression and tumor metastasis by inhibiting EMT involved PTEN/PI3K/AKT signaling pathway in lung cancer cells.

Identification of Symmetrical RNA Editing Events in the Mitochondria of Salvia miltiorrhiza by Strand-specific RNA Sequencing

Salvia miltiorrhiza is one of the most widely-used medicinal plants. Here, we systematically analyzed the RNA editing events in its mitochondria. We developed a pipeline using REDItools to predict RNA editing events from stand-specific RNA-Seq data. The predictions were validated using reverse transcription, RT-PCR amplification and Sanger sequencing experiments. Putative sequences motifs were characterized. Comparative analyses were carried out between S. miltiorrhiza, Arabidopsis thaliana and Oryza sativa. We discovered 1123 editing sites, including 225 “C to U” sites in the protein-coding regions. Fourteen of sixteen (87.5%) sites were validated. Three putative DNA motifs were identified around the predicted sites. The nucleotides on both strands at 115 of the 225 sites had undergone RNA editing, which we called symmetrical RNA editing (SRE). Four of six these SRE sites (66.7%) were experimentally confirmed. Re-examination of strand-specific RNA-Seq data from A. thaliana and O. sativa identified 327 and 369 SRE sites respectively. 78, 20 and 13 SRE sites were found to be conserved among A. thaliana, O. sativa and S. miltiorrhiza respectively. This study provides a comprehensive picture of RNA editing events in the mitochondrial genome of S. miltiorrhiza. We identified SREs for the first time, which may represent a universal phenomenon.

Structural and Functional Insights into WRKY3 and WRKY4 Transcription Factors to Unravel the WRKY-DNA (W-Box) Complex Interaction in Tomato (Solanum lycopersicum L.). A Computational Approach

The WRKY transcription factors (TFs), play crucial role in plant defense response against various abiotic and biotic stresses. The role of WRKY3 and WRKY4 genes in plant defense response against necrotrophic pathogens is well-reported. However, their functional annotation in tomato is largely unknown. In the present work, we have characterized the structural and functional attributes of the two identified tomato WRKY transcription factors, WRKY3 (SlWRKY3), and WRKY4 (SlWRKY4) using computational approaches. Arabidopsis WRKY3 (AtWRKY3: NP_178433) and WRKY4 (AtWRKY4: NP_172849) protein sequences were retrieved from TAIR database and protein BLAST was done for finding their sequential homologs in tomato. Sequence alignment, phylogenetic classification, and motif composition analysis revealed the remarkable sequential variation between, these two WRKYs. The tomato WRKY3 and WRKY4 clusters with Solanum pennellii showing the monophyletic origin and evolution from their wild homolog. The functional domain region responsible for sequence specific DNA-binding occupied in both proteins were modeled [using AtWRKY4 (PDB ID:1WJ2) and AtWRKY1 (PDBID:2AYD) as template protein structures] through homology modeling using Discovery Studio 3.0. The generated models were further evaluated for their accuracy and reliability based on qualitative and quantitative parameters. The modeled proteins were found to satisfy all the crucial energy parameters and showed acceptable Ramachandran statistics when compared to the experimentally resolved NMR solution structures and/or X-Ray diffracted crystal structures (templates). The superimposition of the functional WRKY domains from SlWRKY3 and SlWRKY4 revealed remarkable structural similarity. The sequence specific DNA binding for two WRKYs was explored through DNA-protein interaction using Hex Docking server. The interaction studies found that SlWRKY4 binds with the W-box DNA through WRKYGQK with Tyr⁴⁰⁸, Arg⁴⁰⁹, and Lys⁴¹⁹ with the initial flanking sequences also get involved in binding. In contrast, the SlWRKY3 made interaction with RKYGQK along with the residues from zinc finger motifs. Protein-protein interactions studies were done using STRING version 10.0 to explore all the possible protein partners involved in associative functional interaction networks. The Gene ontology enrichment analysis revealed the functional dimension and characterized the identified WRKYs based on their functional annotation.

Genome-wide Identification and Structural, Functional and Evolutionary Analysis of WRKY Components of Mulberry

Mulberry is known to be sensitive to several biotic and abiotic stresses, which in turn have a direct impact on the yield of silk, because it is the sole food source for the silk worm. WRKYs are a family of transcription factors, which play an important role in combating various biotic and abiotic stresses. In this study, we identified 54 genes with conserved WRKY motifs in the Morus notabilis genome. Motif searches coupled with a phylogenetic analysis revealed seven sub-groups as well as the absence of members of Group Ib in mulberry. Analyses of the 2K upstream region in addition to a gene ontology terms enrichment analysis revealed putative functions of mulberry WRKYs under biotic and abiotic stresses. An RNA-seq-based analysis showed that several of the identified WRKYs have shown preferential expression in the leaf, bark, root, male flower, and winter bud of M. notabilis. Finally, expression analysis by qPCR under different stress and hormone treatments revealed genotype-specific responses. Taken together, our results briefs about the genome-wide identification of WRKYs as well as their differential response to stresses and hormones. Importantly, these data can also be utilized to identify potential molecular targets for conferring tolerance to various stresses in mulberry.

Map-based cloning, identification and characterization of the w gene controlling white immature fruit color in cucumber (Cucumis sativus L.)

A single-nucleotide insertion resulted in a premature stop codon that is responsible for white immature fruit color in cucumber. Despite our previous progress in the mapping of the gene controlling white color in immature cucumber fruit and the identification of candidate genes, the specific gene that governs chlorophyll metabolism and its regulatory mechanism remains unknown. Here, we generated a mapping population consisting of 9497 F2 plants to delimit the controlling gene to an 8.2-kb physical interval that defines a sole candidate gene, APRR2. Sequencing the full-length DNA and cDNA of APRR2 allowed for identification of an allele, aprr2, encoding a truncated 101-amino acid protein due to a frameshift mutation and a premature stop codon. Gene structure prediction indicated that these 101 residues are located in a domain necessary for the function of the protein. The expression patterns of APRR2 were entirely consistent with the visual changes in green color intensity during fruit development. A microscopic observation of the fruit pericarp revealed fewer chloroplasts and a lower chloroplast chlorophyll storage capacity in Q24 (white) than in Q30 (green). A single-base insertion in the white color gene w, which leads to a premature stop codon, is hypothesized to have disabled the function of this gene in chlorophyll accumulation and chloroplast development. These findings contribute to basic research and the genetic improvement of fruit color.

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.

Genome-Wide Identification and Characterization of WRKY Gene Family in Peanut

WRKY, an important transcription factor family, is widely distributed in the plant kingdom. Many reports focused on analysis of phylogenetic relationship and biological function of WRKY protein at the whole genome level in different plant species. However, little is known about WRKY proteins in the genome of Arachis species and their response to salicylic acid (SA) and jasmonic acid (JA) treatment. In this study, we identified 77 and 75 WRKY proteins from the two wild ancestral diploid genomes of cultivated tetraploid peanut, Arachis duranensis and Arachis ipaënsis, using bioinformatics approaches. Most peanut WRKY coding genes were located on A. duranensis chromosome A6 and A. ipaënsis chromosome B3, while the least number of WRKY genes was found in chromosome 9. The WRKY orthologous gene pairs in A. duranensis and A. ipaënsis chromosomes were highly syntenic. Our analysis indicated that segmental duplication events played a major role in AdWRKY and AiWRKY genes, and strong purifying selection was observed in gene duplication pairs. Furthermore, we translate the knowledge gained from the genome-wide analysis result of wild ancestral peanut to cultivated peanut to reveal that gene activities of specific cultivated peanut WRKY gene were changed due to SA and JA treatment. Peanut WRKY7, 8 and 13 genes were down-regulated, whereas WRKY1 and 12 genes were up-regulated with SA and JA treatment. These results could provide valuable information for peanut improvement.

Plastoquinone and Ubiquinone in Plants: Biosynthesis, Physiological Function and Metabolic Engineering

Plastoquinone (PQ) and ubiquinone (UQ) are two important prenylquinones, functioning as electron transporters in the electron transport chain of oxygenic photosynthesis and the aerobic respiratory chain, respectively, and play indispensable roles in plant growth and development through participating in the biosynthesis and metabolism of important chemical compounds, acting as antioxidants, being involved in plant response to stress, and regulating gene expression and cell signal transduction. UQ, particularly UQ₁₀, has also been widely used in people's life. It is effective in treating cardiovascular diseases, chronic gingivitis and periodontitis, and shows favorable impact on cancer treatment and human reproductive health. PQ and UQ are made up of an active benzoquinone ring attached to a polyisoprenoid side chain. Biosynthesis of PQ and UQ is very complicated with more than thirty five enzymes involved. Their synthetic pathways can be generally divided into two stages. The first stage leads to the biosynthesis of precursors of benzene quinone ring and prenyl side chain. The benzene quinone ring for UQ is synthesized from tyrosine or phenylalanine, whereas the ring for PQ is derived from tyrosine. The prenyl side chains of PQ and UQ are derived from glyceraldehyde 3-phosphate and pyruvate through the 2-C-methyl-D-erythritol 4-phosphate pathway and/or acetyl-CoA and acetoacetyl-CoA through the mevalonate pathway. The second stage includes the condensation of ring and side chain and subsequent modification. Homogentisate solanesyltransferase, 4-hydroxybenzoate polyprenyl diphosphate transferase and a series of benzene quinone ring modification enzymes are involved in this stage. PQ exists in plants, while UQ widely presents in plants, animals and microbes. Many enzymes and their encoding genes involved in PQ and UQ biosynthesis have been intensively studied recently. Metabolic engineering of UQ₁₀ in plants, such as rice and tobacco, has also been tested. In this review, we summarize and discuss recent research progresses in the biosynthetic pathways of PQ and UQ and enzymes and their encoding genes involved in side chain elongation and in the second stage of PQ and UQ biosynthesis. Physiological functions of PQ and UQ played in plants as well as the practical application and metabolic engineering of PQ and UQ are also included.

Genome-wide analysis, molecular cloning and expression profiling reveal tissue-specifically expressed, feedback-regulated, stress-responsive and alternatively spliced novel genes involved in gibberellin metabolism in Salvia miltiorrhiza

BACKGROUND

Gibberellin (GA), a classical phytohormone, plays significant roles in plant growth and development. It shares the important intermediate diphosphate precursor, GGPP, with the main lipophilic bioactive components, diterpenoid tanshinones in Salvia miltiorrhiza Bunge, one of the most important Traditional Chinese Medicine materials and an emerging model medicinal plant. Analysis of GA metabolism and regulation may help to demonstrate the biological functions of GAs and the crosstalk between GA metabolism and tanshinone biosynthesis in S. miltiorrhiza. However, genes involved in the conversion of ent-kaurene to GAs have not been systematically studied.

RESULTS

Through genome-wide prediction and molecular cloning, twenty two candidate gibberellin metabolism pathway genes were systematically identified for the first time. It includes a SmKO, two SmKAOs, six SmGA20oxs, two SmGA3oxs and eleven SmGA2oxs, of which twenty genes are novel. The deduced proteins showed sequence conservation and divergence. Gibberellin metabolism pathway genes exhibited tissue-specific expression patterns and responded differentially to exogenous GA3 treatment, indicating differential regulation of gibberellin metabolism in different tissue types in S. miltiorrhiza. SmKAO1, SmKAO2, SmGA2ox2, and SmGA2ox4-SmGA2ox7 were significantly up-regulated; SmGA20ox2, SmGA3ox1, SmGA2ox1, SmGA2ox8, SmGA2ox10 and SmGA2ox11 were significantly down-regulated; while the responses of many other genes varied among different tissue-types and time-points of GA3 treatment, suggesting the complexity of feedback regulation. Tanshinone biosynthesis-related genes, such as SmCPS1 and SmKSL1, were up-regulated in response to GA3 treatment. Among the 22 identified genes, nine responded to yeast extract and Ag(+)-treatment in S. miltiorrhiza hairy roots. Moreover, tissue-specifically expressed splice variants were identified for SmKO, SmGA20ox3, SmGA2ox3 and SmGA2ox11, of which SmKOv1, SmGA20ox3v and SmGA2ox11v1 were GA3-responsive, suggesting the importance of alternative splicing in regulating GA metabolism.

CONCLUSIONS

The results show tissue-specifically expressed, feedback-regulated, stress-responsive and alternatively spliced novel genes and reveal multiple layer regulation of GA metabolism and crosstalk between gibberellin metabolism and tanshinone biosynthesis in S. miltiorrhiza.

Promoter analysis of the sweet potato ADP-glucose pyrophosphorylase gene IbAGP1 in Nicotiana tabacum

KEY MESSAGE

The IbAGP1 gene of sweet potato ( Ipomoea batatas ) encodes the sucrose-inducible small subunit of ADP-glucose pyrophosphorylase. Through expression analysis of 5'-truncations and synthetic forms of the IbAGP1 promoter in transgenic tobacco, we show that SURE-Like elements and W-box elements of the promoter contribute to the sucrose inducibility of this gene. Sweet potato (Ipomoea batatas) contains two genes (IbAGP1 and IbAGP2) encoding the catalytically active small subunits of ADP-glucose pyrophosphorylase, an enzyme with an important role in regulating starch synthesis in higher plants. Previous studies have shown that IbAGP1 is expressed in the storage roots, leaves, and stem tissues of sweet potato, and its transcript is strongly induced by applying sucrose exogenously to detached leaves. To investigate the tissue-specific expression of the IbAGP1 promoter, a series of 5'-truncated promoters extending from bases -1913, -1598, -1298, -1053, -716, and -286 to base +75 were used to drive the expression of the β-glucuronidase reporter gene (GUS) in tobacco plants (Nicotiana tabacum). Histochemical and fluorometric GUS assays showed that (1) GUS expression driven by the longest fragment (1989 bp) of the IbAGP1 promoter was detected in vegetative tissues (roots, stems, leaves), (2) fragments extending to -1053 or beyond retained strong GUS expression in roots, stems, and leaves, whereas further 5'-deletions resulted in considerable reduction in GUS activity, and (3) the series of 5'-truncated promoters responded differently to exogenously applied sucrose. The 1989-bp IbAGP1 promoter contains five sequences (two AATAAAA, one AATAAAAAA, and two AATAAATAAA) that are similar to sucrose-responsive elements (SURE). These SURE-Like sequences are found at nucleotide positions -1273, -1239, -681, -610, and -189. Moreover, putative W-box elements are found at positions -1985, -1434, -750, and -578. Synthetic promoters containing tandem repeats of the 4X SURE-Like or 4X W-box upstream from a minimal CaMV35S promoter-GUS fusion showed significant expression in transgenic tobacco in response to exogenous sucrose. These results show that SURE-Like elements and W-box elements of the IbAGP1 promoter contribute to the sucrose inducibility of this gene.

De Novo Assembly and Characterization of the Transcriptome of the Chinese Medicinal Herb, Gentiana rigescens

Gentiana rigescens is an important medicinal herb in China. The main validated medicinal component gentiopicroside is synthesized in shoots, but is mainly found in the plant's roots. The gentiopicroside biosynthetic pathway and its regulatory control remain to be elucidated. Genome resources of gentian are limited. Next-generation sequencing (NGS) technologies can aid in supplying global gene expression profiles. In this study we present sequence and transcript abundance data for the root and leaf transcriptome of G. rigescens, obtained using the Illumina Hiseq2000. Over fifty million clean reads were obtained from leaf and root libraries. This yields 76,717 unigenes with an average length of 753 bp. Among these, 33,855 unigenes were identified as putative homologs of annotated sequences in public protein and nucleotide databases. Digital abundance analysis identified 3306 unigenes differentially enriched between leaf and root. Unigenes found in both tissues were categorized according to their putative functional categories. Of the differentially expressed genes, over 130 were annotated as related to terpenoid biosynthesis. This work is the first study of global transcriptome analyses in gentian. These sequences and putative functional data comprise a resource for future investigation of terpenoid biosynthesis in Gentianaceae species and annotation of the gentiopicroside biosynthetic pathway and its regulatory mechanisms.