Defining reference genes in Oryza sativausing organ, development, biotic and abiotic transcriptome datasets
© Narsai et al; licensee BioMed Central Ltd. 2010
Received: 19 November 2009
Accepted: 31 March 2010
Published: 31 March 2010
Reference genes are widely used to normalise transcript abundance data determined by quantitative RT-PCR and microarrays. However, the approaches taken to define reference genes can be variable. Although Oryza sativa (rice) is a widely used model plant and important crop specie, there has been no comprehensive analysis carried out to define superior reference genes.
Analysis of 136 Affymetrix transcriptome datasets comprising of 373 genome microarrays from studies in rice that encompass tissue, developmental, abiotic, biotic and hormonal transcriptome datasets identified 151 genes whose expression was considered relatively stable under all conditions. A sub-set of 12 of these genes were validated by quantitative RT-PCR and were seen to be stable under a number of conditions. All except one gene that has been previously proposed as a stably expressed gene for rice, were observed to change significantly under some treatment.
A new set of reference genes that are stable across tissue, development, stress and hormonal treatments have been identified in rice. This provides a superior set of reference genes for future studies in rice. It confirms the approach of mining large scale datasets as a robust method to define reference genes, but cautions against using gene orthology or counterparts of reference genes in other plant species as a means of defining reference genes.
The analysis of gene expression, or more correctly transcript abundance, is widely carried out in a variety of laboratories in various disciplines. Northern blotting, quantitative RT-PCR (QRT-PCR) and microarray approaches are commonly used to assess transcript abundance. All these approaches need a standard or reference for comparison, so that the changes observed can be attributed to a biological process rather than an artefact of the particular technique used [1, 2]. The use of northern blotting often involves the use of equal RNA (total or mRNA) loading as a reference point. Although this can lead to errors, the variability of many steps in northern blotting means that northern blots are generally only used to assess large changes in transcript abundance. In contrast, microarray analysis assesses the transcript abundance of tens of thousands of genes, thus it has required the application of statistical methods to normalise the distribution of signals and also requires correction for large samples sets, so called false discovery rate correction [3, 4]. For QRT-PCR analysis, house-keeping or reference genes can be used as a standard and by definition; the transcript abundance of this gene should not change under the experimental conditions being studied.
The definition of reference genes is important as the use of common sets of reference genes by scientists allows direct comparisons between studies. The benefits of comparing transcripts abundance datasets between a variety of studies is best exemplified with microarray studies, where the predominant use of a single robust platform for studies in Arabidopsis thaliana has led to the development of a number of databases where in silico or digital northern analyses can be carried out. Thus, databases such as Genevestigator  and the Botany Array Resource (BAR)  are just two examples that provide a valuable resource for researchers to obtain information of transcript abundance patterns for genes of interest.
QRT-PCR is often used to validate transcriptome data obtained from array studies or is used in more directed studies where the transcript abundance of a limited number of genes is analysed. Increasingly large scale studies encompassing several hundred to thousands of genes are also analysed by QRT-PCR and represent an important resource to the scientific community, e.g. expression profiling of transcription factors [7–9]. Thus, accurate reference genes are required to interpret such data. In an Arabidopsis study that defined stably expressed genes under a wide variety of conditions and organs, a "superior set" of reference genes were identified that are widely used in QRT-PCR studies in Arabidopsis . An alternative approach to define reference genes is the use of various statistical tests that essentially rank the variability of transcripts abundances for sets of genes that are analysed . Bestkeeper , Norm-Finder  and geNORM  are examples of such widely used programs, albeit their use is limited to some extent in studies with plants .
A variety of studies in different plant species have defined reference genes . Many studies selected a number of potential reference genes based on what is used in other plant species, and tested changes in transcript abundance, using statistical algorithms outlined above to test for variations in different organs or environmental conditions, to determine their suitability as reference genes [14–17]. All these studies have defined reference genes, but the limited number of conditions tested and the lack of genome wide searches for superior reference genes means that these sets may not represent the best reference genes under a wide variety of conditions. The ability of software programs to define variations in gene expression is limited by the input data. However, it is desirable to define reference genes that are stable in transcript abundance under as many conditions as possible and analysing as many genes in the genome as possible.
Oryza sativa (rice) represents an important model plant  and as a crop, provides 21% of the calorie needs of the world's population (and up to ~75% for the population of south east Asia . As such, it is the focus of intense research by a wide variety of researchers. One of the fundamental problems facing researchers carrying out gene expression studies is the use of control or reference genes that should not change, preferably under all experimental conditions. Reference genes in rice have been proposed by testing commonly used reference genes in plants and orthologues of reference genes that have been defined as in Arabidopsis [7, 20]. It is unclear under how many different parameters these genes are appropriate reference genes and also if superior reference genes could be defined using a genome wide approach as previously carried out in Arabidopsis .
In order to define suitable reference genes in rice in an objective manner, a similar procedure to that used to define reference genes in Arabidopsis was undertaken . We collated 373 Affymetrix genome arrays from rice that encompassed tissue, abiotic, biotic and hormonal parameters to define a set of 151 probesets that were stably expressed under all conditions. Of these, 12 genes were chosen as reference genes and validated using QRT-PCR, for different tissues and under stress. In this way, a superior set of reference genes for rice was identified that are suitable for organ, development and stress based experiments.
Results and Discussion
Selection of transcriptome datasets
Overview of experiments involving 373 Affymetrix rice genome microarrays used for the global analysis in this study.
Dry seed and aerobic germination (up to 24 h) cv. Amaroo
Dry and germinating seed
Dry seed and anaerobic germination (up to 24 h) and switch conditions cv. Amaroo
Aerobic and anaerobic grown coleoptiles cv. Nipponbare
Embryo, endosperm, leaf and root from 7-d seedling, 10-d seedling cv. Zhonghua
Embryo, endosperm, leaf and root from 7-d seedling, 10-d seedling
Stigma, Ovary+7 single arrays cv. Nipponbare
Stigma, ovary+7 single arrays
Mature leaf, young leaf, semi apical meristem, inflorescence, seed cv. IR64
Mature leaf, young leaf, semi apical meristem, inflorescence, seed
Drought, salt, cold stress cv. IR64
Heat stress cv. Zhonghua
Salt stress on 2 cultivars; indica, FL478 (salt tolerant), indica, IR29 (salt sensitive)
Crown and growing point
Salt stress on 4 cultivars; japonica, m103 (salt sensitive), indica, IR29 (salt sensitive), japonica, Agami (salt tolerant), indica, IR63731 (salt tolerant)
Panicle initiation stage
Salt stress on root using 4 cultivars; FL478 (salt tolerant), IR29 (salt sensitive), IR63731 (salt tolerant), Pokkali (salt tolerant)
Fe and P treatments cv. Nipponbare
Arsenate treatment cv. Azucena
Physical stress at roots tips cv. Bala
S.Hermonthica plant parasite infection cv. Nipponbare (resistant), IAC165 (susceptible)
M.grisea blast fungus infection cv. Nipponbare
Rice stripe virus infection cv. WuYun3, KT95-418
Infection with bacteria X.Oryzae pv. oryzicola and oryzae cv. Nipponbare
Cytokinin treatment on root and leaf cv. Nipponbare
Root, 2-week old seedlings
Indole-3-actetic acid and benzyl aminopurine treatment cv. IR64
Global analysis of transcriptome datasets
Selection of reference genes
The GC-RMA normalised data for all microarrays with publically available CEL files (331 microarrays; Table 1) was used to calculate the mean, standard deviation (SD) and coefficient of variance (CV; CV = SD/mean) for all 7,922 probesets, where a low CV is indicative of lower variation. This was followed by selection process undertaken to determine which of these genes were suitable as reference genes (Figure 1). Only 151 of the 7,922 probesets were defined as stably expressed across the developmental, stress and/or entire dataset (Figure 1).
The list of reference genes for rice, defined in this and previous studies.
Probe Set Identifier
Nucleic acid binding protein
Tumor protein homolog
Tumor protein homolog
No TIGR identifier
Eukaryotic initiation factor 5C
Translation factor SUI1
SKP1-like protein 1A
Splicing factor U2af
; commonly used reference gene
; commonly used reference gene
; commonly used reference gene
; commonly used reference gene
; commonly used reference gene
; commonly used reference gene
; commonly used reference gene
Nucleotide tract-binding protein
AP-2 complex subunit
Closer analysis of these 12 genes reveals that the genes encoding, a 3-phosphoinositide-dependent protein kinase-1 (LOC_Os06g48970.1) and a nucleic acid binding protein (LOC_Os06g11170.1) showed stable, moderate expression levels across the stress microarray set (Genes 1-2 in Table 2; Figure 2B). While the genes encoding a tumor protein homolog (LOC_Os11g43900.1) and translation initiation factor SUI1 (LOC_Os07g34589.1) showed stable expression across the developmental and entire microarray sets respectively (Genes 4, 7 in Table 2; Figure 2B). As would be expected, it can be seen that many of these stably expressed genes are involved in core cellular functions such as mRNA splicing and translation initiation (Genes 1-12 denoted by blue asterisks in Figure 2A; 2B; Table 2).
Validation of reference genes in quantitative RT-PCR in tissue and stress samples
In order to confirm stable expression of the reference genes identified in this study primers were designed to 26 genes,12 stably expressed genes identified in this study and 14 previously suggested reference genes (Table 1, Additional file 1, Table S1). The stability of transcript abundance of these genes was analysed by QRT-PCR across 15 different samples from a variety of developmental (dry seed, imbibed seed, leaf and roots from young and old plants) and stress treated tissues (shoots from cold treated and heat treated young seedlings over time; Materials and methods). High quality total RNA was isolated from these samples and reverse transcribed to generate cDNA. The same cDNA pool from each of the samples was used to measure the transcript abundance by QRT-PCR, with melt curve analysis for each gene confirming primer specificity.
The geNORM v3.5 software was used to analyse the expression stability for the reference genes analysed by QRT-PCR from the 12 tissue samples (Additional file 1, Table S1) . This software allows calculation of a gene stability measure (M) value for all the genes analysed, where genes with the lowest M value shown the most stable expression (Figure 3A). Authors of the geNORM software suggest using the 3 most stable genes (3 lowest M values) as the most appropriate reference genes . It can be seen that even when commonly or previously suggested reference genes and the novel reference genes from this study are analysed together, all 3 of the most stable genes are the novel reference genes identified in this study (Figure 3A). It is important to note that this M value is only calculated based on data from the limited number of samples that were analysed by QRT-PCR, thus not representing the wide variety of tissues/treatments analysed by microarrays. Therefore, in order to visualise the variation in expression across in the microarrays in parallel, the CV values for each gene was also plotted with the M values, where a lower CV value indicates greater stability. In this way, the most stable genes were identified as those with both low M and CV values. In this combined analysis, the 12 genes chosen all outperformed previously used reference genes, particularly in terms of having a lower CV (Figure 3A), the genes indicated with a black diamond all had lower CV values as indicated by the bar graph, with a gene encoding a nucleic acid binding protein (LOC_Os06g11170.1) apparently the most stable (Figure 3A).
To further test the stability of the reference genes defined in this study, the expression of the 12 novel reference genes defined in this study were analysed independently by geNORM for the samples analysed by QRT-PCR (Figure 3B and 3C). Overall, it can be seen that the most stable genes had low M values as well and low CV values, indicating stable expression (Figure 3B). Furthermore, the geNORM pair-wise analysis to determine the number of control genes recommended for use in normalisation , revealed that 2 or even one gene is stable enough for accurate normalisation, however 2 genes is recommended for more robust normalisation (V < 0.15; Additional file 2, Figure S1) . Using QRT-PCR analysis, we also compared the expression of these 12 reference genes to 3 heat or cold responsive genes including, an Apetala type transcription factor (AP2), a heat shock responsive factor (HSF-82) and alternative oxidase (AOX) over time under i) cold or ii) heat conditions respectively (Figure 3C). It can be observed that under cold treatment, all 12 reference genes show very stable expression over time (Figure 3Ci). Similarly, despite slight variation of some genes under heat conditions, it is evidenced that overall, these genes are also stably expressed over time following heat treatment (Figure 3Cii). In addition, the observed induction of AP2 and HSF-82 under cold and heat treatment, confirmed the success of the respective treatments (Figure 3C). Furthermore, comparison of this induction (at 3 h) to the induction observed from the analogous microarray data, showed that normalisation of the QRT-PCR data using the reference genes defined in this study resulted in comparable increases to those seen using the microarray data (Figure 3C).
Comparison to previous studies and other expression platforms
A large-scale study of reference genes in Arabidopsis revealed superior reference genes using Affymetrix microarray data . Using the Inparanoid orthologue output  for Arabidopsis and rice, it was seen that only 15 rice orthologues of the 30 novel Arabidopsis reference genes were also expressed across all the microarrays in this study and 3 of these were randomly selected for further analysis by QRT-PCR (Genes 24-26; Table 2). Notably, only 1 gene (LOC_Os03g05290.1) encoding an aquaporin TIP protein, was seen to be stably expressed i.e. one of the 151 stably expressed in this study (red asterisk only; Figure 2A). It may be noted that the overall CV values are higher in this study compared to the CV values calculated in the Arabidopsis study . The main reason for this is likely to be due to significant differences in the variability of the input data from both studies. That is, the Arabidopsis reference gene study used microarray data generated from only 7 studies using a large number of microarrays each e.g. 237 microarrays in the single developmental study , whilst this study involved analysis of microarrays from 20 studies carried out in different laboratories, using between 4 and 60 microarrays in each.
Previous studies in rice have examined reference genes using QRT-PCR analysis, however these only involved analysis of a small number of commonly used reference genes such as Actin, Actin1, alpha and beta tubulin, polyubiquitin, ubiquitin, GAPDH and elongation factor 1 in up to 25 samples, under a limited range of conditions [7, 20]. Analysis of these genes in the context of this study (Genes 13-20; Table 2) revealed that some of these were not detected as expressed in one or more tissue/stress microarray experiments, notably, this included Actin1 (LOC_Os05g36290.1; Gene 14 in Table 2) which was not expressed in all 3 biological replicates of the semi apical meristem (GSE6901) (Figure 2C). Similarly, a recent study in rice defined a set of 248 stably expressed genes across 40 developmental tissues that were analysed using Yale/BGI oligonucleotide microarrays . Only 61 of these genes were found to be expressed across all the microarrays analysed in this study, nevertheless 3 of these were randomly selected for further analysis by QRT-PCR (Genes 21-23; Table 2). Notably, one of the 61 genes (LOC_Os07g02340.1) encoding an "expressed protein" was also found to fulfil all the criteria outlined in Figure 1, and showed stable expression across all the samples analysed in the present study (Gene 8 in Table 2; denoted by red and blue asterisk in Figure 2A and 2B).
The use of the large datasets of rice microarray data has provided identification of sets of genes that are stably expressed under a wide variety of parameters. Although microarray platforms were not designed to be quantitative, direct comparison of over 1000 QRT-PCR assays with microarray data has revealed a high degree of correlation . This is consistent with the use of microarray data to define superior reference genes as outlined here, and previously in Arabidopsis . Based on these principles, we suggest the use of one or more of the novel reference genes presented in this study for the normalisation of rice microarray or QRT-PCR data. However although the reference genes identified in this study are stable under a wide variety of parameters, such as developmental, tissue and various stresses, it is essential that each study validate the stability of the selected reference gene(s) to achieve the systematic validation of reference genes that is required to compare different studies .
Analyses of all publically available rice microarrays
To compile the entire publically available Affymetrix rice microarray (as at 1st August 2009), all experiments containing CEL files were downloaded from the Gene Expression Omnibus within the National Centre for Biotechnology Information database or from the MIAME ArrayExpress database http://www.ebi.ac.uk/arrayexpress/. The GSE or EXP numbers for the respective rice studies are shown in Table 1. There was a total of 373 microarrays for which there was either MAS5.0 data available, thus all of these were used for present/absent determination in defining the list of 7,922 probesets expressed in all microarrays. However of these 373 microarrays, 7 had no biological replicates and 35 did not have available CEL files, thus the remaining 331 microarrays were used to carry out further normalisation (GC-RMA) and calculation of mean, standard deviation and coefficient of variance (CV). This allowed analysis of 117 tissues/conditions, with a minimum of 2 biological replicates. The 117 included 41 organ/developmental tissues, 65 samples within abiotic and biotic stress experiments and 11 samples within hormone treatment experiments.
All raw intensity CEL files were imported into Avadis 4.3 (Strand Genomics) and the standard MAS5.0 normalisation was first carried out in order to determine present/absent/marginal calls for each probeset. For all 331 microarrays with available CEL files (and carried out for biological replicates), GC-RMA normalisation was carried out. The mean expression, SD and CV (=SD/mean) was then calculated for each of the 7,922 probesets across the developmental set, stress set and entire dataset (which included the hormone experiments). On the basis of CV cut-offs, the list of 151 probesets was generated (Figure 1). The averaged log GC-RMA normalised values for these 151 probesets, across the developmental tissues, stress and hormone treatment experiments were hierarchically clustered using average linkage on Euclidean distance. The clustering analysis and heatmap generation was carried out using Partek Genomics Suite, version 6.3 (Partek). For the Agilent microarray comparison, data was retrieved under the accession GSE8518 from the Gene Expression Omnibus within the National Centre for Biotechnology Information database.
Analysis of orthologues
The InParanoid: Eukaryotic Orthologue Groups database (version 7.0) was used to analyse all orthologues between rice and Arabidopsis . The orthologous group file was downloaded for the whole-genome comparison of rice versus Arabidopsis. This produced information for orthologues identified by TIGR identifiers for rice and AGIs for Arabidopsis.
Stress treatments, tissue collection and RNA isolation
In order to analyse the expression of all the genes in Table 2, a selection of tissues were collected across different developmental stages/tissues and under different stress conditions in wild type rice, cv. Amaroo. In order to analyse different developmental tissues; embryos were extracted from dry seed, seeds imbibed for 24 h with oxygen gas (24 h A), seeds imbibed for 24 h in the absence of oxygen gas i.e. in the presence of nitrogen gas (24 h N), seeds imbibed for 24 h under nitrogen gas and switched to oxygen gas for 3 h (27 NA), leaf and root tissues from 2-week old seedlings and 3 month old plants. Furthermore, to examine the effects of abiotic stress, 2-week old seedlings were transferred to 4°C and 42°C for cold and heat treatment respectively over a 9 h time course, whilst the controls remained at a constant temperature of 30°C. All 15 tissue samples were analysed using three biological replicates, the RNA was isolated using the Qiagen RNeasy Plant RNA isolation kit and DNase treated using both the Qiagen on-column DNase digestion as well as the Ambion Turbo DNase treatment exactly as carried out in Howell et al., 2009 .
Details of the primer sequences and amplicon lengths for each of the genes are shown in Additional file 1, Table S1. The transcript abundance for each gene was analysed using the SYBR green I master (Roche, Sydney) with the Roche LC480. Each sample was analysed in biological triplicate, using individual plants and treatments to test for reproducibility. Following RNA isolation each of the samples was quantitated using a Nanodrop spectrophotometer. This provided the following information for each sample: concentration (ng/μl), the absorbance (A) in nm at 230, 260 and 280, the A230/A260 and A260/A280 ratios. Using this information the RNA yield and purity was calculated to ensure that they all had no significant impurities between samples that may affect reverse transcription and/or amplification during QRT-PCR. 1 μg of total RNA was reverse transcribed using the Bio-Rad® (Sydney) iScript reverse transcription kit, according manufacturer's instructions. In parallel for each sample, another 1 μg of RNA was used for the same reverse transcription reaction, with the exception of the addition of the reverse transcriptase enzyme (no RT samples). Following this, the Qiagen® PCR purification kit was used according to manufacturer's instructions on all samples (RT and "no RT" samples). This purified cDNA was diluted 1 in 10 with nuclease-free water and 1 μl was used for QRT-PCR analysis. For the no RT samples, no dilution was carried out and 1 μl was used in the same manner as the diluted cDNA for QRT-PCR analysis, this enables the detection of any genomic DNA contamination.
coefficient of variance
rice massively parallel signature sequencing
elongation factor 1 delta
heat shock factor
This work was supported by an Australian Research Council Centre of Excellence Grant CEO561495.
- Guenin S, Mauriat M, Pelloux J, Van Wuytswinkel O, Bellini C, Gutierrez L: Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J Exp Bot. 2009, 60 (2): 487-493. 10.1093/jxb/ern305.PubMedView ArticleGoogle Scholar
- Gutierrez L, Mauriat M, Guenin S, Pelloux J, Lefebvre JF, Louvet R, Rusterucci C, Moritz T, Guerineau F, Bellini C, Van Wuytswinkel O: The lack of a systematic validation of reference genes: a serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnol J. 2008, 6 (6): 609-618. 10.1111/j.1467-7652.2008.00346.x.PubMedView ArticleGoogle Scholar
- Nettleton D: A discussion of statistical methods for design and analysis of microarray experiments for plant scientists. Plant Cell. 2006, 18 (9): 2112-2121. 10.1105/tpc.106.041616.PubMedPubMed CentralView ArticleGoogle Scholar
- Quackenbush J: Microarray data normalization and transformation. Nat Genet. 2002, 32 (Suppl): 496-501. 10.1038/ng1032.PubMedView ArticleGoogle Scholar
- Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136 (1): 2621-2632. 10.1104/pp.104.046367.PubMedPubMed CentralView ArticleGoogle Scholar
- Toufighi K, Brady SM, Austin R, Ly E, Provart NJ: The Botany Array Resource: e-Northerns, Expression Angling, and promoter analyses. Plant J. 2005, 43 (1): 153-163. 10.1111/j.1365-313X.2005.02437.x.PubMedView ArticleGoogle Scholar
- Caldana C, Scheible WR, Mueller-Roeber B, Ruzicic S: A quantitative RT-PCR platform for high-throughput expression profiling of 2500 rice transcription factors. Plant Methods. 2007, 3: 7. 10.1186/1746-4811-3-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Czechowski T, Bari RP, Stitt M, Scheible WR, Udvardi MK: Real-time RT-PCR profiling of over 1400 Arabidopsis transcription factors: unprecedented sensitivity reveals novel root- and shoot-specific genes. Plant J. 2004, 38 (2): 366-379. 10.1111/j.1365-313X.2004.02051.x.PubMedView ArticleGoogle Scholar
- Kakar K, Wandrey M, Czechowski T, Gaertner T, Scheible WR, Stitt M, Torres-Jerez I, Xiao Y, Redman JC, Wu HC, Cheung F, Town CD, Udvardi MK: A community resource for high-throughput quantitative RT-PCR analysis of transcription factor gene expression in Medicago truncatula. Plant Methods. 2008, 4: 18. 10.1186/1746-4811-4-18.PubMedPubMed CentralView ArticleGoogle Scholar
- 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 Physiol. 2005, 139 (1): 5-17. 10.1104/pp.105.063743.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26 (6): 509-515. 10.1023/B:BILE.0000019559.84305.47.PubMedView ArticleGoogle Scholar
- 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 Res. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- 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 Biol. 2002, 3 (7): RESEARCH0034. 10.1186/gb-2002-3-7-research0034.PubMedPubMed CentralView ArticleGoogle Scholar
- Barsalobres-Cavallari CF, Severino FE, Maluf MP, Maia IG: Identification of suitable internal control genes for expression studies in Coffea arabica under different experimental conditions. BMC Mol Biol. 2009, 10: 1. 10.1186/1471-2199-10-1.PubMedPubMed CentralView ArticleGoogle Scholar
- Exposito-Rodriguez M, Borges AA, Borges-Perez A, Perez JA: Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant Biol. 2008, 8: 131. 10.1186/1471-2229-8-131.PubMedPubMed CentralView ArticleGoogle Scholar
- Hu R, Fan C, Li H, Zhang Q, Fu YF: Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR. BMC Mol Biol. 2009, 10: 93. 10.1186/1471-2199-10-93.PubMedPubMed CentralView ArticleGoogle Scholar
- Reid KE, Olsson N, Schlosser J, Oenq 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 Biol. 2006, 6: 27-37. 10.1186/1471-2229-6-27.PubMedPubMed CentralView ArticleGoogle Scholar
- Flavell R: Role of model plant species. Methods Mol Biol. 2009, 513: 1-18. full_text.PubMedView ArticleGoogle Scholar
- Fitzgerald MA, McCouch SR, Hall RD: Not just a grain of rice: the quest for quality. Trends Plant Sci. 2009, 14 (3): 133-139. 10.1016/j.tplants.2008.12.004.PubMedView ArticleGoogle Scholar
- 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. Biochem Biophys Res Commun. 2006, 345 (2): 646-651. 10.1016/j.bbrc.2006.04.140.PubMedView ArticleGoogle Scholar
- Narsai R, Howell KA, Carroll A, Ivanova A, Millar AH, Whelan J: Defining core metabolic and transcriptomic responses to oxygen availability in rice embryos and young seedlings. Plant Physiol. 2009, 151 (1): 306-322. 10.1104/pp.109.142026.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiao Y, Tausta SL, Gandotra N, Sun N, Liu T, Clay NK, Ceserani T, Chen M, Ma L, Holford M, Zhang HY, Zhao H, Deng XW, Nelson T: A transcriptome atlas of rice cell types uncovers cellular, functional and developmental hierarchies. Nat Genet. 2009, 41 (2): 258-263. 10.1038/ng.282.PubMedView ArticleGoogle Scholar
- Remm M, Storm CE, Sonnhammer EL: Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001, 314 (5): 1041-1052. 10.1006/jmbi.2000.5197.PubMedView ArticleGoogle Scholar
- Mosquera G, Giraldo MC, Khang CH, Coughlan S, Valent B: Interaction transcriptome analysis identifies Magnaporthe oryzae BAS1-4 as Biotrophy-associated secreted proteins in rice blast disease. Plant Cell. 2009, 21 (4): 1273-1290. 10.1105/tpc.107.055228.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Aken O, Zhang B, Carrie C, Uggalla V, Paynter E, Giraud E, Whelan J: Defining the Mitochondrial Stress Response in Arabidopsis thaliana. Mol Plant. 2009, 2 (6): 1310-1324. 10.1093/mp/ssp053.PubMedView ArticleGoogle Scholar
- Canales RD, Luo Y, Willey JC, Austermiller B, Barbacioru CC, Boysen C, Hunkapiller K, Jensen RV, Knight CR, Lee KY, Ma Y, Maqsodi B, Papallo A, Peters EH, Poulter K, Ruppel PL, Samaha RR, Shi L, Yang W, Zhang L, Goodsaid FM: Evaluation of DNA microarray results with quantitative gene expression platforms. Nat Biotechnol. 2006, 24 (9): 1115-1122. 10.1038/nbt1236.PubMedView ArticleGoogle Scholar
- Howell KA, Narsai R, Carroll A, Ivanova A, Lohse M, Usadel B, Millar AH, Whelan J: Mapping metabolic and transcript temporal switches during germination in rice highlights specific transcription factors and the role of RNA instability in the germination process. Plant Physiol. 2009, 149 (2): 961-980. 10.1104/pp.108.129874.PubMedPubMed CentralView ArticleGoogle Scholar
- Lasanthi-Kudahettige R, Magneschi L, Loreti E, Gonzali S, Licausi F, Novi G, Beretta O, Vitulli F, Alpi A, Perata P: Transcript profiling of the anoxic rice coleoptile. Plant Physiol. 2007, 144 (1): 218-231. 10.1104/pp.106.093997.PubMedPubMed CentralView ArticleGoogle Scholar
- Xue LJ, Zhang JJ, Xue HW: Characterization and expression profiles of miRNAs in rice seeds. Nucleic Acids Res. 2009, 37 (3): 916-930. 10.1093/nar/gkn998.PubMedPubMed CentralView ArticleGoogle Scholar
- Li M, Xu W, Yang W, Kong Z, Xue Y: Genome-wide gene expression profiling reveals conserved and novel molecular functions of the stigma in rice. Plant Physiol. 2007, 144 (4): 1797-1812. 10.1104/pp.107.101600.PubMedPubMed CentralView ArticleGoogle Scholar
- Jain M, Nijhawan A, Arora R, Agarwal P, Ray S, Sharma P, Kapoor S, Tyagi AK, Khurana JP: F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol. 2007, 143 (4): 1467-1483. 10.1104/pp.106.091900.PubMedPubMed CentralView ArticleGoogle Scholar
- Hu W, Hu G, Han B: Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Science. 2009, 176: 583-590. 10.1016/j.plantsci.2009.01.016.PubMedView ArticleGoogle Scholar
- Walia H, Wilson C, Condamine P, Liu X, Ismail AM, Zeng L, Wanamaker SI, Mandal J, Xu J, Cui X, Close TJ: Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage. Plant Physiol. 2005, 139 (2): 822-835. 10.1104/pp.105.065961.PubMedPubMed CentralView ArticleGoogle Scholar
- Walia H, Wilson C, Zeng L, Ismail AM, Condamine P, Close TJ: Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage. Plant Mol Biol. 2007, 63 (5): 609-623. 10.1007/s11103-006-9112-0.PubMedPubMed CentralView ArticleGoogle Scholar
- Zheng L, Huang F, Narsai R, Wu J, Giraud E, He F, Cheng L, Wang F, Wu P, Whelan J, Shou H: Physiological and transcriptome analysis of iron and phosphorus interaction in rice seedlings. Plant Physiol. 2009, 151 (1): 262-74. 10.1104/pp.109.141051.PubMedPubMed CentralView ArticleGoogle Scholar
- Norton GJ, Lou-Hing DE, Meharg AA, Price AH: Rice-arsenate interactions in hydroponics: whole genome transcriptional analysis. J Exp Bot. 2008, 59 (8): 2267-2276. 10.1093/jxb/ern097.PubMedPubMed CentralView ArticleGoogle Scholar
- Swarbrick PJ, Huang K, Liu G, Slate J, Press MC, Scholes JD: Global patterns of gene expression in rice cultivars undergoing a susceptible or resistant interaction with the parasitic plant Striga hermonthica. New Phytol. 2008, 179 (2): 515-529. 10.1111/j.1469-8137.2008.02484.x.PubMedView ArticleGoogle Scholar
- Ribot C, Hirsch J, Balzergue S, Tharreau D, Notteghem JL, Lebrun MH, Morel JB: Susceptibility of rice to the blast fungus, Magnaporthe grisea. J Plant Physiol. 2008, 165 (1): 114-124. 10.1016/j.jplph.2007.06.013.PubMedView ArticleGoogle Scholar
- Hirose N, Makita N, Kojima M, Kamada-Nobusada T, Sakakibara H: Overexpression of a type-A response regulator alters rice morphology and cytokinin metabolism. Plant Cell Physiol. 2007, 48 (3): 523-539. 10.1093/pcp/pcm022.PubMedView ArticleGoogle Scholar
- Jain M, Khurana JP: Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. Febs J. 2009, 276 (11): 3148-3162. 10.1111/j.1742-4658.2009.07033.x.PubMedView ArticleGoogle Scholar
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.