Longitudinal genetic dynamics of weaning index and implications for cow-calf production efficiency

In beef cattle, the selection for higher weights at young ages has been questioned with the argument that this criterion may increase the adult weight of cows, resulting in higher costs. Therefore, selection criteria should be employed to increase weights at young ages with minimal impact on the adult weight of cows. Additionally, the relationship between measures of cow production efficiency and other well-established selection criteria in breeding programs remains poorly understood. The objective of this study was to longitudinally evaluate the relationship between the weaning index (WIndex) as a measure of efficiency and growth traits of the cows. Possible changes over time in WIndex due to selection applied for yearling weight (YW) were also investigated. The WIndex was proposed to maximize genetic response in the weaning weight of the calf while maintaining genetic gain in BW of the cow at zero. A random regression model was adopted to estimate correlations between WIndex, BW, hip height (HH), and body condition score (BCS) using records of Nelore cows from three lines. Genetic trends were calculated for the control line (NeC) and lines selected for greater YW (NeS and NeT). The age of 3 years was the most critical for the weaning efficiency of the cows. At this stage, young cows are still growing and wean lighter calves than their adult counterparts. The genetic correlation estimates between WIndex and BW (-0.58 to 0.04), HH (-0.05 to -0.34), and BCS (-0.51 to -0.17) were close to zero or negative. BW and HH were strongly correlated genetically across all ages (0.73-0.76). In general, HH exhibited a weak and negative genetic relationship with BCS. The genetic correlation between BW and BCS was stronger for advanced ages (0.45-0.68). In lines selected for YW, important increases in WIndex were observed. However, NeS has been selected since the 1980s until the present for YW, and thus, it showed a more pronounced trend of increasing BW and, consequently, a more modest trend of increasing WIndex compared to NeT. In contrast, WIndex exhibited a trend close to zero for NeC. In this context, monitoring HH and BCS can be useful to avoid losses in the weaning efficiency of cows. Furthermore, we suggest that one way to mitigate efficiency losses in calf production could involve stabilizing the BW of cows and increasing the weaning weight of calves using the WIndex.

Genome-wide scan reveals genomic regions and candidate genes underlying direct and maternal effects of preweaning calf mortality in Nellore cattle

We conducted analysis to estimate genetic parameters and to identify genomic regions and candidate genes affecting direct and maternal effects of preweaning calf mortality (PWM) in Nellore cattle. Phenotypic records of 67,196 animals, and 8443 genotypes for 410,936 SNPs were used. Analysis were performed through the weighted single-step GBLUP approach and considering a threshold animal model via Bayesian Inference. Direct and maternal heritability estimates were of 0.2143 ± 0.0348 and 0.0137 ± 0.0066, respectively. The top 10 genomic regions accounted for 13.61 and 14.23% of the direct and maternal additive genetic variances and harbored a total of 63 and 91 positional candidate genes, respectively. Two overlapping regions on BTA2 were identified for both direct and maternal effects. Candidate genes are involved in biological mechanisms i.e. embryogenesis, immune response, feto-maternal communication, circadian rhythm, hormone alterations, myometrium adaptation, and milk secretion, which are critical for the successful calf growth and survival during preweaning period.

Genome-enabled prediction of reproductive traits in Nellore cattle using parametric models and machine learning methods

