==== Front Environ Entomol Environ Entomol ee Environmental Entomology 0046-225X 1938-2936 Oxford University Press US 37170880 10.1093/ee/nvad035 nvad035 Pest Management AcademicSubjects/SCI01382 Interactions between integrated pest management, pollinator supplementation, and normalized difference vegetation index in pumpkin, Cucurbita maxima (Cucurbitales: Cucurbitaceae), production https://orcid.org/0000-0001-8394-1284 Waithaka Nduta A Department of Agricultural Science and Technology, Kenyatta University, P.O. Box 43844-00100, Nairobi, Kenya Apiculture and Beneficial Insects Research Institute (ABIRI), Kenya Agricultural and Livestock Research Organization (KALRO), P.O. Box 32-30403, Marigat, Kenya https://orcid.org/0000-0002-4983-4546 Kasina Muo Apiculture and Beneficial Insects Research Institute (ABIRI), Kenya Agricultural and Livestock Research Organization (KALRO), P.O. Box 32-30403, Marigat, Kenya https://orcid.org/0000-0003-3399-1616 Samita Namikoye E Department of Agricultural Science and Technology, Kenyatta University, P.O. Box 43844-00100, Nairobi, Kenya Guantai Mary M Kenya Plant Health Inspectorate Service (KEPHIS), P.O. Box 19164-00501, Nairobi, Kenya https://orcid.org/0000-0003-4476-3467 Omuse Evanson R International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya https://orcid.org/0000-0002-8463-8807 Toukem Nadia K Department of Zoology and Entomology, University of Pretoria, Hatfield, 0028 Pretoria, South Africa https://orcid.org/0000-0002-8603-6332 Lattorff H Michael G International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya Abdel-Rahman Elfatih M International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya https://orcid.org/0000-0003-3808-9920 Adan Marian International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya Mohamed Samira A International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya Dubois Thomas International Centre for Insect Physiology and Ecology (icipe), P.O. Box 30772-00100, Nairobi, Kenya Schmidt-Jeffris Rebecca Subject Editor Corresponding author, mail: muo.kasina@kalro.org Present address: Department of Chemistry, University of Nairobi, P.O. Box 30197-00100, Nairobi, Kenya 6 2023 12 5 2023 12 5 2023 52 3 416425 09 12 2022 11 3 2023 12 4 2023 © The Author(s) 2023. Published by Oxford University Press on behalf of Entomological Society of America. 2023 https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Sustainable production of pumpkin (Cucurbita maxima Duchesne) partly relies on integrated pest management (IPM) and pollination services. A farmer-managed field study was carried out in Yatta and Masinga Sub-Counties of Machakos County, Kenya, to determine the effectiveness of a recommended IPM package and its interaction with stingless bee colonies (Hypotrigona sp.) for pollinator supplementation (PS). The IPM package comprised Lynfield traps with cuelure laced with the organophosphate malathion, sprays of Metarhizium anisopliae (Mechnikoff) Sorokin isolate ICIPE 69, the most widely used fungal biopesticide in sub-Saharan Africa, and protein baits incorporating spinosad. Four treatments—IPM, PS, integrated pest and pollinator management (which combined IPM and PS), and control—were replicated 4 times. The experiment was conducted in 600 m2 farms in 2 normalized difference vegetation index (NDVI) classes during 2 growing seasons (October 2019–March 2020 and March–July 2020). Fruits showing signs of infestation were incubated for emergence, fruit fly trap catches were counted weekly, and physiologically mature fruits were harvested. There was no effect of IPM, PS, and NDVI on yield across seasons. This study revealed no synergistic effect between IPM and PS in suppressing Tephritid fruit fly population densities and damage. Hypotrigona sp. is not an efficient pollinator of pumpkin. Therefore, we recommend testing other African stingless bees in pumpkin production systems for better pollination services and improved yields. Graphical Abstract Graphical Abstract Cucurbita maxima integrated pest and pollinator management (IPPM) Hypotrigona sp Zeugodacus cucurbitae German Federal Ministry of Economic Cooperation and Development Deutsche Gesellschaft für Internationale Zusammenarbeit 10.13039/501100011099 International Agricultural Research 17.7860.4–001 Norwegian Agency for Development Cooperation, the Section for Research, Innovation, and Higher Education RAF-3058 KEN-18/0005 Swedish International Development Cooperation Agency 10.13039/100004441 Swiss Agency for Development and Cooperation 10.13039/100009131 Federal Democratic Republic of Ethiopia; and the Government of the Republic of Kenya Government of the Republic of Kenya 10.13039/100013987 ==== Body pmcIntroduction Pumpkins, Cucurbita maxima Duchesne (Cucurbitales: Cucurbitaceae), produce fruits and leaves that are rich in proteins, carbohydrates, and oils containing unsaturated fatty acids (Indrianingsih et al. 2019). Their seeds contain minerals such as nitrogen, phosphorous, sodium, potassium, calcium, magnesium, iron, copper, zinc, and manganese (Nawirska et al. 2009). In Kenya, pumpkins are mainly grown by smallholder farmers for household use, while the surplus is traded for cash income (Machakos County Integrated Development Plan 2019). Pumpkins, therefore, provide affordable nutritious food to low-income families and curb food insecurity (Rődlach 2011). Production in the country is still low compared with neighboring countries, with the area under pumpkin production at 1,293 ha and 1,380 ha in 2017 and 2018, respectively (Horticultural Crops Directorate 2018). In Kenya, pumpkin production mainly occurs in semiarid areas, which comprise the largest area of arable land in the country (Ndengwa et al. 2016). Machakos County is a major production area, where households prefer fruit to leaf consumption (Ndegwa et al. 2016). Pumpkin production in Kenya is seriously constrained by pests, including the silverleaf whitefly, Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) (Messelink et al. 2008), and several species of Tephritid fruit flies (Diptera: Tephritidae) (De Meyer et al. 2012, Kambura et al. 2018). Tephritid fruit flies seem to be the most destructive pests of pumpkin fruits, particularly the melon fruit fly, Zeugodacus cucurbitae Coquillett (Diptera: Tephritidae), which can cause up to 100% yield loss without adequate control measures (Mwatawala et al. 2006, Deguine et al. 2012, Badii et al. 2015, Doorenweerd et al. 2018). In Kenya, this pest has been reported to cause very high yield losses to cucurbits, particularly pumpkin fruits (Kambura et al. 2016). Zeugodacus cucurbitae occurrence is directly influenced by abiotic factors, particularly temperature, which affect its development, survival, and distribution (Fletcher 1989, Aluja et al. 2014, Hill et al. 2016). Other important Tephritid fruit fly species associated with cucurbits include the greater pumpkin fruit fly, Dacus bivittatus Bigot; the lesser pumpkin fly, Dacus ciliatus Loew; the jointed pumpkin fly, Dacus vertebratus Bezzi; and the tomato fruit fly, Dacus punctatifrons Karsch (Okolle and Ntonifor 2005, Kambura et al. 2016, Layodé et al. 2019). The host preference of the different fruit fly species affecting cucurbits varies according to the environment, availability of host plants, and other competing pests (Charlery et al. 2017). Crop production, including pumpkin, is highly influenced by the agricultural landscape, which affects the distribution, diversity, and abundance of pest and beneficial arthropods because of their interaction with abiotic factors such as temperature, altitude, and rainfall (O’Rourke and Jones 2011). The normalized difference vegetation index (NDVI) has been widely used to understand landscape and vegetation phenology. Recent studies have reported the use of NDVI as a proxy for green biomass, which is related to canopy photosynthesis (Wang et al. 2017) and vegetation biodiversity (Beck et al. 2007). Vegetation biodiversity at landscape scales is key as it influences population densities of both pollinators and natural enemies, which use available resources in alternative habitats for their development (Shackelford et al. 2013). The abundance of natural enemies particularly may in turn influence pest population in the ecosystem (Copeland et al. 2006). Studies have reported how NDVI can predict the distribution and abundance of insect pests and pollinators and monitor long-term landscape structure (Abdel-Rahman et al. 2017, Toukem et al. 2020, 2022). For example, Toukem et al. (2020) reported a high density of Tephritid fruit flies infesting avocado, Persea americana Miller (Laurales: Lauraceae), in low NDVI. In this study, NDVI was used as a proxy for landscape vegetation. Tephritid fruit flies in cucurbit production systems are usually controlled by chemical insecticides (Dhillion et al. 2005), which have negative effects on human, animal, and environmental health, including pollinators, and may cause the development of insecticide resistance (Cloyd and Bethke 2015, Nobre et al. 2019). Furthermore, misuse of insecticides may lead to surpassing allowed residue limits and raising trade challenges (Badii et al. 2015, Mebdoua 2018). To reduce reliance on chemical pesticides, integrated pest management (IPM) approaches are being researched and deployed for Tephritid fruit fly control in Africa, largely on mango, Mangifera indica Linnaeus (Sapindales: Anacardiaceae) (Korir et al. 2015, Muriithi et al. 2016, 2020). The use of cultural practices, particularly field sanitation whereby infested fruits are removed, is commonly used for fruit fly management in Africa and can significantly reduce Tephritid fruit fly populations and minimize yield losses (Niassy et al. 2022). Sterile insect technique, which suppresses the pest population by maintaining a barrier of sterile male flies, has been used successfully in Japan against the melon fruit fly (Ito et al. 2003). Male annihilation technique consists of the deployment of high-density trapping stations with male-specific lures combined with an insecticide. The use of cuelure (the pheromone 4-(p-acetoxyphenyl)-2-butanone, which is highly attractive to male Z. cucurbitae) laced with fipronil, has been used successfully against Z. cucurbitae and fruit fly species in genera Bactrocera (Spafford et al. 2018). Among the biopesticides used in sub-Saharan Africa, Metarhizium anisopliae (Mechnikoff) Sorokin (Hypocreales: Clavicipitaceae) isolate ICIPE 69 (Real IPM, Thika, Kenya) is widely used as a fungal biopesticide (Akutse et al. 2020) and has been found to be highly potent against Z. cucurbitae (Onsongo et al. 2019). Protein baits, which are mainly female attractants, entail the application of food baits, combined with insecticides in localized spots against fruit flies, and have been effective in the management of Z. cucurbitae (Shikano et al. 2022). IPM against Z. cucurbitae has been developed and implemented across the Pacific regions (Vargas et al. 2015), but there are few reports on IPM practices for fruit fly management in pumpkin production systems in East Africa. Pumpkin is a pollination-dependent crop (Vidal 2010), with the honeybee Apis mellifera Linnaeus (Hymenoptera: Apidae) reported as the most important pollinator for cucurbits in Kenya (Njoroge et al. 2010). Often, there is a pollination deficit in Africa, as demonstrated for avocado (Sagwe et al. 2021), which can be solved through pollination supplementation with honeybees (Sagwe et al. 2021). The study area in Machakos is characterized as a semiarid area with limited natural habitats (Ndengwa et al. 2016). Limitations in pumpkin pollination in the area may be attributed to the reduced natural habits, which have been reported to reduce species richness and diversity of pollinators (Carvalheiro et al. 2010). However, African stingless bees (Hymenoptera: Apidae) might offer a better alternative for cucurbit pollination, especially for smallholders because of their ecological adaptability, easy domestication, floral constancy, perennial colonies, and polylectic nature (Heard 1999). Hypotrigona sp. was the preferred pollinator for the present study. This parasitic wasp is widely distributed across the tropics, lives in large colonies (Kiatoko et al. 2012), is a generalist in terms of food and nesting resources, and has been reported to pollinate other cucurbits such as the greenhouse cucumber Cucumis sativus Linnaeus (Cucurbitales: Cucurbitaceae), yielding 90% seed germination (Kiatoko et al. 2021b). The genus Hypotrigona contains several species (Hypotrigona gribodoi Magretti, Hypotrigona araujoi Michener, and Hypotrigona ruspolii Magretti in East Africa and Hypotrigona penna Eardley in West Africa) (Ndungu et al. 2017), which are reported to offer good pollination services in cucurbits (Kiatoko et al. 2021b). Studies carried out in Kenya have shown that Hypotrigona sp. successfully pollinate cucurbits such as cucumber and Galia muskmelon Cucumis melo L. var. reticulatus ser. (Cucurbitales: Cucurbitaceae) (Kiatoko et al. 2021a, 2021b). Other important African stingless bees that can be exploited for cucurbit pollination include Meliponula ferruginea Lepeletier (Hymenoptera: Apidae) and Meliponula bocandei Spinola (Hymenoptera: Apidae). These bees have been reported to successfully pollinate coffee Coffea canephora Pierre ex A. Froehner (Gentianales: Rubiaceae) in Uganda, with pollination efficiency of 98% and 95% fruit set, respectively (Munyuli 2010). Integrated pest and pollination management (IPPM) is a holistic strategy that synergizes pest management practices and pollination services in an ecologically and economically sustainable manner (Biddinger and Rajotte 2015, Merle et al. 2021) especially in pollinator-dependent crops that are highly affected by damaging insect pests, such as cucurbits (Leach et al. 2022). IPPM focuses on managing pests while simultaneously mitigating the negative impact of control methods on pollinator health (Biddinger and Rajotte 2015). The aim of this study was, therefore, to determine the effect of IPM and pollinator supplementation (PS), and their interaction, on Tephritid fruit fly population densities and damage, and pumpkin yield and assess whether landscape vegetation influences these effects. Materials and Methods Study Site The study was carried out in farmer-managed pumpkin fields in Yatta (1°24ʹ S, 37°59ʹ E) and Masinga Sub-Counties (0°58’ S, 37°36’ E) in Machakos County within low and medium normalized difference vegetation index (NDVI) classes (Fig. 1). Machakos County covers a total area of 6,208 km2 and lies at an elevation of 790–1,594 m above sea level. The County receives an average annual rainfall of between 500 and 1,250 mm, with October–November as the main rainy season (Machakos County Integrated Development Plan 2019). Fig. 1. Distribution of pumpkin farms in Yatta and Masinga Sub-Counties, Machakos County. Computation of NDVI Rainfall data from the Climate Hazard Group Infrared Precipitation with Station (CHIRPS, University of California, Santa Barbara, USA) were used at a 5 × 5 km spatial resolution to determine the exact dry and wet seasons in 2019 for the study area. Sentinel-2 (S2) satellite images were then utilized to compute the NDVI in the cloud-based Google Earth Engine platform (GEE, Alphabet, Mountain View, USA) according to Gorelick et al. (2017). One hundred and ninety-five S2 images were acquired for both the wet (February–May, 99 images) and dry (June–August, 96 images) seasons. We then created a composite image for each season using the multi-date S2 median reflectance pixel values at 10-m spatial resolution. A detailed description of the S2 image processing steps is found in Adan et al. (2021). Using Equation (1), we calculated the NDVI for both dry and wet composite images based on red and near-infrared bands. NDVI =(NIR−red)(NIR+red) (1) where NIR and red are reflectance values at near-infrared and red bands 8 and 4 of S2, respectively. The wet and dry season NDVI values were reclassified into 3 vegetation productivity classes (low, medium, and high) using the unsupervised K-means clustering method (Adan et al. 2021). The ranges of NDVI values for each class for the dry season were as follows: 0.320–0.077 (low), 0.077–0.133 (medium), and 0.133–0.619 (high), while the NDVI ranges for the wet season were as follows: 0.373–0.158 (low), 0.158–0.249 (medium), and 0.249–0.730 (high). Subsequently, a multiseason composite classified map was created by combining the outputs of the wet and dry season classifications based on the cluster values for both seasons (Adan et al. 2021). Experimental Layout Two experiments were carried out during the 2 pumpkin growing seasons (October 2019–March 2020 and March–July 2020). The first and second seasons are characterized by periods of long and short rains, respectively. The experiment consisted of 2 factors—(i) PS (absence or presence of managed Hypotrigona sp. colonies) and (ii) IPM (absence or presence of IPM)—conducted as a full-factorial design, yielding 4 treatments (control, PS, IPM, and IPPM, with IPPM representing the presence of both PS and IPM). A total of 32 farmer fields were included in the experiments. Each treatment was replicated 4 times in separate farmer fields, both in low and medium NDVI classes, yielding a total of 8 farmer fields per treatment (Fig. 1). Suitable farms were determined using a baseline survey carried out during July–August 2019. Treatments were allocated to pumpkin farmers based on their willingness to pay for the IPM, PS, or IPPM. Those willing to pay for both IPM and PS were randomly selected to host the IPPM treatment; those willing to pay for either IPM or PS only were randomly selected to host the IPM and PS treatments, respectively; and those not willing to pay for either package were randomly selected to host the control treatment. A minimal distance was kept among farms, depending on the treatment option, based on criteria that were developed and described by Adan et al. (2021). Specifically, (i) the distance between IPPM and PS farms was at least 1.5–3.0 km, (ii) the distance between IPM and control farms was at least 0.5 km, and (iii) the distance between IPPM or PS farms was at least 3.5 km away from either IPM or control. During the experiments, all plots were farmer-managed. Therefore, each farmer carried out similar agronomic practices such as land preparation, planting, weeding, and top-dressing with fertilizer independently. No plots were sprayed with insecticides during the entire period of the study. Planting and Application of IPM and PS Packages In each field, plots measuring 20 × 30 m used for growing pumpkins were demarcated. Each farmer was supplied with 150 pumpkin seeds (Dora F1 Hybrid, Safari Seeds, Nairobi, Kenya), which were sown directly in the soil at a spacing of 2 × 2 m and a depth of 3–4 cm. Three weeks from sowing, the crops were top-dressed with 15 g/plant of 17-17-17 NPK fertilizer (DMBL Ruiru, Nairobi, Kenya). The IPM package for this study comprised the following: (i) 2 Lynfield traps fitted with cuelure and cotton balls laced with the organophosphate malathion (Farmtrack Consulting, Nairobi, Kenya) were placed along a transect line in the middle of the farm running along the longest side of the farm. The traps were placed 50 m from each other and installed 8 weeks after crop germination (at the onset of flowering) (Toukem et al. 2020). A metallic stand was anchored 0.3 m into the ground and was used to hold the traps 1 m above the ground. The lures and malathion cotton balls in the traps were replaced 6 weeks after installation. (ii) One hundred milliliters protein bait made from brewer’s yeast (Fruit Fly Mania, Kenya Biologics, Makuyu, Kenya), laced with 0.2-ml spinosad (Tracer, Dow Chemical East Africa, Nairobi, Kenya) was added to a hand spray pump containing 900 ml of water, thoroughly mixed and applied weekly from the onset of flowering to harvesting. (iii) The fungal biopesticide Metarhizium anisopliae isolate ICIPE 69 (Mazao Campaign, Real IPM) was applied once at the onset of flowering by adding 10 ml of the product in 15 liters of water contained in a knapsack. Soil was drenched with the mixture on the entire plot. Unlike the IPM and IPPM plots where 2 Lynfield traps were installed, only one Lynfield trap baited with cuelure and malathion was installed in the control and PS plots to monitor fruit fly populations. Hypotrigona sp. bee colonies for PS treatment were obtained from stingless beekeeping farms in Mwingi, Kenya. The colonies were contained in stingless bee-specific hives, which were small rectangle wooden boxes (30.45 × 8.50 × 7.59 cm). They were constructed by local carpenters using local wood, with all spaces sealed except for a tiny hole (10 mm diameter), which the bees used to create their entrance/exit. The top cover was held with hinges and tightly fixed to avoid open spaces. Stingless bees were collected from the wild and allowed to settle in the hives before transferring them to the trial sites. During transportation, care was taken to avoid workers escaping by sealing the entrance/exit holes with mud (Kiatoko, personal communication). Beehives were installed in farms hosting IPPM and PS treatments at a density of 2 hives/farm at the onset of flowering. The hives were hung on trees within each farm at the edge of the long side of the plot (no separation distance was maintained since the stingless bees do not attack people) and kept throughout the study period. Data Collection Field data collection began 8 weeks after pumpkin germination and continued throughout the flowering and fruiting crop stages for 8 weeks. Fruit flies caught in the Lynfield traps were collected weekly in 75-milliliter plastic bottles and transported for further analysis to the entomology laboratory at National Sericulture Research Centre (NSRC), Kenya Agricultural and Livestock Research Organization (KALRO), Thika, Nairobi. Contents of the plastic bottles were emptied into 85 mm (diameter) × 15 mm (height) plastic Petri dishes containing 70% ethyl alcohol. Fruit flies were morphologically identified to species level according to the keys described by White (2006). The number of each species was recorded per sample. If in doubt, confirmatory morphological identification of the fruit flies was done at the International Centre of Insect Physiology and Ecology (icipe), Nairobi, Kenya by R. Copeland. Pumpkin fruits that had detached from the plant and had symptoms of fruit fly infestation (presence of oviposition marks and larvae) were collected weekly from the field, placed in 30 × 21 cm brown paper bags, and transported to NSRC. Damaged fruits were counted, weighed individually, and incubated. Incubation containers varied based on the size of the fruits: small fruits were incubated in 2-liter plastic containers, medium-sized fruits were incubated in 11-liter plastic containers, and large fruits were incubated in 20-liter plastic containers. All containers were lined with 305 × 305 mm serviettes and the rims were covered with 0.1 µm (pore size) muslin nets held in place with rubber bands. The containers were checked daily, and emerged puparia were collected, counted, and placed in 85 (diameter) × 15 mm (height) Petri dishes for emergence. Emerged adults were identified based on the above procedures and their counts recorded. On a weekly basis, the health of the hives was confirmed by verifying adequate numbers of worker bees leaving and entering the hives. Absconded hives were replaced. Sixteen weeks after crop germination (at maturity), all pumpkin fruits were harvested once in all plots. The number of healthy and damaged harvested fruits/plant and the weight/fruit were recorded. Data Analysis Fruit fly abundance was expressed as daily catches/trap/farm before analysis. Infestation of pumpkin fruits was expressed as the number of puparia per kilogram prior to analysis. Yield was expressed in kilogram per hectare. Trap catches for Z. cucurbitae, D. bivittatus, and D. punctatifrons were [log10(x + 1)] transformed to obtain a normal distribution with equal variance and subjected to a linear mixed-effect model to analyze the effects of IPM, PS, NDVI, and interactions, with the farm included as a random factor. However, very low catches of Bactrocera dorsalis Hendel (Diptera: Tephritidae), D. vertebratus, and D. ciliatus were obtained from trap catches; hence, they were not considered for subsequent analyses. The performance of the models was assessed by inspecting residual error distribution and checking for autocorrelation. A correlation structure was included in the model when there was evidence of autocorrelation (Durbin Watson test statistic < 1 or > 3). Wald chi-square tests were used to test fixed effects in the models, and when significant, means were separated using a Tukey HSD test. We performed all mixed models using the lme4 package (Bates et al. 2015), with the lmer function. However, the number of puparia per kilogram did not assume a normal distribution and was analyzed using the aligned rank transformation analysis of variance (ANOVA) to assess the effects of IPM, PS, NDVI, and interactions. The direct effect of IPM, PS, NDVI, and their interactions on the yield and number of fruits per plant were assessed using a linear model. All analysis was performed in R software version 4.0.5 (R Core Team 2020), and the significance was assessed at 5%. Results Fruit Fly Trap Catches Three fruit fly taxa belonging to the genera Zeugodacus, Dacus, and Bactrocera were identified from the trap catches. Fruit flies were identified as Z. cucurbitae, D. punctatifrons, D. bivittatus, D. vertebratus, D. ciliatus, and B. dorsalis. During the first season, there was no significant effect of either IPM, PS, or NDVI, or their interaction (χ2 ≤ 0.98, df = 1, 246, P ≥ 0.37) on the daily catches of Z. cucurbitae, which averaged 5.76 ± 0.42 catches/trap/day (Table 1). There was a main effect of PS in the second season (χ2 = 5.19, df = 1, 246, P = 0.023), with densities without PS (12.51 ± 1.17 catches/trap/day) higher than densities with PS (9.84 ± 0.90 catches/trap/day). The mean daily catches of Z. cucurbitae were not significantly different between farms with and without IPM (χ2 = 0.83, df = 1, 246, P = 0.36), or across NDVI (χ2 = 2.16, df = 1, 246, P = 0.14). In the second season, a significant interaction occurred among IPM, PS, and NDVI (χ2 = 4.82, df = 1, 246, P = 0.028), with the lowest densities recorded on IPPM farms in medium NDVI (5.41 ± 1.07 catches/trap/day) and the highest densities on control farms in low NDVI (19.47 ± 3.32 catches/trap/day). Table 1. Model estimates (standard error) from linear mixed-effects models fitted to estimate effects of pest control (IPM) (no, yes), PS (no, yes), and NDVI (low, medium), on Z. cucurbitae, D. bivittatus and D. punctatifrons abundance in Machakos County, Kenya Z. cucurbitae D. bivittatus D. punctatifrons Season 1 Season 2 Season 1 Season 2 Season 1 Season 2 (Intercept) 5.78*** 12.16*** 1.43*** 1.71*** 1.96*** 1.29*** (0.33) (0.26) (0.07) (0.12) (0.11) (0.05) IPM (yes) −0.87 −0.71 −0.84 −0.86 −0.97 −0.99 (0.47) (0.37) (0.09) (0.17) (0.15) (0.08) Stingless bee introduction (yes) −0.75 −0.43* −0.83 −0.71 −0.87 −0.86 (0.47) (0.37) (0.09) (0.17) (0.15) (0.08) NDVI (medium) −0.63 −0.58 −0.73** −0.74 1.00 −0.96 (0.47) (0.37) (0.09) (0.17) (0.15) (0.08) IPM (yes) × stingless bee introduction (yes) 1.02 2.76 1.58** 1.81* 1.11 1.08 (0.66) (0.53) (0.13) (0.24) (0.21) (0.11) IPM (yes) × NDVI (medium) 1.83 1.25 1.26 1.21 1.37 −0.99 (0.66) (0.53) (0.13) (0.24) (0.21) (0.11) Stingless bee introduction (yes) × NDVI (medium) 1.56 3.13* 1.28 1.81* 1.40 1.23 (0.66) (0.53) (0.13) (0.24) (0.21) (0.11) IPM (yes) × stingless bee introduction (yes) × NDVI (medium) −0.64 −0.20* −0.59** −0.43* −0.74 −0.82 (0.94) (0.74) (0.19) (0.34) (0.30) (0.15) Estimates were tested with the Wald χ2 test and appreciated at the 5% level of significance. Significance code: ***P < 0.001, **P < 0.01, *P < 0.05. In the first season, the mean daily catches of D. bivittatus differed significantly across NDVI (χ2 = 10.87, df = 1, 246, P < 0.0001), with densities in low NDVI (0.43 ± 0.05 catches/trap/day) higher than densities in medium NDVI (0.10 ± 0.01 catches/trap/day) (Table 1). However, no significant effect of IPM (χ2 = 3.28, df = 1, 246, P = 0.070) or PS (χ2 = 3.76, df = 1, 246, P = 0.053) was observed. A significant interaction among IPM, PS, and NDVI (χ2 = 7.87, df = 1, P = 0.005) was also reported. In low NDVI, IPPM farms recorded the highest densities (0.73 ± 0.16 catches/trap/day), while the lowest densities were reported on control farms in medium NDVI (0.06 ± 0.24 catches/trap/day). In the second season, there was a main effect of PS on mean daily catches of D. bivittatus (χ2 = 4.08, df = 1, 246, P = 0.043) with higher densities on PS farms (0.64 ± 0.08 catches/trap/day) than on farms without PS (0.52 ± 0.05 catches/trap/day). There was a significant interaction among IPM, PS, and NDVI (χ2 = 6.20, df = 1, 246, P = 0.013) on the mean daily catches of D. bivittatus, with highest densities on IPPM farms in low NDVI (1.13 ± 0.22 catches/trap/day) and lowest densities on PS farms in low NDVI (0.25 ± 0.07 catches/trap/day). During the first season, there was no significant effect of IPM, PS, NDVI, or their interaction on the daily catches of D. punctatifrons (χ2 ≤ 1.30, df = 1, 246, P ≥ 0.25), which averaged 1.38 ± 0.09 catches/trap/day (Table 1). In the second season, there was a main effect of PS on the daily catches of D. punctatifrons (χ2 = 4.09, df = 1, 246, P = 0.043) with higher densities on farms without PS (0.29 ± 0.03 catches/trap/day) than on PS farms (0.22 ± 0.05 catches/trap/day). However, no significant effects of IPM (χ2 = 0.00, df = 1, 246, P = 0.95) or NDVI (χ2 = 0.35, df = 1, 246, P = 0.56) were observed. There was no significant interaction among IPM, PS, and NDVI on D. punctatifrons catches in the second season (χ2 = 1.81, df = 1, 246, P = 0.18). Pumpkin Infestation Zeugodacus cucurbitae, D. bivittatus, D. vertebratus, and D. ciliatus were recovered from incubated pumpkin fruits. Their population densities were lumped together since their emergences were very low. In the first season, the number of puparia per kilogram of fruit did not differ significantly across IPM (F = 0.28, df = 1, 248, P = 0.60), PS (F = 0.24, df = 1, 248, P = 0.62), or NDVI (F = 0.21, df = 1, 248, P = 0.65) (Table 2). In the second season, IPM significantly (F = 4.88, df = 1, 248, P = 0.028) influenced pumpkin fruit infestation. There was, however, no main effect of NDVI (F = 3.57, df = 1, 248, P = 0.060) and PS (F = 0.16, df = 1, 248, P = 0.69) on infestation. A significant interaction among IPM, PS, and NDVI was observed (F = 5.81, df = 1, 248, P = 0.017), with lowest infestation on IPM farms in low NDVI and highest infestation on IPPM farms in low NDVI. Table 2. Interaction of IPM (no, yes), PS (no, yes), and NDVI (low, and medium) on Tephritid fruit fly infestation index using recovered puparia/kg of infested pumpkin fruits across 2 growing seasons in Machakos County, Kenya Puparia/kg (infestation index) IPM PS NDVI N Season 1 Season 2 No No Low 32 155.60 ± 88.87a 292.09 ± 142.59ab No Yes Low 32 110.91 ± 35.32a 131.38 ± 57.86ab Yes No Low 32 137.59 ± 78.20a 29.13 ± 19.39a Yes Yes Low 32 331.33 ± 207.76a 369.49 ± 248.94b No No Medium 32 17.63 ± 12.35a 92.80 ± 45.61ab No Yes Medium 32 320.79 ± 153.66a 71.73 ± 38.78ab Yes No Medium 32 180.71 ± 84.50a 174.42 ± 76.60b Yes Yes Medium 32 249.47 ± 168.16a 57.05 ± 34.43ab Different small letters in the same column indicate a significant difference at α = 0.05 according to the Tukey HSD test. Number of Fruits and Yield During the first season, there was no significant effect of IPM (F = 1.46, df = 1, 7.73, P = 0.24), PS (F = 0.40, df = 1, 7.73, P = 0.53), NDVI (F = 2.96, df = 1, 7.73, P = 0.10), or their interaction (F = 0.11, df = 1, 7.73, P = 0.75) on the number of fruits/plant, which averaged 1.02 ± 0.32 fruits/plant (Fig. 2). In the second season, IPM significantly (F = 10.12, df = 1, 6.56, P = 0.004) influenced the number of fruits/plant, with more fruits on IPM farms (1.58 ± 1.13 fruits/plant) than on farms without IPM (0.22 ± 0.08 fruits/plant). However, the number of fruits per plant did not differ significantly across NDVI (F = 0.29, df = 1, 6.56, P = 0.60) or PS (F = 0.13, df = 1, 6.56, P = 0.72). There was no significant interaction among IPM, PS, and NDVI on the number of fruits/plant (F = 2.22, df = 1, 6.56, P = 0.15). Fig. 2. Interaction between IPM (no, yes), PS (no, yes), and NDVI (low, medium), on the number of fruits/plant harvested during the second season in Machakos County, Kenya. In the first season, there was no main effect of IPM, PS, NDVI, and their interaction (F ≤ 1.54, df = 1, 6.57, P ≥ 0.23) on yield, which averaged 475.31 ± 133.17 kg/ha across treatments. During the second season, IPM, PS, NDVI, and their interaction had no significant effect on the yield (F ≤ 0.29, df = 1, 9.10, P ≥ 0.60), which averaged 180.78 ± 59.85 kg/ha. Discussion The sustainability of pumpkin production relies in part on IPM and PS to suppress key Tephritid fruit fly pests and augment pollination, ultimately increasing yield. During the study, Z. cucurbitae was reported as the key fruit fly species from the cuelure traps with minimal catches recorded for D. bivittatus and D. punctatifrons. These results are in line with findings from earlier studies, which reported high attraction of fruit flies in the genera Bactrocera and Zeugodacus to cuelure traps, with minimal attraction to fruit flies in the genus Dacus (Gnanvossou et al. 2017, Royer et al. 2020). In the first season, there was no effect of the treatments on the daily catches of Z. cucurbitae, while in the second season, Z. cucurbitae catches were higher in farms with no PS than in farms with PS, implying that the presence of stingless bee colonies did not increase the abundance of Z. cucurbitae. Contrary to our study, Toukem et al. (2022) reported an increase in B. dorsalis catches on avocado farms with honeybee colonies, which was attributed to the smell of honey sugar and fermentation of beebread, which attracted fruit flies (Lee et al. 2015). Furthermore, in the second season, Z. cucurbitae catches were highest in control farms across low NDVI, presumably due to fewer natural enemies (Copeland et al. 2006) and as reported by Mwatawala et al. (2006) for Tanzanian lowland regions. Moreover, Toukem et al. (2020) reported a higher abundance of Tephritid fruit flies on avocado in low NDVI. We found the lowest daily catches of Z. cucurbitae in IPPM farms in medium NDVI, but only during the second season, illustrating a synergistic effect. Unlike in our study, Toukem et al. (2022) reported no synergistic effect of IPM and PS in reduction of Tephritid fruit fly densities in avocado production systems, however, independently, IPM reduced fruit fly pest abundance. Findings from our study reported a main effect of NDVI on the catches of D. bivittatus, with the highest catches observed across low NDVI. This confirms the importance of presence of vegetation to reduce cucurbit-infesting fruit fly populations. Although the daily catches of D. bivittatus were generally lower than those of Z. cucurbitae, their population densities were higher in IPPM farms and lower across PS farms. Also, D. punctatifrons populations were higher in farms without colonies during the second season. Presumably, D. bivittatus and D. punctatifrons adults could have been attracted to the smell of honey sugar and fermentation of beebread from the Hypotrigona sp. colonies (Mohammad et al. 2021). Our study reports emergence of Z. cucurbitae, D. bivittatus, D. vertebratus, and D. ciliatus from incubated pumpkin fruits in very low numbers. The low emergence could be attributed to the natural control by Psyttalia sp., a fruit fly parasitoid, which was documented emerging from the incubated infested pumpkin fruits. Although present in trap catches, D. punctatifrons was not among the fruit fly species that emerged from incubated fruits. However, other studies have reported the emergence of D. punctatifrons from cucurbits, in particular watermelon (Layodé et al. 2019). Surprisingly, in the first season, highest infestation was reported in IPPM farms and lowest infestation in control farms, possibly due to high attraction of adult Tephritid fruit flies to the smell of honey and fermenting beebread from the Hypotrigona sp. colonies in IPPM farms, which may have resulted in successful mating and oviposition. In the second season, lowest infestation level was reported in IPM farms in low NDVI, indicating that the IPM package reduced the number of fruit flies. During the first season, the number of fruits per plant was not different across treatments, whereas in the second season, the number of fruits per plant was higher in farm plots with IPM than those without, illustrating the effectiveness of the IPM package to combat Tephritid fruit flies in cucurbits, as illustrated by Githiomi et al. (2019) and Midingoyi et al. (2019) for mango, Pecenka et al. (2021) for watermelon, and Toukem et al. (2022) for avocado. However, a higher number of fruits/plant reported in IPM farms in the second season did not translate to higher yields across the treatments. We, therefore, recommend for the IPM package to be fine-tuned for improved pumpkin yields. Introduction of Hypotrigona sp. hives did not significantly increase the number of pumpkin fruits/plant in both seasons, indicating that Hypotrigona sp. may not be an effective pumpkin flower pollinator. In contrast to this finding, a previous study reported effective pollination of C. sativus using Hypotrigona gribodoi yielding 90% seed germination (Kiatoko et al. (2021b). Hypotrigona sp. has a small body size of 2- to 3-mm body length (Heard, 1999), and the number of pollen attached to the body of stingless bees is positively correlated to body weight (Pangestika et al. 2017). Since a total of 1,500–2,000 pollen grains/flower are required for complete pumpkin pollination (Vidal et al. 2010), sufficient pollination may not have taken place. Additionally, the short period of anthesis (6–12 h) in pumpkin may have limited duration for pollen retrieval by stingless bees (Vidal et al. 2010). Although spinosad-based protein baits used in the study are widely incorporated in IPM programs as a replacement for organophosphates for the control of fruit flies, they have been reported to have toxic effects on Africanized honeybees and stingless bees such as Plebeia lucii Moore (Apidae: Meliponini) (Lopes et al. 2018, Marques et al. 2020). However, a recent study by Toukem et al. (2022) in avocado production systems in Kenya reported contradictory results, with abundant honeybee populations documented in IPM farms incorporating spinosad-based protein baits. Additionally, Koech et al. (2023) reported low spinosad contamination on pollen collected by honeybees in avocado farms across Muranga landscapes. The contamination was below the recommended EU limits and therefore not detrimental to honeybee health. Yield in the first season was 2.6 times higher than yield in the second season, reflecting the importance of sufficient rain during the growing season. In the second season, fruit fly densities were much higher (12.39 ± 0.78 catches/trap/day) than in the first season (7.51 ± 0.45 catches/trap/day). Frequent heavy rains encountered during the first season may have attributed to the reduction in fruit fly densities, as demonstrated by Ye and Liu (2005). A previous study has indicated that fruit fly populations in avocado are affected by landscape and time factors, with population densities generally higher in lower NDVI compared to higher NDVI (Toukem et al. 2020). We reported similar findings, but only for Z. cucurbitae and D. bivittatus populations, mainly in the second season. Presumably, when high pest populations are present, lack of vegetation translates in a reduced presence of natural enemies (Bianchi et al. 2006). This study revealed no synergistic effect between IPM and PS in suppressing Tephritid fruit fly population densities and damage. This was especially evident during the main growing season with low Tephritid fruit fly pest densities. However, in season 2, with high Tephritid fruit fly pest densities, Z. cucurbitae population densities were reduced in IPPM, but this reduction did not translate in reduced fruit damage or increased yield. Hypotrigona sp. may not be an efficient pollinator of pumpkin and we recommend testing other African stingless bees such as M. ferruginea and M. bocandei in pumpkin production systems for better pollination services and improved yields. Acknowledgments This study was carried out with financial support from the German Federal Ministry of Economic Cooperation and Development (BMZ) through icipe to KALRO. The funds were administered through the Deutsche Gesellschaft für Internationale Zusammenarbeit (GIZ) Fund for International Agricultural Research (FIA), grant number 17.7860.4–001; the Norwegian Agency for Development Cooperation, the Section for Research, Innovation, and Higher Education, grant number RAF-3058 KEN-18/0005; the Swedish International Development Cooperation Agency (Sida); the Swiss Agency for Development and Cooperation (SDC); the Australian Centre for International Agricultural Research (ACIAR); the Federal Democratic Republic of Ethiopia; and the Government of the Republic of Kenya. The views expressed herein do not necessarily reflect the official opinion of the donors. Author Contributions Nduta A. Waithaka (Conceptualization-Equal, Formal analysis-Supporting, Investigation-Lead, Methodology-Equal, Writing – original draft-Lead, Writing – review & editing-Equal), Muo Kasina (Conceptualization-Equal, Methodology-Equal, Project administration-Equal, Resources-Equal, Supervision-Equal, Writing – review & editing-Lead), Namikoye Samita (Conceptualization-Equal, Methodology-Equal, Supervision-Equal, Writing – review & editing-Equal), Mary M. Guantai (Conceptualization-Equal, Methodology-Equal, Supervision-Equal, Writing – review & editing-Equal), Evanson R. Omuse (Formal analysis-Equal, Writing – review & editing-Equal), Nadia K. Toukem (Formal analysis-Equal, Writing – review & editing-Equal), H. Michael Lattorff (Conceptualization-Equal, Methodology-Equal, Writing – review & editing-Equal), Elfatih Abdel-Rahman (Investigation-Equal, Methodology-Equal, Writing – review & editing-Equal), Marian Adam (Methodology-Equal, Writing – review & editing-Equal), Thomas Dubois (Conceptualization-Equal, Methodology-Equal, Project administration-Lead, Resources-lead, Writing – review & editing-Lead). Conflict of interest The authors declare no conflict of interest. ==== Refs References Abdel-Rahman  EM, Landmann  T, Kyalo  R, Ong’amo  G, Mwalusepo  S, Sulieman  S, Le Ru  BP.  Predicting stem borer density in maize using RapidEye data and generalized linear models. Int J Appl Earth Obs. 2017:57 :61–74. 10.1016/j.jag.2016.12.008 Adan  M, Abdel-Rahman  EM, Gachoki  S, Muriithi  BW, Lattorff  HMG, Kerubo  V, Landmann  T, Mohamed  SA, Tonnang  HEZ, Dubois  T.  Use of earth observation satellite data to guide the implementation of integrated pest and pollinator management (IPPM) technologies in an avocado production system. Remote Sens Appl Soc Environ. 2021:23 :10566. 10.1016/j.rsase.2021.100566 Akutse  KS, Subramanian  S, Maniania  KN, Dubois  T, Ekesi  S.  Biopesticide research and product development in Africa for sustainable agriculture and food security—experiences from the International Centre of Insect Physiology and Ecology (icipe). Front Sustain Food Syst. 2020:4 :563016. 10.3389/fsufs.2020.563016 Aluja  M, Birke  A, Ceymann  M, Guillen  L, Arrigoni  E, Baumgartner  D, Pasacio- Villafan  C, Samietz  JR.  Agroecosystem resilience to an invasive insect species that could expand its geographical range in response to global climate change. Agric Ecosyst Environ. 2014:186 :54–63. 10.1016/j.agee.2014.01.017. Badii  KB, Billah  MK, Nyarko  G.  Review of pest status, economic impact and management of fruit-infesting flies (Diptera: Tephritidae) in Africa. Afr J Agric Res. 2015:10 :1488–1498. 10.5897/AJAR2014.9278 Bates  D, Maechler  M, Bolker  B, Walker  S.  Fitting linear mixed effects models using lme4. J Stat Softw. 2015:67 :1–48. 10.18637/jss.v067.i01 Beck  PSA, Jonsson  P, Høgda  KA, Karlsen  SR, Eklundh  L, Skidmore  AK.  A ground-validated NDVI dataset for monitoring vegetation dynamics and mapping phenology in Fennoscandia and the Kola peninsula. Int J Rem Sens. 2007:28 :4311–4330. 10.1080/01431160701241936 Bianchi  FJJA, Booij  CJH, Tscharntke  T.  Sustainable pest regulation in agricultural landscapes: a review on landscape composition, biodiversity and natural pest control. Proc R Soc B. 2006:273 :1715–1727. 10.1098/rspb.2006.3530 Biddinger  DJ, Rajotte  EG.  Integrated pest and pollinator management—adding a new dimension to an accepted paradigm. Curr Opin Insect Sci. 2015:10 :204–209. 10.1016/j.cois.2015.05.012 29588010 Carvalheiro  LG, Seymour  CL, Veldtman  R, Nicolson  SW.  Pollination services decline with distance from natural habitat even in biodiversity rich areas. J Appl Ecol. 2010:47 :810–820. 10.1111/j.1365-2664.2010.01829.x Charlery de la Masselière  M, Ravingé  V, Facon  B, Lefeuvre  P, Massol  F, Quilici  S, Duyck  PF.  Changes in phytophagous insect host ranges following the invasion of their community: long- term data for fruit flies. Ecol Evol. 2017:7 :5181–5190. 10.1002/ece3.2968 28770058 Cloyd  AR, Bethke  JA.  Impact of neonicotinoid insecticides on natural enemies in greenhouses and interior spaces environment. Pest Manag Sci. 2015:67 :3–9. 10.1002/ps.2015 Copeland  RS, Wharton  RA, Luke  Q, De Meyer  M, Lux  S, Zenx  N, Machera  P, Okumu  M.  Geographic distribution, host fruit and parasitoids of African fruit fly pests Ceratitis anonae, Ceratitis cosyra, Ceratitis fasciventris and Ceratitis rosa (Diptera: Tephritidae) in Kenya. Ann Entomol Soc Am. 2006:99 :261–278. 10.1603/0013-8746(2006)099[0261:GDHFAP]2.0.CO;2 De Meyer  M, Quilici  S, Franck  A, Chadhouliati  AC, Issimaila  MA, Youssoufa  MA, Abdoul-Karime  AL, Barbet  A, Attié  M, White  IM.  Records of frugivorous fruit flies (Diptera: Tephritidae: Dacini) from the Comoro archipelago. Afr Invertebr. 2012:53 :69–77. 10.5733/afin.053.0104 Deguine  JP, Atiama-Nurbel  T, Douraguia  E, Chiroleu  F, Quilici  S.  Species diversity within a community of the cucurbit fruit flies Bactrocera cucurbitae, Dacus ciliatus, and Dacus demmerezi roosting in corn borders near cucurbit production areas of Reunion Island. J Insect Sci. 2012:12 :32. 10.1673/031.012.3201 22958393 Dhillion  M, Singh  R, Naresh  J, Sharma  H.  The melon fruit fly, Bactrocera cucurbitae: a review of its biology and management. J Insect Sci. 2005:5 :40. 10.1093/jis/5.1.40 17119622 Doorenweerd  C, Leblanc  L, Norrbom  AL, Jose  MS, Rubinoff  D.  A global checklist of the 932 fruit fly species in the tribe Dacini (Diptera: Tephritidae). Zookeys. 2018:730 :17–54. 10.3897/zookeys.730.21786 Fletcher  BS.  Temperature-development rate relationships of the immature stages and adults of tephritid fruit flies, In: Robinson  AS, Hooper  G, editors. Fruit flies: their biology, natural enemies and control, vol. 3B . Amsterdam (The Netherlands): Elsevier; 1989. p. 273–289. Githiomi  C, Muriithi  B, Irungu  P, Mwungu  CM, Diiro  G, Affognon  H, Mburu  J, Ekesi  S.  Economic analysis of spillover effects of an integrated pest management (IPM) strategy for suppression of mango fruit fly in Kenya. Food Pol. 2019:84 :121–132. 10.1016/J.FOODPOL.2019.03.006 Gnanvossou  D, Hanna  R, Goergen  G, Salifu  D, Tanga  MC, Mohamed  SA, Ekesi  S.  Diversity and seasonal abundance of tephritid fruit flies in three agroecosystems in Benin, West Africa. J Appl Entomol. 2017:00 :1–12. 10.1111/jen.12429 Gorelick  N, Hancher  M, Dixon  M, Ilyushchenko  S, Thau  D, Moore  R.  Google earth engine: planetary-scale geospatial analysis for everyone. Remote Sens Environ. 2017:202 :18–27. 10.1016/j.rse.2017.06.031 Heard  TA.  The role of stingless bees in crop pollination. Annu Rev Entomol. 1999:44 :183–206. 10.1146/annurev.ento.44.1.183 15012371 Hill  MP, Bertelsmeier  C, Clusella-Trullas  S, Garnas  J, Robertson  MP, Terblanche  JS.  Predicted decrease in global climate suitability masks regional complexity of invasive fruit fly species response to climate change. Biol Invasions. 2016:18 :1105–1119. 10.1007/s10530-016-1078-5 (HCD) Horticultural Crops Directorate. 2018. Horticultural annual validated report for 2017–2018: agriculture and food authority (AFA), government of Kenya. Nairobi (Kenya): HCD. Indrianingsih, AW, Rosyida  VT, Apriyana  W, Hayati  SN, Nisa  K, Darsh  C, Kusumaningrum  A, Ratih  D, Indirayati  N.  Comparisons of antioxidant activities of two varieties of pumpkin Cucurbita moschata and Cucurbita maxima extracts.  IOP Conf Ser: Earth Environ Sci.  2019:251 :012021. 10.1088/1755-1315/251/1/012021 Ito  Y, Kakinohana  H, Yamagishi  M, Kohama  T.  Eradication of the melon fly, Bactrocera cucurbitae, from Okinawa, Japan, by means of the sterile insect technique, with special emphasis on the role of basic studies. J Asia-Pac Entomol. 2003:6 :119–129. 10.1016/S1226-8615(08)60177-6 Kambura, C, D.  Kilalo, J.  Muthomi, S.  Ekesi. 2016. Diversity and host preference of tephritid fruit fly species infesting cucurbit and major horticultural crops grown in the lower coastal, Kenya [MS thesis]. [Nairobi (Kenya)]: University of Nairobi. Kambura  C, Tanga  CM, Kilalo  D, Muthomi  J, Salifu  D, Rwomushana  I, Mohamed  SA, Ekesi  S.  Composition, host range and host suitability of vegetable infesting tephritids on cucurbits cultivated in Kenya. Afr Entomol. 2018:26 :379–397. 10.4001/003.026.0379 Kiatoko, N, Mueke  JM, Raina  KS, Muli  E.  2012. Distribution, behavioural biology, rearing and pollination efficiencies of five stingless bee species (Apidae: Meliponinae) in Kakamega forest, Kenya [PhD thesis]. [Nairobi (Kenya)]: Kenyatta University. Kiatoko  N, Pozo  MI, Van Oystaeyen  A, Langevelde  FV, Wäckers  F, Kumar  RS, Hundt  B, Jaramillo  J.  Effective pollination of greenhouse galia musk melon (Cucumis melo L. var. reticulatus ser.) by afrotropical stingless bee species. J Apic Res. 2021a:61 :664–674. 10.1080/00218839.2021.2021641 Kiatoko  N, Pozo  MI, Van Oystaeyen  A, Musonye  M, Kika  J, Wäckers  F, Langevelde  FV, Hundt  B, Jaramillo  J.  African endemic stingless bees as an efficient alternative pollinator to honey bees in greenhouse cucumber Cucumis sativus L. J Apic Res. 2021b:61 :1–13. 10.1080/00218839.2021.2013421 Koech  SJ, Karanja  RHN, Kurgat  JK, Mokaya  HO, Dubois  T, Lattorff  HMG.  Pesticide contamination and their botanical sources in pollen loads collected by honeybees in Kenya: a spatio-temporal context. Agric Ecosyst Environ. 2023:343 :108264. 10.1016/j.agee.2022.108264 Korir  JK, Affognon  HD, Ritho  CN, Kingori  WS, Irungu  P, Mohamed  SA, Ekesi  S.  Grower adoption of an integrated pest management package for management of mango-infesting fruit flies (Diptera: Tephritidae) in Embu, Kenya. Int J Trop Insect Sci. 2015:35 :80–89. 10.1017/S1742758415000077 Layodé  BFR, Onzo  A, Karlsson  MF.  Watermelon-infesting tephritidae fruit fly guild and parasitism by Psyttalia phaeostigma (Hymenoptera: Braconidae). Int J Trop Insect Sci. 2019:40 (1 ):157–166. 10.1007/s42690-019-00066-x Leach  A, Pecenka  J, Kaplan  I.  Does IPPM bear fruits? Evaluating reduced-risk insecticide programmes on pests, pollinators and marketable yield. J Appl Ecol. 2022:59 :2993–3002. 10.1111/1365-2664.14294 Lee  FJ, Rusch  DB, Stewart  FJ, Mattila  HR, Newton  IL.  Saccharide breakdown and fermentation by the honeybee gut microbiome. Environ Microbiol. 2015:17 (3 ):796–815. 10.1111/1462-2920.12526 24905222 Lopes  MP, Fernandes  KM, Tomé  HVV, Gonçalves  WG, Miranda  FR, Serrão  JE, Martins  GF.  Spinosad-mediated effects on the walking ability, midgut and malpigian tubules of Africanized honeybee workers. Pest Manag Sci. 2018:74 (6 ):1311–1318. 10.1002/ps.4815 29194936 Machakos CIDP. 2019. Machakos County Integrated Development Plan II (2018–2022). Machakos (Kenya):  County Government of Machakos. Marques  RD, Lima  MAP, Marques  RD, Bernardes  RC.  A spinosad-based formulation reduces the survival and alters the behavior of the stingless bee Plebeia lucii. Neutrop Entomol. 2020:49 :578–585. 10.1007/s13744-020-00766-x Mebdoua, S.  2018. Pesticide residues in fruits and vegetables. In: Mérillon  JM, Ramawat  K, editors. Bioactive molecules in food. Reference series in phytochemistry. Cham (Switzerland): Springer; 2018. p. 1715–1753. Merle  I, Hipólito  J, Requier  F.  Towards integrated pest and pollinator management in tropical crops. Curr Opin Insect Sci. 2021:50 :100866. 10.1016/j.cois.2021.12.006 34971783 Messelink  GJ, Van Maanen  R, Van Steenpaal  SEF, Janssen  A.  Biological control of thrips and whiteflies by a shared predator: two pests are better than one. Biol Control. 2008:44 :372–379. 10.1016/j.biocontrol.2007.10.017 Midingoyi  SKG, Kassie  M, Muriithi  B, Diiro  G, Ekesi  S.  Do farmers and the environment benefit from adopting integrated pest management practices? Evidences from Kenya. J Agric Econ. 2019:70 :452–470. 10.1111/1477-9552.12306 Mohammad  SM, Ab-Rashid  NKM, Zawawi  N.  Stingless bee-collected pollen (beebread): chemical and microbiology properties and health benefits. Molecules. 2021:26 :957. 10.3390/molecules26040957 33670262 Munyuli  TMB.  2010. Pollinator biodiversity and economics of pollination services in Uganda [PhD dissertation]. [Kampala (Uganda)]: Makerere University. Muriithi  BW, Affognon  HD, Diiro  GM, Kingori  SW, Tanga  CM, Nderitu  PW, Mohamed  SA, Ekesi  S.  Impact assessment of integrated pest management (IPM) strategy for suppression of mango-infesting fruit flies in Kenya. Crop Prot. 2016:81 :20–29. 10.1016/J.CROPRO.2015.11.014 Muriithi  BW, Gathogo  NG, Diiro  GM, Mohamed  SA, Ekesi  S.  Potential adoption of integrated pest management strategy for suppression of mango fruit flies in East Africa: an ex ante and ex post analysis in Ethiopia and Kenya. Agriculture. 2020:10 (7 ):278. 10.3390/agriculture10070278 Mwatawala  MW, De Meyer  M, Makundi  R, Maerere  AP.  Seasonality and host utilization of the invasive fruit fly, Bactrocera invadens (Diptera: Tephritidae) in central Tanzania. J Appl Entomol. 2006:130 :530–537. 10.1111/j.1439-0418.2006.01099.x Nawirska  A, Figiel  A, Kucharska  AZ, Sokół-Łętowska  A, Biesiada  A.  Drying kinetics and quality parameters of pumpkin slices dehydrated using different methods. J Food Eng. 2009:94 (1 ):14–20. 10.1016/J.JFOODENG.2009.02.025 Ndengwa, R, Bett  E, Muchungi  A, Nguthi  NF.  2016. Socio-economic factors influencing smallholder pumpkin production, consumption and marketing in Eastern and Central Kenya regions [MS thesis]. [Nairobi (Kenya)]: Kenyatta University. Ndungu  NN, Kiatoko  N, Sole  CL, Pirk  CWW, Abdullahi  AY, Raina  SK, Masiga  DK.  Resolving taxonomic ambiguity and cryptic speciation of Hypotrigona species through morphometrics and DNA barcoding. J Apic Res. 2017:57 :354–363. 10.1080/00218839.2018.1426348 Niassy  S, Murithii  B, Omuse  ER, Kimathi  E, Tonnang  H, Ndlela  S, Mohamed  S, Ekesi  S.  Insight on fruit fly IPM technology uptake and barriers to scaling in Africa. Sustainability. 2022:14 (5 ):2954. 10.3390/su14052954 Njoroge  GN, Gemmil  B, Bussmann  R, Newton  LE, Ngumi  VM.  Diversity and efficiency of wild pollinators of watermelon Citrullus lanatus (Thunb.) Mansf.) at Yatta Kenya. J Appl Hortic. 2010:12 :35–41. 10.37855/jah.2010.v12i01.07 Nobre  T, Gomes  L, Rei  FT.  A Re-evaluation of olive fruit fly organophosphate-resistant ace alleles in Iberia, and field-testing population effects after in-practice dimethoate use. Insects. 2019:10 (8 ):232. 10.3390/insects10080232 31374903 O’Rourke  ME, Jones  LE.  Analysis of landscape-scale insect pest dynamics and pesticide use: an empirical and modelling study. Ecol Appl. 2011:21 (8 ):3199–3210. 10.1890/10-1180.1 Okolle  JN, Ntoifor  NN.  Field ovipositional behavior and laboratory studies on development of Dacus punctatifrons (Diptera: Tephritidae) on tomato. Insect Sci. 2005:12 :393–398. 10.1111/j.1005-295X.2005.00049.x Onsongo  SK, Gichimu  BM, Akutse  KS, Dubois  T, Mohamed  SA.  Performance of three isolates of Metarhizium anisopliae and their virulence against Zeugodacus cucurbitae under different temperature regimes, with global extrapolation of their efficiency. Insects. 2019:10 (9 ):270. 10.3390/insects10090270 31454931 Pangestika  NW, Atmowidi  T, Kahono  S.  Pollen load for flower constancy of three species of stingless bees (Hymenoptera: Apidae: Meliponinae). Trop Life Sci Res. 2017:28 (2 ):179–187. 10.21315/tlsr2017.28.2.13 28890769 Pecenka  RJ, Ingwell  LL, Foster  ER, Krupke  HC, Kaplan  I.  IPM reduces insecticide applications by 95% while maintaining or enhancing crop yields through wild pollinator conservation. Proc Natl Acad Sci USA. 2021:118 :1–10. 10.1073/pnas.2108429118 R Core Team. R: a language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing; 2020. Rődlach  A.  Aids are in the food: Zimbabwean’s association between nutrition and HIV/AIDS and their potential for addressing food insecurity and HIV/AIDS. Ann Anthropol Pract. 2011:35 :219–237. 10.1111/J.2153-9588.2011.01076.X Royer  JE, Tan  KH, Mayer  DG.  Comparative trap catches of male Bactrocera, Dacus and Zeugodacus fruit flies (Diptera: Tephritidae) with four floral phenylbutanoid lures (anisyl acetone, cuelure, raspberry ketone, and zingerone) in Queensland, Australia. Environ Entomol. 2020:49 (4 ):815–822. 10.1093/ee/nvaa056 32514581 Sagwe  RN, Peters  MK, Dubois  T, Steffan-Dewenter  I, Lattorff  HMG.  Pollination supplementation mitigates pollination deficits in smallholder avocado (Persea americana Mill.) production systems in Kenya. Basic Appl Ecol. 2021:56 :392–400. 10.1016/j.baae.2021.08.013 Shackelford  G, Steward  PR, Benton  TG, Kunin  WE, Potts  SG, Biesmeijer  JC, Sait  SM.  Comparison of pollinators and natural enemies: a meta-analysis of landscape and local effects of abundance and richness in crops. Biol Rev. 2013:88 (4 ):1002–1021. 10.1111/brv.12040 23578337 Shikano  I, Gutierrez-Coarite  R, Streit  C, Perez  E, Fujitani  E, Mau  RFL.  Field tests of three alternative insecticides with protein bait for the development of an insecticide rotation program to control melon flies, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae). Insects. 2022:13 (7 ):629. 10.3390/insects13070629 35886806 Spafford  H, Chou  MY, Mau  RFL, Vargas  RI.  Suppression of female melon fly, Zeugodacus cucurbitae, with cue‐lure and fipronil bait stations through horizontal insecticide transfer. Entomol Exp Appl. 2018:166 :94–101. 10.1111/eea.12638 Toukem, NK, Yusuf  AA, Dubois  T, Abdel-Rahman  EM, Adan  MS, Mohamed  SA.  Landscape vegetation productivity influences population dynamics of key pests in small avocado farms in Kenya. Insects. 2020:11 :424. 10.3390/insects11070424 32659943 Toukem  NK, Mohamed  SA, Yusuf  AA, Lattorff  HMG, Copeland  RS, Dubois  T.  Interactions between integrated pest management, pollinator introduction, and landscape context on avocado Persea americana productivity. Entomol Gen. 2022:42 :579576–579587. 10.1127/entomologia/2022/1365 Vargas  RI, Piñero  JC, Leblanc  L.  An overview of pest species of Bactrocera fruit flies (Diptera: Tephritidae) and the integration of biopesticides with other biological approaches for their management with a focus on the Pacific region. Insects. 2015:6 (2 ):297–318. 10.3390/insects6020297 26463186 Vidal  MDG, Jong  DD, Wien  HC, Morse  RA.  Pollination and fruit set in pumpkin (Cucurbita pepo) by honeybees. Braz J Bot. 2010:33 :107–113. 10.1590/S0100-84042010000100010 Wang  S, Zhang  L, Huang  C, Qiao  N.  An NDVI-based vegetation phenology is improved to be more consistent with photosynthesis dynamics through applying a light use efficiency model over boreal high-latitude forests. Rem Sens. 2017:9 :695. 10.3390/rs9070695 White  IM.  Taxonomy of the Dacina (Diptera: Tephritidae) of Africa and the Middle East. Afr Entomol Mem. 2006:2 :1–156. Ye  H, Liu  J.  Population dynamics of Bactrocera dorsalis (Diptera: Tephritidae) in Xishuangbanna of Southern Yunnan. J Appl Ecol. 2005:16 :1330–1334.