Varroa destructor shapes the unique viral landscapes of the honey bee populations of the Azores archipelago.

The worldwide dispersal of the ectoparasitic mite Varroa destructor from its Asian origins has fundamentally transformed the relationship of the honey bee (Apis mellifera) with several of its viruses, via changes in transmission and/or host immunosuppression. The extent to which honey bee-virus relationships change after Varroa invasion is poorly understood for most viruses, in part because there are few places in the world with several geographically close but completely isolated honey bee populations that either have, or have not, been exposed long-term to Varroa, allowing for separate ecological, epidemiological, and adaptive relationships to develop between honey bees and their viruses, in relation to the mite's presence or absence. The Azores is one such place, as it contains islands with and without the mite. Here, we combined qPCR with meta-amplicon deep sequencing to uncover the relationship between Varroa presence, and the prevalence, load, diversity, and phylogeographic structure of eight honey bee viruses screened across the archipelago. Four viruses were not detected on any island (ABPV-Acute bee paralysis virus, KBV-Kashmir bee virus, IAPV-Israeli acute bee paralysis virus, BeeMLV-Bee macula-like virus); one (SBV-Sacbrood virus) was detected only on mite-infested islands; one (CBPV-Chronic bee paralysis virus) occurred on some islands, and two (BQCV-Black queen cell virus, LSV-Lake Sinai virus,) were present on every single island. This multi-virus screening builds upon a parallel survey of Deformed wing virus (DWV) strains that uncovered a remarkably heterogeneous viral landscape featuring Varroa-infested islands dominated by DWV-A and -B, Varroa-free islands naïve to DWV, and a refuge of the rare DWV-C dominating the easternmost Varroa-free islands. While all four detected viruses investigated here were affected by Varroa for one or two parameters (usually prevalence and/or the Richness component of ASV diversity), the strongest effect was observed for the multi-strain LSV. Varroa unambiguously led to elevated prevalence, load, and diversity (Richness and Shannon Index) of LSV, with these results largely shaped by LSV-2, a major LSV strain. Unprecedented insights into the mite-virus relationship were further gained from implementing a phylogeographic approach. In addition to enabling the identification of a novel LSV strain that dominated the unique viral landscape of the easternmost islands, this approach, in combination with the recovered diversity patterns, strongly suggests that Varroa is driving the evolutionary change of LSV in the Azores. This study greatly advances the current understanding of the effect of Varroa on the epidemiology and adaptive evolution of these less-studied viruses, whose relationship with Varroa has thus far been poorly defined.

The hidden RNA viruses in Blattodea (cockroaches and termites).

The insect order Blattodea (cockroaches and termites) has drawn substantial research attention for their dietary habits and lifestyle of living with or around humans. In the present study, we focused on the discovery of RNA viruses hidden in Blattodea insects using the publicly available RNA sequencing datasets. Overall, 136 distinctive RNA viruses were identified from 36 Blattodea species, of which more than 70 % were most closely related to the invertebrate-associated viral groups within Picornavirales, Sobelivirales, Bunyaviricetes, Jingchuvirales, Durnavirales, Lispiviridae, Orthomyxoviridae, Permutotetraviridae, Flaviviridae and Muvirales. Several viruses were associated with pathogens of vertebrates (Paramyxoviridae), plants (Tymovirales), protozoa (Totiviridae), fungi (Narnaviridae) and bacteria (Norzivirales). Collectively, 93 complete or near-complete viral genomes were retrieved from the datasets, and several viruses appeared to have remarkable temporal and spatial distributions. Interestingly, the newly identified Periplaneta americana dicistrovirus displayed a remarkable distinct bicistronic genome arrangement from the well-recognized dicistroviruses with the translocated structural and non-structural polyprotein encoding open reading frames over the genome. These results significantly enhance our knowledge of RNA virosphere in Blattodea insects, and the novel genome architectures in dicistroviruses and other RNA viruses may break our stereotypes in the understanding of the genomic evolution and the emergence of potential novel viral species.

Identification and Genome Characterization of a Novel Nege-like Virus Isolated from Aphids (Aphis gossypii) in Yunnan Province.

