First-lactation Holstein Cows Research Articles

Estimates of genetic parameters for rumination time, feed efficiency, and methane production traits in first-lactation Holstein cows

The large-scale recording of traits such as feed efficiency and methane emissions for use in genetic improvement programs is complex, costly, and time-consuming. Therefore, heritable traits that can be continuously recorded in dairy herds and are correlated to feed efficiency and methane emission traits could provide useful information for genetic evaluation. Rumination time has been suggested to be associated with feed efficiency, methane production (methane emission in g/day), and production traits at the phenotypic level. Therefore, the objective of this study was to investigate the genetic relationships among rumination time, feed efficiency, methane and production traits using 7,358 records from 656 first lactation Holstein cows. The estimated heritabilities were moderate for rumination time (0.45 ± 0.14), methane production (0.36 ± 0.12), milk yield (0.40 ± 0.08), fat yield (0.29 ± 0.06), protein yield (0.32 ± 0.07), and energy corrected milk (0.28 ± 0.07), while low and non-significant for feed efficiency (0.15 ± 0.07), which was defined as the residual of the multiple linear regression of DMI on ECM and MBW. A favorable negative genetic correlation was estimated between rumination time and methane production (-0.53 ± 0.24), while a positive favorable correlation was estimated between rumination time and energy corrected milk (0.49 ± 0.11). The estimated genetic correlation of rumination time with feed efficiency (-0.01 ± 0.17) was not significantly different from zero but showed a trend of a low correlation with dry matter intake (0.21 ± 0.13, P = 0.11). These results indicate that rumination time is genetically associated with methane production and milk production traits, but high standard errors indicate that further analyses should be conducted to verify these findings when more data for rumination time, methane production and feed efficiency become available.

Elucidating the factors and consequences of the severity of rumen acidosis in first-lactation Holstein cows during transition and early lactation.

First-lactation cows are particularly prone to subacute ruminal acidosis (SARA) during transition. Besides common risk factors of SARA, such as feeding of starch-rich diets, an individual severity of SARA in cows has been recently evidenced. Yet, the factors that play a role in SARA severity have not been elucidated. The main goal of this research was to evaluate the factors of SARA severity in first-lactation cows during transition and early lactation, which go beyond high-grain feeding, and to explore their impact on behavior, health, and fermentation in the rumen and hindgut. Twenty-four first-lactation Holstein cows with the same feeding regime were used starting from 3wk before the expected calving day until 10wk postpartum. Cows received a close-up diet (32% concentrate) until calving and were then transitioned to a lactation diet (60% concentrate) within 1 week. The SARA severity was assessed by cluster analysis of several rumen pH metrics, which revealed exceptionally longer and more severe SARA in cows denominated as high (n = 9), as compared to moderate (n = 9) and low (n = 6) SARA severity cows (P < 0.01). The logistic analysis showed that the length of close-up feeding, age at parturition, and the level of dry matter intake (DMI) were the main factors that influenced the cows' odds for high SARA severity (each P ≤ 0.01). Moreover, the ANOVA hinted differences in the metabolic activity of the ruminal microbiome to promote SARA severity, as indicated by highest ruminal propionate proportions (P = 0.05) in high SARA severity cows, also with similar DMI. The distinct SARA severity was marginally reflected in behavior and there were no effects of SARA severity or high-grain feeding on blood inflammation markers, which peaked at parturition regardless of SARA severity (P < 0.01). Still, ongoing high-grain feeding increased liver enzyme concentrations from 6wk postpartum on, compared to weeks before (P < 0.01), yet irrespectively of SARA severity. In conclusion, first-lactation cows differed in SARA severity under the same feeding regime, which was ascribed to management factors and differences in ruminal fermentation. Further research is warranted to validate these findings and to understand the mechanisms behind differences in the metabolic function of rumen microbiome, in particular in terms of evaluating markers for various SARA severity, as well as to evaluate potential long-term effects on health, performance, fertility, and longevity of dairy cows.

A Million-Cow Genome-Wide Association Study of Three Fertility Traits in U.S. Holstein Cows.

A genome-wide association study (GWAS) of the daughter pregnancy rate (DPR), cow conception rate (CCR), and heifer conception rate (HCR) using 1,001,374-1,194,736 first-lactation Holstein cows and 75,140-75,295 SNPs identified 7567, 3798, and 726 additive effects, as well as 22, 27, and 25 dominance effects for DPR, CCR, and HCR, respectively, with log10(1/p) > 8. Most of these effects were new effects, and some new effects were in or near genes known to affect reproduction including GNRHR, SHBG, and ESR1, and a gene cluster of pregnancy-associated glycoproteins. The confirmed effects included those in or near the SLC4A4-GC-NPFFR2 and AFF1 regions of Chr06 and the KALRN region of Chr01. Eleven SNPs in the CEBPG-PEPD-CHST8 region of Chr18, the AFF1-KLHL8 region of Chr06, and the CCDC14-KALRN region of Chr01 with sharply negative allelic effects and dominance values for the recessive homozygous genotypes were recommended for heifer culling. Two SNPs in and near the AGMO region of Chr04 that were sharply negative for HCR and age at first calving, but slightly positive for the yield traits could also be considered for heifer culling. The results from this study provided new evidence and understanding about the genetic variants and genome regions affecting the three fertility traits in U.S. Holstein cows.

