Despite the popularity of zebrafish as a research model, its sex determination (SD) mechanism is still unknown. Most cytogenetic studies failed to find dimorphic sex chromosomes and no primary sex determining switch has been identified even though the assembly of zebrafish genome sequence is near to completion and a high resolution genetic map is available. Recent publications suggest that environmental factors within the natural range have minimal impact on sex ratios of zebrafish populations. The primary aim of this study is to find out more about how sex is determined in zebrafish. Using classical breeding experiments, we found that sex ratios across families were wide ranging (4.8% to 97.3% males). On the other hand, repeated single pair crossings produced broods of very similar sex ratios, indicating that parental genotypes have a role in the sex ratio of the offspring. Variation among family sex ratios was reduced after selection for breeding pairs with predominantly male or female offspring, another indication that zebrafish sex is regulated genetically. Further examinations by a PCR-based “blind assay" and array comparative genomic hybridization both failed to find universal sex-linked differences between the male and female genomes. Together with the ability to increase the sex bias of lines by selective breeding, these data suggest that zebrafish is unlikely to utilize a chromosomal sex determination (CSD) system. Taken together, our study suggests that zebrafish sex is genetically determined with limited, secondary influences from the environment. As we have not found any sign for CSD in the species, we propose that the zebrafish has a polygenic sex determination system.

In this review, we provide a detailed overview of studies on the elusive sex determination (SD) and gonad differentiation mechanisms of zebrafish (Danio rerio). We show that the data obtained from most studies are compatible with polygenic sex determination (PSD), where the decision is made by the allelic combinations of several loci. These loci are typically dispersed throughout the genome, but in some teleost species a few of them might be located on a preferential pair of (sex) chromosomes. The PSD system has a much higher level of variation of SD genotypes both at the level of gametes and the sexual genotype of individuals, than that of the chromosomal sex determination systems. The early sexual development of zebrafish males is a complicated process, as they first develop a ‘juvenile ovary’, that later undergoes a transformation to give way to a testis. To date, three major developmental pathways were shown to be involved with gonad differentiation through the modulation of programmed cell death. In our opinion, there are more pathways participating in the regulation of zebrafish gonad differentiation/transformation. Introduction of additional powerful large-scale genomic approaches into the analysis of zebrafish reproduction will result in further deepening of our knowledge as well as identification of additional pathways and genes associated with these processes in the near future.

As complete absence of germ cells leads to sterile males in zebrafish, we explored the relationship between primordial germ cell (PGC) number and sexual development. Our results revealed dimorphic proliferation of PGCs in the early zebrafish larvae, marking the beginning of sexual differentiation. We applied morpholino-based gene knockdown and cell transplantation strategies to demonstrate that a threshold number of PGCs is required for the stability of ovarian fate. Using histology and transcriptomic analyses, we determined that zebrafish gonads are in a meiotic ovarian stage at 14 days postfertilization and identified signaling pathways supporting meiotic oocyte differentiation and eventual female fate. The development of PGC-depleted gonads appears to be restrained and delayed, suggesting that PGC number may directly regulate the variability and length of gonadal transformation and testicular differentiation in zebrafish. We propose that gonadal transformation may function as a developmental buffering mechanism to ensure the reproductive outcome.

During development all zebrafish males first develop a “juvenile ovary” that later degenerates and transforms into a testis. In this study, individuals undergoing gonadal transformation were identified from a vas::egfp transgenic line and used for gene expression analysis of anti‐Müllerian hormone (amh), ovarian aromatase (cyp19a1a) and 11β‐hydroxylase (cyp11b, also known as P450_11β) by real‐time polymerase chain reaction and in situ hybridization. In the “normal (i.e., nontransforming) juvenile ovary” cyp19a1a was expressed around the oocytes, but cyp11b and amh could not be detected. During gonadal transformation cyp19a1a was down‐regulated and could not be detected anymore; in contrast amh was up‐regulated and highly expressed at similar regions where cyp19a1a had been expressed earlier. Furthermore, the normalized transcript levels of cyp19a1a and amh showed a reciprocal picture, i.e., the higher was the level of amh, the lower was that of cyp19a1a. Expression of cyp11b was also up‐regulated but later than amh, and its localization was not related to the position of degenerating oocytes. These data indicate that amh is a candidate gene down‐regulating cyp19a1a, leading to “juvenile ovary‐to‐testis” transformation. Whereas, cyp11b or its product, 11‐ketotestosterone, is unlikely to be the inducer of zebrafish gonad transformation, as proposed earlier for some protogynous hermaphroditic fish species. Developmental Dynamics 236:1329–1338, 2007.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species’ native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics. We describe the genome assembly of Asian seabass (Lates calcarifer), a marine teleost with aquaculture relevance. Though >500 eukaryotic genome sequences are available in public repositories, the majority are highly fragmented with incomplete assemblies, which explains why considerable effort and resources are often spent to improve their quality after publication. In our study, we employed long read sequencing combined with genetic and optical mapping, and syntenic information to produce a chromosomal level assembly. The largely continuous genome assembly will be useful for comparative genomics and offers an opportunity to look into regions less explored such as tandem repeats (the core component of centromeres and telomeres). In addition, population structure of the species was analysed based on low-coverage genome sequence information from 61 individuals representing diverse geographic locations stretching from North-Western India across South-East Asia and Australia to Papua New Guinea.