The Nme family, previously known as Nm23 or NDPK, is involved in various molecular processes including tumor metastasis and some members of the family, but not all, exhibit a Nucleoside Diphosphate Kinase (NDPK) activity. Ten genes are known in humans, in which some members have been extensively studied. In non-mammalian species, the Nme protein family has received, in contrast, far less attention. The picture of the vertebrate Nme family remains thus incomplete and orthology relationships with mammalian counterparts were only partially characterized. The present study therefore aimed at characterizing the Nme gene repertoire in vertebrates with special interest for teleosts, and providing a comprehensive overview of the Nme gene family evolutionary history in vertebrates.

Pedomorphism is the retention of ancestrally juvenile traits by adults in a descendant taxon. Despite its importance for evolutionary change, there are few examples of a molecular basis for this phenomenon. Notothenioids represent one of the best described species flocks among marine fishes, but their diversity is currently threatened by the rapidly changing Antarctic climate. Notothenioid evolutionary history is characterized by parallel radiations from a benthic ancestor to pelagic predators, which was accompanied by the appearance of several pedomorphic traits, including the reduction of skeletal mineralization that resulted in increased buoyancy.

The molecular genetic mechanisms of sex determination are not known for most vertebrates, including zebrafish. We identified a mutation in the zebrafish fancl gene that causes homozygous mutants to develop as fertile males due to female-to-male sex reversal. Fancl is a member of the Fanconi Anemia/BRCA DNA repair pathway. Experiments showed that zebrafish fancl was expressed in developing germ cells in bipotential gonads at the critical time of sexual fate determination. Caspase-3 immunoassays revealed increased germ cell apoptosis in fancl mutants that compromised oocyte survival. In the absence of oocytes surviving through meiosis, somatic cells of mutant gonads did not maintain expression of the ovary gene cyp19a1a and did not down-regulate expression of the early testis gene amh; consequently, gonads masculinized and became testes. Remarkably, results showed that the introduction of a tp53 (p53) mutation into fancl mutants rescued the sex-reversal phenotype by reducing germ cell apoptosis and, thus, allowed fancl mutants to become fertile females. Our results show that Fancl function is not essential for spermatogonia and oogonia to become sperm or mature oocytes, but instead suggest that Fancl function is involved in the survival of developing oocytes through meiosis. This work reveals that Tp53-mediated germ cell apoptosis induces sex reversal after the mutation of a DNA-repair pathway gene by compromising the survival of oocytes and suggests the existence of an oocyte-derived signal that biases gonad fate towards the female developmental pathway and thereby controls zebrafish sex determination.

Genomic resources for hundreds of species of evolutionary, agricultural, economic, and medical importance are unavailable due to the expense of well-assembled genome sequences and difficulties with multigenerational studies. Teleost fish provide many models for human disease but possess anciently duplicated genomes that sometimes obfuscate connectivity. Genomic information representing a fish lineage that diverged before the teleost genome duplication (TGD) would provide an outgroup for exploring the mechanisms of evolution after whole-genome duplication. We exploited massively parallel DNA sequencing to develop meiotic maps with thrift and speed by genotyping F(1) offspring of a single female and a single male spotted gar (Lepisosteus oculatus) collected directly from nature utilizing only polymorphisms existing in these two wild individuals. Using Stacks, software that automates the calling of genotypes from polymorphisms assayed by Illumina sequencing, we constructed a map containing 8406 markers. RNA-seq on two map-cross larvae provided a reference transcriptome that identified nearly 1000 mapped protein-coding markers and allowed genome-wide analysis of conserved synteny. Results showed that the gar lineage diverged from teleosts before the TGD and its genome is organized more similarly to that of humans than teleosts. Thus, spotted gar provides a critical link between medical models in teleost fish, to which gar is biologically similar, and humans, to which gar is genomically similar. Application of our F(1) dense mapping strategy to species with no prior genome information promises to facilitate comparative genomics and provide a scaffold for ordering the numerous contigs arising from next generation genome sequencing.

