Metagenomic Analysis of Viromes of Aedes Mosquitoes across India.

Metagenomic analysis of Aedes aegypti and Ae. albopictus mosquitoes from diverse geographical regions of India revealed the presence of several insect viruses of human interest. Most abundant reads found in Ae. aegypti mosquitoes were of Phasi Charoen-like virus (PCLV), Choristoneura fumiferana granulovirus (CfGV), Cell fusing agent virus (CFAV), and Wenzhou sobemo-like virus 4 (WSLV4), whereas WSLV4 and CfGV constituted the highest percentage of reads in Ae. albopictus viromes. Other reads that were of low percentage included Hubei mosquito virus 2 (HMV2), Porcine astrovirus 4 (PAstV4), and Wild Boar astrovirus (WBAstV). PCLV and CFAV, which were found to be abundant in Ae. aegypti viromes were absent in Ae. albopictus viromes. Among the viromes analyzed, Ae. aegypti sampled from Pune showed the highest percentage (79.82%) of viral reads, while Ae. aegypti mosquitoes sampled from Dibrugarh showed the lowest percentage (3.47%). Shamonda orthobunyavirus (SHAV), African swine fever virus (ASFV), Aroa virus (AROAV), and Ilheus virus (ILHV), having the potential to infect vertebrates, including humans, were also detected in both mosquito species, albeit with low read numbers. Reads of gemykibivirus, avian retrovirus, bacteriophages, herpesviruses, and viruses infecting protozoans, algae, etc., were also detected in the mosquitoes. A high percentage of reads in the Ae. albopictus mosquito samples belonged to unclassified viruses and warrant further investigation. The data generated in the present work may not only lead to studies to explain the influence of these viruses on the replication and transmission of viruses of clinical importance but also to find applications as biocontrol agents against pathogenic viruses.

CmHem, a hemolin-like gene identified from Cnaphalocrocis medinalis, involved in metamorphosis and baculovirus infection.

Background: As a member of the immunoglobulin superfamily, hemolins play a vital role in insect development and defense against pathogens. However, the innate immune response of hemolin to baculovirus infection varies among different insects. Methods and results: In this study, the hemolin-like gene from a Crambidae insect, Cnaphalocrocis medinalis, CmHem was cloned, and its role in insect development and baculovirus infection was analyzed. A 1,528 bp contig as potential hemolin-like gene of C. medinalis was reassembled from the transcriptome. Further, the complete hemolin sequence of C. medinalis (CmHem) was cloned and sequenced. The cDNA of CmHem was 1,515 bp in length and encoded 408 amino acids. The deduced amino acid of CmHem has relatively low identities (41.9-62.3%) to various insect hemolins. However, it contains four Ig domains similarity to other insect hemolins. The expression level of CmHem was the highest in eggs, followed by pupae and adults, and maintained a low expression level at larval stage. The synthesized siRNAs were injected into mature larvae, and the CmHem transcription decreased by 51.7%. Moreover, the abdominal somites of larvae became straightened, could not pupate normally, and then died. Infection with a baculovirus, C. medinalis granulovirus (CnmeGV), the expression levels of CmHem in the midgut and fat body of C. medinalis significantly increased at 12 and 24 h, respectively, and then soon returned to normal levels. Conclusions: Our results suggested that hemolin may be related to the metamorphosis of C. medinalis. Exposure to baculovirus induced the phased expression of hemolin gene in the midgut and fat body of C. medinalis, indicated that hemolin involved in the immune recognition of Crambidae insects to baculovirus.

Developments in the classification and nomenclature of arthropod-infecting large DNA viruses that contain pif genes.

Viruses of four families of arthropod-specific, large dsDNA viruses (the nuclear arthropod large DNA viruses, or NALDVs) possess homologs of genes encoding conserved components involved in the baculovirus primary infection mechanism. The presence of such homologs encoding per os infectivity factors (pif genes), along with their absence from other viruses and the occurrence of other shared characteristics, suggests a common origin for the viruses of these families. Therefore, the class Naldaviricetes was recently established, accommodating these four families. In addition, within this class, the ICTV approved the creation of the order Lefavirales for three of these families, whose members carry homologs of the baculovirus genes that code for components of the viral RNA polymerase, which is responsible for late gene expression. We further established a system for the binomial naming of all virus species in the order Lefavirales, in accordance with a decision by the ICTV in 2019 to move towards a standardized nomenclature for all virus species. The binomial species names for members of the order Lefavirales consist of the name of the genus to which the species belongs (e.g., Alphabaculovirus), followed by a single epithet that refers to the host species from which the virus was originally isolated. The common names of viruses and the abbreviations thereof will not change, as the format of virus names lies outside the remit of the ICTV.

