Acessibilidade / Reportar erro

Seleção de genótipos de trigo para rendimento de grãos e qualidade de panificação em ensaios multiambientes1 1 Este Trabalho é parte da dissertação de mestrado do primeiro autor.

Selecting wheat genotypes for yield and baking quality in multi-environment trials

Resumos

O objetivo deste estudo foi selecionar linhagens de trigo para qualidade de panificação e para rendimento de grãos (RG), em ensaios multiambientes, e identificar ambientes de teste ideais para a seleção. Foram avaliados 39 genótipos de trigo, sendo 18 linhagens e nove cultivares testemunhas, em 2010, e 11 linhagens e quatro cultivares testemunhas, em 2011, no Estado do Paraná. Algumas linhagens e locais não foram repetidos entre os dois anos. O delineamento experimental utilizado foi o de blocos casualizados, com três repetições. A metodologia, em gráfico biplot GGE (Genotype and Genotype-by-Environment), foi utilizada para a avaliação de genótipos e ambientes de teste. Na safra 2010, os locais de Nova Fátima, Pitangueiras e Ventania foram discriminantes (maximizaram a expressão do RG e diferenciação entre genótipos) e representativos do conjunto de ambientes avaliados. Em 2011, Rolândia e Astorga destacaram-se como ambientes ideais para seleção simultânea para RG e qualidade de panificação. As linhagens BIO-08528 e BIO-08228 foram classificadas como genótipos ideais para o RG e a concentração proteica dos grãos (CPG), respectivamente. Na safra 2011, as linhagens BIO-10161 e BIO-10141 foram superiores para RG e para qualidade de panificação, respectivamente. O teste de sedimentação SDS e a CPG correlacionaram-se entre si (r = 0,61**) e foram moderadamente associados com a força de glúten (r = 0,49** e 0,74**), indicando que podem ser utilizados na seleção indireta para qualidade de panificação.

Triticum aestivum L.; estabilidade; interação genótipo-ambiente; GGE biplot


The objective of this study was to select wheat lines for baking quality and grain yield in multi-environment trials in order to identify ideal test environments for selection. A total of 39 wheat genotypes, 18 inbred lines and 9 checks were evaluated in 2010, and other 11 inbred lines and 4 checks in 2011, in the state of Paraná. Some lines and locations were not repeated between the two years. The experimental design was a randomized block design with three replications. GGE biplot analysis (Genotype and Genotype-by-environment) was used to evaluate the genotypes and test environments. The locations Nova Fátima, Pitangueiras and Ventania were discriminative (maximized the expression of yield potential and the differentiation among genotypes) and representatives of the environments evaluated in the 2010 crop season. Rolândia and Astorga stood out as ideal environments to select better genotypes for baking quality in 2011. The inbred lines BIO-08528 and BIO-08228 were classified as ideal genotypes for grain yield (GY) and grain protein concentration (GPC), respectively. However, in 2011 the inbred lines BIO-10161 and BIO-10141 were superior for baking quality and GY, respectively. SDS sedimentation test and GPC were correlated each other (r= 0,61**) and moderately associated with the gluten strength (r=0,49 and** 0,74**), indicating that both traits can be used for indirect selection of baking quality.

Triticumaestivum L.; stability; genotype-by-environment interaction; GGE biplot


Introdução

O Brasil importa aproximadamente metade do trigo que consome, o que se deve, principalmente, a questões econômicas que tornam essa cultura pouco atrativa para os agricultores. Por isso, melhorar a produtividade e a qualidade de panificação é fundamental para tornar a triticultura brasileira mais competitiva. O grande desafio dos melhoristas de trigo consiste em produzir, de forma contínua, novos cultivares que combinem elevada qualidade de panificação e rendimento de grãos, atendendo, simultaneamente, as exigências dos triticultores e do mercado consumidor.

A ocorrência de interação entre genótipo e ambientes (G×A) limita a obtenção de ganhos genéticos expressivos e dificulta a seleção e recomendação de cultivares (Ahmadi et al., 2012Ahmadi J, Mohammadi A & Mirak N (2012) Targeting promising bread wheat (Triticumaestivum L.) lines for cold climate growing environments using AMMI and SREG GGE Biplot analyses. Journal of Agricultural Science and Technology, 14:645-657.; Silva & Benin, 2012Silva RR & Benin G (2012) Análises Biplot: conceitos, interpretações e aplicações. Ciência Rural, 42:1404-1412.). Sendo assim, é necessário realizar avaliações mais precisas, conduzindo ensaios em vários locais e por anos, a fim de identificar genótipos adaptados, estáveis e com características agronômicas adequadas.

Linhagens são testadas, anualmente, em múltiplos locais, por várias razões, entre as quais se destacam: I) identificar genótipos com ampla adaptabilidade com base no desempenho médio em multiambientes; II) selecionar materiais com adaptabilidade específica, com base no desempenho médio em determinados ambientes de um grupo de ambientes. Entretanto, a necessidade de extensiva avaliação em ensaios de competição implica maiores custos de manutenção de uma rede experimental.

