Abstract:
The objective of this work was to identify core locations in two soybean macroregions in Brazil for the evaluation and selection of soybean (Glycine max) lineages. Twenty-two cultivars were tested in 23 locations in these two macroregions (MR1 and MR2), during four years, from 2012 to 2015. Trials were conducted in a randomized complete block design with three replicates. All analyses were performed using the GGEbiplot software. The genotype main effects plus genotype × location interaction + genotype main effects plus genotype × environment interaction and genotypic effect vs. genotype x environment effect analyses were used to identify core locations, i.e., locations with high representativeness and consistency of results. Chapada, in the state of Rio Grande do Sul, and Maracaju, in the state of Mato Grosso do Sul, were the core locations in MR1 and MR2, respectively. These locations were the most representative and consistent over the years, and the genotypic effect explained a high proportion of phenotypic variance.
Index terms:
Glycine max; GGL + GGE; soybean macroregion; test locations
Resumo:
O objetivo deste trabalho foi identificar locais-chave de teste em duas macrorregiões sojícolas no Brasil, para avaliação e seleção de linhagens de soja (Glycine max). Foram testadas 22 cultivares em 23 locais, nessas duas macrorregiões (MR1 e MR2), durante quatro anos, de 2012 a 2015. Os ensaios foram conduzidos em delineamento de blocos ao acaso, com três repetições. Todas as análises foram realizadas com uso do programa GGEbiplot. As análises efeito principal de genótipo mais interação genótipo × local + efeito principal de genótipo mais interação genótipo × ambiente e efeito genotípico vs. efeito do genótipo x ambiente foram utilizadas para identificar locais de teste, isto é, locais com alta representatividade e consistência de resultados. Chapada, no Rio Grande do Sul, e Maracaju, no Mato Grosso do Sul, foram os principais locais em MR1 e MR2, respectivamente. Estes locais foram os mais representativos e consistentes ao longo dos anos, e o efeito genotípico explicou uma alta proporção da variância fenotípica.
Termos para indexação:
Glycine max; GGL + GGE; macrorregião sojícola; locais de teste
Introduction
In plant breeding programs, it is essential to conduct trials at locations that are representative of the target environment. This representativeness must be consistent over the years and present a sufficiently high genetic correlation with the target environment (Dia et al., 2016DIA, M.; WEHNER, T.C.; HASSELL, R.; PRICE, D.S.; BOYHAN, G.E.; OLSON, S.; KING, S.; DAVIS, A.R.; TOLLA, G.E.; BERNIER, J.; JUAREZ, B. Value of locations for representing mega-environments and for discriminating yield of watermelon in the U.S. Crop Science, v.56, p.1726-1735, 2016. DOI: https://doi.org/10.2135/cropsci2015.11.0698.
https://doi.org/10.2135/cropsci2015.11.0...
; Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
). Furthermore, locations and selection must be efficient in differentiating superior genotypes (Qin et al., 2015QIN, J.; XU, R.; LI, H.; YANG, C.; LIU, D.; LIU, Z.; ZHANG, L.; LU, W.; FRETT, T.; CHEN, P.; ZHANG, M.; QIU, L. Evaluation of productivity and stability of elite summer soybean cultivars in multi-environment trials. Euphytica, v.206, p.759-773, 2015. DOI: https://doi.org/10.1007/s10681-015-1513-1.
https://doi.org/10.1007/s10681-015-1513-...
; Krishnamurthy et al., 2017KRISHNAMURTHY, S.L.; SHARMA, P.C.; SHARMA, D.K.; RAVIKIRAN, K.T.; SINGH, Y.P.; MISHRA, V.K.; BURMAN, D.; MAJI, B.; MANDAL, S.; SARANGI, S.K.;GAUTAM, R.K.; SINGH, P.K.; MANOHARA, K.K.; MARANDI, B.C.; PADMAVATHI, G.; VANVE, P.B.; PATIL, K.D.; THIRUMENI, S.; VERMA, O.P.; KHAN, A.H.; TIWARI, S.; GEETHA, S.; SHAKILA, M.; GILL, R.; YADAV, V.K.; ROY, S.K.B.; PRAKASH, M.; BONIFACIO, J.; ISMAIL, A.; GREGORIO, G.B.; SINGH, R.K. Identification of mega-environments and rice genotypes for general and specific adaptation to saline and alkaline stresses in India. Scientific Reports, v.7, art.7968, 2017. DOI: https://doi.org/10.1038/s41598-017-08532-7.
https://doi.org/10.1038/s41598-017-08532...
). When a location combines both representativeness and the ability to discriminate genotypes, it is called a core location (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). The identification of a core location is fundamental in generations with low seed availability, such as in segregating populations and recently selected lines, i.e., in progenies from a single plant and in preliminary trials (Wu et al., 2013WU, J.; BONDALAPATI, K.; GLOVER, K.; BERZONSKY, W.; JENKINS, J.N.; MCCARTY, J.C. Genetic analysis without replications: model evaluation and application in spring wheat. Euphytica, v.190, p.447-458, 2013. DOI: https://doi.org/10.1007/s10681-012-0835-5.
https://doi.org/10.1007/s10681-012-0835-...
). Some of the difficulties faced by soybean breeders when choosing the best location for the selection of superior genotypes include both limited seed availability in early generations (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.) and resources (Mi et al., 2011MI, X.; WEGENAST, T.; UTZ, H.F.; DHILLON, B.S.; MELCHINGER, A.E. Best linear unbiased prediction and optimum allocation of test resources in maize breeding with doubled haploids. Theoretical and Applied Genetics, v.123, p.1-10, 2011. DOI: https://doi.org/10.1007/s00122-011-1561-4.
https://doi.org/10.1007/s00122-011-1561-...
); because of the latter, trials are usually conducted in a few locations - usually in only one -, which has direct implications on the effectiveness of the breeding program.
Macroregions for soybean adaptation are widely used to identify core locations for the selection and evaluation of new lineages (Kaster & Farias, 2012KASTER, M.; FARIAS, J.R.B. Regionalização dos testes de valor de cultivo e uso e da indicação de cultivares de soja: terceira aproximação. Londrina: Embrapa Soja, 2012. Available at: <Available at: http://ainfo.cnptia.embrapa.br/digital/bitstream/item/54939/1/Doc-330-OL1.pdf
>. Accessed on: Feb. 14 2018.
http://ainfo.cnptia.embrapa.br/digital/b...
). Within these regions, breeders identify representative locations to conduct their trials (Yan et al., 2010YAN, W.; FRÉGEAU-REID, J.; PAGEAU, D.; MARTIN, R.; MITCHELL-FETCH, J.; ETIENNE, M.; ROWSELL, J.; SCOTT, P.; PRICE, M.; HAAN, B. de; CUMMISKEY, A. LAJEUNESSE, J.; DURAND, J.; SPARRY, E. Identifying essential test locations for oat breeding in eastern Canada. Crop Science, v.50, p.504-515, 2010. DOI: https://doi.org/10.2135/cropsci2009.03.0133.
https://doi.org/10.2135/cropsci2009.03.0...
). It should be noted that the task of identifying a single location is an important factor (Das et al., 2019DAS, A.; PARIHAR, A.K.; SAXENA, D.; SINGH, D.; SINGHA, K.D.; KUSHWAHA, K.P.S.; CHAND, R.; BAL, R.S.; CHANDRA S.; GUPTA, S. Deciphering genotype-by-environment interaction for targeting test environments and rust resistant genotypes in field pea (Pisum sativum L.). Frontiers in Plant Science, v.10, art.825, 2019. DOI: https://doi.org/10.3389/fpls.2019.00825.
https://doi.org/10.3389/fpls.2019.00825...
), requiring a theoretical and practical foundation to avoid mistakes in selection. Misconceptions in selecting the best lineages in the preliminary testing phase may result in the elimination of highly desirable lineages and in a reduction in the efficiency of the breeding program. Up to date, no known work has been performed to identify the key test locations for soybean breeding in the South of Brazil.
The objective of this work was to identify core locations in two soybean macroregions in Brazil for the evaluation and selection of soybean lineages.
Materials and Methods
Data from soybean [Glycine max (L.) Merr.] value for cultivation and use (VCU) trials carried out in the 2012-2013, 2013-2014, 2014-2015, and 2015-2016 crop seasons were used. Trials were conducted in 23 locations, 11 in the first and 12 in the second evaluated soybean macroregion (MR1 and MR2, respectively) (Figure 1 and Table 1). MR1, located in the South of Brazil, is characterized by a higher altitude and a mesothermal climate without dry seasons, classified as Cfa or Cfb, according to Köppen, with temperatures in the coldest month ranging from -3 to 18ºC (Alvares et al., 2013ALVARES, C.A.; STAPE, J.L.; SENTELHAS, P.C.; GONÇALVES, J.L. de M.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. DOI: https://doi.org/10.1127/0941-2948/2013/0507.
https://doi.org/10.1127/0941-2948/2013/0...
); this macroregion covers the states of Rio Grande do Sul and Santa Catarina, as well as the mid-southern and southeastern regions of the state of Paraná and the south of the state of São Paulo. In contrast, MR2, which consists of western and northern Paraná, southern and western São Paulo, and southern Mato Grosso do Sul, is considered a transitional region between Cfa and Cwa climates (Alvares et al., 2013ALVARES, C.A.; STAPE, J.L.; SENTELHAS, P.C.; GONÇALVES, J.L. de M.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. DOI: https://doi.org/10.1127/0941-2948/2013/0507.
https://doi.org/10.1127/0941-2948/2013/0...
), presenting higher predominant temperatures in summer and dry winters (Kaster & Farias, 2012KASTER, M.; FARIAS, J.R.B. Regionalização dos testes de valor de cultivo e uso e da indicação de cultivares de soja: terceira aproximação. Londrina: Embrapa Soja, 2012. Available at: <Available at: http://ainfo.cnptia.embrapa.br/digital/bitstream/item/54939/1/Doc-330-OL1.pdf
>. Accessed on: Feb. 14 2018.
http://ainfo.cnptia.embrapa.br/digital/b...
).
Brazilian map with trial locations (black dots) in the two evaluated macroregions (MR1 in dark gray and MR2 in light gray) of soybean (Glycine max) adaptation in Southern Brazil. States: MS, Mato Grosso do Sul; MG, Minas Gerais; SP, São Paulo; PR, Paraná; SC, Santa Catarina; and RS, Rio Grande do Sul.
Locations (municipalities) of the value for cultivation and use trials in the two evaluated macroregions (MR1 and MR2) of soybean (Glycine max) adaptation in Brazil.
A total of 22 genotypes were tested, including 8 cultivars and 14 lineages (Table 2). It should be pointed out that, as locations and genotypes were not the same in all studied years, the dataset is highly imbalanced.
Experiments were conducted in a randomized complete block design with three replicates. Plots consisted of four rows with a length of 5 m, and row and plot spacing were of 0.5 m. Sowing density was standardized for all genotypes at 30 seeds per square meter. Basic fertilization consisted of a mineral fertilizer with 7 kg ha-1 N, 70 kg ha-1 P2O5, and 70 kg ha-1 K2O. Harvest was performed using a plot harvester when plants reached harvest maturity at the R8 phenological stage. Both central rows of each plot (5 m2) were harvested, and seed moisture content was routinely adjusted to 13% to obtain grain yield (kg ha-1).
The statistical analysis for the identification of core locations was performed using the GGEbiplot software (Yan, 2001YAN, W. GGEbiplot - a Windows application for graphical analysis of multienvironment trial data and other types of two-way data. Agronomy Journal, v.93, p.1111-1118, 2001. DOI: https://doi.org/10.2134/agronj2001.9351111x.
https://doi.org/10.2134/agronj2001.93511...
). The genotype main effects plus genotype × location interaction (GGL) + genotype main effects plus genotype × environment interaction effect (GGE) analysis was used to identify core locations in each of the two previously defined macroregions. In this analysis, the location on the biplot was defined by the mean of both PC1 and PC2 over the tested years. The data were scaled using a heritability-adjusted environment interaction effect (HA-GGE) biplot, in which data are scaled by standard deviation and adjusted heritability (Yan & Holland, 2010YAN, W.; HOLLAND, J.B. A heritability-adjusted GGE biplot for test environment evaluation. Euphytica, v.171, p.355-369, 2010. DOI: https://doi.org/10.1007/s10681-009-0030-5.
https://doi.org/10.1007/s10681-009-0030-...
). The cosine of the angle between the vector of the location and the average environment axis (AEA) line corresponds to the genetic correlation (rg) between them (Yan & Holland, 2010YAN, W.; HOLLAND, J.B. A heritability-adjusted GGE biplot for test environment evaluation. Euphytica, v.171, p.355-369, 2010. DOI: https://doi.org/10.1007/s10681-009-0030-5.
https://doi.org/10.1007/s10681-009-0030-...
; Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). Therefore, the smaller the angle between location and average environment, the more representative is the test location. Vector length indicates the consistency of the results over the years, i.e., the representativeness of the location; moreover, when the biplot explanation is high, vector length is proportional to the squared root of heritability (h) (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.).
A linear map was built to facilitate the identification of similarity among environments. This way, it was possible to identify the patterns of the genotype × location interaction and genotype x environment (GE) interaction. GE was defined as follows: GE = GL + GY + GLY, where GL is the genotype × location interaction, GY is the genotype × year interaction, and GLY is the genotype × location × year interaction (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
). Therefore, the GE interaction increases when the year effect is high. The genotypic effect (G) vs. GE analysis allows for inferences about the capability of the environment (considered a combination of location + year) for the selection of superior genotypes (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). The vector analysis of the environments, which also allows for the construction of a linear map, shows the distance between environments; the proximity among environments indicates a positive genetic correlation among them (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). In addition, the position of environments on the linear map makes it possible to identify patterns in relation to GL and GE. Therefore, if environments are mainly grouped by location and not year on the map, GL effects will surpass those of GE and vice versa (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.).
For all analyses, the following parameters of the GGEbiplot software (Yan, 2001YAN, W. GGEbiplot - a Windows application for graphical analysis of multienvironment trial data and other types of two-way data. Agronomy Journal, v.93, p.1111-1118, 2001. DOI: https://doi.org/10.2134/agronj2001.9351111x.
https://doi.org/10.2134/agronj2001.93511...
, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.) were used: “Transform = 0”, meaning data were not transformed; “Scaling = 2”, data scaled by standard deviation (SD-scaled) and by adjusted heritability (h-weighted); and “Data centering = 2”, data centering by genotype + genotype x environment interaction (G + GE) and singular value partitioning (“SVP = 2”), focusing on the environment.
Results and Discussion
According to the GGL + GGE analysis, Chapada is the most representative test location because it was highly correlated with the average environment, indicating that the obtained results must be consistent over the years (long vector) (Figure 2). Furthermore, it was observed that this test location discriminates genotypes mainly through G (Yan, 2015YAN, W. Mega-environment analysis and test location evaluation based on unbalanced multiyear data. Crop Science, v.55, p.113-122, 2015. DOI: https://doi.org/10.2135/cropsci2014.03.0203.
https://doi.org/10.2135/cropsci2014.03.0...
).
Biplot for genotype main effects plus genotype × location interaction (GGL) + genotype main effects plus genotype × environment interaction (GGE), using data from the value for cultivation and use trials in the 2012/2013-2015/2016 crop seasons in macroregion 1 of soybean (Glycine max) adaption in Brazil. Genotypes are represented by +. Environments consist of locations (municipalities) in the Southern region of the country: SFA, São Francisco de Assis; GUA, Guarapuava; VAC, Vacaria; ABL, Abelardo Luz; RLZ, Realeza; PSO, Palma Sola; CHA, Chapada; CPN, Campos Novos; STO, Santo Augusto; ERC, Erechim; and STC, Santa Cruz do Sul.
In the G vs. GE analysis based on several environments, the interpretation of the biplot may be limited by location overlapping (Figure 3). However, more intelligible results for MR1 and MR2 are presented in Tables 3 and 4, respectively. The average environment coordinate index on the x-axis (AEC_X) is integrated for environmental evaluations and is useful even in cases of a low explanation of the biplot (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). Therefore, the higher the AEC_X value, the higher is the potential of the environment to select superior genotypes. The AEC y-axis (AEC_Y) column refers to the potential of a location for the selection of genotypes through G or GE. The higher the modular AEC_Y value, the lower the potential of the location to select genotypes by G. Therefore, environments characterized by AEC_Y values between 0.30 and -0.30 were considered adequate for selecting genotypes through G, whereas environments with AEC_Y above these values were considered adequate for selecting genotypes based on GE. In summary, the most desirable environment will be characterized by a high AEC_X value and a close-to-zero AEC_Y value.
Biplot displaying the ability of each environment to select according to genotypic effect (G) vs. genotype x environment effects (GE) (A) and environment vector view of the GGE biplot in a linear map (B), using data from the value for cultivation and use trials in the 2012/2013-2015/2016 crop seasons in macroregion 1 of soybean (Glycine max) adaptation in Brazil. Environments consist of location (Brazilian municipality) combined with different experimental years: STC, Santa Cruz do Sul; GUA, Guarapuava; SFA, São Francisco de Assis; RLZ, Realeza; VAC, Vacaria; CHA, Chapada; ABL, Abelardo Luz; PSO, Palma Sola; STO, Santo Augusto; CPN, Campos Novos; and ERC, Erechim.
Numerical values of the genotypic effect (G) vs. genotype x environment effect (GE) analysis for locations with the ability to select superior genotypes (AEC_X), stable by G and unstable by GE (AEC_Y), as well as vector length and representativeness of the environment correlated with the average environment axis (AEA), for macroregion 1 (MR1) of soybean (Glycine max) adaptation in Brazil.
Numerical values of the genotypic effect (G) vs. genotype x environment effect (GE) analysis for locations with the ability to select superior (AEC_X) genotypes, stable by G and unstable by GE (AEC_Y), as well as vector length and representativeness of the environment correlated with the average environment axis (AEA), for macroregion 2 (MR2) of soybean (Glycine max) adaptation in Brazil.
The GGL + GGE analysis also indicated that Chapada and Palma Sola are the most representative locations in MR1, since both had a high rg with an elevated average environment and a genetic correlation between them (Figure 2). The G vs. GE analysis revealed that Campos Novos, in 2014, presented the highest potential for selection of superior genotypes, i.e., the highest rgh (Figure 3 A and Table 3); however, this environment is less representative of MR1, compared with Chapada and Palma Sola (Figure 2). The Chapada environment was placed to the right of the AEC line across years, which indicates its effectiveness as a site for selecting superior genotypes. This environment showed the best performance among these locations in 2015, but it exceeded optimal vector length in 2014, which indicates that it would not be an effective site for selecting genotypes by G. Also in 2014, Palma Solo had the third highest potential for selecting superior genotypes. The linear correlation analysis (linear mapping) showed that Chapada, in 2012 and 2015, was highly correlated and close on the linear map, indicating its similarity over the years (Figure 3 B). However, most environments were mainly grouped by year and not by location, which indicates that GE dominates over GL.
GE consists mainly in the GL, GY, and GLY interactions (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
). In this way, the identification of locations that have a part of GL that can be predictable, such as Chapada, is important to confirm their efficiency for the selection of superior genotypes. The predictable GL mainly occurs when the location factor is dominated by stable, predictable, or controllable physical factors. The controllable physical factors consist of characteristics as altitude, latitude, soil properties, and established routine managements; the latter includes planting dates, planting density, irrigation, fertilization, and herbicide and fungicide applications (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
).
The São Francisco de Assis and Santa Cruz do Sul locations are the less representative of MR1, with a lower genetic association with the average environment and short vectors in the GGL + GGE analysis (Figure 2). Furthermore, Santa Cruz do Sul, in 2015, and Guarapuava, in 2013, showed negative AEC_X values, which indicate a negligible potential for selecting superior genotypes. Therefore, Chapada, Palma Sola, and Campos Novos are the best candidates for a core location within MR1, but Chapada stands out, indicating that it can be considered the core location for this macroregion. According to Yan et al. (2011)YAN, W.; PAGEAU, D.; FRÉGEAU-REID, J.; DURAND, J. Assessing the representativeness and repeatability of test locations for genotype evaluation. Crop Science, v.51, p.1603-1610, 2011. DOI: https://doi.org/10.2135/cropsci2011.01.0016.
https://doi.org/10.2135/cropsci2011.01.0...
, test locations as Chapada are of type I, with a high representativity (small acute angles with the AEA) and repeatability (small acute angles between years within a location), indicating that a high selection intensity can be applied. Locations such as Palma Sola and Campos Novos, with a low-to-moderate representativity and a high repeatability, represent type II(a), i.e., they are useful for culling unstable genotypes, and a low selection intensity should be adopted during plant selection. Therefore, due to these characteristics, only Chapada is appropriate for plant selection in preliminary trials.
In MR2, the GGL + GGE analysis indicated that Maracaju is the most representative location, i.e., it had the highest rg (Figure 4). In 2013 and 2012, the environment showed the highest potential for identifying superior genotypes, since it was located to the right of the AEA line (higher rgh). In addition, selection in 2013 seems to be mostly dependent on G, because the environmental vector is shorter towards AEC than towards AEA (Figure 5 A and Table 4). The linear map showed that this environment was closely related in 2012, 2013, and 2015, representing a high rg, and that the GL effect partially dominates over that of GE (Figure 5 B). Similarly to Chapada, the predictable GL in Maracaju is important for choosing this core location, where a repeatable GE occurs across years (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
).
Biplot for genotype main effects plus genotype × location interaction (GGL) + genotype main effects plus genotype × environment interaction (GGE), using data from the value for cultivation and use trials in the 2012/2013-2015/2016 crop seasons for macroregion 2 of soybean (Glycine max) adaption in Brazil. Genotypes are represented by +. Environments consist of locations (municipalities) in the Southern region of the country: PLT, Palotina; UBI, Ubiratã; SJI, São Jorge do Ivaí; CCV, Cascavel; PPO, Ponta Porã; LON, Londrina; CMO, Cândido Mota; SID, Sidrolândia; MCJ, Maracaju; DOU, Dourados; PRB, Perobal; BRA, Brasilândia do Sul.
Biplot displaying the ability of the environment to select for genotypic effect (G) vs. genotype x environment effects (GE) (A) and environment vector view of the GGE biplot in a linear map (B), using data from the value for cultivation and use trials in the 2012/2013-2015/2016 crop seasons in microregion 2 of soybean (Glycine max) adaptation in Brazil. Environments consist of location (Brazilian municipality) combined with different experimental years: CCV, Cascavel; PRB, Perobal; SID, Sidrolândia; PLT, Palotina; UBI, Ubiratã; SJI, São Jorge do Ivaí; PPO, Ponta Porã; BRA, Brasilândia do Sul; CMO, Cândido Mota; LON, Londrina; DOU, Dourados; and MCJ, Maracaju.
According to the G vs. GE analysis, the following environments had negative AEC_X values and were located on the left of the AEC line: Cascavel in 2013, Perobal in 2014, Palotina in 2014, Sidrolândia in 2014, Ubiratã in 2014, Ponta Porã in 2014, São Jorge do Ivaí in 2014, and Brasilândia do Sul in 2014 (Figure 5 A). Moreover, Palotina, Ubiratã, Cascavel, São Jorge do Ivaí, Sidrolândia, and Perobal showed short vectors in the GGL + GGE analysis, which indicates inconsistency of results over the years and a low heritability. It was also observed that Brasilândia do Sul and Palotina are scarcely representative environments, and that Palotina, Ubiratã, Cascavel, São Jorge do Ivaí, Sidrolândia, and Perobal present short vectors with inconsistent results throughout the tested years. Therefore, the obtained results are indicative that Maracaju is the most suitable to be designated as a core location in MR2 for evaluating segregating populations. This means that this environment is representative of the macroregion, i.e., has a high rg with the average environment and is characterized by a long vector, which indicates consistent results over the years. In addition, in 2012 and 2013, Maracaju exhibited the highest potential as a site for the efficient selection of superior genotypes. In this way, this location is suitable to conduct trials with no-repeatable evaluation in locations or years, as done for progeny lines, commonly in one or two locations in a single year, with no or few replicates (Yan et al., 2011YAN, W.; PAGEAU, D.; FRÉGEAU-REID, J.; DURAND, J. Assessing the representativeness and repeatability of test locations for genotype evaluation. Crop Science, v.51, p.1603-1610, 2011. DOI: https://doi.org/10.2135/cropsci2011.01.0016.
https://doi.org/10.2135/cropsci2011.01.0...
).
In summary, Chapada and Maracaju were found to be the closest to a core location for MR1 and MR2, respectively. Both are desirable environments in most of the studied years, according to the G vs. GE analysis (Figures 2 and 4). Environments with short vectors towards AEC are desirable for selecting superior genotypes through G (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.), as observed in Chapada in 2015 and 2013 and in Maracaju in 2012 and 2013. However, locations to the left side of the AEC line present negative rgh and are inadequate for selecting superior genotypes (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). In practice, the identification of core locations, such as Chapada and Maracaju, can improve selection accuracy and provide a better allocation of resources. This is important due to the limited availability of resources in plant breeding programs, which requires breeders to choose carefully the correct number of replicates and environments for each genotype (González-Barrios et al., 2019GONZÁLEZ-BARRIOS, P.; DÍAZ-GARCÍA, L.; GUTIÉRREZ, L. Mega-environmental design: using genotype × environment interaction to optimize resources for cultivar testing. Crop Science, v.59, p.1899-1915, 2019. DOI: https://doi.org/10.2135/cropsci2018.11.0692.
https://doi.org/10.2135/cropsci2018.11.0...
). Moreover, the repeatability of the genotype x environment interaction between years, which was verified mainly in Maracaju, is essential to increase the predictive ability of the environment (Lado et al., 2016LADO, B.; GONZÁLES BARRIOS, P.; QUINCKE, M.; SILVA, P.; GUTIÉRREZ, L. Modeling genotype × environment interaction for genomic selection with unbalanced data from a wheat breeding program. Crop Science, v.56, p.2165-2179, 2016. DOI: https://doi.org/10.2135/cropsci2015.04.0207.
https://doi.org/10.2135/cropsci2015.04.0...
).
Another important use of GGE biplots, besides identifying the best test locations, is the possibility of also identifying those that are inappropriate for trials (Krishnamurthy et al., 2017KRISHNAMURTHY, S.L.; SHARMA, P.C.; SHARMA, D.K.; RAVIKIRAN, K.T.; SINGH, Y.P.; MISHRA, V.K.; BURMAN, D.; MAJI, B.; MANDAL, S.; SARANGI, S.K.;GAUTAM, R.K.; SINGH, P.K.; MANOHARA, K.K.; MARANDI, B.C.; PADMAVATHI, G.; VANVE, P.B.; PATIL, K.D.; THIRUMENI, S.; VERMA, O.P.; KHAN, A.H.; TIWARI, S.; GEETHA, S.; SHAKILA, M.; GILL, R.; YADAV, V.K.; ROY, S.K.B.; PRAKASH, M.; BONIFACIO, J.; ISMAIL, A.; GREGORIO, G.B.; SINGH, R.K. Identification of mega-environments and rice genotypes for general and specific adaptation to saline and alkaline stresses in India. Scientific Reports, v.7, art.7968, 2017. DOI: https://doi.org/10.1038/s41598-017-08532-7.
https://doi.org/10.1038/s41598-017-08532...
). In this sense, environments on the left side of the AEC line on the G vs. GE analysis are considered ineffective for selection and, therefore, undesirable, and may be replaced by other environments with a higher potential for selecting superior genotypes (Yan, 2014YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.). This was the case of Santa Cruz do Sul in 2015 and Guarapava in 2013 in MR1, and of Cascavel in 2013, Perobal in 2014, Ponta Porão in 2014, Brasilândia do Sul in 2014, Ubiratã in 2014, and São Jorge do Ivaí in 2014 in MR2 (Figures 3 and 5). Nonrepresentative locations with short vectors (inconsistent results) in the GGL + GGE analysis may also be replaced. São Francisco de Assis and Santa Cruz do Sul, in MR1, and Brasilândia do Sul, Palotina, Ubiratã, Cascavel, São Jorge do Ivaí, Sidrolândia, and Perobal, in MR2, had short vectors and an elevated angle (low rg with the average environment); therefore, they do not represent the target environment by the GGL + GGE analysis. Moreover, all these environments showed inconsistent results across years, indicating that they cannot be designated as core locations (Yan, 2015).
In both macroregions, the environment effects are mainly ranked by year and not by location, leading to a higher correlation due to GE. This occurs because the year effect is naturally randomized and unpredictable; therefore, GY and GLY are not reproducible (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
). Minimizing environmental effects is essential for plant breeding, because environmental control decreases the effect of GE over GL, besides enabling a higher heritability and experimental precision. This was observed for Maracaju, which kept a similar position on the linear map, especially in 2012, 2013, and 2015. A portion of GL may repeat with a decreased year effect, which is related to location effects such as sowing date, planting density, irrigation, soil fertility, and other crop management factors. These characteristics minimize the effects of the interactions involving the year factor (Yan, 2016YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
https://doi.org/1770.2135/cropsci2015.06...
). Therefore, the success of any breeding effort greatly depends on the high representativeness and potential for the selection of superior genotypes of each core location, as well as on the homogeneity of the experimental area and on an efficient crop management practice (Lucio & Sari, 2017LUCIO, A.D.; SARI, B.G. Planning and implementing experiments and analyzing experimental data in vegetable crops: problems and solutions. Horticultura Brasileira, v.35, p.316-327, 2017. DOI: https://doi.org/10.1590/s0102-053620170302.
https://doi.org/10.1590/s0102-0536201703...
).
The obtained results are indicative that the identification of core locations for the evaluation and selection of lineages in preliminary trials, with a large number of lines and a low seed availability, is of great importance for plant breeding.
Conclusion
The municipalities of Chapada, in the state of Rio Grande do Sul, and of Maracaju, in the state of Mato Grosso do Sul, Brazil, are the best core locations for macroregions 1 and 2 of soybean (Glycine max) adaptation in the country, respectively.
References
- ALVARES, C.A.; STAPE, J.L.; SENTELHAS, P.C.; GONÇALVES, J.L. de M.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. DOI: https://doi.org/10.1127/0941-2948/2013/0507.
» https://doi.org/10.1127/0941-2948/2013/0507 - DAS, A.; PARIHAR, A.K.; SAXENA, D.; SINGH, D.; SINGHA, K.D.; KUSHWAHA, K.P.S.; CHAND, R.; BAL, R.S.; CHANDRA S.; GUPTA, S. Deciphering genotype-by-environment interaction for targeting test environments and rust resistant genotypes in field pea (Pisum sativum L.). Frontiers in Plant Science, v.10, art.825, 2019. DOI: https://doi.org/10.3389/fpls.2019.00825.
» https://doi.org/10.3389/fpls.2019.00825 - DIA, M.; WEHNER, T.C.; HASSELL, R.; PRICE, D.S.; BOYHAN, G.E.; OLSON, S.; KING, S.; DAVIS, A.R.; TOLLA, G.E.; BERNIER, J.; JUAREZ, B. Value of locations for representing mega-environments and for discriminating yield of watermelon in the U.S. Crop Science, v.56, p.1726-1735, 2016. DOI: https://doi.org/10.2135/cropsci2015.11.0698.
» https://doi.org/10.2135/cropsci2015.11.0698 - GONZÁLEZ-BARRIOS, P.; DÍAZ-GARCÍA, L.; GUTIÉRREZ, L. Mega-environmental design: using genotype × environment interaction to optimize resources for cultivar testing. Crop Science, v.59, p.1899-1915, 2019. DOI: https://doi.org/10.2135/cropsci2018.11.0692.
» https://doi.org/10.2135/cropsci2018.11.0692 - KASTER, M.; FARIAS, J.R.B. Regionalização dos testes de valor de cultivo e uso e da indicação de cultivares de soja: terceira aproximação. Londrina: Embrapa Soja, 2012. Available at: <Available at: http://ainfo.cnptia.embrapa.br/digital/bitstream/item/54939/1/Doc-330-OL1.pdf >. Accessed on: Feb. 14 2018.
» http://ainfo.cnptia.embrapa.br/digital/bitstream/item/54939/1/Doc-330-OL1.pdf - KRISHNAMURTHY, S.L.; SHARMA, P.C.; SHARMA, D.K.; RAVIKIRAN, K.T.; SINGH, Y.P.; MISHRA, V.K.; BURMAN, D.; MAJI, B.; MANDAL, S.; SARANGI, S.K.;GAUTAM, R.K.; SINGH, P.K.; MANOHARA, K.K.; MARANDI, B.C.; PADMAVATHI, G.; VANVE, P.B.; PATIL, K.D.; THIRUMENI, S.; VERMA, O.P.; KHAN, A.H.; TIWARI, S.; GEETHA, S.; SHAKILA, M.; GILL, R.; YADAV, V.K.; ROY, S.K.B.; PRAKASH, M.; BONIFACIO, J.; ISMAIL, A.; GREGORIO, G.B.; SINGH, R.K. Identification of mega-environments and rice genotypes for general and specific adaptation to saline and alkaline stresses in India. Scientific Reports, v.7, art.7968, 2017. DOI: https://doi.org/10.1038/s41598-017-08532-7.
» https://doi.org/10.1038/s41598-017-08532-7 - LADO, B.; GONZÁLES BARRIOS, P.; QUINCKE, M.; SILVA, P.; GUTIÉRREZ, L. Modeling genotype × environment interaction for genomic selection with unbalanced data from a wheat breeding program. Crop Science, v.56, p.2165-2179, 2016. DOI: https://doi.org/10.2135/cropsci2015.04.0207.
» https://doi.org/10.2135/cropsci2015.04.0207 - LUCIO, A.D.; SARI, B.G. Planning and implementing experiments and analyzing experimental data in vegetable crops: problems and solutions. Horticultura Brasileira, v.35, p.316-327, 2017. DOI: https://doi.org/10.1590/s0102-053620170302.
» https://doi.org/10.1590/s0102-053620170302 - MI, X.; WEGENAST, T.; UTZ, H.F.; DHILLON, B.S.; MELCHINGER, A.E. Best linear unbiased prediction and optimum allocation of test resources in maize breeding with doubled haploids. Theoretical and Applied Genetics, v.123, p.1-10, 2011. DOI: https://doi.org/10.1007/s00122-011-1561-4.
» https://doi.org/10.1007/s00122-011-1561-4 - QIN, J.; XU, R.; LI, H.; YANG, C.; LIU, D.; LIU, Z.; ZHANG, L.; LU, W.; FRETT, T.; CHEN, P.; ZHANG, M.; QIU, L. Evaluation of productivity and stability of elite summer soybean cultivars in multi-environment trials. Euphytica, v.206, p.759-773, 2015. DOI: https://doi.org/10.1007/s10681-015-1513-1.
» https://doi.org/10.1007/s10681-015-1513-1 - WU, J.; BONDALAPATI, K.; GLOVER, K.; BERZONSKY, W.; JENKINS, J.N.; MCCARTY, J.C. Genetic analysis without replications: model evaluation and application in spring wheat. Euphytica, v.190, p.447-458, 2013. DOI: https://doi.org/10.1007/s10681-012-0835-5.
» https://doi.org/10.1007/s10681-012-0835-5 - YAN, W. Analysis and handling of G × E in a practical breeding program. Crop Science, v.56, p.2106-2118, 2016. DOI: https://doi.org/1770.2135/cropsci2015.06.0336.
» https://doi.org/1770.2135/cropsci2015.06.0336 - YAN, W. Crop variety trials: data management and analysis. Chichester: Wiley-Blackwell, 2014. 360p.
- YAN, W. GGEbiplot - a Windows application for graphical analysis of multienvironment trial data and other types of two-way data. Agronomy Journal, v.93, p.1111-1118, 2001. DOI: https://doi.org/10.2134/agronj2001.9351111x.
» https://doi.org/10.2134/agronj2001.9351111x - YAN, W. Mega-environment analysis and test location evaluation based on unbalanced multiyear data. Crop Science, v.55, p.113-122, 2015. DOI: https://doi.org/10.2135/cropsci2014.03.0203.
» https://doi.org/10.2135/cropsci2014.03.0203 - YAN, W.; FRÉGEAU-REID, J.; PAGEAU, D.; MARTIN, R.; MITCHELL-FETCH, J.; ETIENNE, M.; ROWSELL, J.; SCOTT, P.; PRICE, M.; HAAN, B. de; CUMMISKEY, A. LAJEUNESSE, J.; DURAND, J.; SPARRY, E. Identifying essential test locations for oat breeding in eastern Canada. Crop Science, v.50, p.504-515, 2010. DOI: https://doi.org/10.2135/cropsci2009.03.0133.
» https://doi.org/10.2135/cropsci2009.03.0133 - YAN, W.; HOLLAND, J.B. A heritability-adjusted GGE biplot for test environment evaluation. Euphytica, v.171, p.355-369, 2010. DOI: https://doi.org/10.1007/s10681-009-0030-5.
» https://doi.org/10.1007/s10681-009-0030-5 - YAN, W.; PAGEAU, D.; FRÉGEAU-REID, J.; DURAND, J. Assessing the representativeness and repeatability of test locations for genotype evaluation. Crop Science, v.51, p.1603-1610, 2011. DOI: https://doi.org/10.2135/cropsci2011.01.0016.
» https://doi.org/10.2135/cropsci2011.01.0016
Publication Dates
-
Publication in this collection
28 Nov 2019 -
Date of issue
2019
History
-
Received
17 Aug 2018 -
Accepted
08 Oct 2019