The rates of mortality and involuntary culling of dairy calves and replacement heifers have great economic implications on the dairy cattle industry around the world. The main objectives of this study were: (1) to obtain population parameters of mortality and involuntary culling rates of dairy calves and replacement heifers; and, (2) to investigate the factors affecting mortality and involuntary culling rates in Chinese Holstein cattle. Two datasets containing records of birth, calving, and culling events from 142,833 Holstein cattle born between 1991 and 2018 were used in this study. The population parameters were obtained using dataset 1, which consisted of dairy calves and replacement heifers that died or were involuntarily culled. Three survival traits were defined in dataset 2, which consisted of females born from 1999 to 2018. A binomial logistic regression was used to analyze the risk factors on the survival traits. The mortality rate of dairy calves and replacement heifers from day 3 to 60, 61 to 365, and 366 to first calving was 5.5%, 7.4%, and 8.7%, and an unfavorable increasing trend was observed. Health events associated with digestive and respiratory or circulatory systems were the main death reasons. Herd-birth year, birth season, and dam parity had significant effects on survival traits. The results from this study will help farmers to better manage calves and replacement heifers and highlight the need to include survival traits in dairy calves and replacement heifers as part of national genetic evaluation schemes.

Genetic selection for improved climatic resilience is paramount to increase the long-term sustainability of high-producing dairy cattle, especially in face of climate change. Various physiological indicators, such as rectal temperature (RT), respiration rate score (RR), and drooling score (DS), can be used to genetically identify animals with more effective coping mechanisms in response to heat stress events. In this study, we investigated genetic parameters for RT, RR (score from 1-3), and DS (score from 1-3). Furthermore, we assessed the genetic relationship among these indicators and other economically important traits for the dairy cattle industry. After data editing, 59,265 (RT), 30,290 (RR), and 30,421 (DS) records from 13,592 lactating Holstein cows were used for the analyses. Variance components were estimated based on a multiple-trait repeatability animal model. The heritability ± standard error estimate for RT, RR, and DS was 0.06 ± 0.01, 0.04 ± 0.01, and 0.02 ± 0.01, respectively, whereas their repeatability was 0.19, 0.14, and 0.14, respectively. Moderate genetic correlations of RR with RT and DS (0.26 ± 0.11 and 0.25 ± 0.16) and nonsignificant correlation between RT and DS (-0.11 ± 0.14) were observed. Furthermore, the approximate genetic correlations between RT, RR, and DS with 12 production, 29 conformation, 5 fertility and reproduction, 5 health, and 9 longevity-indicator traits were assessed. In general, the approximate genetic correlations calculated were low to moderate. In summary, 3 physiological indicators of heat stress response were measured in a large number of animals and shown to be lowly heritable. There is a value in developing a selection index including all the 3 indicators to improve heat tolerance in dairy cattle. All the unfavorable genetic relationships observed between heat tolerance and other economically important traits can be accounted for in a selection index to enable improved climatic resilience while also maintaining or increasing productivity in Holstein cattle.

Inflammaging is well understood in the study of humans; however it is rarely reported for dairy cows. To understand the changing pattern of the gut microbiota, inflammatory status and milk production performance during the aging process in cows, we grouped 180 cows according to their lactation period: L1 (n = 60, 1st lactation), L3 (n = 60, 3rd lactation), and L5+ (n = 60, at least 5th lactation) and analyzed their milk components and daily milk yields to evaluate the changing pattern of milk production. The microbiota was analyzed using high-throughput sequencing of amplicons of 16S rRNA, which also allowed us to predict the functions of microbes and then study the changing pattern of the ruminal and fecal microbiota. Serum cytokines were measured to study the progress of inflammaging in the cows. We found that old cows (L5+) suffered from a long-term and low-level chronic inflammation, as indicated by significantly higher levels of inflammatory cytokines IL-10, TNF-α, and TGF-β in the L5+ group (p 2), showed a similar change pattern during the aging process, both in the rumen and in feces, and across the six farms. Beneficial bacteria, like Bacteroidaceae, Eubacterium, and Bifidobacterium, displayed lower abundance in the feces of the L5+ group (LEfSe, LDA>2). Reconstruction of the fecal bacteria community indicated transformation of the fermenting pattern of older cows’ (L5+) feces microbiota, with increased functions related the protein metabolism and fewer functions related to carbohydrate and lipid metabolism compared with those in L1 (p < 0.05). Finally, the connections among these changing patterns were revealed using redundancy analysis and network analysis. The results support the hypothesis of prolonging a cows’ productive life and improve dairy cow milk productive performances by manipulating the gut microbiota.