Mammalian genomes can vary substantially in haploid chromosome number even within a small taxon (e.g., 3-40 among deer alone); in contrast, teleost fish genomes are stable (24-25 in 58% of teleosts), but we do not yet understand the mechanisms that account for differences in karyotype stability. Among perciform teleosts, platyfish (Xiphophorus maculatus) and medaka (Oryzias latipes) both have 24 chromosome pairs, but threespine stickleback (Gasterosteus aculeatus) and green pufferfish (Tetraodon nigroviridis) have just 21 pairs. To understand the evolution of teleost genomes, we made a platyfish meiotic map containing 16,114 mapped markers scored on 267 backcross fish. We tiled genomic contigs along the map to create chromosome-length genome assemblies. Genome-wide comparisons of conserved synteny showed that platyfish and medaka karyotypes remained remarkably similar with few interchromosomal translocations but with numerous intrachromosomal rearrangements (transpositions and inversions) since their lineages diverged ∼120 million years ago. Comparative genomics with platyfish shows how reduced chromosome numbers in stickleback and green pufferfish arose by fusion of pairs of ancestral chromosomes after their lineages diverged from platyfish ∼195 million years ago. Zebrafish and human genomes provide outgroups to root observed changes. These studies identify likely genome assembly errors, characterize chromosome fusion events, distinguish lineage-independent chromosome fusions, show that the teleost genome duplication does not appear to have accelerated the rate of translocations, and reveal the stability of syntenies and gene orders in teleost chromosomes over hundreds of millions of years.

Perchlorate is a ubiquitous environmental contaminant that has widespread endocrine disrupting effects in vertebrates, including threespine stickleback (Gasterosteus aculeatus). The target of perchlorate is thyroid tissue where it induces changes in the organization, activation, and morphology of thyroid follicles and surrounding tissues. To test the hypothesis that some phenotypes of perchlorate toxicity are not mediated by thyroid hormone, we chronically exposed stickleback beginning at fertilization to perchlorate (10, 30, 100ppm) or control water with and without supplementation of either iodide or thyroxine (T4). Stickleback were sampled across a one-year timespan to identify potential differences in responses to treatment combinations before and after sexual maturation. We found that most thyroid histomorphological phenotypes induced by perchlorate (follicle proliferation, reduced follicle area (adults only), colloid depletion, thyrocyte hypertrophy (subadults only)) were significantly ameliorated by exogenous iodide supplementation. In contrast, treatment with exogenous T4 did not correct any of the thyroid-specific histopathologies induced by perchlorate. Whole-body thyroid hormone concentrations were not significantly affected by perchlorate exposure; however, supplementation with iodide and T4 significantly increased T4 concentrations. This study also revealed an increased erythrocyte area in the thyroid region of perchlorate-exposed adults, while lipid droplet number increased in perchlorate-exposed subadults. Increased erythrocyte area was ameliorated by both iodide and T4, while neither supplement was able to correct lipid droplet number. Our finding on lipid droplets indicates that exposure to perchlorate in early development may have obesogenic effects.

STAT1 mediates response to interferons and regulates immunity, cell proliferation, apoptosis, and sensitivity of Fanconi Anemia cells to apoptosis after interferon signaling; the roles of STAT1 in embryos, however, are not understood. To explore embryonic functions of STAT1, we investigated stat1b, an unstudied zebrafish co-ortholog of human STAT1. Zebrafish stat1a encodes all five domains of the human STAT1-alpha splice form but, like the human STAT1-beta splice variant, stat1b lacks a complete transactivation domain; thus, two unlinked zebrafish paralogs encode protein forms translated from two splice variants of a single human gene, as expected by sub-functionalization after genome duplication. Phylogenetic and conserved synteny studies showed that stat1b and stat1a arose as duplicates in the teleost genome duplication (TGD) and clarified the evolutionary origin of STAT1, STAT2, STAT3, STAT4, STAT5A, STAT5B and STAT6 by tandem and genome duplication. RT-PCR revealed maternal expression of stat1a and stat1b. In situ hybridization detected stat1b but not stat1a expression in embryonic hematopoietic tissues. Morpholino knockdown of stat1b, but not stat1a, decreased expression of the myeloid and granulocyte markers spi and mpo and increased expression of the hematopoietic progenitor marker scl, the erythrocyte marker gata1, and hemoglobin. These results suggest that zebrafish Stat1b promotes myeloid development at the expense of erythroid development.

