A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu)

Sex Determination and Differentiation in Teleost: Roles of Genetics, Environment, and Brain

The fish reproductive system is a complex biological system. Nonetheless, reproductive organ development is conserved, which starts with sex determination and then sex differentiation. The sex of a teleost is determined and differentiated from bipotential primordium by genetics, environmental factors, or both. These two processes are species-specific. There are several prominent genes and environmental factors involved during sex determination and differentiation. At the cellular level, most of the sex-determining genes suppress the female pathway. For environmental factors, there are temperature, density, hypoxia, pH, and social interaction. Once the sexual fate is determined, sex differentiation takes over the gonadal developmental process. Environmental factors involve activation and suppression of various male and female pathways depending on the sexual fate. Alongside these factors, the role of the brain during sex determination and differentiation remains elusive. Nonetheless, GnRH III knockout has promoted a male sex-biased population, which shows brain involvement during sex determination. During sex differentiation, LH and FSH might not affect the gonadal differentiation, but are required for regulating sex differentiation. This review discusses the role of prominent genes, environmental factors, and the brain in sex determination and differentiation across a few teleost species.

Heterogeneous Histories of Recombination Suppression on Stickleback Sex Chromosomes

How consistent are the evolutionary trajectories of sex chromosomes shortly after they form? Insights into the evolution of recombination, differentiation, and degeneration can be provided by comparing closely related species with homologous sex chromosomes. The sex chromosomes of the threespine stickleback (Gasterosteus aculeatus) and its sister species, the Japan Sea stickleback (G. nipponicus), have been well characterized. Little is known, however, about the sex chromosomes of their congener, the blackspotted stickleback (G. wheatlandi). We used pedigrees to obtain experimentally phased whole genome sequences from blackspotted stickleback X and Y chromosomes. Using multispecies gene trees and analysis of shared duplications, we demonstrate that Chromosome 19 is the ancestral sex chromosome and that its oldest stratum evolved in the common ancestor of the genus. After the blackspotted lineage diverged, its sex chromosomes experienced independent and more extensive recombination suppression, greater X-Y differentiation, and a much higher rate of Y degeneration than the other two species. These patterns may result from a smaller effective population size in the blackspotted stickleback. A recent fusion between the ancestral blackspotted stickleback Y chromosome and Chromosome 12, which produced a neo-X and neo-Y, may have been favored by the very small size of the recombining region on the ancestral sex chromosome. We identify six strata on the ancestral and neo-sex chromosomes where recombination between the X and Y ceased at different times. These results confirm that sex chromosomes can evolve large differences within and between species over short evolutionary timescales.

Expanding the classical paradigm: what we have learnt from vertebrates about sex chromosome evolution

Until recently, the field of sex chromosome evolution has been dominated by the canonical unidirectional scenario, first developed by Muller in 1918. This model postulates that sex chromosomes emerge from autosomes by acquiring a sex-determining locus. Recombination reduction then expands outwards from this locus, to maintain its linkage with sexually antagonistic/advantageous alleles, resulting in Y or W degeneration and potentially culminating in their disappearance. Based mostly on empirical vertebrate research, we challenge and expand each conceptual step of this canonical model and present observations by numerous experts in two parts of a theme issue of Phil. Trans. R. Soc. B. We suggest that greater theoretical and empirical insights into the events at the origins of sex-determining genes (rewiring of the gonadal differentiation networks), and a better understanding of the evolutionary forces responsible for recombination suppression are required. Among others, crucial questions are: Why do sex chromosome differentiation rates and the evolution of gene dose regulatory mechanisms between male versus female heterogametic systems not follow earlier theory? Why do several lineages not have sex chromosomes? And: What are the consequences of the presence of (differentiated) sex chromosomes for individual fitness, evolvability, hybridization and diversification? We conclude that the classical scenario appears too reductionistic. Instead of being unidirectional, we show that sex chromosome evolution is more complex than previously anticipated and principally forms networks, interconnected to potentially endless outcomes with restarts, deletions and additions of new genomic material. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

A sex-linked SNP mutation in amhr2 is responsible for male differentiation in obscure puffer (Takifugu obscurus)

BACKGROUND

Anti-Mullerian hormone receptor type II (Amhr2) is a key receptor of Amh signaling in regulating gonad development. The amhr2 gene has been identified in numerous species, including a few teleost fishes. However, the roles of Amhr2 in Amh signaling in fish are poorly studied.

METHODS AND RESULTS

In this study, an amhr2 homolog from obscure puffer (Takifugu obscurus) was identified, and its molecular characteristics were systematically analyzed. Expression analysis revealed that amhr2 was highly expressed in the gonads of adult pufferfish and significantly upregulated during sex differentiation. Significantly, a sex-linked SNP site was verified in obscure puffer amhr2. Females exhibited a homozygous genotype (C/C), while males possessed a heterozygous genotype (C/G), resulting in an amino acid variation (His/Asp384) in the kinase domain of Amhr2. Then, the functions of the different Amhr2 genotypes were further investigated. The male genotype protein (Amhr2^D384) showed an enhanced ability to interact with the type I receptor (Bmpr1a) compared to the female genotype (Amhr2^H384). The phosphorylation levels of Smads and activity of the target gene (id3) induced by the male genotype were also much higher than those induced by the female genotype. These results confirmed that the male genotype had an enhanced effect on the Amh signaling pathway compared with the female genotype.

CONCLUSIONS

This study provides direct experimental evidence for the roles of different Amhr2 genotypes in pufferfish and suggests that amhr2 is responsible for male sex differentiation in obscure puffer.

Characterization of a Sex-Determining Region and Its Genomic Context via Statistical Estimates of Haplotype Frequencies in Daughters and Sons Sequenced in Pools

