Identification of a candidate sex determination gene in Culaea inconstans suggests convergent recruitment of an Amh duplicate in two lineages of stickleback

The fourspine stickleback (Apeltes quadracus) has an XY sex chromosome system with polymorphic inversions on both X and Y chromosomes

Teleost fish are well-known for possessing a diversity of sex chromosomes and for undergoing frequent turnovers of these sex chromosomes. However, previous studies have mainly focused on variation between species, while comparatively little attention has been given to sex chromosome polymorphisms within species, which may capture early stages of sex chromosome changes. To better understand the evolution of sex chromosomes, we used the fourspine stickleback (Apeltes quadracus) as a model organism. Previous cytogenetic studies suggested that females of this species possessed a ZW heteromorphic sex chromosome system. However, genetic crosses and our whole-genome sequencing of three geographically distinct wild populations revealed that A. quadracus has an XY sex chromosome on chromosome 23. This chromosome has not previously been identified as a sex chromosome in any other stickleback species, indicating a recent sex chromosome turnover. We also identified two genes - rxfp2a and zar1l - as novel candidate sex determination genes. Notably, we observed inversions on both the X and Y chromosomes in different populations, resulting in distinctive patterns of differentiation between the X and Y chromosomes across populations. The new sex chromosome and intraspecies inversion polymorphisms observed in A. quadracus provide an excellent system for future work assessing the relative fitness effects of the inversions, which will enable testing theoretical models about the drivers of sex chromosome evolution and turnover.

The Cause-Effect Model of Master Sex Determination Gene Acquisition and the Evolution of Sex Chromosomes

The canonical model of vertebrate sex chromosome evolution predicts a one-way trend toward degradation. However, most sex chromosomes in lower vertebrates are homomorphic. Recent progress in studies of sex determination has resulted in the discovery of more than 30 master sex determination (MSD) genes, most of which are from teleost fish. An analysis of MSD gene acquisition, recombination suppression, and sex chromosome-specific sequences revealed correlations in the modes of MSD gene acquisition and the evolution of sex chromosomes. Sex chromosomes remain homomorphic with MSD genes acquired by simple mutations, gene duplications, allelic variations, or neofunctionalization; in contrast, they become heteromorphic with MSD genes acquired by chromosomal inversion, fusion, and fission. There is no recombination suppression with sex chromosomes carrying MSD genes gained through simple mutations. In contrast, there is extensive recombination suppression with sex chromosomes carrying MSD genes gained through chromosome inversion. There is limited recombination suppression with sex chromosomes carrying MSD genes gained through transposition or translocation. We propose a cause-effect model that predicts sex chromosome evolution as a consequence of the acquisition modes of MSD genes, which explains the evolution of sex chromosomes in various vertebrates. A key factor determining the trend of sex chromosome evolution is whether non-homologous regions are created during the acquisition of MSD genes. Chromosome inversion creates inversely homologous but directly non-homologous sequences, which lead to recombination suppression but retain recombination potential. Over time, recurrent recombination in the inverted regions leads to the formation of strata and may cause the degradation of sex chromosomes. Depending on the nature of deletions in the inverted regions, sex chromosomes may evolve with dosage compensation, or the selective retention of haplo-insufficient genes may be used as an alternative strategy.

Repetitive DNAs and differentiation of the ZZ/ZW sex chromosome system in the combtail fish Belontia hasselti (Perciformes: Osphronemidae)

BACKGROUND

Java combtail fish Belontia hasselti (Cuvier, 1831), a member of the Osphronemidae family, inhabits lakes and rivers throughout Southeast Asia and Sri Lanka. Previous cytogenetic research revealed it possesses a diploid chromosome number of 48 chromosomes with a female-heterogametic ZZ/ZW sex chromosome system, where the W chromosome is distinguishable as the only metacentric element in the complement. Female-heterogametic sex chromosome systems seem to be otherwise surprisingly rare in the highly diverse order Perciformes and, therefore, B. hasselti provides an important comparative model to evolutionary studies in this teleost lineage. To examine the level of sex chromosome differentiation in B. hasselti and the contribution of repetitive DNAs to this process we combined bioinformatic analyses with chromosomal mapping of selected repetitive DNA classes, and comparative genomic hybridization.

RESULTS

By providing the first satellitome study in Perciformes, we herein identified 13 satellite DNA monomers in B. hasselti, suggesting a very low diversity of satDNA in this fish species. Using fluorescence in situ hybridization, we revealed detectable clusters on chromosomes only for four satellite DNA monomers. Together with the two mapped microsatellite motifs, the repeats primarily accumulated on autosomes, with no distinct clusters located on the sex chromosomes. Comparative genomic hybridization showed no region with accumulated female-specific or enriched repeats on the W chromosome. Telomeric repeats terminated all chromosomes, and no additional interstitial sites were detected.

CONCLUSION

