Assessment of associations between transition diseases and reproductive performance of dairy cows using survival analysis and decision tree algorithms

This study aimed to evaluate the associations between transition cow conditions and diseases TD with fertility in Holstein cows, and to compare analytic methods for doing so. Kaplan-Meier, Cox proportional hazard, and decision tree models were used to analyze the associations of TD with the pregnancy risk at 120 and 210 DIM from a 1-year cohort with 1946 calvings from one farm. The association between TD and fertility was evaluated as follows: 1 cows with TD whether complicated with another TD or not TD-all, versus healthy cows, and 2 cows with uncomplicated TD TD-single, versus cows with multiple TD TD+; complicated cases, versus healthy cows. The occurrence of twins, milk fever, retained placenta, metritis, ketosis, displaced abomasum, and clinical mastitis were recorded. Using Kaplan-Meier models, in primiparous cows the 120 DIM pregnancy risk was 62% (95% CI: 57-67 %) for healthy animals. This was not significantly different for TD-single (58%; 95% CI: 51-66 %) but was reduced for TD+ (45%; 95% CI: 33-60 %). Among healthy primiparous cows, 80% (95% CI: 75-84 %) were pregnant by 210 DIM, but pregnancy risk at that time was reduced for primiparous cows with TD-single (72%; 95% CI: 65-79 %) and TD+ (62%; 95% CI: 49-75 %). In healthy multiparous cows, the 120 DIM pregnancy risk was 53% (95% CI: 49-56 %), which was reduced for TD-single (36%; 95% CI: 31-42 %) and TD+ (30%; 95% CI: 24-38 %). The 210 DIM pregnancy risk for healthy multiparous cows was 70% (95% CI: 67-72 %), being higher than the 210 DIM pregnancy risk for multiparous cows with TD-single (47%; 95% CI: 42-53 %) or TD+ (46%; 95% CI: 38-54 %). Cows with TD-all presented similar pregnancy risk estimates as for TD+. Cox proportional hazards regressions provided similar magnitudes of effects as the Kaplan-Meier estimates. Survival analysis and decision tree models identified parity as the most influential variable affecting fertility. Both modeling techniques concurred that TD + had a greater effect than TD-single on the probability of pregnancy at 120 and 210 DIM. Decision trees for individual TD identified that displaced abomasum affected fertility at 120 DIM in primiparous while metritis was the most influential TD at 120 and 210 DIM for multiparous cows. The data were too sparse to assess multiple interactions in multivariable Cox proportional hazard models for individual TD. Machine learning helped to explore interactions between individual TD to study their hierarchical effect on fertility, identifying conditional relationships that merit further investigation

Effect of Dietary Phosphate Deprivation on Red Blood Cell Parameters of Periparturient Dairy Cows

Postparturient hemoglobinuria is a sporadic disease characterized by intravascular hemolysis and hemoglobinuria in early lactating dairy cows. The condition has empirically been associated with phosphorus (P) deficiency or hypophosphatemia; however, the exact etiology remains obscure. This paper summarizes two controlled studies investigating the effect of P deprivation during the transition period. In Study I, 36 late pregnant dairy cows were randomly assigned to either a diet with low, or adequate, P content from four weeks before calving to four weeks after calving. In Study II, 30 late pregnant dairy cows were again assigned to either a diet with low, or adequate, P for the last four weeks before calving only. Pronounced hypophosphatemia developed during periods of restricted P supply. In early lactation, a subtle decline of the red blood cell count occurred independently of the dietary P supply. In Study I, anemia developed in 11 cows on deficient P supply, which was associated with hemoglobinuria in five cases. Neither erythrocyte total P content nor osmotic resistance of erythrocytes were altered by dietary P deprivation. Restricted dietary P supply, particularly in early lactation, may lead to postparturient hemoglobinuria, but more frequently causes clinically inapparent hemolysis and anemia in cows

Between and within-herd variation in blood and milk biomarkers in Holstein cows in early lactation