The distribution and characteristic of two transposable elements in the genome of Cydia pomonella granulovirus and codling moth.

Baculoviruses are capable to acquire insect host transposable elements (TEs) in their genomes and are hypothesized as possible vectors of insect transposons between Lepidopteran species. Here, we investigated the host origin of two TEs, namely the Tc1/mariner-like element TCp3.2 and a 0.7 kbp insertion sequence (IS07), found in the genome of different isolates of Cydia pomonella granulovirus (CpGV), a member of the Betabaculovirus genus. The sequences of both TEs were searched for in the full genome sequence database of codling moth (CM, Cydia pomonella L.). A total of eleven TCp3.2 TE copies and 76 copies of the IS07 fragments were identified in the CM genome. These TEs were distributed over the 22 autosomes and the Z chromosome (chr1) of CM, except chr6, chr12, chr16, chr23, chr27 and the W chromosome (chr29). TCp3.2 copies with two transposase genes in opposite direction, representing a novel feature, were identified on chr10 and chr18. The TCp3.2 transposase was characterized by DD41D motif of classic Tc1/mariner transposons, consisting of DNA-binding domain, catalytic domain and nuclear localization signal (NLS). Transcription analyses of uninfected and CpGV-infected CM larvae suggested a doubling of the TCp3.2 transposase transcription rate in virus infected larvae. Furthermore, IS07 insertion into the CpGV genome apparently added new transcription initiation sites to the viral genome. The global analysis of the distribution of two TEs in the genome of CM addressed the influx of mobile TEs from CM to CpGV, a genetic process that contributes to the population diversity of baculoviruses.

Quantitative proteomic analysis of extracellular vesicles in response to baculovirus infection of a Trichoplusia ni cell line.

Due to its outstanding suitability to produce complex biopharmaceutical products including virus-like particles and subunit vaccines, the baculovirus/insect cell expression system has developed into a highly popular production platform in the biotechnological industry. For high productivity, virus-cell communication and an efficient spreading of the viral infection are crucial, and, in this context, extracellular vesicles (EVs) might play a significant role. EVs are small particles, utilized by cells to transfer biologically active compounds such as proteins, lipids as well as nucleic acids to recipient cells for intracellular communication. Studies in mammalian cells showed that the release of EVs is altered in response to infection with many viruses, ultimately either limiting or fostering infection spreading. In this study we isolated and characterized EVs, from both uninfected and baculovirus infected Tnms42 insect cells. Via quantitative proteomic analysis we identified more than 3000 T. ni proteins in Tnms42 cell derived EVs, of which more than 400 were significantly differentially abundant upon baculovirus infection. Subsequent gene set enrichment analysis revealed a depletion of proteins related to the extracellular matrix in EVs from infected cultures. Our findings show a significant change of EV protein cargo upon baculovirus infection, suggesting a major role of EVs as stress markers. Our study might serve in designing new tools for process monitoring and control to further improve biopharmaceutical production within the baculovirus/insect cell expression system.

Regulatory Mechanisms, Protein Expression and Biological Activity of Photolyase Gene from Spodoptera littoralis Granulovirus Genome.

One of the most important factor that affects the efficient using of baculoviruses as a biopesticide is their sensitivity to UV irradiation. In this study, a photolyase gene (phr) of 1.4 kbp DNA fragment was cloned and characterized from Spodoptera littoralis granulovirus, an Egyptian isolate (SpliGV-EG1). A sequence of 466 amino acid were deduced when the gene was completely sequenced with a predicted molecular mass of ~ 55 kDa. Transcriptional regulation analyses revealed that phr transcripts were detected early at 6-h post-infection (hpi) and remained detectable until 72 hpi, suggesting their transcriptional regulation from a putative early promoter motif. An approximately ~ 55 kDa protein fragment was expressed from phr-induced bacterial culture and detected by SDS-PAGE and western blotting. In addition, direct exposure to UV irradiation resulted in a twofold decrease in SpliGV-EG1 occlusion bodies activation compared with Spodoptera littoralis nucleopolyhedrovirus (SpliNPV) occlusion bodies which decreased with about 129-fold after exposure to UV irradiation based on median lethal concentration value (LC50). The obtained results suggested that the presence of photolyase gene possibly alters the inactivation of SpliGV-EG1-occluded bodies by UV irradiation. These results support the role and application of the photolyase protein to improve the damaged DNA repair mechanism as well as resistance of SpliGV to UV light inactivation.