Genetic selection for resilience is essential to improve the long-term sustainability of the dairy cattle industry, especially the ability of cows to maintain their level of production when exposed to environmental disturbances. Recording of daily milk yield provides an opportunity to develop resilience indicators based on milk losses and fluctuations in daily milk yield caused by environmental disturbances. In this context, our study aimed to explore milk loss traits and measures of variability in daily milk yield, including log-transformed standard deviation of milk deviations (Lnsd), lag-1 autocorrelation (Ra), and skewness of the deviations (Ske), as indicators of general resilience in dairy cows. The unperturbed dynamics of milk yield as well as milk loss were predicted using an iterative procedure of lactation curve modeling. Milk fluctuations were defined as a period of at least 10 successive days of negative deviations in which milk yield dropped at least once below 90% of the expected values. Genetic parameters of these indicators and their genetic correlation with economically important traits were estimated using single-trait and bivariate animal models and 8,935 lactations (after quality control) from 6,816 Chinese Holstein cows. In general, cows experienced an average of 3.73 environmental disturbances with a milk loss of 267 kg of milk per lactation. Each fluctuation lasted for 19.80 ± 11.46 days. Milk loss traits are heritable with heritability estimates ranging from 0.004 to 0.061. The heritabilities differed between Lnsd (0.135–0.250), Ra (0.008–0.058), and Ske (0.001–0.075), with the highest heritability estimate of 0.250 ± 0.020 for Lnsd when removing the first and last 10 days in milk in a lactation (Lnsd2). Based on moderate to high genetic correlations, lower Lnsd2 is associated with less milk losses, better reproductive performance, and lower disease incidence. These findings indicate that among the variables evaluated, Lnsd2 is the most promising indicator for breeding for improved resilience in Holstein cattle.

Due to the complexity of longevity trait in dairy cattle, two groups of trait definitions are widely used to measure longevity, either covering the full lifespan or representing only a part of it to achieve an early selection. Usually, only one group of longevity definition is used in breeding program for one population, and genetic studies on the comparisons of two groups of trait definitions are scarce. Based on the data of eight traits well representing the both groups of trait definitions, the current study investigated genetic parameters and genetic architectures of longevity in Holsteins. Heritabilities and correlations of eight longevity traits were estimated using single-trait and multi-trait animal models, with the data from 103,479 cows. Among the cows with phenotypes, 2,630 cows were genotyped with the 150K-SNP panel. A single-trait fixed and random Circuitous Probability Unification model was performed to detect candidate genes for eight longevity traits. Generally, all eight longevity traits had low heritabilities, ranging from 0.038 for total productive life and herd life to 0.090 for days from the first calving to the end of first lactation or culling. High genetic correlations were observed among the traits within the same definition group: from 0.946 to 0.997 for three traits reflecting full lifespan and from 0.666 to 0.997 for five traits reflecting partial productive life. Genetic correlations between two groups of traits ranged from 0.648 to 0.963, and increased gradually with the extension of lactations number regarding the partial productive life traits. A total of 55 SNPs located on 25 chromosomes were found genome-wide significantly associated with longevity, in which 12 SNPs were associated with more than one trait, even across traits of different definition groups. This is the first study to investigate the genetic architecture of longevity representing both full and the partial lifespan simultaneously, which will assist the selection of an appropriate trait definition for genetic improvement of longevity. Because of high genetic correlations with the full lifespan traits and higher heritability, the partial productive life trait measured as the days from the first calving to the end of the third lactation or culling could be a good alternative for early selection on longevity. The candidate genes identified by this study, such as RPRM, GRIA3, GTF2H5, CA5A, CACNA2D1, FGF10, and DNAJA3, could be used to pinpoint causative mutations for longevity and further benefit the genomic improvement of longevity in dairy cattle.

Background: Staphylococcus aureus (S. aureus) mastitis is one of the most difficult diseases to treat in lactating dairy cows worldwide. S. aureus with different lineages leads to different host immune responses. Long non-coding RNAs (lncRNAs) are reported to be widely involved in the progress of inflammation. However, no research has identified stable lncRNAs among different S. aureus strain infections. In addition, folic acid (FA) can effectively reduce inflammation, and whether the inflammatory response caused by S. aureus can be reduced by FA remains to be explored. Methods: lncRNA transcripts were identified from Holstein mammary gland tissues infected with different concentrations of S. aureus (in vivo) and mammary alveolar cells (Mac-T cells, in vitro) challenged with different S. aureus strains. Differentially expressed (DE) lncRNAs were evaluated, and stable DE lncRNAs were identified in vivo and in vitro. On the basis of the gene sequence conservation and function conservation across species, key lncRNAs with the function of potentially immune regulation were retained for further analysis. The function of FA on inflammation induced by S. aureus challenge was also investigated. Then, the association analysis between these keys lncRNA transcripts and hematological parameters (HPs) was carried out. Lastly, the knockdown and overexpression of the important lncRNA were performed to validate the gene function on the regulation of cell immune response. Results: Linear regression analysis showed a significant correlation between the expression levels of lncRNA shared by mammary tissue and Mac-T cells (P < 0.001, R² = 0.3517). lncRNAs PRANCR and TNK2–AS1 could be regarded as stable markers associated with bovine S. aureus mastitis. Several HPs could be influenced by SNPs around lncRNAs PRANCR and TNK2–AS1. The results of gene function validation showed PRANCR regulates the mRNA expression of SELPLG and ITGB2 within the S. aureus infection pathway and the Mac-T cells apoptosis. In addition, FA regulated the expression change of DE lncRNA involved in toxin metabolism and inflammation to fight against S. aureus infection. Conclusions: The remarkable association between SNPs around these two lncRNAs and partial HP indicates the potentially important role of PRANCR and TNK2–AS1 in immune regulation. Stable DE lncRNAs PRANCR and TNK2–AS1 can be regarded as potential targets for the prevention of bovine S. aureus mastitis. FA supplementation can reduce the negative effect of S. aureus challenge by regulating the expression of lncRNAs.