An Innate Preference of Bumblebees for Volatile Organic Compounds Emitted by Phaseolus vulgaris Plants Infected With Three Different Viruses

Cucumber mosaic virus (CMV)-infected tomato (Solanum lycopersicum L.) plants emit volatile organic compounds (VOCs) attractive to bumblebees (Bombus terrestris L.), which are important tomato pollinators, but which do not transmit CMV. We investigated if this effect was unique to the tomato-CMV pathosystem. In two bean (Phaseolus vulgaris L.) cultivars, infection with the potyviruses bean common mosaic virus (BCMV) or bean common mosaic necrosis virus (BCMNV), or with the cucumovirus CMV induced quantitative changes in VOC emission detectable by coupled gas chromatography–mass spectrometry. In free-choice olfactometry assays bumblebees showed an innate preference for VOC blends emitted by virus-infected non-flowering bean plants and flowering CMV-infected bean plants, over VOCs emitted by non-infected plants. Bumblebees also preferred VOCs of flowering BCMV-infected plants of the Wairimu cultivar over non-infected plants, but the preference was not significant for BCMV-infected plants of the Dubbele witte cultivar. Bumblebees did not show a significant preference for VOCs from BCMNV-infected flowering bean plants but differential conditioning olfactometric assays showed that bumblebees do perceive differences between VOC blends emitted by flowering BCMNV-infected plants over non-infected plants. These results are consistent with the concept that increased pollinator attraction may be a virus-to-host payback, and show that virus-induced changes in bee-attracting VOC emission is not unique to one virus-host combination.


INTRODUCTION
Plants emit volatile organic compounds (VOCs) that function as signals in plant-plant communication, defense, and attraction of pollinators, seed dispersers, and other beneficial organisms (Dudareva and Pichersky, 2008). Insects, and in particular bees, are the most important group of plant pollinators (Knight et al., 2005;Ghazoul, 2006;Klein et al., 2007). With respect to insect-pollinated plants, pollinator visitation (or artificial pollination) is required to obtain maximal seed and fruit production (Klein et al., 2007;Klatt et al., 2014). Bees possess over 130 distinct antennal olfactory receptors and exploit olfactory cues when foraging and pollinating flowers (Chittka and Raine, 2006). The olfactory cues are important in differentiating between plant species, and in signaling the appropriate phenological stage for pollinator visits (Chen et al., 2009;Hossaert-McKey et al., 2010;Song et al., 2014). Indeed, bees exhibit innate preferences for certain floral odors (Raguso, 2008), and VOCs have been shown to act as long-distance flower location signals for pollinators (Kunze and Gumbert, 2001;Burger et al., 2010;Suchet et al., 2011) and can trigger an immediate pollinator decision to land on a flower (Lunau, 1992). However, the complex effects of non-floral VOCs or floral/non-floral VOC combinations on bee foraging and other insect behaviors are still not well understood (Junker et al., 2017), and in particular the influence of microorganisms (including viruses) on plant VOC emission is relatively under-investigated (Beck and Vanette, 2017).
Phaseolus vulgaris L. (common bean) is a plant that benefits from insect pollination with increased seed production (Ibarra-Perez et al., 1999;Kingha et al., 2012;Ramos et al., 2018;Elisante et al., 2020) and is also host to a number of viruses that can reduce seed yield. Bean common mosaic virus (BCMV) and bean common mosaic necrosis virus (BCMNV) are closely related potyviruses that are among the most prevalent viruses infecting legumes and both result in serious crop loss (Morales, 2006). Cucumber mosaic virus (CMV) is a cucumovirus with a very wide host range (Yoon et al., 2019) that causes serious epidemics and losses in cultivated common bean (Morales, 2006;Thompson et al., 2015;Worrall et al., 2015). Aphids transmit CMV, BCMV, and BCMNV, particles of which are carried between hosts on the stylets (probing mouthparts) of these insects, but the viruses are also transmitted efficiently through seed of infected bean plants (Morales, 2006;Worrall et al., 2015;Yoon et al., 2019). Viruses induce changes in plant metabolism that influence insect vector behavior (Carr et al., 2018;Mauck et al., 2016). For example, CMV, BCMV, and BCMNV induce biochemical changes in common bean that alter the feeding behavior of aphid vectors (Wamonje et al., 2020), and CMV induces alterations in the emission of VOCs that attract aphid vectors to infected hosts (Mauck et al., 2010). Such virus-induced changes in VOC emission by infected hosts are thought to contribute to an increase the rate of aphid-mediated virus transmission (Mauck et al., 2016;Carr et al., 2018;Donnelly et al., 2019). Currently, there is a paucity of data about virus-induced changes on other community members, such as pollinators.
There is no evidence that CMV, BCMV, or BCMNV are transmitted by bees, or by other pollinators. Indeed, comparatively few plant viruses are vectored by pollinators (Hull, 2014). It was therefore surprising to discover that bumblebees (Bombus terrestris) could perceive differences in the VOC blends emitted by uninfected and CMV-infected plants of tomato (Solanum lycopersicum), and Arabidopsis thaliana (Groen et al., 2016). In the same study it was also found that bumblebees enhanced seed production by CMV-infected tomato plants to a greater extent than for mock-inoculated plants (Groen et al., 2016). Mathematical modeling indicated that if a pollinator bias in favor of virus-infected plants in wild host populations occurred, it could outweigh the benefits of resistance, and natural selection would instead favor inheritance by plants of alleles for virus susceptibility, rather than those conditioning resistance (Groen et al., 2016). In this study we investigated whether the innate preference of bumblebees for VOC blends emitted by virus-infected plants is peculiar only to the previously investigated tomato-CMV pathosystem, or if this can occur in other virus-host combinations, such as bean plants infected with CMV, BCMV, or BCMNV.

