Breakup of a homeobox cluster after genome duplication in teleosts

Spatial Colinear but Broken Temporal Expression of Duplicated ParaHox Genes in Asexually Reproducing Annelids, Nais communis and Pristina longiseta

ParaHox genes are key developmental regulators involved in the patterning of the digestive tract along the anteroposterior axis and the development of the nervous system. Most studies have focused on the function of these genes in embryogenesis, while their expression patterns in postembryonic development often remain unknown. In this study, we identified for the first time all ParaHox orthologs in two naidid oligochaetes, N. communis and P. longiseta, and described their expression patterns during normal growth and fission in these animals. We showed that Gsx and Cdx are presented by two paralogs, while Xlox is a single copy gene in both species. Using whole-mount in situ hybridization, we also found that orthologs, except for the Xlox gene, have similar activity patterns with minor differences in details, while the expression patterns of paralogs can differ significantly. However, all these genes are involved in axial patterning and/or in tissue remodeling during growth and asexual reproduction in naidids. Moreover, during paratomic fission, these genes are expressed with spatial colinearity but temporal colinearity is broken. The results of this study may be evidence of the functional diversification of duplicated genes and suggest involvement of the ParaHox genes in whole-body patterning during growth and asexual reproduction in annelids.

Genome of the four-finger threadfin Eleutheronema tetradactylum (Perciforms: Polynemidae)

Background

Teleost fish play important roles in aquatic ecosystems and aquaculture. Threadfins (Perciformes: Polynemidae) show a range of interesting biology, and are of considerable importance for both wild fisheries and aquaculture. Additionally, the four-finger threadfin Eleutheronema tetradactylum is of conservation relevance since its populations are considered to be in rapid decline and it is classified as endangered. However, no genomic resources are currently available for the threadfin family Polynemidae.

Results

We sequenced and assembled the first threadfin fish genome, the four-finger threadfin E. tetradactylum. We provide a genome assembly for E. tetradactylum with high contiguity (scaffold N50 = 56.3 kb) and high BUSCO completeness at 96.5%. The assembled genome size of E. tetradactylum is just 610.5 Mb, making it the second smallest perciform genome assembled to date. Just 9.07–10.91% of the genome sequence was found to consist of repetitive elements (standard RepeatMasker analysis vs custom analysis), making this the lowest repeat content identified to date for any perciform fish. A total of 37,683 protein-coding genes were annotated, and we include analyses of developmental transcription factors, including the Hox, ParaHox, and Sox families. MicroRNA genes were also annotated and compared with other chordate lineages, elucidating the gains and losses of chordate microRNAs.

Conclusions

The four-finger threadfin E. tetradactylum genome presented here represents the first available genome sequence for the ecologically, biologically, and commercially important clade of threadfin fish. Our findings provide a useful genomic resource for future research into the interesting biology and evolution of this valuable group of food fish.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07145-1.

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.

Mutation of amphioxus Pdx and Cdx demonstrates conserved roles for ParaHox genes in gut, anus and tail patterning

BACKGROUND

The homeobox genes Pdx and Cdx are widespread across the animal kingdom and part of the small ParaHox gene cluster. Gene expression patterns suggest ancient roles for Pdx and Cdx in patterning the through-gut of bilaterian animals although functional data are available for few lineages. To examine evolutionary conservation of Pdx and Cdx gene functions, we focus on amphioxus, small marine animals that occupy a pivotal position in chordate evolution and in which ParaHox gene clustering was first reported.

RESULTS

Using transcription activator-like effector nucleases (TALENs), we engineer frameshift mutations in the Pdx and Cdx genes of the amphioxus Branchiostoma floridae and establish mutant lines. Homozygous Pdx mutants have a defect in amphioxus endoderm, manifest as loss of a midgut region expressing endogenous GFP. The anus fails to open in homozygous Cdx mutants, which also have defects in posterior body extension and epidermal tail fin development. Treatment with an inverse agonist of retinoic acid (RA) signalling partially rescues the axial and tail fin phenotypes indicating they are caused by increased RA signalling. Gene expression analyses and luciferase assays suggest that posterior RA levels are kept low in wild type animals by a likely direct transcriptional regulation of a Cyp26 gene by Cdx. Transcriptome analysis reveals extensive gene expression changes in mutants, with a disproportionate effect of Pdx and Cdx on gut-enriched genes and a colinear-like effect of Cdx on Hox genes.

CONCLUSIONS

These data reveal that amphioxus Pdx and Cdx have roles in specifying middle and posterior cell fates in the endoderm of the gut, roles that likely date to the origin of Bilateria. This conclusion is consistent with these two ParaHox genes playing a role in the origin of the bilaterian through-gut with a distinct anus, morphological innovations that contributed to ecological change in the Cambrian. In addition, we find that amphioxus Cdx promotes body axis extension through a molecular mechanism conserved with vertebrates. The axial extension role for Cdx dates back at least to the origin of Chordata and may have facilitated the evolution of the post-anal tail and active locomotion in chordates.

Comparative studies on duplicated tdrd7 paralogs in teleosts: Molecular evolution caused neo-functionalization

The third-round whole genome duplication (3R-WGD) event occurred in the stem lineage of teleost during evolution, and is considered to be responsible for the biological diversification of ray-finned fishes. TUDOR domain containing protein 7 (Tdrd7), which belongs to the Tudor family proteins has been widely discussed in mammals. However, information about this gene in teleost is still lacking. In this study, two teleost tdrd7 genes (tdrd7a and tdrd7b) were identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). Through genomic structure, phylogenetic, synteny analysis and online bioinformatic mining of tdrd7 duplications in other selected species, we confirmed that tdrd7a/7b were originated from the teleost-specific 3R-WGD. The tdrd7a is specific to teleost except for spotted gar. The tdrd7a showed a higher molecular evolution rate than tdrd7b with longer branch-length in the phylogenetic tree and multiple positively selected sites. Interestingly, it showed gonad specific expression pattern in adult tissues and germ cell specific distribution in embryos and gonads. Its 3'-untranslated region (3'UTR) labeled eGFP/DsRED could visualize primordial germ cells (PGCs) in zebrafish embryos. The tdrd7b did not show similar tissue and cell type specificity. These characteristic differences between the duplicated tdrd7 paralogues suggest that tdrd7a and tdrd7b have undergone neofunctionalization in Japanese flounder. Our results provide novel insight into the evolution and functional diversification of teleost tdrd7 genes deserving further investigations.

A TALE of shrimps: Genome-wide survey of homeobox genes in 120 species from diverse crustacean taxa

