SHORT COMMUNICATION

 

First data on the prevalence and distribution of pathogens in bumblebees (Bombus terrestris and Bombus pascuorum) from Spain

 

Clara Jabal-Uriel

Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF). Centro Agrario de Marchamalo, Marchamalo, 19180 Guadalajara, Spain.

Raquel Martín-Hernández

Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF). Centro Agrario de Marchamalo, Marchamalo, 19180 Guadalajara, Spain.

Fundación Parque Científico y Tecnológico de Albacete, Instituto de Recursos Humanos para la Ciencia y la Tecnología (INCRECYT), Albacete, Spain.

Concepción Ornosa

Universidad Complutense de Madrid, Facultad de Biología, Dept. Zoología y Antropología Física, 28040 Madrid, Spain.

Mariano Higes

Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF). Centro Agrario de Marchamalo, Marchamalo, 19180 Guadalajara, Spain.

Eduardo Berriatúa

Universidad de Murcia, Facultad de Veterinaria, Dept. Sanidad Animal, 30100 Murcia, Spain.

Pilar De la Rúa

Universidad de Murcia, Facultad de Veterinaria, Dept. Zoología y Antropología Física, 30100 Murcia, Spain.

 

Abstract

Bumblebees provide pollination services not only to wildflowers but also to economically important crops. In the context of the global decline of pollinators, there is an increasing interest in determining the pathogen diversity of bumblebee species. In this work, wild bumblebees of the species Bombus terrestris and Bombus pascuorum from northern and southern Spain were molecularly screened to detect and estimate prevalence of pathogens. One third of bumblebees were infected: while viruses only infected B. pascuorum, B. terrestris was infected by Apicystis bombi, Crithidia bombi and Nosema bombi. Ecological differences between host species might affect the success of the pathogens biological cycle and consequently infection prevalence. Furthermore, sex of the bumblebees (workers or males), sampling area (north or south) and altitude were important predictors of pathogen prevalence. Understanding how these factors affect pathogens distribution is essential for future conservation of bumblebee wild populations.

Additional key words: pollinators; pathogen dispersion; PCR; Apicystis bombi; Crithidia bombi; Nosema bombi.

Abbreviations used: AKI (Acute bee paralysis, Kashmir bee and Israeli acute paralysis viruses complex); BQCV (Black Queen Cell Virus); DWV (Deformed Wing Virus); LSV (Lake Sinai Virus); PCR (Polymerase Chain Reaction) RT-PCR (Real Time PCR).

Authors’ contributions: Conceived and designed the experiments: CO, RMH, MH and PDLR. Performed the experiments: CJU and RMH. Performed the statistical analyses: CJU and EB. Contributed taxonomic analysis: CO. Wrote the paper: CJU, RMH, CO, MH, EB, PDLR.

Citation: Jabal-Uriel, C.; Martín-Hernández, R.; Ornosa, C.; Higes, M.; Berriatúa, E.; De la Rúa, P. (2017). Short communication: First data on the prevalence and distribution of pathogens in bumblebees (Bombus terrestris and Bombus pascuorum) from Spain. Spanish Journal of Agricultural Research, Volume 15, Issue 1, e05SC01. https://doi.org/10.5424/sjar/2017151-9998

Received: 25 May 2016. Accepted: 8 Feb 2017.

Copyright © 2017 INIA. This is an open access article distributed under the terms of the Creative Commons Attribution (CC-by) Spain 3.0 License.