This study aimed to assess the predictive ability of different machine learning (ML) methods for genomic prediction of reproductive traits in Nellore cattle. The studied traits were age at first calving (AFC), scrotal circumference (SC), early pregnancy (EP) and stayability (STAY). The numbers of genotyped animals and SNP markers available were 2342 and 321 419 (AFC), 4671 and 309 486 (SC), 2681 and 319 619 (STAY) and 3356 and 319 108 (EP). Predictive ability of support vector regression (SVR), Bayesian regularized artificial neural network (BRANN) and random forest (RF) were compared with results obtained using parametric models (genomic best linear unbiased predictor, GBLUP, and Bayesian least absolute shrinkage and selection operator, BLASSO). A 5-fold cross-validation strategy was performed and the average prediction accuracy (ACC) and mean squared errors (MSE) were computed. The ACC was defined as the linear correlation between predicted and observed breeding values for categorical traits (EP and STAY) and as the correlation between predicted and observed adjusted phenotypes divided by the square root of the estimated heritability for continuous traits (AFC and SC). The average ACC varied from low to moderate depending on the trait and model under consideration, ranging between 0.56 and 0.63 (AFC), 0.27 and 0.36 (SC), 0.57 and 0.67 (EP), and 0.52 and 0.62 (STAY). SVR provided slightly better accuracies than the parametric models for all traits, increasing the prediction accuracy for AFC to around 6.3 and 4.8% compared with GBLUP and BLASSO respectively. Likewise, there was an increase of 8.3% for SC, 4.5% for EP and 4.8% for STAY, comparing SVR with both GBLUP and BLASSO. In contrast, the RF and BRANN did not present competitive predictive ability compared with the parametric models. The results indicate that SVR is a suitable method for genome-enabled prediction of reproductive traits in Nellore cattle. Further, the optimal kernel bandwidth parameter in the SVR model was trait-dependent, thus, a fine-tuning for this hyper-parameter in the training phase is crucial.

Genomic reaction norm models exploiting genotype × environment interaction on sexual precocity indicator traits in Nellore cattle

Brazilian beef cattle are raised predominantly on pasture in a wide range of environments. In this scenario, genotype by environment (G×E) interaction is an important source of phenotypic variation in the reproductive traits. Hence, the evaluation of G×E interactions for heifer's early pregnancy (HP) and scrotal circumference (SC) traits in Nellore cattle, belonging to three breeding programs, was carried out to determine the animal's sensitivity to the environmental conditions (EC). The dataset consisted of 85 874 records for HP and 151 553 records for SC, from which 1800 heifers and 3343 young bulls were genotyped with the BovineHD BeadChip. Genotypic information for 826 sires was also used in the analyses. EC levels were based on the contemporary group solutions for yearling body weight. Linear reaction norm models (RNM), using pedigree information (RNM_A) or pedigree and genomic information (RNM_H), were used to infer G×E interactions. Two validation schemes were used to assess the predictive ability, with the following training populations: (a) forward scheme-dataset was split based on year of birth from 2008 for HP and from 2011 for SC; and (b) environment-specific scheme-low EC (-3.0 and -1.5) and high EC (1.5 and 3.0). The inclusion of the H matrix in RNM increased the genetic variance of the intercept and slope by 18.55 and 23.00% on average respectively, and provided genetic parameter estimates that were more accurate than those considering pedigree only. The same trend was observed for heritability estimates, which were 0.28-0.56 for SC and 0.26-0.49 for HP, using RNM_H, and 0.26-0.52 for SC and 0.22-0.45 for HP, using RNM_A. The lowest correlation observed between unfavorable (-3.0) and favorable (3.0) EC levels were 0.30 for HP and -0.12 for SC, indicating the presence of G×E interaction. The G×E interaction effect implied differences in animals' genetic merit and re-ranking of animals on different environmental conditions. SNP marker-environment interaction was detected for Nellore sexual precocity indicator traits with changes in effect and variance across EC levels. The RNM_H captured G×E interaction effects better than RNM_A and improved the predictive ability by around 14.04% for SC and 21.31% for HP. Using the forward scheme increased the overall predictive ability for SC (20.55%) and HP (11.06%) compared with the environment-specific scheme. The results suggest that the inclusion of genomic information combined with the pedigree to assess the G×E interaction leads to more accurate variance components and genetic parameter estimates.

Reaction norm for yearling weight in beef cattle using single-step genomic evaluation

When the environment on which the animals are raised is very diverse, selecting the best sires for different environments may require the use of models that account for genotype by environment interaction (G × E). The main objective of this study was to evaluate the existence of G × E for yearling weight (YW) in Nellore cattle using reaction norm models with only pedigree and pedigree combined with genomic relationships. Additionally, genomic regions associated with each environment gradient were identified. A total of 67,996 YW records were used in reaction norm models to calculate EBV and genomic EBV. The method of choice for genomic evaluations was single-step genomic BLUP (ssGBLUP). Traditional and genomic models were tested on the ability to predict future animal performance. Genetic parameters for YW were obtained with the average information restricted maximum likelihood method, with and without adding genomic information for 5,091 animals. Additive genetic variances explained by windows of 200 adjacent SNP were used to identify genomic regions associated with the environmental gradient. Estimated variance components for the intercept and the slope in traditional and genomic models were similar. In both models, the observed changes in heritabilities and genetic correlations for YW across environments indicate the occurrence of genotype by environment interactions. Both traditional and genomic models were capable of identifying the genotype by environment interaction; however, the inclusion of genomic information in reaction norm models improved the ability to predict animals' future performance by 7.9% on average. The proportion of genetic variance explained by the top SNP window was 0.77% for the regression intercept (BTA5) and 0.82% for the slope (BTA14). Single-step GBLUP seems to be a suitable model to predict genetic values for YW in different production environments.

