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 ( ) were 0.07 ± 0.02, 0.04 ± 0.01, 0.06 ± 0.01, 0.05 ± 0.01, 0.5 ± 0.01, 0.09 ± 0.02, 0.07 ± 0.02 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. For F:P, in addition to the maternal imprinting effect, a significant paternal imprinting component was also detected with a 7% contribution to phenotypic variance of F:P. Estimates of direct heritability ( ) were 0.29 ± 0.02, 0.17 ± 0.01, 0.22 ± 0.02, 0.11 ± 0.01, 0.18 ± 0.02, 0.22 ± 0.02, 0.15 ± 0.04 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. Maternal imprinting correlations (rmi) were in a wide range between -0.75 ± 0.15 (P%-SCS) and 0.95 ± 0.11 (MY-ADM). Additive genetic correlations (ra) ranged between -0.54 ± 0.05 (MY-P%) and 0.99 ± 0.01 (MY-ADM) and phenotypic correlations (rp) ranged from -0.30 ± 0.01 (MY-F%) to 0.93 ± 0.01 (MY-ADM). The Spearman's correlation between additive breeding values including and excluding imprinting effects deviated from unity especially for top-ranked animals implying re-ranking of top animals following the inclusion of imprinting effects in the model. Since including imprinting effects in the model resulted in better data fit and re-ranking of top animals, including these effects in the genetic evaluation models for milk production traits was recommended.