Genome-Wide Association Study to Identify QTL for Carcass Traits in Korean Hanwoo Cattle

This study aimed to identify genetic associations with carcass traits in Hanwoo cattle using a genome-wide association study. A total of 9302 phenotypes were analyzed, and all animals were genotyped using the Illumina Bovine 50K v.3 SNP chip. Heritabilities for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) were estimated as 0.42, 0.36, 0.36, and 0.47, respectively, using the GBLUP model, and 0.47, 0.37, 0.36, and 0.42, respectively, using the Bayes B model. We identified 129 common SNPs using DGEBV and 118 common SNPs using GEBV on BTA6, BTA13, and BTA14, suggesting their potential association with the traits of interest. No common SNPs were found between the GBLUP and Bayes B methods when using residuals as a response variable in GWAS. The most promising candidate genes for CWT included SLIT2, PACRGL, KCNIP4, RP1, XKR4, LYN, RPS20, MOS, FAM110B, UBXN2B, CYP7A1, SDCBP, NSMAF, TOX, CA8, LAP3, FAM184B, and NCAPG. For EMA, the genes IBSP, LAP3, FAM184B, LCORL, NCAPG, SLC30A9, and BEND4 demonstrated significance. Similarly, CYP7B1, ARMC1, PDE7A, and CRH were associated with BF, while CTSZ, GNAS, VAPB, and RAB22A were associated with MS. This finding offers valuable insights into genomic regions and molecular mechanisms influencing Hanwoo carcass traits, aiding efficient breeding strategies.

Strategies of preserving genetic diversity while maximizing genetic response from implementing genomic selection in pulse breeding programs

KEY MESSAGE

Genomic selection maximizes genetic gain by recycling parents to germplasm pool earlier and preserves genetic diversity by restricting the number of fixed alleles and the relationship in pulse breeding programs. Using a stochastic computer simulation, we investigated the benefit of optimization strategies in the context of genomic selection (GS) for pulse breeding programs. We simulated GS for moderately complex to highly complex traits such as disease resistance, grain weight and grain yield in multiple environments with a high level of genotype-by-environment interaction for grain yield. GS led to higher genetic gain per unit of time and higher genetic diversity loss than phenotypic selection by shortening the breeding cycle time. The genetic gain obtained from selecting the segregating parents early in the breeding cycle (at F1 or F2 stages) was substantially higher than selecting at later stages even though prediction accuracy was moderate. Increasing the number of F1 intercross (F1i) families and keeping the total number of progeny of F1i families constant, we observed a decrease in genetic gain and increase in genetic diversity, whereas increasing the number of progeny per F1i family while keeping a constant number of F1i families increased the rate of genetic gain and had higher genetic diversity loss per unit of time. Adding 50 F2 family phenotypes to the training population increased the accuracy of genomic breeding values (GEBVs) and genetic gain per year and decreased the rate of genetic diversity loss. Genetic diversity could be preserved by applying a strategy that restricted both the percentage of alleles fixed and the average relationship of the group of selected parents to preserve long-term genetic improvement in the pulse breeding program.

Genotype by heat stress interactions for production and functional traits in dairy cows from an across-generation perspective

