Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.     

Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.     

Estimation of genetic parameters and single-step genome-wide association studies for milk urea nitrogen in Holstein cattle

The main objectives of this study were to estimate genetic parameters for milk urea nitrogen (MUN) in Holstein cattle and to conduct a single-step (ss)GWAS to identify candidate genes associated with MUN. Phenotypic measurements from 24,435 Holstein cows were collected from March 2013 to July 2019 in 9 dairy farms located in the Beijing area, China. A total of 2,029 cows were genotyped using the Illumina 150K Bovine Bead Chip, containing 121,188 SNP. A single-trait repeatability model was used to evaluate the genetic background of MUN. We found that MUN is a trait with low heritability (0.06 ± 0.004) and repeatability (0.12). Considering similar milk production levels, a lower MUN concentration indicates higher nitrogen digestibility. The genetic correlations between MUN and milk yield, net energy concentration, fat percentage, protein percentage, and lactose percentage were positive and ranged from 0.02 to 0.26. The genetic correlation between MUN and somatic cell score (SCS) was negative (?0.18), indicating that animals with higher MUN levels tend to have lower SCS. Both ssGWAS and pathway enrichment analyses were used to explore the genetic mechanisms underlying MUN. A total of 18 SNP (located on BTA11, BTA12, BTA14, BTA17, and BTA18) were found to be significantly associated with MUN. The genes CFAP77, CAMSAP1, CACNA1B, ADGRB1, FARP1, and INTU are considered to be candidate genes for MUN. These candidate genes are associated with important biological processes such as protein and lipid metabolism and binding to specific proteins. This set of candidate genes, metabolic pathways, and their functions provide a better understanding of the genomic architecture and physiological mechanisms underlying MUN in Holstein cattle.     

On the genomic regions associated with milk lactose in Fleckvieh cattle

Lactose is a sugar uniquely found in mammals' milk and it is the major milk solid in bovines. Lactose yield (LY, kg/d) is responsible for milk volume, whereas lactose percentage (LP) is thought to be more related to epithelial integrity and thus to udder health. There is a paucity of studies that have investigated lactose at the genomic level in dairy cows. This paper aimed to improve our knowledge on LP and LY, providing new insights into the significant genomic regions affecting these traits. A genome-wide association study for LP and LY was carried out in Fleckvieh cattle by using bulls' deregressed estimated breeding values of first lactation as pseudo-phenotypes. Heritabilities of first-lactation test-day LP and LY estimated using linear animal models were 0.38 and 0.25, respectively. A total of 2,854 bulls genotyped with a 54K SNP chip were available for the genome-wide association study; a linear mixed model approach was adopted for the analysis. The significant SNP of LP were scattered across the whole genome, with signals on chromosomes 1, 2, 3, 7, 12, 16, 18, 19, 20, 28, and 29; the top 4 significant SNP explained 4.90% of the LP genetic variance. The signals were mostly in regions or genes with involvement in molecular intra- or extracellular transport; for example, CDH5, RASGEF1C, ABCA6, and SLC35F3. A significant region within chromosome 20 was previously shown to affect mastitis or somatic cell score in cattle. As regards LY, the significant SNP were concentrated in fewer regions (chromosomes 6 and 14), related to mastitis/somatic cell score, immune response, and transport mechanisms. The 5 most significant SNP for LY explained 8.45% of genetic variance and more than one-quarter of this value has to be attributed to the variant within ADGRB1. Significant peaks in target regions remained even after adjustment for the 2 most significant variants previously detected on BTA6 and BTA14. The present study is a prelude for deeper investigations into the biological role of lactose for milk secretion and volume determination, stressing the connection with genes regulating intra- or extracellular trafficking and immune and inflammatory responses in dairy cows. Also, these results improve the knowledge on the relationship between lactose and udder health; they support the idea that LP and its derived traits are potential candidates as indicators of udder health in breeding programs aimed to enhance cows' resistance to mastitis.     

Short communication: Confirmation of candidate causative variants on milk composition and cheesemaking properties in Montbéliarde cows

