Evolutionary mutant models for human disease

Model Organisms Facilitate Rare Disease Diagnosis and Therapeutic Research

Efforts to identify the genetic underpinnings of rare undiagnosed diseases increasingly involve the use of next-generation sequencing and comparative genomic hybridization methods. These efforts are limited by a lack of knowledge regarding gene function, and an inability to predict the impact of genetic variation on the encoded protein function. Diagnostic challenges posed by undiagnosed diseases have solutions in model organism research, which provides a wealth of detailed biological information. Model organism geneticists are by necessity experts in particular genes, gene families, specific organs, and biological functions. Here, we review the current state of research into undiagnosed diseases, highlighting large efforts in North America and internationally, including the Undiagnosed Diseases Network (UDN) (Supplemental Material, File S1) and UDN International (UDNI), the Centers for Mendelian Genomics (CMG), and the Canadian Rare Diseases Models and Mechanisms Network (RDMM). We discuss how merging human genetics with model organism research guides experimental studies to solve these medical mysteries, gain new insights into disease pathogenesis, and uncover new therapeutic strategies.

Transcriptome assembly and candidate genes involved in nutritional programming in the swordtail fish Xiphophorus multilineatus

BACKGROUND

Nutritional programming takes place in early development. Variation in the quality and/or quantity of nutrients in early development can influence long-term health and viability. However, little is known about the mechanisms of nutritional programming. The live-bearing fish Xiphophorus multilineatus has the potential to be a new model for understanding these mechanisms, given prior evidence of nutritional programming influencing behavior and juvenile growth rate. We tested the hypotheses that nutritional programming would influence behaviors involved in energy homeostasis as well gene expression in X. multilineatus.

METHODS

We first examined the influence of both juvenile environment (varied in nutrition and density) and adult environment (varied in nutrition) on behaviors involved in energy acquisition and energy expenditure in adult male X. multilineatus. We also compared the behavioral responses across the genetically influenced size classes of males. Males stop growing at sexual maturity, and the size classes of can be identified based on phenotypes (adult size and pigment patterns). To study the molecular signatures of nutritional programming, we assembled a de novo transcriptome for X. multilineatus using RNA from brain, liver, skin, testis and gonad tissues, and used RNA-Seq to profile gene expression in the brains of males reared in low quality (reduced food, increased density) and high quality (increased food, decreased density) juvenile environments.

RESULTS

We found that both the juvenile and adult environments influenced the energy intake behavior, while only the adult environment influenced energy expenditure. In addition, there were significant interactions between the genetically influenced size classes and the environments that influenced energy intake and energy expenditure, with males from one of the four size classes (Y-II) responding in the opposite direction as compared to the other males examined. When we compared the brains of males of the Y-II size class reared in a low quality juvenile environment to males from the same size class reared in high quality juvenile environment, 131 genes were differentially expressed, including metabolism and appetite master regulator agrp gene.

DISCUSSION

Our study provides evidence for nutritional programming in X. multilineatus, with variation across size classes of males in how juvenile environment and adult diet influences behaviors involved in energy homeostasis. In addition, we provide the first transcriptome of X. multilineatus, and identify a group of candidate genes involved in nutritional programming.

Evolution of gremlin 2 in cetartiodactyl mammals: gene loss coincides with lack of upper jaw incisors in ruminants

Understanding the processes that give rise to genomic variability in extant species is an active area of research within evolutionary biology. With the availability of whole genome sequences, it is possible to quantify different forms of variability such as variation in gene copy number, which has been described as an important source of genetic variability and in consequence of phenotypic variability. Most of the research on this topic has been focused on understanding the biological significance of gene duplication, and less attention has been given to the evolutionary role of gene loss. Gremlin 2 is a member of the DAN gene family and plays a significant role in tooth development by blocking the ligand-signaling pathway of BMP2 and BMP4. The goal of this study was to investigate the evolutionary history of gremlin 2 in cetartiodactyl mammals, a group that possesses highly divergent teeth morphology. Results from our analyses indicate that gremlin 2 has experienced a mixture of gene loss, gene duplication, and rate acceleration. Although the last common ancestor of cetartiodactyls possessed a single gene copy, pigs and camels are the only cetartiodactyl groups that have retained gremlin 2. According to the phyletic distribution of this gene and synteny analyses, we propose that gremlin 2 was lost in the common ancestor of ruminants and cetaceans between 56.3 and 63.5 million years ago as a product of a chromosomal rearrangement. Our analyses also indicate that the rate of evolution of gremlin 2 has been accelerated in the two groups that have retained this gene. Additionally, the lack of this gene could explain the high diversity of teeth among cetartiodactyl mammals; specifically, the presence of this gene could act as a biological constraint. Thus, our results support the notions that gene loss is a way to increase phenotypic diversity and that gremlin 2 is a dispensable gene, at least in cetartiodactyl mammals.

Cavefish and the basis for eye loss

