Ambulacrarian prototypical Hox and ParaHox gene complements of the indirect-developing hemichordate <Emphasis Type="Italic">Balanoglossus simodensis</Emphasis>

The Hox genes and its evolutionary sister, the ParaHox genes, are widely distributed among animals. Although it has been expected that hemichordates and echinoderms have a single set of Hox genes and most likely a single set of ParaHox genes, it is not known whether the ortholog of Hox8 is absent in hemichordates, and in turn, consensus view about Hox/ParaHox gene complements in hemichordates has not been established. In this study, we isolated either complete or nearly complete coding sequences of 12 Hox genes, including the ortholog of the Hox8 that has not been reported in the previous studies, and three ParaHox genes from the recently discovered indirect-developing acorn worm, Balanoglossus simodensis. Our data suggest that the ancestral hemichordate had intact complements of ambulacrarian prototypical Hox and ParaHox genes, consisting of 12 and three members, respectively.     

Evolution of <Emphasis Type="Italic">Hox3</Emphasis> and <Emphasis Type="Italic">ftz</Emphasis> in arthropods: insights from the crustacean <Emphasis Type="Italic">Daphnia pulex</Emphasis>

The Drosophila melanogaster genes zerknüllt (zen) and fushi tarazu (ftz) are members of the Hox gene family whose roles have changed significantly in the insect lineage and thus provide an opportunity to study the mechanisms underlying the functional evolution of Hox proteins. We have studied the expression of orthologs of zen (DpuHox3) and ftz (Dpuftz) in the crustacean Daphnia pulex (Branchiopoda), both of which show a dynamic expression pattern. DpuHox3 is expressed in a complex pattern in early embryogenesis, with the most anterior boundary of expression lying at the anterior limit of the second antennal segment as well as a ring of expression around the embryo. In later embryos, DpuHox3 expression is restricted to the mesoderm of mandibular limb buds. Dpuftz is first expressed in a ring around the embryo following the posterior limit of the mandibular segment. Later, Dpuftz is restricted to the posterior part of the mandibular segment. This is the first report of expression of a Hox3 ortholog in a crustacean, and together with Dpuftz data, the results presented here show that Hox3 and ftz have retained a Hox-like expression pattern in crustaceans. This is in accordance with the proposed model of Hox3 and ftz evolution in arthropods and allows a more precise pinpointing of the loss of ftz “Hox-like behaviour”: in the lineage between the Branchiopoda and the basal insect Thysanura.     

The duplication of the <Emphasis Type="Italic">Hox</Emphasis> gene clusters in teleost fishes

Higher teleost fishes, including zebrafish and fugu, have duplicated their Hox genes relative to the gene inventory of other gnathostome lineages. The most widely accepted theory contends that the duplicate Hox clusters orginated synchronously during a single genome duplication event in the early history of ray-finned fishes. In this contribution we collect and re-evaluate all publicly available sequence information. In particular, we show that the short Hox gene fragments from published PCR surveys of the killifish Fundulus heteroclitus, the medaka Oryzias latipes and the goldfish Carassius auratus can be used to determine with little ambiguity not only their paralog group but also their membership in a particular cluster. Together with a survey of the genomic sequence data from the pufferfish Tetraodon nigroviridis we show that at least percomorpha, and possibly all eutelosts, share a system of 7 or 8 orthologous Hox gene clusters. There is little doubt about the orthology of the two teleost duplicates of the HoxA and HoxB clusters. A careful analysis of both the coding sequence of Hox genes and of conserved non-coding sequences provides additional support for the “duplication early” hypothesis that the Hox clusters in teleosts are derived from eight ancestral clusters by means of subsequent gene loss; the data remain ambiguous, however, in particular for the HoxC clusters. Assuming the “duplication early” hypothesis we use the new evidence on the Hox gene complements to determine the phylogenetic positions of gene-loss events in the wake of the cluster duplication. Surprisingly, we find that the resolution of redundancy seems to be a slow process that is still ongoing. A few suggestions on which additional sequence data would be most informative for resolving the history of the teleostean Hox genes are discussed. Supplemental material is available at http://www.bioinf.uni-leipzig.de/Publications/SUPPLEMENTS/04-006/.     

The <Emphasis Type="Italic">C. elegans Hox</Emphasis> gene <Emphasis Type="Italic">ceh-13</Emphasis> regulates cell migration and fusion in a non-colinear way. Implications for the early evolution of <Emphasis Type="Italic">Hox</Emphasis>clusters

