Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies

Over the past 50 years, many millions of European honey bee (Apis mellifera) colonies have died as the ectoparasitic mite, Varroa destructor, has spread around the world. Subsequent studies have indicated that the mite's association with a group of RNA viral pathogens (Deformed Wing Virus, DWV) correlates with colony death. Here, we propose a phenomenon known as superinfection exclusion that provides an explanation of how certain A. mellifera populations have survived, despite Varroa infestation and high DWV loads. Next-generation sequencing has shown that a non-lethal DWV variant 'type B' has become established in these colonies and that the lethal 'type A' DWV variant fails to persist in the bee population. We propose that this novel stable host-pathogen relationship prevents the accumulation of lethal variants, suggesting that this interaction could be exploited for the development of an effective treatment that minimises colony losses in the future.The ISME Journal advance online publication, 27 October 2015; doi:10.1038/ismej.2015.186.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Detection and identification of species-specific bacteria associated with synanthropic mites

Internal bacterial communities of synanthropic mites Acarus siro, Dermatophagoides farinae, Lepidoglyphus destructor, and Tyrophagus putrescentiae (Acari: Astigmata) were analyzed by culturing and culture-independent approaches from specimens obtained from laboratory colonies. Homogenates of surface-sterilized mites were used for cultivation on non-selective agar and DNA extraction. Isolated bacteria were identified by sequencing of the 16S rRNA gene. PCR amplified 16S rRNA genes were analyzed by terminal restriction fragment length polymorphism analysis (T-RFLP) and cloning sequencing. Fluorescence in situ hybridization using universal bacterial probes was used for direct bacterial localization. T-RFLP analysis of 16S rRNA gene revealed distinct species-specific bacterial communities. The results were further confirmed by cloning and sequencing (284 clones). L. destructor and D. farinae showed more diverse communities then A. siro and T. putrescentiae. In the cultivated part of the community, the mean CFUs from four mite species ranged from 5.2 × 102 to 1.4 × 103 per mite. D. farinae had significantly higher CFUs than the other species. Bacteria were located in the digestive and reproductive tract, parenchymatical tissue, and in bacteriocytes. Among the clones, Bartonella-like bacteria occurring in A. siro and T. putresecentiae represented a distinct group related to Bartonellaceae and to Bartonella-like symbionts of ants. The clones of high similarity to Xenorhabdus cabanillasii were found in L. destructor and D. farinae, and one clone related to Photorhabdus temperata in A. siro. Members of Sphingobacteriales cloned from D. farinae and A. siro clustered with the sequences of “Candidatus Cardinium hertigii” and as a separate novel cluster.