Next Article in Journal
Willow (Salix babylonica) Extracts Can Act as Biostimulants for Enhancing Salinity Tolerance of Maize Grown in Soilless Culture
Next Article in Special Issue
Comparative Genome Analysis of Old World and New World TYLCV Reveals a Biasness toward Highly Variable Amino Acids in Coat Protein
Previous Article in Journal
Applicability of Point- and Polygon-Based Vegetation Monitoring Data to Identify Soil, Hydrological and Climatic Driving Forces of Biological Invasions—A Case Study of Ailanthus altissima, Elaeagnus angustifolia and Robinia pseudoacacia
Previous Article in Special Issue
Melatonin as a Possible Natural Anti-Viral Compound in Plant Biocontrol
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 4
10.3390/plants12040857
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera

by
Ge Zhang
1,2,†,
Zhuo Zhang
3,†,
Qionglian Wan
1,2,
Huijie Zhou
1,2,
Mengting Jiao
1,2,
Hongying Zheng
1,2,
Yuwen Lu
1,2,
Shaofei Rao
1,2,
Guanwei Wu
1,2,
Jianping Chen
1,2,
Fei Yan
1,2,
Jiejun Peng
1,2,* and
Jian Wu
1,2,*
1
State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agroproducts, Institute of Plant Virology, Ningbo University, Ningbo 315211, China
2
Key Laboratory of Biotechnology in Plant Protection of MARA and Zhejiang Province, Institute of Plant Virology, Ningbo University, Ningbo 315211, China
3
Hunan Plant Protection Institute, Hunan Academy of Agricultural Sciences, Changsha 410125, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2023, 12(4), 857; https://doi.org/10.3390/plants12040857
Submission received: 17 January 2023 / Revised: 11 February 2023 / Accepted: 13 February 2023 / Published: 14 February 2023
(This article belongs to the Special Issue Plant Virus Disease Control)
Browse Figures
Review Reports Versions Notes

Abstract

:
Quantitative real-time PCR (RT-qPCR) is a widely used method for studying alterations in gene expression upon infections caused by diverse pathogens such as viruses. Positive-sense single-stranded (ss(+)) RNA viruses form a major part of all known plant viruses, and some of them are damaging pathogens of agriculturally important crops. Analysis of gene expression following infection by ss(+) RNA viruses is crucial for the identification of potential anti-viral factors. However, viral infections are known to globally affect gene expression and therefore selection and validation of reference genes for RT-qPCR is particularly important. In this study, the expression of commonly used reference genes for RT-qPCR was studied in Nicotiana benthamiana following single infection by 11 ss(+) RNA viruses, including five tobamoviruses, four potyviruses, one potexvirus and one polerovirus. Stability of gene expression was analyzed in parallel by four commonly used algorithms: geNorm, NormFinder, BestKeeper, and Delta CT, and RefFinder was finally used to summarize all the data. The most stably expressed reference genes differed significantly among the viruses, even when those viruses were from the same genus. Our study highlights the importance of the selection and validation of reference genes upon different viral infections.

1. Introduction

Infections of positive-sense single-stranded RNA (ss(+) RNA) viruses are major threats to agricultural production [1]. To systemically infect plant cells, ss(+) RNA viruses cooperate with host factors to manipulate host genes, thereby establishing a suitable microenvironment for viral replication and trafficking. Therefore, analysis of gene expression upon viral infection is critical for the identification of potential anti-viral factors. To study virus–host interactions, fluorescence-based quantitative real-time PCR (RT-qPCR) is often used to quantify the transcript levels of target genes [2,3]. The use of qPCR to accurately quantify transcript abundance is largely influenced by the amount of initiating material, RNA quality, reverse transcriptase efficiency, amplification efficiency, etc. [4]. To normalize these differences, the use of stably expressed reference genes is required [5]. The internal reference gene is an internal reaction control whose sequence is different from the target gene. To serve as a reliable reference for RT-qPCR a gene must meet several important criteria, of which the most important is that its expression level is not affected by any of the experimental treatments [6]. In addition, its expression should vary minimally in different tissues as well as under different physiological conditions. Actin, GAPDH and EF1α play critical roles in life activities, and their expression is largely independent of external factors such as abiotic stresses [7]. Therefore, these genes are widely used as references in RT-qPCR to analyze gene expression [8]. In Nicotiana benthamiana (N. benthamiana), the stability of several reference genes has been tested in developmentally distinct tissues and during abiotic stresses [9]. Moreover, previous studies have also analyzed the stability of certain reference genes during infection with a variety of viruses, including necroviruses, benyviruses, hordeiviruses, potexviruses, and tobacco rattle virus [10,11]. However, the stability of expression of most reference genes in host plants during infection with most of the reported plant viruses has not been well studied.
It is known that the expression of some genes that can be used as references under certain conditions may show significantly altered expression under other conditions [7]. To systemically infect hosts, viruses replicate and assemble in host cells using host components to produce new progeny viruses and initiate a new round of life cycle [12]. In this process, plant viruses typically interfere with and engage the host’s cellular pathways, and a large portion of the most widely used reference genes are components of these cellular pathways [13]. In addition, different viral infestations may affect different pathways, and the extent to which these pathways are affected varies depending on the time of viral infestation and the type of cells infected. Some ss(+) RNA viruses may hijack host metabolic enzymes and housekeeping proteins to facilitate their replication. For instance, tomato bushy stunt virus (TBSV) can bind GAPDH, a commonly used reference gene that plays critical roles in glycolysis [14], membrane fusion, vascular bundles, nuclear RNA output, DNA replication, DNA repair [15] and RNA stability [16]. By binding GAPDH to the 3′ untranslated region (UTR) of the TBSV negative-sense strand during replication the synthesis of the positive-sense chain is stimulated by increasing the stability of template [17,18]. However, the stability of GAPDH expression during TBSV infection has not been carefully analyzed in many host plants. Another example is that tobacco mosaic virus (TMV) can bind eEF1A, a highly conserved and abundant housekeeping protein in eukaryotes [19], to maintain the activity of viral RdRp [20]. Moreover, TMV is known to induce the expression of EF1A [21]. Therefore, it is particularly important to select and validate RT-qPCR reference genes to analyze the expression of host genes upon viral infection [22].
Owing to its ability to be infected by a wide variety of plant viruses, N. benthamiana has been widely used as an experimental model for studying plant–virus interactions [23]. Plant virus vector-based expression systems, the construction of virus-induced gene silencing (VIGS) systems, and agrobacterium infiltration are all developed using N. benthamiana plants [24]. In this study, N. benthamiana infection was therefore used as a model to analyze the stability of commonly used reference genes in the context of infections caused by 11 ss(+) RNA viruses, including five tobamoviruses (TMV, tomato brown rugose fruit virus (ToBRFV), pepper mild mottle virus (PMMoV), tobacco mild green mosaic virus (TMGMV) and tomato mottle mosaic virus (ToMMV)), four potyviruses (turnip mosaic virus (TuMV), chilli ringspot virus (ChiRSV), chilli veinal mottle virus (ChiVMV) and pepper mottle virus (PePMoV)), one potexvirus (potato virus X (PVX)) and one polerovirus (pod pepper vein yellows virus (PoPeVYV)). Five commonly used algorithms, geNorm, NormFinder, BestKeeper, Delta CT, and RefFinder were adopted to identify the most stably expressed reference genes following infection by the different viruses. The results show that the most promising reference genes can even be different for viruses of the same genus and identify useful reference genes for future studies of gene expression following infection by these 11 viruses. Our analysis also revealed that the selection and validation of reference genes in the context of different viral infections are necessary.