Bayesian genome-wide association analysis for body weight in farmed Atlantic salmon (Salmo salar L.)

We performed a genome-wide association study to detect markers associated with growth traits in Atlantic salmon. The analyzed traits included body weight at tagging (BWT) and body weight at 25 months (BW25M). Genotypes of 4662 animals were imputed from the 50K SNP chip to the 200K SNP chip using fimpute software. The markers were simultaneously modeled using Bayes C to identify genomic regions associated with the traits. We identified windows explaining a maximum of 3.71% and 3.61% of the genetic variance for BWT and BW25M respectively. We found potential candidate genes located within the top ten 1-Mb windows for BWT and BW25M. For instance, the vitronectin (VTN) gene, which has been previously reported to be associated with cell growth, was found within one of the top ten 1-Mb windows for BWT. In addition, the WNT1-inducible-signaling pathway protein 3, melanocortin 2 receptor accessory protein 2, myosin light chain kinase, transforming growth factor beta receptor type 3 and myosin light chain 1 genes, which have been reported to be associated with skeletal growth in humans, growth stimulation during the larval stage in zebrafish, body weight in pigs, feed conversion in chickens and growth rate of sheep skeletal muscle respectively, were found within some of the top ten 1-Mb windows for BW25M. These results indicate that growth traits are most likely controlled by many variants with relatively small effects in Atlantic salmon. The genomic regions associated with the traits studied here may provide further insight into the functional regions underlying growth traits in this species.

Genetic analysis of feet and leg conformation traits in Nelore cattle

Feet and leg conformation scores are important traits in beef cattle because they encompass a wide range of locomotion disorders that can lead to productive and reproductive losses. Thus, the study of feet and legs in beef cattle is essential for evaluating possible responses to selection focusing on minimizing economic losses caused by the occurrence of feet and leg problems. The aim of this study was to estimate variance components for feet and leg conformation traits in Nelore cattle. The data set contained records of approximately 300,000 animals that were born between 2000 and 2013. These animals belonged to the commercial beef cattle breeding program of the CRV Lagoa (). Feet and legs were evaluated by assigning visual scores at 2 different time points: feet and leg evaluated as a binary trait (FL1), measured at yearling (about 550 d of age) to identify whether (or not) an animal has feet and leg defects, and feet and leg score (FL2), ranging from 1 (less desirable) to 5 (more desirable) was assigned to the top 20% of animals according to the selection index adopted by the beef cattle breeding program, which was measured 2 to 5 mo after the yearling evaluation. The FL1 and FL2 traits were analyzed together with yearling weight (YW). The (co)variance components and breeding values were estimated by Bayesian inference using 2-trait animal models. The posterior means (standard errors) of the heritabilities for FL1, FL2, and YW were 0.18 (0.04), 0.39 (0.07), and 0.47 (0.01), respectively. The results indicate that the incidence of feet and leg problems in this population might be reduced by selection. The genetic correlation between FL1 and FL2 (-0.47) was moderate and negative as expected because the classification score that holds up each trait has opposite numerical values. The genetic trends estimated for FL1 and FL2 (-0.042 and 0.021 genetic standard deviations per year, respectively) were favorable and they indicate that the independent culling strategy for feet and leg problems promotes favorable changes and contributes to the genetic progress of these traits in the population under study.

Genetic control of residual variance of yearling weight in Nellore beef cattle

