Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Evolution of Vertebrate Hormones and Their Receptors: Insights from Non-Osteichthyan Genomes

Homeostatic control and reproductive functions of humans are regulated at the molecular levels largely by peptide hormones secreted from endocrine and/or neuroendocrine cells in the central nervous system and peripheral organs. Homologs of those hormones and their receptors function similarly in many vertebrate species distantly related to humans, but the evolutionary history of the endocrine system involving those factors has been obscured by the scarcity of genome DNA sequence information of some taxa that potentially contain their orthologs. Focusing on non-osteichthyan vertebrates, namely jawless and cartilaginous fishes, this article illustrates how investigating genome sequence information assists our understanding of the diversification of vertebrate gene repertoires in four broad themes: (a) the presence or absence of genes, (b) multiplication and maintenance of paralogs, (c) differential fates of duplicated paralogs, and (d) the evolutionary timing of gene origins.

Evolution of calcitonin/calcitonin gene-related peptide family in chordates: Identification of CT/CGRP family peptides in cartilaginous fish genome

The calcitonin (CT)/CT gene-related peptide (CGRP) family is a peptide gene family that is widely found in bilaterians. CT, CGRP, adrenomedullin (AM), amylin (AMY), and CT receptor-stimulating peptide (CRSP) are members of the CT/CGRP family. In mammals, CT is involved in calcium homeostasis, while CGRP and AM primarily function in vasodilation. AMY and CRSP are associated with anorectic effects. Diversification of the molecular features and physiological functions of the CT/CGRP family in vertebrate lineages have been extensively reported. However, the origin and diversification mechanisms of the vertebrate CT/CGRP family of peptides remain unclear. In this review, the molecular characteristics of CT/CGRP family peptides and their receptors, along with their major physiological functions in mammals and teleosts, are introduced. Furthermore, novel candidates of the CT/CGRP family in cartilaginous fish are presented based on genomic information. The CT/CGRP family peptides and receptors in urochordates and cephalochordates, which are closely related to vertebrates, are also described. Finally, a putative evolutionary scenario of the CT/CGRP family peptides and receptors in chordates is discussed.

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

Evolutionary History of TOPIIA Topoisomerases in Animals

TOPIIA topoisomerases are required for the regulation of DNA topology by DNA cleavage and re-ligation and are important targets of antibiotic and anticancer agents. Humans possess two TOPIIA paralogue genes (TOP2A and TOP2B) with high sequence and structural similarity but distinct cellular functions. Despite their functional and clinical relevance, the evolutionary history of TOPIIA is still poorly understood. Here we show that TOPIIA is highly conserved in Metazoa. We also found that TOPIIA paralogues from jawed and jawless vertebrates had different origins related with tetraploidization events. After duplication, TOP2B evolved under a stronger purifying selection than TOP2A, perhaps promoted by the more specialized role of TOP2B in postmitotic cells. We also detected genetic signatures of positive selection in the highly variable C-terminal domain (CTD), possibly associated with adaptation to cellular interactions. By comparing TOPIIA from modern and archaic humans, we found two amino acid substitutions in the TOP2A CTD, suggesting that TOP2A may have contributed to the evolution of present-day humans, as proposed for other cell cycle-related genes. Finally, we identified six residues conferring resistance to chemotherapy differing between TOP2A and TOP2B. These six residues could be targets for the development of TOP2A-specific inhibitors that would avoid the side effects caused by inhibiting TOP2B. Altogether, our findings clarify the origin, diversification and selection pressures governing the evolution of animal TOPIIA.

Ferroptosis: regulation by competition between NRF2 and BACH1 and propagation of the death signal

Ferroptosis is triggered by a chain of intracellular labile iron-dependent peroxidation of cell membrane phospholipids. Ferroptosis is important not only as a cause of ischaemic and neurodegenerative diseases but also as a mechanism of cancer suppression, and a better understanding of its regulatory mechanism is required. It has become clear that ferroptosis is finely controlled by two oxidative stress-responsive transcription factors, NRF2 (NF-E2-related factor 2) and BACH1 (BTB and CNC homology 1). NRF2 and BACH1 inhibit and promote ferroptosis, respectively, by activating or suppressing the expression of genes in the major regulatory pathways of ferroptosis: intracellular labile iron metabolism, the GSH (glutathione) -GPX4 (glutathione peroxidase 4) pathway and the FSP1 (ferroptosis suppressor protein 1)-CoQ (coenzyme Q) pathway. In addition to this, NRF2 and BACH1 control ferroptosis through the regulation of lipid metabolism and cell differentiation. This multifaceted regulation of ferroptosis by NRF2 and BACH1 is considered to have been acquired during the evolution of multicellular organisms, allowing the utilization of ferroptosis for maintaining homeostasis, including cancer suppression. In terms of cell-cell interaction, it has been revealed that ferroptosis has the property of propagating to surrounding cells along with lipid peroxidation. The regulation of ferroptosis by NRF2 and BACH1 and the propagation phenomenon could be used to realize anticancer cell therapy in the future. In this review, these points will be summarized and discussed.

Alternative Pre-mRNA Splicing of the Mu Opioid Receptor Gene, OPRM1: Insight into Complex Mu Opioid Actions

Most opioid analgesics used clinically, including morphine and fentanyl, as well as the recreational drug heroin, act primarily through the mu opioid receptor, a class A Rhodopsin-like G protein-coupled receptor (GPCR). The single-copy mu opioid receptor gene, OPRM1, undergoes extensive alternative splicing, creating multiple splice variants or isoforms via a variety of alternative splicing events. These OPRM1 splice variants can be categorized into three major types based on the receptor structure: (1) full-length 7 transmembrane (TM) C-terminal variants; (2) truncated 6TM variants; and (3) single TM variants. Increasing evidence suggests that these OPRM1 splice variants are pharmacologically important in mediating the distinct actions of various mu opioids. More importantly, the OPRM1 variants can be targeted for development of novel opioid analgesics that are potent against multiple types of pain, but devoid of many side-effects associated with traditional opiates. In this review, we provide an overview of OPRM1 alternative splicing and its functional relevance in opioid pharmacology.

Molecular Phylogenetic Analysis of the AIG Family in Vertebrates

Androgen-inducible genes (AIGs), which can be regulated by androgen level, constitute a group of genes characterized by the presence of the AIG/FAR-17a domain in its protein sequence. Previous studies on AIGs demonstrated that one member of the gene family, AIG1, is involved in many biological processes in cancer cell lines and that ADTRP is associated with cardiovascular diseases. It has been shown that the numbers of AIG paralogs in humans, mice, and zebrafish are 2, 2, and 3, respectively, indicating possible gene duplication events during vertebrate evolution. Therefore, classifying subgroups of AIGs and identifying the homologs of each AIG member are important to characterize this novel gene family further. In this study, vertebrate AIGs were phylogenetically grouped into three major clades, ADTRP, AIG1, and AIG-L, with AIG-L also evident in an outgroup consisting of invertebrsate species. In this case, AIG-L, as the ancestral AIG, gave rise to ADTRP and AIG1 after two rounds of whole-genome duplications during vertebrate evolution. Then, the AIG family, which was exposed to purifying forces during evolution, lost or gained some of its members in some species. For example, in eutherians, Neognathae, and Percomorphaceae, AIG-L was lost; in contrast, Salmonidae and Cyprinidae acquired additional AIG copies. In conclusion, this study provides a comprehensive molecular phylogenetic analysis of vertebrate AIGs, which can be employed for future functional characterization of AIGs.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

Immunohistochemical detection of hepatocyte nuclear factor-4α in vertebrates

Hepatocyte nuclear factor-4α (HNF4α) presents in multiple isoforms generated using alternative promoter (P1 and P2) and splicing. Neither conservation of tissue distribution of HNF4α isoforms, nor presence of alternative promoter usage is known. In this study, to detect the expression of HNF4α in some species of animals, we have applied monoclonal antibodies against P1 (K9218) and P2 (H6939) promoter-driven and P1/P2 promoter-driven H1415 HNF4α for immunohistochemistry and western blot analysis. Antibody K9218 was observed in the hepatocytes, proximal tubules of the kidney, and epithelial cells in the mucosa of the small intestine and colon of rats, chicken, and tortoise, whereas antibody H6939 signal were detected in the stomach, pancreas, bile duct, and pancreatic duct of human and rats. The signal for antibody K9218 was recognized in tissues of a wide range of mammals, bird, reptile, amphibian, and fish as well. Antibody H1415 displayed a positive reaction in hepatocytes and intestinal epithelial cells in chicken and tortoise, whereas the bile duct, mucosal epithelial cells in the stomach, or pancreas in these animals were negative. Western blotting showed the binding of the antibody with HNF4α protein from each animal. The sequence of human HNF4α was 100% identical to murine and rat HNF4α, 88.9% to chicken, 77.8% to Xenopus HNF4α, and 81.5% to medaka. However, the specific part of human and invertebrate Drosophila HNF4 shares only 14.8% sequence identity. This antibody is useful for detecting HNF4α isoforms in a wide range of vertebrates, and suggests many insights into animal evolution.

Function and Evolution of Nuclear Receptors in Environmental-Dependent Postembryonic Development

Nuclear receptors (NRs) fulfill key roles in the coordination of postembryonal developmental transitions in animal species. They control the metamorphosis and sexual maturation in virtually all animals and by that the two main environmental-dependent developmental decision points. Sexual maturation and metamorphosis are controlled by steroid receptors and thyroid receptors, respectively in vertebrates, while both processes are orchestrated by the ecdysone receptor (EcR) in insects. The regulation of these processes depends on environmental factors like nutrition, temperature, or photoperiods and by that NRs form evolutionary conserved mediators of phenotypic plasticity. While the mechanism of action for metamorphosis and sexual maturation are well studied in model organisms, the evolution of these systems is not entirely understood and requires further investigation. We here review the current knowledge of NR involvement in metamorphosis and sexual maturation across the animal tree of life with special attention to environmental integration and evolution of the signaling mechanism. Furthermore, we compare commonalities and differences of the different signaling systems. Finally, we identify key gaps in our knowledge of NR evolution, which, if sufficiently investigated, would lead to an importantly improved understanding of the evolution of complex signaling systems, the evolution of life history decision points, and, ultimately, speciation events in the metazoan kingdom.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Chemical Screening of Nuclear Receptor Modulators

Nuclear receptors are ligand-inducible transcriptional factors that control multiple biological phenomena, including proliferation, differentiation, reproduction, metabolism, and the maintenance of homeostasis. Members of the nuclear receptor superfamily have marked structural and functional similarities, and their domain functionalities and regulatory mechanisms have been well studied. Various modulators of nuclear receptors, including agonists and antagonists, have been developed as tools for elucidating nuclear receptor functions and also as drug candidates or lead compounds. Many assay systems are currently available to evaluate the modulation of nuclear receptor functions, and are useful as screening tools in the discovery and development of new modulators. In this review, we cover the chemical screening methods for nuclear receptor modulators, focusing on assay methods and chemical libraries for screening. We include some recent examples of the discovery of nuclear receptor modulators.

Comparative Analysis of the Organization of the Catecholaminergic Systems in the Brain of Holostean Fishes (Actinopterygii/Neopterygii)

