Reasons for disposal and cull cow value of Holstein cows compared with Holstein, Jersey, Montbeliarde, Normande, and Viking Red crossbred cows.

Cull cows from dairy farms in the United States account for 10% of the beef market; however, few studies have evaluated the effects of crossbreeding on reasons for disposal and cull value of dairy cows. The objective of the study was to compare reasons for disposal and cull cow value of Holstein and crossbred cows (n = 1,292) in an experimental dairy herd at the University of Minnesota West Central Research and Outreach Center, Morris, Minnesota. Cows were Holstein (n = 272), 1964 genetic control Holstein (n = 161), 3-breed crossbred cows (n = 538) composed of the Montbéliarde, Viking Red, and Holstein (MVH) breeds, and 3-breed crossbred cows (n = 321) composed of the Normande, Jersey, and Viking Red (NJV) breeds. The 1964 genetic control Holsteins cows originated from a design initiated at the University of Minnesota for comparison of the1964 Holstein and contemporary Holstein cows selected for production across time. Records spanned from January 2010 to December 2023. Reasons for disposal were recorded in PCDart Herd Management Software, and cull value and body weight were from receipts from livestock cull markets. For all cows, the primary reason for disposal was for reproduction (44.4%), mastitis (18%), other reasons (14.6%), died (10.3%), dairy purposes (8.4%), and low production (4.3%). Independent variables for statistical analysis of cull value were the fixed effects of body weight at time of culling, DIM at culling (0 to 49 DIM, 50 to 99 DIM, 100 to 149 DIM, 150 to 199 DIM, 200 to 249 DIM, 250 to 299 DIM, 300 to 350 DIM, and 350+ DIM), year (2010 to 2023), season (spring, summer, autumn, winter), parity (1, 2, 3, 4, 5+), primary reason for culling, breed group and sire breed group nested within breed group. Least squares means for bodyweight at culling was 558 kg for Holsteins, 543 kg for 1964 Holsteins, 551 kg for MVH crossbred cow, and 499 kg for NJV crossbred cows. For the analysis across time, 2014 had the highest cull value ($1,126.98) and 2020 had the lowest gross cull value ($515.21). Cows culled greater than 300 DIM had the highest cull value and cows culled less than 50 DIM had the lowest cull value. Cows culled during the spring and summer had higher cull value compared with cows culled during the autumn and winter. The Holstein cows ($730.04) had lower gross value compared with 1964 Holsteins ($804.38) cows and MVH ($767.39) cows. The NJV ($771.39) cows were not different from crossbreds sired by Montbeliarde, Viking Red, and Holstein bulls. Results for comparisons of breed group are from one experimental herd, so inferences to the wider dairy cow population should be undertaken with caution. In summary, dairy producers may receive greater cull value from crossbred cows compared with Holstein cows.

Lactation curves of Montbéliarde-sired and Viking Red-sired crossbred cows and their Holstein herdmates in commercial dairies.

Lactation curves were estimated for Montbéliarde (MO) × Holstein (HO) and Viking Red (VR) × HO 2-breed crossbred cows and for MO × VR/HO and VR × MO/HO 3-breed crossbred cows and their HO herdmates from test-day observations in 7 high-performance herds that participated in a designed study. Cows calved from 2010 to 2017. Test-day observations from milk recording were used to fit the lactation curves of cows in their first 3 lactations. Lactations of cows were required to have at least 250 DIM and to have at least 6 test days ≤265 DIM. Lactation curves from random regression (RR) were compared for 305-d production (kg), peak production (kg), peak day of production, and production from 4 to 103 DIM (kg), from 104 to 205 DIM (kg), and from 206 to 305 DIM (kg) for milk, fat, and protein. Also, the persistency of production was compared. First-lactation versus second- and third-lactation cows were analyzed separately for both the 2-breed and 3-breed crossbred cows and their respective HO herdmates. Legendre polynomial RR had the best goodness of fit for the lactation curves compared with Ali-Schaeffer and Wilmink RR from the test-day observations of milk, fat, and protein production. For fluid milk production of first-lactation cows, the MO × HO 2-breed crossbreds were not different from their HO herdmates for any of the lactation-curve characteristics, except persistency. However, the VR × HO 2-breed crossbreds had less fluid milk production compared with their HO herdmates. For first lactation, the MO × HO 2-breed crossbreds had more persistency of milk, fat, and protein production compared with their HO herdmates. The first-lactation MO × VR/HO 3-breed crossbreds had more persistency of fluid milk production compared with their HO herdmates. For second and third lactations, both the MO × HO and the VR × HO 2-breed crossbreds had higher fat production compared with their HO herdmates. Furthermore, the MO × HO 2-breed crossbreds had more protein production (kg) in all 3 periods of lactation compared with their HO herdmates. Crossbred cows may have advantages over HO cows for persistency of production in high-performance herds.

Genetic parameters and association of national evaluations with breeding values for health traits in US organic Holstein cows.