The homeodomain-containing proteins are an important group of transcription factors found in most eukaryotes including animals, plants and fungi. Homeobox genes are responsible for a wide range of critical developmental and physiological processes, ranging from embryonic development, innate immune homeostasis to whole-body regeneration. With continued fascination on this key class of proteins by developmental and evolutionary biologists, multiple efforts have thus far focused on the identification and characterization of homeobox orthologs from key model organisms in attempts to infer their evolutionary origin and how this underpins the evolution of complex body plans. Despite their importance, the genetic complement of homeobox genes has yet been described in one of the most valuable groups of animals representing economically important food crops. With crustacean aquaculture being a growing industry worldwide, it is clear that systematic and cross-species identification of crustacean homeobox orthologs is necessary in order to harness this genetic circuitry for the improvement of aquaculture sustainability. Using publicly available transcriptome data sets, we identified a total of 4183 putative homeobox genes from 120 crustacean species that include food crop species, such as lobsters, shrimps, crayfish and crabs. Additionally, we identified 717 homeobox orthologs from 6 other non-crustacean arthropods, which include the scorpion, deer tick, mosquitoes and centipede. This high confidence set of homeobox genes will now serve as a key resource to the broader community for future functional and comparative genomics studies.

Lampreys, the jawless vertebrates, contain only two ParaHox gene clusters

ParaHox genes (Gsx, Pdx, and Cdx) are an ancient family of developmental genes closely related to the Hox genes. They play critical roles in the patterning of brain and gut. The basal chordate, amphioxus, contains a single ParaHox cluster comprising one member of each family, whereas nonteleost jawed vertebrates contain four ParaHox genomic loci with six or seven ParaHox genes. Teleosts, which have experienced an additional whole-genome duplication, contain six ParaHox genomic loci with six ParaHox genes. Jawless vertebrates, represented by lampreys and hagfish, are the most ancient group of vertebrates and are crucial for understanding the origin and evolution of vertebrate gene families. We have previously shown that lampreys contain six Hox gene loci. Here we report that lampreys contain only two ParaHox gene clusters (designated as α- and β-clusters) bearing five ParaHox genes (Gsxα, Pdxα, Cdxα, Gsxβ, and Cdxβ). The order and orientation of the three genes in the α-cluster are identical to that of the single cluster in amphioxus. However, the orientation of Gsxβ in the β-cluster is inverted. Interestingly, Gsxβ is expressed in the eye, unlike its homologs in jawed vertebrates, which are expressed mainly in the brain. The lamprey Pdxα is expressed in the pancreas similar to jawed vertebrate Pdx genes, indicating that the pancreatic expression of Pdx was acquired before the divergence of jawless and jawed vertebrate lineages. It is likely that the lamprey Pdxα plays a crucial role in pancreas specification and insulin production similar to the Pdx of jawed vertebrates.

"Holostei versus Halecostomi" Problem: Insight from Cytogenetics of Ancient Nonteleost Actinopterygian Fish, Bowfin Amia calva

Bowfin belongs to an ancient lineage of nonteleost ray-finned fishes (actinopterygians) and is the only extant survivor of a once diverged group, the Halecomorphi or Amiiformes. Owing to the scarcity of extant nonteleost ray-finned lineages, also referred as "living fossils," their phylogenetic interrelationships have been the target of multiple hypotheses concerning their sister group relationships. Molecular and morphological data sets have produced controversial results; bowfin is considered as either the sister group to genome-duplicated teleosts (together forming the group of Halecostomi) or to gars (Lepisosteiformes; together forming the group of Holostei). However, any detailed cytogenetic analysis of bowfin chromosomes has never been performed to address this issue. Here we examined bowfin chromosomes by conventional (Giemsa-staining, C-banding, base-specific fluorescence and silver staining) and molecular (FISH with rDNA probes) cytogenetic protocols. We identified diploid chromosome number 2n = 46 with a middle-sized submetacentric chromosome pair as the major ribosomal DNA-bearing (45S rDNA), GC-positive and silver-positive element. The minor rDNA (5S rDNA) sites were localized in the pericentromeric region of one middle-sized acrocentric chromosome pair. Comparison with available cytogenetic data of other nonteleost actinopterygians (bichirs, sturgeons, gars) and teleost species including representative of basally branching lineages showed bowfin chromosomal characteristics more similar to the teleost type than to any other nonteleosts. Particularly striking differences were identified between bowfin and gars, the latter of which were found to mimic mammalian AT/GC genomic organisation. Such conclusion however contradicts the most recent phylogenomic results and raises the question what states are ancestral and what are derived.

Evolution history of duplicated smad3 genes in teleost: insights from Japanese flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication (WGD) during deuterosome evolution, a third genome duplication occurred in the ray-fined fish lineage and is considered to be responsible for the teleost-specific lineage diversification and regulation mechanisms. As a receptor-regulated SMAD (R-SMAD), the function of SMAD3 was widely studied in mammals. However, limited information of its role or putative paralogs is available in ray-finned fishes. In this study, two SMAD3 paralogs were first identified in the transcriptome and genome of Japanese flounder (Paralichthys olivaceus). We also explored SMAD3 duplication in other selected species. Following identification, genomic structure, phylogenetic reconstruction, and synteny analyses performed by MrBayes and online bioinformatic tools confirmed that smad3a/3b most likely originated from the teleost-specific WGD. Additionally, selection pressure analysis and expression pattern of the two genes performed by PAML and quantitative real-time PCR (qRT-PCR) revealed evidence of subfunctionalization of the two SMAD3 paralogs in teleost. Our results indicate that two SMAD3 genes originate from teleost-specific WGD, remain transcriptionally active, and may have likely undergone subfunctionalization. This study provides novel insights to the evolution fates of smad3a/3b and draws attentions to future function analysis of SMAD3 gene family.

Origin and evolution of GATA2a and GATA2b in teleosts: insights from tongue sole, Cynoglossus semilaevis

Background. Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication occurred in the lineage of teleost fish and is considered to be responsible for much of the biological diversification within the lineage. GATA2, a member of GATA family of transcription factors, is an important regulator of gene expression in hematopoietic cell in mammals, yet the role of this gene or its putative paralogs in ray-finned fishes remains relatively unknown.

Methods. In this study, we attempted to identify GATA2 sequences from the transcriptomes and genomes of multiple teleosts using the bioinformatic tools MrBayes, MEME, and PAML. Following identification, comparative analysis of genome structure, molecular evolution rate, and expression by real-time qPCR were used to predict functional divergence of GATA2 paralogs and their relative transcription in organs of female and male tongue soles (Cynoglossus semilaevis).

Results. Two teleost GATA2 genes were identified in the transcriptomes of tongue sole and Japanese flounder (Paralichthysolivaceus). Synteny and phylogenetic analysis confirmed that the two genes likely originated from the teleost-specific genome duplication . Additionally, selection pressure analysis predicted these gene duplicates to have undergone purifying selection and possible divergent new functions. This was supported by differential expression pattern of GATA2a and GATA2b observed in organs of female and male tongue soles.

