Genomic survey of the ectoparasitic mite Varroa destructor, a major pest of the honey bee Apis mellifera

Virus replication in the honey bee parasite, Varroa destructor

IMPORTANCE

The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

The Impact of Whole Genome Duplication on the Evolution of the Arachnids

The proliferation of genomic resources for Chelicerata in the past 10 years has revealed that the evolution of chelicerate genomes is more dynamic than previously thought, with multiple waves of ancient whole genome duplications affecting separate lineages. Such duplication events are fascinating from the perspective of evolutionary history because the burst of new gene copies associated with genome duplications facilitates the acquisition of new gene functions (neofunctionalization), which may in turn lead to morphological novelties and spur net diversification. While neofunctionalization has been invoked in several contexts with respect to the success and diversity of spiders, the overall impact of whole genome duplications on chelicerate evolution and development remains imperfectly understood. The purpose of this review is to examine critically the role of whole genome duplication on the diversification of the extant arachnid orders, as well as assess functional datasets for evidence of subfunctionalization or neofunctionalization in chelicerates. This examination focuses on functional data from two focal model taxa: the spider Parasteatoda tepidariorum, which exhibits evidence for an ancient duplication, and the harvestman Phalangium opilio, which exhibits an unduplicated genome. I show that there is no evidence that taxa with genome duplications are more successful than taxa with unduplicated genomes. I contend that evidence for sub- or neofunctionalization of duplicated developmental patterning genes in spiders is indirect or fragmentary at present, despite the appeal of this postulate for explaining the success of groups like spiders. Available expression data suggest that the condition of duplicated Hox modules may have played a role in promoting body plan disparity in the posterior tagma of some orders, such as spiders and scorpions, but functional data substantiating this postulate are critically missing. Spatiotemporal dynamics of duplicated transcription factors in spiders may represent cases of developmental system drift, rather than neofunctionalization. Developmental system drift may represent an important, but overlooked, null hypothesis for studies of paralogs in chelicerate developmental biology. To distinguish between subfunctionalization, neofunctionalization, and developmental system drift, concomitant establishment of comparative functional datasets from taxa exhibiting the genome duplication, as well as those that lack the paralogy, is sorely needed.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

First detection of Lake Sinai virus in the Czech Republic: a potential member of a new species

Lake Sinai virus (LSV) is one of over 20 honey bee viruses. Variants of LSV have been classified as members of two officially recognised species, Lake Sinai virus 1 and Lake Sinai virus 2. However, there are currently a limited number of whole-genome sequences, and the genetic variability of the virus indicates that additional species may need to be established. Extracted nucleic acid of 209 honey bee samples was screened by PCR for 11 honey bee viruses. LSV was the third most abundant virus (36.9% of positive samples), after Apis mellifera filamentous virus (72.2%) and deformed wing virus (52.5%). LSV-positive samples were analyzed further by PCR with primers targeting the region encoding the viral RNA-dependent RNA polymerase. Subsequently, the PCR products were sequenced, and the resulting sequences were used for a first round of phylogenetic analysis. Based on those results, several isolates were selected for whole-genome sequencing, and the complete genome sequences were used for additional phylogenetic analysis. The results indicated the presence of at least three genetically distinct groups of LSV in the Czech Republic, the most prevalent one being related to LSV 2 but too dissimilar to be considered a member of the same species. Two sequences of a major LSV strain cluster native to the Czech Republic were determined, representing the first Czech LSV strains published to date.

Viral communities in the parasite Varroa destructor and in colonies of their honey bee host (Apis mellifera) in New Zealand

The parasitic mite Varroa destructor is a leading cause of mortality for Western honey bee (Apis mellifera) colonies around the globe. We sought to confirm the presence and likely introduction of only one V. destructor haplotype in New Zealand, and describe the viral community within both V. destructor mites and the bees that they parasitise. A 1232 bp fragment from mitochondrial gene regions suggests the likely introduction of only one V. destructor haplotype to New Zealand. Seventeen viruses were found in bees. The most prevalent and abundant was the Deformed wing virus A (DWV-A) strain, which explained 95.0% of the variation in the viral community of bees. Black queen cell virus, Sacbrood virus, and Varroa destructor virus 2 (VDV-2) played secondary roles. DWV-B and the Israeli acute paralysis virus appeared absent from New Zealand. Ten viruses were observed in V. destructor, with > 99.9% of viral reads from DWV-A and VDV-2. Substantially more variation in viral loads was observed in bees compared to mites. Where high levels of VDV-2 occurred in mites, reduced DWV-A occurred in both the mites and the bees co-occurring within the same hive. Where there were high loads of DWV-A in mites, there were typically high viral loads in bees.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Genomics and Proteomics of Apis mellifera Filamentous Virus Isolated from Honeybees in China

Apis mellifera filamentous virus (AmFV) is a large DNA virus that is endemic in honeybee colonies. The genome sequence of the AmFV Swiss isolate (AmFV CH–C05) has been reported, but so far very few molecular studies have been conducted on this virus. In this study, we isolated and purified AmFV (AmFV CN) from Chinese honeybee (Apis mellifera) colonies and elucidated its genomics and proteomics. Electron microscopy showed ovoid purified virions with dimensions of 300–500 ​× ​210–285 ​nm, wrapping a 3165 ​× ​40 ​nm filamentous nucleocapsid in three figure-eight loops. Unlike AmFV CH–C05, which was reported to have a circular genome, our data suggest that AmFV CN has a linear genome of approximately 493 ​kb. A total of 197 ORFs were identified, among which 36 putative genes including 18 baculoviral homologs were annotated. The overall nucleotide similarity between the CN and CH–C05 isolates was 96.9%. Several ORFs were newly annotated in AmFV CN, including homologs of per os infectivity factor 4 (PIF4) and a putative integrase. Phylogenomic analysis placed AmFVs on a separate branch within the newly proposed virus class Naldaviricetes. Proteomic analysis revealed 47 AmFV virion-associated proteins, of which 14 had over 50% sequence coverage, suggesting that they are likely to be main structural proteins. In addition, all six of the annotated PIFs (PIF-0–5) were identified by proteomics, suggesting that they may function as entry factors in AmFV infection. This study provides fundamental information regarding the molecular biology of AmFV.

•

The AmFV CN contains a 493 ​kb linear genome encoding 197 ORFs.

•

Proteomics revealed 14 putative major structural proteins.

•

AmFV belongs to the class Naldaviricetes but not the order Lefavirales.

De novo transcriptome sequencing of the northern fowl mite, Ornithonyssus sylviarum, shed light on parasitiform poultry mites evolution and its chemoreceptor repertoires

