Sex dimorphic expression of five dmrt genes identified in the Atlantic cod genome. The fish-specific dmrt2b diverged from dmrt2a before the fish whole-genome duplication

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.

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

Genome-wide identification, phylogeny and expressional profile of the Dmrt gene family in Chinese sturgeon (Acipenser sinensis)

Chinese sturgeon Dmrt gene family was identified and characterized for the first time. A total of 5 putative Dmrt genes were identified. The gene structure, conserved protein domain and the phylogenetic relationship of Dmrt gene family were systematically analyzed. The expressed profile of Chinese sturgeon Dmrt genes in gonad, pituitary and hypothalamus in the male and female were investigated. The results indicated that the accumulation of Dmrt genes was involved in different tissues, and the expression profile also differed among each Dmrt genes. ASDmrt1A, ASDmrt2, ASDmrt3, and ASDmrtA1 were highly expressed in the testis in comparison with other tissue. This result showed that ASDmrt1A, ASDmrt2, ASDmrt3, and ASDmrtA1 played an important role in the development of testicle, and may be useful tool in distinguishing between male and female of Chinese sturgeon. Our study will provide a basis for additional analyses of Chinese sturgeon Dmrt genes. This systematic analysis provided a foundation for further functional characterization of Dmrt genes with an aim of study of Chinese sturgeon Dmrt gene family.

Genome-wide investigation of the DMRT gene family sheds new insight into the regulation of sex differentiation in spotted knifejaw (Oplegnathus punctatus) with fusion chromosomes (Y)

The role of the DMRT family in male sex determination and differentiation is significant, but its regulatory role in spotted knifejaw with Y fusion chromosomes remains unclear. Through genome-wide scanning, transcriptome analysis, qPCR, FISH, and RNA interference (RNAi), we investigated the DMRT family and the dmrt1-based sex regulation network. Seven DMRTs were identified (DMRT1/2 (2a,2b)/6, DMRT4/5, DMRT3), and dmrt gene dispersion among chromosomes is possibly driven by three whole-genome duplications. Transcriptome analysis enriched genes were associated with sex regulation and constructed a network associated with dmrt1. qPCR and FISH results showed the expression dimorphism of sex-related genes in dmrt-related regulatory networks. RNAi experiments indicated a distinct sex regulation mode in spotted knifejaw. Dmrt1 knockdown upregulated male-related genes (sox9a, sox9b, dmrt1, amh, amhr2) and hsd11b2 expression, which is critical for androgen synthesis. Amhr2 is located on the heterozygous chromosome (Y) and is specifically localized in primary spermatocytes, and is extremely upregulated after dmrt1 knockdown which suggested besides the important role of dmrt1 in male differentiation, the amhr2 along with amhr2/amh system, also play important regulatory roles in maintaining high expression of the hsd11b2 and male differentiation. This study aims to further investigate sex regulatory mechanisms in species with fusion chromosomes.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Full-length gonad transcriptome analysis of Amur sturgeon Dmrt family genes: identification, characterization, and expression patterns during gonadal differentiation

The regulatory mechanisms that govern sex differentiation in sturgeon are still poorly understood. The doublesex and Mab-3-related transcription factor (Dmrt) gene family is known for its extensive roles in sex determination and differentiation across vertebrates. This study aimed to identify new members of sturgeon Dmrt family genes and core actors in the gonadal differentiation of Amur sturgeon. A full-length gonad transcriptome database was exploited to identify Dmrt gene orthologs. Analyses of phylogenetic relationships and selection pressure were performed, and tissue expression profiles and spatiotemporal expression patterns in gonads were then analyzed using real-time PCR. In total, five Dmrt family genes were identified from the full-length gonad transcriptome, including Dmrt2, DmrtA1, DmrtA2, DmrtB1a, and DmrtB1b. Phylogenetic analysis showed that these genes were clustered into clades corresponding to the doublesex/Mav-3 (DM) genes of vertebrates. Furthermore, the analysis of evolutionary selective pressure indicated that DmrtB1a and DmrtB1b were subject to positive selection, suggesting the existence of adaptive evolution in sturgeon. The extensive tissue expression profiling of each Dmrt family gene revealed typical characteristics. Remarkably, according to a spatiotemporal expression pattern analysis, in later stages, DmrtB1b expression increased in testes and was significantly higher in testes than in ovaries at 24 months after hatching (M) and 36 M. This study provides a genetic resource of full-length Dmrt family genes and increases the understanding of Dmrt functions in sex differentiation in sturgeon.