Animals have colonized the entire world from rather moderate to the harshest environments, some of these so extreme that only few animals are able to survive. Cave environments present such a challenge and obligate cave animals have adapted to perpetual darkness by evolving a multitude of traits. The most common and most studied cave characteristics are the regression of eyes and the overall reduction in pigmentation. Studying these traits can provide important insights into how evolutionary forces drive convergent and regressive adaptation. The blind Mexican cavefish (Astyanax mexicanus) has emerged as a useful model to study cave evolution owing to the availability of genetic and genomic resources, and the amenability of embryonic development as the different populations remain fertile with each other. In this review, we give an overview of our current knowledge underlying the process of regressive and convergent evolution using eye degeneration in cavefish as an example.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Xiphophorus and Medaka Cancer Models

Besides recently developed zebrafish cancer models, other fish species have been employed for many years as cancer models in laboratory studies. Two models, namely in Xiphophorus and medaka have proven useful in providing important clues to cancer etiology. Medaka is a complementary model to zebrafish in many areas of research since it offers similar resources and experimental tools. Xiphophorus provides the advantages of a natural ("evolutionary mutant") model with established genetics. Xiphophorus hybrids can develop spontaneous and radiation or carcinogen induced cancers. This chapter describes the tumor models in both species, which mainly focus on melanoma, and summarizes the main findings and future research directions.

Embryogenesis and early skeletogenesis in the antarctic bullhead notothen, Notothenia coriiceps

BACKGROUND

Environmental temperature influences rates of embryonic development, but a detailed staging series for vertebrate embryos developing in the subzero cold of Antarctic waters is not yet available from fertilization to hatching. Given projected warming of the Southern Ocean, it is imperative to establish a baseline to evaluate potential effects of changing climate on fish developmental dynamics.

RESULTS

We studied the Bullhead notothen (Notothenia coriiceps), a notothenioid fish inhabiting waters between -1.9 and +2 °C. In vitro fertilization produced embryos that progressed through cleavage, epiboly, gastrulation, segmentation, organogenesis, and hatching. We compared morphogenesis spatially and temporally to Zebrafish and medaka. Experimental animals hatched after about 6 months to early larval stages. To help understand skeletogenesis, we analyzed late embryos for expression of sox9 and runx2, which regulate chondrogenesis, osteogenesis, and eye development. Results revealed that, despite their prolonged developmental time course, N. coriiceps embryos developed similarly to those of other teleosts with large yolk cells.

CONCLUSIONS

Our studies set the stage for future molecular analyses of development in these extremophile fish. Results provide a foundation for understanding the impact of ocean warming on embryonic development and larval recruitment of notothenioid fish, which are key factors in the marine trophic system. Developmental Dynamics 245:1066-1080, 2016. © 2016 Wiley Periodicals, Inc.

Cichlid fishes as a model to understand normal and clinical craniofacial variation

We have made great strides towards understanding the etiology of craniofacial disorders, especially for 'simple' Mendelian traits. However, the facial skeleton is a complex trait, and the full spectrum of genetic, developmental, and environmental factors that contribute to its final geometry remain unresolved. Forward genetic screens are constrained with respect to complex traits due to the types of genes and alleles commonly identified, developmental pleiotropy, and limited information about the impact of environmental interactions. Here, we discuss how studies in an evolutionary model - African cichlid fishes - can complement traditional approaches to understand the genetic and developmental origins of complex shape. Cichlids exhibit an unparalleled range of natural craniofacial morphologies that model normal human variation, and in certain instances mimic human facial dysmorphologies. Moreover, the evolutionary history and genomic architecture of cichlids make them an ideal system to identify the genetic basis of these phenotypes via quantitative trait loci (QTL) mapping and population genomics. Given the molecular conservation of developmental genes and pathways, insights from cichlids are applicable to human facial variation and disease. We review recent work in this system, which has identified lbh as a novel regulator of neural crest cell migration, determined the Wnt and Hedgehog pathways mediate species-specific bone morphologies, and examined how plastic responses to diet modulate adult facial shapes. These studies have not only revealed new roles for existing pathways in craniofacial development, but have identified new genes and mechanisms involved in shaping the craniofacial skeleton. In all, we suggest that combining work in traditional laboratory and evolutionary models offers significant potential to provide a more complete and comprehensive picture of the myriad factors that are involved in the development of complex traits.

Perchlorate Exposure Reduces Primordial Germ Cell Number in Female Threespine Stickleback

Perchlorate is a common aquatic contaminant that has long been known to affect thyroid function in vertebrates, including humans. More recently perchlorate has been shown to affect primordial sexual differentiation in the aquatic model fishes zebrafish and threespine stickleback, but the mechanism has been unclear. Stickleback exposed to perchlorate from fertilization have increased androgen levels in the embryo and disrupted reproductive morphologies as adults, suggesting that perchlorate could disrupt the earliest stages of primordial sexual differentiation when primordial germ cells (PGCs) begin to form the gonad. Female stickleback have three to four times the number of PGCs as males during the first weeks of development. We hypothesized that perchlorate exposure affects primordial sexual differentiation by reducing the number of germ cells in the gonad during an important window of stickleback sex determination at 14-18 days post fertilization (dpf). We tested this hypothesis by quantifying the number of PGCs at 16 dpf in control and 100 mg/L perchlorate-treated male and female stickleback. Perchlorate exposure from the time of fertilization resulted in significantly reduced PGC number only in genotypic females, suggesting that the masculinizing effects of perchlorate observed in adult stickleback may result from early changes to the number of PGCs at a time critical for sex determination. To our knowledge, this is the first evidence of a connection between an endocrine disruptor and reduction in PGC number prior to the first meiosis during sex determination. These findings suggest that a mode of action of perchlorate on adult reproductive phenotypes in vertebrates, including humans, such as altered fecundity and sex reversal or intersex gonads, may stem from early changes to germ cell development.