Genome analysis of Psilogramma increta granulovirus and its intrapopulation diversity.

The complete genome of Psilogramma increta granulovirus (PsinGV), isolated from P. increta (Lepidoptera: Sphingidae), was ultra-deep sequenced with a Novaseq PE150 platform and de novo assembled and annotated. The PsinGV genome is a circular double-stranded DNA, 103,721 bp in length, with a G+C content of 33.0%, the third lowest G+C content in present sequenced baculoviruses. It encodes 123 putative open reading frames, including 38 baculovirus core genes, 42 lepidopteran baculovirus conserved genes, 38 betabaculovirus conserved genes, and 5 genes unique to PsinGV. Meanwhile, 3 homologous repeated regions with the core sequence TTGCAA and 3 direct repeated sequences were identified within the PsinGV genome. Kimura two-parameters distance analysis confirmed that Psilogramma increta granulovirus is a representative of a prospective new species of the genus Betabaculovirus. Phylogenetic analysis based on the baculovirus core genes showed that PsinGV is closely related to Choristoneura fumiferana granulovirus, Clostera anastomosis granulovirus-B, and Erinnyis ello granulovirus. These four species therefore share a common ancestor residing in the Betabaculovirus genus. The genome of the PsinGV isolate contained two p10 copies: p10 and p10-2. Phylogenetic reconstruction of P10 suggests a transfer event of the p10-2 gene from an alphabaculovirus to the aforementioned common ancestor. Analysis of genomic diversity showed that 203 intrahost variants, including 182 single nucleotide variants and 21 short insertions/deletions, are present within the PsinGV isolate. Meanwhile, allele frequency indicated that the isolate contains three major genotypes, implying the archived isolate consists of several P. increta carcasses killed by PsinGV with different genotypes.

MAPK-mediated transcription factor GATAd contributes to Cry1Ac resistance in diamondback moth by reducing PxmALP expression.

The benefits of biopesticides and transgenic crops based on the insecticidal Cry-toxins from Bacillus thuringiensis (Bt) are considerably threatened by insect resistance evolution, thus, deciphering the molecular mechanisms underlying insect resistance to Bt products is of great significance to their sustainable utilization. Previously, we have demonstrated that the down-regulation of PxmALP in a strain of Plutella xylostella (L.) highly resistant to the Bt Cry1Ac toxin was due to a hormone-activated MAPK signaling pathway and contributed to the resistance phenotype. However, the underlying transcriptional regulatory mechanism remains enigmatic. Here, we report that the PxGATAd transcription factor (TF) is responsible for the differential expression of PxmALP observed between the Cry1Ac susceptible and resistant strains. We identified that PxGATAd directly activates PxmALP expression via interacting with a non-canonical but specific GATA-like cis-response element (CRE) located in the PxmALP promoter region. A six-nucleotide insertion mutation in this cis-acting element of the PxmALP promoter from the resistant strain resulted in repression of transcriptional activity, affecting the regulatory performance of PxGATAd. Furthermore, silencing of PxGATAd in susceptible larvae reduced the expression of PxmALP and susceptibility to Cry1Ac toxin. Suppressing PxMAP4K4 expression in the resistant larvae transiently recovered both the expression of PxGATAd and PxmALP, indicating that the PxGATAd is a positive responsive factor involved in the activation of PxmALP promoter and negatively regulated by the MAPK signaling pathway. Overall, this study deciphers an intricate regulatory mechanism of PxmALP gene expression and highlights the concurrent involvement of both trans-regulatory factors and cis-acting elements in Cry1Ac resistance development in lepidopteran insects.

Establishment, verification and application of rapid detection of baculovirus infectious titer by flow cytometry.

Titer detection of baculovirus usually is time-consuming. It is important to establish a rapid detection method for baculovirus titer. In this report, Staining of cells with a fluorescently labeled anti-gp64 antibody allows for identification of infected insect cells. By inoculating cultures with a series of log dilutions of virus, and staining of the cultures 13-22 h post inoculation, the ratio of infected to un-infected insect cells can be determined by flow cytometry(FCM). Statistical analysis of the percentage of infected cells in the virus dilution series enables accurate infectious titer determination. The culture time, cell growth state, the concentration of GP64-APC antibody and the concentration of inactivated FBS in diluent were optimized. The generality, repeatability and intermediate precision of the method were verified. The FCM method has the advantages of simplicity, accuracy, low cost and good repeatability.

