search for


Gibberellins enhance plant growth and ginsenoside content in Panax ginseng
J Plant Biotechnol 2021;48:186-192
Published online September 30, 2021
© 2021 The Korean Society for Plant Biotechnology.

Chang Pyo Hong ・Gwi Yeong Jang ・Hojin Ryu

Theragen Bio Co., Ltd, Suwon 16229, Republic of Korea
Department of Herbal Crop Research, NIHHS, RDA, Eumseong 27709, Republic of Korea
Department of Biology, Chungbuk National University, Cheongju 28644, Republic of Korea
Department of Biological Sciences and Biotechnology, Chungbuk National University, Cheongju 28644, Republic of Korea
Correspondence to: e-mail:
These authors contributed equally to this work as the first authors
Received September 4, 2021; Revised September 23, 2021; Accepted September 23, 2021.
cc This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
The roots of Korean ginseng (Panax ginseng) have a long history of usage as a medicinal drug. Ginsenosides, a group of triterpenioid saponins in ginseng, have been reported to show important pharmacological effects. Many studies have attempted to identify the ginsenoside synthesis pathways of P. ginseng and to increase crop productivity. Recent studies have shown that exogenous gibberellin (GA) treatments promote storage root secondary growth by integration of the modulating cambium stem cell homeostasis with a secondary cell wall-related gene network. However, the dynamic regulation of ginsenoside synthesis related genes and their contents by external signaling cues has been rarely evaluated. In this study, we confirmed that GA treatment not only enhanced the secondary growth of P. ginseng storage roots, but also significantly enriched the terpenoid biosynthesis process in RNA-seq analysis. Consistently, we also found that the expression of most genes involved in the ginsenoside synthesis pathways, including those encoding methylerythritol-4-phosphate (MEP) and mevalonate (MVA), and the saponin content in both leaves and roots was increased by exogenous GA application. These results can be used in future development of biotechnology for ginseng breeding and enhancement of saponin content.
Keywords : Panax ginseng, GA, Hormone, Ginsenoside, Terpenoid biosynthesis

Korean ginseng (P. ginseng) is a medicinal herb that contains a lot of functional substances such as saponins and is recognized worldwide for its medicinal efficacy (Baeg and So 2013; Hong et al. 2018; Kaneko and Nakanishi 2004). Ginsenoside, a major physiologically active substances in ginseng, is attracting attention for its excellent pharmacological action. It has been reported to have various effects such as anti-cancer, anti-fatigue, and anti-diabetes etc. (Lee and Kim 2014). The recently advanced integration of genomics and analytical chemistry research has increased our understanding of ginsenoside biosynthesis in plants (Bergman et al. 2019; Hong et al. 2018; Zhao et al. 2014). Ginsenosides are composed of several types of triterpenoid saponins. Terpenoids are synthesized from the precursor isopentenyl pyrophosphate (IPP), which is biosynthesized in both cytoplasmic MVA and plastid MEP pathways. (Dubey et al. 2003; Kim et al. 2018; Zhao et al. 2014). The 2,3-oxidosqualene produced via the MVA and MEP pathways, is biosynthesized to ginsenosides via cyclization, hydroxylation, and finally glycosylation (Kim et al. 2018; Zhao et al. 2014). Similar to most plants, the 2,3-oxidosqualene is cyclized into cycloartenol, α-, β-amyrin or lupeol, to be further converted to phytosterols and pentacyclic triterpenoids in P. ginseng (Benveniste 2004; Kim et al. 2018). Then, the major dammarane-type saponins [protopanaxadiol (PPD)/pro- topanaxatriol (PPT)] are synthesized by a specific cyclization reaction activated by a set of cytochrome P450 enzymes. It is further glycosylated via several types of UDP-glycosyltransferases (UGTs) to synthesize ginsenosides. Most of genes involved in the ginsenoside synthesis pathway have been identified through recent studies (Benveniste 2004; Kim et al. 2018).

The secondary growth driving from cambium stem cells is a major factor for the yield of P. ginseng. It has been recently reported that secondary growth of storage root in P. ginseng is regulated by brassinosteroid (BR) and gibberellin (GA) (Hwang et al. 2020; Lee et al. 2018). BR signaling activated by prolonged high ambient temperature inhibited the cell division of cambium stem cells in the tap roots of P. ginseng by yet unknown mechanisms (Hwang et al. 2020). Our recent study shows that exogenous GA treatments promoted the secondary growth of ginseng storage roots, and the gene network related to secondary cell wall synthesis is critically involved in this process (Hong et al. 2021). The RNA-seq analysis of GA-treated ginseng roots showed that the terms related to secondary cell wall synthesis formed a close network with genes related to cell division, cell elongation, jasmonic acid (JA) and nitrogen assimilation for promoting secondary growth of the ginseng roots (Hong et al. 2021). In addition, it was recently reported that genes related to canonical GA signaling pathway are evolutionary conserved in the genome of P. ginseng (Hong et al. 2021; Jo et al. 2017; Kim et al. 2021). These results indicate that hormones are essentially involved in the secondary growth of the storage root of P. ginseng.

Recently, as the concept of well-being and consumer demand for eco-friendly agricultural products increase, research on the development of eco-friendly ginseng cultivation technology is being widely conducted. The well-being trend is expected to emerge as one of the big issues in the future ginseng market. With these demands, recent studies of P. ginseng have been focused on the mechanism of pharmacological action in animal model systems and manual development for eco-friendly cultivation (Bergman et al. 2019; Lee et al. 2020; Ratan et al. 2021; Wang et al. 2020). In addition, research is underway to measure the change in ginsenoside content of ginseng sprouts caused by various external environmental factors in the greenhouse system (Li et al. 2016; Huang et al. 2013; Rahimi et al. 2016; Zhang et al. 2021). However, genetic and physiological studies related to the increase in the yield of P. ginseng and the synthesis of saponin content have been rarely conducted. In this study, we found that the secondary growth and ginsenoside content was increased in GA-treated ginseng roots. It was also consistently confirmed that the expression of genes involved in ginsenoside synthesis through RNA-seq analysis was significantly increased by GA treatment in the ginseng roots. Finally, our results also showed that exogenous GA treatment in roots increased plant growth and saponin content in both roots and leaves.

Materials and Methods

Plant growth and RNA-seq analysis

One-year-old P. ginseng seedlings (Yunpoong, kindly supplied by National Institute of Horticultural and Herbal Science) were grown on the ginseng cultivation soil medium. After two weeks, the developed ginseng seedlings were treated 10 μM GA3 and DMSO (mock treatment) once a week for eight weeks with a soaking method. The GA-treated ginseng plants were sampled, and the patterns of shoot primary growth and main root secondary growth were further analyzed.

