An Integrated Proteomics Reveals Pathological Mechanism of Honeybee (Apis cerena) Sacbrood Disease

Life-history stage determines the diet of ectoparasitic mites on their honey bee hosts

Ectoparasitic mites of the genera Varroa and Tropilaelaps have evolved to exclusively exploit honey bees as food sources during alternating dispersal and reproductive life history stages. Here we show that the primary food source utilized by Varroa destructor depends on the host life history stage. While feeding on adult bees, dispersing V. destructor feed on the abdominal membranes to access to the fat body as reported previously. However, when V. destructor feed on honey bee pupae during their reproductive stage, they primarily consume hemolymph, indicated by wound analysis, preferential transfer of biostains, and a proteomic comparison between parasite and host tissues. Biostaining and proteomic results were paralleled by corresponding findings in Tropilaelaps mercedesae, a mite that only feeds on brood and has a strongly reduced dispersal stage. Metabolomic profiling of V. destructor corroborates differences between the diet of the dispersing adults and reproductive foundresses. The proteome and metabolome differences between reproductive and dispersing V. destructor suggest that the hemolymph diet coincides with amino acid metabolism and protein synthesis in the foundresses while the metabolism of non-reproductive adults is tuned to lipid metabolism. Thus, we demonstrate within-host dietary specialization of ectoparasitic mites that coincides with life history of hosts and parasites.

Discovery of SNP Molecular Markers and Candidate Genes Associated with Sacbrood Virus Resistance in Apis cerana cerana Larvae by Whole-Genome Resequencing

Sacbrood virus (SBV) is a significant problem that impedes brood development in both eastern and western honeybees. Whole-genome sequencing has become an important tool in researching population genetic variations. Numerous studies have been conducted using multiple techniques to suppress SBV infection in honeybees, but the genetic markers and molecular mechanisms underlying SBV resistance have not been identified. To explore single nucleotide polymorphisms (SNPs), insertions, deletions (Indels), and genes at the DNA level related to SBV resistance, we conducted whole-genome resequencing on 90 Apis cerana cerana larvae raised in vitro and challenged with SBV. After filtering, a total of 337.47 gigabytes of clean data and 31,000,613 high-quality SNP loci were detected in three populations. We used ten databases to annotate 9359 predicted genes. By combining population differentiation index (FST) and nucleotide polymorphisms (π), we examined genome variants between resistant (R) and susceptible (S) larvae, focusing on site integrity (INT < 0.5) and minor allele frequency (MAF < 0.05). A selective sweep analysis with the top 1% and top 5% was used to identify significant regions. Two SNPs on the 15th chromosome with GenBank KZ288474.1_322717 (Guanine > Cytosine) and KZ288479.1_95621 (Cytosine > Thiamine) were found to be significantly associated with SBV resistance based on their associated allele frequencies after SNP validation. Each SNP was authenticated in 926 and 1022 samples, respectively. The enrichment and functional annotation pathways from significantly predicted genes to SBV resistance revealed immune response processes, signal transduction mechanisms, endocytosis, peroxisomes, phagosomes, and regulation of autophagy, which may be significant in SBV resistance. This study presents novel and useful SNP molecular markers that can be utilized as assisted molecular markers to select honeybees resistant to SBV for breeding and that can be used as a biocontrol technique to protect honeybees from SBV.

Sacbrood Virus: A Growing Threat to Honeybees and Wild Pollinators

Sacbrood virus (SBV) is one of the many viruses that infect both the Western honeybee (Apis mellifera) and the Eastern honeybee (Apis cerana). Recently, the interspecies transmission of SBV has been discovered, especially among wild pollinators. This newly discovered evolutionary occurrence regarding SBV indicates a much wider host range than previously believed, causing further concern about the future sustainability of agriculture and the resilience of ecosystems. Over the past few decades, vast numbers of studies have been undertaken concerning SBV infection in honeybees, and remarkable progress has been made in our understanding of the epidemiology, pathogenesis, transmission, and manifestations of SBV infection in honeybees and other pollinators. Meanwhile, some methods, including Chinese medicine, have been established to control and prevent sacbrood disease in A. cerana in Asian countries. In this review, we summarize the existing knowledge of SBV and address the gaps in the knowledge within the existing literature in the hope of providing future directions for the research and development of management strategies for controlling the spread of this deadly disease.

Identification of Immune Response to Sacbrood Virus Infection in Apis cerana Under Natural Condition

