Models for genetic evaluation of Nelore cattle mature body weight

Records of 18,770 Nelore animals, born from 1975 to 2002, in 8 herds participating in the Nelore Cattle Breeding Program, were analyzed to estimate genetic parameters for mature BW. The mature BW were analyzed as a single BW taken closest to 4.5 yr of age for each cow in the data file, considering BW starting from 2 (W2Y_S), 3 (W3Y_S), or 4 (W4Y_S) yr of age or as repeated records, including all BW starting from 2 (W2Y_R), 3 (W3Y_R), or 4 (W4Y_R) yr of age. The variance components were estimated by restricted maximum likelihood, fitting univariate and bivariate animal models, including weaning weight. The heritability estimates were 0.29, 0.34, 0.36, 0.41, 0.44, and 0.46 for W2Y_S, W3Y_S, W4Y_S, W2Y_R, W3Y_R, and W4Y_R, respectively. The repeatability estimates for W2Y_R, W3Y_R, and W4Y_R were 0.59, 0.64, and 0.72, respectively. Larger accuracy values associated with the EBV were obtained in the repeated records models. The results indicated the bivariate repeated records model as the most appropriate for analyzing mature BW.

Genetic associations between age at first calving and heifer body weight and scrotal circumference in Nelore cattle

P>Age at first calving (AFC) measures the entry of heifers into the beef cattle production system. This trait can be used as a selection criterion for earlier reproductive performance. Using data from Nelore cattle participating in the `Program for Genetic Improvement of the Nelore Breed` (PMGRN-Nelore Brazil), bi-trait analyses were performed using the restricted maximum likelihood method, based on an AFC animal model and the following traits: female body weight adjusted to 365 (BW365) and 450 (BW450) days of age, and male scrotal circumference adjusted to 365 (SC365), 450 (SC450), 550 (SC550) and 730 (SC730) days of age. The heritability estimates for AFC ranged from 0.02 +/- 0.02 to 0.04 +/- 0.02. The estimates of additive direct heritabilities (with standard error) for BW365, BW450, SC365, SC450, SC550 and SC730 were 0.36 +/- 0.07, 0.38 +/- 0.07, 0.48 +/- 0.07, 0.65 +/- 0.07, 0.64 +/- 0.07 and 0.42 +/- 0.07, respectively, and the genetic correlations with AFC were -0.38, -0.33, 0.10, -0.13, -0.13 and 0.06, respectively. In the herds studied, selection for SC365, SC450, SC550 or SC730 should not cause genetic changes in AFC. Selection based on BW365 or BW450 would favor smaller AFC breeding values. However, the low magnitude of direct heritability estimates for AFC in these farms indicates that changes in phenotypical expression depend mostly on non-genetic factors.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

Genetic associations between accumulated productivity, and reproductive and growth traits in Nelore cattle

