Multiple trait genetic evaluation of clinical mastitis in three dairy cattle breeds

Exploring the impact of biosecurity measures on somatic cell score in mountain dairy farms considering the CLASSYFARM system

Biosecurity plays a crucial role in preventing the introduction of infectious diseases to a herd as well as the spread of diseases within or between animals and herds. In particular, biosecurity measures are crucial for maintaining animal health and reducing the need for the application of antibiotic substances for fighting the rising antibiotic resistance. The object of this cross-sectional study was to investigate the presence of biosecurity measures and their association with milk quality parameters, with a special focus on somatic cell count (SCC) - an indicator for udder health -, in small-scale mountain dairy farms. Therefore, for the very first time, the CLASSYFARM system was considered, which is a computer platform integrated into the Italian national veterinary portal, that processes a significant amount of data from various sources collected in the field or from other information systems (e.g. animal welfare, health status, biosecurity, antimicrobial use, slaughterhouse information). A total of 169 dairy farms were included in the study. Biosecurity measures, based on 15 questions required in the CLASSYFARM welfare assessment protocol, as well as information about husbandry systems, milking systems and pasture practices were gathered and combined with milk yield data, provided by the South Tyrolean dairy association. Farms only scored 44.00 points on average in a scoring system from 0.00 to 100.00 points that was be able to summarize 15 different biosecurity measures in one index. Our results show a clear negative correlation (-0.713) between the biosecurity index and somatic cell score (SCS) indicating that a higher level of biosecurity, which reflects the presence of biosecurity measures within a farm, is associated with lower SCC levels. Furthermore, we found significant correlations between SCS and milk production (-0.629), confirming that udder health is linked to higher milk production. Fat, protein, and the fat-to-protein ratio showed a positive correlation with SCS (0.281, 0.146, 0.106), likely to be caused by a concentration shift effect (dilution effect). Husbandry system, breed, milking system, and pasture practices seem to have an impact as well, but the main factor was the biosecurity score. This study highlights the importance of implementing biosecurity measures for ensuring animal health and thus productivity and quality in milk production, even in small-scale farms, which are characterized by limited structure availability and smaller herds compared to big dairy enterprises in the lowlands.

A comprehensive meta-analysis of genetic parameters for resilience and productivity indicator traits in Holstein cattle

Selection for resilience indicator (RIND) traits in Holstein cattle is becoming an important breeding objective as the worldwide population is expected to be exposed to increased environmental stressors due to both climate change and changing industry standards. However, genetic correlations between RIND and productivity indicator (PIND) traits, which are already being selected for and have the most economic value, are often unfavorable. As a result, it is necessary to fully understand these genetic relationships when incorporating novel traits into selection indices, so that informed decisions can be made to fully optimize selection for both groups of traits. In the past 2 decades, there have been many estimates of RIND traits published in the literature, albeit in small populations. To provide valuable pooled summary estimates, a random-effects meta-analysis was conducted for heritability and genetic correlation estimates for PIND and RIND traits in worldwide Holstein cattle. In total, 926 heritability estimates for 9 PIND and 27 RIND traits, along with 362 estimates of genetic correlation (PIND × RIND traits) were collected. Resilience indicator traits were grouped into the following subgroups: Metabolic Diseases, Hoof Health, Udder Health, Fertility, Heat Tolerance, Longevity, and Other. Pooled estimates of heritability for PIND traits ranged from 0.201 ± 0.05 (energy-corrected milk) to 0.377 ± 0.06 (protein content), while pooled estimates of heritability for RIND traits ranged from 0.032 ± 0.02 (incidence of lameness, incidence of milk fever) to 0.497 ± 0.05 (measures of body weight). Pooled estimates of genetic correlations ranged from -0.360 ± 0.25 (protein content vs. milk acetone concentration) to 0.535 ± 0.72 (measures of fat-to-protein ratio vs. milk acetone concentration). Additionally, out of 243 potential genetic correlations between PIND and RIND traits that could have been reported, only 40 had enough published estimates to implement the meta-analysis model. Our results confirmed that the interactions between PIND and RIND traits are complex, and all relationships should be evaluated when incorporating novel traits into selection indices. This study provides a valuable reference for breeders looking to incorporate RIND traits for Holstein cattle into selection indices.

Genetic aspects of immunoglobulins and cyclophilin A in milk as potential indicators of mastitis resistance in Holstein cows

Mastitis is one of the most frequent and costly diseases affecting dairy cattle. Natural antibodies (immunoglobulins) and cyclophilin A (CyPA), the most abundant member of the family of peptidyl prolyl cis/trans isomerases, in milk may serve as indicators of mastitis resistance in dairy cattle. However, genetic information for CyPA is not available, and knowledge on the genetic and nongenetic relationships between these immune-related traits and somatic cell score (SCS) and milk yield in dairy cattle is sparse. Therefore, we aimed to comprehensively evaluate whether immune-related traits consisting of 5 Ig classes (IgG, IgG1, IgG2, IgA, and IgM) and CyPA in the test-day milk of Holstein cows can be used as genetic indicators of mastitis resistance by evaluating the genetic and nongenetic relationships with SCS in milk. The nongenetic factors affecting immune-related traits and the effects of these traits on SCS were evaluated. Furthermore, the genetic parameters of immune-related traits according to health status and genetic relationships under different SCS environments were estimated. All immune-related traits were significantly associated with SCS and directly proportional. Additionally, evaluation using a classification tree revealed that IgA, IgG2, and IgG were associated with SCS levels. Genetic factor analyses indicated that heritability estimates were low for CyPA (0.08) but moderate for IgG (0.37), IgA (0.44), and IgM (0.44), with positive genetic correlations among Ig (0.25-0.96). We also evaluated the differences in milk yield and SCS of cows between the low and high groups according to their sires' estimated breeding value for immune-related traits. In the high group, IgA had a significantly lower SCS in milk at 7 to 30 d compared with that in the low group. Furthermore, the Ig in milk had high positive genetic correlations between healthy and infected conditions (0.82-0.99), suggesting that Ig in milk under healthy conditions could interact with those under infected conditions, owing to the genetic ability based on the level of Ig in milk. Thus, Ig in milk are potential indicators for the genetic selection of mastitis resistance. However, because only the relationship between immune-related traits and SCS was investigated in this study, further study on the relationship between clinical mastitis and Ig in milk is needed before Ig can be used as an indicator of mastitis resistance.

