Identidade de modelos não lineares para comparar curvas de crescimento de bovinos da raça Tabapuã

O objetivo deste trabalho foi avaliar o uso de teste de identidade de modelos não lineares, na comparação de curvas de crescimento de bovinos da raça Tabapuã, de cinco regiões de produção do Nordeste do Brasil. Foram analisados dados de peso de 3.695 machos e 4.236 fêmeas, originários das regiões Maranhão, Gado Algodão, Mata Agreste, Sertão e Itapetinga-Valadares. Após ajuste do modelo Brody, aplicou-se o teste da razão de verossimilhança, com aproximação de qui-quadrado, para avaliar a igualdade de parâmetros de curvas de crescimento entre as regiões. O modelo reduzido, com igualdade da taxa de maturidade, com 14 parâmetros, foi o mais adequado para descrever o crescimento dos animais. As curvas de crescimento dos machos apresentaram taxas de maturidade em comum nas regiões de produção Gado Algodão e Mata Agreste, Maranhão e Itapetinga-Valadares, e Sertão. Para as fêmeas, as regiões com taxas de maturidade em comum foram: Mata Agreste e Sertão, Maranhão e Itapetinga-Valadares, e Gado Algodão. A utilização de uma única curva não é adequada para descrever o crescimento de bovinos da raça Tabapuã nas regiões estudadas.

Estimação de parâmetros genéticos para peso em diferentes idades para animais da raça Tabapuã

Foram estimados parâmetros genéticos para pesos do nascimento aos 570 dias de idade para 35.308 animais da raça Tabapuã, nascidos entre 1975 e 2000, pertencentes à ABCZ (Associação Brasileira de Criadores de Zebu) sob três modelos distintos. O modelo 1 incluiu o efeito genético aditivo direto como aleatório, além dos efeitos fixos de grupo de contemporâneos, definido pelas variáveis: proprietário, rebanho, criador, rebanho do criador, sexo, condição de criação, ano e mês de nascimento, ano e mês da pesagem, e os efeitos linear e quadrático de idade do bezerro e idade da vaca ao parto como covariáveis. O modelo 2 compreendeu, além dos efeitos supracitados, o efeito de ambiente permanente materno. O modelo 3 constou dos efeitos genético aditivo direto e materno e de ambiente permanente materno (aleatórios) e os mesmos incluídos no modelo 1 (fixos). de acordo com o teste de razão de verossimilhança, o modelo 3 foi o mais adequado para ajustar os efeitos estudados. As estimativas de herdabilidade direta foram baixas a moderadas (0,08 a 0,26), decrescendo do nascimento às idades subseqüentes, com picos aos 90 e aos 180 dias de idade. Aos 345 dias de idade, ocorreu novo aumento nas estimativas, com menor oscilação entre as estimativas subseqüentes até 570 dias de idade. As estimativas de herdabilidade materna foram baixas, sendo maiores no período da desmama.

Estimação de parâmetros genéticos para peso do nascimento aos 550 dias de idade para animais da raça Tabapuã utilizando-se modelos de regressão aleatória

Foram utilizados 21.762 registros de peso do nascimento aos 550 dias de idade de 4.221 animais para estimativa das funções de covariância empregando modelos de regressão aleatória. Os modelos incluíram, como aleatórios, os efeitos genéticos aditivo direto e materno, de ambiente permanente de animal e de ambiente permanente materno e, como fixos, os efeitos de grupo contemporâneo, a idade da vaca ao parto (linear e quadrático) e o polinômio ortogonal de Legendre da idade do animal (regressão cúbica), como covariáveis. As variâncias residuais foram modeladas por uma função de variâncias com ordens de 2 a 6. Análises com polinômios ortogonais de diversas ordens foram realizadas para os efeitos genético aditivo direto, genético aditivo materno, de ambiente permanente de animal e de ambiente permanente materno. Os modelos foram comparados pelos critérios de informação Bayesiano de Schwarz (BIC) e Akaike (AIC). O melhor modelo indicado por todos os critérios foi o que considerou o efeito genético aditivo direto ajustado por um polinômio cúbico, o efeito genético materno ajustado por um polinômio quadrático, o efeito de ambiente permanente de animal ajustado por polinômio quártico e o efeito de ambiente permanente materno ajustado por polinômio linear. As estimativas de herdabilidade para o efeito direto foram maiores no início e no final do período estudado, com valores de 0,28 ao nascimento, 0,21 aos 240 dias e 0,24 aos 550 dias de idade. As estimativas de herdabilidade materna foram maiores aos 160 dias de idade (0,10) que nas demais fases do crescimento. As correlações genéticas variaram de moderadas a altas, diminuindo conforme o aumento da distância entre as idades. Maior eficiência na seleção para peso pode ser obtida considerando os pesos pós-desmama, período em que as estimativas de variância genética e herdabilidade foram superiores.