Living holosteans, comprising 8 species of bowfins and gars, form a small monophyletic group of actinopterygian fishes, which are currently considered as the sister group to the enormously numerous teleosts and have largely been neglected in neuroanatomical studies. We have studied the catecholaminergic (CAergic) systems by means of antibodies against tyrosine hydroxylase (TH) and dopamine (DA) in the brain of representative species of the 3 genera included in the 2 orders of holostean fishes: Amia calva (Amiiformes) and Lepisosteus platyrhincus, Lepisosteus oculatus, and Atractosteus spatula (Lepisosteiformes). Different groups of TH/DA-immunoreactive (ir) cells were observed in the olfactory bulb, subpallium, and preoptic area of the telencephalon. Hypothalamic groups were labeled in the suprachiasmatic nucleus, tuberal (only in A. calva), retrotuberal, and retromamillary areas; specifically, the paraventricular organ showed only DA immunoreactivity. In the diencephalon, TH/DA-ir groups were detected in the prethalamus, posterior tubercle, and pretectum. In the caudal hindbrain, the solitary tract nucleus and area postrema presented TH/DA-ir cell groups, and also the spinal cord and the retina. Only in A. calva, particular CAergic cell groups were observed in the habenula, the mesencephalic tegmentum, and in the locus coeruleus. Following a neuromeric analysis, the comparison of these results with those obtained in other classes of fishes and tetrapods shows many common traits of CAergic systems shared by most vertebrates and in addition highlights unique features of actinopterygian fishes.

Duplication of toll-like receptor 22 in teleost fishes

The TLRs of teleost fishes have distinct features and are highly diverse, but the duplication characteristics and expression patterns of the tlr22 gene remain unclear. Here, we identified paralogous tlr22 genes in 13 teleost fishes by screening available fish genomic resources and using molecular cloning. We then conducted comprehensive bioinformatics analyses and investigated spatiotemporal differences in the expression patterns of the tlr22 genes in G. eckloni. The results indicated that more than three paralogous tlr22 genes were possessed by some teleost fishes. Of these, tlr22c is specific to some subfamilies of the Cyprinidae (e.g., Barbinae, Cyprininae, Schizothoracinae, and Leuciscinae). Phylogenetic and syntenic analyses showed that the paralogous tlr22 genes originated from two single-gene duplication events. Molecular clock calculations dated the two gene duplication events at 49.5 and 39.3 MYA, which is before the common carp-specific genome duplication event and well after the fish-specific genome duplication. Gene duplication of tlr22 was followed by gene loss or pseudogene events in certain lineages. Spatiotemporal expression differences between the three duplicated tlr22 genes from G. eckloni suggested that these genes diverged functionally after gene duplication.

The evolutionary road to invertebrate thyroid hormone signaling: Perspectives for endocrine disruption processes

Thyroid hormones (THs) are the only iodine-containing hormones that play fundamental roles in chordates and non-chordates. The chemical nature, mode of action and the synthesis of THs are well established in mammals and other vertebrates. Although thyroid-like hormones have been detected in protostomes and non-chordate deuterostomes, TH signaling is poorly understood as compared to vertebrates, particularly in protostomes. Therefore, the central objective of this article is to review TH system components and TH-induced effects in non-vertebrate chordates, non-chordate deuterostomes and protostomes based on available genomes and functional information. To accomplish this task, we integrate here the available knowledge on the THs signaling across non-vertebrate chordates, non-chordate deuterostomes and protostomes by considering studies encompassing TH system components and physiological actions of THs. We also address the possible interactions of thyroid disrupting chemicals and their effects in protostomes and non-chordate deuterostomes. Finally, the perspectives on current and future challenges are discussed.

Evolution of vertebrate spinal cord patterning

The vertebrate spinal cord is organized across three developmental axes, anterior-posterior (AP), dorsal-ventral (DV), and medial-lateral (ML). Patterning of these axes is regulated by canonical intercellular signaling pathways: the AP axis by Wnt, fibroblast growth factor, and retinoic acid (RA), the DV axis by Hedgehog, Tgfβ, and Wnt, and the ML axis where proliferation is controlled by Notch. Developmental time plays an important role in which signal does what and when. Patterning across the three axes is not independent, but linked by interactions between signaling pathway components and their transcriptional targets. Combined this builds a sophisticated organ with many different types of cell in specific AP, DV, and ML positions. Two living lineages share phylum Chordata with vertebrates, amphioxus, and tunicates, while the jawless fish such as lampreys, survive as the most basally divergent vertebrate lineage. Genes and mechanisms shared between lampreys and other vertebrates tell us what predated vertebrates, while those also shared with other chordates tell us what evolved early in chordate evolution. Between these lie vertebrate innovations: genetic and developmental changes linked to evolution of new morphology. These include gene duplications, differences in how signals are received, and new regulatory connections between signaling pathways and their target genes.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

Cartography of rhodopsin-like G protein-coupled receptors across vertebrate genomes

We conduct a cartography of rhodopsin-like non-olfactory G protein-coupled receptors in the Ensembl database. The most recent genomic data (releases 90–92, 90 vertebrate genomes) are analyzed through the online interface and receptors mapped on phylogenetic guide trees that were constructed based on a set of ~14.000 amino acid sequences. This snapshot of genomic data suggest vertebrate genomes to harbour 142 clades of GPCRs without human orthologues. Among those, 69 have not to our knowledge been mentioned or studied previously in the literature, of which 28 are distant from existing receptors and likely new orphans. These newly identified receptors are candidates for more focused evolutionary studies such as chromosomal mapping as well for in-depth pharmacological characterization. Interestingly, we also show that 37 of the 72 human orphan (or recently deorphanized) receptors included in this study cluster into nineteen closely related groups, which implies that there are less ligands to be identified than previously anticipated. Altogether, this work has significant implications when discussing nomenclature issues for GPCRs.

Transport, Metabolism, and Function of Thyroid Hormones in the Developing Mammalian Brain