Both blood- and milk-based biomarkers have been analysed for decades in research settings, although often only in one herd, and without focus on the variation in the biomarkers that are specifically related to herd or diet. Biomarkers can be used to detect physiological imbalance and disease risk and may have a role in precision livestock farming (PLF). For use in PLF, it is important to quantify normal variation in specific biomarkers and the source of this variation. The objective of this study was to estimate the between- and within-herd variation in a number of blood metabolites (β-hydroxybutyrate (BHB), non-esterified fatty acids, glucose and serum IGF-1), milk metabolites (free glucose, glucose-6-phosphate, urea, isocitrate, BHB and uric acid), milk enzymes (lactate dehydrogenase and N-acetyl-β-D-glucosaminidase (NAGase)) and composite indicators for metabolic imbalances (Physiological Imbalance-index and energy balance), to help facilitate their adoption within PLF. Blood and milk were sampled from 234 Holstein dairy cows from 6 experimental herds, each in a different European country, and offered a total of 10 different diets. Blood was sampled on 2 occasions at approximately 14 days-in-milk (DIM) and 35 DIM. Milk samples were collected twice weekly (in total 2750 samples) from DIM 1 to 50. Multilevel random regression models were used to estimate the variance components and to calculate the intraclass correlations (ICCs). The ICCs for the milk metabolites, when adjusted for parity and DIM at sampling, demonstrated that between 12% (glucose-6-phosphate) and 46% (urea) of the variation in the metabolites’ levels could be associated with the herd-diet combination. Intraclass Correlations related to the herd-diet combination were generally higher for blood metabolites, from 17% (cholesterol) to approximately 46% (BHB and urea). The high ICCs for urea suggest that this biomarker can be used for monitoring on herd level. The low variance within cow for NAGase indicates that few samples would be needed to describe the status and potentially a general reference value could be used. The low ICC for most of the biomarkers and larger within cow variation emphasises that multiple samples would be needed - most likely on the individual cows - for making the biomarkers useful for monitoring. The majority of biomarkers were influenced by parity and DIM which indicate that these should be accounted for if the biomarker should be used for monitoring

Potential of milk mid-infrared spectra to predict nitrogen use efficiency of individual dairy cows in early lactation

Improving nitrogen use efficiency (NUE) at both the individual cow and the herd level has become a key target in dairy production systems, for both environmental and economic reasons. Cost-effective and large-scale phenotyping methods are required to improve NUE through genetic selection and by feeding and management strategies. The aim of this study was to evaluate the possibility of using mid-infrared (MIR) spectra of milk to predict individual dairy cow NUE during early lactation. Data were collected from 129 Holstein cows, from calving until 50 d in milk, in 3 research herds (Denmark, Ireland, and the UK). In 2 of the herds, diets were designed to challenge cows metabolically, whereas a diet reflecting local management practices was offered in the third herd. Nitrogen intake (kg/d) and nitrogen excreted in milk (kg/d) were calculated daily. Nitrogen use efficiency was calculated as the ratio between nitrogen in milk and nitrogen intake, and expressed as a percentage. Individual daily values for NUE ranged from 9.7 to 81.7%, with an average of 36.9% and standard deviation of 10.4%. Milk MIR spectra were recorded twice weekly and were standardized into a common format to avoid bias between apparatus or sampling periods. Regression models predicting NUE using milk MIR spectra were developed on 1,034 observations using partial least squares or support vector machines regression methods. The models were then evaluated through (1) a cross-validation using 10 subsets, (2) a cow validation excluding 25% of the cows to be used as a validation set, and (3) a diet validation excluding each of the diets one by one to be used as validation sets. The best statistical performances were obtained when using the support vector machines method. Inclusion of milk yield and lactation number as predictors, in combination with the spectra, also improved the calibration. In cross-validation, the best model predicted NUE with a coefficient of determination of cross-validation of 0.74 and a relative error of 14%, which is suitable to discriminate between low- and high-NUE cows. When performing the cow validation, the relative error remained at 14%, and during the diet validation the relative error ranged from 12 to 34%. In the diet validation, the models showed a lack of robustness, demonstrating difficulties in predicting NUE for diets and for samples that were not represented in the calibration data set. Hence, a need exists to integrate more data in the models to cover a maximum of variability regarding breeds, diets, lactation stages, management practices, seasons, MIR instruments, and geographic regions. Although the model needs to be validated and improved for use in routine conditions, these preliminary results showed that it was possible to obtain information on NUE through milk MIR spectra. This could potentially allow large-scale predictions to aid both further genetic and genomic studies, and the development of farm management tools