Plant-Emitted Volatile Organic Compound Collection and Analysis
Headspace VOCs from flowering plants of the Wairimu variety that had been mock-inoculated or infected with BCMV, BCMNV, or CMV, were analyzed by GC-MS. VOCs were collected from flowering plants between 1 and 2 weeks after the first flowers opened (22-26 days post-inoculation) by dynamic headspace trapping and analyzed by gas chromatography-mass spectrometry (GC-MS) as previously described (Beale et al., 2006;Webster et al., 2008;Groen et al., 2016). Briefly, headspace volatiles were collected from plants over a period of 24 h onto Porapak Q filters [50 mg, 60/80 mesh size, Supelco (Sigma-Aldrich)]. The plants were contained in a 1.0 liter bell jar clamped to two semi-circular metal plates with a hole in the center to accommodate the stem. Charcoal-filtered air was pumped in at the bottom of the container at a rate of 800 ml.min −1 and drawn out through the Porapak Q filter at the top, at a rate of 750 ml.min −1 . Leaf fresh weight and dry weight were measured to enable normalization of the volatile abundance. Trapped organic chemicals were eluted from the Porapak Q filter with diethyl ether for analysis by gas chromatography coupled to mass spectrometry (GC-MS). Volatiles were separated on a capillary GC column (ZB-5MSi, Thermo Scientific, United Kingdom). The injection volume (splitless) was 1 µl, the injector temperature was 200 • C, and helium was used as the carrier gas at a constant flow rate of 2.6 ml min −1 in an oven maintained at 30 • C for 5 min and then programmed at 15 • C.min −1 to 230 • C. The column was directly coupled to a mass spectrometer (ISQ LT, Thermo Scientific, United Kingdom) with a MS transfer line temperature of 240 • C. Ionization was by electron impact with an ion source temperature of 250 • C in positive ionization. Mass ions were detected between 30 and 650 m/z. Data were collected using Xcalibur software (Thermo Scientific). VOCs were identified by comparing spectra with those curated by the National Institute of Standards and Technology 1 and by comparison with an authentic standard. Abundant VOCs were quantified using standard curves generated with pure benzaldehyde, linalool, ocimene, and pinene, and others quantified as pinene equivalents.
To determine if virus infection causes overall changes in the VOC blends emitted by plants, principal component analysis was performed with MetaboAnalyst 5.0 (Chong et al., 2019) on the quantified VOCs after Pareto scaling. Multivariate analysis of variance (MANOVA), using the Wilks' lambda test statistic, was carried out with the Real Statistics Resource Pack software, Release 7.6 (Zaiontz, 2020). Significant results from MANOVAs were followed up with univariate one-way ANOVAs to assess which VOC(s) differed between groups using R (RStudio Team, 2020). Significant univariate ANOVAs were then followed up with Tukey's HSD post hoc test to allow detailed comparisons between individual pairs of groups.