Combined single-step evaluation of functional longevity of dairy cows including correlated traits

BACKGROUND

For years, multiple trait genetic evaluations have been used to increase the accuracy of estimated breeding values (EBV) using information from correlated traits. In France, accurate approximations of multiple trait evaluations were implemented for traits that are described by different models by combining the results of univariate best linear unbiased prediction (BLUP) evaluations. Functional longevity (FL) is the trait that has most benefited from this approach. Currently, with many single-step (SS) evaluations, only univariate FL evaluations can be run. The aim of this study was to implement a "combined" SS (CSS) evaluation that extends the "combined" BLUP evaluation to obtain more accurate genomic (G) EBV for FL when information from five correlated traits (somatic cell score, clinical mastitis, conception rate for heifers and cows, and udder depth) is added.

RESULTS

GEBV obtained from univariate SS (USS) evaluations and from a CSS evaluation were compared. The correlations between these GEBV showed the benefits of including information from correlated traits. Indeed, a CSS evaluation run without any performances on FL showed that the indirect information from correlated traits to evaluate FL was substantial. USS and CSS evaluations that mimic SS evaluations with data available in 2016 were compared. For each evaluation separately, the GEBV were sorted and then split into 10 consecutive groups (deciles). Survival curves were calculated for each group, based on the observed productive life of these cows as known in 2021. Regardless of their genotyping status, the worst group of heifers based on their GEBV in 2016 was well identified in the CSS evaluation and they had a substantially shorter herd life, while those in the best heifer group had a longer herd life. The gaps between groups were more important for the genotyped than the ungenotyped heifers, which indicates better prediction of future survival.

CONCLUSIONS

A CSS evaluation is an efficient tool to improve FL. It allows a proper combination of information on functional traits that influence culling. In contrast, because of the strong selection intensity on young bulls for functional traits, the benefit of such a "combined" evaluation of functional traits is more modest for these males.

Molecular regulatory mechanism of key LncRNAs in subclinical mastitic cows with folic acid supplementation

BACKGROUND

Folic acid is a water-soluble B vitamin (B9), which is closely related to the body's immune and other metabolic pathways. The folic acid synthesized by rumen microbes has been unable to meet the needs of high-yielding dairy cows. The incidence rate of subclinical mastitis in dairy herds worldwide ranged between 25%~65% with no obvious symptoms, but it significantly causes a decrease in lactation and milk quality. Therefore, this study aims at exploring the effects of folic acid supplementation on the expression profile of lncRNAs, exploring the molecular mechanism by which lncRNAs regulate immunity in subclinical mastitic dairy cows.

RESULTS

The analysis identified a total of 4384 lncRNA transcripts. Subsequently, differentially expressed lncRNAs in the comparison of two groups (SF vs. SC, HF vs. HC) were identified to be 84 and 55 respectively. Furthermore, the weighted gene co-expression network analysis (WGCNA) and the KEGG enrichment analysis result showed that folic acid supplementation affects inflammation and immune response-related pathways. The two groups have few pathways in common. One important lncRNA MSTRG.11108.1 and its target genes (ICAM1, CCL3, CCL4, etc.) were involved in immune-related pathways. Finally, through integrated analysis of lncRNAs with GWAS data and animal QTL database, we found that differential lncRNA and its target genes could be significantly enriched in SNPs and QTLs related to somatic cell count (SCC) and mastitis, such as MSTRG.11108.1 and its target gene ICAM1, CXCL3, GRO1.

CONCLUSIONS

For subclinical mastitic cows, folic acid supplementation can significantly affect the expression of immune-related pathway genes such as ICAM1 by regulating lncRNAs MSTRG.11108.1, thereby affecting related immune phenotypes. Our findings laid a ground foundation for theoretical and practical application for feeding folic acid supplementation in subclinical mastitic cows.

Genetical analysis of mastitis and reproductive traits in first-parity Holstein cows using standard and structural equation modelling

The present study aimed to investigate the causal relationships between clinical mastitis and some reproductive traits, including success at first insemination (SFI), the number of inseminations to pregnancy (INS), the interval from calving to first service (CTFS), first and last service interval (IFL), and open days (OD) in first-parity Holstein cows. For this purpose, the records of 58 281 first parity Holstein cows were analysed. These data sets were collected from 17 large dairy herds from 2008 to 2017. Recursive Mixed Models (RMMs) were applied and compared with the estimations under Standard Mixed Models. Then, one trivariate and three bivariate Gaussian-threshold models were used for the analyses. Recursive models were applied, considering that clinical mastitis can influence fertility traits. Mastitis is considered a covariate for the reproductive traits to determine their causal relationship. The results of this study indicated that causal effects of mastitis on SFI (on the observed scale, %), CTFS, IFL, OD, and INS were -5.7%, 3.3 days, 12.27 days, seven days, and 0.26 services, respectively. The estimated structural coefficients of the recursive models in the first parity imply that mastitis significantly lengthened the fertility interval and decreased the conception rate. In addition, genetic, residual, and phenotypic correlations between mastitis and the reproductive traits under both models were statistically significant. Results of genetic correlations between mastitis and fertility traits suggest that more incidence of mastitis during lactation is related to the delays in the heat show and pregnancy rate after insemination. In summary, considering the causal effects under RMMs may be advantageous to comprehend complicated relationships between complex traits better.

