key: cord-0696916-2gsso3p5
authors: Maroniche, Guillermo A; Sagadín, Mónica; Mongelli, Vanesa C; Truol, Graciela A; del Vas, Mariana
title: Reference gene selection for gene expression studies using RT-qPCR in virus-infected planthoppers
date: 2011-06-16
journal: Virol J
DOI: 10.1186/1743-422x-8-308
sha: 80a8911f3ccc06ddfc7581569014032487688283
doc_id: 696916
cord_uid: 2gsso3p5

BACKGROUND: Planthoppers not only severely affect crops by causing mechanical damage when feeding but are also vectors of several plant virus species. The analysis of gene expression in persistently infected planthoppers might unveil the molecular basis of viral transmission. Quantitative real-time RT-PCR (RT-qPCR) is currently the most accurate and sensitive method used for quantitative gene expression analysis. In order to normalize the resulting quantitative data, reference genes with constant expression during the experimental procedures are needed. RESULTS: Partial sequences of the commonly used reference genes actin (ACT), α1-tubulin (TUB), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), elongation factor 1 alpha (EF1A), ribosomal protein S18 (RPS18) and polyubiquitin C (UBI) from Delphacodes kuscheli, a planthopper capable of persistently transmitting the plant fijivirus Mal de Río Cuarto virus (MRCV), were isolated for the first time. Specific RT-qPCR primers were designed and the expression stability of these genes was assayed in MRCV-infective and naïve planthoppers using geNorm, Normfinder and BestKeeper tools. The overall analysis showed that UBI, followed by 18S and ACT, are the most suitable genes as internal controls for quantitative gene expression studies in MRCV-infective planthoppers, while TUB and EF1A are the most variable ones. Moreover, EF1A was upregulated by MRCV infection. CONCLUSIONS: A RT-qPCR platform for gene expression analysis in the MRCV-infected planthopper vector Delphacodes kuscheli was developed. Our work is the first report on reference gene selection in virus-infected insects, and might serve as a precedent for future gene expression studies on MRCV and other virus-planthopper pathosystems.

Insect genomics is a relatively new field of research. Since the report of the fruit fly genome sequence in 2000 [1] , several other insect genomes have been completed including three mosquitoes species, silkworm, honeybee, red flour beetle, pea aphid, three closely related parasitoid wasps and human body louse [2] [3] [4] [5] . A considerable number of other insect genomes are in progress [6] . In parallel, high throughput tools for functional studies in insects are being developed. For example, Gateway-based vectors applicable for the subcellular localization analysis of proteins in cultured Bombyx mori cells became available [7] . Also important, new insect promoters for gene expression studies are being isolated and characterized [8] and an insect two-hybrid system to analyze protein interactions in cultured insect cells has been reported [9] . Planthoppers (order Hemiptera, superfamily Fulgoroidea) are severe pests of plants because of their sucking damage and ability to transmit at least 18 phytopathogenic viruses [10] . For example, the brown planthopper Nilaparvata lugens is a serious pest of rice in tropical Asia and propagatively transmits Rice ragged stunt virus (RRSV, Oryzavirus genus of the Reoviridae family), Rice grassy stunt virus (RGSV, genus Tenuivirus) [11] and Nilaparvata lugens reovirus (NLRV, Fijivirus genus of the Reoviridae family) [12] . In turn, the small brown planthopper Laodelphax striatellus transmits Rice stripe virus (RSV, genus Tenuivirus) and Rice black streak dwarf virus (RBSDV, Fijivirus genus of the Reoviridae family) [13, 14] and the planthopper Peregrinus maidis transmits Maize mosaic virus (MMV: Nucleorhabdovirus genus of the Rhabdoviridae family) [15] . Another planthopper species, Delphacodes kuscheli, is a natural vector that propagatively transmits Mal de Río Cuarto virus (MRCV, Fijivirus genus of the Reoviridae family). MRCV causes the most important maize disease in Argentina and its genome sequence has recently been completed [16] [17] [18] [19] .

Although no planthopper genonome sequencing project has yet been announced, great progress has been made with the sequencing of more than 37,000 Nilaparvata lugens ESTs originated from various tissues [20] and of 85,526 unigenes from different developmental stages, sexes and wing forms using short-read sequencing technology combined with a tag-based digital gene expression system [21] . This information will certainly contribute to the unveiling of molecular events underlying virus transmission by insects. In particular, very little is known regarding the effects of persistent virus infection on planthoppers transcriptome and new robust and accurate methods for gene expression studies will be required. Most recently, massive Laodelphax striatellus parallel pyrosequencing-based transcriptome analyses was performed in both naïve and RSV-infected insects and several planthopper genes diferentially expressed during virus infection were identified [22] .