Discussion. Our results indicate that two GATA2 genes originating from the first teleost-specific genome duplication have remained transcriptionally active in some fish species and have likely undergone neofunctionalization. This knowledge provides novel insights into the evolution of the teleost GATA2 genes and constituted important groundwork for further research on the GATA gene family.

TCF/Lef regulates the Gsx ParaHox gene in central nervous system development in chordates

BACKGROUND

The ParaHox genes play an integral role in the anterior-posterior (A-P) patterning of the nervous system and gut of most animals. The ParaHox cluster is an ideal system in which to study the evolution and regulation of developmental genes and gene clusters, as it displays similar regulatory phenomena to its sister cluster, the Hox cluster, but offers a much simpler system with only three genes.

RESULTS

Using Ciona intestinalis transgenics, we isolated a regulatory element upstream of Branchiostoma floridae Gsx that drives expression within the central nervous system of Ciona embryos. The minimal amphioxus enhancer region required to drive CNS expression has been identified, along with surrounding sequence that increases the efficiency of reporter expression throughout the Ciona CNS. TCF/Lef binding sites were identified and mutagenized and found to be required to drive the CNS expression. Also, individual contributions of TCF/Lef sites varied across the regulatory region, revealing a partial division of function across the Bf-Gsx-Up regulatory element. Finally, when all TCF/Lef binding sites are mutated CNS expression is not only abolished, but a latent repressive function is also unmasked.

CONCLUSIONS

We have identified a B. floridae Gsx upstream regulatory element that drives CNS expression within transgenic Ciona intestinalis, and have shown that this CNS expression is dependent upon TCF/Lef binding sites. We examine the evolutionary and developmental implications of these results, and discuss the possibility of TCF/Lef not only as a regulator of chordate Gsx, but as a deeply conserved regulatory factor controlling all three ParaHox genes across the Metazoa.

Molecular characterization, origin, and evolution of teleost p68 gene family: Insights from Japanese flounder, Paralichthys olivaceus

Two rounds of whole-genome duplication occurred in the common ancestor of vertebrates. Later, a third round genome duplication occurred in the teleost fishes. As a prototype member of DEAD-box RNA helicases, the function of p68 helicase in development has been well investigated in human, however, limited information is available regarding the regulatory function of this gene in the development of teleosts. In this study, being an important farmed fish in North China, Japanese flounder (Paralichthys olivaceus) was used as model fish to investigate the role of p68 gene in teleost development. Two p68 genes were first identified from Japanese flounder. Molecular characterization of them was performed by analyzing the exon-intron boundaries. Then, we confirmed that such two teleost p68 genes originated from teleost-specific genome duplication through phylogenetic and synteny analyses. Additionally, comparative analyses of amino acid sequences, variation in selective pressure, and expression profiles of p68 genes revealed probable sub-functionalization fate of teleost p68 genes after the duplication. Therefore, this study supplements the evolutionary properties of teleost p68 gene family and provides the groundwork for further studying the regulatory function of p68 genes in the development of teleosts.

Cdx ParaHox genes acquired distinct developmental roles after gene duplication in vertebrate evolution

Background

The functional consequences of whole genome duplications in vertebrate evolution are not fully understood. It remains unclear, for instance, why paralogues were retained in some gene families but extensively lost in others. Cdx homeobox genes encode conserved transcription factors controlling posterior development across diverse bilaterians. These genes are part of the ParaHox gene cluster. Multiple Cdx copies were retained after genome duplication, raising questions about how functional divergence, overlap, and redundancy respectively contributed to their retention and evolutionary fate.

Results

We examined the degree of regulatory and functional overlap between the three vertebrate Cdx genes using single and triple morpholino knock-down in Xenopus tropicalis followed by RNA-seq. We found that one paralogue, Cdx4, has a much stronger effect on gene expression than the others, including a strong regulatory effect on FGF and Wnt genes. Functional annotation revealed distinct and overlapping roles and subtly different temporal windows of action for each gene. The data also reveal a colinear-like effect of Cdx genes on Hox genes, with repression of Hox paralogy groups 1 and 2, and activation increasing from Hox group 5 to 11. We also highlight cases in which duplicated genes regulate distinct paralogous targets revealing pathway elaboration after whole genome duplication.

Conclusions

Despite shared core pathways, Cdx paralogues have acquired distinct regulatory roles during development. This implies that the degree of functional overlap between paralogues is relatively low and that gene expression pattern alone should be used with caution when investigating the functional evolution of duplicated genes. We therefore suggest that developmental programmes were extensively rewired after whole genome duplication in the early evolution of vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0165-x) contains supplementary material, which is available to authorized users.

Comparative Evolution of Duplicated Ddx3 Genes in Teleosts: Insights from Japanese Flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication event occurred in the stem lineage of ray-finned fishes. This teleost-specific genome duplication is thought to be responsible for the biological diversification of ray-finned fishes. DEAD-box polypeptide 3 (DDX3) belongs to the DEAD-box RNA helicase family. Although their functions in humans have been well studied, limited information is available regarding their function in teleosts. In this study, two teleost Ddx3 genes were first identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). We confirmed that the two genes originated from teleost-specific genome duplication through synteny and phylogenetic analysis. Additionally, comparative analysis of genome structure, molecular evolution rate, and expression pattern of the two genes in Japanese flounder revealed evidence of subfunctionalization of the duplicated Ddx3 genes in teleosts. Thus, the results of this study reveal novel insights into the evolution of the teleost Ddx3 genes and constitute important groundwork for further research on this gene family.

Did homeobox gene duplications contribute to the Cambrian explosion?

The Cambrian explosion describes an apparently rapid increase in the diversity of bilaterian animals around 540-515 million years ago. Bilaterian animals explore the world in three-dimensions deploying forward-facing sense organs, a brain, and an anterior mouth; they possess muscle blocks enabling efficient crawling and burrowing in sediments, and they typically have an efficient 'through-gut' with separate mouth and anus to process bulk food and eject waste, even when burrowing in sediment. A variety of ecological, environmental, genetic, and developmental factors have been proposed as possible triggers and correlates of the Cambrian explosion, and it is likely that a combination of factors were involved. Here, I focus on a set of developmental genetic changes and propose these are part of the mix of permissive factors. I describe how ANTP-class homeobox genes, which encode transcription factors involved in body patterning, increased in number in the bilaterian stem lineage and earlier. These gene duplications generated a large array of ANTP class genes, including three distinct gene clusters called NK, Hox, and ParaHox. Comparative data supports the idea that NK genes were deployed primarily to pattern the bilaterian mesoderm, Hox genes coded position along the central nervous system, and ParaHox genes most likely originally specified the mouth, midgut, and anus of the newly evolved through-gut. It is proposed that diversification of ANTP class genes played a role in the Cambrian explosion by contributing to the patterning systems used to build animal bodies capable of high-energy directed locomotion, including active burrowing.