Effects of High Concentrate-Induced Subacute Ruminal Acidosis Severity on Claw Health in First-Lactation Holstein Cows.

This study aimed to evaluate the effects of diet-induced subacute rumen acidosis (SARA) severity during transition and the early lactation period on claw health in 24 first-lactation Holstein heifers. All heifers were fed a 30% concentrate (in dry matter) close-up ration three weeks before calving, then switched to a high-concentrate ration (60% dry matter), which was fed until the 70th day in milk (DIM) to induce SARA. Thereafter, all cows were fed the same post-SARA ration with around 36% concentrate in dry matter. Hoof trimming was performed before calving (visit 1), at 70 (visit 2) and at 160 DIM (visit 3). All claw lesions were recorded, and a Cow Claw Score (CCS) was calculated for each cow. Locomotion scores (LCS 1-5) were assessed at two-week intervals. Intraruminal sensors for continuous pH measurements were used to determine SARA (pH below 5.8 for more than 330 min in 24 h). The cluster analysis grouped the cows retrospectively into light (≤11%; n = 9), moderate (>11-<30%; n = 7), and severe (>30%; n = 8) SARA groups, based on the percentage of days individual cows experienced SARA. Statistically significant differences were found between SARA groups light and severe in terms of lameness incidence (p = 0.023), but not for LCS and claw lesion prevalence. Further, the analysis of maximum likelihood estimates revealed that for each day experiencing SARA, the likelihood of becoming lame increased by 2.52% (p = 0.0257). A significant increase in white line lesion prevalence was observed between visits 2 and 3 in the severe SARA group. The mean CCS in severe SARA group cows were higher at each visit compared to cows in the other two groups, but without statistical significance. Overall, this is the first study indicating that first-lactation cows fed a similar high-concentrate diet but with a higher severity of SARA tended to have poorer claw health, albeit with only partial statistical evidence.

503 Late-Breaking: Using Random Regression Models to Estimate Genetic Parameters for Milk Production Traits under Different Levels of Heat Stress in Canadian Holstein Cattle

Abstract Genetically selecting animals that are heat tolerant can help to maintain production efficiency of the dairy industry under thermal challenges. During the summer months in Canada, dairy cows could be under thermal stress. For instance, in Ontario, cows are expected to be 105 days under thermal stress. The objective of this study was to estimate genetic parameters for milk production traits under heat stress in Canadian Holstein cows. A total of 604,050 test-day records for milk, fat, and protein yields from 78,516 first-lactation Holstein cows were used in this study. The production data was combined with daily maximum temperature-humidity index (THI; ranging from 30 to 83) calculated based on meteorological records from public weather stations. In total, 58 weather stations were located within a maximum distance of 20 km from each of the 1,284 herds included in this study. Single-trait random regressions on THI, fitting linear splines to describe the random additive genetic and permanent environmental effects, were used in the genetic analyses. Based on preliminary analyses, the knots of the spline functions were set to THI = 70, 57, and 58 for milk, fat, and protein yields, respectively. Therefore, the models assumed a linear effect of heat stress beyond the THI threshold, allowing to estimate the additive genetic effect under thermal comfort and at different levels of heat stress. The genetic correlations estimated between thermal comfort and heat stress were -0.09, -0.33, and -0.21 for milk, fat, and protein yields, respectively. This indicates an antagonistic relationship between level of production and heat tolerance. However, small differences in heritability estimates were observed above the THI thresholds (they ranged from 0.20 to 0.22 for milk, 0.18 to 0.15 for fat, and 0.14 to 0.15 for protein yields). This study demonstrates the possibility of genetically selecting for more heat-tolerant animals in Canada.

Single-step genomic best linear unbiased predictor genetic parameter estimations and genome-wide associations for milk fatty acid profiles, interval from calving to first insemination, and ketosis in Holstein cattle

