Genetic parameters for production,health, fertility and longevity traits in dairy cows

Milk production, fertility, longevity and health records, were extracted from databases of two milk recording organisations in the United Kingdom for the first three lactations of the Holstein–Friesian breed. These included data related to health events (mastitis and lameness), voluntarily recorded on a proportion of farms. The data were analysed to calculate disease incidence levels and to estimate genetic parameters for health traits and their relationships with production and other functional traits. The resulting dataset consisted of 124 793 lactations from 75 137 animals of 1586 sires, recorded in 2434 herds. Incidence of health events increased with parity. The overall incidence of mastitis (MAS) and lameness (LAM), defined as binary traits, were 17% and 16%, respectively. Heritability estimates for MAS and LAM were 0.04 and 0.02, respectively, obtained from repeatability linear sire models. Heritability estimates of mastitis and lameness as count traits were slightly higher, 0.05 and 0.03, respectively. Genetic correlations were obtained by bivariate analyses of all pair-wise combinations between milk 305-day yield (MY), protein 305-day yield (PY), fat 305-day yield (FY), lactation average log_e transformed lactation average somatic cell count (SCS), calving interval (CI), days to first service (DFS), non-return at 56 days (NR56), number of inseminations (NINS), mastitis (MAS), number of mastitis episodes (NMAS), lameness (LAM), number of lameness episodes (NLAM) and lifespan score (LS). As expected, MAS was correlated most strongly with SCS (0.69), which supports the use of SCS as an indicator trait for mastitis. Genetic correlations between MAS and yield and fertility traits were of similar magnitude ranging from 0.27 to 0.33. Genetic correlations between MAS with LAM and LS were 0.38 and −0.59, respectively. Not all genetic correlations between LAM and other traits were significant because of fewer numbers of lameness records. LAM had significant genetic correlations with MY (0.38), PY (0.28), CI (0.35), NINS (0.38) and LS (−0.53). The heritability estimates of mastitis and lameness were low; therefore, genetic gain through direct selection alone would be slow, yet still positive and cumulative. Direct selection against mastitis and lameness as additional traits should reduce incidence of both diseases, and simultaneously improve fertility and longevity. However, both health traits had antagonistic relationships with production traits, thus genetic gain in production would be slower.     

Estimates of correlations among yield traits and somatic cell score with different models to adjust for bovine somatotropin effects on Holstein dairy cows

Records of Holstein cows from the Dairy Records Processing Center at Raleigh, NC were edited to obtain three data sets: 65,720 first, 50,694 second, and 65,445 later lactations. Correlations among yield traits and somatic cell score were estimated with three different models: 1) bovine somatotropin (bST) administration ignored, 2) bST administration as a fixed effect and 3) administration of bST as part of the contemporary group (herd-year-month-bST). Heritability estimates ranged from 0.13 to 0.17 for milk, 0.12 to 0.20 for fat, 0.14 to 0.16 for protein yields, and 0.08 to 0.09 for somatic cell score. Estimates were less for later than first lactations. Estimates of genetic correlations among yields ranged from 0.35 to 0.85 with no important differences between estimates with the 3 models. Estimates for lactation 2 agreed with estimates for lactation 1. Estimates of genetic correlations for later lactations were generally greater than for lactations 1 and 2 except between milk and protein yields. Estimates of genetic correlations between yields and somatic cell score were mostly negative or small (-0.45 to 0.11). Estimates of environmental correlations among yield traits were similar with all models (0.77 to 0.97). Estimates of environmental correlations between yields and somatic cell score were negative (-0.22 to -0.14). Estimates of phenotypic correlations among yield traits ranged from 0.70 to 0.95. Estimates of phenotypic correlations between yields and somatic cell score were small and negative. For all three data sets and all traits, no important differences in estimates of genetic parameters were found for the two models that adjusted for bST and the model that did not.     

Genetic parameters of a biological lactation model: early lactation and secretion rate traits of dairy heifers

Early lactation parameters are difficult to estimate from commercial dairy records due to the small number of records available before the peak of production. A biological model of lactation was used with weekly milk records from a single Holstein herd to estimate these early lactation parameters and the secretion rate of milk from the average cell throughout lactation. A genetic analysis of the lactation curve parameters, calculated curve characteristics and secretion rate traits was undertaken. Early lactation traits were found to have little genetic variation and effectively zero heritability. Secretion rate traits for milk, protein, lactose and water were all moderately heritable and highly genetically correlated (>0.87) but fat secretion rate had lower genetic correlations with the other secretion rates. A similar pattern of correlations was seen between total lactation yield traits for fat, protein, lactose and water. The genetic correlations between the lactation curve traits and the secretion rate traits were calculated. Total milk yield, peak yield and maximum secretion potential were all highly correlated with milk, lactose and water secretion rates but less so with fat and protein secretion rates. In particular, fat secretion rate had a moderate to low genetic correlation with these lactation curve traits. Persistency of lactation was highly correlated with fat and protein secretion rates, more persistent lactations being associated with lower rates of secretion of these milk components. Similar levels of heritability were found, where trait genetic parameters were directly equivalent to those derived from the same dataset by random regression methods. However, by using a biological model of lactation to analyse lactation traits new insights into the biology of lactation are possible and ways to select cows on a range of lactation traits may be achieved.     

中国荷斯坦牛IL8基因遗传多态性与泌乳性状以及体细胞评分的关联

以上海某奶牛场30个公牛家系的610头中国荷斯坦牛为试验材料,采用聚合酶链式反应-单链构象多态性(PCR-SSCP)技术对Interleukin-8(IL8)基因的遗传多态性进行了分析,采用混合动物模型分析了IL8基因突变位点与测定日产奶量、测定日乳脂率、测定日乳蛋白率、305d校正产奶量、305d乳脂量、305d乳蛋白量及测定日体细胞评分7个性状的相关性,寻找可用于生产实际的分子标记。共检测到KK、KA和AA3种基因型,频率分别为0.187、0.451和0.362,等位基因K和A的频率分别为0.412和0.588。该位点突变对测定日产奶量、305d乳蛋白量、305d校正产奶量和305d乳脂量以及体细胞评分影响达到极显著水平(P0.01),对测定日乳蛋白率的影响达到显著水平(P0.05),对测定日乳脂率影响不显著(P0.05)。多重比较表明:KK基因型对测定日产奶量、305d校正产奶量、305d乳蛋白量和305d乳脂量极显著高于AA和KA基因型(P0.01)。KK基因型的体细胞评分(SCS)最小二乘均值极显著低于KA、AA基因型(P0.01)。对于测定日的乳蛋白率AA基因型显著低于KA、KK型(P0.05)。IL8基因遗传突变对中国荷斯坦牛泌乳性状和乳房炎抗性有较大的遗传效应,可用于中国荷斯坦牛的分子标记辅助选择。     

Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) with production traits in Canadian dairy cattle

Associations of two alleles of the bovine major histocompatibility complex DRB3 gene (BoLA-DRB3) with lowered somatic cell score (SCS) and occurrence of disease (BoLA-DRB3.2* 16 and *23, respectively) have previously been documented. The objective of this study was to evaluate potential relationships between BoLA-DRB3 alleles with production traits, namely 305-day milk, milk fat and milk protein yield, in a population of Canadian dairy cattle (Holstein, n = 835 and Jersey, n = 66) over the course of two lactations. No significant associations were detected between BoLA alleles and production traits in Jerseys. In Holsteins, alleles *16 and *23 also did not show associations with production traits but allele *8 was significantly associated with increased 305-day milk, fat and protein yields in the previous lactation (the lactation prior to immunization with a gram negative core antigen vaccine), and with increased protein production in the subsequent (with reference to the time of immunization) lactation. Allele *22 was associated with decreased milk and protein yield in both previous and subsequent lactations. Therefore, it can be concluded that increasing or decreasing the frequency of BoLA alleles *16 and *23 to reduce SCS or increase resistance to mastitis in this population would not have adverse effects on production in this population, and that certain BoLA alleles (*8 and *22) are associated with altered production traits in Canadian Holsteins.     

Genetic parameters for serial, automatically recorded milkability and its relationship to udder health in dairy cattle