These data collectively indicate a low degree of sex chromosome differentiation in B. hasselti despite their considerable heteromorphy. Possible mechanisms that may underlie this pattern are discussed.

Sex chromosome turnover and biodiversity in fishes

The impact of sex chromosomes and their turnover in speciation remains a subject of ongoing debate in the field of evolutionary biology. Fishes are the largest group of vertebrates, and they exhibit unparalleled sexual plasticity, as well as diverse sex-determining (SD) genes, sex chromosomes, and sex-determination mechanisms. This diversity is hypothesized to be associated with the frequent turnover of sex chromosomes in fishes. Although it is evident that amh and amhr2 are repeatedly and independently recruited as SD genes, their relationship with the rapid turnover of sex chromosomes and the biodiversity of fishes remains unknown. We summarize the canonical models of sex chromosome turnover and highlight the vital roles of gene mutation and hybridization with empirical evidence. We revisit Haldane's rule and the large X-effect and propose the hypothesis that sex chromosomes accelerate speciation by multiplying genotypes via hybridization. By integrating recent findings on the turnover of SD genes, sex chromosomes, and sex-determination systems in fish species, this review provides insights into the relationship between sex chromosome evolution and biodiversity in fishes.

Sex chromosome turnover in hybridizing stickleback lineages

Recent discoveries of sex chromosome diversity across the tree of life have challenged the canonical model of conserved sex chromosome evolution and evoked new theories on labile sex chromosomes that maintain less differentiation and undergo frequent turnover. However, theories of labile sex chromosome evolution lack direct empirical support due to the paucity of case studies demonstrating ongoing sex chromosome turnover in nature. Two divergent lineages (viz. WL & EL) of nine-spined sticklebacks (Pungitius pungitius) with different sex chromosomes (linkage group [LG] 12 in the EL, unknown in the WL) hybridize in a natural secondary contact zone in the Baltic Sea, providing an opportunity to study ongoing turnover between coexisting sex chromosomes. In this study, we first identify an 80 kbp genomic region on LG3 as the sex-determining region (SDR) using whole-genome resequencing data of family crosses of a WL population. We then verify this region as the SDR in most other WL populations and demonstrate a potentially ongoing sex chromosome turnover in admixed marine populations where the evolutionarily younger and homomorphic LG3 sex chromosome replaces the older and heteromorphic LG12 sex chromosome. The results provide a rare glimpse of sex chromosome turnover in the wild and indicate the possible existence of additional yet undiscovered sex chromosome diversity in Pungitius sticklebacks.

Transforming growth factor-β (TGF-β): A master signal pathway in teleost sex determination

Sex determination and differentiation in fish has always been a hot topic in genetic breeding of aquatic animals. With the advances in next-generation sequencing (NGS) in recent years, sex chromosomes and sex determining genes can be efficiently identified in teleosts. To date, master sex determination genes have been elucidated in 114 species, of which 72 species have sex determination genes belonging to TGF-β superfamily. TGF-β is the only signaling pathway that the largest proportion of components, which including ligands (amhy, gsdfy, gdf6), receptors (amhr, bmpr), and regulator (id2bby), have opportunity recognized as a sex determination gene. In this review, we focus on the recent studies about teleost sex-determination genes within TGF-β superfamily and propose several hypotheses on how these genes regulate sex determination process. Differing from other reviews, our review specifically devotes significant attention to all members of the TGF-β signal pathway, not solely the sex determination genes within the TGF-β superfamily. However, the functions of the paralogous genes of TGF superfamily are still needed ongoing research. Further studies are required to more accurately interpret the molecular mechanism of TGF-β superfamily sex determination genes.

Heterogeneous genomic architecture of skeletal armour traits in sticklebacks

Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonizing and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to 3 moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the 4 QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and 3 moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.

SLRfinder: A method to detect candidate sex-linked regions with linkage disequilibrium clustering

Despite their critical roles in genetic sex determination, sex chromosomes remain unknown in many non-model organisms, especially those having recently evolved sex-linked regions (SLRs). These evolutionarily young and labile sex chromosomes are important for understanding early sex chromosome evolution but are difficult to identify due to the lack of Y/W degeneration and SLRs limited to small genomic regions. Here, we present SLRfinder, a method to identify candidate SLRs using linkage disequilibrium (LD) clustering, heterozygosity and genetic divergence. SLRfinder does not rely on specific sequencing methods or a specific type of reference genome (e.g., from the homomorphic sex). In addition, the input of SLRfinder does not require phenotypic sexes, which may be unknown from population sampling, but sex information can be incorporated and is necessary to validate candidate SLRs. We tested SLRfinder using various published datasets and compared it to the local principal component analysis (PCA) method and the depth-based method Sex Assignment Through Coverage (SATC). As expected, the local PCA method could not be used to identify unknown SLRs. SATC works better on conserved sex chromosomes, whereas SLRfinder outperforms SATC in analysing labile sex chromosomes, especially when SLRs harbour inversions. Power analyses showed that SLRfinder worked better when sampling more populations that share the same SLR. If analysing one population, a relatively larger sample size (around 50) is needed for sufficient statistical power to detect significant SLR candidates, although true SLRs are likely always top-ranked. SLRfinder provides a novel and complementary approach for identifying SLRs and uncovering additional sex chromosome diversity in nature.

