Agronomic Traits of Seven Accessions of Cenchrus purpureus under Rainfed Conditions in the Tropic Region of Veracruz, México

1. Introduction

Fodder production for cattle feeding occupies a third of the Earth’s surface, which is important for its’ contribution to food security [1]. 25% of the total protein, 18% of the total calories [2], and 34% of the beef [3] consumed at a global level are generated in that space.

The main use of agricultural land in Latin America is cattle rearing [4]; in Mexico, animal husbandry occupies more than 55.9% of the land area, and 29.3% of the total agricultural land area [5]. The largest farming system is focused on meat production, which is concentrated in 10 Mexican states including Veracruz [6]. The Huasteca Alta Veracruzana (HAV) region is representative of the cattle rearing activity. This region concentrates more than 162,000 cattle production units in 3.7 million hectares, with a total of 5.8 million bovines [7]. In HAV, cattle feeding is based on fodder grasses and pastures, which are considered pillars of animal husbandry in Mexico and all over the world [8] due to their efficiency for the net production of proteins for human consumption [9]. Fodder quality and production rates have been affected by climate change, and this has fostered the development of research that seeks to determine the impact of the interaction between the modified climatological conditions and the different areas of animal husbandry, but this has not been satisfactorily achieved so far [10,11,12]. The efficient use of grazing lands is key for the sustainability and productivity of the cattle farming sector; however, to achieve this, commonly used fodder grasses must be evaluated to identify those that were able to adapt to the current climate conditions. It has been proved that improved grasses have a greater biomass production potential, which makes them viable options to replace local species, or those considered native or traditional [13].

Cenchrus is a highly diversified genus [14], with around 145 extant species [15]. Its productive potential in conditions of adequate irrigation and fertilization has been proven [16], as in the case of Elephant grass, a species characterized by its’ large biomass production [17]. This is a species native to Africa and it has been introduced successfully to several tropical and subtropical countries due to its high capability to adapt to different agro-ecological systems, in which it has been used as fodder crop [18]. However, weather conditions in the tropics comprise well-defined rain and drought periods, which affect crop productivity [19].

Fodder production and quality must be carefully evaluated and, based on this evaluation, take proper management decisions to improve the productivity of cattle farming units, and to adapt to the conditions provoked by climate change [20]. These modified environmental conditions (climatic variability) have not been thoroughly investigated, nor their global impact evaluated [12]; in the case of the HAV, the months classified as abnormally dry have increased to 10 from 2016 to 2021, and the region suffered from extreme drought in some years of this period. [21]. Thus, this study aimed to evaluate the agronomic performance of seven accessions of Cenchrus purpureus under the current climate conditions of the HAV.

2. Materials and Methods

2.1. Description of the Experimental Site

2.1.1.

The research was carried out at the experimental field of the Instituto Tecnológico Superior de Tantoyuca, in the municipality of Tantoyuca, state of Veracruz, Mexico, located in the HAV region (Figure 1), between parallels 97° 59’ and 98° 24’W, and 21°06’ and 21°40’N, at an altitude up to 300 masl. This region has a warm subhumid climate, with rainy summers, a mean annual temperature between 22 and 26°C, average annual rainfall between 1100 and 1300 mm, and a relative humidity of 44% [22]. The characteristics of the soil in the site where the pre-study experiment was performed are the following: argillaceous texture, a pH of 7.8, light salinity, 5 % nitrogen content, 50 ppm phosphorus, 49 ppm potassium and organic matter (>5.1 %), according to the analysis performed at the soil laboratory of the Universidad Veracruzana.

2.1.2. Evaluated Germplasm and Experiment Design

The vegetal material used for this study belongs to the Cenchrus purpureus species, and the evaluated accessions, taken as treatments, were King grass, Maralfalfa., Elephant grass, Merkeron, Purple grass, CT-115 and Taiwan, obtained from the Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias at the La Posta experimental field, in the town of Paso del Toro, Veracruz. A randomized block experimental design with seven treatments and four replications, which were evaluated throughout a 12-month period under rainfed conditions.

2.1.3. Experiment Site and Sowing

