Genetic relationship of discrete-time survival with fertility and production in dairy cattle using bivariate models

Bivariate analyses of functional longevity in dairy cattle measured as survival to next lactation (SURV) with milk yield and fertility traits were carried out. A sequential threshold-linear censored model was implemented for the analyses of SURV. Records on 96 642 lactations from 41 170 cows were used to estimate genetic parameters, using animal models, for longevity, 305 d-standardized milk production (MY305), days open (DO) and number of inseminations to conception (INS) in the Spanish Holstein population; 31% and 30% of lactations were censored for DO and INS, respectively. Heritability estimates for SURV and MY305 were 0.11 and 0.27 respectively; while heritability estimates for fertility traits were lower (0.07 for DO and 0.03 for INS). Antagonist genetic correlations were estimated between SURV and fertility (-0.78 and -0.54 for DO and INS, respectively) or production (-0.53 for MY305), suggesting reduced functional longevity with impaired fertility and increased milk production. Longer days open seems to affect survival more than increased INS. Also, high productive cows were more problematic, less functional and more liable to being culled. The results suggest that the sequential threshold model is a method that might be considered at evaluating genetic relationship between discrete-time survival and other traits, due to its flexibility.     

Genetic parameters for production,health, fertility and longevity traits in dairy cows

Milk production, fertility, longevity and health records, were extracted from databases of two milk recording organisations in the United Kingdom for the first three lactations of the Holstein–Friesian breed. These included data related to health events (mastitis and lameness), voluntarily recorded on a proportion of farms. The data were analysed to calculate disease incidence levels and to estimate genetic parameters for health traits and their relationships with production and other functional traits. The resulting dataset consisted of 124 793 lactations from 75 137 animals of 1586 sires, recorded in 2434 herds. Incidence of health events increased with parity. The overall incidence of mastitis (MAS) and lameness (LAM), defined as binary traits, were 17% and 16%, respectively. Heritability estimates for MAS and LAM were 0.04 and 0.02, respectively, obtained from repeatability linear sire models. Heritability estimates of mastitis and lameness as count traits were slightly higher, 0.05 and 0.03, respectively. Genetic correlations were obtained by bivariate analyses of all pair-wise combinations between milk 305-day yield (MY), protein 305-day yield (PY), fat 305-day yield (FY), lactation average log_e transformed lactation average somatic cell count (SCS), calving interval (CI), days to first service (DFS), non-return at 56 days (NR56), number of inseminations (NINS), mastitis (MAS), number of mastitis episodes (NMAS), lameness (LAM), number of lameness episodes (NLAM) and lifespan score (LS). As expected, MAS was correlated most strongly with SCS (0.69), which supports the use of SCS as an indicator trait for mastitis. Genetic correlations between MAS and yield and fertility traits were of similar magnitude ranging from 0.27 to 0.33. Genetic correlations between MAS with LAM and LS were 0.38 and −0.59, respectively. Not all genetic correlations between LAM and other traits were significant because of fewer numbers of lameness records. LAM had significant genetic correlations with MY (0.38), PY (0.28), CI (0.35), NINS (0.38) and LS (−0.53). The heritability estimates of mastitis and lameness were low; therefore, genetic gain through direct selection alone would be slow, yet still positive and cumulative. Direct selection against mastitis and lameness as additional traits should reduce incidence of both diseases, and simultaneously improve fertility and longevity. However, both health traits had antagonistic relationships with production traits, thus genetic gain in production would be slower.     

Genetic parameters for milk urea concentration and milk traits in Polish Holstein-Friesian cows

Milk urea concentration (MU) used by dairy producers for management purposes can be affected by selection for milk traits. To assess this problem, genetic parameters for MU in Polish Holstein-Friesian cattle were estimated for the first three lactations. The genetic correlation of MU with milk production traits, lactose percentage, fat to protein ratio (FPR) and somatic cell score (SCS) were computed with two 5-trait random regression test-day models, separately for each lactation. Data used for estimation (159,044 daily observations) came from 50 randomly sampled herds. (Co)variance components were estimated with the Bayesian Gibbs sampling method. The coefficient of variation for MU in all three parities was high (40–41 %). Average daily heritabilities of MU were 0.22 for the first parity and 0.21 for the second and third lactations. Average genetic correlations for different days in milk in the first three lactations between MU and other traits varied. They were small and negative for protein percentage (from ?0.24 to ?0.11) and for SCS (from ?0.14 to ?0.09). The weakest genetic correlation between MU and fat percentage, and between MU and lactose percentage were observed (from ?0.10 to 0.10). Negative average genetic correlation with the fat to protein ratio was observed only in the first lactation (?0.14). Genetic correlations with yield traits were positive and ranged from low to moderate for protein (from 0.09 to 0.33), fat (from 0.16 to 0.35) and milk yield (from 0.20 to 0.42). These results suggest that the selection on yield traits and SCS tends to increase MU slightly.     

