Calcium-Dependent Protein Kinase Genes in Glycyrrhiza Uralensis Appear to be Involved in Promoting the Biosynthesis of Glycyrrhizic Acid and Flavonoids under Salt Stress

Diverse Physiological Roles of Flavonoids in Plant Environmental Stress Responses and Tolerance

Flavonoids are characterized as the low molecular weight polyphenolic compounds universally distributed in planta. They are a chemically varied group of secondary metabolites with a broad range of biological activity. The increasing amount of evidence has demonstrated the various physiological functions of flavonoids in stress response. In this paper, we provide a brief introduction to flavonoids' biochemistry and biosynthesis. Then, we review the recent findings on the alternation of flavonoid content under different stress conditions to come up with an overall picture of the mechanism of involvement of flavonoids in plants' response to various abiotic stresses. The participation of flavonoids in antioxidant systems, flavonoid-mediated response to different abiotic stresses, the involvement of flavonoids in stress signaling networks, and the physiological response of plants under stress conditions are discussed in this review. Moreover, molecular and genetic approaches to tailoring flavonoid biosynthesis and regulation under abiotic stress are addressed in this review.

Expression profiling of the mitogen-activated protein kinase gene family reveals their diverse response pattern in two different salt-tolerant Glycyrrhiza species

BACKGROUND

Mitogen-activated protein kinases (MPKs) play important role in response to environmental stress as crucial signal receptors or sensors. Our previous study indicated that salt stress acts as a positive factor to stimulate the production of pharmacodynamic metabolites in the medicinal plant Glycyrrhiza uralensis. Currently, little is known about the MPK gene family and their functions in the medicinal plant G. uralensis.

OBJECTIVE

Identification, comprehensive bioinformatic analysis, expression profiling, and response pattern under salt stress of the G. uralensis GuMPK gene family.

METHODS

Genome-wide investigation and expression profiling of the MPK gene family in G. uralensis, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, promoter cis-acting element, and expression pattern under salt stress in two different salt-tolerant Glycyrrhiza species were performed.

RESULTS

A total of 20 G. uralensis GuMPK genes were identified and categorized into five groups, and had conserved gene structure and motif distribution. Expression profiling of GuMPK genes suggested their potentially diverse functions in plant growth and in response to phytohormones and environmental stress, particularly GuMPK1, 2, 5, and 10 as key components for G. uralensis in response to abiotic stress. Further expression analysis under NaCl treatment in two different salt-tolerant Glycyrrhiza species displayed the MPKs' different response patterns, emphasizing the role of MPK2, 5, 7, and 16 as potentially crucial genes for Glycyrrhiza to respond to salt stress.

CONCLUSION

Our results provide a genome-wide identification and expression profiling of MPK gene family in G. uralensis, and establish the foundation for screening key responsive genes and understanding the potential function and regulatory mechanism of GuMPKs in salt responsiveness.

The Calcium-Dependent Protein Kinase TaCDPK27 Positively Regulates Salt Tolerance in Wheat

As essential calcium ion (Ca2+) sensors in plants, calcium-dependent protein kinases (CDPKs) function in regulating the environmental adaptation of plants. However, the response mechanism of CDPKs to salt stress is not well understood. In the current study, the wheat salt-responsive gene TaCDPK27 was identified. The open reading frame (ORF) of TaCDPK27 was 1875 bp, coding 624 amino acids. The predicted molecular weight and isoelectric point were 68.905 kDa and 5.6, respectively. TaCDPK27 has the closest relationship with subgroup III members of the CDPK family of rice. Increased expression of TaCDPK27 in wheat seedling roots and leaves was triggered by 150 mM NaCl treatment. TaCDPK27 was mainly located in the cytoplasm. After NaCl treatment, some of this protein was transferred to the membrane. The inhibitory effect of TaCDPK27 silencing on the growth of wheat seedlings was slight. After exposure to 150 mM NaCl for 6 days, the NaCl stress tolerance of TaCDPK27-silenced wheat seedlings was reduced, with shorter lengths of both roots and leaves compared with those of the control seedlings. Moreover, silencing of TaCDPK27 further promoted the generation of reactive oxygen species (ROS); reduced the activities of superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT); aggravated the injury to photosystem II (PS II); and increased programmed cell death (PCD) in wheat leaves under NaCl treatment, confirming that the TaCDPK27-silenced seedlings exhibited more NaCl injury than control seedlings. Taken together, the decrease in NaCl tolerance in TaCDPK27-silenced seedlings was due to excessive ROS accumulation and subsequent aggravation of the NaCl-induced PCD. TaCDPK27 may be essential for positively regulating salt tolerance in wheat seedlings.

