Abstract
The health and sustainability of honeybee populations are essential for global food production and ecological balance. More than one-third of agricultural crops depend on pollination, making honeybees indispensable to the global economy. However, colonies worldwide are increasingly threatened by the gut parasite Nosema ceranae, whose recurrent outbreaks cause severe productivity losses and economic damage to the apicultural and agricultural sectors. Despite its prevalence, the nonlinear mechanisms responsible for the persistence and resurgence of Nosema ceranae remain poorly understood. This study aims to develop a mathematical model that interpret the complex epidemiological behaviour of Nosema ceranae under realistic biological conditions, specifically focusing on renewal of colony-level susceptibility and the strained resources available for control and management. A nonlinear Susceptible-Infected-Recovered-Susceptible (SIRS) model incorporating renewal of colony-level susceptibility and resource saturation is formulated and analysed using stability theory and numerical bifurcation techniques. Analytical derivations identify the conditions for transcritical and Hopf bifurcations, while numerical continuation analysis and time series simulations validate the theoretical predictions and reveal the rich dynamical structure. The analysis uncovers two critical thresholds governing colony dynamics. A forward (transcritical) bifurcation marks the transition from disease eradication to endemic persistence, while a Hopf bifurcation arising from resource limitations induces sustained oscillations representing recurrent infection waves. The coexistence of stable equilibria and periodic orbits highlights bistability, indicating that small perturbations or change in control interventions can trigger large amplitude outbreaks. These nonlinear feedbacks provide a mechanistic explanation for the cyclical prevalence of Nosema ceranae observed in distinct colonies. The proposed framework establishes, for the first time, how the interplay between partial recovery and limited resource availability can drive complex epidemic patterns in honeybee colonies. Beyond its theoretical contribution, the study provides actionable insight for the agriculture industry: efficient allocation of control resources and timely interventions are essential to prevent recurrent epidemics that threaten pollination services and agricultural productivity. The results bridge mathematical modelling and ecological management, offering a predictive foundation for mitigating the economic and environmental impact of Nosema ceranae on global food security.
Data availability
Simulation data and codes will be shared upon reasonable request
References
Hung, K.-L.J., Kingston, J. M., Albrecht, M., Holway, D. A. & Kohn, J. R. The worldwide importance of honey bees as pollinators in natural habitats. Proc. R. Soc. B: Biol. Sci. 285(1870), 20172140 (2018).
Khalifa, S. A. M. et al. Overview of bee pollination and its economic value for crop production. Insects 12 (8), 688 (2021).
Babin, A. et al. Distribution of infectious and parasitic agents among three sentinel bee species across european agricultural landscapes. Sci. Rep. 14(1), 3524 (2024).
Goulson, D., Nicholls, E., Botías, C. & Rotheray, E. L. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 347(6229), 1255957 (2015).
Aizen, M. A. & Harder, L. D. The global stock of domesticated honey bees is growing slower than agricultural demand for pollination. Curr. Biol. 19(11), 915–918 (2009).
Genersch, Elke. Honey bee pathology: current threats to honey bees and beekeeping. Appl. Microbiol. Biotechnol. 87(1), 87–97 (2010).
Collison, E., Hird, H., Cresswell, J. & Tyler, C. Interactive effects of pesticide exposure and pathogen infection on bee health-a critical analysis. Biol. Rev. 91(4), 1006–1019 (2016).
El-Sayed, A. S. A. et al. Biological control of nosemosis in apis mellifera l. with acacia nilotica extract. Sci. Rep. 14(1), 28340 (2024).
Fries, Ingemar. Nosema ceranae in european honey bees (apis mellifera). J. Invertebrate Pathol. 103, S73–S79 (2010).
Martín-Hernández, R. et al. Nosema ceranae in apis mellifera: a 12 years postdetection perspective. Environ. Microbiol. 20(4), 1302–1329 (2018).
Park, D. et al. Uncovering the novel characteristics of asian honey bee, apis cerana, by whole genome sequencing. BMC Genom. 16(1), 1 (2015).
Chen, Y. P. & Huang, Z. Y. Nosema ceranae, a newly identified pathogen of apis mellifera in the usa and asia. Apidologie 41(3), 364–374 (2010).
Koeniger, N., Koeniger, G. & Delfinado-Baker, M. Observations on mites of the asian honeybee species (apis cerana, apis dorsata, apis florea). Apidologie 14(3), 197–204 (1983).
Kairo, G. et al. Nosema ceranae, fipronil and their combination compromise honey bee reproduction via changes in male physiology. Sci. Rep. 7(1), 8556 (2017).
Holt, H. L. & Grozinger, C. M. Approaches and challenges to managing nosema (microspora: Nosematidae) parasites in honey bee (hymenoptera: Apidae) colonies. J. Econ. Entomol. 109(4), 1487–1503 (2016).
Fisher, R. A., Gollan, B. & Helaine, S. Persistent bacterial infections and persister cells. Nat. Rev. Microbiol. 15(8), 453–464 (2017).
Barrett, M. P., Kyle, D. E., Sibley, L. D., Radke, J. B. & Tarleton, R. L. Protozoan persister-like cells and drug treatment failure. Nat. Rev. Microbiol. 17(10), 607–620 (2019).
Hopkins, S. R. et al. Environmental persistence of the world’s most burdensome infectious and parasitic diseases. Front. Public Health 10, 892366 (2022).
Gisder, S., Schüler, V., Horchler, L. L., Groth, D. & Genersch, E. 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 (2017).
Jabal-Uriel, C., Albarracín, V. N., Calatayud, J., Higes, M. & Martín-Hernández, R. Age and season effect the timing of adult worker honeybee infection by nosema ceranae. Front. Cell. Infect. Microbiol. 11, 823050 (2022).
Mendoza, Y. et al. Nosema ceranae winter control: study of the effectiveness of different fumagillin treatments and consequences on the strength of honey bee (hymenoptera: Apidae) colonies. J. Econ. Entomol. 110(1), 1–5 (2017).
Burnham, A. J. Scientific advances in controlling nosema ceranae (microsporidia) infections in honey bees (apis mellifera). Front. Vet. Sci. 6, 79 (2019).
Betti, M. I., Wahl, L. M. & Zamir, M. Effects of infection on honey bee population dynamics: a model. PloS one 9(10), e110237 (2014).
Laomettachit, T., Liangruksa, M., Termsaithong, T., Tangthanawatsakul, A. & Duangphakdee, O. A model of infection in honeybee colonies with social immunity. PloS one 16(2), e0247294 (2021).
Jabal-Uriel, C. et al. Epidemiology of the microsporidium nosema ceranae in four mediterranean countries. Insects 13(9), 844 (2022).
Schüler, V. et al. Significant, but not biologically relevant: Nosema ceranae infections and winter losses of honey bee colonies. Commun. Biol. 6(1), 229 (2023).
Kermack, W. O. & McKendrick, A. G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A, Contain. Pap. Math. Phys. Character 115(772), 700–721 (1927).
Mohd, M. H. & Sulayman, F. Unravelling the myths of r0 in controlling the dynamics of covid-19 outbreak: a modelling perspective. Chaos, Solitons Fractals 138, 109943 (2020).
Salman, A. M., Ahmed, I., Mohd, M. H., Jamiluddin, M. S. & Dheyab, M. A. Scenario analysis of covid-19 transmission dynamics in malaysia with the possibility of reinfection and limited medical resources scenarios. Comput. Biol. Med. 133, 104372 (2021).
Salman, A. M. & Mohd, M. H. Dynamical systems analysis of a reaction-diffusion sirs model with optimal control for the covid-19 spread. Comput. Methods Biomech. Biomed. Eng. 2024, 1–18 (2024).
Salman, A. M. et al. An optimal control of sirs model with limited medical resources and reinfection problems. Malays. J. Fundam. Appl. Sci. 18(3), 332–342 (2022).
MacInnis, C. I., Luong, L. T. & Pernal, S. F. A tale of two parasites: responses of honey bees infected with nosema ceranae and lotmaria passim. Sci. Rep. 13(1), 22515 (2023).
Cronin, A. L., Molet, M., Doums, C., Monnin, T. & Peeters, C. Recurrent evolution of dependent colony foundation across eusocial insects. Annu. Rev. Entomol. 58, 37–55 (2013).
Friedman, D. A., Johnson, B. R. & Linksvayer, T. A. Distributed physiology and the molecular basis of social life in eusocial insects. Hormones Behav. 122, 104757 (2020).
Feng, Tao & Kang, Yun. Foraging dynamics in social insect colonies: mechanisms of backward bifurcations and impacts of stochasticity. Math. Biosci. 384, 109436 (2025).
Opachaloemphan, C., Yan, H., Leibholz, A., Desplan, C. & Reinberg, D. Recent advances in behavioral (epi) genetics in eusocial insects. Annu. Rev. Genet. 52(1), 489–510 (2018).
Zhou, L. & Fan, M. Dynamics of an sir epidemic model with limited medical resources revisited. Nonlinear Anal. Real World Appl. 13(1), 312–324 (2012).
Muhammad, N. & Eberl, H. J. Two routes of transmission for nosema infections in a honeybee population model with polyethism and time-periodic parameters can lead to drastically different qualitative model behavior. Commun. Nonlinear Sci. Numer. Simul. 84, 105207 (2020).
Funding
The authors acknowledge support from the Fundamental Research Grant Scheme with Project Code: FRGS/1/2022/STG06/USM/02/1 by the Ministry of Higher Education, Malaysia (MOHE).
Author information
Authors and Affiliations
Contributions
A. M. S., Contributed to the conception and design of the study, wrote the original version of the article, and performed the analysis and simulations. M. H. M., Supervised the study, contributed to its conception and design, and assisted in writing the original version of the article and funding. A. K. K. A., Contributed to reviewing, and editing, and assisted with the biological concept in the original manuscript. K. M. A. M., Contributed to the writing, validation, investigation, formal analysis, and funding.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
About this article
Cite this article
Salman, A.M., Mohd, M.H., Almansoori, A.K.K. et al. Nonlinear dynamics of Nosema ceranae and the fragile resilience of honeybee colonies under environmental strain.
Sci Rep (2026). https://doi.org/10.1038/s41598-026-45351-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-026-45351-1
Keywords
Nosema ceranae dynamics- Mathematical ecology of pollinators
- Bifurcation and stability analysis
- Colony collapse and resilience
- Eco-epidemiological modelling
Source: Ecology - nature.com
