Preference for feed bins shared with cows of the same or different parity and relationships with feeding behavior and feed efficiency.

Our objectives were to (1) evaluate cows' preferences for visiting feed bins limited to either same- versus mixed-parity social interactions, depending on their parity; (2) examine the effect of parity and bin social dynamic type on competition behavior and feeding patterns, and (3) investigate cow-level relationships between feed bunk competition behavior, feeding patterns, and feed efficiency. Twenty-eight primiparous and 28 multiparous (2.4 ± 0.6 lactations) lactating Holstein cows (127.8 ± 30.1 and 145.3. ± 10.4 DIM, respectively) were housed in a freestall pen with 28 roughage intake control bins (2:1 stocking density). Each cow was assigned to 2 bins, including 1 shared with 3 other cows of the same parity (SM) and 1 with 3 cows of mixed parities (MX, 50% primiparous and 50% multiparous). Feed bunk competition was recorded via video in the first hour after morning feed delivery for 2 d, and feeding patterns were recorded from 24-h roughage intake control data. Residual feed intake was calculated as the difference between predicted and observed dry matter intake after accounting for known energy sinks. Based on the first visit to the feed bunk after fresh feed delivery, multiparous cows tended to prefer the MX bin compared with the SM one; cows showed no other overall preference for bin type based on number of visits. Over time, multiparous cows remained consistent in their magnitude of preference for visiting each bin type, but involvement in competition was not consistent over time. Primiparous cows tended to be involved in more total competitive contacts and ate faster at the SM bin compared with the MX one. Those primiparous cows who visited the SM bin more often within the first hour after morning feed delivery tended to be less feed efficient. Multiparous cows initiated more successful replacements after a displacement at the MX versus SM bin, with no difference in feeding patterns between bin types. Regardless of parity or bin type, visiting the bunk sooner after feed delivery was correlated with involvement in more competitive interactions and more time eating within the first 30 min. Consuming more feed during a longer first visit to the bunk after fresh feed delivery was correlated with being less feed efficient. Overall, when given the choice of feeding from bins shared with cows of the same or mixed parities at a 2:1 stocking density, primiparous cows showed differences in behavior between those bin types, with implications for feed efficiency; these effects are perhaps an unintended consequence of compensatory strategies to avoid direct competition with multiparous cows.

Long-term growth, feed efficiency, enteric methane emission, and blood metabolite responses to in utero hyperthermia in Holstein heifers.

Dairy producers are experiencing production and animal welfare pressures from the increasing frequency and severity of heat stress events due to global climate change. Offspring performance during the preweaning and lactating periods is compromised when exposed to heat stress during late gestation (in utero). However, knowledge of the lingering effects of in utero heat stress on yearling dairy heifers is limited. Herein, we investigated the long-term effects of in utero heat stress on heifer growth, feed efficiency, and enteric methane emissions in postpubertal heifers. During the last 56 d of gestation, 38 pregnant cows carrying heifer calves were exposed to either heat stress (IUHT; n = 17) or artificial cooling (IUCL; n = 21). At 18 ± 1 mo of age, the resulting IUCL and IUHT heifers were enrolled in the present 63-d study. Heifers were blocked by weight and randomly assigned to 3 pens with Calan gates. Body weights were recorded on 3 consecutive days at the start and end of the trial and used to calculate ADG. Body condition score, hip width, body length, and chest girth were measured at the start and end of the study. All heifers were fed a TMR comprised of 46.6% oatlage, 44.6% grass/alfalfa haylage, 7.7% male-sterile corn silage, 0.3% urea, and 0.8% mineral/vitamin supplement (on a DM basis). The TMR and refusal samples were obtained daily, composited weekly, and dried to calculate DMI. During the study, each pen had access to a GreenFeed unit for 8 ± 1d to measure CH4 and CO2 gas fluxes. During the last 3 d of measuring CH4 and CO2 fluxes, fecal samples were collected, composited by animal, dried, and analyzed to calculate NDF, OM, and DM digestibility. On the last day of fecal sampling, blood samples were also collected via coccygeal venipuncture, and GC time-of-flight MS analysis was performed. Residual feed intake (RFI; predicted DMI - observed DMI), and feed conversion efficiency (FCE; DMI/ADG) were calculated to estimate feed efficiency. No differences were found in initial or final BW, hip width, chest girth, or BCS; however, IUCL heifers were longer in body length compared with IUHT heifers. Dry matter intake, ADG, RFI, and FCE were similar between IUHT and IUCL heifers. In utero heat-stressed and IUCL heifers produced similar amounts of CH4 and CO2, and no differences were found in the number of GreenFeed visits or latency to approach the GreenFeed. The concentrations of 6 blood metabolites involved in lipogenic pathways were different between in utero treatments. In conclusion, in utero heat stress does not seem to have long-term effects on feed efficiency or methane emissions during the postpubertal growing phase; however, IUCL heifers maintained a body-length advantage over their IUHT counterparts and differed in concentrations of several candidate metabolites that encourage further exploration of their potential function in key organs, such as the liver and mammary gland.