Genetic sex determination in three closely related hydrothermal vent gastropods, including one species with intersex individuals

Molluscs have undergone many transitions between separate sexes and hermaphroditism, which is of interest in studying the evolution of sex determination and differentiation. Here, we combined multi-locus genotypes obtained from restriction site-associated DNA (RAD) sequencing with anatomical observations of the gonads of three deep-sea hydrothermal vent gastropods of the genus Alviniconcha living in the southwest Pacific. We found that all three species (Alviniconcha boucheti, Alviniconcha strummeri, and Alviniconcha kojimai) share the same male-heterogametic XY sex-determination system but that the gonads of XX A. kojimai individuals are invaded by a variable proportion of male reproductive tissue. The identification of Y-specific RAD loci (found only in A. boucheti) and the phylogenetic analysis of three sex-linked loci shared by all species suggested that X-Y recombination has evolved differently within each species. This situation of three species showing variation in gonadal development around a common sex-determination system provides new insights into the reproductive mode of poorly known deep-sea species and opens up an opportunity to study the evolution of recombination suppression on sex chromosomes and its association with mixed or transitory sexual systems.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

Homeology of sex chromosomes in Amazonian Harttia armored catfishes supports the X-fission hypothesis for the X1X2Y sex chromosome system origin

The Neotropical monophyletic catfish genus Harttia represents an excellent model to study karyotype and sex chromosome evolution in teleosts. Its species split into three phylogenetic clades distributed along the Brazilian territory and they differ widely in karyotype traits, including the presence of standard or multiple sex chromosome systems in some members. Here, we investigate the chromosomal rearrangements and associated synteny blocks involved in the origin of a multiple X1X2Y sex chromosome system present in three out of six sampled Amazonian-clade species. Using 5S and 18S ribosomal DNA fluorescence in situ hybridization and whole chromosome painting with probes corresponding to X1 and X2 chromosomes of X1X2Y system from H. punctata, we confirm previous assumptions that X1X2Y sex chromosome systems of H. punctata, H. duriventris and H. villasboas represent the same linkage groups which also form the putative XY sex chromosomes of H. rondoni. The shared homeology between X1X2Y sex chromosomes suggests they might have originated once in the common ancestor of these closely related species. A joint arrangement of mapped H. punctata X1 and X2 sex chromosomes in early diverging species of different Harttia clades suggests that the X1X2Y sex chromosome system may have formed through an X chromosome fission rather than previously proposed Y-autosome fusion.

Pelvic spine reduction affects diet but not gill raker morphology in two polymorphic brook stickleback (Culaea inconstans) populations

Pelvic spine polymorphism occurs in several species in the stickleback family (Gasterosteidae). Given the similar phenotypic polymorphisms in multiple stickleback species, we sought to determine the extent of parallelism in the ecological correlates of pelvic spine reduction. Based on a metabarcoding analysis of brook stickleback gut contents in two polymorphic populations, we found that significant diet differences were associated with pelvic spine reduction, but we found no clear or consistent trend supporting a tendency for benthic feeding in pelvic-reduced brook sticklebacks. These results contrast with those found in threespine sticklebacks where pelvic spine reduction is often associated with a benthic diet. Hence, we found non-parallel consequences of spine polymorphism across species. Furthermore, a difference in gill raker morphology has been frequently observed between ecomorphs with different diets in many fish species. However, we found no evidence of any difference in gill raker morphology associated with pelvic spine polymorphism in brook sticklebacks.

Rapid turnover and evolution of sex-determining regions in Sebastes rockfishes

Nature has evolved a wealth of sex determination (SD) mechanisms, driven by both genetic and environmental factors. Recent studies of SD in fishes have shown that not all taxa fit the classic paradigm of sex chromosome evolution and diverse SD methods can be found even among closely related species. Here, we apply a suite of genomic approaches to investigate sex-biased genomic variation in eight species of Sebastes rockfish found in the northeast Pacific Ocean. Using recently assembled chromosome-level rockfish genomes, we leverage published sequence data to identify disparate sex chromosomes and sex-biased loci in five species. We identify two putative male sex chromosomes in S. diaconus, a single putative sex chromosome in the sibling species S. carnatus and S. chrysomelas, and an unplaced sex determining contig in the sibling species S. miniatus and S. crocotulus. Our study provides evidence for disparate means of sex determination within a recently diverged set of species and sheds light on the diverse origins of sex determination mechanisms present in the animal kingdom.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.