Genetic parameters of colostrum and calf serum antibodies in Swedish dairy cattle

Background

A sufficient IgG content in the colostrum is essential for the newborn calf, as it provides passive immunity which substantially affects the probability of survival during rearing. Failure of passive transfer (FPT) occurs when a calf does not absorb enough antibodies from the colostrum and is defined by an IgG concentration in calf serum lower than 10 g/L. Apart from delayed access to colostrum, FPT can be due to a low production of IgG in the mother or poor IgG absorption by the calf. The aim of this study was to estimate the genetic background of antibody levels and indicator traits for antibodies in the colostrum and calf serum, and their correlation with milk production.

Results

Colostrum data were available for 1340 dairy cows with at least one calving and calf serum data were available for 886 calves from these cows. Indicator traits for antibody concentrations were estimated using refractometry (a digital Brix refractometer for colostrum and an optical refractometer for serum), and enzyme-linked immunosorbent assays (ELISA) were used to determine the levels of total IgG and natural antibodies (NAb) of various antibody isotypes in the colostrum and calf serum. Colostrum traits had heritabilities ranging from 0.16 to 0.31 with repeatabilities ranging from 0.21 to 0.55. Brix percentages had positive genetic correlations with all colostrum antibody traits including total IgG (0.68). Calf serum antibody concentrations had heritabilities ranging from 0.25 to 0.59, with a significant maternal effect accounting for 17 to 27% of the variance. When later in life calves produced their first lactation, the lactation average somatic cell score was found to be negatively correlated with NAb levels in calf serum.

Conclusions

Our results suggest that antibody levels in the colostrum and calf serum can be increased by means of selection.

Estimation of milkability breeding values and variance components for Italian Holstein

•

The national genetic evaluation for milkability in Italian Holstein was revised.

•

A stricter data editing and a different approach to the phenotype, from ratio to single observations, were applied.

•

A different statistical model was used, changing from a linear to a threshold model.

•

Higher genomic reliability was achieved with the revised model.

•

The revised model provides more reliable breeding values for decision-making at the farm level.

The importance of milkability as a trait is growing because of the need to efficiently use labor and machinery; therefore, it is crucial to update the statistical model for the trait to improve the accuracy of the estimated breeding values, and thus provide a more accurate tool for decision-making at the farm level. In the Italian Holstein Friesian cattle population, milkability is recorded twice a year by the milk recording system as a binary trait (slow, coded as 2, or not slow, coded as 1). Data consisted of 7,862,371 records from 2,945,249 cows collected between 2004 and 2021. A single-trait threshold animal model with repeated measures was used, with parity, days in milk class, calving season, and regression of production (fat + protein grams) within days in milk class as fixed effects and herd-year-season of recording, permanent environment, and animal as random effects. The results for heritability and repeatability were 0.275 and 0.5, estimated with the Gibbs sampler THRGIBBS1F90. Genomic validation, carried out using genotyped proven bulls born before 2009 as the training set, gave a result of 0.386 for reliability. The genetic correlations of this trait confirmed that both extremes of the estimated breeding value must be treated cautiously, because correlations with important traits such as mastitis resistance, body condition score, and teat length are unfavorable.

Associations between differential somatic cell count and milk yield, quality, and technological characteristics in Holstein cows

The aim of this study was to investigate the associations between differential somatic cell count (DSCC) and milk quality and udder health traits, and for the first time, between DSCC and milk coagulation properties and cheesemaking traits in a population of 1,264 Holstein cows reared in northern Italy. Differential somatic cell count represents the combined proportions of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in the total somatic cell count (SCC), with macrophages (MAC) making up the remaining proportion. The milk traits investigated in this study were milk yield (MY), 8 traits related to milk composition and quality (fat, protein, casein, casein index, lactose, urea, pH, and milk conductivity), 9 milk coagulation traits [3 milk coagulation properties (MCP) and 6 curd firming (CF) traits], 7 cheesemaking traits, 3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits. A linear mixed model was fitted to explore the associations between SCS combined with DSCC and the aforementioned milk traits. An additional model was run, which included DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the percentage of PMN-LYM and MAC by SCC in the milk for each cow in the data set. The unfavorable association between SCS and milk quality and technological traits was confirmed. Increased DSCC was instead associated with a linear increase in MY, casein index, and lactose proportion and a linear decrease in milk fat and milk conductivity. Accordingly, DSCC was favorably associated with all MCP and CF traits (with the exception of the time needed to achieve maximum, CF), particularly with rennet coagulation time, and it always displayed linear relationships. Differential somatic cell count was also positively associated with the recovery of milk nutrients in the curd (protein, fat, and energy), which increased linearly with increasing DSCC. The PMN-LYM count was rarely associated with milk traits, even though the pattern observed confirmed the results obtained when both SCS and DSCC were included in the model. The MAC count, however, showed the opposite pattern: MY, casein index, and lactose percentage decreased and milk conductivity increased with an increasing MAC count. No significant association was found between PMN-LYM count and MCP, CF, CY, and REC traits, whereas MAC count was unfavorably associated with MCP, CF traits, some CY traits, and all REC traits. Our results showed that the combined information derived from SCS and DSCC might be useful to monitor milk quality and cheesemaking-related traits.

Polymorphisms in TLR4 Gene Associated With Somatic Cell Score in Water Buffaloes (Bubalus bubalis)