DMRT Transcription Factors in the Control of Nervous System Sexual Differentiation

Sexual phenotypic differences in the nervous system are one of the most prevalent features across the animal kingdom. The molecular mechanisms responsible for sexual dimorphism throughout metazoan nervous systems are extremely diverse, ranging from intrinsic cell autonomous mechanisms to gonad-dependent endocrine control of sexual traits, or even extrinsic environmental cues. In recent years, the DMRT ancient family of transcription factors has emerged as being central in the development of sex-specific differentiation in all animals in which they have been studied. In this review, we provide an overview of the function of Dmrt genes in nervous system sexual regulation from an evolutionary perspective.

Genome-Wide Identification, Phylogeny, and Expression Profile of the Dmrt (Doublesex and Mab-3 Related Transcription Factor) Gene Family in Channel Catfish (Ictalurus punctatus)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM domain related to sex determination and differentiation, which has been systematically described in various teleost fish, but less in channel catfish (Ictalurus punctatus), an important global aquaculture species in the US and China. In this study, seven Dmrt genes from channel catfish genome were identified and analyzed using bioinformatics methods. Seven IpDmrt genes were distributed unevenly across five chromosomes. Synteny analysis revealed that Dmrt1, Dmrt2a, Dmrt2b, Dmrt3, Dmrt4, and Dmrt5 were relatively conserved in teleost fish. Tissue distribution analysis showed that IpDmrt1, IpDmrt2b, IpDmrt5, and IpDmrt6 exhibited sexually dimorphic expression patterns and, among them, IpDmrt1 and IpDmrt6 had high expression levels in the testes, while IpDmrt2b and IpDmrt5 had more significant expression levels in the ovaries than in other tissues. After 17β-estradiol treatment, IpDmrt2b and IpDmrt5 were significantly up regulated, while the expression of IpDmrt1 and IpDmrt6 was significantly repressed in XY channel catfish ovaries compared with XX channel catfish ovaries. The present study provides a comprehensive insight into the Dmrt gene family of channel catfish. The results suggest that IpDmrt1 and IpDmrt6 may play an important role in testis differentiation/development, while IpDmrt2b and IpDmrt5 are critical in ovary development in this species.

Functional Differentiation of BMP7 Genes in Zebrafish: bmp7a for Dorsal-Ventral Pattern and bmp7b for Melanin Synthesis and Eye Development

Bone morphogenetic protein 7 (BMP7) belongs to the transforming growth factor β (TGF-β) family, which not only induces cartilage and bone formation, but also regulates eye development and melanoma tumorigenesis in mammals. In teleosts, BMP7 differentiates into two subtypes, bmp7a and bmp7b, which have clearly differentiated structures. To fully understand the functional differentiation of bmp7a and bmp7b in fish species, we successfully constructed bmp7a and bmp7b gene deletion mutants in zebrafish using CRISPR/Cas9-mediated gene editing technology. Our results showed that bmp7a mutation caused abnormal development of the embryo’s dorsal-ventral pattern that led to death; bmp7b mutation induced growth inhibition and increased melanin production in the skin and eye of mutants. Histological analysis revealed that melanin in the retina of the eyes in bmp7b mutants increased, and behavioral observation showed that the vision and sensitivity to food of the mutants were reduced. Transcriptome analysis of the skin and eye tissues showed that the expression changes of wnt7ba and gna14 in bmp7b mutants might promote the increase of melanin. Additionally, the eye transcriptome analysis indicated that changes in the structure of the eyes in bmp7b mutants led to defects in phototransduction, and seven DEGs (rgs9a, rgs9b, rcvrn2, guca1d, grk1b, opn1mw4, and gc2) were identified as key candidate genes that affected the photonic response of the eyes. The study revealed the functional differentiation of bmp7a and bmp7b in teleosts and the first report about the inhibitory effect of bmp7b on melanogenesis may provide useful information for the future research on human melanoma-related diseases.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Comparative Genomics Studies on the dmrt Gene Family in Fish

