Heritabilities and genetic and phenotypic correlations for milk production and fertility traits of spring-calved once-daily or twice-daily milking cows in New Zealand.

The objectives of this study were to estimate the genetic and phenotypic correlations and heritabilities for milk production and fertility traits in spring-calved once-daily (OAD) milking cows for the whole season in New Zealand and compare those estimates with twice-daily (TAD) milking cows. Data used in the study consisted of 69,252 first parity cows from the calving seasons 2015-2016 to 2017-2018 in 113 OAD and 531 TAD milking herds. Heritability estimates for production and fertility traits were obtained through single-trait animal models, and estimates of genetic and phenotypic correlations were obtained through bivariate animal models. Heritability estimates of production traits varied from 0.26 to 0.61 in OAD and from 0.13 to 0.63 in TAD. Heritability estimates for fertility traits were low in both OAD and TAD milking cow populations, and estimates were consistent (OAD: 0.01 to 0.10 and TAD: 0.01 to 0.08) across milking regimens. Estimates of phenotypic and genetic correlations among production traits were consistent across populations. In both populations, phenotypic correlations between milk production and fertility traits were close to zero, and most of the genetic correlations were antagonistic. In OAD milking cows, genetic correlations of milk and lactose yields with the start of mating to conception, 6-wk in-calf, not-in-calf, and 6-wk calving rate were close to zero. Interval from first service to conception was negatively genetically correlated with milk and lactose yields in OAD milking cows. Protein percentage was positively genetically correlated with 3-wk and 6-wk submission, 3-wk in-calf, 6-wk in-calf, first service to conception, 3-wk calving, and 6-wk calving rate in the TAD milking cow population, but these correlations were low in the OAD milking cow population. Further studies are needed to understand the relationship of protein percentage and fertility traits in the OAD milking system. The phenotypic correlations between fertility traits were similar in OAD and TAD milking populations. Genetic correlations between fertility traits were strong (≥0.70) in cows milked TAD, but genetic correlations varied from weak to strong in cows milked OAD. Further research is required to evaluate the interaction between genotype by milking regimen for fertility traits in terms of sire selection in the OAD milking cow population.

Estimation of genetic parameters and individual and maternal breed, heterosis, and recombination loss effects for production and fertility traits of spring-calved cows milked once daily or twice daily in New Zealand.

The objectives of this study were to estimate genetic parameters and individual and maternal breed, heterosis, and recombination loss effects for milk production and fertility traits of Holstein Friesian (F), Jersey (J), and crossbred Holstein Friesian and Jersey (F × J) cows milked once daily (OAD) or twice daily (TAD) in New Zealand. Data on 278,776 lactations from 30,217 OAD and 170,680 TAD milking cows across 644 spring-calving herds were available. Genetic parameters and individual and maternal breed, heterosis, and recombination loss estimates were obtained from univariate animal models. Heritability and repeatability estimates for milk production, milk composition, and fertility traits were consistent across the milking frequencies. Heritability estimates for yields of milk, fat, protein, and lactose varied between 0.21 and 0.29 in OAD and TAD. Heritability estimates for fertility traits ranged from 0.01 to 0.08 in both populations, and estimates were slightly greater in TAD than OAD milking cows. In both milking populations, individual breed effects for yields were in favor of F cows; however, maternal breed effects for yields were in favor of J dams. Jersey cows were more fertile than the F cows in both milking populations, but maternal breed effects for fertility traits were in favor of F dams. Individual heterosis effects were favorable for all traits and were consistent across milking regimens. Crossbred F × J cows had significantly shorter intervals from start of mating to first service and from start of mating to conception, and a higher proportion of 3-wk submission, 3-wk in calf, and 3-wk calving relative to the average of purebred F and J cows. Recombination loss effects were not always unfavorable for production and fertility traits, but most estimates were small with larger standard errors. Favorable maternal heterosis effects were associated with production traits in both milking systems, but maternal heterosis effects were less likely to influence reproductive performance.

