This study was performed to explore the relationship between morphological characters and nrDNA ITS data belonging to two subspecies of Primula vulgaris Huds. (Primulaceae) along the altitudinal gradient over NE Anatolia. There were no significant molecular differences between P. vulgaris subsp. vulgaris and P. vulgaris subsp. sibthorpii (Hoffmanns) W.W.Sm. & Forrest related to the flower colour and altitudinal gradient. However, Maximum parsimony (MP) analysis showed that all the examined populations fall in two distinct clusters confirming that the analyzed sample corresponds to two subspecies. It was also determined that the presence or absence of leaf pubescence does not vary along the altitudinal gradient within each subspecies.