Skip to main content


We’d like to understand how you use our websites in order to improve them. Register your interest.

Assessment of reference gene stability in Rice stripe virus and Rice black streaked dwarf virus infection rice by quantitative Real-time PCR



Stably expressed reference gene(s) normalization is important for the understanding of gene expression patterns by quantitative Real-time PCR (RT-qPCR), particularly for Rice stripe virus (RSV) and Rice black streaked dwarf virus (RBSDV) that caused seriously damage on rice plants in China and Southeast Asia.


The expression of fourteen common used reference genes of Oryza sativa L. were evaluated by RT-qPCR in RSV and RBSDV infected rice plants. Suitable normalization reference gene(s) were identified by geNorm and NormFinder algorithms.


UBQ 10 + GAPDH and UBC + Actin1 were identified as suitable reference genes for RT-qPCR normalization under RSV and RBSDV infection, respectively. When using multiple reference genes, the expression patterns of OsPRIb and OsWRKY, two virus resistance genes, were approximately similar with that reported previously. Comparatively, by using single reference gene (TIP41-Like), a weaker inducible response was observed.


We proposed that the combination of two reference genes could obtain more accurate and reliable normalization of RT-qPCR results in RSV- and RBSDV-infected plants. This work therefore sheds light on establishing a standardized RT-qPCR procedure in RSV- and RBSDV-infected rice plants, and might serve as an important point for discovering complex regulatory networks and identifying genes relevant to biological processes or implicated in virus.


Rice viral diseases are major threats to rice production and have been distributed worldwide across regions depending on rice cultivation [1]. Two of the most prevalent rice viruses are RSV and RBSDV, which were transmitted by a small brown planthopper (SBPH, Laodelphax striatellus Fallen) [24]. When infected with RSV at the seedling stage, normally, rice plants grow poorly and often develop folded and twisted leaves, with the central leaves yellowing and withering; and plant growth may terminate and ultimately the plant will die [57]. In China, rice stripe is very serious, especially in Jiangsu province, where about 0.6 M ha per year of rice were infected by RSV during the period of 2000 to 2003, increasing to 1 M ha in 2004. In heavily infected fields, rice yield is reduced by 30–50 %, and in some of the most severely infected fields, no harvest is possible [8]. In RBSDV infected rice always develops stunted stems, dark green, twisted leaves, and white waxy swellings along veins on the abaxial surface of the leaves [3, 9, 10]. The disease caused severely damage on rice in most parts of eastern China with due to widespread release of susceptible cultivars. Since the infection damage was very severe in China and Southeast Asia, the understanding of the responses of rice to viral infection, especially gene expression analysis, is very important for developing strategies for disease control [8, 11].

