Genetic polymorphisms of TRAPPC9 and CD4 genes and their association with milk production and mastitis resistance phenotypic traits in Chinese Holstein

The present study was designed to evaluate the association of polymorphisms in bovine trafficking protein particle complex subunit 9 (TRAPPC9) and cluster of differentiation 4 (CD4) genes with milk production and mastitis resistance phenotypic traits in a different cattle population. Three single nucleotide polymorphisms (SNPs) (SNP1 Position: Chr14:2484891, SNP2 (rs110017379), SNP3 Position: Chr14:2525852) in bovine TRAPPC9 and one SNP (Position: Chr5:104010752) in CD4 were screened through Chinese Cow's SNPs Chip-I (CCSC-I) and genotyped in a population of 312 Chinese Holsteins (156: Mastitis, 156: Healthy). The results were analyzed using the general linear model in SAS 9.4. Our analysis revealed that milk protein percentage, somatic cell count (SCC), somatic cell score (SCS), serum cytokines interleukin 6 (IL-6) and interferon-gamma (IFN-γ) were significantly (P < 0.05) associated with at least one or more identified SNPs of TRAPPC9 and CD4 genes. Furthermore, the expression status of SNPs in CD4 and TRAPPC9 genes were verified through RT-qPCR. The expression analysis showed that genotypes GG in SNP3 of TRAPPC9 and TT genotype in SNP4 of CD4 showed higher expression level compared to other genotypes. The GG genotype in SNP2 and TT genotype in SNP3 of TRAPPC9 were associated with higher bovine milk SCC and lower IL6. Altogether, our findings suggested that the SNPs of TRAPPC9 and CD4 genes could be useful genetic markers in selection for milk protein improvement and mastitis resistance phenotypic traits in dairy cattle. The CCSC-I used in current study is proposed to be validate in different and large population of dairy cattle not only in China but also in other countries. Moreover, our analyses recommended that besides SCC and SCS, the association of genetic markers could also be considered with the serum cytokines (IL-6, IFN-γ) while selecting genetically mastitis resistance dairy cattle.

Polymorphisms of the IL-17A Gene Influence Milk Production Traits and Somatic Cell Score in Chinese Holstein Cows

The cow’s milk production characteristics are a significant economic indicator in the livestock industry. Serum cytokines such as interleukin-17 (IL-17) may be potential indicators for bovine mastitis concerning the milk somatic cell count (SCC) and somatic cell score (SCS). The current study aims to find previously undiscovered single nucleotide polymorphisms in the bovine (IL-17A) gene and further investigates their associations with milk production traits in Chinese Holstein cows. Twenty Chinese Holstein cows were randomly chosen from six farms in Jiangsu Province, China. The DNA was extracted from selected samples of bloods for PCR amplification Sequence analyses were used to find SNPs in the bovine (IL-17A) gene. The discovered five SNPs are g-1578A>G, g-1835G>A, and g-398T>A in the 5′UTR; g3164T>C and g3409G>C in the exon region. The genotyping of Holstein cows (n = 992) was performed based on Sequenom Mass ARRAY and SNP data. The connection between SNPs, milk production variables, and the somatic cell score was investigated using the least-squares method. Based on the results, SNP g-398T>A had a significant linkage disequilibrium with g3164T>C. SNPs were found to have significant (p < 0.05) correlations with the test-day milk yield. In conclusion, IL-17A affects cow’s milk production traits significantly.

Genetic parameters and association of national evaluations with breeding values for health traits in US organic Holstein cows

Among other regulations, organic cows in the United States cannot receive antibiotics and preserve their organic status, emphasizing the importance of prevention of illness and benefit of high genetic merit for disease resistance. At the same time, data underlying national genetic evaluations primarily come from conventional cows, drawing concern to the possibility of a genotype by environment interaction whereby the value of a genotype varies depending on the environment, and potentially limits the relevance of these evaluations to organic cows. The objectives of this study were to characterize the genetics of and determine the presence of genotype by environment interaction for health traits in US organic dairy cows. Individual cow health data were obtained from 16 US Department of Agriculture certified organic dairy farms from across the United States that used artificial insemination and maintained detailed records. Data were obtained for the following traits: died, lameness, mastitis, metabolic diseases (displaced abomasum, ketosis, and milk fever), reproductive diseases (abortion, metritis, and retained placenta), transition health events (any health event occurring 21 d before or after parturition), and all health events. Binary phenotypes (1 = diseased, 0 = otherwise) for 38,949 lactations on 19,139 Holstein cows were used. Genotypes from 2,347 cows with 87.5% or greater Holstein breed-based representation were incorporated into single-step multitrait threshold animal models that included stayability (1 = completed lactation, 0 = otherwise). Gibbs sampling was used. Genomic predicted transmitting abilities (gPTA) from national genetic evaluations were obtained for sires for production, fitness, health, and conformation traits. We approximated genetic correlations for sires using these gPTA and our estimated breeding values. We also regressed health phenotypes on cow estimated breeding values and sire gPTA. Heritabilities (± standard error) ranged from 0.03 ± 0.01 (reproductive diseases) to 0.11 ± 0.03 (metabolic diseases). Most genetic correlations among health traits were positive, though the genetic correlation between metabolic disease and mastitis was -0.42 ± 0.17. Approximate genetic correlations between disease resistance for our health trait categories and disease resistance for the nationally-evaluated health traits generally carried the expected sign with the strongest correlation for mastitis (0.72 ± 0.084). Regression coefficients carried the expected sign and were mostly different from zero, indicating that evaluations from primarily conventional herd data predicted health on organic farms. In conclusion, use of national evaluations for health traits should afford genetic improvement for health in US organic herds.

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.

