# Honey Bee Larval Hemolymph as a Source of Key Nutrients and Proteins Offers a Promising Medium for Varroa destructor Artificial Rearing

Vincent Piou, Caroline Vilarem, Solène Blanchard, Jean-Marc Strub, Fabrice Bertile, Michel Bocquet, Karim Arafah, Philippe Bulet, Angélique Vétillard

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Vincent Piou, Caroline Vilarem, Solène Blanchard, Jean-Marc Strub, Fabrice Bertile, et al.. Honey Bee Larval Hemolymph as a Source of Key Nutrients and Proteins Offers a Promising Medium for Varroa destructor Artificial Rearing. International Journal of Molecular Sciences, 2023, 24 (15), pp.1-19. 10.3390/ijms241512443 . hal-04186784

HAL Id: hal-04186784

https://cnam.hal.science/hal-04186784v1

Submitted on 24 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Article

# Honey Bee Larval Hemolymph as a Source of Key Nutrients and Proteins Offers a Promising Medium for *Varroa destructor* Artificial Rearing

Vincent Piou<fnref n="1" />, Caroline Vilarem<fnref n="1" />,<fnref n="2" /> ID, Solène Blanchard<fnref n="1" />, Jean-Marc Strub<fnref n="3" /> ID, Fabrice Bertile<fnref n="3" /> ID,
Michel Bocquet<fnref n="4" />, Karim Arafah<fnref n="5" />, Philippe Bulet<fnref n="5" />,<fnref n="6" /> ID and Angélique Vétillard<fnref n="1" />,<fnref n="7" />,*

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Piou, V.; Vilarem, C.; Blanchard, S.; Strub, J.-M.; Bertile, F.; Bocquet, M.; Arafah, K.; Bulet, P.; Vétillard, A. Honey Bee Larval Hemolymph as a Source of Key Nutrients and Proteins Offers a Promising Medium for *Varroa destructor* Artificial Rearing. *Int. J. Mol. Sci.* **2023**, *24*, 12443.
https://doi.org/10.3390/ijms241512443

Academic Editors: Michael Stear and
Katarzyna Donskow-Lysoniewska

Received: 12 July 2023
Revised: 28 July 2023
Accepted: 29 July 2023
Published: 4 August 2023

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** *Varroa destructor*, a major ectoparasite of the Western honey bee *Apis mellifera*, is a widespread pest that damages colonies in the Northern Hemisphere. Throughout their lifecycle, *V. destructor* females feed on almost every developmental stage of their host, from the last larval instar to the adult. The parasite is thought to feed on hemolymph and fat body, although its exact diet and nutritional requirements are poorly known. Using artificial Parafilm™ dummies, we explored the nutrition of *V. destructor* females and assessed their survival when fed on hemolymph from bee larvae, pupae, or adults. We compared the results with mites fed on synthetic solutions or filtered larval hemolymph. The results showed that the parasites could survive for several days or weeks on different diets. Bee larval hemolymph yielded the highest survival rates, and filtered larval plasma was sufficient to maintain the mites for 14 days or more. This cell-free solution therefore theoretically contains all the necessary nutrients for mite survival. Because some bee proteins are known to be hijacked without being digested by the parasite, we decided to run a proteomic analysis of larval honey bee plasma to highlight the most common proteins in our samples. A list of 54 proteins was compiled, including several energy metabolism proteins such as Vitellogenin, Hexamerin, or Transferrins. These molecules represent key nutrient candidates that could be crucial for *V. destructor* survival.

**Keywords:** mite; artificial dummy; hemolymph; parasite; protein; nutrition; insect physiology

## 1. Introduction

Each organism comes with its load of parasites and pathogens, either as a result of a long coevolution process or because of new contacts between two species [1,2]. These parasites can have deleterious impacts on ecosystems or animal and human health [3–5]. Most of the time, parasitic species are studied in the context of a two- or even three-way relationship with their host and with the pathogens they transmit [6,7]. Although challenging, studying such harmful species separately in laboratory conditions represents a complementary and more controlled approach [8,9]. The difficulty in focusing on parasites alone is to isolate them and keep them away from their host. Artificial feeding systems are often necessary to successfully maintain and study parasites separately over prolonged periods [10]. Many protocols exist for ticks or hematophagous insects and have enlightened our understanding of parasite disease transmission and susceptibility to treatments [11–14].

<footnote n="1">¹ Laboratoire Evolution et Diversité Biologique, UMR5174, CNRS-Université de Toulouse III-IRD—Université Paul Sabatier, 31077 Toulouse, France; vincent.piou@univ-tlse3.fr (V.P.); solene.blanchard@univ-tlse3.fr (S.B.)</footnote>

<footnote n="2">² M2i Biocontrol-Entreprise SAS, 46140 Parnac, France</footnote>

<footnote n="3">³ Laboratoire de Spectrométrie de Masse Bio-Organique, Département des Sciences Analytiques, Institut Pluridisciplinaire Hubert Curien, UMR 7178 (CNRS-UdS), 67037 Strasbourg, France; fbertile@unistra.fr (F.B.)</footnote>

<footnote n="4">⁴ Apimedia, 74370 Pringy, Annecy, France</footnote>

<footnote n="5">⁵ Plateforme BioPark d'Archamps, 74160 Archamps, France</footnote>

<footnote n="6">⁶ Institute pour l'Avancée des Biosciences, CR Université Grenoble Alpes, Inserm U1209, CNRS UMR 5309, 38000 Grenoble, France</footnote>

<footnote n="7">⁷ Conservatoire National des Arts et Métiers (CNAM), Unité Métabiot, 22440 Ploufragan, France</footnote>

<footnote>* Correspondence: angelique.vetillard@lecnam.net</footnote>
Since standardized laboratory media and rearing environments can be artificially modified, these methods are also unrivaled for exploring the feeding biology of parasitic pests.

Because of its impact on the Western honey bee (Apis mellifera), beekeeping, and agriculture, the ectoparasite *Varroa destructor* is one of the main pests studied throughout the world [15–17]. Since it is an obligate parasite, it is only possible to isolate the mite from its host and to maintain it in the laboratory for limited periods. Until now, most rearing methods have relied on easy-to-make Parafilm™ membranes but on complex synthetic feeding solutions [18–20]. Recently, the importance of 3D round-shaped substrates in triggering mite feeding behavior under artificial conditions was highlighted [21]. The use of a phosphate-buffered saline (PBS) solution also offers a nutrient-deprived control solution that can be employed as a basic medium to be supplemented in *V. destructor* studies [21]. Although the exact composition of the natural, balanced diet remains unknown, testing solutions prepared from honey bee tissue samples can yield better survival rates [22]. Thus, they could represent a reliable alternative to synthetic solutions for laboratory rearing. Altogether, these artificial feeding protocols offer promising prospects for exploring the physiology, feeding behavior, and nutritional requirements of *V. destructor*.

Unlike many other bee parasites, *V. destructor* infests both adult and immature stages of its host throughout its cycle [16,17]. Female mites have to puncture the cuticle of the larval, pupal, or adult host to get nutrients from hemolymph and fat body [22–25]. However, we do not know whether the necessary nutrients come mostly from lysed hemocytes or fat body cells or directly from molecules circulating in the hemolymph. The precise nature of the proteins, carbohydrates, vitamins, or lipids essential to the mite's survival and reproduction has never been determined either. In addition to being a potential source of amino acids necessary for protein synthesis, undigested honey bee proteins could act as key nutrients in the mite's physiology, as they have been found inside *V. destructor* hemocoel and reproductive organs. The hijacking of such molecules could be of prime importance for parasite physiology [24,26]. Determining the nutritional requirements of *V. destructor*, especially in terms of proteins, is thus crucial to better understanding the parasite and rearing it under artificial conditions.

Following improvements in rearing protocols, we experimentally explored the nu-
trition and feeding requirements of mites under laboratory conditions. In this study, we
focused on hemolymph, as it is usually easier to collect in great quantities compared to fat
body. Hemolymph from three different stages of honey bee development was tested to
assess its nutritional potential as a natural food source. Larval, pupal, or adult hemolymph
was included in 3D rounded Parafilm™ dummies to assess the consumption of food and
the survival of mites. The idea was to investigate if the nutritional value of hemolymph
could vary throughout honey bee development and impact the mite. To further assess
whether essential nutrients came from lysed cells or from molecules circulating in the
hemolymph, we tested treated larval hemolymph or synthetic media as feeding solutions.
Finally, the protein content of filtered cell-free larval hemolymph—used successfully in
our rearing bioassays—was analyzed by off-gel proteomics in order to identify protein
candidates that could be key nutrients for *V. destructor*.

2. Results

2.1. *Bioassay 1: Effect of Hemolymph from Different Bee Developmental Stages on V. destructor Survival*

Since untreated hemolymph can quickly become unsuitable for *V. destructor* feeding
after a few days at 34 °C (Figure S1), freshly collected larval, pupal, or adult hemolymph
was heated at 65 °C for 7 min (as in [20]). This process prevents contamination and limits
clotting or melanization once hemolymph is included into feeding Parafilm™ dummies
 displays in Figure 1).
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Artificial feeding chambers used to maintain mature female mites. Chambers consist of Petri dishes (Ø 5 cm) covered with Parafilm™. Two dummies containing an artificial feeding solution were placed in the chamber before 10 to 12 females *V. destructor* were inserted.</figcaption></figure>

In our artificial conditions, survival reached more than 70% of still alive mites at the end of the first week (day 7) in each of the three rearing conditions (Figure 2). The parasite was significantly more likely to stay alive for a longer period when fed on larval hemolymph compared to pupal or adult hemolymph (log rank test df = 2, $\chi^2$ = 42.9, *p* < 0.001). The median survival even reached 23 days in the case of larval hemolymph. Pupal and adult hemolymph solutions do not result in significant survival differences, and the median survival equals 11 days (Figure 2).

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Survival curves and confidence intervals of mites fed on hemolymph of larvae (Hemlar, red curve), pupae (Hempu, light blue curve), or adult bees (Hemadu, green curve). The dashed black line shows the median survival for the three groups. Colored dashed lines surrounding the survival curves represent the 95% confidence interval. *p*-value resulting from the log rank test is shown in the top right corner of the graph along with pairwise *p*-values comparing the different feeding groups.</figcaption></figure>

## 2.2. *Bioassay 2: Feeding on Treated Larval Hemolymph Has No Detrimental Effect on V. destructor Survival*

