Phylogenetic Dating and Characterization of Gene Duplications in Vertebrates: The Cartilaginous Fish Reference

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

The evolution and immunomodulatory role of Zc3h12 proteins in zebrafish (Danio rerio)

Zc3h12 family is an important RNA-binding protein family regulating mRNA of inflammatory cytokines in mammals. However, there are few studies on their post-transcriptional level regulation of inflammatory cytokines in fish. Here, we investigated the evolution of zebrafish Zc3h12 family and explored their immunomodulatory role. Phylogenetic and syntenic analysis indicated the number of zc3h12 family members had increased ranging from a single member in invertebrates to a single copy of four members in mammals. As the most evolutionarily diverse group of vertebrates, the number of zc3h12 family members was more complex and diverse in the teleost, each member experienced different fates and followed different rules in multiple rounds of whole-genome duplication events. Thereinto, zebrafish contained three zc3h12 genes, among which zc3h12aa and zc3h12ab were duplicated from the same gene. Zebrafish Zc3h12 family could recognize the 3'-UTR regions of inflammatory cytokines through binding to the specific RNA secondary structure and negatively regulate their expression. Deletion of either Zc3h12 domains or mutation of the key amino acid in RNAase domain attenuated their modulatory effect, suggesting both domain and RNAase activity are important to the immunomodulatory role. These results elucidated the evolution of Zc3h12 family and uncovered Zc3h12-mediated post-transcriptional regulation of cytokines in zebrafish.

Genetic dissection of steroid-hormone modulated social behavior: Novel paralogous genes are a boon for discovery

Research across species has led to important discoveries on the functions of steroid hormones in the regulation of behavior. However, like in many fields, advancements in transgenic and mutagenic technology allowed mice to become the premier genetic model for conducting many experiments to understand how steroids control social behavior. Since there has been a general lack of parallel methodological developments in other species, many of the findings cannot be generalized. This is especially the case for teleost fish, in which a whole-genome duplication produced novel paralogs for key steroid hormone signaling genes. In this review, we summarize technical advancements over the history of the field of neuroendocrinology that have led to important insights in our understanding of the control of social behavior by steroids. We demonstrate that early mouse genetic models to understand these mechanisms suffered from several issues that were remedied by more precise transgenic technological advancements. We then highlight the importance of CRISPR/Cas9 gene editing tools that will in time bridge the gap between mice and non-traditional model species for understanding principles of steroid hormone action in the modulation of social behavior. We specifically highlight the role of teleost fish in bridging this gap because they are 1) highly genetically tractable and 2) provide a novel advantage in achieving precise genetic control. The field of neuroendocrinology is entering a new "gene editing revolution" that will lead to novel discoveries about the roles of steroid hormones in the regulation and evolutionary trajectories of social behavior.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Physiological relevance of nitric oxide in ovarian functions: An overview

Nitric oxide (NO, nitrogen monoxide), a short-lived, free radical carrying an unpaired electron, is one of the smallest molecules synthesized in the biological system. In addition to its role in angiogenesis, neuronal function and inflammatory response, NO has wide-spread significance in regulation of ovarian function in vertebrates. Based on tissue-specific expression, three different nitric oxide synthase (NOS) isoforms, neuronal (nNOS) or NOS1, inducible (iNOS) or NOS2 and endothelial (eNOS) or NOS3 have been identified. While expression of both inducible (iNOS) and constitutive NOS (eNOS) isoforms varies considerably in the ovary at various stages of follicular growth and development, selective binding of NO with proteins containing heme moieties have significant influence on ovarian steroidogenesis. Besides, NO modulation of ovulatory response suggests physiological significance of NO/NOS system in mammalian ovary. Compared to the duality of NO action on follicular development, steroidogenesis and meiotic maturation in mammalian models, participation of NO/NOS system in teleost ovary is less investigated. Genes encoding nos1 and nos2 have been identified in fish; however, presence of nos3 is still ambiguous. Interestingly, two distinct nos2 genes, nos2a and nos2b in zebrafish, possibly arose through whole genome duplication. Differential expression of major NOS isoforms in catfish ovary, NO inhibition of meiosis resumption in Anabas testudineus follicle-enclosed oocytes and NO/sGC/cGMP modulation of oocyte maturation in zebrafish are some of the recent advancements. The present overview is an update on the advancements made and shortfalls still remaining in NO/NOS modulation of intercellular communication in teleost vis-à-vis mammalian ovary.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BACKGROUND

Gene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.

RESULTS

Based on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.

CONCLUSIONS

Here, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects-i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements-suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

Role of elasmobranchs and holocephalans in understanding peptide evolution in the vertebrates: Lessons learned from gonadotropin releasing hormone (GnRH) and corticotropin releasing factor (CRF) phylogenies

The cartilaginous fishes (Class Chondrichthyes) comprise two morphologically distinct subclasses; Elasmobranchii and Holocephali. Evidence indicates early divergence of these subclasses, suggesting monophyly of their lineage. However, such a phylogenetic understanding is not yet developed within two highly conserved peptide lineages, GnRH and CRF. Various GnRH forms exist across the Chondrichthyes. Although 4-7 immunoreactive forms have been described in Elasmobranchii, only one has been elucidated in Holocephali. In contrast, Chondrichthyan CRF phylogeny follows a pattern more consistent with vertebrate evolution. For example, three forms are expressed within the lamprey, with similar peptides present within the genome of the Callorhinchus milii, a holocephalan. Although these findings are consistent with recent evidence regarding the phylogenetic age of Chondrichthyan lineages, CRF evolution in vertebrates remains elusive. Assuming that the Elasmobranchii and Holocephali are part of a monocladistic clade within the Chondrichthyes, we interpret the findings of GnRH and CRF to be products of their respective lineages.

Evolutionary history of glucose-6-phosphatase encoding genes in vertebrate lineages: towards a better understanding of the functions of multiple duplicates

Background

Glucose-6-phosphate (G6pc) is a key enzyme involved in the regulation of the glucose homeostasis. The present study aims at revisiting and clarifying the evolutionary history of g6pc genes in vertebrates.

Results