RNA-Seq reads of DMSO- and GA-treated ginseng roots, of which three replicates have deposited with the accession number SAMN12240097 and SAMN122731128in NCBI Short Read Archive database (Hong et al. 2021), were reanalyzed for this study. The RNA-seq analysis was performed as described in our previous work (Lee et al., 2021). Genes differentially expressed between each replicate of DMSO- and GA-treated root samples were identified using DESeq2 (Love et al. 2014), based on cutoff values of q < 0.05 and FC ≥ 1.5. GO enrichment analysis of DEGs homologous to Arabidopsis genes (TAIR 10 release) was performed using DAVID, with an EASE score cutoff of < 0.01. Genes enriched under specific GO terms were selected, and their expression patterns were visualized as a heatmap using MeV ( Expression levels of genes were shown as Z-score of TPM values. To test if a particular gene set (i.e., Terpenoid biosynthesis proecess) was enriched, GSEA (Subramanian et al. 2005; was applied to a background dataset comprising transcripts with TPM> 0.3, which balances the numbers of false positives and false negatives, in either DMSO- or GA-treated root samples.

Measurement of ginsenoside composition

The dried ginseng roots and leaves were ground using a hammer mill (Micro hammer cutter mill type-3, Culatti AG, Zurich, Switzerland), and these powders were used for analysis. A 50 mg of sample (ginseng roots and leaves powder) was extracted with 2 mL of 70% ethanol solution using a sonicator at room temperature for 1 hr, and the extraction process was repeated three times. The extract was diluted to 10-fold dilution with distilled water, loaded on to Diaion HP-20 open column (20 mm × 10 mm), sequentially eluted with distilled water and methanol, and the methanol fraction was used to HPLC analysis of ginsenoside composition.

Ginsenoside standards (Rg1, Re, Rb1, Rg2(S), Rc, Rg2(R), Rb2, Rb3, F1, Rg6, F4, Rk3, Rh4, Rg3(S), Rg3(R), protopanaxatrio, Rk1, Compound K, Rg5, Rh2, and protopanaxadiol) were purchased from Chengdu BiopurifyPhytochemicals Ltd. (Chengdu, Sichuan, China). HPLC-grade water and acetonitrile (ACN) were purchased from J.T. Baker (Phillipsburg, NJ, USA). Diaion HP20 resin was purchased from Sigma Aldrich (USA). All other chemicals used were of reagent grade.

The ginsenoside composition was analyzed using a high- performance liquid chromatography (HPLC) with a UV-visible detector (HPLC system: 1200 series, Agilent technologies, Santa Clara, CA, USA; column: Kinetex XB-C18, 100 × 4.6 mm, 2.6 µm, Phenomenex, Torrance, CA, USA). The mobile phase of the analytical system consisted of ACN (A) and water (B), using the following gradient: 0 min (18% A), 0–5 min (18~22% A), 5–15 min (22~22% A), 15–19 min (22~35% A), 19–25 min (35~35% A), 25-31 min (35-50% A), 31-51 min (50~65% A), 51-55 min (65~100% A), and 55-63 min (100~100% A). The flow rate, detection wavelength, and injection volume were set at 1.0 mL/min, 203 nm, and 10 µL, respectively.

Results and Discussion

Our previous study showed that exogenous GA treatments greatly enhanced the secondary growth of the tap roots of P. ginseng (Hong et al. 2021). To further investigate the physiological effects of GA in plant growth and secondary metabolites in P. ginseng, we reanalyzed the phenotypes and RNA-seq data of GA3-treated roots of 1-year-old P. ginseng. As previously reported, exogenous application of GA significantly promoted both stem growth and tap root grith of ginseng compared with DMSO treatment (Figure 1A, B and C). GA treatment increased the length of the ginseng stem and the area of the root by about 75% and 50%, respectively (Fig. 1B and C). To investigate the relationship between plant growth promotion and the accumulation of secondary metabolites, we next conducted that the GO term related to the metabolic process was identified in the previously reported RNA-seq data of GA3-treated ginseng roots (Hong et al. 2021). As present in Figure. 1D, the gene set enrichment analysis (GSEA) of differentially expressed genes (DEGs) revealed that the upregulated genes in GA treatment were significantly enriched under the functional categories of terpenoid biosynthesis process (FDR : 0.0; Fig. 1D and Table 1). Terpenoid is a representative plant secondary metabolite involved in the synthesis of critical for substances in medicine uses (Bergman et al. 2019; Lee and Kim 2014; Benveniste 2004). These results indicate that GA treatment has a great influence on the promotion of plant growth and the synthesis pathway of secondary metabolites of P. ginseng.

Expression of genes involved in triterpene saponin biosynthetic pathway in Panax ginseng

Fig. 1. Exogenous GA treatment affects plant growth and the expression pattern of terpenoid biosynthesis-related genes. (A) GA treatment enhances primary and secondary growth of shoot and storage roots, respectively. One-year-old P. ginseng plants were treated with DMSO as a mock control or 10 μM GA3 once a week for 8 weeks. Scale bar = 2 cm. (B, C) Measurements of shoot length (B) and root grith (C). Error bars represent the standard error (n > 18). Different lowercase letters indicate statistically significant differences (P < 0.05; one-way ANOVA, followed by Tukey’s multiple range test). (D) Enrichment plot for the terpenoid biosynthesis process and an expression heatmap of 19 genes related to this pathway (FDR = 0.0) were obtained from the RNA-seq data of the roots of GA3-treated P. ginseng (Hong et al. 2021)

Since the synthesis of terpenoids is used as a precursor for the synthesis of ginsenoside in ginseng, we next investigated the expression patterns of genes involved in saponin biosynthetic pathways (Bergman et al. 2019; Kim et al. 2018). It was confirmed that most of the expression of AACT (Acetyl-CoA C-acetyltransferase), HMGS (Hydroxymethylglutaryl-CoA synthase), HMGR (Hydroxymethylglutaryl-CoA reductase), MVK (Mevalonate kinase) and PMK (Phosphomevalonate kinase) genes involved in the MVA pathway occurring in the cytoplasm was increased by GA treatment (Fig. 2A). In addition, we found that the expression of most genes involved in another MEP (occurred in plastid) pathway for IPP synthesis was highly maintained by GA treatment (Fig. 2A). IPP synthesized by either the MEP or MVA pathways is converted to triterpene and then synthesized to 2,3-oxidosqualene by several SQEs (Squalene epoxidase). Twelve SQE genes are known in the genome of ginseng, and the expression of 8 SQE genes was increased by GA treatment (Fig. 2B). The 2,3-oxidosqualene was then cyclized to β-amyrin by the b-AS (Squalene cyclase) family, but only 3 genes out of a total of 8 gene families were upregulated in GA-treated ginseng roots. However, the expression of genes involved in the synthesis of dammaranediol, a cyclic compound additionally synthesized by 2,3-oxidosqualene, and the synthesis of major dammarane-type sapogenin [protopanaxadiol (PPD)/protopanaxatriol (PPT)] were mostly upregulated by exogenous GA3 (Fig. 2B). Finally, we confirmed that most UDP-glycosyltransferases that finally synthesize ginsenoside by further glycosylation of these precursors was also highly expressed in the GA-treated ginseng (Fig. 2B). Taken together, these results suggest that the terpenoid synthesis pathway increased by GA treatment would positively control the synthesis of various ginsenosides.

Fig. 2. Transcriptome analysis of the saponin biosynthetic pathways of GA-treated ginseng roots. The ginsenoside biosynthetic pathway and expression levels of its related genes in P. ginseng are presented. The expression heatmaps of genes involved in ginsenoside biosynthetic pathways (A; MEP and MVA pathways, B; downstream pathway) were obtained from the RNA-seq data of the roots of GA3-treated P. ginseng (Hong et al. 2021)

We next measure the content of various types of ginsenosides in GA-treated ginseng plants. Typical chromatograms of ginsenosides for P. ginseng roots and leaves were initially measured (Fig. 3). Thereafter, the contents of each ginsenoside were comparatively analyzed in the leaves and roots. As shown in Figure. 4, the ginsenoside contents were generally measured to be high in the GA3-treated ginseng samples. In particular, the ginsenoside contents were higher in leave samples than in roots. Re and Rb2 types were measured to be high in leaves, and most ginsenosides except for F1 and Rg3(R) were increased by about 10% in GA3-treated ginseng root samples (Fig. 4). However, the content of Rb3 did not show a significant difference in both leaves and roots (Fig. 4). Taken together, these results indicate that GA treatment can increase the expression and content of genes involved in the synthesis of ginsenoside in ginseng.

Fig. 3. Typical chromatograms of P. ginseng roots and leaves treated with or without GA3. Samples: A, ginsenoside standards; B, roots (GA3); C, roots (Con.); D, leaves (GA3); E, leaves (Con.). Ginsenoside standards: 1, Rg1; 2, Re; 3, Rb1; 4, Rg2(S); 5, Rc; 6, Rg2(R); 7, Rb2; 8, Rb3; 9, F1; 10, Rg6; 11, F4; 12, Rk3; 13, Rh4; 14, Rg3(S); 15, Rg3(R); 16, protopanaxatriol; 17, Rk1; 18, Compound K; 19, Rg5; 20, Rh2; 21, protopanaxadiol

Fig. 4. Ginsenoside content of P. ginseng roots and leaves after GA treatment. The ginsenoside content was measured. Error bars indicate SEM (n = 3; *p < 0.05; **p < 0.01; Student’s t-test). ns: not significant, nd: not detected

In this study, we confirmed once again that the primary and secondary growth of P. ginseng was promoted by exogenous GA treatments. Our recent study demonstrates that genes involved in GA signaling in the ginseng genome are evolutionarily well conserved and functional genomic analysis (Kim et al. 2021; Hong et al. 2021). In addition, terms related to the terpenoid synthesis pathway among DEG whose expression is increased by GA treatment are very significantly related, showing that the expression of genes related to saponin synthesis in ginseng is specifically increased. Several studies have demonstrated that transcription factors activated by JA signaling directly reduced GA synthesis and signaling (Heinrich et al. 2013; Cheng et al. 2009). JA is well known as a representative substance that enhances saponin synthesis in the invitro culture of ginseng. This study describes the effect of short-term GA treatment using one-year-old ginseng seedlings. However, several experimental results suggest that GA would enhance saponin synthesis through direct interaction with JA as well as promoting plant growth. These results can be expected to promote the growth of ginseng seedlings and increase saponin in a short period of time, and furthermore, they can be used as important clues in the breeding of P. ginseng.


This research was supported by Chungbuk National University Korea National University Development Project (2021).

Conflict of Interest Disclosures

All Authors have read the manuscript and declared that they have no conflict of interest.

  1. Baeg IH, So SH (2013) The world ginseng market and the ginseng (Korea). J Ginseng Res 37(1): 1-7
    Pubmed KoreaMed CrossRef
  2. Benveniste P (2004) Biosynthesis and accumulation of sterols. Annu Rev Plant Biol 55: 429-457
    Pubmed CrossRef
  3. Bergman ME, Davis B, Phillips MA (2019) Medically Useful Plant Terpenoids: Biosynthesis, Occurrence, and Mechanism of Action. Molecules 24(21)
    Pubmed KoreaMed CrossRef
  4. Cheng H, Song S, Xiao L, Soo HM, Cheng Z, Xie D, Peng J (2009) Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLoS Genet 5(3)
    Pubmed KoreaMed CrossRef
  5. Heinrich M, Hettenhausen C, Lange T, Wunsche H, Fang J, Baldwin IT, Wu J (2013) High levels of jasmonic acid antagonize the biosynthesis of gibberellins and inhibit the growth of Nicotiana attenuata stems. Plant J 73(4): 591-606
    Pubmed CrossRef
  6. Hong CP, Kim J, Lee J, Yoo S-i, Bae W, Geem KR, Yu J, Jang I-b, Jo IH, Cho H, Shim D, Ryu H (2021) Gibberellin signaling promotes the secondary growth of storaage roots in Panax ginseng. Int J Mol Sci 22: 8694
    Pubmed KoreaMed CrossRef
  7. Hong J, Kim H, Ryu H (2018) Identification of ABSCISIC ACID (ABA) signaling related genes in Panax ginseng. J Plant Biotechnol (45): 306-314
  8. Huang C, Qian ZG, Zhong JJ (2013) Enhancement of ginsenoside biosynthesis in cell cultures of Panax ginseng by N,N'- dicyclohexylcarbodiimide elicitation. J Biotechnol 165(1): 30-36
    Pubmed CrossRef
  9. Hwang H, Lee HY, Ryu H, Cho H (2020) Functional Characterization of BRASSINAZOLE-RESISTANT 1 in Panax Ginseng (PgBZR1) and Brassinosteroid Response during Storage Root Formation. Int J Mol Sci 21(24)
    Pubmed KoreaMed CrossRef
  10. Jo IH, Lee J, Hong CE, Lee DJ, Bae W, Park SG, Ahn YJ, Kim YC, Kim JU, Lee JW, Hyun DY, Rhee SK, Hong CP, Bang KH, Ryu H (2017) Isoform Sequencing Provides a More Comprehensive View of the Panax ginseng Transcriptome. Genes (Basel) 8(9)
    Pubmed KoreaMed CrossRef
  11. Kaneko H, Nakanishi K (2004) Proof of the mysterious efficacy of ginseng: basic and clinical trials: clinical effects of medical ginseng, korean red ginseng: specifically, its anti-stress action for prevention of disease. J Pharmacol Sci 95(2): 158-162
    Pubmed CrossRef
  12. Kim NH, Jayakodi M, Lee SC, Choi BS, Jang W, Lee J, Kim HH, Waminal NE, Lakshmanan M, van Nguyen B, Lee YS, Park HS, Koo HJ, Park JY, Perumal S, Joh HJ, Lee H, Kim J, Kim IS, Kim K, Koduru L, Kang KB, Sung SH, Yu Y, Park DS, Choi D, Seo E, Kim S, Kim YC, Hyun DY, Park YI, Kim C, Lee TH, Kim HU, Soh MS, Lee Y, In JG, Kim HS, Kim YM, Yang DC, Wing RA, Lee DY, Paterson AH, Yang TJ (2018) Genome and evolution of the shade-requiring medicinal herb Panax ginseng. Plant Biotechnol J 16(11): 1904-1917
    Pubmed KoreaMed CrossRef
  13. Kim J, Shin WR, Kim YH, Shim D, Ryu H (2021) Functional characterization of gibberellin signaling related genes in Panax ginseng. . J Plant Biotech 48(3): 148-155
  14. Lee CH, Kim JH (2014) A review on the medicinal potentials of ginseng and ginsenosides on cardiovascular diseases. J Ginseng Res 38(3): 161-166
    Pubmed KoreaMed CrossRef
  15. Lee DK, Park S, Long NP, Min JE, Kim HM, Yang E, Lee SJ, Lim J, Kwon SW (2020) Research Quality-Based Multivariate Modeling for Comparison of the Pharmacological Effects of Black and Red Ginseng. Nutrients 12(9)
    Pubmed KoreaMed CrossRef
  16. Lee J, Shim D, Moon S, Kim H, Bae W, Kim K, Kim YH, Rhee SK, Hong CP, Hong SY, Lee YJ, Sung J, Ryu H (2018) Genome- wide transcriptomic analysis of BR-deficient Micro-Tom reveals correlations between drought stress tolerance and brassinosteroid signaling in tomato. Plant Physiol Biochem 127: 553-560
    Pubmed CrossRef
  17. Lee J, Kim H, Park SG, Hwang H, Yoo Si, Bae W, Kim E, Kim J, Lee HY, Heo TY (2021) Brassinosteroid‐BZR1/2‐WAT1 module determines the high level of auxin signalling in vascular cambium during wood formation. New Phytologist 230(4): 1503-1516
    Pubmed CrossRef
  18. Li J, Liu S, Wang J, Li J, Liu D, Li J, Gao W (2016) Fungal elicitors enhance ginsenosides biosynthesis, expression of functional genes as well as signal molecules accumulation in adventitious roots of Panax ginseng C. A. Mey. J Biotechnol 239: 106-114
    Pubmed CrossRef
  19. Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome biology 15(12): 1-21
    Pubmed KoreaMed CrossRef
  20. Rahimi S, Kim YJ, Sukweenadhi J, Zhang D, Yang DC (2016) PgLOX6 encoding a lipoxygenase contributes to jasmonic acid biosynthesis and ginsenoside production in Panax ginseng. J Exp Bot 67(21): 6007-6019
    Pubmed KoreaMed CrossRef
  21. Ratan ZA, Haidere MF, Hong YH, Park SH, Lee JO, Lee J, Cho JY (2021) Pharmacological potential of ginseng and its major component ginsenosides. J Ginseng Res 45(2): 199-210
    Pubmed KoreaMed CrossRef
  22. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences 102(43): 15545-15550
    Pubmed KoreaMed CrossRef
  23. Wang C, Liu J, Deng J, Wang J, Weng W, Chu H, Meng Q (2020) Advances in the chemistry, pharmacological diversity, and metabolism of 20(R)-ginseng saponins. J Ginseng Res 44(1): 14-23
    Pubmed KoreaMed CrossRef
  24. Zhang Y, Jin T, Ryu G, Gao Y (2021) Effects of screw configuration on chemical properties and ginsenosides content of extruded ginseng. Food Sci Nutr 9(1): 251-260
    Pubmed KoreaMed CrossRef
  25. Zhao S, Wang L, Liu L, Liang Y, Sun Y, Wu J (2014) Both the mevalonate and the non-mevalonate pathways are involved in ginsenoside biosynthesis. Plant Cell Rep 33(3): 393-400
    Pubmed CrossRef

September 2022, 49 (3)
Full Text(PDF) Free

Social Network Service

Cited By Articles

Funding Information
  • CrossMark
  • Crossref TDM