Variants in the mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) infected with Oryctes rhinoceros nudivirus in oil palm and coconut plantations.

The CRB (coconut rhinoceros beetle) haplotype was classified into CRB-S and CRB-G, based on the presence of single nucleotide polymorphisms (SNPs) in the mitochondrial cox1 gene. Mitochondrial genomes (mitogenomes) are the most widely used genetic resources for molecular evolution, phylogenetics, and population genetics in relation to insects. This study presents the mitogenome CRB-G and CRB-S which were collected in Johor, Malaysia. The mitogenome of CRB-G collected from oil palm plantations in 2020 and 2021, and wild coconut palms in 2021 was 15,315 bp, 15,475 bp, and 17,275 bp, respectively. The CRB-S was discovered in coconut and oil palms in 2021, and its mitogenome was 15,484 bp and 17,142 bp, respectively. All the mitogenomes have 37 genes with more than 99% nucleotide sequence homology, except the CRB-G haplotype collected from oil palm in 2021 with 89.24% nucleotide sequence homology. The mitogenome of Johor CRBs was variable in the natural population due to its elevated mutation rate. Substitutions and indels in cox1, cox2, nad2 and atp6 genes were able to distinguish the Johor CRBs into two haplotypes. The mitogenome data generated in the present study may provide baseline information to study the infection and relationship between the two haplotypes of Johor CRB and OrNV in the field. This study is the first report on the mitogenomes of mixed haplotypes of CRB in the field.

Conserved Viral Transcription Plays a Key Role in Virus-Like Particle Production of the Parasitoid Wasp Venturia canescens.

Nudiviruses are large double-stranded DNA viruses related to baculoviruses known to be endogenized in the genomes of certain parasitic wasp species. These wasp-virus associations allow the production of viral particles or virus-like particles that ensure wasp parasitism success within lepidopteran hosts. Venturia canescens is an ichneumonid wasp belonging to the Campopleginae subfamily that has endogenized nudivirus genes belonging to the Alphanudivirus genus to produce "virus-like particles" (Venturia canescens virus-like particles [VcVLPs]), which package proteic virulence factors. The main aim of this study was to determine whether alphanudivirus gene functions have been conserved following endogenization. The expression dynamics of alphanudivirus genes was monitored by a high throughput transcriptional approach, and the functional role of lef-4 and lef-8 genes predicted to encode viral RNA polymerase components was investigated by RNA interference. As described for baculovirus infections and for endogenized nudivirus genes in braconid wasp species producing bracoviruses, a transcriptional cascade involving early and late expressed alphanudivirus genes could be observed. The expression of lef-4 and lef-8 was also shown to be required for the expression of alphanudivirus late genes allowing correct particle formation. Together with previous literature, the results show that endogenization of nudiviruses in parasitoid wasps has repeatedly led to the conservation of the viral RNA polymerase function, allowing the production of viruses or viral-like particles that differ in composition but enable wasp parasitic success. IMPORTANCE This study shows that endogenization of a nudivirus genome in a Campopleginae parasitoid wasp has led to the conservation, as for endogenized nudiviruses in braconid parasitoid wasps, of the viral RNA polymerase function, required for the transcription of genes encoding viral particles involved in wasp parasitism success. We also showed for the first time that RNA interference (RNAi) can be successfully used to downregulate gene expression in this species, a model in behavioral ecology. This opens the opportunity to investigate the function of genes involved in other traits important for parasitism success, such as reproductive strategies and host choice. Fundamental data acquired on gene function in Venturia canescens are likely to be transferable to other parasitoid wasp species used in biological control programs. This study also renders possible the investigation of other nudivirus gene functions, for which little data are available.

The Complete Genome of Chelonus insularis Reveals Dynamic Arrangement of Genome Components in Parasitoid Wasps That Produce Bracoviruses.

