Evolutionarily diverse origins of deformed wing viruses in western honey bees.

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

The first steps toward a global pandemic: Reconstructing the demographic history of parasite host switches in its native range.

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real-time, and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after being brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploiting the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

A toolkit for studying Varroa genomics and transcriptomics: preservation, extraction, and sequencing library preparation.

BACKGROUND: The honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity. RESULTS: For both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4 °C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms. CONCLUSIONS: Ethanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

When European meets African honeybees (Apis mellifera L.) in the tropics: Morphological changes related to genetics in Mauritius Island (South-West Indian Ocean).

The previous genetic characterization of the honeybee population of Mauritius Island (Indian Ocean) revealed an ongoing process of hybridization between the first established African subspecies Apis mellifera unicolor and recently imported European subspecies (A. m. ligustica, A. m. carnica and A. m. mellifera). This context offers the rare opportunity to explore the influence of hybridization between African and European honeybees on phenotypic traits out of the case largely studied of the Africanized honeybee (hybrid between A. m. scutellata from South Africa and European subspecies). We thus conducted geometric morphometric analyses on forewings of 283 workers genetically characterized at 14 microsatellite loci to evaluate (1) if the morphological variability coincides well with the neutral genetic variability, (2) if hybrids exhibited rather parental, intermediate or transgressive traits, and (3) to test if fluctuating asymmetry (FA) of size and shape, as a measure of developmental stability, was elevated in hybrids (due to genetic stress) and/or European bees (due to unsuitable environment) compared to African bees. A strong concordance was found between morphological variability and neutral genetic variability, especially for wing shape, based on partial least-square analyses (PLS). However, on average, the morphology of hybrids was more similar to the African bees, potentially reflecting the dynamics and direction of introgression. Significant FA for wing size as well as wing shape was detected, suggesting the overall presence of stress during the development of the studied individuals. In contrast, the asymmetry levels do not differ according to the ancestry (African, European or hybrid) of the individuals. Therefore, if ongoing hybridization contributed to increasing the genetic and phenotypic diversity of the populations and influences its adaptive potential, developmental stressors could not be identified and their evolutionary consequences remain uncertain.

Varroa destructor: A Complex Parasite, Crippling Honey Bees Worldwide.

The parasitic mite, Varroa destructor, has shaken the beekeeping and pollination industries since its spread from its native host, the Asian honey bee (Apis cerana), to the naïve European honey bee (Apis mellifera) used commercially for pollination and honey production around the globe. Varroa is the greatest threat to honey bee health. Worrying observations include increasing acaricide resistance in the varroa population and sinking economic treatment thresholds, suggesting that the mites or their vectored viruses are becoming more virulent. Highly infested weak colonies facilitate mite dispersal and disease transmission to stronger and healthier colonies. Here, we review recent developments in the biology, pathology, and management of varroa, and integrate older knowledge that is less well known.

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.

Autosomal and Mitochondrial Adaptation Following Admixture: A Case Study on the Honeybees of Reunion Island.

The honeybee population of the tropical Reunion Island is a genetic admixture of the Apis mellifera unicolor subspecies, originally described in Madagascar, and of European subspecies, mainly A. m. carnica and A. m. ligustica, regularly imported to the island since the late 19th century. We took advantage of this population to study genetic admixing of the tropical-adapted indigenous and temperate-adapted European genetic backgrounds. Whole genome sequencing of 30 workers and 6 males from Reunion, compared with samples from Europe, Madagascar, Mauritius, Rodrigues, and the Seychelles, revealed the Reunion honeybee population to be composed on an average of 53.2 ± 5.9% A. m. unicolor nuclear genomic background, the rest being mainly composed of A. m. carnica and to a lesser extent A. m. ligustica. In striking contrast to this, only 1 out of the 36 honeybees from Reunion had a mitochondrial genome of European origin, suggesting selection has favored the A. m. unicolor mitotype, which is possibly better adapted to the island's bioclimate. Local ancestry was determined along the chromosomes for all Reunion samples, and a test for preferential selection for the A. m. unicolor or European background revealed 15 regions significantly associated with the A. m. unicolor lineage and 9 regions with the European lineage. Our results provide insights into the long-term consequences of introducing exotic specimen on the nuclear and mitochondrial genomes of locally adapted populations.