Reducing milking frequency in early lactation improved the energy status but reduced milk yield during the whole lactation of primiparous Holstein cows consuming a total mixed ration and pasture.

To investigate the immediate and long-term performance effects of milking frequency during early lactation of primiparous dairy cows consuming a total mixed ration and pasture, 20 Holstein cows were assigned in a randomized block design to either once-daily (1×) or twice-daily (2×) milking during the first 8 wk of lactation (treatment period). After the treatment period, all cows were milked 2× until wk 43 of lactation. Cows were fed a total mixed ration (approximately 15 kg of DM/cow per day) and allowed to graze an oat pasture (Avena sativa). Dry matter intake was 19.1 kg of DM/cow per day on average and was not affected by treatments. Milk yield was 40% lower in cows milked 1× during the treatment period, and a carryover effect existed until wk 21 of lactation, resulting in a final reduction of 15% of milk yield in the whole lactation. Milk lactose concentration decreased, whereas fat and protein concentrations increased for cows milked 1×. Mobilization of energy reserves during the treatment period occurred in both groups, but cows milked 1× showed greater body condition score and greater backfat thickness. In conclusion, milking 1× during the first 8 wk of lactation resulted in immediate and carryover negative effects on milk and milk solid yield without affecting feed intake, resulting in the improved energy status of primiparous dairy cows.

Genetic parameters of milk production traits in response to a short once-daily milking period in crossbred Holstein × Normande dairy cows.

Despite its potential utility for predicting cows' milk yield responses to once-daily milking (ODM), the genetic basis of cow milk trait responses to ODM has been scarcely if ever described in the literature, especially for short ODM periods. This study set out to (1) estimate the genetic determinism of milk yield and composition during a 3-wk ODM period, (2) estimate the genetic determinism of milk yield responses (i.e., milk yield loss upon switching cows to ODM and milk yield recovery upon switching them back to twice-daily milking; TDM), and (3) seek predictors of milk yield responses to ODM, in particular using the first day of ODM. Our trial used 430 crossbred Holstein × Normande cows and comprised 3 successive periods: 1 wk of TDM (control), 3 wk of ODM, and 2 wk of TDM. Implementing ODM for 3 wk reduced milk yield by 27.5% on average, and after resuming TDM cows recovered on average 57% of the milk lost. Heritability estimates in the TDM control period and 3-wk ODM period were, respectively, 0.41 and 0.35 for milk yield, 0.66 and 0.61 for milk fat content, 0.60 and 0.80 for milk protein content, 0.66 and 0.36 for milk lactose content, and 0.20 and 0.15 for milk somatic cell score content. Milk yield and composition during 3-wk ODM and TDM periods were genetically close (within-trait genetic correlations between experimental periods all exceeding 0.80) but were genetically closer within the same milking frequency. Heritabilities of milk yield loss observed upon switching cows to ODM (0.39 and 0.34 for milk yield loss in kg/d and %, respectively) were moderate and similar to milk yield heritabilities. Milk yield recovery (kg/d) upon resuming TDM was a trait of high heritability (0.63). Because they are easy to measure, TDM milk yield and composition and milk yield responses on the first day of ODM were investigated as predictors of milk yield responses to a 3-wk ODM to easily detect animals that are well adapted to ODM. Twice-daily milking milk yield and composition were found to be partly genetically correlated with milk yield responses but not closely enough for practical application. With genetic correlations of 0.98 and 0.96 with 3-wk ODM milk yield losses (in kg/d and %, respectively), milk yield losses on the first day of ODM proved to be more accurate in predicting milk yield responses on longer term ODM than TDM milk yield.

Once-daily milking during late lactation in pasture-fed dairy cows has minor effects on feed intake, body condition score gain, and hepatic gene expression.

