Validation of internal control for gene expression study in soybean by quantitative real-time PCR
© Jian et al; licensee BioMed Central Ltd. 2008
Received: 06 November 2007
Accepted: 23 June 2008
Published: 23 June 2008
Normalizing to housekeeping gene (HKG) can make results from quantitative real-time PCR (qRT-PCR) more reliable. Recent studies have shown that no single HKG is universal for all experiments. Thus, a suitable HKG should be selected before its use. Only a few studies on HKGs have been done in plants, and none in soybean, an economically important crop. Therefore, the present study was conducted to identify suitable HKG(s) for normalization of gene expression in soybean.
All ten HKGs displayed a wide range of Ct values in 21 sample pools, confirming that they were variably expressed. GeNorm was used to determine the expression stability of the HGKs in seven series sets. For all the sample pools analyzed, the stability rank was ELF1B, CYP2 > ACT11 > TUA > ELF1A > UBC2 > ACT2/7 > TUB > G6PD > UBQ10. For different tissues under the same developmental stage, the rank was ELF1B, CYP2 > ACT2/7 > UBC2 > TUA > ELF1A > ACT11 > TUB > G6PD > UBQ10. For the developmental stage series, the stability rank was ACT2/7, TUA > ELF1A > UBC2 > ELF1B > TUB > CYP2 > ACT11 > G6PD > UBQ10. For photoperiodic treatments, the rank was ACT11, ELF1B > CYP2 > TUA > ELF1A > UBC2 > ACT2/7 > TUB > G6PD > UBQ10. For different times of the day, the rank was ELF1A, TUA > ELF1B > G6PD > CYP2 > ACT11 > ACT2/7 > TUB > UBC2 > UBQ10. For different cultivars and leaves on different nodes of the main stem, the ten HKGs' stability did not differ significantly. ΔCt approach and 'Stability index' were also used to analyze the expression stability in all 21 sample pools. Results from ΔCt approach and geNorm indicated that ELF1B and CYP2 were the most stable HKGs, and UBQ10 and G6PD the most variable ones. Results from 'Stability index' analysis were different, with ACT11 and CYP2 being the most stable HKGs, and ELF1A and TUA the most variable ones.
Our data suggests that HKGs are expressed variably in soybean. Based on the results from geNorm and ΔCt analysis, ELF1B and CYP2 could be used as internal controls to normalize gene expression in soybean, while UBQ10 and G6PD should be avoided. To achieve accurate results, some conditions may require more than one HKG to be used for normalization.
Gene expression analysis is becoming much more prevalent since it promotes our understanding of biological processes. Compared with the traditional methods for transcript analysis including Northern blot, RNase protection analysis, in situ hybridization and semi-RT-PCR, the fluorescence-based qRT-PCR has recently been considered as the most reliable method for the detection of mRNA  because of its high sensitivity, no post-PCR processing , and wide dynamic range , which allows a straightforward comparison between RNAs that differ widely in their abundance. Furthermore, it is easy to use, allows high throughput production of data and eliminates the need for radioactive isotopes . Moreover, it is especially suitable when only a small number of cells are available. Although qRT-PCR is frequently used due to these advantages, some disadvantages may include variations between samples which may differ in the amount and quality of starting material, RNA preparation, cDNA synthesis, dilutions and pipetting. Normalizing a target gene to the HKGs makes qRT-PCR reliable by minimizing the variations.
The HKGs, which are referred to as internal controls or reference genes, are presumed to have constant expression level among different tissues and at all developmental stages, regardless of the experimental conditions or treatments. Additionally, the HKG and target gene should have similar transcript levels to avoid analytical problems . Commonly used HKGs are cellular maintenance genes, which regulate basic and ubiquitous cellular functions , such as components of the cytoskeleton (actins), glycolytic pathway (glyceraldehyde-3-phosphate dehydrogenase (GAPDH)), protein folding (cyclophilin), synthesis of ribosome subunits (rRNA), electron transporter (succinate dehydrogenase complex, SDH), protein degradation (ubiquitin), etc. These genes are supposed to have constant expression levels between different samples, and are frequently used as a normalizer without proper validation. However, recent studies show that the transcriptional levels of these HKGs are not always stable, and that no single HKG has a constant level under all experimental conditions [8–10]. A recent study even suggests that such a 'foolproof' gene does not exist . The reason for this expression variability may be that the HKGs not only take part in the basic cell metabolism but also participate in other cellular process [12, 13]. Therefore, selecting a suitable HKG(s) which has a constant expression level in certain experimental conditions for normalization is crucial for getting accurate results in gene expression studies.
Recently, many procedures have been constructed to find the best suitable HKG(s) in a set of samples, such as geNorm , NormFinder , ΔCt approach  and 'Stability index' . For example, using geNorm, YWHAZ, GAPD and SDHA were found to be the most stable HKGs across the examined embryonic stages in bovine pre-implantation embryos, while the commonly used ACTB was variably expressed . By comparing the expression results of the non-stimulated tissues and leucocytes from Atlantic salmon (Salmo salar L.) using the Normfinder program, it was shown that EF1-alpha was the most stably expressed gene . Using the ΔCt approach, GAPDH was found to be the most suitable HKG for expression studies in reticulocytes while the commonly used B2M should be avoided . UBQ and TUA were selected as reference genes to normalize gene expression in a single female poplar hybrid clone (P. trichocarpa × P. deltoidies) using the 'Stability index' .
Nevertheless, many studies on HKGs selection refer to human or animal tissue. As far as is known, only a few have been focused on plants such as rice [19–21], poplar , potato  and Arabidopsis thaliana . Moreover, there is no report on soybean (Glycine max [L.] Merr.), a very important crop and a model plant in the early studies of photoperiodism [24, 25].
Photoperiod controls several responses throughout the plant life cycle, such as germination , flowering induction , post-flowering development [27, 28], maturity, dormancy , and yield formation [30, 31]. The photoperiodic control of flowering in Arabidopsis, a long-day (LD) plant, is a hot topic attracting many scientists to enter this field and helping to better understand the processes involved [32–35]. However, there is little known about the mechanism in soybean, a typical short-day (SD) plant. Thus, the understanding of some key genes' expression patterns will help illuminate the mechanism involved in this process. Different cultivars of soybean may have different sensitivity to photoperiod; studies of the expression pattern of key genes in the photoperiodic pathway may also help elucidate what leads to the varietal difference.
In the present study, the expression profiles of ten HKGs, including ACT11, ACT2/7, TUA, TUB, UBC2, CYP2, G6PD, ELF1A, ELF1B and UBQ10, were studied during the development of Zigongdongdou (ZGGG), a late-maturing soybean (Glycine max [L.] Merr) cultivar, under LD and SD conditions. The expression patterns of the ten HKGs were also detected in soybean cultivar Heihe No. 27 (HH27), Zhonghuang No. 24 (ZH24) and Suinong No. 14 (SN14). GeNorm, ΔCt approach and 'Stability index' were used to assess the value of ten HKGs as suitable internal control(s) for soybean gene expression studies.
Expression profiling of HKGs
'Stability index' assay
Summary of statistics measuring stability of HKG expression
The ΔCt approach was employed by Silver et al. to select the most suitable HKG in reticulocytes. If the ΔCt value between two HKGs does not change when analyzed in all the samples, it means either both genes have stable expression patterns or they are co-regulated among those samples. However, the fluctuation in ΔCt means that at least one of them was variably expressed. Introduction of a third, fourth, or fifth gene into the comparisons will provide more information on which pairs show less variability and hence which gene(s) has stable expression among the samples tested. It is easy for large panels of genes to be compared with one another and then selected or eliminated on the basis of ΔCt. Ultimately, an appropriate HKG can be selected for a particular experimental system .
Relative expression of GmBFT at different developmental stages
Relative expression of GmBFT using different HKG(s) as internal controls (mean ± StdDev)
0.47 ± 0.35
1.08 ± 0.13
0.92 ± 0.18
2.59 ± 0.09
0.85 ± 0.10
1.18 ± 0.23
0.25 ± 0.12
0.95 ± 0.34
1.16 ± 0.13
2.56 ± 0.35
1.24 ± 0.25
0.97 ± 0.21
Relative expression of GmBFT under different photoperiodic treatments
Relative expression of GmBFT using different HKG(s) as internal controls (mean ± StdDev)
2.59 ± 0.45
0.98 ± 0.09
0.75 ± 0.12
0.87 ± 0.08
1.32 ± 0.15
1.00 ± 0.15
2.78 ± 0.34
1.15 ± 0.04
1.13 ± 0.08
0.11 ± 0.12
0.54 ± 0.12
0.63 ± 0.11
qRT-PCR has become a powerful tool for gene expression analysis because of its high throughput, sensitivity and accuracy [1, 4]. The use of suitable HKGs to normalize the variation made by RNA preparation, cDNA synthesis or PCR processing would make the results more reliable. In order to select suitable HKG(s) for normalization, many procedures such as geNorm , Normfinder , ΔCt approach  and 'Stability index'  have been used. Since all methods mentioned above are based on the Ct value, which is determined mostly by the quantity of cDNA , the prerequisite for selecting a set of reliable HKG(s) is based on the equal input cDNA when doing qRT-PCR. In early studies, the most commonly used method to measure the input cDNA was by a spectrophotometer such as NanoDrop ND-1000 [20, 40, 41]. Considering the importance of input cDNA, more than one method should be used to detect its quantity and quality to ensure the cDNA equality for each PCR run in order to get reliable results for HKG(s) expression stability analysis. In the present study, cDNA was verified by measurements on both the ND-1000 and SMA3000 spectrophotometers to ensure the equality of cDNA in the PCR reactions and reliability of the results when using Ct values for analysis.
When gene expression stability in soybean was analyzed by geNorm, the most stable genes in the seven series were different as shown in Figure 3. In all seven series analyzed, ELF1B was the most stable HKG. UBQ10 and G6PD were the most variable ones, so these genes should be avoided as internal controls when doing gene expression studies in soybean. Our findings were in accordance with the result that UBQ10 exhibited the least stable expression in different tissues or cell types at different developmental stages in rice . Similarly, in the development of grape berry, UBQ10 was not the recommended HKG for normalization . However, in an earlier study in Arabidopsis, UBQ10 showed highly stable expression . An ubiquitin tag is not only used to mark particular proteins for proteolytic elimination, but can also have non-proteolytic functions  which may lead to the variable expression of ubiquitin in different plants. G6PD was suggested to be an inappropriate internal control in qRT-PCR studies of estrogens effects in fish  and it also showed significant differences in expression between malignant and nonmalignant pairs (at least p < 0.04) of human bladder cancer . There are few published works using G6PD as an internal control. This may be because G6PD not only acts as a component of the glycolytic pathway but also participates in other processes as well. Thus, the expression profile of G6PD might change according to the corresponding experimental conditions. In the present study, TUA was found to be one of the most stable HKGs with the lowest M value in the developmental series in soybean. This result was consistent with an earlier study in poplar with TUA as one of the most stable HKGs . Other most commonly used HKGs, like TUB, displayed an unacceptably high variable expression pattern limiting its use as an internal control except in the different cultivar series where all the other HKGs showed relatively stable expression. Taken together, these results suggested that these HKGs were regulated differently in each plant species and may exhibit differential expression patterns. Therefore, a HKG with stable expression in one organism may be not suitable to normalize gene expression in another organism under a given set of conditions and thus needs to be validated before its use.
To further verify the suitability of HKGs selected in the present study, GmBFT expression levels were detected at different developmental stages and under different photoperiodic conditions in soybean (Figure 5). In unifoliate leaves, GmBFT expression after a 25-day SD treatment was significantly higher than after a 1-day SD treatment, and was also higher than after a 25-day LD treatment. The case was similar in the shoot tips of soybean using HKGs selected in corresponding conditions as internal controls. The fact that this result is in accordance with earlier work by Sun et al.  means that the HKGs identified in this study are suitable at various development stages and photoperiodic conditions. The relative expression of GmBFT was also analyzed using UBQ10 and G6PD as internal controls (Figure 5). No significant expression difference of GmBFT was observed in unifoliate leaves at different developmental stages or under different photoperiodic conditions. The result was similar in shoot tips. Obviously, UBQ10 and G6PD are not suitable HKGs to normalize gene expression in soybean under such conditions.
The ΔCt approach  and 'Stability index'  were also used to analyze gene expression stability in all 21 sample pools in soybean to compare the accuracy of these three methods. Results obtained from the ΔCt method were very similar to that from the geNorm analysis. ELF1B was the most stable HKG followed by CYP2, while UBQ10 was the most variable HKG followed by G6PD. When 'Stability index' was used to measure the stability of HKG(s), results changed, with ACT11 being the most stable HKG followed by CYP2, while TUA and ELF1A became the least stable ones in all the 21 sample pools in soybean. An explanation might be that the 'Stability index'did not take the PCR efficiency into account, which played an important role in the data analysis. Thus, in order to get a reliable result by only comparing Ct deviation for the individual genes in the tested tissue samples, the HKGs analyzed should have similar PCR efficiency.
Although gene expression stability analyzed by the ΔCt approach was similar to that by geNorm, it was still not the first choice to get accurate normalization especially for the research which can get enough samples. The ΔCt method can be used to detect the most stable HKG as geNorm provides; however, it could not provide the number of HKGs necessary for accurate normalization. Some studies may require more than one HKG to be included, and using the single most stable HKG for normalization might not get the most accurate result. Indeed, as shown in Figure 3A, when all 21 sample pools were analyzed by geNorm, the geometric average of four HKGs (ELF1B, CYP2, ACT11 and TUA) was recommended as the normalization factor to get accurate results in soybean. However, the analysis by ΔCt method only indicated that the most stable HKG was ELF1B, and it could not be used to find out the optimal number of HKGs. Thus, gene expression accuracy might be undermined if the result is only based on ELF1B. It is true that the ΔCt method is useful for validating the most stable HKG in some specific tissue samples or cell types for which it is difficult to obtain enough material, such as reticulocytes . However, for material where RNA samples are easy to obtain, such as plants, geNorm is recommended because it is easy to determine the optimal number of stable HKGs for accurate normalization.
Selection of suitable HKG(s) is necessary for accurate gene expression, but it is quite expensive and time-consuming. To avoid the additional expense and labor of using multiple internal control genes, a potential strategy suggested by Brunner et al.  was to design a PCR primer pair which could amplify two or more members of a control gene family. However, for two members of the same HKG family, their expression pattern and stability might vary. Therefore, results may be questionable for gene expression based on these two HKGs. For example, UBQ5 was one of the most suitable HKGs in a given set of tissue samples in rice, while UBQ10 expressed variably . A similar situation was observed for the actin gene family in the developmental stage series in the present study. ACT2/7 was stably expressed, while ACT11 showed variable profiling. Thus, using primer pair sets to amplify two or more members of a HKG family for qRT-PCR may not be recommendable. After all, the accuracy of the results should always be given first priority.
The photoperiod plays an important role throughout the life cycle of soybean [27, 28, 46]. The sensitivity of different cultivars of soybean to photoperiod is quite variable. Understanding the mechanisms involved in this process will be beneficial in the molecular breeding of soybean. In the present study, the most stable HKG in the different cultivar series under the same photoperiodic conditions and developmental stage was UBC2, whereas ACT2/7 was the most variable one. The difference of average expression stability among the ten HKGs in this series was smaller than that in the other series. Thus, the gene expression would be similar using any of the ten HKGs as internal controls when studying different cultivar. It would be helpful to understand the mechanisms involved in differential sensitivities to photoperiod for different cultivars.
A large number of studies have been carried out concerning the validation of HKG(s) in many different tissue samples and cell types. However, there is no correlative report in soybean, a SD dicot. Our data showed the variable expression profiles of ten commonly used HKGs in different tissue samples and under different photoperiodic conditions in soybean. Based on geNorm and ΔCt methods, ELF1B and CYP2 appears to be the most suitable HKGs to normalize gene expression during the development of soybean, while UBQ10 and G6PD seems to be unsuitable as reference genes. Under some conditions, more than one HKG should be used as internal controls to normalize gene expression in soybean in order to get the most reliable results.
Sample collection and RNA extraction
Soybean tissues used for gene expression analysis
Time after photoperiodic treatment for seeding
Total RNA was extracted using Trizol (Invitrogen) according to the Manufacture's Instruction with little modification. One more chloroform extraction step was added to the RNA extraction process. RNA was quantified by the absorbance at OD260 using NanDrop ND-1000 spectrophotometer. The absorbance ratio at OD260/280 and OD260/230 were used to assess the purity of all the RNA samples. Only RNA samples with OD260/280 ratio (protein contamination) between 1.8–2.0 and OD260/230 (organic pollutant) higher than 2.0 was used for the further analysis. RNA integrity was verified by 2% agar gel electrophoresis and ethidium bromide staining. The samples with 25S/18S ribosomal RNA between 1.5–2.0 and absence of smears were used for the following experiment.
cDNA synthesis and quantification
Before cDNA synthesis, 5 μg total RNA was treated with RQ1 RNase-free DNase (Promega) according to the Manufacture's Instruction to ensure no DNA contamination, and then cDNA synthesis was carried out with the purified RNA using the SuperScript III First-Strand Synthesis System (Invitrogen) following the instruction. The RT reaction was performed using Mastercycler Gradient (Eppendorf). Briefly, 1 μg RNA, 50 μM oligodT(20) and 10 mM dNTP mix were added together to incubate at 65°C for 5 min, then placed on ice for at least 1 min. After that, 2 μl 10 × RT buffer, 1 μl 25 mM MgCl2, 2 μl 0.1 M DTT, 40 U RNaseOUT and 200 U SuperScript III were added and then incubated at 50°C for 50 min. The RT reaction was terminated by incubating at 85°C for 5 min and the residual RNA was removed by incubated at 37°C for 20 min with the addition of 1 μl RNaseH. After cDNA was synthesized, it was used as the template for PCR amplification. This amplification was made using primer pair sets which span an intron to detect DNA contamination.
cDNA was 2× diluted before quantification and the quantity and quality of input cDNA were determined by the SMA3000 and NanDrop ND-1000 spectrophotometer to make sure the cDNA amount for each PCR run is equal. For each method, the measurement was done in duplicate. The slope of the regression line of the concentrations measured with both methods did not differ indicating the equality of cDNA measurement with both methods for the qRT-PCR reaction (Yconc.ND-1000 = 1.0444 × Xconc.SMA3000+0.5443; n = 21; r = 0.9763). To reduce the system error, all the cDNA was diluted to about 2.3 ng/μl, so there were 20 ng/8.8 μl cDNA for the real-time RT-PCR reaction. All the cDNA were stored at -20°C until PCR.
Selection of soybean sequences, primer design and PCR optimization
Description of soybean genes for qRT-PCR
Soybase Accession Number
Arabidopsis homolog locusb
Arabidopsis locus description
BLASTX score/E Value
Cytoskeletal structural protein
Cytoskeletal structural protein
Glucose metabolic process
Eukaryotic elongation factor 1-beta
Ubiquitin-conjugating enzyme E2
Ubiquitin-dependent protein catabolic process
Eukaryotic elongation factor 1-alpha
Structural constituent of cytoskeleton
Structural constituent of cytoskeleton
Protein binding, protein modification process
Brother of FT and TFL1 protein
Phosphatidylethanola -mine binding
Primers and amplicons for each of the 10 HKGs and GmBFT
Forward Primer Sequence [5'-3']
Reverse Primer Sequence [5'-3']
Amplicon Length (bp)
qRT-PCR was conducted on ABI PRISM 7000 Sequence Detection System using power SYBR Green Mix (Applied Biosystems, USA). Each reaction was run in a 20 μl volume which contained 8.8 μl cDNA equal to 20 ng, 10 μl 2 × power SYBR mix, 0.6 μl each primer to a final concentration of 300 nM. All the reactions were performed as the following conditions: 2 min at 50°C, 10 min at 95°C, and 40 cycles of 10 s at 95°C, and 1 min at 60°C in 96-well optical reaction plates (Applied Biosystems, USA). To verify the specificity of the amplicon for each primer pair, a melting curve was made from 60°C to 95°C at the end of each PCR run and all the ten primer pairs amplified a single product. The PCR efficiencies showed in Table 6 for each gene was determined with the slope of a liner regression model. Each cDNA sample pool was bulked and then used as the PCR template in a range of 50, 25, 10, 5, and 2 ng . The corresponding real-time PCR efficiencies were calculated according to the equation: E = 10-1/slope .
Expression levels of the ten HKGs in all the sample pools were determined by the number of cycles (Ct) needed for the amplification related fluorescence to reach a specific threshold level of detection . The raw Ct value obtained from ABI 7000 after each PCR run was converted into relative quantities using the PCR efficiencies for each gene according to the requirement of geNorm software [11, 49] to calculate gene expression stability (M). The expression stability of the ten HKGs was also determined by the 'Stability index'  and ΔCt approach  to compare the three methods in all the 21 sample pools.
quantitative real-time PCR
- ACT11 :
- ACT2/7 :
- G6PD :
- ELF1B :
eukaryotic elongation factor 1-beta
- UBC2 :
ubiquitin-conjugating enzyme E2
- ELF1A :
eukaryotic elongation factor 1-alpha
- TUB :
- TUA :
- CYP2 :
- UBQ10 :
Heihe No. 27
Zhonghuang No. 24
Suinong No. 14
- GmBFT :
Glycine max brother of FT and TFL1.
We thank Hongbo Sun and Shikui Song for suggestions in performing the experiment and Jiantian Leng and Ying Wang for assistance in data analysis. We especially thank Charles H. Leseberg and Tore Brembu for reading the manuscript and their suggestions for the writing process. This work was supported by the National Natural Science Foundation of China (30471054), National High Technology Research and Development Program of China (2007AA10Z133) and Beijing Natural Science Foundation (5042019).
- Bustin S: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of molecular endocrinology 2000, 25: 169-193. 10.1677/jme.0.0250169.View ArticlePubMed
- Ginzinger DG: Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Experimental hematology 2002, 30(6):503-512. 10.1016/S0301-472X(02)00806-8.View ArticlePubMed
- Garson JA, Grant PR, Ayliffe U, Ferns RB, Tedder RS: Real-time PCR quantitation of hepatitis B virus DNA using automated sample preparation and murine cytomegalovirus internal control. Journal of virological methods 2005, 126(1–2):207-213. 10.1016/j.jviromet.2005.03.001.View ArticlePubMed
- Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochemical and biophysical research communications 2004, 313(4):856-862. 10.1016/j.bbrc.2003.11.177.View ArticlePubMed
- Gal AB, Carnwath JW, Dinnyes A, Herrmann D, Niemann H, Wrenzycki C: Comparison of real-time polymerase chain reaction and end-point polymerase chain reaction for the analysis of gene expression in preimplantation embryos. Reproduction, fertility, and development 2006, 18(3):365-371. 10.1071/RD05012.View ArticlePubMed
- Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, Kristiansen G, Jung K: In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC molecular biology 2007, 8: 47. 10.1186/1471-2199-8-47.PubMed CentralView ArticlePubMed
- de Kok JB, Roelofs RW, Giesendorf BA, Pennings JL, Waas ET, Feuth T, Swinkels DW, Span PN: Normalization of gene expression measurements in tumor tissues: comparison of 13 endogenous control genes. Laboratory investigation; a journal of technical methods and pathology 2005, 85(1):154-159.View ArticlePubMed
- 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. Journal of biotechnology 1999, 75(2–3):291-295. 10.1016/S0168-1656(99)00163-7.View ArticlePubMed
- Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. Journal of biochemical and biophysical methods 2000, 46(1–2):69-81. 10.1016/S0165-022X(00)00129-9.View ArticlePubMed
- Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Analytical biochemistry 2002, 309(2):293-300. 10.1016/S0003-2697(02)00311-1.View ArticlePubMed
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome biology 2002, 3(7):RESEARCH0034. 10.1186/gb-2002-3-7-research0034.PubMed CentralView ArticlePubMed
- Singh R, Green MR: Sequence-specific binding of transfer RNA by glyceraldehyde-3-phosphate dehydrogenase. Science 1993, 259(5093):365-368. 10.1126/science.8420004.View ArticlePubMed
- Ishitani R, Sunaga K, Hirano A, Saunders P, Katsube N, Chuang DM: Evidence that glyceraldehyde-3-phosphate dehydrogenase is involved in age-induced apoptosis in mature cerebellar neurons in culture. Journal of neurochemistry 1996, 66(3):928-935.View ArticlePubMed
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer research 2004, 64(15):5245-5250. 10.1158/0008-5472.CAN-04-0496.View ArticlePubMed
- Silver N, Best S, Jiang J, Thein SL: Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC molecular biology 2006, 7: 33. 10.1186/1471-2199-7-33.PubMed CentralView ArticlePubMed
- Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC plant biology 2004, 4: 14. 10.1186/1471-2229-4-14.PubMed CentralView ArticlePubMed
- Goossens K, Van Poucke M, Van Soom A, Vandesompele J, Van Zeveren A, Peelman LJ: Selection of reference genes for quantitative real-time PCR in bovine preimplantation embryos. BMC developmental biology 2005, 5: 27. 10.1186/1471-213X-5-27.PubMed CentralView ArticlePubMed
- Ingerslev HC, Pettersen EF, Jakobsen RA, Petersen CB, Wergeland HI: Expression profiling and validation of reference gene candidates in immune relevant tissues and cells from Atlantic salmon ( Salmo salar L.). Molecular immunology 2006, 43(8):1194-1201. 10.1016/j.molimm.2005.07.009.View ArticlePubMed
- Ding J, Jia J, Yang L, Wen H, Zhang C, Liu W, Zhang D: Validation of a rice specific gene, sucrose phosphate synthase, used as the endogenous reference gene for qualitative and real-time quantitative PCR detection of transgenes. Journal of agricultural and food chemistry 2004, 52(11):3372-3377. 10.1021/jf049915d.View ArticlePubMed
- Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochemical and biophysical research communications 2006, 345(2):646-651. 10.1016/j.bbrc.2006.04.140.View ArticlePubMed
- Kim BR, Nam HY, Kim SU, Kim SI, Chang YJ: Normalization of reverse transcription quantitative-PCR with housekeeping genes in rice. Biotechnology letters 2003, 25(21):1869-1872. 10.1023/A:1026298032009.View ArticlePubMed
- Nicot N, Hausman JF, Hoffmann L, Evers D: Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. Journal of experimental botany 2005, 56(421):2907-2914. 10.1093/jxb/eri285.View ArticlePubMed
- Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR: Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis . Plant physiology 2005, 139(1):5-17. 10.1104/pp.105.063743.PubMed CentralView ArticlePubMed
- Garner W, Alland H: Effect of the relative length of day and night and other factors of the environment on growth and reproduction in plants. J Agric Res 1920, 18: 553-606.
- Borthwick H, Parker M: Influence of photoperiods upon the differentiation of meristems and the blossoming of Biloxi soybeans. Bot Gaz 1938, 99: 825-839. 10.1086/334749.View Article
- Munir J, Dorn LA, Donohue K, Schmitt J: The effect of maternal photoperiod on seasonal dormancy in Arabidopsis thaliana (Brassicaceae). Am J Bot 2001, 88(7):1240-1249. 10.2307/3558335.View ArticlePubMed
- Han T, Wu C, Tong Z, Mentreddy RS, Tan K, Gai J: Postflowering photoperiod regulates vegetative growth and reproductive development of soybean. Enviromental and Experimental Botany 2006, 55: 120-129. 10.1016/j.envexpbot.2004.10.006.View Article
- Wu C, Ma Q, Yam KM, Cheung MY, Xu Y, Han T, Lam HM, Chong K: In situ expression of the GmNMH7 gene is photoperiod-dependent in a unique soybean ( Glycine max [L.] Merr.) flowering reversion system. Planta 2006, 223(4):725-735. 10.1007/s00425-005-0130-y.View ArticlePubMed
- Masuda J, Urakawa T, Ozaki Y, Okubo H: Short photoperiod induces dormancy in Lotus ( Nelumbo nucifera ). Annals of botany 2006, 97(1):39-45. 10.1093/aob/mcj008.PubMed CentralView ArticlePubMed
- Kantolic AG, Slafer GA: Photoperiod sensitivity after flowering and seed number determination in indeterminate soybean cultivars. Field Crops Res 2001, 72: 109-118. 10.1016/S0378-4290(01)00168-X.View Article
- Kantolic AG, Slafer GA: Development and seed number in indeterminate soybean as affected by timing and duration of exposure to long photoperiods after flowering. Annals of botany 2007, 99(5):925-933. 10.1093/aob/mcm033.PubMed CentralView ArticlePubMed
- Putterill J, Robson F, Lee K, Simon R, Coupland G: The CONSTANS gene of Arabidopsis promotes flowering and encodes a protein showing similarities to zinc finger transcription factors. Cell 1995, 80(6):847-857. 10.1016/0092-8674(95)90288-0.View ArticlePubMed
- Teper-Bamnolker P, Samach A: The flowering integrator FT regulates SEPALLATA3 and FRUITFULL accumulation in Arabidopsis leaves. Plant Cell 2005, 17(10):2661-2675. 10.1105/tpc.105.035766.PubMed CentralView ArticlePubMed
- Corbesier L, Vincent C, Jang S, Fornara F, Fan Q, Searle I, Giakountis A, Farrona S, Gissot L, Turnbull C, Coupland G: FT protein movement contributes to long-distance signaling in floral induction of Arabidopsis . Science 2007, 316(5827):1030-1033. 10.1126/science.1141752.View ArticlePubMed
- Jung JH, Seo YH, Seo PJ, Reyes JL, Yun J, Chua NH, Park CM: The GIGANTEA-Regulated MicroRNA172 Mediates Photoperiodic Flowering Independent of CONSTANS in Arabidopsis . Plant Cell 2007, 19(9):2736-2748. 10.1105/tpc.107.054528.PubMed CentralView ArticlePubMed
- Toegel S, Huang W, Piana C, Unger FM, Wirth M, Goldring MB, Gabor F, Viernstein H: Selection of reliable reference genes for qPCR studies on chondroprotective action. BMC molecular biology 2007, 8: 13. 10.1186/1471-2199-8-13.PubMed CentralView ArticlePubMed
- Sun H, Liu Y, Hu P, Hou W, Wu C, Cao D, Han T: Cloning and characterization of GmBFT, a soybean BFT homologue encoding the phosphatidylethanolamine-binding protein. DNA Sequcence 2008. (accepted)
- Kobayashi Y, Kaya H, Goto K, Iwabuchi M, Araki T: A pair of related genes with antagonistic roles in mediating flowering signals. Science 1999, 286(5446):1960-1962. 10.1126/science.286.5446.1960.View ArticlePubMed
- Walker NJ: Tech. Sight. A technique whose time has come. Science 2002, 296(5567):557-559. 10.1126/science.296.5567.557.View ArticlePubMed
- Ohl F, Jung M, Xu C, Stephan C, Rabien A, Burkhardt M, Nitsche A, Kristiansen G, Loening SA, Radonic A, Jung K: Gene expression studies in prostate cancer tissue: which reference gene should be selected for normalization? Journal of molecular medicine (Berlin, Germany) 2005, 83(12):1014-1024.View Article
- Steinau M, Rajeevan MS, Unger ER: DNA and RNA references for qRT-PCR assays in exfoliated cervical cells. J Mol Diagn 2006, 8(1):113-118. 10.2353/jmoldx.2006.050088.PubMed CentralView ArticlePubMed
- Reid KE, Olsson N, Schlosser J, Peng F, Lund ST: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC plant biology 2006, 6: 27. 10.1186/1471-2229-6-27.PubMed CentralView ArticlePubMed
- Hochstrasser M: Evolution and function of ubiquitin-like protein-conjugation systems. Nature cell biology 2000, 2(8):E153-157. 10.1038/35019643.View ArticlePubMed
- Filby AL, Tyler CR: Appropriate 'housekeeping' genes for use in expression profiling the effects of environmental estrogens in fish. BMC molecular biology 2007, 8: 10. 10.1186/1471-2199-8-10.PubMed CentralView ArticlePubMed
- Ohl F, Jung M, Radonic A, Sachs M, Loening SA, Jung K: Identification and validation of suitable endogenous reference genes for gene expression studies of human bladder cancer. The Journal of urology 2006, 175(5):1915-1920. 10.1016/S0022-5347(05)00919-5.View ArticlePubMed
- Han T, Wang J: Studies on the post-flowering photoperiodic responses in soybean. Acta Botanica Sinica 1995, 37(11):863-869.
- Soybean EST database[http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=soybean]
- TAIR BLAST[http://www.arabidopsis.org/Blast/]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.