Genetic parameters of differential somatic cell count, milk composition, and cheese-making traits measured and predicted using spectral data in Holstein cows

Evidence of additive genetic variation for major milk proteins in dairy cows: A meta-analysis

In the past, there have been reports of genetic parameters for milk proteins in various dairy cattle populations. The high variability among genetic parameter estimates has been caused by this. This study aimed to use a random-effects meta-analysis model to compile published estimates of genetic parameter for major milk proteins of α-lactalbumin, β-lactoglobulin, sum of whey proteins, casein, αs1-casein, αs2-casein, β-casein, and κ-casein in dairy cows. The study used a total of 140 heritability and 256 genetic correlation estimates from 23 papers published between 2004 and 2022. The estimated range of milk protein heritability is from 0.284 (for α-lactalbumin in milk) to 0.596 (for sum of whey proteins). The genetic correlation estimates between casein and milk yield, milk fat and protein percentages were -0.461, 0.693, and 0.976, respectively (p < 0.05). The genetic correlation estimates between milk proteins expressed as a percentage of milk were significant and varied from 0.177 (between β-lactoglobulin and κ-casein) to 0.892 (between αs1-casein and αs2-casein). Moderate-to-high heritability estimates for milk proteins and their low genetic associations with milk yield and composition indicated the possibility for improving milk proteins in a genetic selection plan with negligible correlated effects on production traits in dairy cows.

Genetic parameters for novel mastitis traits defined by combining test-day somatic cell score and differential somatic cell count in the first lactation of Japanese Holsteins

In this study, we aimed to improve current udder health genetic evaluations by addressing the limitations of monthly sampled somatic cell score (SCS) for distinguishing cows with robust innate immunity from those susceptible to chronic infections. The objectives were to (1) establish novel somatic cell traits by integrating SCS and the differential somatic cell count (DSCC), which represents the combined proportion of polymorphonuclear leukocytes and lymphocytes in somatic cells and (2) estimate genetic parameters for the new traits, including their daily heritability and genetic correlations with milk production traits and SCS, using a random regression test-day model (RRTDM). We derived 3 traits, termed ML_SCS_DSCC, SCS_4_DSCC_65_binary, and ML_SCS_DSCC_binary, by using milk loss (ML) estimates at corresponding SCS and DSCC levels, thresholds established in previous studies, and a threshold established from milk loss estimates, respectively. Data consisted of test-day records collected during January 2021 through March 2022 from 265 herds in Hokkaido, Japan. From these records, we extracted records between 7 to 305 d in milk (DIM) in the first lactation to fit the RRTDM. The model included the random effect of herd-test-day, the fixed effect of year-month, fixed lactation curves nested with calving age groups, and random regressions with Legendre polynomials of order 3 for additive genetic and permanent environmental effects. The analysis was performed using Gibbs sampling with Gibbsf90+ software. The averages (ranges) of the daily heritability estimates over lactation were 0.086 (0.075-0.095) for SCS, 0.104 (0.073-0.127) for ML_SCS_DSCC, 0.137 (0.014-0.297) for SCS_4_DSCC_65_binary, and 0.138 (0.115-0.185) for ML_SCS_DSCC_binary; the heritability curve for SCS_4_DSCC_65_binary was erratic. Genetic correlations within the trait decreased as the DIM interval widened, especially for those integrating DSCC, indicating that these traits should be analyzed using RRTDM rather than repeatability models. The averages (ranges) of genetic correlations with milk yield over lactation were 0.01 (-0.22 to 0.28) for SCS, -0.05 (-0.40 to 0.13) for ML_SCS_DSCC, -0.08 (-0.17 to 0.09) for SCS_4_DSCC_65_binary, and -0.08 (-0.22 to 0.27) for ML_SCS_DSCC_binary. Compared with SCS, the newly defined traits exhibited slightly stronger negative genetic correlations with milk yield. Especially in late lactation stages, the genetic correlation between ML_SCS_DSCC and milk yield was significantly below zero, with a posterior median of -0.40. Furthermore, the new traits showed positive correlations with SCS, having estimates varying from 0.68 to 0.85 for ML_SCS_DSCC, 0.14 to 0.47 for SCS_4_DSCC_65_binary, and 0.61 to 0.66 for ML_SCS_DSCC_binary, depending on DIM. Considering that ML_SCS_DSCC and ML_SCS_DSCC_binary have relatively high heritability (compared with SCS) and favorable genetic correlations with milk production traits and SCS, their incorporation into breeding programs appears promising. Nevertheless, their genetic relationships with (sub)clinical mastitis require further investigation.