The widely used method to measure transcript abundance is RT-qPCR compared to reverse transcription-polymerase chain reaction (RT-PCR) and northern blot [1214]. Besides being a powerful tool, RT-qPCR suffers from certain pitfalls, most important being the normalization with a reference gene [1517]. In recent years, the reference genes, such as those encoding actin (Actin), tubulin (TUB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and 18S rRNA, are often separately chosen for the normalization in RT-qPCR because of their constant expression levels in living organisms [16]. Nevertheless, different studies sometimes proved different or even opposing results of these reference genes, and it was demonstrated that the transcript levels of these genes actually vary under different experimental conditions [1821]. For example, different expression levels normalized by a different reference gene could be approximately 100 folds [22]. Furthermore, Myzus persicae’s actin and GAPDH protein were found to interaction with Beet western yellows virus in vitro [23] and some A.pisumwere’s genes (Actin and GAPDH) considered to be potentially related to the transmission of Peaenation mosaic virus and Soybean dwarf virus [24]. Thus, it is important and necessary to select suitable reference gene(s) for different experimental paradigms, particularly in RSV and RBSDV infection conditions which those appropriate internal reference(s) were not identified [2527].

In this study, we reported the validation of reference genes to identify the most suitable internal control gene(s) for the normalization of RT-qPCR data upon viral infection in rice plants. Using statistical algorithms geNorm and Norm Finder [28, 29], the stability of 14 candidate reference genes (Actin, UBC, 18S rRNA, EF-1α, UBQ 5, GAPDH, α-TUB, β-TUB, eIF-4α, Actin1, UBQ 10, TIP41-like, EXP and Os AOC) was examined and compared. Two best reference genes were identified more stably expressed than traditional ones in RSV- and RBSDV-infected treatments. Our results further indicated that the combination of these two reference genes provides a good starting point for gene expression analysis in rice viral infection plants by RT-qPCR.


Rice infectivity assay

We characterized the phenotype of RSV- and RBSDV-infected rice plants, and the symptoms were allowed to develop under the controlled environmental conditions. RSV-infected rice developed folded and twisted leaves, with the central leaves yellowing and withering. Meanwhile in RBSDV-infected rice plants, it developed stunted stems and white waxy swellings along veins on the abaxial surface of the leaves (Fig. 1a). RSV and RBSDV were respectively detected in inoculated rice, using RT-PCR (Fig. 1b).

Fig. 1

RSV and RBSDV infection in rice plants. a Disease symptoms developed on rice during virus infection. From left to right, non-inoculated plant, RSV- and RBSDV-infected rice plants. b Agarose gel (2 %) showing RSV and RBSDV in inoculated rice (1–3 means three replications)

Identification of candidate reference genes

In order to evaluate the expression stability of reference genes, the candidates should be identified first. In our test, primers were designed for the fourteen commonly used reference genes, identified by BLASTN and TBLASTN searches at National Center for Biotechnology Information (NCBI). The gene names, accession numbers, gene description and primer sequences were all provided in Table 1. First, melting curve for four representative genes (Dissociation curves for all other genes with single peak were not shown) and agarose gel analysis confirmed that each primer set gave a single amplified product of desired size (Fig. 2a and b). Therefore, these genes were selected for RT-qPCR validation. It is also important to accurately quantitate the quality of RNA before reverse transcription. The concentration and quality of isolated RNA were determined using the NanoDrop 2000 spectrophotometer. Agarose gel electrophoresis assay also confirmed the integrity of RNA samples (Fig. 2c).

Table 1 Candidate reference genes and their primer sequences used in this study
Fig. 2

The specificity of real time PCR amplification. a Dissociation curves for four representative genes with single peak obtained from three technical replicates. b Agarose gel (2 %) showing amplification of a specific PCR product of expected size for each reference gene tested in the study; c Agarose gel showing the quality of RNA (0, 7, 14 and 21 means days after RSV- and RBSDV-infection)

To give an overview picture of the relative abundance of candidate reference genes, the calculated cycle threshold (Ct) values were determined for each gene across all the tested RSV- and RBSDV-infected samples (Fig. 3a). RT-qPCR analysis indicated that 14 candidate reference genes exhibited different levels of abundance. The Ct values ranged from 14 to 30, with the most lying between 22 and 28. 18S rRNA and UBQ 10 were more abundantly transcribed than others, while β-TUB was the least expressed gene across all the tested samples.

Fig. 3

Expression levels of candidate reference genes. a Average cycle threshold (Ct) values for the 14 candidate reference genes used in this study; b Comparison of the expression levels of candidate reference genes in different treatments (Con: virus-free SBPHs infection; RSV: RSV SBPHs infection; RBSDV: RBSDV SBPHs infection). 14-day-old rice seedlings were inoculated with or without viruliferous nymphs (RSV and RBSDV) for 3 days. Total RNA was extracted from RSV- and RBSDV-infected seedlings, respectively. Values were given in the form of RT-qPCR quantification cycle numbers across all tested samples. Bars indicate standard error of the mean

Because the transcript levels of reference genes actually varied under different experimental conditions [18, 21], we analyzed if virus infection (RSV and RBSDV infection) altered the expression of any of the 14 candidate genes. The Ct values obtained for each gene were compared in RSV- and RBSDV-infected against virus-free treatment. In this study, our results showed that each candidate reference gene expressed constantly in different treatments (Fig. 3b). Variation in transcript quantity of RSV- and RBSDV-infected samples revealed that each tested gene approximately exhibited the similar transcript levels. Therefore, above 14 references genes were selected in the subsequent analyses.

Comparison of the expression stability of housekeeping genes by geNorm and NormFinder

The expression stability of abovementioned candidate reference genes could be assessed by different software programs. To date, the most popular and useful method is geNorm algorithm [22, 28]. The geNorm is a statistical algorithm which determines the gene stability measure (M) of all the genes under investigation, based on the geometric averaging of multiple control genes and mean pairwise variation of a gene from all other control genes in a given set of samples [28]. It relies on the principle that the expression ratio of two ideal internal control genes is identical in all the samples, regardless of the experimental condition and cell-type. Genes with the lowest M values have the most stable expression.

In our experimental conditions, geNorm analysis showed that GAPDH and UBQ 10 were ranked as the best reference genes in RSV-infected plants, followed by EXP, while UBC and 18S rRNA were the least stable genes (Fig. 4a). When considering RBSDV-infected treatments, the average expression stability value (M) of Actin1 and UBC were two most stable genes, followed by UBQ 10, and those of β-TUB and TIP41-like were the highest two (Fig. 4b).

Fig. 4

Average expression stability values (M) and pairwise variation (V) analyses of candidate reference genes under virus-infected conditions by geNorm. 14 candidate reference genes were amplified in cDNA samples from RSV- and RBSDV-infected seedlings. A lower M value indicated more stable expression. Mean expression stability following stepwise exclusion of the least stable gene across samples from RSV-infected treatments (a) or RBSDV-infected samples (b), respectively. The pairwise variation (V n/n+1) measured the effect of adding additional reference genes on the normalization factor (The dash line denotes 0.15 cut-off V value) for these treatments (c and d). Calculations were performed as described in “Materials and Methods” section

In general, it is not sufficient by using only one most stable reference gene to obtain an accurate and reliable result. Therefore, the next question is how many reference genes should be included for RT-qPCR normalization. The geNorm software also calculated the pairwise variations (V n/n+1 ) between two sequential normalization factors to determine the necessity of adding further reference gene(s). As suggested by Vandesompele et al. [28], 0.15 is a cutoff V value, below which the inclusion of an additional reference gene is not required. But this proposed value can not be taken as an absolute rule and might depend on the data. It is advisable to add additional reference genes to the normalization factor until the added gene has no significant effect [28, 30].

The pairwise variation analyses showed that all of V values were less than 0.15 in this set of samples (Fig. 4c and d). Moreover, the V 3/4 values were higher than that of V 2/3 and V 4/5. Although V 4/5 values were much lower than that of V 3/4 and even lower values were obtained by adding more reference genes, considering practical applications, two reference genes was optimal in our experimental conditions.

Different algorithm methods sometimes may produce different results from the same data set. Hence, all of the data were reassessed by NormFinder to avoid introducing unnecessary bias. NormFinder is a mathematical model which performs separate analysis of sample subgroups. It estimates intra- and inter-group variations and combines both results into a consistent value for each investigated gene [29]. Interestingly, we found that the ranking generated by this approach (Fig. 5) was approximately similar with those determined by geNorm (Fig. 4a and b). UBC and 18S rRNA were sill ranked higher than other reference genes in RSV-infected samples. In the RBSDV-infected samples, β-TUB and TIP41-Like were also ranking the highest. Although in RSV-infected samples, the rank order of EXP and GAPDH were slightly altered in NormFinder analysis in comparison with those in geNorm analyses, the stability of UBQ 10 + GAPDH was high enough for reliable normalization compared to those of UBC and 18S rRNA. Also, the similar result was obtained in the RBSDV-treated samples, showing that Actin1 + UBC were suitable for normalization in RT-qPCR analysis.

Fig. 5

The ranking of candidate reference genes based on stability values calculated by NormFinder under virus-infected conditions. 14 candidate reference genes were amplified in cDNA samples from RSV- and RBSDV-infected seedlings. Relative quantifications of RSV-infected (a) or RBSDV-infected (b) seedlings, were performed respectively, as described in “Materials and Methods” section

Comparison of single and multiple reference gene(s) in quantitative Real-time PCR normalization

In order to demonstrate the usefulness of the above validated reference genes in RT-qPCR, two important genes in virus resistant response, OsPR1b [3133] and OsWRKY [3437], were selected and applied as target genes (Fig. 6). In RSV- and RBSDV-infected samples, UBQ 10 + GAPDH and UBC + Actin1 were used as multiple reference genes, respectively. Meanwhile, TIP41-Like was used as single reference gene in our experimental conditions. As expected, in comparison with single reference gene, the normalization of OsPR1b using UBQ 10 + GAPDH and UBC + Actin1, respectively, resulted in significant increase in transcript levels under virus infection conditions (Fig. 6a and b). The result was similar with previous studies [38, 39]. Meanwhile, the expression levels of OsWRKY were progressively increased in the early stage of infection and then decreased (Fig. 6c and d). By contrast, when TIP41-Like was used as single reference gene, the increasing transcript patterns of OsPR1b and OsWRKY were differentially lower than those normalized by multiple reference genes.

Fig. 6

Relative expression levels of OsPR1b and OsWRKY using single or multiple reference gene(s) for normalization during RSV- (a and c) and RBSDV- (b and d) infection. 14-day-old rice seedlings were inoculated with viruliferous nymphs (RSV and RBSDV) for 3 days. Total RNA was extracted from RSV- and RBSDV-infected seedlings, respectively. TIP41-Like was used as a single reference gene, while UBQ 10 + GAPDH and UBC + Actin1 were used multiple reference genes under RSV- and RBSDV-infection plants, respectively, in our experimental conditions. Bars indicate standard error of the mean

We also detect the virus gene expression in plants during viral infection. According to the results, the expression patterns of the virus gene are consistent with the expression level of virus resistant gene (Fig. 7). The expression level normalized by multiple reference genes were better than those by single reference, therefore multiple reference genes were suitable for normalization in RT-qPCR analysis.

Fig. 7

Relative expression levels of virus genes (CP and SP for RSV, P5-1 and P9-1 for RBSDV) using single or multiple reference gene(s) for normalization during RSV- (a and c) and RBSDV- (b and d) infection. 14-day-old rice seedlings were inoculated with viruliferous nymphs (RSV and RBSDV) for 3 days. Total RNA was extracted from RSV- and RBSDV-infected seedlings, respectively. UBQ 10 + GAPDH and UBC + Actin1 were used multiple reference genes under RSV- and RBSDV-infection plants, respectively


RSV and RBSDV are two significant rice viruses that threat rice production. To understand the mechanism of viral transmission and discover some resistance genes, a reliable quantitation method is particularly needed. RT-qPCR is still the most commonly used technique. However, an accurate and reliable gene expression analyzed by RT-qPCR highly requires stably expressed reference genes for normalization. In fact, no one gene can act as a universal reference under different experimental conditions, and the normalization of gene expression with a single reference gene can usually lead to relatively errors [28, 30]. Thus, two or more stably reference genes for normalizing are essential.

In this study, the suitable reference genes for normalizing gene expression in RSV- and RBSDV-infected rice plants were identified. We first assessed the integrity of RNA samples (Fig. 2c). Meanwhile, the specificity of the RT-qPCR primer pairs was confirmed by agarose gel electrophoresis (Fig. 2b) and melting curves analysis (Fig. 2a). The relative abundance of candidate reference genes was further identified (Fig. 3). The results of these indicated that the candidate reference genes could be used for the subsequent analyses.

To date, the commonly used methods to assess the stability of reference genes are geNorm [28], NormFinder [29] and BestKeeper [40]. As BestKeeper cannot test more than 10 candidates, it was not used in this study. By using geNorm and NormFinder, algorithms, 14 reference genes (Table 1) were evaluated. Our findings revealed that UBQ 10 and GAPDH were overall the most stable genes in RSV-infected rice (Figs. 4a and 5a), and Actin1 and UBC ranked the best candidate genes under RBSDV infection (Figs. 4b and 5b). Interestingly, the suitable reference genes for RSV- and RBSDV-infected rice plants were different, it is indicated that the expression of the reference genes can vary under given situations [21]. Therefore, it is vital to choose suitable reference gene(s) for normalization of gene expression. Moreno et al. [41] showed that in different Cassava brown streak virus (CBSV)-infected tissues, PP2A, UBQ10 and GTPb appeared to be the most stable genes. In infected tomato plants, GAPDH and UBQ indicated as the most appropriate internal standards both in leaves and root tissues, and Actin was also stably expressed in the infected plants [42]. When infected with four commonly known tomato viral pathogens, ACT, CAC and EF1-α were considered as the most suitable reference genes in the studies of host–virus interactions [43]. Therefore, the selected multiple reference genes (UBQ 10 + GAPDH for RSV, and UBC + Actin1 for RBSDV) could be used for the normalization of gene expression pattern in RSV- and RBSDV-infected rice plants.

Normally, EF1-α and 18S rRNA were used as the reference genes in RT-qPCR experiments. Previous results showed that EF1-α was expressed stably in potato during biotic and abiotic stress [44], and 18S rRNA was identified to be suitable for normalisation in Barley yellow dwarf virus-infected cereals [45]. In our study, the performance of above two genes was dissatisfactory. These were similar to the earlier studies in Cicer arietinum and virus-infected tomato, showing the considerable alterations in the transcript levels of EF1-α and 18S rRNA [42, 46]. Thus, the two traditional reference genes might not be the optimal choices for quantitating transcript level in rice during RSV- and RBSDV-infection. Meanwhile, TUB was widely used as reference gene for gene expression analysis. However, due to the potential regulation in various physiological states, the suitability as internal control has been questioned [45, 47]. In our experimental conditions, α-TUB and β-TUB (in particularly) were not stable either, especially in RBSDV-infected plants. Taken together, our data demonstrated that α-TUB and β-TUB were unsuitable for the normalization of gene expression levels in viral infected rice.

To further investigate the suitability of the selected reference genes, the expression levels of two virus resistant response genes, the gene expression of OsPR1b and OsWRKY, was compared by using single and multiple reference gene(s) for normalization. In the earlier studies [38, 39], for example, the expression of OsPR1b gene was strongly induced by virus infection. Our results further showed that when using TIP41-Like gene for normalization, the induction patterns of OsPR1b and OsWRKY transcripts were not stronger than those using two reference genes (Fig. 6). Thus, the application of multiple reference genes is a better choice for gene expression analysis under virus infection situation.

In conclusion, we evaluated 14 candidate reference genes, and identified UBQ 10 + GAPDH and UBC + Actin1 as the most stably expressed reference genes in RSV- and RBSDV-infected rice, respectively. These genes will enable more accurate and reliable normalization of RT-qPCR analysis for gene expression studies to get insight into complex regulatory networks and will most probably lead to the identification of genes relevant to new biological processes, in RSV- and RBSDV-infected plants.


Plant materials, virus isolates and inoculation assay

Rice (Oryza sativa L. cv. Nipponbare) seeds were used throughout the experiments. After the disinfection with 0.1 % HgCl2 for 1 h and thorough washing with RO (reverse-osmosis) water, seeds were soaked overnight in RO water and incubation at 25 °C for 1 day. Seedlings were then grown on plastic chambers using Kimura B nutrient solution [48], with a 14/10 h (day/night) regimes at 28 ± 1 °C.

Rice plants infected with RSV and RBSDV, were collected from Jianhu county, Jiangsu Province in July 2014. Young instar nymphs of SBPHs were fed RSV- and RBSDV-infected rice plants for 2 days to acquire the virus, respectively. The virus was maintained by SBPHs in an insect-rearing room at a temperature of 25 °C. Virus-free SBPHs were also used for the control inoculation. Viruliferous or virus-free SBPHs were reared on rice seedlings (Oryza sativa L. cv. Wuyujing No. 3) in glass vessels at 22 °C under alternating photoperiods of 14 h of light and 10 h of dark. 14-day-old seedlings were inoculated with 10 viruliferous nymphs per plant and were kept in a growth chamber. After the incubation for 3 days, planthoppers were removed and plants were transferred to field. Virus-free SBPHs were also used for control inoculation. Leaf tissues were collected after 7, 14 and 21 days of virus treatment and immediately frozen in liquid nitrogen and stored at −80 °C until further analyses.

According to previous report [49], RT-PCR was carried out to detect RSV and RBSDV in rice plants. Primer RSRB-R (5′-CCYATCACAAASAAATMAAAAT-3′) paired with primer RSV-F (5′-AGATCCAGAGAGAGTCACGGAAG-3′) was used to amplify a specific 1114 bp fragment for detection of RSV. Primers RSRB-R and RBSDV-F (5′-GTTCAAAGACAATACACTCAAAA-3′) were used to amplify a 414 bp product, which was specific for RBSDV.

Preparation and quantification of DNA-free total RNA

Total RNA was extracted by Trizol reagent (Invitrogen, Gaithersburg, MD, USA) according to the manufacturer’s instructions. The RNA was dissolved in DNase-treated distilled water. The concentration and quality of isolated RNA were determined using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). Only the RNA samples with 260/280 ratio (an indication of protein contamination) between 1.9 and 2.1 and 260/230 ratio (an indication of reagent contamination) greater than 2.0, were used for the analyses [50, 51]. The integrity of RNA samples was assessed by agarose gel electrophoresis. All RNA samples were adjusted to the same concentration and measured again to homogenize RNA for the subsequent experiments.

