Introduction

Insect pollination is one of the imperative ecosystem services for more than 80% of known flowering plant species, many of which are part of our human diet1,2. These pollinators differentially improve crop quality/quantity and the availability of viable superior seeds3. Studies show that 75% of leading insect pollinator-dependent crops could account for 35% of global food production4. When the world was disaggregated by its economic regions, 59.8% of crops in developing regions were found to depend directly on these pollinators5. From a financial perspective, at the global scale, 51% (n = 46) of crops with a direct insect pollinator dependency accounted for 39% (625 billion euros) of the world's production value (€1618 trillion)6. The share from insect pollinators was estimated at €153 billion, while the agriculture sector's vulnerability to insect pollinator losses was about 9.8%6.

Food production constraints and economic losses due to declines of suitable insect pollinators could increase our dependence on already extensively consumed carbohydrate-dense staples (cereals, tubers, pulses), many of which are wind pollinated5,7,8,9. Insufficient consumption of a diverse number of nutrient-rich foods increases our vulnerability to non-communicable diseases10,11,12,13. Annual global deaths of up to 691,000 due to inadequate intake of insect-pollinator-dependent fruits, vegetables and nuts, which provide essential vitamins and oils, have been postulated from modelling studies14. Disease and death burdens would be higher in upper–middle–income countries14, while production losses would be higher in low-income countries, several of which are in Africa14,15.

More evidence is, however, still required to better establish the rate of species declines and extinction. This will allow scientists to fine-tune insect pollinator dependency studies and project high-risk crops. Nevertheless, a few global studies confirm that insect species, including taxa with insect pollinators like bees and butterflies, are declining faster than those increasing16. Model projections further depict more than 50% loss in suitable habitats globally for up to 49% of insect species and 44% of plant species if warming levels of 3.2 °C are attained17. Insect extinction risks of about 10% have also been predicted18, with the primary triggers being climate change, land use, pests, pathogens, and pollutants (pesticides and fertilizers)19. In Africa, most studies on insect declines are from farm-level studies on honey bees. Inputs like pesticides have been flagged as a threat to honey bee colonies in Rwanda20 and Kenya21. However, honeybees cannot account for the spatiotemporal behaviour of all species to stressors18.

As more tailored evidence on insect species decline is generated, global concerns and the existence of triggers that can affect plant-pollinator interactions19,22 are sufficient to drive preliminary impact assessment studies on crop-insect pollinator dependency. Results from the latter studies could provide feedback to impel insect pollinator monitoring actions. Unfortunately, in Africa and for several local crop varieties, insights on the percentage dependency of crops on insect pollinators are still lagging. While this remains a research gap, data can be inferred from available publications4 to support downstream studies. This has been the case in an increasing number of countries outside Africa, such as India23, Belgium24, Latin America25, China26, Pakistan27, and on a global scale4,5. However, only some studies have focused on African countries like Benin28, Burkina Faso29, Ethiopia30, Morroco31, Malawi, Nigeria, Tanzania, and Uganda10. Results show a high dependence of 46% (Malawi) to 76% (Tanzania).

To further exploit existing datasets and to expand knowledge on the dependency of the agriculture sector on insect pollinators in Africa, we elected to work in Rwanda, where a comprehensive study is yet to be undertaken. This is despite the country's high dependency on agriculture and the presence of crops like passion fruits, macadamias, mangoes, avocadoes, tomatoes, beans, eggplants, soybeans, coffee, groundnuts, etc., that can benefit from insect pollinators32. These crops are also essential for food security, health, nutrition and economic prosperity in Rwanda33,34. For instance, export revenues of up to $465.4M (2018/2019) have been generated by the sector, with marked contributions from coffee (~ $69M; insect-pollinated) and tea (~ $88M; wind-pollinated)4,33. More than half of the sector comprises household farmers32,35. In the Huye district of Rwanda, such farmers could lose up to 11% of farm income and 20% of marketed production share in the complete absence of insect pollinators11. Yet, like several African countries, Rwanda's Strategic Plan for Agriculture Transformation (PSTA4) 2018-202437 does not mention a crop-insect pollinator strategy. This, therefore, limits actions on and investments in insect pollination as an agronomic input.

