Rodrigo Junqueira Pereira

Milk yield persistency in Brazilian Gyr cattle based on a random regression model

With the objective of evaluating measures of milk yield persistency, 27,000 test-day milk yield records from 3362 first lactations of Brazilian Gyr cows that calved between 1990 and 2007 were analyzed with a random regression model. Random, additive genetic and permanent environmental effects were modeled using Legendre polynomials of order 4 and 5, respectively. Residual variance was modeled using five classes. The average lactation curve was modeled using a fourth-order Legendre polynomial. Heritability estimates for measures of persistency ranged from 0.10 to 0.25. Genetic correlations between measures of persistency and 305-day milk yield (Y305) ranged from -0.52 to 0.03. At high selection intensities for persistency measures and Y305, few animals were selected in common. As the selection intensity for the two traits decreased, a higher percentage of animals were selected in common. The average predicted breeding values for Y305 according to year of birth of the cows had a substantial annual genetic gain. In contrast, no improvement in the average persistency breeding value was observed. We conclude that selection for total milk yield during lactation does not identify bulls or cows that are genetically superior in terms of milk yield persistency. A measure of persistency represented by the sum of deviations of estimated breeding value for days 31 to 280 in relation to estimated breeding value for day 30 should be preferred in genetic evaluations of this trait in the Gyr breed, since this measure showed a medium heritability and a genetic correlation with 305-day milk yield close to zero. In addition, this measure is more adequate at the time of peak lactation, which occurs between days 25 and 30 after calving in this breed.

Breeding value accuracy estimates for growth traits using random regression and multi-trait models in Nelore cattle.

We quantified the potential increase in accuracy of expected breeding value for weights of Nelore cattle, from birth to mature age, using multi-trait and random regression models on Legendre polynomials and B-spline functions. A total of 87,712 weight records from 8144 females were used, recorded every three months from birth to mature age from the Nelore Brazil Program. For random regression analyses, all female weight records from birth to eight years of age (data set I) were considered. From this general data set, a subset was created (data set II), which included only nine weight records: at birth, weaning, 365 and 550 days of age, and 2, 3, 4, 5, and 6 years of age. Data set II was analyzed using random regression and multi-trait models. The model of analysis included the contemporary group as fixed effects and age of dam as a linear and quadratic covariable. In the random regression analyses, average growth trends were modeled using a cubic regression on orthogonal polynomials of age. Residual variances were modeled by a step function with five classes. Legendre polynomials of fourth and sixth order were utilized to model the direct genetic and animal permanent environmental effects, respectively, while third-order Legendre polynomials were considered for maternal genetic and maternal permanent environmental effects. Quadratic polynomials were applied to model all random effects in random regression models on B-spline functions. Direct genetic and animal permanent environmental effects were modeled using three segments or five coefficients, and genetic maternal and maternal permanent environmental effects were modeled with one segment or three coefficients in the random regression models on B-spline functions. For both data sets (I and II), animals ranked differently according to expected breeding value obtained by random regression or multi-trait models. With random regression models, the highest gains in accuracy were obtained at ages with a low number of weight records. The results indicate that random regression models provide more accurate expected breeding values than the traditionally finite multi-trait models. Thus, higher genetic responses are expected for beef cattle growth traits by replacing a multi-trait model with random regression models for genetic evaluation. B-spline functions could be applied as an alternative to Legendre polynomials to model covariance functions for weights from birth to mature age.

Inbreeding depression in Zebu cattle traits

The productivity of herds may be negatively affected by inbreeding depression, and it is important to know how intense is this effect on the livestock performance. We performed a comprehensive analysis involving five Zebu breeds reared in Brazil to estimate inbreeding depression in productive and reproductive traits. Inbreeding depression was estimated for 13 traits by including the individual inbreeding rate as a linear covariate in the standard genetic evaluation models. For all breeds and for almost all traits (no effect was observed on gestation length), the performance of the animals was compromised by an increase in inbreeding. The average inbreeding depression was -0.222% and -0.859% per 1% of inbreeding for linear regression coefficients scaled on the percentage of mean (βm ) and standard deviation (βσ ), respectively. The means for βm (and βσ ) were -0.269% (-1.202%) for weight/growth traits and -0.174% (-0.546%) for reproductive traits. Hence, inbreeding depression is more pronounced in weight/growth traits than in reproductive traits. These findings highlight the need for the management of inbreeding in the respective breeding programmes of the breeds studied here.