The total meat yield in a beef cattle production cycle is economically very important and depends on the number of calves born per year or birth season, being directly related to reproductive potential. Accumulated Productivity (ACP) is an index that expresses a cow`s capacity to give birth regularly at a young age and to wean animals of greater body weight. Using data from cattle participating in the ""Program for Genetic Improvement of the Nelore Breed"" (PMGRN - Nelore Brasil), bi-trait analyses were performed using the Restricted Maximum Likelihood method based on an ACP animal model and the following traits: age at first calving (AFC), female body weight adjusted for 365 (BW365) and 450 (BW450) days of age, and male scrotal circumference adjusted for 365 (SC365), 450 (SC450), 550 (SC550) and 730 (SC730) days of age. Median estimated ACP heritability was 0.19 and the genetic correlations with AFC, BW365, BW450, SC365, SC450, SC550 and SC730 were 0.33, 0.70, 0.65, 0.08, 0.07, 0.12 and 0.16, respectively. ACP increased and AFC decreased over time, revealing that the selection criteria genetically improved these traits. Selection based on ACP appears to favor the heaviest females at 365 and 450 days of age who showed better reproductive performance as regards AFC. Scrotal circumference was not genetically associated with ACP. (C) 2007 Elsevier B.V. All rights reserved.

Genetic and environmental factors affecting ultrasound measures of longissimus muscle area and backfat thickness in Nelore cattle

The aim of the present study was to evaluate the genetic and environmental factors affecting records of longissimus muscle area (LMA) and back fat thickness (BF) obtained between the 12th and 13th ribs, and rump fat thickness (RF) between the hook and pin bones, measured by real-time ultrasound in Nelore cattle. Also, weight records of 22,778 animals born from 1998 to 2003, in ten farms across six Brazilian states were used. Carcass traits as measured by ultrasound of the live animal were recorded from 2002 to 2004 in 2590 males and females with ages varying from 450 to 599 days. Fixed models including farm, year and season of birth, sex and type of feed effects, and the covariates age of dam (AOD) and age of animal at measurement were used to study the effect of environmental factors on these traits. The genetic parameters for LMA, BF and RF were estimated with two and three-trait animal models with 120-day weights using a restricted maximum likelihood method. All environmental effects significantly affected carcass traits, with the exception of year of birth for BF and RF and AOD for LMA. The heritability estimates for LMA, BF and RF were 0.35, 0.51 and 0.39, respectively. Standard errors obtained in one-trait analyses were from 0.07 to 0.09. Genetic correlation estimates between LMA and the two traits of subcutaneous fat were low (close to zero) and 0.74 between BF and RF, indicating that the selection for LMA should not cause antagonism in the genetic improvement of subcutaneous fat measured by real-time ultrasound. (C) 2007 Elsevier B.V. All fights reserved.

Multi-trait and random regression mature weight heritability and breeding value estimates in Nelore cattle

Mature weight breeding values were estimated using a multi-trait animal model (MM) and a random regression animal model (RRM). Data consisted of 82 064 weight records from 8 145 animals, recorded from birth to eight years of age. Weights at standard ages were considered in the MM. All models included contemporary groups as fixed effects, and age of dam (linear and quadratic effects) and animal age as covariates. In the RRM, mean trends were modelled through a cubic regression on orthogonal polynomials of animal age and genetic maternal and direct and maternal permanent environmental effects were also included as random. Legendre polynomials of orders 4, 3, 6 and 3 were used for animal and maternal genetic and permanent environmental effects, respectively, considering five classes of residual variances. Mature weight (five years) direct heritability estimates were 0.35 (MM) and 0.38 (RRM). Rank correlation between sires' breeding values estimated by MM and RRM was 0.82. However, selecting the top 2% (12) or 10% (62) of the young sires based on the MM predicted breeding values, respectively 71% and 80% of the same sires would be selected if RRM estimates were used instead. The RRM modelled the changes in the (co) variances with age adequately and larger breeding value accuracies can be expected using this model.

Immunological responses and cytokine gene expression analysis to Cooperia punctata infections in resistant and susceptible Nelore cattle

Cellular and humoral immune response, as well as cytokine gene expression, was assessed in Nelore cattle with different degrees of resistance to Cooperia punctata natural infection. One hundred cattle (male, weaned, 11-12 months old), kept together on pasture, were evaluated. Faecal and blood samples were collected for parasitological and immunological assays. Based on nematode faecal egg counts (FEC) and worm burden, the seven most resistant and the eight most susceptible animals were selected. Tissue samples of the small intestine were collected for histological quantification of inflammatory cells and analysis of cytokine gene expression (IL-2, IL-4, IL-8, IL-1 2p35, IL-13, TNF-alpha, IFN-gamma, MCP-1, MCP-2, and MUC- 1) using real-time RT-PCR. Mucus samples were also collected for IgA levels determination. Serum IgG1 mean levels against C. punctata antigens were higher in the resistant group, but significant differences between groups were only observed 14 days after the beginning of the experiment against infective larvae (1-3) and 14 and 84 days against adult antigens. The resistant group also presented higher IgA levels against C. punctata (L3 and adult) antigens with significant difference 14 days after the beginning of the trial (P < 0.05). In the small-intestine mucosa, levels of IgA anti-L3 and anti-adult C. punctata were higher in the resistant group, compared with the susceptible group (P < 0.05). Gene expression of both T(H)2 cytokines (IL-4 and IL-13) in the resistant group and T(H)1 cytokines (IL-2, IL-1 2p35, IFN-gamma and MCP-1) in the susceptible group was up-regulated. Such results suggested that immune response to C. punctata was probably mediated by TH2 cytokines in the resistant group and by T(H)1 cytokines in the susceptible group. (C) 2008 Elsevier B.V. All rights reserved.

Comparison of different nonlinear functions to describe Nelore cattle growth

This work aims to compare different nonlinear functions for describing the growth curves of Nelore females. The growth curve parameters, their (co) variance components, and environmental and genetic effects were estimated jointly through a Bayesian hierarchical model. In the first stage of the hierarchy, 4 nonlinear functions were compared: Brody, Von Bertalanffy, Gompertz, and logistic. The analyses were carried out using 3 different data sets to check goodness of fit while having animals with few records. Three different assumptions about SD of fitting errors were considered: constancy throughout the trajectory, linear increasing until 3 yr of age and constancy thereafter, and variation following the nonlinear function applied in the first stage of the hierarchy. Comparisons of the overall goodness of fit were based on Akaike information criterion, the Bayesian information criterion, and the deviance information criterion. Goodness of fit at different points of the growth curve was compared applying the Gelfand`s check function. The posterior means of adult BW ranged from 531.78 to 586.89 kg. Greater estimates of adult BW were observed when the fitting error variance was considered constant along the trajectory. The models were not suitable to describe the SD of fitting errors at the beginning of the growth curve. All functions provided less accurate predictions at the beginning of growth, and predictions were more accurate after 48 mo of age. The prediction of adult BW using nonlinear functions can be accurate when growth curve parameters and their (co) variance components are estimated jointly. The hierarchical model used in the present study can be applied to the prediction of mature BW in herds in which a portion of the animals are culled before adult age. Gompertz, Von Bertalanffy, and Brody functions were adequate to establish mean growth patterns and to predict the adult BW of Nelore females. The Brody model was more accurate in predicting the birth weight of these animals and presented the best overall goodness of fit.

