Clara Jabal-Uriel1
Verónica N. Albarracín2Joaquín Calatayud3Mariano Higes1
Raquel Martín-Hernández1,4*- 1Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental (CIAPA), Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal (IRIAF), Consejería de Agricultura de la Junta de Comunidades de Castilla-La Mancha, Marchamalo, Spain
- 2Facultad de Agronomía y Zootecnia de la Universidad Nacional de Tucumán, Tucumán, Argentina
- 3Departamento de Biología, Geología, Física y Química inorgánica, Universidad Rey Juan Carlos, Madrid, Spain
- 4Instituto de Recursos Humanos para la Ciencia y la Tecnología (INCRECYT – ESF/EC-FSE), Fundación Parque Científico y Tecnológico de Castilla – La Mancha, Albacete, Spain
The microsporidia Nosema ceranae is an intracellular parasite of honeybees’ midgut, highly prevalent in Apis mellifera colonies for which important epidemiological information is still unknown. Our research aimed at understanding how age and season influence the onset of infection in honeybees and its development in the colony environment. Adult worker honeybees of less than 24h were marked and introduced into 6 different colonies in assays carried out in spring and autumn. Bees of known age were individually analyzed by PCR for Nosema spp. infection and those resulting positive were studied to determine the load by Real Time-qPCR. The age of onset and development of infection in each season was studied on a total of 2401 bees and the probability and the load of infection for both periods was established with two statistical models. First N. ceranae infected honeybees were detected at day 5 post emergence (p.e.; spring) and at day 4 p.e. (autumn) and in-hive prevalence increased from that point onwards, reaching the highest mean infection on day 18 p.e. (spring). The probability of infection increased significantly with age in both periods although the age variable better correlated in spring. The N. ceranae load tended to increase with age in both periods, although the age-load relationship was clearer in spring than in autumn. Therefore, age and season play an important role on the probability and the development of N. ceranae infection in honeybees, bringing important information to understand how it spreads within a colony.
Introduction
The European honeybee, Apis mellifera, currently faces many stressors that affect colony viability such as pesticides, climate change, loss of biodiversity, pests and pathogens (vanEngelsdorp and Meixner, 2010). Among the latter, the microsporidia Nosema ceranae has a detrimental impact on infected colonies (Higes et al., 2008; Emsen et al., 2020). Initially described in the Asian bee Apis cerana (Fries et al., 1996), N. ceranae was first reported in 2006 in the European honeybee A. mellifera in Europe (Higes et al., 2006) and one year later in Asia (Huang et al., 2007). It has since been reported in other parts of the globe (Klee et al., 2007; Chen et al., 2008; Williams et al., 2008; Giersch et al., 2009), becoming a very widespread and prevalent parasite of honeybees.
N. ceranae, which has been proposed to be reclassified as Vairimorpha ceranae (Tokarev et al., 2020) even though we maintain the original name, is an obligate intracellular parasite that is acquired by ingesting spores from contaminated food or water, during cleaning duties or through trophallaxis (reviewed in Martín-Hernández et al., 2018). After intake, the spores reach the ventriculum where they extrude the polar filament to introduce the sporoplasm into the bees’ epithelial cells in which the parasite completes its cell cycle and reproduces, extending infection throughout the organ and producing severe tissue damage (Higes et al., 2007; García-Palencia et al., 2010; Maiolino et al., 2014; Panek et al., 2018).
At the individual level, N. ceranae infection induces immune modulation (Antúnez et al., 2009; Li et al., 2018), it alters polyethism, promoting infected bees to forage prematurely (Dussaubat et al., 2013; Lecocq et al., 2016) and shortening their lifespan (Higes et al., 2007; Goblirsch et al., 2013). At the colony level, infection reduces colony strength, i. e. population size and honey production (Higes et al., 2010; Botías et al., 2013; Bravo et al., 2014; Simeunovic et al., 2014), that finally can lead to colony collapse (Higes et al., 2008; Higes et al., 2010; Villa et al., 2013; Bekele et al., 2015; Adjlane and Haddad, 2016).
This microsporidia parasitizes colonies all year round, although the levels of infection fluctuate across seasons (Higes et al., 2008; Gisder et al., 2010; Copley et al., 2012; Traver et al., 2012; Gisder et al., 2017; Emsen et al., 2020). While all honeybee castes (queen, drones, workers) can be infected (Higes et al., 2009; Traver and Fell, 2011; Martín-Hernández et al., 2012), the prevalence differs among them, being drones and workers the most infected (Traver and Fell, 2011; Martín-Hernández et al., 2012).
In terms of worker bees, N. ceranae parasitize all ages even though the prevalence of infection varies depending on age (Smart and Sheppard, 2012; Jack et al., 2016; Li et al., 2017). In fact, infection is usually more prevalent in foragers than in younger bees that perform in-hive duties (Smart and Sheppard, 2012; Botías et al., 2013; Jack et al., 2016). Thus, the structure of the colony might play an important role in the dynamics of pathogen transmission as the network in which the colony is organized shields younger bees from infection, which has been proposed as a system of organizational immunity (Naug and Smith, 2007; Baracchi and Cini, 2014). Nevertheless, nurse bees are also susceptible to infection (Smart and Sheppard, 2012; Jack et al., 2016) and indeed, N. ceranae appears to be capable of infecting larvae (Eiri et al., 2015; Benvau and Nieh, 2017) even under natural conditions (Urbieta-Magro et al., 2019a). However, there is still little knowledge regarding the age at which honeybees begin to be infected by N. ceranae in natural conditions, as the information available so far is limited to groups of bees of the same age, such as recently emerged bees, nurses or foragers (Smart and Sheppard, 2012; Jack et al., 2016; Li et al., 2017). Therefore, the objectives of this work were to determine at what age honeybees are infected by N. ceranae under field conditions and how age and season influence the development of such infections as this information could help to understand the epidemiology of the infection within the colony.
Material and Methods
Experimental Design
This study was carried out in late spring (May-June 2019) and autumn (October 2019) on experimental colonies located at the Centro de Investigación Apícola y Agroambiental (CIAPA), in Marchamalo, Spain (40° 40´55,77” N; 3° 12´ 32,72” W) and it is shown in the Figure 1. Meteorological conditions during the experiment can been seen in Figure S1.
Figure 1 Experimental design of spring and autumn trials. Capped brood frames were taken from a donor apiary and placed overnight in an incubator and the following day, the newly emerged bees (NEB) were marked and placed in the selected recipient colonies. Marked bees were collected every day and processed in the laboratory. The midgut was removed from each bee and molecular analysis were carried out to detect N. ceranae infection and whenever the samples were positive, RT-qPCR was conducted to quantify the parasitic load.
Spring Experiment
This assay was performed during the end of May and the first fortnight of June.
To avoid genetic homogeneity, eight capped worker brood combs were taken from 8 A. mellifera iberiensis colonies located at an experimental apiary 16 km away from the CIAPA the day before the experiment started. Those colonies had been previously tested and they were negative or had a very low prevalence of N. ceranae infection, and they did not show any clinical signs of other diseases. The combs were immediately moved to the laboratory, removing all the remaining adult bees on them, and placed in an incubator (Memmert® IPP 500) at 34 (± 1) °C overnight.
The next day (day 0), the newly emerged adult bees, all of them of less than 24h old, were allotted in three groups of approximately 300 individuals each and marked with an assigned color of enamel paint on the thorax (Posca PC-5M, Mitsubishi Pencil Co.) before they were introduced into three recipient colonies naturally infected with N. ceranae (300 marked bees per colony). In addition, 15 newly emerged bees were collected to check for the absence of N. ceranae and they were frozen at -20°C until analysis. To make it easier for the new bees to be accepted in the recipient colonies, sugar syrup was sprayed upon them prior to their release. The whole procedure was repeated 3 days later and again, 300 marked bees of each group were released into another three different N. ceranae infected colonies. Therefore, a total of 1800 one-day-old marked worker bees were introduced into six recipient colonies. In this way, we had 2 groups of 3 hives, overlapping each other so we were able to collect daily data from at least 3 colonies for the whole studied period (in order to be able to collect samples during working days). During 21 days, 15 marked bees were collected manually every morning, beginning at the same time, from each of the six colonies starting the day after introduction (day 1) and they were taken to the laboratory for analysis.
The recipient colonies were analyzed to determine the basal level of infection. To do this, prior to the introduction of newly emerged bees on day 0, adult worker bees from those recipient colonies were collected after brushing a comb with no brood. From each colony, 25 bees were analyzed individually (honeybee per honeybee) to determine how many of them were infected by Nosema spp.
Autumn Experiment
This assay was carried out in October following the same experimental design as described above. In this case, 6 capped brood combs were taken from 6 colonies in the same donor apiary to avoid genetic homogeneity, taken to the lab and maintained under the same conditions. The newly emerged bees were marked as indicated above, again in two different batches, this time separated by 2 days. In this case, the absence of infection was determined in 15 newly emerged bees sampled each day of marking and release (day 0), such that 30 bees were analyzed in total. Regarding the recipient colonies, three colonies used for the spring trial remained alive and they were used again in this assay. The other three colonies died (colonies 1, 3 and 6) and they were replaced by three other naturally infected colonies (colonies 7, 8 and 9). During this trial, the marked bees were collected over 22 days as indicated and Nosema spp. infection of all the recipient colonies was again determined as described above (25 worker bees/colony) at the day of introduction.
Molecular Detection of Nosema spp.
Once the marked bees were collected each morning, they were transported to the laboratory and processed immediately. The gut of each honeybee was carefully removed by pulling the last segment of the abdomen with sterile tweezers in a laminar flow cabinet (Telstar AV–30/70). Each gut was cut with sterile dissection scissors, separating the ventriculum and rectum by a cut at half of the ileum. The first part, containing the ventriculum and half of the ileum, was used in this study. Each sample was placed in one well of a 96-well plate (Qiagen®) containing 250 µL of sterile PBS buffer and four 2 mm glass beads (Sigma®) for homogenization. The tissues were homogenized for 2 minutes at 30 Hz (TyssueLyser II, QIAgen®). A well was left with PBS buffer and the reagents alone as a control for DNA extraction between the bees from different colonies and days.
In both the spring and autumn trials, the presence of Nosema spp. was checked in newly emerged bees collected on day 0 and those adult bees collected to determine the basal infection in the recipient colonies, processing the bees individually. As such, the full abdomen of each honeybee was placed in a well of a 96-well plate (Qiagen ®) containing 250 µL of nuclease-free water (Sigma®) and four 2 mm glass beads (Sigma®), and homogenized as indicated above.
DNA was extracted as described previously (Urbieta-Magro et al., 2019b) transferring 50 µL of the homogenate to 50 µL of a Tris-HCl lysis solution (10 nM Tris-HCl [pH 8.0], 1 nM EDTA - TE). All the DNA samples obtained were kept frozen (-20°C) until further analysis. To assess the presence of Nosema spp., conventional triplex PCR was performed on all the samples using gelified plates (BioTools®) (Martín-Hernández et al., 2012) in a Mastercycler® ep gradient S (Eppendorf) and the resulting amplicons were analyzed in a QIAxcel Advanced System (Qiagen®) following the protocol described elsewhere (Urbieta-Magro et al., 2019b), including an internal control of A. mellifera DNA. In addition, the aforementioned DNA extraction controls, non-template controls (NTCs), and a positive control of N. ceranae and N. apis were included in the reaction plates.
Quantification of the N. ceranae Load Per Bee
The bees of known age (marked bees) that were detected positive for Nosema spp. were afterwards analyzed by real time quantitative-PCR (RT-qPCR) to detect the polar tubule protein 3 (PTP3) gene and determine the N. ceranae load per honeybee using the primers and conditions as previously described in Urbieta-Magro et al. (2019b). This was performed in 384-well plates in a Roche LightCycler® 480 thermocycler in a final reaction volume of 20 µL. The amplification cycles were analyzed with the LightCycler® 480 software v1.5.1 (Roche Diagnostics GmbH, Basel, Switzerland) and the crossing point (Cp) was recorded in all the samples using the second maximum derivative method for its calculation based on the standard curve. For each sample, there were two replicates in the same kinetic qPCR run (intra-assay variation), and negative and N. ceranae positive controls were analyzed in parallel. The parasite load (pg/µL) was quantified in all samples relative to the specific synthetic oligonucleotides (gBlocks®, IDT DNA Technologies, Coralville, Iowa, USA) for the N. ceranae-PTP3 gene fragment (Urbieta-Magro et al., 2019b). The standard curve for quantification was prepared following the manufacturer’s protocol at an initial concentration of 10 ng/μL (in TE) and using serial dilutions up to 1 × 10−14 ng/μL. In addition, the GAPDH gene of A. mellifera was analysed in those samples to check for DNA integrity using the primers and probe at the same concentration as described in Martín-Hernández et al. (2017), with 5 µL of the DNA template in a final reaction of 20 µL using a Roche LightCycler® 480 thermocycler. All RT-qPCR programmes consisted of an initial 10 min denaturation step at 95°C, and 45 cycles of 10 seconds at 95°C and 30 seconds of annealing at 60 °C, 1 second at 72°C and ended by a cooling step at 40°C for 30 seconds.
Statistical Analyses
The probability of Nosema spp. infection in bees was analyzed using a Generalized Linear Mixed Model fit by maximum likelihood (Laplace approximation) with binomial distribution and logit link function. Hive and sampling date were included as random variables to account for the potential effects of colony idiosyncrasies and daily climatic variations. As fixed variables we included age and its quadratic term, season and the basal infection of the recipient colonies in the moment of introduction. Moreover, an interaction between age and season was also included. Starting from this full model, we followed a backward model selection based on the second order variant of Akaike Information Criteria (AICc). In case two models were equivalent (i. e. ΔAICc ≤ 2) the model with less parameters was chosen, following a rule of parsimony. We calculated the model’s goodness of fit following Nakagawa and Schielzeth (2013) as implemented in the R package MuMIn (Barton, 2020). A Linear Mix Model fit by maximum likelihood was used for the model of N. ceranae load. The dependent variable was the log transformed load and the fixed variables were age and season. We included the same fixed and random variables as explained above, using also a backward model selection based on AICc. Calculations were made in R (R 4. 1. 0. http://cran.r-project.org/) using the lme4 package (Bates et al., 2015) to fit mixed models and using the AICcmodavg package (Mazzerolle, 2019) to compute the AICc scores. In addition, the loads were analyzed using SPSS software for Windows (Version 12.0), performing an ANOVA and Games-Howell post hoc analysis to correlate the N. ceranae loads among the different ages (in days) of the bees during and within both seasons.
Results
Infection With Nosema spp.
All DNA extraction controls and NTCs did not show any amplification product, indicating the absence of cross-contamination during sample processing and analysis. Significantly, N. apis was not detected in any honeybee throughout the study.
None of the 15 newly emerged bees analyzed were infected by N. ceranae on the day of introduction (day 0) in spring. However, a single honeybee of the 30 bees analyzed in autumn was positive for N. ceranae, which represents a 3.3% of newly emerged bees infected in autumn and a 2.2% among all the newly emerged bees (spring and autumn) analyzed on day 0. The basal infection in the recipient colonies ranged from 8-20% in spring and 4-28% in autumn (Table 1).
Table 1 Daily percentage (and absolute number) of bees infected per day detected positive for N. ceranae infection in each colony during spring and autumn assays.
A total of 2401 bees (1067 in spring and 1334 in autumn) of known age (from day 1 until day 22 post-emergence - p.e.) were collected and analyzed in this study. Of those collected in spring, 35% were positive for N. ceranae when analyzed by conventional triplex-PCR whereas 24% of those collected in autumn were positive (Table 1 and Figure 2).
Figure 2 Mean percentage of N. ceranae in spring and autumn. The percentage in spring (solid line) starts one day later than in autumn (dotted line) and it remains higher except on days 6 and 13 p.e. On the last day of autumn, the percentage of infection drops dramatically.
In spring the first N. ceranae infected bees were detected at day 5 p.e. in the three colonies analyzed (33%, 7%, 7% respectively; Table 1 and Figure S2). In autumn, a single honeybee was positive in one colony on day 4 p.e. (Table 1) and 4 infected bees were detected in two colonies on the following day. From that point onwards, infection was detected in most of the colonies in spring and autumn, and the percentage of bees infected by N. ceranae generally increased (Figure 2, Figure S2). In spring, the highest mean infection was observed on day 18 p.e. and it remained above 75% until day 21 p.e. This was the general trend seen in all the colonies except for colonies 2 and 4, where a decrease was observed at the end of the study period (Table 1). In fact, this latter colony was that with the lowest infection in spring (max. 55%). Two colonies reached 100% infection, colony 2 on days 17 and 18 p.e., and colony 6 on day 15 p.e. However, only one honeybee was found in this latter colony that was positive on day 15, which was the actually last day when any marked honeybee was observed.
In autumn, the mean infection increased until day 21 p.e. when it reached its highest peak (70%), which was followed by a fall to 28% on day 22 p.e. (Figure 2). In general, the mean relative infection was lower in this season than in spring, a trend observed in all the colonies, and no colony reached 100% infection. Moreover, colony 4 was again the colony with the lowest level of infection (max. 40%).
The basal infection by N. ceranae of the recipient colonies did not have a significant effect on the probability of infection of the marked bees (p>0.05), therefore this parameter was not selected in the best model (Table 2). Indeed, the earliest detection of N. ceranae infection was seen in colony 4, associated with 8% basal infection. According to the observations, the best model based on the AICc included the age of the bee, its quadratic term and the interaction of both with the season (Table 2 and Table S1). The overall model fit was R2m = 0.47; R2c = 0.53. On average, the probability of infection was higher in spring (least square mean = 0.364, SE = 0.564) than in autumn (mean = 0.170; SE = 0.564). Moreover, although the probability of infection increased significantly with age in both periods, the effect of age was stronger in spring than in autumn (Figure 3).
Table 2 Results of model selection based on AICc for the probability of infection by N. ceranae and the N. ceranae load in adult bees.
Figure 3 Fitted probability of infection by N. ceranae in spring and autumn by GLMM. There is an exponential growth of infection in spring until the onset of foraging, whereas in autumn this probability remains lower. The shadow area represents 95% confidence interval.
N. ceranae Load in the Infected Bees
This parameter was analyzed in 666 bees which have been determined positive by conventional PCR, as out of the 695 bees that were positive for N. ceranae infection, 6 wells broke and the DNA was poorly preserved in 23 samples (GAPDH not detected), excluding these from the subsequent analysis. All the extraction controls and NTCs gave a negative result, indicating the absence of any contamination or non-specific DNA amplification.
The load of N. ceranae-PTP3 gene was calculated using a standard curve and the limit of quantification (LOQ) was established at 2.5 × 10-6 pg/µL. After analysis, 30 samples were found to be below the LOQ and therefore, they were also removed from the subsequent calculations. As such, the N. ceranae load was calculated for 357 bees collected in spring and 279 bees collected in autumn (636 bees in total). The average concentration of the N. ceranae-PTP3 DNA (pg/µL) was assessed in both seasons (Figure 4) and no significant differences were evident when the load on a given day was compared between them, except between the bees at 21 days p.e., that had higher load in spring (Games-Howell, p < 0.05). This specific day had a significantly higher load than other days in spring but also, to some other of autumn. The N. ceranae-PTP3 load in autumn did not show significant differences between the days of sampling (Games-Howell; p > 0.05), whereas in spring the bees from days 8 and 9 p.e. had significantly lower loads than those bees from days 14, 17, 19, 20 and 21 p.e. (Games-Howell, p < 0.05).
Figure 4 Boxplots showing the N. ceranae-PTP3 load of picograms in logarithmic scale (log pgr) of the infected bees on each day in spring and autumn. The line represents the median while the box represents the 50% of the observations and the whiskers reach the 1.5 x interquartile range. Outliers are also shown as dots.
The best model to explain the N. ceranae load in the bees included the age of the honeybee interacting with season (Table 2 and Table S2). The overall model fit was R2m = 0.12; R2c = 0.16. On average, the load was higher in spring (least square mean = 0.028, SE = 1.48) than in autumn (least square mean = 0.007, SE = 1.52). The N. ceranae load in the bees tended to increase with age in both periods, although the age-load relationship was steeper in spring than in autumn when it was almost flat (Figure 5).
Figure 5 Prediction of N. ceranae-PTP3 load of picograms in logarithmic scale (log pgr) in infected bees in spring and autumn. The load in spring is lower during the initial days and it increases steeply, whereas in autumn the trend remains flatter across the days.
Discussion
Honeybee colonies have evolved as high-level units of biological organization, assembling unified societies of lower-level units. In fact, a honeybee colony is considered a single living entity formed by a society of many thousands of individuals that functions as a unified whole or superorganism (Seeley, 2010). These superorganisms have developed social immunity by coordinating behavioral interactions among nest mates, which results in colony-level immune responses (reviewed in Evans and Spivak, 2010). Among these collective defenses, social organization has been proposed to represent a barrier to pathogen transmission (Evans and Spivak, 2010). Indeed, social organization leads to organizational immunity (Naug and Smith, 2007) whereby the nurse bees and larvae are in the more isolated part of the colony to reduce their contact with pathogens.
However, despite this complex colony organization, our data indicate that adult bees begin to be infected by N. ceranae as young as 4 (autumn) or 5 days (spring) after their emergence. The number of bees infected increases in a time dependent manner, first moderately until day 9 p.e. and then rising remarkably from this point until the age when bees become foragers. To the best of our knowledge, this is the first report confirming the age of infection of worker bees at a colony level, covering every day from their emergence until the onset of foraging and comparing two different seasons. Other reports investigating groups of bees of a range of ages have confirmed older bees to be those with the higher prevalence of N. ceranae infection (Higes et al., 2008; Meana et al., 2010; Smart and Sheppard, 2012; Jack et al., 2016; Li et al., 2017). Our data show that infection at early ages was not anecdotal as a total of 12 bees from 5 colonies were infected by day 5 p.e. Therefore, despite the comb compartmentalization to reduce the overlap of bees of different ages (Naug, 2008; Baracchi and Cini, 2014), our data clearly demonstrate that young bees are exposed to N. ceranae spores. In this study, a single newly emerged bee was infected by N. ceranae and just in autumn and no other infected-bee was detected until day 4 p.e. Therefore, this low prevalence in newly emerged bees did not have a major influence on the results. In fact, of the 695 positive bees only 15 could have been infected from when they emerged (2.2% infected in day 0). This may not be surprising since N. ceranae infected capped brood has been reported in the same location which might acquire the infection during larvae stage, probably due to the food provide by nurse bees (Urbieta-Magro et al., 2019a). Due to the absence of treatments to control Nosema infection and the high prevalence of N. ceranae in Spain, procurement of Nosema-free colonies from which to obtain such a large number of newly emerged bees is currently not possible.
Some differences were found in the development of the infection when compared spring and autumn, which was higher in spring and with an earlier increase in the probability of infection (Figure 3). Similar results have been reported previously when comparing bees of unknown age collected in spring, summer and autumn (Mulholland et al., 2012; Emsen et al., 2020). In our study, we focused on spring rather than in summer due to the harsh conditions (high average temperature and low relative humidity in the experimental site, central Spain) during that season since the time required for the collection of samples could influence the viability of colonies during the experiment.
There was considerable variability in the mean N. ceranae load per bee, consistent with previous results (Mulholland et al., 2012). In the first days of infection in spring, the N. ceranae-PTP3 load was significantly lower than at the end but in autumn there were no differences over the days studied (Figure 4). By contrast, no significant differences were found when comparing bees of the same age between spring and autumn, except on day 21 p.e. The different patterns observed between spring and autumn could be related to colony size and the physiological features of the bees in those seasons. In spring, colonies were larger and all the colonies except one had a shallow super, while on the contrary, the colonies were just the nest in autumn. This might influence the distribution of the honeybee population in the colony and therefore, the possibility of contact among bees of different ages. By contrast, previously reported physiological differences between spring and winter bees might explain the different patterns found here. Winter (or diutinus) bees have swollen hypopharyngeal glands, increased fat body size storing proteins and fat, higher vitellogenin concentration and a considerably extended worker lifespan (Morse, 1985; Ricigliano et al., 2018; Dostálková et al., 2021), and they also display differential microbiota (Kešnerová et al., 2020) and different content of physiologically essential, electrolytic, enzymatic of toxic bioelements (Ptaszyńska et al., 2018). These bees exhibit as well a more intense immune response with stronger expression of antimicrobial genes and more intense antimicrobial activity (Dostálková et al., 2021). This might enhance the bee’s resistance to infection, reducing the number of bees infected compared to spring.
The size of the dataset collected here after analysing more than two thousand bees (from 9 colonies) allowed us to test the effect of age and season on the risk of N. ceranae infection (i) and the parasite load (ii). For both models, age was correlated with season, and the effects of age were stronger in spring than in autumn, therefore, an old honeybee in spring has a higher probability of being infected and it will have a higher parasite load (Figures 3, 5), what is in accordance with other studies where prevalence of N. ceranae infection was found higher in spring than in autumn (Gisder et al., 2017; Emsen et al., 2020). Regarding the probability of infection by N. ceranae (i), our model explains 47% of the cases considering only the variables age and season. Other factors like the colony and the sample date enhanced the model fit (53%) but they were included as random effects as some information that might have helped explain the probability of infection was not available. These factors include the genetic background of the colonies, given that nuclear and mitochondrial variability among A. mellifera iberiensis colonies in the same apiary have been published (Henriques et al., 2021b), to some unpredictable daily events (e.g. the time required to sample each colony increased as the number of bees dropped throughout the study, weather conditions, etc.). Similarly, the model of parasite load (ii) explains 12% of the variability with age and season. The high percentage of unexplained variability could be due to factors like individual honeybee genetics, differential microbiota (Engel et al., 2015; reviewed in Zheng et al., 2018), coinfection with other pathogens (Forsgren and Fries, 2010; Costa et al., 2011; Doublet et al., 2015; Zheng et al., 2015; Ptaszyńska et al., 2021) or the nutraceutical properties of food (Nanetti et al., 2021), which might differ among the colonies.
The establishment of infection in very young bees and the different progression of infection between seasons might have important consequences on the posterior development and viability of the colony. Indeed, N. ceranae infection plays an important role on honeybee polyethism (Dussaubat et al., 2013; Lecocq et al., 2016; Natsopoulou et al., 2016; Mayack et al., 2021) and colony dynamics (Higes et al., 2008; Goblirsch et al., 2013; Wolf et al., 2016; Biganski et al., 2017).
It would be interesting to further investigate these factors that might influence the probability of infection and parasitic load. As a matter of fact, new tools such SNPs (single nucleotide polymorphism) assays are being developed to improve our understanding of immune gene diversity in honeybees (Henriques et al., 2021a). Both the quality and quantity of the feed the bees receive and the presence of nutraceutical substances that may be present in it could modify the development of the infection. Indeed, high-quality or high-diversity pollen have been reported to increase the survival of N. ceranae infected honeybees (Di Pasquale et al., 2013; Ferguson et al., 2018; Castelli et al., 2020) despite that it could increase the parasite load (Jack et al., 2016). Certainly, there are some differences in pollen quality according to season, being generally more abundant and diverse in spring than in autumn. Although we did not collect pollen data in this study, the diversity of plants available for foraging in autumn in the experimental area (Asteraceae and Cruciferae, Dr. González Porto, personal communication) was adequate to cover the needs of the colony and did not require additional supplementation. On the other hand, naturally derived bioactive compounds that are gathered by bees from plants could have beneficial effects on health (nutraceutical effect), as recently shown for N. ceranae infection (Nanetti et al., 2021), and they may vary among colonies as foragers select different plants (González-Porto, personal communication). With regard to coinfection with other pathogens, the interaction of N. ceranae have has been investigated mainly with N. apis (Forsgren and Fries, 2010) and deformed wing virus (Costa et al., 2011; Doublet et al., 2015; Zheng et al., 2015). Hence, exploring the effects of other viruses on N. ceranae infection could be of interest.
Finally, the bee’s microbiota might also influence the age of infection or even its development, although this has not yet been well stablished. Newly emerged bees have very few or no gut bacteria but if they are exposed to natural hive conditions, as in this study, a stable community dominated by core microbial species will be established by days 4 to 6 p. e. (Powell et al., 2014; reviewed in Kwong and Moran, 2016). Nevertheless, some non-core bacterial species (Engel et al., 2015; reviewed Zheng et al., 2018) or even variants of core gut symbionts (Rubanov et al., 2019) might be present at variable levels among individuals or found in the colony but not in every individual and the presence of these bacteria even differ between seasons (Kešnerová et al., 2020), possibly explaining the individual variation observed here. Besides that, N. ceranae infection modulates the gut microbiota by modifying the relative abundance of some bacterial species. Indeed, positive associations have been reported between this microsporidia and some major bacteria such as Gilliamella. apicola, Snodgrasella alvi and Bifidobacteria spp. (Huang et al., 2018; Rubanov et al., 2019; Zhang et al., 2019; Paris et al., 2020). Indeed, it has been suggested that associations between N. ceranae and normal gut microbiota could sustain host survival and benefit the pathogen, as bees with an enhanced natural microbiota have higher parasite loads but higher survival, increasing disease transmission, (Zhang et al., 2021). Therefore, the study of the interaction between the parasite and the microbiota highly deserves future research.
In summary, age and season play an important role on the probability and the dynamics of N. ceranae infection in bees. Although some other factors might influence this, as indicated above, this is a good starting point to understand how N. ceranae spreads within a colony.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
The study was designed by RM-H, VA, and MH. The experiments (marking bees, collection of samples from colonies and subsequent molecular analysis) were done by CJ-U and VA. The statistical analysis was performed by JC and CJ-U. RM-H and CJ-U wrote the manuscript, which was revised and approved by all the other authors (VA, JC, and MH).
Funding
This work was funded by the Consejería de Educación, Cultura y Deportes, of the Junta de Castilla – La Mancha (European Regional development Fund) project No. SBPLY/19/180501/000334. INCRECYT program funded by ESF/EC (Fondo Social Europeo). CJ-U contract funded by the Ministerio de Asuntos Económicos y Transformación Digital (grant no. BES-2017-080176, RTA-2015-0003-CO3-01).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors wish to thank D. Aguado, V. Albendea, J. Almagro, C. Botías, T. Corrales, M. Gajero, J. García and C. Uceta, of the Honeybee Pathology Laboratory for their technical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.823050/full#supplementary-material
References
Adjlane, N., Haddad, N. (2016). Effect of Some Honeybee Diseases on Seasonal Mortality of Apis Mellifera Intermissa in Algeria Apiaries. Proc. Zool. Soc 71, 83–87. doi: 10.1007/s12595-016-0188-5
Antúnez, K., Martín-Hernández, R., Prieto, L., Meana, A., Zunino, P., Higes, M. (2009). Immune Suppression in the Honey Bee (Apis Mellifera) Following Infection by Nosema Ceranae (Microsporidia). Environ. Microbiol. 11, 2284–2290. doi: 10.1111/j.1462-2920.2009.01953.x
Baracchi, D., Cini, A. (2014). A Socio-Spatial Combined Approach Confirms a Highly Compartmentalised Structure in Honeybees. Ethology 120, 1167–1176. doi: 10.1111/eth.12290
Barton, B. (2020) MuMIn: Multi-Model Interference. R Package Version 1.43.17. Available at: https://CRAN.R-project.org/package=MuMIn (Accessed July 26, 2021).
Bates, D., Mächler, M., Bolker, B. M., Walker, S. C. (2015). Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Softw. 67, 1–48. doi: 10.18637/jss.v067.i01
Bekele, A. Z., Mor, S. K., Phelps, N. B. D., Goyal, S. M., Armién, A. G. (2015). A Case Report of Nosema Ceranae Infection in Honey Bees in Minnesota, USA. Vet. Q. 35, 48–50. doi: 10.1080/01652176.2014.981766
Benvau, L. R., Nieh, J. C. (2017). Larval Honey Bees Infected With Nosema Ceranae Have Increased Vitellogenin Titers as Young Adults. Sci. Rep. 7, 1–8. doi: 10.1038/s41598-017-14702-4
Biganski, S., Kurze, C., Müller, M. Y., Moritz, R. F. A. (2017). Social Response of Healthy Honeybees Towards Nosema Ceranae-Infected Workers: Care or Kill? Apidologie 49, 325–334. doi: 10.1007/s13592-017-0557-8
Botías, C., Martín-Hernández, R., Barrios, L., Meana, A., Higes, M. (2013). Nosema Spp. Infection and its Negative Effects on Honey Bees (Apis Mellifera Iberiensis) at the Colony Level. Vet. Res. 44, 1–15. doi: 10.1186/1297-9716-44-25
Bravo, J., Carbonell, V., Valdebenito, J. T., Figueroa, C., Valdovinos, C. E., Martín-Hernández, R., et al. (2014). Identification of Nosema Ceranae in the Valparaíso District, Chile. Arch. Med. Vet. 46, 487–491. doi: 10.4067/S0301-732X2014000300021
Castelli, L., Branchiccela, B., Garrido, M., Invernizzi, C., Porrini, M., Romero, H., et al. (2020). Impact of Nutritional Stress on Honeybee Gut Microbiota, Immunity, and Nosema Ceranae Infection. Microb. Ecol. 80, 908–919. doi: 10.1007/s00248-020-01538-1
Chen, Y., Evans, J. D., Smith, I. B., Pettis, J. S. (2008). Nosema Ceranae is a Long-Present and Wide-Spread Microsporidian Infection of the European Honey Bee (Apis Mellifera) in the United States. J. Invertebr. Pathol. 97, 186–188. doi: 10.1016/j.jip.2007.07.010
Copley, T. R., Chen, H., Giovenazzo, P., Houle, E., Jabaji, S. H. (2012). Prevalence and Seasonality of Nosema Species in Québec Honey Bees. Can. Entomol. 144, 577–588. doi: 10.4039/tce.2012.46
Costa, C., Tanner, G., Lodesani, M., Maistrello, L., Neumann, P. (2011). Negative Correlation Between Nosema Ceranae Spore Loads and Deformed Wing Virus Infection Levels in Adult Honey Bee Workers. J. Invertebr. Pathol. 108, 224–225. doi: 10.1016/j.jip.2011.08.012
Di Pasquale, G., Salignon, M., Le Conte, Y., Belzunces, L. P., Decourtye, A. (2013). Influence of Pollen Nutrition on Honey Bee Health: Do Pollen Quality and Diversity Matter? PloS One. 8, 8, 1–13. doi: 10.1371/journal.pone.0072016
Dostálková, S., Dobeš, P., Kunc, M., Hurychová, J., Škrabišová, M., Petřivalský, M., et al. (2021). Winter Honeybee (Apis Mellifera) Populations Show Greater Potential to Induce Immune Responses Than Summer Populations After Immune Stimuli. J. Exp. Biol. 224, 1–12. doi: 10.1242/jeb.232595
Doublet, V., Natsopoulou, M. E., Zschiesche, L., Paxton, R. J. (2015). Within-Host Competition Among the Honey Bees Pathogens Nosema Ceranae and Deformed Wing Virus is Asymmetric and to the Disadvantage of the Virus. J. Invertebr. Pathol. 124, 31–34. doi: 10.1016/j.jip.2014.10.007
Dussaubat, C., Maisonnasse, A., Crauser, D., Beslay, D., Costagliola, G., Soubeyrand, S., et al. (2013). Flight Behavior and Pheromone Changes Associated to Nosema Ceranae Infection of Honey Bee Workers (Apis Mellifera) in Field Conditions. J. Invertebr. Pathol. 113, 42–51. doi: 10.1016/j.jip.2013.01.002
Eiri, D. M., Suwannapong, G., Endler, M., Nieh, J. C. (2015). Nosema Ceranae can Infect Honey Bee Larvae and Reduces Subsequent Adult Longevity. PloS One 10, 1–17. doi: 10.1371/journal.pone.0126330
Emsen, B., de la Mora, A., Lacey, B., Eccles, L., Kelly, P. G., Medina-Flores, C. A., et al. (2020). Seasonality of Nosema Ceranae Infections and Their Relationship With Honey Bee Populations, Food Stores, and Survivorship in a North American Region. Vet. Sci. 7, 1–14. doi: 10.3390/VETSCI7030131
Engel, P., Bartlett, K. D., Moran, N. A. (2015). The Bacterium Frischella Perrara Causes Scab Formation in the Gut of its Honeybee Host. MBio 6, 1–8. doi: 10.1128/mBio.00193-15
Evans, J. D., Spivak, M. (2010). Socialized Medicine: Individual and Communal Disease Barriers in Honey Bees. J. Invertebr. Pathol. 103, S62–S72. doi: 10.1016/j.jip.2009.06.019
Ferguson, J. A., Northfield, T. D., Lach, L. (2018). Honey Bee (Apis Mellifera) Pollen Foraging Reflects Benefits Dependent on Individual Infection Status. Microb. Ecol. 76, 482–491. doi: 10.1007/s00248-018-1147-7
Forsgren, E., Fries, I. (2010). Comparative Virulence of Nosema Ceranae and Nosema Apis in Individual European Honey Bees. Vet. Parasitol. 170, 212–217. doi: 10.1016/j.vetpar.2010.02.010
Fries, I., Feng, F., Da Silva, A., Slemenda, S. B., Pieniazek, N. J. (1996). Nosema Ceranae N. Sp. (Microspora, Nosematidae), Morphological and Molecular Characterization of a Microsporidian Parasite of the Asian Honey Bee Apis Cerana (Hymenoptera, Apidae). Eur. J. Protistol. 32, 356–365. doi: 10.1016/S0932-4739(96)80059-9
García-Palencia, P., Martín-Hernández, R., González-Porto, A. V., Marin, P., Meana, A., Higes, M. (2010). Natural Infection by Nosema Ceranae Causes Similar Lesions as in Experimentally Infected Caged-Worker Honey Bees (Apis Mellifera). J. Apic. Res. 49, 278–283. doi: 10.3896/IBRA.1.49.1.08
Giersch, T., Berg, T., Galea, F., Hornitzky, M. (2009). Nosema Ceranae Infects Honey Bees (Apis Mellifera) and Contaminates Honey in Australia. Apidologie 40, 117–123. doi: 10.1051/apido/2008065
Gisder, S., Hedtke, K., Möckel, N., Frielitz, M. C., Linde, A., Genersch, E. (2010). Five-Year Cohort Study of Nosema Spp. In Germany: Does Climate Shape Virulence and Assertiveness of Nosema Ceranae? Appl. Environ. Microbiol. 76, 3032–3038. doi: 10.1128/AEM.03097-09
Gisder, S., Schüler, V., Horchler, L. L., Groth, D., Genersch, E. (2017). Long-Term Temporal Trends of Nosema Spp. Infection Prevalence in Northeast Germany: Continuous Spread of Nosema Ceranae, an Emerging Pathogen of Honey Bees (Apis Mellifera), But No General Replacement of Nosema Apis. Front. Cell. Infect. Microbiol. 7, 301. doi: 10.3389/fcimb.2017.00301
Goblirsch, M., Huang, Z. Y., Spivak, M. (2013). Physiological and Behavioral Changes in Honey Bees (Apis Mellifera) Induced by Nosema Ceranae Infection. PloS One 8, 1–8. doi: 10.1371/journal.pone.0058165
Henriques, D., Lopes, A. R., Chejanovsky, N., Dalmon, A., Higes, M., Jabal-Uriel, C., et al. (2021a). A SNP Assay for Assessing Diversity in Immune Genes in the Honey Bee (Apis Mellifera L.). Sci. Rep. 11, 1–12. doi: 10.1038/s41598-021-94833-x
Henriques, D., Lopes, A. R., Chejanovsky, N., Dalmon, A., Higes, M., Jabal-Uriel, C., et al. (2021b). Mitochondrial and Nuclear Diversity of Colonies of Varying Origins: Contrasting Patterns Inferred From the Intergenic Trnaleu-Cox2 Region and Immune SNPs. J. Apic. Res. 1–4. doi: 10.1080/00218839.2021.2010940
Higes, M., García-Palencia, P., Martín-Hernández, R., Meana, A. (2007). Experimental Infection of Apis Mellifera Honeybees With Nosema Ceranae (Microsporidia). J. Invertebr. Pathol. 94, 211–217. doi: 10.1016/j.jip.2006.11.001
Higes, M., Martín-Hernández, R., Botías, C., Bailón, E. G., González-Porto, A. V., Barrios, L., et al. (2008). How Natural Infection by Nosema Ceranae Causes Honeybee Colony Collapse. Environ. Microbiol. 10, 2659–2669. doi: 10.1111/j.1462-2920.2008.01687.x
Higes, M., Martín-Hernández, R., García-Palencia, P., Marín, P., Meana, A. (2009). Horizontal Transmission of Nosema Ceranae (Microsporidia) From Worker Honeybees to Queens (Apis Mellifera). Environ. Microbiol. Rep. 1, 495–498. doi: 10.1111/j.1758-2229.2009.00052.x
Higes, M., Martín-Hernández, R., Martínez-Salvador, A., Garrido-Bailón, E., González-Porto, A. V., Meana, A., et al. (2010). A Preliminary Study of the Epidemiological Factors Related to Honey Bee Colony Loss in Spain. Environ. Microbiol. Rep. 2, 243–250. doi: 10.1111/j.1758-2229.2009.00099.x
Higes, M., Martín, R., Meana, A. (2006). Nosema Ceranae, a New Microsporidian Parasite in Honeybees in Europe. J. Invertebr. Pathol. 92, 93–95. doi: 10.1016/j.jip.2006.02.005
Huang, W. F., Jiang, J. H., Chen, Y. W., Wang, C. H. (2007). A Nosema Ceranae Isolate From the Honeybee Apis Mellifera. Apidologie 38, 30–37. doi: 10.1051/apido:2006054
Huang, S. K., Ye, K. T., Huang, W. F., Ying, B. H., Su, X., Lin, L. H., et al. (2018). Influence of Feeding Type and Nosema Ceranae Infection on the Gut Microbiota of Apis Cerana Workers. mSystems 3, 1–11. doi: 10.1128/msystems.00177-18
Jack, C. J., Lucas, H. M., Webster, T. C., Sagili, R. R. (2016). Colony Level Prevalence and Intensity of Nosema Ceranae in Honey Bees (Apis Mellifera L.). PloS One 11, 1–20. doi: 10.1371/journal.pone.0163522
Kešnerová, L., Emery, O., Troilo, M., Liberti, J., Erkosar, B., Engel, P. (2020). Gut Microbiota Structure Differs Between Honeybees in Winter and Summer. ISME. J. 14, 801–814. doi: 10.1038/s41396-019-0568-8
Klee, J., Besana, A. M., Genersch, E., Gisder, S., Nanetti, A., Tam, D. Q., et al. (2007). Widespread Dispersal of the Microsporidian Nosema Ceranae, an Emergent Pathogen of the Western Honey Bee, Apis Mellifera. J. Invertebr. Pathol. 96, 1–10. doi: 10.1016/j.jip.2007.02.014
Kwong, W. K., Moran, N. A. (2016). Gut Microbial Communities of Social Bees. Nat. Rev. Microbiol. 14, 374–384. doi: 10.1038/nrmicro.2016.43
Lecocq, A., Jensen, A. B., Kryger, P., Nieh, J. C. (2016). Parasite Infection Accelerates Age Polyethism in Young Honey Bees. Sci. Rep. 6, 1–11. doi: 10.1038/srep22042
Li, W., Chen, Y., Cook, S. C. (2018). Chronic Nosema Ceranae Infection Inflicts Comprehensive and Persistent Immunosuppression and Accelerated Lipid Loss in Host Apis Mellifera Honey Bees. Int. J. Parasitol. 48, 433–444. doi: 10.1016/j.ijpara.2017.11.004
Li, W., Evans, J. D., Li, J., Su, S., Hamilton, M., Chen, Y. (2017). Spore Load and Immune Response of Honey Bees Naturally Infected by Nosema Ceranae. Parasitol. Res. 116, 3265–3274. doi: 10.1007/s00436-017-5630-8
Maiolino, P., Iafigliola, L., Rinaldi, L., De Leva, G., Restucci, B., Martano, M. (2014). Histopathological Findings of the Midgut in European Honey Bee (Apis Mellifera L.) Naturally Infected by Nosema Spp. Vet. Med. Anim. Sci. 2:4. doi: 10.7243/2054-3425-2-4
Martín-Hernández, R., Bartolomé, C., Chejanovsky, N., Le Conte, Y., Dalmon, A., Dussaubat, C., et al. (2018). Nosema Ceranae in Apis Mellifera: A 12 Years Postdetection Perspective. Environ. Microbiol. 20, 1302–1329. doi: 10.1111/1462-2920.14103
Martín-Hernández, R., Botías, C., Bailón, E. G., Martínez-Salvador, A., Prieto, L., Meana, A., et al. (2012). Microsporidia Infecting Apis Mellifera: Coexistence or Competition. Is Nosema Ceranae Replacing Nosema Apis? Environ. Microbiol. 14, 2127–2138. doi: 10.1111/j.1462-2920.2011.02645.x
Martín-Hernández, R., Higes, M., Sagastume, S., Juarranz, Á., Dias-Almeida, J., Budge, G. E., et al. (2017). Microsporidia Infection Impacts the Host Cell’s Cycle and Reduces Host Cell Apoptosis. PloS One 12, 1–17. doi: 10.1371/journal.pone.0170183
Mayack, C., Broadrup, R. L., Schick, S. J., Eppley, E. J., Khan, Z., MacHerone, A. (2021). Increased Alarm Pheromone Component is Associated With Nosema Ceranae Infected Honeybee Colonies. R. Soc Open Sci. 8, 1–6. doi: 10.1098/rsos.210194
Mazzerolle, M. J. (2019) AICcmodavg: Model Selection and Multimodel Inference Based on (Q)AIC©. R Package Version 2.2-2. Available at: https://cran.r-project.org/package=AICcmodavg (Accessed July 26, 2021).
Meana, A., Martín-Hernández, R., Higes, M. (2010). The Reliability of Spore Counts to Diagnose Nosema Ceranae Infections in Honey Bees. J. Apic. Res. 49, 212–214. doi: 10.3896/IBRA.1.49.2.12
Morse, R. (1985). The Illustrated Encyclopedia of Beekeeping. Eds. Morse, R., Hoper, T. (Poole, Dorset: Blandford Press).
Mulholland, G. E., Traver, B. E., Johnson, N. G., Fell, R. D. F. (2012). Individual Variability of Nosema Ceranae Infections in Apis Mellifera Colonies. Insects 3, 1143–1155. doi: 10.3390/insects3041143
Nakagawa, S., Schielzeth, H. (2013). A General and Simple Method for Obtaining R2 From Generalized Linear Mixed-Effects Models. Methods Ecol. Evol. 4, 133–142. doi: 10.1111/j.2041-210x.2012.00261.x
Nanetti, A., Ugolini, L., Cilia, G., Pagnotta, E., Malaguti, L., Cardaio, I., et al. (2021). Seed Meals From Brassica Nigra and Eruca Sativa Control Artificial Nosema Ceranae Infections in Apis Mellifera. Microorganisms 9, 949. doi: 10.3390/microorganisms9050949
Natsopoulou, M. E., McMahon, D. P., Paxton, R. J. (2016). Parasites Modulate Within-Colony Activity and Accelerate the Temporal Polyethism Schedule of a Social Insect, the Honey Bee. Behav. Ecol. Sociobiol. 70, 1019–1031. doi: 10.1007/s00265-015-2019-5
Naug, D. (2008). Structure of the Social Network and its Influence on Transmission Dynamics in a Honeybee Colony. Behav. Ecol. Sociobiol. 62, 1719–1725. doi: 10.1007/s00265-008-0600-x
Naug, D., Smith, B. (2007). Experimentally Induced Change in Infectious Period Affects Transmission Dynamics in a Social Group. Proc. R. Soc B. Biol. Sci. 274, 61–65. doi: 10.1098/rspb.2006.3695
Panek, J., Paris, L., Roriz, D., Mone, A., Dubuffet, A., Delbac, F., et al. (2018). Impact of the Microsporidian Nosema Ceranae on the Gut Epithelium Renewal of the Honeybee, Apis Mellifera. J. Invertebr. Pathol. 159, 121–128. doi: 10.1016/j.jip.2018.09.007
Paris, L., Peghaire, E., Moné, A., Diogon, M., Debroas, D., Delbac, F., et al. (2020). Honeybee Gut Microbiota Dysbiosis in Pesticide/Parasite Co-Exposures is Mainly Induced by Nosema Ceranae. J. Invertebr. Pathol. 172, 107348. doi: 10.1016/j.jip.2020.107348
Powell, J. E., Martinson, V. G., Urban-Mead, K., Moran, N. A. (2014). Routes of Acquisition of the Gut Microbiota of the Honey Bee Apis Mellifera. Appl. Environ. Microbiol. 80, 7378–7387. doi: 10.1128/AEM.01861-14
Ptaszyńska, A. A., Gancarz, M., Hurd, P. J., Borsuk, G., Wiącek, D., Nawrocka, A., et al. (2018). Changes in the Bioelement Content of Summer and Winter Western Honeybees (Apis Mellifera) Induced by Nosema Ceranae Infection. PloS One 13, 1–18. doi: 10.1371/journal.pone.0200410
Ptaszyńska, A. A., Latoch, P., Hurd, P. J., Polaszek, A., Michalska-Madej, J., Grochowalski, >Ł., et al. (2021). Amplicon Sequencing of Variable 16s rRNA From Bacteria and ITS2 Regions From Fungi and Plants, Reveals Honeybee Susceptibility to Diseases Results From Their Forage Availability Under Anthropogenic Landscapes. Pathogens 10, 1–17. doi: 10.3390/pathogens10030381
Ricigliano, V. A., Mott, B. M., Floyd, A. S., Copeland, D. C., Carroll, M. J., Anderson, K. E. (2018). Honey Bees Overwintering in a Southern Climate: Longitudinal Effects of Nutrition and Queen Age on Colony-Level Molecular Physiology and Performance. Sci. Rep. 8, 1–11. doi: 10.1038/s41598-018-28732-z
Rubanov, A., Russell, K. A., Rothman, J. A., Nieh, J. C., McFrederick, Q. S. (2019). Intensity of Nosema Ceranae Infection is Associated With Specific Honey Bee Gut Bacteria and Weakly Associated With Gut Microbiome Structure. Sci. Rep. 9, 1–8. doi: 10.1038/s41598-019-40347-6
Seeley, T. D. (2010). Honeybee Democracy (Princeton, New Jersey: Princeton University Press). doi: 10.1515/9781400835959
Simeunovic, P., Stevanovic, J., Cirkovic, D., Radojicic, S., Lakic, N., Stanisic, L., et al. (2014). Nosema Ceranae and Queen Age Influence the Reproduction and Productivity of the Honey Bee Colony. J. Apic. Res. 53, 545–554. doi: 10.3896/IBRA.1.53.5.09
Smart, M. D., Sheppard, W. S. (2012). Nosema Ceranae in Age Cohorts of the Western Honey Bee (Apis Mellifera). J. Invertebr. Pathol. 109, 148–151. doi: 10.1016/j.jip.2011.09.009
Tokarev, Y. S., Huang, W. F., Solter, L. F., Malysh, J. M., Becnel, J. J., Vossbrinck, C. R. (2020). A Formal Redefinition of the Genera Nosema and Vairimorpha (Microsporidia: Nosematidae) and Reassignment of Species Based on Molecular Phylogenetics. J. Invertebr. Pathol. 169, 107279. doi: 10.1016/j.jip.2019.107279
Traver, B. E., Fell, R. D. (2011). Nosema Ceranae in Drone Honey Bees (Apis Mellifera). J. Invertebr. Pathol. 107, 234–236. doi: 10.1016/j.jip.2011.05.016
Traver, B. E., Williams, M. R., Fell, R. D. (2012). Comparison of Within Hive Sampling and Seasonal Activity of Nosema Ceranae in Honey Bee Colonies. J. Invertebr. Pathol. 109, 187–193. doi: 10.1016/j.jip.2011.11.001
Urbieta-Magro, A., Higes, M., Meana, A., Barrios, L., Martín-Hernández, R. (2019a). Age and Method of Inoculation Influence the Infection of Worker Honey Bees (Apis Mellifera) by Nosema Ceranae. Insects 10, 1–17. doi: 10.3390/insects10120417
Urbieta-Magro, A., Higes, M., Meana, A., Gómez-Moracho, T., Rodríguez-García, C., Barrios, L., et al. (2019b). The Levels of Natural Nosema Spp. Infection in Apis Mellifera Iberiensis Brood Stages. Int. J. Parasitol. 49, 657–667. doi: 10.1016/j.ijpara.2019.04.002
vanEngelsdorp, D., Meixner, M. D. (2010). A Historical Review of Managed Honey Bee Populations in Europe and the United States and the Factors That may Affect Them. J. Invertebr. Pathol. 103, S80–S95. doi: 10.1016/j.jip.2009.06.011
Villa, J. D., Bourgeois, A. L., Danka, R. G. (2013). Negative Evidence for Effects of Genetic Origin of Bees on Nosema Ceranae, Positive Evidence for Effects of Nosema Ceranae on Bees. Apidologie 44, 511–518. doi: 10.1007/s13592-013-0201-1
Williams, G. R., Shafer, A. B. A., Rogers, R. E. L., Shutler, D., Stewart, D. T. (2008). First Detection of Nosema Ceranae, a Microsporidian Parasite of European Honey Bees (Apis Mellifera), in Canada and Central USA. J. Invertebr. Pathol. 97, 189–192. doi: 10.1016/j.jip.2007.08.005
Wolf, S., Nicholls, E., Reynolds, A. M., Wells, P., Lim, K. S., Paxton, R. J., et al. (2016). Optimal Search Patterns in Honeybee Orientation Flights are Robust Against Emerging Infectious Diseases. Sci. Rep. 6, 1–10. doi: 10.1038/srep32612
Zhang, Y., Lu, X., Huang, S., Zhang, L., Su, S., Huang, W. F. (2019). Nosema Ceranae Infection Enhances Bifidobacterium Spp. Abundances in the Honey Bee Hindgut. Apidologie 50, 353–362. doi: 10.1007/s13592-019-00644-5
Zhang, Y., Su, M., Wang, L., Huang, S., Su, S., Huang, W. F. (2021). Vairimorpha (Nosema) Ceranae Infection Alters Honey Bee Microbiota Composition and Sustains the Survival of Adult Honey Bees. Biol. (Basel). 10, 1–14. doi: 10.3390/biology10090905
Zheng, H. Q., Gong, H. R., Huang, S. K., Sohr, A., Hu, F. L., Chen, Y. P. (2015). Evidence of the Synergistic Interaction of Honey Bee Pathogens Nosema Ceranae and Deformed Wing Virus. Vet. Microbiol. 177, 1–6. doi: 10.1016/j.vetmic.2015.02.003
Keywords: Microsporidia, Nosema ceranae, Apis mellifera, honeybee, natural infection, age of infection
Citation: Jabal-Uriel C, Albarracín VN, Calatayud J, Higes M and Martín-Hernández R (2022) Age and Season Effect the Timing of Adult Worker Honeybee Infection by Nosema ceranae. Front. Cell. Infect. Microbiol. 11:823050. doi: 10.3389/fcimb.2021.823050
Received: 26 November 2021; Accepted: 29 December 2021;
Published: 28 January 2022.
Edited by:
Giovanni Cilia, Council for Agricultural and Economics Research (CREA), ItalyReviewed by:
Christopher Mayack, Sabancı University, TurkeyDaniele Alberoni, University of Bologna, Italy
Copyright © 2022 Jabal-Uriel, Albarracín, Calatayud, Higes and Martín-Hernández. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Raquel Martín-Hernández, rmhernandez@jccm.es