The Novel Agrotis ipsilon Nora Virus Confers Deleterious Effects to the Fitness of Spodoptera frugiperda (Lepidoptera: Noctuidae)

Abstract

In the present study, we identified a novel, positive-sense single-stranded RNA virus in the Chinese black cutworm, Agrotis ipsilon. It has a genome length of 11,312 nucleotides, excluding the poly(A) tails, and contains five open reading frames. The ORF2 encodes the conserved domains of RNA helicase and RNA-dependent RNA polymerase, while ORF4 and 5 encode three viral proteins. Herein, the A. ipsilon virus was clustered with a Helicoverpa armigera Nora virus and was thus provisionally named “Agrotis ipsilon Nora virus” (AINV). AINV was successfully transmitted into a novel host, Spodoptera frugiperda, through injection, causing a stable infection. This found the possibility of horizontal AINV transmission among moths belonging to the same taxonomic family. Nonetheless, AINV infection was deleterious to S. frugiperda and mainly mediated by antiviral and amino acid metabolism-related pathways. Furthermore, the infection significantly increased the S. frugiperda larval period but significantly reduced its moth eclosion rate. It suggests that AINV is probably to be a parasitic virus of S. frugiperda.

Introduction

The second-generation sequencing technology has expanded our understanding of the diversity of insect viruses (Xu et al., 2014; Liu et al., 2017; Shi et al., 2018; Yang et al., 2019). Recently, picorna-like symbiotic viruses are commonly identified in insects (Shi et al., 2016; Yang et al., 2016; Cholleti et al., 2018). Taxonomically, picorna-like viruses belong to the order Picornavirales, which share some properties with members of the family Picornaviridae (Le Gall et al., 2008). Generally, picorna-like viruses possess a (+) ssRNA genome (positive-sense, single-stranded RNA), and their translated polyproteins are cleaved into structural and non-structural proteins. Their pathogenicity is diverse, ranging from lethal to symbiotic (de Miranda et al., 2010; Jakubowska et al., 2014).

The Nora virus was first reported in Drosophila melanogaster and then proposed to represent a new picorna-like virus family (Habayeb et al., 2006). Subsequent studies have revealed that Nora viruses interfere with the D. melanogaster RNAi system (van Mierlo et al., 2012, 2014). Recently, additional members of the Nora virus clade have been identified in arthropods through high-throughput sequencing technologies (Shi et al., 2016). Agrotis ipsilon (Hufnagel) is a worldwide pest that causes significant damage to vegetables and grains. Despite its interactions with its endosymbionts being thoroughly studied (Shi et al., 2013), its viral landscape remains ambiguous. Herein, we identified a novel Nora virus in Chinese A. ipsilon through RNA sequencing. The virus was provisionally named as “Agrotis ipsilon Nora virus” (AINV). The complete AINV genome was subsequently determined, and its phylogenetic position and gene expression patterns were uncovered among the A. ipsilon developmental stages. AINV was also introduced into wild Spodoptera frugiperda, an invasive pest in China (Sun et al., 2021), to evaluate its effects on fitness and gene expressions.

Materials and Methods

Insect Rearing

Colonies of A. ipsilon and S. frugiperda were established with caterpillars collected in Shanghai, China, in 2018 and Zhengzhou, China, in 2020, respectively. The caterpillars were fed with artificial diets at 25°C (±1°C) under a 14/10h light/dark cycle. Adult moths were fed with 10% honey water.

Virus Detection in A. ipsilon Using RNA-Seq

