Alternative subclinical mastitis traits for genetic evaluation in dairy cattle

Alternative Traits for Genetic Evaluation of Mastitis Based on Lifetime Merit

Genetic selection has achieved little progress in reducing mastitis incidence. Mastitis traits are problematic due to the lack of sensitivity of the data and reliance on clinical diagnosis, often missing subclinical cases, and/or on monthly somatic cell count (SCC) measurements. The current measure for mastitis is the lactation average of the somatic cells score (LSCS). We studied two datasets: (1) 148 heifers divided into non-intramammary infected, sub-clinically infected and clinical mastitis groups; (2) data from 89,601 heifers from Israeli Holsteins through the same period divided into "udder healthy" (UH) and "non-healthy" (UNH) by a threshold of SCC 120,000 cells/mL in all nine monthly milk recordings. In study 1, non-infected heifers had significantly (p < 0.05) more partum, production days and overall lifetime milk production compared to clinical and sub-clinically infected. In study 2, UH heifers (20.3%) had significantly higher (p < 0.01) lifetime milk, production days, and lactations. Subdividing datasets by sires, the same analyses detected differences in percentages of UH daughters between the sire groups. Lifetime milk production correlated (r = +0.83, p < 0.001) with udder health status. SCC threshold of less than 120,000 cells/mL during all first lactation measurements indicated healthy udder, providing a valuable insight that this dichotomous trait is advantageous for calculating lifetime net-merit index (NM$) over LSCS.

High-throughput analysis of lncRNA in cows with naturally infected Staphylococcus aureus mammary gland

LncRNA (long non-coding RNA) is an RNA molecule with a length between 200 and 100,000 nt. It does not encode proteins and is involved in a variety of intracellular processes, becoming a research hotspot of genetics. To identify key lncRNAs associated with dairy mastitis, we collected mammary epithelial tissue samples of Normal disease-free Holstein cows (HCN) and unhealthy Holstein cows with Staphylococcus aureus (HCU) and performed RNA sequencing (RNA-seq) on the samples. A total of 270 differentially expressed lncRNAs and 500 differentially expressed mRNAs were identified by high-throughput sequencing and bioinformatics analysis. Furthermore, Hydrolase activity is the most enriched in GO, and ErbB signaling pathway is significantly enriched in KEGG. In addition, through qPCR validation of 5 candidate lncRNAs in HCN and HCU, four differentially expressed lncRNAs MSTRG.498, MSTRG57.1, MSTRG.41.1 and MSTRG 124.1 were confirmed to have significant differentially expressed in cow mastitis. Also, lncRNA MSTRG.498 and its target gene, SMC4, might directly or indirectly play a role in cow mastitis. The regulatory network of lncRNA-miRNA-mRNA has been inferred from a bioinformatics perspective, which may assist understand the underlying molecular mechanism of lncRNAs involved in regulating mastitis in cows. Our findings will provide meaningful resources for further research on the regulatory function of lncRNAs in cow mastitis.

Genomic-based genetic parameters for milkability traits derived from automatic milking systems in North American Holstein cattle

