Genesis and evolution of the Evx and Mox genes and the extended Hox and ParaHox gene clusters

New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore

BACKGROUND

Early during onychophoran development and prior to the formation of the germ band, a posterior tissue thickening forms the posterior pit. Anterior to this thickening forms a groove, the embryonic slit, that marks the anterior-posterior orientation of the developing embryo. This slit is by some authors considered the blastopore, and thus the origin of the endoderm, while others argue that the posterior pit represents the blastopore. This controversy is of evolutionary significance because if the slit represents the blastopore, then this would support the amphistomy hypothesis that suggests that a slit-like blastopore in the bilaterian ancestor evolved into protostomy and deuterostomy.

RESULTS

In this paper, we summarize our current knowledge about endoderm and mesoderm development in onychophorans and provide additional data on early endoderm- and mesoderm-determining marker genes such as Blimp, Mox, and the T-box genes.

CONCLUSION

We come to the conclusion that the endoderm of onychophorans forms prior to the development of the embryonic slit, and thus that the slit is not the primary origin of the endoderm. It is thus unlikely that the embryonic slit represents the blastopore. We suggest instead that the posterior pit indeed represents the lips of the blastopore, and that the embryonic slit (and surrounding tissue) represents a morphologically superficial archenteron-like structure. We conclude further that both endoderm and mesoderm development are under control of conserved gene regulatory networks, and that many of the features found in arthropods including the model Drosophila melanogaster are likely derived.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

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.

Jellyfish genomes reveal distinct homeobox gene clusters and conservation of small RNA processing

The phylum Cnidaria represents a close outgroup to Bilateria and includes familiar animals including sea anemones, corals, hydroids, and jellyfish. Here we report genome sequencing and assembly for true jellyfish Sanderia malayensis and Rhopilema esculentum. The homeobox gene clusters are characterised by interdigitation of Hox, NK, and Hox-like genes revealing an alternate pathway of ANTP class gene dispersal and an intact three gene ParaHox cluster. The mitochondrial genomes are linear but, unlike in Hydra, we do not detect nuclear copies, suggesting that linear plastid genomes are not necessarily prone to integration. Genes for sesquiterpenoid hormone production, typical for arthropods, are also now found in cnidarians. Somatic and germline cells both express piwi-interacting RNAs in jellyfish revealing a conserved cnidarian feature, and evidence for tissue-specific microRNA arm switching as found in Bilateria is detected. Jellyfish genomes reveal a mosaic of conserved and divergent genomic characters evolved from a shared ancestral genetic architecture.

Solving Classification Problems for Large Sets of Protein Sequences with the Example of Hox and ParaHox Proteins

Phylogenetic methods are key to providing models for how a given protein family evolved. However, these methods run into difficulties when sequence divergence is either too low or too high. Here, we provide a case study of Hox and ParaHox proteins so that additional insights can be gained using a new computational approach to help solve old classification problems. For two (Gsx and Cdx) out of three ParaHox proteins the assignments differ between the currently most established view and four alternative scenarios. We use a non-phylogenetic, pairwise-sequence-similarity-based method to assess which of the previous predictions, if any, are best supported by the sequence-similarity relationships between Hox and ParaHox proteins. The overall sequence-similarities show Gsx to be most similar to Hox2–3, and Cdx to be most similar to Hox4–8. The results indicate that a purely pairwise-sequence-similarity-based approach can provide additional information not only when phylogenetic inference methods have insufficient information to provide reliable classifications (as was shown previously for central Hox proteins), but also when the sequence variation is so high that the resulting phylogenetic reconstructions are likely plagued by long-branch-attraction artifacts.

The origin of the Hox/ParaHox genes, the Ghost Locus hypothesis and the complexity of the first animal

A key aim in evolutionary biology is to deduce ancestral states to better understand the evolutionary origins of clades of interest and the diversification process(es) that has/have elaborated them. These ancestral deductions can hit difficulties when undetected loss events are misinterpreted as ancestral absences. With the ever-increasing amounts of animal genomic sequence data, we are gaining a much clearer view of the preponderance of differential gene losses across animal lineages. This has become particularly clear with recent progress in our understanding of the origins of the Hox/ParaHox developmental control genes relative to the earliest branching lineages of the animal kingdom: the sponges (Porifera), comb jellies (Ctenophora) and placozoans (Placozoa). These reassessments of the diversity and complexity of developmental control genes in the earliest animal ancestors need to go hand-in-hand with complementary advances in comparative morphology, phylogenetics and palaeontology to clarify our understanding of the complexity of the last common ancestor of all animals. The field is currently undergoing a shift from the traditional consensus of a sponge-like animal ancestor from which morphological and molecular elaboration subsequently evolved, to a scenario of a more complex animal ancestor, with subsequent losses and simplifications in various lineages.