Total RNA was extracted from A. ipsilon samples using the TRIzol reagent (Invitrogen, Carlsbad, CA, United States) following the manufacturer’s protocol. The samples were collected at four developmental stages, i.e., egg and larvae (n=30), pupae (n=30), and adults (n=30) in both male and female insects. RNA quality and quantity were assessed on an Agilent 2,100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States) and RNase-free agarose gel electrophoresis. The mRNA in the RNA samples was enriched by removing rRNA using the Ribo-ZeroTM Magnetic Kit (Epicentre, Madison, WI, United States). The mRNA libraries were subsequently prepared by BGI (Shenzhen, China) and sequenced on an Illumina Hiseq 2000 platform. De novo transcriptome assemblies were performed using the SOAPdenovo software with default parameters (Xie et al., 2014). Protein sequences of sRNA viruses (txid 439,488), dsRNA viruses (txid 35,325), and Delta viruses (txid 39,759) were retrieved from the NCBI RefSeq database (Reference Sequences) and employed as subjects in local blastx searches to uncover possible viral fragments in assembled unigenes. The non-virus originating hits were filtered by performing online blastx searches against the NCBI nr database (non-redundant protein sequences). An e-value threshold of 1×10⁻⁵ was used in these searches.

AINV Genome Amplification and Phylogenetic Analysis

Full-length cDNA sequences of the AINV were obtained using the 3′ and 5′ rapid amplification of cDNA ends (RACE) system (Life Technologies, Carlsbad, CA, United States), following the manufacturer’s instructions. The sequences were subsequently verified through amplification and sequencing using specific primers. Total RNA for specific RT-PCR was extracted using the TaKaRa MiniBEST Universal RNA Extraction Kit (Takara, Dalian, China), followed by RT-PCR amplification of the AINV fragments using PrimeScript^™ One-Step RT-PCR Kit Ver.2 (Dye Plus; Takara). The gene-specific primers used herein were listed in Table 1. The standard genetic code was employed to predict the open reading frames of the AINV using the NCBI online ORF finder program.¹ Conserved domains within the ORFs were subsequently predicted using the NCBI conserved domain database v3.17,² with threshold of 1×10⁻².

The RNA-dependent RNA polymerase (RdRp) protein sequences from taxa within Picornavirales and Caliciviridae were retrieved to uncover the phylogenetic position of AINV. The sequences were aligned using the MUSCLE program in MEGA 7.0 (Kumar et al., 2016) and then trimmed using trimAl to remove the poorly aligned regions (Capella-Gutierrez et al., 2009). Phylogenetic analysis was performed in IQ-TREE 1.6.6 (Nguyen et al., 2015), followed by resampling 1,000 ultrafast bootstraps to assess the support for each node. The substitution model based on the Bayesian information criterion in ModelFinder (Kalyaanamoorthy et al., 2017) was selected to improve the accuracy of phylogenetic estimates.

AINV Transmission and Quantification in a Novel Host

An AINV infected liquid was formulated following the methods described by Xu et al. (2020). The AINV infected fourth A. ipsilon larvae were ground in liquid nitrogen and then homogenized with 1ml PBS buffer (0.01M, pH 7.4). The homogenate was then centrifuged at 6500×g for 15min at 4°C, followed by filtration of the liquid supernatant through 0.2μm syringe filters (Pall Corporation, New York, United States). The filtrate (10μl) was subsequently injected into the third S. frugiperda larvae using a Hamilton Microliter syringe (705N) to introduce the AINV into the host. Control hosts were injected with an equal volume of sterile water. The S. frugiperda larvae were reared in room, followed by detection of AINV infections in the moths using AINV specific primers, AINV_72F and AINV_44R (Table 1). The S. frugiperda larvae were reared in room, followed by detection of AINV infections in the moths using AINV specific primers, AINV_72F and AINV_44R (Table 1). AINV titers in S. frugiperda were quantified based on the ORF2 (polyprotein) copies and then normalized with S. frugiperda β-actin gene copies. The qPCR reactions were performed in a CFX384 Touch^™ Real-Time PCR Detection System (Bio-Rad, United States) using the SYBR Green method. The cDNA used in qPCR reactions was synthesized using the PrimeScript^™ RT reagent Kit with gDNA Eraser (Takara). The reactions were carried out in a 10μl volume containing 5 ul of 2×TB Green^® Premix Ex Taq^™ (Takara), 0.4μl of each primer (10μm), 1μl of cDNA, and 3.2μl sterile water. The qPCR primers used are listed in Table 1. AINVQF and AINVQR2 amplified a 102bp fragment of the AINV, while SFACTINF and SFACTINR amplified a 156bp fragment of the S. frugiperda β-actin gene. The qPCR reactions for quantifying AINV titers in infected male and female S. frugiperda larvae were done in triplicate. Four replicates of AINV infected larvae were used in these qPCRs.

