Honey Bee Colonies Research Articles

Toxicological Impact Of Lambda-Cyhalothrin On Apis Mellifera: Comparative Analysis Under Semi-Field And Field Conditions

Abstract Lambda-cyhalothrin is the most widely used insecticide worldwide. This study investigates the impact of lambda-cyhalothrin on honeybees (Apis mellifera) colonies under semi-field and field conditions utilizing the formulated product (Karate 5% at a concentration of 0.85 ml/l). Considering the application of lambda-cyhalothrin in onion fields, the foraging activities reduced by 27.89 and 29.28% under semi-field and field conditions, respectively. The number of dead workers in the comparison group was three-fold higher than that in the control group under both semi-field and field conditions. Under semi-field and field conditions, reductions were observed in the brood (22.6 and 24.01%), pollen (94.46 and 41.51%), and nectar (69.27 and 44.50%) areas. The effect of lambda-cyhalothrin on the strength of A. mellifera under semi-field conditions resulted in a significant decrease until the 14th day (2.20 bee frames/hive) and in field condition it was on the 7th day (4.79 bee frames/hive). This highlights the non-bee-safe nature of lambda-cyhalothrin when applied during flowering, and spraying should be avoided to reduce its contamination of nectar and pollen. Overall, our results show the significance of conducting field trials to accurately assess the risk of pesticides, depicting a more realistic exposure scenario under actual colony conditions.

Molecular Detection and Phylogenetic Relationships of Honey Bee-Associated Viruses in Bee Products.

In the last few years, the isolation and amplification of DNA or RNA from the environment (eDNA/eRNA) has proven to be an alternative and non-invasive approach for molecular identification of pathogens and pests in beekeeping. We have recently demonstrated that bee pollen and bee bread represent suitable biological material for the molecular identification of viral RNA. In the present study, we extracted total RNA from different bee products (pollen, n = 25; bee bread, n = 17; and royal jelly, n = 15). All the samples were tested for the presence of six of the most common honey bee-associated viruses-Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Chronic bee paralysis virus (CBPV), Sacbrood virus (SBV), Kashmir bee virus (KBV), and Black queen cell virus (BQCV)-using a reverse transcription polymerase chain reaction (RT-PCR). We successfully detected six records of DWV (10.5%, 6/57), four of ABPV (7.0%, 4/57), three of Israeli acute paralysis virus (IAPV) (5.3%, 3/57), and two of BQCV (3.5%, 2/57). Using ABPV primers, we also successfully detected the presence of IAPV. The obtained viral sequences were analyzed for phylogenetic relationships with the highly similar sequences (megablast) available in the GenBank database. The Bulgarian DWV isolates revealed a high homology level with strains from Syria and Turkey. Moreover, we successfully detected a DWV strain B for the first time in Bulgaria. In contrast to DWV, the ABPV isolates formed a separate clade in the phylogenetic tree. BQCV was closely grouped with Russian isolates, while Bulgarian IAPV formed its own clade and included a strain from China. In conclusion, the present study demonstrated that eRNA can be successfully used for molecular detection of honey bee-associated viruses in bee products. The method can assist the monitoring of the health status of honey bee colonies at the local, regional, and even national levels.

Evaluating the efficiency of supplementary feeding as a management strategy for enhancing honeybee (Apis mellifera L.) colony growth and productivity

Sustaining honeybee colonies is challenging during dearth periods as their metabolic functions are reduced due to limited foraging activities. The experiment used honeybee colonies of Apis mellifera, and five different low-cost supplementary foods—sugar, banana, pumpkin, maize flour, and rice flour syrups—were introduced as treatments. Every box for each treatment received a daily 300-ml supplementary food syrup consisting of a specific amount of feeding materials along with 100 g of brown sugar and 20 g of honey. The amount of food consumed was assessed on the second day following the supplementation. Supplemental food with low-cost feeding materials significantly impacts the growth and strength of the colonies. Results revealed significant impacts on colony growth and strength, with all supplements contributing to food consumption over 78%. Despite variations in brood and pollen cells, all feeding supplements showcased efficiency in supporting honeybee feeding, indicating their potential utility in mitigating the challenges during the dearth period. Notably, pumpkin syrup emerged as the best supplement, offering cost-effectiveness compared to sugar and banana syrups, and it could reduce sugar syrup costs by 50% while enhancing brood, honey, and pollen cell production by 71.36%, 108.36%, and 58.73%, respectively. The findings of the economic analysis revealed that the cost of feeding materials was the highest for sugar syrup ($1.89), followed by banana ($0.91), pumpkin ($0.83), maize ($0.53), and rice ($0.53). This study suggests that supplementing honeybee colonies with low-cost feeding materials can positively impact colony growth and strength during dearth periods and advance the beekeeper’s decision as a cost-effective alternative to traditional sugar syrup.

MSPB: a longitudinal multi-sensor dataset with phenotypic trait measurements from honey bees

We present a one-year-long multi-sensor dataset collected from honey bee colonies (Apis mellifera) with rich phenotypic measurements. Data were collected non-stop from April 2020 to April 2021 from 53 hives located at two apiaries in Québec, Canada. The sensor data included audio features, temperature, and relative humidity. The phenotypic measurements contained beehive population, number of brood cells (eggs, larva and pupa), Varroa destructor infestation levels, defensive and hygienic behaviors, honey yield, and winter mortality. Our study is amongst the first to combine a wide variety of phenotypic trait measurements annotated by apicultural science experts with multi-sensor data, which facilitate a broader scope of analysis. We first summarize the data collection procedure, sensor data pre-processing steps, and data composition. We then provide an overview of the phenotypic data distribution as well as a visualization of the sensor data patterns. Lastly, we showcase several hive monitoring applications based on sensor data analysis and machine learning, such as winter mortality prediction, hive population estimation, and the presence of an active and laying queen.

Effect of artificial sugar supplement on the lifespan and learning memory ability of honey bee (Apis cerana) worker bee offspring.

Honeybees maintain their growth and reproduction mainly by collecting nutrients from nectar-source plants. Apis cerana, a unique species of honeybee in China, is capable of sporadically collecting nectar. In traditional beekeeping, sugar syrup or a honey-water solution must be artificially fed to bees to supplement their diet during rainy weather or nectar-deficient periods. In this study, 2 groups of honeybee colonies were each fed sugar syrup or a honey-water solution, and a third group consisting of colonies that were allowed to naturally forage without any dietary supplement was used as the control. The effects of the 2 sugar sources on A. cerana worker bee offspring were compared. The results showed that the sugar source affected the lifespan and learning memory of the worker bee offspring. The lifespan, learning memory ability, and expression of related genes in the sugar syrup group were significantly lower than those in the honey-water solution and natural nectar foraging groups (P < 0.05). A honey-water solution supplement was more beneficial to the healthy development of worker bee offspring than a sugar syrup supplement when the colonies lacked dietary resources. These findings provide a theoretical basis that can guide beekeepers in choosing the appropriate dietary supplements for honeybees.

Prevalence and molecular analysis of some important viruses in honey bee colonies in Türkiye: the status of multiple infections.

In this study, seven bee viruses of significant importance for bee health in Türkiye were investigated using one-step RT-PCR. For this purpose, larvae from 1183 hives and adult bees from 1196 hives were sampled from 400 apiaries in 40 provinces. The prevalence of viral infections in hives was as follows: acute bee paralysis virus (ABPV), 6.4%; black queen cell virus (BQCV), 77%; chronic bee paralysis virus (CBPV), 3.2%; deformed wing virus (DWV), 63.8%; Israel acute bee paralysis virus (IAPV), 7%; Kashmir bee virus (KBV), 2.7%; sacbrood virus (SBV), 49.7%. Moreover, 50 different combinations of viral infections were identified in the hives. While dual infections (36.1%) were the most common in hives, triple infections with BQCV, DWV, and SBV were found to have the highest prevalence (22.1%). At least one viral infection was detected in all of the apiaries tested. Phylogenetic analysis showed that the isolates from this study generally exhibited the highest similarity to previously reported Turkish isolates. When similarity ratios and the locations and types of amino acid mutations were analyzed, it was observed that the isolates from our study exhibited high similarity to isolates from various countries, including China, the United Kingdom, Syria, and Germany.

