No abstract

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records.

The monitoring of population structure, inbreeding and other related parameters has great potential to prevent major losses of genetic diversity in populations of Zebu cattle in Brazil. Therefore, the objective of the present study was to investigate the structure and genetic diversity of Brazilian Zebu cattle breeds by pedigree analysis. The national pedigree file of the seven Brazilian Zebu breeds was used, which included all registered animals (12,290,243) born between 1938 and 2012: Brahman, Gir, Guzerá, Indubrasil, Nelore, Sindi, and Tabapuã. Almost all breeds studied undergo expansion in their census which, however, is not accompanied by the maintenance of genetic diversity. Problems were encountered in all breeds, but most of them can currently be considered less important. Using the calculation method considered as the most accurate for pedigree analysis when some population substructure exists, all breeds, except Sindi, had effective population size greater than 100. The most common problems were the presence of tight bottlenecks in the pedigree, which were mainly due to the intensive use of few animals as parents and the high degree of population subdivision. The use of a wider range of sires is therefore recommended. However, most Zebu breeds can deal with breeding programs using high selection intensities. Greater care should be taken in the case of the Indubrasil breed since its census was reduced drastically over the last few years, a fact favoring the occurrence of serious problems related to inbreeding. Although Sindi presents problems due to subdivision and possesses a relatively small census compared to other Zebu breeds, this population would have a promising future if its breeding policy were revised.

The present study had the following objectives: to compare random regression models (RRM) considering the time-dependent (days in milk, DIM) and/or temperature × humidity-dependent (THI) covariate for genetic evaluation; to identify the effect of genotype by environment interaction (G×E) due to heat stress on milk yield; and to quantify the loss of milk yield due to heat stress across lactation of cows under tropical conditions. A total of 937,771 test-day records from 3603 first lactations of Brazilian Holstein cows obtained between 2007 and 2013 were analyzed. An important reduction in milk yield due to heat stress was observed for THI values above 66 (-0.23 kg/day/THI). Three phases of milk yield loss were identified during lactation, the most damaging one at the end of lactation (-0.27 kg/day/THI). Using the most complex RRM, the additive genetic variance could be altered simultaneously as a function of both DIM and THI values. This model could be recommended for the genetic evaluation taking into account the effect of G×E. The response to selection in the comfort zone (THI ≤ 66) is expected to be higher than that obtained in the heat stress zone (THI > 66) of the animals. The genetic correlations between milk yield in the comfort and heat stress zones were less than unity at opposite extremes of the environmental gradient. Thus, the best animals for milk yield in the comfort zone are not necessarily the best in the zone of heat stress and, therefore, G×E due to heat stress should not be neglected in the genetic evaluation.

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.