Correction: SNV discovery and functional candidate gene identification for milk composition based on whole genome resequencing of Holstein bulls with extremely high and low breeding values

Single-step genomic best linear unbiased predictor genetic parameter estimations and genome-wide associations for milk fatty acid profiles, interval from calving to first insemination, and ketosis in Holstein cattle

Milk fatty acids (FA) have been suggested as biomarkers for early-lactation metabolic diseases and for female fertility status. The aim of the present study was to infer associations between FA, the metabolic disorder ketosis (KET), and the interval from calving to first insemination (ICF) genetically and genomically. In this regard, we focused on a single-step genomic BLUP approach, allowing consideration of genotyped and ungenotyped cows simultaneously. The phenotypic data set considered 38,375 first-lactation Holstein cows, kept in 45 large-scale co-operator herds from 2 federal states in Germany. The calving years for these cows were from 2014 to 2017. Concentrations in milk from the first official milk recording test-day for saturated, unsaturated (UFA), monounsaturated (MUFA), polyunsaturated, palmitic, and stearic (C18:0) FA were determined via Fourier-transform infrared spectroscopy. Ketosis was defined as a binary trait according to a veterinarian diagnosis key, considering diagnoses within a 6-wk interval after calving. A subset of 9,786 cows was genotyped for 40,989 SNP markers. Variance components and heritabilities for all Gaussian distributed FA and for ICF, and for binary KET were estimated by applying single-step genomic BLUP single-trait linear and threshold models, respectively. Genetic correlations were estimated in series of bivariate runs. Genomic breeding values for the single-step genomic BLUP estimations were dependent traits in single-step GWAS. Heritabilities for FA were moderate in the range from 0.09 to 0.20 (standard error = 0.02-0.03), but quite small for ICF (0.08, standard error = 0.01) and for KET (0.05 on the underlying liability scale, posterior standard deviation = 0.02). Genetic correlations between KET and UFA, MUFA, and C18:0 were large (0.74 to 0.85, posterior standard deviation = 0.14-0.19), and low positive between KET and ICF (0.17, posterior standard deviation = 0.22). Genetic correlations between UFA, MUFA, and C18:0 with ICF ranged from 0.34 to 0.46 (standard error = 0.12). In single-step GWAS, we identified a large proportion of overlapping genomic regions for the different FA, especially for UFA and MUFA, and for saturated and palmitic FA. One identical significantly associated SNP was identified for C18:0 and KET on BTA 15. However, there was no genomic segment simultaneously significantly affecting all trait categories ICF, FA, and KET. Nevertheless, some of the annotated potential candidate genes DGKA, IGFBP4, and CXCL8 play a role in lipid metabolism and fertility mechanisms, and influence production diseases in early lactation. Genetic and genomic associations indicate that Fourier-transform infrared spectroscopy FA concentrations in milk from the first official test-day are valuable predictors for KET and for ICF.