Safety assessment of high doses of vaporized oxalic acid on honey bee worker health and queen quality

IntroductionThe honey bee ectoparasitic mite, Varroa destructor, is one of the main causes of honey bee colony loss worldwide. Synthetic acaricides are the most commonly used strategy for varroa control, however, resistance to these acaricides has emerged. Consequently, the use of organic acids for varroa control is gaining more interest among beekeepers. For example, oxalic acid (OA) is a natural compound that has been shown to be an effective acaricide against varroa mites, however, the potential toxicity of OA to adult bees and queens is poorly understood. The objective of the study was to evaluate the toxicity of incremental doses of vaporized OA on honey bee workers and queens.MethodsWe exposed 32 colonies to incremental doses (0, 5, 10 or 20 g per colony) of vaporized OA once per week over four consecutive weeks and we monitored the acute and long-term toxicity. We investigated the short-term effects of OA administration by evaluating adult bee mortality, brood production, and population size. Next, we evaluated the long-term effects of OA application on both worker bees and queens. Regarding workers, we investigated their ability to rear new queens. As for queens, we measured acceptance, performance, and sperm quality.ResultsWe found that colonies treated with 20 g OA (20 times the label dose) had a statistically significant increase in worker bee mortality, with a non-significant, 23% decrease in brood relative to controls. No significant differences were observed in queen performance nor sperm quality.DiscussionWe found that repeated application of vaporized OA, at up to 20 times the label dose, had no significant short-term nor long-term, negative effects on colony or queen health, with the exception of a short-term increase in adult bee mortality in the 20 g OA-treated group. The results of this study support the safety of higher-than-label doses of OA for honey bee colonies. The observed increase in adult bee mortality in the 20 g OA dose group in this study suggests that OA doses should be maintained below 20 g per brood chamber.

Effect of Media on the Mycelial Growth of Honey Bee Fungal Pathogen &lt;i&gt;Ascosphaera apis&lt;/i&gt; Mating Types

The fungal pathogen Ascosphaera apis significantly impacts the honey bee Apis mellifera by causing larval mortality and reducing colony productivity. Understanding the growth dynamics of A. apis under different media conditions is crucial for advancing research on its biology and pathogenicity. This study evaluates the effect of eight artificial media namely Potato Dextrose Agar (PDA), Carrot Dextrose Agar (CDA), Beetroot Dextrose Agar (BDA), Tender Coconut Water (TCW) media, Sabouraud Dextrose Agar (SDA), Sabouraud Dextrose Agar Yeast (SDAY), Sabouraud Maltose Agar Yeast (SMAY), and Malt Extract Agar (MEA) on the mycelial growth and reproductive structure formation of A. apis's two mating types. Results indicated significant differences in mycelial growth across media. SDA (75.40 mm dia) and SDAY (74.83 mm dia) supported maximum growth for mating type 1 (MAT1) and, SDA (72.72 mm dia) and BDA (68.17 mm dia) yielded the maximum growth for mating type 2 (MAT2). The TCW resulted in the least growth for MAT1 (44.19 mm dia) and MAT2 (41.25 mm dia). The growth patterns showed that media supplemented with sugars generally enhanced A. apis growth, whereas non-sugar media did not. Furthermore, reproductive structures formed faster on PDA and CDA than on other media. The study confirms that media composition, particularly sugar content, critically influences the growth and development of A. apis, providing insights that could help to develop better management strategies for controlling chalkbrood disease in honey bee colonies.