A 1,222 m² area was used to create 28 plots of 25 m² (5 m x 5 m), distributed in 4 blocks (7 plots p/block), and a distance of 2 m between plots and blocks. Sowing of the assessed accessions was done with two stems (phytomers with three internodes) per sowing point, introducing two internodes in the seedbed, and placing the stems at a 45° angle. The distance between crop rows and sowing points was 80 cm, obtaining a density of 36 tillers in the 25 m² of every experimental unit. The seeded ground was only prepared by primary tillage (plowing and harrowing). Seeding was done in September 2018, and the cut of homogenization was made on July 2019, at 15 cm above soil level.

2.1.4. Sampling and Cutting Frequency (CF)

A sample area of 2 m² (1 m x 2 m) was obtained leaving 2 m of border of each of the experimental units to avoid the edge effect. The cutting frequency was determined considering a mean intercepted radiation rate of ≥ 95 in the four replications of each treatment [24]. Four consecutive cuts were performed from July 2019 to April 2020.

2.1.5. Evaluated Variables

Green matter yield (GMY) and dry matter yield (DMY). Cuts were done at 15 cm above soil level for GMY purposes, and all the biomass of the usable sample (2 m²) was weighed immediately after cutting, using a digital scale with a minimum precision of 1 g 5 whole plants (leaf and stem) were randomly selected to determine the DMY to obtain a subsample of 300 g, which was dried in a forced air lab oven at 65 °C until achieving a constant weight [25]. The obtained results were extrapolated to 10,000 m² (ha).

Leaf:stem ratio (LSR) and accumulated dry leaf matter rate (ADLMR). 5 plants were selected from the biomass obtained from the sample, and the leaves were separated from the stem. Once separated, the stems were all placed in a forced air lab oven at 65 °C until achieving a constant weight. The LSR was estimated by dividing the dry weight of the leaf samples between the dry weight of the stem samples. The ADLMR was calculated dividing the dry leaf matter yield between the number of days passed from cutting period to cutting period. To determine the values of the GMY, DMY, LSR and ADLMR variables per season (Dry and Rainy), the value of each variable per accession obtained in the cutting of each period were added up, and this was added up to the measurements of the variables of each accession obtained within a year to calculate the annual production [26].

Leaf area index (LAI) and height (H). The LAI was determined with five randomly selected plants taken from the sample; afterwards, the width and length of the fourth leaf of each selected plant were measured to use the method proposed by Herrera [27]. The HI was determined using a measuring tape, calibrated in meters, centimeters and millimeters, to measure five plants of each sample at every cutting period, from soil level to the last leaf.

2.1.6. Climate Variables

Some climate variables that are important for biomass production, such as average rainfall (mm), relative humidity (%), global radiation (W/m²), mean temperature (°C) and reference evapotranspiration (mm), were taken during the development of the research activities (Figure 2). These data were obtained from the weather station located at the experimental field of the Instituto Tecnológico Superior de Tantoyuca. For the purposes of this research, “dry season” is defined as the months of the year in which reference evapotranspiration is greater than the average rainfall because, in these conditions, crops are in a state of stress due to a lack of moisture, and this period is considered key to determine the need for irrigation [28,29].

2.1.7. Statistical Analysis

Analysis of variance (ANOVA). The Statistica v. 10 software was used for statistical analysis; the general linear model (GLM) was as follows:

Where:

Yijk = the quantitative response variable of the i-th accession and the j-th replication.

µ = Overall mean

αi = fixed effect of the i-th C. purpureus accession

= effect of the j-th block

εijk = random error associated to each observation, where ij∼NI(0,σ²), normality was analyzed with the Shapiro-Wilk test; and homoscedasticity, with Bartlett’s test (p < 0.05). The means were compared with the Fisher method (p < 0.05).

Multivariate analysis. Using the studied response variables from each cutting performed, and the monthly average climatological variables of the research site, a multiple linear analysis was executed to ascertain the weight climatological variables have on the assessed accessions. The MLR model is as follows:

$${{\rm Y}}_{\iota }={\beta }_0+{\beta }_1{X}_1+{\beta }_2{X}_2+{\beta }_3{X}_3+{\beta }_4{X}_4+\mathcal{e}\mathfrak{i}$$

Where:

${{\rm Y}}_{\iota }$= Dependent variable of the i-th accession

${\beta }_0=Intercept$

$X_1=Accumulatedprecipitation$

$X_2=Relativehumidity$

$X_3=Globalradiation$

$X_3=Globalradiation$

$\mathcal{e}\mathfrak{i}=Randomcomponent$