Serial measurements of three milkability traits from two commercial dairy farms in Germany were used to estimate heritabilities and breeding values (BVs). Overall, 6352 cows in first, second and third lactations supplied 2 188 810 records based on daily values recorded from 1998 to 2003. Only the records between day 8 and day 305 after calving were considered. The estimated genetic correlations between different parities within the three milkability traits ranged from rg = 0.88 to 0.98, i.e. they were sufficiently high to warrant a repeatability model. The resulting estimated heritability coefficients were h2 = 0.42 for average milk flow, h2 = 0.56 for maximum milk flow and h2 = 0.38 for milking time. We analysed the genetic correlation between milkability and somatic cell score (SCS) and between milkability and the liability to mastitis, respectively, as the optimum milk flow for udder health is not well defined. There were 66 146 records with information on somatic cell count. Furthermore, 23 488 days of medical treatment for udder diseases were available, resulting in 2 600 302 days of observation in total. Heritabilities for the liability to mastitis, estimated with a test-day threshold model, were h2 = 0.19 and h2 = 0.13, depending on the data-recording period (first 50 days of lactation and first 305 days of lactation, respectively). With respect to the relationship between milkability and udder health, the results indicated a slight and linear correlation insofar as one can assume: the higher the milk flow, the worse the udder health. For this reason, bulls and cows with high BVs for milk flow should be excluded from breeding to avoid a deterioration of udder health. The establishment of a special data-recording scheme for functional traits such as milkability and mastitis on commercial dairy farms may be possible according to these results.     

Merging and characterising phenotypic data on conventional and rare traits from dairy cattle experimental resources in three countries

This study set out to demonstrate the feasibility of merging data from different experimental resource dairy populations for joint genetic analyses. Data from four experimental herds located in three different countries (Scotland, Ireland and the Netherlands) were used for this purpose. Animals were first lactation Holstein cows that participated in ongoing or previously completed selection and feeding experiments. Data included a total of 60 058 weekly records from 1630 cows across the four herds; number of cows per herd ranged from 90 to 563. Weekly records were extracted from the individual herd databases and included seven traits: milk, fat and protein yield, milk somatic cell count, liveweight, dry matter intake and energy intake. Missing records were predicted with the use of random regression models, so that at the end there were 44 weekly records, corresponding to the typical 305-day lactation, for each cow. A total of 23 different lactation traits were derived from these records: total milk, fat and protein yield, average fat and protein percentage, average fat-to-protein ratio, total dry matter and energy intake and average dry matter intake-to-milk yield ratio in lactation weeks 1 to 44 and 1 to 15; average milk somatic cell count in lactation weeks 1 to 15 and 16 to 44; average liveweight in lactation weeks 1 to 44; and average energy balance in lactation weeks 1 to 44 and 1 to 15. Data were subsequently merged across the four herds into a single dataset, which was analysed with mixed linear models. Genetic variance and heritability estimates were greater (P < 0.05) than zero for all traits except for average milk somatic cell count in weeks 16 to 44. Proportion of total phenotypic variance due to genotype-by-environment (sire-by-herd) interaction was not different (P > 0.05) from zero. When estimable, the genetic correlation between herds ranged from 0.85 to 0.99. Results suggested that merging experimental herd data into a single dataset is both feasible and sensible, despite potential differences in management and recording of the animals in the four herds. Merging experimental data will increase power of detection in a genetic analysis and augment the potential reference population in genome-wide association studies, especially of difficult-to-record traits.     

Genetic parameters for milk mineral content and acidity predicted by mid-infrared spectroscopy in Holstein–Friesian cows

The aim of the present study was to estimate genetic parameters for calcium (Ca), phosphorus (P) and titratable acidity (TA) in bovine milk predicted by mid-IR spectroscopy (MIRS). Data consisted of 2458 Italian Holstein−Friesian cows sampled once in 220 farms. Information per sample on protein and fat percentage, pH and somatic cell count, as well as test-day milk yield, was also available. (Co)variance components were estimated using univariate and bivariate animal linear mixed models. Fixed effects considered in the analyses were herd of sampling, parity, lactation stage and a two-way interaction between parity and lactation stage; an additive genetic and residual term were included in the models as random effects. Estimates of heritability for Ca, P and TA were 0.10, 0.12 and 0.26, respectively. Positive moderate to strong phenotypic correlations (0.33 to 0.82) existed between Ca, P and TA, whereas phenotypic weak to moderate correlations (0.00 to 0.45) existed between these traits with both milk quality and yield. Moderate to strong genetic correlations (0.28 to 0.92) existed between Ca, P and TA, and between these predicted traits with both fat and protein percentage (0.35 to 0.91). The existence of heritable genetic variation for Ca, P and TA, coupled with the potential to predict these components for routine cow milk testing, imply that genetic gain in these traits is indeed possible.     

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.     

