Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 17691.
Published online 2017 December 18. doi:  10.1038/s41598-017-17970-2
PMCID: PMC5735132

Bee pollination increases yield quantity and quality of cash crops in Burkina Faso, West Africa


Mutualistic biotic interactions as among flowering plants and their animal pollinators are a key component of biodiversity. Pollination, especially by insects, is a key element in ecosystem functioning, and hence constitutes an ecosystem service of global importance. Not only sexual reproduction of plants is ensured, but also yields are stabilized and genetic variability of crops is maintained, counteracting inbreeding depression and facilitating system resilience. While experiencing rapid environmental change, there is an increased demand for food and income security, especially in sub-Saharan communities, which are highly dependent on small scale agriculture. By combining exclusion experiments, pollinator surveys and field manipulations, this study for the first time quantifies the contribution of bee pollinators to smallholders’ production of the major cash crops, cotton and sesame, in Burkina Faso. Pollination by honeybees and wild bees significantly increased yield quantity and quality on average up to 62%, while exclusion of pollinators caused an average yield gap of 37% in cotton and 59% in sesame. Self-pollination revealed inbreeding depression effects on fruit set and low germination rates in the F1-generation. Our results highlight potential negative consequences of any pollinator decline, provoking risks to agriculture and compromising crop yields in sub-Saharan West Africa.


Insect pollination of flowering plants is a process of significant importance in terrestrial environments, and it provides vital ecosystem services for human well-being1,2 such as crop production. About 75% of agricultural crop species rely, to some degree, on animal pollination, and about one-third benefit from cross-pollination by developing higher fruit quantity and/or quality35. Among other insects, bees are known as one of the most important groups of pollinators worldwide6,7; a decrease of this pollination service could potentially reduce yields by ca. 40%3. However, non-bee insects have been shown to be important pollinators, too8. Land-use change and rapid habitat transformation during recent decades are seen as important drivers of insect pollinator declines2,9, thereby increasing the risk of future pollination deficits in areas of high, and increasing, pollination demands10. Areas where food production mostly depends on animal pollination are also those from which the fewest data are available11,12; this is due to lack of infrastructure and funding in many areas of the world, particularly in developing countries. The same areas are often poorly buffered against disruption of ecosystem service provision from whatever cause, and for this reason, effects of any ecological incidents on human well-being could be more severe here than elsewhere13.

Burkina Faso’s economic development largely depends on agriculture, with cotton as main export product (over 60% of the total export volume), making it the main cotton producer in sub-Saharan Africa14,15. Sesame ranks third among the top ten export commodities of the country. According to FAO16 cotton is sold as lint (265,500 tonnes in 2014) and seed (12,589 tonnes in 2014), sesame is sold as seed (321,837 tonnes production on 506,095 ha area harvested in 2014). Around 80% of the population relies on subsistence farming, as it is true for West Africa in general17, with the majority of agricultural production on properties of less than 2 ha up to 6 ha. A standard field of a Burkinabé farmer is of ca. 1 ha in size and a farmer usually cultivates a diversity of cotton, maize, sorghum, millet, sesame, etc. on each 1 ha field18.

Cotton and sesame are known to be able to self-pollinate, though the levels of selfing and outcrossing differ in the literature19,20. Additionally, only few studies have investigated the quality of the offspring in terms of germination rate, survival and growth performance in comparison to individuals resulting from self-pollination. Furthermore, different crop varieties can exhibit variable levels of pollinator dependence and possible yield deficits, hence an investigation of potential impacts of variety on pollination is important21,22.

Despite the variable level and mode of pollination, bees have been reported to be important pollinators of cotton and sesame2325. Some systematic assessments of pollinator species and their relative pollination service contribution to the yield quantity and quality of these two important cash crops have been carried out26,27, but estimations of the relative economic benefits of bee pollination have rarely been made. Most studies assessing the contribution of insect pollinators to yield quantity and quality were conducted in Europe and North America25,28,29. Studies in Africa mostly concentrate on South African fruit30,31 and seed crops3234. Additionally, the majority of studies are based on large-scale farms, while smallholding crop systems with fields smaller than 2 ha have received less attention in ecosystem-service research35,36. Following FAO, the insect pollination economic value (IPEV) for West Africa is estimated as 5.6 Billion US Dollars, which is the highest on the African continent along with one of the highest vulnerability rates in terms of pollination (i.e., IPEV/total production economic value)11,37.

In order to assess the supply of pollination services provided by insects to upland cotton and sesame in Burkina Faso, this study utilizes a novel combination of pollinator effectiveness assays, pollinator dependence manipulations and quality measures of F1-generation individuals from selfing or outcross-pollination. We have (1) determined the level of self-compatibility and the selfing rate of the crop varieties; (2) determined the pollinator species of cotton and sesame, and quantified their relative effectiveness; (3) quantified the demand for insect pollination services; (4) investigated signs of inbreeding depression after self-pollination; and (5) estimated the economic benefits of pollinators to smallholder production.


Study area and study sites

Field work was conducted from June to November 2015 in the southwest of Burkina Faso. This period covers the rainy season and outset of the dry season, which equally mark the flowering and harvesting season of cotton and sesame in Burkina Faso.

We selected the variety FK37 of conventional upland cotton (Gossypium hirsutum, Malvaceae), which is the most commonly used variety in Burkina Faso with more than 50% of all cotton varieties cultivated. In the study region of southwest Burkina Faso, this variety was cultivated on 100% of all cotton fields in 2015 (79% in 2017, depending on seed availability; local cotton producer society, SOFITEX, pers. comm.). For sesame (Sesamum indicum, Pedaliaceae), we selected the variety S-42, which is most widely cultivated, and covers over 80% of the cultivated areas in Burkina Faso. Three sites for exclusion experiments and pollinator monitoring in cotton and sesame fields were chosen: Dano (11°08'56.566″ N, 003°03'36.446″ W), Bontioli (10°48'26.393″ N, 003°04'39.564″ W) and Nazinga (11°06'34.998″ N, 001°29'07.181″ W), at elevations between 271 and 448 m a.s.l. All sites are located in the south-Sudanian vegetation zone and are characterized by a mosaic of agricultural land, small villages and fragments of near-natural savanna woodlands. The average proportion land cover (%) of grass-, shrub- and tree savannas of the sites amounts to 72.98% ± 18.15%. Average crop land cover (farmland and fallows) amounts to 16.86% ± 18.59%38,39. Smallholder agricultural fields were ca. 1 ha each in size, typically embedded in a matrix of other small fields of wind-pollinated cereal crops such as maize, millet and sorghum, or entomophilous crops such as beans and groundnut, with a maximum distance of 500 m (mean 81 m ± 101 m SD) to the next near natural savanna fragment.