Bumblebee Free-Choice and Differential Conditioning Olfactometry Assays
Buff-tailed bumblebees (Bombus terrestris audax L.) were sourced from Koppert Biological Systems (Berkel en Rodenrijs, Netherlands), or Syngenta-Bioline (Leicester, United Kingdom). Free-choice preference assays were performed to identify any innate preference of bumblebees for VOCs emitted by virusinfected versus mock-inoculated bean plants. Assays were carried out using flight arenas (Whitney et al., 2009) and procedures previously described with 8 cm diameter/15 cm height cylindrical towers enclosing infected or mock-inoculated plants (Groen et al., 2016). To obscure plants, but allow VOC diffusion, towers were roofed with muslin overlaid with plastic mesh, supporting cups containing 30% (w/v) sucrose to incentivize visitation (Supplementary Figure 1). Between experiments, arenas and tower roofs were cleaned with 30% (v/v) ethanol to remove scent marks. Five mock-inoculated and five virus-infected bean plants were concealed under towers randomly placed in the arena. Thirty paint-marked bees were released one at a time (Supplementary Figure 1). The first 10 choices of each bee were recorded either as acceptance or rejection of a tower. A similar choice test was carried out to assess bee preference for volatiles from flowering versus non-flowering plants, except the first 10 choices of 40 bees was recorded.
Free choice assays were analyzed using binomial generalized linear mixed models, fitted using mixed_model in the GLMMadaptive package (Rizopoulos, 2021). The identity of each bee was used as a random effect in fitting the models, allowing differences between bees to be accounted for. Wald tests assessing whether the fixed effect (i.e., the intercept) in the fitted model differed significantly from 0 were used to determine whether bumblebee preference significantly differed from an expected frequency distribution of 50:50 for either mock-inoculated or virus-infected bean plants. To account for multiple comparisons, p-values were subjected to the Holm-Bonferroni correction.
Differential conditioning ("learning curve") assays determined bumblebees' ability to discriminate between VOCs emitted by infected and mock-inoculated plants, regardless of instinctive preference (Dyer and Chittka, 2004;Whitney et al., 2009;Groen et al., 2016). Cups on towers concealing five flowering mockinoculated plants and five flowering BCMNV-infected plants contained 30% (w/v) sucrose "reward" or 0.12% (w/v) quinine hemisulphate solution "punishment, " respectively, and the first 100 choices of each bee were recorded (Groen et al., 2016). Bumblebees cannot distinguish quinine from sucrose except by taste (Whitney et al., 2009). Statistical analysis of data from differential conditioning assays were analyzed using binomial mixed-effects logistic regression (Foster et al., 2014;Groen et al., 2016), again using mixed_model in the GLMMadaptive package (Rizopoulos, 2021). The identity of individual bees controlled a random intercept, allowing differences between bees to be accounted for in the analysis. Whether or not bees "learnt" was tested using a likelihood ratio test comparing the fit of the model in which the probability of selecting the sucrose reward depended on how many choices had been made by an individual bee with a simpler model in which the probability was independent of the choice (but depended on the bee in question).

Viruses Induce Quantitative Changes in Volatile Organic Compounds Emitted by Flowering Bean Plants
Virus infection decreased the quantity of VOCs emitted by approximately half compared with emission levels of uninfected plants ( Figure 1A). This might be expected since all three viruses cause stunting (Supplementary Figure 2) and the VOC emission rate at the tissue-mass level for plants infected with BCMV, BCMNV, and CMV was not significantly different from mockinoculated plants ( Figure 1B). Principle component analysis on the emitted VOCs showed that the scores from mock-inoculated plants form a discrete cluster while the scores from the virusinfected plants group together (Figure 2A). The scores and loadings biplot explains which variables have the largest effect on each sample ( Figure 2B). (E)-Ocimene, linalool and 4-hexen-1-ol had a strong influence on PC1 (that is mainly responsible for separating the mock and virus-inoculated plant scores), while benzaldehyde had a stronger influence on PC2.
Benzaldehyde, (E)-ocimene, 4-hexen-1-ol acetate, linalool and D-limonene dominated the VOC blend from the mockinoculated bean plants ( Figure 2C). These VOCs, with the exception of 4-hexen-1-ol acetate, are known to be bee attractants (Granero et al., 2005;Knudsen et al., 2006;Dötterl and Vereecken, 2010;Klatt et al., 2013;Krug et al., 2018). Despite stunting, virus-infected plant benzaldehyde-output was maintained at a level equivalent to uninfected plants ( Figure 2C) but (E)-ocimene and 4-hexen-1-ol acetate emissions were, in contrast, disproportionately lower than their decreased tissue mass could account for Figures 2C,D. The lower emission rates of (E)-ocimene and 4-hexen-1-ol acetate were the main cause for the overall decrease in total VOC emission from virusinfected plants (Figure 1) and supports their importance in PCA (Figures 2A,B). BCMNV-infected plants were the most severely stunted (Supplementary Figure 2) and on a tissue per dry weight basis produced greater levels of several VOCs compared to other treatments ( Figure 2D: linalool, dodecane, α-pinene, tetradecane, and tridecane). In most cases this resulted in BCMNV-infected whole plants matching or exceeding the VOC output of mockinoculated plants. Having established that virus infection altered the VOC emissions from bean plants, we went on to test whether this affected bumblebee behavior.