Doublesex and mab-3-related transcription factor (dmrt) genes are widely distributed across various biological groups and play critical roles in sex determination and neural development. Here, we applied bioinformatics methods to exam cross-species changes in the dmrt family members and evolutionary relationships of the dmrt genes based on genomes of 17 fish species. All the examined fish species have dmrt1–5 while only five species contained dmrt6. Most fish harbored two dmrt2 paralogs (dmrt2a and dmrt2b), with dmrt2b being unique to fish. In the phylogenetic tree, 147 DMRT are categorized into eight groups (DMRT1–DMRT8) and then clustered in three main groups. Selective evolutionary pressure analysis indicated purifying selections on dmrt1–3 genes and the dmrt1–3–2(2a) gene cluster. Similar genomic conservation patterns of the dmrt1–dmrt3–dmrt2(2a) gene cluster with 20-kb upstream/downstream regions in fish with various sex-determination systems were observed except for three regions with remarkable diversity. Synteny analysis revealed that dmrt1, dmrt2a, dmrt2b, and dmrt3–5 were relatively conserved in fish during the evolutionary process. While dmrt6 was lost in most species during evolution. The high conservation of the dmrt1–dmrt3–dmrt2(2a) gene cluster in various fish genomes suggests their crucial biological functions while various dmrt family members and sequences across fish species suggest different biological roles during evolution. This study provides a molecular basis for fish dmrt functional analysis and may serve as a reference for in-depth phylogenomics.

Genome-wide investigation of Dmrt gene family in large yellow croaker (Larimichthys crocea)

The Dmrt (Doublesex and Mab-3 related transcription factor) gene family is a class of crucial transcription factors characterized by a conserved DM (Doublesex/Mab-3) domain. Previous researches indicate this gene family is involved in various physiological processes, especially in sex determination/differentiation and gonad development. Despite the vital roles of the Dmrt gene family in physiological processes, the comprehensive characterization and analysis of the dmrt genes in large yellow croaker (Larimichthys crocea), one of the most commercially important marine fish in China, have not been described. In this study, we performed the first genome-wide systematic analysis of L. crocea dmrt genes through the bioinformatics method. A total of seven members of the Dmrt gene family including Lcdmrt1, Lcdmrt2a, Lcdmrt2b, Lcdmrt3, Lcdmrt4, Lcdmrt5, and Lcdmrt6 were excavated based on the genome data of L. crocea. Further analysis revealed that the dmrt genes of L. crocea were distributed unevenly across four chromosomes. There were three dmrt genes (Lcdmrt1, Lcdmrt2a, and Lcdmrt3) on 3rd chromosome, one (Lcdmrt6) on 13th chromosome, one (Lcdmrt4) on 14th chromosome, two on (Lcdmrt5 and Lcdmrt2b) 17th chromosome. The gene structure analysis indicated that the number of introns of different dmrt genes of L. crocea had some differences: Lcdmrt1 had four introns, Lcdmrt2a, Lcdmrt2b, and Lcdmrt6 had two introns, Lcdmrt3, Lcdmrt4, and Lcdmrt5 had only one intron. The expression pattern analysis with published gonad transcriptome datasets and further confirmed by qRT-PCR revealed that these members of the Dmrt gene family except for Lcdmrt4 were all sexually dimorphic and preferred expressing in testis. Furthermore, the expression pattern analysis also revealed that the expression level of Lcdmrt1 and Lcdmrt6 was significantly higher than that of other members, suggesting that these two genes may play a more important role in testis. Overall, our studies provide a comprehensive insight into the Dmrt gene family members and a basis for the further study of their biological functions in L. crocea.

Behavioral Characterization of dmrt3a Mutant Zebrafish Reveals Crucial Aspects of Vertebrate Locomotion through Phenotypes Related to Acceleration

Vertebrate locomotion is orchestrated by spinal interneurons making up a central pattern generator. Proper coordination of activity, both within and between segments, is required to generate the desired locomotor output. This coordination is altered during acceleration to ensure the correct recruitment of muscles for the chosen speed. The transcription factor Dmrt3 has been proposed to shape the patterned output at different gaits in horses and mice.

