Wide spectrum and high frequency of genomic structural variation, including transposable elements, in large double-stranded DNA viruses

A qPCR Assay for the Quantification of Selected Genotypic Variants of Spodoptera frugiperda Multiple Nucleopolyhedrovirus (Baculoviridae)

Alphabaculoviruses are lethal dsDNA viruses of Lepidoptera that have high genetic diversity and are transmitted in aggregates within proteinaceous occlusion bodies. This mode of transmission has implications for their efficacy as biological insecticides. A Nicaraguan isolate of Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV-NIC) comprising nine genotypic variants has been the subject of considerable study due to the influence of variant interactions on the insecticidal properties of mixed-variant occlusion bodies. As part of a systematic study on the replication and transmission of variant mixtures, a tool for the accurate quantification of a selection of genotypic variants was developed based on the quantitative PCR technique (qPCR). First, primer pairs were designed around a region of high variability in four variants named SfNic-A, SfNic-B, SfNic-C and SfNic-E to produce amplicons of 103-150 bp. Then, using cloned purified amplicons as standards, amplification was demonstrated over a dynamic range of 108-101 copies of each target. The assay was efficient (mean ± SD: 98.5 ± 0.8%), reproducible, as shown by low inter- and intra-assay coefficients of variation (<5%), and specific to the target variants (99.7-100% specificity across variants). The quantification method was validated on mixtures of genotype-specific amplicons and demonstrated accurate quantification. Finally, mixtures of the four variants were quantified based on mixtures of budded virions and mixtures of DNA extracted from occlusion-derived virions. In both cases, mixed-variant preparations compared favorably to total viral genome numbers by quantification of the polyhedrin (polh) gene that is present in all variants. This technique should prove invaluable in elucidating the influence of variant diversity on the transmission and insecticidal characteristics of this pathogen.

Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

Spatially Segregated Transmission of Co-Occluded Baculoviruses Limits Virus-Virus Interactions Mediated by Cellular Coinfection during Primary Infection

The occlusion bodies (OBs) of certain alphabaculoviruses are polyhedrin-rich structures that mediate the collective transmission of tens of viral particles to the same insect host. In addition, in multiple nucleopolyhedroviruses, occlusion-derived virions (ODVs) form nucleocapsid aggregates that are delivered to the same host cell. It has been suggested that, by favoring coinfection, this transmission mode promotes evolutionarily stable interactions between different baculovirus variants. To quantify the joint transmission of different variants, we obtained OBs from cells coinfected with two viral constructs, each encoding a different fluorescent reporter, and used them for inoculating Spodoptera exigua larvae. The microscopy analysis of midguts revealed that the two reporter genes were typically segregated into different infection foci, suggesting that ODVs show limited ability to promote the co-transmission of different virus variants to the same host cell. However, a polyhedrin-deficient mutant underwent inter-host transmission by exploiting the OBs of a fully functional virus and re-acquired the lost gene through recombination, demonstrating cellular coinfection. Our results suggest that viral spatial segregation during transmission and primary infection limits interactions between different baculovirus variants, but that these interactions still occur within the cells of infected insects later in infection.

viralFlye: assembling viruses and identifying their hosts from long-read metagenomics data

Although the use of long-read sequencing improves the contiguity of assembled viral genomes compared to short-read methods, assembling complex viral communities remains an open problem. We describe the viralFlye tool for identification and analysis of metagenome-assembled viruses in long-read assemblies. We show it significantly improves viral assemblies and demonstrate that long-reads result in a much larger array of predicted virus-host associations as compared to short-read assemblies. We demonstrate that the identification of novel CRISPR arrays in bacterial genomes from a newly assembled metagenomic sample provides information for predicting novel hosts for novel viruses.

Nucleopolyhedrovirus Coocclusion Technology: A New Concept in the Development of Biological Insecticides