Estimação de funções de covariância para características de crescimento da raça Tabapuã utilizando modelos de regressão aleatória

Foram utilizados 35.732 registros de peso do nascimento aos 660 dias de idade de 8.458 animais da raça Tabapuã para estimar funções de covariância utilizando modelos de regressão aleatória sobre polinômios de Legendre. Os modelos incluíram: como aleatórios, os efeitos genético aditivo direto, materno, de ambiente permanente de animal e materno; como fixos, os efeitos de grupo de contemporâneo; como covariáveis, a idade do animal à pesagem e a idade da vaca ao parto (linear e quadrática); e sobre a idade à pesagem, polinômio ortogonal de Legendre (regressão cúbica) foi considerado para modelar a curva média da população. O resíduo foi modelado considerando sete classes de variância e os modelos foram comparados pelos critérios de informação Bayesiano de Schwarz e Akaike. O melhor modelo apresentou ordens 4, 3, 6, 3 para os efeitos genético aditivo direto e materno, de ambiente permanente de animal e materno, respectivamente. As estimativas de covariância e herdabilidades, obtidas utilizando modelo bicaracter, e de regressão aleatória foram semelhantes. As estimativas de herdabilidade para o efeito genético aditivo direto, obtidas com o modelo de regressão aleatória, aumentaram do nascimento (0,15) aos 660 dias de idade (0,45). Maiores estimativas de herdabilidade materna foram obtidas para pesos medidos logo após o nascimento. As correlações genéticas variaram de moderadas a altas e diminuíram com o aumento da distância entre as pesagens. A seleção para maiores pesos em qualquer idade promove maior ganho de peso do nascimento aos 660 dias de idade.

Tendência Genética dos Efeitos Direto e Materno sobre os Pesos à Desmama e Pós-Desmama de Bovinos da Raça Tabapuã no Brasil

Os dados são relativos aos pesos de animais da raça Tabapuã, nascidos no período de 1959 a 1996 em várias regiões brasileiras. As observações foram analisadas com o objetivo de avaliar as mudanças genéticas aditivas diretas e maternas dos pesos padronizados para 205 (P205), 365 (P365) e 550 (P550) dias de idade. As estimativas dos componentes de (co) variância utilizadas no cálculo dos valores genéticos foram obtidas pelo método de máxima verossimilhança restrita livre de derivadas (REML), usando o aplicativo MTDFREML, cujo modelo continha os efeitos aleatórios aditivo direto, materno e de ambiente permanente, além dos efeitos fixos de grupo de contemporâneos (unidade da federação, fazenda, sexo, estação e ano de nascimento do animal) e a covariável idade da vaca ao parto (efeitos linear e quadrático). As tendências genéticas dos efeitos genéticos direto e materno foram estimadas pela regressão, ponderada, das médias anuais dos valores genéticos dos animais. As tendências genéticas dos efeitos direto no período estudado foram 0,134 ; 0,207, e 0,276 kg/ano, para P205, P365 e P550, respectivamente. Ainda para os três pesos, na mesma ordem, as estimativas das tendências genéticas maternas foram 0,019; -0,011; e -0,022 kg/ano. em virtude da variação genética existente, os resultados observados estão bem aquém das mudanças possíveis.

Genotype by environment interaction in different birth seasons for weight at 240, 365 and 450 days of age in Tabapuã cattle

