Next Article in Journal
L-Type Amino Acid Transporter 1 Regulates Cancer Stemness and the Expression of Programmed Cell Death 1 Ligand 1 in Lung Cancer Cells
Next Article in Special Issue
Genetic and Enzymatic Characteristics of CYP2A13 in Relation to Lung Damage
Previous Article in Journal
Biocompatibility and Cytotoxicity of Gold Nanoparticles: Recent Advances in Methodologies and Regulations
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 20
10.3390/ijms222010953
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification of CYP72A Gene Family and Expression Patterns Related to Jasmonic Acid Treatment and Steroidal Saponin Accumulation in Dioscorea zingiberensis

by
Lixiu Hou
,
Xincheng Yuan
,
Song Li
,
Yi Li
,
Zihao Li
and
Jiaru Li
*
State Key Laboratory of Hybrid Rice, Department of Plant Science, College of Life Sciences, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(20), 10953; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222010953
Submission received: 15 September 2021 / Revised: 5 October 2021 / Accepted: 8 October 2021 / Published: 11 October 2021
(This article belongs to the Special Issue Cytochrome P450 (CYP))
Browse Figures
Review Reports Versions Notes

Abstract

:
Dioscorea zingiberensis is a medicinal herb containing a large amount of steroidal saponins, which are the major bioactive compounds and the primary storage form of diosgenin. The CYP72A gene family, belonging to cytochromes P450, exerts indispensable effects on the biosynthesis of numerous bioactive compounds. In this work, a total of 25 CYP72A genes were identified in D. zingiberensis and categorized into two groups according to the homology of protein sequences. The characteristics of their phylogenetic relationship, intron–exon organization, conserved motifs and cis-regulatory elements were performed by bioinformatics methods. The transcriptome data demonstrated that expression patterns of DzCYP72As varied by tissues. Moreover, qRT-PCR results displayed diverse expression profiles of DzCYP72As under different concentrations of jasmonic acid (JA). Likewise, eight metabolites in the biosynthesis pathway of steroidal saponins (four phytosterols, diosgenin, parvifloside, protodeltonin and dioscin) exhibited different contents under different concentrations of JA, and the content of total steroidal saponin was largest at the dose of 100 μmol/L of JA. The redundant analysis showed that 12 DzCYP72As had a strong correlation with specialized metabolites. Those genes were negatively correlated with stigmasterol and cholesterol but positively correlated with six other specialized metabolites. Among all DzCYP72As evaluated, DzCYP72A6, DzCYP72A16 and DzCYP72A17 contributed the most to the variation of specialized metabolites in the biosynthesis pathway of steroidal saponins. This study provides valuable information for further research on the biological functions related to steroidal saponin biosynthesis.

1. Introduction

Cytochromes P450 (CYPs) are a superfamily of heme-containing enzymes that mainly function as monooxygenases in all kingdoms of life. The genes of CYPs in the plant kingdom are usually larger than those in the kingdoms of animals and microorganisms, accounting for approximately 1% of the protein-coding genes [1]. In plants, CYPs are involved in the biosynthesis of diverse specialized metabolites, such as phytohormones, fatty acids and flavonoids, playing a crucial role in plant growth and development [2].
CYPs can be divided into many subfamilies, and CYPs belonging to the same subfamilies have the same catalytic functions in organisms. The CYPs are mainly composed of single-family clans and multiple-family clans based on phylogenetic classification [3]. The CYP72A gene subfamily, an important component of the multiple-family clan in CYPs, catalyzes numerous crucial reactions in the biosynthesis pathway of many specialized metabolites. Gibberellins (GAs), especially GAs with low bioactivity, exert positive influences on seed dormancy. AtCYP72A9 has been confirmed to encode bioactive GA 13-hydroxylase in Arabidopsis thaliana, which contributes to the accumulation of such low-bioactivity GAs [4]. Moreover, a great number of CYP72A genes play indispensable roles in the biosynthesis of medicinal metabolites. In the Medicago genus, MtCYP72A67 encodes a key enzyme that catalyzes hydroxylation at the C-2 position in the hemolytic sapogenin biosynthesis [5]. GsCYP72As have been shown to be involved in the biosynthesis of triterpenoid saponins in Gleditsia sinensis [6]. In addition, the expression of some CYP72A genes can be largely influenced by biotic or abiotic stresses. The methylation level of LpCYP72A161 is closely connected with the response to temperature stress in ryegrass [7].
Dioscorea zingiberensis, a perennial vine, is mainly utilized as a source for the production of diosgenin and has attracted much more interest on account of its pharmaceutical value. Diosgenin is an indispensable precursor to steroidal hormones, and it also shows pharmacological activities against many cancers, such as liver cancer and pancreatic cancer [8,9]. However, natural diosgenin is not abundant in D. zingiberensis, and it will immediately be converted into steroidal saponins by combining with aglycone structures [10]. Steroidal saponins are not only the primary storage form of diosgenin but also the main bioactive compounds in D. zingiberensis [11]. Steroidal saponins exhibited beneficial activities toward decreasing the risk of hyperlipidemia and cardiovascular diseases [12,13]. Previous studies have showed that the diosgenin was biosynthesized by the mevalonate acid (MVA) pathway. Acetyl coenzyme A, the initial substrate in the MVA pathway, can be subsequently converted to phytosterols and diosgenin after being catalyzed via a variety of enzymes. In fenugreek, it was proposed that the CYP72A gene family played a crucial role in the diosgenin biosynthesis and might participate in the biosynthesis pathway of phytosterols [14]. Moreover, enzymes encoded by the CYP72A gene family in legumes were reported to have a catalytic activity in the formation of triterpenoids, such as oleanolic acid, ursolic acid and betulinic acid [15,16]. Nevertheless, research on D. zingiberensis has mainly been focused on developing pharmacological activities [17]. Hence, it is urgent to explore the correlation between the CYP72A gene family and diosgenin biosynthesis.
Considering the importance of CYP72A proteins in many medicinal plants, the objective of this study was to analyze the characters of the CYP72A family using bioinformatics based on the genome of D. zingiberensis, including phylogenetic relationship, gene structure organization, conserved motif analysis, etc. In the meantime, specialized metabolites in the biosynthesis pathway of steroidal saponins (phytosterols, diosgenin and steroidal saponins) were determined, and qRT-PCR was applied to examine the expression difference of the CYP72A gene family under jasmonic acid (JA) treatment. Moreover, important CYP72A genes were screened by investigating the correlation between gene family and steroidal saponin biosynthesis. The results obtained will provide valuable insights into the function of DzCYP72As and identification of gene resources for breeding of D. zingiberensis.

2. Results

2.1. Identification of the CYP72A Proteins in D. zingiberensis

The published CYP72A protein data of A. thaliana were used as the reference sequence against the genomic information of D. zingiberensis to identify DzCYP72A proteins, and 25 sequences were identified as CYP72A proteins (Table S1). Such identified CYP72As were mapped on three chromosomes (Chr1, Chr4 and Chr9) and assigned as DzCYP72A1–DzCYP72A25 based on their genomic location (Table 1). The DzCYP72As were unevenly distributed across the three chromosomes, and three genes of the family were located on unknown chromosomes. The majority of DzCYP72A proteins were localized on Chr1 (n = 14, 56%). By contrast, Chr9 contained seven DzCYP72A genes (28%), and Chr4 contained only one (4%).
The detailed information and physicochemical properties of each DzCYP72A protein were predicted by the ExPASy online tool. The amino acids, molecular weight, isoelectric point, instability index, aliphatic index and grand average of hydropathicity are exhibited in Table 1. The lengths of the DzCYP72A proteins varied from 173 (DzCYP72A19) to 525 (DzCYP72A14) amino acids, with a theoretical molecular weight ranging from 19.63 (DzCYP72A19) to 60.09 (DzCYP72A25) kD. The predicted isoelectric point (pI) values of DzCYP72A proteins ranged from 6.29 (DzCYP72A19) to 9.44 (DzCYP72A14). DzCYP72A15 displayed the better thermostability while DzCYP72A7 showed relatively poor thermostability. According to the predicted results, all CYP72As were a kind of hydrophilic protein (GRAVY < 0), and most of them were unstable proteins. DzCYP72A15 was the only stable protein as its instability index was only 37.99. The subcellular localization of DzCYP72A proteins was predicted by four different online tools, and the predicted results are displayed in Table 2. Most proteins were presumed to be located in membranous organelles, such as the chloroplast, mitochondrion and endoplasmic reticulum.