Bumblebees Possess an Innate Preference for Volatile Organic Compounds Emitted by Virus-Infected Bean Plants
In free-choice olfactometry assays bumblebees displayed an innate preference for VOC blends emitted by CMV-infected and BCMV-infected flowering plants of the Wairimu cultivar, as compared to the blend emitted by mock-inoculated plants ( Figure 3A). The preferences for VOCs produced by plants infected with CMV or BCMV were statistically significant. However, bees showed no preference for VOC blends emitted by flowering Wairimu plants that had been infected with BCMNV. In another cultivar of bean, Dubbele witte, bees showed a preference for VOC blends emitted by CMV-infected flowering plants but no preference for those emitted by BCMV-or BCMNV-infected plants. When bumblebees were presented with a choice between non-infected and virus-infected non-flowering plants, they preferred the odors from all virus-infected plants over healthy counterparts. This was the case for both cultivars ( Figure 3B). Because the bee preference for VOCs released from virus-infected plants was more apparent when plants were not flowering, we also presented bumblebees with feeding towers concealing five non-flowering and five flowering non-infected plants ( Figure 3C). Bees did not show a preference between odors from flowering plants over those from non-flowering plants.

Bumblebees Perceive Bean Common Mosaic Necrosis Virus-Induced Changes in Volatile Organic Compound Emission by Flowering Plants but Do Not Find Them Intrinsically Attractive or Repellent
Bumblebees did not exhibit a preference for VOCs emitted by BCMNV-infected flowering plants over those emitted by mockinoculated flowering plants ( Figure 3B). This prompted us to investigate whether this was due to an inability of bumblebees to perceive any differences between the VOC blends emitted by BCMNV-infected plants versus non-infected plants, once they begin flowering. We used differential conditioning assays to determine if bumblebees can distinguish between stimuli (in this case VOC blends) by monitoring the insects' ability to learn to associate one VOC blend with a sucrose reward and the other stimulus with a quinine "punishment." Bumblebees learned rapidly to associate the reward and punishment with VOCs emitted by mock-inoculated and BCMNV-infected plants, respectively, and did so for flowering plants of both bean cultivars (Figure 4). This indicates that bumblebees can distinguish between VOC blends of BCMNV-infected and mock-inoculated flowering plants (Figure 4), despite not demonstrating a significant preference for either VOC blend in free-choice assays ( Figure 3B).  or 40 (C) naïve bumblebees, and relates to the Wald test of whether the intercept in the fitted binomial generalized linear mixed model was zero, i.e., whether there is any evidence -after controlling for systematic differences between individuals -that bees did anything other than choose at random, with all p-values subjected to a Holm-Bonferroni correction to account for multiple testing.

