1. Introduction
Cassava (
Manihot esculenta Crantz) is a highly significant root crop, widely utilized for human consumption, animal feed, and industrial applications [
1]. The starchy storage roots of cassava have emerged as a crucial source of dietary energy in sub-Saharan Africa (SSA), offering higher returns per unit of input compared to other staple crops [
1,
2]. Cassava exhibits exceptional resilience, thriving in nutrient-poor soils and displaying superior yield performance in comparison to other root and tuber crops [
2,
3]. Nevertheless, the response of cassava genotypes to diverse environmental factors (such as soil and climate) and biotic stresses can vary significantly [
3,
4].
In Sierra Leone, the lack of comprehensive data regarding the genetic variation within the cassava breeding population hampers the development of superior cassava genotypes. The generation of such data involves the use of genetic markers. Genetic markers serve as an effective tool for evaluating genetic divergence within the breeding population [
5]. Preliminary studies often utilize genetic markers, including agro-morphological markers, due to their efficiency and simplicity in assessing the diversity of the cassava germplasm [
6]. Certain morphological traits have proven to accurately reflect the diversity perceived by farmers [
7,
8]. Additionally, Elias et al. [
9] reported that morphological traits exhibit heritable genetic variation. Morphological traits, both quantitative and qualitative, have been widely employed in the systematic identification of genotypes, species, and genera in various crops [
10]. Qualitative traits are typically governed by a few major genes, resulting in distinct observable differences that facilitate genotype differentiation and identification. On the other hand, quantitative traits are influenced by numerous minor genes with complex inheritance patterns. These traits are more susceptible to environmental effects during the developmental stage of the crop. Genetic diversity studies based solely on morphological traits can be constrained by environmental factors and genotype–environment interactions [
11]. These limitations may hinder the accurate identification of duplicate genotypes through morphological classification techniques alone. To overcome these challenges, Collard et al. [
11] highlighted the use of molecular markers, which enable the detection of genetic distinctions among closely related genotypes.
Assessing the extent of genetic variation within the breeding population is essential for identifying valuable genetic divergence that is crucial for improving cassava populations [
12]. Knowledge regarding the diversity of cassava accessions and their agronomic traits could be useful for breeding programs to develop high-yielding and disease-resistant varieties with improved culinary traits [
13]. Therefore, the characterization of the cassava germplasm can help to identify and promote the use of diverse cassava varieties with different agronomic and culinary traits. This can improve the resilience of cassava production systems in Sierra Leone by reducing the risk of crop failure due to pest and disease outbreaks or changes in environmental conditions. The agro-morphological and culinary trait characterization of the cassava germplasm can provide breeders with information on the genetic diversity of cassava varieties in Sierra Leone. This information can be used to identify desirable traits, such as high yield, disease resistance, and culinary qualities, which can be incorporated into new breeding programs. This study assessed the genetic variability and relationships within 103 cassava genotypes using agro-morphological and culinary markers.
3. Results
3.1. Correlation Analysis, Principal Component Analysis and Genetic Relationships among Cassava Accessions Based on 25 Qualitative Agro-Morphological Traits
The results of the correlation matrices of 25 qualitative characteristics revealed a significant and positive correlation between color of leaf veins and leaf color (r = 0.22 *), color of leaf veins and petiole color (r = 0.62 ***), leaf color and color of apical leaves (r = 0.39 ***), leaf color and petiole color (r = 0.43 ***), leaf color and leaf retention (r = 0.22 *), color of apical leaves and leaf retention (r = 0.23 *), shape of center leaflet and color of stem epidermis (r = 0.23 *), leaf margin and presence of foliar scars (r = 0.29 **), leaf margin and color of end branches (r = 0.24 *), petiole color and leaf retention (r = 0.24 *), root taste and presence of fruit (r = 0.25 *), color of stem epidermis and plant shape (r = 0.21 *), leaf retention and stipule margin (r = 0.24 *), height at first branching and presence of fruit (r = 0.32 **), height at first branching and seed (r = 0.31 **), plant shape and presence of fruit (r = 0.47 ***), plant shape and seed (r = 0.44 ***), and presence of fruit and seeds (r = 0.86 ***) (
Table 3). Conversely, significant and negative correlations were noted between leaf color and the shape of the center leaflet (r = −0.30 ***), leaf color and branching habit (r = 0.20 *), the shape of the center leaflet and petiole color (r = −0.42 ***), the shape of the center leaflet and root constriction (r = −0.29 **), petiole color and color of stem exterior (r = −0.25 *), root shape and color of root cortex (r = −0.24 *), color of stem cortex and leaf retention (r = −0.27 **), color of stem exterior and presence of seeds (r = −0.24 *), orientation of petiole and plant shape (r = −0.23 *), stipule margin and root constriction (r = −0.19 *), and branching habit and plant shape (r = −0.53 ***).
The eigenvalues and percentage variations of the principal components are presented in
Table 4. Ten principal components (PCs), which accounted for 66.744% of the total variation among the genotypes were identified. The first PC axis, with an eigenvalue of 2.971, accounted for 11.885% of the total variation, while the second, third, and fourth PCs, with eigenvalues of 2.611, 1.774, and 1.651, accounted for 10.444%, 7.098%, and 6.604% of the total variation, respectively. The fifth, sixth, and seventh PCs, with eigenvalues of 1.577, 1.501, and 1.177, accounted for 6.309%, 6.004%, and 4.707% of the total variation, respectively, while the eighth, ninth, and tenth PCs, with eigenvalues of 1.167, 1.156, and 1.1, accounted for 4.669%, 4.625%, and 4.4% of the total variation, respectively. The first principal component was positively associated with root taste, orientation of petiole, branching habit, presence of fruit, and presence of seeds. The second PC, with the second highest loading factor, was positively associated with color of leaf vein, leaf color, color of apical leaf, and petiole color. The third PC was associated with root taste, color of apical leaf, leaf retention, and orientation of petiole, while the fourth PC was associated with leaf margin, root cortex, and color of end branches. Traits that were positively associated with the fifth PC were leaf color, color of apical leaf, color of root cortex, color of stem exterior, and presence of fruit and seeds, while leaf retention, color of end branches, and stipule margin were the traits was associated with the sixth PC and the seventh PC was positively associated mainly with storage root characteristics (root taste and external color of storage). The eighth and ninth PCs were positively associated with root taste, root shape, external color of storage and leaf margin, color of root cortex and leaf retention, respectively, while petiole color and prominence of foliar scale were the only traits that were associated positively with the tenth PC.
The hierarchical classification of 103 cassava genotypes based on qualitative traits grouped the genotypes into three classes with similar characteristics as a function of the variable. The genetic similarity for the 25 qualitative traits ranged from zero to two, with a mean similarity of 2.0. The cassava genotypes were grouped into three distinct clusters with similarities of 1.4. Clusters I and III have a higher number of genotypes, with 40 and 34, respectively, while 29 individuals were in cluster II (
Figure 1).
3.2. Descriptive Statistics, Correlation Analysis, Principal Component Analysis, and Genetic Relationships among Cassava Accessions Based on 14 Quantitative Agro-Morphological Traits
The cassava genotypes showed variability for fourteen quantitative morphological traits assessed (
Table 5). The ranges of the values produced were 1.0 to 5.0 cm for distance between leaf scars, 0.0 to 165.9 cm for height at first branching, 0.16 to 0.68 for harvest index, 0.0 to 20.0 for level of branching, 5.7 cm to 21.8 cm for length of leaf lobes, 1.0 cm to 5.0 cm for length of stipule, 1.0 to 24.0 for number of commercial roots, 3.0 to 9.0 for number of leaf lobes, 2.0 to 53.0 for number of storage roots, 41.0 cm to 346 cm for plant height, 4.7 cm to 37.3 cm for petiole length, 15.0% to 36.0% for root dry matter content, 1.0 kg to 19.0 kg for root yield per plant, and 1.5 cm to 6.0 cm for width of leaf lobes (
Table 5). The coefficients of variation varied from 14.7% (root dry matter content) to 81.2% (level of branching). Based on the 14 quantitative characteristics/traits, 11 had high (CV > 20%) coefficients of variation (
Table 5). Only three characteristics had low variation, including the root dry matter content (14.8%), the distance between leaf scars (18.7%), and the number of leaf lobes (17.8%).
Phenotypic correlations among 14 quantitative traits revealed that harvest index is the principal trait that is significantly and positively correlated with root yield per plant (r = 0.40 ***) and number of storage roots (r = 0.24 *) (
Table 6). The root yield per plant is highly significant and positively correlated with the number of storage roots (r = 0.76 ***) and the number of commercial roots (r = 0.68 ***). The petiole length is significantly and positively correlated with the length of leaf lobes (r = 0.22 *) and the number of leaf lobes (r = 0.26 **). The length of leaf lobes is significantly and positively correlated with the width of leaf lobes (r = 0.66 ***) and the number of leaf lobes (r = 0.41 ***); the width of leaf lobes is also positively correlated with plant height (r = 0.20 *).
The number of storage roots was significantly and positively correlated with the number of commercial roots (r = 0.84 ***), while the plant height was positively correlated with the number of leaf lobes (r = 0.24 *), and the height at first branching was positively correlated with the level of branching (r = 0.25 *). Conversely, significant and negative correlations were noted between the harvest index and the distance between leaf scars (r = 0.20 *), length of leaf lobes, and level of branching (r = 0.21 *), and between the width of leaf lobes and the level of branching (r = 0.21 *).
The principal component analysis of quantitative agro-morphological traits revealed that the six main principal components accounted for 67.72% of the total variation among the genotypes (
Table 7). The first factorial plane contains 18.34% of the variance. The traits that significantly correlated with axis 1 were harvest index (43.0%), root yield per plant (86.0%), number of storage roots (90.0%), and number of commercial roots (85.0%). The variables that were significantly correlated with axis 2 were petiole length (43.0%), length of leaf lobes (83.0%), width of leaf lobes (70.0%), plant height (40.0%), and number of leaf lobes (53.0%).
The traits that were significantly related to axis 3 were petiole length (44.0%), height at first branching (35.0%), and level of branching (30.0%). The variables significantly correlated to axis 4 were plant height (30.0%), height at first branching (66.0%), and level of branching (61.0%). The variables significantly related to axis 5 were harvest index (54.0%), number of leaf lobes (44.0%), and length of stipule (32.0%), and the traits that were significantly correlated to axis 6 were root dry matter content (54.0%) and length of stipule (44.0%) (
Table 7).
Hierarchical classification of 14 quantitative agro-morphological traits grouped the genotypes into three classes with similar characteristics as a function of the variable (
Figure 2). The genetic similarity for the even quantitative traits ranged from zero to two, with a mean similarity of 1.8. Cluster I contains 53 genotypes, cluster II contains six genotypes, and cluster III contains 44 genotypes.
3.3. Culinary Traits and Marketable Leaf Assessment of 15 Selected Cassava Genotypes
The mean cooking time (CT) ranged from 15.2 to 30.0 min (
Table 8). The genotypes that had shorter and better cooking times were Cookson (15.2 min) and Butter cassava (15.4 min), while those with the longest cooking times were Pink lady (30.0 min), SLICASS 4 (30.0 min), and SLICASS 6 (30 min). The cooking percentage was included as an important culinary trait due to the large variation that is normally observed between the cooking time of cassava varieties. The mean cooking percentage (CP) ranged from 27.9% to 96.4%. The mean percentages were classified into two groups. The group with lowest membership included the genotypes SLICASS 4 (27.9%), SLICASS 6 (30.2%), and SLICASS 7 (58.3%). The means of the other groups of genotypes ranging from 63.2% to 96.4% were Pink lady (63.2%), Soja color (70.3%), Kendemeh (73.4%), Tangagboi (76.2%), Mende tangai (78.4%), Yaa kanu (80.2%), Nikaneh (80.2%), Tapiyoka (82.0%), Cocoa cassada (83.2%), Ndiamonyamalo (88.3%), Butter cassava (95.0%), and Cookson (96.4%).
Mean scores for surface appearance, mealiness, taste, and aroma ranged from 2.0 to 4.75, and for texture, from 2.0 to 4.50, indicating the acceptable mean threshold (
Table 8). Comparing the means for sensory characteristics, genotypes that had the highest mean scores for surface appearance were Cookson (4.50 ± 0.50), Butter cassava (4.50 ± 0.58), and Cocoa cassava (4.50 ± 0.50), closely followed by Ndiamonyamalo (4.3 ± 0.50), Nikaneh (4.3 ± 0.58), Tapiyoka (4.3 ± 0.55),Tangagboi (4.25 ± 0.58), Mende Tangai (4.0 ± 0.65), Yaa kanu (4.0 ± 0.50), and SLICASS 7 (4.0 ± 0.50), which clearly indicate good acceptability by consumers for these traits, while other genotypes like Kendemeh (3.5 ± 0.58), Pink lady (3.5 ± 0.82), SLICASS 6 (3.5 ± 0.58) and Soja color (3.5 ± 0.50) were fairly acceptable to consumers. SLICASS 4 has the lowest mean score (2.0 ± 0.40) for surface appearance, thus indicating poor acceptability to consumers/farmers (
Table 8).
The results for mealiness revealed that Cookson has the highest mean score (4.75 ± 0.80), followed by Butter cassava (4.50 ± 0.50), Tapiyoka (4.3 ± 0.50), Nikaneh (4.3 ± 0.58), Ndiamonyamalo (4.25 ± 0.48), Cocoa cassava (4.25 ± 0.80), Tangagboi (4.0 ± 0.40), Mende tangai (4.0 ± 0.65), Yaa kanu (4.0 ± 0.60) and Kendemeh (4.0 ± 0.63), thus indicating good acceptability of the genotypes by consumers/farmers for mealiness. Pink lady (3.5 ± 0.82), SLICASS 6 (3.5 ± 0.58), and Soja color (3.0 ± 0.50) were fairly acceptable to farmers. The genotype that recorded the lowest mean score in terms of mealiness was SLICASS 4 (2.0 ± 0.82), indicating poor acceptability by consumers of the genotype for this trait. Cookson was highly rated for its taste and thus records the highest mean score (4.75 ± 0.50), closely followed by Butter cassava (4.50 ± 0.43), Ndiamonyamalo (4.50 ± 0.50), Cocoa cassava (4.50 ± 0.58), Tapiyoka (4.50 ± 0.80),Tangagboi (4.50 ± 0.50), Nikaneh (4.3 ± 0.58), Yaa kanu (4.25 ± 0.50), Kendemeh (4.25 ± 0.50), and Mende tangai (4.0 ± 0.65), indicating good acceptability to farmers, while Pink lady (0.5 ± 0.65) and SLICASS 7 (3.0 ± 0.45) were fairly acceptable to farmers for this trait. However, SLICASS 4 (2.0 ± 0.50), SLICASS 6 (2.0 ± 0.50), and Soja color (2.0 ± 0.50) were rated poorly for taste quality. Sensory evaluation on the aroma of the cooked roots of the genotypes revealed that Cookson (4.75 ± 0.50) and Butter cassava (4.75 ± 0.50) scored the highest mean followed by Ndiamonyamalo (4.3 ± 0.50), Tapiyoka (4.3 ± 0.50), Cocoa cassava (4.25 ± 0.58), Tangagboi (4.25 ± 0.58), Nikaneh (4.25 ± 0.50), Yaa kanu (4.25 ± 0.82), Mende tangai (4.0 ± 0.65), Kendemeh (4.0 ± 0.58), and SLICASS 7 (4.0 ± 0.60), which indicate good acceptability by consumers, while Soja color (3.0 ± 0.65), Pink lady (3.0 ± 0.50), and SLICASS 6 (3.0 ± 0.50) were fairly acceptable to consumers in terms of aroma. However, SLICASS 4 was poorly accepted by consumers, with a mean score of 2.0 ± 0.82. Cookson and Butter cassava scored the highest means of 4.50 ± 0.58 and 4.50 ± 0.50, respectively, in terms their texture, followed by Ndiamonyamalo (4.3 ± 0.50), Cocoa cassava (4.25 ± 0.58), Tangagboi (4.25 ± 0.58), Mende Tangai (4.0 ± 0.65), Nikaneh (4.0 ± 0.50) and Tapiyoka (4.0 ± 0.60), which indicate good acceptability to consumers. Other genotypes like Yaa kanu (3.5 ± 0.58), Kendemeh (3.5 ± 0.50), Pink lady (3.0 ± 0.58), SLICASS 7 (3.0 ± 0.50), and SLICASS 6 (3.0 ± 0.50) were fairly acceptable to consumers. SLICASS 4 (2.0 ± 0.82) and Soja color (2.0 ± 0.82) were poorly accepted.
3.4. Marketable Leaf Assessment
The five elite varieties selected by farmers based on the desired leaf traits for selling and consumption included Tangaigboi, Coco cassava, Cookson, Butter cassava, and Ndiamonyamalo, followed by the moderately preferred varieties Tapiyoka, Yaa kanu, Nikaneh, Mende tangai, and Kendemeh (
Table 9). However, the leaves of the varieties SLICASS 4, SLICASS 6, SLICASS 7, Pink lady, and Soja color are not usually consumed and sold.
4. Discussion
Significant variation was noted from the analysis of the qualitative morphological traits (root taste, external color of storage roots, color of root pulp, ease of peeling, color of leaf veins, lobe margins, leaf color, color of apical leaves, shape of central leaflets, leaf color, color of root cortex, petiole color, root cortex, color of stem exterior, presence of fruit and seeds, leaf retention, color of end branches and stipule margin, and external color of storage) among the studied genotypes. The most representative and distinctive trait was color, possibly due to the fact that most of the genotypes exhibited white root pulp and dark brown external storage roots. The aboveground leaf traits of the studied genotypes were green leaf veins, light green leaves, light green apical leaves, smooth lobe margins, and elliptic-lanceolate center leaflets. The leaf attributes play significant roles in cultivar identification and, more significantly, in cassava selection by farmers for leafy vegetable markets. The findings agree with the work of Karim et al. [
13], who reported that leaf attributes are important for cassava identification and selection for leafy vegetable markets in Sierra Leone, where cassava leaves are widely consumed. The findings also agree with Agre et al. [
20], who reported that farmers use the color of the leaves and stems to identify their cassava cultivars.
Principal component analysis is a powerful data reduction technique utilized to reduce large numbers of correlated variables to a small number that is independent and very useful. The PCA revealed the traits that contributed most to the variation present in the cassava germplasm. Qualitative traits with the highest positive contribution to the first PCA included root shape, orientation of petioles, branching habits, and presence of fruit and seeds. The findings of the study revealed the usefulness of these traits for genotype identification and genetic diversity studies in cassava. The high positive values of phenotypic correlations indicate that the positive selection of one of these traits results in increases in the other trait, whereas the reverse is also true. Traits with high correlation coefficients are among the key traits often considered relevant for the selection of varieties for the genetic improvement of the crop. The findings concur with the view that PCA results reveal a trend towards higher genetic values for some traits with significant strong positive or negative contributions to the observed genetic variability compared to the contributions of those with weak correlation coefficients [
21]. Thus, the PCA results of the present study demonstrated a pattern of higher phenotypic values for some traits with significant strong positive or negative contributions to the observed phenotypic variability compared to the contributions of those with weak correlation coefficients. The trait variations reflected among the genotypes could be attributable to their different genetic backgrounds.
The clustering based on the similarity index of the qualitative traits in this study grouped the 103 cassava accessions into three clusters. Cluster I contained the accessions characterized by green apical leaves, light green leaves, green leaf veins, yellowish-green petioles, umbrella plant shapes, and average leaf retention. Cluster II included accessions with elliptic-lanceolate shaped leaflets, light brown stem epidermises, light brown stem exteriors, ease of peeling, and sweet root taste. Cluster III included accessions based on the absence of fruit, absence of seeds, dichotomous branching habits, horizontal petiole orientation, dark brown external storage roots, light green leaves, and ease of peeling. A similar study conducted by Karim et al. [
13] using 102 cassava accessions identified five distinct groups. In this research, the first two principal components explained 22.33% of the total cumulative variance for the qualitative traits. This result underscored the findings of Karim et al. [
13] and Afonso et al. [
22], who found 31.18% and 32.56% of the genetic variance in the first factorial plane, respectively. It can also be explained by the fact that the variance distribution is associated with the nature and number of characteristics used in the analysis and focuses on the first principal components.
The 103 cassava accessions also showed variability for the 14 quantitative morphological traits assessed. In order to appreciate the variability of each trait, all the quantitative traits recorded were subjected to descriptive statistical analysis (minimum, maximum, mean, variance, standard deviation, and coefficient of variation). This study revealed high variability and higher coefficients of variation for most of the studied traits, which indicated the presence of high heterogeneity within the characterized population, which could be exploited for future breeding. These results are similar to those reported by Agre et al. [
23] and Karim [
24]. Agre et al. [
23] found higher coefficients of variation for 75% of the traits studied.
The analysis of the first six principal components explained 67.72% of the overall variability in the quantitative analysis. Principal component I presented yield and yield attribute traits such as harvest index, root yield per plant, number of storage roots, and number of commercial roots, while principal components II, III, and IV presented growth attribute traits such as petiole length, length of leaflets, width of leaflets, plant height, height at first branching, and level of branching. Principal components V and VI presented both yield and growth attribute traits such as harvest index, root dry matter content, and length of stipule. The quantitative traits with the highest positive contribution to the distinguishing accessions in the first PCA included harvest index, root yield per plant, number of storage roots, and number of commercial roots. These traits are among the key traits often considered for selection of varieties for the genetic improvement of cassava. This study agreed with the work of Karim et al. [
13], who observed that the first principal component presented yields and attribute traits such as harvest index, average yield per plant, number of storage roots, root dry matter content, and starch content, and they were the highest positive contributors to the principal component analysis. The cluster analysis of the 14 quantitative agro-morphological traits also grouped the accessions into three groups. Cluster I comprised accessions with high numbers of storage roots, numbers of commercial roots, root yields per plant, harvest indexes, and numbers of leaf lobes; cluster II included accessions based on high petiole length and root dry matter content; while cluster III accessions exhibited high numbers of commercial roots, numbers of storage roots, root yields per plant, width of leaf lobes, and harvest indexes. The results of this study demonstrate the relevance of aboveground and root attributes in characterizing cassava accessions. This study also demonstrated the usefulness of Fukuda et al.’s [
14] agro-morphological descriptor in identifying variability and reducing dimensionality in the trait set. In this study, the 25 qualitative and 14 quantitative trait sets sufficiently discriminated the 103 genotypes into distinct cluster groups.
Root tenderness after cooking is a relevant attribute for cassava since dry matter content in cassava sometimes plays an important role. The results revealed that only the Cookson and butter cassava accessions recorded cooking times below 20 min, which is considered acceptable for cassava according to the work of De Pedri et al. [
25]. These authors further reported that the easier the cooking, the better the taste of the cooked product, indicating the relevance of dry matter content for acceptable cassava taste. Due to the importance of cooking sweet cassava and the variations in the time of cooking that usually exist between cassava accessions, it is important to consider the cooking percentage of each accession. Among all the accessions evaluated for this trait, Cookson and butter cassava had higher and better cooking percentages than other accessions. This result implied that the shorter the cooking time, the higher the cooking percentage. Thus, Cookson and butter cassava had shorter cooking times and higher cooking percentages than other accessions. SLICASS 4 and SLICASS 6 had low cooking percentages during the 30 min maximum cooking time set. This result confirms the work of Santos et al. [
26], who reported wide variations in the cooking percentages of 11 cassava clones evaluated at different harvesting times. The sensory evaluation results indicated that 60% of the total accessions evaluated had good acceptable cooking qualities (surface appearance, mealiness, taste, aroma, and texture). Only one accession (SLICASS 4) was poorly rated by farmers for all the culinary traits. This variety is an improved variety that was developed by the Sierra Leone Agricultural Institute (SLARI). The variety is high yielding but lacks acceptable cooking qualities (mealiness, taste, and texture). The findings imply the significance of the incorporation of desired end-user traits in cassava breeding programs. Similar research conducted by Mugalavai et al. [
27] revealed that 18.0% of the 51 accessions evaluated had acceptable cooking qualities, such as texture, taste, and mealiness.
5. Conclusions
This study successfully determined the extent of genetic divergence within the cassava germplasm of Sierra Leone using morphological markers. It also provides vital information to help cassava scientists to make informed decisions for parental selection for morphological and culinary traits based on genetic diversity. The useful genetic variability in root dry matter content, number of storage roots, number of commercial roots, root yield per plant, harvest index, length of leaflets, width of leaflets, height at first branching, distance between leaf scars, and level of branching that were identified could be exploited for the genetic improvement of the crop and its conservation. The root shape, orientation of petiole, branching habit, presence of fruit and seeds, and color attributes of various qualitative traits studied contributed most to the differentiation of the genotypes. The agro-morphological and culinary traits were useful in distinguishing and selecting cassava genotypes, respectively. Both approaches should therefore be used for genetic diversity studies of cassava.