Search

Background: Genome-wide association studies (GWAS) were performed at the sequence level to identify candidate mutations that affect the expression of six major milk proteins in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) dairy cattle. Whey protein (α-lactalbumin and β-lactoglobulin) and casein (αs1, αs2, β, and κ) contents were estimated by mid-infrared (MIR) spectrometry, with medium to high accuracy (0.59 ≤ R 2  ≤ 0.92), for 848,068 test-day milk samples from 156,660 cows in the first three lactations. Milk composition was evaluated as average test-day measurements adjusted for environmental effects. Next, we genotyped a subset of 8080 cows (2967 MON, 2737 NOR, and 2306 HOL) with the BovineSNP50 Beadchip. For each breed, genotypes were first imputed to high-density (HD) using HD single nucleotide polymorphisms (SNPs) genotypes of 522 MON, 546 NOR, and 776 HOL bulls. The resulting HD SNP genotypes were subsequently imputed to the sequence level using 27 million high-quality sequence variants selected from Run4 of the 1000 Bull Genomes consortium (1147 bulls). Within-breed, multi-breed, and conditional GWAS were performed., Results: Thirty-four distinct genomic regions were identified. Three regions on chromosomes 6, 11, and 20 had very significant effects on milk composition and were shared across the three breeds. Other significant effects, which partially overlapped across breeds, were found on almost all the autosomes. Multi-breed analyses provided a larger number of significant genomic regions with smaller confidence intervals than within-breed analyses. Combinations of within-breed, multi-breed, and conditional analyses led to the identification of putative causative variants in several candidate genes that presented significant protein-protein interactions enrichment, including those with previously described effects on milk composition (SLC37A1, MGST1, ABCG2, CSN1S1, CSN2, CSN1S2, CSN3, PAEP, DGAT1, AGPAT6) and those with effects reported for the first time here (ALPL, ANKH, PICALM)., Conclusions: GWAS applied to fine-scale phenotypes, multiple breeds, and whole-genome sequences seems to be effective to identify candidate gene variants. However, although we identified functional links between some candidate genes and milk phenotypes, the causality between candidate variants and milk protein composition remains to be demonstrated. Nevertheless, the identification of potential causative mutations that underlie milk protein composition may have immediate applications for improvements in cheese-making.

In the context of the PhénoFinLait project, a genome-wide analysis was performed to detect quantitative trait loci (QTL) that affect milk protein composition estimated using mid-infrared spectrometry in the Montbéliarde (MO), Normande (NO), and Holstein (HO) French dairy cattle breeds. The 6 main milk proteins (α-lactalbumin, β-lactoglobulin, and αS1-, αS2-, β-, and κ-caseins) expressed as grams per 100g of milk (% of milk) or as grams per 100g of protein (% of protein) were estimated in 848,068 test-day milk samples from 156,660 cows. Genotyping was performed for 2,773 MO, 2,673 NO, and 2,208 HO cows using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). Individual test-day records were adjusted for environmental effects and then averaged per cow to define the phenotypes analyzed. Quantitative trait loci detection was performed within each breed using a linkage disequilibrium and linkage analysis approach. A total of 39 genomic regions distributed on 20 of the 29 Bos taurus autosomes (BTA) were significantly associated with milk protein composition at a genome-wide level of significance in at least 1 of the 3 breeds. The 9 most significant QTL were located on BTA2 (133 Mbp), BTA6 (38, 47, and 87 Mbp), BTA11 (103 Mbp), BTA14 (1.8 Mbp), BTA20 (32 and 58 Mbp), and BTA29 (8 Mbp). The BTA6 (87 Mbp), BTA11, and BTA20 (58 Mbp) QTL were found in all 3 breeds, and they had highly significant effects on κ-casein, β-lactoglobulin, and α-lactalbumin, expressed as a percentage of protein, respectively. Each of these QTL explained between 13% (BTA14) and 51% (BTA11) of the genetic variance of the trait. Many other QTL regions were also identified in at least one breed. They were located on 14 additional chromosomes (1, 3, 4, 5, 7, 15, 17, 19, 21, 22, 24, 25, 26, and 27), and they explained 2 to 8% of the genetic variance of 1 or more protein composition traits. Concordance analyses, performed between QTL status and sequence-derived polymorphisms from 13 bulls, revealed previously known causal polymorphisms in LGB (BTA11) and GHR (BTA20 at 32 Mbp) and excluded some other previously described mutations. These results constitute a first step in identifying causal mutations and using routinely collected mid-infrared predictions in future genomic selection programs to improve bovine milk protein composition., (Copyright © 2016 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

In 2010, a routine genetic evaluation on occurrence of clinical mastitis in three main dairy cattle breeds-- Montbéliarde (MO), Normande (NO) and Holstein (HO)--was implemented in France. Records were clinical mastitis events reported by farmers to milk recording technicians and the analyzed trait was the binary variable describing the occurrence of a mastitis case within the first 150 days of the first three lactations. Genetic parameters of clinical mastitis were estimated for the three breeds. Low heritability estimates were found: between 2% and 4% depending on the breed. Despite its low heritability, the trait exhibits genetic variation so efficient genetic improvement is possible. Genetic correlations with other traits were estimated, showing large correlations (often>0.50, in absolute value) between clinical mastitis and somatic cell score (SCS), longevity and some udder traits. Correlation with milk yield was moderate and unfavorable (ρ=0.26 to 0.30). High milking speed was genetically associated with less mastitis in MO (ρ=-0.14) but with more mastitis in HO (ρ=0.18). A two-step approach was implemented for routine evaluation: first, a univariate evaluation based on a linear animal model with permanent environment effect led to pre-adjusted records (defined as records corrected for all non-genetic effects) and associated weights. These data were then combined with similar pre-adjusted records for others traits in a multiple trait BLUP animal model. The combined breeding values for clinical mastitis obtained are the official (published) ones. Mastitis estimated breeding values (EBV) were then combined with SCSs EBV into an udder health index, which receives a weight of 14.5% to 18.5% in the French total merit index (ISU) of the three breeds. Interbull genetic correlations for mastitis occurrence were very high (ρ=0.94) with Nordic countries, where much stricter recording systems exist reflecting a satisfactory quality of phenotypes as reported by the farmers. They were lower (around 0.80) with countries supplying SCS as a proxy for the international evaluation on clinical mastitis.