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

Evidence of circulating recombinants between deformed wing virus and Varroa destructor virus-1 in honey bee colonies in Türkiye

Deformed wing virus (DWV), which is an important honey bee virus transmitted by Varroa destructor (V. destructor), causes colony losses in honey bee colonies. This study aimed to investigate the prevalence and genetic diversity of DWV in honey bees in Türkiye and to determine the role of V. destructor in the transmission of the genetic variants of DWV. Honey bee samples were collected from 62 apiaries, by simple random sampling, during March 2022 and April 2023. The presence of V. destructor in collected bee samples was examined using a stereo microscope. Real-time RT-PCR was used for the detection of DWV-A and DWV-B (Varroa destructor virus-1 (VDV-1)) viruses. Genetic characterisation of the positive samples was conducted by sequencing polyprotein genomic region. Considering the V. destructor infestation rate of 3% as relevant, out of the 62 apiaries examined, 17 (27.4%) were positive. However, DWV-A and VDV-1 specific RNA was not detected in V. destructor samples. VDV-1 specific RNA was detected in 6.5% (4/62) of the apiaries, whereas DWV-A was not detected in the sampled apiaries. Phylogenetic analysis showed that isolates detected in this study were located in a separate cluster from previously characterised DWV-A and VDV-1 isolates. According to RDP4 and GARD analyses, DWV-VDV-1 recombination breakpoints were detected in field isolates. To the best our knowledge, this is the first report of the presence of VDV-1-DWV recombinants in Türkiye. Further studies are needed to determine the impact of VDV-1-DWV recombinants and their virological and antigenic properties.

Recombination analysis on the receptor switching event of MERS-CoV and its close relatives: implications for the emergence of MERS-CoV

BACKGROUND

PlMERS-CoV is a coronavirus known to cause severe disease in humans, taxonomically classified under the subgenus Merbecovirus. Recent findings showed that the close relatives of MERS-CoV infecting vespertillionid bats (family Vespertillionidae), named NeoCoV and PDF-2180, use their hosts' ACE2 as their entry receptor, unlike the DPP4 receptor usage of MERS-CoV. Previous research suggests that this difference in receptor usage between these related viruses is a result of recombination. However, the precise location of the recombination breakpoints and the details of the recombination event leading to the change of receptor usage remain unclear.

METHODS

We used maximum likelihood-based phylogenetics and genetic similarity comparisons to characterise the evolutionary history of all complete Merbecovirus genome sequences. Recombination events were detected by multiple computational methods implemented in the recombination detection program. To verify the influence of recombination, we inferred the phylogenetic relation of the merbecovirus genomes excluding recombinant segments and that of the viruses' receptor binding domains and examined the level of congruency between the phylogenies. Finally, the geographic distribution of the genomes was inspected to identify the possible location where the recombination event occurred.

RESULTS

Similarity plot analysis and the recombination-partitioned phylogenetic inference showed that MERS-CoV is highly similar to NeoCoV (and PDF-2180) across its whole genome except for the spike-encoding region. This is confirmed to be due to recombination by confidently detecting a recombination event between the proximal ancestor of MERS-CoV and a currently unsampled merbecovirus clade. Notably, the upstream recombination breakpoint was detected in the N-terminal domain and the downstream breakpoint at the S2 subunit of spike, indicating that the acquired recombined fragment includes the receptor-binding domain. A tanglegram comparison further confirmed that the receptor binding domain-encoding region of MERS-CoV was acquired via recombination. Geographic mapping analysis on sampling sites suggests the possibility that the recombination event occurred in Africa.

CONCLUSION

Together, our results suggest that recombination can lead to receptor switching of merbecoviruses during circulation in bats. These results are useful for future epidemiological assessments and surveillance to understand the spillover risk of bat coronaviruses to the human population.

Nationwide genomic surveillance reveals the prevalence and evolution of honeybee viruses in China

BACKGROUND

The economic and environmental value of honeybees has been severely challenged in recent years by the collapse of their colonies worldwide, often caused by outbreaks of infectious diseases. However, our understanding of the diversity, prevalence, and transmission of honeybee viruses is largely obscure due to a lack of large-scale and longitudinal genomic surveillance on a global scale.

RESULTS

We report the meta-transcriptomic sequencing of nearly 2000 samples of the two most important economic and widely maintained honeybee species, as well as an associated ectoparasite mite, collected across China during 2016-2019. We document the natural diversity and evolution of honeybee viruses in China, providing evidence that multiple viruses commonly co-circulate within individual bee colonies. We also expanded the genomic data for 12 important honeybee viruses and revealed novel genetic variants and lineages associated with China. We identified more than 23 novel viruses from the honeybee and mite viromes, with some exhibiting ongoing replication in their respective hosts. Together, these data provide additional support to the idea that mites are an important reservoir and spill-over host for honeybee viruses.

CONCLUSIONS

Our data show that honeybee viruses are more widespread, prevalent, and genetically diverse than previously realized. The information provided is important in mitigating viral infectious diseases in honeybees, in turn helping to maintain sustainable productive agriculture on a global scale. Video Abstract.

Comparative Genomics and Characterization of SARS-CoV-2 P.1 (Gamma) Variant of Concern From Amazonas, Brazil

BACKGROUND

P.1 lineage (Gamma) was first described in the State of Amazonas, northern Brazil, in the end of 2020, and has emerged as a very important variant of concern (VOC) of SARS-CoV-2 worldwide. P.1 has been linked to increased infectivity, higher mortality, and immune evasion, leading to reinfections and potentially reduced efficacy of vaccines and neutralizing antibodies.

METHODS

The samples of 276 patients from the State of Amazonas were sent to a central referral laboratory for sequencing by gold standard techniques, through Illumina MiSeq platform. Both global and regional phylogenetic analyses of the successfully sequenced genomes were conducted through maximum likelihood method. Multiple alignments were obtained including previously obtained unique human SARS-CoV-2 sequences. The evolutionary histories of spike and non-structural proteins from ORF1a of northern genomes were described and their molecular evolution was analyzed for detection of positive (FUBAR, FEL, and MEME) and negative (FEL and SLAC) selective pressures. To further evaluate the possible pathways of evolution leading to the emergence of P.1, we performed specific analysis for copy-choice recombination events. A global phylogenomic analysis with subsampled P.1 and B.1.1.28 genomes was applied to evaluate the relationship among samples.

RESULTS

Forty-four samples from the State of Amazonas were successfully sequenced and confirmed as P.1 (Gamma) lineage. In addition to previously described P.1 characteristic mutations, we find evidence of continuous diversification of SARS-CoV-2, as rare and previously unseen P.1 mutations were detected in spike and non-structural protein from ORF1a. No evidence of recombination was found. Several sites were demonstrated to be under positive and negative selection, with various mutations identified mostly in P.1 lineage. According to the Pango assignment, phylogenomic analyses indicate all samples as belonging to the P.1 lineage.

CONCLUSION

P.1 has shown continuous evolution after its emergence. The lack of clear evidence for recombination and the positive selection demonstrated for several sites suggest that this lineage emergence resulted mainly from strong evolutionary forces and progressive accumulation of a favorable signature set of mutations.

Deformed wing virus variant shift from 2010 to 2016 in managed and feral UK honey bee colonies

Deformed wing virus (DWV) has been linked to the global decline of honey bees. DWV exists as three master variants (DWV-A, DWV-B, and DWV-C), each with differing outcomes for the honey bee host. Research in the USA showed a shift from DWV-A to DWV-B between 2010 to 2016 in honey bee colonies. Likewise, in the UK, a small study in 2007 found only DWV-A, whereas in 2016, DWV-B was the most prevalent variant. This suggests a shift from DWV-A to DWV-B might have occurred in the UK between 2007 and 2016. To investigate this further, data from samples collected in 2009/10 (n = 46) were compared to existing data from 2016 (n = 42). These samples also allowed a comparison of DWV variants between Varroa-untreated (feral) and Varroa-treated (managed) colonies. The results revealed that, in the UK, DWV-A was far more prevalent in 2009/10 (87%) than in 2016 (43%). In contrast, DWV-B was less prevalent in 2009/10 (76%) than in 2016 (93%). Regardless if colonies had been treated for Varroa (managed) or not (feral), the same trend from DWV-A to DWV-B occurred. Overall, the results reveal a decrease in DWV-A and an increase in DWV-B in UK colonies.

Ten Years of Deformed Wing Virus (DWV) in Hawaiian Honey Bees (Apis mellifera), the Dominant DWV-A Variant Is Potentially Being Replaced by Variants with a DWV-B Coding Sequence

The combination of Deformed wing virus (DWV) and Varroa destructor is arguably one of the greatest threats currently facing western honey bees, Apis mellifera. Varroa's association with DWV has decreased viral diversity and increased loads of DWV within honey bee populations. Nowhere has this been better studied than in Hawaii, where the arrival of Varroa progressively led to the dominance of the single master variant (DWV-A) on both mite-infested Hawaiian Islands of Oahu and Big Island. Now, exactly 10 years following the original study, we find that the DWV population has changed once again, with variants containing the RdRp coding sequence pertaining to the master variant B beginning to co-dominate alongside variants with the DWV-A RdRp sequence on the mite-infested islands of Oahu and Big Island. In speculation, based on other studies, it appears this could represent a stage in the journey towards the complete dominance of DWV-B, a variant that appears better adapted to be transmitted within honey bee colonies.

