Complete genome sequence and genetic characterization of a novel segmented RNA virus infecting Nilaparvata lugens.

The brown planthopper (BPH), Nilaparvata lugens, is a significant agricultural pest capable of long-distance migration and transmission of viruses that cause severe disease in rice. In this study, we identified a novel segmented RNA virus in a BPH, and this virus exhibited a close relationship to members of a recently discovered virus lineage known as "quenyaviruses" within the viral kingdom Orthornavirae. This newly identified virus was named "Nilaparvata lugens quenyavirus 1" (NLQV1). NLQV1 consists of five positive-sense, single-stranded RNAs, with each segment containing a single open reading frame (ORF). The genomic characteristics and phylogenetic analysis support the classification of NLQV1 as a novel quenyavirus. Notably, all of the genome segments of NLRV contained the 5'-terminal sequence AUCUG. The characteristic virus-derived small interfering RNA (vsiRNA) profile of NLQV1 suggests that the antiviral RNAi pathway of the host BPH was activated in response to virus infection. These findings represent the first documented report of quenyaviruses in planthoppers, contributing to our understanding of quenyaviruses and expanding our knowledge of insect-specific viruses in planthoppers.

Molecular and biological characterization of a bunyavirus infecting the brown planthopper (Nilaparvata lugens).

A negative-strand symbiotic RNA virus, tentatively named Nilaparvata lugens Bunyavirus (NLBV), was identified in the brown planthopper (BPH, Nilaparvata lugens). Phylogenetic analysis indicated that NLBV is a member of the genus Mobuvirus (family Phenuiviridae, order Bunyavirales). Analysis of virus-derived small interfering RNA suggested that antiviral immunity of BPH was successfully activated by NLBV infection. Tissue-specific investigation showed that NLBV was mainly accumulated in the fat-body of BPH adults. Moreover, NLBV was detected in eggs of viruliferous female BPHs, suggesting the possibility of vertical transmission of NLBV in BPH. Additionally, no significant differences were observed for the biological properties between NLBV-infected and NLBV-free BPHs. Finally, analysis of geographic distribution indicated that NLBV may be prevalent in Southeast Asia. This study provided a comprehensive characterization on the molecular and biological properties of a symbiotic virus in BPH, which will contribute to our understanding of the increasingly discovered RNA viruses in insects.

Cross-kingdom RNA interference mediated by insect salivary microRNAs may suppress plant immunity.

Communication between insects and plants relies on the exchange of bioactive molecules that traverse the species interface. Although proteinic effectors have been extensively studied, our knowledge of other molecules involved in this process remains limited. In this study, we investigate the role of salivary microRNAs (miRNAs) from the rice planthopper Nilaparvata lugens in suppressing plant immunity. A total of three miRNAs were confirmed to be secreted into host plants during insect feeding. Notably, the sequence-conserved miR-7-5P is specifically expressed in the salivary glands of N. lugens and is secreted into saliva, distinguishing it significantly from homologues found in other insects. Silencing miR-7-5P negatively affects N. lugens feeding on rice plants, but not on artificial diets. The impaired feeding performance of miR-7-5P-silenced insects can be rescued by transgenic plants overexpressing miR-7-5P. Through target prediction and experimental testing, we demonstrate that miR-7-5P targets multiple plant genes, including the immune-associated bZIP transcription factor 43 (OsbZIP43). Infestation of rice plants by miR-7-5P-silenced insects leads to the increased expression of OsbZIP43, while the presence of miR-7-5P counteracts this upregulation effect. Furthermore, overexpressing OsbZIP43 confers plant resistance against insects which can be subverted by miR-7-5P. Our findings suggest a mechanism by which herbivorous insects have evolved salivary miRNAs to suppress plant immunity, expanding our understanding of cross-kingdom RNA interference between interacting organisms.

A pipeline contributes to efficient identification of salivary proteins in short-headed planthopper, Epeurysa nawaii.

