Early steps in the evolution of multicellularity: deep structural and functional homologies among homeobox genes in sponges and higher metazoans

The sponge homeobox gene EmH-3 had not been attributed to any homeobox family. Comparative promoter and homeodomain sequence analyses suggest that it is related to the Hox11 gene, which belongs to the Tlx homeobox family. Hox11 is highly expressed in proliferating progenitor cells, but expression is downregulated during cell differentiation. Using reporter gene methodology, we monitored function of the sponge EmH-3 promoter transfected into human erythroleukemia K562 cells. These cells express the Tlx/Hox11 gene constitutively, and downregulate its expression upon differentiation. The same pattern of expression and downregulation was observed for the sponge reporter construct. We propose that Tlx/Hox11 genes have structural and functional homologies conserved in phylogenetically distant groups, that represent a deep homology in the regulation of cell proliferation, commitment and differentiation.     

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.     

Isolation of Hox genes from the scyphozoan Cassiopeia xamachana: implications for the early evolution of Hox genes.

The isolation of Hox genes from two cnidarian groups, the Hydrozoa and Anthozoa, has sparked hypotheses on the early evolution of Hox genes and a conserved role for these genes for defining a main body axis in all metazoan animals. We have isolated the first five Hox genes, Scox-1 to Scox-5, from the third cnidarian class, the Scyphozoa. For all but one gene, we report full-length homeobox plus flanking sequences. Four of the five genes show close relationship to previously reported Cnox-1 genes from Hydrozoa and Anthozoa. One gene, Scox-2, is an unambiguous homologue of Cnox-2 genes known from Hydrozoa, Anthozoa, and also Placozoa. Based on sequence similarity and phylogenetic analyses of the homeobox and homeodomain sequences of known Hox genes from cnidarians, we suggest the presence of at least five distinct Hox gene families in this phylum, and conclude that the last common ancestor of the Recent cnidarian classes likely possessed a set of Hox genes representing three different families, the Cnox-1, Cnox-2, and Cnox-5 families. The data presented are consistent with the idea that multiple duplication events of genes have occurred within one family at the expense of conservation of the original set of genes, which represent the three ancestral Hox gene families.     

Isolation of Hox cluster genes from insects reveals an accelerated sequence evolution rate

Among gene families it is the Hox genes and among metazoan animals it is the insects (Hexapoda) that have attracted particular attention for studying the evolution of development. Surprisingly though, no Hox genes have been isolated from 26 out of 35 insect orders yet, and the existing sequences derive mainly from only two orders (61% from Hymenoptera and 22% from Diptera). We have designed insect specific primers and isolated 37 new partial homeobox sequences of Hox cluster genes (lab, pb, Hox3, ftz, Antp, Scr, abd-a, Abd-B, Dfd, and Ubx) from six insect orders, which are crucial to insect phylogenetics. These new gene sequences provide a first step towards comparative Hox gene studies in insects. Furthermore, comparative distance analyses of homeobox sequences reveal a correlation between gene divergence rate and species radiation success with insects showing the highest rate of homeobox sequence evolution.     

Molecular cloning and developmental expression of Tlx (Hox11) genes in zebrafish (Danio rerio)

Tlx (Hox11) genes are orphan homeobox genes that play critical roles in the regulation of early developmental processes in vertebrates. Here, we report the identification and expression patterns of three members of the zebrafish Tlx family. These genes share similar, but not identical, expression patterns with other vertebrate Tlx-1 and Tlx-3 genes. Tlx-1 is expressed early in the developing hindbrain and pharyngeal arches, and later in the putative splenic primordium. However, unlike its orthologues, zebrafish Tlx-1 is not expressed in the cranial sensory ganglia or spinal cord. Two homologues of Tlx-3 were identified: Tlx-3a and Tlx-3b, which are both expressed in discrete regions of the developing nervous system, including the cranial sensory ganglia and Rohon-Beard neurons. However, only Tlx-3a is expressed in the statoacoustic cranial ganglia, enteric neurons and non-neural tissues such as the fin bud and pharyngeal arches and Tlx-3b is only expressed in the dorsal root ganglia.     