Negeviruses are insect-specific enveloped RNA viruses that exhibit a wide geographic distribution. A novel nege-like virus, tentatively named Aphis gossypii nege-like virus (AGNLV, GenBank: OR880429.1), was isolated from aphids (Aphis gossypii) in Lijiang City, Yunnan, China. AGNLV has a genome sequence of 9258 nt (excluding the polyA tail) encoding three open reading frames (ORFs). ORF1 (7149 nt) encodes a viral methyltransferase, a viral RNA helicase, and an RNA-dependent RNA polymerase. ORF2 (1422 nt) encodes a DiSB-ORF2_chro domain and ORF3 encodes an SP24 domain. The genome sequence of AGNLV shares the highest nucleotide identity of 60.0% and 59.5% with Wuhan house centipede virus 1 (WHCV1) and Astegopteryx formosana nege-like virus (AFNLV), respectively. Phylogenetic analysis based on the RNA-dependent RNA polymerase shows that AGNLV is clustered with other negeviruses and nege-like viruses discovered in aphids, forming a distinct "unclassified clade". Interestingly, AGNLV only encodes three ORFs, whereas AFNLV and WHCV1 have four ORFs. Structure and transmembrane domain predictions show the presence of eight alpha helices and five transmembrane helices in the AGNLV ORF3. Translational enhancement of the AGNLV 5' UTR was similar to that of the 5' UTR of plant viruses. Our findings provide evidence of the diversity and structure of nege-like viruses and are the first record of such a virus from a member of the genus Aphis.

Known and novel viruses in Belgian honey bees: yearly differences, spatial clustering, and associations with overwintering loss.

In recent years, managed honey bee colonies have been suffering from an increasing number of biotic and abiotic stressors, resulting in numerous losses of colonies worldwide. A pan-European study, EPILOBEE, estimated the colony loss in Belgium to be 32.4% in 2012 and 14.8% in 2013. In the current study, absolute viral loads of four known honey bee viruses (DWV-A, DWV-B, AmFV, and BMLV) and three novel putative honey bee viruses (Apis orthomyxovirus 1, apthili virus, and apparli virus) were determined in 300 Flemish honey bee samples, and associations with winter survival were determined. This revealed that, in addition to the known influence of DWV-A and DWV-B on colony health, one of the newly described viruses (apthili virus) shows a strong yearly difference and is also associated with winter survival. Furthermore, all scrutinized viruses revealed significant spatial clustering patterns, implying that despite the limited surface area of Flanders, local virus transmission is paramount. The vast majority of samples were positive for at least one of the seven investigated viruses, and up to 20% of samples were positive for at least one of the three novel viruses. One of those three, Apis orthomyxovirus 1, was shown to be a genuine honey bee-infecting virus, able to infect all developmental stages of the honey bee, as well as the Varroa destructor mite. These results shed light on the most prevalent viruses in Belgium and their roles in the winter survival of honey bee colonies. IMPORTANCE: The western honey bee (Apis mellifera) is a highly effective pollinator of flowering plants, including many crops, which gives honey bees an outstanding importance both ecologically and economically. Alarmingly high annual loss rates of managed honey bee colonies are a growing concern for beekeepers and scientists and have prompted a significant research effort toward bee health. Several detrimental factors have been identified, such as varroa mite infestation and disease from various bacterial and viral agents, but annual differences are often not elucidated. In this study, we utilize the viral metagenomic survey of the EPILOBEE project, a European research program for bee health, to elaborate on the most abundant bee viruses of Flanders. We complement the existing metagenomic data with absolute viral loads and their spatial and temporal distributions. Furthermore, we identify Apis orthomyxovirus 1 as a potentially emerging pathogen, as we find evidence for its active replication honey bees.

Landscape structure affects temporal dynamics in the bumble bee virome: Landscape heterogeneity supports colony resilience.