Reverse transcription and quantitative Real-time PCR assay

cDNA was synthesized from 2 μg of total RNA using an oligo (dT) primer and M-MLV reverse transcriptase (BioTeke, Beijing, China). RT-qPCR was performed using the SsoFastTM Eva Green® Supermix (Bio-Rad) with the Bio-Rad iQ5 RT-qPCR system. Combined with all internal control genes used recently [25, 27, 28], we further selected 14 genes as the candidate reference genes in this study (Table 1). Additionally, the specific primers 5′-ACGCCTTCACGGTCCATAC-3′ and 5′-AAACAGAAAGAAACAGAGGGAGTAC-3′ were used for OsPR1b (AK107926); 5′-TCAGTGGAGAAGCGGGTGGTG-3′ and 5′-GGGTGGTTGTGCTCGAAGGAG-3′ were used for OsWRKY (EF143611). The efficiency and specificity of all the primers were checked by melting curve analysis, similar to previous report [51].

Using the RSV CP, SP gene and RBSDV P5-1, P9-1 gene, the dynamics of viral infection in rice plants was measured by reverse-transcription real-time PCR. The primer and products were shown in Table 2.

Table 2 Real-time PCR primers of virus genes (RSV and RBSDV) used in this study

Each quantitative Real-time PCR was performed using 0.5 μL cDNA, 10 μL of SsoFastTM Eva Green® Supermix and 0.2 μM forward and reverse primers were used in a total volume of 20 μL. All tubes were subjected to denaturation for 10 min at 95 °C, followed by 40 cycles of 95 °C for 10 s, and 56 °C for 20 s. SYBR Green absorbance was detected at 56 °C. All reactions were conducted in triplicate. Amplicon dissociation curves (melting curves), were recorded after cycle 40 by heating from 60 °C to 95 °C at a ramp speed of 1.9 °C ⁄min.