"Wrecks of Ancient Life": Genetic Variants Vetted by Natural Selection

The Genetics Society of America’s George W. Beadle Award honors individuals who have made outstanding contributions to the community of genetics researchers and who exemplify the qualities of its namesake as a respected academic, administrator, and public servant. The 2015 recipient is John Postlethwait. He has made groundbreaking contributions in developing the zebrafish as a molecular genetic model and in understanding the evolution of new gene functions in vertebrates. He built the first zebrafish genetic map and showed that its genome, along with that of distantly related teleost fish, had been duplicated. Postlethwait played an integral role in the zebrafish genome-sequencing project and elucidated the genomic organization of several fish species. Postlethwait is also honored for his active involvement with the zebrafish community, advocacy for zebrafish as a model system, and commitment to driving the field forward.

Fish-T1K (Transcriptomes of 1,000 Fishes) Project: large-scale transcriptome data for fish evolution studies

Ray-finned fishes (Actinopterygii) represent more than 50 % of extant vertebrates and are of great evolutionary, ecologic and economic significance, but they are relatively underrepresented in ‘omics studies. Increased availability of transcriptome data for these species will allow researchers to better understand changes in gene expression, and to carry out functional analyses. An international project known as the “Transcriptomes of 1,000 Fishes” (Fish-T1K) project has been established to generate RNA-seq transcriptome sequences for 1,000 diverse species of ray-finned fishes. The first phase of this project has produced transcriptomes from more than 180 ray-finned fishes, representing 142 species and covering 51 orders and 109 families. Here we provide an overview of the goals of this project and the work done so far.

Evolution by gene loss

The recent increase in genomic data is revealing an unexpected perspective of gene loss as a pervasive source of genetic variation that can cause adaptive phenotypic diversity. This novel perspective of gene loss is raising new fundamental questions. How relevant has gene loss been in the divergence of phyla? How do genes change from being essential to dispensable and finally to being lost? Is gene loss mostly neutral, or can it be an effective way of adaptation? These questions are addressed, and insights are discussed from genomic studies of gene loss in populations and their relevance in evolutionary biology and biomedicine.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

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.

Developmental basis of phenotypic integration in two Lake Malawi cichlids

Background

Cichlid fishes from the Rift Lakes of East Africa have undergone the most spectacular adaptive radiations in vertebrate history. Eco-morphological adaptations in lakes Victoria, Malawi and Tanganyika have resulted in a vast array of skull shapes and sizes, yet primary axes of morphological variation are conserved in all three radiations, prominently including the size of the preorbital region of the skull. This conserved pattern suggests that development may constrain the trajectories of cichlid head morphological evolution.

Results

Here, we (1) present a comparative analysis of adult head morphology in two sand-dweller cichlids from Lake Malawi with preorbital size differences representative of the main axis of variation among the three lakes and (2) analyze the ontogeny of shape and size differences by focusing on known developmental modules throughout the head. We find that (1) developmental differences between the two species correlate with known developmental modules; (2) differences in embryonic cartilage development result in phenotypically integrated changes among all bones derived from a single cartilage, while differences in dermal bone development tend to influence isolated regions within a bone; and lastly (3) species-specific morphologies appear in the embryo as subtle differences, which become progressively amplified throughout ontogeny. We propose that this amplification takes place at skeletal growth zones, the locations and shapes of which are patterned during embryogenesis.

Conclusions

This study is the most anatomically comprehensive analysis of the developmental differences underlying cichlid skull evolution in the Rift Lakes of East Africa. The scale of our analysis reveals previously unnoticed correlations between developmental modules and patterns of phenotypic integration. We propose that the primary axes of morphological variation among East African cichlid adaptive radiations are constrained by the hierarchical modularity of the teleost head skeleton.

Natural bone fragmentation in the blind cave-dwelling fish, Astyanax mexicanus: candidate gene identification through integrative comparative genomics

