The miR-9b microRNA mediates dimorphism and development of wing in aphids.

Wing dimorphism is a phenomenon of phenotypic plasticity in aphid dispersal. However, the signal transduction for perceiving environmental cues (e.g., crowding) and the regulation mechanism remain elusive. Here, we found that aci-miR-9b was the only down-regulated microRNA (miRNA) in both crowding-induced wing dimorphism and during wing development in the brown citrus aphid Aphis citricidus We determined a targeted regulatory relationship between aci-miR-9b and an ABC transporter (AcABCG4). Inhibition of aci-miR-9b increased the proportion of winged offspring under normal conditions. Overexpression of aci-miR-9b resulted in decline of the proportion of winged offspring under crowding conditions. In addition, overexpression of aci-miR-9b also resulted in malformed wings during wing development. This role of aci-miR-9b mediating wing dimorphism and development was also confirmed in the pea aphid Acyrthosiphon pisum The downstream action of aci-miR-9b-AcABCG4 was based on the interaction with the insulin and insulin-like signaling pathway. A model for aphid wing dimorphism and development was demonstrated as the following: maternal aphids experience crowding, which results in the decrease of aci-miR-9b. This is followed by the increase of ABCG4, which then activates the insulin and insulin-like signaling pathway, thereby causing a high proportion of winged offspring. Later, the same cascade, "miR-9b-ABCG4-insulin signaling," is again involved in wing development. Taken together, our results reveal that a signal transduction cascade mediates both wing dimorphism and development in aphids via miRNA. These findings would be useful in developing potential strategies for blocking the aphid dispersal and reducing viral transmission.

Fusion dsRNA in targeting salivary protein genes enhance the RNAi-based aphid control.

RNA interference (RNAi) is a promising tool for pest control and relies on sequence-specific gene silencing. Salivary proteins are cooperatively secreted into plants to guarantee the feeding of aphids; thus they have potential to develop as selective targets for RNAi-based pest control strategy. For this purpose, we firstly analyzed 18 salivary proteomes of various aphid species, and these salivary proteins can be mainly categorized into seven functional groups. Secondly, we created a work-flow for fusion dsRNA design that can target multiple genes but were selectively safe to beneficial insects. Based on this approach, seven fusion dsRNAs were designed to feed the green peach aphid, which induced a significant reduction in aphid fitness. Among them, ingestion of dsperoxidase induced the highest mortality in aphids, which was also significantly higher than that of traditional dsRNAs in targeting three peroxidases separately. In addition, dsperoxidase-fed green peach aphids triggered the highest H2O2 content of host plants as well as the attraction to natural enemies (ladybeetle and parasitic wasp) but repellent to other control aphids. Our results indicate that the fusion dsRNA design approach can improve aphid control capacity, and the fusion dsRNA targeting salivary protein-encoding genes can enhance the direct and indirect defenses of host plants, thus providing a new strategy for RNAi-based aphid control.

The Diversity of Viral Community in Invasive Fruit Flies (Bactrocera and Zeugodacus) Revealed by Meta-transcriptomics.

RNA viruses are extremely diverse and rapidly evolving in various organisms. Our knowledge on viral evolution with interacted hosts in the manner of ecology is still limited. In the agricultural ecosystem, invasive insect species are posing a great threat to sustainable crop production. Among them, fruit flies (Diptera: Tephritidae Bactrocera and Zeugodacus) are destructive to fruits and vegetables, which are also closely related and often share similar ecological niches. Thus, they are ideal models for investigating RNA virome dynamics in host species. Using meta-transcriptomics, we found 39 viral sequences in samples from 12 fly species. These viral species represented the diversity of the viromes including Dicistroviridae, negev-like virus clades, Thika virus clades, Solemoviridae, Narnaviridae, Nodaviridae, Iflaviridae, Orthomyxoviridae, Bunyavirales, Partitiviridae, and Reoviridae. In particular, dicistrovirus, negev-like virus, orthomyxovirus, and orbivirus were common in over four of the fly species, which suggests a positive interaction between fly viromes that exist under the same ecological conditions. For most of the viruses, the virus-derived small RNAs displayed significantly high peaks in 21 nt and were symmetrically distributed throughout the viral genome. These results suggest that infection by these viruses can activate the host's RNAi immunity. Our study provides RNA virome diversity and evidence on their infection activity in ecologically associated invasive fruit fly species, which could help our understanding of interactions between complex species and viruses.

