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

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation, however their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of eleven Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using HCR (hybridization chain reaction) probes for in-situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox genes expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

Genomic and transcriptomic survey of bryozoan Hox and ParaHox genes with emphasis on phylactolaemate bryozoans

BACKGROUND

Bryozoans are mostly sessile aquatic colonial invertebrates belonging to the clade Lophotrochozoa, which unites many protostome bilaterian phyla such as molluscs, annelids and brachiopods. While Hox and ParaHox genes have been extensively studied in various lophotrochozoan lineages, investigations on Hox and ParaHox gene complements in bryozoans are scarce.

RESULTS

Herein, we present the most comprehensive survey of Hox and ParaHox gene complements in bryozoans using four genomes and 35 transcriptomes representing all bryozoan clades: Cheilostomata, Ctenostomata, Cyclostomata and Phylactolaemata. Using similarity searches, phylogenetic analyses and detailed manual curation, we have identified five Hox genes in bryozoans (pb, Dfd, Lox5, Lox4 and Post2) and one ParaHox gene (Cdx). Interestingly, we observed lineage-specific duplication of certain Hox and ParaHox genes (Dfd, Lox5 and Cdx) in some bryozoan lineages.

CONCLUSIONS

The bryozoan Hox cluster does not retain the ancestral lophotrochozoan condition but appears relatively simple (includes only five genes) and broken into two genomic regions, characterized by the loss and duplication of serval genes. Importantly, bryozoans share the lack of two Hox genes (Post1 and Scr) with their proposed sister-taxon, Phoronida, which suggests that those genes were missing in the most common ancestor of bryozoans and phoronids.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

The gastropod Lottia peitaihoensis as a model to study the body patterning of trochophore larvae

The body patterning of trochophore larvae is important for understanding spiralian evolution and the origin of the bilateral body plan. However, considerable variations are observed among spiralian lineages, which have adopted varied strategies to develop trochophore larvae or even omit a trochophore stage. Some spiralians, such as patellogastropod mollusks, are suggested to exhibit ancestral traits by producing equal-cleaving fertilized eggs and possessing "typical" trochophore larvae. In recent years, we developed a potential model system using the patellogastropod Lottia peitaihoensis (= Lottia goshimai). Here, we introduce how the species were selected and establish sources and techniques, including gene knockdown, ectopic gene expression, and genome editing. Investigations on this species reveal essential aspects of trochophore body patterning, including organizer signaling, molecular and cellular processes connecting the various developmental functions of the organizer, the specification and behaviors of the endomesoderm and ectomesoderm, and the characteristic dorsoventral decoupling of Hox expression. These findings enrich the knowledge of trochophore body patterning and have important implications regarding the evolution of spiralians as well as bilateral body plans.

Shell field morphogenesis in the polyplacophoran mollusk Acanthochitona rubrolineata

BACKGROUND

The polyplacophoran mollusks (chitons) possess serially arranged shell plates. This feature is unique among mollusks and believed to be essential to explore the evolution of mollusks as well as their shells. Previous studies revealed several cell populations in the dorsal epithelium (shell field) of polyplacophoran larvae and their roles in the formation of shell plates. Nevertheless, they provide limited molecular information, and shell field morphogenesis remains largely uninvestigated.

RESULTS

In the present study, we investigated shell field development in the chiton Acanthochitona rubrolineata based on morphological characteristics and molecular patterns. A total of four types of tissue could be recognized from the shell field of A. rubrolineata. The shell field comprised not only the centrally located, alternatively arranged plate fields and ridges, but also the tissues surrounding them, which were the precursors of the girdle and we termed as the girdle field. The girdle field exhibited a concentric organization composed of two circularly arranged tissues, and spicules were only developed in the outer circle. Dynamic engrailed expression and F-actin (filamentous actin) distributions revealed relatively complicated morphogenesis of the shell field. The repeated units (plate fields and ridges) were gradually established in the shell field, seemingly different from the manners used in the segmentation of Drosophila or vertebrates. The seven repeated ridges also experienced different modes of ontogenesis from each other. In the girdle field, the presumptive spicule-formation cells exhibited different patterns of F-actin aggregations as they differentiate.

CONCLUSIONS

These results reveal the details concerning the structure of polyplacophoran shell field as well as its morphogenesis. They would contribute to exploring the mechanisms of polyplacophoran shell development and molluscan shell evolution.

Brachiopod and mollusc biomineralisation is a conserved process that was lost in the phoronid-bryozoan stem lineage

BACKGROUND

Brachiopods and molluscs are lophotrochozoans with hard external shells which are often believed to have evolved convergently. While palaeontological data indicate that both groups are descended from biomineralising Cambrian ancestors, the closest relatives of brachiopods, phoronids and bryozoans, are mineralised to a much lower extent and are comparatively poorly represented in the Palaeozoic fossil record. Although brachiopod and mollusc shells are structurally analogous, genomic and proteomic evidence indicates that their formation involves a complement of conserved, orthologous genes. Here, we study a set of genes comprised of 3 homeodomain transcription factors, one signalling molecule and 6 structural proteins which are implicated in mollusc and brachiopod shell formation, search for their orthologs in transcriptomes or genomes of brachiopods, phoronids and bryozoans, and present expression patterns of 8 of the genes in postmetamorphic juveniles of the rhynchonelliform brachiopod T. transversa.

RESULTS

Transcriptome and genome searches for the 10 target genes in the brachiopods Terebratalia transversa, Lingula anatina, Novocrania anomala, the bryozoans Bugula neritina and Membranipora membranacea, and the phoronids Phoronis australis and Phoronopsis harmeri resulted in the recovery of orthologs of the majority of the genes in all taxa. While the full complement of genes was present in all brachiopods with a single exception in L. anatina, a bloc of four genes could consistently not be retrieved from bryozoans and phoronids. The genes engrailed, distal-less, ferritin, perlucin, sp1 and sp2 were shown to be expressed in the biomineralising mantle margin of T. transversa juveniles.