g6pc duplications happened by successive rounds of whole genome duplication that occurred during vertebrate evolution. g6pc duplicated before or around Osteichthyes/Chondrichthyes radiation, giving rise to g6pca and g6pcb as a consequence of the second vertebrate whole genome duplication. g6pca was lost after this duplication in Sarcopterygii whereas both g6pca and g6pcb then duplicated as a consequence of the teleost-specific whole genome duplication. One g6pca duplicate was lost after this duplication in teleosts. Similarly one g6pcb2 duplicate was lost at least in the ancestor of percomorpha. The analysis of the evolution of spatial expression patterns of g6pc genes in vertebrates showed that all g6pc were mainly expressed in intestine and liver whereas teleost-specific g6pcb2 genes were mainly and surprisingly expressed in brain and heart. g6pcb2b, one gene previously hypothesised to be involved in the glucose intolerant phenotype in trout, was unexpectedly up-regulated (as it was in liver) by carbohydrates in trout telencephalon without showing significant changes in other brain regions. This up-regulation is in striking contrast with expected glucosensing mechanisms suggesting that its positive response to glucose relates to specific unknown processes in this brain area.

Conclusions

Our results suggested that the fixation and the divergence of g6pc duplicated genes during vertebrates’ evolution may lead to adaptive novelty and probably to the emergence of novel phenotypes related to glucose homeostasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3727-1) contains supplementary material, which is available to authorized users.

Evolutionary diversification of retinoic acid receptor ligand-binding pocket structure by molecular tinkering

Whole genome duplications (WGDs) have been classically associated with the origin of evolutionary novelties and the so-called duplication-degeneration-complementation model describes the possible fates of genes after duplication. However, how sequence divergence effectively allows functional changes between gene duplicates is still unclear. In the vertebrate lineage, two rounds of WGDs took place, giving rise to paralogous gene copies observed for many gene families. For the retinoic acid receptors (RARs), for example, which are members of the nuclear hormone receptor (NR) superfamily, a unique ancestral gene has been duplicated resulting in three vertebrate paralogues: RARα, RARβ and RARγ. It has previously been shown that this single ancestral RAR was neofunctionalized to give rise to a larger substrate specificity range in the RARs of extant jawed vertebrates (also called gnathostomes). To understand RAR diversification, the members of the cyclostomes (lamprey and hagfish), jawless vertebrates representing the extant sister group of gnathostomes, provide an intermediate situation and thus allow the characterization of the evolutionary steps that shaped RAR ligand-binding properties following the WGDs. In this study, we assessed the ligand-binding specificity of cyclostome RARs and found that their ligand-binding pockets resemble those of gnathostome RARα and RARβ. In contrast, none of the cyclostome receptors studied showed any RARγ-like specificity. Together, our results suggest that cyclostome RARs cover only a portion of the specificity repertoire of the ancestral gnathostome RARs and indicate that the establishment of ligand-binding specificity was a stepwise event. This iterative process thus provides a rare example for the diversification of receptor-ligand interactions of NRs following WGDs.

Membrane-associated collagens with interrupted triple-helices (MACITs): evolution from a bilaterian common ancestor and functional conservation in C. elegans

BACKGROUND

Collagens provide structural support and guidance cues within the extracellular matrix of metazoans. Mammalian collagens XIII, XXIII and XXV form a unique subgroup of type II transmembrane proteins, each comprising a short N-terminal cytosolic domain, a transmembrane domain and a largely collagenous ectodomain. We name these collagens as MACITs (Membrane-Associated Collagens with Interrupted Triple-helices), and here investigate their evolution and conserved properties. To date, these collagens have been studied only in mammals. Knowledge of the representation of MACITs in other extant metazoans is lacking. This question is of interest for understanding structural/functional relationships in the MACIT family and also for insight into the evolution of MACITs in relation to the secreted, fibrillar collagens that are present throughout the metazoa.

RESULTS

MACITs are restricted to bilaterians and are represented in the Ecdysozoa, Hemichordata, Urochordata and Vertebrata (Gnathostomata). They were not identified in available early-diverging metazoans, Lophotrochozoa, Echinodermata, Cephalochordata or Vertebrata (Cyclostomata). Whereas invertebrates encode a single MACIT, collagens XIII/XXIII/XXV of jawed vertebrates are paralogues that originated from the two rounds of en-bloc genome duplication occurring early in vertebrate evolution. MACITs have conserved domain architecture in which a juxta-membrane furin-cleavage site and the C-terminal 34 residues are especially highly conserved, whereas the cytoplasmic domains are weakly conserved. To study protein expression and function in a metazoan with a single MACIT gene, we focused on Caenorhabditis elegans and its col-99 gene. A col-99 cDNA was cloned and expressed as protein in mammalian CHO cells, two antibodies against COL-99 protein were generated, and a col-99-bearing fosmid gene construct col-99::egfp::flag was used to generate transgenic C. elegans lines. The encoded COL-99 polypeptide is 85 kDa in size and forms a trimeric protein. COL-99 is plasma membrane-associated and undergoes furin-dependent ectodomain cleavage and shedding. COL-99 is detected in mouth, pharynx, body wall and the tail, mostly in motor neurons and muscle systems and is enriched at neuromuscular junctions.

CONCLUSIONS

Through identification of MACITs in multiple metazoan phyla we developed a model for the evolution of MACITs. The experimental data demonstrate conservation of MACIT molecular and cellular properties and tissue localisations in the invertebrate, C. elegans.

Neofunctionalization of Androgen Receptor by Gain-of-Function Mutations in Teleost Fish Lineage

Steroid hormone receptor family provides an example of evolution of diverse transcription factors through whole-genome duplication (WGD). However, little is known about how their functions have been evolved after the duplication. Teleosts present a good model to investigate an accurate evolutionary history of protein function after WGD, because a teleost-specific WGD (TSGD) resulted in a variety of duplicated genes in modern fishes. This study focused on the evolution of androgen receptor (AR) gene, as two distinct paralogs, ARα and ARβ, have evolved in teleost lineage after TSGD. ARα showed a unique intracellular localization with a higher transactivation response than that of ARβ. Using site-directed mutagenesis and computational prediction of protein-ligand interactions, we identified two key substitutions generating a new functionality of euteleost ARα. The substitution in the hinge region contributes to the unique intracellular localization of ARα. The substitution on helices 10/11 in the ligand-binding domain possibly modulates hydrogen bonds that stabilize the receptor-ligand complex leading to the higher transactivation response of ARα. These substitutions were conserved in Acanthomorpha (spiny-rayed fish) ARαs, but not in an earlier branching lineage among teleosts, Japanese eel. Insertion of these substitutions into ARs from Japanese eel recapitulates the evolutionary novelty of euteleost ARα. These findings together indicate that the substitutions generating a new functionality of teleost ARα were fixed in teleost genome after the divergence of the Elopomorpha lineage. Our findings provide a molecular explanation for an adaptation process leading to generation of the hyperactive AR subtype after TSGD.

The evolution of tenascins and fibronectin