Genetic diversity and differentiation among insular honey bee populations in the southwest Indian Ocean likely reflect old geographical isolation and modern introductions.

With globalization the Western honey bee has become a nearly cosmopolitan species, but it was originally restricted to the Old World. This renowned model of biodiversity has diverged into five evolutionary lineages and several geographic "subspecies." If Apis mellifera unicolor is indubitably an African subspecies endemic to Madagascar, its relationship with honey bees from three archipelagos in the southwest Indian Ocean (SWIO) hotspot of biodiversity is misunderstood. We compared recent mtDNA diversity data to an original characterization of the nuclear diversity from honey bees in the Mascarenes and Comoros archipelagos, using 14 microsatellites, but also additional mtDNA tRNALeu-cox2 analysis. Our sampling offers the most comprehensive dataset for the SWIO populations with a total of 3,270 colonies from 10 islands compared with 855 samples from Madagascar, 113 from Africa, and 138 from Europe. Comprehensive mitochondrial screening confirmed that honey bees from La Réunion, Mauritius, and Comoros archipelagos are mainly of African origin (88.1% out of 2,746 colonies) and that coexistence with European lineages occurs only in the Mascarenes. PCA, Bayesian, and genetic differentiation analysis showed that African colonies are not significantly distinct on each island, but have diversified among islands and archipelagos. FST levels progressively decreased in significance from European and African continental populations, to SWIO insular and continental populations, and finally among islands from the same archipelago. Among African populations, Madagascar shared a nuclear background with and was most closely related to SWIO island populations (except Rodrigues). Only Mauritius Island presented clear cytoplasmic disequilibrium and genetic structure characteristic of an admixed population undergoing hybridization, in this case, between A. m. unicolor and A. m. ligustica, A. m. carnica and A. m. mellifera-like individuals. Finally, global genetic clustering analysis helped to better depict the colonization and introduction pattern of honey bee populations in these archipelagos.

Large-scale mitochondrial DNA analysis of native honey bee Apis mellifera populations reveals a new African subgroup private to the South West Indian Ocean islands.

BACKGROUND: The South West Indian Ocean (SWIO) archipelagos and Madagascar constitute a hotspot of biodiversity with a high rate of endemism. In this area, the endemic subspecies A. m. unicolor has been described in Madagascar. It belongs to the African lineage, one of the four described evolutionary lineages in honey bees. Despite a long beekeeping tradition and several recorded European introductions, few studies have been carried out on the diversity and proportion of honey bee subspecies. In order to identify and define which evolutionary lineages and potential sub-lineages are present in the SWIO, the COI-COII intergenic region and the ND2 gene of the mtDNA were sequenced in honey bee colonies from three archipelagos. An extensive sampling (n = 1184 colonies) was done in the Mascarene (La Réunion, Mauritius, Rodrigues), Seychelles (Mahé, Praslin, La Digue) and Comoros (Grande Comore, Mohéli, Anjouan, Mayotte) archipelagos. Islands genetic diversity was compared to newly sampled populations from Madagascar, continental African and European populations. RESULTS: African lineage haplotypes were found in all islands (except for Rodrigues). Madagascar, Comoros and Seychelles had 100% of A lineage, 95.5% in La Réunion and 56.1% in Mauritius. Among all African colonies detected in the SWIO, 98.1% (n = 633) of COI-COII haplotypes described the presence of the subspecies A. M. unicolor. Both genetic markers revealed i) a new private AI mitochondrial group shared by the SWIO archipelagos and Madagascar distant from continental populations; ii) the private African haplotypes for each island suggested diversity radiation in the archipelagos; iii) the detection of the Comoros archipelago as a possible contact area between insular and continental African populations. The exotic European C and M lineages were only detected in the Mascarene archipelago, but striking differences of proportion were observed among islands. Merely 4.6% of European colonies were found in La Réunion whereas Mauritius cumulated 44%. Here, among the 84 observed COI-COII haplotypes, 50 were newly described including 13 which were private to the SWIO archipelagos and Madagascar. Similarly, 24 of the 34 found ND2 haplotypes were novel which included six haplotypes particular to the SWIO populations. CONCLUSION: A new African subgroup was described in the SWIO region with mitochondrial genetic evidence that A. m. unicolor is the indigenous subspecies of the archipelagos surrounding Madagascar.