Sex chromosomes are generally derived from a pair of autosomes that have acquired a locus controlling sex. Sex chromosomes may evolve reduced recombination around this locus and undergo a long process of molecular divergence. At that point, the original loci controlling sex may be difficult to pinpoint. This difficulty has affected many model species from mammals to birds to flies, which present highly diverged sex chromosomes. Identifying sex-controlling loci is easier in species with molecularly similar sex chromosomes. Here we aimed at pinpointing the sex-determining region (SDR) of Armadillidium vulgare, a terrestrial isopod with female heterogamety (ZW females and ZZ males) and whose sex chromosomes appear to show low genetic divergence. To locate the SDR, we assessed single-nucleotide polymorphism (SNP) allele frequencies in F1 daughters and sons sequenced in pools (pool-seq) in several families. We developed a Bayesian method that uses the SNP genotypes of individually sequenced parents and pool-seq data from F1 siblings to estimate the genetic distance between a given genomic region (contig) and the SDR. This allowed us to assign more than 43 Mb of contigs to sex chromosomes, and to demonstrate extensive recombination and very low divergence between these chromosomes. By taking advantage of multiple F1 families, we delineated a very short genomic region (∼65 kb) that presented no evidence of recombination with the SDR. In this short genomic region, the comparison of sequencing depths between sexes highlighted female-specific genes that have undergone recent duplication, and which may be involved in sex determination in A. vulgare.

Seadragon genome analysis provides insights into its phenotype and sex determination locus

The iconic phenotype of seadragons includes leaf-like appendages, a toothless tubular mouth, and male pregnancy involving incubation of fertilized eggs on an open "brood patch." We de novo-sequenced male and female genomes of the common seadragon (Phyllopteryx taeniolatus) and its closely related species, the alligator pipefish (Syngnathoides biaculeatus). Transcription profiles from an evolutionary novelty, the leaf-like appendages, show that a set of genes typically involved in fin development have been co-opted as well as an enrichment of transcripts for potential tissue repair and immune defense genes. The zebrafish mutants for scpp5, which is lost in all syngnathids, were found to lack or have deformed pharyngeal teeth, supporting the hypothesis that the loss of scpp5 has contributed to the loss of teeth in syngnathids. A putative sex-determining locus encoding a male-specific amhr2y gene shared by common seadragon and alligator pipefish was identified.

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis

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.

A duplicated amh is the master sex-determining gene for Sebastes rockfish in the Northwest Pacific

Teleost fish are the most diverse group of vertebrates and provide opportunities to study the evolution of sex determination (SD) systems. Using genomic and functional analyses, we identified a male-specific duplication of anti-Müllerian hormone (amh) gene as the male master sex-determining (MSD) gene in Sebastes schlegelii. By resequencing 10 males and 10 females, we characterized a 5 kb-long fragment in HiC_Scaffold_12 as a male-specific region, which contained an amh gene (named amhy). We then demonstrated that amhy is a duplication of autosomal amh that was later translocated to the ancestral Y chromosome. amha and amhy shared high-nucleotide identity with the most significant difference being two insertions in intron 4 of amhy. Furthermore, amhy overexpression triggered female-to-male sex reversal in S. schlegelii, displaying its fundamental role in driving testis differentiation. We developed a PCR assay which successfully identified sexes in two species of northwest Pacific rockfish related to S. schlegelii. However, the PCR assay failed to distinguish the sexes in a separate clade of northeast Pacific rockfish. Our study provides new examples of amh as the MSD in fish and sheds light on the convergent evolution of amh duplication as the driving force of sex determination in different fish taxa.

The replaceable master of sex determination: bottom-up hypothesis revisited

Different group of vertebrates and invertebrates demonstrate an amazing diversity of gene regulations not only at the top but also at the bottom of the sex determination genetic network. As early as 1995, based on emerging findings in Drosophila melanogaster and Caenorhabditis elegans, Wilkins suggested that the evolution of the sex determination pathway evolved from the bottom to the top of the hierarchy. Based on our current knowledge, this review revisits the 'bottom-up' hypothesis and applies its logic to vertebrates. The basic operation of the determination network is through the dynamics of the opposing male and female pathways together with a persistent need to maintain the sexual identity of the cells of the gonad up to the reproductive stage in adults. The sex-determining trigger circumstantially acts from outside the genetic network, but the regulatory network is not built around it as a main node, thus maintaining the genetic structure of the network. New sex-promoting genes arise either through allelic diversification or gene duplication and act specially at the sex-determination period, without integration into the complete network. Due to this peripheral position the new regulator is not an indispensable component of the sex-determining network and can be easily replaced. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Epigenetic mechanisms in sex determination and in the evolutionary transitions between sexual systems

The hypothesis that epigenetic mechanisms of gene expression regulation have two main roles in vertebrate sex is presented. First, and within a given generation, by contributing to the acquisition and maintenance of (i) the male or female function once during the lifetime in individuals of gonochoristic species; and (ii) the male and female function in the same individual, either at the same time in simultaneous hermaphrodites, or first as one sex and then as the other in sequential hermaphrodites. Second, if environmental conditions change, epigenetic mechanisms may have also a role across generations, by providing the necessary phenotypic plasticity to facilitate the transition: (i) from one sexual system to another, or (ii) from one sex-determining mechanism to another. Furthermore, if the environmental change lasts enough time, epimutations could facilitate assimilation into genetic changes that stabilize the new sexual system or sex-determining mechanism. Examples supporting these assertions are presented, caveats or difficulties and knowledge gaps identified, and possible ways to test this hypothesis suggested.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Lessons from an unusual vertebrate sex-determining gene