Por esta razão, identificar ambientes que facilitem a seleção de genótipos com interação específica contribui para a otimização da rede de ensaios e maximiza a eficiência de seleção (Navabi et al. 2006Navabi A, Yanga RC, Helmb J & Spaner DM (2006) Can Spring Wheat-Growing Mega environments in the Northern Great Plains Be Dissected for Representative Locations or Niche-Adapted Genotypes? Crop Science, 46:1107-1116.; Yan et al. 2011Yan W, Pageau D, Frégeau-Reid JA & Durand J (2011) Assessing the representativeness and repeatability of test locations for genotype evaluation. Crop Science, 51:1603-1610.; Ahmadi et al., 2012Ahmadi J, Mohammadi A & Mirak N (2012) Targeting promising bread wheat (Triticumaestivum L.) lines for cold climate growing environments using AMMI and SREG GGE Biplot analyses. Journal of Agricultural Science and Technology, 14:645-657.). Embora um ambiente ideal possa não existir, ele pode ser estimado e utilizado como referência para seleção de genótipos em ensaios multiambientes (Kaya et al., 2006Kaya Y, Akçura M & Taner S (2006) GGE-Biplot Analysis of Multi-Environment Yield Trials in Bread Wheat. Turkish Journal of Agriculture and Forestry, 30:325-337.), assim como para a escolha de locais com desempenhos próximos ao ideal.

Na presença de interação de natureza complexa, ou seja, quando ocorre inversão no ranqueamento dos genótipos entre os ambientes de teste, a análise GGE (Genotype and Genotype-by-Environment) destaca-se por gerar inferências quanto ao desempenho de genótipos e ambientes, facilitando a visualização dos resultados em saídas gráficas (Silva et al., 2011Silva RR, Benin G, Silva GO, Marchioro VS, Almeida JL & Matei G (2011) Adaptabilidade e estabilidade de cultivares de trigo em diferentes épocas de semeadura, no Paraná. Pesquisa Agropecuária Brasileira, 46:1439-1447.; Benin et al., 2012Benin G, Pinnow C, Silva CL, Pagliosa ES, Beche E, Bornhofen E, Munaro LB & Silva RR (2012) Análises biplot na avaliação de cultivares de trigo em diferentes níveis de manejo. Bragantia, 71:28-36.; Silva & Benin, 2012Silva RR & Benin G (2012) Análises Biplot: conceitos, interpretações e aplicações. Ciência Rural, 42:1404-1412.). Além disso, o modelo GGE biplotidentifica ambientes representativos e discriminativos e seleciona os genótipos mais adaptados e estáveis a ambientes específicos (Asfaw et al., 2009Asfaw A, Alemayehu F, Gurum F & Atnaf M (2009) AMMI and GGE Biplot analysis for matching varieties onto soybean production environments in Ethiopia. Scientific Researchand Essay, 4:1322-1330.; Yan, 2010Yan W (2010) Optimal use of biplots in analysis of multi-location variety test data. Acta Agronomica Sinica, 36:1805-1819.).

O objetivo deste estudo foi identificar linhagens de trigo promissoras e ambientes de testes ideais, que facilitem a seleção de genótipos de trigo com superior rendimento de grãos e qualidade de panificação.

Material e Métodos

Os ensaios foram conduzidos nos anos agrícolas de 2010 e 2011, no Estado do Paraná, totalizando 11 locais de cultivo (Tabela 1). Na safra agrícola 2010, 27 genótipos de trigo (18 linhagens e nove cultivares testemunhas) foram testados em nove locais. Já em 2011, 15 genótipos (11 linhagens e quatro cultivares) foram avaliados em quatro locais de cultivo, representativos das regiões de VCU 2 e 3. A rede experimental pertencia ao programa de melhoramento de trigo da empresa Biotrigo Genética.

Tabela 1: Identificação
das linhagens, cultivares testemunhas e informações de dois locais de ensaios avaliados nas safras agrícolas de 2010 e 2011 nas regiões de Valor de Cultivo e Uso (VCU) 2 e 3

Os ensaios foram conduzidos em delineamento experimental de blocos ao acaso, com três repetições. Cada parcela constituiu-se de sete linhas de 5,0 metros de comprimento, com espaçamento entre linhas de 0,17 m, formando uma área de 6,0 m², sendo consideradas como área útil da parcela apenas as cinco linhas centrais. A semeadura foi realizada na primeira e segunda quinzena de maio, nas regiões de VCU2 e VCU3, respectivamente. A densidade de semeadura adotada foi de 350 sementes viáveis por m2.

A adubação de base utilizada foi de 30 kg N ha-1, 60 kg P2O5 ha-1,60 kg K2O ha-1 e mais 70 kg ha-1 de nitrogênio na forma de ureia (45% de N), aplicados em cobertura no início do perfilhamento (Z 2.2 da escala Zadocks). Os controles de plantas daninhas, de pragas e de doenças foram realizados de acordo com as recomendações técnicas para a cultura do trigo. O rendimento de grãos (RG) foi determinado a partir do rendimento das parcelas, corrigindo-se para a umidade de 13% (base úmida) e convertido para kg ha-1.