Among other regulations, organic cows in the United States cannot receive antibiotics and preserve their organic status, emphasizing the importance of prevention of illness and benefit of high genetic merit for disease resistance. At the same time, data underlying national genetic evaluations primarily come from conventional cows, drawing concern to the possibility of a genotype by environment interaction whereby the value of a genotype varies depending on the environment, and potentially limits the relevance of these evaluations to organic cows. The objectives of this study were to characterize the genetics of and determine the presence of genotype by environment interaction for health traits in US organic dairy cows. Individual cow health data were obtained from 16 US Department of Agriculture certified organic dairy farms from across the United States that used artificial insemination and maintained detailed records. Data were obtained for the following traits: died, lameness, mastitis, metabolic diseases (displaced abomasum, ketosis, and milk fever), reproductive diseases (abortion, metritis, and retained placenta), transition health events (any health event occurring 21 d before or after parturition), and all health events. Binary phenotypes (1 = diseased, 0 = otherwise) for 38,949 lactations on 19,139 Holstein cows were used. Genotypes from 2,347 cows with 87.5% or greater Holstein breed-based representation were incorporated into single-step multitrait threshold animal models that included stayability (1 = completed lactation, 0 = otherwise). Gibbs sampling was used. Genomic predicted transmitting abilities (gPTA) from national genetic evaluations were obtained for sires for production, fitness, health, and conformation traits. We approximated genetic correlations for sires using these gPTA and our estimated breeding values. We also regressed health phenotypes on cow estimated breeding values and sire gPTA. Heritabilities (± standard error) ranged from 0.03 ± 0.01 (reproductive diseases) to 0.11 ± 0.03 (metabolic diseases). Most genetic correlations among health traits were positive, though the genetic correlation between metabolic disease and mastitis was -0.42 ± 0.17. Approximate genetic correlations between disease resistance for our health trait categories and disease resistance for the nationally-evaluated health traits generally carried the expected sign with the strongest correlation for mastitis (0.72 ± 0.084). Regression coefficients carried the expected sign and were mostly different from zero, indicating that evaluations from primarily conventional herd data predicted health on organic farms. In conclusion, use of national evaluations for health traits should afford genetic improvement for health in US organic herds.

Birth traits of Holstein calves compared with Holstein, Jersey, Montbéliarde, Normande, and Viking Red-sired crossbred calves.

Holstein (HO) calves, 3-breed crossbred calves of Montbéliarde, Viking Red, and HO (MVH), and 3-breed crossbred calves of Normande, Jersey, and Viking Red (NJV) were compared for gestation length (GL), calf weight at birth (CW), calving difficulty (CD), and stillbirth (SB) in 2 research herds at the University of Minnesota. Calves were born from January 2009 to December 2019. For the St. Paul and Morris herds, HO calves (n = 1,121) were compared with MVH calves (n = 1,393) from primiparous and multiparous cows. For the single herd analysis at Morris, HO calves (n = 476), MVH calves (n = 922), and NJV calves (n = 405) were compared from primiparous and multiparous cows. Primiparous and multiparous births were analyzed separately because multiparous cows had multiple births, and CD and SB are likely different traits for primiparous and multiparous cows. Statistical analysis of GL, CW, CD, and SB included fixed effects of sex of calf, herd, breed group of calf, and year-season of calving. For the St. Paul and Morris herds, HO calves from primiparous (278 d) and multiparous (279 d) HO cows had shorter GL compared with MVH calves from primiparous (280 d) and multiparous (282 d) crossbred cows. The HO calves (39.4 and 43.2 kg, respectively) from primiparous and multiparous HO cows had lower CW compared with MVH calves (40.3 and 44.3 kg, respectively) from primiparous and multiparous crossbred cows. Calving difficulty and SB were not different for HO and MVH calves from primiparous and multiparous cows. For the single herd analysis at Morris, HO calves (278 and 279 d, respectively) from primiparous and multiparous HO cows had shorter GL compared with MVH calves (281 and 282 d, respectively) and NJV calves (282 and 282 d, respectively) from primiparous and multiparous crossbred cows. The CW of HO calves (38.6 and 42.0 kg, respectively) from primiparous and multiparous HO cows was lower compared with MVH calves (39.7 and 42.9 kg, respectively), but higher compared with NJV calves (35.1 and 38.0 kg, respectively) from primiparous and multiparous crossbred cows. Calving difficulty and SB did not differ for HO, MVH, and NJV calves from primiparous and multiparous cows. The longer GL for crossbred calves and higher CW for MVH calves did not increase CD and SB for primiparous and multiparous cows. Dairy producers may implement 3-breed rotational crossbreeding systems that include the HO, Jersey, Normande, Montbéliarde, and Viking Red breeds, and some breeds may increase GL and CW without an increase in CD and SB.

Comparison of 3-breed rotational crossbreds of Montbéliarde, Viking Red, and Holstein with Holstein cows fed 2 alternative diets for dry matter intake, production, and residual feed intake.

The objective of this study was to compare 3-breed rotational crossbred (CB) cows of the Montbéliarde, Viking Red, and Holstein (HO) breeds with HO cows fed 2 alternative diets for dry matter intake (DMI), fat plus protein production (CFP), body weight (BW), body condition score (BCS), feed efficiency, and residual feed intake (RFI) from 46 to 150 days in milk (DIM) during first lactation. The CB cows (n = 17) and HO cows (n = 19) calved from September 2019 to March 2020. Cows were fed either a traditional total mixed ration diet (TRAD) or a higher fiber, lower starch total mixed ration diet (HFLS). The HFLS had 21% more corn silage, 47% more alfalfa hay, 44% less corn grain, and 43% less corn gluten feed than the TRAD. The 2 diets were analyzed for dry matter content, crude protein, forage digestibility, starch, and net energy for lactation. The BW and BCS were recorded once weekly. Daily milk, fat, and protein production were estimated from twice monthly milk recording with random regression. Measures of efficiency were CFP per kilogram of DMI and DMI per kilogram of BW. The RFI from 46 to 150 DIM was the residual error from regression of DMI on milk energy, metabolic BW, and the energy required for change in BW. Statistical analysis of all variables included the fixed effects of diet, breed group, and the interaction of diet and breed group. The CB cows fed HFLS had less DMI (-12%) and lower DMI/BW (-14%) compared with the HO cows fed TRAD. For CFP, CB and HO cows were not different when fed TRAD or HFLS. Furthermore, the CB cows fed HFLS had higher BW (+50 kg) compared with HO cows fed HFLS. The CB cows fed TRAD had higher BCS than HO cows fed TRAD and HO cows fed HFLS (+0.46 and +0.62, respectively). The HO cows fed TRAD had more DMI (+14%) and lower CFP per kilogram of DMI (-12%) compared with the HO cows fed HFLS. In addition, mean RFI from 46 to 150 DIM was lower and more desirable for CB cows fed HFLS (-120.0 kg) compared with HO cows fed TRAD (85.3 kg). Dairy producers may feed either TRAD or HFLS to CB cows without loss of CFP.