Vertebrate locomotion is orchestrated by spinal interneurons making up a central pattern generator. Proper coordination of activity, both within and between segments, is required to generate the desired locomotor output. This coordination is altered during acceleration to ensure the correct recruitment of muscles for the chosen speed. The transcription factor Dmrt3 has been proposed to shape the patterned output at different gaits in horses and mice. Here, we characterized dmrt3a mutant zebrafish, which showed a strong, transient, locomotor phenotype in developing larvae. During beat-and-glide swimming, mutant larvae showed fewer and shorter movements with decreased velocity and acceleration. Developmental compensation likely occurs as the analyzed behaviors did not differ from wild-type at older larval stages. However, analysis of maximum swim speed in juveniles suggests that some defects persist within the mature locomotor network of dmrt3a mutants. Our results reveal the pivotal role Dmrt3 neurons play in shaping the patterned output during acceleration in vertebrates.

Loss of Dmrt5 Affects the Formation of the Subplate and Early Corticogenesis

Dmrt5 (Dmrta2) and Dmrt3 are key regulators of cortical patterning and progenitor proliferation and differentiation. In this study, we show an altered apical to intermediate progenitor transition, with a delay in SP neurogenesis and premature birth of Ctip2⁺ cortical neurons in Dmrt5^−/− mice. In addition to the cortical progenitors, DMRT5 protein appears present in postmitotic subplate (SP) and marginal zone neurons together with some migrating cortical neurons. We observed the altered split of preplate and the reduced SP and disturbed radial migration of cortical neurons into cortical plate in Dmrt5^−/− brains and demonstrated an increase in the proportion of multipolar cells in primary neuronal cultures from Dmrt5^−/− embryonic brains. Dmrt5 affects cortical development with specific time sensitivity that we described in two conditional mice with slightly different deletion time. We only observed a transient SP phenotype at E15.5, but not by E18.5 after early (Dmrt5^lox/lox;Emx1^Cre), but not late (Dmrt5^lox/lox;Nestin^Cre) deletion of Dmrt5. SP was less disturbed in Dmrt5^lox/lox;Emx1^Cre and Dmrt3^−/− brains than in Dmrt5^−/− and affects dorsomedial cortex more than lateral and caudal cortex. Our study demonstrates a novel function of Dmrt5 in the regulation of early SP formation and radial cortical neuron migration.

Summary Statement

Our study demonstrates a novel function of Dmrt5 in regulating marginal zone and subplate formation and migration of cortical neurons to cortical plate.

Sex-dependent dominance maintains migration supergene in rainbow trout

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

Genome-Wide identification of doublesex and Mab-3-Related transcription factor (DMRT) genes in nile tilapia (oreochromis niloticus)

Doublesex and Mab-3-related transcription factor (DMRT) gene family is extensively known for its contribution in sex determination and differentiation across phyla. Here we report the identification of five DM (doublesex and mab-3) domain genes in the Nile tilapia which includes DMRT1, DMRTa2, DMRT2a, DMRT2b and DMRT3a. The full-length sequence of DMRT genes ranges from 3526 (DMRTA2) to 1471bp (DMRT1) which encode putative proteins series from 469 to 372 amino acids. All the DMRT proteins contained at least one conserved DNA-binding DM domain. Sub-cellular localization and gene ontology revealed DMRT1 protein is maximum localized in nuclear region and gene ontology analysis showed the molecular function of 48.2%, biological process 43.6% and cellular component 25%. Chromosomal location and synteny analysis displayed that DMRT genes mostly cluster linkage group 12. Altogether, our findings provide vital genomic information for future studies of biochemical, physiological, and phylogenetic studies on DMRT genes in teleost.

Independent evolution for sex determination and differentiation in the DMRT family in animals

Some DMRT family genes including arthropod dsx, nematode mab-3, and vertebrate dmrt1 are involved in sex determination and/or differentiation in bilaterian animals. Although there have been some reports about evolutionary analyses of the family by using its phylogenetic trees, it is still undecided as to whether these three sex determination-related genes share orthologous relationships or not. To clarify this question, we analyzed evolutional relationships among the family members in various bilaterians by using not only phylogenetic tree analysis, but also synteny analysis. We found that only four genes, dmrt2a/2b, dmrt3, dmrt4/5 and dmrt93B were commonly present in invertebrate bilateria. The syntenies of dmrt2a/2b-dmrt3 and dmrt4/5-dmrt93B are conserved before and after two rounds of whole genome duplication in the ancestral vertebrate. Importantly, this indicates that dmrt1 must have appeared in the common vertebrate ancestor. In addition, dmrt1, dsx, or mab-3 formed each different cluster at a distance in our phylogenetic tree. From these findings, we concluded that the three sex determination-related genes, dmrt1, dsx, and mab-3 have no orthologous relationships, and suggested independent evolution for sex determination and differentiation in the DMRT gene family. Our results may supply clues about why sex-determining systems have diverged during animal evolution.

