Search

The olfactory system is an important component of insect behavior and is vital for survival and reproduction. However, the genomic characterization and molecular basis of the olfactory response of Megalurothrips usitatus remain relatively unknown. RNA sequencing-built developmental transcriptomes of nymphs, pupae, and adult M. usitatus were examined in order to establish the sequence-based background of M. usitatus olfactory responses. A total of 56,669 unigenes were annotated using GO, NR, Pfam, eggNOG, SwissProt, and KEGG. The number of differentially expressed genes between pupae and nymphs, males and nymphs, and females and nymphs were 10,498, 9,235, and 10,964, respectively. One odorant-binding protein (MusiOBP1) and one chemosensory protein (MusiCSP1) were selected from the transcriptome, and their full-length sequences were obtained using RACE PCR. The relative expression of MusiOBP1 was the highest in primordial females, whereas the relative expression of MusiCSP1 was the highest in primordial pupae. The strongest binding ability to the odor-binding protein MusiOBP1 was observed for β-citronellol. 3-Hydroxy-2-methyl-4-pyrone showed the strongest binding affinity to MusiCSP1. Our analysis suggests that MusiOBP1 and MusiCSP1 may play significant roles in mediating M. usitatus host recognition. This research will improve our knowledge of odorant-binding proteins and chemosensory proteins, which will in turn improve our understanding of insect olfactory systems.

Henosepilachna vigintioctopunctata is an emerging horticultural pest of Solanaceae plants in Asia. Here, we employed transcriptome sequencing of three female and three male H. vigintioctopunctata libraries to identify sex-related genes. We recorded 281 599 294 bp of clean reads, and de novo assembly generated 85 206 unigenes, with N50 of 1 311 bp and average length of 763 bp. We found 5 002 genes highly expressed in males and 2 179 genes highly expressed in females. Expression profiles of six unigenes specific to females, six unigenes specific to males, and ten unigenes common to both sexes were verified using semi-quantitive RT-PCR; results showed expression levels of these 22 unigenes were relatively consistent with the transcriptome analysis. Using all unigenes as references, we identified 6 657 putative simple sequence repeats. The transcriptome sequences and gene expression profiles of female and male H. vigintioctopunctata help explore genetic mechanisms of sexual dimorphism.

Abstract Background While virus-vector-host interactions have been a major focus of both basic and applied ecological research, little is known about how different levels of plant defense interact with prior herbivory to affect these relationships. We used genetically-modified strains of tomato (Solanum lycopersicum) varying in the jasmonic acid (JA) plant defense pathways to explore how plant defense and prior herbivory affects a plant virus (tomato yellow leaf curl virus, ‘TYLCV’), its vector (the whitefly Bemisia tabaci MED), and the host. Results Virus-free MED preferred low-JA over high-JA plants and had lower fitness on high-JA plants. Viruliferous MED preferred low-JA plants but their survival was unaffected by JA levels. While virus-free MED did not lower plant JA levels, viruliferous MED decreased both JA levels and the expression of JA-related genes. Infestation by viruliferous MED reduced plant JA levels. In preference tests, neither virus-free nor viruliferous MED discriminated among JA-varying plants previously exposed to virus-free MED. However, both virus-free and viruliferous MED preferred low-JA plant genotypes when choosing between plants that had both been previously exposed to viruliferous MED. The enhanced preference for low-JA genotypes appears linked to the volatile compound neophytadiene, which was found only in whitefly-infested plants and at concentrations inversely related to plant JA levels. Conclusions Our findings illustrate how plant defense can interact with prior herbivory to affect both a plant virus and its whitefly vector, and confirm the induction of neophytadiene by MED. The apparent attraction of MED to neophytadiene may prove useful in pest detection and management.