The TALE face of Hox proteins in animal evolution

Hox genes are major regulators of embryonic development. One of their most conserved functions is to coordinate the formation of specific body structures along the anterior-posterior (AP) axis in Bilateria. This architectural role was at the basis of several morphological innovations across bilaterian evolution. In this review, we traced the origin of the Hox patterning system by considering the partnership with PBC and Meis proteins. PBC and Meis belong to the TALE-class of homeodomain-containing transcription factors and act as generic cofactors of Hox proteins for AP axis patterning in Bilateria. Recent data indicate that Hox proteins acquired the ability to interact with their TALE partners in the last common ancestor of Bilateria and Cnidaria. These interactions relied initially on a short peptide motif called hexapeptide (HX), which is present in Hox and non-Hox protein families. Remarkably, Hox proteins can also recruit the TALE cofactors by using specific PBC Interaction Motifs (SPIMs). We describe how a functional Hox/TALE patterning system emerged in eumetazoans through the acquisition of SPIMs. We anticipate that interaction flexibility could be found in other patterning systems, being at the heart of the astonishing morphological diversity observed in the animal kingdom.

The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis

Homeobox genes play crucial roles for the development of multicellular eukaryotes. We have generated a revised list of all homeobox genes for Caenorhabditis elegans and provide a nomenclature for the previously unnamed ones. We show that, out of 103 homeobox genes, 70 are co-orthologous to human homeobox genes. 14 are highly divergent, lacking an obvious ortholog even in other Caenorhabditis species. One of these homeobox genes encodes 12 homeodomains, while three other highly divergent homeobox genes encode a novel type of double homeodomain, termed HOCHOB. To understand how transcription factors regulate cell fate during development, precise spatio-temporal expression data need to be obtained. Using a new imaging framework that we developed, Endrov, we have generated spatio-temporal expression profiles during embryogenesis of over 60 homeobox genes, as well as a number of other developmental control genes using GFP reporters. We used dynamic feedback during recording to automatically adjust the camera exposure time in order to increase the dynamic range beyond the limitations of the camera. We have applied the new framework to examine homeobox gene expression patterns and provide an analysis of these patterns. The methods we developed to analyze and quantify expression data are not only suitable for C. elegans, but can be applied to other model systems or even to tissue culture systems.

How Hox genes can shed light on the place of echinoderms among the deuterostomes

Background

The Hox gene cluster ranks among the greatest of biological discoveries of the past 30 years. Morphogenetic patterning genes are remarkable for the systems they regulate during major ontogenetic events, and for their expressions of molecular, temporal, and spatial colinearity. Recent descriptions of exceptions to these colinearities are suggesting deep phylogenetic signal that can be used to explore origins of entire deuterostome phyla. Among the most enigmatic of these deuterostomes in terms of unique body patterning are the echinoderms. However, there remains no overall synthesis of the correlation between this signal and the variations observable in the presence/absence and expression patterns of Hox genes.

Results

Recent data from Hox cluster analyses shed light on how the bizarre shift from bilateral larvae to radial adults during echinoderm ontogeny can be accomplished by equally radical modifications within the Hox cluster. In order to explore this more fully, a compilation of observations on the genetic patterns among deuterostomes is integrated with the body patterning trajectories seen across the deuterostome clade.

Conclusions

Synthesis of available data helps to explain morphogenesis along the anterior/posterior axis of echinoderms, delineating the origins and fate of that axis during ontogeny. From this, it is easy to distinguish between ‘seriality’ along echinoderm rays and true A/P axis phenomena such as colinearity within the somatocoels, and the ontogenetic outcomes of the unique translocation and inversion of the anterior Hox class found within the Echinodermata. An up-to-date summary and integration of the disparate lines of research so far produced on the relationship between Hox genes and pattern formation for all deuterostomes allows for development of a phylogeny and scenario for the evolution of deuterostomes in general, and the Echinodermata in particular.

The origin of patterning systems in bilateria-insights from the Hox and ParaHox genes in Acoelomorpha