The northern fowl mite (NFM), Ornithonyssus sylviarum, and the poultry red mite (PRM), Dermanyssus gallinae, are the most serious pests of poultry, both of which have an expanding global prevalence. Research on NFM has been constrained by a lack of genomic and transcriptomic data. Here, we report and analyze the first global transcriptome data across all mite live stages and sexes. A total of 28,999 unigenes were assembled, of which 19,750 (68.10%) were annotated using seven functional databases. The biological function of these unigenes was classified using the GO, KOG, and KEGG databases. To gain insight into the chemosensory receptor-based system of parasitiform mites, we furthermore assessed the gene repertoire of gustatory receptors (GRs) and ionotropic receptors (IRs), both of which encode putative ligand-gated ion channel proteins. While these receptors are well characterized in insect model species, our understanding of chemosensory detection in mites and ticks is in its infancy. To address this paucity of data, we identified 9 IR/iGluRs and 2 GRs genes by analyzing transcriptome data in the NFM, while 9 GRs and 41 IR/iGluRs genes were annotated in the PRM genome. Taken together, the transcriptomic and genomic annotation of these two species provide a valuable reference for studies of parasitiform mites and also help to understand how chemosensory gene family expansion/contraction events may have been reshaped by an obligate parasitic lifestyle compared with their free-living closest relatives. Future studies should include additional species to validate this observation and functional characterization of the identified proteins as a step forward in identifying tools for controlling these poultry pests.

Integrated Pest Management Control of Varroa destructor (Acari: Varroidae), the Most Damaging Pest of (Apis mellifera L. (Hymenoptera: Apidae)) Colonies

Varroa destructor is among the greatest biological threats to western honey bee (Apis mellifera L.) health worldwide. Beekeepers routinely use chemical treatments to control this parasite, though overuse and mismanagement of these treatments have led to widespread resistance in Varroa populations. Integrated Pest Management (IPM) is an ecologically based, sustainable approach to pest management that relies on a combination of control tactics that minimize environmental impacts. Herein, we provide an in-depth review of the components of IPM in a Varroa control context. These include determining economic thresholds for the mite, identification of and monitoring for Varroa, prevention strategies, and risk conscious treatments. Furthermore, we provide a detailed review of cultural, mechanical, biological, and chemical control strategies, both longstanding and emerging, used against Varroa globally. For each control type, we describe all available treatments, their efficacies against Varroa as described in the primary scientific literature, and the obstacles to their adoption. Unfortunately, reliable IPM protocols do not exist for Varroa due to the complex biology of the mite and strong reliance on chemical control by beekeepers. To encourage beekeeper adoption, a successful IPM approach to Varroa control in managed colonies must be an improvement over conventional control methods and include cost-effective treatments that can be employed readily by beekeepers. It is our intention to provide the most thorough review of Varroa control options available, ultimately framing our discussion within the context of IPM. We hope this article is a call-to-arms against the most damaging pest managed honey bee colonies face worldwide.

The Odorant-Binding Proteins of the Spider Mite Tetranychus urticae

Spider mites are one of the major agricultural pests, feeding on a large variety of plants. As a contribution to understanding chemical communication in these arthropods, we have characterized a recently discovered class of odorant-binding proteins (OBPs) in Tetranychus urticae. As in other species of Chelicerata, the four OBPs of T. urticae contain six conserved cysteines paired in a pattern (C1-C6, C2-C3, C4-C5) differing from that of insect counterparts (C1-C3, C2-C5, C4-C6). Proteomic analysis uncovered a second family of OBPs, including twelve members that are likely to be unique to T. urticae. A three-dimensional model of TurtOBP1, built on the recent X-ray structure of Varroa destructor OBP1, shows protein folding different from that of insect OBPs, although with some common features. Ligand-binding experiments indicated some affinity to coniferyl aldehyde, but specific ligands may still need to be found among very large molecules, as suggested by the size of the binding pocket.

Adaptive population structure shifts in invasive parasitic mites, Varroa destructor

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors governing host-parasite interactions. Even though invasive parasites are considered of major biological importance, little is known about their adaptative potential when infesting the new hosts. Here, the genetic diversification of Varroa destructor, a novel parasite of Apis mellifera originating from Asia, was investigated using population genetics to determine how the genetic structure of the parasite changed in distinct European populations of its new host. To do so, mites infesting two categories of hosts in four European regions were compared: (a) adapted hosts surviving through means of natural selection, thereby expected to impose strong selective pressure on the mites, and (b) treated host populations, surviving mite infestations because acaricides are applied, therefore characterized by a relaxed selection imposed by the host on the mites. Significant genetic divergence was found across regions, partially reflecting the invasion pattern of V. destructor throughout Europe and indicating local adaptation of the mite to the host populations. Additionally, varying degrees of genotypic changes were found between mites from adapted and treated colonies. Altogether, these results indicate that V. destructor managed to overcome the genetic bottlenecks following its introduction in Europe and that host-mediated selection fostered changes in the genetic structure of this mite at diverse geographic scales. These findings highlight the potential of parasites to adapt to their local host populations and confirm that adaptations developed within coevolutionary dynamics are a major determinant of population genetic changes.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging bees.

A genome-wide screening for RNAi pathway proteins in Acari

BACKGROUND

RNA interference (RNAi) is a highly conserved, sequence-specific gene silencing mechanism present in Eukaryotes. Three RNAi pathways are known, namely micro-RNA (miRNA), piwi-interacting RNA (piRNA) and short interfering RNA (siRNA). However, little knowledge exists about the proteins involved in these pathways in Acari. Moreover, variable successes has been obtained in gene knockdown via siRNA pathway in their functional genomics and management. We hypothesized that the clue may be in the variability of the composition and the efficacy of siRNA machinery among Acari.

RESULTS

Both comparative genomic analyses and domain annotation suggest that all the analyzed species have homologs of putative core proteins that mediate cleaving of targeted genes via the three RNAi pathways. We identified putative homologs of Caenorhabditis elegans RNA-dependent RNA polymerase (RdRP) protein in all species though no secondary Argonaute homologs that operate with this protein in siRNA amplification mechanism were found, suggesting that the siRNA amplification mechanism present in Acari may be distinct from that described in C. elegans. Moreover, the genomes of these species do not encode homologs of C. elegans systemic RNAi defective-1 (Sid-1) protein that mediate silencing of the mRNA target throughout the treated organisms suggesting that the phenomena of systemic RNAi that has been reported in some Acari species probably occur through a different mechanism. However, homologs of putative RNAi spreading defective-3 (Rsd-3) protein and scavenger receptors namely Eater and SR-CI that mediate endocytosis cellular update of dsRNA in C. elegans and Drosophila melanogaster were found in Acari genomes. This result suggests that cellular dsRNA uptake in Acari is endocytosis-dependent. Detailed phylogenetic analyses of core RNAi pathway proteins in the studied species revealed that their evolution is compatible with the proposed monophyletic evolution of this group.