So far, very few sex-determining genes have been identified in vertebrates and most of them, the so-called 'usual suspects', evolved from genes which fulfil essential functions during sexual development and are thus already tightly linked to the process that they now govern. The single exception to this 'usual suspects' rule in vertebrates so far is the conserved salmonid sex-determining gene, sdY (sexually dimorphic on the Y chromosome), that evolved from a gene known to be involved in regulation of the immune response. It is contained in a jumping sex locus that has been transposed or translocated into different ancestral autosomes during the evolution of salmonids. This special feature of sdY, i.e. being inserted in a 'jumping sex locus', could explain how salmonid sex chromosomes remain young and undifferentiated to escape degeneration. Recent knowledge on the mechanism of action of sdY demonstrates that it triggers its sex-determining action by deregulating oestrogen synthesis that is a conserved and crucial pathway for ovarian differentiation in vertebrates. This result suggests that sdY has evolved to cope with a pre-existing sex differentiation regulatory network. Therefore, 'limited options' for the emergence of new master sex-determining genes could be more constrained by their need to tightly interact with a conserved sex differentiation regulatory network rather than by being themselves 'usual suspects', already inside this sex regulatory network. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Evolution of master sex determiners: TGF-β signalling pathways at regulatory crossroads

To date, more than 20 different vertebrate master sex-determining genes have been identified on different sex chromosomes of mammals, birds, frogs and fish. Interestingly, six of these genes are transcription factors (Dmrt1- or Sox3- related) and 13 others belong to the TGF-β signalling pathway (Amh, Amhr2, Bmpr1b, Gsdf and Gdf6). This pattern suggests that only a limited group of factors/signalling pathways are prone to become top regulators again and again. Although being clearly a subordinate member of the sex-regulatory network in mammals, the TGF-β signalling pathway made it to the top recurrently and independently. Facing this rolling wave of TGF-β signalling pathways, this review will decipher how the TGF-β signalling pathways cope with the canonical sex gene regulatory network and challenge the current evolutionary concepts accounting for the diversity of sex-determining mechanisms. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

A novel c.1759T>G variant in follicle-stimulating hormone-receptor gene is concordant with male determination in the flathead grey mullet (Mugil cephalus)

Various master key regulators (MKRs) that control a binary switch of sex determination (SD) have been found in fish; these provide an excellent model for the study of vertebrate genetic SD. The SD region in flathead grey mullet has been previously mapped to a 1 Mbp region harboring 27 genes, of which one is follicle-stimulating hormone receptor (fshr). Although this gene is involved in gonad differentiation and function, it has not been considered as an MKR of SD. We systematically investigated polymorphism in mullet fshr using DNA shotgun sequences, and compared them between males and females. Capable of encoding nonconservative amino acid substitutions, c.1732G>A and c.1759T>G exhibited association with sex on a population level (N = 83; P ≤ 6.7 × 10-19). Hence, 1732 A and 1759 G represent a male-specific haplotype of the gene, designated as "fshry." Additional flanking SNPs showed a weaker degree of association with sex, delimiting the SD critical region to 143 nucleotides on exon 14. Lack of homozygotes for fshry, and the resulting divergence from Hardy-Weinberg equilibrium (N = 170; P ≤ 3.9 × 10-5), were compatible with a male heterogametic model (XY/XX). Capable of replacing a phenylalanine with valine, c.1759T>G alters a conserved position across the sixth transmembrane domain of vertebrate FSHRs. Amino acid substitutions in this position in vertebrates are frequently associated with constant receptor activation and consequently with FSH/FSHR signaling alteration; thus, indicating a potential role of fshr as an MKR of SD.

Might Gene Duplication and Neofunctionalization Contribute to the Sexual Lability Observed in Fish?

Sex determination and differentiation varies widely across vertebrates, but is most dramatically diverse in fishes. Among fishes sex reversal and sex change are observed in 41 teleost families spanning 7 orders. These sex-changing fish perhaps highlight better than any other system that sex determination is not the narrow and fixed construct we once thought, but a plastic trait that is better viewed as a reaction norm. However, while this stunning transformation is increasingly understood, a fundamental question arises, which is why some fish species have retained this inherent plasticity in sexual fate, while others have not? Here, we explore our current understanding of sex change in fish, some of the factors that permit and constrain sex reversal, and posit that gene duplication and neofunctionalization contribute to the sexual lability observed in fish.

No evidence for genetic sex determination in Daphnia magna

Mechanisms of sex determination (SD) differ widely across the tree of life. In genotypic sex determination (GSD), genetic elements determine whether individuals are male or female, while in environmental sex determination (ESD), external cues control the sex of the offspring. In cyclical parthenogens, females produce mostly asexual daughters, but environmental stimuli such as crowding, temperature or photoperiod may cause them to produce sons. In aphids, sons are induced by ESD, even though GSD is present, with females carrying two X chromosomes and males only one (X0 SD system). By contrast, although ESD exists in Daphnia, the two sexes were suggested to be genetically identical, based on a 1972 study on Daphnia magna (2n=20) that used three allozyme markers. This study cannot, however, rule out an X0 system, as all three markers may be located on autosomes. Motivated by the life cycle similarities of Daphnia and aphids, and the absence of karyotype information for Daphnia males, we tested for GSD (homomorphic sex chromosomes and X0) systems in D. magna using a whole-genome approach by comparing males and females of three genotypes. Our results confirm the absence of haploid chromosomes or haploid genomic regions in D. magna males as well as the absence of sex-linked genomic regions and sex-specific single-nucleotide polymorphisms. Within the limitations of the three studied populations here and the methods used, we suggest that our results make the possibility of genetic differences among sexes in the widely used Daphnia model system very unlikely.

Environmental Cues and Mechanisms Underpinning Sex Change in Fish