Summary: Three DMRT family genes, vertebrate dmrt1, arthropod dsx and nematode mab-3, involved in sex determination and primary sex differentiation have no orthologous relationships, indicating independent evolution in bilaterian animals.

Phylogenetic analysis and embryonic expression of panarthropod Dmrt genes

Background

One set of the developmentally important Doublesex and Male-abnormal-3 Related Transcription factors (Dmrt) is subject of intense research, because of their role in sex-determination and sexual differentiation. This likely non-monophyletic group of Dmrt genes is represented by the Drosophila melanogaster gene Doublesex (Dsx), the Caenorhabditis elegans Male-abnormal-3 (Mab-3) gene, and vertebrate Dmrt1 genes. However, other members of the Dmrt family are much less well studied, and in arthropods, including the model organism Drosophila melanogaster, data on these genes are virtually absent with respect to their embryonic expression and function.

Results

Here we investigate the complete set of Dmrt genes in members of all main groups of Arthropoda and a member of Onychophora, extending our data to Panarthropoda as a whole. We confirm the presence of at least four families of Dmrt genes (including Dsx-like genes) in Panarthropoda and study their expression profiles during embryogenesis. Our work shows that the expression patterns of Dmrt11E, Dmrt93B, and Dmrt99B orthologs are highly conserved among panarthropods. Embryonic expression of Dsx-like genes, however, is more derived, likely as a result of neo-functionalization after duplication.

Conclusions

Our data suggest deep homology of most of the panarthropod Dmrt genes with respect to their function that likely dates back to their last common ancestor. The function of Dsx and Dsx-like genes which are critical for sexual differentiation in animals, however, appears to be much less conserved.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0322-0) contains supplementary material, which is available to authorized users.

Expression profiles of dmrts and foxls during gonadal development and sex reversal induced by 17α-methyltestosterone in the orange-spotted grouper

The orange-spotted grouper, Epinephelus coioides, is a marine protogynous hermaphrodite fish of commercial importance. There are many examples of sex change species among marine fish, but the molecular basis for the sex change is still unknown. Gonadal expression patterns of the dmrts and foxls genes in E. coioides have pointed to sexual dimorphism in this species and it has been shown that mRNA levels of dmrts and foxls to vary significantly during reproduction cycles. The steroid 17α-methyltestosterone was used to induce sex reversal in these fish, during which dmrts and foxls levels changed significantly and subsequently reverted to normal when 17α-methyltestosterone was withdrawn. Interestingly, the expression of dmrt2b and dmrt3 was not affected by this steroid. We speculate that the role of foxl2 in reproduction may be conserved via regulation of early differentiation of the ovary by the hypothalamus-pituitary-gonad axis, and dmrt2 may have a significant role in premature ovarian differentiation and maintenance in E. coioides. dmrt1 and foxl3 played a role in the development of the testes and are believed to be potential male regulatory genes.

Endodermal pouch-expressed dmrt2b is important for pharyngeal cartilage formation

Pharyngeal pouches, a series of outpocketings derived from the foregut endoderm, are essential for craniofacial skeleton formation. However, the molecular mechanisms underlying endodermal pouch-regulated head cartilage development are not fully understood. In this study, we find that zebrafish dmrt2b, a gene encoding Doublesex- and Mab-3-related transcription factor, is specifically expressed in endodermal pouches and required for normal pharyngeal cartilage development. Loss of dmrt2b doesn't affect cranial neural crest (CNC) specification and migration, but leads to prechondrogenic condensation defects by reducing cxcl12b expression after CNC cell movement into the pharyngeal arches. Moreover, dmrt2b inactivation results in reduced proliferation and impaired differentiation of CNC cells. We also show that dmrt2b suppresses crossveinless 2 expression in endodermal pouches to maintain BMP/Smad signaling in the arches, thereby facilitating CNC cell proliferation and chondrogenic differentiation. This work provides insight into how transcription factors expressed in endodermal pouches regulate pharyngeal skeleton development through tissue-tissue interactions.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Genomic characterization of the Atlantic cod sex-locus