Data processing

The expression stability of the candidate reference genes were analyzed by using geNorm [28] and NormFinder [29]. Expression levels were assessed based on the number of amplification cycles needed to reach a specific threshold (cycle threshold; Ct) in the exponential phase of RT-qPCR. For both programs, raw Ct values of each gene were converted into relative quantities before inputting into software. The relative expression levels of corresponding genes were calculated relative to the maximum abundance in different samples. The highest relative expression for each gene was set to 1.0. The geNorm algorithm [28] derives a stability measure (M). Via a stepwise exclusion of the least stable gene, it creates a stability ranking. It also estimates the number of genes required to calculate a robust normalization factors, and performs a stepwise analysis (more stable to less stable genes) to calculate the pairwise variation (V n/n+1 ) between two sequential normalization factors containing an increasing number of genes. NormFinder algorithm [29] used an ANOVA-based model to estimate intra- and inter-group variation. It combines these results to provide a direct measurement of the variation in the expression for each gene [52]. Statistical significance of Ct differences between treatments was calculated by the Mann–Whitney t test using the GraphPad Prism 5 software.


  1. 1.

    Hibino H. Biology and epidemiology of rice viruses. Annu Rev Phytopathol. 1996;34:249–74.

  2. 2.

    Kuribayashi K. On the relationship between rice stripe disease and Delphacodes striatella Fallen. J Plant Prot. 1931;18:565–71. 636–640.

  3. 3.

    Shikata E, Kitagawa Y. Rice black-streaked dwarf virus: its properties, morphology and intracellular localization. Virology. 1977;77:826–42.

  4. 4.

    Falk BW, Tsai JH. Biology and molecular biology of viruses in the genus Tenuivirus. Annu Rev Phytopathol. 1998;36:139–63.

  5. 5.

    Kiso A, Yamamoto T. Infection and symptom in rice stripe disease with special reference to disease-specific protein other than virus. Rev Plant Prot Res. 1973;6:75–100.

  6. 6.

    Shinkai A. Studies on insect transmission of rice virus diseases in Japan. Bull Nat Inst Agric Sci Ser C. 1962;14:1–12.

  7. 7.

    Iida TT, Shinkai A. Transmission of dwarf, yellow dwarf, stripe and black-streaked dwarf. In: RF Chandler, editor. The virus diseases of the rice plant. International Rice Research Institute, the Johns Hopkins Press: Baltimore, MD; 1969. p. 125–129.

  8. 8.

    Sun DZ, Jiang L. Research on the inheritance and breeding of rice stripe resistance. Chin Agric Sci Bull. 2006;22:318–22.

  9. 9.

    Fang S, Yu J, Feng J, Han C, Li D, Liu Y. Identification of rice black-streaked dwarf fijivirus in maize with rough dwarf disease in China. Arch Virol. 2001;146:167–70.

  10. 10.

    Isogai M, Uyeda I, Lee BC. Detection and assignment of proteins encoded by rice black streaked dwarf fijivirus S7, S8, S9 and S10. J Gen Virol. 1998;79:1487–94.

  11. 11.

    Dong GK, Wang EG, Luo GL, Lin LW, Guan MP, Zhang ZD, et al. Occurrence of rice black-streaked dwarf disease in late season hybrid rice and its control strategy. Acta Agric Zhejiangensis. 1999;11:364–7.

  12. 12.

    Gachon C, Mingam A, Charrier B. Real-time PCR: what relevance to plant studies? J Exp Bot. 2004;55:1445–54.

  13. 13.

    Bustin SA, Benes V, Nolan T, Pfaffl MW. Quantitative real-time RT-PCR–a perspective. J Mol Endocrinol. 2005;34:597–601.

  14. 14.

    Dekkers BJW, Willems L, Bassel GW, van Bolderen-Veldkamp RP, Ligterink W, Hilhorst HWM, et al. Identification of reference genes for RT-qPCR expression analysis in Arabidopsis and tomato seeds. Plant Cell Physiol. 2012;53:28–37.

  15. 15.

    Bustin SA, Nolan T. Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004;15:155–66.

  16. 16.

    Huggett J, Dheda K, Bustin S, Zumla A. Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005;6:279–84.

  17. 17.

    Skern R, Frost P, Nilsen F. Relative transcript quantification by quantitative PCR: roughly right or precisely wrong? BMC Mol Biol. 2005;6:10.

  18. 18.

    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:5–17.

  19. 19.

    Lee PD, Sladek R, Greenwood CMT, Hudson TJ. Control genes and variability: absence of ubiquitous reference transcripts in diverse mammalian expression studies. Genome Res. 2002;12:292–7.

  20. 20.

    Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A. Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004;313:856–62.

  21. 21.

    Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, et al. Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999;75:291–5.

  22. 22.

    Gutierrez L, Mauriat M, Guénin S, Pelloux J, Lefebvre JF, Louvet R, et al. 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:609–18.

  23. 23.

    Seddas P, Boissinot S, Strub JM, Van Dorsselaer A, Van Regenmortel MH, Pattus F. Rack-1, GAPDH3, and actin: proteins of Myzus persicae potentially involved in the transcytosis of beet western yellows virus particles in the aphid. Virology. 2004;325:399–412.

  24. 24.

    Tamborindeguy C, Monsion B, Brault V, Hunnicutt L, Ju HJ, Nakabachi A, et al. genomic analysis of transcytosis in the pea aphid, Acyrthosiphon pisum, a mechanism involved in virus transmission. Insect Mol Biol. 2010;19:259–72.

  25. 25.

    Brenndörfer M, Boshart M. Selection of reference genes for mRNA quantification in Trypanosoma brucei. Mol Biochem Parasitol. 2010;172:52–5.

  26. 26.

    Lord JC, Hartzer K, Toutges M, Oppert B. Evaluation of quantitative PCR reference genes for gene expression studies in Tribolium castaneum after fungal challenge. J Microbiol Methods. 2010;80:219–21.

  27. 27.

    Radonic A, Thulke S, Bae HG, Muller M, Siegert W, Nitsche A. Reference gene selection for quantitative real-time PCR analysis in virus infected cells: SARS corona virus, Yellow fever virus, Human Herpesvirus-6, Camelpox virus and Cytomegalo virus infections. Virol J. 2006;2:7.

  28. 28.

    Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:7. resrarch0034–resrarch0034.11.

  29. 29.

    Andersen CL, Jensen JL, Ørntoft 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:5245–50.

  30. 30.

    Gutierrez L, Mauriat M, Pelloux J, Bellini C, Van Wuytswinkel O. Towards a systematic validation of references in real-time RT-PCR. Plant Cell. 2008;20:1734–5.

  31. 31.

    Hao ZN, Wang LP, He YP, Liang JG, Tao RX. Expression of defense genes and activities of antioxidant enzymes in rice resistance to rice stripe virus and small brown planthopper. Plant Physiol Biochem. 2011;49:744–51.

  32. 32.

    Satoh K, Kondoh H, De Leon T, Macalalad RJA, Cabunagan RC, Cabauatan PC, et al. Gene expression responses to Rice tungro spherical virus in susceptible and resistant near-isogenic rice plants. Virus Res. 2013;171:111–20.

  33. 33.

    Lee KJ, Kim K. The rice serine/threonine protein kinase OsPBL1 (ORYZA SATIVA ARABIDOPSIS PBS1- LIKE 1) is potentially involved in resistance to rice stripe disease. Plant Growth Regul. 2015, 1–9. doi: 10.1007/s10725-015-0036-z.

  34. 34.

    Qiu D, Xiao J, Ding X, Xiong M, Cai M, Cao Y, et al. OsWRKY13 mediates rice disease resistance by regulating defense-related genes in salicylate- and jasmonate-dependent signaling. Mol Plant Microbe Interact. 2007;20:492–9.

  35. 35.

    Bai FW, Yan J, Qu ZC, Zhang HW, Xu J, Ye MM, et al. Phylogenetic analysis reveals that a dwarfing disease on different cereal crops in China is due to rice black streaked dwarf virus (RBSDV). Virus Genes. 2002;25:201–6.

  36. 36.

    Catoni M, Miozzi L, Fiorilli V, Lanfranco L, Accotto GP. Comparative analysis of expression profiles in shoots and roots of tomato systemically infected by Tomato spotted wilt virus reveals organ-specific transcriptional responses. Mol Plant Microbe Interact. 2009;22:1504–13.

  37. 37.

    Tao Z, Liu H, Qiu D, Zhou Y, Li X, Xu C, et al. A pair of allelic WRKY genes play opposite roles in rice–bacteria interactions. Plant Physiol. 2009;151:936–48.

  38. 38.

    Satoh K, Kondoh H, Sasaya T, Shimizu T, Choi IR, Omura T, et al. Selective modification of rice (Oryza sativa) gene expression by rice stripe virus infection. J Gen Virol. 2010;91:294–305.

  39. 39.

    Zheng WJ, Ma L, Zhao JM, Li ZQ, Sun FY, Lu XU. Comparative transcriptome analysis of two rice varieties in response to rice stripe virus and small brown planthoppers during early interaction. PLoS One. 2013;8:e82126.

  40. 40.

    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:509–15.

  41. 41.

    Moreno I, Gruissem W, Vanderschuren H. Reference genes for reliable potyvirus quantitation in cassava and analysis of Cassava brown streak virus load in host varieties. J Virol Methods. 2011;177:49–54.

  42. 42.

    Mascia T, Santovito E, Gallitelli D, Cillo F. Evaluation of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in infected tomato plants. Mol Plant Pathol. 2010;11:805–16.

  43. 43.

    Wieczorek P, Wrzesińska B, Obrępalska-Stęplowska A. Assessment of reference gene stability influenced by extremely divergent disease symptoms in Solanum lycopersicum L. J Virol Methods. 2013;194:161–8.

  44. 44.

    Wan H, Zhao Z, Qian C, Sui Y, Malik AA, Chen J. Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber. Anal Biochem. 2010;399:257–61.

  45. 45.

    Jarošová J, Kundu J. Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR. BMC Plant Biol. 2010;10:146.

  46. 46.

    Castro P, Roman B, Rubio J, Die JV. Selection of reference genes for expression studies in Cicer arietinum L.: analysis of cyp81E3 gene expression against Ascochyta rabiei. Mol Breed. 2012;29:261–74.

  47. 47.

    Mafra V, Kubo KS, Alves-Ferreira M, Ribeiro-Alves M, Stuart RM, Boava LP, et al. Reference genes for accurate transcript normalization in citrus genotypes under different experimental conditions. PLoS One. 2012;7:e31263.

  48. 48.

    Ma JF, Goto S, Tamai K, Ichii M. Role of root hairs and lateral roots in silicon uptake by rice. Plant Physiol. 2011;127:1773–80.

  49. 49.

    Li S, Wang X, Xu JX, Ji YH, Zhou YJ. A simplified method for simultaneous detection of Rice stripe virus and Rice black-streaked dwarf virus in insect vector. J Virol Methods. 2015;211:32–5.

  50. 50.

    Nolan T, Hands RE, Bustin SA. Quantification of mRNA using real-time RT-PCR. Nat Protoc. 2006;1:1559–82.

  51. 51.

    Kumar K, Muthamilarasan M, Prasad M. Reference genes for quantitative real-time PCR analysis in the model plant foxtail millet (Setaria italica L.) subjected to abiotic stress conditions. Plant Cell Tiss Org. 2013;115:13–22.

  52. 52.

    Chi XY, Hu RB, Yang QL, Zhang XW, Pan LJ, Chen N, et al. Validation of reference genes for gene expression studies in peanut by quantitative real-time RT-PCR. Mol Genet Genomics. 2012;287:167–76.