The assessed accessions were typified with the obtained agronomic characteristics and yield indicators [30], through multivariate analysis (K-means clustering), Dendrogram (Graph clustering) was done through complete Euclidean distance transform, amalgamation was done through the unweighted pair-group method, and normalized through the 1-Pearson correlation coefficient [Correlation distance/maximum distance]*100.

Lastly, Principal Component Analysis were used to observe the possible grouping of the agronomic and yield variables, as well as the accession grouping [31].

3. Results

The Taiwan accession showed a greater LAI in the rainy season (p < 0.05) than King grass, Elephant, and Purple grass accessions, and similar to the rest. Regarding the ADLMR, the Merkeron grass accession showed the highest value, statistically different (p < 0.05) only to the Taiwan and King grass accessions, thus demonstrating its capability to accumulate dry matter in the leaves; consequently, this accession showed a LSR of 0.88, the greatest in all the accessions, showing a statistical difference (p < 0.05) with the CT-115 accession, which had the lowest value (0.58). CT-115 was the accession with the greatest capability to accumulate dry matter on the stems, which is made patent by the fact that it is one of the tallest accessions, only lower than King grass and Maralfalfa. An important characteristic from the point of view of productivity is the cutting frequency (CF), that is, the time (days) it takes for the plant to reach maturity and, consequently, harvest time. Merkeron, Elephant and Purple grasses were harvested at intervals ∼60 days, while the rest of the accessions took ∼80 days (p < 0.05) (Table 1).

CT-115 showed the lowest LAI off all the accessions (p < 0.05) during the dry season; on the other hand, Elephant grass had the highest ADLMR, showing a statistical difference (p < 0.05) with Merker grass, Taiwan A-144, King grass and CT-115. The highest LSR was observed in Purple grass (p < 0.05); nevertheless, this accession had the lowest H, showing notable differences with Elephant grass and Taiwan. There was little variation in CF, with Elephant grass having the lowest value, like King grass and CT-115 (p > 0.05) and showing statistical differences with the rest of the accessions (p < 0.05). Merkeron grass was the accession with the most days of cutting frequency of all (p < 0.05) (Table 2).

Table 3 shows that the LAI showed no statistical difference (p > 0.05) among the accessions. However, it is important to highlight the ADLMR variable in Elephant grass, which had a value of 228.51 kg/ha/day, the highest in all treatments, showing statistical differences (p < 0.05) with Merkeron grass, Taiwan, King grass and CT-115 and a greater capability to deposit dry matter on the leaves. Notwithstanding the fact that Elephant grass showed and adequate LSR, it was surpassed by Purple grass, which was statistically different (p < 0.05) from King grass, CT-115 and Maralfalfa. The H obtained was not statistically different (p > 0.05) among the assessed accessions. The shortest CF was 71 days, corresponding to Elephant grass, while the longest corresponded to the Taiwan accession, showing statistical differences (p < 0.05).

The differences in the climate parameters between the dry and rainy seasons (Figure 2) had an influence on the agronomic behavior of the assessed accessions. In the rainy season, LAI, ADLMR, H and CF noteworthy values were reported. In dry season, on the other hand, the LSR was greatest (p < 0.05).

The GMY of the assessed accessions was not statistically different during the rainy season; nevertheless, the DMY results indicate that CT-115 was quantitatively superior to all the other accessions, and statistically different (p < 0.05) from Elephant and Merkeron grass. Conversely, the GMY of Maralfalfa was greatest in the dry season, and it was statistically different from CT-115 and Taiwan (p < 0.05). In turn, Elephant grass had the highest accumulation of dry matter (t/ha), followed by Purple grass and Maralfalfa (p > 0.05) in the same period. The annual cumulative green matter yield was greatest in Maralfalfa, and statistically different (p < 0.05) from Taiwan, while Elephant grass had the greatest amount of accumulated dry matter of all the accessions after a year of assessment (Table 4).

Table 5 shows the average GMY and DMY per cutting period. The capability of Elephant grass of producing dry matter (DM) under the particular climatological conditions of the HAV is evident. There are statistical differences between GMY and DMY (p > 0.05) in both seasons. The greatest values for these variables were reported in the rainy season.