Phylogenetic relationships of the marine Haplosclerida (Phylum Porifera) employing ribosomal (28S rRNA) and mitochondrial (cox1, nad1) gene sequence data

The systematics of the poriferan Order Haplosclerida (Class Demospongiae) has been under scrutiny for a number of years without resolution. Molecular data suggests that the order needs revision at all taxonomic levels. Here, we provide a comprehensive view of the phylogenetic relationships of the marine Haplosclerida using many species from across the order, and three gene regions. Gene trees generated using 28S rRNA, nad1 and cox1 gene data, under maximum likelihood and Bayesian approaches, are highly congruent and suggest the presence of four clades. Clade A is comprised primarily of species of Haliclona and Callyspongia, and clade B is comprised of H. simulans and H. vansoesti (Family Chalinidae), Amphimedon queenslandica (Family Niphatidae) and Tabulocalyx (Family Phloeodictyidae), Clade C is comprised primarily of members of the Families Petrosiidae and Niphatidae, while Clade D is comprised of Aka species. The polyphletic nature of the suborders, families and genera described in other studies is also found here.     

Hox genes from the Polystomatidae (Platyhelminthes, Monogenea)

Hox genes form a multigenic family that play a fundamental role during the early stages of development. They are organised in a single cluster and share a 60 amino acid conserved sequence that corresponds to the DNA binding domain, i.e. the homeodomain. Sequence conservation in this region has allowed investigators to explore Hox diversity in the metazoan lineages. Within parasitic flatworms only homeobox sequences of parasite species from the Cestoda and Digenea have been reported. In the present study we surveyed species of the Polyopisthocotylea (Monogenea) in order to clarify Hox identification and diversification processes in the neodermatan lineage. From cloning of degenerative PCR products of the central region of the homeobox, we report one ParaHox and 25 new Hox sequences from 10 species of the Polystomatidae and one species of the Diclidophoridae, which extend Hox gene diversity from 46 to 72 within Neodermata. Hox sequences from the Polyopisthocotylea were annotated and classified from sequence alignments and Bayesian inferences of 178 Hox, ParaHox and related gene families recovered from all available parasitic platyhelminths and other bilaterian taxa. Our results are discussed in the light of the recent Hox evolutionary schemes. They may provide new perspectives to study the transition from turbellarians to parasitic flatworms with complex life-cycles and outline the first steps for evolutionary developmental biological approaches within platyhelminth parasites.     

Evidence for 14 homeobox gene clusters in human genome ancestry

The arrangement of Hox genes into physical clusters is fundamental to the patterning of animal body plans. Other homeobox genes are often described as dispersed, with only occasional examples of linkage reported, such as the amphioxus ParaHox and Drosophila 93D/E clusters. This clustering is unlikely to be the derived condition, as the genes of the ParaHox and 93D/E clusters are phylogenetically widespread. To assess whether clustering is retained in mammals, and to infer its history, we considered the distribution of ANTP superclass homeobox genes in human and mouse genomes. We postulate four ancient arrays of ANTP superclass genes in animal genomes, denoted 'extended Hox' (Hox, Evx and Mox), NKL (including NK1, NK3, NK4, Lbx, Tlx, Emx, Vax, Hmx, NK6, Msx), ParaHox (Cdx, Xlox, Gsx) and EHGbox (En, HB9, Gbx). Each of these duplicated in the ancestry of the human genome to yield four Hox, four NKL, four ParaHox and at least two EHGbox clusters or arrays. Two of the human NKL clusters (four in mouse) have subsequently been split by chromosome rearrangement, as has one human EHGbox array. We date all cluster duplications to early chordate evolution and infer that three clusters (Hox, NKL, EHGbox) resided on the same chromosome before duplication.     

The freshwater planarian Dugesia (G.) tigrina contains a great diversity of homeobox genes

