Abstract

Cassava is a significant contributor to food security and an income source for smallholder farmers in southern Ethiopia. However, little research effort has been done so far based on designing field experiment samples for the biochemical composition of cassava accession at the country level. The study was conducted to assess the biochemical composition of cassava accessions in southwest Ethiopia. Flour samples from the storage roots of 64 cassava accessions were collected and were run in duplicates. Data on 13 biochemical characters were collected and analyzed using standard methods. The analysis of variance showed significant to very highly significant differences among the tested accessions for biochemical composition. The flour moisture ranged from 4.83–10.11%, dry matter (89.89-95.17%), organic matter (86.71–92.65%), ash (2.1–3.96%), fiber (1.14–3.00%), fat (0.26-1.4%), crude protein (1.28-2.86%), starch (65.1–74.2%), carbohydrate (81.29–87.94%), energy (341.44–367.61 kcal/100g DM), and cyanide (1.67–3.14). The highest GCV = 29.54% was shown for crude fat, followed by GCV = 16.94% for crude fiber, and GCV = 16.11% for tannin, whereas, among the characters, dry matter was observed to be the lowest (GCV = 0.84%). The GAM ranged from protein 0.30% to 54.94% for fat, while heritability ranged from flour moisture and dry matter (17.29%) to 84.88% for cyanide. The first five principal components explained 80.1% of the total variation, with PC I accounting for 37%, PC II 15.4%, PC III 11.6%, PC IV 8.4%, and PC V 8.20% of the total variation. This study found the presence of high biochemical variability among the tested accessions’ roots and could be used to select accessions with desirable biochemical composition in future breeding work.

1. Introduction

Cassava (Manihot esculenta Crantz) is a lucrative crop in terms of food calories generated per unit of area and time, much exceeding other staple crops [1]. Cassava produces more than 250,000 food calories per hectare per day, compared to 200,000 for maize, 176,000 for rice, and 110,000 for wheat [1]. Nonetheless, cassava, like some crops, has poisonous and antinutritional components. Of particular concern are the cyanogenic glucosides of cassava, which include linamarin and lotaustralin at concentrations ranging from 150 to 300 mg/kg in peeled root or 300 to 900 mg/kg in a dry matter [2, 3]. These are hydrolyzed by the plant’s endogenous enzyme, linamarase, to release free cyanide. Cyanide suppresses cellular respiration in all aerobic organisms by limiting mitochondrial electron transport and preventing oxygen intake [3].

The hazardous hydrogen cyanide concentration of cassava compromises its dietary value [4]. The presence of cyanide in cassava has been reported as a source of concern for human and animal consumption [5]. The amount of these antinutritional and potentially harmful glycosides varies greatly depending on variety, climate, and cultural factors [1, 6]. Cyanide content is used as a criterion for selecting cassava cultivars and their flour for use in the food industry and breeding programs [7]. Cassava cultivars are divided into two classes based on their cyanide content: high-HCN or bitter cassava with more than 100 ppm of cyanogenic equivalents and low-HCN or sweet cassava with less than 50 ppm [8]. Some authors categorize cassava cultivars as bitter when the cyanide is spread throughout the tuberous root at levels greater than 100 mg/kg fresh root weight. The other (sweet/cool) cultivars are those in which cyanide is contained primarily in their peels at a low level. The sweet cassava can be cooked and eaten raw, whereas the bitter cassava requires processing before consumption [8].

As a result of inherent characteristics that differ from one cassava genotype to the other, cultivars play a vital role in the creation of a wide range of food products. Proximate composition (proteins, lipids, fiber, ash, and moisture), starch yield, and dry matter content are examples of such parameters [9]. These authors indicated that proximate, dry matter, starch yield, and cyanide content are some of the most important quality indicators in the food sector when selecting raw materials. As a result, determining the starch and other chemical components of cassava cultivars might be helpful in selecting cassava cultivars for various food formulations, processing, and ultimately industrial applications [10]. Cultivar, geographical location, plant maturity stage, and environmental variables are elements that influence chemical composition in plants [2, 11]. The biochemical composition of cassava roots varies based on cultural methods such as root pruning, age at harvest, maturity at harvest, storage environment, area of growth, and postharvest practices [10].

