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.

Consistency of dry matter intake in Holstein cows: Heritability estimates and associations with feed efficiency.

Resilience can be defined as the capacity to maintain performance or bounce back to normal functioning after a perturbation, and studying fluctuations in daily feed intake may be an effective way to identify resilient dairy cows. Our goal was to develop new phenotypes based on daily dry matter intake (DMI) consistency in Holstein cows, estimate genetic parameters and genetic correlations with feed efficiency and milk yield consistency, and evaluate their relationships with production, longevity, health, and reproduction traits. Data consisted of 397,334 daily DMI records of 6,238 lactating Holstein cows collected from 2007 to 2022 at 6 research stations across the United States. Consistency phenotypes were calculated based on the deviations from expected daily DMI for individual cows during their respective feeding trials, which ranged from 27 to 151 d in duration. Expected values were derived from different models, including simple average, quadratic and cubic quantile regression with a 0.5 quantile, and locally estimated scatterplot smoothing (LOESS) regression with span parameters 0.5 and 0.7. We then calculated the log of variance (log-Var-DMI) of daily deviations for each model as the consistency phenotype. Consistency of milk yield was also calculated, as a reference, using the same methods (log-Var-Milk). Genetic parameters were estimated using an animal model, including lactation, days in milk and cohort as fixed effects, and animal as random effect. Relationships between log-Var-DMI and traits currently considered in the US national genetic evaluation were evaluated using Spearman's rank correlations between sires' breeding values. Heritability estimates for log-Var-DMI ranged from 0.11 ± 0.02 to 0.14 ± 0.02 across models. Different methods (simple average, quantile regressions, and LOESS regressions) used to calculate log-Var-DMI yielded very similar results, with genetic correlations ranging from 0.94 to 0.99. Estimated genetic correlations between log-Var-DMI and log-Var-Milk ranged from 0.51 to 0.62. Estimated genetic correlations between log-Var-DMI and feed efficiency ranged from 0.55 to 0.60 with secreted milk energy, from 0.59 to 0.63 with metabolic body weight, and from 0.26 to 0.31 with residual feed intake (RFI). Relationships between log-Var-DMI and the traits in the national genetic evaluation were moderate and positive correlations with milk yield (0.20 to 0.21), moderate and negative correlations with female fertility (-0.07 to -0.20), no significant correlations with health and longevity, and favorable correlations with feed efficiency (-0.23 to -0.25 with feed saved and 0.21 to 0.26 with RFI). We concluded that DMI consistency is heritable and may be an indicator of resilience. Cows with lower variation in the difference between actual and expected daily DMI (more consistency) may be more effective in maintaining performance in the face of challenges or perturbations, whereas cows with greater variation in observed versus expected daily DMI (less consistency) are less feed efficient and may be less resilient.

Impact of parity differences on residual feed intake estimation in Holstein cows.

Residual feed intake (RFI) has been used as a measure of feed efficiency in farm animals. In lactating dairy cattle, RFI is typically obtained as the difference between dry matter intake observations and predictions from regression on known energy sinks, and effects of parity, days in milk, and cohort. The impact of parity (lactation number) on the estimation of RFI is not well understood, so the objectives of this study were to (1) evaluate alternative RFI models in which the energy sinks (metabolic body weight, body weight change, and secreted milk energy) were nested or not nested within parity, and (2) estimate variance components and genetic correlations for RFI across parities. Data consisted of 72,474 weekly RFI records of 5,813 lactating Holstein cows collected from 2007 to 2022 in 5 research stations across the United States. Estimates of heritability, repeatability, and genetic correlations between weekly RFI for parities 1, 2, and 3 were obtained using bivariate repeatability animal models. The nested RFI model showed better goodness of fit than the nonnested model, and some partial regression coefficients of dry matter intake on energy sinks were heterogeneous between parities. However, the Spearman's rank correlation between RFI values calculated from nested and nonnested models was equal to 0.99. Similarly, Spearman's rank correlation between the RFI breeding values from these 2 models was equal to 0.98. Heritability estimates for RFI were equal to 0.16 for parity 1, 0.19 for parity 2, and 0.22 for parity 3. Repeatability estimates for RFI across weeks within parities were high, ranging from 0.51 to 0.57. Spearman's rank correlations of sires' breeding values were 0.99 between parities 1 and 2, 0.91 between parities 1 and 3, and 0.92 between parities 2 and 3. We conclude that nesting energy sinks within parity when computing RFI improves model goodness of fit, but the impact on the estimated breading values appears to be minimal.