Genome-wide association for metabolic clusters in early-lactation Holstein dairy cows

The aim of this study was to detect the genomic region or regions associated with metabolic clusters in early-lactation Holstein cows. This study was carried out in 2 experiments. In experiment I, which was carried out on 105 multiparous Holstein cows, animals were classified through k-means clustering on log-transformed and standardized concentrations of blood glucose, insulin-like growth factor I, free fatty acids, and β-hydroxybutyrate at 14 and 35 d in milk (DIM), into metabolic clusters, either balanced (BAL) or other (OTR). Forty percent of the animals were categorized in the BAL group, and the remainder were categorized as OTR. The cows were genotyped for a total of 777,962 SNP. A genome-wide association study was performed, using a case-control approach through the GEMMA software, accounting for population structure. We found 8 SNP (BTA11, BTA23, and BTAX) associated with the predicted metabolic clusters. In experiment II, carried out on 4,267 second-parity Holstein cows, milk samples collected starting from the first week until 50 DIM were used to determine Fourier-transform mid-infrared (FT-MIR) spectra and subsequently to classify the animals into the same metabolic clusters (BAL vs. OTR). Twenty-eight percent of the animals were categorized in the BAL group, and the remainder were classified in the OTR category. Although daily milk yield was lower in BAL cows, we found no difference in daily fat- and protein-corrected milk yield in cows from the BAL metabolic cluster compared with those in the OTR metabolic cluster. In the next step, a single-step genomic BLUP was used to identify the genomic region(s) associated with the predicted metabolic clusters. The results revealed that prediction of metabolic clusters is a highly polygenic trait regulated by many small-sized effects. The region of 36,258 to 36,295 kb on BTA27 was the highly associated region for the predicted metabolic clusters, with the closest genes to this region (ANK1 and miR-486) being related to hematopoiesis, erythropoiesis, and mammary gland development. The heritability for metabolic clustering was 0.17 (SD 0.03), indicating that the use of FT-MIR spectra in milk to predict metabolic clusters in early-lactation across a large number of cows has satisfactory potential to be included in genetic selection programs for modern dairy cows

Comparison of the transcriptome in circulating leukocytes in early lactation between primiparous and multiparous cows provides evidence for age-related changes

