Identification and functional analysis of the origins of DNA replication in the Cydia pomonella granulovirus genome

A Renewed Appreciation of Helicoverpa armigera Nucleopolyhedrovirus BJ (Formerly Helicoverpa assulta Nucleopolyhedrovirus) with Whole Genome Sequencing

Helicoverpa assulta is a pest that causes severe damage to tobacco, pepper and other cash crops. A local strain of HearNPV-BJ (formerly Helicoverpa assulta nucleopolyhedrovirus (HeasNPV-DJ0031)) was isolated from infected H. assulta larvae in Beijing, which had been regarded as a new kind of baculovirus in previous studies. Describing the biological characteristics of the strain, including its external morphology, internal structure and the pathological characteristics of the infection of various cell lines, can provide references for the identification and function of the virus. HearNPV-BJ virion was defined as a single-nucleocapsid nucleopolyhedrovirus by scanning electron microscopy. QB-Ha-E-5 (H. armigera) and BCIRL-Hz-AM1 (H. zea) cell lines were sensitive to HearNPV-BJ. Undoubtedly modern developed sequencing technology further facilitates the increasing understanding of various strains. The whole genome sequence of the HearNPV-BJ was sequenced and analyzed. The HearNPV-BJ isolate genome was 129, 800 bp nucleotides in length with a G + C content of 38.87% and contained 128 open reading frames (ORFs) encoding predicted proteins of 50 or over 50 amino acids, 67 ORFs in the forward orientation and 61 ORFs in the reverse orientation, respectively. The genome shared 99% sequence identity with Helicoverpa armigera nucleopolyhedrovirus C1 strain (HearNPV-C1), and 103 ORFs had very high homology with published HearNPV sequences. Two bro genes and three hrs were found to be dispersed along the HearNPV-BJ genome. Three of the highest homologs, ORFs with HearNPV, were smaller due to the earlier appearance of the stop codon with unknown functions. P6.9 of HearNPV-BJ, a structural protein, is distinctly different from that of Autographa californica nucleopolyhedrovirus (AcMNPV); its homology with the corresponding gene in HearNPV-C1 was 93.58%. HearNPV-BJ contains 38 core genes identified in other baculoviruses, and phylogenetic analysis indicates HearNPV-BJ belongs to Alphabaculovirus Group II, same as HearNPV-C1. The resulting data provide a better understanding of virion structure, gene function and character of infection. By supplementing the whole-genome sequencing data and Kimura-2 model index, there is more evidence to indicate that HearNPV-BJ may be a variant of Helicoverpa armigera nucleopolyhedrovirus, which also deepens our understanding of the virus species demarcation criteria.

Genome sequence and organization of the Mythimna (formerly Pseudaletia) unipuncta granulovirus Hawaiian strain

Purified occlusion bodies (OBs) of Mythimna (formerly Pseudaletia) unipuncta (the true armyworm) granulovirus Hawaiian strain (MyunGV-A) were observed, showing typical GV morphological characteristics under scanning and transmission electron microscopy (EM). The genome of MyunGV-A was completely sequenced and analysed. The genome is 176,677 bp in size, with a G+C content of 39.79%. It contains 183 open reading frames (ORFs) encoding 50 or more amino acids with minimal overlap. Comparison of MyunGV-A with TnGV, XcGV, and HearGV genomes revealed extensive sequence similarity and collinearity, and the four genomes contain the same nine homologous regions (hrs) with conserved structures and locations. Three unique genes, 12 baculovirus repeated ORF (bro), 2 helicase, and 3 enhancin genes, were identified. In particular, two repeated genes (ORF39 and 49) are present in the genome, in reverse and complementarily orientations. Twenty-four OB proteins were identified from the putative protein database of MyunGV-A. In addition, MyunGV-A belongs to the Betabaculovirus group and is most closely related to TnGV (99% amino acid identity) according to a phylogenetic tree based on the combined amino acid sequences of 38 core gene contents.

Genome sequence analysis of a Helicoverpa armigera single nucleopolyhedrovirus (HearNPV-TR) isolated from Heliothis peltigera in Turkey

The entire genome of Helicoverpa armigera single nucleopolyhedrovirus (HearNPV-TR) was sequenced, and compared to genomes of other existing isolates. HearNPV-TR genome is 130.691 base pairs with a 38.9% G+C content and has 137 open reading frames (ORFs) of ≥ 150 nucleotides. Five homologous repeated sequences (hrs) and two baculovirus repeated ORFs (bro-a and bro-b) were identified. Phylogenetic analysis showed that HearNPV-TR is closer to HaSNPV-C1, HaSNPV-G4, HaSNPV-AU and HasNPV. However, there are significant differences in hr3, hr5 regions and in bro-a gene. Pairwise Kimura-2 parameter analysis of 38 core genes sequences of HearNPV-TR and other Helicoverpa NPVs showed that the genetic distances for these sequences were below 0.015 substitutions/site. Genomic differences as revealed by restriction profiles indicated that hr3, hr5 regions and bro-a gene may play a role in the virulence of HearNPV-TR.