Genomic prediction of blood biomarkers of metabolic disorders in Holstein cattle using parametric and nonparametric models

BACKGROUND

Metabolic disturbances adversely impact productive and reproductive performance of dairy cattle due to changes in endocrine status and immune function, which increase the risk of disease. This may occur in the post-partum phase, but also throughout lactation, with sub-clinical symptoms. Recently, increased attention has been directed towards improved health and resilience in dairy cattle, and genomic selection (GS) could be a helpful tool for selecting animals that are more resilient to metabolic disturbances throughout lactation. Hence, we evaluated the genomic prediction of serum biomarkers levels for metabolic distress in 1353 Holsteins genotyped with the 100K single nucleotide polymorphism (SNP) chip assay. The GS was evaluated using parametric models best linear unbiased prediction (GBLUP), Bayesian B (BayesB), elastic net (ENET), and nonparametric models, gradient boosting machine (GBM) and stacking ensemble (Stack), which combines ENET and GBM approaches.

RESULTS

The results show that the Stack approach outperformed other methods with a relative difference (RD), calculated as an increment in prediction accuracy, of approximately 18.0% compared to GBLUP, 12.6% compared to BayesB, 8.7% compared to ENET, and 4.4% compared to GBM. The highest RD in prediction accuracy between other models with respect to GBLUP was observed for haptoglobin (hapto) from 17.7% for BayesB to 41.2% for Stack; for Zn from 9.8% (BayesB) to 29.3% (Stack); for ceruloplasmin (CuCp) from 9.3% (BayesB) to 27.9% (Stack); for ferric reducing antioxidant power (FRAP) from 8.0% (BayesB) to 40.0% (Stack); and for total protein (PROTt) from 5.7% (BayesB) to 22.9% (Stack). Using a subset of top SNPs (1.5k) selected from the GBM approach improved the accuracy for GBLUP from 1.8 to 76.5%. However, for the other models reductions in prediction accuracy of 4.8% for ENET (average of 10 traits), 5.9% for GBM (average of 21 traits), and 6.6% for Stack (average of 16 traits) were observed.

CONCLUSIONS

Our results indicate that the Stack approach was more accurate in predicting metabolic disturbances than GBLUP, BayesB, ENET, and GBM and seemed to be competitive for predicting complex phenotypes with various degrees of mode of inheritance, i.e. additive and non-additive effects. Selecting markers based on GBM improved accuracy of GBLUP.

Observational study on the associations between milk yield, composition, and coagulation properties with blood biomarkers of health in Holstein cows