The relative quantification of mRNA levels by RT-qPCR is currently an extensively used technique, in which reliable quantification depends on the use of one or more stably expressed endogenous genes, usually housekeeping genes, as internal controls. However, because housekeeping gene expression is not always stable and might vary depending on the samples or treatments, it is necessary to first study the stability of several endogenous gene expression in order to select suitable internal references [23] . Thus, stability analysis of reference genes in diverse conditions and organisms are rapidly increasing [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] . Even though the behavior of insect housekeeping genes has been studied in bacterially challenged bees [30] and in Tribolium beetles infected with fungus [28] , this kind of information is lacking for virus-infected insects. In addition, there is no information currently available on reference gene stability in Delphacodes kuscheli or any other planthopper for accurate quantification of mRNA levels by RT-qPCR. Our study aimed to examine the stability of potential reference genes expression upon MRCV infection of single D. kuscheli planthoppers and to rank them according to their reliability as internal controls. As a result, a robust RT-qPCR platform was set up for analyzing changes on gene expression upon virus infection in individual planthoppers.

The RC MRCV isolate from the endemic disease area (County of Río Cuarto, Province of Córdoba, Argentina) was maintained on oat (Avena sativa L.) plants through successive transmissions by D. kuscheli as described previously [35] . To obtain MRCV-infective insects, 5 III instar D. kuscheli nymphs were fed on a MRCV-infected oat plant for two days. Next, the insects were moved to a cage and fed on uninfected wheat (Triticum aestivum cv Prointa federal) seedlings (displaying one fully developed leaf) for 19 days to achieve the optimum latency period [35] . Subsequently, each insect was placed on an individual uninfected wheat seedling for 24 h (inoculation period) and then fixed in absolute EtOH until use. Finally, after 25 days, the MRCV transmission to the wheat plants was individually determined by monitoring the presence of virus symptoms and by DAS-ELISA.

Total RNA was extracted from each planthopper using the Trizol basic protocol (Invitrogen, CA, USA) with some modifications. First, previous to the extraction the sample was pulverized by freezing in liquid nitrogen and homogenized in a 1.5 ml microcentrifuge tube, and 800 μl of Trizol were added. Second, during the isopropanol precipitation 20 μg of glycogen (Sigma, USA) was added as a carrier. RNA concentration and purity were measured with a spectrophotometer (Thermo Scientific NanoDrop™ 1000, USA) and the integrity was checked by agarose gel electrophoresis. Finally, cDNA was synthesized using Superscript III and random primers from 500 ng of DNAse I (Invitrogen, CA, USA)-treated total RNA according to the manufacturer's instructions.

For the isolation of Delphacodes kuscheli candidate gene sequences, degenerate oligonucleotides were designed from conserved regions of publicly available insect sequences of the actin (ACT), α-1-tubulin (TUB), elongation factor 1-α (EF1A), ribosomal protein S18 (RPS18), gliceraldehyde-3-phosphate dehydrogenase (GAPDH) and polyubiquitin C (UBI) genes. The using D. kuscheli total RNA as template. The resulting amplified fragments were purified using QIAquick gel extraction kit (QIAGEN, Germany), cloned into pGemT Easy (Promega, USA) and sequenced. All Delphacodes kuscheli sequences were deposited in the GenBank database (Table 1) . Finally, real time RT-qPCR specific oligonucleotides were designed from D. kuscheli sequences using Primer Express 3 software (Applied Biosystems, USA). The S18 rRNA primers were designed from the Laodelphax striatellus sequence [GenBank: AB085211]. The sequences of primers used for RT-qPCR analysis are listed on Table 1 .

RT-qPCR experiments were carried out using the Quan-tiTec SYBR Green PCR kit (QIAGEN, Germany) according to the manufacturer's instructions, in an ABI7500 Real Time PCR System (Applied Biosystems, USA) with a standard program (1 min elongation time at 60°C). Reactions were performed in a 25 μl final volume reaction, using primers in a final concentration of 200 nM and 1 μl of a 1/10 dilution of the cDNA as template. No template was added to negative control reactions. Output results were exported as raw data, and introduced to the LinReg software [36] for baseline correction and PCR efficiency calculations. Resulting Ct values were used for expression stability analysis using the Microsoft Excel based tools geNorm 3.5 [37] , Normfinder 0.953 [38] and Bestkeeper v1 [39] according to the developer's instructions. Statistical significance of Ct differences between treatments was calculated by the Mann-Whitney t test using the Graphpad Prism 5 software.