In a previous study, we identified candidate causative variants located in 24 functional candidate genes for milk protein and fatty acid composition in Montbéliarde, Normande, and Holstein cows. We designed these variants on the custom part of the EuroG10K BeadChip (Illumina Inc., San Diego, CA), which is routinely used for genomic selection analyses in French dairy cattle. To validate the effects of these candidate variants on milk composition and to estimate their effects on cheesemaking properties, a genome-wide association study was performed on milk protein, fatty acid and mineral composition, as well as on 9 cheesemaking traits (3 laboratory cheese yields, 5 coagulation traits, and milk pH). All the traits were predicted from midinfrared spectra in the Montbéliarde cow population of the Franche-Comté region. A total of 194 candidate variants located in 24 genes and 17 genomic regions were imputed on 19,862 cows with phenotypes and genotyped with either the BovineSNP50 (Illumina Inc.) or the EuroG10K BeadChip. We then tested the effect of each SNP in a mixed linear model including random polygenic effects estimated with a genomic relationship matrix. We confirm here the effects of candidate causative variants located in 17 functional candidate genes on both cheesemaking properties and milk composition traits. In each candidate gene, we identified the most plausible causative variant: 4 are missense in the ALPL, SLC26A4, CSN3, and SCD genes, 7 are located in 5′UTR (AGPAT6), 3′ untranslated region (GPT), or upstream (CSN1S1, CSN1S2, PAEP, DGAT1, and PICALM) regions, and 6 are located in introns of the SLC37A1, MGST1, CSN2, BRI3BP, FASN, and ANKH genes.     

Short communication: Replication of genome-wide association studies for milk production traits in Chinese Holstein by an efficient rotated linear mixed model

Milk is regarded as an important nutrient for humans, and Chinese Holstein cows provide high-quality milk for billions of Chinese people. Therefore, detecting quantitative trait nucleotides (QTN) or candidate genes for milk production traits in Chinese Holstein is important. In this study, we performed genome-wide association studies (GWAS) in a Chinese Holstein population of 6,675 cows and 71,633 SNP using deregressed proofs (DRP) as phenotypes to replicate our previous study in a population of 1,815 cows and 39,163 SNP using estimated breeding values (EBV) as phenotypes. The associations between 3 milk production traits—milk yield (MY), fat percentage (FP), and protein percentage (PP)—and the SNP were determined by using an efficient rotated linear mixed model, which benefits from linear transformations of genomic estimated values and Eigen decomposition of the genomic relationship matrix algorithm. In total, we detected 94 SNP that were significantly associated with one or more milk production traits, including 7 SNP for MY, 76 for FP, and 36 for PP; 87% of these SNP were distributed across Bos taurus autosomes 14 and 20. In total, 83 SNP were found to be located within the reported quantitative trait loci (QTL) regions, and one novel segment (between 1.41 and 1.49 Mb) on chromosome 14 was significantly associated with FP, which could be an important candidate QTL region. In addition, the detected intervals were narrowed down from the reported regions harboring causal variants. The top significant SNP for the 3 traits was ARS-BFGL-NGS-4939, which is located within the DGAT1 gene. Five detected genes (CYHR1, FOXH1, OPLAH, PLEC, VPS28) have effects on all 3 traits. Our study provides a suite of QTN, candidate genes, and a novel QTL associated with milk production traits, and thus forms a solid basis for genomic selection and molecular breeding for milk production traits in Chinese Holstein.     

Genome-wide association study of normal and atypical progesterone profiles in Holstein-Friesian dairy cows

Female fertility has a major role in dairy production and affects the profitability of dairy cattle. The genetic progress obtained by traditional selection can be slow because of the low heritability of classical fertility traits. Endocrine fertility traits based on progesterone concentration in milk have higher heritability and more directly reflect the cow's own reproductive physiology. The aim of our study was to identify genomic regions for 7 endocrine fertility traits in dairy cows by performing a genome-wide association study with 54,000 SNP. The next step was to fine-map targeted genomic regions with significant SNP using imputed sequences to identify potential candidate genes associated with the normal and atypical progesterone profiles. The association between a SNP and a phenotype was assessed by a single SNP analysis, using a linear mixed model that included a random polygenic effect. Phenotypes and genotypes were available for 1,126 primiparous and multiparous Holstein-Friesian cows from research herds in Ireland, the Netherlands, Sweden, and the United Kingdom. In total, 44 significant SNP associated with 7 endocrine fertility traits were identified on Bos taurus autosome (BTA) 1–4, 6, 8–9, 11–12, 14–17, 19, 21–24, and 29. Three chromosomes, BTA8, BTA17, and BTA23, were imputed from 54,000 SNP genotypes to the whole-genome sequence level with Beagle version 4.1. The fine-mapping identified several significant associations with delayed cyclicity, cessation of cyclicity, commencement of luteal activity, and inter-ovulatory interval. These associations may contribute to an index of markers for genetic improvement of fertility. Several potential candidate genes reported to affect reproduction were also identified in the targeted genomic regions. However, due to high linkage disequilibrium, it was not possible to identify putative causal genes or polymorphisms for any of the regions.     