Saliva, an oral secretion primarily originating from salivary glands (SGs), exert critical roles in the ongoing evolutionary interaction between insects and plants. However, identifying insect salivary components poses challenges due to the tiny size of insects, low secretion amounts, and the propensity for degradation after secretion. In this study, we developed a transcriptome-based approach to comprehensively analyze the salivary proteins of the short-headed planthopper, Epeurysa nawaii, a species with unique feeding habits on bamboo. A total of 165 salivary proteins were identified, with 114 secretory genes highly and specifically expressed in SGs. Consistent with most phloem-feeding insects, digestive enzymes, calcium-binding proteins, oxidoreductases, and a few previously reported salivary effectors were ubiquitously distributed in E. nawaii saliva. However, we also identified a substantial portion of salivary proteins exhibiting taxonomy specificity, including 60 E. nawaii-specific and 62 Delphacidae-specific proteins. These taxonomy-restricted proteins potentially play a role in insect adaptation to specific host plants. Our study provides an efficient pipeline for salivary protein identification and serves as a valuable resource for the functional characterization of effectors.

Salivary proteins potentially derived from horizontal gene transfer are critical for salivary sheath formation and other feeding processes.

Herbivorous insects employ an array of salivary proteins to aid feeding. However, the mechanisms behind the recruitment and evolution of these genes to mediate plant-insect interactions remain poorly understood. Here, we report a potential horizontal gene transfer (HGT) event from bacteria to an ancestral bug of Eutrichophora. The acquired genes subsequently underwent duplications and evolved through co-option. We annotated them as horizontal-transferred, Eutrichophora-specific salivary protein (HESPs) according to their origin and function. In Riptortus pedestris (Coreoidea), all nine HESPs are secreted into plants during feeding. The RpHESP4 to RpHESP8 are recently duplicated and found to be indispensable for salivary sheath formation. Silencing of RpHESP4-8 increases the difficulty of R. pedestris in probing the soybean, and the treated insects display a decreased survivability. Although silencing the other RpHESPs does not affect the salivary sheath formation, negative effects are also observed. In Pyrrhocoris apterus (Pyrrhocoroidea), five out of six PaHESPs are secretory salivary proteins, with PaHESP3 being critical for insect survival. The PaHESP5, while important for insects, no longer functions as a salivary protein. Our results provide insight into the potential origin of insect saliva and shed light on the evolution of salivary proteins.

Comparative transcriptomic analysis of salivary glands between the zoophytophagous Cyrtorhinus lividipennis and the phytozoophagous Apolygus lucorum.

BACKGROUND: Saliva plays a crucial role in shaping the feeding behavior of insects, involving processes such as food digestion and the regulation of interactions between insects and their hosts. Cyrtorhinus lividipennis serves as a predominant natural enemy of rice pests, while Apolygus lucorum, exhibiting phytozoophagous feeding behavior, is a destructive agricultural pest. In this study, a comparative transcriptome analysis, incorporating the published genomes of C.lividipennis and A.lucorum, was conducted to reveal the role of salivary secretion in host adaptation. RESULTS: In contrast to A.lucorum, C.lividipennis is a zoophytophagous insect. A de novo genome analysis of C.lividipennis yielded 19,706 unigenes, including 16,217 annotated ones. On the other hand, A.lucorum had altogether 20,111 annotated genes, as obtained from the published official gene set (20,353 unigenes). Functional analysis of the top 1,000 salivary gland (SG)-abundant genes in both insects revealed that the SG was a dynamically active tissue engaged in protein synthesis and secretion. Predictions of other tissues and signal peptides were compared. As a result, 94 and 157 salivary proteins were identified in C.lividipennis and A.lucorum, respectively, and were categorized into 68 and 81 orthogroups. Among them, 26 orthogroups were shared, potentially playing common roles in digestion and detoxification, including several venom serine proteases. Furthermore, 42 and 55 orthogroups were exclusive in C.lividipennis and A.lucorum, respectively, which were exemplified by a hyaluronidase in C.lividipennis that was associated with predation, while polygalacturonases in A.lucorum were involved in mesophyll-feeding patterns. CONCLUSIONS: Findings in this study provide a comprehensive insight into saliva secretions in C.lividipennis and A.lucorum via a transcriptome approach, reflecting the intricate connections between saliva secretions and feeding behaviors. It is found that conserved salivary secretions are involved in shaping the overlapping feeding patterns, while a plethora of unique salivary secretions may drive the evolution of specific feeding behaviors crucial for their survival. These results enhance our understanding of the feeding mechanisms in different insects from the perspective of saliva and contribute to future environmentally friendly pest control by utilizing predatory insects.