This work thus aims to elucidate if the government of Rwanda should include insect pollinators as an agronomic input in its subsequent strategic plan for agriculture transformation. We hypothesized that insect pollinators are important for Rwanda's agriculture and would bring a valuable monetary contribution to the economy. Our study had the objective of mapping the insect pollinator dependency of 71 crops that are cultivated in Rwanda, evaluating the total economic value of crop production and that due to insect pollination for 32 representative crops using the bio-economic approach and FAO crop data for the year 2014 and 2020, and finally, evaluating whether production from 21 representative crops could meet consumption demands before and after the complete loss of insect pollinators. We found evidence for Rwanda to take strategic and urgent actions on insect pollination in its crop production subsector, including stringent insect pollinator conservation practices. Our findings also strengthen the need for other countries to commence such assessments and take action.

Methodology

In this study, we assumed that all pollinators behave similarly (completely absent or completely present). This is irrespective of which crops they pollinate or how many pollinators can pollinate a given crop. Secondly, we used the mean pollinator dependency range/ratio of crops and not their actual pollinator dependency values. Data manipulation, analysis and visualization steps were performed using the R platform and suitable packages (Table_5_SuppInfo)36.

Data extraction

Crop price and agriculture variables (yield, production and area harvested) data for the period 2010–2020 were downloaded from the Food and Agricultural Organization website37 in 2021 and then updated in February 2023, just before the completion of this study. These were downloaded as single files and then merged into one data set. Additional crops mentioned in the grey and published literature were used to enrich this data set. This included crops captured in the Agricultural Household Survey by the National Institute of Statistics of Rwanda32, the National Agricultural Export Development Board annual report33 and the Baseline Report on Rwanda Horticulture38. We extracted distinct crops from this data set, which resulted in a total of 71 representative crops. This data set was used to map the insect pollinator dependency of these crops. Crop type information was extracted from the Indicative Crop Classification (I.C.C.) list39.

Crop-insect pollinator dependency mapping

The crop-insect pollinator dependency classification model by Klein et al.4 was used to classify the above crops. The insect pollinator dependency depicts the extent to which a plant depends on insect pollination for improved agronomic outputs (yields) or the production quantity the plant would lose in the complete absence of suitable insect pollinators. These dependencies are represented by lower and upper limit values between 0 and 100%, referred to as the insect pollinator dependency ratio, rate or range of a crop.

The model classifies plants into three main groups: direct, indirect, and none insect pollinator-dependent. Based on the above values, direct insect pollinator-dependent crops are further divided into four subgroups: essential (> 90–100%), great (> 40–90%), modest (> 10–40%), and little (> 0–10%). Indirect insect pollinator-dependent crops are divided into two subgroups: crops whose seed yields increase in the presence of insect pollinators and we consume their resultant vegetative parts (increase seed production, ISP), and crops whose seed production increases only during breeding, but we do not consume their resulting vegetative parts (increase breeding, IB). The latter crops reproduce through vegetative propagation. None insect pollinator-dependent crops do not require insect pollinators (wind-pollinated). Crops like tree tomatoes, stevia, moringa, and leeks with an unknown insect pollinator dependency were excluded.

The above classification model was further confirmed against country-specific crop dependency data sets from Brazil40, Argentina41, India23 and Burkina Faso29. Where discrepancies occurred, we used information from the most recent publications or the African-specific context (Table_1_SuppInfo). Stenchly et al.29 used a different classification for the indirect dependent group. They classified carrots, taro, onions and cabbages as directly dependent on insect pollinators, given that we consume their vegetative parts, whose production is influenced by seed quality/quantity achieved during insect pollination42,43.

The insect pollinator dependency table was generated after assigning an insect pollinator dependency (mean ratio, subgroup, group), crop name, species name, and crop type to each of the 71 crops in our data set (Table_1_SuppInfo). Crops with a different dependency value defined by Stenchly et al.29 were marked with an asterisk (*).

Evaluation of the crops' distribution

The percentage distribution of a given pollinator-dependent group, subgroup or crop type was computed using data in the crop-insect pollinator-dependency table (Table_1_SuppInfo). Firstly, we aggregated (subaggregated) crops based on their pollinator-dependency groups (direct, indirect, none), subgroups (essential, great, modest, little, IB, ISP, none) or crop types (fruits, vegetables, nuts/oilseeds, stimulants, pulses, tubers, cereals, sugars, spices/herbs). For a given number of crops in a subaggregate (nsubaggregate) and the total number of crops in the aggregate (naggregate), the percentage distribution of the subaggregate is computed as in Eq. (1).