Milk fatty acids (FA) have been suggested as biomarkers for early-lactation metabolic diseases and for female fertility status. The aim of the present study was to infer associations between FA, the metabolic disorder ketosis (KET), and the interval from calving to first insemination (ICF) genetically and genomically. In this regard, we focused on a single-step genomic BLUP approach, allowing consideration of genotyped and ungenotyped cows simultaneously. The phenotypic data set considered 38,375 first-lactation Holstein cows, kept in 45 large-scale co-operator herds from 2 federal states in Germany. The calving years for these cows were from 2014 to 2017. Concentrations in milk from the first official milk recording test-day for saturated, unsaturated (UFA), monounsaturated (MUFA), polyunsaturated, palmitic, and stearic (C18:0) FA were determined via Fourier-transform infrared spectroscopy. Ketosis was defined as a binary trait according to a veterinarian diagnosis key, considering diagnoses within a 6-wk interval after calving. A subset of 9,786 cows was genotyped for 40,989 SNP markers. Variance components and heritabilities for all Gaussian distributed FA and for ICF, and for binary KET were estimated by applying single-step genomic BLUP single-trait linear and threshold models, respectively. Genetic correlations were estimated in series of bivariate runs. Genomic breeding values for the single-step genomic BLUP estimations were dependent traits in single-step GWAS. Heritabilities for FA were moderate in the range from 0.09 to 0.20 (standard error = 0.02-0.03), but quite small for ICF (0.08, standard error = 0.01) and for KET (0.05 on the underlying liability scale, posterior standard deviation = 0.02). Genetic correlations between KET and UFA, MUFA, and C18:0 were large (0.74 to 0.85, posterior standard deviation = 0.14-0.19), and low positive between KET and ICF (0.17, posterior standard deviation = 0.22). Genetic correlations between UFA, MUFA, and C18:0 with ICF ranged from 0.34 to 0.46 (standard error = 0.12). In single-step GWAS, we identified a large proportion of overlapping genomic regions for the different FA, especially for UFA and MUFA, and for saturated and palmitic FA. One identical significantly associated SNP was identified for C18:0 and KET on BTA 15. However, there was no genomic segment simultaneously significantly affecting all trait categories ICF, FA, and KET. Nevertheless, some of the annotated potential candidate genes DGKA, IGFBP4, and CXCL8 play a role in lipid metabolism and fertility mechanisms, and influence production diseases in early lactation. Genetic and genomic associations indicate that Fourier-transform infrared spectroscopy FA concentrations in milk from the first official test-day are valuable predictors for KET and for ICF.

A Large-Scale Genome-Wide Association Study of Epistasis Effects of Production Traits and Daughter Pregnancy Rate in U.S. Holstein Cattle.

Epistasis is widely considered important, but epistasis studies lag those of SNP effects. Genome-wide association study (GWAS) using 76,109 SNPs and 294,079 first-lactation Holstein cows was conducted for testing pairwise epistasis effects of five production traits and three fertility traits: milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FPC), protein percentage (PPC), and daughter pregnancy rate (DPR). Among the top 50,000 pairwise epistasis effects of each trait, the five production traits had large chromosome regions with intra-chromosome epistasis. The percentage of inter-chromosome epistasis effects was 1.9% for FPC, 1.6% for PPC, 10.6% for MY, 29.9% for FY, 39.3% for PY, and 84.2% for DPR. Of the 50,000 epistasis effects, the number of significant effects with log10(1/p) ≥ 12 was 50,000 for FPC and PPC, and 10,508, 4763, 4637 and 1 for MY, FY, PY and DPR, respectively, and A × A effects were the most frequent epistasis effects for all traits. Majority of the inter-chromosome epistasis effects of FPC across all chromosomes involved a Chr14 region containing DGAT1, indicating a potential regulatory role of this Chr14 region affecting all chromosomes for FPC. The epistasis results provided new understanding about the genetic mechanism underlying quantitative traits in Holstein cattle.

The hidden cost of disease: II. Impact of the first incidence of lameness on production and economic indicators of primiparous dairy cows

Lameness is a persistent and underreported health and welfare problem in the dairy industry, resulting in reduced cow performance and profitability as well as early culling. The study objectives were (1) to quantify the impact of the first instance of lameness, at different stages of lactation, on production and economic performance, and (2) to further quantify the impacts of the first instance of lameness when only cows that remain in the herd for at least 100 d in milk (DIM) and those that remain for 305 DIM are included in the analysis. A retrospective longitudinal study was conducted using pre-existing data from animal health records and Dairy Herd Improvement Association records. Data were edited based on selected inclusion criteria, yielding a data set containing records from 15,159 first-lactation Holstein cows from 120 herds with year of first calving between 2003 and 2014. Lame cows were assigned to 1 of 4 groups based on when in the lactation the first event of lameness occurred: transition (1-21 DIM), early lactation (22-100 DIM), mid-lactation (101-200 DIM), or late lactation (201+ DIM). Mid- and late-lactation lame cows were also stratified by cumulative milk yield before the lameness event. Healthy cows (i.e., no recorded lameness event) were randomly assigned for each lactation stage, with mid-lactation healthy and late-lactation healthy cows similarly stratified. Production performance (cumulative milk, fat, and protein yield) and economic performance [milk value, margin over feed cost (MOFC), and gross profit] were analyzed using a mixed model with herd as a random effect. Cumulative milk yields were 811 to 1,290 kg lower for lame cows than for healthy cows, with milk component yields undergoing similar reductions. Reductions in milk yield contributed to losses in milk value (-Can$527 to -Can$1,083; -US$419 to -US$862) and MOFC (-Can$510 to -Can$774; -US$406 to -US$616). Higher losses were reported using gross profit (-Can$753 to -Can$1,052; -US$599 to -US$837), which includes all lameness-related costs. Production and performance losses were smaller when 100 DIM and 305 DIM thresholds were applied (i.e., exclusion of cows culled before 100 and 305 DIM, respectively), however, mid- and late-lactation lame cows maintained high levels of significant losses for all 6 variables analyzed. Lameness also led to higher levels of culling, masking losses for transition and early-lactation lame cows in the 305-DIM analysis. Increasing producer understanding of the costs associated with lameness not only serves to provide insight to producers for more informed culling decisions, but may also help producers weigh the costs of adopting new methods and technologies targeted at reducing on-farm lameness.

The hidden cost of disease: I. Impact of the first incidence of mastitis on production and economic indicators of primiparous dairy cows

