Next Article in Journal
Altered Fatty Acid Metabolism-Related Gene Expression in Liver from Morbidly Obese Women with Non-Alcoholic Fatty Liver Disease
Previous Article in Journal
Comparing Biofouling Control Treatments for Use on Aquaculture Nets
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 15
Issue 12
10.3390/ijms151222155
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Identification and Evaluation of Reliable Reference Genes for Quantitative Real-Time PCR Analysis in Tea Plant (Camellia sinensis (L.) O. Kuntze)

by
Xinyuan Hao
1,2,3,4,
David P. Horvath
2,
Wun S. Chao
2,
Yajun Yang
1,3,4,*,
Xinchao Wang
3,4 and
Bin Xiao
1
1
College of Horticulture, Northwest A&F University, Yangling 712100, Shaanxi, China
2
USDA-Agricultural Research Service, Biosciences Research Lab, Sunflower and Plant Biology Research Unit, 1605 Albrecht Blvd N, Fargo, ND 58102, USA
3
Tea Research Institute, Chinese Academy of Agricultural Sciences, No. 9 South Meiling Road, Hangzhou 310008, Zhejiang, China
4
National Center for Tea Improvement, Key Laboratory of Tea Biology and Resources Utilization, Ministry of Agriculture, No. 9 South Meiling Road, Hangzhou 310008, Zhejiang, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2014, 15(12), 22155-22172; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms151222155
Submission received: 12 October 2014 / Revised: 24 November 2014 / Accepted: 24 November 2014 / Published: 2 December 2014
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
Reliable reference selection for the accurate quantification of gene expression under various experimental conditions is a crucial step in qRT-PCR normalization. To date, only a few housekeeping genes have been identified and used as reference genes in tea plant. The validity of those reference genes are not clear since their expression stabilities have not been rigorously examined. To identify more appropriate reference genes for qRT-PCR studies on tea plant, we examined the expression stability of 11 candidate reference genes from three different sources: the orthologs of Arabidopsis traditional reference genes and stably expressed genes identified from whole-genome GeneChip studies, together with three housekeeping gene commonly used in tea plant research. We evaluated the transcript levels of these genes in 94 experimental samples. The expression stabilities of these 11 genes were ranked using four different computation programs including geNorm, Normfinder, BestKeeper, and the comparative ∆CT method. Results showed that the three commonly used housekeeping genes of CsTUBULIN1, CsACINT1 and Cs18S rRNA1 together with CsUBQ1 were the most unstable genes in all sample ranking order. However, CsPTB1, CsEF1, CsSAND1, CsCLATHRIN1 and CsUBC1 were the top five appropriate reference genes for qRT-PCR analysis in complex experimental conditions.

1. Introduction

Gene expression analysis is an important approach, providing insight into the genetic and developmental mechanisms in biological research. Quantitative real-time polymerase chain reaction (qRT-PCR) has become a critical tool for rapid and reliable quantification of gene transcript levels due to its simplicity, specificity, and sensitivity [1,2,3]. However, a successful qRT-PCR assay relies on many factors, including the integrity of purified RNA, reverse transcription efficiency, primer design, detection chemistry selection, and data analysis [4,5]. Because experimental variation is unavoidable, an accurate method of normalization is essential in qRT-PCR assays [6]. Normalizing results with a stably-expressed internal reference gene is a simple and popular method for controlling error in qRT-PCR, and is generally preferred over other normalization strategies based on sample size, total RNA quantification, and spiked control molecules [7]. Therefore, the selection and validation of endogenous reference genes for qRT-PCR in each species is necessary.
Often actin, ubiquitin, or 18S rRNA, or other housekeeping genes are used without any evidence that such genes are stably expressed across the experimental conditions described. These genes were reasonable choices for hybridization-based expression analyses where the probes for these housekeeping genes often hybridized to multiple gene family members and thus were representative of the overall gene expression. However, we now know from RNAseq studies that individual gene family members can have very different expression values in regards to both tissue specificity and response to environmental changes. Because qRT-PCR is often gene specific, it is no longer acceptable to arbitrarily choose a common housekeeping gene as an internal control unless there is evidence that it is likely to show stable expression across the samples of any given experiment. Indeed, several published results examining the suitability of various commonly used control genes indicate substantial variation in expression of commonly used control genes, but have identified several genes that appear to be more appropriate for use as controls for qRT-PCR [8,9]. However, an examination of such results suggests that the best of control genes can vary between species and even between experimental conditions. Thus, an evidence-based choice of control genes is needed for any given gene expression analysis using qRT-PCR.
Tea is the most consumed beverage in the world aside from water [10], and has many beneficial health effects [11,12,13]. Tea plant (Camellia sinensis), is an important cash crop, and is grown commercially in about 30 different countries. To improve the quality and productivity of tea plant, an increasing number of researchers are investigating the physiology, biochemistry and metabolism of tea plant using molecular biotechnology [14]. Singh et al. [15] evaluated 26S rRNA genes by non-quantitative RT-PCR in tea plant together with six other plant species. Gohain et al. [16] validated the suitability of primers that amplify 18S rRNAs, 26s rRNAs, ribulose-1,5-bisphosphates (RuBPs) and actins as potential reference genes, but only examined transcript levels in tea leaves. To date, no comprehensive study on endogenous reference genes for qRT-PCR analysis has been done in tea plant. Surprisingly, various members of the housekeeping gene families encoding actins, tubulins, glyceraldehyde 3-phosphate dehydrogenases (GAPDHs), 18S rRNAs and 26S rRNAs are the only reference genes currently used in gene expression analysis of tea plant [17,18,19,20,21,22,23,24]. Due to the increased sensitivity and dynamic range of real-time PCR over traditional quantization techniques, an increasing number of studies have documented the poor stability of commonly known housekeeping genes like β-actins, GAPDHs and 18S rRNAs [25,26]. To avoid erroneous results, the expression stability of chosen reference gene needs to be validated depending on individual samples/experiments rather than relying on previously published materials unless validation of the chosen reference genes was demonstrated [6]. Recently, studies on response of tea plant to cold, seasonal dormancy, and tissue specific gene expression have become subjects of interest [17,19,22,23]. Thus, the validation of highly stable reference genes under these conditions is imperative.
To identify suitable reference genes for normalizing qRT-PCR tea plant gene expression analysis specifically suitable for the expression analysis in several different experimental conditions, we examined the expression stability of 11 candidate reference genes from three different sources: the orthologs of traditional Arabidopsis reference gene family members encoding GAPDHs, elongation factors (EFs), ubiquitins (UBQs) and ubiquitin carrier proteins (UBCs), stably expressed genes identified from whole-genome GeneChip studies that encode clathrins, MONENSIN SENSITIVITY1 (MON1 or SAND), tonoplast intrinsic proteins (TIPs) and polypyrimidine tract-binding proteins (PTBs) [27], and three housekeeping gene currently used in tea plant research that encode 18S rRNAs, actins and tubulins. We evaluated the transcript levels of these genes together with two tea plant target genes SHORT VEGETATIVE PHASE 1 (CsSVP1) and FLOWERING LOCUS T 1 (CsFT1) in 94 experimental samples consisting of diurnal variety (24 samples), different organs (18 samples), different duration of cold stress in leaves and shoots (each contains 18 samples), and different time points after auxin antagonist auxinole treatment (eight samples) and lanolin treatment (eight samples). Accurate quantification of cDNA synthesis products was performed using radioactivity detection of [α-32P]-dCTP incorporation. The expression stability of candidate reference genes was ranked using comparative ∆CT, BestKeeper, NormFinder and geNorm software. We determined that genes encoding 18S rRNAs, actins and tubulins, currently used in tea plant studies were the most unstable reference genes, while specific genes CsPTB1, CsEF1, CsSAND1, CsCLATHRIN1 and CsUBC1 were the most appropriate reference genes for qRT-PCR analysis in tea plant under complex experimental conditions.