CONCLUSIONS

The gene expression patterns we recovered indicate that while mineralised shells in brachiopods and molluscs are structurally analogous, their formation builds on a homologous process that involves a conserved complement of orthologous genes. Losses of some of the genes related to biomineralisation in bryozoans and phoronids indicate that loss of the capacity to form mineralised structures occurred already in the phoronid-bryozoan stem group and supports the idea that mineralised skeletons evolved secondarily in some of the bryozoan subclades.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

Various evaluations of the last common bilaterian ancestor (lcba) currently suggest that it resembled either a microscopic, non-segmented motile adult; or, on the contrary, a complex segmented adult motile urbilaterian. These fundamental inconsistencies remain largely unexplained. A majority of multidisciplinary data regarding sedentary adult ancestral bilaterian organization is overlooked. The sedentary-pelagic model is supported now by a number of novel developmental, paleontological and molecular phylogenetic data: (1) data in support of sedentary sponges, in the adult stage, as sister to all other Metazoa; (2) a similarity of molecular developmental pathways in both adults and larvae across sedentary sponges, cnidarians, and bilaterians; (3) a cnidarian-bilaterian relationship, including a unique sharing of a bona fide Hox-gene cluster, of which the evolutionary appearance does not connect directly to a bilaterian motile organization; (4) the presence of sedentary and tube-dwelling representatives of the main bilaterian clades in the early Cambrian; (5) an absence of definite taxonomic attribution of Ediacaran taxa reconstructed as motile to any true bilaterian phyla; (6) a similarity of tube morphology (and the clear presence of a protoconch-like apical structure of the Ediacaran sedentary Cloudinidae) among shells of the early Cambrian, and later true bilaterians, such as semi-sedentary hyoliths and motile molluscs; (7) recent data that provide growing evidence for a complex urbilaterian, despite a continuous molecular phylogenetic controversy. The present review compares the main existing models and reconciles the sedentary model of an urbilaterian and the model of a larva-like lcba with a unified sedentary(adult)-pelagic(larva) model of the lcba.

Uniaxial orientation of β-chitin nanofibres used as an organic framework in the scales of a hot vent snail

Organisms use various forms and orientations of chitin nanofibres to make structures with a wide range of functions, from insect wings to mussel shells. Lophotrochozoan animals such as snails and annelid worms possess an ancient ‘biomineralization toolkit’, enabling them to flexibly and rapidly evolve unique hard parts. The scaly-foot snail is a gastropod endemic to deep-sea hydrothermal vents, unique in producing dermal sclerites used as sites of sulfur detoxification. Once considered to be strictly proteinaceous, recent data pointed to the presence of chitin in these sclerites, but direct evidence is still lacking. Here, we show that β-chitin fibres (approx. 5% of native weight) are indeed the building framework, through a combination of solid-state nuclear magnetic resonance spectroscopy, wide-angle X-ray diffraction, and electron microscopy. The fibres are uniaxially oriented, likely forming a structural basis for column-like channels into which the scaly-foot snail is known to actively secrete sulfur waste—expanding the known function of chitinous hard parts in animals. Our results add to the existing evidence that animals are capable of modifying and co-opting chitin synthesis pathways flexibly and rapidly, in order to serve novel functions during their evolution.

The genome of the poecilogonous annelid Streblospio benedicti

Streblospio benedicti is a common marine annelid that has become an important model for developmental evolution. It is the only known example of poecilogony (where two distinct developmental modes occur within a single species) that is due to a heritable difference in egg size. The dimorphic developmental programs and life-histories exhibited in this species depend on differences within the genome, making it an optimal model for understanding the genomic basis of developmental divergence. Studies using S. benedicti have begun to uncover the genetic and genomic principles that underlie developmental uncoupling, but until now they have been limited by the lack of availability of genomic tools. Here we present an annotated chromosomal-level genome assembly of S. benedicti generated from a combination of Illumina reads, Nanopore long reads, Chicago and Hi-C chromatin interaction sequencing, and a genetic map from experimental crosses. At 701.4 Mb, the S. benedicti genome is the largest annelid genome to date that has been assembled to chromosomal scaffolds. The complete genome of S. benedicti is valuable for functional genomic analyses of development and evolution, as well as phylogenetic comparison within the annelida and the Lophotrochozoa. Despite having two developmental modes, there is no evidence of genome duplication or substantial gene number expansions. Instead, lineage specific repeats account for much of the expansion of this genome compared to other annelids.

Co-option of the limb patterning program in cephalopod eye development

Background

Across the Metazoa, similar genetic programs are found in the development of analogous, independently evolved, morphological features. The functional significance of this reuse and the underlying mechanisms of co-option remain unclear. Cephalopods have evolved a highly acute visual system with a cup-shaped retina and a novel refractive lens in the anterior, important for a number of sophisticated behaviors including predation, mating, and camouflage. Almost nothing is known about the molecular-genetics of lens development in the cephalopod.

Results

Here we identify the co-option of the canonical bilaterian limb patterning program during cephalopod lens development, a functionally unrelated structure. We show radial expression of transcription factors SP6-9/sp1, Dlx/dll, Pbx/exd, Meis/hth, and a Prdl homolog in the squid Doryteuthis pealeii, similar to expression required in Drosophila limb development. We assess the role of Wnt signaling in the cephalopod lens, a positive regulator in the developing Drosophila limb, and find the regulatory relationship reversed, with ectopic Wnt signaling leading to lens loss.