Effect of AINV Infection on S. frugiperda

Six fitness indices, including larval period, pupal period, pupal weight, pupation rate, eclosion rate, and sex ratio, were measured in AINV positive and negative S. frugiperda populations to uncover the effects of AINV infection on the novel host. The GraphPad Prism 8.3.0 software was then employed to analyze and visualize the data for the six indices.

RNA-seq was also employed to reveal the effects of AINV infection on the gene expression of female S. frugiperda adults. The first day hatched moths were subjected to RNA-seq, and three replicates of AINV positive and negative females were used following the aforementioned sequencing procedure. The HISAT2 (2.0.4; Kim et al., 2015) alignment program was used to map the obtained reads to the S. frugiperda published genome (ZJU_Sfru_1.0, GCF_011064685.1). Then, the mapped reads were assembled using the StringTie (Version: 1.3.4) sequence assembler (Pertea et al., 2015), followed by merging of the transcripts using gffcompare (Version: 0.9.8).³ The ballgown R package⁴ was subsequently used to estimate the expression levels of all transcripts and to detect the differentially expressed mRNAs by calculating FPKM (Fragments Per Kilobase Million; FPKM = [total exon fragments/mapped reads(millions)×exon length(kB)]). The edgeR R package (Robinson et al., 2010) was employed to filter the differentially expressed genes (DEGs) at a threshold absolute log₂ (FC, fold change)≥1 and false discovery rate (FDR)<0.05. The DEGs were shown with volcano plot, which was visualized using the ggplot2, ggrepel, gridExtra, and ggthemes R packages. Principal component analysis (PCA) of the samples was further determined with the gene expression matrix. The analysis was then visualized using the stats R package. GO and KEGG enrichment results of the DEGs were visualized using the clusterProfiler R package (Yu et al., 2012), while the heat maps were visualized using TBtools (Chen et al., 2020).

Results

RNA Virus Detection in A. ipsilon

There were 62,933 unigenes within the A. ipsilon transcriptome. Two virus-originated fragments were detected as: unigene-23,909 (5,239 nucleotides) and unigene-23,910 (6,052 nucleotides). Blastp searches revealed that both unigenes were similar to the protein sequences of a Nora virus isolated in Helicoverpa armigera (HaNV, accession no: MK033133). These findings suggested that the unigenes probably originated from a Nora virus member and was tentatively named AINV. Supplementary Figure S1 showed the alignment of the two unigenes. The overlapping region between unigene-23,909 and unigene-23,910 was subsequently amplified, and the sequences were manually verified using specific primers (Table 1). Furthermore, following the methods described by Yang et al., (2019), we enriched the virus particles in AINV positive population and then observed them with a transmission electron microscope. The result was provided in Supplementary Figure S3.

AINV Genome Characterization, Phylogenetic Analysis, and Quantification in A. ipsilon

The RACE alignments (data not shown) indicated that the complete genome of AINV was 11,312 nucleotides in length, excluding the poly(A) tails (Figure 1A). The locations of unigene-23,909 and unigene-23,910 on AINV genome were position 25 to position 5,263, and position 5,261 to position 11,312, respectively. The genome shared 85.69% nucleotide identity with HaNV with a 96% coverage. The ORF finder results suggested that AINV had a (+) ssRNA genome containing five ORFs. The conserved replication polyprotein domains were in ORF2. Comparison of the replication polyproteins of AINV and HaNV revealed a 14.60% genetic divergence between them. In the same line, sequences homologous to the three structural proteins of Drosophila Nora viruses, named viral protein (VP) 4A, VP4B, and VP4C (Ekstrom et al., 2011), were also detected in ORF4 and ORF5 of the AINV.

Figure 1

