- 1Plant Genetics and Rhizosphere Processes Laboratory, TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech, Gembloux, Belgium
- 2Department of Crop Science, School of Agriculture, University of Rwanda, Musanze, Rwanda
- 3Department of Research, Rwanda Agriculture and Animal Resources Development Board, Huye, Rwanda
- 4Integrated and Urban Plant Pathology Laboratory, TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech, Gembloux, Belgium
- 5Departement des Sciences Agronomiques, Université Catholique de Bukavu, Bukavu, République Démocratique du Congo
- 6Unité Santé des Plantes et Forêts, Centre wallon de Recherches agronomiques, Gembloux, Belgium
- 7Numerical Genetics, Genomics and Modeling Laboratory, TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech, Gembloux, Belgium
- 8Phytopathology Laboratory, University of Louvain, Louvain-la-Neuve, Belgium
- 9Department of Public Health, School of Health Sciences, Mount Kenya University, Kigali, Rwanda
- 10Modelisation and Development, Gembloux Agro-Bio Tech, University of Liège, Gembloux, Belgium
- 11Agroecology Lab, Université Libre de Bruxelles (ULB), Brussels, Belgium
- 12Tropical Crop Improvement Laboratory, Department of Biosystems, Katholieke Universiteit Leuven, Heverlee, Belgium
Cassava (Manihot esculenta Crantz) is a vital crop in Rwanda where it ranks as the third most consumed staple. However, cassava productivity remains below its yield potential due to several constraints, including important viral diseases, such as cassava brown streak disease (CBSD). Because various factors can be addressed to mitigate the impact of viral diseases, it is essential to identify routes of virus contamination in the cassava agrosystems from the seed system to farmer's practices and knowledge. The present study aimed at (1) assessing the current cassava seed system and farmers' practices and their knowledge of the biotic constraints to cassava production, (2) determining the status of CBSD as well as critical factors associated with its spread through the seed system channels, and (3) determining factors that influence cassava productivity in Rwanda. A cross-sectional study was carried out from May to September 2019 in 13 districts of Rwanda. A total of 130 farmers and cassava fields were visited, and the incidence and severity of CBSD were evaluated. CBSD was detected in all cassava-producing districts. The highest field incidence of CBSD was recorded in the Nyanza district (62%; 95% CI = 56–67%) followed by the Bugesera district (60%; 95% CI = 54–65%), which recorded the highest severity score of 3.0 ± 0.6. RT-PCR revealed the presence of CBSD at the rate of 35.3%. Ugandan cassava brown streak virus was predominant (21.5%) although cassava brown streak virus was 4% and mixed infection was 10%. An informal cassava seed system was dominant among individual farmers, whereas most cooperatives used quality seeds. Cassava production was found to be significantly influenced by the use of fertilizer, size of the land, farming system, cassava viral disease, and type of cassava varieties grown (p < 0.001). Disease management measures were practiced by a half of participants only. Factors found to be significantly associated with CBSD infection (p < 0.05) were the source of cuttings, proximity to borders, age of cassava, and knowledge of CBSD transmission and management.
Introduction
Cassava (Manihot esculenta Cranzt) ranks as the sixth most important food crop worldwide and the fourth after rice, maize, and wheat among developing and emerging countries (Otekunrin and Sawicka, 2019; Saranraj et al., 2019). In Rwanda, cassava is the third most important crop after banana and sweet potato (Night et al., 2011). Because of its importance in several tropical regions and its relatively good performance on marginal lands under suboptimal climatic conditions (Burns et al., 2010), cassava is recognized as a subsistence crop to overcome food insecurity for the fast-growing population in areas prone to important climatic changes (El-Sharkawy, 2004; Chavez et al., 2005; Lobell et al., 2008; Burns et al., 2010). Although cassava plays an important role as a food security crop in sub-Saharan Africa, it is also used as a cash crop in various cassava-growing regions (Spencer and Ezedinma, 2017; Munganyinka et al., 2018).
The yield potential of cassava under optimum conditions is about 90 tons of fresh roots per hectare, which is equivalent to 30 tons of cassava dry matter per hectare (El-Sharkawy, 2004). More than half (61%) of cassava production is taking place in sub-Saharan Africa; however, cassava yield in tropical countries is still far below its production potential. Indeed, in 2017, the world cassava yield was about 11.08 tons of fresh roots per hectare, and the top cassava producer (Nigeria) had an average yield of 8.75 tons per hectare, followed by the Democratic Republic of Congo with 8.14 tons per hectare (FAO, 2019; Otekunrin and Sawicka, 2019). Cassava production in Rwanda varied between 3,000 and 3,701 Mt of fresh roots per year from 2015 to 2018 with a reported average yield of about 14.5 tons per hectare (FAO, 2018; Rwanda Agricultural Board1). Despite its resilience under adverse environmental conditions, the production of cassava remains constrained by several abiotic and biotic factors. The former includes postharvest deterioration, infertile soils, planting unimproved traditional varieties, and inadequate farming practices, whereas the latter includes green mites, mealy bug, cassava bacterial blight, and viral diseases (Bull et al., 2011; Kombate et al., 2017).
As a consequence of viral diseases and the lack of resistant varieties, cassava yields have drastically decreased in many countries. Cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) are the most economically important cassava diseases, causing yield losses of more than US$1 billion yearly globally (Legg and Thresh, 2003; Legg et al., 2006; Patil et al., 2015; Rey and Vanderschuren, 2017).
CBSD has so far only been reported in sub-Saharan Africa. CBSD is particularly devastating because it negatively impacts cassava tuberous roots both quantitatively and qualitatively, causing important economic losses to African farmers (Mohammed et al., 2012). For decades, CMD has been managed through dissemination of resistant varieties, but unfortunately, the distributed CMD-resistant varieties were found to be sensitive to CBSD in Rwanda and in many other African countries (Legg et al., 2001; Thresh and Cooter, 2005; Bua, 2017; Nyirahorana et al., 2017). CBSD is caused by two species of single-stranded RNA viruses of the family Potyviridae, Genus Ipomovirus; Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) (Mbanzibwa et al., 2009). In Rwanda, CBSD was first reported in 2009 in the Southern province (Muhanga district). CBSD has since spread to reach most cassava-producing regions in the south (Kamonyi, Ruhango, Nyanza, Gisagara and Muhanga Districts) and east of the country (Bugesera, Nyagatare, Gatsibo, Kirehe, Kayonza and Ngoma Districts) (Nyirahorana et al., 2017). A study conducted in Rwanda in 2014 reported a distribution of CBSD incidence as follows: 74.2% UCBSV, 15.3% CBSV, and 10.5% mixed infection (Munganyinka et al., 2018).
The CBSVs are transmitted by either the whitefly Bemisia tabaci and/or exchange of infected planting materials between farmers. Plant pathologists and extension services recognize the importance of establishing a disease-free seed system to mitigate the spread of CBSD (Mbanzibwa et al., 2009; Patil et al., 2015; Maruthi et al., 2017). In Rwanda, formal distribution of clean planting seeds usually involves the whole production chain from the Rwanda Agricultural and Animal Resources Development Board (RAB) where researchers produce basic clean seeds. Basics seeds are then distributed to seed multipliers across different regions for further multiplication, and before their dissemination to farmers, a quality seed certification agency is involved to ensure the quality of the planting materials (Broek and van den Byakweli, 2014; Andrade-Piedra et al., 2016). Quality seed refers to the seed preferred by farmers and consumers with good health (virus free), genetic purity, appropriate physiological age, and physical quality (Andrade-Piedra et al., 2016).
The concomitant existence of an informal seed system involves production and dissemination of seeds by farmers usually without quality seed certification, which tends to favor disease dissemination and maintenance in cassava agrosystems (Muthoni and Nyamongo, 2008). Because of its prominence in most cassava agrosystems worldwide, the informal seed system plays a major role in the rapid spread of diseases (Andrade-Piedra et al., 2016; Maruthi et al., 2017).
Various management measures are applied to reduce CBSD impact, including investigation and breeding of CBSD-tolerant varieties (Masumba et al., 2017; Anjanappa et al., 2018; Sheat et al., 2019; Shirima et al., 2020) as well as the dissemination of disease-free planting material to farmers (Alicai et al., 2016). Conversely, farmers who reuse cuttings from their own fields will not escape the disease as this tends to maintain 30–50% of infection, especially in CBSD hot spots (Rwegasira and Rey, 2011; Patil et al., 2015). Unfortunately, the seed system sustainability remains fragile and needs to be consolidated to provide appropriate healthy planting cassava material to farmers. Despite the emergence of CBSD in cassava fields in Rwanda, there has been limited information about the CBSD distribution and factors associated with its transmission as well as factors that affect cassava productivity in the country.
The present study aimed at assessing the impact of farmers' practices and their knowledge of the biotic constraints on cassava production as well as determining the status of CBSD and the critical factors associated with its spread through the seed system channels.
Materials and Methods
Study Area
The study was conducted in 13 cassava-growing districts of Rwanda during 2019. Districts in the southern and eastern provinces are considered to be major cassava-producing areas. In the south, five districts were surveyed, namely, Gisagara, Nyanza, Ruhango, Muhanga, and Kamonyi, whereas in the east, six districts were surveyed, namely, Bugesera, Nyagatare, Kayonza, Gatsibo, Kirehe, and Ngoma. In addition, two districts from the western and Northern provinces, Nyamasheke and Gakenke, respectively, were included in the study (Figure 1).
Figure 1. Map of Rwanda showing surveyed districts. The red dots represent the location of the cassava fields assessed in the study.
Farmer and Field Selection
A multistage sampling method was applied to select cassava farmers and fields. In the first stage, 13 districts, representing both major and minor cassava-growing areas, were selected. In the second stage, according to the information provided by district agronomists and RAB, five main cassava-growing sectors were purposely selected within each district. Sectors in which cassava is marginally cropped were not kept in the selection because they may not be representative of current limitations encountered by cassava growers. In the third stage, as cassava farmers are classified into two categories (either individual farmers or belonging to cooperatives) (Miklyaev et al., 2021), two farmers per sector, one individual farmer and one farmer belonging to a cooperative, were selected from a sampling frame provided by sector agronomists using a simple random-sampling approach, making a total of 10 farmers per district and 130 interviewed farmers for the 13 districts surveyed. Selected participants were always heads of household, either a man or woman, depending on their availability. In case both were available, simple random sampling was applied to select one of them.
Furthermore, for each participant, a field with cassava plants older than 6 months was also visited for disease evaluation. The distance between two cassava fields visited was around 10 km. Within the selected fields, 30 plants were selected for leaf and stem CBSD symptom examination including five plants of the two diagonals and five of the four sides (Rwegasira et al., 2011). The field incidence per district was recorded as the percentage of symptomatic plants out of the total examined. The 10 plants examined at the two diagonals were further pooled and used for CBSVs indexing by RT-PCR.
Farmer Interviews
Primary data used in the study were collected using a structured questionnaire (Supplementary Datasheet 1), semistructured interviews and observations on five key subject areas (such as demography, social economics, agronomy, seed accessibility, and availability factors and disease aspect) relevant to cassava production. For farmers, quality seed was defined as seed certified by a seed-quality inspection agency (Rwanda Agriculture Board) that may be obtained from a professional seed multiplier. District and sector agriculture extension officers liaising with local community leaders were involved in the mobilization of farmers. Permission to conduct research in the area was sought from the administration of the study area (district and sector agronomists) through official communication by RAB authorities. Participants were told the purpose of the research and that participation was voluntary. Oral consent was given before starting the interview and field visit. All records were identified by study identification number to keep participant privacy and confidentiality.
Cassava productivity among participants was considered as cassava yield (kg fresh matter per hectare) in the present study, and income generated from cassava was estimated without considering price variation between districts. The price of cassava was estimated as the average of the cassava root prices recorded during the survey, which was 90 Rwandan francs (Rwf) per kilogram (Max = 95 Rwf, Min = 85 Rwf, std = 4.3 Rwf). Therefore, the income was estimated as the yield times the estimated price of cassava. To learn the factors influencing cassava productivity, the benefit generated from cassava per hectare was calculated as the total income minus the cost of production. Variables collected during the survey are listed in Supplementary Table 1. Biotic variables observed and measured in the visited fields are listed in Supplementary Table 2.
Disease Severity Assessment
A 1–5 CBSD symptom scale (Gondwe et al., 2003) was used to measure the degree of severity of CBSD aerial symptoms in the fields. The scale used was 1 = no apparent symptoms; 2 = slight leaf feathery chlorosis with no stem lesions; 3 = pronounced leaf feathery chlorosis, mild stem lesions, and no dieback; 4 = severe leaf feathery chlorosis, severe stem lesions, and no dieback; and 5 = defoliation, severe stem lesions, and dieback (Gondwe et al., 2003). The average degree of severity was calculated omitting the score of 1, which represent asymptomatic plants to provide a true picture of the severity in the fields assessed (Sseruwagi et al., 2004). An average of disease severity per district was calculated based on the observation of 30 × 10 = 300 plants.
Sampling Test Materials for RT-PCR
Samples were collected from May to September 2019. In each cassava field assessed, 10 cuttings from 10 plants examined along the two diagonals were collected per field and established in the screenhouse. In total 1,300 plants were grown in the screenhouse from 130 fields. Five-month-old plants from collected cuttings were used for molecular analysis. In the laboratory, 20 leaf samples resulting from growth of cuttings collected on the same field were pooled (lower and middle leaves were used per plant); thus, a total of 130 samples were tested for CBSVs using RT-PCR.
Molecular Analysis
RNA Extraction
Total RNA was extracted from ~0.2 g cassava leaf using the cetyltrimethylammonium bromide protocol previously described (Abarshi et al., 2010).
cDNA Synthesis and RT-PCR
Synthesis of cDNA was done using a ProtoScript II Reverse Transcriptase kit (BioLabs, UK) following the manufacturer's instructions. Briefly, a Master Mix containing 1 μl d(T)23 (50 mM), 2 μl of buffer, 1 μl of 0.1 mM DTT, 0.5 μl Protoscript II RT, 0.5 dNTP Mix, and 3 μl of nuclease-free water was prepared. Then, 2 μl of RNA template was subsequently added, making 10 μl per reaction. The reaction mixture was incubated in a PCR thermocycler at 42°C for 1 h for primer annealing and cDNA synthesis, followed by 20 min at 65°C for inactivation of the ProtoScript II Reverse Transcriptase. The resulting cDNA samples were stored at −20°C.
The synthesized cDNA was subjected to polymerase chain reaction using a Taq G2 Hot Start Master Mix from Promega. The primer pair F: 5′-CCTCCATCWCATGCTATAGACA-3′ and R: 5′-GGATATGGAGAAAGRKCTCC-3′ that amplifies ~703 bp of CBSV and ~800 bp of UCBSV isolates was used (Elegba, 2018). The 10 μl PCR reaction contained 5 μl G2 Mix, 0.4 μl each primer (0.4 μM final concentration), 1 μl cDNA and the volume was brought to 10 μl with nuclease-free water. PCR conditions were as follows: predenaturation at 95°C for 2 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 30″ and elongation at 72°C for 50″ and final elongation at 72°C for 5 min.
An internal control gene from cassava called Manihot esculenta protein phosphatase 2 A (MePP2A) was detected in parallel using a pair of primers F: 5′-TGCAAGGCTCACACTTTCATC-3′ and R: 5′-CTGAGCGTAAAGCAGGGAAG-3′ that amplifies 150 bp of MePP2A to ensure the accuracy of the PCR results by ruling out any false negative results (Moreno et al., 2011). PCR amplification was checked by loading 10 μl of PCR products in 1% (w/v) agarose gel stained with Gel red in 1X Tris-acetate-EDTA (TAE) buffer for 1 h at 200 V to allow the separation of amplicons from the two isolates. The PCR products were visualized on a UV gel documentation system and photographed.
Data Analysis
Raw data were transferred into Statistical Package for Social Sciences (SPSS) version 22 for analysis. Frequency and proportions for categorical variables were computed to describe the basic attributes of the respondents (farmers) as well as the occurrence of cassava infection (defined as the presence of CBSV or UCBSV after PCR) in sampled fields. A bivariate analysis with chi square tests was used to determine factors associated with cassava virus infections (categorical variable). Then, multivariate logistic regression analyses were performed by considering together all significant factors during bivariate analysis. The goodness of fit was assessed using Hosmer–Lemeshow test. In all statistical tests, differences were considered statistically significant at p < 0.05.
Similarly, descriptive statistics using mean and standard deviation for continuous variables were used. After testing for normality, cassava benefit was not normally distributed, and Sqrt transformation was used to ensure the normality of the cassava benefit distribution. Then, an ANOVA test was used to compare means to estimate the effect of sociodemographic, economic, and agronomic variables on cassava productivity. Multiple linear regression analyses were performed to assess the independent factors affecting cassava productivity. Ten factors were considered together, and upon fitting all against cassava productivity using multiple linear regression and specifying a “stepwise” method, five factors remained independently affecting cassava productivity. To validate the regression analysis, the histogram and P-P plot presented in Supplementary Figure 3 were used to show the reasonable normality of data. Furthermore, multi-collinearity problems were assessed using tolerance and variance inflation factors (VIFs) as well as Durbin Watson was used to check for autocorrelation that the residuals from linear regressions are independent.
Results
Sociodemographic Characteristics of Participants, Source of Cassava Cuttings, Trust, and Disease Management
Among the 130 farmers interviewed, a majority of them were individuals (58.5%), men (64.6%), and married people (66.9%). The percentages of men and women in the participants also reflects the gender distribution for cash-crop growers in Rwanda (Gender Monitoring Office, 2017; Munganyinka et al., 2018). Most of the participants (72.3%) had attended primary school. The respondents' age was in three categories with 45 years and above being slightly higher (36.2%) than the other two categories (Supplementary Table 3). The majority of participants (62.1%) reported that they obtain planting material from their own fields and use the same materials over many seasons. A minority of interviewed farmers (25.1%) acquired planting materials from seed multipliers every season. Among those, a large proportion (86%) were cooperative members (Figures 2A,B).
Figure 2. Sources of cassava planting materials, trust, and disease management among farmers in Rwanda during 2019. (A) Sources of planting materials, (B) categories of farmers vs. sources of seeds, (C) level of trust in seed quality among participants, (D) disease management methods applied by participants (NGO, non-government organization).
Although most farmers use seeds from their fields, about 45% questioned their quality and feared that their cassava fields may succumb to diseases. A gap in cassava viral disease management was noted through the study as 50.3% of interviewed farmers took no action to control viral disease and only 19.5% declared roguing out infected plants from their field (Figures 2C,D).
CBSD Incidence and Severity Score of Observed Symptoms
Field incidence of CBSD and severity of aerial symptoms were evaluated in the 13 districts. The highest incidence (62%; 95% CI = 56–67%) was recorded in the Nyanza district, and the lowest (12%; 95% CI = 8–16%) was observed in Gakenke district. Disease mean severity scores varied from district to district, ranging from 3.0 ± 0.6 recorded in Bugesera to 2 ± 0.2 recorded in Nyamasheke. Table 1 shows field incidence and mean severity score of observed CBSD symptoms in the examined fields.
Table 1. Field incidence, severity, and frequency of cassava plants showing aerial CBSD symptoms in Rwanda, 2019.
Molecular diagnostics was performed on a total of 130 samples collected from 13 districts. Samples were analyzed by RT-PCR for detection of CBSV and UCBSV. The overall incidence of CBSVs was 35.3%. Among the positive samples, 61% (28/46) were UCBSV whereas 11% (5/46) were CBSV, and 28% (13/46) had mixed infection of both CBSV and UCBSV (Figure 3), Supplementary Figure 1 shows RT-PCR detection of CBSVs in field samples.
Figure 3. Incidence of CBSD based on RT-PCR in different districts of Rwanda, 2019. The highest incidence was recorded from Nyanza, followed by Gisagara and Bugesera districts. UCBSV, Ugandan cassava brown streak virus; CBSV, Cassava brown streak virus.
All 13 districts surveyed were found to be affected by CBSVs based on RT-PCR results, and the highest incidence (60%) was recorded in the south, in Nyanza district followed by Gisagara and Bugesera districts, both displaying an incidence of 50%. A single infection of UCBSV was found in all districts except Muhanga, whereas a single infection of CBSV occurred in Muhanga, Bugesera, Gatsibo, and Gakenke. Mixed infections were recorded in most districts except in Ngoma, Gatsibo, Gakenke, and Nyamasheke (Figure 3).
Bivariate Analysis of Factors Associated With CBSD Incidence
Fields From Farmers Working in Cooperatives Display Lower CBSD Infection Rate
Using data collected during the survey, we tested whether sociodemographic characteristics influence the level of cassava infection. Bivariate analysis (using Chi square test) of sociodemographic factors of 130 farmers (for whom cassava fields were visited) stratified by CBSV infection revealed that there is a significant association between category of respondents and cassava infection, and individual farmers had more infected fields than farmers in cooperatives (p = 0.023) (Supplementary Table 3). Farmer's age also significantly influenced the level of cassava infection of their field (p = 0.043). All the other sociodemographic factors had no influence on levels of cassava infection.
Fields Established With Planting Material From Seed Multipliers Have a Lower Probability to Be CBSD Infected
We further investigated the link between cassava seed accessibility and CBSV infection. Our analysis showed that the source of cassava cuttings and proximity to the border had a significant impact on cassava infection with p values of 0.001 and 0.021, respectively (Table 2). Farmers who used seeds from their own field were more likely to have infected fields (60.9%) than those who got seeds from seed multipliers (10.9%). It was also noted that farmers near the country's border had more infected fields (54.3%) (Table 2).
A majority of cassava fields surveyed were more than 8 months old (56.2%). Most participants grew improved varieties (65.4%) and had access to extension services (83.1%). In fact, the analysis revealed that there was a significant association between age of the plant and cassava infection (p < 0.001) where the plants aged <8 months (57.1%) were significantly more likely to be healthy than to be infected (19.5%). Furthermore, a significant association was noted between use of fertilizers and cassava infection (p = 0.002) as the farmers using fertilizers (68%) were significantly more likely to have healthy than infected fields (39%) (Table 3).
Farmers' Awareness of Cassava Viral Diseases Is Associated With Lower CBSD Incidence
Although all farmers were aware of at least one cassava viral disease's existence, 34.6% of them did not know the symptoms of cassava viral diseases. Symptoms of CMD were easily recognized by 31.5% followed by 20.8% who recognized both CMD and CBSD (Supplementary Figure 2 shows symptoms of CMD and CBSD). Farmers who were not aware that the viruses can be transmitted had more infected fields (65.2%) (p < 0.001), and likewise, those who did not know disease management had more infections in their fields (67.4%) compared with those who knew the management techniques (p < 0.001) (Supplementary Table 4).
Multivariate Analysis of Risk Factors Associated With CBSD
Eight factors that showed significant association (p < 0.05) during bivariate analysis (including source of cassava cuttings, proximity to the border, age of the plants, use of fertilizers, category of respondents, age of farmers, knowledge of cassava viral disease transmission, and knowledge of cassava disease management) were considered together in a multivariable analysis to identify the variables associated with cassava infections. Upon fitting the factors using multiple logistic regression and specifying the “backward conditional” method with removal at p < 0.05, five factors remained in the final analysis as shown in Table 4. After testing the goodness of fit using the Hosmer–Lemeshow test, the chi-square value was 4.80 with degree of freedom 6, and p value was 0.570, which indicates that the fitted model was adequate.
Farmers who use cuttings from their own fields or from other farmers' fields had a more than seven-fold higher risk than those who used cuttings from a seed multiplier (p < 0.05). Respondents located near the border had a four times higher risk of having CBSVs than those located far away (p < 0.05). Cassava plants <8 months old had fewer risks of infection compared with the older ones (p < 0.05). Similarly, those who were not aware of the disease transmission and management had nearly a three times higher risk of having the infected plants (p < 0.05) (Table 4).
Effect of Sociodemographic Variables on Cassava Productivity
Although diseases are among the main constraints to cassava yield and production (FAO, 2015; Rey and Vanderschuren, 2017), socioeconomic factors might also be associated with suboptimal cassava yields. Taking advantage of the farmer survey, we performed an ANOVA to assess the effect of different sociodemographic variables on cassava productivity. The demographic variables that significantly affected cassava productivity were district, level of education and farmer category. The district of Ruhango had a significantly higher average of cassava benefit per hectare than other districts (p < 0.001). The mean of cassava benefit was significantly higher among respondents with a secondary level of education compared with illiterate farmers (p < 0.001). Similarly, cooperative farmers had significantly more average cassava benefit than individual farmers (p < 0.001) (Supplementary Table 5).
Effect of Economic and Agronomic Variables on Cassava Productivity
During bivariate analysis using ANOVA, all agronomic and economic factors tested showed that they significantly influence cassava productivity (p < 0.001). The mean cassava benefit was significantly higher among farmers with fields displaying less viral disease (after PCR), larger size of land, growing improved varieties, practicing monoculture, applying fertilizers, and using seeds from seed multipliers (Table 5).
Multiple Linear Regression for Cassava Productivity
A multiple linear method was used to determine which variables had the most significant impact on cassava productivity. The results show that 5 out of 10 factors predict cassava productivity. The 10 variables considered in the model were all those variables significant at the bivariate analysis using the ANOVA test (see Table 5; Supplementary Table 5). Upon fitting these 10 factors against the dependent variables using multiple linear regression and specifying the stepwise method, five variables (using fertilizer, size of the land, farming system, cassava viral disease, and type of cassava varieties grown) remained affecting the productivity of cassava (Table 6). There were no collinearity issues found in this study between the different outcome and independent variables as tolerance was above 0.1 and VIFs were below 10. Durbin Watson was also used to check for autocorrelation that the residuals from linear regressions are independent. Durbin Watson with zero indicated positive autocorrelation, and four indicated negative autocorrelation while around two indicated that the residuals are uncorrelated.
Using fertilizer was the main positive effect in the first place and explained 57% (R2 = 0.57) of the changes in cassava benefit. In the second model, size of the land is added, which led to 67% (R2 = 0.67) of the changes in the cassava benefit. In the third model, farming system is added, which led to 69% (R2 = 0.69). In the fourth model, viral disease led to 71% (R2 = 0.71), and type of cassava varieties grown in the fifth model led to 72% (R2 = 0.72) variation in cassava benefit (Table 6).
Discussion
The present study used a comprehensive cross-country survey to assess the current cassava seed system, farmers' practices, and their knowledge of the biotic constraints, the status of CBSD, and critical factors associated with its spread throughout the cassava seed system as well as factors that influence cassava productivity in Rwanda.
The current findings confirmed the occurrence of CBSD (both CBSV and UCBSV) in Rwanda. The disease was found in all 13 districts surveyed, indicating that it has spread out in all major cassava-growing regions, including Kirehe and Nyagatare, where CBSVs were not detected in previous studies (Munganyinka et al., 2018).
In our survey, we found that districts located near the border displayed a higher rate of CBSD incidence. The highest field incidences and severities were recorded in the three districts, namely, Nyanza, Bugesera, and Gisagara, bordering Burundi. In an earlier study carried out by Munganyinka et al. (2014), Nyanza and Gisagara districts also displayed the highest CBSD incidence, which confirms them as hot spots for CBSD (Munganyinka et al., 2018). This observation might correspond to the informal movement of cassava cuttings across countries that leads to the importation of infected cuttings or use of genetic material that is more susceptible to CBSVs. Furthermore, the high CBSD incidence in those districts could be because since its first report in 2009, the virus could have flourished in those areas season after season due to the relatively warm environments that favor proliferation of whitefly vectors (Campo et al., 2011).
A CBSD survey performed in Burundi previously reported an average incidence and severity of 15.3 and 2.3%, respectively (Bigirimana et al., 2011). A decade ago, UCBSV was the only viral species that was associated with the disease in Burundi although it was already present in Tanzania (Rwegasira et al., 2011).
Based on RT-PCR diagnostics, the overall CBSV incidence was found to be less than CBSV incidences reported elsewhere in East-Central Africa (Kenya, Tanzania, Malawi, Zambia) (Rwegasira et al., 2011; Mbewe et al., 2014; Koima et al., 2018; Mulenga et al., 2018). This difference might be due to a later introduction of CBSVs in Rwanda. UCBSV was prevalent across the country, indicating that it is the commonest cause of CBSD in the country. Similar findings were reported in the survey conducted in 2014 among major cassava-growing regions (Munganyinka et al., 2018). It should be noted that our study revealed an increased rate of mixed infections from 10% in 2014 to 28% in 2019, highlighting the dynamism of disease spread over time either through vectors or exchange of unhealthy cuttings among farmers.
Despite farmers being aware of the importance of quality seed to optimize output, it was observed that most farmers have difficulty identifying quality seed as was reported earlier (Minot et al., 2007).
Our present study indicates that a majority of farmers (76.9%) use informal ways to get cuttings for free from either their own fields or from neighbors. Because the supplied planting material often suffers from a lack of quality control, farmers are more likely to plant virus-infected cuttings, which can lead to low yield and reduced benefits. The informal seed system might be further maintained by the lack of knowledge about cassava viral diseases as identified in the survey, highlighting the need to increase farmers' awareness of using quality seed and also to mobilize the private sector to invest in the commercial cassava seed business. Previous studies have already highlighted the need to promote farmers' awareness (Chipeta et al., 2016; Nduwumuremyi et al., 2016; Bentley et al., 2017) as farmers using cassava planting materials from appropriate sources (research institutions, NGOs, etc.) appear to have fields with reduced CBSD infection (Gondwe et al., 2003).
Despite differences in CBSV incidence between districts, farmers who used quality seeds kept having lower CBSV infection compared with those reusing seeds from their own field or from other farmers. This observation also highlights the importance of the human factor (transport and exchange of unhealthy cuttings), contributing to the propagation and dissemination of CBSD (Patil et al., 2015; Maruthi et al., 2017).
The ultimate goal of cassava farming is to optimize yields for food security and income. Taking into account that cassava is a climate-resilient crop essential in fighting food insecurity for the fast-growing population (El-Sharkawy, 2004; Chavez et al., 2005; Burns et al., 2010), crop-intensification programs have been launched in Rwanda and other African countries to improve agricultural production (FAO, 2016). Despite that cassava production is a profitable investment in Rwanda (Gasangwa, 2013), its production still remains constrained by several factors.
The present study indicates that high yield and income generation from cassava logically increase with the use of fertilizers, the size of the land, monoculture, and improved cassava varieties, whereas incidence of viral diseases is associated with a decrease in the income. Although the use of fertilizers is known to increase yields of cassava storage roots (Munyahali et al., 2017), cassava continues to be seen as a resilient crop with limited or no requirement for fertilization. Therefore, promoting adequate use of fertilizers in cassava should remain a priority in future crop-intensification programs.
The survey also indicates that cooperatives nearly double their benefits as compared with individual farmers. This observation could be explained by the fact that most cooperatives exploit bigger land and better comply with good agriculture practices (including the use of fertilizer, improved quality seeds, and practicing monoculture) than individuals. A recent assessment of the benefits associated with cultivating improved high-yielding varieties shows that yields can be increased more than four-fold as compared with the cultivation of local varieties in Africa (Khonje et al., 2015). Thus, farmers should be encouraged to group into cooperatives and to practice good agriculture practices (GAP) to meet the high yield and food security objectives of crop intensification programs.
Following the emergence of CBSD in 2009, efforts made by governmental entities to combat cassava viral diseases, have helped reduce their impact on cassava production. However, CBSD remains a major challenge to cassava production despite past and ongoing efforts to breed for virus-resistant varieties, to distribute clean planting material, and to promote GAP (Catholic Relief Services, 2012). Although time and resource constraints prevented us from performing a more extensive survey, the random approach operated in the cassava-growing sectors makes our survey representative of the cassava-growing areas.
Our study reveals that viral diseases remain a constraint to cassava productivity with a disease prevalence that has increased to 35.3% in Rwanda. Therefore, there is a need to continue efforts to introgress virus-resistance traits into farmer-preferred varieties and to establish a cassava seed system enabling sustainable and affordable supply of clean planting material to farmers. The strengthening of the cassava seed system also requires the development of important capacities for virus diagnostics (Mukasa, 2015; Wossen et al., 2020). To sustain the implementation of GAP, there is also a need to increase farmers' awareness of cassava diseases and the immediate benefit of using quality seeds.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
Ethical review and approval was not required for this study with human participants, in accordance with the local legislation and institutional requirements.
Author Contributions
HV, CBr, CBu, NG, CN, and JB contributed to the conception and design of the study. CN, JB, AN, VB, SM, CT, and HV contributed to the methodology design. HV, SM, CBr, CT validated the study. MH, YB, CT, CN, and JB performed the statistical analysis. CN and JB wrote the first original draft. YK, YB, CT, LL, SM, EM, AN, and HV revised and edited the manuscript. HV provided resources and supervision. All authors contributed to the article and approved the submitted version.
Funding
The study was carried out as part of the ARES (Académie de Recherche et d'Enseignement Supérieur, Fédération Wallonie-Bruxelles, Belgium) funded project iCARE (Improved Cassava Virus Resistance mitigation strategies and development of a disease-free seed system).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fsufs.2021.699655/full#supplementary-material
Supplementary Figure 1. Detection of CBSV and UCBSV by reverse transcription-PCR in samples collected from (A) Nyanza and (B) Bugesera districts. Expected PCR product sizes were CBSV (703 bp) and UCBSV (800 bp). From left to right, 100 bp DNA ladder; C+ correspond to positive control for both CBSV and UCBSV, whereas C– corresponds to negative control; lanes 1–10 represent pool of samples from fields 1–10; faster DNA ladder. (C) Detection of the internal control MePP2A for samples in (A); (D) detection of the internal control MePP2A for samples in (B); expected PCR product sizes was 150 bp. From left to right, 100 bp DNA ladder; Lanes 1–10: samples; internal control was detected in all tested samples.
Supplementary Figure 2. (A) CBSD symptoms on cassava leaves; (B) CBSD symptoms on cassava stem; (C) CBSD symptoms on cassava roots; (D) CMD symptoms on cassava leaves.
Supplementary Figure 3. Histogram and P-P plot normality tests for cassava productivity linear regression analysis.
Supplementary Table 1. List of collected variables during the survey.
Supplementary Table 2. List of biotic variables recorded or measured during the survey of the 130 visited fields.
Supplementary Table 3. Socio-demographic factors stratified by cassava infection.
Supplementary Table 4. Disease awareness variables stratified by cassava infection.
Supplementary Table 5. Effects of demographic variables on cassava productivity.
Supplementary Datasheet 1. Questionnaire.
Footnotes
1. ^Available online at: https://www.newtimes.co.rw/business/what-holding-back-cassava-yield.
References
Abarshi, M. M., Mohammed, I. U., Wasswa, P., Hillocks, R. J., Holt, J., Legg, J. P., et al. (2010). Optimization of diagnostic RT-PCR protocols and sampling procedures for the reliable and cost-effective detection of Cassava brown streak virus. J. Virol. Methods 163, 353–359. doi: 10.1016/j.jviromet.2009.10.023
Alicai, T., Ndunguru, J., Sseruwagi, P., Tairo, F., Okao-okuja, G., Nanvubya, R., et al. (2016). A rapidly evolving genome: implications for virus speciation, variability, diagnosis and host resistance. Sci. Rep. 6, 1–14. doi: 10.1038/srep36164
Andrade-Piedra, J., Bentley, J., Almekinders, C., Walsh, K., and Thiele, G. (2016). Case Studies of Roots, Tubers and Banana Seed Systems. RTB Working Paper N° 2016-3. Lima: CGIAR Research Program on Roots, Tubers and Banana (RTB), 244. doi: 10.4160/23096586RTBWP20163
Anjanappa, R. B., Mehta, D., Okoniewski, M. J., Szabelska-Ber?sewicz, A., Gruissem, W., and Vanderschuren, H. (2018). Molecular insights into Cassava brown streak virus susceptibility and resistance by profiling of the early host response. Mol. Plant Pathol. 19, 476–489. doi: 10.1111/mpp.12565
Bentley, J., Olanrewaju, A., Madu, T., Olaosebikan, O., Abdoulaye, T., Assfaw Wossen, T., et al. (2017). Cassava Farmers' Preferences for Varieties and Seed Dissemination System in Nigeria: Gender and Regional Perspectives. IITA Monograph, Ibadan: IITA
Bigirimana, S., Barumbanze, P., Ndayihanzamaso, P., Shirima, R., and Legg, J. P. (2011). First report of cassava brown streak disease and associated Ugandan cassava brown streak virus in Burundi. New Dis. Rep. 24:26. doi: 10.5197/j.2044-0588.2011.024.026
Broek, J., and van den Byakweli, J. M. (2014). Exploratory Study on Rwanda's Seed Sector: Options for Dutch Support. CDI Report CDI-14-002. Wageningen: Centre for Development Innovtion, Wageningen UR (University and Research Centre).
Bua, B. (2017). Farmers' knowledge of cassava streak virus disease in selected districts of Central Uganda. Int. J. Environ. Agric. Biotechnol. 2, 230–234. doi: 10.22161/ijeab/2.1.29
Bull, S. E., Ndunguru, J., Gruissem, W., Beeching, J. R., and Vanderschuren, H. (2011). Cassava: constraints to production and the transfer of biotechnology to African laboratories. Plant Cell Rep. 30, 779–787. doi: 10.1007/s00299-010-0986-6
Burns, A., Gleadow, R., Cliff, J., Zacarias, A., and Cavagnaro, T. (2010). Cassava: the drought, war and famine crop in a changing world. Sustainability 2, 3572–3607. doi: 10.3390/su2113572
Campo, B. V. H., Hyman, G., and Bellotti, A. (2011). Threats to cassava production: known and potential geographic distribution of four key biotic constraints. Food Security 3, 329–345. doi: 10.1007/s12571-011-0141-4
Catholic Relief Services (2012). Final Report on the Great Lakes Cassava Initiative (Vol. 1). Catholic Relief Services. doi: 10.20622/jlta.1.0_toc1
Chavez, A. L., Sanchez, T., Jaramillo, G., Bedoya, J. M., and Echeverry, J. (2005). Variation of quality traits in cassava roots evaluated in landraces and improved clones. Euphytica 143, 125–133. doi: 10.1007/s10681-005-3057-2
Chipeta, M. M., Shanahan, P., Melis, R., Sibiya, J., and Ibrahim, R. M. (2016). Farmers' knowledge of cassava brown streak disease and its management in Malawi. Int. J. Pest Manage. 62, 175–184. doi: 10.1080/09670874.2016.1167268
Elegba, W. (2018). Engineering Cassava Mosaic Disease (CMD) Resistanc in a Ghanaian Cassava Cultivar (Doctoral dissertation). ETH Zurich.
El-Sharkawy, M. A. (2004). Cassava biology and physiology. Plant Mol. Biol. 53, 621–641. doi: 10.1007/s11103-005-2270-7
FAO (2016). The state of food and agriculture: climate change, agriculture and food security. Rome: Food and Agriculture Organization of the United Nations.
FAO (2018). “Food outlook—biannual report on global food markets—November 2018,” in Global Information and Early Warning System on Food and Agriculture. Geneva: FAO. doi: 10.1044/leader.PPL.19102014.18
Gasangwa, E. (2013). Socio-Economic Study of Cassava Value Chain in Rwanda. Sunnyvale, CA: Lambert Academic Publishing.
Gender Monitoring Office (2017). Republic of Rwanda, Gender and Agriculture. Gender Monitoring Office
Gondwe, F. M. T., Mahungu, N. M., Hillocks, R. J., Raya, M. D., Moyo, C. C., Soko, M. M., et al. (2003). Economic Losses Experienced by Small-scale Farmers in Malawi due to Cassava Brown Streak Virus Disease. London: Natural Resources International Limited.
Khonje, M., Mkandawire, P., Manda, J., and Alene, A. D. (2015). “Analysis of adoption and impacts of improved cassava varieties in Zambia,” Paper presented at the 29th Triennial Conference of the International Association of Agricultural Economists (IAAE) (Milan).
Koima, I. N., Orek, C. O., and Nguluu, S. N. (2018). Distribution of cassava mosaic and cassava brown streak diseases in agro-ecological zones of lower Eastern Kenya. Int. J. Innov. Sci. Res. Technol. 3, 391–399.
Kombate, K., Dossou-Aminon, I., Dansi, A., Adjatin, R. A., Kpemoua, K., Dassou, G. A., et al. (2017). Cassava (Manihot esculenta Crantz) production constraints, farmers' preference criteria and diversity management in Togo. Int. J. Curr. Microbiol. Appl. Sci. 6, 3328–3240. doi: 10.20546/ijcmas.2017.606.391
Legg, J. P., Okao-Okuja, G., Mayala, R., and Muhinyuza, J.-B. (2001). Spread into Rwanda of the severe cassava mosaic virus disease pandemic and the associated Uganda variant of East African cassava mosaic virus (EACMV-Ug). Plant Pathol. 50:796. doi: 10.1046/j.1365-3059.2001.00619.x
Legg, J. P., Owor, B., Sseruwagi, P., and Ndunguru, J. (2006). Cassava mosaic virus disease in East and Central Africa: epidemiology and management of a regional pandemic. Adv. Virus Res. 67, 355–418. doi: 10.1016/S0065-3527(06)67010-3
Legg, J. P., and Thresh, J. M. (2003). “Cassava virus diseases in Africa,” in Proceedings of a Conference on Plant Virology in Sub-Saharan Africa (Ibadan), 517–522.
Lobell, B. D., Burke, B. M., Tebaldi, C., Mastrandrea, D. M., Falcon, P. W., and Naylor, L. R. (2008). Needs for food security in 2030 region. Science 319, 607–610. doi: 10.1126/science.1152339
Maruthi, M. N., Legg, J. P., Jeremiah, S. C., and Mohammed, I. U. (2017). The role of the whitefly, Bemisia tabaci (Gennadius), and farmer practices in the spread of cassava brown streak ipomoviruses. J. Phytopathol. 165, 707–717. doi: 10.1111/jph.12609
Masumba, E. A., Kapinga, F., Mkamilo, G., Salum, K., Kulembeka, H., Rounsley, S., et al. (2017). QTL associated with resistance to cassava brown streak and cassava mosaic diseases in a bi-parental cross of two Tanzanian farmer varieties, Namikonga and Albert. Theor. Appl. Genet. 130, 2069–2090. doi: 10.1007/s00122-017-2943-z
Mbanzibwa, D. R., Tian, Y. P., Tugume, A. K., Mukasa, S. B., Tairo, F., Kyamanywa, S., et al. (2009). Genetically distinct strains of Cassava brown streak virus in the Lake Victoria basin and the Indian Ocean coastal area of East Africa. Arch. Virol. 1996, 353–359. doi: 10.1007/s00705-008-0301-9
Mbewe, W., Kumar, P. L., Changadeya, W., Ntawuruhunga, P., and Legg, J. (2014). Diversity, distribution and effects on cassava cultivars of cassava brown streak viruses in Malawi. J. Phytopathol. 163, 433–443. doi: 10.1111/jph.12339
Miklyaev, M., Jenkins, G., and Shobowale, D. (2021). Sustainability of agricultural crop policies in rwanda: an integrated cost–benefit analysis. Sustainability (Switzerland) 13, 1–22. doi: 10.3390/su13010048
Minot, N., Smale, M., Kling, J., and Eicher, C. (2007). Seed development programs in Sub-Saharan Africa: a review of experiences seed development programs in sub-saharan Africa. Nairobi: Rockefeller Foundation.
Mohammed, I. U., Abarshi, M. M., Muli, B., Hillocks, R. J., and Maruthi, M. N. (2012). The symptom and genetic diversity of cassava brown streak viruses infecting cassava in East Africa. Adv. Virol. 2012, 1–10. doi: 10.1155/2012/795697
Moreno, I., Gruissem, W., and Vanderschuren, H. (2011). Reference genes for reliable potyvirus quantitation in cassava and analysis of cassava brown streak virus load in host varieties. J. Virol. Methods 177, 49–54. doi: 10.1016/j.jviromet.2011.06.013
Mukasa, S. (2015). Developing a Community-Based Cassava Seed System for Increased Productivity and Market Linkages in Uganda. Cassava CARP Annual Report 2015.
Mulenga, R. M., Makulu, M., and Laura, M. (2018). Cassava brown streak disease and Ugandan cassava brown streak virus reported for the first time in Zambia. Plant Dis. 102, 1410–1418. doi: 10.1094/PDIS-11-17-1707-RE
Munganyinka, E., Ateka, E. M., Kihurani, A. W., Kanyange, M. C., and Tairo, F. (2018). Cassava brown streak disease in Rwanda, the associated viruses and disease phenotypes. Plant Pathol. 67, 377–387. doi: 10.1111/ppa.12789
Munyahali, W., Pypers, P., Swennen, R., Walangululu, J., Vanlauwe, B., and Merckx, R. (2017). Responses of cassava growth and yield to leaf harvesting frequency and NPK fertilizer in South Kivu, Democratic Republic of Congo. Field Crops Res. 214, 194–201. doi: 10.1016/j.fcr.2017.09.018
Muthoni, J., and Nyamongo, D. O. (2008). Seed systems in Kenya and their relationship to on-farm conservation of food crops. J. New Seeds 9, 330–342. doi: 10.1080/15228860802492273
Nduwumuremyi, A., Melis, R., Shanahan, P., and Asiimwe, T. (2016). Participatory appraisal of preferred traits, production constraints and postharvest challenges for cassava farmers in Rwanda. Food Security 8, 375–388. doi: 10.1007/s12571-016-0556-z
Night, G., Asiimwe, P., Gashaka, G., Nkezabahizi, D., Legg, P. J., Okao-Okuja, G., et al. (2011). Occurrence and distribution of cassava pests and diseases in Rwanda. Agric. Ecosyst. Environ. 140, 492–497. doi: 10.1016/j.agee.2011.01.014
Nyirahorana, C., Mburu, D. M., Mulyungi, P., Ntaganira, E., Nduguru, J., Sseruwagi, P., et al. (2017). Drivers behind adoption of cassava brown streak disease control measures in Rwanda. Int. J. Sci. Technol. Res. 6, 113–117.
Otekunrin, O. A., and Sawicka, B. (2019). Cassava, a 21st century staple crop: how can nigeria harness its enormous trade potentials? Acta Sci. Agric. 3, 194–202. doi: 10.31080/ASAG.2019.03.0586
Patil, B. L., Legg, J. P., Kanju, E., and Fauquet, C. M. (2015). Cassava brown streak disease: a threat to food security in Africa. J. General Virol. 96, 956–968. doi: 10.1099/jgv.0.000014
Rey, C., and Vanderschuren, H. (2017). Cassava mosaic and brown streak diseases: current perspectives and beyond. Annu. Rev. Virol. 4, 429–452. doi: 10.1146/annurev-virology-101416-041913
Rwegasira, G. M., Momanyi, G., Rey, M. E. C., Kahwa, G., and Legg, J. P. (2011). Widespread occurrence and diversity of cassava brown streak virus (Potyviridae: Ipomovirus) in Tanzania. Virology 101, 1159–1167.
Rwegasira, G. M., and Rey, H. N. (2011). Approaches to diagnosis and detection of cassava brown streak virus(Potiviridae, Ipomovirus) in field grown cassava crop. Afric. J. Food Agric. Nutr. Dev. 11, 4739–4756. doi: 10.1094/PHYTO-11-10-0297
Saranraj, P., Behera, S. S., and Ray, R. C. (2019). Traditional Foods From Tropical Root and Tuber Crops. Sawston: Woodhead Publishing. doi: 10.1016/B978-0-12-814887-7.00007-1
Sheat, S., Fuerholzner, B., Stein, B., and Winter, S. (2019). Resistance against cassava brown streak viruses from africa in cassava germplasm from South America. Front. Plant Sci. 10:567. doi: 10.3389/fpls.2019.00567
Shirima, R. R., Legg, J. P., Maeda, D. G., Tumwegamire, S., Mkamilo, G., Mtunda, K., et al. (2020). Genotype by environment cultivar evaluation for cassava brown streak disease resistance in Tanzania. Virus Res. 286:198017. doi: 10.1016/j.virusres.2020.198017
Spencer, D. S. C., and Ezedinma, C. (2017). Cassava Cultivation in Sub-Saharan Africa. Cambridge: Burleigh Dodds Science Publishing Limited, 123–148. doi: 10.19103/as.2016.0014.06
Sseruwagi, P., Sserubombwe, W. S., Legg, J. P., Ndunguru, J., and Thresh, J. M. (2004). Methods of surveying the incidence and severity of cassava mosaic disease and whitefly vector populations on cassava in Africa: a review. Virus Res. 100, 129–142. doi: 10.1016/j.virusres.2003.12.021
Keywords: cassava, seed system, CBSD, field survey, cassava production, Rwanda
Citation: Nyirakanani C, Bizimana JP, Kwibuka Y, Nduwumuremyi A, Bigirimana VdP, Bucagu C, Lassois L, Malice E, Gengler N, Massart S, Bragard C, Habtu M, Brostaux Y, Thonar C and Vanderschuren H (2021) Farmer and Field Survey in Cassava-Growing Districts of Rwanda Reveals Key Factors Associated With Cassava Brown Streak Disease Incidence and Cassava Productivity. Front. Sustain. Food Syst. 5:699655. doi: 10.3389/fsufs.2021.699655
Received: 23 April 2021; Accepted: 24 September 2021;
Published: 03 December 2021.
Edited by:
Katherine Steele, Bangor University, United KingdomReviewed by:
Makeshkumar Thangaraj, Central Tuber Crops Research Institute (ICAR), IndiaHale Tufan, Cornell University, United States
Copyright © 2021 Nyirakanani, Bizimana, Kwibuka, Nduwumuremyi, Bigirimana, Bucagu, Lassois, Malice, Gengler, Massart, Bragard, Habtu, Brostaux, Thonar and Vanderschuren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chantal Nyirakanani, Q2hhbnRhbC5OeWlyYWthbmFuaSYjeDAwMDQwO3VsaWVnZS5iZQ==; Jean Pierre Bizimana, SmVhbi5CaXppbWFuYSYjeDAwMDQwO3VsaWVnZS5iZQ==; Hervé Vanderschuren, aGVydmUudmFuZGVyc2NodXJlbiYjeDAwMDQwO2t1bGV1dmVuLmJl; aGVydmUudmFuZGVyc2NodXJlbiYjeDAwMDQwO3VsaWVnZS5iZQ==
†These authors have contributed equally to this work