2. Results

2.1. Assessment of Primer Specificity and Amplification Efficiency

A total of 28 reference genes and their forward and reverse primers used for RT-qPCR were obtained from published literature. The primer sequences of all genes and information on the length of amplification products are provided in Table S1. We used RT-qPCR melting curve analysis, 1% agarose gel electrophoresis and ethidium bromide staining, and PCR product sequencing to verify the specificity of the primers. Among all the 28 selected candidate genes in Table S1, genes in black letters indicate that the melting curve obtained by RT-qPCR was single-peaked and only one band was visualized after ethidium bromide staining (Figure 1A and Figure S1), indicating high specificity of the amplification product; genes in red letters indicate the presence of double peaks in the melting curve or a red curve (Figure 1B), revealing the presence of non-specific amplification for that primer or low transcript levels of the target genes; genes in blue letters indicate that, as revealed by Sanger sequencing of the RT-qPCR products after cloning into pGEM®-T vector, the obtained sequence of at least three clones did not match the expected target gene sequences. We finally selected 13 genes from the 28 candidates for the subsequent experiments. The amplification efficiencies of these 13 primer pairs were analyzed. The results showed that the amplification efficiencies ranged from 91.9% to 110% (Table 1), which were all within the acceptable range of 90–110%. Furthermore, the standard curves demonstrated good linear relationships (R2 > 0.9874) between the cycle threshold (Ct) values and the log-transformed copy numbers of 13 tested reference genes (Table 1). The above results indicate that the primers of these 13 pairs of candidate reference genes were suitable for following RT-qPCR experiments.

2.2. Confirmation of Viral Infection and Sampling

At 6 days after inoculation, plants infected with six viruses (TMV, ToBRFV, PMMoV, TMGMV, ToMMV, and TuMV) had obvious symptoms (curved leaves, yellow spotting, etc.) in the systemic leaves, whereas it took 15 days for plants infected with the other five viruses to show typical symptoms. Symptoms induced by infections of representative viruses from four genera are shown in Figure S2A. The data of RT-PCR confirmation of viral infection are shown in Figure S2B.

2.3. Expression Analysis of the 13 Selected Reference Genes by RT-qPCR

To analyze the differential expression of the 13 candidate reference genes in N. benthamiana under infection of each of the 11 viruses, RNA samples were prepared using leaf samples collected at either 6 or 15 dpi depending on the time symptoms appeared (see above). RNA integrity was analyzed by 5% urea polyacrylamide gel electrophoresis. After ethidium bromide staining, ribosomal RNA bands with sharp edges were observed, indicating high RNA integrity (Figure S3). RNA concentrations were determined and only the ones with a 260/280 ratio close to 2.0 and a 260/230 ratio within the range of 2.0–2.2 were used in RT-qPCR. Prior to reverse transcription, genomic DNA was digested by incubating the RNA samples with DNase I for 2 h at 37 °C. Melting curves of the RT-qPCR products were checked to make sure only a single sharp peak was obtained. Ct values of three biological replicates of each virus-infected group and the Mock infection group were displayed by box line plots to represent the variations (Figure 2).
Among the tobamoviruses, the greatest variation in gene expression was in GAPDH after infection by TMV, TMGMV, ToBRFV and ToMMV. Following PMMoV infection, Ct values of gene AGO2 varied the most. Gene Lip showed the smallest variations in Ct values after infection by TMV or ToBRFV; Ct values of gene Tspan varied the least in the case of TMGMV and ToMMV single infection; for PMMoV infection, the gene with the least variation in Ct value was KLC.
Among the four potyviruses TuMV, ChiRSV, ChiVMV, and PePMoV, the genes that varied the most in Ct values were GAPDH, AGO2, AGO2, and GADPH, respectively; the genes with the least variations in Ct values were RdR6, UBC, Lip, and Lip, respectively.
In the case of PVX infection (genus Potexvirus), the Ct value of gene Lip varied the most and the gene Tspan varies the least and for PoPeVYV (genus Polerovirus) the gene with the biggest variation in Ct value was L23 and the gene with the smallest variation was UBC.

2.4. Stability Analysis of the 13 Selected Reference Genes

To further evaluate the expression stability of the 13 candidate reference genes, four different algorithms (geNorm, NormFinder, BestKeeper, and Delta-CT) were used to calculate expression stabilities individually. Finally, a comprehensive analysis platform RefFinder was used to integrate the values obtained from the four algorithms to identify the one with the highest stability upon the infection of each virus.
As shown in Table 2, the different algorithms were not consistent with one another in identifying the most and least stable genes for each virus. When integrated in RefFinder (Table 3), the top four reference genes for TMV were Lip, eIF4A, EF1α, and L23 while GAPDH, AGO2, and ACT were the least stable; since the M value of GAPDH was greater than 1.5 it is very unsuitable as a reference in this regard. For PMMoV infection, RefFinder identified the top four stable reference genes as RdR6, UBC, Tspan, and ACT while the least stable were AGO2, PGK, and L23, none of which were suitable as a reference gene. For TMGMV, F-BOX, RdR6, Lip, and ACT were the top four stable reference genes while GAPDH, AGO2, UBC, and eIF4A were not suitable as reference genes. For ToBRFV, the top four stable reference genes were F-BOX, UBC, EF1α, and PGK, while AGO2, GAPDH, and L23 were the three least stable genes. For ToMMV, RdR6, Tspan, Lip, and eIF4A were the top four stable genes while the least stable were GAPDH, ACT, AGO2, and L23.
Among the potyviruses, RefFinder identified PGK, RdR6, eIF4A, and Tspan as the top four stable reference genes for TuMV, while the most unstable were GAPDH, AGO2, ACT, and UBC (Table 3). For ChiRSV, Tspan, UBC, EF1α, and RdR6 were the top four stable reference genes, while the most unstable were AGO2, L23, ACT, and GAPDH. For ChiVMV, Tspan, eIF4A, UBC, and F-BOX were the four most stable genes while the most unstable were AGO2, GAPDH, L23, and RdR6). For PePMoV, the four algorithms gave very similar results with Lip, eIF4A, KLC, and Tspan being the most stable (Table 2); the most unstable were AGO2, GAPDH, ACT, and RdR6 (Table 3).
For PVX, RefFinder identified F-BOX, Tspan, ACT, and PGK as the top four stable candidate reference genes, while the most unstable were Lip, KLC, EF1α, and RdR6. For PoPeVYV, the four most stable genes were Tspan, F-BOX, eIF4A, and PGK, with RdR6, AGO2, L23, and Lip being the most unstable.