Nucleopolyhedroviruses (NPV, Baculoviridae) that infect lepidopteran pests have an established record as safe and effective biological insecticides. Here, we describe a new approach for the development of NPV-based insecticides. This technology takes advantage of the unique way in which these viruses are transmitted as collective infectious units, and the genotypic diversity present in natural virus populations. A ten-step procedure is described involving genotypic variant selection, mixing, coinfection and intraspecific coocclusion of variants within viral occlusion bodies. Using two examples, we demonstrate how this approach can be used to produce highly pathogenic virus preparations for pest control. As restricted host range limits the uptake of NPV-based insecticides, this technology has recently been adapted to produce custom-designed interspecific mixtures of viruses that can be applied to control complexes of lepidopteran pests on particular crops, as long as a shared host species is available for virus production. This approach to the development of NPV-based insecticides has the potential to be applied across a broad range of NPV-pest pathosystems.

Assessing the Impact of a Viral Infection on the Expression of Transposable Elements in the Cabbage Looper Moth (Trichoplusia ni)

Most studies of stress-induced transposable element (TE) expression have so far focused on abiotic sources of stress. Here, we analyzed the impact of an infection by the AcMNPV baculovirus on TE expression in a cell line (Tnms42) and midgut tissues of the cabbage looper moth (Trichoplusia ni). We find that a large fraction of TE families (576/636 in Tnms42 cells and 503/612 in midgut) is lowly expressed or not expressed at all [≤ 4 transcripts per million (TPM)] in the uninfected condition (median TPM of 0.37 in Tnms42 and 0.46 in midgut cells). In the infected condition, a total of 62 and 187 TE families were differentially expressed (DE) in midgut and Tnms42 cells, respectively, with more up- (46) than downregulated (16) TE families in the former and as many up- (91) as downregulated (96) TE families in the latter. Expression log2 fold changes of DE TE families varied from -4.95 to 9.11 in Tnms42 cells and from -4.28 to 7.66 in midgut. Large variations in expression profiles of DE TEs were observed depending on the type of cells and on time after infection. Overall, the impact of AcMNPV on TE expression in T. ni is moderate but potentially sufficient to affect TE activity and genome architecture. Interestingly, one host-derived TE integrated into AcMNPV genomes is highly expressed in infected Tnms42 cells. This result shows that virus-borne TEs can be expressed, further suggesting that they may be able to transpose and that viruses may act as vectors of horizontal transfer of TEs in insects.

Generation of Variability in Chrysodeixis includens Nucleopolyhedrovirus (ChinNPV): The Role of a Single Variant

The mechanisms generating variability in viruses are diverse. Variability allows baculoviruses to evolve with their host and with changes in their environment. We examined the role of one genetic variant of Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) and its contribution to the variability of the virus under laboratory conditions. A mixture of natural isolates (ChinNPV-Mex1) contained two genetic variants that dominated over other variants in individual larvae that consumed high (ChinNPV-K) and low (ChinNPV-E) concentrations of inoculum. Studies on the ChinNPV-K variant indicated that it was capable of generating novel variation in a concentration-dependent manner. In cell culture, cells inoculated with high concentrations of ChinNPV-K produced OBs with the ChinNPV-K REN profile, whereas a high diversity of ChinNPV variants was recovered following plaque purification of low concentrations of ChinNPV-K virion inoculum. Interestingly, the ChinNPV-K variant could not be recovered from plaques derived from low concentration inocula originating from budded virions or occlusion-derived virions of ChinNPV-K. Genome sequencing revealed marked differences between ChinNPV-K and ChinNPV-E, with high variation in the ChinNPV-K genome, mostly due to single nucleotide polymorphisms. We conclude that ChinNPV-K is an unstable genetic variant that is responsible for generating much of the detected variability in the natural ChinNPV isolates used in this study.

Paleovirology of the DNA viruses of eukaryotes

Paleovirology is the study of ancient viruses and how they have coevolved with their hosts. An increasingly detailed understanding of the diversity, origins, and evolution of the DNA viruses of eukaryotes has been obtained through the lens of paleovirology in recent years. Members of multiple viral families have been found integrated in the genomes of eukaryotes, providing a rich fossil record to study. These elements have extended our knowledge of exogenous viral diversity, host ranges, and the timing of viral evolution, and are revealing the existence of entire new families of eukaryotic integrating dsDNA viruses and transposons. Future work in paleovirology will continue to provide insights into antiviral immunity, viral diversity, and potential applications, and reveal other secrets of the viral world.

Monitoring insect transposable elements in large double-stranded DNA viruses reveals host-to-virus and virus-to-virus transposition

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

A holistic perspective on herpes simplex virus (HSV) ecology and evolution

Herpes simplex viruses (HSV) cause chronic infection in humans that are characterized by periodic episodes of mucosal shedding and ulcerative disease. HSV causes millions of infections world-wide, with lifelong bouts of viral reactivation from latency in neuronal ganglia. Infected individuals experience different levels of disease severity and frequency of reactivation. There are two distantly related HSV species, with HSV-1 infections historically found most often in the oral niche and HSV-2 infections in the genital niche. Over the last two decades, HSV-1 has emerged as the leading cause of first-episode genital herpes in multiple countries. While HSV-1 has the highest level of genetic diversity among human alpha-herpesviruses, it is not yet known how quickly the HSV-1 viral population in a human host adapts over time, or if there are population bottlenecks associated with viral reactivation and/or transmission. It is also unknown how the ecological environments in which HSV infections occur influence their evolutionary trajectory, or that of co-occurring viruses and microbes. In this review, we explore how HSV accrues genetic diversity within each new infection, and yet maintains its ability to successfully infect most of the human population. A holistic examination of the ecological context of natural human infections can expand our awareness of how HSV adapts as it moves within and between human hosts, and reveal the complexity of these lifelong human-virus interactions. These insights may in turn suggest new areas of exploration for other chronic pathogens that successfully evolve and persist among their hosts.

Genomic diversity in a population of Spodoptera frugiperda nucleopolyhedrovirus

Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) represents a strong candidate to develop environmental-friendly pesticides against the fall armyworm (Spodoptera frugiperda), a widespread pest that poses a severe threat to different crops around the world. To date, SfMNPV genomic diversity of different isolates has been mainly studied by means of restriction pattern analyses and by sequencing of the egt region. Here, the genomic diversity present inside an isolate of SfMNPV was explored using high-throughput sequencing for the first time. We identified 704 intrahost single nucleotide variants, from which 184 are nonsynonymous mutations distributed among 82 different coding sequences. We detected several structural variants affecting SfMNPV genome, including two previously reported deletions inside the egt region. A comparative analysis between polymorphisms present in different SfMNPV isolates and our intraisolate diversity data suggests that coding regions with higher genetic diversity are associated with oral infectivity or unknown functions. In this context, through molecular evolution studies we provide evidence of diversifying selection acting on sf29, a putative collagenase which could contribute to the oral infectivity of SfMNPV. Overall, our results contribute to deepen our understanding of the coevolution between SfMNPV and the fall armyworm and will be useful to improve the applicability of this virus as a biological control agent.

Mixtures of Insect-Pathogenic Viruses in a Single Virion: towards the Development of Custom-Designed Insecticides

Alphabaculoviruses (Baculoviridae) are pathogenic DNA viruses of Lepidoptera that have applications as the basis for biological insecticides and expression vectors in biotechnological processes. These viruses have a characteristic physical structure that facilitates the transmission of groups of genomes. We demonstrate that coinfection of a susceptible insect by two different alphabaculovirus species results in the production of mixed-virus occlusion bodies containing the parental viruses. This occurred between closely related and phylogenetically more distant alphabaculoviruses. Approximately half the virions present in proteinaceous viral occlusion bodies produced following coinfection of insects with a mixture of two alphabaculoviruses contained both viruses, indicating that the viruses coinfected and replicated in a single cell and were coenveloped within the same virion. This observation was confirmed by endpoint dilution assay. Moreover, both viruses persisted in the mixed-virus population by coinfection of insects during several rounds of insect-to-insect transmission. Coinfection by viruses that differed in genome size had unexpected results on the length of viral nucleocapsids, which differed from those of both parental viruses. These results have unique implications for the development of alphabaculoviruses as biological control agents of insect pests.IMPORTANCE Alphabaculoviruses are used as biological insecticides and expression vectors in biotechnology and medical applications. We demonstrate that in caterpillars infected with particular mixtures of viruses, the genomes of different baculovirus species can be enveloped together within individual virions and occluded within proteinaceous occlusion bodies. This results in the transmission of mixed-virus populations to the caterpillar stages of moth species. Once established, mixed-virus populations persist by coinfection of insect cells during several rounds of insect-to-insect transmission. Mixed-virus production technology opens the way to the development of custom-designed insecticides for control of different combinations of caterpillar pest species.

Population structure of Cydia pomonella granulovirus isolates revealed by quantitative analysis of genetic variation

Genetic diversity of viruses is driven by genomic mutations and selection through its host, resulting in differences in virulence as well as host responses. For baculoviruses, which are naturally occurring pathogens of insects and which are frequently sprayed on hundred thousands to millions of hectares as biocontrol agents of insect pests, the phenomenon of virus-host co-evolution is of particular scientific interest and economic importance because high virulence of baculovirus products is essential and emergence of host resistance needs to be avoided as much as possible. In the present study, the population structure of twenty isolates of the Cydia pomonella granulovirus (CpGV), including twelve isolates from different geographic origins and eight commercial formulations, were studied on the basis of next-generation sequencing data and by analyzing the distribution of single nucleotide polymorphisms (SNPs). An entirely consensus sequence-free quantitative SNP analysis was applied for the identification of 753 variant SNP sites being specific for single as well as groups of CpGV isolates. Based on the quantitative SNP analysis, homogenous, heterogenous as well as mixed isolates were identified and their proportions of genotypes were deciphered, revealing a high genetic diversity of CpGV isolates from around the world. Based on hierarchical clustering on principal components (HCPC), six distinct isolate/group clusters were identified, representing the proposed main phylogenetic lineages of CpGV but comprising full genome information from virus mixtures. The relative location of different isolates in HCPC reflected the proportion of variable compositions of different genotypes. The established methods provide novel analysis tools to decipher the molecular complexity of genotype mixtures in baculovirus isolates, thus depicting the population structure of baculovirus isolates in a more adequate form than consensus based analyses.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

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.

Bacsnp: Using Single Nucleotide Polymorphism (SNP) Specificities and Frequencies to Identify Genotype Composition in Baculoviruses

Natural isolates of baculoviruses (as well as other dsDNA viruses) generally consist of homogenous or heterogenous populations of genotypes. The number and positions of single nucleotide polymorphisms (SNPs) from sequencing data are often used as suitable markers to study their genotypic composition. Identifying and assigning the specificities and frequencies of SNPs from high-throughput genome sequencing data can be very challenging, especially when comparing between several sequenced isolates or samples. In this study, the new tool “bacsnp”, written in R programming langue, was developed as a downstream process, enabling the detection of SNP specificities across several virus isolates. The basis of this analysis is the use of a common, closely related reference to which the sequencing reads of an isolate are mapped. Thereby, the specificities of SNPs are linked and their frequencies can be used to analyze the genetic composition across the sequenced isolate. Here, the downstream process and analysis of detected SNP positions is demonstrated on the example of three baculovirus isolates showing the fast and reliable detection of a mixed sequenced sample.

Using Direct RNA Nanopore Sequencing to Deconvolute Viral Transcriptomes

The genomes of DNA viruses encode deceptively complex transcriptomes evolved to maximize coding potential within the confines of a relatively small genome. Defining the full range of viral RNAs produced during an infection is key to understanding the viral replication cycle and its interactions with the host cell. Traditional short-read (Illumina) sequencing approaches are problematic in this setting due to the difficulty of assigning short reads to individual RNAs in regions of transcript overlap and to the biases introduced by the required recoding and amplification steps. Additionally, different methodologies may be required to analyze the 5' and 3' ends of RNAs, which increases both cost and effort. The advent of long-read nanopore sequencing simplifies this approach by providing a single assay that captures and sequences full length RNAs, either in cDNA or native RNA form. The latter is particularly appealing as it reduces known recoding biases whilst allowing more advanced analyses such as estimation of poly(A) tail length and the detection of RNA modifications including N6 -methyladenosine. Using herpes simplex virus (HSV)-infected primary fibroblasts as a template, we provide a step-by-step guide to the production of direct RNA sequencing libraries suitable for sequencing using Oxford Nanopore Technologies platforms and provide a simple computational approach to deriving a high-quality annotation of the HSV transcriptome from the resulting sequencing data. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Productive infection of primary fibroblasts with herpes simplex virus Support Protocol: Cell passage and plating of primary fibroblasts Basic Protocol 2: Preparation and sequencing of dRNA-seq libraries from virus-infected cells Basic Protocol 3: Processing, alignment, and analysis of dRNA-seq datasets.