The considerable increase in the production capacity of individual cows owing to both selective breeding and innovations in the dairy sector has posed challenges to management practices in terms of maintaining the nutritional and metabolic health status of dairy cows. In this observational study, we investigated the associations between milk yield, composition, and technological traits and a set of 21 blood biomarkers related to energy metabolism, liver function or hepatic damage, oxidative stress, and inflammation or innate immunity in a population of 1,369 high-yielding Holstein-Friesian dairy cows. The milk traits investigated in this study included 4 production traits (milk yield, fat yield, protein yield, daily milk energy output), 5 traits related to milk composition (fat, protein, casein, and lactose percentages and urea), 11 milk technological traits (5 milk coagulation properties and 6 curd-firming traits). All milk traits (i.e., production, composition, and technological traits) were analyzed according to a linear mixed model that included the days in milk, the parity order, and the blood metabolites (tested one at a time) as fixed effects and the herd and date of sampling as random effects. Our findings revealed that milk yield and daily milk energy output were positively and linearly associated with total cholesterol, nonesterified fatty acids, urea, aspartate aminotransferase, γ-glutamyl transferase, total bilirubin, albumin, and ferric-reducing antioxidant power, whereas they were negatively associated with glucose, creatinine, alkaline phosphatase, total reactive oxygen metabolites, and proinflammatory proteins (ceruloplasmin, haptoglobin, and myeloperoxidase). Regarding composition traits, the protein percentage was negatively associated with nonesterified fatty acids and β-hydroxybutyrate (BHB), while the fat percentage was positively associated with BHB, and negatively associated with paraoxonase. Moreover, we found that the lactose percentage increased with increasing cholesterol and albumin and decreased with increasing ceruloplasmin, haptoglobin, and myeloperoxidase. Milk urea increased with an increase in cholesterol, blood urea, nonesterified fatty acids, and BHB, and decreased with an increase in proinflammatory proteins. Finally, no association was found between the blood metabolites and milk coagulation properties and curd-firming traits. In conclusion, this study showed that variations in blood metabolites had strong associations with milk productivity traits, the lactose percentage, and milk urea, but no relationships with technological traits of milk. Specifically, increasing levels of proinflammatory and oxidative stress metabolites, such as ceruloplasmin, haptoglobin, myeloperoxidase, and total reactive oxygen metabolites, were shown to be associated with reductions in milk yield, daily milk energy output, lactose percentage, and milk urea. These results highlight the close connection between the metabolic and innate immunity status and production performance. This connection is not limited to specific clinical diseases or to the transition phase but manifests throughout the entire lactation. These outcomes emphasize the importance of identifying cows with subacute inflammatory and oxidative stress as a means of reducing metabolic impairments and avoiding milk fluctuations.

Uncovering genetic parameters and environmental influences on fertility, milk production, and quality in autochthonous Reggiana cattle

Reggiana is a local cattle breed from northern Italy known for its rusticity and profitability, due to the production of branded Parmigiano Reggiano cheese. To ensure the persistence of such profitability in the long term, an adequate breeding program is required. To this aim, in the present study we estimate the genetic parameters of the main productive and reproductive traits, and we evaluate the effect of genotype by environment interaction (GxE) on these traits using 2 environmental covariates: (1) productivity and (2) temperature-humidity index (THI). Milk, fat, protein, and casein yield were considered as daily production traits, whereas protein, fat, casein percentage, casein index, and somatic cell score were considered as milk quality traits. Finally, reproductive traits such as the number of inseminations, days open, calving interval, and calving-to-first-insemination interval were evaluated. Reggiana cattle produce an average of 19 kg of milk per day with 3.7% fat and 3.4% protein content and have excellent fertility parameters. Compared with other breeds, they have slightly lower heritability for production and quality for production traits (e.g., 0.12 [0.09; 0.15] for milk yield), but similar heritability for fertility traits. Milk, protein, and fat daily yields are highly correlated but negatively correlated with the percentage of protein, fat, and casein, whereas fertility traits have an unfavorable genetic correlation with daily production traits. When considering productivity, a consistent amount of variability due to GxE was observed for all daily production traits, somatic cell count, and casein index. A modest amount of GxE was observed for fertility parameters, while the percentage of solid content showed almost no GxE effect. A similar situation occurred when considering the THI, but no GxE interaction was observed for reproduction traits. In conclusion, this study provides useful information for the implementation of accurate selection plans in this local breed, accounting for environmental plasticity measured through the consistent GxE interaction observed.

Novel insights into the associations between immune cell population distribution in mammary glands and milk minerals in Holstein cows