Genetic parameters of calf morbidity and stayability for US organic Holstein calves.

The objectives of this study were to estimate genetic parameters of calf health in organic US Holstein calves. Calves were born on farms across the United States from 2006 to 2019. Three calf health traits were evaluated in the study: calf respiratory disease until 365 d of age, calf scours until 60 d of age, and heifer stayability until 365 d of age. For respiratory disease and scours, animals were assigned a phenotype of 0 if they were healthy and a phenotype of 1 if they were diseased. For stayability, animals were assigned a phenotype of 0 if they were removed from the herd by 365 d of age and 1 if they remained in the herd at 365 d of age. Genetic parameters were estimated from threshold models that included the fixed effects of mean, year-season of birth, and dam age (respiratory disease and scours only) as well as the random effects of herd-year of birth and additive genetics. Heritability estimates were 0.100, 0.075, and 0.085 for respiratory disease, scours, and stayability, respectively. Solutions for estimated breeding values for respiratory disease and scours were transformed from disease risk to disease resistance by reversing the signs before calculating genetic correlations such that higher values of scours, respiratory disease, and stayability were favored. There was a moderate favorable genetic correlation estimate between respiratory disease resistance and stayability of 0.675. However, genetic correlation estimates between respiratory disease resistance and scours resistance (0.148) and between scours resistance and stayability (0.165) were low. Estimated breeding value correlations between calf health traits and other traits evaluated nationally were generally low in magnitude. The strongest correlation estimates were with longevity, particularly between stayability and heifer livability (0.217) and between stayability and cow livability (0.288); respiratory disease resistance was also favorably correlated with heifer (0.190) and cow (0.178) livability. Correlations with cow health traits were generally low and unfavorable. Linear models including the random effect of herd-by-sire indicated that herd-by-sire accounted for approximately 2% of phenotypic variance for scours and stayability, which may indicate a genotype by environment interaction effect for these traits. In conclusion, there is significant genetic variation in organic calf health, and there was evidence of genotype by environment interaction.

Genetic parameters of passive transfer of immunity for US organic Holstein calves.

Passive transfer of immunity is important for calf health and survival. The objectives of this study were to estimate genetic parameters for calf passive transfer of immunity through producer-recorded serum total protein (STP) and to determine associations with other routinely evaluated traits in organic Holstein calves (n = 16,725) that were born between July 2013 to June 2018; a restricted subset (n = 7,518) of calves with known Holstein maternal grandsires was analyzed separately. Producers measured STP on farm, and STP was extracted from farm management software. Failure of passive transfer of immunity (FPT) was declared for calves with STP ≤5.2 g/dL. Calves that had the opportunity to reach 1 yr of age were recorded as either staying in the herd or leaving the herd (STAY365). Univariate and threshold models were fitted for STP and FPT, respectively, and included the fixed effects of herd-year-month of birth, calf age in days at STP measurement, dam age in years, and random effects of animal and birthdate within herd. Model effects for STAY365 included the fixed effects of herd-year-month of birth and random effects of animal and birthdate within herd. Multivariate analyses of STP with FPT or STAY365 were conducted to determine the genetic correlation between traits and STP was also regressed on gestation length. Heritability estimates of STP were 0.06 and 0.08 for full and restricted data, respectively. Heritability estimates for FPT were 0.04 and 0.06 for full and restricted data, respectively. The genetic correlation between STP and FPT was near unity. Heritability estimates for STAY365 ranged from 0.08 to 0.11 with genetic correlation estimates between STP and STAY365 ranging from 0.19 and 0.25. Approximate genetic correlations were estimated for sires (n = 302 and n = 256 for full and restricted data, respectively) with at least 10 daughters for STP and predicted transmitting abilities for health, calving traits, and production. Positive approximate genetic correlations were estimated for STP with cow livability, productive life, net merit dollars, and milk yield; favorable approximate genetic correlations were observed for daughter and sire calving ease, and sire stillbirth. Longer gestation length was associated with reduced STP genetically and phenotypically. These results suggest that passive transfer as measured through STP is heritable and favorably correlated with current measures of health, calving, and production.

Evaluation of solar photovoltaic systems to shade cows in a pasture-based dairy herd.

