All Metazoan nuclear genomes underwent a continuous process of both complete and partial genetic material gain and loss. The forces modulating these events are also subject to the strict interaction between nuclear and mitochondrial (mt) genome. In this context we investigate the evolution of nuclear genes encoding proteins which target the mitochondrion, with a particular attention to genes involved in oxidative phosphorylation (OXPHOS), one of the most ancient and conserved functions. To examine thoroughly the evolutionary strategies that preserve OXPHOS and coordinate the two cellular genomes, a comparative analysis has been carried out for 78 OXPHOS gene families in several Metazoa (insects, tunicates, fishes and mammals). We demonstrate that the duplication rate of OXPHOS genes increases passing from invertebrates to vertebrates consistently with the total increase in genome size, but all species are prone to negatively select OXPHOS duplicates compared to the general trend of nuclear gene families. These results are consistent with the ‘balance hypothesis’ and, at least in insects, the expression of duplicate genes is low and strongly testis-biased.