Hox and ParaHox genes constitute two families of developmental regulators that pattern the Anterior–Posterior body axis in all bilaterians. The members of these two groups of genes are usually arranged in genomic clusters and work in a coordinated fashion, both in space and in time. While the mechanistic aspects of their action are relatively well known, it is still unclear how these systems evolved. For instance, we still need a proper model of how the Hox and ParaHox clusters were assembled over time. This problem is due to the shortage of information on gene complements for many taxa (mainly basal metazoans) and the lack of a consensus phylogenetic model of animal relationships to which we can relate our new findings. Recently, several studies have shown that the Acoelomorpha most probably represent the first offshoot of the Bilateria. This finding has prompted us, and others, to study the Hox and ParaHox complements in these animals, as well as their activity during development. In this review, we analyze how the current knowledge of Hox and ParaHox genes in the Acoelomorpha is shaping our view of bilaterian evolution.

Expression of the pair-rule gene homologs runt, Pax3/7, even-skipped-1 and even-skipped-2 during larval and juvenile development of the polychaete annelid Capitella teleta does not support a role in segmentation

Background

Annelids and arthropods each possess a segmented body. Whether this similarity represents an evolutionary convergence or inheritance from a common segmented ancestor is the subject of ongoing investigation.

Methods

To investigate whether annelids and arthropods share molecular components that control segmentation, we isolated orthologs of the Drosophila melanogaster pair-rule genes, runt, paired (Pax3/7) and eve, from the polychaete annelid Capitella teleta and used whole mount in situ hybridization to characterize their expression patterns.

Results

When segments first appear, expression of the single C. teleta runt ortholog is only detected in the brain. Later, Ct-runt is expressed in the ventral nerve cord, foregut and hindgut. Analysis of Pax genes in the C. teleta genome reveals the presence of a single Pax3/7 ortholog. Ct-Pax3/7 is initially detected in the mid-body prior to segmentation, but is restricted to two longitudinal bands in the ventral ectoderm. Each of the two C. teleta eve orthologs has a unique and complex expression pattern, although there is partial overlap in several tissues. Prior to and during segment formation, Ct-eve1 and Ct-eve2 are both expressed in the bilaterial pair of mesoteloblasts, while Ct-eve1 is expressed in the descendant mesodermal band cells. At later stages, Ct-eve2 is expressed in the central and peripheral nervous system, and in mesoderm along the dorsal midline. In late stage larvae and adults, Ct-eve1 and Ct-eve2 are expressed in the posterior growth zone.

Conclusions

C. teleta eve, Pax3/7 and runt homologs all have distinct expression patterns and share expression domains with homologs from other bilaterians. None of the pair-rule orthologs examined in C. teleta exhibit segmental or pair-rule stripes of expression in the ectoderm or mesoderm, consistent with an independent origin of segmentation between annelids and arthropods.

Extensive chordate and annelid macrosynteny reveals ancestral homeobox gene organization

Genes with the homeobox motif are crucial in developmental biology and widely implicated in the evolution of development. The Antennapedia (ANTP)-class is one of the two major classes of animal homeobox genes, and includes the Hox genes, renowned for their role in patterning the anterior-posterior axis of animals. The origin and evolution of the ANTP-class genes are a matter of some debate. A principal guiding hypothesis has been the existence of an ancient gene Mega-cluster deep in animal ancestry. This hypothesis was largely established from linkage data from chordates, and the Mega-cluster hypothesis remains to be seriously tested in protostomes. We have thus mapped ANTP-class homeobox genes to the chromosome level in a lophotrochozoan protostome. Our comparison of gene organization in Platynereis dumerilii and chordates indicates that the Mega-cluster, if it did exist, had already been broken up onto four chromosomes by the time of the protostome-deuterostome ancestor (PDA). These results not only elucidate an aspect of the genome organization of the PDA but also reveal high levels of macrosynteny between P. dumerilii and chordates. This implies a very low rate of interchromosomal genome rearrangement in the lineages leading to P. dumerilii and the chordate ancestor since the time of the PDA.

Evolution of Hox complexes

Recent years have seen a plethora of ideas and hypotheses, and lots of debate, about the origin and evolution of the Hox gene cluster. Here I will attempt to summarize these hypotheses, identify their strengths and weaknesses and highlight the types of new data that may lead to further resolution of the competing ideas. The major theme is that Hox genes originated very early in animal evolution and extensive independent duplications occurred in major lineages. Duplications however have not been the only route to change in the composition and structure of the Hox cluster, as extensive gene losses have occurred as well. Indeed it is gene loss that is one of the main obstacles in our understanding of the origin and evolution of Hox clusters. Matters should be improved with wider taxon sampling along with a clearer understanding of how duplicated genes evolve.