Bracoviruses (BVs) are endogenized nudiviruses in parasitoid wasps of the microgastroid complex (family Braconidae). Microgastroid wasps have coopted nudivirus genes to produce replication-defective virions that females use to transfer virulence genes to parasitized hosts. The microgastroid complex further consists of six subfamilies and ∼50,000 species but current understanding of BV gene inventories and organization primarily derives from analysis of two wasp species in the subfamily Microgastrinae (Microplitis demolitor and Cotesia congregata) that produce M. demolitor BV (MdBV) and C. congregata BV (CcBV). Notably, several genomic features of MdBV and CcBV remain conserved since divergence of M. demolitor and C. congregata ∼53 million years ago (MYA). However, it is unknown whether these conserved traits more broadly reflect BV evolution, because no complete genomes exist for any microgastroid wasps outside the Microgastrinae. In this regard, the subfamily Cheloninae is of greatest interest because it diverged earliest from the Microgastrinae (∼85 MYA) after endogenization of the nudivirus ancestor. Here, we present the complete genome of Chelonus insularis, which is an egg-larval parasitoid in the Cheloninae that produces C. insularis BV (CinsBV). We report that the inventory of nudivirus genes in C. insularis is conserved but are dissimilarly organized compared to M. demolitor and C. congregata. Reciprocally, CinsBV proviral segments share organizational features with MdBV and CcBV but virulence gene inventories exhibit almost no overlap. Altogether, our results point to the functional importance of a conserved inventory of nudivirus genes and a dynamic set of virulence genes for the successful parasitism of hosts. Our results also suggest organizational features previously identified in MdBV and CcBV are likely not essential for BV virion formation. IMPORTANCE Bracoviruses are a remarkable example of virus endogenization, because large sets of genes from a nudivirus ancestor continue to produce virions that thousands of wasp species rely upon to parasitize hosts. Understanding how these genes interact and have been coopted by wasps for novel functions is of broad interest in the study of virus evolution. This work characterizes bracovirus genome components in the parasitoid wasp Chelonus insularis, which together with existing wasp genomes captures a large portion of the diversity among wasp species that produce bracoviruses. Results provide new information about how bracovirus genome components are organized in different wasps while also providing additional insights on key features required for function.

Confirmation of Oryctes rhinoceros nudivirus infections in G-haplotype coconut rhinoceros beetles (Oryctes rhinoceros) from Palauan PCR-positive populations.

Coconut rhinoceros beetle (CRB), Oryctes rhinoceros, is a pest of palm trees in the Pacific. Recently, a remarkable degree of palm damage reported in Guam, Hawaii, Papua New Guinea and Solomon Islands has been associated with a particular haplotype (clade I), known as "CRB-G". In the Palau Archipelago, both CRB-G and another haplotype (clade IV) belonging to the CRB-S cluster coexist in the field. In this study, more than 75% of pheromone trap-captured adults of both haplotypes were Oryctes rhinoceros nudivirus (OrNV)-positive by PCR. No significant difference in OrNV prevalence between the haplotypes was detected. In PCR-positive CRB-G tissue specimens from Palau, viral particles were observed by electron microscopy. Hemocoel injection of CRB larvae with crude virus homogenates from these tissues resulted in viral infection and mortality. OrNV isolated from Palauan-sourced CRB was designated as OrNV-Palau1. Both OrNV-Palau1 and OrNV-X2B, a CRB biological control isolate released in the Pacific, were propagated using the FRI-AnCu-35 cell line for production of inoculum. However, the OrNV-Palau1 isolate exhibited lower viral production levels and longer larval survival times compared to OrNV-X2B in O. rhinoceros larvae. Full genome sequences of the OrNV-Palau1 and -X2B isolates were determined and found to be closely related to each other. Altogether these results suggest CRB adults in Palau are infected with a less virulent virus, which may affect the nature and extent of OrNV-induced pathology in Palauan populations of CRB.

Identification and Full Characterisation of Two Novel Crustacean Infecting Members of the Family Nudiviridae Provides Support for Two Subfamilies.

Multiple enveloped viruses with rod-shaped nucleocapsids have been described, infecting the epithelial cell nuclei within the hepatopancreas tubules of crustaceans. These bacilliform viruses share the ultrastructural characteristics of nudiviruses, a specific clade of viruses infecting arthropods. Using histology, electron microscopy and high throughput sequencing, we characterise two further bacilliform viruses from aquatic hosts, the brown shrimp (Crangon crangon) and the European shore crab (Carcinus maenas). We assembled the full double stranded, circular DNA genome sequences of these viruses (~113 and 132 kbp, respectively). Comparative genomics and phylogenetic analyses confirm that both belong within the family Nudiviridae but in separate clades representing nudiviruses found in freshwater and marine environments. We show that the three thymidine kinase (tk) genes present in all sequenced nudivirus genomes, thus far, were absent in the Crangon crangon nudivirus, suggesting there are twenty-eight core genes shared by all nudiviruses. Furthermore, the phylogenetic data no longer support the subdivision of the family Nudiviridae into four genera (Alphanudivirus to Deltanudivirus), as recently adopted by the International Committee on Taxonomy of Viruses (ICTV), but rather shows two main branches of the family that are further subdivided. Our data support a recent proposal to create two subfamilies within the family Nudiviridae, each subdivided into several genera.