Tenascins are extracellular matrix glycoproteins that act both as integrin ligands and as modifiers of fibronectin-integrin interactions to regulate cell adhesion, migration, proliferation and differentiation. In tetrapods, both tenascins and fibronectin bind to integrins via RGD and LDV-type tripeptide motifs found in exposed loops in their fibronectin-type III domains. We previously showed that tenascins appeared early in the chordate lineage and are represented by single genes in extant cephalochordates and tunicates. Here we have examined the genomes of the coelacanth Latimeria chalumnae, the elephant shark Callorhinchus milii as well as the lampreys Petromyzon marinus and Lethenteron japonicum to learn more about the evolution of the tenascin gene family as well as the timing of the appearance of fibronectin during chordate evolution. The coelacanth has 4 tenascins that are more similar to tetrapod tenascins than are tenascins from ray-finned fishes. In contrast, only 2 tenascins were identified in the elephant shark and the Japanese lamprey L. japonicum. An RGD motif exposed to integrin binding is observed in tenascins from many, but not all, classes of chordates. Tetrapods that lack this RGD motif in tenascin-C have a similar motif in the paralog tenascin-W, suggesting the potential for some overlapping function. A predicted fibronectin with the same domain organization as the fibronectin from tetrapods is found in the sea lamprey P. marinus but not in tunicates, leading us to infer that fibronectin first appeared in vertebrates. The motifs that recognize LDV-type integrin receptors are conserved in fibronectins from a broad spectrum of vertebrates, but the RGD integrin-binding motif may have evolved in gnathostomes.

Evolution of retinoic acid receptors and retinoic acid signaling

Retinoic acid (RA) is a vitamin A-derived morphogen controlling important developmental processes in vertebrates, and more generally in chordates, including axial patterning and tissue formation and differentiation. In the embryo, endogenous RA levels are controlled by RA synthesizing and degrading enzymes and the RA signal is transduced by two retinoid receptors: the retinoic acid receptor (RAR) and the retinoid X receptor (RXR). Both RAR and RXR are members of the nuclear receptor superfamily of ligand-activated transcription factors and mainly act as heterodimers to activate the transcription of target genes in the presence of their ligand, all-trans RA. This signaling pathway was long thought to be a chordate innovation, however, recent findings of gene homologs involved in RA signaling in the genomes of a wide variety of non-chordate animals, including ambulacrarians (sea urchins and acorn worms) and lophotrochozoans (annelids and mollusks), challenged this traditional view and suggested that the RA signaling pathway might have a more ancient evolutionary origin than previously thought. In this chapter, we discuss the evolutionary history of the RA signaling pathway, and more particularly of the RARs, which might have experienced independent gene losses and duplications in different animal lineages. In sum, the available data reveal novel insights into the origin of the RA signaling pathway as well as into the evolutionary history of the RARs.

Nuclear Receptors, RXR, and the Big Bang

Isolation of genes encoding the receptors for steroids, retinoids, vitamin D, and thyroid hormone and their structural and functional analysis revealed an evolutionarily conserved template for nuclear hormone receptors. This discovery sparked identification of numerous genes encoding related proteins, termed orphan receptors. Characterization of these orphan receptors and, in particular, of the retinoid X receptor (RXR) positioned nuclear receptors at the epicenter of the "Big Bang" of molecular endocrinology. This Review provides a personal perspective on nuclear receptors and explores their integrated and coordinated signaling networks that are essential for multicellular life, highlighting the RXR heterodimer and its associated ligands and transcriptional mechanism.

Evolution of melanocortin receptors in cartilaginous fish: melanocortin receptors and the stress axis in elasmobranches

There is general agreement that the presence of five melanocortin receptor genes in tetrapods is the result of two genome duplications that occurred prior to the emergence of the gnathostomes, and at least one local gene duplication that occurred early in the radiation of the ancestral gnathostomes. Hence, it is assumed that representatives from the extant classes of gnathostomes (i.e., Chondrichthyes, Actinopterygii, Sarcopterygii) should also have five paralogous melanocortin genes. Current studies on cartilaginous fishes indicate that while there is evidence for five paralogous melanocortin receptor genes in this class, to date all five paralogs have not been detected in the genome of a single species. This mini-review will discuss the ligand selectivity properties of the melanocortin-3 receptor of the elephant shark (subclass Holocephali) and the ligand selectivity properties of the melanocortin-3 receptor, melanocortin-4 receptor, and the melanocortin-5 receptor of the dogfish (subclass Elasmobranchii). The potential relationship of these melanocortin receptors to the hypothalamus/pituitary/interrenal axis will be discussed.

Body plan convergence in the evolution of skates and rays (Chondrichthyes: Batoidea)

Skates, rays and allies (Batoidea) comprise more than half of the species diversity and much of the morphological disparity among chondrichthyan fishes, the sister group to all other jawed vertebrates. While batoids are morphologically well characterized and have an excellent fossil record, there is currently no consensus on the interrelationships of family-level taxa. Here we construct a resolved, robust and time-calibrated batoid phylogeny using mitochondrial genomes, nuclear genes, and fossils, sampling densely across taxa. Data partitioning schemes, biases in the sequence data, and the relative informativeness of each fossil are explored. The molecular phylogeny is largely congruent with morphology crownward in the tree, but the branching orders of major batoid groups are mostly novel. Body plan convergence appears to be widespread in batoids. A depressed, rounded pectoral disk supported to the snout tip by fin radials, common to skates and stingrays, is indicated to have been derived independently by each group, while the long, spiny rostrum of sawfishes similarly appears to be convergent with that of sawsharks, which are not batoids. The major extant batoid lineages are inferred to have arisen relatively rapidly from the Late Triassic into the Jurassic, with long stems followed by subsequent radiations in each group around the Cretaceous/Tertiary boundary. The fossil record indicates that batoids were affected with disproportionate severity by the end-Cretaceous extinction event.

Molecular characterization and transcriptional regulation of the sodium-dependent vitamin C transporter genes (slc23a1 and slc23a2) in a teleost fish, the Senegalese sole (Solea senegalensis)