The circadian clock of teleost fish: a comparative analysis reveals distinct fates for duplicated genes

The circadian clock is a central oscillator that coordinates endogenous rhythms. Members of six gene families underlie the metabolic machinery of this system. Although this machinery appears to correspond to a highly conserved genetic system in metazoans, it has been recognized that vertebrates possess a more diverse gene inventory than that of non-vertebrates. This difference could have originated in the two successive rounds of whole-genome duplications that took place in the common ancestor of the group. Teleost fish underwent an extra event of whole-genome duplication, which is thought to have provided an abundance of raw genetic material for the biological innovations that facilitated the radiation of the group. In this study, we assessed the relative contributions of whole-genome duplication and small-scale gene duplication to generate the repertoire of genes associated with the circadian clock of teleost fish. To achieve this goal, we annotated genes from six gene families associated with the circadian clock in eight teleost fish species, and we reconstructed their evolutionary history by inferring phylogenetic relationships. Our comparative analysis indicated that teleost species possess a variable repertoire of genes related to the circadian clock gene families and that the actual diversity of these genes has been shaped by a variety of phenomena, such as the complete deletion of ohnologs, the differential retention of genes, and lineage-specific gene duplications. From a functional perspective, the subfunctionalization of two ohnolog genes (PER1a and PER1b) in zebrafish highlights the power of whole-genome duplications to generate biological diversity.

Genomic organisation of the seven ParaHox genes of coelacanths

Human and mouse genomes contain six ParaHox genes implicated in gut and neural patterning. In coelacanths and cartilaginous fish, an additional ParaHox gene exists—Pdx2—that dates back to the genome duplications in early vertebrate evolution. Here we examine the genomic arrangement and flanking genes of all ParaHox genes in coelacanths, to determine the full complement of these genes. We find that coelacanths have seven ParaHox genes in total, in four chromosomal locations, revealing that five gene losses occurred soon after vertebrate genome duplication. Comparison of intergenic sequences reveals that some Pdx1 regulatory regions associated with development of pancreatic islets are older than tetrapods, that Pdx1 and Pdx2 share few if any conserved non-coding elements, and that there is very high sequence conservation between coelacanth species.

Intact cluster and chordate-like expression of ParaHox genes in a sea star

Background

The ParaHox genes are thought to be major players in patterning the gut of several bilaterian taxa. Though this is a fundamental role that these transcription factors play, their activities are not limited to the endoderm and extend to both ectodermal and mesodermal tissues. Three genes compose the ParaHox group: Gsx, Xlox and Cdx. In some taxa (mostly chordates but to some degree also in protostomes) the three genes are arranged into a genomic cluster, in a similar fashion to what has been shown for the better-known Hox genes. Sea urchins possess the full complement of ParaHox genes but they are all dispersed throughout the genome, an arrangement that, perhaps, represented the primitive condition for all echinoderms. In order to understand the evolutionary history of this group of genes we cloned and characterized all ParaHox genes, studied their expression patterns and identified their genomic loci in a member of an earlier branching group of echinoderms, the asteroid Patiria miniata.

Results

We identified the three ParaHox orthologs in the genome of P. miniata. While one of them, PmGsx is provided as maternal message, with no zygotic activation afterwards, the other two, PmLox and PmCdx are expressed during embryogenesis, within restricted domains of both endoderm and ectoderm. Screening of a Patiria bacterial artificial chromosome (BAC) library led to the identification of a clone containing the three genes. The transcriptional directions of PmGsx and PmLox are opposed to that of the PmCdx gene within the cluster.

Conclusions

The identification of P. miniata ParaHox genes has revealed the fact that these genes are clustered in the genome, in contrast to what has been reported for echinoids. Since the presence of an intact cluster, or at least a partial cluster, has been reported in chordates and polychaetes respectively, it becomes clear that within echinoderms, sea urchins have modified the original bilaterian arrangement. Moreover, the sea star ParaHox domains of expression show chordate-like features not found in the sea urchin, confirming that the dynamics of gene expression for the respective genes and their putative regulatory interactions have clearly changed over evolutionary time within the echinoid lineage.

Identification of an intact ParaHox cluster with temporal colinearity but altered spatial colinearity in the hemichordate Ptychodera flava

BACKGROUND

ParaHox and Hox genes are thought to have evolved from a common ancestral ProtoHox cluster or from tandem duplication prior to the divergence of cnidarians and bilaterians. Similar to Hox clusters, chordate ParaHox genes including Gsx, Xlox, and Cdx, are clustered and their expression exhibits temporal and spatial colinearity. In non-chordate animals, however, studies on the genomic organization of ParaHox genes are limited to only a few animal taxa. Hemichordates, such as the Enteropneust acorn worms, have been used to gain insights into the origins of chordate characters. In this study, we investigated the genomic organization and expression of ParaHox genes in the indirect developing hemichordate acorn worm Ptychodera flava.

RESULTS

We found that P. flava contains an intact ParaHox cluster with a similar arrangement to that of chordates. The temporal expression order of the P. flava ParaHox genes is the same as that of the chordate ParaHox genes. During embryogenesis, the spatial expression pattern of PfCdx in the posterior endoderm represents a conserved feature similar to the expression of its orthologs in other animals. On the other hand, PfXlox and PfGsx show a novel expression pattern in the blastopore. Nevertheless, during metamorphosis, PfXlox and PfCdx are expressed in the endoderm in a spatially staggered pattern similar to the situation in chordates.

CONCLUSIONS

Our study shows that P. flava ParaHox genes, despite forming an intact cluster, exhibit temporal colinearity but lose spatial colinearity during embryogenesis. During metamorphosis, partial spatial colinearity is retained in the transforming larva. These results strongly suggest that intact ParaHox gene clustering was retained in the deuterostome ancestor and is correlated with temporal colinearity.

The cytochrome P450 genesis locus: the origin and evolution of animal cytochrome P450s

The neighbourhoods of cytochrome P450 (CYP) genes in deuterostome genomes, as well as those of the cnidarians Nematostella vectensis and Acropora digitifera and the placozoan Trichoplax adhaerens were examined to find clues concerning the evolution of CYP genes in animals. CYP genes created by the 2R whole genome duplications in chordates have been identified. Both microsynteny and macrosynteny were used to identify genes that coexisted near CYP genes in the animal ancestor. We show that all 11 CYP clans began in a common gene environment. The evidence implies the existence of a single locus, which we term the 'cytochrome P450 genesis locus', where one progenitor CYP gene duplicated to create a tandem set of genes that were precursors of the 11 animal CYP clans: CYP Clans 2, 3, 4, 7, 19, 20, 26, 46, 51, 74 and mitochondrial. These early CYP genes existed side by side before the origin of cnidarians, possibly with a few additional genes interspersed. The Hox gene cluster, WNT genes, an NK gene cluster and at least one ARF gene were close neighbours to this original CYP locus. According to this evolutionary scenario, the CYP74 clan originated from animals and not from land plants nor from a common ancestor of plants and animals. The CYP7 and CYP19 families that are chordate-specific belong to CYP clans that seem to have originated in the CYP genesis locus as well, even though this requires many gene losses to explain their current distribution. The approach to uncovering the CYP genesis locus overcomes confounding effects because of gene conversion, sequence divergence, gene birth and death, and opens the way to understanding the biodiversity of CYP genes, families and subfamilies, which in animals has been obscured by more than 600 Myr of evolution.