The LG substitution model with unequal amino acid frequencies (+F) and rate variation among sites (+R10) was selected in the phylogenetic analysis. The results strongly supported the monophyly of the Nora virus clade (UFBoot=99; Supplementary Figure S2). Furthermore, AINV shared a sister relationship with HaNV (UFBoot=100) within the Nora virus clade (Figure 1B) and further clustered with Nora viruses detected in Drosophila and Odonata (UFBoot=78). These findings confirmed that AINV was a new Nora virus member.

Both unigenes exhibited the same expression pattern in A. ipsilon, with higher expression levels in larvae and pupae than in eggs and adults (Figure 1C). These findings suggested that AINV titers varied among the developmental stages of A. ipsilon, with the larvae and adults having the highest and lowest AINV titers, respectively.

AINV Novel Host Transmission

The AINV liquid was injected into 45 third S. frugiperda larvae, while another 45 larvae were injected with sterile water to act as the control. The adults of these larvae were then harvested and subjected to AINV detection using specific primers. Notably, all the adults in the AINV liquid injected group were positive for AINV, but none was positive in the controls. These findings suggested that AINV had a highly accurate horizontal transmission rate even between different hosts with the liquid injection method. Further quantification of the AINV based on the ORF2 polyprotein copies and subsequent normalization with S. frugiperda β-actin gene copies revealed the highest AINV titers in the larvae (1.77±0.32, n=4). In addition, male (1.56±0.05) and female (1.57±0.07) S. frugiperda had almost similar AINV titers. Similarly, there was no significant difference in AINV titer between the larvae and adults (F=1.086, p=0.388; Figure 2).

Figure 2

Effects of AINV on Fitness and Gene Expression of S. frugiperda

AINV infection significantly increased the larval (p<0.05) and pupal growth periods and decreased the pupal weight (Figure 3A). There were no significant differences in the pupae rate and sex ratio between the AINV positive and negative populations (Figures 3B,C). However, the AINV negative population had a higher pupae rate than the positive population. Similarly, the moth eclosion rate was significantly decreased in the AINV positive population (p<0.05; Figure 3D). Moreover, we also observed that the AINV positive females laid about average 30 eggs that are much less than the AINV negative females, which can averagely lay more than 500 eggs under room conditions. It indicated that AINV infection would significantly affect the reproductive capacity of S. frugiperda females.

Figure 3

PCA results of AINV positive and negative adult subjected to RNA-seq revealed that the AINV positive and negative replicates were well distinguished in females (Figure 4A), but not in the males (Figure 5A). This finding suggested that AINV infection strongly changed the gene expression pattern in the S. frugiperda females, nor the males. Hence, in this study, we mainly uncovered the effects of AINV on the gene expression patterns in S. frugiperda females. However, the analyses of males were provided in Figure 5.

Figure 4

Figure 5

Based on the log₂FC>1 or log₂FC<−1 and FDR<0.05 threshold, 627 unigenes were filtered as DEGs. Among them, 26 were significantly downregulated while 601 were upregulated in AINV positive females compared to the AINV negative females (Figure 4B). The enriched GO terms of the DEGs included a response to virus infection, iron-binding, amino acid metabolism and transport, and insect metamorphosis and mating process (FDR<0.05; Figure 4C). Similarly, the enriched KEGG pathways of the DEGs included the lysosome, ECM-receptor interaction, amino acids biosynthesis and metabolism, and longevity regulating pathways (FDR<0.05; Figure 4D). Moreover, all the DEGs grouped into these KEGG pathways were upregulated in AINV positive females except DEG118270764 (Figure 6).

Figure 6

Discussion