Conclusion

This regulatory divergence suggests that duplication of SP6-9 in cephalopods may mediate the co-option of the limb patterning program. Thus, our study suggests that this program could perform a more universal developmental function in radial patterning and highlights how canonical genetic programs are repurposed in novel structures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01182-2.

The evolution of the metazoan Toll receptor family and its expression during protostome development

BACKGROUND

Toll-like receptors (TLRs) play a crucial role in immunity and development. They contain leucine-rich repeat domains, one transmembrane domain, and one Toll/IL-1 receptor domain. TLRs have been classified into V-type/scc and P-type/mcc TLRs, based on differences in the leucine-rich repeat domain region. Although TLRs are widespread in animals, detailed phylogenetic studies of this gene family are lacking. Here we aim to uncover TLR evolution by conducting a survey and a phylogenetic analysis in species across Bilateria. To discriminate between their role in development and immunity we furthermore analyzed stage-specific transcriptomes of the ecdysozoans Priapulus caudatus and Hypsibius exemplaris, and the spiralians Crassostrea gigas and Terebratalia transversa.

RESULTS

We detected a low number of TLRs in ecdysozoan species, and multiple independent radiations within the Spiralia. V-type/scc and P-type/mcc type-receptors are present in cnidarians, protostomes and deuterostomes, and therefore they emerged early in TLR evolution, followed by a loss in xenacoelomorphs. Our phylogenetic analysis shows that TLRs cluster into three major clades: clade α is present in cnidarians, ecdysozoans, and spiralians; clade β in deuterostomes, ecdysozoans, and spiralians; and clade γ is only found in spiralians. Our stage-specific transcriptome and in situ hybridization analyses show that TLRs are expressed during development in all species analyzed, which indicates a broad role of TLRs during animal development.

CONCLUSIONS

Our findings suggest that a clade α TLR gene (TLR-Ca) and a clade β/γ TLR gene (TLR-Cβ/γ) were already present in the cnidarian-bilaterian common ancestor. However, although TLR-Ca was conserved in cnidarians, TLR-Cβ/γ was lost during the early evolution of these taxa. Moreover, TLR-Cβ/γ duplicated to generate TLR-Cβ and TLR-Cγ in the lineage to the last common protostome-deuterostome ancestor. TLR-Ca, TLR-Cβ and TLR-Cγ further expanded generating the three major TLR clades. While all three clades radiated in several spiralian lineages, specific TLRs clades have been presumably lost in other lineages. Furthermore, the expression of the majority of these genes during protostome ontogeny suggests a likely role in development.

HES and Mox genes are expressed during early mesoderm formation in a mollusk with putative ancestral features

The mesoderm is considered the youngest of the three germ layers. Although its morphogenesis has been studied in some metazoans, the molecular components underlying this process remain obscure for numerous phyla including the highly diverse Mollusca. Here, expression of Hairy and enhancer of split (HES), Mox, and myosin heavy chain (MHC) was investigated in Acanthochitona fascicularis, a representative of Polyplacophora with putative ancestral molluscan features. While AfaMHC is expressed throughout myogenesis, AfaMox1 is only expressed during early stages of mesodermal band formation and in the ventrolateral muscle, an autapomorphy of the polyplacophoran trochophore. Comparing our findings to previously published data across Metazoa reveals Mox expression in the mesoderm in numerous bilaterians including gastropods, polychaetes, and brachiopods. It is also involved in myogenesis in molluscs, annelids, tunicates, and craniates, suggesting a dual role of Mox in mesoderm and muscle formation in the last common bilaterian ancestor. AfaHESC2 is expressed in the ectoderm of the polyplacophoran gastrula and later in the mesodermal bands and in putative neural tissue, whereas AfaHESC7 is expressed in the trochoblasts of the gastrula and during foregut formation. This confirms the high developmental variability of HES gene expression and demonstrates that Mox and HES genes are pleiotropic.

Molecular and morphological analysis of the developing nemertean brain indicates convergent evolution of complex brains in Spiralia

BACKGROUND

The brain anatomy in the clade Spiralia can vary from simple, commissural brains (e.g., gastrotrichs, rotifers) to rather complex, partitioned structures (e.g., in cephalopods and annelids). How often and in which lineages complex brains evolved still remains unclear. Nemerteans are a clade of worm-like spiralians, which possess a complex central nervous system (CNS) with a prominent brain, and elaborated chemosensory and neuroglandular cerebral organs, which have been previously suggested as homologs to the annelid mushroom bodies. To understand the developmental and evolutionary origins of the complex brain in nemerteans and spiralians in general, we investigated details of the neuroanatomy and gene expression in the brain and cerebral organs of the juveniles of nemertean Lineus ruber.

RESULTS

In the juveniles, the CNS is already composed of all major elements present in the adults, including the brain, paired longitudinal lateral nerve cords, and an unpaired dorsal nerve cord, which suggests that further neural development is mostly related with increase in the size but not in complexity. The ultrastructure of the juvenile cerebral organ revealed that it is composed of several distinct cell types present also in the adults. The 12 transcription factors commonly used as brain cell type markers in bilaterians show region-specific expression in the nemertean brain and divide the entire organ into several molecularly distinct areas, partially overlapping with the morphological compartments. Additionally, several of the mushroom body-specific genes are expressed in the developing cerebral organs.

CONCLUSIONS

The dissimilar expression of molecular brain markers between L. ruber and the annelid Platynereis dumerilii indicates that the complex brains present in those two species evolved convergently by independent expansions of non-homologous regions of a simpler brain present in their last common ancestor. Although the same genes are expressed in mushroom bodies and cerebral organs, their spatial expression within organs shows apparent differences between annelids and nemerteans, indicating convergent recruitment of the same genes into patterning of non-homologous organs or hint toward a more complicated evolutionary process, in which conserved and novel cell types contribute to the non-homologous structures.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