There is evidence for genetic variability in residual variance of livestock traits, which offers the potential for selection for increased uniformity of production. Different statistical approaches have been employed to study this topic; however, little is known about the concordance between them. The aim of our study was to investigate the genetic heterogeneity of residual variance on yearling weight (YW; 291.15 ± 46.67) in a Nellore beef cattle population; to compare the results of the statistical approaches, the two-step approach and the double hierarchical generalized linear model (DHGLM); and to evaluate the effectiveness of power transformation to accommodate scale differences. The comparison was based on genetic parameters, accuracy of EBV for residual variance, and cross-validation to assess predictive performance of both approaches. A total of 194,628 yearling weight records from 625 sires were used in the analysis. The results supported the hypothesis of genetic heterogeneity of residual variance on YW in Nellore beef cattle and the opportunity of selection, measured through the genetic coefficient of variation of residual variance (0.10 to 0.12 for the two-step approach and 0.17 for DHGLM, using an untransformed data set). However, low estimates of genetic variance associated with positive genetic correlations between mean and residual variance (about 0.20 for two-step and 0.76 for DHGLM for an untransformed data set) limit the genetic response to selection for uniformity of production while simultaneously increasing YW itself. Moreover, large sire families are needed to obtain accurate estimates of genetic merit for residual variance, as indicated by the low heritability estimates (<0.007). Box-Cox transformation was able to decrease the dependence of the variance on the mean and decreased the estimates of genetic parameters for residual variance. The transformation reduced but did not eliminate all the genetic heterogeneity of residual variance, highlighting its presence beyond the scale effect. The DHGLM showed higher predictive ability of EBV for residual variance and therefore should be preferred over the two-step approach.

Accuracies of genomic prediction of feed efficiency traits using different prediction and validation methods in an experimental Nelore cattle population

Animal feeding is the most important economic component of beef production systems. Selection for feed efficiency has not been effective mainly due to difficult and high costs to obtain the phenotypes. The application of genomic selection using SNP can decrease the cost of animal evaluation as well as the generation interval. The objective of this study was to compare methods for genomic evaluation of feed efficiency traits using different cross-validation layouts in an experimental beef cattle population genotyped for a high-density SNP panel (BovineHD BeadChip assay 700k, Illumina Inc., San Diego, CA). After quality control, a total of 437,197 SNP genotypes were available for 761 Nelore animals from the Institute of Animal Science, Sertãozinho, São Paulo, Brazil. The studied traits were residual feed intake, feed conversion ratio, ADG, and DMI. Methods of analysis were traditional BLUP, single-step genomic BLUP (ssGBLUP), genomic BLUP (GBLUP), and a Bayesian regression method (BayesCπ). Direct genomic values (DGV) from the last 2 methods were compared directly or in an index that combines DGV with parent average. Three cross-validation approaches were used to validate the models: 1) YOUNG, in which the partition into training and testing sets was based on year of birth and testing animals were born after 2010; 2) UNREL, in which the data set was split into 3 less related subsets and the validation was done in each subset a time; and 3) RANDOM, in which the data set was randomly divided into 4 subsets (considering the contemporary groups) and the validation was done in each subset at a time. On average, the RANDOM design provided the most accurate predictions. Average accuracies ranged from 0.10 to 0.58 using BLUP, from 0.09 to 0.48 using GBLUP, from 0.06 to 0.49 using BayesCπ, and from 0.22 to 0.49 using ssGBLUP. The most accurate and consistent predictions were obtained using ssGBLUP for all analyzed traits. The ssGBLUP seems to be more suitable to obtain genomic predictions for feed efficiency traits on an experimental population of genotyped animals.

Genotype × environment interaction for age at first calving, scrotal circumference, and yearling weight in Nellore cattle using reaction norms in multitrait random regression models