Ever since the discovery of thyroid hormone deficiency as the primary cause of cretinism in the second half of the 19th century, the crucial role of thyroid hormone (TH) signaling in embryonic brain development has been established. However, the biological understanding of TH function in brain formation is far from complete, despite advances in treating thyroid function deficiency disorders. The pleiotropic nature of TH action makes it difficult to identify and study discrete roles of TH in various aspect of embryogenesis, including neurogenesis and brain maturation. These challenges notwithstanding, enormous progress has been achieved in understanding TH production and its regulation, their conversions and routes of entry into the developing mammalian brain. The endocrine environment has to adjust when an embryo ceases to rely solely on maternal source of hormones as its own thyroid gland develops and starts to produce endogenous TH. A number of mechanisms are in place to secure the proper delivery and action of TH with placenta, blood-brain interface, and choroid plexus as barriers of entry that need to selectively transport and modify these hormones thus controlling their active levels. Additionally, target cells also possess mechanisms to import, modify and bind TH to further fine-tune their action. A complex picture of a tightly regulated network of transport proteins, modifying enzymes, and receptors has emerged from the past studies. TH have been implicated in multiple processes related to brain formation in mammals-neuronal progenitor proliferation, neuronal migration, functional maturation, and survival-with their exact roles changing over developmental time. Given the plethora of effects thyroid hormones exert on various cell types at different developmental periods, the precise spatiotemporal regulation of their action is of crucial importance. In this review we summarize the current knowledge about TH delivery, conversions, and function in the developing mammalian brain. We also discuss their potential role in vertebrate brain evolution and offer future directions for research aimed at elucidating TH signaling in nervous system development.

Identification and Evolution of TGF-β Signaling Pathway Members in Twenty-Four Animal Species and Expression in Tilapia

Transforming growth factor β (TGF-β) signaling controls diverse cellular processes during embryogenesis as well as in mature tissues of multicellular animals. Here we carried out a comprehensive analysis of TGF-β pathway members in 24 representative animal species. The appearance of the TGF-β pathway was intrinsically linked to the emergence of metazoan. The total number of TGF-β ligands, receptors, and smads changed slightly in all invertebrates and jawless vertebrates analyzed. In contrast, expansion of the pathway members, especially ligands, was observed in jawed vertebrates most likely due to the second round of whole genome duplication (2R) and additional rounds in teleosts. Duplications of TGFB2, TGFBR2, ACVR1, SMAD4 and SMAD6, which were resulted from 2R, were first isolated. Type II receptors may be originated from the ACVR2-like ancestor. Interestingly, AMHR2 was not identified in Chimaeriformes and Cypriniformes even though they had the ligand AMH. Based on transcriptome data, TGF-β ligands exhibited a tissue-specific expression especially in the heart and gonads. However, most receptors and smads were expressed in multiple tissues indicating they were shared by different ligands. Spatial and temporal expression profiles of 8 genes in gonads of different developmental stages provided a fundamental clue for understanding their important roles in sex determination and reproduction. Taken together, our findings provided a global insight into the phylogeny and expression patterns of the TGF-β pathway genes, and hence contribute to the greater understanding of their biological roles in the organism especially in teleosts.

The Calcitonin/Calcitonin Gene-Related Peptide Family in Invertebrate Deuterostomes

Calcitonin (CT)/CT gene-related peptide (CGRP) family peptides (CT/CGRP family peptides) including CT, CGRP, adrenomedullin, amylin, and CT receptor-stimulating peptide have been identified from various vertebrates and perform a variety of important physiological functions. These peptides bind to two types of receptors including CT receptor (CTR) and CTR-like receptor (CLR). Receptor recognition of CT/CGRP family peptides is determined by the heterodimer between CTR/CLR and receptor activity-modifying protein (RAMP). Comparative studies of the CT/CGRP family have been exclusively performed in vertebrates from teleost fishes to mammals and strongly manifest that the CGRP family system containing peptides, their receptors, and RAMPs was derived from a common ancestor. In addition, CT/CGRP family peptides and their receptors are also identified and inferred from various invertebrate species. However, the evolutionary process of the CT/CGRP family from invertebrates to vertebrates remains enigmatic. In this review, I principally summarize the CT/CGRP family peptides and their receptors in invertebrate deuterostomes, highlighting the study of invertebrate chordates including ascidians and amphioxi. The CT/CGRP family peptide that shows similar molecular structure and function with that of vertebrate CT has been identified from ascidian, Ciona intestinalis. Amphioxus, Branchiostoma floridae also possessed three CT/CGRP family peptides, one CTR/CLR receptor, and three RAMP-like proteins. The molecular function of the receptor complex formed by amphioxus CTR/CLR and a RAMP-like protein was clarified. Moreover, CT/CGRP family peptides have been identified in the superphylum Ambulacraria, which is close to Chordata. Finally, this review provides potential hypotheses of the evolution of CGRP family peptides and their receptors from invertebrates to vertebrates.

Lamprey metamorphosis: Thyroid hormone signaling in a basal vertebrate