Estimates of genetic parameters for feeding behavior traits and their associations with feed efficiency in Holstein cows.

Residual feed intake (RFI) is commonly used to measure feed efficiency but individual intake recording systems are needed. Feeding behavior may be used as an indicator trait for feed efficiency using less expensive precision livestock farming technologies. Our goal was to estimate genetic parameters for feeding behavior and the genetic correlations with feed efficiency in Holstein cows. Data consisted of 75,877 daily feeding behavior records of 1,328 mid-lactation Holstein cows in 31 experiments conducted from 2009 to 2020 with an automated intake recording system. Feeding behavior traits included number of feeder visits per day, number of meals per day, duration of each feeder visit, duration of each meal, total duration of feeder visits, intake per visit, intake per meal [kg of dry matter (DM)], feeding rate per visit, and feeding rate per meal (kg of DM per min). The meal criterion was estimated as 26.4 min, which means that any pair of feeder visits separated by less than 26.4 min were considered part of the same meal. The statistical model included lactation and days in milk as fixed effects, and experiment-treatment, animal, and permanent environment as random effects. Genetic parameters for feeding behavior traits were estimated using daily records and weekly averages. Estimates of heritability for daily feeding behavior traits ranged from 0.09 ± 0.02 (number of meals; mean ± standard error) to 0.23 ± 0.03 (feeding rate per meal), with repeatability estimates ranging from 0.23 ± 0.01 (number of meals) to 0.52 ± 0.02 (number of feeder visits). Estimates of heritability for weekly averages of feeding behavior traits ranged from 0.19 ± 0.04 (number of meals) to 0.32 ± 0.04 (feeding rate per visit), with repeatability estimates ranging from 0.46 ± 0.02 (duration of each meal) to 0.62 ± 0.02 (feeding rate per visit and per meal). Most of the feeding behavior measures were strongly genetically correlated, showing that with more visits or meals per day, cows spend less time in each feeder visit or meal with lower intake per visit or meal. Weekly averages for feeding behavior traits were analyzed jointly with RFI and its components. Number of meals was genetically correlated with milk energy (0.48), metabolic body weight (-0.27), and RFI (0.19). Duration of each feeder visit and meal were genetically correlated with milk energy (0.43 and 0.44, respectively). Total duration of feeder visits per day was genetically correlated with DM intake (0.29), milk energy (0.62), metabolic body weight (-0.37), and RFI (0.20). Intake per visit and meal were genetically correlated with DM intake (0.63 and 0.87), milk energy (0.47 and 0.69), metabolic body weight (0.47 and 0.68), and RFI (0.31 and 0.65). Feeding rate was genetically correlated with DM intake (0.69), metabolic body weight (0.67), RFI (0.47), and milk energy (0.21). We conclude that measures of feeding behavior could be useful indicators of dairy cow feed efficiency, and individual cows that eat at a slower rate may be more feed efficient.

Gene mapping, gene-set analysis, and genomic prediction of postpartum blood calcium in Holstein cows.

