Epidemiological assessment of cassava mosaic disease in Burkina Faso

Received: 6 April 2021 | Accepted: 5 August 2021Β 

DOI: 10.1111/ppa.13459Β Β 

ORIGINAL ARTICLEΒ 

Epidemiological assessment of cassava mosaic disease inΒ  Burkina FasoΒ 

Monique Soro1,2,3,4 | Fidèle Tiendrébéogo1,3,4 | Justin S. Pita1,2 |  Edwig T. Traoré3,4,5 | Koussao Somé3,4,6 | Ezechiel B. Tibiri3,4 | James B. Néya3,4 |  J. Musembi Mutuku1 | Jacques Simporé5 | Daouda Koné2,7 

1Central and West African VirusΒ  Epidemiology (WAVE), PΓ΄le scientifiqueΒ  et d’innovation de Bingerville, Université  FΓ©lix HouphouΓ«t-Boigny (UFHB),Β  Bingerville, Ivory CoastΒ 

2Laboratoire de Biotechnologie, AgricultureΒ  et Valorisation des Ressources Biologiques,Β  UFR Biosciences, UniversitΓ© FΓ©lixΒ  HouphouΓ«t-Boigny, Abidjan, Ivory Coast 3Laboratoire de Virologie et deΒ Β 

Biotechnologies VΓ©gΓ©tales, Institut deΒ  l’Environnement et de Recherches AgricolesΒ  (INERA), Ouagadougou, Burkina Faso 4Laboratoire Mixte International Patho-Bios,Β  IRD-INERA, Ouagadougou, Burkina Faso 5Laboratoire de Biologie MolΓ©culaire et deΒ  GΓ©nΓ©tique (LABIOGENE), UniversitΓ© JosephΒ  Ki-Zerbo, Ouagadougou, Burkina Faso 6Laboratoire de GΓ©nΓ©tique et deΒ Β 

Biotechnologies VΓ©gΓ©tales, Institut deΒ  l’Environnement et de Recherches AgricolesΒ  (INERA), Ouagadougou, Burkina Faso 7Centre d’Excellence Africain sur leΒ  Changement Climatique, la Biodiversité  et l’Agriculture Durable (WASCAL/CEA CCBAD, UniversitΓ© FΓ©lix HouphouΓ«t Boigny), PSI-UniversitΓ© FΓ©lix HouphouΓ«t Boigny, Abidjan, Ivory CoastΒ 

CorrespondenceΒ 

FidΓ¨le TiendrΓ©bΓ©ogo, Laboratoire deΒ  Virologie et de Biotechnologies VΓ©gΓ©tales,Β  Institut de l’Environnement et deΒ  Recherches Agricoles (INERA), 01 BP 476Β  Ouagadougou 01, Burkina Faso.Β 

Emails: fidelet@gmail.com; fidele. tiendrebeogo@wave-center.orgΒ 

Funding informationΒ 

Bill & Melinda Gates Foundation; TheΒ  United Kingdom Foreign, CommonwealthΒ  & Development Office, Grant/AwardΒ  Number: OPP1082413Β 

AbstractΒ 

Surveys were conducted in 2016 and 2017 across the main cassava-growing regionsΒ  of Burkina Faso to assess the status of cassava mosaic disease (CMD) and to determineΒ  the virus strains causing the disease, using field observation and phylogenetic analy sis. CMD incidence varied between regions and across years but was lowest in Hauts Bassins (6.0%, 2016 and 5.4%, 2017) and highest in Centre-Sud (18.5%, 2016) and inΒ  Boucle du Mouhoun (51.7%, 2017). The lowest CMD severity was found in Est regionΒ  (2.0) for both years and the highest in Sud-Ouest region (3.3, 2016) and Centre-SudΒ  region (2.8, 2017). The CMD infection was primarily associated with contaminatedΒ  cuttings in all regions except in Hauts-Bassins, where whitefly-borne infection wasΒ  higher than cuttings-borne infection in 2016. PCR screening of 687 samples coupledΒ  with sequence analysis revealed the presence of African cassava mosaic-like (ACMV like) viruses and East African cassava mosaic-like (EACMV-like) viruses as single infec tions at 79.5% and 1.1%, respectively. Co-infections of ACMV-like and EACMV-likeΒ  viruses were detected in 19.4% of the tested samples. In addition, 86.7% of the sam ples positive for EACMV-like virus were found to be positive for East African cassavaΒ  mosaic Cameroon virus (EACMCMV). Phylogenetic analysis revealed the segregationΒ  of cassava mosaic geminiviruses (CMGs) from Burkina Faso into three clades specificΒ  to ACMV, African cassava mosaic Burkina Faso virus (ACMBFV), and EACMCMV, con firming the presence of these viruses. The results of this study show that EACMCMVΒ  occurrence may be more prevalent in Burkina Faso than previously thought.Β 

KEYWORDSΒ 

African cassava mosaic virus (ACMV), cassava mosaic geminiviruses, East African cassavaΒ  mosaic Cameroon virus (EACMCMV), geminiviruses characterization, geminivirusesΒ  distribution

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium,Β  provided the original work is properly cited.Β 

Β© 2021 The Authors. Plant Pathology published by John Wiley & Sons Ltd on behalf of British Society for Plant Pathology.Β 

Plant Pathology. 2021;70:2207–2216. wileyonlinelibrary.com/journal/ppa | 2207Β 

2208 |    SORO et al.Β 

1  | INTRODUCTIONΒ 

Cassava (Manihot esculenta, Euphorbiaceae), which originates fromΒ  Latin America, is a major source of food for more than 700 millionΒ  people in tropical and subtropical developing countries and en hances food security in these countries (Ntawuruhunga et al., 2013;Β  Patil & Fauquet, 2009; Saediman et al., 2016). Cassava is a stapleΒ  food crop in the sub-Saharan region of Africa and consequently aΒ  source of income for many processors and traders (NtawuruhungaΒ  et al., 2013). The high calorie yield per hectare (250 kcal/ha/day),Β  drought tolerance, hardiness in stressful environments, and flexibil ity of harvesting time are the major advantages of this crop com pared to many other crops (Byju & Suja, 2020; El-Sharkawy, 2004;Β  Pushpalatha & Gangadharan, 2020). In Burkina Faso, cassava wasΒ  introduced by farmers decades ago from Ghana and Ivory CoastΒ  (CΓ΄te d’Ivoire) (Guira et al., 2017). It has long been cultivated aroundΒ  vegetable gardens for domestic consumption. Formerly consideredΒ  as a neglected crop, cassava has become a cash crop since the for mal introduction of improved varieties from IITA in 2003 (DabirΓ© &Β  Belem, 2003).Β 

