Epidemiology of a major honey bee pathogen, deformed wing virus: potential worldwide replacement of genotype A by genotype B

The western honey bee (Apis mellifera) is of major economic and ecological importance, with elevated rates of colony losses in temperate regions over the last two decades thought to be largely caused by the exotic ectoparasitic mite Varroa destructor and deformed wing virus (DWV), which the mite transmits. DWV currently exists as two main genotypes: the formerly widespread DWV-A and the more recently described and rapidly expanding DWV-B. It is an excellent system to understand viral evolution and the replacement of one viral variant by another. Here we synthesise published results on the distribution and prevalence of DWV-A and -B over the period 2008-2021 and present novel data for Germany, Italy and the UK to suggest that (i) DWV-B has rapidly expanded worldwide since its first description in 2004 and (ii) that it is potentially replacing DWV-A. Both genotypes are also found in wild bee species. Based on a simple mathematical model, we suggest that interference between viral genotypes when co-infecting the same host is key to understanding their epidemiology. We finally discuss the consequences of genotype replacement for beekeeping and for wild pollinator species.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids

Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.

A56 Visualization of recombination in deformed wing virus infecting bees

Honey bees suffer increasing colony mortality worldwide, partially caused by the spread of viral pathogens. Among these pathogens, deformed wing virus (DWV) is one of the major, widespread viruses of honey bees resulting in wing deformities and weakening colonies. DWV can be found in honey bees, bumble bees, and other wild bees as three major genotypes named DWV-A, -B (also named Varroa destructor virus 1), and -C. Various recombinants of DWV-A and -B have been previously found in honey bees, some of which have been suggested to have higher virulence over non-recombinant, parental virus. In most of these cases, recombinants were only shown as consensus sequences from previous assemblies and alignments and may not reflect the biological reality of all variants present within a host bee. It is therefore important to build a method of recombinant detection and quantification within mixed infections in single-host individuals, including both parental and various recombinant genomes, so as to evaluate the relevance of recombinants for viral genome evolution and the impact on hosts. Here, we propose to visualize and quantify these recombinants using next-generation sequencing data to better understand how these genomes evolve within bees. Our method will be performed directly from raw sequence reads from various datasets (including field and lab experiments as well as screening of public databases) in order to obtain an overview of DWV recombination in various in vivo and in vitro conditions. Recombination of viral genomes is a key point for virus evolution. The detection and quantification of recombination will facilitate analysis of the determinants of recombination and help in understanding the routes by which new viral variants emerge. The emergence of new (more virulent) recombinant viruses can result from acquisition of new capabilities, such as escape from host immunity or increased transmission rates. Recombination can also lead to adaptation to new environments and new hosts by a change in cell tropism, allowing cross-species transmission, which may be particularly relevant for bumble bees and wild bees infected by honey bee-derived DWV.

The Two Prevalent Genotypes of an Emerging Infectious Disease, Deformed Wing Virus, Cause Equally Low Pupal Mortality and Equally High Wing Deformities in Host Honey Bees

Deformed wing virus (DWV) is an emerging infectious disease of the honey bee (Apis mellifera) that is considered a major cause of elevated losses of honey bee colonies. DWV comprises two widespread genotypes: the originally described genotype A, and genotype B. In adult honey bees, DWV-B has been shown to be more virulent than DWV-A. However, their comparative effects on earlier host developmental stages are unknown. Here, we experimentally inoculated honey bee pupae and tested for the relative impact of DWV-A versus DWV-B on mortality and wing deformities in eclosing adults. DWV-A and DWV-B caused similar, and only slightly elevated, pupal mortality (mean 18% greater mortality than control). Both genotypes caused similarly high wing deformities in eclosing adults (mean 60% greater wing deformities than control). Viral titer was high in all of the experimentally inoculated eclosing adults, and was independent of wing deformities, suggesting that the phenotype 'deformed wings' is not directly related to viral titer or viral genotype. These viral traits favor the emergence of both genotypes of DWV by not limiting the reproduction of its vector, the ectoparasitic Varroa destructor mite, in infected pupae, and thereby facilitating the spread of DWV in honey bees infested by the mite.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Discovery of parvovirus-related sequences in an unexpected broad range of animals