3. Discussion

RT-qPCR has become one of the important methods for rapid, sensitive, and quantitative comparison of gene expression levels [25]. Many experimental factors, such as the quantity of starting material, the quality of RNA, and the efficiency of reverse transcriptase, can affect the efficacy of RT-qPCR [4,26]. Therefore, accurate normalization, which relies on the selection of housekeeping reference genes, is crucial. Unfortunately, there are no universal references that can be used under diverse conditions to accurately reflect alterations in gene expression [27]. In the case of different viral infections that can globally affect gene expression, it is particularly important to select and validate suitable internal references to identify the targets of viruses and potential anti-viral factors [28]. Therefore, a systematic evaluation of reference genes must be performed before the analysis of gene expression upon viral infections [29].
Interestingly, our comprehensive data show that different reference genes should be used for different viruses, even when those viruses are members of the same genus. For example, while Lip, eIF4A, EF1α, and L23 appeared to be most suitable as reference genes for infection by TMV, those most suitable for PMMoV (another tobamovirus) were RdR6, UBC, Tspan, and ACT. Thus, there is no single gene appropriate for use with all tobamoviruses. Similar conclusions were obtained after the analysis of gene expression in the context of infections by four different potyviruses. Our study therefore suggests that reference genes need to be evaluated and selected on a case-by-case basis in the context of viral infections; they cannot be simply adopted from previous studies using a healthy plant or a different virus.
Our conclusions are also supported by previous studies. F-BOX is stably expressed in a variety of plant species under different biotic or abiotic stresses. In Arabidopsis thaliana, F-BOX is one of the genes that exhibits the most stable expression levels when plants are infected by viruses [30] or when they are under metal stress [31]. Analysis of gene expression in N. benthamiana plants infected with five RNA plant viruses also revealed that the F-BOX was the most stable reference gene [10]. In our study, F-BOX was also the most suitable reference gene in N. benthamiana plants infected with TMGMV, ToBRFV, and PVX. However, F-BOX expression levels varied a lot in plants infected by TMV or TuMV. Tetraspanins are a family of small transmembrane proteins that play crucial roles in intracellular transport, cellular signal transduction, and cell migration [32]. More than 33 Tspans, such as Tspan7, span9, CD9, CD81, and CD151, are known to be involved in the process of virus infestation in human cells [33]. Tetraspanins in A. thaliana facilitate CMV infection through specific membrane-recognition motifs [34]. However, in our study Tspan (Tetraspanin-20) displayed stable expression in N. benthamiana plants infected by ToMMV, ChiRSV, ChiVMV or PVX, suggesting that Tspan plays no significant role in infection by these viruses. GAPDH, AGO2, and actins are three of the most widely used reference genes in various plant species [22,35], but they have altered expression under various biotic or abiotic stresses [26,30,36]. Consistently, our study also revealed the altered expression of GAPDH, AGO2, and actins upon viral infections, suggesting that these are not suitable as reference genes in studies of viral infection. Interestingly, a previous study reported that L25 and EF-1α can be used as stable reference genes for gene expression analysis in both developmental and stress-treated samples of N. benthamiana plants [9]. The complexity of gene expression analysis under biotic stress treatments, especially viral infections, was highlighted by that study and our results.

4. Materials and Methods

4.1. Plant Materials and Virus Inoculation

N. benthamiana plants were grown in a growth chamber with a 16-h light/8-h dark cycle at 24 °C and 50% relative humidity. Infectious clones of 11 viruses, including TMV, PMMoV, TMGMV, ToBRFV, ToMMV, TuMV, ChiRSV, ChiVMV, PePMoV, PVX, and PoPeVYV were all from our laboratory. The information of genome size, backbone vector, and promoter were listed in Table S2. All viruses were agro-inoculated into N. benthamiana at the five-true-leaf stage according to published methods [37]. Three biological replicates (three plants per replicate) were included for each virus and mock control. Viral infection was confirmed by disease symptoms.

4.2. Total RNA Isolation and cDNA Synthesis

Depending on the onset of symptoms (see Results Section 2.2, above), leaf samples were collected from the first two fully expanded systemic leaves of N. benthamiana plants infected with TMV, ToBRFV, PMMoV, TMGMV, ToMMV, and TuMV at 6 days post inoculation (dpi) and from those infected with the remaining five viruses and 15 dpi. Total RNA was extracted using TRIzolTM Reagent (Invitrogen, Shanghai, China) according to the manufacturer’s instructions. RNA concentration was measured using a Nanodrop ONE spectrophotometer (Nanodrop Technologies, Rockland, DE, USA). RNA integrity was analyzed by 5% urea polyacrylamide gel electrophoresis, according to standard procedures [38]. First Strand cDNA was synthesized by reverse transcription (RT) using EasyScript®. All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (Transgen, Beijing, China) with 1 μg total RNA as template. All procedures were completed following the manufacturer’s protocol. Thermal conditions for RT were: 15 min at 37 °C, and then 5 s at 85 °C to inactivate the enzyme.

4.3. Selection of Candidate Reference Genes