As one of the most basal living vertebrates, lampreys represent an excellent model system to study the evolution of thyroid hormone (TH) signaling. The lamprey hypothalamic-pituitary-thyroid and reproductive axes overlap functionally. Lampreys have 3 gonadotropin-releasing hormones and a single glycoprotein hormone from the hypothalamus and pituitary, respectively, that regulate both the reproductive and thyroid axes. TH synthesis in larval lampreys takes place in an endostyle that transforms into typical vertebrate thyroid tissue during metamorphosis; both the endostyle and follicular tissue have all the typical TH synthetic components found in other vertebrates. Furthermore, lampreys also have the vertebrate suite of peripheral regulators including TH distributor proteins (THDPs), deiodinases and TH receptors (TRs). Although at the molecular level the components of the lamprey thyroid system are ancestral to other vertebrates, their functions have been largely conserved. TH signaling as it relates to lamprey metamorphosis represents a particularly interesting phenomenon. Unlike other metamorphosing vertebrates, lamprey THs increase throughout the larval period, peak prior to metamorphosis and decline rapidly at the onset of metamorphosis; patterns of deiodinase activity are consistent with these increases and declines. Moreover, goitrogens (which suppress TH levels) initiate precocious metamorphosis, and exogenous TH treatment blocks goitrogen-induced metamorphosis and disrupts natural metamorphosis. Despite this clear physiological difference, TH action via TRs is consistent with higher vertebrates. Based on observations that TRs are upregulated in a tissue-specific fashion during morphogenesis and the finding that lamprey TRs upregulate genes via THs in a fashion similar to higher vertebrates, we propose the following hypothesis for further testing. THs have a dual role in lampreys where high TH levels promote larval feeding and growth and then at the onset of metamorphosis TH levels decrease rapidly; at this time the relatively low TH levels function via TRs in a fashion similar to that of other metamorphosing vertebrates.

BACH transcription factors in innate and adaptive immunity

BTB and CNC homology (BACH) proteins are transcriptional repressors of the basic region leucine zipper (bZIP) transcription factor family. Recent studies indicate widespread roles of BACH proteins in controlling the development and function of the innate and adaptive immune systems, including the differentiation of effector and memory cells of the B and T cell lineages, CD4+ regulatory T cells and macrophages. Here, we emphasize similarities at a molecular level in the cell-type-specific activities of BACH factors, proposing that competitive interactions of BACH proteins with transcriptional activators of the bZIP family form a common mechanistic theme underlying their diverse actions. The findings contribute to a general understanding of how transcriptional repressors shape lineage commitment and cell-type-specific functions through repression of alternative lineage programmes.

CRH peptide evolution occurred in three phases: Evidence from characterizing sea lamprey CRH system members

The corticotropin releasing hormone (CRH) system, which includes the CRH family of peptides, their receptors (CRHRs) and a binding protein (CRHBP), has been strongly conserved throughout vertebrate evolution. The identification of invertebrate homologues suggests this system evolved over 500 million years ago. However, the early vertebrate evolution of the CRH system is not understood. Current theory indicates that agnathans (hagfishes and lampreys) are monophyletic with a conservative evolution over the past 500million years and occupy a position at the root of vertebrate phylogeny. We isolated the cDNAs for three CRH family members, two CRHRs and a CRHBP from the sea lamprey, Petromyzon marinus. Two of the CRH peptides are related to the CRH/urotensin-1 (UI) lineage, whereas the other is a urocortin (Ucn) 3 orthologue. The predicted amino acid identity of CRH and UI is 61% but they possess distinct motifs indicative of each peptide, suggesting an early divergence of the two genes. Based on our findings we propose the CRH peptides evolved in at least 3 distinct phases. The first occurring prior to the agnathans gave rise to the CRH/UI-like and Ucn2/3-like paralogous lineages. The second was a partial sub-genomic duplication of the ancestral CRH/UI-like gene, but not the Ucn2/3-like gene, giving rise to the CRH and UI (Ucn) lineages. The third event which resulted in the appearance of Ucn2 and Ucn3 must have occurred after the evolution of the cartilaginous fishes. Interestingly, unlike other vertebrate CRHRs, we were unable to classify our two P. marinus receptors (designated CRHRα and CRHRβ) as either type 1 or type 2, indicating that this split evolved later in vertebrate evolution. A single CRHBP gene was found suggesting that either this gene has not been affected by the vertebrate genome duplications or there have been a series of paralogous gene deletions. This study suggests that P. marinus possess a functional CRH system that differs from that of the gnathostomes and may represent a model for the earliest functioning CRH system in vertebrates.

Evolutionary and Functional Diversification of the Vitamin D Receptor-Lithocholic Acid Partnership

The evolution, molecular behavior, and physiological function of nuclear receptors are of particular interest given their diverse roles in regulating essential biological processes. The vitamin D receptor (VDR) is well known for its canonical roles in calcium homeostasis and skeletal maintenance. Additionally, VDR has received an increased amount of attention due to the discovery of numerous non-calcemic functions, including the detoxification of lithocholic acid. Lithocholic acid is a toxic metabolite of chenodeoxycholic acid, a primary bile acid. The partnership between the VDR and lithocholic acid has been hypothesized to be a recent adaptation that evolved to mediate the detoxification and elimination of lithocholic acid from the gut. This partnership is speculated to be limited to higher vertebrates (birds and mammals), as lower vertebrates do not synthesize the parent compound of lithocholic acid. However, the molecular functions associated with the observed insensitivity of basal VDRs to lithocholic acid have not been explored. Here we characterize canonical nuclear receptor functions of VDRs from select species representing key nodes in vertebrate evolution and span a range of bile salt phenotypes. Competitive ligand binding assays revealed that the receptor's affinity for lithocholic acid is highly conserved across species, suggesting that lithocholic acid affinity is an ancient and non-adaptive trait. However, transient transactivation assays revealed that lithocholic acid-mediated VDR activation might have evolved more recently, as the non-mammalian receptors did not respond to lithocholic acid unless exogenous coactivator proteins were co-expressed. Subsequent functional assays indicated that differential lithocholic acid-mediated receptor activation is potentially driven by differential protein-protein interactions between VDR and nuclear receptor coregulator proteins. We hypothesize that the vitamin D receptor-lithocholic acid partnership evolved as a by-product of natural selection on the ligand-receptor partnership between the vitamin D receptor and the native VDR ligand: 1α,25-dihydroxyvitamin D3, the biologically active metabolite of vitamin D3.