Whole-genome duplication and the functional diversification of teleost fish hemoglobins

Subsequent to the two rounds of whole-genome duplication that occurred in the common ancestor of vertebrates, a third genome duplication occurred in the stem lineage of teleost fishes. This teleost-specific genome duplication (TGD) is thought to have provided genetic raw materials for the physiological, morphological, and behavioral diversification of this highly speciose group. The extreme physiological versatility of teleost fish is manifest in their diversity of blood–gas transport traits, which reflects the myriad solutions that have evolved to maintain tissue O₂ delivery in the face of changing metabolic demands and environmental O₂ availability during different ontogenetic stages. During the course of development, regulatory changes in blood–O₂ transport are mediated by the expression of multiple, functionally distinct hemoglobin (Hb) isoforms that meet the particular O₂-transport challenges encountered by the developing embryo or fetus (in viviparous or oviparous species) and in free-swimming larvae and adults. The main objective of the present study was to assess the relative contributions of whole-genome duplication, large-scale segmental duplication, and small-scale gene duplication in producing the extraordinary functional diversity of teleost Hbs. To accomplish this, we integrated phylogenetic reconstructions with analyses of conserved synteny to characterize the genomic organization and evolutionary history of the globin gene clusters of teleosts. These results were then integrated with available experimental data on functional properties and developmental patterns of stage-specific gene expression. Our results indicate that multiple α- and β-globin genes were present in the common ancestor of gars (order Lepisoteiformes) and teleosts. The comparative genomic analysis revealed that teleosts possess a dual set of TGD-derived globin gene clusters, each of which has undergone lineage-specific changes in gene content via repeated duplication and deletion events. Phylogenetic reconstructions revealed that paralogous genes convergently evolved similar functional properties in different teleost lineages. Consistent with other recent studies of globin gene family evolution in vertebrates, our results revealed evidence for repeated evolutionary transitions in the developmental regulation of Hb synthesis.

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.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

eGOB: eukaryotic Gene Order Browser

A large number of genomes have been sequenced, allowing a range of comparative studies. Here, we present the eukaryotic Gene Order Browser with information on the order of protein and non-coding RNA (ncRNA) genes of 74 different eukaryotic species. The browser is able to display a gene of interest together with its genomic context in all species where that gene is present. Thereby, questions related to the evolution of gene organization and non-random gene order may be examined. The browser also provides access to data collected on pairs of adjacent genes that are evolutionarily conserved.

AVAILABILITY

eGOB as well as underlying data are freely available at http://egob.biomedicine.gu.se

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

tore.samuelsson@medkem.gu.se.

Homeodomain subtypes and functional diversity

The homeodomain is a protein domain of about 60 amino acids that is encoded by homeobox genes. The homeodomain is a DNA binding domain, and hence homeodomain proteins are essentially transcription factors (TFs). They have been shown to play major roles in many developmental processes of animals, as well as fungi and plants. A primary function of homeodomain proteins is to regulate the expression of other genes in development and differentiation. Thousands of homeobox genes have been identified, and they can be grouped into many different classes. Often other conserved protein domains are found linked to a homeodomain. Several particular types of homeobox genes are organized into chromosomal clusters. The best-known cluster, the HOX cluster, is found in all bilaterian animals. Tetrapods contain four HOX clusters that arose through duplication in early vertebrate evolution. The genes in these clusters are called Hox genes. Lower chordates, insects and nematodes tend to have only one HOX cluster. Of particular interest is that many of the HOX cluster genes function in the process of pattern formation along the anterior-posterior body axis. Many other types of homeodomain proteins play roles in the determination of cell fates and cell differentiation. Homeobox genes thus perform key roles for all aspects of the development of an organism.

Retinoic acid-dependent establishment of positional information in the hindbrain was conserved during vertebrate evolution

Zebrafish hoxb1b is expressed during epiboly in the posterior neural plate, with its anterior boundary at the prospective r4 region providing a positional cue for hindbrain formation. A similar function and expression is known for Hoxa1 in mice, suggesting a shared regulatory mechanism for hindbrain patterning in vertebrate embryos. To understand the evolution of the regulatory mechanisms of key genes in patterning of the central nervous system, we examined how hoxb1b transcription is regulated in zebrafish embryos and compared the regulatory mechanisms between mammals and teleosts that have undergone an additional genome duplication. By promoter analysis, we found that the expression of the reporter gene recapitulated hoxb1b expression when driven in transgenic embryos by a combination of the upstream 8.0-kb DNA and downstream 4.6-kb DNA. Furthermore, reporter expression expanded anteriorly when transgenic embryos were exposed to retinoic acid (RA) or LiCl, or injected with fgf3/8 mRNA, implicating the flanking DNA examined here in the responsiveness of hoxb1b to posteriorizing signals. We further identified at least two functional RA responsive elements in the downstream DNA that were shown to be major regulators of early hoxb1b expression during gastrulation, while the upstream DNA, which harbors repetitive sequences with apparent similarity to the autoregulatory sequence of mouse Hoxb1, contributed only to later hoxb1b expression, during somitogenesis. Possible implications in vertebrate evolution are discussed based on these findings.

Conservation of ParaHox genes' function in patterning of the digestive tract of the marine gastropod Gibbula varia

BACKGROUND

Presence of all three ParaHox genes has been described in deuterostomes and lophotrochozoans, but to date one of these three genes, Xlox has not been reported from any ecdysozoan taxa and both Xlox and Gsx are absent in nematodes. There is evidence that the ParaHox genes were ancestrally a single chromosomal cluster. Colinear expression of the ParaHox genes in anterior, middle, and posterior tissues of several species studied so far suggest that these genes may be responsible for axial patterning of the digestive tract. So far, there are no data on expression of these genes in molluscs.

RESULTS