Integrative analysis of the pharmaceutical active ingredient and transcriptome of the aerial parts of Glycyrrhiza uralensis under salt stress reveals liquiritin accumulation via ABA-mediated signaling

The aerial parts of Glycyrrhiza uralensis supply substantial raw material for the extraction of active pharmaceutical ingredients comprehensively utilized in many industries. Our previous study indicated that salt stress increased the content of active ingredients. However, the regulatory mechanism remains unclear. In this study, RNA-sequencing (RNA-seq) of the aerial parts of G. uralensis treated with 150 mM NaCl for 0, 2, 6, and 12 h was performed to identify the key genes and metabolic pathways regulating pharmacological active component accumulation. The main active component detection showed that liquiritin was the major ingredient and exhibited more than a ten-fold significant increase in the 6 h NaCl treatment. Temporal expression analysis of the obtained 4245 differentially expressed genes (DEGs) obtained by RNA-seq revealed two screened profiles that included the significant up-regulated DEGs (UDEGs) at different treatment points. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of these UDEGs identified phenylpropanoid metabolism and flavonoid biosynthesis as the most significantly enriched pathways in 2 h treated materials. Interestingly, the carotenoid biosynthesis pathway that is related to ABA synthesis was also discovered, and the ABA content was significantly promoted after 6 h NaCl treatment. Following ABA stimulation, the content of liquiritin demonstrated a significant and immediate increase after 2 h treatment, with the corresponding consistent expression of genes involved in the pathways of ABA signal transduction and flavonoid biosynthesis, but not in the pathway of glycyrrhizic acid biosynthesis. Our study concludes that salt stress might promote liquiritin accumulation through the ABA-mediated signaling pathway, and provides effective reference for genetic improvement and comprehensive utilization of G. uralensis.

Combined Analysis of Pharmaceutical Active Ingredients and Transcriptomes of Glycyrrhiza uralensis Under PEG6000-Induced Drought Stress Revealed Glycyrrhizic Acid and Flavonoids Accumulation via JA-Mediated Signaling

Glycyrrhiza uralensis contains many secondary metabolites with a wide range of pharmacological activities. Drought stress acts as a positive regulator to stimulate the production of medicinal active component in G. uralensis, however, the underlying mechanism remains unclear. The aim of this work is to investigate the accumulation and regulatory mechanism of pharmaceutical active ingredients in G. uralensis under drought stress. The materials of the aerial and underground parts of G. uralensis seedlings treated by 10% PEG6000 for 0, 2, 6, 12, and 24 h were used for RNA sequencing and determination of phytohormones and pharmaceutical active ingredients. PEG6000, ibuprofen (IBU), and jasmonic acid (JA) were utilized to treat G. uralensis seedlings for content detection and gene expression analysis. The results showed that, the contents of glycyrrhizic acid, glycyrrhetinic acid, and flavonoids (licochalcone A, glabridin, liquiritigenin, isoliquiritigenin, and liquiritin) were significantly accumulated in G. uralensis underground parts under drought stress. Kyoto Encyclopedia of Genes and Genomes analysis of the transcriptome data of drought-treated G. uralensis indicated that up-regulated differentially expressed genes (UDEGs) involved in glycyrrhizic acid synthesis in the underground parts and flavonoids synthesis in both aerial and underground parts were significantly enriched. Interestingly, the UDEGs participating in jasmonic acid (JA) signal transduction in both aerial and underground parts were discovered. In addition, JA content in both aerial and underground parts under drought stress showed the most significantly accumulated. And drought stress stimulated the contents of JA, glycyrrhizic acid, and flavonoids, coupled with the induced expressions of genes regulating the synthesis and transduction pathway. Moreover, In PEG6000- and JA-treated G. uralensis, significant accumulations of glycyrrhizic acid and flavonoids, and induced expressions of corresponding genes in these pathways, were observed, while, these increases were significantly blocked by JA signaling inhibitor IBU. JA content and expression levels of genes related to JA biosynthesis and signal transduction were also significantly increased by PEG treatment. Our study concludes that drought stress might promote the accumulation of pharmaceutical active ingredients via JA-mediated signaling pathway, and lay a foundation for improving the medicinal component of G. uralensis through genetic engineering technology.

Molecular cloning and functional characterization of UGTs from Glycyrrhiza uralensis flavonoid pathway

