Bayesian analysis of the association between casein complex haplotype variants and milk yield, composition, and curve shape parameters in murciano-granadina goats

Considering casein haplotype variants rather than SNPs may maximize the understanding of heritable mechanisms and their implication on the expression of functional traits related to milk production. Effects of casein complex haplotypes on milk yield, milk composition, and curve shape parameters were used using a Bayesian inference for ANOVA. We identified 48 single nucleotide polymorphisms (SNPs) present in the casein complex of 159 unrelated individuals of diverse ancestry, which were organized into 86 haplotypes. The Ali and Schaeffer model was chosen as the best fitting model for milk yield (Kg), protein, fat, dry matter, and lactose (%), while parabolic yield-density was chosen as the best fitting model for somatic cells count (SCC × 103 sc/mL). Peak and persistence for all traits were computed respectively. Statistically significant differences (p < 0.05) were found for milk yield and components. However, no significant difference was found for any curve shape parameter except for protein percentage peak. Those haplotypes for which higher milk yields were reported were the ones that had higher percentages for protein, fat, dry matter, and lactose, while the opposite trend was described by somatic cells counts. Conclusively, casein complex haplotypes can be considered in selection strategies for economically important traits in dairy goats

Software-automatized individual lactation model fitting, peak and persistence and Bayesian criteria comparison for milk yield genetic studies in Murciano-Granadina goats

SPSS model syntax was defined and used to evaluate the individual performance of 49 linear and non-linear models to fit the lactation curve of 159 Murciano-Granadina goats selected for genotyping analyses. Lactation curve shape, peak and persistence were evaluated for each model using 3107 milk yield controls with an average of 3.78 ± 2.05 lactations per goat. Best fit (Adjusted R2) values (0.47) were reached by the five-parameter logarithmic model of Ali and Schaeffer. Three main possibilities were detected: non-fitting (did not converge), standard (Adjusted R2 over 75%) and atypical curves (Adjusted R2 below 75%). All the goats fitted for 38 models. The ability to fit different possible functional forms for each goat, which progressively increased with the number of parameters comprised in each model, translated into a higher sensitivity to explaining curve shape individual variability. However, for models for which all goats fitted, only moderate increases in explanatory and predictive potential (AIC, AICc or BIC) were found. The Ali and Schaeffer model reported the best fitting results to study the genetic variability behind goat milk yield and perhaps enhance the evaluation of curve parameters as trustable future selection criteria to face the future challenges offered by the goat dairy industry

Non-parametric association analysis of additive and dominance effects of casein complex SNPs on milk content and quality in Murciano-Granadina goats

Goat milk casein proteins (αS1, αS2, β and κ) are encoded by four loci (CSN1S1, CSN1S2, CSN2 and CSN3, respectively) clustered within 250 kb in chromosome 6. In this study, 159 Murciano-Granadina goats were genotyped for 48 SNPs within the entire casein region. Phenotypes on milk yield and components were obtained from 2,594 dairy registries. Additive and dominance effects on milk composition and quality were studied using non-parametric tests and principal component analysis to prevent SNPs multicollinearity. Two deletions in exon 4 (CSN1S1 and CSN3), one in exon 7 (CSN2) and one in exon 15 (CSN1S2) have been found at frequencies ranging from 0.12 to 0.50. Bonferroni-corrected significant SNP additive and dominance effects were found for milk yield, fat, protein, dry matter and lactose, and somatic cells. Exons 15 and 7 were significantly associated with milk yield and components except for lactose and somatic cells, while exon 4 was significantly associated with milk yield and components except for protein and dry matter. SNPs' associations with somatic cells were less frequent and weaker than those with milk yield and components. As caseins increase, somatic cells decrease, reducing milk enzymatic activity and consumption suitability. Hence, including molecular information in breeding schemes may promote production efficiency, as selecting against undesirable alleles could prevent the compromises derived from their dominance effects