Vitamin C (ascorbic acid, AA) is an antioxidant that acts as a free radical scavenger and cofactor for several important enzymatic reactions, thus being important for normal cellular functions, growth and development. Accumulation of AA in cells depends on two types of sodium-dependent vitamin C transporters (SVCTs), designed as SVCT1 and SVCT2. In human, they are the products of SLC23A1 and SLC23A2 genes, respectively. In the present work, the molecular cloning of the cDNAs corresponding to slc23a1 and slc23a2 in a teleost fish, the Senegalese sole (Solea senegalensis Kaup, 1858) is first described. Sequence analysis of the predicted polypeptides revealed a conserved topology with those of mammals with important motifs involved in structure and function, being also present in svct1 and svct2. Phylogenetic analyses including a range of vertebrate SVCTs suggest that both transporters are the result of an ancient gene duplication event that occurred prior to the divergence of tetrapods and teleosts, which took place 450 million years ago. Expression profiles in juvenile tissues and during larval development were analyzed using a real-time PCR approach. In juvenile fish, slc23a1 was strongly expressed in intestine, whereas slc23a2 exhibited a widespread distribution in tissues. Transcripts of both genes were detected at early developmental stages, probably representing mRNAs of maternal origin. A possible regulation by their own substrate was detected after first uptakes of AA from diet in both genes. During metamorphosis, both slc23a1 and slc23a2 were down-regulated, the former in a thyroid hormone (TH) dependent way. This pattern coincided with a significant reduction in the AA content of larvae during metamorphosis. These results are interpreted in a physiological context of general reduction in the metabolism of metamorphic larvae. Data presented here provide the first step toward a better understanding of the physiological role of SVCTs in teleost fish.

The two tryptophans of β2-microglobulin have distinct roles in function and folding and might represent two independent responses to evolutionary pressure

Background

We have recently discovered that the two tryptophans of human β2-microglobulin have distinctive roles within the structure and function of the protein. Deeply buried in the core, Trp95 is essential for folding stability, whereas Trp60, which is solvent-exposed, plays a crucial role in promoting the binding of β2-microglobulin to the heavy chain of the class I major histocompatibility complex (MHCI). We have previously shown that the thermodynamic disadvantage of having Trp60 exposed on the surface is counter-balanced by the perfect fit between it and a cavity within the MHCI heavy chain that contributes significantly to the functional stabilization of the MHCI. Therefore, based on the peculiar differences of the two tryptophans, we have analysed the evolution of β2-microglobulin with respect to these residues.

Results

Having defined the β2-microglobulin protein family, we performed multiple sequence alignments and analysed the residue conservation in homologous proteins to generate a phylogenetic tree. Our results indicate that Trp60 is highly conserved, whereas some species have a Leu in position 95; the replacement of Trp95 with Leu destabilizes β2-microglobulin by 1 kcal/mol and accelerates the kinetics of unfolding. Both thermodynamic and kinetic data fit with the crystallographic structure of the Trp95Leu variant, which shows how the hydrophobic cavity of the wild-type protein is completely occupied by Trp95, but is only half filled by Leu95.

Conclusions

We have established that the functional Trp60 has been present within the sequence of β2-microglobulin since the evolutionary appearance of proteins responsible for acquired immunity, whereas the structural Trp95 was selected and stabilized, most likely, for its capacity to fully occupy an internal cavity of the protein thereby creating a better stabilization of its folded state.

Molecular phylogeny and functional genomics of beta-galactoside alpha2,6-sialyltransferases that explain ubiquitous expression of st6gal1 gene in amniotes

Sialyltransferases are key enzymes in the biosynthesis of sialoglycoconjugates that catalyze the transfer of sialic residue from its activated form to an oligosaccharidic acceptor. β-Galactoside α2,6-sialyltransferases ST6Gal I and ST6Gal II are the two unique members of the ST6Gal family described in higher vertebrates. The availability of genome sequences enabled the identification of more distantly related invertebrates' st6gal gene sequences and allowed us to propose a scenario of their evolution. Using a phylogenomic approach, we present further evidence of an accelerated evolution of the st6gal1 genes both in their genomic regulatory sequences and in their coding sequence in reptiles, birds, and mammals known as amniotes, whereas st6gal2 genes conserve an ancestral profile of expression throughout vertebrate evolution.

Molecular characterization and gene expression of thyrotropin receptor (TSHR) and a truncated TSHR-like in Senegalese sole

Thyroid hormones (THs) play a key role in larval development, growth and metamorphosis in flatfish. Their synthesis is tightly regulated by the hypothalamic-pituitary-thyroid axis. Thyroid-stimulating hormone receptor (TSHR) is a key protein in the control of thyroid function stimulating TH synthesis after binding its ligand, the thyrotropin. In teleost fish, numerous reports have associated the TSHR with gametogenesis. However, little information about its role during larval development is available. In this study, we report the cloning of two different cDNAs with high similarity to TSHR. Phylogenetic analysis clustered both cDNAs separately. One of them (referred to TSHR) grouped with TSHR orthologs in tetrapods and teleost fish and possessed the three typical conserved domains and regulatory motifs. The second receptor (referred to as TSHRtr-like) represented a novel truncated cDNA bearing the extracellular and part of the transmembrane domain. TSHRtr-like orthologs were only found in teleosts, which suggests that it could have appeared after fish-specific 3R genome duplication. Expression profiles of both genes are analyzed in juvenile tissues and during larval development using a real-time PCR approach. In juvenile fish, TSHR and TSHRtr-like are expressed ubiquitously although transcript levels varied between organs. In both cases, the highest mRNAs levels are detected in brain. During larval development, both genes are expressed to a high level during the first stages (2-3days after hatching) reducing progressively their abundance in the whole larvae during metamorphosis. This reduction in mRNA abundance is more accentuated for the TSHRtr-like gene. To evaluate the possible regulation of both receptors by T4 during sole metamorphosis, larvae are exposed to the goitrogen thiourea (TU). Only TSHRtr-like modifies its expression, increasing its transcripts at 11days after treatment. Moreover, adding exogenous T4 hormone to TU-treated larvae restores the TSHRtr-like steady-state levels similar to the untreated control. Overall, these results demonstrate the existence of two thyrotropin receptors differentially regulated by THs in teleosts.

Evolution of the nitric oxide synthase family in metazoans

Nitric oxide (NO) is essential to many physiological functions and operates in several signaling pathways. It is not understood how and when the different isoforms of nitric oxide synthase (NOS), the enzyme responsible for NO production, evolved in metazoans. This study investigates the number and structure of metazoan NOS enzymes by genome data mining and direct cloning of Nos genes from the lamprey. In total, 181 NOS proteins are analyzed from 33 invertebrate and 63 vertebrate species. Comparisons among protein and gene structures, combined with phylogenetic and syntenic studies, provide novel insights into how NOS isoforms arose and diverged. Protein domains and gene organization--that is, intron positions and phases--of animal NOS are remarkably conserved across all lineages, even in fast-evolving species. Phylogenetic and syntenic analyses support the view that a proto-NOS isoform was recurrently duplicated in different lineages, acquiring new structural configurations through gains and losses of protein motifs. We propose that in vertebrates a first duplication took place after the agnathan-gnathostome split followed by a paralog loss. A second duplication occurred during early tetrapod evolution, giving rise to the three isoforms--I, II, and III--in current mammals. Overall, NOS family evolution was the result of multiple gene and genome duplication events together with changes in protein architecture.

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.