Genetic associations between carcass traits measured by real-time ultrasound and scrotal circumference and growth traits in Nelore cattle

The aim of the present study was to evaluate the genetic correlations among real-time ultrasound carcass, BW, and scrotal circumference (SC) traits in Nelore cattle. Carcass traits, measured by real-time ultrasound of the live animal, were recorded from 2002 to 2004 on 10 farms across 6 Brazilian states on 2,590 males and females ranging in age from 450 to 599 d. Ultrasound records of LM area (LMA) and backfat thickness (BF) were obtained from cross-sectional images between the 12th and 13th ribs, and rump fat thickness (RF) was measured between the hook and pin bones over the junction between gluteus medius and biceps femoris muscles. Also, BW (n = 22,778) and SC ( n = 5,695) were recorded on animals born between 1998 and 2003. The BW traits were 120, 210, 365, 450, and 550-d standardized BW (W120, W210, W365, W450, and W550), plus BW (WS) and hip height (HH) on the ultrasound scanning date. The SC traits were 365-, 450-, and 550-d standardized SC (SC365, SC450, and SC550). For the BW and SC traits, the database used was from the Nelore Breeding Program-Nelore Brazil. The genetic parameters were estimated with multivariate animal models and REML. Estimated genetic correlations between LMA and other traits were 0.06 (BF), -0.04 ( RF), 0.05 (HH), 0.58 (WS), 0.53 (W120), 0.62 (W210), 0.67 (W365), 0.64 ( W450 and W550), 0.28 (SC365), 0.24 (SC450), and 0.00 ( SC550). Estimated genetic correlations between BF and with other traits were 0.74 ( RF), -0.32 (HH), 0.19 (WS), -0.03 (W120), -0.10 (W210), 0.04 (W365), 0.01 (W450), 0.06 ( W550), 0.17 (SC365 and SC450), and -0.19 (SC550). Estimated genetic correlations between RF and other traits were -0.41 (HH), -0.09 (WS), -0.13 ( W120), -0.09 ( W210), -0.01 ( W365), 0.02 (W450), 0.03 (W550), 0.05 ( SC365), 0.11 ( SC450), and -0.18 (SC550). These estimates indicate that selection for carcass traits measured by real-time ultrasound should not cause antagonism in the genetic improvement of SC and BW traits. Also, selection to increase HH might decrease subcutaneous fat as correlated response. Therefore, to obtain animals suited to specific tropical production systems, carcass, BW, and SC traits should be considered in selection programs.