Para a análise de qualidade de panificação, o volume de grãos das três repetições de cada genótipo foi misturado e homogeneizado para obtenção de uma amostra composta. As amostras foram analisadas, para qualidade de panificação, no Laboratório de Cereais do Centro de Pesquisa em Alimentação da Universidade de Passo Fundo. As características viscoelásticas da farinha de trigo foram mensuradas em alveógrafo Chopin, modelo NG, conforme a AACC (1995)AACC - American Association of Cereal Chemists (1995) Approved methods. 9ª ed. Saint Paul, AACC., método no 54-30. Esse teste avalia a capacidade de absorção de água pelas proteínas formadoras de glúten e a capacidade de retenção do gás carbônico. Entre os parâmetros avaliados na alveografia, apenas a força de glúten (W) foi considerado para análise; esse caractere representa o trabalho mecânico necessário para expandir a massa até a ruptura e é expresso em 10-4 Joules. A concentração proteica dos grãos (CPG) foi determinada por espectrometria de infravermelho próximo (NIR), com aparelho Infratec Foss 1240, conforme o método nº 38-12 da AACC (1995)AACC - American Association of Cereal Chemists (1995) Approved methods. 9ª ed. Saint Paul, AACC.. O volume de sedimentação na presença de SDS (dodecil sulfato de sódio) foi determinado de acordo com o método descrito por Peña & Amaya (1985Peña RJ & Amaya A (1985) Rapid estimation of gluten quality in bread wheat, durum wheat and triticale using the sodium dodecyl sulfate (SDS) - sedimentation test: grain quality methods. Mexico, International Maize and Wheat Improvement Center. 12p.). O rendimento de proteína (RP) foi obtido a partir da multiplicação do rendimento de grãos pela concentração proteica dos grãos, de acordo com a expressão: 𝑅𝑃=𝑅𝐺𝑘𝑔ℎ𝑎−1×𝐶𝑃𝐺(%).

Por causa da prática usual de descarte de linhagens inferiores, e adição de novas nos anos seguintes, grupos diferentes de genótipos foram avaliados em 2010 e 2011. Dessa forma, fez-se necessário analisar individualmente os dados de cada ano.

A análise AMMI1(Additive Main Effects and Multiplicative Interaction) foi utilizada para estimar a contribuição relativa dos efeitos principais de ambiente, genótipo e da interação genótipo × ambiente. Essa análise considera somente o primeiro eixo de cada componente principal (CP1) e o rendimento de cada genótipo, em cada ambiente (Yij ), e foi obtida pela equação: 𝑌𝑖𝑗= 𝑌𝑖+ (𝐼𝑃𝐶𝐴1)𝑖(𝐼𝑃𝐶𝐴1)𝑗, em que Yi é a média do genótipo em todos os ambientes; IPCA1i é a pontuação do genótipo i no primeiro eixo principal da interação (𝐼𝑃𝐶𝐴1), estimado pela análise de componentes principais; IPCA1j é a pontuação do ambiente j, sobre o mesmo eixo (𝐼𝑃𝐶𝐴1). Desta forma, o rendimento de cada genótipo em cada ambiente (Yij ), é representado como uma função linear simples da pontuação ambiental (𝐼𝑃𝐶𝐴1)𝑗, isto é, uma equação reta de intercepção (IPCA1)i (Camargo-Buitrago et al., 2011Camargo-Buitrago I, Intire EQ & Gordón-Mendoza R (2011) Identificación de mega-ambientes para potenciar el uso de genotipos superiores de arroz em Panamá. Pesquisa Agropecuária Brasileira, 46:1061-1069.).

A capacidade discriminativa, a representatividade e a associação entre locais de teste foram avaliadas por meio da análise GGE (Genotype and Genotype-by-Environment). Essa análise considera apenas os efeitos de genótipo (G) e da interação genótipo × ambiente (G×E) (Yan et al., 2000Yan W K, Hunt LA, Sheng QL & Szlavnics Z (2000) Cultivar evaluation and mega-environment investigation based on the GGE biplot. Crop Science, 40:597-605.). A interpretação das associações foi realizada pelo cosseno do ângulo entre vetores, o qual equivale ao coeficiente de correlação. Foram adotados os seguintes critérios: a associação é positiva se o ângulo entre os vetores for < 90º, é negativa se o ângulo for > 90º, e é nula se o ângulo entre os vetores for igual a 90º (Yan & Tinker, 2006Yan W & Tinker A (2006) Biplot analysis of multi environment trial data: principles and applications. Canadian Journal of Plant Science, 86:623-645.). A notação Transform = 0 indica que os dados não foram transformados. No Scaling=1 (SD), os valores são divididos pelo desvio padrão de cada ambiente, conferindo a cada ambiente importância similar. A notação Centering refere-se ao modelo utilizado, em que Centering=2 contém o efeito de GGE (genótipo + interação genótipo × ambiente). O SVP, ou "singular value partitioning", é uma técnica matemática para decomposição de matrizes. O SVP = 1, com foco no desempenho genotípico, foi utilizado para a análise da associação entre os caracteres avaliados, e o SVP = 2, com foco no ambiente, foi utilizado para as demais inferências testadas (Yan et al., 2000Yan W K, Hunt LA, Sheng QL & Szlavnics Z (2000) Cultivar evaluation and mega-environment investigation based on the GGE biplot. Crop Science, 40:597-605.). Os coeficientes de correlação de Pearson foram apresentados em forma gráfica, utilizando-se o software Sigmaplot 11.0. As análises em gráfico biplot foram realizadas com o software estatístico GGE biplot (Yan, 2001Yan W (2001) GGE biplot - A Windows application for graphical analysis of multi-environment trial data and other types of two-way data. Agronomy Journal, 93:1111-1118.).

