Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

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.

Chromosomal rearrangements and the first indication of an ♀X1 X1 X2 X2 /♂X1 X2 Y sex chromosome system in Rineloricaria fishes (Teleostei: Siluriformes)

Rineloricaria is the most diverse genus within the freshwater fish subfamily Loricariinae, and it is widely distributed in the Neotropical region. Despite limited cytogenetic data, records from southern and south-eastern Brazil suggest a high rate of chromosomal rearrangements in this genus, mirrored in remarkable inter- and intraspecific karyotype variability. In the present work, we investigated the karyotype features of Rineloricaria teffeana, an endemic representative from northern Brazil, using both conventional and molecular cytogenetic techniques. We revealed different diploid chromosome numbers (2n) between sexes (33♂/34♀), which suggests the presence of an ♀X₁ X₁ X₂ X₂ /♂X₁ X₂ Y multiple sex chromosome system. The male-limited Y chromosome was the largest and the only biarmed element in the karyotype, implying Y-autosome fusion as the most probable mechanism behind its origination. C-banding revealed low amounts of constitutive heterochromatin, mostly confined to the (peri)centromeric regions of most chromosomes (including the X₂ and the Y) but also occupying the distal regions of a few chromosomal pairs. The chromosomal localization of the 18S ribosomal DNA (rDNA) clusters revealed a single site on chromosome pair 4, which was adjacent to the 5S rDNA cluster. Additional 5S rDNA loci were present on the autosome pair 8, X₁ chromosome, and in the presumed fusion point on the Y chromosome. The probe for telomeric repeat motif (TTAGGG)_n revealed signals of variable intensities at the ends of all chromosomes except for the Y chromosome, where no detectable signals were evidenced. Male-to-female comparative genomic hybridization revealed no sex-specific or sex-biased repetitive DNA accumulations, suggesting a presumably low level of neo-Y chromosome differentiation. We provide evidence that rDNA sites might have played a role in the formation of this putative multiple sex chromosome system and that chromosome fusions originate through different mechanisms among different Rineloricaria species.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.