Prevalence of a Novel Bunyavirus in Tea Tussock Moth Euproctis pseudoconspersa (Lepidoptera: Lymantriidae).

Euproctis pseudoconspersa is a major pest of tea plants, and also causes a skin rash on workers in tea plantations. Research on virus could provide fundamental insights for classification, genetic diversity, evolution, and host-virus interaction mechanisms. Here, we identified a novel RNA virus, Euproctis pseudoconspersa bunyavirus (Phenuiviridae), and found that it is widely distributed in field populations of E. pseudoconspersa. The replication of virus in E. pseudoconspersa was indicated by Tag-PCR. These results contribute to the classification of bunyaviruses and provide insight into the diversity of commensal E. pseudoconspersa bunyavirus and the host.

Genome-Wide Characterization and Identification of Long Non-Coding RNAs during the Molting Process of a Spider Mite, Panonychus citri.

Molting is essential for arthropods to grow. As one of the important arthropod pests in agriculture, key spider mite species (Tetranychus and Panonychus) can normally molt three times from the larva to adult stage within a week. This physiological strategy results in the short lifecycle of spider mites and difficulties in their control in the field. Long non-coding RNAs (lncRNAs) regulate transcriptional editing, cellular function, and biological processes. Thus, analysis of the lncRNAs in the spider mite molting process may provide new insights into their roles in the molting mechanism. For this purpose, we used high-throughput RNA-seq to examine the expression dynamics of lncRNAs and mRNAs in the molting process of different development stages in Panonychus citri. We identified 9199 lncRNAs from 18 transcriptomes. Analysis of the lncRNAs suggested that they were shorter and had fewer exons and transcripts than mRNAs. Among these, 356 lncRNAs were differentially expressed during three molting processes: late larva to early protonymph, late protonymph to early deutonymph, and late deutonymph to early adult. A time series profile analysis of differentially expressed lncRNAs showed that 77 lncRNAs were clustered into two dynamic expression profiles (Pattern a and Pattern c), implying that lncRNAs were involved in the molting process of spider mites. Furthermore, the lncRNA-mRNA co-expression networks showed that several differentially expressed hub lncRNAs were predicted to be functionally associated with typical molting-related proteins, such as cuticle protein and chitin biosynthesis. These data reveal the potential regulatory function of lncRNAs in the molting process and provide datasets for further analysis of lncRNAs and mRNAs in spider mites.

The RNA Virome and Its Dynamics in an Invasive Fruit Fly, Bactrocera dorsalis, Imply Interactions Between Host and Viruses.

The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important invasive agricultural insect pest with a wide host range, and has spread around the world over the last century. This evolutionary trait may have arisen primarily from interactions between B. dorsalis and other invertebrates that share the same ecological niches. The invasive behavior of B. dorsalis also frequently exposes them to diverse species of viruses. Thereby, RNA viromes may be useful microbial markers to understand the ecological evolution of B. dorsalis as well as to investigate virus-host interactions. Here, we reported eight novel RNA viruses in B. dorsalis of a lab colony, including four positive-strand RNA viruses, two negative-strand RNA viruses, and two double-stranded RNA viruses using high-throughput sequencing technology. Analysis of the virus-derived small RNAs suggested that most of these viruses may be active and trigger the host antiviral RNAi responses. The viruses were also detected in various geographical populations of B. dorsalis, implying that there is a strong association between the viromes and host. In addition, these viruses infected specific fly tissues, predominately the central nervous system and gut. Furthermore, we explored the dynamics of the viruses when hosts were exposed to short- or long-term stressors, which showed that titers of some viruses were responsively altered in the stressed B. dorsalis. The discovery of these viruses may enrich our understanding of the species diversity of RNA viruses and also provide information on viruses in association with host adaptation in insects.