The aim of this study was to evaluate the effect of genotype × environment interaction (G×E) on age at first calving (AFC), scrotal circumference (SC), and yearling weight (YW) and to estimate genetic correlations between these traits in Nellore cattle using reaction norms in multitrait random regression models. In this study, 28,871, 41,386, and 89,152 records of Nellore cattle for AFC, SC, and YW, respectively, were used. The data were obtained from farms located in the north, northeast, midwest, and southeast regions of Brazil that participate in the DeltaGen Breeding Program. Environmental levels were defined as a function of contemporary groups, that is, animals born in the same herd and year, from the same management group (from birth to yearling), and of the same sex. Postweaning weight gain was used as a criterion to evaluate the environmental conditions for all traits. For reaction norm analyses, residual variances were modeled with homogeneous and heterogeneous classes. The model for SC and YW included the fixed effects of contemporary group and age of the animal as a covariate as well as random direct additive genetic and residual effects. The same model, excluding the covariate age of the animal, was used for AFC. The heritability estimates were low to high for AFC (0.09 to 0.50), high for SC (0.51 to 0.67), and moderate to high for YW (0.33 to 0.71). The genetic correlations (within each trait) along the environmental levels varied from -0.27 to 1.0 for AFC, from 0.73 to 1.0 for SC, and from 0.26 to 1.0 for YW. The genetic correlations between different traits in different environments varied from -0.14 to -0.60 between AFC and SC, from -0.05 to -0.32 between AFC and YW, and from -0.05 to 0.72 between YW and SC. The genetic correlations have had different magnitudes for AFC, SC, and YW, which could indicate the presence of G×E. The present results should support researchers and farmers in defining selection criteria to improve growth traits and sexual precocity. Our results suggest that animals for breeding have to be selected in the same environment and management conditions as their progeny will be reared.

Genetic parameters for an alternative criterion to improve productive longevity of Nellore cows

Number of calvings at 53 mo (NC53) was proposed as an alternative selection criterion to improve productive longevity of Nellore cows. This study was carried out to estimate variance components for NC53 by assuming different models, so that the potential for using this selection criterion to improve fertility of Nellore cows could be assessed. Genetic correlations between NC53, number of calvings at 89 mo (NC89), and 2 selection indexes used in this breed were also estimated. The NC53 trait is moderately heritable (posterior mean heritability ≈ 0.17) and selecting for this criterion could improve productive longevity of Nellore cows. Greater response to selection is expected by fitting a threshold animal model for this trait, rather than a linear animal model. Greater accuracy of prediction for this criterion could be achieved by fitting a threshold-linear model, considering this trait and a selection index composed by traits evaluated at weaning and long-yearling.

Analysis of genetic correlations of hip height with selection indices and mature weight in Nelore cattle

Body size is directly related to the productive and reproductive performance of beef cattle raised under free-range conditions. In an attempt to better plan selection criteria, avoiding extremes in body size, this study estimated the heritabilities and genetic correlations of yearling hip height (YH) and mature hip height (MH) with selection indices obtained at weaning (WI) and yearling (YI) and mature weight (MW). Data from 102,373 Nelore animals born between 1984 and 2010, which belong to 263 farms that participate in genetic evaluation programmes of beef cattle conducted in Brazil and Paraguay, were used. The (co)variance components and genetic parameters were estimated by Bayesian inference in multi-trait analysis using an animal model. The mean heritabilities for YH, MH and MW were 0.56 ± 0.06, 0.47 ± 0.02 and 0.42 ± 0.02, respectively. The genetic correlation of YH with WI (0.13 ± 0.01) and YI (0.11 ± 0.01) was practically zero, whereas a higher correlation was observed with MW (0.22 ± 0.03). Positive genetic correlations of medium magnitude were estimated between MH and WI and YI (0.23 ± 0.01 and 0.43 ± 0.02, respectively). On the other hand, a high genetic correlation (0.68 ± 0.03) was observed between the indicator traits of mature body size (MH and MW). Considering the top 20 % of sire (896 sires) in terms of breeding values for the yearling index, the rank sire correlations between breeding values for MH and MW was 0.62. In general, the results indicate that selection based on WI and YI should not lead to important changes in YH. However, an undesired correlated response in mature cow height is expected, particularly when selection is performed using YI. Therefore, changes in the body structure of Nelore females can be obtained when MH and MW is used as a selection criterion for cows.

Evaluation of mature cow weight: genetic correlations with traits used in selection indices, correlated responses, and genetic trends in Nelore cattle

