Selection of Reference Genes for Quantitative Real-Time PCR Normalization in Panax ginseng at Different Stages of Growth and in Different Organs Jing Liu 1 , Qun Wang 1 , Minying Sun 1 , Linlin Zhu 1 , Michael Yang 2 , Yu Zhao 1 * 1 Traditional Chinese Medicine and Biotechnology Research and Development Center, Changchun University of Traditional Chinese Medicine, Changchun, Jilin, People’s Republic of China, 2 College of Dental Medicine, Columbia University, New York, United States of America Abstract Quantitative real-time reverse transcription PCR (qRT-PCR) has become a widely used method for gene expression analysis; however, its data interpretation largely depends on the stability of reference genes. The transcriptomics of Panax ginseng, one of the most popular and traditional ingredients used in Chinese medicines, is increasingly being studied. Furthermore, it is vital to establish a series of reliable reference genes when qRT-PCR is used to assess the gene expression profile of ginseng. In this study, we screened out candidate reference genes for ginseng using gene expression data generated by a high-throughput sequencing platform. Based on the statistical tests, 20 reference genes (10 traditional housekeeping genes and 10 novel genes) were selected. These genes were tested for the normalization of expression levels in five growth stages and three distinct plant organs of ginseng by qPCR. These genes were subsequently ranked and compared according to the stability of their expressions using geNorm, NormFinder, and BestKeeper computational programs. Although the best reference genes were found to vary across different samples, CYP and EF-1a were the most stable genes amongst all samples. GAPDH/30S RPS20, CYP/60S RPL13 and CYP/QCR were the optimum pair of reference genes in the roots, stems, and leaves. CYP/60S RPL13, CYP/eIF-5A, aTUB/V-ATP, eIF-5A/SAR1, and aTUB/pol IIa were the most stably expressed combinations in each of the five developmental stages. Our study serves as a foundation for developing an accurate method of qRT-PCR and will benefit future studies on gene expression profiles of Panax Ginseng. Citation: Liu J, Wang Q, Sun M, Zhu L, Yang M, et al. (2014) Selection of Reference Genes for Quantitative Real-Time PCR Normalization in Panax ginseng at Different Stages of Growth and in Different Organs. PLoS ONE 9(11): e112177. doi:10.1371/journal.pone.0112177 Editor: Robert W. Dettman, Northwestern University, United States of America Received November 29, 2013; Accepted October 13, 2014; Published November 13, 2014 Copyright: ß 2014 Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This work was supported by the National Natural Science Foundation of China (no. 81373937, URL: http://www.nsfc.gov.cn/Portal0/default152.htm) and the National Key Technology R&D Program of China (no. 2011BAI03B01-04, URL: http://www.863.gov.cn/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * Email: [email protected]Introduction Ginseng (Panax ginseng C.A. Meyer) is a perennial herb and is well-known for its adaptogenic and restorative properties. It has been widely used in traditional Chinese medicine and Western herbal medicine [1,2]. Ginseng root, the most commonly used part of the plant, contains ginsenosides that are major bioactive constituents with complex and multiple pharmacological effects [3,4]. Ginseng leaf-stem extract also contains numerous important bioactive components [5,6]. A recent report demonstrated that American ginseng leaf contains similar pharmacologically active ingredients in higher quantity than found in ginseng root [7]. Research has shown that ginseng leaf-stem may as well be a valuable source of ginsenosides as ginseng root [8]. From germination to withering, the stages of growth of ginseng can be generally classified into the leaf-expansion period (LP), the flowering stage (FS), the green fruit stage (GFS), the red fruit stage (RFS), the root growing after fruit stage (RGS), and the withering stage [9]. In recent years, the research focus has expanded considerably towards elucidating the gene expression of ginseng at different developmental stages. Various researchers have high- lighted the genetic aspects of ginseng, including the marker gene identification or authentication, genes that confer resistance to environmental and biological stresses, the regulatory factors of its growth and development, and key enzymes involved in the ginsenoside biosynthetic pathway [7,8,10–15]. qRT-PCR has been widely used as a powerful technique to quantify the expression levels of transcripts. The accuracy of qRT- PCR largely depends on the stability of the reference gene(s) applied to data normalization [16]. A series of presumably stable expressed genes have been used as internal references. Some of the best known and most frequently used reference transcripts, often referred to as housekeeping genes [17], include actin (ACT), tubulin (TUB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), polyubiquitin (UBQ), and translational initiation factor (eIF). These have been extensively used as reference genes in different organisms because of their stable and uniform expression patterns [18]. However, these references have shown a significant variance when tested across species and under a broad range of experimental tests [19]. Failure to use suitable reference genes may deflect gene expression profiles and lead to misguiding results [16]. So far, there have been no reports on the use of such genes in ginseng. Therefore, it is essential to determine appropriate reference genes in order to undertake genetic engineering studies in ginseng. PLOS ONE | www.plosone.org 1 November 2014 | Volume 9 | Issue 11 | e112177
14
Embed
Selection of Reference Genes for Quantitative Real-Time ...
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Selection of Reference Genes for Quantitative Real-TimePCR Normalization in Panax ginseng at Different Stagesof Growth and in Different OrgansJing Liu1, Qun Wang1, Minying Sun1, Linlin Zhu1, Michael Yang2, Yu Zhao1*
1 Traditional Chinese Medicine and Biotechnology Research and Development Center, Changchun University of Traditional Chinese Medicine, Changchun, Jilin, People’s
Republic of China, 2 College of Dental Medicine, Columbia University, New York, United States of America
Abstract
Quantitative real-time reverse transcription PCR (qRT-PCR) has become a widely used method for gene expression analysis;however, its data interpretation largely depends on the stability of reference genes. The transcriptomics of Panax ginseng,one of the most popular and traditional ingredients used in Chinese medicines, is increasingly being studied. Furthermore, itis vital to establish a series of reliable reference genes when qRT-PCR is used to assess the gene expression profile ofginseng. In this study, we screened out candidate reference genes for ginseng using gene expression data generated by ahigh-throughput sequencing platform. Based on the statistical tests, 20 reference genes (10 traditional housekeeping genesand 10 novel genes) were selected. These genes were tested for the normalization of expression levels in five growth stagesand three distinct plant organs of ginseng by qPCR. These genes were subsequently ranked and compared according to thestability of their expressions using geNorm, NormFinder, and BestKeeper computational programs. Although the bestreference genes were found to vary across different samples, CYP and EF-1a were the most stable genes amongst allsamples. GAPDH/30S RPS20, CYP/60S RPL13 and CYP/QCR were the optimum pair of reference genes in the roots, stems,and leaves. CYP/60S RPL13, CYP/eIF-5A, aTUB/V-ATP, eIF-5A/SAR1, and aTUB/pol IIa were the most stably expressedcombinations in each of the five developmental stages. Our study serves as a foundation for developing an accuratemethod of qRT-PCR and will benefit future studies on gene expression profiles of Panax Ginseng.
Citation: Liu J, Wang Q, Sun M, Zhu L, Yang M, et al. (2014) Selection of Reference Genes for Quantitative Real-Time PCR Normalization in Panax ginseng atDifferent Stages of Growth and in Different Organs. PLoS ONE 9(11): e112177. doi:10.1371/journal.pone.0112177
Editor: Robert W. Dettman, Northwestern University, United States of America
Received November 29, 2013; Accepted October 13, 2014; Published November 13, 2014
Copyright: � 2014 Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricteduse, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the National Natural Science Foundation of China (no. 81373937, URL: http://www.nsfc.gov.cn/Portal0/default152.htm)and the National Key Technology R&D Program of China (no. 2011BAI03B01-04, URL: http://www.863.gov.cn/). The funders had no role in study design, datacollection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
The samples, processed according to the Illumina kit instruc-
tions, were prepared for the transcriptome analysis. Protocols for
the cDNA library construction, sequencing, assembly, and gene
expression level analysis have been previously described by Baojin
Yao [30]. Based on the RPKM values, the estimated gene
expression was used directly for comparing the differences in gene
expressions between samples. Distinct sequences were used for the
BLAST search and annotated against the NCBI nr database using
an E-value cut-off of 1025 [31].
Using the Illumina sequencing platform, we generated more
than 39 million high-quality sequencing reads for each sample.
After clustering via the TGICL software, more than 80,000
unigenes were produced in every database. Unigene sequences
were aligned by BLASTX to four common protein databases (Nr,
Swiss-Prot, KEGG, and COG; e-value ,0.00001). Simultaneous-
ly, we obtained the highest sequence similarities Unigenes along
with their protein functional annotations.
Selection of candidate reference genes for normalizationOn analyzing the existing databases, ten commonly used
housekeeping genes were selected as endogenous control genes.
Based on the calculated statistical values of the coefficient of
variation (CV = SD/Mean) and the maximum fold change
(MFC = Max(RPKM)/Min(RPKM)) [32], we obtained ten novel
reference genes from the 15 databases.
Table 1. Cont.
Root Stem Leaf LP FS GFS RFS RGS
CV(%) 12.65 8.03 19.65 18.90 19.53
MFC 1.39 1.21 1.49 1.46 1.50
Gene 6-PG 30SRPS20 SAR1
Mean 105.83 88.88 82.35
CV(%) 14.48 10.38 19.23
MFC 1.49 1.34 1.50
Gene ARF QCR
Mean 119.60 166.70
CV(%) 17.00 10.92
MFC 1.41 1.24
Gene V-ATP pol IIa
Mean 177.53 128.55
CV(%) 19.70 15.33
MFC 1.47 1.39
Notes: Descriptive statistics of the candidate genes based on the coefficient of variance (CV) and the maximum fold change (MFC). In total, 10 untraditional referencegenes were screened, which had the CV less than 20% and MFC less than 1.5. LP, leaf-expansion period; FS, the flower stage; GFS, the green fruit stage; RFS, the red fruitstages; RGS, the root growing after fruit stage.doi:10.1371/journal.pone.0112177.t001
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 3 November 2014 | Volume 9 | Issue 11 | e112177
In total, 20 candidate reference genes were selected, including
10 housekeeper reference genes (ACT1, GAPDH, UBQ,
18SrRNA, eIF-5A, aTUB, bTUB, CYP, F-box, and EF-1a) and
F-box, and EF-1a). Besides 18S rRNA, the RPKM value
distribution of the remaining nine housekeeping genes was in the
Figure 1. RPKM value distribution of 20 candidate reference genes. LP, leaf-expansion period; FS, the flower stage; GFS, the green fruit stage;RFS, the red fruit stages; RGS, the root growing after fruit stage.doi:10.1371/journal.pone.0112177.g001
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 4 November 2014 | Volume 9 | Issue 11 | e112177
Figure 2. Specificity of primer pairs for RT-qPCR amplification. Agarose gel (2%) electrophoresis showing amplification of a specific PCRproduct of the expected size for each gene (M:DL1000 DNA Marker).doi:10.1371/journal.pone.0112177.g002
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 5 November 2014 | Volume 9 | Issue 11 | e112177
three vegetative organs at one growth stage (Table 1-a). The CV
values of the housekeeping genes were found to vary from 3.06%
to 88.21%, while the MFC values ranged from 1.06 to 8.76. In
order to screen more stable reference genes, we set the threshold
values for CV to ,20% and MFC to ,1.5. Additionally, ten novel
genes (CDP, 6-PG, 30S RPS20, 60S RPL13, V-ATP, pol IIa,
ARF, QCR, SAR1, and TCTP) were screened as candidate
reference genes (Table 1-b). Screening of the potential reference
genes was based on the statistical tests (CV and MFC), which
reflected the RPKM values of stably and moderately or highly
expressed genes among all the databases. RPKM value expression
abundance ratios are presented in Figure 1. To determine the
distribution of transcript populations of 20 candidate reference
genes in three vegetative organs of ginseng during the five stages of
growth, the quantity of transcript for each gene was estimated as a
ratio relative to the sum of the 20 transcript populations. The
results clearly revealed a fluctuation in the relative magnitude of
RPKM values and the ratios, thus indicating that all of the 20
genes did not exhibit stable expression patterns. A summary of the
sequence information for the 20 ginseng candidate reference genes
is presented in Table 2.
Validating the expression levels of candidate referencegenes by qRT-PCR
By reverse transcription PCR, the specificity of the primers used
for candidate reference genes was verified. A single band for each
gene was revealed through electrophoresis, without primer-dimers
or non-specific amplification (Figure 2).
Based on SYBR Green detection, qRT-PCR analysis was
employed to evaluate the stability of the expressions of the 20
candidate reference genes in different organs and different
developmental stages of P. ginseng. The samples were divided
into fifteen groups comprising of three organs (roots, stems, and
leaves) and five developmental stages. The Ct values of the
reference genes of each group were then used to compare the
various degrees of expression.
Statistical data analysisThe gene expression data were analyzed by Ct value, geNorm,
NormFinder, and BestKeeper applets to obtain the expression
stability of 20 candidate reference genes.
With a higher gene expression, a smaller Ct value was obtained,
and vice versa. Figure 3 shows a relatively broad range of Ct
values for all the 20 putative reference genes. The highest Ct value
was 26.40 (bTUB), while the lowest was 15.06 (18S rRNA). Ct
values of the remaining genes were distributed between 19 and 24.
On comparing the Ct values of the 20 candidate reference genes,
the expression level of each reference gene was found to differ,
with respect to the developmental stage or the organ under study.
The expression patterns of the 20 reference genes displayed
irregular variation; this may be attributed to change in the level of
reference gene expression abundance with the cell type and the
analysis under different experimental conditions in ginseng
requires careful selection of reliable reference genes.
Based on the expression stability of the genes and the
assumption that two ideal reference genes should not vary with
each other under different test conditions [34], geNorm ranked the
best out of the three data analysis applications used. geNorm
computes the average pair-wise variation of a given candidate
reference gene with all the other genes and assigns a score of its
expression stability (M) to each gene. Stepwise exclusion of genes
with the highest M values (indicating the least stable expressions)
before recalculation finally reveals the two most stable candidate
genes [35]. After calculating the pair-wise variation Vn/n+1,
geNorm selects the optimal number of control genes. The cut-off
value is usually set to a default value of 0.15 [34]. Gene expression
stability and ranking of 20 candidate reference genes, as calculated
by geNorm using nine sets of samples, are presented in Figure 4.
Analyses of all fifteen samples revealed that the CYP and EF-1acombination showed the lowest M value (0.31), while 30S RPS20
showed the highest M value (0.89). Among the different organs,
GAPDH/30S RPS20, CYP/60S RPL13, and CYP/QCR were
the most stably expressed gene combinations in roots, stems, and
Figure 3. RT-qPCR CT values for the candidate reference genes (n = 3). Expression date displayed as CT values for each reference gene in allginseng samples. A line across the box is depicted as the median. The box indicates the 25th and 75th percentiles. Whiskers represent the maximumand minimum values.doi:10.1371/journal.pone.0112177.g003
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 6 November 2014 | Volume 9 | Issue 11 | e112177
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 7 November 2014 | Volume 9 | Issue 11 | e112177
leaves, respectively; while 18S rRNA, UBQ, and TCTP were the
least stably expressed. Among the five developmental stages under
SAR1, and aTUB/pol IIa were the most stably expressed
combination, respectively, and 30S RPS20 was the least stably
expressed gene in all the five stages. Based on these observations,
CYP was evidently the most stably expressed gene and may be
considered as the most suitable reference gene for the analyses of
gene expressions in P. ginseng. Furthermore, the addition of a
third reference gene would not have significantly increased the
statistical reliability of this calculation, as V2/3 = 0.033 or V3/
4 = 0.041 (in roots) was significantly below the default cut-off value
of 0.15 (Figure 5). Although the pair-wise variation for all the
samples (V2/3) was estimated as 0.145, it was still less than the
limiting value. Hence, our study showed that two reference genes
were sufficient to normalize gene expression for all the samples of
P. ginseng.
The NormFinder algorithm uses a model-based approach to
evaluate modifications amongst the reference gene expression
levels [36]. Similar to the geNorm method, NormFinder imparts a
score of expression stability (M) to each gene, which is negatively
correlated with the stability of gene expression [37]. In addition,
NormFinder can determine the estimated inter- and intra-group
variances [38]. The calculated values generated by NormFinder
are shown in Table 3. In the final outcome, CYP, QCR, and
aTUB show the most stable expression levels for the total samples,
stems, leaves, and the five developmental stages, while 30SRPS20
and TCTP were observed to be less stable. In the roots, GAPDH
and V-ATP were the most stably expressed genes with values of
0.013 and 0.110, while 18S rRNA was the least stable.
Nevertheless, EF-1a and eIF-5A were found to be in the forefront
of the rankings. The results of NormFinder and geNorm were
almost consistent.
The stability of the candidate reference gene expression was also
analyzed using BestKeeper, an Excel-based tool. In this analysis,
the average Ct value of every single reaction is applied to analyze
the stability of each candidate reference gene [25]. Rankings of the
candidate reference genes are based on their pair-wise correlation
with this index value, which is indicated by the Pearson correlation
coefficient (r) [35]. BestKeeper calculates the standard deviation
(SD) and the coefficient of variation (CV) based on the Ct values.
The most stable reference genes exhibit the lowest CV and SD
(CV6SD) [39]. Because the maximum number of genes analyzed
by this algorithm is 10 [40], the candidate genes that rank lower in
the previous analyses are generally ruled out. The ranking of the
genes revealed through BestKeeper analysis is presented in
Table 4. These results were mostly consistent with those obtained
Figure 4. Gene expression stability and ranking of 20 candidate reference genes as caluculated by geNorm. The stability value (M) wasdetermined by assessing the mean pairwise variations of all genes; the least stable gene (the highest M value) was excluded, and the M value wasrecalculated until the most stable pair was selected.doi:10.1371/journal.pone.0112177.g004
Figure 5. Determination of the optimal number of reference genes required for effective normalization. The geNorm programcalculated an NF and used the variable V to determine pairwise variation (Vn/Vn+1) between two sequential NFs (NFn and NFn+1). Additional genesare included when V exceeds the cutoff value, which is typically set at 0.15 but is not always achievable. The number of reference genes is deemedoptimal when the lowest possible V value is achieved, at which point it is unnecessary to include additional genes in the normalization strategy.doi:10.1371/journal.pone.0112177.g005
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 8 November 2014 | Volume 9 | Issue 11 | e112177
Ta
ble
3.
Ran
kin
go
fca
nd
idat
ere
fere
nce
ge
ne
sin
ord
er
of
the
ire
xpre
ssio
nst
abili
tyas
calc
ula
ted
by
No
rmFi
nd
er
soft
war
e.
Ra
nk
To
tal
Ro
ot
Ste
mL
ea
fL
PF
SG
FS
RF
SR
GS
1C
YP
GA
PD
HQ
CR
CY
PC
YP
aTU
BaT
UB
aTU
BQ
CR
Mv
alu
e0
.22
70
.01
30
.09
50
.04
00
.11
40
.03
80
.02
10
.01
30
.03
2
2Q
CR
V-A
TP
AR
FQ
CR
60
SR
PL1
3C
YP
V-A
TP
F-b
ox
V-A
TP
Mv
alu
e0
.25
50
.11
00
.10
50
.04
60
.14
60
.06
20
.02
10
.01
30
.03
2
3e
IF-5
A3
0SR
PS2
0EF
-1a
18
SrR
NA
QC
RU
BQ
eIF
-5A
18
SrR
NA
aTU
B
Mv
alu
e0
.26
80
.11
40
.13
00
.07
50
.15
70
.07
00
.06
20
.11
80
.03
3
4EF
-1a
CY
P3
0SR
PS2
0EF
-1a
EF-1a
V-A
TP
QC
Rb
TU
Bp
ol
IIa
Mv
alu
e0
.29
30
.14
80
.13
40
.16
90
.24
30
.09
40
.14
30
.16
50
.09
3
5A
CT
16
-PG
GA
PD
HG
AP
DH
eIF
-5A
eIF
-5A
EF-1a
eIF
-5A
bT
UB
Mv
alu
e0
.30
10
.15
80
.18
60
.18
30
.28
20
.10
50
.14
70
.19
50
.10
9
6V
-AT
PSA
R1
6-P
GA
RF
V-A
TP
EF-1a
po
lIIa
SAR
1G
AP
DH
Mv
alu
e0
.30
60
.18
20
.19
00
.19
10
.30
30
.18
20
.20
40
.22
20
.12
5
7G
AP
DH
EF-1a
60
SRP
L13
SAR
1A
RF
18
SrR
NA
AC
T1
GA
PD
HSA
R1
Mv
alu
e0
.34
20
.20
50
.25
60
.21
00
.31
20
.18
80
.24
30
.30
00
.20
9
8A
RF
eIF
-5A
CY
PC
DP
18
SrR
NA
QC
R6
-PG
QC
RA
RF
Mv
alu
e0
.36
20
.21
10
.25
80
.22
60
.37
00
.27
30
.29
10
.31
40
.23
7
9U
BQ
AR
FA
CT
13
0SR
PS2
0SA
R1
AC
T1
CY
PC
YP
60
SRP
L13
Mv
alu
e0
.36
80
.21
60
.28
30
.23
20
.40
60
.32
00
.30
70
.33
90
.25
8
10
po
lIIa
F-b
ox
V-A
TP
AC
T1
AC
T1
AR
FU
BQ
AC
T1
UB
Q
Mv
alu
e0
.38
30
.22
40
.28
70
.25
90
.40
80
.32
90
.35
20
.35
10
.30
0
11
SAR
1b
TU
BU
BQ
eIF
-5A
aTU
Bp
ol
IIab
TU
BEF
-1a
CY
P
Mv
alu
e0
.38
60
.25
90
.29
00
.26
00
.46
70
.37
50
.39
80
.36
80
.30
3
12
aTU
BaT
UB
po
lIIa
TC
TP
po
lIIa
bT
UB
F-b
ox
60
SRP
L13
EF-1a
Mv
alu
e0
.40
50
.27
70
.32
00
.30
10
.50
30
.39
10
.44
10
.39
10
.30
8
13
bT
UB
UB
Qe
IF-5
A6
-PG
GA
PD
HG
AP
DH
GA
PD
HV
-AT
PC
DP
Mv
alu
e0
.42
70
.34
40
.35
20
.31
20
.52
50
.40
50
.49
10
.40
50
.33
6
14
F-b
ox
AC
T1
SAR
1aT
UB
6-P
G6
-PG
AR
FU
BQ
AC
T1
Mv
alu
e0
.44
30
.34
40
.35
20
.33
60
.59
30
.45
60
.53
40
.50
60
.36
7
15
18
SrR
NA
QC
RF-
bo
xp
ol
IIab
TU
BF-
bo
xT
CT
PA
RF
F-b
ox
Mv
alu
e0
.45
30
.35
10
.41
20
.33
80
.63
40
.50
90
.59
90
.56
80
.37
1
16
60
SRP
L13
po
lIIa
aTU
Bb
TU
BU
BQ
60
SRP
L13
SAR
1p
ol
IIa1
8Sr
RN
A
Mv
alu
e0
.47
60
.38
60
.54
20
.41
80
.66
30
.52
80
.60
90
.71
50
.42
6
17
6-P
GC
DP
CD
PF-
bo
xC
DP
SAR
11
8Sr
RN
A6
-PG
TC
TP
Mv
alu
e0
.51
90
.41
10
.55
20
.43
20
.67
80
.54
70
.67
10
.71
90
.55
0
18
CD
PT
CT
P1
8Sr
RN
A6
0SR
PL1
3F-
bo
xC
DP
60
SRP
L13
CD
Pe
IF-5
A
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 9 November 2014 | Volume 9 | Issue 11 | e112177
using geNorm, including the total samples, roots, stems, leaves and
LP.
In summary, CYP and EF-1a were demonstrated to be the best
reference genes under all the treatment conditions. In addition,
GAPDH and V-ATP showed the highest CV6SD values
(0.3760.08 and 0.4660.11, respectively) in the roots. However,
ACT1 and QCR were the most stable reference genes in FS, and
their CV6SD values were 0.2960.06 and 0.4460.11, respective-
ly, which slightly differed between geNorm and NormFinder.
Discussion
Selection of suitable reference genes is a crucial pre-condition to
a successful gene expression study based on qRT-PCR. Using
inaccurate reference genes can lead to conflicting results,
particularly when the variations in the rate of transcription
between sample groups are small [41]. Herein, we have described
a systematic analysis involving the stability of mRNA expression of
candidate genes for data normalization in qPCR experiments
using different developmental stages and the three vegetative
organs of Panax ginseng. Investigation of 20 candidate reference
genes by Ct value, geNorm, NormFinder, and BestKeeper applets
led to the identification of the best reference genes for differential
gene expression analyses at different developmental stages and
various organs of ginseng. In qRT-PCR analysis, certain
housekeeping genes (such as, ACT, UBQ, F-box) are considered
stably expressed in different environmental conditions and are
commonly employed as reference gene(s) [36]. The analysis data
revealed certain changes in the mRNA gene expression levels in
majority of the traditional housekeeping genes of ginseng under
different treatment conditions; therefore, these genes could not be
considered as ideal ginseng reference genes. However, a stable
reference gene is essential for genetic engineering studies in
ginseng. To the best of our knowledge, this is the first report on the
identification and validation of suitable reference genes for qRT-
PCR analysis of ginseng.
An ‘‘ideal’’ reference gene(s) should be continually transcribed
in all cell types and organs. Additionally, its RNA transcription
level should be relatively constant in response to the internal and
external stimulations [42]. For example, during housekeeping gene
selection for qRT-PCR normalization in potato, it was found that
the expression of EF-1a was not influenced by cold, salt, or late
blight stressors [29]. In the analysis of reference genes for
Arabidopsis, EF-1a was relatively stable in different organs [43].
However, under nutrition deficiency or abiotic stress, the stability
of EF-1a was poor [44]. Selected as the appropriate reference gene
in cucumber, the CYP gene was the most stable gene under cold
and heat stress treatments; nevertheless it was less stable in various
other tissues [45]. Based on our statistical analyses using Ct value,
geNorm, NormFinder, and BestKeeper applets, the mRNA
expression level of CYP, a traditional housekeeping gene, was
found to be the most stable in different organs and developmental
stages, and was followed by EF-1a (Table 5). Furthermore, out of
the 10 novel reference genes, it was interesting to note that QCR
was relatively stable in all the experimental samples.
Although the results of all the three applets were reasonable,
they were not found to be completely consistent. However, this
variation was not surprising, since the three software applications
are based on different calculation algorithms [25]. geNorm is
known to be a more effective and feasible algorithm for ensuring
the optimal stability of reference genes, whereas NormFinder and
BestKeeper are best applied for assessing the quality of the gene
rankings obtained by geNorm [26,46,47]. The results of the
geNorm analysis have been satisfactorily accepted by many
Ta
ble
3.
Co
nt.
Ra
nk
To
tal
Ro
ot
Ste
mL
ea
fL
PF
SG
FS
RF
SR
GS
Mv
alu
e0
.63
50
.46
70
.56
60
.43
60
.69
50
.92
80
.76
50
.78
00
.57
4
19
TC
TP
60
SRP
L13
bT
UB
V-A
TP
TC
TP
TC
TP
CD
PT
CT
P6
-PG
Mv
alu
e0
.67
20
.50
60
.56
80
.46
10
.94
91
.16
60
.79
50
.83
20
.66
5
20
30
SRP
S20
18
SrR
NA
TC
TP
UB
Q3
0S
RP
S20
30
SRP
S20
30
SRP
S20
30
SRP
S20
30
SRP
S20
Mv
alu
e1
.07
00
.59
10
.81
30
.48
41
.21
61
.24
01
.24
21
.33
71
.16
6
No
tes:
LP,
leaf
-exp
ansi
on
pe
rio
d;
FS,
the
flo
we
rst
age
;G
FS,
the
gre
en
fru
itst
age
;R
FS,
the
red
fru
itst
age
s;R
GS,
the
roo
tg
row
ing
afte
rfr
uit
stag
e.
do
i:10
.13
71
/jo
urn
al.p
on
e.0
11
21
77
.t0
03
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 10 November 2014 | Volume 9 | Issue 11 | e112177
Ta
ble
4.
Ran
kin
go
fca
nd
idat
ere
fere
nce
ge
ne
sin
ord
er
of
the
ire
xpre
ssio
nst
abili
tyas
calc
ula
ted
by
Be
stK
ee
pe
rso
ftw
are
.
Ra
nk
To
tal
Ro
ot
Ste
mL
ea
fL
PF
SG
FS
RF
SR
GS
1C
YP
GA
PD
HC
YP
CY
P6
0SR
PS2
0A
CT
1C
YP
CY
PC
YP
CV
%±
SD
0.8
56
0.1
80
.376
0.0
80
.716
0.1
50
.536
0.1
10
.416
0.0
90
.296
0.0
61
.046
0.2
20
.026
0.0
00
.636
0.1
3
2EF
-1a
V-A
TP
V-A
TP
QC
REF
-1a
QC
REF
-1a
EF-1a
EF-1a
CV
%±
SD
1.0
76
0.2
60
.466
0.1
10
.736
0.1
70
.576
0.1
50
.436
0.1
00
.446
0.1
11
.306
0.3
30
.406
0.1
00
.656
0.1
6
3e
IF-5
AC
YP
EF-1a
EF-1a
CY
PA
RF
po
lIIa
SAR
1b
TU
B
CV
%±
SD
1.7
56
0.3
60
.476
0.1
00
.736
0.1
70
.576
0.1
40
.636
0.1
30
.446
0.1
11
.716
0.4
30
.806
0.1
91
.116
0.2
8
4A
RF
30
SRP
S20
AR
FSA
R1
QC
REF
-1a
V-A
TP
ase
eIF
-5A
GA
PD
H
CV
%±
SD
1.9
06
0.4
50
.616
0.1
51
.126
0.2
71
.176
0.2
90
.796
0.2
00
.756
0.1
91
.786
0.4
10
.846
0.1
71
.456
0.3
2
5SA
R1
EF-1a
AC
T1
18
SrR
NA
AC
T1
V-A
TP
aTU
Bb
TU
BA
RF
CV
%±
SD
1.9
86
0.4
70
.656
0.1
61
.196
0.4
01
.276
0.2
21
.826
0.3
71
.006
0.2
31
.906
0.4
61
.106
0.2
91
.616
0.3
8
6Q
CR
SAR
16
-PG
AR
F1
8Sr
RN
AC
YP
eIF
-5A
18
SrR
NA
QC
R
CV
%±
SD
1.9
96
0.5
11
.016
0.2
41
.256
0.3
01
.276
0.3
12
.246
0.3
81
.136
0.2
42
.036
0.4
21
.446
0.2
41
.656
0.4
2
7V
-AT
Pe
IF-5
AG
AP
DH
30
SRP
S20
V-A
TP
eIF
-5A
bT
UB
GA
PD
HaT
UB
CV
%±
SD
2.2
06
0.5
21
.056
0.2
11
.346
0.3
11
.716
0.3
72
.326
0.5
41
.186
0.2
42
.286
0.5
91
.466
0.3
21
.696
0.4
0
8aT
UB
6-P
GQ
CR
GA
PD
HA
RF
UB
QQ
CR
aTU
Bp
ol
IIa
CV
%±
SD
2.3
66
0.5
61
.176
0.2
71
.46
0.3
61
.866
0.4
12
.436
0.5
71
.776
0.3
92
.626
0.6
71
.826
0.4
51
.746
0.4
5
9G
AP
DH
F-b
ox
60
SRP
L13
eIF
-5A
eIF
-5A
aTU
BA
CT
1F-
bo
xV
-AT
P
CV
%±
SD
2.3
96
0.5
31
.246
0.2
81
.926
0.4
82
.086
0.4
42
.506
0.5
21
.876
0.4
52
.916
0.6
21
.996
0.4
52
.096
0.4
9
10
AC
T1
AR
F3
0SR
PS2
0C
DP
SAR
11
8Sr
RN
A6
-PG
QC
RSA
R1
CV
%±
SD
2.5
16
0.5
21
.316
0.3
11
.926
0.4
82
.216
0.5
32
.516
0.6
03
.016
0.5
03
.066
0.7
32
.276
0.5
62
.586
0.6
2
No
tes:
LP,l
eaf
-exp
ansi
on
pe
rio
d;F
S,th
efl
ow
er
stag
e;G
FS,t
he
gre
en
fru
itst
age
;RFS
,th
ere
dfr
uit
stag
es;
RG
S,th
ero
ot
gro
win
gaf
ter
fru
itst
age
.De
scri
pti
vest
atis
tics
of
10
can
did
ate
ge
ne
sb
ase
do
nth
eco
eff
icie
nt
of
vari
ance
(CV
)an
dst
and
ard
de
viat
ion
(SD
)o
fth
eir
Ct
valu
es
we
red
ete
rmin
ed
usi
ng
the
wh
ole
dat
ase
t.R
efe
ren
ceg
en
es
we
reid
en
tifi
ed
asth
em
ost
stab
leg
en
es,
i.e.t
ho
sew
ith
the
low
est
coe
ffic
ien
to
fva
rian
cean
dst
and
ard
de
viat
ion
(CV
%6
SD).
do
i:10
.13
71
/jo
urn
al.p
on
e.0
11
21
77
.t0
04
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 11 November 2014 | Volume 9 | Issue 11 | e112177
Ta
ble
5.
Stab
ility
ran
kin
go
f2
0ca
nd
idat
ere
fere
nce
ge
ne
su
sin
gg
eN
orm
,No
rmfi
nd
er
and
Be
stke
ep
er.
To
tal
Ro
ot
Ste
mL
ea
fL
PF
SG
FS
RF
SR
GS
GN
BG
NB
GN
BG
NB
GN
BG
NB
GN
BG
NB
GN
B
CY
P1
11
54
31
81
11
11
13
12
69
91
49
17
11
1
EF
-1a
14
29
75
83
34
43
37
23
64
41
12
51
12
81
22
QC
R4
26
15
15
—*
61
71
22
42
44
82
68
81
28
10
51
6
eIF
-5A
33
34
87
14
13
—1
01
19
10
49
15
73
56
15
41
81
8—
AC
T1
75
10
16
14
—9
98
12
10
—6
10
56
91
71
09
13
10
—1
41
4—
V-A
TP
86
73
22
71
02
17
19
—1
15
75
45
11
34
14
13
—6
29
GA
PD
H1
17
91
11
45
65
58
17
14
—1
31
3—
17
7—
97
73
64
AR
F5
84
69
10
52
49
66
96
87
10
31
11
5—
15
15
—9
85
UB
Q1
29
—8
13
—1
31
1—
20
20
—1
41
5—
83
81
51
4—
10
14
—1
01
0—
po
lII
a1
41
0—
14
16
—1
21
2—
13
15
—1
31
2—
14
11
—5
16
31
71
6—
14
8
SA
R1
61
15
12
66
11
14
—7
74
88
10
12
17
—1
36
—1
63
11
71
0
aT
UB
91
28
10
12
—1
61
3—
14
14
—7
11
—9
19
11
57
18
13
7
bT
UB
10
13
—1
11
1—
17
19
—1
61
6—
15
18
—1
11
2—
16
47
64
54
53
F-b
ox
15
14
—7
10
91
51
5—
19
7—
15
17
—1
61
5—
10
2—
82
91
51
5—
18
SrR
NA
13
15
—2
02
0—
18
18
—3
35
59
61
07
10
14
3—
33
61
61
6—
60
SR
PL
13
17
16
—1
91
9—
17
91
81
8—
13
11
71
6—
18
12
—1
11
9—
12
9—
6-P
G1
61
7—
13
58
36
51
11
3—
12
13
—1
51
4—
81
71
01
61
2—
19
19
—
CD
P1
91
8—
17
17
—1
91
7—
88
10
16
16
—1
81
8—
19
18
—1
91
7—
13
13
—
TC
TP
18
19
—1
81
8—
20
20
—1
51
2—
19
19
—1
91
9—
12
19
—1
81
8—
17
17
—
30
SR
PS
20
20
20
—1
34
10
41
06
97
20
20
—2
02
0—
20
20
—2
02
0—
20
20
—
No
tes:
LP,l
eaf
-exp
ansi
on
pe
rio
d;F
S,th
efl
ow
er
stag
e;G
FS,t
he
gre
en
fru
itst
age
;RFS
,th
ere
dfr
uit
stag
es;
RG
S,th
ero
ot
gro
win
gaf
ter
fru
itst
age
.G,g
eN
orm
soft
war
e;N
,No
rmfi
nd
er
soft
war
e;B
,Be
stke
ep
er
soft
war
e.*
me
ans
Ith
asn
ot
be
en
test
ing
by
Be
stke
ep
er.
do
i:10
.13
71
/jo
urn
al.p
on
e.0
11
21
77
.t0
05
Reference Gene Selection in Panax ginseng
PLOS ONE | www.plosone.org 12 November 2014 | Volume 9 | Issue 11 | e112177
researchers [21-26,48,49]. In the present study, the two top ranked
reference genes for the total samples, roots, leaves, and the
developmental stage, LP, obtained through geNorm were consis-
tent with the ranking of NormFinder and BestKeeper. However,
the two best ranked reference genes in the stems and other
developmental stages (FS, GFS, RFS, and RGS), as analyzed by
geNorm, were slightly different from the results produced by
NormFinder or BestKeeper; interestingly, the genes were still top-
ranked. Our data showed that CYP and EF-1a were the most
stable reference genes among all the samples. Meanwhile, different
types of samples revealed their own best reference genes amongst
the 20 selected candidate reference genes. In the different
vegetative organs of ginseng, GAPDH and 30SRPS20 were the
best reference genes found in the roots; CYP and 60SRPL13 were
the top-ranked reference genes in the stems; and CYP and QCR
were the best reference genes in the leaves. In different
developmental stages of ginseng, CYP/60SRPL13, CYP/eIF-5A,
aTUB/V-ATP, eIF-5A/SAR1, and aTUB/pol IIa were the most
stably expressed combinations in LP, FS, GFS, RFS, and RGS,
respectively. Their CV and MFC values were relatively low.
Although 30SRPS20 was the least stable among the 20 candidate
reference genes in all five developmental stages, it ranked high in
the roots, as determined by geNorm, NormFinder, and Best-
Keeper.
Taken together, we identified 20 potential reference genes from
15 P. ginseng samples (different organs and developmental stages)
for the normalization of qRT-PCR data. CYP and EF-1a were the
most suitable reference genes in ginseng, as evaluated by the three
software applications.
Conclusion
Gene transcription studies using real-time quantitative reverse
and multiple actions. Biochem Pharmacol 58(11): 1685–1693.
5. Hou JP (1977) The chemical constituents of ginseng plants. Comp Med East and
West 5: 123–145.
6. Yip TT, Lau CN, But PP, Kong YC (1985) Quantitative analysis of ginsenosides
in fresh Panax ginseng. Am J Chin Med 13: 77–88.
7. Wang HW, Peng DC, Xie JT(2009)Ginseng leaf-stem: Bioactive constituents
and pharmacological functions. Chin Med 4: 20.
8. Jackson CJC, Dini JP, Lavandier C, Faulkner H (2003) Ginsenoside content of
North American ginseng (Panax quinquefolius L. Araliaceae) in relation to plantdevelopment and growing locations. J Ginseng Research 27: 135–140.
9. Chen WY, Chen CB, Li YB (2010) Ginseng development practical manual.
Changchun: Jilin Science and Technology Bureau.
10. Sathiyaraj G, Srinivasan S, Subramanium S (2010) Polygalacturonase inhibitingprotein: isolation, developmental regulation and pathogen related expression in