Comparative and evolutionary analyses reveal conservation and divergence of the notch pathway in lophotrochozoa

Lophotrochozoan species exhibit wide morphological diversity; however, the molecular basis underlying this diversity remains unclear. Here, we explored the evolution of Notch pathway genes across 37 metazoan species via phylogenetic and molecular evolutionary studies with emphasis on the lophotrochozoans. We displayed the components of Notch pathway in metazoans and found that Delta and Hes/Hey-related genes, as well as their functional domains, are duplicated in lophotrochozoans. Comparative transcriptomics analyses allow us to pinpoint sequence divergence of multigene families in the Notch signalling pathway. We identified the duplication mechanism of a mollusc-specific gene, Delta2, and found it displayed complementary expression throughout development. Furthermore, we found the functional diversification not only in expanded genes in the Notch pathway (Delta and Hes/Hey-related genes), but also in evolutionary conservative genes (Notch, Presenilin, and Su(H)). Together, this comprehensive study demonstrates conservation and divergence within the Notch pathway, reveals evolutionary relationships among metazoans, and provides evidence for the occurrence of developmental diversity in lophotrochozoans, as well as a basis for future gene function studies.

FGF signaling acts on different levels of mesoderm development within Spiralia

FGF signaling is involved in mesoderm induction in members of deuterostomes (e.g. tunicates, hemichordates), but not in flies and nematodes, in which it has a role in mesoderm patterning and migration. However, we need comparable studies in other protostome taxa in order to decipher whether this mesoderm-inducing function of FGF extends beyond the lineage of deuterostomes. Here, we investigated the role of FGF signaling in mesoderm development in three species of lophophorates, a clade within the protostome group Spiralia. Our gene expression analyses show that the mesodermal molecular patterning is conserved between brachiopods and phoronids, but the spatial and temporal recruitment of transcription factors differs significantly. Moreover, the use of the inhibitor SU5402 demonstrates that FGF signaling is involved in different steps of mesoderm development, as well as in morphogenetic movements of gastrulation and axial elongation. Our findings suggest that the mesoderm-inducing role of FGF extends beyond the group of deuterostomes.

The Iron-Responsive Genome of the Chiton Acanthopleura granulata

Molluscs biomineralize structures that vary in composition, form, and function, prompting questions about the genetic mechanisms responsible for their production and the evolution of these mechanisms. Chitons (Mollusca, Polyplacophora) are a promising system for studies of biomineralization because they build a range of calcified structures including shell plates and spine- or scale-like sclerites. Chitons also harden the calcified teeth of their rasp-like radula with a coat of iron (as magnetite). Here we present the genome of the West Indian fuzzy chiton Acanthopleura granulata, the first from any aculiferan mollusc. The A. granulata genome contains homologs of many genes associated with biomineralization in conchiferan molluscs. We expected chitons to lack genes previously identified from pathways conchiferans use to make biominerals like calcite and nacre because chitons do not use these materials in their shells. Surprisingly, the A. granulata genome has homologs of many of these genes, suggesting that the ancestral mollusc may have had a more diverse biomineralization toolkit than expected. The A. granulata genome has features that may be specialized for iron biomineralization, including a higher proportion of genes regulated directly by iron than other molluscs. A. granulata also produces two isoforms of soma-like ferritin: one is regulated by iron and similar in sequence to the soma-like ferritins of other molluscs, and the other is constitutively translated and is not found in other molluscs. The A. granulata genome is a resource for future studies of molluscan evolution and biomineralization.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

The Scaly-foot Snail genome and implications for the origins of biomineralised armour

The Scaly-foot Snail, Chrysomallon squamiferum, presents a combination of biomineralised features, reminiscent of enigmatic early fossil taxa with complex shells and sclerites such as sachtids, but in a recently-diverged living species which even has iron-infused hard parts. Thus the Scaly-foot Snail is an ideal model to study the genomic mechanisms underlying the evolutionary diversification of biomineralised armour. Here, we present a high-quality whole-genome assembly and tissue-specific transcriptomic data, and show that scale and shell formation in the Scaly-foot Snail employ independent subsets of 25 highly-expressed transcription factors. Comparisons with other lophotrochozoan genomes imply that this biomineralisation toolkit is ancient, though expression patterns differ across major lineages. We suggest that the ability of lophotrochozoan lineages to generate a wide range of hard parts, exemplified by the remarkable morphological disparity in Mollusca, draws on a capacity for dynamic modification of the expression and positioning of toolkit elements across the genome.

A technical review and guide to RNA fluorescence in situ hybridization

RNA-fluorescence in situ hybridization (FISH) is a powerful tool to visualize target messenger RNA transcripts in cultured cells, tissue sections or whole-mount preparations. As the technique has been developed over time, an ever-increasing number of divergent protocols have been published. There is now a broad selection of options available to facilitate proper tissue preparation, hybridization, and post-hybridization background removal to achieve optimal results. Here we review the technical aspects of RNA-FISH, examining the most common methods associated with different sample types including cytological preparations and whole-mounts. We discuss the application of commonly used reagents for tissue preparation, hybridization, and post-hybridization washing and provide explanations of the functional roles for each reagent. We also discuss the available probe types and necessary controls to accurately visualize gene expression. Finally, we review the most recent advances in FISH technology that facilitate both highly multiplexed experiments and signal amplification for individual targets. Taken together, this information will guide the methods development process for investigators that seek to perform FISH in organisms that lack documented or optimized protocols.