Design, synthesis and biological activities of echinopsine derivatives containing acylhydrazone moiety.

Based on the broad-spectrum biological activities of echinopsine and acylhydrazones, a series of echinopsine derivatives containing acylhydrazone moieties have been designed, synthesized and their biological activities were evaluated for the first time. The bioassay results indicated that most of the compounds showed moderate to good antiviral activities against tobacco mosaic virus (TMV), among which echinopsine (I) (inactivation activity, 49.5 ± 4.4%; curative activity, 46.1 ± 1.5%; protection activity, 42.6 ± 2.3%) and its derivatives 1 (inactivation activity, 44.9 ± 4.6%; curative activity, 39.8 ± 2.6%; protection activity, 47.3 ± 4.3%), 3 (inactivation activity, 47.9 ± 0.9%; curative activity, 43.7 ± 3.1%; protection activity, 44.6 ± 3.3%), 7 (inactivation activity, 46.2 ± 1.6%; curative activity, 45.0 ± 3.7%; protection activity, 41.7 ± 0.9%) showed higher anti-TMV activity in vivo at 500 mg/L than commercial ribavirin (inactivation activity, 38.9 ± 1.4%; curative activity, 39.2 ± 1.8%; protection activity, 36.4 ± 3.4%). Some compounds exhibited insecticidal activities against Plutella xylostella, Mythimna separate and Spodoptera frugiperda. Especially, compounds 7 and 27 displayed excellent insecticidal activities against Plutella xylostell (mortality 67 ± 6% and 53 ± 6%) even at 0.1 mg/L. Additionally, most echinopsine derivatives exhibited high fungicidal activities against Physalospora piricola and Sclerotinia sclerotiorum.

Natural Coinfection between Novel Species of Baculoviruses in Spodoptera ornithogalli Larvae.

Spodoptera ornithogalli (Guenée) (Lepidoptera: Noctuidae) is an important pest in different crops of economic relevance in America. For its control, strategies that include chemicals are usually used; so, the description of entomopathogens would be very useful for the formulation of biopesticides. In this regard, two different baculoviruses affecting S. ornithogalli were isolated in Colombia, with one of them being an NPV and the other a GV. Ultrastructural, molecular, and biological characterization showed that both isolates possess the 38 core genes and are novel species in Baculoviridae, named as Spodoptera ornithogalli nucleopolyhedrovirus (SporNPV) and Spodoptera ornithogalli granulovirus (SporGV). The bioassays carried out in larvae of S. ornithogalli and S. frugiperda showed infectivity in both hosts but being higher in the first. In addition, it was observed that SporGV potentiates the insecticidal action of SporNPV (maximum value in ratio 2.5:97.5). Both viruses are individually infective but coexist in nature, producing mixed infections with a synergistic effect that improves the performance of the NPV and enables the transmission of the GV, which presents a slowly killing phenotype.

Cross-Resistance of the Codling Moth against Different Isolates of Cydia pomonella Granulovirus Is Caused by Two Different but Genetically Linked Resistance Mechanisms.

Cydia pomonella granulovirus (CpGV) is a widely used biological control agent of the codling moth. Recently, however, the codling moth has developed different types of field resistance against CpGV isolates. Whereas type I resistance is Z chromosomal inherited and targeted at the viral gene pe38 of isolate CpGV-M, type II resistance is autosomal inherited and targeted against isolates CpGV-M and CpGV-S. Here, we report that mixtures of CpGV-M and CpGV-S fail to break type II resistance and is expressed at all larval stages. Budded virus (BV) injection experiments circumventing initial midgut infection provided evidence that resistance against CpGV-S is midgut-related, though fluorescence dequenching assay using rhodamine-18 labeled occlusion derived viruses (ODV) could not fully elucidate whether the receptor binding or an intracellular midgut factor is involved. From our peroral and intra-hemocoel infection experiments, we conclude that two different (but genetically linked) resistance mechanisms are responsible for type II resistance in the codling moth: resistance against CpGV-M is systemic whereas a second and/or additional resistance mechanism against CpGV-S is located in the midgut of CpR5M larvae.