Reference genes and the forward and reverse primers were obtained from published literature [10,39,40,41,42,43,44], or designed in this study. Primer sequences and the references are listed in Table S1. Primer specificity was evaluated using RT-qPCR and PCR amplification, in which the cDNA sample was derived from an RNA sample from a healthy N. benthamiana plant. The amplified PCR products were visualized by ethidium bromide staining after 1% agarose gel electrophoresis. PCR products were also purified and cloned into blunt zero T-vector (Transgen, Beijing, China) and sequenced. The amplification efficiency of each pair of primers was analyzed using the built-in software of the Roche LightCycler480 instrument and EXCEL. Standard curves were generated from a five-fold dilution (initial concentration = 14.6 ng/μL) series of one sample.

4.4. RT-PCR Confirmation of Viral Infection and RT-qPCR

RT-PCR was performed to confirm the existence of viral genome RNA in systemic leaves. Random primers were used to synthesize cDNA through RT. Primer sequences of PCR were listed in Table S3. RT-qPCR was performed on a LightCycler480 instrument (Roche, Switzerland) with 384-well PCR plates. Reaction mixtures were composed of 18 μL of ChamQ Universal YBR QPCR Master Mix (Vazyme, Nanjing, China), 0.6 μL of each primer (10 μM), 6 μL of cDNA template, and 10.8 μL of RNase free water. Three biological replicates were included for each experiment, and three technical replicates were set for each biological replicate. Thermal cycles of RT-qPCR were: 95 °C for 30 s, then 40 cycles of 15 s at 95 °C, 15 s at 55 °C, and 30 s at 72 °C. Melting curves were prepared through the following thermal conditions: 95 °C for 10 s, 60 °C for 30 s, and 95 °C for 15 s. Amplification efficiencies were calculated using the built-in software of the Roche LightCycler 480 instrument.

4.5. Evaluation of the Expression Stability of Reference Genes Using GeNorm, NormFinder, Delta-CT, BestKeeper, and RefFinder

To analyze gene expression stability, four algorithms (geNorm [45], NormFinder [46], BestKeeper [47], and Delta-CT [48]) and an online comprehensive analysis platform RefFinder [49] were used following manufacturer’s instructions.
For geNorm and NormFinder, the raw Ct values were converted to the data input format (2−ΔCt) required by the software. Moreover, ΔCt = each Ct value—the lowest Ct value. In contrast, raw Ct values were directly used as input data in BestKeeper and RefFinder analysis. The output of geNorm software analysis is the expression stability value M, which represents the average pairwise variation between a single gene and other internal reference genes. With M = 1.5 as the cutoff value, a low M value indicates stable gene expression. In addition, geNorm also provides pairwise variation values (Vn/n+1) to represent the optimal number of reference genes for normalization. When Vn/n+1 < 0.15, the most suitable number of reference genes was n. When Vn/n+1 > 0.15, the most appropriate number of reference genes is n + 1. Normfinder provides the stability value SV to evaluate the differences in the expression of candidate reference genes and the variations among sample groups. Like the M value, a low SV value indicates stable gene expression. The Delta-CT method identifies the potential reference genes by comparing the relative expression between gene pairs. According to the lowest coefficient of variation (CV) and the lowest relative standard deviation (SD), BestKeeper software was used to determine the most stable reference genes. When SD > 1, the expression of the reference gene was unstable. RefFinder is a web-based comprehensive reference gene analysis platform. It integrates the values obtained using algorithms geNorm, NormFinder, BestKeeper, and Delta-CT to obtain the optimal internal reference genes.

5. Conclusions

In conclusion, our study identified optimal RT-qPCR reference genes for gene expression analysis in N. benthamiana plants upon single infections by 11 ss(+) RNA viruses from four genera. Our data also highlighted the importance of the selection and validation of reference genes for gene expression analysis in the context of viral infections.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12040857/s1, Figure S1: Analysis of RT-qPCR product of 28 genes through 1% agarose gel electrophoresis and ethidium bromide staining; Figure S2: Symptoms induced by infections of 11 viruses from four genera (A) and RT-PCR confirmation of viral infection (B); Figure S3: Analysis of RNA integrity through 5% urea-PAGE gel electrophoresis and ethidium bromide staining; Table S1: Information of the 28 selected potential RT-qPCR reference genes and primer sequences; Table S2: Information of expression vectors of the 11 viruses used in this study; Table S3: Sequences of primers used in RT-PCR to confirm viral infection.

Author Contributions

J.P. and J.W. generated the idea; G.Z. and Z.Z. performed experiments, G.Z., Z.Z, Q.W., H.Z. (Huijie Zhou), M.J., H.Z. (Hongying Zheng), Y.L., S.R., G.W., J.C. and F.Y. carried out the analysis. G.Z. and Z.Z prepared the first draft of the manuscript. J.P. and J.W. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China awarded to J.W. [Grant No. 32272483] and Natural Science Foundation of Zhejiang Province awarded to J.P. [Grant No. LY22C140001].

Data Availability Statement

The data presented in this article are available on request from the corresponding authors.

Acknowledgments

We are grateful to Mike Adams for careful reading and revising the manuscript.

Conflicts of Interest

There authors have no conflict of interest.