Mastitis is a highly prevalent disease, which negatively affects cow performance, profitability, welfare, and longevity. The objectives of this study were (1) to quantify the impact of the first instance of mastitis, at different stages of lactation, on production and economic performance, and (2) to further quantify the impact of the first instance of mastitis when only cows that remain in the herd for at least 100 d in milk (DIM) and those that remain for 305 DIM are included in the analysis. A retrospective longitudinal study was conducted using data from existing animal health record files and Dairy Herd Improvement records. After editing based on selected inclusion criteria and completeness of health records, data consisted of records from first-lactation Holstein cows, from 120 herds, that calved for the first time between 2003 and 2014, inclusive. Mastitic cows were assigned to 1 of 4 groups based on when in the lactation the first event of mastitis occurred: transition (1-21 DIM), early lactation (22-100 DIM), mid lactation (101-200 DIM), or late lactation (201+ DIM). Mid-lactation and late-lactation mastitic cows were also stratified by cumulative milk yield before the mastitis event. Healthy cows (i.e., no recorded mastitis event) were randomly assigned for each lactation stage, with mid-lactation healthy and late-lactation healthy cows similarly stratified. Production performance (cumulative milk, fat, and protein yield) and economic performance [milk value, margin over feed cost (MOFC), and gross profit] were analyzed using a mixed model with herd as a random effect. Significant losses in cumulative milk yield (-382 to -989 kg) and correspondingly lower fat and protein yields were found in mastitic cows, with transition and late-lactation mastitic cows having the highest losses. Drops in production translated to significant reductions in cumulative milk value (-Can$287 to -Can$591; -US$228 to -US$470), MOFC (-Can$243 to -Can$540; -US$193 to -US$429), and gross profit (-Can$649 to -Can$908; -US$516 to -US$722) for mastitic cows at all stages. Differences between mastitic and healthy cows in the early lactation and transition stages remained for all variables in the 100-DIM analysis, but, aside from gross profit, were nonsignificant in the 305-DIM analysis. Gross profit accounted for all costs associated with mastitis and thus continued to be lower for mastitic cows at all stages, even in the 305-DIM analysis in which culled cows were omitted (-Can$485 to -Can$979; -US$386 to -US$779). The research reflects the performance implications of mastitis, providing more information upon which the producer can make informed culling decisions and maximize both herd profitability and cow longevity.

9 Pregnancy loss in Holstein lactating recipient cows diagnosed pregnant by pregnancy-associated glycoprotein test in blood

The objective of this study was to evaluate the pregnancy loss (PL) between Day 30 (P30) and Day 80 (P80) of pregnancy in lactating Holstein recipients that received an invivo- (flush) or invitro-produced (IVF) embryo. The recipient cows were located at Maddox Dairy in Riverdale, CA, USA, a Holstein herd that milks 3500 cows with a 305-day mature-equivalent milk production of 12 800 kg. First-lactation cows were enrolled in a Presynch-Ovsynch oestrus-synchronization program and scheduled to be artificially bred for the first time at 80 days after calving or to receive an embryo 7 or 8 days after the expected heat (recipients). The data from 590 pregnancies (1045 embryo transfers) from embryos transferred between January 2018 and March 2019 was analysed. Only grade 1 embryos (from morula to hatched blastocyst stage) produced invivo or invitro from Holstein donor heifers, lactating or dry cows, were transferred fresh (invivo or invitro) or frozen-thawed (invivo), and pregnancy rates are shown in Table 1. The ruminant trophoblast produces pregnancy-associated glycoproteins (PAG) that can be detected by enzyme-linked immunoassay (ELISA) in the blood of pregnant cattle as early as 28 days after insemination. Various dairy herds in the USA have been using this test to supplement or replace the use of transrectal ultrasonography for early pregnancy diagnosis. Blood was sampled on P30 after expected heat day (23 days after embryo transfer) from the recipient cows and sent to IDEXX for the PAG Bovine Pregnancy Test, which was reconfirmed on P80 of pregnancy by transrectal ultrasonography. Pregnancy loss was considered to have occurred when a cow was pregnant on P30 but not pregnant on P80. The variable PL was analysed by binary logistic regression in the MINITAB program, and the model included effects of donor status (heifer vs. milk vs. dry) and embryo type. The total PL was 12.2% (72/590) and the details are shown in Table 1. No effect of donor status (P=0.80) was detected. However, there was effect of embryo type (P=0.004). The IVF embryos had a PL of 18.0% compared with 9.5% for the invivo-produced embryos. Further research should be performed to study heifer embryos PL, because currently more embryos are being produced from very young donor cows and sires due to intensive use of genomic testing. In conclusion, there is a higher PL in lactating dairy recipients receiving IVF fresh embryos compared with fresh or frozen invivo-produced embryos. Table 1.Embryo transfer (ET) pregnancy rates (PR, P30) in first-lactation Holstein cows (top) and pregnancy loss (PL) from Day 30 (P30) to 80 (P80) in first-lactation recipient Holstein cows (bottom) Item Heifer donor Lactating donor Dry donor Total ET P30 PR% ET P30 PR% ET P30 PR% ET P30 PR% Invivo - fresh 6 3 50.0 43 25 58.1 329 206 62.6 378 234 61.9 Invivo - frozen 75 38 50.7 221 129 58.4 296 167 56.4 IVF - fresh 123 61 49.6 91 42 46.2 157 86 54.8 371 189 50.9 Total 129 64 49.6 209 105 50.2 707 421 59.5 1045 590 56.5 P30 P80 PL% P30 P80 PL% P30 P80 PL% P30 P80 PL% Invivo - fresh 3 3 0.0 25 24 4.0 206 186 9.7 234 213 9.0 Invivo - frozen 38 35 7.9 129 115 10.9 167 150 10.2 IVF - fresh 61 48 21.3 42 35 16.7 86 72 16.3 189 155 18.0 Total 64 51 20.3 105 94 10.5 421 373 11.4 590 518 12.2

