Genome analysis and phylogenetic characterization of two deformed wing virus strains from Apis cerana in Vietnam

Thu, Ha T.; Lien, Nguyen T.K.; Lanh, Pham T.; Duong, Bui T.T.; Hoa, Nguyen T.; Phuoc, Man H.; Thai, Pham H.; Quyen, Dong Van

BACKGROUND: Deformed wing virus (DWV) is a virulent virus that causes honeybee disease. DWV can exist as a latent infection in honeybees, outbreak into epidemics, and cause serious damage to beekeeping cross the world, including Vietnam. METHODS: The two DWV strains circulating in Vietnamese honeybee, Apis cerana, were first isolated from adult honeybees in North Vietnam (DWV-NVN) and South Vietnam (DWV-SVN). Their complete nucleotide sequences were determined, aligned, and compared with other DWV strains. RESULTS: The two Vietnamese DWV strains comprised 10,113 bp and contained a large single open reading frame (ORF) of 2,893 amino acids, initiating at nucleotide 1,130 and terminating at nucleotide 9,812. Multiple nucleotide sequence alignment between these two DWV-VN strains and DWV strains in A. mellifera was performed. The DWV-VN strains showed a low genetic identity (from 91.4% to 92.0%) with almost of these strains, but lower identities (89.2% and 89.4%) with UK2 and (89.6%) with the China2 strain. Low identities (91.7% and 91.9%) were also observed between the China3 strain (in A. cerana) and the DWV-VN strains, respectively. The deduced amino acid sequence alignment showed high genetic similarities (97.0%–97.9%) when the USA1, Chile, Italy1, France, UK1, UK2, Japan, Korea2, China1, China2 and China3 strains were compared to the DWV-VN strains. This ratio was 96.7% and 96.8% when the Korea1 strain was compared to the DWV-SVN and DWV-NVN strains, respectively. Numerous amino acid substitutions were identified in the L, VP3, and RdRp sequences. Notably, we observed six substitutions positioned at amino acids 27 (E > I), 98 (S > T), 120 (A > V), 153 (M > T), 170 (D > F), and 174 (Y > F) in the L protein, two amino acid changes at positions 980 (S > A) and 1032 (E > T) in VP3, and one amino acid change at position 2627 (R > C) unique to the DWV-VN strains. Phylogenetic analysis based on complete genome sequences, RdRp sequences and Simplot analysis indicated that there was a significant difference between DWV-VN strains in A. cerana and DWV strains in A. mellifera. The results suggested that the genetic variations of the DWV-VN strains in A. cerana help them to adapt geographical conditions and may lead to change the viral pathogenicity of DWV-VN strains.