The Case of the Missing Ancient Fungal Polyploids

Polyploidy-the increase in the number of whole chromosome sets-is an important evolutionary force in eukaryotes. Polyploidy is well recognized throughout the evolutionary history of plants and animals, where several ancient events have been hypothesized to be drivers of major evolutionary radiations. However, fungi provide a striking contrast: while numerous recent polyploids have been documented, ancient fungal polyploidy is virtually unknown. We present a survey of known fungal polyploids that confirms the absence of ancient fungal polyploidy events. Three hypotheses may explain this finding. First, ancient fungal polyploids are indeed rare, with unique aspects of fungal biology providing similar benefits without genome duplication. Second, fungal polyploids are not successful in the long term, leading to few extant species derived from ancient polyploidy events. Third, ancient fungal polyploids are difficult to detect, causing the real contribution of polyploidy to fungal evolution to be underappreciated. We consider each of these hypotheses in turn and propose that failure to detect ancient events is the most likely reason for the lack of observed ancient fungal polyploids. We examine whether existing data can provide evidence for previously unrecognized ancient fungal polyploidy events but discover that current resources are too limited. We contend that establishing whether unrecognized ancient fungal polyploidy events exist is important to ascertain whether polyploidy has played a key role in the evolution of the extensive complexity and diversity observed in fungi today and, thus, whether polyploidy is a driver of evolutionary diversifications across eukaryotes. Therefore, we conclude by suggesting ways to test the hypothesis that there are unrecognized polyploidy events in the deep evolutionary history of the fungi.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERβ) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

A histological atlas of the tissues and organs of neotenic and metamorphosed axolotl

Axolotl (Ambystoma Mexicanum) has been emerging as a promising model in stem cell and regeneration researches due to its exceptional regenerative capacity. Although it represents lifelong lasting neoteny, induction to metamorphosis with thyroid hormones (THs) treatment advances the utilization of Axolotl in various studies. It has been reported that amphibians undergo anatomical and histological remodeling during metamorphosis and this transformation is crucial for adaptation to terrestrial conditions. However, there is no comprehensive histological investigation regarding the morphological alterations of Axolotl organs and tissues throughout the metamorphosis. Here, we reveal the histological differences or resemblances between the neotenic and metamorphic axolotl tissues. In order to examine structural features and cellular organization of Axolotl organs, we performed Hematoxylin & Eosin, Luxol-Fast blue, Masson's trichrome, Alcian blue, Orcein and Weigart's staining. Stained samples from brain, gallbladder, heart, intestine, liver, lung, muscle, skin, spleen, stomach, tail, tongue and vessel were analyzed under the light microscope. Our findings contribute to the validation of the link between newly acquired functions and structural changes of tissues and organs as observed in tail, skin, gallbladder and spleen. We believe that this descriptive work provides new insights for a better histological understanding of both neotenic and metamorphic Axolotl tissues.

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.

Comparative analysis of zebrafish bone morphogenetic proteins 2, 4 and 16: molecular and evolutionary perspectives

BMP2, BMP4 and BMP16 form a subfamily of bone morphogenetic proteins acting as pleiotropic growth factors during development and as bone inducers during osteogenesis. BMP16 is the most recent member of this subfamily and basic data regarding protein structure and function, and spatio-temporal gene expression is still scarce. In this work, insights on BMP16 were provided through the comparative analysis of structural and functional data for zebrafish BMP2a, BMP2b, BMP4 and BMP16 genes and proteins, determined from three-dimensional models, patterns of gene expression during development and in adult tissues, regulation by retinoic acid and capacity to activate BMP-signaling pathway. Structures of Bmp2a, Bmp2b, Bmp4 and Bmp16 were found to be remarkably similar; with residues involved in receptor binding being highly conserved. All proteins could activate the BMP-signaling pathway, suggesting that they share a common function. On the contrary, stage- and tissue-specific expression of bmp2, bmp4 and bmp16 suggested the genes might be differentially regulated (e.g. different transcription factors, enhancers and/or regulatory modules) but also that they are involved in distinct physiological processes, although with the same function. Retinoic acid, a morphogen known to interact with BMP-signaling during bone formation, was shown to down-regulate the expression of bmp2, bmp4 and bmp16, although to different extents. Taxonomic and phylogenetic analyses indicated that bmp16 diverged before bmp2 and bmp4, is not restricted to teleost fish lineage as previously reported, and that it probably arose from a whole genomic duplication event that occurred early in vertebrate evolution and disappeared in various tetrapod lineages through independent events.

Identification of Thyroid Hormones and Functional Characterization of Thyroid Hormone Receptor in the Pacific Oyster Crassostrea gigas Provide Insight into Evolution of the Thyroid Hormone System

Thyroid hormones (THs) play important roles in development, metamorphosis, and metabolism in vertebrates. During the past century, TH functions were regarded as a synapomorphy of vertebrates. More recently, accumulating evidence has gradually convinced us that TH functions also occur in invertebrate chordates. To date, however, TH-related studies in non-chordate invertebrates have been limited. In this study, THs were qualitatively detected by two reliable methods (HPLC and LC/MS) in a well-studied molluscan species, the Pacific oyster Crassostrea gigas. Quantitative measurement of THs during the development of C. gigas showed high TH contents during embryogenesis and that oyster embryos may synthesize THs endogenously. As a first step in elucidating the TH signaling cascade, an ortholog of vertebrate TH receptor (TR), the most critical gene mediating TH effects, was cloned in C. gigas. The sequence of CgTR has conserved DNA-binding and ligand-binding domains that normally characterize these receptors. Experimental results demonstrated that CgTR can repress gene expression through binding to promoters of target genes and can interact with oyster retinoid X receptor. Moreover, CgTR mRNA expression was activated by T4 and the transcriptional activity of CgTR promoter was repressed by unliganded CgTR protein. An atypical thyroid hormone response element (CgDR5) was found in the promoter of CgTR, which was verified by electrophoretic mobility shift assay (EMSA). These results indicated that some of the CgTR function is conserved. However, the EMSA assay showed that DNA binding specificity of CgTR was different from that of the vertebrate TR and experiments with two dual-luciferase reporter systems indicated that l-thyroxine, 3,3′,5-triiodothyronine, and triiodothyroacetic acid failed to activate the transcriptional activity of CgTR. This is the first study to functionally characterize TR in mollusks. The presence of THs and the functions of CgTR in mollusks contribute to better understanding of the evolution of the TH system.

