Cumulative t-link threshold models for the genetic analysis of calving ease scores

In this study, a hierarchical threshold mixed model based on a cumulative t-link specification for the analysis of ordinal data or more, specifically, calving ease scores, was developed. The validation of this model and the Markov chain Monte Carlo (MCMC) algorithm was carried out on simulated data from normally and t4 (i.e. a t-distribution with four degrees of freedom) distributed populations using the deviance information criterion (DIC) and a pseudo Bayes factor (PBF) measure to validate recently proposed model choice criteria. The simulation study indicated that although inference on the degrees of freedom parameter is possible, MCMC mixing was problematic. Nevertheless, the DIC and PBF were validated to be satisfactory measures of model fit to data. A sire and maternal grandsire cumulative t-link model was applied to a calving ease dataset from 8847 Italian Piemontese first parity dams. The cumulative t-link model was shown to lead to posterior means of direct and maternal heritabilities (0.40 ± 0.06, 0.11 ± 0.04) and a direct maternal genetic correlation (-0.58 ± 0.15) that were not different from the corresponding posterior means of the heritabilities (0.42 ± 0.07, 0.14 ± 0.04) and the genetic correlation (-0.55 ± 0.14) inferred under the conventional cumulative probit link threshold model. Furthermore, the correlation (> 0.99) between posterior means of sire progeny merit from the two models suggested no meaningful rerankings. Nevertheless, the cumulative t-link model was decisively chosen as the better fitting model for this calving ease data using DIC and PBF

Relação da amostragem do solo com a variabilidade das estimativas de estoque de carbono.

Resumo: Com o objetivo de estabelecer protocolos para a amostragem de solos agrícolas com foco na determinação precisa do seu estoque de carbono (EC), informação fundamental para o estabelecimento da referência inicial (baseline) nos projetos e auditorias sobre sequestro de Carbono (C), o presente trabalho investigou uma questão central nessa amostragem: qual o tamanho da amostra (número de pontos de coleta) necessário para uma estimativa com precisão aceitável? Essa é uma indagação realmente importante, pois uma amostra muito grande pode significar desperdício de recursos, enquanto uma amostra muito pequena pode reduzir a utilidade dos resultados ou mesmo torná-los inúteis para o objetivo pretendido. Com base em um conjunto de dados amostrais de campo, oriundos de propriedades agrícolas localizadas em diferentes estados brasileiros, foi possível, utilizando a técnica de bootstrap, simular o erro amostral para diversos tamanhos de amostra. A coleta e análise de amostras estratificadas por profundidade, que têm seus estoques de carbono somados para compor o estoque final de 0 a 30 centímetros (EC30), se mostrou mais eficiente, no sentido de proporcionar o mesmo coeficiente de variação (CV = 10%) com menor tamanho de amostra (n=9), do que se a escolha do tamanho de amostra fosse baseada na variabilidade do estoque de carbono em cada camada separadamente, para depois somá-las

Uso de função de pedotransferência para estimar a densidade do solo em áreas de produção de grãos no Brasil.

A função de pedotransferência obtida com os dados do procarbono permite estimar com boa eficiência a densidade do solo e resulta em baixo erro no estoque de carbono

Use of linear mixed models for genetic evaluation of gestation length and birth weight allowing for heavy-tailed residual effects