A variety of sex determination mechanisms can be observed in evolutionary divergent teleosts. Sex determination is genetic in Atlantic cod (Gadus morhua), however the genomic location or size of its sex-locus is unknown. Here, we characterize the sex-locus of Atlantic cod using whole genome sequence (WGS) data of 227 wild-caught specimens. Analyzing more than 55 million polymorphic loci, we identify 166 loci that are associated with sex. These loci are located in six distinct regions on five different linkage groups (LG) in the genome. The largest of these regions, an approximately 55 Kb region on LG11, contains the majority of genotypes that segregate closely according to a XX-XY system. Genotypes in this region can be used genetically determine sex, whereas those in the other regions are inconsistently sex-linked. The identified region on LG11 and its surrounding genes have no clear sequence homology with genes or regulatory elements associated with sex-determination or differentiation in other species. The functionality of this sex-locus therefore remains unknown. The WGS strategy used here proved adequate for detecting the small regions associated with sex in this species. Our results highlight the evolutionary flexibility in genomic architecture underlying teleost sex-determination and allow practical applications to genetically sex Atlantic cod.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

A Novel Female-Specific and Sexual Reproduction-Associated Dmrt Gene Discovered in the Stony Coral, Euphyllia ancora

Transcription factors encoded by the Dmrt gene family regulate multiple aspects of animal reproduction. Most studies investigating the Dmrt gene family were conducted in model organisms from bilateral species, with a particular emphasis on gene function in male sex determination. It is still unclear whether the E. ancora Dmrt (EaDmrt) genes found in basal metazoans such as cnidarians share similar characteristics with orthologs in other metazoans. In this study, seven full Dmrt gene transcript sequences for a gonochoric coral, Euphyllia ancora (phylum: Cnidaria; class: Anthozoa), were obtained through transcriptome data mining, RT-PCR analysis, rapid amplification of cDNA ends, and sequencing. These EaDmrts were subjected to quantitative assays measuring temporal and tissue-specific expression. Results demonstrated a unique gene expression pattern for EaDmrtE, which is enriched in female germ cells during the spawning season. Based on the phylogenetic analyses performed across the homologous Dmrt genes in metazoans, we found that the female-specific EaDmrtE gene is not related to the DM1 gene of Acropora spp. coral nor to Dmrt1 of vertebrates, which are involved in sexual reproduction, especially in sex determination (vertebrate Dmrt1). Additionally, high levels of EaDmrtE transcripts detected in unfertilized mature eggs are retained in newly formed zygotes but decrease during embryonic development. We suggest that the newly discovered gene may play a role in oogenesis and early embryogenesis as a maternal factor in corals. Therefore, the sexual reproduction-associated Dmrt gene(s) should have arisen in cnidarians and might have evolved multiple times in metazoans.

Effects of crude oil exposure and elevated temperature on the liver transcriptome of polar cod (Boreogadus saida)

Petroleum-related activities in the Arctic have raised concerns about the adverse effects of potential oil spill on the environment and living organisms. Polar cod plays a key role in the Arctic marine ecosystem and is an important species for monitoring oil pollution in this region. We examined potential interactions of oil pollution and global warming by analysing liver transcriptome changes in polar cod exposed to crude oil at elevated temperature. Adult males and females were kept at high (11°C) or normal (4°C) temperature for 5 days before exposure to mechanically dispersed crude oil for 2 days followed by recovery in clean sea water for 11 days at the two temperatures. Genome-wide microarray analysis of liver samples revealed numerous differentially expressed genes induced by uptake of oil as confirmed by increased levels of bile polycyclic aromatic hydrocarbon (PAH) metabolites. The hepatic response included genes playing important roles in xenobiotic detoxification and closely related biochemical processes, but also of importance for protein stress response, cell repair and immunity. Though magnitude of transcriptome responses was similar at both temperatures, the upregulated expression of cyp1a1 and several chaperone genes was much stronger at 11°C. Most gene expression changes returned to basal levels after recovery. The microarray results were validated by qPCR measurement of eleven selected genes representing both known and novel biomarkers to assess exposure to anthropogenic threats on polar cod. Strong upregulation of the gene encoding fibroblast growth factor 7 is proposed to protect the liver of polar fish with aglomerular kidneys from the toxic effect of accumulated biliary compounds. The highly altered liver transcriptome patterns after acute oil exposure and recovery suggests rapid responses in polar cod to oil pollutants and the ability to cope with toxicity in relatively short time.

Isolation of doublesex- and mab-3-related transcription factor 6 and its involvement in spermatogenesis in tilapia

The dmrt6 gene has been isolated from tetrapods and recently from a coelacanth, Latimeria chalumnae. Its evolutionary history and exact function remain unclear. In the present study, dmrt6 was isolated from Perciformes (five cichlids and stickleback), Siluriformes (southern catfish), and Lepisosteiformes (spotted gar). Syntenic and phylogenetic analyses indicated that dmrt6 experienced gene transposition after the divergence of teleosts from other bony fish as gene loci surrounding dmrt6 were conserved among teleosts (but was completely different from gene loci surrounding dmrt6 in tetrapods and spotted gar), while these gene loci were conserved among nonteleost species. Real-time PCR and in situ hybridization revealed that dmrt6 was highly expressed in the XY gonads from 90 days after hatching (dah) onward and was observed exclusively in spermatocytes of the testes in tilapia. Dmrt6 knockout by CRISPR/Cas9 resulted in fewer spermatocytes, down-regulated Cyp11b2 in testes, and consequently produced a lower level of serum 11-ketotestosterone (11-KT) in Dmrt6-deficient XY fish compared with the XY control at 120 dah. From 150 to 180 dah, spermatogenesis gradually recovered, and cyp11b2 expression and serum 11-KT level were restored to the same levels as those of the XY control fish. In addition, a Dmrt6 mutation was observed in genomic DNA of sperm of G0 mutant fish and F1 fish. Taken together, our data suggest that dmrt6 also exists in bony fish. Its absence in most fish genomes was probably due to incomplete sequencing and/or secondary loss. The dmrt6 gene is highly expressed in spermatocytes and is involved in spermatogenesis in tilapia.

Molecular cloning of Pcc-dmrt1s and their specific expression patterns in Pengze crucian carp (Carassius auratus var. Pengze) affected by 17α-methyltestosterone

Dmrt1, an important transcription factor associated with testicular differentiation, is conserved among teleost, which could also be detected in ovaries. In the present study, three isoforms of Pcc-dmrt1s (Pcc-dmrt1a, Pcc-dmrt1b and Pcc-dmrt1c) resulting from alternative splicing of the dmrt1 gene were cloned and characterized in the triploid gynogenetic fish, the Pengze crucian carp. Their mRNA expression profiling was investigated in juvenile developmental stages, tissues of the adult fish, and the juveniles under 84.2 ng/L 17α-methyltestosterone (MT) treatments. Results showed that their putative proteins shared high identities to Dmrt1 in cyprinid fish species. Gene expression profiling in the developmental stages showed that all the three target genes had a highest/lowest expression at 56/40 days post-hatching (dph), respectively. The period of 40 dph appeared to be a key time during the process of the ovary development of Pengze crucian carp. The tissue distribution results indicated that Pcc-dmrt1s were predominantly expressed in hepatopancreas, brain, spleen and ovary of the female fish. MT significantly increased the mRNA expression of Pcc-dmrt1a (all 4-week exposures) and Pcc-dmrt1b (except for week 2), while repressed Pcc-dmrt1c transcripts at all exposure period except for week 2. MT extremely significant repressed cyp19a1a transcripts for 1 week. The present study indicated that MT could influence the ovary development of Pengze crucian carp by disturbing gene expressions of Pcc-dmrt1s and cyp19a1a. Furthermore, the present study will be of great significance to broaden the understanding of masculinizing pathway during ovary development in gynogenetic teleost.

Genome-wide microarray analysis of Atlantic cod (Gadus morhua) oocyte and embryo

Background

Regulation of gene expression plays a central role in embryonic development. Early stages are controlled by gametic transcripts, which are subsequently substituted with transcripts from the genome of the zygote. Transcriptomic analyses provide an efficient approach to explore the temporal gene expression profiles in embryos and to search for the developmental regulators. We report a study of early Atlantic cod development that used a genome-wide oligonucleotide microarray to examine the composition and putative roles of polyadenylated transcripts.

Results