2.2. Phylogenetic Analysis and Multiple Sequence Alignment of DzCYP72A Genes

The CYP72A proteins in Asparagus officinalis, Solanum lycopersicum, Dioscorea rotundata and Oryza sativa were identified in the same way as DzCYP72A proteins. The phylogenetic relationship among 25 DzCYP72As, 6 AoCYP72As, 13 SlCYP72As, 6 DrCYP72As, 9 OsCYP72As and 9 AtCYP72As was analyzed by MEGA7 based on the aligned results of protein sequences (Figure 1). In summary, all CYP72A proteins could be divided into two categories: monocotyledons and dicotyledons. In the meantime, the DzCYP72A family was categorized into two clades. A total of eight DzCYP72As, together with two DrCYP72As and five AoCYP72As were categorized as clade Ⅰ. Clade Ⅱ contained the remaining 17 DzCYP72As, 6 OsCYP72As and 1 AoCYP72A. Clade Ⅲ comprised four DrCYP72As and three OsCYP72As. All AtCYP72As (n = 9) and the majority of SlCYP72As (n = 12) were categorized as clade Ⅳ. Multiple alignment of DzCYP72A proteins was performed using DNAMAN8 software (Figure 2). The DzCYP72A protein sequence in clade Ⅱ had higher homology than that in clade Ⅰ. All DzCYP72A proteins could be divided into two clusters, which kept consistence with the phylogenetic results.

2.3. Gene Structure and Motif Analysis of DzCYP72As

The exon/intron composition of DzCYP72As displayed the structural diversity and complexity of this gene family (Figure 3). The exon number of DzCYP72As in clade Ⅰ ranged from one to six. DzCYP72A18 had the largest gene length, and it had the maximum number of exons. In contrast, the gene length of DzCYP72A19 was the shortest, and it had only two exons. All DzCYP72As in clade Ⅱ contained five exons, and DzCYP72A5 was intron free. The coding sequence of DzCYP72As is listed in Table S2.
In the meanwhile, the conserved motifs of DzCYP72A proteins were further analyzed by the MEME online search tool. A total of 10 motifs were found, and the detailed information is displayed in Figure 4. Based on the protein sequence of all DzCYP72As, we found that all DzCYP72A proteins contained motif 1, implying that motif 1 might be one of conserved motifs among all CYP72A proteins in D. zingiberensis. The motif number of all DzCYP72A proteins was in the range of 4–10, and motif 1 and 2 were the common composition in clade Ⅰ. The motif 1, 2, 8 and 10 were fundamental components of the DzCYP72A proteins in clade Ⅱ. The DzCYP72As in clade Ⅱ contained more motifs than those in clade Ⅰ, suggesting that such diverse motifs may be strongly correlated with corresponding functions.

2.4. Analysis of Cis-Regulatory Elements in the Promoters of DzCYP72As

Cis-regulatory elements (CREs) in the promoter region exert a crucial influence on gene functions. In order to further conduct research on the genetic functions and regulatory mechanisms of the DzCYP72A gene family, the 2 kb upstream sequences of 25 DzCYP72As were uploaded to the PlantCARE online tool. A total of 10 representative CREs (light responsive, auxin responsive, abscisic acid responsive, MYB binding site, MeJA responsive, low-temperature responsive, anaerobic induction, gibberellin responsive, salicylic acid responsive and defense and stress responsive) were visualized on each gene and are shown in Figure 5. Among all CREs, the proportion of light-responsive elements was the largest, accounting for 47%. The MeJA-responsive CREs were the second most abundant and stood at 15%. The proportions of anaerobic induction, abscisic acid responsive and MYB binding site were similar, with 9%, 8% and 7%, respectively. The proportions of auxin-responsive, gibberellin-responsive and salicylic-acid-responsive CREs were at the same level, approximately making up 3%. Additionally, the CRE related to defense and stress responsiveness occupied the least proportion, only 2%. The detailed information on CREs is exhibited in Table S3.

2.5. Expression Profiles of CYP72As in D. zingiberensis

In order to further explore the functions of DzCYP72A genes, we analyzed the transcriptome data and investigated the tissue-specific expression patters of each DzCYP72A gene. The expression patterns of DzCYP72A genes varied by tissues (Figure 6A). DzCYP72A23 and DzCYP72A15 displayed the highest transcript level in both the leaf and the stem, whereas DzCYP72A15 exhibited the highest transcript level in the rhizome. DzCYP72A4 and DzCYP72A5 showed the lowest expression level in the leaf, while DzCYP72A14 represented the lowest expression level in the stem. In addition, both DzCYP72A25 and DzCYP72A19 exhibited the lowest expression level in the rhizome. These results implied that DzCYP72As may have indispensable influences on the growth and development of D. zingiberensis.
Numerous studies have illustrated that JA and its derivates exerted positive effects on the biosynthesis of saponins [18,19,20]. Moreover, the DzCYP72A gene family played a key role in the biosynthesis of steroidal saponins and the analysis of cis-regulatory elements suggested that DzCYP72As had a strong correlation with JA. Thus, the primers of DzCYP72As were designed (Table S4) and qRT-PCR was applied to analyze the transcriptional expressions of DzCYP72As under treatment with different concentrations of JA (Figure 6B). Almost all DzCYP72As displayed diverse expression patterns under all treatments, but it was difficult to obtain acceptable qRT-PCR results for gene expression analysis of DzCYP72A18, DzCYP72A19 and DzCYP72A25 in rhizomes on account of their low expression abundance and amplification efficiency, which cohered with the transcriptome data.
Some genes, such as DzCYP72A24, DzCYP72A13 and DzCYP72A7, responded rapidly and displayed higher expression levels at a relatively low concentration of JA (25 μmol/L). The transcriptional expressions of DzCYP72A13, DzCYP72A14, DzCYP72A16, DzCYP72A20 and DzCYP72A24 showed significant upregulation. The majority of DzCYP72A genes exhibited the highest expression level under JA with 100 μmol/L.

2.6. Effects of JA Concentration on the Specialized Metabolites in D. zingiberensis

Phytosterols, such as cholesterol, campesterol, stigmasterol and β-sitosterol, were reported as the intermediates in the biosynthesis pathway of steroidal saponins [21]. Therefore, the contents of specialized metabolites, including phytosterols, diosgenin and steroidal saponins, were determined to investigate the impacts of JA on steroidal saponin biosynthesis. It has been found that JA treatment exerted significant effects on the contents of bioactive compounds compared to those in untreated rhizomes (Table 3). Both dioscin and protodeltonin responded rapidly to low concentrations of JA, exhibiting the highest content of 38.25 μg/g and 15.69 mg/g, respectively, at JA with 50 μmol/L. The accumulation of parvifloside showed a positive relationship with the higher concentrations of JA and was approximately 1.6-fold increase over that in the control group at a dose of 100 μmol/L. Therefore, the total saponins obtained the maximum content (69.26 mg/g) at JA with 100 μmol/L, which was 1.8-fold higher than in the untreated group. Likewise, the content of diosgenin also positively responded to the JA and obtained the highest yield (647.18 µg/g) at a dose of 100 μmol/L. However, the content of cholesterol showed no significant difference under different concentrations of JA. In contrast, the yield of campesterol and β-sitosterol varied with different concentrations of JA, and the maximum yields of these two metabolites all occurred at a dose of 100 μmol/L. Stigmasterol obtained the maximum content at a dose of 50 μmol/L of JA, and its accumulation was decreased by relatively higher concentrations of JA.