The correlation analysis achieved to determine the degree of association between the DMY with others agronomic variables indicated positive correlations with LAI, H, ADLMR, and negative correlation with LSR. Nevertheless, it was found that there is little association with the GMY variable (r = 0.147), with a value of p= 0.083. In the case of the H and LSR variables, the association was low (r = 0.192 and 0.262, respectively). There is an intermediate value association with the LAI variable (r = 0.572). There was a high degree of association with the ADLMR variable (r = 0.889), and a value of p= 0.001 (Table 6).

The MLR implemented for the creation of statistical models to describe the association degree between climatological and agronomic variables of C. purpureus in the HAV is presented in Table 7. From all the performed analysis, the one that obtained a medium to high value in the goodness-of-fit indicators (r= 0.73 and r²a=0.52) was the GMY variable, where the Forward method included all the climatological variables. The goodness-of-fit indicators for the DMY mathematical model (r = 0.22 and r²a=0.03) do not indicate a good fit of the data, and likewise with the LAI, LSR and ADLMR. The MLR model (Forward method) for H variable showed a medium-high degree of association (r= 0.70 and r²a = 0.48), only including temperature and global radiation.

The implemented cluster analysis formed three accession groups according to their agronomic descriptors. Group 1 (G1) includes Elephant and Purple grass, and it is characterized by having the highest DMY, LSR and ADLMR values. Group 2 (G2), comprising Maralfalfa and Merkeron grass, stands out for its high productive capability with regards to GMY. Finally, Group 3 (G3), made up by CT-115, King grass and Taiwan, was characterized by the lowest records of all the indicators used for this analysis (Table 8 and Figure 4).

Figure 3. Dendrogram, clustering different C. purpureus accessions according to their agronomic and yield descriptors in the Huasteca Veracruzana region. (a). Dry season. (b). Rainy season.

Figure 3. Dendrogram, clustering different C. purpureus accessions according to their agronomic and yield descriptors in the Huasteca Veracruzana region. (a). Dry season. (b). Rainy season.

The principal component analysis, completed with all the estimated variables (rainy and dry seasons) of the seven C. purpureus accessions, allowed the creation of two groups of orthogonally independent variables. The first group (Factor 1) corresponds to those known as agronomic characteristics, and the second group (Factor 2) is made up of estimated yield variables (Figure 4).

Figure 4. Grouping of C. purpureus accessions according to two principals’ components: Factor 1. Agronomic characteristics and Factor 2. Yield parameters. (a) and (c). Dry season. (b) and (d). Rainy season.

Figure 4. Grouping of C. purpureus accessions according to two principals’ components: Factor 1. Agronomic characteristics and Factor 2. Yield parameters. (a) and (c). Dry season. (b) and (d). Rainy season.

The analysis implemented in the dry and rainy seasons (Figure 4), formed three groups of Cenchrus purpureus accessions, a) outstanding agronomic and yield characteristics, b) outstanding yield characteristics, and c) unremarkable agronomic and yield characteristics. The accessions integrating each group in both seasons are not the same; their genetic characteristics and interactions with the environment define their performance.

4. Discussion

The LAI is an important parameter for the development and production of biomass in a grass crop, because it partially depends on their photosynthetic capability [31]. Taiwan was the accession that had the highest LAI of the assessed accessions. The different climatological conditions in the rainy and dry seasons were decisive in the modification of this agronomic characteristic. The LAI value of CT-115 and Maralfalfa doubled in the rainy season with respect to the dry season (Table 1 and Table 2), which was also reported by Maranhão et al. [32], upon evaluation of a species of the Cenchrus genus in eastern Brazil along three climatologically defined periods (rainy, transitional, and dry) and different cutting ages. In their research, they found that LAI was two times greater in the rainy season, in comparison with the dry season; additionally, according to their results and analysis, they determine that there is a positive correlation with age. Seemingly, the high LAI favor apical growth due to the shading effect caused by the elongation of the stem, as mentioned by these authors, which resembles the results of the research presented in Table 1.

The ADLMR of the Taiwan accession on the rainy season was greater than reported by Vinay-Vadillo et al. [33], in the central region of Veracruz state. Despite the application of different doses of fertilizers, the values that they obtained were lower than those reported in this study. The sowing density remained the same; thus, the greater accumulation of DM/ha, per day, may be attributed to the longer CF, and to a greater acidity of the soil where the experiment was carried out, which negatively affected yield. ADLMR was led by Elephant grass, Purple grass and Maralfalfa (Table 3); these results were corroborated by Reyes-Pérez et al. [16], who reported that these same accessions have a greater potential for this production indicator and have the capability of increasing it with the application of fertilizers, and a greater cutting age.