Udder health has a crucial role in sustainable milk production, and various reports have pointed out that changes in udder condition seem to affect milk mineral content. The somatic cell count (SCC) is the most recognized indicator for the determination of udder health status. Recently, a new parameter, the differential somatic cell count (DSCC), has been proposed for a more detailed evaluation of intramammary infection patterns. Specifically, the DSCC is the combined proportions of polymorphonuclear neutrophils and lymphocytes (PMN-LYM) on the total SCC, with macrophages (MAC) representing the remainder proportion. In this study, we evaluated the association between DSCC in combination with SCC on a detailed milk mineral profile in 1,013 Holstein-Friesian cows reared in 5 herds. An inductively coupled plasma-optical emission spectrometry was used to quantify 32 milk mineral elements. Two different linear mixed models were fitted to explore the associations between the milk mineral elements and first, the DSCC combined with SCC, and second, DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the proportion of PMN-LYM and MAC by SCC. We observed a significant positive association between SCC and milk Na, S, and Fe levels. Differential somatic cell count showed an opposite behavior to the one displayed by SCC, with a negative association with Na and positive association with K milk concentrations. When considering DSCC as count, Na and K showed contrasting behavior when associated with PMN-LYM or MAC counts, with decreasing of Na content and increasing K when associated with increasing PMN-LYM counts, and increasing Na and decreasing K when associated with increasing MAC count. These findings confirmed that an increase in SCC is associated with altered milk Na and K amounts. Moreover, MAC count seemed to mirror SCC patterns, with the worsening of inflammation. Differently, PMN-LYM count exhibited patterns of associations with milk Na and K contents attributable more to LYM than PMN, given the nonpathological condition of the majority of the investigated population. An interesting association was observed for milk S content, which increased with increasing of inflammatory conditions (i.e., increased SCC and MAC count) probably attributable to its relationship with milk proteins, especially whey proteins. Moreover, milk Fe content showed positive associations with the PMN-LYM population, highlighting its role in immune regulation during inflammation. Further studies including individuals with clinical condition are needed to achieve a comprehensive view of milk mineral behavior during udder health impairment.

Heritability and genetic correlations of total and differential somatic cell count with milk yield and composition traits in Italian Simmental cows

Costs of production have deeply increased each year in the last decades, breeders are continuously looking for more cost effective and more efficient ways to produce milk. Despite the major signs of progress in productivity, it is fundamental to optimize rather than maximize the performances of the dairy cows. Mastitis is still a highly prevalent disease in the dairy sector which causes several economic losses and environmental effect. Its accurate and early diagnosis is crucial to improve profitability of dairy cows and contribute to a more sustainable dairy industry. Among mastitis reduction strategies, there is the urgent need to implement breeding objectives to select cows displaying mastitis resistance by investigating the genetic mechanisms at the base of the inflammatory response. Therefore, in this study we aimed to further understand the genetic background of the differential somatic cell count (DSCC), which provides thorough insights on the actual inflammatory status of the mammary glands. The objectives of this study were to estimate on a cohort of 20,215 Italian Simmental cows over a 3-yr period: (1) the heritability and repeatability values of somatic cell score (SCS) and DSCC, (2) the genetic and phenotypic correlations between these 2 traits and milk production and milk composition traits, (3) the heritability and repeatability values of SCS and DSCC within class of udder health status. Heritability was low both for SCS (0.06) and DSCC (0.08), whereas the repeatability values for these traits were 0.43 and 0.36, suggesting that the magnitude of cow permanent environmental effect for these traits is remarkable. The genetic and phenotypic correlation of SCS with DSCC was 0.612 and 0.605, respectively. Because both significantly differed from the unit, we must consider those traits as different ones. This latter aspect corroborates the need to consider the DSCC as a further indicator of inflammatory status which might be implemented in the Simmental breed genetic evaluation. It is worthy to mention that heritability estimates for SCS and DSCC were the highest in healthy cows compared with the other udder health classes. This implies that when the udder health status changes, it is most likely due to environmental factors rather than aspects related to the animal's genetics. In contrast, the highest additive genetic variance and heritability found for SCS and DSCC in the healthy group might reveal the potential to further implement breeding strategies to select for healthier animals.

High-concentrate diet elevates histone lactylation mediated by p300/CBP through the upregulation of lactic acid and induces an inflammatory response in mammary gland of dairy cows