To identify potential pattern control and cell determination and/or differentiation genes in the freshwater planarian Dugesial (G.) tigrina, we searched for homeobox genes of different types in the genome of this primitive metazoan. We applied two basic approaches: 1) Screening the cDNA library with degenerate oligonucleotides corresponding to the most conserved amino acid sequence from helix-3 of the homeodomain of each family; and 2) PCR amplification of genomic DNA or cDNA, using two sets of degenerated oligonucleotides corresponding to helices 1 and 3 of the homeodomain or two specific domains of the POU family. Using the first strategy we have identified and characterized two tissue-specific cell determination and/or differentiation NK-type homeobox genes. Using the second strategy we have identified several homeobox genes that belong to the HOM/Hox, paired (prd) or POU families.     

Expression of a homologue of the fushi tarazu (ftz) gene in a cirripede crustacean

In Metazoa, Hox genes control the identity of the body parts along the anteroposterior axis. In addition to this homeotic function, these genes are characterized by two conserved features: They are clustered in the genome, and they contain a particular sequence, the homeobox, encoding a DNA-binding domain. Analysis of Hox homeobox sequences suggests that the Hox cluster emerged early in Metazoa and then underwent gene duplication events. In arthropods, the Hox cluster contains eight genes with a homeotic function and two other Hox-like genes, zerknullt (zen)/Hox3 and fushi tarazu (ftz). In insects, these two genes have lost their homeotic function but have acquired new functions in embryogenesis. In contrast, in chelicerates, these genes are expressed in a Hox-like pattern, which suggests that they have conserved their ancestral homeotic function. We describe here the characterization of Diva, the homologue of ftz in the cirripede crustacean Sacculina carcini. Diva is located in the Hox cluster, in the same position as the ftz genes of insects, and is not expressed in a Hox-like pattern. Instead, it is expressed exclusively in the central nervous system. Such a neurogenic expression of ftz has been also described in insects. This study, which provides the first information about the Hoxcluster in Crustacea, reveals that it may not be much smaller than the insect cluster. Study of the Diva expression pattern suggests that the arthropod ftz gene has lost its ancestral homeotic function after the divergence of the Crustacea/Hexapoda clade from other arthropod clades. In contrast, the function of ftz during neurogenesis is well conserved in insects and crustaceans.     

Small inverted repeats drive mitochondrial genome evolution in Lake Baikal sponges

Demosponges, the largest and most diverse class in the phylum Porifera, possess mitochondrial DNA (mtDNA) markedly different from that in other animals. Although several studies investigated evolution of demosponge mtDNA among major lineages of the group, the changes within these groups remain largely unexplored. Recently we determined mitochondrial genomic sequence of the Lake Baikal sponge Lubomirskia baicalensis and described proliferation of small inverted repeats (hairpins) that occurred in it since the divergence between L. baicalensis and the most closely related cosmopolitan freshwater sponge Ephydatia muelleri. Here we report mitochondrial genomes of three additional species of Lake Baikal sponges: Swartschewskia papyracea, Rezinkovia echinata and Baikalospongia intermedia morpha profundalis (Demospongiae, Haplosclerida, Lubomirskiidae) and from a more distantly related freshwater sponge Corvomeyenia sp. (Demospongiae, Haplosclerida, Metaniidae). We use these additional sequences to explore mtDNA evolution in Baikalian sponges, paying particular attention to the variation in the rates of nucleotide substitutions and the distribution of hairpins, abundant in these genomes. We show that most of the changes in Lubomirskiidae mitochondrial genomes are due to insertion/deletion/duplication of these elements rather than single nucleotide substitutions. Thus inverted repeats can act as an important force in evolution of mitochondrial genome architecture and be a valuable marker for population- and species-level studies in this group. In addition, we infer (((Rezinkovia+Lubomirskia)+Swartschewskia)+Baikalospongia) phylogeny for the family Lubomirskiidae based on the analysis of mitochondrial coding sequences from freshwater sponges.     