How to make a minimal genome for synthetic minimal cell

As a key focus of synthetic biology, building a minimal artificial cell has given rise to many discussions. A synthetic minimal cell will provide an appropriate chassis to integrate functional synthetic parts, devices and systems with functions that cannot generally be found in nature. The design and construction of a functional minimal genome is a key step while building such a cell/chassis since all the cell functions can be traced back to the genome. Kinds of approaches, based on bioinformatics and molecular biology, have been developed and proceeded to derive essential genes and minimal gene sets for the synthetic minimal genome. Experiments about streamlining genomes of model bacteria revealed genome reduction led to unanticipated beneficial properties, such as high electroporation efficiency and accurate propagation of recombinant genes and plasmids that were unstable in other strains. Recent achievements in chemical synthesis technology for large DNA segments together with the rapid development of the whole-genome sequencing, have transferred synthesis of genes to assembly of the whole genomes based on oligonucleotides, and thus created strong preconditions for synthesis of artificial minimal genome. Here in this article, we review briefly the history and current state of research in this field and summarize the main methods for making a minimal genome. We also discuss the impacts of minimized genome on metabolism and regulation of artificial cell.

Are Hox genes ancestrally involved in axial patterning? Evidence from the hydrozoan Clytia hemisphaerica (Cnidaria)

BACKGROUND

The early evolution and diversification of Hox-related genes in eumetazoans has been the subject of conflicting hypotheses concerning the evolutionary conservation of their role in axial patterning and the pre-bilaterian origin of the Hox and ParaHox clusters. The diversification of Hox/ParaHox genes clearly predates the origin of bilaterians. However, the existence of a "Hox code" predating the cnidarian-bilaterian ancestor and supporting the deep homology of axes is more controversial. This assumption was mainly based on the interpretation of Hox expression data from the sea anemone, but growing evidence from other cnidarian taxa puts into question this hypothesis.

METHODOLOGY/PRINCIPAL FINDINGS

Hox, ParaHox and Hox-related genes have been investigated here by phylogenetic analysis and in situ hybridisation in Clytia hemisphaerica, an hydrozoan species with medusa and polyp stages alternating in the life cycle. Our phylogenetic analyses do not support an origin of ParaHox and Hox genes by duplication of an ancestral ProtoHox cluster, and reveal a diversification of the cnidarian HOX9-14 genes into three groups called A, B, C. Among the 7 examined genes, only those belonging to the HOX9-14 and the CDX groups exhibit a restricted expression along the oral-aboral axis during development and in the planula larva, while the others are expressed in very specialised areas at the medusa stage.

CONCLUSIONS/SIGNIFICANCE

Cross species comparison reveals a strong variability of gene expression along the oral-aboral axis and during the life cycle among cnidarian lineages. The most parsimonious interpretation is that the Hox code, collinearity and conservative role along the antero-posterior axis are bilaterian innovations.

Genomic organization and expression demonstrate spatial and temporal Hox gene colinearity in the lophotrochozoan Capitella sp. I

Hox genes define regional identities along the anterior–posterior axis in many animals. In a number of species, Hox genes are clustered in the genome, and the relative order of genes corresponds with position of expression in the body. Previous Hox gene studies in lophotrochozoans have reported expression for only a subset of the Hox gene complement and/or lack detailed genomic organization information, limiting interpretations of spatial and temporal colinearity in this diverse animal clade. We studied expression and genomic organization of the single Hox gene complement in the segmented polychaete annelid Capitella sp. I. Total genome searches identified 11 Hox genes in Capitella, representing 11 distinct paralog groups thought to represent the ancestral lophotrochozoan complement. At least 8 of the 11 Capitella Hox genes are genomically linked in a single cluster, have the same transcriptional orientation, and lack interspersed non-Hox genes. Studying their expression by situ hybridization, we find that the 11 Capitella Hox genes generally exhibit spatial and temporal colinearity. With the exception of CapI-Post1, Capitella Hox genes are all expressed in broad ectodermal domains during larval development, consistent with providing positional information along the anterior–posterior axis. The anterior genes CapI-lab, CapI-pb, and CapI-Hox3 initiate expression prior to the appearance of segments, while more posterior genes appear at or soon after segments appear. Many of the Capitella Hox genes have either an anterior or posterior expression boundary coinciding with the thoracic–abdomen transition, a major body tagma boundary. Following metamorphosis, several expression patterns change, including appearance of distinct posterior boundaries and restriction to the central nervous system. Capitella Hox genes have maintained a clustered organization, are expressed in the canonical anterior–posterior order found in other metazoans, and exhibit spatial and temporal colinearity, reflecting Hox gene characteristics that likely existed in the protostome–deuterostome ancestor.