Bioinformatic Characterization of Genes and Proteins Involved in Blood Clotting in Lampreys

Lampreys and hagfish are the earliest diverging of extant vertebrates and are obvious targets for investigating the origins of complex biochemical systems found in mammals. Currently, the simplest approach for such inquiries is to search for the presence of relevant genes in whole genome sequence (WGS) assemblies. Unhappily, in the past a high-quality complete genome sequence has not been available for either lampreys or hagfish, precluding the possibility of proving gene absence. Recently, improved but still incomplete genome assemblies for two species of lamprey have been posted, and, taken together with an extensive collection of short sequences in the NCBI trace archive, they have made it possible to make reliable counts for specific gene families. Particularly, a multi-source tactic has been used to study the lamprey blood clotting system with regard to the presence and absence of genes known to occur in higher vertebrates. As was suggested in earlier studies, lampreys lack genes for coagulation factors VIII and IX, both of which are critical for the "intrinsic" clotting system and responsible for hemophilia in humans. On the other hand, they have three each of genes for factors VII and X, participants in the "extrinsic" clotting system. The strategy of using raw trace sequence "reads" together with partial WGS assemblies for lampreys can be used in studies on the early evolution of other biochemical systems in vertebrates.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.

Evolution of retinoic acid receptors in chordates: insights from three lamprey species, Lampetra fluviatilis, Petromyzon marinus, and Lethenteron japonicum

Background

Retinoic acid (RA) signaling controls many developmental processes in chordates, from early axis specification to late organogenesis. The functions of RA are chiefly mediated by a subfamily of nuclear hormone receptors, the retinoic acid receptors (RARs), that act as ligand-activated transcription factors. While RARs have been extensively studied in jawed vertebrates (that is, gnathostomes) and invertebrate chordates, very little is known about the repertoire and developmental roles of RARs in cyclostomes, which are extant jawless vertebrates. Here, we present the first extensive study of cyclostome RARs focusing on three different lamprey species: the European freshwater lamprey, Lampetra fluviatilis, the sea lamprey, Petromyzon marinus, and the Japanese lamprey, Lethenteron japonicum.

Results

We identified four rar paralogs (rar1, rar2, rar3, and rar4) in each of the three lamprey species, and phylogenetic analyses indicate a complex evolutionary history of lamprey rar genes including the origin of rar1 and rar4 by lineage-specific duplication after the lamprey-hagfish split. We further assessed their expression patterns during embryonic development by in situ hybridization. The results show that lamprey rar genes are generally characterized by dynamic and highly specific expression domains in different embryonic tissues. In particular, lamprey rar genes exhibit combinatorial expression domains in the anterior central nervous system (CNS) and the pharyngeal region.

Conclusions

Our results indicate that the genome of lampreys encodes at least four rar genes and suggest that the lamprey rar complement arose from vertebrate-specific whole genome duplications followed by a lamprey-specific duplication event. Moreover, we describe a combinatorial code of lamprey rar expression in both anterior CNS and pharynx resulting from dynamic and highly specific expression patterns during embryonic development. This ‘RAR code’ might function in regionalization and patterning of these two tissues by differentially modulating the expression of downstream effector genes during development.

Electronic supplementary material

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

Molecular cloning, functional characterization, and evolutionary analysis of vitamin D receptors isolated from basal vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC50 values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D3 acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D3 throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.

Current advances on ABC drug transporters in fish

Most members of the large ATP-binding cassette (ABC) gene family are transporters involved in substrate translocation across biological membranes. In eukaryotes, ABC proteins functioning as drug transporters are located in the plasma membrane and mediate the cellular efflux of a wide range of organic chemicals, with some transporters also transporting certain metals. As the enhanced expression of ABC drug transporters can confer multidrug resistance (MDR) to cancers and multixenobiotic resistance (MXR) to organisms from polluted habitats, these ABC family members are also referred to as MDR or MXR proteins. In mammals, ABC drug transporters show predominant expression in tissues involved in excretion or constituting internal or external body boundaries, where they facilitate the excretion of chemicals and their metabolites, and limit chemical uptake and penetration into "sanctuary" sites of the body. Available knowledge about ABC proteins is still limited in teleost fish, a large vertebrate group of high ecological and economic importance. Using transport activity measurements and immunochemical approaches, early studies demonstrated similarities in the tissue distribution of ABC drug transporters between teleosts and mammals, suggesting conserved roles of the transporters in the biochemical defence against toxicants. Recently, the availability of teleost genome assemblies has stimulated studies of the ABC family in this taxon. This review summarises the current knowledge regarding the genetics, functional properties, physiological function, and ecotoxicological relevance of teleostean ABC transporters. The available literature is reviewed with emphasis on recent studies addressing the tissue distribution, substrate spectrum, regulation, physiological function and phylogenetic origin of teleostean ABC transporters.

Thyroid hormone and retinoid X receptor function and expression during sea lamprey (Petromyzon marinus) metamorphosis