The combined use of solar photovoltaics and agriculture may provide farmers with an alternative source of income and reduce heat stress in dairy cows. The objective of this study was to determine the effects on grazing cattle under shade from a solar photovoltaic system. The study was conducted at the University of Minnesota West Central Research and Outreach Center in Morris, Minnesota on a grazing dairy. Twenty-four crossbred cows were randomly assigned to 2 treatment groups (shade or no shade) from June to September in 2019. The replicated (n = 4) treatment groups of 6 cows each were provided shade from a 30-kW photovoltaic system. Two groups of cows had access to shade in paddocks, and 2 groups of cows had no shade in paddocks. All cows were located in the same pasture during summer. Behavior observations and milk production were evaluated for cows during 4 periods of summer. Boluses and an eartag sensor monitored internal body temperature, activity, and rumination on all cows, respectively. Independent variables were the fixed effects of breed, treatment group, coat color, period, and parity, and random effects were replicate group, date, and cow. No differences in fly prevalence, milk production, fat and protein production, or drinking bouts were observed between the treatment groups. Shade cows had more ear flicks (11.4 ear flicks/30 s) than no-shade cows (8.6 ear flicks/30 s) and had dirtier bellies and lower legs (2.2 and 3.2, respectively) than no-shade cows (1.9 and 2.9, respectively). During afternoon hours, shade cows had lower respiration rates (66.4 breaths/min) than no-shade cows (78.3 breaths/min). From 1200 to 1800 h and 1800 to 0000 h, shade cows had lower body temperature (39.0 and 39.2°C, respectively) than no-shade cows (39.3 and 39.4°C, respectively). Furthermore, between milking times (0800 and 1600 h), the shade cows had lower body temperature (38.9°C) than no-shade cows (39.1°C). Agrivoltaics incorporated into pasture dairy systems may reduce the intensity of heats stress in dairy cows and increase well-being of cows and the efficiency of land use.

Evaluation of the RumiWatch system as a benchmark to monitor feeding and locomotion behaviors of grazing dairy cows.

Direct visual observation is a common method for validation of animal behavior technologies; however, visual observations are time consuming and subject to human error. The objective of this study was to evaluate the RumiWatch system (Itin and Hoch GmbH, Liestal, Switzerland), which is composed of a noseband sensor and a pedometer, for monitoring feeding and locomotion behaviors of grazing dairy cows, to determine its accuracy for use as a benchmark in validation studies. The study was conducted at the University of Minnesota West Central Research and Outreach Center in Morris, Minnesota, from May to June 2018. Two experiments were conducted and validated: (1) feeding and locomotion behaviors and (2) rumination cycle and grazing bites. Lactating crossbred dairy cows (n = 12) were offered pasture for 22 h/d, and cows were milked twice daily. Visual observations were recorded by 3 observers with the Pocket Observer app (Noldus Information Technology, Leesburg, VA). The first experiment determined agreement for visual observations and the RumiWatch noseband sensor and pedometer from 144 h of feeding and locomotion behaviors. The second experiment determined agreement for visual observations and the RumiWatch noseband sensor from 17.75 h of rumination cycle and grazing bites. Pearson correlations evaluated associations for visual observations, and the RumiWatch noseband sensor and pedometer and were 0.84 for rumination, 0.76 for grazing, 0.39 for drinking, 0.57 for other activities, 0.83 for standing, 0.91 for lying, and 0.38 for walking. Correlations for visual observations and rumination cycle and grazing bites were -0.13 and 0.47, respectively. The RumiWatch system evaluated rumination, grazing, standing, and lying behaviors with high precision and accuracy, and the RumiWatch system may be used as a benchmark instead of visual observation to validate animal behavior technologies.

Herd life, lifetime production, and profitability of Viking Red-sired and Montbéliarde-sired crossbred cows compared with their Holstein herdmates.

The first 2 generations from a 3-breed rotation of the Viking Red (VR), Montbéliarde (MO), and Holstein (HO) breeds were compared with their HO herdmates in high-performance commercial herds in Minnesota. The designed study enrolled pure HO females in 2008 to initiate a comparison of 3-breed rotational crossbreds with their HO herdmates. Sires of cows were proven artificial insemination bulls selected for high genetic merit in each of the 3 breeds. The first-generation cows calved for a first time from 2010 to 2014 and had 376 VR × HO and 358 MO × HO crossbreds to compare with their 640 HO herdmates. The second-generation cows calved for a first time from 2012 to 2014 and had 109 VR × MO/HO and 117 MO × VR/HO crossbreds to compare with their 250 HO herdmates. Collection of data ceased on December 31, 2017, and all cows studied had the opportunity for 45 mo in the herd after first calving. Production of milk, fat, and protein (kg) during lifetimes of cows was estimated from test-day observations with best prediction. The lifetime profit function included revenue and cost. Revenue was from production, calves, and slaughter of cull cows. Costs included feed cost during lactation, lactating overhead cost, dry cow cost (including feed cost during dry periods), replacement cost, health treatment cost, insemination cost, fertility hormone cost, pregnancy diagnosis cost, hoof trimming cost, and carcass disposal cost. For individual cows with herd life longer than 45 mo after first calving, survival of cows was projected beyond 45 mo after first calving to estimate herd life, production, and profitability. The 2-breed crossbreds had +158 d longer herd life and the 3-breed crossbreds had +147 d longer herd life compared with their respective HO herdmates. Also, 12.4% of the 2-breed crossbreds died up to 45 mo after first calving compared with 16.3% of their HO herdmates. Furthermore, approximately 29% of both the 2-breed and 3-breed crossbreds lived beyond 45 mo after first calving compared with approximately 18% of their respective HO herdmates. On a lifetime basis, the 2-breed and 3-breed crossbreds provided +$122 and +$134, respectively, more cull cow revenue compared with their HO herdmates. For lifetime replacement cost, the 2-breed crossbreds did not differ from their HO herdmates; however, the 3-breed crossbreds had -$28 less lifetime replacement cost compared with their HO herdmates because of their younger age at first calving. The combined 2-breed crossbreds had +$0.473 (+13%) more daily profit (ignoring potential differences for feed efficiency) and the combined 3-breed crossbreds had +$0.342 (+9%) more daily profit compared with their respective HO herdmates. This resulted in +$173 more profit/cow annually for the combined 2-breed crossbreds and +$125 more profit/cow annually for the combined 3-breed crossbreds compared with their respective HO herdmates.