Resultados e Discussão

Os locais de Ventania (7.514 kg ha-1), Nova Fátima (6.588 kg ha-1) e Cafelândia (5.371 kg ha-1) apresentaram rendimentos de grãos (RG) superiores à média geral dos locais dos ensaios conduzidos na safra agrícola de 2010 (Tabela 2). A concentração proteica dos grãos (CPG) foi superior em Astorga, Palotina, Apucarana e Cafelândia. Os maiores valores para o rendimento de proteína (RP) foram observados em Nova Fátima, Ventania, Cafelândia, Apucarana e Astorga. Quanto à relação RG/RP, a qual indica o volume de grãos produzido para se obter um kg de proteína, os menores valores foram verificados em Palotina, Astorga e Apucarana. Na safra agrícola de 2011, Arapongas e Rolândia destacaram-se quanto ao RG (5.295 e 4.950 kg ha-1, respectivamente) e RP (603 e 592, respectivamente), sendo que Astorga e Apucarana mostraram os menores valores para a relação RG/RP. Estes locais também maximizaram a expressão da maioria dos parâmetros indicativos de qualidade de panificação. Os caracteres força de glúten (W) e CPG tiveram médias de maior magnitude em Apucarana e o caractere SDS apresentou a melhor média em Astorga.

Tabela 2:
Médias2, amplitude de variação3 e contribuição relativa dos efeitos principais (EP)4 de Genótipo (G), Ambiente (E) e interação (G×E), com base na análise AMMI1, de quatro caracteres agronômicos mensurados em 28 genótipos de trigo avaliados nas regiões de VCU 2 e 3 em 2010 e seis caracteres mensurados em 15 genótipos de trigo avaliados na região de VCU 3 na safra 2011

Considerando-se a contribuição relativa da variação proporcionada pelo desdobramento dos efeitos principais de genótipo (G), ambiente (E) e da interação G×E, por meio da análise AMMI1, verificou-se que o efeito principal de ambiente influenciou em maior magnitude o desempenho dos caracteres avaliados, corroborando os resultados reportados por Asfaw et al. (2009Asfaw A, Alemayehu F, Gurum F & Atnaf M (2009) AMMI and GGE Biplot analysis for matching varieties onto soybean production environments in Ethiopia. Scientific Researchand Essay, 4:1322-1330.) e Ahmadi et al. (2012Ahmadi J, Mohammadi A & Mirak N (2012) Targeting promising bread wheat (Triticumaestivum L.) lines for cold climate growing environments using AMMI and SREG GGE Biplot analyses. Journal of Agricultural Science and Technology, 14:645-657.). Em 2010, o efeito principal de ambiente explicou 78,92, 70,08, 69,32 e 71,28% da variação dos componentes RG, CPG, RP e RG/RP, respectivamente, sendo que o efeito de G e da interação G×E foram relativamente menos importantes (Tabela 2). Da mesma forma, em 2011 o efeito de ambiente foi responsável pela maior variação do caráter rendimento de grãos (57,91%), em comparação com o efeito de genótipo (18,17%) e com a interação G×E (13,4%). Para os demais caracteres, o efeito principal de ambiente representou aproximadamente 13, 45, 11, 45 e 27% da variação total do RP, RG/RP, SDS, CPG e W, respectivamente. A maior participação do efeito de genótipo em 2011 deve-se à proximidade entre os locais de teste, o que implicou a redução da importância relativa do ambiente sobre os genótipos avaliados.

O genótipo ideal é definido, na análise GGE, com base no alto rendimento de grãos e na estabilidade (Yan & Kang, 2003Yan W & Kang MS (2003) GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomists. 1ª ed. Florida,CRC Press. 288p.). Este "genótipo" é definido graficamente pelo vetor de maior comprimento, no PC1, e sem projeções, em PC2, ou seja, deve estar mais próximo do menor círculo concêntrico central.

Para rendimento de grãos, em 2010, a linhagem BIO 08528 foi a que mais se aproximou do "genótipo idealizado", seguida por BIO 08527, Mirante, Quartzo, TBIO Itaipu, BIO 07221, BIO 08537, BIO 08544 e BIO 08545 (Figura 1), enquanto o cultivar TBIO Ivaí e as linhagens BIO 08264, BIO 08228, BIO 08033, CD 104, BRS 208 e Marfim destacaram-se como os mais próximos ao genótipo ideal para CPG (Figura 1). No ano agrícola de 2011, as linhagens BIO 10161 e BIO 10152 foram as mais próximas do "genótipo ideal" para rendimento de grãos (Figura 2). Da mesma forma, as linhagens BIO 10120, BIO 10122, BIO 10123, BIO 10141 e o cultivar Marfim classificaram-se como mais próximas ao genótipo ideal, para parâmetros indicativos de qualidade de panificação, indicando que reúnem características de elevado rendimento de grãos e qualidade de panificação.