Possible Spillover of Pathogens between Bee Communities Foraging on the Same Floral Resource

Simple Summary

Floral resource availability is one of the keys to preserving the health of bee communities. However, flowers also present a risk of pathogen transmission, as infected pollinators could deposit pathogens while foraging, exposing other pollinators to infection via the consumption of contaminated nectar or pollen. Here, we studied, over time, the prevalence of seven viruses in bee communities that share the same small surface of floral resource in order to assess the risk of virus spillover. In total, 2057 bee specimens from 30 species were caught, identified and checked for the presence of viruses. Specimens from the Halictidae family were the dominant wild bees. The prevalence of viruses was quite high: at least one virus was detected in 78% of the samples, and co-infections were frequent. The genetic diversity of the viruses was also investigated to look for the possible association of geographic origin or host with shared ancestry.

Abstract

Viruses are known to contribute to bee population decline. Possible spillover is suspected from the co-occurrence of viruses in wild bees and honey bees. In order to study the risk of virus transmission between wild and managed bee species sharing the same floral resource, we tried to maximize the possible cross-infections using Phacelia tanacetifolia, which is highly attractive to honey bees and a broad range of wild bee species. Virus prevalence was compared over two years in Southern France. A total of 1137 wild bees from 29 wild bee species (based on COI barcoding) and 920 honey bees (Apis mellifera) were checked for the seven most common honey bee RNA viruses. Halictid bees were the most abundant. Co-infections were frequent, and Sacbrood virus (SBV), Black queen cell virus (BQCV), Acute bee paralysis virus (ABPV) and Israeli acute paralysis virus (IAPV) were widespread in the hymenopteran pollinator community. Conversely, Deformed wing virus (DWV) was detected at low levels in wild bees, whereas it was highly prevalent in honey bees (78.3% of the samples). Both wild bee and honey bee virus isolates were sequenced to look for possible host-specificity or geographical structuring. ABPV phylogeny suggested a specific cluster for Eucera bees, while isolates of DWV from bumble bees (Bombus spp.) clustered together with honey bee isolates, suggesting a possible spillover.

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

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.

Comparison of coxsackievirus A12 genome isolated from patients with different symptoms in Yunnan, Southwest China

Aim: Coxsackievirus A12 (CVA12) mainly causes hand–foot–mouth disease (HFMD) or herpangina (HA). Thus, the genomic characteristics of CVA12 isolates from Southwest China, especially the discrepancy between CVA12-HFMD and CVA12-HA were analyzed. Patients & methods: The full-length genome sequences of CVA12-HFMD and CVA12-HA were obtained and phylogenetic, nucleotide mutation, amino acid substitution and recombinant analyses were performed. Results: All CVA12 were closest to the Netherlands (2013) and have possibly recombined in the capsid coding region (P1) with other HEV-A. In the coding region, 54 base mutation result in 11 nonsynonymous mutations and four of them were identical mutations between CVA12-HFMD and CVA12-HA. Conclusion: Whether the presence of four consistent nonsynonymous mutation sites affect the virulence of the CVA12 deserves further study.

The honeybee (Apis mellifera) developmental state shapes the genetic composition of the deformed wing virus-A quasispecies during serial transmission

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it.

Phylogenetic Analysis of Belgian Small Ruminant Lentiviruses Supports Cross Species Virus Transmission and Identifies New Subtype B5 Strains

Small ruminant lentiviruses (SRLV) are a group of highly divergent viruses responsible for global and fatal infections in sheep and goats. Since the current phylogenetic classification of these viruses was proposed in 2004, it nowadays consists out of 5 genotypes and 28 subtypes. In support of our national SRLV control program, we performed the genetic characterization of SRLV strains circulating in the Belgian sheep and goat population. Fourteen sheep and 9 goat strains were sequenced in the gag-pol and pol regions using the method described by Shah. Most SRLV strains from sheep and goats belonged to prototype A1 and B1 subtypes, respectively. We, however, also found indications for cross-species transmission of SRLV strains between sheep and goats and vice versa, and identified a new subtype designated as B5. An in-depth analysis of the current SRLV phylogeny revealed that many subtypes have been defined over the years based on limited sequence information. To keep phylogeny as a useful tool, we advocate to apply more rigorous sequencing standards to ensure the correct classification of current and new emerging strains. The genetic characterization of Belgian SRLV strains will help in the development of appropriate diagnostic tools to assist the national control program.

Dynamic evolution in the key honey bee pathogen deformed wing virus: Novel insights into virulence and competition using reverse genetics

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population.

Deformed wing virus, a key pathogen of honey bees, shows rapid diversification after genetic bottlenecks; a novel reverse-genetic system provides insights into the forces that shape virus diversity, suggesting that virus quasi-species diversification may be driven by selection of genotypes capable of evading host RNAi defences.