Milking cows once daily (1×) is a management practice occasionally used during mid/late lactation in pasture-based systems. It has been postulated that 1× milking will reduce dry matter intake (DMI) and increase body condition score (BCS) gain; however, this has not been quantified. Lactating, pregnant Holstein-Friesian dairy cows (n=52) were allocated to either 1× or twice-daily (2×) milking in mid-January (summer, 175d in milk). To obtain accurate DMI measurements, cows underwent 4 periods in a Calan gate indoor feeding facility, interspersed with grazing outdoors. Milk production, body weight (BW), and BCS were recorded 2 wk before treatment start (-2 wk) and weekly thereafter. Blood variables were recorded at -2 wk and weekly when indoors. Liver was biopsied at -2, 2, and 10 wk, and hepatic gene expression measured using quantitative PCR. Milking cows 1× tended to lower DMI (17.8 vs. 18.2 kg of dry matter), but increased BCS gain (0.36 vs. 0.13 BCS units) and BW (546 vs. 533 kg) at wk 12 relative to 2×. The greater BCS and BW of cows milked 1× compared with 2× were reflected in lower plasma concentrations of nonesterified fatty acids and lower transcription of genes involved in the oxidation of fatty acids, indicating reduced release and processing of fatty acids. Cows milked 1× produced 20% less milk, and although milk fat and protein concentrations were increased relative to cows milked 2×, yields of fat and protein were 14 and 17% less, respectively. The reduction in milk production with 1× milking (14.1 vs. 16.8 kg/cow per d energy-corrected milk) was accompanied by increases in blood concentrations of glucose and insulin, with a concurrent decrease in the transcription of the insulin receptor and gluconeogenic genes. These results indicate a coordinated response to reduce glucose production due to decreased mammary demand. Expression of 2 genes linked to inflammation and adipokine signaling was reduced in cows milked 1× and may indicate a lower inflammatory state in the liver of cows milked 1× in late lactation. No effect was found of milking frequency during late lactation on milk production in the subsequent lactation. In summary, although 1× milking tended to reduce DMI and increase BCS in late lactation, these effects were lower than what is commonly supposed in pasture-based dairy systems. The modest BCS gains need to be considered with the reduced milk production when adopting 1× milking as a management strategy.

Effects of DGAT1 K232A polymorphism and milking frequency on milk composition and spontaneous lipolysis in dairy cows.

Milk spontaneous lipolysis (SL) of milk triglycerides is induced by the lipoprotein lipase, a milk native enzyme, and leads to an accumulation of partial glycerides and free fatty acids that are responsible for the deterioration of the taste of milk products. The gene coding for diacylglycerol acyltransferase 1 (DGAT1), an enzyme implicated in triglycerides synthesis, has an important polymorphic site at the K232A locus. This gene is well known to modulate milk composition. No data are available on the effects of DGAT1 on SL. Thus, a trial was carried out to evaluate the effects of DGAT1 K232A polymorphism on milk SL upon milking frequency variations [once- (ODM) and twice-daily milking (TDM)]. Twenty-one cows were divided into 3 groups according their DGAT1 K232A genotype: 8 cows had the KK genotype of DGAT1 (KK cows), 8 had the KA genotype (KA cows), and 5 had the AA genotype (AA cows). The trial consisted in 3 successive periods: 3 wk of TDM, 3 of ODM, and 4 of TDM. Samples were collected for fat and protein contents, SL, fatty acid, and protein profiles determinations. The KK cows presented higher fat and protein contents, lower milk production, and higher κ-casein percentage. No significant difference in fatty acid composition was noted between groups. The SL was twice as high for KK cows in TDM situations (1.13 vs. 0.59 and 0.63mEq/100g of fat, respectively, for KK, KA, and AA cows during the first period of TDM, and 0.46 vs. 0.25 and 0.21mEq/100g of fat, respectively, for KK, KA, and AA during the second period of TDM). The SL remained lower in TDM2 than in TDM1. During ODM, no difference in SL was found between groups and SL remained below 0.2mEq/100g of fat. These results demonstrate the existence of a correlation between DGAT1 genotypes and spontaneous lipolysis, in interaction with an environmental factor, milking frequency, although it has not been possible to clarify the causal mechanism at this stage.

Lactation traits associated with short- and long-term once-daily milking performance in New Zealand crossbred dairy cattle.

The main objectives of this study were to establish the relative value of milk yields under twice-daily milking (TDM) as a predictor of yield and yield loss under once-daily milking (ODM), and to understand the role of residual milk and udder storage capacity-related traits in regulating yield and yield loss during ODM. A Holstein-Friesian × Jersey crossbred herd was established over 2 seasons (years), as 2 individual cohorts on the same farm, managed on a pasture-based system over 4 lactations. Short-term (1-wk) ODM studies, with a starting total of 690 cows, were undertaken in mid- and late-lactation in lactation 2 and in mid-lactation in lactation 3 for each cohort. A 10-wk study of ODM performance began in mid-lactation in lactation 3, whereas lactation 4 was a full-lactation assessment of ODM. In the short-term studies, milk yield under ODM was well predicted (R(2)=0.7 to 0.8 in 5 of 6 studies) by the daily yield under TDM in the week before ODM. Yield loss (kg/d) increased with increasing milk yield and with increasing somatic cell count (SCC), although predictions were relatively poor (R(2)=0.09 to 0.30). Yield loss (%) decreased with increasing TDM yield in 3 of the 6 studies and was positively correlated with SCC during ODM. Nevertheless, ODM yield loss, in absolute or percentage terms, was a poorly repeatable trait in grazing cows. Part of the variation in yield loss percentage (30%) was positively associated with residual milk (%), measured pretrial, during measurement of functional udder capacity in lactation 3. Total production (kg of milk) over the full-lactation ODM study in lactation 4 was correlated with total production in the 10-wk trial in lactation 3 (r=0.72 and 0.63 for cohorts 1 and 2, respectively). Identifying the highest- and lowest-producing 10% of animals during the full lactation of ODM indicated that poor production was associated with high yields of residual milk (measured in lactation 3) and, conversely, high production was associated with low yields of residual milk, relative to the other 80% of animals. These same "high" and "low" production groups from lactation 4 had similar differences in performance in the earlier short-term studies and a larger or smaller percentage yield loss associated with the residual milk measurement. Breeding strategies for ODM may benefit, therefore, from greater emphasis on selecting for a low residual milk fraction to optimize milking performance. Nevertheless, the level of milk production under TDM is a strong phenotypic predictor of milk production under ODM.

Reducing milking frequency during nutrient restriction has no effect on the hepatic transcriptome of lactating dairy cattle.

The objective of this study was to investigate if a reduced milking frequency altered the effect of dietary energy restriction on the hepatic transcriptome of grazing dairy cows during early lactation. Multiparous Holstein-Friesian and Holstein-Friesian × Jersey cows (n = 120) were milked twice daily (2×) from calving until 34 ± 6 days in milk (mean ± SD). Cows were then allocated to one of four treatments in a 2 × 2 factorial arrangement. Treatments consisted of two milking frequencies [2× or once daily (1×)] and two feeding levels for 3 wk: adequately fed (AF) or underfed (UF, 60% of AF). Liver tissue was biopsied from 12 cows per treatment after 3 wk of treatment, and the hepatic transcriptome was profiled with an Agilent 4 × 44k bovine microarray. Over 2,900 genes were differentially expressed in response to the energy restriction; however, no effects resulted from changes to milking frequency. This may indicate that after 3 wk of 1× milking, any changes to the liver transcriptome that may have occurred earlier have returned to normal. After 3 wk of energy restriction, gene expression patterns indicate that glucose-sparing pathways were activated, and gluconeogenesis was increased in UF cows. Genes involved in hepatic stress were upregulated in response to the energy restriction indicative of the pressure energy restriction places on liver function. Other pathways upregulated included "cytoskeletal remodeling," indicating that a 3 wk energy restriction resulted in molecular changes to assist tissue remodeling. Overall, 1× milking does not modify the hepatic transcriptome changes that occur in response to an energy restriction.

Once-daily milking during a feed deficit decreases milk production but improves energy status in early lactating grazing dairy cows.