Chinese sacbrood virus (CSBV) is a serious threat to eastern honeybees (Apis cerana), especially larvae. However, the pathological mechanism of this deadly disease remains unclear. Here, we employed mRNA and small RNA (sRNA) transcriptome approach to investigate the microRNAs (miRNAs) and small interfering RNAs (siRNAs) expression changes of A. cerana larvae infected with CSBV under natural condition. We found that serine proteases involved in immune response were down-regulated, while the expression of siRNAs targeted to serine proteases were up-regulated. In addition, CSBV infection also affected the expression of larvae cuticle proteins such as larval cuticle proteins A1A and A3A, resulting in increased susceptibility to CSBV infection. Together, our results provide insights into sRNAs that they are likely to be involved in regulating honeybee immune response.

Genes and Proteomes Associated With Increased Mutation Frequency and Multidrug Resistance of Naturally Occurring Mismatch Repair-Deficient Salmonella Hypermutators

The emergence of antibiotic-resistant Salmonella through mutations led to mismatch repair (MMR) deficiency that represents a potential hazard to public health. Here, four representative MMR-deficient Salmonella hypermutator strains and Salmonella Typhimurium LT2 were used to comprehensively reveal the influence of MMR deficiency on antibiotic resistance among Salmonella. Our results indicated that the mutation frequency ranged from 3.39 × 10^–4 to 5.46 × 10^–2 in the hypermutator. Mutation sites in MutS, MutL, MutT, and UvrD of the four hypermutators were all located in the essential and core functional regions. Mutation frequency of the hypermutator was most highly correlated with the extent of mutation in MutS. Mutations in MMR genes (mutS, mutT, mutL, and uvrD) were correlated with increased mutation in antibiotic resistance genes, and the extent of antibiotic resistance was significantly correlated with the number of mutation sites in MutL and in ParC. The number of mutation sites in MMR genes and antibiotic resistance genes exhibited a significant positive correlation with the number of antibiotics resisted and with expression levels of mutS, mutT, and mutL. Compared to Salmonella Typhimurium LT2, a total of 137 differentially expressed and 110 specifically expressed proteins were identified in the four hypermutators. Functional enrichment analysis indicated that the proteins significantly overexpressed in the hypermutators primarily associated with translation and stress response. Interaction network analysis revealed that the ribosome pathway might be a critical factor for high mutation frequency and multidrug resistance in MMR-deficient Salmonella hypermutators. These results help elucidate the mutational dynamics that lead to hypermutation, antibiotic resistance, and activation of stress response pathways in Salmonella.

Identification of the Novel Host Protein Interacting With the Structural Protein VP1 of Chinese Sacbrood Virus by Yeast Two-Hybrid Screening

Chinese sacbrood virus (CSBV) is the major cause and lead to the collapse of Apis cerana colonies. VP1, the structural protein of CSBV, shows the highest variation in the amino acid sequences among proteins from different CSBV strains as well as exhibits excellent immunogenicity. However, its function with host protein still remains unclear. To clarify its function with host protein, we screened out host cellular proteins that interact with VP1 using the membrane protein yeast two-hybrid system. In addition, we verified interactions between heat shock protein 70 cognate 5 (Hsp70-c5) and VP1 using glutathione S-transferase (GST) pull-down and co-immunoprecipitation assays. VP1 and Hsp70-c5 were colocalized in the cytoplasm and nucleus. Using western blot and real-time polymerase chain reaction (PCR), Hsp70-c5 expression in CSBV-infected larvae was upregulated compared with that in healthy larvae. We observed that when we silenced Hsp70-c5, VP1 expression was significantly downregulated. These results demonstrate that Hsp70-c5 is involved in at least one stage(s) of the viral life cycle.

Metabolomics-based biomarker discovery for bee health monitoring: A proof of concept study concerning nutritional stress in Bombus terrestris