Detrimental effect of selection for milk yield on genetic tolerance to heat stress in purebred Zebu cattle: Genetic parameters and trends

In an attempt to determine the possible detrimental effects of continuous selection for milk yield on the genetic tolerance of Zebu cattle to heat stress, genetic parameters and trends of the response to heat stress for 86,950 test-day (TD) milk yield records from 14,670 first lactations of purebred dairy Gir cows were estimated. A random regression model with regression on days in milk (DIM) and temperature-humidity index (THI) values was applied to the data. The most detrimental effect of THI on milk yield was observed in the stage of lactation with higher milk production, DIM 61 to 120 (-0.099kg/d per THI). Although modest variations were observed for the THI scale, a reduction in additive genetic variance as well as in permanent environmental and residual variance was observed with increasing THI values. The heritability estimates showed a slight increase with increasing THI values for any DIM. The correlations between additive genetic effects across the THI scale showed that, for most of the THI values, genotype by environment interactions due to heat stress were less important for the ranking of bulls. However, for extreme THI values, this type of genotype by environment interaction may lead to an important error in selection. As a result of the selection for milk yield practiced in the dairy Gir population for 3 decades, the genetic trend of cumulative milk yield was significantly positive for production in both high (51.81kg/yr) and low THI values (78.48kg/yr). However, the difference between the breeding values of animals at high and low THI may be considered alarming (355kg in 2011). The genetic trends observed for the regression coefficients related to general production level (intercept of the reaction norm) and specific ability to respond to heat stress (slope of the reaction norm) indicate that the dairy Gir population is heading toward a higher production level at the expense of lower tolerance to heat stress. These trends reflect the genetic antagonism between production and tolerance to heat stress demonstrated by the negative genetic correlation between these components (-0.23). Monitoring trends of the genetic component of heat stress would be a reasonable measure to avoid deterioration in one of the main traits of Zebu cattle (i.e., high tolerance to heat stress). On the basis of current genetic trends, the need for future genetic evaluation of dairy Zebu animals for tolerance to heat stress cannot be ruled out.

Comparison of factor-analytic and reduced rank models for test-day milk yield in Gyr dairy cattle (Bos indicus).

We analyzed 46,161 monthly test-day records of milk production from 7453 first lactations of crossbred dairy Gyr (Bos indicus) x Holstein cows. The following seven models were compared: standard multivariate model (M10), three reduced rank models fitting the first 2, 3, or 4 genetic principal components, and three models considering a 2-, 3-, or 4-factor structure for the genetic covariance matrix. Full rank residual covariance matrices were considered for all models. The model fitting the first two principal components (PC2) was the best according to the model selection criteria. Similar phenotypic, genetic, and residual variances were obtained with models M10 and PC2. The heritability estimates ranged from 0.14 to 0.21 and from 0.13 to 0.21 for models M10 and PC2, respectively. The genetic correlations obtained with model PC2 were slightly higher than those estimated with model M10. PC2 markedly reduced the number of parameters estimated and the time spent to reach convergence. We concluded that two principal components are sufficient to model the structure of genetic covariances between test-day milk yields.

Effect of pregnancy on the genetic evaluation of dairy cattle.

We investigated the effect of stage of pregnancy on estimates of breeding values for milk yield and milk persistency in Gyr and Holstein dairy cattle in Brazil. Test-day milk yield records were analyzed using random regression models with or without the effect of pregnancy. Models were compared using residual variances, heritabilities, rank correlations of estimated breeding values of bulls and cows, and number of nonpregnant cows in the top 200 for milk yield and milk persistency. The estimates of residual variance and heritabilities obtained with the models with or without the effect of pregnancy were similar for the two breeds. Inclusion of the effect of pregnancy in genetic evaluation models for these populations did not affect the ranking of cows and sires based on their predicted breeding values for 305-day cumulative milk yield. In contrast, when we examined persistency of milk yield, lack of adjustment for the effect of pregnancy overestimated breeding values of nonpregnant cows and cows with a long days open period and underestimated breeding values of cows with a short days open period. We recommend that models include the effect of days of pregnancy for estimation of adjustment factors for the effect of pregnancy in genetic evaluations of Dairy Gyr and Holstein cattle.

Genetics of heat tolerance for milk yield and quality in Holsteins