Stage-specific gene expression during spermatogenesis in the dogfish (Scyliorhinus canicula)

In the dogfish testis, the cystic arrangement and polarization of germ cell stages make it possible to observe all stages of spermatogenesis in a single transverse section. By taking advantage of the zonation of this organ, we have used suppressive subtractive libraries construction, real-time PCR, and in situ hybridization to identify 32 dogfish genes showing differential expressions during spermatogenesis. These include homologs of genes already known to be expressed in the vertebrate testis, but found here to be specifically expressed either in pre-meiotic and/or meiotic zones (ribosomal protein S8, high-mobility group box 3, ubiquitin carboxyl-terminal esterase L3, 20beta-hydroxysteroid dehydrogenase, or cyclophilin B) or in post-meiotic zone (speriolin, Soggy, zinc finger protein 474, calreticulin, or phospholipase c-zeta). We also report, for the first time, testis-specific expression patterns for dogfish genes coding for A-kinase anchor protein 5, ring finger protein 152, or F-box only protein 7. Finally, the study highlights the differential expression of new sequences whose identity remains to be assessed. This study provides the first molecular characterization of spermatogenesis in a chondrichthyan, a key species to gain insight into the evolution of this process in gnathostomes.

Nuclear receptors CAR and PXR: Molecular, functional, and biomedical aspects

Nuclear receptors (NRs) are ligand-activated transcription factors sharing a common evolutionary history and having similar sequence features at the protein level. Selective ligand(s) for some NRs is not known, therefore these NRs have been named "orphan receptors". Whenever ligands have been recognized for any of the orphan receptor, it has been categorized and grouped as "adopted" orphan receptor. This group includes the constitutive androstane receptor (CAR) and the pregnane X receptor (PXR). They function as sensors of toxic byproducts derived from endogenous metabolites and of exogenous chemicals, in order to enhance their elimination. This unique function of CAR and PXR sets them apart from the steroid hormone receptors. The broad response profile has established that CAR and PXR are xenobiotic sensors that coordinately regulate xenobiotic clearance in the liver and intestine via induction of genes involved in drug and xenobiotic metabolism. In the past few years, research has revealed new and mostly unsuspected roles for CAR and PXR in modulating hormone, lipid, and energy homeostasis as well as cancer and liver steatosis. The purpose of this review is to highlight the structural and molecular bases of CAR and PXR impact on human health, providing information on mechanisms through which diet, chemical exposure, and environment ultimately impact health and disease.

The amphioxus genome and the evolution of the chordate karyotype

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Molecular characterization, phylogeny, and expression of c-type and g-type lysozymes in brill (Scophthalmus rhombus)

Lysozymes are key proteins of the innate immune system against bacterial infections. In this study we report the molecular cloning and characterization of the c-type and g-type lysozymes in brill (Scophthalmus rhombus). Catalytic and other conserved residues required for functionality were identified. Phylogenetic analysis revealed distinct evolutionary histories for each lysozyme type. Expression profiles of both lysozyme genes were studied in juvenile tissues using a real-time PCR approach. c-Type lysozyme was expressed mainly in stomach and liver, whereas the g-type was detected in all tissues with highest mRNA levels observed in the spleen. Induction experiments revealed that g-type transcripts increased significantly in head kidney after lipopolysaccharide (25- and 23-fold at 12 and 24h, respectively) and Photobacterium damselae subsp. piscicida (17-fold at 24h) treatments. In contrast, no induction was observed for c-type lysozyme. All these data suggest that g-type lysozyme is involved in the response against bacterial infections, whereas c-type lysozyme may also play a role in digestion.

Phylogenomic analyses of KCNA gene clusters in vertebrates: why do gene clusters stay intact?

Background

Gene clusters are of interest for the understanding of genome evolution since they provide insight in large-scale duplications events as well as patterns of individual gene losses. Vertebrates tend to have multiple copies of gene clusters that typically are only single clusters or are not present at all in genomes of invertebrates. We investigated the genomic architecture and conserved non-coding sequences of vertebrate KCNA gene clusters. KCNA genes encode shaker-related voltage-gated potassium channels and are arranged in two three-gene clusters in tetrapods. Teleost fish are found to possess four clusters. The two tetrapod KNCA clusters are of approximately the same age as the Hox gene clusters that arose through duplications early in vertebrate evolution. For some genes, their conserved retention and arrangement in clusters are thought to be related to regulatory elements in the intergenic regions, which might prevent rearrangements and gene loss. Interestingly, this hypothesis does not appear to apply to the KCNA clusters, as too few conserved putative regulatory elements are retained.

Results

We obtained KCNA coding sequences from basal ray-finned fishes (sturgeon, gar, bowfin) and confirmed that the duplication of these genes is specific to teleosts and therefore consistent with the fish-specific genome duplication (FSGD). Phylogenetic analyses of the genes suggest a basal position of the only intron containing KCNA gene in vertebrates (KCNA7). Sistergroup relationships of KCNA1/2 and KCNA3/6 support that a large-scale duplication gave rise to the two clusters found in the genome of tetrapods. We analyzed the intergenic regions of KCNA clusters in vertebrates and found that there are only a few conserved sequences shared between tetrapods and teleosts or between paralogous clusters. The orthologous teleost clusters, however, show sequence conservation in these regions.

Conclusion

The lack of overall conserved sequences in intergenic regions suggests that there are either other processes than regulatory evolution leading to cluster conservation or that the ancestral regulatory relationships among genes in KCNA clusters have been changed together with their regulatory sites.

Functional evolution of the pregnane X receptor

The pregnane X receptor (PXR; NR1I2) is a nuclear hormone receptor (NR) that transcriptionally regulates genes encoding transporters and drug-metabolising enzymes in the liver and intestine. PXR activation leads to enhanced metabolism and elimination of xenobiotics and endogenous compounds such as hormones and bile salts. Relative to other vertebrate NRs, PXR has the broadest specificity for ligand activators by virtue of a large, flexible ligand-binding cavity. In addition, PXR has the most extensive sequence diversity across vertebrate species in the ligand-binding domain of any NR, with significant pharmacological differences between human and rodent PXRs, and especially marked divergence between mammalian and nonmammalian PXRs. The unusual properties of PXR complicate the use of in silico and animal models to predict in vivo human PXR pharmacology. Research into the evolutionary history of the PXR gene has also provided insight into the function of PXR in humans and other animals.