Genetic parameters for milk urea concentration and milk traits in Polish Holstein-Friesian cows

Milk urea concentration (MU) used by dairy producers for management purposes can be affected by selection for milk traits. To assess this problem, genetic parameters for MU in Polish Holstein-Friesian cattle were estimated for the first three lactations. The genetic correlation of MU with milk production traits, lactose percentage, fat to protein ratio (FPR) and somatic cell score (SCS) were computed with two 5-trait random regression test-day models, separately for each lactation. Data used for estimation (159,044 daily observations) came from 50 randomly sampled herds. (Co)variance components were estimated with the Bayesian Gibbs sampling method. The coefficient of variation for MU in all three parities was high (40–41 %). Average daily heritabilities of MU were 0.22 for the first parity and 0.21 for the second and third lactations. Average genetic correlations for different days in milk in the first three lactations between MU and other traits varied. They were small and negative for protein percentage (from ?0.24 to ?0.11) and for SCS (from ?0.14 to ?0.09). The weakest genetic correlation between MU and fat percentage, and between MU and lactose percentage were observed (from ?0.10 to 0.10). Negative average genetic correlation with the fat to protein ratio was observed only in the first lactation (?0.14). Genetic correlations with yield traits were positive and ranged from low to moderate for protein (from 0.09 to 0.33), fat (from 0.16 to 0.35) and milk yield (from 0.20 to 0.42). These results suggest that the selection on yield traits and SCS tends to increase MU slightly.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

Identification of SNPs in interferon gamma, interleukin-22, and their receptors and associations with health and production-related traits in Canadian Holstein bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

Heritability and repeatability of milk lactose and its relationships with traditional milk traits,somatic cell score and freezing point in Holstein cows

Lactose percentage (LP) in milk is currently determined in most herd-testing schemes, and globally, it is usually routinely recorded in the framework of the official milk recording procedures. However, few studies have investigated the phenotypic and genetic variability of this component. Data used in the present paper consisted of 59 811 test-day records from 4355 Holstein cows in 266 herds. Heritabilities of LP and lactose yield (LY) were estimated through single-trait repeatability animal models, whereas genetic and phenotypic correlations of LP and LY with milk composition and production traits, somatic cell score and milk freezing point were estimated using bivariate models. Fixed effects included in the analyses were herd-test-date, season of calving, parity, stage of lactation and the interaction between parity and stage of lactation. Random effects were animal additive genetic, within and across lactation permanent environment and the residual. Lactation curves of LP and LY increased from parturition to the peak of lactation and decreased thereafter, mirroring the typical curve of milk yield. Lactose percentage was greater in first- than later-parity cows. Heritabilities of LP and LY were 0.43±0.03 and 0.14±0.02, respectively, and LP and protein percentage were the most repeatable traits. Genetic correlations (r_a) of LP with somatic cell score, LY and milk freezing point were −0.22±0.08, 0.28±0.08 and −0.46±0.05, respectively. Genetic relationships of LY with milk yield (r_a=0.97±0.00), fat percentage (r_a=−0.71±0.06), protein percentage (r_a=−0.57±0.06) and protein yield (r_a=0.64±0.06) were moderate to strong. Results suggest that milk LP could be considered in breeding strategies to accelerate the gain of correlated low heritable traits. Further research is needed to evaluate the feasibility of including LP in the selection index of Italian Holstein population to address country-specific needs and market demands.     

Impact of bovine somatotropin on ranking for genetic value of dairy sires for milk yield traits and somatic cell score

