The objective of this study was to investigate genetic variation and genotype by environment (G × E) interactions for fertility (including age at first calving and calving interval), somatic cell score (SCS), and milk production traits for Iranian Holsteins. Different environments were defined based on the climatic zones (cold, semi-cold, and moderate) and considering the herd location. Data were collected between 2003 and 2018 by the National Animal Breeding Center of Iran (Karaj). Variance and covariance components and genetic correlations were estimated using 2 different models, which were analyzed using Bayesian methods. For both models, performance of traits in each climate were considered as different traits. Fertility traits were analyzed using a trivariate model. Furthermore, SCS and production traits were analyzed using trivariate random regression models (records in different climate zones considered as different traits). For the fertility traits, the largest estimates of heritability were observed in cold climate. Fertility performance was always better in cold environment. Genetic correlations between climatic zones ranged from 0.85 to 0.94. For daily measurements of SCS and production traits, heritability ranged from 0.031 to 0.037 and 0.069 to 0.209, respectively. Genetic variances were the highest in the semi-cold and moderate climates for the SCS and production traits, respectively. Furthermore, across the studied climates, 305-d genetic correlation ranged from 0.756 to 0.884 for SCS and from 0.925 to 0.957 for the production traits. The structure of genetic correlation within each climate indicated a negative correlation between early and late lactation for SCS, especially in the cold climate and for milk production in the moderate climate. For fat percentage, in all climatic zones, the lowest genetic correlations were observed between early and mid-lactation. In addition, for protein production in the cold climate, a negative correlation was observed between early and late lactation. Results indicated heterogeneous variance components for all the studied traits across various climatic zones. Estimated genetic correlations for SCS revealed that the genetic expression of animals may vary by climatic zone. Results indicated the existence of G × E interaction due to the climatic condition, only for SCS. Therefore, in Iranian Holsteins, the effect of G × E interactions should not be neglected, especially for SCS, as different sires might be optimal for use in different climatic zones.