Obtainment of MRCV-infective planthoppers and sample processing

The capacity of MRCV to replicate in plant phloem as well as in its planthopper vector tissues [35] suggests that the virus has the ability to regulate gene expression in both hosts. MRCV replication in plant phloem cells (where it produces severe cytopathic effects) and in planthopper vector tissues (where it establishes a persistent propagative infection) probably involves a differential regulation of gene expression in each case. To study these still unknown gene regulation processes in insects, reliable reference genes must be selected as internal controls to perform gene expression analysis. Consequently, planthoppers undoubtedly able to transmit MRCV are needed to analyze variations on potential reference gene expression levels with respect to naïve insects. One-to-one transmission assays were performed to obtain MRCV-infective insects as described in Methods. As a result, only 4 out of 50 independently evaluated insects were capable of successfully transmitting MRCV to wheat. This transmission efficiency was similar to the reported previously for this kind of assays [35] . Only planthoppers able to transmit MRCV to wheat plants (infective insects) were kept for subsequent RT-qPCR studies along with control naïve insects subjected to the same procedure but always fed on uninfected plants.

In parallel, four different methods for insect storage were tested: immediate freezing in liquid nitrogen List of D. kuscheli selected candidate genes names, assigned function, GenBank database accession codes, designed RT-qPCR primer sequences, length of the amplified fragment (A) and PCR efficiency as calculated by LinReg (E). ND: not determined.

followed by storage at -70°C or stored for 10 days in absolute ethanol, acetone or Trizol. Next, RNA extraction was performed from each planthopper using a Trizol protocol adapted to the small size of the samples (see Methods). After measuring the resulting RNA concentration and analyzing its integrity, ethanol storage was chosen for further experiments because it was a simpler procedure and gave rise to good RNA quality (Additional file 1).

Isolation of partial sequences of candidate Delphacodes kuscheli genes for the design of qPCR primers

Seven housekeeping genes implicated in diverse cellular processes and commonly used as internal references for RT-qPCR in animals were selected as candidates for expression stability analysis. Partial sequences of Delphacodes kuscheli genes homologous to these candidate genes were obtained by RT-PCR with degenerated oligonucleotide primers (see Methods). The amplified actin (ACT), α-1-tubulin (TUB), elongation factor 1-α (EF1A), ribosomal protein S18 (RPS18), glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and polyubiquitin C (UBI) D. kuscheli sequences showed 95%, 86%, 92%, 92%, 90% and 84% identity to Nilaparvata lugens ESTs sequences, respectively, and were deposited in the Gen-Bank database (Table 1) . Next, RT-qPCR oligonucleotides were designed for each of the D. kuscheli isolated gene fragments. Owed to the high conservation of ribosomal RNA sequences, the 18S rRNA primers were based on an available Laodelphax striatellus sequence [GenBank: AB085211]. The designed RT-qPCR primers ( Table 1) were able to efficiently amplify a single product of the expected size and did not show primer dimers formation in RT-PCR experiments (Additional file 2). In addition, a primer amplification performance test by real time RT-qPCR showed that all primer combinations yield a single amplicon and the absence of primer dimers formation in the dissociation step (data not shown). The PCR amplification efficiency of each set of primers was evaluated using the LinReg software which performs, for each individual sample, a baseline correction and a PCR efficiency calculation from the linear portion of the amplification curve slope by linear regression and then calculates a mean PCR efficiency for each primer pair [36] . The calculated PCR efficiency means ranged from 82% for UBI and TUB to 55% for 18S ( Table 1 ). The 18S primer pair low PCR efficiency might have been caused by incomplete annealing due to mismatches, given that their design was based on a related species sequence.

Expression stability analysis of the candidate genes and optimum number of genes for normalization A real time RT-qPCR assay was carried out to examine the expression stability of the 7 candidate genes by the following procedure. Total RNA was extracted from the previously obtained MRCV-infective planthoppers along with control (naïve) individuals, and its concentration was precisely quantified with a spectrophotometer. As a mean for normalization of the sample input, equal amounts of total RNA (500 ng) from each planthopper were treated with DNAse I to remove possible genomic DNA contamination and then used for cDNA synthesis by retrotranscription as detailed previously. Parallel reactions without RT enzyme were performed and used for RT-PCR amplification as an additional control to check that no contaminating genomic DNA was present (Additional file 3). Finally, RT-qPCR experiments were carried as described in Methods. The full sample set including both infective and naïve planthoppers was subjected to each gene expression analysis to avoid between-run variations and three independent technical replicates were performed for each sample in all the experiments. The raw data obtained from each experiment was processed using the LinReg software for baseline correction and the resulting cycle threshold (Ct) values were used in candidate gene expression stability analysis by three different methods using geNorm [37] , Normfinder [38] and Bestkeeper [40] tools.

GeNorm is a Visual Basic Application (VBA) for Microsoft Excel that ranks the genes under study by assigning a stability index M to each of them based on the average pair-wise variation of a gene compared to the rest of the studied genes [37] . According to the ranking of stability generated, ACT was the most stable gene with M= 0.4738, followed by UBI, GADPH, 18S and RPS18 with similar M coefficients. On the other hand, EF1A and TUB were listed as the most unstable genes with M > 1.2 ( Figure 1A ). Next, by performing a stepwise exclusion of the least stable gene and recalculating the M value for the remaining genes, geNorm detected that the most stable pair of reference genes was ACT/GAPDH (Figure 2A ) which displayed the first and third place in the individual stability ranking. Finally, starting with the two most stable genes, the program calculates the pair-wise variation V of two consecutive normalization factors (NF) that result from stepwise introduction of another gene (the one with the lowest M remaining) [37] . This analysis showed that if we consider the generally adopted threshold of V = 0.15 [37, 41, 42] , the inclusion of a third gene does not result in any appreciable improvement of the normalization factor. Moreover, the V values resulting from adding a third, fourth and fifth gene were low (V < 0.015) and the threshold was only reached when adding the two most unstable genes ( Figure 2B ). Thus, these results indicated that only the 2 most stable genes are needed for a reliable normalization.

The Normfinder VBA program analyzes gene expression stability by means of a mathematical model, focusing in the inter-and intra-group variation [38] . This approach ranked UBI as the most stable gene (stability value = 0.1754) followed by 18S, ACT and RPS18 ( Figure 1B) . This is in accordance to geNorm results, which situated UBI as the second most stable gene after ACT. However, Normfinder ranked GAPDH as the most unstable gene, whereas geNorm considered it the third most stable. This difference, due to the singular processing and assumptions of each statistical approach, indicates that it is strongly recommended to carry out more than one analysis when studying expression stability of putative reference genes.

Finally, BestKeeper Excel tool, analyzes each gene's expression variability by calculating the Ct set standard deviation (SD) and coefficient of variance (CV) and then by pair-wise comparison calculates the correlation between the genes and with the Bestkeeper index [40] . This application ranked UBI as the least variable gene with a Ct standard deviation (SD) of 0.31 (1.24 fold change) followed by 18S with SD = 0.39 (1.31 fold change) and indicated that TUB, EF1A and GAPDH were the most variable genes ( Figure 1C ). When pairwise comparing the calculated BestKeeper index against each of the candidate genes, UBI and 18S also showed the higher correlation, with R > 0.89 and p = 0.003 ( Figure 1C ). It is worthy noticing that according to this analysis, none of the genes can be considered inconsistent since all of the SD values were lower than 1 [40] .

BestKeeper and Normfinder yielded comparable results, selecting UBI and 18S as the two most stable genes and GAPDH as the least stable one ( Figure 1B  and 1C) . On the other hand, the geNorm ranking was mostly dissimilar, positioning GAPDH as the third most stable gene ( Figure 1A ). Nevertheless, with the exception of this conflictive gene, all three approaches clearly separated the 7 candidate genes in two distinct groups according to their stability indexes: one group of a low variability and high stability composed of UBI, ACT, 18S and RPS18 and a second one more variable and unstable composed of EF1A and TUB (Figure 1 ). Taken the three method's output together, the results indicate that UBI is the most reliable reference gene for quantitative analysis of MRCV-infective planthopper mRNA levels by RT-qPCR. If two (or more) genes are used for normalization, 18S and ACT are the most suitable choices as a second reference. In agreement, the ACT gene has been recently employed as an internal control for the Tomato spotted wilt virus (TSWV) quantification in individual thrip vectors, although in this case the validation of its expression stability was not previously showed [43] .

Finally, the overall cycle threshold values (Ct) obtained for the different genes were compared. As shown in Figure 3A , the Ct values of the candidate genes were distributed from 15 to 34, in which the lowest Ct (and highest expression) corresponded to 18S and the highest Ct to GAPDH. The high expression level of the 18S rRNA is expected as it is one of the most abundant RNA species in the cell. Moreover, ACT second lowest Ct is in agreement with the usually abundant actin mRNA in cells [20] . It is important to note that the four genes consensually detected as the most reliable (UBI, 18S, ACT and RPS18) covered a relatively wide Ct range (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) . It is convenient to count with several reference genes of diverse Ct values, because the use of internal controls with a Ct close to that of the gene whose expression is being studied is preferable.

Because the expression of genes classically considered as housekeeping has been found to be altered upon virus infection [22, 44] , we analyzed if MRCV infection altered the expression of any of the 7 candidate genes. The Ct values obtained for each gene were compared in MRCV-infective against naïve planthoppers and the statistical significance of Ct differences between treatments was calculated by the Mann-Whitney t test using the Graphpad Prism 5 software. Interestingly, EF1A expression was significantly higher (p < 0.05) in MRCV-infective planthoppers ( Figure 3B ). This fact explains the instability of EF1A gene expression previously detected when using a pool of infective and naïve insects ( Figure  1 ). In agreement, EF1A gene expression has been reported to be 97.8 times more abundant in RSVinfected versus naïve small brown planthoppers [22] . It would be interesting to extend this study to a higher number of Delphacodes kuscheli genes, including antiviral genes of the innate immunity and RNAi pathways, to further establish if their expression is altered during persistent MRCV replication in its planthopper vector.

Insect genomics is a growing area of research. In particular planthoppers transmit many important viral diseases and the molecular basis underling insect-virus interactions are relatively poorly studied. Here, a RT-qPCR platform useful for gene expression analysis in virus infected planthoppers was developed. The evaluation of several D. kuscheli housekeeping genes expression stability showed that the UBI gene is the most reliable as internal reference control for mRNA quantification in MRCV-infective planthoppers. On the other hand, the EF1A gene was upregulated in MRCV-infected planthoppers so it shouldn't be used as an internal control. Our work is the first report on reference gene selection in virus-infected insects, and might serve as a precedent for future gene expression studies on virusplanthopper pathosystems. Particularly, these results will surely contribute to the study of gene expression profiles in D. kuscheli in the context of MRCV infection. The RT-qPCR platform developed along this work will also be useful for normalizing viral titers in insect vectors, contributing to the already developed platforms for viral quantification in vectors and to new molecular studies of transmission and epidemiology [45] [46] [47] [48] [49] [50] . Finally, the study of gene expression patterns in virus-transmitting arthropods will not only improve our knowledge of the virus-vector pathosystems but might also impact in the design of new strategies for virus control.

Additional file 1: Analysis of different sample storage conditions. Four storage conditions were assayed: three individual samples were frozen in liquid nitrogen and kept at -70°C or stored during 10 days in absolute ethanol, acetone or Trizol. Next total RNA was extracted and analyzed by agarose gel electrophoresis. Ribosomal RNA is indicated with an arrow.

Additional file 2: Analysis of the performance of the designed primers in regular RT-PCR reactions. Amplification products of the different candidate genes in five samples (1-5) after regular RT-PCR using the specific primers listed on Table 1 . (-) stands for control RT-PCR reactions where no cDNAs were added. DNA markers are indicated to the right.

Additional file 3: RT-PCR control reactions to assay for genomic DNA contamination on the D. kuscheli total RNA samples. RT-PCR amplification of RPS18 performed in individual samples by using cDNA synthesized with or without the RT enzyme. A control RT-PCR reaction with no cDNA was carried out (-).

List of abbreviations RT-qPCR: real-time reverse transcription polymerase chain reaction.

The genome sequence of Drosophila melanogaster

Genome sequence of the pea aphid Acyrthosiphon pisum

The promise of insect genomics. Pest Management Science

The genome of the model beetle and pest Tribolium castaneum

Functional and evolutionary insights from the genomes of three parasitoid Nasonia species

How can planthopper genomics be useful for planthopper management? In Planthoppers: new threats to the sustainability of intensive rice production systems in Asia

Construction of Gateway-based Destination Vectors for Detecting Subcellular Localization of Proteins in the Silkworm, Bombyx mori

Construction of gene expression systems in insect cell lines using promoters from the silkworm, Bombyx mori

Analysis of protein interactions with two-hybrid system in cultured insect cells

Insect vector interactions with persistently transmitted viruses

Insect-borne viruses of rice

A reovirus in the brown planthopper, Nilaparvata lugens

Maize rough dwarf and related viruses

Studies on insects transmission of rice virus diseases in Japan

Maize fine streak virus, a New Leafhopper-Transmitted Rhabdovirus

Sequence analysis of genome segments S4 and S8 of Mal de Rio Cuarto virus (MRCV): evidence that the virus should be a separate Fijivirus species

Hopp HE, del Vas M: Sequence and phylogenetic analysis of genome segments S1, S2, S3 and S6 of Mal de Rio Cuarto virus, a newly accepted Fijivirus species

Sequence analysis of genome segments S5 and S10 of Mal de Rio Cuarto virus (Fijivirus, Reoviridae)

Sequencing of the bicistronic genome segments S7 and S9 of Mal de Rio Cuarto virus (Fijivirus, Reoviridae) completes the genome of this virus

Annotated ESTs from various tissues of the brown planthopper Nilaparvata lugens: a genomic resource for studying agricultural pests

Transcriptome Analysis of the Brown Planthopper Nilaparvata lugens

Massively parallel pyrosequencing-based transcriptome analyses of small brown planthopper (Laodelphax striatellus), a vector insect transmitting rice stripe virus (RSV)

Real-time RT-PCR normalisation; strategies and considerations

Selection of reference genes for expression studies with fish myogenic cell cultures

Selection of reference genes for mRNA quantification in Trypanosoma brucei

Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process

Selection of reference genes for real-time polymerase chain reaction analysis in tissues from Bombus terrestris and Bombus lucorum of different ages

Evaluation of quantitative PCR reference genes for gene expression studies in Tribolium castaneum after fungal challenge

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 Cytomegalovirus infections

Differential protein expression in the honey bee head after a bacterial challenge

Selection of reference genes for quantitative real-time PCR expression studies in the apomictic and sexual grass Brachiaria brizantha

Identification and validation of housekeeping genes in brains of the desert locust Schistocerca gregaria under different developmental conditions

Selection and validation of reference genes for real-time RT-PCR studies in the non-model species Delomys sublineatus, an endemic Brazilian rodent

Selection of Reference Genes in Transcription Analysis of Gene Expression of the Mandarin Fish, Siniperca chuasti

Transmission features of Mal de Río Cuarto virus in wheat by its planthopper vector Delphacodes kuscheli

Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data

Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes

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

Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper-Excel-based tool using pair-wise correlations

Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper -Excel-based tool using pair-wise correlations

Evaluation of putative reference genes for gene expression normalization in soybean by quantitative real-time RT-PCR

Selection of reliable reference genes for gene expression studies in peach using real-time PCR

Variation in Tomato spotted wilt virus Titer in Frankliniella occidentalis and Its Association with Frequency of Transmission

Transcriptomic analysis of intestinal genes following acquisition of pea enation mosaic virus by the pea aphid Acyrthosiphon pisum

Cytochrome c oxidase mRNA as an internal control for detection of Potato virus Y and Potato leafroll virus from single aphids by a co-amplification RT-PCR assay

A one-step real time RT-PCR assay for quantifying rice stripe virus in rice and in the small brown planthopper (Laodelphax striatellus Fallen)

Detecting Rice stripe virus (RSV) in the small brown planthopper (Laodelphax striatellus) with high specificity by RT-PCR

Detection of wheat dwarf virus (WDV) in wheat and vector leafhopper (Psammotettix alienus Dahlb.) by real-time PCR

Fereres Castiel A: Estimation of vector propensity for Lettuce mosaic virus based on viral detection in single aphids

The detection of Tomato spotted wilt virus (TSWV) in individual thrips using real time fluorescent RT-PCR (TaqMan)

Reference gene selection for gene expression studies using RT-qPCR in virus-infected planthoppers

Authors' contributions GAM performed the stability analysis and expression profiles of the candidate genes. MS maintained the MRCV inocula and performed one-toone transmission assays to obtain MRCV-infective insects. VCM designed degenerated primers and amplified, cloned and analyzed the sequence of D. kuscheli reference genes. GAT supervised the virus transmission experiments and the insect rearing. MdV and GAM participated in the conception and design of the studies and wrote the manuscript. All authors read and approved the final manuscript.

The authors declare that they have no competing interests.