“…Since the discovery of the Forkhead domain in the fruit fly Drosophila melanogaster (Weigel et al 1989), Fox genes have been studied in a wide range of traditional developmental systems, mostly vertebrates (Mazet et al 2003; Hannenhalli & Kaestner 2009; Lee & Frasch 2004; Jackson et al 2010; Golson & Kaestner 2016). The initial description of 15 Fox gene classes in chordates, each identified by a letter (Kaestner et al 2000), and the establishment of a unified nomenclature facilitated phylogenetic analyses and comparisons with other major invertebrate clades, such as hemichordates (Fritzenwanker et al 2014), molluscs (Mei et al ., 2014; Wu et al ., 2020), platyhelminthes (Pascual-Carreras et al 2021), panarthropods (Schomburg et al 2022), cnidarians (Magie et al 2005), and sponges (Larroux et al 2006), as well as animal outgroups (Larroux et al 2008), ultimately uncovering a complex evolutionary history for this large family of transcription factors. Today, Fox genes are classified into 26 classes belonging to two major clades (Larroux et al 2008; Mazet et al 2003; Hannenhalli & Kaestner 2009; Benayoun et al 2011; Kaestner et al 2000; Pascual-Carreras et al 2021), where ancestral duplication events (e.g., the former class foxQ split into foxQ1 and foxQ2 , foxN into foxN1/4 and foxN2/3 , foxL into foxL1 and foxL2 , and foxJ into foxJ1 and foxJ2/3 ), gene innovations (e.g., foxR and foxS are unique of vertebrates, foxT is a novelty of panarthropods), expansions and losses (e.g., foxAB in vertebrates, foxQ2 in tetrapods, and foxAB , foxE , foxH , foxI in arthropods) are common (Paps et al 2012; Schomburg et al 2022; Mazet et al 2003; Wotton & Shimeld 2006; Tu et al 2006).…”