Statistical tests between competing hypotheses of Hox cluster evolution

The Hox genes encode transcription factors that play vital roles in the anterior-posterior patterning of all bilaterian phyla studied to date. Additionally, the gain of Hox genes by duplication has been widely implicated as a driving force in the evolution of animal body plans. Because of this, reconstructing the evolution of the Hox cluster has been the focus of intense research interest. It has been commonly assumed that an ancestral four-gene ProtoHox cluster was duplicated early in animal evolution to give rise to the Hox and ParaHox clusters. However, this hypothesis has recently been called into question, and a number of alternative hypotheses of Hox and ParaHox gene evolution have been proposed. Here, we present the first statistical comparisons of current hypotheses of Hox and ParaHox gene evolution. We use two statistical methods that represent two different approaches to the treatment of phylogenetic uncertainty. In the first method, we estimate the maximum-likelihood tree for each hypothesis and compare these trees to one another using a parametric bootstrapping approach. In the second method, we use Bayesian phylogenetics to estimate the posterior distribution of trees, then we calculate the support for each hypothesis from this distribution. The results of both methods are largely congruent. We find that we are able to reject five out of the eight current hypotheses of Hox and ParaHox gene evolution that we consider. We conclude that the ProtoHox cluster is likely to have contained either three or four genes but that there is insufficient phylogenetic signal in the homeodomains to distinguish between these alternatives.

A novel conserved evx1 enhancer links spinal interneuron morphology and cis-regulation from fish to mammals

Spinal interneurons are key components of locomotor circuits, driving such diverse behaviors as swimming in fish and walking in mammals. Recent work has linked the expression of evolutionarily conserved transcription factors to key features of interneurons in diverse species, raising the possibility that these interneurons are functionally related. Consequently, the determinants of interneuron subtypes are predicted to share conserved cis-regulation in vertebrates with very different spinal cords. Here, we establish a link between cis-regulation and morphology of spinal interneurons that express the Evx1 homeodomain transcription factor from fish to mammals. Using comparative genomics, and complementary transgenic approaches, we have identified a novel enhancer of evx1, that includes two non-coding elements conserved in vertebrates. We show that pufferfish evx1 transgenes containing this enhancer direct reporter expression to a subset of spinal commissural interneurons in zebrafish embryos. Pufferfish, zebrafish and mouse evx1 downstream genomic enhancers label selectively Evx1(+) V0 commissural interneurons in chick and rat embryos. By dissecting the zebrafish evx1 enhancer, we identify a role for a 25 bp conserved cis-element in V0-specific gene expression. Our findings support the notion that spinal interneurons shared between distantly related vertebrates, have been maintained in part via the preservation of highly conserved cis-regulatory modules.

Testing gene set enrichment for subset of genes: Sub-GSE

Background

Many methods have been developed to test the enrichment of genes related to certain phenotypes or cell states in gene sets. These approaches usually combine gene expression data with functionally related gene sets as defined in databases such as GeneOntology (GO), KEGG, or BioCarta. The results based on gene set analysis are generally more biologically interpretable, accurate and robust than the results based on individual gene analysis. However, while most available methods for gene set enrichment analysis test the enrichment of the entire gene set, it is more likely that only a subset of the genes in the gene set may be related to the phenotypes of interest.

Results

In this paper, we develop a novel method, termed Sub-GSE, which measures the enrichment of a predefined gene set, or pathway, by testing its subsets. The application of Sub-GSE to two simulated and two real datasets shows Sub-GSE to be more sensitive than previous methods, such as GSEA, GSA, and SigPath, in detecting gene sets assiated with a phenotype of interest. This is particularly true for cases in which only a fraction of the genes in the gene set are associated with the phenotypes. Furthermore, the application of Sub-GSE to two real data sets demonstrates that it can detect more biologically meaningful gene sets than GSEA.

Conclusion

We developed a new method to measure the gene set enrichment. Applications to two simulated datasets and two real datasets show that this method is sensitive to the associations between gene sets and phenotype. The program Sub-GSE can be downloaded from .

Global control regions and regulatory landscapes in vertebrate development and evolution

