Genetic parameters of total milk yield and factors describing the shape of lactation curve in dairy buffaloes

The objective of this study was to apply factor analysis to describe lactation curves in dairy buffaloes in order to estimate the phenotypic and genetic association between common latent factors and cumulative milk yield. A total of 31 257 monthly test-day milk yield records from buffaloes belonging to herds located in the state of São Paulo were used to estimate two common latent factors, which were then analysed in a multi-trait animal model for estimating genetic parameters. Estimates of (co)variance components for the two common latent factors and cumulated 270-d milk yield were obtained by Bayesian inference using a multiple trait animal model. Contemporary group, number of milkings per day (two levels) and age of buffalo cow at calving (linear and quadratic) as covariate were included in the model as fixed effects. The additive genetic, permanent environmental and residual effects were included as random effects. The first common latent factor (F1) was associated with persistency of lactation and the second common latent factor (F2) with the level of production in early lactation. Heritability estimates for Fl and F2 were 0.12 and 0.07, respectively. Genetic correlation estimates between El and F2 with cumulative milk yield were positive and moderate (0.63 and 0.52). Multivariate statistics employing factor analysis allowed the extraction of two variables (latent factors) that described the shape of the lactation curve. It is expected that the response to selection to increase lactation persistency is higher than the response obtained from selecting animals to increase lactation peak. Selection for higher total milk yield would result in a favourable correlated response to increase the level of production in early lactation and the lactation persistency.

Genetic parameters for milk yield, age at first calving and interval between first and second calving in milk buffaloes

Genetic parameters for the relation between the traits of milk yield (MY), age at first calving (AFC) and interval between first and second calving (IBFSC) were estimated in milk buffaloes of the Murrah breed. In the study, data of 1578 buffaloes at first lactation, with calvings from 1974 to 2006 were analyzed. The MTDFREML system was used in the analyses with models for the MY, IBFSC traits which included the fixed effects of herd-year-season of calving, linear and quadratic terms of calving age as covariate and the random animal effects and error. The model for AFC consisted of the herd-year-season fixed effects of calving and the random effects of animal and error. Heritability estimates MY, AFC and IBFSC traits were 0.20, 0.07 and 0.14, respectively. Genetic and phenotypic correlations between the traits were: MY and AFC = -0.12 and -0.15, MY and IBFSC = 0.07 and 0.30, AFC and IBFSC = 0.35 and 0.37, respectively. Genetic correlation between MY and AFC traits showed desirable negative association, suggesting that the daughters of the bulls with high breeding value for MY could be physiological maturity to a precocious age. Genetic correlation between MY and IBFSC showed that the selection of the animals that increased milk yield is also those that tend to intervals of bigger calving.

Variances of direct genetic effects, maternal genetic effects, and cytoplasmic inheritance effects for milk yield, fat yield, and fat percentage

Milk yield, fat yield, and fat percentage during the first three lactations were studied using New York Holsteins that were milked twice daily over a 305-d, mature equivalent lactation. Those data were used to estimate variances from direct and maternal genetic effects, cytoplasmic effects, sire by herd interaction, and cow permanent environmental effects. Cytoplasmic line was traced to the last female ancestor using DHI records from 1950 through 1991. Records were 138,869 lactations of 68,063 cows calving from 1980 through 1991. Ten random samples were based on herd code. Samples averaged 4926 dams and 2026 cytoplasmic lines. Model also included herd-year-seasons as fixed effects and genetic covariance for direct-maternal effects. Mean estimates of the effects of maternal genetic variances and direct-maternal covariances, as fractions of phenotypic variances, were 0.008 and 0.007 for milk yield, 0.010 and 0.010 for fat yield, and 0.006 and 0.025 for fat percentage, respectively. Average fractions of variance from cytoplasmic line were 0.011, 0.008, and 0.009 for milk yield, fat yield, and fat percentage. Removal of maternal genetic effects and covariance for maternal direct effects from the model increased the fraction of direct genetic variance by 0.014, 0.021, and 0.046 for milk yield, fat yield, and fat percentage; little change in the fraction was due to cytoplasmic line. Exclusion of cytoplasmic effects from the model increased the ratio of additive direct genetic variance to phenotypic variance by less than 2%. Similarly, when sire by herd interaction was excluded, the ratio of direct genetic variance to phenotypic variance increased 1% or less.