Genomic sequencing of Troides aeacus nucleopolyhedrovirus (TraeNPV) from golden birdwing larvae (Troides aeacus formosanus) to reveal defective Autographa californica NPV genomic features

Background

The golden birdwing butterfly (Troides aeacus formosanus) is a rarely observed species in Taiwan. Recently, a typical symptom of nuclear polyhedrosis was found in reared T. aeacus larvae. From the previous Kimura-2 parameter (K-2-P) analysis based on the nucleotide sequence of three genes in this isolate, polh, lef-8 and lef-9, the underlying virus did not belong to any known nucleopolyhedrovirus (NPV) species. Therefore, this NPV was provisionally named “TraeNPV”. To understand this NPV, the nucleotide sequence of the whole TraeNPV genome was determined using next-generation sequencing (NGS) technology.

Results

The genome of TraeNPV is 125,477 bp in length with 144 putative open reading frames (ORFs) and its GC content is 40.45%. A phylogenetic analysis based on the 37 baculoviral core genes suggested that TraeNPV is a Group I NPV that is closely related to Autographa californica nucleopolyhedrovirus (AcMNPV). A genome-wide analysis showed that TraeNPV has some different features in its genome compared with other NPVs. Two novel ORFs (Ta75 and Ta139), three truncated ORFs (pcna, he65 and bro) and one duplicated ORF (38.7 K) were found in the TraeNPV genome; moreover, there are fewer homologous regions (hrs) than there are in AcMNPV, which shares eight hrs within the TraeNPV genome. TraeNPV shares similar genomic features with AcMNPV, including the gene content, gene arrangement and gene/genome identity, but TraeNPV lacks 15 homologous ORFs from AcMNPV in its genome, such as ctx, host cell-specific factor 1 (hcf-1), PNK/PNL, vp15, and apsup, which are involved in the auxiliary functions of alphabaculoviruses.

Conclusions

Based on these data, TraeNPV would be clarified as a new NPV species with defective AcMNPV genomic features. The precise relationship between TraeNPV and other closely related NPV species were further investigated. This report could provide comprehensive information on TraeNPV for evolutionary insights into butterfly-infected NPV.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5713-2) contains supplementary material, which is available to authorized users.

Genome sequence analysis and organization of the Hyphantria cunea granulovirus (HycuGV-Hc1) from Turkey

The fall webworm (Hyphantria cunea) impacts a wide variety of crops and cultivated broadleaf plant species. The pest is native to North America, was introduced to Europe and has since spread further as far as central Asia. Despite several attempts to control its distribution, the pest continues to spread causing damage all over the world. A naturally occurring baculovirus, Hyphantria cunea granulovirus (HycuGV-Hc1), isolated from the larvae of H. cunea in Turkey appears to have a potential as microbial control agent against this pest. In this report we describe the complete genome sequence and organization of the granulovirus isolate (HycuGV-Hc1) that infects the larval stages and compare it to other baculovirus genomes. The HycuGV-Hc1 genome is a circular double-stranded DNA of 114,825 bp in size with a nucleotide distribution of 39.3% G + C. Bioinformatics analysis predicted 132 putative open reading frames of (ORFs) ≥ 150 nucleotides. There are 24 ORFs with unknown function. Seven homologous repeated regions (hrs) and two bro genes (bro-1 and bro-2) were identified in the genome. Comparison to other baculovirus genomes, HycuGV-Hc1 revealed some differences in gene content and organization. Gene parity plots and phylogenetics confirmed that HycuGV-Hc1 is a Betabaculovirus and is closely related to Plutella xylostella granulovirus. This study expands our knowledge on the genetic variation of HycuGV isolates and provides further novel knowledge on the nature of granuloviruses.

Complete genome sequence analysis of the Malacosoma neustria nucleopolyhedrovirus from Turkey

The lackey moth, Malacosoma neustria (Linnaeus, 1758), a worldwide pest, causes extensive economic losses particularly on hazelnut, plum, oak, poplar, and willow trees. A baculovirus, Malacosoma neustria nucleopolyhedrovirus (ManeNPV-T2), has been isolated from the larvae collected in Turkey and appears to have a potential as a microbial control agent. In this study, we describe the complete genome sequence of ManeNPV-T2 and compare it to other sequenced baculovirus genomes. The ManeNPV-T2 genome is a circular double-stranded DNA molecule of 130,202 bp, has 38.2% G + C, and is predicted to contain 131 putative open reading frames (ORFs) each with a coding capacity of more then 50 amino acids. There are 27 ORFs with unknown function of which 6 are unique to ManeNPV-T2. Eleven homologous regions (hrs) and two bro genes (bro-a and bro-b) were identified in the genome. There are two homologues of chaB and nicotinamide riboside kinase-1 genes, separated from themselves with a few nucleotides. Additionally, ac145, thought to be per os infectivity factor (pif) gene, is also found as two homologues. All 38 core genes are found in the ManeNPV-T2 genome. The phylogenetic tree of ManeNPV-T2 in relation to 50 other baculoviruses whose genomes have been completely sequenced showed ManeNPV-T2 to be closely related to the group II NPVs. This study expands our knowledge on baculoviruses, describes the characterization ManeNPV, and ultimately contributes to the registration of this virus as a microbial pesticide.