We isolated the complete coding sequences of the three Gibbula varia ParaHox genes, and then tested their expression in larval and postlarval development. In Gibbula varia, the ParaHox genes participate in patterning of the digestive tract and are expressed in some cells of the neuroectoderm. The expression of these genes coincides with the gradual formation of the gut in the larva. Gva-Gsx patterns potential neural precursors of cerebral ganglia as well as of the apical sensory organ. During larval development this gene is involved in the formation of the mouth and during postlarval development it is expressed in the precursor cells involved in secretion of the radula, the odontoblasts. Gva-Xolx and Gva-Cdx are involved in gut patterning in the middle and posterior parts of digestive tract, respectively. Both genes are expressed in some ventral neuroectodermal cells; however the expression of Gva-Cdx fades in later larval stages while the expression of Gva-Xolx in these cells persists.

CONCLUSIONS

In Gibbula varia the ParaHox genes are expressed during anterior-posterior patterning of the digestive system. This colinearity is not easy to spot during early larval stages because the differentiated endothelial cells within the yolk permanently migrate to their destinations in the gut. After torsion, Gsx patterns the mouth and foregut, Xlox the midgut gland or digestive gland, and Cdx the hindgut. ParaHox genes of Gibbula are also expressed during specification of cerebral and ventral neuroectodermal cells. Our results provide additional support for the ancestral complexity of Gsx expression and its ancestral role in mouth patterning in protostomes, which was secondarily lost or simplified in some species.

Evolution of Hox gene clusters in gnathostomes: insights from a survey of a shark (Scyliorhinus canicula) transcriptome

It is now well established that there were four Hox gene clusters in the genome of the last common ancestor of extant gnathostomes. To better understand the evolution of the organization and expression of these genomic regions, we have studied the Hox gene clusters of a shark (Scyliorhinus canicula). We sequenced 225,580 expressed sequence tags from several embryonic cDNA libraries. Blast searches identified corresponding transcripts to almost all the HoxA, HoxB, and HoxD cluster genes. No HoxC transcript was identified, suggesting that this cluster is absent or highly degenerate. Using Hox gene sequences as probes, we selected and sequenced seven clones from a bacterial artificial chromosome library covering the complete region of the three gene clusters. Mapping of cDNAs to these genomic sequences showed extensive alternative splicing and untranslated exon sharing between neighboring Hox genes. Homologous noncoding exons could not be identified in transcripts from other species using sequence similarity. However, by comparing conserved noncoding sequences upstream of these exons in different species, we were able to identify homology between some exons. Some alternative splicing variants are probably very ancient and were already coded for by the ancestral Hox gene cluster. We also identified several transcripts that do not code for Hox proteins, are probably not translated, and all but one are in the reverse orientation to the Hox genes. This survey of the transcriptome of the Hox gene clusters of a shark shows that the high complexity observed in mammals is a gnathostome ancestral feature.

Parallel retention of Pdx2 genes in cartilaginous fish and coelacanths

The Pdx1 or Ipf1 gene encodes an important homeodomain-containing protein with key roles in pancreas development and function. Mutations in human PDX1 are implicated in developmental defects and disease of the pancreas. Extensive research, including genome sequencing, has indicated that Pdx1 is the only member of its gene family in mammals, birds, amphibians, and ray-finned fish, and with the exception of teleost fish, this gene forms part of the ParaHox gene cluster along with Gsx1 and Cdx2. The ParaHox cluster, however, is a remnant of a 4-fold genome duplication; the three other ParaHox paralogues lack a Pdx-like gene in all vertebrate genomes examined to date. We have used bacterial artificial chromosome cloning and synteny analysis to show that the ancestor of living jawed vertebrates in fact had more ParaHox genes, including two Pdx genes (Pdx1 and Pdx2). Surprisingly, the two Pdx genes have been retained in parallel in two quite distantly related lineages, the cartilaginous fish (sharks, skates, and chimeras) and the Indonesian coelacanth, Latimeria menadoensis. The Pdx2 gene has been lost independently in ray-finned fish and in tetrapods.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Lineage-specific co-evolution of the Egf receptor/ligand signaling system

Background

The epidermal growth factor receptor (Egfr) with its numerous ligands has fundamental roles in development, cell differentiation and physiology. Dysfunction of the receptor-ligand system contributes to many human malignancies. Consistent with such various tasks, the Egfr gene family has expanded during vertebrate evolution as a consequence of several rounds of whole genome duplication. Of particular interest is the effect of the fish-specific whole genome duplication (FSGD) on the ligand-receptor system, as it has supplied this largest group of vertebrates with additional opportunities for sub- and/or neofunctionalization in this signaling system.

Results

We identified the predicted components of the Egf receptor-ligand signaling system in teleost fishes (medaka, platyfish, stickleback, pufferfishes and zebrafish). We found two duplicated egfr genes, egfra and egfrb, in all available teleost genomes. Surprisingly only one copy for each of the seven Egfr ligands could be identified in most fishes, with zebrafish hbegf being the only exception. Special focus was put on medaka, for which we more closely investigated all Egf receptors and Egfr ligands. The different expression patterns of egfra, egfrb and their ligands in medaka tissues and embryo stages suggest differences in role and function. Preferential co-expression of different subsets of Egfr ligands corroborates the possible subfunctionalization and specialization of the two receptors in adult tissues. Bioinformatic analyses of the ligand-receptor interface between Egfr and its ligands show a very weak evolutionary conservation within this region. Using in vitro analyses of medaka Egfra, we could show that this receptor is only activated by medaka ligands, but not by human EGF. Altogether, our data suggest a lineage-specific Egfr/Egfr ligand co-evolution.

Conclusions

Our data indicate that medaka Egfr signaling occurs via its two copies, Egfra and Egfrb, each of them being preferentially coexpressed with different subsets of Egfr ligands. This fish-specific occurrence of Egf receptor specialization offers unique opportunities to study the functions of different Egf receptor-ligand combinations and their biological outputs in vertebrates. Furthermore, our results strongly support the use of homologous ligands in future studies, as sufficient cross-specificity is very unlikely for this ligand/receptor system.

Features of the ancestral bilaterian inferred from Platynereis dumerilii ParaHox genes

Background

The ParaHox gene cluster is the evolutionary sister to the Hox cluster. Whilst the role of the Hox cluster in patterning the anterior-posterior axis of bilaterian animals is well established, and the organisation of vertebrate Hox clusters is intimately linked to gene regulation, much less is known about the more recently discovered ParaHox cluster. ParaHox gene clustering, and its relationship to expression, has only been described in deuterostomes. Conventional protostome models (Drosophila melanogaster and Caenorhabditis elegans) are secondarily derived with respect to ParaHox genes, suffering gene loss and cluster break-up.

Results

We provide the first evidence for ParaHox gene clustering from a less-derived protostome animal, the annelid Platynereis dumerilii. Clustering of these genes is thus not a sole preserve of the deuterostome lineage within Bilateria. This protostome ParaHox cluster is not entirely intact however, with Pdu-Cdx being on the opposite end of the same chromosome arm from Pdu-Gsx and Pdu-Xlox. From the genomic sequence around the P. dumerilii ParaHox genes the neighbouring genes are identified, compared with other taxa, and the ancestral arrangement deduced.

