This study addresses the phylogenetic relationship between Pistacia species by amplified fragment length polymorphism (AFLP). The plant materials of this study consisted of a total of 44 accessions belonging to P. vera, P. eurycarpa, P. khinjuk, all subspecies of P. atlantica (atlantica, mutica, kurdica and cabulica), three unknown genotypes and three accessions, proposed to be hybrid from P. eurycarpa x P. atlantica. The accessions were from Iran, Turkey, USA and Syria. Six AFLP primer combinations produced a total of 475 fragments, with average of 79.16 fragments per primer pair, of which, 336 bands were polymorphic. Unweighted pair group method based on arithmetic average (UPGMA) analysis was performed on jaccard's similarity coefficient matrix and also average similarity of each species. According to the results, two main clusters were developed and P. vera, P. eurycarpa, P. atlantica (subsp. atlantica, kurdica, mutica, cabulica) and the hybrid genotypes located in the first main cluster. P. khinjuk accessions from Iran and USA localized in second main cluster. The hybrid accessions located between eurycarpa and atlantica species and their hybrid nature between these two species were confirmed. One of the unknown accessions clustered with the hybrid ones and the two other were grouped closely with P. Khinjuk. According to this study, the closest species to P. vera was Eurycarpa group, followed by P. atlantica. UPGMA analysis separated P. atlantica subsp. mutica and cabulica from P. atlantica and P. eurycarpa. Subspecies mutica and cabulica were two closest genotypes; hence, P. atlantica subsp. mutica could be classified as a distinct species as P. mutica and the cabulica as a subspecies of P. mutica. This study revealed that P. eurycarpa is synonym for P. atlantica subsp. kurdica and should be considered distinct from P. atlantica; however, P. atlantica showed a closer genetic similarity to P. eurycarpa than the other species.