Co-option of a non-retroviral endogenous viral element in planthoppers.

Non-retroviral endogenous viral elements (nrEVEs) are widely dispersed throughout the genomes of eukaryotes. Although nrEVEs are known to be involved in host antiviral immunity, it remains an open question whether they can be domesticated as functional proteins to serve cellular innovations in arthropods. In this study, we found that endogenous toti-like viral elements (ToEVEs) are ubiquitously integrated into the genomes of three planthopper species, with highly variable distributions and polymorphism levels in planthopper populations. Three ToEVEs display exon‒intron structures and active transcription, suggesting that they might have been domesticated by planthoppers. CRISPR/Cas9 experiments revealed that one ToEVE in Nilaparvata lugens, NlToEVE14, has been co-opted by its host and plays essential roles in planthopper development and fecundity. Large-scale analysis of ToEVEs in arthropod genomes indicated that the number of arthropod nrEVEs is currently underestimated and that they may contribute to the functional diversity of arthropod genes.

Horizontally Transferred Salivary Protein Promotes Insect Feeding by Suppressing Ferredoxin-Mediated Plant Defenses.

Herbivorous insects such as whiteflies, planthoppers, and aphids secrete abundant orphan proteins to facilitate feeding. Yet, how these genes are recruited and evolve to mediate plant-insect interaction remains unknown. In this study, we report a horizontal gene transfer (HGT) event from fungi to an ancestor of Aleyrodidae insects approximately 42 to 190 million years ago. BtFTSP1 is a salivary protein that is secreted into host plants during Bemisia tabaci feeding. It targets a defensive ferredoxin 1 in Nicotiana tabacum (NtFD1) and disrupts the NtFD1-NtFD1 interaction in plant cytosol, leading to the degradation of NtFD1 in a ubiquitin-dependent manner. Silencing BtFTSP1 has negative effects on B. tabaci feeding while overexpressing BtFTSP1 in N. tabacum benefits insects and rescues the adverse effect caused by NtFD1 overexpression. The association between BtFTSP1 and NtFD1 is newly evolved after HGT, with the homologous FTSP in its fungal donor failing to interact and destabilize NtFD1. Our study illustrates the important roles of horizontally transferred genes in plant-insect interactions and suggests the potential origin of orphan salivary genes.

Diversity of RNA viruses in agricultural insects.

Recent advancements in next-generation sequencing (NGS) technology and bioinformatics tools have revealed a vast array of viral diversity in insects, particularly RNA viruses. However, our current understanding of insect RNA viruses has primarily focused on hematophagous insects due to their medical importance, while research on the viromes of agriculturally relevant insects remains limited. This comprehensive review aims to address the gap by providing an overview of the diversity of RNA viruses in agricultural pests and beneficial insects within the agricultural ecosystem. Based on the NCBI Virus Database, over eight hundred RNA viruses belonging to 39 viral families have been reported in more than three hundred agricultural insect species. These viruses are predominantly found in the insect orders of Hymenoptera, Hemiptera, Thysanoptera, Lepidoptera, Diptera, Coleoptera, and Orthoptera. These findings have significantly enriched our understanding of RNA viral diversity in agricultural insects. While further virome investigations are necessary to expand our knowledge to more insect species, it is crucial to explore the biological roles of these identified RNA viruses within insects in future studies. This review also highlights the limitations and challenges for the effective virus discovery through NGS and their potential solutions, which might facilitate for the development of innovative bioinformatic tools in the future.

Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection.

The whitefly Bemisia tabaci is one of the most destructive pests worldwide, and causes tremendous economic losses. Tobacco Nicotiana tabacum serves as a model organism for studying fundamental biological processes and is severely damaged by whiteflies. Hitherto, our knowledge of how tobacco perceives and defends itself against whiteflies has been scare. In this study, we analyze the gene expression patterns of tobacco in response to whitefly infestation. A total of 244 and 2417 differentially expressed genes (DEGs) were identified at 12 h and 24 h post whitefly infestation, respectively. Enrichment analysis demonstrates that whitefly infestation activates plant defense at both time points, with genes involved in plant pattern recognition, transcription factors, and hormonal regulation significantly upregulated. Notably, defense genes are more intensely upregulated at 24 h post infestation than at 12 h, indicating an increased immunity induced by whitefly infestation. In contrast, genes associated with energy metabolism, carbohydrate metabolism, ribosomes, and photosynthesis are suppressed, suggesting impaired plant development. Taken together, our study provides comprehensive insights into how plants respond to phloem-feeding insects, and offers a theoretical basis for better research on plant-insect interactions.