Genome-wide association analysis and pathways enrichment for lactation persistency in Canadian Holstein cattle

Lactation persistency (LP), defined as the rate of declining milk yield after milk peak, is an economically important trait for dairy cattle. Improving LP is considered a good alternative method for increasing overall milk production because it does not cause the negative energy balance and other health issues that cows experience during peak milk production. However, little is known about the biology of LP. A genome-wide association study (GWAS) and pathway enrichment were used to explore the genetic mechanisms underlying LP. The GWAS was performed using a univariate regression mixed linear model on LP data of 3,796 cows and 44,100 single nucleotide polymorphisms (SNP). Eight and 47 SNP were significantly and suggestively associated with LP, respectively. The 2 most important quantitative trait loci regions for LP were (1) a region from 106 to 108 Mb on Bos taurus autosome (BTA) 5, where the most significant SNP (ARS-BFGL-NGS-2399) was located and also formed a linkage disequilibrium block with 3 other SNP; and (2) a region from 29.3 to 31.3 Mb on BTA 20, which contained 3 significant SNP. Based on physical positions, MAN1C1, MAP3K5, HCN1, TSPAN9, MRPS30, TEX14, and CCL28 are potential candidate genes for LP because the significant SNP were located in their intronic regions. Enrichment analyses of a list of 536 genes in 0.5-Mb flanking regions of significant and suggestive SNP indicates that synthesis of milk components, regulation of cell apoptosis processes and insulin, and prolactin signaling pathways are important for LP. Upstream regulators relevant for LP positional candidate genes were prolactin (PRL), peroxisome proliferator-activated receptor gamma (PPARG), and Erb-B2 receptor tyrosine kinase 2 (ERBB2). Several networks related to cellular development, proliferation and death were significantly enriched for LP positional candidate genes. In conclusion, this study detected several SNP, genes, and interesting regions for fine mapping and validation of candidate genes and SNP for potential use in selection for improved LP. This study also provided further insights on the biology of LP which will help to prioritize selected candidate genes for functional validation and application.     

Phenotypic and genetic effects of season on milk production traits in dairy cattle in the Netherlands

Milk production systems in several countries show considerable differences between seasons. For example, in the Netherlands, cows are kept inside and fed silage in winter, whereas they are on pasture in summer. The differences between seasons affect milk yield and composition and might influence the genetic background of milk production traits. The objective of this study was to estimate phenotypic and genetic effects of season on milk production traits. For this purpose, 19,286 test-day milk production records of 1,800 first-parity Dutch Holstein-Frisian cows were available, and these cows were genotyped using a 50K SNP panel. Phenotypic effects of season were significant for all milk production traits. Effects of season were large for milk fat yield, fat content, and protein content. Genetic correlations between milk production traits in different seasons showed that genotype by season interaction effects were relatively small for most milk production traits. The genetic background of protein content and lactose content seems to be sensitive to seasonal effects. Furthermore, the genetic correlations between spring and autumn differed significantly from unity for almost all milk production traits. A genome-wide association study for genotype by season interaction identified chromosomal regions on BTA3, BTA14, BTA20, and BTA25 that showed genotype by season interaction effects, including a region containing DGAT1, which showed interaction effects for fat content and protein content.     

Short communication: Validation of somatic cell score–associated loci identified in a genome-wide association study in German Holstein cattle