Background  

Hox genes play a central role in axial patterning during animal development. They are clustered in the genome and specify cell fate in sequential domains along the anteroposterior (A-P) body axis in a conserved order that is co-linear with their relative genomic position. In the soil worm Caenorhabditis elegans, this striking rule of co-linearity is broken by the anterior Hox gene ceh-13, which is located between the two middle Hox paralogs, lin-39 and mab-5, within the loosely organized nematode Hox cluster. Despite its evolutionary and developmental significance, the functional consequence of this unusual genomic organization remains unresolved.     

Isolation and expression analysis of a poriferan <Emphasis Type="Italic">Antp</Emphasis>-class <Emphasis Type="Italic">Bar-/Bsh-</Emphasis>like homeobox gene

A novel non-Hox Antp-class gene (BarBsh-Hb) was isolated from the marine sponge Halichondria sp. This gene shares high sequence identity with eumetazoan genes from the Bsh and Bar gene families and can be distinguished from other non-Hox Antp-class genes by diagnostic residues. We also present an alignment of all known (full-length) poriferan non-Hox Antp-class genes. Maximum likelihood methods were employed to estimate phylogenetic relationships among non-Hox genes and BarBsh-Hb. We employed RT-PCR techniques to look at expression across different developmental stages (larval to rhagon). BarBsh-Hb product was present in newly released larvae, but expression was not detected 8–16 h post-release. Expression of BarBsh-Hb was detected in later-stage (>16 h post-release), free-swimming larvae until they settled and attached to the substratum, after which expression was down-regulated. In a separate set of experiments, low levels of expression were observed in normal adult tissue and disaggregated adult tissue, but BarBsh-Hb expression increased during tissue re-aggregation. These data increase the number of non-Hox homeobox genes identified in sponges and provide evidence of regulation of this non-Hox gene during sponge development. While the Bar and Bsh genes play important roles in the development of nervous tissue—especially visual systems—in metazoans, the specific role(s) BarBsh-Hb play(s) in sponge development is unclear and deserves greater attention.Edited by C. Desplan     

The<Emphasis Type="Italic"> Trox-2</Emphasis> Hox/ParaHox gene of<Emphasis Type="Italic"> Trichoplax</Emphasis> (Placozoa) marks an epithelial boundary

Hox and ParaHox genes are implicated in axial patterning of cnidarians and bilaterians, and are thought to have originated by tandem duplication of a single ProtoHox gene followed by duplication of the resultant gene cluster. It is unclear what the ancestral role of Hox/ParaHox genes was before the divergence of Cnidaria and Bilateria, or what roles the postulated ProtoHox gene(s) played. Here we describe the full coding region, spatial expression and function of Trox-2, the single Hox/ParaHox-type gene identified in Trichoplax adhaerens (phylum Placozoa) and either a candidate ProtoHox or a ParaHox gene. Trox-2 is expressed in a ring around the periphery of Trichoplax, in small cells located between the outer margins of the upper and lower epithelial cell layers. Inhibition of Trox-2 function, either by uptake of morpholino antisense oligonucleotides or by RNA interference, causes complete cessation of growth and binary fission. We speculate that Trox-2 functions within a hitherto unrecognized population of possibly multipotential peripheral stem cells that contribute to differentiated cells at the epithelial boundary of Trichoplax.Edited by D. Tautz     

<Emphasis Type="Italic">Hox</Emphasis> cluster duplication in the basal teleost <Emphasis Type="Italic">Hiodon alosoides</Emphasis> (Osteoglossomorpha)

Large-scale—even genome-wide—duplications have repeatedly been invoked as an explanation for major radiations. Teleosts, the most species-rich vertebrate clade, underwent a “fish-specific genome duplication” (FSGD) that is shared by most ray-finned fish lineages. We investigate here the Hox complement of the goldeye (Hiodon alosoides), a representative of Osteoglossomorpha, the most basal teleostean clade. An extensive PCR survey reveals that goldeye has at least eight Hox clusters, indicating a duplicated genome compared to basal actinopterygians. The possession of duplicated Hox clusters is uncoupled to species richness. The Hox system of the goldeye is substantially different from that of other teleost lineages, having retained several duplicates of Hox genes for which crown teleosts have lost at least one copy. A detailed analysis of the PCR fragments as well as full length sequences of two HoxA13 paralogs, and HoxA10 and HoxC4 genes places the duplication event close in time to the divergence of Osteoglossomorpha and crown teleosts. The data are consistent with—but do not conclusively prove—that Osteoglossomorpha shares the FSGD. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