Conclusion

We relate the organisation of the ParaHox genes to their expression, and from comparisons with other taxa hypothesise that a relatively complex pattern of ParaHox gene expression existed in the protostome-deuterostome ancestor, which was secondarily simplified along several invertebrate lineages. Detailed comparisons of the gene content around the ParaHox genes enables the reconstruction of the genome surrounding the ParaHox cluster of the protostome-deuterostome ancestor, which existed over 550 million years ago.

Overlapping functions of Cdx1, Cdx2, and Cdx4 in the development of the amphibian Xenopus tropicalis

Using Xenopus tropicalis, we present the first analysis of the developmental effects that result from knocking down the function of the three Cdx genes present in the typical vertebrate genome. Knockdowns of individual Cdx genes lead to a similar range of posterior defects; compound Cdx knockdowns result in increasingly severe posterior truncations, accompanied by posterior shifts and reduction of 5' Hox gene expression. We provide evidence that Cdx and Wnt3A genes are components of a positive feedback loop operating in the posterior axis. We show that Cdx function is required during later, but not early stages of development, for correct regional specification of the endoderm and morphogenesis of the gut. Our results support the hypothesis that during amphibian development the overall landscape of Cdx activity in the embryo is more important than the specific function of individual Cdx proteins.

ParaHox cluster evolution--hagfish and beyond

Abstract The ParaHox genes comprise three Hox-related homeobox gene families, found throughout the animals. They were first discovered in the invertebrate chordate amphioxus, where they are tightly clustered. In this paper we carry out a comparative review of ParaHox gene cluster organization among the deuterostomes, and discuss how the recently published hagfish ParaHox clusters fit into current theories about the evolution of this group of genes.