Factors affecting somatic cell counts and their relations with milk and milk constituent yield in buffaloes

Data concerning daily milk yield (MY), percentage of milk fat (%F), protein (%P), lactose (%LT), and total solids (%TS), and somatic cell counts (SCC) for a herd of 222 Murrah buffalo reared in the state of São Paulo, Brazil, were collected monthly from 1997 to 2000 in order to study the factors affecting SCC and their relation to milk production and constituents during lactation. SCC decreased in the second month of lactation and increased thereafter, up to the ninth month of lactation. The interaction of month of lactation x order of calving was significant. Mean MY observed during the first month of lactation was 6.87 kg, which increased to 7.65 kg during the second month, and then decreased until the ninth month of lactation (3.83 kg). During the different months of lactation, %F, %P, %LT, and %TS ranged from 6.28 to 8.38%, 4.05 to 4.59%, 4.96 to 5.34%, and 16.94 to 18.55%, respectively. Calving year, calving order, and order of month of lactation significantly affected MY, %F, %P, %LT, and %TS. The regression coefficients of transformed SCC on MY and %LT were negative and significant during all months of lactation, showing that milk and lactose yield decreased with increased transformed SCC, causing losses to buffalo milk producers.

Factors that cause genotype by environment interaction and use of a multiple-trait herd-cluster model for milk yield of Holstein cattle from Brazil and Colombia

Descriptive herd variables (DVHE) were used to explain genotype by environment interactions (G x E) for milk yield (MY) in Brazilian and Colombian production environments and to develop a herd-cluster model to estimate covariance components and genetic parameters for each herd environment group. Data consisted of 180,522 lactation records of 94,558 Holstein cows from 937 Brazilian and 400 Colombian herds. Herds in both countries were jointly grouped in thirds according to 8 DVHE: production level, phenotypic variability, age at first calving, calving interval, percentage of imported semen, lactation length, and herd size. For each DVHE, REML bivariate animal model analyses were used to estimate genetic correlations for MY between upper and lower thirds of the data. Based on estimates of genetic correlations, weights were assigned to each DVHE to group herds in a cluster analysis using the FASTCLUS procedure in SAS. Three clusters were defined, and genetic and residual variance components were heterogeneous among herd clusters. Estimates of heritability in clusters 1 and 3 were 0.28 and 0.29, respectively, but the estimate was larger (0.39) in Cluster 2. The genetic correlations of MY from different clusters ranged from 0.89 to 0.97. The herd-cluster model based on DVHE properly takes into account G x E by grouping similar environments accordingly and seems to be an alternative to simply considering country borders to distinguish between environments.

Genetic parameters for milk yield, age at first calving and interval between first and second calving in milk Murrah buffaloes

In the present study, data of 1,578 first lactation females, calving from 1985 to 2006 were analysed with the purpose of estimating genetic parameters for milk yield (MY), age at first calving (AFC) and interval between first and second calving (IBFSC) in dairy buffaloes of the Murrah breed in Brazil, Heritability estimates for MY, AFC and IBFSC traits were 0.20, 0.07 and 0.14, respectively. Genetic correlations between MY and AFC and IBFSC were -0.12 and 0.07, respectively, while the corresponding phenotypic correlations were -0.15 and 0.30, respectively. Genetic and phenotypic correlations between AFC and IBFSC were 0.35 and 0.37, respectively. Genetic correlation between MY and AFC showed desirable negative association, suggesting that daughters of the bulls with high breeding values for MY could reach physiological mature at a precocious age. Genetic correlation between MY and IBFSC, showed that the selection for milk production could result in the increase of calving intervals.

Test-day model for milk yield of dairy buffaloes in colombia

The test-day model is the preferred method for genetic evaluations in dairy cattle. For this study, 28372 test-day records of 1220 lactations from 1997 to 2009 were used. The (co)variance components for milk in test-day were estimated using a Uni and multiple-traits repeated animal model with the Restricted Maximum Likelihood method (REML). The Contemporary Group (herd, year, and season of parity) and the age of parity (linear and quadratic) fixed effects, and the additive genetic, permanent environmental, and residual random effects were included in the model. The heritabilities ranged between 0.06 and 0.45 during lactation. The genetic correlations were greater than 0.93. In conclusion, the test-day model is appropriate for the genetic evaluation of dairy buffaloes in Colombia.

Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle

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. © 2013 American Dairy Science Association.

