Clustered brachiopod Hox genes are not expressed collinearly and are associated with lophotrochozoan novelties

Temporal collinearity is often regarded as the force preserving Hox clusters in vertebrate genomes. Studies that combine genomic and gene expression data in invertebrates would allow generalizing this observation across all animals, but are scarce, particularly within Lophotrochozoa (e.g., snails and segmented worms). Here, we use two brachiopod species -Terebratalia transversa, Novocrania anomala- to characterize the complement, cluster and expression of their Hox genes. T. transversa has an ordered, split cluster with ten genes (lab, pb, Hox3, dfd, scr, lox5, antp, lox4, post2, post1), while N. anomala has nine (missing post1). Our in situ hybridization, qPCR and stage specific transcripto