Fishes are the only vertebrates that undergo sex change during their lifetime, but even within this group, a unique reproductive strategy is displayed by only 1.5% of the teleosts. This lability in alternating sexual fate is the result of the simultaneous suppression and activation of opposing male and female networks. Here, we provide a brief review summarizing recent advances in our understanding of the environmental cues that trigger sex change and their perception, integration, and translation into molecular cascades that convert the sex of an individual. We particularly focus on molecular events underpinning the complex behavioral and morphological transformation involved in sex change, dissecting the main molecular players and regulatory networks that shape the transformation of one sex into the opposite. We show that histological changes and molecular pathways governing gonadal reorganization are better described than the neuroendocrine basis of sex change and that, despite important advances, information is lacking for the majority of hermaphrodite species. We highlight significant gaps in our knowledge of how sex change takes place and suggest future research directions.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy

Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal

The yellowfin seabream Acanthopagrus latus is the economically most important Sparidae fish in the northern South China Sea. As euryhaline fish, they are perfect model for investigating osmoregulatory mechanisms in teleosts. Moreover, the reproductive biology of hermaphrodites has long been intriguing; however, little information is known about the molecular pathways underlying their sex change. Here, we report a chromosome level reference genome of A. latus generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The draft genome of yellowfin seabream was 806 Mb, with 732 Mb scaffolds anchored on 24 chromosomes. The contig N50 and scaffold N50 were 2.6 Mb and 30.17 Mb, respectively. The assembly is of high integrity and includes 92.23% universal single-copy orthologues based on benchmarking universal single-copy orthologs (BUSCO) analysis. A total of 19,631 protein-coding genes were functionally annotated in the reference genome. Moreover, ARRDC3 and GSTA gene families which related to osmoregulation underwent an extensive expansion in two euryhaline sparids fish genomes compared to other teleost genomes. Moreover, integrating sex-specific transcriptome analyses, several genes related to the transforming growth factor beta (TGF-β) signalling pathway involved in sex differentiation and development. This genomic resource will not only be valuable for studying the osmoregulatory mechanisms in estuarine fish and sex determination in hermaphrodite vertebrate species, but also provide useful genomic tools for facilitating breeding of the yellowfin seabream.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BACKGROUND

Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.

RESULTS

The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.

CONCLUSIONS

The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Genetic Dissection of a Precocious Phenotype in Male Tiger Pufferfish (Takifugu rubripes) using Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di)

The novel non-targeted PCR-based genotyping system, namely Genotyping by Random Amplicon Sequencing, Direct (GRAS-Di), is characterized by the simplicity in library construction and robustness against DNA degradation and is expected to facilitate advancements in genetics, in both basic and applied sciences. In this study, we tested the utility of GRAS-Di for genetic analysis in a cultured population of the tiger pufferfish Takifugu rubripes. The genetic analyses included family structure analysis, genetic map construction, and quantitative trait locus (QTL) analysis for the male precocious phenotype using a population consisting of four full-sib families derived from a genetically precocious line. An average of 4.7 million raw reads were obtained from 198 fish. Trimmed reads were mapped onto a Fugu reference genome for genotyping, and 21,938 putative single-nucleotide polymorphisms (SNPs) were obtained. These 22 K SNPs accurately resolved the sibship and parent-offspring pairs. A fine-scale linkage map (total size: 1,949 cM; average interval: 1.75 cM) was constructed from 1,423 effective SNPs, for which the allele inheritance patterns were known. QTL analysis detected a significant locus for testes weight on Chr_14 and three suggestive loci on Chr_1, Chr_8, and Chr_19. The significant QTL was shared by body length and body weight. The effect of each QTL was small (phenotypic variation explained, PVE: 3.1-5.9%), suggesting that the precociousness seen in the cultured pufferfish is polygenic. Taken together, these results indicate that GRAS-Di is a practical genotyping tool for aquaculture species and applicable for molecular breeding programs, such as marker-assisted selection and genomic selection.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Differential expression analysis and identification of sex-related genes by gonad transcriptome sequencing in estradiol-treated and non-treated Ussuri catfish Pseudobagrus ussuriensis

The Ussuri catfish (Pseudobagrus ussuriensis) has an XX/XY sex determination system but its sex determination gene(s) remain unknown. To better understand the molecular sex determination mechanism, transcriptome analysis was conducted to obtain sex-related gene expression profiles. Transcriptome analyses were made of male and female developing/differentiating gonads by high-throughput RNA sequencing, including gonads from fish given an estradiol-induced sex reversal treatment. A total of 81,569 unigenes were assembled and 39,904 were significantly matched to known unique proteins by comparison with public databases. Twenty specifically expressed and 142 differentially expressed sex-related genes were extracted from annotated data by comparing the treatment groups. These genes are involved in spermatogenesis (e.g., Dnali1, nectin3, klhl10, mybl1, Katnal1, Eno4, Mns1, Spag6, Tsga10, Septin7), oogenesis (e.g., Lagr5, Fmn2, Npm2, zar1, Fbxo5, Fbxo43, Prdx4, Nrip1, Lfng, Atrip), gonadal development/differentiation (e.g., Cxcr4b, Hmgb2, Cftr, Ch25h, brip1, Prdm9, Tdrd1, Star, dmrt1, Tut4, Hsd17b12a, gdf9, dnd, arf1, Spata22), and estradiol response (e.g., Mmp14, Lhcgr, vtg1, vtg2, esr2b, Piwil1, Aifm1, Hsf1, gdf9). Dmrt1 and gdf9 may play an essential role in sex determination in P. ussuriensis. The expression patterns of six random genes were validated by quantitative real-time PCR, which confirmed the reliability and accuracy of the RNA-seq results. These data provide a valuable resource for future studies of gene expression and for understanding the molecular mechanism of sex determination/differentiation and gonadal development/differentiation (including hormone-induced sexual reversal) in Ussuri catfish. This has the potential to assist in producing monosex Ussuri catfish to increase aquacultural productivity.