At all sites farmers were requested to continue their usual farming practice during the study period: fertilizers were applied at the beginning of the sowing season, insecticides and fungicides were irregularly applied depending on the infestation rate and financial resources of the farmers. Weeds were removed manually. The management practices at all sites were alike. Only at the study site of Bontioli, a neighboring farmer kept 5 beehives of managed honeybees in the vicinity of about 2 km to our study fields. The beehives were not moved towards the fields.

Breeding system, seed germination and sapling performance

To investigate the level of self-pollination, self-compatibility and potential signs of inbreeding depression in the studied crop varieties, pollination experiments were conducted following Dafni40 (see Supplementary Table S1).

Cotton and sesame both have one-day flowers. All manipulations were carried out on the only day of flower anthesis, which ensured stigma receptivity and thus successful hand pollination experiments. Flowers usually open around 7:00 a.m., and pollen is released 30 to 60 min later. For our experiments, 50 plants of cotton and 50 plants of sesame were randomly selected from one field each. Five different pollination treatments were conducted per individual plant. For each treatment, we used two flowers per plant (i.e. 10 flowers per plant were treated). In total we investigated 100 flowers per treatment. Each individual plant was subjected to all pollination treatments rather than only to one treatment. This approach was made to avoid assessing the impact of the individual plant on reproductive performance rather than the impact of the treatment itself. For bagging of buds and flowers we used perforated pollination bags (hole diameter: 0.5 mm). Emasculation (removal of anthers), if required, was done in the morning between 7:00 and 7:30 a.m. before anthers released pollen. Due to high fragility of sesame flowers and difficult access to the anthers, the treatments hand-crossing (HC), hand-selfing (HS) and spontaneous selfing (SS) were carried out in a greenhouse of the Botanical Garden, University of Rostock, Germany under simulated field conditions (day/night: 32 °C/25 °C, 12 h light, 12 h dark). The open pollination treatments control (C) and natural outcrossing (OC) were conducted with the same seed variety in Burkina Faso since they required the presence of natural pollinators.

When fruits were mature, we recorded whether flowers set fruit, and if so, the number and weight of seeds per fruit. Mean fruit set, overall mean number of seeds per fruit and mean seed weight were then calculated for each treatment. Fruit set is calculated as number of fruits divided by the number of flowers.Two indices related to selfing were calculated: 1) the selfing rate “S”41 calculated as:


where Px are seeds resulting from open cross-pollination (OC), Pc are seeds resulting from open pollination (C) and Ps are seeds resulting from self-pollination (SS), and 2) the index of self-incompatibility “ISI”42,43 as the average fruit set by hand self-pollination (HS) divided by the average fruit set by hand cross-pollination (HC).

Inbreeding depression for fruit set, seed production, seed weight, and seed germination were calculated by the inbreeding depression index44 as the relative performance of cross types (RP) following the formula:


where Ws is the mean performance of selfed progeny (SS) and Wo is the mean performance of out-crossed progeny (OC). This index varies from −1 to 1. Positive values indicate that outcrossed progeny outperform selfed progeny (i.e. selfed progeny suffers from inbreeding depression), negative values indicate the opposite case44.

To evaluate seed germination, 270–370 cotton seeds from each pollination treatment (corresponding to the number of available seeds, total number of seeds: 1650) were placed in Petri dishes (10 seeds per dish) on filter paper moistened with purified water and maintained in a climate conditioning chamber at temperatures of 32 °C (day) and 25 °C (night) with artificial lighting and a photoperiod of 12 h from 6:00 a.m.–6:00 p.m.45.

Each day, we counted germinated seeds until all seeds had either germinated or decayed.

The same protocol applied to sesame, but with placing the seeds (210–400 sesame seeds from each pollination treatment, total number of seeds: 2249) on substrates (pots) in the greenhouse under similar conditions. All Petri dishes and pots were manually randomized every day (i.e. array in climate chamber and greenhouse, respectively) to avoid impacts of varying microclimatic and light conditions.

After germination, 50 cotton seedlings per treatment were transferred to pots filled with standard substrate (sand-turf mixture) in the greenhouse where the F1-individuals were cultivated. For sesame, seedlings remained in their pots where they already had germinated. After 35 days of cultivation, height and number of leaves per sapling were measured.

Pollinator species and their effectiveness

To assess pollinating species and their ability to pollinate cotton and sesame flowers, we determined the flower-visiting species and their pollination effectiveness in terms of fruit set and fruit quality. At all three study sites, 100 flowers each of cotton and sesame in bud stage were bagged with pollination bags and marked. In total we assessed flower visitors and their effectiveness as pollinators on 300 flowers of cotton and 300 flowers of sesame. Every day at 7:00 a.m., we checked if flowers had opened. If open, the bag was removed and we waited at a slight distance for the first pollinator to enter the focal flower. When leaving the flower, the pollinator was caught, etherized and stored in 70% ethanol. Thereon, the flower was bagged again after the insect visit. Care was taken that weather conditions were similar (sunny to bright, no rain) on all sampling days. Each pollinator specimen was numbered according to the visited flower and then pinned and identified to genus or species (voucher specimens are held at the Royal Belgian Institute of Natural Sciences, Brussels). At the end of the season, we recorded whether each visited flower set fruit, and if so the quality parameters such as the number of intact and non-intact seeds, seed weight, and additionally for cotton, fiber weight (lint weight) per fruit. For each visitor species, we calculated mean values and standard deviations for all parameters.

Pollinator dependence

To investigate to which extent the focal cash crop species benefit from insect-mediated outcross pollination in comparison to simulated total absence of pollinators, pollinator exclusion experiments were conducted on 11 fields of cotton and sesame each: 4 fields of cotton and sesame each were selected at a maximum distance of 1 km to the nearest semi-natural savanna habitat in “Nazinga” and “Bontioli”. However, at “Dano”, 3 fields per crop species were chosen due to infrastructural restrictions.

On each field, 25 plants were selected randomly. Three pollination treatments were conducted per plant: (1) natural open pollination, OPEN; (2) spontaneous self-pollination, SELF, and (3) natural outcross-pollination, CROSS (see Supplementary Table S1). The natural open pollination (control) includes self- and cross-pollination. Due to the selfing ability of both crop species, we added the outcrossing treatment to investigate the pure contribution of insect-mediated pollen transfer from other conspecifics (cross-pollination). Each pollination treatment was applied to two flowers per plant. In total we had 275 plants of cotton and 275 plants of sesame (100 plants in Nazinga and 100 plants in Bontioli due to 25 plants on 4 fields, and 75 plants in Dano due to 25 plants on 3 fields). Since we used two flowers per plant, the experiments were executed with a total of 550 flowers per crop and treatment.

Cotton plants usually produce only one open flower per day. In rare cases, where two flowers were open on the same day, one of the flowers was bagged and used for an exclusion treatment to prevent self-pollination that occurs from pollinator movement between flowers on the same plant (geitonogamy). In sesame, inflorescences appear to bare several open flowers, most of which are old flowers from the day before, and not being receptive to pollen anymore. We only used flowers that were produced during the peak of the flowering season from mid-July until mid-August (i.e. not the very first and the very last flowers produced by a plant), when the plants had reached their full height, to avoid impacts of potential resource limitation on fruit set in early growth and late fruit ripening periods. At the end of the season, we recorded whether each treated flower set fruit, and the quality parameters: 1) number of intact seeds, 2) number of non-intact seeds, 3) seed weight, additionally 4a) fiber weight (lint weight) per fruit for cotton and 4b) fruit weight for sesame.

Economic valuation of pollination – yield gap

Economic benefits of pollination services to producers were calculated for each crop by comparing the fruit quantity from open natural pollination (mixed self- and cross-pollination), outcross-pollination and pollinator-exclusion treatments (spontaneous self-pollination)46. Cotton yield is determined as the weight of cotton fibers (lint weight) and seeds together, with the pericarp removed. Hence, we summed up fiber and seed weight recorded for each fruit of each treatment from the pollinator exclusion experiments, calculated the mean values and the difference of yield resulting from outcrossing and selfing to the natural open pollination (control). In sesame, only the seeds are sold. Hence, we used the seed weight data for economic valuation. Yield gap, defined here as the difference in crop yield between natural open pollination and self-pollination (pollinator exclusion), was calculated as kg*ha−1 based on average total yield quantity of the study fields (data obtained from the local cotton producer society SOFITEX and the local sesame farmer community). Based on the market prices in Burkina Faso in 2015, yield gaps and economic benefits through insect pollination were valued by multiplying the yield quantity [kg*ha−1] by the lowest national market price per kg (cotton: 0.38 USD or 225 Francs CFA, sesame: 0.43 USD or 250 Francs CFA)47 to state realistic values. The farmers´ expenses for seeds, fertilizers and pesticides per hectare (cotton: 140.87 USD or 82875 Francs CFA, sesame: 34.00 USD or 20000 Francs CFA)47 were subtracted. Human labor was not monetarily considered in this calculation.

Data analysis

We fitted generalized linear mixed models (GLMER) to analyse the effect of pollination treatment on fruit set with a binomial data distribution (fruit set = 1, no fruit set = 0). To analyse and compare the effect of pollination treatment on fruit quality parameters of cotton and sesame, we used linear mixed-effects models (LMER) considering in both models field and plant as random grouping factors. Multiple comparisons were conducted using Tukey procedure with P-value  0.05 as threshold for significance. All quantity and quality parameter data refer to field measures at site in Burkina Faso. Only seed germination of sesame seeds deriving from self-pollination refers to the seeds that were produced in greenhouses in Germany. Differences in group means of seed germination after treatment (manually randomized during data collection) were analysed using Kruskall-Wallis One-Way-ANOVA on Ranks due to non-normal distribution of data and pairwise multiple comparison procedures (Tukey Test or Dunn´s respectively) with P-value of  0.05 as the threshold of significance. Data of the follow-up seedling performance experiment were analysed using LMER with Petri dish (cotton) or pot (sesame) as random grouping factor. Linear mixed-effects models were further used to calculate pollinator effectiveness by separating four bee species that visited cotton flowers most frequently from wasps (exclusively Scoliidae species) and from remaining visitors. The dataset of the latter was pooled due to low visitor numbers (less than ten visits per species). The effect of single pollinator species on fruit set was analysed using GLMER and quality parameters, namely seed weight, fibre weight and number of non-intact seeds, were analysed using LMER with field as random grouping factor. For sesame, no statistical test regarding pollinator efficiency was feasible due to inhomogeneity of visitor numbers. Statistical analyses and in part figure production were conducted in R version 3.3.048 with additional functions provided by the R packages MASS49, lme450, lmerTest51 and multcomp52.

Data availability Statement

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.


Breeding system, seed germination and sapling performance

In both cotton and sesame, pollination treatments revealed a mixed mating system. The cotton variety used in this study has a high ability to set fruit without pollinators via spontaneous selfing (61%), whereas this ability was lower in sesame (20%). The selfing rate ‘S’ accounts for 0.143 in cotton and 0.016 in sesame (i.e. ca. 14% of cotton seed production and ca. 2% of sesame seed production per fruit derive from self-pollination). The hand pollination experiments revealed an index of self-compatibility (ISI) of 0.91 in cotton and 0.71 in sesame. Hence, both crops are partially self-compatible (ISI of 1 = fully self-compatible).

The relative performance of cross types RP (inbreeding depression index) in cotton was determined as 0.197 for fruit set, 0.259 for seed production, 0.206 for seed weight and 0.347 for seed germination rate. In sesame, RP amounted to 0.369 for fruit set, 0.45 for seed production, 0.71 for seed weight and 0.97 for seed germination rate (positive values from 0 to 1 indicate inbreeding depression of selfed progeny, negative values to −1 would indicate outbreeding depression). The germination rate of control seeds was 61.35% in cotton and 26.29% in sesame (C vs. OC, cotton: ANOVA, F(4,160) = 0.536, ns.; sesame: Q = 4.753, P < 0.05), whereas the germination rate of seeds resulting from outcross-pollination was highest in both crops (cotton: 72.16%, sesame: 62.16%). Seeds resulting from spontaneous self-pollination germinated least (cotton: 47.09%, sesame: 1.81%; OC vs. SS: cotton: ANOVA, F(4,160) = 0.643, P < 0.001; sesame: Q = 7.918, P < 0.05).

