Análise Bayesiana da área de olho do lombo e da espessura de gordura obtidas por ultrassom e suas associações com outras características de importância econômica na raça Nelore

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

Análise da curva de crescimento de bovinos da raça Nelore utilizando funções não-lineares em análises Bayesianas: Selma Forni. -

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Modelos de antedependência para avaliação genética de características longitudinais

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

Estudo da interação genótipo x ambiente sobre a produção de leite em rebanhos da raça Parodo-Suiço no Brasil, utilizando inferência Bayesiana

Verifica a presença de heterogeneidade de variâncias sobre a produção de leite na primeira lactação de fêmeas da raça Pardo-Suíço e, seu impacto sobre a avaliação genética dos reprodutores, utilizando a inferência Bayesiana por meio de amostrador de Gibbs, foram utilizados 2981 registros referentes às produções de leite e idade da vaca ao parto, em primeiras lactações de vacas da raça Pardo-Suíço, distribuídos em 62 rebanhos. Os registros foram provenientes do serviço de controle leiteiro da Associação Brasileira de Criadores de Gado Pardo Suíço, com os partos ocorridos entre os anos de 1980 a 2002. Foram estabelecidas duas classes de desvio-padrão fenotípico para produção de leite. Posteriormente, os dados foram analisados desconsiderando e considerando as classes de desvio-padrão. As médias observadas e desvio-padrão para produção de leite nas classes de alto e baixo desviopadrão e em análise geral foram iguais a 5802,02 ± 1929,96, 4844,37 ± 1592,99, 5373,47 ± 1849,13, respectivamente. As médias posteriores para os componentes de variâncias foram maiores na classe de alto desvio-padrão. A herdabilidade obtida na classe de alto desviopadrão foi próxima do valor observado na análise geral e inferior ao valor encontrado na classe de baixo desvio-padrão fenotípico. A correlação genética para a produção de leite entre as classes de desvio-padrão foi igual a 0,48. As correlações de Pearson e Spearman entre os valores genéticos para a produção de leite obtidos na análise geral, com os valores obtidos em cada classe de desvio-padrão foram todas maiores que 0,80, quando se considerou todos os reprodutores. Porém, refinando a amostra de reprodutores verifica-se que, as correlações diminuem em magnitude. Existindo uma maior variabilidade nos rebanhos presentes na classe de alto desvio-padrão e, o impacto dessa heterogeneidade de variância sobre a avaliação genética de reprodutores, é pequeno, pois a fonte principal dessa heterogeneidade é decorrente de fatores genéticos confirmando a presença de heterogeneidade de variâncias.

Estudo genético quantitativo das características andrológicas de touros jovens da raça Nelore

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

Genetic Parameters for Productive and Reproductive Traits for Milk Buffalo in Brazil

Knowing the genetic parameters of productive and reproductive traits in milking buffaloes is essential for planning and implementing of a program genetic selection. In Brazil, this information is still scarce. The objective of this study was to verify the existence of genetic variability in milk yield of buffaloes and their constituents, and reproductive traits for the possibility of application of the selection. A total of 9,318 lactations records from 3,061 cows were used to estimate heritabilities for milk yield (MY), fat percentage (%F), protein percentage (%P), length of lactation (LL), age of first calving (AFC) and calving interval (CI) and the genetic correlations among traits MY, %F and %P. The (co) variance components were estimated using multiple-trait analysis by Bayesian inference method, applying an animal model, through Gibbs sampling. The model included the fixed effects of contemporary groups (herd-year and calving season), number of milking (2 levels), and age of cow at calving as (co) variable (quadratic and linear effect). The additive genetic, permanent environmental, and residual effects were included as random effects in the model. Estimated heritability values for MY, % F, % P, LL, AFC and CI were 0.24, 0.34, 0.40, 0.09, 0.16 and 0.05, respectively. The genetic correlation estimates among MY and % F, MY and % P and % F and % P were -0.29, -0.18 and 0.25, respectively. The production of milk and its constituents showed enough genetic variation to respond to a selection program. Negative estimates of genetic correlations between milk production and its components suggest that selection entails a reduction in the other.

Multiple-trait genomic evaluation for milk yield and milk quality traits using genomic and phenotypic data in buffalo in Brazil

The use of markers distributed all long the genome may increase the accuracy of the predicted additive genetic value of young animals that are candidates to be selected as reproducers. In commercial herds, due to the cost of genotyping, only some animals are genotyped and procedures, divided in two or three steps, are done in order to include these genomic data in genetic evaluation. However, genomic evaluation may be calculated using one unified step that combines phenotypic data, pedigree and genomics. The aim of the study was to compare a multiple-trait model using only pedigree information with another using pedigree and genomic data. In this study, 9,318 lactations from 3061 buffaloes were used, 384 buffaloes were genotyped using a Illumina bovine chip (Illumina Infinium (R) bovineHD BeadChip). Seven traits were analyzed milk yield (MY), fat yield (FY), protein yield (PY), lactose yield (LY), fat percentage (F%), protein percentage (P%) and somatic cell score (SCSt). Two analyses were done: one using phenotypic and pedigree information (matrix A) and in the other using a matrix based in pedigree and genomic information (one step, matrix H). The (co) variance components were estimated using multiple-trait analysis by Bayesian inference method, applying an animal model, through Gibbs sampling. The model included the fixed effects of contemporary groups (herd-year-calving season), number of milking (2 levels), and age of buffalo at calving as (co) variable (quadratic and linear effect). The additive genetic, permanent environmental, and residual effects were included as random effects in the model. The heritability estimates using matrix A were 0.25, 0.22, 0.26, 0.17, 0.37, 0.42 and 0.26 and using matrix H were 0.25, 0.24, 0.26, 0.18, 0.38, 0.46 and 0.26 for MY, FY, PY, LY, % F, % P and SCCt, respectively. The estimates of the additive genetic effect for the traits were similar in both analyses, but the accuracy were bigger using matrix H (superior to 15% for traits studied). The heritability estimates were moderated indicating genetic gain under selection. The use of genomic information in the analyses increases the accuracy. It permits a better estimation of the additive genetic value of the animals.