Recently, we identified 6 genomic loci affecting daughter yield deviations (DYD) for somatic cell score (SCS) in a genome-wide association study (GWAS) performed with German Holstein bulls. In the current study, we tested if these loci were associated with SCS in cows using their own performance data. The study was performed with 1,412 German Holstein cows, of which 483 were daughters of 71 bulls that had been used in the GWAS. We tested 10 single nucleotide polymorphisms (SNP) representing 6 genomic regions that were associated with DYD for SCS in bulls. All tested SNP were significant in cows. Seven of them, located on Bos taurus autosomes (BTA) 6, 13, and 19, had the same direction of effect as those previously reported in the bull population. The most significant associations were detected on BTA6 and BTA19, accounting for 1.8% of the total genetic variance. The major allele of the 2 SNP on BTA6 and the minor allele of the 2 SNP on BTA19 were favorable for lower SCS. The differences between the homozygous genotype classes were up to 15,000 cells/mL. The verification of SNP associated with SCS in this study provides further evidence for the functional role of the linked genomic regions for immune response and contributes to identification of causative mutations. In particular, SNP with minor frequency of the favorable allele possess high potential to reduce SCS in German Holstein cattle by selection.     

Fine-mapping of BTA17 using imputed sequences for associations with de novo synthesized fatty acids in bovine milk

The aim of this study was to fine-map a genomic region associated with milk fatty acids (FA) on Bos taurus autosome (BTA) 17. This genomic region has been discovered with 50,000 (50k) single nucleotide polymorphisms (SNP) imputed to 777,000 (777k) SNP. In this study, high-density genotypes were imputed to whole-genome sequences level to identify candidate gene(s) associated with milk FA composition on BTA17. Phenotypes and genotypes were available for 1,640 cows sampled in winter, and for 1,581 cows sampled in summer. Phenotypes consisted of gas chromatography measurements in winter and in summer milk samples of 6 individual FA and the indicator of de novo synthesis, C6:0–C14:0. Genotypes consisted of imputed 777k SNP, and 89 sequenced ancestors of the population of genotyped cows. In addition, 450 whole-genome sequences from the 1,000 Bull Genome Consortium were available. Using 495 Holstein-Friesian sequences as a reference population, the 777k SNP genotypes of the cows were imputed to sequence level. We then applied single-variant analyses with an animal model, and identified thousands of significant associations with C6:0, C8:0, C10:0, C12:0, C14:0, and C6:0–C14:0. For C8:0 in summer milk samples, the genomic region located between 29 and 34 Mbp on BTA17 revealed a total of 646 significant associations. The most significant associations [–log₁₀(P-value) = 7.82] were 8 SNP in perfect linkage disequilibrium. After fitting one of these 8 SNP as a fixed effect in the model, and re-running the single-variant analyses, no further significant associations were found for any of the 6 FA or C6:0–C14:0. These findings suggest that one polymorphism underlying this QTL on BTA17 influences multiple de novo synthesized milk FA. Thirteen genes in the QTL region were identified and analyzed carefully. Six out of the 8 SNP that showed the strongest associations were located in the La ribonucleoprotein domain family, member 1B (LARP1B) gene, and we suggest LARP1B as a primary candidate gene. Another gene of interest for this QTL region might be PKL4. None of these suggested candidate genes have previously been associated with milk fat synthesis or milk FA composition.     

Genome-wide association study for selected cheese-making properties in Dual-Purpose Belgian Blue cows