Comparison of the expression of medaka (Oryzias latipes) pitx genes with other vertebrates shows high conservation and a case of functional shuffling in the pituitary

With the availability of an increasing number of whole genome sequences in chordates, exhaustive comparisons of multigene families become feasible. Relationships of orthology/paralogy can not only be inferred from sequence similarity but also by comparing synteny conservation on chromosomes. More accurate scenarios for gene and expression domain gain or loss can now be proposed. Here, we take benefit from the recent release of the medaka (Oryzias latipes) genome to analyse the orthology relationships and expression patterns of the three different sub-families of the pitx homeobox genes belonging to the paired class. They are involved in a wide variety of developmental processes and have pleiotropic expression patterns, especially in the case of the pitx2 sub-family. The emerging picture is a strong conservation of expression domains, suggesting that most functions have been present in the common ancestor of actinopterygians and sarcopterygians. Almost all pitx genes are expressed in anterior placodes in all species studied so far, including medaka. It has previously been shown that in mammals, pitx1 and 2 are expressed in the pituitary. Interestingly we demonstrate here that only pitx3 is expressed in medaka pituitary. It will be interesting to analyze what are the corresponding changes in the regulatory elements of pitx genes.

The plasticity of immunoglobulin gene systems in evolution

The mechanism of recombination-activating gene (RAG)-mediated rearrangement exists in all jawed vertebrates, but the organization and structure of immunoglobulin (Ig) genes, as they differ in fish and among fish species, reveal their capability for rapid evolution. In systems where there can exist 100 Ig loci, exon restructuring and sequence changes of the constant regions led to divergence of effector functions. Recombination among these loci created hybrid genes, the strangest of which encode variable (V) regions that function as part of secreted molecules and, as the result of an ancient translocation, are also grafted onto the T-cell receptor. Genomic changes in V-gene structure, created by RAG recombinase acting on germline recombination signal sequences, led variously to the generation of fixed receptor specificities, pseudogene templates for gene conversion, and ultimately to Ig sequences that evolved away from Ig function. The presence of so many Ig loci in fishes raises interesting questions not only as to how their regulation is achieved but also how successive whole-locus duplications are accommodated by a system whose function in other vertebrates is based on clonal antigen receptor expression.

Evolution of axis specification mechanisms in jawed vertebrates: insights from a chondrichthyan

The genetic mechanisms that control the establishment of early polarities and their link with embryonic axis specification and patterning seem to substantially diverge across vertebrates. In amphibians and teleosts, the establishment of an early dorso-ventral polarity determines both the site of axis formation and its rostro-caudal orientation. In contrast, amniotes retain a considerable plasticity for their site of axis formation until blastula stages and rely on signals secreted by extraembryonic tissues, which have no clear equivalents in the former, for the establishment of their rostro-caudal pattern. The rationale for these differences remains unknown. Through detailed expression analyses of key development genes in a chondrichthyan, the dogfish Scyliorhinus canicula, we have reconstructed the ancestral pattern of axis specification in jawed vertebrates. We show that the dogfish displays compelling similarities with amniotes at blastula and early gastrula stages, including the presence of clear homologs of the hypoblast and extraembryonic ectoderm. In the ancestral state, these territories are specified at opposite poles of an early axis of bilateral symmetry, homologous to the dorso-ventral axis of amphibians or teleosts, and aligned with the later forming embryonic axis, from head to tail. Comparisons with amniotes suggest that a dorsal expansion of extraembryonic ectoderm, resulting in an apparently radial symmetry at late blastula stages, has taken place in their lineage. The synthesis of these results with those of functional analyses in model organisms supports an evolutionary link between the dorso-ventral polarity of amphibians and teleosts and the embryonic-extraembryonic organisation of amniotes. It leads to a general model of axis specification in gnathostomes, which provides a comparative framework for a reassessment of conservations both among vertebrates and with more distant metazoans.

Unexpected diversity of aryl hydrocarbon receptors in non-mammalian vertebrates: insights from comparative genomics

Ligand-activated receptors are well-known targets of environmental chemicals that disrupt endocrine signaling. Genomic approaches are providing new opportunities to understand the comparative biology and molecular evolution of these receptors. One example of this is the aryl hydrocarbon receptor (AHR), a basic-helix-loop-helix (bHLH)-Per-Arnt-Sim (PAS) transcription factor through which planar aromatic hydrocarbons cause altered gene expression and toxicity. In contrast to humans and other mammals, which possess a single AHR, teleosts such as the Atlantic killifish (Fundulus heteroclitus) have at least two AHRs (AHR1 and AHR2). Analysis of sequenced genomes has revealed additional, unexpected AHR diversity in non-mammalian vertebrates, including the chicken Gallus gallus (three predicted AHR genes), bony fishes such as the pufferfish Takifugu (formerly Fugu) rubripes (five AHR genes) and zebrafish Danio rerio (three AHR genes), and cartilaginous fishes such as the spiny dogfish Squalus acanthias (three AHR genes). In contrast, invertebrates appear to possess single AHRs that do not bind typical ligands of vertebrate AHRs. We suggest that AHR diversity in vertebrates arose through both gene and whole-genome duplications combined with lineage-specific gene loss, and that sensitivity to the developmental toxicity of planar aromatic hydrocarbons may have had its origin in the evolution of the ligand-binding capacity of the AHR in the chordate lineage. Comparative molecular and genomic studies are providing new insights into AHR diversity and function in non-mammalian species, revealing additional complexity in mechanisms by which environmental chemicals interfere with receptor-dependent signaling.

En bloc duplications, mutation rates, and densities of amino acid changes clarify the evolution of vertebrate alpha-1,3/4-fucosyltransferases