After ensuring that storage at -20 °C did not have any impact on parasites' survival (Figure S2), larval hemolymph was treated differently and inserted in Parafilm™ dummies to feed female mites (Bioassay 2 in Table 1). More specifically, filtration at 0.2 µm led to a cell-free hemolymph preparation, while heating at 65 °C for 7 min limited melanization [20] and released the cell contents into the medium. Eukaryotic cells exposed to high temperatures (above 55 °C) indeed suffer instant necrosis after membrane rupture [27,28]. Solutions were prepared as described in Table 1 and stored at -20 °C.
<table><caption>**Table 1.** Details of the conditions tested in three different bioassays. Bee developmental stage and hemolymph treatment (heating, filtration, storage) are shown in the table, along with abbreviations used in the graphs and sample size N. In rare cases, there is a slight difference between the initial sample size (N₀) and the sample sizes after one or two weeks because some mites escaped from their chambers in the course of our experiment (Figure S2).</caption><thead><tr><th>Bioassay—<br>Factor<br>Tested</th><th>Condition</th><th>Treatment Applied to Feeding Solution<br>Heated<br>(65 °C, 7 min)</th><th>Filtered<br>(0.2 µm)</th><th>Storage<br>(-20 °C)</th><th>Abbreviation</th><th>Number of<br>Mites Tested<br>(N)</th></tr></thead><tbody><tr><td rowspan="3">1—Origin of<br>hemolymph</td><td>Larva</td><td>Yes</td><td>No</td><td>No</td><td>Hemlar</td><td>64</td></tr><tr><td>Pupa</td><td>Yes</td><td>No</td><td>No</td><td>HemPu</td><td>62</td></tr><tr><td>Adult</td><td>Yes</td><td>No</td><td>No</td><td>HemAdu</td><td>39</td></tr><tr><td rowspan="4">2—Treatment<br>of larval<br>hemolymph</td><td>Stored</td><td>Yes</td><td>No</td><td>Yes</td><td>Stored-<br>Hemlar</td><td>90</td></tr><tr><td>Filtered<br>Heated then<br>filtered</td><td>No</td><td>Yes</td><td>Yes</td><td>Hemlar-F</td><td>44</td></tr><tr><td>Filtered then<br>heated</td><td>Yes (firstly)</td><td>Yes (secondly)</td><td>Yes</td><td>Hemlar-HF</td><td>57</td></tr><tr><td>Yes (secondly)</td><td>Yes (firstly)</td><td>Yes</td><td>Yes</td><td>Hemlar-FH</td><td>54</td></tr><tr><td rowspan="4">3—Synthetic<br>diet</td><td>Negative</td><td>No</td><td>No</td><td>No</td><td>Negative<br>control</td><td>24</td></tr><tr><td>PBS</td><td>No</td><td>No</td><td>No</td><td>PBS control</td><td>31</td></tr><tr><td>PBS + Sugar<br>(100 mg/mL)</td><td>No</td><td>Yes</td><td>No</td><td>Sugar</td><td>92</td></tr><tr><td>PBS + Sugar<br>(100 mg/mL)<br>+ Yeast extract<br>(40 mg/mL)</td><td>No</td><td>Yes</td><td>No</td><td>SugarYE</td><td>31</td></tr></tbody></table>

In these conditions, mites survived for seven days in more than 77% of cases regardless of the hemolymph treatment. Although slight differences can be observed in Figure 3, no significance was reached between the four conditions (GLM, df = 3, $\chi^2$ = 2.06, *p* = 0.56). On day 14, the proportion of alive mites remained above 50%. More precisely, it ranged from 51.1% [CI95: 35.8–66.3] for heated cell-free hemolymph (Hemlar-HF) to 70.9% [CI95: 60.1–80.2] in the case of only heated hemolymph (Stored-Hemlar). Again, no significant difference between the four conditions was noticed at day 14 (GLM df = 3, $\chi^2$ = 5.98, *p* = 0.11).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Survival of mites fed on larval hemolymph under different treatments (in % ± CI95). Prior to its insertion in Parafilm™ dummies, hemolymph was heated (Stored-Hemlar), filtered (Hemlar-F), or both heated and filtered (Hemlar-FH for the filtered then heated condition and Hemlar-HF for the heated then filtered condition). No significant difference was highlighted either on day 7 (dark blue) or on day 14 (light blue).</figcaption></figure>
## 2.3. Bioassay 3: Artificial Feeding with Synthetic Diet Impacts V. destructor Survival

To mimic hemolymph, PBS solutions supplemented with sugars or sugars and yeast extracts were tested (Table 1). Both glucose and fructose are major hemolymph carbohydrates detected at a concentration of approximately 5 mg/mL each in adult hemolymph, in addition to the glucose disaccharide trehalose (30 mg/mL) [29]. The concentrations used here (50 mg/mL), higher than those observed in natural adult bee hemolymph, were chosen based on preliminary results (Figure S3). In our second feeding condition, yeast extract was added to the carbohydrate solution to simulate amino acid and vitamin uptakes [30]. Again, concentrations of yeast extracts (40 mg/mL) were chosen to be included within a suitable range determined during preliminary experiments (Figure S3). These two test conditions were compared to starved mites, PBS-fed or hemolymph-fed mites.

The first step was to ensure that mites were fed in such artificial conditions. By coloring the feeding media, we managed to verify that after 24 h, the feeding rate was over 95% in every condition (Figures 4 and 5). No significant difference was detected between bee larval hemolymph and sugar-containing PBS (Biased reduced GLM: df = 3, $\chi^2$ = 0.26, *p* = 0.97).

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Pictures showing the coloration of *V. destructor* guts through the cuticle (upper panels) or after dissection (lower panels). (A,D): Mites that feed on a natural uncolored diet; (B,E): Mites that feed on a synthetic Blue-FCF colored diet; (C,F): Mite that feed on Blue-FCF colored larval hemolymph.</figcaption></figure>

Survival was highly impacted by the type of artificial solution used to feed *V. destructor* (log rank test df = 4, $\chi^2$ = 257, *p* < 0.001). Mortality was the highest for starved mites (negative control), and a slight tendency to reduce mortality was associated with the presence of PBS-filled dummies. The supplementation with glucose and fructose (Sugar) was sufficient to improve the survival of artificially fed mites compared to sugar-free PBS-fed or starved mites (PBS control and negative control, respectively; Figure 6). The median survival was around five days for mites fed with PBS supplemented with sugar or six days when yeast extract was added to the sugar supplementation (Figure 6; Sugar and SugarYE, respectively). The small increase associated with the addition of yeast extracts was significant. However, whatever the synthetic solution tested, survival was always significantly lower than in the case of hemolymph-fed mites (Figure 6).
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Feeding success of mites fed on synthetic solutions in comparison with mites fed on larval hemolymph (Stored-Hemlar). Synthetic solutions consisted of PBS (control group), PBS supplemented with glucose and fructose (Sugar), or with glucose, fructose, and yeast extract (SugarYE). The coloration of the gut was checked 24 h after the beginning of rearing. No significant difference was detected between the four groups (Larval hemolymph N = 30; PBS control N = 31; Sugar N = 92; SugarYE N = 31).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Survival curves and confidence intervals of mites fed synthetic solutions or larval hemolymph. Negative control (yellow curve) = starved mites without feeding dummies; PBS Control (blue curve) = mites with PBS-filled Parafilm™ dummies; Sugar (red curve) = mites fed on a sugar-supplemented PBS solution included in dummies; SugarYE (green curve) = mites fed on sugar and yeast extract supplemented PBS solution included in dummies. Colored dashed lines surrounding the survival curves represent the 95% confidence interval. The *p*-value resulting from the log rank test is shown in the bottom right corner of the graph along with pairwise *p*-values comparing the different feeding groups.</figcaption></figure>
## 2.4. Protein Content of Filtered Hemolymph

The protein content of filtered larval hemolymph (similar to the Hemlar-F condition in Table 1) was analyzed by off-gel proteomics. Four different larval samples collected from three different colonies were used (Table S1). A total of 1118 proteins were detected, matching mostly insect proteins (93.5%) and, more precisely, Hymenoptera species (91.1%, Figure 7A). The small proportion remaining (6.5%) can be imputed to the presence of *V. destructor* and bee pathogens such as *Nosema* spp. or viruses, with 5.2% of proteins matching Acari species and 1.3% of proteins matching mostly unicellular species. Among the 1045 insect proteins, we focused on a set of 88 proteins that were consistently detected in our samples from three different colonies (Figure 7B). Once uncharacterized proteins and redundant identifications were processed, a shortlist of 54 proteins remained. This set likely consists of the most frequent and abundant proteins naturally present in bee larval hemolymph (Table 2). The 54 proteins could thus be relevant candidate nutrients required for *V. destructor* physiology. Immunity and development are the most represented functions, with 39% of the identified filtered hemolymph proteins being involved in each of these functions. Besides tegument and extracellular matrix molecules, lipid transport and energy metabolism proteins such as Vitellogenin, Apolipophorin, Hexamerin, and Transferrin are also well represented (8 of the 54 proteins; Table 2).

<figure><img src="image_12.png" /><figcaption>**Figure 7.** (A) Number of proteins detected in filtered larval hemolymph and systematic origin of the matching protein sequence. (B) Number of identified proteins that were found only once or shared between two, three, or four hemolymph samples from our three different colonies.</figcaption></figure>