Background: Previous studies have identified many immune pathways which are consistently altered in humans and model organisms as they age. Dairy cows are often culled at quite young ages due to an inability to cope adequately with metabolic and infectious diseases, resulting in reduced milk production and infertility. Improved longevity is therefore a desirable trait which would benefit both farmers and their cows. This study analysed the transcriptome derived from RNA-seq data of leukocytes obtained from Holstein cows in early lactation with respect to lactation number. Results: Samples were divided into three lactation groups for analysis: i) primiparous (PP, n = 53), ii) multiparous in lactations 2–3 (MP 2–3, n = 121), and iii) MP in lactations 4–7 (MP > 3, n = 55). Leukocyte expression was compared between PP vs MP > 3 cows with MP 2–3 as background using DESeq2 followed by weighted gene co-expression network analysis (WGCNA). Seven modules were significantly correlated (r ≥ 0.25) to the trait lactation number. Genes from the modules which were more highly expressed in either the PP or MP > 3 cows were pooled, and the gene lists subjected to David functional annotation cluster analysis. The top three clusters from modules more highly expressed in the PP cows all involved regulation of gene transcription, particularly zinc fingers. Another cluster included genes encoding enzymes in the mitochondrial beta-oxidation pathway. Top clusters up-regulated in MP > 3 cows included the terms Glycolysis/Gluconeogenesis, C-type lectin, and Immunity. Differentially expressed candidate genes for ageing previously identified in the human blood transcriptome up-regulated in PP cows were mainly associated with T-cell function (CCR7, CD27, IL7R, CAMK4, CD28), mitochondrial ribosomal proteins (MRPS27, MRPS9, MRPS31), and DNA replication and repair (WRN). Those up-regulated in MP > 3 cows encoded immune defence proteins (LYZ, CTSZ, SREBF1, GRN, ANXA5, ADARB1). Conclusions: Genes and pathways associated with lactation number in cows were identified for the first time to date, and we found that many were comparable to those known to be associated with ageing in humans and model organisms. We also detected changes in energy utilization and immune responses in leukocytes from older cows

Genetic parameters for milk urea and its relationship with milk yield and compositions in Holstein dairy cows

The aim was to estimate genetic parameters for milk urea (MU) concentration and its relationship with milk yield and compositions in Holstein dairy Cows. Edited data were 90,594 test-day records of milk yield and composition collected during 2015 to 2018 on 13,737 lactations obtained from 7,850 Holstein cows in 50 herds. Random regression test-day model was used to estimate genetic parameters. (Co)variance components were estimated with the Bayesian Gibbs sampling method using a single chain of 400,000 iterates. The first 50,000 iterates of each chain were regarded as a burn-in period. Mean (SD) of MU was 23.03 (5.99) and 22.41 (5.74) mg/dl in primiparous and multiparous cows, respectively. Average heritability estimates for daily MU was 0.33 (SD = 0.02) ranged 0.29 to 0.36 and 0.32 (SD = 0.03) ranged 0.27 to 0.34, respectively, for primiparous and multiparous cows. The mean (SD) genetic correlation between MU and milk yield, fat yield, protein yield, lactose yield, fat percentage, protein percentage, lactose percentage, and somatic cell score was, respectively, -0.02 (0.03), -0.02 (0.01), 0.01 (0.04), 0.01 (0.03), 0.00 (0.07), -0.03 (0.04), 0.00 (0.01), -0.11 (0.06) in primiparous cows. The corresponding values in multiparous cows were -0.01 (0.02), -0.01 (0.03), -0.04 (0.04), -0.04 (0.04), 0.04 (0.04), 0.04 (0.07), -0.03 (0.09), 0.06 (0.11), respectively. The results indicate that selection on MU is possible with no effect on milk yield or compositions, however, relationships between MU and other important traits such as longevity, metabolic diseases, and fertility are needed

Validation of Dairy Cow Bodyweight Prediction Using Traits Easily Recorded by Dairy Herd Improvement Organizations and Its Potential Improvement Using Feature Selection Algorithms

Knowing the body weight (BW) of a cow at a specific moment or measuring its changes through time is of interest for management purposes. The current work aimed to validate the feasibility of predicting BW using the day in milk, parity, milk yield, and milk mid-infrared (MIR) spectrum from a multiple-country dataset and reduce the number of predictors to limit the risk of over-fitting and potentially improve its accuracy. The BW modeling procedure involved feature selections and herd-independent validation in identifying the most interesting subsets of predictors and then external validation of the models. From 1849 records collected in 9 herds from 360 Holstein cows, the best performing models achieved a root mean square error (RMSE) for the herd-independent validation between 52 ± 2.34 kg to 56 ± 3.16 kg, including from 5 to 62 predictors. Among these models, three performed remarkably well in external validation using an independent dataset (N = 4067), resulting in RMSE ranging from 52 to 56 kg. The results suggest that multiple optimal BW predictive models coexist due to the high correlations between adjacent spectral points