Morphological and Molecular Identification and Mating Type Detection of Chalkbrood Fungal Pathogen &lt;i&gt;Ascosphaera apis&lt;/i&gt; in &lt;i&gt;Apis mellifera&lt;/i&gt; L. in Southern India

Chalkbrood Ascosphaera apis is a fungal brood pathogen that exerts considerable biotic stress on honey bees worldwide. Chalkbrood was noticed at Coimbatore in Apis mellifera colonies for the first time. Precise identification of the fungal pathogen and their mating types is crucial to develop effective disease management strategies. Hence, the present study aimed to isolate and identify fungal pathogen morphologically and genomically. Scanning Electron Microscopic (SEM) analysis was done and confirmed the essential morphological characteristics of the isolated fungus such as spore cysts (51-86 µm in diameter), spore balls (8-15 µm in diameter), and ascospores (1.7-2.6 µm in length). Molecular characterization using internal transcribed sequence (ITS)-PCR of the fungal isolate indicated 99-100% sequence similarity to A. apis. In addition, the multiplex PCR assay was performed and the mating types MAT1-1 and MAT1-2 were successfully detected and named TNAU CBD MAT1 and TNAU CBD MAT2, respectively.

Phenotypic and genotypic antimicrobial resistance patterns in honey bee (Apis mellifera L.) bacterial symbionts.

Antimicrobial resistance (AMR) is a major global public health problem. Nevertheless, the knowledge of the factors driving the spread of resistance among environmental microorganisms is limited, and few studies have been performed worldwide. Honey bees (Apis mellifera L.) have long been considered bioindicators of environmental pollution and more recently also of AMR. In this study, 53 bacterial strains isolated from the body surface of honey bees at three ontogenetic stages, collected from ten different geographic locations, were tested for their phenotypic and genotypic resistance to eight classes of the most widely used antimicrobials in human and veterinary medicine. Results showed that 83% of the strains were resistant to at least one antimicrobial and 62% were multidrug-resistant bacteria, with a prevalence of resistance to nalidixic acid, cefotaxime, and aztreonam. A high percentage of isolates harbouring at least one antimicrobial gene was also observed (85%). The gene encoding resistance to colistin mcr-1 was the most abundant, followed by those for tetracycline tetM and tetC. Geographical features influenced the distribution of these traits more than bacterial species or bee stage, supporting the use of honey bee colonies and their associated bacteria as indicators to monitor environmental resistance. This approach can improve the scientific understanding of this global threat by increasing data collection capacity.

Evaluation of a new plant-based formulation to control Varroa mite (Varroa destructor) in honey bee (Apis mellifera) colonies

Evaluation of a new plant-based formulation to control Varroa mite (Varroa destructor) in honey bee (Apis mellifera) colonies

First record of the sap beetle Phenolia picta invading Eastern honey bee (Apis cerana) (Hymenoptera, Apidae) colonies

In Asia, the sap beetle Phenolia (Lasiodites) picta (Macleay) (Coleoptera, Nitidulidae) (hereafter P. picta) is endemic to the southeastern Palaearctic and northeastern Oriental regions. It is commonly found in orchards and plantations, where it feeds on ripe or rotten fruits. Here, P. picta was recorded for the first time to invade the hive of the Eastern honey bee (Apis cerana). Both DNA barcode and morphological data were used to identify the beetle, with the DNA barcoding comparison to the P. picta sequence in the BOLD library reaching 100% similarity. Additionally, the body sizes were measured, with the males averaging 7.87 ± 0.21 mm in length and the females averaging 8.27 ± 0.15 mm. The beetle’s body colour is primarily brownish with yellowish to pale orange spots on the elytra. P. picta invades the honey bee hive, feeds on honey and bee pupae, and tunnels through the comb, which causes the honey bees to abandon their hive. Such invasions could have severe impacts on honey bee colonies. This report provides essential information for beekeepers to use when identifying P. picta, which could be useful for preventing and controlling parasitic beetles.

