Prospects, challenges and perspectives in harnessing natural selection to solve the 'varroa problem' of honey bees.

Honey bees, Apis mellifera, of European origin are major pollinators of crops and wild flora. Their endemic and exported populations are threatened by a variety of abiotic and biotic factors. Among the latter, the ectoparasitic mite Varroa destructor is the most important single cause behind colony mortality. The selection of mite resistance in honey bee populations has been deemed a more sustainable solution to its control than varroacidal treatments. Because natural selection has led to the survival of some European and African honey bee populations to V. destructor infestations, harnessing its principles has recently been highlighted as a more efficient way to provide honey bee lineages that survive infestations when compared with conventional selection on resistance traits against the parasite. However, the challenges and drawbacks of harnessing natural selection to solve the varroa problem have only been minimally addressed. We argue that failing to consider these issues could lead to counterproductive results, such as increased mite virulence, loss of genetic diversity reducing host resilience, population collapses or poor acceptance by beekeepers. Therefore, it appears timely to evaluate the prospects for the success of such programmes and the qualities of the populations obtained. After reviewing the approaches proposed in the literature and their outcomes, we consider their advantages and drawbacks and propose perspectives to overcome their limitations. In these considerations, we not only reflect on the theoretical aspects of host-parasite relationships but also on the currently largely neglected practical constraints, that is, the requirements for productive beekeeping, conservation or rewilding objectives. To optimize natural selection-based programmes towards these objectives, we suggest designs based on a combination of nature-driven phenotypic differentiation and human-directed selection of traits. Such a dual strategy aims at allowing field-realistic evolutionary approaches towards the survival of V. destructor infestations and the improvement of honey bee health

Advances and perspectives in selecting resistance traits against the parasitic mite Varroa destructor in honey bees.

In spite of the implementation of control strategies in honey bee (Apis mellifera) keeping, the invasive parasitic mite Varroa destructor remains one of the main causes of colony losses in numerous countries. Therefore, this parasite represents a serious threat to beekeeping and agro-ecosystems that benefit from the pollination services provided by honey bees. To maintain their stocks, beekeepers have to treat their colonies with acaricides every year. Selecting lineages that are resistant to infestations is deemed to be a more sustainable approach. Over the last three decades, numerous selection programs have been initiated to improve the host-parasite relationship and to support honey bee survival in the presence of the parasite without the need for acaricide treatments. Although resistance traits have been included in the selection strategy of honey bees, it has not been possible to globally solve the V. destructor problem. In this study, we review the literature on the reasons that have potentially limited the success of such selection programs. We compile the available information to assess the relevance of selected traits and the potential environmental effects that distort trait expression and colony survival. Limitations to the implementation of these traits in the field are also discussed. Improving our knowledge of the mechanisms underlying resistance to V. destructor to increase trait relevance, optimizing selection programs to reduce environmental effects, and communicating selection outcomes are all crucial to efforts aiming at establishing a balanced relationship between the invasive parasite and its new host

The ectoparasitic mite Tropilaelaps mercedesae (Acari, Laelapidae) as a vector of honeybee viruses

Abstract.: The ectoparasitic mites Varroa destructor and Tropilaelaps mercedesae share life history traits and both infect honeybee colonies, Apis mellifera. Since V. destructor is a biological vector of several honeybee viruses, we here test whether T. mercedesae can also be infected and enable virus replication. In Kunming (China), workers and T. mercedesae mites were sampled from three A. mellifera colonies, where workers were exhibiting clinical symptoms of deformed wing virus (DWV). We analysed a pooled bee sample (15 workers) and 29 mites for the presence of Deformed wing virus (DWV), Black queen cell virus (BQCV), Sacbrood virus (SBV), Kashmir bee virus (KBV), Acute bee paralysis virus (ABPV), and Chronic bee paralysis virus (CBPV). Virus positive samples were analysed with a qPCR. Only DWV +RNA was found but with a high titre of up to 108 equivalent virus copies per mite and 106 per bee. Moreover, in all DWV positive mites (N= 12) and in the bee sample virus-RNA was also detected using RT-PCR and tagged RT-PCR, strongly suggesting virus replication. Our data show for the first time that T. mercedesae may be a biological vector of DWV, which would open a novel route of virus spread in A. mellifer

Paenibacillus melissococcoides sp. nov., isolated from a honey bee colony affected by European foulbrood disease.

