Whole-genome sequence of Enterobacter hormaechei, isolate jjbc recovered from the gut of Plutella xylostella feeding on cabbage.

We present the whole-genome sequence of Enterobacter hormaechei (previously Enterobacter cloacae) obtained from long and short reads. It is a dominant gut symbiont of the notorious crop pest Plutella xylostella, highly prevalent in lepidopteran midguts and a useful model for the evolution of resistance to antimicrobials.

Transcriptome analysis of CpGV in midguts of type II resistant codling moth larvae and identification of contaminant infections by SNP mapping of RNA-Seq data.

Various isolates of the Cydia pomonella granulovirus (CpGV) are used as insect pest control agents against codling moth (CM, Cydia pomonella L.), a predominant pest in apple orchards. Three different types (I-III) of dominantly inherited field resistance of CM larvae to CpGV have been recently identified. In this study, transcription of virus genes in midgut cells of type II-resistant CM larvae infected with different CpGV isolates, i.e., CpGV-M and CpGV-S (both prone to type II resistance) as well as CpGV-E2 (breaking type II resistance) was determined by strand-specific RNA sequencing (RNA-Seq) at an early infection stage (72 h post infection). Based on principal component analysis of read counts and the quantitative distribution of single nucleotide polymorphisms (SNPs) in the RNA-Seq data, a bioinformatics analysis pipeline was developed for an a posteriori identification of the infective agents. We report that (i) identification of infective agent is crucial, especially in in vivo infection experiments, when activation of covert virus infections is a possibility, (ii) no substantial difference between CpGV-M and CpGV-S transcription was found in type II-resistant CM larvae despite a different resistance mechanism, (iii) the transcription level of CpGV-M and CpGV-S was much lower than that of CpGV-E2, and (iv) orf59 (sod), orf89 (pif-6), orf92 (p18), and orf137 (lef-10) were identified as significantly downregulated genes in resistance-prone isolates CpGV-M and CpGV-S. For type II resistance of CM larvae, we conclude that CpGV-M and CpGV-S are both able to enter midgut cells, but viral transcription is significantly impaired at an early stage of infection compared to the resistance-breaking isolate CpGV-E2. IMPORTANCE: CpGV is a highly virulent pathogen of codling moth, and it has been developed into one of the most successful commercial baculovirus biocontrol agents for pome fruit production worldwide. The emergence of field resistance in codling moth to commercial CpGV products is a threat toward the sustainable use of CpGV. In recent years, different types of resistance (type I-III) were identified. For type II resistance, very little is known regarding the infection process. By studying the virus gene expression patterns of different CpGV isolates in midguts of type II-resistant codling moth larvae, we found that the type II resistance mechanism is most likely based on intracellular factors rather than a receptor component. By applying SNP mapping of the RNA-Seq data, we further emphasize the importance of identifying the infective agents in in vivo experiments when activation of a covert infection cannot be excluded.

Pathogens of the oak processionary moth Thaumetopoea processionea: Developing a user-friendly bioassay system and metagenome analyses for microorganisms.

The oak processionary moth (OPM) Thaumetopoea processionea is a pest of oak trees and poses health risks to humans due to the urticating setae of later instar larvae. For this reason, it is difficult to rear OPM under laboratory conditions, carry out bioassays or examine larvae for pathogens. Biological control targets the early larval instars and is based primarily on commercial preparations of Bacillus thuringiensis ssp. kurstaki (Btk). To test the entomopathogenic potential of other spore-forming bacteria, a user-friendly bioassay system was developed that (i) applies bacterial spore suspensions by oak bud dipping, (ii) targets first instar larvae through feeding exposure and (iii) takes into account their group-feeding behavior. A negligible mortality in the untreated control proved the functionality of the newly established bioassay system. Whereas the commercial Btk HD-1 strain was used as a bioassay standard and confirmed as being highly efficient, a Bacillus wiedmannii strain was ineffective in killing OPM larvae. Larvae, which died during the infection experiment, were further subjected to Nanopore sequencing for a metagenomic approach for entomopathogen detection. It further corroborated that B.wiedmannii was not able to infect and establish in OPM, but identified potential insect pathogenic species from the genera Serratia and Pseudomonas.

Haplotype determination of the Bombyx mori nucleopolyhedrovirus by Nanopore sequencing and linkage of single nucleotide variants.

Naturally occurring isolates of baculoviruses, such as the Bombyx mori nucleopolyhedrovirus (BmNPV), usually consist of numerous genetically different haplotypes. Deciphering the different haplotypes of such isolates is hampered by the large size of the dsDNA genome, as well as the short read length of next generation sequencing (NGS) techniques that are widely applied for baculovirus isolate characterization. In this study, we addressed this challenge by combining the accuracy of NGS to determine single nucleotide variants (SNVs) as genetic markers with the long read length of Nanopore sequencing technique. This hybrid approach allowed the comprehensive analysis of genetically homogeneous and heterogeneous isolates of BmNPV. Specifically, this allowed the identification of two putative major haplotypes in the heterogeneous isolate BmNPV-Ja by SNV position linkage. SNV positions, which were determined based on NGS data, were linked by the long Nanopore reads in a Position Weight Matrix. Using a modified Expectation-Maximization algorithm, the Nanopore reads were assigned according to the occurrence of variable SNV positions by machine learning. The cohorts of reads were de novo assembled, which led to the identification of BmNPV haplotypes. The method demonstrated the strength of the combined approach of short- and long-read sequencing techniques to decipher the genetic diversity of baculovirus isolates.

Advancing pathogen surveillance by nanopore sequencing and genotype characterization of Acheta domesticus densovirus in mass-reared house crickets.

Rapid and reliable detection of pathogens is crucial to complement the growing industry of mass-reared insects, in order to safeguard the insect colonies from outbreak of diseases, which may cause significant economic loss. Current diagnostic methods are mainly based on conventional PCR and microscopic examination, requiring prior knowledge of disease symptoms and are limited to identifying known pathogens. Here, we present a rapid nanopore-based metagenomics approach for detecting entomopathogens from the European house cricket (Acheta domesticus). In this study, the Acheta domesticus densovirus (AdDV) was detected from diseased individuals using solely Nanopore sequencing. Virus reads and genome assemblies were obtained within twenty-four hours after sequencing. Subsequently, due to the length of the Nanopore reads, it was possible to reconstruct significantly large parts or even the entire AdDV genome to conduct studies for genotype identification. Variant analysis indicated the presence of three AdDV genotypes within the same house cricket population, with association to the vital status of the diseased crickets. This contrast provided compelling evidence for the existence of non-lethal AdDV genotypes. These findings demonstrated nanopore-based metagenomics sequencing as a powerful addition to the diagnostic tool kit for routine pathogen surveillance and diagnosis in the insect rearing industry.

A practical guide and Galaxy workflow to avoid inter-plasmidic repeat collapse and false gene loss in Unicycler's hybrid assemblies.

Generating complete, high-quality genome assemblies is key for any downstream analysis, such as comparative genomics. For bacterial genome assembly, various algorithms and fully automated pipelines exist, which are free-of-charge and easily accessible. However, these assembly tools often cannot unambiguously resolve a bacterial genome, for example due to the presence of sequence repeat structures on the chromosome or on plasmids. Then, a more sophisticated approach and/or manual curation is needed. Such modifications can be challenging, especially for non-bioinformaticians, because they are generally not considered as a straightforward process. In this study, we propose a standardized approach for manual genome completion focusing on the popular hybrid assembly pipeline Unicycler. The provided Galaxy workflow addresses two weaknesses in Unicycler's hybrid assemblies: (i) collapse of inter-plasmidic repeats and (ii) false loss of single-copy sequences. To demonstrate and validate how to detect and resolve these assembly errors, we use two genomes from the Bacillus cereus group. By applying the proposed pipeline following an automated assembly, the genome sequence quality can be significantly improved.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176.

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

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.

Identification of a new nucleopolyhedrovirus isolated from the olive leaf moth, Palpita vitrealis, from two locations in Egypt.

The olive leaf moth (jasmine moth), Palpita vitrealis (Lepidoptera: Crambidae), is an important insect pest of olives in several Mediterranean countries. A new alphabaculovirus was isolated from diseased larvae of P. vitrealis in Egypt, first in Giza in spring 2005 and again in Marsa Matrouh in 2019.The larvae exhibited typical symptoms of a baculovirus infection. Light and scanning electron microscopy studies revealed polyhedral occlusion bodies. Transmission electron microscopy of ultrathin sections of purified OBs revealed virions with multiple embedded nucleocapsids. The identity of the two virus isolates was confirmed by sequencing the partial polyhedrin and lef-8 genes, and sequence comparison suggested a relationship to group I alphabaculoviruses. Therefore, this virus was termed Palpita vitrealis nucleopolyhedrovirus (PaviNPV). Whole genome sequencing of the PaviNPV isolate from Giza (Gz05) revealed a genome of 117,533 bp, 131 open reading frames (ORFs) and four homologous repeat (hr) regions. Phylogenetic reconstruction and genetic distance analyses using 38 core genes indicated that PaviNPV should be considered to belong to a novel species within the genus Alphabaculovirus. In bioassays, PaviNPV was highly virulent against second-instar larvae of P. vitrealis. The study reports a novel baculovirus that might have potential as a biological control agent of the olive leaf moth.

Viruses of the Fall Armyworm Spodoptera frugiperda: A Review with Prospects for Biological Control.

The fall armyworm (FAW), Spodoptera frugiperda, is a native pest species in the Western hemisphere. Since it was first reported in Africa in 2016, FAW has spread throughout the African continent and is now also present in several countries in Asia as well as Australia. The invasion of FAW in these areas has led to a high yield reduction in crops, leading to huge economic losses. FAW management options in the newly invaded areas are limited and mainly rely on the use of synthetic pesticides. Since there is a risk of resistance development against pesticides in addition to the negative environmental and human health impacts, other effective, sustainable, and cost-efficient control alternatives are desired. Insect pathogenic viruses fulfil these criteria as they are usually effective and highly host-specific with no significant harmful effect on beneficial insects and non-target organisms. In this review, we discuss all viruses known from FAW and their potential to be used for biological control. We specifically focus on baculoviruses and describe the recent advancements in the use of baculoviruses for biological control in the native geographic origin of FAW, and their potential use in the newly invaded areas. Finally, we identify current knowledge gaps and suggest new avenues for productive research on the use of viruses as a biopesticide against FAW.

Genome Sequence of a Spodoptera frugiperda Multiple Nucleopolyhedrovirus Isolated from Fall Armyworm (Spodoptera frugiperda) in Nigeria, West Africa.

We report the entire genome sequence of an isolate of Spodoptera frugiperda multiple nucleopolyhedrovirus from Nigeria, West Africa. The genome is 132,710 bp long and contains 144 open reading frames. The GC content is 40.3% and, based on baculovirus species demarcation criteria, the isolate belongs to the species Spodoptera frugiperda multiple nucleopolyhedrovirus.

Patterns in Genotype Composition of Indian Isolates of the Bombyx mori Nucleopolyhedrovirus and Bombyx mori Bidensovirus.

The mulberry silkworm, Bombyx mori (L.), is a model organism of lepidopteran insects with high economic importance. The viral diseases of the silkworm caused by Bombyx mori nucleopolyhedrovirus (BmNPV) and Bombyx mori bidensovirus (BmBDV) inflict huge economic losses and significantly impact the sericulture industry of India and other countries. To understand the distribution of Indian isolates of the BmNPV and to investigate their genetic composition, an in-depth population structure analysis was conducted using comprehensive and newly developed genomic analysis methods. The seven new Indian BmNPV isolates from Anantapur, Dehradun, Ghumarwin, Jammu, Kashmir, Mysore and Salem grouped in the BmNPV clade, and are most closely related to Autographa californica multiple nucleopolyhedrovirus and Rachiplusia ou multiple nucleopolyhedrovirus on the basis of gene sequencing and phylogenetic analyses of the partial polh, lef-8 and lef-9 gene fragments. The whole genome sequencing of three Indian BmNPV isolates from Mysore (-My), Jammu (-Ja) and Dehradun (-De) was conducted, and intra-isolate genetic variability was analyzed on the basis of variable SNP positions and the frequencies of alternative nucleotides. The results revealed that the BmNPV-De and BmNPV-Ja isolates are highly similar in their genotypic composition, whereas the population structure of BmNPV-My appeared rather pure and homogenous, with almost no or few genetic variations. The BmNPV-De and BmNPV-Ja samples further contained a significant amount of BmBDV belonging to the Bidnaviridae family. We elucidated the genotype composition within Indian BmNPV and BmBDV isolates, and the results presented have broad implications for our understanding of the genetic diversity and evolution of BmNPV and co-occurring BmBDV isolates.

Monitoring Insect Transposable Elements in Large Double-Stranded DNA Viruses Reveals Host-to-Virus and Virus-to-Virus Transposition.

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

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.

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

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

Molecular and morphological characterisation of a novel microsporidian species, Tubulinosema suzukii, infecting Drosophila suzukii (Diptera: Drosophilidae).

