The marker effects of a single-step random regression model for 4 test-day traits in German Holsteins.

The single-step genomic model has become the golden standard for routine evaluation in livestock species, such as Holstein dairy cattle. The single-step genomic model with direct estimation of marker effects has been proven to be efficient in accurately accounting for millions of genotype records. For diverse applications including frequent genomic evaluation updates on a weekly basis, estimates of the marker effects from the single-step evaluations play a central role in genomic prediction. In this study we focused on exploring the marker effect estimates from the single-step evaluation. Phenotypic, genotypic, and pedigree data were taken from the official evaluation for German dairy breeds in April 2021. A multilactation random regression test-day model was applied to more than 242 million test-day records separately for 4 traits: milk, fat, and protein yields, and somatic cell scores (SCS). Approximately one million genotyped Holstein animals were considered in the single-step genomic evaluations including ∼21 million animals in pedigree. Deregressed multiple across-country breeding values of Holstein bulls having daughters outside Germany were integrated into the national test-day data to increase the reliability of genomic breeding values. To assess the stability and bias of the marker effects of the single-step model, test-day records of the last 4 yr were deleted, and the integrated bulls born in the last 4 yr were truncated from the complete phenotypic dataset. Estimates of the marker effects were shown to be highly correlated, with correlations ∼0.9, between the full and truncated evaluations. Regression slope values of the marker-effect estimates from the full on the truncated evaluations were all close to their expected value, being ∼1.03. Calculated using random regression coefficients of the marker effect estimates, drastically different shapes of the genetic lactation curve were seen for 2 markers on chromosome 14 for the 4 test-day traits. The contribution of individual chromosomes to the total additive genetic variances seemed to follow the polygenic inheritance mode for protein yield and SCS. However, chromosome 14 was found to make an exceptionally large contribution to the total additive genetic variance for milk and fat yields because of markers near the major gene DGAT1. For the first lactation test-day traits, we obtained ∼0 correlations of chromosomal direct genomic values between any pair of the chromosomes; no spurious correlations were found in our analysis, thanks to the large reference population. For trait milk yield, chromosomal direct genomic values appeared to have a large variation in the between-lactation correlations among the chromosomes, especially between first and second or third lactations. The optimal features of the random regression test-day model and the single-step marker model allowed us to track the differences in the shapes of genetic lactation curves down to the individual markers. Furthermore, the single-step random regression test-day model enabled us to better understand the inheritance mode of the yield traits and SCS (e.g., variable chromosomal contributions to the total additive genetic variance and to the genetic correlations between lactations).

Estimation of the threshold for heat stress and genetic features for milk yield in Mehsana buffaloes in India.

Heat stress is one of the primary environmental factors that harm both the productivity and health of buffaloes. The current study was conducted to estimate the threshold of temperature humidity index (THI)1 and genetic features for milk yield of first-lactation Mehsana buffaloes using an univariate repeatability test-day model. The data included 130,475 first lactation test-day milk yield (FLTDMY) records of 13,887 Mehsana buffaloes and the daily temperature and humidity. The statistical model included herd test day as fixed effects, days-in-milk (DIM) classes, age of the animal, as well as random factors such as the additive genetic effect (AGE) of animal in general conditions (intercept), AGE of the buffaloes subjected to heat stress (slope), permanent environmental effect of animal in general conditions (intercept), permanent environmental effect of animal under heat stress conditions (slope) and random residual effect. It was expected that the general effects and the heat-tolerance effects would be correlated, represented by the present investigation's repeatability models. The variance components of FLTDMY in the present study were computed using the REML method. The threshold for THI was 78. At the THI below the threshold, the heritability estimated for the FLTDMY trait was 0.29, and the additive genetic variance (AGV) for heat stress conditions was 0. At THI of 83, AGV for heat stress conditions was highest for FLTDMY. The genetic correlation of general AGE to heat-tolerant AGE was -0.40. The results indicated that a consistent selection for milk production, avoiding the thermal tolerance, may diminish the thermal tolerance capacity of Mehsana buffaloes.

Genetic relationships between weight loss in early lactation and daily milk production throughout lactation in Holstein cows.