Genetic correlations of selection indices and the traits considered in these indices with mature weight (MW) of Nelore females and correlated responses were estimated to determine whether current selection practices will result in an undesired correlated response in MW. Genetic trends for weaning and yearling indices and MW were also estimated. Data from 612,244 Nelore animals born between 1984 and 2010, belonging to different beef cattle evaluation programs from Brazil and Paraguay, were used. The following traits were studied: weaning conformation (WC), weaning precocity (WP), weaning muscling (WM), yearling conformation (YC), yearling precocity (YP), yearling muscling (YM), weaning and yearling indices, BW gain from birth to weaning (BWG), postweaning BW gain (PWG), scrotal circumference (SC), and MW. The variance and covariance components were estimated by Bayesian inference in a multitrait analysis, including all traits in the same analysis, using a nonlinear (threshold) animal model for visual scores and a linear animal model for the other traits. The mean direct heritabilities were 0.21±0.007 (WC), 0.22±0.007 (WP), 0.20±0.007 (WM), 0.43±0.005 (YC), 0.40±0.005 (YP), 0.40±0.005 (YM), 0.17±0.003 (BWG), 0.21±0.004 (PWG), 0.32±0.001 (SC), and 0.44±0.018 (MW). The genetic correlations between MW and weaning and yearling indices were positive and of medium magnitude (0.30±0.01 and 0.31±0.01, respectively). The genetic changes in weaning index, yearling index, and MW, expressed as units of genetic SD per year, were 0.26, 0.27, and 0.01, respectively. The genetic trend for MW was nonsignificant, suggesting no negative correlated response. The selection practice based on the use of sires with high final index giving preference for those better ranked for yearling precocity and muscling than for conformation generates only a minimal correlated response in MW.

Índice de seleção bioeconômico para fêmeas de corte da raça nelore

A eficiência reprodutiva da fêmea Nelore foi descrita com base na precocidade sexual, permanência produtiva no rebanho (NP), produtividade materna (PM) e custo de mantença estimado (CM). A combinação dessas características deu origem ao índice bioeconômico retorno maternal (RMat), que estima o retorno em quilos de peso vivo produzidos por uma vaca em um ano. Em adição, incluiu-se a composição do peso produzido, adicionando à PM os escores de conformação, precocidade e musculatura a desmama, compondo o biótipo do bezerro. Foram consideradas precoces as fêmeas cuja idade ao primeiro parto foi inferior a 30 meses. A NP foi expressa pelo número de partos até 53 meses de idade. O CM foi calculado em função do consumo estimado de matéria seca da vaca. O RMat médio estimado foi 62,02±24,12 kg/vaca/ano. As estimativas da variância genética aditiva e residual do RMat, usando a metodologia da máxima verossimilhança restrita, sob um modelo animal unicaracterística, foram 195,35 e 242,96, respectivamente. A herdabilidade estimada para Rmat foi 0,45±0,02, indicando que o índice é herdável e pode ser aplicado na seleção para eficiência reprodutiva. A NP foi o componente de principal variação do RMat. Touros selecionados com base no RMat apresentaram filhas mais eficientes.

Genetic effects on preweaning weight gain of Nelore-Hereford calves according to different models and estimation methods1

Additive and nonadditive genetic effects on preweaning weight gain (PWG) of a commercial crossbred population were estimated using different genetic models and estimation methods. The data set consisted of 103,445 records on purebred and crossbred Nelore-Hereford calves raised under pasture conditions on farms located in south, southeast, and middle west Brazilian regions. In addition to breed additive and dominance effects, the models including different epistasis covariables were tested. Models considering joint additive and environment (latitude) by genetic effects interactions were also applied. In a first step, analyses were carried out under animal models. In a second step, preadjusted records were analyzed using ordinary least squares (OLS) and ridge regression (RR). The results reinforced evidence that breed additive and dominance effects are not sufficient to explain the observed variability in preweaning traits of Bos taurus x Bos indicus calves, and that genotype x environment interaction plays an important role in the evaluation of crossbred calves. Data were ill-conditioned to estimate the effects of genotype x environment interactions. Models including these effects presented multicolinearity problems. In this case, RR seemed to be a powerful tool for obtaining more plausible and stable estimates. Estimated prediction error variances and variance inflation factors were drastically reduced, and many effects that were not significant under ordinary least squares became significant under RR. Predictions of PWG based on RR estimates were more acceptable from a biological perspective. In temperate and subtropical regions, calves with intermediate genetic compositions (close to 1/2 Nelore) exhibited greater predicted PWG. In the tropics, predicted PWG increased linearly as genotype got closer to Nelore.