Field Survey
Surveys were conducted throughout the regions indicated in Figure 1.They were carried out in September and October of 2015 and December 2016 to January 2017, using a harmonized sampling protocol developed by the West African Virus Epidemiology (WAVE) for root and tuber project head quartered at the University of Félix Houphouët-Boigny, Abidjan, Côte d’Ivoire. During each survey, coordinates of fields visited were captured using a GPS device. The minimum distance between two fields where samples were collected was 10 Km and in areas where cassava production was not widespread, the distance between two fields varied depending on how fields were encountered. In each field transect walks along two diagonals were made, and disease incidence and severity were assessed on 15 plants within each transect. Thus, a total of 30 plants were assessed in each field for disease incidence and severity. Disease incidence was measured as a percentage of number of plants infected per field based on the sampled plants while disease severity was assessed based on a scale ranging from 1 to 5 as defined by the International Institute of Tropical Agriculture (IITA). Using this scale, 1 represents absence of infection; 2: mild infection; 3: moderate infection; 4: severe infection and 5 represents very severe infection (IITA 1990). Data were collected via an electronic device application developed by the University of Cambridge, UK. The data were then uploaded into a central repository at the University of Cambridge, UK for analysis. From each field, leaf samples showing symptoms of cassava mosaic virus infection were collected; these were kept in herbarium pressers and labeled. In addition, stem cuttings of infected plants were also collected and labeled. The leaf samples were stored in the laboratory under ambient conditions while the stem cuttings were maintained in insect proof screenhouses prior to their laboratory analysis. Whitefly Bemisia tabaci vectors were collected using an aspirator from 5 plants/field and the total number counted and recorded. The whiteflies were preserved in 70% ethanol, labeled and preserved at the laboratory in a freezer. The diversity of cassava genotypes encountered on the fields were differentiated using the leaf petiole colour as most of the farm owners could not readily give the names of the varieties they had planted.
Laboratory Diagnostics
Genomic DNA was extracted from samples based on cetyl trimethylammonium bromide (CTAB) method (Dellaporta et al., 1983).
Polymerase chain reaction (PCR)
PCR was carried out in a 25 µl reaction consisting of the following reagents: PCR-grade water; 5 µl of 5x standard buffer; 2.5 µl 5% Tween-20; 0.25 µl of 10 mM dNTPs; 0.25 µl each of Forward and Reverse primers; 5 U Taq polymerase and 5 µl of DNA template. The list of primers used for the analysis can be found in Table 1.
Thermocycling conditions
Amplification conditions included a first PCR cycle comprising denaturation at 94 ˚C for 5 min, annealing of primers at 52 oC for 1 min and elongation at 72 oC for 2 min. This initial amplification cycle was followed by 35 cycles of 1 min at 94 oC, 1 min at 52 oC and 2 min at 72 oC. At the end of the reaction, a final elongation step was achieved at 72 oC for 10 min and the reaction hold at 4o.
PCR products were separated by electrophoresis on a 1.5 % agarose gel containing ethidium bromide (1 lg⁄10 ml) under a constant current of 100 V. Visualization of the amplified bands under UV light was done using a Vilber Lourmat (infinityTM) gel documentation imaging system.
Sequencing
Where the expected bands on the agarose gel were slightly different rom the expected band, recovered PCR products were prepared for sequencing. Sequencing was carried out by GENEWIZ, Inc (South Plainfield, NJ) using Applied Biosystems Big Dye version 3.1. The reactions were then run through an Applied Biosystem 3730xl DNA Analyzer for DNA sequencing.
Data Analysis
For our analysis, we made use of data-analysis software, namely a multi-dimensional data cube referred to, within the WAVE programme, as the WAVE Cube. This tool enabled us to interrogate and view our survey data from different perspectives to present our extracted results as tables and graphics.
Nucleotide sequences obtained from GENEWIZ were used as query and BLAST searched with ACMV and EACMV strains/isolate nucleotide sequences stored in NCBI GenBank database to detect their corresponding strains/isolates. Nucleotide sequence of samples that yielded above 80% identity with those of ACMV and EACMV isolates in the NCBI GenBank database repository were considered similar (Fondong et al., 2000).
Patristic distances between pairs of sequences analyses were conducted using the Maximum Composite Likelihood model (Tamura et al., 2004). This analysis involved 56 nucleotide sequences. All ambiguous positions were removed for each sequence pair (pairwise deletion option). Evolutionary analyses were conducted in MEGA X (Kumar et al., 2018).