Diverse Host Immune Responses of Different Geographical Populations of the Coconut Rhinoceros Beetle to Oryctes Rhinoceros Nudivirus (OrNV) Infection.

Incursions of the coconut rhinoceros beetle (CRB), Oryctes rhinoceros, into different islands in the South Pacific have been detected in recent years. It has been suggested that this range expansion is related to an O. rhinoceros haplotype reported to show reduced susceptibility to the well-established classical biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). Our understanding of the genetic characteristics which distinguish the population of O. rhinoceros that has recently established in Solomon Islands from other well-established populations across the region is very limited. Here, we hypothesized that the recently established O. rhinoceros population should have greater innate immune responses when challenged by OrNV than those of well-established and native O. rhinoceros populations. We used the RNA sequencing (RNA-Seq) approach to generate gene expression profiles of midgut tissue from OrNV-infected and noninfected individuals collected in the Solomon Islands (recent incursion), Papua New Guinea and Fiji (previously established), and the Philippines (within the native range). The collections included individuals from each of the three major mitochondrial lineages (CRB-G, CRB-PNG, and CRB-S) known to the region, allowing us to explore the specific responses of each haplotype to infection. Although insects from the Philippines and Solomon Islands that were tested belong to the same mitochondrial lineage (CRB-G), their overall responses to infection were different. The number of differentially expressed genes between OrNV-infected and noninfected wild-caught individuals from the four different locations varied from 148 to 252. Persistent OrNV infection caused a high level of induced antimicrobial activity and immune responses in O. rhinoceros, but the direction and magnitude of the responses were population specific. The insects tested from the Solomon Islands displayed extremely high expression of genes which are known to be involved in immune responses (e.g. coleoptericin, cecropin, and serpin). These variations in the host immune system among insects from different geographical regions might be driven by variations in the virulence of OrNV isolates, and this requires further investigation. Overall, our current findings support the importance of immunity in insect pest incursion and an expansion of the pest's geographic range. IMPORTANCE Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress coconut rhinoceros beetle (CRB) in the Pacific Islands. Recently a new wave of CRB incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles (CRB-G). Our comparative analysis of OrNV-infected and noninfected CRBs revealed that specific sets of genes were induced by viral infection in the beetles. This induction was much stronger in beetles collected from the Solomon Islands, a newly invaded country, than in individuals collected from within the beetle's native range (the Philippines) or from longer-established populations in its exotic range (Fiji and Papua New Guinea [PNG]). Beetles from the Philippines and the Solomon Islands that were tested in this study all belonged to the CRB-G haplotype, but the country-specific responses of the beetles to OrNV infection were different.

Nudivirus Sequences Identified from the Southern and Western Corn Rootworms (Coleoptera: Chrysomelidae).

Analysis of pooled genomic short read sequence data revealed the presence of nudivirus-derived sequences from U.S. populations of both southern corn rootworm (SCR, Diabrotica undecimpunctata howardi Barber) and western corn rootworm (WCR, Diabrotica virgifera virgifera LeConte). A near complete nudivirus genome sequence was assembled from sequence data for an SCR population with relatively high viral titers. A total of 147,179 bp was assembled from five contigs that collectively encode 109 putative open reading frames (ORFs) including 20 nudivirus core genes. In contrast, genome sequence recovery was incomplete for a second nudivirus from WCR, although sequences derived from this virus were present in three geographically dispersed populations. Only 48,989 bp were assembled with 48 putative ORFs including 13 core genes, representing about 20% of a typical nudivirus genome. Phylogenetic analysis indicated that both corn rootworm nudiviruses grouped with the third known nudivirus of beetles, Oryctes rhinoceros nudivirus in the genus Alphanudivirus. On the basis of phylogenetic and additional analyses, we propose further taxonomic separation of nudiviruses within Alphanudivirus and Betanudivirus into two subfamilies and five genera. Identification of nudivirus-derived sequences from two species of corn rootworm highlights the diversity of viruses associated with these agricultural insect pests.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