The Nora virus is a single-stranded, positive-sense RNA virus that was first reported in two Drosophila species: D.melanogaster and Drosophila simulans (Habayeb et al., 2006). Thereafter, metagenomic analyses have revealed novel Nora-like viruses in Drosophila (van Mierlo et al., 2014; Medd et al., 2018), Nasonia (Oliveira et al., 2010), honey bee (Remnant et al., 2017), and diverse environmental samples (Shi et al., 2016). Recently, Nora-like virus infections were also detected in two lepidopteran pests: S. exigua and H. armigera (Yang et al., 2019). Herein, AINV clustered with the H. armigera Nora virus with robust supports and grouped into the Nora virus clade suggesting that A. ipsilon is the third lepidopteran natural host of Nora virus. Nora virus shares a conserved picornavirus-like helicase-protease-replicase (H-P-Rep) cassette and has a unique genome organization pattern (Habayeb et al., 2006; Ekstrom et al., 2011). The taxonomic position of the Nora virus is further confirmed by the recent structure analyses of its capsid proteins (Laurinmaki et al., 2020). The virus has a relatively large RNA genome, comprising more than 11,000 nucleotides in length (Habayeb et al., 2006; Laurinmaki et al., 2020). AINV has 11,312 nucleotides and which include five ORFs. The ORF2 of AINV encodes the picornavirus-like replicative cassette, while ORF4 and ORF5 encode the three structural proteins. The AINV ORFs have short overlaps, indicating that ribosomal frameshifting maybe involved in the expression of its proteins (Dreher and Miller, 2006).

Nora virus can cause persistent infections in Drosophila, and these infections are not affected by the Drosophila RNAi pathways (Habayeb et al., 2009b). The Drosophila Nora virus is horizontally transmitted through the fecal-oral route, causing mild effects on the fitness of infected flies (Habayeb et al., 2009a). A recent study postulated that geotaxis dysfunction is a phenotypic hallmark of Drosophila Nora virus infection, which causes the infected flies to significantly lose their climbing ability (Rogers et al., 2020). Notably, the Nora virus is horizontally and vertically transmitted within H. armigera populations, with unclear pathogenic effects on the infected moths (Yang et al., 2019). Similarly, AINV stably infects the A. ipsilon population, with no clear negative effects on the fitness of A. ipsilon (data not shown). As such, AINV is probably a mutualistic virus in A. ipsilon. Novel partiti-like viruses can horizontally transfer among congener lepidopteran hosts and subsequently cause deleterious effects on the new hosts (Xu et al., 2020). These findings strongly suggest that the novel viruses are potential biological resources for pest management. S. frugiperda is a destructive crop pest discovered in China in 2018 (Sun et al., 2021). Herein, AINV was injected into S. frugiperda and established stable infections, indicating that AINV can horizontally transfer among the Noctuidae hosts. Our previous study indicates that relative to the oral transmitted pathway, microinjection has its own advantage in the horizontal transmission of novel insect viruses (Xu et al., 2020). Hence, in this study, we establish the AINV infection S. frugiperda population by microinjection. However, we also notice that Nora virus can horizontally transmitted in its original host with perfect transmitted rate, by sharing the virus contamination food (Yang et al., 2019). It implies that AINV is probably can also be transmitted among A. ipsilon individuals through food-borne transmission route. AINV titers in A. ipsilon and S. frugiperda varied among its host developmental stages, with the larvae having the highest titer. AINV infection also had some side effects on the fitness of its new host, S. frugiperda. It significantly extended the S. frugiperda larval period but significantly reduced its moth eclosion rate. Moreover, after AINV infection, S. frugiperda pupal weight loss, it indicates that AINV and S. frugiperda lack a long co-evolutionary history which causes AINV to be a parasitic virus of S. frugiperda.

Previous studies postulate that the Drosophila Nora virus structural protein has RNAi suppressive activity in its natural host (van Mierlo et al., 2012, 2014). Herein, RNA-seq analysis revealed a higher expression of Argonaute-2 gene (AGO2) in AINV positive S. frugiperda females than the negative ones (Fold change=1.47, FDR>0.05), indicating that AINV stimulates the RNAi antiviral defense in the new host. Nora virus infection upregulates immune genes in Drosophila (Cordes et al., 2013; Lopez et al., 2018). GO and KEGG enrichments herein revealed no immune-related pathway in response to AINV infection. However, some potential antiviral pathways, such as ko04142 (lysosome), were uncovered (Du and Jin, 2017). Notably, ko04512 (glycine, serine, and threonine metabolism; Ding et al., 2019) and ko00260 (ECM-receptor interaction; Guo et al., 2015) have been reported to respond to animal virus infections. Herein, the DEGs are mainly enriched in the amino acid-related pathways, indicating that AINV infection affects the amino acid metabolism in S. frugiperda.