Records of Holstein cows were used to examine how different models account for the effect of bovine somatotropin (bST) treatment on genetic evaluation of dairy sires for yield traits and somatic cell score. Data set 1 included 65,720 first-lactation records. Set 2 included 50,644 second-lactation records. Set 3 included 45,505 records for lactations three, four and five. Estimated breeding values (EBV) of sires were with three different animal models. With Model 1, bST administration was ignored. With Model 2, bST administration was used as a fixed effect. With Model 3, administration of bST was used to define the contemporary group (herd-year-month of calving-bST). Correlations for EBV of 1,366 sires with treated daughters between pairs of the three models were calculated for milk, fat and protein yields and somatic cell score for the three data sets. Correlations for EBV of sires between pairs of models for all traits ranged from 0.971 to 0.999. Fractions of sires with bST-treated progeny selected in common (top 10 to 15%) were 0.94 and usually greater for all pairs of models for all traits and parities. For this study, the method of statistical adjustment for bST treatment resulted in a negligible effect on genetic evaluations of sires when some daughters were treated with bST and suggests that selection of sires to produce the next generation of sires and cows might not be significantly affected by how the effect of bST is modeled for prediction of breeding values for milk, fat and protein yields and somatic cell score.     

Bayesian inference of genetic parameters for test-day milk yield, milk quality traits, and somatic cell score in Burlina cows

The aim of the study was to infer (co)variance components for daily milk yield, fat and protein contents, and somatic cell score (SCS) in Burlina cattle (a local breed in northeast Italy). Data consisted of 13576 monthly test-day records of 666 cows (parities 1 to 8) collected in 10 herds between 1999 and 2009. Repeatability animal models were implemented using Bayesian methods. Flat priors were assumed for systematic effects of herd test date, days in milk, and parity, as well as for permanent environmental, genetic, and residual effects. On average, Burlina cows produced 17.0 kg of milk per day, with 3.66 and 3.33% of fat and protein, respectively, and 358000 cells per mL of milk. Marginal posterior medians (highest posterior density of 95%) of heritability were 0.18 (0.09–0.28), 0.28 (0.21–0.36), 0.35 (0.25–0.49), and 0.05 (0.01–0.11) for milk yield, fat content, protein content, and SCS, respectively. Marginal posterior medians of genetic correlations between the traits were low and a 95% Bayesian confidence region included zero, with the exception of the genetic correlation between fat and protein contents. Despite the low number of animals in the population, results suggest that genetic variance for production and quality traits exists in Burlina cattle.     

Genetic correlations between yield traits or days open measured in cows and semen production traits measured in bulls

We used a bivariate animal model to investigate the genetic correlations between yield traits or days open (DO) as characters measured in cows and semen production traits as characters measured in bulls. Lactation records of 305-day milk, fat, and protein yields, and DO, from 386 809 first-lactation Holstein cows in Hokkaido, Japan, that calved between 2008 and 2014 were used. Semen production records were collected between 2005 and 2014 and included volume per ejaculate (VOL), sperm concentration (CON), number of sperm per ejaculate (NUM), progressive motility index of sperm (MOT), and MOT after freeze-thawing (A-MOT). Number of sperm per ejaculate was log-transformed into a NUM score (NUMS). A total of 30 373 semen production records from 1196 bulls were obtained. The pedigree file used for analysing the records was involving 885 345 animals. Heritability was estimated for VOL (0.42), CON (0.12), NUMS (0.37), MOT (0.08), and A-MOT (0.11). Weak and negative genetic correlations were recorded between yield traits measured in cows and VOL, CON or NUMS measured in bulls. Moderate and negative genetic correlations were obtained between DO and MOT (–0.42) or A-MOT (–0.43). Selection focused on MOT or A-MOT measured in bulls may therefore improve DO measured in cows.     

Genetics of milking characteristics in dairy cows