CONCLUSIONS

Our analyses have revealed the potential activity of all three pathways in Acari. Still, much experimental work remains to be done to confirm the mechanisms behind these pathways in particular those that govern systemic/parental RNAi and siRNA amplification in Acari. Disclosure of these mechanisms will facilitate the development of new and specific management tools for the harmful species and enrichment of the beneficial species.

Genome-Wide Identification of Long Non-coding RNAs in the Gravid Ectoparasite Varroa destructor

Long non-coding RNAs (lncRNAs) emerge as critical regulators across a wide variety of biological functions in living organisms. However, to date, no systematic characterization of lncRNAs has been investigated in the ectoparasitic mite Varroa destructor, the most severe biotic threat to honey bees worldwide. Here, we performed an initial genome-wide identification of lncRNAs in V. destructor via high-throughput sequencing technology and reported, for the first time, the transcriptomic landscape of lncRNAs in the devastating parasite. By means of a lncRNA identification pipeline, 6,645 novel lncRNA transcripts, encoded by 3,897 gene loci, were identified, including 2,066 sense lncRNAs, 2,772 lincRNAs, and 1,807 lncNATs. Compared with protein-coding mRNAs, V. destructor lncRNAs are shorter in terms of full length, as well as of the ORF length, contain less exons, and express at lower level. GO term and KEGG pathway enrichment analyses of the lncRNA target genes demonstrated that these predicted lncRNAs may be potentially responsible for the regulatory functions of cellular and biological progresses in the reproductive phase of V. destructor. To our knowledge, this is the first catalog of lncRNA profile in the parasitiformes species, providing a valuable resource for genetic and genomic studies. Understanding the characteristics and features of lncRNAs in V. destructor would promote sustainable parasite control.

Diversity and Global Distribution of Viruses of the Western Honey Bee, Apis mellifera

In the past centuries, viruses have benefited from globalization to spread across the globe, infecting new host species and populations. A growing number of viruses have been documented in the western honey bee, Apis mellifera. Several of these contribute significantly to honey bee colony losses. This review synthetizes the knowledge of the diversity and distribution of honey-bee-infecting viruses, including recent data from high-throughput sequencing (HTS). After presenting the diversity of viruses and their corresponding symptoms, we surveyed the scientific literature for the prevalence of these pathogens across the globe. The geographical distribution shows that the most prevalent viruses (deformed wing virus, sacbrood virus, black queen cell virus and acute paralysis complex) are also the most widely distributed. We discuss the ecological drivers that influence the distribution of these pathogens in worldwide honey bee populations. Besides the natural transmission routes and the resulting temporal dynamics, global trade contributes to their dissemination. As recent evidence shows that these viruses are often multihost pathogens, their spread is a risk for both the beekeeping industry and the pollination services provided by managed and wild pollinators.

Genomic insights into mite phylogeny, fitness, development, and reproduction

BACKGROUND

Predatory mites (Acari: Phytoseiidae) are the most important beneficial arthropods used in augmentative biological pest control of protected crops around the world. However, the genomes of mites are far less well understood than those of insects and the evolutionary relationships among mite and other chelicerate orders are contested, with the enigmatic origin of mites at one of the centres in discussion of the evolution of Arachnida.

RESULTS

We here report the 173 Mb nuclear genome (from 51.75 Gb pairs of Illumina reads) of the predatory mite, Neoseiulus cucumeris, a biocontrol agent against pests such as mites and thrips worldwide. We identified nearly 20.6 Mb (~ 11.93% of this genome) of repetitive sequences and annotated 18,735 protein-coding genes (a typical gene 2888 bp in size); the total length of protein-coding genes was about 50.55 Mb (29.2% of this assembly). About 37% (6981) of the genes are unique to N. cucumeris based on comparison with other arachnid genomes. Our phylogenomic analysis supported the monophyly of Acari, therefore rejecting the biphyletic origin of mites advocated by other studies based on limited gene fragments or few taxa in recent years. Our transcriptomic analyses of different life stages of N. cucumeris provide new insights into genes involved in its development. Putative genes involved in vitellogenesis, regulation of oviposition, sex determination, development of legs, signal perception, detoxification and stress-resistance, and innate immune systems are identified.

CONCLUSIONS

Our genomics and developmental transcriptomics analyses of N. cucumeris provide invaluable resources for further research on the development, reproduction, and fitness of this economically important mite in particular and Arachnida in general.

Characterisation of three novel β-1,3 glucanases from the medically important house dust mite Dermatophagoides pteronyssinus (airmid)

The European house dust mite, Dermatophagoides pteronyssinus is a major source of airborne allergens worldwide and is found in half of European homes. Interactions between microbes and house dust mites (HDM) are considered important factors that allow them to persist in the home. Laboratory studies indicate the European HDM, D. pteronyssinus is a mycophagous mite, capable of utilising a variety of fungi for nutrients, however specific mycolytic digestive enzymes are unknown. Our previous work identified a number of putative glycosyl hydrolases present in the predicted proteome of D. pteronyssinus airmid and validated the expression of 42 of these. Of note, three GH16 proteins with predicted β-1,3 glucanase activity were found to be consistently present in the mite body and excretome. Here, we performed an extensive bioinformatic, proteomic and biochemical study to characterize three-novel β-1,3 glucanases from this medically important house dust mite. The genes encoding novel β-1,3 glucanases designated Glu1, Glu2 and Glu3 were identified in D. pteronyssinus airmid, each exhibited more than 59% amino acid identity to one another. These enzymes are encoded by Glu genes present in a tri-gene cluster and protein homologs are found in other acari. The patchy phyletic distribution of Glu proteins means their evolutionary history remains elusive, however horizontal gene transfer cannot be completely excluded. Recombinant Glu1 and Glu2 exhibit hydrolytic activity toward laminarin, pachyman and barley glucan. Excreted β-1,3 glucanase activity was increased in response to D. pteronyssinus airmid feeding on baker's yeast. Active β-1,3 glucanases are expressed and excreted in the faeces of D. pteronyssinus airmid indicating they are digestive enzymes capable of breaking down β-1,3 glucans of fungi present in house dust.

