Whole genome duplication (WGD) generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. After WGD, retained ohnologues can undergo subfunctionalisation to partition ancestral functions and/or neofunctionalisation, where one copy assumes a new function. We previously found that there had been a WGD in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues were produced by WGD versus tandem duplications. Understanding this is key to determining the relative contributions of tandem duplications and WGD to arachnopulmonate genome evolution. Here we characterised the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that the majority of retained duplicate homeobox genes in spiders likely originated from WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE, in all eight species. Consistently, we noticed one degenerated copy of each cluster in terms of gene content and organisation while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalisation between the duplicated NK genes in the spider Parasteatoda tepidariorum to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.