RNA virome screening in diverse but ecologically related citrus pests reveals potential virus-host interactions.

As an evergreen ecosystem, citrus orchards have specialized pest species and stable ecological homeostasis; thus, they provide an ideal model for investigating RNA viromes in diverse but ecologically related species. For this purpose, we collected specialized citrus pests from three classes of invertebrates, Insecta, Arachnida, and Gastropoda and we constructed two kinds of libraries (RNA and small RNA) for the pests by deep sequencing. In total, six virus-derived sequences were identified, including four Picornavirales, one Jingchuvirales and one Nidovirales. The picornavirus-derived small RNAs showed significant small RNA peaks and symmetric distribution patterns along the genome, which suggests these viruses infected the hosts and triggered host antiviral immunity RNA interference. Screening of virus-derived sequences in multiple species of citrus pests (n = 10 per species) showed that Eotetranychus kankitus picorna-like virus and Tetranychus urticae mivirus may be present in multiple pests. Our investigation in citrus pests confirmed that RNA viruses revealed by metagenomics could impact host immunity (e.g. RNAi). An approach with parallel deep sequencing of RNAs and small RNAs is useful not only for viral discoveries but also for understanding virus-host interactions of ecologically related but divergent pest species.

First Report of a Mesonivirus and Its Derived Small RNAs in an Aphid Species Aphis citricidus (Hemiptera: Aphididae), Implying Viral Infection Activity.

We report a new positive-sense single-stranded RNA (ss RNA+) virus from the brown citrus aphid Aphis citricidus. The 20,300 nucleotide (nt)-long viral genome contains five open-reading frames and encodes six conserved domains (TM2, 3CLpro, TM3, RdRp, Zm, and HEL1). Phylogenetic analysis and amino acid sequence analysis revealed this virus might belong to an unassigned genus in the family Mesoniviridae. The presence of the virus was also confirmed in the field population. Importantly, analysis of the virus-derived small RNAs showed a 22-nt peak, implying that viral infection triggers the small interfering RNA pathway as antiviral immunity in aphids. This is the first report of a mesonivirus in invertebrates other than mosquitoes.

RNAi of the nuclear receptor HR3 suggests a role in the molting process of the spider mite Panonychus citri.

Ecdysteroids regulate molting in arthropods by binding to heterodimers of the ecdysone receptor and retinoid-X-receptor, homologous to the ultraspiracle protein, to induce the expression of downstream signal response genes including the nuclear receptor HR3. However, the detailed expression dynamics of HR3 during molting in spider mites are not yet clear. In this study, the full length of PcHR3 was retrieved based on the genome of citrus red mite, Panonychus citri. The open reading frame is 1707 bp encoding 568 amino acids, which contains a DNA binding domain and a ligand binding domain. Then, the expression pattern of PcHR3 was analyzed throughout the development of the deutonymph by RT-qPCR. The result showed that PcHR3 was mainly transcribed in the late deutonymph stage, when the deutonymph was at least 24 h old and motionless, the critical point at which the mites started molting. Transcription reached the highest level in 32-h-old deutonymphs and decreased by 36 h, where the mites remained in a quiescent state. Further silencing of PcHR3 by leaf-disc-based delivery of dsRNA to 8-h-old deutonymph mites, resulted in retarded development and death of 58% of deutonymphs. In summary, we suggest that PcHR3 regulates the latter stages of molting in P. citri.

Characterization of a new bunyavirus and its derived small RNAs in the brown citrus aphid, Aphis citricidus.

High-throughput sequencing is widely used for virus discovery, and many RNA viruses have been discovered and identified. A new negative-sense single-stranded RNA virus was identified in the brown citrus aphid and named Aphis citricidus bunyavirus. The genome consists of large (7037 nt), medium (3462 nt), and small (1163 nt) segments. Phylogenetic analysis and amino acid sequences identities of this virus with other bunyaviruses suggest that it is a new species belonging to the family Phenuiviridae. The small interfering RNA pathway could be involved against the infection of this virus in brown citrus aphid as supported by the viral derived small RNAs. The discovery of this virus illustrates the diversity of RNA viruses and contributes to the classification of bunyaviruses.