Cassava’s nutritional value is determined by the plant portion (root or leaves), cultivar, age, geographical location, and environmental conditions [11, 12]. Cassava storage roots are a good source of energy, with carbohydrate content ranging from 32% to 35% on a fresh weight (FW) basis and from 80% to 91% on a dry matter (DM) basis [13]. Starch content varies with genotypes, in which improved, landraces, account for 73–85% of dry root weight [14]. The high starch content (18–24% amylose and 70% amylopectin) allows easy digestion [15]. Bitter cultivars have low levels of glucose, fructose, sucrose, and maltose [12], whereas sweet cultivars have a sucrose concentration of more than 17% [16]. The lipid content of cassava root varies between 0.1% and 0.3%of fresh weight [17, 18]. The protein content is low, ranging between 1% and 3%of dry matter and 1.5 mg/100 g of fresh mass [19]. Cassava cultivars with high moisture and reducing sugar contents could be employed as raw materials in the ethanol, organic acids, lactic bacteria, and biofuel industries [18]. Fresh cassava root, on the other hand, is very perishable due to its high moisture content (33–72%) [11, 20] and has a short postharvest life of fewer than 72 hours. Cassava is then processed into shelf-stable primary products such as flour, chips, and pellets shortly after harvest.

Many smallholder farmers in Ethiopia grow a variety of root and tuber crops due to the country’s diversified agroecologies and favorable conditions [21]. In Ethiopia, Cassava has great adaptation and growth performance in a variety of agroecologies with varying productivity [22]. It is an important food crop in Southern Ethiopia, providing a significant amount of the family’s daily meals and serving as a key source of carbohydrates. Currently, farmers are growing cassava as a food security crop and a significant source of household income. Despite the fact that cassava has a wide range of culinary applications, the majority of its products are consumed in Ethiopia by boiling the root and flour [23]. There has also been limited research information on the chemical composition of cassava accessions and knowledge in biochemical compositions. Therefore, the present study aimed to assess the biochemical composition n of cassava accessions collected from various areas and provide information for utilization and conservation.

2. Materials and Methods

2.1. Description of the Experimental Site

The experiment was conducted in the postharvest and animal nutrition laboratory of Jimma University College of Agriculture and Veterinary Medicine (JUCAVM) in the year 2021.

2.2. Experimental Material

The study used 64 cassava accessions among which 15 were provided by the International Institute of Tropical Agriculture (IITA), Nigeria. The remaining 49 accessions were obtained from Jimma Agricultural Research Center (5) and Hawassa Agricultural Research Center (44) accessions (Table 1).

2.3. Samples Collection and Preparation

Root samples of 0.5 kg fresh weight each sixty-four cassava accessions were collected from the cassava field experiment located at Tarcha. Cassava roots were gathered when they were 18 months old and brought to the laboratory within three hours. To remove sand and other dirt particles, the samples were cleaned, peeled, and washed with room temperature water. A sample was then processed using three typical processing procedures (chipping, grating, and soaking). The peeled root was soaked for 48 hours in plastic drums filled with water, after which the cassava roots were removed, cut into pieces, and oven dried for 48 hours at 60°C. The dried chips were subsequently crushed into fine-powder cassava flour using an electric grinder [24, 25]. The flour was sieved through a 1 mm sieve, measured, placed into an airtight plastic bag, and kept in the refrigerator until the analysis was done. Triplicate analyses were performed on each sample of flour.

2.4. Biochemical Analysis

The biochemical analysis was performed with cassava flour, and triplicate samples were obtained for each accession in the postharvest and animal nutrition laboratory of Jimma University College of Agriculture and Veterinary Medicine (JUCAVM). Standard procedures [2628] were used to determine the moisture content, dry matter content, ash, organic matter, crude fiber, crude fat, starch, tannin, crude protein, saponin, carbohydrate, cyanide, and total energy contents of the flour. The following are the details of the analyses.

2.4.1. Determination of Moisture Content

The moisture content of flour was determined using a typical standard analytical method [28], and a sample of 100 g duplicate flour was weighed in aluminum dishes and oven dried for three days at 65°C. The dried sample was then weighed after cooling at room temperature in a desiccator. The moisture content of flour was estimated by converting the weight loss due to drying to percent flour moisture content as follows: percentage flour moisture = (weight of moisture evaporated/weight of flour sample) × 100.

2.4.2. Determination of Dry Matter Content

The flour dry matter content (DM) was estimated by collecting a representative duplicate sample of 100 g (W1), oven drying it at 65°C for 72 hours, weighing it (W2), and expressing the value in percentage [29, 30]. The percentage of dry matter content was calculated as follows: % DM = (W2/W1) × 100 or % DM = 100 − % moisture content.

2.4.3. Determination Total of Ash Content

The guidance of [27] was followed to determine the ash content. The crucibles were rinsed and dried for 30 minutes at 105°C in a hot air oven (SM9053). These were weighed after cooling in desiccators. To get rid of the smoke, a five-gram sample was burned inside a heater in a fume cupboard. The sample was then placed in a preheated muffle furnace (SM9080) and heated to 550°C until a light grey ash appeared, after which the crucibles were cooled in desiccators and weighed. The ash content was determined as follows:

2.4.4. Determination of Organic Matter Content

By removing the percent ash from the % total dry matter, the organic matter content was calculated and reported as a percentage. The organic matter content was determined as follows: organic matter content (%) = % DM − % ash.

2.4.5. Determination of Crude Fiber Content

Nonenzymatic gravimetric technique was used to determine the crude fiber content [27,920.168]. After defeating the two-gram sample with petroleum ether, 200 ml of a solution containing 1.25 g of H2SO4 per 100 ml was heated in reflux for 30 minutes [31]. Finally, using a fluted funnel and linen, the solution was filtered. After 5 washings with hot water and trapping the particles with a two-food muslin cloth, the sample was quantitatively transferred back to the flask and boiled for 30 minutes in 200 ml of 1.25 g of carbonate-free NaOH per 100 ml, then washed before being transferred to a weighed Gooch crucible and dried in the oven at 105°C for three hours. After cooling in desiccators, it was reweighed. The percentage of crude fiber was estimated as follows:

2.4.6. Determination of Crude Fat

According to [27], the crude fat was calculated using a fat extractor with an automated control unit (FOSS Soxtec 2055). A one-gram sample was weighed into the thimble, which was then put into the extraction unit with its mouth plugged with defatted cotton wool. Eighty milliliters of petroleum ether were poured into each cup, which was then aligned with their respective thimbles at 135°C. For 30 minutes, each extraction and rinsing was performed, then the sample was aerated for 15 minutes, and the crude fat % was determined as follows:where  = weight of sample,  = weight of empty cup, and  = weight of cup with the extracted oil.

2.4.7. Determination of Crude Protein

The technique 988.05 micro-Kjeldahl method of nitrogen analysis was used to determine the crude protein content of the samples [28]. About 0.3g of sample was measured using an analytical balance (Model: ABJ220-4M, Australia); 1g of K2SO4 and CuSO4 catalyst mixture, and 5 ml of sulfuric acid were added to each digestion flask (Kjeldahl flask KF250, Germany), which contained the mixture of sample and catalysts. The solution (0.3 gram sample + 1 g K2SO4 and copper sulfate + 5 ml H2SO4) was immediately placed in a digestion flask and heated to around 420°C for 3–4 hours, until it became clear.

The digested sample was then passed to the distillation equipment, where 25 ml of 40% (w/v) NaOH was continuously added until the solution became murky, indicating that the solution had become alkaline. The ammonia was collected into a 200 ml conical flask containing 25 ml of 4% boric acid with mixed methyl red indicator solution after the solutions were steam distilled. Then, with the delivery tube below the acid level, distillation was carried out into the boric acid solution in the receiver flask. The pink color solution in the receiver flask turned green as the distillation progressed, showing that ammonia was present. Distillation was maintained until the flask’s contents reached the appropriate level. Titrations of the green color solution against 0.1 N HCl solutions were then performed. The green color changed to a reddish pink tone near the conclusion, indicating that all of the nitrogen bound as ammonium borate had been released as ammonium chloride. To get a reddish color, the distillate was titrated with 0.1 N sulfuric acid. Finally, using the formula below, the percentage of nitrogen content was calculated.where VA = volume (ml) of the HCl solution consumed in the sample titration, VB = volume (ml) of the standard solution used in the sample blank titration, N = Normality of hydrochloric used which was 0.1 N, and W = weight of the sample (g).

The crude protein content was estimated using the formula equation (5).

2.4.8. Starch Determination

The wet approach was used to extract starch [32]. Roots were washed, peeled, and diced into 1 cm cubes before being pulverized for 5 minutes in a high-powered blender (Model KING, Osaka, Japan). The pulp was suspended in ten times its volume of water, stirred for five minutes, and then filtered through a double-fold cotton towel. The top liquid was decanted and discarded after the filtrate had been allowed to settle for 2 hours. The sediment was restirred for 5 minutes after the addition of water. The starch from the filtrate was allowed to settle after another round of filtration. The sediment (starch) was sun dried for 24 hours after decanting the top liquid, then kept and converted to a percentage by mass. %starch = (weight of starch isolated/weight of sample root) × 100.