The number of dairy farms adopting automatic milking systems (AMS) has considerably increased around the world aiming to reduce labor costs, improve cow welfare, increase overall performance, and generate a large amount of daily data, including production, behavior, health, and milk quality records. In this context, this study aimed to (1) estimate genomic-based variance components for milkability traits derived from AMS in North American Holstein cattle based on random regression models; and (2) derive and estimate genetic parameters for novel behavioral indicators based on AMS-derived data. A total of 1,752,713 daily records collected using 36 milking robot stations and 70,958 test-day records from 4,118 genotyped Holstein cows were used in this study. A total of 57,600 SNP remained after quality control. The daily-measured traits evaluated were milk yield (MY, kg), somatic cell score (SCS, score unit), milk electrical conductivity (EC, mS), milking efficiency (ME, kg/min), average milk flow rate (FR, kg/min), maximum milk flow rate (FRM, kg/min), milking time (MT, min), milking failures (MFAIL), and milking refusals (MREF). Variance components and genetic parameters for MY, SCS, ME, FR, FRM, MT, and EC were estimated using the AIREMLF90 software under a random regression model fitting a third-order Legendre orthogonal polynomial. A threshold Bayesian model using the THRGIBBS1F90 software was used for genetically evaluating MFAIL and MREF. The daily heritability estimates across days in milk (DIM) ranged from 0.07 to 0.28 for MY, 0.02 to 0.08 for SCS, 0.38 to 0.49 for EC, 0.45 to 0.56 for ME, 0.43 to 0.52 for FR, 0.47 to 0.58 for FRM, and 0.22 to 0.28 for MT. The estimates of heritability (± SD) for MFAIL and MREF were 0.02 ± 0.01 and 0.09 ± 0.01, respectively. Slight differences in the genetic correlations were observed across DIM for each trait. Strong and positive genetic correlations were observed among ME, FR, and FRM, with estimates ranging from 0.94 to 0.99. Also, moderate to high and negative genetic correlations (ranging from -0.48 to -0.86) were observed between MT and other traits such as SCS, ME, FR, and FRM. The genetic correlation (± SD) between MFAIL and MREF was 0.25 ± 0.02, indicating that both traits are influenced by different sets of genes. High and negative genetic correlations were observed between MFAIL and FR (-0.58 ± 0.02) and MFAIL and FRM (-0.56 ± 0.02), indicating that cows with more MFAIL are those with lower FR. The use of random regression models is a useful alternative for genetically evaluating AMS-derived traits measured throughout the lactation. All the milkability traits evaluated in this study are heritable and have demonstrated selective potential, suggesting that their use in dairy cattle breeding programs can improve dairy production efficiency in AMS.

Impact of different search strategies on the results of a meta-analysis for mastitis in dairy cattle

Abstract Poorly formulated search strategies can have great influences on the results of a meta-analysis, since it directly impacts the amount and adherence to the theme of the works used for study, therefore, the formulation of a consistent and functional search strategy is essential for the review system to achieve its goals. The objective of this work was to study the impact of different search strategies in a systematic review by performing a meta-analysis to estimate heritability for the mastitis trait in dairy cattle. Once the search strategies were defined, the searches carried out in the Web of science, Scopus, Scielo and Pubmed databases returned 921 studies from which, after going through the identification, selection, eligibility and inclusion processes, 25 studies were selected. Withdrawals from selected articles, 26 heritability estimates were used in the meta-analysis. A random effect model was used, with all analyzes performed by the R program, through the Metafor package. The estimates obtained through the combined statistics of studies for mastitis, presented values of low magnitude (0.05 and 0.06). The effects of search strategies have a significant impact on the meta-analysis estimates produced.

Impacto de diferentes estratégias de busca sob os resultados de uma metanálise para a característica mastite em gado de leite

Resumo Estratégias de buscas mal formuladas podem apresentar grandes influências nos resultados de uma metanálise, uma vez que impacta diretamente na quantidade e aderência ao tema dos trabalhos utilizados para estudo, portanto, a formulação de uma estratégia de busca consistente e funcional é fundamental para que a revisão sistemática atinja seus objetivos. O objetivo deste trabalho foi estudar o impacto de diferentes estratégias de busca em uma revisão sistemática por meio da realização de uma metanálise para estimação de herdabilidade para a característica mastite em gado de leite. Uma vez definidas as estratégias de busca, as pesquisas realizadas nas bases de dados Web of science, Scopus, Scielo e Pubmed retornaram 921 estudos dos quais, após passarem pelos processos de identificação, seleção, elegibilidade e inclusão, 25 estudos foram selecionados. Retiradas dos artigos selecionados, 26 estimativas de herdabilidades foram utilizadas na realização da metanálise. Utilizou-se um modelo de efeito aleatório, sendo todas as análises realizadas pelo programa R, por meio do pacote Metafor. As estimativas obtidas através da estatística combinada de estudos para mastite, apresentou valores de baixa magnitude (0,05 e 0,06). Os efeitos das estratégias de busca têm impacto significativo nas estimativas de metanálise produzidas.