Considering the importance of the diseases affecting the productive performance of animals in the dairy industry worldwide, it is necessary to implement tools that help to control and limit the occurrence of such diseases. As the increased somatic cell counts (SCC) are a direct expression of the inflammatory process, they are candidates to become the usual parameter for assessing udder health regarding milk quality and for monitoring mastitis incidences. Toll-Like Receptors are membrane proteins that play a key role in immunity, recognizing pathogens and, subsequently, activating immune responses. The present study was conducted to identify single nucleotide polymorphisms in the TLR4 gene of buffaloes and to analyze its associations with somatic cell counts. DNA samples of 120 Murrah buffaloes were used. The whole coding region of the TLR4 gene was amplified by polymerase chain reaction reactions and sequenced for polymorphism scanning. A total of 13 polymorphisms were identified for the sequenced regions of the TLR4, most of which are in the coding region. The association with the somatic cell score was highly significant (p < 0.001) for all identified polymorphisms of TLR4 gene (g.54621T>A, g.54429G>T, g.54407T>A, g.46616C>A, g.46613T>G, g.46612A>G, g.46611C>A, g.46609T>G, g.46541C>G, g.46526C>A, g.46516T>C, g.46376C>T, g.46372T>C). Therefore, it is suggested that the markers of the TLR4 gene can be used as molecular markers for mastitis resistance in buffaloes, due to their association with somatic cell counts.

Confirmed effects of candidate variants for milk production, udder health, and udder morphology in dairy cattle

Background

Over the last years, genome-wide association studies (GWAS) based on imputed whole-genome sequences (WGS) have been used to detect quantitative trait loci (QTL) and highlight candidate genes for important traits. However, in general this approach does not allow to validate the effects of candidate mutations or determine if they are truly causative for the trait(s) in question. To address these questions, we applied a two-step, within-breed GWAS approach on 15 traits (5 linked with milk production, 2 with udder health, and 8 with udder morphology) in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) cattle. We detected the most-promising candidate variants (CV) using imputed WGS of 2515 MON, 2203 NOR, and 6321 HOL bulls, and validated their effects in three younger populations of 23,926 MON, 9400 NOR, and 51,977 HOL cows.

Results

Bull sequence-based GWAS detected 84 QTL: 13, 10, and 30 for milk production traits; 3, 0, and 2 for somatic cell score (SCS); and 8, 2 and 16 for udder morphology traits, in MON, NOR, and HOL respectively. Five genomic regions with effects on milk production traits were shared among the three breeds whereas six (2 for production and 4 for udder morphology and health traits) had effects in two breeds. In 80 of these QTL, 855 CV were highlighted based on the significance of their effects and functional annotation. The subsequent GWAS on MON, NOR, and HOL cows validated 8, 9, and 23 QTL for production traits; 0, 0, and 1 for SCS; and 4, 1, and 8 for udder morphology traits, respectively. In 47 of the 54 confirmed QTL, the CV identified in bulls had more significant effects than single nucleotide polymorphisms (SNPs) from the standard 50K chip. The best CV for each validated QTL was located in a gene that was functionally related to production (36 QTL) or udder (9 QTL) traits.

Conclusions

Using this two-step GWAS approach, we identified and validated 54 QTL that included CV mostly located within functional candidate genes and explained up to 6.3% (udder traits) and 37% (production traits) of the genetic variance of economically important dairy traits. These CV are now included in the chip used to evaluate French dairy cattle and can be integrated into routine genomic evaluation.

Identification of novel molecular markers of mastitis caused by Staphylococcus aureus using gene expression profiling in two consecutive generations of Chinese Holstein dairy cattle

Background

Mastitis in dairy cows caused by Staphylococcus aureus is a major problem hindering economic growth in dairy farms worldwide. It is difficult to prevent or eliminate due to its asymptomatic nature and long persistence of infection. Although transcriptomic responses of bovine mammary gland cells to pathogens that cause mastitis have been studied, the common responses of peripheral blood leukocytes to S. aureus infection across two consecutive generations of dairy cattle have not been investigated.

Methods

In the current study, RNA-Seq was used to profile the transcriptomes of peripheral blood leukocytes sampled from S. aureus-infected mothers and their S. aureus-infected daughters, and also healthy non-infected mothers and their healthy daughters. Differential gene expression was evaluated as follows: 1) S. aureus-infected cows versus healthy non-infected cows (S vs. H, which include all the mothers and daughters), 2) S. aureus-infected mothers versus healthy non-infected mothers (SM vs. HM), and 3) S. aureus-infected daughters versus healthy non-infected daughters (SMD vs. HMD).

Results

Analysis of all identified expressed genes in the four groups (SM, SMD, HM, and HMD) showed that EPOR, IL9, IFNL3, CCL26, IL26 were exclusively expressed in both the HM and HMD groups, and that they were significantly (P <  0.05) enriched for the cytokine-cytokine receptor interaction pathway. A total of 17, 13 and 10 differentially expressed genes (DEGs) (FDR P_adj. < 0.1 and |FC| > 1.2) were detected in the three comparisons, respectively. DEGs with P <  0.05 and |FC| > 2 were used for functional enrichment analyses. For the S vs. H comparison, DEGs detected included CCL20, IL13 and MMP3, which are associated with the IL-17 signaling pathway. In the SM vs. HM and SMD vs. HMD comparisons, five (BLA-DQB, C1R, C2, FCGR1A, and KRT10) and six (BLA-DQB, C3AR1, CFI, FCAR, FCGR3A, and LOC10498484) genes, respectively, were involved in the S. aureus infection pathway.

Conclusions

Our study provides insights into the transcriptomic responses of bovine peripheral blood leukocytes across two generations of cattle naturally infected with S. aureus. The genes highlighted in this study could serve as expression biomarkers for mastitis and may also contain sequence variation that can be used for genetic improvement of dairy cattle for resilience to mastitis.

Genotype by environment interaction for somatic cell score in Holstein cattle of southern Brazil via reaction norms

