Quantitative proteomics reveals divergent responses in Apis mellifera worker and drone pupae to parasitization by Varroa destructor

A blood test to monitor bee health across a European network of agricultural sites of different land-use by MALDI BeeTyping mass spectrometry

There are substantial concerns about impaired honey bee health and colony losses due to several poorly understood factors. We used MALDI profiling (MALDI BeeTyping®) analysis to investigate how some environmental and management factors under field conditions across Europe affected the honey bee haemolymph peptidome (all peptides in the circulatory fluid), as a profile of molecular markers representing the immune status of Apis mellifera. Honey bees were exposed to a range of environmental stressors in 128 agricultural sites across eight European countries in four biogeographic zones, with each country contributing eight sites each for two different cropping systems: oilseed rape (OSR) and apple (APP). The full haemolymph peptide profiles, including the presence and levels of three key immunity markers, namely the antimicrobial peptides (AMPs) Apidaecin, Abaecin and Defensin-1, allowed the honey bee responses to environmental variables to be discriminated by country, crop type and site. When considering just the AMPs, it was not possible to distinguish between countries by the prevalence of each AMP in the samples. However, it was possible to discriminate between countries on the amounts of the AMPs, with the Swedish samples in particular expressing high amounts of all AMPs. A machine learning model was developed to discriminate the haemolymphs of bees from APP and OSR sites. The model was 90.6 % accurate in identifying the crop type from the samples used to build the model. Overall, MALDI BeeTyping® of bee haemolymph represents a promising and cost-effective "blood test" for simultaneously monitoring dozens of peptide markers affected by environmental stressors at the landscape scale, thus providing policymakers with new diagnostic and regulatory tools for monitoring bee health.

Varroa destructor parasitism and Deformed wing virus infection in honey bees are linked to peroxisome-induced pathways

The ectoparasitic mite Varroa destructor transmits and triggers viral infections that have deleterious effects on honey bee colonies worldwide. We performed a manipulative experiment in which worker bees collected at emergence were exposed to Varroa for 72 h, and their proteomes were compared with those of untreated control bees. Label-free quantitative proteomics identified 77 differentially expressed A. mellifera proteins (DEPs). In addition, viral proteins were identified by orthogonal analysis, and most importantly, Deformed wing virus (DWV) was found at high levels/intensity in Varroa-exposed bees. Pathway enrichment analysis suggested that the main pathways affected included peroxisomal metabolism, cyto-/exoskeleton reorganization, and cuticular proteins. Detailed examination of individual DEPs revealed that additional changes in DEPs were associated with peroxisomal function. In addition, the proteome data support the importance of TGF-β signaling in Varroa-DWV interaction and the involvement of the mTORC1 and Hippo pathways. These results suggest that the effect of DWV on bees associated with Varroa feeding results in aberrant autophagy. In particular, autophagy is selectively modulated by peroxisomes, to which the observed proteome changes strongly corresponded. This study complements previous research with different study designs and suggests the importance of the peroxisome, which plays a key role in viral infections.

Omics-based analysis of honey bee (Apis mellifera) response to Varroa sp. parasitisation and associated factors reveals changes impairing winter bee generation

The extensive annual loss of honey bees (Apis mellifera L.) represents a global problem affecting agriculture and biodiversity. The parasitic mite Varroa destructor, associated with viral co-infections, plays a key role in this loss. Despite years of intensive research, the complex mechanisms of Varroa - honey bee interaction are still not fully defined. Therefore, this study employed a unique combination of transcriptomic, proteomic, metabolomic, and functional analyses to reveal new details about the effect of Varroa mites and naturally associated factors, including viruses, on honey bees. We focused on the differences between Varroa parasitised and unparasitised ten-day-old worker bees collected before overwintering from the same set of colonies reared without anti-mite treatment. Supplementary comparison to honey bees collected from colonies with standard anti-Varroa treatment can provide further insights into the effect of a pyrethroid flumethrin. Analysis of the honey bees exposed to mite parasitisation revealed alterations in the transcriptome and proteome related to immunity, oxidative stress, olfactory recognition, metabolism of sphingolipids, and RNA regulatory mechanisms. The immune response and sphingolipid metabolism were strongly activated, whereas olfactory recognition and oxidative stress pathways were inhibited in Varroa parasitised honey bees compared to unparasitised ones. Moreover, metabolomic analysis confirmed the depletion of nutrients and energy stores, resulting in a generally disrupted metabolism in the parasitised workers. The combined omics-based analysis conducted on strictly parasitised bees revealed the key molecular components and mechanisms underlying the detrimental effects of Varroa sp. and its associated pathogens. This study provides the theoretical basis and interlinked datasets for further research on honey bee response to biological threats and the development of efficient control strategies against Varroa mites.

Varroa destructor: how does it harm Apis mellifera honey bees and what can be done about it?

Since its migration from the Asian honey bee (Apis cerana) to the European honey bee (Apis mellifera), the ectoparasitic mite Varroa destructor has emerged as a major issue for beekeeping worldwide. Due to a short history of coevolution, the host–parasite relationship between A. mellifera and V. destructor is unbalanced, with honey bees suffering infestation effects at the individual, colony and population levels. Several control solutions have been developed to tackle the colony and production losses due to Varroa, but the burden caused by the mite in combination with other biotic and abiotic factors continues to increase, weakening the beekeeping industry. In this synthetic review, we highlight the main advances made between 2015 and 2020 on V. destructor biology and its impact on the health of the honey bee, A. mellifera. We also describe the main control solutions that are currently available to fight the mite and place a special focus on new methodological developments, which point to integrated pest management strategies for the control of Varroa in honey bee colonies.