<table><caption>**Table 2.** List of proteins detected in each of the four filtered larval hemolymph samples from three different colonies. The protein name and the species from the matching protein sequence are indicated, along with the consistency of detection and the probable function described in the literature.</caption><thead><tr><th>Protein Name</th><th>Species</th><th>Detected in All or Part of Technical Replicates</th><th>Function</th><th>Reference</th></tr></thead><tbody><tr><td>A-agglutinin anchorage subunit isoform X1</td><td>Apis mellifera</td><td>All</td><td>Cell-cell adhesion; immunity</td><td>Rowley and Ratcliffe (1980) [31]</td></tr><tr><td>Abaecin</td><td>Apis cerana</td><td>11/12</td><td>Immunity</td><td>Casteels et al., (1990) [32]; Plua and sokol (2020) [33]</td></tr><tr><td>AF4/FMR2 family member 4</td><td>Apis florea/Nomia melanderi</td><td>7/12</td><td>Cell growth/identity; development</td><td>Wittwer et al., (2001) [34]</td></tr><tr><td>Apidaecin</td><td>Apis cerana</td><td>All</td><td>Immunity</td><td>Plua and sokol (2020) [33]</td></tr><tr><td>Apolipophorins</td><td>Apis mellifera</td><td>All</td><td>Lipid metabolism; energy storage; immunity</td><td>Cabri et al., (2018) [35]; Kim and Jin (2014) [36]</td></tr><tr><td>Apolipoprotein D</td><td>Apis mellifera/Apis cerana</td><td>All</td><td>Lipid metabolism; energy storage</td><td>Chan and Foster (2008) [37]</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Protein Name</th><th>Species</th><th>Detected in All or Part of Technical Replicates</th><th>Function</th><th>Reference</th></tr></thead><tbody><tr><td>Beta-galactosidase</td><td>Apis mellifera</td><td>All</td><td>Carbohydrates<br>metabolism</td><td>Peng (1980) [38];<br>Ricigliano et al., (2017) [39]</td></tr><tr><td>Bifunctional methylenetetrahydrofolate<br>dehydrogenase/cyclohydrolase,<br>mitochondrial<br>BMP and activin<br>membrane-bound<br>inhibitor</td><td>Bombus vancouverensis</td><td>8/12</td><td>Oxidative stress;<br>development; aging</td><td>Tremblay et al., (1995) [40];<br>Yu et al., (2015) [41]</td></tr><tr><td>chitinase-like protein<br>EN03 isoform X1</td><td>Apis cerana</td><td>All</td><td>Development</td><td>Upadhyay et al.,<br>(2017) [42]</td></tr><tr><td>Chymotrypsin inhibitor</td><td>Apis mellifera</td><td>All</td><td>Development</td><td>Li et al., (2010) [43]<br>Kim et al., (2013) [44];<br>Corral-Rodriguez et al.,<br>(2009) [45]; Rhoads et al.,<br>(2000) [46]</td></tr><tr><td>Class E basic<br>helix-loop-helix protein 22</td><td>Apis mellifera</td><td>10/12</td><td>Development<br>(transcription factor)</td><td>Ledent and Vervoort<br>(2001) [47]; Wan et al.,<br>(2016) [48]; Wang et al.,<br>(2008) [49]</td></tr><tr><td>Collagen alpha-1(IV) chain</td><td>Apis mellifera</td><td>All</td><td>Extracellular matrix;<br>development</td><td>Pastor-Pareja et al.,<br>(2011) [50]; Sutherland<br>et al., (2013) [51]</td></tr><tr><td>Cubilin</td><td>Apis cerana</td><td>11/12</td><td>Endocytosis</td><td>Zhang et al., (2013) [52]</td></tr><tr><td>Cuticular protein 2<br>precursor</td><td>Apis mellifera</td><td>All</td><td>Tegument protein;<br>development (preparation<br>of ecdysis)</td><td>Hopkins et al., (2000) [53];<br>Soares et al., (2007) [54]</td></tr><tr><td>Endocuticle structural<br>glycoprotein SgAbd-1</td><td>Pseudomyrmex gracilis</td><td>11/12</td><td>Tegument protein</td><td>Micas et al., (2016) [55]</td></tr><tr><td>FABP-like protein</td><td>Apis cerana</td><td>All</td><td>Lipid transport/uptake;<br>immunity; development</td><td>Caccia et al., (2012) [56];<br>Chen et al., (2022) [57];<br>Cheng et al., (2013) [58]</td></tr><tr><td>Fibrillin-2</td><td>Apis mellifera</td><td>10/12</td><td>Extracellular matrix</td><td>Piha-Gossack et al.,<br>(2012) [59]</td></tr><tr><td>Fibroin heavy chain<br>isoform X1</td><td>Apis mellifera</td><td>9/12</td><td>Silk protein</td><td>Sutherland et al.,<br>(2006) [60]</td></tr><tr><td>Flexible cuticle protein<br>12-like isoform X1</td><td>Apis florea</td><td>All</td><td>Tegument protein;<br>development</td><td>Rebers et al., (1988) [61]</td></tr><tr><td>Floculation protein FLO11<br>isoform X1</td><td>Apis mellifera</td><td>All</td><td>Extracellular matrix;<br>development</td><td>Zhao et al., (2020) [62]</td></tr><tr><td>Glucose dehydrogenase<br>[FAD, quinone]</td><td>Apis mellifera</td><td>All</td><td>Development; immunity</td><td>Cox-Foster et al., (1990,<br>1994) [63,64]</td></tr><tr><td>Glycine-rich cell wall<br>structural protein 1<br>isoform X1</td><td>Apis mellifera</td><td>All</td><td>Cuticle protein;<br>development; response to<br>stress</td><td>Zhang et al., (2008) [65];<br>Zhong et al., (2005) [66]</td></tr><tr><td>Hexamerin</td><td>Apis mellifera</td><td>All</td><td>Storage protein</td><td>Martin et al., (2010) [67]</td></tr><tr><td>Inactive serine protease<br>scarface</td><td>Apis mellifera</td><td>11/12</td><td>Development</td><td>Contreras et al., (2021) [68]</td></tr><tr><td>Larval-specific very high<br>density lipoprotein<br>precursor</td><td>Apis mellifera</td><td>All</td><td>Lipid transport/storage</td><td>Shipman et al., (1987) [69]</td></tr><tr><td>Leukocyte elastase<br>inhibitor</td><td>Apis mellifera</td><td>All</td><td>Coagulation/immunity</td><td>Kim et al., (2013) [44];<br>Corral-Rodriguez et al.,<br>(2009) [45]; Rhoads et al.,<br>(2000) [46]</td></tr><tr><td>Interferon-related<br>developmental regulator<br>1-like</td><td>Apis mellifera</td><td>All</td><td>Development; immunity</td><td>Arockiaraj et al.,<br>(2014) [70]; Hoffmann<br>et al., (1996) [71]; Stanifer<br>et al., (2019) [72]</td></tr><tr><td>Lysozyme</td><td>Apis mellifera</td><td>All</td><td>Immunity</td><td>Al-Ghamdi et al.,<br>(2021) [73]</td></tr><tr><td>Neurofilament heavy<br>polypeptide</td><td>Apis mellifera</td><td>All</td><td>Neural cytoskeleton;<br>neural development</td><td>Bezabih et al., (2017) [74];<br>Petzold (2005) [75]</td></tr><tr><td>OBP13</td><td>Apis mellifera</td><td>All</td><td>Olfaction; transport<br>protein in larvae</td><td>Forêt and Maleszka<br>(2006) [76]</td></tr><tr><td>Olorant binding protein<br>14 precursor</td><td>Apis mellifera</td><td>All</td><td>Olfaction; transport<br>protein in larvae</td><td>Forêt and Maleszka<br>(2006) [76]</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Protein Name</th><th>Species</th><th>Detected in All or Part of Technical Replicates</th><th>Function</th><th>Reference</th></tr></thead><tbody><tr><td>Omega-conotoxin-like protein 1</td><td>Apis mellifera</td><td>All</td><td>Immunity/melanization</td><td>Bloch et Cohen (2014) [77]</td></tr><tr><td>Peptidyl-prolyl cis-trans isomerase B precursor</td><td>Apis mellifera</td><td>All</td><td>Protein folding; development, cell differentiation; oxidative stress and immunity</td><td>Wang and Heitmann (2005) [78]; Yoon et al., (2022) [79]</td></tr><tr><td>Peritrophin-1</td><td>Apis mellifera</td><td>All</td><td>Tegument protein (of the peritrophic membrane)</td><td>Park et al., (2016) [80]</td></tr><tr><td>Phenoloxidase subunit A3</td><td>Apis mellifera</td><td>All</td><td>Melanization/immunity</td><td>Wilson-Rich et al., (2008) [81]</td></tr><tr><td>Phenoloxidase-activating factor 2 isoform X1</td><td>Apis mellifera</td><td>All</td><td>Melanization/immunity</td><td>Wilson-Rich et al., (2008) [81]</td></tr><tr><td>Plexin domain-containing protein 2 isoform X1</td><td>Apis florea</td><td>All</td><td>Clotting; development; immunity</td><td>Miller-Delaney et al., (2011) [82]; Thibord et al., (2019) [83]</td></tr><tr><td>Odorant receptor 43a-like isoform X1</td><td>Vollenhovia emeryi</td><td>5/12</td><td>Olfaction</td><td>Liu et al., (2020) [84]</td></tr><tr><td>Prisilkin-39 isoform X1 Probable G-protein coupled receptor Mth-like 10</td><td>Apis cerana</td><td>8/12</td><td>Silk matrix protein</td><td>Jung et al., (2021) [85]</td></tr><tr><td>Protein D2</td><td>Pseudomyrmex gracilis</td><td>6/12</td><td>Oxidative stress and longevity</td><td>Liu et al., (2021) [86]</td></tr><tr><td>Protein mesh isoform X1</td><td>Apis mellifera</td><td>All</td><td>Brain development</td><td>Jørgensen (1983) [87]</td></tr><tr><td>Putative acyl-CoA-binding protein Putative cyclin-dependent serine/threonine-protein kinase</td><td>Apis mellifera</td><td>10/12</td><td>Cell-cell adhesion; development</td><td>Jonusaite et al., (2020) [88]</td></tr><tr><td>Secapin-2 precursor</td><td>Apis mellifera</td><td>All</td><td>Lipid metabolism</td><td>Majerowicz et al., (2016) [89]</td></tr><tr><td>Secapin-3 precursor</td><td>Apis mellifera</td><td>All</td><td>Regulation of transcription; cell division</td><td>Zhao et al., (2018) [90]</td></tr><tr><td>Serine protease inhibitor 3 Thymosin beta-a</td><td>Apis mellifera</td><td>All</td><td>Immunity; venom</td><td>Al-Naggar et al., (2023) [91]; Doublet et al., (2017)</td></tr><tr><td>Titin homolog, partial</td><td>Apis cerana</td><td>11/12</td><td>Immunity; development</td><td>Al-Naggar et al., (2023) [91]; Doublet et al., (2017) [92]</td></tr><tr><td>Transcriptional regulator ATRX homolog</td><td>Bombus terrestris</td><td>10/12</td><td>Muscle protein</td><td>Shakeel et al., (2019) [93]</td></tr><tr><td>Transferrin</td><td>Apis mellifera</td><td>All</td><td>Gene regulation; cell division</td><td>Zhang et al., (2011) [94]</td></tr><tr><td>Venom carbohydrate-rich protein precursor</td><td>Apis mellifera</td><td>All</td><td>Iron transport; immunity; energy metabolism</td><td>Hooper and Thuma (2005) [95]</td></tr><tr><td>Vitellogenin-6-like</td><td>Apis dorsata</td><td>10/12</td><td>Venom component allergen</td><td>Lopez-Falcon et al., (2014) [96]</td></tr><tr><td></td><td></td><td></td><td>Immunity; energy metabolism; lipid transport; development</td><td>Geiser et al., (2012) [97]; Rodriguez Garcia et al., (2021) [98]</td></tr><tr><td></td><td></td><td></td><td></td><td>Peiren et al., (2006) [99]</td></tr><tr><td></td><td></td><td></td><td></td><td>Leipart et al., (2022a,b) [100]</td></tr></tbody></table>

### 3. Discussion

