The clustered Hox genes show a conserved role in patterning the body axis of bilaterian metazoans. Increasingly, a broader phylogenetic sampling of non-model system organisms is being examined to detect a correlation, if any, between Hox gene evolution, and body plan innovations. To assess how Hox gene expression and function evolve with changing cluster arrangements, we must be able to reliably assign gene orthologies between Hox genes. Recent evidence suggests that a four-gene proto-Hox cluster duplicated to form the precursor of the present cluster and an additional sister-cluster, the ParaHox group. Here, phylogenetic methods are used to determine Hox-gene orthologies and to infer probable clustering events leading to the current bilaterian Hox complement. This analysis supports the ParaHox hypothesis and gives first confirmation that ind (intermediate neuroblasts defective) is an anterior ParaHox ortholog from protostomes. This analysis supports a proto-Hox cluster of four genes in which the central-class member of the ParaHox cluster may have been lost. It is also proposed here that ancestral diploblasts had central-class members of both Hox and ParaHox clusters. Primitive Hox gene ancestors are estimated by phylogenetic methods and found to have no strong affinity to any particular class of extant Hox members.