Reproduction of Distinct Varroa destructor Genotypes on Honey Bee Worker Brood

Honey bees play important roles in pollination for many crops and wild plants, but have been facing great threats posed by various pathogens and parasites. Among them, Varroa destructor, an obligate ectoparasite of honey bees, is considered the most damaging. Within the last century, V. destructor shifted from the original host, the Asian honey bee Apis cerana to the new host, the European honey bee A. mellifera. However, the reproduction of Varroa mites, especially of different haplotypes in the two hosts, is still largely unknown. In this study, we first investigated the existing Varroa haplotypes in local colonies in southern China, and then compared the reproduction of different haplotypes on the worker brood of both the original and new hosts by artificial inoculation. We confirmed that there are two haplotypes of V. destructor in southern China, one is the Korea haplotype and the other is the China haplotype, and the two types parasitized different honey bee species. Although Varroa females from A. mellifera (Korea haplotype) are able to reproduce on the worker brood of both honey bee species, they showed better reproductive performance in the new host A. mellifera with significantly higher fecundity (number of offspring per mother mite) and reproductive rate (number of adult daughters per mother mite), suggesting that this parasite gains higher fitness after host shift. The data further showed that a short stay of Varroa females inside the A. cerana worker cells decreased their fecundity and especially the reproductive rate in a time-dependent manner, suggesting that the A. cerana worker cells may inhibit Varroa reproduction. In contrast, Varroa mites derived from A. cerana colonies (China haplotype) were entirely sterile in A. mellifera worker cells during two sequential inoculations, while the control mites from A. mellifera colonies (Korea haplotype) reproduced normally. In addition, all the infertile mites were found to defecate on the abdomen of bee pupae. We have revealed that two haplotypes of V. destructor exhibit differential reproduction on the worker brood of the original and new host honey bees, providing novel insights into the diversity and complexity of the reproduction of V. destructor.

In silico chemical library screening and experimental validation of novel compounds with potential varroacide activities

The mite Varroa destructor is an ectoparasite and has been identified as a major cause of worldwide honey bee colony losses. The use of yearly treatments for the control of varroosis is the most common answer to prevent collapses of honey bee colonies due to the mite. However, the number of effective acaricides is small and the mite tends to become resistant to these few active molecules. In this study, we have been looking for a new original varroacide treatment inhibiting selectively Varroa destructor AChE (vdAChE) with respect to Apis mellifera AChE (amAChE). To do this an original drug design methodology was used applying virtual screening of the CERMN chemolibrary, starting from a vdAChE homology sequence model. By combining the in silico screening with in vitro experiments, two promising compounds were found. In vitro tests of AChE inhibition for both species have confirmed good selectivity toward the mite vdAChE. Moreover, an in vivo protocol was performed and highlighted a varroacide activity without acute consequences on honey bee survival. The two compounds discovered have the potential to become new drug leads for the development of new treatments against the mite varroa. The method described here clearly shows the potential of a drug-design approach to develop new solutions to safeguard honey bee health.

Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities.

Deformed wing virus type A, a major honey bee pathogen, is vectored by the mite Varroa destructor in a non-propagative manner

Honey bees, the primary managed insect pollinator, suffer considerable losses due to Deformed wing virus (DWV), an RNA virus vectored by the mite Varroa destructor. Mite vectoring has resulted in the emergence of virulent DWV variants. The basis for such changes in DWV is poorly understood. Most importantly, it remains unclear whether replication of DWV occurs in the mite. In this study, we exposed Varroa mites to DWV type A via feeding on artificially infected honey bees. A significant, 357-fold increase in DWV load was observed in these mites after 2 days. However, after 8 additional days of passage on honey bee pupae with low viral loads, the DWV load dropped by 29-fold. This decrease significantly reduced the mites’ ability to transmit DWV to honey bees. Notably, negative-strand DWV RNA, which could indicate viral replication, was detected only in mites collected from pupae with high DWV levels but not in the passaged mites. We also found that Varroa mites contain honey bee mRNAs, consistent with the acquisition of honey bee cells which would additionally contain DWV replication complexes with negative-strand DWV RNA. We propose that transmission of DWV type A by Varroa mites occurs in a non-propagative manner.

Population genetics of ectoparasitic mites suggest arms race with honeybee hosts

The ectoparasitic mite, Varroa destructor, is the most severe biotic threat to honeybees (Apis mellifera) globally, usually causing colony death within a few years without treatments. While it is known that a few A. mellifera populations survive mite infestations by means of natural selection, the possible role of mite adaptations remains unclear. To investigate potential changes in mite populations in response to host adaptations, the genetic structure of V. destructor in the mite-resistant A. mellifera population on Gotland, Sweden, was studied. Spatio-temporal genetic changes were assessed by comparing mites collected in these colonies, as well as from neighboring mite-susceptible colonies, in historic (2009) and current (2017/2018) samples. The results show significant changes in the genetic structure of the mite populations during the time frame of this study. These changes were more pronounced in the V. destructor population infesting the mite-resistant honeybee colonies than in the mite-susceptible colonies. These results suggest that V. destructor populations are reciprocating, in a coevolutionary arms race, to the selection pressure induced by their honeybee host. Our data reveal exciting new insights into host-parasite interactions between A. mellifera and its major parasite.

Varroa chemosensory proteins: some are conserved across Arthropoda but others are arachnid specific

The tight synchronization between the life cycle of the obligatory parasitic mite Varroa destructor (Varroa) and its host, the honeybee, is mediated by honeybee chemical stimuli. These stimuli are mainly perceived by a pit organ located on the distal part of the mite's foreleg. In the present study, we searched for Varroa chemosensory molecular components by comparing transcriptomic and proteomic profiles between forelegs from different physiological stages, and rear legs. In general, a comparative transcriptomic analysis showed a clear separation of the expression profiles between the rear legs and the three groups of forelegs (phoretic, reproductive and tray-collected mites). Most of the differentially expressed transcripts and proteins in the mite's foreleg were previously uncharacterized. Using a conserved domain approach, we identified 45 transcripts with known chemosensory domains belonging to seven chemosensory protein families, of which 14 were significantly upregulated in the mite's forelegs when compared to rear legs. These are soluble and membrane bound proteins, including the somewhat ignored receptors of degenerin/epithelial Na+ channels and transient receptor potentials. Phylogenetic clustering and expression profiles of the putative chemosensory proteins suggest their role in chemosensation and shed light on the evolution of these proteins in Chelicerata.