The Complete Genome Sequence of a Second Distinct Betabaculovirus from the True Armyworm, Mythimna unipuncta

The betabaculovirus originally called Pseudaletia (Mythimna) sp. granulovirus #8 (MyspGV#8) was examined by electron microscopy, host barcoding PCR, and determination of the nucleotide sequence of its genome. Scanning and transmission electron microscopy revealed that the occlusion bodies of MyspGV#8 possessed the characteristic size range and morphology of betabaculovirus granules. Barcoding PCR using cytochrome oxidase I primers with DNA from the MyspGV#8 collection sample confirmed that it had been isolated from the true armyworm, Mythimna unipuncta (Lepidoptera: Noctuidae) and therefore was renamed MyunGV#8. The MyunGV#8 genome was found to be 144,673 bp in size with a nucleotide distribution of 49.9% G+C, which was significantly smaller and more GC-rich than the genome of Pseudaletia unipuncta granulovirus H (PsunGV-H), another M. unipuncta betabaculovirus. A phylogeny based on concatenated baculovirus core gene amino acid sequence alignments placed MyunGV#8 in clade a of genus Betabaculovirus. Kimura-2-parameter nucleotide distances suggested that MyunGV#8 represents a virus species different and distinct from other species of Betabaculovirus. Among the 153 ORFs annotated in the MyunGV#8 genome, four ORFs appeared to have been obtained from or donated to the alphabaculovirus lineage represented by Leucania separata nucleopolyhedrovirus AH1 (LeseNPV-AH1) during co-infection of Mythimna sp. larvae. A set of 33 ORFs was identified that appears only in other clade a betabaculovirus isolates. This clade a-specific set includes an ORF that encodes a polypeptide sequence containing a CIDE_N domain, which is found in caspase-activated DNAse/DNA fragmentation factor (CAD/DFF) proteins. CAD/DFF proteins are involved in digesting DNA during apoptosis.

The Complete Genome Sequence of Plodia Interpunctella Granulovirus: Evidence for Horizontal Gene Transfer and Discovery of an Unusual Inhibitor-of-Apoptosis Gene

The Indianmeal moth, Plodia interpunctella (Lepidoptera: Pyralidae), is a common pest of stored goods with a worldwide distribution. The complete genome sequence for a larval pathogen of this moth, the baculovirus Plodia interpunctella granulovirus (PiGV), was determined by next-generation sequencing. The PiGV genome was found to be 112, 536 bp in length with a 44.2% G+C nucleotide distribution. A total of 123 open reading frames (ORFs) and seven homologous regions (hrs) were identified and annotated. Phylogenetic inference using concatenated alignments of 36 baculovirus core genes placed PiGV in the “b” clade of viruses from genus Betabaculovirus with a branch length suggesting that PiGV represents a distinct betabaculovirus species. In addition to the baculovirus core genes and orthologues of other genes found in other betabaculovirus genomes, the PiGV genome sequence contained orthologues of the bidensovirus NS3 gene, as well as ORFs that occur in alphabaculoviruses but not betabaculoviruses. While PiGV contained an orthologue of inhibitor of apoptosis-5 (iap-5), an orthologue of inhibitor of apoptosis-3 (iap-3) was not present. Instead, the PiGV sequence contained an ORF (PiGV ORF81) encoding an IAP homologue with sequence similarity to insect cellular IAPs, but not to viral IAPs. Phylogenetic analysis of baculovirus and insect IAP amino acid sequences suggested that the baculovirus IAP-3 genes and the PiGV ORF81 IAP homologue represent different lineages arising from more than one acquisition event. The presence of genes from other sources in the PiGV genome highlights the extent to which baculovirus gene content is shaped by horizontal gene transfer.

Genome of Cnaphalocrocis medinalis Granulovirus, the First Crambidae-Infecting Betabaculovirus Isolated from Rice Leaffolder to Sequenced

Cnaphalocrocis medinalis is a major pest of rice in South and South-East Asia. Insecticides are the major means farmers use for management. A naturally occurring baculovirus, C. medinalis granulovirus (CnmeGV), has been isolated from the larvae and this has the potential for use as microbial agent. Here, we described the complete genome sequence of CnmeGV and compared it to other baculovirus genomes. The genome of CnmeGV is 112,060 base pairs in length, has a G+C content of 35.2%. It contains 133 putative open reading frames (ORFs) of at least 150 nucleotides. A hundred and one (101) of these ORFs are homologous to other baculovirus genes including 37 baculovirus core genes. Thirty-two (32) ORFs are unique to CnmeGV with no homologues detected in the GeneBank and 53 tandem repeats (TRs) with sequence length from 25 to 551 nt intersperse throughout the genome of CnmeGV. Six (6) homologous regions (hrs) were identified interspersed throughout the genome. Hr2 contains 11 imperfect palindromes and a high content of AT sequence (about 73%). The unique ORF28 contains a coiled-coil region and a zinc finger-like domain of 4–50 residues specialized by two C2C2 zinc finger motifs that putatively bound two atoms of zinc. ORF21 encoding a chit-1 protein suggesting a horizontal gene transfer from alphabaculovirus. The putative protein presents two carbohydrate-binding module family 14 (CBM_14) domains rather than other homologues detected from betabaculovirus that only contains one chit-binding region. Gene synteny maps showed the colinearity of sequenced betabaculovirus. Phylogenetic analysis indicated that CnmeGV grouped in the betabaculovirus, with a close relation to AdorGV. The cladogram obtained in this work grouped the 17 complete GV genomes in one monophyletic clade. CnmeGV represents a new crambidae host-isolated virus species from the genus Betabaculovirus and is most closely relative of AdorGV. The analyses and information derived from this study will provide a better understanding of the pathological symptoms caused by this virus and its potential use as a microbial pesticide.

