Genetic analysis of health traits and their associations with longevity, fertility, production, and conformation traits in Holstein cattle

Health traits have high economic values in dairy cattle breeding, which can cause considerable financial loss through involuntary culling. In this study, fourteen health traits were analysed, including five composite health traits: reproductive disorders, udder health (UH), digestive disorders, metabolic disorders, locomotory diseases (LD), and nine independent health traits: gestation disorders and peripartum disorders, irregular estrus cycle and sterility, metritis (ME), mastitis (MA), abomasal displacement (AD), enteritis (EN), and ketosis, claw diseases (CD), laminitis complex. This study analysed variance components for health traits through both single and bivariate repeatability animal models. All health traits showed low heritability, ranging from 0.001 to 0.025. Most of the health traits in five categories showed negative genetic correlations, ranging from -0.012 (CD and EN) to -0.634 (ME and EN). Strong positive genetic correlations appeared within the same category, ranging from 0.469 (EN and AD) to 0.994 (UH and MA, LD and CD). Furthermore, approximate genetic correlations were evaluated between health traits and routinely collected traits (longevity, fertility, production, and conformation). In general, the low to moderate approximate genetic correlations were estimated between health traits and routinely collected traits. The estimated correlations between health traits and longevity, fertility, production, and conformation traits could provide an indirect reference for disease-resistance breeding in Holstein cattle.

Genetic analysis of longevity and their associations with fertility traits in Holstein cattle

The increase of longevity is intended to reduce involuntary culling rates, not extend the life span, and it reflects the ability of animals to successfully cope with the environment and disease during production. Sire model, animal model and repeatability animal models were used to estimate the (co) variance components of longevity and fertility traits. Six longevity and thirteen fertility traits were analysed, including herd life (HL), productive life (PL), number of days between first calving and the end of first lactation or culling (L1); number of days between first calving and the end of the second lactation or culling (L2); number of days between first calving and the end of the third lactation or culling (L3); number of days between first calving and the end of the fourth lactation or culling (L4); age at first service, age at first calving (AFC), the interval from first to last inseminations in heifer (IFLh), conception rate of first insemination in heifer, days open (DO), calving interval, gestation length, interval from calving to first insemination (ICF), interval from first to last inseminations in cow (IFLc), conception rate of first insemination in cow, calving ease (CE), birth weight, and calf survival. The estimated heritabilities (±SE) were 0.018 (±0.003), 0.015 (±0.003), 0.049 (±0.004), 0.025 (±0.003), 0.009 (±0.002) and 0.011 (±0.002) for HL, PL, L1, L2, L3 and L4, respectively. Strong correlations were appeared in HL and PL; the genetic and phenotypic correlation coefficients were 0.998 and 0.985, respectively. There were high genetic and phenotypic correlations which were observed in L1 and L2, L2 and L3, L3 and L4, respectively. All fertility traits of heifer showed medium to high heritability, while the cow showed low heritability. All heifer fertility traits had low genetic associations with longevity traits, ranging from -0.018 (L2 and IFLh) to 0.257 (L3 and AFC). Most of the fertility traits showed negative correlations with longevity traits in different parities, and we recommend DO, ICF, IFLc and CE as indirect indicators of longevity traits in dairy cows, but we also need to take into account the differences between parities.

Reduction in Mortality of Calves with Bovine Respiratory Disease in Detection with Influenza C and D Virus