After calving, high-yielding dairy cows mobilize body reserves for energy, sometimes to the detriment of health and fertility. This study aimed to estimate the genetic correlation between body weight loss until nadir and daily milk production (MY24) in first- (L1) and second-lactation (L2) Holstein cows. The data set included 859,020 MY24 records and 570,651 daily raw body weight (BWr) phenotypes from 3,989 L1 cows, and 665,361 MY24 records and 449,449 BWr phenotypes from 3,060 L2 cows, recorded on 36 French commercial farms equipped with milking robots that included an automatic weighing platform. To avoid any bias due to change in digestive content, BWr was adjusted for variations in feed intake, estimated from milk production and BWr. Adjusted body weight was denoted BW. The genetic parameters of BW and MY24 in L1 and L2 cows were estimated using a 4-trait random regression model. In this model, the random effects were fitted by second-order Legendre polynomials on a weekly basis from wk 1 to 44. Nadir of BW was found to be earlier than reported in the literature, at 29 d in milk, and BW loss from calving to nadir was also lower than generally assumed, close to 29 kg. To estimate genetic correlations between body weight loss and production, we defined BWL5 as the loss of weight between wk 1 and 5 after calving. Genetic correlations between BWL5 and MY24 ranged from -0.26 to 0.05 in L1 and from -0.11 to 0.10 in L2, according to days in milk. These moderate to low values suggest that it may be possible to select for milk production without increasing early body mobilization.

Predicting daily milk yield for primiparous cows using data of within-herd relatives to capture genotype-by-environment interactions.

This study develops and illustrates a hybrid k-medoids, random forest, and support vector regression (K-R-S) approach for predicting the lactation curves of individual primiparous cows within a targeted environment using monthly milk production data from their dams and paternal siblings. The model simulation and evaluation were based on historical test-day (TD) milk production data from 2010 to 2016 for 260 Wisconsin dairy farms. Data from older paternal siblings and dams were used to create family units (n = 6,400) of individual calves, from which their future performance was predicted. Test-day milk yield (MY) records from 2010 to 2014 were used for model training, whereas monthly milk production records of Holstein calves born in 2014 were used for model evaluation. The K-R-S hybrid approach was used to generate MY predictions for 5 randomly selected batches of 320 primiparous cows, which were used to evaluate model performance at the individual cow level by cross-validation. Across all 5 batches, the mean absolute error and the root mean square error of the K-R-S predictions were lower (by 24.2 and 23.4%, respectively) than that of the mean daily MY of paternal siblings. The K-R-S predictions of TD MY were closer to actual values 74.2 ± 2.0% of the time, as compared with means of paternal siblings'. The correlation between actual TD MY and K-R-S predictions was greater (0.34 ± 0.01) than the correlation between the actual yield and the mean of paternal siblings (0.08 ± 0.01). The results of this study demonstrate the effectiveness of the K-R-S hybrid approach for predicting future first-lactation MY of dairy calves in management applications, such as milk production forecasting or decision-support simulation, using only monthly TD yields of within-herd relatives and in the absence of detailed genomic data.

Genomic prediction and genetic correlations estimated for milk production and fatty acid traits in Walloon Holstein cattle using random regression models.

The aims of this study were to: (1) estimate genetic correlation for milk production traits (milk, fat and protein yields and fat and protein contents) and fatty acids (FA: C16:0, C18:1 cis-9, LCFA, SFA, and UFA) over days in milk, (2) investigate the performance of genomic predictions using single-step GBLUP (ssGBLUP) based on random regression models (RRM), and (3) identify the optimal scaling and weighting factors to be used in the construction of the H matrix. A total of 302 684 test-day records of 63.875 first lactation Walloon Holstein cows were used. Positive genetic correlations were found between milk yield and fat and protein yield (rg from 0.46 to 0.85) and between fat yield and milk FA (rg from 0.17 to 0.47). On the other hand, negative correlations were estimated between fat and protein contents (rg from -0.22 to -0.59), between milk yield and milk FA (rg from -0.22 to -0.62), and between protein yield and milk FA (rg from -0.11 to -0.19). The selection for high fat content increases milk FA throughout lactation (rg from 0.61 to 0.98). The test-day ssGBLUP approach showed considerably higher prediction reliability than the parent average for all milk production and FA traits, even when no scaling and weighting factors were used in the H matrix. The highest validation reliabilities (r2 from 0.09 to 0.38) and less biased predictions (b1 from 0.76 to 0.92) were obtained using the optimal parameters (i.e., ω = 0.7 and α = 0.6) for the genomic evaluation of milk production traits. For milk FA, the optimal parameters were ω = 0.6 and α = 0.6. However, biased predictions were still observed (b1 from 0.32 to 0.81). The findings suggest that using ssGBLUP based on RRM is feasible for the genomic prediction of daily milk production and FA traits in Walloon Holstein dairy cattle.

