Breast cancer is a heterogeneous disease that is affected by ethnicity of patients. According to hormone receptor status and gene expression profiling, breast cancers are classified into four molecular subtypes, each showing distinct clinical behavior. Lack of sufficient data on molecular subtypes of breast cancer in Iran, prompted us to investigate the prevalence and the clinicopathological features of each subtype among Iranian women. A total of 428 women diagnosed with breast cancer from 2002 to 2011 were included and categorized into four molecular subtypes using immunohistochemistry. Prevalence of each subtype and its association with patients' demographics and tumor characteristics, such as size, grade, lymph-node involvement and vascular invasion, were investigated using Chi-square, analysis of variance and multivariate logistic regression. Luminal A was the most common molecular subtype (63.8%) followed by Luminal B (8.4%), basal-like (15.9%) and HER-2 (11.9%). Basal-like and HER-2 subtypes were mostly of higher grades while luminal A tumors were more of grade 1 (P<0.001). Vascular invasion was more prevalent in HER-2 subtype, and HER-2 positive tumors were significantly associated with vascular invasion (P=0.013). Using muti-variate analysis, tumor size greater than 5 cm and vascular invasion were significant predictors of 3 or more nodal metastases. Breast cancer was most commonly diagnosed in women around 50 years of age and the majority of patients had lymph node metastasis at the time of diagnosis. This points to the necessity for devising an efficient screening program for breast cancer in Iran. Further, prospective surveys are suggested to evaluate prognosis of different subtypes in Iranian patients.