Estimation of genetic parameters for predicted nitrogen use efficiency and losses in early lactation of Holstein cows

The objective of this study was to estimate genetic parameters of predicted N use efficiency (PNUE) and N losses (PNL) as proxies of N use and loss for Holstein cows. Furthermore, we have assessed approximate genetic correlations between PNUE, PNL, and dairy production, health, longevity, and conformation traits. These traits are considered important in many countries and are currently evaluated by the International Bull Evaluation Service (Interbull). The values of PNUE and PNL were obtained by using the combined milk mid-infrared (MIR) spectrum, parity, and milk yield–based prediction equations on test-day MIR records with days in milk (DIM) between 5 and 50 d. After editing, the final data set comprised 46,163 records of 21,462 cows from 154 farms in 5 countries. Each trait was divided into primiparous and multiparous (including second to fifth parity) groups. Genetic parameters and breeding values were estimated by using a multitrait (2-trait, 2-parity classes) repeatability model. Herd-year-season of calving, DIM, age of calving, and parity were used as fixed effects. Random effects were defined as parity (within-parity permanent environment), nongenetic cow (across-parity permanent environment), additive genetic animal, and residual effects. The estimated heritability of PNUE and PNL in the first and later parity were 0.13, 0.12, 0.14, and 0.13, and the repeatability values were 0.49, 0.40, 0.55, and 0.43, respectively. The estimated approximate genetic correlations between PNUE and PNL were negative and high (from −0.89 to −0.53), whereas the phenotypic correlations were also negative but relatively low (from −0.45 to −0.11). At a level of reliability of more than 0.30 for all novel traits, a total of 504 bulls born after 1995 had also publishable Interbull multiple-trait across-country estimated breeding values (EBV). The approximate genetic correlations between PNUE and the other 30 traits of interest, estimated as corrected correlations between EBV of bulls, ranged from −0.46 (udder depth) to 0.47 (milk yield). Obtained results showed the complex genetic relationship between efficiency, production, and other traits: for instance, more efficient cows seem to give more milk, which is linked to deeper udders, but seem to have lower health, fertility, and longevity. Additionally, the approximate genetic correlations between PNL, lower values representing less loss of N, and the 30 other traits, were from −0.32 (angularity) to 0.57 (direct calving ease). Even if further research is needed, our results provided preliminary evidence that the PNUE and PNL traits used as proxies could be included in genetic improvement programs in Holstein cows and could help their management

Genome-wide association for milk production and lactation curve parameters in Holstein dairy cows

The aim of this study was to identify genomic regions associated with 305-day milk yield and lactation curve parameters on primiparous (n = 9,910) and multiparous (n = 11,158) Holstein cows. The SNP solutions were estimated using a weighted single-step genomic BLUP approach and imputed high-density panel (777k) genotypes. The proportion of genetic variance explained by windows of 50 consecutive SNP (with an average of 165 Kb) was calculated, and regions that accounted for more than 0.50% of the variance were used to search for candidate genes. Estimated heritabilities were 0.37, 0.34, 0.17, 0.12, 0.30 and 0.19, respectively, for 305-day milk yield, peak yield, peak time, ramp, scale and decay for primiparous cows. Genetic correlations of 305-day milk yield with peak yield, peak time, ramp, scale and decay in primiparous cows were 0.99, 0.63, 0.20, 0.97 and −0.52, respectively. The results identified three windows on BTA14 associated with 305-day milk yield and the parameters of lactation curve in primi- and multiparous cows. Previously proposed candidate genes for milk yield supported by this work include GRINA, CYHR1, FOXH1, TONSL, PPP1R16A, ARHGAP39, MAF1, OPLAH and MROH1, whereas newly identified candidate genes are MIR2308, ZNF7, ZNF34, SLURP1, MAFA and KIFC2 (BTA14). The protein lipidation biological process term, which plays a key role in controlling protein localization and function, was identified as the most important term enriched by the identified genes