Virus spillovers from managed honey bees, Apis mellifera, are thought to contribute to the decline of wild pollinators, including bumble bees. However, data on the impact of such viruses on wild pollinators remain scarce, and the influence of landscape structure on virus dynamics is poorly understood. In this study, we deployed bumble bee colonies in an agricultural landscape and studied changes in the bumble bee virome during field placement under varying habitat composition and configuration using a multiscale analytical framework. We estimated prevalence of viruses and viral loads (i.e. number of viral genomic equivalent copies) in bumble bees before and after placing them in the field using next generation sequencing and quantitative PCR. The results show that viral loads and number of different viruses present increased during placement in the field and that the virus composition of the colonies shifted from an initial dominance of honey bee associated viruses to a higher number (in both viral loads and number of viruses present) of bumble bee associated viruses. Especially DWV-B, typical for honey bees, drastically decreased after the time in the field. Viral loads prior to placing colonies in the field showed no effect on colony development, suggesting low impacts of these viruses in field settings. Notably, we further demonstrate that increased habitat diversity results in a lower number of different viruses present in Bombus colonies, while colonies in areas with well-connected farmland patches decreased in their total viral load after field placement. Our results emphasize the importance of landscape heterogeneity and connectivity for wild pollinator health and that these influences predominate at fine spatial scales.

Naturally occurring viruses of Drosophila reduce offspring number and lifespan.

Drosophila remains a pre-eminent insect model system for host-virus interaction, but the host range and fitness consequences of the drosophilid virome are poorly understood. Metagenomic studies have reported approximately 200 viruses associated with Drosophilidae, but few isolates are available to characterize the Drosophila immune response, and most characterization has relied on injection and systemic infection. Here, we use a more natural infection route to characterize the fitness effects of infection and to study a wider range of viruses. We exposed laboratory Drosophila melanogaster to 23 naturally occurring viruses from wild-collected drosophilids. We recorded transmission rates along with two components of female fitness: survival and the lifetime number of adult offspring produced. Nine different viruses transmitted during contact with laboratory D. melanogaster, although for the majority, rates of transmission were less than 20%. Five virus infections led to a significant decrease in lifespan (D. melanogaster Nora virus, D. immigrans Nora virus, Muthill virus, galbut virus and Prestney Burn virus), and three led to a reduction in the total number of offspring. Our findings demonstrate the utility of the Drosophila model for community-level studies of host-virus interactions, and suggest that viral infection could be a substantial fitness burden on wild flies.

Development of a multiplex real-time quantitative reverse-transcription polymerase chain reaction for the detection of four bee viruses.

Viruses in the families Dicistroviridae and Iflaviridae are among the main threats to western honey bees (Apis mellifera) and native bee species. Polymerase chain reaction (PCR) is the gold standard for pathogen detection in bees. However, high throughput screening for bee virus infections in singleplex PCR reactions is cumbersome and limited by the high quantities of sample RNA required. Thus, the development of a sensitive and specific multiplex PCR detection method for screening for multiple viruses simultaneously is necessary. Here, we report the development of a one-step multiplex reverse-transcription quantitative polymerase chain reaction (RT-qPCR) assay to detect four viruses commonly encountered in pollinator species. The optimized multiplex RT-qPCR protocol described in this study allows simultaneous detection of two dicistroviruses (Israeli acute paralysis virus and Black queen cell virus) and two iflaviruses (Sacbrood virus and Deformed wing virus) with high efficiency and specificity comparable to singleplex detection assays. This assay provides a broad range of detection and quantification, and the results of virus quantification in this study are similar to those performed in other studies using singleplex detection assays. This method will be particularly useful for data generation from small-bodied insect species that yield low amounts of RNA.

Lack of abundant core virome in Culex mosquitoes from a temperate climate region despite a mosquito species-specific virome.

In arthropod-associated microbial communities, insect-specific viruses (ISVs) are prevalent yet understudied due to limited infectivity outside their natural hosts. However, ISVs might play a crucial role in regulating mosquito populations and influencing arthropod-borne virus transmission. Some studies have indicated a core virome in mosquitoes consisting of mostly ISVs. Employing single mosquito metagenomics, we comprehensively profiled the virome of native and invasive mosquito species in Belgium. This approach allowed for accurate host species determination, prevalence assessment of viruses and Wolbachia, and the identification of novel viruses. Contrary to our expectations, no abundant core virome was observed in Culex mosquitoes from Belgium. In that regard, we caution against rigidly defining mosquito core viromes and encourage nuanced interpretations of other studies. Nonetheless, our study identified 45 viruses of which 28 were novel, enriching our understanding of the mosquito virome and ISVs. We showed that the mosquito virome in this study is species-specific and less dependent on the location where mosquitoes from the same species reside. In addition, because Wolbachia has previously been observed to influence arbovirus transmission, we report the prevalence of Wolbachia in Belgian mosquitoes and the detection of several Wolbachia mobile genetic elements. The observed prevalence ranged from 83% to 92% in members from the Culex pipiens complex.IMPORTANCECulex pipiens mosquitoes are important vectors for arboviruses like West Nile virus and Usutu virus. Virome studies on individual Culex pipiens, and on individual mosquitoes in general, have been lacking. To mitigate this, we sequenced the virome of 190 individual Culex and 8 individual Aedes japonicus mosquitoes. We report the lack of a core virome in these mosquitoes from Belgium and caution the interpretation of other studies in this light. The discovery of new viruses in this study will aid our comprehension of insect-specific viruses and the mosquito virome in general in relation to mosquito physiology and mosquito population dynamics.