Genotype by environment interaction for Holstein cattle populations using autoregressive and within- and across-country multi-trait reaction norms test-day models

The progenies of international bulls in diverse climatic conditions and management levels may lead to different expressions of their genetic potential resulting in a re-ranking of these bulls. Therefore, evaluate the presence of genotype by environment interaction (G × E) within and across countries is important to guide the decision-making on alternative selection strategies. Thus, a two-step reaction norm (RN) approach was used to investigate the presence of G × E in Portuguese and Brazilian Holstein cattle. In step 1, we performed a within-country genetic evaluation using an autoregressive model to obtain precorrected phenotypes and environmental gradients (herd test-day solutions, HTD levels). In step 2, the precorrected phenotypes were considered as two distinct traits in a bi-trait RN model to estimate variance components across HTD levels, genetic correlation between HTD levels in Portugal and Brazil, and RN of the estimated breeding values. Additionally, the genetic correlation between countries using a bi-trait random regression (RR) sire model was obtained. In step 1, genetic additive variance for milk yield (MY) in Portugal was 14.1% higher than in Brazil. For somatic cell score (SCS), the genetic additive variance in Portugal was 12.7% lower than in Brazil. Although similar heritability estimates for SCS were observed in both countries, MY heritabilities were 0.31 for Portugal and 0.23 for Brazil. Genetic correlations (SD) between both countries obtained using RR sire model were 0.78 (0.051) for MY and 0.75 (0.062) for SCS. In step 2, MY genetic correlations among HTD levels within countries were higher than 0.94 for Portugal and 0.98 for Brazil. Somatic cell score genetic correlations among HTD levels ranged from 0.70 to 0.99 for Portugal and from 0.84 to 0.99 for Brazil. The average (SD) of genetic correlation estimates between Portuguese and Brazilian HTD levels were 0.74 (0.009) for MY and 0.57 (0.060) for SCS. These results suggest the presence of G × E for MY and SCS of Holstein cattle between both countries. Although there was no indication of G × E between Brazilian herd environments, the low genetic correlation for SCS indicates potential re-ranking of bulls between extreme environmental gradient in Portugal. Overall, the results of this study evidence the importance of national and international genetic evaluation systems to assist dairy farmers in the selection of the best genotypes to obtain the expected returns from investments in imported semen and to realize genetic progress in dairy populations under local environmental conditions.     

Genetic correlations among milk yield,morphology, performance test traits and somatic cells in dual-purpose Rendena breed

Selection in native local breeds needs great carefulness due to the small population size and the risk of inbreeding. Furthermore, most breeds are dual-purpose, and milk and beef attitudes are antagonistic. For preservation purposes functional traits need to be considered. Focusing on the small local Rendena cattle, this study aimed to analyse the genetic correlations among milk, beef and udder health traits and the response to selection predicted under different scenarios. The study considered milk, fat and protein yields (MY), factor scores for udder volume (UV), conformation (UC) and muscularity obtained from type traits scored on primiparous cows, and performance test traits (PT) measured on young bulls at test station: average daily gain, in vivo SEUROP fleshiness, in vivo dressing percentage. Somatic cell score (SCS) was considered as a functional trait, with a possibility of restricting its genetic gain to zero. The study considered 281 497 MY test-day data collected on 16 974 cows, and data from linear type evaluation on 11 992 primiparous cows for factor scores. The PT data were recorded on 1428 young bulls, and SCS obtained from cell counts at milk recording. Bi-trait restricted maximum likelihood animal model analyses were performed to assess genetic parameters. Heritability varied from 0.157 (fat) to 0.442 (dressing percentage). Udder volume and MY resulted positively genetically correlated (average correlation 0.427), whereas the low-negative genetic correlation between MY and UC (−0.141) suggested a negative impact of milk gain on udder form. Beef traits of factor muscularity and PT showed medium-high favourable genetic correlations (from 0.357 to 0.984), excluding a null correlation between daily gain and muscularity. The genetic correlation MY v. muscularity was unfavourable (−0.328 on average), whereas null correlations were found in MY v. PT, apart from fat v. dressing percentage (−0.151). Somatic cell score showed low unfavourable correlations with protein (0.111) and UV (0.092), and favourable correlations with UC (−0.193). Response to selection in different scenarios indicated a good balanced gain for milk and beef when standardized economic weights of 0.66 and 0.34 are given to the two attitudes, and SCS genetic gain is restricted. Current genetic trends (MY and PT increasing, but muscularity lessening) reflect a stronger selection for milk, suggesting a slight progressive change towards a milk conformation. Aiming to preserve the dual-purpose characteristics of a breed, proper breeding policies taking into account the genetic relationships among traits and including functional traits should be applied in local dual-purpose populations.     