A Large-Scale Genome-Wide Association Study in U.S. Holstein Cattle.

Genome-wide association study (GWAS) is a powerful approach to identify genomic regions and genetic variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. We conducted a large-scale GWAS using 294,079 first-lactation Holstein cows and identified new additive and dominance effects on five production traits, three fertility traits, and somatic cell score. Four chromosomes had the most significant SNP effects on the five production traits, a Chr14 region containing DGAT1 mostly had positive effects on fat yield and negative effects on milk and protein yields, the 88.07–89.60 Mb region of Chr06 with SLC4A4, GC, NPFFR2, and ADAMTS3 for milk and protein yields, the 30.03–36.67 Mb region of Chr20 with C6 and GHR for milk yield, and the 88.19–88.88 Mb region with ABCC9 as well as the 91.13–94.62 Mb region of Chr05 with PLEKHA5, MGST1, SLC15A5, and EPS8 for fat yield. For fertility traits, the SNP in GC of Chr06, and the SNPs in the 65.02–69.43 Mb region of Chr01 with COX17, ILDR1, and KALRN had the most significant effects for daughter pregnancy rate and cow conception rate, whereas SNPs in AFF1 of Chr06, the 47.54–52.79 Mb region of Chr07, TSPAN4 of Chr29, and NPAS1 of Chr18 had the most significant effects for heifer conception rate. For somatic cell score, GC of Chr06 and PRLR of Chr20 had the most significant effects. A small number of dominance effects were detected for the production traits with far lower statistical significance than the additive effects and for fertility traits with similar statistical significance as the additive effects. Analysis of allelic effects revealed the presence of uni-allelic, asymmetric, and symmetric SNP effects and found the previously reported DGAT1 antagonism was an extreme antagonistic pleiotropy between fat yield and milk and protein yields among all SNPs in this study.

Phenotypic relationships, genetic parameters, genome-wide associations, and identification of potential candidate genes for ketosis and fat-to-protein ratio in German Holstein cows

Energy demand for milk production in early lactation exceeds energy intake, especially in high-yielding Holstein cows. Energy deficiency causes increasing susceptibility to metabolic disorders. In addition to several blood parameters, the fat-to-protein ratio (FPR) is suggested as an indicator for ketosis, because a FPR >1.5 refers to high lipolysis. The aim of this study was to analyze phenotypic, quantitative genetic, and genomic associations between FPR and ketosis. In this regard, 8,912 first-lactation Holstein cows were phenotyped for ketosis according to a veterinarian diagnosis key. Ketosis was diagnosed if the cow showed an abnormal carbohydrate metabolism with increased content of ketone bodies in the blood or urine. At least one entry for ketosis in the first 6 wk after calving implied a score = 1 (diseased); otherwise, a score = 0 (healthy) was assigned. The FPR from the first test-day was defined as a Gaussian distributed trait (FPRgauss), and also as a binary response trait (FPRbin), considering a threshold of FPR = 1.5. After imputation and quality controls, 45,613 SNP markers from the 8,912 genotyped cows were used for genomic studies. Phenotypically, an increasing ketosis incidence was associated with significantly higher FPR, and vice versa. Hence, from a practical trait recording perspective, first test-day FPR is suggested as an indicator for ketosis. The ketosis heritability was slightly larger when modeling the pedigree-based relationship matrix (pedigree-based: 0.17; SNP-based: 0.11). For FPRbin, heritabilities were larger when modeling the genomic relationship matrix (pedigree-based: 0.09; SNP-based: 0.15). For FPRgauss, heritabilities were almost identical for both pedigree and genomic relationship matrices (pedigree-based: 0.14; SNP-based: 0.15). Genetic correlations between ketosis with FPRbin and FPRgauss using either pedigree- or genomic-based relationship matrices were in a moderate range from 0.39 to 0.71. Applying genome-wide association studies, we identified the specific SNP rs109896020 (BTA 5, position: 115,456,438 bp) significantly contributing to ketosis. The identified potential candidate gene PARVB in close chromosomal distance is associated with nonalcoholic fatty liver disease in humans. The most important SNP contributing to FPRbin was located within the DGAT1 gene. Different SNP significantly contributed to ketosis and FPRbin, indicating different mechanisms for both traits genomically.

Effects of partial replacement of soybean meal with other protein sources in diets of lactating cows