Conclusion

This study reports a new Nora virus infecting A. ipsilon and provisionally names it AINV. Its genome has 11,312 nucleotides, which include five ORFs. AINV was successfully transmitted into a novel host, S. frugiperda, through injection, causing stable infection, suggesting horizontal AINV transmission among moths of the same taxonomic family. Furthermore, the fitness measurements and RNA-seq analysis showed that the AINV infection was deleterious to S. frugiperda and mainly mediated by antiviral and amino acid metabolism-related pathways.

Funding

This work was supported by the Fund for Distinguished Young Scholars from the Henan Academy of Agricultural Sciences (grant no. 2020JQ05) and the National Natural Science Foundation of China (31772520; 31702057; and 31601897).

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  Capella-GutierrezS.Silla-MartinezJ. M.GabaldonT. (2009). Trimal: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics25, 1972–1973. doi: 10.1093/bioinformatics/btp348

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  ChenC.ChenH.ZhangY.ThomasH. R.FrankM. H.HeY.et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant13, 1194–1202. doi: 10.1016/j.molp.2020.06.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  CholletiH.HayerJ.FafetineJ.BergM.BlomstromA. L. (2018). Genetic characterization of a novel picorna-like virus in Culex spp. mosquitoes from Mozambique. Virol. J.15, 1–10. doi: 10.1186/s12985-018-0981-z

  • CrossRef
  • Google Scholar

• 4

  CordesE. J.Licking-MurrayK. D.CarlsonK. A. (2013). Differential gene expression related to Nora virus infection of Drosophila melanogaster. Virus Res.175, 95–100. doi: 10.1016/j.virusres.2013.03.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  de MirandaJ. R.DainatB.LockeB.CordoniG.BerthoudH.GauthierL.et al. (2010). Genetic characterization of slow bee paralysis virus of the honeybee (Apis mellifera L.). J. Gen. Virol.91, 2524–2530. doi: 10.1099/vir.0.022434-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  DingF.ChengJ.FuY.ChenT.LiB.JiangD.et al. (2019). Early transcriptional response to DNA virus infection in Sclerotinia sclerotiorum. Viruses11:278. doi: 10.3390/v11030278

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  DreherT. W.MillerW. A. (2006). Translational control in positive strand RNA plant viruses. Virology344, 185–197. doi: 10.1016/j.virol.2005.09.031

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  DuZ. Q.JinY. H. (2017). Comparative transcriptome and potential antiviral signaling pathways analysis of the gills in the red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus (WSSV). Genet. Mol. Biol.40, 168–180. doi: 10.1590/1678-4685-gmb-2016-0133

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  EkstromJ. O.HabayebM. S.SrivastavaV.KieselbachT.WingsleG.HultmarkD. (2011). Drosophila Nora virus capsid proteins differ from those of other picorna-like viruses. Virus Res.160, 51–58. doi: 10.1016/j.virusres.2011.05.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  GuoR.WangS.XueR.CaoG.HuX.HuangM.et al. (2015). The gene expression profile of resistant and susceptible Bombyx mori strains reveals cypovirus-associated variations in host gene transcript levels. Appl. Microbiol. Biotechnol.99, 5175–5187. doi: 10.1007/s00253-015-6634-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  HabayebM. S.CanteraR.CasanovaG.EkströmJ.-O.AlbrightS.HultmarkD. (2009a). The drosophila Nora virus is an enteric virus, transmitted via feces. J. Invertebr. Pathol.101, 29–33. doi: 10.1016/j.jip.2009.02.003

  • CrossRef
  • Google Scholar