Selection and hybridization shaped the rapid spread of African honey bee ancestry in the Americas

Recent biological invasions offer 'natural' laboratories to understand the genetics and ecology of adaptation, hybridization, and range limits. One of the most impressive and well-documented biological invasions of the 20th century began in 1957 when Apis mellifera scutellata honey bees swarmed out of managed experimental colonies in Brazil. This newly-imported subspecies, native to southern and eastern Africa, both hybridized with and out-competed previously-introduced European honey bee subspecies. Populations of scutellata-European hybrid honey bees rapidly expanded and spread across much of the Americas in less than 50 years. We use broad geographic sampling and whole genome sequencing of over 300 bees to map the distribution of scutellata ancestry where the northern and southern invasions have presently stalled, forming replicated hybrid zones with European bee populations in California and Argentina. California is much farther from Brazil, yet these hybrid zones occur at very similar latitudes, consistent with the invasion having reached a climate barrier. At these range limits, we observe genome-wide clines for scutellata ancestry, and parallel clines for wing length that span hundreds of kilometers, supporting a smooth transition from climates favoring scutellata-European hybrid bees to climates where they cannot survive winter. We find no large effect loci maintaining exceptionally steep ancestry transitions. Instead, we find most individual loci have concordant ancestry clines across South America, with a build-up of somewhat steeper clines in regions of the genome with low recombination rates, consistent with many loci of small effect contributing to climate-associated fitness trade-offs. Additionally, we find no substantial reductions in genetic diversity associated with rapid expansions nor complete dropout of scutellata ancestry at any individual loci on either continent, which suggests that the competitive fitness advantage of scutellata ancestry at lower latitudes has a polygenic basis and that scutellata-European hybrid bees maintained large population sizes during their invasion. To test for parallel selection across continents, we develop a null model that accounts for drift in ancestry frequencies during the rapid expansion. We identify several peaks within a larger genomic region where selection has pushed scutellata ancestry to high frequency hundreds of kilometers past the present cline centers in both North and South America and that may underlie high-fitness traits driving the invasion.

Varroa destructor parasitism has a greater effect on proteome changes than the deformed wing virus and activates TGF-β signaling pathways

Honeybee workers undergo metamorphosis in capped cells for approximately 13 days before adult emergence. During the same period, Varroa mites prick the defenseless host many times. We sought to identify proteome differences between emerging Varroa-parasitized and parasite-free honeybees showing the presence or absence of clinical signs of deformed wing virus (DWV) in the capped cells. A label-free proteomic analysis utilizing nanoLC coupled with an Orbitrap Fusion Tribrid mass spectrometer provided a quantitative comparison of 2316 protein hits. Redundancy analysis (RDA) showed that the combination of Varroa parasitism and DWV clinical signs caused proteome changes that occurred in the same direction as those of Varroa alone and were approximately two-fold higher. Furthermore, proteome changes associated with DWV signs alone were positioned above Varroa in the RDA. Multiple markers indicate that Varroa activates TGF-β-induced pathways to suppress wound healing and the immune response and that the collective action of stressors intensifies these effects. Furthermore, we indicate JAK/STAT hyperactivation, p53-BCL-6 feedback loop disruption, Wnt pathway activation, Wnt/Hippo crosstalk disruption, and NF-κB and JAK/STAT signaling conflict in the Varroa–honeybee–DWV interaction. These results illustrate the higher effect of Varroa than of DWV at the time of emergence. Markers for future research are provided.

Hypoxic conditioned medium derived from bone marrow mesenchymal stromal cells protects against ischemic stroke in rats

In recent years, studies have shown that the secretome of bone marrow mesenchymal stromal cells (BMSCs) contains many growth factors, cytokines, and antioxidants, which may provide novel approaches to treat ischemic diseases. Furthermore, the secretome may be modulated by hypoxic preconditioning. We hypothesized that conditioned medium (CM) derived from BMSCs plays a crucial role in reducing tissue damage and improving neurological recovery after ischemic stroke and that hypoxic preconditioning of BMSCs robustly improves these activities. Rats were subjected to ischemic stroke by middle cerebral artery occlusion and then intravenously administered hypoxic CM, normoxic CM, or Dulbecco modified Eagle medium (DMEM, control). Cytokine antibody arrays and label-free quantitative proteomics analysis were used to compare the differences between hypoxic CM and normoxic CM. Injection of normoxic CM significantly reduced the infarct area and improved neurological recovery after stroke compared with administering DMEM. These outcomes may be associated with the attenuation of apoptosis and promotion of angiogenesis. Hypoxic preconditioning significantly enhanced these therapeutic effects. Fourteen proteins were significantly increased in hypoxic CM compared with normoxic CM as measured by cytokine arrays. The label-free quantitative proteomics analysis revealed 163 proteins that were differentially expressed between the two groups, including 107 upregulated proteins and 56 downregulated proteins. Collectively, our results demonstrate that hypoxic CM protected brain tissue from ischemic injury and promoted functional recovery after stroke in rats and that hypoxic CM may be the basis of a potential therapy for stroke patients.