The aim of this study was to analyze time-lagged heat stress (HS) effects during late gestation on genetic co(variance) components in dairy cattle across generations for production, female fertility, and health traits. The data set for production and female fertility traits considered 162,492 Holstein Friesian cows from calving years 2003 to 2012, kept in medium-sized family farms. The health data set included 69,986 cows from calving years 2008 to 2016, kept in participating large-scale co-operator herds. Production traits were milk yield (MKG), fat percentage (fat%), and somatic cell score (SCS) from the first official test-day in first lactation. Female fertility traits were the nonreturn rate after 56 d (NRR56) in heifers and the interval from calving to first insemination (ICFI) in first-parity cows. Health traits included clinical mastitis (MAST), digital dermatitis (DD), and endometritis (EM) in the early lactation period in first-parity cows. Meteorological data included temperature and humidity from public weather stations in closest herd distance. The HS indicator was the temperature-humidity index (THI) during dams' late gestation, also defined as in utero HS. For the genetic analyses of production, female fertility, and health traits in the offspring generation, a sire-maternal grandsire random regression model with Legendre polynomials of order 3 for the production and of order 2 for the fertility and health traits on prenatal THI, was applied. All statistical models additionally considered a random maternal effect. THI from late gestation (i.e., prenatal climate conditions), influenced genetic parameter estimates in the offspring generation. For MKG, heritabilities and additive genetic variances decreased in a wave-like pattern with increasing THI. Especially for THI >58, the decrease was very obvious with a minimal heritability of 0.08. For fat% and SCS, heritabilities increased slightly subjected to prenatal HS conditions at THI >67. The ICFI heritabilities differed marginally across THI [heritability (h²) = 0.02-0.04]. For NRR56, MAST, and DD, curves for heritabilities and genetic variances were U-shaped, with largest estimates at the extreme ends of the THI scale. For EM, heritability increased from THI 25 (h² = 0.13) to THI 71 (h² = 0.39). The trait-specific alterations of genetic parameters along the THI gradient indicate pronounced genetic differentiation due to intrauterine HS for NRR56, MAST, DD, and EM, but decreasing genetic variation for MKG and ICFI. Genetic correlations smaller than 0.80 for NRR56, MAST, DD, and EM between THI 65 with corresponding traits at remaining THI indicated genotype by environment interactions. The lowest genetic correlations were identified when considering the most distant THI. For MKG, fat%, SCS, and ICFI, genetic correlations throughout were larger than 0.80, disproving concerns for any genotype by environment interactions. Variations in genetic (co)variance components across prenatal THI may be due to epigenetic modifications in the offspring genome, triggered by in utero HS. Epigenetic modifications have a persistent effect on phenotypic responses, even for traits recorded late in life. However, it is imperative to infer the underlying epigenetic mechanisms in ongoing molecular experiments.

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.

Suitability of traits related to aggression and handleability for integration into pig breeding programmes: Genetic parameters and comparison between Gaussian and binary trait specifications

Changes in husbandry systems as well as consumers' increasing demands for animal welfare lead to increasing importance of traits such as handleability and aggressiveness in pigs. However, before using such novel traits for selection decisions, information on genetic parameters for these traits for the specific population is required. Therefore, weight gain and behaviour-related traits were recorded in 1004 pigs (814 Pietrain x German Landrace crossbred, 190 German Landrace purebred) at different ages. Behaviour indicators and tests were assessed and conducted, respectively under commercial farm conditions and included scoring of skin lesions (twice) and behaviour during backtests (twice), injections (once), handling (twice) and weighing (three times). Since behaviour scores often exhibit suboptimal statistical properties for parametric analyses, variance components were estimated using an animal model assuming a normal (Gaussian, GA; all traits) and additionally a binary distribution of variables (BI; using a logit-link function for all behaviour traits). Heritabilities for behavioural traits ranged from 0.02 ± 0.04 (finishing pig handling test; BI) to 0.36 ± 0.08 (backtest 2; GA) suggesting that some of the traits are potentially useful for genetic selection (e.g. finishing pig weighing test: h2 (GA) = 0.20 ± 0.07). Only minor differences were observed for results from binary and Gaussian analyses of the same traits suggesting that either approach might yield valid results. However, four-fold cross-validation using correlations between breeding values of a sub-set of animals for the sample trait finishing pig weighing score indicated slight superiority of the logit model (r = 0.85 ± 0.04 vs. r = 0.77 ± 0.03). Generally, only weak to moderate associations were found between behavioural reactions to the same test at different ages (rp ≤ 0.11 for weighing at different ages; rp = 0.30 but rg (GA) = 0.84 ± 0.11 for the backtests) as well as between reactions to different tests. Therefore, for inclusion of behaviour traits into breeding programmes, and considering high labour input required for some tests such as the backtest, it is recommended to assess behaviour during situations that are relevant and identical to practical conditions, while the use of indicator traits generally does not appear to be a very promising alternative.

Random forest estimation of genomic breeding values for disease susceptibility over different disease incidences and genomic architectures in simulated cow calibration groups

A simulation study was conducted to investigate the performance of random forest (RF) and genomic BLUP (GBLUP) for genomic predictions of binary disease traits based on cow calibration groups. Training and testing sets were modified in different scenarios according to disease incidence, the quantitative-genetic background of the trait (h(2)=0.30 and h(2)=0.10), and the genomic architecture [725 quantitative trait loci (QTL) and 290 QTL, populations with high and low levels of linkage disequilibrium (LD)]. For all scenarios, 10,005 SNP (depicting a low-density 10K SNP chip) and 50,025 SNP (depicting a 50K SNP chip) were evenly spaced along 29 chromosomes. Training and testing sets included 20,000 cows (4,000 sick, 16,000 healthy, disease incidence 20%) from the last 2 generations. Initially, 4,000 sick cows were assigned to the testing set, and the remaining 16,000 healthy cows represented the training set. In the ongoing allocation schemes, the number of sick cows in the training set increased stepwise by moving 10% of the sick animals from the testing set to the training set, and vice versa. The size of the training and testing sets was kept constant. Evaluation criteria for both GBLUP and RF were the correlations between genomic breeding values and true breeding values (prediction accuracy), and the area under the receiving operating characteristic curve (AUROC). Prediction accuracy and AUROC increased for both methods and all scenarios as increasing percentages of sick cows were allocated to the training set. Highest prediction accuracies were observed for disease incidences in training sets that reflected the population disease incidence of 0.20. For this allocation scheme, the largest prediction accuracies of 0.53 for RF and of 0.51 for GBLUP, and the largest AUROC of 0.66 for RF and of 0.64 for GBLUP, were achieved using 50,025 SNP, a heritability of 0.30, and 725 QTL. Heritability decreases from 0.30 to 0.10 and QTL reduction from 725 to 290 were associated with decreasing prediction accuracy and decreasing AUROC for all scenarios. This decrease was more pronounced for RF. Also, the increase of LD had stronger effect on RF results than on GBLUP results. The highest prediction accuracy from the low LD scenario was 0.30 from RF and 0.36 from GBLUP, and increased to 0.39 for both methods in the high LD population. Random forest successfully identified important SNP in close map distance to QTL explaining a high proportion of the phenotypic trait variations.

Assessing the impact of natural service bulls and genotype by environment interactions on genetic gain and inbreeding in organic dairy cattle genomic breeding programs

The objective of the present study was to compare genetic gain and inbreeding coefficients of dairy cattle in organic breeding program designs by applying stochastic simulations. Evaluated breeding strategies were: (i) selecting bulls from conventional breeding programs, and taking into account genotype by environment (G×E) interactions, (ii) selecting genotyped bulls within the organic environment for artificial insemination (AI) programs and (iii) selecting genotyped natural service bulls within organic herds. The simulated conventional population comprised 148 800 cows from 2976 herds with an average herd size of 50 cows per herd, and 1200 cows were assigned to 60 organic herds. In a young bull program, selection criteria of young bulls in both production systems (conventional and organic) were either 'conventional' estimated breeding values (EBV) or genomic estimated breeding values (GEBV) for two traits with low (h 2=0.05) and moderate heritability (h 2=0.30). GEBV were calculated for different accuracies (r mg), and G×E interactions were considered by modifying originally simulated true breeding values in the range from r g=0.5 to 1.0. For both traits (h 2=0.05 and 0.30) and r mg⩾0.8, genomic selection of bulls directly in the organic population and using selected bulls via AI revealed higher genetic gain than selecting young bulls in the larger conventional population based on EBV; also without the existence of G×E interactions. Only for pronounced G×E interactions (r g=0.5), and for highly accurate GEBV for natural service bulls (r mg>0.9), results suggests the use of genotyped organic natural service bulls instead of implementing an AI program. Inbreeding coefficients of selected bulls and their offspring were generally lower when basing selection decisions for young bulls on GEBV compared with selection strategies based on pedigree indices.