Parasite and virus dynamics in the honeybee Apis mellifera unicolor on a tropical island recently invaded by Varroa destructor.

In La Réunion, the established honeybee subspecies Apis mellifera unicolor, an endemic subspecies of African lineage, is facing considerable challenges. Since the introduction of the Varroa destructor mite in 2017 high colony losses have been recorded. We investigated the dynamics of V. destructor and two viruses, the Deformed Wing Virus (DWV), known to be transmitted by the mite, and the Chronic Bee Paralysis Virus (CBPV), in A. m. unicolor. Colonies from two apiaries located at 300 and 900 m a.s.l were monitored twice for one year without any acaricide treatment. The brood area, V. destructor infestation rates, DWV and CBPV prevalence and load were recorded monthly. A. m. unicolor maintained brood rearing throughout the year. Varroa destructor infestation resulted in high colony mortality (up to 85 %) and high phoretic mite rates (up to 52 mites per hundred bees). The establishment of DWV in colonies occurred after that of V. destructor and the mite infestation rate had a significant effect on the virus prevalence and load. CBPV appeared only transiently throughout the surveys. The data showed that, in tropical colonies with permanent brood rearing, V. destructor and DWV can reach high levels, but are still subject to seasonal variations that appear to be influenced by environmental conditions. This suggests that beekeeping practices could be adapted by favouring sites and periods for transhumance or acaricide treatment.

Cell fusing agent virus rarely transmits vertically in artificially infected laboratory-colonized Aedes aegypti mosquitoes.

BACKGROUND: Vertical transmission (VT) of arboviruses (arthropod-borne viruses) can serve as an essential link in the transmission cycle during adverse environmental conditions. The extent of VT among mosquito-borne arboviruses can vary significantly among different virus families and even among different viruses within the same genus. For example, orthobunyaviruses exhibit a higher VT rate than orthoflaviviruses and alphaviruses. Mosquitoes are also the natural hosts of a large number of insect-specific viruses (ISV) that belong to several virus families, including Bunyaviridae, Flaviviridae, and Togaviridae. Cell fusing agent virus (CFAV), an insect-specific orthoflavivirus, displays higher VT rates than other dual-host orthoflaviviruses, such as Zika and dengue viruses. High VT rates require establishment of stabilized infections in the germinal tissues of female vectors. To delve deeper into understanding the mechanisms governing these differences in VT rates and the establishment of stabilized infections, the ovary infection patterns and VT of Zika virus (ZIKV) and CFAV were compared. METHODS: Laboratory colonized Aedes aegypti females were infected with either ZIKV or CFAV by intrathoracic injection. Ovary infection patterns were monitored by in situ hybridization using virus-specific probes, and VT was determined by detecting the presence of the virus among the progeny, using a reverse-transcription quantitative polymerase chain reaction (PCR) assay. RESULTS: Both ZIKV and CFAV infect mosquito ovaries after intrathoracic injection. Infections then become widespread following a non-infectious blood meal. VT rates of ZIKV are similar to previously reported results (3.33%). CFAV, on the contrary transmits vertically very rarely. VT was not observed in the first gonotrophic cycle following intrathoracic injection, and only rarely in the second gonotrophic cycle. VT of CFAV is mosquito population independent, since similar results were obtained with Aedes aegypti collected from two different geographic locations. CONCLUSIONS: Although CFAV infects mosquito ovaries, the occurrence of VT remains infrequent in artificially infected Ae. aegypti, despite the observation of high VT rates in field-collected mosquitoes. These results suggest that infections of insect-specific viruses are stabilized in mosquitoes by some as yet unidentified mechanisms.

Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing.

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Transcriptome Analysis Reveals a Diverse Range of Novel Viruses in Australian Sugarcane Soldier Fly (Inopus flavus) Larvae.

