Phylogenetic analysis of forkhead transcription factors in the Panarthropoda

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.