Numerous vertebrates have four alpha-1,3/4-fucosyltransferase genes (FUT9, FUT7, FUT4, and FUT Lewis) belonging to the same family. Until now, studies on the evolution of this family have mainly focused on Lewis genes but how the other alpha-1,3/4-fucosyltransferases have emerged from a common ancestor is not well known. In order to define the respective roles of duplications and mutations, we have compared amino acid sequences representative of bony fish (Takifugu rubripes), amphibians (Xenopus laevis), birds (Gallus gallus), and mammals (Bos taurus). The FUT tree has two fundamental branches, each split into two subfamilies. We found evidence for two duplication events, dated around 710-760 Myr and 590-640 Myr, respectively, compatible with the hypothesis of two rounds of whole genome duplications in chordate genomes, before the emergence of bony vertebrates. Based on the Homo sapiens (human) physical map, we identified blocks of paralogues belonging to regions of FUT9 (6q16), FUT4 (11q21), FUT7 (9q34), and FUT Lewis (19p13) and to a region on HSA1p that is devoid of any FUT. In zebrafish (Danio rerio), an orthologue region of HSA1 harbors an FUT9 specific to bony fish, showing that duplications are not restricted to a single FUT gene but involve blocks of paralogues. In addition, sets of genes within each block clarify the order of duplication events and, as a result, the order of alpha-1,3/4-fucosyltransferase gene emergence. We have also determined the mutation rates and the density of amino acid changes along protein sequences in each alpha-1,3/4-fucosyltransferase subfamily during the main vertebrate transitions. After the emergence of tetrapods, the mutation rate of FUT9 decreased dramatically, suggesting the early acquisition of a crucial fucosyltransferase activity in the first stages of development. The FUT7 mutation rate, which in tetrapod ancestors is about half that in amniote ancestors, may be related to the role of this gene in immune systems. In contrast to other subfamilies, we found a constant mutation rate in FUT Lewis and a rather homogeneous amino acid density change, independently of the vertebrate transition, suggesting that hitherto Lewis epitopes have dispensable functions.

Neofunctionalization in vertebrates: the example of retinoic acid receptors

Understanding the role of gene duplications in establishing vertebrate innovations is one of the main challenges of Evo-Devo (evolution of development) studies. Data on evolutionary changes in gene expression (i.e., evolution of transcription factor-cis-regulatory elements relationships) tell only part of the story; protein function, best studied by biochemical and functional assays, can also change. In this study, we have investigated how gene duplication has affected both the expression and the ligand-binding specificity of retinoic acid receptors (RARs), which play a major role in chordate embryonic development. Mammals have three paralogous RAR genes—RARα, β, and γ—which resulted from genome duplications at the origin of vertebrates. By using pharmacological ligands selective for specific paralogues, we have studied the ligand-binding capacities of RARs from diverse chordates species. We have found that RARβ-like binding selectivity is a synapomorphy of all chordate RARs, including a reconstructed synthetic RAR representing the receptor present in the ancestor of chordates. Moreover, comparison of expression patterns of the cephalochordate amphioxus and the vertebrates suggests that, of all the RARs, RARβ expression has remained most similar to that of the ancestral RAR. On the basis of these results together, we suggest that while RARβ kept the ancestral RAR role, RARα and RARγ diverged both in ligand-binding capacity and in expression patterns. We thus suggest that neofunctionalization occurred at both the expression and the functional levels to shape RAR roles during development in vertebrates.

In eukaryotic organisms, each gene is a stretch of DNA composed of control regions that bind transcription factors and coding regions that transcribe the mRNA that is later translated into proteins. At the molecular level, changes in control regions can affect the time and place at which a protein is synthesized, whereas changes in the coding region can alter the protein's function. Retinoic acid receptors (RARs) are chordate-specific transcription factors which, upon binding the natural morphogen retinoic acid, bind to and activate transcription from target genes. Here, the authors show how the ligand specificity of RARs has changed during vertebrate evolution in parallel with changes in expression. Through functional characterization of the RARs from several vertebrates, the chordate amphioxus, and the reconstructed ancestral RAR sequence, the authors show that of the three vertebrate RARs, RARβ has retained the ancestral characteristics in terms of both function and expression, while RARα and γ have evolved by acquiring new functions, both new binding specificity and new expression patterns. Thus both types of evolution have been important in the diversification of vertebrate RARs.

Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates

Background

Evolution of the deuterostome lineage was accompanied by an increase in systematic complexity especially with regard to highly specialized tissues and organs. Based on the observation of an increased number of paralogous genes in vertebrates compared with invertebrates, two entire genome duplications (2R) were proposed during the early evolution of vertebrates. Most glycolytic enzymes occur as several copies in vertebrate genomes, which are specifically expressed in certain tissues. Therefore, the glycolytic pathway is particularly suitable for testing theories of the involvement of gene/genome duplications in enzyme evolution.

Results

We assembled datasets from genomic databases of at least nine vertebrate species and at least three outgroups (one deuterostome and two protostomes), and used maximum likelihood and Bayesian methods to construct phylogenies of the 10 enzymes of the glycolytic pathway. Through this approach, we intended to gain insights into the vertebrate specific evolution of enzymes of the glycolytic pathway. Many of the obtained gene trees generally reflect the history of two rounds of duplication during vertebrate evolution, and were in agreement with the hypothesis of an additional duplication event within the lineage of teleost fish. The retention of paralogs differed greatly between genes, and no direct link to the multimeric structure of the active enzyme was found.

Conclusion

The glycolytic pathway has subsequently evolved by gene duplication and divergence of each constituent enzyme with taxon-specific individual gene losses or lineage-specific duplications. The tissue-specific expression might have led to an increased retention for some genes since paralogs can subdivide the ancestral expression domain or find new functions, which are not necessarily related to the original function.

Exploring developmental, functional, and evolutionary aspects of amphioxus sensory cells

Amphioxus has neither elaborated brains nor definitive sensory organs, so that the two may have evolved in a mutually affecting manner and given rise to the forms seen in extant vertebrates. Clarifying the developmental and functional aspects of the amphioxus sensory system is thus pivotal for inferring the early evolution of vertebrates. Morphological studies have identified and classified amphioxus sensory cells; however, it is completely unknown whether the morphological classification makes sense in functional and evolutionary terms. Molecular markers, such as gene expression, are therefore indispensable for investigating the developmental and functional aspects of amphioxus sensory cells. This article reviews recent molecular studies on amphioxus sensory cells. Increasing evidence shows that the non-neural ectoderm of amphioxus can be subdivided into molecularly distinct subdomains by the combinatorial code of developmental cues involving the RA-dependent Hox code, suggesting that amphioxus epithelial sensory cells developed along positional information. This study focuses particularly on research involving the molecular phylogeny and expression of the seven-transmembrane, G protein-coupled receptor (GPCR) genes and discusses the usefulness of this information for characterizing the sensory cells of amphioxus.

Evolution and function of the NR1I nuclear hormone receptor subfamily (VDR, PXR, and CAR) with respect to metabolism of xenobiotics and endogenous compounds