In Australia, Soldier flies (Inopus spp.) are economically significant pests of sugarcane that currently lack a viable management strategy. Despite various research efforts, the mechanisms underlying the damage caused by soldier fly larvae remain poorly understood. Our study aims to explore whether this damage is associated with the transmission of plant viruses during larval feeding. We also explore the larval transcriptome to identify any entomopathogenic viruses with the potential to be used as biocontrol agents in future pest management programs. Seven novel virus sequences are identified and characterised using de novo assembly of RNA-Seq data obtained from salivary glands of larvae. The novel virus sequences belong to different virus families and are tentatively named SF-associated anphevirus (SFaAV), SF-associated orthomyxo-like virus (SFaOV), SF-associated narna-like virus (SFaNV), SF-associated partiti-like virus (SFaPV), SF-associated toti-like virus (SFaTV-1 and SFaTV-2) and SF-associated densovirus (SFaDV). These newly identified viruses are more likely insect-associated viruses, as phylogenetic analyses show that they cluster with other insect-specific viruses. Small RNA analysis indicates prominent peaks at both 21 nt and 26-29 nt, suggesting the activation of host siRNA and piwiRNA pathways. Our study helps to improve understanding of the virome of soldier flies and could identify insect viruses for deployment in novel pest management strategies.

Chronic bee paralysis virus exploits host antimicrobial peptides and alters gut microbiota composition to facilitate viral infection.

The significance of gut microbiota in regulating animal immune response to viral infection is increasingly recognized. However, how chronic bee paralysis virus (CBPV) exploits host immune to disturb microbiota for its proliferation remains elusive. Through histopathological examination, we discovered that the hindgut harbored the highest level of CBPV, and displayed visible signs of damages. The metagenomic analysis showed that a notable reduction in the levels of Snodgrassella alvi and Lactobacillus apis, and a significant increase in the abundance of the opportunistic pathogens such as Enterobacter hormaechei and Enterobacter cloacae following CBPV infection. Subsequent co-inoculation experiments showed that these opportunistic pathogens facilitated the CBPV proliferation, leading to accelerated mortality in bees and exacerbation of bloated abdomen symptoms after CBPV infection. The expression level of antimicrobial peptide (AMP) was found to be significantly up-regulated by over 1000 times in response to CBPV infection, as demonstrated by subsequent transcriptome and quantitative real-time PCR investigations. In particular, through correlation analysis and a bacteriostatic test revealed that the AMPs did not exhibit any inhibitory effect against the two opportunistic pathogens. However, they did demonstrate inhibitory activity against S. alvi and L. apis. Our findings provide different evidence that the virus infection may stimulate and utilize the host's AMPs to eradicate probiotic species and facilitate the proliferation of opportunistic bacteria. This process weakens the intestinal barrier and ultimately resulting in the typical bloated abdomen.

Impact of symbiotic insect-specific viruses on mosquito vector competence for arboviruses.

Mosquitoes are vectors for arboviruses, such as dengue, Zika, and Chikungunya. Symbiotic interactions can affect the intrinsic ability of mosquitoes to acquire and transmit arboviruses, referred to as vector competence. Insect-specific viruses (ISVs) are commonly found in symbiotic associations with mosquitoes in the wild and can affect many aspects of mosquito biology. Here, we review current knowledge on the effects of symbiotic ISV-mosquito interactions on vector competence. We discuss potential mechanisms underlying these interactions and their implications for shaping new biological control strategies. Finally, we highlight the need for field data analyzing the circulation of ISVs in mosquitoes associated with mechanistic studies in the laboratory.

Negevirus Piura Suppresses Zika Virus Replication in Mosquito Cells.