Both influenza C (ICV) and influenza D (IDV) viruses were recently included as bovine respiratory disease (BRD) causes, but their role in BRD has not been evaluated. Therefore, the mortality and reproductive performances of BRD calves with different isolated viruses were determined in this study. Data on 152 BRD calves with bovine viral diarrhoea virus (BVDV), bovine respiratory syncytial virus (BRSV), bovine coronavirus (BCoV), bovine parainfluenza virus 3 (BPIV-3), ICV, or IDV from nasal swab samples using real-time rt-PCR were used. The general data and respiratory signs were recorded immediately, and thereafter, the data on dead or culling calves due to BRD and reproductive performance were collected. The percentages of the BRD calves were 71.7%, 52.6%, 40.8%, 10.5%, 68.4%, and 65.8% for BVDV, BRSV, BCoV, BPIV-3, ICV, and IDV, respectively. Mucous secretion (OR = 4.27) and age ≤ 6 months (OR =14.97) had higher risks of mortality than those with serous secretion and older age. The calves with IDV had lower risks of culling than those without IDV (OR = 0.19). This study shows that most viral infections in BRD calves are a combination of viruses with BVDV, ICV, and IDV. In addition, IDV might have a role in reducing the severity of BRD calves.

An Appropriate Genetic Approach for Improving Reproductive Traits in Crossbred Thai-Holstein Cattle under Heat Stress Conditions

Thailand is a tropical country affected by global climate change and has high temperatures and humidity that cause heat stress in livestock. A temperature−humidity index (THI) is required to assess and evaluate heat stress levels in livestock. One of the livestock types in Thailand experiencing heat stress due to extreme climate change is crossbred dairy cattle. Genetic evaluations of heat tolerance in dairy cattle have been carried out for reproductive traits. Heritability values for reproductive traits are generally low (<0.10) because environmental factors heavily influence them. Consequently, genetic improvement for these traits would be slow compared to production traits. Positive and negative genetic correlations were found between reproductive traits and reproductive traits and yield traits. Several selection methods for reproductive traits have been introduced, i.e., the traditional method, marker-assisted selection (MAS), and genomic selection (GS). GS is the most promising technique and provides accurate results with a high genetic gain. Single-step genomic BLUP (ssGBLUP) has higher accuracy than the multi-step equivalent for fertility traits or low-heritability traits.

Genomic prediction of milk-production traits and somatic cell score using single-step genomic best linear unbiased predictor with random regression test-day model in Thai dairy cattle

Cow genotypes are expected to improve the accuracy of genomic estimated breeding values (GEBV) for young bulls in relatively small populations such as Thai Holstein-Friesian crossbred dairy cattle in Thailand. The objective of this study was to investigate the effect of cow genotypes on the predictive ability and individual accuracies of GEBV for young dairy bulls in Thailand. Test-day data included milk yield (n = 170,666), milk component traits (fat yield, protein yield, total solids yield, fat percentage, protein percentage, and total solids percentage; n = 160,526), and somatic cell score (n = 82,378) from 23,201, 82,378, and 13,737 (for milk yield, milk component traits, and SCS, respectively) cows calving between 1993 and 2017, respectively. Pedigree information included 51,128; 48,834; and 32,743 animals for milk yield, milk component traits, and somatic cell score, respectively. Additionally, 876, 868, and 632 pedigreed animals (for milk yield, milk component traits, and SCS, respectively) were genotyped (152 bulls and 724 cows), respectively, using Illumina Bovine SNP50 BeadChip. We cut off the data in the last 6 yr, and the validation animals were defined as genotyped bulls with no daughters in the truncated set. We calculated GEBV using a single-step random regression test-day model (SS-RR-TDM), in comparison with estimated breed value (EBV) based on the pedigree-based model used as the official method in Thailand (RR-TDM). Individual accuracies of GEBV were obtained by inverting the coefficient matrix of the mixed model equations, whereas validation accuracies were measured by the Pearson correlation between deregressed EBV from the full data set and (G)EBV predicted with the reduced data set. When only bull genotypes were used, on average, SS-RR-TDM increased individual accuracies by 0.22 and validation accuracies by 0.07, compared with RR-TDM. With cow genotypes, the additional increase was 0.02 for individual accuracies and 0.06 for validation accuracies. The inflation of GEBV tended to be reduced using cow genotypes. Genomic evaluation by SS-RR-TDM is feasible to select young bulls for the longitudinal traits in Thai dairy cattle, and the accuracy of selection is expected to be increased with more genotypes. Genomic selection using the SS-RR-TDM should be implemented in the routine genetic evaluation of the Thai dairy cattle population. The genetic evaluation should consider including genotypes of both sires and cows.