Whole-Genome Re-sequencing and Transcriptome Reveal Oogenesis-Related Genes in Autotetraploid Carassius auratus

Oogenesis involves a series of biochemical and physiological transformations and numerous regulated genes. The autotetraploid Carassius auratus (4nRR) originated from whole-genome duplication of Carassius auratus red var. (RCC), which produces diploid eggs through pairing of diploid-like chromosome during female meiosis. To explore the molecular mechanisms underlying oogenesis in 4nRR, we used the Illumina sequencing platform to characterize the ovaries of 4nRR and RCC. Transcriptome and whole-genome re-sequencing were performed to uncover the key genes and potential genetic mutations related to oogenesis. Each sample produced paired-end reads in the range of 66.97 to 98.36 million via Illumina HiSeq™ 2500. After comparing of the transcriptome profiles between the 4nRR and RCC, we uncovered 8562 differentially expressed genes (DEGs). The DEGs were enriched in oogenesis-related processes, including oogenesis, oocyte development, ubiquitin-mediated proteolysis, the signaling pathways of MAPK and calcium, and oocyte meiosis as investigated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Additionally, whole-genome re-sequencing revealed 34,058,834 SNPs and 6,153,711 InDels, including 6,677,638 non-synonymous variations (SNPs) and 706,210 frame-shift InDels in the 8510 DEGs of 4nRR fish. Subsequently, whole-genome re-sequencing and transcriptomatic analyses revealed the genes that participate in oogenesis associated processes. Specifically, genes involved in ubiquitin-mediated proteolysis (SMURF1, UBE2I), calcium transport (CALM3, CAMK4), and meiosis (MAPK3, GRB2, CPEB1, CCNB2, YWHAE) were related to oogenesis in 4nRR. These findings enrich our understanding of oogenesis in the autopolyploid fish.

Long-term persistence of supernumerary B chromosomes in multiple species of Astyanax fish

BACKGROUND

Eukaryote genomes frequently harbor supernumerary B chromosomes in addition to the "standard" A chromosome set. B chromosomes are thought to arise as byproducts of genome rearrangements and have mostly been considered intraspecific oddities. However, their evolutionary transcendence beyond species level has remained untested.

RESULTS

Here we reveal that the large metacentric B chromosomes reported in several fish species of the genus Astyanax arose in a common ancestor at least 4 million years ago. We generated transcriptomes of A. scabripinnis and A. paranae 0B and 1B individuals and used these assemblies as a reference for mapping all gDNA and RNA libraries to quantify coverage differences between B-lacking and B-carrying genomes. We show that the B chromosomes of A. scabripinnis and A. paranae share 19 protein-coding genes, of which 14 and 11 were also present in the B chromosomes of A. bockmanni and A. fasciatus, respectively. Our search for B-specific single-nucleotide polymorphisms (SNPs) identified the presence of B-derived transcripts in B-carrying ovaries, 80% of which belonged to nobox, a gene involved in oogenesis regulation. Importantly, the B chromosome nobox paralog is expressed > 30× more than the A chromosome paralog. This indicates that the normal regulation of this gene is altered in B-carrying females, which could potentially facilitate B inheritance at higher rates than Mendelian law prediction.

CONCLUSIONS

Taken together, our results demonstrate the long-term survival of B chromosomes despite their lack of regular pairing and segregation during meiosis and that they can endure episodes of population divergence leading to species formation.

When and how do sex-linked regions become sex chromosomes?

The attention given to heteromorphism and genetic degeneration of "classical sex chromosomes" (Y chromosomes in XY systems, and the W in ZW systems that were studied first and are best described) has perhaps created the impression that the absence of recombination between sex chromosomes is inevitable. I here argue that continued recombination is often to be expected, that absence of recombination is surprising and demands further study, and that the involvement of selection in reduced recombination is not yet well understood. Despite a long history of investigations of sex chromosome pairs, there is a need for more quantitative approaches to studying sex-linked regions. I describe a scheme to help understand the relationships between different properties of sex-linked regions. Specifically, I focus on their sizes (differentiating between small regions and extensive fully sex-linked ones), the times when they evolved, and their differentiation, and review studies using DNA sequencing in nonmodel organisms that are providing information about the processes causing these properties.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable 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 specifically designed to search for and visualize 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 the 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-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic 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 nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

Molecular characterization and expression analysis of anti-Müllerian hormone in common carp (Cyprinus carpio) adult testes

Anti-Müllerian hormone (Amh) is a member of the transforming growth factor-β (Tgf-β) superfamily required in the regression of Müllerian ducts during gonadal sex differentiation of higher vertebrates. Teleost fish lack Müllerian ducts, but identified Amh orthologs have been shown to exert crucial functions during sex determination and differentiation of several species of teleosts. However, the function of Amh during gametogenesis in adult fish remains poorly investigated. Therefore, to expand present knowledge on the role of Amh in teleosts, the present study aimed to isolate and clone full-length amh cDNA in the common carp, Cyprinus carpio, and examine its expression levels throughout the male reproductive cycle and in response to different hormone treatments of testicular explants. Molecular cloning and characterization showed that the common carp Amh precursor amino acid sequence shared common features to other fish Amh precursors, including a conserved C-terminus (Tgf-β domain) and a double proteolytic cleavage site (R-X-X-R-X-X-R) upstream to the Tgf-β domain. Expression analysis showed amh dimorphic expression in the adult gonads with higher expression in the testes than ovaries. In testes, amh mRNA was detected in Sertoli cells contacting different types of germ cells, although the expression was greatest in Sertoli cells associated with type A undifferentiated spermatogonia. Expression analysis during the reproductive cycle showed that amh transcripts were down-regulated during the developing phase, which is characterized by an increased proliferation of type A undifferentiated spermatogonia and Sertoli cells and appearance of spermatocytes (meiosis) in the testes. Furthermore, ex vivo experiments showed that a 7 day exposure to Fsh or estrogens was required to decrease amh mRNA levels in common carp testicular explants. In summary, this study provided information on the molecular characterization and transcript abundance of amh in common carp adult testes. Altogether, these data will be useful for further investigations on sex determination and differentiation in this species, and also to improved strategies for improved carp aquaculture, such as inhibiting precocious maturation of males.

