Longitudinal data analysis of polymorphisms in the κ-casein and β-lactoglobulin genes shows differential effects along the trajectory of the lactation curve in tropical dairy goats.

The κ-casein (CSN-3) and β-lactoglobulin (BLG) genes are extensively polymorphic in ruminants. Several association studies have estimated the effects of polymorphisms in these genes on milk yield, milk composition, and cheese-manufacturing properties. Usually, these results are based on production integrated over the lactation curve or on cross-sectional studies at specific days in milk (DIM). However, as differential expression of milk protein genes occurs over lactation, the effect of the polymorphisms may change over time. In this study, we fitted a mixed-effects regression model to test-day records of milk yield and milk quality traits (fat, protein, and total solids yields) from Colombian tropical dairy goats. We used the well-characterized A/B polymorphisms in the CSN-3 and BLG genes. We argued that this approach provided more efficient estimators than cross-sectional designs, given the same number and pattern of observations, and allowed exclusion of between-subject variation from model error. The BLG genotype AA showed a greater performance than the BB genotype for all traits along the whole lactation curve, whereas the heterozygote showed an intermediate performance. We observed no such constant pattern for the CSN-3 gene between the AA homozygote and the heterozygote (the BB genotype was absent from the sample). The differences among the genotypic effects of the BLG and the CSN-3 polymorphisms were statistically significant during peak and mid lactation (around 40-160 DIM) for the BLG gene and only for mid lactation (80-145 DIM) for the CSN-3 gene. We also estimated the additive and dominant effects of the BLG locus. The locus showed a statistically significant additive behavior along the whole lactation trajectory for all quality traits, whereas for milk yield the effect was not significant at later stages. In turn, we detected a statistically significant dominance effect only for fat yield in the early and peak stages of lactation (at about 1-45 DIM). The longitudinal analysis of test-day records allowed us to estimate the differential effects of polymorphisms along the lactation curve, pointing toward stages that could be affected by the gene.