Animals that colonize dark and nutrient-poor subterranean environments evolve numerous extreme phenotypes. These include dramatic changes to the craniofacial complex, many of which are under genetic control. These phenotypes can demonstrate asymmetric genetic signals wherein a QTL is detected on one side of the face but not the other. The causative gene(s) underlying QTL are difficult to identify with limited genomic resources. We approached this task by searching for candidate genes mediating fragmentation of the third suborbital bone (SO3) directly inferior to the orbit of the eye. We integrated positional genomic information using emerging Astyanax resources, and linked these intervals to homologous (syntenic) regions of the Danio rerio genome. We identified a discrete, approximately 6 Mb, conserved region wherein the gene causing SO3 fragmentation likely resides. We interrogated this interval for genes demonstrating significant differential expression using mRNA-seq analysis of cave and surface morphs across life history. We then assessed genes with known roles in craniofacial evolution and development based on GO term annotation. Finally, we screened coding sequence alterations in this region, identifying two key genes: transforming growth factor β3 (tgfb3) and bone morphogenetic protein 4 (bmp4). Of these candidates, tgfb3 is most promising as it demonstrates significant differential expression across multiple stages of development, maps close (<1 Mb) to the fragmentation critical locus, and is implicated in a variety of other animal systems (including humans) in non-syndromic clefting and malformations of the cranial sutures. Both abnormalities are analogous to the failure-to-fuse phenotype that we observe in SO3 fragmentation. This integrative approach will enable discovery of the causative genetic lesions leading to complex craniofacial features analogous to human craniofacial disorders. This work underscores the value of cave-dwelling fish as a powerful evolutionary model of craniofacial disease, and demonstrates the power of integrative system-level studies for informing the genetic basis of craniofacial aberrations in nature.

Innate immune responses to gut microbiota differ between oceanic and freshwater threespine stickleback populations

Animal hosts must co-exist with beneficial microbes while simultaneously being able to mount rapid, non-specific, innate immune responses to pathogenic microbes. How this balance is achieved is not fully understood, and disruption of this relationship can lead to disease. Excessive inflammatory responses to resident microbes are characteristic of certain gastrointestinal pathologies such as inflammatory bowel disease (IBD). The immune dysregulation of IBD has complex genetic underpinnings that cannot be fully recapitulated with single-gene-knockout models. A deeper understanding of the genetic regulation of innate immune responses to resident microbes requires the ability to measure immune responses in the presence and absence of the microbiota using vertebrate models with complex genetic variation. Here, we describe a new gnotobiotic vertebrate model to explore the natural genetic variation that contributes to differences in innate immune responses to microbiota. Threespine stickleback, Gasterosteus aculeatus, has been used to study the developmental genetics of complex traits during the repeated evolution from ancestral oceanic to derived freshwater forms. We established methods to rear germ-free stickleback larvae and gnotobiotic animals monoassociated with single bacterial isolates. We characterized the innate immune response of these fish to resident gut microbes by quantifying the neutrophil cells in conventionally reared monoassociated or germ-free stickleback from both oceanic and freshwater populations grown in a common intermediate salinity environment. We found that oceanic and freshwater fish in the wild and in the laboratory share many intestinal microbial community members. However, oceanic fish mount a strong immune response to residential microbiota, whereas freshwater fish frequently do not. A strong innate immune response was uniformly observed across oceanic families, but this response varied among families of freshwater fish. The gnotobiotic stickleback model that we have developed therefore provides a platform for future studies mapping the natural genetic basis of the variation in immune response to microbes.

Summary: We developed a gnotobiotic threespine stickleback system to identify the cellular and genetic bases for variation in inflammatory responses to microbiota.

Molecular basis and drug sensitivity of the delayed rectifier (IKr) in the fish heart

Fishes are increasingly used as models for human cardiac diseases, creating a need for a better understanding of the molecular basis of fish cardiac ion currents. To this end we cloned KCNH6 channel of the crucian carp (Carassius carassius) that produces the rapid component of the delayed rectifier K(+) current (IKr), the main repolarising current of the fish heart. KCNH6 (ccErg2) was the main isoform of the Kv11 potassium channel family with relative transcript levels of 98.9% and 99.6% in crucian carp atrium and ventricle, respectively. KCNH2 (ccErg1), an orthologue to human cardiac Erg (Herg) channel, was only slightly expressed in the crucian carp heart. The native atrial IKr and the cloned ccErg2 were inhibited by similar concentrations of verapamil, terfenadine and KB-R7943 (P>0.05), while the atrial IKr was about an order of magnitude more sensitive to E-4031 than ccErg2 (P<0.05) suggesting that some accessory β-subunits may be involved. Sensitivity of the crucian carp atrial IKr to E-4031, terfenadine and KB-R7943 was similar to what has been reported for the Herg channel. In contrast, the sensitivity of the crucian carp IKr to verapamil was approximately 30 times higher than the previously reported values for the Herg current. In conclusion, the cardiac IKr is produced by non-orthologous gene products in fish (Erg2) and mammalian hearts (Erg1) and some marked differences exist in drug sensitivity between fish and mammalian Erg1/2 which need to be taken into account when using fish heart as a model for human heart.

Development of the Synarcual in the Elephant Sharks (Holocephali; Chondrichthyes): Implications for Vertebral Formation and Fusion