Differentiating cells interact with their extracellular environment over time. Chondrocytes embed themselves in a proteoglycan (PG)-rich matrix, then undergo a developmental transition, termed "maturation," when they express ihh to induce bone in the overlying tissue, the perichondrium. Here, we ask whether PGs regulate interactions between chondrocytes and perichondrium, using zebrafish mutants to reveal that cartilage PGs inhibit chondrocyte maturation, which ultimately dictates the timing of perichondral bone development. In a mutagenesis screen, we isolated a class of mutants with decreased cartilage matrix and increased perichondral bone. Positional cloning identified lesions in two genes, fam20b and xylosyltransferase1 (xylt1), both of which encode PG synthesis enzymes. Mutants failed to produce wild-type levels of chondroitin sulfate PGs, which are normally abundant in cartilage matrix, and initiated perichondral bone formation earlier than their wild-type siblings. Primary chondrocyte defects might induce the bone phenotype secondarily, because mutant chondrocytes precociously initiated maturation, showing increased and early expression of such markers as runx2b, collagen type 10a1, and ihh co-orthologs, and ihha mutation suppressed early perichondral bone in PG mutants. Ultrastructural analyses demonstrated aberrant matrix organization and also early cellular features of chondrocyte hypertrophy in mutants. Refining previous in vitro reports, which demonstrated that fam20b and xylt1 were involved in PG synthesis, our in vivo analyses reveal that these genes function in cartilage matrix production and ultimately regulate the timing of skeletal development.

Retinoic acid signaling is important for patterning the central nervous system, paired appendages, digestive tract, and other organs. To begin to investigate retinoic acid signaling in zebrafish, we determined orthologies between zebrafish and tetrapod retinoic acid receptors (Rars) and examined the expression patterns of rar genes during embryonic development. Analysis of phylogenies and conserved syntenies showed that the three cloned zebrafish rar genes include raraa and rarab, which are co-orthologs of tetrapod Rara, and rarg, which is the zebrafish ortholog of tetrapod Rarg. We did not, however, find an ortholog of Rarb. RNA in situ hybridization experiments showed that rarab and rarg, are maternally expressed. Zygotic expression of raraa occurs predominantly in the hindbrain, lateral mesoderm, and tailbud. Zygotic expression of rarab largely overlaps that of raraa, except that in later stages rarab is expressed more broadly in the brain and in the pectoral fin bud and pharyngeal arches. Zygotic expression of zebrafish rarg also overlaps the other two genes, but it is expressed more strongly in the posterior hindbrain beginning in late somitogenesis as well as in neural crest cells in the pharyngeal arches. Thus, these three genes have largely overlapping expression patterns and a few gene-specific expression domains. Knowledge of these expression patterns will guide the interpretation of the roles these genes play in development.

We present seven new vertebrate homologs of the prolactin-releasing hormone receptor (PRLHR) and show that these are found as two separate subtypes, PRLHR1 and PRLHR2. Analysis of a number of vertebrate sequences using phylogeny, pharmacology, and paralogon analysis indicates that the PRLHRs are likely to share a common ancestry with the neuropeptide Y (NPY) receptors. Moreover, a micromolar level of NPY was able to bind and inhibit completely the PRLH-evoked response in PRLHR1-expressing cells. We suggest that an ancestral PRLH peptide started coevolving with a redundant NPY binding receptor, which then became PRLHR, approximately 500 million years ago. The PRLHR1 subtype was shown to have a relatively high evolutionary rate compared to receptors with fixed peptide preference, which could indicate a drastic change in binding preference, thus supporting this hypothesis. This report suggests how gene duplication events can lead to novel peptide ligand/receptor interactions and hence spur the evolution of new physiological functions.