Provided with 3D tube-shaped Parafilm™ dummies filled with nutritional liquids, V. de-structer females are able to feed under laboratory conditions. Specifically, in flat Petri dishes (Ø 5 cm), the parasite can detect odorless dummies and pierce the Parafilm™ membrane to ingest a liquid diet, whether synthetic or natural. This finding agrees with our previous work showing that under artificial conditions and deprived of olfactive cues, the parasite's detection of a food source relies on shape-related information [21]. Using Parafilm™ dummies filled with the appropriate feeding solution, we were thus able to maintain mites for weeks. Although rearing methods have already been described [18–20,22], the capacity to feed a high number of mites systematically and easily under laboratory conditions is new and promising. Furthermore, the experimental design we developed in this work allowed us to test a variety of diets either completely synthetic or of natural origin.
Larval, pupal, or adult honey bee hemolymph has long been described as the only component of the *V. destructor* diet in natural conditions, until a recent study showed that adult fat body could also be targeted by the mite [22]. In our artificial conditions, the consumption of hemolymph collected from larvae, pupae, or adults was associated with a high survival rate after seven days of rearing. Pure undiluted heated hemolymph is thus sufficient to maintain mites for a week or more. After the first week, the mites survived longer when fed on larval rather than pupal and adult hemolymph. The composition of insect hemolymph varies between developmental stages [101–103]. Under our experimental conditions, larval hemolymph could represent an energetically richer food source more suitable to maintain *V. destructor*, whereas pupal and adult hemolymph could be less nutritional. In a natural environment, this lack of nutrients in the adult hemolymph could be compensated for by the consumption of fat body, as described by Ramsey et al. [22]. Since the hemolymph we collected from adult heads is unlikely to contain fat body cells [104], the nutrient deficiency could lead to premature death of *V. destructor* in our experiment. However, our findings do not completely match the work of Ramsey et al. [22], as mites fed on untreated adult hemolymph through a flat Parafilm™ membrane did not survive for more than two days in their study. This divergence could be due to our experimental designs, which differ both in the shape of the Parafilm™ membrane used and in the treatment applied to hemolymph. Indeed, in our study, adult hemolymph was heated for seven minutes at 65 °C to limit early melanization or contamination of the food, which could lead to survival differences in fed mites. This hypothesis, however, seems unlikely, as a survival rate of 80% [CI95: 59.3–93.2] was later observed on day 7 when using untreated adult hemolymph (Table S2). On the other hand, compared to flat Parafilm™ membranes, the 3D rounded shape of the feeding membrane used in our study had an arrestant effect on the mite [21]. This might facilitate the finding of food sources or even stimulate feeding behavior, which could explain the differences with previous studies [22,105].

In any case, the honey bee larva was the most suitable stage for sampling hemolymph as a food source for *V. destructor*. Indeed, in addition to the easy collection of large volumes with only a few individuals, parasites fed on larval hemolymph survived for more than 2 weeks in our study. This longevity is similar to what has been observed when the mite is kept on a natural host under laboratory conditions [106]. This means that female mites obtain the nutrients they need from larval hemolymph, either directly from cell-free hemolymph or from crude hemolymph. We showed that filtered unheated larval hemolymph (i.e., plasma) was as effective as heated hemolymph to feed parasites for weeks. As a result, free-floating hemocytes in plasma appear unnecessary for *V. destructor* nutrition, and plasma should therefore contain the vast majority of nutrients required for mite survival. Studying honey bee hemolymph plasma molecular composition is therefore a priority for understanding the physiological and nutritional needs of *V. destructor*.

Honey bee hemolymph is known to contain a variety of carbohydrates, lipids, free amino acids, and proteins [102,107,108]. Simple carbohydrates such as glucose, fructose, and trehalose are among the most abundant [29,109–111]. A mere fructose and glucose solution can thus simulate an extremely simplified version of bee cell-free hemolymph. In our bioassays, the use of such solutions led to a non-negligible survival gain compared with our control PBS solution, reaching levels similar to or higher than those observed with standard complex synthetic solutions [18–20]. Glucose and fructose thus represent two of the most important nutrients sought by the parasite. In addition to carbohydrates, other metabolites, such as fatty acids or carboxylic acids, were detected in adult bee hemolymph [107] along with several abundant immune or lipid transport proteins and a variety of free amino acids [101,107,112,113]. In honey bee juvenile stages, studies are scarce and show that the protein content can vary between mobile feeding larvae and immobile fasting and metamorphosing pupae [101,102,114]. In our study, protein supplementation was tested through the addition of yeast extracts, which mostly consist of free amino acids and polypeptides [30]. The addition of yeast extracts to the sugar solution significantly improved the survival of artificially fed mites. Nonetheless, survival rates were far from
the results obtained with filtered hemolymph. This highlights that some crucial nutrients remain unknown and hidden in the complexity of the hemolymph, perhaps in the form of whole bee proteins. In *V. destructor*, some proteins obtained from the parasitized bee can indeed act as essential nutrients on their own and not through digestion. Honey bee Vitellogenin has long been known to be found in the parasite's reproductive organs and hemocoel. More recently, a dozen additional bee proteins were found in the mite after it fed on its host [26]. The amino acid and polypeptides mimicked by the addition of yeast extract thus could not be sufficient for the mite survival, and specific honey bee proteins absent from these artificial feeding solutions could be required by the parasite.

The exact protein content of bee larval cell-free hemolymph was studied in our experiment, and a list of 54 omnipresent proteins was highlighted in each of our samples. Although additional interesting candidate molecules may be missing from this list, these commonly found proteins likely include many of the most frequent and abundant nutrients in bee larval plasma. Unsurprisingly, as hemolymph came from the last larval instar before metamorphosis, the proteins involved in insect development were highly represented [115]. Immune factors were also abundant in our samples, which is understandable, as these proteins are secreted directly into the hemolymph, mostly by hemocytes and fat body tissue [116]. The infectious background of larvae and physical wounds inflicted on the bees during hemolymph collection could also have triggered an immune/inflammatory reaction [114,117]. Nutritionally, the proteins involved in energy metabolism or lipid transport and storage could be particularly interesting. Such proteins might also be present as digested peptides in insect hemolymph, which cannot be discriminated against by our proteomic analyses. In any case, *Varroa destructor* females start their energy-consuming reproductive phase as soon as they enter a fifth instar larva cell, and the first 12 h are crucial for the initiation of reproduction in terms of olfactive signals [118]. The food ingested during the first hours after infestation could also be essential, as it represents the initial energy input used during reproduction. In this case, Vitellogenin, Hexamerin, Apolipophorins, or Transferrins seem to be relevant candidates for the physiological requirements of mites. Recently, Ramsey et al. [26] showed that these honey bee proteins can be directly hijacked by the parasite and accumulate in the females' body and eggs. The effect of such proteins in artificial feeding media remains to be elucidated. In addition, proteins are not the only nutrients targeted by *V. destructor*, and future metabolomics analyses should shed light on lipids, carbohydrates, or vitamins that may be of particular importance to *V. destructor*.

In conclusion, this study brings us closer to understanding mites' physiological requirements. From a methodological point of view, the protocol presented here could provide a sustainable procedure for maintaining mites under laboratory conditions, even in the absence of bees. So far, the mean longevity of the mites has reached 3 weeks when fed larval hemolymph, with a maximum of 56 days observed twice. The results obtained with filtered hemolymph also highlight that it is theoretically possible to create a synthetic solution that enables the maintenance of mites for days. The characterization of the nutrients required in such a solution, although fastidious, would be of prime interest for *Varroa* research. Compared with existing artificial feeding protocols for *V. destructor*, which rely on complex synthetic solutions made up of a mixture of cell culture media [18–21], we described here a simplified sugar solution that achieves similar survival rates, and shed light on several candidate proteins that need to be tested under artificial conditions. This new feeding protocol will be a helpful tool for research about many diverse aspects of *V. destructor* biology and could give rise to new questions about the nutrition of other hematophagous parasites.

It will also allow the development of sustainable rearing techniques without any living hosts, opening a path to new studies under entirely controlled conditions. This could lead to new opportunities in exploring biological control methods against this parasite and reducing its threat to honey bees.
# 4. Materials and Methods

## 4.1. Biological Material

### 4.1.1. Mites and Bees

Our studies were conducted according to the European ethics laws for scientific research currently in force (Directive 2010/63/EU of the European Parliament and the Council of 22 September 2010 on the protection of animals used for scientific purposes). Eight infested honey bee colonies were maintained on the university campus (Albi, France). *Varroa destructor* adults were collected during spring, summer, and autumn 2022 from sealed brood frames 8 to 12 days postcapping. Only mature females were used as males and juvenile stages are unable to pierce the bee cuticle and depend on their mother to feed. Young and not fully melanized adult females were also excluded.

### 4.1.2. Hemolymph Collection

Regardless of the honey bee stage, hemolymph collection was carried out on ice to maintain low temperature and avoid rapid melanization and proteolysis. Larval hemolymph was sampled by puncturing larvae with sharp tweezers or entomological pins. A drop of hemolymph was collected using a micropipette and the volume was transferred into a 1.5 mL microcentrifuge tube. About 30 larvae are needed to collect around 1 mL of larval hemolymph. Pupal hemolymph was sampled following a similar process. The abdomen was punctured and the thorax was pressed slightly to collect a droplet of hemolymph with a micropipette. At this stage, the hemolymph is full of suspended fat body cells [26,119], so the collected samples were not as clear and pure as larvae or adults' hemolymph. Around 30 pupae are necessary to gather 400 µL of pupal hemolymph.

To sample adult hemolymph, freshly freeze-killed worker bees from brood frames were used, and the method described by Borsuk et al., 2017 [104] was followed. Around 160 bees were kept on ice, and their antennae were removed. The bee thorax and abdomen were slightly pressed, and a droplet of hemolymph was collected from the site of the removed antenna using a micropipette. The droplets were transferred into a 1.5 mL microcentrifuge tube kept on ice until a sufficient volume was acquired. Around 160 bees allow the collection of 400 µL of hemolymph.

## 4.2. Artificial Rearing

### 4.2.1. Hemolymph as a Feeding Solution

During preliminary testing, pure untreated larval hemolymph was fed to the mites. However, due to risks of clotting and contamination, untreated hemolymph can quickly become unsuitable for *V. destructor* feeding after a few days at 34 °C (Figure S1). To avoid contamination and limit clotting or melanization, freshly collected larval, pupal, or adult hemolymph was heated at 65 °C for 7 min (as in [20]) before being included into Parafilm™ dummies. In our first bioassay, we explored the impact that the bee developmental stage used as a source of hemolymph could have on *V. destructor* survival (Table 1).

In a second bioassay, different treatments were applied to larval hemolymph. More specifically, filtration at 0.2 µm and heating at 65 °C for 7 min were conducted in order to obtain cell-free hemolymph with limited risks of melanization and contamination [20]. Heating also released the cell contents into the medium, as eukaryotic cells exposed to high temperature (above 55 °C) suffer instant necrosis after membrane rupture [27,28]. Solutions were prepared as described in Table 1 and stored at -20 °C. We also ensured that cold storage of our nutritive solutions was harmless by comparing the outcome of artificial feeding using freshly collected or frozen hemolymph (Figure S2).

### 4.2.2. Synthetic Diet

As alternatives to hemolymph, synthetic solutions were fed to mites in a third bioassay (Table 1), and PBS was used as control. In the first feeding condition, a mix of fructose (50 mg/mL) and glucose (50 mg/mL) was added to PBS. These carbohydrates are nutrients naturally present within the hemolymph in the form of mono- or disaccharides. More
precisely, in addition to trehalose (30 mg/mL), both glucose and fructose were detected at a concentration of approximately 5 mg/mL in adult hemolymph [29]. The concentrations used here, higher than those observed in natural adult bee hemolymph, were chosen based on preliminary results (Figure S3). In our second feeding condition, yeast extract was added to the carbohydrate solution to simulate proteins and vitamins uptakes [30]. Again, concentrations of yeast extracts (40 mg/mL) were chosen to be included within a suitable range determined during the preliminary experiments (Figure S3). In these two conditions, as well as in the PBS control, feeding solutions were colored using 0.5% artificial Blue-FCF food dye (Vahiné, France) to assess the feeding success of the mites. In all cases, the synthetic feeding solutions were filtered at 0.2 µm to limit the risk of microbial contamination.