Three novel RNA viruses in the spider mite Tetranychus urticae and their possible interactions with the host RNA interference response.

As a polyphagous herbivore, the two-spotted spider mite Tetranychus urticae is engaged with various plant hosts and interacts with diverse organisms that share the same ecological niche. Thus, T. urticae faces frequent challenges from viral infections. However, the RNA viruses of T. urticae are still unknown. Here, we constructed two libraries (~8 Gb for RNA and ~10 Mb for small RNA) from a strain of T. urticae using deep sequencing, and identified three novel RNA viruses from the families Kitaviridae, Dicistroviridae, and Chuviridae. Among them, the Kitaviridae and Dicistroviridae viruses presented a possible interaction pattern with the host RNA interference pathway.

Expression dynamics of key ecdysteroid and juvenile hormone biosynthesis genes imply a coordinated regulation pattern in the molting process of a spider mite, Tetranychus urticae.

In insects, the ecdysteroid 20-hydroxyecdysone coordinates with juvenile hormone (JH) to regulate the process of molting, development and metamorphosis; however, this interaction is still unclear in the mites. In this study, we investigated the gene related to ecdysteroid and JH biosynthesis pathways, including four ecdysteroid and 11 JH biosynthesis genes. We examined their expression patterns during molting of different developmental stages of the two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae), an important agricultural pest that feeds on more than 1100 plant species. The expression of ecdysteroid biosynthesis Halloween genes exhibited a positive zigzag-like pattern, with a peak after 8 h of molting and a drop 8 h after entering each quiescent stage. In contrast, JH biosynthesis genes expression displayed a negative zigzag-like pattern, with a peak at 8 h after entering each quiescent stage and a drop after 8 h of each molting. These opposite patterns imply that ecdysteroid and JH expression is coordinated during the developmental transition. Our data provide an initial perspective on the co-expression of ecdysteroid and JH biosynthesis genes to regulate this important developmental process in the two-spotted spider mite.

Silencing of Two Insulin Receptor Genes Disrupts Nymph-Adult Transition of Alate Brown Citrus Aphid.

Insulin receptors play key roles in growth, development, and polymorphism in insects. Here, we report two insulin receptor genes (AcInR1 and AcInR2) from the brown citrus aphid, Aphis (Toxoptera) citricidus. Transcriptional analyses showed that AcInR1 increased during the nymph-adult transition in alate aphids, while AcInR2 had the highest expression level in second instar nymphs. AcInR1 is important in aphid development from fourth instar nymphs to adults as verified by dsRNA feeding mediated RNAi. The silencing of AcInR1 or/and AcInR2 produced a variety of phenotypes including adults with normal wings, malformed wings, under-developed wings, and aphids failing to develop beyond the nymphal stages. Silencing of AcInR1 or AcInR2 alone, and co-silencing of both genes, resulted in 73% or 60%, and 87% of aphids with problems in the transition from nymph to normal adult. The co-silencing of AcInR1 and AcInR2 resulted in 62% dead nymphs, but no mortality occurred by silencing of AcInR1 or AcInR2 alone. Phenotypes of adults in the dsInR1 and dsInR2 were similar. The results demonstrate that AcInR1 and AcInR2 are essential for successful nymph-adult transition in alate aphids and show that RNAi methods may be useful for the management of this pest.

Abamectin treatment affects glutamate decarboxylase expression and induces higher GABA levels in the citrus red mite, Panonychus citri.