During the course of evolution, many genes that control the development of metazoan body plans were co-opted to exert novel functions, along with the emergence or modification of structures. Gene amplification and/or changes in the cis-regulatory modules responsible for the transcriptional activity of these genes have certainly contributed in a major way to evolution of gene functions. In some cases, these processes led to the formation of groups of adjacent genes that appear to be controlled by both global and shared mechanisms.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications

Background

The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes.

Results

We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons.

Conclusion

There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular – but possibly clusters of genes more generally – might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Ancient complexity of the non-Hox ANTP gene complement in the anthozoan Nematostella vectensis: implications for the evolution of the ANTP superclass

The origin and evolution of ANTP superclass genes has raised controversial discussions. While recent evidence suggests that a true Hox cluster emerged after the cnidarian bilaterian split, the origin of the ANTP superclass as a whole remains unclear. Based on analyses of bilaterian genomes, it seems very likely that clustering has once been a characteristic of all ANTP homeobox genes and that their ancestors have emerged through several series of cis-duplications from the same genomic region. Since the diploblastic Cnidaria possess orthologs of some non-Hox ANTP genes, at least some steps of the expansion of this hypothetical homeobox gene array must have occurred in the last common ancestor of both lineages--but it is unknown to what extent. By screening the unassembled Nematostella genome, we have identified unambiguous orthologs to almost all non-Hox ANTP genes which are present in Bilateria--with the exception of En, Tlx and (possibly) Vax. Furthermore, Nematostella possesses ANTP genes that are missing in some bilaterian lineages, like the rough gene or NK7. In addition, several ANTP homeobox gene families have been independently duplicated in Nematostella. We conclude that the last cnidarian/bilaterian ancestor already harboured the almost full complement of non-Hox ANTP genes before the Hox system evolved.

Pre-bilaterian origins of the Hox cluster and the Hox code: evidence from the sea anemone, Nematostella vectensis

BACKGROUND

Hox genes were critical to many morphological innovations of bilaterian animals. However, early Hox evolution remains obscure. Phylogenetic, developmental, and genomic analyses on the cnidarian sea anemone Nematostella vectensis challenge recent claims that the Hox code is a bilaterian invention and that no "true" Hox genes exist in the phylum Cnidaria.

METHODOLOGY/PRINCIPAL FINDINGS

Phylogenetic analyses of 18 Hox-related genes from Nematostella identify putative Hox1, Hox2, and Hox9+ genes. Statistical comparisons among competing hypotheses bolster these findings, including an explicit consideration of the gene losses implied by alternate topologies. In situ hybridization studies of 20 Hox-related genes reveal that multiple Hox genes are expressed in distinct regions along the primary body axis, supporting the existence of a pre-bilaterian Hox code. Additionally, several Hox genes are expressed in nested domains along the secondary body axis, suggesting a role in "dorsoventral" patterning.

CONCLUSIONS/SIGNIFICANCE

A cluster of anterior and posterior Hox genes, as well as ParaHox cluster of genes evolved prior to the cnidarian-bilaterian split. There is evidence to suggest that these clusters were formed from a series of tandem gene duplication events and played a role in patterning both the primary and secondary body axes in a bilaterally symmetrical common ancestor. Cnidarians and bilaterians shared a common ancestor some 570 to 700 million years ago, and as such, are derived from a common body plan. Our work reveals several conserved genetic components that are found in both of these diverse lineages. This finding is consistent with the hypothesis that a set of developmental rules established in the common ancestor of cnidarians and bilaterians is still at work today.

Ancient linkage of a POU class 6 and an anterior hox-like gene in cnidaria: implications for the evolution of homeobox genes

Linkage analyses in metazoan genomes suggest two ancestral arrays for the majority of homeobox genes. The related homeobox genes and chromosomal regions that are dispersed in extant species derived possibly from only two single common ancestor regions. One proposed ancestral array, designated as ANTP mega-array, contains most of the ANTP class homeobox genes; the second, named the contraHox super-paralogon, would consist of the classes PRD, POU, LIM, CUT, prospero, TALE and SIX. Here, we report the tight linkage of a POU class 6 gene to an anterior Hox-like gene in the hydrozoan Eleutheria dichotoma and discuss its possible significance for the evolution of homeobox genes. POU class 6 genes also seem to be ancestrally linked to the HoxC and A clusters in vertebrates, despite POU homeobox genes belonging to the contraHox paralogon. Hence, the much tighter linkage of a POU class 6 gene to an anterior Hox-like gene in a cnidarian is possibly the evolutionary echo of an ancestral genomic region from which most metazoan homeobox classes emerged.