Host and rumen microbiome contributions to feed efficiency traits in Holstein cows.

It is now widely accepted that dairy cow performance is influenced by both the host genome and rumen microbiome composition. The contributions of the genome and the microbiome to the phenotypes of interest are quantified by heritability (h2) and microbiability (m2), respectively. However, if the genome and microbiome are included in the model, then the h2 reflects only the contribution of the direct genetic effects quantified as direct heritability (hd2), and the holobiont effect reflects the joint action of the genome and the microbiome, quantified as the holobiability (ho2). The objectives of this study were to estimate h2, hd2,m2, and ho2 for dry matter intake, milk energy, and residual feed intake; and to evaluate the predictive ability of different models, including genome, microbiome, and their interaction. Data consisted of feed efficiency records, SNP genotype data, and 16S rRNA rumen microbial abundances from 448 mid-lactation Holstein cows from 2 research farms. Three kernel models were fit to each trait: one with only the genomic effect (model G), one with the genomic and microbiome effects (model GM), and one with the genomic, microbiome, and interaction effects (model GMO). The model GMO, or holobiont model, showed the best goodness-of-fit. The hd2 estimates were always 10% to 15% lower than h2 estimates for all traits, suggesting a mediated genetic effect through the rumen microbiome, and m2 estimates were moderate for all traits, and up to 26% for milk energy. The ho2 was greater than the sum of hd2 and m2, suggesting that the genome-by-microbiome interaction had a sizable effect on feed efficiency. Kernel models fitting the rumen microbiome (i.e., models GM and GMO) showed larger predictive correlations and smaller prediction bias than the model G. These findings reveal a moderate contribution of the rumen microbiome to feed efficiency traits in lactating Holstein cows and strongly suggest that the rumen microbiome mediates part of the host genetic effect.

Investigating relationships between the host genome, rumen microbiome, and dairy cow feed efficiency using mediation analysis with structural equation modeling.

The rumen microbiome is crucial for converting feed into absorbable nutrients used for milk synthesis, and the efficiency of this process directly impacts the profitability and sustainability of the dairy industry. Recent studies have found that the rumen microbial composition explains part of the variation in feed efficiency traits, including dry matter intake, milk energy, and residual feed intake. The main goal of this study was to reveal relationships between the host genome, rumen microbiome, and dairy cow feed efficiency using structural equation models. Our specific objectives were to (i) infer the mediation effects of the rumen microbiome on feed efficiency traits, (ii) estimate the direct and total heritability of feed efficiency traits, and (iii) calculate the direct and total breeding values of feed efficiency traits. Data consisted of dry matter intake, milk energy, and residual feed intake records, SNP genotype data, and 16S rRNA rumen microbial abundances from 448 mid-lactation Holstein cows from 2 research farms. We implemented structural equation models such that the host genome directly affects the phenotype (GP → P) and the rumen microbiome (GM → P), while the microbiome affects the phenotype (M → P), partially mediating the effect of the host genome on the phenotype (G → M → P). We found that 7 to 30% of microbes within the rumen microbial community had structural coefficients different from zero. We classified these microbes into 3 groups that could have different uses in dairy farming. Microbes with heritability <0.10 but significant causal effects on feed efficiency are attractive for external interventions. On the other hand, 2 groups of microbes with heritability ≥0.10, significant causal effects, and genetic covariances and causal effects with the same or opposite sign to feed efficiency are attractive for selective breeding, improving or decreasing the trait heritability and response to selection, respectively. In general, the inclusion of the different microbes in genomic models tends to decrease the trait heritability rather than increase it, ranging from -15% to +5%, depending on the microbial group and phenotypic trait. Our findings provide more understanding to target rumen microbes that can be manipulated, either through selection or management interventions, to improve feed efficiency traits.

Genetic relationships between behavioral traits and feed efficiency traits in lactating Holstein cows.