The objective of this study was to evaluate the effect of genotype by environment interaction (GEI) on the weight of Tabapuã cattle at 240 (W240), 365 (W365) and 450 (W450) days of age. In total, 35,732 records of 8,458 Tabapuã animalswhich were born in the state of Bahia, Brazil, from 1975 to 2001, from 167 sires and 3,707 dams, were used. Two birth seasons were tested as for the environment effect: the dry (D) and rainy (R) ones. The covariance components were obtainedby a multiple-trait analysis using Bayesian inference, in which each trait was considered as being different in each season. Covariance components were estimated by software gibbs2f90. As for W240, the model was comprised of contemporary groups and cow age (in classes) as fixed effects; animal and maternal genetic additive, maternal permanent environmental and residual were considered as random effects. Concerning W365 and W450, the model included only the contemporary aged cow groups as fixed effects and the genetic additive and residual effects of the animal as the random ones. The GEI was assessed considering the genetic correlation, in which values below 0.80 indicated the presence of GEI. Regarding W365 and W450, the GEI was found in both seasons. As for post-weaning weight (W240), the effect of such interaction was not observed. ©2012 Sociedade Brasileira de Zootecnia.

Resposta sorológica de bezerras da raça Tabapuã recém-vacinadas com amostra B19 de Brucella abortus

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise de pedigree e estimativas de parâmetros genéticos para características de importância econômica em bovinos Tabapuã

Pós-graduação em Genética e Melhoramento Animal - FCAV

Propriedades físicas de latossolo e argissolo em função de práticas de manejo aplicadas na entrelinha da cultura da seringueira (Hevea brasiliensis)

Avaliou-se a influência das práticas de manejo (gradagem, Pueraria phaseoloides e roçadeira) nas entrelinhas da cultura da seringueira (Hevea brasiliensis), plantada em 1992, sobre a densidade do solo e a macro e microporosidade de solos do Planalto Paulista. Foram retiradas amostras nos anos de 1998 e 1999, após seis anos consecutivos de manejo, nas profundidades de 0-10, 10-20, 20-30 e 30-40 cm do Latossolo Vermelho distrófico textura argilosa A moderado caulinítico hipoférrico relevo plano e Argissolo Vermelho-Amarelo distrófico abrúptico, Tb, A moderado textura areia/média fase floresta tropical subperenifólia e relevo suave ondulado, localizados, respectivamente, nos municípios de Jaboticabal e Tabapuã (SP). Para cada profundidade, foram retiradas oito amostras por tratamento, por experimento, em cada ano, totalizando 394 amostras. Os atributos físicos do solo (densidade do solo e macroporosidade) foram avaliados de acordo com as práticas de manejo aplicadas. Dentre os sistemas de manejos aplicados na entrelinha da seringueira, a roçadeira foi o que provocou maior compactação do solo em todas as profundidades, evidenciada pelos elevados valores de densidade do solo e reduzida macroporosidade, sendo este efeito mais pronunciado nas camadas superficiais dos solos.

Modelagem geoestatística das incertezas da distribuição espacial do fósforo disponível no solo, em área de cana-de-açúcar

O fósforo destaca-se como um dos nutrientes limitantes ao desenvolvimento da cultura da cana-de-açúcar em solos brasileiros. Esse elemento apresenta grande variabilidade espacial, coordenada pelos atributos que regem as reações de adsorção e dessorção. Estimativas espaciais são conduzidas por meio de interpolações geoestatísticas para a caracterização dessa variabilidade. No entanto, tais estimativas apresentam incertezas inerentes ao procedimento que estão associadas à estrutura de variabilidade do atributo em estudo e à configuração amostral da área. Dessa forma, avaliar a incerteza das predições associada à distribuição espacial do fósforo disponível (Plábil) é importante para otimizar o uso dos fertilizantes fosfatados. O objetivo deste estudo foi avaliar o desempenho da simulação sequencial gaussiana (SSG) e da krigagem ordinária (KO) na modelagem da incerteza das predições do fósforo disponível. Uma malha amostral com 626 pontos foi instalada em uma área experimental de 200 hectares de cana-de-açúcar no município de Tabapuã, São Paulo. Foram geradas 200 realizações por meio do algoritmo da SSG. As realizações da SSG reproduziram as estatísticas e a distribuição dos dados amostrais. A estatística G (0,81) indicou boa proximidade entre as frações dos valores simulados e as dos observados. As realizações da SSG preservaram a variabilidade espacial do Plábil, sem o efeito de suavização obtido pelo mapa da KO. A acurácia na reprodução do variograma dos dados amostrais, obtida pelas realizações da SSG foi, em média, 240 vezes maior que obtida por meio da KO. O mapa de incertezas, obtido por meio da KO, apresentou menor variação na área de estudo do que por SSG. Dessa forma, a avaliação das incertezas, pela SSG, evidenciou-se mais informativa, podendo ser utilizada para definir e delimitar, de forma mais precisa, as áreas de manejo específico.