Honey bee predisposition of resistance to ubiquitous mite infestations

Host-parasite co-evolution history is lacking when parasites switch to novel hosts. This was the case for Western honey bees (Apis mellifera) when the ectoparasitic mite, Varroa destructor, switched hosts from Eastern honey bees (Apis cerana). This mite has since become the most severe biological threat to A. mellifera worldwide. However, some A. mellifera populations are known to survive infestations, largely by suppressing mite population growth. One known mechanism is suppressed mite reproduction (SMR), but the underlying genetics are poorly understood. Here, we take advantage of haploid drones, originating from one queen from the Netherlands that developed Varroa-resistance, whole exome sequencing and elastic-net regression to identify genetic variants associated with SMR in resistant honeybees. An eight variants model predicted 88% of the phenotypes correctly and identified six risk and two protective variants. Reproducing and non-reproducing mites could not be distinguished using DNA microsatellites, which is in agreement with the hypothesis that it is not the parasite but the host that adapted itself. Our results suggest that the brood pheromone-dependent mite oogenesis is disrupted in resistant hosts. The identified genetic markers have a considerable potential to contribute to a sustainable global apiculture.

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host-parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid-20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural-selection-based breeding programmes have resulted in the evolution of Varroa-resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high-density genome-wide association analysis in a Varroa-resistant honey bee population, we identify an ecdysone-induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone-related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone-linked resistance genes to Varroa's meals and reproduction. If Varroa co-opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.

Multiple Virus Infections in Western Honeybee (Apis mellifera L.) Ejaculate Used for Instrumental Insemination

Instrumental insemination of Apis mellifera L. queens is a widely employed technique used in honeybee breeding that enables the effective control of mating. However, drone semen represents a potential source of honeybee viruses. In this study, 43 semen doses collected from apparently healthy drones, and consequently used in instrumental insemination, were analysed using PCR or RT-PCR to detect the presence of viral genome of 11 honeybee viruses. In 91% of samples, viral infection was detected. The survey revealed genomes of five viruses, namely Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Black queen cell virus (BQCV), Sacbrood virus (SBV), and A. mellifera filamentous virus (AmFV) in 84%, 19%, 14%, 2%, and 67% of samples, respectively. Single infection (30% of samples) as well as multiple infection (61% of samples) of two, three or four pathogens were also evaluated. To the best of our knowledge, this is the first study describing the presence of the BQCV and SBV genome sequence in drone ejaculate. Phylogenetic analysis of BQCV partial helicase gene sequence revealed the high similarity of nucleotide sequence of described Czech strains, which varied from 91.4% to 99.6%. The findings of our study indicate the possibility of venereal transmission of BQCV and SBV.

Association of Varroa destructor females in multiply infested cells of the honeybee Apis mellifera

The genetic diversity of Varroa destructor (Anderson & Trueman) is limited outside its natural range due to population bottlenecks and its propensity to inbreed. In light of the arms race between V. destructor and its honeybee (Apis mellifera L.) host, any mechanism enhancing population admixture of the mite may be favored. One way that admixture can occur is when two genetically dissimilar mites coinvade a brood cell, with the progeny of the foundresses admixing. We determined the relatedness of 393 pairs of V. destructor foundresses, each pair collected from a single bee brood cell (n = five colonies). We used six microsatellites to identify the genotypes of mites coinvading a cell and calculated the frequency of pairs with different or the same genotypes. We found no deviation from random coinvasion, but the frequency of cells infested by mites with different genotypes was high. This rate of recombination, coupled with a high transmission rate of mites, homogenized the allelic pool of mites within the apiary.

Haem Biology in Metazoan Parasites - 'The Bright Side of Haem'

Traditionally, host haem has been recognized as a cytotoxic molecule that parasites need to eliminate or detoxify in order to survive. However, recent evidence indicates that some lineages of parasites have lost genes that encode enzymes involved specifically in endogenous haem biosynthesis. Such lineages thus need to acquire and utilize haem originating from their host animal, making it an indispensable molecule for their survival and reproduction. In multicellular parasites, host haem needs to be systemically distributed throughout their bodies to meet the haem demands in all cell and tissue types. Host haem also gets deposited in parasite eggs, enabling embryogenesis and reproduction. Clearly, a better understanding of haem biology in multicellular parasites should elucidate organismal adaptations to obligatory blood-feeding.

Transcriptome profiling of the honeybee parasite Varroa destructor provides new biological insights into the mite adult life cycle

Background

The parasite Varroa destructor represents a significant threat to honeybee colonies. Indeed, development of Varroa infestation within colonies, if left untreated, often leads to the death of the colony. Although its impact on bees has been extensively studied, less is known about its biology and the functional processes governing its adult life cycle and adaptation to its host. We therefore developed a full life cycle transcriptomic catalogue in adult Varroa females and included pairwise comparisons with males, artificially-reared and non-reproducing females (10 life cycle stages and conditions in total).

Results

Extensive remodeling of the Varroa transcriptome was observed, with an upregulation of energetic and chitin metabolic processes during the initial and final phases of the life cycle (e.g. phoretic and post-oviposition stages), whereas during reproductive stages in brood cells genes showing functions related to transcriptional regulation were overexpressed. Several neurotransmitter and neuropeptide receptors involved in behavioural regulation, as well as active compounds of salivary glands, were also expressed at a higher level outside the reproductive stages. No difference was detected between artificially-reared phoretic females and their counterparts in colonies, or between females who failed to reproduce and females who successfully reproduced, indicating that phoretic individuals can be reared outside host colonies without impacting their physiology and that mechanisms underlying reproductive failure occur before oogenesis.

Conclusions

We discuss how these new findings reveal the remarkable adaptation of Varroa to its host biology and notably to the switch from living on adults to reproducing in sealed brood cells. By spanning the entire adult life cycle, our work captures the dynamic changes in the parasite gene expression and serves as a unique resource for deciphering Varroa biology and identifying new targets for mite control.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4668-z) contains supplementary material, which is available to authorized users.

Genetics and physiology of Varroa mites

Varroa destructor is the primary biological threat to domesticated honey bee colonies in much of the world, impacting host fitness both directly and by transmitting RNA viruses. Genomic, proteomic, and functional-genetic resources provide a framework for Varroa biology. When coupled with physiological analyses of development, host finding, and reproduction, these resources reveal general traits of arthropods and offer new strategies for mite control. Efforts to develop novel controls are focused on efficacy, efficient delivery, and the avoidance of both host impacts and the swift evolution of resistance by mites.