Our knowledge of the genetic diversity and host ranges of viruses is fragmentary. This is particularly true for the Parvoviridae family. Genetic diversity studies of single stranded DNA viruses within this family have been largely focused on arthropod- and vertebrate-infecting species that cause diseases of humans and our domesticated animals: a focus that has biased our perception of parvovirus diversity. While metagenomics approaches could help rectify this bias, so too could transcriptomics studies. Large amounts of transcriptomic data are available for a diverse array of animal species and whenever this data has inadvertently been gathered from virus-infected individuals, it could contain detectable viral transcripts. We therefore performed a systematic search for parvovirus-related sequences (PRSs) within publicly available transcript, genome and protein databases and eleven new transcriptome datasets. This revealed 463 PRSs in the transcript databases of 118 animals. At least 41 of these PRSs are likely integrated within animal genomes in that they were also found within genomic sequence databases. Besides illuminating the ubiquity of parvoviruses, the number of parvoviral sequences discovered within public databases revealed numerous previously unknown parvovirus-host combinations; particularly in invertebrates. Our findings suggest that the host-ranges of extant parvoviruses might span the entire animal kingdom.

Comparative Analysis of Transcriptomes from Secondary Reproductives of Three Reticulitermes Termite Species

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.

The effects of quisqualate and nocodazole on the organization of MAP2 and neurofilaments in spinal cord neurons in vitro

The relationship between microtubules, neurofilaments and microtubule-associated protein (MAP)2 was investigated in spinal cord neurons grown for up to 14 days in vitro. Neurons were labelled using antibodies against MAP2, neurofilaments and tubulin, and immunofluorescence analyzed by confocal microscopy. A well-structured network of neurofilaments and microtubules was observed in unstimulated cultures. MAP2 staining was poorly structured but became more filamentous following depolymerization of microtubules with nocodazole. Double-staining experiments suggested that MAP2 was now closely associated with neurofilaments in cell bodies and dendrites. Stimulation of cultures with excitatory amino acids increased the resistance of the microtubular cytoskeleton to depolymerization by nocodazole. Again double-labelling experiments demonstrated an increased association between neurofilaments and MAP2 immunofluorescence. Previous results suggested that the stability of the neuronal cytoskeleton could be modulated by glutamate receptors acting through an increased binding of MAP2 to microtubules. From the results presented here, we further suggest that cross-linking of neurofilaments to microtubules may also play a role in this process.

Reorganization of the Cytoskeleton in Rat Neurons Following Stimulation With Excitatory Amino Acids In Vitro

The state of neuronal microtubule polymerization is influenced by microtubule-associated proteins such as MAP2, which is specifically localized within neuronal dendrites and cell bodies. We have demonstrated that stimulation of spinal cord or cortical neurons in vitro with excitatory amino acids results in a dramatic modification of the neuronal cytoskeleton as monitored with antibodies against MAP2 and tubulin. Stimulation of cultures with glutamate receptor agonists induced a reorganization of MAP2 immunoreactivity into a distinctive network of bundles within certain neuronal cell bodies and their proximal neurites. The effect was not abolished by depolymerizing drugs such as nocodazole, or protein synthesis inhibitors. The effect was dependent upon the entry of sodium following depolarization and was not associated with neuronal damage. We suggest that in neurons the state of the neuronal cytoskeleton can be modulated by glutamate receptor activation acting through MAP2.

Effect of excitatory amino acids on microtubule-associated proteins in cultured cortical and spinal neurones

The effect of excitatory amino acid stimulation on the cytoskeleton of cultured spinal cord and cortical neurons was monitored with antibodies against microtubule-associated proteins tau and MAP2. In unstimulated cultures tau-1 immunoreactivity was restricted to axon-like processes. Stimulation with glutamate (0.1-1 mM) or N-methyl-D-aspartate (NMDA) (0.1 mM) resulted in a dramatic increase in the intensity of tau labelling in axons and the appearance of staining within a proportion of neuronal cell bodies and dendrites. Quisqualate or kainate stimulation resulted only in an increase in tau immunoreactivity within axons. The NMDA mediated events were calcium dependent and the effects of all excitatory amino acids could be blocked by specific antagonists. In contrast, following stimulation with excitatory amino acids, MAP2-immunoreactivity was associated with filaments which formed a complex network within the cell body. This suggests that the different excitatory amino acid receptor subtypes can have differential effects on the neuronal cytoskeleton.