Test-day somatic cell score, fat-to-protein ratio and milk yield as indicator traits for sub-clinical mastitis in dairy cattle

Test-day (TD) records of milk, fat-to-protein ratio (F:P) and somatic cell score (SCS) of first-lactation Canadian Holstein cows were analysed by a three-trait finite mixture random regression model, with the purpose of revealing hidden structures in the data owing to putative, sub-clinical mastitis. Different distributions of the data were allowed in 30 intervals of days in milk (DIM), covering the lactation from 5 to 305 days. Bayesian analysis with Gibbs sampling was used for model inferences. Estimated proportion of TD records originated from cows infected with mastitis was 0.66 in DIM from 5 to 15 and averaged 0.2 in the remaining part of lactation. Data from healthy and mastitic cows exhibited markedly different distributions, with respect to both average value and the variance, across all parts of lactation. Heterogeneity of distributions for infected cows was also apparent in different DIM intervals. Cows with mastitis were characterized by smaller milk yield (down to -5 kg) and larger F:P (up to 0.13) and SCS (up to 1.3) compared with healthy contemporaries. Differences in averages between healthy and infected cows for F:P were the most profound at the beginning of lactation, when a dairy cow suffers the strongest energy deficit and is therefore more prone to mammary infection. Residual variances for data from infected cows were substantially larger than for the other mixture components. Fat-to-protein ratio had a significant genetic component, with estimates of heritability that were larger or comparable with milk yield, and was not strongly correlated with milk and SCS on both genetic and environmental scales. Daily milk, F:P and SCS are easily available from milk-recording data for most breeding schemes in dairy cattle. Fat-to-protein ratio can potentially be a valuable addition to SCS and milk yield as an indicator trait for selection against mastitis.     

Application of multiple‐trait finite mixture model to test‐day records of milk yield and somatic cell score of Canadian Holsteins

Multiple‐trait (MT) finite mixture random regression (MIX) model was applied using Bayesian methods to first lactation test‐day (TD) milk yield and somatic cell score (SCS) of Canadian Holsteins, allowing for heterogeneity of distributions with respect to days in milk (DIM) in lactation. The assumption was that the associations between patterns of variation in these traits and mastitis would allow revealing the hidden structure in the data distribution because of unknown health status of cows. The MIX model assumed separate means and residual co‐variance structures for two components in four intervals of lactation, in addition to fitting the fixed effect of herd‐test‐day, and fixed and random regressions with Legendre polynomials. Results indicated that the mixture model was superior to standard MT model, as supported by the Bayes factor. Approximately 20% of TD records were classified as originated from cows with a putative, sub‐clinical form of mastitis. The proportion of records from mastitic cows was the largest at the beginning of lactation. The MIX model exhibited different distributions of data from healthy and infected cows in different parts of lactation. Records from sick cows were characterized by larger (smaller) means for SCS (milk) and larger variances. Residual, and daily genetic and environmental correlations between milk and SCS were smaller from the MIX model when compared with MT estimates. Heritabilities of both traits differed significantly among records from healthy, sick and MT model estimates. Both models fitted milk records from healthy cows relatively well. The ability of the MT model in handling SCS records, measured by model residuals, was low, but improved substantially, however, where the data were allowed to be separated into two components in the MIX parameterization. Correlations between estimated breeding values (EBV) for sires from both models were very high for cumulative milk yield (>0.99) and slightly lower (0.95 in the interval from 5 to 45 DIM) for daily SCS. EBV for SCS from MT and MIX models were weakly correlated with posterior probability of sub‐clinical mastitis on the phenotypic scale.     

Genetic analysis of fat‐to‐protein ratio,milk yield and somatic cell score of Holstein cows in Japan in the first three lactations by using a random regression model

We estimated the genetic parameters of fat‐to‐protein ratio (FPR) and the genetic correlations between FPR and milk yield or somatic cell score in the first three lactations in dairy cows. Data included 3 079 517 test‐day records of 201 138 Holstein cows in Japan from 2006 to 2011. Genetic parameters were estimated with a multiple‐trait random regression model in which the records within and between parities were treated as separate traits. The phenotypic values of FPR increased soon after parturition and peaked at 10 to 20 days in milk, then decreased slowly in mid‐ and late lactation. Heritability estimates for FPR yielded moderate values. Genetic correlations of FPR among parities were low in early lactation. Genetic correlations between FPR and milk yield were positive and low in early lactation, but only in the first lactation. Genetic correlations between FPR and somatic cell score were positive in early lactation and decreased to become negative in mid‐ to late lactation. By using these results for genetic evaluation it should be possible to improve energy balance in dairy cows.     

Genetic evaluation of Australian dairy cattle for somatic cell scores using multi-trait random regression test-day model

