مریم دوستی توسه چالکی

 In dairy cows, milk yield changes over the lactation period. Lactation curve modelling could be of importance from nutritional management and genetic selection point of view. The trajectory shape of the lactation curve has inclining and declining slopes which determine total amount of milk yield during the lactation. A great number of linear and nonlinear mathematical models have been so far utilized to describe the shape of the lactation curve in dairy cows. Among the models, mechanistic functions are considered to be more accurate than empirical models in terms of taking account of biological mechanisms undertaken in mammary gland to produce milk. To our knowledge, no research has been carried out to use a complex mechanistic model for describing the lactation curve of Iranian dairy cows. Based on this, the present research aimed to apply a five-parameter nonlinear mechanistic model for mathematical description of the shape of the lactation curve in Iranian first-parity dairy cows.

Animal Breeding Centre of Iran provided the data used in this research. Initial data in separate files were merged and edited by the SPSS software. Final data comprised 5596039 milk test day records from 821153 first-parity cows distributed in 579 herds of 26 provinces over the country. The cows were the progeny of 7957 sires and 530394 dams and calved between 1996 and 2020. The cows were categorized into two groups based on the age of first calving: <=25 or >25 months. Also, two groups of cows were created based on the proportion of Holstein gene inheritance: <100% (grade cows) or 100% (pure Holstein cows). A five-parameter nonlinear mechanistic model was applied to mathematically describe the shape of the lactation curve. The fitted mechanistic model had five parameters including MSmax (maximum milk secretion potential of the lactation), GR (relative proliferation rate of secretory cell number during early lactation), MSLmax (maximum secretion loss), NOD (proportion of parenchyma cell dead at parturition), and DR (relative decline rate in cell number) which their estimates were obtained for different seasons and ages of calving, as well as different genotypes (grade and pure Holstein) using NLIN procedure of the SAS software. Based on the estimated parameters of the model in each year of calving, phenotypic trends were calculated using the SPSS software.

The results indicated that Iranian first-parity dairy cows reach to peak of milk production during the third month of the lactation curve. Based on standard deviation as well as coefficient of variation (CV) of milk test day milk records, maximum variation was observed in the last month of the lactation period. The findings also indicated that cows calved in autumn had the greatest MSmax and GR while cows calved in spring and summer had the lowest magnitude for these parameters (P<0.05). The minimum magnitude of MSLmax and NOD were detected for the cows calved in summer (P<0.05). Regarding the DR parameter of the model, cows calved in spring were observed to have a minimum value. Maximum MSmax, Gr, and DR values ​​were determined for cows calving up to 25 months of age, while maximum MSLmax and NOD values ​​were observed for cows of later ages (P<0.05). Compared to pure Holstein cows, grade cows had higher MSLmax and NOD values, but other parameters were found to be greater in the pure Holstein cows (P<0.05), suggesting that the maximum milk secretion potential during lactation in pure Holstein cows is expected to be higher than that of grade cows. Annual phenotypic change trends were found to be 0.376 kg/year (R2=0.9), 0.00009869 (R2=071), -0.146 kg/year (R2=0.59), -0.011 (R2=0.85), and 0.001 (R2=0.85) for MSmax, GR, MSLmax, NOD, and DR, respectively (P<0.0001).     

This study found that the parameters of the mechanistic model fitted to the milk test day records of Iranian primiparous dairy cows are significantly influenced by the age and season of calving as well as the genotype of the cow. Cows calving in summer are expected to have more persistency. Moreover, pure Holsteins are more persistent than grade cows. Positive and negative annual trends have been detected for MSmax and MSLmax parameters, respectively, during 1996-2020 indicating a favorable increase in maximum milk secretion potential of the lactation and also a favorable decrease in maximum secretion loss over that period.