$${\text{Percentage distribution}}\;\left( \% \right) = 100\left( {n_{subaggregate} /n_{aggregate} } \right)$$
(1)

For instance, to compute the percentage distribution of the direct insect pollinator-dependent group, we divide the number of direct insect pollinator-dependent crops (nsubaggregate) by the number of crops in the main aggregate/entire data set (naggregate) and multiply by 100.

Determination of the economic value of insect pollinators to crop production

The bioeconomic approach was applied to determine the economic/monetary value of insect pollinators to crop production6. It takes into account three input parameters: total crop production quantity (Q), farm gate price (P) and the insect pollinator dependency value (D). Hence, for a given crop (i), where i Є [1: I], the economic value (EV) is the product of the price (Pi) and production quantity (Qi) as in Eq. (2). The economic value due to insect pollinators (IPEV) for a crop (i) is the product of its EV and the mean of its dependency ratio (Di) as in Eq. (3). The final results are expressed in the currency of the price.

$${\text{Economic value of crop production }}\left( {{\text{EV}}} \right) = \mathop \sum \limits_{i = 1}^{I} \left( {P_{i} *Q_{i} } \right)$$
(2)
$${\text{Economic value of insect pollination }}\left( {IPEV} \right) = \mathop \sum \limits_{i = 1}^{I} \left( {P_{i} *Q_{i} *D_{i} } \right)$$
(3)

\({Q}_{i}\) = total agricultural production (tons) of a crop (i), \({P}_{i}\)($/tons) = the corresponding market value of the crop, \({D}_{i}\)= the crop's pollination dependence rate (%), which is represented by its mean dependency ratio. This includes essential crops (0.95), great (0.65), modest (0.25), and little (0.05), indirect and none pollinator-dependent crops (0.00) (Table_1_SuppInfo).

To calculate the EVs, we downloaded the FAO data set on the value of agricultural production from 2014 to 2020 as a single file. It contains crop data and the calculated economic value (U.S. dollars; $; Eq. 2) for several crops. We combined this with the data set on crop price, yield, production, and area harvested (Sect. 2.1) to easily compare different variables and obtain additional crop names, price and production data. Where applicable, we used the price and production information to calculate the EV of crops as outlined in Eq. (2). Crops lacking price and production data, or an EV, were eventually excluded, resulting in a final data set of 32 crops.

Using data in our crop insect pollinator-dependency table (Table_1_SuppInfo), we assigned an insect pollinator dependency (mean ratio, group and subgroup) primarily from Klein et al.4 as a new data set or from Stenchly et al.29, as a separate data set. We added the species name and crop type to these crops before calculating the IPEV (Eq. 3). All downstream steps were applied to both data sets. However, for visualization, the detailed outputs from the data set of the former authors are often displayed and described in the results.

For visualization, we aggregated crops based on their insect pollinator dependency group and crop type before computing their total EV, IPEV, area harvested and production for each year. Then, we displayed the results for 2014 and 2020 (Table_2_SuppInfo). The annual outputs mean and standard deviation for the entire study period (2014–2020) were also displayed based on their insect pollinator dependency groups when the classification model of Klein et al.4 (Table_3_SuppInfo) and Stenchly et al.29 (Table_4_SuppInfo) was applied.

Determination of the vulnerability of the agricultural sector to insect pollinators

We calculated the vulnerability of the crop production subsector (CVR) if faced with a complete pollinator loss6 by dividing the IPEV sum by the EV sum for any crop category. The result is multiplied by 100 and expressed as a percentage, as in Eq. (4). In our study, the IPEV for indirect and wind-pollinated crops is zero when Di = 0. Hence, this vulnerability solely depends on direct insect-pollinated crops when Di > 0.

$${\text{Crop Vulnerability Ratio}}\;\left( {\% {\text{CVR}}} \right) = \frac{IPEV}{{EV}} = \frac{{\mathop \sum \nolimits_{i = 1}^{I} \left( {P_{i} *Q_{i} *D_{i} } \right)}}{{\mathop \sum \nolimits_{i = 1}^{I} \left( {P_{i} *Q_{i} } \right)}}*100$$
(4)

