SciELO - Scientific Electronic Library Online

 
vol.40 issue1Estimate of post-harvest losses in potato (Solanum tuberosum): A case-study at Distrito de Moamba, Bloco I (Moamba-Sede)Hydroretentor polymers in soybean in Latosol conditions author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista de Ciências Agrárias

Print version ISSN 0871-018X

Rev. de Ciências Agrárias vol.40 no.1 Lisboa Mar. 2017

https://doi.org/10.19084/RCA16013 

ARTIGO

Divergência genética entre genótipos de milho (Zea mays L.) em ambientes distintos

Genetic divergence among corn (Zea mays L.) genotypes in distinct environments

Maicon Nardino1*, Diego Baretta1, Ivan R. Carvalho1, Diego N. Follmann2, Mauricio Ferrari1, Alan J. de Pelegrin1, Vinícius J. Szareski3, Valmor A. Konflanz4 e Velci Q. de Souza5

 

1 Departamento de Agronomia, Centro de Genômica e Fitomelhoramento, Faculdade de Agronomia Eliseu Maciel – Universidade Federal de Pelotas, C.P. 354, 96010-900, Capão do Leão, Brasil; *E-mail: nardinomn@gmail.com

2 Programa de Pós-Graduação em Agronomia, Universidade Federal de Santa Maria, Departamento de Fitotecnia, CEP 97105-900, Santa Maria, Rio Grande do Sul, Brasil;

3 Programa de Pós Graduação em Ciência e Tecnologia de Sementes, Departamento de Fitotecnia, Faculdade de Agronomia Eliseu Maciel, C.P. 354, 96010-900, Capão do Leão, Brasil;

4 KSP Sementes e Pesquisas Ltda, CEP 85504-050, Pato Branco – PR, Brasil;

5 Universidade Federal do Pampa, Campus Dom Pedrito – RS, C.P. 96450-000, Dom Pedrito, Brasil.

 

RESUMO

O objetivo deste trabalho foi caracterizar morfologicamente a dissimilaridade genética entre 25 genótipos de milho (Zea mays L.), visando a seleção dos mais divergentes e de maior potencial produtivo, por meio da análise multivariada, em cinco ambientes de cultivo do Sul do Brasil. As linhagens endogâmicas dos dois grupos heteróticos foram provenientes do banco de germoplasma da empresa KSP Sementes Ltda. Utilizaram-se 15 linhagens como genitores femininos e oito linhagens como genitores masculinos. Os ensaios foram conduzidos na safra agrícola 2011/2012 em cinco ambientes, com delineamento de blocos completos ao acaso com três repetições por local. Foram avaliados 15 caracteres de interesse agronômico. A divergência genética foi avaliada por meio do método de agrupamento de Tocher e UPGMA. As estimativas da D2 indicam os pares de genótipos 6 e 22 como mais distantes geneticamente, e 1 e 14 os pares mais similares. Foram identificados nove grupos divergentes em ambos os métodos de agrupamento. Os caracteres profundidade de grãos, diâmetro do sabugo, massa de sabugo, diâmetro de espiga e peso de 100 grãos foram os principais determinantes na quantificação da divergência genética. Para o melhoramento de milho, os genótipos com maiores potenciais são 5, 15, 12, 20, 9, 25 e 10.

Palavra-chave: Análise multivariada, dissimilaridade genética, desempenho agronômico

 

ABSTRACT

The aim of this study was to morphologically characterize the genetic dissimilarity among 25 maize (Zea mays L.) genotypes in order to select the most divergent and higher yield potential, by multivariate analysis, in five cultivation environments in southern Brazil. The inbred lines of the two heterotic groups were from the genebank company KSP Seeds Ltda. Fifteen inbred lines were used as female parents and eight inbred lines as male parents. The tests were conducted during the harvest 2011/2012 in five environments, in randomized complete block design with three replications per location. We evaluated 15 characters of agronomic interest. Genetic divergence was assessed by the clustering method of Tocher and UPGMA. Estimates of D2 indicated genotype pairs 6 and 22 as genetically more distant, and 1 and 14 the most similar pairs. Nine different groups were identified in both clustering methods. The depth of grain, cob diameter, cob weight, ear diameter and weight of 100 grains characters were the main determinants in the quantification of genetic divergence. For breeding purposes, genotypes 5, 15, 12, 20, 9, 25 and 10 are those of higher potential.