Bee pollinators are exposed to multiple natural and anthropogenic stressors. Understanding the effects of a single stressor in the complex environmental context of antagonistic/synergistic interactions is critical to pollinator monitoring and may serve as early warning system before a pollination crisis. This study aimed to methodically improve the diagnosis of bee stressors using a simultaneous untargeted and targeted metabolomics-based approach. Analysis of 84 Bombus terrestris hemolymph samples found 8 metabolites retained as potential biomarkers that showed excellent discrimination for nutritional stress. In parallel, 8 significantly altered metabolites, as revealed by targeted profiling, were also assigned as candidate biomarkers. Furthermore, machine learning algorithms were applied to the above-described two biomarker sets, whereby the untargeted eight components showed the best classification performance with sensitivity and specificity up to 99% and 100%, respectively. Based on pathway and biochemistry analysis, we propose that gluconeogenesis contributed significantly to blood sugar stability in bumblebees maintained on a low carbohydrate diet. Taken together, this study demonstrates that metabolomics-based biomarker discovery holds promising potential for improving bee health monitoring and to identify stressor related to energy intake and other environmental stressors.

Understanding the genetic mechanism of resistance in aphid-treated alfalfa (Medicago sativa L.) through proteomic analysis

To minimize dependency on chemical pesticides, plant breeders are trying to emphasize on important agricultural pests for the development of pest resistant cultivars. However, the molecular approach and associated genetic tools conferring resistance have not been widely studied. In the current study, proteomic analysis of two of the alfalfa cultivars viz. a resistant (R) (Zhongmu-1) and a susceptible (S) (WL343), with (+ A) and without (- A) aphids rearing were carried out. Results indicated that 325 differentially expressed proteins (DEPs) up-regulated while 319 down-regulated with a pattern of R + A/R - A plants, whereas 371 up- and 583 down-regulated DEPs were identified in the S + A/S - A plants. Total number of DEPs found in (S + A/S - A) was around 19.7% greater than that of (R + A/R - A), whereas, the down-regulated DEPs of susceptible variety was 11.6% higher than the resistant cultivar. Applying the KEGG analysis, 96 and 142 DEPs were portrayed to 15 and 10 substantively augmented pathways for Zhongmu-1 and WL343, respectively. We also found that two of the shared pathways (carbon metabolism and pyruvate metabolism) are linking to important traits conferring resistance in alfalfa. Most importantly, the specific role of linoleic acid metabolism was found to be associated with jasmonic acid, flavonoid biosynthesis, and terpenoid backbone biosynthesis that might have been associated with the insect-resistant material synthesis in the resistant alfalfa cultivar. Our study suggested that both alfalfa cultivars (R, S) could govern protein expression through carbon and pyruvate metabolism. But only the resistant alfalfa cultivar (Zhongmu-1) can tune protein expression via linoleic acid metabolism and terpenoid backbone biosynthesis to induce the defensive protein expressions (e.g., jasmonic acid and flavonoid biosynthesis along with terpenoid backbone biosynthesis), to enhance plant defense capacity.

The Wisdom of Honeybee Defenses Against Environmental Stresses

As one of the predominant pollinator, honeybees provide important ecosystem service to crops and wild plants, and generate great economic benefit for humans. Unfortunately, there is clear evidence of recent catastrophic honeybee colony failure in some areas, resulting in markedly negative environmental and economic effects. It has been demonstrated that various environmental stresses, including both abiotic and biotic stresses, functioning singly or synergistically, are the potential drivers of colony collapse. Honeybees can use many defense mechanisms to decrease the damage from environmental stress to some extent. Here, we synthesize and summarize recent advances regarding the effects of environmental stress on honeybees and the wisdom of honeybees to respond to external environmental stress. Furthermore, we provide possible future research directions about the response of honeybees to various form of stressors.

Proteomic analysis of trochophore and veliger larvae development in the small abalone Haliotis diversicolor

Background

Haliotis diversicolor is commercially important species. The trochophore and veliger are distinct larval stages in gastropod development. Their development involves complex morphological and physiological changes. We studied protein changes during the embryonic development of H. diversicolor using two dimensional electrophoresis (2-DE) and label-free methods, tandem mass spectrometry (MS/ MS), and Mascot for protein identification.

Results

A total of 150 2-DE gel spots were identified. Protein spots showed upregulation of 15 proteins and downregulation of 28 proteins as H. diversicolor developed from trochophore to veliger larvae. Trochophore and veliger larvae were compared using a label-free quantitative proteomic approach. A total of 526 proteins were identified from both samples, and 104 proteins were differentially expressed (> 1.5 fold). Compared with trochophore larvae, veliger larvae had 55 proteins upregulated and 49 proteins downregulated. These differentially expressed proteins were involved in shell formation, energy metabolism, cellular and stress response processes, protein synthesis and folding, cell cycle, and cell fate determination. Compared with the 5 protein (fructose-bisphosphate aldolase, 14–3-3ε, profilin, actin-depolymerizing factor (ADF)/cofilin) and calreticulin) expression patterns, the mRNA expression exhibited similar patterns except gene of fructose-bisphosphate aldolase.