Determining extracellular vesicles properties and miRNA cargo variability in bovine milk from healthy cows and cows undergoing subclinical mastitis

Background

Subclinical mastitis, the inflammation of the mammary gland lacking clinical symptoms, is one of the most prevalent and costly diseases in dairy farming worldwide. Milk microRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) have been proposed as potential biomarkers of different mammary gland conditions, including subclinical mastitis. However, little is known about the robustness of EVs analysis regarding sampling time-point and natural infections. To estimate the reliability of EVs measurements in raw bovine milk, we first evaluated changes in EVs size and concentration using Tunable Resistive Pulse Sensing (TRPS) during three consecutive days of sampling. Then, we analysed daily differences in miRNA cargo using small RNA-seq. Finally, we compared milk EVs differences from naturally infected udder quarters with their healthy adjacent quarters and quarters from uninfected udders, respectively.

Results

We found that the milk EV miRNA cargo was very stable over the course of three days regardless of the health status of the quarter, and that infected quarters did not induce relevant changes in milk EVs of adjacent healthy quarters. Chronic subclinical mastitis induced changes in milk EV miRNA cargo, but neither in EVs size nor concentration. We observed that the changes in immunoregulatory miRNAs in quarters with chronic subclinical mastitis were cow-individual, however, the most upregulated miRNA was bta-miR-223-3p across all individuals.

Conclusions

Our results showed that the miRNA profile and particle size characteristics remained constant throughout consecutive days, suggesting that miRNAs packed in EVs are physiological state-specific. In addition, infected quarters were solely affected while adjacent healthy quarters remained unaffected. Finally, the cow-individual miRNA changes pointed towards infection-specific alterations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08377-z.

Milk Somatic Cell Count and Polymorphonuclear Cells in Healthy Quarters of Cows That Underwent Blanket and Selective Dry Therapy: An Italian Case Study

The incidence of mastitis increases with parity in dairy cattle and multiparous cows are often treated at drying off to limit the risk of udder health issues and support mammary gland tissues recovery. Milk somatic cells count (SCC, cells/mL) comprises different white blood cells fractions and is worldwide used to monitor and genetically improve udder health. Nevertheless, only certain SCC fractions increase when an udder inflammation occurs. Considering that antibiotic use for preventive purposes will be forbidden in 2022, we compared two different dry therapy protocols, blanket (BDCT) and selective (SDCT), on different SCC fractions in healthy quarters milk. Multiparous Holstein cows were enrolled in a randomized controlled trial and SCC, neutrophils, macrophages, lymphocytes, polymorphonuclear cells (PMN) and differential somatic cell count (DSCC) recorded after the experimental drying off were available. Significant differences were observed between the two protocols, with more favorable parameters in BDCT than SDCT cows. Results showed that moving from BDCT to SDCT is expected to significantly increase some SCC fractions, such as PMN, in healthy quarters. The baseline SCC level at the onset of lactation was greater in cows under SDCT than BDCT. Although not significant, clinical mastitis prevalence was numerically lower in BDCT (7.32%) than SDCT (8.62%). In this study we referred to a limited number of cows, but still findings will be useful to improve the knowledge on the impact of SDCT on milk SCC fractions in healthy quarters.

Relationship of daily total somatic cell output with somatic cell concentration and clinical mastitis

•

Measures of somatic cells in milk typically reflect concentration of cells but not total cell output.

•

Milk yield has a dilution effect, and elevated concentration of somatic cells in early and late lactation partly reflects lower yield.

•

Consideration of both somatic cell concentration and total daily somatic cell output may improve mastitis detection.

