Ancient homeobox gene loss and the evolution of chordate brain and pharynx development: deductions from amphioxus gene expression

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Gene Loss may have Shaped the Cnidarian and Bilaterian Hox and ParaHox Complement

Hox and ParaHox transcription factors are important for specifying cell fates along the primary body axes during the development of most animals. Within Cnidaria, much of the research on Hox/ParaHox genes has focused on Anthozoa (anemones and corals) and Hydrozoa (hydroids) and has concentrated on the evolution and function of cnidarian Hox genes in relation to their bilaterian counterparts. Here we analyze together the full complement of Hox and ParaHox genes from species representing all four medusozoan classes (Staurozoa, Cubozoa, Hydrozoa, and Scyphozoa) and both anthozoan classes (Octocorallia and Hexacorallia). Our results show that Hox genes involved in patterning the directive axes of anthozoan polyps are absent in the stem leading to Medusozoa. For the first time, we show spatial and temporal expression patterns of Hox and ParaHox genes in the upside-down jellyfish Cassiopea xamachana (Scyphozoa), which are consistent with diversification of medusozoan Hox genes both from anthozoans and within medusozoa. Despite unprecedented taxon sampling, our phylogenetic analyses, like previous studies, are characterized by a lack of clear homology between most cnidarian and bilaterian Hox and Hox-related genes. Unlike previous studies, we propose the hypothesis that the cnidarian-bilaterian ancestor possessed a remarkably large Hox complement and that extensive loss of Hox genes was experienced by both cnidarian and bilaterian lineages.

Hox Gene Collinearity May Be Related to Noether Theory on Symmetry and Its Linked Conserved Quantity

Hox Gene Collinearity (HGC) is a fundamental property that controls the development of many animal species, including vertebrates. In the Hox gene clusters, the genes are located in a sequential order Hox1, Hox2, Hox3, etc., along the 3’ to 5’ direction of the cluster in the chromosome. During Hox cluster activation, the Hox genes are expressed sequentially in the ontogenetic units D1, D2, D3, etc., along the anterior–posterior axis (A-P) of the early embryo. This collinearity, first observed by E.B. Lewis, is surprising because the spatial collinearity of these structures (Hox clusters and embryos) correlates entities that differ by about four orders of magnitude. Biomolecular mechanisms alone cannot explain such correlations. Long-range physical interactions, such as diffusion or electric attractions, should be involved. A biophysical model (BM) was formulated, which, in alignment with the biomolecular processes, successfully describes the existing vertebrate genetic engineering data. One hundred years ago, Emmy Noether made a fundamental discovery in mathematics and physics. She proved, rigorously, that a physical system obeying a symmetry law (e.g., rotations or self-similarity) is followed by a conserved physical quantity. It is argued here that HGC obeys a ‘primitive’ self-similarity symmetry. In this case, the associated primitive conserved quantity is the irreversibly increasing ‘ratchet’-like Hoxgene ordering where some genes may be missing. The genes of a vertebrate Hox clusterare located along a finite straight line. The same order follows the ontogenetic unitsof the vertebrate embryo. Therefore, HGC is a manifestation of a primitive Noether Theory (NT). NT may be applied to other than the vertebrate case, for instance, to animals with a circular topological symmetry. For example, the observed abnormal Hox gene ordering of the echinoderm Hox clusters may be reproduced by a double-strand break of the circular Hox gene ordering and its subsequent incorporation in the flanking chromosome.

Expression of NK genes that are not part of the NK cluster in the onychophoran Euperipatoides rowelli (Peripatopsidae)

Background

NK genes are a group of homeobox transcription factors that are involved in various molecular pathways across bilaterians. They are typically divided into two subgroups, the NK cluster (NKC) and NK-linked genes (NKL). While the NKC genes have been studied in various bilaterians, corresponding data of many NKL genes are missing to date. To further investigate the ancestral roles of NK family genes, we analyzed the expression patterns of NKL genes in the onychophoran Euperipatoides rowelli.

Results

The NKL gene complement of E. rowelli comprises eight genes, including BarH, Bari, Emx, Hhex, Nedx, NK2.1, vax and NK2.2, of which only NK2.2 was studied previously. Our data for the remaining seven NKL genes revealed expression in different structures associated with the developing nervous system in embryos of E. rowelli. While NK2.1 and vax are expressed in distinct medial regions of the developing protocerebrum early in development, BarH, Bari, Emx, Hhex and Nedx are expressed in late developmental stages, after all major structures of the nervous system have been established. Furthermore, BarH and Nedx are expressed in distinct mesodermal domains in the developing limbs.

Conclusions

Comparison of our expression data to those of other bilaterians revealed similar patterns of NK2.1, vax, BarH and Emx in various aspects of neural development, such as the formation of anterior neurosecretory cells mediated by a conserved molecular mechanism including NK2.1 and vax, and the development of the central and peripheral nervous system involving BarH and Emx. A conserved role in neural development has also been reported from NK2.2, suggesting that the NKL genes might have been primarily involved in neural development in the last common ancestor of bilaterians or at least nephrozoans (all bilaterians excluding xenacoelomorphs). The lack of comparative data for many of the remaining NKL genes, including Bari, Hhex and Nedx currently hampers further evolutionary conclusions. Hence, future studies should focus on the expression of these genes in other bilaterians, which would provide a basis for comparative studies and might help to better understand the role of NK genes in the diversification of bilaterians.

Electronic supplementary material

The online version of this article (10.1186/s12861-019-0185-9) contains supplementary material, which is available to authorized users.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

The Nodal signaling pathway controls left-right asymmetric development in amphioxus

Background