DISCUSSION
We found that naïve bumblebees have an innate preference for VOC blends emitted by common bean plants infected with three different viruses over non-infected plants. It was previously established that CMV-infected tomato plants are more attractive to bumblebees than uninfected plants, and benefit more from bumblebee-mediated pollination (Groen et al., 2016). Since FIGURE 4 | Bumblebees can distinguish between VOC blends emitted by flowering non-infected and BCMNV-infected bean plants despite having no innate preference for either. Using differential conditioning with a sucrose reward and quinine punishment bumblebees learned to distinguish between VOC blends emitted by BCMNV-infected and mock-inoculated flowering common bean. The learning curve indicates the overall ability to distinguish between plant-emitted volatiles analyzed after 100 choices per bee using 10 bees. In flowering cv. Wairimu (A), bumblebees initially chose the sucrose reward at a mean accuracy of 50% in their first 10 choices increasing to about 85% by choices 90 to 100 meaning bumblebees learnt to identify sucrose rewards based on the association with volatiles from flowering mock-inoculated plants. Similar results were obtained with flowering cv. Dubbele witte (B), where bumblebees' first 10 choices had a mean accuracy of 60% for sucrose increasing to a mean of c. 75% by choices 90 to 100. Data are shown pooled over all bees (n = 10) into successive groups of 10 choices, with error bars showing 95% binomial confidence intervals for the proportion of correct choices. The white curve shows the fitted binomial generalized linear mixed model, with blue borders showing 95% confidence intervals on the fitted response. The χ 2 statistic and significance level for the likelihood ratio test assessing whether bees can learn are given at the bottom left of each panel.
CMV is not transmitted by bumblebees, it was hypothesized that this could be a "payback" or compensation to susceptible hosts (Groen et al., 2016). Mathematical modeling further suggested that in wild populations of plants that were dependent upon pollinators for reproductive success, increased pollinator visitation to virus-infected plants might outweigh the benefits of resistance and favor the persistence of alleles for susceptibility in the population (Groen et al., 2016). Our current work shows that virus-induced attraction of bumblebees is not unique or peculiar to the interaction of CMV and tomato, but extends to interactions of CMV and two other viruses (neither of which are pollinatorvectored), the potyviruses BCMV and BCMNV, with common bean. These observations in two different plant species suggest that modification by viruses of the emission of bee-perceivable plant VOCs may have arisen more than once. Common bean and tomato have very different bumblebee pollination mechanisms. Tomato and several other solanaceous plants exhibit "buzz" or "pepper pot" pollination, in which foraging bumblebees grasp the fused anther cone and vibrate it to release pollen. Excess pollen provides bumblebees with a proteinrich reward but no nectar reward is provided (Glover et al., 2004). In contrast, the papilionate floral morphology of beans and many other legumes requires, almost exclusively, Hymenopteran pollinators with sufficient weight and strength to mechanically "trip" the hull and wing petals to expose the stamens and pistil. This enables bees to carry out pollination, and access nectar (a rich sucrose source) as the predominant reward (Leppik, 1966;Córdoba and Cocucci, 2011;Bailes et al., 2018). CMV is a cucumovirus with a broad host range encompassing many plant families (Yoon et al., 2019). Contrastingly, BCMV and BCMNV are legume pathogens that belong to one of the largest family of plant-infecting viruses, the Potyviridae (Ivanov et al., 2014), and are unrelated to CMV, which lies within the family Bromoviridae (Palukaitis and García-Arenal, 2003). Thus, we have shown that infection by diverse viruses can alter host VOCs to be more attractive to bees, and that the effect occurs in distantly related plants with distinct pollination mechanisms: bean (this study), and tomato (Groen et al., 2016).
Bumblebees are generalist pollinators. This is made possible by their ability to use several sensory channels to perceive multimodal cues for flower detection such as visual and olfactory stimuli (Katzenberger et al., 2013;Glover, 2014;Lawson et al., 2017) coupled to their learning capacity that allows effective association of particular floral features with nectar and pollen rewards (Cnaani et al., 2006;Gomez et al., 2008;Eisenhardt, 2014;Konzmann and Lunau, 2014). In this study, we have discovered an innate preference of bees for odors emitted from virus-infected bean plants in the absence of visual stimuli. Although this is consistent with the idea that viruses may benefit their plant hosts by encouraging pollinator visits, further studies will show how bee foraging behavior is affected by the combination of visual and olfactory cues. In multimodal experiments, Telles et al. (2017) found that flower color had the strongest effect on bumblebee foraging but that familiar floral scents enhanced the ability of bees to locate flowers with a novel color. Virus-infected plants are often stunted making their flowers less conspicuous, and it is conceivable that odor cues may be important to make them more evident to foragers.
BCMV and CMV induced changes in volatile emission by flowering bean plants that rendered the VOC blend more attractive to bumblebees. Work with yellow-faced bumblebees (Bombus vosnesenskii Radoszkowsk: Byers et al., 2014) and black bean aphids (Aphis fabae Scopoli: Webster et al., 2010) showed that neither are attracted to plants by individual VOCs but rather by specific VOC combinations. Yellow-faced bumblebee workers preferred a combination of limonene, myrcene and (E)-ocimene over any single compound (Byers et al., 2014) while volatiles within a blend that function as host cues to the black bean aphid become non-host cues (or repellent) when presented individually (Webster et al., 2010). Our data on VOC emission from virusinfected bean plants supports the idea that the innate olfactory preferences of insects are driven by VOC blends rather than by single VOCs. Total VOC emission from virus-infected whole plants was roughly half that from healthy plants due in part to the negative effect of virus infection on plant growth, but emission of individual VOCs was also differentially affected at both whole plant and tissue dry weight levels. Two of the most abundant compounds in the VOC blend emitted by flowering bean plants, (E)-ocimene and 4-hexen-1-ol acetate, were significantly lower in all virus-infected plants. Diminishing the proportion of these VOCs in the blend from flowering bean plants is unlikely to account for the innate attraction of bumblebees per se, as the insects found the VOCs from mock-inoculated and BCMNVinfected plants equally attractive in olfactometry tests. The similar trajectories of the loading vectors for the monoterpenes (E)ocimene and D-linalool (Figure 2) indicate that these VOC outputs are positively correlated and illustrate their closely related biosynthetic pathways derived from geranyl diphosphate (Dudareva and Pichersky, 2000). D-Linalool, which is widespread among floral scents (Knudsen et al., 2006), influences various diurnal (Dötterl and Vereecken, 2010), and nocturnal bees (Krug et al., 2018) and other nocturnal visitors, including moths or bats (Dobson, 2006). However, less D-linalool was emitted by BCMV-infected and CMV-infected plants relative to healthy and BCMNV-infected plants, which is puzzling. However, variation in a few key VOCs in a floral scent blend can profoundly benefit plant competitiveness and reproductive success (Parachnowitsch et al., 2012), which leads us to speculate that (as is the situation for aphid-bean interactions: Webster et al., 2010) it is the blend of VOCs, rather than the amount of specific VOCs, that is most important in driving olfactory attraction.
VOCs emitted by all virus-infected non-flowering plants attracted bumblebees indicating that changes in volatile emission by vegetative tissues provide some olfactory cues in beanbumblebee interactions. This is consistent with work with CMVinfected tomato plants where bumblebees showed similar innate preferences for flowering and non-flowering plant VOCs (Groen et al., 2016). Although it may seem counterintuitive that bees are attracted to vegetative tissue, Kárpáti et al. (2013) demonstrated that for Manduca sexta, conspecific leaf and flower-derived volatiles combined were more attractive than the flower blend alone from the nectar sources in Nicotiana attenuata and Datura wrightii, suggesting this could be a strategy for foraging insects to optimize food location using odor cues. It is also known that leaf volatiles act in tandem with floral volatiles to attract pollinators from long and short distances. Studies in dwarf palm (Chamaerops humilis) showed that green leaf VOCs attract pollinators over longer distances than floral scents (Caissard et al., 2004;Dufaÿ et al., 2004), and green leaf VOCs may also reinforce floral visual clues (Kunze and Gumbert, 2001;Kulahci et al., 2008). Also, a recent paper has shown that bumblebees can accelerate flowering by interacting with and wounding leaves (Pashalidou et al., 2020). It will be interesting to see if this wounding behavior is influenced by plant VOCs, and if it would be affected by a plant's infection status, since in nature plants are affected by multiple biotic and abiotic stresses (Gruden et al., 2020).
Finally, revealing further insights into the effects of virus infection on the emission of pollinator-perceivable VOCs by plants may be useful in devising strategies to increase pollination services for crops, and for understanding how pathogens or other stresses may affect pollinated wild plants. Around 87.5% of flowering plants depend upon animal pollinators including many crops (Ollerton et al., 2011), and populations of pollinators, in particular those of bees, are under threat (Potts et al., 2010;Godfray et al., 2015). Studies such as ours on bee-attracting VOCs may aid in improving pollinator service for crops through plant breeding or other methods.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
NM, JCr, AM, FW, and BG conceived the work. NM and JCu collected the data. AM and NC carried out statistical analysis. NM and AM conducted the data analysis. NM, AM, and JCr drafted the manuscript. All authors contributed to manuscript editing and gave final approval for submission.

FUNDING
This work was supported by the Biotechnology and Biological Sciences Research Council (BBSRC-SCPRID Grant BB/J011762/1 and BBSRC-GCRF Grant BB/P023223/1), the Leverhulme Trust (F/09 741F), Royal Society (FCG\R1\201005), and the Cambridge-Africa ALBORADA Research Fund. FW was supported by a Royal Society-FLAIR Fellowship (FLR\R1\190462). NM was supported by studentships from the Schlumberger Foundation, Cambridge Trust, and Magdalene College Cambridge.

ACKNOWLEDGMENTS
We thank Adrienne Pate for expert technical assistance and Sanjie Jiang, and Niels Groen for useful discussions.