Analyses of the genetic relationships between lactose,somatic cell score,and growth traits in Simmental cattle

Lactose and somatic cell score (SCS) are major economic traits of milk. However, for many countries, they are typically not directly considered in the national genetic evaluation of Simmental cattle. This study aimed to estimate the genetic relationships between lactose, SCS, and growth traits of Simmental cattle to provide information for the national genetic evaluation of the selection of traits of this cattle population. The data of 1781 animals with 6519 records obtained over a period of 41 years (1975–2016) were collected from Xinjiang Hutubi Farm, China. The analyzed traits included 305 days of milk yield (305MY), milk fat percentage (MFP), milk protein percentage (MPP), milk lactose percentage (MLP), total solids (TS), SCS, body height (BH), body length (BL), chest girth (CG), abdominal circumference (AC), rump width (RW), rump length (RL), leg circumference (LC), and cannon circumference (CC). The multiple-trait repeatability model was adopted to estimate (co)variance components using the average information-restricted maximum likelihood method implemented using the DMU statistical package. The heritability estimates for milk components and growth traits ranged from 0.09 (SCS) to 0.51 (BH). Genetic correlations for milk components ranged from 0.03 ± 0.14 (MFP and MLP) to 0.81 ± 0.08 (MFP and MPP). Genetic correlation between MLP and SCS was moderate and negative (− 0.50 ± 0.15) compared with that among other traits. Genetic correlations between the milk components and growth traits ranged from 0.00 ± 0.07 (305MY and RW) to − 0.64 ± 0.15 (MLP and BL). Genetic correlations of BL, LC, RW, and RL with MLP were moderate to high and negative ranging from − 0.39 to − 0.64. Somatic cell score showed the highest correlation with BL (0.41) followed by LC (0.21). An increase in MLP would result in an increase in 305MY or TS and a decrease in BL, LC, RW, and RL. Additionally, a decrease in SCS would occur with the selection of increased MLP and reduced BL. We conclude that selection based on easily and inexpensively measured growth traits could improve the milk quality from Simmental cattle.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Genes influencing milk production traits predominantly affect one of four biological pathways

In this study we introduce a method that accounts for false positive and false negative results in attempting to estimate the true proportion of quantitative trait loci that affect two different traits. This method was applied to data from a genome scan that was used to detect QTL for three independent milk production traits, Australian Selection Index (ASI), protein percentage (P%) and fat percentage corrected for protein percentage (F% – P%). These four different scenarios are attributed to four biological pathways: QTL that (1) increase or decrease total mammary gland production (affecting ASI only); (2) increase or decrease lactose synthesis resulting in the volume of milk being changed but without a change in protein or fat yield (affecting P% only); (3) increase or decrease protein synthesis while milk volume remains relatively constant (affecting ASI and P% in the same direction); (4) increase or decrease fat synthesis while the volume of milk remains relatively constant (affecting F% – P% only). The results indicate that of the positions that detected a gene, most affected one trait and not the others, though a small proportion (2.8%) affected ASI and P% in the same direction.     

Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) with production traits in Canadian dairy cattle