Genetic selection for milking speed is feasible. The existence of a correlation structure between milking speed and milk yield, however, necessitates a selection strategy to increase milking speed with no repercussion on genetic merit for milk yield. Residual milking duration (RMD) and residual milking duration including somatic cell score (RMDS), defined as the residuals from a regression model of milking duration on milk yield or milk yield plus somatic cell score (SCS) have been advocated. The objective of this study was to undertake a first ever genetic analysis of these novel traits. Data on electronically recorded milking duration and other milking characteristics from 235 005 test-day records on 74 608 cows in 1075 Irish dairy herds were available. Variance components for the milking characteristic traits were estimated using animal linear mixed models and covariances with other performance traits, including udder-related type traits, were estimated using sire models. The heritability of milking duration, RMD and RMDS was 0.20, 0.22 and 0.18, respectively. There were little differences in the heritability of RMD or RMDS when defined using genetic regression. The genetic standard deviation of RMDS defined on the phenotypic or genetic level was 36.8 s and 37.6 s, respectively, clearly indicating considerable exploitable genetic variation in milking duration independent of both milk yield and SCS. The genetic correlation between phenotypically derived RMDS and milk yield was favourable (−0.43), but RMDS was unfavourably genetically correlated with SCS (−0.30); the genetic correlations with both traits when RMDS was defined at a genetic level were zero. RMDS defined at the phenotypic level was negatively (i.e. unfavourable) genetically correlated (−0.35; s.e. = 0.15) with mastitis; however, when defined using genetic regression, shorter RMDS was not associated with greater expected incidence of mastitis. RMDS, defined at the genetic level, is a useful heritable trait with ample genetic variation for inclusion in a national breeding strategy without influencing genetic gain in either milk yield or udder health.     

Genetic and phenotypic parameters of lactations longer than 305 days (extended lactations)

Test-day milk yield and somatic cell count data over extended lactation (lactation to 540-600 days) were analysed considering part lactations as different traits and fitting random regression (RR) models. Data on Australian Jersey and Holstein Friesian (HF) were used to demonstrate the shape of the lactation curve and data on HF were used for genetic study. Test-day data from about 100 000 cows that calved between 1998 and 2005 were used for this study. In all analyses, a sire model was used.When part lactations were considered as different traits, protein yield early in the lactation (e.g. first 2 months) had a genetic correlation of about 0.8 with protein yield produced after 300 days of lactation. Genetic correlations between lactation stages that are adjacent to each other were high (0.9 or more) within parity. Across parities, genetic correlations were high for both protein and milk yield if they are within the same stage of lactation. Phenotypic correlations were lower than genetic correlations. Heritability of milk-yield traits estimated from the RR model varied from 0.15 at the beginning of the lactation to as high as 0.37 by the 4th month of lactation. All genetic correlations between different days in milk were positive, with the highest correlations between adjacent days in milk and decreasing correlations with increasing time-span. The pattern of genetic correlations between milk yield in the second 300 days (301 to 600 days of lactation) do not markedly differ from the pattern in the first 300 days of lactation. The lowest estimated genetic correlation was 0.15 between milk yield on days 45 and 525 of lactation. The result from this study shows that progeny of bulls with high estimated breeding values for yield traits and those that produce at a relatively high level in the first few months are the most likely candidates for use in herds favouring extended lactations.     

Relationships between conception rate in Holstein heifers and cows and milk yield at various stages of lactation

We investigated the relationships between conception rates (CRs) at first service in Japanese Holstein heifers (i.e. animals that had not yet had their first calf) and cows and their test-day (TD) milk yields. Data included records of artificial insemination (AI) for heifers and cows that had calved for the first time between 2000 and 2008 and their TD milk yields at 6 through 305 days in milk (DIM) from first through third lactations. CR was defined as a binary trait for which first AI was a failure or success. A threshold-linear animal model was applied to estimate genetic correlations between CRs of heifers or cows and TD milk yield at various lactation stages. Two-trait genetic analyses were performed for every combination of CR and TD milk yield by using the Bayesian method with Gibbs sampling. The posterior means of the heritabilities of CR were 0.031 for heifers, 0.034 for first-lactation cows and 0.028 for second-lactation cows. Heritabilities for TD milk yield increased from 0.324 to 0.433 with increasing DIM but decreased slightly after 210 DIM during first lactation. These heritabilities from the second and third lactations were higher during late stages of lactation than during early stages. Posterior means of the genetic correlations between heifer CR and all TD yields were positive (range, 0.082 to 0.287), but those between CR of cows and milk yields during first or second lactation were negative (range, −0.121 to −0.250). Therefore, during every stage of lactation, selection in the direction of increasing milk yield may reduce CR in cows. The genetic relationships between CR and lactation curve shape were quite weak, because the genetic correlations between CR and TD milk yield were constant during the lactation period.