Origin and Molecular Evolution of Receptor Tyrosine Kinases with Immunoglobulin-Like Domains

Receptor tyrosine kinases (RTKs) are involved in the control of fundamental cellular processes in metazoans. In vertebrates, RTK could be grouped in distinct classes based on the nature of their cognate ligand and modular composition of their extracellular domain. RTK with immunoglobulin-like domains (IG-like RTK) encompass several RTK classes and have been found in early metazoans, including sponges. Evolution of IG-like RTK is characterized by extended molecular and functional diversification, which prompted us to study their evolutionary history. For that purpose, a nonredundant data set including annotated protein sequences of IG-like RTK (n = 85) was built, representing 19 species ranging from sponges to humans. Phylogenetic trees were generated from alignment of conserved regions using maximum likelihood approach. Molecular phylogeny strongly suggests that IG-like RTK diversification occurred according to a complex scenario. In particular, we propose that specific cis duplications of a common ancestor to both platelet-derived growth factor receptor (class III) and vascular endothelial growth factor receptor (class V) families preceded two trans duplications. In contrast, other IG-like RTK genes, like Musk and PTK7, apparently did not evolve by duplications, whereas fibroblast growth factor receptors (class IV) evolved through two rounds of trans duplications. The proposed model of IG-like RTK evolution is supported by high bootstrap values and by the clustering of genes encoding class III and class V RTKs at specific chromosomal locations in mouse and human genomes.

Gene duplication, gene loss and evolution of expression domains in the vertebrate nuclear receptor NR5A (Ftz-F1) family

Fushi tarazu factor 1 (Ftz-F1, NR5A) is a zinc-finger transcription factor that belongs to the nuclear receptor superfamily and regulates genes that are involved in sterol and steroid metabolism in gonads, adrenals, liver and other tissues. To understand the evolutionary origins and developmental genetic relationships of the Ftz-F1 genes, we have cloned four homologous Ftz-f1 genes in zebrafish, called ff1a, ff1b, ff1c and ff1d. These four genes have different temporal and spatial expression patterns during development, indicating that they have distinct mechanisms of genetic regulation. Among them, the ff1a expression pattern is similar to mammalian Nr5a2, while the ff1b pattern is similar to that of mammalian Nr5a1. Genetic mapping experiments show that these four ff1 genes are located on chromosome segments conserved between the zebrafish and human genomes, indicating a common ancestral origin. Phylogenetic and conserved synteny analysis show that ff1a is the orthologue of NR5A2, and that ff1b and ff1d genes are co-orthologues of NR5A1 that arose by a gene-duplication event, probably a whole-genome duplication, in the ray-fin lineage, and each gene is located next to an NR6A1 co-orthologue as in humans, showing that the tandem duplication occurred before the divergence of human and zebrafish lineages. ff1c does not have a mammalian counterpart. Thus we have characterized the phylogenetic relationships, expression patterns and chromosomal locations of these Ftz-F1 genes, and have demonstrated their identities as NR5A genes in relation to the orthologous genes in other species.

The genesis and evolution of homeobox gene clusters

Once called the 'Rosetta stone' of developmental biology, the homeobox continues to fascinate both evolutionary and developmental biologists. The birth of the homeotic, or Hox, gene cluster, and its subsequent evolution, has been crucial in mediating the major transitions in metazoan body plan. Comparative genomics studies indicate that the more recently discovered ParaHox and NK clusters were linked to the Hox cluster early in evolution, and that together they constituted a 'megacluster' of homeobox genes that conspicuously contributed to body-plan evolution.

Pdx1-related homeodomain transcription factors are distinctly expressed in mouse adult pancreatic islets

Complex gene networks are responsible for the proper operation of the endocrine pancreas. A central member of such networks, the homeodomain transcription factor Pdx1, belongs to the ParaHox gene cluster, an array of Hox-like homeobox genes. With a combination of mRNA in situ hybridisation and immunodetection, we have found that the rest of ParaHox cluster genes, Cdx1, Cdx2/3, and Cdx4, and Gsh1 and Gsh2, are all expressed in specific islet cell types of the endocrine pancreas. To our knowledge, this is the first report that locates ParaHox genes other than Pdx1 and Cdx2/3 in a place as to be involved in the pancreatic transcriptional regulatory networks, potentially regulating glucagon-insulin homeostasis.

No more than 14: the end of the amphioxus Hox cluster

