Abstract Reports of honey bee population decline has spurred many national efforts to understand the extent of the problem and to identify causative or associated factors. However, our collective understanding of the factors has been hampered by a lack of joined up trans-national effort. Moreover, the impacts of beekeeper knowledge and beekeeping management practices have often been overlooked, despite honey bees being a managed pollinator. Here, we established a standardised active monitoring network for 5 798 apiaries over two consecutive years to quantify honey bee colony mortality across 17 European countries. Our data demonstrate that overwinter losses ranged between 2% and 32%, and that high summer losses were likely to follow high winter losses. Multivariate Poisson regression models revealed that hobbyist beekeepers with small apiaries and little experience in beekeeping had double the winter mortality rate when compared to professional beekeepers. Furthermore, honey bees kept by professional beekeepers never showed signs of disease, unlike apiaries from hobbyist beekeepers that had symptoms of bacterial infection and heavy Varroa infestation. Our data highlight beekeeper background and apicultural practices as major drivers of honey bee colony losses. The benefits of conducting trans-national monitoring schemes and improving beekeeper training are discussed.

Citation: Jacques A, Laurent M, EPILOBEE Consortium, Ribière-Chabert M, Saussac M, Bougeard S, et al. (2017) A pan-European epidemiological study reveals honey bee colony survival depends on beekeeper education and disease control. PLoS ONE 12(3): e0172591. https://doi.org/10.1371/journal.pone.0172591 Editor: Nicolas Chaline, Universidade de São paulo, BRAZIL Received: July 12, 2016; Accepted: February 7, 2017; Published: March 9, 2017 Copyright: © 2017 Jacques et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: The data belong to each Member State and the European Commission. They are available upon request to the Member State. Each Member State has a list of people who worked in EPILOBEE. The request for data should be addressed to them. The authors of this study may be contacted at eurl.bee@anses.fr. Funding: The financers of EPILOBEE were the European Commission, the Member States and the EU Reference Laboratory. The statistical analysis was co-funded by EFSA together with the EURL for honey bee health. GB acknowledges support from Fera and the Welsh Government (and MB, RG, GM, MP AW and IW in the Epilobee consortium) but did not have any additional role in the study design, data analysis, decision to publish, or preparation of the manuscript. The specific role of these authors rare articulated in author contribution section. Competing interests: This commercial affiliation Fera does not alter our adherence to PLOS ONE policies on sharing data and materials.

Introduction Honey bees are highly effective pollinators with an annual global contribution to crop productivity of € 147 million [1]. Recent decades have seen heightened concern about honey bee colony mortaility across the United States [2, 3], Asia [4] and Europe [5]. Whilst the global number of managed colonies has risen by about 45% over the last 60 years [6, 7], the seemingly unpredicable loss of honey bee colonies exacerbates the shortage of pollinators leading to concerns that pollination deficits may limit crop production [8]. Indeed in the US, it was shown that the quantity and quality of pollination services haved declined through time [9]. Honey bees are subject to many interacting pressures including pests, pathogens, pesticides and climate change—for a review see Vanbergen et al.[10]. National studies suggest different drivers of poor honey bee colony health depending on the geography. Pathogens were linked to poor honey bee health in the UK [11] and in Germany [12], pathogens and pesticides in Italy [13], while this was not found to be the case in Africa [14]. A first and recent Europe-wide experiment clearly indicated the presence of genotype-environment interactions originating from specific local adaptation of the honey bee populations [15]. Indeed, the most important result of this study was a significantly higher survivorship of the local genotypes compared to the non-local ones [16]. In the US, the diagnosed causes of overwintering mortality were different according to the beekeeper operation type: backyard beekeepers generally identified “manageable” factors (e.g., starvation, weak colony in the fall), while commercial beekeepers suggested environmental factors including pesticides [2]. In South-Africa, the effect of migratory beekeeping practices was found significant on colony losses, migratory beekeepers losing on average more colonies than stationary beekeepers [17]. For pollinators in general, other studies have shown the negative impact of land-use intensification [18], and interacting negative effects of human-induced changes on climate [19]. Despite all these studies, our collective understanding of honey bee health has been hampered by a lack of collaborative trans-national efforts following common protocols and the impacts of beekeeper knowledge and beekeeping management practices has often been overlooked, despite honey bees being a managed pollinator. The objective of the present study was to identify the key risk factors surrounding honey bee colony mortality through data from the first surveillance programme based on randomly selected participants and deploying standardised methods to monitor honey bee colony health, pests, diseases and management practices across 17 European countries.