Associations of two alleles of the bovine major histocompatibility complex DRB3 gene (BoLA-DRB3) with lowered somatic cell score (SCS) and occurrence of disease (BoLA-DRB3.2* 16 and *23, respectively) have previously been documented. The objective of this study was to evaluate potential relationships between BoLA-DRB3 alleles with production traits, namely 305-day milk, milk fat and milk protein yield, in a population of Canadian dairy cattle (Holstein, n = 835 and Jersey, n = 66) over the course of two lactations. No significant associations were detected between BoLA alleles and production traits in Jerseys. In Holsteins, alleles *16 and *23 also did not show associations with production traits but allele *8 was significantly associated with increased 305-day milk, fat and protein yields in the previous lactation (the lactation prior to immunization with a gram negative core antigen vaccine), and with increased protein production in the subsequent (with reference to the time of immunization) lactation. Allele *22 was associated with decreased milk and protein yield in both previous and subsequent lactations. Therefore, it can be concluded that increasing or decreasing the frequency of BoLA alleles *16 and *23 to reduce SCS or increase resistance to mastitis in this population would not have adverse effects on production in this population, and that certain BoLA alleles (*8 and *22) are associated with altered production traits in Canadian Holsteins.     

Bayesian segregation analysis of milk flow in Swiss dairy cattle using Gibbs sampling

Segregation analyses with Gibbs sampling were applied to investigate the mode of inheritance and to estimate the genetic parameters of milk flow of Swiss dairy cattle. The data consisted of 204 397, 655 989 and 40 242 lactation records of milk flow in Brown Swiss, Simmental and Holstein cattle, respectively (4 to 22 years). Separate genetic analyses of first and multiple lactations were carried out for each breed. The results show that genetic parameters especially polygenic variance and heritability of milk flow in the first lactation were very similar under both mixed inheritance (polygenes + major gene) and polygenic models. Segregation analyses yielded very low major gene variances which favour the polygenic determinism of milk flow. Heritabilities and repeatabilities of milk flow in both Brown Swiss and Simmental were high (0.44 to 0.48 and 0.54 to 0.59, respectively). The heritability of milk flow based on scores of milking ability in Holstein was intermediate (0.25). Variance components and heritabilities in the first lactation were slightly larger than those estimates for multiple lactations. The results suggest that milk flow (the quantity of milk per minute of milking) is a relevant measurement to characterise the cows milking ability which is a good candidate trait to be evaluated for a possible inclusion in the selection objectives in dairy cattle.     

Genome‐wide association studies for milk production traits in Valle del Belice sheep using repeated measures

Genome‐wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test‐day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.     

荷斯坦牛Nramp1基因遗传多态性及其与乳房炎相关性的研究

利用PCR-SSCP技术检测了344头中国荷斯坦牛Nramp1基因exon 11的基因多态性, 并分析了其不同基因型与乳房炎及产奶量性状的关系。结果表明: 实验群体发现3种基因型AA、AB、BB, 其中A等位基因为优势等位基因, 等位基因频率为0.767, 而B等位基因频率则为0.233。经χ2适合性检验, 群体处于Hardy-Weinberg平衡状态(P>0.05)。测序结果显示: 扩增片段分别在200 bp(C/G)和254 bp(T/G)存在碱基突变, 并导致了氨基酸改变, 分别为丙氨酸替换为脯氨酸(Ala356Pro)、亮氨基酸替换为蛋氨酸(Leu374Met)。通过构建最小二乘线性模型, 进行Nramp1基因多态性与产奶量、体细胞评分(SCS)的相关性分析表明, AA型个体的SCS最小二乘均值显著低于BB﹑AB型(P<0.05), 而AA型﹑AB个体的产奶量最小二乘均值显著高于BB型(P<0.01, P<0.05), AA基因型可作为乳房炎抗性的优良基因型。因此, 可将Nramp1作为奶牛乳房炎候选基因应用于分子标记辅助选择育种。     

An analysis of the genetic relationship between udder health and udder conformation traits in South African Jersey cows

A multi-trait animal model was used to estimate genetic parameters among lactation somatic cell score (SCS) and udder-type traits in South African Jersey cattle, through restricted maximum likelihood (REML) procedures. Data comprised records on 18 321 Jersey cows in 470 herds, collected through the National Milk Recording Scheme from 1996 to 2002. Average SCS in the first three lactations (SCS1, SCS2 and SCS3) were considered as different traits and the udder-type traits were fore udder attachment (FUA), rear udder height (RUH), rear udder width (RUW), udder cleft (UC), udder depth (UD), fore teat placement (FTP), rear teat placement (RTP) and fore teat length (FTL). Heritability estimates for the respective lactation SCS were 0.07 ± 0.01, 0.11 ± 0.01 and 0.11 ± 0.02. Udder-type traits had heritability estimates ranging from 0.14 ± 0.01 for UD to 0.30 ± 0.02 for FTL. Genetic correlations between SCS and udder-type traits ranged from -0.003 ± 0.07 between FUA and SCS3 to -0.50 ± 0.07 between UD and SCS3. Slow genetic progress is expected when selection is applied independently on SCS and udder-type traits, due to the generally low heritabilities. Tightly attached shallow udders with narrowly placed rear teats are associated with low SCS in the Jersey population.     

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