A chromosome-level genome assembly of Stenchaetothrips biformis and comparative genomic analysis highlights distinct host adaptations among thrips.

Insects have a limited host range due to genomic adaptation. Thysanoptera, commonly known as thrips, occupies distinct feeding habitats, but there is a lack of comparative genomic analyses and limited genomic resources available. In this study, the chromosome-level genome of Stenchaetothrips biformis, an oligophagous pest of rice, is assembled using multiple sequencing technologies, including PacBio, Illumina short-reads, and Hi-C technology. A 338.86 Mb genome is obtained, consisting of 1269 contigs with a contig N50 size of 381 kb and a scaffold N50 size of 18.21 Mb. Thereafter, 17,167 protein-coding genes and 36.25% repetitive elements are annotated. Comparative genomic analyses with two other polyphagous thrips, revealing contracted chemosensory-related and expanded stress response and detoxification gene families in S. biformis, potentially facilitating rice adaptation. In the polyphagous thrips species Frankliniella occidentalis and Thrips palmi, expanded gene families are enriched in metabolism of aromatic and anthocyanin-containing compounds, immunity against viruses, and detoxification enzymes. These expansion gene families play crucial roles not only in adapting to hosts but also in development of pesticide resistance, as evidenced by transcriptome results after insecticides treatment. This study provides a chromosome-level genome assembly and lays the foundation for further studies on thrips evolution and pest management.

Long-wave opsin involved in body color plastic development in Nilaparvata lugens.

BACKGROUND: As one of the components of visual photopigments in photoreceptor cells, opsin exhibits different spectral peaks and plays crucial roles in visual function. Besides, it is discovered to evolve other functions despite color vision. However, research on its unconventional function is limited nowadays. With the increase in genome database numbers, various numbers and types of opsins have been identified in insects due to gene duplications or losses. The Nilaparvata lugens (Hemiptera) is a rice pest known for its long-distance migration capability. In this study, opsins were identified in N. lugens and characterized by genome and transcriptome analyses. Meanwhile, RNA interference (RNAi) was carried out to investigate the functions of opsins, and then the Illumina Novaseq 6000 platform-based transcriptome sequencing was performed to reveal gene expression patterns. RESULTS: Four opsins belonging to G protein-coupled receptors were identified in the N. lugens genome, including one long-sensitive opsin (Nllw) together with two ultraviolet-sensitive opsins (NlUV1/2) and an additional new opsin with hypothesized UV peak sensitivity (NlUV3-like). A tandem array of NlUV1/2 on the chromosome suggested the presence of a gene duplication event, with similar exons distribution. Moreover, as revealed by spatiotemporal expression, the four opsins were highly expressed in eyes with age-different expression levels. Besides, RNAi targeting each of the four opsins did not significantly affect the survival of N. lugens in phytotron, but the silencing of Nllw resulted in the melanization of body color. Further transcriptome analysis revealed that silencing of Nllw resulted in up-regulation of a tyrosine hydroxylase gene (NlTH) and down-regulation of an arylalkylamine-N-acetyltransferases gene (NlaaNAT) in N. lugens, demonstrating that Nllw is involved in body color plastic development via the tyrosine-mediated melanism pathway. CONCLUSIONS: This study provides the first evidence in a Hemipteran insect that an opsin (Nllw) takes part in the regulation of cuticle melanization, confirming a cross-talk between the gene pathways underlying the visual system and the morphological differentiation in insects.

Prevalence of Reversed Genome Organizations for Viruses in the Family Iflaviridae, Order Picornavirales.