Conclusion

Our results provide insight into novel aspects of protein function in shell formation, torsion, and nervous system development, and muscle system differentiation in H. diversicolor larvae. “Quality control” proteins were identified to be involved in abalone larval development.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4203-7) contains supplementary material, which is available to authorized users.

A comparison of biological characteristics of three strains of Chinese sacbrood virus in Apis cerana

We selected and sequenced the entire genomes of three strains of Chinese sacbrood virus (CSBV): LNQY-2008 (isolated in Qingyuan, Liaoning Province), SXYL-2015 (isolated in Yulin, Shanxi Province), and JLCBS-2014 (isolated in Changbaishan, Jilin Province), by VP1 amino acid (aa) analysis. These strains are endemic in China and infect Apis cerana. Nucleotide sequences, deduced amino acid sequences, genetic backgrounds, and other molecular biological characteristics were analysed. We also examined sensitivity of these virus strains to temperature, pH, and organic solvents, as well as to other physicochemical properties. On the basis of these observations, we compared pathogenicity and tested cross-immunogenicity and protective immunity, using antisera raised against each of the three strains. Our results showed that compared with SXYL-2015, LNQY-2008 has a 10-aa deletion and 3-aa deletion (positions 282–291 and 299–301, respectively), whereas JLCBS-2014 has a 17-aa deletion (positions 284–300). However, the three strains showed no obvious differences in physicochemical properties or pathogenicity. Moreover, there was immune cross-reactivity among the antisera raised against the different strains, implying good protective effects of such antisera. The present study should significantly advance the understanding of the pathogenesis of Chinese sacbrood disease, and offers insights into comprehensive prevention and treatment of, as well as possible protection from, the disease by means of an antiserum.

Phenotypic, histological and proteomic analyses reveal multiple differences associated with chloroplast development in yellow and variegated variants from Camellia sinensis

Leaf colour variation is observed in several plants. We obtained two types of branches with yellow and variegated leaves from Camellia sinensis. To reveal the mechanisms that underlie the leaf colour variations, combined morphological, histological, ionomic and proteomic analyses were performed using leaves from abnormal branches (variants) and normal branches (CKs). The measurement of the CIE-Lab coordinates showed that the brightness and yellowness of the variants were more intense than the CKs. When chloroplast profiles were analysed, HY1 (branch with yellow leaves) and HY2 (branch with variegated leaves) displayed abnormal chloroplast structures and a reduced number and size compared with the CKs, indicating that the abnormal chloroplast development might be tightly linked to the leaf colour variations. Moreover, the concentration of elemental minerals was different between the variants and the CKs. Furthermore, DEPs (differentially expressed proteins) were identified in the variants and the CKs by a quantitative proteomics analysis using the label-free approach. The DEPs were significantly involved in photosynthesis and included PSI, PSII, cytochrome b6/f complex, photosynthetic electron transport, LHC and F-type ATPase. Our results suggested that a decrease in the abundance of photosynthetic proteins might be associated with the changes of leaf colours in tea plants.

Proteome Analysis of the Hemolymph, Mushroom Body, and Antenna Provides Novel Insight into Honeybee Resistance against Varroa Infestation

Varroa destructor has been identified as a major culprit responsible for the losses of millions of honeybee colonies. Varroa sensitive hygiene (VSH) is a suite of behaviors from adult bees to suppress mite reproduction by uncapping and/or removing mite infested pupae from a sealed brood. Despite the efforts to elucidate the molecular underpinnings of VSH, they remain largely unknown. We investigated the proteome of mushroom bodies (MBs) and antennae of adult bees with and without VSH from a stock selected for VSH based on their response to artificially Varroa-infected brood cells by near-infrared camera observation. The pupal hemolymph proteome was also compared between the VSH-line and the line that was not selected for VSH. The identified 8609 proteins in the hemolymph, MBs, and antennae represent the most depth coverage of the honeybee proteome (>55%) to date. In the hemolymph, the VSH-line adapts a unique strategy to boost the social immunity and drive pupal organogenesis by enhancing energy metabolism and protein biosynthesis. In MBs, the up-regulated proteins implicated in neuronal sensitivity suggest their roles to promote the execution of VSH by activation of synaptic vesicles and calcium channel activities. In antennae, the highly expressed proteins associated with sensitivity of olfactory senses and signal transmissions signify their roles by inputting a strong signal to the MBs for initiating VSH. These observations illustrate that the enhanced social immunities and olfactory and neuronal sensitivity play key roles in the combat against Varroa infestation. The identified candidate markers may be useful for accelerating marker-associated selection for VSH to aid in resistance to a parasite responsible for decline in honeybee health.

Codon Optimization, Expression in Escherichia coli, and Immunogenicity of Recombinant Chinese Sacbrood Virus (CSBV) Structural Proteins VP1, VP2, and VP3

Chinese sacbrood virus (CSBV) is a small RNA virus family belonging to the genus Iflavirus that causes larval death, and even the collapse of entire bee colonies. The virus particle is spherical, non-enveloped, and its viral capsid is composed of four proteins, although the functions of the structural proteins are unclear. In this study, we used codon recoding to express the recombinant proteins VP1, VP2, and VP3 in Escherichia coli. SDS-PAGE analysis and Western blotting revealed that the target genes were expressed at high levels. Mice were then immunized with the purified, recombinant proteins, and antibody levels and lymphocyte proliferation were analyzed by ELISA and the MTT assay, respectively. The results show that the recombinant proteins induced high antibody levels and promoted lymphocyte proliferation. Polyclonal antibodies directed against these proteins will aid future studies of the molecular pathogenesis of CSBV.

In-depth phosphoproteomic analysis of royal jelly derived from western and eastern honeybee species

The proteins in royal jelly (RJ) play a pivotal role in the nutrition, immune defense, and cast determination of honeybee larvae and have a wide range of pharmacological and health-promoting functions for humans as well. Although the importance of post-translational modifications (PTMs) in protein function is known, investigation of protein phosphorylation of RJ proteins is still very limited. To this end, two complementary phosphopeptide enrichment materials (Ti(4+)-IMAC and TiO2) and high-sensitivity mass spectrometry were applied to establish a detailed phosphoproteome map and to qualitatively and quantitatively compare the phosphoproteomes of RJ produced by Apis mellifera ligustica (Aml) and Apis cerana cerana (Acc). In total, 16 phosphoproteins carrying 67 phosphorylation sites were identified in RJ derived from western bees, and nine proteins phosphorylated on 71 sites were found in RJ produced by eastern honeybees. Of which, eight phosphorylated proteins were common to both RJ samples, and the same motif ([S-x-E]) was extracted, suggesting that the function of major RJ proteins as nutrients and immune agents is evolutionary preserved in both of these honeybee species. All eight overlapping phosphoproteins showed significantly higher abundance in Acc-RJ than in Aml-RJ, and the phosphorylation of Jelleine-II (an antimicrobial peptide, TPFKLSLHL) at S(6) in Acc-RJ had stronger antimicrobial properties than that at T(1) in Aml-RJ even though the overall antimicrobial activity of Jelleine-II was found to decrease after phosphorylation. The differences in phosphosites, peptide abundance, and antimicrobial activity of the phosphorylated RJ proteins indicate that the two major honeybee species employ distinct phosphorylation strategies that align with their different biological characteristics shaped by evolution. The phosphorylation of RJ proteins are potentially driven by the activity of extracellular serine/threonine protein kinase FAM20C-like protein (FAM20C-like) through the [S-x-E] motif, which is supported by evidence that mRNA and protein expression of FAM20C-like protein kinase are both found in the highest level in the hypopharyngeal gland of nurse bees. Our data represent the first comprehensive RJ phosphorylation atlas, recording patterns of phosphorylated RJ protein abundance and antibacterial activity of some RJ proteins in two major managed honeybee species. These data constitute a firm basis for future research to better understand the biological roles of each RJ protein for honeybee biology and human health care.

Proteomic analysis of Apis cerana and Apis mellifera larvae fed with heterospecific royal jelly and by CSBV challenge