The objective of this study was to investigate the effect of milking frequency (MF) at 2 feeding levels (FL) on milk production, body condition score, and metabolic indicators of energy status in grazing dairy cows during early lactation. Multiparous Holstein-Friesian and Holstein-Friesian × Jersey cows (n=120) grazed pasture and were milked twice daily (2×) from calving until 34 ± 6 d in milk (mean ± standard deviation). Cows were then allocated to 1 of 4 treatments in a 2 × 2 factorial arrangement. Treatments consisted of 2 FL: adequately fed [AF; 14.3 kg dry matter intake (DMI)/cow per d] or underfed (UF; 8.3 kg of DMI/cow per d) and 2 MF: 2× or once daily (1×). Treatments were imposed for 3 wk. After the treatment period, all cows were offered a generous pasture allowance (grazing residuals >1,600 kg of dry matter/ha) and milked 2×. During the 3-wk treatment period, we observed an interaction between FL and MF for energy-corrected milk (ECM), such that the decrease due to 1× milking was greater in AF than in UF cows (20 and 14% decrease, respectively). No interactions were found posttreatment. Cows previously UF produced 7% less ECM than AF cows during wk 4 to 12; however, no subsequent effect was observed of the previous underfeeding. Cows previously milked 1× produced 5% less ECM during wk 4 to 12, and differences remained during wk 13 to 23. During the 3-wk treatment period, UF cows lost 0.2 body condition score units (1-10 scale) and this was not affected by 1× milking. During the treatment period, UF cows had lower plasma glucose, insulin, and insulin-like growth factor I, and greater nonesterified fatty acids and ß-hydroxybutyrate concentrations than AF cows. Cows milked 1× had greater plasma glucose, insulin, and insulin-like growth factor I, and lower nonesterified fatty acids and ß-hydroxybutyrate concentrations compared with cows milked 2×. In conclusion, energy status was improved by 1× milking; however, when UF cows were milked 1×, milk production was reduced by more than underfeeding alone. The immediate and residual responses to 1× milking need to be considered when using this management strategy during a feed deficit.

Response profiles of dairy cows to a single 24-h milking interval in relation with milk proteolysis, udder expansion and immune traits.

An extended milking interval of 24 h (24-h milking interval (24h-MI)) constitutes the acute phase of cow adaptation to once-daily milking (ODM). A recent trial including 724 24h-MI challenges demonstrated that milk yield responses to this acute phase of ODM are highly variable (from+22% to -52% of milk yield when switching to the 24h-MI, mean=-25.3%) and that factors such as stage of lactation parity and milk yield level influenced cows' responses but did not account for all individual variability. Additional traits related to physiological, immune and behavioural adaptation were measured on a subset (96 observations) of this data set. This study aimed to determine (1) the relationship of these traits with cows' milk yield responses, (2) their ability - combined with previously identified traits - to help predict milk yield responses to 24h-MI (adaptive profiles). The 24h-MI challenge consisted of three successive periods: one control week of twice-daily milking (cTDM), one single day of 24h-MI and then 13 days of TDM (pTDM). Milk yield responses to the 24h-MI (corrected for effects of stage of lactation, parity, milk yield level and milk yield) were related to physiological traits measured during cTDM (milk flow rate, presence or absence of interleukin-8) and to their changes during the 24h-MI (absolute increase in milk flow rate and relative udder distension). Analysis of associations between milk yield responses, stage of lactation, parity, milk yield level, proteolysis, udder expansion and immune traits found three adaptive cow profile clusters. Cows in cluster 1 had a less compliant udder than cows in cluster 2, and they lost more milk during the 24h-MI than cluster-2 and cluster-3 cows. After resuming twice daily-milking (TDM), cluster-2 cows fully recovered the milk they had lost during the 24h-MI. On the opposite, cluster-3 cows did not recover the milk they lost, likely due to udder inflammation during cTDM, as suggested by elevated concentrations of interleukin-8 in their milk. These results combining new traits with stage of lactation, parity and milk yield level constitute a first step towards predicting individual cow responses to a 24h-MI.