Objective

The objective of this study was to evaluate the genetic behavior of a population of Holstein cattle in response to the variation of environmental temperature by analyzing the effects of genotype by environment interaction (GEI) through reaction norms for the somatic cell score (SCS).

Methods

Data was collected for 67,206 primiparous cows from the database of the Paraná Holstein Breeders Association in Brazil, with the aim of evaluating the temperature effect, considered as an environmental variable, distinguished under six gradients, with the variation range found being 17°C to 19.5°C, over the region. A reaction norm model was adopted utilizing the fourth order under the Legendre polynomials, using the mixed models of analysis by the restricted maximum likelihood method by the WOMBAT software. Additionally, the genetic behavior of the 15 most representative bulls was assessed, in response to the changes in the temperature gradient.

Results

A mean score of 2.66 and a heritability variation from 0.17 to 0.23 was found in the regional temperature increase. The correlation between the environmental gradients proved to be higher than 0.80. Distinctive genetic behaviors were observed according to the increase in regional temperature, with an observed increase of up to 0.258 in the breeding values of some animals, as well as a reduction in the breeding of up to 0.793, with occasional reclassifications being observed as the temperature increased.

Conclusion

Non-relevant GEI for SCS were observed in Holstein cattle herds of southern Brazil. Thus, the inclusion of the temperature effect in the model of genetic evaluation of SCS for the southern Brazilian Holstein breed is not required.

Population genetic structure and maternal lineage of South African crossbred Nguni cattle using the cytochrome b gene in mtDNA

The Nguni cattle breed predominates South Africa and is endowed with traits favourable against environmental stressors such as heat stress and resistance to diseases. Interventions to improve production have led to the erosion of the genetic integrity of local breeds and the introduction of exotic breeds has proved ineffective as they fail to perform well due to different climatic conditions and production systems. In this study, the genetic structure and genetic lineage of Nguni crossbreds from 6 populations were assessed using the mitochondrial cytochrome b gene. Twelve polymorphic sites were detected resulting in 11 haplotypes with haplotype and nucleotide diversities of 0.550 ± 0.135 and 0.0019 ± 0.0011, respectively. Only 2 of the 6 populations displayed recent population expansion events, whereas the majority adhered to neutral evolution. The basal haplotype contained approximately 60% of the studied populations and there were four unique haplotypes that were revealed. A possible Nguni descript haplotype was uncovered, and this haplotype was found in all populations but was however devoid of individuals from around the world. The genetic structure of the populations was rather low (average pairwise F_ST = 0.066 and Slatkins F_ST = 0.094), and approximately 96% of the total genetic variation was accounted for by differences within populations. Phylogenetic analyses supported the clustering of all the samples within the Bos taurus clade and no Bos indicus haplotype was detected. Furthermore, no intermediate haplotype of taurine and indicine was detected. Overall, the maternal lineage of the crossbreds points to a taurine origin and the low genetic diversity depicts the retention of the Nguni genetic pool and possibly its superior adaptive traits.

Characterization of Holstein and Normande whole milk miRNomes highlights breed specificities

The concept of milk as a healthy food has opened the way for studies on milk components, from nutrients to microRNAs, molecules with broad regulatory properties present in large quantities in milk. Characterization of these components has been performed in several species, such as humans and bovine, depending on the stages of lactation. Here, we have studied the variation in milk microRNA composition according to genetic background. Using high throughput sequencing, we have characterized and compared the milk miRNomes of Holstein and Normande cattle, dairy breeds with distinct milk production features, in order to highlight microRNAs that are essential for regulation of the lactation process. In Holstein and Normande milk, 2,038 and 2,030 microRNAs were identified, respectively, with 1,771 common microRNAs, of which 1,049 were annotated and 722 were predicted. The comparison of the milk miRNomes of two breeds allowed to highlight 182 microRNAs displaying significant differences in the abundance. They are involved in the regulation of lipid metabolism and mammary morphogenesis and development, which affects lactation. Our results provide new insights into the regulation of molecular mechanisms involved in milk production.

Genetic and Epigenetic Regulation of Immune Response and Resistance to Infectious Diseases in Domestic Ruminants

Infectious diseases are the outcome of complex interactions between the host, pathogen, and environment. After exposure to a pathogen, the host immune system uses various mechanisms to remove the pathogen. However, environmental factors and characteristics of pathogens can compromise the host immune responses and subsequently alter the outcome of infection. In this article, genetic and epigenetic factors that shape the individual variation in mounting protective responses are reviewed. Different approaches that have been used by researchers to investigate the genetic regulation of immunity in ruminants and various sources of genetic information are discussed.

Genetic correlations of fighting ability with somatic cells and longevity in cattle

The success in competitions may be stressful for animals and costly in terms of immune functions and longevity. Focusing on Aosta Chestnut and Aosta Black Pied cattle, selected for their fighting ability in traditional competitions, this study investigated the genetic relationships of fighting ability with udder health traits (somatic cell score and two threshold traits for somatic cells), longevity (length of productive life and number of calvings) and test-day milk, fat and protein yield. Herdbook information and phenotypic records that have been routinely collected for breeding programs in 16 years were used for the abovementioned traits. Data belonged to 9328 cows and 19 283 animals in pedigree. Single-trait animal model analyses were run using a Gibbs sampling algorithm to estimate the variance components of traits, and bivariate analyses were then performed to estimate the genetic correlations. Moderate positive genetic correlations (ra) were found for fighting ability with somatic cell score (ra=0.255), suggesting that greater fighting ability is genetically related to a detriment in udder health, in agreement with the theory. The high positive genetic correlation between fighting ability and longevity (average ra=0.669) suggests that the economic importance of fighting ability (the winning cows get an higher price at selling) had probably masked the true genetic covariances. The genetic correlation between milk yield traits and fighting ability showed large intervals, but the negative values (average ra=-0.121) agreed with previous research. This study is one of the few empirical studies on genetic correlations for the competitive success v. immune functions and longevity traits. The knowledge of the genetic correlations among productive and functional traits of interest, including fighting ability, is important in animal breeding for a sustainable genetic improvement.