Chinese honeybee Apis cerana (Ac) is one of the major Asian honeybee species for local apiculture. However, Ac is frequently damaged by Chinese sacbrood virus (CSBV), whereas Apis mellifera (Am) is usually resistant to it. Heterospecific royal jelly (RJ) breeding in two honeybee species may result in morphological and genetic modification. Nevertheless, knowledge on the resistant mechanism of Am to this deadly disease is still unknown. In the present study, heterospecific RJ breeding was conducted to determine the effects of food change on the larval mortality after CSBV infection at early larval stage. 2-DE and MALDI-TOF/TOF MS proteomic technology was employed to unravel the molecular event of the bees under heterospecific RJ breeding and CSBV challenge. The change of Ac larval food from RJC to RJM could enhance the bee resistance to CSBV. The mortality rate of Ac larvae after CSBV infection was much higher when the larvae were fed with RJC compared with the larvae fed with RJM. There were 101 proteins with altered expressions after heterospecific RJ breeding and viral infection. In Ac larvae, 6 differential expression proteins were identified from heterospecific RJ breeding only, 21 differential expression proteins from CSBV challenge only and 7 differential expression proteins from heterospecific RJ breeding plus CSBV challenge. In Am larvae, 17 differential expression proteins were identified from heterospecific RJ breeding only, 26 differential expression proteins from CSBV challenge only and 24 differential expression proteins from heterospecific RJ breeding plus CSBV challenge. The RJM may protect Ac larvae from CSBV infection, probably by activating the genes in energy metabolism pathways, antioxidation and ubiquitin-proteasome system. The present results, for the first time, comprehensively descript the molecular events of the viral infection of Ac and Am after heterospecific RJ breeding and are potentially useful for establishing CSBV resistant populations of Ac for apiculture.

Hemolymph proteome changes during worker brood development match the biological divergences between western honey bees (Apis mellifera) and eastern honey bees (Apis cerana)

Background

Hemolymph plays key roles in honey bee molecule transport, immune defense, and in monitoring the physiological condition. There is a lack of knowledge regarding how the proteome achieves these biological missions for both the western and eastern honey bees (Apis mellifera and Apis cerana). A time-resolved proteome was compared using two-dimensional electrophoresis-based proteomics to reveal the mechanistic differences by analysis of hemolymph proteome changes between the worker bees of two bee species during the larval to pupal stages.

Results

The brood body weight of Apis mellifera was significantly heavier than that of Apis cerana at each developmental stage. Significantly, different protein expression patterns and metabolic pathways were observed in 74 proteins (166 spots) that were differentially abundant between the two bee species. The function of hemolymph in energy storage, odor communication, and antioxidation is of equal importance for the western and eastern bees, indicated by the enhanced expression of different protein species. However, stronger expression of protein folding, cytoskeletal and developmental proteins, and more highly activated energy producing pathways in western bees suggests that the different bee species have developed unique strategies to match their specific physiology using hemolymph to deliver nutrients and in immune defense.

Conclusions

Our disparate findings constitute a proof-of-concept of molecular details that the ecologically shaped different physiological conditions of different bee species match with the hemolymph proteome during the brood stage. This also provides a starting point for future research on the specific hemolymph proteins or pathways related to the differential phenotypes or physiology.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-563) contains supplementary material, which is available to authorized users.

Identification of Himetobi P virus in the small brown planthopper by deep sequencing and assembly of virus-derived small interfering RNAs

Profiling and assembly of virus-derived small interfering RNAs (siRNAs) using next-generation sequencing technologies have been very useful for identification and diagnosis of a number of plant and invertebrate viruses. In this work, we have conducted high-throughput pyrosequencing and bioinformatic analysis of the small brown planthopper (SBPH, Laodelphax striatellus), and these analyses unexpectedly showed that the Himetobi P virus (HiPV) was present in our laboratory cultures. HiPV was also found to infect our brown planthopper (BPH, Nilaparvata lugens) and the white-backed planthopper (WBPH, Sogatella furcifera) cultures. The majority of the HiPV-derived siRNAs (Hd-siRNAs) were 21 and 22 nucleotides in length and nearly two-thirds of the siRNAs originated from the HiPV genomic RNA strand. The Hd-siRNAs were evenly distributed across the genome and this indicates that the HiPV genome contributes uniformly to production of Hd-siRNAs. Although HiPV infection appeared to be innocuous to the SBPH, alterations of gene expressions involved in reproduction, cytoskeleton structure and defense responses such as RNA interference pathways (RNAi) genes were observed. Furthermore, we demonstrated that silencing Agronaute 2 in L. striatellus enhanced HiPV accumulation, and this observation provides evidence for the existence of RNAi defenses against HiPV in the SBPH.