The protein nutrition of dairy cows is of great importance because of its direct influence on milk production, reproductive efficiency, and feeding cost. Eight first-lactation Holstein cows were randomly assigned to two contemporary 4 × 4 Latin squares in a 2 × 2 factorial design to evaluate the effects of replacing soybean meal with yeast-derived microbial protein (YMP) as a protein source (0% or 1.5% of dry matter (DM)) and its combination with slow-release urea (SRU; 0% or 0.75% of DM) on DM intake and milk production and composition, as well as blood parameters and nitrogen balance. Each experimental period lasted 28 days, with 21 days of adaptation and 7 days of data collection. The diets were formulated to attend the nutritional recommendations of the National Research Council and consisted of 49% forage (47% corn silage and 2% Tifton hay) and 51% concentrate, with 16.8% CP and 1.6 Mcal net energy for lactation/kg DM. For diets without YMP, the inclusion of SRU decreased DM intake, milk production as well as N intake and balance, but did not affect efficiency of production, milk composition or most of blood parameters. On the contrary, for diets with YMP, DM intake and milk production were increased by inclusion of SRU, while minor effects were observed for milk efficiency and composition, blood parameters as well as N intake, excretion and balance. When diets with SRU were compared, the inclusion of YMP increased DM intake, 4% fat-corrected milk, and N intake and balance (P<0.05), with no differences in milk production (kg/day), milk energy, efficiency of milk production or most of the blood parameters. For diets without SRU, YMP inclusion decreased DM intake, milk production, milk energy, N intake, fecal N and N balance (P<0.05), with no effects on milk efficiency and composition, or most of blood parameters. In conclusion, the use of YMP, SRU or both as partial substitutes of soybean meal in the diet of lactating cows has no negative effects on productivity parameters.

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.

Heritabilities and genetic correlations in the same traits across different strata of herds created according to continuous genomic, genetic, and phenotypic descriptors

The most common approach in dairy cattle to prove genotype by environment interactions is a multiple-trait model application, and considering the same traits in different environments as different traits. We enhanced such concepts by defining continuous phenotypic, genetic, and genomic herd descriptors, and applying random regression sire models. Traits of interest were test-day traits for milk yield, fat percentage, protein percentage, and somatic cell score, considering 267,393 records from 32,707 first-lactation Holstein cows. Cows were born in the years 2010 to 2013, and kept in 52 large-scale herds from 2 federal states of north-east Germany. The average number of genotyped cows per herd (45,613 single nucleotide polymorphism markers per cow) was 133.5 (range: 45 to 415 genotyped cows). Genomic herd descriptors were (1) the level of linkage disequilibrium (r2) within specific chromosome segments, and (2) the average allele frequency for single nucleotide polymorphisms in close distance to a functional mutation. Genetic herd descriptors were the (1) intra-herd inbreeding coefficient, and (2) the percentage of daughters from foreign sires. Phenotypic herd descriptors were (1) herd size, and (2) the herd mean for nonreturn rate. Most correlations among herd descriptors were close to 0, indicating independence of genomic, genetic, and phenotypic characteristics. Heritabilities for milk yield increased with increasing intra-herd linkage disequilibrium, inbreeding, and herd size. Genetic correlations in same traits between adjacent levels of herd descriptors were close to 1, but declined for descriptor levels in greater distance. Genetic correlation declines were more obvious for somatic cell score, compared with test-day traits with larger heritabilities (fat percentage and protein percentage). Also, for milk yield, alterations of herd descriptor levels had an obvious effect on heritabilities and genetic correlations. By trend, multiple trait model results (based on created discrete herd classes) confirmed the random regression estimates. Identified alterations of breeding values in dependency of herd descriptors suggest utilization of specific sires for specific herd structures, offering new possibilities to improve sire selection strategies. Regarding genomic selection designs and genetic gain transfer into commercial herds, cow herds for the utilization in cow training sets should reflect the genomic, genetic, and phenotypic pattern of the broad population.

Genetic analysis of groups of mid-infrared predicted fatty acids in milk

The objective of this study was to investigate genetic variability of mid-infrared predicted fatty acid groups in Canadian Holstein cattle. Genetic parameters were estimated for 5 groups of fatty acids: short-chain (4 to 10 carbons), medium-chain (11 to 16 carbons), long-chain (17 to 22 carbons), saturated, and unsaturated fatty acids. The data set included 49,127 test-day records from 10,029 first-lactation Holstein cows in 810 herds. The random regression animal test-day model included days in milk, herd-test date, and age-season of calving (polynomial regression) as fixed effects, herd-year of calving, animal additive genetic effect, and permanent environment effects as random polynomial regressions, and random residual effect. Legendre polynomials of the third degree were selected for the fixed regression for age-season of calving effect and Legendre polynomials of the fourth degree were selected for the random regression for animal additive genetic, permanent environment, and herd-year effect. The average daily heritability over the lactation for the medium-chain fatty acid group (0.32) was higher than for the short-chain (0.24) and long-chain (0.23) fatty acid groups. The average daily heritability for the saturated fatty acid group (0.33) was greater than for the unsaturated fatty acid group (0.21). Estimated average daily genetic correlations were positive among all fatty acid groups and ranged from moderate to high (0.63-0.96). The genetic correlations illustrated similarities and differences in their origin and the makeup of the groupings based on chain length and saturation. These results provide evidence for the existence of genetic variation in mid-infrared predicted fatty acid groups, and the possibility of improving milk fatty acid profile through genetic selection in Canadian dairy cattle.