The evaluation of dairy cow feed efficiency using residual feed intake accounts for known energy sinks. However, behavioral traits may also contribute to the variation in feed efficiency. Our objective was to estimate the heritability and repeatability of behavioral traits and their genetic correlations with feed efficiency and its components in lactating Holstein cows. The first data set consisted of 36,075 daily rumination and lying time records collected using a SMARTBOW ear tag accelerometer (Zoetis, Parsippany, NJ) and 6,371 weekly feed efficiency records of 728 cows from the University of Wisconsin-Madison. The second data set consisted of 59,155 daily activity records, measured as number of steps, recorded by pedometers (AfiAct; S.A.E. Afikim, Kibbutz Afikim, Israel), and 8,626 weekly feed efficiency records of 635 cows from the University of Florida. Feed efficiency and its components included dry matter intake, change in body weight, metabolic body weight, secreted milk energy, and residual feed intake. The statistical models included the fixed effect of cohort, lactation number, and days in milk, and the random effects of animal and permanent environment. Heritability estimates for behavioral traits using daily records were 0.19 ± 0.06 for rumination and activity, and 0.37 ± 0.07 for lying time. Repeatability estimates for behavioral traits using daily data ranged from 0.56 ± 0.02 for activity to 0.62 ± 0.01 for lying time. Both heritability and repeatability estimates were larger when weekly records instead of daily records were used. Rumination and activity had positive genetic correlations with residual feed intake (0.40 ± 0.19 and 0.31 ± 0.22, respectively) while lying time had a negative genetic correlation with this residual feed intake (-0.27 ± 0.11). These results indicate that more efficient cows tend to spend more time lying and less time active. Additionally, less efficient cows tend to eat more and therefore also tend to ruminate longer. Overall, sensor-based behavioral traits are heritable and genetically correlated with feed efficiency and its components and, therefore, they could be used as indicators to identify feed efficient cows within the herd.

Dry matter intake in US Holstein cows: Exploring the genomic and phenotypic impact of milk components and body weight composite.

Large datasets allow estimation of feed required for individual milk components or body maintenance. Phenotypic regressions are useful for nutrition management, but genetic regressions are more useful in breeding programs. Dry matter intake records from 8,513 lactations of 6,621 Holstein cows were predicted from phenotypes or genomic evaluations for milk components and body size traits. The mixed models also included DIM, age-parity subclass, trial date, management group, and BW change during 28- and 42-d feeding trials in mid lactation. Phenotypic regressions of DMI on milk (0.014 ± 0.006), fat (3.06 ± 0.01), and protein (4.79 ± 0.25) were much less than corresponding genomic regressions (0.08 ± 0.03, 11.30 ± 0.47, and 9.35 ± 0.87, respectively) or sire genomic regressions multiplied by 2 (0.048 ± 0.04, 6.73 ± 0.94, and 4.98 ± 1.75). Thus, marginal feed costs as fractions of marginal milk revenue were higher from genetic than phenotypic regressions. According to the ECM formula, fat production requires 69% more DMI than protein production. In the phenotypic regression, it was estimated that protein production requires 56% more DMI than fat. However, the genomic regression for the animal showed a difference of only 21% more DMI for protein compared with fat, whereas the sire genomic regressions indicated approximately 35% more DMI for fat than protein. Estimates of annual maintenance in kilograms DMI/kilograms BW per lactation were similar from phenotypic regression (5.9 ± 0.14), genomic regression (5.8 ± 0.31), and sire genomic regression multiplied by 2 (5.3 ± 0.55) and are larger than those estimated by the National Academies for Science, Engineering, and Medicine based on NEL equations. Multiple regressions on genomic evaluations for the 5 type traits in body weight composite (BWC) showed that strength was the type trait most associated with BW and DMI, agreeing with the current BWC formula, whereas other traits were less useful predictors, especially for DMI. The Net Merit formula used to weight different genetic traits to achieve an economically optimal overall selection response was revised in 2021 to better account for these estimated regressions. To improve profitability, breeding programs should select smaller cows with negative residual feed intake that produce more milk, fat, and protein.

Long-Term Effects of Pre-Weaning Individual or Pair Housing of Dairy Heifer Calves on Subsequent Growth and Feed Efficiency.

Our objective in this exploratory study was to evaluate the long-term impacts of pre-weaning social isolation vs. contact on subsequent growth and feed efficiency of Holstein heifers. As pre-weaned calves, 41 heifers were housed individually (n = 15 heifers) or in pairs (n = 13 pairs; 26 heifers). At 18 months of age, heifers were blocked by body weight and randomly assigned to one of three pens within a block (six to eight heifers per pen; six pens total), with original pairs maintained. Body weight (BW), hip height and width, and chest girth were measured at the start and end of the study. Each pen was given 3 days of access to a GreenFeed greenhouse gas emissions monitor to assess potential physiological differences between treatments in enteric methane emissions or behavioral differences in propensity to approach a novel object. During the 9-week study, heifers were fed a common diet containing 62.3% male-sterile corn silage, 36.0% haylage, 0.7% urea, and 1.0% mineral (DM basis). To calculate daily feed intake, as-fed weights and refusals were recorded for individual heifers using Calan gates. Feed samples were collected daily, composited by week, and dried to calculate dry matter intake (DMI). Feed refusal and fecal samples were collected on 3 consecutive days at 3 timepoints, composited by heifer, dried, and analyzed to calculate neutral detergent fiber (NDF), organic matter (OM), and DM digestibility. Feed efficiency was calculated as feed conversion efficiency (FCE; DMI/average daily gain [ADG]) and residual feed intake (RFI; observed DMI-predicted DMI). Paired and individually housed heifers did not differ in DMI, ADG, FCE, or RFI. Although no differences were found in initial or final hip height, hip width, or chest girth, heifers which had been pair-housed maintained a greater BW than individually housed heifers during the trial. Methane production, intensity, and yield were similar between treatments. Pre-weaning paired or individual housing did not impact the number of visits or latency to approach the GreenFeed; approximately 50% of heifers in each treatment visited the GreenFeed within 8 h of exposure. Digestibility of OM, DM, and NDF were also similar between housing treatments. In conclusion, pre-weaning pair housing had no adverse effects on growth, feed efficiency, or methane emissions at 18 to 20 months of age.