A Saliva Protein of Varroa Mites Contributes to the Toxicity toward Apis cerana and the DWV Elevation in A. mellifera

Varroa destructor mites express strong avoidance of the Apis cerana worker brood in the field. The molecular mechanism for this phenomenon remains unknown. We identified a Varroa toxic protein (VTP), which exhibited toxic activity toward A. cerana worker larvae, in the saliva of these mites, and expressed VTP in an Escherichia coli system. We further demonstrated that recombinant VTP killed A. cerana worker larvae and pupae in the absence of deformed-wing virus (DWV) but was not toxic to A. cerana worker adults and drones. The recombinant VTP was safe for A. mellifera individuals, but resulted in elevated DWV titers and the subsequent development of deformed-wing adults. RNAi-mediated suppression of vtp gene expression in the mites partially protected A. cerana larvae. We propose a modified mechanism for Varroa mite avoidance of worker brood, due to mutual destruction stress, including the worker larvae blocking Varroa mite reproduction and Varroa mites killing worker larvae by the saliva toxin. The discovery of VTP should provide a better understanding of Varroa pathogenesis, facilitate host-parasite mechanism research and allow the development of effective methods to control these harmful mites.

Chemosensing of honeybee parasite, Varroa destructor: Transcriptomic analysis

Chemosensing is a primary sense in nature, however little is known about its mechanism in Chelicerata. As a model organism we used the mite Varroa destructor, a key parasite of honeybees. Here we describe a transcriptomic analysis of two physiological stages for the Varroa foreleg, the site of primary olfactory organ. The transcriptomic analysis revealed transcripts of chemosensory related genes belonging to several groups. These include Niemann-Pick disease protein, type C2 (NPC2), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and odorant binding proteins (OBP). However, no insect odorant receptors (ORs) and odorant co-receptors (ORcos) were found. In addition, we identified a homolog of the most ancient IR co-receptor, IR25a, in Varroa as well as in other members of Acari. High expression of this transcript in the mite's forelegs, while not detectable in the other pairs of legs, suggests a function for this IR25a-like in Varroa chemosensing.

A Varroa destructor protein atlas reveals molecular underpinnings of developmental transitions and sexual differentiation

