Masatoshi Nakamoto

Japanese flounder (Paralichthys olivaceus) is an important marine fish species of both fisheries and aquaculture in Northeast Asia. The commercial interest for all-female progenies due to several sex-related traits has prompted basic research on the mechanisms of sex determination in this species. By conducting a linkage analysis of the sex-determining locus, we initially identified 12 microsatellite markers linked to sex in 11 scaffolds, whose localization was restricted to a specific region of linkage group 9. Sequence analysis of this region identified 181 genes based on the UniProt database annotations. Among them, the amh gene was considered a potential candidate for sex determination because this gene is known to have taken over the role of sex determination in many teleosts. An in-depth sequence analysis of both the coding and non-coding regions of amh in XX and XY individuals detected nine SNPs linked with maleness. However, because these substitutions were synonymous, the upstream and downstream regions of amh were also investigated and a male-specific variant with deletions in the promoter region was detected. This truncated Y-specific amh variant was named amhy, and the amh shared by both sexes was named amhx. The association analysis using both females and males of the genotypic sex inferred by the presence/absence of amhy found complete association with phenotypic sex and genotype. Gene expression analysis in larvae derived from a single-pair progeny by quantitative real-time PCR detected amhy transcripts in the larval trunks between 20 and 100 days after hatching only in XY larvae. Localization of amhy by in situ hybridization was detected in presumptive Sertoli cells of XY gonads. Expression of amhx was almost undetectable in both XX and XY genotypes. Loss of Amh function by CRISPR-Cas9 induced male-to-female sex reversal, indicating that this gene was necessary for the masculinization of XY individuals. In conclusion, the complete linkage of amhy with males, its early expression in XY gonads before testicular differentiation, and the induction of sex reversal by loss-of-function mutation support the view that amhy is the sex-determining gene in this species.

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

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.