Possible functional co-operation of palindromes hr3 and hr4 in the genome of Cydia pomonella granulovirus affects viral replication capacity

After previous studies had shown that natural transposon insertion between the two homologous regions hr3 and hr4 of the genome of the Mexican (M) strain of Cydia pomonella granulovirus (CpGV-M) resulted in a loss of viral competitiveness, the function of these homologous regions was investigated. A CpGV-based bacmid (CpBAC) was constructed and mutants with deleted hr3 and hr4 palindromes (CpBAChr3/hr4KO) and a construct (CpBAChr3-kan-hr4) with physically separated hr3 and hr4 repeats were generated to investigate their involvement in in vivo replication. Based on median lethal concentration (LC50) and median survival time (ST50) of the mutant viruses vCpBAChr3/hr4KO and vCpBAChr3-kan-hr4 it was found that the infectivity of both mutants for codling moth Cydia pomonella L. (Lep.: Tortricidae) larvae was not influenced compared with the parental virus vCpBAC. Co-infection experiments with vCpBAChr3-kan-hr4 and vCpBAC using different virus ratios revealed that vCpBAChr3-kan-hr4 was efficiently out-competed by vCpBAC during in vivo replication. These findings suggested that the separation of hr3 and hr4 resulted in a replication disadvantage of the mutant similar to the observation made in previous co-infection experiments using the transposon-carrying mutant CpGV-MCp5 and WT CpGV-M. It was concluded that the palindromes hr3 and hr4 may play a non-essential but co-functional role in the replication of CpGV-M.

Baculovirus resistance in codling moth is virus isolate-dependent and the consequence of a mutation in viral gene pe38

The baculovirus Cydia pomonella granulovirus (CpGV) is widely applied as a biocontrol agent of codling moth. After field resistance of codling moth populations had been observed against the commercially used Mexican (M) isolate of CpGV, infection experiments of larvae of the resistant codling moth strain CpRR1 showed that several other naturally occurring CpGV isolates (I12, S, E2, and I07) from different geographic origins are still infectious to resistant CpRR1. Whole-genome sequencing and phylogenetic analyses of these geographic CpGV variants revealed that their genomes share only a single common difference from that of CpGV-M, which is a mutation coding for a repeat of 24 nucleotides within the gene pe38; this mutation results in an additional repeat of eight amino acids that appears to be inserted to PE38 of CpGV-M only. Deletion of pe38 from CpGV-M totally abolished virus infection in codling moth cells and larvae, demonstrating that it is an essential gene. When the CpGV-M deletion mutant was repaired with pe38 from isolate CpGV-S, which originated from the commercial product Virosoft and is infectious for the resistant codling moth strain CpRR1, the repaired CpGV-M mutant was found to be fully infectious for CpRR1. Repair using pe38 from CpGV-M restored infectivity for the virus in sensitive codling moth strains, but not in CpRR1. Therefore, we conclude that CpGV resistance of codling moth is directed to CpGV-M but not to other virus isolates. The viral gene pe38 is not only essential for the infectivity of CpGV but it is also the key factor in overcoming CpGV resistance in codling moth.

Homologous region 1 of Plutella xylostella granulovirus functions as an enhancer for early gene expression

Homologous regions (hrs) contained in genomes of baculoviruses have been shown to function as origins of viral DNA replication in alphabaculoviruses and betabaculoviruses, and as enhancers for early gene expression in alphabaculoviruses. The hr sequences of betabaculoviruses differ substantially from the ones of alphabaculoviruses. The enhancing property of betabaculovirus hrs has not been reported. In this study, transient assays were performed to investigate the effects of Plutella xylostella granulovirus (PlxyGV) hr1 on early gene expression of the virus. It was shown that hr1 stimulated reporter gene expression from the promoters of four early genes--ie1, dnapol, lef1, and lef9--independent of additional viral gene expression. The PlxyGV ie1 was shown to repress reporter gene expression from all four early gene promoters in a Trichoplusia ni cell line, both in the presence and absence of hr.

The genome and occlusion bodies of marine Penaeus monodon nudivirus (PmNV, also known as MBV and PemoNPV) suggest that it should be assigned to a new nudivirus genus that is distinct from the terrestrial nudiviruses

BACKGROUND