We computed the vulnerability ratio for each crop type and the annual vulnerabilities using the total economic values generated in Sect. 2.4 (Tables 2, 3, 4_SuppInfo).

Determination of overproduction based on consumption demands

To understand whether production can meet consumption demands without insect pollinators, we used the food balance sheet data from FAO for Rwanda for 2014–20206, which contains consumption and production data. Crops and their products (e.g. maize and its products) were first aligned to the main crop (maize). Processed foods (e.g. infant food), crop types (e.g. cereals and others) or miscellaneous were excluded as they do not represent a single crop. This resulted in a final data set of 21 crops. These crops were assigned an insect pollinator dependency (mean ratio, group and subgroup), species name and crop type as defined in our crop-insect pollinator dependency table (Table_1_SuppInfo).

For overproduction before pollination loss for a specific crop (i), where i Є [1: I], its consumption (Ci) is subtracted from its total production quantity (Qi) and then divided by its consumption as in Eq. (5). For overproduction after pollination loss, we obtained the proportion of production that is not dependent on insect pollinators (1-Di), multiplied this proportion by its total production quantity (Qi), subtracted its consumption (Ci) and divided the result by its consumption as in Eq. (6). These outputs were expressed as a percentage of consumption.

$${\text{Overproduction before pollination loss }}\left( {\% consumption} \right) = \frac{{\mathop \sum \nolimits_{i = 1}^{I} \left( {Q_{i} - C_{i} } \right)}}{{\mathop \sum \nolimits_{i = 1}^{I} \left( {C_{i} } \right)}}$$
(5)
$${\text{Overproduction after pollination loss }}\left( {\% consumption} \right) = \frac{{\mathop \sum \nolimits_{i = 1}^{I} \left( {Q_{i} \left( {1 - D_{i} } \right) - C_{i} } \right)}}{{\mathop \sum \nolimits_{i = 1}^{I} \left( {C_{i} } \right)}}$$
(6)

\({\mathrm{Q}}_{\mathrm{i}}\) = total agriculture production (tons) of a given crop i, \({\mathrm{C}}_{\mathrm{i}}\) (tons) = consumption of crop i, \({\mathrm{D}}_{\mathrm{i}}\)= the pollinator dependence rate of crop i (%), \((1-{\mathrm{D}}_{\mathrm{i}})\)= portion of crop i without insect pollination.

Determination of relative changes

Relative changes between two years of interest were calculated using the method described by Azien et al.5 and information in Table 2_SuppInfo. For this, one of the years was set as a reference year. Therefore, to compute the relative change in a variable (x) between two years, the value of a variable (x) in the reference year \((x_{{t_{reference} }} )\) is subtracted from that of the year of interest \((x_{t} )\) before dividing the output by the value of variable \((x)\) of the reference year \((x_{{t_{reference} }} )\) as in Eq. (7). In our study, 2014 was set as the reference year. Our variables included yield, production, area harvested, production value, economic value and consumption demand before or after the complete loss of insect pollination for any given crop. The final outputs are expressed as percentages.

$${\text{Relative percent change}}\;\left( {\% {\text{RC}}} \right) = \frac{{ (x_{{\dot{t}}} - x_{{t_{reference} }} )100}}{{x_{{t_{reference} }} }}$$
(7)

With \(x_{{\dot{t}}}\) = value of a variable \((x)\) for a given crop in year t, \(x_{treference}\) = value of variable \((x)\) for a given crop in the reference year.

Missing data imputation

The Multivariate Imputation by Chained Equations (MICE) package in R was used to fill in missing data44. MICE has the advantage over other methods of filling in missing values many times based on the observed values of a given sample and the relationship in the data set45. The MICE imputation was used to fill the production value of green beans (2014 and 2015) and green chillies/peppers (2014).

Statistical analysis

Statistical differences in production value between the three insect pollinator-dependent groups were calculated using the non-parametric Kruskal–Wallis test for each year. This is recommended for non-normally distributed data, especially if the researcher chooses not to manipulate the data through log transformation or removing outliers, for instance46. This test was selected over the analysis of variance (ANOVA) since the data for specific years (2014, 2015, 2016, 2020) arguably satisfied the test for normality, even though they satisfied the homogeneity of variance assumption. Secondly, not all crop types are represented in each insect pollinator-dependent group—this narrowed the number of independent variables that could be included.

The Kruskal–Wallis test was performed using its built-in function in RStudio. We defined two hypotheses, null and alternative, and then used the probability (p) value to statistically interpret the results47. For the null hypothesis, H0 implies that the means of the dependent variable (production value) is not significantly different (p > 0.05) among the independent variables (pollinator-dependent groups). Otherwise, we accept the alternative hypothesis H1 if at least one group mean differs significantly (p < 0.05) from the others. Lastly, the dunn_test function in RStudio was used to obtain information about which two combinations of the insect pollinator-dependent subgroups were statistically different.

Results

Dependence of edible crops on insect pollinators

Of the edible crops grown in Rwanda, 62% directly depend on insect pollinators

We found that 62% (n = 44) of representative edible crops (n = 71) that are cultivated in Rwanda are directly dependent on insect pollinators (Fig. 1, table). Of these, the majority (14 crops) have a modest dependency, 12 crops a little dependency, 11 crops a great dependency and seven crops an essential dependency (Table_1_SuppInfo). Crops that did not require pollinators constituted the second largest group (21.1%, n = 15), followed by 16.9% (12 crops) with an indirect dependency. Among the latter group, the majority (n = 9) are pollinated by insects during crop breeding, whereas for the others (n = 4), insect pollination is vital for increased seed production.

Figure 1
figure 1

Dependence of 71 representative edible crops cultivated in Rwanda on insect pollinators. The crop-insect pollinator dependency ratio was inferred mainly from Klein et al.4 and other country-specific publications (Table_1_SuppInfo). The table in the chart depicts crops in each insect pollinator dependency group (direct, indirect, and none). The bar graph and associated digits represent the percentage distribution of each crop type.

Fruits, vegetables, and oilseeds substantially contribute to the direct insect pollinator dependency

The most representative crops were fruits, vegetables, nuts/oilseeds, and spices/herbs. The first three contributed the most to the overall dependence of the direct insect pollinator group. All represented nuts/oils have a direct dependency, while up to 77% (n = 17) of fruits and 58.3% (n = 7) of vegetables also have a direct dependency (Fig. 1, white bars). However, most nuts/oilseeds (75%) are at the lower end of the dependency grid of insect pollinators (modest – little), in contrast to 58.8% (n = 10) of fruits that occupy the higher end of the dependency grid (great-essential). For vegetables, there is a comparable distribution among the insect pollinator-dependent groups.

Among the least represented crop types, a majority have a direct dependency. We found that direct insect pollination is vital for 75% (n = 6) of spices/herbs. Of the six pulses, 83.4% occupy the lower end of the insect pollinator dependency grid. Of the two stimulants, there was an equal share, with coffee having a modest dependency and tea having no dependency. Tubers populated the indirect insect pollinator group, predominantly for breeding purposes. All representative sugar crops and cereals are wind-pollinated (Table_1_SuppInfo).

Economic value of crop production and insect pollination

The total economic value of crop production increased in Rwanda from $2.551 billion (2014) to 2.788 billion (2020)

The total economic value of crop production from the 32 representative edible crops in Rwanda was in 2014 estimated at $ 2.551 billion (Table 1), with a production of 6.6 million tons across an estimated area of 2.1 million ha (Table_2_SuppInfo). Of this, the direct insect pollinator-dependent group had the highest number of representative crops (43.8%). This accounted for 20.4% ($521.5 million) of the economic value. Furthermore, its production (14.1% of the total) was slightly lower than that of the none-dependent group (15.7%). This former group also occupied the second-largest harvested area (31.8%).

Table 1 Total economic value of production (EV) in Rwanda in the year 2014 and 2020, share from insect pollination (IPEV), percentage relative change (RC) in the year 2020 against the reference year 2014, and vulnerability ratio (%VR).

The indirect insect pollinator-dependent group, which had 34.4% of representative crops, accounted for the highest total economic value (64.4%; $1.643 billion), production (70.2%), and harvested area (46.0%). The wind-pollinated group had the lowest number of representative crops (21.9%), made the lowest contribution to the economic value (15.1% = $386.2 million) and further occupied the smallest harvested area (22.2%).