Viruses in the order Picornavirales possess a positive-strand RNA genome that encodes structural proteins (SPs) and nonstructural proteins (NSPs). According to the recent report of the International Committee on Taxonomy of Viruses (ICTV), there are 8 families in Picornavirales, and monopartite picornaviruses in each family exhibit distinct types of genome organizations with rearranged genes coding for SPs and NSPs, namely, TypeI (5'-SPs-NSPs-3') and TypeII (5'-NSPs-SPs-3'). In the present study, 2 iflaviruses with the 2 genome types were unexpectedly identified in a damselfly host species, suggesting that these 2 genome types coexisted in the same host species, and the families of order Picornavirales might be more complex than previously thought. The consequent systematic homologous screening with all the publicly available picornaviruses successfully revealed a considerable number of candidates rearranged genome types of picornaviruses in various families of Picornavirales. Subsequently, phylogenetic trees were reconstructed based on RNA dependent RNA polymerase and coat protein, which evidently confirmed the prevalence of the 10 typeII iflaviruses in the Iflaviridae family. This suggests that genome types may not be relevant to viral taxonomy in this family. However, candidate picornaviruses with reversed genome types in the Secoviridae and Dicistroviridae families require further investigation. All in all, as the number of newly discovered viruses increases, more viruses with non-canonical genome arrangements will be uncovered, which can expand our current knowledge on the genome complexity and evolution of picornaviruses. IMPORTANCE Monopartite viruses in the order Picornavirales exhibit distinct genome arrangement of nonstructural proteins and structural proteins for each of the 8 families. Recent studies indicated that at least 4 ifla-like viruses possessed reversed genome organization in the family Iflaviridae, raising the possibility that this phenomenon may commonly present in different families of picornaviruses. Since we discovered 2 iflaviruses with exchanged structural and nonstructural proteins simultaneously in the damselfly, a systematic screening was subsequently performed for all of the current available picornaviruses (1,543 candidates). The results revealed 10 picornaviruses with reversed genome organization in the family Iflaviridae, implying that this phenomenon might prevalence in the order Picornavirales. These results will contribute to a better understanding for the future study on the genome complexity and taxonomy of picornaviruses.

Planthopper salivary sheath protein LsSP1 contributes to manipulation of rice plant defenses.

Salivary elicitors secreted by herbivorous insects can be perceived by host plants to trigger plant immunity. However, how insects secrete other salivary components to subsequently attenuate the elicitor-induced plant immunity remains poorly understood. Here, we study the small brown planthopper, Laodelphax striatellus salivary sheath protein LsSP1. Using Y2H, BiFC and LUC assays, we show that LsSP1 is secreted into host plants and binds to salivary sheath via mucin-like protein (LsMLP). Rice plants pre-infested with dsLsSP1-treated L. striatellus are less attractive to L. striatellus nymphs than those pre-infected with dsGFP-treated controls. Transgenic rice plants with LsSP1 overexpression rescue the insect feeding defects caused by a deficiency of LsSP1 secretion, consistent with the potential role of LsSP1 in manipulating plant defenses. Our results illustrate the importance of salivary sheath proteins in mediating the interactions between plants and herbivorous insects.

High-throughput sequencing yields a complete mitochondrial genome of the rice thrips, Stenchaetothrips biformis (Thysanoptera: Thripidae).

Rice thrips, Stenchaetothrips biformis (Bagnall, 1913), are one of the destructive pests of rice. Here, the complete mitochondrial genome of S. biformis was sequenced using high-throughput sequencing. The mitogenome is 15,359 bp long with an A + T content of 76.94%, which contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs), 2 ribosomal RNA genes (rRNAs) and 2 putative control regions (CRs). The phylogenetic analysis showed that S. biformis is closely related to Thrips imaginis and Thrips palmi. This new mitochondrial genome data can be better used to provide a basis for studies of the mitochondrial evolution of Thysanoptera.

Complete genome analysis of a novel chuvirus from a southern green stink bug (Nezara viridula).