The synarcual is a structure incorporating multiple elements of two or more anterior vertebrae of the axial skeleton, forming immediately posterior to the cranium. It has been convergently acquired in the fossil group ‘Placodermi’, in Chondrichthyes (Holocephali, Batoidea), within the teleost group Syngnathiformes, and to varying degrees in a range of mammalian taxa. In addition, cervical vertebral fusion presents as an abnormal pathology in a variety of human disorders. Vertebrae develop from axially arranged somites, so that fusion could result from a failure of somite segmentation early in development, or from later heterotopic development of intervertebral bone or cartilage. Examination of early developmental stages indicates that in the Batoidea and the ‘Placodermi’, individual vertebrae developed normally and only later become incorporated into the synarcual, implying regular somite segmentation and vertebral development. Here we show that in the holocephalan Callorhinchus milii, uniform and regular vertebral segmentation also occurs, with anterior individual vertebra developing separately with subsequent fusion into a synarcual. Vertebral elements forming directly behind the synarcual continue to be incorporated into the synarcual through growth. This appears to be a common pattern through the Vertebrata. Research into human disorders, presenting as cervical fusion at birth, focuses on gene misexpression studies in humans and other mammals such as the mouse. However, in chondrichthyans, vertebral fusion represents the normal morphology, moreover, taxa such Leucoraja (Batoidea) and Callorhinchus (Holocephali) are increasingly used as laboratory animals, and the Callorhinchus genome has been sequenced and is available for study. Our observations on synarcual development in three major groups of early jawed vertebrates indicate that fusion involves heterotopic cartilage and perichondral bone/mineralised cartilage developing outside the regular skeleton. We suggest that chondrichthyans have potential as ideal extant models for identifying the genes involved in these processes, for application to human skeletal heterotopic disorders.

The rise of Astyanax cavefish

Numerous animals have invaded subterranean caverns and evolved remarkably similar features. These features include loss of vision and pigmentation, and gains in nonvisual sensation. This broad convergence echoes smaller-scale convergence, in which members of the same species repeatedly evolve the same cave-associated phenotypes. The blind Mexican tetra of the Sierra de El Abra region of northeastern Mexico has a complex origin, having recurrently colonized subterranean environments through numerous invasions of surface-dwelling fish. These colonizations likely occurred ∼1-5 MYa. Despite evidence of historical and contemporary gene flow between cave and surface forms, the cave-associated phenotype appears to remain quite stable in nature. This model system has provided insight to the mechanisms of phenotypic regression, the genetic basis for constructive trait evolution, and the origin of behavioral novelties. Here, we document the rise of this model system from its discovery by a Mexican surveyor in 1936, to a powerful system for cave biology and contemporary genetic research. The recently sequenced genome provides exciting opportunities for future research, and will help resolve several long-standing biological problems. Developmental Dynamics 244:1031-1038, 2015. © 2015 Wiley Periodicals, Inc.

Gene expression responses of threespine stickleback to salinity: implications for salt-sensitive hypertension

Despite recent success with genome-wide association studies (GWAS), identifying hypertension (HTN)-susceptibility loci in the general population remains difficult. Here, we present a novel strategy to address this challenge by studying salinity adaptation in the threespine stickleback, a fish species with diverse salt-handling ecotypes. We acclimated native freshwater (FW) and anadromous saltwater (SW) threespine sticklebacks to fresh, brackish, and sea water for 30 days, and applied RNA sequencing to determine the gene expression in fish kidneys. We identified 1844 salt-responsive genes that were differentially expressed between FW sticklebacks acclimated to different salinities and/or between SW and FW sticklebacks acclimated to full-strength sea water. Significant overlap between stickleback salt-responsive genes and human genes implicated in HTN was detected (P < 10⁻⁷, hypergeometric test), suggesting a possible similarity in genetic mechanisms of salt handling between threespine sticklebacks and humans. The overlapping genes included a newly discovered HTN gene—MAP3K15, whose expression in FW stickleback kidneys decreases with salinity. These also included genes located in the GWAS loci such as AGTRAP-PLOD1 and CYP1A1-ULK3, which contain multiple potentially causative genes contributing to HTN susceptibility that need to be prioritized for study. Taken together, we show that stickleback salt-responsive genes provide valuable information facilitating the identification of human HTN genes. Thus, threespine sticklebacks may be used as a model, complementary to existing animal models, in human HTN research.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

Divergence in skeletal mass and bone morphology in antarctic notothenioid fishes

Although notothenioid fishes lack swim bladders, some species live temporarily or permanently in the water column. Given its relatively high density, skeletal mass is a key determinant of buoyancy. Notothenioids have reduced skeletal ossification, but there is little quantitative data on the phylogenetic distribution of this trait. We obtained dry skeletal masses for 54 specimens representing 20 species from six notothenioid families. Although comparative data are sparse, notothenioid skeletons comprise a smaller percentage of body mass, <3.5%, than those of three non-notothenioid perciforms. With relatively high skeletal mass, the non-Antarctic Bovichtus diacanthus is similar in skeletal mass to some non-notothenioids. Eleginops maclovinus, the non-Antarctic sister group of the Antarctic clade, has a relatively light skeleton (<2% of body mass) similar to many species in the Antarctic clade. Low skeletal mass is therefore a synapomorphy shared by Eleginops plus the Antarctic clade. We provide gross, histological, and micro-CT documentation of the structure and location of bone and cartilage in skulls, pectoral girdles, and vertebrae, with emphasis on the bovichtid B. diacanthus, the eleginopsid E. maclovinus, and the channichthyid Chaenodraco wilsoni. In Eleginops and the Antarctic clade, most bone is spongy and most species have persisting cartilage in the skull and appendicular skeleton. We also measured the relative size of the notochordal canal in adult vertebral centra of 38 species representing all eight families. There is considerable interspecific variation in this pedomorphic trait and all species show an ontogenetic reduction in the relative size of the canal. However, large persisting canals are present in adults of the Antarctic clade, especially in the nototheniids Pleuragramma and Aethotaxis and in a number of bathydraconid and channichthyid genera.