A novel, facultatively anaerobic, Gram-stain-positive, motile, endospore-forming bacterium of the genus Paenibacillus, designated strain 2.1T, was isolated from a colony of Apis mellifera affected by European foulbrood disease in Switzerland. The rod-shaped cells of strain 2.1T were 2.2–6.5 μm long and 0.7–1.1 μm wide. Colonies of strain 2.1T were orange-pigmented under oxic growth conditions on solid basal medium at 35–37 °C. Strain 2.1T showed catalase and cytochrome c oxidase activity. Its polar lipid profile consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, aminophospholipid and phospholipid. The only respiratory quinone was menaquinone 7, and the major cellular fatty acids were anteiso-C15 : 0, anteiso-C17 : 0, iso-C15 : 0, iso-C17 : 0 and palmitic acid (C16 : 0), which is consistent with other members of the genus Paenibacillus. The G+C content of the genomic DNA of strain 2.1T was 53.3 mol%. Phylogenetic analyses based on the 16S rRNA gene sequence similarity showed that strain 2.1T was closely related to Paenibacillus dendritiformis LMG 21716T (99.7 % similarity) and Paenibacillus thiaminolyticus DSM 7262T (98.8 %). The whole-genome average nucleotide identity between strain 2.1T and the type strains of P. dendritiformis and P. thiaminolyticus was 92 and 91 %, respectively, and thus lower than the 95 % threshold value for delineation of genomic prokaryotic species. Based on the results of phylogenetic, genomic, phenotypic and chemotaxonomic analyses we propose the name Paenibacillus melissococcoides sp. nov. for this novel Paenibacillus species. The type strain is 2.1T (=CCOS 2000T=DSM 113619T=LMG 32539T)

Assembly and analysis of the genome sequence of the yeast Brettanomyces naardenensis CBS 7540

Brettanomyces naardenensis is a spoilage yeast with potential for biotechnological applications for production of innovative beverages with low alcohol content and high attenuation degree. Here, we present the first annotated genome of B. naardenensis CBS 7540. The genome of B. naardenensis CBS 7540 was assembled into 76 contigs, totaling 11,283,072 nucleotides. In total, 5168 protein-coding sequences were annotated. The study provides functional genome annotation, phylogenetic analysis, and discusses genetic determinants behind notable stress tolerance and biotechnological potential of B. naardenensis

Genetic characterization of slow bee paralysis virus of the honeybee (Apis mellifera L.)

Complete genome sequences were determined for two distinct strains of slow bee paralysis virus (SBPV) of honeybees (Apis mellifera). The SBPV genome is approximately 9 5 kb long and contains a single ORF flanked by 5'- and 3'-UTRs and a naturally polyadenylated 3' tail, with a genome organization typical of members of the family Iflaviridae The two strains, labelled `Rothamsted' and 'Harpenden', are 83% identical at the nucleotide level (94% identical at the amino acid level), although this variation is distributed unevenly over the genome. The two strains were found to co-exist at different proportions in two independently propagated SBPV preparations The natural prevalence of SBPV for 847 colonies in 162 apiaries across five European countries was <2%, with positive samples found only in England and Switzerland, in colonies with variable degrees of Varroa infestatio

Predictive Markers of Honey Bee Colony Collapse

Across the Northern hemisphere, managed honey bee colonies, Apis mellifera, are currently affected by abrupt depopulation during winter and many factors are suspected to be involved, either alone or in combination. Parasites and pathogens are considered as principal actors, in particular the ectoparasitic mite Varroa destructor, associated viruses and the microsporidian Nosema ceranae. Here we used long term monitoring of colonies and screening for eleven disease agents and genes involved in bee immunity and physiology to identify predictive markers of honeybee colony losses during winter. The data show that DWV, Nosema ceranae, Varroa destructor and Vitellogenin can be predictive markers for winter colony losses, but their predictive power strongly depends on the season. In particular, the data support that V. destructor is a key player for losses, arguably in line with its specific impact on the health of individual bees and colonies

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

GLADX: An Automated Approach to Analyze the Lineage-Specific Loss and Pseudogenization of Genes

A well-established ancestral gene can usually be found, in one or multiple copies, in different descendant species. Sometimes during the course of evolution, all the representatives of a well-established ancestral gene disappear in specific lineages; such gene losses may occur in the genome by deletion of a DNA fragment or by pseudogenization. The loss of an entire gene family in a given lineage may reflect an important phenomenon, and could be due either to adaptation, or to a relaxation of selection that leads to neutral evolution. Therefore, the lineage-specific gene loss analyses are important to improve the understanding of the evolutionary history of genes and genomes. In order to perform this kind of study from the increasing number of complete genome sequences available, we developed a unique new software module called GLADX in the DAGOBAH framework, based on a comparative genomic approach. The software is able to automatically detect, for all the species of a phylum, the presence/absence of a representative of a well-established ancestral gene, and by systematic steps of re-annotation, confirm losses, detect and analyze pseudogenes and find novel genes. The approach is based on the use of highly reliable gene phylogenies, of protein predictions and on the analysis of genomic mutations. All the evidence associated to evolutionary approach provides accurate information for building an overall view of the evolution of a given gene in a selected phylum. The reliability of GLADX has been successfully tested on a benchmark analysis of 14 reported cases. It is the first tool that is able to fully automatically study the lineage-specific losses and pseudogenizations. GLADX is available at http://ioda.univ-provence.fr/IodaSite/gladx/