Figura 1:
Plotagem dos escores dos componentes principais, segundo o modelo GGE biplot, quanto a identificação do genótipo ideal para os caracteres rendimento de grãos (RG), concentração proteica dos grãos (CPG), rendimento de proteína (RP) e relação RG/RP de 28 genótipos de trigo avaliados em 9 locais de ensaio, na safra 2010.

Figura 2:
Plotagem dos escores dos componentes principais, segundo o modelo GGE biplot, quanto à identificação do genótipo ideal para os caracteres rendimento de grãos (RG), rendimento de proteína (RP), teste de sedimentação em SDS (SDS), relação RG/RP, concentração proteica dos grãos (CPG) e forca de glúten (W) de 15 genótipos de trigo.

Um ambiente de teste ideal deve discriminar os melhores genótipos e ser representativo de todos os ambientes de teste, ou seja, deve apresentar um alto escore de PC1 (maior poder de discriminação de genótipos) e pontuação nula no PC2 (maior capacidade de representar todos os outros ambientes). Da mesma forma, quanto maior o comprimento do vetor de um local de teste, maior a sua capacidade discriminativa (diferenciação de genótipos) e quanto menor o ângulo do vetor de local (em relação ao eixo de seta única), maior a representatividade deste local, em relação ao conjunto de locais avaliados (Yan & Kang, 2003Yan W & Kang MS (2003) GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomists. 1ª ed. Florida,CRC Press. 288p.; Yan & Tinker, 2006Yan W & Tinker A (2006) Biplot analysis of multi environment trial data: principles and applications. Canadian Journal of Plant Science, 86:623-645.). Ambientes discriminantes e representativos favorecem a seleção de cultivares com adaptação ampla, enquanto ambientes discriminantes, mas não representativos, podem ser úteis para descartar genótipos instáveis, da mesma forma que ambientes não discriminantes e não representativos são pouco úteis em redes experimentais de ensaios e devem ser descartados.

No ano de 2010, os locais de Pitangueiras, Nova Fátima e Ventania foram classificados como ambientes ideais de teste para rendimento de grãos e, além disso, foram simultaneamente os ambientes mais discriminantes e representativos (Figura 3). Provavelmente isso se deve à condição de irrigação (pivô central), em Pitangueiras e em Nova Fátima, e aos elevados níveis de fertilidade, em Ventania, que maximizaram a expressão do potencial de rendimento dos genótipos. Ventania mostrou-se o local menos discriminativo para CPG e para a relação RG/RP e os locais de Campo Mourão e Faxinal apresentaram informações redundantes para rendimento de grãos. A identificação e o descarte de locais redundantes reduzem os custos de pesquisa com avaliações de campo e mão de obra, sem perda de informação genotípica (Pacheco et al., 2009Pacheco RM, Duarte JB, Souza PIM, Silva AS & Nunes Junior J (2009) Key locations for soybean genotype assessment in Central Brazil. Pesquisa Agropecuaria Brasileira, 44:478-486.; Zhe et al., 2010Zhe Y, Lauer JG, Borges R & Leon N (2010) Effects of genotype × environment interaction on agronomic traits in soybean. Crop Science, 50:696-702.). Cafelândia, Faxinal e Apucarana foram classificados como ambientes ideais para avaliar CPG; porém, com alta relação RG/RP, por causa do baixo rendimento de grãos (RG) obtido nesses locais, uma vez que a CPG e o RG são inversamente associados (Bogard et al., 2010Bogard M, Allard V, Brancourt-Hulmel M, Heumez E, Machet JM, Jeuffroy MH, Gate P, Martre P & Le Gouis J (2010) Deviation from the grain protein concentration-grain yield negative relationship is highly correlated to post-anthesis N uptake in winter wheat. Journal of Experimental Botany, 61:4303-4312.; Gaju et al., 2011Gaju O, Allard V, Martre P, Snape JW, Heumez E, Le Gouis J, Moreau D, Grifiths S, Orford S, Hubbart S, Foulkes MJ & Bogard M (2011) Identification of traits to improve the nitrogen-use efficiency of wheat genotypes. Field Crops Research, 123:139-152.). Assim, fica evidente a necessidade de se identificarem genótipos que maximizem o aproveitamento de estímulos específicos de ambiente, a fim de incrementar o RG e a qualidade de panificação, simultaneamente.

avaliados em 4 locais da região de VCU3 no Estado do Paraná, na safra 2011.

Figura 3:
Plotagem dos escores dos componentes principais, segundo o modelo GGE biplot, quanto a identificação de ambiente ideal para os caracteres rendimento de grãos (RG), concentração proteica dos grãos (CPG), rendimento de proteína (RP) e relação RG/RP de 28 genótipos de trigo avaliados em 9 locais de ensaio, na safra 2010.