The rise and fall of the ancient northern pike master sex-determining gene

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

Establishment of a new sex-determining allele driven by sexually antagonistic selection

The turnover of sex-determining loci has repeatedly occurred in a number of species, rather than having a diverged pair of sex chromosomes. We model the turnover process by considering a linked locus under sexually antagonistic selection. The entire process of a turnover may be divided into two phases, which are referred to as the stochastic and deterministic phases. The stochastic phase is when a new sex-determining allele just arises and is still rare and random genetic drift plays an important role. In the deterministic phase, the new allele further increases in frequency by positive selection. The theoretical results currently available are for the deterministic phase, which demonstrated that a turnover of a newly arisen sex-determining locus could benefit from selection at a linked locus under sexually antagonistic selection, by assuming that sexually antagonistic selection works in a form of balancing selection. In this work, we provide a comprehensive theoretical description of the entire process from the stochastic phase to the deterministic phase. In addition to balancing selection, we explore several other modes of selection on the linked locus. Our theory allows us make a quantitative argument on the rate of turnover and the effect of the mode of selection at the linked locus. We also performed simulations to explore the pattern of polymorphism around the new sex-determining locus. We find that the pattern of polymorphism is informative to infer how selection worked through the turnover process.

Effects of the Thyroid Endocrine System on Gonadal Sex Ratios and Sex-Related Gene Expression in the Pufferfish Takifugu rubripes

To examine the effect and mechanism of thyroid hormone on gonadal sex differentiation, Takifugu rubripes larvae were treated with goitrogen (methimazole, MET, 1000 g/g), and thyroxine (T4, 2nM) from 25 to 80 days after hatching (dah). Gonadal histology and sex ratios of fish were then determined at 80 dah. MET treatment induced masculinization, but T4 treatment did not induce feminization in T. rubripes larvae. Transcriptomic analysis of gonads at 80 dah was then conducted. Among the large number of differentially expressed genes between the groups, the expression of foxl2, cyp19a1a, and dmrt1 was altered. The expression of foxl2, cyp19a1a, dmrt1 and gsdf at 25, 40, 55 days after treatment (dat) was further analyzed by qPCR. MET treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 in genetic females (p < 0.05). Additionally, T4 treatment induced an increase in the expression of cyp19a1a in genetic XY gonads only at 25 dat. However, the increase in cyp19a1a expression did not continue to 40 and 55 dat. This may explain why feminization of larvae was not found in the T4-treated group. Thus, the present study provides the first evidence that MET treatment causes masculinization in teleost fish. The effects of MET-induced masculinization in T. rubripes may act primarily via suppression of the expression of foxl2 and cyp19a1a, and stimulation of the expression of dmrt1. Moreover, the effects of higher concentrations of T4 or different concentrations of T3, on sex differentiation require further testing.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Disagreement in FST estimators: A case study from sex chromosomes

Sewall Wright developed FST for describing population differentiation and it has since been extended to many novel applications, including the detection of homomorphic sex chromosomes. However, there has been confusion regarding the expected estimate of FST for a fixed difference between the X- and Y-chromosome when comparing males and females. Here, we attempt to resolve this confusion by contrasting two common FST estimators and explain why they yield different estimates when applied to the case of sex chromosomes. We show that this difference is true for many allele frequencies, but the situation characterized by fixed differences between the X- and Y-chromosome is among the most extreme. To avoid additional confusion, we recommend that all authors using FST clearly state which estimator of FST their work uses.

A frog with three sex chromosomes that co-mingle together in nature: Xenopus tropicalis has a degenerate W and a Y that evolved from a Z chromosome

In many species, sexual differentiation is a vital prelude to reproduction, and disruption of this process can have severe fitness effects, including sterility. It is thus interesting that genetic systems governing sexual differentiation vary among-and even within-species. To understand these systems more, we investigated a rare example of a frog with three sex chromosomes: the Western clawed frog, Xenopus tropicalis. We demonstrate that natural populations from the western and eastern edges of Ghana have a young Y chromosome, and that a male-determining factor on this Y chromosome is in a very similar genomic location as a previously known female-determining factor on the W chromosome. Nucleotide polymorphism of expressed transcripts suggests genetic degeneration on the W chromosome, emergence of a new Y chromosome from an ancestral Z chromosome, and natural co-mingling of the W, Z, and Y chromosomes in the same population. Compared to the rest of the genome, a small sex-associated portion of the sex chromosomes has a 50-fold enrichment of transcripts with male-biased expression during early gonadal differentiation. Additionally, X. tropicalis has sex-differences in the rates and genomic locations of recombination events during gametogenesis that are similar to at least two other Xenopus species, which suggests that sex differences in recombination are genus-wide. These findings are consistent with theoretical expectations associated with recombination suppression on sex chromosomes, demonstrate that several characteristics of old and established sex chromosomes (e.g., nucleotide divergence, sex biased expression) can arise well before sex chromosomes become cytogenetically distinguished, and show how these characteristics can have lingering consequences that are carried forward through sex chromosome turnovers.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

Transcriptome-wide single nucleotide polymorphisms related to electric organ discharge differentiation among African weakly electric fish species