SPATIAL UNCERTAINTY OF NUTRIENT LOSS BY EROSION IN SUGARCANE HARVESTING SCENARIOS

The assessment of spatial uncertainty in the prediction of nutrient losses by erosion associated with landscape models is an important tool for soil conservation planning. The purpose of this study was to evaluate the spatial and local uncertainty in predicting depletion rates of soil nutrients (P, K, Ca, and Mg) by soil erosion from green and burnt sugarcane harvesting scenarios, using sequential Gaussian simulation (SGS). A regular grid with equidistant intervals of 50 m (626 points) was established in the 200-ha study area, in Tabapuã, São Paulo, Brazil. The rate of soil depletion (SD) was calculated from the relation between the nutrient concentration in the sediments and the chemical properties in the original soil for all grid points. The data were subjected to descriptive statistical and geostatistical analysis. The mean SD rate for all nutrients was higher in the slash-and-burn than the green cane harvest scenario (Student’s t-test, p<0.05). In both scenarios, nutrient loss followed the order: Ca>Mg>K>P. The SD rate was highest in areas with greater slope. Lower uncertainties were associated to the areas with higher SD and steeper slopes. Spatial uncertainties were highest for areas of transition between concave and convex landforms.

Predição de desempenho com uso de sistemas de exigências nutricionais e características de carcaça de tourinhos zebuínos

O objetivo deste trabalho foi avaliar a predição de desempenho pelos sistemas CNCPS 5.0, NRC 2000 e BR-CORTE, bem como o ganho de peso diário (GPD) e as características de carcaça em tourinhos zebuínos. Foram utilizados 44 tourinhos: 19 Nelore PO, 7 Nelore LA, 10 Tabapuã PO e 8 Guzerá PO, com peso corporal médio inicial de 266, 236, 222 e 219 kg e idade inicial média de 9, 10, 8 e 8 meses, respectivamente. O período experimental foi de 112 dias. O consumo alimentar individual foi obtido com o uso de indicadores (LIPE, óxido crômico e FDAi). Informações sobre área de olho de lombo (AOL), espessura de gordura subcutânea (EGS) e espessura de gordura na garupa (P8) foram obtidas por ultrassonografia. Para comparação do GPD predito e observado, foi utilizada análise de regressão linear. Os valores de GPD, AOL, EGS e P8 não diferiram entre os grupos genéticos. Para os sistemas CNCPS 5.0 e NRC 2000, os menores valores de GPD preditos foram estimados com base na disponibilidade de energia; para BR-CORTE, foram baseados na disponibilidade de proteína. Os sistemas NRC 2000, CNCPS 5.0 e BR-CORTE superestimaram o GPD, e não se mostraram adequados para predição do desempenho de tourinhos zebuínos.

Modelagem e quantificação da incerteza espacial do potássio disponível no solo por simulações estocásticas

O objetivo deste trabalho foi avaliar o desempenho da simulação sequencial gaussiana (SSG) e da simulação sequencial indicatriz (SSI) na modelagem da incerteza das predições do K disponível em área de cana-de-açúcar, e comparar as simulações com o método já consagrado de krigagem ordinária (KO). Uma malha amostral com 626 pontos foi instalada em área de 200 ha, no Município de Tabapuã, em São Paulo. As simulações reproduziram a variabilidade dos dados amostrais de K disponível, enquanto a KO superestimou os baixos teores de K e subestimou os altos. O mapa de desvio-padrão obtido a partir da KO mostrou menor variação ao longo da área de estudo, quando comparado aos mapas obtidos a partir das simulações. A SSI obteve acurácia 22% superior à obtida pela SSG, na modelagem da função de distribuição condicional do K. As simulações apresentam maior eficiência que a KO para modelar incerteza na distribuição espacial do K. A SSI apresenta melhor desempenho que a SSG na estimativa dos teores de K disponível, em área de cana-de-açúcar.

