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

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. Correspondence should be addressed to Pilar De la Rúa: pdelarua@um.es 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). 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). 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. 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. Copyright © 2017 INIA. This is an open access article distributed under the terms of the Creative Commons Attribution (CC-by) Spain 3.0 License. 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. 15(1), e05SC01, 6 pages (2017) eISSN: 2171-9292 https://doi.org/10.5424/sjar/2017151-9998


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).

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).
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).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 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.

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 bum blebee 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 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.24.1 14.9-33.31.2 0.0-3.6 0.0 0.0-0.00.0 0.0-0.0 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.