Hox gene expression during development of the phoronid Phoronopsis harmeri

Background

Phoronida is a small group of marine worm-like suspension feeders, which together with brachiopods and bryozoans form the clade Lophophorata. Although their development is well studied on the morphological level, data regarding gene expression during this process are scarce and restricted to the analysis of relatively few transcription factors. Here, we present a description of the expression patterns of Hox genes during the embryonic and larval development of the phoronid Phoronopsis harmeri.

Results

We identified sequences of eight Hox genes in the transcriptome of Ph. harmeri and determined their expression pattern during embryonic and larval development using whole mount in situ hybridization. We found that none of the Hox genes is expressed during embryonic development. Instead their expression is initiated in the later developmental stages, when the larval body is already formed. In the investigated initial larval stages the Hox genes are expressed in the non-collinear manner in the posterior body of the larvae: in the telotroch and the structures that represent rudiments of the adult worm. Additionally, we found that certain head-specific transcription factors are expressed in the oral hood, apical organ, preoral coelom, digestive system and developing larval tentacles, anterior to the Hox-expressing territories.

Conclusions

The lack of Hox gene expression during early development of Ph. harmeri indicates that the larval body develops without positional information from the Hox patterning system. Such phenomenon might be a consequence of the evolutionary intercalation of the larval form into an ancestral life cycle of phoronids. The observed Hox gene expression can also be a consequence of the actinotrocha representing a “head larva”, which is composed of the most anterior body region that is devoid of Hox gene expression. Such interpretation is further supported by the expression of head-specific transcription factors. This implies that the Hox patterning system is used for the positional information of the trunk rudiments and is, therefore, delayed to the later larval stages. We propose that a new body form was intercalated to the phoronid life cycle by precocious development of the anterior structures or by delayed development of the trunk rudiment in the ancestral phoronid larva.

Unravelling spiral cleavage

Snails, earthworms and flatworms are remarkably different animals, but they all exhibit a very similar mode of early embryogenesis: spiral cleavage. This is one of the most widespread developmental programs in animals, probably ancestral to almost half of the animal phyla, and therefore its study is essential for understanding animal development and evolution. However, our knowledge of spiral cleavage is still in its infancy. Recent technical and conceptual advances, such as the establishment of genome editing and improved phylogenetic resolution, are paving the way for a fresher and deeper look into this fascinating early cleavage mode.

Dorsoventral decoupling of Hox gene expression underpins the diversification of molluscs

In contrast to the Hox genes in arthropods and vertebrates, those in molluscs show diverse expression patterns with differences reported among lineages. Here, we investigate 2 phylogenetically distant molluscs, a gastropod and a polyplacophoran, and show that the Hox expression in both species can be divided into 2 categories. The Hox expression in the ventral ectoderm generally shows a canonical staggered pattern comparable to the patterns of other bilaterians and likely contributes to ventral patterning, such as neurogenesis. The other category of Hox expression on the dorsal side is strongly correlated with shell formation and exhibits lineage-specific characteristics in each class of mollusc. This generalized model of decoupled dorsoventral Hox expression is compatible with known Hox expression data from other molluscan lineages and may represent a key characteristic of molluscan Hox expression. These results support the concept of widespread staggered Hox expression in Mollusca and reveal aspects that may be related to the evolutionary diversification of molluscs. We propose that dorsoventral decoupling of Hox expression allowed lineage-specific dorsal and ventral patterning, which may have facilitated the evolution of diverse body plans in different molluscan lineages.

A developmental perspective on the evolution of the nervous system

The evolution of nervous systems in animals has always fascinated biologists, and thus multiple evolutionary scenarios have been proposed to explain the appearance of neurons and complex neuronal centers. However, the absence of a robust phylogenetic framework for animal interrelationships, the lack of a mechanistic understanding of development, and a recapitulative view of animal ontogeny have traditionally limited these scenarios. Only recently, the integration of advanced molecular and morphological studies in a broad range of animals has allowed to trace the evolution of developmental and neuronal characters on a better-resolved animal phylogeny. This has falsified most traditional scenarios for nervous system evolution, paving the way for the emergence of new testable hypotheses. Here we summarize recent progress in studies of nervous system development in major animal lineages and formulate some of the arising questions. In particular, we focus on how lineage analyses of nervous system development and a comparative study of the expression of neural-related genes has influenced our understanding of the evolution of an elaborated central nervous system in Bilateria. We argue that a phylogeny-guided study of neural development combining thorough descriptive and functional analyses is key to establish more robust scenarios for the origin and evolution of animal nervous systems.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

Sclerite-bearing annelids from the lower Cambrian of South China

Cambrian annelids are strikingly diverse and reveal important details of annelid character acquisition. Their contribution, however, to a wider understanding of the evolution of the trochozoans (encompassing the annelids as well as such groups as the brachiopods and molluscs) remains limited. Thus the early annelids had been linked to a variety of cataphract Cambrian metazoans, notably Wiwaxia and the halkieriids, but recent work assigns such fossils to stem-group molluscs. Here we report two new annelids from the Lower Cambrian Chengjiang Lagerstätte, South China. Ipoliknus avitus n. gen., n. sp. is biramous with neurochaetae and notochaetae, but significantly also bears dorsal spinose sclerites and dorso-lateral dentate sclerites. Adelochaeta sinensis n. gen., n. sp. is unique amongst Cambrian polychaetes in possessing the rod-like supports of the parapodia known as aciculae. This supports phylogenetic placement of Adelochaeta as sister to some more derived aciculate Palaeozoic taxa, but in contrast Ipoliknus is recovered as the most basal of the stem-group annelids. Sclerites and chaetae of I. avitus are interpreted respectively as the remnants and derivatives of a once more extensive cataphract covering that was a characteristic of more primitive trochozoans. The two sets of chaetae (noto- and neurochaetae) and two sets of sclerites (spinose and dentate) suggest that in a pre-annelid an earlier and more complete scleritome may have consisted of four zones of sclerites. Other cataphract taxa from the Lower Palaeozoic show a variety of scleritome configurations but establishing direct links with such basal annelids as Ipoliknus at present must remain conjectural.