Glycyrrhiza uralensis Fisch., a well-known medicinal plant, contains flavonoids including liquiritigenin and isoliquiritigenin, and their corresponding glycoside liquiritin and isoliquiritin. Although some genes encoding UDP-glycosyltransferases (UGTs) have been functionally characterized in G. uralensis, other UGTs mechanisms of glycosylation remain to be elucidated. Against this background the aim of the present study included cloning and characterization of two full-length cDNA clones of GuUGT isoforms from the UGT multigene family. These included GuUGT2 (NCBI acc. MK341791) and GuUGT3 (NCBI acc. MK341793) with an ORF of 1473 and 1332 bp, respectively. Multiple alignments and phylogenetic analysis revealed GuUGTs protein of Glycine max had a high homology to that of other plants. Meanwhile, quantitative real-time PCR was performed to detect the transcript levels of GuUGTs in different tissues. The results indicated that GuUGTs was more expressed in roots as compared to the leaves, and significantly up-regulated upon NaCl stress. The recombinant protein was heterologous expressed in Escherichia coli and exhibited a high level of UGT activity, catalyzing formation of isoliquiritin and liquiritin from isoliquiritigenin and liquiritigenin. The key residues of GuUGT2 for liquiritigenin glycosylation (Asn223), isoliquiritigenin (Asp272) were predicted by molecular docking and residue scanning based on simulated mutations. These results could serve as an important reference to understand the function of the UGT family. In addition, the identification of GuUGT2 and GuUGT3 provides a foundation for future studies of flavonoid biosynthesis in G. uralensis.

Genome-Wide Characterization and Expression Analysis of bZIP Gene Family Under Abiotic Stress in Glycyrrhiza uralensis

bZIP gene family is one of the largest transcription factor families. It plays an important role in plant growth, metabolic, and environmental response. However, complete genome-wide investigation of bZIP gene family in Glycyrrhiza uralensis remains unexplained. In this study, 66 putative bZIP genes in the genome of G. uralensis were identified. And their evolutionary classification, physicochemical properties, conserved domain, functional differentiation, and the expression level under different stress conditions were further analyzed. All the members were clustered into 13 subfamilies (A–K, M, and S). A total of 10 conserved motifs were found in GubZIP proteins. Members from the same subfamily shared highly similar gene structures and conserved domains. Tandem duplication events acted as a major driving force for the evolution of bZIP gene family in G. uralensis. Cis-acting elements and protein–protein interaction networks showed that GubZIPs in one subfamily are involved in multiple functions, while some GubZIPs from different subfamilies may share the same functional category. The miRNA network targeting GubZIPs showed that the regulation at the transcriptional level may affect protein–protein interaction networks. We suspected that domain-mediated interactions may categorize a protein family into subfamilies in G. uralensis. Furthermore, the tissue-specific gene expression patterns of GubZIPs were analyzed using the public RNA-seq data. Moreover, gene expression level of 66 bZIP family members under abiotic stress treatments was quantified by using qRT-PCR. The results of this study may serve as potential candidates for functional characterization in the future.

Comparative transcriptome analysis of leaves during early stages of chilling stress in two different chilling-tolerant brown-fiber cotton cultivars

Chilling stress generates significant inhibition of normal growth and development of cotton plants and lead to severe reduction of fiber quality and yield. Currently, little is known for the molecular mechanism of brown-fiber cotton (BFC) to respond to chilling stress. Herein, RNA-sequencing (RNA-seq)-based comparative analysis of leaves under 4°C treatment in two different-tolerant BFC cultivars, chilling-sensitive (CS) XC20 and chilling-tolerant (CT) Z1612, was performed to investigate the response mechanism. A total of 72650 unigenes were identified with eight commonly used databases. Venn diagram analysis identified 1194 differentially expressed genes (DEGs) with significant up-regulation in all comparison groups. Furthermore, enrichment analyses of COG and KEGG, as well as qRT-PCR validation, indicated that 279 genes were discovered as up-regulated DEGs (UDEGs) with constant significant increased expression in CT cultivar Z1612 groups at the dimensions of both each comparison group and treatment time, locating in the enriched pathways of signal transduction, protein and carbohydrate metabolism, and cell component. Moreover, the comprehensive analyses of gene expression, physiological index and intracellular metabolite detections, and ascorbate antioxidative metabolism measurement validated the functional contributions of these identified candidate genes and pathways to chilling stress. Together, this study for the first time report the candidate key genes and metabolic pathways responding to chilling stress in BFC, and provide the effective reference for understanding the regulatory mechanism of low temperature adaptation in cotton.