Short communication: Association of foot and leg conformation and body weight with claw disorders in Spanish Holstein cows.

Data of first-lactation Holstein cows including claw disorders, foot and leg type traits, and structure and capacity type traits were used to study the phenotypic association of conformation and body weight with claw health status. The claw disorders studied were digital and interdigital dermatitis, sole ulcer, white line disease, and an overall claw disorder trait, indicating the presence of any of the 6 claw disorders recorded by the Spanish I-SAP program. Body weight was estimated indirectly with the Von Bertalanffy equation using live weight, which was also estimated from body depth, stature, and chest width. Cows with poor scores for foot and leg traits were more likely to have a claw disorder, with the exception of rear leg rear view and foot and leg composite that did not show any association with the studied disorders. The heavier the cow was, the higher was the probability of having sole ulcer, white line disease, or overall claw disorder trait, but digital and interdigital dermatitis, as an infectious disorder, did not show any association with body weight. Therefore, it is recommended that the preventive trimming routine be improved, which means checking each cow at least once per lactation and trimming if necessary, to achieve balanced weight-bearing for heavier cows and cows with poor feet and leg classification.

Production and physiological responses of heat-stressed lactating dairy cattle to conductive cooling

The objective of this research was to test the effectiveness of conductive cooling in alleviating heat stress of lactating dairy cows. A conductive cooling system was built with waterbeds (Dual Chamber Cow Waterbeds, Advanced Comfort Technology Inc., Reedsburg, WI) modified to circulate chilled water. The experiment lasted 7 wk. Eight first-lactation Holstein cows producing 34.4±3.7kg/d of milk at 166±28 d in milk were used in the study. Milk yield, dry matter intake (DMI), and rectal temperature were recorded twice daily, and respiration rate was recorded 5 times per day. During wk 1, the cows were not exposed to experimental heat stress or conductive cooling. For the remaining 6 wk, the cows were exposed to heat stress from 0900 to 1700h each day. During these 6 wk, 4 of the 8 cows were cooled with conductive cooling (experimental cows), and the other 4 were not cooled (control cows). The study consisted of 2 thermal environment exposures (temperature-humidity index mean ± standard deviation of 80.7±0.9 and 79.0±1.0) and 2 cooling water temperatures (circulating water through the water mattresses at temperatures of 4.5°C and 10°C). Thus, a total of 4 conductive cooling treatments were tested, with each treatment lasting 1 wk. During wk 6, the experimental and control cows were switched and the temperature-humidity index of 79.0±1.0 with 4.5°C cooling water treatment was repeated. During wk 7, waterbeds were placed directly on concrete stalls without actively cooling the water. Least squares means and P-values for the different treatments were calculated with multivariate mixed models. Conductively cooling the cows with 4.5°C water decreased rectal temperature by 1.0°C, decreased respiration rate by 18 breaths/min, increased milk yield by 5%, and increased DMI by 14% compared with the controls. When the results from the 2 cooling water temperatures (4.5°C and 10°C circulating water) were compared, we found that the rectal temperature from 4.5°C cooling water was 0.3°C lower than the rectal temperature with 10°C cooling water, but the other measurements (respiration rate, milk production, and DMI) did not show a statistically significant difference between the cooling water temperatures. Placing waterbeds on concrete stalls without additional cooling did not have a measurable effect in alleviating the heat stress of the cows.

Genetic analysis of heat stress effects on yield traits, udder health, and fatty acids of Walloon Holstein cows

Genetic parameters that considered tolerance for heat stress were estimated for production, udder health, and milk composition traits. Data included 202,733 test-day records for milk, fat, and protein yields, fat and protein percentages, somatic cell score (SCS), 10 individual milk fatty acids (FA) predicted by mid-infrared spectrometry, and 7 FA groups. Data were from 34,468 first-lactation Holstein cows in 862 herds in the Walloon region of Belgium and were collected between 2007 and 2010. Test-day records were merged with daily temperature-humidity index (THI) values based on meteorological records from public weather stations. The maximum distance between each farm and its corresponding weather station was 21km. Linear reaction norm models were used to estimate the intercept and slope responses of 23 traits to increasing THI values. Most yield and FA traits had phenotypic and genetic declines as THI increased, whereas SCS, C18:0, C18:1 cis-9, and 4 FA groups (unsaturated FA, monounsaturated FA, polyunsaturated FA, and long-chain FA) increased with THI. Moreover, the latter traits had the largest slope-to-intercept genetic variance ratios, which indicate that they are more affected by heat stress at high THI levels. Estimates of genetic correlations within trait between cold and hot environments were generally high (>0.80). However, lower estimates (<=0.67) were found for SCS, fat yield, and C18:1 cis-9, indicating that animals with the highest genetic merit for those traits in cold environments do not necessarily have the highest genetic merit for the same traits in hot environments. Among all traits, C18:1 cis-9 was the most sensitive to heat stress. As this trait is known to reflect body reserve mobilization, using its variations under hot conditions could be a very affordable milk biomarker of heat stress for dairy cattle expressing the equilibrium between intake and mobilization under warm conditions.