Hox cluster duplication in the basal teleost Hiodon alosoides (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.

Do cnidarians have a ParaHox cluster? Analysis of synteny around a Nematostella homeobox gene cluster

The Hox gene cluster is renowned for its role in developmental patterning of embryogenesis along the anterior-posterior axis of bilaterians. Its supposed evolutionary sister or paralog, the ParaHox cluster, is composed of Gsx, Xlox, and Cdx, and also has important roles in anterior-posterior development. There is a debate as to whether the cnidarians, as an outgroup to bilaterians, contain true Hox and ParaHox genes, or instead the Hox-like gene complement of cnidarians arose from independent duplications to those that generated the genes of the bilaterian Hox and ParaHox clusters. A recent whole genome analysis of the cnidarian Nematostella vectensis found conserved synteny between this cnidarian and vertebrates, including a region of synteny between the putative Hox cluster of N. vectensis and the Hox clusters of vertebrates. No syntenic region was identified around a potential cnidarian ParaHox cluster. Here we use different approaches to identify a genomic region in N. vectensis that is syntenic with the bilaterian ParaHox cluster. This proves that the duplication that gave rise to the Hox and ParaHox regions of bilaterians occurred before the origin of cnidarians, and the cnidarian N. vectensis has bona fide Hox and ParaHox loci.

Intestinal differentiation in zebrafish requires Cdx1b, a functional equivalent of mammalian Cdx2

BACKGROUND & AIMS

The ParaHox transcription factor Cdx2 is an essential determinant of intestinal phenotype in mammals throughout development, influencing gut function, homeostasis, and epithelial barrier integrity. Cdx2 expression demarcates the zones of intestinal stem cell proliferation in the adult gut, with deregulated expression implicated in intestinal metaplasia and cancer. However, in vivo analysis of these prospective roles has been limited because inactivation of Cdx2 in mice leads to preimplantation embryonic lethality. We used the zebrafish, a valuable model for studying gut development, to generate a system to further understanding of the role of Cdx2 in normal intestinal function and in disease states.

METHODS

We isolated and characterized the zebrafish cdx1b ortholog and analyzed its function by antisense morpholino gene knockdown.

RESULTS

We showed that zebrafish Cdx1b replaces the role of Cdx2 in gut development. Evolutionary studies have indicated that the zebrafish cdx2 loci were lost following the genome-wide duplication event that occurred in teleosts. Zebrafish Cdx1b is expressed exclusively in the developing intestine during late embryogenesis and regulates intestinal cell proliferation and terminal differentiation.

CONCLUSIONS

This work established an in vivo system to explore further the activity of Cdx2 in the gut and its impact on processes such as inflammation and cancer.

SynBlast: assisting the analysis of conserved synteny information

Motivation

In the last years more than 20 vertebrate genomes have been sequenced, and the rate at which genomic DNA information becomes available is rapidly accelerating. Gene duplication and gene loss events inherently limit the accuracy of orthology detection based on sequence similarity alone. Fully automated methods for orthology annotation do exist but often fail to identify individual members in cases of large gene families, or to distinguish missing data from traceable gene losses. This situation can be improved in many cases by including conserved synteny information.

Results

Here we present the SynBlast pipeline that is designed to construct and evaluate local synteny information. SynBlast uses the genomic region around a focal reference gene to retrieve candidates for homologous regions from a collection of target genomes and ranks them in accord with the available evidence for homology. The pipeline is intended as a tool to aid high quality manual annotation in particular in those cases where automatic procedures fail. We demonstrate how SynBlast is applied to retrieving orthologous and paralogous clusters using the vertebrate Hox and ParaHox clusters as examples.

Software

The SynBlast package written in Perl is available under the GNU General Public License at .

Comparative genomics of Lbx loci reveals conservation of identical Lbx ohnologs in bony vertebrates

Background

Lbx/ladybird genes originated as part of the metazoan cluster of Nk homeobox genes. In all animals investigated so far, both the protostome genes and the vertebrate Lbx1 genes were found to play crucial roles in neural and muscle development. Recently however, additional Lbx genes with divergent expression patterns were discovered in amniotes. Early in the evolution of vertebrates, two rounds of whole genome duplication are thought to have occurred, during which 4 Lbx genes were generated. Which of these genes were maintained in extant vertebrates, and how these genes and their functions evolved, is not known.

Results

Here we searched vertebrate genomes for Lbx genes and discovered novel members of this gene family. We also identified signature genes linked to particular Lbx loci and traced the remnants of 4 Lbx paralogons (two of which retain Lbx genes) in amniotes. In teleosts, that have undergone an additional genome duplication, 8 Lbx paralogons (three of which retain Lbx genes) were found. Phylogenetic analyses of Lbx and Lbx-associated genes show that in extant, bony vertebrates only Lbx1- and Lbx2-type genes are maintained. Of these, some Lbx2 sequences evolved faster and were probably subject to neofunctionalisation, while Lbx1 genes may have retained more features of the ancestral Lbx gene. Genes at Lbx1 and former Lbx4 loci are more closely related, as are genes at Lbx2 and former Lbx3 loci. This suggests that during the second vertebrate genome duplication, Lbx1/4 and Lbx2/3 paralogons were generated from the duplicated Lbx loci created during the first duplication event.

Conclusion

Our study establishes for the first time the evolutionary history of Lbx genes in bony vertebrates, including the order of gene duplication events, gene loss and phylogenetic relationships. Moreover, we identified genetic hallmarks for each of the Lbx paralogons that can be used to trace Lbx genes as other vertebrate genomes become available. Significantly, we show that bony vertebrates only retained copies of Lbx1 and Lbx2 genes, with some Lbx2 genes being highly divergent. Thus, we have established a base on which the evolution of Lbx gene function in vertebrate development can be evaluated.

ParaHox gene expression in larval and postlarval development of the polychaete Nereis virens (Annelida, Lophotrochozoa)

Background

Transcription factors that encode ANTP-class homeobox genes play crucial roles in determining the body plan organization and specification of different organs and tissues in bilaterian animals. The three-gene ParaHox family descends from an ancestral gene cluster that existed before the evolution of the Bilateria. All three ParaHox genes are reported from deuterostomes and lophotrochozoans, but not to date from any ecdysozoan taxa, and there is evidence that the ParaHox genes, like the related Hox genes, were ancestrally a single chromosomal cluster. However, unlike the Hox genes, there is as yet no strong evidence that the ParaHox genes are expressed in spatial and temporal order during embryogenesis.

Results

We isolated fragments of the three Nereis virens ParaHox genes, then used these as probes for whole-mount in situ hybridization in larval and postlarval worms. In Nereis virens the ParaHox genes participate in antero-posterior patterning of ectodermal and endodermal regions of the digestive tract and are expressed in some cells in the segment ganglia. The expression of these genes occurs in larval development in accordance with the position of these cells along the main body axis and in postlarval development in accordance with the position of cells in ganglia along the antero-posterior axis of each segment. In none of these tissues does expression of the three ParaHox genes follow the rule of temporal collinearity.

Conclusion

In Nereis virens the ParaHox genes are expressed during antero-posterior patterning of the digestive system (ectodermal foregut and hindgut, and endodermal midgut) of Nereis virens. These genes are also expressed during axial specification of ventral neuroectodermal cell domains, where the expression domains of each gene are re-iterated in each neuromere except for the first parapodial segment. These expression domains are probably predetermined and may be directed on the antero-posterior axis by the Hox genes, whose expression starts much earlier during embryogenesis. Our results support the hypothesis that the ParaHox genes are involved in antero-posterior patterning of the developing embryo, but they do not support the notion that these genes function only in the patterning of endodermal tissues.

ParaHox genes in pancreatic cell cultures: effects on the insulin promoter regulation

The gene encoding PDX1 (pancreatic duodenum homeobox 1), the main transcription factor regulating the glucose-dependent transactivation of the insulin promoter in pancreatic β-cells, clusters with two closely related homeobox genes (Gsh1 and Cdx2/3), all of them belonging to the ParaHox gene family. The ParaHox gene evolutionary history in the vertebrate lineage involved duplications of the cluster and subsequent loss of some members, so that eventually, the human and murine genomes contain only 6 ParaHox genes. The crucial role of PDX1 in pancreas development, beta-cell formation and insulin transcription regulation has long been established. There is some data on CDX2/3 function in α-cells, but remarkably, nothing is known on the role of the other ParaHox genes, which are also expressed in the endocrine pancreas. Homeobox transcription factors that belong to the same family show high conservation of the homeodomain and share similar target sites and oligomeric partners, and thus may act redundantly, synergistically or antagonistically on the same promoters. Therefore, we explored the effects of the Parahox proteins (GSH1, GSH2, CDX1, CDX2/3 and CDX4) on the regulation of the insulin promoter in transfected α- and β- cultured cell lines at different glucose concentrations and compared them to those of PDX1. Noticeably, several ParaHox transcription factors are able to transactivate or inhibit the insulin promoter, depending on the cell type and glucose concentration, thus suggesting their possible participation in the regulation of similar target genes, such as insulin, either by silencing or activating them, in the absence of PDX1.

Cdx4 is required in the endoderm to localize the pancreas and limit beta-cell number

Cdx transcription factors have crucial roles in anteroposterior patterning of the nervous system and mesoderm. Here we focus on the role of cdx4 in patterning the endoderm in zebrafish. We show that cdx4 has roles in determining pancreatic beta-cell number, directing midline convergence of beta-cells during early pancreatic islet formation, and specifying the anteroposterior location of foregut organs. Embryos deficient in cdx4 have a posteriorly shifted pancreas, liver and small intestine. The phenotype is more severe with knockdown of an additional Cdx factor, cdx1a. We show that cdx4 functions within the endoderm to localize the pancreas. Morpholino knockdown of cdx4 specifically in the endoderm recapitulates the posteriorly shifted pancreas observed in cdx4 mutants. Conversely, overexpression of cdx4 specifically in the endoderm is sufficient to shift the pancreas anteriorly. Together, these results suggest a model in which cdx4 confers posterior identity to the endoderm. Cdx4 might function to block pancreatic identity by preventing retinoic acid (RA) signal transduction in posterior endoderm. In support of this, we demonstrate that in cdx4-deficient embryos treated with RA, ectopic beta-cells are located well posterior to the normal pancreatic domain.

Retroviral enhancer detection insertions in zebrafish combined with comparative genomics reveal genomic regulatory blocks - a fundamental feature of vertebrate genomes

A large-scale enhancer detection screen was performed in the zebrafish using a retroviral vector carrying a basal promoter and a fluorescent protein reporter cassette. Analysis of insertional hotspots uncovered areas around developmental regulatory genes in which an insertion results in the same global expression pattern, irrespective of exact position. These areas coincide with vertebrate chromosomal segments containing identical gene order; a phenomenon known as conserved synteny and thought to be a vestige of evolution. Genomic comparative studies have found large numbers of highly conserved noncoding elements (HCNEs) spanning these and other loci. HCNEs are thought to act as transcriptional enhancers based on the finding that many of those that have been tested direct tissue specific expression in transient or transgenic assays. Although gene order in hox and other gene clusters has long been known to be conserved because of shared regulatory sequences or overlapping transcriptional units, the chromosomal areas found through insertional hotspots contain only one or a few developmental regulatory genes as well as phylogenetically unrelated genes. We have termed these regions genomic regulatory blocks (GRBs), and show that they underlie the phenomenon of conserved synteny through all sequenced vertebrate genomes. After teleost whole genome duplication, a subset of GRBs were retained in two copies, underwent degenerative changes compared with tetrapod loci that exist as single copy, and that therefore can be viewed as representing the ancestral form. We discuss these findings in light of evolution of vertebrate chromosomal architecture and the identification of human disease mutations.