Transcription Profile and Genomic Variations of Oryctes Rhinoceros Nudivirus in Coconut Rhinoceros Beetles.

Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress the coconut rhinoceros beetle (Oryctes rhinoceros) in Southeast Asia and the Pacific Islands. A new wave of O. rhinoceros incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. In this study, chronically infected beetles were collected from Philippines, Fiji, Papua New Guinea (PNG), and the Solomon Islands (SI). RNA sequencing (RNA-seq) was performed to investigate the global viral gene expression profiles and for comparative genomic analysis of structural variations. Maximum likelihood phylogenic analysis indicated that OrNV strains from the SI and Philippines are closely related, while OrNV strains from PNG and Fiji formed a distinct adjacent clade. We detected several polymorphic sites with a frequency higher than 35% in 892 positions of the viral genome. Nonsynonymous mutations were detected in several hypothetical proteins and 15 nudivirus core genes, such as gp034, lef-8, lef-4, and vp91 We found limited evidence of variation in viral gene expression among geographic populations. Only a few genes, such as gp01, gp022, and gp107, were differentially expressed among different strains. Additionally, small RNA sequencing from the SI population suggested that OrNV is targeted by the host RNA interference (RNAi) response with abundant 21-nucleotide small RNAs. Some of these genomic changes are specific to the geographic population and could be related to particular phenotypic characteristics of the strain, such as viral pathogenicity or transmissibility, and this requires further investigation.IMPORTANCE Oryctes rhinoceros nudivirus has been an effective biocontrol agent against the coconut rhinoceros beetle in Southeast Asia and the Pacific Islands for decades. The recent outbreak of these beetles in many South Pacific islands has had a significant impact on livelihoods in the region. It has been suggested that the resurgence and spread of the pest are related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. We examined viral genomic and transcriptional variations in chronically infected beetles from different geographical populations. A high number of polymorphic sites among several geographical strains of OrNV were identified, but potentially only a few of these variations in the genome are involved in functional changes and can potentially alter the typical function. These findings provide valuable resources for future studies to improve our understanding of the OrNV genetic variations in different geographic regions and their potential link to virus pathogenicity.

Endogenous viral elements integrated into the genome of the soybean aphid, Aphis glycines.

Sequence analysis of the genomic DNA isolated from four biotypes of the soybean aphid, Aphis glycines (AG), revealed that in addition to the commonly observed retrovirus-related retrotransposons, viral sequences derived from multiple RNA and DNA viruses have integrated into the genome. Notably, sequences of more than 60 nudiviral genes were identified from de novo assembled DNA contigs, and mapped to assembled genomic scaffolds of AG, indicating that an ancient nudivirus, named Aphis glycines endogenous nudivirus (AgENV), had integrated into the AG genome. Furthermore, sequences derived from a similar endogenous nudivirus, Melanaphis sacchari endogenous nudivirus (MsENV), were identified from the genomic scaffolds of the sugarcane aphid, Melanaphis sacchari. Analysis of transcriptome and small RNA sequence data derived from AG did not provide evidence for transcription of the integrated AgENV genes. Hence, the genes of AgENV may be present as pseudogenes. Phylogenetic analysis based on nudivirus core genes indicated that these aphid ENVs belong to the genus Alphanudivirus.

Chalcid wasp paleoviruses bridge the evolutionary gap between bracoviruses and nudiviruses.

Polydnaviruses are obligate mutualists of parasitoid wasps and are divided into two genera, Bracovirus and Ichnovirus. Bracoviruses are thought to originate from a single integration of an ancestral nudivirus into the ancestor of microgastroid complex ~100 million years ago. However, all the known nudiviruses are only distantly related to bracoviruses, and much remains obscure about the origin of bracoviruses. Here we employ a paleovirological method to screen endogenous nudivirus-like elements across arthropods. Interestingly, we identify many endogenous nudivirus-like elements within the genome of Eurytoma brunniventris, a species of the Chalcidoidea superfamily. Among them, we find 14 core gene sequences are likely to be derived from a betanudivirus (designated EbrENV-ß), suggesting that betanudivirus has been circulating in parasitoid wasps. Phylogenomic analysis suggests that EbrENV-ß is the known closest relative of bracoviruses. Synteny analyses show the order of core genes is not well conserved between EbrENV-ß and nudiviruses, revealing the dynamic nature of the evolution of nudivirus genome structures. Our findings narrow down the evolutionary gap between bracoviruses and nudiviruses and provide novel insights into the origin and evolution of polydnaviruses.

ICTV Virus Taxonomy Profile: Nudiviridae.

Members of the family Nudiviridae are large dsDNA viruses with distinctive rod-shaped nucleocapsids and circular genomes of 96-232 kbp. Nudiviruses have been identified from a diverse range of insects and crustaceans and are closely related to baculoviruses. This is a summary of the International Committee on Taxonomy of Viruses Report on the taxonomy of the family Nudiviridae, which is available at ictv.global/report/nudiviridae.