The onset of lactation results in a sudden irreversible loss of Ca for colostrum and milk synthesis. Some cows are unable to quickly adapt to this demand and succumb to clinical hypocalcemia, whereas a larger proportion of cows develop subclinical hypocalcemia that predisposes them to other peripartum diseases. The objective of this study was to perform a comprehensive genomic analysis of blood total Ca concentration in periparturient Holstein cows. We first performed a genomic scan and a subsequent gene-set analysis to identify candidate genes, biological pathways, and molecular mechanisms affecting postpartum Ca concentration. Then, we assessed the prediction of postpartum Ca concentration using genomic information. Data consisted of 7,691 records of plasma or serum concentrations of Ca measured in the first, second, and third day after parturition of 959 primiparous and 1,615 multiparous cows that calved between December 2015 and June 2020 in 2 dairy herds. All cows were genotyped with 80k SNPs. The statistical model included lactation (1 to 5+), calf category (male, females, twins), and day as fixed effects, and season-treatment-experiment, animal, and permanent environmental as random effects. Model predictive ability was evaluated using 10-fold cross-validation. Heritability and repeatability estimates were 0.083 (standard error = 0.017) and 0.444 (standard error = 0.028). The association mapping identified 2 major regions located on Bos taurus autosome (BTA)6 and BTA16 that explained 1.2% and 0.7% of additive genetic variance of Ca concentration, respectively. Interestingly, the region on BTA6 harbors the GC gene, which encodes the vitamin D binding protein, and the region on BTA16 harbors LRRC38, which is actively involved in K transport. Other sizable peaks were identified on BTA5, BTA2, BTA7, BTA14, and BTA9. These regions harbor genes associated with Ca channels (CACNA1S, CRACR2A), K channels (KCNK9), bone remodeling (LRP6), and milk production (SOCS2). The gene-set analysis revealed terms related to vitamin transport, calcium ion transport, calcium ion binding, and calcium signaling. Genomic predictions of phenotypic and genomic estimated breeding values of Ca concentration yielded predictive correlations up to 0.50 and 0.15, respectively. Overall, the present study contributes to a better understanding of the genetic basis of postpartum blood Ca concentration in Holstein cows. In addition, the findings may contribute to the development of novel selection and management strategies for reducing periparturient hypocalcemia in dairy cattle.

Genomic Analysis of Visceral Fat Accumulation in Holstein Cows.

Visceral fat is related to important metabolic processes, including insulin sensitivity and lipid mobilization. The goal of this study was to identify individual genes, pathways, and molecular processes implicated in visceral fat deposition in dairy cows. Data from 172 genotyped Holstein cows classified at slaughterhouse as having low (n = 77; omental fold < 5 mm in thickness and minimum fat deposition in omentum) or high (n = 95; omental fold ≥ 20 mm in thickness and marked fat deposition in omentum) omental fat were analyzed. The identification of regions with significant additive and non-additive genetic effects was performed using a two-step mixed model-based approach. Genomic scans were followed by gene-set analyses in order to reveal the genetic mechanisms controlling abdominal obesity. The association mapping revealed four regions located on BTA19, BTA20 and BTA24 with significant additive effects. These regions harbor genes, such as SMAD7, ANKRD55, and the HOXB family, that are implicated in lipolysis and insulin tolerance. Three regions located on BTA1, BTA13, and BTA24 showed marked non-additive effects. These regions harbor genes MRAP, MIS18A, PRNP and TSHZ1, that are directly implicated in adipocyte differentiation, lipid metabolism, and insulin sensitivity. The gene-set analysis revealed functional terms related to cell arrangement, cell metabolism, cell proliferation, cell signaling, immune response, lipid metabolism, and membrane permeability, among other functions. We further evaluated the genetic link between visceral fat and two metabolic disorders, ketosis, and displaced abomasum. For this, we analyzed 28k records of incidence of metabolic disorders from 14k cows across lactations using a single-step genomic BLUP approach. Notably, the region on BTA20 significantly associated with visceral fat deposition was also associated with the incidence of displaced abomasum. Overall, our findings suggest that visceral fat deposition in dairy cows is controlled by both additive and non-additive effects. We detected at least one region with marked pleiotropic effects affecting both visceral fat accumulation and displaced abomasum.

Genomic Study of Babesia bovis Infection Level and Its Association With Tick Count in Hereford and Braford Cattle.

