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:pollinatorspathogen dispersionPCRApicystis bombiCrithidia bombiNosema bombiAbbreviations 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)
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
).
Author´s 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.
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.
Introduction
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 methods
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).
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.
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 discussion
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).
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 bumble-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
., 2013b) 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.
Acknowledgments
We very much appreciate the English editing by Carolyn Daher.
CameronSALozierJDStrangeJPKochJBCordesNSolterLFGriswoldTL2011Patterns of widespread decline in North American bumble bees.108662667https://doi.org/10.1073/pnas.1014743108ChantawannakulPWardLBoonhamNBrownM2006A scientific note on the detection of honeybee viruses using real-time PCR (TaqMan) in Varroa mites collected from a Thai honeybee (Apis mellifera) apiary.916973https://doi.org/10.1016/j.jip.2005.11.001FanthamHBPorterA1914The morphology, biology and economic importance of Nosema bombi, N. sp., parasitic in various humble bees (Bombus spp.).8623638https://doi.org/10.1080/00034983.1914.11687667FrancisRMKrygerP2012Single assay detection of Acute Bee Paralysis Virus, Kashmir Bee Virus and Israeli Acute Paralysis Virus.56137146https://doi.org/10.2478/v10289-012-0014-xGallot-LavalléeMSchmid-HempelRVandameRVergaraCHSchmid-HempelP2016Large scale patterns of abundance and distribution of parasites in Mexican bumblebees.1337382https://doi.org/10.1016/j.jip.2015.12.004GamboaVRavoetJBrunainMSmaggheGMeeusIFigueroaJRia-o D. de GraafDC2015Bee pathogens found in Bombus atratus from Colombia: A case study.1293639https://doi.org/10.1016/j.jip.2015.05.013GenerschEYueCFriesIde MirandaJR2006Detection of Deformed wing virus, a honey bee viral pathogen, in bumble bees (Bombus terrestris and Bombus pascuorum) with wing deformities.916163https://doi.org/10.1016/j.jip.2005.10.002GoulsonD2010
Bumblebees: behaviour, ecology, and conservation. Oxford University Press, Oxford, UK.
https://doi.org/10.1017/cbo9780511778230.025GoulsonDWhitehornPFowleyM2012Influence of urbanisation on the prevalence of protozoan parasites of bumblebees.378389https://doi.org/10.1111/j.1365-2311.2011.01334.xGoulsonDNichollsEBotíasCRotherayEL2015Bee declines driven by combined stress from parasites, pesticides, and lack of flowers.34762291255957https://doi.org/10.1126/science.1255957GraystockPYatesKDarvillBGoulsonDHughesWHO2013Emerging dangers: deadly effects of an emergent parasite in a new pollinator host.114114119https://doi.org/10.1016/j.jip.2013.06.005GraystockPGoulsonDHughesWHO2014The relationship between managed bees and the prevalence of parasites in bumblebees.2e522https://doi.org/10.7717/peerj.522GraystockPGoulsonDHughesWHO2015Parasites in bloom: flowers aid dispersal and transmission of pollinator parasites within and between bee species.28220151371https://doi.org/10.1098/rspb.2015.1371Martín-HernándezRMeanaAPrietoLMartínez-SalvadorAGarrido-BailónEHigesM2007Outcome of colonization of Apis mellifera by Nosema ceranae.7363316338https://doi.org/10.1128/AEM.00270-07MeeusIde GraafDCJansKSmaggheG2010Multiplex PCR detection of slowly-evolving trypanosomatids and neogregarines in bumblebees using broad-range primers.109107115MeeusIBrownMJFde GraafDCSmaggheG2011Effects of invasive parasites in bumble bee declines.25662671https://doi.org/10.1111/j.1523-1739.2011.01707.xMurrayTEFitzpatrickUBrownMJFPaxtonRJ2008Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs.9653666https://doi.org/10.1007/s10592-007-9394-zNinyerolaMPonsXRoureJM2007Objective air temperature mapping for the Iberian Peninsula using spatial interpolation and GIS.2712311242https://doi.org/10.1002/joc.1462OrnosaC1996Una 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).20259260OrnosaCOrtiz-SánchezFJ2004Hymenoptera: Apoidea I. In: Fauna Ibérica23556Ramos MA et al. (eds.).Ortiz-SánchezFJ1992. Introducción de Bombus terrestris terrestris(Linneaus1758
)
en el sur de España para la polinización de cultivos de invernadero (Hymenoptera, Apidae).16247248PloquinEFHerreraJMObesoJR2013Bumblebee community homogenization after uphill shifts in montane areas of northern Spain.17316491660https://doi.org/10.1007/s00442-013-2731-7PottsSGImperatriz-FonsecaVNgoHTAizenMABiesmeijerJCBreezeTDDicksLVGaribaldiLAHillRSteeleJVanbergenAJ2016Safeguarding pollinators and their values to human well-being.540220229https://doi.org/10.1038/nature20588RasmontPFranzénMLecocqTHarpkeARobertsSPMBiesmeijerKCastroLCederbergBDvorákLFitzpatrickÚetal.2015Climatic risk and distribution atlas of European bumblebees.10(Spec Issue)246RavoetJMaharramovJMeeusIDe SmetLWenseleersTSmaggheGDe GraafDC2013Comprehensive bee pathogen screening in Belgium reveals Crithidia mellificae as a new contributory factor to winter mortality.88e72443https://doi.org/10.1371/journal.pone.0072443ScrivenJJWoodallLCGoulsonD2013Nondestructive DNA sampling from bumblebee faeces.13
(
2
):
225229https://doi.org/10.1111/1755-0998.12036SinghRLevittALRajotteEGHolmesECOstiguyNvanEngelsdorpDLipkinWIdePamphilisCWTothALCox-FosterDL2010RNA viruses in hymenopteran pollinators: Evidence of inter-taxa virus transmission via pollen and potential impact on non-Apis hymenopteran species.5e14357https://doi.org/10.1371/journal.pone.0014357TamuraKStecherGPetersonDFilipskiAKumarS2013MEGA6: Molecular Evolutionary Genetics Analysis version 6.0.3027252729https://doi.org/10.1093/molbev/mst197VelthuisHHWvan DoornA2006A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination.37421451https://doi.org/10.1051/apido:2006019WhitehornPRTinsleMCBrownMJFDarvillBGoulsonD2011Genetic diversity, parasite prevalence and immunity in wild bumblebees.27811951202https://doi.org/10.1098/rspb.2010.1550WilliamsPHOsborneJL2009Bumblebee vulnerability and conservation world-wide.40367387https://doi.org/10.1051/apido/2009025