The LSR results of the evaluated accessions showed greater values in the dry season, in comparison with the rainy season. Purple grass stood out from the rest with regards to this variable, showing a value of 1.03 (at 80 cutting age). This result does not coincide with the results reported by Vinay-Vadillo et al. [33], who evaluated this parameter at 45 days regrowth, however, they found that, upon harvest at the 60th day, LSR values decreased in the cultivars submitted to irrigation with and without nitrogen fertilization, showing that the stem proportions in this grass increase with age. The LSR decreased due to the increase of the stem proportion, mainly in the rainy season, which favors apical growth. This was also described by Habte et al. [26] upon evaluating 45 genotypes of C. purpureus in the dry and rainy seasons of Bishoftu, Ethiopia, at 1890 masl. However, it has been shown [354] that sowing density plays an important role on LSR of C. purpureus cultivars, in which values of 1.6 were obtained when the number of plants per hectare was 11,350. For the purposes of this research, plant density per hectare in the research units was established at 15,625, which could partially explain the lowers values found in the assessed accessions. The increase of the stem proportion in this genus reduces fodder quality in terms of their neutral detergent fiber content, as indicated by Ledea-Rodríguez et al. [35].

In the case of King grass harvested at 80 days, the LSR was 0.78 and 0.89 in the rainy and dry season (Table 1 and Table 2). Seemingly, this result is related to an increased availability of water, because it differs to the reports of Reyes-Pérez et al. [16], who reported a LSR value of 1.0 of this accession at 45 days of regrowth, with and without fertilizers, at different doses and concentrations. However, when harvested after 60 days, the value of this parameter oscillated between 0.66 and 0.70, which led to the conclusion that there is a negative correlation between the LSR and yield, which contrasts with the results obtained in this study, which found values of r= 0.262 and r²= 0.069 between the LRS and DMY (Table 6).

King grass is referred to in this study as one of the accessions with the greatest H value per cutting period (2.00 m). This result differs with the reports of Reyes-Pérez et al. [16], who, upon application of fertilizers at several concentrations, managed to increase this agronomic characteristic, surpassing their control treatment (without fertilizer) at 45 and 60 days. The heights obtained from the non-fertilized accessions that they assessed were generally lower than those reported in this study, which suggests a favorable response to this kind of agricultural practices. On their part, Cruz et al. [36], upon evaluating new clones generated from CT-115, reported that the height parameters surpassed those of the parent plants despite having been established in moderately saline soils, and in drought conditions. This is similar to what is reported in this study, where the height of the assessed accessions differs numerically (p > 0.05). Plant height in different C. purpureus genotypes partially depends on regrowth days; thus, the older the plant is, the taller it is [37], which could explain the mentioned difference.

The accessions with the greatest average GMY’s per cutting period were Maralfalfa and CT-115 (243 and 220 kg/ha, respectively), at CF of 83 days. These results are higher to those presented by Reyes-Pérez et al. [16], who assessed five C. purpureus accessions in a humid subtropical region of South America. They established short cutting periods, and used a greater sowing density, to which the difference is attributed. Nevertheless, these same authors coincide in that Maralfalfa obtained the highest value for this variable, demonstrating its capability to accumulate green matter Moreover, they reported that, assessed conditions of irrigation without fertilizers, CT-115 showed the lowest GMY when harvested at 45 days of regrowth (40.63 t/ha). This result is lower than that obtained with this same accession in this study, despite the fact it was assessed in unirrigated conditions. The difference in results may be related to the cutting frequency because the latter was above 40 days. Encountered yields were lower in the dry season, showing significance in them (p < 0.05). This same effect was reported by Habte et al. [26] and Tulu et al. [38], who associated the high GMY and DMY to the months of highest rainfall. Moreover, they found no significant differences between the genotypes they assessed, unlike the reported results here (Table 5).

The DMY is a desirable attribute in fodder crops, because there are nutrients in this fraction, which could be beneficial to ruminants. The results presented showed that some accessions have a greater capability to deposit dry matter, and the variation of the climatological variables between the rainy and dry seasons did not affect the behavior of this indicator in any considerable degree, which is the case for Elephant grass, Purple grass and Maralfalfa (Table 4). However, it was observed that DMY of the rest of the assessed accessions was affected by the climatological conditions of the rainy and dry seasons (p < 0.05). This indicates that these genotypes are susceptible to climate fluctuations, as demonstrated in an experiment carried out by Ledea-Rodríguez et al. [39] in Granma province, Cuba, where they evaluated three new C. purpureus accessions (CT-601, CT-603 and CT-605), conceived to tolerate drought conditions, in two periods (rainy and lightly rainy); however, they demonstrated their lower DM accumulation capabilities in the period of least rainfall, as reported in this study. Taiwan showed greater DMY in the rainy season, duplicating the amount shown in the dry season. Vinay-Vadillo et al. [34] evaluated the effect of fertilization on DMY of the Taiwan cultivar during the rainy season, which led them to the conclusion that the application of agrochemical agents helps to duplicate this variable; however, their yields do not surpass those reported in this study, possibly due to the shortened CF. The photosynthetic and DM accumulation capabilities show variation within the Cenchrus genus, the higher productive yields of some accessions in the seasons of most abundant rainfall are due to their intrinsic potential and the higher amount of nutrients dissolved in the ground because of the greater humidity, as illustrated by Habte et al. [26], upon evaluation of the DMY and GMY of several genotypes of this genus in similar conditions.

Among the assessed accessions, DMY per cutting, per year, was greatest in Purple grass and Maralfalfa (Table 5). In this respect, Reyes-Pérez et al. [16] mention that these cultivars showed the highest DMY, above king grass, CT-115 and Elephant grass, and they highlighted that fertilization and longer CF are factors that increase the deposit of DM [37,40,41] in all the assessed genotypes. This explains the greater DMY found in this study.

Most of the evaluated agronomic variables showed a low association degree, it was positive, low, and not significant for DMY and GMY (Table 6). With this result, it is possible to infer that the assessed accessions accumulate variable amounts of moisture on the stem, which occupies a greater proportion of the plant (Table 3). On the contrary, Habte et al. [26] found high positive correlations (0.86) between the total amounts of DMY and GMY, accumulated in several C. purpureus cultivars. Seemingly, the higher proportion of leaves (LSR 2.41 and 5.87 in the rainy and dry seasons, respectively) boosted this correlation. Moreover, these authors coincide in the results of the correlation between DMY and LSR (negative correlation), indicating that, the higher LSR is, higher is the total proportion of leaves in the plant and, consequently, the amount of DM represented by the stem is lower.

The data that best fit the MLR model correspond to GMY variable, which is shown by their goodness indicators (Table 7); Araujo et al. [42], report similar values in these indicators for several species of C4 grass crops. P. maximum (Tanzania and Mombasa), grass species with an erect growth habit, like those of the C. genus, reached r² values between 0.78 and 0.87 in the univariate linear regression analysis in the ADLMR and ambient temperature variables. The step-by-step execution with the Forward method used in this investigation for the multivariate regression model included all the used climatological variables; from these, it can be observed that temperature has a negative angular value on biomass production. Andrade et al. [43] found that temperature is related to enzymatic physiological processes, such as photosynthesis, and other processes related to transport through cellular membranes. The production behavior of the assessed accessions may have been affected by the extreme climatological conditions under which this research was done (Figure 2).

Clustering analysis used ADLMR variable as the most significant in the configuration of the groups. This variable demonstrated importance for its significant value in the research developed by Ray et al. [30], in which they assessed 11 drought-resistant clones of C. purpureus. In their results, they reported that the clustering analysis formed three groups, which they describe by their biomass production yield as: promising, medium and low. The clustering and characteristics of the groups formed by these authors coincides with the reports of this study, and this confirms the genetic and productive variability that exists within the C. genus, as cited by Olivera-Castro et al. [14].

5. Conclusions

In general, the assessed accessions showed variations on their agronomic behavior and yield indicators, which are attributed to changes in the climatological variables between the rainy and dry seasons. Despite this, Elephant gras and Maralfalfa achieved the highest DMY and GMY in both seasons, even though there were minimal differences between them, proving their genetic capability to maintain a stable production in different environmental conditions. It is necessary to simultaneously evaluate the agronomic characteristics, yield and quality indicators of the assessed accessions in order to identify those who show the best levels in these attributes.