Genome-wide association for milk urea concentration in Dual-Purpose Belgian Blue cows.

The objectives of this study were to estimate genetic parameters and identify genomic regions associated with milk urea concentration (MU) in Dual-Purpose Belgian Blue (DPBB) cows. The data were 29,693 test-day records of milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FP), protein percentage (PP) and MU collected between 2014 and 2020 on 2498 first parity cows (16,935 test-day records) and 1939 second-parity cows (12,758 test-day records) from 49 herds in the Walloon Region of Belgium. Data of 28,266 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), on 1699 animals (639 males and 1060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method using a single chain of 100,000 iterations after a burn-in period of 20,000. SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by windows of 25 consecutive SNPs (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. The mean (SD) of MU was 22.89 (10.07) and 22.35 (10.07) mg/dl for first and second parity, respectively. The mean (SD) heritability estimates for daily MU were 0.18 (0.01) and 0.22 (0.02), for first and second parity, respectively. The mean (SD) genetic correlations between daily MU and MY, FY, PY, FP and PP were -0.05 (0.09), -0.07 (0.11), -0.03 (0.13), -0.05 (0.08) and -0.03 (0.11) for first parity, respectively. The corresponding values estimated for second parity were 0.02 (0.10), -0.02 (0.09), 0.02 (0.08), -0.08 (0.06) and -0.05 (0.05). The genome-wide association analyses identified three genomic regions (BTA2, BTA3 and BTA13) associated with MU. The identified regions showed contrasting results between parities and among different stages within each parity. This suggests that different groups of candidate genes underlie the phenotypic expression of MU between parities and among different lactation stages within a parity. The results of this study can be used for future implementation and use of genomic evaluation to reduce MU in DPBB cows.

Effect of a calcium-energy supplement drink at calving on lactation performance: Milk yield and composition, odds to reach a next lactation, and calving interval.

Supplementation of Ca products to cows after calving is common in calving protocols. This study evaluated the effect of a Ca-energy drink voluntarily consumed on milk yield and composition, odds to reach a next lactation, and calving interval. This prospective randomized study included a blinded placebo and was conducted in 10 commercial dairy farms that included 504 Holstein dairy cows. Cows were blocked within farm by calving sequence and parity (primiparous or multiparous). Within each block of 2 animals, cows were randomly assigned to 1 of 2 treatments: a Ca-energy supplement drink (CAE, n = 255) providing 45 g of Ca and other components (dextrose, lactose, protein, fat, other minerals and vitamins), a placebo (i.e., 100 g of cellulose and 20 g of dextrose; CON, n = 249), both strictly offered to the animals for voluntary consumption. Treatments were offered mixed in 20 L of water within 3 h after calving. Milk data were analyzed using 2 approaches. The first, most classical, evaluated the effect of the treatments on observed milk data, whereas the second approach evaluated the effect on milk residuals (i.e., the difference between observed milk data and a prediction made by a herd test-day model). Eighty-one percent of the CAE cows fully consumed the treatment, whereas only 50% of CON cows did. No differences were detected for observed milk yield, nor for composition in multiparous cows. The only production effect observed on multiparous cows was a treatment by time interaction for milk fat yield, reflecting greater yield for CAE cows between 100 and 150 d in milk only. However, primiparous cows receiving CAE had increased milk (+0.8 kg/d) and component yields (i.e., +40 g/d of protein) compared with CON cows. These effects were more evident when milk and milk components residuals data were analyzed (i.e., +1.5 kg/d for milk yield and +57 g/d of protein). This was achieved with a herd test-day model that allowed milk and milk components data to be adjusted for environmental and genetic factors (i.e., farm effect, time effect, age at calving, parity, stage of lactation, breeding value). The treatment had no effect on the probability of reaching the next lactation (i.e., 72% of CAE cows had a next calving against 69% in CON). Primiparous cows receiving CAE had a longer calving interval compared with CON cows. At 400 d after the application of the treatment, 65% of CAE primiparous cows had a next calving, whereas 81% of CON primiparous cows had calved already. The supplementation of the tested oral Ca-energy solution at calving did not increase the probability to reach a next lactation for neither primiparous or multiparous, but positively influenced milk yield and milk component yields for primiparous.

Autoregressive single-step test-day model for genomic evaluations of Portuguese Holstein cattle.

The multiple-lactation autoregressive test-day (AR) model is the adopted model for the national genetic evaluation of dairy cattle in Portugal. Under this model, animals' permanent environment effects are assumed to follow a first-order autoregressive process over the long (auto-correlations between parities) and short (auto-correlations between test-days within lactation) terms. Given the relevance of genomic prediction in dairy cattle, it is essential to include marker information in national genetic evaluations. In this context, we aimed to evaluate the feasibility of applying the single-step genomic (G)BLUP to analyze milk yield using the AR model in Portuguese Holstein cattle. In total, 11,434,294 test-day records from the first 3 lactations collected between 1994 and 2017 and 1,071 genotyped bulls were used in this study. Rank correlations and differences in reliability among bulls were used to compare the performance of the traditional (A-AR) and single-step (H-AR) models. These 2 modeling approaches were also applied to reduced data sets with records truncated after 2012 (deleting daughters of tested bulls) to evaluate the predictive ability of the H-AR. Validation scenarios were proposed, taking into account young and proven bulls. Average EBV reliabilities, empirical reliabilities, and genetic trends predicted from the complete and reduced data sets were used to validate the genomic evaluation. Average EBV reliabilities for H-AR (A-AR) using the complete data set were 0.52 (0.16) and 0.72 (0.62) for genotyped bulls with no daughters and bulls with 1 to 9 daughters, respectively. These results showed an increase in EBV reliabilities of 0.10 to 0.36 when genomic information was included, corresponding to a reduction of up to 43% in prediction error variance. Considering the 3 validation scenarios, the inclusion of genomic information improved the average EBV reliability in the reduced data set, which ranged, on average, from 0.16 to 0.26, indicating an increase in the predictive ability. Similarly, empirical reliability increased by up to 0.08 between validation tests. The H-AR outperformed A-AR in terms of genetic trends when unproven genotyped bulls were included. The results suggest that the single-step GBLUP AR model is feasible and may be applied to national Portuguese genetic evaluations for milk yield.

Variance Components and Genetic Parameters Estimated for Fat and Protein Content in Individual Months of Lactation: The Case of Tsigai Sheep.

The objective of this study was to assess variance components and genetic parameters for fat and protein content in Tsigai sheep using multivariate animal models in which fat and protein content in individual months of lactation were treated as different traits, and univariate models in which fat and protein content were treated as repeated measures of the same traits. Test day measurements were taken between the second and the seventh month of lactation. The fixed effects were lactation number, litter size and days in milk. The random effects were animal genetic effect and permanent environmental effect of ewe. The effect of flock-year-month of test day measurement was fitted either as a fixed (FYM) or random (fym) effect. Heritabilities for fat content were estimated between 0.06 and 0.17 (FYM fitted) and between 0.06 and 0.11 (fym fitted). Heritabilities for protein content were estimated between 0.15 and 0.23 (FYM fitted) and between 0.10 and 0.18 (fym fitted). For fat content, variance ratios of permanent environmental effect of ewe were estimated between 0.04 and 0.11 (FYM fitted) and between 0.02 and 0.06 (fym fitted). For protein content, variance ratios of permanent environmental effect of ewe were estimated between 0.13 and 0.20 (FYM fitted) and between 0.08 and 0.12 (fym fitted). The proportion of phenotypic variance explained by fym effect ranged from 0.39 to 0.43 for fat content and from 0.25 to 0.36 for protein content. Genetic correlations between individual months of lactation ranged from 0.74 to 0.99 (fat content) and from 0.64 to 0.99 (protein content). Fat content heritabilities estimated with univariate animal models roughly corresponded with heritability estimates from multivariate models: 0.13 (FYM fitted) and 0.07 (fym fitted). Protein content heritabilities estimated with univariate animal models also corresponded with heritability estimates from multivariate models: 0.18 (FYM fitted) and 0.13 (fym fitted).

Short communication: Reliability of single-step genomic BLUP breeding values by multi-trait test-day model analysis.

The purpose of our study was to develop an approximation procedure to estimate reliabilities of single-step genomic BLUP breeding values in a test-day model for routine evaluation of milk yield in a dairy cattle population. Input data consisted of 20,220,047 first-, second-, and third-lactation test-day milk yield records of 1,126,102 Czech Holstein cows (each lactation being considered a separate trait), with 1,844,679 animals in the pedigree file and with genomic data from 2,236 bulls. Evaluation was according to a multi-lactation model. The procedure was based on the effective number of records per animal from milk recording as well as from genomic and pedigree relationships. Traits were analyzed individually, and genetic covariances among traits were subsequently taken into account. The use of genomic information increased average reliability in young bulls from 0.276 to 0.505, but increased reliability in proven bulls only from 0.828 to 0.855. The reliabilities of genomic breeding values in multi-trait evaluation for first, second and third lactations, respectively, averaged 0.652, 0.673, and 0.633 for young bulls and 0.907, 0.894, and 0.852 for proven bulls. For an index combining all 3 lactations, the average reliability of a single-step genomic BLUP prediction was 0.712 and 0.925 for younger and proven bulls, respectively. Increased reliability due to genotyping in the population of all genotyped and nongenotyped animals was very small (<0.01) because of the small proportion of genotyped animals in the population.

Single-step genomic evaluation using multitrait random regression model and test-day data.

The objectives of this study were to evaluate the feasibility of use of the test-day (TD) single-step genomic BLUP (ssGBLUP) using phenotypic records of Nordic Red Dairy cows. The critical point in ssGBLUP is how genomically derived relationships (G) are integrated with population-based pedigree relationships (A) into a combined relationship matrix (H). Therefore, we also tested how different weights for genomic and pedigree relationships affect ssGBLUP, validation reliability, and validation regression coefficients. Deregressed proofs for 305-d milk, protein, and fat yields were used for a posteriori validation. The results showed that the use of phenotypic TD records in ssGBLUP is feasible. Moreover, the TD ssGBLUP model gave considerably higher validation reliabilities and validation regression coefficients than the TD model without genomic information. No significant differences were found in validation reliability between the different TD ssGBLUP models according to bootstrap confidence intervals. However, the degree of inflation in genomic enhanced breeding values is affected by the method used in construction of the H matrix. The results showed that ssGBLUP provides a good alternative to the currently used multi-step approach but there is a great need to find the best option to combine pedigree and genomic information in the genomic matrix.

Genetic Analysis of Milk Yield in First-Lactation Holstein Friesian in Ethiopia: A Lactation Average vs Random Regression Test-Day Model Analysis.

The development of effective genetic evaluations and selection of sires requires accurate estimates of genetic parameters for all economically important traits in the breeding goal. The main objective of this study was to assess the relative performance of the traditional lactation average model (LAM) against the random regression test-day model (RRM) in the estimation of genetic parameters and prediction of breeding values for Holstein Friesian herds in Ethiopia. The data used consisted of 6,500 test-day (TD) records from 800 first-lactation Holstein Friesian cows that calved between 1997 and 2013. Co-variance components were estimated using the average information restricted maximum likelihood method under single trait animal model. The estimate of heritability for first-lactation milk yield was 0.30 from LAM whilst estimates from the RRM model ranged from 0.17 to 0.29 for the different stages of lactation. Genetic correlations between different TDs in first-lactation Holstein Friesian ranged from 0.37 to 0.99. The observed genetic correlation was less than unity between milk yields at different TDs, which indicated that the assumption of LAM may not be optimal for accurate evaluation of the genetic merit of animals. A close look at estimated breeding values from both models showed that RRM had higher standard deviation compared to LAM indicating that the TD model makes efficient utilization of TD information. Correlations of breeding values between models ranged from 0.90 to 0.96 for different group of sires and cows and marked re-rankings were observed in top sires and cows in moving from the traditional LAM to RRM evaluations.

Association of total mixed ration particle fractions retained on the Penn State Particle Separator with milk, fat, and protein yield lactation curves at the cow level.

As part of a larger project aiming to develop management evaluation tools based on results from test-day (TD) models, the objective of this study was to examine the effect of physical composition of total mixed rations (TMR) tested quarterly from March 2006 through December 2008 on milk, fat, and protein yield curves for 25 herds in Ragusa, Sicily. A random regression sire-maternal grandsire model was used to estimate variance components for milk, fat, and protein yields fitted on a full data set, including 241,153 TD records from 9,809 animals in 42 herds recorded from 1995 through 2008. The model included parity, age at calving, year at calving, and stage of pregnancy as fixed effects. Random effects were herd × test date, sire and maternal grandsire additive genetic effect, and permanent environmental effect modeled using third-order Legendre polynomials. Model fitting was carried out using ASREML. Afterward, for the 25 herds involved in the study, 9 particle size classes were defined based on the proportions of TMR particles on the top (19-mm) and middle (8-mm) screen of the Penn State Particle Separator. Subsequently, the model with estimated variance components was used to examine the influence of TMR particle size class on milk, fat, and protein yield curves. An interaction was included with the particle size class and days in milk. The effect of the TMR particle size class was modeled using a ninth-order Legendre polynomial. Lactation curves were predicted from the model while controlling for TMR chemical composition (crude protein content of 15.5%, neutral detergent fiber of 40.7%, and starch of 19.7% for all classes), to have pure estimates of particle distribution not confounded by nutrient content of TMR. We found little effect of class of particle proportions on milk yield and fat yield curves. Protein yield was greater for sieve classes with 10.4 to 17.4% of TMR particles retained on the top (19-mm) sieve. Optimal distributions different from those recommended may reflect regional differences based on climate and types and quality of forages fed.

Random Regression Model for Genetic Evaluation and Early Selection in the Iranian Holstein Population.

The objective of this study was to use a model to predict breeding values for sires and cows at an early stage of the first lactation of cows and progeny groups in the Iranian Holstein population to enable the early selection of sires. An additional objective was to estimate genetic and phenotypic parameters associated with this model. The accuracy of predicted breeding values was investigated using cross-validation based on sequential genetic evaluations emulating yearly evaluation runs. The data consisted of 2,166,925 test-day records from 456,712 cows calving between 1990 and 2015. (Co)-variance components and breeding values were estimated using a random regression test-day model and the average information (AI) restricted maximum likelihood method (REML). Legendre polynomial functions of order three were chosen to fit the additive genetic and permanent environmental effects, and a homogeneous residual variance was assumed throughout lactation. The lowest heritability of daily milk yield was estimated to be just under 0.14 in early lactation, and the highest heritability of daily milk yield was estimated to be 0.18 in mid-lactation. Cross-validation showed a highly positive correlation of predicted breeding values between consecutive yearly evaluations for both cows and sires. Correlation between predicted breeding values based only on records of early lactation (5-90 days) and records including late lactation (181-305 days) were 0.77-0.87 for cows and 0.81-0.94 for sires. These results show that we can select sires according to their daughters' early lactation information before they finish the first lactation. This can be used to decrease generation interval and to increase genetic gain in the Iranian Holstein population.

Genetic Parameter Estimation and Genome-Wide Association Study-Based Loci Identification of Milk-Related Traits in Chinese Holstein.

Accurately estimating the genetic parameters and revealing more genetic variants underlying milk production and quality are conducive to the genetic improvement of dairy cows. In this study, we estimate the genetic parameters of five milk-related traits of cows-namely, milk yield (MY), milk fat percentage (MFP), milk fat yield (MFY), milk protein percentage (MPP), and milk protein yield (MPY)-based on a random regression test-day model. A total of 95,375 test-day records of 9,834 cows in the lower reaches of the Yangtze River were used for the estimation. In addition, genome-wide association studies (GWASs) for these traits were conducted, based on adjusted phenotypes. The heritability, as well as the standard errors, of MY, MFP, MFY, MPP, and MPY during lactation ranged from 0.22 ± 0.02 to 0.31 ± 0.04, 0.06 ± 0.02 to 0.15 ± 0.03, 0.09 ± 0.02 to 0.28 ± 0.04, 0.07 ± 0.01 to 0.16 ± 0.03, and 0.14 ± 0.02 to 0.27 ± 0.03, respectively, and the genetic correlations between different days in milk (DIM) within lactations decreased as the time interval increased. Two, six, four, six, and three single nucleotide polymorphisms (SNPs) were detected, which explained 5.44, 12.39, 8.89, 10.65, and 7.09% of the phenotypic variation in MY, MFP, MFY, MPP, and MPY, respectively. Ten Kyoto Encyclopedia of Genes and Genomes pathways and 25 Gene Ontology terms were enriched by analyzing the nearest genes and genes within 200 kb of the detected SNPs. Moreover, 17 genes in the enrichment results that may play roles in milk production and quality were selected as candidates, including CAMK2G, WNT3A, WNT9A, PLCB4, SMAD9, PLA2G4A, ARF1, OPLAH, MGST1, CLIP1, DGAT1, PRMT6, VPS28, HSF1, MAF1, TMEM98, and F7. We hope that this study will provide useful information for in-depth understanding of the genetic architecture of milk production and quality traits, as well as contribute to the genomic selection work of dairy cows in the lower reaches of the Yangtze River.

Development of genetic evaluation for milk production traits of Holsteins in Japan.

The procedure used for the genetic evaluation of dairy cattle in Japan has developed from a lactation sire-MGS model to a multiple-lactation random regression test-day animal model. Genetic evaluation of Holstein bulls in Japan began in 1989 with the use of field-style progeny testing; dairy herd improvement program data from all over Japan were used, along with a sire and maternal grandsire model. In 1993, an animal model was introduced to estimate breeding values for yield and type traits. A random regression test-day model was first applied in 2010. In the business of breeding dairy cattle, it is very important to users that estimated breeding values are reliable and stable among subsequent routine evaluations. With experience in the genetic evaluation of dairy cattle in Japan, Japanese researchers have found ways to improve the stability of estimated breeding values. These modifications involve changes in data editing, development of evaluation models, changes to the structures of unknown-parent groups, awareness of the problems of predicting lactation yield from partial test-day records, and adjustment for heterogeneity within herd variances. Here, I introduce developments in, and our experiences with, the genetic evaluation of yield traits of Holstein cattle in Japan.

Effects of stage of pregnancy on variance components, daily milk yields and 305-day milk yield in Holstein cows, as estimated by using a test-day model.

Pregnancy and calving are elements indispensable for dairy production, but the daily milk yield of cows decline as pregnancy progresses, especially during the late stages. Therefore, the effect of stage of pregnancy on daily milk yield must be clarified to accurately estimate the breeding values and lifetime productivity of cows. To improve the genetic evaluation model for daily milk yield and determine the effect of the timing of pregnancy on productivity, we used a test-day model to assess the effects of stage of pregnancy on variance component estimates, daily milk yields and 305-day milk yield during the first three lactations of Holstein cows. Data were 10 646 333 test-day records for the first lactation; 8 222 661 records for the second; and 5 513 039 records for the third. The data were analyzed within each lactation by using three single-trait random regression animal models: one model that did not account for the stage of pregnancy effect and two models that did. The effect of stage of pregnancy on test-day milk yield was included in the model by applying a regression on days pregnant or fitting a separate lactation curve for each days open (days from calving to pregnancy) class (eight levels). Stage of pregnancy did not affect the heritability estimates of daily milk yield, although the additive genetic and permanent environmental variances in late lactation were decreased by accounting for the stage of pregnancy effect. The effects of days pregnant on daily milk yield during late lactation were larger in the second and third lactations than in the first lactation. The rates of reduction of the 305-day milk yield of cows that conceived fewer than 90 days after the second or third calving were significantly (P<0.05) greater than that after the first calving. Therefore, we conclude that differences between the negative effects of early pregnancy in the first, compared with later, lactations should be included when determining the optimal number of days open to maximize lifetime productivity in dairy cows.

Domestic estimated breeding values and genomic enhanced breeding values of bulls in comparison with their foreign genomic enhanced breeding values.

Estimated breeding values (EBVs) and genomic enhanced breeding values (GEBVs) for milk production of young genotyped Holstein bulls were predicted using a conventional BLUP - Animal Model, a method fitting regression coefficients for loci (RRBLUP), a method utilizing the realized genomic relationship matrix (GBLUP), by a single-step procedure (ssGBLUP) and by a one-step blending procedure. Information sources for prediction were the nation-wide database of domestic Czech production records in the first lactation combined with deregressed proofs (DRP) from Interbull files (August 2013) and domestic test-day (TD) records for the first three lactations. Data from 2627 genotyped bulls were used, of which 2189 were already proven under domestic conditions. Analyses were run that used Interbull values for genotyped bulls only or that used Interbull values for all available sires. Resultant predictions were compared with GEBV of 96 young foreign bulls evaluated abroad and whose proofs were from Interbull method GMACE (August 2013) on the Czech scale. Correlations of predictions with GMACE values of foreign bulls ranged from 0.33 to 0.75. Combining domestic data with Interbull EBVs improved prediction of both EBV and GEBV. Predictions by Animal Model (traditional EBV) using only domestic first lactation records and GMACE values were correlated by only 0.33. Combining the nation-wide domestic database with all available DRP for genotyped and un-genotyped sires from Interbull resulted in an EBV correlation of 0.60, compared with 0.47 when only Interbull data were used. In all cases, GEBVs had higher correlations than traditional EBVs, and the highest correlations were for predictions from the ssGBLUP procedure using combined data (0.75), or with all available DRP from Interbull records only (one-step blending approach, 0.69). The ssGBLUP predictions using the first three domestic lactation records in the TD model were correlated with GMACE predictions by 0.69, 0.64 and 0.61 for milk yield, protein yield and fat yield, respectively.

Genetic parameters for milk, fat and protein yields in Murrah buffaloes (Bubalus bubalis Artiodactyla, Bovidae).

The objective of the present study was to estimate genetic parameters for test-day milk, fat and protein yields and 305-day-yields in Murrah buffaloes. 4,757 complete lactations of Murrah buffaloes were analyzed. Co-variance components were estimated by the restricted maximum likelihood method. The models included additive direct genetic and permanent environmental effects as random effects, and the fixed effects of contemporary group, milking number and age of the cow at calving as linear and quadratic covariables. Contemporary groups were defined by herd-year-month of test for test-day yields and by herd-year-season of calving for 305-day yields. The heritability estimates obtained by two-trait analysis ranged from 0.15 to 0.24 for milk, 0.16 to 0.23 for protein and 0.13 to 0.22 for fat, yields. Genetic and phenotypic correlations were all positive. The observed population additive genetic variation indicated that selection might be an effective tool in changing population means in milk, fat and protein yields.

Genetic parameters for milk, fat and protein yields in Murrah buffaloes (Bubalus bubalis Artiodactyla, Bovidae)

The objective of the present study was to estimate genetic parameters for test-day milk, fat and protein yields and 305-day-yields in Murrah buffaloes. 4,757 complete lactations of Murrah buffaloes were analyzed. Co-variance components were estimated by the restricted maximum likelihood method. The models included additive direct genetic and permanent environmental effects as random effects, and the fixed effects of contemporary group, milking number and age of the cow at calving as linear and quadratic covariables. Contemporary groups were defined by herd-year-month of test for test-day yields and by herd-year-season of calving for 305-day yields. The heritability estimates obtained by two-trait analysis ranged from 0.15 to 0.24 for milk, 0.16 to 0.23 for protein and 0.13 to 0.22 for fat, yields. Genetic and phenotypic correlations were all positive. The observed population additive genetic variation indicated that selection might be an effective tool in changing population means in milk, fat and protein yields.