Fourier transform infrared milk spectral data are routinely used for milk quality control and have been revealed to be driven by genetics. This study aimed to (1) estimate heritability for 1,060 wavenumbers in the infrared region from 5,008 to 925 cm −1 , (2) estimate genomic correlations between wavenumbers with increased heritability, and (3) compare results between Danish Holstein and Danish Jersey cows. For Danish Holstein, 3,275 cows and 19,656 milk records were available. For Danish Jersey, 3,408 cows and 20,228 milk records were available. We used a hierarchical mixed model, with a Bayesian approach. Heritability of individual wavenumbers ranged from 0.00 to 0.31 in Danish Holstein, and from 0.00 to 0.30 in Danish Jersey. Genomic correlation was calculated between 15 selected wavenumbers, and varied from weak to very strong, in both Danish Holstein and Danish Jersey (0.03 to 0.97, and −0.11 to −0.97). Within the 15 selected wavenumbers, a subdivision into 2 groups of wavenumbers was observed, where genomic correlations were negative between groups, and positive within groups. Heritability and genomic correlations were higher in Danish Holstein compared with Danish Jersey, but followed a similar pattern in both breeds. Breed differences were most pronounced in the mid-infrared region that interacts with lactose and the spectral region that interacts with protein.In conclusion, heritability for individual wavenumbers of Fourier transform milk spectra was moderate, and strong genomic correlations were observed between wavenumbers across the spectrum. Heritability and genomic correlations were higher in Danish Holstein, with the strongest breed differences showing in spectral regions interacting with protein or lactose. bers in breeding programs, or use genetic correlations between wavenumbers and a trait of interest.The aims of this study were to (1) calculate heritability (h 2 ) for 1,060 wavenumbers in the MIR region from 5,008 to 925 cm −1 , (2) calculate genomic correlations between transmittance values of wavenumbers with high h 2 , and (3) compare results between Danish Holstein (DH) and Danish Jersey (DJ) cows.
MATERIALS AND METHODS
Study PopulationPhenotypes. Milk records were provided by the Danish milk recording organization (RYK, Aarhus, Denmark). The study population consisted of 3,275 DH cows from 354 farms, and 3,408 DJ cows from 175 farms. On 4 farms, both DH and DJ cows were present. From farms with only DH, 1 to 112 cows were sampled, and from farms with only DJ, 1 to 217 cows were sampled. Milk records were sent to Eurofins-Steins laboratory (Vejen, Denmark) for FT-IR spectral analyses using MilkoScan FT+ (Foss, Hillerød, Denmark) and for analysis of the infrared region 5,008 to 925 cm −1 , including 1,060 individual wavenumbers.Visual inspection of spectra was performed using PLSR Hotelling T 2 versus Q-residual plots and leverage versus studentized residual plots, and outliers were removed. Milk records with outlying fat percentage (fat%), and protein percentage (protein%) were excluded fr...