This study aimed to estimate genetic parameters and identify genomic region(s) associated with selected cheese-making properties (CMP) in Dual-Purpose Belgian Blue (DPBB) cows. Edited data were 46,301 test-day records of milk yield, fat percentage, protein percentage, casein percentage, milk calcium content (CC), coagulation time (CT), curd firmness after 30 min from rennet addition (a30), and milk titratable acidity (MTA) collected from 2014 to 2020 on 4,077 first-parity (26,027 test-day records), and 3,258 second-parity DPBB cows (20,274 test-day records) distributed in 124 herds in the Walloon Region of Belgium. Data of 28,266 SNP, located on 29 Bos taurus autosomes (BTA) of 1,699 animals were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 25 consecutive SNPs (with an average size of ～2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Heritability estimates for the included CMP ranged from 0.19 (CC) to 0.50 (MTA), and 0.24 (CC) to 0.41 (MTA) in the first and second parity, respectively. The genetic correlation estimated between CT and a30 varied from ?0.61 to ?0.41 and from ?0.55 to ?0.38 in the first and second lactations, respectively. Negative genetic correlations were found between CT and milk yield and composition, while those estimated between curd firmness and milk composition were positive. Genome-wide association analyses results identified 4 genomic regions (BTA1, BTA3, BTA7, and BTA11) associated with the considered CMP. The identified genomic regions showed contrasting results between parities and among the different stages of each parity. It suggests that different sets of candidate genes underlie the phenotypic expression of the considered CMP between parities and lactation stages of each parity. The findings of this study can be used for future implementation and use of genomic evaluation to improve the cheese-making traits in DPBB cows.     

Longitudinal genome-wide association studies of milk production traits in Holstein cattle using whole-genome sequence data imputed from medium-density chip data

Longitudinal traits, such as milk production traits in dairy cattle, are featured by having phenotypic values at multiple time points, which change dynamically over time. In this study, we first imputed SNP chip (50–100K) data to whole-genome sequence (WGS) data in a Chinese Holstein population consisting of 6,470 cows. The imputation accuracies were 0.88 to 0.97 on average after quality control. We then performed longitudinal GWAS in this population based on a random regression test-day model using the imputed WGS data. The longitudinal GWAS revealed 16, 39, and 75 quantitative trait locus regions associated with milk yield, fat percentage, and protein percentage, respectively. We estimated the 95% confidence intervals (CI) for these quantitative trait locus regions using the logP drop method and identified 581 genes involved in these CI. Further, we focused on the CI that covered or overlapped with only 1 gene or the CI that contained an extremely significant top SNP. Twenty-eight candidate genes were identified in these CI. Most of them have been reported in the literature to be associated with milk production traits, such as DGAT1, HSF1, MGST1, GHR, ABCG2, ADCK5, and CSN1S1. Among the unreported novel genes, some also showed good potential as candidate genes, such as CCSER1, CUX2, SNTB1, RGS7, OSR2, and STK3, and are worth being further investigated. Our study provided not only new insights into the candidate genes for milk production traits, but also a general framework for longitudinal GWAS based on random regression test-day model using WGS data.     

Fine mapping of a quantitative trait locus for bovine milk fat composition on Bos taurus autosome 19

A major quantitative trait locus (QTL) for milk fat content and fatty acids in both milk and adipose tissue has been detected on Bos taurus autosome 19 (BTA19) in several cattle breeds. The objective of this study was to refine the location of the QTL on BTA19 for bovine milk fat composition using a denser set of markers. Opportunities for fine mapping were provided by imputation from 50,000 genotyped single nucleotide polymorphisms (SNP) toward a high-density SNP panel with up to 777,000 SNP. The QTL region was narrowed down to a linkage disequilibrium block formed by 22 SNP covering 85,007 bp, from 51,303,322 to 51,388,329 bp on BTA19. This linkage disequilibrium block contained 2 genes: coiled-coil domain containing 57 (CCDC57) and fatty acid synthase (FASN). The gene CCDC57 is minimally characterized and has not been associated with bovine milk fat previously, but is expressed in the mammary gland. The gene FASN has been associated with bovine milk fat and fat in adipose tissue before. This gene is a likely candidate for the QTL on BTA19 because of its involvement in de novo fat synthesis. Future studies using sequence data of both CCDC57 and FASN, and eventually functional studies, will have to be pursued to assign the causal variant(s).     

Genome-wide association study identifies loci influencing natural antibody titers in milk of Dutch Holstein-Friesian cattle

Natural antibodies (NAb) are produced without any antigenic stimulation as a part of the innate immune system and provide a first line of defense against pathogens. Hence, they may be a useful trait when estimating an animal's potential immune competence and in selection for disease resistance. The aim of this study was to identify genomic regions associated with different NAb traits in milk and potentially describe candidate genes. Milk samples from 1,695 first-lactation Holstein Friesian cows with titer measurements for keyhole limpet hemocyanin, lipopolysaccharide, lipoteichoic acid, and peptidoglycan-binding total NAb and isotypes IgG₁, IgM, and IgA were used. Genome-wide association study analyses were performed using imputed 777K SNP genotypes, accounting for relationships using pedigree information. Functional enrichment analysis was performed on the significantly associated genomic regions to look for candidate genes. For IgM NAb, significant associations (false discovery rate <0.05) were found on Bos taurus autosome (BTA) 17, 18, and 21 with candidate genes related to immunoglobulin structure and early B cell development. For IgG₁, associations were found on BTA3, and we confirmed a quantitative trait loci on BTA21 previously reported for IgG NAb in serum. Our results provide new insights into the regulation of milk NAb that will help unravel the complex relationship between milk immunoglobulins and disease resistance in dairy cattle.     

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.     

Genome-wide scan for bovine milk-fat composition. I. Quantitative trait loci for short- and medium-chain fatty acids

A genome-wide scan was performed to identify quantitative trait loci (QTL) for short- and medium-chain fatty acids (expressed in wt/wt %). Milk samples were available from 1,905 cows from 398 commercial herds in the Netherlands, and milk-fat composition was measured by gas chromatography. DNA was available from 7 of the paternal half-sib families: 849 cows and their 7 sires. A genetic map was constructed comprising 1,341 SNP and 2,829 cM, with an average information content of 0.83. Multimarker interval mapping was used in an across-family regression on corrected phenotypes for the 7 half-sib families. Four QTL were found: on Bos taurus autosome (BTA) 6, a QTL was identified for C6:0 and C8:0; on BTA14, a QTL was identified for fat percentage, all odd-chain fatty acids, and C14:0, C16:0, C16:1, and their unsaturation indices; on BTA19, a QTL affected C14:0; and on BTA26, a QTL was identified for the monounsaturated fatty acids and their unsaturation indices. The QTL explained 3 to 19% of phenotypic variance. Furthermore, 49 traits with suggestive evidence for linkage were found on 21 chromosomes. Additional analyses revealed that the QTL on BTA14 was most likely caused by a mutation in DGAT1, whereas the QTL on BTA26 was most likely caused by a mutation in the SCD1 gene. Quantitative trait loci that affect specific fatty acids might increase the understanding of physiological processes regarding fat synthesis and the position of the causal genes.     

Gene mapping,gene-set analysis,and genomic prediction of postpartum blood calcium in Holstein cows

The onset of lactation results in a sudden irreversible loss of Ca for colostrum and milk synthesis. Some cows are unable to quickly adapt to this demand and succumb to clinical hypocalcemia, whereas a larger proportion of cows develop subclinical hypocalcemia that predisposes them to other peripartum diseases. The objective of this study was to perform a comprehensive genomic analysis of blood total Ca concentration in periparturient Holstein cows. We first performed a genomic scan and a subsequent gene-set analysis to identify candidate genes, biological pathways, and molecular mechanisms affecting postpartum Ca concentration. Then, we assessed the prediction of postpartum Ca concentration using genomic information. Data consisted of 7,691 records of plasma or serum concentrations of Ca measured in the first, second, and third day after parturition of 959 primiparous and 1,615 multiparous cows that calved between December 2015 and June 2020 in 2 dairy herds. All cows were genotyped with 80k SNPs. The statistical model included lactation (1 to 5+), calf category (male, females, twins), and day as fixed effects, and season-treatment-experiment, animal, and permanent environmental as random effects. Model predictive ability was evaluated using 10-fold cross-validation. Heritability and repeatability estimates were 0.083 (standard error = 0.017) and 0.444 (standard error = 0.028). The association mapping identified 2 major regions located on Bos taurus autosome (BTA)6 and BTA16 that explained 1.2% and 0.7% of additive genetic variance of Ca concentration, respectively. Interestingly, the region on BTA6 harbors the GC gene, which encodes the vitamin D binding protein, and the region on BTA16 harbors LRRC38, which is actively involved in K transport. Other sizable peaks were identified on BTA5, BTA2, BTA7, BTA14, and BTA9. These regions harbor genes associated with Ca channels (CACNA1S, CRACR2A), K channels (KCNK9), bone remodeling (LRP6), and milk production (SOCS2). The gene-set analysis revealed terms related to vitamin transport, calcium ion transport, calcium ion binding, and calcium signaling. Genomic predictions of phenotypic and genomic estimated breeding values of Ca concentration yielded predictive correlations up to 0.50 and 0.15, respectively. Overall, the present study contributes to a better understanding of the genetic basis of postpartum blood Ca concentration in Holstein cows. In addition, the findings may contribute to the development of novel selection and management strategies for reducing periparturient hypocalcemia in dairy cattle.     

Whole-genome association study for milk protein composition in dairy cattle

Our objective was to perform a genome-wide association study for content in bovine milk of α_S1-casein (α_S1-CN), α_S2-casein (α_S2-CN), β-casein (β-CN), κ-casein (κ-CN), α-lactalbumin (α-LA), β-lactoglobulin (β-LG), casein index, protein percentage, and protein yield using a 50K single nucleotide polymorphism (SNP) chip. In total, 1,713 Dutch Holstein-Friesian cows were genotyped for 50,228 SNP and a 2-step association study was performed. The first step involved a general linear model and the second step used a mixed model accounting for all family relationships. Associations with milk protein content and composition were detected on 20 bovine autosomes. The main genomic regions associated with milk protein composition or protein percentage were found on chromosomes 5, 6, 11, and 14. The number of chromosomal regions showing significant (false discovery rate <0.01) effects ranged from 3 for β-CN and 3 for β-LG to 12 for α_S2-CN. A genomic region on Bos taurus autosome (BTA) 6 was significantly associated with all 6 major milk proteins, and a genomic region on BTA 11 was significantly associated with the 4 caseins and β-LG. In addition, regions were detected that only showed a significant effect on one of the milk protein fractions: regions on BTA 13 and 22 with effects on α_S1-CN; regions on BTA 1, 9, 10, 17, 19, and 28 with effects on α_S2-CN; a region on BTA 6 with an effect on β-CN; regions on BTA 13 and 21 with effects on κ-CN; regions on BTA 1, 5, 9, 16, 17, and 26 with effects on α-LA; and a region on BTA 24 with an effect on β-LG. The proportion of genetic variance explained by the SNP showing the strongest association in each of these genomic regions ranged from <1% for α_S1-CN on BTA 22 to almost 100% for casein index on BTA 11. Variation associated with regions on BTA 6, 11, and 14 could in large part but not completely be explained by known protein variants of β-CN (BTA 6), κ-CN (BTA 6), and β-LG (BTA 11) or DGAT1 variants (BTA 14). Our results indicate 3 regions with major effects on milk protein composition, in addition to several regions with smaller effects involved in the regulation of milk protein composition.     

Genome-wide association study and functional analyses for clinical and subclinical ketosis in Holstein cattle

Ketosis is one of the most frequent metabolic diseases in high-yielding dairy cows and is characterized by high concentrations of ketone bodies in blood, urine, and milk, causing high economic losses. The search for polymorphic genes, whose alleles have different effects on resistance to developing the disease, is of extreme importance to help select less susceptible animals. The aims of this study were to identify genomic regions associated with clinical and subclinical ketosis (β-hydroxybutyrate concentration) in North American Holstein dairy cattle and to investigate these regions to identify candidate genes and metabolic pathways associated with these traits. To achieve this, a GWAS was performed for 4 traits: clinical ketosis lactation 1, clinical ketosis lactation 2 to 5, subclinical ketosis lactation 1, and subclinical ketosis lactation 2 to 5. The estimated breeding values from 77,277 cows and 7,704 bulls were deregressed and used as pseudophenotypes in the GWAS. The top-20 genomic regions explaining the largest proportion of the genetic variance were investigated for putative genes associated with the traits through functional analyses. Regions of interest were identified on chromosomes 2, 5, and 6 for clinical ketosis lactation 1; 3, 6, and 7 for clinical ketosis lactation 2 to 5; 1, 2, and 12 for subclinical ketosis lactation 1; and 20, 11, and 25 for subclinical ketosis lactation 2 to 5. The highlighted genes potentially related to clinical and subclinical ketosis included ACAT2 and IGF1. Enrichment analysis of the list of candidate genes for clinical and subclinical ketosis showed molecular functions and biological processes involved in fatty acid metabolism, lipid metabolism, and inflammatory response in dairy cattle. Several genomic regions and SNPs related to susceptibility to ketosis in dairy cattle that were previously described in other studies were confirmed. The novel genomic regions identified in this study aid to characterize the most important genes and pathways that explain the susceptibility to clinical and subclinical ketosis in dairy cattle.