Varroa destructor is the most economically damaging honey bee pest, weakening colonies by simultaneously parasitizing bees and transmitting harmful viruses. Despite these impacts on honey bee health, surprisingly little is known about its fundamental molecular biology. Here, we present a Varroa protein atlas crossing all major developmental stages (egg, protonymph, deutonymph, and adult) for both male and female mites as a web-based interactive tool (http://foster.nce.ubc.ca/varroa/index.html). We used intensity-based label-free quantitation to find 1,433 differentially expressed proteins across developmental stages. Enzymes for processing carbohydrates and amino acids were among many of these differences as well as proteins involved in cuticle formation. Lipid transport involving vitellogenin was the most significantly enriched biological process in the foundress (reproductive female) and young mites. In addition, we found that 101 proteins were sexually regulated and functional enrichment analysis suggests that chromatin remodeling may be a key feature of sex determination. In a proteogenomic effort, we identified 519 protein-coding regions, 301 of which were supported by two or more peptides and 169 of which were differentially expressed. Overall, this work provides a first-of-its-kind interrogation of the patterns of protein expression that govern the Varroa life cycle and the tools we have developed will support further research on this threatening honey bee pest.

Draft genome of the honey bee ectoparasitic mite, Tropilaelaps mercedesae, is shaped by the parasitic life history

The number of managed honey bee colonies has considerably decreased in many developed countries in recent years and ectoparasitic mites are considered as major threats to honey bee colonies and health. However, their general biology remains poorly understood. We sequenced the genome of Tropilaelaps mercedesae, the prevalent ectoparasitic mite infesting honey bees in Asia, and predicted 15 190 protein-coding genes that were well supported by the mite transcriptomes and proteomic data. Although amino acid substitutions have been accelerated within the conserved core genes of two mites, T. mercedesae and Metaseiulus occidentalis, T. mercedesae has undergone the least gene family expansion and contraction between the seven arthropods we tested. The number of sensory system genes has been dramatically reduced, but T. mercedesae contains all gene sets required to detoxify xenobiotics. T. mercedesae is closely associated with a symbiotic bacterium (Rickettsiella grylli-like) and Deformed Wing Virus, the most prevalent honey bee virus. T. mercedesae has a very specialized life history and habitat as the ectoparasitic mite strictly depends on the honey bee inside a stable colony. Thus, comparison of the genome and transcriptome sequences with those of a tick and free-living mites has revealed the specific features of the genome shaped by interaction with the honey bee and colony environment. Genome and transcriptome sequences of T. mercedesae, as well as Varroa destructor (another globally prevalent ectoparasitic mite of honey bee), not only provide insights into the mite biology, but may also help to develop measures to control the most serious pests of the honey bee.

Changes in the Bacteriome of Honey Bees Associated with the Parasite Varroa destructor, and Pathogens Nosema and Lotmaria passim

The honey bee, Apis mellifera, is a globally important species that suffers from a variety of pathogens and parasites. These parasites and pathogens may have sublethal effects on their bee hosts via an array of mechanisms, including through a change in symbiotic bacterial taxa. Our aim was to assess the influence of four globally widespread parasites and pathogens on the honey bee bacteriome. We examined the effects of the ectoparasitic mite Varroa destructor, the fungal pathogens Nosema apis and Nosema ceranae, and the trypanosome Lotmaria passim. Varroa was detected by acaricidal treatment, Nosema and L. passim by PCR, and the bacteriome using MiSeq 16S rRNA gene sequencing. Overall, the 1,858,850 obtained sequences formed 86 operational taxonomic units (OTUs) at 3 % dissimilarity. Location, time of year, and degree of infestation by Varroa had significant effects on the composition of the bacteriome of honey bee workers. Based on statistical correlations, we found varroosis more important factor than N. ceranae, N. apis, and L. passim infestation influencing the honey bee bacteriome and contributing to the changes in the composition of the bacterial community in adult bees. At the population level, Varroa appeared to modify 20 OTUs. In the colonies with high Varroa infestation levels (varroosis), the relative abundance of the bacteria Bartonella apis and Lactobacillus apis decreased. In contrast, an increase in relative abundance was observed for several taxa including Lactobacillus helsingborgensis, Lactobacillus mellis, Commensalibacter intestini, and Snodgrassella alvi. The results showed that the "normal" bacterial community is altered by eukaryotic parasites as well as displaying temporal changes and changes associated with the geographical origin of the beehive.

Seasonal cycle of inbreeding and recombination of the parasitic mite Varroa destructor in honeybee colonies and its implications for the selection of acaricide resistance

Varroa destructor is the most devastating parasite of the Western honeybee, Apis mellifera. In the light of the arm race opposing the host and its parasite, the population dynamics and genetic diversity of these organisms are key parameters. However, the life cycle of V. destructor is characterized by extreme inbreeding due to full sibling mating in the host brood cells. We here present an equation reflecting the evolution of inbreeding in such a clonal system, and compare our predictions with empirical data based on the analysis of seven microsatellite markers. This comparison revealed that the mites perform essentially incestuous mating in the beginning of the brood season. However, this pattern changes with the development of mite infestation. Despite the fact that the overall level of genetic diversity of the mites remained low through the season, multiple inbred lineages were identified in the mites we sampled in June. As a response to the decrease of brood availability and the increase of the parasite population in parallel in the colonies, these lineages recombined towards the end of the season as mites co-infest brood cells. Our results suggest that the ratio of the number of mite per brood cell in the colony determines the genetic structure of the populations of V. destructor. This intracolonial population dynamics has great relevance for the selection of acaricide resistance in V. destructor. If chemical treatments occur before the recombination phase, inbreeding will greatly enhance the fixation of resistance alleles at the colony level.

Two novel viruses associated with the Apis mellifera pathogenic mite Varroa destructor

Varroa destructor infestation of Apis mellifera colonies carries and/or promotes replication of honey bee viruses like the Deformed wing virus, the Varroa destructor virus-1, the Acute bee paralysis virus, the Israeli acute bee paralysis virus and the Kashmir bee virus that have been well described and characterized; but viruses exclusively associated with Varroa were not found. To look for viruses that may associate with- or infect V. destructor we performed deep sequencing (RNA-seq) of RNA extracted from honey bees and mites in Varroa-infested untreated colonies. Comparative bioinformatic analysis of the two separate contig-assemblies generated from the sequences' reads annotated using Blastx enabled identification of new viruses unique to Varroa and absent in A. mellifera: an Iflavirus and a virus with homology to Ixodes scapularis associated virus 2, that we named Varroa destructor virus 2 (VDV-2) and 3(VDV-3), respectively. We validated these findings sequencing the mite- and honey bee-viromes and in separate mites and honey bees randomly sampled. The complete genomes of VDV-2 and VDV-3 bear 9576 nucleotides and 4202 nucleotides, respectively. Phylogenetic analysis of VDV-3 suggests that it belongs to a new group of viruses. Our results open venues for investigating the pathogenicity of these V. destructor viruses.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Comparison of Varroa destructor and Worker Honeybee Microbiota Within Hives Indicates Shared Bacteria

The ectoparasitic mite Varroa destructor is a major pest of the honeybee Apis mellifera. In a previous study, bacteria were found in the guts of mites collected from winter beehive debris and were identified using Sanger sequencing of their 16S rRNA genes. In this study, community comparison and diversity analyses were performed to examine the microbiota of honeybees and mites at the population level. The microbiota of the mites and honeybees in 26 colonies in seven apiaries in Czechia was studied. Between 10 and 50 Varroa females were collected from the bottom board, and 10 worker bees were removed from the peripheral comb of the same beehive. Both bees and mites were surface sterilized. Analysis of the 16S rRNA gene libraries revealed significant differences in the Varroa and honeybee microbiota. The Varroa microbiota was less diverse than was the honeybee microbiota, and the relative abundances of bacterial taxa in the mite and bee microbiota differed. The Varroa mites, but not the honeybees, were found to be inhabited by Diplorickettsia. The relative abundance of Arsenophonus, Morganella, Spiroplasma, Enterococcus, and Pseudomonas was higher in Varroa than in honeybees, and the Diplorickettsia symbiont detected in this study is specific to Varroa mites. The results demonstrated that there are shared bacteria between Varroa and honeybee populations but that these bacteria occur in different relative proportions in the honeybee and mite bacteriomes. These results support the suggestion of bacterial transfer via mites, although only some of the transferred bacteria may be harmful.

Genome Sequencing of the Phytoseiid Predatory Mite Metaseiulus occidentalis Reveals Completely Atomized Hox Genes and Superdynamic Intron Evolution

Metaseiulus occidentalis is an eyeless phytoseiid predatory mite employed for the biological control of agricultural pests including spider mites. Despite appearances, these predator and prey mites are separated by some 400 Myr of evolution and radically different lifestyles. We present a 152-Mb draft assembly of the M. occidentalis genome: Larger than that of its favored prey, Tetranychus urticae, but considerably smaller than those of many other chelicerates, enabling an extremely contiguous and complete assembly to be built-the best arachnid to date. Aided by transcriptome data, genome annotation cataloged 18,338 protein-coding genes and identified large numbers of Helitron transposable elements. Comparisons with other arthropods revealed a particularly dynamic and turbulent genomic evolutionary history. Its genes exhibit elevated molecular evolution, with strikingly high numbers of intron gains and losses, in stark contrast to the deer tick Ixodes scapularis Uniquely among examined arthropods, this predatory mite's Hox genes are completely atomized, dispersed across the genome, and it encodes five copies of the normally single-copy RNA processing Dicer-2 gene. Examining gene families linked to characteristic biological traits of this tiny predator provides initial insights into processes of sex determination, development, immune defense, and how it detects, disables, and digests its prey. As the first reference genome for the Phytoseiidae, and for any species with the rare sex determination system of parahaploidy, the genome of the western orchard predatory mite improves genomic sampling of chelicerates and provides invaluable new resources for functional genomic analyses of this family of agriculturally important mites.

A Toolbox for Quantitative Gene Expression in Varroa destructor: RNA Degradation in Field Samples and Systematic Analysis of Reference Gene Stability

Varroa destructor is the major pest of Apis mellifera and contributes to the global honey bee health crisis threatening food security. Developing new control strategies to combat Varroa will require the application of molecular biology, including gene expression studies by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Both high quality RNA samples and suitable stable internal reference genes are required for accurate gene expression studies. In this study, ten candidate genes (succinate dehydrogenase (SDHA), NADH dehydrogenase (NADH), large ribsosmal subunit, TATA-binding protein, glyceraldehyde-3-phosphate dehydrogenase, 18S rRNA (18S), heat-shock protein 90 (HSP90), cyclophilin, α-tubulin, actin), were evaluated for their suitability as normalization genes using the geNorm, Normfinder, BestKeeper, and comparative ΔCq algorithims. Our study proposes the use of no more than two of the four most stable reference genes (NADH, 18S, SDHA and HSP90) in Varroa gene expression studies. These four genes remain stable in phoretic and reproductive stage Varroa and are unaffected by Deformed wing virus load. When used for determining changes in vitellogenin gene expression, the signal-to-noise ratio (SNR) for the relatively unstable genes actin and α-tubulin was much lower than for the stable gene combinations (NADH + HSP90 +18S; NADH + HSP90; or NADH). Using both electropherograms and RT-qPCR for short and long amplicons as quality controls, we demonstrate that high quality RNA can be recovered from Varroa up to 10 days later stored at ambient temperature if collected into RNAlater and provided the body is pierced. This protocol allows the exchange of Varroa samples between international collaborators and field sample collectors without requiring frozen collection or shipping. Our results make important contributions to gene expression studies in Varroa by proposing a validated sampling protocol to obtain high quality Varroa RNA and the validation of suitable reference genes for expression studies in this globally important pest.

Comparison of tau-fluvalinate, acrinathrin, and amitraz effects on susceptible and resistant populations of Varroa destructor in a vial test

The parasitic mite Varroa destructor is a major pest of the western honeybee, Apis mellifera. The development of acaricide resistance in Varroa populations is a global issue. Discriminating concentrations of acaricides are widely used to detect pest resistance. Two methods, using either glass vials or paraffin capsules, are used to screen for Varroa resistance to various acaricides. We found the glass vial method to be useless for testing Varroa resistance to acaridices, so we developed a polypropylene vial bioassay. This method was tested on tau-fluvalinate-, acrinathrin-, and amitraz-resistant mite populations from three apiaries in Czechia. Acetone was used as a control and technical grade acaricide compounds diluted in acetone were applied to the polypropylene vials. The solutions were spread on the vial surface by rolling the vial, and were then evaporated. Freshly collected Varroa females were placed in the vials and the mortality of the exposed mites was measured after 24 h. The Varroa populations differed in mortality between the apiaries and the tested compounds. Mites from the Kyvalka site were resistant to acrinathrin, tau-fluvalinate, and amitraz, while mites from the Postrizin site were susceptible to all three acaricides. In Prelovice apiary, the mites were susceptible to acrinathrin and amitraz, but not to tau-fluvalinate. The calculated discriminating concentrations for tau-fluvalinate, acrinathrin, and amitraz were 0.66, 0.26 and 0.19 µg/mL, respectively. These results indicate that polyproplyne vial tests can be used to determine discriminating concentrations for the early detection of acaricide resistant Varroa. Finally, multiple-resistance in Kyvalka may indicate metabolic resistance.

Identification and gene-silencing of a putative odorant receptor transcription factor in Varroa destructor: possible role in olfaction

The ectoparasitic mite Varroa destructor is one of the major threats to apiculture. Using a behavioural choice bioassay, we determined that phoretic mites were more successful in reaching a bee than reproductive mites, suggesting an energy trade-off between reproduction and host selection. We used both chemo-ecological and molecular strategies to identify the regulation of the olfactory machinery of Varroa and its association with reproduction. We focused on transcription regulation. Using primers designed to the conserved DNA binding region of transcription factors, we identified a gene transcript in V. destructor homologous to the pheromone receptor transcription factor (PRTF) gene of Pediculus humanus corporis. Quantitative PCR (qPCR) revealed that this PRTF-like gene transcript is expressed in the forelegs at higher levels than in the body devoid of forelegs. Subsequent comparative qPCR analysis showed that transcript expression was significantly higher in the phoretic as compared to the reproductive stage. Electrophysiological and behavioural studies revealed a reduction in the sensitivity of PRTF RNA interference-silenced mites to bee headspace, consistent with a reduction in the mites' ability to reach a host. In addition, vitellogenin expression was stimulated in PRTF-silenced mites to similar levels as found in reproductive mites. These data shed light upon the regulatory mechanism of host chemosensing in V. destructor.

Transcriptome analysis of the synganglion from the honey bee mite, Varroa destructor and RNAi knockdown of neural peptide targets

Varroa mites (Varroa destructor) and the viruses that they transmit are one of the major contributing factors to the global honey bee crisis. Gene products within the nervous system are the targets of all the insecticides currently used to control Varroa but there is a paucity of transcriptomic data available for Varroa neural tissues. A cDNA library from the synganglia ("brains") of adult female Varroa was constructed and 600 ESTs sequenced and analysed revealing several current and potential druggable targets. Contigs coding for the deformed wing virus (DWV) variants V. destructor virus-1 (VDV-1) and the recombinant (VDV-1DVD) were present in the synganglion library. Negative-sense RNA-specific PCR indicated that VDV-1 replicates in the Varroa synganglion and all other tissues tested, but we could not detect DWV replicating in any Varroa tissue. Two neuropeptides were identified in the synganlion EST library: a B-type allatostatin and a member of the crustacean hyperglycaemic hormone (CHH) superfamily. Knockdown of the allatostatin or the CHH-like gene by double-stranded RNA-interference (dsRNAi) resulted in 85% and 55% mortality, respectively, of Varroa. Here, we present the first transcriptomic survey in Varroa and demonstrate that neural genes can be targeted by dsRNAi either for genetic validation of putative targets during drug discovery programmes or as a potential control measure in itself.

Natural Larval Diet Differently Influences the Pattern of Developmental Changes in DNA 5-Methylcytosine Levels in Apis mellifera Queens as Compared with Workers and Drones

The principal mechanism of gene activation/silencing is DNA 5-methylcytosine methylation. This study was aimed at determining global DNA methylation levels in larvae, prepupae, pupae, and 1-day-old adults of Apis mellifera queens, workers and drones. The Imprint Methylated DNA Quantification Kit MDQ1 was used. Percentages of DNA 5-methylcytosine were low and relatively similar in the larvae of all the castes until 4th day of larval development (3-5%). However, they were higher in the drone and worker larvae than in the queen larvae. Generally, the developmental patterns of changes in the DNA methylation levels were different in the queens in comparison with the drones and workers. While methylation increased in the queens, it decreased in the drones and workers. Methylated DNA methylcytosine percentages and weights in the queen prepupae (15%, 9.18 ng) and pupae (21%, 10.74 ng) were, respectively, three and four times higher than in the worker/drone brood of the same age (2.5-4%, 0.03-0.07 ng). Only in the queens, after a substantial increase, did DNA methylation decrease almost twice between the pupal stage and queen emergence (from 21% and 10.74 ng to 12% and 6.78 ng). This finding seems very interesting, particularly for experimental gerontology.