We investigated the interaction between the insect-specific virus, Piura virus (PIUV), and the arbovirus Zika virus (ZIKV) in Aedes albopictus cells. We performed coinfection experiments in C6/36 cells. Piura virus (Cor 33 strain, Colombia) and ZIKV (PRVABC58 strain, Puerto Rico) were co-inoculated into C6/36 cells using two multiplicity of infection (MOI) combinations: 0.1 for both viruses and 1.0 for ZIKV, 0.1 for PIUV. Wells were infected in triplicate with either PIUV and ZIKV coinfection, ZIKV-only, or PIUV-only. Mock infected cells served as control wells. The cell suspension was collected daily 7 days post-infection. Zika virus load was titrated by TCID50 on Vero 76 cells. The ZIKV-only infection and PIUV and ZIKV coinfection experiments were also quantified by RT-qPCR. We also investigated whether ZIKV interfered in the PIUV replication. PIUV suppressed the replication of ZIKV, resulting in a 10,000-fold reduction in ZIKV titers within 3 days post-infection. PIUV viral loads were not reduced in the presence of ZIKV. We conclude that, when concurrently infected, PIUV suppresses ZIKV in C6/36 cells while ZIKV does not interfere in PIUV replication.

Nearly complete genome sequences of the first two identified Colorado potato beetle viruses.

The Colorado potato beetle is one of the most devastating potato pests in the world. However, its viral pathogens, which might have potential in pest control, have remained unexplored. With high-throughput sequencing of Colorado potato beetle samples derived from prepupal larvae which died from an unknown infection, we have identified two previously unknown RNA viruses and assembled their nearly complete genome sequences. The subsequent genetic and phylogenetic analysis demonstrated that the viruses, tentatively named Leptinotarsa iflavirus 1 and Leptinotarsa solinvi-like virus 1, are the novel representatives of the Iflaviridae and Solinviviridae families, respectively. To the best of our knowledge, these are the first sequencing-confirmed insect viruses derived from Colorado potato beetle samples. We propose that Leptinotarsa iflavirus 1 may be associated with a lethal disease in the Colorado potato beetle.

Identification of long noncoding RNAs of silkworm at the early stage of Bombyx mori bidensovirus infection.

Bombyx mori bidensovirus (BmBDV) is one of the most important pathogens of silkworm. It mainly infects midgut cells of silkworm and causes losses to the sericulture industry. Long noncoding RNAs (lncRNAs) have been reported to play an important role in the regulation of antiviral immune response in silkworm. To explore whether lncRNAs are involved in BmBDV infection and immune response of silkworm, we performed a comparative transcriptome analysis to identify the lncRNAs and mRNAs between the BmBDV infected and noninfected silkworm larvae at the early stage. A total of 16,069 genes and 974 candidate lncRNAs were identified, among which 142 messenger RNA (mRNAs) and four lncRNAs were differentially expressed (DE). Target gene prediction revealed that 142 DEmRNAs were coexpressed with four DElncRNAs, suggesting that the expression of mRNA is mainly affected through trans-regulation activities. A regulatory network of DElncRNAs and DEmRNAs was constructed, showing that many genes targeted by different DElncRNAs are involved in metabolism and immunity, which implies that these genes and lncRNAs play an important role in the replication of BmBDV. Our results will help us to improve our understanding of lncRNA-mediated regulatory roles in BmBDV infection, providing a new perspective for further exploring the interaction between host and BmBDV.

Insect-specific RNA virus affects the stylet penetration activity of brown citrus aphid (Aphis citricidus) to facilitate its transmission.

Sap-sucking insects often transmit plant viruses but also carry insect viruses, which infect insects but not plants. The impact of such insect viruses on insect host biology and ecology is largely unknown. Here, we identified a novel insect-specific virus carried by brown citrus aphid (Aphis citricidus), which we tentatively named Aphis citricidus picornavirus (AcPV). Phylogenetic analysis discovered a monophyletic cluster with AcPV and other unassigned viruses, suggesting that these viruses represent a new family in order Picornavirales. Systemic infection with AcPV triggered aphid antiviral immunity mediated by RNA interference, resulting in asymptomatic tolerance. Importantly, we found that AcPV was transmitted horizontally by secretion of the salivary gland into the feeding sites of plants. AcPV influenced aphid stylet behavior during feeding and increased the time required for intercellular penetration, thus promoting its transmission among aphids with plants as an intermediate site. The gene expression results suggested that this mechanism was linked with transcription of salivary protein genes and plant defense hormone signaling. Together, our results show that the horizontal transmission of AcPV in brown citrus aphids evolved in a manner similar to that of the circulative transmission of plant viruses by insect vectors, thus providing a new ecological perspective on the activity of insect-specific viruses found in aphids and improving the understanding of insect virus ecology.