Genetic parameters for mastitis incidence and its indicators based on somatic cell score for Holsteins in Hokkaido, Japan

The objectives of this study were to estimate the heritability of mastitis incidence and genetic correlations between the mastitis and the somatic cell score (SCS) statistics, and to compare the practicability between different models. We used test-day records with the mastitis incidence and SCS collected from Holstein cows calving from 1988 to 2015 in Hokkaido, Japan. As indicators of mastitis, the average SCS (avSCS), the standard deviation of SCS (sdSCS), and the maximum SCS (maxSCS) were calculated using test-day records up to the first 305 days in milk within a lactation. We compared a four-trait repeatability animal model (MTRP) with a four-trait multiple-lactation animal model (MTML). The heritability for mastitis was equal to or lower than 0.05 in all the models. Genetic correlations between lactations with MTML within the same trait were positive and close to 1. With MTRP, the estimated genetic correlations of the mastitis incidence with avSCS, sdSCS, and maxSCS were 0.66, 0.79, and 0.82, respectively. A joint evaluation with SCS statistics is expected to give an extra reliability for mastitis because of high and positive genetic correlations among the traits.

A longitudinal prospective cohort study investigating the association of premilking stimulation and teat-end shape on milking characteristics and teat tissue condition in dairy cows

Background

Premilking udder preparation is essential for harvesting high-quality milk as gently, completely, and quickly as possible. The associations between characteristics such as teat-end shape and premilking stimulation on milking characteristics and machine milking-induced changes to the teat tissue condition have not been rigorously investigated. The primary objective was to investigate the interactive effects of manual premilking stimulation (i.e., preparation lag time) and teat-end shape on total milk yield, two-minute milk yield, milking unit-on time, and time in low milk flow rate. Our secondary objective was to study the association of manual premilking stimulation and changes to the teat tissue condition after machine milking (i.e., short-term changes). In a longitudinal prospective cohort study, 384 milking observations from 129 cows were analysed. Holstein cows were housed in sand-bedded free-stall pens, fed a total mixed ration, and milked 3 times a day. Cows were classified by teat-end shape into 1 of 3 categories: pointed, flat, or round. Individual cow milking characteristics were recorded with electronic on-farm milk meters. The duration of manual stimulation, preparation lag time, and presence of short-term changes were documented for each milking observation. General linear mixed models were used to study the interactive effects of preparation lag time and teat-end shape on milking characteristics.

Results

There was an interaction between preparation lag time and teat-end shape for two-minute milk yield and time in low milk flow rate. The preparation lag time effect was modified by teat-end shape, while no interaction was observed for total milk yield or milking unit-on time. A generalized linear mixed model revealed that preparation lag time was associated with short-term changes in teat tissue condition, where the odds of short-term changes decreased as preparation lag time increased.

Conclusions

In summary, cows with different teat-end shapes may require different premilking stimulation regimens. Increasing preparation lag time benefits teat tissue condition during machine milking. Further research is warranted to optimize individual premilking stimulation in dairy cows.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-1803-2) contains supplementary material, which is available to authorized users.

Detection of quantitative trait loci and putative causal variants affecting somatic cell score in dairy sheep by using a 50K SNP chip and whole-genome sequencing

