Baculovirus Molecular Evolution via Gene Turnover and Recurrent Positive Selection of Key Genes

Sequencing, Analysis and Organization of the Complete Genome of a Novel Baculovirus Calliteara abietis Nucleopolyhedrovirus (CaabNPV)

Baculoviridae, a virus family characterized by a single large double stranded DNA, encompasses the majority of viral bioinsecticides, representing a highly promising and environmentally friendly pesticide approach to insect control. This study focuses on the characterization of a baculovirus isolated from larvae of Calliteara abietis (Erebidae, Lymantriidae) collected in Mongolian pinaceae forests. This new isolate was called Calliteara abietis nucleopolyhedrovirus (CaabNPV). CaabNPV exhibits an irregular polyhedron shape, and significant variation in the diameter of its occlusion bodies (OBs) was observed. Nucleotide distance calculations confirmed CaabNPV as a novel baculovirus. The CaabNPV genome spans 177,161 bp with a G+C content of 45.12% and harbors 150 potential open reading frames (ORFs), including 38 core genes. A comprehensive genomic analysis categorizes CaabNPV within Group II alphabaculovirus, revealing a close phylogenetic relationship with Alphabaculovirus orleucostigmae (OrleNPV). Additionally, repeat sequence analysis identified three highly repetitive sequences consisting of 112 bp repeat units, known as homologous regions (hrs). This research contributes valuable insights into CaabNPV's phylogenetic placement, genomic structure, and its potential applications in insect biocontrol.

Investigation on Pathological Aspects, Mode of Transmission, and Tissue Tropism of Antheraea proylei Nucleopolyhedrovirus Infecting Oak Tasar Silkworm

The temperate oak tasar silkworm, Antheraea proylei, is frequently infested with Antheraea proylei nucleopolyhedrovirus (AnprNPV) causing tiger band disease. This disease is one of the key factors that obstructs production and productivity of oak tasar sericulture. The current study aimed to investigate the pathogenicity of AnprNPV, its mode of transmission, and detection of AnprNPV in different tissues. Transmission electron micrographs of AnprNPV showed single rod-shaped bodies and occlusion derived virus (ODV) enclosed within multiple envelopes. The infecting AnprNPV displayed tissue tropism with higher copy numbers detected in the insect fat body and ovary. The virus was observed to multiply in all developmental stages of the silkworm such as egg, larva, pupa, and moth, confirming its ability to spread throughout the silkworm lifecycle. Baculovirus isolated from infected A. proylei showed cross-infectivity in other Saturniidae wild silkworm species such as Antheraea pernyi, A. frithi, and Samia ricini, widening their probable host range for infection. Baculoviruses generally display a horizontal mode of transmission, mainly through ingestion of occlusion bodies (OBs); however, the present study revealed a trans-ovum vertical mode of transmission in addition to a horizontal mode. The observations made in this study aid a detailed understanding of the tiger band disease and its causative pathogen AnprNPV, which will support future studies and disease management in oak tasar sericulture.

Genomic analysis of two Chinese isolates of hyphantria cunea nucleopolyhedrovirus reveals a novel species of alphabaculovirus that infects hyphantria cunea drury (lepidoptera: arctiidae)

Background

Baculoviruses act as effective biological control agents against the invasive pest Hyphantria cunea Drury. In this study, two Chinese Hyphantria cunea nucleopolyhedrovirus (HycuNPV) isolates, HycuNPV-BJ and HycuNPV-HB, were deep sequenced and compared with the Japanese isolate, HycuNPV-N9, to determine whole-genome level diversity and evolutionary history.

Results

The divergence of the phylogenetic tree and the K2P distances based on 38 core-gene concatenated alignment revealed that two Chinese HycuNPV isolates were a novel species of Alphabaculovirus that infected Hyphantria cunea in China. The gene contents indicated significant differences in the HycuNPV genomes between the Chinese and Japanese isolates. The differences included gene deletions, acquisitions and structural transversions, but the main difference was the high number of single nucleotide polymorphisms (SNPs). In total, 10,393 SNPs, corresponding to approximately 8% of the entire HycuNPV-N9 genome sequence, were detected in the aligned reads. By analyzing non-synonymous variants, we found that hotspot mutation-containing genes had mainly unknown functions and most were early expressing genes. We found that the hycu78 gene which had early and late promoter was under positive selection. Biological activity assays revealed that the infectivity of HycuNPV-HB was greater than that of HycuNPV-BJ, and the killing speed of HycuNPV-HB was faster than that of HycuNPV-BJ. A comparison of molecular genetic characteristics indicated that the virulence differences between the two isolates were affected by SNP and structural variants, especially the homologous repeat regions.

Conclusions