A PCR survey of <Emphasis Type="Italic">Hox</Emphasis> genes in the myzostomid <Emphasis Type="Italic">Myzostoma cirriferum</Emphasis>

Using degenerate primers, we were able to identify seven Hox genes for the myzostomid Myzostoma cirriferum. The recovered fragments belong to anterior class (Mci_lab, Mci_pb), central class (Mci_Dfd, Mci_Lox5, Mci_Antp, Mci_Lox4), and posterior class (Mci_Post2) paralog groups. Orthology assignment was verified by phylogenetic analyses and presence of diagnostic regions in the homeodomain as well as flanking regions. The presence of Lox5, Lox4, and Post2 supports the inclusion of Myzostomida within Lophotrochozoa. We found signature residues within flanking regions of Lox5, which are also found in annelids, but not in Platyhelminthes. As such the available Hox genes data of myzostomids support an annelid relationship. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genesis and evolution of the <Emphasis Type="Italic">Evx</Emphasis> and <Emphasis Type="Italic">Mox</Emphasis>genes and the extended Hox and ParaHox gene clusters

Background

Hox and ParaHox gene clusters are thought to have resulted from the duplication of a ProtoHox gene cluster early in metazoan evolution. However, the origin and evolution of the other genes belonging to the extended Hox group of homeobox-containing genes, that is, Mox and Evx, remains obscure. We constructed phylogenetic trees with mouse, amphioxus and Drosophila extended Hox and other related Antennapedia-type homeobox gene sequences and analyzed the linkage data available for such genes.

Results

We claim that neither Mox nor Evx is a Hox or ParaHox gene. We propose a scenario that reconciles phylogeny with linkage data, in which an Evx/Mox ancestor gene linked to a ProtoHox cluster was involved in a segmental tandem duplication event that generated an array of all Hox-like genes, referred to as the 'coupled' cluster. A chromosomal breakage within this cluster explains the current composition of the extended Hox cluster (with Evx, Hox and Mox genes) and the ParaHox cluster.

Conclusions

Most studies dealing with the origin and evolution of Hox and ParaHox clusters have not included the Hox-related genes Mox and Evx. Our phylogenetic analyses and the available linkage data in mammalian genomes support an evolutionary scenario in which an ancestor of Evx and Mox was linked to the ProtoHox cluster, and that a tandem duplication of a large genomic region early in metazoan evolution generated the Hox and ParaHox clusters, plus the cluster-neighbors Evx and Mox. The large 'coupled' Hox-like cluster EvxHox/MoxParaHox was subsequently broken, thus grouping the Mox and Evx genes to the Hox clusters, and isolating the ParaHox cluster.

     

Isolation and sequence analysis of <Emphasis Type="Italic">Sox</Emphasis> genes from lizard <Emphasis Type="Italic">Eremias</Emphasis> multiocellata

The Sox (SRY-related high-mobility-group box) family of genes shares a conserved HMG box and is involved in a diverse range of developmental processes and sex determination in vertebrates. Twenty Sox genes are present in the genomes of humans and mice, but far less is known about the Sox gene family in reptiles. Using two pairs of highly degenerate primers designed from a multiple alignment of Sox amino acid sequences in several species, different positive clones were obtained from male and female Eremias multiocellata, a viviparous lizard which is subject to TSD (temperature-dependent sex determination). These clones were sequenced and identified. They are members of the SoxB (Sox2, Sox14), SoxC (Sox11, Sox12) and SoxE (Sox9a, Sox9b, Sox10) groups. No sex-specific differences were observed. Based on the amino acid sequence similarities, the phylogenetic analysis was carried out and these genes clustered with their orthologues. In addition, we found the gene duplication in E. multiocellata, it may be a mechanism to produce new functional genes.     

<Emphasis Type="Italic">Hox</Emphasis> gene survey in the chaetognath<Emphasis Type="Italic"> Spadella cephaloptera</Emphasis>: evolutionary implications

We present the isolation of six Hox genes in the chaetognath Spadella cephaloptera. We identified one member of the paralogy group 3, four median genes and a mosaic gene that shares features of both median and posterior classes ( SceMedPost). Several hypotheses may account for the presence of a mosaic Hox gene in this animal. Here we propose that SceMedPost may represent an ancestral gene, which has not diverged totally into a posterior or a median one. This hypothesis has interesting implications for the reconstruction of the evolutionary history of Hox genes and suggests that Chaetognatha lineage divergence could predate the deuterostome/protostome split. Such a phylogenetic position is considered in the light of their embryological and morphological characters.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