Download references


This study was financially supported by grants from the National Key Basic Research and Development Program (973 Program) of China (2013CBA01403), the National Natural Science Foundation (31101412), the Special Fund for Agro-scientific Research in the Public Interest (201303021), the Jiangsu Agricultural Scientific Self-innovation Fund ([15]1053) and the Jiangsu Province Science and Technology Support Project (BE2013301).

Author information



Corresponding author

Correspondence to Tong Zhou.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: TZ, YZ, WS. Performed the experiments: PF, RL, FS, YL, LD. Analysed the data: PF. Wrote the paper: PF, WS, TZ, YZ. All authors read and approved the final manuscript.

Additional file

Additional file 1: Figure S1.

Relative expression levels of OsPR1b and OsWRKY using single or multiple reference gene(s) for normalization during RSV- (A and C) and RBSDV- (B and D) infection. 14-day-old rice seedlings were inoculated with viruliferous nymphs (RSV and RBSDV) for 3 days. Total RNA was extracted from RSV- and RBSDV-infected seedlings, respectively. In our experimental conditions, UBQ 10 + GAPDH and UBC + Actin1 were used multiple reference genes under RSV- and RBSDV-infection plants. (DOC 303 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fang, P., Lu, R., Sun, F. et al. Assessment of reference gene stability in Rice stripe virus and Rice black streaked dwarf virus infection rice by quantitative Real-time PCR. Virol J 12, 175 (2015).

Download citation


  • Gene expression
  • RT-qPCR
  • Reference genes
  • RSV- and RBSDV-infected rice