Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2014 May 8;9(5):e97038.
doi: 10.1371/journal.pone.0097038. eCollection 2014.

Sequencing and validation of reference genes to analyze endogenous gene expression and quantify yellow dwarf viruses using RT-qPCR in viruliferous Rhopalosiphum padi

Keke Wu  1 , Wenwen Liu  2 , Thithi Mar  2 , Yan Liu  2 , Yunfeng Wu  3 , Xifeng Wang  2
Affiliations

Sequencing and validation of reference genes to analyze endogenous gene expression and quantify yellow dwarf viruses using RT-qPCR in viruliferous Rhopalosiphum padi

Keke Wu et al. PLoS One. .

Abstract

The bird cherry-oat aphid (Rhopalosiphum padi), an important pest of cereal crops, not only directly sucks sap from plants, but also transmits a number of plant viruses, collectively the yellow dwarf viruses (YDVs). For quantifying changes in gene expression in vector aphids, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a touchstone method, but the selection and validation of housekeeping genes (HKGs) as reference genes to normalize the expression level of endogenous genes of the vector and for exogenous genes of the virus in the aphids is critical to obtaining valid results. Such an assessment has not been done, however, for R. padi and YDVs. Here, we tested three algorithms (GeNorm, NormFinder and BestKeeper) to assess the suitability of candidate reference genes (EF-1α, ACT1, GAPDH, 18S rRNA) in 6 combinations of YDV and vector aphid morph. EF-1α and ACT1 together or in combination with GAPDH or with GAPDH and 18S rRNA could confidently be used to normalize virus titre and expression levels of endogenous genes in winged or wingless R. padi infected with Barley yellow dwarf virus isolates (BYDV)-PAV and BYDV-GAV. The use of only one reference gene, whether the most stably expressed (EF-1α) or the least stably expressed (18S rRNA), was not adequate for obtaining valid relative expression data from the RT-qPCR. Because of discrepancies among values for changes in relative expression obtained using 3 regions of the same gene, different regions of an endogenous aphid gene, including each terminus and the middle, should be analyzed at the same time with RT-qPCR. Our results highlight the necessity of choosing the best reference genes to obtain valid experimental data and provide several HKGs for relative quantification of virus titre in YDV-viruliferous aphids.

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Box and whisker plots of Cq values for the 4 candidate reference genes in 6 experimental groups (one of three viruses in either winged or wingless adults of Rhopalosiphum padi).
Each box shows the lower 25th and upper 75th percentiles with median Cq values; the whiskers mark the lower 5th and upper 95th percentiles of the Cq values in each data set. Experimental groups from left to right: WYDV-GPV in winged adult, BYDV-PAV in winged adult, BYDV-GAV in winged adult, WYDV-GPV and wingless adult, BYDV-PAV and wingless adult, and BYDV-GAV and wingless adult.
Figure 2
Figure 2. Determination of optimal reference genes for normalization by GeNorm.
The optimal number of reference genes was determined by pairwise variation V (A) between two sequential normalization factors containing an increasing number of reference genes (n and n +1 reference genes). The pairwise variation begins with the first two and three (V2/3) most stable candidate reference genes (B), terminating with the addition of the least stable one, here V3/4. The cut-off V value is set at 0.15 by the program GeNorm, below which an additional reference gene does not need to be included for calculating the normalization factor.
Figure 3
Figure 3. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-GAV-viruliferous winged adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 4
Figure 4. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-GAV-viruliferous wingless adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 5
Figure 5. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-PAV-viruliferous winged adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 6
Figure 6. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-PAV-viruliferous wingless adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 7
Figure 7. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-GPV-viruliferous winged adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 8
Figure 8. Mean relative expression values (± SE, n = 3) for endogenous aphid genes ago-1a-1 (A), ago-1a-2 (B), ago-1a-3 (C) and dcr1 (D) in BYDV-GPV-viruliferous wingless adults of Rhopalosiphum padi after different virus-feeding durations.
A–D, Expression values for each gene or each region of one gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 9
Figure 9. Mean relative titre of YDVs (± SE, n = 3) in viruliferous winged adults of Rhopalosiphum padi after different virus-feeding durations.
Virus titre is illustrated by relative expression values of the CP (A, C and E) and RTD (B, D and F) gene of YDVs. Expression values for the CP or RTD gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. A and B: WYDV-GPV; C and D: BYDV-PAV; E and F: BYDV-GAV. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
Figure 10
Figure 10. Mean relative titre of YDVs (± SE, n = 3) in viruliferous wingless adults of Rhopalosiphum padi after different virus-feeding durations.
Virus titre is illustrated by relative expression values of the CP (A, C and E) and RTD (B, D and F) gene of YDVs. Expression values for the CP or RTD gene at each duration were normalized with the reference gene(s) selected by GeNorm and then compared with a one-way ANOVA (p) among these durations with each normalization condition. A and B: WYDV-GPV; C and D: BYDV-PAV; E and F: BYDV-GAV. 2re = 2 best reference genes; 3re = 3 best reference genes; 4re = all 4 reference genes; bad2re = 2 least stable reference genes; EF-1α = EF-1α as the reference gene; 18S = 18S rRNA as the reference gene.
See this image and copyright information in PMC

References

    1. Hesler LS (2005) Resistance to Rhopalosiphum padi L. (Homoptera: Aphididae) in three triticale accessions. J Econ Entomol 98: 603–610. - PubMed
    1. Descamps LR, Sanchez C (2011) Population growth of Rhopalosiphum padi L. (Homoptera: Aphididae) on different cereal crops from the semiarid pampas of Argentina under laboratory conditions. Chilean Journal of Agricultural Research 71: 390–394.
    1. Hadi BAR, Flanders KL, Bowen KL, Murphy JF, Blount AR (2012) Survey of Barley Yellow Dwarf Virus and Cereal Yellow Dwarf Virus on Three Perennial Pasture Grasses in Florida. Journal of Entomological Science 47: 35–43.
    1. Krueger EN, Beckett RJ, Gray SM, Miller WA (2013) The complete nucleotide sequence of the genome of Barley yellow dwarf virus-RMV reveals it to be a new Polerovirus distantly related to other yellow dwarf viruses. Front Microbiol 4: 205. - PMC - PubMed
    1. Zhang W, Cheng Z, Xu L, Wu M, Waterhouse P, et al. (2009) The complete nucleotide sequence of the barley yellow dwarf GPV isolate from China shows that it is a new member of the genus Polerovirus. Arch Virol 154: 1125–1128. - PubMed

Publication types

LinkOut - more resources