CpGV-M Replication in Type I Resistant Insects: Helper Virus and Order of Ingestion Are Important.

The genetic diversity of baculoviruses provides a sustainable agronomic solution when resistance to biopesticides seems to be on the rise. This genetic diversity promotes insect infection by several genotypes (i.e., multiple infections) that are more likely to kill the host. However, the mechanism and regulation of these virus interactions are still poorly understood. In this article, we focused on baculoviruses infecting the codling moth, Cydia pomonella: two Cydia pomonella granulovirus genotypes, CpGV-M and CpGV-R5, and Cryptophlebia peltastica nucleopolyhedrovirus (CrpeNPV). The influence of the order of ingestion of the virus genotypes, the existence of an ingestion delay between the genotypes and the specificity of each genotype involved in the success of multiple infection were studied in the case of Cydia pomonella resistance. To obtain a multiple infection in resistant insects, the order of ingestion is a key factor, but the delay for ingestion of the second virus is not. CrpeNPV cannot substitute CpGV-R5 to allow replication of CpGV-M.

High Water Density at Non-Ice-Binding Surfaces Contributes to the Hyperactivity of Antifreeze Proteins.

Antifreeze proteins (AFPs) can bind to ice nuclei thereby inhibiting their growth and their hydration shell is believed to play a fundamental role. Here, we use molecular dynamics simulations to characterize the hydration shell of four moderately-active and four hyperactive AFPs. The local water density around the ice-binding-surface (IBS) is found to be lower than that around the non-ice-binding surface (NIBS) and this difference correlates with the higher hydrophobicity of the former. While the water-density increase (with respect to bulk) around the IBS is similar between moderately-active and hyperactive AFPs, it differs around the NIBS, being higher for the hyperactive AFPs. We hypothesize that while the lower water density at the IBS can pave the way to protein binding to ice nuclei, irrespective of the antifreeze activity, the higher density at the NIBS of the hyperactive AFPs contribute to their enhanced ability in inhibiting ice growth around the bound AFPs.

Bacmid Expression of Granulovirus Enhancin En3 Accumulates in Cell Soluble Fraction to Potentiate Nucleopolyhedrovirus Infection.

Enhancins are metalloproteinases that facilitate baculovirus infection in the insect midgut. They are more prevalent in granuloviruses (GVs), constituting up to 5% of the proteins of viral occlusion bodies (OBs). In nucleopolyhedroviruses (NPVs), in contrast, they are present in the envelope of the occlusion-derived virions (ODV). In the present study, we constructed a recombinant Autographa californica NPV (AcMNPV) that expressed the Trichoplusia ni GV (TnGV) enhancin 3 (En3), with the aim of increasing the presence of enhancin in the OBs or ODVs. En3 was successfully produced but did not localize to the OBs or the ODVs and accumulated in the soluble fraction of infected cells. As a result, increased OB pathogenicity was observed when OBs were administered in mixtures with the soluble fraction of infected cells. The mixture of OBs and the soluble fraction of Sf9 cells infected with BacPhEn3 recombinant virus was ~3- and ~4.7-fold more pathogenic than BacPh control OBs in the second and fourth instars of Spodoptera exigua, respectively. In contrast, when purified, recombinant BacPhEn3 OBs were as pathogenic as control BacPh OBs. The expression of En3 in the soluble fraction of insect cells may find applications in the development of virus-based insecticides with increased efficacy.

Gene expression patterns of Cydia pomonella granulovirus in codling moth larvae revealed by RNAseq analysis.

The Cydia pomonella granulovirus (CpGV) has been used as a biological control agent of codling moth (Cydia pomonella), a severe global pest on pome fruit. Despite the economic importance, our knowledge of its molecular biology is still limited and a detailed picture of its gene expression is still missing. Here, we sequenced the transcriptome of codling moth larvae infected with the Mexican isolate CpGV-M and analyzed the expression of viral genes at 12, 48, and 96 h post infection (hpi). The results showed that two genes (p6.9 and pp31/39K) related to DNA binding of virus production, were highly expressed at 48 and 96 hpi. From 48 to 96 hpi, the expression of genes associated with virus replication and dissemination decreased, whereas the expression of genes related to infectious virion production and per os infectivity increased. This study provides a comprehensive view of CpGV gene expression patterns in host larvae.

Transcriptome of Cydia pomonella granulovirus in susceptible and type I resistant codling moth larvae.