Penaeus monodon nudivirus (PmNV) is the causative agent of spherical baculovirosis in shrimp (Penaeus monodon). This disease causes significant mortalities at the larval stage and early postlarval (PL) stage and may suppress growth and reduce survival and production in aquaculture. The nomenclature and classification status of PmNV has been changed several times due to morphological observation and phylogenetic analysis of its partial genome sequence. In this study, we therefore completed the genome sequence and constructed phylogenetic trees to clarify PmNV's taxonomic position. To better understand the characteristics of the occlusion bodies formed by this marine occluded virus, we also compared the chemical properties of the polyhedrin produced by PmNV and the baculovirus AcMNPV (Autographa californica nucleopolyhedrovirus).

RESULTS

We used next generation sequencing and traditional PCR methods to obtain the complete PmNV genome sequence of 119,638 bp encoding 115 putative ORFs. Phylogenetic tree analysis showed that several PmNV genes and sequences clustered with the non-occluded nudiviruses and not with the baculoviruses. We also investigated the characteristics of PmNV polyhedrin, which is a functionally important protein and the major component of the viral OBs (occlusion bodies). We found that both recombinant PmNV polyhedrin and wild-type PmNV OBs were sensitive to acid conditions, but unlike the baculoviral OBs, they were not susceptible to alkali treatment.

CONCLUSIONS

From the viral genome features and phylogenetic analysis we conclude that PmNV is not a baculovirus, and that it should be assigned to the proposed Nudiviridae family with the other nudiviruses, but into a distinct new genus (Gammanudivirus).

Identification of a conserved non-protein-coding genomic element that plays an essential role in Alphabaculovirus pathogenesis

Highly homologous sequences 154-157 bp in length grouped under the name of "conserved non-protein-coding element" (CNE) were revealed in all of the sequenced genomes of baculoviruses belonging to the genus Alphabaculovirus. A CNE alignment led to the detection of a set of highly conserved nucleotide clusters that occupy strictly conserved positions in the CNE sequence. The significant length of the CNE and conservation of both its length and cluster architecture were identified as a combination of characteristics that make this CNE different from known viral non-coding functional sequences. The essential role of the CNE in the Alphabaculovirus life cycle was demonstrated through the use of a CNE-knockout Autographa californica multiple nucleopolyhedrovirus (AcMNPV) bacmid. It was shown that the essential function of the CNE was not mediated by the presumed expression activities of the protein- and non-protein-coding genes that overlap the AcMNPV CNE. On the basis of the presented data, the AcMNPV CNE was categorized as a complex-structured, polyfunctional genomic element involved in an essential DNA transaction that is associated with an undefined function of the baculovirus genome.

Complete genomic sequences and comparative analysis of Mamestra brassicae nucleopolyhedrovirus isolated in Korea

Mamestra brassicae nucleopolyhedrovirus-K1 (MabrNPV-K1) was isolated from naturally infected M. brassicae (Lepidoptera: Noctuidae) larvae in Korea. The full genome sequences of MabrNPV-K1 were determined, analysed and compared to those of other baculoviruses. The MabrNPV-K1 genome consisted of 152,710 bp and had an overall G + C content of 39.9%. Computer-assisted analysis predicted 158 open reading frames (ORFs) of 150 nucleotides or greater that showed minimal overlap. Two inhibitor of apoptosis (iap) and six baculovirus repeated ORFs were interspersed in the MabrNPV-K1 genome. The unique MabrNPV-K1 ORF133 was identified in the MabrNPV-K1 genome that was not previously reported in baculoviruses. The gene content and arrangement in MabrNPV-K1 had the highest similarity with those of Helicoverpa armigera MNPV (HearMNPV) and Mamestra configurata NPV-B (MacoNPV-B), and their shared homologous genes were 99% collinear. The MabrNPV-K1 genome contained four homologous repeat regions (hr1, hr2, hr3 and hr4) that accounted for 3.3% of the genome. The genomic positions of the four MabrNPV-K1 hr regions were conserved among those of HearMNPV and MacoNPV-B. The gene parity plot, percent identity of the gene homologues and a phylogenetic analysis suggested that these three viruses are closely related not only to each other but also to the same virus strains rather than different virus species.

Determination and analysis of the genome sequence of Spodoptera littoralis multiple nucleopolyhedrovirus

The Spodoptera littoralis multiple nucleopolyhedrovirus (SpliMNPV), a pathogen of the Egyptian cotton leaf worm S. littoralis, was subjected to sequencing of its entire DNA genome and bioassay analysis comparing its virulence to that of other baculoviruses. The annotated SpliMNPV genome of 137,998 bp was found to harbor 132 open reading frames and 15 homologous repeat regions. Four unique genes not present in SpltMNPV were identified, as were 14 genes that were absent or translocated by comparison. Bioassay analysis of experimentally infected Spodoptera frugiperda revealed an extended killing time for SpliMNPV as compared to S. frugiperda MNPV (SfMNPV), but a level of mortality similar to that caused by infection with SfMNPV and superior to that of Autographa californica MNPV (AcMNPV). Although extensive similarity was observed between the genome structure and predicted translation products of SpliMNPV and Spodoptera litura MNPV (SpltMNPV), genetic distances between isolates of SpliMNPV and SpltMNPV suggest that they are in fact different species of genus Alphabaculovirus.

Genome of Epinotia aporema granulovirus (EpapGV), a polyorganotropic fast killing betabaculovirus with a novel thymidylate kinase gene