Insights into the Evolution of Shells and Love Darts of Land Snails Revealed from Their Matrix Proteins

Over the past decade, many skeletal matrix proteins that are possibly related to calcification have been reported in various calcifying animals. Molluscs are among the most diverse calcifying animals and some gastropods have adapted to terrestrial ecological niches. Although many shell matrix proteins (SMPs) have already been reported in molluscs, most reports have focused on marine molluscs, and the SMPs of terrestrial snails remain unclear. In addition, some terrestrial stylommatophoran snails have evolved an additional unique calcified character, called a “love dart,” used for mating behavior. We identified 54 SMPs in the terrestrial snail Euhadra quaesita, and found that they contain specific domains that are widely conserved in molluscan SMPs. However, our results also suggest that some of them possibly have evolved independently by domain shuffling, domain recruitment, or gene co-option. We then identified four dart matrix proteins, and found that two of them are the same proteins as those identified as SMPs. Our results suggest that some dart matrix proteins possibly have evolved by independent gene co-option from SMPs during dart evolution events. These results provide a new perspective on the evolution of SMPs and “love darts” in land snails.

Hox gene expression in postmetamorphic juveniles of the brachiopod Terebratalia transversa

BACKGROUND

Hox genes encode a family of homeodomain containing transcription factors that are clustered together on chromosomes of many Bilateria. Some bilaterian lineages express these genes during embryogenesis in spatial and/or temporal order according to their arrangement in the cluster, a phenomenon referred to as collinearity. Expression of Hox genes is well studied during embryonic and larval development of numerous species; however, relatively few studies focus on the comparison of pre- and postmetamorphic expression of Hox genes in animals with biphasic life cycle. Recently, the expression of Hox genes was described for embryos and larvae of Terebratalia transversa, a rhynchonelliformean brachiopod, which possesses distinct metamorphosis from planktonic larvae to sessile juveniles. During premetamorphic development, T. transversa does not exhibit spatial collinearity and several of its Hox genes are recruited for the morphogenesis of novel structures. In our study, we determined the expression of Hox genes in postmetamorphic juveniles of T. transversa in order to examine metamorphosis-related changes of expression patterns and to test whether Hox genes are expressed in the spatially collinear way in the postmetamorphic juveniles.

RESULTS

Hox genes are expressed in a spatially non-collinear manner in juveniles, generally showing similar patterns as ones observed in competent larvae: genes labial and post1 are expressed in chaetae-related structures, sex combs reduced in the shell-forming epithelium, whereas lox5 and lox4 in dorso-posterior epidermis. After metamorphosis, expression of genes proboscipedia, hox3, deformed and antennapedia becomes restricted to, respectively, shell musculature, prospective hinge rudiments and pedicle musculature and epidermis.

CONCLUSIONS

All developmental stages of T. transversa, including postmetamorphic juveniles, exhibit a spatial non-collinear Hox genes expression with only minor changes observed between pre- and postmetamorphic stages. Our results are concordant with morphological observation that metamorphosis in rhynchonelliformean brachiopods, despite being rapid, is rather gradual. The most drastic changes in Hox gene expression patterns observed during metamorphosis could be explained by the inversion of the mantle lobe, which relocates some of the more posterior larval structures into the anterior edge of the juveniles. Co-option of Hox genes for the morphogenesis of novel structures is even more pronounced in postmetamorphic brachiopods when compared to larvae.

Staggered Hox expression is more widespread among molluscs than previously appreciated

Hox genes are expressed along the anterior-posterior body axis in a colinear fashion in the majority of bilaterians. Contrary to polyplacophorans, a group of aculiferan molluscs with conserved ancestral molluscan features, gastropods and cephalopods deviate from this pattern by expressing Hox genes in distinct morphological structures and not in a staggered fashion. Among conchiferans, scaphopods exhibit many similarities with gastropods, cephalopods and bivalves, however, the molecular developmental underpinnings of these similar traits remain unknown. We investigated Hox gene expression in developmental stages of the scaphopod Antalis entalis to elucidate whether these genes are involved in patterning morphological traits shared by their kin conchiferans. Scaphopod Hox genes are predominantly expressed in the foot and mantle but also in the central nervous system. Surprisingly, the scaphopod mid-stage trochophore exhibits a near-to staggered expression of all nine Hox genes identified. Temporal colinearity was not found and early-stage and late-stage trochophores, as well as postmetamorphic individuals, do not show any apparent traces of staggered expression. In these stages, Hox genes are expressed in distinct morphological structures such as the cerebral and pedal ganglia and in the shell field of early-stage trochophores. Interestingly, a re-evaluation of previously published data on early-stage cephalopod embryos and of the gastropod pre-torsional veliger shows that these developmental stages exhibit traces of staggered Hox expression. Considering our results and all gene expression and genomic data available for molluscs as well as other bilaterians, we suggest a last common molluscan ancestor with colinear Hox expression in predominantly ectodermal tissues along the anterior-posterior axis. Subsequently, certain Hox genes have been co-opted into the patterning process of distinct structures (apical organ or prototroch) in conchiferans.

Downregulation of HOXA3 in lung adenocarcinoma and its relevant molecular mechanism analysed by RT-qPCR, TCGA and in silico analysis