Detection of Genomic Regions with Pleiotropic Effects for Growth and Carcass Quality Traits in the Rubia Gallega Cattle Breed

Simple Summary

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. We have developed a ssGWAS by backsolving the SNP effects after implementing a ssGBLUP. The results showed an apparent heterogeneity of the additive genetic variance across the genome. Some of the genomic regions explaining the most of this additive variance were shared across traits, indicating the presence of pleiotropic effects, which were reflected in their genetic correlations.

Abstract

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. In recent years, genomic evaluation has been implemented by using a ssGBLUP (single-step Genomic Best Linear Unbiased Prediction). This procedure can reparameterized to perform ssGWAS (single-step Genome Wide Association Studies) by backsolving the SNP (single nucleotide polymorphisms) effects. Therefore, the objective of this study was to identify genomic regions associated with the genetic variability in growth and carcass quality traits. We implemented a ssGBLUP by using a database that included records for Birth Weight (BW-327,350 records-), Weaning Weight (WW-83,818-), Cold Carcass Weight (CCW-91,621-), Fatness (FAT-91,475-) and Conformation (CON-91,609-). The pedigree included 464,373 individuals, 2449 of which were genotyped. After a process of filtering, we ended up using 43,211 SNP markers. We used the GBLUP and SNPBLUP model equivalences to obtain the effects of the SNPs and then calculated the percentage of variance explained by the regions of the genome between 1 Mb. We identified 7 regions of the genome for CCW; 8 regions for BW, WW, FAT and 9 regions for CON, which explained the percentage of variance above 0.5%. Furthermore, a number of the genome regions had pleiotropic effects, located at: BTA1 (131–132 Mb), BTA2 (1–11 Mb), BTA3 (32–33 Mb), BTA6 (36–38 Mb), BTA16 (24–26 Mb), and BTA 21 (56–57 Mb). These regions contain, amongst others, the following candidate genes: NCK1, MSTN, KCNA3, LCORL, NCAPG, and RIN3.

Genetic correlations between fertility traits and milk composition and fatty acids in Holstein-Friesian, Brown Swiss, and Simmental cattle using recursive models

This study aimed to investigate the genetic and putative causal relationships between fertility traits [i.e., days open (DO) and calving rate (CR)] and milk quality, composition, and fatty acid contents (milk composition traits) in Holstein-Friesian, Brown Swiss, and Simmental cattle, using recursive models within a Bayesian framework. Trivariate animal models were run, each including one fertility trait, one milk composition trait, and, in all models, milk yield. The DO and CR data were merged with the test days closest to the insemination date for milk composition traits. After editing, 16,468 to 23,424 records for Holstein-Friesian, 23,424 to 46,660 for Brown Swiss, and 26,105 to 35,574 for Simmental were available for the analyses. Recursive animal models were applied to investigate the possible causal influences of milk composition traits on fertility and the genetic relationships among these traits. The results suggested a potential cause-and-effect relationship between milk composition traits and fertility traits, with the first trait influencing the latter. We also found greater recursive effects of milk composition traits on DO than on CR, the latter with some putative differences among breeds in terms of sensitivity. For instance, the putative causal effects of somatic cell score on CR (on the observed scale, %) varied from -0.96 to -1.39%, depending on the breed. Concerning fatty acids, we found relevant putative effects of C18:0 on CR, with estimates varying from -7.8 to -9.9%. Protein and casein percentages, and short-chain fatty acid showed larger recursive effects on CR, whereas fat, protein, and casein percentages, C16:0, C18:0, and long-chain fatty acid had larger effects on DO. The results obtained suggested that these milk traits could be considered as effective indicators of the effects of animal metabolic and physiological status on fertility.

Genetic parameters for fertility traits assessed in herds divergent in milk energy output in Holstein-Friesian, Brown Swiss, and Simmental cattle