Assessing the Efficacy of Organic and Chemical Treatments for Varroa Mite Control in Macedonian Apiaries

Varroa destructor infestations pose a significant global threat to honeybee colonies, affecting bee populations and honey production. The purpose of this study is to assess the performance of Apis mellifera L. honeybee colonies under various control methods with a particular focus on critical parameters like mite dewdrop, honey production and colony growth. The study took place in a normal apiary in the Pelister district for a year from April 2021 to March 2022. There were three treatment groups: Group A got Amitraz (a chemical treatment), group M got formic acid (an organic treatment) and group K got nothing (a control group) to see how the treatments affected the bees. Regular monitoring assessed parameters including dropped mites and Varroa infestation levels with honey production evaluated seasonally. The results showed that both chemical and organic treatments successfully decreased the levels of Varroa mite infestation compared to the control group, while also noting differences in colony growth and honey production. The Amitraz treatment did prevent the mites from dropping, but it also slowed down the young's development. On the other hand, the formic acid treatment kept the brood's development healthy, suggesting that it might be a better way to get rid of varroa. This study shows how important it is to think about how treatments affect important colony parameters. It also calls for more research like longterm observations, to fully understand how treatments affect the health, development and longevity of bee families which will eventually lead to the creation of more effective combined pest control methods.

Effects of Vitamin Supplements in a Pollen Substitute on Some Characteristics of Bee Nucleus Colonies

This study aimed to evaluate the effects of vitamin supplementation in a pollen substitute on the performance of honey bee colonies (Apis mellifera mellifera). Twenty nucleus colonies housed in Langstroth-Ruth hives were randomly selected and fed a diet consisting of corn gluten, sugar, and brewer’s yeast residuals, yielding 21.0% crude protein, 7.1% crude fat, and 1.7% crude ash. This diet was enriched with a vitamin mixture at doses of 1, 2, and 3 grams. Control group I received a sugar-water syrup (1:1, w/v), while control group II received the basic diet without vitamin supplementation. The three experimental groups were given the diet with varying amounts of the vitamin mixture. Results indicated that experimental group II demonstrated colony strength comparable to experimental group I but significantly higher than the control group by 14.6-49.4%, and 34.3% higher than experimental group III (P&lt;0.05). Queen bee egg yield increased significantly, with control group II showing a 27.5% increase, experimental group I showing 52.1% increase, experimental group II showing 67.6% increase, and experimental group III showing 28.5% increase on July 28th. Similar trends were observed on August 28th, with respective increases of 30.2%, 32.5%, 37.4%, and 14.7% compared to control group I. Additionally, honey yields for experimental group II were significantly higher by 25.8-57.9% compared to the control groups and 46.8% higher than experimental group III (P&lt;0.05). These findings suggest that vitamin-enriched pollen substitutes positively impact colony strength, queen bee egg laying, and honey production, underscoring the potential benefits of such supplementation in beekeeping practices.

Diversity of plant pollen sources, microbial communities, and phenolic compounds present in bee pollen and bee bread.

The pollination of several crops, as well as wild plants, depends on honeybees. To get the nutrients required for growth and survival, honeybee colonies are dependent on pollen supply. Bee pollen (BP) is partially packed in honeycomb cells and processed into beebread (BB) by microbial metabolism. The composition of pollen is highly variable and is mainly dependent on ecological habitat, geographical origin, honey plants, climatic conditions, and seasonal variations. Although there are important differences between the BP and the BB, little comparative chemical and microbiological data on this topic exists in the literature, particularly for samples with the same origin. In this study, BP and BB pollen samples were collected from two apiaries located in the Campania and Molise regions of Southern Italy. Phenolic profiles were detected via HPLC, while antioxidant activity was determined by ABTS·+ and DPPH· assay. The next-generation sequencing (NGS) based on RNA analysis of 16S (rRNA) and internal transcribed spacer (ITS2) regions were used to investigate the microbial community (bacteria and fungi) and botanical origin of the BP and BB. Chemical analysis showed a higher content of flavonols in BP (rutin, myricetin, quercetin, and kaempferol), while in BB there was a higher content of phenolic acids. The NGS analysis revealed that the microbial communities and pollen sources are dependent on the geographical location of apiaries. In addition, diversity was highlighted between the microbial communities present in the BP and BB samples collected from each apiary.