Influence of dengue virus serotypes on the abundance of Aedes aegypti insect-specific viruses (ISVs).

A comprehensive understanding of the virome in mosquito vectors is crucial for assessing the potential transmission of viral agents, designing effective vector control strategies, and advancing our knowledge of insect-specific viruses (ISVs). In this study, we utilized Oxford Nanopore Technologies metagenomics to characterize the virome of Aedes aegypti mosquitoes collected in various regions of Colombia, a country hyperendemic for dengue virus (DENV). Analyses were conducted on groups of insects with previous natural DENV infection (DENV-1 and DENV-2 serotypes), as well as mosquito samples that tested negative for virus infection (DENV-negative). Our findings indicate that the Ae. aegypti virome exhibits a similar viral composition at the ISV family and species levels in both DENV-positive and DENV-negative samples across all study sites. However, differences were observed in the relative abundance of viral families such as Phenuiviridae, Partitiviridae, Flaviviridae, Rhabdoviridae, Picornaviridae, Bromoviridae, and Virgaviridae, depending on the serotype of DENV-1 and DENV-2. In addition, ISVs are frequently found in the core virome of Ae. aegypti, such as Phasi Charoen-like phasivirus (PCLV), which was the most prevalent and showed variable abundance in relation to the presence of specific DENV serotypes. Phylogenetic analyses of the L, M, and S segments of the PCLV genome are associated with sequences from different regions of the world but show close clustering with sequences from Brazil and Guadeloupe, indicating a shared evolutionary relationship. The profiling of the Ae. aegypti virome in Colombia presented here improves our understanding of viral diversity within mosquito vectors and provides information that opens the way to possible connections between ISVs and arboviruses. Future studies aimed at deepening our understanding of the mechanisms underlying the interactions between ISVs and DENV serotypes in Ae. aegypti could provide valuable information for the design of effective vector-borne viral disease control and prevention strategies.IMPORTANCEIn this study, we employed a metagenomic approach to characterize the virome of Aedes aegypti mosquitoes, with and without natural DENV infection, in several regions of Colombia. Our findings indicate that the mosquito virome is predominantly composed of insect-specific viruses (ISVs) and that infection with different DENV serotypes (DENV-1 and DENV-2) could lead to alterations in the relative abundance of viral families and species constituting the core virome in Aedes spp. The study also sheds light on the identification of the genome and evolutionary relationships of the Phasi Charoen-like phasivirus in Ae. aegypti in Colombia, a widespread ISV in areas with high DENV incidence.

Genetic diversity of RNA viruses infecting invertebrate pests of rice.

Invertebrate species are a natural reservoir of viral genetic diversity, and invertebrate pests are widely distributed in crop fields. However, information on viruses infecting invertebrate pests of crops is limited. In this report, we describe the deep metatranscriptomic sequencing of 88 invertebrate samples covering all major invertebrate pests in rice fields. We identified 296 new RNA viruses and 13 known RNA viruses. These viruses clustered within 31 families, with many highly divergent viruses constituting potentially new families and genera. Of the identified viruses, 13 RNA viruses clustered within the Fiersviridae family of bacteriophages, and 48 RNA viruses clustered within families and genera of mycoviruses. We detected known rice viruses in novel invertebrate hosts at high abundances. Furthermore, some novel RNA viruses have genome structures closely matching to known plant viruses and clustered within genera of several plant virus species. Forty-five potential insect pathogenic RNA viruses were detected in invertebrate species. Our analysis revealed that host taxonomy plays a major role and geographical location plays an important role in structuring viral diversity. Cross-species transmission of RNA viruses was detected between invertebrate hosts. Newly identified viral genomes showed extensive variation for invertebrate viral families or genera. Together, the large-scale metatranscriptomic analysis greatly expands our understanding of RNA viruses in rice invertebrate species, the results provide valuable information for developing efficient strategies to manage insect pests and virus-mediated crop diseases.