### 4.2.3. Artificial Feeding Chambers and Bioassays

Rearing chambers consisted of Petri dishes (Ø 5 cm) covered with stretched Parafilm™ to limit the effect of static electricity observed on plastic surfaces [120]. The chambers were sterilized under UV light for 20 min. Under a fume hood, two tube-shaped dummies were fashioned from 2 × 2 cm Parafilm™ pieces fully stretched to a thickness of around 16 µm according to Posada-Florez et al. [19]. First, 90 microliters of feeding medium were deposited at the center of the stretched Parafilm™. The Parafilm™ was then carefully folded to create a tube-shaped dummy (Figure 1). In each chamber, 10 to 12 mites from the same colony can be followed at once. Under all conditions, the feeding solutions were renewed twice a week. The survival of mites was recorded daily for every condition through simple observation under a stereomicroscope.

In the first and second bioassays, parasite survival was investigated over 56 and 21 days, respectively. In the case of synthetic diet (Bioassay 3), the survival was also recorded daily, and the feeding success of mites was assessed after 24 h using the dye included in the medium. The blue dye can be seen throughout the parasite's cuticle, although dissection can be performed to confirm it [121] (Figure 4). An additional sample of 30 mites was also fed on colored larval hemolymph to allow comparison with synthetic diets.

## 4.3. *Proteomic Analyses*

### 4.3.1. Sample Processing and Bottom-up Proteomics

The filtered bee larval hemolymph samples (100 µL) were treated with an equal volume of pure hexafluoroisopropanol surfactant (HFIP, Sigma Aldrich, Saint-Quentin-Fallavier, France) and incubated at 4 °C for 4 h while shaking. The HFIP was totally evaporated by centrifugation under vacuum (Labconco, Kansas City, MO, USA), and the hemolymph proteins were reduced, alkylated, and digested as previously described [122,123]. The samples were then stored at -20 °C and dried out prior to proteomics analyses. Lyophilized samples were resuspended in formic acid 0.1% in water. Protein concentration was determined using Pierce 660-nm Protein Assay Reagent (Thermofisher, Rockford, IL, USA). Samples were analyzed on a nanoUPLC system (nanoAcquity, Waters, Milford, MA, USA) coupled to a quadrupole-Orbitrap hybrid mass spectrometer (Q-Exactive plus, Thermo Scientific, San Jose, CA, USA). The system was fully controlled using XCalibur software (v3.0.63; Thermo Fisher Scientific, Waltham, MA, USA). The UPLC system was equipped with a Symmetry C₁₈ precolumn (300 Å, 20 × 0.18 mm, 5 µm particle size, Waters, Milford, MA, USA) and an ACQUITY UPLC® BEH C₁₈ separation column (130 Å, 75 µm × 200 mm, 1.7 µm particle size, Waters). The solvent system consisted of formic acid 0.1% in water (solvent A) and 0.1% formic acid in acetonitrile (solvent B). Peptides (500 ng for each sample) were first trapped for 3 min at 5 µL/min with 99% A and 1% B. Elution was carried out at 60 °C at a flow rate of 350 nL/min, using a linear gradient from 2 to 35% B in 126 min, then 70% B in 4 min and equilibration was carried out for 14 min at 2% B. To minimize carry-over, a column wash (50% acetonitrile for 20 min) was performed between each sample in addition to a solvent blank.
The Q-Exactive Plus was operated in positive ion mode with the source temperature set to 250 °C and spray voltage to 1.8 kV. Full-scan MS spectra (300–1800 *m/z*) were acquired at a resolution of 140,000 at *m/z* 200, with a maximum injection time of 50 ms and an AGC target value of 3 × 10⁶ charges. The lock-mass option was enabled (polysiloxane, 445.12002 *m/z*). Up to 10 most intense peptides (at least doubly charged) per full scan were isolated using a 2 *m/z* window and they were fragmented using higher energy collisional dissociation (normalized collision energy of 27 eV, and dynamic exclusion of already fragmented precursors set to 60 s). MS/MS spectra (200–2000 *m/z*) were acquired at a resolution of 17,500 at *m/z* 200, with a maximum injection time of 100 ms and an automatic gain control (AGC) target value of 1 × 10⁵. The peptide match selection option was turned on. Peak intensities and retention times of reference peptides were monitored in daily.

### 4.3.2. Proteomics Data Processing

Proteome Discoverer 2.5 (Thermo Fisher Scientific, Waltham, MA, USA) was used for bottom-up proteomics sequencing and to quantify the proteins based on the intensities of the precursor ions of unique and razor peptides. A processing workflow and a consensus workflow were used to characterize and validate/quantify the proteins in each analyzed sample, respectively. A Fasta file protein database containing entries from NCBI gathered 1,532,988 sequences of Hymenoptera, 20,214 protein sequences from *Aethina tumida*, 14,703 sequences from *Tropilaelaps* species, 59,574 sequences from *Varroa* species, 7578 entries from bee viruses, 32,221 entries from *Crithidia* and *Lotmaria* species, and 115,707 entries from mites. The trypsin enzyme, a precursor mass tolerance of 20 ppm, a fragment mass tolerance of 0.5 Da, and the chemical modifications (Oxidation/+15.995 Da as dynamic on methionine and tryptophan residues and pyridylethyl/+105.058 Da as static on cysteine) were selected for subsequent data processing. The datasets were aligned to extract for each individual file, the LC-MS/MS mapping features. The settings were as follows: A maximum retention time (RT) shift of 10 min., a mass tolerance of 10 ppm, and coarse parameter tuning were set. Protein annotations were successfully validated when the q-value scores were calculated below the false discovery rate (0.05).

### 4.4. Statistical Analyses

Results were analyzed using R.4.0.4 [124] and graphs were generated using the ggplot2 package [125]. When mortality on day 7 and 14 was compared, binary data (survival 0 or 1) were analyzed using generalized linear models (GLM) with binomial distribution. The feeding success after 24 h was analyzed following the same process.

Survival curves were fitted to the data to follow mortality over time. Data were analyzed using the survfit package in R and log rank tests. If significance was reached, pairwise comparisons with Benjamini-Hochberg corrections were further computed.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241512443/s1.

**Author Contributions:** Conceptualization, V.P., C.V., S.B. and A.V.; methodology, V.P., C.V., S.B., J.-M.S., F.B., M.B., K.A. and P.B.; formal analysis, V.P., K.A. and M.B.; investigation, V.P.; resources, V.P., C.V., S.B., J.-M.S., F.B., M.B., K.A., P.B. and A.V.; data curation, K.A., M.B. and P.B.; writing—original draft preparation, V.P.; writing—review and editing, V.P., C.V., S.B., J.-M.S., F.B., M.B., K.A. and P.B.; visualization, V.P.; supervision, A.V. and P.B.; funding acquisition, A.V. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Region Occitanie (ADEME grant No. 2082C0061). We thank the Association Plateforme BioPark d'Archamps (PBA, France) for supporting part of this research through its research and development program.

**Institutional Review Board Statement:** No ethical approval required.

**Informed Consent Statement:** Not applicable.
**Data Availability Statement:** All data are available upon request. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org, accessed on 10 July 2023) via the PRIDE [126] partner repository with the dataset identifier PXD043487.

**Acknowledgments:** We are thankful to ADA Occitanie for the biological material.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Chalkowski, K.; Lepczyk, C.A.; Zohdy, S. *Parasite Ecology of Invasive Species: Conceptual Framework and New Hypotheses*. *Trends Parasitol.* **2018**, *34*, 655–663. [CrossRef]

2. Chapman, N.C.; Colin, T.; Cook, J.; Da Silva, C.R.B.; Gloag, R.; Hogendoorn, K.; Howard, S.R.; Remnant, E.J.; Roberts, J.M.K.; Tierney, S.M.; et al. The final frontier: Ecological and evolutionary dynamics of a global parasite invasion. *Biol. Lett.* **2023**, *5*, 20220589. [CrossRef]

3. Insolia, L.; Molinari, R.; Rogers, S.R.; Williams, G.R.; Chiaromonte, F.; Calovi, M. Honey bee colony loss linked to parasites, pesticides and extreme weather across the United States. *Sci. Rep.* **2022**, *12*, 20787. [CrossRef]

4. Grupe II, A.C.; Quandt, C.A. A growing pandemic: A review of *Nosema* parasites in globally distributed domesticated and native bees. *PLoS Pathog.* **2020**, *16*, e1008580. [CrossRef]

5. Marcogliese, D.J.; Pietrock, M. Combined effects of parasites and contaminants on animal health: Parasites do matter. *Trends Parasitol.* **2011**, *27*, 123–130. [CrossRef]

6. Giuffre, C.; Lubkin, S.R.; Tarpy, D.R. Does viral load alter behavior of the bee parasite *Varroa destructor*? *PLoS ONE* **2019**, *14*, e0217975. [CrossRef]

7. Thaduri, S.; Stephan, J.G.; de Miranda, J.R.; Locke, B. Disentangling host-parasite-pathogen interactions in a varroa-resistant honeybee population reveals virus tolerance as an independent, naturally adapted survival mechanism. *Sci. Rep.* **2019**, *9*, 6221. [CrossRef]

8. Mazier, D.; Beaudoin, R.L.; Mellouk, S.; Druilhe, P.; Texier, B.; Trosper, J.; Miltgen, F.; Landau, I.; Paul, C.; Brandicourt, O.; et al. Complete Development of Hepatic Stages of *Plasmodium falciparum* In Vitro. *Science* **1985**, *227*, 440–442. [CrossRef]

9. Fivelman, Q.L.; McRobert, L.; Sharp, S.; Taylor, C.J.; Saeed, M.; Swales, C.A.; Sutherland, C.J.; Baker, D.A. Improved synchronous production of Plasmodium falciparum gametocytes in vitro. *Mol. Biochem. Parasitol.* **2007**, *154*, 119–123. [CrossRef]

10. Waladde, S.M.; Young, A.S.; Morzaria, S.P. Artificial feeding of ixodid ticks. *Parasitol. Today* **1996**, *12*, 272–278. [CrossRef]

11. Montes, C.; Cuadrillero, C.; Vilella, D. Maintenance of a laboratory colony of *Cimex lectularius* (Hemiptera: Cimicidae) using an artificial feeding technique. *J. Med. Entomol.* **2002**, *39*, 675–679. [CrossRef]

12. Bonnet, S.; Liu, X. Laboratory artificial infection of hard ticks: A tool for the analysis of tick-borne pathogen transmission. *Acarologia* **2012**, *52*, 453–464. [CrossRef]

13. Bonnet, S.; Jouglin, M.; Malandrin, L.; Becker, C.; Agoulon, A.; L'Hostis, M.; Chauvin, A. Transstadial and transovarial persistence of *Babesia divergens* DNA in *Ixodes ricinus* ticks fed on infected blood in a new skin-feeding technique. *Parasitology* **2007**, *134*, 197–207. [CrossRef]

14. Kröber, T.; Guerin, P.M. An in vitro feeding assay to test acaricides for control of hard ticks. *Pest Manag. Sci.* **2008**, *63*, 17–22. [CrossRef]

