Determination of the number of harvests to select elite sugarcane genotypes in NE Brazil
Received: December 14, 2017
Accepted: January 08, 2018
Published: January 20, 2018
Genet.Mol.Res. 17(1): gmr16039869
Sixty-eight sugarcane genotypes were evaluated to determine the number of harvests required to select new genotypes in the sugarcane belt in the state of Pernambuco, northeast Brazil. Using a random block design, the following parameters were evaluated: (i) polarizable sugars (POL) per hectare, in metric tons (TPH); (ii) culm productivity per hectare (TCH); (iii) fiber content (FIB); (iv) adjusted percent POL (PCC); (v) soluble solids (BRIX); (vi) total recoverable sugars (ATR); and (viii) metric tons of ATR per hectare (ATR t.ha-1). A variance analysis and genetic parameter estimation were carried out. Means were analyzed using the Scott and Knott test. The repeatability coefficient and the number of harvests were determined using analysis of variance, principal component analysis, and structured data analyses. The best cultivars were SP79-1011, RB952692, RB952675, RB813804, SP78-4764, RB952522, RB952511, RB953265, RB952754, RB952875, SP80-1816, RB763710, and RB892575. Two evaluations are enough to select elite genotypes in the early experimental stages, reducing the time-to-market of these cultivars by three years under the edaphic and climatic conditions in the sugarcane belt evaluated
Sugarcane holds considerable importance in Brazil’s economy, accounting for approximately 2% of the country’s gross domestic product (Biosev, 2013). Recent biotechnological breakthroughs have afforded to explore the biological potential of sugarcane more fully. For example, current research has shown that the plant species may be used to produce biodiesel, in a process based on a genetically modified yeast in a partnership between the company Amyris and the São Martinho Group (Única, 2013). Also, Coelho et al. (2001) described a microbiological process to produce a sugarcane biopolymer with interesting healing properties both for domestic animals and for humans with skin lesions or submitted to surgery. However, sugar and ethanol remain the core sugarcane products today due to the best return on interest they afford.
As with any industrial process, quality of sugarcane raw material is essential in the generation of energy as well as the production of sugar, ethanol, and any other item based on this plant species. In this sense, the quality of sugarcane cultivars is the starting point in estimating how successful an enterprise will be, and is the most important and least expensive technological variable from the producer’s standpoint (Barbosa and Silveira, 2012).
Several factors determine sugarcane productivity when manufacturing sugar, ethanol, and other sugarcane-based products. These parameters include the management, planning, and implementation of up-to-date agricultural technologies, in addition to using agricultural inputs reasonably. In this sense, more productive and resistant cultivars with favorable traits developed by genetic improvement programs are central to obtain higher yields at lower costs and better quality. For Silva (2008), the ideal sugarcane cultivar is the one that affords high productivity even under harsh environmental conditions.
More specifically in Brazil, the state of Pernambuco (PE) stands out as the second largest sugarcane producer in the country’s northeast region, with approximately 14.90 million metric tons produced to meet the demands of the sugar and energy industries (CONAB, 2013). The main hurdle to increasing productivity is the interaction genotype × environment, which is influenced by specific soil characteristics such as altitude, for instance, and more particularly so by irregular rainfall patterns marked by long droughts.
In their pioneering study, Koffler et al. (1986) characterized the sugarcane belt in PE. The authors systematically gathered a whole set of relevant environmental information with a view to widening the horizons of future genetic improvement research on sugarcane. With that in mind, they split the sugarcane belt in PE into five micro-regions, namely North Woodlands (NW), South Woodlands (SW), Mid Belt (MB), North Coast (NC), and South Coast (SC), and described the geology, geomorphology, climate, hydrology, natural plant cover, soils, and ecological zoning of each.
Subsequently, edaphic and climatic aspects of these micro-regions were investigated considering the suitability to grow sugarcane. It was found that the occurrence of distinctive environments is behind the specific limitations affecting the interaction of soil and climate factors with sugarcane productivity. In genetic improvement terms this means that the agro-industrial performance of a given sugarcane cultivar in a region may not match that in another, and that the environment may influence the expression of desired traits in plants. Interestingly, climate characteristics also vary across harvest seasons, and may shorten the lifespan of ratoons. As a result, sugarcane plantations require renewing at comparatively short intervals, increasing productivity losses.
For this reason, the samples used in cultivar assays at the final stages of an experiment are based on the means of three to four harvests (Koffler et al., 1986). These experiments compare the performance of new clones with that of cultivars more widely grown (Ferreira et al., 2005), identifying high-productivity and high-longevity cultivars in ratoons, as carried out by Melo et al. (2009). The authors assessed the agronomic performance of sugarcane genotypes in SC (PE) using both univariate and multivariate analysis of variance (ANOVA) of experiments and relevant genetic parameters, managing to characterize genotypes that performed better both in the field and in industrial processes. However, these analyses do not suffice to identify the genotypes with high longevity potential in ratoons.
The objective of genotype selection is to increase a genotype’s agroindustry performance throughout its lifetime, and to make sure that the good performance of an individual plant in certain traits reflects the genetic potential of the cultivar as a whole (Cruz and Regazzi, 2001). More specifically in sugarcane, the expectation lies in more long-lived ratoons, that is, high productivity in the plant crop and the 1st, 2nd, and 3rd ratoon crops. In this context, a feasible alternative in the identification of more long-lived genotypes includes the use of the repeatability coefficient (r). The coefficient is useful to determine the number of harvests required to select the best genotypes, significantly reducing costs—since cultivar selection assays are costly and time-consuming. In addition,r also helps shorten the time-to-market of new cultivars from the 12 to 15 years on average required Currently, several methods are available to estimate r, such as ANOVA, principal component analysis (PCA) (Abeywardena, 1972) and structured data analysis (SA) (Mansour et al., 1981).
Ferreira et al. (2005) evaluated sugarcane genotypes grown in the production belts in the Brazilian states of São Paulo, Minas Gerais, and Goiás. Based on r, the authors discovered that three harvests are enough when selecting new genotypes with predictability degree above 80% of actual values. Santos et al. (2004) evaluated sugarcane genotypes in the state of Alagoas, Brazil, and concluded that the five harvests afford the selection of genotypes with predictability above 80% of actual values.
These studies are especially relevant in the context of sugarcane production in PE, Brazil, since no research covering the whole state’s sugarcane belt has been published. So, the objectives of the present study were to evaluate the agro-industrial performance of 68 sugarcane genotypes and to determine the number of harvests needed to select elite genotypes for commercial growth of the plant species in the sugarcane belt in PE, Brazil.
Materials and Methods
General description of the experiments
The experiments were conducted in plantations managed by the sugarcane mills participating in the Sugarcane Genetic Improvement Program, Federal Rural University of Pernambuco, Brazil, which is a affiliate of the Inter-university Network for the Development of the Sugar and Energy Sector (PMGCA/UFRPE/RIDESA): Usina São José, Usina Trapiche, Usina Central Olho d’Água, Usina Pumaty, and Usina Petribú. According to the classification system proposed by Koffler (1986), these mills are located in the sugarcane producer micro-regions North Woodlands (NW), South Woodlands (SW), Mid Belt (MB), North Coast (NC), and South Coast (SC), respectively, in the state of PE, Brazil.
A four-repeat randomized block design was used. A quadrat was defined as five 8-m lines 1 m apart with sugarcane grown using two culms placed next to each other, though the end of the first culm was next to the head of the second (Stolf, 1986). Soil pH adjustments and fertilization were carried out following the procedures used by the sugarcane mills cited above. The parameters evaluated were (i) polarizable sugars (POL) per hectare, in metric tons (TPH); (ii) culm productivity per hectare (TCH); (iii) fiber content (FIB); (iv) adjusted percent POL (PCC); (v) soluble solids (BRIX); (vi) total recoverable sugars (ATR); and (viii) metric tons of ATR per hectare (ATR t.ha-1).
Productivity per area unit (TCH) was estimated weighing all culms of a quadrat (kg) according to the formula:
TCH = total sugarcane weight (kg) × 10/quadrat used area (m2)
In turn, TPH was calculated using the formula: TCH × PCC/100
Also, ATR t.ha-1 was calculated using:
ATR t.ha-1 = TCH × ATR/1000.
BRIX was obtained using a refractometer and a homogenized sample of sugarcane juice obtained from 10 culms randomly collected in the quadrat. FIB, PCC, and ATR were calculated according to Fernandes (2003).
Genotypes evaluated and number of measurements
Four measurements were carried out (plant crop and the 1st, 2nd, and 3rd ratoon crops) in micro-region NC using the genotypes RB72454, RB763710, RB813804, RB952517, RB952522, RB952571, RB952675, RB952681, RB952692, RB952749, RB952754, RB952796, RB952826, and SP7910-11. In micro-region SC, measurements were carried out using the genotypes RB813804, RB952517, RB952522, RB952692, RB952754, RB952796, SP79-1011, RB952511, RB952904, RB953002, RB953133, RB953265, and SP78-4764. In NW the genotypes used were RB72454, RB763710, RB813804, RB952517, CP851491, RB942898, RB952675, RB952514, RB952609, RB952875, RB952884, RB953155, and SP79-1011. In SW the genotypes used were RB952692, RB952597, RB952517, RB952522, RB952571, RB952675, RB952681, RB952609, RB952749, RB952754, SP79-1011, RB952511, SP78-4764, and SP80-1816. In micro-region MB, three measurements were carried out (plant crop and 1st and 2nd ratoon crops) using the genotypes IAC851491, RB763710, RB813804, RB892575, CP851491, RB953214, RB953270, RB953245, SP79-1011, RB953265, RB943339, and RB953206.
Genetic and statistical analyses
Multivariate ANOVA was carried out using the statistical model proposed by Cruz (2006a): Yijk = µ + (b/c)jk + gi + ck + gcik + εijk
where Yijk = ith genotype evaluated in the jth block in the kth harvest
µ = overall mean value
(b/c)jk = effect of block j on the k harvest
gi = effect of treatment (genotype) i
ck = effect of harvest k
gcik = effect of the interaction genotype i and harvest k
εijk = random error associated with the ijk value
So, means (μ) and genotype (g) were considered effect variables, while block (b), harvest (c), interaction genotype-harvest (gc), and experimental error (ε) were considered the random variables.
The results of the multivariate ANOVA of experiments were obtained using the formula shown in Table 1.
|Harvests (H)||c – 1||MSH/MSB|
|Genotypes (G)||g – 1||MSG/MSGH|
|Interaction (G x H)||(c – 1)(g –1)||MSGH/MSR|
|Residual||(g-1)(r – 1)c|
|ℓ = g/(g-1)|
Table 1. Expected mean squares to obtain the results of the multivariate ANOVA of experiments conducted in the sugarcane belt in PE, Brazil, considering the interaction genotype-harvest (sugarcane harvest cycles).
Means were grouped using the Scott and Knott (1974) test at 5% probability. The genetic parameters were estimated according to Cruz (2006a):
Genetic variance component:
Genotype-harvest interaction variance component:
Genetic variance coefficient:
Hartley's F max was calculated before the multivariate ANOVA to evaluate the homogeneity of residuals and, according to Gomes (1990), all variables with a highest-to-lowest error variance ratio above 7 were adjusted for degrees of freedom of the mean residual error (n) and for the interaction between factors (n’) of the multivariate analysis, as described by Cochran (1974).
Yij = μ + gi + aj + εij
Yij = Estimated r for the ith genotype in the jth harvest (time or space)
μ = Overall mean
gi = Effect of the ith genotype under the influence of last harvest (i = 1, 2, …p)
aj = Effect of the temporary harvest on the jth measurement (j = 1, 2, … η)
εij= experimental error based on the temporary effects on harvest in the jth measurement of the ith genotype. The ANOVA arrangement is shown in Table 2.
Table 2. Expected mean squares to obtain the results of estimated r values of experiments based on a two-way ANOVA in experiments conducted in the sugarcane belt in the state of Pernambuco, Brazil.
After the ANOVA, r was obtained using the formula:
Therefore, the number of repeats required to predict the real value of individuals was obtained using the formula:
Next, r was estimated using PCA based on a phenotypic variance and multivariate matrix according to the statistical model introduced by Cruz and Regazzi (2001):
Yij = + gi + aj + εij
However, here the covariance matrix was obtained as:
The main eigenvalue was obtained with:
The corresponding eigenvalue was obtained using the formula:
The r estimator was obtained using:
In turn, the SA used to estimate r was based on the parametric matrix of correlations between genotypes in each analysis pair, and was calculated using:
This estimator of r is the arithmetic mean of the phenotypical correlations between genotypes considering each pair of measurements.
The determination coefficient was obtained using the formula:
While the number of measurements required to predict the real individual value was:
All genetic and statistical analyses were carried out using the Genes (Cruz, 2006b) software.
Results and Discussion
The multivariate ANOVA of experiments revealed significant differences in TPH, TCH, FIB, BRIX, and ATR t. ha-1 across all genotypes grown in the five sugarcane micro-regions in PE, Brazil (Table 3).
|Environments||Traits||Genotype||Harvest||G x H||Residual||Mean||CV(%)||H|
NC: North Coast, SC: South Coast, NW: North Woodlands, SW: South Woodlands, MB: Mid Belt.
G × H: Interaction genotype – harvest, CV (%): Coefficient of experimental variation, H: Hartley's F max test
Table 3. Multivariate ANOVA in experiments conducted in the sugarcane belt in the state of Pernambuco, Brazil, considering the parameters polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable sugars (ATR); and metric tons of ATR per hectare (ATR t. ha-1).
This indicates the existence of high genetic variability to be exploited between the genotypes analyzed concerning these traits, which are considered the most important factors in sugarcane production (Bastos et al., 2003). For PCC, statistically significant differences were observed in all regions, except SC, indicating that it is not possible to select genotypes with high sugar levels in this micro-region. The same finding was observed for the parameter ATR.
For the harvest cycles considered (plant crop and the 1st, 2nd, and 3rd ratoon crops, Table 3), significant differences were observed for all traits, except FIB in MB. For Rosse et al. (2002), these results indicate that harvest cycles represent environmental contrasts due to climate factors and therefore affect the traits evaluated in the present study. Silva (2008) investigated the interaction genotype-environment in sugarcane plantations in the state of São Paulo (SP), Brazil, and underline the fact that the higher mean square values based on site, as compared with other mean squares, signal the existence of considerable differences in yield potential of genotypes in the locations investigated. In the present study a large difference was observed in mean square values for harvest cycles in comparison with mean squares for genotype and the interaction genotype-harvest. Therefore, it may be said that genotypes differ as to yield potential in different harvest cycles, pointing to the need for continuous improvement in long-lived elite genetic material, that is, traits should remain essentially constant throughout harvest cycles.
The coefficients of variation varied from low to medium, according to the classification system proposed by Gomes (1990). Interestingly, in micro-region NW coefficient values were low (< 10) for all traits (Table 3).
The Scott and Knott test (Table 4) shows the emergence of elite genotypes throughout the sugarcane belt in PE, Brazil. In micro-region NC, for instance, genotypes SP79-1011, RB952692, and RB952675 stood out in terms of TPH and TCH. Genotypes RB813804, RB952749, RB952796, RB952754, RB952826, and RB952681 had interesting FIB values. Concerning PCC, except for RB952517 and RB952754, all other genotypes were classified into group ‘a’, with high potential for improvement in this richness trait. This was also valid for BRIX and ATR. For ATR t. ha-1, genotypes SP79-1011, RB952692, RB952675, and RB813804 presented the best performance, forming an elite group.
|Genotypes||TPH t.ha-1||TCH t.ha-1||FIB%||PCC%||BRIX%||ATR kg/t||ATR t.ha-1|
NC: North Coast, SC: South Coast, NW: North Woodlands, SW: South Woodlands, MB: Mid Belt. Means followed by the same letter in a column indicate that values belong to the same group in the Scott and Knott test at 5% probability.
Table 4. Groupings of mean values of the traits polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable sugars (ATR); and metric tons of ATR per hectare (ATR t. ha-1).
In micro-region SC, genotypes RB813804, SP78-4764, and RB962675 exhibited good TPH and TCH values. Concerning the traits FIB, PCC, BRIX, and ATR, genotypes exhibited essentially equivalent potential and were not classified into distinct groups. As for ATR t. ha-1, RB813804, SP7847-64, RB952522, RB952692, SP-79-1011, RB952511, RB953265, and RB952754 exhibited high potential, and were classified into group ‘a’.
In micro-region NW, genotypes RB952875, RB952675, and RB813804 presented the best performance, falling in group ‘a’ in the parameters TCH and ATR t. ha-1. Therefore, improvement programs in NW should be focused on the commercial value of genotype RB952675, since the TPH value was significant. Similarly, important findings in terms of energy generation were observed for genotypes RB952884, RB952609, RB952514, and CP85-1491, which were classified into group ‘a’ concerning the trait FIB, with values over 15%. These values indicate that these genotypes should be exploited specifically in the generation of bioenergy. However, for BRIX, genotypes were not included in any specific group, showing that they present the same potential.
In micro-region SW, genotypes SP78-4764, SP80-1816, RB952511, and SP79-1011 presented excellent performance in TCH and ATR t. ha-1, and were classified into group ‘a’. Similarly, SP78-4764 and RB952511 had high TPH value and were also included in the same group. Genotype RB952609 stood out due to the very high level of FIB, which rendered it the best genotype to be included in improvement programs addressing biomass and bioenergy production. Concerning the parameter PCC, no elite genotype groups were formed. BRIX values of genotypes SP78-4764, SP80-1816, SP79-1011, and RB952681 were high and prompted inclusion in group ‘a’. For ATR, SP79-1011, SP80-1816, and RB952681 also demonstrated superiority and were classified into group ‘a’.
In micro-region MB, RB813804, RB763710, RB892575, RB953265, and SP7910-11 were included in group ‘a’, with high performance for TCH and ATR t. ha-1. Genotypes RB763710 and RB892575 were included in group ‘a’ due to the good TPH values. Concerning FIB levels, genotype RB953245 exhibited the best performance. As for PCC and BRIX, genotypes RB813804, RB953265, SP79-1011, RB943365, and RB953206 were included in group ‘a’, and were superior to the other genotypes. In turn, genotypes RB813804, RB953265, SP79-1011, RB943365, and RB953206 presented the best performance and were included in group ‘a’.
Of the genetic parameters evaluated (Table 5), TPH, TCH, and ATR t. ha-1 presented high genetic variance, superior to the variance of the interaction genotype-harvest in all micro-regions evaluated.
|Environments||Traits||H2||CVg||CVg / CVe|
NC: North Coast, SC: South Coast, NW: North Woodlands, SW: South Woodlands, MB: Mid Belt.
Table 5. Genetic parameters associated with the traits polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable
Genetic variance component.
Interaction genotype-harvest variance component.
h2: Genotypic determination at mean level.
CVg: Genetic variance coefficient.
CVg / CVe: b index.
Such a result is highly desirable, since, according to Dutra Filho (2008a), it indicates that the expression of this important production aspect is mostly due to the genetic effects, suggesting the feasibility to select elite genotypes in the sugarcane belt in PE, Brazil. Concerning ATR t.ha-1, it should be emphasized that the results obtained are highly significant, pointing to the possibility to increase sugarcane productivity significantly in PE and, therefore, to secure better return on investment, since sugarcane prices are defined based on ATR t.ha-1 in Brazil.
For the traits PCC, BRIX, and ATR t.ha-1, selecting elite genetic material is possible only in micro-regions NC, NW, and SW, where genetic variance was higher than the variance in the interaction genotype-harvest. Concerning FIB, only micro-region MB affords to select elite genotypes (Table 5).
The mean heritability coefficient values were high (> 75%) for TPH, THC, and ATR t.ha-1 in all micro-regions (Table 5). For Falconer (1987), such a scenario is highly favorable in the selection of elite genotypes based on these traits, since it indicates the robustness of phenotypical value as an indicator of genotypic value. These results confirm previous findings by Dutra Filho et al. (2008b), in a study that evaluated the progeny of sugarcane in NL at the early improvement stage, and endorse the results published by Melo et al. (2006), who assessed sugarcane genotypes from NW. We also observed that the mean heritability coefficient for TPH was higher than that obtained by Souza et al. (2012) in MB. According to Melo et al. (2009), these interesting heritability coefficients translate as an important aspect in genetic improvement programs, since they clearly show that it is possible to obtain significant genetic advantages in the selection of elite genotypes considering these traits. Since heritability is defined as the transferrable part of the total genetic heritability, Gonçalves et al. (2007) highlight the fact that it is possible to improve production using genotypes with higher performance in these traits obtained crossing the appropriate genetic materials. Mean heritability coefficients were also high for the traits FIB, PCC, BRIX, and ATR t.ha-1 considering the edaphic and climatic conditions in NC, NW, and SW in PE, Brazil. This indicates that elite genotypes should be exploited based on these traits in these micro-regions.
The selection of elite genotypes in the sugarcane belt in PE should be based on TPH, TCH, and ATR t.ha-1. These traits presented the highest genetic variability, with genetic variation coefficients above 10% in all micro-regions (Table 5). According to Oliveira et al. (2008), genetic variation coefficients above 10% are considered high. The ratio of the genetic variation coefficient to the experimental variation coefficient (index b) were above 1 for the traits TPH, TCH, and ATR t.ha-1 (Table 5), showing that these traits are essential in the expression of genetic variability between the genotypes evaluated in the sugarcane belt in PE, Brazil, and lending strength to the hypothesis that selection may become more effective.
Also, r (Table 6) was higher than 0.5 for TPH, TCH, and ATR t.ha-1 across the sugarcane belt in PE, in the three methodologies used. However, in NC and SC, values of r were above 0.75. In NW, r values were considered excellent using the three methodologies, with values over 0.85. For Santos (2004), these results indicate superior genetic control in the expression of these traits in the genotypes evaluated. Consequently, these traits are more evenly expressed across the harvest cycles of sugarcane in the region, which increases ratoon longevity. Cruz and Regazzi (2011) also underline the fact that the higher the r value, the lower the number of measurements required to predict the actual value of genotypes.
NC: North Coast, SC: South Coast, NW: North Woodlands, SW: South Woodlands, MB: Mid Belt.
r: Repeatability coefficient
R²: Determination coefficient
Table 6. Repeatability coefficient (r) values for the traits polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable sugars (ATR); and metric tons of ATR per hectare (ATR t.ha-1) of sugarcane genotypes evaluated in experiments conducted in the sugarcane belt in the state of Pernambuco (PE), Brazil.
Similarly, the values of the traits FIB, PCC, BRIX, and ATR t.ha-1 were regularly repeated in the genotypes evaluated in NC and NW, PE using the three methodologies. In SW, FIB and BRIX also exhibited consistent repeatability, in the three methods. However, in MB only the traits PCC, BRIX, and ATR t.ha-1 were regularly repeated, in PCA.
The estimated number of repeated measurements required to select genotypes with 80%, 90%, and 95% predictability of actual values are shown in Table 7.
NC: North Coast, SC: South Coast, NW: North Woodlands, SW: South Woodlands, MB: Mid Belt.
ANOVA: Analysis of variance; PCA: Principal component analysis; SA: Structured data analyses
Table 7. Number of measurements required to select superior sugarcane genotypes considering the traits polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable sugars (ATR); and metric tons of ATR per hectare (ATR t.ha-1) of sugarcane genotypes evaluated in experiments conducted in the sugarcane belt in the state of Pernambuco (PE), Brazil.
Resende (2002) concluded that determination coefficients above 80% afford appropriate predictive power for the actual value of an individual and, consequently, the number of measurements required for selection. Therefore, it becomes clear that the traits TPH, TCH, and ATR t. ha-1 require no more than two evaluations to select elite genotypes intended for commercial sugarcane plantations in the sugarcane belt in PE, with predictability of 90% of actual values for each trait in NC and NW, using the three methodologies. Similarly, predictability above 80% was observed for SC and SW using the three methodologies, and for MB using PCA and SA.
Concerning the trait FIB, two assessments are enough to select elite genotypes with 80% predictability of actual values for growth in micro-regions NC, NW, and SW, using the three methodologies. In micro-regions SC and MB at least four evaluations, or four harvests, are required to select genotypes with the same 80% predictability level, which is economically unfeasible (Table 7).
For PCC, two evaluations are enough to select genotypes with 80% predictability of actual values using PCA and SA in NC and NW, respectively. For micro-region MB, two evaluations are required to select genotypes with 80% predictability of actual values in a PCA. For micro-region SC, more than four evaluations are necessary (Table 7).
In turn, for the trait BRIX, two evaluations suffice to select genotypes with 80% predictability of actual values in micro-regions NC and SW. For micro-regions NW and MB, two and three evaluations are required, respectively, to select genotypes with 80% predictability, using PCA (Table 7).
For ATR t.ha-1, three evaluations are necessary to select genotypes with 80% predictability of actual values in micro-region NC, in the three methods used, while for micro-region NW two evaluations are required using the three methods. For SW, three evaluations should be carried out, using PCA. For MB, two assessments are necessary using the same method (Table 7).
In estimating the number of measurements, the PCA was more effective than ANOVA and SA to select elite genotypes considering PCC in micro-regions SW and MB. The same was observed for the trait BRIX in micro-regions NW and MB, and for ATR in SW and MB. For Cruz and Regazzi (2001), the experimental error component is more effectively eliminated in the PCA than in the ANOVA. Despite that, the number of assessments required to select superior genotypes with 80% predictability of the actual value in SC, in the three methodologies used, is much higher than four.
In this sense, it is important to evaluate the genotypic stabilization of materials evaluated under the edaphic and climatic conditions in micro-region SC (Table 8).
|Trait||Correlation / No of assessments||ANOVA||R2||CP||R2|
Table 8. Genotypic stabilization of sugarcane genotypes considering the traits polarizable sugars (POL) per hectare, in metric tons (TPH); culm productivity per hectare (TCH); fiber content (FIB); adjusted percent POL (PCC); soluble solids (BRIX); total recoverable sugars (ATR); and metric tons of ATR per hectare (ATR t.ha-1) evaluated in experiments conducted in the sugarcane belt in the state of Pernambuco (PE), Brazil.
Genotypic stabilization is obtained based on the correlations between all measurements carried out. Concerning FIB, for instance, the correlation between the values obtained in plant crop and the 1st and 2nd harvests was negative in ANOVA and very low in PCA. This explains the low R² value, which was zero using ANOVA and only 13.69 using PCA. This explains the need for a significantly higher number of evaluations required to predict the actual value of genotypes. The poor correlation between the measurements the present study indicates the severe climatic variations between harvest years, which negatively influence the regularity of expression of these traits in micro-region SC. However, this was not observed for TPH, when the high correlations between measurements both in ANOVA and in PCA were high, increasing R² values and reducing the number of measurements required to predict the actual value of genotypes.
However, as previously demonstrated for the estimation of genetic traits the selection of superior genotypes in the commercial growth of sugarcane in the sugarcane belt in PE, Brazil should be carried out based on the traits TPH, TCH, and ATR t.ha-1. Also, considering r values, the traits were regularly expressed in the genotypes evaluated, and two measurements are enough to select superior genotypes with 80% predictability of actual values. Pedroso et al. (2011) evaluated sugarcane families grown in the state of Minas Gerais, Brazil, and obtained high repeatability values, concluding that the selection of these families based on two harvests and of clones during the 2nd harvest may be an interesting option in improvement programs, increasing the efficiency in the production of cultivars. Therefore, the results of the present study afford to reduce the time-to-market of new cultivars by three years in the sugarcane belt in PE, Brazil.
The use of two evaluations in the beginning of experiments was enough to select superior genotypes in the sugarcane belt in PE, Brazil. It is possible to reduce the time-to-market of new sugarcane cultivars by three harvest years in the sugarcane belt in PE, Brazil, which affords to decrease labour costs significantly. The best cultivars in micro-region NC were SP79-1011, RB952692, RB952675, and RB813804. For micro-region SC, the best cultivars were RB813804, SP78-4764, RB952675, RB952522, RB952692, SP79-1011, RB952511, RB953265, and RB952754. In NW, the best cultivars were RB952675, RB813804, RB952875, and RB952675, while in SW the best cultivars were SP78-4764, SP80-1816, RB952511, and RB952675. For MB, the best cultivars were RB813804, RB763710, RB892975, RB953265, and SP79-1011.
Conflicts of interest
The authors declare that they have no conflicts of interest.
CAPES - BRAZIL by grant of scholarship.
RIDESA - Interuniversity Network for the Development of the Sucroenergetic Sector for the granting of genotypes for research.
About the Authors
J.A. Dutra Filho
Unidade Academic de Ciências Agrárias, Brazil
- [email protected]
- Abeywardena V (1972). An application of principal component analysis in genetics. J. Gen. 16(01): 27-51. https://doi.org/10.5772/38049
- Barbosa MHP, Silveira LCI (2012). Melhoramento Genético e Recomendação de cultivares. In: Cana-de-açúcar Bioenergia, Açúcar e Etanol – Tecnologias e Perspectivas (Santos F, Borém A and Caldas C, eds.). UFV, Viçosa, 313-353.
- Barbosa MHP, Resende MD, Silveira LCI, Peternelli AL (2005). Estratégias de melhoramento genético da cana-de-açúcar em universidades. In: IX Simpósio sobre seleção recorrente. UFLA, Lavras.
- Bastos IT, Barbosa MHP, Cruz CD, Burniquist WL, et al. (2003). Análise dialélica em clones de cana-de-açúcar. Bragantia 62(02): 199-206.
- Biosev, A. Louis Dreyfus Commodities Company. Available at [http://www.biosev.com]. Accessed April 4, 2013.
- Cochran WG (1974). The combination of Estimates form Different Experiments. Biometrics 10(01): 101-129. https://doi.org/10.2307/3001666
- Coelho COC, Carrazoni PG, Monteiro VLC, Melo FAD, et al. (2001). Biopolímero produzido a partir da cana-de-áçucar para cicatrização cutânea. Acta Cir. Bras. 17(8): 11-13. https://doi.org/10.1590/s0102-86502002000700003
- CONAB, Companhia Nacional de Abastecimento. Available at [http://www.conab.gov.br]. Accessed April 4, 2013.
- Cruz CD (2006a). Programa Genes: Estatística experimental e matrizes. UFV, Viçosa, Brazil.
- Cruz CD (2006b). Programa Genes: Biometria. UFV, Viçosa, Brazil.
- Cruz CD and Regazzi A (2001). Modelos biométricos aplicados ao melhoramento genético .ed. rev. UFV, Viçosa, Brazil.
- Dos Santos MSM, Madalena JA, Soares L, Ferreira PV, et al. (2004). Repetibilidade de características agroindustriais em cana-de-açúcar Pesq. Agropec. Bras. 39(04): 301-306. https://doi.org/10.1590/s0100-204x2004000400001
- Dutra Filho JA, Bastos GQ, Resende LV, Simões Neto DE, et al. (2011b). Avaliação agroindustrial e dissimilaridade genética em progênies e variedades RB de cana-de-açúcar. Agropec. Téc. 32(1): 55-61.
- Dutra Filho JA, T Melo LJO, Resende LV, Anunciação Filho CJ, et al. (2011a). Aplicação de técnicas multivariadas no estudo da divergência genética em cana-de-açúcar. Rev. Ciênc. Agron. 42(1): 185-192. https://doi.org/10.1590/s1806-66902011000100023
- Falconer DS (1987). Introdução a genética quantitativa. UFV, Viçosa, Brazil.
- Fernandes AC (2003). Cálculos na agroindustria da cana-de-açúcar. 2. ed. EME, Piracicaba, Brazil.
- Ferreira A, Barbosa MHP, Cruz CD, Hoffmann HP, et al. (2005). Repetibilidade e número de colheitas para a seleção de clones de cana-de-açúcar. Pesq. Agropec. Bras. 40(8): 761-767. https://doi.org/10.1590/s0100-204x2005000800005
- Gomes PF (1990). Curso de estatística experimental, 13. ed. USP, Piracicaba.
- Gonçalves GM, Viana AP, Bezeera Neto FV, Pereira MG, et al. (2007). Seleção e herdabilidade na produção de ganhos genéticos em maracujá amarelo. Pesq. Agropec. Bras. 42(02): 193-198. https://doi.org/10.1590/s0100-204x2007000200007
- Koffler NF, Lima JFWF, De Lacerda MF, Santana JF, et al. (1986). Caracterização edafo-climática das regiões canavieiras do Brasil: PERNAMBUCO. IAA/PLANALSUCAR, Piracicaba, Brazil.
- Mansour H, Nordheim EV, Ruledge JJ (1981). Estimators of repeatability, Theor. Appl. Gen. 60(03): 151-156. https://doi.org/10.1007/bf00264520
- Melo LJOT, Oliveira FJ, Bastos GQ, Anunciação Filho CJ, et al. (2009). Desempenho agroindustrial de cultivares de cana-de-açúcar na zona da mata litoral sul de Pernambuco. Ciênc. Agrotecnol. 33(03): 684-691. https://doi.org/10.1590/s1413-70542009000300004
- Melo LJTO, Oliveira FJ, Bastos GQ, Anunciação Filho CJ, et al. (2006). Interação genótipo × ciclos de colheita de cana-de-açúcar da Zona da mata Norte de Pernambuco. Bragantia 65(02): 197-205. https://doi.org/10.1590/s0006-87052006000200002
- Oliveira RA, Daros E, Bespalhok-Filho JC, Zambon JLC, et al. (2008). Seleção de famílias de cana-de-açúcar via modelos mistos. Sci. Agrária 09(03): 269-274. https://doi.org/10.5380/rsa.v9i3.11564
- Pedroso CA, Barbosa MHP, Silva FL, Resende MDV, et al. (2011). Repeatability of full-sib sugarcane families across harvests and the efficiency of early selection. Euphytica 182(3): 423-430. https://doi.org/10.1007/s10681-011-0521-z
- Resende MDV (2002). Genética biométrica e estatística no melhoramento de plantas perenes. Embrapa Informação Tecnológica, Brasília, Brazil.
- Rosse LN, Vencovsky R, Ferreira DF. (2002). Comparação de métodos de regressão para avaliar a estabilidade fenotípica em cana-de-açúcar. Pesq. Agropec. Bras. 37(1): 25-32. https://doi.org/10.1590/s0100-204x2002000100004
- Scott AJ and Knott M (1974). A cluster analysis method for grouping means in the analysis of variance. Biometrics 30(3): 507-512. https://doi.org/10.2307/2529204
- Silva MA (2008). Interação genótipo × ambiente e estabilidade fenotípica de cana-de-açúcar em ciclo de cana de ano. Bragantia 67(1): 109-117. https://doi.org/10.1590/s0006-87052008000100013
- Souza PHN, Bastos GQ, Anunciação Filho CG, Dutra Filho JA, et al. (2012). Avaliação de genótipos de cana-de-açúcar para início de safra na Microrregião Centro de Pernambuco. Rev. Ceres 59(5): 427-432. https://doi.org/10.1590/s0034-737x2012000500013
- Stolf R (1986). Metodologia de avaliação de falhas nas linhas de cana-de-açúcar. STAB, Piracicaba 4: 22-36.
- Única, União da Indústria da Cana-de-Açúcar. Available at [http://www.unica.com.br]. Accessed April 4, 2013.
- Share This