High-concentrate diet can cause metabolic diseases, such as subacute ruminal acidosis (SARA), and secondary mastitis. To investigate the effect of SARA induced by high-concentrate diet on the lysine lactylation (Kla) and inflammatory responses in the mammary gland of dairy cows and the mechanism between them, we selected twelve mid-lactation Holstein cows with similar body conditions for modelling. They were randomly divided into two groups, fed a low-concentrate diet (LC) and a high-concentrate diet (HC) for 21 days. Our results showed that high-concentrate diet feeding significantly reduced ruminal pH, and the pH was below 5.6 for more than 3 h per day, indicating successful induction of the SARA model. Lactic acid concentrations in mammary gland and plasma were higher in the HC group than that in the LC group. HC diet feeding significantly up-regulated the expression levels of the Pan Kla, H3K18la, p300/CBP and monocarboxylate transporter 1 (MCT1) in the mammary gland. In addition, the mRNA expression levels of inflammatory factors were significantly regulated, including IL-1β, IL-1α, IL-6, IL-8, SAA3, and TNF-α, while the anti-inflammatory factor IL-10 was down-regulated. The mammary gland of HC group was structurally disorganized with incomplete glandular vesicles, with a large number of detached mammary epithelial cells and inflammatory cells infiltration. The up-regulation of TLR4, TNF-α, p-p65, and p-IκBα indicated that the TLR4/NF-κB signaling pathway was activated. In conclusion, this study found that HC diet feeding can induce SARA and increase the concentration of lactic acid in mammary gland and plasma. Then, lactic acid could be transported into cells by MCT1 and up-regulate the expression level of histone lactylation mediated by p300/CBP, and subsequently promote the activation of TLR4/NF-κB signaling pathway, ultimately causing inflammatory responses in the mammary gland.

Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by β-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CV_r) and herd-out (CV_h). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R² value from 5% to 75% for CV_r and 2% to 139% for CV_h, whereas the stacking ensemble increased the R² value from 4% to 70% for CV_r and 4% to 150% for CV_h. Considering the best model, with the CV_r scenario, good prediction accuracies were obtained for glucose (R² = 0.81), urea (R² = 0.73), albumin (R² = 0.75), total reactive oxygen metabolites (R² = 0.79), total thiol groups (R² = 0.76), ceruloplasmin (R² = 0.74), total proteins (R² = 0.81), globulins (R² = 0.87), and Na (R² = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

Integrating on-farm and genomic information improves the predictive ability of milk infrared prediction of blood indicators of metabolic disorders in dairy cows

BACKGROUND

Blood metabolic profiles can be used to assess metabolic disorders and to evaluate the health status of dairy cows. Given that these analyses are time-consuming, expensive, and stressful for the cows, there has been increased interest in Fourier transform infrared (FTIR) spectroscopy of milk samples as a rapid, cost-effective alternative for predicting metabolic disturbances. The integration of FTIR data with other layers of information such as genomic and on-farm data (days in milk (DIM) and parity) has been proposed to further enhance the predictive ability of statistical methods. Here, we developed a phenotype prediction approach for a panel of blood metabolites based on a combination of milk FTIR data, on-farm data, and genomic information recorded on 1150 Holstein cows, using BayesB and gradient boosting machine (GBM) models, with tenfold, batch-out and herd-out cross-validation (CV) scenarios.

RESULTS

The predictive ability of these approaches was measured by the coefficient of determination (R²). The results show that, compared to the model that includes only FTIR data, integration of both on-farm (DIM and parity) and genomic information with FTIR data improves the R² for blood metabolites across the three CV scenarios, especially with the herd-out CV: R² values ranged from 5.9 to 17.8% for BayesB, from 8.2 to 16.9% for GBM with the tenfold random CV, from 3.8 to 13.5% for BayesB and from 8.6 to 17.5% for GBM with the batch-out CV, and from 8.4 to 23.0% for BayesB and from 8.1 to 23.8% for GBM with the herd-out CV. Overall, with the model that includes the three sources of data, GBM was more accurate than BayesB with accuracies across the CV scenarios increasing by 7.1% for energy-related metabolites, 10.7% for liver function/hepatic damage, 9.6% for oxidative stress, 6.1% for inflammation/innate immunity, and 11.4% for mineral indicators.

CONCLUSIONS

Our results show that, compared to using only milk FTIR data, a model integrating milk FTIR spectra with on-farm and genomic information improves the prediction of blood metabolic traits in Holstein cattle and that GBM is more accurate in predicting blood metabolites than BayesB, especially for the batch-out CV and herd-out CV scenarios.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [α_S1-casein (CN), α_S2-CN, β-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (β-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (r_a) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Impact of somatic cell count combined with differential somatic cell count on milk protein fractions in Holstein cattle

Udder health in dairy herds is a very important issue given its implications for animal welfare and the production of high-quality milk. Somatic cell count (SCC) is the most widely used means of assessing udder health status. However, differential somatic cell count (DSCC) has recently been proposed as a new and more effective means of evaluating intramammary infection dynamics. Differential SCC represents the combined percentage of polymorphonuclear neutrophils and lymphocytes (PMN-LYM) in the total SCC, with macrophages (MAC) accounting for the remaining proportion. The aim of this study was to evaluate the association between SCC and DSCC and the detailed milk protein profile in a population of 1,482 Holstein cows. A validated reversed-phase HPLC method was used to quantify 4 caseins (CN), namely α_S1-CN, α_S2-CN, κ-CN, and β-CN, and 3 whey protein fractions, namely β-lactoglobulin, α-lactalbumin, and lactoferrin, which were expressed both quantitatively (g/L) and qualitatively (as a percentage of the total milk nitrogen content, %N). A linear mixed model was fitted to explore the associations between somatic cell score (SCS) combined with DSCC and the protein fractions expressed quantitatively and qualitatively. We ran an additional model that included DSCC expressed as PMN-LYM and MAC counts, obtained by multiplying the percentages of PMN-LYM and MAC by SCC for each cow in the data set. When the protein fractions were expressed as grams per liter, SCS was significantly negatively associated with almost all the casein fractions and positively associated with the whey protein α-lactalbumin, while DSCC was significantly associated with α_S1-CN, β-CN, and α-lactalbumin, but in the opposite direction to SCS. We observed the same pattern with the qualitative data (i.e., %N), confirming opposite effects of SCS and DSCC on milk protein fractions. The PMN-LYM count was only slightly associated with the traits of concern, although the pattern observed was the same as when both SCS and DSCC were included in the model. The MAC count, however, generally had a greater impact on many casein fractions, in particular decreasing both β-CN content (g/L) and proportion (%N), and exhibited the opposite pattern to the PMN-LYM count. Our results show that information obtained from both SCS and DSCC may be useful in assessing milk quality and protein fractions. They also demonstrate the potential of MAC count as a novel udder health trait.

Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak

Simple Summary

Increasing milk fat rate has a good effect on the milk quality of cattleyak. SNPs can help us find potential molecular markers for the milk fat traits of cattleyak, and they can be screened according to molecular markers when they are young. It provides a reference for cultivating high milk fat cattle population in the future. The results of this study suggest that the SORBS1 gene polymorphism is closely related to the milk fat traits of cattleyak, which could be used as a candidate genetic marker for milk fat trait selection in cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak. It has a certain reference value for the research and improvement of milk quality.

Abstract

This study aimed to find the SNPs in the SORBS1 gene of cattleyak, analyze the relationship between its polymorphisms and the milk fat traits, and find potential molecular markers for the milk fat traits of cattleyak. The polymorphism of the SORBS1 gene in 350 cattleyak from Hongyuan County (Sichuan, China) were detected by PCR and DNA sequencing, and the correlation between these SNPs and the milk production traits of cattleyak was analyzed. The results showed that there were nine SNPs in the CDS and their adjacent non-coding regions of the SORBS1 gene, and all SNPs have three genotypes. The correlation analysis found that the genotypes with superior milk fat traits in the other eight alleles were homozygous genotypes with a high genotype frequency except the g.96284 G > A (c.3090 G > A) (p < 0.05). However, at locus g.96284 G > A, the milk fat percentage, monounsaturated fatty acids (MUFAs), polyunsaturated fatty acids (PUFAs) and saturated fatty acids (SFAs) of the GA genotype were significantly higher than that of GG and AA genotypes (p < 0.05). Among these SNPs, three SNPs (g.6256 C > T (c.298 C > T), g.24791 A > G (c.706 A > G) and g.29121 A > G (c.979 A > G)) caused the amino acids change. The genotypes of the three SNPs consist of three haplotypes and four diplotypes. The amino acid mutation degree of diplotype H1–H1 (CCAAAA) was the highest, and its milk fat percentage, MUFAs, PUFAs and SFAs were also the highest (p < 0.05). Taken together, we found nine SNPs in the SORBS1 gene that are closely related to the milk fat traits of cattleyak. Moreover, the mutation of amino acids caused by SNPs had positive effects on the milk fat traits of cattleyak. H1-H1 is the dominant diplotype which significantly related to the milk fat traits of cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak.