Closing the yield gap of soybean ( Glycine max (L.) Merril) in Southern Africa: a case of Malawi, Zambia, and Mozambique

300%. Although these technologies enhance soybean yields, they are not fully embraced by smallholder farmers. Hence, re ﬁ ning and bundling them in a digital advisory tool will enhance the availability of the correct information to smallholder farmers at the right time and improve soybean yields per unit area.


Introduction
Soybean (Glycine max (L.) Merril) is a legume crop that improves soil fertility (Agomoh et al., 2021), human and animal nutrition (Ali et al., 2020), and income (Tufa et al., 2019). Its demand continuously rises globally, while its production in Southern Africa is low, contributing to 1% of the global output (Murithi et al., 2016). This is against a regional demand of two million metric tonnes (TechnoServe, 2011), leading to significant net importation. Commercial production majorly contributes to 1% of the global output and dominates in the region, with South Africa producing about 65% of total regional produce, followed by Zambia, Malawi, and Zimbabwe (TechnoServe, 2011). Contribution by smallholder farmers (owning farms that are less or equal to 2 ha (AGRA, 2017)) is limited as the crop is marginally profitable due to low yields, in both quality and quantity, resulting from poor agronomic management practices, poor or low input use, poor soil fertility, inappropriate variety choices specific to agroecological zones, and erratic rainfall patterns and amounts (Van Vugt et al., 2017). Further exacerbating these low yields is the focus of most governments in the region on staple food crops like maize (Grant et al., 2012;Jakobsen and Westengen, 2022), leaving minimal support for the soybean value chain.
Soybean production in Sub-Saharan Africa (SSA) is increasing (Figure 1), a fact confirmed also by Khojely et al. (2018), yet this is not due to a yield increase per hectare but to increased land for production, i.e., the hectareage under soybean. Indeed, an increase in land area for soybean will dominate its future production (Foyer et al., 2019). Our analysis here indicated that in the past 10 years, from 2010 to 2020, soybean production increase was parallel to the harvested area, while yield declined in Zambia and Malawi (Figure 1). This was also reported by Meyer et al. (2018). Improved soybean market prices are the drivers of increased production with an increase in land area (Nsomba et al., 2022) due to food and feed demand, especially in the poultry industry (Ncube et al., 2017); the need to diversify from maize (Mhango et al., 2012;Renwick et al., 2021), particularly to combat climate change (Siamabele, 2021) and improve smallholder farmers' nutrition (Kolapo, 2011); improve soil fertility through biological nitrogen fixation (Ciampitti et al., 2021); break pest and disease cycles-pest control in integrated pest management such as fall armyworm and stemborer (Hailu et al., 2018;Matova et al., 2020) among others. Such land expansion for more production will gradually be unsustainable as the population grows; the land size per household reduces (Headey and Jayne, 2014;Muyanga and Jayne, 2014), and hence yield per unit area requires improvement (Mellor, 2014).
In 2020, FAOSTAT reported soybean yields in Malawi, Zambia, and Mozambique as 0.95, 1.41, and 1.67 t ha −1 , respectively (FAO, 2023). Comparing this with FAOSTAT Africa's average production in the same year indicates that soybean yield from Malawi is lower by 39.8%, while those of Zambia and Mozambique are above by 6.1% and 25.9%, respectively. This is an indication that the reasons behind the yield gap are country-specific and could further narrow to agroecology, households, and individual farms. In comparison to the global average yield, the yield gap widens in all three countries; the average soybean yield in the world was 195.2%, 99.0%, and 67.8% higher than in Malawi, Zambia, and Mozambique, respectively. This shows that globally, some countries are narrowing their yield gap, inching closer to the potential yield, while most countries in SSA lag. Identifying these countries that are narrowing the yield gap and their strategies and adapting them to Sub-Saharan Africa conditions and site-specifics is important. A further comparison of the FAOSTAT soybean yield in 2020 with the attainable soybean yield under breeding trials at the International Institute of Tropical (IITA) in those three countries indicated a yield gap of 276.3%, 117.1%, and 125.0% in Malawi, Zambia, and Mozambique, respectively. This is an even wider yield gap than a comparison with either Africa's or the global average yield. To accelerate increased production by smallholder farmers in the region and narrow the yield gap, research in development could play an essential role through the use of decision-support tools (Chivenge et al., 2022) and digital platforms. Ezui (2020) indicated a yield boost to smallholder farmers who used the AKILIMO App (Chernet and Pypers, 2020) in cassava production in Tanzania, Nigeria, and Rwanda. This shows the importance of digital platforms in enhancing crop yields. In soybean, such a digital tool does not exist for smallholder farmers in Sub-Saharan Africa. A tool that could address yield-limiting challenges in soybean and bottlenecks such as site-specific timing of planting, poor soil fertility/site-specific nutrient recommendation, suitable varieties, and disease/pest control (Van Vugt et al., 2017) and provide a wellcoordinated and informed agronomy advisory and market services. Furthermore, poor access to and use of improved agronomic strategies that address production constraints is a major challenge among the region's smallholder farmers. Approaches that take agronomy to scale and provide location-and context-specific decision support to farmers facing poor extension services could provide a viable route to enhancing soybean yields and production in the region. Moreover, improved soybean varieties and yield-enhancing agronomic practices developed and disseminated so far have had less adoption, and their scaling is required (Tufa et al., 2019). Therefore, the study: (i) assessed soybean production and yield gaps in Southern African countries-Malawi, Zambia, and Mozambique; (ii) determined site-specific yieldenhancing practices; and (iii) modeled through APSIM the potential soybean yields for Malawi, Zambia, and Mozambique under the right sowing date, and suitable variety. Later, these objectives will lead to the development of a soybean agronomic advisory tool that can assist smallholder farmers in getting the right information for suitable planting dates, crop variety choice, and site-specific fertilizer recommendations.

Literature search
The literature search, selection criteria of the articles, and metaanalysis followed the steps described by Hansen et al. (2022). The search engines used were Web of Science, Google Scholar, Google, and Scopus, and the treatments, i.e., management practices are described in Table 1. The search strings included the following keywords: (1) soybean: TOPIC: (soybean OR soybean crop) AND TOPIC: (grain yield) AND TOPIC: (Africa OR Sub-Saharan Africa OR Southern Africa OR Malawi OR Zambia OR Mozambique); (2) for nitrogen application: TOPIC: (nitrogen fertilization OR nitrogen) AND TOPIC: (soybean grain yield OR yield) AND TOPIC: (Africa OR Sub-Saharan Africa OR Southern Africa OR Malawi OR Zambia OR Mozambique); (3) for phosphorus application: TOPIC: (phosphorus fertilization OR phosphorus) AND TOPIC: (soybean grain yield OR yield) AND TOPIC: (Africa OR Sub-Saharan Africa OR Southern Africa OR Malawi OR Zambia OR Mozambique); (4) for potassium: TOPIC: (potassium fertilization OR potassium) AND TOPIC: (soybean grain yield OR yield) AND TOPIC: (Africa OR Sub-Saharan Africa OR Southern Africa OR Malawi OR Zambia OR Mozambique); (5) for inoculation: TOPIC: (biological nitrogen fixation OR BNF OR nitrogen fixation OR fixed nitrogen OR nitrogen derived from atmosphere (ndfa)) AND TOPIC: (soybean   Table S1). Data from these articles were cleaned, curated, and subjected to different analyses. Comparisons among the countries were done by finding the averages of production, harvested area, and yield; for Malawi and Zambia, it was from 2011 to 2020, whereas for Mozambique, it was from 2016 to 2020. Comparing the average yields between Malawi and Zambia was done from 2011 to 2020, while between Malawi and Mozambique, for both countries, it was conducted from 2016 to 2020.

Meta-analysis
A meta-analysis was conducted using Comprehensive Meta-Analysis V3 (Borenstein et al., 2003) on the effects of (1) P or K, (2) inoculants, (3) a combination of either P or K with inoculants, and (4) NPK and inoculants on soybean yield in Malawi, Zambia, Mozambique, and Sub-Saharan Africa. Analysis of the effect size of phosphorus application on soybean grain yield, heterogeneity, and publication bias tests were conducted. Furthermore, using SAS software (SAS, 2002), the means of the nitrogen, phosphorus, potassium, inoculant, and their combinations were compared using Fisher's least significant difference (LSD) method at a 95% confidence level. These yield data were also subjected to analysis of variance using the Proc Glm statement.

Crop modeling methodology Site
The modeling component of the study was based on the breeding trials that were undertaken at Chipata in Zambia, Chitedze in Malawi, and Angonia in Mozambique. The countries were selected because they have a significant proportion of smallholder soybean-growing farmers. Chipata district is in Eastern Province, Zambia, at an altitude of 1,070 M above sea level (m.a.s.l) at 13°70′ S and 32°61′ E. The soils are predominantly clay-loam textured, which are moderately acidic with low phosphorus content. The area receives 850-1,050 mm of rainfall per season (Thierfelder et al., 2013). Angonia (14°72′ S and 34°37′ E) is located within the north-western district of Mozambique within Tete province and at an altitude of over 1,000 m.a.s.l. The region receives at least 850 mm of rainfall per season, and the mean annual temperature ranges between 15°C and 23°C. The predominant soils are Lixisols and Luvisols but are nutrientdepleted (Nyagumbo et al., 2016;Inacio et al., 2022). Chitedze agriculture research station (13°98′ S and 33°67′ E) is at an altitude of 1,145 m.a.s.l. The predominant soils are ferruginous Latosols, which are deep and freely draining. The pH ranges from 4.5 to 6.0. Chitedze has unimodal rainfall ranging from 700 to 1,100 mm per season. The mean annual temperature ranges from 18°C to 22°C (Ivy et al., 2017).

Crop model
The APSIM Next-Generation model, known as APSIM NextGen , was calibrated based on the biophysical data from the breeding trials that were undertaken in Malawi, Mozambique, and Zambia. The APSIM NextGen was developed to improve execution speed, cross-platform development, model construction and visualization, and manager script flexibility (Holzworth et al., 2018). This improvement allows model developers to run larger simulations faster on multiple operating systems and simulate complex farming systems on temporal and spatial scales (Holzworth et al., 2018). It uses the APSIM Plant Modelling Framework (PMF) described by Brown et al. (2014), which allows the model developer to choose from a library of commonly used functions and algorithms for plant modeling. These are subsequently configured into a model description using the eXtensible Markup Language (XML). The APSIM Soybean model simulates crop development by pooling management and biophysical modules. The APSIM Soybean modules mimic crop growth upon input of certain data, for instance, crop management data, weather data, and soil data.

Model calibration
The APSIM Soybean model was calibrated based on the biophysical data collected from on-station IITA breeding trials. For effective and reliable simulations, the APSIM crop model needs to be calibrated for leaf area index (LAI), biomass, maturity, anthesis and emergence dates, and grain yield. In Sub-Saharan Africa, there are generally challenges due to inadequate data on soybean field research. The study, therefore, used one season of data but replicated across three sites in Malawi, Mozambique, and Zambia. The model was therefore calibrated based on one season's data. In this study, most of the parameters needed were not collected, so the APSIM crop model was calibrated only for maturity, anthesis, and grain yield. In addition, in breeding trials, each line only appears once per year per location; it is therefore not possible to have data for at least three seasons for each variety. The model was also calibrated with daily rainfall, solar radiation, minimum temperature, and maximum temperature data which were downloaded from NASA (Supplementary Table S2). The soil data were extracted from the ISRIC global soil database (Table 2) (Batjes, 2022). The datasets were used since there was no weather or soil data of satisfactory quality. Crop model performance was evaluated using the root mean square error (RMSE) (Moriasi et al., 2007), which compares measured and observed values of the different parameters. In this study, all the RMSE values for all the parameters, varieties, and locations were less than 30%, which is acceptable (Table 3). Specifically, RMSE values between 20% and 30% are considered "fair," 10% an 20% are considered "good," and 0% and 10% are considered "excellent" (Moriasi et al., 2007). The soybean yields were also validated with key informants.

Agronomic scenarios
The study coupled agronomic scenarios with the calibrated APSIM Soybean model. Soybean yield simulations were undertaken under different varieties, planting dates, and different phosphorus fertilizer application rates. The study evaluated long-duration varieties (grouped as taking > 120 days to physiological maturity) (Sc Serenade, Sc Sentinile, and Sc Safari) and short-duration varieties (days to physiological maturity is < 120) (MRI Dina). The attributes of these varieties are in Table 3. The study also evaluated the effect of different planting dates on yields from 1 November to 28 February, with a single planting every week in a similar fashion as done by Nyagumbo et al. (2017).

Soybean production in Malawi, Zambia, and Mozambique
Soybean production and harvested area in Malawi, Zambia, and Mozambique have been increasing from 2011 to 2020 (Figure 1). From 2011 to 2020, production in Malawi increased by 137.9%. This was parallel to harvested area, which increased by 150.5% while the yield decreased by 5.0%. In Zambia, the production and the harvested area increased by 154.7% and 252.2%, respectively. Conversely, the yield declined by 27.7%. In Mozambique, even though the data were only available from 2016 to 2020, there was an increase in production, harvested area, and yield per area by 59.6%, 18.4%, and 34.8%, respectively. However, the production and the harvested area were lower than that of Zambia and Malawi; both were highest in Zambia. Considering the available data from Mozambique from 2016 to 2020, despite the average yield declining in Zambia (1,583.3 kg ha −1 ), it was more than that of Mozambique (1,394.2 kg ha −1 ) and Malawi (972.6 kg ha −1 ). Malawi produced less per area, 85.5% and 43.4% lower than Zambia and Mozambique, respectively (Figure 1).

Yield gap
The actual yield of soybean reported by FAOSTATS for 2020 (FAO, 2023) and Solidaridad in 2020 and the attainable yield from a variety of evaluation trials across the three countries (Malawi, Zambia, and Mozambique) by IITA in 2020 were compared. The attainable yield in Malawi (3,565.3 kg ha −1 ) was 2.6 and 3.8 times more than Solidaridad's and FAOSTAT''s (FAO, 2023) yields, respectively ( Figure 2). In Zambia, the attainable yield (3,051.3 kg ha −1 ) was 2.5 and 2.2 times higher than Solidaridad's and FAOSTAT's (FAO, 2023) yields, respectively. Whereas in Mozambique, the attainable yield (3,750.7 kg ha −1 ) was 2.1 and 2.3 times more than the actual yields from Solidaridad's and FAOSTAT's (FAO, 2023), respectively ( Figure 2). The potential yield calculated by the APSIM model was 3.2% and 19.6% higher than attainable in Malawi and Zambia, respectively. Interestingly, it was 20.0% lower than the attainable yield in Mozambique (Figure 2). Comparing it with the actual FAOSTAT yields in Malawi, Zambia, and Mozambique, it was 288.4%, 195.2%, and 46.7% higher, respectively, and with Solidaridad's actual yields, it was 171.9%, 159.7%, and 42.4% more in Malawi, Zambia, and Mozambique, respectively (Figure 2).

Yield increase per unit area constraints Nutrients
The mean effect size of the study meta-analysis was 0.363 with a 95% confidence interval of 0.279 to 0.473, and the Z-value was −7.511 with p < 0.001. Both the Tau 2 (an estimate of the standard deviation of the distribution of true effect sizes) and I 2 (also known as Higgins I 2 , a statistical method to analyze heterogeneity and a The variety attributes are from Ngonga (2020) and Seed Co (2023).
measure for the proportion of observed variance that reflects real differences in effects (Borenstein et al., 2009)) were 0.332% and 0.000%, respectively (Figure 3). Phosphorus significantly increased the grain yield of soybeans ( Figure 4). Applying 22 kg ha −1 of phosphorus to soybean in Malawi increased yields by 52.4% ( Figure 4A), while in Zambia, an application of 20 kg ha −1 led to a 6.6% yield increase ( Figure 4B). In Mozambique, 40 kg ha −1 of P improved soybean yield by 42.6% ( Figure 4C). Overall, in Sub-Saharan Africa, the application of 21-30 kg ha −1 of P raised soybean yield by 48.2% ( Figure 4D). However, there was a negligible decline of 1% if more than 30 kg ha −1 of P.
Inoculating soybeans with rhizobia significantly improved their grain yield. By only inoculating soybean, the yield increased by 23.8% ( Figure 5A), whereas, coupling this with either 40 kg ha −1 of P or 60 kg ha −1 of K boosted the yields by 89.1% ( Figure 5B) and 26.0% ( Figure 5C), respectively. Furthermore, combining the inoculant with N, P, and K improved yield by 32.7% compared to not inoculating but applying P and K ( Figure 5D).

Planting dates and varieties
Simulations of different soybean varieties and planting dates impacted yields (Figures 6, 7). Variety Sc Serenade had the highest yields regardless of the planting dates over the 30 years across different locations and countries, as well as for the selected season of 2020/2021. However, the pattern was more pronounced in Angonia ( Figure 6A). This was similar to Chipata and Chitedze over the 30 years for the first half of the season, where from January till the end of the season there were no notable differences between all the four varieties (Figures 6, 8). The pattern was similar for the selected case study years, where there was a notable yield decline in the two locations after January (Figures 6, 8).
For Angonia, over the 30 years, the optimal planting window was remarkably wide from 12 December to 23 January, with Sc Serenade consistently being the highest-yielding variety ( Figure 6A). During the case study year, 2020-2021, the optimal planting date range was 2-23 January across all cultivars. However, for Sc Serenade and Sc Sentinile, the optimal planting date was 9 January, with the highest yield of 3,600 and 2,600 kg ha −1 , respectively. Also, the optimal planting date for Safari and MRI Dina was 16 January, achieving the highest yields of 2,600 had 2,500 kg ha −1 , respectively ( Figure 6B).

FIGURE 2
Actual soybean grain yield according to FAOSTAT 2020 and Solidaridad in 2020, the attainable yield from soybean variety breeding trials at the International Institute Tropical of Agriculture (IITA) in 2020, and the APSIM model's simulated potential yield.  At Chipata in Zambia, over the 30 years, the optimal planting period was from 1 November to 31 December, with the variety SC Serenade having high yields of 3,600-4,200 kg ha −1 . The yields of all four varieties began to decline, and the yield difference narrows after 2 January, with yields dropping to as low as 400 kg ha −1 when soybean was planted as late as 27 February ( Figure 8A). The pattern was similar at Chitedze in Malawi ( Figure 7A). For the season 2020/2021 at Good Nature in Zambia, the optimal planting window was from 28 November to 2 January, which led to the highest yields of 3,990 to 4,050 kg ha −1 . After 2 January, yields began to drop until the end of the season, reaching 1,200 kg ha −1 planted on 27 February ( Figure 8B). The pattern was similar in Malawi, Chitedze ( Figure 7B).

Discussion
The soybean yield gap in three of the Southern African countries covering an area locally referred to as the Chinyanja Triangle is wide, with an attainable yield of above 3.5 t ha −1 in Malawi and Mozambique, against an actual of 0.9 and 1.7 t ha −1 in those countries, respectively. In Zambia, the attainable is above 3.1 t ha −1 , against an actual of 1.4 t ha −1 . Also, the world's average soybean yield is 2.1 times that of Africa's and 2.8, 1.9, and 1.7 times that of Malawi, Zambia, and Mozambique, respectively (FAO, 2023). This huge difference between the attainable and actual yield, or the actual yields from various agroecologies and regions -the yield gap (Sadras et al., 2015) is attributed to a myriad of factors. These include the use of poor-performing varieties, limited or no fertilizer application, especially phosphorus, and a lack of or poor usage of inoculants in soils without a previous history of soybean production (Khojely et al., 2018). Interestingly, these factors are continuously being addressed through research ( Gbegbelegbe et al., 2019), from the release of high-yielding soybean varieties through breeding to proper crop, disease, pest, and soil fertility management, yet the yields of smallholder farmers are still low. The overall significant positive effect, as indicated by the metaanalysis, showed the importance of nutrient application in closing the yield gap in soybean. We revealed that inoculant application improves yields by 23.8%, application of between 20 and 30 kg P ha −1 increases yields by 48.2%, and a combination of P and K elevates the yields by 89.1%. These are also corroborated by some of the studies in Supplementary Table S1. Despite indications that these improved technologies increased soybean yield, most smallholder farmers are yet to adopt them (Tufa et al., 2019). However, there is hope, as early adopters lead the others in cementing that these technologies increase yields (Tufa et al., 2021). There are many reasons cited for low adoption, including lack or low ability to purchase inputs, lack or minimum skills and knowledge, variability in climate, poor access and availability of inputs, especially seeds, performance of the technology, poor soil fertility, gender and cultural practices (Odame et al., 2013), etc. Although poor extension services of the new improved technologies feature prominently too. Therefore, there is a need to improve extension services. In this era of mobile technology, e-extension is becoming popular, for example, the use of the Akilimo App by smallholder farmers producing cassava to enhance yields (Chernet and Pypers, 2020;Ezui, 2020). This can be developed for soybeans too to improve extension, farmers' skills, management, and production. However, technologies required for yield improvement and sustainability must be gathered, refined, and developed into a one-stop application tool that can be easily used by soybean smallholder farmers in Sub-Saharan Africa, and crop modeling is an essential component of that.
Evaluation of the contribution of different farm management systems to the total yield output can potentially inform farmers on the choice and potential of strategies that can be used to close the yield gap. This also enhances decision-making and resource allocation. The actual average soybean yields are relatively low. The modeling outputs showed that the use of different planting dates notably changed soybean crop yield from the current yield of an average of 1 t ha −1 to as high as 3.3 t ha −1 , which is 1.1 t ha −1 lower than the attainable yield. This shows the value of the right timing of planting. Previous studies in maize research showed at least a 5% yield decline for every week's delay in planting (Shumba et al., 1992;Nyagumbo et al., 2017). This pattern was also realized in this study, where there was at least 6% soybean loss per week's delay in sowing. Incorporation of other crop management practices such as suitable fertilizer quantity and timing, inoculant, and varieties could possibly further narrow the soybean yield gap. Factors such as fertilizer quantity and timing of planting have been identified as some of the major edaphic and socioeconomic yield gaps explaining factors in smallholder farming systems (Beza et al., 2017).
The use of different varieties significantly increased yields, thus providing further opportunities for narrowing the soybean yield gap. Most smallholder farmers use their own seed (sometimes called recycled seed), which reduces vigor, thus attaining yields of less than 1.5 t ha −1 . From the model simulations, the use of certified short-and long-duration varieties will lead to yield increases of up to 300% and 250%, respectively. The use of different varieties to narrow the soybean yield gap is, however, dependent on the right time of planting and location-specific, as the study showed that some varieties consistently performed well in high-rainfall areas such as Angonia, Mozambique. Angonia receives about 800-2,000 mm per season (Nyagumbo et al., 2016;Inacio et al., 2022). In environments such as Chitedze and Chipata in Zambia and Malawi, respectively, the effect of genotype diminishes more as rainfall starts to decrease and planting is late in the season. For example, at Chitedze (Malawi), which receives 700-1,100 mm per season (Ivy et al., 2017), SC Seranade and MRI Dina, whose yields are as high as 3,600-4,200 and 3,500-3,800 kg ha −1 , respectively, when planted early on 12 December, both drop to as low as 400 kg ha −1 when planted late on 27 February. Yield gap evaluation and yield decomposition, as in rice (Saito et al., 2017), wheat, and maize (Silva, 2017), are indeed important for soybean production in Sub-Sahara Africa to enhance targeted intervention and technology deployment.

Conclusion
Soybean production increased, although it was through an expansion in the cropped area rather than a yield increase per hectare. The average soybean yield in the world was 195.2%, 99.0%, and 67.8% higher than in Malawi, Zambia, and Mozambique, respectively. Also, the actual smallholder farmers' yield was 3.8, 2.2, and 2.3 times lower than the attainable yield in Malawi, Zambia, and Mozambique, respectively -a wide yield gap. Some site-specific yield-enhancing practices, including the application of (1) between 20 and 30 kg P ha −1 , (2) other nutrients like K and starter N, and (3) inoculants, improved soybean yields in Sub-Saharan Africa. Moreover, improved high-yielding soybean varieties preferred by most farmers in the region have been identified, and suitable management practices are being coupled with them to produce best-fit solutions. Furthermore, as shown by the APSIM model, suitable sowing dates have been identified for specific varieties in respective countries and agroecological zones (12 December to 23 January in Angonia and 1 November to 31 December in Chipata were the optimal sowing windows for SC Serenade). Altogether, proper management, soil fertility improvement strategies, and suitable sowing dates and varieties bundled in a digital tool that is easily available and usable by soybean smallholder farmers could enhance e-extension services and improve soybean yields per unit area.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.

Author contributions
JO, SM, JM, PM, and IN conceived and designed the manuscript. JO, SM, JM, PM, and IN collected, synthesized, analyzed, and presented the data. JO, SM, JM, and PM drafted the manuscript, which was reviewed by AC, MC, GC, and SK-B. JO edited and revised the manuscript's final version. All authors contributed to the article and approved the submitted version.

Funding
This study was funded under the OneCGIAR initiative, Excellence in Agronomy.