Recent studies have indicated that homeobox A3 (HOXA3) functions as a carcinogen in colon cancer and the methylation level of HOXA3 is significantly increased in lung adenocarcinoma (LUAD) tissues. However, at least to the best of our knowledge, few studies to date have been performed on HOXA3 in non-small cell lung cancer (NSCLC). Therefore, further studies on HOXA3 expression in NSCLC and the potential regulatory mechanisms are urgently required. In this study, HOXA3 expression in 55 tissues of cases of NSCLC and corresponding non-lung cancer tissues was detected by reverse transcription-quantitative PCR (RT-qPCR). In addition, the clinical significance of HOXA3 expression in NSCLC was evaluated using the Cancer Genome Atlas (TCGA) database. Bioinformatics analysis was then performed to elucidate the potential molecular mechanisms of action of HOXA3. Furthermore, the potential target microRNAs (miRNAs or miRs) of HOXA3 were predicted using miRWalk2.0. Based on Gene Expression Omnibus (GEO) and TGCA databases, standardized mean difference (SMD) and sROC methods were used for meta-analyses of the expression of potential target miRNAs of HOXA3 in NSCLC to evaluate their association with HOXA3. The results revealed that the HOXA3 expression levels in NSCLC, LUAD and lung squamous cell carcinoma (LUSC) were 0.1130±0.1398, 0.1295±0.16890 and 0.0906±0.0846, respectively. These values were all decreased compared with the normal tissues (0.1877±0.1975, 0.2337±0.2405 and 0.1249±0.0873, respectively, P<0.05). The TCGA database also revealed the low expression trend of HOXA3. The downregulation of HOXA3 may play an important role in the progression and the poor prognosis of LUAD. The TCGA database also suggested that HOXA3 in LUAD and LUSC tissues exhibited certain mutational levels. In addition, the methylation levels in the NSCLC, LUAD and LUSC tissues significantly increased [NSCLC: fold change (FC), 1.3226; P<0.001; LUAD: FC, 1.2712; P<0.001; and LUSC: FC, 1.3786; P<0.001]. According to the analyses using the Kyoto Encyclopedia of Genes and Genomes (KEGG), we found that the co-expression HOXA3 genes were mainly associated with the focal adhesion signalling pathway and the ECM-receptor interaction signalling pathway. Furthermore, the predicted miRNA, miR-372-3p, exhibited a high expression in both the NSCLC and LUAD tissues (P<0.05). On the whole, the findings of this study indicate that low HOXA3 expression may play a certain role in LUAD; however, its association with LUSC still requires further investigation. HOXA3 function may be achieved through different pathways or target miRNAs. However, the specific underlying mechanisms need to be confirmed through various functional studies.

Lophotrochozoan Zic Genes

Lophotrochozoa is a sister taxon of Ecdysozoa in the Protostomia that includes mollusks, annelids, brachiopods, and platyhelminths. Recent studies have clarified the structure, expression, and roles of lophotrochozoan Zic family genes. Zic genes in oligochaete annelid Tubifex tubifex (freshwater sludge worm) and polychaete annelid Capitella teleta (bristle worm) are commonly expressed in a subset of developing brain and mesoderm derivatives. The latter includes the naïve mesoderm and the associated chaetal sacs in each body segment, although the segmentation processes differ between the two species. Furthermore, in brachiopod Terebratalia transversa (lamp shell), Zic is expressed in the anterior ectodermal domains and mesodermal derivatives, including those associated with the chaetal sacs. This result suggests the common involvement of Zic genes in the development of chaetae, a lophotrochozoan novelty acquired in the course of evolution. In addition, the highly simplified lophotrochozoan Dicyema acuticephalum (dicyemid mesozoan, a cephalopod endoparasite), which lost its gut, nervous system, and muscles during evolution, expresses its Zic genes in hermaphroditic gonads, highlighting the role of Zic genes in germ cell development. The role of Zic in head regeneration was revealed in studies on platyhelminth Schmidtea mediterranea (freshwater planarian). Planarian Zic expression was induced in a subpopulation of neoblasts that includes adult pluripotent stem cells. It is needed for head regeneration and production of an anterior signaling center. Suppression of Wnt-β-catenin signaling underlies Zic-mediated head regeneration, reminiscent of Wnt-β-catenin suppression by vertebrate Zic genes. Taken together, studies on the lophotrochozoan Zic genes are essential to understanding not only the roles of these genes in body plan evolution but also the molecular mechanism underlying adult stem cell regulation.

An ancient FMRFamide-related peptide-receptor pair induces defence behaviour in a brachiopod larva

Animal behaviour often comprises spatially separated sub-reactions and even ciliated larvae are able to coordinate sub-reactions of complex behaviours (metamorphosis, feeding). How these sub-reactions are coordinated is currently not well understood. Neuropeptides are potential candidates for triggering larval behaviour. However, although their immunoreactivity has been widely analysed, their function in trochozoan larvae has only been studied for a few cases. Here, we investigate the role of neuropeptides in the defence behaviour of brachiopod larvae. When mechanically disturbed, the planktonic larvae of Terebratalia transversa protrude their stiff chaetae and sink down slowly. We identified endogenous FLRFamide-type neuropeptides (AFLRFamide and DFLRFamide) in T. transversa larvae and show that the protrusion of the chaetae as well as the sinking reaction can both be induced by each of these peptides. This also correlates with the presence of FLRFamidergic neurons in the apical lobe and adjacent to the trunk musculature. We deorphanized the AFLRFamide/DFLRFamide receptor and detected its expression in the same tissues. Furthermore, the ability of native and modified FLRFamide-type peptides to activate this receptor was found to correspond with their ability to trigger behavioural responses. Our results show how FLRFamide-type neuropeptides can induce two coherent sub-reactions in a larva with a simple nervous system.

Hox and Wnt pattern the primary body axis of an anthozoan cnidarian before gastrulation

Hox gene transcription factors are important regulators of positional identity along the anterior–posterior axis in bilaterian animals. Cnidarians (e.g., sea anemones, corals, and hydroids) are the sister group to the Bilateria and possess genes related to both anterior and central/posterior class Hox genes. Here we report a previously unrecognized domain of Hox expression in the starlet sea anemone, Nematostella vectensis, beginning at early blastula stages. We explore the relationship of two opposing Hox genes (NvAx6/NvAx1) expressed on each side of the blastula during early development. Functional perturbation reveals that NvAx6 and NvAx1 not only regulate their respective expression domains, but also interact with Wnt genes to pattern the entire oral–aboral axis. These findings suggest an ancient link between Hox/Wnt patterning during axis formation and indicate that oral–aboral domains are likely established during blastula formation in anthozoan cnidarians.

Hox genes regulate anterior–posterior axis formation but their role in cnidarians is unclear. Here, the authors disrupt Hox genes NvAx1 and NvAx6 in the starlet sea anemone, Nematostella vectensis, showing antagonist function in patterning the oral–aboral axis and a link to Wnt signaling.

Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates

Hox genes exert fundamental roles for proper regional specification along the main rostro-caudal axis of animal embryos. They are generally expressed in restricted spatial domains according to their position in the cluster (spatial colinearity)-a feature that is conserved across bilaterians. In jawed vertebrates (gnathostomes), the position in the cluster also determines the onset of expression of Hox genes (a feature known as whole-cluster temporal colinearity (WTC)), while in invertebrates this phenomenon is displayed as a subcluster-level temporal colinearity. However, little is known about the expression profile of Hox genes in jawless vertebrates (cyclostomes); therefore, the evolutionary origin of WTC, as seen in gnathostomes, remains a mystery. Here, we show that Hox genes in cyclostomes are expressed according to WTC during development. We investigated the Hox repertoire and Hox gene expression profiles in three different species-a hagfish, a lamprey and a shark-encompassing the two major groups of vertebrates, and found that these are expressed following a whole-cluster, temporally staggered pattern, indicating that WTC has been conserved during the past 500 million years despite drastically different genome evolution and morphological outputs between jawless and jawed vertebrates.

Nemertean and phoronid genomes reveal lophotrochozoan evolution and the origin of bilaterian heads

Nemerteans (ribbon worms) and phoronids (horseshoe worms) are closely related lophotrochozoans-a group of animals including leeches, snails and other invertebrates. Lophotrochozoans represent a superphylum that is crucial to our understanding of bilaterian evolution. However, given the inconsistency of molecular and morphological data for these groups, their origins have been unclear. Here, we present draft genomes of the nemertean Notospermus geniculatus and the phoronid Phoronis australis, together with transcriptomes along the adult bodies. Our genome-based phylogenetic analyses place Nemertea sister to the group containing Phoronida and Brachiopoda. We show that lophotrochozoans share many gene families with deuterostomes, suggesting that these two groups retain a core bilaterian gene repertoire that ecdysozoans (for example, flies and nematodes) and platyzoans (for example, flatworms and rotifers) do not. Comparative transcriptomics demonstrates that lophophores of phoronids and brachiopods are similar not only morphologically, but also at the molecular level. Despite dissimilar head structures, lophophores express vertebrate head and neuronal marker genes. This finding suggests a common origin of bilaterian head patterning, although different heads evolved independently in each lineage. Furthermore, we observe lineage-specific expansions of innate immunity and toxin-related genes. Together, our study reveals a dual nature of lophotrochozoans, where conserved and lineage-specific features shape their evolution.

Possible co-option of engrailed during brachiopod and mollusc shell development

In molluscs, two homeobox genes, engrailed (en) and distal-less (dlx), are transcription factors that are expressed in correlation with shell development. They are expressed in the regions between shell-forming and non-shell-forming cells, likely defining the boundaries of shell-forming fields. Here we investigate the expression of two transcription factors in the brachiopod Lingula anatina We find that en is expressed in larval mantle lobes, whereas dlx is expressed in larval tentacles. We also demonstrate that the embryonic shell marker mantle peroxidase (mpox) is specifically expressed in mantle lobes. Our results suggest that en and mpox are possibly involved in brachiopod embryonic shell development. We discuss the evolutionary developmental origin of lophotrochozoan biomineralization through independent gene co-option.

Orthrozanclus elongata n. sp. and the significance of sclerite-covered taxa for early trochozoan evolution

Orthrozanclus is a shell-bearing, sclerite covered Cambrian organism of uncertain taxonomic affinity, seemingly representing an intermediate between its fellow problematica Wiwaxia and Halkieria. Attempts to group these slug-like taxa into a single ‘halwaxiid’ clade nevertheless present structural and evolutionary difficulties. Here we report a new species of Orthrozanclus from the early Cambrian Chengjiang Lagerstätte. The scleritome arrangement and constitution in this material corroborates the link between Orthrozanclus and Halkieria, but not with Wiwaxia — and calls into question its purported relationship with molluscs. Instead, the tripartite construction of the halkieriid scleritome finds a more compelling parallel in the camenellan tommotiids, relatives of the brachiopods and phoronids. Such a phylogenetic position would indicate the presence of a scleritome in the common ancestor of the three major trochozoan lineages, Mollusca, Annelida and Brachiozoa. On this view, the absence of fossil Ediacaran sclerites is evidence against any ‘Precambrian prelude’ to the explosive diversification of these phyla in the Cambrian, c. 540–530 million years ago.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.