Effects of local domestication warrant attention in honey bee population genetics.

Honey bees, Apis mellifera, have for millennia been managed and exploited by humans and introduced into most suitable regions worldwide. However, given the lack of records for many introduction events, treating A. mellifera populations as native would predictably bias genetic studies regarding origin and evolution. Here, we used the Dongbei bee, a well-documented population, introduced beyond the natural distribution range approximately 100 years ago, to elucidate the effects of local domestication on animal population genetic analyses. Strong domestication pressure was detected in this population, and the genetic divergence between Dongbei bee and its ancestral subspecies was found to have occurred at the lineage level. Results of phylogenetic and time divergence analyses could consequently be misinterpreted. Proposing new subspecies or lineages and performing analyses of origin should thus strive to eliminate anthropogenic effects. We highlight the need for definitions of landrace and breed in honey bee sciences and make preliminary suggestions.

Acute and chronic viruses mediated by an ectoparasite targeting different developmental stages of honeybee (Apis mellifera and Apis cerana) brood.

The health of the western honeybee, Apis mellifera, the most crucial pollinator, has been challenged globally over the past decades. An ectoparasitic mite, Varroa destructor, together with the viruses it vectored, is generally regarded as the vital pathogenic agent. Although the poor health status of A. mellifera compared to its eastern counterpart, Apis cerana, has been broadly identified, the underlying mechanism remains poorly understood and comparison between susceptible and resistant hosts will potentially ameliorate this predicament. Here, we investigated the impacts of two widespread viruses-deformed wing virus type A (DWV-A) and Israeli acute paralysis virus (IAPV), mediated by V. destructor mite, on the capped developing honeybee brood, in the absence of adult workers, of A. mellifera and A. cerana, with positive and negative controls. Our results demonstrated that the endogenous viruses imposed limited damage on the hosts even if the brood was wounded. In contrast, the exogenous viruses introduced by ectoparasites triggered variable mortality of the infested brood between host species. Intriguingly, death causes of both honeybee species presented a similar trend: the acute IAPV generally causes morbidity and mortality of late larvae, while the chronic DWV-A typically leads to brood mortality during and after pupation. Notably, the susceptible immature A. cerana individuals, supported by higher observed mortality and a lower virus tolerance, serve the interests of the colony and foster the overall survival of a resistant honeybee superorganism. These results improve our understanding of the interactions between viruses carried by ectoparasites and their developing hosts, and the novel insight of weak individuals fostering strong colonies may promote breeding efforts to mitigate the indefensible colony losses globally.

Investigation of circular RNAs in an ectoparasitic mite Varroa destructor (Acarina: Varroidae) of the honey bee.

Circular RNAs (circRNAs) are a large class of non-protein-coding transcripts that are involved in a diverse spectrum of regulatory mechanisms across a broad range of biological processes. To date, however, few studies on circRNAs have investigated their role in the biology of invertebrate parasites. The ectoparasitic mite Varroa destructor is perceived as the principal biotic threat towards global honey bee health. This parasite cannot be sustainably controlled partially due to the lack of knowledge about its basic molecular biology. In this paper, we unveil the circRNA profile of V. destructor for the first time and report the sources, distribution, and features of the identified circRNAs. Exonic, intronic, exon-intron, and intergenic circRNAs were discovered and exon-intron circRNAs were the most abundant within the largest spliced length. Three hundred and eighty-six (8.3%) circRNAs were predicted to possess translational potential. Eleven circRNAs, derived from six parental genes, exhibited strong bonds with miRNAs as sponges, suggesting an efficient post-transcriptional regulation. GO term and KEGG pathway enrichment analyses of the parental genes of the identified circRNAs showed that these non-coding RNAs were mainly engaged in protein processing, signal transduction, and various metabolism processes. To our knowledge, this is the first catalog of a circRNA profile of parasitiformes species, which reveals the prevalence of circRNAs in the parasite and provides biological insights for future genetic studies on this ubiquitous parasitic mite.