Key words: genetic dissimilarity, multivariate analysis, agronomic performance.

 

INTRODUÇÃO

O milho (Zea mays L.) destaca-se economicamente como uma importante cultura, empregada para fonte de alimento, fibras, combustível e rações. O cereal é usado como ingrediente em uma inúmera lista de produtos manufaturados e empregados na nutrição da população mundial (Hallauer et al., 2010). O Brasil se destaca como um dos maiores produtores de milho do mundo, superado apenas pelos Estados Unidos e China (Faostat, 2014). De acordo com levantamento realizado pela Conab, na safra de 2013/2014, a produção brasileira de milho foi de 78,96 milhões de toneladas em 15,49 milhões de hectares, com uma produtividade média de 5.095 kg ha-1 (Conab, 2014). Entretanto, algumas regiões brasileiras atingem produtividades, superiores a 8.000 kg por hectare (Dotto et al., 2010).

A avaliação da divergência genética é amplamente utilizada pelos melhoristas de milho para seleção de genitores. Este método visa selecionar materiais mais promissores, diminuindo custos e tempo necessário para realização de vários cruzamentos, muitas vezes desnecessários. Dessa forma, esforços são concentrados em combinações híbridas entre genitores mais divergentes (Rinaldi et al., 2007). O conhecimento da divergência genética fornece parâmetros para a correta escolha de genitores que, ao serem cruzados, resultem em elevado efeito heterótico nas progênies, maximizando as chances de obtenção de genótipos superiores em gerações segregantes (Rotili et al., 2012).

A divergência genética tem sido avaliada por meio de técnicas biométricas, baseadas na quantificação da heterose, destacando-se as análises dialélicas, ou por processos preditivos. Nos dialelos, quando o número de genitores é elevado, a obtenção de material experimental pode ser impraticável. Por dispensarem a obtenção prévia das combinações híbridas, os métodos preditivos tem merecido considerável ênfase. Na predição de divergência genética, vários métodos de análise multivariada podem ser aplicados, entre os quais citam-se a análise de componentes principais, variáveis canônicas e os métodos aglomerativos. Os métodos aglomerativos dependem de medidas de dissimilaridades estimadas previamente, como a distância euclidiana ou a distância generalizada de Mahalanobis. Contudo, a distância de Mahalanobis tem sido amplamente utilizada, por levar em consideração as matrizes de covariância residuais estimadas a partir de ensaios experimentais com repetição (Cruz et al., 2012).

A divergência genética tem sido estudada em várias culturas com base em caracteres morfoagronômicos, moleculares e ambos, visando a seleção de genitores para formação de combinações híbridas e/ou formação de novas populações segregantes oriundas de cruzamentos com genótipos mais divergentes (Rinaldi et al., 2007; Ganesan et al., 2010; Vogt et al., 2010; Quintal et al., 2012; Rotili et al., 2012; Silva et al., 2012; Kage et al., 2013). Quanto mais divergentes forem os genitores, maior é a variabilidade resultante na população segregante, e maior é a probabilidade de reagrupar alelos em novas combinações favoráveis (Barbieri et al., 2005).

O objetivo deste trabalho foi caracterizar morfologicamente a dissimilaridade genética entre 25 genótipos de milho, visando a seleção dos mais divergentes e de maior potencial produtivo, por meio de análise multivariada, em cinco ambientes de cultivo no Sul do Brasil.

 

MATERIAL E MÉTODOS

Os cruzamentos artificiais entre as linhagens endogâmicas foram conduzidos na estação de pesquisa da empresa KSP Sementes Ltda. no município de Clevelândia-PR no ano agrícola de 2010/2011. Para condução do trabalho foram utilizadas 15 linhagens endogâmicas como genitores femininos, que apresentavam características de elevado rendimento de grão e oito como genitores masculinos, com características de rusticidade, formando dois grupos heteróticos distintos. Os grupos heteróticos foram cruzados artificialmente conforme o modelo dialélico parcial, resultando consequentemente, em 25 combinações híbridas. O germoplasma utilizado foi oriundo do programa de melhoramento genético da empresa KSP Sementes Ltda., com sede no município de Pato Branco - PR. As sementes híbridas obtidas dos cruzamentos foram colhidas manualmente, secadas e preparadas para a sementeira dos ensaios nos cinco locais.