The total economic value of crop production in 2020 increased by 9.3%, equating to $2.788 billion (Table 1). This was the second-highest economic value attained between 2014 and 2020 from representative crops (Table_3_SuppInfo). There was also a 33.1% increase in relative production (8.8 million tons) and a 5.8% increase in harvested area (2.2 million ha), despite the pressures imposed by COVID-19. Additionally, the economic value of the direct insect pollinator-dependent group decreased slightly by 3.4%, resulting in an estimated value of 503.5 million. Nevertheless, this was still higher than that of the none-dependent group ($414.6 million; relative increase = 7.3%) but lower than that of the indirect dependent group, which still made the highest contribution of 58.5% ($1.869 billion; relative increase = 13.8%).

The direct insect pollinator group had a 4.5% increase in production (966,356 tons) and a 25% increase in its area harvested (826,475 ha) in 2020. Its production remained slightly less than that of the non-insect pollinator-dependent group (1 million tons) and far less than that of the indirect-dependent group (Table_2_SuppInfo). This latter group accounted for up to 77.5% (6.8 million tons) of the overall production, even though its area harvested had declined by 15% (811,063 ha).

At the crop type level, indirect insect pollinator-dependent fruits and tubers contributed the most to the total economic value (Table 1). The respective amounts were $879.1 million and $737.2 million in 2014, before decreasing slightly in fruits (0.8%) but increasing in 2020 by 27.8% in tubers. Cereals, pulses (direct) and vegetables (direct) followed suit, with 130–354 million.

The economic value of crops due to insect pollination declined by 35.8% in 2020

The share of the total economic value of production due to insect pollination was in 2014 at an estimated $186.2 million (Table 1). This translated into a 7.3% vulnerability of the crop production subsector to insect pollination. However, this share declined in 2020 by 35.8% (Table 1), with a reduced vulnerability to insect pollinators of 4.3%. Stimulants (76.6%), fruits (71.9%) and vegetables (38.3%), represented mainly by coffee, mangoes, and pumpkins, contributed to this decline (Figure_2_SuppInfo).

Amongst the crop types, vegetables (direct) accounted for the highest economic value due to insect pollination in both study years (Table 1). This value was in 2014, $153 million, before declining by 38.3% ($94.5 million) in 2020. They were also the most vulnerable crop type in 2014 (78.2%) and 2020 (68.7%). The pulses had an estimated value of 12.2–15.2 million. Unlike vegetables, pulses were the least vulnerable to insect pollinators (5%). Fruits, estimated at $10.8 million in 2014, were the second most vulnerable crop type (38.7%). They also experienced one of the highest declines (71.9%) in 2020. The stimulant, coffee, was the fourth highest contributor in 2014 ($5.4 million) but experienced a 76.6% decline in 2020, resulting in nuts/oilseeds taking this spot. Stimulants also had the same vulnerability (25%) as spices, which made the least contribution in both years (Table 1).

When we recalculated the IPEV using the classification model of Stenchly et al.29, the total monetary value was estimated at $219.91 million in 2014 and increased to $235.24 in 2020. The sector's vulnerability was estimated at 8.62% in 2014 and 8.44% in 2020 (Table 4_SuppInfo).

The insect pollinator-dependent group had the lowest production per hectare of harvested area but the highest production value per ton

We also observed that in terms of production quantity per hectare of harvested area, the direct insect pollinator-dependent crops were producing, on average, the lowest quantity per hectare (977,737.69 tons/733,687.5 ha = 1.33 tons/ha) over the 7-year (2014–2020) period (Table 3_SuppInfo). The wind-pollinated group had a slightly higher production per hectare of harvested area (934,986.5 tons /514,047.43 ha = 1.82 tons/ha), while the indirect dependency group generated the highest quantity (5,509,515.06 tons/875,269.86 ha = 6.29 tons/ha).