The baculovirus Cydia pomonella granulovirus (CpGV) is a biocontrol agent used worldwide against the codling moth (CM), Cydia pomonella L., a severe pest in organic and integrated pome fruit production. Its successful application is increasingly challenged by the occurrence of CM populations resistant to commercial CpGV products. Whereas three types (I-III) of CpGV resistance have been identified, type I resistance compromising the efficacy of CpGV-M, the so-called Mexican isolate of CpGV, is assumed to be the most widely distributed resistance type in Central Europe. Despite the wide use of CpGV products as biocontrol agents, little information is available on gene-expression levels in CM larvae. In this study, the in vivo transcriptome of CpGV-M infecting susceptible (CpS) and resistant (CpRR1) CM larvae was analysed at 24, 48, 72, 96 and 120 hours post infection in the midgut and fat body tissue by using a newly developed microarray covering all ORFs of the CpGV genome. According to their transcript abundance, the CpGV genes were grouped into four temporal clusters to which groups of known and unknown function could be assigned. In addition, sets of genes differentially expressed in the midgut and fat body were found in infected susceptible CpS larvae. For the resistant CpRR1 larvae treated with CpGV-M, viral entry in midgut cells could be confirmed from onset but a significantly reduced gene expression, indicating that type I resistance is associated with a block of viral gene transcription and replication.

Granulovirus PK-1 kinase activity relies on a side-to-side dimerization mode centered on the regulatory αC helix.

The life cycle of Baculoviridae family insect viruses depends on the viral protein kinase, PK-1, to phosphorylate the regulatory protein, p6.9, to induce baculoviral genome release. Here, we report the crystal structure of Cydia pomenella granulovirus PK-1, which, owing to its likely ancestral origin among host cell AGC kinases, exhibits a eukaryotic protein kinase fold. PK-1 occurs as a rigid dimer, where an antiparallel arrangement of the αC helices at the dimer core stabilizes PK-1 in a closed, active conformation. Dimerization is facilitated by C-lobe:C-lobe and N-lobe:N-lobe interactions between protomers, including the domain-swapping of an N-terminal helix that crowns a contiguous ß-sheet formed by the two N-lobes. PK-1 retains a dimeric conformation in solution, which is crucial for catalytic activity. Our studies raise the prospect that parallel, side-to-side dimeric arrangements that lock kinase domains in a catalytically-active conformation could function more broadly as a regulatory mechanism among eukaryotic protein kinases.

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.

Comparative genomic analysis of three geographical isolates from China reveals high genetic stability of Plutella xylostella granulovirus.

In this study, the genomes of three Plutella xylostella granulovirus (PlxyGV) isolates, PlxyGV-W and PlxyGV-Wn from near Wuhan and PlxyGV-B from near Beijing, China were completely sequenced and comparatively analyzed to investigate genetic stability and diversity of PlxyGV. PlxyGV-W, PlxyGV-B and PlxyGV-Wn consist of 100,941bp, 100,972bp and 100,999bp in length with G + C compositions of 40.71-40.73%, respectively, and share nucleotide sequence identities of 99.5-99.8%. The three individual isolates contain 118 putative protein-encoding ORFs in common. PlxyGV-W, PlxyGV-B and PlxyGV-Wn have ten, nineteen and six nonsynonymous intra isolate nucleotide polymorphisms (NPs) in six, fourteen and five ORFs, respectively, including homologs of five DNA replication/late expression factors and two per os infectivity factors. There are seventeen nonsynonymous inter isolate NPs in seven ORFs between PlxyGV-W and PlxyGV-B, seventy three nonsynonymous NPs in forty seven ORFs between PlxyGV-W and PlxyGV-Wn, seventy seven nonsynonymous NPs in forty six ORFs between PlxyGV-B and PlxyGV-Wn. Alignment of the genome sequences of nine PlxyGV isolates sequenced up to date shows that the sequence homogeneity between the genomes are over 99.4%, with the exception of the genome of PlxyGV-SA from South Africa, which shares a sequence identity of 98.6-98.7% with the other ones. No events of gene gain/loss or translocations were observed. These results suggest that PlxyGV genome is fairly stable in nature. In addition, the transcription start sites and polyadenylation sites of thirteen PlxyGV-specific ORFs, conserved in all PlxyGV isolates, were identified by RACE analysis using mRNAs purified from larvae infected by PlxyGV-Wn, proving the PlxyGV-specific ORFs are all genuine genes.