In cotton, the saplings of the F1-generation resulting from self-pollination showed a trend towards reduced height and number of leaves in comparison to outcrossed saplings (OC vs. SS: height: LMER, t 26 = −1.8, P = 0.08; number of leaves: LMER, t 23 = −1.7, P = 0.09). Whereas outcrossed saplings had a mean height of 33.20 cm (SD = 12.49) and mean number of leaves of 5.89 (SD = 2.08), saplings from self-pollination were smaller in height (mean: 23.93, SD = 17.76) and had fewer leaves (mean: 4.47, SD = 3.09).

In sesame, no significant difference in the performance of the offspring was found between self-pollination and outcrossing treatments.

Pollinator species and their effectiveness

The diversity of flower visitor species differed between cotton and sesame. We recorded the semi-wild (few colonies are managed, while the majority is feral) honeybee Apis mellifera, twenty-six wild bee species and four species of wasps as flower visitors in cotton. Six wild bee species did not initiate fruit set; all other species are considered as pollinators. Flowers of sesame were visited by Apis mellifera and four wild bee species, but only the honeybees and two wild bee species initiated fruit set and are thus considered as pollinators. No other animal taxa or insects were observed visiting the cotton and sesame flowers during the effectiveness experiments (Table 1).

Table 1
Taxonomic list of Hymenoptera visitor species to the flowers of cotton (Gossypium hirsutum, Malvaceae) and sesame (Sesamum indicum, Pedaliaceae), their number (n) of visits and their ability to initiate fruit set. Crosses indicate the initiation of fruit ...

The effectiveness assays revealed that not all pollinators were equally effective, and the relative abundance of different pollinator species varied strongly. In cotton, the most abundant pollinator species was the honeybee Apis mellifera with 150 visits followed by the native wild bee species Tetralonia fraterna (long-horned bee; 59 visits), the small carpenter bee Ceratina spec.1 (15 visits) and the stingless bee Hypotrigona gribodoi (11 visits). All other twenty-three wild bee species and four Scoliid wasp species were infrequent flower-visitors (Table 1). Only pollination by honeybees significantly increased fruit set, whereas all other pollinator species had no significant effect on fruit set in cotton (Table 2). However, their effectiveness differed regarding other quality parameters: fruits resulting from pollination by T. fraterna had significantly fewer non-intact seeds (LMER, t 59 = −3.43, P = 0.002) than fruits after pollination by honeybees (LMER, t 150 = 7.94, P = 0.013) and the other pollinator species. In comparison to fruits after pollination by the other wild bee and wasp species, seed weight and fibre weight of fruits significantly increased after pollination by honeybees (seed weight: LMER, t 150 = 13.70, P = 0.004; fibre weight: LMER, t 150 = 13.57, P = 0.005) and by T. fraterna (seed weight: LMER, t 59 = 2.67, P = 0.008; fibre weight: LMER, t 59 = 2.83, P = 0.005; Table 2, see Supplementary Table S2).

Table 2
Most frequent pollinator species of cotton in Burkina Faso and their effectiveness in terms of fruit set and quality parameters (n = 300 flowers; N = number of pollinator visits). Apis mellifera is a semi-wild honeybee ...

In sesame, honeybees were by far the most abundant (290 visits) and most effective pollinators in terms of seed- and fruit weight. The other bee pollinator species were infrequent visitors. Despite initiating fruit set, fruit and seed quality decreased. After pollination by the Mmason bee Chalicodoma mephistolica the number of non-intact seeds increased, and seed and fruit weight were decreased. Pollination by the stingless bee Hypotrigona gribodoi was not effective (see Supplementary Table S2).

Pollinator dependency

Pollinator exclusion (spontaneous self-pollination, SELF) significantly decreased fruit set and fruit quality in both cotton and sesame. In cotton, natural pollination (OPEN) resulted in a mean fruit set of 48.73%, fruit set from self-pollination due to pollinator exclusion accounted for 37.45% (GLMER, z = −4.1, P < 0.001). In sesame, mean fruit set from natural pollination (OPEN) was 45.82%, in contrast to fruit set from self-pollination, which was 20.36% (GLMER, z = −9.3, P < 0.001). Similar results were found for all quality parameters (Fig. 1). In particular, the economically important parameter fibre weight in cotton (seed weight in sesame) significantly decreased by ca. 33% (ca. 59% in sesame) when pollinators were excluded (LMER, cotton: t 639 = −8.8, P < 0.001; sesame: t 486 = −9.8, P < 0.001). Fruits resulting from self-pollination contained eight times (cotton) or three times (sesame) as many non-intact seeds as fruits of the natural open pollination treatment (cotton: LMER, t 668 = 26.9, P < 0.001; sesame: LMER, t 427 = 11.4, P < 0.001). Outcross-pollination did not result in more abundant mean fruit set in both crops in comparison to natural pollination (cotton: 46.36%, sesame: 32.36%). However, cross-pollination by bees significantly increased fruit quality. In cotton, fruits resulting from outcross-pollination had significantly fewer non-intact seeds (LMER, t 641 = −3.2, P < 0.01) and increased seed weight (LMER, t 613 = 8.5, P < 0.001) and fibre weight (LMER, t 615 = 7.1, P < 0.001) in comparison to the natural open pollination treatment (OPEN). In sesame, seed weight from outcross-pollination was ca. 41% higher than from mixed pollination of the control treatment (LMER, t 504 = 8.6, P < 0.001; Fig. 1, see Supplementary Table S3).

Figure 1
Fruit quality parameters of cotton (a) and sesame (b) resulting from pollinator dependency experiments in Burkina Faso, 2015. 550 flowers (2 flowers per treatment of 275 plants) from 11 fields per crop species were subjected to each pollination treatment: ...

Economic valuation of bee pollination

The economic impact of bee pollination on producer profits was significantly beneficial in both crop species, and a loss of pollination service by bees would cause severe yield gaps.

The national average yield of cotton in Burkina Faso in 201547 accounts for 1043 kg*ha−1. The mean yield of our study fields amounted to 953.91 kg*ha−1 (SD = 474.04 kg*ha−1; local cotton producer society Sofitex, unpublished data). This value corresponds to the control treatment. According to the national market prices of the Burkina Faso 2015–2016 cotton campaign, this translates into earnings of 364.39 USD per ha (SD = 181.08 USD per ha) were famers expenses of 140.87 USD per ha are still to be subtracted.