BACKGROUND

Epinotia aporema (Lepidoptera: Tortricidae) is an important pest of legume crops in South America. Epinotia aporema granulovirus (EpapGV) is a baculovirus that causes a polyorganotropic infection in the host larva. Its high pathogenicity and host specificity make EpapGV an excellent candidate to be used as a biological control agent.

RESULTS

The genome of Epinotia aporema granulovirus (EpapGV) was sequenced and analyzed. Its circular double-stranded DNA genome is 119,082 bp in length and codes for 133 putative genes. It contains the 31 baculovirus core genes and a set of 19 genes that are GV exclusive. Seventeen ORFs were unique to EpapGV in comparison with other baculoviruses. Of these, 16 found no homologues in GenBank, and one encoded a thymidylate kinase. Analysis of nucleotide sequence repeats revealed the presence of 16 homologous regions (hrs) interspersed throughout the genome. Each hr was characterized by the presence of 1 to 3 clustered imperfect palindromes which are similar to previously described palindromes of tortricid-specific GVs. Also, one of the hrs (hr4) has flanking sequences suggestive of a putative non-hr ori. Interestingly, two more complex hrs were found in opposite loci, dividing the circular dsDNA genome in two halves. Gene synteny maps showed the great colinearity of sequenced GVs, being EpapGV the most dissimilar as it has a 20 kb-long gene block inversion. Phylogenetic study performed with 31 core genes of 58 baculoviral genomes suggests that EpapGV is the baculovirus isolate closest to the putative common ancestor of tortricid specific betabaculoviruses.

CONCLUSIONS

This study, along with previous characterization of EpapGV infection, is useful for the better understanding of the pathology caused by this virus and its potential utilization as a bioinsecticide.

Genetic and biological variation among nucleopolyhedrovirus isolates from the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae)

A PCR-based method was used to identify and distinguish among 40 uncharacterized nucleopolyhedrovirus (NPV) isolates from larvae of the moth Spodoptera frugiperda that were part of an insect virus collection. Phylogenetic analysis was carried out with sequences amplified from two strongly conserved loci (polh and lef-8) from the 40 isolates in the collection and from eight previously studied S. frugiperda NPV (SfMNPV) isolates. To further distinguish these isolates, analysis was also carried out with sequences from two less-conserved loci, hr4 and hr5. Phylogenetic inference from the sequence data could distinguish among several of the individual isolates and between different groups of isolates from Georgia (USA) and Colombia, South America. A stronger degree of bootstrap support for the phylogenetic trees was obtained with the hr4 and hr5 homologous repeat sequences. Sequencing and phylogenetic analysis detected a relatively high degree of larva-to-larva sequence divergence occurring among isolates of SfMNPV collected from the same field in Missouri, USA. Restriction endonuclease analysis of viral DNA from larvae infected with five isolates from Georgia, Missouri, Louisiana, Florida (USA), and Colombia allowed for further comparison with other previously reported isolates of SfMNPV. Bioassays with these five geographically distinct isolates detected minor differences in virulence. This study highlights the use of PCR to rapidly distinguish and characterize large numbers of historical baculovirus isolates from the same host using minimal quantities of material, and the use of sequences from homologous repeat regions to distinguish closely related isolates of the same NPV species.

Genomic sequencing and analyses of Lymantria xylina multiple nucleopolyhedrovirus

BACKGROUND

Outbreaks of the casuarina moth, Lymantria xylina Swinehoe (Lepidoptera: Lymantriidae), which is a very important forest pest in Taiwan, have occurred every five to 10 years. This moth has expanded its range of host plants to include more than 65 species of broadleaf trees. LyxyMNPV (L. xylina multiple nucleopolyhedrovirus) is highly virulent to the casuarina moth and has been investigated as a possible biopesticide for controlling this moth. LdMNPV-like virus has also been isolated from Lymantria xylina larvae but LyxyMNPV was more virulent than LdMNPV-like virus both in NTU-LY and IPLB-LD-652Y cell lines. To better understand LyxyMNPV, the nucleotide sequence of the LyxyMNPV DNA genome was determined and analysed.

RESULTS

The genome of LyxyMNPV consists of 156,344 bases, has a G+C content of 53.4% and contains 157 putative open reading frames (ORFs). The gene content and gene order of LyxyMNPV were similar to those of LdMNPV, with 151 ORFs identified as homologous to those reported in the LdMNPV genome. Two genes (Lyxy49 and Lyxy123) were homologous to other baculoviruses, and four unique LyxyMNPV ORFs (Lyxy11, Lyxy19, Lyxy130 and Lyxy131) were identified in the LyxyMNPV genome, including a gag-like gene that was not reported in baculoviruses. LdMNPV contains 23 ORFs that are absent in LyxyMNPV. Readily identifiable homologues of the gene host range factor-1 (hrf-1), which appears to be involved in the susceptibility of L. dispar to NPV infection, were not present in LyxyMNPV. Additionally, two putative odv-e27 homologues were identified in LyxyMNPV. The LyxyMNPV genome encoded 14 bro genes compared with 16 in LdMNPV, which occupied more than 8% of the LyxyMNPV genome. Thirteen homologous regions (hrs) were identified containing 48 repeated sequences composed of 30-bp imperfect palindromes. However, they differed in the relative positions, number of repeats and orientation in the genome compared to LdMNPV.

CONCLUSION

The gene parity plot analysis, percent identity of the gene homologues and a phylogenetic analysis suggested that LyxyMNPV is a Group II NPV that is most closely related to LdMNPV but with a highly distinct genomic organisation.

André Paillot (1885-1944): his work lives on

André Paillot (1885-1944) is one of the true founders of insect pathology. His applied and basic research culminated in the first identification granulovirus infection in caterpillars of the cabbage white, Pieris brassicae, and in the characterization of many other insect diseases. This paper reviews the life-time achievements of one of the most remarkable invertebrate pathologists in the first half of the last century.

Structural divergence among genomes of closely related baculoviruses and its implications for baculovirus evolution

Baculoviruses are members of a large, well-characterized family of dsDNA viruses that have been identified from insects of the orders Lepidoptera, Hymenoptera, and Diptera. Baculovirus genomes from different virus species generally exhibit a considerable degree of structural diversity. However, some sequenced baculovirus genomes from closely related viruses are structurally very similar and share overall nucleotide sequence identities in excess of 95%. This review focuses on the comparative analysis of partial and complete nucleotide sequences from two groups of closely related baculoviruses with broad host ranges: (a) group I multiple nucleopolyhedroviruses (MNPVs) from a cluster including Autographa californica (Ac)MNPV, Rachiplusia ou (Ro)MNPV, and Plutella xylostella (Plxy)MNPV; and (b) granuloviruses (GVs) from a cluster including Xestia c-nigrum (Xecn)GV and Helicoverpa armigera (Hear)GV. Even though the individual viruses in these clusters share high nucleotide sequence identities, a significant degree of genomic rearrangement (in the form of insertions, deletions, and homologous recombination resulting in allelic replacement) is evident from alignments of their genomes. These observations suggest an important role for recombination in the early evolution and biological characteristics of baculoviruses of these two groups.

Diversity and evolution of the Cydia pomonella granulovirus

Eight new field isolates of Cydia pomonella granulovirus (CpGV) originating in Iran and Georgia and one English CpGV isolate were analysed for restriction fragment length polymorphisms (RFLPs) and by partial genome amplification and sequencing. According to the observed RFLPs, most of the predominant genotypes of these isolates could be assigned to those present in previously found isolates originating from Mexico (CpGV-M), England (CpGV-E) and Russia (CpGV-R). We suggest that these isolates should be designated genome A, B and C types, respectively. A fourth genome type was identified in three isolates and is designated D type. The isolates with A, B and D type genomes contained four open reading frames (ORFs) (ORF63-ORF66) not present in C type genomes. The lack of these ORFs in other granuloviruses suggests that the C type genome is evolutionarily ancestral to the other genome types. The B and D type genomes contained an additional insertion of a non-protein coding region of 0.7 kb, which was at different genome locations. Analysis of the partial gene sequences of late expression factor 8 (lef-8), lef-9 and polyhedrin/granulin (polh/gran) genes revealed single nucleotide polymorphisms (SNPs) that corresponded to the RFLP types. Phylogenetic analyses based on these SNPs corroborated the proposed ancestry of the C type genome. C type viruses were also less virulent to neonate codling moth larvae than the other virus types. In conclusion, the known diversity of CpGV isolates can be described by four major genome types, which appear to exist in different isolates as genotype mixtures.

Genomic sequence and biological characterization of a nucleopolyhedrovirus isolated from the summer fruit tortrix, Adoxophyes orana

Adoxophyes orana nucleopolyhedrovirus (AdorNPV) was isolated from overwintering larvae from an orchard in the UK. The nucleotide sequence of the AdorNPV DNA genome was determined and analysed. The genome contains 111724 bp and has a G+C content of 35.0 mol%. The analysis predicted 121 ORFs of 150 nt or larger. Of these putative genes, 118 were homologous to genes identified previously in the Adoxophyes honmai nucleopolyhedrovirus (AdhoNPV) genome (83.3-100 % aa identity), and three AdorNPV ORFs were unique. There were four small homologous regions that consisted of a similar core sequence and at the same relative positions in the genome as AdhoNPV, but they differed in the number of repeats and orientation. Some genes that have been reported to have major roles in baculovirus biology were either absent or truncated in the AdorNPV genome. These included chitinase, which is involved in the liquefaction of the host, and the C-terminal of the ecdysteroid UDP-glucosyltransferase (egt) protein, which was truncated by 149 aa compared with AdhoNPV, with essential amino acids absent. The AdorNPV genome encoded two inhibitor of apoptosis (iap) genes compared with three in AdhoNPV and three bro genes compared with four in AdhoNPV. The susceptibility of A. orana larvae to AdorNPV was evaluated in laboratory bioassays using inoculation by microdroplet feeding and applied dose assays. LD50 for neonates was 56 occlusion bodies rising to 2.3x10(4) for fifth instar larvae. Median survival time values using an LD80 dose were 8.8 days for neonates and 7.0 days for fifth instar larvae.