Tropical and sub-tropical climates are characterized by high temperature and humidity, during at least part of the year. Consequently, heat stress is common in Holstein cattle and productive and reproductive losses are frequent. Our objectives were as follows: (1) to quantify losses in production and quality of milk due to heat stress; (2) to estimate genetic correlations within and between milk yield (MY) and milk quality traits; and (3) to evaluate the trends of genetic components of tolerance to heat stress in multiple lactations of Brazilian Holstein cows. Thus, nine analyses using two-trait random regression animal models were carried out to estimate variance components and genetic parameters over temperature-humidity index (THI) values for MY and milk quality traits (three lactations: MY×fat percentage (F%), MY×protein percentage (P%) and MY×somatic cell score (SCS)) of Brazilian Holstein cattle. It was demonstrated that the effects of heat stress can be harmful for traits related to milk production and milk quality of Holstein cattle even though most herds were maintained in a modified environment, for example, with fans and sprinklers. For MY, the effect of heat stress was more detrimental in advanced lactations (-0.22 to -0.52 kg/day per increase of 1 THI unit). In general, the mean heritability estimates were higher for lower THI values and longer days in milk for all traits. In contrast, the heritability estimates for SCS increased with increasing THI values in the second and third lactation. For each trait studied, lower genetic correlations (different from unity) were observed between opposite extremes of THI (THI 47 v. THI 80) and in advanced lactations. The genetic correlations between MY and milk quality trait varied across the THI scale and lactations. The genotype×environment interaction due to heat stress was more important for MY and SCS, particularly in advanced lactations, and can affect the genetic relationship between MY and milk quality traits. Selection for higher MY, F% or P% may result in a poor response of the animals to heat stress, as a genetic antagonism was observed between the general production level and specific ability to respond to heat stress for these traits. Genetic trends confirm the adverse responses in the genetic components of heat stress over the years for milk production and quality. Consequently, the selection of Holstein cattle raised in modified environments in both tropical and sub-tropical regions should take into consideration the genetic variation in heat stress.

Contemporary group alternatives for genetic evaluation of milk yield in small populations of dairy cattle

In the present study, different random regression models, focussed on the nature of the effect of the contemporary group, fixed or random, were compared for genetic evaluation of test-day milk yield in dairy Gir breed, whose herds are characterised by relatively reduced numbers of cows. Cows were assigned to the same contemporary group if they were tested in the same herd, year and month. In one of the evaluated models, the contemporary group was investigated as a fixed effect and a clustering procedure was adopted to group herd-year subclasses in order to, at the end of the process, all contemporary groups contained at least three cows. The models were compared by the quality of fit, the accuracy of the predicted breeding values and their estimates of genetic parameters. The clustering procedure did not significantly improve the accuracy of predicted breeding values. Moreover, the fit to the data for this model was negatively affected. Therefore, this strategy should not be further implemented. The model including the herd-calving period effect and contemporary group treated as random effect showed similar characteristics to its equivalent, where the contemporary group was modelled as a fixed effect. However, the fit to the data for this model was slightly worse. Thus, the results suggest a random regression model including the herd-calving period effect and the fixed effect of contemporary group for the genetic evaluations of production traits in dairy Gir cattle. These findings could be extended to small dairy cattle populations whose herds are small-sized.

Fatores ambientais e parâmetros genéticos para características reprodutivas de fêmeas Nelore na região meio norte do Brasil

The study was carried out to evaluate the age-at-first-calving and calving interval (IPP and IDP) in Nelore cows, using IPP and IDP data from 23,066 females reared in a grazing regime in the States of Maranhao and Piaui from 1976 to 2012. Data were analyzed statistically using the SAS (Statistical Analysis System) program. The mean and coefficient of variation for the age-at-first-calving was 1252.77 ± 125 days (10.41 months; 10.56%) and, for the calving interval, values were 591.73 ± 149.91 days (12.49 months; 38.74%). Uncharacteristic analyses were used to estimate genetic parameters for heritability values, obtaining values of approximately 0.026 and 0.006 for IPP and IDP, respectively. The repeatability estimate for the IDP characteristic, as calculated by uncharacteristic analysis, had a generally low value of 0.025. Similarly, for trait analysis, heritability values were 0.03 and 0.06 for IPP and IDP, respectively, and the estimated repeatability coefficient was 0.14. The genetic correlations between the reproductive characteristics (IPP and IDP) observed in this study were equal to 0.06. The reproductive performance of the herd studied, in tropical climate conditions, can be considered to be satisfactory, in accordance with the parameters obtained by other authors in different regions.