In Africa, cassava production is negatively affected by two mainΒ  viral diseases: cassava brown streak disease (CBSD) and cassava mo saic disease (CMD). CMD is a major constraint to cassava production,Β  which causes tuber yield losses estimated at $2.7 billion annuallyΒ  (Patil & Fauquet, 2009). CMD is caused by distinct cassava mosaicΒ  geminiviruses (CMGs) (family Geminiviridae, genus Begomovirus) andΒ  naturally transmitted by the whitefly Bemisia tabaci (Hemiptera:Β  Aleyrodidae) (Ally et al., 2019; Legg et al., 2002; MacFadyen et al.,Β  2018). CMD is also widely disseminated by infected stem cuttings,Β  used for vegetative propagation (Bock & Woods, 1983; FondongΒ  et al., 2000). CMD is endemic in Africa, with nine distinct CMGΒ  species officially recognized by the International Committee onΒ  Taxonomy of Viruses (ICTV; https://talk.ictvonline.org/ictv-reporΒ 

ts/ictv_online_report/ssdna-viruses/w/geminiviridae/479/membe r-species-begomovirus): African cassava mosaic Burkina Faso virusΒ  (ACMBFV; TiendrΓ©bΓ©ogo et al., 2012), African cassava mosaic virusΒ  (ACMV; Stanley & Gay, 1983), East African cassava mosaic CameroonΒ  virus (EACMCMV; Fondong et al., 2000), East African cassava mo saic Kenya virus (EACMKV; Bull et al., 2006), East African cassavaΒ  mosaic Malawi virus (EACMMV; Zhou et al., 1998), East AfricanΒ  cassava mosaic virus (EACMV; Pita, Fondong, SangarΓ©, Otim-Nape,Β  et al., 2001), East African cassava mosaic Zanzibar virus (EACMZV;Β  Maruthi et al., 2004), cassava mosaic Madagascar virus (CMMGV),Β  and South African cassava mosaic virus (SACMV; Berrie et al., 2001).Β 

In West Africa, the presence of ACMV and EACMV was reportedΒ  in Ivory Coast (Pita, Fondong, SangarΓ©, Kokora, et al., 2001; ToualyΒ  et al., 2014), Ghana (Torkpo et al., 2017), and Nigeria (AbubakarΒ  et al., 2019; Ariyo et al., 2005; Eni et al., 2021; Ogbe et al., 2006). TheΒ  presence of EACMCMV was also reported in Ivory Coast and NigeriaΒ  (Ariyo et al., 2005; Eni et al., 2021; Pita, Fondong, SangarΓ©, Kokora,Β  et al., 2001). In previous studies, the presence of cassava mosaic disΒ 

ease (CMD) has been reported in some localities in Burkina Faso.Β  Indeed, the presence of ACMV was reported in 1995 using the tripleΒ Β 

antibody sandwich-ELISA method with cross-reacting monoclonalΒ  antibodies to ACMV (KonatΓ© et al., 1995). The molecular featuresΒ  of an ACMV-like virus (ACMBFV, whose Rep protein gene and in tergenic region differ from ACMV) was described and the pres ence of EACMV-UG variant was reported around OuagadougouΒ  (TiendrΓ©bΓ©ogo et al., 2009, 2012). Since then, the real status of CMDΒ  and its epidemiological parameters such as the incidence and sever ity of the disease, the whitefly abundance, and the mode of infectionΒ  remain unclear. To overcome this knowledge gap, we conducted forΒ  the first time georeferenced surveys in the main cassava productionΒ  areas in Burkina Faso.Β 

2  |  MATERIALS AND METHODSΒ 

2.1  |  Cassava mosaic disease status assessmentΒ 

Surveys were conducted in 2016 and 2017 in eight major cassava growing regions of Burkina Faso (Figure 1). The number of fieldsΒ  sampled within a region depended on the number of cassava growing localities and the availability of cassava fields at 3–6 monthsΒ  after planting (MAP). Our field sampling protocol was a modificationΒ  of one previously described (Sseruwagi et al., 2004) and has beenΒ  adopted by 10 countries in Central and West Africa to harmonize ef forts at surveillance and monitoring of these transboundary patho gens of high economic importance. Briefly, in each field, 30 cassavaΒ  plants were assessed randomly along two diagonals to form an β€œX”  pattern. Then each selected plant was assessed visually for the pres ence or absence of CMD symptoms (leaf mosaic, leaf distortion, andΒ  stunted growth) and the number of whiteflies settling on the leaves,Β  and if infected, we determined whether the mode of infection wasΒ  through cuttings or whitefly transmission. The whitefly popula tion was estimated by counting the number of whiteflies on theΒ  top five fully expanded leaves. The mode of infection in each plantΒ  was determined based on the location of the leaves with symptomsΒ  as previously described by Sseruwagi et al. (2004). According toΒ  these authors, from 3 to 6 MAP it is possible to distinguish betweenΒ  cutting-borne and whitefly-borne infections. Symptoms appearingΒ  only on upper leaves were taken to have resulted from whitefly transmitted infection, whereas plants that showed symptoms eitherΒ  only on the lower leaves or on all leaves were taken as having beenΒ  infected through cassava cuttings. CMD symptom severity was re corded using a scale from 1 (no symptoms) to 5 (very severe symp toms) (Terry, 1975). We acknowledge that the severity level dependsΒ  on the variety, climate, crop management, and mainly the time atΒ  which the infection occurred. To minimize the effects of these vari ables on our data, we sampled fields within the same locations thatΒ  were within the 3–6 MAP age. The CMD incidence was calculatedΒ  as the percentage of plants with symptoms in relation to the numberΒ  of plants assessed.Β 

A total of 237 leaf samples from 212 plants with symptoms andΒ  25 without symptoms in 2016, and 450 leaf samples from 240 plantsΒ  with symptoms and 210 without symptoms in 2017 were collectedΒ 

   SORO et al. |  2209

FIGURE 1 Map of Burkina Faso showing the regions and localities where surveys were done in 2016 and 2017 [Colour figure can beΒ  viewed at wileyonlinelibrary.com]Β 

TABLE 1 Primer pairs used for the amplification of ACMV-like virus, EACMV-like virus, and EACMCMV ExpectedΒ Β 

Primer Sequence (5′–3β€²) Target regionΒ 

size (bp) Virus species ReferenceΒ 

JSP 001 ATGTCGAAGCGACCAGGAGAT DNA-A (CP) 783 ACMV-like virus Pita, Fondong, SangarΓ©, Otim-Nape, et al.Β  (2001) JSP 002 TGTTTATTAATTGCCAATACTΒ 

JSP 001 ATGTCGAAGCGACCAGGAGAT DNA-A (CP) 780 EACMV-like virus Pita, Fondong, SangarΓ©, Otim-Nape, et al.Β  (2001) JSP 003 CCTTTATTAATTTGTCACTGCΒ 

VNF031/F GGATACAGATAGGGTTCCCAC DNA-A (AC2/Β 

560 EACMCMV Fondong et al. (2000)Β 

VNF032/R GACGAGGACAAGAATTCCAAT

AC3)Β 

for laboratory analysis using PCR and Sanger sequencing. The globalΒ  positioning system (GPS) coordinates were recorded for each field.Β 

2.2  |  Molecular characterization of CMGsΒ 

Total DNA was extracted from cassava leaves using the CTAB pro tocol as previously described (Permingeat et al., 1998). The concen tration of DNA in each sample was determined using a NanoDropΒ  2000 spectrophotometer (Thermo Fisher Scientific) and adjustedΒ  to 150 ng/ΞΌl. We previously discovered that the most problematicΒ  CMGs in smallholder cassava production in Burkina Faso were ACMVΒ  and a variant of the East African cassava mosaic virus (EACMV), theΒ Β 

EACMV-Uganda variant (TiendrΓ©bΓ©ogo et al., 2009). Because theΒ  current status of the incidence and severity of these two CMD causing viruses is unknown in Burkina Faso, we surveyed the wholeΒ  country and focused our surveys on ACMV and EACMV. The ex tracted DNA was subjected to PCR using the specific primers listedΒ  in Table 1 for the detection of ACMV-like virus (JSP001/JSP002)Β  and EACMV-like virus (JSP001/JSP003). The samples positive forΒ  EACMV-like virus were subjected to another round of PCR usingΒ  specific primers for the detection of EACMCMV (VNF031/VNF032;Β  Table 1). The PCR mix was prepared in a final volume of 25 ΞΌl usingΒ  20.9 ΞΌl of molecular biology grade water, 2.5 Β΅l of 10Γ— reactionΒ  buffer, 0.5 Β΅l of 10 mM dNTPs, 0.5 Β΅l of 10 Β΅M of each primer, 0.1 Β΅lΒ  of 5 U/Β΅l of Maximo Taq DNA polymerase (GeneON), and 150 ngΒ Β 

2210 |    SORO et al.Β 

DNA template of each sample. The DNA amplification was carriedΒ  out in a SimpliAmp thermal cycler (Life Technologies Holdings PteΒ  Ltd). The PCR temperature profile was set at 94Β°C for 4 or 5 min forΒ  initial denaturation, followed by 35 cycles of amplification at 94Β°CΒ  for 45 or 60 s, 55Β°C for 45 or 60 s, and 72Β°C for 55 or 60 s (dependingΒ  on primers). The final elongation step was performed at 72Β°C for 7Β  or 10 min. PCR-amplified products were subjected to 1% agarose gelΒ  electrophoresis, stained with ethidium bromide. The electrophoresisΒ  was performed at 100 V and the gel was visualized using a CompactΒ  Digimage System, UVDI series (MS major science). PCR products ofΒ  40 ACMV-like positive samples (randomly selected from the regions)Β  were directly sequenced in both forward and reverse orientationsΒ  using the Sanger method at Inqaba Biotec company (South Africa)Β  to determine their identity. PCR products of 15 EACMCMV posi tive samples were also subjected to sequencing in both forward andΒ  reverse orientations to confirm their identity.Β 

2.3  |  Statistical analysisΒ 

Data analysis was performed using the R software v. 3.6.1 (RΒ  Development Core Team). The normality of the variables was deter mined using the Shapiro–Wilk test. When the variable was not dis tributed according to the normal distribution, the generalized linearΒ  model was used. The difference in the number of whiteflies per plantΒ  between regions and the difference in the severity score of CMDΒ  between regions in the same year were assessed using the general ized linear model and Tukey’s pairwise mean comparison test. TheΒ  difference in the number of whiteflies per plant and the differenceΒ  in the severity score of CMD between 2016 and 2017 were assessedΒ  using Wilcoxon test with continuity correction. A test of pairwiseΒ  comparison of proportions was used based on a G-test with correc tion of BY (Benjamini & Yekutieli, 2001) to compare the incidencesΒ  of CMD between regions. The map of Burkina Faso showing the re gions where surveys were done in 2016 and 2017 was developedΒ  using QGIS software v. 2.18.26 (https://qgis.org/downloads/).Β 

2.4  |  Phylogenetic analysisΒ 

The amplicon sequences were trimmed and assembled de novo usingΒ  Geneious v. 8.1.7 (Biomatters Ltd) software. Consensus sequenceΒ  obtained from forward and reverse sequences for each sample wasΒ  subjected to BLASTn in NCBI for preliminary species assignment andΒ  subsequently for pairwise sequence comparison (Bao et al., 2014).Β  The sequences were aligned with representative isolates of begoΒ 

moviruses using ClustalW alignment method in MEGA X softwareΒ  (Kumar et al., 2018). The sequences of 25 out of 40 ACMV-like virusΒ  positive samples and six out of 15 EACMCMV positive samples wereΒ  used for phylogenetic tree construction. The maximum-likelihoodΒ  (ML) method with general time reversible (GTR) model (as the bestΒ  fit model for substitution pattern description) was used for phyloΒ 

genetic trees construction using FastTree v. 2.1.9 (Price et al., 2010)Β  with bootstrap replicates of 1000. The tree was visualized and ed ited using FigTree v. 1.4.4 (http://tree.bio.ed.ac.uk/software/figtr ee/).Β 

3  |  RESULTSΒ 

3.1  |  CMD distribution in 2016 and 2017Β 

In 2016, CMD symptoms were found in 84.0% (42/50) of surveyedΒ  localities, with the lowest proportion (57.1%, 8/14) of infectedΒ  fields in the region of Hauts-Bassins. Cassava fields in 65.9%Β  (29/44) of localities showed CMD symptoms in 2017, with the lowΒ 

est proportion (45.5%, 5/11) in Centre-Est region (Table 2). TheΒ  proportion of localities where CMD-affected cassava fields wereΒ  found varied significantly between 2016 and 2017 (p < 0.05).Β  Indeed, compared to 2016, no cassava fields were found to haveΒ  CMD symptoms in the provinces of Nahouri (Centre-Sud region)Β  and Bougouriba (Sud-Ouest region) in 2017. Typical CMD sympΒ 

toms observed across farmers’ fields included distinctive leaf mo saic symptoms often associated with leaf distortion and reductionΒ Β 

TABLE 2 Proportion of localities in Burkina Faso where typical cassava mosaic disease (CMD) symptoms were found on cassava plants inΒ  2016 and 2017Β 

2016 2017Β 

RegionΒ 

SurveyedΒ  localitiesΒ 

Localities withΒ  CMDΒ 

Localities withΒ  CMD (%)Β 

SurveyedΒ  localitiesΒ 

Localities withΒ  CMDΒ 

Localities withΒ  CMD (%)Β 

Boucle du Mouhoun 5 4 80.0 2 2 100.0 Cascades 11 10 91.9 7 7 100.0 Centre-Est 1 1 100.0 11 5 45.4 Centre-Ouest 6 6 100.0 6 5 83.3 Centre-Sud 8 8 100.0 4 2 50.0 Est 2 2 100.0 2 2 100.0 Hauts-Bassins 14 8 57.1 8 4 50.0 Sud-Ouest 3 3 100.0 4 2 50.0 Total 50 42 84.0 44 29 65.9

   SORO et al. |  2211Β 

(Figure 2b–e), as well as an overall stunted appearance of the af fected plants.Β 

3.2  |  Incidence and symptom severity of CMD inΒ  2016 and 2017Β 

The CMD incidence in 2016 varied significantly from that observed inΒ  2017 (p ≀ 0.001). In 2016, the overall CMD incidence across the sur veyed fields in Burkina Faso was 11.3% (216/1920) and ranged fromΒ  6.0% (36/600) in Hauts-Bassins region to 18.5% (50/270) in Centre Sud region. In 2016, the difference between the lowest incidenceΒ  and the highest incidence was highly significant (p ≀ 0.001). For theΒ  2017 survey, the overall CMD incidence was 18.9% (329/1720). TheΒ  lowest incidence in 2017 was observed in the Hauts-Bassins regionΒ  with 5.4% (21/390), whereas the Boucle du Mouhoun region wasΒ  observed to have the highest incidence 51.7% (30/60, p ≀ 0.001;Β  Figure 3a).Β 

The mean CMD symptom severity score was 2.9 and the rangeΒ  was from 2.0 (Est region) to 3.3 (Sud-Ouest region) in 2016. In 2017,Β  the mean CMD symptom severity score was 2.5 with the lowest se verity in Est region (2.0) and the highest severity in the Centre-SudΒ  region (2.8) (Figure 3b). Significant differences (p < 0.05) were foundΒ  between CMD symptom severity scores in 2016 and 2017 in threeΒ  regions (Cascades, Centre-Ouest, and Hauts-Bassins). In most of theΒ  regions, no significant difference was found between the proportionΒ  of plants with different CMD symptom severity scores in 2016. InΒ  2017, the proportion of plants with CMD symptom severity score 2.0Β Β 

was higher than the other symptom severity scores in most of the re gions, whereas in the region of Centre-Sud the proportion of plantsΒ  with CMD symptom severity score 3 was the highest (Figure 3c).Β 

3.3  |  Adult whitefly distribution andΒ Β 

mode of infectionΒ 

Determination of whitefly counts and distribution was conductedΒ  at the time of the CMD incidence and severity survey to ensure theΒ  parameters that might affect the epidemiology of CMD in the fieldΒ  were the same, and the plants were of similar age. The methods haveΒ  been harmonized across 10 West and Central African countries,Β  Sierra Leone, Ivory Coast, Burkina Faso, Ghana, Nigeria, Benin, Togo,Β  Cameroon, Gabon, and Democratic Republic of Congo, to ensure theΒ  data is comparable. The adult whitefly counts were very low in 2016Β  and a similar trend was observed in 2017, with a mean of 0.1 and 0.7Β  per plant in 2016 and 2017, respectively. The highest mean whiteflyΒ  count was observed in the region of Boucle du Mouhoun (1.08) inΒ  2016. In 2017, the highest mean whitefly count (2.7) was observedΒ  in the regions of Boucle du Mouhoun and Cascades. In most of theΒ  regions, the mean number of whiteflies per plant was higher in 2017Β  than in 2016 (Figure 3d).Β 

When the number of plants with symptoms infected throughΒ  cuttings or by whitefly transmission were compared, a prepon derance of cutting-borne infections was detected in 2016 (83.3%,Β  180/216) and 2017 (88.8%, 292/329). The exception was Hauts Bassins region, where greater whitefly-borne infections wereΒ Β 

FIGURE 2 Symptoms of cassava mosaic disease observed on infected cassava plants during the surveys, using a scale from 1 (noΒ  symptoms) to 5 (very severe symptoms). (a) = 1, (b) = 2, (c) = 3, (d) = 4, (e) = 5 [Colour figure can be viewed at wileyonlinelibrary.com]

2212 |    SORO et al.

FIGURE 3 Epidemiological assessment of cassava mosaic disease (CMD) in Burkina Faso. (a) CMD incidence (percentage of plants withΒ  symptoms). (b) Severity of CMD (mean CMD severity score of plants with symptoms). (c) Proportion of plants with different CMD severityΒ  scores in 2016 and 2017. (d) Mean whitefly counts in 2016 and 2017. (e) Proportion of plants with symptoms infected by cutting or whiteflyΒ  in 2016 and 2017. The bars represent the standard error. Bars sharing the same lower case letters are not significantly different betweenΒ  regions in 2016 and those sharing the same upper case letters are not significantly different between regions in 2017 [Colour figure can beΒ  viewed at wileyonlinelibrary.com]

recorded in 2016 (Figure 3e). It is likely that a number of factors in cluding climatic conditions, infection status, or indeed cassava vari ety differences affected whitefly counts (Mugerwa et al., 2021). TheΒ  challenge is that during our survey years, whitefly pressure was notΒ  strong enough to explain the significant differences in disease inci dence observed between regions and years. Also, considering thatΒ  high incidence is associated with seeding using infected cuttings, weΒ  are not able to without doubt correlate whitefly numbers with theΒ  disease incidence or severity.Β 

3.4  |  CMGs detected by PCR in cassavaΒ  leaf samplesΒ 

A total of 687 cassava leaf samples were collected from 452 plantsΒ  with symptoms and 235 plants without symptoms in 2016 and 2017Β  for PCR analysis. Among the samples having observable symptoms,Β  4.0% (18/452) tested negative for ACMV-like virus and EACMV-likeΒ  virus. On the other hand, 2.1% (5/235) of symptomless samplesΒ  tested positive for ACMV-like virus. Approximately 63.9% (439/687)Β Β 

   SORO et al. |  2213 TABLE 3 PCR results obtained from samples collected during 2016 and 2017 surveys in eight main cassava-growing regions of BurkinaΒ  FasoΒ 

ACMV-like virus singleΒ Β 

infectionΒ 

PositiveΒ Β 

EACMV-like virusΒ Β 

single infection Mixed infectionΒ 

Region Tested samplesΒ 

samplesΒ 

n % n % n %Β 

Boucle du Mouhoun 30 26 19 73.1 ad 0 0.0 7 26.9 ab Cascades 202 132 107 81.1 ab 3 2.3 22 16.6 bc Centre-Est 90 43 43 100.0 c 0 0.0 0 0.0 d Centre-Ouest 95 70 62 88.6 ab 0 0.0 8 11.4 bc Centre-Sud 84 68 38 55.9 d 1 1.5 29 42.6 a Est 35 18 18 100.0 c 0 0.0 0 0.0 d Hauts-Bassins 113 54 48 88.9 ab 0 0.0 0 11.1 bc Sud-Ouest 38 28 14 50.0 d 1 3.6 13 46.4 a Total 687 439 349 79.5 5 1.1 85 19.4 Note: Percentages followed by the same letters are not significantly different between regions.

of collected samples tested positive for CMGs. Among the positiveΒ  samples, the single ACMV-like virus infection was by far the mostΒ  frequent, accounting for 79.5% (349/439) of all infection, followedΒ  by mixed infections of ACMV-like virus and EACMV-like virus withΒ  19.4% (85/439), and single infection of EACMV-like virus with 1.1%Β  (5/439). The single infection of ACMV-like virus was predominantΒ  in all surveyed regions, with the highest proportion (100%) in theΒ  Centre-Est and Est regions. The mixed infection occurred in theΒ  remaining six regions, with the highest proportions in Centre-SudΒ  (42.6%, 29/68) and Sud-Ouest (46.4%, 13/28) regions. The single in fection of EACMV-like virus was found in the regions of Sud-OuestΒ  (3.6%, 1/28), Cascades (2.3%, 3/132) and Centre-Sud (1.5%, 1/68)Β  but no significant difference was found between these proportionsΒ  (Table 3). Of the 90 EACMV-like virus positive samples (single andΒ  mixed infections), 86.7% (78/90) tested positive for EACMCMVΒ  using the primer pair VNF031/VNF032.Β 

3.5  |  CMGs identity confirmed by sequencingΒ 

A search for related sequences in the GenBank database (NCBI,Β  BLASTN) showed that the sequences of the 40 samples that testedΒ  positive for the ACMV-like virus were most closely related to ACMVΒ  and ACMBFV. Indeed, they shared the highest nucleotide identityΒ  (98%–99%) with ACMV isolates from Ghana (MG250119, MG250156,Β  MG250088), Ivory Coast (AF259894), Burkina Faso (FM877473),Β  and Nigeria (MH251339), and with ACMBFV isolates from BurkinaΒ  Faso (HE616777, HE616779, HE616780, HE616781). The sequencesΒ  of the 15 samples that tested positive for EACMCMV were mostΒ  closely related to the EACMCMV and shared the highest nucleotideΒ  identities (97%–98%) with isolates from Ghana (MG250164), IvoryΒ  Coast (AF259896), Nigeria (EU685319, EU685326), and MadagascarΒ  (KJ887944). The ML phylogenetic tree inferred from alignment ofΒ  coat protein (CP) gene sequences from Burkina Faso (25 ACMV like virus and six EACMCMV) and other CMGs confirmed that theΒ Β 

sequences from Burkina Faso are phylogenetically associated withΒ  ACMV, ACMBFV, or EACMCMV (Figure 4).Β 

4  | DISCUSSIONΒ 

In general, the mean CMD severity scores recorded in the studyΒ  areas were moderate in both 2016 and 2017. This could be due toΒ  similarities in factors affecting disease establishment across theΒ  country. The significant difference observed between the proporΒ 

tion of localities with CMD-affected fields in 2016 and 2017 couldΒ  be explained by better awareness by farmers of the risk of CMDΒ  transmission via infected cuttings through outreach programmesΒ  initiated in 2016.Β 

The relatively lower incidence recorded in Burkina Faso couldΒ  be explained by the fact that the intensification of cassava produc tion is a more recent phenomenon in the country compared to otherΒ  African countries (Guira et al., 2017; Legg et al., 2006), coupled withΒ  CMD awareness and the adoption of good farming practices byΒ  farmers.Β 

Our results showed that cases of CMD transmitted by cassavaΒ  cuttings were more prevalent as compared to cases resulting fromΒ  whitefly transmission. This phenomenon appears to be widespreadΒ  in sub-Saharan Africa (Chikoti et al., 2013; Mulenga et al., 2016;Β  Mwatuni et al., 2015; Torkpo et al., 2018; Zinga et al., 2013). TheΒ  high incidence of cutting-borne infection is probably due to farmΒ 

ers’ inability to select virus-free cassava cuttings when planting. TheΒ  very low incidences of whitefly-borne infections observed in BurkinaΒ  Faso is consistent with the low counts of whiteflies observed in theΒ  cassava fields in both years under study. It is notable that althoughΒ  the mean whitefly counts in the Hauts-Bassins region was less thanΒ  one per plant in 2016, a higher proportion of whitefly-borne infecΒ 

tions were recorded from the region during the same period. TheseΒ  results can be interpreted as suggesting that the rate of whitefly borne infection is not always correlated with whitefly abundance,Β Β 

2214 |    SORO et al.

FIGURE 4 Maximum-likelihood phylogenetic tree obtained from alignment of partial nucleotide sequences of coat protein (CP) genesΒ  of African cassava mosaic-like viruses (ACMV-like) and East African cassava mosaic-like viruses (EACMV-like). The names of the sequencesΒ  characterized in this study are in red. The horizontal scale indicates the genetic distance [Colour figure can be viewed at wileyonlinelibrary. com]

as was recently reported by Eni et al. (2021). Although our resultsΒ  showed that whiteflies may not be a key factor in the epidemiologyΒ  of CMD in our study area, it would be interesting to conduct otherΒ  field experiments using CMD-free planting material, in differentΒ  localities and at different times of the year, to determine the roleΒ  played by whiteflies in CMD epidemiology.Β 

This study shows the presence of ACMV-like viruses in the eightΒ  cassava-growing regions and EACMV-like viruses in six of them,Β  occurring as single or mixed infections in CMD-affected cassavaΒ Β 

plants in Burkina Faso. This is probably due to the exchange ofΒ  planting material between Burkina Faso and the neighbouringΒ  countries such as Ivory Coast, Togo, and Ghana where ACMV-likeΒ  viruses and EACMV-like viruses have been also reported (AdjataΒ  et al., 2009; Torkpo et al., 2017; Toualy et al., 2014). ACMV-like viΒ 

ruses were the predominant CMGs species in each cassava-growingΒ  region as the majority of CMD resulted from single ACMV-like virusΒ  infections. The predominance of single ACMV-like virus infection inΒ  West Africa has previously been reported (Abubakar et al., 2019;Β Β 

   SORO et al. |  2215Β 

Ogbe et al., 2006; Pita, Fondong, SangarΓ©, Kokora, et al., 2001;Β  Toualy et al., 2014) accompanied by a low distribution of EACMV like virus single infections (Ariyo et al., 2005; Ogbe et al., 2006;Β  Toualy et al., 2014). Our current work confirms that the situationΒ  has not changed. In addition, we discovered that most EACMV like virus isolates occurred as mixed infections with ACMV. OverΒ  86% of the EACMV-like virus positive samples were found to haveΒ  East African cassava mosaic Cameroon virus (EACMCMV). TheseΒ  results show that EACMCMV occurrence may be more prevalentΒ  in Burkina Faso than previously thought. Our analysis confirms thatΒ  the CMG isolates obtained from Burkina Faso samples are phyloge netically associated with ACMV-like viruses (ACMV and ACMBFV)Β  and EACMCMV. We propose further analysis, such as the use ofΒ  specific primers for each CMG species or next-generation sequenc ing, to resolve the issue of the occurrence of CMG species andΒ  strains in Burkina Faso.Β 

We detected the occurrence of CMGs in symptomless samplesΒ  (2.1%), which shows that the viruses can be latent in the plantsΒ  without manifesting symptoms. Therefore, the use of symptom less cassava landraces as an option to manage CMD could inad vertently result in increased cutting-borne transmission becauseΒ  they may harbour CMGs. We propose that the use of certifiedΒ  virus-free cuttings for the establishment of new cassava fields willΒ  be crucial for fighting the transmission of CMD. In the absenceΒ  of certified virus-free cuttings, the training of farmers on how toΒ  select healthy cuttings for the new planting season and on useΒ  of in-field diagnostic applications will be crucial to bring downΒ  the incidence or transmission of these viruses of high economicΒ  importance.Β 

ACKNOWLEDGEMENTSΒ 

This work was supported, in whole or in part, by the Bill & MelindaΒ  Gates Foundation and The United Kingdom Foreign, CommonwealthΒ  & Development Office (FCDO) under grant number OPP1082413 toΒ  the Central and West African Virus Epidemiology (WAVE) ProgramΒ  for root and tuber cropsβ€”through a subgrant from UniversitΓ© FΓ©lixΒ  HouphouΓ«t-Boigny (UFHB) to the Institut de l’Environnement et deΒ  Recherches Agricoles (INERA). Under the grant conditions of theΒ  Foundation, a Creative Commons Attribution 4.0 Generic LicenseΒ  has already been assigned to the Author Accepted Manuscript ver sion that might arise from this submission. We thank Adja Ndiaye forΒ  reviewing the manuscript.Β 

CONFLICT OF INTERESTΒ 

The authors declare that they have no conflict of interest.Β 

DATA AVAILABILITY STATEMENTΒ 

The data that support the findings of this study are available fromΒ  the corresponding author upon reasonable request.Β 

ORCIDΒ 

Monique Soro https://orcid.org/0000-0003-0459-5884 Fidèle Tiendrébéogo https://orcid.org/0000-0002-3619-3268 

REFERENCESΒ 

Abubakar, M., Singh, D. & Keta, J.N. (2019) Cassava mosaic disease andΒ  associated geminiviruses in Bauchi state, Nigeria: occurrence andΒ  distribution. American Journal of Plant Biology, 4, 85–90.Β 

Adjata, K.D., Muller, E., Peterschmi, M., Traore, O. & Gumedzoe, Y.M.D.Β  (2009) Molecular evidence for the association of a strain of UgandaΒ  variant of East African cassava mosaic virus to symptom severity inΒ  cassava (Manihot esculenta Crantz) fields in Togo. American JournalΒ  of Biochemistry and Biotechnology, 5, 196–201.Β 

Ally, H.M., Hamss, H.E., Simiand, C., Maruthi, M.N., Colvin, J., Omongo,Β  C.A. et al. (2019) What has changed in the outbreaking populationsΒ  of the severe crop pest whitefly species in cassava in two decades?Β  Scientific Reports, 9, 14796.Β 

Ariyo, O.A., Koerbler, M., Dixon, A.G.O., Atiri, G.I. & Winter, S. (2005)Β  Molecular variability and distribution of cassava mosaic begomovi ruses in Nigeria. Journal of Phytopathology, 153, 226–231.Β 

Bao, Y., Chetvernin, V. & Tatusova, T. (2014) Improvements to pairwiseΒ  sequence comparison (PASC): a genome-based web tool for virusΒ  classification. Archives of Virology, 159, 3293–3304.Β 

Benjamini, Y. & Yekutieli, D. (2001) The control of the false discoveryΒ  rate in multiple testing under dependency. Annals of Statistics, 29,Β  1165–1188.Β 

Berrie, L.C., Rybicki, E.P. & Rey, M.E.C. (2001) Complete nucleotide se quence and host range of South African cassava mosaic virus: fur ther evidence for recombination amongst begomoviruses. JournalΒ  of General Virology, 82, 53–58.Β 

Bock, K.R. & Woods, R.D. (1983) Etiology of African cassava mosaic dis ease. Plant Disease, 67, 994–995.Β 

Bull, S.E., Briddon, R.W., Sserubombwe, W.S., Ngugi, K., Markham, P.G. &Β  Stanley, J. (2006) Genetic diversity and phylogeography of cassavaΒ  mosaic viruses in Kenya. Journal of General Virology, 87, 3053–3065.Β 

Byju, G. & Suja, G. (2020) Mineral nutrition of cassava. Advances inΒ  Agronomy, 159, 169–235.Β 

Chikoti, P.C., Ndunguru, J., Melis, R., Tairo, F., Shanahan, P. & Sseruwagi,Β  P. (2013) Cassava mosaic disease and associated viruses inΒ  Zambia: occurrence and distribution. International Journal of PestΒ  Management, 59, 63–72.Β 

DabirΓ©, R. & Belem, J. (2003) Les plantes Γ  tubercules et racines auΒ  Burkina Faso. WASNET, 8, 12–16.Β 

El-Sharkawy, M.A. (2004) Cassava biology and physiology. PlantΒ  Molecular Biology, 56, 481–501.Β 

Eni, A.O., Efekemo, O.P., Onile-ere, O.A. & Pita, J.S. (2021) South WestΒ  and North Central Nigeria: assessment of cassava mosaic diseaseΒ  and field status of African cassava mosaic virus and East AfricanΒ  cassava mosaic virus. Annals of Applied Biology, 178, 466–479.Β 

Fondong, V.N., Pita, J.S., Rey, M.E.C., De Kochko, A., Beachy, R.N. &Β  Fauquet, C.M. (2000) Evidence of synergism between African cas sava mosaic virus and a new double-recombinant geminivirus infect ing cassava in Cameroon. Journal of General Virology, 81, 287–297.Β 

Guira, F., Some, K., Kabore, D., Sawadogo-Lingani, H., Traore, Y. &Β  Savadogo, A. (2017) Origins, production, and utilization of cassavaΒ  in Burkina Faso, a contribution of a neglected crop to householdΒ  food security. Food Science and Nutrition, 5, 415–423.Β 

KonatΓ©, G., Barro, N., Fargette, D., Swanson, M.M. & Harrison, B.D.Β  (1995) Occurrence of whitefly-transmitted geminiviruses in cropsΒ  in Burkina Faso and their serological detection and differentiation.Β  Annals of Applied Biology, 126, 121–130.Β 

Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. (2018) MEGA X: mo lecular evolutionary genetics analysis across computing platforms.Β  Molecular Biology and Evolution, 35, 1547–1549.Β 

Legg, J.P., French, R., Rogan, D., Okao-Okuja, G. & Brown, J.K. (2002)Β  A distinct Bemisia tabaci (Gennadius) (Hemiptera: Sternorrhyncha:Β  Aleyrodidae) genotype cluster is associated with the epidemic ofΒ  severe cassava mosaic virus disease in Uganda. Molecular Ecology,Β  11, 1219–1229.

2216 |    SORO et al.Β 

Legg, J.P., Owor, B., Sseruwagi, P. & Ndunguru, J. (2006) Cassava mosaicΒ  virus disease in East and Central Africa: epidemiology and man agement of a regional pandemic. Advances in Virus Research, 67,Β  355–418.Β 

MacFadyen, S., Paull, C., Boykin, L.M., De Barro, P., Maruthi, M.N.,Β  Otim, M. et al. (2018) Cassava whitefly, Bemisia tabaci (Gennadius)Β  (Hemiptera: Aleyrodidae) in East African farming landscapes: a re view of the factors determining abundance. Bulletin of EntomologicalΒ  Research, 108, 565–582.Β 

Maruthi, M.N., Seal, S., Colvin, J., Briddon, R.W. & Bull, S.E. (2004) EastΒ  African cassava mosaic Zanzibar virus – a recombinant begomo virus species with a mild phenotype. Archives of Virology, 149,Β  2365–2377.Β 

Mugerwa, H., Colvin, J., Alicai, T., Omongo, C.A., Kabaalu, R., Visendi, P.Β  et al. (2021) Genetic diversity of whitefly (Bemisia spp.) on crop andΒ  uncultivated plants in Uganda: implications for the control of thisΒ  devastating pest species complex in Africa. Journal of Pest Science,Β  94, 1307–1330.Β 

Mulenga, R.M., Legg, J.P., Ndunguru, J., Miano, D.W., Mutitu, E.W.,Β  Chikoti, P.C. et al. (2016) Survey, molecular detection, and charac terization of geminiviruses associated with cassava mosaic diseaseΒ  in Zambia. Plant Disease, 100, 1379–1387.Β 

Mwatuni, F., Ateka, E., Karanja, L., Mwaura, S. & Obare, I. (2015)Β  Distribution of cassava mosaic geminiviruses and their associ ated DNA satellites in Kenya. American Journal of ExperimentalΒ  Agriculture, 9, 1–12.Β 

Ntawuruhunga, P., Dixon, A.G.O., Kanju, E., Ssemakula, G., Okechukwu,Β  R.U., Whyte, J.B.A. et al. (2013) Successful innovations and les sons learnt in cassava improvement and deployment by IITA in theΒ  Eastern African Region. African Journal of Root and Tuber Crops, 10,Β  41–54.Β 

Ogbe, F.O., Dixon, A.G.O., Hughes, J.D’A., Alabi, O.J. & Okechukwu,Β  R. (2006) Status of cassava begomoviruses and their new naturalΒ  hosts in Nigeria. Plant Disease, 90, 548–553.Β 

Patil, B.L. & Fauquet, C.M. (2009) Cassava mosaic geminiviruses: ac tual knowledge and perspectives. Molecular Plant Pathology, 10,Β  685–701.Β 

Permingeat, H.R., Romagnoli, M.V., Juliana, I. & Vallejos, R.H. (1998) AΒ  simple method for isolating DNA of high yield and quality from cot ton (Gossypium hirsutum L.) leaves. Plant Molecular Biology Reporter,Β  16, 89.Β 

Pita, J.S., Fondong, V.N., SangarΓ©, A., Kokora, R.N.N. & Fauquet, C.M.Β  (2001) Genomic and biological diversity of the African cassavaΒ  geminiviruses. Euphytica, 120, 115–125.Β 

Pita, J.S., Fondong, V.N., SangarΓ©, A., Otim-Nape, G.W., Ogwal, S. &Β  Fauquet, C.M. (2001) Recombination, pseudorecombination andΒ  synergism of geminiviruses are determinant keys to the epidemicΒ  of severe cassava mosaic disease in Uganda. Journal of GeneralΒ  Virology, 82, 655–665.Β 

Price, M.N., Dehal, P.S. & Arkin, A.P. (2010) FastTree 2 – approximatelyΒ  maximum-likelihood trees for large alignments. PLoS One, 5, e9490.Β 

Pushpalatha, R. & Gangadharan, B. (2020) Is cassava (Manihot esculenta Crantz) a climate β€œsmart” crop? A review in the context of bridgingΒ  future food demand gap. Tropical Plant Biology, 13, 201–211.Β 

Saediman, H., Limi, M.A., Rosmawaty, A.P. & Indarsyih, Y. (2016) CassavaΒ  consumption and food security status among cassava growingΒ  households in southeast Sulawesi. Pakistan Journal of Nutrition, 15,Β  1008–1016.Β 

Sseruwagi, P., Sserubombwe, W.S., Legg, J.P., Ndunguru, J. & 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 Research, 100, 129–142.Β 

Stanley, J. & Gay, M. (1983) Nucleotide sequence of cassava latent virusΒ  DNA. Nature, 301, 260–262.Β 

Terry, E.R. (1975) Description and evaluation of cassava mosaic diseaseΒ  in Africa. In: Terry, E.R. & MacIntyre, R. (Eds.) The international ex change and testing of cassava germplasm in Africa. Ibadan, Nigeria:Β  IITA, pp. 53–54.Β 

TiendrΓ©bΓ©ogo, F., Lefeuvre, P., Hoareau, M., Harimalala, M.A., De Bruyn,Β  A., Villemot, J. et al. (2012) Evolution of African cassava mosaicΒ  virus by recombination between bipartite and monopartite bego moviruses. Virology Journal, 9, 1–7.Β 

TiendrΓ©bΓ©ogo, F., Lefeuvre, P., Hoareau, M., TraorΓ©, V.S.E., Barro, N.,Β  Reynaud, B. et al. (2009) Occurrence of East African cassava mo saic virus -Uganda (EACMV-UG) in Burkina Faso. Plant Pathology,Β  58, 783.Β 

Torkpo, S.K., Gafni, Y., Danquah, E.Y. & Offei, S.K. (2018) Incidence andΒ  severity of cassava mosaic disease in farmer’s fields in Ghana.Β  Ghana Journal of Agricultural Science, 53, 61.Β 

Torkpo, S.K., Offei, K., Danquah, E.Y. & Gafni, Y. (2017) Status of cassavaΒ  mosaic begomoviruses in farmers’ fields in Ghana. AIMS AgricultureΒ  and Food, 2, 279–289.Β 

Toualy, M.N.Y., Akinbade, S.A., Koutoua, S., Diallo, H.A. & Lava, P.K.Β  (2014) Incidence and distribution of cassava mosaic begomovirusesΒ  in CΓ΄te d’Ivoire. International Journal of Agronomy and AgriculturalΒ  Research, 4, 131–139.Β 

Zhou, X., Robinson, D.J. & Harrison, B.D. (1998) Types of variationΒ  in DNA-A among isolates of East African cassava mosaic virusΒ  from Kenya, Malawi and Tanzania. Journal of General Virology, 79,Β  2835–2840.Β 

Zinga, I., Chiroleu, F., Legg, J., Lefeuvre, P., Komba, E.K., Semballa, S.Β  et al. (2013) Epidemiological assessment of cassava mosaic diseaseΒ  in Central African Republic reveals the importance of mixed viral in fection and poor health of plant cuttings. Crop Protection, 44, 6–12.Β 

How to cite this article: Soro, M., TiendrΓ©bΓ©ogo, F., Pita, J.S.,Β  TraorΓ©, E.T., SomΓ©, K., Tibiri, E.B., et al (2021) EpidemiologicalΒ  assessment of cassava mosaic disease in Burkina Faso. PlantΒ  Pathology, 70, 2207–2216. https://doi.org/10.1111/ppa.13459