Covariance functions for body weight from birth to maturity in Nellore cows

The objective of this study was to estimate (co)variance functions using random regression models on Legendre polynomials for the analysis of repeated measures of BW from birth to adult age. A total of 82,064 records from 8,145 females were analyzed. Different models were compared. The models included additive direct and maternal effects, and animal and maternal permanent environmental effects as random terms. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of animal age (cubic regression) were considered as random co-variables. Eight models with polynomials of third to sixth order were used to describe additive direct and maternal effects, and animal and maternal permanent environmental effects. Residual effects were modeled using 1 (i.e., assuming homogeneity of variances across all ages) or 5 age classes. The model with 5 classes was the best to describe the trajectory of residuals along the growth curve. The model including fourth- and sixth-order polynomials for additive direct and animal permanent environmental effects, respectively, and third-order polynomials for maternal genetic and maternal permanent environmental effects were the best. Estimates of (co) variance obtained with the multi-trait and random regression models were similar. Direct heritability estimates obtained with the random regression models followed a trend similar to that obtained with the multi-trait model. The largest estimates of maternal heritability were those of BW taken close to 240 d of age. In general, estimates of correlation between BW from birth to 8 yr of age decreased with increasing distance between ages.

Validation of body condition score as a predictor of subcutaneous fat in Nelore (Bos indicus) cows

The aim of the present study was to determine the relationship among body weight (BW), body condition score (BCS) and rump fat thickness (RFAT) measured by ultrasonography, and validate the relationship between BCS and RFAT over the time. Two hundred sixty and six Nelore cows had their BW, BCS and RFAT evaluated at five different moments during the production cycle: M1) weaning: M2) parturition, M3) 42 days post-partum; M4) 82 days postpartum and M5) 112 days post-partum. A BCS value was attributed for each cow following a I to 5 points scale. Ultrasonographic images for RFAT measurement were obtained using a 3.5 MHz linear transducer. Images were immediately analyzed as soon as they were formed and frozen. Body condition scores and ultrasound measurements were collected on the same day by a single trained technician. The relationship between BCS and RFAT values was investigated by regression models. The analysis of similarity among the five obtained models was performed using the proc MIXED from SAS and the correlations among variables were analyzed with proc CORR from SAS. The BCS was able to predict RFAT in Nelore cows in all different moments evaluated. Also, it was shown that BCS presented high correlation (r=0.82 to 0.93) and relationship (R(2) = 0.73 to 0.92) with RFAT. However, both BCS and RFAT showed low correlation (r=0.37 to 0.50) and relationship (R(2) = 0.13 to 0.25) with BW. The BCS classification by visual method using a 1 to 5 point scale, was able to predict the RFAT in Nelore cows over the time. (C) 2008 Elsevier B.V. All rights reserved.

Causas de variação nos preços de bovinos nelore elite no Brasil

Este trabalho teve por objetivo estudar as causas de variação nos preços de bovinos da raça nelore pertencentes a rebanhos de seleção, os quais foram comercializados em leilões, para verificar as influências das avaliações genéticas e dos julgamentos de exterior sobre esses preços. Para tanto, foram computados os preços de venda de 426 bovinos da referida raça em 12 leilões ocorridos em diversas localidades brasileiras (regiões Centro-Oeste, Norte e Sudeste), entre os anos de 2002 e 2005. O valor médio foi de R$ 3.325,49, sendo o mínimo de R$ 1.400,00 e o máximo de R$ 10.500,00. Esses dados foram digitados juntamente com outras informações que eram apresentadas nos catálogos dos leilões. As informações registradas incluíram o sexo de cada animal, o nome do leilão e as DEPs informadas nos catálogos. Além da avaliação da influência das informações dos catálogos, também foi avaliada a influência das informações dos reprodutores, pais dos animais vendidos nos leilões, envolvendo suas DEPs publicadas em um sumário de reprodutores da raça e as pontuações de suas progênies em julgamentos. Os métodos estatísticos aplicados foram análises de variâncias e análises de agrupamento (método K-médias). Como resultado, foi observado que animais com superioridade genética em características relacionadas a desempenho ponderal, considerando-se os efeitos diretos e maternos, foram valorizados ao serem comercializados nos leilões. Em contra-partida, a pontuação dos reprodutores nos julgamentos não teve influência significativa sobre os preços médios de venda de suas progênies nos leilões.