The Antarctic icefishes (family Channichthyidae, suborder Notothenioidei) constitute the only vertebrate taxon that fails to produce red blood cells. These fishes can be paired with closely related, but erythrocyte-producing, notothenioids to discover erythropoietic genes via representational difference analysis. Using a B30.2-domain-encoding DNA probe so derived from the hematopoietic kidney (pronephros) of a red-blooded Antarctic rockcod, Notothenia coriiceps, we discovered a related, novel gene, bloodthirsty (bty), that encoded a 547-residue protein that contains sequential RING finger, B Box, coiled-coil, and B30.2 domains. bty mRNA was expressed by the pronephric kidney of N. coriiceps at a steady-state level 10-fold greater than that found in the kidney of the icefish Chaenocephalus aceratus. To test the function of bty, we cloned the orthologous zebrafish gene from a kidney cDNA library. Whole-mount in situ hybridization of zebrafish embryos showed that bty mRNA was present throughout development and, after the mid-blastula transition, was expressed in the head and in or near the site of primitive erythropoiesis in the tail just prior to red cell production. One- to four-cell embryos injected with two distinct antisense morpholino oligonucleotides (MOs) targeted to the 5'-end of the bty mRNA failed to develop red cells, whereas embryos injected with 4- and 5-bp mismatch control MOs produced wild-type quantities of erythrocytes. The morphant phenotype was rescued by co-injection of synthetic bty mRNA containing an artificial 5'-untranslated region (UTR) with the antisense MO that bound the 5'-UTR of the wild-type bty transcript. Furthermore, the expression of genes that mark terminal erythroid differentiation was greatly reduced in the antisense-MO-treated embryos. We conclude that bty is likely to play a role in differentiation of the committed red cell progenitor.

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

In the frigid, oxygen-rich Southern Ocean (SO), Antarctic icefishes (Channichthyidae; Notothenioidei) evolved the ability to survive without producing erythrocytes and hemoglobin, the oxygen-transport system of virtually all vertebrates. Here, we integrate paleoclimate records with an extensive phylogenomic dataset of notothenioid fishes to understand the evolution of trait loss associated with climate change. In contrast to buoyancy adaptations in this clade, we find relaxed selection on the genetic regions controlling erythropoiesis evolved only after sustained cooling in the SO. This pattern is seen not only within icefishes but also occurred independently in other high-latitude notothenioids. We show that one species of the red-blooded dragonfish clade evolved a spherocytic anemia that phenocopies human patients with this disease via orthologous mutations. The genomic imprint of SO climate change is biased toward erythrocyte-associated conserved noncoding elements (CNEs) rather than to coding regions, which are largely preserved through pleiotropy. The drift in CNEs is specifically enriched near genes that are preferentially expressed late in erythropoiesis. Furthermore, we find that the hematopoietic marrow of icefish species retained proerythroblasts, which indicates that early erythroid development remains intact. Our results provide a framework for understanding the interactions between development and the genome in shaping the response of species to climate change.

The retina has a very high energy demand but lacks an internal blood supply in most vertebrates. Here we explore the hypothesis that oxygen diffusion limited the evolution of retinal morphology by reconstructing the evolution of retinal thickness and the various mechanisms for retinal oxygen supply, including capillarization and acid-induced haemoglobin oxygen unloading. We show that a common ancestor of bony fishes likely had a thin retina without additional retinal oxygen supply mechanisms and that three different types of retinal capillaries were gained and lost independently multiple times during the radiation of vertebrates, and that these were invariably associated with parallel changes in retinal thickness. Since retinal thickness confers multiple advantages to vision, we propose that insufficient retinal oxygen supply constrained the functional evolution of the eye in early vertebrates, and that recurrent origins of additional retinal oxygen supply mechanisms facilitated the phenotypic evolution of improved functional eye morphology.

Ancient bony fishes had heterocercal tails, like modern sharks and sturgeons, with asymmetric caudal fins and a vertebral column extending into an elongated upper lobe. Teleost fishes, in contrast, developed a homocercal tail characterized by two separate equal-sized fin lobes and the body axis not extending into the caudal fin. A similar heterocercal-to-homocercal transition occurs during teleost ontogeny, although the underlying genetic and developmental mechanisms for either transition remain unresolved. Here, we investigated the role of hox13 genes in caudal fin formation as these genes control posterior identity in animals. Analysis of expression profiles of zebrafish hox13 paralogs and phenotypes of CRISPR/Cas9-induced mutants showed that double hoxb13a and hoxc13a mutants fail to form a caudal fin. Furthermore, single mutants display heterocercal-like morphologies not seen since Mesozoic fossil teleosteomorphs. Relaxation of functional constraints after the teleost genome duplication may have allowed hox13 duplicates to neo- or subfunctionalize, ultimately contributing to the evolution of a homocercal tail in teleost fishes.

The mangrove killifish Kryptolebias marmoratus, and its close relative Kryptolebias hermaphroditus, are the only vertebrate species known to reproduce by self-fertilization due to functional ovotestis development. To improve our understanding of their genomes, we constructed a genetic map. First, a single F1 fish was made by artificial fertilization between K. marmoratus and K. hermaphroditus strains. F2 progeny were then obtained by self-fertilization of the F1 fish. We used RAD-seq to query genomic DNAs from the two parental strains, the F1 individual and 49 F2 progeny. Results identified 9904 polymorphic RAD-tags (DNA markers) that mapped to 24 linkage groups, corresponding to the haploid chromosome number of these species. The total length of the map was 1248 cM, indicating that about one recombination occurred for each of the 24 homologous chromosome pairs in each meiosis. Markers were not evenly distributed along the chromosomes: in all chromosomes, many markers (> 8% of the total markers for each chromosome) mapped to chromosome tips. Centromeres suppress recombination, and this uneven distribution is probably due to the species' acrocentric chromosomes. Mapped marker sequences were compared to genomic sequences of medaka and platyfish, the next most closely related species with sequenced genomes that are anchored to genetic maps. Results showed that each mangrove killifish chromosome corresponds to a single chromosome of both platyfish and medaka, suggesting strong conservation of chromosomes over 100 million years of evolution. Our genetic map provides a framework for the K. marmoratus/K. hermaphroditus genome sequence and an important resource for understanding the biology of hermaphroditism.

Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species.

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

Antarctic fish have adapted to the freezing waters of the Southern Ocean. Representative adaptations to this harsh environment include a constitutive heat shock response and the evolution of an antifreeze protein in the blood. Despite their adaptations to the cold, genome-wide studies have not yet been performed on these fish due to the lack of a sequenced genome. Notothenia coriiceps, the Antarctic bullhead notothen, is an endemic teleost fish with a circumpolar distribution and makes a good model to understand the genomic adaptations to constant sub-zero temperatures.

Several attributes intuitively considered to be typical mammalian features, such as complex behavior, live birth and malignant disease such as cancer, also appeared several times independently in lower vertebrates. The genetic mechanisms underlying the evolution of these elaborate traits are poorly understood. The platyfish, X. maculatus, offers a unique model to better understand the molecular biology of such traits. We report here the sequencing of the platyfish genome. Integrating genome assembly with extensive genetic maps identified an unexpected evolutionary stability of chromosomes in fish, in contrast to in mammals. Genes associated with viviparity show signatures of positive selection, identifying new putative functional domains and rare cases of parallel evolution. We also find that genes implicated in cognition show an unexpectedly high rate of duplicate gene retention after the teleost genome duplication event, suggesting a hypothesis for the evolution of the behavioral complexity in fish, which exceeds that found in amphibians and reptiles.