Beyond the zebrafish: diverse fish species for modeling human disease

In recent years, zebrafish, and to a lesser extent medaka, have become widely used small animal models for human diseases. These organisms have convincingly demonstrated the usefulness of fish for improving our understanding of the molecular and cellular mechanisms leading to pathological conditions, and for the development of new diagnostic and therapeutic tools. Despite the usefulness of zebrafish and medaka in the investigation of a wide spectrum of traits, there is evidence to suggest that other fish species could be better suited for more targeted questions. With the emergence of new, improved sequencing technologies that enable genomic resources to be generated with increasing efficiency and speed, the potential of non-mainstream fish species as disease models can now be explored. A key feature of these fish species is that the pathological condition that they model is often related to specific evolutionary adaptations. By exploring these adaptations, new disease-causing and disease-modifier genes might be identified; thus, diverse fish species could be exploited to better understand the complexity of disease processes. In addition, non-mainstream fish models could allow us to study the impact of environmental factors, as well as genetic variation, on complex disease phenotypes. This Review will discuss the opportunities that such fish models offer for current and future biomedical research.

Deficiency of zebrafish fgf20a results in aberrant skull remodeling that mimics both human cranial disease and evolutionarily important fish skull morphologies

The processes that direct skull remodeling are of interest to both human-oriented studies of cranial dysplasia and evolutionary studies of skull divergence. There is increasing awareness that these two fields can be mutually informative when natural variation mimics pathology. Here we describe a zebrafish mutant line, devoid of blastema (dob), which does not have a functional fgf20a protein, and which also presents cranial defects similar to both adaptive and clinical variation. We used geometric morphometric methods to provide quantitative descriptions of the effects of the dob mutation on skull morphogenesis. In combination with "whole-mount in situ hybridization" labeling of normal fgf20a expression and assays for osteoblast and osteoclast activity, the results of these analyses indicate that cranial dysmorphologies in dob zebrafish are generated by aberrations in post-embryonic skull remodeling via decreased osteoblasotgenesis and increased osteoclastogenesis. Mutational effects include altered skull vault geometries and midfacial hypoplasia that are consistent with key diagnostic signs for multiple human craniofacial syndromes. These phenotypic shifts also mimic changes in the functional morphology of fish skulls that have arisen repeatedly in several highly successful radiations (e.g., damselfishes and East-African rift-lake cichlids). Our results offer the dob/fgf20a mutant as an experimentally tractable model with which to examine post-embryonic skull development as it relates to human disease and vertebrate evolution.

Evolution of the eye transcriptome under constant darkness in Sinocyclocheilus cavefish

In adaptating to perpetual darkness, cave species gradually lose eyes and body pigmentation and evolve alternatives for exploring their environments. Although troglodyte features evolved independently many times in cavefish, we do not yet know whether independent evolution of these characters involves common genetic mechanisms. Surface-dwelling and many cave-dwelling species make the freshwater teleost genus Sinocyclocheilus an excellent model for studying the evolution of adaptations to life in constant darkness. We compared the mature retinal histology of surface and cave species in Sinocyclocheilus and found that adult cavefish showed a reduction in the number and length of photoreceptor cells. To identify genes and genetic pathways that evolved in constant darkness, we used RNA-seq to compare eyes of surface and cave species. De novo transcriptome assemblies were developed for both species, and contigs were annotated with gene ontology. Results from cave-dwelling Sinocyclocheilus revealed reduced transcription of phototransduction and other genes important for retinal function. In contrast to the blind Mexican tetra cavefish Astyanax mexicanus, our results on morphologies and gene expression suggest that evolved retinal reduction in cave-dwelling Sinocyclocheilus occurs in a lens-independent fashion by the reduced proliferation and downregulation of transcriptional factors shown to have direct roles in retinal development and maintenance, including cone-rod homeobox (crx) and Wnt pathway members. These results show that the independent evolution of retinal degeneration in cavefish can occur by different developmental genetic mechanisms.

Quantitative genetic analysis of retinal degeneration in the blind cavefish Astyanax mexicanus