Polymorphisms in the ghrelin gene and their associations with milk yield and quality in water buffaloes

Ghrelin is a gastrointestinal hormone that acts in releasing growth hormone and influences the body general metabolism. It has been proposed as a candidate gene for traits such as growth, carcass quality, and milk production of livestock because it influences feed intake. In this context, the aim of this study was to verify the existence of polymorphisms in the ghrelin gene and their associations with milk, fat and protein yield, and percentage in water buffaloes (Bubalus bubalis). A group of 240 animals was studied. Five primer pairs were used and 11 single nucleotide polymorphisms (SNP) were found in the ghrelin gene by sequencing. The animals were genotyped for 8 SNP by PCR-RFLP. The SNP g.960G>A and g.778C>T were associated with fat yield and the SNP g.905T>C was associated with fat yield and percentage and protein percentage. These SNP are located in intronic regions of DNA and may be in noncoding RNA sites or affect transcriptional efciency. The ghrelin gene in buffaloes influences milk fat and protein synthesis. The polymorphisms observed can be used as molecular markers to assist selection. © 2013 American Dairy Science Association.

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal. © 2013 American Dairy Science Association.

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.

Application of Reactions Norms in Study of Genotype Environmental Interaction for Milk Yield of Buffaloes

The aim of this study was to estimate genetic parameters for milk yield (MY) in buffaloes using reaction norms. Model included the additive direct effect as random and contemporary group (herd and year of birth) were included as fixed effects and cow age classes (linear) as covariables. The animal additive direct random effect was modeled through linear Legendre polynomials on environment gradient (EG) standardized means. Mean trends were taken into account by a linear regression on Legendre polynomials of environmental group means. Residual variance was modeled trough 6 heterogeneity classes (EG). These classes of residual variance was formed : EG1: mean = 866,93 kg (621,68 kg-1011,76 kg); EG2: mean = 1193,00 kg (1011,76 kg-1251,49 kg); EG3: mean = 1309,37 kg (1251,49 kg -1393,20 kg); EG4: mean = 1497,59 kg (1393,20 kg-1593,53 kg); EG5: mean = 1664,78 kg (1593,53 kg -1727,32kg) e EG6: mean = 1973,85 kg (1727,32 kg -2422,19 kg).(Co) variance functions were estimated by restricted maximum likelihood (REML) using the GIBBS3F90 package. The heritability estimates for MY raised as the environmental gradient increased, varying from 0.20 to 0.40. However, in intermediate to favorable environments, the heritability estimates obtained with Considerable genotype-environment interaction was found for MY using reaction norms. For genetic evaluation of MY is necessary to consider heterogeneity of variances to model the residual variance.

Estimates of genetic parameters for total milk yield over multiple ages in Brazilian Murrah buffaloes using different models

The objective of this study was to estimate variance components and genetic parameters for accumulated 305-day milk yield (MY305) over multiple ages, from 24 to 120 months of age, applying random regression (RRM), repeatability (REP) and multi-trait (MT) models. A total of 4472 lactation records from 1882 buffaloes of the Murrah breed were utilized. The contemporary group (herd-year-calving season) and number of milkings (two levels) were considered as fixed effects in all models. For REP and RRM, additive genetic, permanent environmental and residual effects were included as random effects. MT considered the same random effects as did REP and RRM with the exception of permanent environmental effect. Residual variances were modeled by a step function with 1, 4, and 6 classes. The heritabilities estimated with RRM increased with age, ranging from 0.19 to 0.34, and were slightly higher than that obtained with the REP model. For the MT model, heritability estimates ranged from 0.20 (37 months of age) to 0.32 (94 months of age). The genetic correlation estimates for MY305 obtained by RRM (L23.res4) and MT models were very similar, and varied from 0.77 to 0.99 and from 0.77 to 0.99, respectively. The rank correlation between breeding values for MY305 at different ages predicted by REP, MT, and RRM were high. It seems that a linear and quadratic Legendre polynomial to model the additive genetic and animal permanent environmental effects, respectively, may be sufficient to explain more parsimoniously the changes in MY305 genetic variation with age.