The Hox gene cluster has been a key paradigm for a generation of developmental and evolutionary biologists. Since its discovery in the mid-1980's, the identification, genomic organization, expression, colinearity, and regulation of Hox genes have been immediate targets for study in any new model organism, and metazoan genome projects always refer to the structure of the particular Hox cluster(s). Since the early 1990's, it has been dogma that vertebrate Hox clusters are composed of thirteen paralogous groups. Nonetheless, we showed that in the otherwise prototypical cephalochordate amphioxus (Branchiostoma floridae), the Hox cluster contains a fourteenth Hox gene, and very recently, a 14^th Hox paralogous group has been found in the coelacanth and the horn shark, suggesting that the amphioxus cluster was anticipating the finding of Hox 14 in some vertebrate lineages. In view of the pivotal place that amphioxus occupies in vertebrate evolution, we thought it of considerable interest to establish the limits of its Hox gene cluster, namely resolution of whether more Hox genes are present in the amphioxus cluster (e.g., Hox 15). Using two strategies, here we report the completion and characterization of the Hox gene content of the single amphioxus Hox cluster, which encompasses 650 kb from Hox1 to Evx. Our data have important implications for the primordial Hox gene cluster of chordates: the prototypical nature of the single amphioxus Hox cluster makes it unlikely that additional paralogous groups will be found in any chordate lineage. We suggest that 14 is the end.

Hox, ParaHox, ProtoHox: facts and guesses

The Hox gene cluster has captivated the imagination of evolutionary and developmental biologists worldwide. In this review, the origin of the Hox and ParaHox gene clusters by duplication of a ProtoHox gene cluster, and the changes in their gene numbers in major Metazoan Transitions are reviewed critically. Re-evaluation of existing data and recent findings in Cnidarians, Acoels, and critical stages of vertebrate evolution suggest alternative scenarios for the origin, structure, and changes in Hox gene numbers in relevant events of Metazoan evolution. I discuss opposing views and propose that (i) the ProtoHox cluster had only two genes, and not four as commonly believed: a corollary is that the origin of Bilaterians was coincident with the invention of new Hox and ParaHox gene classes, which may have facilitated such a transition; (ii) the ProtoHox cluster duplication was a cis duplication event, rather than a trans duplication event, as previously suggested, and (iii) the ancestral vertebrate cluster possessed 14 Hox genes, and not the 13 generally assumed. These hypotheses could be verified or refuted in the near future, but they may help critical discussion of the evolution of the Hox/ParaHox family in the metazoan kingdom.

Ciona intestinalis Hox gene cluster: Its dispersed structure and residual colinear expression in development

Ascidians, belonging to the subphylum Urochordata, the earliest branch from the lineage to the vertebrates, exhibit a prototypical morphogenesis of chordates in the larval development, although they subsequently metamorphose into adults with a unique body structure. Recent draft genome analysis of the ascidian Ciona intestinalis has identified 9 Hox genes, which, however, have been located on five scaffolds. Similarly, expression patterns of Ciona Hox genes are largely unknown, although some data have been available for a few Hox member genes. Thus, the cluster structure and colinearity of Hox genes are still an enigma in C. intestinalis. To address these issues, we used fluorescence in situ hybridization and whole-mount in situ hybridization techniques and examined the genomic organization and spatiotemporal expression of all Hox as well as extended Hox member genes (Evx and Mox) of C. intestinalis. We found that seven of nine Ciona Hox genes are located on a single chromosome with some ordering exceptions, and the other genes, including Evx and Mox, are located on three other chromosomes. Some Ciona Hox genes, if not all, exhibited spatially coordinated expression within the larval central nervous system and the gut of the juvenile. In light of these observations, we suggest that the cluster organization and colinearity of the Hox genes are under dispersing and disintegrating conditions in C. intestinalis.

The Hox gene complement of acoel flatworms, a basal bilaterian clade

Several molecular data sets suggest that acoelomorph flatworms are not members of the phylum Platyhelminthes but form a separate branch of the Metazoa that diverged from all other bilaterian animals before the separation of protostomes and deuterostomes. Here we examine the Hox gene complement of the acoel flatworms. In two distantly related acoel taxa, we identify only three distinct classes of Hox gene: an anterior gene, a posterior gene, and a central class gene most similar to genes of Hox classes 4 and 5 in other Bilateria. Phylogenetic analysis of these genes, together with the acoel caudal homologue, supports the basal position of the acoels. The similar gene sets found in two distantly related acoels suggest that this reduced gene complement may be ancestral in the acoels and that the acoels may have diverged from other bilaterians before elaboration of the 8- to 10-gene Hox cluster that characterizes most bilaterians.