Complete genome sequence of Oryctes rhinoceros nudivirus isolated from the coconut rhinoceros beetle in Solomon Islands.

Oryctes rhinoceros nudivirus (OrNV) has been an effective biocontrol agent against the insect pest Oryctes rhinoceros (Coleoptera: Scarabaeidae) for decades, but there is evidence that resistance could be evolving in some host populations. We detected OrNV infection in O. rhinoceros from Solomon Islands and used Oxford Nanopore Technologies (ONT) long-read sequencing to determine the full length of the virus genomic sequence isolated from an individual belonging to a mitochondrial lineage (CRB-G) that was previously reported as resistant to OrNV. The complete circular genome of the virus consisted of 125,917 nucleotides, 1.698 bp shorter than the originally-described full genome sequence of Ma07 strain from Malaysia. We found 130 out of 139 previously annotated ORFs (seven contained interrupted/non-coding sequences, two were identified as duplicated versions of the existing genes), as well as a putatively inverted regions containing four genes. These results demonstrate the usefulness of a long-read sequencing technology for resolving potential structural variations when describing new virus isolates. While the Solomon Islands isolate exhibited 99.41 % nucleotide sequence identity with the originally described strain, we found several genes, including a core gene (vlf-1), that contained multiple amino acid insertions and/or deletions as putative polymorphisms of large effect. Our complete annotated genome sequence of a newly found isolate in Solomon Islands provides a valuable resource to help elucidate the mechanisms that compromise the efficacy of OrNV as a biocontrol agent against the coconut rhinoceros beetle.

A Deep Learning Approach for Detecting Copy Number Variation in Next-Generation Sequencing Data.

Copy number variants (CNV) are associated with phenotypic variation in several species. However, properly detecting changes in copy numbers of sequences remains a difficult problem, especially in lower quality or lower coverage next-generation sequencing data. Here, inspired by recent applications of machine learning in genomics, we describe a method to detect duplications and deletions in short-read sequencing data. In low coverage data, machine learning appears to be more powerful in the detection of CNVs than the gold-standard methods of coverage estimation alone, and of equal power in high coverage data. We also demonstrate how replicating training sets allows a more precise detection of CNVs, even identifying novel CNVs in two genomes previously surveyed thoroughly for CNVs using long read data.

The first clawed lobster virus Homarus gammarus nudivirus (HgNV n. sp.) expands the diversity of the Nudiviridae.

Viral diseases of crustaceans are increasingly recognised as challenges to shellfish farms and fisheries. Here we describe the first naturally-occurring virus reported in any clawed lobster species. Hypertrophied nuclei with emarginated chromatin, characteristic histopathological lesions of DNA virus infection, were observed within the hepatopancreatic epithelial cells of juvenile European lobsters (Homarus gammarus). Transmission electron microscopy revealed infection with a bacilliform virus containing a rod shaped nucleocapsid enveloped in an elliptical membrane. Assembly of PCR-free shotgun metagenomic sequencing produced a circular genome of 107,063 bp containing 97 open reading frames, the majority of which share sequence similarity with a virus infecting the black tiger shrimp: Penaeus monodon nudivirus (PmNV). Multiple phylogenetic analyses confirm the new virus to be a novel member of the Nudiviridae: Homarus gammarus nudivirus (HgNV). Evidence of occlusion body formation, characteristic of PmNV and its closest relatives, was not observed, questioning the horizontal transmission strategy of HgNV outside of the host. We discuss the potential impacts of HgNV on juvenile lobster growth and mortality and present HgNV-specific primers to serve as a diagnostic tool for monitoring the virus in wild and farmed lobster stocks.

Common themes in three independently derived endogenous nudivirus elements in parasitoid wasps.

Endogenous Viral Elements (EVEs) are remnants of viral genomes that are permanently integrated into the genome of another organism. Parasitoid wasps have independently acquired nudivirus-derived EVEs in three lineages. Each parasitoid produces virions or virus-like particles (VLPs) that are injected into hosts during parasitism to function in subversion of host defenses. Comparing the inventory of nudivirus-like genes in different lineages of parasitoids can provide insights into the importance of each encoded function in virus or VLP production and parasitism success. Comparisons revealed the following conserved features: first, retention of genes encoding a viral RNA polymerase and infectivity factors; second, loss of the ancestral DNA polymerase gene; and third, signatures of viral ancestry in patterns of gene retention.