Funding: Spanish Ministry of Economy and Competitiveness and FEDER funds (BIOBOMBUS (CGL2012-34897); INIA, Madrid (E-RTA2014-00003-C03 and 01); Regional Government of Murcia (19908/GERM/2015, Fundación Séneca). PDLR is presently member of and receives support from COST Action FA1307 (Sustainable pollination in Europe: joint research on bees and other pollinators (SUPER-B; http://www.cost.eu/COST_Actions/fa/Actions/FA1307).

Competing interests: None of the authors of this paper has a financial or personal relationship with other people or organizations that could inappropriately influence or bias the content of the paper. The authors declare that they have no competing interests.

Correspondence should be addressed to Pilar De la Rúa: pdelarua@um.es


 

CONTENTS

Abstract

Introduction

Material and methods

Results and discussion

Acknowledgments

References

IntroductionTop

Pollination is important for ecological processes and worldwide agricultural productivity (Potts et al., 2016). Several causes of global pollinator decline have been described, including climate change, pathogens, habitat loss and pesticides (Goulson et al., 2015).

Among pollinators, bumblebees (Bombus Latreille, 1802) provide pollination services not only to wildflowers but also to economically important crops; for this, some species are reared and commercialized for greenhouse production (Ortiz-Sánchez, 1992; Ornosa, 1996; Velthuis & van Doorn, 2006). Several pathogens (Crithidia Léger, 1902: Trypanosomatidae, Apicystis (Lipa & Triggiani, 1996): Neogregarine, Nosema Nägueli, 1857: Microsporidia, and viruses) are related to the decline of bumblebee species, and their prevalence has been monitored elsewhere but not in the Iberian Peninsula (Williams & Osborne, 2009; Cameron et al., 2011; Meeus et al., 2011; Gallot-Lavallée et al., 2016). To address this knowledge gap, samples of the most widely distributed species Bombus terrestris (Linnaeus, 1758) and Bombus pascuorum (Scopoli, 1763) were collected to map pathogen prevalence. We hypothesized that infection rates differ between northern and southern Spain, between species and between sexes.

Material and methodsTop

Worker and male bumblebees were sampled from 19 locations in Spain: either from a southern mountain region (Sierra Nevada National Park, Granada) at altitudes between 1041 and 2752 m a.s.l., where conservation and protection measures of fauna and flora are implemented, or in northern human-altered landscapes (mainly edges of agricultural fields, river banks and urban gardens) between 63 and 1031 m a.s.l (Table 1).


Table 1. Sampling information of Spanish bumblebees (M = male and F = female worker). Mean temperatures during August at each location were taken from Ninyerola et al. (2007).

Bombus terrestris and B. pascuorum species remain abundant in the study area, and are not considered threatened according to the International Union for Conservation of Nature (Rasmont et al., 2015). Individuals were sampled with nets while foraging on flowers during August 2013 with dry weather (above 18 ºC) and clear weather conditions, and preserved in 100% ethanol at – 4 ºC.

Every sampled bumblebee was dissected to remove the genitalia and only the remaining abdominal tissues (gut and fat body) were used for DNA extraction following the Chelex method (Scriven et al., 2013). RNA was isolated with RNase-Free DNase kit (Quiagen) and transformed in cDNA with QuantiTec reverse transcription kit (Quiagen) following manufacturer’s instructions.

Bumblebee species identification was confirmed through morphological and sequence analyses (Murray et al., 2008).

To amplify the fragment of the 18S rDNA gene of Crithidia spp. and Apicystis spp., the reaction mixture consisted of a triplex PCR reaction including Apidae primers (as DNA extraction and amplification control for Bombus species) (Meeus et al., 2010). To amplify Nosema bombi (Fantham & Porter, 1914) the reaction mixture consisted of a duplex PCR also including the Apidae primers. For Nosema apis (Zander 1909) and Nosema ceranae (Fries et al., 1996) BioTools amplification plates were used (Martín-Hernández et al., 2007), consisting of a prepared microplate with gelled enzymes and primers to which only water and the DNA sample need to be added. To amplify viruses, the reaction mixtures were performed following Ravoet et al. (2013) for Lake Sinai Virus (LSV), Francis & Kryger (2012) for AKI (Acute bee paralysis, Kashmir bee and Israeli acute paralysis viruses complex) and Chantawannakul et al. (2006) for Deformed Wing Virus (DWV) and Black Queen Cell Virus (BQCV). For detecting BQCV and DWV virus, cDNA amplification was performed with a RT-PCR, using as internal control Apidae primers from Meeus et al. (2010). Detailed information about primer and probe sequences are provided in (Table 2).


Table 2. Primers (F: forward; R: reverse; T: probe), oligonucleotide sequences of the primers, annealing temperature (Ta), expected size of the fragments (bp: base pairs) used for PCR and RT-PCR amplification of the genes from the different organisms.

PCR reactions were run for 10 min at 95 ºC and followed by 35 cycles of denaturing at 95 ºC for 30 s, specific annealing temperature for 30 s, extending at 72 ºC for 45 s, and final extension of 7 min at 72 ºC. Positive (DNA of the pathogen) and negative (water instead of DNA) controls were included in each PCR reaction.

To verify that amplicons (those from the bumblebees for species confirmation, and from the pathogens) corresponded to the target organisms, positive samples were sequenced and compared with sequences uploaded in the NCBI through the BLAST tool in MEGA v6 (Tamura et al., 2013).

Prevalence of infection was estimated as a percentage of the number of infected divided by the total number of individuals. Proportions were compared using Yates-corrected chi-squared test, or when appropriate, Fisher´s exact test. Medians were compared with the non-parametric Kruskall-Wallis test. The relationship between infection with any pathogen or with Crithidia spp. and Apicystis spp. and explanatory variables was further explored using mixed logistic regression analysis to correct for spatial autocorrelation. Infection status (infected or non-infected) was the response variable. Explanatory variables including sex (female worker or male), environmental temperature, altitude, geographical distribution (north and south) and bumblebee species (B. terrestris and B. pascuorum) were fitted as fixed effects, and location was included as a random effect. A backward selection strategy was used to select explanatory variables, starting with a model that included all variables. Due to the high correlation between explanatory variables (for example altitude and sampling area), different combinations of models with two variables were finally compared using Akaike´s Information Criteria to select the model with the lowest value. Models were estimated using the maximum likelihood estimation method and significance was set at p<0.05 for a double-sided test.

Results and discussionTop

A total of 115 bumblebees were sampled: 83 B. terrestris and 32 B. pascuorum. Bumblebee species abundance differed between regions: 75% B. terrestris came from the south, while 68% B. pascuorum were from the north (p<0.05). The sex ratio was similar for B. terrestris individuals, and 72% of B. pascuorum individuals were females (p<0.05).

Alignment of the amplicons with sequences FN546181 and FN546182 (Meeus et al., 2010) confirmed the presence of Crithidia bombi (Lipa & Triggiani 1988) and Apicystis bombi (Liu et al., 1974). Pathogens distribution is described in (Table 3). Both protozoa were present in 37% B. terrestris and in none B. pascuorum. The estimated prevalence (95% CI) for B. terrestris was 15.7% (7.8-23.5%), for C. bombi and 24.1% (14.9-33.3%) for A. bombi, respectively. Two individuals were infected by both pathogens. Prevalence of both protozoans was significantly higher in males than in female workers, and A. bombi was more frequent in the south and C. bombi in lower altitude areas (p<0.05). Logistic regression models corroborated the relationship between protozoa infection, sex and area (north and south) or altitude. The risk of infection was greatest in middle range altitude (p<0.05).


Table 3. Pathogen prevalence: percentage of PCR-positive individuals and 95% confidence interval (CI) in bumblebees from Spain per species, area (N = northern and S = southern provinces) and sex (F = female worker and M = male). BQCV refers to Black Queen Cell Virus and DWV to Deformed Wing Virus.

Of the three Nosema species, only one B. terrestris male from the north was infected with N. bombi. Viruses were only found in B. pascuorum and with low prevalence: two workers were infected with DWV (one southern and one northern), and three were infected with BQCV (one male and two workers from the north). No bumblebees were infected with LSV or AKI.

Overall, pathogens (protozoa, N. bombi and bum­ble-bee viruses) were detected in 34% of sampled bumblebees; 16% of B. pascuorum and 37% of B. terrestris (p<0.05). Pathogen prevalence was significantly greater in males than in females (p<0.05), and did not differ according to environmental temperature or altitude, although it was numerically greatest in B. terrestris males, in low altitudes and warmer places (p>0.05). There was no evidence for significant variation in infection risk according to location.

This study shows that pathogen prevalence in two bumblebee species in Spain varies depending on host and pathogen species, sex and area, although in general the level of pathogens was low. Prevalence of infection was higher in B. terrestris than in B. pascuorum, however, the former had protozoa and a microsporidium infection, and the latter only carried viruses.

A. bombi showed the highest prevalence and appeared exclusively in southern B. terrestris. This neogregarine has been recently detected in other Bombus species (Gamboa et al., 2015), but its dispersion has been reported mainly due to managed bumblebees such as B. terrestris (Graystock et al., 2014). In Spain, the extensive use of this species for greenhouse tomato pollination may have contributed to the spread of this pathogen.

The second most common parasite in B. terrestris was C. bombi as in UK (Goulson et al., 2012). Trypanosomatids have rapid adult-adult transmission linearly related to host density (Goulson et al., 2012), occurring either in sympatric or allopatric populations (Meeus et al., 2011). They can be transmitted by direct contact or indirectly by dispersal onto flowers (Graystock et al., 2015). B. terrestris has larger nests than B. pascuorum (Goulson et al., 2012), which could facilitate faster transmission. In this study, males showed higher infection prevalence than workers, which may be explained by biological reasons: males live mainly outside the nests, and are more exposed to pathogen spillover from other pollinator species (Goulson, 2010).

Regarding microsporidia, N. bombi was the only species detected in B. terrestris. This pathogen has also been found in the UK (Goulson et al., 2012), but in Spain has not yet been reported. Notably, N. ceranae was not present in our sample, despite its high prevalence in bumblebees from UK (Graystock et al., 2013) and in honeybees in Spain (Martín-Hernández et al., 2007).

As in other studies (Singh et al., 2010; Meeus et al., 2011), LSV and AKI-complex viruses did not appear in Spanish bumblebees albeit BQCV and DWV have been found with a low prevalence only in B. pascuorum. DWV is highly widespread among bumblebees in Europe, and has been found in B. terrestris and B. pascuorum living nearby honeybee colonies in Germany (Genersch et al., 2006). Given that we did not account for honeybee colony density in the sampling areas, a cause-effect of this low prevalence cannot be assessed.

Two individual bumblebees were infected by more than one pathogen, suggesting individual differences regarding susceptibility to parasitic infection, or that pathogens can act synergistically (Whitehorn et al., 2011). The higher prevalence of pathogens in B. terrestris might reflect a different resistance and susceptibility to infections between species. Alternatively, protozoa virulence in B. terrestris may be lower compared to B. pascuorum (Goulson et al., 2012). Moreover, B. terrestris emerges early in the spring, whereas B. pascuorum emerges later (Ornosa & Ortiz-Sánchez, 2004; Ploquin et al., 2013). We conclude that this ecological difference between host species might affect the success of the parasite’s biological cycle and consequently infection prevalence.

AcknowledgmentsTop

We very much appreciate the English editing by Carolyn Daher.


ReferencesTop

Cameron SA, Lozier JD, Strange JP, Koch JB, Cordes N, Solter LF, Griswold TL, 2011. Patterns of widespread decline in North American bumble bees. Proc Natl Acad Sci USA 108: 662-667. https://doi.org/10.1073/pnas.1014743108

Chantawannakul P, Ward L, Boonham N, Brown M, 2006. A scientific note on the detection of honeybee viruses using real-time PCR (TaqMan) in Varroa mites collected from a Thai honeybee (Apis mellifera) apiary. J Invert Pathol 91: 69-73. https://doi.org/10.1016/j.jip.2005.11.001

Fantham HB, Porter A, 1914. The morphology, biology and economic importance of Nosema bombi, N. sp., parasitic in various humble bees (Bombus spp.). Ann Trop Med Parasitol 8: 623-638. https://doi.org/10.1080/00034983.1914.11687667

Francis RM, Kryger P, 2012. Single assay detection of Acute Bee Paralysis Virus, Kashmir Bee Virus and Israeli Acute Paralysis Virus. J Apic Science 56: 137-146. https://doi.org/10.2478/v10289-012-0014-x

Gallot-Lavallée M, Schmid-Hempel R, Vandame R, Vergara CH, Schmid-Hempel P, 2016. Large scale patterns of abundance and distribution of parasites in Mexican bumblebees. J Invert Pathol 133: 73-82. https://doi.org/10.1016/j.jip.2015.12.004

Gamboa V, Ravoet J, Brunain M, Smagghe G, Meeus I, Figueroa J, Ria-o D. de Graaf DC, 2015. Bee pathogens found in Bombus atratus from Colombia: A case study. J Invert Pathol 129: 36-39. https://doi.org/10.1016/j.jip.2015.05.013

Genersch E, Yue C, Fries I, de Miranda JR, 2006. Detection of Deformed wing virus, a honey bee viral pathogen, in bumble bees (Bombus terrestris and Bombus pascuorum) with wing deformities. J Invert Pathol 91: 61-63. https://doi.org/10.1016/j.jip.2005.10.002

Goulson D, 2010. Bumblebees: behaviour, ecology, and conservation. Oxford University Press, Oxford, UK. https://doi.org/10.1017/cbo9780511778230.025

Goulson D, Whitehorn P, Fowley M, 2012. Influence of urbanisation on the prevalence of protozoan parasites of bumblebees. Ecol Entomol 37: 83-89. https://doi.org/10.1111/j.1365-2311.2011.01334.x

Goulson D, Nicholls E, Botías C, Rotheray EL, 2015. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 347(6229): 1255957. https://doi.org/10.1126/science.1255957

Graystock P, Yates K, Darvill B, Goulson D, Hughes WHO, 2013. Emerging dangers: deadly effects of an emergent parasite in a new pollinator host. J Invert Pathol 114: 114-119. https://doi.org/10.1016/j.jip.2013.06.005

Graystock P, Goulson D, Hughes WHO, 2014. The relationship between managed bees and the prevalence of parasites in bumblebees. PeerJ 2: e522. https://doi.org/10.7717/peerj.522

Graystock P, Goulson D, Hughes WHO, 2015. Parasites in bloom: flowers aid dispersal and transmission of pollinator parasites within and between bee species. Proc R Soc B 282: 20151371. https://doi.org/10.1098/rspb.2015.1371

Martín-Hernández R, Meana A, Prieto L, Martínez-Salvador A, Garrido-Bailón E, Higes M 2007. Outcome of colonization of Apis mellifera by Nosema ceranae. Appl Environ Microbiol 73: 6331-6338. https://doi.org/10.1128/AEM.00270-07

Meeus I, de Graaf DC, Jans K, Smagghe G, 2010. Multiplex PCR detection of slowly-evolving trypanosomatids and neogregarines in bumblebees using broad-range primers. J Appl Microbiol 109: 107-115.

Meeus I, Brown MJF, de Graaf DC, Smagghe G 2011. Effects of invasive parasites in bumble bee declines. Conserv Biol 25: 662-671. https://doi.org/10.1111/j.1523-1739.2011.01707.x

Murray TE, Fitzpatrick U, Brown MJF, Paxton RJ, 2008. Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs. Conserv Genet 9: 653-666. https://doi.org/10.1007/s10592-007-9394-z

Ninyerola M, Pons X, Roure JM, 2007. Objective air temperature mapping for the Iberian Peninsula using spatial interpolation and GIS. Int J Climatol 27: 1231-1242. https://doi.org/10.1002/joc.1462

Ornosa C, 1996. Una nota de atención sobre la introducción artificial de subespecies foráneas de abejorros polinizadores en la Península Ibérica (Hymenoptera, Apidae, Bombinae). Bol Soc Esp Entomol 20: 259-260.

Ornosa C, Ortiz-Sánchez FJ, 2004. Hymenoptera: Apoidea I. In: Fauna Ibérica, vol. 23. Ramos MA et al. (eds.). Museo Nacional de Ciencias Naturales, CSIC. Madrid. 556 pp.

Ortiz-Sánchez FJ, 1992. Introducción de Bombus terrestris terrestris (Linneaus, 1758) en el sur de España para la polinización de cultivos de invernadero (Hymenoptera, Apidae). Bol Soc Esp Entomol 16: 247-248.

Ploquin EF, Herrera JM, Obeso JR, 2013. Bumblebee community homogenization after uphill shifts in montane areas of northern Spain. Oecologia 173: 1649-1660. https://doi.org/10.1007/s00442-013-2731-7

Potts SG, Imperatriz-Fonseca V, Ngo HT, Aizen MA, Biesmeijer JC, Breeze TD, Dicks LV, Garibaldi LA, Hill R, Steele J, Vanbergen AJ (2016) Safeguarding pollinators and their values to human well-being. Nature 540: 220–229. https://doi.org/10.1038/nature20588

Rasmont P, Franzén M, Lecocq T, Harpke A, Roberts SPM, Biesmeijer K, Castro L, Cederberg B, Dvorák L, Fitzpatrick Ú, et al., 2015. Climatic risk and distribution atlas of European bumblebees. Biorisk 10 (Spec Issue), 246 pp.

Ravoet J, Maharramov J, Meeus I, De Smet L, Wenseleers T, Smagghe G, De Graaf DC, 2013. Comprehensive bee pathogen screening in Belgium reveals Crithidia mellificae as a new contributory factor to winter mortality. PLoS One 8(8): e72443. https://doi.org/10.1371/journal.pone.0072443

Scriven JJ, Woodall LC, Goulson D, 2013. Nondestructive DNA sampling from bumblebee faeces. Mol Ecol Res 13 (2): 225-229. https://doi.org/10.1111/1755-0998.12036

Singh R, Levitt AL, Rajotte EG, Holmes EC, Ostiguy N, vanEngelsdorp D, Lipkin WI, dePamphilis CW, Toth AL, Cox-Foster DL, 2010. RNA viruses in hymenopteran pollinators: Evidence of inter-taxa virus transmission via pollen and potential impact on non-Apis hymenopteran species. PLoS ONE 5: e14357. https://doi.org/10.1371/journal.pone.0014357

Tamura K, Stecher G, Peterson D, Filipski A, Kumar S, 2013. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30: 2725-2729. https://doi.org/10.1093/molbev/mst197

Velthuis HHW, van Doorn A, 2006. A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie 37: 421-451. https://doi.org/10.1051/apido:2006019

Whitehorn PR, Tinsle, MC, Brown MJF, Darvill B, Goulson D, 2011. Genetic diversity, parasite prevalence and immunity in wild bumblebees. Proc R Soc Lond 278: 1195-1202. https://doi.org/10.1098/rspb.2010.1550

Williams PH, Osborne JL, 2009. Bumblebee vulnerability and conservation world-wide. Apidologie 40: 367-387. https://doi.org/10.1051/apido/2009025