The genomes of the two Chinese HycuNPV isolates were characterized, they belonged to a novel species of Alphabaculovirus that infected Hyphantria cunea in China. We inferred that the loss or gain of genetic material in the HycuNPV-HB and HycuNPV-BJ genomes resulted in new important adaptive capabilities to the H. cunea host. These results extend the current understanding of the genetic diversity of HycuNPV and will be useful for improving the applicability of this virus as a biological control agent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08604-7.

Role of Diversity and Recombination in the Emergence of Chilli Leaf Curl Virus

Chilli leaf curl virus (ChiLCV), (Genus Begomovirus, family Geminiviridae) and associated satellites pose a serious threat to chilli production, worldwide. This study highlights the factors accountable for genetic diversity, recombination, and evolution of ChiLCV, and associated chilli leaf curl alphasatellite (ChiLCA) and chilli leaf curl betasatellite (ChiLCB). Phylogenetic analysis of complete genome (DNA-A) sequences of 132 ChiLCV isolates from five countries downloaded from NCBI database clustered into three major clades and showed high population diversity. The dN/dS ratio and Tajima D value of all viral DNA-A and associated betasatellite showed selective control on evolutionary relationships. Negative values of neutrality tests indicated purified selection and an excess of low-frequency polymorphism. Nucleotide diversity (π) for C4 and Rep genes was higher than other genes of ChiLCV with an average value of π = 18.37 × 10⁻² and π = 17.52 × 10⁻² respectively. A high number of mutations were detected in TrAP and Rep genes, while ChiLCB has a greater number of mutations than ChiLCA. In addition, significant recombination breakpoints were detected in all regions of ChiLCV genome, ChiLCB and, ChiLCA. Our findings indicate that ChiLCV has the potential for rapid evolution and adaptation to a range of geographic conditions and could be adopted to infect a wide range of crops, including diverse chilli cultivars.

Differential insecticidal properties of Spodoptera frugiperda multiple nucleopolyhedrovirus isolates against corn-strain and rice-strain fall armyworm, and genomic analysis of three isolates

The fall armyworm, Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae) is a destructive crop pest native to North, Central, and South America that recently has spread to Africa and Asia. Isolates of Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) have the potential to be developed as low-risk biopesticides for management of fall armyworm, and a commercially available formulation has been developed for control of fall armyworm in North and South America. In this study, the virulence (LC₅₀ and LT₅₀) of several SfMNPV isolates towards larvae of both corn-strain and rice-strain fall armyworm was assessed. Bioassays with corn-strain larvae revealed that the isolates could be organized into fast-killing (LT₅₀ < 56 h post-infection) and slow-killing (LT₅₀ > 68 h post-infection) groups. Rice-strain larvae exhibited narrower ranges of susceptibility to baculovirus infection and of survival times in bioassays with different isolates. Two SfMNPV isolates with rapid speeds of kill (SfMNPV-459 from Colombia and SfMNPV-1197 from Georgia, USA) along with an isolate that killed corn-strain at relatively low concentrations (SfMNPV-281 from Georgia) were selected for the complete determination of their genome sequences. The SfMNPV-1197 genome sequence shared high sequence identity with genomes of a Nicaraguan isolate, while SfMNPV-281 formed a separate clade with a USA and a Brazilian isolate in phylogenetic trees. The SfMNPV-459 sequence was more divergent with the lowest genome sequence identities in pairwise alignments with other sequenced SfMNPV genomes, and was not grouped reliably with either the 1197 clade or the 281 clade. SfMNPV-459 contained homologs of two ORFs that were unique to another Colombian isolate, but these isolates were not placed in the same clade in phylogenetic trees. This study identifies isolates with superior properties for control of fall armyworm and adds to our knowledge of the genetics of SfMNPV.

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.

Mosaic genome evolution and phylogenetics of Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) and virulence of seven new isolates from the Brazilian states of Minas Gerais and Mato Grosso

