HOX gene complement and expression in the planarian Schmidtea mediterranea

Background: Freshwater planarians are well known for their regenerative abilities. Less well known is how planarians maintain spatial patterning in long-lived adult animals or how they re-pattern tissues during regeneration. HOX genes are good candidates to regulate planarian spatial patterning, yet the full complement or genomic clustering of planarian HOX genes has not yet been described, primarily because only a few have been detectable by in situ hybridization, and none have given morphological phenotypes when knocked down by RNAi. Results: Because the planarian Schmidtea mediterranea (S. med) is unsegmented, appendage-less, and morphologically simple, it has been proposed that it may have a simplified HOX gene complement. Here we argue against this hypothesis and show that S. med has a total of 13 HOX genes, which represent homologs to all major axial categories, and can be detected by whole-mount in situ hybridization using a highly-sensitive method. In addition, we show that planarian HOX genes do not cluster in the genome, yet 5/13 have retained aspects of axially-restricted expression. Finally, we confirm HOX gene axial expression by RNA-deep-sequencing 6 anterior-to-posterior “zones” of the animal, which we provide as a dataset to the community to discover other axially-restricted transcripts. Conclusions: Freshwater planarians have an unappreciated HOX gene complexity, with all major axial categories represented. However, we conclude based on adult expression patterns that planarians have a derived body plan and their asexual lifestyle may have allowed for large changes in HOX expression from the last common ancestor between arthropods, flatworms, and vertebrates. Using our in situ method and axial zone RNAseq data, it should be possible to further understand the pathways that pattern the anterior-posterior axis of adult planarians. 10 animals were immobilized on peltier coolers and carefully amputated into 6 Anterior-Posterior zones, as illustrated in Figure 2 of the manuscript. RNA was extracted from each zone, and deep sequenced to a depth >52 million reads.