Variance and covariance components and genetic parameters for fat and protein yield of first-lactation Holstein cows using random regression models

Summary Background: the genetic parameters of the lactation curve in dairy cattle can be analyzed as longitudinal data using Random Regression Models (RRM). Objective: to estimate the (co) variance components and genetic parameters for fat (F) and protein (P) yield in first lactation Holstein cows of Antioquia (Colombia) by RRM based on Legendre polynomials. Methods: monthly F and P records (9,479) from 1,210 first-lactation Holstein cows were used. Twenty-two and 24 RRM were used for F and P, respectively, with different orthogonal Legendre-polynomial orders to estimate the fixed-curve population coefficients and predict direct-genetic additive and permanent environment effects. The models considered homogeneous and heterogeneous residual variances of 5, 7, and 10 classes. Results: the best fit for F was the fourth order model for the population fixed-curve and the additive genetic effect, and the third order for the permanent environment with seven heterogeneous variances. The best fit for P was the fifth order model for the population fixed-curve and the additive genetic and permanent environmental effects with five heterogeneous variances. The variance for the animals' genetic, phenotypic, permanent environment, and residual effects for both F and P decreased as lactation progressed. F and P heritabilities were between 0.13 and 0.38, and 0.12 and 0.32, respectively. Conclusion: first-birth animals can be selected in Antioquia for F and P characteristics. Selection should be done preferably at the beginning of lactation since they reach the highest heritability values at this time. Keywords: genetic evaluation, heritability, milk quality, solids in milk. Resumen Antecedentes: los parámetros genéticos de la curva de lactancia en ganado de leche pueden ser analizados como datos longitudinales usando Modelos de Regresión Aleatoria (RRM). Objetivo: estimar mediante RRM basados en polinomios de Legendre componentes de (co) varianza y parámetros genéticos para producción de grasa (F) y proteína (P) láctea en vacas Holstein de primera lactancia de Antioquia (Colombia). Métodos: se incluyeron 9.479 registros mensuales de F y P pertenecientes a 1.210 vacas Holstein de primera lactancia. Para F y P se usaron 22 y 24 RRM respectivamente, con diferentes órdenes de polinomio ortogonal de Legendre para estimar los coeficientes de la curva fija de la población, la predicción de los efectos genético aditivo directo y del ambiente permanente. Los modelos consideraron varianzas residuales homogéneas y heterogéneas de 5, 7 y 10 clases. Resultados: para F, el mejor modelo fue el de cuarto orden para la curva fija de la población y el efecto genético aditivo, y de tercer orden para el ambiente permanente con siete varianzas heterogéneas. Para P, el modelo que presentó mejor ajuste fue el de quinto orden para la curva fija de la población, el efecto genético aditivo y el ambiente permanente y 5 varianzas heterogéneas. Para ambas características las varianzas genética aditiva directa, fenotípica, de ambiente permanente y residual disminuyeron a medida que avanzaba la lactancia. Las heredabilidades para F y P estuvieron entre 0,13 y 0,38, y entre 0,12 y 0,32, respectivamente. Conclusión: es posible realizar la selección para F y P en animales de primer parto en el departamento de Antioquia, preferiblemente al inicio de la lactancia, ya que ambas características presentan heredabilidades altas en esta etapa. Palabras clave: calidad de leche, evaluación genética, heredabilidad, sólidos lácteos. Resumo Antecedentes: os parâmetros genéticos da produção de leite podem ser estimados usando Modelos de Regressão Aleatória (RRM). Objetivo: estimar por RRM com base em polinômios de Legendre os componentes de variância e covariância e os parâmetros genéticos para produção de gordura (F) e proteína (P) em vacas leiteiras de primeira lactação da raça holandesa em rebanhos de Antioquia, Colômbia. Métodos: foram avaliadas 9.479 registros mensais de F e P pertencentes a 1.210 vacas. Para F e P foram usados 22 e 24 RRM, respectivamente, com diferentes ordens de polinomiais ortogonais de Legendre para estimar os coeficientes da curva fixa da população, os efeitos genéticos aditivo direito, do ambiente permanente e residual. Os modelos consideraram variâncias residuais homogêneas e heterogêneas de 5, 7 e 10 classes. Resultados: para F, o melhor modelo foi o de quarta ordem para a curva fixa da população e o efeito genético aditivo, e de terceira ordem para o ambiente permanente com sete variâncias heterogêneas. Para P, o modelo que forneceu o melhor ajuste foi do quinto ordem para a curva fixa da população, o efeito genético aditivo e de ambiente permanente, e com 5 variâncias heterogêneas. Para ambas as características as variâncias genética aditiva direita, fenotípica e residual diminuíram no tempo. As herdabilidades para F e P ficaram entre 0,13 e 0,38 e entre 0,12 e 0,32, respectivamente. Conclusão: é possível fazer a seleção de animais da raça holandesa para F e P, de preferência no início da lactação, pois as duas características têm altas herdabilidades. Palavras chave: avaliação genética, herdabilidade, qualidade do leite, sólidos lácteos.