Fluralaner, a systemic pesticide, was originally registered with the US Food and Drug Administration in 2014 under the trade name Bravecto for flea treatment for pets. As a GABA antagonist, the footprint of fluralaner has expended beyond medical and veterinary pests in recent years. In this study, we examined the acute toxicity of fluralaner against three pests of Henosepilachna vigintioctopunctata, Megalurothrips usitatus, and Phyllotreta striolata in the Solanaceae, Fabaceae, and Cruciferae families, respectively, and the sublethal impact of fluralaner on Propylaea japonica, a widely distributed predatory ladybeetle. Based on LC50, fluralaner was effective against H. vigintioctopunctata (0.098 mg a.i. L−1 for the second instar larvae), M. usitatus (0.134 mg a.i. L−1 for adult females), and P. striolata (0.595 mg a.i. L−1 for adults). For P. japonica, however, fluralaner was substantially less effective (1.177 mg a.i. L−1 for the third instar larvae). Furthermore, the LC10 and LC30 of P. japonica were also consistently higher than the LC50 of the three pests. In addition, we did not observe any significant impacts of fluralaner at LC10 and LC30 on the life history traits, including body weight, developmental time, pre-oviposition period, and fecundity of P. japonica. Based on our results from acute toxicities and sublethal impacts, fluralaner is effective against vegetable pests, while potentially friendly to P. japonica when employed as a biological control agent.

The gut bacteria of insects positively influence the physiology of their host, however, the dynamics of this complicated ecosystem are not fully clear. To improve our understanding, we characterized the gut prokaryotic of Henosepilachna vigintioctopunctata that fed on two host plants, Solanum melongena (referred to as QZ hereafter) and Solanum nigrum (referred to as LK hereafter), by sequencing the V3-V4 hypervariable region of the 16S rRNA gene using the Illumina MiSeq system. The results revealed that the gut bacterial composition varied between specimens that fed on different host plants. The unweighted pair group method with arithmetic mean analyses and principal coordinate analysis showed that the bacterial communities of the LK and QZ groups were distinct. Four phyla (Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria) were present in all H. vigintioctopunctata gut samples. It is noteworthy that bacteria of the phylum Cyanobacteria were only found in the LK group, with a low relative abundance. Proteobacteria and Enterobacteriaceae were the predominant phylum and family, respectively, in both the LK and QZ groups. Linear discriminant analysis effect size (LEfSe) analyses showed that the QZ group enriched the Bacilli class and Lactococcus genus; while the LK group enriched the Alphaproteobacteria class and Ochrobactrum genus. PICRUSt analysis showed that genes predicted to be involved in xenobiotic biodegradation and metabolism, metabolism of other amino acids, signaling molecules, and interaction were significantly higher in the QZ group. Genes predicted to be involved in the metabolism of cofactors and vitamins were significantly higher in the LK group. Furthermore, the complexity of the network structure and the modularity were higher in the LK group than in the QZ group. This is the first study to characterize the gut bacteria of H. vigintioctopunctat, our results demonstrate that the two host plants tested had a considerable impact on bacterial composition in the gut of H. vigintioctopunctata and that the bacterial communities were dominated by relatively few taxa.

Asian citrus psyllid Diaphorina citri Kuwayama is an important economic pest of citrus, as it transmits Candidatus Liberibacter asiaticus, the causative agent of huanglongbing. In this study, we used RNA-seq to identify novel genes and provide the first high-resolution view of the of D. citri transcriptome throughout development. The transcriptomes of D. citri during eight developmental stages, including the egg, five instars, and male and female adults were sequenced. In total, 115 million clean reads were obtained and assembled into 354,726 unigenes with an average length of 925.65 bp and an N50 length of 1733 bp. Clusters of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses were conducted to functionally annotate the genes. Differential expression analysis highlighted developmental stage-specific expression patterns. Furthermore, two trehalase genes were characterized with lower expression in adults compared to that in the other stages. The RNA interference (RNAi)-mediated suppression of the two trehalase genes resulted in significantly high D. citri mortality. This study enriched the genomic information regarding D. citri. Importantly, these data represent the most comprehensive transcriptomic resource currently available for D. citri and will facilitate functional genomics studies of this notorious pest.

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.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a momentous technique for quantifying expression levels of the targeted genes across various biological processes. Selection and validation of appropriate reference genes for RT-qPCR analysis are a pivotal precondition for reliable expression measurement. Henosepilachna vigintioctopunctata is one of the most serious insect pests that attack Solanaceae plants in Asian countries. Recently, the transcriptomes of H. vigintioctopunctata were sequenced, promoting gene functional studies of this insect pest. Unfortunately, the reference genes for H. vigintioctopunctata have not been selected and validated. Here, a total of 7 commonly used reference genes, namely, Actin, GAPDH, RPL13, RPL6, RPL32, RPS18, and ATPB, were selected and assessed for suitability under four experimental conditions, namely, developmental stage, tissue, temperature, and host plant, using RefFinder, which integrates four different analytical tools (Normfinder, geNorm, the ΔCt method, and BestKeeper). The results displayed that RPL13 and RPS18 were the best suitable reference genes for each experimental condition. The relative transcript levels of 2 target genes, lov and TBX1, varied greatly according to normalization with the two most- and least-suited reference genes. Our results will be helpful for improving the accuracy of the RT-qPCR analysis for future functional investigations of target gene expression in H. vigintioctopunctata.

Diverse pathogens, plant hosts, insect vectors, and non-vector herbivores coexist and interact in natural systems. An example is the cooccurrence of insects Bemisia tabaci Q and Frankliniella occidentalis and the pathogens tomato yellow leaf curl virus (TYLCV) and tomato spotted wilt virus (TSWV) on the same plant. In addition, both TYLCV and TSWV are persistently transmitted in these insect species. However, TSWV reduces the fitness of B. tabaci Q; therefore, we investigated whether TSWV affects the transmission of TYLCV to tomato. Both TYLCV and TSWV are persistently transmitted. Although B. tabaci Q cannot transmit TSWV, we found that this insect species is able to acquire and retain this virus serotype, indicating that the effects of TSWV on TYLCV transmission in the current study result from effects on the vector. The acquisition, retention, and transmission of TYLCV by B. tabaci Q were reduced when the insect vector contained TSWV. Additionally, the TYLCV acquisition and transmission by B. tabaci Q were reduced when the host plant was inoculated with TSWV before TYLCV or simultaneously with TYLCV. We also found that F. occidentalis fecundity and transmission of TSWV were reduced when F. occidentalis contained TYLCV. Our findings are consistent with the hypothesis that persistently transmitted viruses can restrict the transmission of other viruses by affecting their insect vectors.

RNA interference (RNAi) techniques have emerged as powerful tools in the development of novel management strategies for the control of insect pests, such as Henosepilachna vigintioctopunctata, which is a major solanaceous pest in Asia. Our results showed that levels of HvSnf7 expression were greater in larval midguts than in other tissues. Silencing of HvSnf7 led to greater H. vigintioctopunctata mortality rates and appeared to be time- and partially dose-dependent. Bacterially expressed dsHvSnf7 that was applied to detached plant leaves caused 98, 88, and 60% mortality in 1st and 3rd instars, and adults after 10, 12, and 14 d, respectively; when applied to living plants, bacterially expressed dsHvSnf7 led to mortality in 1st and 3rd instars, with no effect on adults. Bacterially expressed dsHvSnf7 led to improved plant protection against H. vigintioctopunctata. Ultrastructural changes caused by HvSnf7-RNAi in larval midguts showed extensive loss of cellular contents that indicate loss of membrane integrity. This study indicate that HvSnf7 potentially can be used as RNAi target gene for controlling of H. vigintioctopunctata.

The sweet potato whitefly, Bemisia tabaci, is one of the most invasive insect pests worldwide. The two most destructive whitefly cryptic species are MEAM1/B and MED/Q. Given that MED/Q has replaced MEAM1/B in China and the invasion of MED/Q has coincided with the outbreak of tomato yellow leaf curl virus (TYLCV), we hypothesize that pre-infestation with viruliferous B. tabaci will affect the subsequent host preferences. To test this hypothesis, we (1) conducted bioassays to compare the host preference of viruliferous and non-viruliferous MEAM1/B and MED/Q, respectively, on plants pre-infested with viruliferous and non-viruliferous MEAM1/B and MED/Q; (2) profiled plant volatiles using GC-MS; and (3) functionally characterized chemical cues could potentially modulate B. tabaci-TYLCV-tomato interactions, including ρ-cymene, thujene and neophytadiene, using a Y-tube olfactometer. As a result, plants pre-infested with MEAM1/B whiteflies carrying TYLCV or not, did not attract more or less B or Q whiteflies. Plants pre-infested with non-viruliferous MED/Q resisted MEAM1/B but did not affect MED/Q. However, plants pre-infested with viruliferous MED/Q attracted more whiteflies. Feeding of viruliferous MED/Q reduced the production of ρ-cymene, and induced thujene and neophytadiene. Functionally analyses of these plant volatiles show that ρ-cymene deters while neophytadiene recruits whiteflies. These combined results suggest that pre-infestation with viruliferous MED/Q promotes the subsequent whitefly infestation and induces plant volatile neophytadiene which recruits whiteflies.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique commonly used in molecular biology to analyze RNA expression. The selection of suitable reference genes for data normalization is a precondition for credible measurements of gene expression levels using RT-qPCR. Propylea japonica is one of the most common pests of many crop systems throughout East Asia, and has often been used in the testing of non-target impacts during environmental risk assessments of genetically engineered plants. The present study assessed the suitability of nine frequently used reference genes for comparisons of P. japonica gene expression. Expression stability was compared across developmental stages, sex, a range of tissues, and following exposure to different temperatures. Data were analyzed using RefFinder, which integrated the results obtained using NormFinder, geNorm, BestKeeper, and the ΔCt method. This led to the identification of unique sets of reference genes for each experimental condition: ribosomal protein S18 (RPS18) and elongation factor 1 α (EF1A) for developmental stage comparisons, RPS18 and EF1A for sex comparisons, EF1A and ribosomal protein L4 for tissue comparisons, and RPS18 and EF1A for analyses of temperature-mediated effects. These reference genes will help to enhance the accuracy of RT-qPCR analyses of P. japonica gene expression. This work represents an initial move towards building a standardized system for RT-qPCR analysis of P. japonica, providing a basis for the ecological risk assessment of RNAi-based insect control products.

Sugar transporters (STs) play pivotal roles in the growth, development, and stress responses of phloem-sucking insects, such as the whitefly, Bemisia tabaci. In this study, 137 sugar transporters (STs) were identified based on analysis of the genome and transcriptome of B. tabaci MEAM1. B. tabaci MEAM1 encodes a larger number of STs than other selected insects. Phylogenetic and molecular evolution analysis showed that the 137 STs formed three expanded clades and that the genes in Sternorrhyncha expanded clades had accelerated rates of evolution. B. tabaci sugar transporters (BTSTs) were divided into three groups based on their expression profiles across developmental stages; however, no host-specific BTST was found in B. tabaci fed on different host plants. Feeding of B. tabaci adults with feeding diet containing dsRNA significantly reduced the transcript level of the target genes in B. tabaci and mortality was significantly improved in B. tabaci fed on dsRNA compared to the control, which indicates the sugar transporters may be used as potential RNAi targets for B. tabaci bio-control. These results provide a foundation for further studies of STs in B. tabaci.

The whitefly (Bemisia tabaci), an important invasive pest that causes severe damage to crops worldwide, has developed resistance to a variety of insecticides. Carboxylesterases (COEs) are important multifunctional enzymes involved in the growth, development, and xenobiotic metabolism of insects. However, systematic studies on the COEs of B. tabaci are scarce. Here, 42 putative COEs in different functional categories were identified in the Mediterranean species of B. tabaci (B. tabaci MED) based on a genome database and neighbor-joining phylogeny. The expression patterns of the COEs were affected by the development of B. tabaci. The expression levels of six COEs were positively correlated with the concentration of imidacloprid to which B. tabaci adults were exposed. The mortality of B. tabaci MED adults fed dsBTbe5 (67.5%) and dsBTjhe2 (58.4%) was significantly higher than the adults fed dsEGFP (41.1%) when treated with imidacloprid. Our results provide a basis for functional research on COEs in B. tabaci and provide new insight into the imidacloprid resistance of B. tabaci.

Henosepilachna vigintioctopunctata (Coleoptera: Coccinellidae) is a major pest affecting Solanaceae plants in Asian countries. In this study, we sequenced the ovary and testis transcriptomes of H. vigintioctopunctata to identify gonad-related genes. Comparison of the unigene sequences in ovary and testis libraries identified 1,421 and 5,315 ovary- and testis-specific genes, respectively. Among the ovary-specific genes, we selected the RC2-like and PSHS-like genes to investigate the effects of gene silencing on the mortality, percentage infertility, pre-oviposition period, fecundity, daily number of eggs laid, and hatching rate in female adults. Although the percentage mortality and infertility of females did not differ significantly among dsRNA treatments, fecundity was significantly reduced in the dsRC2-like and dsPSHS-like treatment groups. Moreover, the pre-oviposition period was markedly prolonged in response to dsPSHS-like treatment. This is the first reported RNA sequencing of H. vigintioctopunctata. The transcriptome sequences and gene expression profiles of the ovary and testis libraries will provide useful information for the identification of gonad-related genes in H. vigintioctopunctata and facilitate further research on the reproductive biology of this species. Moreover, the gonad-specific genes identified may represent candidate target genes for inhibiting the population growth of H. vigintioctopunctata.

The development of genetically-engineered plants that employ RNA interference (RNAi) to suppress invertebrate pests opens up new avenues for insect control. While this biotechnology shows tremendous promise, the potential for both non-target and off-target impacts, which likely manifest via altered mRNA expression in the exposed organisms, remains a major concern. One powerful tool for the analysis of these un-intended effects is RT-qPCR, a technique for quantifying gene expression using a suite of reference genes for normalization. The seven-spotted ladybeetle Coccinella septempunctata, a commonly-used predator in both classical and augmentative biological controls, is a model surrogate species used in the environmental risk assessment (ERA) of plant incorporated protectants (PIPs). Here, we assessed the suitability of eight reference gene candidates for the normalization and analysis of C. septempunctata v-ATPase A gene expression under both biotic and abiotic conditions. Five computational tools with distinct algorisms, geNorm, Normfinder, BestKeeper, the ΔCt method, and RefFinder, were used to evaluate the stability of these candidates. As a result, unique sets of reference genes were recommended respectively for experiments involving different developmental stages, tissues, and ingested dsRNAs. By providing a foundation for standardized RT-qPCR analysis in C. septempunctata, our work improves the accuracy and replicability of the ERA of PIPs involving RNAi transgenic plants.

RNAi-based genetically engineered (GE) crops for the management of insect pests are likely to be commercialized by the end of this decade. Without a workable framework for conducting the ecological risk assessment (ERA) and a standardized ERA protocol, however, the utility of RNAi transgenic crops in pest management remains uncertain. The overall goal of this study is to assess the risks of RNAi-based GE crops on a non-target soil micro-arthropod, Sinella curviseta, which could be exposed to plant-protected dsRNAs deposited in crop residues. Based on the preliminary research, we hypothesized that insecticidal dsRNAs targeting at the western corn rootworm, Diabrotica virgifera virgifera, a billion-dollar insect pest, has no adverse impacts on S. curviseta, a soil decomposer. Following a tiered approach, we tested this risk hypothesis using a well-designed dietary RNAi toxicity assay. To create the worst-case scenario, the full-length cDNA of v-ATPase subunit A from S. curviseta were cloned and a 400 bp fragment representing the highest sequence similarity between target pest and non-target arthropods was selected as the template to synthesize insecticidal dsRNAs. Specifically, 10-day old S. curviseta larvae were subjected to artificial diets containing v-ATPase A dsRNAs from both D. v. virgifera (dsDVV) and S. curviseta (dsSC), respectively, a dsRNA control, β-glucuronidase, from plant (dsGUS), and a vehicle control, H2O. The endpoint measurements included gene expression profiles, survival, and life history traits, such as developmental time, fecundity, hatching rate, and body length. Although S. curviseta larvae developed significantly faster under the treatments of dsDVV and dsSC than the vehicle control, the combined results from both temporal RNAi effect study and dietary RNAi toxicity assay support the risk hypothesis, suggesting that the impacts of ingested arthropod-active dsRNAs on this representative soil decomposer are negligible.

Interspecific competition between herbivores is widely recognized as an important determinant of community structure. Although researchers have identified a number of factors capable of altering competitive interactions, few studies have addressed the influence of neighboring plant species. If adaptation to/ epigenetic effects of an herbivore's natal host plant alter its performance on other host plants, then interspecific herbivore interactions may play out differently in heterogeneous and homogenous plant communities. We tested wether the natal host plant of a whitefly population affected interactions between the Middle-east Asia Minor 1 (MEAM1) and Mediterranean (MED) cryptic species of the whitefly Bemisia tabaci by rearing the offspring of a cabbage-derived MEAM1 population and a poinsettia-derived MED population together on three different host plants: cotton, poinsettia, and cabbage. We found that MED dominated on poinsettia and that MEAM1 dominated on cabbage, results consistent with previous research. MED also dominated when reared with MEAM1 on cotton, however, a result at odds with multiple otherwise-similar studies that reared both species on the same natal plant. Our work provides evidence that natal plants affect competitive interactions on another plant species, and highlights the potential importance of neighboring plant species on herbivore community composition in agricultral systems.

In China, Bemisia tabaci Q (commonly known as biotype Q) has rapidly displaced B (commonly known as biotype B) in the past 6 years. The mechanisms underlying such phenomenon have been studied extensively in recent years; however, we have not come to a definitive conclusion yet. In the present study, the differences in host suitability between B and Q whitefly adults to five host plants (cabbage, cotton, cucumber, poinsettia, and tomato) were evaluated based on their respective feeding behaviors using a direct-current electrical penetration graph (DC-EPG) system. Pair-wise comparisons of B. tabaci B and Q feeding on each of the five host plants clearly indicate that Q feeds better than B on tomato, cotton and poinsettia, while B feeds better than Q on cabbage and cucumber. The EPG parameters related to both phloem and non-phloem phases confirm that cabbage and cucumber are best suited to B, while tomato, cotton, and poinsettia are best suited to Q. Our present results support the contention that host suitability and adult feeding behavior contribute to the competitive displacement of biotype B by biotype Q. The discrepancy between field (previous studies) and laboratory results (this study), however, suggests that 1) whitefly displacement is apparently contributed by multiple factors; and 2) factor(s) other than the host plant suitability may play a vital role in dictating the whitefly biotypes in the field.

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.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a powerful technique to quantify gene expression. To facilitate gene expression study and obtain accurate results, normalization relative to stably expressed reference genes is crucial. The monarch butterfly, Danaus plexippus (L.), is one of the most recognized insect species for its spectacular annual migration across North America. Besides its great voyages, D. plexippus has drawn attention to its role as a bio-indicator, ranging from genetically modified organisms (GMOs) to natural ecosystems. In this study, nine reference genes from D. plexippus genome were selected as the candidate reference genes. The expression profiles of these candidates under various biotic and abiotic conditions were evaluated using the four readily available computational programs, BestKeeper, Normfinder, geNorm, and ΔCt method, respectively. Moreover, RefFinder, a web-based computational platform integrating the four above mentioned algorisms, provided a comprehensive ranking of the stability of these reference genes. As a result, a suite of reference genes were recommended for each experimental condition. Specifically, elongation factor 1α (EF1A) and ribosomal protein 49 (RP49) were the most stable reference genes, respectively, under biotic (development, tissue, and sex) and abiotic (photoperiod, temperature, and dietary RNAi) conditions. With the recent release of a 273-million base pair draft genome, results from this study allow us to establish a standardized RT-qPCR analysis and lay a foundation for the subsequent genomic and functional genomic research in D. plexippus, a major bio-indicator and an emerging model for migratory animals.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable, rapid, and reproducible technique for measuring and evaluating changes in gene expression. To facilitate gene expression studies and obtain more accurate RT-qPCR data, normalization relative to stable reference genes is required. In this study, expression profiles of seven candidate reference genes, including β-actin (Actin), elongation factor 1 α (EF1A), glyceralde hyde-3-phosphate dehydro-genase (GAPDH), cyclophilins A (CypA), vacuolar-type H+-ATPase (ATPase), 28S ribosomal RNA (28S), and 18S ribosomal RNA (18S) from Hippodamia convergens were investigated. H. convergens is an abundant predatory species in the New World, and has been widely used as a biological control agent against sap-sucking insect pests, primarily aphids. A total of four analytical methods, geNorm, Normfinder, BestKeeper, and the ΔCt method, were employed to evaluate the performance of these seven genes as endogenous controls under diverse experimental conditions. Additionally, RefFinder, a comprehensive evaluation platform integrating the four above mentioned algorithms, ranked the overall stability of these candidate genes. A suite of reference genes were specifically recommended for each experimental condition. Among them, 28S, EF1A, and CypA were the best reference genes across different development stages; GAPDH, 28S, and CypA were most stable in different tissues. GAPDH and CypA were most stable in female and male adults and photoperiod conditions, 28S and EF1A were most stable under a range of temperatures, Actin and CypA were most stable under dietary RNAi condition. This work establishes a standardized RT-qPCR analysis in H. convergens. Additionally, this study lays a foundation for functional genomics research in H. convergens and sheds light on the ecological risk assessment of RNAi-based biopesticides on this non-target biological control agent.

Quantitative real-time PCR (qRT-PCR) is a reliable and reproducible technique for measuring mRNA expression. To facilitate gene expression studies and obtain more accurate qRT-PCR analysis, normalization relative to stable housekeeping genes is mandatory. In this study, ten housekeeping genes, including beta-actin (Actin) , elongation factor 1 α (EF1A) , glyceralde hyde-3-phosphate dehydrogenase (GAPDH) , ribosomal protein L13 (RPL13) , ribosomal protein 49 (RP49) , α-tubulin (Tubulin) , vacuolar-type H+-ATPase (v-ATPase) , succinate dehydrogenase subunit A (SDHA) , 28S ribosomal RNA (28S) , and 18S ribosomal RNA (18S) from the two-spotted spider mite, Tetranychus urticae, were selected as the candidate reference genes. Four algorithms, geNorm, Normfinder, BestKeeper, and the ΔCt method, were used to evaluate the performance of these candidates as endogenous controls across different developmental stages. In addition, RefFinder, which integrates the above-mentioned software tools, provided the overall ranking of the stability/suitability of these candidate reference genes. Among them, PRL13 and v-ATPase were the two most stable housekeeping genes across different developmental stages. This work is the first step toward establishing a standardized qRT-PCR analysis in T. urticae following the MIQE guideline. With the recent release of the T. urticae genome, results from this study provide a critical piece for the subsequent genomics and functional genomics research in this emerging model system.

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for measuring and evaluating gene expression during variable biological processes. To facilitate gene expression studies, normalization of genes of interest relative to stable reference genes is crucial. The western flower thrips Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae), the main vector of tomato spotted wilt virus (TSWV), is a destructive invasive species. In this study, the expression profiles of 11 candidate reference genes from nonviruliferous and viruliferous F. occidentalis were investigated. Five distinct algorithms, geNorm, NormFinder, BestKeeper, the ΔCt method, and RefFinder, were used to determine the performance of these genes. geNorm, NormFinder, BestKeeper, and RefFinder identified heat shock protein 70 (HSP70), heat shock protein 60 (HSP60), elongation factor 1 α, and ribosomal protein l32 (RPL32) as the most stable reference genes, and the ΔCt method identified HSP60, HSP70, RPL32, and heat shock protein 90 as the most stable reference genes. Additionally, two reference genes were sufficient for reliable normalization in nonviruliferous and viruliferous F. occidentalis. This work provides a foundation for investigating the molecular mechanisms of TSWV and F. occidentalis interactions.

To facilitate gene expression study and obtain accurate qRT-PCR analysis, normalization relative to stable expressed housekeeping genes is required. In this study, expression profiles of 11 candidate reference genes, including actin (Actin), elongation factor 1 α (EF1A), TATA-box-binding protein (TATA), ribosomal protein L12 (RPL12), β-tubulin (Tubulin), NADH dehydrogenase (NADH), vacuolar-type H+-ATPase (v-ATPase), succinate dehydrogenase B (SDHB), 28S ribosomal RNA (28S), 16S ribosomal RNA (16S), and 18S ribosomal RNA (18S) from the pea aphid Acyrthosiphon pisum, under different developmental stages and temperature conditions, were investigated. A total of four analytical tools, geNorm, Normfinder, BestKeeper, and the ΔCt method, were used to evaluate the suitability of these genes as endogenous controls. According to RefFinder, a web-based software tool which integrates all four above-mentioned algorithms to compare and rank the reference genes, SDHB, 16S, and NADH were the three most stable house-keeping genes under different developmental stages and temperatures. This work is intended to establish a standardized qRT-PCR protocol in pea aphid and serves as a starting point for the genomics and functional genomics research in this emerging insect model.

BackgroundAccurate evaluation of gene expression requires normalization relative to the expression of reliable reference genes. Expression levels of "classical" reference genes can differ, however, across experimental conditions. Although quantitative real-time PCR (qRT-PCR) has been used extensively to decipher gene function in the sweetpotato whitefly Bemisia tabaci, a world-wide pest in many agricultural systems, the stability of its reference genes has rarely been validated.ResultsIn this study, 15 candidate reference genes from B. tabaci were evaluated using two Excel-based algorithms geNorm and Normfinder under a diverse set of biotic and abiotic conditions. At least two reference genes were selected to normalize gene expressions in B. tabaci under experimental conditions. Specifically, for biotic conditions including host plant, acquisition of a plant virus, developmental stage, tissue (body region of the adult), and whitefly biotype, ribosomal protein L29 was the most stable reference gene. In contrast, the expression of elongation factor 1 alpha, peptidylprolyl isomerase A, NADH dehydrogenase, succinate dehydrogenase complex subunit A and heat shock protein 40 were consistently stable across various abiotic conditions including photoperiod, temperature, and insecticide susceptibility.ConclusionOur finding is the first step toward establishing a standardized quantitative real-time PCR procedure following the MIQE (Minimum Information for publication of Quantitative real time PCR Experiments) guideline in an agriculturally important insect pest, and provides a solid foundation for future RNA interference based functional study in B. tabaci.

BACKGROUND: The whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), is one of the most widely distributed agricultural pests. In recent years, B. tabaci Q has invaded China, and Q has displaced B in many areas now. In a number of regions of the world, invasion by B and/or Q has been followed by outbreaks of tomato yellow leaf curl virus (TYLCV). Our previous study showed TYLCV directly and indirectly modified the feeding behavior of B. tabaci in favor of Q rather than B. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we quantified the salicylic acid (SA) titers and relative gene expression of SA in tomato leaves that were infested with viruliferous or non-viruliferous B and Q. We also measured the impacts of exogenous SA on the performance of B and Q, including the effects on ovary development. SA titer was always higher in leaves that were infested with viruliferous B than with viruliferous Q, whereas the SA titer did not differ between leaves infested with non-viruliferous B and Q. The relative gene expression of SA signaling was increased by feeding of viruliferous B but was not increased by feeding of viruliferous Q. The life history traits of B and Q were adversely affected on SA-treated plants. On SA-treated plants, both B and Q had lower fecundity, shorter longevity, longer developmental time and lower survival rate than on untreated plants. Compared with whiteflies feeding on control plants, those feeding on SA-treated plants had fewer oocytes and slower ovary development. On SA-treated plants, viruliferous B had fewer oocytes than viruliferous Q. CONCLUSIONS/SIGNIFICANCE: These results indicate that TYLCV tends to induce SA-regulated plant defense against B but SA-regulated plant defense against Q was reduced. In other words, Q may have a mutualistic relationship with TYLCV that results in the reduction of the plant's defense response.

BACKGROUND: The sweetpotato whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), is one of the most widely distributed agricultural pests. Although it has developed resistance to many registered insecticides including the neonicotinoid insecticide thiamethoxam, the mechanisms that regulate the resistance are poorly understood. To understand the molecular basis of thiamethoxam resistance, "omics" analyses were carried out to examine differences between resistant and susceptible B. tabaci at both transcriptional and translational levels. RESULTS: A total of 1,338 mRNAs and 52 proteins were differentially expressed between resistant and susceptible B. tabaci. Among them, 11 transcripts had concurrent transcription and translation profiles. KEGG analysis mapped 318 and 35 differentially expressed genes and proteins, respectively, to 160 and 59 pathways (p

While every individual of Bemisia tabaci (Hemiptera: Aleyrodidae) harbors the primary symbiont (P-symbiont) Portiera, the infection frequencies of the six secondary symbionts (S-symbionts) including Hamiltonella, Arsenophonus, Cardinium, Wolbachia, Rickettsia and Fritschea vary greatly among different populations. To characterize the factors influencing the infection dynamics of the six S-symbionts in B. tabaci, gene-specific PCR were conducted to screen for the presence of the P-symbiont Portiera and the six S-symbionts in 61 (17 B and 44 Q biotypes) field populations collected from different plant species and locations in China. All individuals of the 61 populations hosted the P-symbiont Portiera, but none of them harbored Arsenophonus and Fritschea. The presence and infection rates of Hamiltonella, Cardinium, Rickettsia, Wolbachia and their co-infections Rickettsia + Hamiltonella (RH), Rickettsia + Cardinium (RC), Hamiltonella + Cardinium (HC) and Rickettsia + Hamiltonella + Cardinium (RHC) varied significantly among the 61 field populations; and the observed variations can be explained by biotypes, sexes, host plants and geographical locations of these field populations. Taken together, at least three factors including biotype, host plant and geographical location affect the infection dynamics of S-symbionts in B. tabaci.