A microsporidium showing morphological characteristics typical of a Tubulinosema species was discovered in Drosophila suzukii. All developmental stages were diplokaryotic and grew in direct contact with the host cell cytoplasm. Spores from fresh preparations were ovoid to slightly pyriform and measured 4.29 × 2.47 µm in wet mount preparations. The spore wall consisted of a 125 nm thick endospore covered by a double layered exospore of 39 nm and 18 nm. The polar filament measured 67 µm in length, was slightly anisofilar and was arranged in ten coils in one or rarely two rows. The two posterior coils were 95 nm in diameter while the anterior coils were 115 nm in diameter. Early developmental stages were surrounded by electron-dense, 35.3 nm diameter, surface ornaments scattered over the membrane. Tubular elements with diameters of approximately 75 nm were seen attaching to the periphery of meronts and sporonts. Tissues infected included fat body, midgut and muscle. A 1915 bp rDNA fragment, covering the small subunit (SSU), the internal transcribed spacer (ITS) and the 5' end of the large subunit ribosomal DNA, was amplified by PCR and sequenced. Phylogenetic analyses of the SSU rDNA fragment revealed closest relationship to Tubulinosema pampeana (Host: Bombus atratus, South America) and Tubulinosema loxostegi (Host: Loxostege sticticalis, ubiquitous), but using the complete dataset of SSU-ITS-LSU rDNA genes revealed T. hippodamiae (Host: Hippodamiae convergens) as the most closely related species. Based on the morphological and genetic features a new species, Tubulinosema suzukii sp. nov., is proposed for this microsporidium isolated from D. suzukii.

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

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

Novel Diversity and Virulence Patterns Found in New Isolates of Cydia pomonella Granulovirus from China.

Cydia pomonella granulovirus (CpGV) is successfully used worldwide as a biocontrol agent of the codling moth (CM) (Cydia pomonella). The occurrence of CM populations with different modes of resistance against commercial CpGV preparations in Europe, as well as the invasiveness of CM in China, threatening major apple production areas there, requires the development of new control options. Utilizing the naturally occurring genetic diversity of CpGV can improve such control strategies. Here, we report the identification of seven new CpGV isolates that were collected from infected CM larvae in northwest China. Resistance testing using a discriminating CpGV concentration and the determination of the median lethal concentration (LC50) were performed to characterize their levels of virulence against susceptible and resistant CM larvae. The isolates were further screened for the presence of the 2 × 12-bp-repeat insertion in CpGV gene pe38 (open reading frame 24 [ORF24]), which was shown to be the target of type I resistance. It was found that three isolates, CpGV-JQ, -KS1, and -ZY2, could break type I resistance, although delayed mortality was observed in the infection process. All isolates followed the pe38 model of breaking type I resistance, except for CpGV-WW, which harbored the genetic factor but failed to overcome type I resistance. However, CpGV-WW was able to overcome type II and type III resistance. The bioassay results and sequencing data of pe38 support previous findings that pe38 is the major target for type I resistance. The new isolates show some distinct virulence characteristics when infection of different CM strains is considered.IMPORTANCE CpGV is a highly virulent pathogen of the codling moth (CM). It is registered and widely applied as a biocontrol agent in nearly all apple-growing countries worldwide. The emergence of CpGV resistance and the increasing lack of chemical control options require improvements to current control strategies. Natural CpGV isolates, as well as resistance-breaking isolates selected in resistant CM strains, have provided resources for improved resistance-breaking CpGV products. Here, we report novel CpGV isolates collected in China, which have new resistance-breaking capacities and may be an important asset for future application in the biological control of codling moths.

Single nucleotide polymorphism (SNP) frequencies and distribution reveal complex genetic composition of seven novel natural isolates of Cydia pomonella granulovirus.

The co-evolution between baculoviruses and their insect hosts results in selection of virus populations. To explore this phenomenon at the molecular level, seven natural isolates of Cydia pomonella granulovirus (CpGV) collected from orchards in northwest China were studied using Illumina next generation sequencing (NGS). A total of 540 genome positions with single nucleotide polymorphisms (SNPs) were detected in comparison with known CpGV isolates. New members of previously defined phylogenetic genome groups A, D and E of CpGV, as well as two novel phylogenetic lines, termed genome group F and G, were identified. Combining SNP frequency distribution with the prevalence of genome group-specific SNPs, revealed that six isolates of CpGV were mixtures of different ratios of at least two genotypes, whereas only one isolate, CpGV-WW, was genetically highly homogeneous. This study significantly extends our current understanding of the genetic diversity of CpGV and opens new lines of application of this virus.