A novel chuvirus from a southern green stink bug (Nezara viridula) was identified by RNA sequencing in this study and was tentatively named "Ningbo southern green stink bug chuvirus 1" (NBSGSBV-1). The complete genome sequence of NBSGSBV-1 consists of 11,375 nucleotides, and the genome was found to be circular by 'around-the-genome' reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing. Three open reading frames (ORFs) were predicted in the NBSGSBV-1 genome, encoding a large polymerase protein (L protein), a glycoprotein (G protein), and a nucleocapsid protein (N protein). A phylogenetic tree was constructed based on all of the currently available RNA-dependent RNA polymerase amino acid sequences of viruses of the family Chuviridae, and NBSGSBV-1 was found to cluster together with Sanya chuvirus 2 and Hubei odonate virus 11, indicating that NBSGSBV-1 might belong to the genus Odonatavirus. Five conserved sites were identified in the L proteins of NBSGSBV-1 and other chuviruses. The abundance and characteristics of the NBSGSBV-1-derived small interfering RNAs suggested that NBSGSBV-1 actively replicates in the host insect. To the best of our knowledge, this is the first report of a chuvirus identified in a member of the insect family Pentatomidae. The discovery and characterization of NBSGSBV-1 will help us to understand the diversity of chuviruses in insects.

Complete genome analysis of a novel iflavirus from the spotted lanternfly Lycorma delicatula.

The spotted lanternfly (Lycorma delicatula) is an invasive pest that causes serious economic losses in fruit and wood production. Here, we identified a novel iflavirus named "Lycorma delicatula iflavirus 1" (LDIV1), in a spotted lanternfly. The full genome sequence of LDIV1 is 10,222 nt in length and encodes a polyprotein containing a picornavirus capsid-protein-domain-like domain, a cricket paralysis virus capsid superfamily domain, an RNA helicase domain, a peptidase C3 superfamily domain, and an RNA-dependent RNA polymerase (RdRp) domain. LDIV1 replicates in the host insect and activates small interfering RNA (siRNA)-based host antiviral immunity. Phylogenetic analysis demonstrated that LDIV1 is most closely related to an unspecified member of the order Picornavirales, with 61.7% sequence identity in the RdRp region and 57.6% sequence identity in the coat protein region, and thus meets the demarcation criteria for new species in the genus Iflavirus. To the best of our knowledge, LDIV1 is the first virus discovered in L. delicatula.

Characterization of Two Novel Insect-Specific Viruses Discovered in the Green Leafhopper, Cicadella viridis.

Insect-specific viruses (ISV) are one of the most promising agents for the biological control of insects. The green leafhopper, Cicadella viridis (Linnaeus), is an significant pest in agriculture, and causes economic losses to fruit trees, willows, and field crops. As a representative species of the large family Cicadellidae, ISVs in C. viridis have never been studied, to date. In this study, the full genome sequences of two ISVs, named Cicadella viridis iflavirus1 (CvIfV1), and Cicadella viridis nido-like virus 1 (CvNiLV1), were revealed using a metatranscriptomic approach. A homology search and phylogenetic analysis indicated that CvIfV1 is a new member in the family Iflaviridae (genus Iflavirus) with a typical iflavirus genome organization, whereas CvNiLV1 belongs to the unclassified clade/family of the order Nidovirales. In addition, analysis of virus-derived small interfering RNAs (vsiRNAs) was performed to investigate the antiviral RNA interference (RNAi) response of C. viridis. The vsiRNAs exhibit typical patterns produced by host siRNA-mediated antiviral immunity, including a preference of 21-nt vsiRNAs derived equally from the sense and antisense genomic strands, and a strong A/U bias in the 5'-terminus of the viral genomes. Our study provides valuable information for ISVs in leafhoppers for the first time, which might prove useful in the control of C. viridis in future.

Draft Genome Sequence of a Novel Iflavirus from Leafhoppers (Exitianus capicola) in Iraq.

The draft genome sequence of a novel virus was determined from a leafhopper insect (Exitianus capicola) via transcriptomic sequencing and rapid amplification of cDNA ends. Results of the sequence alignment suggested that the new virus is a member of the genus Iflavirus, family Iflaviridae. Thus, it was proposed to be named Exitianus capicola iflavirus 2 (ECIV2). ECIV2 has a single positive-strand RNA genome of 7,821 nucleotides [excluding the poly(A) tail] containing a single open reading frame (ORF) with typical iflavirus conserved domains, including a picornavirus capsid protein-like domain, a cricket paralysis virus capsid-like domain, a helicase, a peptidase, and an RNA-dependent RNA polymerase domain.