Phylogenetic relationship between species of Allium section Cepa and A. roylei (section Rhizirideum) have been inferred from nuclear DNA variation (RAPDs; nDNA dataset) and from morphological, pollen epidermis texture, chromosomal and chemical variation (supranuclear dataset). These sets were complemented with data, taken from the literature, on cpDNA variation and crossability. The trees produced with the supranuclear, nDNA and cpDNA datasets were compared by using the topology of the most parsimonious tree of one dataset as the constraint for the construction of a most parsimonious tree of another dataset. The accuracy of the trees were evaluated by calculating several Consistency and Incongruence Indices. The constrained tree of supranuclear-nDNA dataset showed the highest index values. The tree topologies of the supranuclear and cpDNA datasets were the least similar. The cpDNA tree and crossability dendrograms were identical. The most important difference between the nDNA-supranuclear trees and the cpDNA-crossability trees pertains to the position of Allium roylei, which is much closer to the clade A. cepal/A. vavilovii in the cpDNA tree than in the nDNA tree. This difference is considered to be the result of chloroplast capture from one species to another after an introgression event. A shorter distance between species inferred from a cpDNA tree from a nDNA or comparable tree might be indicative for the level of crossability.