Phylogenetic and recombination analyses of two deformed wing virus strains from different honeybee species in China

Deformed wing virus (DWV) is one of many viruses that infect honeybees and has been extensively studied because of its close association with honeybee colony collapse that is induced by Varroa destructor. However, virus genotypes, sequence characteristics, and genetic variations of DWV remain unknown in China.

Two DWV strains were isolated from Jinzhou and Qinhuangdao cities in China, and were named China1-2017 (accession number: MF770715) and China2-2018 (accession number: MH165180), respectively, and their complete genome sequences were analyzed. To investigate the phylogenetic relationships of the DWV isolates, a phylogenetic tree of the complete open reading frame (ORF), structural protein VP1, and non-structural protein 3C+RdRp of the DWV sequences was constructed using the MEGA 5.0 software program. Then, the similarity and recombinant events of the DWV isolated strains were analyzed using recombination detection program (RDP4) software and genetic algorithm for recombination detection (GARD).

The complete genomic analysis showed that the genomes of the China1-2017 and China2-2018 DWV strains consisted of 10,141 base pairs (bp) and 10,105 bp, respectively, and contained a single, large ORF (China1-2017: 1,146-9,827 bp; China2-2018: 1,351-9,816 bp) that encoded 2,894 amino acids. The sequences were compared with 20 previously reported DWV sequences from different countries and with sequences of two closely related viruses, Kakugo virus (KV) and V. destructor virus-1. Multiple sequence comparisons revealed a nucleotide identity of 84.3-96.7%, and identity of 94.7-98.6% in amino acids between the two isolate strains and 20 reference strains. The two novel isolates showed 96.7% nucleotide identity and 98.1% amino acid identity. The phylogenetic analyses showed that the two isolates belonged to DWV Type A and were closely related to the KV-2001 strain from Japan. Based on the RDP4 and GARD analyses, the recombination of the China2-2018 strain was located at the 4,266-7,507 nt region, with Korea I-2012 as an infer unknown parent and China-2017 as a minor parent, which spanned the entire helicase ORF. To the best of our knowledge, this is the first study to the complete sequence of DWV isolated from Apis cerana and the possible DWV recombination events in China. Our findings are important for further research of the phylogenetic relationship of DWVs in China with DWV strains from other countries and also contribute to the understanding of virological properties of these complex DWV recombinants.