In a comparative analysis of genome sequences from isolates of the baculovirus Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) from Brazil and Guatemala, we identified a subset of isolates possessing chimeric genomes. We identified six distinct phylogenetically incongruous regions (PIRs) dispersed in the genomes, of between 279 and 3345 bp in length. The individual PIRs possessed high sequence similarity among the affected ChinNPV isolates but varied in coverage in some instances. The donor for four of the PIRs implicated in horizontal gene transfer (HGT) was identified as Trichoplusia ni single nucleopolyhedrovirus (TnSNPV), an alphabaculovirus closely related to ChinNPV, or another unknown but closely related virus. BLAST searches of the other two PIRs returned only ChinNPV sequences, but HGT from an unknown donor baculovirus cannot be excluded. Although Chrysodeixis includens and Trichoplusia ni are frequently co-collected from soybean fields in Brazil, pathogenicity data suggest that natural coinfection of C. includens larvae with ChinNPV and TnSNPV is probably uncommon. Additionally, since the chimeric ChinNPV genomes with tracts of TnSNPV sequence were restricted to a single monophyletic lineage of closely related isolates, a model of progressive restoration of the native DNA sequence by recombination with ChinNPV possessing a fully or partially non-chimeric genome is reasonable. However, multiple independent HGT from TnSNPV to ChinNPV during the evolution of these isolates cannot be excluded. Mortality data suggest that the ChinNPV isolates with chimeric genomes are not significantly different in pathogenicity towards C. includens when compared to most other ChinNPV isolates. Exclusion of the PIRs prior to phylogenetic analysis had a large impact on the topology of part of the maximum-likelihood tree, revealing a homogenous clade of three isolates (IB, IC and ID) from Paraná state in Brazil collected in 2006, together with an isolate from Guatemala collected in 1972 (IA), comprising the lineage uniquely affected by HGT from TnSNPV. The other 10 Brazilian ChinNPV isolates from Paraná, Mato Grosso, and Minas Gerais states showed higher variability, where only three isolates from Paraná state formed a monophyletic group correlating with geographical origin.

Recurrent evolution of high virulence in isolated populations of a DNA virus

Hosts and viruses are constantly evolving in response to each other: as a host attempts to suppress a virus, the virus attempts to evade and suppress the host's immune system. Here, we describe the recurrent evolution of a virulent strain of a DNA virus, which infects multiple Drosophila species. Specifically, we identified two distinct viral types that differ 100-fold in viral titer in infected individuals, with similar differences observed in multiple species. Our analysis suggests that one of the viral types recurrently evolved at least four times in the past ~30,000 years, three times in Arizona and once in another geographically distinct species. This recurrent evolution may be facilitated by an effective mutation rate which increases as each prior mutation increases viral titer and effective population size. The higher titer viral type suppresses the host-immune system and an increased virulence compared to the low viral titer type.

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.

Biology and Interaction of the Natural Occurrence of Distinct Monopartite Begomoviruses Associated With Satellites in Capsicum annum From India

Chili (Capsicum annuum L.) is an important vegetable and spice crop of tropical and sub-tropical regions. Chili plants showing upward leaf curling, leaf crinkling, and leaf yellowing symptoms, collected from Sikar district of Rajasthan, India, were found to be associated with begomovirus and satellite molecules. The presence of virus was confirmed by PCR using virus-specific primer. The full-length genomic DNA-A of three begomovirus (MM-1, CS-1 and RV-1) and two satellites (MM-2 and MM-3) were cloned which was identified from single symptomatic chili plant. The genome organization of isolated three viruses is similar to those of other Old World monopartite begomoviruses. The comparison of the sequences and closest phylogenetic relationships for the begomoviruses, betasatellite and alphasatellite DNAs revealed that MM-1 was designated as DNA-A of Chili leaf curl virus (ChiLCV), CS-1 is considered to be a new distinct species of Tomato leaf curl Gujrat virus (ToLCGV) whereas RV-1 as a new strain of Cotton leaf curl Multan virus (CLCuMuV). The DNA-A component of ChiLCV showed 8.6%, ToLCGV of 16.6% and CLCuMuV of 7.7% average evolutionary divergence, concomitantly, the betasatellite and alphasatellite molecule had 9.9% and 5.9% overall sequence divergence, respectively. Interestingly, most of the begomoviruses were found to be intra-species recombinants. The dN/dS ratio and Tajima D value of all viral DNA-A component and their associated betasatellite showed their selective control on evolutionary relationships. The nucleotide substitution rates were determined for the DNA-A genomes of ChiLCV (7.22 × 10^–4 substitutions site^–1 year^–1), CLCuMuV (1.49 × 10^–4 substitutions site^–1 year^–1), ToLCGV (7.47 × 10^–4 substitutions site^–1 year^–1), the genome of associated ChiLCB (4.20 × 10^–4 substitutions site^–1 year^–1) and CLCuMuA (1.49 × 10^–4 substitutions site^–1 year^–1). Agro-inoculation studies indicate that the presence of DNA betasatellite induce severe symptoms in N. benthamiana and chili, suggesting prerequisite association for typical disease development.

Identification of Loci Associated with Enhanced Virulence in Spodoptera litura Nucleopolyhedrovirus Isolates Using Deep Sequencing