In this study, we aimed to investigate differences in the genetics of fertility traits (heritability of traits and correlations between traits in divergent environments) in dairy cows of different production levels defined on the basis of the herd-average daily milk energy output (herd-dMEO). Data were obtained from Holstein-Friesian (n = 37,359 for fertility traits, 381,334 for dMEO), Brown Swiss (n = 79,638 for fertility traits, 665,697 for dMEO), and Simmental cows (n = 63,048 for fertility traits, 448,445 for dMEO) reared in northeastern Italy. Fertility traits under study were interval from calving to first service, interval from first service to conception, days open, calving interval, calving rate, and nonreturn rate at d 56. We classified herds into low and high productivity based on the herd-average dMEO (inferred using mixed effects models). We estimated genetic parameters using Bayesian bivariate animal models, where expressions of a phenotype in the low and high dMEO herds were taken as being different-albeit correlated-traits. Fertility traits were more favorable in Simmental than in Holstein-Friesian cows, whereas for all traits, Holstein-Friesian had the highest estimates of intraherd heritability [ranging from 0.021 (0.006-0.038) to 0.126 (0.10-0.15)] and Simmental the lowest [ranging from 0.008 (0.001-0.017) to 0.101 (0.08-0.12)]. The genetic correlations between fertility traits and dMEO were moderate and unfavorable, ranging, in absolute values, from 0.527 (0.37-0.68) to 0.619 (0.50-0.73) in Holstein-Friesian; from 0.339 (0.20-0.47) to 0.556 (0.45-0.66) in Brown Swiss; and from 0.340 (0.10-0.60) to 0.475 (0.33-0.61) in Simmental cattle. The only exception was the nonreturn rate at d 56, which had weak genetic correlations with dMEO in all 3 breeds. The herd correlations between fertility and dMEO tended to be modest and favorable and the residual correlations modest and variable. The heritability of fertility traits tended to be greater in the low dMEO than in the high dMEO herds in the case of the Holstein-Friesians, but not in the case of the Brown Swiss or Simmentals. The additive genetic correlations between fertility traits in the low and high dMEO herds were always lower than 1 [0.329 (-0.17 to 0.85) to 0.934 (0.86 to 0.99)] for all traits considered in all breeds. The correlation was particularly low for the threshold characters and the interval from first service to conception in Holstein-Friesian, suggesting that the relative performances of genotypes vary significantly between herds of different dMEO levels. Although there was large variability in the estimates, results might support making separate genetic evaluations of fertility in the different herd production groups. Our results also indicate that Simmental, a dual-purpose breed, has higher fertility and lower environmental sensitivity than Holstein-Friesian, with Brown Swiss being intermediate.

Genetic associations between milk fat-to-protein ratio, milk production and fertility in the first two lactations of Thai Holsteins dairy cattle

The aims of this study were to estimate, simultaneously, the genetic parameters of test-day milk fat-to-protein ratio (FPR), test-day milk yield (MY), and days-open (DO) in the first two lactations of Thai Holsteins. A total of 76 194 test-day production records collected from 8874 cows with 8674 DO records between 2001 and 2011 from different lactations were treated as separated traits. The estimates of heritability for test-day FPR in the first lactation showed an increasing trend, whereas the estimates in the second lactation showed a U-shape trend. Genetic correlations for FPR-DO and MY-DO showed a decreasing trend along days in milk (DIM) in both lactations, whereas genetic correlations for FPR-MY increased along DIM in the first lactation but decreased in the second lactation. Genetic correlations of FPR between consecutive DIM were moderate to high, which showed the effectiveness of simultaneous analyses. Selection of FPR in the early stage has no adverse effect on MY and DO for the first lactation but has a negative effect on MY and positive effect on DO for the second lactation. This study showed that genetic improvement of the energy balance using FPR, MY and DO with multi-trait test day model could be applied in a Thailand dairy cattle breeding program.