Genomic sequence analysis of the Illinois strain of the Agrotis ipsilon multiple nucleopolyhedrovirus

The Agrotis ipsilon multiple nucleopolyhedrovirus (AgipMNPV) is a group II nucleopolyhedrovirus (NPV) from the black cutworm, A. ipsilon, with potential as a biopesticide to control infestations of cutworm larvae. The genome of the Illinois strain of AgipMNPV was completely sequenced. The AgipMNPV genome was 155,122 nt in size and contained 163 open reading frames (ORFs), including 61 ORFs found among all lepidopteran baculoviruses sequenced to date. Phylogenetic inference placed AgipMNPV in a clade with group II NPVs isolated from larvae of Agrotis and Spodoptera species. Though closely related to the Agrotis segetum NPV (AgseNPV), AgipMNPV was found to be missing 15 ORFs present in the AgseNPV genome sequence, including two of the three AgseNPV enhancin genes. Remarkably few polymorphisms were identified in the AgipMNPV sequence even though an uncloned field isolate of this virus was sequenced. A genotype characterized by a 128-bp deletion in the ecdysteroid UDP-glucosyltransferase gene (egt) was identified in the AgipMNPV field isolate and among clonal isolates of AgipMNPV. The deletion in egt was not associated with differences in budded virus or occluded virus production among AgipMNPV clones in cell culture.

The origins of replication of granuloviruses

The genomes of eight granuloviruses (GVs), have been analyzed for the presence of homologous regions (hrs) that may act as origins of replication. Thirteen 74-76-bp palindromes within 11 hrs have previously been identified in the Cydia pomonella GV (CpGV) genome and found to replicate in an infection-dependent DNA replication assay. We report a further palindrome within one of the hrs, which was found to replicate, bringing the total to 14 palindromes. We also report imperfect palindromes, with similar 13-bp end sequences to the CpGV palindromes, within the Adoxophyes orana GV, Cryptophlebia leucotreta GV (CrleGV), Choristoneura occidentalis GV and Phthorimaea operculella GV genomes. No hrs were detected in Agrotis segetum GV, and no additional hrs or palindromes, other than those published, were detected in the Plutella xylostella GV and Xestia c-nigrum GV genomes. Several putative hrs from the GVs were tested for replication in C. pomonella cells using a CpGV-dependent replication assay. Two CrleGV hrs were found to replicate at a low level.

A bacmid approach to the genetic manipulation of granuloviruses

A Cydia pomonella granulovirus (CpGV) bacmid has been constructed, which allows rapid and efficient production of recombinant baculoviruses in Escherichia coli. An 8.6kbp bacterial DNA cassette derived from the AcMNPV Bac-to-Bac system was ligated into a unique PacI restriction site within an intergenic region flanking the DNA ligase gene of the CpGV genome. The CpGV bacmids produced in E. coli were transfected into a CpGV-permissive C. pomonella cell line and the transfected cells fed to larvae to amplify the virus. The enhanced green fluorescent protein (EGFP) gene under the constitutive Drosophila heat-shock promoter was transposed into the mini-attTn7 transposition site, using a modified pFASTBAC donor plasmid, to generate a recombinant CpGV bacmid which caused infected larvae to glow under UV light. Targeted homologous recombination was also achieved in a recombinant proficient E. coli strain (BJ5183). A chloramphenicol acetyl transferase (CAT) gene replaced the cathepsin (v-cath) gene in the bacmid to produce a v-cath-deletion mutant. This is the first published report of a granulovirus bacmid, which will allow easy manipulation of the CpGV genome, enabling future studies on granulovirus genes and biology.

Genomic sequence analysis of a granulovirus isolated from the Old World bollworm, Helicoverpa armigera

The genome of a granulovirus isolated from the Old World bollworm, Helicoverpa armigera, was completely sequenced. The size of the Helicoverpa armigera granulovirus (HearGV) genome is 169,794 nt containing 179 open reading frames (ORFs), making it the second largest baculovirus genome analyzed to date. The genomes of HearGV and the Xestia c-nigrum GV (XecnGV) exhibit extensive sequence similarity and co-linearity, with both genomes containing the same nine homologous regions (hrs) with conserved structure and locations and sharing 167 open reading frames (ORFs). Phylogenetic inference and pairwise analysis of Kimura-2-parameter nucleotide distances for the lef-8, lef-9, and granulin genes indicate that HearGV is part of a cluster of granuloviruses typified by XecnGV. The HearGV genome contains all 62 ORFs found in common among other fully sequenced lepidopteran baculovirus genomes, as well as seven ORFs unique to HearGV. In addition, HearGV and XecnGV genomes share 20 ORFs not found among other baculovirus genomes sequenced to date. In addition to possessing ten ORFs with sequence similarity to baculovirus repeated ORFs (bro), the HearGV genome contains members of two other gene families with homologues in ascovirus, nucleopolyhedrovirus, and entomopoxvirus genomes. Alignment of the HearGV and XecnGV genome sequences revealed that HearGV is missing approximately 16.6 kbp of XecnGV-homologous sequence and contains approximately 8.2 kbp of sequence not found in the XecnGV genome.