The NR1I subfamily of nuclear hormone receptors includes the 1,25-(OH)(2)-vitamin D(3) receptor (VDR; NR1I1), pregnane X receptor (PXR; NR1I2), and constitutive androstane receptor (CAR; NR1I3). PXR and VDR are found in diverse vertebrates from fish to mammals while CAR is restricted to mammals. Current evidence suggests that the CAR gene arose from a duplication of an ancestral PXR gene, and that PXR and VDR arose from duplication of an ancestral gene, represented now by a single gene in the invertebrate Ciona intestinalis. Aside from the high-affinity effects of 1,25-(OH)(2)-vitamin D(3) on VDRs, the NR1I subfamily members are functionally united by the ability to bind potentially toxic endogenous compounds with low affinity and initiate changes in gene expression that lead to enhanced metabolism and elimination (e.g., induction of cytochrome P450 3A4 expression in humans). The detoxification role of VDR seems limited to sensing high concentrations of certain toxic bile salts, such as lithocholic acid, whereas PXR and CAR have the ability to recognize structurally diverse compounds. PXR and CAR show the highest degree of cross-species variation in the ligand-binding domain of the entire vertebrate nuclear hormone receptor superfamily, suggesting adaptation to species-specific ligands. This review examines the insights that phylogenetic and experimental studies provide into the function of VDR, PXR, and CAR, and how the functions of these receptors have expanded to evolutionary advantage in humans and other animals.

Phylogenomic analysis of vertebrate thrombospondins reveals fish-specific paralogues, ancestral gene relationships and a tetrapod innovation

BACKGROUND

Thrombospondins (TSPs) are evolutionarily-conserved, extracellular, calcium-binding glycoproteins with important roles in cell-extracellular matrix interactions, angiogenesis, synaptogenesis and connective tissue organisation. Five TSPs, designated TSP-1 through TSP-5, are encoded in the human genome. All but one have known roles in acquired or inherited human diseases. To further understand the roles of TSPs in human physiology and pathology, it would be advantageous to extend the repertoire of relevant vertebrate models. In general the zebrafish is proving an excellent model organism for vertebrate biology, therefore we set out to evaluate the status of TSPs in zebrafish and two species of pufferfish.

RESULTS

We identified by bioinformatics that three fish species encode larger numbers of TSPs than vertebrates, yet all these sequences group as homologues of TSP-1 to -4. By phylogenomic analysis of neighboring genes, we uncovered that, in fish, a TSP-4-like sequence is encoded from the gene corresponding to the tetrapod TSP-5 gene. Thus, all TSP genes show conservation of synteny between fish and tetrapods. In the human genome, the TSP-1, TSP-3, TSP-4 and TSP-5 genes lie within paralogous regions that provide insight into the ancestral genomic context of vertebrate TSPs.

CONCLUSION

A new model for TSP evolution in vertebrates is presented. The TSP-5 protein sequence has evolved rapidly from a TSP-4-like sequence as an innovation in the tetrapod lineage. TSP biology in fish is complicated by the presence of additional lineage- and species-specific TSP paralogues. These novel results give deeper insight into the evolution of TSPs in vertebrates and open new directions for understanding the physiological and pathological roles of TSP-4 and TSP-5 in humans.

Duplication events and the evolution of segmental identity

Duplication of genes, genomes, or morphological structures (or some combination of these) has long been thought to facilitate evolutionary change. Here we focus on studies of the teleost fishes to consider the conceptual similarities in the evolutionary potential of these three different kinds of duplication events. We review recent data that have confirmed the occurrence of a whole-genome duplication event in the ray-finned fish lineage, and discuss whether this event may have fuelled the radiation of teleost fishes. We then consider the fates of individual duplicated genes, from both a theoretical and an experimental viewpoint, focusing on our studies of teleost Hox genes and their functions in patterning the segmented hindbrain. Finally, we consider the duplication of morphological structures, once again drawing on our experimental studies of the hindbrain, which have revealed that experimentally induced duplicated neurons can produce functionally redundant neural circuits. We posit that the availability of duplicated material, independent of its nature, can lead to functional redundancy, which in turn enables evolutionary change.

Organisation of the lamprey (Lampetra fluviatilis) embryonic brain: insights from LIM-homeodomain, Pax and hedgehog genes

To investigate the embryonic development of the central nervous system of the lamprey Lampetra fluviatilis, we have isolated and analysed the expression patterns of members of the LIM-homeodomain, Pax, Hedgehog and Nkx2.1 families. Using degenerate RT-PCR, single representatives of Lhx1/Lhx5, Lhx2/Lhx9, Pax3/Pax7 and Hedgehog families could be isolated in L. fluviatilis. Expression analysis revealed that the lamprey forebrain presents a clear neuromeric pattern. We describe the existence of 4 embryonic diencephalic prosomeres whose boundaries can be identified by the combined and relative expressions of LfPax37, LfLhx15 and LfLhx29. This suggests that the embryonic lamprey and gnathostome forebrain are patterned in a highly similar manner. Moreover, analysis of the LfHh gene, which is expressed in the hypothalamus, zona limitans intrathalamica and floor plate, reveals the possible molecular origin of this neuromeric brain pattern. By contrast, LfHh and LfNkx2.1 expressions suggest major differences in patterning mechanisms of the ventral telencephalon when compared to gnathostomes. In summary, our findings highlight a neuromeric organisation of the embryonic agnathan forebrain and point to the possible origin of this organisation, which is thus a truly vertebrate character. They also suggest that Hh/Shh midline signalling might act as a driving force for forebrain evolution.

Molecular Evolution of the Ankyrin Gene Family

Ankyrins are membrane adaptor molecules that play important roles in coupling integral membrane proteins to the spectrin-based cytoskeleton network. Human mutations of ankyrin genes lead to severe genetic diseases such as fatal cardiac arrhythmias and hereditary spherocytosis. To elucidate the evolutionary history of ankyrins, we have identified novel ankyrin sequences in insect, fish, frog, chicken, dog, and chimpanzee genomes and explored the phylogenetic relationships of the ankyrin gene family. Our data demonstrate that duplication of ankyrin genes occurred at two different stages. The first duplication resulted from an independent evolution event specific in Arthropoda after its divergence from Chordata. Following the separation from Urochordata, expansion of ankyrins in vertebrates involved ancestral genome duplications. We did not find evidence of coordinated arrangements of gene families of ankyrin-associated membrane proteins on paralogous chromosomes. In addition, evolution of the 24 ANK-repeats strikingly correlated with the exon boundary sites of ankyrin genes, which might have occurred before its duplication in vertebrates. Such correlation is speculated to bring functional diversity and complexity. Moreover, based on the phylogenetic analysis of the ANK-repeat domain, we put forward a novel model for the putative primordial ankyrin that contains the fourth six-ANK-repeat subdomain and the spectrin-binding domain. These findings will provide guides for future studies concerning structure, function, evolutionary origins of ankyrins, and possibly other cytoskeletal proteins.