Bovine babesiosis is a tick-borne disease caused by intraerythrocytic protozoa and leads to substantial economic losses for the livestock industry throughout the world. Babesia bovis is considered the most pathogenic species, which causes bovine babesiosis in Brazil. Genomic data could be used to evaluate the viability of improving resistance against B. bovis infection level (IB) through genomic selection, and, for that, knowledge of genetic parameters is needed. Furthermore, genome-wide association studies (GWAS) could be conducted to provide a better understanding of the genetic basis of the host response to B. bovis infection. No previous work in quantitative genetics of B. bovis infection was found. Thus, the objective of this study was to estimate the genetic correlation between IB and tick count (TC), evaluate predictive ability and applicability of genomic selection, and perform GWAS in Hereford and Braford cattle. The single-step genomic best linear unbiased prediction method was used, which allows the estimation of both breeding values and marker effects. Standard phenotyping was conducted for both traits. IB quantifications from the blood of 1,858 animals were carried using quantitative PCR assays. For TC, one to three subsequent tick counts were performed by manually counting adult female ticks on one side of each animal's body that was naturally exposed to ticks. Animals were genotyped using the Illumina BovineSNP50 panel. The posterior mean of IB heritability, estimated by the Bayesian animal model in a bivariate analysis, was low (0.10), and the estimations of genetic correlation between IB and TC were also low (0.15). The cross-validation genomic prediction accuracy for IB ranged from 0.18 to 0.35 and from 0.29 to 0.32 using k-means and random clustering, respectively, suggesting that genomic predictions could be used as a tool to improve genetics for IB, especially if a larger training population is developed. The top 10 single nucleotide polymorphisms from the GWAS explained 5.04% of total genetic variance for IB, which were located on chromosomes 1, 2, 5, 6, 12, 17, 18, 16, 24, and 26. Some candidate genes participate in immunity system pathways indicating that those genes are involved in resistance to B. bovis in cattle. Although the genetic correlation between IB and TC was weak, some candidate genes for IB were also reported in tick infestation studies, and they were also involved in biological resistance processes. This study contributes to improving genetic knowledge regarding infection by B. bovis in cattle.

Selection signatures in candidate genes and QTL for reproductive traits in Nellore heifers.

The identification of selection signature genes may help to detect genomic regions that underwent artificial selection and contributed to phenotypic diversity. The aim of this study, therefore, was to detect selection signatures in candidate genes and quantitative trait locus (QTL) for reproductive traits in a Nellore population being selected for sexual precocity. A total of 2035 Nellore heifers, sourced from breeding programs focused on sexual precocity, were used. Candidate genes and some specific QTL related to reproductive traits were chosen based on published literature and Animal QTL databases, respectively, for investigation whether these regions were affected by selection. Selection signature DNA sequences were detected in the selected regions using the extended haplotype homozygosity (EHH) and relative extended haplotype homozygosity (REHH) methods. From 22,241 single nucleotide polymorphisms (SNPs) located in the candidate genes and QTL, 17,312 SNPs generated 2756 haplotype blocks. A total of 7518 EHH tests were analyzed using haplotypes with a frequency of more than 25%, for which there were 39 tests that were significant for REHH (P<0.01). Selection signature DNA sequences were detected that contained several QTLs for important reproductive traits in cattle, suggesting that reproductive traits may have been affected by selection for sexual precocity in this population. Forty-six genes were located in the selection signature regions, whereas 24 genes participated in important biological processes or pathways that may underlie sexual precocity. These results indicate there are possible molecular mechanisms related to sexual precocity in the Nellore breed.

Genome Association Study for Visual Scores in Nellore Cattle Measured at Weaning.

BACKGROUND: Genome-wide association studies (GWAS) are utilized in cattle to identify regions or genetic variants associated with phenotypes of interest, and thus, to identify design strategies that allow for the increase of the frequency of favorable alleles. Visual scores are important traits of cattle production in Brazil because they are utilized as selection criteria, helping to choose more harmonious animals. Despite its importance, there are still no studies on the genome association for these traits. This study aimed to identify genome regions associated with the traits of conformation, precocity and muscling, based on a visual score measured at weaning. RESULTS: Bayesian approaches with BayesC and Bayesian LASSO were utilized with 2873 phenotypes of Nellore cattle for a GWAS. The animals were genotyped with Illumina BovineHD BeadChip, and a total of 309,865 SNPs were utilized after quality control. In the analyses, phenotype and deregressed breeding values were utilized as dependent variables; a threshold model was utilized for the former and a linear model for the latter. The association criterion was the percentage of genetic variance explained by SNPs found in 1 Mb-long windows. The Bayesian approach BayesC was better adjusted to the data because it could explain a larger phenotypic variance for both dependent variables. CONCLUSIONS: There were no large effects for the visual scores, indicating that they have a polygenic nature; however, regions in chromosomes 1, 3, 5, 7, 14, 15, 16, 19, 20 and 23 were identified and explained a large part of the genetic variance.