Validation of reference genes for gene expression studies in Aphis glycines (Hemiptera: Aphididae)

Evaluation of reference genes for quantitative expression analysis in Mylabris sibirica (Coleoptera, Meloidae)

Mylabris sibirica is a hypermetamorphic insect whose adults feed on oilseed rape. However, due to a shortage of effective and appropriate endogenous references, studies on molecular functional genes in Mylabris sibirica, have been tremendously limited. In this study, ten internal reference genes (ACT, ARF1, AK, EF1α, GAPDH, α-TUB, RPL6, RPL13, RPS3 and RPS18) were tested and assessed under four selected treatments including adult ages, adult tissues, temperatures, and sex by RT-qPCR based on five methods (Ct value, geNorm, NormFinder, BestKeeper and RefFinder). Our findings showed that RPL6 and RPL13 were the most optimal internal reference gene combination for gene expression during various adult ages and under diverse temperatures; The combination of RPL6 and RPS18 was recommended to test gene transcription levels under different adult tissues. AK and RPL6 were the best reference genes in male and female adults. RPL6 and RPL13 were the most appropriate reference gene pair to estimate gene expression levels under four different tested backgrounds. The relative transcript levels of a uridine diphosphate (UDP)-N-acetylglucosamine-pyrophosphorylase (MsUAP), varied greatly according to normalization with the two most- and least-suited reference genes. This study will lay the basis for further molecular physiology and biochemistry studies in M. sibirica, such as development, reproduction, sex differentiation, cold and heat resistance.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Nanoparticle-Shielded dsRNA Delivery for Enhancing RNAi Efficiency in Cotton Spotted Bollworm Earias vittella (Lepidoptera: Nolidae)

The spotted bollworm Earias vittella (Lepidoptera: Nolidae) is a polyphagous pest with enormous economic significance, primarily affecting cotton and okra. However, the lack of gene sequence information on this pest has a significant constraint on molecular investigations and the formulation of superior pest management strategies. An RNA-seq-based transcriptome study was conducted to alleviate such limitations, and de novo assembly was performed to obtain transcript sequences of this pest. Reference gene identification across E. vittella developmental stages and RNAi treatments were conducted using its sequence information, which resulted in identifying transcription elongation factor (TEF), V-type proton ATPase (V-ATPase), and Glyceraldehyde -3-phosphate dehydrogenase (GAPDH) as the most suitable reference genes for normalization in RT-qPCR-based gene expression studies. The present study also identified important developmental, RNAi pathway, and RNAi target genes and performed life-stage developmental expression analysis using RT-qPCR to select the optimal targets for RNAi. We found that naked dsRNA degradation in the E. vittella hemolymph is the primary reason for poor RNAi. A total of six genes including Juvenile hormone methyl transferase (JHAMT), Chitin synthase (CHS), Aminopeptidase (AMN), Cadherin (CAD), Alpha-amylase (AMY), and V-type proton ATPase (V-ATPase) were selected and knocked down significantly with three different nanoparticles encapsulated dsRNA conjugates, i.e., Chitosan-dsRNA, carbon quantum dots-dsRNA (CQD-dsRNA), and Lipofectamine-dsRNA conjugate. These results demonstrate that feeding nanoparticle-shielded dsRNA silences target genes and suggests that nanoparticle-based RNAi can efficiently manage this pest.

Selection of Reference Genes for Normalization of qRT‒PCR Analysis in the Soybean Aphid Aphis glycines Matsumura (Hemiptera: Aphididae)

The soybean aphid Aphis glycines Matsumura is a predominant insect pest in Asia and North America and causes great losses to soybean. The release of genome data for A. glycines will facilitate gene function research in the future. However, suitable reference genes for A. glycines under various experimental conditions are scarce. To search for appropriate reference genes for A. glycines, nine candidate reference genes, including Act, α-Tub, β-Tub, RPS12, RPS18, RPL5, RPL27, EF1α, and Fer, were tested under six experimental conditions to evaluate their suitability for use in the normalization of qRT‒PCR data. Results showed that EF1α and RPS12 were optimal for the developmental stages of A. glycines, RPS18 and RPS12 were appropriate for wing dimorphism, β-Tub and RPS18 were suitable for different tissues and RPL5, and α-Tub could be used for normalization at different temperatures. β-Tub and EF1α could be proposed as reference genes for insecticide treatment, and RPL5 and RPS12 were found to be the most stable reference genes in different photoperiods. The results provide appropriate reference genes for analyzing gene expression in A. glycines and contribute to future research on the molecular physiology and biochemistry of A. glycines.

Assessment of Suitable Reference Genes for qRT-PCR Normalization in Eocanthecona furcellata (Wolff)

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is a widely used tool for measuring gene expression; however, its accuracy relies on normalizing the data to one or more stable reference genes. Eocanthecona furcellata (Wolff) is a polyphagous predatory natural enemy insect that preferentially feeds on more than 40 types of agricultural and forestry pests, such as those belonging to the orders Lepidoptera, Coleoptera, and Hymenoptera. However, to our knowledge, the selection of stable reference genes has not been reported in detail thus far. In this study, nine E. furcellata candidate reference genes (β-1-TUB, RPL4, RPL32, RPS17, RPS25, SDHA, GAPDH2, EF2, and UBQ) were selected based on transcriptome sequencing results. The expression of these genes in various samples was examined at different developmental stages, in the tissues of male and female adults, and after temperature and starvation treatments. Five algorithms were used, including ΔCt, geNorm, NormFinder, BestKeeper, and RefFinder, to evaluate reference gene expression stability. The results revealed that the most stable reference genes were RPL32 and RPS25 at different developmental stages; RPS17, RPL4, and EF2 for female adult tissue samples; RPS17 and RPL32 for male adult tissue samples; RPS17 and RPL32 for various temperature treatments of nymphs; RPS17 and RPS25 for nymph samples under starvation stress; and RPS17 and RPL32 for all samples. Overall, we obtained a stable expression of reference genes under different conditions in E. furcellata, which provides a basis for future molecular studies on this organism.

Selection of Reference Genes for RT-qPCR Analysis of Wing Dimorphism in English Grain Aphid, Sitobion avenae (Hemiptera: Aphididae)

The English grain aphid, Sitobion avenae (Fabricius), exhibits classic and dramatic phenotypic plasticity in wing development. Both genetic and environmental inputs contribute to the wing polyphenism in aphids, an extreme form of phenotypic plasticity in which a single genotype produces discrete winged and wingless morphs. Validated reference genes are needed to accurately normalize temporal and spatial variation in gene expression estimates by RT-qPCR. In this research, the stability of 11 candidate reference genes selected from S. avenae transcriptomes was evaluated under an array of abiotic and biotic conditions relevant to wing development. RefFinder, a comprehensive software integrating rankings from delta Ct, BestKeeper, NormFinder, and geNorm, offered a series of reference genes for every experimental condition. Overall, helicase (HEL) and ubiquitin ribosomal protein S27A fusion protein (RpS27) are suited for most of the conditions examined in this study, although exceptions do exist. Specifically, NADH dehydrogenase (Ap-NADH) and 28S ribosomal RNA (28S) are recommended for insecticide and antibiotic treatments, while ribosomal RNA L14 (RPL14) and 18S ribosomal RNA (18S) are selected for density treatment, respectively. This study provides a suite of reference genes to investigate the wing polyphenism in S. avenae, and is important for application of RT-qPCR in future experiments of novel tactics to control aphids.

Identification of suitable reference genes for expression profiling studies using qRT-PCR in an important insect pest, Maruca vitrata

BACKGROUND

Maruca vitrata is one of the potential insect pests that cause devastating losses to legume cultivation worldwide. Gene functional studies facilitate dissecting the molecular mechanisms underlying the infection process and enable devising appropriate molecular strategies to control this insect pest. Expression profiling using quantitative real-time PCR (qRT-PCR) provides insights into the functional characterization of target genes; however, ideal reference genes should be deployed in such studies to nullify the background variation and improve the accuracy of target gene expression. An ideal reference gene should have a stable expression across developmental stages, biological conditions, tissues, or experimental conditions.

METHODS AND RESULTS

Given this, the stability of eight candidate reference genes was evaluated in M. vitrata at different developmental stages, diets, and sexes by qRT-PCR method, and the data was analyzed using four independent algorithms, namely GeNorm, NormFinder, BestKeeper, and ΔCt, and one comprehensive algorithm, RefFinder.

CONCLUSION

The analysis showed that RP49 and RPL13 were the best suitable reference genes for studying target gene expression at different developmental stages. Further, the study identified RP49 and RPL24, and GAPDH and RPL24 as the ideal reference genes in M. vitrata fed with different diets and sexes, respectively. The reference genes reported in the present study will ensure the accuracy of target gene expression, and thus, will serve as an important resource for gene functional studies in M. vitrata.

Selection and Validation of Reference Genes For qRT-PCR Analysis of Rhopalosiphum padi (Hemiptera: Aphididae)

Rhopalosiphum padi (L.) (Hemiptera: Aphididae) is an important cosmopolitan pest in cereal crops. Reference genes can significantly affect qRT-PCR results. Therefore, selecting appropriate reference genes is a key prerequisite for qRT-PCR analyses. This study was conducted to identify suitable qRT-PCR reference genes in R. padi. We systematically analyzed the expression profiles of 11 commonly used reference genes. The ΔCt method, the BestKeeper, NormFinder, geNorm algorithms, and the RefFinder online tool were used to evaluate the suitability of these genes under diverse experimental conditions. The data indicated that the most appropriate sets of reference genes were β-actin and GAPDH (for developmental stages), AK and TATA (for populations), RPS18 and RPL13 (for tissues), TATA and GAPDH (for wing dimorphism), EF-1α and RPS6 (for antibiotic treatments), GAPDH and β-actin (for insecticide treatments), GAPDH, TATA, RPS18 (for starvation-induced stress), TATA, RPS6, and AK (for temperatures), and TATA and GAPDH (for all conditions). Our study findings, which revealed the reference genes suitable for various experimental conditions, will facilitate the standardization of qRT-PCR programs, while also improving the accuracy of qRT-PCR analyses, with implications for future research on R. padi gene functions.

Reference gene validation in Eotetranychus sexmaculatus (Acari: Tetranychidae) feeding on mite-susceptible and mite-resistant rubber tree germplasms

Reliable reference genes are quite important in calculating gene transcript levels by using real-time quantitative reverse transcription-PCR (RT-qPCR). Eotetranychus sexmaculatus is known as a dangerous mite causing significant yield reduction of rubber tree latex; however, selection of appropriate reference genes for validation of target gene expression in E. sexmaculatus has not been conducted yet. In the present study, nine candidate reference genes were analyzed for their expression stability in different life stages of E. sexmaculatus by using common algorithms including comparative ΔCq method, geNorm, NormFinder and BestKeeper. In addition, a comprehensive analysis software (RefFinder) was used to assign an overall final rank for each candidate gene. The results showed that β-actin and β-TUB were the best two reference genes and were subjected to evaluate expression of two protective enzyme genes (EsCu/ZnSOD and EsCAT1) in E. sexmaculatus. We found that the expression of EsCu/ZnSOD and EsCAT1 in E. sexmaculatus feeding on mite-resistant rubber tree germplasm was significantly lower compared with those feeding on mite-susceptible germplasm. These results will facilitate research in revealing molecular mechanisms underlying rubber tree resistance to the spider mite.

Diverse patterns of constitutive and inducible overexpression of detoxifying enzyme genes among resistant Aphis glycines populations

Understanding the mechanisms of pyrethroid resistance is essential to the effective management of pesticide resistance in Aphis glycines Matsumura (Hemiptera: Aphididae). We mined putative detoxifying enzyme genes in the draft genome sequence of A. glycines for cytochrome oxidase P450 (CYP), glutathione-S-transferase (GST) and esterases (E4 and carboxylesterases-CES). Aphids from clonal populations resistant to pyrethroids from three sites in Minnesota, USA, were screened against a diagnostic LC99 concentration of either λ-cyhalothrin or bifenthrin and detoxifying enzyme genes expression in survivors was analyzed by qPCR. Their expression profiles were compared relative to a susceptible clonal population. We found 61 CYP (40 full-length), seven GST (all full-length), seven E4 (five full-length) and three CES (two full-length) genes, including 24 possible pseudogenes. The detoxifying enzymes had different expression profiles across resistant aphid populations, possibly reflecting differences in the genetic background and pyrethroid selection pressures as the number of constitutively overexpressed detoxifying enzyme genes was correlated with the level of resistance. Our findings will strengthen the understanding of the pyrethroid resistance mechanisms in A. glycines.

Identification and Evaluation of Reference Genes for Normalization of Gene Expression in Developmental Stages, Sexes, and Tissues of Diaphania caesalis (Lepidoptera, Pyralidae)

Diaphania caesalis (Walker) is an important boring insect mainly distributed in subtropical and tropical areas and attacked tropical woody grain crops, such as starchy plants of Artocarpus. Quantitative real-time polymerase chain reaction (qRT-PCR) is a powerful approach for investigating target genes expression profiles at the transcriptional level. However, the identification and selection of internal reference genes, which is often overlooked, is the most vital step before the analysis of target gene expression by qRT-PCR. So far, the reliable internal reference genes under a certain condition of D. caesalis have not been investigated. Therefore, this study evaluated the expression stability of eight candidate reference genes including ACT, β-TUB, GAPDH, G6PDH, RPS3a, RPL13a, EF1α, and EIF4A in different developmental stages, tissues and sexes using geNorm, NormFinder and BestKeeper algorithms. To verify the stability of the recommended internal reference genes, the expression levels of DcaeOBP5 were analyzed under different treatment conditions. The results indicated that ACT, RPL13a, β-TUB, RPS3a, and EF1α were identified as the most stable reference genes for further studies on target gene expression involving different developmental stages of D. caesalis. And ACT and EIF4A were recommended as stable reference genes for different tissues. Furthermore, ACT, EF1α, and RPS3a were ranked as the best reference genes in different sexes based on three algorithms. Our research represents the critical first step to normalize qRT-PCR data and ensure the accuracy of expression of target genes involved in phylogenetic and physiological mechanism at the transcriptional level in D. caesalia.

Evaluation and Validation of Reference Genes for Quantitative Real-Time PCR in Helopeltis theivora Waterhouse (Hemiptera: Miridae)

Helopeltis theivora Waterhouse is a predominant sucking pest in many tropic economic crops, such as tea, cocoa and coffee. Quantitative real-time PCR (qRT-PCR) is one of the most powerful tools to analyze the gene expression level and investigate the mechanism of insect physiology at transcriptional level. Gene expression studies utilizing qRT-PCR have been applied to numerous insects so far. However, no universal reference genes could be used for H. theivora. To obtain accurate and reliable normalized data in H. theivora, twelve candidate reference genes were examined under different tissues, developmental stages and sexes by using geNorm, NormFinder, BestKeeper, Delta Ct and RefFinder algorithms, respectively. The results revealed that the ideal reference genes differed across the treatments, and the consensus rankings generated from stability values provided by these programs suggested a combination of two genes for normalization. To be specific, RPS3A and Actin were the best suitable reference genes for tissues, RPL13A and GAPDH were suitable for developmental stages, EF1α and RPL13A were suitable for sexes, and RPL13A and RPS3A were suitable for all samples. This study represents the first systematic analysis of reference genes for qRT-PCR experiments in H. theivora, and the results can provide a credible normalization for qRT-PCR data, facilitating transcript profiling studies of functional genes in this insect.

Screening of reference genes using real-time quantitative PCR for gene expression studies in Neoseiulus barkeri Hughes (Acari: Phytoseiidae)

A stable reference gene is a key prerequisite for accurate assessment of gene expression. At present, the real-time reverse transcriptase quantitative polymerase chain reaction has been widely used in the analysis of gene expression in a variety of organisms. Neoseiulus barkeri Hughes (Acari: Phytoseiidae) is a major predator of mites on many important economically crops. Until now, however, there are no reports evaluating the stability of reference genes in this species. In view of this, we used GeNorm, NormFinder, BestKeeper, and RefFinder software tools to evaluate the expression stability of 11 candidate reference genes in developmental stages and under various abiotic stresses. According to our results, β-ACT and Hsp40 were the top two stable reference genes in developmental stages. The Hsp60 and Hsp90 were the most stable reference genes in various acaricides stress. For alterations in temperature, Hsp40 and α-TUB were the most suitable reference genes. About UV stress, EF1α and α-TUB were the best choice, and for the different prey stress, β-ACT and α-TUB were best suited. In normal conditions, the β-ACT and α-TUB were the two of the highest stable reference genes to respond to all kinds of stresses. The current study provided a valuable foundation for the further analysis of gene expression in N. barkeri.

Evaluation of reference genes for real-time quantitative PCR analysis in southern corn rootworm, Diabrotica undecimpunctata howardi (Barber)

Quantitative reverse transcription PCR (RT-qPCR) is one of the most efficient, reliable and widely used techniques to quantify gene expression. In this study, we evaluated the performance of six southern corn rootworm, Diabrotica undecimpunctata howardi (Barber), housekeeping genes (HKG), β-actin (Actin), β-tubulin (Tubulin), elongation factor 1 alpha (EF1α), glyceraldehyde-3 phosphate dehydrogenase (GAPDH), 40 S ribosomal protein S9 (RpS9) and ubiquitin-conjugating protein (Ubi), under different experimental conditions such as developmental stage, exposure of neonate and adults to dsRNA, exposure of adults to different temperatures, different 3^rd instar larva tissues, and neonate starvation. The HKGs were analyzed with four algorithms, including geNorm, NormFinder, BestKeeper, and delta-CT. Although the six HKGs showed a relatively stable expression pattern among different treatments, some variability was observed. Among the six genes, EF1α exhibited the lowest Ct values for all treatments while Ubi exhibited the highest. Among life stages and across treatments, Ubi exhibited the least stable expression pattern. GAPDH, Actin, and EF1α were among the most stable HKGs in the majority of the treatments. This research provides HKG for accurate normalization of RT-qPCR data in the southern corn rootworm. Furthermore, this information can contribute to future genomic and functional genomic research in Diabrotica species.

Robust reference gene design and validation for expression studies in the large milkweed bug, Oncopeltus fasciatus, upon cardiac glycoside stress

Despite its history as a developmental and evolutionary model organism, gene expression analysis in the large milkweed bug, Oncopeltus fasciatus, has rarely been explored using quantitative real-time PCR. The strength of this method depends greatly on the endogenous controls used for normalization, which are lacking for the milkweed bug system. Here, to fill in this gap in our knowledge, we validated the stability of a set of ten candidate reference genes identified from the O. fasciatus transcriptome, and did so upon exposure to a dietary toxin, a cardiac glycoside, and across four different exposure periods. To increase robustness against gDNA contaminants, genome resources were used to design intron-bridging primers. A comprehensive stability validation by the Bestkeeper, Normfinder, geNorm and comparative ΔCt methods identified ef1a and tubulin as the most stable genes across treatments and time points, whereas 18S rRNA was the most unstable. However, accounting for the temporal scale indicated that time point confined normalizers might enable higher quantification accuracy for treatment comparison. Overall this study demonstrates: (i) a robust RT-qPCR primer design approach is possible for non-model organisms where genome annotation is often incomplete, and (ii) the importance of detailed reference gene stability exploration in multifactorial experimental designs.

Reference genes for gene expression studies by RT-qPCR in Brevipalpus yothersi (Acari: Tenuipalpidae), the mite vector of citrus leprosis virus

Quantitative reverse transcription PCR (RT-qPCR) is a high-throughput method to analyze the transcriptional expression of genes. Currently, no reference genes have been described for evaluating gene expression in Brevipalpus yothersi, the false spider mite, a polyphagous that act as vector of the citrus leprosis virus C (CiLV-C), an important citrus disease. This study aimed to identify the most stable reference genes in B. yothersi. The RT-qPCR expression data for selected genes were evaluated from three conditions: different developmental stages, plant hosts and acquisition of CiLV-C. To analyze the stability of the candidate reference genes we used ΔCq method, GeNorm, NormFinder, BestKeeper and RefFinder. Ubiq and GAPDH are best suited for normalizing gene expression data in viruliferous and non-viruliferous mites. Ubiq, EF1α and GAPDH are the most stable for different developmental stages. RPL13 and RPL32 are the best reference genes for approaches to B. yothersi in different host plants. Considering all the experimental conditions, Ubiq, EF1α, and GAPDH were the most stable genes. Here we developed an accurate and comprehensive RT-qPCR strategy for use in B. yothersi gene expression analysis. These results will improve the understanding of the biology of the false spider mites and their role as virus vectors.

Selection of Reference Genes for Optimal Normalization of Quantitative Real-Time Polymerase Chain Reaction Results for Diaphorina citri Adults

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae), can cause direct damage to citrus trees and is the main vector for the devastating disease, citrus greening disease or huanglongbing. Most molecular studies on this important insect pest use real-time reverse-transcription quantitative polymerase chain reaction (RT-qPCR) to quantify gene expression, including analyzing molecular basis for insecticide resistance in field populations. One critical factor to cause inaccuracy in RT-qPCR results is the lack of appropriate internal reference genes for optimal data normalization. In this study, the expression levels of 10 selected reference genes were evaluated in different tissue samples of psyllid adults and in the insects treated with different temperatures and insecticides. Data were analyzed using different computational algorithms, including Delta Ct, BestKeeper, NormFinder, geNorm, and RefFinder. According to our results, at least two reference genes should be used for the normalization of RT-qPCR data in this insect. The best choices of reference genes for different samples are as follows: ACT1 and Ferritin for different tissue samples, RPS20 and Ferritin for samples treated with different temperatures, TBP and EF1α for samples treated with imidacloprid, and Ferritin and TBP for samples treated with beta-cypermethrin. The reference genes identified in this study should be useful for future studies to analyze the expression patterns of target genes, especially for genes linked with temperature adaptability and insecticide resistance in this insect species in the future.

Comparison of Reliable Reference Genes Following Different Hormone Treatments by Various Algorithms for qRT-PCR Analysis of Metasequoia

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is the most sensitive technique for evaluating gene expression levels. Choosing appropriate reference genes for normalizing target gene expression is important for verifying expression changes. Metasequoia is a high-quality and economically important wood species. However, few systematic studies have examined reference genes in Metasequoia. Here, the expression stability of 14 candidate reference genes in different tissues and following different hormone treatments were analyzed using six algorithms. Candidate reference genes were used to normalize the expression pattern of FLOWERING LOCUS T and pyrabactin resistance-like 8. Analysis using the GrayNorm algorithm showed that ACT2 (Actin 2), HIS (histone superfamily protein H3) and TATA (TATA binding protein) were stably expressed in different tissues. ACT2, EF1α (elongation factor-1 alpha) and HIS were optimal for leaves treated with the flowering induction hormone solution, while Cpn60β (60-kDa chaperonin β-subunit), GAPDH (glyceraldehyde-3-phosphate dehydrogenase) and HIS were the best reference genes for treated buds. EF1α, HIS and TATA were useful reference genes for accurate normalization in abscisic acid-response signaling. Our results emphasize the importance of validating reference genes for qRT-PCR analysis in Metasequoia. To avoid errors, suitable reference genes should be used for different tissues and hormone treatments to increase normalization accuracy. Our study provides a foundation for reference gene normalization when analyzing gene expression in Metasequoia.

Selection of Reference Genes for the Normalization of RT-qPCR Data in Gene Expression Studies in Insects: A Systematic Review

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for quantifying expression levels of targeted genes during various biological processes in numerous areas of clinical and biological research. Selection of appropriate reference genes for RT-qPCR normalization is an elementary prerequisite for reliable measurements of gene expression levels. Here, by analyzing datasets published between 2008 and 2017, we summarized the current trends in reference gene selection for insect gene expression studies that employed the most widely used SYBR Green method for RT-qPCR normalization. We curated 90 representative papers, mainly published in 2013–2017, in which a total of 78 insect species were investigated in 100 experiments. Furthermore, top five journals, top 10 frequently used reference genes, and top 10 experimental factors have been determined. The relationships between the numbers of the reference genes, experimental factors, analysis tools on the one hand and publication date (year) on the other hand was investigated by linear regression. We found that the more recently the paper was published, the more experimental factors it tended to explore, and more analysis tools it used. However, linear regression analysis did not reveal a significant correlation between the number of reference genes and the study publication date. Taken together, this meta-analysis will be of great help to researchers that plan gene expression studies in insects, especially the non-model ones, as it provides a summary of appropriate reference genes for expression studies, considers the optimal number of reference genes, and reviews the average number of experimental factors and analysis tools per study.

Host Plant Infection by Soybean Mosaic Virus Reduces the Fitness of Its Vector, Aphis glycines (Hemiptera: Aphididae)

Coevolutionary interactions between pathogens and their insect vectors can dramatically impact the fitness of herbivorous insects and patterns of plant disease transmission. Soybean mosaic virus (SMV) is a common disease in soybean production worldwide. Infected seed is the primary source of inoculum in fields and the virus is secondarily spread among plants by the soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), in a nonpersistent manner. In this study, we compared the biological fitness of A. glycines colonizing both SMV-infected and uninfected soybean plants. Aphids feeding on SMV-infected soybean seedlings were significantly smaller and lighter than those feeding on uninfected plants across life stages. SMV infection caused delayed development of aphid nymphs on soybean seedlings, but this was more than compensated by a reduction in the pre-reproductive period of apterous adults. The fecundity of A. glycines was reduced when feeding on SMV-infected seedlings, resulting in a lower reproductive rate, a longer generation time, and a slower population doubling time. A smaller proportion of aphid offspring developed into alatae when feeding on SMV-infected soybean seedling, and these took longer to mature than their counterparts on uninfected plants. We infer that SMV infection has significantly negative effects on the biological performance of A. glycines, which may be consistent with the long-term coevolution of SMV, soybean, and A. glycines in the transmission cycle of SMV.

Evaluation of potential reference genes for real-time qPCR analysis in a biparental beetle, Lethrus apterus (Coleoptera: Geotrupidae)

Hormones play an important role in the regulation of physiological, developmental and behavioural processes. Many of these mechanisms in insects, however, are still not well understood. One way to investigate hormonal regulation is to analyse gene expression patterns of hormones and their receptors by real-time quantitative polymerase chain reaction (RT-qPCR). This method, however, requires stably expressed reference genes for normalisation. In the present study, we evaluated 11 candidate housekeeping genes as reference genes in samples of Lethrus apterus, an earth-boring beetle with biparental care, collected from a natural population. For identifying the most stable genes we used the following computational methods: geNorm, NormFinder, BestKeeper, comparative delta Ct method and RefFinder. Based on our results, the two body regions sampled (head and thorax) differ in which genes are most stably expressed. We identified two candidate reference genes for each region investigated: ribosomal protein L7A and RP18 in samples extracted from the head, and ribosomal protein L7A and RP4 extracted from the muscles of the thorax. Additionally, L7A and RP18 appear to be the best reference genes for normalisation in all samples irrespective of body region. These reference genes can be used to study the hormonal regulation of reproduction and parental care in Lethrus apterus in the future.

Selection and validation of reference genes for quantitative real-time PCR analysis under different experimental conditions in the leafminer Liriomyza trifolii (Diptera: Agromyzidae)

Liriomyza trifolii is a highly-invasive leafmining insect that causes significant damage to vegetables and horticultural crops worldwide. Relatively few studies have quantified gene expression in L. trifolii using real-time quantitative PCR (RT-qPCR), which is a reliable and sensitive technique for measuring gene expression. RT-qPCR requires the selection of reference genes to normalize gene expression data and control for internal differences between samples. In this study, nine housekeeping genes from L. trifolii were selected for their suitability in normalizing gene expression using geNorm, Normfinder, BestKeeper, the ΔCt method and RefFinder. HSP21.7, which encodes heat shock protein 21.7, was used as a target gene to validate the expression of candidate reference genes. Results indicated that ACTIN and 18S were optimal for developmental stage and low temperature, TUB and 18S showed the most stable expression for sex, and GAPDH and ACTIN were the best reference genes for monitoring gene expression at high temperature. Selection and validation of appropriate reference genes are critical steps in normalizing gene expression levels, which improve the accuracy and quality of expression data. Results of this study provide vital information on reference genes and is valuable in developing a standardized RT-qPCR protocol for functional genomics research in L. trifolii.

Efficacy of RNA interference knockdown using aerosolized short interfering RNAs bound to nanoparticles in three diverse aphid species

RNA interference (RNAi) has emerged as a promising method for validating gene function; however, its utility in nonmodel insects has proven problematic, with delivery methods being one of the main obstacles. This study investigates a novel method of RNAi delivery in aphids, the aerosolization of short interfering RNA (siRNA)-nanoparticle complexes. By using nanoparticles as a siRNA carrier, the likelihood of cellular uptake is increased, when compared to methods previously used in insects. To determine the efficacy of this RNAi delivery system, siRNAs were aerosolized with and without nanoparticles in three aphid species: Acyrthosiphon pisum, Aphis glycines and Schizaphis graminum. The genes targeted for knockdown were carotene dehydrogenase (tor), which is important for pigmentation in Ac. pisum, and branched chain-amino acid transaminase (bcat), which is essential in the metabolism of branched-chain amino acids in all three aphid species. Overall, we observed modest gene knockdown of tor in Ac. pisum and moderate gene knockdown of bcat in Ap. glycines along with its associated phenotype. We also determined that the nanoparticle emulsion significantly increased the efficacy of gene knockdown. Overall, these results suggest that the aerosolized siRNA-nanoparticle delivery method is a promising new high-throughput and non-invasive RNAi delivery method in some aphid species.

Selection and validation of reference genes for qRT-PCR analysis during biological invasions: The thermal adaptability of Bemisia tabaci MED

The Bemisia tabaci Mediterranean (MED) cryptic species has been rapidly invading to most parts of the world owing to its strong ecological adaptability, which is considered as a model insect for stress tolerance studies under rapidly changing environments. Selection of a suitable reference gene for quantitative stress-responsive gene expression analysis based on qRT-PCR is critical for elaborating the molecular mechanisms of thermotolerance. To obtain accurate and reliable normalization data in MED, eight candidate reference genes (β-act, GAPDH, β-tub, EF1-α, GST, 18S, RPL13A and α-tub) were examined under various thermal stresses for varied time periods by using geNorm, NormFinder and BestKeeper algorithms, respectively. Our results revealed that β-tub and EF1-α were the best reference genes across all sample sets. On the other hand, 18S and GADPH showed the least stability for all the samples studied. β-act was proved to be highly stable only in case of short-term thermal stresses. To our knowledge this was the first comprehensive report on validation of reference genes under varying temperature stresses in MED. The study could expedite particular discovery of thermotolerance genes in MED. Further, the present results can form the basis of further research on suitable reference genes in this invasive insect and will facilitate transcript profiling in other invasive insects.

Selection of reference genes for expression analysis in the potato psyllid, Bactericera cockerelli

The selection of reference genes is a crucial step for quantitative real-time PCR analyses and increasingly the use of more than one reference gene for accurate and reliable normalization is being recommended. In this study, a set of six genes was selected and their stability was assessed in different life stages and female organs of Bactericera cockerelli harbouring or not the bacterial pathogen 'Candidatus Liberibacter solanacearum' (Lso) haplotype B. The stability of each gene was determined using the BestKeeper, NormFinder and GeNorm programs. These analyses identified elongation factor-1a, ribosomal protein subunit L5 and ribosomal protein subunit 18 as the most stable genes to analyse gene expression during the insect life stages irrespective of Lso presence; Lso haplotype B only affected their respective ranking. By contrast, a common set of normalizers could not be found amongst the different female organs tested (bacteriomes, alimentary canals and reproductive organs).

Quantitative RT-PCR Gene Evaluation and RNA Interference in the Brown Marmorated Stink Bug

The brown marmorated stink bug (Halyomorpha halys) has emerged as one of the most important invasive insect pests in the United States. Functional genomics in H. halys remains unexplored as molecular resources in this insect have recently been developed. To facilitate functional genomics research, we evaluated ten common insect housekeeping genes (RPS26, EF1A, FAU, UBE4A, ARL2, ARP8, GUS, TBP, TIF6 and RPL9) for their stability across various treatments in H. halys. Our treatments included two biotic factors (tissues and developmental stages) and two stress treatments (RNAi injection and starvation). Reference gene stability was determined using three software algorithms (geNorm, NormFinder, BestKeeper) and a web-based tool (RefFinder). The qRT-PCR results indicated ARP8 and UBE4A exhibit the most stable expression across tissues and developmental stages, ARL2 and FAU for dsRNA treatment and TBP and UBE4A for starvation treatment. Following the dsRNA treatment, all genes except GUS showed relatively stable expression. To demonstrate the utility of validated reference genes in accurate gene expression analysis and to explore gene silencing in H. halys, we performed RNAi by administering dsRNA of target gene (catalase) through microinjection. A successful RNAi response with over 90% reduction in expression of target gene was observed.

Identification and Validation of Reference Genes for the Normalization of Gene Expression Data in qRT-PCR Analysis in Aphis gossypii (Hemiptera: Aphididae)

To obtain accurate and reliable results from quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) analysis, it is necessary to select suitable reference genes as standards for normalizing target gene expression data. QRT-PCR is a popular analytical methodology for studying gene expression and it has been used widely in studies of Aphis gossypii Glover in recent years. However, there is absence of study on the stability of the expression of reference genes in A. gossypii. In this study, eight commonly used candidate reference genes, including 18S, 28S, β-ACT, GAPDH, EF1α, RPL7, α-TUB, and TBP, were evaluated under various experimental conditions to assess their suitability for use in the normalization of qRT-PCR data. The optimal number of reference genes was determined using the geNorm program, and the suitability of particular reference genes was empirically validated by performing normalizations of expression data for the HSP70 gene. The results showed the most suitable combinations of reference genes for the different experimental conditions. For experiments based on divergent developmental stages, EF1α, β-ACT, and RPL7 are the optimal reference gene combination, both EF1α and β-ACT are the optimal combination used in the experiments of different geographical populations, whereas for experiments of the temperature changes, the combination of GAPDH and RPL7 is optimal, both 18S and β-ACT are an optimal combination for feeding assay experiments. These research results should be useful for the selection of the suitable reference genes to obtain reliable qRT-PCR data in the gene expression study of A. gossypii.

Recommended Reference Genes for Quantitative PCR Analysis in Soybean Have Variable Stabilities during Diverse Biotic Stresses

For real-time reverse transcription-PCR (qRT-PCR) in soybean, reference genes in different tissues, developmental stages, various cultivars, and under stress conditions have been suggested but their usefulness for research on soybean under various biotic stresses occurring in North-Central U.S. is not known. Here, we investigated the expression stabilities of ten previously recommended reference genes (ABCT, CYP, EF1A, FBOX, GPDH, RPL30, TUA4, TUB4, TUA5, and UNK2) in soybean under biotic stress from Bean pod mottle virus (BPMV), powdery mildew (PMD), soybean aphid (SBA), and two-spotted spider mite (TSSM). BPMV, PMD, SBA, and TSSM are amongst the most common pest problems on soybean in North-Central U.S. and other regions. Reference gene stability was determined using three software algorithms (geNorm, NormFinder, BestKeeper) and a web-based tool (RefFinder). Reference genes showed variability in their expression as well as stability across various stressors and the best reference genes were stress-dependent. ABCT and FBOX were found to be the most stable in soybean under both BPMV and SBA stress but these genes had only minimal to moderate stability during PMD and TSSM stress. Expression of TUA4 and CYP was found to be most stable during PMD stress; TUB4 and TUA4 were stable under TSSM stress. Under various biotic stresses on soybean analyzed, GPDH expression was found to be consistently unstable. For all biotic stressors on soybean, we obtained pairwise variation (V2/3) values less than 0.15 which suggested that combined use of the two most stable reference genes would be sufficient for normalization. Further, we demonstrated the utility of normalizing the qRT-PCR data for target genes using the most stable reference genes validated in current study. Following of the recommendations from our current study will enable an accurate and reliable normalization of qRT-PCR data in soybean under biotic stress.

Temperature and Development Impacts on Housekeeping Gene Expression in Cowpea Aphid, Aphis craccivora (Hemiptera: Aphidiae)

Quantitative real-time PCR (qRT-PCR) is a powerful technique to quantify gene expression. To standardize gene expression studies and obtain more accurate qRT-PCR analysis, normalization relative to consistently expressed housekeeping genes (HKGs) is required. In this study, ten candidate HKGs including elongation factor 1 α (EF1A), ribosomal protein L11 (RPL11), ribosomal protein L14 (RPL14), ribosomal protein S8 (RPS8), ribosomal protein S23 (RPS23), NADH-ubiquinone oxidoreductase (NADH), vacuolar-type H+-ATPase (ATPase), heat shock protein 70 (HSP70), 18S ribosomal RNA (18S), and 12S ribosomal RNA (12S) from the cowpea aphid, Aphis craccivora Koch were selected. Four algorithms, geNorm, Normfinder, BestKeeper, and the ΔCt method were employed to evaluate the expression profiles of these HKGs as endogenous controls across different developmental stages and temperature regimes. Based on RefFinder, which integrates all four analytical algorithms to compare and rank the candidate HKGs, RPS8, RPL14, and RPL11 were the three most stable HKGs across different developmental stages and temperature conditions. This study is the first step to establish a standardized qRT-PCR analysis in A. craccivora following the MIQE guideline. Results from this study lay a foundation for the genomics and functional genomics research in this sap-sucking insect pest with substantial economic impact.

Abiotic and biotic stressors causing equivalent mortality induce highly variable transcriptional responses in the soybean aphid

Environmental stress affects basic organismal functioning and can cause physiological, developmental, and reproductive impairment. However, in many nonmodel organisms, the core molecular stress response remains poorly characterized and the extent to which stress-induced transcriptional changes differ across qualitatively different stress types is largely unexplored. The current study examines the molecular stress response of the soybean aphid (Aphis glycines) using RNA sequencing and compares transcriptional responses to multiple stressors (heat, starvation, and plant defenses) at a standardized stress level (27% adult mortality). Stress-induced transcriptional changes showed remarkable variation, with starvation, heat, and plant defensive stress altering the expression of 3985, 510, and 12 genes, respectively. Molecular responses showed little overlap across all three stressors. However, a common transcriptional stress response was identified under heat and starvation, involved with up-regulation of glycogen biosynthesis and molecular chaperones and down-regulation of bacterial endosymbiont cellular and insect cuticular components. Stressor-specific responses indicated heat affected expression of heat shock proteins and cuticular components, whereas starvation altered a diverse set of genes involved in primary metabolism, oxidative reductive processes, nucleosome and histone assembly, and the regulation of DNA repair and replication. Exposure to host plant defenses elicited the weakest response, of which half of the genes were of unknown function. This study highlights the need for standardizing stress levels when comparing across stress types and provides a basis for understanding the role of general vs. stressor specific molecular responses in aphids.

Selection and evaluation of reference genes for expression studies with quantitative PCR in the model fungus Neurospora crassa under different environmental conditions in continuous culture

Neurospora crassa has served as a model organism for studying circadian pathways and more recently has gained attention in the biofuel industry due to its enhanced capacity for cellulase production. However, in order to optimize N. crassa for biotechnological applications, metabolic pathways during growth under different environmental conditions must be addressed. Reverse-transcription quantitative PCR (RT-qPCR) is a technique that provides a high-throughput platform from which to measure the expression of a large set of genes over time. The selection of a suitable reference gene is critical for gene expression studies using relative quantification, as this strategy is based on normalization of target gene expression to a reference gene whose expression is stable under the experimental conditions. This study evaluated twelve candidate reference genes for use with N. crassa when grown in continuous culture bioreactors under different light and temperature conditions. Based on combined stability values from NormFinder and Best Keeper software packages, the following are the most appropriate reference genes under conditions of: (1) light/dark cycling: btl, asl, and vma1; (2) all-dark growth: btl, tbp, vma1, and vma2; (3) temperature flux: btl, vma1, act, and asl; (4) all conditions combined: vma1, vma2, tbp, and btl. Since N. crassa exists as different cell types (uni- or multi-nucleated), expression changes in a subset of the candidate genes was further assessed using absolute quantification. A strong negative correlation was found to exist between ratio and threshold cycle (CT) values, demonstrating that CT changes serve as a reliable reflection of transcript, and not gene copy number, fluctuations. The results of this study identified genes that are appropriate for use as reference genes in RT-qPCR studies with N. crassa and demonstrated that even with the presence of different cell types, relative quantification is an acceptable method for measuring gene expression changes during growth in bioreactors.

RNA-Seq reveals a xenobiotic stress response in the soybean aphid, Aphis glycines, when fed aphid-resistant soybean

BACKGROUND

While much recent research has expanded our understanding of the molecular interactions between aphids and their host plants, it is lacking for the soybean aphid, Aphis glycines. Since its North American invasion, A. glycines has become one of the most damaging insect pests on this important crop. Five soybean genes for host plant resistance to A. glycines have been identified, but populations of A. glycines have already adapted to overcome these resistance genes. Understanding the molecular interactions between resistant soybean and A. glycines can provide clues to its adaptation mechanisms. Here, we used RNA-Sequencing to compare and contrast A. glycines gene expression when fed resistant (Rag1) and susceptible soybean.

RESULTS

Combining results from a previous A. glycines transcriptome, we generated 64,860 high quality transcripts, totaling 41,151,086 bases. Statistical analysis revealed 914 genes with significant differential expression. Most genes with higher expression in A. glycines on resistant plants (N = 352) were related to stress and detoxification such as cytochrome P450s, glutathione-S-transferases, carboxyesterases, and ABC transporters. A total of 562 genes showed lower transcript abundance in A. glycines on resistant plants. From our extensive transcriptome data, we also identified genes encoding for putative salivary effector proteins (N = 73). Among these, 6 effector genes have lower transcript abundance in A. glycines feeding on resistant soybean.

CONCLUSIONS

Overall, A. glycines exhibited a pattern typical of xenobiotic challenge, thereby validating antibiosis in Rag1, presumably mediated through toxic secondary metabolites. Additionally, this study identified many A. glycines genes and gene families at the forefront of its molecular interaction with soybean. Further investigation of these genes in other biotypes may reveal adaptation mechanisms to resistant plants.

Validation of reference housekeeping genes for gene expression studies in western corn rootworm (Diabrotica virgifera virgifera)

Quantitative Real-time PCR (qRT-PCR) is a powerful technique to investigate comparative gene expression. In general, normalization of results using a highly stable housekeeping gene (HKG) as an internal control is recommended and necessary. However, there are several reports suggesting that regulation of some HKGs is affected by different conditions. The western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae), is a serious pest of corn in the United States and Europe. The expression profile of target genes related to insecticide exposure, resistance, and RNA interference has become an important experimental technique for study of western corn rootworms; however, lack of information on reliable HKGs under different conditions makes the interpretation of qRT-PCR results difficult. In this study, four distinct algorithms (Genorm, NormFinder, BestKeeper and delta-CT) and five candidate HKGs to genes of reference (β-actin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; β-tubulin; RPS9, ribosomal protein S9; EF1a, elongation factor-1α) were evaluated to determine the most reliable HKG under different experimental conditions including exposure to dsRNA and Bt toxins and among different tissues and developmental stages. Although all the HKGs tested exhibited relatively stable expression among the different treatments, some differences were noted. Among the five candidate reference genes evaluated, β-actin exhibited highly stable expression among different life stages. RPS9 exhibited the most similar pattern of expression among dsRNA treatments, and both experiments indicated that EF1a was the second most stable gene. EF1a was also the most stable for Bt exposure and among different tissues. These results will enable researchers to use more accurate and reliable normalization of qRT-PCR data in WCR experiments.

Identification and validation of reference genes for quantitative real-time PCR in Drosophila suzukii (Diptera: Drosophilidae)

To accurately evaluate gene expression levels and obtain more accurate quantitative real-time RT-PCR (qRT-PCR) data, normalization relative to reliable reference gene(s) is required. Drosophila suzukii, is an invasive fruit pest native to East Asia, and recently invaded Europe and North America, the stability of its reference genes have not been previously investigated. In this study, ten candidate reference genes (RPL18, RPS3, AK, EF-1β, TBP, NADH, HSP22, GAPDH, Actin, α-Tubulin), were evaluated for their suitability as normalization genes under different biotic (developmental stage, tissue and population), and abiotic (photoperiod, temperature) conditions. The three statistical approaches (geNorm, NormFinder and BestKeeper) and one web-based comprehensive tool (RefFinder) were used to normalize analysis of the ten candidate reference genes identified α-Tubulin, TBP and AK as the most stable candidates, while HSP22 and Actin showed the lowest expression stability. We used three most stable genes (α-Tubulin, TBP and AK) and one unstably expressed gene to analyze the expression of P-glycoprotein in abamectin-resistant and sensitive strains, and the results were similar to reference genes α-Tubulin, TBP and AK, which show good stability, while the result of HSP22 has a certain bias. The three validated reference genes can be widely used for quantification of target gene expression with qRT-PCR technology in D.suzukii.

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

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.

Reduction in fecundity and shifts in cellular processes by a native virus on an invasive insect

Pathogens and their vectors have coevolutionary histories that are intricately intertwined with their ecologies, environments, and genetic interactions. The soybean aphid, Aphis glycines, is native to East Asia but has quickly become one of the most important aphid pests in soybean-growing regions of North America. In this study, we used bioassays to examine the effects of feeding on soybean infected with a virus it vectors (Soybean mosaic virus [SMV]) and a virus it does not vector (Bean pod mottle virus [BPMV]) have on A. glycines survival and fecundity. The genetic underpinnings of the observed changes in fitness phenotype were explored using RNA-Seq. Aphids fed on SMV-infected soybean had transcriptome and fitness profiles that were similar to that of aphids fed on healthy control plants. Strikingly, a significant reduction in fecundity was seen in aphids fed on BPMV-infected soybean, concurrent with a large and persistent downregulation of A. glycines transcripts involved in regular cellular activities. Although molecular signatures suggested a small regulatory RNA pathway defense response was repressed in aphids feeding on infected plants, BPMV did not appear to be replicating in the vector. These results suggest that incompatibilities with BPMV or the effects of BPMV infection on soybean caused A. glycines to allot available energy resources to survival rather than reproduction and other core cellular processes. Ultimately, the detrimental impacts to A. glycines may reflect the short tritrophic evolutionary histories between the insect, plant, and virus.

Microbiome diversity of Aphis glycines with extensive superinfection in native and invasive populations

Associations among insects and microbes can lead to beneficial or parasitic interactions. Using 454 sequencing of 16S RNA genes, we compared microbiome diversity and abundance among field-collected (F) and laboratory-reared (L) populations of the invasive soybean aphid (Aphis glycines), a pest of soybean. Additionally, we screened A. glycines populations from native (Japan, South Korea and China) and invasive regions (North America) to broadly determine the microbiome diversity. Our results suggested that Arsenophonus (relative abundance of 54.6%), Buchnera (38.7%) and Wolbachia (3.7%) were the major bacteria associated with A. glycines. Arsenophonus was the most abundant in F populations but was significantly reduced in L populations; additional bacteria species also had lower relative abundances in L populations. Native and invasive populations were largely similar in bacteria communities and revealed substantial superinfection of Arsenophonus and Wolbachia. The lone exception was a lack of Arsenophonus in A. glycines from Japan. Divergent selection pressures among natural and laboratory populations were inferred as factors driving the differential bacterial communities observed. Our results will allow for improved comparative aphid-symbiont research and broaden our understanding of the interactions among insects, endosymbionts and their environments.

Selection of reference genes from two leafhopper species challenged by phytoplasma infection, for gene expression studies by RT-qPCR

BACKGROUND

Phytoplasmas are phloem-limited phytopathogenic wall-less bacteria and represent a major threat to agriculture worldwide. They are transmitted in a persistent, propagative manner by phloem-sucking Hemipteran insects. For gene expression studies based on mRNA quantification by RT-qPCR, stability of housekeeping genes is crucial. The aim of this study was the identification of reference genes to study the effect of phytoplasma infection on gene expression of two leafhopper vector species. The identified reference genes will be useful tools to investigate differential gene expression of leafhopper vectors upon phytoplasma infection.

RESULTS

The expression profiles of ribosomal 18S, actin, ATP synthase β, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and tropomyosin were determined in two leafhopper vector species (Hemiptera: Cicadellidae), both healthy and infected by "Candidatus Phytoplasma asteris" (chrysanthemum yellows phytoplasma strain, CYP). Insects were analyzed at three different times post acquisition, and expression stabilities of the selected genes were evaluated with BestKeeper, geNorm and Normfinder algorithms. In Euscelidius variegatus, all genes under all treatments were stable and could serve as reference genes. In Macrosteles quadripunctulatus, BestKeeper and Normfinder analysis indicated ATP synthase β, tropomyosin and GAPDH as the most stable, whereas geNorm identified reliable genes only for early stages of infection.

CONCLUSIONS

In this study a validation of five candidate reference genes was performed with three algorithms, and housekeeping genes were identified for over time transcript profiling of two leafhopper vector species infected by CYP. This work set up an experimental system to study the molecular basis of phytoplasma multiplication in the insect body, in order to elucidate mechanisms of vector specificity. Most of the sequences provided in this study are new for leafhoppers, which are vectors of economically important plant pathogens. Phylogenetic indications were also drawn from sequence analysis of these genes.

Molecular characterization and expression analysis of soluble trehalase gene in Aphis glycines, a migratory pest of soybean

In insects, the enzyme trehalase plays a crucial role in energy metabolism, chitin synthesis and possibly during plant-insect interactions. We have characterized a soluble trehalase gene (Tre-1) from cDNA of Aphis glycines, a serious migratory pest of soybean. The full-length cDNA of Tre-1 in A. glycines (AyTre-1) was 2550 bp long with an open reading frame of 1770 bp that encoded for a 589 amino acid residues protein. Sequence assessment and phylogenetic analysis of the putative protein suggested that the selected cDNA belongs to soluble trehalase group. Quantitative PCR (qPCR) analysis in different tissues and developmental stages revealed peak mRNA levels of AyTre-1 in the gut (compared with other tissues assayed) and highest expression in the second instar compared with the other developmental stages assayed. Interestingly, a significantly increased expression of AyTre-1 (1.9-fold, P < 0.05) was observed in the alate morphs compared with that in apterate morphs. However, there was no significant difference in AyTre-1 expression in A. glycines-nymphs fed with resistant and susceptible plants. Expression patterns identified in this study provide a platform to investigate the role of AyTre-1 in physiological activities such as flight and feeding in A. glycines. The characterization of soluble trehalase gene may help to develop novel strategies to manage A. glycines using trehalase inhibitors and using RNA interference for knock-down of AyTre-1 expression.

Validation of reference genes for expression analysis by quantitative real-time PCR in Leptinotarsa decemlineata (Say)

Background

L. decemlineata is an exotic invasive insect pest, and invaded in Xinjiang Uygur autonomous region in China in the 1990s from Kazakhstan. It is a notorious defoliator of potato throughout most of the northern Xinjiang in current, and often causes extremely large yield losses of potato.

Results

The expression stability of nine L. decemlineata house-keeping genes (Actin, ACT1 and ACT2; ADP-ribosylation factor, ARF1 and ARF4; TATA box binding protein, TBP1 and TBP2; ribosomal protein RP4 and RP18; translation elongation factor 1α EF1α) was evaluated by quantitative real-time polymerase chain reaction (qRT-PCR) in seven developmental stages, three larval tissues and two insecticide treatments. The results were analyzed using three software programs: geNorm, NormFinder and BestKeeper. Although there was no consistent ranking observed among the house-keeping genes across the samples, the overall analysis revealed that RP18, RP4, ARF1, and ARF4 were the four most stable house-keeping genes. In contrast, ACT1 and ACT2, two of the most widely used reference genes, had the least stability. Our results suggest that the combined use of the four most stably expressed genes may produce optimal normalization for qRT-PCR.

Conclusions

The expression stability of the house-keeping genes varies among different developing stages, in different tissues and under different experimental conditions. Our results will enable a more accurate and reliable normalization of qRT-PCR data in L. decemlineata.

Core RNAi Machinery and Sid1, a Component for Systemic RNAi, in the Hemipteran Insect, Aphis glycines

RNA interference (RNAi) offers a novel tool to manage hemipteran pests. For the success of RNAi based pest control in the field, a robust and systemic RNAi response is a prerequisite. We identified and characterized major genes of the RNAi machinery, Dicer2 (Dcr2), Argonaute2 (Ago2), and R2d2 in Aphis glycines, a serious pest of soybean. The A. glycines genome encodes for at least one copy of Dcr2, R2d2 and Ago2. Comparative and molecular evolution analyses (dN/dS) showed that domain regions of encoded proteins are highly conserved, whereas linker (non-domain) regions are diversified. Sequence homology and phylogenetic analyses suggested that the RNAi machinery of A. glycines is more similar to that of Tribolium casteneum as compared to that of Drosophila melanogaster. We also characterized Sid1, a major gene implicated in the systemic response for RNAi-mediated gene knockdown. Through qPCR, Dcr2, R2d2, Ago2, and Sid1 were found to be expressed at similar levels in various tissues, but higher expression of Dcr2, R2d2, and Ago2 was seen in first and second instars. Characterization of RNAi pathway and Sid1 in A. glycines will provide the foundation of future work for controlling one of the most important insect pests of soybean in North America.

Characterization of a chitin synthase encoding gene and effect of diflubenzuron in soybean aphid, Aphis glycines

Chitin synthases are critical enzymes for synthesis of chitin and thus for subsequent growth and development in insects. We identified the cDNA of chitin synthase gene (CHS) in Aphis glycines, the soybean aphid, which is a serious pest of soybean. The full-length cDNA of CHS in A. glycines (AyCHS) was 5802 bp long with an open reading frame of 4704 bp that encoded for a 1567 amino acid residues protein. The predicted AyCHS protein had a molecular mass of 180.05 kDa and its amino acid sequence contained all the signature motifs (EDR, QRRRW and TWGTR) of chitin synthases. The quantitative real-time PCR (qPCR) analysis revealed that AyCHS was expressed in all major tissues (gut, fat body and integument); however, it had the highest expression in integument (~3.5 fold compared to gut). Interestingly, the expression of AyCHS in developing embryos was nearly 7 fold higher compared to adult integument, which probably is a reflection of embryonic molts in hemimetabolus insects. Expression analysis in different developmental stages of A. glycines revealed a consistent AyCHS expression in all stages. Further, through leaf dip bioassay, we tested the effect of diflubenzuron (DFB, Dimilin ®), a chitin-synthesis inhibitor, on A. glycines' survival, fecundity and body weight. When fed with soybean leaves previously dipped in 50 ppm DFB solution, A. glycines nymphs suffered significantly higher mortality compared to control. A. glycines nymphs feeding on diflubenzuron treated leaves showed a slightly enhanced expression (1.67 fold) of AyCHS compared to nymphs on untreated leaves. We discussed the potential applications of the current study to develop novel management strategies using chitin-synthesis inhibitors and using RNAi by knocking down AyCHS expression.