• 12

  HabayebM. S.EkengrenS. K.HultmarkD. (2006). Nora virus, a persistent virus in drosophila, defines a new picorna-like virus family. J. Gen. Virol.87, 3045–3051. doi: 10.1099/vir.0.81997-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  HabayebM. S.EkströmJ. O.HultmarkD. (2009b). Nora virus persistent infections are not affected by the RNAi machinery. PLoS One4:e5731. doi: 10.1371/journal.pone.0005731

  • CrossRef
  • Google Scholar

• 14

  JakubowskaA. K.D'AngioloM.Gonzalez-MartinezR. M.Millan-LeivaA.CarballoA.MurilloR.et al. (2014). Simultaneous occurrence of covert infections with small RNA viruses in the lepidopteran Spodoptera exigua. J. Invertebr. Pathol.121, 56–63. doi: 10.1016/j.jip.2014.06.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  KalyaanamoorthyS.MinhB. Q.WongT. K. F.von HaeselerA.JermiinL. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods14, 587–589. doi: 10.1038/nmeth.4285

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods12, 357–360. doi: 10.1038/nmeth.3317

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  KumarS.StecherG.TamuraK. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol.33, 1870–1874. doi: 10.1093/molbev/msw054

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  LaurinmakiP.ShakeelS.EkstromJ. O.MohammadiP.HultmarkD.ButcherS. J. (2020). Structure of Nora virus at 2.7 A resolution and implications for receptor binding, capsid stability and taxonomy. Sci. Rep.10, 1–11. doi: 10.1038/s41598-020-76613-1

  • CrossRef
  • Google Scholar

• 19

  Le GallO.ChristianP.FauquetC. M.KingA. M.KnowlesN. J.NakashimaN.et al. (2008). Picornavirales, a proposed order of positive-sense single-stranded RNA viruses with a pseudo-T = 3 virion architecture. Arch. Virol.153, 715–727. doi: 10.1007/s00705-008-0041-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  LiuS. J.ChenY. T.BonningB. C. (2017). RNA virus discovery in insects. Curr. Opin. Insect. Sci.8, 54–61. doi: 10.1016/j.cois.2014.12.005

  • CrossRef
  • Google Scholar

• 21

  LiuW. Z.XieY. B.MaJ. Y.LuoX. T.NieP.ZuoZ. X.et al. (2015). IBS: an illustrator for the presentation and visualization of biological sequences. Bioinformatics31, 3359–3361. doi: 10.1093/bioinformatics/btv362

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  LopezW.PageA. M.CarlsonD. J.EricsonB. L.CserhatiM. F.GudaC.et al. (2018). Analysis of immune-related genes during Nora virus infection of Drosophila melanogaster using next generation sequencing. AIMS. Microbiol.4, 123–139. doi: 10.3934/microbiol.2018.1.123

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  MeddN. C.FellousS.WaldronF. M.XuerebA.NakaiM.CrossJ. V.et al. (2018). The virome of Drosophila suzukii, an invasive pest of soft fruit. Virus. Evol.4:vey009. doi: 10.1093/ve/vey009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  NguyenL. T.SchmidtH. A.von HaeselerA.MinhB. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol.32, 268–274. doi: 10.1093/molbev/msu300

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  OliveiraD. C.HunterW. B.NgJ.DesjardinsC. A.DangP. M.WerrenJ. H. (2010). Data mining cDNAs reveals three new single stranded RNA viruses in Nasonia (hymenoptera: Pteromalidae). Insect Mol. Biol.19(Suppl 1), 99–107. doi: 10.1111/j.1365-2583.2009.00934.x

  • CrossRef
  • Google Scholar

• 26

  PerteaM.PerteaG. M.AntonescuC. M.ChangT. C.MendellJ. T.SalzbergS. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol.33, 290–295. doi: 10.1038/nbt.3122

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  RemnantE. J.ShiM.BuchmannG.BlacquiereT.HolmesE. C.BeekmanM.et al. (2017). A diverse range of novel RNA viruses in geographically distinct honey bee populations. J. Virol.91, e00158–e00117. doi: 10.1128/JVI.00158-17

  • CrossRef
  • Google Scholar