Identification of genomic regions and genes associated with subclinical ketosis in periparturient dairy cows.

Subclinical ketosis (SCK) is a prevalent metabolic disorder that occurs during the transition to lactation period. It is defined as a high blood concentration of ketone bodies (beta-hydroxybutyric acid f ≥ 1.2 mmol/L) within the first few weeks of lactation, and often presents without clinical signs. SCK is mainly caused by negative energy balance (NEB). The objective of this study is to identify single nucleotide polymorphisms (SNPs) associated with SCK using genome-wide association studies (GWAS), and to predict the biological functions of proximal genes using gene-set enrichment analysis (GSEA). Blood samples were collected from 112 Holstein cows between 5 and 18 days postpartum to determine the incidence of SCK. Genomic DNA extracted from both SCK and healthy cows was examined using the Illumina Bovine SNP50K BeadChip for genotyping. GWAS revealed 194 putative SNPs and 163 genes associated with those SNPs. Additionally, GSEA showed that the genes retrieved by Database for Annotation, Visualization, and Integrated Discovery (DAVID) belonged to calcium signaling, starch and sucrose, immune network, and metabolic pathways. Furthermore, the proximal genes were found to be related to germ cell and early embryo development. In summary, this study proposes several feasible SNPs and genes associated with SCK through GWAS and GSEA. These candidates can be utilized in selective breeding programs to reduce the genetic risk for SCK and subfertility in high-performance dairy cows.

Behavioral consistency of competitive behaviors and feeding patterns in lactating dairy cows across stocking densities at the feed bunk.

Introduction: High feed bunk stocking densities can differentially impact individual dairy cows' competitive behaviors, feeding patterns, and feed efficiency. Our objective was to manipulate feed bunk stocking densities to evaluate intra-individual behavioral consistency across stocking densities and quantify associations with feed efficiency and production. Methods: Thirty-two primiparous (130.7 ± 29.0 days in milk, DIM) and 32 multiparous (111.3 ± 38.3 DIM) lactating Holstein cows were housed with 32 roughage intake control (RIC) bins. Each cow was assigned to share 8 bins with others of the same parity and similar body weight (16 cows/block; 2:1 feed bunk stocking density except during tests). Competition and feeding patterns were evaluated via video and RIC data, respectively, at 3 stocking densities (1:1, 2:1, 4:1 cows/bin) during 1-h tests (2 tests/stocking density; 6 tests/cow) following 2 h feed deprivation. Residual feed intake (RFI) was calculated across the 45-d study as the difference between observed and predicted dry matter intake (DMI) after accounting for known energy sinks. Linear mixed models were used to evaluate the overall impact of test stocking density on competition and feeding patterns. To evaluate intra-individual consistency between stocking densities, individual stability statistic (ISS) scores were computed. Correlational relationships were determined between RFI and ISS scores. Results and dicsussion: Cows displayed the most competitive behaviors at 2:1 stocking density (p < 0.0001) but experienced the highest rate of contacts per minute of eating time at 4:1 (1:1 vs. 2:1 vs. 4:1: 0.09 vs. 0.95 vs. 1.60 contacts/min; p < 0.0001). Feeding patterns were modulated as stocking density increased; eating rate increased (0.16 vs. 0.18 vs. 0.22 kg/min; p < 0.001) as eating time (40.3 vs. 28.2 vs. 14.6 min; p < 0.001) and DMI decreased (6.3 vs. 5.1 vs. 3.0 kg; p < 0.001). As stocking density doubled, individuals remained consistent (p = 0.018) in time spent near others actively eating and tended to remain consistent in competition behavior and feeding patterns (0.053 ≤ p ≤ 0.094). Between 2:1 and 4:1, cows with higher DMI and milk production were more consistent in first-visit DMI and duration. Feed efficiency was not associated with behavioral consistency across the tests (p ≥ 0.14). Nonetheless, feed bunk stocking density has behavioral implications which should be considered in nutritional management decisions.