References

  1. Chauhan, P.; Singla, K.; Rajbhar, M.; Singh, A.; Das, N.; Kumar, K. A systematic review of conventional and advanced approaches for the control of plant viruses. J. Appl. Biol. Biotechnol. 2019, 7, 89–98. [Google Scholar]
  2. Bustin, S.A. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J. Mol. Endocrinol. 2000, 25, 169–193. [Google Scholar] [CrossRef]
  3. Wang, T.; Brown, M.J. mRNA Quantification by Real Time TaqMan Polymerase Chain Reaction: Validation and Comparison with RNase Protection. Anal. Biochem. 1999, 269, 198–201. [Google Scholar] [CrossRef]
  4. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef]
  5. Mallona, I.; Lischewski, S.; Weiss, J.; Hause, B.; Egea-Cortines, M. Validation of reference genes for quantitative real-time PCR during leaf and flower development in Petunia hybrida. BMC Plant Biol. 2010, 10, 4. [Google Scholar] [CrossRef]
  6. Chervoneva, I.; Li, Y.; Schulz, S.; Croker, S.; Wilson, C.; Waldman, S.A.; Hyslop, T. Selection of optimal reference genes for normalization in quantitative RT-PCR. BMC Bioinform. 2010, 11, 253. [Google Scholar] [CrossRef] [PubMed]
  7. Thellin, O.; Zorzi, W.; Lakaye, B.; De Borman, B.; Coumans, B.; Hennen, G.; Grisar, T.; Igout, A.; Heinen, E. Housekeeping genes as internal standards: Use and limits. J. Biotechnol. 1999, 75, 291–295. [Google Scholar] [CrossRef] [PubMed]
  8. Eisenberg, E.; Levanon, E.Y. Human housekeeping genes, revisited. Trends Genet. 2013, 29, 569–574. [Google Scholar] [CrossRef]
  9. Schmidt, G.W.; Delaney, S.K. Stable internal reference genes for normalization of real-time RT-PCR in tobacco (Nicotiana tabacum) during development and abiotic stress. Mol. Genet. Genom. 2010, 283, 233–241. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, D.; Shi, L.; Han, C.; Yu, J.; Li, D.; Zhang, Y. Validation of Reference Genes for Gene Expression Studies in Virus-Infected Nicotiana benthamiana Using Quantitative Real-Time PCR. PLoS ONE 2012, 7, e46451. [Google Scholar] [CrossRef]
  11. Rotenberg, D.; Thompson, T.S.; German, T.L.; Willis, D.K. Methods for effective real-time RT-PCR analysis of virus-induced gene silencing. J. Virol. Methods 2006, 138, 49–59. [Google Scholar] [CrossRef]
  12. Inoue, T.; Tsai, B. How Viruses Use the Endoplasmic Reticulum for Entry, Replication, and Assembly. Cold Spring Harb. Perspect. Biol. 2013, 5, a013250. [Google Scholar] [CrossRef]
  13. Hyodo, K.; Okuno, T. Host factors used by positive-strand RNA plant viruses for genome replication. J. Gen. Plant Pathol. 2014, 80, 123–135. [Google Scholar] [CrossRef]
  14. Sirover, M.A. New nuclear functions of the glycolytic protein, glyceraldehyde-3-phosphate dehydrogenase, in mammalian cells. J. Cell. Biochem. 2005, 95, 45–52. [Google Scholar] [CrossRef]
  15. Sirover, M.A. New insights into an old protein: The functional diversity of mammalian glyceraldehyde-3-phosphate dehy-drogenase. Biochim. Biophys. Acta 1999, 1432, 159–184. [Google Scholar] [CrossRef]
  16. Bonafé, N.; Gilmore-Hebert, M.; Folk, N.L.; Azodi, M.; Zhou, Y.; Chambers, S.K. Glyceraldehyde-3-Phosphate Dehydrogenase Binds to the AU-Rich 3′ Untranslated Region of Colony-Stimulating Factor–1 (CSF-1) Messenger RNA in Human Ovarian Cancer Cells: Possible Role in CSF-1 Posttranscriptional Regulation and Tumor Phenotype. Cancer Res. 2005, 65, 3762–3771. [Google Scholar] [CrossRef]
  17. Wang, R.Y.; Nagy, P.D. Tomato bushy stunt virus co-opts the RNA-binding function of a host metabolic enzyme for viral ge-nomic RNA synthesis. Cell Host Microbe 2008, 3, 178–187. [Google Scholar] [CrossRef]
  18. Wang, X.; Ahlquist, P. Filling a GAP(DH) in Asymmetric Viral RNA Synthesis. Cell Host Microbe 2008, 3, 124–125. [Google Scholar] [CrossRef] [PubMed]
  19. Gross, S.R.; Kinzy, T.G. Translation elongation factor 1A is essential for regulation of the actin cytoskeleton and cell morphology. Nat. Struct. Mol. Biol. 2005, 12, 772–778. [Google Scholar] [CrossRef]
  20. Nishikiori, M.; Dohi, K.; Mori, M.; Meshi, T.; Naito, S.; Ishikawa, M. Membrane-Bound Tomato Mosaic Virus Replication Proteins Participate in RNA Synthesis and Are Associated with Host Proteins in a Pattern Distinct from Those That Are Not Membrane Bound. J. Virol. 2006, 80, 8459–8468. [Google Scholar] [CrossRef]
  21. Shin, R.; Park, C.J.; Paek, K.H. A translation elongation factor 1A (CaEF1A) gene from hot pepper (Capsicum annuum L.) is induced by the tobacco mosaic virus and by wounding. J. Plant Biol. 2001, 44, 199–204. [Google Scholar] [CrossRef]
  22. Jarošová, J.; Kundu, J.K. 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. [Google Scholar] [CrossRef] [PubMed]
  23. Goodin, M.M.; Zaitlin, D.; Naidu, R.A.; Lommel, S.A. Nicotiana benthamiana: Its History and Future as a Model for Plant–Pathogen Interactions. Mol. Plant-Microbe Interact. 2008, 21, 1015–1026. [Google Scholar] [CrossRef] [PubMed]
  24. Bally, J.; Jung, H.; Mortimer, C.; Naim, F.; Philips, J.G.; Hellens, R.; Bombarely, A.; Goodin, M.M.; Waterhouse, P.M. The Rise and Rise of Nicotiana benthamiana: A Plant for All Reasons. Annu. Rev. Phytopathol. 2018, 56, 405–426. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, S.; Fernald, R.D.; Bylund, J.B.; Trinh, L.T.; Awgulewitsch, C.P.; Paik, D.T.; Jetter, C.; Jha, R.; Zhang, J.; Nolan, K.; et al. Comprehensive Algorithm for Quantitative Real-Time Polymerase Chain Reaction. J. Comput. Biol. 2005, 12, 1047–1064. [Google Scholar] [CrossRef]
  26. Wan, H.; Zhao, Z.; Qian, C.; Sui, Y.; Malik, A.A.; 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–261. [Google Scholar] [CrossRef]
  27. Guénin, S.; Mauriat, M.; Pellouxn, 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, 487–493. [Google Scholar] [CrossRef]
  28. Huang, Y.-W.; Hu, C.-C.; Lin, N.-S.; Hsu, Y.-H. Unusual roles of host metabolic enzymes and housekeeping proteins in plant virus replication. Curr. Opin. Virol. 2012, 2, 676–682. [Google Scholar] [CrossRef]
  29. Zhou, L.; Meng, J.Y.; Ruan, H.Y.; Yang, C.L.; Zhang, C.Y. Expression stability of candidate RT-qPCR housekeeping genes in Spodoptera frugiperda (Lepidoptera: Noctuidae). Arch. Insect Biochem. Physiol. 2021, 108, e21831. [Google Scholar]
  30. Lilly, S.T.; Drummond, R.S.M.; Pearson, M.N.; MacDiarmid, R.M. Identification and Validation of Reference Genes for Normalization of Transcripts from Virus-Infected Arabidopsis thaliana. Mol. Plant-Microbe Interact. 2011, 24, 294–304. [Google Scholar] [CrossRef]
  31. Remans, T.; Smeets, K.; Opdenakker, K.; Mathijsen, D.; Vangronsveld, J.; Cuypers, A. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 2008, 227, 1343–1349. [Google Scholar] [CrossRef] [PubMed]
  32. Hemler, M.E. Tetraspanin functions and associated microdomains. Nat. Rev. Mol. Cell Biol. 2005, 6, 801–811. [Google Scholar] [CrossRef] [PubMed]
  33. Florin, L.; Lang, T. Tetraspanin Assemblies in Virus Infection. Front. Immunol. 2018, 9, 1140. [Google Scholar] [CrossRef]
  34. Zhu, T.; Sun, Y.; Chen, X. Arabidopsis Tetraspanins Facilitate Virus Infection via Membrane-Recognition GCCK/RP Motif and Cysteine Residues. Front. Plant Sci. 2022, 13, 805633. [Google Scholar] [CrossRef]
  35. Kumar, V.; Sharma, R.; Trivedi, P.; Vyas, G.K.; Khandelwal, V. Traditional and novel references towards systematic normalization of qRT-PCR data in plants. Aust. J. Crop. Sci. 2011, 5, 1455–1468. [Google Scholar]
  36. Ding, S.-W. RNA-based antiviral immunity. Nat. Rev. Immunol. 2010, 10, 632–644. [Google Scholar] [CrossRef]
  37. Helderman, T.A.; Deurhof, L.; Bertran, A.; Richard, M.M.S.; Kormelink, R.; Prins, M.; Joosten, M.H.A.J.; Burg, H.A. Members of the ribosomal protein S6 (RPS6) family act as pro-viral factor for tomato spotted wilt orthotospovirus infectivity in Nicotiana benthamiana. Mol. Plant Pathol. 2021, 23, 431–446. [Google Scholar] [CrossRef] [PubMed]
  38. Summer, H.; Grämer, R.; Dröge, P. Denaturing urea polyacrylamide gel electrophoresis (Urea PAGE). J. Vis. Exp. 2009, 32, e1485. [Google Scholar]
  39. Ismayil, A.; Haxim, Y.; Wang, Y.; Li, H.; Qian, L.; Han, T.; Chen, T.; Jia, Q.; Yihao Liu, A.; Zhu, S.; et al. Cotton Leaf Curl Multan virus C4 protein suppresses both tran-scriptional and post-transcriptional gene silencing by interacting with SAM synthetase. PLoS Pathog. 2018, 14, e1007282. [Google Scholar] [CrossRef]
  40. Pombo, M.A.; Ramos, R.N.; Zheng, Y.; Fei, Z.; Martin, G.B.; Rosli, H.G. Transcriptome-based identification and validation of reference genes for plant-bacteria interaction studies using Nicotiana benthamiana. Sci. Rep. 2019, 9, 1–10. [Google Scholar] [CrossRef]
  41. Li, F.; Zhang, C.; Li, Y.; Wu, G.; Hou, X.; Zhou, X.; Wang, A. Beclin1 restricts RNA virus infection in plants through suppression and degradation of the viral polymerase. Nat. Commun. 2018, 9, 1–17. [Google Scholar] [CrossRef]
  42. Wang, Y.; Gong, Q.; Wu, Y.; Huang, F.; Ismayil, A.; Zhang, D.; Li, H.; Gu, H.; Ludman, M.; Fátyol, K.; et al. A calmodulin-binding transcription factor links calcium signaling to antiviral RNAi defense in plants. Cell Host Microbe 2021, 29, 1393–1406. [Google Scholar] [CrossRef] [PubMed]
  43. Grech-Baran, M.; Witek, K.; Poznański, J.T.; Grupa-Urbańska, A.; Malinowski, T.; Lichocka, M.; Jones, J.D.; Hennig, J. The Rysto immune receptor recognises a broadly conserved feature of potyviral coat proteins. New Phytol. 2022, 235, 1179–1195. [Google Scholar] [CrossRef]
  44. Cui, W.; Wang, S.; Han, K.; Zheng, E.; Ji, M.; Chen, B.; Wang, X.; Chen, J.; Yan, F. Ferredoxin 1 is downregulated by the accumulation of abscisic acid in an ABI5-dependent manner to facilitate rice stripe virus infection in Nicotiana benthamiana and rice. Plant J. 2021, 107, 1183–1197. [Google Scholar] [CrossRef]
  45. Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quan-titative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, Research0034. [Google Scholar] [CrossRef] [PubMed]
  46. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. 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–5250. [Google Scholar] [CrossRef] [PubMed]
  47. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. 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–515. [Google Scholar] [CrossRef] [PubMed]
  48. Silver, N.; Best, S.; Jiang, J.; Thein, S.L. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 2006, 7, 33. [Google Scholar] [CrossRef]
  49. Xie, F.; Xiao, P.; Chen, D.; Xu, L.; Zhang, B. miRDeepFinder: A miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 2012, 80, 75–84. [Google Scholar] [CrossRef] [PubMed]
Plants 12 00857 g001 550
Figure 1. Melting curves of the RT-qPCR products of 28 reference genes tested in this study. (A) Melting curves of 13 genes showing a single peak. (B) Melting curves of the remaining 13 genes showing multiples peaks or red curves. Red curves indicate low RNA abundance in one or more replicates.
Figure 1. Melting curves of the RT-qPCR products of 28 reference genes tested in this study. (A) Melting curves of 13 genes showing a single peak. (B) Melting curves of the remaining 13 genes showing multiples peaks or red curves. Red curves indicate low RNA abundance in one or more replicates.
Plants 12 00857 g001
Plants 12 00857 g002 550
Figure 2. RT-qPCR Ct values of 13 candidate genes evaluated in the context of single infections of 11 viruses in N. benthamiana. Values are given as the mean Ct of three technical replicates of three biological replicates of each viral infection and Mock infection group. In total, six values are presented for each candidate reference gene. The same Mock group was used in all cases. Ct values are inversely proportional to the amount of template. The three lines of the box represent the 25th quartiles, median, and 75th quartiles. The whiskers represent minimum and maximum values.
Figure 2. RT-qPCR Ct values of 13 candidate genes evaluated in the context of single infections of 11 viruses in N. benthamiana. Values are given as the mean Ct of three technical replicates of three biological replicates of each viral infection and Mock infection group. In total, six values are presented for each candidate reference gene. The same Mock group was used in all cases. Ct values are inversely proportional to the amount of template. The three lines of the box represent the 25th quartiles, median, and 75th quartiles. The whiskers represent minimum and maximum values.
Plants 12 00857 g002
Table
Table 1. Amplification efficiency of 13 pairs of primers.
Table 1. Amplification efficiency of 13 pairs of primers.
GenesAmplicon Length (bp)E (%)R2
NbeIF4A135107.40.996
NbLip10399.30.9963
NbL2311099.40.9951
NbKLC81106.40.996
NbEF1α1351100.9874
NbPGK105105.50.9988
NbTspan10494.80.9961
NbRdR6134103.30.9960
NbF-BOX12596.10.9994
NbUBC10092.50.9956
NbACT1451080.9914
NbAGO212191.90.9923
NbGAPDH12596.10.9907
E: PCR amplification efficiency; R2: correlation coefficient.
Table
Table 2. Expression stability ranking of 13 candidate reference genes in the context of single infections by 11 viruses, as calculated by four algorithms.
Table 2. Expression stability ranking of 13 candidate reference genes in the context of single infections by 11 viruses, as calculated by four algorithms.
VirusesGeneAlgorithms
GeNormNormFinderBestKeeperDelta-CT
MRSVRCVSDRSDR
TMVNbeIF4A0.3960.0911.83 0.37 61.03 1
NbLip0.1110.1950.17 0.05 11.06 3
NbL230.4370.0912.30 0.49 71.06 2
NbKLC0.3450.1641.40 0.34 51.06 4
NbEF1α0.1110.2870.46 0.08 21.08 5
NbPGK0.4680.0932.79 0.62 91.10 6
NbTspan0.2130.3380.70 0.16 31.11 7
NbRdR60.2640.4690.74 0.19 41.15 8
NbF-BOX0.5290.2763.14 0.81 101.23 9
NbUBC0.59101.00102.46 0.55 81.42 10
NbACT0.73111.12116.57 1.39 111.69 11
NbAGO21.07123.08128.50 2.24 123.11 12
NbGAPDH1.57134.301316.87 4.00 134.31 13
PMMoVNbeIF4A0.27 60.25 62.66 0.53 60.52 7
NbLip0.32 80.42 93.01 0.79 70.55 8
NbL230.36 100.51 114.37 0.88 120.59 10
NbKLC0.41 120.40 80.48 0.11 10.65 12
NbEF1α0.29 70.28 74.12 0.67 110.49 5
NbPGK0.38 110.54 124.51 0.93 130.62 11
NbTspan0.10 10.09 41.43 0.30 30.50 6
NbRdR60.10 10.05 11.36 0.34 20.49 4
NbF-BOX0.22 50.15 52.47 0.60 50.46 3
NbUBC0.18 30.05 12.35 0.52 40.42 1
NbACT0.20 40.05 13.12 0.57 80.43 2
NbAGO20.62 131.73 133.72 0.97 91.74 13
NbGAPDH0.34 90.43 103.89 0.76 100.56 9
TMGMVNbeIF4A0.53 80.39 63.81 0.76 80.96 7
NbLip0.48 60.21 32.02 0.53 40.86 3
NbL230.17 10.49 95.54 1.13 111.00 9
NbKLC0.57 90.47 81.08 0.25 10.91 6
NbEF1α0.23 30.25 55.75 0.96 120.91 5
NbPGK0.17 10.45 75.43 1.13 100.97 8
NbTspan0.67 101.01 101.09 0.22 21.18 10
NbRdR60.41 50.15 22.51 0.64 70.83 1
NbF-BOX0.51 70.11 12.20 0.53 50.84 2
NbUBC0.74 111.09 112.02 0.43 31.26 11
NbACT0.37 40.23 44.50 0.83 90.89 4
NbAGO20.84 121.47 122.30 0.61 61.54 12
NbGAPDH1.16 132.94 1314.26 3.11 132.96 13
ToBRFVNbeIF4A0.19 30.33 73.14 0.64 100.80 6
NbLip0.48 90.29 60.33 0.10 10.82 7
NbL230.40 70.66 114.10 0.89 110.97 9
NbKLC0.44 80.16 41.03 0.25 30.79 5
NbEF1α0.32 50.04 11.71 0.29 60.74 2
NbPGK0.13 10.28 52.80 0.62 80.78 4
NbTspan0.55 100.62 91.17 0.26 40.99 10
NbRdR60.59 110.63 100.89 0.23 21.00 11
NbF-BOX0.13 10.06 32.08 0.53 70.74 1
NbUBC0.30 40.04 11.54 0.35 50.75 3
NbACT0.36 60.44 83.04 0.62 90.88 8
NbAGO21.06 132.38 136.33 1.70 122.40 13
NbGAPDH0.82 122.12 129.52 2.07 132.14 12
ToMMVNbeIF4A0.46 70.13 23.53 0.71 70.90 2
NbLip0.20 10.82 100.45 0.13 21.09 9
NbL230.62 100.40 65.36 1.12 111.10 10
NbKLC0.29 40.54 70.84 0.20 40.98 5
NbEF1α0.34 50.27 52.18 0.36 50.92 3
NbPGK0.53 80.15 44.07 0.87 80.96 4
NbTspan0.20 10.72 90.54 0.12 11.03 8
NbRdR60.39 60.13 11.97 0.50 60.88 1
NbF-BOX0.57 90.14 33.61 0.92 100.99 6
NbUBC0.26 30.67 80.78 0.17 31.03 7
NbACT0.71 110.86 117.27 1.47 121.36 11
NbAGO20.87 121.89 123.46 0.89 91.93 12
NbGAPDH1.29 133.62 1316.14 3.70 133.63 13
TuMVNbeIF4A0.17 40.10 20.83 0.16 50.64 1
NbLip0.30 70.43 90.78 0.22 40.78 9
NbL230.39 90.25 72.47 0.50 90.77 7
NbKLC0.26 60.14 51.10 0.26 60.70 5
NbEF1α0.22 50.11 41.70 0.27 70.71 6
NbPGK0.09 10.05 10.37 0.08 20.64 2
NbTspan0.12 30.18 60.45 0.10 30.67 4
NbRdR60.09 10.10 30.29 0.07 10.66 3
NbF-BOX0.36 80.27 81.80 0.45 80.77 8
NbUBC0.47 100.88 112.67 0.56 101.04 10
NbACT0.55 110.76 104.45 0.87 111.05 11
NbAGO20.70 121.61 124.46 1.13 121.64 12
NbGAPDH0.96 132.37 1310.29 2.20 132.38 13
ChiRSVNbeIF4A0.15 10.35 91.96 0.39 80.52 7
NbLip0.33 80.34 70.48 0.14 20.56 8
NbL230.47 120.62 122.87 0.57 110.70 12
NbKLC0.37 90.34 81.50 0.35 60.58 9
NbEF1α0.15 10.24 51.96 0.31 90.48 5
NbPGK0.23 40.25 61.51 0.32 70.47 3
NbTspan0.22 30.04 11.07 0.23 40.43 1
NbRdR60.29 60.12 30.79 0.19 30.48 4
NbF-BOX0.31 70.17 41.23 0.29 50.49 6
NbUBC0.27 50.11 20.36 0.08 10.47 2
NbACT0.44 110.53 112.92 0.50 120.68 11
NbAGO20.61 131.38 134.28 1.09 131.41 13
NbGAPDH0.40 100.52 102.59 0.47 100.66 10
ChiVMVNbeIF4A0.11 30.16 41.10 0.22 40.35 2
NbLip0.28 100.20 60.56 0.16 10.45 10
NbL230.17 60.32 101.44 0.28 90.43 8
NbKLC0.25 90.20 50.94 0.22 30.44 9
NbEF1α0.19 70.15 31.87 0.30 110.39 5
NbPGK0.08 10.29 91.38 0.29 80.40 6
NbTspan0.12 40.14 21.13 0.25 50.34 1
NbRdR60.33 120.32 110.88 0.21 20.51 11
NbF-BOX0.21 80.13 11.23 0.29 60.39 4
NbUBC0.08 10.24 71.30 0.29 70.38 3
NbACT0.15 50.28 81.63 0.28 100.40 7
NbAGO20.49 131.33 134.32 1.10 131.34 13
NbGAPDH0.31 110.46 122.15 0.39 120.54 12
PePMoVNbeIF4A0.07 10.04 10.49 0.10 30.72 2
NbLip0.07 10.04 10.28 0.08 10.71 1
NbL230.34 70.50 81.94 0.41 70.89 8
NbKLC0.12 30.07 30.48 0.11 20.73 3
NbEF1α0.41 80.45 72.71 0.46 100.89 7
NbPGK0.29 60.41 61.64 0.35 60.84 6
NbTspan0.14 40.08 40.68 0.15 40.74 4
NbRdR60.53 100.63 102.50 0.62 81.00 10
NbF-BOX0.22 50.13 50.90 0.22 50.77 5
NbUBC0.48 90.59 92.52 0.55 90.96 9
NbACT0.59 110.83 113.16 0.64 111.10 11
NbAGO21.07 132.36 138.10 2.06 132.39 13
NbGAPDH0.84 122.19 127.35 1.54 122.23 12
PVXNbeIF4A0.43 60.50 72.15 0.42 70.67 7
NbLip0.67 130.68 132.46 0.70 90.78 13
NbL230.38 50.45 42.15 0.43 60.65 5
NbKLC0.65 120.67 123.59 0.82 120.76 12
NbEF1α0.62 110.61 104.61 0.78 130.72 10
NbPGK0.50 80.22 22.65 0.53 100.56 2
NbTspan0.13 10.48 60.64 0.14 10.64 4
NbRdR60.54 90.54 82.41 0.58 80.72 9
NbF-BOX0.13 10.42 30.85 0.20 20.61 3
NbUBC0.19 30.59 90.98 0.21 30.70 8
NbACT0.47 70.10 12.09 0.38 50.52 1
NbAGO20.23 40.66 111.17 0.30 40.76 11
NbGAPDH0.58 100.47 53.21 0.61 110.66 6
PoPeVYVNbeIF4A0.10 10.15 31.22 0.25 70.26 4
NbLip0.24 100.31 91.17 0.34 60.37 9
NbL230.22 90.34 121.91 0.40 130.38 10
NbKLC0.27 110.32 101.08 0.25 50.39 12
NbEF1α0.15 50.17 61.53 0.27 120.28 6
NbPGK0.10 10.15 51.23 0.26 80.26 3
NbTspan0.11 30.06 10.88 0.19 40.24 1
NbRdR60.31 130.38 131.48 0.37 110.42 13
NbF-BOX0.17 60.07 20.60 0.15 20.24 2
NbUBC0.19 80.17 60.58 0.13 10.28 7
NbACT0.19 70.15 40.79 0.15 30.27 5
NbAGO20.29 120.33 111.30 0.36 100.38 11
NbGAPDH0.14 40.20 81.29 0.25 90.29 8
M: expression stability value; SV: stability value; CV: coefficient of variation; SD: standard deviation; R: Ranking.
Table
Table 3. Gene expression stability ranked using RefFinder.
Table 3. Gene expression stability ranked using RefFinder.
GenesTobamovirusPotyvirusPotexvirusPolerovirus
TMVPMMoVTMGMVToBRFVToMMVTuMVChiRSVChiVMVPePMoVPVXPoPeVYV
NbeIF4A271074362273
NbLip1935378611310
NbL2341181110812117611
NbKLC567676983129
NbEF1α386355377118
NbPGK812548155644
NbTspan639102411421
NbRdR6712812410101013
NbF-BOX95119974512
NbUBC10211261023955
NbACT1144912111191136
NbAGO2121312131112131313812
NbGAPDH13101312131310121297
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, G.; Zhang, Z.; Wan, Q.; Zhou, H.; Jiao, M.; Zheng, H.; Lu, Y.; Rao, S.; Wu, G.; Chen, J.; et al. Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera. Plants 2023, 12, 857. https://doi.org/10.3390/plants12040857

AMA Style

Zhang G, Zhang Z, Wan Q, Zhou H, Jiao M, Zheng H, Lu Y, Rao S, Wu G, Chen J, et al. Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera. Plants. 2023; 12(4):857. https://doi.org/10.3390/plants12040857

Chicago/Turabian Style

Zhang, Ge, Zhuo Zhang, Qionglian Wan, Huijie Zhou, Mengting Jiao, Hongying Zheng, Yuwen Lu, Shaofei Rao, Guanwei Wu, Jianping Chen, and et al. 2023. "Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera" Plants 12, no. 4: 857. https://doi.org/10.3390/plants12040857

APA Style

Zhang, G., Zhang, Z., Wan, Q., Zhou, H., Jiao, M., Zheng, H., Lu, Y., Rao, S., Wu, G., Chen, J., Yan, F., Peng, J., & Wu, J. (2023). Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera. Plants, 12(4), 857. https://doi.org/10.3390/plants12040857

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop