INTRODUCTION
Cassava (Manihot esculenta Crantz) is of great importance worldwide. It is a staple food for more than one billion people in more than 100 countries around the world (Albuquerque et al., 2012). Its cultivation is widely developed in tropical Africa, Asia and Latin America, mainly in underdeveloped and developing countries. In regions such as Northeast Brazil, countries such as Ghana and Nigeria and some parts of Indonesia, cassava provides about 70% of the daily calorie consumption of the population (Nassar et al., 2006).
Its vast cultivation is mainly due to its genetic characteristics, which include rusticity, tolerance to drought and acidic soils, high adaptive capacity and low production cost. It has become one of the main sources of carbohydrates in the world (Bester et al., 2021). According to the latest survey by the Food and Agriculture Organization of the United Nations, global production of cassava roots was 291 million tons in an area of approximately 26.3 million hectares, with an increase of 55 and 66%, respectively, since 2000 (FAO, 2019).
Through genetic improvement, cassava crop has made great progress, such as an expressive increase in productive potential and nutritional root quality through new cultivars. Bester et al. (2021), in a study on a 30-year historical series, reported an increase in cassava productivity in Brazil. Species diversity is a resource that can be used as a source of genetic variation for incorporating genes into new varieties (Pádua, 2018). Many genes of great importance have not yet been found in this vast germplasm bank, whether they are linked to increased productivity, starch content and quality or tolerance to drought, pests and diseases, representing an important line of research (Oliveira et al., 2016).
Selection strategies were created, including the Restricted Maximum Likelihood (REML) and Best Linear Unbiased Prediction (BLUP) methods. These methods are used to estimate parameters and predict genetic values, without the environmental effects on phenotype, in which only the studied genetic characteristics remain. They are important methods in guiding breeding programs (Pimentel et al., 2014). Thus, the mixed model makes it possible to obtain variance components and genetic parameters through the restricted maximum likelihood (Baretta et al., 2016) and can help in the selection and prediction of genotypes through the best non-targeted linear predictor (BLUP), indicating the genetic value and the new predicted mean (Carvalho et al., 2016).
For the importance of nutritional quality and the need to demonstrate the effects of the cultivation environment on these characteristics, another method of analysis called Annicchiarico can be used, which demonstrates the stability of genotypes in a given environment. The results reflected a confidence index for each genotype that showed which genotypes were superior and classified environments as favorable or unfavorable for the relevant characteristic (Cruz et al., 2014).
This study aimed to highlight the behavior of cassava cultivars when subjected to different densities and biostimulants at planting and to select superior cultivars based on nutritional and productive attributes using the multivariate approach.
MATERIALS AND METHODS
This study was carried out in the municipality of Augusto Pestana/RS (Brazil). The soil was classified as Typical Distroferric Red Latosol (Santos et al., 2018). According to the Köppen climate characterization, the climate of the region was Cfa type.
The experiment design used randomized blocks with three replications in a 3×2×2 triple factorial scheme (three cultivars, two planting densities and with or without plant regulator). The treatments consisted of three cassava cultivars: FEPAGRO-RS 13 Vassourinha, BRS CS01, and Iapar - 19 Pioneira in four environments, namely: density 10 (10 buds per linear meter) with biostimulator, density 10 without biostimulator, density 20 with biostimulator (20 buds per linear meter) and density 20 without biostimulator. The experiment consisted of 36 plots of 6 linear meters, according to the treatments, with a spacing of one meter between lines.
The soil was prepared with subsoiling and opening of the planting lines. At the same time, the fertilization and soil correction were carried out according to the soil analysis for the area following the technical recommendations for cassava crops (CQFSRS, 2004). The crop was planted at the end of September 2019 by manually transplanting stems that were deposited in a horizontal position at the bottom of the furrows.
The cassava harvest was carried out 9 months after the transplanting the cuttings. The start was done manually where the roots were collected, cleaned and sent to the UNIJUÍ plant production laboratory, where productivity measurements were taken (t ha-1), and drying was done in an oven (105±5°C). After drying, the samples were ground in a willey mill and sent to the Food Microbiology Laboratory of the Federal University of Pelotas, where the following variables were measured:
Mineral material (MM, %): the Mineral material content was obtained after incineration of the sample in a muffle furnace at a temperature of 550°C, following the method described by the Association of Official Analytical Chemises (AOAC, 2005); lipid content (oil, %): the lipid content was determined with solvent extraction using a Soxhlet extractor with ether as the reagent, following the method described by Soxhlet (1879); protein content (P, %): the protein determination was based on nitrogen with the Kjeldah digestion process. The organic matter was broken down, and the existing nitrogen was finally transformed into ammonia, introducing an empirical factor of 6.25 to transform the number of grams of nitrogen into the number of grams of protein, according to the AOAC procedures; starch content (S, %): The tuber samples were immersed in water overnight, then ground in a blender, and the starch was filtered through a fine mesh (opening size: 350 μm), immersed in water to separate the husks from the suspension, and allowed to settle overnight. After decanting, the sedimented starch was rinsed several times with distilled water, pressed onto a clean muslin cloth and dried at 45°C in a hot air oven (D-37520, Thermo Fischer Scientific, Pretoria, South Africa) for 12 h. The dry starch was stored at 4°C until analysis, a maximum period of one week; this methodology was adapted from the article by Oyeyinka et al. (2019).
The data were analyzed using the R software and submitted to analysis of variance (ANOVA) to detect the presence or absence of interaction between the factors. Then, based on this information, the means comparison test was performed for the variables with the Tukey method (P<0.05). Afterwards, the method of Annicchiarico (1992) was used, which was applied to the starch, protein, lipid, mineral material and productivity, following the methodology proposed by Cruz et al. (2014). To estimate the variance components and genetic parameters (REML) for half-sibling progenies and meet the assumptions of the experiment, the genotypic variance (σ²G), residual variance (σ²E), individual phenotypic variance (σ²F), broad-sense heritability for total genotypic effects (Ĥ²g), mean genotype heritability (Ĥ²mg), accuracy for genotype selection (řgǧ), genotypic coefficient of variation (CVg %), residual coefficient of variation (CVe %) and rate of the coefficient of variation (CVr) were estimated according to the methodology proposed by Ramalho et al. (2012). Subsequently, the Deviance analysis (LRT) was performed at P<0.05 probability with the chi-square test (X²) to identify the significance of the characteristic. BLUP (Best Linear Unbiased Predictor) estimates were used to obtain the components of the means. Finally, the multi-trait genotype-ideotype distance index (MGIDI) was used to identify genotypes that combined high productive performance with high nutritional quality.
RESULTS AND DISCUSSION
In the analysis of variance (Tab. 1), the cultivar × density interaction exhibited a significant effect only for starch, whereas the protein variable exhibited an effect with the cultivar × density × biostimulator interaction, and the density × biostimulator exhibited a significant effect on protein. The source of the cultivar variation showed a significant effect on the variables starch, protein and lipid. This result is similar to that found by Teixeira (2017), who, when evaluating nineteen varieties of table cassava, observed a difference between all characteristics, including starch, protein and lipids, with the exception of mineral material. On the other hand, in terms of density, only a significant effect was observed for starch, and the source of the biostimulator variation had an effect on the lipid content.
Variation factor | DF | Starch (%) | Protein (%) | Mineral material (%) | Productivity (t ha-1) | Lipid (%) |
---|---|---|---|---|---|---|
Block | 2 | 3.1 | 0.1 | 1.00·10-05 | 2862.8 | 9.00·10-05 |
Cultivar | 2 | 14.5* | 0.8* | 1.00·10-05 | 9755.2 | 0.20* |
Density | 1 | 11.6* | 0.007 | 0 | 5536.4 | 0.04 |
Biostimulator | 1 | 2.4 | 0.2 | 0 | 130.2 | 0.08* |
Cultivar × Density | 2 | 10.5* | 0.4 | 0 | 2185.8 | 0.06 |
Cultivar × Bioestimulator | 2 | 7.6 | 3.4 | 1.00·10-05 | 7695.4 | 0.01 |
Density × Bioestimulator | 1 | 7.9 | 0.2* | 0 | 1159.4 | 0.10 |
Cultivar × Density × Bioestimulador | 2 | 2.8 | 3.7* | 0 | 643.9 | 0.02 |
Residue | 22 | 1.4 | 3.7 | 0 | 3212.7 | 0.01 |
Total | 33 |
* Statistically significant at P<0.05.
The descriptive analysis (Fig. 1) showed that cultivar BRS CS01 had the highest productivity, and cultivar Iapar 19 - Pioneira exhibited the greatest variation between treatments for lipid content in the roots, cultivar FEPAGRO - RS 13 Vassourinha stood out with little variation between treatments, with Iapar 19 - Pioneira showing the lowest content, followed by a large variation between treatments. For starch content, there was no significant difference between cultivars and treatments. In terms of protein content, cultivar Iapar 19 - Pioneira stood out with the highest averages, and FEPAGRO - RS 13 Vassourinha presented the lowest content. For the mineral material content, cultivar BRS CS01 had the highest productivity; however, it did not differ from cultivar Iapar 19 - Pioneira. FEPAGRO - RS 13 Vassourinha had the lowest average.
For starch content (Tab. 2), density 10 had no statistical difference between the cultivars; however, at density 20, cultivar FEPAGRO - RS 13 Vassourinha showed the highest starch content, followed by BRS CS01 and Iapar 19 - Pioneira. When the spacing between plants was reduced, there was an influence on starch content. The shape of the root system may have caused this modification since FEPAGRO - RS 13 Vassourinha has shorter roots than the other two cultivars, resulting in less competition between plants. For biostimulator use, there was no difference between the cultivars; however, genotype FEPAGRO - RS 13 Vassourinhas had an average of 9.96 without it. For the interaction, density 10 had no effect from the biostimulator. Density 20 with the biostimulator presented the highest average, 7.26.
Density | FEPAGRO-RS 13 Vassourinha | BRS CS01 | Iapar 19 - Pioneira |
---|---|---|---|
10 | 8.78 aA | 7.39 aA | 8.38 aA |
20 | 8.92 aA | 7.12 bA | 5.10 bB |
Bioestimulator | |||
With | 7.34 aB | 7.13 aA | 7.20 aA |
Without | 9.96 aA | 7.39 bA | 6.28 bB |
Density | Bioestimulator | ||
With | Without | ||
10 | 6.46 aB | 6.92 aA | |
20 | 7.26 aA | 6.84 bB |
Means followed by the same lowercase letter in the row and uppercase in the column do not differ from each other by Tukey test (P<0.05).
When considering the protein content for density 10 with the biostimulator (Tab. 3), cultivar Iapar 19 - Pioneira had the highest average, 3.02. On the other hand, without the biostimulator, cultivars FEPAGRO - RS 13 Vassourinha and BRS CS01 presented the highee means, 2.84 and 2.6, respectively. At density 20, only cultivar FEPAGRO - RS 13 Vassourinha differed in the treatment with the biostimulator, presenting the lowest average. However, the others did not differ statistically. Comparing the different densities, there was a difference for FEPAGRO - RS 13 Vassourinha with the biostimulator, where the treatment with 20 buds per linear meter was superior. The same was repeated in the BRS CS01 cultivar. Cultivar Iapar 19 - Pioneira only differed in the treatment without biostimulator use, where density 20 presented the best average.
Density | FEPAGRO-RS 13 Vassourinha | BRS CS01 | Iapar 19 - Pioneira | |||
---|---|---|---|---|---|---|
With | Without | With | Without | With | Without | |
10 | 0.92 bAβ | 2.84 aAα | 1.53 bBβ | 2.60 aAα | 3.02 aAα | 0.85 bBβ |
20 | 1.29 bAα | 1.54 aBα | 2.36 aAα | 2.02 aAα | 2.11 aBα | 2.24 aAα |
Averages followed by the same lowercase letter compares cultivars for each density and hormone, uppercase compares hormones within the same cultivar and density, and Greek letter compares density within hormone, do not differ from each other by Tukey test (P<0.05).
For the lipid content (Tab. 4), cultivar FEPAGRO - RS 13 Vassourinha expressed the highest average at density 10 and at density of 20 buds per linear meter, the cultivars that stood out were FEPAGRO - RS 13 Vassourinha and BRS CS01. Cultivar Iapar 19 - Pioneira differed between the two densities, where density 10 presented the highest average. For the interaction between density and biostimulator use, the worst average was obtained in the treatment without biostimulator at density 20.
Density | FEPAGRO-RS 13 Vassourinha | BRS CS01 | Iapar 19 - Pioneira |
---|---|---|---|
10 | 0.7 aA | 0.58 bA | 0.56 aA |
20 | 0.64 aA | 0.66 aA | 0.37 bB |
Density | Bioestimulator | ||
With | Without | ||
10 | 0.62 aA | 0.63 aA | |
20 | 0.66 aA | 0.46 bB |
Means followed by the same lowercase letter in the row and uppercase in the column do not differ from each other by Tukey test (P<0.05).
The production environments that were efficient were favorable for the variable’s percentage of mineral material, starch, protein, lipid and productivity. The phenotypic stability analysis was used using the Annicchiarico method (Tab. 5) to identify which cassava genotype and cultivation environments showed greater efficiency. Density 20 with the biostimulant was favorable to the characteristics mineral and lipid material, along with density 10 with and without biostimulator for the latter. For starch and protein, only density 10 without the biostimulator was favorable. In terms of productivity, the favorable environment was density 10 with and without biostimulator.
Mineral material (%) | |||
---|---|---|---|
Cultivation environments | Mean | Environmental index | Classification |
Density 10* with biostimulant | 0.0215 | -0.000397 | U |
Density 10 without biostimulant | 0.0217 | -0.000208 | U |
Density 20 with biostimulant | 0.0227 | 0.000817 | F |
Density 20 without biostimulant | 0.0217 | -0.000212 | U |
Starch (%) | |||
Density 10 with biostimulant | 7.46 | -0.161 | U |
Density 10 without biostimulant | 8.92 | 1.30 | F |
Density 20 with biostimulant | 7.26 | -0.361 | U |
Density 20 without biostimulant | 6.84 | -0.776 | U |
Protein (%) | |||
Density 10 with biostimulant | 1.82 | -0.122 | U |
Density 10 without biostimulant | 2.10 | 0.152 | F |
Density 20 with biostimulant | 1.92 | -0.0190 | U |
Density 20 without biostimulant | 1.93 | -0.0107 | U |
Lipid (%) | |||
Density 10 with biostimulant | 0.622 | 0.0275 | F |
Density 10 without biostimulant | 0.635 | 0.0399 | F |
Density 20 with biostimulant | 0.663 | 0.0686 | F |
Density 20 without biostimulant | 0.459 | -0.136 | U |
Productivity (t ha-1) | |||
Density 10 with biostimulant | 0.622 | 4.82 | F |
Density 10 without biostimulant | 0.635 | 20.0 | F |
Density 20 with biostimulant | 0.663 | -8.63 | U |
Density 20 without biostimulant | 0.459 | -16.2 | U |
* 10 cassava buds per linear meter and 20 yolks per linear meter. Favorable (F), Unfavorable (U).
Therefore, in terms of favorable and unfavorable environments (Tab. 6) for percentage of mineral material, starch, lipid, protein and productivity, the best cultivar was BRS CS01, and the lowest performance was in FEPAGRO-RS 13 Vassourinha. The same pattern was repeated for mineral material and productivity for unfavorable environments, different from the results found for starch and lipid content, where the best was FEPAGRO-RS 13 Vassourinha, and Iapar -19 Pioneer was inferior. For protein, in unfavorable environments, the genotype with the best performance was Iapar 19 Pioneira. FEPAGRO-RS 13 Vassourinha showed lower results.
Favorable | Unfavorable |
---|---|
Mineral material (%) | |
BRS CS01 | BRS CS01 |
Iapar-19 Pioneira | Iapar-19 Pioneira |
FEPAGR-RS 13 Vassourinha | FEPAGR-RS 13 Vassourinha |
Starch (%) | |
BRS CS01 | FEPAGR-RS 13 Vassourinha |
Iapar-19 Pioneira | BRS CS01 |
FEPAGR-RS 13 Vassourinha | Iapar-19 Pioneira |
Protein (%) | |
BRS CS01 | Iapar-19 Pioneira |
Iapar-19 Pioneira | BRS CS01 |
FEPAGR-RS 13 Vassourinha | FEPAGR-RS 13 Vassourinha |
Lipid (%) | |
BRS CS01 | FEPAGR-RS 13 Vassourinha |
Iapar-19 Pioneira | BRS CS01 |
FEPAGR-RS 13 Vassourinha | Iapar-19 Pioneira |
Productivity (t ha-1) | |
BRS CS01 | BRS CS01 |
Iapar-19 Pioneira | Iapar-19 Pioneira |
FEPAGR-RS 13 Vassourinha | FEPAGR-RS 13 Vassourinha |
The genotype selected by the multi-trait genotype-ideotype distance index was FEPAGRO - RS 13 Vassourinha (Fig. 2), which was close to the cut-off point, or red line, which indicated the number of genotypes selected according to the selection pressure, which suggests that this genotype may present ideal characteristics for the study. Table 7 shows that the selection index for the five main components were retained; however, the analyzed data can be 100% explained by two factors when considering the first five main components and cultivar FEPAGRO - RS 13 Vassourinha, which was selected by the MGIDI index. The first factor explained the content of protein, mineral material, lipids and starch, and the second factor explained productivity.
Cultivar FEPAGRO - RS 13 Vassourinha showed an increase in starch content and productivity, which allows it to be evaluated in table 7 in relation to the other cultivars. The heritability factor showed that the lipid variable exhibited the highest genetic heritability, while protein was the lowest. The MGIDI provided desired (positive) gains for lipid and starch characteristics, and the other variables were all negative.
Factor | Variables | Mean | Sense | Heritability | MGIDI* |
---|---|---|---|---|---|
FA 1 | Protein | 1.94 | Decreasing | 0.00105 | -1.69·10-05 |
FA 1 | Mineral material | 0.0219 | Decreasing | 0.643 | -2.04 |
FA 1 | Lipid | 0.595 | Decreasing | 0.921 | 14.7 |
FA 1 | Starch | 7.62 | Increasing | 0.805 | 10.5 |
FA 2 | Productivity | 2.78 | Increasing | 0.615 | -0.22 |
*multi-trait genotype-ideotype distance index. Produtivity
Variance components and genetic parameters (REML) were estimated for three cassava genotypes cultivated in four environments in Rio Grande do Sul (Tab. 8). The phenotypic magnitude of the characteristics was related to effects attributed to the environment, which was part of the genetic variation. Thus, by establishing a joint relationship between individual phenotypic variance (σ²F) and genotypic variance (σ²G), it was possible to show that productivity, lipid, starch, protein and mineral material content were determined at 11.7, 49.3, 25.6, 0.00877 and 13.1% by genetic effects, respectively.
Components of variance REML1 | Characters | ||||
---|---|---|---|---|---|
Productivity (t ha-1) | Lipid (%) | Starch (%) | Protein (%) | Mineral material (%) | |
σ²G | 500 | 0.01 | 0.973 | 7.08·10-05 | 0.000000572 |
σ²F | 4,258 | 0.03 | 3.8 | 0.80 | 0.00000438 |
Ĥ²g | 0.12 | 0.49 | 0.26 | 8.77·10-05 | 0.13 |
Ĥ²mg | 0.62 | 0.92 | 0.80 | 0.00105 | 0.64 |
řgǧ | 0.78 | 0.96 | 0.90 | 0.03 | 0.80 |
CVg (%) | 8.04 | 19.1 | 13.00 | 0.43 | 3.46 |
CVe (%) | 22.1 | 19.3 | 22.1 | 46.3 | 8.92 |
CVr | 0.36 | 0.99 | 0.59 | 0.00936 | 0.39 |
1σ²G: genotypic variance; σ²E: residual variance; σ²F: individual phenotypic variance; Ĥ²g: broad-sense heritability for total genotypic effects; Ĥ²mg: heritability of the genotype mean; řgǧ: accuracy for genotype selection; CVg (%): genotypic coefficient of variation; CVe (%): residual coefficient of variation; CVr: rate of the coefficient of variation.
Individual phenotypic variance (σ²F) contributed to the yield and starch characteristics, where σ²F takes into account the genotype × environment interaction. Removing σ²G influenced the environment for the characteristic. Broad-sense heritability is related to the percentage of genetic variance existing within the phenotypic variance (Ramalho et al., 2012). Several recent studies have focused on heritability estimates for maize (Ferrari et al., 2022), wheat (Carvalho et al., 2019), rice (Facchinello et al., 2021), white oat (Rosa et al., 2021) and soybean (Barbosa et al., 2021), contributing to the genetic gain of these crops.
In this context, broad-sense heritability for the total genotypic effects without interference of genotypes × environments interaction (Ĥ² g) showed el formato dehigher magnitudes for lipids (0.49) and starch (0.26). Mineral material (0,13), productivity (0,12) and protein (8,77e-05) had the lowest Ĥ²g. The broad-sense heritability of the genotype mean (Ĥ² mg) was high for lipids (0.92) and starch (0.80). High accuracy values demonstrated high experiment precision, with efficiency in the methods of selection and genetic increment of the characteristics (Costa et al., 2000). High accuracies (řgǧ>0.75) were obtained for all characteristics, except for protein, which showed low accuracy (0.03).
The genotypic variation coefficient (CVg) expressed greater magnitudes for lipids (19.1%), starch (13%) and productivity (8.04%), indicating a greater genetic variability of genotypes for these traits. The residual coefficient of variation (CVe) was higher for protein (46.3%) because this characteristic was influenced by the interaction between genotypes x environments.
In general, the optimal procedure for verifying genetic values is the best unbiased linear predictor (BLUP) as it allows understanding and selection of promising genotypes with information that represents the true genetic value and minimizes the interference of estimates by environmental effects (Borges et al., 2010).
Figure 3 shows that cultivar BRS CS01 had higher genetic value for productivity. For lipids and starch, cultivar FEPAGRO - RS 13 Vassourinha showed the highest genetic values. For protein content of the roots, there was no genetic difference between the cultivars. On the other hand, for mineral material, cultivars Iapar 19 - Pioneira and BRS CS01 showed the highest averages.
The methods used by breeding to select certain genetic characteristics show the ideotype with a designated function, either industrial or for both human and animal food. The better methods include descriptive analysis, Annicchiarico, MGIDI index, unbiased linear predictor, REML, and BLUP, which, through desired characteristics productivity and starch content, indicate the best genotype.
CONCLUSION
Cultivar BRS CS01 showed the highest yield and concentration of mineral material, genotype FEPAGRO - RS 13 Vassourinha had the highest lipid content, and Iapar 19 - Pioneira showed the highest protein concentrations.
The test of comparison of means and MGIDI index showed that cultivar FEPAGRO - RS 13 Vassourinha had the most starch and was the ideal cultivar based on multi-characteristics.
Density 10 with the biostimulator was favorable for productivity and lipids, whereas density 10 without the biostimulator was favorable for starch, lipids, proteins and productivity. Density 20 with the biostimulator was favorable for lipids.
The lipid and starch characteristics exhibited greater genetic contributions, with a heritability of 49 and 26%, respectively