Hox and paraHox genes in bivalve molluscs

In this study, we sought the presence and analysed the sequences of the Hox and ParaHox genes in bivalve molluscs. The clustered Hox genes play a central role in anterior-posterior axial patterning in bilaterian metazoa, whereas the ParaHox gene cluster is a paralogue (evolutionary sister) of the Hox cluster.Using polymerase chain reaction (PCR)-based approaches, we isolated nine different sequences in five species belonging to three of the main bivalve subclasses: Ensis ensis and Tapes philippinarum (Heterodonta), Pecten maximus and Mytilus galloprovincialis (Pteriomorphia), and Yoldia eightsi (Protobranchia). Comparison with the Hox and ParaHox genes of other bilaterians, particularly lophotrochozoans, allowed us to attribute six of these sequences to the Hox gene cluster (one to paralog group [PG] 3 class, and five to the central class), two to the ParaHox cluster and one to the Gbx gene family.The results of our investigation seem to indicate that homeotic Hox and ParaHox gene clusters are homogeneous for both presence and characteristics in molluscs.     

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in <Emphasis Type="Italic">Solanum lycopersicum</Emphasis>

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A–E) and five subclades (A1–A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.     

Functional comparison of the Hoxa 4, Hoxa 10, and Hoxa 11 homeoboxes

A number of models attempt to explain the functional relationships of Hox genes. The functional equivalence model states that mammalian Hox-encoded proteins are largely functionally equivalent, and that Hox quantity is more important than Hox quality. In this report, we describe the results of two homeobox swaps. In one case, the homeobox of Hoxa 11 was replaced with that of the very closely related Hoxa 10. Developmental function was assayed by analyzing the phenotypes of all possible allele combinations, including the swapped allele, and null alleles for Hoxa 11 and Hoxd 11. This chimeric gene provided wild-type function in the development of the axial skeleton and male reproductive tract, but served as a hypomorph allele in the development of the appendicular skeleton, kidneys, and female reproductive tract. In the other case, the Hoxa 11 homeobox was replaced with that of the divergent Hoxa 4 gene. This chimeric gene provided near recessive null function in all tissues except the axial skeleton, which developed normally. These results demonstrate that even the most conserved regions of Hox genes, the homeoboxes, are not functionally interchangeable in the development of most tissues. In some cases, developmental function tracked with the homeobox, as previously seen in simpler organisms. Homeoboxes with more 5' cluster positions were generally dominant over more 3' homeoboxes, consistent with phenotypic suppression seen in Drosophila. Surprisingly, however, all Hox homeoboxes tested did appear functionally equivalent in the formation of the axial skeleton. The determination of segment identity is one of the most evolutionarily ancient functions of Hox genes. It is interesting that Hox homeoboxes are interchangeable in this process, but are functionally distinct in other aspects of development.     

Primary structure, developmentally regulated expression and potential duplication of the zebrafish homeobox gene ZF-21.

We report the molecular cloning and characterization of a cDNA derived from a zebrafish gene (ZF-21) related to the mouse homeobox containing gene Hox2.1. Interesting information about the differential conservation of various domains was gained from comparisons between the putative protein sequences from ZF-21 (275 amino acids) and Hox2.1 (279 aa). A separate DNA binding domain including the ZF-21 homeodomain and 36 additional flanking residues is completely identical to the C-terminal part of Hox2.1. As a consequence, these two mouse and zebrafish proteins must have identical DNA binding properties. A lower level of sequence identity between the N-terminal coding regions of ZF-21 and Hox2.1 reduces the total protein homology to 81%. However, short stretches of perfect homology in these N-terminals suggests that the essential biochemical functions are the same. As expected for true homologues, the ZF-21 and Hox2.1 genes also share extensive similarities with respect to non-coding sequences and temporal expression during embryogenesis. The finding of a potential ZF-21 duplication is discussed in relation to functional and evolutionary aspects of vertebrate homeobox genes.