In 2010, a routine genetic evaluation on occurrence of clinical mastitis in three main dairy cattle breeds – Montbéliarde (MO), Normande (NO) and Holstein (HO) – was implemented in France. Records were clinical mastitis events reported by farmers to milk recording technicians and the analyzed trait was the binary variable describing the occurrence of a mastitis case within the first 150 days of the first three lactations. Genetic parameters of clinical mastitis were estimated for the three breeds. Low heritability estimates were found: between 2% and 4% depending on the breed. Despite its low heritability, the trait exhibits genetic variation so efficient genetic improvement is possible. Genetic correlations with other traits were estimated, showing large correlations (often>0.50, in absolute value) between clinical mastitis and somatic cell score (SCS), longevity and some udder traits. Correlation with milk yield was moderate and unfavorable (ρ=0.26 to 0.30). High milking speed was genetically associated with less mastitis in MO (ρ=−0.14) but with more mastitis in HO (ρ=0.18). A two-step approach was implemented for routine evaluation: first, a univariate evaluation based on a linear animal model with permanent environment effect led to pre-adjusted records (defined as records corrected for all non-genetic effects) and associated weights. These data were then combined with similar pre-adjusted records for others traits in a multiple trait BLUP animal model. The combined breeding values for clinical mastitis obtained are the official (published) ones. Mastitis estimated breeding values (EBV) were then combined with SCSs EBV into an udder health index, which receives a weight of 14.5% to 18.5% in the French total merit index (ISU) of the three breeds. Interbull genetic correlations for mastitis occurrence were very high (ρ=0.94) with Nordic countries, where much stricter recording systems exist reflecting a satisfactory quality of phenotypes as reported by the farmers. They were lower (around 0.80) with countries supplying SCS as a proxy for the international evaluation on clinical mastitis.     

Three novel SNPs of the bovine Tf gene in Chinese native cattle and their associations with milk production traits

Transferrin (Tf) is a β-globulin protein that transports iron ions in mammalian cells. It contributes to innate immunity to microbial pathogens, primarily by limiting microbial access to iron. Thus, polymorphisms present in bovine Tf could potentially underlie inherited differences in mastitis resistance and milk production traits. We detected three novel single-nucleotide polymorphisms of the Tf gene in Chinese native cattle by screening for genetic variation of Tf in 751 individuals of three Chinese cattle breeds, namely China Holstein, Luxi Yellow and Bohai Black, using PCR-RFLP and DNA sequencing techniques. The three new SNPs, g.-1748G>A ss250608649, g.13942T>C ss250608650, and g.14037A>G ss250608651, had allele frequencies of 85.9, 86.3 and 92.5%, 64.5, 73.3 and 65.0%, and 67.6, 73.7 and 60.0%, respectively. SNP g.-1748G>A was located in the 5' flanking region of Tf. SNP g.14037A>G was located in intron 8 of Tf. SNP g.13942T>C, located in exon 8 of Tf, was a synonymous mutation (TTA > CTA), encoding a leucine (326 aa) in the Tf protein. Associations of the Tf SNPs with milk traits were also analyzed. Significant (P < 0.05) relationships among the Tf polymorphisms, somatic cell scores (SCS), and milk productive traits were observed. Cows with genotypes TT (g.13942T>C), GG (g.-1748G>A) and AG (g.14037A>G) had a lower SCS and higher protein levels and 305-day milk yield. Nineteen combinations of different haplotypes from the three SNPs were identified in Chinese Holstein cattle. The haplotype combination ATA/GCA, GCA/GCA and GCG/ GTA was dominant in cows with a lower SCS, a higher protein level and a higher 305-day milk yield, respectively. Moreover, the gene expression level of Tf was higher in mastitis-affected mammary tissues than in normal mammary tissues. These results suggest that the Tf gene affects milk production, as well as mastitis-resistance traits, in Chinese Holsteins.     