African weakly electric fish of the mormyrid genus Campylomormyrus generate pulse-type electric organ discharges (EODs) for orientation and communication. Their pulse durations are species-specific and elongated EODs are a derived trait. So far, differential gene expression among tissue-specific transcriptomes across species with different pulses and point mutations in single ion channel genes indicate a relation of pulse duration and electrocyte geometry/excitability. However, a comprehensive assessment of expressed Single Nucleotide Polymorphisms (SNPs) throughout the entire transcriptome of African weakly electric fish, with the potential to identify further genes influencing EOD duration, is still lacking. This is of particular value, as discharge duration is likely based on multiple cellular mechanisms and various genes. Here we provide the first transcriptome-wide SNP analysis of African weakly electric fish species (genus Campylomormyrus) differing by EOD duration to identify candidate genes and cellular mechanisms potentially involved in the determination of an elongated discharge of C. tshokwe. Non-synonymous substitutions specific to C. tshokwe were found in 27 candidate genes with inferred positive selection among Campylomormyrus species. These candidate genes had mainly functions linked to transcriptional regulation, cell proliferation and cell differentiation. Further, by comparing gene annotations between C. compressirostris (ancestral short EOD) and C. tshokwe (derived elongated EOD), we identified 27 GO terms and 2 KEGG pathway categories for which C. tshokwe significantly more frequently exhibited a species-specific expressed substitution than C. compressirostris. The results indicate that transcriptional regulation as well cell proliferation and differentiation take part in the determination of elongated pulse durations in C. tshokwe. Those cellular processes are pivotal for tissue morphogenesis and might determine the shape of electric organs supporting the observed correlation between electrocyte geometry/tissue structure and discharge duration. The inferred expressed SNPs and their functional implications are a valuable resource for future investigations on EOD durations.

Locating the Sex Determining Region of Linkage Group 12 of Guppy (Poecilia reticulata)

Despite over 100 years of study, the location of the fully sex-linked region of the guppy (Poecilia reticulata) carrying the male-determining locus, and the regions where the XY pair recombine, remain unclear. Previous population genomics studies to determine these regions used small samples from recently bottlenecked captive populations, which increase the false positive rate of associations between individuals' sexes and SNPs. Using new data from multiple natural populations, we show that a recently proposed candidate for this species' male-determining gene is probably not completely sex-linked, leaving the maleness factor still unidentified. Variants in the chromosome 12 region carrying the candidate gene sometimes show linkage disequilibrium with the sex-determining factor, but no consistently male-specific variant has yet been found. Our genetic mapping with molecular markers spread across chromosome 12 confirms that this is the guppy XY pair. We describe two families with recombinants between the X and Y chromosomes, which confirm that the male-determining locus is in the region identified by all previous studies, near the terminal pseudo-autosomal region (PAR), which crosses over at a very high rate in males. We correct the PAR marker order, and assign two unplaced scaffolds to the PAR. We also detect a duplication, with one copy in the male-determining region, explaining signals of sex linkage in a more proximal region.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction

Background

Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes.

Results

We generated 96.7 × 10⁹ 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 10⁹ paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females.

Conclusions

The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y.

Preferential Mapping of Sex-Biased Differentially-Expressed Genes of Larvae to the Sex-Determining Region of Flathead Grey Mullet (Mugil cephalus)

Flathead gray mullet (Mugil cephalus) is a cosmopolitan mugilid species popular in fishery and aquaculture with an economic preference for all-female population. However, it displays neither sexual dimorphisms nor heteromorphic sex chromosomes. We have previously presented a microsatellite-based linkage map for this species locating a single sex determination region (SDR) on linkage group 9 (LG9) with evidence for XX/XY sex determination (SD) mechanism. In this work, we refine the critical SDR on LG9, and propose positional- and functional- candidate genes for SD. To elucidate the genetic mechanism of SD, we assembled and compared male and female genomic sequences of 19 syntenic genes within the putative SDR on mullet's LG9, based on orthology to tilapia's LG8 (tLG8) physical map. A total of 25 sequence-based markers in 12 genes were developed. For all markers, we observed association with sex in at least one of the two analyzed M. cephalus full-sib families, but not in the wild-type population. Recombination events were inferred within families thus setting the SDR boundaries to a region orthologous to ∼0.9 Mbp with 27 genes on tLG8. As the sexual phenotype is evident only in adults, larvae were assigned into two putative sex-groups according to their paternal haplotypes, following a model of XY/XX SD-system. A total of 107 sex-biased differentially expressed genes in larvae were observed, of which 51 were mapped to tLG8 (48% enrichment), as compared to 5% in random control. Furthermore, 23 of the 107 genes displayed sex-specific expression; and 22 of these genes were positioned to tLG8, indicating 96% enrichment. Of the 27 SDR genes, BCCIP, DHX32A, DOCK1, and FSHR (GTH-RI) are suggested as positional and functional gene candidates for SD.

Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758)

BACKGROUND

Goldfish is an important model for various areas of research, including neural development and behavior and a species of significant importance in aquaculture, especially as an ornamental species. It has a male heterogametic (XX/XY) sex determination system that relies on both genetic and environmental factors, with high temperatures being able to produce female-to-male sex reversal. Little, however, is currently known on the molecular basis of genetic sex determination in this important cyprinid model. Here we used sequencing approaches to better characterize sex determination and sex-chromosomes in an experimental strain of goldfish.

RESULTS