Materials and methods Building a surveillance network The selection of Member States taking part to EPILOBEE, the methodology to randomly pick up the apiaries and the beekeepers, the evaluation of the protocols set up in each Member State are described in details in Chauzat et al. [20] and in the guidelines published by the European Reference Laboratory for Honey bee Health (EURL) [21, 22]. The EPILOBEE protocol is described briefly below. Surveillance was implemented during two consecutive years, between autumn 2012 and summer 2014. Three visits were set up in each Member State: before winter (autumn visit: V1), after winter (spring visit: V2) and during the beekeeping season (summer visit: V3). At each visit, beekeeping practices and clinical signs of the main honey bee diseases were recorded during field inspections using a standardised questionnaire. If colonies exhibited clinical signs of a disease, samples were collected for subsequent laboratory analyses. The main honey bee diseases clinically investigated were those listed for notification for intra-EU trade and import rules or for national eradication programmes at the European level [23, 24]: the fungal disease Nosemosis; the parasitic disease varroosis; the two main bacterial diseases affecting honey bee brood; the American foulbrood (AFB) and the European foulbrood (EFB); and a viral disease caused by the chronic bee paralysis virus (CBPV). An apiary was considered to be suffering from a disease after its diagnosis in one or more colonies. All case definitions were agreed between the EURL and the Member States (20). An agreed sampling protocol was circulated by the EURL after consultation with all participating Member States and is available upon request to the EURL [20]. Subsequently, each Member State organised training, managed the visits and stored the national data in an online database. For the second year, at least one third of the total beekeepers from the first year was renewed and for comparison the same methodology was utilised. Population representativeness was achieved through random sampling of apiaries (primary units) and bee colonies (secondary units) in each Member State considered to be representative of the national beekeeping population [20]. England & Wales (taken as one Member State) did not take part in the second year of the program. Study population studied and retained variables To calculate the colony mortality rates, only the apiaries with three consecutive visits and consistent data needed for mortality calculation were included. Data consistency was tested through several checks applied to the dataset including some editing to retain as much data as possible for analysis [25]. The consistency between various variables was checked; for exemple, the number of colonies randomly selected at the first visit and information reported at the second visit; the number of colonies randomly selected at the first visit and the size of the apiary; or the number of colonies owned by the beekeeper and the size of the apiary. After quality checking, 2 332 apiaries (out of 3 053) were retained for the analysis in the first year and 2 426 apiaries (out of 2 745) for the second year. Apiary was designated as the epidemiological unit. Among the available 138 variables in 12 tables, 36 were selected for the statistical analysis after a Delphi-like selection. In summary, these 36 retained variables were related to the beekeeper (age, activity and experience in beekeeping), the operation type (type of production i.e. honey, pollen, queens), the operation scale (the number of colonies owned by the beekeeper, the number of colonies in the visited apiary), the husbandry (honey bee subspecies [26], swarms and queens produced or bought), the disease state, whether depopulation had been observed in the apiary (clinical signs and mortality of colonies or honey bees observed before and during site visits) and the landscape surrounding the apiary. They were used as explanatory variables with the seasonal and the winter mortalities as response variables. With the seasonal mortality as response variable, the previous winter mortality was included as a supplementary explanatory variable. Statistical methods First, data sets from each year were considered independently and the 17 Member States were clustered according to honey bee colony mortality. Second, for each year considered again independently, the links between the seasonal and the winter mortalities were explored. To work with a complete dataset, missing data were completed using the data imputation method based on multiple correspondence analysis. Finally, risk factors of mortality were sought on the overall study including both years of data by means of generalised linear models (see below). Clustering the Member States according to yearly mortality Yearly honey bee colony mortality was studied pooling together the winter mortality and the seasonal mortality on 2 332 apiaries for the first year and 2 426 apiaries for the second year. The 17 Member States were clustered according to mortality patterns (i.e. groups with high mortality rates versus lower mortality rates) using vectors composed of two values (for the winter mortality and for the seasonal mortality rates). The 17 vectors were subsequently analysed through the hierarchical clustering of observations on the main principal components function (HCPC) from the FactoMineR R package [27]. This resulted in producing groups of Member States with similar annual mortalities for the first and second year separately. Link between the winter and the seasonal mortalities This specific analysis was performed for EPILOBEE first year and second year separately. Spearman correlation coefficients were calculated taking into account two variables for each year: a first one corresponding to winter mortality rates and a second one corresponding to seasonal mortality rates. Risk factors related to honey bee colony mortality (both years combined) Considering the missing data in data sets from each year, 1 139 (out of 2 332) apiaries should have been kept in the analysis for the first year, and 1 020 (out of 2 426) for the second year. To take into acount a complete dataset, an imputation method was performed to process the missing data separately for each year, based on multiple correspondence analysis (MCA) [28] through the functions estim_ncpMCA and imputeMCA from the missMDA R package [29]. After imputation, data sets from each year were merged and the new global dataset was used in the following analysis. The winter and seasonal honey bee colony mortalities (response variables) were considered as count numbers with quasi-Poisson distributions. In the subsequent analysis, the models used were the quasi-Poisson generalised linear models. Calculations of p-values (associated to the likelihood ratio) were implemented for the winter mortality and the seasonal mortality for each explanatory variable (glm function R package). The 0.20 threshold was used to select the variables kept in the study: the ones with a p-value superior to the threshold were removed from the study, the others variables being incorporated in the subsequent steps. To avoid multicollinearity in Poisson regressions, variables from each theme (e.g., diseases, management) were summarized in a single categorical latent variable through MCA and classification using the FactoMineR package [30] as detailed in Jacques et al. [25]. A unique new synthetic explanatory variable resulting from clustering and composed of seven categories was included in the generalised linear model as a fixed effect. The Country and the year variables were included as random effects. For each mortality, a mixed multivariate Poisson regression model was performed using the glmer function from the lme4 R package with a log link [31]. Thereafter, contrasts were used to compare mortalities of each cluster of the synthetic explanatory variable to the other ones. At the end of the classification, the seven clusters of the synthetic variable were built according to eight criteria: the typology of the beekeeper (hobby, part time or professional), the size of the apiaries and of the colonies according to four classes: size 1 = operation having less than 50 colonies in total and apiary having less than 20 colonies; size 2 = operation having less than 50 colonies and apiary having between 21 and 50 colonies; size 3 = operation having between 51 and 300 colonies and apiary having more than 50 colonies; and size 4 = operation with more than 300 colonies and apiaries having more than 50 colonies. The age of the beekeeper was recorded in four groups (less than 30 years of age, between 30 and 45 years, between 45 and 65 years and more than 65 years of age). The involvement in a cooperative treatment against Varroa and the observation of clinical diseases in colonies visited and the subsequent confirmation in laboratory were also considered in the synthetic variable. The actions implemented at the operation level to improve the quality of production were reported under the section Apiary management: if the beekeeper produced his/her own queens and swarms or whether he/she bought some; if the management of the bee population by acquiring swarms and queens was intended to increase the production, to compensate the bad health conditions or to maintain the livestock; which subspecies of honey bees were used in the colonies visited; which environment surrounded the apiary under study according to the beekeeper; what kind of main production was targeted by the beekeeper; if the colonies were merged; and which honeyflow was targeted by the migration. The bee subspecies were reported by the beekeepers. There was no morphological nor molecular analysis conducted on honeybees to assess the subspecies. The environment surrounded the apiary under study was also reported by the beekeepers. No landscape measures were recorded to assess this criterion. The section on the beekeeper background gathered the experience in beekeeping (training, qualification, how long he/she is a beekeeper) and some practices (using a apiarist book, be part of a beekeeping organisation). The last section noted the presence of any health events in the apiary before the start of the project. For seasonal mortality only, a supplementary criteria was added to account for the winter colony mortality previously recorded.

