<p>The aim was to estimate the relative contribution of imprinting effects from both paternal and maternal sides to phenotypic variation in milk production traits including 305 days milk yield (MY), average daily milk production (ADM), fat percentage (F%), protein percentage (P%), 305 days fat yield (FY), 305 days protein yield (PY), ratio of fat percentage to protein percentage (F:P) and somatic cell score (SCS) in Iranian Holstein cows. To do this, each trait was analysed with a series of four animal models, which were identical for fixed and additive genetic effects but differed for combinations of paternal and maternal imprinting effects. The log-likelihood ratio test (LRT) and Akaike's information criteria (AIC) were used to select the best model for each trait. Correlations between traits due to additive and imprinting effects were estimated by bivariate analyses. For all traits studied, fitting the imprinting effect led to a better data fit. Also, it resulted in a noticeable decrease in additive genetic variance from 8% (SCS) to 28% (F:P). A significant maternal imprinting effect was detected on all traits studied. Estimates of maternal imprinting heritability (<span></span><math>