2. Results

To identify stable reference genes for gene expression studies in tea plant, 11 candidate reference genes and two target genes were investigated by qRT-PCR using six series of samples. First, the primer specificity was detected based on melting curve analysis and agarose gel electrophoresis analysis using pooled cDNAs as template. Results showed that at 53 °C annealing temperature, all primer pairs amplified a single major peak with all samples. Only CsUBQ1 contained a second minor shoulder peak at 53 °C or lower temperatures (Supplemental Figure S1). After 45 cycles of amplification, the PCR products were examined on a 1.5% agarose gel, and all the tested amplicons produced a single visible band of the expected size (Supplemental Figure S2). We also checked the primer specificity performance at a wide range of annealing temperature from 50 to 59 °C with intervals of 3 °C. Melting curve analysis indicated that all the designed primers worked well at all of the tested annealing temperature, except for CsUBQ1 which was not specific when amplified at an annealing temperature of 53 °C (Supplemental Table S1). The efficiency of PCR amplification, defined as percentage of PCR product increase per cycle, was determined for each gene using two-fold serial dilution of the pooled cDNAs (Table 1). Resulted showed that our primers which amplified Cs18S rRNA1, CsACTIN1, CsCLATHRIN1, CsEF1, CsSAND1, CsTIP1 and CsUBQ1 had excellent amplification efficiency of between 98% and 106%, while CsGAPDH1, CsTUBULIN1, CsUBC1 and CsPTB1 had less but still acceptable efficiencies of between 92% and 95%.
Table
Table 1. Description of tea plant reference genes for qRT-PCR.
Table 1. Description of tea plant reference genes for qRT-PCR.
Reference GeneGeneBank Accession NumberArabidopsis OrthologArabidopsis Locus DescriptionForward/Reverse Primer Sequence (5'–3')Amplicon Size (bp)qRT-PCR Efficiency (%)
Cs18S rRNA1AY563528.1At3g4176818S ribosomal RNATGGCCTTCGGGATCGGAGTAAT/GCTTTCGCAGTTGTTCGTCTTTCAT108101
CsACTIN1KA280216.1At5g09810Actin 2/7TGGGCCAGAAAGATGCTTATGTAGG/ATGCCAGATCTTTTCCATGTCATCC118103
CsCLATHRIN1KA291473.1AT5G46630Clathrin adaptor complex subunitTAGAGCGGGTAGTGGAGACCTCGTT/TACCAAAGCCGGCTCGTATGAGATT129102
CsEF1KA280301.1AT5G60390Elongation factor 1 alphaTTGGACAAGCTCAAGGCTGAACG/ATGGCCAGGAGCATCAATGACAGT11098
CsGAPDH1KA295375.1AT1G13440Glyceraldehyde-3-phosphate dehydrogenaseTTTTTGGCCTTAGGAACCCAGAGG/GGGCAGCAGCCTTATCCTTATCAGT10793
CsSAND1KM057790AT2G28390SAND family proteinTCCAATTGCCCCCTTAATGACTCA/GTAAGGGCAGGCAAACACCAGGTA109106
CsTIP1KA283116.1AT4G34270TIP41-like proteinGGGTGCTTATGAGATTGAGGGACAC/ATATGCCGCAGAATCAGATGGGTAT148101
CsTUBULIN1DQ444294.1XM_002871860Tubulin alpha-3TGCCGCTAATAACTTTGCCAGAGG/GCCACCGCCAACAGCGTTAGTA14192
CsUBQ1HM003234.1AT4G05320UbiquitinGCCGGAAAACAGCTTGAAGATGG/AGGACGGCTCAATAATCCCACCAC11198
CsUBC1KA281185.1AT4G27960Ubiquitin-conjugating enzymeTGCTGGTGGGGTTTTTCTTGTTACC/AAGGCATATGCTCCCATTGCTGTTT12492
CsPTB1GAAC01052498.1AT3G01150Polypyrimidine tract-binding proteinTGACCAAGCACACTCCACACTATCG/TGCCCCCTTATCATCATCCACAA10795

2.1. Expression Profiles of Candidate Reference Genes

Based on radioactive tracer normalized and quantified cDNA template, gene expression analysis of the 11 candidate reference gene exhibited a narrow CT range among all the experimental series (Figure 1). The CT value ranged from 13.1 to 24.4, while the majority of these values were distributed between 20 and 23. The transcript of Cs18S rRNA1 was the most abundant (with lowest CT value of 13.1) housekeeping gene, while the CsSAND1 gene produced the least abundant transcript with a CT value of 24.4. Each reference genes had varied expression ranges across all studied cDNA samples. CsPTB1, CsSAND1 and CsUBC1 had the lowest transcript variations among studied reference genes, while CsUBQ1, CsACTIN1 and CsTUBULIN1 showed the highest expression variations. For target genes, CsSVP1 gene had a mean CT value of 24.4, while CsFT1 had a very high CT value of 36.6. CsFT1 gene also had a more dynamic range of transcript accumulation compared to SVP1 and the various reference genes. Most of the reference genes studied had smaller expression variation compared to the variation of target genes.
Ijms 15 22155 g001 1024
Figure 1. Average cycle threshold (CT) values for 11 candidate reference genes. The filled diamond symbol indicates the mean of CT values. The bars indicate standard deviation.
Figure 1. Average cycle threshold (CT) values for 11 candidate reference genes. The filled diamond symbol indicates the mean of CT values. The bars indicate standard deviation.
Ijms 15 22155 g001
Expression profiles of candidate reference genes and target genes are shown in Figure 2. In each series of experiments, the first sampling time point or organ was chosen as the control. Differences in gene expression relative to the chosen control were presented as ratios of log2 transformed relative expression values. Results showed that the target gene of CsFT1 had significant expression fluctuation among the six experimental series, while another target gene of CsSVP1 had relatively stable expression. Except for CsACTIN1 and CsTUBULIN1, most of the candidate reference genes had relatively more stable expression in different experimental series compared with the target genes. During the diurnal expression in leaves, CsTUBULIN1was significantly down-regulated at 9:00 am and up-regulated at 6:00 pm, while all of the other candidate reference genes had relatively stable expression. CsACTIN1 and CsTUBULIN1 had clear differential transcript abundance in different organs, but the expressions of CsSAND1, CsTIP1 and CsPTB1 were consistent across the different organs tested. CsACTIN1 and CsTUBULIN1 were highly induced in leaves after cold treatment. CsUBQ1, CsCLATHRIN1 and CsTIP1 also had significant expression variations at various time points during the cold treatment. In shoots, no significant expression fluctuation was observed with most reference genes at 1 day of cold and short photoperiod treatment, except for the substantial repression of CsUBQ1, CsCLATHRIN1 and CsTIP1, and the observed increase in accumulation of the CsTUBULIN1 transcript. Following auxinole (an auxin action inhibitor) and lanolin treatment, the expression of CsACTIN1, CsGAPDH1, CsUBQ1, CsTIP1 and Cs18S rRNA1 varied substantially in shoots. Among all the candidate reference genes, CsACTIN1 and CsTUBULIN1 were the most unstable, while CsPTB1 seemed to be most stable in all of the experimental series.
Ijms 15 22155 g002 1024
Figure 2. Expression profile of candidate reference genes and target genes under different experimental conditions. Six experimental series including 94 samples were used to examine gene expression. The fold difference is designated as log2 value. Red indicates up-regulated genes and green indicates down-regulated genes as compared with first sampling point/tissue. Bars at the bottom indicate the range of transcript changes in log2 value. The range of transcript changes for each experimental series is provided inside the parenthesis.
Figure 2. Expression profile of candidate reference genes and target genes under different experimental conditions. Six experimental series including 94 samples were used to examine gene expression. The fold difference is designated as log2 value. Red indicates up-regulated genes and green indicates down-regulated genes as compared with first sampling point/tissue. Bars at the bottom indicate the range of transcript changes in log2 value. The range of transcript changes for each experimental series is provided inside the parenthesis.
Ijms 15 22155 g002