15. Amiri, E.; Waiker, P.; Rueppell, O.; Manda, P. Using manual and computer-based text-mining to uncover research trends for *Apis mellifera*. *Vet. Sci.* **2020**, *7*, 61. [CrossRef]

16. Traynor, K.S.; Mondet, F.; de Miranda, J.R.; Techer, M.; Kowallik, V.; Oddie, M.A.Y.; Chantawannakul, P.; McAfee, A. *Varroa destructor*: A complex parasite, crippling honey bees worldwide. *Trends Parasitol.* **2020**, *36*, 592–606. [CrossRef]

17. Vilarem, C.; Piou, V.; Vogelweith, F.; Vétillard, A. *Varroa destructor* from the laboratory to the field: Control, biocontrol and IPM perspectives—A review. *Insects* **2021**, *12*, 800. [CrossRef]

18. Bruce, W.A.; Chiesa, F.; Marchetti, S.; Griffiths, D.A. Laboratory feeding of *Varroa jacobsoni* Oudemans on natural and artificial diets (Acari: Varroidae). *Apidologie* **1988**, *19*, 209–218. [CrossRef]

19. Posada-Florez, F.; Ryabov, E.V.; Heerman, M.C.; Chen, Y.; Evans, J.D.; Sonenshine, D.E.; Cook, S.C. *Varroa destructor* mites vector and transmit pathogenic honey bee viruses acquired from an artificial diet. *PLoS ONE* **2020**, *15*, e0242688. [CrossRef]

20. Tabart, J.; Colin, M.E.; Carayon, J.L.; Tene, N.; Payre, B.; Vetillard, A. Artificial feeding of *Varroa destructor* through a chitosan membrane: A tool for studying the host-microparasite relationship. *Exp. Appl. Acarol.* **2013**, *61*, 107–118. [CrossRef]

21. Piou, V.; Vilarem, C.; Blanchard, S.; Armengaud, C.; Heeb, P.; Vétillard, A. A foraging enigma: How does *Varroa destructor* find a food source in artificial conditions? *Parasite*, 2023; under review.

22. Ramsey, S.D.; Ochoa, R.; Bauchan, G.; Gulbronson, C.; Mowery, J.D.; Cohen, A.; Lim, D.; Joklik, J.; Cicero, J.M.; Ellis, J.D.; et al. *Varroa destructor* feeds primarily on honey bee fat body tissue and not hemolymph. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 1792–1801. [CrossRef]

23. Annoscia, D.; Brown, S.P.; Di Prisco, G.; De Paoli, E.; Del Fabbro, S.; Frizzera, D.; Zanni, V.; Galbraith, D.A.; Caprio, E.; Grozinger, C.M.; et al. Haemolymph removal by *Varroa mite* destabilizes the dynamical interaction between immune effectors and virus in bees, as predicted by Volterra's model. *Proc. R. Soc. B Biol. Sci.* **2019**, *268*, 20190331. [CrossRef]
24. Tewarson, N.C.; Engels, W.; Tewarson, N.C. Undigested Uptake of Non-Host Proteins by *Varroa jacobsoni*. *J. Apic. Res.* **1982**, *21*, 222–225. [CrossRef]

25. Weinberg, K.P.; Madel, G. The Influence of the Mite *Varroa jacobsoni* Oud. on the Protein Concentration and the Haemolymph Volume of the Brood of Worker Bees and Drones of the Honey Bee *Apis mellifera* L. *Apidologie* **1985**, *16*, 421–436. [CrossRef]

26. Ramsey, S.D.; Cook, S.C.; Gulbronson, C.; Vanengelsdorp, D.; Evans, J.; Posada, F.; Sonenshine, D. Kleptocytosis: A Novel Parasitic Strategy for Accelerated Reproduction via Host Protein Stealing in *Varroa destructor*. *bioRxiv*, 2022; preprint.

27. Janko, C.; Munoz, L.; Chaurio, R.; Maueröder, C.; Berens, C.; Lauber, K.; Herrmann, M. Navigation to the graveyard-induction of various pathways of necrosis and their classification by flow cytometry. *Methods Mol. Biol.* **2013**, *1004*, 3–15. [CrossRef]

28. Yan, G.; Elbadawi, M.; Efferth, T. Multiple cell death modalities and their key features (Review). *World Acad. Sci. J.* **2020**, *2*, 39–48. [CrossRef]

29. Blatt, J.; Roces, F. Haemolymph sugar levels in foraging honeybees (*Apis mellifera carnia*): Dependence on metabolic rate and in vivo measurement of maximal rates of trehalose synthesis. *J. Exp. Biol.* **2001**, *204*, 2709–2716. [CrossRef]

30. Tomé, D. Yeast Extracts: Nutritional and Flavoring Food Ingredients. *ACS Food Sci. Technol.* **2021**, *1*, 487–494. [CrossRef]

31. Rowley, A.F.; Ratcliffe, N.A. Insect erythrocyte agglutinins. In Vitro opsonization experiments with Clitumnus extradentatus and Periplaneta americana haemocytes. *Immunology* **1980**, *40*, 483–492.

32. Casteels, P.; Ampe, C.; Riviere, L.; Van Damme, J.; Elicone, C.; Fleming, M.; Jacobs, F.; Tempst, P. Isolation and characterization of abaecin, a major antibacterial response peptide in the honeybee (*Apis mellifera*). *Eur. J. Biochem.* **1990**, *187*, 381–386. [CrossRef]

33. Pluta, P.; Sokól, R. Changes in the expression of antimicrobial peptide genes in honey bees (*Apis mellifera*) under the influence of various pathogens. *Ann. Parasitol.* **2020**, *66*, 457–465.

34. Wittwer, F.; van der Straten, A.; Keleman, K.; Dickson, B.J.; Hafen, E. Lilliputian: An AF4/FMR2-related protein that controls cell identity and cell growth. *Development* **2001**, *128*, 791–800. [CrossRef]

35. Cabbri, R.; Ferlizza, E.; Nanetti, A.; Monari, E.; Andreani, G.; Galuppi, R.; Isani, G. Biomarkers of nutritional status in honeybee haemolymph: Effects of different biotechnical approaches for *Varroa destructor* treatment and wintering phase. *Apidologie* **2018**, *49*, 606–618. [CrossRef]

36. Kim, B.Y.; Jin, B.R. Apolipophorin III from honeybees (*Apis cerana*) exhibits antibacterial activity. *Comp. Biochem. Physiol. Part B Biochem. Mol. Biol.* **2015**, *182*, 6–13. [CrossRef]

37. Waladde, S.M.; Rice, M.J. The Sensory Basis of Tick Feeding Behaviour. In *Physiology of Ticks*; Obenchain, F.D., Galun, R., Eds.; Pergamon Press Ltd.: Amsterdam, The Netherlands, 1982; pp. 71–118.

38. Peng, Y.S. Activity of β-Galactosidase in the Midgut of the Honeybee. *J. Apic. Res.* **1980**, *19*, 105–111. [CrossRef]

39. Ricigliano, V.A.; Fitz, W.; Copeland, D.C.; Mott, B.M.; Maes, P.; Floyd, A.S.; Dockstader, A.; Anderson, K.E. The impact of pollen consumption on honey bee (*Apis mellifera*) digestive physiology and carbohydrate metabolism. *Arch. Insect Biochem. Physiol.* **2017**, *96*, e21406. [CrossRef]

40. Tremblay, G.B.; Sohi, S.S.; Retnakaran, A.; MacKenzie, R.E. NAD-dependent methylenetetrahydrofolate dehydrogenase-methenyltetrahydrofolate cyclohydrolase is targeted to the cytoplasm in insect cell lines. *FEBS Lett.* **1995**, *368*, 177–182. [CrossRef]

41. Yu, S.; Jang, Y.; Paik, D.; Lee, E.; Park, J.J. Nmdmc overexpression extends Drosophila lifespan and reduces levels of mitochondrial reactive oxygen species. *Biochem. Biophys. Res. Commun.* **2015**, *465*, 845–850. [CrossRef]

42. Upadhyay, A.; Moss-Taylor, L.; Kim, M.J.; Ghosh, A.C.; O'Connor, M.B. TGF-β Family Signaling in Drosophila. *Cold Spring Harb. Perspect. Biol.* **2017**, *9*, a022152. [CrossRef]

43. Li, J.; Wu, J.; Rundassa, D.B.; Song, F.; Zheng, A.; Fang, Y. Differential protein expression in honeybee (*Apis mellifera* L.) larvae: Underlying caste differentiation. *PLoS ONE* **2010**, *5*, e13455. [CrossRef]

44. Kim, B.Y.; Lee, K.S.; Wan, H.; Zou, F.M.; Choi, Y.S.; Yoon, H.J.; Kwon, H.W.; Je, Y.H.; Jin, B.R. Anti-elastolytic activity of a honeybee (*Apis cerana*) chymotrypsin inhibitor. *Biochem. Biophys. Res. Commun.* **2013**, *430*, 144–149. [CrossRef]

45. Corral-Rodríguez, M.Á.; Macedo-Ribeiro, S.; Barbosa Pereira, P.J.; Fuentes-Prior, P. Tick-derived Kunitz-type inhibitors as antihemostatic factors. *Insect Biochem. Mol. Biol.* **2009**, *39*, 579–595. [CrossRef]

46. Rhoads, M.L.; Fetterer, R.H.; Hill, D.E.; Urban, J.F. *Trichuris suis*: A secretory chymotrypsin/elastase inhibitor with potential as an immunomodulator. *Exp. Parasitol.* **2000**, *95*, 36–44. [CrossRef]

47. Ledent, V.; Vervoort, M. The basic helix-loop-helix protein family: Comparative genomics and phylogenetic analysis. *Genome Res.* **2001**, *11*, 754–770. [CrossRef]

48. Wan, P.J.; Yuan, S.Y.; Wang, W.X.; Chen, X.; Lai, F.X.; Fu, Q. A genome-wide identification and analysis of the basic helix-loop-helix transcription factors in brown planthopper, *Nilaparvata lugens*. *Genes* **2016**, *7*, 100. [CrossRef]

49. Wang, Y.; Chen, K.; Yao, Q.; Wang, W.; Zhu, Z. The basic helix-loop-helix transcription factor family in the honey bee, *Apis mellifera*. *J. Insect Sci.* **2008**, *8*, 40. [CrossRef]

50. Pastor-Pareja, J.C.; Xu, T. Shaping Cells and Organs in *Drosophila* by Opposing Roles of Fat Body-Secreted Collagen IV and Perlecan. *Dev. Cell* **2011**, *21*, 245–256. [CrossRef]

51. Sutherland, T.D.; Peng, Y.Y.; Trueman, H.E.; Weisman, S.; Okada, S.; Walker, A.A.; Sriskantha, A.; White, J.F.; Huson, M.G.; Werkmeister, J.A.; et al. A new class of animal collagen masquerading as an insect silk. *Sci. Rep.* **2013**, *3*, 2–7. [CrossRef]