Our results confirmed that sex determination in goldfish is a mix of environmental and genetic factors and that its sex determination system is male heterogametic (XX/XY). Using reduced representation (RAD-seq) and whole genome (pool-seq) approaches, we characterized sex-linked polymorphisms and developed male specific genetic markers. These male specific markers were used to distinguish sex-reversed XX neomales from XY males and to demonstrate that XX female-to-male sex reversal could even occur at a relatively low rearing temperature (18 °C), for which sex reversal has been previously shown to be close to zero. We also characterized a relatively large non-recombining region (~ 11.7 Mb) on goldfish linkage group 22 (LG22) that contained a high-density of male-biased genetic polymorphisms. This large LG22 region harbors 373 genes, including a single candidate as a potential master sex gene, i.e., the anti-Mullerian hormone gene (amh). However, no sex-linked polymorphisms were detected in the coding DNA sequence of the goldfish amh gene.

CONCLUSIONS

These results show that our goldfish strain has a relatively large sex locus on LG22, which is likely the Y chromosome of this experimental population. The presence of a few XX males even at low temperature also suggests that other environmental factors in addition to temperature could trigger female-to-male sex reversal. Finally, we also developed sex-linked genetic markers, which will be important tools for future research on sex determination in our experimental goldfish population. However, additional work would be needed to explore whether this sex locus is conserved in other populations of goldfish.

Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution

BACKGROUND

Heteromorphic sex chromosomes have evolved repeatedly across diverse species. Suppression of recombination between X and Y chromosomes leads to degeneration of the Y chromosome. The progression of degeneration is not well understood, as complete sequence assemblies of heteromorphic Y chromosomes have only been generated across a handful of taxa with highly degenerate sex chromosomes. Here, we describe the assembly of the threespine stickleback (Gasterosteus aculeatus) Y chromosome, which is less than 26 million years old and at an intermediate stage of degeneration. Our previous work identified that the non-recombining region between the X and the Y spans approximately 17.5 Mb on the X chromosome.

RESULTS

We combine long-read sequencing with a Hi-C-based proximity guided assembly to generate a 15.87 Mb assembly of the Y chromosome. Our assembly is concordant with cytogenetic maps and Sanger sequences of over 90 Y chromosome BAC clones. We find three evolutionary strata on the Y chromosome, consistent with the three inversions identified by our previous cytogenetic analyses. The threespine stickleback Y shows convergence with more degenerate sex chromosomes in the retention of haploinsufficient genes and the accumulation of genes with testis-biased expression, many of which are recent duplicates. However, we find no evidence for large amplicons identified in other sex chromosome systems. We also report an excellent candidate for the master sex-determination gene: a translocated copy of Amh (Amhy).

CONCLUSIONS

Together, our work shows that the evolutionary forces shaping sex chromosomes can cause relatively rapid changes in the overall genetic architecture of Y chromosomes.

Comparative transcriptome analysis on four types of gonadal tissues of blotched snakehead (Channa maculata)

Blotched snakehead (Channa maculata) is an economically important freshwater fish in China, of which males grow much faster than females. To illuminate the molecular mechanism of sex differentiation and gonad development, RNA-Sequencing was performed to identify sex-related genes and pathway in gonads of 6-month-old normal XX females (XX-F), normal XY males (XY-M), XY sex reversal females (XY-F) and YY super-males (YY-M). The analysis showed that many differentially expressed genes (DEGs) had similar expression patterns in XY-F and XX-F, which were different from XY-M and YY-M. qRT-PCR indicated that Amh, Dmrt1, and Sox9 had relatively high expression in testes of XY-M and YY-M. Taking Amh as an example, there was a relative fold change of 1.0 in XX-F, 2.1 fold change in XY-F, 36.1 fold change in XY-M, and 26.0 fold change in YY-M. Cyp19a1a, Figla, and Foxl2 were highly expressive in ovaries of XX-F and XY-F. Taking Figla as an example, there was a relative fold change of 557 in XX-F, 304.5 fold change in XY-F, 5.6 fold change in XY-M, and 4.4 fold change in YY-M. KEGG analysis revealed many DEGs distributed in pathways related to sex differentiation, steroid hormone synthesis and growth, etc. Significant variation and trends in relative expression levels tested by qRT-PCR were consistent with those recorded by RNA-Sequencing. This is the first time that transcriptome of snakehead has been investigated systematically and in an integrated way. Large quantities of candidate genes involved in sex differentiation, gonad development and growth dimorphism were identified. The study provides useful resources for understanding sex differentiation and growth dimorphism, potentially assisting mono-sex production of snakehead in aquaculture.

Gonadal Transcriptome Analysis of Sex-Related Genes in the Protandrous Yellowfin Seabream (Acanthopagrus latus)

Yellowfin seabream (Acanthopagrus latus), a protandrous hermaphroditic fish, is a good model for studying the mechanism of sex reversal. However, limited knowledge is known about the genetic information related to reproduction and sex differentiation in this species. Here, we performed de novo transcriptome sequencing analysis of the testis, ovotestis, and ovary to identify sex-related genes in yellowfin seabream. The results assembled 71,765 unigenes in which 16,126 and 17,560 unigenes were differentially expressed in the ovotestis and ovary compared to the testis, respectively. The most differentially expressed gene (DEG)-enriched Kyoto Encyclopedia of Genes and Genomes and GO pathways were closely associated with the synthesis of sex steroid hormones. Functional analyses identified 55 important sex-related DEGs, including 32 testis-biased DEGs (dmrt1, amh, and sox9, etc.), 20 ovary-biased DEGs (cyp19a, foxl2, and wnt4, etc.), and 3 ovotestis-biased DEGs (lhb, dmrt2, and foxh1). Furthermore, the testis-specific expression of dmrt1 and the brain-pituitary-ovary axis expression of foxl2 were characterized, suggesting that they might play important roles in sex differentiation in yellowfin seabream. Our present work provided an important molecular basis for elucidating the mechanisms underlying sexual transition and reproductional regulation in yellowfin seabream.