2.4.9. Determination of Total Cyanide Contents

The AOAC method [26] was used to quantify total cyanide concentration. Cassava flour (20 g) was placed in a 1-liter distillation flask with 100 ml water and let to stand for 2 hours to hydrolyze the bound glucosides and liberate hydrocyanic acid. Following that, 100 mL of distilled water was added to the slurry in the distillation flask (autolysis), and the steam distillate was collected in 20 ml of 0.01 N AgNO3 that had been acidified with 1 ml of HNO3. The distillation was carried out for another 40 minutes with rapid boiling until approximately 150 ml of the distillate was collected twice. The distillates were pooled after being filtered through Gooch with some water. The excess of AgNO3 was then measured using ferric alum indicator in the combined filtrate and washings with 0.02 N KSCN. After the addition of 0.02 N KSCN solutions, the titration reached its end point when a faint crimson tint appeared. The following equation was used to compute the amount of HCN present in the sample.Volume (ml) of AgNO3 consumed to complex CN- = 20 − 2V of the titer; V is the volume of the titrant consumed 1 ml 0.01 NAgNO3 = 0.27 mg HCN.

2.4.10. Determination of Condensed Tannin

The technique adopted by Udosen [33] was used to determine the tannin concentration. Each sample was weighed and placed in a centrifuge tube with 2 ml of distilled water. It was centrifuged at 1500 rpm for 10 minutes. The supernatant (extract) was then dispersed after the centrifuge samples were poured into a beaker. In the beaker, 1 ml of NaCO3 and Folin–Denis reagent was added and allowed to settle. A spectrophotometer was used to take the measurements. The following formula was used to determine the amount of tannin in the sample.where An = absorbance test sample, As = absorbance of standard sample, C = concentration of the standard solution, W = weight of the sample, Vf = total filtrate volume, and Va = volume of filtrate analyzed.

2.4.11. Determination of Condensed Saponin

The method described by Obadoni and Ochuko [34] for determining saponin content was used. In 200 ml of 20% ethanol, 20 grams of each ground sample were distributed. At around 55°C, the suspension was cooked for 4 hours in a hot water bath with constant stirring. The residue was re-extracted with 200 ml of 20% ethanol after the mixture was filtered. Over a water bath at around 90°C, the mixed extracts were reduced to 40 ml. The concentrate was poured into a 250 ml separator funnel, along with 20 ml of diethyl ether, and rapidly agitated. The aqueous layer was kept, while the ether layer was thrown away. It was necessary to repeat the purifying procedure, and n-butanol was added in the amount of 60 ml. The n-butanol and extract mixture was rinsed twice with 10 ml of 5% sodium chloride aqueous solution. In a water bath, the residual solution was heated to around 90°C. The samples were dried in an oven at 100°C until they reached a consistent weight. The saponin content was calculated in percentage and changed to mg/100 g.

2.4.12. Determination of Carbohydrate and Energy

Using arithmetic difference, the carbohydrate content of the sample was estimated [35, 36]. With 4 calories for 1 g of carbohydrates, 4 calories for 1 g of proteins, and 9 calories for 1 g of crude fat, the energy value was estimated using Atwater and Rosa’s thermal coefficients (1899). Using the formula below, the available carbohydrate (CHO) and energy value were calculated: CHO = [100 − (% moisture + % crude protein + % crude fat + crude fiber + % ash)] [35, 36].where CHO = carbohydrate, CP = crude protein, and CF = crude fat.

2.5. Data Analysis

Data were subjected to one-way analysis of variance using the complete randomized design (CRD) by using the Statistical Analysis System (SAS) package [37] and R software package [38] to determine the presence of variation among the accessions for various biochemical content. The data were standardized by dividing each variable by its range, then clustered using the unweighted pair group method of arithmetic mean (UPGMA) and pairwise generalized square distances (D2) between clusters, with principal component analysis (PCA) used to group the accessions based on biochemical character and assess correlations between principal components and the parameters measured.

2.5.1. Variance and Covariance Analysis

Phenotypic and genotypic variances and covariances were calculated according to the methods suggested by [39, 40].(a)Genotypic variance component(1)Genotypic variance () = (MSg − MSe)/r, where Genotypic mean square = MSg, error mean square = MSe, and replication = r(2)Genotypic coefficient of variation (GCV) = , where is the grand mean value of the trait(b)Environmental variance component (on genotypic mean basis) σ2e = MSe/r(c)Phenotypic variance component(1)Phenotypic variance (σ2p) =  + σ2e(2)Phenotypic coefficient of variation (PCV) = (d)CovarianceCovariance was calculated using the following formula:(1)Genotypic covariance (Cov A) = rAσXAσYA, where rA = correlation of breeding values arises from two sources, σXA = standard deviation of trait x, and σyA = standard deviation of trait y(2)Environmental covariance (CovE) = rEσXEσyE, where rE = correlation of environmental deviations, σXE = environmental standard deviation of trait x, and σyE = environmental standard deviation of trait y(3)Phenotypic covariance (CovP) = Genotypic covariance (CovA) + environmental covariance (CovE)

The genetic advance with one cycle of selection, supposing the selection intensity of 5% was estimated as proposed by [41].where σp = the standard deviation of the character phenotypic and k = the standardized selection differential at 5% selection intensity (2.063).

The extent of the estimated genetic advance of different traits under selection, GAM, has been calculated by the following formula:where genetic advance as a percent of the mean = GAM, genetic advance = GA, and  = grand mean for the character.

2.5.2. Correlation Analysis

Correlation coefficient (r) was calculated using the standard procedure suggested by [42].

The genotypic correlation coefficient between traits x and y was calculated as follows:where  = genotypic correlation coefficient between traits x and y,  = genotypic covariance between trait x and y,  = genotypic variance of trait x, and  = genotypic variance of trait y.

Genotypic correlation coefficients were tested for their significance using the formula suggested by Robertson [43], using the t-table at () degrees at 5% and 1% levels of alpha.

3. Results and Discussion

3.1. Analysis of Variance

The mean squares due to accessions were very highly significant () for ash, fiber, fat, protein, cyanide, tannin, and saponin content, whereas organic matter, starch, carbohydrate, and total energy were observed to have highly significant () difference among the tested accessions. Similarly, flour moisture and dry matter content showed significant () differences among accessions (Table 2). This implies the presence of high genetic variability of these characters among the tested accessions, and the observed variability is a good possibility for developing cassava cultivars having desirable biochemical characters.

3.2. Mean and Range Value of Biochemical Characters

The descriptive statistics of the cassava accessions root based on biochemical compositions are presented in Table 3, and wide ranges of variation were identified for most of the biochemical characters studied. The moisture percentage of cassava root flour varied greatly among the 64 accessions, ranging from 4.83% to 10.11%, with a mean value of 7.11 ± 0.16. Accessions G6, G13, G37, and G64 had high moisture levels, whereas accessions G8, G18, G28, and G42 had low moisture contents, among others (Tables 1 and 3). The variation in the moisture content of cassava root flour might own to differences in genetic makeup for solute constituents. According to this study, G8, G18, G28, and G42 were more suitable for long-term storage of their roots than those with a high moisture value. This finding agrees with those reported by Adejumo [44] and Baah et al. [45], in which the moisture content of cassava root flour ranged from 6.68% to 10.96%. The values of dry matter ranged from 89.89% to 95.17% with a mean of 92.89 ± 0.16%. Similarly, organic matter varied from 86.71% to 92.65% with a mean of 89.82 ± 0.16% (Table 3). Among the tested accessions, G8 and G18 had higher dry matter and organic matter contents, while G13 and G27 had lower values (Tables 1 and 3). This implies that accessions G8 and G18 have good eating quality and are more storable than others because of their high dry matter content. With a mean value of 3.07 ± 0.05, the ash content ranged from 2.10% to 3.96%. Among the accessions studied, the G51 had the lowest ash level (2.10%), and the G30 had the highest (3.96%). These ash values were greater than those reported by Rojas et al. [46], which ranged from 1.5% to 2.7%, but were similar to those reported by Offor et al. [47], which ranged from 2.3% to 3.6%, with a mean value of 2.65%. This indicates that high ash content among the examined accessions results in high mineral content, and in contrast, low ash content indicates low mineral concentration in cassava [4].