Additive genetic relationship of longevity with fertility and production traits in Nellore cattle based on bivariate models

Survival or longevity is an economically important trait in beef cattle. The main inconvenience for its inclusion in selection criteria is delayed recording of phenotypic data and the high computational demand for including survival in proportional hazard models. Thus, identification of a longevity-correlated trait that could be recorded early in life would be very useful for selection purposes. We estimated the genetic relationship of survival with productive and reproductive traits in Nellore cattle, including weaning weight (WW), post-weaning growth (PWG), muscularity (MUSC), scrotal circumference at 18 months (SC18), and heifer pregnancy (HP). Survival was measured in discrete time intervals and modeled through a sequential threshold model. Five independent bivariate Bayesian analyses were performed, accounting for cow survival and the five productive and reproductive traits. Posterior mean estimates for heritability (standard deviation in parentheses) were 0.55 (0.01) for WW, 0.25 (0.01) for PWG, 0.23 (0.01) for MUSC, and 0.48 (0.01) for SC18. The posterior mean estimates (95% confidence interval in parentheses) for the genetic correlation with survival were 0.16 (0.13-0.19), 0.30 (0.25-0.34), 0.31 (0.25-0.36), 0.07 (0.02-0.12), and 0.82 (0.78-0.86) for WW, PWG, MUSC, SC18, and HP, respectively. Based on the high genetic correlation and heritability (0.54) posterior mean estimates for HP, the expected progeny difference for HP can be used to select bulls for longevity, as well as for post-weaning gain and muscle score.

Genetic analysis of average annual productivity of Nellore breeding cows (COWPROD)

With the aim of estimating the coefficient of heritability of average annual productivity of Nellore cows (COWPROD), a data set from 24,855 animals with known pedigree was analyzed. COWPROD is defined as the amount (in kilograms) of weaned calves produced yearly by one cow during her remaining time in herd ignoring a fixed period of 365 days. COWPROD was calculated regarding three standards: a) based on the post-weaning weight from the calves ignoring any kind of adjustment (COWPROD_NAJ), b) adjusted weight for the fixed effects (COWPROD_AJFIX) and c) adjusted weight for the fixed effects and for the genetic merit of the sire (COWPROD_AJFIN). The obtained heritabilities were 0.15, 0.15 and 0.16 for COWPROD_NAJ, COWPROD_AJFIX and COWPROD_AJFIN, respectively. A complete set composed of 105,158 COWPROD records on 130,740 animals in pedigree was also analyzed for predicting the genetic merit of all animals in the data set and for the calculation of the genetic, phenotypic and residual trends. Ranking correlation was high for the adjusted and non-adjusted data, yet, for some of the animals, the difference among the genetic values was large. This would be an indication that it would be better to work always with the adjusted weaning weights. The genetic trend was positive, but was of small magnitude (0.26% of the trait average) and the residual trend was negative as a consequence of the large intensification of the production system, which has been occurring in the last years in the farms studied. The phenotypic trend was also negative and intermediate between the genetic and the residual ones.

mRNA profile of Nellore calves after primary infection with Haemonchus placei

Haemonchus parasites are responsible for many losses in animal production. However, few studies are available, especially of zebu cattle. In this study, we investigated mRNA differences of immune response genes in naive Nellore calves infected with Haemonchus placei, relating these differences to patterns of cellular infiltrate. Calves were infected with 15,000 H. placei 13 larvae and after 7 days lymph node and abomasum tissues were collected. IL-2, IL-4, IL-8, IL-12, IL-13, IFN-gamma, MCP-1, lysozyme, pepsinogen and TNF-alpha genes were evaluated by qPCR. Mast cells, eosinophils and globular leukocytes were counted by abomasum histology. In the infected group, IL-4, IL-13 and TNF-alpha were up-regulated in the abomasal lymph node. In the abomasum, IL-13 increased and TNF-alpha was down-regulated (p < 0.05). No differences were detected for mast cells and eosinophil counts in abomasal tissue (p > 0.05). We conclude that for this infection time, there was Th2 polarization but that cellular infiltrate in abomasal tissue takes longer to develop. (C) 2010 Elsevier B.V. All rights reserved.