<p>Abstract</p> <p>Background</p> <p>The distribution of residual effects in linear mixed models in animal breeding applications is typically assumed normal, which makes inferences vulnerable to outlier observations. In order to mute the impact of outliers, one option is to fit models with residuals having a heavy-tailed distribution. Here, a Student's-<it>t </it>model was considered for the distribution of the residuals with the degrees of freedom treated as unknown. Bayesian inference was used to investigate a bivariate Student's-<it>t </it>(BS<it>t</it>) model using Markov chain Monte Carlo methods in a simulation study and analysing field data for gestation length and birth weight permitted to study the practical implications of fitting heavy-tailed distributions for residuals in linear mixed models.</p> <p>Methods</p> <p>In the simulation study, bivariate residuals were generated using Student's-<it>t </it>distribution with 4 or 12 degrees of freedom, or a normal distribution. Sire models with bivariate Student's-<it>t </it>or normal residuals were fitted to each simulated dataset using a hierarchical Bayesian approach. For the field data, consisting of gestation length and birth weight records on 7,883 Italian Piemontese cattle, a sire-maternal grandsire model including fixed effects of sex-age of dam and uncorrelated random herd-year-season effects were fitted using a hierarchical Bayesian approach. Residuals were defined to follow bivariate normal or Student's-<it>t </it>distributions with unknown degrees of freedom.</p> <p>Results</p> <p>Posterior mean estimates of degrees of freedom parameters seemed to be accurate and unbiased in the simulation study. Estimates of sire and herd variances were similar, if not identical, across fitted models. In the field data, there was strong support based on predictive log-likelihood values for the Student's-<it>t </it>error model. Most of the posterior density for degrees of freedom was below 4. Posterior means of direct and maternal heritabilities for birth weight were smaller in the Student's-<it>t </it>model than those in the normal model. Re-rankings of sires were observed between heavy-tailed and normal models.</p> <p>Conclusions</p> <p>Reliable estimates of degrees of freedom were obtained in all simulated heavy-tailed and normal datasets. The predictive log-likelihood was able to distinguish the correct model among the models fitted to heavy-tailed datasets. There was no disadvantage of fitting a heavy-tailed model when the true model was normal. Predictive log-likelihood values indicated that heavy-tailed models with low degrees of freedom values fitted gestation length and birth weight data better than a model with normally distributed residuals.</p> <p>Heavy-tailed and normal models resulted in different estimates of direct and maternal heritabilities, and different sire rankings. Heavy-tailed models may be more appropriate for reliable estimation of genetic parameters from field data.</p

The dog as an animal model for DISH?

Diffuse idiopathic skeletal hyperostosis (DISH) is a systemic disorder of the axial and peripheral skeleton in humans and has incidentally been described in dogs. The aims of this retrospective radiographic cohort study were to determine the prevalence of DISH in an outpatient population of skeletally mature dogs and to investigate if dogs can be used as an animal model for DISH. The overall prevalence of canine DISH was 3.8% (78/2041). The prevalence of DISH increased with age and was more frequent in male dogs, similar to findings in human studies. In the Boxer breed the prevalence of DISH was 40.6% (28/69). Dog breeds represent closed gene pools with a high degree of familiar relationship and the high prevalence in the Boxer may be indicative of a genetic origin of DISH. It is concluded that the Boxer breed may serve as an animal model for DISH in humans

Analysis of alanine aminotransferase in various organs of soybean (Glycine max) and in dependence of different nitrogen fertilisers during hypoxic stress

Alanine aminotransferase (AlaAT) catalyses the reversible conversion of pyruvate and glutamate into alanine and oxoglutarate. In soybean, two subclasses were identified, each represented by two highly similar members. To investigate the role of AlaAT during hypoxic stress in soybean, changes in transcript level of both subclasses were analysed together with the enzyme activity and alanine content of the tissue. Moreover, the dependency of AlaAT activity and gene expression was investigated in relation to the source of nitrogen supplied to the plants. Using semi-quantitative PCR, GmAlaAT genes were determined to be highest expressed in roots and nodules. Under normal growth conditions, enzyme activity of AlaAT was detected in all organs tested, with lowest activity in the roots. Upon waterlogging-induced hypoxia, AlaAT activity increased strongly. Concomitantly, alanine accumulated. During re-oxygenation, AlaAT activity remained high, but the transcript level and the alanine content decreased. Our results show a role for AlaAT in the catabolism of alanine during the initial period of re-oxygenation following hypoxia. GmAlaAT also responded to nitrogen availability in the solution during waterlogging. Ammonium as nitrogen source induced both gene expression and enzyme activity of AlaAT more than when nitrate was supplied in the nutrient solution. The work presented here indicates that AlaAT might not only be important during hypoxia, but also during the recovery phase after waterlogging, when oxygen is available to the tissue again

Evaluating alternate models to estimate genetic parameters of calving traits in United Kingdom Holstein-Friesian dairy cattle

<p>Abstract</p> <p>Background</p> <p>The focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.</p> <p>Methods</p> <p>Data from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (−maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.</p> <p>Results and discussion</p> <p>On average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (−maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (−maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.</p> <p>Conclusions</p> <p>For the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (−maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.</p