Na safra agrícola de 2011, Rolândia classificou-se como ambiente ideal (representativo e discriminante) para RG e CPG, enquanto Apucarana destacou-se como ambiente ideal para os parâmetros indicativos de qualidade de panificação SDS e W (Figura 4). Arapongas destacou-se como o ambiente menos representativo para o conjunto de caracteres avaliados, sendo útil apenas para a seleção de genótipos especificamente adaptados. Por terem sido testados apenas quatro locais, em 2011, não foi possível identificar ambientes redundantes.

Figura 4:
Plotagem dos escores dos componentes principais, segundo o modelo GGE biplot, quanto a identificação de ambiente ideal para os caracteres rendimento de grãos (RG), rendimento de proteína (RP), teste de sedimentação em SDS (SDS), relação RG/RP, concentração proteica dos grãos (CPG) e forca de glúten (W) de 15 genótipos de trigo avaliados em 4 locais da região de VCU3 no Estado do Paraná, na safra 2011.

As associações entre os caracteres avaliados com os genótipos e os locais de cultivo são apresentadas na Figura 5. O cosseno de um ângulo entre os vetores de quaisquer dois caracteres representa a correlação entre esses caracteres, sendo que quanto menor for o ângulo entre dois vetores, maior é a correlação entre eles (Yan & Kang, 2003Yan W & Kang MS (2003) GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomists. 1ª ed. Florida,CRC Press. 288p.). Na safra 2010, a CPG teve sua expressão maximizada (foi positivamente associada) nas condições ambientais de Palotina, Astorga e Apucarana (Figura 5a), com destaque para o cultivar Marfim e para a linhagem BIO 08264 (Figura 5b). Da mesma forma, Nova Fátima e Ventania associaram-se positivamente com os caracteres RG e RP (Figura 5a), com destaque para as linhagens BIO 08528, BIO08527 e BIO 08544 (Figura 5b). Inferências a cerca do rendimento de grãos podem ser extraídas de apenas um desses locais, em vista da associação positiva existente, conduzindo a informações seguras do conjunto de genótipos avaliados. Kaya et al. (2006Kaya Y, Akçura M & Taner S (2006) GGE-Biplot Analysis of Multi-Environment Yield Trials in Bread Wheat. Turkish Journal of Agriculture and Forestry, 30:325-337.) e Jalata (2011Jalata Z (2011) GGE-biplot Analysis of Multi-environment Yield Trials of Barley (Hordeiumvulgare L.) Genotypes in Southeastern Ethiopia Highlands. International Journal of Plant Breeding and Genetics, 5:59-75.) salientam que pode ser aplicada seleção indireta entre ambientes de testes correlacionados, obtendo-se o mesmo nível de informação de ambiente e economia, com a otimização do número de ensaios multiambientes.

Figura 5:
Plotagem dos escores dos componentes principais, segundo o modelo GGE biplot, quanto a associação entre 4 caracteres agronômicos e 9 locais de ensaios de trigo (a) e entre caracteres e genótipos (b), na safra 2010 e associação entre 8 caracteres agronômicos e 4 locais de ensaios de trigo (c) e entre caracteres e genótipos (d), na safra 2011.

Em 2010, os desempenhos dos parâmetros indicativos de qualidade de panificação CPG, SDS e W foram maximizados, nos locais de Astorga e Apucarana (Figura 5c), com destaque para as linhagens BIO10141 e BIO10120 (Figura 5d). Em 2011, os caracteres RG e CPG foram negativamente correlacionados (r = -0,53**) (Figura 6a), provavelmente por causa da competição entre carbono e nitrogênio por energia (Munier-Jolain et al., 2005Munier-Jolain NG & Salon C (2005) Are the carbon costs of seed production related to the quantitative and qualitative performance? An appraisal for legumes and other crops. Plant, Cell and Environment, 28:1388-1395.; Xu et al., 2012Xu G, Fan X & Miller AJ (2012) Plant Nitrogen Assimilation and Use Efficiency. Annual Review of Plant Biology, 63:153-82.). Outros estudos também confirmam esta relação negativa entre RG e CPG (Bogard et al., 2010Bogard M, Allard V, Brancourt-Hulmel M, Heumez E, Machet JM, Jeuffroy MH, Gate P, Martre P & Le Gouis J (2010) Deviation from the grain protein concentration-grain yield negative relationship is highly correlated to post-anthesis N uptake in winter wheat. Journal of Experimental Botany, 61:4303-4312.; Gaju et al., 2011Gaju O, Allard V, Martre P, Snape JW, Heumez E, Le Gouis J, Moreau D, Grifiths S, Orford S, Hubbart S, Foulkes MJ & Bogard M (2011) Identification of traits to improve the nitrogen-use efficiency of wheat genotypes. Field Crops Research, 123:139-152.). Entretanto, associação positiva e de elevada magnitude foi evidenciada entre os caracteres RG e RP (0,89**) (Figura 6b). Embora ocorra redução da CPG em detrimento da elevação no RG, há maior rendimento de proteína associado com o aumento do RG, quando a observação é extraída de uma unidade de área, pois ocorre uma compensação do menor valor de CPG, em virtude do maior volume de grãos.

Figura 6:
Estimativas de correlação de Pearson entre 3 caracteres agronômicos avaliados em 27 genótipos de trigo conduzidos em 9 locais de testes no Estado do Paraná, na safra 2010 (a; b) e correlação entre 4 caracteres agronômicos avaliados em 14 genótipos de trigo conduzidos em 4 locais de teste no Estado do Paraná, safra 2011 (c; d; e; f).

A força de glúten (W) foi positivamente correlacionada com a CPG (r = 0,49**) (Figura 6c) e com o teste de sedimentação em SDS (r=0,74**) (Figura 6d), corroborando outros estudos similares (Clarke et al., 2010Clarke FR, Clarke JM, Ames NA, Knox RE & Ross JR (2010) Gluten index compared with SDS-sedimentation volume for early generation selection for gluten strength in durum wheat. Canadian Journal of Plant Science, 90:01-11.; Li et al., 2013Li Y, Wu Y, Espinosa NH & Peña RJ (2013) The influence of drought and heat stress on the expression of end-use quality parameters of common wheat. Journal of Cereal Science, 57:73-78.). A CPG e o SDS (r = 0,61**) também foram positivamente associados entre si (Figura 6e). Esses resultados indicam que tanto a CPG como o SDS podem ser usados de forma eficiente na seleção indireta para estimar valores de W (Al-Saleh e Brennan, 2012Al-Saleh A & Brenna CS (2012) Bread wheat quality: some physical, chemical and rheological characteristics of Syrian and English bread wheat samples. Foods, 1:03-17.). Clarke et al (2010)Clarke FR, Clarke JM, Ames NA, Knox RE & Ross JR (2010) Gluten index compared with SDS-sedimentation volume for early generation selection for gluten strength in durum wheat. Canadian Journal of Plant Science, 90:01-11.observaram valores de correlação de 0,62 a 0,96, entre SDS e W, e de 0,52, entre SDS e CPG, e sugeriram que o teste de sedimentação SDS pode ser usado com eficácia na seleção de genótipos de trigo com elevada qualidade de panificação, pois esse teste exige baixo volume de amostra, menor necessidade de mão de obra e é relativamente barato, em comparação com as análises convencionais de W e CPG. Segundo Quick & Donnelly (1980Quick JS & Donnelly BJ (1980) A rapid test for estimating durum wheat gluten quality. Crop Science, 20:816-818.), a seleção para força de glúten com base no teste SDS pode ser 84% bem sucedida, em relação às análises tradicionais de qualidade de panificação. Além disso, um laboratorista bem treinado consegue realizar mais de 200 análises diárias.

Neste estudo, foi possível identificar ambientes discriminativos e representativos para a seleção de genótipos de trigo, com rendimentos superiores de grãos e com qualidade superior de panificação. As linhagens BIO-08228 e BIO-10161 que se destacaram com elevado rendimento de grãos e qualidade de panificação foram lançadas comercialmente como cultivares Tbio Alvorada e Tbio Noble, respectivamente. Os resultados deste estudo indicam que a CPG e SDS podem ser utilizados, de forma eficaz, na seleção de genótipos com superior força de glúten, podendo complementar ou até mesmo substituir testes alveográficos tradicionais. Porém, visando à obtenção de informações mais consistentes quanto à qualidade de panificação e a ambientes de teste, sugere-se que um maior número de locais e de anos sejam avaliados.

Conclusões

As linhagens BIO-08228, BIO-10161, BIO-08528 e BIO-10141 destacaram-se simultaneamente como genótipos ideais para o rendimento de grãos e a qualidade de panificação.

Os locais de teste Pitangueiras, Nova Fátima, Rolândia e Astorga foram ambientes representativos e discriminantes do rendimento de grãos e de parâmetros indicativos da qualidade de panificação.

A seleção para qualidade de panificação pode ser feita indiretamente, por meio da concentração proteica dos grãos ou teste de sedimentação em SDS.

Referências

  • AACC - American Association of Cereal Chemists (1995) Approved methods. 9ª ed. Saint Paul, AACC.
  • Ahmadi J, Mohammadi A & Mirak N (2012) Targeting promising bread wheat (Triticumaestivum L.) lines for cold climate growing environments using AMMI and SREG GGE Biplot analyses. Journal of Agricultural Science and Technology, 14:645-657.
  • Al-Saleh A & Brenna CS (2012) Bread wheat quality: some physical, chemical and rheological characteristics of Syrian and English bread wheat samples. Foods, 1:03-17.
  • Asfaw A, Alemayehu F, Gurum F & Atnaf M (2009) AMMI and GGE Biplot analysis for matching varieties onto soybean production environments in Ethiopia. Scientific Researchand Essay, 4:1322-1330.
  • Benin G, Pinnow C, Silva CL, Pagliosa ES, Beche E, Bornhofen E, Munaro LB & Silva RR (2012) Análises biplot na avaliação de cultivares de trigo em diferentes níveis de manejo. Bragantia, 71:28-36.
  • Bogard M, Allard V, Brancourt-Hulmel M, Heumez E, Machet JM, Jeuffroy MH, Gate P, Martre P & Le Gouis J (2010) Deviation from the grain protein concentration-grain yield negative relationship is highly correlated to post-anthesis N uptake in winter wheat. Journal of Experimental Botany, 61:4303-4312.
  • Camargo-Buitrago I, Intire EQ & Gordón-Mendoza R (2011) Identificación de mega-ambientes para potenciar el uso de genotipos superiores de arroz em Panamá. Pesquisa Agropecuária Brasileira, 46:1061-1069.
  • Clarke FR, Clarke JM, Ames NA, Knox RE & Ross JR (2010) Gluten index compared with SDS-sedimentation volume for early generation selection for gluten strength in durum wheat. Canadian Journal of Plant Science, 90:01-11.
  • Gaju O, Allard V, Martre P, Snape JW, Heumez E, Le Gouis J, Moreau D, Grifiths S, Orford S, Hubbart S, Foulkes MJ & Bogard M (2011) Identification of traits to improve the nitrogen-use efficiency of wheat genotypes. Field Crops Research, 123:139-152.
  • Jalata Z (2011) GGE-biplot Analysis of Multi-environment Yield Trials of Barley (Hordeiumvulgare L.) Genotypes in Southeastern Ethiopia Highlands. International Journal of Plant Breeding and Genetics, 5:59-75.
  • Kaya Y, Akçura M & Taner S (2006) GGE-Biplot Analysis of Multi-Environment Yield Trials in Bread Wheat. Turkish Journal of Agriculture and Forestry, 30:325-337.
  • Li Y, Wu Y, Espinosa NH & Peña RJ (2013) The influence of drought and heat stress on the expression of end-use quality parameters of common wheat. Journal of Cereal Science, 57:73-78.
  • Munier-Jolain NG & Salon C (2005) Are the carbon costs of seed production related to the quantitative and qualitative performance? An appraisal for legumes and other crops. Plant, Cell and Environment, 28:1388-1395.
  • Navabi A, Yanga RC, Helmb J & Spaner DM (2006) Can Spring Wheat-Growing Mega environments in the Northern Great Plains Be Dissected for Representative Locations or Niche-Adapted Genotypes? Crop Science, 46:1107-1116.
  • Pacheco RM, Duarte JB, Souza PIM, Silva AS & Nunes Junior J (2009) Key locations for soybean genotype assessment in Central Brazil. Pesquisa Agropecuaria Brasileira, 44:478-486.
  • Peña RJ & Amaya A (1985) Rapid estimation of gluten quality in bread wheat, durum wheat and triticale using the sodium dodecyl sulfate (SDS) - sedimentation test: grain quality methods. Mexico, International Maize and Wheat Improvement Center. 12p.
  • Quick JS & Donnelly BJ (1980) A rapid test for estimating durum wheat gluten quality. Crop Science, 20:816-818.
  • Silva RR & Benin G (2012) Análises Biplot: conceitos, interpretações e aplicações. Ciência Rural, 42:1404-1412.
  • Silva RR, Benin G, Silva GO, Marchioro VS, Almeida JL & Matei G (2011) Adaptabilidade e estabilidade de cultivares de trigo em diferentes épocas de semeadura, no Paraná. Pesquisa Agropecuária Brasileira, 46:1439-1447.
  • Yan W, Pageau D, Frégeau-Reid JA & Durand J (2011) Assessing the representativeness and repeatability of test locations for genotype evaluation. Crop Science, 51:1603-1610.
  • Yan W (2001) GGE biplot - A Windows application for graphical analysis of multi-environment trial data and other types of two-way data. Agronomy Journal, 93:1111-1118.
  • Yan W (2010) Optimal use of biplots in analysis of multi-location variety test data. Acta Agronomica Sinica, 36:1805-1819.
  • Yan W & Kang MS (2003) GGE biplot analysis: A graphical tool for breeders, geneticists, and agronomists. 1ª ed. Florida,CRC Press. 288p.
  • Yan W K, Hunt LA, Sheng QL & Szlavnics Z (2000) Cultivar evaluation and mega-environment investigation based on the GGE biplot. Crop Science, 40:597-605.
  • Yan W & Tinker A (2006) Biplot analysis of multi environment trial data: principles and applications. Canadian Journal of Plant Science, 86:623-645.
  • Xu G, Fan X & Miller AJ (2012) Plant Nitrogen Assimilation and Use Efficiency. Annual Review of Plant Biology, 63:153-82.
  • Zhe Y, Lauer JG, Borges R & Leon N (2010) Effects of genotype × environment interaction on agronomic traits in soybean. Crop Science, 50:696-702.

Datas de Publicação

  • Publicação nesta coleção
    Jul-Aug 2015

Histórico

  • Recebido
    26 Nov 2013
  • Aceito
    14 Jul 2015
Universidade Federal de Viçosa Av. Peter Henry Rolfs, s/n, 36570-000 Viçosa, Minas Gerais Brasil, Tel./Fax: (55 31) 3612-2078 - Viçosa - MG - Brazil
E-mail: ceres@ufv.br