Discussion Our results demonstrate that overwinter losses ranged between 2% and 32%, and that high summer losses were likely to follow high winter losses. Multivariate Poisson regression models revealed that hobbyist beekeepers with small apiaries and little experience in beekeeping had double the winter mortality rate when compared to professional beekeepers. Furthermore, honey bees kept by professional beekeepers never showed signs of disease during the visits, unlike apiaries from hobbyist beekeepers that had often symptoms of bacterial infection and/or heavy Varroa infestation observed during the visits. Our results show that the main factors protecting honey bee colonies are beekeeper background and practices. More efforts are needed in beekeeper training to promote good beekeeping practices and achieve early identification of clinical signs of disease. During the first year of the study there was no clear geographical pattern for the distribution of yearly mortality rates (including both winter and seasonal rates). Conversely, during the second year, the Member States with similar mortality rates were distributed according to geographical patterns, grouped in patches. The winter 2012–2013 was particularly long and cold across Europe. This was not the case for the following winter (2013–2014) which was rather mild [32]. Mild winters are known to favour the build-up of Varroa mite populations. According to a longitudinal study run in Europe, the colony location strongly influence autumn mite loads and viral prevalence [15]. Multi-year longitudinal monitoring of colony health can identify region-specific risk factors associated with colony mortality [33]. Once these factors have been identified, they can inform on management and research priorities including pesticide uses [33]. It has been recently suggested that different honey bee populations may have developed their own specific resistance mechanisms tailored to match the challenge of the environment they are located in, including the locally prevailing combination of pathogens and pathogen variants [16]. Similar results were shown in England and Wales for AFB exhibiting significant spatial aggregation at distances from 10 to 30 km [34]. In the field of ecotoxicology, sublethal effects of neonicotinoid pesticides were also modified in magnitude by environmental interactions specific to the landscape and time of exposure events [35]. These observations advocate the use of honey bee health surveys to be conducted at the landscape scale, and to include environmental factors. However, the distribution of the Member States with similar mortalities observed during the first year of the study without any geographical pattern may indicate that weather could be an overwhelming factor that can mask regional effects. In the US data showed that beekeepers in northern states lose more colonies than those in southern states [2, 33]. Therefore the role of climate in winter colony mortality could be further elucidated by repeating the surveys over serveral years. The mortality rates reported in EPILOBEE were relatively low compared to rates reported in the US which are consistenly around one-third for each winter [3]. It is not known whether these discrepancies could be attributed to sampling issues, analysis, biological pressures, landscapes or management factors, or any combination of the above. Commercial beekeeper in the US typically transport their hives thousands of miles per year, which is exceedingly rare in Europe. Although some studies have already highlighted the stress induced by these transportations, how these factors influence the various mortality rates remains largely unknown at colony level [36, 37]. We also present the first evidence of a relationship between winter and seasonal mortalities. Although winter colony mortality is commonly used to quantify losses in a standardized way, the colony mortality during a full year should also be taken into consideration to estimate population level trends. Colony losses are highly dependant on beekeeper management which includes how colonies are taken care (disease detection) but also which honey flow is targeted and therefore potential specific pesticide exposure [38]. Thus, future surveillance studies should encompass information on land use and pesticide analyses. The statistical analysis of our data from the first year has shown that size of the operation, apiary and the clinical disease (varroosis, AFB and nosemosis) observed during the previous autumn were significant risk indicators of overwinter honey bee colony survival [20]. In the present study, our findings find a similar result when considering two study years. The study variables were grouped into seven summary clusters. The clusters with the highest and the lowest mortality were distinguished by occurrence of diseases and beekeeper background. When the data set was extended to include both years, only prior observations of varroosis for the winter mortality and AFB for the summer mortality were found to be significant risk factors. Beekeepers with better knowledge on disease detection and management (specifically for varroosis and AFB), applying earlier prophylactic measures and good beekeeping practices (e.g. preparation of colonies for winter) had lower mortality rates than others (Tables 1 & 2). Recently, historical data have been used to explore the link between honey bee colony mortality and other risk factors such as changes in the political and socioeconomic systems [39]. Using FAO data, it was shown that the most dramatic decline in the number of colonies in Europe was statistically associated with the collapse of the socialistic regimes in the Eastern European countries after 1989. In contrast, the arrival of the parasitic mite V. destructor in the early eighties in Europe had no detectable effects on the number of managed colonies in the full European data set. This does not show that Varroa does not kill honey bee colonies. However, it shows how beekeepers efficiently adapt their operation to comply with the challenges set by pests even as lethal as V. destructor [39]. Our results are in accordance with the latter study showing that beekeepers tend to have practices that compensate for the colony losses. Indeed, during both years of the study, the number of swarms bought was positively linked with both mortalities (the more swarms were bought, the higher the colony mortality–results not shown). The statistical link was also significant between merged colonies and colony losses. This might reflect the additional work needed by a beekeeper to overcome colony losses: beekeepers got new swarms and merged colonies to maintain their livestock at an acceptable level in quantity and quality. In future studies, the assessment of beekeeper additional work to compensate for losses or weakness of colonies should be better quantified. The limits of our descriptive epidemiological protocol should be taken into account when drawing conclusions. Indeed, any hypotheses expressed in this paper should be fully studied in dedicated experimental protocols to confirm the risk factors and clarify any potential causality [33]. Our study offers observational evidence to suggest the importance of beekeeper training and education. These results must be seen as preliminary until confirmed by direct experimental means.