First metagenomic analysis of virome in Uzbekistan honey bee (Apis mellifera): Investigating basic information on honey bee viruses

Honey bees are economically important insects. However, they face multiple biotic and abiotic stresses, such as diseases, pesticides, climate change, and pests, which cause the loss of honey bee colonies worldwide. Among these factors, viruses have been identified as the major cause of colony loss. Research on honey bee viruses in Uzbekistan is limited. This study investigated the viruses affecting honey bees in Uzbekistan. Virome analysis was conducted for each sample using high-throughput sequencing and bioinformatics. Nine honey bee viruses have been identified: the acute bee paralysis virus, aphid lethal paralysis virus, Apis rhabdovirus 1 and 2, black queen cell virus, deformed wing virus, Lake Sinai virus 10, sacbrood virus, and Hubei partiti-like virus 34. Additionally, 15 plant viruses were identified, 7 of which were novel. This study is the first virome analysis of Uzbekistan honey bees and provides a foundation for understanding the viruses affecting honey bees and plants in Uzbekistan.

A machine learning-based multiclass classification model for bee colony anomaly identification using an IoT-based audio monitoring system with an edge computing framework

Approximately one-third of human food is derived from flowering plants, and 80% of those plants require honeybees for pollination. However, climate change is affecting bees as well as the beekeeping industry. Having continuous information on the condition of honey bee colonies would be helpful in studying new diseases such as colony collapse disorder, as well as in developing novel beekeeping tools to improve the health management of bee colonies. In this study, an audio monitoring system is established based on Internet of Thing (IoT) technology to record the sounds produced by bees. A colony status classification model coupled with a machine learning algorithm is used to analyze the sound data of bee colonies under various treatments to classify the colony conditions, such as queenless, virus-infected, and pesticide-contaminated. With an edge computing framework, the classification is done by the monitoring system, which greatly reduces data transmission delay and cloud computing burden. The field test results show that the accuracy andF1-score for the classification model for the beehive under various treatments reach 98.8% and 0.98, respectively.The proposed IoT-based audio monitoring system for beehives can assist beekeepers and users in managing bee colonies more effectively by providing automatic, accurate, and real-time monitoring information regarding the current colony conditions.

Effects of landscape variation on thermoregulation and performance in Apis mellifera honey bee colonies: insights from mtDNA haplotypes

In a year-long research study, we delved into the thermoregulation and performance of twenty-one honey bee colonies situated in three distinct landscape settings: 1- Urban, 2- Agricultural (AG), and 3- Non-agricultural (Non-AG) environments. All colonies, consistent in size, were equipped with sensors designed to record inner hive temperature and humidity. Additionally, colonies were placed on precision electronic scales for weight monitoring. Varroa mite infestations were documented at seven different time points throughout the experiment. Colonies were genetically analyzed for subspecies and haplotype identification using the in silico DraI mtDNA COI-COII test. Our results unveiled an overall significant (p < 0.001) weight gain (5.76 kg) in colonies located within the urban setting, in stark contrast to both AG and Non-AG locations. Notably, colonies situated in the natural refuge center (Non-AG) experienced a weight loss of −0.05 kg, signaling a lack of sustainable forage resources in such a setting. Similarly, the inner hive temperature was significantly (p < 0.001) higher in the urban area (33.65 °C) with a significantly lower relative humidity (61.2%) compared to both other locations. Among the locations, the temperature variable displayed the strongest significant correlation (p < 0.001) with a positive correlation coefficient (r = 0.95). Two distinct mtDNA haplotypes were identified: C1 and C2j, characterizing the ligustica and carnica maternal origin, respectively. These haplotypes consistently demonstrated significant differences in their thermoregulation and weight gain, irrespective of the landscape compositions. The C1 haplotype displayed significantly higher weight gain alongside lower inner hive temperature and humidity levels compared to the C2j haplotype. No differences in mite infestation were observed between the various locations or haplotypes. Our data underscores the substantial impact of landscape composition on honey bee performance. It also unveils significant variations in thermoregulation and adaptation between haplotypes, shedding light on the multifaceted relationship between environmental factors, genetics, and bee colony dynamics.