The range of crude fat was 0.26% to 1.40%, with a mean value of 0.83 ± 0.03%. Accessions G34 and G35 had the highest crude fat values, while G47 had the lowest (Tables 1 and 3). Findings of the current study differ from those of Sarkiyayi and Agar [7], who reported high-fat values (3.92% for sweet cassava varieties and 3.82% for bitter cassava varieties), while Manano et al. [10] found low crude fat values both for improved (0.48% to 0.63%) and local (0.39% to 0.48%) varieties. The difference in crude fat between accessions may be due to genetic differences rather than environmental factors. The fiber content of the various cassava accessions roots ranged from 1.14% to 3% on a dry weight basis, with a mean value of 2.19 ± 0.05 (Table 3). Roots from the accessions G38 and G41 exhibited the highest fiber value (Table 1), while those from the accessions G54 and G 52 were identified as having the lowest value fiber content. Similar trends are shown by Manano et al. [10], who reported cassava varieties with fiber contents of 1.06% to 1.2%. The results, therefore, suggest that the accessions G38 and G41 are good sources of dietary fiber. This could be due to the high crude fiber content in these accessions. The starch content of the cassava root on a dry weight basis varied from 65.1% to 74.2% with a mean value of 72.11 ± 0.24% across the examined accessions. The highest starch content was found in accessions G11 and G46, while the lowest was found in accession G6 (Tables 1 and 3). This finding is consistent with that of Chijindu and Boateng [48], who reported 76.82% starch content for cassava varieties typically grown in Ghana, and Chiwona-Karltun et al. [49], who reported starch content of cassava varieties to vary from 54.7% to 76. 34%. The results of the current study indicate that the starch content is greater than 60%, and these accessions could potentially be used for various commercial products such as starch, alcohol, and glucose.

The crude protein and carbohydrate percentages ranged from 1.28% to 2.46% and 81.29% to 87.84% with a mean of 1.95 ± 0.03% and 84.84 ± 0.17%, respectively (Table 3). Accessions G10, G28, and G42 had high crude protein values, whereas G27, G33, and G50 had low crude protein value. These findings are similar to those of Rojas et al. [46], who reported crude protein levels ranging from 1.5% to 2.8%, with an average of 2%, and Nyirendah et al. [20], who observed values ranging from 1.2% to 3.5% in six cassava varieties. However, Manano et al. [10] found a lower crude protein content, ranging from 0.74% to 1.5%, and carbohydrate content ranging from 83.86% to 91.33%. The study suggests that crude protein content variation in cassava accessions may be due to genetic differences rather than environmental factors.

The amount of cyanide in each sample ranged from 1.67 mg/100 g to 3.14 mg/100 g with a mean performance of 2.47 ± 0.04 mg/100 g, and the amount of tannin in each sample ranged from 0.17 mg/100 g to 0.36 mg/100 g, with a mean value of 0.25 ± 0.01 mg/100 g. Saponin concentrations ranged from 1.74 mg/100 g to 2.84 mg/100 g, and the average was the performance of 2.29 ± 0.04 mg/100 g (Table 3). In comparison to other accessions, G38, G39, G40, and G63 had high cyanide concentrations, whereas G9, 29, and G51 had low cyanide content in the dry weight base (Tables 1 and 3), suggesting that accessions, G9, G29, and G51 might be suitable for the food industry. In line with this work, Charles et al. [50] reported 0.8 to 2.9 mg/100 g cyanide content on a dry weight basis. Nyakaisiki [51] reported cyanide contents ranging between 2.8 and 5.3 mg/100 g on a fresh weight basis. But Manano et al. [10] reported elevated values, ranging from 3 mg/100 g to 80 mg/100 g cyanide contents on a dry weight basis.

3.3. Variances, Heritability, Coefficient of Variation, and Genetic Advance

All characters evaluated were given their corresponding variance components and coefficients of variation values (Table 4). Total energy (11.31), followed by carbohydrate (0.78) and organic matter (0.65), had the highest genotypic variation (), among the evaluated characters, while tannin (0.002) had the lowest genotypic variation (). Similarly, phenotypic variances (σ2p) ranged from 0.002 for tannin to 60.59 for total energy. Crude fat had the largest genotypic coefficient of variation (GCV = 29.54%), followed by crude fiber content (GCV = 16.94%) and tannin (GCV = 16.11%), while dry matter had the lowest (GCV = 0.84%). The phenotypic coefficients of variation (PCV) for dry matter content ranged from 2.02% to 32.74% for crude fat content (Table 4). According to Deshmukh et al. [52], PCV and GCV values greater than 20% are considered high, whereas values between 10% and 20% are considered medium, and values less than 10% are considered low. Based on this classification, crude fat content was the only character that met the high range of GCV (29.54%) and PCV (32.74%) values; tannin, cyanide, ash, saponin, flour moisture, and crude protein were categorized under the medium range of GCV and PCV values. The current study revealed that the biochemical composition of the tested accessions was highly variable. As a result, it is suggested that selecting a character of interest could be effective if high and medium PCV and GCV values are taken into account.

The expected genetic progress and estimated heritability in the broad sense for the investigated characters are reported in Table 4. Heritability in the broad sense in the current study ranged from 17.29% for flour moisture and dry matter percent to 84.88% for cyanide, while genetic advance as a percentage means ranged from 0.30% for crude protein to 54.94% for crude fat. According to Bhateria et al. [53], cyanide (84.88%), crude fat (81.32%), crude fiber (77.32%), ash (72.13%), tannin (66.67%), saponin (60.96%), and crude protein (54.94%) have high heritability, but the rest of the characters fall under low heritability (Table 4). This suggests that effective genetic progress can be made if some of these characters with high heritability broad sense estimates are viewed as selection criteria for the character of interest in the future cassava improvement program. Crude fat had the largest genetic advance (54.94%), followed by crude fiber (30.68%), tannin (27.11%), cyanide (26.73%), and ash (21.21%) (Table 4). High heritability along with high genetic advance as a percent of the mean was estimated for crude fat, crude fiber, tannin, cyanide, and ash (Table 4). It is hypothesized that these characters with high broad-sense heritability and high genetic progress as a percentage of the mean can be improved by using phenotypic expression to direct the selection of the character of interest.

3.4. Correlation among Character

The coefficients of person correlation for thirteen characters are shown in Table 5. The correlation results revealed that flour moisture content had highly significant () strong negative correlation with dry matter (r = −1), organic matter (r = −0.95), starch content (r = −0.89), carbohydrate content (r = −0.85), and total energy content (r = −0.85) but a highly significant and medium negative correlation with crude protein (r = −0.42). Similarly, dry matter demonstrated a strong positive and very significant relationship with organic matter (r = 0.95), starch content (r = 0.89), carbohydrate content (r = 0.85), and total energy content (r = 0.91), and it was observed as a highly significant and medium-positive correlation with crude protein (r = 0.42) (Table 5). These results agree with the study by Adebola et al. [54], who reported that moisture content was strongly and negatively correlated with dry matter content (r = −0.91), crude protein (r = −0.64), and carbohydrate (r = −0.67). These authors found that dry matter contents correlated with crude protein (r = −0.52) and carbohydrate (r = −0.63). In the other study, Mégnanou et al. [55] discovered that moisture content was strongly correlated with dry matter (r = −1), carbohydrate (r = −0.985), starch (r = −0.733), and energy value (r = −0.971), while dry matter was correlated with carbohydrate (r = 0.985), starch (r = 0.733), and energy value (r = 0.971).

Organic matter content was observed to have a highly significant and strong positive correlation with starch content (r = 0.91), carbohydrate content (r = 0.91), and total energy content (r = 0.95); organic matter and ash were significant () and inversely linked (r = −0.3). Likewise, starch content showed a highly significant and strong positive correlation with carbohydrate content (r = 0.99) and total energy content (r = 0.89). Crude protein was moderate, positive, and highly significantly correlated with starch (r = 46) and total energy content (r = 0.42). Tannin and saponin had a strong link with ash (r = 0.26) from the antinational component, while saponin had a substantial correlation with crude fat (r = 0.43) (Table 5).

Many characters have been found to have strong relationships in the current investigation. However, a few attributes are correlated with antinutritional factors, and the majority of these correlations are nonsignificant (Table 5). The biochemical characters show a stronger and more positive relationship with the parameters tested, implying that any rise in one biochemical character will raise many characters. Furthermore, the observed relationships between the biochemical characters suggest that improving one character will lead to improvements in the others.

3.5. Principal Component Analysis

The presence of genetic variation in the accessions was revealed by principal component analysis (PCA), with the first five principal components with eigenvalues >1 ranging from 1.004 to 4.808, contributing to a total variation as elevated as 80.01% (Table 6). The first principal component explained 37% of the total variation and was positively associated with organic matter (0.450), dry matter (0.440), total energy (0.439), and carbohydrate (0.423), while flour moisture content (−0.44) contributed negatively (Table 6). Principal component II was correlated with crude fat (0.513), saponin (0.461), ash (0.444), and fiber content (0.376) and explained 15.4% of the total variation, whereas principal component III was correlated with starch content (0.649), tannin (0.407), cyanide (−0.401), and crude protein (0.343) and explained 11.6% of the total variation (Table 6). With increased crude fat, saponin, ash, and fiber content, the principal component of the second increases. Therefore, it implies that with an increase in crude fat, saponin, ash, and fiber contents tend to increase too. The principal component of the second may be regarded as a measure of the content of crude fat, saponin, ash, and fiber. Similarly, the third principal component increases simultaneously with increasing starch content, tannin, and crude protein but decreases cyanide content. This indicates that cassava genotypes containing high starch, tannin, and crude protein content tend to have decreased cyanide content.

The fourth principal component (PC IV) accounted for 8.4% of total variability and was mainly associated with crude protein (0.551), fiber (0.414), ash (−0.430), and tannin (−0.464), whereas the 7.7% total variation for the fifth principal component (PC V) was positively correlated with the fiber (0.542), cyanide (0.461), and tannin (0.346) but was negatively correlated with saponin (−0.477) (Table 6). The VI, VII, and VIII principal components accounted for 6.5%, 5.4%, and 4.4% of the total variation, correspondingly, and these principal components were strongly correlated with cyanide (0.596) and crude protein (0.546) for PC VI, ash (−0.608) and tannin (0.511) for VII, and starch (0.599) for PVIII (Table 6).

3.6. Cluster Analysis

It is essential to group accessions according to their similarity or difference. In the present study, a total of 64 cassava accessions were clustered into six diverse groups based on thirteen biochemical characters (Figure 1 and Table 7). Among the clusters, cluster I was the largest, comprising 30 (46.87%) accessions, of which 6 accessions were collected from the International Institute for Tropical Agriculture, 2 accessions from Jimma Agricultural Research Center, and 22 accessions from Hawassa Agricultural Research Center (Figure 1 and Table 7). The clusters II, III, and IV were represented by 20, 6, and 6 accessions within diverse origin sources, respectively, and the whole clusters contributed about 50.01% of the total variations. Finally, clusters V and VI were represented by 1 accession separately and contributed 3.2% of the total variation (Figure 1 and Table 7).

The cluster means for different biochemical characters revealed that wide variations were observed among the cluster means for biochemical characters (Table 8). The highest flour moisture content was shown in cluster mean IV (9.58%), followed by cluster mean II (7.89%) and cluster mean I (6.55%), while cluster mean III (5.41%) was the lowest flour moisture content among the clusters (Table 8). This indicates that cluster means III accessions are better suited for long-term root storage than cluster means IV and II, which have higher moisture values. In their clusters, accessions from cluster means I and III produced the highest levels of dry matter (93.45% and 94.59%) and organic matter (90.38% and 91.71%), respectively (Table 8). Likewise, the maximum cluster mean was shown for ash and fiber in cluster II; for crude protein, starch, carbohydrate, and total energy in cluster III; for cyanide in IV; for tannin in II; and for saponin in III. In the cluster means, the lowest antinutritional was observed for cyanide in cluster III and for tannin and saponin in IV (Table 8). This suggests that genotypes from cluster mean III can be employed for the development of variety in terms of high nutritional compositions.

4. Conclusion

Among the tested accessions, G8 and G18 are the best accessions for long-term root storage and eating quality because of their high dry matter and organic matter levels. The highest starch concentration is recorded in G11 and G46, showing that these accessions might be used for making a variety of commercial products such as starch, alcohol, and glucose. Accessions G9, G29, and G51 had lower cyanide values in the dry weight base than other accessions, indicating that they are appropriate for the food industry.

Along with the medium to high heritability, GCV and genetic advance as a percent of the mean value are shown for ash, crude fiber, crude, cyanide, tannin, and saponin. Thus, these characters can be used as good criteria for the future cassava nutritional composition improvement program. Accessions from cluster mean III have high dry matter, organic matter, crude protein, starch, carbohydrate, total energy, and low antinutritional factors. This is an indication that accessions belonging to cluster mean III could be used as an initial population for further study in biochemical composition improvement. In conclusion, the present study showed the presence of high biochemical variability among the tested accessions and could be used to select accessions with a desirable biochemical composition in future breeding work.

Data Availability

The data sets used to support the findings of this study are available at the hands of the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

BB carried out the fieldwork and laboratory work, performed the statistical analysis, and drafted the manuscript. The rest of the authors coordinated the study, supervised fieldwork, and laboratory work, and contributed to the write-up of the manuscript. All authors have read and approved the final manuscript for submission.

Acknowledgments

The authors would like to thank Jimma and Wolita Sodo Universities for financing the research work. The authors are also pleased to acknowledge Jimma Hawassa Agricultural Research Centers and International Institute for Tropical Agriculture for providing planting materials used for the study. Jimma University and Wolita Sodo University financed all the research work.