52. Zhang, F.; Zhao, Y.; Chao, Y.; Muir, K.; Han, Z. Cubilin and amnionless mediate protein reabsorption in *Drosophila* nephrocytes. *J. Am. Soc. Nephrol.* **2013**, *24*, 209–216. [CrossRef]
53. Hopkins, T.L.; John Krchma, L.; Ahmad, S.A.; Kramer, K.J. Pupal cuticle proteins of *Manduca sexta*: Characterization and profiles during sclerotization. *Insect Biochem. Mol. Biol.* **2000**, *30*, 19–27. [CrossRef]

54. Soares, M.P.M.; Elias-Neto, M.; Simões, Z.L.P.; Bitondi, M.M.G. A cuticle protein gene in the honeybee: Expression during development and in relation to the ecdysteroid titer. *Insect Biochem. Mol. Biol.* **2007**, *37*, 1272–1282. [CrossRef]

55. Micas, A.F.D.; Ferreira, G.A.; Laure, H.J.; Rosa, J.C.; Bitondi, M.M.G. Proteins of the integumentary system of the honeybee, *Apis mellifera*. *Arch. Insect Biochem. Physiol.* **2016**, *93*, 3–24. [CrossRef]

56. Caccia, S.; Grimaldi, A.; Casartelli, M.; Falabella, P.; de Eguileor, M.; Pennacchio, F.; Giordana, B. Functional analysis of a fatty acid binding protein produced by *Aphidius ervi* teratocytes. *J. Insect Physiol.* **2012**, *58*, 621–627. [CrossRef]

57. Chen, J.; Zou, X.; Zhu, W.; Duan, Y.; Merzendorfer, H.; Zhao, Z.; Yang, Q. Fatty acid binding protein is required for chitin biosynthesis in the wing of *Drosophila melanogaster*. *Insect Biochem. Mol. Biol.* **2022**, *149*, 103845. [CrossRef]

58. Cheng, L.; Jin, X.K.; Li, W.W.; Li, S.; Guo, X.N.; Wang, J.; Gong, Y.N.; He, L.; Wang, Q. Fatty Acid Binding Proteins FABP9 and FABP10 Participate in Antibacterial Responses in Chinese Mitten Crab, *Eriocheir sinensis*. *PLoS ONE* **2013**, *8*, e54053. [CrossRef]

59. Piha-Gossack, A.; Sossin, W.; Reinhardt, D.P. The evolution of extracellular fibrillins and their functional domains. *PLoS ONE* **2012**, *7*, e33560. [CrossRef]

60. Sutherland, T.D.; Campbell, P.M.; Weisman, S.; Trueman, H.E.; Sriskantha, A.; Wanjura, W.J.; Haritos, V.S. A highly divergent gene cluster in honey bees encodes a novel silk family. *Genome Res.* **2006**, *16*, 1414–1421. [CrossRef]

61. Rebers, J.E.; Riddiford, L.M. Structure and expression of a Manduca sexta larval cuticle gene homologous to Drosophila cuticle genes. *J. Mol. Biol.* **1988**, *203*, 411–423. [CrossRef]

62. Zhao, X.; Zhang, J.; Yang, J.; Niu, N.; Zhang, J.; Yang, Q. Mucin family genes are essential for the growth and development of the migratory locust, *Locusta migratoria*. *Insect Biochem. Mol. Biol.* **2020**, *123*, 103404. [CrossRef]

63. Cox-Foster, D.L.; Schonbaum, C.P.; Murtha, M.T.; Cavener, D.R. Developmental expression of the glucose dehydrogenase gene in *Drosophila melanogaster*. *Genetics* **1990**, *124*, 873–880. [CrossRef]

64. Cox-Foster, D.L.; Stehr, J.E. Induction and localization of FAD-glucose dehydrogenase (GLD) during encapsulation of abiotic implants in Manduca sexta larvae. *J. Insect Physiol.* **1994**, *40*, 235–249. [CrossRef]

65. Zhang, J.; Goyer, C.; Pelletier, Y. Environmental stresses induce the expression of putative glycine-rich insect cuticular protein genes in adult *Leptinotarsa decemlineata* (Say). *Insect Mol. Biol.* **2008**, *17*, 209–216. [CrossRef]

66. Zhong, Y.S.; Mita, K.; Shimada, T.; Kawasaki, H. Glycine-rich protein genes, which encode a major component of the cuticle, have different developmental profiles from other cuticle protein genes in *Bombyx mori*. *Insect Biochem. Mol. Biol.* **2006**, *36*, 99–110. [CrossRef]

67. Martins, J.R.; Nunes, F.M.F.; Cristino, A.S.; Simões, Z.L.P.; Bitondi, M.M.G. The four hexamerin genes in the honey bee: Structure, molecular evolution and function deduced from expression patterns in queens, workers and drones. *BMC Mol. Biol.* **2010**, *11*, 23. [CrossRef]

68. Contreras, E.G.; Glavic, Á.; Brand, A.H.; Sierralta, J.A. The serine protease homolog, scarface, is sensitive to nutrient availability and modulates the development of the *Drosophila* blood-brain barrier. *J. Neurosci.* **2021**, *41*, 6430–6448. [CrossRef]

69. Ryan, R.O.; Law, J.H.; Shipman, B.A.; Schmidt, J.O. Purification and Properties of a Very High Density Lipoprotein from the Hemolymph of the Honeybee *Apis mellifera*. *Biochemistry* **1987**, *26*, 1885–1889. [CrossRef]

70. Arockiaraj, J.; Easwvaran, S.; Vanaraja, P.; Singh, A.; Othman, R.Y.; Bhassu, S. First report on interferon related developmental regulator-1 from Macrobrachium rosenbergii: Bioinformatic analysis and gene expression. *Fish Shellfish Immunol.* **2012**, *32*, 929–933. [CrossRef]

71. Hoffmann, J.A.; Reichhart, J.M.; Hetru, C. Innate immunity in higher insects. *Curr. Opin. Immunol.* **1996**, *8*, 8–13. [CrossRef]

72. Stanifer, M.L.; Pervolaraki, K.; Boulant, S. Differential regulation of type I and type III interferon signaling. *Int. J. Mol. Sci.* **2019**, *20*, 1445. [CrossRef]

73. Al-Ghamdi, A.A.; Al-Ghamdi, M.S.; Ahmed, A.M.; Mohamed, A.S.A.; Shaker, G.H.; Ansari, M.J.; Dorrah, M.A.; Khan, K.A.; Ayaad, T.H. Immune investigation of the honeybee *Apis mellifera* jemenitica broods: A step toward production of a bee-derived antibiotic against the American foulbrood. *Saudi J. Biol. Sci.* **2021**, *28*, 1528–1538. [CrossRef]

74. Bezabih, G.; Cheng, H.; Han, B.; Feng, M.; Xue, Y.; Hu, H.; Li, J. Phosphoproteome Analysis Reveals Phosphorylation Underpin- nings in the Brains of Nurse and Forager Honeybees (*Apis mellifera*). *Sci. Rep.* **2017**, *7*, 1973. [CrossRef]

75. Petzold, A. Neurofilament phosphoforms: Surrogate markers for axonal injury, degeneration and loss. *J. Neurol. Sci.* **2005**, *233*, 183–198. [CrossRef]

76. Forêt, S.; Maleszka, R. Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (*Apis mellifera*). *Genome Res.* **2006**, *16*, 1404–1413. [CrossRef]

77. Bloch, G.; Cohen, M. The expression and phylogenetics of the Inhibitor Cysteine Knot peptide OCLP1 in the honey bee *Apis mellifera*. *J. Insect Physiol.* **2014**, *65*, 1–8. [CrossRef]

78. Wang, P.; Heitman, J. The cyclophilins. *Genome Biol.* **2005**, *6*, 226. [CrossRef]

79. Yoon, J.S.; Kim, S.S.; Ha, J.; Kang, I.; Choe, W. Cyclophilin B, a molecule chaperone, promotes adipogenesis in 3T3-L1 preadipocytes via AKT/mTOR pathway. *Int. J. Mol. Med.* **2023**, *51*, 6. [CrossRef]

80. Park, M.J.; Kim, B.Y.; Jin, B.R. Molecular characterization of a chitin-binding protein with the peritrophin-A domain from the Asiatic honeybee *Apis cerana*. *J. Asia Pac. Entomol.* **2016**, *19*, 963–968. [CrossRef]
81. Wilson-Rich, N.; Dres, S.T.; Starks, P.T. The ontogeny of immunity: Development of innate immune strength in the honey bee (*Apis mellifera*). *J. Insect Physiol.* **2008**, *54*, 1392–1399. [CrossRef]

82. Miller-Delaney, S.F.C.; Lieberam, I.; Murphy, P.; Mitchell, K.J. Plxdc2 Is a Mitogen for Neural Progenitors. *PLoS ONE* **2011**, *6*, e14565. [CrossRef]

83. Thibord, F.; Hardy, L.; Ibrahim-Kosta, M.; Saut, N.; Pulcrano-Nicolas, A.S.; Goumidi, L.; Civelek, M.; Eriksson, P.; Deleuze, J.F.; Le Goff, W.; et al. A Genome Wide Association Study on plasma FV levels identified PLXDC2 as a new modifier of the coagulation process. *J. Thromb. Haemost.* **2019**, *17*, 1808–1814. [CrossRef]

84. Liu, Y.; Beaurepaire, A.; Rogers, C.W.; Lopez, D.; Evans, J.D.; Straub, L.; Neumann, P.; Cook, S.C.; Huang, Q. Gene expression and functional analyses of odorant receptors in small hive beetles (*Aethina tumida*). *Int. J. Mol. Sci.* **2020**, *21*, 4582. [CrossRef]

85. Jung, D.; Lee, J.; Park, T.Y.; Yang, Y.J.; Cha, H.J. Diverse silk and silk-like proteins derived from terrestrial and marine organisms and their applications. *Acta Biomater.* **2021**, *136*, 56–71. [CrossRef]

86. Liu, N.; Li, T.; Wang, Y.; Liu, S. G-protein coupled receptors (Gpcrs) in insects—A potential target for new insecticide development. *Molecules* **2021**, *26*, 2993. [CrossRef]

87. Jørgensen, O.S. D2-protein and D3-protein as markers for synaptic turnover and concentration. *J. Neural Transm. Suppl.* **1983**, *18*, 245–255.

88. Jonusaite, S.; Beyenbach, K.W.; Meyer, H.; Paululat, A.; Izumi, Y.; Furuse, M.; Rodan, A.R. The septate junction protein Mesh is required for epithelial morphogenesis, ion transport, and paracellular permeability in the *Drosophila* Malpighian tubule. *Am. J. Physiol. Cell Physiol.* **2020**, *318*, C675–C694. [CrossRef]

89. Majerowicz, D.; Hannibal-Bach, H.K.; Castro, R.S.C.; Bozaquel-Morais, B.L.; Alves-Bezerra, M.; Grillo, L.A.M.; Masuda, C.A.; Færgeman, N.J.; Knudsen, J.; Gondim, K.C. The ACBP gene family in *Rhodnius prolixus*: Expression, characterization and function of RpACBP-1. *Insect Biochem. Mol. Biol.* **2016**, *72*, 41–52. [CrossRef]