• 28

  RobinsonM. D.McCarthyD. J.SmythG. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics26, 139–140. doi: 10.1093/bioinformatics/btp616

  • CrossRef
  • Google Scholar

• 29

  RogersA.ToweryL.McCownA.CarlsonK. A. (2020). Impaired geotaxis as a novel phenotype of Nora virus infection of Drosophila melanogaster. Scientifica2020:6. doi: 10.1155/2020/1804510

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  ShiM.LinX. D.ChenX.TianJ. H.ChenL. J.LiK.et al. (2018). The evolutionary history of vertebrate RNA viruses. Nature556, 197–202. doi: 10.1038/s41586-018-0012-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  ShiM.LinX. D.TianJ. H.ChenL. J.ChenX.LiC. X.et al. (2016). Redefining the invertebrate RNA virosphere. Nature540, 539–543. doi: 10.1038/nature20167

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  ShiW.XieS.ChenX.SunS.ZhouX.LiuL.et al. (2013). Comparative genomic analysis of the microbiome [corrected] of herbivorous insects reveals eco-environmental adaptations: biotechnology applications. PLoS Genet.9:e1003131. doi: 10.1371/annotation/91a25db3-8127-42c7-baa0-ce398a2857a6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  SunX. X.HuC. X.JiaH. R.WuQ. L.ShenX. J.ZhaoS. Y.et al. (2021). Case study on the first immigration of fall armyworm, Spodoptera frugiperda invading into China. J. Integr. Agr.20, 664–672. doi: 10.1016/S2095-3119(19)62839-X

  • CrossRef
  • Google Scholar

• 34

  van MierloJ. T.BronkhorstA. W.OverheulG. J.SadanandanS. A.EkstromJ. O.HeestermansM.et al. (2012). Convergent evolution of argonaute-2 slicer antagonism in two distinct insect RNA viruses. PLoS Path.8:e1002872. doi: 10.1371/journal.ppat.1002872

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  van MierloJ. T.OverheulG. J.ObadiaB.van CleefK. W.WebsterC. L.SalehM. C.et al. (2014). Novel drosophila viruses encode host-specific suppressors of RNAi. PLoS Path.10:e1004256. doi: 10.1371/journal.ppat.1004256

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  XieY.WuG.TangJ.LuoR.PattersonJ.LiuS.et al. (2014). SOAPdenovo-trans: de novo transcriptome assembly with short RNA-Seq reads. Bioinformatics30, 1660–1666. doi: 10.1093/bioinformatics/btu077

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  XuP.LiuY.GrahamR. I.WilsonK.WuK. (2014). Densovirus is a mutualistic symbiont of a global crop pest (Helicoverpa armigera) and protects against a baculovirus and Bt biopesticide. PLoS Path.10:e1004490. doi: 10.1371/journal.ppat.1004490

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  XuP.YangL.YangX.LiT.GrahamR. I.WuK.et al. (2020). Novel partiti-like viruses are conditional mutualistic symbionts in their normal lepidopteran host, African armyworm, but parasitic in a novel host Fall armyworm. PLoS Path.16:e1008467. doi: 10.1371/journal.ppat.1008467

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  YangY. T.NaiY. S.LeeS. J.LeeM. R.KimS.KimJ. S. (2016). A novel picorna-like virus, Riptortus pedestris virus-1 (RiPV-1), found in the bean bug, R. pedestris, after fungal infection. J. Invertebr. Pathol.141, 57–65. doi: 10.1016/j.jip.2016.11.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  YangX.XuP.YuanH.GrahamR. I.WilsonK.WuK. (2019). Discovery and characterization of a novel picorna-like RNA virus in the cotton bollworm Helicoverpa armigera. J. Invertebr. Pathol.160, 1–7. doi: 10.1016/j.jip.2018.11.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS16, 284–287. doi: 10.1089/omi.2011.0118

  • Pubmed Abstract
  • CrossRef
  • Google Scholar