Genetic analysis of resilience indicators based on milk yield records in different lactations and at different lactation stages

Resilience is the ability of cows to cope with disturbances, such as pathogens or heat waves. To breed for improved resilience, it is important to know whether resilience genetically changes throughout life. Therefore, the aim was to perform a genetic analysis on 2 resilience indicators based on data from 3 periods of the first lactation (d 11-110, 111-210, and 211-340) and the first 3 full lactations, and to estimate genetic correlations with health traits. The resilience indicators were the natural log-transformed variance (LnVar) and lag-1 autocorrelation (r_auto) of daily deviations in milk yield from an expected lactation curve. Low LnVar and r_auto indicate low variability in daily milk yield and quick recovery, and were expected to indicate good resilience. Data of 200,084 first, 155,784 second, and 89,990 third lactations were used. Heritabilities were similar based on different lactation periods (0.12-0.15 for LnVar, 0.05-0.06 for r_auto). However, the heritabilities of the resilience indicators based on full first lactation were higher than those based on lactation periods (0.20 for LnVar, 0.08 for r_auto), due to lower residual variances. Heritabilities decreased from 0.20 in full lactation 1 to 0.19 in full lactation 3 for LnVar and from 0.08 to 0.06 for r_auto. For LnVar, as well as for r_auto, the strongest genetic correlation between lactation periods was between period 2 and 3 (0.97 for LnVar, 0.96 for r_auto) and the weakest between period 1 and 3 (0.81 for LnVar, 0.65 for r_auto). Similarly, for both traits the genetic correlation between full lactations was strongest between lactations 2 and 3 (0.99 for LnVar, 0.95 for r_auto) and weakest between lactations 1 and 3 (0.91 for LnVar, 0.71 for r_auto). For LnVar, genetic correlations with resilience-related traits, such as udder health, ketosis, and longevity, adjusted for correlations with milk yield, were almost always favorable (-0.59 to 0.02). In most cases these genetic correlations were stronger based on full lactations than on lactation periods. Genetic correlations were similar across full lactations, but the correlation with udder health increased substantially from -0.31 in lactation 1 to -0.51 in lactation 3. For r_auto, genetic correlations with resilience-related traits were always favorable in lactation period 1 and in most full lactations, but not in the other lactation periods. However, correlations were weak (-0.27 to 0.15). Therefore, as a resilience indicator for breeding, LnVar is preferred over r_auto. A multitrait index based on estimated breeding values for LnVar in lactations 1, 2, and 3 is recommended to improve resilience throughout the lifetime of a cow.

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.

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.

Polymorphisms in the promoter region of the bovine lactoferrin gene influence milk somatic cell score and milk production traits in Chinese Holstein cows

Lactoferrin is an iron-binding protein found in cow's milk that plays an important role in preventing mastitis caused by intramammary infection. In this study, 20 Chinese Holstein cows were selected randomly for PCR amplification and sequencing of the bovine lactoferrin gene promoter region and used for SNP discovery in the region between nucleotide positions -461 to -132. Three SNPs (-270T>C, -190G>A and -156A>G) were identified in bovine lactoferrin, then Chinese Holstein cows (n=866) were genotyped using Sequenom MassARRAY (Sequenom Inc., San Diego, CA) based on the previous SNP information in this study, and the associations between SNPs or haplotype and milk somatic cell score (SCS) and production traits were analyzed by the least squares method in the GLM procedure of SAS. SNPs -270T>C and -156A>G showed close linkage disequilibrium (r(2)=0.76). The SNP -190G>A showed a significant association with SCS, and individuals with genotype GG had higher SCS than genotypes AG and AA. Associations were found between the SNPs -270T>C and -190G>A with SCS and the milk composition. The software MatInspector revealed that these SNPs were located within several potential transcription factor binding sites, including NF-κB p50, KLF7 and SP1, and may alter gene expression, but further investigation will be required to elucidate the biological and practical relevance of these SNPs.