Genetic parameters for stayability of Holsteins in US organic herds.

The objectives of this study were to estimate genetic parameters for stayability in US organic Holstein dairy cows and estimate genetic correlations with nationally evaluated traits of interest. Stayability is the binary trait for success or failure to remain in the herd until a given time point. We used birth, calving, and cull dates from 16 USDA certified organic farms recommended by industry personnel as herds maintaining individual cow records and using artificial insemination. Stayability at 5 time points was assigned based on the presence of a calving date for each parity up to 5 (STAY1 to STAY5). We also considered livebirth (vs. stillbirth), stayability from a successful first calving to second calving (STAY12), stayability from a successful second calving to third calving (STAY23), and stayability as a repeated measure encompassing STAY1 to STAY5. In total, 44,995 females were used in this study. Ninety-six percent were born alive and of these, 64% reached first parity. Animals with Holstein sires and no other identified breed for 3 generations on the maternal side were included. Heritabilities for stayability to each parity on the underlying scale were estimated using a threshold model with the fixed effect of herd and the random effects of animal and herd-year-season of birth. Genetic correlations were estimated among livebirth, STAY1, STAY12, and STAY23 with a 4-trait linear model with fixed herd-year-season of birth and random effects of animal, dam of the calf (livebirth), and herd calving date (STAY12 and STAY23). Heritabilities for stayability ranged from 0.07 to 0.15 and was 0.08 for the direct effect of livebirth and 0.06 for the maternal effect of livebirth. The repeatability for stayability was 0.60. Genetic correlations ranged from 0.11 between livebirth and STAY1 to 0.83 between STAY12 and STAY23. Excluding livebirth, stayability to all time points was significantly correlated with productive life and with cow livability. In general, stayability was positively associated with milk yield and negatively associated with fat percent and stillbirth. In conclusion, stayability in organic herds is heritable and positively associated with nationally evaluated longevity traits suggesting that improvement for stayability in organic herds can be achieved with current national evaluations for longevity.

Genetic parameters and genomic regions associated with horn fly resistance in organic Holstein cattle.

Horn flies (Haematobia irritans [L.]) contribute to major economic losses of pastured cattle operations, particularly in organic herds because of limitations on control methods that can be used. The objectives of this research were to determine if resistance to horn flies is a heritable trait in organic Holstein cattle, determine associations with yield traits, and to detect genomic regions associated with fly infestation. Observations of fly load were recorded from 1,667 pastured Holstein cows, of which 640 were genotyped, on 13 organic dairies across the United States. Fly load score was determined using a 0 to 4 scale based on fly coverage from chine to loin on one side of the body, with 0 indicating few to no flies and 4 indicating high infestation. The scoring system was validated by counting flies from photographs taken at the time of scoring from 252 cows. To mitigate the effect of our data structure on potential selection bias effects on genetic parameter estimates, survival to subsequent lactations of scored animals and herd-mates that had been culled before the trial was accounted for as the trait stayability. Genetic parameters were estimated using single-step genomic analysis with 3-trait mixed models that included fly score, stayability, and a third phenotype. Model effects differed by variable, but fixed effects generally included a contemporary group, scorer, parity, and stage of lactation; random effects included animal, permanent environment, and residual error. A genome-wide association study was performed by decomposing estimated breeding values into marker effects to detect significant genomic regions associated with fly score. The rank correlation between the subjective fly score and the objective count was 0.79. The average heritability of fly score (± standard error) estimated across multiple models was 0.25 ± 0.04 when a known Holstein maternal grandsire was required and 0.19 ± 0.03 when only a known Holstein sire was required. Genetic correlation estimates with yield traits were moderately positive, but a greater fly load was associated with reduced yield after accounting for genetic merit. Lower fly loads were associated with white coat coloration; a significant genomic region on Bos taurus autosome 6 was identified that contains the gene KIT, which was the most plausible candidate gene for fly resistance because of its role in coat pattern and coloration. The magnitude of heritable variation in fly infestation is similar to other traits included in selection programs, suggesting that producers can select for resistance to horn flies.

Investigating the cow skin and teat canal microbiomes of the bovine udder using different sampling and sequencing approaches.

There is a need for standardized, efficient, and practical sampling methods to support large population-based studies of the internal and external epithelial microbiomes of the bovine udder. The primary objective of this study was to evaluate different sampling devices for the isolation of microbial DNA originating from the internal and external teat epithelium. Secondary objectives were to survey and compare the microbial diversity of external and teat canal epithelial microbiomes using amplicon and shotgun metagenomic sequencing approaches. To address these objectives, we enrolled a convenience sample of 24 Holstein dairy cows and collected samples from the external epithelium at the base of udder, the external teat barrel epithelium, the external teat apex epithelium, and the teat canal epithelium. Extracted DNA was quantified and subjected to PCR amplification of the V4 hypervariable region of the 16S rRNA gene and sequenced on the Illumina MiSeq platform (Illumina Inc., San Diego, CA). A subset of samples was subjected to a shallow shotgun metagenomic assay on the Illumina HiSeq platform. For samples collected from the external teat epithelium, we found that gauze squares consistently yielded more DNA than swabs, and Simpson's reciprocal index of diversity was higher for gauze than for swabs. The teat canal epithelial samples exhibited significantly lower diversity than the external sampling locations, but there were no significant differences in diversity between teat apex, teat barrel, and base of the udder samples. There were, however, differences in the microbial distribution and abundances of specific bacteria across external epithelial surfaces. The proportion of shotgun sequence reads classified as Bos taurus was highly variable between sampling locations, ranging from 0.33% in teat apex samples to 99.91% in teat canal samples. These results indicate that gauze squares should be considered for studying the microbiome of the external epithelium of the bovine udder, particularly if DNA yield must be maximized. Further, the relative proportion of host to non-host DNA present in samples collected from the internal and external teat epithelium should be considered when designing studies that utilize shotgun metagenomic sequencing.

Fertility and 305-day production of Viking Red-, Montbéliarde-, and Holstein-sired crossbred cows compared with Holstein cows during their first 3 lactations in Minnesota dairy herds.

Three generations of crossbred cows from a 3-breed rotation of the Viking Red (VR), Montbéliarde (MO), and Holstein (HO) breeds were compared with HO herdmates in 7 high-performance, commercial dairy herds in Minnesota. The designed study was initiated in 2008 with the enrollment of 3,550 HO females. Sires of cows were proven artificial insemination bulls and were high-ranking for genetic merit within each of the VR, MO, and HO breeds. The first generation of cows calved a first time from 2010 to 2017 and consisted of 644 VR × HO and 616 MO × HO 2-breed crossbreds and their 1,405 HO herdmates. The second generation calved a first time from 2012 to 2017 and consisted of 615 VR × MO/HO and 568 MO × VR/HO crossbreds and their 1,462 HO herdmates. The third generation calved a first time from 2014 to 2017 and was composed of 466 HO × VR/MO/HO and HO × MO/VR/HO crossbreds combined and their 736 HO herdmates. Collection of data ceased on December 31, 2017, and for the duration of study, many of the 2-breed and 3-breed crossbreds and their HO herdmates had the opportunity to complete at least 3 lactations, whereas the HO-sired crossbreds and their HO herdmates had the opportunity to complete 1 lactation. The 305-d actual (not mature equivalent) production of milk, fat, and protein was estimated from test-day observations with best prediction and cows with <305 DIM were projected to 305 d. The 2-breed and 3-breed crossbred cows had superiority over their HO herdmates for all of the fertility traits measured. The 2-breed crossbreds had -9, -17, and -15 d fewer days open than their HO herdmates during first, second, and third lactation, respectively. Also, the 3-breed crossbreds had -15, -19, and -20 d fewer days open than their HO herdmates during first, second, and third lactation, respectively. Cows in these herds had young ages at first calving of 22 to 23 mo across the breed groups, and the 3-breed crossbreds had significantly younger ages at first, second, and third calving than their HO herdmates. The 2-breed crossbreds had +2% higher fat plus protein production (kg) during first lactation, but did not differ during second and third lactation, from their HO herdmates. The 3-breed crossbreds had -3% to -4% lower fat plus protein production (kg) than their HO herdmates in each of their first 3 lactations. During first lactation, the HO-sired crossbreds did not differ for fat plus protein production (kg) from their HO herdmates.

Exploring physiological and genetic variation of digital cushion thickness in Holstein and Jersey cows and bulls.

The objectives of our study were to expand phenotypic characterization of digital cushion thickness (DCT) to the Jersey breed and include mature bulls and to identify breed-specific quantitative trait loci (QTL) for DCT within the Jersey or Holstein breeds and common QTL across breeds and sexes to better understand biological regulation and guide future use of marker-assisted selection. In a cohort of 698 cows and 85 bulls (Holstein and Jersey) from 8 farms in New York State, data were collected on DCT, body condition score (BCS), withers height, and sacral height. All animals underwent sonographic examination of the digital cushion evaluated at the sole ulcer site for the right front and hind feet. Linear mixed models were conducted on DCT separately for cows and bulls with fixed effects of time point, breed, age, digit, and BCS group. The models included random effects to control for the random subset of animals per farm, repeated measures, and multiple measurements from each animal. The phenotypic results indicated that DCT varied by time point, breed, age, digit, and BCS group for cows and by breed, age, digit, and BCS group for bulls. For the genotypic study, 616 cow DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA), whereas 76 bull DNA samples were genotyped on different platforms ranging from 5K to 150K. Multiple genome-wide association studies were conducted to highlight pertinent phenotyping parameters and genetic markers for genomic selection. Data were separated into 8 data sets based on different combinations of breed and sex. Each data set was assessed for quality of markers and samples before conducting genome-wide association studies for DCT, testing the inheritance models and genetic variation of digit, foot, and average thickness. Ten markers passed the Bonferroni correction threshold and 9 passed false discovery rate from 10 genome-wide association studies using a combination of the covariates breed, sex, genotyping batch plate, age, BCS, withers height, and sacral height. Of the 43 candidate genes, 8 novel biologically plausible genes were identified on Bos taurus autosomes 3, 4, 7, and 9: SFRS18 and LRRFIP1 function in fat deposition, whereas AHR, BZW2, EFNA5, USP45, and VAV3 effect bone growth, and SOSTDC1 is related to epidermal keratinocyte function. The genetic markers associated with DCT in this study were explored for variation between cows and bulls within and across breeds for their potential use in marker-assisted selection.

Health treatment cost, stillbirth, survival, and conformation of Viking Red-, Montbéliarde-, and Holstein-sired crossbred cows compared with pure Holstein cows during their first 3 lactations.