2.2. Stability Ranking of Candidate Reference Genes

Four different computation programs geNorm, Normfinder, BestKeeper, and comparative ∆CT were used to rank the candidate reference genes based on CT values of 94 samples (Supplemental Table S2) from different experimental series. The stability ranking results are shown in Table 2. Additionally, the analytical tool RefFinder also gave a recommended comprehensive ranking for each gene based on weight assignment (Table 2, details were presented in Supplemental Figure S3). Based on this analysis, CsEF1, CsGAPDH1, CsCLATHRIN1 and CsSAND1 were the top four reference genes during the diurnal expression analysis. In contrast, CsACTIN1 and CsTUBULIN1 were the worst reference genes. In different organs, CsCLATHRIN1, CsPTB1 and CsSAND1 were most stable, while CsUBQ1, CsACTIN1 and CsTUBULIN1 were the most unstable genes. In leaves with cold and short photoperiod treatment, the three most stably expressed genes were CsCLATHRIN1, CsGAPDH1 and CsUBQ1 and the four least stable were Cs18S rRNA1, CsACTIN1, CsTIP1 and CsTUBULIN1. Likewise, in shoots with cold and short photoperiod treatment, CsUBC1, CsPTB1, CsGAPDH1 and CsCLATHRIN1 were the best reference genes and had very similar scores. On the contrary, CsACTIN1 and CsTUBULIN1 were again the least stably expressed. After auxinole treatment, CsPTB1, CsCLATHRIN1 and CsUBC1 expression were most stable and CsTIP1, CsUBQ1, CsGAPDH1, Cs18S rRNA1 and CsACTIN1 were relatively less stably expressed. In shoots with lanolin treatment, CsUBC1, CsPTB1, CsCLATHRIN1, and CsEF1 showed good expression stability, while CsGAPDH1, Cs18S rRNA1, CsUBQ1 and CsACTIN1 were not as stable. When all the sample data was included, the order of top ranked reference genes are CsPTB1, CsEF1, CsSAND1, CsCLATHRIN1, and CsUBC1 and the order of worst ranked genes are CsACTIN1, CsTUBULIN1, Cs18S rRNA1 and CsUBQ1.

3. Discussion

In this study, 11 candidate reference genes were selected from currently used reference genes in tea research, orthologs of Arabidopsis genes including superior traditional reference genes, and those identified as stably expressed in Arabidopsis whole-genome GeneChip studies [16,27]. The stability of these candidate reference genes were evaluated by qRT-PCR based on relative cDNA quantification, which is a reliable method to evaluate reference gene stability [28]. Results indicated that the primers had very good specificity through a wide range of annealing temperatures (Supplemental Figures S1 and S2, Table S1). In addition, qRT-PCR for each tested reference gene showed high efficiency (Table 1). These observations suggest our primer design for analyses of reference gene stability in this study is reliable, and, with the exception of the primers for CsUBQ1, were likely gene-specific.
Table
Table 2. Stability ranking of candidate reference genes.
Table 2. Stability ranking of candidate reference genes.
Analysis ToolRanking Order (1 being the most stable, 11 being the least stable)
1234567891011
Diurnal Expression Ranking Order
∆CTCsEF1CsGAPDH1CsSAND1CsCLATHRIN1CsTIP1CsUBQ1CsUBC1CsPTB1Cs18S rRNA1CsACTIN1CsACTIN1
BestKeeperCs18S rRNA1CsPTB1CsEF1CsSAND1CsTIP1CsGAPDH1CsUBC1CsCLATHRIN1CsACTIN1CsUBQ1CsACTIN1
NormfinderCsEF1CsGAPDH1CsUBQ1CsCLATHRIN1CsSAND1CsTIP1CsPTB1CsUBC1Cs18S rRNA1CsACTIN1CsACTIN1
GenormCsGAPDH1|CsCLATHRIN1CsSAND1CsEF1CsTIP1CsUBC1CsUBQ1CsPTB1Cs18S rRNA1CsACTIN1CsACTIN1
Recommended comprehensive rankingCsEF1CsGAPDH1CsCLATHRIN1CsSAND1Cs18S rRNA1CsTIP1CsPTB1CsUBQ1CsUBC1CsACTIN1CsACTIN1
Different Organs Ranking Order
∆CTCsCLATHRIN1CsEF1CsGAPDH1CsPTB1CsSAND1CsTIP1CsUBC1CsUBQ1Cs18S rRNA1CsACTIN1CsACTIN1
BestKeeperCs18S rRNA1CsPTB1CsSAND1CsTIP1CsCLATHRIN1CsEF1CsUBC1CsGAPDH1CsUBQ1CsACTIN1CsACTIN1
NormfinderCsCLATHRIN1CsUBC1CsGAPDH1CsEF1CsTIP1CsPTB1CsSAND1CsUBQ1Cs18S rRNA1CsACTIN1CsACTIN1
GenormCsPTB1|CsSAND1CsEF1CsGAPDH1CsCLATHRIN1CsUBQ1CsTIP1CsUBC1Cs18S rRNA1CsACTIN1CsACTIN1
Recommended comprehensive rankingCsCLATHRIN1CsPTB1CsSAND1CsEF1CsGAPDH1Cs18S rRNA1CsUBC1CsTIP1CsUBQ1CsACTIN1CsACTIN1
Leaves with Cold and Short Photoperiod Treatment Ranking Order
∆CTCsCLATHRIN1CsGAPDH1CsPTB1CsUBQ1CsSAND1CsUBC1CsACTIN1CsEF1Cs18S rRNA1CsTIP1CsACTIN1
BestKeeperCsEF1CsSAND1Cs18S rRNA1CsGAPDH1CsUBC1CsPTB1CsCLATHRIN1CsACTIN1CsTIP1CsUBQ1CsACTIN1
NormfinderCsCLATHRIN1CsGAPDH1CsUBQ1CsPTB1CsSAND1CsUBC1CsACTIN1CsEF1Cs18S rRNA1CsTIP1CsACTIN1
GenormCsCLATHRIN1|CsUBQ1CsGAPDH1CsPTB1CsUBC1CsSAND1CsEF1CsACTIN1Cs18S rRNA1CsTIP1CsACTIN1
Recommended comprehensive rankingCsCLATHRIN1CsGAPDH1CsUBQ1CsPTB1CsSAND1CsEF1CsUBC1Cs18S rRNA1CsACTIN1CsTIP1CsACTIN1
Shoots with Cold and Short Photoperiod Treatment Ranking Order
∆CTCsUBC1CsGAPDH1CsPTB1CsCLATHRIN1CsUBQ1CsEF1CsSAND1Cs18S rRNA1CsTIP1CsACTIN1CsACTIN1
BestKeeperCsPTB1CsGAPDH1CsUBC1CsCLATHRIN1CsEF1Cs18S rRNA1CsACTIN1CsSAND1CsTIP1CsUBQ1CsACTIN1
NormfinderCsUBC1CsGAPDH1CsPTB1CsCLATHRIN1CsUBQ1CsEF1Cs18S rRNA1CsSAND1CsTIP1CsACTIN1CsACTIN1
GenormCsUBC1|CsCLATHRIN1CsPTB1CsTIP1CsSAND1CsGAPDH1CsUBQ1CsEF1Cs18S rRNA1CsACTIN1CsACTIN1
Recommended comprehensive rankingCsUBC1CsPTB1CsGAPDH1CsCLATHRIN1CsEF1CsUBQ1CsSAND1CsTIP1Cs18S rRNA1CsACTIN1CsACTIN1
Auxin Antagonist Auxinole Treatment Ranking Order
∆CTCsPTB1CsCLATHRIN1CsUBC1CsSAND1CsEF1CsACTIN1CsUBQ1CsTIP1CsGAPDH1CsACTIN1Cs18S rRNA1
BestKeeperCsPTB1CsUBC1CsCLATHRIN1CsACTIN1CsEF1Cs18S rRNA1CsSAND1CsACTIN1CsGAPDH1CsTIP1CsUBQ1
NormfinderCsPTB1CsCLATHRIN1CsUBC1CsSAND1CsEF1CsACTIN1CsUBQ1CsTIP1CsGAPDH1CsACTIN1Cs18S rRNA1
GenormCsPTB1|CsCLATHRIN1CsUBC1CsACTIN1CsSAND1CsTIP1CsEF1CsUBQ1CsGAPDH1CsACTIN1Cs18S rRNA1
Recommended comprehensive rankingCsPTB1CsCLATHRIN1CsUBC1CsSAND1CsACTIN1CsEF1CsTIP1CsUBQ1CsGAPDH1Cs18S rRNA1CsACTIN1
Lanolin Treatment Ranking Order
∆CTCsPTB1CsUBC1CsCLATHRIN1CsEF1CsSAND1CsTIP1CsACTIN1CsGAPDH1Cs18S rRNA1CsUBQ1CsACTIN1
BestKeeperCsEF1CsUBC1CsPTB1CsCLATHRIN1CsSAND1CsTIP1CsACTIN1Cs18S rRNA1CsGAPDH1CsUBQ1CsACTIN1
NormfinderCsPTB1CsUBC1CsCLATHRIN1CsSAND1CsEF1CsTIP1CsACTIN1CsGAPDH1Cs18S rRNA1CsUBQ1CsACTIN1
GenormCsUBC1|CsCLATHRIN1CsEF1CsPTB1CsTIP1CsSAND1CsACTIN1CsGAPDH1Cs18S rRNA1CsUBQ1CsACTIN1
Recommended comprehensive rankingCsUBC1CsPTB1CsCLATHRIN1CsEF1CsSAND1CsTIP1CsACTIN1CsGAPDH1Cs18S rRNA1CsUBQ1CsACTIN1
All Samples Ranking Order
∆CTCsEF1CsCLATHRIN1CsGAPDH1CsPTB1CsUBC1CsSAND1CsTIP1CsUBQ1Cs18S rRNA1CsACTIN1CsACTIN1
BestKeeperCsPTB1CsSAND1CsUBC1CsTIP1CsEF1CsGAPDH1CsCLATHRIN1Cs18S rRNA1CsUBQ1CsACTIN1CsACTIN1
NormfinderCsCLATHRIN1CsEF1CsGAPDH1CsUBC1CsPTB1CsUBQ1CsTIP1CsSAND1Cs18S rRNA1CsACTIN1CsACTIN1
GenormCsPTB1|CsSAND1CsUBC1CsTIP1CsEF1CsGAPDH1CsCLATHRIN1CsUBQ1Cs18S rRNA1CsACTIN1CsACTIN1
Recommended comprehensive rankingCsPTB1CsEF1CsSAND1CsCLATHRIN1CsUBC1CsGAPDH1CsTIP1CsUBQ1Cs18S rRNA1CsACTIN1CsACTIN1
The stability of candidate reference genes was validated under different experimental series, including diurnal expression in leaves, different organs, leaves and shoots with cold and short photoperiod treatment, auxinole and lanolin treatment. The experimental conditions were chosen based on the current focuses in physiological and molecular studies of tea plant. The expressions of two target genes CsSVP1 and CsFT1 were examined as negative controls since they were expected to show variation in expression in most of the sample treatments tested. In Arabidopsis, SVP plays an important role in the response of plants to ambient temperature changes [29]. FT, as a member of the small phosphatidylethanolamine-binding protein family, is known to have a tissue-specific gene expression and can be affected by light, photoperiod, temperature and hormone [30]. As expected, the transcript accumulation profile of CsSVP1 was significantly suppressed in shoots after cold and short photoperiod treatment. Likewise, CsFT1 gene had significant and expected expression fluctuation across the entire experimental series (Figure 2). The observed expected expression patterns of the target genes suggests that the relative quantitation of the cDNAs were accurate and that there were no obvious differences resulting from altered amplification efficiencies between samples in any given experiment. Thus, the stability rankings of 11 candidate reference genes is also likely reliable. Due to the differences in statistical algorithms, the four computational programs did not place the order of top ranked genes equally. In each experimental series, a recommend comprehensive ranking was given based on the ranking results from the four analysis programs. Unsurprisingly, no single gene showed stable transcript accumulation across all treatments tested. Consequently, this suggests that reference gene choice should be tested for any given experimental design to insure accurate quantitation of target RNAs. CsUBC1 was the most stable reference gene in both shoots with cold and short photoperiod treatment and lanolin treatment series but performed poorly in diurnal expression series. Above results further indicated that it is difficult to identify a single universally applicable reference gene covering all conditions.
Although no single gene showed consistent expression transcript accumulation over all treatments, several reference genes performed better than the others, suggesting that these genes might be a good choice for reference genes under different experimental settings in tea plants. According to the resulting overall ranking, CsPTB1, CsEF1, CsSAND1, CsCLATHRIN1 and CsUBC1 were the top five ideal reference genes (Table 2). CsPTB1 was chosen as a candidate reference gene because it was stably expressed in Arabidopsis across multiple Affymetrix ATH1 whole-genome GeneChip studies [27]. This gene was also subsequently identified as a stably expressed gene in leafy spurge suggesting it is a good first choice in other plant species [8]. CsPTB1 gene encodes polyprimidine tract-binding protein which is involved in the regulation of pre-mRNA alternative splicing, internal ribosomal entry site-mediated translation and mRNA localizing/sorting [31]. Though the CsPTB1 gene only ranked first at auxinole treatment series, it was ranked among the four top genes during all the experimental series except diurnal expression. Therefore, we suggest CsPTB1 to be the best reference gene for tea plant in experiments not involving diurnal expression changes. CsSAND1, CsCLATHRIN1, and CsEF1 genes were the next in ranking for top reference genes. CsSAND1 was first identified in animal studies and had been validated as a good reference gene in Arabidopsis, leafy spurge, citrus, and grape [8,32,33,34,35,36]. Clathrins are important scaffolding proteins in plant endocytosis and also showed stable character in other reference gene selection studies [37,38,39,40]. CsEF1, which encodes a protein that plays an important role in translation and cytoskeletal rearrangements [41], was perhaps the most superior traditional reference gene [27], and similar genes demonstrated stable expression in different tissue of Brassica rapa [42], during developmental process in Populus [43], and under the effects of different biotic factors (development stage, tissue, and strain) in Plutella xylostella [44]. In our study, CsEF1 gene was the most stable reference gene in diurnal expression series and also had good ranking in other experimental serials.
CsACTIN1, CsTUBULIN1 and Cs18S rRNA1 performed poorly among the 11 selected candidate genes in the experimental series of this study. This is likely due to the fact that these genes are members of larger gene families, and it has been demonstrated that individual members of these genes families have tissue specific and/or differential expression in response to environmental cues [45,46,47]. Thus, our results regarding the unstable expression of individual members of these gene families is not surprising and likely not unique to studies in tea plant. Indeed, the poor stability of CsACTIN1, CsTUBULIN1 and Cs18S rRNA1, which were considered as ideal internal controls in qRT-PCR previously, have been repeatedly reported recently [25,26,48,49]. Wong and Medrano [6] pointed out that the stability of these genes has to be validated under one’s own experimental condition before using them as reference genes. Several previous studies have used CsACTIN1 and Cs18S rRNA1 genes in tea plant as controls to examine various expression patterns of other genes [19,24]. Even compared with the target gene CsSVP1, the expression variation of CsACTIN1 and CsTUBULIN1 was greater in most of the experimental treatments. Therefore, these genes are not reliable reference genes for gene expression research in tea or most other systems. Our observation suggests that previous studies which used these common tea plant reference genes may need to re-examine their work in light of our findings, and our results should provide a guideline for future reference gene choice for gene expression studies in tea plant.

4. Experimental Section

4.1. Plant Materials and Treatments

Three-year-old potted small leaf tea plant (Camellia sinensis (L.) O. Kuntze) was grown in a greenhouse under a 16 h photoperiod at 23–25 °C for one year. Regular fertilization and pest management were applied. The fourth to sixth intact mature leaves away from the terminal growing point were sampled at different time of a day (3:00 am, 6:00 am, 9:00 am, 12:00 pm, 3:00 pm, 6:00 pm 9:00 pm and 12:00 am) for validation of reference genes during diurnal cycles. Different organs, including mature leaves (the fourth to sixth intact leaves away from the terminal growing point), whole flowers at half boom, shoots (including a bud and a leaf), five-month-old stems (generated in spring), axillary buds (taken from one-year-old branches) and roots, were collected separately at 9:00 pm. For analysis of transcript accumulation during cold stress, greenhouse-grown tea plants were moved to a growth chamber set for a 10 h photoperiod at 10 °C. Both mature leaves (the fourth to sixth intact leaves away from the terminal growing point) and shoots (including a bud and a leaf) were collected after 0 day, 1 day, 2 day, 3 day, 4 day, 1 week, 2 week, 3 week and 4 week cold treatment. In addition, an auxin antagonist (auxinole, kindly provided by Ken-ichiro Hayashi, Okayama University of Science, Okayama, Japan [50,51]) was applied to the growing apical buds of greenhouse-grown tea plants as a 2 mM solution in lanolin. Auxinole preparation and treatment were followed the description by Iino [52]. Same number of apical buds was treated with lanolin as control treatment. Shoots were harvested after 2 h, 1 day, 2 day, and 8 day auxinole and lanolin treatment. Two to three biological replicates were collected from each treatment. All the samples were immediately frozen in liquid nitrogen and stored at −80 °C until extraction of RNA.

4.2. Total RNA Extraction and cDNA Template Preparation

Total RNA extractions were performed following the pine tree extraction method [53]. The extracted RNAs were quantified using a NanoDrop® ND-100 Spectrophotometer (Thermo Fisher Scientific, Waltham, MS, USA) and the quantity and quality of RNAs was verified by 1% gel electrophoresis. For each sample, 5 µg RNA was used for cDNA synthesis following DNase treatment which was carried out using DNase I kit (Invitrogen, Carlsbad, CA, USA). cDNA synthesis was carried out using SuperScript® III First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA). 0.4 µCi of [α-32P]-dCTP (PerkinElmer, Waltham, MS, USA) was added with dNTP mix during cDNA synthesis. To filter out free [α-32P]-dCTP, the final cDNA product was purified using a QIAquick® PCR Purification Kit (QIAGEN, Valencia, CA, USA) and eluted with 50 µL final dilution of water. Quantitation of the resulting cDNA was determined by liquid scintillation as follows: 2 µL of eluted cDNA was mixed with 10 mL of scintillation cocktail (MP BIOMEDICALS, Irvine, CA, USA) and incorporated radioactivity was detected on Beckman LS 6000 Scintillation Counter (GMI, Albertville, MN, USA). The quantity of all the cDNA products was normalized based on the incorporated radioactivity of the least concentrated sample and then all samples were diluted 20-fold with ultrapure water. The final cDNA templates were stored in −80 °C for qRT-PCR analysis.

4.3. Primer Design and Quantitative Real-Time PCR (qRT-PCR)

All candidate reference genes were identified by BLASTX of the Arabidopsis reference genes against a tea plant EST database developed from assembly of the combined data from the small read archive (SRA) for tea plant at NCBI (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/) (accession numbers: SRX020193, SRX147828, SRX153144, SRX202319, SRX149415). The sequences with high sequence similarity to each candidate reference genes were aligned separately to highlight potential areas of polymorphism. The primer pairs were designed by PrimerSelect of Lasergene 8 (DNASTAR, Inc., Madison, WI, USA) following the recommendations by Udvardi et al. [54], and primer pairs were selected from sequence region with the fewest polymorphisms. The primer pair information is provided in Table 1. The amplification specificity of designed primers was confirmed via standard qRT-PCR procedures under different annealing temperature (50, 53, 56 and 59 °C). PCR amplification efficiency was calculated using a two-fold serial dilution of the cocktail of cDNA templates following the description by Chao et al. [8]. All PCR products were examined by melting curve analysis and agarose gel electrophoresis.
Gene expression was examined by qRT-PCR on a LightCycle® 480 Real-Time PCR System (Roche, Indianapolis, IN, USA). For real-time PCR reaction, 4 µL cDNA template was added to a 10 µL PCR reaction mixture containing 5 µL of 2× LightCycle® 480 SYBR Green I Master (Roche, Indianapolis, IN, USA) and 0.5 µL of each primer (20 pmol). Thermal cycling was performed with the program setting of 10 min pre-incubation at 95 °C, 45 cycles amplification of 20 s at 95 °C, 10 s at appropriate annealing temperatures (56 °C for CsACTIN1, Cs18S rRNA1, CsEF1, CsGAPDH1, CsTUBULIN1, CsUBC1 and CsPTB1 and 59 °C for CsUBQ1, CsCLATHRIN1, CsSAND1 and CsTIP1), and 35 s at 70 °C, and followed by dissociation analysis with a temperature ramp step added to the end of the amplification with an initial temperature of 65 °C and a final temperature of 97 °C. Three technical replicates were carried out for each reaction.
To validate the normalization effect of cDNA templates and the stability of candidate reference genes, the expressions of two target genes CsSVP1 (Accession number: KM057787, forward primer: 5'-ACATTACGGCGAGGCAAGTGAC-3', reverse primer: 5'-CGGTGGCGGAGAAGATGATG-3') and CsFT1 (Accession number: KC149523, forward primer: 5'-CCTTTACAAGGTCTATCTCTCTCAGG-3', 5'-TATATGATGATGATGATGTGA-3'), were included as negative controls since it was expected that these two genes would likely have differential transcript accumulation in many of the experimental treatments we were testing. We tested the expressions of these target genes under the same qRT-PCR conditions except with 54 and 60 °C annealing temperature, respectively. The average mean and the fluctuation range of the cycle threshold (CT) values between candidate reference genes and target genes were compared. Heat-maps of gene expression profile were performed following the description by Chao et al. [8].
The stabilities of candidate reference genes were determined based on a web-based analysis tool RefFinder (http://www.leonxie.com/referencegene.php). RefFinder is a comprehensive tool which consists of four computational programs geNorm, Normfinder, BestKeeper, and the comparative ∆CT method. The expression stability of each tested candidate reference gene was compared and ranked by each program on the basis of their cycle threshold (CT) values. Furthermore, based on the rankings from each program, an individual gene was assigned an appropriate weight and finally given an overall final ranking by calculating the geometric mean of their weights.

5. Conclusions

We evaluated 11 candidate reference genes under different experimental conditions using reliable qRT-PCR assay based on accurate cDNA normalization. Using the web-based analytical tool RefFinder, CsPTB1, CsEF1, CsSAND1, CsCLATHRIN1 and CsUBC1 appeared to be the top five most stable reference genes for all the sample pools. However, variability in even these top performing reference genes suggests that reference gene choice needs to be evaluated for each experimental condition. In addition, CsUBQ1, Cs18S rRNA1, CsACTIN1 and CsTUBULIN1 seemed to be unsuitable as internal control among most experimental conditions we tested. Our study highlights problems with the currently used reference gene in tea plant for expression analysis under different conditions. Our results should prove useful for future works on expression analysis in tea plant using qRT-PCR.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/15/12/22155/s1.

Acknowledgments

This work was supported by Natural Science Foundation of Zhejiang Province (Z3100473), National Natural Science Foundation of China (31170650), Earmarked Fund for China Agriculture Research System (CARS-23) and Major Project for New Agricultural Varieties Breeding of Zhejiang Province (2012C2905-3). But these funding bodies have no roles in design, analysis and manuscript preparation. Xinyuan Hao acknowledges the fellowship from the China Scholarship Council (CSC). The authors wish to thank Laura Kelley and Brant Bigger for their technical assistance during this study.

Author Contributions

Xinyuan Hao, David P. Horvath, Wun S. Chao, Yajun Yang, Xinchao Wang and Bin Xiao conceived and designed the experiments. Xinyuan Hao performed the experiments and analyzed the data. Xinyuan Hao wrote the paper. Xinyuan Hao, David P. Horvath, Wun S. Chao, Yajun Yang, Xinchao Wang and Bin Xiao revised and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bustin, S.A.; Benes, V.; Nolan, T.; Pfaffl, M.W. Quantitative real-time RT-PCR—A perspective. J. Mol. Endocrinol. 2005, 34, 597–601. [Google Scholar] [CrossRef] [PubMed]
  2. Ginzinger, D.G. Gene quantification using real-time quantitative PCR: An emerging technology hits the mainstream. Exp. Hematol. 2002, 30, 503–512. [Google Scholar] [CrossRef] [PubMed]
  3. Klein, D. Quantification using real-time PCR technology: Applications and limitations. Trends Mol. Med. 2002, 8, 257–260. [Google Scholar] [CrossRef] [PubMed]
  4. Pfaffl, M.W. A-Z of Quantitative PCR; International University Line: La Jolla, CA, USA, 2004; pp. 87–112. [Google Scholar]
  5. Nolan, T.; Hands, R.E.; Bustin, S.A. Quantification of mRNA using real-time RT-PCR. Nat. Protoc. 2006, 1, 1559–1582. [Google Scholar] [CrossRef] [PubMed]
  6. Wong, M.L.; Medrano, J.F. Real-time PCR for mRNA quantitation. Biotechniques 2005, 39, 75–85. [Google Scholar] [CrossRef] [PubMed]
  7. Huggett, J.; Dheda, K.; Bustin, S.; Zumla, A. Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6, 279–284. [Google Scholar] [CrossRef] [PubMed]
  8. Chao, W.S.; Doğramaci, M.; Foley, M.E.; Horvath, D.P.; Anderson, J.V. Selection and validation of endogenous reference genes for qRT-PCR analysis in leafy spurge (Euphorbia esula). PLoS One 2012, 7, e42839. [Google Scholar] [CrossRef] [PubMed]
  9. Thellin, O.; Zorzi, W.; Lakaye, B.; de Borman, B.; Coumans, B.; Hennen, G.; Grisar, T.; Igout, A.; Heinen, E. Housekeeping genes as internal standards: Use and limits. J. Biotechnol. 1999, 75, 291–295. [Google Scholar] [CrossRef] [PubMed]
  10. Graham, H.N. Green tea composition, consumption, and polyphenol chemistry. Prev. Med. 1992, 21, 334–350. [Google Scholar] [CrossRef] [PubMed]
  11. Yang, C.S.; Landau, J.M. Effects of tea consumption on nutrition and health. J. Nutr. 2000, 130, 2409–2412. [Google Scholar] [PubMed]
  12. Fujiki, H.; Suganuma, M.; Imai, K.; Nakachi, K. Green tea: Cancer preventive beverage and/or drug. Cancer Lett. 2002, 188, 9–13. [Google Scholar] [CrossRef] [PubMed]
  13. Basu, A.; Sanchez, K.; Leyva, M.J.; Wu, M.; Betts, N.M.; Aston, C.E.; Lyons, T.J. Green tea supplementation affects body weight, lipids, and lipid peroxidation in obese subjects with metabolic syndrome. J. Am. Coll. Nutr. 2010, 29, 31–40. [Google Scholar] [CrossRef] [PubMed]
  14. Mondal, T.K.; Bhattacharya, A.; Laxmikumaran, M.; Ahuja, P.S. Recent advances of tea (Camellia Sinensis) biotechnology. Plant Cell Tissue Organ Cult. 2004, 76, 195–254. [Google Scholar] [CrossRef]
  15. Singh, K.; Raizada, J.; Bhardwaj, P.; Ghawana, S.; Rani, A.; Singh, H.; Kaul, K.; Kumar, S. 26s rRNA-based internal control gene primer pair for reverse transcription-polymerase chain reaction-based quantitative expression studies in diverse plant species. Anal. Biochem. 2004, 335, 330–333. [Google Scholar] [CrossRef] [PubMed]
  16. Gohain, B.; Bandyopadhyay, T.; Borchetia, S.; Bharalee, R.; Gupta, S.; Bhorali, P.; Agarwala, N.; Das, S. Identification and validation of stable reference genes in camellia species. J. Biotechnol. Pharm. Res. 2011, 2, 9–18. [Google Scholar]
  17. Wang, Y.; Jiang, C.-J.; Li, Y.-Y.; Wei, C.-L.; Deng, W.-W. CsICE1 and csCBF1: Two transcription factors involved in cold responses in Camellia sinensis. Plant Cell. Rep. 2012, 31, 27–34. [Google Scholar] [CrossRef] [PubMed]
  18. Singh, K.; Kumar, S.; Ahuja, P.S. Differential expression of histone H3 gene in tea (Camellia sinensis (L.) O. Kuntze) suggests its role in growing tissue. Mol. Biol. Rep. 2009, 36, 537–542. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, X.-C.; Zhao, Q.-Y.; Ma, C.-L.; Zhang, Z.-H.; Cao, H.-L.; Kong, Y.-M.; Yue, C.; Hao, X.-Y.; Chen, L.; Ma, J.-Q. Global transcriptome profiles of Camellia sinensis during cold acclimation. BMC Genomics 2013, 14. [Google Scholar] [CrossRef]
  20. Das, A.; Das, S.; Mondal, T.K. Identification of differentially expressed gene profiles in young roots of tea [Camellia sinensis (L.) O. Kuntze] subjected to drought stress using suppression subtractive hybridization. Plant Mol. Biol. Rep. 2012, 30, 1088–1101. [Google Scholar] [CrossRef]
  21. Wang, X.; Hao, X.; Ma, C.; Cao, H.; Yue, C.; Wang, L.; Zeng, J.; Yang, Y. Identification of differential gene expression profiles between winter dormant and sprouting axillary buds in tea plant (Camellia sinensis) by suppression subtractive hybridization. Tree Genet. Genomes 2014, 10, 1149–1159. [Google Scholar] [CrossRef]
  22. Wang, L.; Li, X.; Zhao, Q.; Jing, S.; Chen, S.; Yuan, H. Identification of genes induced in response to low-temperature treatment in tea leaves. Plant Mol. Biol. Rep. 2009, 27, 257–265. [Google Scholar] [CrossRef]
  23. Paul, A.; Kumar, S. Responses to winter dormancy, temperature, and plant hormones share gene networks. Funct. Integr. Genomics 2011, 11, 659–664. [Google Scholar] [CrossRef] [PubMed]
  24. Shi, C.-Y.; Yang, H.; Wei, C.-L.; Yu, O.; Zhang, Z.-Z.; Jiang, C.-J.; Sun, J.; Li, Y.-Y.; Chen, Q.; Xia, T. Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics 2011, 12. [Google Scholar] [CrossRef] [PubMed]
  25. Bustin, S.A. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol. 2000, 25, 169–193. [Google Scholar] [CrossRef] [PubMed]
  26. Dheda, K.; Huggett, J.F.; Bustin, S.A.; Johnson, M.A.; Rook, G.; Zumla, A. Validation of housekeeping genes for normalizing rna expression in real-time PCR. Biotechniques 2004, 37, 112–119. [Google Scholar] [PubMed]
  27. Czechowski, T.; Stitt, M.; Altmann, T.; Udvardi, M.K.; Scheible, W.-R. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005, 139, 5–17. [Google Scholar] [CrossRef] [PubMed]
  28. Ailenberg, M.; Silverman, M. Semiquantitative RT-PCR: Enhancement of assay accuracy and reproducibility. Biotechniques 1997, 22, 630–636. [Google Scholar] [PubMed]
  29. Lee, J.H.; Yoo, S.J.; Park, S.H.; Hwang, I.; Lee, J.S.; Ahn, J.H. Role of SVP in the control of flowering time by ambient temperature in Arabidopsis. Genes Dev. 2007, 21, 397–402. [Google Scholar] [CrossRef] [PubMed]
  30. Pin, P.; Nilsson, O. The multifaceted roles of FLOWERING LOCUS T in plant development. Plant Cell Environ. 2012, 35, 1742–1755. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, S.; Okamoto, T. Involvement of polypyrimidine tract-binding protein (ptb)-related proteins in pollen germination in Arabidopsis. Plant Cell Phys. 2009, 50, 179–190. [Google Scholar] [CrossRef]
  32. Cottage, A.; Edwards, Y.J.; Elgar, G. SAND, a new protein family: From nucleic acid to protein structure and function prediction. Comp. Funct. Genomics 2001, 2, 226–235. [Google Scholar] [CrossRef] [PubMed]
  33. Poteryaev, D.; Spang, A. A role of SAND-family proteins in endocytosis. Biochem. Soc. Trans. 2005, 33, 606–608. [Google Scholar] [CrossRef] [PubMed]
  34. Remans, T.; Smeets, K.; Opdenakker, K.; Mathijsen, D.; Vangronsveld, J.; Cuypers, A. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 2008, 227, 1343–1349. [Google Scholar] [CrossRef] [PubMed]
  35. Mafra, V.; Kubo, K.S.; Alves-Ferreira, M.; Ribeiro-Alves, M.; Stuart, R.M.; Boava, L.P.; Rodrigues, C.M.; Machado, M.A. Reference genes for accurate transcript normalization in citrus genotypes under different experimental conditions. PLoS One 2012, 7, e31263. [Google Scholar] [CrossRef] [PubMed]
  36. Reid, K.E.; Olsson, N.; Schlosser, J.; Peng, F.; Lund, S.T. An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6. [Google Scholar] [CrossRef] [PubMed]
  37. Migocka, M.; Papierniak, A. Identification of suitable reference genes for studying gene expression in cucumber plants subjected to abiotic stress and growth regulators. Mol. Breed. 2011, 28, 343–357. [Google Scholar] [CrossRef]
  38. Artico, S.; Nardeli, S.M.; Brilhante, O.; Grossi-de-Sa, M.F.; Alves-Ferreira, M. Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data. BMC Plant Biol. 2010, 10. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, L.; Zhong, H.-Y.; Kuang, J.-F.; Li, J.-G.; Lu, W.-J.; Chen, J.-Y. Validation of reference genes for RT-qPCR studies of gene expression in banana fruit under different experimental conditions. Planta 2011, 234, 377–390. [Google Scholar] [CrossRef] [PubMed]
  40. Demidenko, N.V.; Logacheva, M.D.; Penin, A.A. Selection and validation of reference genes for quantitative real-time PCR in buckwheat (Fagopyrum esculentum) based on transcriptome sequence data. PLoS One 2011, 6, e19434. [Google Scholar] [CrossRef] [PubMed]
  41. Shiina, N.; Gotoh, Y.; Kubomura, N.; Iwamatsu, A.; Nishida, E. Microtubule severing by elongation factor 1 alpha. Science 1994, 266, 282–285. [Google Scholar] [CrossRef] [PubMed]
  42. Qi, J.; Yu, S.; Zhang, F.; Shen, X.; Zhao, X.; Yu, Y.; Zhang, D. Reference gene selection for real-time quantitative polymerase chain reaction of mRNA transcript levels in chinese cabbage (Brassica rapa L. Ssp. Pekinensis.). Plant Mol. Biol. Rep. 2010, 28, 597–604. [Google Scholar] [CrossRef]
  43. Xu, M.; Zhang, B.; Su, X.; Zhang, S.; Huang, M. Reference gene selection for quantitative real-time polymerase chain reaction in Populus. Anal. Biochem. 2011, 408, 337–339. [Google Scholar] [CrossRef] [PubMed]
  44. Fu, W.; Xie, W.; Zhang, Z.; Wang, S.; Wu, Q.; Liu, Y.; Zhou, X.; Zhou, X.; Zhang, Y. Exploring valid reference genes for quantitative real-time PCR analysis in Plutella xylostella (lepidoptera: Plutellidae). Int. J. Biol. Sci. 2013, 9, 792–802. [Google Scholar] [CrossRef] [PubMed]
  45. An, Y.Q.; McDowell, J.M.; Huang, S.; McKinney, E.C.; Chambliss, S.; Meagher, R.B. Strong, constitutive expression of the Arabidopsis ACT2/ACT8 actin subclass in vegetative tissues. Plant J. 1996, 10, 107–121. [Google Scholar] [CrossRef] [PubMed]
  46. Bustin, S. Quantification of mrna using real-time reverse transcription PCR (RT-PCR): Trends and problems. J. Mol. Endocrinol. 2002, 29, 23–39. [Google Scholar] [CrossRef] [PubMed]
  47. Ludwig, S.R.; Oppenheimer, D.G.; Silflow, C.D.; Snustad, D.P. Characterization of the α-tubulin gene family of Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 1987, 84, 5833–5837. [Google Scholar] [CrossRef] [PubMed]
  48. Tong, Z.; Gao, Z.; Wang, F.; Zhou, J.; Zhang, Z. Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol. Biol. 2009, 10. [Google Scholar] [CrossRef] [PubMed]
  49. Nicot, N.; Hausman, J.-F.O.; Hoffmann, L.; Evers, D. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J. Exper. Bot. 2005, 56, 2907–2914. [Google Scholar] [CrossRef]
  50. Hayashi, K.-I.; Neve, J.; Hirose, M.; Kuboki, A.; Shimada, Y.; Kepinski, S.; Nozaki, H. Rational design of an auxin antagonist of the SCFTIR1 auxin receptor complex. ACS Chem. Biol. 2012, 7, 590–598. [Google Scholar] [CrossRef] [PubMed]
  51. Hayashi, K.-I.; Tan, X.; Zheng, N.; Hatate, T.; Kimura, Y.; Kepinski, S.; Nozaki, H. Small-molecule agonists and antagonists of F-box protein-substrate interactions in auxin perception and signaling. Proc. Natl. Acad. Sci. USA 2008, 105, 5632–5637. [Google Scholar] [CrossRef] [PubMed]
  52. Iino, M. Gravitropism and phototropism of maize coleoptiles: Evaluation of the Cholodny-Went theory through effects of auxin application and decapitation. Plant Cell Physiol. 1995, 36, 361–367. [Google Scholar]
  53. Chang, S.; Puryear, J.; Cairney, J. A simple and efficient method for isolating RNA from pine trees. Plant Mol. Biol. Rep. 1993, 11, 113–116. [Google Scholar] [CrossRef]
  54. Udvardi, M.K.; Czechowski, T.; Scheible, W.-R. Eleven golden rules of quantitative RT-PCR. Plant Cell 2008, 20, 1736–1737. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Hao, X.; Horvath, D.P.; Chao, W.S.; Yang, Y.; Wang, X.; Xiao, B. Identification and Evaluation of Reliable Reference Genes for Quantitative Real-Time PCR Analysis in Tea Plant (Camellia sinensis (L.) O. Kuntze). Int. J. Mol. Sci. 2014, 15, 22155-22172. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms151222155

AMA Style

Hao X, Horvath DP, Chao WS, Yang Y, Wang X, Xiao B. Identification and Evaluation of Reliable Reference Genes for Quantitative Real-Time PCR Analysis in Tea Plant (Camellia sinensis (L.) O. Kuntze). International Journal of Molecular Sciences. 2014; 15(12):22155-22172. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms151222155

Chicago/Turabian Style

Hao, Xinyuan, David P. Horvath, Wun S. Chao, Yajun Yang, Xinchao Wang, and Bin Xiao. 2014. "Identification and Evaluation of Reliable Reference Genes for Quantitative Real-Time PCR Analysis in Tea Plant (Camellia sinensis (L.) O. Kuntze)" International Journal of Molecular Sciences 15, no. 12: 22155-22172. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms151222155

Article Metrics

Back to TopTop