Dairy cattle in a temperate climate: the effects of weather on milk yield and composition depend on management

A better understanding of how livestock respond to weather is essential to enable farming to adapt to a changing climate. Climate change is mainly expected to impact dairy cattle through heat stress and an increase in the frequency of extreme weather events. We investigated the effects of weather on milk yield and composition (fat and protein content) in an experimental dairy herd in Scotland over 21 years. Holstein Friesian cows were either housed indoors in winter and grazed over the summer or were continuously housed. Milk yield was measured daily, resulting in 762 786 test day records from 1369 individuals, and fat and protein percentage were sampled once a week, giving 89 331 records from 1220 cows/trait. The relative influence of 11 weather elements, measured from local outdoor weather stations, and two indices of temperature and humidity (THI), indicators of heat stress, were compared using separate maximum likelihood models for each element or index. Models containing a direct measure of temperature (dry bulb, wet bulb, grass or soil temperature) or a THI provided the best fits to milk yield and fat data; wind speed and the number of hours of sunshine were most important in explaining protein content. Weather elements summarised across a week’s timescale from the test day usually explained milk yield and fat content better than shorter-scale (3 day, test day, test day −1) metrics. Then, examining a subset of key weather variables using restricted maximum likelihood, we found that THI, wind speed and the number of hours of sunshine influenced milk yield and composition. The shape and magnitude of these effects depended on whether animals were inside or outside on the test day. The milk yield of cows outdoors was lower at the extremes of THI than at average values, and the highest yields were obtained when THI, recorded at 0900 h, was 55 units. Cows indoors decreased milk yield as THI increased. Fat content was lower at higher THIs than at intermediate THIs in both environments. Protein content decreased as THI increased in animals kept indoors and outdoors, and the rate of decrease was greater when animals were outside than when they were inside. Moderate wind speeds appeared to alleviate heat stress. These results show that milk yield and composition are impacted at the upper extreme of THI under conditions currently experienced in Scotland, where animals have so far experienced little pressure to adapt to heat stress.     

Genetic correlations between production and disease traits during first lactation in Holstein cows

The aim of this study was to estimate genetic correlations between milk yield, somatic cell score (SCS), mastitis, and claw and leg disorders (CLDs) during first lactation in Holstein cows by using a threshold–linear random regression test-day model. We used daily records of milk, fat and protein yields; somatic cell count (SCC); and mastitis and CLD incidences from 46 771 first-lactation Holstein cows in Hokkaido, Japan, that calved between 2000 and 2009. A threshold animal model for binary records (mastitis and CLDs) and linear animal model for yield traits were applied in our multiple trait analysis. For both liabilities and yield traits, additive genetic effects were used as random regression on cubic Legendre polynomials of days on milk. The highest positive genetic correlations between yields and disease incidences (0.36 for milk and mastitis, 0.56 for fat and mastitis, 0.24 for protein and mastitis, 0.32 for milk and CLD, 0.44 for fat and CLD and 0.31 for protein and CLD) were estimated at about the time of peak milk yield (36 to 65 days in milk). Selection focused on early lactation yield may therefore increase the risk of mastitis and CLDs. The positive genetic correlations of SCS with mastitis or CLD incidence imply that selection to reduce SCS in the early stages of lactation would decrease the incidence of both mastitis and CLD.     

Genetic relationship among somatic cell score and some milking traits in Holstein-Friesian primiparous cows milked by an automated milking system