The analyses were carried out in unfertilized oocytes, newly fertilized oocytes and embryos at the stages of mid-blastula transition and segmentation. Numerous genes transcribed in oocytes are involved in multiple aspects of cell maintenance and protection, including metabolism, signal perception and transduction, RNA processing, cell cycle, defense against pathogens and DNA damage. Transcripts found in unfertilized oocytes also encoded a large number of proteins implicated in cell adherence, tight junction and focal adhesion, suggesting high complexity in terms of structure and cellular interactions in embryos prior to midblastula transition (MBT). Prezygotic transcripts included multiple regulators that are most likely involved in developmental processes that take place long after fertilization, such as components of ErbB, hedgehog, notch, retinoid, TGFb, VEGF and Wnt signaling pathways, as well as transcripts involved in the development of nervous system. The major event of MBT was the activation of a large group of histones and other genes that modify chromatin structure preceding massive gene expression changes. A hallmark of events observed during segmentation was the induction of multiple transcription factors, including a large group of homeobox proteins in pace with decay of a large fraction of maternal transcripts. Microarray analyses detected a suite of master developmental regulators that control differentiation and maintenance of diverse cell lineages.

Conclusions

Transcriptome profiling of the early stages in Atlantic cod revealed the presence of transcripts involved in patterning and development of tissues and organs long before activation of the zygotic genome. The switch from maternal to zygotic developmental programs is associated with large-scale modification of chromosomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-594) contains supplementary material, which is available to authorized users.

Sex determination in annual fishes: Searching for the master sex-determining gene in Austrolebias charrua (Cyprinodontiformes, Rivulidae)

Evolution of sex determination and differentiation in fishes involves a broad range of sex strategies (hermaphroditism, gonochorism, unisexuality, environmental and genetic sex determination). Annual fishes inhabit temporary ponds that dry out during the dry season when adults die. The embryos exhibit an atypical developmental pattern and remain buried in the bottom mud until the next rainy season. To elucidate genomic factors involved in the sex determination in annual fish, we explored the presence of a candidate sex-specific gene related to the cascade network in Austrolebias charrua. All phylogenetic analyses showed a high posterior probability of occurrence for a clade integrated by nuclear sequences (aprox. 900 bp) from both adults (male and female), with partial cDNA fragments of A. charrua from juveniles (male) and the dsx D. melanogaster gene. The expressed fragment was detected from blastula to adulthood stages showing a sexually dimorphic expression pattern. The isolated cDNA sequence is clearly related to dsx D. melanogaster gene and might be located near the top of the sex determination cascade in this species.

Identification of Dmrt genes and their up-regulation during gonad transformation in the swamp eel (Monopterus albus)

The swamp eel is a teleost fish with a characteristic of natural sex reversal and an ideal model for vertebrate sexual development. However, underlying molecular mechanisms are poorly understood. We report the identification of five DM (doublesex and mab-3) domain genes in the swamp eel that include Dmrt2, Dmrt2b, Dmrt3, Dmrt4 and Dmrt5, which encode putative proteins of 527, 373, 471, 420 and 448 amino acids, respectively. Phylogenetic tree showed that these genes are clustered into corresponding branches of the DM genes in vertebrates. Southern blot analysis indicated that the Dmrt1-Dmrt3-Dmrt2 genes are tightly linked in a conserved gene cluster. Notably, these Dmrt genes are up-regulated during gonad transformation. Furthermore, mRNA in situ hybridisation showed that Dmrt2, Dmrt3, Dmrt4 and Dmrt5 are expressed in developing germ cells. These results are evidence that the DM genes are involved in sexual differentiation in the swamp eel.

Characterization of sex determination and sex differentiation genes in Latimeria

Genes involved in sex determination and differentiation have been identified in mice, humans, chickens, reptiles, amphibians and teleost fishes. However, little is known of their functional conservation, and it is unclear whether there is a common set of genes shared by all vertebrates. Coelacanths, basal Sarcopterygians and unique "living fossils", could help establish an inventory of the ancestral genes involved in these important developmental processes and provide insights into their components. In this study 33 genes from the genome of Latimeria chalumnae and from the liver and testis transcriptomes of Latimeria menadoensis, implicated in sex determination and differentiation, were identified and characterized and their expression levels measured. Interesting findings were obtained for GSDF, previously identified only in teleosts and now characterized for the first time in the sarcopterygian lineage; FGF9, which is not found in teleosts; and DMRT1, whose expression in adult gonads has recently been related to maintenance of sexual identity. The gene repertoire and testis-specific gene expression documented in coelacanths demonstrate a greater similarity to modern fishes and point to unexpected changes in the gene regulatory network governing sexual development.