Nodal is an important determinant of the left-right (LR) body axis in bilaterians, specifying the right side in protostomes and non-chordate deuterostomes as opposed to the left side in chordates. Amphioxus represents an early-branching chordate group, rendering it especially useful for studying the character states that predate the origin of vertebrates. However, its anatomy, involving offset arrangement of axial structures, marked asymmetry of the oropharyngeal region, and, most notably, a mouth positioned on the left side, contrasts with the symmetric arrangement of the corresponding regions in other chordates.

Results

We show that the Nodal signaling pathway acts to specify the LR axis in the cephalochordate amphioxus in a similar way as in vertebrates. At early neurula stages, Nodal switches from initial bilateral to the left-sided expression and subsequently specifies the left embryonic side. Perturbation of Nodal signaling with small chemical inhibitors (SB505124 and SB431542) alters expression of other members of the pathway and of left/right-sided, organ-specific genes. Upon inhibition, larvae display loss of the innate alternation of both somites and axons of peripheral nerves and loss of left-sided pharyngeal structures, such as the mouth, the preoral pit, and the duct of the club-shaped gland. Concomitantly, the left side displays ectopic expression of otherwise right-sided genes, and the larvae exhibit bilaterally symmetrical morphology, with duplicated endostyle and club-shaped gland structures.

Conclusions

We demonstrate that Nodal signaling is necessary for establishing the LR embryonic axis and for developing profound asymmetry in amphioxus. Our data suggest that initial symmetry breaking in amphioxus and propagation of the pathway on the left side correspond with the situation in vertebrates. However, the organs that become targets of the pathway differ between amphioxus and vertebrates, which may explain the pronounced asymmetry of its oropharyngeal and axial structures and the left-sided position of the mouth.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-6-5) contains supplementary material, which is available to authorized users.

The logic of gene regulatory networks in early vertebrate forebrain patterning

The vertebrate forebrain or prosencephalon is patterned at the beginning of neurulation into four major domains: the telencephalic, hypothalamic, retinal and diencephalic anlagen. These domains will then give rise to the majority of the brain structures involved in sensory integration and the control of higher intellectual and homeostatic functions. Understanding how forebrain pattering arises has thus attracted the interest of developmental neurobiologists for decades. As a result, most of its regulators have been identified and their hierarchical relationship is now the object of active investigation. Here, we summarize the main morphogenetic pathways and transcription factors involved in forebrain specification and propose the backbone of a possible gene regulatory network (GRN) governing its specification, taking advantage of the GRN principles elaborated by pioneer studies in simpler organisms. We will also discuss this GRN and its operational logic in the context of the remarkable morphological and functional diversification that the forebrain has undergone during evolution.

A genome-wide view of transcription factor gene diversity in chordate evolution: less gene loss in amphioxus?

Previous studies of gene diversity in the homeobox superclass have shown that the Florida amphioxus Branchiostoma floridae has undergone remarkably little gene family loss. Here we use a combined BLAST and HMM search strategy to assess the family level diversity of four other transcription factor superclasses: the Paired/Pax genes, Tbx genes, Fox genes and Sox genes. We apply this across genomes from five chordate taxa, including B. floridae and Ciona intestinalis, plus two outgroup taxa. Our results show scattered gene family loss. However, as also found for homeobox genes, B. floridae has retained all ancient Pax, Tbx, Fox and Sox gene families that were present in the common ancestor of living chordates. We conclude that, at least in terms of transcription factor gene complexity, the genome of amphioxus has experienced remarkable stasis compared to the genomes of other chordates.

Pikaia gracilens Walcott, a stem-group chordate from the Middle Cambrian of British Columbia

The Middle Cambrian Pikaia gracilens (Walcott) has an iconic position as a Cambrian chordate, but until now no detailed description has been available. Here on the basis of the 114 available specimens we review its anatomy, confirm its place in the chordates and explore with varying degrees of confidence its relationships to both extant and extinct chordates and other deuterostomes. The body of Pikaia is fusiform, laterally compressed and possesses about 100 myomeres. The head is small, bilobed and bears two narrow tentacles. There is no evidence for eyes. Apart from a thin dorsal fin (without finrays) and a series of at least nine bilaterally arranged appendages with possible pharyngeal pores at the anterior end, there are no other external features. In addition to the musculature the internal anatomy includes an alimentary canal, the anterior of which forms a prominent lenticular unit that is almost invariably preserved in positive relief. The cavity is interpreted as pharyngeal, implying that the mouth itself was almost terminal. The posterior extension of the gut is unclear although the anus appears to have been terminal. The most prominent internal structure is a reflectively preserved unit, possibly hollow, termed here the dorsal organ. Although formerly interpreted as a notochord its position and size make this less likely. Its original function remains uncertain, but it could have formed a storage organ. Ventral to the dorsal organ a narrower strand of tissue is interpreted as representing the nerve chord and notochord. In addition to these structures, there is also evidence for a vascular system, including a ventral blood vessel. The position of Pikaia in the chordates is largely based on the presence of sigmoidal myomeres, and the more tentative identification of a notochord. In many other respects, Pikaia differs from the expected nature of primitive chordates, especially as revealed in amphioxus and the Cambrian record (including Cathaymyrus, Haikouichthys, Metaspriggina, Myllokunmingia, and Zhongxiniscus). Whilst the possibility that Pikaia is simply convergent on the chordates cannot be dismissed, we prefer to build a scenario that regards Pikaia as the most stem-ward of the chordates with links to the phylogenetically controversial yunnanozoans. This hypothesis has implications for the evolution of the myomeres, notochord and gills. Finally, the wealth of material of Pikaia indicates that, although by definition there must be some sort of taphonomic imprint, the consistency of preservational details allows a reliable reconstruction of the anatomy and does not significantly erode phylogenetically relevant characters.