Alterações na estabilidade de agregados de latossolo e argissolo em função do manejo, na entrelinha da seringueira (Hevea brasiliensis)

Realizou-se um experimento com o propósito de avaliar o efeito de diferentes sistemas de manejo aplicados na entrelinha da seringueira (Hevea brasiliensis), na estabilidade de agregados em água de um Latossolo Vermelho, textura argilosa do Município de Jaboticabal, SP, e um Argissolo Vermelho-Amarelo, textura areia/média do Município de Tabapuã, SP, cujos tratamentos consistiram em: roçadeira, adubo-verde perene (Pueraria phaseoloides) e gradagem. As amostras de solo foram retiradas após sete anos de aplicação desses manejos nas camadas de 0-0,1; 0,1-0,2; 0,2-0,3; e 0,3-,4 m. A estabilidade de agregados foi obtida em peneira com classes de tamanhos de 8-4; 4-2; 2-1; 1-0,5; 0,5-0,125; e <0,125 mm. A Pueraria phaseoloides propiciou diferença significativa em relação aos demais sistemas de manejo na camada superficial (10 cm), quanto à distribuição de agregados, apresentando maior distribuição de agregados maiores no Latossolo Vermelho. O manejo com a grade dessa mesma camada, em ambos os solos, apresentou valores superiores de agregados de tamanho pequeno. A matéria orgânica dos solos apresentou relação direta e significativa com a estabilidade dos agregados maiores e relação inversa e significativa com a estabilidade de agregados menores.

Aplicação de microssatélites (STR) de bovinos em animais de raças zebuínas

O gado zebu (Bos indicus) é predominante em todo o mundo. Atualmente, no Brasil, o rebanho de bovinos e zebuínos é composto por aproximadamente 170 milhões de indivíduos, no qual, 80% são zebuínos ou animais de raças cruzadas com zebu. O principal objetivo deste estudo foi determinar a possibilidade da utilização de marcadores moleculares de bovinos na identificação genética de zebuínos. A confirmação desta possibilidade seria de extrema valia, já que a quantidade de informações sobre o genoma zebuíno é limitada e metodologias específicas de identificação de marcadores genéticos para raças zebuínas são raras. Assim, o DNA de 65 zebus, de 4 raças diferentes (Gir, Tabapuã, Nelore e Brahman), foi extraído a partir de sangue aplicado em papel filtro, utilizando o “kit” DNA IQTM SYSTEM (Promega®). Após a extração de DNA, foi feito um PCR Multiplex utilizando o “kit” StockMarks® (Applied Biosystems®), o qual amplifica um total de 11 loci (BM1824, BM2113, ETH10, ETH225, ETH3, INRA023, SPS115, TGLA122, TGLA126, TGLA227 e TGLA53). Os fragmentos gerados pela PCR foram analisados por eletroforese capilar utilizando o analisador genético ABI Prism 3100 (Applied Biosystems®) e o programa GeneScan® v.3.7 (Applied Biosystems®). A freqüência dos marcadores e os índices de variabilidade genética foram calculados utilizando o programa Microsatellite Toolkit v.3.1 e o equilíbrio de Hardy-Weinberg foi determinado através do programa GENEPOP v.3.4. Através destas técnicas foi possível demonstrar que os marcadores moleculares utilizados para identificação genética de bovinos podem também ser utilizados para a realização da técnica de identificação genética de zebuínos Além disso, foi possível demonstrar a freqüência alélica na população de zebuínos analisada como um todo, assim como foi determinada a freqüência alélica para os diferentes marcadores moleculares dentro das 4 raças estudas. Considerando o total de 100 diferentes alelos identificados entre os 11 STR (Short Tandem Repeats) analisados, foi possível observar que a heterozigosidade esperada (He) foi relativamente alta na população analisada, assim como na análise por raça. Estes dados indicam alto grau de diversidade genética nas diferentes raças. Neste sentido, a raça Brahman apresentou a menor diversidade (0,60) e a raça Tabapuã, a maior (0,69). Contudo, estudos complementares são necessários, a fim de confirmar as freqüências alélicas estabelecidas, uma vez que os resultados encontrados no presente trabalho referem-se a um número relativamente pequeno de zebuínos, quando comparados com o total de zebus encontrados no rebanho brasileiro.