A sementeira foi conduzida na safra agrícola 2011/2012 em cinco locais distintos dos três estados da região Sul do Brasil. No Rio Grande do Sul, o ensaio foi conduzido no Município de Frederico Westphalen, sendo as coordenadas (27°23'47''S, 53°25'35''O, 480 m). Em Santa Catarina, o ensaio foi conduzido no município de Itapiranga, sendo as coordenadas (27º10'10"S, 53º42'44"O, 206 m). No Paraná, os ensaios foram conduzidos em três locais: Pato Branco, com as coordenadas (26°13'44?S, 52°40'15?O, 760 m), Ampére com as coordenadas (25°54'20''S, 53º25'54'' O, 718 m), Clevelândia com as coordenadas (26°21'17''S, 52°28'56'' O, 860 m).

As unidades experimentais foram compostas por duas linhas de cinco metros de comprimento, espaçadas 0,70 metros. As linhas de sementeira foram demarcadas com semeadoura de sistema de sementeira direta, realizando-se na mesma operação a adubação das linhas. As sementes híbridas provenientes dos cruzamentos dos dois grupos heteróticos foram semeadas manualmente. A época de sementeira foi realizada de acordo com o zoneamento agroclimático de cada local. O manejo do solo, e os tratos culturais foram os mesmos para os cinco locais obedecendo aos estádios fenológicos e às necessidades da cultura, onde após a emergência e estabelecimento da mesma, realizou-se o raleio manual, para ajuste do estande de 42 plantas por unidade experimental, equivalente a 60.000 plantas ha-¹.

Os seguintes caracteres agronômicos foram mensurados: altura de planta (AP, em cm), altura de inserção da espiga (AE, em cm), prolificidade (PROL, em unidades), rendimento de grãos de parcela (RG, em kg), diâmetro da espiga (DE, em mm), comprimento da espiga (CE, em cm), número de fileiras por espiga (NF, em unidades), número de grãos por fileira (NC, em unidades), número total de grãos por espiga (NG, em unidades), peso da espiga (PE, em g), profundidade de grãos (PFG, em mm), peso total de grãos da espiga (PG, em g), diâmetro do sabugo (DS, em mm), peso do sabugo (PS, em g) e peso de 100 grãos (MCG, em g).

Na análise dos dados, foi realizada a análise de variância individual e, constatada a homogeneidade de variâncias residuais pelo quociente entre o maior e menor quadrado médio residual ter ficado abaixo de 7 (Gomes e Garcia, 2002), realizou-se a análise de variância conjunta. Os dados foram submetidos à análise de agrupamento de médias pelo teste Scott e Knott (1974) a 5%, para o caráter RG. Na aplicação da técnica de agrupamento, adotou-se a distância generalizada de Mahalanobis (D2), como medida de dissimilaridade, sendo efetuada a importância relativa dos caracteres, estimada pela partição dos componentes de D2, relativos a cada caráter, no total da dissimilaridade observada (Singh, 1981). No estudo da divergência genética, foi empregado o método de otimização de Tocher (Rao, 1952) e análise de agrupamento pelo método das médias aritméticas não-ponderadas (Unweighted pair group method using arithmetic average - UPGMA). Todas as análises foram realizadas através do programa computacional GENES (Cruz, 2013).

 

RESULTADOS E DISCUSSÃO

A variação associada a híbridos, locais de cultivo e interação híbridos x locais de cultivo revelou significância (p<0,05) para todos os caracteres analisados, indicando que o comportamento entre os mesmos diferiu em função das variações ambientais existentes entre os locais de cultivo, possibilitando a identificação de constituições genéticas de melhor desempenho per se sobre os caracteres analisados (Quadro 1). Os coeficientes de variação apresentaram amplitude de 3,13% para DE até 7,69% para PFG, indicando boa precisão experimental segundo a escala proposta por Scapim et al. (1995).

Os resultados da análise de variância multivariada (MANOVA), referente a 15 caracteres avaliados em 25 híbridos de milho em cinco locais de cultivo, encontram-se sumarizados no Quadro 2. A análise MANOVA evidenciou, da mesma forma, diferenças significativas entre os vetores de médias para todos os efeitos principais (híbridos e locais) e para o efeito da interação híbrido x locais. Dessa forma, com a rejeição da hipótese de que os vetores de médias são iguais, justifica-se o uso de outras técnicas multivariadas, objetivando estudos de dissimilaridade genética.

As medidas de dissimilaridade estimadas a partir da distância generalizada de Mahalanobis (D2) relativa aos 25 híbridos nos cinco locais são indicadas no Quadro 3. Dentre os híbridos avaliados, os menores valores de D2 foram observados entre KSP1 e KSP14 (2,3), e a maior distância ocorreu entre os híbridos KSP6 e KSP22 (77,2). A relação entre o maior e o menor valor observado de D2 foi da ordem de 33,6 caracterizando, assim, a presença de variabilidade genética entre os materiais estudados, corroborando os resultados obtidos por Paixão et al. (2008) que estudando seis populações de milho obteve como estimativa de D2 o valor de 39.

As amplitudes dos valores das estimativas de dissimilaridade observadas sugerem que há variabilidade genética entre o conjunto de híbridos analisados, visto que, tais constituições híbridas foram originadas de distintos grupos heteróticos, viabilizando desta maneira, a extração de linhagens endogâmicas e/ou possibilitando ganho heterótico com exploração de híbridos duplos. Segundo Cruz e Carneiro (2006) dois pais não distantes geneticamente entre si tendem a compartilhar muitos genes e alelos em comum, apresentando pouca complementariedade quando os mesmos são cruzados, apresentando baixo vigor como consequência do baixo nível de heterozigosidade alélica. Entretanto, quando dois pais são mais distantes geneticamente, é admitido que eles difiram no número de loci no qual o efeito de dominância é evidente, refletindo uma maior manifestação da heterose.

O método de otimização de Tocher, utilizando as médias dos híbridos nos cinco ambientes (análise conjunta), proporcionou a formação de nove grupos distintos, onde o grupo I formou o maior número de genótipos, representado por 11 híbridos, que juntos incorporaram 44% do conjunto dos 25 híbridos estudados, seguido do grupo II que apresentou quatro híbridos, grupo III, IV e V com dois híbridos e grupos VI, VII, VIII e IX com um híbrido cada (Quadro 4). Pode verificar-se elevada dissimilaridade entre o conjunto de híbridos, explicando a formação de nove grupos. A elevada dissimilaridade genética observada entre os genótipos é principalmente devido à origem e composição genética das linhagens endogâmicas utilizadas, pertencentes aos grupos heteróticos distintos. Os grupos de híbridos formados apresentam menor distância média intragrupo que a distância média intergrupo (Cruz e Carneiro, 2006).

Segundo Simon et al. (2012) em programas de melhoramento genético de plantas com ênfase para produção de híbridos, a escolha de genitores para extração de linhagens, que sejam complementares para um maior conjunto de caracteres agronômicos de interesse é de extrema relevância.

A dissimilaridade intergrupo (Quadro 5) obtida pelo método de otimização de Tocher, permitiu discriminar entre os grupos formados, quais os que se revelaram mais divergentes geneticamente. A menor distância intragrupo foi revelada no grupo III (9,010), em contrapartida a maior distância foi verificada pelo grupo V (12,437). Para os demais grupos não foi possível estimar as distâncias intragrupos, por ambos possuírem somente um representante cada. Por sua vez, as maiores distâncias médias intergrupos observadas foram entre os grupos V e IX (62,374), seguido do grupo III e IX (56,951), enquanto que as menores distâncias intergrupos foram evidenciadas entre o grupo I e III (14,002) e o grupo III e IV (14,177). O grupo I revelou o maior número de genótipos agrupados, apresentando a maior distância média em relação ao grupo IX (33,035), seguido do grupo VIII (23,387) e IV (19,846) e proporcionando dessa maneira, a formação de distintos grupos heróticos. A correta alocação dos grupos heteróticos é fundamental, pois segundo Ferreira et al. (1995) a capacidade de combinação e a heterose são diretamente proporcionais a dissimilaridade genética, e dessa forma, maiores são as chances de alcançar combinações híbridas promissoras.

Nos métodos hierárquicos, os genótipos são agrupados por um processo que se repete por vários níveis, até que se estabeleça um dendrograma ou o diagrama de árvore (Cruz e Carneiro, 2006). Dessa forma, o método UPGMA aliado ao método de otimização de Tocher, confere maior eficiência em discriminar os genótipos quanto as suas distâncias genéticas. O método UPGMA (Figura 1) permitiu a formação de nove grupos distintos, sendo considerado como ponto de corte a 50% do eixo das abscissas, que foram: grupo G1 (KSP1, KSP14, KSP7, KSP19, KSP3, KSP8, KSP16, KSP2 e KSP17); grupo G2 (KSP5, KSP15, KSP9 e KSP20); grupo G3 (KSP4 e KSP6); grupo G4 (KSP10), grupo G5 (KSP11, KSP24 e KSP18); grupo G6 (KSP13); grupo G7 (KSP12, KSP21 e KSP25); grupo G8 (KSP22) e grupo G9 (KSP23). O coeficiente de correlação cofenético, que estima a representatividade dos dados da matriz de dissimilaridade do dendrograma, revelou magnitude de 0,60, indicando que os dados da matriz apresentaram ajuste satisfatório na representação gráfica apresentada pelo dendrograma.

Observa-se que ambos os métodos (Tocher e UPGMA) foram parcialmente concordantes quanto à dissimilaridade genética entre o conjunto de genótipos, com a formação do mesmo número de grupos. Da mesma forma, foi observado que a maior proporção dos genótipos foram alocados no primeiro grupo, em ambos métodos de agrupamento, assim como grupos isolados com apenas um híbrido. Estas diferenças devem-se, segundo Büttow et al. (2010), à distinta maneira de cada método realizar o cálculo da variabilidade genética. A semelhança de ambos os métodos de agrupamento na discriminação de genótipos pela dissimilaridade genética corrobora os estudos realizados por outros autores (Neitzke et al., 2009; Büttow et al., 2010; Simon et al., 2012).

Segundo Carpentieri-Pípolo et al. (2000) a seleção de genótipos para cruzamento com base somente na divergência genética, sem considerar seus rendimentos, pode não ser uma estratégia eficiente para um programa de melhoramento. Dessa forma, além da dissimilaridade genética para seleção de genitores, o desempenho per se dos mesmos, nos ambientes específicos, merece atenção. De acordo com a análise de agrupamento de médias obtidas pelo teste de Scott e Knott (1974) (Quadro 6), os híbridos com maior rendimento de grão (RG) para o local Ampére foram KSP8, KSP24, KSP5, KSP16 e KSP17. Destacaram-se em Clevelândia os híbridos KSP5, KSP15, KSP21, KSP1, KSP10, KSP8, KSP9, KSP20, KSP14, KSP4, KSP24, KSP3 e KSP22 que apresentaram valores de rendimento de grão (RG) superior a 10 Kg por parcela. Em Itapiranga os híbridos KSP12, KSP15 e KSP20 formaram o grupo mais produtivo. Para Pato Branco o híbrido KSP12 obteve o melhor desempenho, com produtividade superior a 11 Kg por parcela. Em Frederico Westphalen os híbridos KSP9, KSP15, KSP13 e KSP19 formaram o grupo mais produtivo. Em relação à média conjunta, envolvendo os cinco ambientes, o caráter RG evidenciou uma variação de 8,82 Kg para o híbrido KSP 5 a 6,75 Kg para o híbrido KSP 18, representando uma amplitude de variação de 30,7%, destacando os híbridos KSP5, KSP15, KSP12, KSP20, KSP9, KSP25 e KSP10 que formaram o grupo de produtividade superior.

Embora a interação genótipo x ambiente tenha sido significativa, verificou-se que o desempenho de um genótipo pode ser superior num determinado ambiente não se mantendo em todos os locais em que a experimentação decorre. Dotto et al. (2010) reporta que, levando em consideração a demanda por mão-de-obra no desenvolvimento de programas de melhoramento de milho, específicos para uma micro-região, ser uma tarefa onerosa, os primeiros cruzamentos necessitam ser amplos. Dessa forma, a melhor estratégia seria utilizar informações provenientes de análises conjuntas dos ambientes, podendo no futuro serem feitas avaliações e seleções das cultivares comerciais para ambientes específicos.

O híbrido KSP22, formando o grupo IX e G8 pelo método de Tocher e UPGMA respectivamente, apesar de destacar-se pela elevada divergência dos demais grupos, não apresentou resultados desejáveis em relação ao seu desempenho per se nos ambientes individuais e de forma conjunta, mostrando não ser um genitor em potencial, visto que, a recomendação de cruzamentos, seja realizada entre genótipos divergentes, e que apresentem superioridade para caracteres de importância agronômica (Cruz et al., 2012).

As estimativas da contribuição relativa dos 15 caracteres para a dissimilaridade genética entre 25 híbridos simples de milho foram avaliadas pelo método de Singh (1981) estão apresentadas no Quadro 7. Os caracteres que mais contribuíram para divergência genética levando em consideração a análise conjunta dos cinco ambientes foram: a profundidade de grão (PGF - 33,33%), o diâmetro do sabugo (DS - 18,22%), o peso do sabugo (PS - 9,62%), o diâmetro da espiga (DE - 9,01%) e o peso de 100 grãos (MCG - 6,66%), totalizando 76,84% da divergência genética entre os híbridos. Já os caracteres prolificidade (PROL - 0,54%), número de grãos por fileira (NC - 1,09%), altura da planta (AP - 1,23%), número de fileiras (NF - 1,37%) e peso da espiga (PE - 1,51%) foram os que menos contribuíram para a divergência genética. Segundo Dotto et al. (2010) é fundamental conhecer a influência das principais características na divergência genética entre os materiais, pois estas irão direcionar os trabalhos de melhoramento, assim como aumentar a eficiência na síntese de híbridos promissores. Azad et al. (2012) observaram, em 30 linhagens de milho, que o diâmetro da espiga foi um dos caracteres com maior contribuição para a divergência genética. Kage et al. (2013) reportaram em um estudo envolvendo 82 linhagens endogâmicas de milho, que a altura da planta e a altura da inserção da espiga foram os caracteres que mais contribuíram na divergência genética com 30,6% e 16,65% respectivamente, sendo os caracteres rendimento de grão (8,70%), número de grãos por fileira (6,14%), número de grãos por espiga (1,96%) e comprimento da espiga (1,48%) os que tiveram menores contribuições. Para Nagalakshmi et al. (2010) a contribuição depende principalmente dos genótipos utilizados no estudo e nas influências ambientais sobre o caráter.

Os resultados revelam perspectivas de trabalhos futuros com a exploração da variabilidade genética entre o conjunto de híbridos avaliados, visando a obtenção de híbridos duplos e/ou formação de futuras populações para extração de linhagens. Genitores que apresentaram destaque podem ser incluídos em cruzamentos dialélicos, por exemplo, para se ter confirmação dos resultados pela análise multivariada, uma vez que são de natureza preditiva, permitindo a estimação de parâmetros como heterose e a capacidade de combinação e maximizando as chances de cruzamentos mais promissores.

 

CONCLUSÃO

Os cruzamentos entre os híbridos KSP5, KSP15, KSP12, KSP20, KSP9, KSP25 e KSP10 podem resultar em novas combinações gênicas por serem divergentes e reunirem desempenho per se desejável.

Os caracteres profundidade de grãos, diâmetro do sabugo, massa de sabugo, diâmetro de espiga e peso de 100 grãos são os principais determinantes na quantificação da divergência genética entre os genótipos de milho.

 

REFERÊNCIAS BIBLIOGRÁFICAS

Azad, M.A.K.; Biswas, B.K.; Alam, N. e Alam, S.S. (2012) - Diversity in maize (Zea mays L.) inbred lines. The Agriculturists, vol. 10, n. 1, p. 64-70. http://dx.doi.org/10.3329/agric.v10i1.11066        [ Links ]

Barbieri, R.L.; Leite, D.L.; Choer, E. e Sinigaglia, C. (2002) - Divergência genética entre populações de cebola com base em marcadores morfológicos. Ciência Rural, vol. 35, n. 2, p. 303-308. http://dx.doi.org/10.1590/S0103-84782005000200009        [ Links ]

Vogt, G.A.; Balbinot Júnior, A.A. e Souza, A.M.D. (2010) - Divergência genética entre cultivares de girassol no planalto norte catarinense. Scientia Agrária, vol. 11, n. 4, p. 307-315. http://dx.doi.org/10.5380/rsa.v11i4.18265        [ Links ]

 

Recebido/received: 2016.01.29

Recebido em versão revista/received in revised form: 2016.05.19

Aceite/accepted: 2016.05.31

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License