The citrus red mite, Panonychus citri, is one of the most economically and globally destructive mite pests of citrus. Acaricide resistance has been a growing problem in controlling this pest. As the main inhibitory neurotransmitter in organisms, γ-aminobutyric acid (GABA) is synthesized from the amino acid glutamate by the action of glutamate decarboxylases (GADs). In the present study, one novel GAD gene, PcGAD, was identified and characterized from P. citri. The opening reading frame of PcGAD contained 1548 nucleotides that encode 515 amino acids. The subsequent spatiotemporal expression pattern by RT-qPCR revealed that the expression levels of PcGAD were significantly higher in larvae than in adults. Challenging with various concentrations of abamectin resulted in the upregulation of PcGAD transcript levels. Furthermore, biochemical characterization indicated that changes in GAD activity coincided with its mRNA levels. High-performance liquid chromatography confirmed that the GABA contents of P. citri increased upon abamectin treatment. The application of abamectin induces PcGAD expression and activates GAD activity, thereby resulting in an increase in GABA content in P. citri, which contributes to the adaptability of the mite to abamectin challenge.

Inheritance, Realized Heritability, and Biochemical Mechanisms of Malathion Resistance in Bactrocera dorsalis (Diptera: Tephritidae).

To better characterize the resistance development and therefore establish effective pest management strategies, this study was undertaken to investigate the inheritance mode and biochemical mechanisms of malathion resistance in the oriental fruit fly, Bactrocera dorsalis (Hendel), which is one of the most notorious pests in the world. After 22 generations of selection with malathion, the malathion-resistant (MR) strain of B. dorsalis developed a 34-fold resistance compared with a laboratory susceptible strain [malathion-susceptible (MS)]. Bioassay results showed that there was no significant difference between the LD50 values of malathion against the progenies from both reciprocal crosses (F(1)-SR and F(1)-RS). The degree of dominance values (D) was calculated as 0.39 and 0.32 for F(1)-RS and F(1)-SR, respectively. The logarithm dosage-probit mortality lines of the F(2) generation and progeny from the backcross showed no clear plateaus of mortality across a range of doses. In addition, Chi-square analysis revealed significant differences between the mortality data and the theoretical expectations. The realized heritability (h(2)) value was 0.16 in the laboratory-selected resistant strain of B. dorsalis. Enzymatic activities identified significant changes of carboxylesterases, cytochrome P450 (general oxidases), and glutathione S-transferases in MR compared with the MS strain of B. dorsalis. Taken together, this study revealed for the first time that malathion resistance in B. dorsalis follows an autosomal, incompletely dominant, and polygenic mode of inheritance and is closely associated with significantly elevated activities of three major detoxification enzymes.

Functional analysis of a chitinase gene during the larval-nymph transition in Panonychus citri by RNA interference.

Chitinases are hydrolytic enzymes that are required for chitin degradation and reconstruction in arthropods. In this study, we report a cDNA sequence encoding a putative chitinase (PcCht1) from the citrus red mite, Panonychus citri. The PcCht1 (564 aa) possessed a signal peptide, a conserver domain, and a chitin-binding domain. Structural and phylogenetic analyses found that PcCht1 had high sequence similarity to chitinases in Tetranychus urticae. Real-time quantitative PCR analyses showed that the transcript levels of PcCht1 peaked periodically in larval and nymph stages. Moreover, significant increase of PcCht1 transcript level in the larvae was observed upon the exposure of diflubenzuron. In contrast, exposures of the larvae to diflubenzuron resulted in the decreased chitin content. Furthermore, through a feeding-based RNA interference approach, we were able to reduce the PcCht1 transcript level by 59.7 % in the larvae, and consequently the treated larvae showed a very low molting rate compared with the control. Our results expanded the understanding of the important role of PcCht1 in the growth and development of P. citri.

Molecular characterization of vitellogenin and its receptor genes from citrus red mite, Panonychus citri (McGregor).