2.7. Correlations among DzCYP72As, Phytosterols, Diosgenin and Steroidal Saponins

The de-trended correspondence analysis (DCA) was performed on phytosterols, diosgenin and steroidal saponins, and then was applied to calculate the maximum length of the gradient axis (LGA). Table S5 shows that the LGA value of DCA was below one. According to previous research, the redundancy analysis (RDA) is suitable for subsequent analyses when the LGA values are <3 [22]. Therefore, RDA was carried out to analyze the correlation between specialized metabolites and DzCYP72As. In the RDA algorithm, a total of 22 DzCYP72As responding to JA were used as the explanatory variables, and 8 metabolites were used as the response variables. In the meantime, a forward direction of the Akaike information criterion (AIC) was applied to further screen explanatory variables. In the first two axes of RDA analysis, these genes constrained 91.46% of variance in eight metabolites, suggesting these two axes could represent the total constrained proportion (Table S6). As shown in Table 4, a total of 12 DzCYP72As had correlations with specialized metabolites of steroidal saponin biosynthesis. However, only 9 DzCYP72As were significantly correlated with those metabolites. Among these, DzCYP72A6, DzCYP72A16 and DzCYP72A17 displayed a highly significant correlation with such bioactive compounds. DzCYP72A1, DzCYP72A3, DzCYP72A9, DzCYP72A11, DzCYP72A14 and DzCYP72A20 showed significant correlations with specialized metabolites. In view of Figure 7, the screened DzCYP72As were positively associated with parvifloside, protodeltonin, dioscin, diosgenin, campesterol and β-sitosterol. In contrast, cholesterol and stigmasterol were negatively correlated with these DzCYP72As.

3. Discussion

The CYP72A gene family plays crucial roles in catalyzing many considerable reactions in Plantae, and many CYP72A genes have been cloned from various plants [4,5]. Nevertheless, genome-wide analysis of the CYP72A gene family has not been performed in D. zingiberensis, an important source of diosgenin. In this study, a total of 25 DzCYP72A genes were identified in D. zingiberensis and assigned as DzCYP72A1–25 on the base of their chromosomal location. The phylogenetic relationship, gene structure, conserved motifs, cis-regulatory elements and expression patterns under JA treatment were analyzed. Meanwhile, eight specialized metabolites related to the biosynthesis of steroidal saponins were determined, and the correlation between DzCYP72As and such metabolites was also investigated. This work provided valuable information for subsequent functional analysis of DzCYP72As.
The phylogenetic tree indicated that DzCYP72A proteins had stronger homology with O. sativa, A. officinalis and D. rotundata than with A. thaliana (Figure 1). D. zingiberensis belongs to monocotyledon and A. thaliana belongs to dicotyledon; therefore, O. sativa, A. officinalis and D. rotundata, belonging to monocotyledon, are inclined to have strong homology with D. zingiberensis. Moreover, the functions of CYP gene subfamilies catalyze considerable reactions in the biosynthesis pathway of numerous specialized metabolites in plants, and previous research illustrated that the CYP72A gene family had a profound correlation with the biosynthesis of saponins [23,24]. Hence, the reason for D. zingiberensis containing more DzCYP72A genes may be the large accumulation of steroidal saponins in the rhizome of D. zingiberensis. Likewise, the large amount of steroidal saponins in A. officinalis may be one of the reasons why all AoCYP72A proteins had strong homology with DzCYP72A proteins. Although D. rotundata also belongs to the Dioscorea genus, it contains abundant starch and only a trace amount of steroidal saponins [25]. Therefore, the large amount of starch may be closely related to the high homology of amino acid sequences between DrCYP72A and OsCYP72A proteins. In addition, S. lycopersicum also belongs to dicotyledon, but one SlCYP72A protein has remarkable homology with DzCYP72A proteins. According to previous research, S. lycopersicum contains abundant steroidal glycoalkaloids, which are also biosynthesized from phytosterols [26]. The similar biosynthesis pathway may be connected to the high homology between SlCYP72A and DzCYP72A proteins.
Similar to the case in other gene families in plants, the expression profiles of DzCYP72As displayed special variations in different tissues, which may be closely related to their functions in D. zingiberensis [27,28]. Many DzCYP72A genes had high expressions in the leaf and stem, but the steroidal saponins mainly accumulated in the rhizomes. According to previous research, the CYP72A gene family mainly participates in the biosynthesis of phytosterols, which are not only the intermediates in the biosynthesis pathway of steroidal saponins but also indispensable components for cell membranes [29,30]. Generally, the leaves have higher metabolism than rhizomes, and rhizomes cultivated for at least 3 years can satisfy the industrial requirements. Therefore, the reason for DzCYP72A genes having a high transcript level might be that more phytosterols are needed for maintaining the functions of cell membranes in leaves under a high metabolic rate.
Abiotic stresses, such as extreme temperature, abscisic acid, salicylic acid, JA and its derivates, exert considerable effects on plant development and the accumulation of specialized metabolisms [31,32]. In this study, the MeJA-responsive CREs were the most abundant among all CREs responsive to phytohormones, indicating that JA and its derivates profoundly correlated with the functions of DzCYP72As. Moreover, previous research demonstrated that JA and its derivates had positive correlations with the biosynthesis of steroidal saponins [18,19,20]. Therefore, most DzCYP72As were upregulated under different concentrations of JA. In most cases, MeJA was used as an elicitor in medicinal crops for enhancing the yields of bioactive compounds [33]. However, MeJA is essentially a non-bioactive compound, but it is a volatile substance, which can be helpful to increase the capacity of JA to easily enter plants via the stomata [34]. In the meantime, JA is more effective than the methylation of JA in stimulating the biosynthesis of saponins [19]. Therefore, JA was applied to the rhizomes of D. zingiberensis in this work.
According to previous research, Ankang of Shaanxi and Shiyan of Hubei are the original cities of D. zingiberensis, and the rhizomes from these two regions contained more steroidal saponins than those from other locations in China [21,35]. Therefore, seedlings from Ankang of Shaanxi were used as materials to investigate the effects of JA on DzCYP72As and specialized metabolites. The content of dioscin in this study was lower than that in mature rhizomes, but the contents of cholesterol, stigmasterol, β-sitosterol and diosgenin were higher than that in mature rhizomes [21]. Solar energy provides the source of energy for plants and is stored via photosynthesis [36]. When a plant reaches natural maturation, the aerial part of D. zingiberensis has already withered, inhibiting photosynthesis and energy storage. Meanwhile, the low metabolic rate and insufficient energy supply may pose negative influences on the biosynthesis of steroidal saponins, which might be one of the reasons why mature rhizomes contained only a trace amount of cholesterol, stigmasterol, β-sitosterol and diosgenin. Cholesterol has been proved to be one of the important compounds that provide carbon skeletons for the biosynthesis of steroidal saponins [37]. Nevertheless, the content of cholesterols had no correlation with JA in this study. In contrast, β-sitosterol had a profound correlation with JA. In the meanwhile, it was found that β-sitosterol can also biosynthesize the diosgenin and steroidal saponins [38]. Thus, the reason for the increase of steroidal saponins might be that a large amount of β-sitosterol participates in the biosynthesis of steroidal saponins under JA treatment. Moreover, β-sitosterol is also reported to be the precursor for stigmasterol, but the content of stigmasterol decreased with increased β-sitosterol in the high concentrations of JA, implying the accumulated β-sitosterol may be used for the biosynthesis of steroidal saponins [39]. In most cases, various bioactive compounds had different pharmacological activities, and parvifloside, protodeltonin and dioscin showed diverse responses to different concentrations of JA, which provided a targeted scientific basis for enhancing those metabolites. In addition, other gene families may also have strong correlations with the biosynthesis of steroidal saponins. Meanwhile, other environmental variables, such as soil salinization and moisture, can also affect expression patterns of DzCYP72A genes and other gene families. Therefore, further research should pay more attention to the correlation between specialized metabolites and other gene families under different environmental conditions.

4. Materials and Methods

4.1. Plant Materials

The seeds of D. zingiberensis were obtained from a plantation of Shaanxi Ankang in October of 2020. D. zingiberensis seeds were cultivated in an environmentally controlled greenhouse with 26 ± 2 ℃ and 16 h light/day. According to the manufacturer’s instruction and previous research, jasmonic acid (JA, Macklin, 98%) was dissolved in 80% ethanol and diluted in 1/2 Hoagland solution to obtain 25, 50, 100 and 200 μmol/L solutions [40,41]. Ethanol was diluted in 1/2 Hoagland solution and used as a control solution. Three-month-old seedlings with similar growth status were watered with enough control solution and JA solutions for six days. The seedlings were watered every three days, and all solutions were freshly prepared just before use. The rhizomes were collected, and every treatment was replicated three times. We finally obtained the JA-treated samples: 0 μmol/L (S1), 25 μmol/L (S2), 50 μmol/L (S3), 100 μmol/L (S4) and 200 μmol/L (S5). Some samples were promptly frozen in liquid nitrogen pending RNA extraction. While the others were freeze dried to constant weight at −80 ℃ and ground to powders with a tissue lyser.

4.2. Standards and Chemical Reagents

Methanol, acetonitrile, chloroform, ethanol and n-hexane were purchased from Thermo Fisher Scientific Inc. (Waltham, MA, USA), and all chemical reagents were of HPLC grade. Dioscin (98%), diosgenin (98%), β-sitosterol (98%), stigmasterol (95%), campesterol (98%), cholesterol (98%) and N-methyl-N-(trimethylsilyl) trifluoroacetamide were purchased from Shanghai Yuanye Bio-Technology Co., Ltd. (Shanghai, China). Parvifloside (95%) and protodeltonin (95%) were obtained from Kunming Institute of Botany, CAS (Kunming, China). Detailed information and the structures of eight compounds can be obtained on the NCBI website (https://pubchem.ncbi.nlm.nih.gov/ accessed on 30 July 2021).

4.3. RNA Extraction and Gene Expression Analysis

The total RNA of rhizomes was extracted using RNAiso plus (Takara, Beijing, China) according to its instructions. RNA quality and concentration were determined with the NanoDrop 2000 (Thermo Scientific, Waltham, MA, USA). Subsequently, 1 μg of total RNA was subjected to a reverse transcription reaction according to HiScript III 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China).
All primers were designed by Primer 5.0 and are listed in Table S4. Quantitative real-time PCR (qRT-PCR) was conducted according to the instructions of 2 × TransStart Green qPCR SuperMix (TransGen, Beijing, China). The expression patterns of the DzCYP72A gene family under JA treatment were analyzed using the CFX96 Real-Time PCR Detection System (Bio-Rad, Los Angeles, CA, USA). DzActin and DzGAPDH were used as internal controls. The gene’s relative expression was normalized by the 2−ΔΔCt method.

4.4. Identification and Screening of CYP72A Family Genes

The published nine CYP72A protein sequences of A. thaliana were downloaded from the database (http://www.p450.kvl.dk/At_cyps/family.shtml#72A accessed on 10 July 2021) [42]. The basic local alignment search tool (BLAST) was used to identify and screen DzCYP72A proteins. In brief, these nine published AtCYP72A protein sequences were used as a query to further screen the CYP72A proteins in D. zingiberensis, and 25 sequences fulfilling requirements were finally identified. The genome information of D. zingiberensis was uploaded to the NCBI database under project PRJNA716093 (unpublished).
The genome information of A. officinalis, S. lycopersicum (TAG 3.2) and O. sativa (V7.0) was downloaded from the Joint Genome Institute (https://phytozome.jgi.doe.gov accessed on 10 July 2021). In the meantime, the genome information of D. rotundata was downloaded from the NCBI (project number: PRJNA695139). The CYP72A protein sequences of these species were confirmed in the same way as for DzCYP72A proteins.

4.5. Bioinformatics Analysis

The physicochemical characteristics of DzCYP72A proteins were predicted by an online tool, ExPASy (https://web.expasy.org/protparam/ accessed on 15 July 2021). The subcellular localization information of DzCYP72A proteins was obtained from the Plant-PLoc website (http://www.csbio.sjtu.edu.cn/bioinf/plant/ accessed on 18 July 2021), the WoLF PSORT website (https://wolfpsort.hgc.jp/ accessed on 18 July 2021), the CELLO website (http://cello.life.nctu.edu.tw/ accessed on 18 July 2021) and the YLoc website (https://abi-services.informatik.uni-tuebingen.de/yloc/webloc.cgi accessed on 18 July 2021) [43]. The CYP72A protein sequences of A. thaliana, D. zingiberensis, A. officinalis, S. lycopersicum, O. sativa and D. rotundata were aligned by MEGA 7.0 and visualized using the online tool iTOL (https://itol.embl.de/ accessed on 23 July 2021). The gene structure of DzCYP72A was calculated by TBtools and visualized by the online website Gene Structure Display Server 2.0 (http://gsds.gao-lab.org/ accessed on 25 July 2021). [44]. The conserved motifs of DzCYP72A proteins were analyzed using the online tool Multiple EM for Motif Elicitation (MEME, https://meme-suite.org/meme/tools/meme accessed on 28 July 2021) and visualized using the TBtools [45,46]. The 2 kb upstream sequences of each DzCYP72A gene were predicted using the PlantCARE software (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ accessed on 30 July 2021). and presented using TBtools [47]. The transcriptome data of D. zingiberensis (leaf, stem and rhizome) were uploaded to NCBI database project SRR15660571, SRR15660569 and SRR15660570, respectively (unpublished).

4.6. Analysis of Specialized Metabolites

About 50 mg freeze-dried rhizome powder was saponified with 2 mL 2 mol/L KOH-ethanol in the water bath for one hour at 80 ℃. The unsaponifiable part was extracted using n-hexane and filtered through a membrane solution filter (0.22 μm). The supernatant was evaporated under vacuum at room temperature. Dried samples were derivatized by 50 μL N-methyl-N-(trimethylsilyl) trifluoroacetamide in room temperature for 30 min. Subsequently, 150 μL n-hexane was added pending the determination of phytosterols.
A 2:1 chloroform/methanol mixture solvent was used to extract diosgenin. Approximately 50 mg freeze-dried rhizome powder was dissolved in 1 mL of the mixture solvent and extracted by ultrasonication for 30 min. Diosgenin was obtained by centrifugation at 12,000× g for 10 min. The extraction of each sample was repeated three times. The supernatant extracted three times was merged and filtered through a membrane solution filter (0.22 μm) pending determination. The contents of phytosterols and diosgenin were detected using the Thermo Trace gas chromatography–mass spectrometry instrument equipped with the TG-5 ms column (30 m × 0.25 mm × 0.25 μm). The operation conditions were performed as described previously [21].
Subsequently, 80% ethanol was used for the extraction of steroidal saponins. Approximately 50 mg freeze-dried rhizome powder was dissolved in 1 mL extraction solution and centrifuged at 12,000× g for 10 min. This powder (50 mg) was extracted three times. The obtained supernatant was merged and filtered through a membrane solution filter (0.22 μm) for subsequent analysis. The steroidal saponins were determined by an ultra-performance liquid chromatography instrument equipped with the Q Exactive hybrid quadrupole mass spectrometer (Thermo Scientific, Waltham, MA, USA) and a reverse-phase C18 column (Thermo, 100 mm × 2.1 mm, 3 μm). The mobile phases, flow rate and temperature of chromatography column were operated as described previously [21]. The gradient program was: 0–2 min, isocratic 12% B; 2–8 min, linear gradient of 12–35% B; 8–20 min, isocratic 45% B; 20–20.1 min, linear gradient of 45–20% B; 20.1–22 min, isocratic 20% B. The data collection was performed by full scan mode (m/z 300–200) and selected ion monitoring mode with the diagnostic ion monitored at m/z 869.48.

4.7. Statistical Analyses

The qualitative information regarding phytosterols, diosgenin and steroidal saponins was collected from mass spectra, and the quantified results were calculated by the peak areas of external standards (Table S7). The mean and standard deviation of each specialized metabolite were determined using SPSS 19.0 software. DCA was performed with the package vegan in R.4.0.2, and RDA was carried out on the online tool Gene Denovo (https://www.omicshare.com/tools/ accessed on 15 August 2021). The correlation between DzCYP72As and steroidal saponin biosynthesis was analyzed using Rstudio.

5. Conclusions

In this study, a total of 25 CYP72A genes were screened and isolated from the genome of D. zingiberensis. The physicochemical characteristics, subcellular localization, phylogenetic analysis, exon–intron organization, motifs, cis-regulatory elements and tissue-specific expression were investigated with diverse bioinformatics methods. The results of qRT-PCR and eight metabolites revealed that DzCYP72As and specialized metabolites displayed significant responses to JA treatment. Moreover, the total steroidal saponins, parvifloside, natural diosgenin campesterol and β-sitosterol were the most abundant at a dose of 100 μmol/L of JA, while protodeltonin, dioscin and stigmasterol were the most abundant at the concentration of 50 μmol/L of JA. The spearman analysis revealed that DzCYP72As have a strong correlation with eight metabolites in the biosynthesis pathway of steroidal saponins. The obtained results provide useful information in relation to the DzCYP72 gene family and steroidal saponins, which is beneficial for further investigations into the evolution and functions of these genes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms222010953/s1.

Author Contributions

Conceptualization, J.L. and L.H.; methodology, L.H. and X.Y.; software, L.H., Y.L. and Z.L.; validation, X.Y., Y.L. and S.L.; writing—original draft preparation, L.H.; writing—review and editing, L.H., S.L. and Z.L.; supervision, J.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 31270345.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All study data are included in the main text and Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Nelson, D.; Werck-Reichhart, D. A P450-centric view of plant evolution. Plant J. 2011, 66, 194–211. [Google Scholar] [CrossRef] [PubMed]
  2. Feldmann, K.A. Cytochrome P450s as genes for crop improvement. Curr. Opin. Plant Biol. 2001, 4, 162–167. [Google Scholar] [CrossRef]
  3. Hansen, C.C.; Nelson, D.R.; Moller, B.L.; Werck-Reichhart, D. Plant cytochrome P450 plasticity and evolution. Mol. Plant 2021, 14, 1244–1265. [Google Scholar] [CrossRef] [PubMed]
  4. He, J.; Chen, Q.; Xin, P.; Yuan, J.; Ma, Y.; Wang, X.; Xu, M.; Chu, J.; Peters, R.J.; Wang, G. CYP72A enzymes catalyse 13-hydrolyzation of gibberellins. Nat. Plants 2019, 5, 1057–1065. [Google Scholar] [CrossRef] [Green Version]
  5. Biazzi, E.; Carelli, M.; Tava, A.; Abbruscato, P.; Losini, I.; Avato, P.; Scotti, C.; Calderini, O. CYP72A67 Catalyzes a Key Oxidative Step in Medicago truncatula Hemolytic Saponin Biosynthesis. Mol. Plant 2015, 8, 1493–1506. [Google Scholar] [CrossRef] [Green Version]
  6. Kuwahara, Y.; Nakajima, D.; Shinpo, S.; Nakamura, M.; Kawano, N.; Kawahara, N.; Yamazaki, M.; Saito, K.; Suzuki, H.; Hirakawa, H. Identification of potential genes involved in triterpenoid saponins biosynthesis in Gleditsia sinensis by transcriptome and metabolome analyses. J. Nat. Med. 2018, 73, 369–380. [Google Scholar] [CrossRef] [Green Version]
  7. Dai, Y.; Mao, P.; Tao, X.; Wang, Y.; Wei, C.-M.; Ma, X.-R. Flexible shift on gene body methylation and transcription of LpCYP72A161 exposed to temperature stress in perennial ryegrass. Environ. Exp. Bot. 2017, 143, 29–37. [Google Scholar] [CrossRef]
  8. Chang, H.-Y.; Kao, M.-C.; Way, T.-D.; Ho, C.-T.; Fu, E. Diosgenin Suppresses Hepatocyte Growth Factor (HGF)-Induced Epithelial–Mesenchymal Transition by Down-regulation of Mdm2 and Vimentin. J. Agric. Food Chem. 2011, 59, 5357–5363. [Google Scholar] [CrossRef]
  9. Zhang, J.; Xie, J.-J.; Zhou, S.-J.; Chen, J.; Hu, Q.; Pu, J.-X.; Lu, J.-L. Diosgenin inhibits the expression of NEDD4 in prostate cancer cells. Am. J. Transl. Res 2019, 11, 3461–3471. [Google Scholar]
  10. Moses, T.; Papadopoulou, K.K.; Osbourn, A. Metabolic and functional diversity of saponins, biosynthetic intermediates and semi-synthetic derivatives. Crit. Rev. Biochem. Mol. Biol. 2014, 49, 439–462. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, X.-X.; Ito, Y.; Liang, J.-R.; Liu, J.-L.; He, J.; Sun, W.-J. Therapeutic effects of total steroid saponin extracts from the rhizome of Dioscorea zingiberensis C.H.Wright in Freund’s complete adjuvant induced arthritis in rats. Int. Immunopharmacol. 2014, 23, 407–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hashidume, T.; Sasaki, K.; Hirata, J.; Kato, M.; Yoshikawa, Y.; Iwasaki, Y.; Arai, H.; Miura, S.; Miyoshi, N. Effects of sanyaku and its constituent diosgenin on the fasted and postprandial hpertriacylglycerolemia in high-fat-diet-fed KK-Ay Mice. J. Agric. Food Chem. 2018, 66, 9968–9975. [Google Scholar] [CrossRef] [PubMed]
  13. Tang, Y.-N.; He, X.-C.; Ye, M.; Huang, H.; Chen, H.-L.; Peng, W.-L.; Zhao, Z.-Z.; Yi, T.; Chen, H.-B. Cardioprotective effect of total saponins from three medicinal species of Dioscorea against isoprenaline-induced myocardial ischemia. J. Ethnopharmacol. 2015, 175, 451–455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ciura, J.; Szeliga, M.; Grzesik, M.; Tyrka, M. Changes in fenugreek transcriptome induced by methyl jasmonate and steroid precursors revealed by RNA-Seq. Genomics 2018, 110, 267–276. [Google Scholar] [CrossRef]
  15. Fukushima, E.O.; Seki, H.; Sawai, S.; Suzuki, M.; Ohyama, K.; Saito, K.; Muranaka, T. Combinatorial Biosynthesis of Legume Natural and Rare Triterpenoids in Engineered Yeast. Plant Cell Physiol. 2013, 54, 740–749. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, J.; Li, Y.; Zhang, Y.; Jia, L.; Sun, L.; Wang, S.; Xiao, J.; Zhan, Y.; Yin, J. Functional identification of five CYP450 genes from birch responding to MeJA and SA in the synthesis of betulinic acid from lupitol. Ind. Crop. Prod. 2021, 167, 113513. [Google Scholar] [CrossRef]
  17. Xue, Z.; Cao, Z.; Jin, M.; Zhang, X.; Wang, X.; Dou, J.; Zhu, Y.; Ito, Y.; Guo, Z. New steroid saponins from Dioscorea zingiberensis yam and their medicinal use against I/R via anti-inflammatory effect. Food Funct. 2021, 12, 8314–8325. [Google Scholar] [CrossRef] [PubMed]
  18. Alsoufi, A.S.M.; Pączkowski, C.; Szakiel, A.; Długosz, M. Effect of jasmonic acid and chitosan on triterpenoid production in Calendula officinalis hairy root cultures. Phytochem. Lett. 2019, 31, 5–11. [Google Scholar] [CrossRef]
  19. Li, J.; Wang, J.; Wu, X.; Liu, D.; Li, J.; Li, J.; Liu, S.; Gao, W. Jasmonic acid and methyl dihydrojasmonate enhance saponin biosynthesis as well as expression of functional genes in adventitious roots of Panax notoginseng F.H. Chen. Biotechnol. Appl. Biochem. 2017, 64, 225–238. [Google Scholar] [CrossRef]
  20. Chaudhary, S.; Chikara, S.K.; Sharma, M.C.; Chaudhary, A.; Alam Syed, B.; Chaudhary, P.S.; Mehta, A.; Patel, M.; Ghosh, A.; Iriti, M. Elicitation of Diosgenin Production in Trigonella foenum-graecum (Fenugreek) Seedlings by Methyl Jasmonate. Int. J. Mol. Sci. 2015, 16, 29889–29899. [Google Scholar] [CrossRef] [Green Version]
  21. Hou, L.; Li, S.; Tong, Z.; Yuan, X.; Xu, J.; Li, J. Geographical variations in fatty acid and steroid saponin biosynthesis in Dioscorea zingiberensis rhizomes. Ind. Crop. Prod. 2021, 170, 113779. [Google Scholar] [CrossRef]
  22. Hill, M.O.; Gauch, H.G. Detrended correspondence analysis: An improved ordination technique. Vegetatio 1980, 42, 47–58. [Google Scholar] [CrossRef]
  23. Abdollahi, F.; Alebrahim, M.T.; Ngov, C.; Lallemand, E.; Zheng, Y.; Villette, C.; Zumsteg, J.; André, F.; Navrot, N.; Werck-Reichhart, D.; et al. Innate promiscuity of the CYP706 family of P450 enzymes provides a suitable context for the evolution of dinitroaniline resistance in weed. New Phytol. 2021, 229, 3253–3268. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, Q.; Khakimov, B.; Cárdenas, P.D.; Cozzi, F.; Olsen, C.E.; Jensen, K.R.; Hauser, T.P.; Bak, S. The cytochrome P450 CYP72A552 is key to production of hederagenin-based saponins that mediate plant defense against herbivores. New Phytol. 2019, 222, 1599–1609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Tamiru, M.; Natsume, S.; Takagi, H.; White, B.; Yaegashi, H.; Shimizu, M.; Yoshida, K.; Uemura, A.; Oikawa, K.; Abe, A.; et al. Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination. BMC Biol. 2017, 15, 1–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Nakayasu, M.; Umemoto, N.; Akiyama, R.; Ohyama, K.; Lee, H.J.; Miyachi, H.; Watanabe, B.; Muranaka, T.; Saito, K.; Sugimoto, Y.; et al. Characterization of C-26 aminotransferase, indispensable for steroidal glycoalkaloid biosynthesis. Plant J. 2021. [Google Scholar] [CrossRef]
  27. Yang, L.; Cao, H.; Zhang, X.; Gui, L.; Chen, Q.; Qian, G.; Xiao, J.; Li, Z. Genome-Wide Identification and Expression Analysis of Tomato ADK Gene Family during Development and Stress. Int. J. Mol. Sci. 2021, 22, 7708. [Google Scholar] [CrossRef] [PubMed]
  28. Xu, D.-B.; Ma, Y.-N.; Qin, T.-F.; Tang, W.-L.; Qi, X.-W.; Wang, X.; Liu, R.-C.; Fang, H.-L.; Chen, Z.-Q.; Liang, C.-Y.; et al. Transcriptome-Wide Identification and Characterization of the JAZ Gene Family in Mentha canadensis L. Int. J. Mol. Sci. 2021, 22, 8859. [Google Scholar] [CrossRef]
  29. Zhang, X.; Lin, K.; Li, Y. Highlights to phytosterols accumulation and equilibrium in plants: Biosynthetic pathway and feedback regulation. Plant Physiol. Biochem. 2020, 155, 637–649. [Google Scholar] [CrossRef]
  30. Li, J.; Zhou, H.; Pan, X.; Li, Z.; Lu, Y.; He, N.; Meng, T.; Yao, C.; Chen, C.; Ling, X. The role of fluconazole in the regulation of fatty acid and unsaponifiable matter biosynthesis in Schizochytrium sp. MYA 1381. BMC Microbiol. 2019, 19, 256. [Google Scholar] [CrossRef]
  31. Song, N.; Cheng, Y.; Peng, W.; Peng, E.; Zhao, Z.; Liu, T.; Yi, T.; Dai, L.; Wang, B.; Hong, Y. Genome-Wide Characterization and Expression Analysis of the SBP-Box Gene Family in Sweet Orange (Citrus sinensis). Int. J. Mol. Sci. 2021, 22, 8918. [Google Scholar] [CrossRef]
  32. Wang, Y.; Zhang, Y.; Liu, Q.; Zhang, T.; Chong, X.; Yuan, H. Genome-Wide Identification and Expression Analysis of BBX Transcription Factors in Iris germanica L. Int. J. Mol. Sci. 2021, 22, 8793. [Google Scholar] [CrossRef]
  33. Maes, L.; Van Nieuwerburgh, F.C.W.; Zhang, Y.S.; Reed, D.W.; Pollier, J.; Casteele, S.R.F.V.; Inze, D.; Covello, P.S.; Deforce, D.L.D.; Goossens, A. Dissection of the phytohormonal regulation of trichome formation and biosynthesis of the antimalarial compound artemisinin in Artemisia annua plant. New Phytol. 2011, 189, 176–189. [Google Scholar] [CrossRef] [Green Version]
  34. Stitz, M.; Gase, K.; Baldwin, I.T.; Gaquerel, E. Ectopic Expression of AtJMT in Nicotiana attenuata: Creating a Metabolic Sink Has Tissue-Specific Consequences for the Jasmonate Metabolic Network and Silences Downstream Gene Expression. Plant Physiol. 2011, 157, 341–354. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, X.; Jin, M.; Tadesse, N.; Dang, J.; Zhou, T.; Zhang, H.; Wang, S.; Guo, Z.; Ito, Y. Dioscorea zingiberensis C. H. Wright: An overview on its traditional use, phytochemistry, pharmacology, clinical applications, quality control, and toxicity. J. Ethnopharmacol. 2018, 220, 283–293. [Google Scholar] [CrossRef]
  36. Jiao, Y.; Lau, O.S.; Deng, X.W. Light-regulated transcriptional networks in higher plants. Nat. Rev. Genet. 2007, 8, 217–230. [Google Scholar] [CrossRef] [PubMed]
  37. Cheng, J.; Chen, J.; Liu, X.; Li, X.; Zhang, W.; Dai, Z.; Lu, L.; Zhou, X.; Cai, J.; Zhang, X.; et al. The origin and evolution of the diosgenin biosynthetic pathway in yam. Plant Commun. 2021, 2, 100079. [Google Scholar] [CrossRef] [PubMed]
  38. Vaidya, K.; Ghosh, A.; Kumar, V.; Chaudhary, S.; Srivastava, N.; Katudia, K.; Tiwari, T.; Chikara, S.K. De Novo Transcriptome Sequencing in Trigonella foenum-graecum L. to Identify Genes Involved in the Biosynthesis of Diosgenin. Plant Genome 2013, 6, 1–11. [Google Scholar] [CrossRef] [Green Version]
  39. Sonawane, P.D.; Pollier, J.; Panda, S.; Szymanski, J.; Massalha, H.; Yona, M.; Unger, T.; Malitsky, S.; Arendt, P.; Pauwels, L.; et al. Plant cholesterol biosynthetic pathway overlaps with phytosterol metabolism. Nat. Plants 2016, 3, 16205. [Google Scholar] [CrossRef] [PubMed]
  40. Su, Y.; Huang, Y.; Dong, X.; Wang, R.; Tang, M.; Cai, J.; Chen, J.; Zhang, X.; Nie, G. Exogenous Methyl Jasmonate Improves Heat Tolerance of Perennial Ryegrass Through Alteration of Osmotic Adjustment, Antioxidant Defense, and Expression of Jasmonic Acid-Responsive Genes. Front. Plant Sci. 2021, 12, 664519. [Google Scholar] [CrossRef]
  41. Paponov, M.; Antonyan, M.; Slimestad, R.; Paponov, I.A. Decoupling of Plant Growth and Accumulation of Biologically Active Compounds in Leaves, Roots, and Root Exudates of Hypericum perforatum L. by the Combination of Jasmonate and Far-Red Lighting. Biomolecules 2021, 11, 1283. [Google Scholar] [CrossRef] [PubMed]
  42. Paquette, S.M.; Bak, S.; Feyereisen, R. Intron–Exon Organization and Phylogeny in a Large Superfamily, the Paralogous Cytochrome P450 Genes of Arabidopsis thaliana. DNA Cell Biol. 2000, 19, 307–317. [Google Scholar] [CrossRef] [PubMed]
  43. Yu, C.S.; Lin, C.J.; Hwang, J.K. Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004, 13, 1402–1406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hu, B.; Jin, J.; Guo, A.-Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2014, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, w202–w208. [Google Scholar] [CrossRef]
  46. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  47. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
Ijms 22 10953 g001 550
Figure 1. Phylogenetic tree of DzCYP72A proteins.
Figure 1. Phylogenetic tree of DzCYP72A proteins.
Ijms 22 10953 g001
Ijms 22 10953 g002 550
Figure 2. Multiple alignment of DzCYP72A proteins.
Figure 2. Multiple alignment of DzCYP72A proteins.
Ijms 22 10953 g002
Ijms 22 10953 g003 550
Figure 3. Gene structure of DzCYP72As.
Figure 3. Gene structure of DzCYP72As.
Ijms 22 10953 g003
Ijms 22 10953 g004 550
Figure 4. Motif pattern of DzCYP72As. (A) The motif composition of DzCYP72A proteins predicted by MEME and plotted in TBtools. The detailed information of ten motifs is exhibited by different colored boxes. The length of every DzCYP72A protein can be measured by the scale at the bottom. (B) Detailed information of ten motifs.
Figure 4. Motif pattern of DzCYP72As. (A) The motif composition of DzCYP72A proteins predicted by MEME and plotted in TBtools. The detailed information of ten motifs is exhibited by different colored boxes. The length of every DzCYP72A protein can be measured by the scale at the bottom. (B) Detailed information of ten motifs.
Ijms 22 10953 g004
Ijms 22 10953 g005 550
Figure 5. Predicted cis-regulatory elements in promoter regions of DzCYPP72As.
Figure 5. Predicted cis-regulatory elements in promoter regions of DzCYPP72As.
Ijms 22 10953 g005
Ijms 22 10953 g006 550
Figure 6. Predicted cis-regulatory elements in promoter regions of DzCYPP72As. (A) Tissue-specific expression profile of DzCYP72As. (B) Expression patterns of DzCYP72As under jasmonic acid treatment. S1: 0 μmol/L; S2: 25 μmol/L; S3: 50 μmol/L; S4: 100 μmol/L; S5: 200 μmol/L. The bar in the upper right corner represents the expression value of log10.
Figure 6. Predicted cis-regulatory elements in promoter regions of DzCYPP72As. (A) Tissue-specific expression profile of DzCYP72As. (B) Expression patterns of DzCYP72As under jasmonic acid treatment. S1: 0 μmol/L; S2: 25 μmol/L; S3: 50 μmol/L; S4: 100 μmol/L; S5: 200 μmol/L. The bar in the upper right corner represents the expression value of log10.
Ijms 22 10953 g006
Ijms 22 10953 g007 550
Figure 7. Spearman analysis of correlation between specialized metabolites and DzCYP72As.
Figure 7. Spearman analysis of correlation between specialized metabolites and DzCYP72As.
Ijms 22 10953 g007
Table
Table 1. Characters of DzCYP72A proteins identified in D. zingiberensis.
Table 1. Characters of DzCYP72A proteins identified in D. zingiberensis.
NameAmino
Acids		MW (kD)pIInstability
Index		Aliphatic
Index		GRAVYGenomic Location
DzCYP72A151859.53 9.42 50.42 92.08 −0.207 Chr1:54408546–54412582
DzCYP72A252059.52 9.09 46.29 93.19 −0.169 Chr1:54419523–54423635
DzCYP72A351658.89 9.01 45.88 92.60 −0.192 Chr1:54482572–54486042
DzCYP72A451959.74 8.98 42.74 92.02 −0.182 Chr1:54491399–54496848
DzCYP72A551959.40 9.23 49.52 94.68 −0.161 Chr1:54620742–54635835
DzCYP72A649956.98 8.77 48.86 91.44 −0.208 Chr1:54694341–54697800
DzCYP72A748655.91 9.06 47.04 90.86 −0.241 Chr1:54703122–54708457
DzCYP72A851859.46 9.16 47.69 92.82 −0.186 Chr1:54817699–54821732
DzCYP72A949356.52 9.03 44.79 92.35 −0.176 Chr1:54828984–54833020
DzCYP72A1051859.28 9.04 46.59 93.92 −0.200 Chr1:54911925–54915459
DzCYP72A1151759.57 9.10 43.23 93.52 −0.162 Chr1:54920693–54924967
DzCYP72A1251859.21 9.27 47.44 95.98 −0.142 Chr1:54948190–54951851
DzCYP72A1351859.29 9.19 46.19 95.98 −0.186 Chr1:54952667–54955365
DzCYP72A1452560.06 9.44 46.68 93.79 −0.204 Chr1:55018259–55022509
DzCYP72A1551959.51 8.23 37.99 97.86 −0.093 Chr4:20462820–20465388
DzCYP72A1651559.04 8.36 45.35 94.45 −0.167 Chr9:1757458–1760130
DzCYP72A1748655.77 8.12 50.32 92.49 −0.271 Chr9:1797953–1806784
DzCYP72A1845151.87 8.40 46.16 91.02 −0.292 Chr9:1957935–1963395
DzCYP72A1917319.63 6.29 42.93 94.68 −0.070 Chr9:1963533–1964154
DzCYP72A2050357.97 9.13 47.24 95.37 −0.203 Chr9:2174147–2176010
DzCYP72A2148756.05 7.27 51.27 93.31 −0.244 Chr9:2181213–2183031
DzCYP72A2248655.86 8.58 50.57 91.30 −0.290 Chr9:2219226–2221041
DzCYP72A2351859.07 9.22 46.65 96.74 −0.121 Unknown
DzCYP72A2451859.36 9.11 46.56 95.23 −0.207 Unknown
DzCYP72A2552560.09 9.41 46.32 93.79 −0.205 Unknown
MW: molecular weight; pI: isoelectric point; GRAVY: grand average of hydropathicity; Chr: chromosome.
Table
Table 2. Predicted subcellular localization of CYP72A proteins in D. zingiberensis.
Table 2. Predicted subcellular localization of CYP72A proteins in D. zingiberensis.
Gene NamePlant-PLocWoLF PSORTCELLOYLoc
CYP72A1chloroplastchlo: 5mito(2.82)chloroplast
CYP72A2mitochondrionchlo: 6mito(2.38)chloroplast
CYP72A3chloroplastchlo: 6mito(2.07)chloroplast
CYP72A4cytoplasmplas: 7mito(2.06)chloroplast
CYP72A5chloroplastchlo: 6mito(2.76)chloroplast
CYP72A6chloroplastchlo: 5mito(1.7)/cyto(1.33)chloroplast
CYP72A7chloroplastcyto: 9mito(2.19)nucleus
CYP72A8mitochondrionchlo: 6mito(2.56)chloroplast
CYP72A9mitochondrionchlo: 7mito(2.08)/cyto(1.39)chloroplast
CYP72A10chloroplastchlo: 8mito(2.15)chloroplast
CYP72A11cytoplasmchlo: 8mito(2.07)chloroplast
CYP72A12cytoplasmchlo: 8mito(2.16)chloroplast
CYP72A13cytoplasmchlo: 5mito(2.59)chloroplast
CYP72A14cytoplasmE.R.: 5mito(2.97)chloroplast
CYP72A15chloroplastchlo: 9mito(1.68)/cyto(1.67)chloroplast
CYP72A16chloroplastchlo: 5cyto(1.71)chloroplast
CYP72A17chloroplastcyto: 7cyto(1.70)/mito(1.21)/nucl(1.14)nucleus
CYP72A18chloroplastcyto: 9cyto(3.09)cytoplasm
CYP72A19cytoplasmcyto: 11cyto(2.26)cytoplasm
CYP72A20chloroplastcyto:11mito(1.99)/cyto(1.42)secreted pathway
CYP72A21chloroplastcyto:7cyto(1.97)chloroplast
CYP72A22chloroplastcyto: 7cyto(1.40)/nucl(1.35)/mito(1.30)nucleus
CYP72A23chloroplastchlo: 5mito(2.21)chloroplast
CYP72A24chloroplastchlo: 5mito(2.63)chloroplast
CYP72A25cytoplasmE.R.: 5mito(2.83)chloroplast
chlo: chloroplast; plas: plastid; E.R.: endoplasmic reticulum; cyto: cytoplasm; mito: mitochondrion; nucl: nucleus.
Table
Table 3. Effects of jasmonic acid on steroidal saponin biosynthesis in D. zingiberensis.
Table 3. Effects of jasmonic acid on steroidal saponin biosynthesis in D. zingiberensis.
Compound Name0 μmol/L25 μmol/L50 μmol/L100 μmol/L200 μmol/L
Parvifloside (mg/g)34.88 ± 0.85 c39.26 ± 1.08 b24.83 ± 1.25 d54.89 ± 0.88 a34.20 ± 0.55 c
Protodeltonin (mg/g)3.81 ± 0.21 d10.75 ± 1.51 bc15.69 ± 1.68 a12.37 ± 0.11 b10.38 ± 0.25 c
Dioscin (µg/g)19.35 ± 1.00 c23.79 ± 0.84 b38.25 ± 1.07 a23.20 ± 1.07 b11.44 ± 0.23 d
Diosgenin (µg/g)404.30 ± 19.97 d434.24 ± 6.37 d595.27 ± 20.86 b647.18 ± 24.21 a478.64 ± 22.59 c
Campesterol (µg/g)78.33 ± 10.92 a76.12 ± 2.06 ab80.53 ± 19.24 a88.13 ± 21.41 a51.1 ± 1.17 b
Stigmasterol (µg/g)128.91 ± 3.92 b117.33 ± 4.64 ab253.94 ± 16.34 a68.93 ± 6.93 c69.68 ± 1.66 c
β-Sitosterol (µg/g)278.62 ± 12.23 bc299.21 ± 16.38 ab281.27 ± 17.31 c322.62 ± 14.77 a226.12 ± 12.55 d
Note: The different lowercase letters represent significant difference at p < 0.05 (ANOVA with least-significant difference).
Table
Table 4. Redundancy analysis of steroidal biosynthesis and DzCYP72As in D. zingiberensis.
Table 4. Redundancy analysis of steroidal biosynthesis and DzCYP72As in D. zingiberensis.
StatisticsRDA1RDA2FPr(>F)
DzCYP72A10.9925 −0.1219 2.6747 0.0500 *
DzCYP72A30.9799 0.1996 2.6286 0.0417 *
DzCYP72A60.9486 0.3164 2.9872 0.0083 **
DzCYP72A90.9873 −0.1586 3.0657 0.0417 *
DzCYP72A110.9947 −0.1026 3.7190 0.0417 *
DzCYP72A130.9960 0.0896 2.5909 0.0667
DzCYP72A140.9838 −0.1794 4.1352 0.0167 *
DzCYP72A160.9993 −0.0386 4.9162 0.0083 **
DzCYP72A170.9997 0.0241 4.3933 0.0083 **
DzCYP72A200.9812 −0.1931 3.8978 0.0167 *
DzCYP72A210.9313 −0.3643 1.9759 0.1333
DzCYP72A220.9250 −0.3800 2.4202 0.1583
Note: * represents significant; ** represents highly significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hou, L.; Yuan, X.; Li, S.; Li, Y.; Li, Z.; Li, J. Genome-Wide Identification of CYP72A Gene Family and Expression Patterns Related to Jasmonic Acid Treatment and Steroidal Saponin Accumulation in Dioscorea zingiberensis. Int. J. Mol. Sci. 2021, 22, 10953. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222010953

AMA Style

Hou L, Yuan X, Li S, Li Y, Li Z, Li J. Genome-Wide Identification of CYP72A Gene Family and Expression Patterns Related to Jasmonic Acid Treatment and Steroidal Saponin Accumulation in Dioscorea zingiberensis. International Journal of Molecular Sciences. 2021; 22(20):10953. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222010953

Chicago/Turabian Style

Hou, Lixiu, Xincheng Yuan, Song Li, Yi Li, Zihao Li, and Jiaru Li. 2021. "Genome-Wide Identification of CYP72A Gene Family and Expression Patterns Related to Jasmonic Acid Treatment and Steroidal Saponin Accumulation in Dioscorea zingiberensis" International Journal of Molecular Sciences 22, no. 20: 10953. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms222010953

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop