Sox5 is involved in germ-cell regulation and sex determination in medaka following co-option of nested transposable elements

Dissecting SOX2 expression and function reveals an association with multiple signaling pathways during embryonic development and in cancer progression

SRY (Sex Determining Region) box 2 (SOX2) is an essential transcription factor that plays crucial roles in activating genes involved in pre- and post-embryonic development, adult tissue homeostasis, and lineage specifications. SOX2 maintains the self-renewal property of stem cells and is involved in the generation of induced pluripotency stem cells. SOX2 protein contains a particular high-mobility group domain that enables SOX2 to achieve the capacity to participate in a broad variety of functions. The information about the involvement of SOX2 with gene regulatory elements, signaling networks, and microRNA is gradually emerging, and the higher expression of SOX2 is functionally relevant to various cancer types. SOX2 facilitates the oncogenic phenotype via cellular proliferation and enhancement of invasive tumor properties. Evidence are accumulating in favor of three dimensional (higher order) folding of chromatin and epigenetic control of the SOX2 gene by chromatin modifications, which implies that the expression level of SOX2 can be modulated by epigenetic regulatory mechanisms, specifically, via DNA methylation and histone H3 modification. In view of this, and to focus further insights into the roles SOX2 plays in physiological functions, involvement of SOX2 during development, precisely, the advances of our knowledge in pre- and post-embryonic development, and interactions of SOX2 in this scenario with various signaling pathways in tumor development and cancer progression, its potential as a therapeutic target against many cancers are summarized and discussed in this article.

Genome-wide identification and expression profiles of sex-related gene families in the Pacific abalone Haliotis discus hannai

In recent years, members of the Dmrt family, TGF-β superfamily and Sox family have been recognized as crucial genes for sex determination/differentiation across diverse animal species. Nevertheless, knowledge regarding the abundance and potential functions of these genes in abalone remains limited. In this study, a total of 5, 10, and 7 members of the Dmrt family, the TGF-β superfamily and the Sox family, respectively, were identified in the Pacific abalone Haliotis discus hannai. Sequence characteristics, phylogenetic relationships and spatiotemporal expression profiles of these genes were investigated. Notably, HdDmrt-04 (Dmrt1/1L-like) emerged as a potential mollusc-specific gene with a preponderance for expression in the testis. Interestingly, none of the TGF-β superfamily members exhibited specific or elevated expression in the gonads, highlighting the need for further investigation into their role in abalone sex differentiation. The Sox proteins in H. discus hannai were categorized into 7 subfamilies: B1, B2, C, D, E, F, and H. Among them, HdSox-07 (SoxH-like) was observed to play a crucial role in testis development, while HdSox-03 (SoxB1-like) and HdSox-04 (SoxC-like) probably cooperate in abalone ovary development. Taken together, the results of the present study suggested that HdDmrt-04 and HdSox-07 can be used as male-specific markers for gonad differentiation in H. discus hannai and imply conservation of their functions across invertebrates and vertebrates. Our findings provide new insights into the evolution and genetic structure of the Dmrt family, the TGF-β superfamily and the Sox family in abalone and pave the way for a deeper understanding of sex differentiation in gastropods.

De Novo Assembly, Characterization, and Comparative Transcriptome Analysis of Mature Male and Female Gonads of Rabbitfish (Siganus oramin) (Bloch & Schneider, 1801)

The rabbitfish, Siganus oramin, is a commercially important table fish in southeastern China. However, there have been few studies on its gonad development and reproduction regulation. Comparative transcriptome analysis was first performed on adult male and female gonads of S. oramin. In total, 47,070 unigenes were successfully assembled and 22,737 unigenes were successfully annotated. Through comparative transcriptome analysis of male and female gonads, a total of 6722 differentially expressed genes were successfully identified, with 3528 upregulated genes and 3154 downregulated genes in the testes. In addition, 39 differentially expressed reproduction-related genes were identified. Finally, quantitative real-time PCR was used to validate the expression levels of several differentially expressed genes. These results provide important data for further studying the function of reproduction-related genes and the molecular mechanism regulating gonad development and reproduction in S. oramin.

Molecular functional characterization of the setdb1 and its potential target gene sox5 illuminate the histone modification-mediated orchestration of gonadal development in Chinese tongue sole (Cynoglossus semilaevis)

SET (SuVar3-9, Enhancer of Zeste, Trithorax) domain bifurcated histone lysine methyltransferase 1, setdb1, is the predominant histone lysine methyltransferase catalyzing H3K9me3. Prior studies have illustrated that setdb1 and H3K9me3 critically regulate sex differentiation and gametogenesis. However, the molecular details by which setdb1 is involved in these processes in fish have been poorly reported. Here, we cloned and characterized the setdb1 ORF (open reading frame) sequence from Chinese tongue sole (Cynoglossus semilaevis). The setdb1 ORF sequence was 3,669 bp, encoding a 1,222-amino-acid protein. Phylogenetic analysis showed that setdb1 was structurally conserved. qRT-PCR revealed that setdb1 had a high expression level in the testes at 12 mpf (months post fertilization). Single-cell RNA-seq data at 24 mpf indicated that setdb1 was generally expressed in spermatogenic cells at each stage except for sperm and was centrally expressed in oogonia. H3K9me3 modification was observed in gonads with the immunofluorescence technique. Furthermore, the overexpression experiment suggested that sox5 was a candidate target of setdb1. sox5 was abundantly expressed in male and pseudomale gonads at 24 mpf. Single-cell RNA-seq data showed that sox5 was mainly expressed in spermatogonia and its expression gradually declined with differentiation. Taken together, our findings imply that setdb1 regulates sox5 transcription in gonads, which provides molecular clues into histone modification-mediated orchestration of sex differentiation and gametogenesis.

Transposable elements acquire time- and sex-specific transcriptional and epigenetic signatures along mouse fetal gonad development

Gonadal sex determination in mice is a complex and dynamic process, which is crucial for the development of functional reproductive organs. The expression of genes involved in this process is regulated by a variety of genetic and epigenetic mechanisms. Recently, there has been increasing evidence that transposable elements (TEs), which are a class of mobile genetic elements, play a significant role in regulating gene expression during embryogenesis and organ development. In this study, we aimed to investigate the involvement of TEs in the regulation of gene expression during mouse embryonic gonadal development. Through bioinformatics analysis, we aimed to identify and characterize specific TEs that operate as regulatory elements for sex-specific genes, as well as their potential mechanisms of regulation. We identified TE loci expressed in a time- and sex-specific manner along fetal gonad development that correlate positively and negatively with nearby gene expression, suggesting that their expression is integrated to the gonadal regulatory network. Moreover, chromatin accessibility and histone post-transcriptional modification analyses in differentiating supporting cells revealed that TEs are acquiring a sex-specific signature for promoter-, enhancer-, and silencer-like elements, with some of them being proximal to critical sex-determining genes. Altogether, our study introduces TEs as the new potential players in the gene regulatory network that controls gonadal development in mammals.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

Regulatory mechanisms of SoxD transcription factors and their influences on male fertility

Members of the SRY-related box (SOX) subfamily D (SoxD) of transcription factors are well conserved among vertebrate species and play important roles in different stages of male reproductive development. In mammals, the SoxD subfamily contains three members: SOX5, SOX6 and SOX13. Here, we describe their implications in testicular development and spermatogenesis, contributing to fertility. We also cover the mechanisms of action of SoxD transcription factors in gene regulation throughout male development. The specificity of activation of target genes by SoxD members depends, in part, on their post-translational modifications and interactions with other partners. Sperm production in adult males requires the coordination in the regulation of gene expression by different members of the SoxD subfamily of transcription factors in the testis. Specifically, the regulation of genes promoting adequate spermatogenesis by SoxD members is discussed in comparison between species.

Evaluating the role of sexual antagonism in the evolution of sex chromosomes: new data from fish

The recent increase in available molecular and genomic data for diverse taxa helps to shed new light on long-standing theories. Research into sex chromosome evolution has particularly benefited from a growing number of studies of fish, motivated by their highly diverse mechanisms of sex determination. Sexual antagonism is regularly cited as an influential force in sex chromosome emergence; however, this so far proves difficult to demonstrate. In this review, we highlight recent developments in the investigation of sexual antagonism in sex chromosome research in fish. We find strong emphasis placed on study-organism specific genomic features and patterns of recombination, rather than evidence for a comprehensive role of sexual antagonism. In this light, we discuss the alternative models of sex chromosome evolution. We conclude that fish represents a key resource for further research, provided attention is given to species-specific effects while simultaneously integrating comparative studies across taxa for a vital and comprehensive understanding of sex chromosome evolution and investigation of proposed models.

Sex steroid profile during oocyte development and maturation in the intertidal worm Marphysa madrasi (Polychaeta: Eunicidae) from the east coast of India

Marphysa madrasi is a commercially valuable maturation diet in crustacean aquaculture. This study presents the first detailed investigation of oogenesis in the intertidal polychaete worm M. madrasi and reports the steroid profile during oocyte growth and development. Oogenesis is extraovarian type I, originating from coelomic epithelial cells, with four stages of development - primary growth, early vitellogenic, late vitellogenic, and maturation. The primary growth phase contains oogonial cells and previtellogenic oocyte clusters in the early, mid, and late stages of development form a dispersed ovary attached to blood vessels. The late previtellogenic oocytes detach from the ovary at the onset of vitellogenesis. The detached oocytes complete vitellogenesis and final maturation in the coelomic fluid as solitary free-floating cells without any connection with follicle cells. The worms display asynchronous reproduction with a heterogeneous population of developing oocytes. Steroid extracts from the polychaete homogenates in different stages of oogenesis were identified by HPLC and confirmed by LC-MS/MS. In M. madrasi, two vertebrate-type steroids, pregnenolone (P5) and 17α-hydroxyprogesterone (17-OHP) were detected and quantified. The P5 levels were low in immature worms but increased significantly by ∼ 8.3-fold in the previtellogenic stage and peaked during oocyte maturation. 17-OHP levels were low in immature worms but gradually increase as the oogenesis progress to the primary growth and early vitellogenic phase, with a significant increase (p < 0.001) during the late vitellogenic phase. Although an increase in the concentration of P5 and 17-OHP during vitellogenesis and maturation of oocytes points to a possible role in reproduction, the absence of other vertebrate-type steroids in the investigated polychaete signifies a plausible uptake of P5 and 17-OHP from the environment.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Mature Gonads in Spinibarbus hollandi

Spinibarbus hollandi is an important commercial aquaculture species in southeastern China, but with long maturity period and low egg laying amount. However, there has been little study of its gonad development and reproductive regulation, which limits aquaculture production. Here, for the first time, gonadal transcriptomes of male and female S. hollandi were analyzed. A total of 167,152 unigenes were assembled, with only 48,275 annotated successfully. After comparison, a total of 21,903 differentially expressed genes were identified between male and female gonads, of which 16,395 were upregulated and 5508 were downregulated in the testis. In addition, a large number of differentially expressed genes participating in reproduction, gonad formation and differentiation, and gametogenesis were screened out and the differential expression profiles of partial genes were further validated using quantitative real-time PCR. These results will provide basic information for further research on gonad differentiation and development in S. hollandi.

Identification of Iguania Ancestral Syntenic Blocks and Putative Sex Chromosomes in the Veiled Chameleon (Chamaeleo calyptratus, Chamaeleonidae, Iguania)

The veiled chameleon (Chamaeleo calyptratus) is a typical member of the family Chamaeleonidae and a promising object for comparative cytogenetics and genomics. The karyotype of C. calyptratus differs from the putative ancestral chameleon karyotype (2n = 36) due to a smaller chromosome number (2n = 24) resulting from multiple chromosome fusions. The homomorphic sex chromosomes of an XX/XY system were described recently using male-specific RADseq markers. However, the chromosomal pair carrying these markers was not identified. Here we obtained chromosome-specific DNA libraries of C. calyptratus by chromosome flow sorting that were assigned by FISH and sequenced. Sequence comparison with three squamate reptiles reference genomes revealed the ancestral syntenic regions in the C. calyptratus chromosomes. We demonstrated that reducing the chromosome number in the C. calyptratus karyotype occurred through two fusions between microchromosomes and four fusions between micro-and macrochromosomes. PCR-assisted mapping of a previously described Y-specific marker indicates that chromosome 5 may be the sex chromosome pair. One of the chromosome 5 conserved synteny blocks shares homology with the ancestral pleurodont X chromosome, assuming parallelism in the evolution of sex chromosomes from two basal Iguania clades (pleurodonts and acrodonts). The comparative chromosome map produced here can serve as the foundation for future genome assembly of chameleons and vertebrate-wide comparative genomic studies.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Analysis of yellow mutant rainbow trout transcriptomes at different developmental stages reveals dynamic regulation of skin pigmentation genes

Yellow mutant rainbow trout (YR), an economically important aquaculture species, is popular among consumers due to its excellent meat quality and attractive appearance. Skin color is a key economic trait for YR, but little is known about the molecular mechanism of skin color development. In this study, YR skin transcriptomes were analyzed to explore temporal expression patterns of pigmentation-related genes in three different stages of skin color development. In total, 16,590, 16,682, and 5619 genes were differentially expressed between fish at 1 day post-hatching (YR1d) and YR45d, YR1d and YR90d, and YR45d and YR90d. Numerous differentially expressed genes (DEGs) associated with pigmentation were identified, and almost all of them involved in pteridine and carotenoid synthesis were significantly upregulated in YR45d and YR90d compared to YR1d, including GCH1, PTS, QDPR, CSFIR1, SLC2A11, SCARB1, DGAT2, PNPLA2, APOD, and BCO2. Interestingly, many DEGs enriched in melanin synthesis pathways were also significantly upregulated, including melanogenesis (MITF, MC1R, SLC45A2, OCA2, and GPR143), tyrosine metabolism (TYR, TYRP1, and DCT), and MAPK signaling (KITA) pathways. Using short time-series expression miner, we identified eight differential gene expression pattern profiles, and DEGs in profile 7 were associated with skin pigmentation. Protein–protein interaction network analysis showed that two modules were related to xanthophores and melanophores. In addition, 1,812,329 simple sequence repeats and 2,011,334 single-nucleotide polymorphisms were discovered. The results enhance our understanding of the molecular mechanism underlying skin pigmentation in YR, and could accelerate the molecular breeding of fish species with valuable skin color traits and will likely be highly informative for developing new therapeutic approaches to treat pigmentation disorders and melanoma.

Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish

Background

Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens.

Results

We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication.

Conclusions

Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01205-y.

The Important Role of Sex-Related Sox Family Genes in the Sex Reversal of the Chinese Soft-Shelled Turtle (Pelodiscus sinensis)

The Chinese soft-shelled turtle Pelodiscus sinensis shows obvious sexual dimorphism. The economic and nutrition value of male individuals are significantly higher than those of female individuals. Pseudo-females which are base to all-male breeding have been obtained by estrogen induction, while the gene function and molecular mechanism of sex reversal remain unclear in P. sinensis. Here, comparative transcriptome analyses of female, male, and pseudo-female gonads were performed, and 14,430 genes differentially expressed were identified in the pairwise comparison of three groups. GO and KEGG analyses were performed on the differentially expressed genes (DEGs), which mainly concentrated on steroid hormone synthesis. Furthermore, the results of gonadal transcriptome analysis revealed that 10 sex-related sox genes were differentially expressed in males vs. female, male vs. pseudo-female, and female vs. pseudo-female. Through the differential expression analysis of these 10 sox genes in mature gonads, six sox genes related to sex reversal were further screened. The molecular mechanism of the six sox genes in the embryo were analyzed during sex reversal after E₂ treatment. In mature gonads, some sox family genes, such as sox9sox12, and sox30 were highly expressed in the testis, while sox1, sox3, sox6, sox11, and sox17 were lowly expressed. In the male embryos, exogenous estrogen can activate the expression of sox3 and inhibit the expression of sox8, sox9, and sox11. In summary, sox3 may have a role in the process of sex reversal from male to pseudo-female, when sox8 and sox9 are inhibited. Sox family genes affect both female and male pathways in the process of sex reversal, which provides a new insight for the all-male breeding of the Chinese soft-shelled turtle.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

A new linkage group involved in sex determination in haplotilapiine cichlids

The high diversity of sex chromosomes and sex determination systems among haplotilapiines suggests that this large cichlid clade is a good model for investigating the evolution of genetics of sex determination. Nonetheless, information about sex determination in this clade remains sparse. The present study reports a microsatellite marker that is closely associated with sex in Xenotilapia rotundiventralis from Lake Tanganyika, Africa. This study is the first to suggest the role of linkage group 17 in sex determination in haplotilapiine cichlids.

Comparative transcriptomic analysis reveals the gonadal development-related gene response to environmental temperature in Mauremys mutica

The Asian yellow pond turtle (Mauremys mutica) displays temperature-dependent sex determination (TSD), in which incubation temperature during embryonic development determines the sexual fate of the individual. However, the mechanism of the sex determination/differentiation of Mauremys mutica remains a mystery. Here, we first analyzed the temperature-specific gonadal transcriptomes of Mauremys mutica prior to gonad formation and gonads during the thermosensitive period. We uncovered a list of candidates that respond to temperature stimuli enriched in several categories, such as heat shock protein family members dnajb6a, dnaja4, hspa8 and hsp90aa1, temperature sensor genes mmp17 and mmp28, and putative novel temperature-responsive genes tmco6, gria3 and eif3f. Notably, striking differences were identified in the expression profiles of genes underlying sexual development, such as tex15, insr, igf1r, cirbp, esr1, dmrt2 and Serpinh1. Moreover, we analyzed the similarity and divergence of the timecourse of gene expression among Mauremys mutica and two other reported TSD turtles (Trachemys scripta and Chrysemys picta). The shared genes revealed the common gonad-specific regulatory mechanisms existing in these three TSD turtles that initiate their sexual development. Therefore, our findings could provide basic data to elucidate the mechanisms of sex determination/differentiation of M. mutica, even contributing to further understanding of these mechanisms in other TSD turtles.

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)'.

How did the guppy Y chromosome evolve?

The sex chromosome pairs of many species do not undergo genetic recombination, unlike the autosomes. It has been proposed that the suppressed recombination results from natural selection favouring close linkage between sex-determining genes and mutations on this chromosome with advantages in one sex, but disadvantages in the other (these are called sexually antagonistic mutations). No example of such selection leading to suppressed recombination has been described, but populations of the guppy display sexually antagonistic mutations (affecting male coloration), and would be expected to evolve suppressed recombination. In extant close relatives of the guppy, the Y chromosomes have suppressed recombination, and have lost all the genes present on the X (this is called genetic degeneration). However, the guppy Y occasionally recombines with its X, despite carrying sexually antagonistic mutations. We describe evidence that a new Y evolved recently in the guppy, from an X chromosome like that in these relatives, replacing the old, degenerated Y, and explaining why the guppy pair still recombine. The male coloration factors probably arose after the new Y evolved, and have already evolved expression that is confined to males, a different way to avoid the conflict between the sexes.

We report new findings concerning the long-studied the guppy XY pair, which has remained somewhat mystifying. We show that it can be understood as a case of a recent sex chromosome turnover event in which an older, highly degenerated Y chromosome was lost, and creation of a new sex chromosome from the ancestral X. This chromosome acquired a male-determining factor, possibly by a mutation in (or a duplication of) a previously X-linked gene, or (less likely) by movement of an ancestral Y-linked maleness factor onto the X. We relate the findings to theoretical models of such events, and argue that the proposed change was free from factors thought to impede such turnovers. The change resulted in the intriguing situation where the X chromosome is old and the Y is much younger, and we discuss some other species where a similar change seems likely to have occurred.

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.

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)'.

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)'.

Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY. By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

Nuclear heme oxygenase-1 improved the hypoxia-mediated dysfunction of blood-spinal cord barrier via the miR-181c-5p/SOX5 signaling pathway

Our previous study demonstrated that adenovirus-delivered GFP nuclear heme oxygenase-1 (nuclear HO-1, NHO-1) fragments lacking 23 amino acids at the C-terminus (Ad-GFP-HO-1C[INCREMENT]23) showed the potential therapeutic effects mediated by its improvement of the blood-spinal cord barrier (BSCB) integrity. However, the NHO-1-mediated molecular mechanism in regulating the BSCB function remains unclear. The BSCB model in vitro was established via a coculture of primary rat brain microvascular endothelial cells (RBMECs) and spinal cord astrocytes on transwell system. NHO-1 markedly reduced the disruption of the BSCB integrity induced by hypoxia. And NHO-1 significantly attenuated the expression of miR-181c-5p, but increased the expression level of SOX5 protein. miR-181c-5p was shown as an essential miRNA for increasing the BSCB permeability under hypoxia condition. Furthermore, we identified that miR-181c-5p could regulate the expression of SOX5 through binding to the 3'-UTR of its mRNA. And the decreased BSCB permeability and upregulation of tight junction (TJ) protein expression induced by NHO-1 could be partly reversed by the inhibition of SOX5 or miR-181c-5p (+). The present study results provide a better understanding of the molecular mechanisms induced by NHO-1 in improving the BSCB integrity, which is associated with the regulation of miR-181c-5p/SOX5/TJ signaling pathway.

Role of sox30 in regulating testicular steroidogenesis of common carp

Expression of transcription factors is crucial for the regulation of steroidogenesis and gonadal development in fish. SRY-related box (SOX) proteins regulate gene expression of various events related to vertebrate reproduction. This study reports the role of sox30 and its influence on sox9a/b in regulating testicular steroidogenesis of the common carp, Cyprinus carpio. Tissue distribution showed predominant expression of sox30 in gonads, while gonadal ontogeny indicated significant dimorphic expression of sox30 from 120 days post hatch. Higher sox30 transcripts during the spawning season, an elevation of sox30 after human chorionic gonadotropin induction, and 11-ketotestosterone (11-KT) treatment authenticate gonadotropin dependency. Treatment of 17α-methyl-di-hydroxy-testosterone to juvenile common carp for mono-sex induction, vis-à-vis elevated sox30 expression. Sox30 protein was detected abundantly in spermatocytes and spermatid/sperm of carp testis. Transient silencing of sox30 using small interfering RNAs decreased sox9a/b expression, lead to downregulation of certain molecule/factor, transcription factor, germ/stem cell marker, and steroidogenesis-related enzyme genes. Serum testosterone and 11-KT decreased significantly upon transient silencing of sox30, in vivo. Concomitantly, a reduction in testicular microsomal 11-β hydroxysteroid dehydrogenase activity was observed. These results demonstrate the influence of sox30 as well as sox9a/b in the regulation of testicular steroidogenesis in common carp.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Transcriptional Inhibition of Sp-IAG by Crustacean Female Sex Hormone in the Mud Crab, Scylla paramamosain

In crustaceans, the regulation of sex differentiation is mediated by insulin-like androgenic hormone (IAG) and crustacean female sex hormone (CFSH). CFSH is reported to inhibit IAG gene (Sp-IAG) expression in the mud crab Scylla paramamosain, but the regulatory mechanism is not well understood. A 2674 bp 5′ flanking Sp-IAG contains many potential transcription factor binding sites. In this study, analysis of serially deleted 5′ flanking Sp-IAG and site-directed mutation (SDM) of transcription factor binding sites of the same gene showed that the promoter activity of reporter vectors with Sox-5-binding site, signal transducers and activators of transcription (STAT)-binding site and activator protein 1 (AP-1)-binding site were significantly higher than that of vectors without these regions, suggesting that they were involved in transcriptional regulation of Sp-IAG expression. The expression analysis of these transcription factor showed that there was no difference in the level of mRNA in Sox-5 and AP-1 in androgenic gland treated with recombinant CFSH, but expression of Sp-STAT was significantly reduced, suggesting that CFSH regulates the expression of Sp-STAT, inhibiting its function to regulate Sp-IAG. Further experiment revealed that RNAi mediated Sp-STAT gene knockdown reduced the expression of Sp-IAG. These results suggested that Sp-CFSH regulates Sp-IAG by inhibiting STAT. This is a pioneering finding on the transcriptional mechanism of IAG gene in crustaceans.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Genomic Access to the Diversity of Fishes

The number of fishes exceeds that of all other vertebrates both in terms of species numbers and in their morphological and phylogenetic diversity. They are an ecologically and economically important group and play an essential role as a resource for humans. This makes the genomic exploration of fishes an important area of research, both from an applied and a basic research perspective. Fish genomes can vary greatly in complexity, which is partially due to differences in size and content of repetitive DNA, a history of genome duplication events and because fishes may be polyploid, all of which complicate the assembly and analysis of genome sequences. However, the advent of modern sequencing techniques now facilitates access to genomic data that permit genome-wide exploration of genetic information even for previously unexplored species. The development of genomic resources for fishes is spearheaded by model organisms that have been subject to genetic analysis and genome sequencing projects for a long time. These offer a great potential for the exploration of new species through the transfer of genomic information in comparative analyses. A growing number of genome sequencing projects and the increasing availability of tools to assemble and access genomic information now move boundaries between model and nonmodel species and promises progress in many interesting but unexplored species that remain to be studied.

Sex and the TEs: transposable elements in sexual development and function in animals

Transposable elements are endogenous DNA sequences able to integrate into and multiply within genomes. They constitute a major source of genetic innovations, as they can not only rearrange genomes but also spread ready-to-use regulatory sequences able to modify host gene expression, and even can give birth to new host genes. As their evolutionary success depends on their vertical transmission, transposable elements are intrinsically linked to reproduction. In organisms with sexual reproduction, this implies that transposable elements have to manifest their transpositional activity in germ cells or their progenitors. The control of sexual development and function can be very versatile, and several studies have demonstrated the implication of transposable elements in the evolution of sex. In this review, we report the functional and evolutionary relationships between transposable elements and sexual reproduction in animals. In particular, we highlight how transposable elements can influence expression of sexual development genes, and how, reciprocally, they are tightly controlled in gonads. We also review how transposable elements contribute to the organization, expression and evolution of sexual development genes and sex chromosomes. This underscores the intricate co-evolution between host functions and transposable elements, which regularly shift from a parasitic to a domesticated status useful to the host.

Evolution and diversity of transposable elements in fish genomes

Transposable elements (TEs) are genomic sequences that can move, multiply, and often form sizable fractions of vertebrate genomes. Fish belong to a unique group of vertebrates, since their karyotypes and genome sizes are more diverse and complex, with probably higher diversity and evolution specificity of TE. To investigate the characteristics of fish TEs, we compared the mobilomes of 39 species, and observed significant variation of TE content in fish (from 5% in pufferfish to 56% in zebrafish), along with a positive correlation between fish genome size and TE content. In different classification hierarchies, retrotransposons (class), long terminal repeat (order), as well as Helitron, Maverick, Kolobok, CMC, DIRS, P, I, L1, L2, and 5S (superfamily) were all positively correlated with fish genome size. Consistent with previous studies, our data suggested fish genomes to not always be dominated by DNA transposons; long interspersed nuclear elements are also prominent in many species. This study suggests CR1 distribution in fish genomes to be obviously regular, and provides new clues concerning important events in vertebrate evolution. Altogether, our results highlight the importance of TEs in the structure and evolution of fish genomes and suggest fish species diversity to parallel transposon content diversification.

A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka

Dmrt1 is a highly conserved transcription factor, which is critically involved in regulation of gonad development of vertebrates. In medaka, a duplicate of dmrt1-acting as master sex-determining gene-has a tightly timely and spatially controlled gonadal expression pattern. In addition to transcriptional regulation, a sequence motif in the 3' UTR (D3U-box) mediates transcript stability of dmrt1 mRNAs from medaka and other vertebrates. We show here that in medaka, two RNA-binding proteins with antagonizing properties target this D3U-box, promoting either RNA stabilization in germ cells or degradation in the soma. The D3U-box is also conserved in other germ-cell transcripts, making them responsive to the same RNA binding proteins. The evolutionary conservation of the D3U-box motif within dmrt1 genes of metazoans-together with preserved expression patterns of the targeting RNA binding proteins in subsets of germ cells-suggest that this new mechanism for controlling RNA stability is not restricted to fishes but might also apply to other vertebrates.

Differential isoform expression and alternative splicing in sex determination in mice

Background

Alternative splicing (AS) may play an important role in gonadal sex determination (GSD) in mammals. The present study was designed to identify differentially expressed isoforms and AS modifications accompanying GSD in mice.

Results

Using deep RNA-sequencing, we performed a transcriptional analysis of XX and XY gonads during sex determination on embryonic days 11 (E11) and 12 (E12). Analysis of differentially expressed genes (DEG) identified hundreds of genes related to GSD and early sex differentiation that may represent good candidates for sex reversal. Expression at time point E11 in males was significantly enriched in RNA splicing and mRNA processing Gene Ontology terms. Differentially expressed isoform analysis identified hundreds of specific isoforms related to GSD, many of which showed no differences in the DEG analysis. Hundreds of AS events were identified as modified at E11 and E12. Female E11 gonads featured sex-biased upregulation of intron retention (in genes related to regulation of transcription, protein phosphorylation, protein transport and mRNA splicing) and exon skipping (in genes related to chromatin repression) suggesting AS as a post-transcription mechanism that controls sex determination of the bipotential fetal gonad.

Conclusion

Our data suggests an important role of splicing regulatory mechanisms for sex determination in mice.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5572-x) contains supplementary material, which is available to authorized users.

Sex biased expression of anti-Mullerian hormone (amh) gene in a live bearing fish, Gambusia holbrooki: Evolutionary implications and potential role in sex differentiation

The amh, a member of transforming growth factor-β (TGF-β) family, is known to play a critical role in vertebrate male sex differentiation, with its paralogue/s evolving to determine sex in few heterogametic (XX/XY) teleosts. However, it remains relatively unexplored in the reproductively unique live bearing teleosts. Therefore, this study comparatively examined the structure and content of G. holbrooki amh as well as characterised its expression. A paralogous Y-specific amh (amhy) was not detected, suggesting an unlikely role in sex determination. Two transcripts (1.4 and 1.5 kb) were detected in adults: the larger (1.5 kb) retaining intron 5, coding for a truncated AMH-N and no TGF-β domain. The small (1.4 kb) transcript, had both domains intact and clustered with members of Poeciliidae. In contrast to other vertebrates, a higher conservation between the N- rather than the C- terminus of amh in Poeciliidae was observed, suggesting an adaptation that may be unique to live bearing teleosts. The amh expression was 6 times higher in brain of both sexes and testis compared with ovaries (p = .001). Intriguingly, female splenic tissues showed 10 times higher expression (p = .006) and such female bias splenic expression has not been reported in any teleosts. Ontogenic expression was 25 times higher in male embryos at gastrulation stage (p = .001), much earlier than those reported in egg-laying teleosts. Such heightened expression in male embryos suggests a repressive role associated with proliferation and migration of primordial germ cells (PGCs) that are known to occur earlier at blastulation in teleosts-potentially influencing gonadal fate.

Diverse and variable sex determination mechanisms in vertebrates

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

A Comparative View on Sex Differentiation and Gametogenesis Genes in Lungfish and Coelacanths

Gonadal sex differentiation and reproduction are the keys to the perpetuation of favorable gene combinations and positively selected traits. In vertebrates, several gonad development features that differentiate tetrapods and fishes are likely to be, at least in part, related to the water-to-land transition. The collection of information from basal sarcopterygians, coelacanths, and lungfishes, is crucial to improve our understanding of the molecular evolution of pathways involved in reproductive functions, since these organisms are generally regarded as “living fossils” and as the direct ancestors of tetrapods. Here, we report for the first time the characterization of >50 genes related to sex differentiation and gametogenesis in Latimeria menadoensis and Protopterus annectens. Although the expression profiles of most genes is consistent with the intermediate position of basal sarcopterygians between actinopterygian fish and tetrapods, their phylogenetic placement and presence/absence patterns often reveal a closer affinity to the tetrapod orthologs. On the other hand, particular genes, for example, the male gonad factor gsdf (Gonadal Soma-Derived Factor), provide examples of ancestral traits shared with actinopterygians, which disappeared in the tetrapod lineage.