Genetic parameters for earnings in Quarter Horse

In this study, we estimated the heritability (h(2)) of earnings in the Quarter Horse in order to evaluate the inclusion of this trait in breeding programs. Records from 14,754 races of 2443 horses from 1978-2009 were provided by Sorocaba Hippodrome, Sao Paulo, Brazil. All ancestors of the registered horses were included in the pedigree file until the 4th generation. Log-transformed performance measures (LPM) were analyzed for animals aged 2, 3, and 4 years and during their entire career. The h(2) estimates were obtained using a multi-trait model and Gibbs sampling that included the effects of sex, year of race, and animal in all analyses. Five analyses were performed: 1 in which LPM was divided by the number of prizes, 1 in which LPM was divided by the number of race starts, and 3 analyses that included the number of prizes, number of race starts, and both (LPM_cNPS) as covariates. Analysis was performed with and without inclusion of the maternal effect. Models were compared based on the deviance information criterion and LPM_cNPS including maternal effects was found to be the best model. The h(2) estimates and standard deviation obtained using model LPM_cNPS were 0.19 +/- 0.08, 0.21 +/- 0.08, 0.22 +/- 0.09, and 0.21 +/- 0.07 for earnings at 2, 3, and 4 years of age and total career, respectively. Our analyses indicate that earnings are subject to selection and can be included in breeding programs to improve the racing performance of Quarter Horses.

Finishing precocity visual score and genetic associations with growth traits in Angus beef cattle

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

Produção de leite, duração da lactação e intervalo de partos em búfalas mestiças Murrah

This study aimed to evaluate the potential for milk production (MP), lactation length (LL) and calving interval (CI), analyze the environmental component affecting these traits, and to estimate the heritability and repeatability for milk production in crossbreds of Murrah buffalo cows in the state of Alagoas, Brazil. Data was composed of 487 observations of MP from 136 lactations recorded between the years of 2000 and 2010. In the analysis of variance for PL, the fixed effects were season (1- October to March, 2 -April to September) and year of the beginning of lactation, calving order and the LL (covariate). For the analysis of LL only the fixed effect of year of the beginning of lactation was included, and finally for the CI analysis, year of the beginning of lactation and calving order. The estimates of covariance were obtained using unicharacteristic analysis by Bayesian inference method, applyingan animal model, through Gibbs sampling. The additive genetic, permanent environment and residual effects were included as random effects in the model. The averages (sd) of MP, LL and CI were 2,218.03 kg (408.18), 282.59 days (39.48) e 422.49 days (91.05), respectively. All the effects included in the models were important (P<0.01). The estimates of heritability and repeatability for PL were 0.29 and 0.69, respectively. The results suggest that there is a moderate genetic variability among individuals for PL, indicating the possibility to obtain gain using selection.

Genetic associations between temperament and performance traits in Nellore beef cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Análise bayesiana para características de carcaça avaliadas por ultrassonografia de bovinos da raça Nelore Mocho, criados em bioma Cerrado

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Parâmetros genéticos para características indicadoras de eficiência reprodutiva e produtiva de ovinos da raça Santa Inês

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

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera

Abstract Background In honeybees, differential feeding of female larvae promotes the occurrence of two different phenotypes, a queen and a worker, from identical genotypes, through incremental alterations, which affect general growth, and character state alterations that result in the presence or absence of specific structures. Although previous studies revealed a link between incremental alterations and differential expression of physiometabolic genes, the molecular changes accompanying character state alterations remain unknown. Results By using cDNA microarray analyses of >6,000 Apis mellifera ESTs, we found 240 differentially expressed genes (DEGs) between developing queens and workers. Many genes recorded as up-regulated in prospective workers appear to be unique to A. mellifera, suggesting that the workers' developmental pathway involves the participation of novel genes. Workers up-regulate more developmental genes than queens, whereas queens up-regulate a greater proportion of physiometabolic genes, including genes coding for metabolic enzymes and genes whose products are known to regulate the rate of mass-transforming processes and the general growth of the organism (e.g., tor). Many DEGs are likely to be involved in processes favoring the development of caste-biased structures, like brain, legs and ovaries, as well as genes that code for cytoskeleton constituents. Treatment of developing worker larvae with juvenile hormone (JH) revealed 52 JH responsive genes, specifically during the critical period of caste development. Using Gibbs sampling and Expectation Maximization algorithms, we discovered eight overrepresented cis-elements from four gene groups. Graph theory and complex networks concepts were adopted to attain powerful graphical representations of the interrelation between cis-elements and genes and objectively quantify the degree of relationship between these entities. Conclusion We suggest that clusters of functionally related DEGs are co-regulated during caste development in honeybees. This network of interactions is activated by nutrition-driven stimuli in early larval stages. Our data are consistent with the hypothesis that JH is a key component of the developmental determination of queen-like characters. Finally, we propose a conceptual model of caste differentiation in A. mellifera based on gene-regulatory networks.