The direct insect pollinator-dependent group had the highest mean production value per ton in 2014 and 2020. In 2014, this value was 769.9 \(\pm\) 91.3$/ton (± standard error) and was significantly greater (p = 0.025) than that of the indirect dependent group (366.2 \(\pm\) 51.5$/ton) but not the none-dependent group (584.0 \(\pm\) 207.5 $/ton; p = 0.18). This value decreased in 2020 by 13.0% and 11.5% in the direct and indirect insect pollinator-dependent groups, respectively, but increased slightly (5.2%) in the none-dependent group (Fig. 2). The value in 2020 remained significantly different (p = 0.010) only between the direct and indirect insect pollinator-dependent groups. Dry peas (none) had the highest production value/ton in both years, followed by groundnuts, coffee, lemon, oranges, mangoes (direct dependencies) and rice (none) in at least one of these years.

Figure 2
figure 2

Production value per ton ($/ton) of representative crops cultivated in Rwanda in 2014 and 2020. The values are represented with respect to their insect pollinator-dependent group (direct, indirect, and none). Grey dots represent the mean production value per ton in each group. Kruskal–Wallis and Dunn tests were used to calculate significant differences between groups using their functions in RStudio.

Production and consumption demands

Insect pollinator-dependent crops could face major constraints in meeting consumption demands

Using data from 2014 to 2020, we found that among the eight representative direct insect pollinator-dependent crops, only coffee, and in selected study years, chillies, groundnuts and dry beans could meet consumption demands (Figure_2_SuppInfo). In groundnuts, this was apparent during the later years of the study (2017–2020), while in chillies, this was between 2015 and 2019. In contrast, oranges, soybeans, lemon/limes, and tomatoes had differential deficits ranging from 30 to 70% during the study period. The latter was affected after the loss of insect pollinators, whereas the others were affected under both situations. On average, over the seven years, oranges and soybeans showed deficits during both pollination situations, while the others had marked deficits (tomatoes and oranges) or declines (chillies and dry beans) only after pollinator loss (Table 2).

Table 2 Mean overproduction from 2014 to 2020 before and after insect pollinator loss (%consumption).

Of the eight indirect insect pollinator-dependent crops, only potatoes, onions, and cassava could not meet consumption demands (Figure_3_SuppInfo and Table 2). While this was marginal in onions and cassava (± 10%), it was pronounced in potatoes during the early study years (2014–2018), with deficits of at least 40%, before recovering to exceed consumption demands (10%) during the later years.

Wheat and rice were the two wind-pollinated crops that could not meet consumption demands (Table 2), with deficits of at least 30% (Figure_4_SuppInfo). Although dry peas initially had a slight deficit (~ − 10%), they recovered in the subsequent years, remaining at marginal levels (~ 10–15%). Maize and sorghum, on the other hand, remained in excess of consumption demands throughout the study period, with excesses of about 30–100% observed.

Discussion

Our study has provided for the first time a more advanced insight into the dependence of the crop production subsector in Rwanda on insect pollinators, the economic value of crop production and the impact on consumption demands before and after the complete loss of suitable insect pollinators. What is unique in our study is the multiyear approach applied to observe fluctuations, if any, using the most recent data.

In our study, we found support for our hypothesis that insect pollinators are important to Rwanda's agriculture and will bring in a valuable monetary contribution. This was despite the absence of the added economic value of 30 insect-pollinated crops (e.g. passion fruits, macadamia, sunflowers, etc.) with no price data6. Losing this monetary contribution has repercussions on farmers' revenue11 and GDP48.

A loss of more than 60% of insect-pollinated crops in Rwanda due to a complete absence of these pollinators would affect household farmers who cultivate fruits (with 67.3% growing avocadoes and 32.7% mangoes), vegetables (with 43.7% growing amaranth, 25.6% tomato, 24.2% eggplant), pulses (88% growing beans), nuts/oilseed (12.3% growing soybean, 6% growing groundnuts) and stimulants (355, 771 farmers)32,49. Production losses will compound food scarcity, shift consumption demands to indirect- or wind-pollinated staples (maize, rice, cassava, beans, etc.), and heightened food insecurity, especially in resource-limited households50.

Even though food choices are not guided by pollination, farming households in Rwanda (80%) rely heavily on their farm produce for consumption10,32. This includes insect-pollinated produce, especially vitamin A-rich vegetables (pumpkin, carrot, squash, amaranth/dodo), fruits (passion fruit, avocado, papaya), but also nuts/oilseed (groundnuts), pulses (dried beans), spices/herbs and coffee34,50. Additionally, their consumption, which increases with availability, can be enhanced by promoting kitchen gardens and capacity-building programs on food and nutrition13,51. These actions align with the government of Rwanda's strategic plan to sustain a nutrition-sensitive agriculture sector which increases dietary diversity and the consumption of nutrient-dense foods34,35. This continues to be important in combating the 35% food insecurity risk, 33% stunting, 45% mortality risk (especially among children under five years), 26% obesity in women, and inadequate consumption of diverse crop types in Rwanda34,52,53.

A focused intervention should, therefore, be on insect-pollinated vegetables and fruits, which, in our study, were the most vulnerable to insect pollinator losses and the highest contributors to the vulnerability of the production subsector. The latter vulnerability has remained slightly lower than that for world agriculture4 but is higher in given years than for East Africa6. Together with nuts/oilseeds (soybean and groundnuts), these vegetables (tomatoes) and fruits (oranges, lemon) are among the crops that struggled to meet consumption demands in the 7-year study period.

Low-vulnerability crops (nuts/oilseeds, pulses, coffee) should also be given attention due to their added nutritional value (oils, proteins, vitamins) and international traction (coffee)14. To close nutrition gaps, biofortification is used to augment the nutrient composition of consumer-preferred crops (e.g. dry beans for iron delivery) in Rwanda9 and globally54. However, if insect pollination is vital to these biofortified crops, appropriate interventions should be prescribed to achieve superior yields. For instance, dry beans can further benefit from yield increases per hectare of up to 46% during open pollination by insects55.

Our findings are consistent with postulations that low-income countries are yet to reap the full benefits of insect pollination5,14,15. Our representative direct insect pollinator-dependent group generated the least production per hectare of area harvested. Given that agricultural land is a limited resource37, Rwanda might struggle to compete in production quantity. However, they can capitalize on improving quality through strategic insect pollinator interventions as an added agronomic input. This will better position their produce in the international market.

Conclusion and recommendations

Given the potential value-addition of crop-insect pollination to food systems, nutrition and health, Rwanda should incorporate relevant insect pollination actions in its strategic plans.

A first recommendation is for Rwanda to increase efforts on beehive establishment, mainly where pollinator-dependent crops are cultivated. This will directly benefit several vital crops like coffee, sunflower, macadamia, mango, avocado, etc., that depend on bees and other pollinators4.

Secondly, policies guiding bee hive establishment should align with those regulating pollutive fertilizers, herbicides and pesticides56. While these are strategic farm inputs in Rwanda35, their negative consequences on bees and other pollinators are a concern20. Insect-friendly farm inputs (organic fertilizers) should be promoted, while pesticides should be carefully analyzed for their impacts on insects and discontinued if harmful.

Thirdly, Rwanda should nurture natural ecosystems, which are pollinator-rich habitats57,58. Insect-pollinated crops could be cultivated close to such vicinities to facilitate interactions. Additionally, Rwanda's land consolidation action35 should promote the establishment of patches of fields with high pollinator attractants (e.g. sunflower or flowering plants with bright colours) close to areas with pollinator-dependent crops58,59.

Above all, advocacy and capacity-building activities should be intensified35. Additional emphasis should be on the adequate handling of approved farm inputs, techniques to mitigate the impacts of harmful farm practices, the importance of pollinators in food systems, monitoring and documenting insect visitors, hand pollination techniques, and climate change mitigation actions (e.g. afforestation)60.

Challenges and future perspectives

This study has also identified gaps, mainly around the availability of data. Research efforts should be geared towards characterizing the insect pollinator dependency of crop cultivars in Rwanda, their nutritional value, and possible health benefits. This will open the gateway for impact-driven crop pollinator studies. A second research effort should focus on monitoring insect pollinator species in different ecosystems to clarify concerns on whether, how and at what rate insect pollinators are declining. Thirdly, national surveys and data collection efforts should capture crop production, area cultivated, area harvested, yield, farm gate price, and consumption demand data more consistently for an increasing number of crops, including minority crops in Rwanda. This data should also be disaggregated per crop and not aggregated under a crop type. These findings and recommendations could contribute to Rwanda's poverty alleviation and food security agenda, and could be translated in other African countries.