Genome-Wide Analysis of MDHAR Gene Family in Four Cotton Species Provides Insights into Fiber Development via Regulating AsA Redox Homeostasis

Monodehydroasorbate reductase (MDHAR) (EC1.6.5.4), a key enzyme in ascorbate-glutathione recycling, plays important roles in cell growth, plant development and physiological response to environmental stress via control of ascorbic acid (AsA)-mediated reduction/oxidation (redox) regulation. Until now, information regarding MDHAR function and regulatory mechanism in Gossypium have been limited. Herein, a genome-wide identification and comprehensive bioinformatic analysis of 36 MDHAR family genes in four Gossypium species, Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, were performed, indicating their close evolutionary relationship. Expression analysis of GhMDHARs in different cotton tissues and under abiotic stress and phytohormone treatment revealed diverse expression features. Fiber-specific expression analysis showed that GhMDHAR1A/D, 3A/D and 4A/D were preferentially expressed in fiber fast elongating stages to reach peak values in 15-DPA fibers, with corresponding coincident observances of MDHAR enzyme activity, AsA content and ascorbic acid/dehydroascorbic acid (AsA/DHA) ratio. Meanwhile, there was a close positive correlation between the increase of AsA content and AsA/DHA ratio catalyzed by MDHAR and fiber elongation development in different fiber-length cotton cultivars, suggesting the potential important function of MDHAR for fiber growth. Following H₂O₂ stimulation, GhMDHAR demonstrated immediate responses at the levels of mRNA, enzyme, the product of AsA and corresponding AsA/DHA value, and antioxidative activity. These results for the first time provide a comprehensive systemic analysis of the MDHAR gene family in plants and the four cotton species and demonstrate the contribution of MDHAR to fiber elongation development by controlling AsA-recycling-mediated cellular redox homeostasis.

CDPK6 phosphorylates and stabilizes MYB30 to promote hyperoside biosynthesis that prolongs the duration of full-blooming in okra

The flowers of okra (Abelmoschus esculentus) open and wilt within only a few hours, and this is accompanied by accumulation of hyperoside, a secondary metabolite in the flavonoid pathway. However, little is known about the relationship between flavonoids and flowering. Here, we found that exogenous application of hyperoside extended the duration of the full-blooming period by more than 3-fold, and this was accompanied by a 14.7-fold increase in the expression of CALCIUM-DEPENDENT PROTEIN KINASE6 (AeCDPK6). Gene expression profiling indicated that the transcription factor AeMYB30 was co-expressed with AeCDPK6, and detailed protein interaction and phosphorylation experiments together with yeast two-hybrid and bimolecular fluorescence complementation assays demonstrated an interaction between AeMYB30 and AeCDPK6. AeCDPK6 specifically phosphorylated AeMYB30S191, leading to increased protein stability and prevention of degradation. Furthermore, AeMYB30 directly bound to the promoter of AeUF3GaT1, a key enzyme in the hyperoside biosynthesis pathway. Analysis of transgenic plants showed that AeCDPK6 was required for the hyperoside-induced phosphorylation of AeMYB30 to enhance its stability and transcriptional activity. Ectopic expression of AeCDPK6 promoted hyperoside accumulation and prolonged the full-blooming period in an AeMYB30-dependent manner. Our results indicate the role of AeCDPK6-AeMYB30 in the molecular mechanism by which hyperoside regulates the period of full blooming in okra, a plant with a short duration of flowering.

Insights on Calcium-Dependent Protein Kinases (CPKs) Signaling for Abiotic Stress Tolerance in Plants

Abiotic stresses are the major limiting factors influencing the growth and productivity of plants species. To combat these stresses, plants can modify numerous physiological, biochemical, and molecular processes through cellular and subcellular signaling pathways. Calcium-dependent protein kinases (CDPKs or CPKs) are the unique and key calcium-binding proteins, which act as a sensor for the increase and decrease in the calcium (Ca) concentrations. These Ca flux signals are decrypted and interpreted into the phosphorylation events, which are crucial for signal transduction processes. Several functional and expression studies of different CPKs and their encoding genes validated their versatile role for abiotic stress tolerance in plants. CPKs are indispensable for modulating abiotic stress tolerance through activation and regulation of several genes, transcription factors, enzymes, and ion channels. CPKs have been involved in supporting plant adaptation under drought, salinity, and heat and cold stress environments. Diverse functions of plant CPKs have been reported against various abiotic stresses in numerous research studies. In this review, we have described the evaluated functions of plant CPKs against various abiotic stresses and their role in stress response signaling pathways.