Conclusion Our results show that the main factors protecting honey bee colonies are beekeeper background and practices. More efforts are needed in beekeeper training to promote good beekeeping practices and achieve early identification of clinical signs of disease. Considerable variation of colony losses exist across different Member States and between years. Climate conditions might have a strong effect on colony mortality during the whole year, requiring long term surveillance study to overcome the weather factor. Data from descriptive survey such as EPILOBEE should be used to set up dedicated protocols to study further targeted hypotheses. The promotion of regional scale studies of local practices should be encouraged. Further to this work, the causes of colony losses should be investigated by conducting studies on specific issues as potential causes of honey bee losses, for example case-control studies that include pesticide analyses and landscape recording.

Acknowledgments The authors wish to thank the beekeepers who partipated in this important field work. The complete membership of the EPILOBEE Consortium Belgium: De Graaf Dirk (Ghent University, Department of Physiology, Laboratory of Zoophysiology); Méroc Estelle (National Reference Laboratory for honeybee diseases CODA-CERVA-VAR); Nguyen Bach Kim (Ulg, Faculté Gembloux Agro-Bio Tech); Roelandt Sophie (National Reference Laboratory for honeybee diseases CODA-CERVA-VAR); Roels Stefan (CODA-CERVA-VAR); Van der Stede Yves (CODA-CERVA-VAR). Denmark: Tonnersen Tina (National Reference Laboratory for honeybee diseases Aarhus University); Kryger Per (NRL Aarhus University). Estonia: Jaarma Kärt (Estonian Veterinary and Food Board); Kuus Merle (Estonian Veterinary and Food Laboratory); Raie Arvi (Estonian Veterinary and Food Board). Finland: Heinikainen Sirpa; Pelkonen Sinikka, Vähänikkilä Nella (Finnish Food Safety Authority EVIRA, Veterinary Bacteriology Research Unit, Kuopio). France: Andrieux Christophe (DDPP du Cantal); Ballis Alexis (Chambre d'Agriculture du Haut-Rhin); Barrieu Guy (DDPP des Bouches du Rhône); Bendali Fatah (Direction Générale de l'Alimentation); Brugoux Corinne (Groupement de Défense Sanitaire du Cantal); Franco Stéphanie (National Reference Laboratory for honeybee diseases, Anses, Sophia Antipolis); Fuentes Anne-Marie (Groupement de Défense Sanitaire de la Drôme); Joel Alain (DDPP Finistère); Layec Yves (Groupement de Défense Sanitaire Apicole du Finistère); Lopez Jacqueline (DDPP Indre et Loire); Lozach André (Groupement de Défense Sanitaire Apicole du Finistère); Malherbe-Duluc Laure (Groupement de Défense Sanitaire Indre et Loire); Mariau Viviane (DDPP Indre et Loire); Meziani Fayçal (Direction Générale de l'Alimentation); Monod Denis (Groupement de Défense Sanitaire Apicole des Bouches du Rhône); Mutel Sébastien (DDCSPP Haut-Rhin); Oesterle Eric (Groupement de Défense Sanitaire Indre et Loire); Orlowski Muriel (DDPP de la Drôme); Petit Manuel (DDPP Finistère); Pillu Patricia (DDPP du Cantal); Poret Florence (Groupement de Défense Sanitaire du Cantal); Viry Alain (Laboratoire d'Analyses du Jura). Berg Stefan (Bavarian State Institute for Viticulture and Horticulture, Bee Research Center, Veitshöchheim); Büchler Ralph (LLH Bieneninstitut Kirchhain); de Craigher Doris (University of Hohenheim, Apicultural State Institute, Stuttgart); Genersch Elke (Institute for Bee Research, Hohen Neuendorf); Kaatz Hans-Hinrich (University of Halle-Wittenberg, Zoology Dept., Halle); Meixner Marina D. (LLH Bieneninstitut Kirchhain); von der Ohe Werner (LAVES Institut für Bienenkunde, Celle); Otten Christoph (Dienstleistungszentrum Ländlicher Raum, Fachzentrum Bienen und Imkerei Mayen); Rosenkranz Peter (University of Hohenheim, Apicultural State Institute, Stuttgart);Schäfer Marc O. (Institute of Infectology, Friedrich-Loeffler-Institut, Greifswald—Insel Riems); Schroeder Annette (University of Hohenheim, Apicultural State Institute, Stuttgart). Greece: Agianiotaki Eirini (Centre of Veterinary Institute of Athens); Arfara Stamatia (Centre of Veterinary Institute of Athens); Boutsini Sofia (Centre of Veterinary Institute of Athens); Giannoulopoulou Marianthi (Regional Veterinary Laboratory of Heraclio); Hondrou Varvara (Regional Veterinary Laboratory of Mytilini); Karipidou Sonia (Regional Veterinary Laboratory of Kozani); Katsaros Dimitrios (Regional Veterinary Laboratory of Chalkis); Katzagiannakis Aristomenis (Regional Veterinary Laboratory of Heraclio); Kiriakopoulos Antonios (Regional Veterinary Laboratory of Mytilini); Oureilidis Konstantinos (Regional Veterinary Laboratory of Kavala); Panteli Aspasia (Centre of Veterinary Institute of Athens); Pantoleon Fotis (Regional Veterinary Laboratory of Tripoli); Papagianni Zoi (Centre of Veterinary Institute of Athens); Papalexiou Eleni (Centre of Veterinary Institute of Athens); Perdikaris Socrates (Ministry of Rural Development and Food); Prapas Athanasios (Centre of Veterinary Institute of Athens); Siana Panagiota (Regional Veterinary Laboratory of Tripoli); Skandalakis Ioannis (Regional Veterinary Laboratory of Chania); Stougiou Despoina (Centre of Veterinary Institute of Athens); Tomazinakis Ioannis (Regional Veterinary Laboratory of Chania); Tsali Eleftheria (Regional Veterinary Laboratory of Larisa); Tseliou Evgenia (Regional Veterinary Laboratory of Kerkyra); Tsiplakidis Achillefs (Regional Veterinary Laboratory of Kavala); Tsompanellis Efstratios (Regional Veterinary Laboratory of Mytilini); Vamvakas Giorgos (Regional Veterinary Laboratory of Kozani); Varvarouta Vasiliki (Regional Veterinary Laboratory of Heraclio); Vourvidis Dimitris (Ministry of Rural Development and Food). Hungary: Dán Ádám, Daróczi Gyöngyi, Láng Mária, Papp Melitta, Paulus Petra–Deákné, Pupp Eszter (National Food Chain Safety Office, Veterinary Diagnostic Directorate -NRL for bee health); Szaló Márta (Ministry of Agriculture, Food Chain Control Department); Tóth Ádám, Zséli Szilvia (NRL for bee health). Italy: Bressan Gianluigi (Ulss22, Bussolengo); Cerrone Anna (IZS del Mezzogiorno); Formato Giovanni (IZS delle Regioni Lazio e Toscana); Granato Anna (IZS delle Venezie); Lavazza Antonio (IZS della Lombardia e dell'Emilia Romagna); Macellari Piero (Asl Umbria); Marcello Paola (ASL Sassari); Ghittino Claudio (IZS dell'Umbria e delle Marche); Maroni Ponti Andrea (Ministero della Salute); Possidente Rosaria (IZS del Piemonte Liguria e Valle d'Aosta); Mutinelli Franco (IZS delle Venezie); Nassuato Claudia (Regione Lombardia); Pintore Antonio (IZS della Sardegna); Ricchiuti Luciano (IZS dell'Abruzzo e del Molise); Ruocco Luigi (Ministero della Salute); Salvaggio Antonio (IZS della Sicilia); Troiano Pasquale (IZS di Puglia e Basilicata); Voltini Barbara (Regione Toscana). Latvia: Avsejenko Jelena, Ciekure Elïna, Deksne Gunita (Institute of Food safety, Animal Health and Environment, "BIOR"); Eglïte Ineta (Latvian Beekeepers Association); Granta Rita (BIOR); Olševski Edvïns (Food and Veterinary Service of the Republic of Latvia); Rodze Ieva (BIOR); Stinka Madara (Food and Veterinary Service of the Republic of Latvia). Lithuania: Sirutkaityte Rasa, Siriukaitis Sigitas (The State Food and Veterinary Service, Animal Health and Welfare department); Poland: Bober Andrzej (National Veterinary Research Institute, Pulawy); Jażdżewski Krzysztof (General Veterinary Inspectorate, Warsaw); Pohorecka Krystyna, Skubida Marta, Zdańska Dagmara (National Veterinary Research Institute, Pulawy). Portugal: Amador Maria Rita Ramos, Freitas Susana, Quintans Sofia (Direção Geral de Alimentacão e Veterinaria); Santos Patricia Tavares (Direção Geral de Alimentacão e Veterinaria). Slovakia: Březinová Nicole, Brtková Andrea, Čuvalová Zuzana, Filipová Miriam, Jurovčiková Júlia, Kantíková Miriam, Kubicová Zuzana (State Veterinary and Food Institute, Dolny Kubin); Papierniková Erika (State Veterinary and Food Administration of the Slovak Republic), Šulejová Lucia (State Veterinary and Food Institute, Dolny Kubin); Toporčák Juraj (The University of veterinary medicine and pharmacy in Košice). Spain: Ares Cenador Carmen Maria (Consejería de Agroganadería y Recursos Autóctonos del Principado de Asturias); Ariza Javier (Consejería de Agricultura de La Junta de Comunidades de Castilla La Mancha); Berná Serna Nieves (GVA, Consejería de la Presidencia de Agricultura, Pesca, Alimentación y Agua); Cabeza Núñez Amparo (Consejería de Agricultura y Pesca de la junta de Andalucía); Casasempere Cascales Jorge (GVA, Consejería de la Presidencia de Agricultura, Pesca, Alimentación y Agua); Cid González Carlos (Subdireción Xeral de Gandería—Consellería do Medio Rural e do Mar—Xunta de Galicia); Corzán Ripoll Jose Manuel (Consejería de Agricultura, Ganadería y Medio Ambiente, Diputación General de Aragón); De Abajo Domingo Miguel Ángel (Consejería de Agricultura y Ganadería de la Junta de Castilla y León); Díaz Rey Roberto (Subdireción Xeral de Gandería—Consellería do Medio Rural e do Mar—Xunta de Galicia); Esteban Royo Ángel (Consejería de Agricultura, Ganadería y Medio Ambiente, Diputación General de Aragón); Fernández Somalo Pilar (SG de Sanidad, Higiene Animal y Trazabilidad del M° Agricultura, Alim. y Medio Ambiente); García Pascualvaca Alejandra (Consejería de Agricultura y Pesca de la junta de Andalucía); González Breña Carlos (Consejería de Agricultura, Desarrollo Rural, Medio Ambiente y Energía, J. de Extremadura); Mínguez Gonzalez Olga (Consejería de Agricultura y Ganadería de la Junta de Castilla y León); Oñate Maria Luisa (Consejería de Agricultura, Ganadería y Medio Ambiente, Diputación General de Aragón); Oteiza Orradre Pedro (Dpto. de Desarrollo Rural, Industria, Empleo y Medio Ambiente, Diputación Foral Navarra); Pérez Cobo Iratxe (SG de Sanidad, Higiene Animal y Trazabilidad del M° Agricultura, Alim. y Medio Ambiente); Plaza Pérez Margarita (Consejería de Agricultura y Agua de la Región de Murcia); Puy Pitarque D. Juan Ramón (Departamento de Agricultura, Pesca y Alimentación del Gobierno Vasco); Riol Guinea Rubén (Consejería de Agricultura y Ganadería de la Junta de Castilla y León); Romero González Luis José (SG de Sanidad e Higiene Animal y Trazabilidad del M° Agricultura, Alim. y Medio Ambiente); Soldevilla Yanguas Jose Fernando (Consejería de Agricultura, Ganadería y Medio Ambiente, Comunidad Autónoma de la Rioja); Soler i Barrasús Mercè (Dept. Agricultura, Ramaderia, Pesca, Alimentación i Medi Natural, Generalitat de Catalunya); Soriano González Mario (Consejería de la Presidencia de Agricultura, Pesca, Alim. y Agua, Generalitat Valenciana); Vigo López Virginia (Consejería de Agricultura, Ganadería, Pesca y Aguas del Gobierno de Canarias); Villarta Rivas José Luis (Consejería de Agricultura de La Junta de Comunidades de Castilla La Mancha). Sweden: Fabricius-Kristiansen Lotta (Swedish Board of Agriculture); Forsgren Eva (Department of Ecology, Swedish University of Agricultural Sciences). UK: Brown Mike, Budge Giles, Grant Ruth, Marris Gay, Powell Michelle, Wattam Andy, Whiting Ilex (The National Bee Unit (NBU)—Food and Environment Research Agency (FERA). EURL staff: Cauquil Laura, Rivière Marie-Pierre (Unit of Honey bee Pathology); Garin Emmanuel (Unit of coordination and support to surveillance).

Author Contributions Conceptualization: MRC PH MPC. Data curation: ML MRC PH MPC. Formal analysis: AJ ML MS SB PH MPC. Funding acquisition: ML MRC SB PH MPC. Investigation: AJ ML MRC MS SB PH MPC. Methodology: AJ ML MS SB PH MPC. Project administration: ML MRC PH MPC. Resources: AJ ML MRC MS SB PH MPC. Software: AJ MS ML SB. Supervision: ML MRC PH MPC. Validation: AJ ML MS SB PH MPC. Visualization: AJ ML MS SB GEB MPC. Writing – original draft: AJ ML MS SB MPC. Writing – review & editing: AJ ML MS SB GEB MPC.