Three generations of crossbreds from a 3-breed rotation of the Viking Red (VR), Montbéliarde (MO), and Holstein (HO) breeds were compared with their HO herdmates in 7 commercial dairy herds in Minnesota. The designed study enrolled 3,550 HO females in 2008 to initiate crossbreeding and a control of pure HO herdmates within each herd. Service sires were high-ranking, proven AI bulls selected for high genetic merit within each of the VR, MO, and HO breeds. Cows in this study calved from 2010 to 2017 and collection of data ended on December 31, 2017. The first generation of cows consisted of 644 VR × HO and 616 MO × HO crossbreds and their 1,405 HO herdmates. The second generation had 615 VR × MO/HO and 568 MO × VR/HO crossbreds and their 1,462 HO herdmates. The third generation had 466 combined HO × VR/MO/HO and HO × MO/VR/HO crossbreds and their 736 HO herdmates. Total health cost was the sum of veterinary treatment cost, pharmaceutical cost, and farm labor cost to treat 16 different health disorders. Conformation traits and body condition score were subjectively scored once during early lactation for each of the first 3 lactations of cows. Total health cost of the 2-breed crossbreds was significantly lower during first (-23%), second (-29%), and third (-21%) lactation compared with their HO herdmates. For the 3-breed crossbreds, total health cost did not differ during first lactation but was -26% lower during both second and third lactation compared with their HO herdmates. The stillbirth rate for calves born to 2-breed crossbred dams (4%) was significantly lower compared with calves born to their HO herdmates (8%) at first calving. Survival from first to third calving (+9%) and first to fourth calving (+11%) was significantly higher for the 2-breed crossbreds compared with their HO herdmates. Also, the 3-breed crossbreds had significantly higher survival to third (+11%) and fourth (+19%) calving compared with their HO herdmates. Across each generation of crossbreeding, the crossbreds had uniformly shorter stature, less angularity, and less body depth compared with their respective HO herdmates. The crossbred cows also had significantly less udder clearance from the hock but significantly more rear teat width and longer teat length compared with their respective HO herdmates. Furthermore, the crossbred cows had higher body condition score compared with their HO herdmates during each of their first 3 lactations.

Validation of an ear tag-based accelerometer system for detecting grazing behavior of dairy cows.

The objective of the study was to develop a grazing algorithm for an ear tag-based accelerometer system (Smartbow GmbH, Weibern, Austria) and to validate the grazing algorithm with data from a noseband sensor. The ear tag has an acceleration sensor, a radio chip, and temperature sensor for calibration and it can monitor rumination and detect estrus and localization. To validate the ear tag, a noseband sensor (RumiWatch, Itin and Hoch GmbH, Liestal, Switzerland) was used. The noseband sensor detects pressure and acceleration patterns, and, with a software program specific to the noseband, pressure and acceleration patterns are used to classify data into eating, ruminating, drinking, and other activities. The study was conducted at the University of Minnesota West Central Research and Outreach Center (Morris, MN) and at Teagasc Animal and Grassland Research and Innovation Centre (Moorepark, Fermoy, Co. Cork, Ireland). During May and June 2017, observational data from Minnesota and Ireland were used to develop the grazing algorithm. During September 2018, data were collected by the ear tag and noseband sensor from 12 crossbred cows in Minnesota for a total of 248 h and from 9 Holstein-Friesian cows in Ireland for a total of 248 h. A 2-sided t-test was used to compare the percentage of grazing and nongrazing time recorded by the ear tag and the noseband sensor. Pearson correlations and concordance correlation coefficients (CCC) were used to evaluate associations between the ear tag and noseband sensor. The percentage of total grazing time recorded by the ear tag and by the noseband sensor was 37.0% [95% confidence interval (CI): 32.1 to 42.0] and 40.5% (95% CI: 35.5 to 45.6), respectively, in Minnesota, and 35.4% (95% CI: 30.6 to 40.2) and 36.9% (95% CI: 32.1 to 41.8), respectively, in Ireland. The ear tag and noseband sensor agreed strongly for monitoring grazing in Minnesota (r = 0.96; 95% CI: 0.94 to 0.97, CCC = 0.95) and in Ireland (r = 0.92; 95% CI: 0.90 to 0.94, CCC = 0.92). The results suggest that there is potential for the ear tag to be used on pasture-based dairy farms to support management decision-making.

Phenotypic characterization and genome-wide association studies of digital cushion thickness in Holstein cows.

The bovine digital cushion is a compression pad between the distal phalanx and sole and has been associated with claw horn disruption lesions. Digital cushion thickness (DCT) is estimated to be moderately heritable. Therefore, the objectives of our study were to examine influences of management and environment on DCT and to identify genetic markers and candidate genes associated with DCT. In a cohort of 502 Holsteins from 5 farms in New York State, DCT and body condition score (BCS) were collected twice, at <137 d prepartum and from 86 to 127 d in milk, corresponding to periods when the digital cushion is thickest and thinnest, respectively, as determined by previous research. Cows underwent sonographic examination of the digital cushion evaluated at the typical sole ulcer site for the right front and hind foot. Linear mixed models were conducted on DCT with the fixed effects of time point, digit, wither height, sacral height, BCS group, and multiple farm system variables separately and included random effects to control for the random subset of cows per farm, repeated measures, and multiple measurements from each cow. The phenotypic results indicated that DCT varied by sample time point, sacral height, parity, digit, BCS group, and wither height. For the genotypic study, 447 DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA). Quality assessment of markers and samples provided a final data set of 431 samples and 579,449 markers. Genome-wide association studies were conducted for DCT testing inheritance models and genetic variation of digit, foot, time point, and average thickness. One marker passed the Bonferroni correction threshold and 26 passed false discovery rate from 4 genome-wide association studies with covariates of sequencing batch plate, parity group, BCS, wither height, and sacral height. Ten candidate genes were identified, with 2 genes on Bos taurus autosomes 24 and 29 involved in biological functions related to the digital cushion: MC4R and DLG2 were related to fat deposition and bone growth, respectively. The genetic markers discovered in this study have the opportunity to be used in breeding programs using genomic selection to select against claw horn disruption lesions and lameness due to associations between the markers and DCT. Further studies on the biologically plausible candidate genes may identify causative genetic variants and how they relate to DCT through gene regulation, expression, structure, or copy number variation.