The production and uptake of yolk protein play an important role in the reproduction of all oviparous organisms. Vitellogenin (Vg) is the precursor of vitellin (Vn), which is the major egg storage protein, and vitellogenin receptor (VgR) is a necessary protein for the uptake of Vg into developing oocytes. In this paper, we characterize the full-length Vg and VgR, PcVg1 and PcVgR, respectively, of the citrus red mite Panonychus citri (McGregor). The PcVg1 cDNA is 5748 nucleotides (nt) with a 5553-nt open reading frame (ORF) coding for 1851 amino acids (aa), and the PcVgR is 6090 nt, containing an intact ORF of 5673 nt coding an expected protein of 1891 aa. The PcVg1 aa sequence shows a typical GLCG domain and several K/RXXR cleavage sites, and PcVgR comprises two ligand-binding domains, two epidermal growth factor (EGF)-like regions containing YWTD motifs, a transmembrane domain, and a cytoplasmic domain. An analysis of the aa sequences and phylogenetics implied that both genes were genetically distinct from those of ticks and insects. The transcriptional profiles determined by real-time quantitative PCR in different developmental stages showed that both genes present the same expressional tendencies in eggs, larvae, nymphs, and adults. This suggested that the biosynthesis and uptake of PcVg occurs coordinately. The strong reproductive capacity of P. citri has been hypothesized as an important factor in its resistance; consequently, understanding the molecular mechanisms regulating Vg and VgR are fundamental for mite control.

The Essential Role of Vitellogenin Receptor in Ovary Development and Vitellogenin Uptake in Bactrocera dorsalis (Hendel).

The vitellogenin receptor (VgR) functions as an essential component in uptaking and transporting vitellogenin (Vg) in female adults, which is involved in ovary development and oviposition. This study aimed to clarify the molecular characteristics and function of VgR in the oriental fruit fly Bactrocera dorsalis (Hendel). Here, we identified the full-length of BdVgR (GenBank Accession No. JX469118), encoding a 1925 residue (aa) protein with a 214.72 kDa molecular mass and several typical motifs of low-density lipoprotein receptor superfamily (LDLR). Phylogenic analysis suggested that BdVgR was evolutionary conserved with other Dipteran VgRs. The expression of BdVgR was exclusively detected in the ovaries rather than head, thorax or other tissues. The developmental expression patterns showed that the signal of BdVgR was detectable in very beginning of adult stage, and positively correlated with the growth rate of ovaries and the expression levels of its ligands. In addition, we also demonstrated that the expression level of BdVgR, and ovary development were significantly suppressed after being injected with BdVgR-targeted dsRNA. Together, all of these results indicated that BdVgR was critical for yolk protein absorption and ovary maturation in B. dorsalis, playing a vital role in female reproduction.

Analysis of reference gene stability after Israeli acute paralysis virus infection in bumblebees Bombus terrestris.

To date, there are no validated internal reference genes for the normalization of RT-qPCR data from virus infection experiments with pollinating insects. In this study we evaluated the stability of five candidate internal reference genes: elongation factor-1-alpha (ELF1α), peptidylprolyl isomerase A (PPIA), 60S ribosomal protein L23 (RPL23), TATA-binding protein (TBP) and polyubiquitin (UBI), in relation to Israeli acute paralysis virus (IAPV) infection of Bombus terrestris. We investigated the stability of these genes: in whole bodies and individual body parts, as well as in whole bodies collected at different time intervals after infection with IAPV. Our data identified PPIA as the single, most-optimal internal reference gene and the combination of PPAI-RPL23-UBI as a fully-sufficient multiple internal reference genes set for IAPV infection experiments in B. terrestris.

Gamma irradiation of pollen and eradication of Israeli acute paralysis virus.

Honeybees and bumblebees are the most important pollinators of agricultural crops. For this purpose honeybees and bumblebees are reared and transported. A pathogen-free status of bees in general, is crucial. Indeed anthropogenic transports of hosts carrying parasites could alter the natural host/pathogen association, inducing an extra pathogenic stress. Therefore the creation of a pathogen-free rearing environment is needed. For bumblebees this is possible, as these species are reared in a closed environment. Although, a link remains between reared bumblebees and the outside bee community, as honeybee-collected pollen is essential food for bumblebee mass rearing. Here we evaluated if gamma irradiation can minimize the risk of this potential route of exposure and can inactivate viral particles present in honeybee-collected pollen. We show that 16.9kGy gamma irradiation induced a 100-1000 fold reduction on the ability of IAPV to cause mortality after injections. This result opens avenues toward rearing pathogen-free bumblebees and towards eliminating the risks of pathogen spillover to native wild bee species.