The salivary gland transcriptome of Varroa destructor reveals suitable targets for RNAi-based mite control.

The mite Varroa destructorAnderson and Trueman (Mesostigmata: Varroidae) has a dramatic impact on beekeeping and is one of the main causes of honey bee colony losses. This ectoparasite feeds on honey bees' liquid tissues, through a wound created on the host integument, determining weight loss and a reduction of lifespan, as well as the transmission of viral pathogens. However, despite its importance, the mite feeding strategy and the host regulation role by the salivary secretions have been poorly explored. Here, we contribute to fill this gap by identifying the salivary components of V. destructor, to study their functional importance for mite feeding and survival. The differential expression analysis identified 30 salivary gland genes encoding putatively secreted proteins, among which only 15 were found to be functionally annotated. These latter include proteins with putative anti-bacterial, anti-fungal, cytolytic, digestive and immunosuppressive function. The three most highly transcribed genes, coding for a chitin-binding domain protein, a Kazal domain serine protease inhibitor and a papain-like cysteine protease were selected to study their functional importance by reverse genetics. Knockdown (90%-99%) by RNA interference (RNAi) of the transcript of a chitin-binding domain protein, likely interfering with the immune reaction to facilitate mite feeding, was associated with a 40%-50% decrease of mite survival. This work expands our knowledge of the host regulation and nutritional exploitation strategies adopted by ectoparasites of arthropods and allows the identification of potential targets for RNAi, paving the way towards the development of new strategies for Varroa mite control.

Complex dynamical behaviors of a honeybee-mite model in parameter plane

The honeybee has a significant impact on ecosystem stability, biodiversity conservation, and pollinating crops. However, during the past few years, there has been a sharp reduction in both the honeybee population and their colonies. It has been found that the parasitic mite Varroa destructor is responsible for the decline in honeybee colonies around the world, which in turn immensely affects the economic growth of a country. To investigate how the dynamics of a system is influenced by the growth of honeybees and the parasitism of mites, we study a nonlinear honeybee-mite population model in a discrete-time setup. We observe that an increase in the value of the queen’s egg-laying rate drives the system towards chaos. However, chaos can be controlled as well if the parasite attachment effect is increased. The intrinsic dynamical properties of the proposed system are investigated with the simultaneous variation of the queen’s egg-laying rate and the mite’s parasite attachment effect by constructing several largest Lyapunov exponent and isoperiodic diagrams. The investigation reveals the existence of several periodic structures in the quasiperiodic and chaotic regimes of the parameter plane, including Arnold tongues, saddle area, spring area, connected shrimp-shaped structure, and connected saddle area. We also find the appearance of a novel ‘jellyfish’-shaped periodic structure. One of the most fascinating findings of this study is the appearance of Arnold tongues along the inner boundary region of another Arnold tongue. In addition, this work also reveals different types of multistability, e.g., the coexistence of two, three, and even four attractors. What is more interesting is that the current analysis unveils the coexistence of five attractors as well, more specifically, four different periodic attractors coexist with the trivial fixed point attractor, which is quite rare in ecological systems. The structures of the basins of these coexisting attractors are either smooth or very complex in nature. Furthermore, the present study also discloses the fact that variation in the initial condition of the system can significantly change the appearance of the periodic structures in the parameter plane.