According to the data obtained in the pollinator experiments, an increase of yield quantity by bee-mediated outcrossing between 27% and 31% for cotton was constituted, which translates into additional earnings between 98.38 USD and 112.96 USD per ha.

On the other hand, the scenario of pollinator loss, and hence pure self-pollination of cotton flowers revealed a decrease of yield quantity between 33% and 43% This translates into a yield gap and net loss between 120.35 USD and 156.69 USD per ha. For sesame, the national average yield in Burkina Faso in 201547 stands at 250–350 kg*ha−1. The mean sesame yield (seed weight) of our study fields amounts to 202.2 kg*ha−1 (SD = 44.26 kg*ha−1; local sesame producers, unpublished data). This value corresponds to the control treatment and equates to 86.95 USD per ha (SD = 19.03 USD per ha), were farmers expenses of 52.95 USD per ha are still to be subtracted.

Outcross-pollination by bees increased yield quantity between 37 and 42%, which translates into additional earnings between 32.17 USD and 36.52 USD per ha.

A lack of pollinators and their service would cause a decrease in yield quantity between 50 and 87% and would lead to a yield gap and net loss between 43.47 USD and 75,64 USD per ha.


The breeding system experiments have revealed that both investigated cotton and sesame varieties are partially self-compatible and facultatively xenogamous. They have a mixed mating system and produce fruits through self- and cross-pollination. The selfing rate is high in cotton and less pronounced in sesame. This is in line with other studies on cotton stating that cotton is mainly self-pollinating although cross-pollination occurs20,23. Cross-pollination is mostly attributed to insect pollination because cotton pollen is too large and heavy to be easily carried by air currents53. Studies on Gossypium hirsutum in Greece and California revealed less than 10% outcross pollination. They both state that cross-pollination of cotton highly depends on pollinator activity and is low in the absence of pollinators54,55. Our study revealed a relatively high rate of outcross-pollination which might be due to high pollinator abundances at our study sites. Furthermore, the cotton fields in our study had a size of about 1 ha and were embedded in a landscape mosaic of various other crop fields and savanna. In contrast, the studies in Greece and California were conducted on large holdings. Large monocultures may decrease pollinator abundance and hence lead to lower cross-pollination56. For sesame, our current finding contradicts a study stating that sesame is predominantly autogamous57, but is in line with another study revealing best performance in sesame in Brazil after cross-pollination58. Nevertheless, crossing rates reported in some studies ranged from less than 10% to 68%59,60. The outcrossing rate of the sesame variety used in this study was 32.36%. Such high variations in levels of pollinator dependence for outcross pollination is typical of crop species with mixed mating systems19 and evidencing the need for breeding system studies and pollinator dependence experiments for each variety and environmental setting.

Despite high selfing rates, self-pollination led to clear signs of inbreeding depression, being highest in both species for number of seeds, and additionally, in cotton, for fruit set and in sesame for seed germination.

We recorded honeybees and a plethora of wild bees as pollinators of cotton and sesame at our study sites, with honeybees and the long-horned wild bee Tetralonia fraterna being most effective in terms of fruit set and fruit quality. Since flower visitor abundances were very uneven, one might argue that abundance might account for pollinator efficiency. This might be true for the two most efficient species recorded. In a previous study based on pantraps Apis mellifera (364 specimens in cotton fields, 142 in sesame) and Tetralonia fraterna (120 specimens in cotton fields) belonged to the top three most abundant bee species in cotton fields assessed in our study area within two sampling years (2014 and 2015) (KStein, data unpublished). The pantrap study revealed that the stingless bee Hypotrigona gribodoi was most abundant in both cotton (4253 specimens) and sesame (165 specimens) fields. In the present study though, this species turned out to be less efficient, indicating that sole abundance does not stand for species efficiency. It might rather be the pollen collecting behaviour and body size that account for bee species efficiency as pollinators. Pantraps do not allow for concluding on a species to be a pollinator, since the insects might have been caught by accident and do not necessarily need to be pollinators. This method can only be used to assess the community of “potential pollinators”. Our finding underlines the importance of direct pollinator observations on the flowers and an assessment of their pollination success as we did with the pollinator efficiency approach.

Many native bee species other than honeybees have been shown to visit crops61 and have been recognized for their role in increasing and stabilizing crop-pollination services62,63. Bees, and in particular wild bees, are known to improve seed set, quality, and commercial value of a variety of crops29,64. A meta-study including 90 studies from five continents found that wild bee communities contributed about the same to crop production as managed honeybees did65. Hence, wild bees complement, in a number of ways, the service provided by honeybees; first biologically, by enhancing the efficacy of honeybee pollination in some cases6, and then economically, by insuring against pollination shortages. Availability of accurate estimates of this value could improve land use planning by quantifying the costs and benefits of conserving habitats for pollinators in agricultural landscapes.

Our research revealed that cross-pollination by honeybees and wild bees significantly increased yield quality in both cotton and sesame, which is in line with the findings of other studies on cotton23 and on sesame58. However, these studies were carried out in Australia and Brazil and our results reveal even higher benefits from bee pollination in terms of fruit quality. To our knowledge, our study is the first from the sub-Saharan West African country of Burkina Faso. Quality significantly increased for the economically important parameters of fibre (lint) and seed weight by two thirds in cotton, and seed weight was tripled in sesame. A higher number of intact seeds and increased seed weight per fruit are desirable features for cotton and sesame from both the commercial and ecological point of view, which is an important indicator of plant reproductive success. A loss of bee pollinators concomitant with self-pollination of the crop species would result in yield gaps between 33% and 43% for cotton and 50% and 87% for sesame, respectively. Furthermore, seeds from outcross-pollination were more viable and vital, and germinated more readily, contrary to seeds from self-pollination. This finding is relevant for the smallholders and cotton companies, since they use the seeds from the previous harvest for next year’s sowing season. A high portion of outcrossed seeds significantly enhances the success of next year’s cultivation. Thus, our study clearly reveals the value of bee pollination to cotton and sesame production in smallholder systems in sub-Saharan Africa. Since it lacks detailed social data on the livelihood of the farmers, the economic evaluation in this study only allows broad conclusions on potentially increased earnings or losses from existing or lacking pollination services in relation to national statistics.

In the context of human population growth and land conversion with natural habitat loss as one of the main drivers of bee decline66,67, and increasing climate unpredictability68 especially for West Africa, targeted establishment and preservation of semi-natural habitats would consistently increase the diversity and abundance of wild pollinators62. The expected agricultural loss in the absence of animal pollination is estimated at 5–7.5% in Burkina Faso56. As with a number of previous studies, estimates of the economic benefits of pollination services are limited by the assumption of constant market prices69. In particular, the benefits reported here may vary depending on the presence of other inputs such as fertilizers and pesticides or ecosystem services70. For our calculation, we have used the market prices of 2nd class quality cotton and the lowest market price of sesame seeds.

Additionally, world cotton production is projected to grow annually by 2.1% over the next ten years (2.3% in West Africa). Besides Brazil, least developed countries in sub-Saharan Africa are also expected to augment their cotton exports, increasing their share of world trade growing from 8 to 9% by 2024. Cotton consumption is limited throughout sub-Saharan Africa, and many countries export virtually their entire production71. In sesame, the annual production growth rate in Burkina Faso is estimated to be 18.4%15. Based on a sesame seeds export quantity of 58,650 tons in 2012, which equates to an export value of 56,730^1000 USD15, the loss of pollination service would lead to a yield loss of 34,498 tons on a national scale and would have severe economic impact on the country.


By combining results from exclusion experiments, pollinator surveys and field manipulations, this study for the first time quantifies the contribution of bee pollinators to smallholders’ cash crop production in Burkina Faso. Measures of pollinator diversity and effectiveness were combined and compared between species to quantify their contribution to cash crop production and economic output in West Africa. This study is also unique in that pollinator species diversity and efficiency were assessed, and that an in-depth assessment was provided, tracking multiple steps of the different modes of pollination from the breeding system to the F1-generation.

As our knowledge on pollination efficiency of different pollinators to different crops grows and consolidates globally72, the concepts used in this study can be applied to better quantify economic impacts of different components of the pollinator community on crop production. Ultimately, this can result in more holistic models of pollination service provision and better forecasting of the risks of pollinator declines68,69.

These findings have implications for the management of insect pollination services in cash crop fields such as cotton and sesame, and highlight potential consequences of any decline in specific taxa, thereby advocating the conservation of bee populations.

Electronic supplementary material


The study was realized within the WASCAL programme, financed by the German Federal Ministry of Education and Research. We thank the “Ministère de la Recherche Scientifique et de l´Innovation”, Burkina Faso, for research permit, Adjima Thiombiano of the University of Ouagadougou for his support, Alain Pauly from the Royal Belgian Institute of Natural Sciences Brussels, for his help regarding determination of the bee species, Adewole Olagoke of the Technische Universität Dresden for proof reading, the student assistants Hannah Parnow-Kloth, Paul Goschnick, Sarah Mueller and Jasper Sehrt for conducting breeding system and germination studies, numerous field assistants and farmers in Burkina Faso, and the local cotton producer society Sofitex for their collaboration. KStein and DC were funded by WASCAL.We acknowledge financial support by Deutsche Forschungsgemeinschaft and Universität Rostock within the funding programme Open Access Publishing.

Author Contributions

Author Contributions

K. Stein, K.E.L. and S.K. designed the study. K. Stein and D.C. collected data in the field. D.C. determined bee species. K. Stein, K. Stenchly and A.L. analysed and plotted output data. S.P. and D.G. contributed towards data sorting and focus of the manuscript and supervised pollination and germination experiments in the greenhouses. K.S. wrote the first draft of the manuscript, and all authors contributed substantially to revisions.


Competing Interests

The authors declare that they have no competing interests.


Electronic supplementary material

Supplementary information accompanies this paper at 10.1038/s41598-017-17970-2.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Garibaldi LA, Aizen MA, Klein AM, Cunningham SA, Harder LD. Global growth and stability of agricultural yield decrease with pollinator dependence. Proc. Natl. Acad. Sci. USA. 2011;108:5909–5914. doi: 10.1073/pnas.1012431108. [PubMed] [Cross Ref]
2. Garibaldi LA, et al. From research to action: enhancing crop yield through wild pollinators. Front. Ecol. Environ. 2014;12:439–447. doi: 10.1890/130330. [Cross Ref]
3. Klein AM, et al. Importance of pollinators in changing landscapes for world crops. Proc. R. Soc. B. 2007;274:303–313. doi: 10.1098/rspb.2006.3721. [PMC free article] [PubMed] [Cross Ref]
4. Aizen MA, Garibaldi LA, Cunningham SA, Klein AM. How much does agriculture depend on pollinators? Lessons from long-term trends in crop production. Ann. Bot.-London. 2009;103:1579–1588. doi: 10.1093/aob/mcp076. [PMC free article] [PubMed] [Cross Ref]
5. Classen A, et al. Complementary ecosystem services provided by pest predators and pollinators increase quantity and quality of coffee yields. Proc. R. Soc. B. 2014;281:20133148. doi: 10.1098/rspb.2013.3148. [PMC free article] [PubMed] [Cross Ref]
6. Greenleaf SS, Kremen C. Wild bees enhance honey bees´ pollination of hybrid sunflower. Proc. Natl. Acad. Sci. USA. 2006;103:13890–13895. doi: 10.1073/pnas.0600929103. [PubMed] [Cross Ref]
7. Winfree R, Williams NM, Dushoff J, Kremen C. Native bees provide insurance against ongoing honeybee losses. Ecol. Lett. 2007;10:1105–1113. doi: 10.1111/j.1461-0248.2007.01110.x. [PubMed] [Cross Ref]
8. Rader R, et al. Non-bee insects are important contributors to global crop pollination. Proc. Natl. Acad. Sci. USA. 2015;113:146–151. doi: 10.1073/pnas.1517092112. [PubMed] [Cross Ref]
9. Ollerton J, Erenler H, Edwards M, Crockett R. Extinctions of aculeate pollinators in Britain and the role of large-scale agricultural changes. Science. 2014;346:1360–1362. doi: 10.1126/science.1257259. [PubMed] [Cross Ref]
10. Lautenbach S, Seppelt R, Liebscher J, Dormann CF. Spatial and temporal trends of global pollination benefit. PLoS ONE. 2012;7:e35954. doi: 10.1371/journal.pone.0035954. [PMC free article] [PubMed] [Cross Ref]
11. Gallai N, Salles J-M, Settele J, Vaissière BE. Economic valuation of the vulnerability of world agriculture confronted with pollinator decline. Ecol. Econ. 2009;68:810–821. doi: 10.1016/j.ecolecon.2008.06.014. [Cross Ref]
12. Archer CR, Pirk CWW, Carvalheiro LG, Nicolson SW. Economic and ecological implications of geographic bias in pollinator ecology in the light of pollinator declines. Oikos. 2014;123:401–407. doi: 10.1111/j.1600-0706.2013.00949.x. [Cross Ref]
13. De Palma A, et al. Predicting bee community responses to land-use changes: Effects of geographic and taxonomic biases. Scientific Reports. 2016;6:31153. doi: 10.1038/srep31153. [PMC free article] [PubMed] [Cross Ref]
14. Kampmann, D. Biodiversity Atlas of West Africa Vol.2 Burkina Faso (ed. Thiombiano, A.) 484–487 (BIOTA, 2010).
15. FAO - Food and Agriculture Organization of the United Nations. Statistical database. (2013).
16. FAO - Food and Agriculture Organization of the United Nations. Statistical database. (2014).
17. GIZ – Deutsche Gesellschaft für Internationale Zusammenarbeit GMbH. Africa, Burkina Faso. (2016).
18. FAO - Food and Agriculture Organization of the United Nations. Cadre de Programmation Pays 2013–2015, Burkina Faso. (2013).
19. Pathirana R. Natural cross-pollination in sesame (Sesamum indicum L.) Plant Breed. 1994;112:167–170. doi: 10.1111/j.1439-0523.1994.tb00665.x. [Cross Ref]
20. Abrol DP. Pollination biology: biodiversity, conservation and agricultural production. Q. Rev. Biol. 2013;88:132.
21. Hudewenz A, Pufal G, Bogeholz A-L, Klein A-M. Cross-pollination benefits differ among oilseed rape varieties. The Journal of Agricultural Science. 2014;152:770–778. doi: 10.1017/S0021859613000440. [Cross Ref]
22. Benjamin FE, Winfree R. Lack of pollinators limits fruit production in commercial blueberry (Vaccinium corymbosum) Environ. Entomol. 2014;43:1574–1583. doi: 10.1603/EN13314. [PubMed] [Cross Ref]
23. Rhodes J. Cotton pollination by honeybees. Aust. J. Exp. Agric. 2002;42:513–518. doi: 10.1071/EA01063. [Cross Ref]
24. Kumar R, Lenin JK. Insect pollinators and effects of cross pollination on yield attributes of sesame (Sesamum indicum L.) Indian Bee J. 2000;62:75–80.
25. Vaissière BE, Moffett JO, Loper GM. Honey bees as pollinators for hybrid cotton seed production on the Texas High Plains. Agron. J. 1984;76:1005–1010. doi: 10.2134/agronj1984.00021962007600060033x. [Cross Ref]
26. Waller D, Vaissière BE, Moffett JO, Martin JH. Comparison of Carpenter Bees (Xylocopa varipuncta Patton) (Hymenoptera: Anthophoridae) and Honey Bees (Apis mellifera L.) (Hymenoptera: Apidae) as pollinators of male-sterile cotton in cages. J. Econ. Entomol. 1985;78:558–561. doi: 10.1093/jee/78.3.558. [Cross Ref]
27. Berger LA, Vaissière BE, Moffett JO, Merritt SJ. Bombus spp. (Hymenoptera: Apidae) as pollinators of male-sterile upland cotton on the Texas High Plains. Environ. Entomol. 1988;17:789–794. doi: 10.1093/ee/17.5.789. [Cross Ref]
28. Bommarco R, Marini L, Vaissière BE. Insect pollination enhances seed yield, quality, and market value in oilseed rape. Oecologia. 2012;169:1025–1032. doi: 10.1007/s00442-012-2271-6. [PubMed] [Cross Ref]
29. Klatt BK, et al. Bee pollination improves crop quality, shelf life and commercial value Proc. R. Soc. B. Biol. Sci. 2014;281:20132440. doi: 10.1098/rspb.2013.2440. [PMC free article] [PubMed] [Cross Ref]
30. Carvalheiro LG, Seymour CL, Nicolson SW, Veldtman R. Creating patches of native flowers facilitates crop pollination in large agricultural fields: mango as a case study. J. Appl. Ecol. 2012;49:1373–1383. doi: 10.1111/j.1365-2664.2012.02217.x. [Cross Ref]
31. 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. doi: 10.1111/j.1365-2664.2010.01829.x. [Cross Ref]
32. Allsopp MH, DeLange WJ, Veldtman R. Valuing insect pollination services with cost of replacement. PloS ONE. 2008;3:e3128. doi: 10.1371/journal.pone.0003128. [PMC free article] [PubMed] [Cross Ref]
33. Melin A, Rouget M, Midgley J, Donaldson JS. Pollination ecosystem services in South African agricultural systems. S. Afr. J. Sci. 2014;110:1–9. doi: 10.1590/sajs.2014/20140078. [Cross Ref]
34. Carvalheiro LG, et al. Natural and within-farmland biodiversity enhances crop productivity. Ecol. Lett. 2011;14:251–9. doi: 10.1111/j.1461-0248.2010.01579.x. [PubMed] [Cross Ref]
35. Steward PR, et al. Pollination and biological control research: are we neglecting two billion smallholders. Agriculture & Food Security. 2014;3:5. doi: 10.1186/2048-7010-3-5. [Cross Ref]
36. Garibaldi LA, et al. Mutually beneficial pollinator diversity and crop yield outcomes in small and large farms. Science. 2016;351:388–391. doi: 10.1126/science.aac7287. [PubMed] [Cross Ref]
37. FAO - Food and Agriculture Organization of the United Nations. Statistical database. http://faostat3.fao.orgAgricultural data/Agricultural production/Crops primary. (2007).
38. Dimobe K, et al. Identification of driving factors of land degradation and deforestation in the Wildlife Reserve of Bontioli (Burkina Faso, WestAfrica) Glob. Ecol. Conserv. 2015;4:559–571. doi: 10.1016/j.gecco.2015.10.006. [Cross Ref]
39. Dimobe K, et al. Spatio-temporal dynamics in land use and habitat fragmentation within a protected area dedicated to tourism in a Sudanian savanna of West Africa. Journal of Landscape Ecology. 2017;10:5–25. doi: 10.1515/jlecol-2017-0011. [Cross Ref]
40. Dafni, A. Pollination ecology: A practical approach, Ch. 8, 36 (Oxford University Press, 1992).
41. Charlesworth D, Charlesworth B. The effect of investment in attractive structures on allocation to male and female functions in plants. Evolution. 1987;41:948–968. doi: 10.1111/j.1558-5646.1987.tb05869.x. [PubMed] [Cross Ref]
42. Zapata TR, Arroyo MTK. Plant reproductive ecology of a secondary deciduous tropical forest in Venezuela. Biotropica. 1978;10:221–230. doi: 10.2307/2387907. [Cross Ref]
43. Schmid S, Schmid VS, Zillikens A, Harter-Marques B, Steiner J. Bimodal pollination system of the bromeliad Aechmea nudicaulis involving hummingbirds and bees. Plant Biol. 2011;13:41–50. doi: 10.1111/j.1438-8677.2010.00348.x. [PubMed] [Cross Ref]
44. Ågren J, Schemske DW. Outcrossing rate and inbreeding depression in two annual monoecious herbs, Begonia hirsuta and B. semiovata. Evolution. 1993;47:125–135. doi: 10.1111/j.1558-5646.1993.tb01204.x. [PubMed] [Cross Ref]
45. Ali MA, Khan IA, Awan SI, Ali S, Niaz S. Genetics of fibre quality traits in cotton (Gossypium hirsutum L.) Aust. J. Crop. Sci. 2008;2:10–17.
46. Garratt MPD, et al. Avoiding a bad apple: Insect pollination enhances fruit quality and economic value. Agric. Ecosyst. Environ. 2014;184:34–40. doi: 10.1016/j.agee.2013.10.032. [PMC free article] [PubMed] [Cross Ref]
47. Banque Nationale de Développement Agricole – BNDA, Bamako, Mali Agricultural Commodity Technical Cards from BNDA_Coton and _Sesame. (2015).
48. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. (2016).
49. Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S. (Springer, 2002).
50. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015;67:1–48. doi: 10.18637/jss.v067.i01. [Cross Ref]
51. Kuznetsova A, Brockhoff PB, Christensen RH. B. lmerTest: Tests in linear mixed effects models. R package version. 2016;2:0–30.
52. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom. J. 2008;50:346–363. doi: 10.1002/bimj.200810425. [PubMed] [Cross Ref]
53. Umbeck PF, et al. Degree of pollen dispersal by insects from a field test of genetically engineered cotton. J. Econ. Entomol. 1991;84:1943–1950. doi: 10.1093/jee/84.6.1943. [Cross Ref]
54. Xanthopoulos FP, Kechagia UE. Natural crossing in cotton (Gossypium hirsutum L.) Aust. J. Agric. Res. 2000;51:9789–983. doi: 10.1071/AR00026. [Cross Ref]
55. Van Deynze AE, Sundstrom FJ, Bradford KJ. Pollen-mediated gene flow in California cotton depends on pollinator activity. Crop Science. 2004;45:1565–1570. doi: 10.2135/cropsci2004.0463. [Cross Ref]
56. Potts SG, et al. Safeguarding pollinators and their values to human well-being. Nature. 2016;540:220–229. doi: 10.1038/nature20588. [PubMed] [Cross Ref]
57. Weiss, E.A. Sesame in Oilseed crops., 282–340 (Longman, 1983).
58. Andrade PB, Freitas BM, Macedo Rocha EE, Lima JA, Rufino LL. Floral biology and pollination requirements of sesame (Sesamum indicum L.) Acta Sci. Anim. Sci. 2014;36:93–99. doi: 10.4025/actascianimsci.v36i1.21838. [Cross Ref]
59. Free, J. B. Insect pollination of crops. (Academic Press, 1993).
60. Sarker AM. Effect of heoneybee pollination on the yield of rapeseed, mustard and sesame. Geobros (Jodhpur) 2004;31:49–51.
61. Delaplane, K. S. & Mayer, D. F. Crop pollination by bees, 215–219 (CABI, 2000).
62. Garibaldi LA, et al. Stability of pollination services decreases with isolation from natural areas despite honey bee visits. Ecol. Lett. 2011;14:1062–1072. doi: 10.1111/j.1461-0248.2011.01669.x. [PubMed] [Cross Ref]
63. Garibaldi LA, et al. Wild pollinators enhance fruit set of crops regardless of honey bee abundance. Science. 2013;339:1230200. doi: 10.1126/science.1230200. [PubMed] [Cross Ref]
64. Hoehn P, Tscharntke T, Tylianakis JM, Steffan-Dewenter I. Functional group diversity of bee pollinators increases crop yield. Proc. R. Soc. B. 2008;275:2283–2291. doi: 10.1098/rspb.2008.0405. [PMC free article] [PubMed] [Cross Ref]
65. Kleijn D, et al. Delivery of crop pollination services is an insufficient argument for wild pollinator conservation. Nat. Commun. 2015;6:7414. doi: 10.1038/ncomms8414. [PMC free article] [PubMed] [Cross Ref]
66. Potts SG, et al. Global pollinator declines: trends, impacts and drivers. Tree. 2010;25:345–353. [PubMed]
67. Winfree R, Bartomeus I, Cariveau DP. Native pollinators in anthropogenic habitats. Annu. Rev. Ecol. Evol. Syst. 2011;42:1–22. doi: 10.1146/annurev-ecolsys-102710-145042. [Cross Ref]
68. Settele J, Bishop J, Potts SG. Climate change impacts on pollination. Nat. Plants. 2016;2:1–3. [PubMed]
69. Melathopoulos AP, Cutler GC, Tyedmers P. Where is the value in valuing pollination ecosystem services to agriculture? Ecol. Econ. 2015;109:59–70. doi: 10.1016/j.ecolecon.2014.11.007. [Cross Ref]
70. Melathopoulos AP, Tyedmers P, Cutler GC. Contextualising pollination benefits: effect of insecticide and fungicide use on fruit set and weight from bee pollination in lowbush blueberry. Ann. Appl. Biol. 2014;165:387–394. doi: 10.1111/aab.12143. [Cross Ref]
71. OECD/FAO in: OECD-FAO Agricultural Outlook 2015. 10.1787/agr_outlook-2015-14-en Cotton (OECD Publishing, 2015).
72. Godfray HCJ, et al. Food security: the challenge of feeding 9 billion people. Science. 2010;327:812–818. doi: 10.1126/science.1185383. [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group