Short communication: Effects of mesh leggings on fly pressure and fly avoidance behaviors of pastured dairy cows.

Ectoparasitic stable flies (Stomoxys calcitrans [L.]), horn flies (Haematobia irritans [L.]), and face flies (Musca autumnalis De Geer) negatively affect dry matter intake, milk production, and health of pastured dairy cows. These flies cause fly avoidance behaviors and are a major welfare concern for dairy producers. The objective of this study was to evaluate the effects of mesh Shoofly Leggins (Stone Manufacturing & Supply, Kansas City, MO) on fly avoidance behaviors and numbers of flies attacking pastured dairy cows. In a crossover design, lactating dairy cows (n = 80) were randomly assigned to groups with and without leggings (Shoofly Leggins worn on all legs). All cows were managed in one group. Cows were observed for 2-wk periods, and then treatments were reversed in the next 2-wk interval. Counts of stable flies, horn flies, and face flies on all cows were recorded twice daily (once in morning per cow: 0930 to 1230 h; and once in the afternoon per cow: 1330 to 1630 h), 3 times per wk on Tuesday, Wednesday, and Thursday of each week during the study period. The total number of flies per cow was greater on cows with leggings than cows without leggings. The number of horn flies per cow was greater on cows in with leggings (26.4 flies/side) compared with cows without leggings (24.1 flies/side). Stable fly numbers were similar for cows with and without leggings (12.8 flies/leg). A random subset of 20 focal cows per group was observed during 5-min intervals to record frequencies of 4 behaviors: leg stomps, head tosses, skin twitches, and tail swishes. Counts of head tosses (2.6 vs. 3.1), skin twitches (20.9 vs. 19.6), and tail swishes (21.3 vs. 19.3) were similar for cows without leggings versus cows with leggings, respectively. However, foot stomps were 39% lower for cows with leggings compared with cows without leggings, and leg stomps were 26% higher in the afternoon than in the morning (2.9 vs. 2.4, respectively). A positive correlation was observed between stable and horn flies and all insect avoidance behaviors. Numbers of stable flies were 1.5 times greater in the afternoon than in the morning. The results of this study indicated that flies were associated with cow fly avoidance behaviors regardless of the use of leggings, but leggings effectively reduced foot stomps by 39%, so their use may provide some relief from stable fly injury to pastured dairy cows.

Three-breed rotational crossbreds of Montbéliarde, Viking Red, and Holstein compared with Holstein cows for feed efficiency, income over feed cost, and residual feed intake.

Rotational 3-breed crossbred cows of Montbéliarde, Viking Red, and Holstein (CB) were compared with Holstein (HO) cows for alternative measures of feed efficiency as well as income over feed cost (IOFC) and residual feed intake (RFI) during the first 150 d of first, second, and third lactations. Primiparous and multiparous CB (n = 63 and n = 43, respectively) and HO (n = 60 and n = 37, respectively) cows were fed the same total mixed ration twice daily with refusals weighed once daily. Feed was analyzed for dry matter content, net energy for lactation, and crude protein content. Body weight (BW) was recorded twice weekly. Daily production of milk, fat, and protein were estimated from monthly test days with best prediction. Measures of efficiency from 4 to 150 d in milk (DIM) were feed conversion efficiency (FCE), defined as fat plus protein production (kg) per kilogram of dry matter intake (DMI); ECM/DMI, defined as kilograms of energy-corrected milk (ECM) per kilogram of DMI; net energy for lactation efficiency (NELE), defined as ECM (kg) per megacalorie of net energy for lactation intake; crude protein efficiency (CPE), defined as true protein production (kg) per kilogram of crude protein intake; and DMI/BW, defined as DMI (kg) per kilogram of BW. The IOFC was defined as revenue from fat plus protein production minus feed cost. The RFI from 4 to 150 DIM for each lactation was the residual error remaining from regression of DMI on milk energy output (Mcal), metabolic BW, and energy required for change in BW (Mcal). Statistical analysis of measures of feed efficiency and RFI for primiparous cows included the fixed effects of year of calving and breed group. For multiparous cows, statistical analysis included breed as a fixed effect and cow as a repeated effect nested within breed group. Primiparous CB cows had higher means for FCE (+5.5%), ECM/DMI (+4.0%), NELE (+4.0%), and CPE (+5.2%) and a lower mean DMI/BW (-5.3%) than primiparous HO cows. Primiparous CB cows ($875) also had higher mean IOFC than primiparous HO cows ($825). In addition, mean RFI from 4 to 150 DIM was significantly lower (more desirable) for primiparous CB cows than HO cows. Likewise, multiparous CB cows had higher means for FCE (+8.2%), ECM/DMI (+5.9%), NELE (+5.8%), and CPE (+8.1%) and a lower mean for DMI/BW (-4.8%) than multiparous HO cows. Multiparous CB cows ($1,296) also had a higher mean for IOFC than multiparous HO cows ($1,208) and a lower mean for RFI from 4 to 150 DIM than HO cows.