Summary A multi-trait (MT) random regression (RR) test day (TD) model has been developed for genetic evaluation of somatic cell scores for Australian dairy cattle, where first, second and third lactations were considered as three different but correlated traits. The model includes herd-test-day, year-season, age at calving, heterosis and lactation curves modelled with Legendre polynomials as fixed effects, and random genetic and permanent environmental effects modelled with Legendre polynomials. Residual variance varied across the lactation trajectory. The genetic parameters were estimated using asreml . The heritability estimates ranged from 0.05 to 0.16. The genetic correlations between lactations and between test days within lactations were consistent with most of the published results. Preconditioned conjugate gradient algorithm with iteration on data was implemented for solving the system of equations. For reliability approximation, the method of Tier and Meyer was used. The genetic evaluation system was validated with Interbull validation method III by comparing proofs from a complete evaluation with those from an evaluation based on a data set excluding the most recent 4 years. The genetic trend estimate was in the allowed range and correlations between the two sets of proofs were very high. Additionally, the RR model was compared to the previous test day model. The correlations of proofs between both models were high (0.97) for bulls with high reliabilities. The correlations of bulls decreased with increasing incompleteness of daughter performance information. The correlations between the breeding values from two consecutive runs were high ranging from 0.97 to 0.99. The MT RR TD model was able to make effective use of available information on young bulls and cows, and could offer an opportunity to breeders to utilize estimated breeding values for first and later lactations.     

Length of lags in responses of milk yield and somatic cell score on test day to heat stress in Holsteins

We used daily records from provincial Japanese weather stations and monthly test‐day records of milk production to investigate the length of the lags in the responses of cows’ milk yield and somatic cell score (SCS) to heat stress (HS). We also investigated the HS thresholds in milk yield and SCS. Data were a total of 17,245,709 test‐day records for milk and SCS in Holstein cows that had calved for the first time between 2000 and 2015, along with weather records from 60 weather stations. Temperature–humidity index (THI) values were estimated by using average daily temperature and average daily relative humidity. Adjusted THI values were calculated by using temperature, relative humidity, wind speed, and solar radiation. The model contained herd, calving year, month of test day, age group, days in milk, and THI as a fixed effect. THIs for each day from 14 days before the test day until the test day were used to represent the HS effects. The HS occurring 3 days, and between 8 and 10 days, before the test day had the greatest effect on the milk yield and SCS, respectively. The threshold THI values for the HS effect were about 60–65 for both traits.     

Genetic parameters for total lactation yields of milk,fat, protein,and somatic cell score in New Zealand dairy goats

The aim of this study was to estimate genetic parameters for lactation yields of milk (MY), fat (FY), protein (PY), and somatic cell score (SCS) of New Zealand dairy goats. The analysis used 64,604 lactation records from 23,583 does, kidding between 2004 and 2017, distributed in 21 flocks and representing 915 bucks. Estimates of genetic and residual (co) variances, heritabilities, and repeatabilities were obtained using a multiple‐trait repeatability animal model. The model included the fixed effects of contemporary group (does kidding in the same flock and year), age of the doe (in years), and as covariates, kidding day, proportion of Alpine, Nubian, Toggenburg, and “unknown” breeds (Saanen was used as the base breed), and heterosis. Random effects included additive animal genetic and doe permanent environmental effects. Estimates of heritabilities were 0.25 for MY, 0.24 for FY, 0.24 for PY, and 0.21 for SCS. The phenotypic correlations between MY, FY, and PY ranged from 0.90 to 0.96, and the genetic correlations ranged from 0.81 to 0.93. These results indicate lactation yield traits exhibit useful heritable variation and that multiple trait selection for these traits could improve milk revenue produced from successive generations of New Zealand dairy goats.     

A robust method for simultaneous estimation of single gene and polygenic effects in dairy cows using externally estimated breeding values as prior information

The aim of this study was to develop a robust method to estimate single gene and random polygenic animal effects simultaneously in a small field dataset with limited pedigree information. The new method was based on a Bayesian approach using additional prior information on the distribution of externally estimated breeding values. The field dataset consisted of 40 269 test‐day records for milk performance traits for 1455 genotyped dairy cows for the 11 bp‐deletion in the coding sequence of the myostatin gene. For all traits, estimated additive effects of the favoured wild‐type allele (‘+’ allele) were smaller when applying the new method in comparison with the application of a conventional mixed inheritance test‐day model. Dominance effects of the myostatin gene showed the same behaviour but were generally lower than additive effects. Robustness of methods was tested using a data‐splitting technique, based on the correlation of estimated breeding values from two samples, with one‐half of the data eliminated randomly from the first sample and the remaining data eliminated from the second sample. Results for 100 replicates showed that the correlation between split datasets when prior information included was higher than the conventional method. The new method led to more robust estimations for genetic effects and therefore has potential for use when only a small number of genotyped animals with field data and limited pedigree information are available.