90. Zhao, G.; Wang, C.; Wang, H.; Gao, L.; Liu, Z.; Xu, B.; Guo, X. Characterization of the CDK5 gene in *Apis cerana cerana* (AccCDK5) and a preliminary identification of its activator gene, AccCDK5r1. *Cell Stress Chaperones* **2018**, *23*, 13–28. [CrossRef]

91. Al-Naggar, Y.; Shafiey, H.; Paxton, R.J. Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (*Apis mellifera*). *Viruses* **2023**, *15*, 1284. [CrossRef]

92. Doublet, V.; Poeschl, Y.; Gogol-Döring, A.; Alaux, C.; Annoscia, D.; Aurori, C.; Barribeau, S.M.; Bedoya-Reina, O.C.; Brown, M.J.F.; Bull, J.C.; et al. Unity in defence: Honeybee workers exhibit conserved molecular responses to diverse pathogens. *BMC Genom.* **2017**, *18*, 207. [CrossRef]

93. Shakeel, M.; Xu, X.; De Mandal, S.; Jin, F. Role of serine protease inhibitors in insect-host-pathogen interactions. *Arch. Insect Biochem. Physiol.* **2019**, *102*, e21556. [CrossRef]

94. Zhang, F.X.; Shao, H.L.; Wang, J.X.; Zhao, X.F. B-Thymosin Is Upregulated By the Steroid Hormone 20-Hydroxyecdysone and Microorganisms. *Insect Mol. Biol.* **2011**, *20*, 519–527. [CrossRef]

95. Hooper, S.L.; Thuma, J.B. Invertebrate muscles: Muscle specific genes and proteins. *Physiol. Rev.* **2005**, *85*, 1001–1060. [CrossRef] [PubMed]

96. López-Falcón, B.; Meyer-Nava, S.; Hernández-Rodríguez, B.; Campos, A.; Montero, D.; Rudiño, E.; Vázquez, M.; Zurita, M.; Valadez-Graham, V. Characterization of the *Drosophila* group ortholog to the amino-terminus of the alpha-thalassemia and mental retardation X-linked (ATRX) vertebrate protein. *PLoS ONE* **2014**, *9*, e113182. [CrossRef] [PubMed]

97. Geiser, D.L.; Winzerling, J.J. Insect transferrins: Multifunctional proteins. *Biochim. Biophys. Acta Gen. Subj.* **2012**, *1820*, 437–451. [CrossRef] [PubMed]

98. Rodríguez-García, C.; Heerman, M.C.; Cook, S.C.; Evans, J.D.; DeGrandi-Hoffman, G.; Banmeke, O.; Zhang, Y.; Huang, S.; Hamilton, M.; Chen, Y.P. Transferrin-mediated iron sequestration suggests a novel therapeutic strategy for controlling *Nosema* disease in the honey bee, *Apis mellifera*. *PLoS Pathog.* **2021**, *17*, e1009270. [CrossRef] [PubMed]

99. Peiren, N.; de Graaf, D.C.; Brunain, M.; Bridts, C.H.; Ebo, D.G.; Stevens, W.J.; Jacobs, F.J. Molecular cloning and expression of icarapin, a novel IgE-binding bee venom protein. *FEBS Lett.* **2006**, *580*, 4895–4899. [CrossRef]

100. Leipart, V.; Ludvigsen, J.; Kent, M.; Sandve, S.; To, T.H.; Árnyasi, M.; Kreibich, C.D.; Dahle, B.; Amdam, G.V. Identification of 121 variants of honey bee Vitellogenin protein sequences with structural differences at functional sites. *Protein Sci.* **2022**, *31*, e4369. [CrossRef]

101. Erban, T.; Harant, K.; Kamler, M.; Markovic, M.; Titera, D. Detailed proteome mapping of newly emerged honeybee worker hemolymph and comparison with the red-eye pupal stage. *Apidologie* **2016**, *47*, 805–817. [CrossRef]

102. Woltedji, D.; Fang, Y.; Han, B.; Feng, M.; Li, R.; Lu, X.; Li, J. Proteome analysis of hemolymph changes during the larval to pupal development stages of honeybee workers (*Apis mellifera ligustica*). *J. Proteome Res.* **2013**, *12*, 5189–5198. [CrossRef]

103. Zhou, L.; Li, H.; Hao, F.; Li, N.; Liu, X.; Wang, G.; Wang, Y.; Tang, H. Developmental Changes for the Hemolymph Metabolome of Silkworm (*Bombyx mori* L.). *J. Proteome Res.* **2015**, *176*, 2331–2347. [CrossRef]

104. Borsuk, G.; Ptaszyńska, A.A.; Olszewski, K.; Domaciuk, M.; Krutmuang, P.; Paleolog, J. A new method for quick and easy hemolymph collection from apidae adults. *PLoS ONE* **2017**, *12*, e0170487. [CrossRef] [PubMed]

105. Jack, C.J.; Dai, P.L.; van Santen, E.; Ellis, J.D. Comparing four methods of rearing *Varroa destructor* in vitro. *Exp. Appl. Acarol.* **2020**, *80*, 463–476. [CrossRef] [PubMed]
106. Egekwu, N.I.; Posada, F.; Sonenshine, D.E.; Cook, S. Using an in vitro system for maintaining *Varroa destructor* mites on *Apis mellifera* pupae as hosts: Studies of mite longevity and feeding behavior. *Exp. Appl. Acarol.* **2018**, *74*, 301–315. [CrossRef]

107. Aliferis, K.A.; Copley, T.; Jabaji, S. Gas chromatography-mass spectrometry metabolite profiling of worker honey bee (*Apis mellifera* L.) hemolymph for the study of Nosema ceranae infection. *J. Insect Physiol.* **2012**, *58*, 1349–1359. [CrossRef]

108. Cremonez, T.M.; De Jong, D.; Bitondi, M.M.G. Quantification of Hemolymph Proteins as a Fast Method for Testing Protein Diets for Honey Bees (Hymenoptera: Apidae). *J. Econ. Entomol.* **1998**, *91*, 1284–1289. [CrossRef]

109. Cournoyer, A.; Plamondon, L.; Bau-Gaudreault, L.; Deschamps, A.; Dubreuil, P.; Benoit-Biancamano, M.O. Effects of *Varroa destructor* on Hemolymph Sugars and Secondary Infections in Honeybees (*Apis mellifera*). *Appl. Sci.* **2022**, *12*, 11630. [CrossRef]

110. Tsao, W.; Shuel, R.W. Studies in the mode of action of royal jelly in honeybee development IX. The carbohydrates and lipids in the haemolymph and the fat body of developing larvae. *Can. J. Zool.* **1973**, *51*, 1139–1148. [CrossRef]

111. Wyatt, G.R.; Kalf, G.F. The chemistry of insect hemolymph II. Trehalose and other carbohydrates. *J. Gen. Physiol.* **1957**, *40*, 833–847. [CrossRef]

112. Micheu, S.; Crailsheim, K.; Leonhard, B. Importance of proline and other amino acids during honeybee flight (*Apis mellifera carnica* Pollmann). *Amino Acids* **2000**, *18*, 157–175. [CrossRef]

113. Wang, D.I.; Moeller, F.E. Comparison of the free amino acid composition in the hemolymph of healthy and *Nosema*-infected female honey bees. *J. Invertebr. Pathol.* **1970**, *15*, 202–206. [CrossRef]

114. Randolt, K.; Gimple, O.; Geissendörfer, J.; Reinders, J.; Prusko, C.; Mueller, M.J.; Albert, S.; Tautz, J.; Beier, H. Immune-related proteins induced in the hemolymph after aseptic and septic injury differ in honey bee worker larvae and adults. *Arch. Insect Biochem. Physiol.* **2008**, *69*, 155–167. [CrossRef] [PubMed]

115. Cao, X.; Wang, Y.; Rogers, J.; Hartson, S.; Kanost, M.R.; Jiang, H. Changes in composition and levels of hemolymph proteins during metamorphosis of *Manduca sexta*. *Insect Biochem. Mol. Biol.* **2020**, *127*, 103489. [CrossRef]

116. Yu, S.; Luo, F.; Xu, Y.; Zhang, Y.; Jin, L.H. *Drosophila* Innate Immunity Involves Multiple Signaling Pathways and Coordinated Communication Between Different Tissues. *Front. Immunol.* **2022**, *13*, 905370. [CrossRef]

117. Vierstraete, E.; Verleyen, P.; Baggerman, G.; D'Hertog, W.; Van Den Bergh, G.V.; Arckens, L.; De Loof, A.; Schoofs, L. A proteomic approach for the analysis of instantly released wound and immune proteins in *Drosophila melanogaster* hemolymph. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 470–475. [CrossRef] [PubMed]

118. Frey, E.; Odemer, R.; Blum, T.; Rosenkranz, P. Activation and interruption of the reproduction of *Varroa destructor* is triggered by host signals (*Apis mellifera*). *J. Invertebr. Pathol.* **2013**, *113*, 56–62. [CrossRef] [PubMed]

119. Cohen, E. Fat body. In *Encyclopedia of Insects, Volume 3*; Resh, V.H., Cardé, R.T., Eds.; Academic Press: London, UK, 2003; p. 1295.

120. Piou, V.; Urrutia, V.; Laffont, C.; Hemptinne, J.L.; Vétillard, A. The nature of the arena surface affects the outcome of host-finding behavior bioassays in *Varroa destructor* (Anderson & Trueman). *Parasitol. Res.* **2019**, *118*, 2935–2943. [CrossRef]

121. Piou, V.; Vilarem, C.; Rein, C.; Sprau, L.; Vétillard, A. Standard Methods for Dissection of *Varroa destructor* Females. *Insects* **2022**, *13*, 37. [CrossRef]

122. Bournonville, L.; Askri, D.; Arafah, K.; Voisin, S.N.; Bocquet, M.; Bulet, P. Unraveling the *Bombus terrestris* Hemolymph, an Indicator of the Immune Response to Microbial Infections, through Complementary Mass Spectrometry Approaches. *Int. J. Mol. Sci.* **2023**, *24*, 4658. [CrossRef]

123. Askri, D.; Straw, E.A.; Arafah, K.; Voisin, S.N.; Bocquet, M.; Brown, M.J.F.; Bulet, P. Parasite and Pesticide Impacts on the Bumblebee (*Bombus terrestris*) Haemolymph Proteome. *Int. J. Mol. Sci.* **2023**, *24*, 5384. [CrossRef]

124. R Core Team. *R: A Language and Environment for Statistical Computing*; R Core Team: Vienna, Austria, 2015.

125. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis*; Springer International Publishing: Cham, Switzerland, 2016; ISBN 9783319242750.

126. Perez-Riverol, Y.; Bai, J.; Bandla, C.; García-Seisdedos, D.; Hewapathirana, S.; Kamatchinathan, S.; Kundu, D.J.; Prakash, A.; Fericks-Zipper, A.; Eisenacher, M.; et al. The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. *Nucleic Acids Res.* **2022**, *7*, 543–552. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.