Molecular analysis of the genus <Emphasis Type="Italic">Anoxybacillus</Emphasis> based on sequence similarity of the genes <Emphasis Type="Italic">recN</Emphasis>, <Emphasis Type="Italic">flaA</Emphasis>, and <Emphasis Type="Italic">ftsY</Emphasis>

Genome predictions based on selected genes would be a very welcome approach for taxonomic studies. We analyzed three genes, recN, flaA, and ftsY, for determining if these genes are useful tools for systematic analyses in the genus Anoxybacillus. The genes encoding a DNA repair and genetic recombination protein (recN), the flagellin protein (flaA), and GTPase signal docking protein (ftsY) were sequenced for ten Anoxybacillus species. The sequence comparisons revealed that recN sequence similarities range between 61% and 99% in the genus Anoxybacillus. Comparisons to other bacterial recN genes indicated that levels of similarity did not differ from the levels within genus Anoxybacillus. These data showed that recN is not a useful marker for the genus Anoxybacillus. A 550–600-bp region of the flagellin gene was amplified for all Anoxybacillus strains except for Anoxybacillus contaminans. The sequence similarity of flaA gene varies between 61% and 76%. Comparisons to other bacterial flagellin genes obtained from GenBank (Bacillus, Pectinatus, Proteus, and Vibrio) indicated that the levels of similarity were lower (3–42%). Based on these data, we concluded that the variability in this single gene makes it a particularly useful marker. Another housekeeping gene ftsY suggested to reflect the G+C (mol/mol) content of whole genome was analyzed for Anoxybacillus strains. A mean difference of 1.4% was observed between the G+C content of the gene ftsY and the G+C content of the whole genome. These results showed that the gene ftsY can be used to represent whole G+C content of the Anoxybacillus species.     

Micro-colinearity between rice, <Emphasis Type="Italic">Brachypodium</Emphasis>, and <Emphasis Type="Italic">Triticum monococcum</Emphasis> at the wheat domestication locus <Emphasis Type="Italic">Q</Emphasis>

Brachypodium, a wild temperate grass with a small genome, was recently proposed as a new model organism for the large-genome grasses. In this study, we evaluated gene content and microcolinearity between diploid wheat (Triticum monococcum), Brachypodium sylvaticum, and rice at a local genomic region harboring the major wheat domestication gene Q. Gene density was much lower in T. monococcum (one per 41 kb) because of gene duplication and an abundance of transposable elements within intergenic regions as compared to B. sylvaticum (one per 14 kb) and rice (one per 10 kb). For the Q gene region, microcolinearity was more conserved between wheat and rice than between wheat and Brachypodium because B. sylvaticum contained two genes apparently not present within the orthologous regions of T. monococcum and rice. However, phylogenetic analysis of Q and leukotriene A-4 hydrolase-like gene orthologs, which were colinear among the three species, showed that Brachypodium is more closely related to wheat than rice, which agrees with previous studies. We conclude that Brachypodium will be a useful tool for gene discovery, comparative genomics, and the study of evolutionary relationships among the grasses but will not preclude the need to conduct large-scale genomics experiments in the Triticeae.     

Sequence analysis of the grey mouse lemur (<Emphasis Type="Italic">Microcebus murinus</Emphasis>) <Emphasis Type="Italic">MHC class II DQ</Emphasis> and <Emphasis Type="Italic">DR</Emphasis> region

We here report the genomic organisation of the grey mouse lemur (Microcebus murinus) MHC class II DQ and DR region based on BAC clone analysis. The sequenced Mimu-MHC haplotype spans 343 kb and encompasses the genes TAP2, DOB, DQB, DQA, DRB, DRA, BTNL2 and a further BTNL gene. The DQ and DR genes of this haplotype are not duplicated. Mimu-DOB is not transcribed and represents a pseudogene due to deletions and premature stop codons. Analysis of BAC clone DNA, a cDNA sample and eight genomic DNA samples suggests that Mimu-DRB, Mimu-DQA and Mimu-DQB are highly polymorphic with the majority of peptide-binding residues being affected by polymorphisms. In contrast, Mimu-DRA is moderately polymorphic, and the variable amino acid positions are not part of the peptide-binding region. Phylogenetic analysis of Mimu-DQA and Mimu-DQB and other primate DQA and DQB genes indicates that duplication of DQA and DQB loci occurred in Anthropoidea after the split from Strepsirrhini.     