The automated milking system provides breeders with a large amount of automatically collected information about each cow in herd that cannot be easily obtained in non-robotised systems. This knowledge can be used in breeding programs improving somatic cell count (SCC) level. The objective of this study was to estimate heritabilities and genetic correlations among test-day (TD) somatic cell score (SCS) and selected milking traits, such as daily milk yield (MY), milking frequency (MF), milking time (MT) and milking speed (MS), attachment time (AT) to single teat cups, electrical conductivity (EC) and milk temperature (MTEMP). Data were collected for 1 899 Polish Holstein-Friesian primiparous cows milked in an automatic milking system. Genetic parameters of the studied traits were estimated using Bayesian method via Gibbs sampling and two-trait random regression animal model with fixed effect of herd x TD, fixed regressions on days in milk (DIM) nested within age at calving by season of calving and RR for additive genetic and permanent environmental effects. Both fixed and RR were fitted with fourth-order Legendre polynomials on DIM. The estimated daily heritabilities were in the following ranges: MY – 0.162–0.338, MF – 0.156–0.444, MT – 0.090–0.320, MS – 0.252–0.665, AT – 0.105–0.394, EC – 0.269–0.466, MTEMP – 0.135–0.304 and SCS – 0.155–0.321. The heritabilities for traits expressed on a 305-d basis were moderate to high: 0.460 for MY, 0.514 for MF, 0.315 for MT, 0.431 for MS, 0.256 for AT, 0.386 for EC, 0.407 for MTEMP and 0.359 for SCS. Genetic correlations between traits on a 305-d basis showed that SCS was most strongly genetically correlated with MTEMP (0.572) and MS (0.480), whereas genetic relationships of SCS with MT (0.221) and EC (− 0.216) were moderate. Phenotypic correlations between traits on a 305-d basis were moderate or low. Somatic cell score was negatively phenotypically correlated with MY, MF and MT, with the highest relationship with MT (− 0.302). The largest positive phenotypic correlations were observed between SCS and MS (0.311) and with MTEMP (0.286). In summary, it is concluded that there is a chance to carry out effective selection for lower SCS and for some other traits, in particular MS and MTEMP. The obtained results are promising enough to conduct further research to evaluate how these traits can be used both to increase the accuracy of genetic evaluations of SCC and to improve udder health.     

Modeling heat stress effects on dairy cattle milk production in a tropical environment using test-day records and random regression models

In tropical environments, dairy cattle production is constrained by several factors, including climate. The seasonal loss of milk due to heat stress is a recurring challenge for many dairy producers. The objective of this study was to detect heat stress thresholds, milk yield loss and individual animal variations using random regression models for dairy cattle from test-day milk records. Data were obtained from the Kenya Livestock Breeders Organization for the years 2000–2017 and merged with weather data. The weather parameters were grid-interpolated solar and meteorological data obtained from the National Aeronautics and Space Administration/Prediction Of Worldwide Energy Resources (NASA/POWER). After editing, the records comprised 49 993, 45 251 and 36 136 test-day records for first, second, and third lactations, respectively, for the four main dairy breeds: Friesian (68.0%), Ayrshire (21.1%), Jersey (7.6%) and Guernsey (3.3%). Variance components were estimated using Restricted Maximum Likelihood in ASReml software. Random regression models with third-order Legendre polynomials were fitted to the average and individual lactation curves and the reaction norms. An extended factor analytic variance structure for the random cow effects was used to estimate (co)variances between days in milk and thermal load. The daily average temperature (TA) and temperature humidity index (THI) were identified as the most suitable thermal load indicators for assessing milk yield losses. Considering a one day lag, the estimated heat stress thresholds were about 22 °C and 69 index units for TA and THI, respectively. Almost no differences were observed for estimated residual variances between the thermal load indicators, indicating there was no better model fit by TA or THI. The heat stress thresholds and milk loss patterns are important for management of dairy production systems in the tropics with climatic conditions similar to this study. Data recording should be improved as a tool to monitor the expected impacts of climate change and mitigation measures.     

Bayesian inference of genetic parameters for test-day milk yield, milk quality traits, and somatic cell score in Burlina cows

The aim of the study was to infer (co)variance components for daily milk yield, fat and protein contents, and somatic cell score (SCS) in Burlina cattle (a local breed in northeast Italy). Data consisted of 13576 monthly test-day records of 666 cows (parities 1 to 8) collected in 10 herds between 1999 and 2009. Repeatability animal models were implemented using Bayesian methods. Flat priors were assumed for systematic effects of herd test date, days in milk, and parity, as well as for permanent environmental, genetic, and residual effects. On average, Burlina cows produced 17.0 kg of milk per day, with 3.66 and 3.33% of fat and protein, respectively, and 358000 cells per mL of milk. Marginal posterior medians (highest posterior density of 95%) of heritability were 0.18 (0.09–0.28), 0.28 (0.21–0.36), 0.35 (0.25–0.49), and 0.05 (0.01–0.11) for milk yield, fat content, protein content, and SCS, respectively. Marginal posterior medians of genetic correlations between the traits were low and a 95% Bayesian confidence region included zero, with the exception of the genetic correlation between fat and protein contents. Despite the low number of animals in the population, results suggest that genetic variance for production and quality traits exists in Burlina cattle.