Selection and validation of reference genes desirable for gene expression analysis by qRT-PCR in MeJA-treated ginseng hairy roots
Authors:
Li Li aff001; Kangyu Wang aff001; Mingzhu Zhao aff001; Shaokun Li aff001; Yue Jiang aff001; Lei Zhu aff001; Jing Chen aff001; Yanfang Wang aff002; Chunyu Sun aff001; Ping Chen aff001; Jun Lei aff001; Meiping Zhang aff001; Yi Wang aff001
Authors place of work:
College of Life Science, Jilin Agricultural University, Changchun, Jilin, China
aff001; Research Center of Ginseng Genetic Resources Development and Utilization, Changchun, Jilin, China
aff002; College of Chinese Medicinal Materials, Jilin Agricultural University, Changchun, Jilin, China
aff003; Key Laboratory of Natural Resources of Changbai Mountain & Functional Molecules (Yanbian University), Ministry of Education, Yanji, Jilin, China
aff004
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0226168
Summary
Ginseng is a valuable herb of traditional Chinese medicine and ginsenosides, the main bioactive components of ginseng, have been proven to have multiple functions in human therapies and health. Methyl jasmonate (MeJA) is an elicitor that has been demonstrated to have a vital influence on ginsenoside biosynthesis. Quantitative real-time polymerase chain reaction (qRT-PCR) has been widely used in quantification of gene expressions. Here, we report the selection and validation of reference genes desirable for normalization of gene expressions quantified by qRT-PCR in ginseng hairy roots treated with MeJA. Twelve reference genes were selected as candidate genes, and their expressions were quantified by qRT-PCR, and analyzed by geNorm, NormFinder and BestKeeper. CYP and EF-1α were shown to be the most stable reference genes in geNorm, CYP was the most stable reference gene in NormFinder, and 18S was the most stable reference gene in BestKeeper. On this basis, we further quantified the relative expression levels of four genes encoding key enzymes that are involved in ginsenoside biosynthesis using CYP and 18S as the reference genes, respectively. Moreover, correlation analysis was performed between the quantified expressions of four genes and the ginsenoside content in MeJA-treated ginseng hairy roots. The results of relative expressions of the four genes quantified using CYP as the reference gene and their significant correlations with the ginsenoside content were better than those using 18S as the reference gene. The CYP gene, hence, was concluded as the most desirable reference gene for quantification of the expressions of genes in MeJA-treated ginseng hairy roots. This finding, therefore, provides information useful for gene research in ginseng, particularly in MeJA-treated ginseng hairy roots, which includes identification and characterization of genes involved in ginsenoside biosynthesis.
Keywords:
Gene expression – Hair – Biosynthesis – Polymerase chain reaction – RNA extraction – Seedlings – Glycosides – Secondary metabolism
1. Introduction
Ginseng (Panax ginseng C.A. Meyer) is a valuable herb of traditional Chinese medicine that has curative effects on several human diseases, including cancerous[1], inflammatory[2], cardiovascular[3] and neurodegenerative diseases[4]. It also has effects on reducing obesity and the regulation of immunity[5]. Nevertheless, ginseng research is facing several challenges. For instance, ginseng must grow at least four years before it could be effectively used in medicine, during which it is likely subjected to a variety of diseases. Wild ginseng is rare and thus, is extremely valuable. Therefore, several tissue culture methods, such as callus, somatic embryogenesis, adventitious roots and hairy roots culture[6] have been developed and used to generate plant materials for ginseng research. Among the plant materials generated through these methods, hairy roots are desirable for numerous aspects of ginseng research since they contain components that are very similar to native ginseng, grow rapidly and are hormone-autotrophic. Ginsenosides are the main bioactive components in ginseng, including protopanaxadiol (PPD)-type, protopanaxatriol (PPT)-type, oleanolic acid-type and ocotillo-type[7]. As the secondary metabolites, ginsenosides also play a vital role in plant defense to biotic and abiotic stresses. More than 150 naturally occurring ginsenosides have been so far identified from Panax species[8], but their types and content vary greatly among tissues and through periods of growth and development. It was reported that the ginsenoside content was the highest in leaves and the lowest in primary roots. The content ratio of PPT-type and PPD-type ginsenoside was about 1.0 in primary roots, while it was 1.37–3.14 in leaves[9]. Moreover, the tissue-specific content and distributions of ginsenosides could be changed by elicitors through controlling growth conditions, for which the elicitors used for such purposes including fungi[10], N, N-dicyclohexylcarbodiimide (DCCD)[11], salicylic acid (SA)[12] and Methyl jasmonate (MeJA)[13]. Among these elicitors, MeJA has been proven to be the most effective and therefore, been used most widely. It has been reported that MeJA has a vital influence on the expression of the genes involved in ginsenoside biosynthesis and ginsenoside content in hairy roots[14, 15], adventitious roots[16, 17] and other tissues[18]. Therefore, these elicitors can be potentially used for increasing ginsenoside production, identifying genes involved in ginsenoside biosynthesis, as well as studying the molecular mechanisms of ginsenoside biosynthesis and regulation. Especially, the MeJA-treated ginseng hairy roots have been widely used in the study of ginsenoside biosynthesis.
Quantitative real-time polymerase chain reaction (qRT-PCR) has been widely used for quantification of gene expressions[19]. The reference gene used for gene expression quantification by qRT-PCR is a vital factor to normalize gene differential expression among samples, but its expression often substantially varies with tissues, treatments and growth conditions[20]. Liu et al.[21] investigated 20 genes as the reference genes for qRT-PCR at different growth stages and in different tissues for ginseng research. Of these, 20 candidate genes studied, 10 were traditional housekeeping genes and 10 were novel genes selected by transcriptome analysis. Later, Wang et al.[22] also conducted a similar study in different tissues of ginseng and ginseng seedlings treated with heat stress. Nevertheless, reference genes that were desirable for ginseng research using hairy roots grown under the stress or elicited condition, especially those treated with MeJA, have not been reported.
Several sorts of software have been used to test the stability of reference genes for gene expression analysis by qRT-PCR, including geNorm, NormFinder, RefFinder and Bestkeeper. They all calculate the stability of reference genes based on the Ct value. Recently, a validation step has been included in the process, following the Ct value calculation[23–28]. Unfortunately, only genes that are sensitive to treatment, such as HSF (heat shock protein), were used for reference gene validation, and the reference gene was believed to be suitable when its expression level significantly changed in different conditions. However, it is essential and desirable to further validate the reference gene by examining the expressions of a selection of targeted genes and the correlation of their expression with traits.
In the present study, we treated the ginseng hairy roots with MeJA, determined the change of its ginsenoside content, selected the reference genes for gene expression analysis by qRT-PCR, and validated the desirable reference gene by correlation analysis between ginsenoside content and the expressions of genes involved in ginsenosides biosynthesis. The reference genes selected and validated in this study, therefore, will facilitate the identification and characterization of genes involved in ginsenoside biosynthesis.
2. Materials and methods
2.1 Plant material and MeJA treatment
Ginseng cultivar ‘Damaya’ was used as plant materials for hairy root inducing in this study. After washed and sterilized, the ginseng seeds were germinated on 1/2 MS medium and cultured at 23°C for 15 days till the seedlings leaf fully unfolded. The seedlings were cut into 0.5 cm segments and dipped into Agrobacterium rhizogenes strain A4 for 15 min. After infection, the explants were dried on sterilized filter papers, and cultured at 23°C in dark on 1/2 MS medium containing 200 mg/L acetosyringone. After 2 days of co-cultivation, the explants were transferred to 1/2 MS medium containing 50 mg/L kanamycin and 500 mg/L cefotaxime to eliminate the Agrobacterium. Putative hairy roots were obtained approximately in a month after selection (50 mg/L kanamycin). Every induced single hairy root was individually further cultured to ensure the homogeneity of genetic background. The fastest-growing hairy root line was used as the materials for MeJA treatment. The ginseng hairy roots for MeJA treatment were cultured in flasks containing 150 ml 1/2 MS liquid medium with an initial inoculation of 1.0 g hairy roots at 22 oC with shaking at 110 rpm. After cultured for 23 days, the hairy roots were treated with 200 μM MeJA that was previously shown to be the most desirable for the MeJA treatment in our laboratory. The hairy roots were harvested after treatment at 6 h, 12 h, 24 h, 48 h, 72 h, 96 h and 120 h, respectively. Three biologically duplicated samples and one blank control with no MeJA treatment were harvested at every time point. A part of hairy root samples was flash-frozen in liquid nitrogen for RNA extraction and the remaining part was dried to constant weight for ginsenoside extraction.
2.2 Ginsenoside extraction and content quantification
Ginsenosides were extracted by the Soxhlet extraction method[29]. Weighed samples (0.5 g) were wrapped up with filter paper and transferred into a 250 ml refluxing-type Soxhlet extractor. The samples were extracted by 100 ml methanol in 90 oC water bath for 12 h, with every sample having three technical replicates. The extracts were collected and steamed to dry, and then the residue was dissolved by 10 ml chromatographic methanol, and filtered for analysis.
The mono-ginsenosides in each sample were separated using the Waters Alliance HPLC, with e2695 Separations Module, and their content was determined using the Waters 2489 Ultraviolet Spectrophotometric Detector (Waters, Milford, MA, USA). A Waters Uxbridge C18 column (4.6 mm × 250 mm, 5 μm) was used to achieve the separation. The mobile phase containing Solvent A (acetonitrile) and Solvent B (water) was set to a gradient elution program as reported by Cong et al.[30]: 0–40 min (18–21% A), 40–42 min (21–26% A), 42–46 min (26–32% A), 46–66 min (32–33.5% A), 66–71 min (33.5–38% A), 71–86 min (38–65% A), 86–91 min (65% A), 91–96 min (65–85% A), 96–103 min (85% A), 103–105 min (85–18% A) and 105–106 min (18% A). The flow rate was 1.0 ml/min, the sample injection volume was 20 μL, the column oven temperature was 35 oC, and the detection wavelength was 203 nm. The content of total ginsenoside was examined by the Vanillin colorimetric method as reported[31].
2.3 Total RNA isolation and cDNA synthesis
Total RNA of ginseng hairy roots was isolated using the TRIpure Reagent Total RNA Extraction Reagent (Bioteke, Beijing, China). The concentration and quality of the RNA were determined by Scandrop 100 (Analytic Jena AG). cDNA synthesis from the RNA was carried out using the PrimeScript™ RT Reagent Kit with gDNA Eraser (Perfect Real Time) (TaKaRa, Dalian, China). The concentration of cDNA was also determined by Scandrop 100 and then diluted to 100 ng/μL for qRT-PCR.
2.4 Candidate reference genes and key enzyme coding genes of ginsenoside biosynthesis for qRT-PCR
Based on the previous studies in ginseng[21, 22], 12 widely used qRT-PCR reference genes were selected as candidate genes, including 18S (18S ribosomal RNA), ACT1 (actin 1), aTUB (tubulin alpha-1 chain), bTUB (beta-tubulin), CYC (cyclophilin ABH-like protein), CYP (cyclophilin), EF-1α (elongation factor 1-alpha), eIF-5A (translational initiation factor eIF-5a), F-box (F-box containing protein), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), IF3G1 (eukaryotic translation initiation factor 3G1), and UBQ (polyubiquitin). The primers of these genes reported by these researchers were used in this study. Four already verified key enzyme genes in mevalonate pathway of ginsenoside biosynthesis were randomly selected as the positive controls for validation of the reference genes by qRT-PCR, including PgFPS[32], (DQ087959.1, F-GGATGATTATCTGGATTGCTTTGG, R-CAGTGCTTTTACTACCAACCAGGAG) PgDDS[33] (AB122080.1, F- CGGAACGATTGACACTATTCTGAC, R- CTGACCCAATCATCGTGCTGT), PgUGT71A27 (UGTPg1)[34] (KM491309.1, F- TGCGTCCGTCTATCCCTAAAG, R- TGATGTCCTGTCCAAGAATCCTAC) and PgCYP716A47 [35] (JN604536.1, F- TTAGGTGATACGGCGGCAG, R-CTGGGGGATGCGTTTTGTAT).
2.5 qRT-PCR
qRT-PCR was carried out with the ABI 7500 Real-Time PCR System (Applied Biosystem, Foster City, CA, USA) and TB Green™ Premix Ex Taq™ (Tli RNase H Plus) (TaKaRa, Dalian, China). The reaction system included 10 μL TB Green™ Premix Ex Taq™ (2X), 0.4 μL Rox Ⅱ, 0.8 μL each primer, 1 μL cDNA and 7 μL RNase-free water. The reaction was performed as follows: pre-denaturation at 95 oC for 30 s; PCR of 40 cycles at 95 oC for 5 s and 60 oC for 34 s; melting curve at 95 oC for 15 s, 60 oC for 1 min, and 95 oC for 15 s. Each biological sample was tested in triplicate. The gene relative expression levels were calculated by the 2-ΔΔCt method.
2.6 Gene expression stability analysis
To validate the expression stability of candidate reference genes, the results of qRT-PCR were analyzed by the geNorm[36], NormFinder[37], and BestKeeper[38] methods, respectively, according to the manufacturers’ instructions. In geNorm and NormFinder, the input data were the gene relative quantities that transferred from the raw Ct results of qRT-PCR, and the output data were the stability value of gene expression (M value). In BestKeeper, the input data were the untransformed Ct results of qRT-PCR, and the output data were a series of results that indicate the expression stability of the genes.
2.7 Correlation analysis between gene expressions and ginsenoside content
The relative expression levels of genes involved in ginsenoside biosynthesis normalized by the candidate reference genes and the mono-ginsenoside content in the same sample were used for the correlation analysis. The Spearman's correlation coefficient was calculated by SPSS (version 23) software to reflect the relationship between gene expressions and ginsenoside content.
3. Results
3.1. Changes of ginsenoside content in the ginseng hairy roots after the MeJA treatment
The content of 14 mono-ginsenosides in the MeJA-treated hairy roots were determined. These 14 mono-ginsenosides included Rb1, Rb2, Rb3, Rc, Rd, Re, Rf, Rg1, Rg2, Rh1, Rh2, F1, F2 and CK. We also quantified the content of glycoside PPT, glycoside PPD and total ginsenoside in the MeJA-treated hairy roots. The content of mono-ginsenoside and total ginsenoside changing with MeJA-treated time were showed as Fig 1. The content of 13 mono-ginsenosides, two glycosides and total ginsenoside were significantly changed with MeJA treatment (CK have not been detected). The data of the control group (0 h) with no MeJA treatment was the average content of blank samples collected at seven-time points after MeJA treatment because the mono-ginsenoside content little changed through the 120-hour period for the untreated samples.
The changing trend of ginsenoside content showed two patterns. Nine of the 13 mono-ginsenosides and total-ginsenoside content showed constantly and significantly increasing trend, while the remaining four mono-ginsenosides and glycoside content fluctuated along with the MeJA treatment time. The constant increase of the mono-ginsenoside content indicated the accumulation of their biosynthesis while the fluctuation of the mono-ginsenoside content might reflect the dynamic status of their synthesis and decomposition as the time of MeJA treatment lasted.
The time points when the significant change of ginsenoside content initially occurred varied from 12 hours to 96 hours. The F1 mono-ginsenoside content significantly increased at 12 h after the treatment. The Re, Rg1, Rh1, Rb1, Rb2, Rc, Rd, and F2 content significantly increased at 24 h after the treatment. The Rb3, Rh2, PPD and total ginsenoside content significantly increased at 48 h after the treatment. The Rf content significantly increased at 72 h after the treatment and the Rg2 content significantly increased at 96 h after the treatment. The time points from which the significant increase of ginsenoside content in the MeJA-treated hairy roots started might suggest the positions and the order of the mono-ginsenoside in the ginsenoside biosynthesis pathway.
Moreover, we summed the content of PPT-type and PPD-type ginsenosides respectively, and analyzed their content ratio at the different time point of the MeJA treatment (Table 1). Before MeJA-treated, the content ratio of PPT-type ginsenosides to PPD-type ginsenosides in ginseng hairy roots was 1.28, but it increased along with the MeJA treatment and reached a maximum of 1.50 at 48 h after treatment. This result showed that the biosynthesis and accumulation of PPT-type ginsenosides were faster than PPD-type in the MeJA-treated ginseng hairy roots.
In addition, the content of F1 in ginseng hairy roots was extremely high after the MeJA treatment. It was 0.599 mg/g (0.06%) in ginseng hairy roots before the MeJA treatment and reached 7.298 mg/g (0.73%) after 120-hour MeJA treatment that was similar to the content of Rb1. This was a new finding that the F1 content could increase to the same level of Rb1 content in ginseng hairy roots under MeJA treatment. This result provides information useful for the production of ginsenoside F1, even though additional research is needed to further confirm this result.
3.2 Selection of reference genes desirable for qRT-PCR in the MeJA-treated ginseng hairy roots
3.2.1 Expression profiles of the candidate reference genes
The expression stabilities of all 12 candidate reference genes were preliminarily displayed by cycle threshold (Ct) values that represent the transcript abundance of the genes in the tested samples. In all samples, the Ct values of the candidate reference genes ranged from 15 to 31 (Fig 2). Among which, the 18S (15.73±0.80) was the most abundant in the MeJA-treated hairy roots, followed by those of CYC (21.24±1.22), 5S (21.98±0.67), ACT1 (24.44±1.07), GAPDH (25.02±1.27), CYP (25.94±1.19), EF-1α (26.70±1.23), IF3G1 (26.70±1.20), bTUB (26.73±1.19), UBQ (27.20±1.24), F-box (27.33±1.26), aTUB (28.02±1.17) and eIF-5A (31.13±1.26).
3.2.2 Stability analysis by geNorm
The pairwise variation and the M value are the two main indicators of reference gene expression stability in geNorm, which represent the suitable gene number for normalization and the most desirable reference gene. Fig 3 showed the analysis results of the 12 candidate reference genes in MeJA-treated hairy roots using geNorm. The pairwise variation values of all gene pairs were smaller than the cut-off value (0.15), suggesting that two genes should be used for qRT-PCR normalization in the hairy roots treated by MeJA. The M values of these genes in descending order were 18S, eIF-5A, IF3G1, CYC, GAPDH, ACT1, bTUB, aTUB, UBQ, F-box, CYP and EF-1α. Since the most suitable reference genes for qRT-PCR normalization should have the minimum M value, the CYP and EF-1α genes were selected for qRT-PCR normalization for the hairy roots treated with MeJA.
3.2.3 Stability analysis by NormFinder
The stability value based on variance analysis was the only indicator of reference gene expression stability in NormFinder. As the M value of geNorm, the smaller the stability value, the more stable the gene and the more suitable the gene for qRT-PCR normalization. The analysis results of the 12 candidate reference genes were shown in Table 2. The most stable reference gene was shown to be CYP, with the value of 0.234; the most unstable genes were shown to be 18S (0.909) and eIF-5A (0.858), and the M values of the remaining genes were between 0.2 and 0.4. This result was in line with the results of geNorm.
3.2.4 Stability analysis by BestKeeper
The standard deviation (SD), coefficient variation (CV) and correlation coefficients are the vital indicators of gene expression stability in BestKeeper. The result of the BestKeeper analysis was shown in Table 3. Only the gene with an SD value below 1.0 was desirable for qRT-PCR, so none of the 12 candidate reference genes met this criterion. Nevertheless, the correlation coefficients of all the 12 genes were quite high and the p-values were all < 0.001. Among the 12 candidate genes, the SD value of 18S was the lowest (SD = 1.25).
3.3 Validation of the most desirable reference genes
3.3.1 The relative expression levels of key enzyme genes involved in ginsenoside biosynthesis normalized by different reference genes
Because the ginseng hairy roots are widely used in the study of ginsenoside biosynthesis, we chose four genes encoding key enzymes in the ginsenoside biosynthesis pathway as the target genes to validate the most desirable reference genes obtained above. In the results of geNorm and NormFinder, the gene of CYP was the most stable reference gene and the 18S was the most unstable gene. However, in the result of BestKeeper, the 18S was the most stable reference gene. Therefore, the CYP and 18S genes were selected as the reference genes to confirm whether they were suitable for normalization of gene expression quantified by qRT-PCR in the MeJA-treated hairy roots.
The qRT-PCR was conducted for the four selected target genes with CYP and 18S as reference genes. The relative expression levels of these four genes were calculated using CYP and 18S as the reference genes respectively (Fig 4). When CYP was used as the reference gene, the relative expression levels of all four genes significantly increased compared with the no-treated hairy roots. The relative expression levels of PgFPS and PgCYP716A47 reached a maximum at 6 h after the MeJA treatment and then, gradually reduced along with the time of MeJA treatment. The relative expression levels of PgDDS and PgUGT71A27 reached their maximum at 48 h and 24 h, respectively, after the MeJA treatment and then gradually decreased as the time of the MeJA treatment lasted. Besides, the trends of these four genes showed good regularity. When 18S was used as the reference gene, the relative expression levels of all four genes significantly increased only partially, and absent in regularity. Therefore, CYP was more desirable as the reference gene for gene expression analysis in ginseng hairy roots treated with MeJA.
3.3.2 The correlation analysis of gene expression levels with ginsenoside content
To further confirm the above observation, correlation analysis was performed between the relative expressions of four key enzyme genes involved in ginsenoside biosynthesis and ginsenoside content in the MeJA-treated ginseng hairy roots (Table 4). When CYP was used as the reference gene, the expressions of PgFPS and PgUGT71A27 were significantly correlated with the content of 11 kinds of ginsenosides; PgDDS and PgCYP716A47 were significantly correlated with the content of 8 kinds of mono-ginsenosides. When 18S was used as the reference gene, the expression of PgFPS was significantly correlated with the content of 8 kinds of mono-ginsenosides; PgDDS was significantly correlated with the content of only 1 kind of mono-ginsenoside; PgUGT71A27 was significantly correlated with the content of 9 kinds of mono-ginsenosides; PgCYP716A47 was not significantly correlated with the content of any kinds of mono-ginsenoside analyzed. This result further verified that CYP was more desirable as a reference gene than 18S for gene expression analysis by qRT-PCR in the MeJA-treated ginseng hairy roots.
4. Discussion
4.1 The changes of different type mono-ginsenoside content under MeJA treatment
Methyl jasmonate as the elicitor of secondary metabolism can induce the accumulation of ginsenosides in ginseng. Kim et al.[39] reported that the ginsenosides of the Rb group accumulated more than that of the Rg group under the MeJA treatment. Kang et al. [40] and Oh et al.[18] reported that MeJA treatment could significantly enhance only PPD-type ginsenosides content. This study showed that the content of all the detected mono-ginsenosides and glycosides increased in hairy roots after the MeJA treatment. Moreover, we showed that the accumulation of the PPT-type ginsenosides was faster than that of the PPD-type ginsenosides under the MeJA treatment, which was different from the results reported before.
There may be two reasons that led to the difference among these studies. The one may be due to the different plant materials used for the studies and the other could be attributed to different mono-ginsenosides detected. The types and content of secondary metabolites as well as synthetic substrates and enzymes[41] may differ in different plant materials. Kim et al.[39] and Kang et al.[40] used the adventitious roots of ginseng for their studies, and Oh et al.[18] used fine root, root body, epidermis, rhizome, stems and leaf for their studies while we used hairy roots as materials in this study. However, we consider the types of detected mono-ginsenosides was the major reason for the observed discrepancy. The numbers of PPD-type mono-ginsenosides and PPT-type mono-ginsenosides detected by Oh et al.[18] were four and four, respectively. Those detected by Kang et al.[40] were sixteen and four, while ours were eight and seven. As more than 150 naturally occurring ginsenosides have been identified from Panax species, even if some of them account for a large proportion, we still cannot conclude that the MeJA treatment influences only certain types of ginsenosides. Moreover, the mechanism of MeJA in secondary metabolism is complicated. Several factors have been identified to play a vital role in the induction of ginsenoside biosynthesis, including JASMONATE ZIM DOMAIN (JAZ) proteins[42], transcription factors (TFs)[43], phytohormones[44] and some enzyme-coding genes[45]. All these factors were shown to have a global influence on secondary metabolisms, except for enzyme-coding genes. However, enzyme-coding genes may also have multiple functions that a single enzyme may catalyze biosynthesis of multiple secondary metabolites[8]. Hence, the influence of MeJA on secondary metabolism is likely to be overall, instead of metabolite-specific. Therefore, we should consider the induction of ginsenosides by MeJA from a global perspective.
4.2 The validation of CYP as the reference gene in ginseng hairy roots under MeJA treatment
By using gene expression quantified with qRT-PCR, followed by correlation analysis with the ginsenoside content variation, the genes involved in ginsenoside biosynthesis can be confirmed, at least preliminarily. Therefore, the stably expressed reference genes used for qRT-PCR have become a significant factor of properly quantifying gene expressions. In fact, several reference genes have previously been used in ginseng. The β-actin gene has been the most commonly used reference gene for all types of ginseng materials, including hairy roots[46], adventitious roots[47], different tissues[48], and transgenic materials[49]. However, it did not perform well in this study. Liu et al.[21] reported that the CYP and EF-1α genes were the most stable reference genes at different growth stages and in different organs of ginseng. Wang et al.[22] found that the EF1-γ/IF3G1 genes were the most stable reference genes in different tissues and that the IF3G1/ACT11 were the most stable reference genes in seedlings grown under heat stress condition. In this study, we investigated the stability of 12 candidate reference genes including 10 traditional housekeeping genes tested by Liu et al.[21] and two other housekeeping genes tested by Wang et al.[22]. Our results showed that the CYP gene was the most stable reference gene that was consistent with the result of Liu et al.[21].
To further confirm whether CYP is the most suitable reference gene in ginseng hairy roots under the MeJA treatment, four genes that encode key enzymes involved in ginsenoside biosynthesis were used as the positive controls. The expression of PgFPS was up-regulated under MeJA treatment as reported by Kim et al.[32]. PgCYP716A47 also responded to MeJA treatment as reported by Han et al.[35]. PgUGT71A27[34] and PgDDS[33] were reported to participate in ginsenoside biosynthesis. Since they have been already showed to relate with ginsenoside biosynthesis, it is expected that their expressions should be also correlated with the change of ginsenoside content. Our results showed that the relative expression levels of the genes used CYP as the reference gene were more significantly correlated with the ginsenoside content than those used 18S as the reference gene. Especially for PgCYP716A47, its expression was significantly correlated with the content of 8 kinds of ginsenosides when CYP was used as the reference gene, while its expression had no significant relevance with the content of any kinds of ginsenosides when 18S was used as the reference gene. These results indicate that our analysis of the candidate reference genes is reliable and the CYP gene can be used as the reference gene for qRT-PCR when studying the genes involved in ginsenoside biosynthesis. This study lays a foundation for further explorations of the genes involved in ginsenoside biosynthesis and the molecular mechanism underlying the process.
Supporting information
S1 Fig [a]
The induction and culture of ginseng hairy roots.
S2 Fig [a]
The chromatogram of standard products and samples.
Zdroje
1. Ahuja A, Kim JH, Kim J-H, Yi Y-S, Cho JY. Functional role of ginseng-derived compounds in cancer. Journal of ginseng research. 2018;42(3):248–254. doi: 10.1016/j.jgr.2017.04.009 29983605
2. Park SB, Park GH, Um Y, Kim HN, Song HM, Kim N, et al. Wood-cultivated ginseng exerts anti-inflammatory effect in LPS-stimulated RAW264. 7 cells. International journal of biological macromolecules. 2018;116:327–334. doi: 10.1016/j.ijbiomac.2018.05.039 29751039
3. Kim J-H. Pharmacological and medical applications of Panax ginseng and ginsenosides: A review for use in cardiovascular diseases. Journal of ginseng research. 2018;42(3):264–269. doi: 10.1016/j.jgr.2017.10.004 29983607
4. Kim KH, Lee D, Lee HL, Kim C-E, Jung K, Kang KS. Beneficial effects of Panax ginseng for the treatment and prevention of neurodegenerative diseases: past findings and future directions. Journal of ginseng research. 2018;42(3):239–247. doi: 10.1016/j.jgr.2017.03.011 29989012
5. Patel S, Rauf A. Adaptogenic herb ginseng (Panax) as medical food: Status quo and future prospects. Biomedicine & Pharmacotherapy. 2017;85:120–127.
6. Adil M, Jeong BR. In vitro cultivation of Panax ginseng C.A. Meyer. Industrial crops and products. 2018;122:239–251.
7. Zhang T, Zhong S, Hou L, Wang Y, Xing X, Guan T, et al. Computational and experimental characterization of estrogenic activities of 20 (S, R)-protopanaxadiol and 20 (S, R)-protopanaxatriol. Journal of Ginseng Research. 2018.
8. Kim Y-J, Zhang D, Yang D-C. Biosynthesis and biotechnological production of ginsenosides. Biotechnology advances. 2015;33(6):717–735.
9. Zhang Y-C, Li G, Jiang C, Yang B, Yang H-J, Xu H-Y, et al. Tissue-specific distribution of ginsenosides in different aged ginseng and antioxidant activity of ginseng leaf. Molecules. 2014;19(11):17381–17399. doi: 10.3390/molecules191117381 25353387
10. Li J, Liu S, Wang J, Li J, Liu D, Li J, et al. Fungal elicitors enhance ginsenosides biosynthesis, expression of functional genes as well as signal molecules accumulation in adventitious roots of Panax ginseng C.A. Mey. Journal of biotechnology. 2016;239:106–114. doi: 10.1016/j.jbiotec.2016.10.011 27746309
11. Huang C, Qian Z-G, Zhong J-J. Enhancement of ginsenoside biosynthesis in cell cultures of Panax ginseng by N, N′-dicyclohexylcarbodiimide elicitation. Journal of biotechnology. 2013;165(1):30–36. doi: 10.1016/j.jbiotec.2013.02.012 23467002
12. Tewari RK, Paek K-Y. Salicylic acid-induced nitric oxide and ROS generation stimulate ginsenoside accumulation in Panax ginseng roots. Journal of plant growth regulation. 2011;30(4):396–404.
13. Vasconsuelo A, Boland R. Molecular aspects of the early stages of elicitation of secondary metabolites in plants. Plant Science. 2007;172(5):861–875.
14. Yu K-W, Gao W-Y, Son S-H, Paek K-Y. Improvement of ginsenoside production by jasmonic acid and some other elicitors in hairy root culture of ginseng (Panax ginseng C.A. Meyer). In Vitro Cellular & Developmental Biology-Plant. 2000;36(5):424–428.
15. Palazón J, Cusidó RM, Bonfill M, Mallol A, Moyano E, Morales C, et al. Elicitation of different Panax ginseng transformed root phenotypes for an improved ginsenoside production. Plant Physiology and Biochemistry. 2003;41(11–12):1019–1025.
16. Ali MB, Yu K-W, Hahn E-J, Paek K-Y. Differential responses of anti-oxidants enzymes, lipoxygenase activity, ascorbate content and the production of saponins in tissue cultured root of mountain Panax ginseng CA Mayer and Panax quinquefolium L. in bioreactor subjected to methyl jasmonate stress. Plant science. 2005;169(1):83–92.
17. Cao H, Nuruzzaman M, Xiu H, Huang J, Wu K, Chen X, et al. Transcriptome analysis of methyl jasmonate-elicited Panax ginseng adventitious roots to discover putative ginsenoside biosynthesis and transport genes. International journal of molecular sciences. 2015;16(2):3035–3057. doi: 10.3390/ijms16023035 25642758
18. Oh JY, Kim Y-J, Jang M-G, Joo SC, Kwon W-S, Kim S-Y, et al. Investigation of ginsenosides in different tissues after elicitor treatment in Panax ginseng. Journal of ginseng research. 2014;38(4):270–277. doi: 10.1016/j.jgr.2014.04.004 25379007
19. Udvardi MK, Czechowski T, Scheible W-R. Eleven golden rules of quantitative RT-PCR. The Plant Cell. 2008;20(7):1736–1737. doi: 10.1105/tpc.108.061143 18664613
20. Chapman JR, Waldenström J. With reference to reference genes: a systematic review of endogenous controls in gene expression studies. PloS one. 2015;10(11):e0141853. doi: 10.1371/journal.pone.0141853 26555275
21. Liu J, Wang Q, Sun M, Zhu L, Yang M, Zhao Y. Selection of reference genes for quantitative real-time PCR normalization in Panax ginseng at different stages of growth and in different organs. PloS one. 2014;9(11):e112177. doi: 10.1371/journal.pone.0112177 25393243
22. Wang M, Lu S. Validation of suitable reference genes for quantitative gene expression analysis in Panax ginseng. Frontiers in plant science. 2016;6:1259. doi: 10.3389/fpls.2015.01259 26793228
23. Chi C, Shen Y, Yin L, Ke X, Han D, Zuo Y. Selection and validation of reference genes for gene expression analysis in Vigna angularis using quantitative real-time RT-PCR. PloS one. 2016;11(12):e0168479. doi: 10.1371/journal.pone.0168479 27992593
24. Hu Y, Deng T, Chen L, Wu H, Zhang S. Selection and Validation of Reference Genes for qRT-PCR in Cycas elongata. PloS one. 2016;11(4):e0154384. doi: 10.1371/journal.pone.0154384 27124298
25. Reddy DS, Bhatnagar-Mathur P, Reddy PS, Cindhuri KS, Ganesh AS, Sharma KK. Identification and validation of reference genes and their impact on normalized gene expression studies across cultivated and wild cicer species. PloS one. 2016;11(2):e0148451. doi: 10.1371/journal.pone.0148451 26863232
26. Chang Y-W, Chen J-Y, Lu M-X, Gao Y, Tian Z-H, Gong W-R, et al. Selection and validation of reference genes for quantitative real-time PCR analysis under different experimental conditions in the leafminer Liriomyza trifolii (Diptera: Agromyzidae). PloS one. 2017;12(7):e0181862. doi: 10.1371/journal.pone.0181862 28746411
27. Dai T-M, Lü Z-C, Liu W-X, Wan F-H. Selection and validation of reference genes for qRT-PCR analysis during biological invasions: The thermal adaptability of Bemisia tabaci MED. PloS one. 2017;12(3):e0173821. doi: 10.1371/journal.pone.0173821 28323834
28. Liu Y, Liu J, Xu L, Lai H, Chen Y, Yang Z, et al. Identification and validation of reference genes for seashore paspalum response to abiotic stresses. International journal of molecular sciences. 2017;18(6):1322.
29. Yue Y, Qiu Z-D, Qu X-Y, Deng A-P, Yuan Y, Huang L-Q. Discoursing on Soxhlet extraction of ginseng using association analysis and scanning electron microscopy. Journal of pharmaceutical analysis. 2018;8(5):312–317. doi: 10.1016/j.jpha.2018.08.003 30345145
30. Yue-Yi C JS, En Y U, et al. Study on transformation of ginsenoside Rg_3 fermented by Monascus purpureus. Chinese Traditional and Herbal Drugs. 2018. (In Chinese)
31. Yu-Nan Z, Zhong-Li W, Jian-Guo D, Lin C, HUANG Y-F. Preparation and quality assessment of high-purity ginseng total saponins by ion exchange resin combined with macroporous adsorption resin separation. Chinese journal of natural medicines. 2014;12(5):382–392. doi: 10.1016/S1875-5364(14)60048-0 24856763
32. Kim O, Bang K, Jung S, Kim Y, Hyun D, Kim S, et al. Molecular characterization of ginseng farnesyl diphosphate synthase gene and its up-regulation by methyl jasmonate. Biologia plantarum. 2010;54(1):47–53.
33. Han JY, Kwon YS, Yang DC, Jung YR, Choi YE. Expression and RNA interference-induced silencing of the dammarenediol synthase gene in Panax ginseng. Plant and cell physiology. 2006;47(12):1653–1662. doi: 10.1093/pcp/pcl032 17088293
34. Jung S-C, Kim W, Park SC, Jeong J, Park MK, Lim S, et al. Two ginseng UDP-glycosyltransferases synthesize ginsenoside Rg3 and Rd. Plant and Cell Physiology. 2014;55(12):2177–2188. doi: 10.1093/pcp/pcu147 25320211
35. Han J-Y, Kim H-J, Kwon Y-S, Choi Y-E. The Cyt P450 enzyme CYP716A47 catalyzes the formation of protopanaxadiol from dammarenediol-II during ginsenoside biosynthesis in Panax ginseng. Plant and cell physiology. 2011;52(12):2062–2073. doi: 10.1093/pcp/pcr150 22039120
36. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome biology. 2002;3(7):research0034. 1.
37. Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer research. 2004;64(15):5245–5250. doi: 10.1158/0008-5472.CAN-04-0496 15289330
38. Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnology letters. 2004;26(6):509–515. doi: 10.1023/b:bile.0000019559.84305.47 15127793
39. Kim Y-S, Hahn E-J, Murthy HN, Paek K-y. Adventitious root growth and ginsenoside accumulation in Panax ginseng cultures as affected by methyl jasmonate. Biotechnology Letters. 2004;26(21):1619–1622. doi: 10.1007/s10529-004-3183-2 15604808
40. Kang KB, Jayakodi M, Lee YS, Nguyen VB, Park H-S, Koo HJ, et al. Identification of candidate UDP-glycosyltransferases involved in protopanaxadiol-type ginsenoside biosynthesis in Panax ginseng. Scientific reports. 2018;8(1):11744. doi: 10.1038/s41598-018-30262-7 30082711
41. Verma N, Shukla S. Impact of various factors responsible for fluctuation in plant secondary metabolites. Journal of Applied Research on Medicinal and Aromatic Plants. 2015;2(4):105–113.
42. Thines B, Katsir L, Melotto M, Niu Y, Mandaokar A, Liu G, et al. JAZ repressor proteins are targets of the SCF COI1 complex during jasmonate signalling. Nature. 2007;448(7154):661. doi: 10.1038/nature05960 17637677
43. Wasternack C, Strnad M. Jasmonates are signals in the biosynthesis of secondary metabolites—Pathways, transcription factors and applied aspects—A brief review. New biotechnology. 2019;48:1–11. doi: 10.1016/j.nbt.2017.09.007 29017819
44. Per TS, Khan MIR, Anjum NA, Masood A, Hussain SJ, Khan NA. Jasmonates in plants under abiotic stresses: Crosstalk with other phytohormones matters. Environmental and experimental botany. 2018;145:104–120.
45. Shen Q, Lu X, Yan T, Fu X, Lv Z, Zhang F, et al. The jasmonate‐responsive AaMYC2 transcription factor positively regulates artemisinin biosynthesis in Artemisia annua. New Phytologist. 2016;210(4):1269–1281. doi: 10.1111/nph.13874 26864531
46. Lu C, Zhao S, Wei G, Zhao H, Qu Q. Functional regulation of ginsenoside biosynthesis by RNA interferences of a UDP-glycosyltransferase gene in Panax ginseng and Panax quinquefolius. Plant physiology and biochemistry. 2017;111:67–76. doi: 10.1016/j.plaphy.2016.11.017 27914321
47. Subramaniyam S, Mathiyalagan R, Natarajan S, Kim Y-J, Jang M-g, Park J-H, et al. Transcript expression profiling for adventitious roots of Panax ginseng Meyer. Gene. 2014;546(1):89–96. doi: 10.1016/j.gene.2014.05.024 24831831
48. Zhang R, Zhu J, Cao H-Z, An Y-R, Huang J-J, Chen X-H, et al. Molecular cloning and expression analysis of PDR1-like gene in ginseng subjected to salt and cold stresses or hormonal treatment. Plant physiology and biochemistry. 2013;71:203–211. doi: 10.1016/j.plaphy.2013.07.011 23968928
49. Park S-B, Chun J-H, Ban Y-W, Han JY, Choi YE. Alteration of Panax ginseng saponin composition by overexpression and RNA interference of the protopanaxadiol 6-hydroxylase gene (CYP716A53v2). Journal of ginseng research. 2016;40(1):47–54. doi: 10.1016/j.jgr.2015.04.010 26843821
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts