Selection of stable reference genes for quantitative real-time PCR in the Varroa mite, Varroa destructor.

To investigate the acaricide toxicity and resistance mechanisms in the Varroa mite, it is essential to understand the genetic responses of Varroa mites to acaricides, which are usually evaluated by transcriptional profiling based on quantitative real-time polymerase chain reaction (qPCR). In this study, to select reference genes showing consistent expression patterns regardless of the acaricide treatment or the type of tissue, Varroa mites treated with each of the three representative acaricides (coumaphos, fluvalinate, and amitraz) were processed for transcriptomic analysis, from which eight genes (NADH dehydrogenase [NADHD], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], eukaryotic translation elongation factor 1 α 1 [eEF1A1], eukaryotic translation elongation factor 2 [eEF2], ribosomal protein L5 [RpL5], Actin, tubulin α-1D chain [α-tubulin], and Rab1) were selected as candidates. The transcription profiles of these genes, depending on the treatment of the three acaricides or across different tissues (cuticle, legs, gut/fat bodies, and synganglion), were analyzed using qPCR with four validation programs, BestKeeper, geNorm, NormFinder, and RefFinder. Following acaricide treatment, eEF1A1 and NADHD showed the least variation in their expression levels, whereas the expression levels of α-tubulin and RpL5 were the most stable across different tissue groups. Rab1/GAPDH and Actin/eEF2 showed the least stable expression patterns following acaricide treatments and across different tissues, respectively, requiring precautions for use. When vitellogenin gene expression was analyzed by different reference genes, its expression profiles varied significantly depending on the reference genes, highlighting the importance of proper reference gene use. Thus, it is recommended using eEF1A1 and NADHD as reference genes for the comparison of the effects of acaricide on the whole body, whereas α-tubulin and RpL5 are recommended for investigating the tissue-specific expression profiles of target genes.

Selection of lethal genes for ingestion RNA interference against western flower thrips, Frankliniella occidentalis, via leaf disc-mediated dsRNA delivery.

The western flower thrips, Frankliniella occidentalis, is a major pest that damages a wide variety of crops and vegetables. Following extensive use of insecticides, it has developed high levels of resistance to almost all groups of insecticides due to its high reproduction rate and short generation time. Therefore, an alternative pest control strategy, such as RNA interference (RNAi)-based control, is essential. To establish an ingestion RNAi-based control, a total of 57 genes involved in various biological processes were selected, and their double-stranded RNAs (dsRNA) were delivered to an insecticide-susceptible strain of F. occidentalis via the leaf disc-feeding method using a bioassay chamber optimized by 3D printing. The mortality of dsRNA-ingested thrips was examined every 24 h until 120 h post-treatment. Of the 57 genes screened, dsRNAs of the Toll-like receptor 6, apolipophorin, coatomer protein subunit epsilon and sorting and assembly machinery component were most lethal when ingested by thrips. The dsRNA-fed thrips showed substantially reduced transcription levels of target genes, demonstrating that the observed mortality was likely due to RNAi. When these genes were tested for ingestion RNAi against an insecticide-resistant strain of F. occidentalis, bioassay results were similar. In conclusion, this study provides the first evidence that ingestion RNAi can be lethal to F. occidentalis, a mesophyll sucking pest, and further suggests that transgenic plants expressing hairpin RNA of these essential genes can be employed to control insecticide-resistant thrips.

Pathway profiles based on gene-set enrichment analysis in the honey bee Apis mellifera under brood rearing-suppressed conditions.

Perturbation of normal behaviors in honey bee colonies by any external factor can immediately reduce the colony's capacity for brood rearing, which can eventually lead to colony collapse. To investigate the effects of brood-rearing suppression on the biology of honey bee workers, gene-set enrichment analysis of the transcriptomes of worker bees with or without suppressed brood rearing was performed. When brood rearing was suppressed, pathways associated with both protein degradation and synthesis were simultaneously over-represented in both nurses and foragers, and their overall pathway representation profiles resembled those of normal foragers and nurses, respectively. Thus, obstruction of normal labor induced over-representation in pathways related with reshaping of worker bee physiology, suggesting that transition of labor is physiologically reversible. In addition, some genes associated with the regulation of neuronal excitability, cellular and nutritional stress and aggressiveness were over-expressed under brood rearing suppression perhaps to manage in-hive stress under unfavorable conditions.

Transcriptome-based identification and characterization of genes commonly responding to five different insecticides in the diamondback moth, Plutella xylostella.

When the 3rd instar larvae of the diamondback moth (DBM), Plutella xylostella, were pretreated with sublethal doses (LC10) and then subsequently exposed to lethal doses (LC50) of chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad via leaf dipping, their tolerance to insecticides was significantly enhanced. To identify genes that commonly respond to the treatment of different insecticides and are responsible for the tolerance enhancement, transcriptomic profiles of larvae treated with sublethal doses of the five insecticides were compared with that of untreated control. A total of 117,181 transcripts with a mean length of 662bp were generated by de novo assembly, of which 35,329 transcripts were annotated. Among them, 125, 143, 182, 215 and 149 transcripts were determined to be up-regulated whereas 67, 45, 60, 60 and 38 genes were down-regulated following treatments with chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad, respectively. Gene ontology (GO) analysis of differentially expressed genes (DEGs) revealed little differences in their GO profiles between treatments with different insecticides except for spinosad. Finally, the DEGs commonly responding to all insecticides were selected for further characterization, and some of their over-transcription levels were confirmed by quantitative PCR. The most notable examples of commonly responding over-transcribed genes were two cytochrome P450 genes (Cyp301a1 and Cyp9e2) and nine cuticular protein genes. In contrast, several genes composing the mitochondrial energy generation system were significantly down-regulated in all treated larvae. Considering the distinct structure and mode of action of the five insecticides tested, the differentially expressed genes identified in this study appear to be involved in general chemical defense at the initial stage of intoxication. Their possible roles in the tolerance/resistance development were discussed.

Characterization of Venom Components and Their Phylogenetic Properties in Some Aculeate Bumblebees and Wasps.

: To identify and compare venom components and expression patterns, venom gland-specific transcriptome analyses were conducted for 14 Aculeate bees and wasps. TPM (transcripts per kilobase million) values were normalized using the average transcription level of a reference housekeeping gene (dimethyladenosine transferase). Orthologous venom component genes across the 14 bee and wasp species were identified, and their relative abundance in each species was determined by comparing normalized TPM values. Based on signal sequences in the transcripts, the genes of novel venom components were identified and characterized to encode potential allergens. Most of the allergens and pain-producing factors (arginine kinase, hyaluronidase, mastoparan, phospholipase A1, phospholipase A2, and venom allergen 5) showed extremely high expression levels in social wasps. Acid phosphatase, neprilysin, and tachykinin, which are known allergens and neurotoxic peptides, were found in the venom glands of solitary wasps more often than in social wasps. In the venom glands of bumblebees, few or no transcripts of major allergens or pain-producing factors were identified. Taken together, these results indicate that differential expression patterns of the venom genes in some Aculeate species imply that some wasps and bumblebee species have unique groups of highly expressed venom components. Some venom components reflected the Aculeate species phylogeny, but others did not. This unique evolution of specific venom components in different groups of some wasps and bumblebee species might have been shaped in response to both ecological and behavioral influences.

Differential expression of acetylcholinesterase 1 in response to various stress factors in honey bee workers.

The honey bee acetylcholinesterase 1 (AmAChE1) has been suggested to be related to stress response as judged from its elevated expression level under brood rearing-suppressed conditions. To further investigate the involvement of AmAChE1 expression in the stress response and its physiological functions, we analyzed altered expression profiles of AmAChE1 induced by diverse stress factors. In addition, transcription profiles of several heat shock protein (Hsp) genes (hsps) and the vitellogenin (Vg) gene (vg) known as general stress markers were investigated as positive references. Among the tested stress conditions, AmAChE1 expression was induced under the brood rearing-suppressed, crowding and heat shock conditions. The hsps, particularly hsp70 and hsp90, responded to seven of nine stress conditions tested, confirming that hsp expression profiles can serve as a general stress marker. Taken together, AmAChE1 expression is not suitable for using as a stress marker due to its limited response. Nevertheless, AmAChE1 expression appears to be connected, at least in part, to heat shock response and other pathways. Considering that AmAChE1 likely regulates the ACh titer particularly in non-neuronal tissues, thereby modulating the signal cascades mediated by mAChR, the AmAChE1 expression profile under different conditions likely provides important information on its physiological roles in honey bees.

Expression of acetylcholinesterase 1 is associated with brood rearing status in the honey bee, Apis mellifera.

Acetylcholinesterase 1 (AmAChE1) of the honey bee, Apis mellifera, has been suggested to have non-neuronal functions. A systematic expression profiling of AmAChE1 over a year-long cycle on a monthly basis revealed that AmAChE1 was predominantly expressed in both head and abdomen during the winter months and was moderately expressed during the rainy summer months. Interestingly, AmAChE1 expression was inhibited when bees were stimulated for brood rearing by placing overwintering beehives in strawberry greenhouses with a pollen diet, whereas it resumed when the beehives were moved back to the cold field, thereby suppressing brood rearing. In early spring, pollen diet supplementation accelerated the induction of brood-rearing activity and the inhibition of AmAChE1 expression. When active beehives were placed in a screen tent in late spring, thereby artificially suppressing brood-rearing activity, AmAChE1 was highly expressed. In contrast, AmAChE1 expression was inhibited when beehives were allowed to restore brood rearing by removing the screen, supporting the hypothesis that brood rearing status is a main factor in the regulation of AmAChE1 expression. Since brood rearing status is influenced by various stress factors, including temperature and diet shortage, our finding discreetly suggests that AmAChE1 is likely involved in the stress response or stress management.