The retina is the light-sensitive tissue of the eye that facilitates vision. Mutations within genes affecting eye development and retinal function cause a host of degenerative visual diseases, including retinitis pigmentosa and anophthalmia/microphthalmia. The characin fish Astyanax mexicanus includes both eyed (surface fish) and eyeless (cavefish) morphs that initially develop eyes with normal retina; however, early in development, the eyes of cavefish degenerate. Since both surface and cave morphs are members of the same species, they serve as excellent evolutionary mutant models with which to identify genes causing retinal degeneration. In this study, we crossed the eyed and eyeless forms of A. mexicanus and quantified the thickness of individual retinal layers among 115 F₂ hybrid progeny. We used next generation sequencing (RAD-seq) and microsatellite mapping to construct a dense genetic map of the Astyanax genome, scan for quantitative trait loci (QTL) affecting retinal thickness, and identify candidate genes within these QTL regions. The map we constructed for Astyanax includes nearly 700 markers assembled into 25 linkage groups. Based on our scans with this map, we identified four QTL, one each associated with the thickness of the ganglion, inner nuclear, outer plexiform, and outer nuclear layers of the retina. For all but one QTL, cavefish alleles resulted in a clear reduction in the thickness of the affected layer. Comparative mapping of genetic markers within each QTL revealed that each QTL corresponds to an approximately 35 Mb region of the zebrafish genome. Within each region, we identified several candidate genes associated with the function of each affected retinal layer. Our study is the first to examine Astyanax retinal degeneration in the context of QTL mapping. The regions we identify serve as a starting point for future studies on the genetics of retinal degeneration and eye disease using the evolutionary mutant model Astyanax.

Epigenetic synthesis: a need for a new paradigm for evolution in a contaminated world

Epigenetics is a perspective, not a set of techniques. Molecular biology and genetics are the dominant disciplines in biology today, but practitioners of these fields have only recently ‘rediscovered’ the importance of the environment. This has led to increasing research into molecular epigenetics and the interface between the environment and gene regulation. Beyond the study of epigenetic mechanisms at the level of the gene, more investigation of epigenetic outcomes at the level of both the individual organism and the evolution of the population is needed.

Impact of silicon-based quantum dots on the antioxidative system in white muscle of Carassius auratus gibelio

Silicon-based quantum dots were intraperitoneally injected in individuals of Carassius auratus gibelio. Their effects on white muscle were investigated by following their distribution and impact on the antioxidative system. The GSH level significantly increased after 1 and 3 days of exposure by, respectively, 85.3 and 25.4%. Seven days later, GSH levels were similar to control concentrations. MDA concentration rose after three days by 46.9% and remained at the same level after 7 days. Protein thiol levels significantly decreased by 6.7 and 8.1% after 3 and 7 days, whereas advanced oxidation protein products increased by 12.7, respectively, 28.1% in the same time intervals. The protein reactive carbonyl groups were raised only after the first day of exposure and returned to the control level later on. SOD specific activity increased up to 48% after 7 days, while CAT activity increased by 328, 176, and 26% after 1, 3, and 7 days of treatment. GST specific activity was up-regulated by 87, 18, and 9%, while GR activity increased by 68, 34, and 9%. G6PD activity was up-regulated by 12, 22, and 50%, whereas GPx activity raised by 75 and 109% compared to control after, respectively, 1, 3, and 7 days. Our results suggest that oxidative stress induced by silicon-based quantum dots was not strong enough to cause permanent damage in the white muscle of crucian carp.

Intestinal synthesis and secretion of bile salts as an adaptation to developmental biliary atresia in the sea lamprey

Bile salt synthesis is a specialized liver function in vertebrates. Bile salts play diverse roles in digestion and signaling, and their homeostasis is maintained by controlling input (biosynthesis) and intestinal conservation. Patients with biliary atresia (i.e., obliteration of the biliary tree) suffer liver fibrosis and cirrhosis. In contrast, sea lamprey thrives despite developmental biliary atresia. We discovered that the sea lamprey adapts to biliary atresia through a unique mechanism of de novo synthesis and secretion of bile salts in intestine after developmental biliary atresia, in addition to known mechanisms, such as the reduction of bile salt synthesis in liver. During and after developmental biliary atresia, expression of cyp7a1 in intestine increased by more than 100-fold (P < 0.001), whereas in liver it decreased by the same magnitude (P < 0.001). Concurrently, bile salt pools changed in similar patterns and magnitudes in these two organs and the composition shifted from C24 bile alcohol sulfates to taurine-conjugated C24 bile acids. In addition, both in vivo and ex vivo experiments showed that aductular sea lamprey secreted taurocholic acid into its intestinal lumen. Our results indicate that the sea lamprey, a jawless vertebrate, may be in an evolutionarily transitional state where bile salt synthesis occurs in both liver and intestine. Understanding the molecular basis of these mechanisms may shed light on the evolution of bile salt synthesis and possible therapy for infant biliary atresia.

Applications for next-generation sequencing in fish ecotoxicogenomics

The new technologies for next-generation sequencing (NGS) and global gene expression analyses that are widely used in molecular medicine are increasingly applied to the field of fish biology. This has facilitated new directions to address research areas that could not be previously considered due to the lack of molecular information for ecologically relevant species. Over the past decade, the cost of NGS has decreased significantly, making it possible to use non-model fish species to investigate emerging environmental issues. NGS technologies have permitted researchers to obtain large amounts of raw data in short periods of time. There have also been significant improvements in bioinformatics to assemble the sequences and annotate the genes, thus facilitating the management of these large datasets.The combination of DNA sequencing and bioinformatics has improved our abilities to design custom microarrays and study the genome and transcriptome of a wide variety of organisms. Despite the promising results obtained using these techniques in fish studies, NGS technologies are currently underused in ecotoxicogenomics and few studies have employed these methods. These issues should be addressed in order to exploit the full potential of NGS in ecotoxicological studies and expand our understanding of the biology of non-model organisms.

Genome evolution and meiotic maps by massively parallel DNA sequencing: spotted gar, an outgroup for the teleost genome duplication

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.

Common mechanism underlies repeated evolution of extreme pollution tolerance

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.

Morphogenesis of the zebrafish jaw: development beyond the embryo

The zebrafish has emerged as an important model for vertebrate development as it relates to human health and disease. Work in this system has provided significant insights into the variety of genetic signals that direct the cellular activities and tissue interactions necessary for proper assembly of the pharyngeal skeleton. Unfortunately our understanding of craniofacial development beyond embryonic stages is far less complete. Stated another way, we know a great deal about the early patterning of the skull, but we know comparatively little about how mature craniofacial shape is determined and maintained over time. Here we propose ways to expand the current molecular genetic paradigm beyond the embryo to gain an understanding of the processes and mechanisms that guide growth and remodeling of mineralized craniofacial, skeletal, and dental tissues. First, we discuss sources of adult mutant phenotypes that can be used to study of postembryonic development. Next, we review salient quantitative methods that are necessary to define complex adult phenotypes. We also discuss how other organismal systems can be used to inform and complement studies in zebrafish. We conclude by discussing the implications for such studies within the context of furthering an understanding of the etiology and pathophysiology of human craniofacial malformations, as well as informing an understanding of adaptive craniofacial variation among natural populations.

Genome enablement of the notothenioidei: genome size estimates from 11 species and BAC libraries from 2 representative taxa

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species flock in the cold Southern Ocean surrounding Antarctica. To enable genome-level studies of these psychrophilic fishes, we estimated the sizes of the genomes of 11 Antarctic species and generated high-quality BAC libraries for 2, the notothen Notothenia coriiceps and the icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families, [e.g., the icefishes (Channichthyidae)], was accompanied by genome expansion. Species (n=6) of the basal family Nototheniidae had C values that ranged between 0.98 and 1.20 pg, whereas those of the icefishes, the notothenioid crown group, were 1.66-1.83 pg (n=4 species). The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprised 12X and 10X coverage of the respective genomes and had average insert sizes of 138 and 168 kb. Greater than 60% of paired BAC ends sampled from each library ( approximately 0.1% of each genome) contained repetitive sequences, and the repetitive element landscapes of the 2 genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The representation and depth of coverage of the libraries were verified by identification of multiple Hox gene contigs: six discrete Hox clusters were found in N. coriiceps and at least five Hox clusters were found in C. aceratus. Given the unusual anatomical and physiological adaptations of the notothenioids, the availability of these BAC libraries sets the stage for expanded analysis of the psychrophilic mode of life.

Antarctic notothenioid fishes: genomic resources and strategies for analyzing an adaptive radiation

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species-flock in the cold Southern Ocean that surrounds Antarctica. To facilitate genome-level studies of the diversification of these fishes, we present estimates of the genome sizes of 11 Antarctic species and describe the production of high-quality bacterial artificial chromosome (BAC) libraries for two, the red-blooded notothen Notothenia coriiceps and the white-blooded icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families (e.g., the crown group Channichthyidae [icefishes]), was accompanied by genome expansion. Six species from the basal family Nototheniidae had C-values between 0.98 and 1.20 pg, a range that is consistent with the genome sizes of proposed outgroups (e.g., percids) of the notothenioid suborder. In contrast, four icefishes had C-values in the range 1.66-1.83 pg. The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprise 12× and 10× coverage of the respective genomes and have average insert sizes of 138 and 168 kb. Paired BAC-end reads representing ∼0.1% of each genome showed that the repetitive element landscapes of the two genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The availability of these high-quality and well-characterized BAC libraries sets the stage for targeted genomic analyses of the unusual anatomical and physiological adaptations of the notothenioids, some of which mimic human diseases. Here we consider the evolution of secondary pelagicism by various taxa of the group and illustrate the utility of Antarctic icefishes as an evolutionary-mutant model of human osteopenia (low-mineral density of bones).

Molecular pedomorphism underlies craniofacial skeletal evolution in Antarctic notothenioid fishes

Background

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.

Results

We compared craniofacial skeletal development in two pelagic notothenioids, Chaenocephalus aceratus and Pleuragramma antarcticum, to that in a benthic species, Notothenia coriiceps, and two outgroups, the threespine stickleback and the zebrafish. Relative to these other species, pelagic notothenioids exhibited a delay in pharyngeal bone development, which was associated with discrete heterochronic shifts in skeletal gene expression that were consistent with persistence of the chondrogenic program and a delay in the osteogenic program during larval development. Morphological analysis also revealed a bias toward the development of anterior and ventral elements of the notothenioid pharyngeal skeleton relative to dorsal and posterior elements.

Conclusions

Our data support the hypothesis that early shifts in the relative timing of craniofacial skeletal gene expression may have had a significant impact on the adaptive radiation of Antarctic notothenioids into pelagic habitats.