Somatic cell count (SCC) measures the concentration of somatic cells in milk and is used as a mastitis diagnostic tool. It is plausible that variation in milk yield could alter the relationship between SCC and mastitis status. Our objective was to evaluate total daily SCC output as a predictor of clinical mastitis. Data included 37,035 test-day records from 4,179 lactations of 1,679 cows and 1,286 mastitis events from an experimental herd. Daily total SCC was derived by multiplying SCC by daily milk yield in milliliters and transformed to daily total somatic cell score (DTSCS) via a log₂ transformation. Milk yield, SCS, and DTSCS were evaluated with mixed models that included the proximity of a mastitis event to the test date and days in milk as the main fixed effects. A second series of logistic regression was conducted that considered mastitis (1 = mastitis occurred during a test interval; 0 = no mastitis) as the dependent variable with milk yield, SCS, and DTSCS the main independent effects. Least squares means for test dates associated with mastitis-free lactations were 2.43 and 2.25 for SCS and DTSCS, respectively. The corresponding values were 5.96 for SCS and 5.66 for DTSCS for the week of a mastitis event. Whereas SCS declined rapidly in early lactation and then increased steadily thereafter, DTSCS was lowest in early lactation and increased by a proportionally smaller amount throughout lactation. Including both SCS and DTSCS in the same model improved the logistic regression model fit over a model with SCS only. Dilution effects from milk yield influence SCS, and consideration of DTSCS in management and genetic selection schemes could improve mastitis detection and resistance.

Genetic Parameters Estimation of Milking Traits in Polish Holstein-Friesians Based on Automatic Milking System Data

The automatic milking system (AMS) provides a large amount of information characterizing the course of each milking cow, which is not available in the conventional system. The aim of our study was to estimate heritability and genetic correlations for milk yield (MY), milking frequency (MF), and speed (MS) for 1713 Polish Holstein-Friesian primiparous cows milked in barns with an AMS. Daily heritability indicators estimated using second-order Legendre polynomials and Random Regression Models showed high variation during lactation, ranging 0.131-0.345 for MY, 0.153-0.322 for MF, and 0.336-0.493 for MS. The rates of genetic correlation between traits ranged: 0.561-0.929 for MY-MF, (-0.255)-0.090 for MF-MS, (-0.174)-0.020 for MY-MS. It is possible to carry out effective selection for milking speed, which provides an opportunity to increase the number of cows per milking robot, and thus increase the profitability of production in the herd. The results proved that selection for milk yield and daily milking frequency is also feasible. The research showed a high, positive genetic correlation between milking frequency and milk yield, which allows us to conclude that preferring breeding cows with a natural tendency to frequent visits to the milking robot should indirectly improve the genetic basis of milking.

Expression analysis of candidate genes for chronic subclinical mastitis in Norwegian Red cattle

Chronic subclinical mastitis (SCM) is characterized by a long-term inflammation in the udder with high somatic cell count (SCC) in milk. Previously, several novel alternative SCM traits for Norwegian Red (NR) cattle have been defined to improve breeding strategies against chronic SCM. Quantitative trait loci and candidate genes affecting chronic SCM in NR have been identified. The aim of this study was to analyze the expression profiles of 14 selected candidate genes (RAD17, ACOT2, ACOT4, FOS, CXCL1, CXCL8, CCNB1, CDK7, TGFB3, SEL1L, STAT4, C6, GLI2, and SLC18A2). Twenty healthy NR cows with official genomic estimated breeding values (GEBV) for lactation average somatic cell scores (LSCS) were selected. Ten cows had high GEBV for LSCS (cows with low probability to have high SCC in milk during lactation) and 10 cows had low GEBV for LSCS (cows with high probability of having high SCC in milk). We isolated RNA from unstimulated peripheral blood mononuclear cells from these. Two out of the 14 analyzed genes showed significantly different results between groups. The group with high GEBV for LSCS displayed significantly higher expression of the CXCL1 gene than the low GEBV group. Grouping by lactation stage revealed significant differential expression of the FOS gene, with higher expression in early lactation (2-3 mo after calving) compared with late lactation (7-8 mo after calving). In addition, flow cytometry was performed on the peripheral blood mononuclear cells samples to analyze if number and type of isolated cells influenced the gene expression in the groups. The results in the current study provide identified genes that can be considered as possible candidate genes for chronic SCM in NR cows.

A genetic study of new udder health indicator traits with data from automatic milking systems

The current study aimed to investigate new udder health traits based on data from automatic milking systems (AMS) for use in routine genetic evaluations. Data were from 77 commercial herds; out of these, 24 had equipment for measuring online cell count (OCC), whereas all had data on electrical conductivity (EC). A total of 4,714 Norwegian Red dairy cows and 2,363,928 milkings were included in the genetic analyses. Electrical conductivity was available on quarter level for each milking, whereas OCC was measured per milking. The AMS traits analyzed were log-transformed online cell count (lnOCC), maximum conductivity (ECmax), mean conductivity (ECmean), elevated mastitis risk (EMR), and log-transformed EMR (lnEMR). In addition, lactation mean somatic cell score (LSCS) was collected from the Norwegian dairy herd recording system. Elevated mastitis risk expresses the probability of a cow having mastitis and was calculated from smoothed lnOCC values according to individual trend and level of the OCC curve. The udder health traits from AMS were analyzed as repeated milkings from 30 to 320 DIM, and LSCS as repeated parities. In addition, both ECmax and lnOCC were analyzed as multiple traits by splitting the lactation into 5 periods. (Co)variance components were estimated from bivariate mixed linear animal models, and investigated traits showed genetic variation. Estimated heritabilities of ECmean, ECmax, and lnEMR were 0.35, 0.23, and 0.12, respectively, whereas EMR and lnOCC both showed heritabilities of 0.09. Heritability varied between periods of lactation, from 0.04 to 0.13 for lnOCC and from 0.12 to 0.27 for ECmax, although standard errors of certain periods were large. Genetic correlations among the AMS traits ranged from 0 to 0.99. The genetic correlations between EC-based traits and OCC-based traits in AMS were 0. Genetic correlations with LSCS were favorable, ranging from 0.37 to 0.80 (±0.11-0.22). The strongest correlation (0.80 ± 0.13) was found between LSCS and lnEMR. Results question the value of ECmax and ECmean as indicators of udder health in genetic evaluations and suggest OCC to be more valuable in this manner. This study demonstrates a potential of using AMS data as additional information on udder health for genetic evaluations, although further investigation is recommended before these traits can be implemented.

Identification of candidate genes affecting chronic subclinical mastitis in Norwegian Red cattle: combining genome-wide association study, topologically associated domains and pathway enrichment analysis

The aim of this study was to identify genes associated with chronic subclinical mastitis (SCM) in Norwegian Red (NR) cattle. Twelve SCM traits defined based on fixed threshold for test-day somatic cell count (SCC) were, together with lactation-average somatic cell score (LSCS) used for association and pathway enrichment analyses. A GWAS was performed on 3795 genotyped NR bulls with 777K SNP data and phenotypic information from 7 300 847 test-day SCC observations from 3 543 764 cows. At 5% chromosome-wide significance level 36 unique SNP were detected to be associated with one or more of the traits. These SNPs were analysed for linked genes using genomic positions of topologically associated domains (TAD). For the SCM traits with SCC >50 000 and >100 000 cells/ml on two test-days in a row and LSCS, the same top significant genes were identified - checkpoint clamp loader component (RAD17) and cyclin B1 (CCNB1). The SCM traits with SCC >250 000, 300 000, 350 000 or 400 000 cells/ml on two test-days in a row and D400 (number of days before the first case with SCC >400 000 cells/ml) displayed similar top significant genes: acyl-CoA thioesterase 2 and 4 (ACOT2; ACOT4). For the traits SCM200_3 (SCC >200 000 cells/ml on three test-days in a row) and SCM150, SCM200 (SCC >150 000; 200 000 cells/ml on two test-days in a row) a group of chemokine (C-X-C motif) ligand genes and the Fos proto-oncogene, AP-1 transcription factor subunit (FOS) gene, were identified. Further functional studies of these identified candidate genes are necessary to clarify their actual role in development of chronic SCM in NR cattle.