Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RADtag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.. . .researchers using standard lines of zebrafish that have long been maintained in laboratories are often plagued by severe sex ratio distortions. . .. C. Lawrence, J. P. Ebersole, and R. V. Kesseli (2008) C ONSIDERING the fundamental importance of sex for species propagation, it is surprising that primary sexdetermining mechanisms are not strongly conserved among animal taxa (Bull 1983;Charlesworth 1996;Ming et al. 2011;Bachtrog et al. 2014). Closely related species or even populations of the same species can have different sexdetermining mechanisms (Takehana et al. 2007;Ross et al. 2009;Kobayashi et al. 2013;Heule et al. 2014;Larney et al. 2014). Zebrafish (Danio rerio) is a popular model for studies of vertebrate development, behavior, physiology, evolution, disease, and human health (Mills et al. 2007;Seth et al. 2013;Braasch et al. 2014;Ota and Kawahara 2014;Wilkinson et al. 2014), but researchers struggle with highly variable and distorted sex ratios, and investigations into the genetic nature of zebrafish sex determination are conflicting. To help understand these issues, we conducted a population genomic study of sex determination in multiple zebrafish strains. Zebrafish exhibit...

Zebrafish (Danio rerio) is an important model for vertebrate development, genomics, physiology, behavior, toxicology, and disease. Additionally, work on numerous Danio species is elucidating evolutionary mechanisms for morphological development. Yet, the relationships of zebrafish and its closest relatives remain unclear possibly due to incomplete lineage sorting, speciation with gene flow, and interspecies hybridization. To clarify these relationships, we first constructed phylogenomic data sets from 30,801 restriction-associated DNA (RAD)-tag loci (483,026 variable positions) with clear orthology to a single location in the sequenced zebrafish genome. We then inferred a well-supported species tree for Danio and tested for gene flow during the diversification of the genus. An approach independent of the sequenced zebrafish genome verified all inferred relationships. Although identification of the sister taxon to zebrafish has been contentious, multiple RAD-tag data sets and several analytical methods provided strong evidence for Danio aesculapii as the most closely related extant zebrafish relative studied to date. Data also displayed patterns consistent with gene flow during speciation and postspeciation introgression in the lineage leading to zebrafish. The incorporation of biogeographic data with phylogenomic analyses put these relationships in a phylogeographic context and supplied additional support for D. aesculapii as the sister species to D. rerio. The clear resolution of this study establishes a framework for investigating the evolutionary biology of Danio and the heterogeneity of genome evolution in the recent history of a model organism within an emerging model genus for genetics, development, and evolution.

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.

The study of sex determination and sex chromosome organisation in non-model species has long been technically challenging, but new sequencing methodologies are now enabling precise and high-throughput identification of sex-specific genomic sequences. In particular, Restriction Site-Associated DNA Sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software designed to specifically search for sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analyzed a published dataset of Japanese medaka, Oryzias latipes , where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyze new RAD-Seq datasets from 15 fish species spanning multiple systematic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in non-model species and outperforms the commonly used RAD-Seq analysis software STACKS.RADSex in speed, resource usage, ease of application, and visualization options. Furthermore, our analysis of new datasets from 15 species provides new insights on sex determination in fish.

The study of sex determination and sex chromosome organisation in non-model species has long been technically challenging, but new sequencing methodologies are now enabling precise and high-throughput identification of sex-specific genomic sequences. In particular, Restriction Site-Associated DNA Sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software designed to specifically search for sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analyzed a published dataset of Japanese medaka, Oryzias latipes , where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyze new RAD-Seq datasets from 15 fish species spanning multiple systematic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in non-model species and outperforms the commonly used RAD-Seq analysis software STACKS.RADSex in speed, resource usage, ease of application, and visualization options. Furthermore, our analysis of new datasets from 15 species provides new insights on sex determination in fish.

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons.

Molecular and cellular mechanisms underlying variation in adult form remain largely unknown. Adult pigment patterns of fishes in the genus Danio, which includes zebrafish, D. rerio, consist of horizontal stripes, vertical bars, spots and uniform patterns, and provide an outstanding opportunity to identify causes of species level variation in a neural crest derived trait. Understanding pigment pattern variation requires quantitative approaches to assess phenotypes, yet such methods have been mostly lacking for pigment patterns. We introduce metrics derived from information theory that describe patterns and pattern variation in Danio fishes. We find that these metrics used singly and in multivariate combinations are suitable for distinguishing general pattern types, and can reveal even subtle phenotypic differences attributable to mutations. Our study provides new tools for analyzing pigment pattern in Danio and potentially other groups, and sets the stage for future analyses of pattern morphospace and its mechanistic underpinnings.

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons.Author SummaryPigment patterns of fishes are diverse and function in a wide range of behaviors. Common pattern themes include stripes and spots, exemplified by the closely related minnows Danio quagga and D. kyathit, respectively. We show that these patterns arise late in development owing to alterations in the development and arrangements of pigment cells. In the closely related model organism zebrafish (D. rerio) single genes can switch the pattern from stripes to spots. Yet, we show that pattern differences between D. quagga and D. kyathit have a more complex genetic basis, depending on multiple genes and interactions between these genes. Our findings illustrate the importance of characterizing naturally occuring genetic variants, in addition to laboratory induced mutations, for a more complete understanding of pigment pattern development and evolution.