Two different genomic regions (ND2 mitochondrial gene and EF1-alpha intron) were PCR amplified, cloned and sequenced for the ten known honey bee species collected within their natural range distribution. DNA sequences were analyzed using parsimony, distance and maximum likelihood methods to investigate phylogenetic relationships within Apis. The phylogenetic analyses strongly supported the basic topology recoverable from morphometric analysis, grouping the honey bees into three major clusters: giant bees (A. dorsata, A. binghami, and A. laboriosa), dwarf bees (A. andreniformis and A. florea), and cavity-nesting bees (A. mellifera, A. cerana, A. koschevnikovi, A. nuluensis, and A. nigrocincta). However, the clade of Asian cavity-nesting bees included paraphyletic taxa. Exemplars of Apis cerana collected from divergent portions of its range were less related to each other than were sympatric A. cerana, A. nuluensis, and A. nigrocincta taxa. Nucleotide sequence divergence between allopatrically distributed western (A. mellifera) and eastern (A. cerana, A. koschevnikovi, A. nigrocincta, and A. nuluensis) cavity-nesting species, around 18% for the mitochondrial gene and 10-15% for the nuclear intron, suggested an earlier divergence for these groups than previously estimated from morphometric and behavioral studies. This latter finding neccessitates reevaluation of the hypothesized origin of extant European, African, and west Asian Apis mellifera. Sequence divergence between A. laboriosa and A. dorsata was consistent with behavioral data and supports the species status of A. laboriosa.