Type III Secretion Homologs Are Present in <Emphasis Type="Italic">Brucella melitensis</Emphasis>, <Emphasis Type="Italic">B. ovis</Emphasis>, and <Emphasis Type="Italic">B. suis</Emphasis> biovars 1, 2, and 3

Protein sequences from characterized type III secretion (TTS) systems were used as probes in silico to identify several TTS gene homologs in the genome sequence of Brucella suis biovar 1 strain 1330. Four of the genes, named flhB, fliP, fliR, and fliF on the basis of greatest homologies to known flagellar apparatus proteins, were targeted in PCR and hybridization assays to determine their distribution among other Brucella nomen species and biovars. The results indicated that flhB, fliP, fliR and fliF are present in Brucella melitensis, Brucella ovis, and Brucella suis biovars 1, 2 and 3. Similar homologos have been reported previously in Brucella abortus. Using RT-PCR assays, we were unable to detect any expression of these genes. It is not yet known whether the genes are the cryptic remnants of a flagellar system or are actively involved in a process contributing to pathogenicity or previously undetected motility, but they are distributed widely in Brucella and merit further study to determine their role. Received: 11 February 2002 / Accepted: 13 June 2002     

Distribution of <Emphasis Type="Italic">Hordoindoline</Emphasis> genes in the genus <Emphasis Type="Italic">Hordeum</Emphasis>

Hordoindoline (Hin) genes, which are known to comprise Hina, Hinb-1, and Hinb-2, are associated with grain hardness in barley. However, the interspecific variation in the Hin genes in the genus Hordeum has not been studied in detail. We examined the variation in Hin genes and used it to infer the phylogenetic relationships between the genes found in two H. vulgare subspecies (cultivated barley and H. vulgare subsp. spontaneum) and 10 wild relatives (H. bogdanii, H. brachyantherum, H. bulbosum, H. chilense, H. comosum, H. marinum, H. murinum, H. patagonicum, H. pusillum, and H. roshevitzii). The Hina and Hinb genes of these species were amplified by PCR. We found two Hinb genes in three wild species (H. bogdanii, H. brachyantherum, and H. roshevitzii) and preliminarily named them Hinb-A and Hinb-B. Cluster analysis showed that the 17 Hinb genes present in Hordeum formed two distinct clusters (named A and B). Seven Hinb genes were included in Cluster-A, and 10 Hinb genes were included in Cluster-B. All Hinb-A genes were included in Cluster-A, while all of the Hinb-B genes were included in Cluster-B. In contrast, the Hinb-1 and Hinb-2 genes in H. vulgare were included in Cluster-B. These results suggest that the Hinb genes duplicated during the early stages of diversification in the genus Hordeum. On the other hand, the Hinb-1 and Hinb-2 genes in H. vulgare seem to have been generated by a duplication of the Hinb gene after the split of the lineages leading to H. vulgare and H. bulbosum.     

Characterization of Maltase Clusters in the Genus <Emphasis Type="Italic">Drosophila</Emphasis>

To reveal evolutionary history of maltase gene family in the genus Drosophila, we undertook a bioinformatics study of maltase genes from available genomes of 12 Drosophila species. Molecular evolution of a closely related glycoside hydrolase, the α-amylase, in Drosophila has been extensively studied for a long time. The α-amylases were even used as a model of evolution of multigene families. On the other hand, maltase, i.e., the α-glucosidase, got only scarce attention. In this study, we, therefore, investigated spatial organization of the maltase genes in Drosophila genomes, compared the amino acid sequences of the encoded enzymes and analyzed the intron/exon composition of orthologous genes. We found that the Drosophila maltases are more numerous than previously thought (ten instead of three genes) and are localized in two clusters on two chromosomes (2L and 2R). To elucidate the approximate time line of evolution of the clusters, we estimated the order and dated duplication of all the 10 genes. Both clusters are the result of ancient series of subsequent duplication events, which took place from 352 to 61 million years ago, i.e., well before speciation to extant Drosophila species. Also observed was a remarkable intron/exon composition diversity of particular maltase genes of these clusters, probably a result of independent intron loss after duplication of intron-rich gene ancestor, which emerged well before speciation in a common ancestor of all extant Drosophila species.