Spodoptera litura is an emerging pest insect in cotton and arable crops in Central Asia. To explore the possibility of using baculoviruses as biological control agents instead of chemical pesticides, in a previous study we characterized a number of S. litura nucleopolyhedrovirus (SpltNPV) isolates from Pakistan. We found significant differences in speed of kill, an important property of a biological control agent. Here we set out to understand the genetic basis of these differences in speed of kill, by comparing the genome of the fast-killing SpltNPV-Pak-TAX1 isolate with that of the slow-killing SpltNPV-Pak-BNG isolate. These two isolates and the SpltNPV-G2 reference strain from China were deep sequenced with Illumina. As expected, the two Pakistani isolates were closely related with >99% sequence identity, whereas the Chinese isolate was more distantly related. We identified two loci that may be associated with the fast action of the SpltNPV-Pak-TAX1 isolate. First, an analysis of rates of synonymous and non-synonymous mutations identified neutral to positive selection on open reading frame (ORF) 122, encoding a viral fibroblast growth factor (vFGF) that is known to affect virulence in other baculoviruses. Second, the homologous repeat region hr17, a putative enhancer of transcription and origin of replication, is absent in SpltNPV-Pak-TAX1 suggesting it may also affect virulence. Additionally, we found there is little genetic variation within both Pakistani isolates, and we identified four genes under positive selection in both isolates that may have played a role in adaptation of SpltNPV to conditions in Central Asia. Our results contribute to the understanding of the enhanced activity of SpltNPV-Pak-TAX1, and may help to select better SpltNPV isolates for the control of S. litura in Pakistan and elsewhere.

Biological and molecular characterization of two Anticarsia gemmatalis multiple nucleopolyhedrovirus clones exhibiting contrasting virulence

Baculovirus natural populations are known to be genetically heterogeneous and such genotypic diversity could have implications in the performance of biocontrol agents. The Anticarsia gemmatalis nucleopolyhedrovirus (AgMNPV) has been widely used to control the velvetbean caterpillar, Anticarsia gemmatalis, in Brazil. In the present work, morphological and molecular analyses as well as the biological activity of AgMNPV genotypes derived from a Brazilian field isolate (AgMNPV-79) were carried out. The existence of genotypic variants in the population was confirmed by DNA restriction analysis. Although difference in virulence was observed among the variants, the most (Ag79-01) and the least (AgL-16) virulent clones do not show any morphological and cytopathological changes when compared to the most studied isolate (AgMNPV-2D). The complete genome analysis of the two viral clones showed the presence of single open reading frames (ORFs) of the pe-38 and he65 genes, which contrasts with the two split ORFs present in the genome of the AgMNPV-2D isolate. The viral clone AgL-16 has many variations in the ie-2 and pe-38 genes, which are transcription regulatory genes responsible for the regulation of viral early gene expression during insect cell infection. Furthermore, other genes showed alterations like the odv-e56, which have an essential role in the maturation and envelopment of the ODVs, and bro-a and bro-b genes which were fused to form a single ORF. For the Ag79-01, although the total number of single nucleotide variants (SNVs) was more prominent in the pe-38 gene, its genome showed very few modifications in comparison to the AgMNPV-2D genome.

Exploitation of receptor tyrosine kinases by viral-encoded growth factors

Receptor tyrosine kinases (RTKs) are essential components of cell communication pathways utilized from the embryonic to adult stages of life. These transmembrane receptors bind polypeptide ligands, such as growth factors, inducing signalling cascades that control cellular processes such as proliferation, survival, differentiation, motility and inflammation. Many viruses have acquired homologs of growth factors encoded by the hosts that they infect. Production of growth factors during infection allows viruses to exploit RTKs for entry and replication in cells, as well as for host and environmental dissemination. This review describes the genetic diversity amongst virus-derived growth factors and the mechanisms by which RTK exploitation enhances virus survival, then highlights how viral ligands can be used to further understanding of RTK signalling and function during embryogenesis, homeostasis and disease scenarios.

The dynamic evolution of Drosophila innubila Nudivirus

Viruses coevolve with their hosts to overcome host resistance and gain the upper hand in the evolutionary arms race. Drosophila innubila nudivirus (DiNV) is a double stranded DNA virus, closely related to Oryctes rhinoceros nudivirus (OrNV) and Kallithea virus. DiNV is the first DNA virus found to naturally infect Drosophila and therefore has the potential to be developed as a model for DNA virus immune defense and host/virus coevolution within its well-studied host system. Here we sequence and annotate the genome of DiNV and identify signatures of adaptation, revealing clues for genes involved in host-parasite coevolution. The circular genome is 155,555bp and contains 107 coding open reading frames (ORFs) and a wealth of AT-rich simple sequence repeats. While synteny is highly conserved between DiNV and Kallithea virus, it drops off rapidly as sequences become more divergent, consistent with rampant rearrangements across nudiviruses. Overall, we show that evolution of DiNV is likely due to adaptation of a very few genes coupled with high gene turnover.