Sea lampreys (Petromyzon marinus) are members of the ancient class Agnatha and undergo a metamorphosis that transforms blind, sedentary, filter-feeding larvae into free-swimming, parasitic juveniles. Thyroid hormones (THs) appear to be important for lamprey metamorphosis, however, serum TH concentrations are elevated in the larval phase, decline rapidly during early metamorphosis and remain low until metamorphosis is complete; these TH fluctuations are contrary to those of other metamorphosing vertebrates. Moreover, thyroid hormone synthesis inhibitors (goitrogens) induce precocious metamorphosis and exogenous TH treatments disrupt natural metamorphosis in P. marinus. Given that THs exert their effects by binding to TH nuclear receptors (TRs) that often act as heterodimers with retinoid X receptors (RXRs), we cloned and characterized these receptors from P. marinus and examined their expression during metamorphosis. Two TRs (PmTR1 and PmTR2) and three RXRs (PmRXRs) were isolated from P. marinus cDNA. Phylogenetic analyses group the PmTRs together on a branch prior to the gnathostome TRα/β split. The three RXRs also group together, but our data indicated that these transcripts are most likely either allelic variants of the same gene locus, or the products of a lamprey-specific duplication event. Importantly, these P. marinus receptors more closely resemble vertebrate as opposed to invertebrate chordate receptors. Functional analysis revealed that PmTR1 and PmTR2 can activate transcription of TH-responsive genes when treated with nanomolar concentrations of TH and they have distinct pharmacological profiles reminiscent of vertebrate TRβ and TRα, respectively. Also similar to other metamorphosing vertebrates, expression patterns of the PmTRs during lamprey metamorphosis suggest that PmTR1 has a dynamic, tissue-specific expression pattern that correlates with tissue morphogenesis and biochemical changes and PmTR2 has a more uniform expression pattern. This TR expression data suggests that THs, either directly or via a metabolite, may function to positively modulate changes at the tissue or organ levels during lamprey metamorphosis. Collectively the results presented herein support the hypothesis that THs have a dual functional role in the lamprey life cycle whereby high levels promote larval feeding, growth and lipogenesis and low levels promote metamorphosis.

Nuclear receptors in nematode development: Natural experiments made by a phylum

The development of complex multicellular organisms is dependent on regulatory decisions that are necessary for the establishment of specific differentiation and metabolic cellular states. Nuclear receptors (NRs) form a large family of transcription factors that play critical roles in the regulation of development and metabolism of Metazoa. Based on their DNA binding and ligand binding domains, NRs are divided into eight NR subfamilies from which representatives of six subfamilies are present in both deuterostomes and protostomes indicating their early evolutionary origin. In some nematode species, especially in Caenorhabditis, the family of NRs expanded to a large number of genes strikingly exceeding the number of NR genes in vertebrates or insects. Nematode NRs, including the multiplied Caenorhabditis genes, show clear relation to vertebrate and insect homologues belonging to six of the eight main NR subfamilies. This review summarizes advances in research of nematode NRs and their developmental functions. Nematode NRs can reveal evolutionarily conserved mechanisms that regulate specific developmental and metabolic processes as well as new regulatory adaptations. They represent the results of a large number of natural experiments with structural and functional potential of NRs for the evolution of the phylum. The conserved and divergent character of nematode NRs adds a new dimension to our understanding of the general biology of regulation by NRs. This article is part of a Special Issue entitled: Nuclear receptors in animal development.

Mu opioids and their receptors: evolution of a concept

Opiates are among the oldest medications available to manage a number of medical problems. Although pain is the current focus, early use initially focused upon the treatment of dysentery. Opium contains high concentrations of both morphine and codeine, along with thebaine, which is used in the synthesis of a number of semisynthetic opioid analgesics. Thus, it is not surprising that new agents were initially based upon the morphine scaffold. The concept of multiple opioid receptors was first suggested almost 50 years ago (Martin, 1967), opening the possibility of new classes of drugs, but the morphine-like agents have remained the mainstay in the medical management of pain. Termed mu, our understanding of these morphine-like agents and their receptors has undergone an evolution in thinking over the past 35 years. Early pharmacological studies identified three major classes of receptors, helped by the discovery of endogenous opioid peptides and receptor subtypes-primarily through the synthesis of novel agents. These chemical biologic approaches were then eclipsed by the molecular biology revolution, which now reveals a complexity of the morphine-like agents and their receptors that had not been previously appreciated.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

Untangling the evolution of Rab G proteins: implications of a comprehensive genomic analysis

BACKGROUND

Membrane-bound organelles are a defining feature of eukaryotic cells, and play a central role in most of their fundamental processes. The Rab G proteins are the single largest family of proteins that participate in the traffic between organelles, with 66 Rabs encoded in the human genome. Rabs direct the organelle-specific recruitment of vesicle tethering factors, motor proteins, and regulators of membrane traffic. Each organelle or vesicle class is typically associated with one or more Rab, with the Rabs present in a particular cell reflecting that cell's complement of organelles and trafficking routes.

RESULTS

Through iterative use of hidden Markov models and tree building, we classified Rabs across the eukaryotic kingdom to provide the most comprehensive view of Rab evolution obtained to date. A strikingly large repertoire of at least 20 Rabs appears to have been present in the last eukaryotic common ancestor (LECA), consistent with the 'complexity early' view of eukaryotic evolution. We were able to place these Rabs into six supergroups, giving a deep view into eukaryotic prehistory.

CONCLUSIONS

Tracing the fate of the LECA Rabs revealed extensive losses with many extant eukaryotes having fewer Rabs, and none having the full complement. We found that other Rabs have expanded and diversified, including a large expansion at the dawn of metazoans, which could be followed to provide an account of the evolutionary history of all human Rabs. Some Rab changes could be correlated with differences in cellular organization, and the relative lack of variation in other families of membrane-traffic proteins suggests that it is the changes in Rabs that primarily underlies the variation in organelles between species and cell types.