This study presents a scan of the ovine genome to identify quantitative trait loci (QTL) influencing the somatic cell score (SCS), a classical indicator of subclinical mastitis in sheep, and a subsequent high-resolution analysis of one of the identified QTL regions based on the analysis of whole-genome sequence data sets. A half-sib commercial population of Churra sheep genotyped with a 50K SNP chip was analyzed using linkage analysis (LA) and combined linkage and linkage disequilibrium analysis (LDLA). By LA, 2 5% chromosome-wide significant QTL on OAR5 and OAR25 and one 5% genome-wide significant QTL on ovine chromosome 20 (OAR20) were detected, whereas 22 significant associations were identified by LDLA. Two of the associations detected by LDLA replicated LA-detected effects (OAR20, OAR25). We compared the detected associations with previously reported QTL in sheep and cattle, and functional candidate genes were identified within the estimated confidence intervals. We then performed a high-resolution analysis of the OAR20 QTL region, the most significant QTL region identified by LA that replicated a QTL previously described in Churra sheep for SCS using microsatellite markers. For that, 2 segregating trios of 2 segregating families for the OAR20 QTL (each including the Qq sire and 2 daughters, QQ and qq) were selected for whole-genome sequencing. The bioinformatic analysis of the 6 sequenced samples performed across the genomic interval considered (14.2-41.7 Mb) identified a total of 227,030 variants commonly identified by 2 independent software packages. For the 3 different concordance tests considered, due to discrepancies regarding the QTL peak in the segregating families, the list of mutations concordant with the QTL segregating pattern was processed to identify the variants identified in immune-related genes that show a moderate/high impact on the encoded protein function. Among a list of 85 missense variants concordant with the QTL segregation pattern that were within candidate immune-related genes, 13 variants distributed across 7 genes [PKHD1, NOTCH4, AGER, ENSOARG00000009395 (HLA-C, Homo sapiens), ENSOARG00000015002 (HLA-B, H. sapiens), MOG, and ENSOARG00000018075 (BoLA, Bos taurus, orthologous to human HLA-A] were predicted to cause deleterious effects on protein function. Future studies should assess the possible associations of the candidate variants identified herein in commercial populations with indicator traits of udder inflammation (SCS, clinical mastitis).

Invited review: Low milk somatic cell count and susceptibility to mastitis

An enduring controversy exists about low milk cell counts and susceptibility to mastitis. The concentration of milk leukocytes, or somatic cell count (SCC), is a well-established direct indicator of mammary gland inflammation that is highly correlated with the presence of a mammary infection. The SCC is also used as a trait for the selection of dairy ruminants less prone to mastitis. As selection programs favor animals with less SCC, and as milk cells contribute to the defense of the mammary gland, the idea that susceptibility to mastitis could possibly be increased in the long term has been put forward and is still widely debated. Epidemiological and experimental studies aimed at relating SCC to susceptibility to mastitis have yielded results that seem contradictory at first sight. Nevertheless, by taking into account the immunobiology of milk and mammary tissue cells and their role in the defense against infection, along with recent studies on SCC-based divergent selection of animals, the issue can be settled. Apparent SCC-linked susceptibility to mastitis is a phenotypic trait that may be linked to immunomodulation but not to selection.

Genome-wide association study for milking speed in French Holstein cows

Using a combination of data from the BovineSNP50 BeadChip SNP array (Illumina, San Diego, CA) and a EuroGenomics (Amsterdam, the Netherlands) custom single nucleotide polymorphism (SNP) chip with SNP pre-selected from whole genome sequence data, we carried out an association study of milking speed in 32,491 French Holstein dairy cows. Milking speed was measured by a score given by the farmer. Phenotypes were yield deviations as obtained from the French evaluation system. They were analyzed with a linear mixed model for association studies. We identified SNP on 22 chromosomes significantly associated with milking speed. As clinical mastitis and somatic cell score have an unfavorable genetic correlation with milking speed, we tested whether the most significant SNP on these 22 chromosomes associated with milking speed were also associated with clinical mastitis or somatic cell score. Nine hundred seventy-one genome-wide significant SNP were associated with milking speed. Of these, 86 were associated with clinical mastitis and 198 with somatic cell score. The most significant association signals for milking speed were observed on chromosomes 7, 8, 10, 14, and 18. The most significant signal was located on chromosome 14 (ZFAT gene). Eleven novel milking speed quantitative trait loci (QTL) were observed on chromosomes 7, 10, 11, 14, 18, 25, and 26. Twelve candidate SNP for milking speed mapped directly within genes. Of these, 10 were QTL lead SNP, which mapped within the genes HMHA1, POLR2E, GNB5, KLHL29, ZFAT, KCNB2, CEACAM18, CCL24, and LHPP. Limited pleiotropy was observed between milking speed QTL and clinical mastitis.

Symposium review: Novel strategies to genetically improve mastitis resistance in dairy cattle

Mastitis is a disease of major economic importance to the dairy cattle sector because of the high incidence of clinical mastitis and prevalence of subclinical mastitis and, consequently, the costs associated with treatment, production losses, and reduced animal welfare. Disease-recording systems compiling data from a large number of farms are still not widely implemented around the world; thus, selection for mastitis resistance is often based on genetically correlated indicator traits such as somatic cell count (SCC), udder depth, and fore udder attachment. However, in the past years, several countries have initiated collection systems of clinical mastitis, based on producers recording data in most cases. The large data sets generated have enabled researchers to assess incidence of this disease and to investigate the genetic background of clinical mastitis itself, as well as its relationships with other traits of interest to the dairy industry. The genetic correlations between clinical mastitis and its previous proxies were estimated more accurately and confirmed the strong relationship of clinical mastitis with SCC and udder depth. New traits deriving from SCC were also studied, with the most relevant findings being associated with mean somatic cell score (SCS) in early lactation, standard deviation of SCS, and excessive test-day SCC pattern. Genetic correlations between clinical mastitis and other economically important traits indicated that selection for mastitis resistance would also improve resistance against other diseases and enhance both fertility and longevity. However, milk yield remains negatively correlated with clinical mastitis, emphasizing the importance of including health traits in the breeding objectives to achieve genetic progress for all important traits. These studies enabled the establishment of new genetic and genomic evaluation models, which are more efficient for selection to mastitis resistance. Further studies that are potential keys for future improvement of mastitis resistance are deep investigation of the bacteriology of mastitis, identification of novel indicator traits and tools for selection, and development of a larger female reference population to improve reliability of genomic evaluations. These cutting-edge studies will result in a better understanding of the genetic background of mastitis resistance and enable a more accurate phenotyping and genetic selection to improve mastitis resistance, and consequently, animal welfare and industry profitability.

Genome-wide methylation analysis reveals differentially methylated loci that are associated with an age-dependent increase in bovine fibroblast response to LPS

BACKGROUND

Differences in DNA methylation are known to contribute to the development of immune-related disorders in humans but relatively little is known about how methylation regulates immune function in cattle. Utilizing whole-transcriptome analyses of bovine dermal fibroblasts, we have previously identified an age and breed-dependent up-regulation of genes within the toll-like receptor 4 (TLR4) pathway that correlates with enhanced fibroblast production of IL-8 in response to lipopolysaccharide (LPS). Age-dependent differences in IL-8 production are abolished by treatment with 5-aza-2-deoxycytidine and Trichostatin A (AZA-TSA), suggesting epigenetic regulation of the innate response to LPS. In the current study, we performed reduced representation bisulfite sequencing (RRBS) on fibroblast cultures isolated from the same animals at 5- and 16-months of age to identify genes that exhibit variable methylation with age. To validate the role of methylation in gene expression, six innate response genes that were hyper-methylated in young animals were assessed by RT-qPCR in fibroblasts from animals at different ages and from different breeds.

RESULTS

We identified 14,094 differentially methylated CpGs (DMCs) that differed between fibroblast cultures at 5- versus 16-months of age. Of the 5065 DMCs that fell within gene regions, 1117 were located within promoters, 1057 were within gene exons and 2891 were within gene introns and 67% were more methylated in young cultures. Transcription factor enrichment of the promoter regions hyper-methylated in young cultures revealed significant regulation by the key pro-inflammatory regulator, NF-κB. Additionally, five out of six chosen genes (PIK3R1, FES, NFATC1, TNFSF13 and RORA) that were more methylated in young cultures showed a significant reduction in expression post-LPS treatment in comparison with older cultures. Two of these genes, FES and NFATC1, were similarly down-regulated in Angus cultures that also exhibit a low LPS response phenotype.

CONCLUSIONS

Our study has identified immune-related loci regulated by DNA methylation in cattle that may contribute to differential cellular response to LPS, two of which exhibit an identical expression profile in both low-responding age and breed phenotypes. Methylation biomarkers of differential immunity may prove useful in developing selection strategies for replacement cows that are less susceptible to severe infections, such as coliform mastitis.

Bivariate threshold models for genetic evaluation of susceptibility to and ability to recover from mastitis in Danish Holstein cows

Mastitis in dairy cows is an unavoidable problem and genetic variation in recovery from mastitis, in addition to susceptibility, is therefore of interest. Genetic parameters for susceptibility to and recovery from mastitis were estimated for Danish Holstein-Friesian cows using data from automatic milking systems equipped with online somatic cell count measuring units. The somatic cell count measurements were converted to elevated mastitis risk, a continuous variable [on a (0-1) scale] indicating the risk of mastitis. Risk values >0.6 were assumed to indicate that a cow had mastitis. For each cow and lactation, the sequence of health states (mastitic or healthy) was converted to a weekly transition: 0 if the cow stayed within the same state and 1 if the cow changed state. The result was 2 series of transitions: one for healthy to diseased (HD, to model mastitis susceptibility) and the other for diseased to healthy (DH, to model recovery ability). The 2 series of transitions were analyzed with bivariate threshold models, including several systematic effects and a function of time. The model included effects of herd, parity, herd-test-week, permanent environment (to account for the repetitive nature of transition records from a cow) plus two time-varying effects (lactation stage and time within episode). In early lactation, there was an increased risk of getting mastitis but the risk remained stable afterwards. Mean recovery rate was 45% per lactation. Heritabilities were 0.07 [posterior mean of standard deviations (PSD) = 0.03] for HD and 0.08 (PSD = 0.03) for DH. The genetic correlation between HD and DH has a posterior mean of -0.83 (PSD = 0.13). Although susceptibility and recovery from mastitis are strongly negatively correlated, recovery can be considered as a new trait for selection.

Genetic trends for fertility, udder health and protein yield in Swedish red cattle estimated with different models

The aim of this study was to estimate and compare genetic trends in Swedish Red cattle using a full multiple-trait (MT) model and trait-group-wise models for female fertility, udder health and protein yield. Field data for maiden heifers from 1989 and cows with a first and second lactation between 1990 and 2007 were included. (Co)variance components were estimated prior to prediction of breeding values. The estimated genetic trends were clearly favourable for protein yield and udder conformation, and in most cases neutral to favourable for clinical mastitis and calving to first insemination. In maiden heifers, the trends were neutral for number of inseminations per service period. Unfavourable genetic trends were estimated for number of inseminations in the first two lactations, but the trends seemed less unfavourable from evaluations within trait groups compared with when using the full MT model. Excluding maiden heifer data affected genetic trends less than using trait-group-wise analyses instead of a full MT model. Unfavourable genetic trends in functional traits may be missed unless the traits are evaluated in a MT model including traits under strong selection.

Differences during the first lactation between cows cloned by somatic cell nuclear transfer and noncloned cows

Lactation performance is dependent on both the genetic characteristics and the environmental conditions surrounding lactating cows. However, individual variations can still be observed within a given breed under similar environmental conditions. The role of the environment between birth and lactation could be better appreciated in cloned cows, which are presumed to be genetically identical, but differences in lactation performance between cloned and noncloned cows first need to be clearly evaluated. Conflicting results have been described in the literature, so our aim was to clarify this situation. Nine cloned Prim' Holstein cows were produced by the transfer of nuclei from a single fibroblast cell line after cell fusion with enucleated oocytes. The cloned cows and 9 noncloned counterparts were raised under similar conditions. Milk production and composition were recorded monthly from calving until 200d in milk. At 67d in milk, biopsies were sampled from the rear quarter of the udder, their mammary epithelial cell content was evaluated, and mammary cell renewal, RNA, and DNA were then analyzed in relevant samples. The results showed that milk production did not differ significantly between cloned and noncloned cows, but milk protein and fat contents were less variable in cloned cows. Furthermore, milk fat yield and contents were lower in cloned cows during early lactation. At around 67 DIM, milk fat and protein yields, as well as milk fat, protein, and lactose contents, were also lower in cloned cows. These lower yields could be linked to the higher apoptotic rate observed in cloned cows. Apoptosis is triggered by insulin-like factor growth binding protein 5 (IGFBP5) and plasminogen activator inhibitor (PAI), which both interact with CSN1S2. During our experiments, CSN1S2 transcript levels were lower in the mammary gland of cloned cows. The mammary cell apoptotic rate observed in cloned cows may have been related to the higher levels of DNA (cytosine-5-)-methyltransferase 1 (DNMT1) transcripts, coding for products that maintain the epigenetic status of cells. We conclude, therefore, that milk production in cloned cows differs slightly from that of noncloned cows. These differences may be due, in part, to a higher incidence of subclinical mastitis. They were associated with differences in cell apoptosis and linked to variations in DNMT1 mRNA. However, milk protein and fat contents were more similar among cloned cows than among noncloned cows.