Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish

Development and Validation of Sex-Specific Markers in Pelodiscus Sinensis Using Restriction Site-Associated DNA Sequencing

The sex of an animal influences its economic traits, especially in species displaying sexual dimorphism. The Chinese soft-shelled turtle, Pelodiscus sinensis, is an economically important aquatic species that shows significant male sexual dimorphism, with a large body size, faster growth, a thick and wide calipash, and lower body fat. In this study, ten male and ten female turtles were subjected to restriction site-associated DNA sequencing (RAD-seq) using the Hi-Seq 4000 sequencing platform to isolate female-specific DNA fragments. We identified 5967 bp and 6532 bp fragments using genome walking. Three female-specific markers designed from these two fragments were confirmed to separate the sexes of Pelodiscus sinensis perfectly. One of the female-specific markers showed dosage association in female and male individuals. Individuals from different populations (n = 296) were used to validate that the female-specific markers could identify the genetic sex of Pelodiscus sinensis with 100% accuracy. The results of the present study demonstrated that RAD-seq was useful to develop sex-related markers in animals, and verified that the sex determination system of Pelodiscus sinensis belonged to the ZZ/ZW heterogametic system. Importantly, the developed markers could lead to a method for sex-controlled breeding in the Chinese soft-shelled turtle.

A High-Density Genetic Linkage Map and QTL Mapping for Sex in Black Tiger Shrimp (Penaeus monodon)

The black tiger shrimp, Penaeus monodon, is important in both fishery and aquaculture and is the second-most widely cultured shrimp species in the world. However, the current strains cannot meet the market needs in various cultural environments, and the genome resources for P. monodon are still lacking. Restriction-site associated DNA sequencing (RADseq) has been widely used in genetic linkage map construction and in quantitative trait loci (QTL) mapping. We constructed a high-density genetic linkage map with RADseq in a full-sib family. This map contained 6524 single nucleotide polymorphisms (SNPs) and 2208 unique loci. The total length was 3275.4 cM, and the genetic distance was estimated to be 1.1 Mb/cM. The sex trait is a dichotomous phenotype, and the same interval was detected as a QTL using QTL mapping and genome-wide association analysis. The most significant locus explained 77.4% of the phenotype variance. The sex locus was speculated to be the same in this species based on the sequence alignments in Mozambique, India, and Hawaii populations. The constructed genetic linkage map provided a valuable resource for QTL mapping, genome assembly, and genome comparison for shrimp. The demonstrated common sex locus is a step closer to locating the underlying gene.

Delayed elimination of paternal mtDNA in the interspecific hybrid of Pelteobagrus fulvidraco and Pelteobagrus vachelli during early embryogenesis

Mitochondrial homoplasmy is essential for normal development, as its heteroplasmy usually leads to abnormal or diseased phenotypes in mammals. So far, diverse mechanisms have been proposed to play roles in ensuring uniparental inheritance of mitochondria in many organisms. In recent years, hybrid yellow catfish from mating female yellow catfish (Pelteobagrus fulvidraco) with male darkbarbel catfish (Pelteobagrus vachelli) has been widely cultured in China due to its fast-growing. However, a high rate of abnormal and defective embryos was observed in the offsprings of hybrid yellow catfish. In this study, we systematically investigated the elimination process of paternal mitochondrial DNA (mtDNA) in yellow catfish and hybrid yellow catfish. The mtDNA contents significantly decreased in the isolated mature sperm compared with the semen. Different from the elimination of paternal mtDNA after fertilization in yellow catfish, paternal mtDNA was retained in the developmental embryos of hybrid yellow catfish as later as gastrula stage, indicating a delay of elimination for paternal mtDNA and mitochondrial heteroplasmy during embryogenesis in hybrid yellow catfish. Altogether, the present study suggests that mitochondrial heteroplasmy may affect embryonic development of hybrid progeny between catfish species.

Abundance of Early Embryonic Primordial Germ Cells Promotes Zebrafish Female Differentiation as Revealed by Lifetime Labeling of Germline

Teleost sex differentiation largely depends on the number of undifferentiated germ cells. Here, we describe the generation and characterization of a novel transgenic zebrafish line, Tg(piwil1:egfp-UTRnanos3)^ihb327Tg, which specifically labels the whole lifetime of germ cells, i.e., from primordial germ cells (PGCs) at shield stage to the oogonia and early stage of oocytes in the ovary and to the early stage of spermatogonia, spermatocyte, and spermatid in the testis. By using this transgenic line, we carefully observed the numbers of PGCs from early embryonic stage to juvenile stage and the differentiation process of ovary and testis. The numbers of PGCs became variable at as early as 1 day post-fertilization (dpf). Interestingly, the embryos with a high amount of PGCs mainly developed into females and the ones with a low amount of PGCs mainly developed into males. By using transient overexpression and transgenic induction of PGC-specific bucky ball (buc), we further proved that induction of abundant PGCs at embryonic stage promoted later ovary differentiation and female development. Taken together, we generate an ideal transgenic line Tg(piwil1:egfp-UTRnanos3)^ihb327Tg which can visualize zebrafish germline for a lifetime, and we have utilized this line to study germ cell development and gonad differentiation of teleost and to demonstrate that the increase of PGC number at embryonic stage promotes female differentiation.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

A novel PDZ domain-containing gene is essential for male sex differentiation and maintenance in yellow catfish (Pelteobagrus fulvidraco)

The sex-determining genes are found to be variable among different fish species. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish species in China with XX/XY sex-determining type. Recently, YY super-male yellow catfish has been successfully produced by combining hormonal-induced sex reversal method with sex chromosome-linked markers. Here, we identified a novel PDZ domain-containing gene in yellow catfish designated as pfpdz1, in whose intron the sex-linked marker was located. The coding sequence of pfpdz1 in Y chromosome was identical to that in X chromosome except a missense SNP (A/T) that changes an amino acid (E8V) in the N-terminal region. Pfpdz1 displayed male-specific expression during sex differentiation. Overexpression of pfpdz1 using additive transgenesis induces XX ovary to differentiate into testis-like tissue, while the targeted inactivation of pfpdz1 in Y chromosome using CRISPR/Cas9-mediated mutagenesis triggers ovarian differentiation. Furthermore, we demonstrated that pfpdz1 initiates testicular differentiation through upregulating expression of amh, dmrt1 and sox9a1, as well as downregulating expression of cyp19a1, foxl2 and wnt4. Our data provide functional evidence that pfpdz1 is significant for male differentiation and maintenance in yellow catfish.

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.

Chromosomal-level assembly of yellow catfish genome using third-generation DNA sequencing and Hi-C analysis

Background

The yellow catfish, Pelteobagrus fulvidraco, belonging to the Siluriformes order, is an economically important freshwater aquaculture fish species in Asia, especially in Southern China. The aquaculture industry has recently been facing tremendous challenges in germplasm degeneration and poor disease resistance. As the yellow catfish exhibits notable sex dimorphism in growth, with adult males about two- to three-fold bigger than females, the way in which the aquaculture industry takes advantage of such sex dimorphism is another challenge. To address these issues, a high-quality reference genome of the yellow catfish would be a very useful resource.

Findings

To construct a high-quality reference genome for the yellow catfish, we generated 51.2 Gb short reads and 38.9 Gb long reads using Illumina and Pacific Biosciences (PacBio) sequencing platforms, respectively. The sequencing data were assembled into a 732.8 Mb genome assembly with a contig N50 length of 1.1 Mb. Additionally, we applied Hi-C technology to identify contacts among contigs, which were then used to assemble contigs into scaffolds, resulting in a genome assembly with 26 chromosomes and a scaffold N50 length of 25.8 Mb. Using 24,552 protein-coding genes annotated in the yellow catfish genome, the phylogenetic relationships of the yellow catfish with other teleosts showed that yellow catfish separated from the common ancestor of channel catfish ∼81.9 million years ago. We identified 1,717 gene families to be expanded in the yellow catfish, and those gene families are mainly enriched in the immune system, signal transduction, glycosphingolipid biosynthesis, and fatty acid biosynthesis.

Conclusions

Taking advantage of Illumina, PacBio, and Hi-C technologies, we constructed the first high-quality chromosome-level genome assembly for the yellow catfish P. fulvidraco. The genomic resources generated in this work not only offer a valuable reference genome for functional genomics studies of yellow catfish to decipher the economic traits and sex determination but also provide important chromosome information for genome comparisons in the wider evolutionary research community.

Effects of Temperature on the Expression of Two Ovarian Differentiation-Related Genes foxl2 and cyp19a1a

Exposure to stress induces a series of responses and influences a wide range of biological processes including sex differentiation in fish. The present work investigated the molecular and physiological response to thermal stress throughout the early development stage covering the whole period of sex differentiation of bluegill, Lepomis macrochirus. Larvae were treated using three temperatures, 17, 24, and 32°C from 6 to 90 days posthatching (dph) in 30-L round tanks. There is no significant difference of the sex ratio and survival among the three temperature groups in the geographic population used in this study. Two ovarian differentiation-related genes foxl2 and cyp19a1a were detected at 7 dph suggesting that these genes have already played a role prior to sex differentiation. The expression of foxl2 reached the peak and was thermosensitive just prior to the onset of ovarian differentiation at 27 dph. Histological examination displayed that the proliferation of germ cells and ovarian differentiation were delayed at the low-temperature treatment (17°C) at 97 dph compared with higher temperatures. In conclusion, the water temperature regulates the sex differentiation of bluegill through modulation of the expression of foxl2 and cyp19a1a. A comparative study of the expression profile of sex differentiation-related genes in species will shed light on the evolution of sex-determination mechanisms and the impact of stress on sex differentiation.

Whole Genome Incorporation and Epigenetic Stability in a Newly Synthetic Allopolyploid of Gynogenetic Gibel Carp

Allopolyploidization plays an important role in speciation, and some natural or synthetic allopolyploid fishes have been extensively applied to aquaculture. Although genetic and epigenetic inheritance and variation associated with plant allopolyploids have been well documented, the relative research in allopolyploid animals is scarce. In this study, the genome constitution and DNA methylation inheritance in a newly synthetic allopolyploid of gynogenetic gibel carp were analyzed. The incorporation of a whole genome of paternal common carp sperm in the allopolyploid was confirmed by genomic in situ hybridization, chromosome localization of 45S rDNAs, and sequence comparison. Pooled sample-based methylation sensitive amplified polymorphism (MSAP) revealed that an overwhelming majority (98.82%) of cytosine methylation patterns in the allopolyploid were inherited from its parents of hexaploid gibel carp clone D and common carp. Compared to its parents, 11 DNA fragments in the allopolyploid were proved to be caused by interindividual variation, recombination, deletion, and mutation through individual sample-based MSAP and sequencing. Contrast to the rapid and remarkable epigenetic changes in most of analyzed neopolyploids, no cytosine methylation variation was detected in the gynogenetic allopolyploid. Therefore, the newly synthetic allopolyploid of gynogenetic gibel carp combined genomes from its parents and maintained genetic and epigenetic stability after its formation and subsequently seven successive gynogenetic generations. Our current results provide a paradigm for recurrent polyploidy consequences in the gynogenetic allopolyploid animals.

Social factors and aromatase gene expression during adult male-to-female sex change in captive leopard grouper Mycteroperca rosacea

Social factors and aromatase gene expression in the leopard grouper Mycteroperca rosacea was studied when captive fish were separated by sex during the reproductive (April-June) and post-reproductive (July-September) seasons. Monosex females, monosex males, and mixed-sex, held in social sextet units were analyzed for sex steroids throughout confinement. At the end of the experiment, the gonad-sex was defined by histology, and gonad and brain aromatase gene expressions were quantified. Only males held in the monosex social units changed sex. Histology showed one male remained unchanged, six were found in a transitional sexual stage, in which two had intersex-predominantly-testes, and four had a more defined intersex ovo-testes pattern, and 11 were immature de novo females (neofemales). Neofemales and most intersex fish did not survive. In spring, 11-ketosterone showed a specific male profile, which suggests that male-to-female sex change was not triggered during the reproductive season. The low steroid levels in summer made it impossible to associate the sex change to a gonad hormonal shift; in September, gonad aromatase gene expression was not significantly different among groups. However, brain aromatase expression in intersex fish was significantly higher than monosex females, mixed-sex females, and neofemale groups. These results suggest that in the absence of female hormonal compounds, and at a time when male gonad steroidogenesis was diminished, the brain mediated male-to-male social-behavioral interactions, including stress, by increasing aromatization, resulting in derived intersex-male, which triggered more aromatization, followed by a sex change.

Experimental evidence of masculinization by continuous illumination in a temperature sex determination teleost (Atherinopsidae) model: is oxidative stress involved?

The present study evaluates the influence of continuous light on phenotypic sex ratios in Chirostoma estor, a temperature sex determination animal model. Relative gene expression levels of 5 day old larvae were performed on two early gonad differentiation genes (sox9 and foxl2), two stress axis activation genes (gcr1 and crf) and four reactive oxygen species (ROS) antagonist effector genes (sod2, ucp2, gsr and cat). Two light treatments were applied from fertilization; control (12L:12D) simulated natural photoperiod and a continuous illumination photoperiod. By the end of the trial (12 weeks after hatching), differentiated and normal gonads were clearly identifiable in both treatments by histological observations. Regarding sex ratio, 73% of phenotypic males were found in continuous illumination compared with 40% in controls. Consistently, the sox9 gene (involved in early testis differentiation) showed an over expression in 64% of the individual larvae analysed compared with foxl2 (ovarian differentiation) suggesting a masculinization tendency in continuous illumination. On the other hand, only 36% of individuals showed the same tendency in the control treatment consistent with phenotypic sex ratios found under normal culture conditions. Relative gene expression results did not show significant difference in sod2, ucp2 and gcr1 levels, but cat, gsr and crf showed significantly higher expression levels in the continuous illumination treatment suggesting that both, the stress axis and ROS response mechanisms were activated at this time. This study suggests, a link between continuous light, oxidative stress and environmental sex determination in vertebrates. However, further research is necessary to describe this possible upstream mechanism that may drive some aspects of sexual plasticity in vertebrates.

Differential expression and functional diversification of diverse immunoglobulin domain-containing protein (DICP) family in three gynogenetic clones of gibel carp

Diverse immunoglobulin (Ig) domain-containing protein (DICP) family is a novel bony fish-specific multi-gene family encoding diversified immune receptors. However, their function and the implication of binding partners remain unknown. In this study, we first identified 28 DICPs from three gibel carp gynogenetic clones and revealed their high variability and clone-specific feature. After crucian carp herpesvirus (CaHV) infection, these DICPs were significantly upregulated in head kidney, kidney and spleen. The up-regulation folds in clone A⁺, F and H were related to the susceptibility to CaHV, progressively increasing from resistant clone to susceptible clone. Overexpression of gibel carp DICPs inhibited interferon (IFN) and viperin promoter-driven luciferase activity. The additions of E. coli extracts and lipid A significantly enhanced the inhibition effect. In addition, gibel carp DICPs can interact with SHP-1 and SHP-2. These findings suggest that gible carp DICPs, as inhibitory receptors, might specifically recognize lipid A, and then interact with SHP-1 and SHP-2 to inhibit the induction of IFN and ISGs.

Origin and transition of sex determination mechanisms in a gynogenetic hexaploid fish

Most vertebrates reproduce sexually, and plastic sex determination mechanisms including genotypic sex determination (GSD) and environmental sex determination (ESD) have been extensively revealed. However, why sex determination mechanisms evolve diversely and how they correlate with diverse reproduction strategies remain largely unclear. Here, we utilize the superiority of a hexaploid gibel carp (Carassius gibelio) that is able to reproduce by unisexual gynogenesis and contains a rare but diverse proportion of males to investigate these puzzles. A total of 2248 hexaploid specimens were collected from 34 geographic wild populations throughout mainland China, in which 24 populations were revealed to contain 186 males with various incidences ranging from 1.2 to 26.5%. Subsequently, the proportion of temperature-dependent sex determination (TSD) was revealed to be positively correlated to average annual temperature in wild populations, and male incidence in lab gynogenetic progenies was demonstrated to increase with the increasing of larval rearing temperature. Meanwhile, extra microchromosomes were confirmed to play genotypic male determination role as previously reported. Thereby, GSD and TSD were found to coexist in gibel carp, and the proportions of GSD were observed to be much higher than that of TSD in sympatric wild populations. Our findings uncover a potential new mechanism in the evolution of sex determination system in polyploid vertebrates with unisexual gynogenesis ability, and also reveal a possible association of sex determination mechanism transition between TSD and GSD and reproduction mode transition between unisexual gynogenesis and bisexual reproduction.

Evolutionary conservation and divergence of Vasa, Dazl and Nanos1 during embryogenesis and gametogenesis in dark sleeper (Odontobutis potamophila)

Germline-specific genes, Vasa, Dazl and Nanos1, have highly conserved functions in germline development and fertility across animal phyla. In this study, the full-length sequences of Opvasa, Opdazl and Opnanos1 were cloned and characterized from the dark sleeper (Odontobutis potamophila). Gonad-specific expression patterns of Opvasa and Opdazl were confirmed in adult tissues by quantitative real-time PCR (qRT-PCR). Different from Opvasa and Opdazl, the expression of Opnanos1 was ubiquitously detected in all examined tissues except for the liver and spleen. Time-course dynamic expressions during embryogenesis were assessed, and all three genes (Opvasa, Opdazl and Opnanos1) persisted at a high level until gastrulation. qRT-PCR and Western blotting analyses revealed that all three genes were highly expressed throughout gametogenesis. In testis, the expressions of all three genes at the mRNA and protein levels were down-regulated during spermatogenesis. In ovary, different expression patterns were found, and all three genes had a differential role in translational regulation during oogenesis. The expressions of Opvasa, Opdazl and Opnanos1 at the mRNA but not the protein level were high in stage IV. Different expression patterns were found in premeiotic gonads treated by HPG axis hormones (HCG and LHRH-A). Immunolocalization analysis demonstrated that in testis, Opvasa, Opdazl and Opnanos1 were detected in spermatogonia and spermatocytes but absent in the meiotic products, such as spermatids and spermatozoa. In ovary, Opvasa, Opdazl and Opnanos1 persisted at a high level throughout oogenesis. These findings indicated that Opvasa, Opdazl and Opnanos1 played an important role in mitotic and early meiotic phases of oogenesis and spermatogenesis, and they functioned as maternal factors in early embryogenesis. Their proteins could be used as three new markers for germ cells during gametogenesis in O. potamophila gonad. Our data laid a good foundation for improving the breeding efficiency of O. potamophila.

Distinct sperm nucleus behaviors between genotypic and temperature-dependent sex determination males are associated with replication and expression-related pathways in a gynogenetic fish

BACKGROUND

Coexistence and transition of diverse sex determination strategies have been revealed in some ectothermic species, but the variation between males caused by different sex determination strategies and the underlying mechanism remain unclear. Here, we used the gynogenetic gibel carp (Carassius gibelio) with both genotypic sex determination (GSD) and temperature-dependent sex determination (TSD) strategies to illustrate this issue.

RESULTS

We found out that males of GSD and TSD in gibel carp had similar morphology, testicular histology, sperm structure and sperm vitality. However, when maternal individuals were mated with males of GSD, sperm nucleus swelling and fusing with the female pronucleus were observed in the fertilized eggs. On the contrary, when maternal individuals were mated with males of TSD, sperm nucleus remained in the condensed status throughout the whole process. Subsequently, semen proteomics analysis unveiled that DNA replication and gene expression-related pathways were inhibited in the sperm from males of TSD compared to males of GSD, and most differentially expressed proteins associated with DNA replication, transcription and translation were down-regulated. Moreover, via BrdU incorporation and immunofluorescence detection, male nucleus replication was revealed to be present in the fertilized eggs by the sperm from males of GSD, but absent in the fertilized eggs by the sperm from males of TSD.

CONCLUSIONS

These findings indicate that DNA replication and gene expression-related pathways are associated with the distinct sperm nucleus development behaviors in fertilized eggs in response to the sperm from males of GSD and TSD. And this study is the first attempt to screen the differences between males determined via GSD and TSD in gynogenetic species, which might give a hint for understanding evolutionary adaption of diverse sex determination mechanisms in unisexual vertebrates.

Screening and characterisation of sex differentiation-related long non-coding RNAs in Chinese soft-shell turtle (Pelodiscus sinensis)

Long non-coding RNAs (lncRNAs) perform distinct functions in various biological processes in mammals, including sex differentiation. However, the roles of lncRNAs in other vertebrates, especially in the Chinese soft-shell turtle (Pelodiscus sinensis), remain to be clarified. In this study, we performed genome-wide analysis of the lncRNA expression profiles in gonad tissues and screened numerous sex-specific lncRNAs in the Chinese soft-shell turtle. Of the 363,310,650 clean reads obtained, 5,994 sequences were typed as lncRNAs, of which 4,463 were novel. A selection of sex-specific lncRNAs (♀ 932, ♂ 449) from female ovaries and male testis were shown to act on target genes in cis and in trans, and most were involved in gonad differentiation based on Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Furthermore, interactions among the differentially expressed lncRNA-mRNAs and protein coding genes were identified by construction of correlation networks. Overall, our systematic analysis of lncRNA expression profiles in gonad tissues revealed numerous sex-specific lncRNAs in P. sinensis. Thereby, these findings provide new insights into the function of lncRNAs in sex differentiation and highlight a group of candidate lncRNAs for future studies.

Diethylstilbestrol arrested spermatogenesis and somatic growth in the juveniles of yellow catfish (Pelteobagrus fulvidraco), a fish with sexual dimorphic growth

In fish, spermatogenesis and somatic growth are mainly regulated by hypothalamic-pituitary-gonadal (HPG) and hypothalamic-pituitary-somatic (HPS) axes, respectively. Xenoestrogens have been reported to impair spermatogenesis in some fishes, and arrest somatic growth in some others, whereas, whether xenoestrogens are capable of disrupting spermatogenesis and somatic growth simultaneously in fish that exhibits sexual dimorphic growth is little known, and the underlying mechanisms remain poorly understood. In this study, male juveniles of yellow catfish (Pelteobagrus fulvidraco), which exhibits a sexual dimorphic growth that favors males, were exposed to diethylstilbestrol (DES) for 28 days. After exposure, DES significantly disrupted the spermatogenesis (decreased gonadal-somatic index (GSI) and germ cell number) and arrested the somatic growth (declined body weight) of the catfish juveniles. Gene expression and plasma steroid analyses demonstrated the suppressed mRNA levels of genes in HPG axis (gnrh-II, fshβ, and lhβ in the brain and dmrt1, sf1, fshr, cyp17a1, cyp19a1a, and cyp11b2 in the testis) and decreased 17β-estrodial (E2) and 11-ketotestosterone (11-KT) levels in plasma. Further analysis revealed the arrested germ cell proliferation (cyclin d1), meiosis (dmc1, sycp3), and enhanced apoptosis (decreased bcl-2 and elevated bax/bcl-2 ratio) in the testis. Besides, DES also suppressed the mRNA levels of genes in HPS axis (ghrh, gh, and prl in the brain and ghr, igf1, igf2a, and igf2b in the liver). The suppressed HPG and HPS axes were thus supposed to disturb spermatogenesis and arrest somatic growth in yellow catfish. The present study greatly extended our understanding on the mechanisms underlying the toxicity of DES on spermatogenesis and somatic growth of fish.

Stat5b Regulates Sexually Dimorphic Gene Expression in Zebrafish Liver

Sexual size dimorphism is an interesting phenomenon occurred in many fish species. Wildtype zebrafish exhibits a significant sexual dimorphism in body size at the adult stage. Previous studies indicated that sexual size dimorphism was eliminated in stat5b-mutated zebrafish. Herein, the comparative transcriptome analysis was conducted to observe the genes and pathways involved in sexual size dimorphism. The number of male-biased and female-biased genes was much less in the liver of stat5b mutant zebrafish than in wildtype. Gene ontology (GO) enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that multiple pathways related to metabolism were affected upon loss of stat5b function. qRT-PCR results also validated that sexually dimorphic expression of a set of genes was lost when stat5b was mutated. Furthermore, the weighted correlation network analysis (WGCNA) detected many candidate genes related to the growth traits and stat5b function, such as greb1, lepr, and igf2b. Our data suggest that stat5b should regulate the sexually dimorphic gene expression in zebrafish liver and add in understanding of the molecular mechanisms underlying sexual size dimorphism in fish species.

An miR-200 Cluster on Chromosome 23 Regulates Sperm Motility in Zebrafish

Besides its well-documented roles in cell proliferation, apoptosis, and carcinogenesis, the function of the p53-microRNA axis has been recently revealed in the reproductive system. Recent studies indicated that miR-200 family members are dysregulated in nonobstructive azoospermia patients, whereas their functions remain poorly documented. The aim of this study was to investigate the function of the miR-200 family on zebrafish testis development and sperm activity. There was no substantial difference in testis morphology and histology between wild-type (WT) and knockout zebrafish with deletion of miR-200 cluster on chromosome 6 (chr6-miR-200-KO) or on chromosome 23 (chr23-miR-200-KO). Interestingly, compared with WT zebrafish, the chr6-miR-200-KO zebrafish had no difference on sperm motility, whereas chr23-miR-200-KO zebrafish showed significantly improved sperm motility. Consistently, ectopic expression of miR-429a, miR-200a, and miR-200b, which are located in the miR-200 cluster on chromosome 23, significantly reduced motility traits of sperm. Several sperm motility-related genes, such as amh, wt1a, and srd5a2b have been confirmed as direct targets of miR-200s on chr23. 17α-ethynylestradiol (EE2) exposure resulted in upregulated expression of p53 and miR-429a in testis and impairment of sperm motility. Strikingly, in p53 mutant zebrafish testis, the expression levels of miR-200s on chr23 were significantly reduced and accompanied by a stimulation of sperm motility. Moreover, the upregulation of miR-429a associated with EE2 treatment was abolished in testis with p53 mutation. And the impairment of sperm activity by EE2 treatment was also eliminated when p53 was mutated. Together, our results reveal that miR-200 cluster on chromosome 23 controls sperm motility in a p53-dependent manner.

Insights into teleost sex determination from the Seriola dorsalis genome assembly

Background

The assembly and annotation of a genome is a valuable resource for a species, with applications ranging from conservation genomics to gene discovery. Genomic resource development is especially important for species in culture, such as the California Yellowtail (Seriola dorsalis), the likely candidate for the establishment of commercial offshore aquaculture production in southern California. Genomic resource development for this species will improve the understanding of sex and other phenotypic traits, and allow for rapid increases in genetic improvement for and economic gain in culture production.

Results

We describe the assembly and annotation of the S. dorsalis genome, and present resequencing data from 45 male and 45 female wild-caught S. dorsalis used to identify a sex-determining region and marker in this species. The genome assembly captured approximately 93% of the total 685 MB genome with an average coverage depth of 180×. Using the assembled genome, resequencing data from the 90 fish were aligned to place boundaries on the sex-determining region. Sex-specific markers were developed based on a female-specific, 61 nucleotide deletion identified in that region. We hypothesize that Estradiol 17-beta-dehydrogenase is the putative sex-determining gene and propose a plausible genetic mechanism for ZW sex determination in S. dorsalis involving a female-specific deletion of a transcription factor binding motif that may be targeted by Sox3.

Conclusions

Understanding the mechanism of sex determination and development of assays to determine sex is critical both for management of wild fisheries and for development of efficient and sustainable aquaculture practices. In addition, this genome assembly for S. dorsalis will be a substantial resource for a variety of future research applications.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4403-1) contains supplementary material, which is available to authorized users.

Sex-specific markers developed by next-generation sequencing confirmed an XX/XY sex determination system in bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix)

Sex-specific markers are powerful tools for identifying sex-determination system in various animals. Bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix) are two of the most important edible fish in Asia, which have a long juvenility period that can lasts for 4-5 years. In this study, we found one sex-specific marker by next-generation sequencing together with bioinformatics analysis in bighead carp. The male-specific markers were used to perform molecular sexing in the progenies of artificial gynogenetic diploids and found all progenies (n = 160) were females. Meanwhile, around 1 : 1 sex ratio was observed in a total of 579 juvenile offspring from three other families. To further extend the male-specific region, we performed genome walking and got a male-specific sequence of 8,661 bp. Five pairs of primers were designed and could be used to efficiently distinguish males from females in bighead carp and silver carp. The development of these male-specific markers and results of their molecular sexing in different populations provide strong evidence for a sex determination system of female homogametry or male heterogametry (XX/XY) in bighead carp and silver carp. To the best of our knowledge, this is the first report of effective sex-specific markers in these two large carp species.

Potential Contributions of miR-200a/-200b and Their Target Gene-Leptin to the Sexual Size Dimorphism in Yellow Catfish

Sexual size dimorphism is the consequence of differential expression of sex-biased genes related to feeding and growth. Leptin is known to regulate energy balance by regulating food intake. In order to investigate the molecular mechanism of sexual size dimorphism in yellow catfish (Pelteobagrus fulvidraco), the expression of leptin (lep) and its functional receptor (lepr) were detected during larval development. Both lep and lepr have lower expression in males than in females during 1–4 weeks post hatching. 17a-Methyltestosterone (MT) treatment resulted in decreased expression of lep and lepr in both male and female larval fish. Interestingly, the mRNA levels of lep and lepr in juvenile male were significantly decreased compared with juvenile female during short-term fasting periods. Lep was predicted to be a potential target of miR-200a and miR-200b that had an opposite expression pattern to lep in male and female larvas. The results of luciferase reporter assay suggested that lep is a target of miR-200a/-200b. Subsequently, male hormone and fasting treatment have opposite effects on the expression of miR-200a/-200b and lep between males and females. In summary, our results suggest that sexual size dimorphism in fish species is probably caused by the sexually dimorphic expression of leptin, which could be negatively regulated by miR-200a/-200b.

Essential roles of stat5.1/stat5b in controlling fish somatic growth

Signal transducer and activator of transcription 5b (STAT5b) has been identified as a key downstream mediator of growth hormone (GH) signaling in somatic growth of mammalian. However, the corresponding homologue gene of Stat5b is unknown in fish species. In this study, we generated loss-of-function mutants in stat5.1 and stat5.2, two stat5 homologues existing in zebrafish. In stat5.1-deficient zebrafish, a significant reduction of body length and body weight was detected in the embryos/larvae and adults compared with the wild-type control fish, and sexual size dimorphism in adult zebrafish was also eliminated. However, the stat5.2-deficient zebrafish displayed a normal developmental phenotype during all lifespan. Chromatin immunoprecipitation combined with deep sequencing (ChIP-seq) method was adopted to further investigate the potential transcriptional targets of Stat5 protein and cast much light upon the biological function of Stat5. We identified more than 800 genes as transcriptional targets of Stat5 during zebrafish embryogenesis. KEGG analysis indicated that the Stat5 target gene network is predominantly linked to the metabolic pathways, neuroactive ligand-receptor interaction and JAK-STAT signaling pathways. Further validation studies suggested that Stat5.1 protein could directly regulate the expression of gh1, and stat5.1-mutated zebrafish showed a reduction of gh1 mRNA level. In the present study, stat5.1 was revealed as the corresponding homologue gene of Stat5b in fish species. Additionally, we found a novel molecular interaction between Stat5.1/Stat5b and GH, and unraveled a positive feedback loop Stat5.1-GH-Stat5.1 which is necessary for somatic growth and body development in zebrafish.

Numerous mitochondrial DNA haplotypes reveal multiple independent polyploidy origins of hexaploids in Carassius species complex

Evolutionary trajectory and occurrence history of polyploidy have been extensively studied in plants, but they remain quite elusive in vertebrates. Here, we sampled and gathered 4,159 specimens of polyploid Carassius species complex including 1,336 tetraploids and 2,823 hexaploids from a large geographic scale (49 localities) across East Asia, and identified a huge number of 427 diverse haplotypes of mitochondrial control region, in which 74 haplotypes with total occurrence frequency up to 75.498% were shared by hexaploids and tetraploids. Significantly, these diverse haplotypes were clustered into four major lineages, and many haplotypes of hexaploids and tetraploids were intermixed in every lineage. Moreover, the evolutionary trajectory and occurrence history of four different lineages were revealed by a simplified time-calibrated phylogenetic tree, and their geographic distribution frequencies and haplotype diversity were also analyzed. Furthermore, lineage C and D were revealed to undergo population expansion throughout mainland China. Therefore, our current data indicate that hexaploids should undergo multiple independent polyploidy origins from sympatric tetraploids in the polyploid Carassius species complex across East Asia.

Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

Spermatogenesis is a fundamental process in male reproductive biology and depends on precise balance between self-renewal and differentiation of male germ cells. However, the regulative factors for controlling the balance are poorly understood. In this study, we examined the roles of amh and dmrt1 in male germ cell development by generating their mutants with Crispr/Cas9 technology in zebrafish. Amh mutant zebrafish displayed a female-biased sex ratio, and both male and female amh mutants developed hypertrophic gonads due to uncontrolled proliferation and impaired differentiation of germ cells. A large number of proliferating spermatogonium-like cells were observed within testicular lobules of the amh-mutated testes, and they were demonstrated to be both Vasa- and PH3-positive. Moreover, the average number of Sycp3- and Vasa-positive cells in the amh mutants was significantly lower than in wild-type testes, suggesting a severely impaired differentiation of male germ cells. Conversely, all the dmrt1-mutated testes displayed severe testicular developmental defects and gradual loss of all Vasa-positive germ cells by inhibiting their self-renewal and inducing apoptosis. In addition, several germ cell and Sertoli cell marker genes were significantly downregulated, whereas a prominent increase of Insl3-positive Leydig cells was revealed by immunohistochemical analysis in the disorganized dmrt1-mutated testes. Our data suggest that amh might act as a guardian to control the balance between proliferation and differentiation of male germ cells, whereas dmrt1 might be required for the maintenance, self-renewal, and differentiation of male germ cells. Significantly, this study unravels novel functions of amh gene in fish.

Expression profile of amh/Amh during bi-directional sex change in the protogynous orange-spotted grouper Epinephelus coioides

Gonadal differentiation is tightly regulated by the initial sex determining gene and the downstream sex-related genes in vertebrates. However, sex change in fish can alter the sexual fate from one sex to the other. Chemical-induced maleness in the protogynous orange-spotted grouper is transient, and a reversible sex change occurs after the chemical treatment is withdrawn. We used these characteristics to study Amh signaling during bi-directional sex change in the grouper. We successfully induced the female-to-male sex change by chemical (aromatase inhibitor, AI, or methyltestosterone, MT) treatment. A dormant gonad (a low proliferation rate of early germ cells and no characteristics of both sexes) was found during the transient phase of reversible male-to-female sex change after the withdrawal of chemical administration. Our results showed that amh (anti-mullerian hormone) and its receptor amhr2 (anti-mullerian hormone receptor type 2) were significantly increased in the gonads during the process of female-to-male sex change. Amh is expressed in the Sertoli cells surrounding the type A spermatogonia in the female-to-male grouper. Male-related gene (dmrt1 and sox9) expression was immediately decreased in MT-terminated males during the reversible male-to-female sex change. However, Amh expression was found in the surrounding cells of type A spermatogonia-like cells during the transient phase of reversible male-to-female sex change. This phenomenon is correlated with the dormancy of type A spermatogonia-like cells. Thus, Amh signaling is suggested to play roles in regulating male differentiation during the female-to-male sex change and in inhibiting type-A spermatogonia-like cell proliferation/differentiation during the reversible male-to-female sex change. We suggest that Amh signaling might play dual roles during bi-directional sex change in grouper.

Gonadal morphogenesis and sex differentiation in cultured Ussuri catfish Tachysurus ussuriensis

The objective of this study was to investigate the optimal developmental time to perform sex reversal in Ussuri catfish Tachysurus ussuriensis, to develop monosex breeding in aquaculture. Systematic observations of gonadal sex differentiation of P. ussiriensis were conducted. The genital ridge formed at 9 days post fertilization (dpf) and germ cells begin to proliferate at 17 dpf. The ovarian cavity began forming on 21 dpf and completed by 25 dpf while presumptive testis remained quiescent. The primary oocytes were at the chromatin nucleolus stage by 30 dpf, the peri-nucleolus stage by 44 dpf and the cortical alveoli stage by 64 dpf. The germinal vesicle migrated towards the animal pole (polarization) at 120 dpf. In presumptive testis, germ cells entered into mitosis and blood vessels appeared in the proximal gonad on 30 dpf. The efferent duct anlage appeared on 36 dpf and formation of seminal lobules with spermatogonia and lobules interstitium occurred at 120 dpf. Therefore, gonadal sex differentiation occurred earlier in females than in males, with the histological differentiation preceding cytologic differentiation in T. ussuriensis. This indicates that undifferentiated gonads directly differentiate into ovary or testis between 17 and 21 dpf and artificial induction of sexual reversal by oral steroid administration must be conducted before 17 dpf.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3995-9) contains supplementary material, which is available to authorized users.

Distinct herpesvirus resistances and immune responses of three gynogenetic clones of gibel carp revealed by comprehensive transcriptomes

Background

Gibel carp is an important aquaculture species in China, and a herpesvirus, called as Carassius auratus herpesvirus (CaHV), has hampered the aquaculture development. Diverse gynogenetic clones of gibel carp have been identified or created, and some of them have been used as aquaculture varieties, but their resistances to herpesvirus and the underlying mechanism remain unknown.

Results

To reveal their susceptibility differences, we firstly performed herpesvirus challenge experiments in three gynogenetic clones of gibel carp, including the leading variety clone A⁺, candidate variety clone F and wild clone H. Three clones showed distinct resistances to CaHV. Moreover, 8772, 8679 and 10,982 differentially expressed unigenes (DEUs) were identified from comparative transcriptomes between diseased individuals and control individuals of clone A⁺, F and H, respectively. Comprehensive analysis of the shared DEUs in all three clones displayed common defense pathways to the herpesvirus infection, activating IFN system and suppressing complements. KEGG pathway analysis of specifically changed DEUs in respective clones revealed distinct immune responses to the herpesvirus infection. The DEU numbers identified from clone H in KEGG immune-related pathways, such as “chemokine signaling pathway”, “Toll-like receptor signaling pathway” and others, were remarkably much more than those from clone A⁺ and F. Several IFN-related genes, including Mx1, viperin, PKR and others, showed higher increases in the resistant clone H than that in the others. IFNphi3, IFI44-like and Gig2 displayed the highest expression in clone F and IRF1 uniquely increased in susceptible clone A⁺. In contrast to strong immune defense in resistant clone H, susceptible clone A⁺ showed remarkable up-regulation of genes related to apoptosis or death, indicating that clone A⁺ failed to resist virus offensive and evidently induced apoptosis or death.

Conclusions

Our study is the first attempt to screen distinct resistances and immune responses of three gynogenetic gibel carp clones to herpesvirus infection by comprehensive transcriptomes. These differential DEUs, immune-related pathways and IFN system genes identified from susceptible and resistant clones will be beneficial to marker-assisted selection (MAS) breeding or molecular module-based resistance breeding in gibel carp.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3945-6) contains supplementary material, which is available to authorized users.

Wider geographic distribution and higher diversity of hexaploids than tetraploids in Carassius species complex reveal recurrent polyploidy effects on adaptive evolution

Polyploidy roles on adaptive evolution and ecological novelty have been extensively studied in plants but remained unclear in vertebrates owing to the rare polyploidy incidences. Here, a huge number of 3105 specimens in Carassius species complex including 2211 hexaploids and 894 tetraploids were sampled from 34 locations through mainland China. And hexaploids had wider geographic distribution than tetraploids especially in the areas with high altitude, high latitude and low annual precipitation. Then, an approximate 1050 bp transferrin (tf) fragments were amplified from all the samples, and 526 tf alleles were identified from a total of 37260 sequences at last. Intriguingly, higher nucleotide diversity of tf alleles in hexaploids than in tetraploids was revealed. Moreover, via phylogenetic analysis of tf alleles, potential origin center of Carassius species complex was deduced to be Yangtze River basin and hexaploids should undergo multiple independent polyploidy origins from sympatric tetraploids. These findings indicate that the hexaploids might possess stronger environmental adaptation and ecological novelty than the tetraploids, which provide an association paradigm of recurrent polyploidy and ecological context in polyploid vertebrates.

Seasonal and daily protandry in a cyprinid fish

In polygynandrous mating systems, in which females limit reproductive success, males can increase their success by investing in courtship. Earlier arrival at the spawning ground compared to when females arrive may increase their opportunities in competitive mating systems. In this study, we used passive telemetry to test whether a male minnow known as the asp, Leuciscus aspius, times its arrival at spawning grounds relative to the arrival of females. Males arrived in a model stream approximately five days earlier than females on average and left four to five days later than females over two years. Both sexes performed a daily migration between a staging ground (standing water, low energy costs) and the fluvial spawning ground (high energy costs). Fish abundance peaked twice a day, with a major peak at sunset and a minor peak at sunrise and with the evening peak abundance for males occurring 1 hour 40 minutes earlier than that of females. The number of females on the spawning ground never exceeded the number of males. While the degree of protandry is hypothesized to be influenced by the operational sex ratio (ranging from 0.5 to 1 in our study), our data did not support this theory.

Dmrt1 is required for primary male sexual differentiation in Chinese soft-shelled turtle Pelodiscus sinensis

In vertebrates, the primary sex-determining signals that initiate sexual development are remarkably diverse, ranging from complete genetic to environmental cues. However, no sex determination-related genes have been functionally identified in reptiles. Here, we characterized a conserved DM domain gene, Dmrt1, in Chinese soft-shelled turtle Pelodiscus sinensis (P. sinensis), which exhibits ZZ/ZW sex chromosomes. Dmrt1 exhibited early male-specific embryonic expression, preceding the onset of gonadal sex differentiation. The expression of Dmrt1 was induced in ZW embryonic gonads that were masculinized by aromatase inhibitor treatment. Dmrt1 knockdown in ZZ embryos by RNA interference resulted in male to female sex reversal, characterized by obvious feminization of gonads, significant down-regulation of testicular markers Amh and Sox9, and remarkable up-regulation of ovarian regulators, Cyp19a1 and Foxl2. Conversely, ectopic expression of Dmrt1 led to largely masculinized genetic females, production of Amh and Sox9, and a decline in Cyp19a1 and Foxl2. These findings demonstrate that Dmrt1 is both necessary and sufficient to initiate testicular development, thereby acting as an upstream regulator of the male pathway in P. sinensis.

Expression pattern and functional analysis of fundc1 in rare minnow (Gobiocypris rarus)

Fundc1 is a mitochondrial outer membrane protein and plays important roles in mitochondria fission and hypoxia-induced mitophagy in mammalian cells. However, there is no relevant report of fundc1 in fish. In the present study, we cloned a 942bp fundc1 cDNA from rare minnow. The cDNA, designated as Grfundc1 cDNA, contains an open reading frame (ORF) of 459bp which encodes a polypeptide of 152 amino acid residues. Comparisons of deduced amino acid sequences demonstrated that Grfundc1 was highly homologous with those of other vertebrates. RT-PCR and real time PCR detection revealed that the transcripts of Grfundc1 were not detectable in the unfertilized eggs and had high levels at blastula and gastrula stages. Whole mount in situ hybridization analysis observed that Grfundc1 was ubiquitously expressed at early stage and later riched in specific regions, such as brain, branchial arch, eye and somite during embryogenesis. Grfundc1 was expressed in all the tissues of rare minnow adult, including brain, liver, gill, eyes, heart, kidney, intestine, muscle, testis and ovary. The expression of Grfundc1 in the brain, gill, heart and eye of rare minnow adult was significantly down-regulated by hypoxia. Similar hypoxic response was observed in the rare minnow embryos at 48hpf following hypoxia exposure. Functional analysis showed that knockdown of Grfundc1 significantly caused defects in the body axis and dorsal neural tissues of rare minnow embryos. These results indicate that Grfundc1 may play important roles in embryogenesis in fish.

MicroRNA-203a regulates fast muscle differentiation by targeting dmrt2a in zebrafish embryos

Dmrt2b (doublesex and mab-3 related transcription factor 2b) has been revealed to be involved in zebrafish slow muscle development. However, the function of dmrt2a, a paralogue gene of dmrt2b, remains unclear during zebrafish muscle development. Here, we demonstrated that knockdown of dmrt2a resulted in severe developmental defects, and caused downregulation of fast muscle marker myhz-2 and upregulation of slow muscle marker myhz-5, respectively. It is known that microRNAs (miRNAs) control many biological events including muscle development. Dmrt2a was predicted to be a target gene of miR-203, which was further verified by luciferase reporter assay, since miR-203a was found to directly reduce the expression of dmrt2a by binding to the seed sequence of its 3'UTR. After miR-203a injection into zebrafish embryos, the expression of dmrt2a was significantly inhibited. Similar to the effect of dmrt2a knockdown, miR-203a overexpression led to downregulation of myhz-2 and upregulation of myhz-5. Our studies indicated that miR-203a directly regulated dmrt2a expression to control fast and slow muscle differentiation, while overexpression of miR-203a or knockdown of dmrt2a will impair fast muscle development and promote slow muscle development.

A novel male-specific SET domain-containing gene setdm identified from extra microchromosomes of gibel carp males

Various genes have been screened on extra chromosomes, but their molecular characterization, expression pattern and biological function are still unclear. Here, we utilized a male-specific sequence of polyploid gibel carp (Carassius gibelio) to identify a novel male-specific SET (Su(var)3-9, Enhancer-of-zeste, Trithorax) domain-containing gene setdm on extra microchromosomes of gibel carp males. And setdm was characterized in molecule and expression aspects, in which its expression was specific to testis, and had relative high transcription during middle/late stages of testis development. Moreover, prominent expression of Setdm in spermatogenic cells was observed in testis through immunofluorescence co-localization analysis. These results suggest that biological function of setdm might be related to testis development and spermatogenesis of gibel carp. Additionally, the homeologous gene setdmf of setdm, was also characterized, and its expression was gonad-specific, in which its expressed product was detected to mainly distribute in gametogenic cells of testis and ovary, and to have dynamic expression pattern similar to that of setdm. Based on the current results, we propose that the novel male-specific setdm on extra microchromosomes might be functional divergence gene of the gonad-specific setdmf. Therefore, these findings will help us to further understand evolutionary fate and functional role of genes on extra microchromosomes.

Identification and profiling of Cyprinus carpio microRNAs during ovary differentiation by deep sequencing

Background

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that regulate gene expression by targeting specific mRNAs. However, the possible role of miRNAs in the ovary differentiation and development of fish is not well understood. In this study, we examined the expression profiles and differential expression of miRNAs during three key stages of ovarian development and different developmental stages in common carp Cyprinus carpio.

Results

A total of 8765 miRNAs were identified, including 2155 conserved miRNAs highly conserved among various species, 145 miRNAs registered in miRBase for common carp, and 6505 novel miRNAs identified in common carp for the first time. Comparison of miRNA expression profiles among the five libraries identified 714 co-expressed and 2382 specific expressed miRNAs. Overall, 150, 628, and 431 specifically expressed miRNAs were identified in primordial gonad, juvenile ovary, and adult ovary, respectively. MiR-6758-3p, miR-3050-5p, and miR-2985-3p were highly expressed in primordial gonad, miR-3544-5p, miR-6877-3p, and miR-9086-5p were highly expressed in juvenile ovary, and miR-154-3p, miR-5307-5p, and miR-3958-3p were highly expressed in adult ovary. Predicted target genes of specific miRNAs in primordial gonad were involved in many reproductive biology signaling pathways, including transforming growth factor-β, Wnt, oocyte meiosis, mitogen-activated protein kinase, Notch, p53, and gonadotropin-releasing hormone pathways. Target-gene prediction revealed upward trends in miRNAs targeting male-bias genes, including dmrt1, atm, gsdf, and sox9, and downward trends in miRNAs targeting female-bias genes including foxl2, smad3, and smad4. Other sex-related genes such as sf1 were also predicted to be miRNA target genes.

Conclusions

This comprehensive miRNA transcriptome analysis demonstrated differential expression profiles of miRNAs during ovary development in common carp. These results could facilitate future exploitation of the sex-regulatory roles and mechanisms of miRNAs, especially in primordial gonads, while the specifically expressed miRNAs represent candidates for studying the mechanisms of ovary determination in Yellow River carp.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3701-y) contains supplementary material, which is available to authorized users.

Potential contributions of heat shock proteins and related genes in sexual differentiation in yellow catfish (Pelteobagrus fulvidraco)

Sex determination and differentiation in ectotherms are very complicated affairs and usually affected by both genetic and environmental factors. Because of their temperature-sensitive expression, heat shock proteins (HSPs) are good candidates for temperature-dependent sex determination (TSD). Similar to most thermosensitive fish species, the male to female ratio increases with temperature in yellow catfish (Pelteobagrus fulvidraco). Yellow catfish is also a type of sexual size dimorphic fish, and the male individuals grow much faster than females of the same age. Therefore, research of sex differentiation in yellow catfish is important in aquiculture. In this attempt, a total of seven HSPs and related genes were identified from transcriptomes of yellow catfish by 454 pyrosequencing and Solexa sequencing that we did previously, including five genes with complete open reading frame (ORF). Phylogenetically, all these genes were compared with their counterparts from other vertebrates. All these genes were sex-biased expressed in gonads. Hspa5, Hip, and Cdc37 were expressed more highly in ovary than in testis, whereas Hsp90α, Hspb2, Hspb8, and Hspbp1 were expressed more highly in testis than in ovary. Additionally, the expression of these genes was assessed after 17α-methyltestosterone (MT) and 17α-ethinylestradiol (EE2) treatment, respectively. Our result showed that working as co-chaperones, these HSPs and related genes may regulate sex steroid receptor activities to influence gonad development in yellow catfish. Our work would help in the understanding of the mechanism of sexual differentiation in teleosts.

Sequential, Divergent, and Cooperative Requirements of Foxl2a and Foxl2b in Ovary Development and Maintenance of Zebrafish

Foxl2 is essential for mammalian ovary maintenance. Although sexually dimorphic expression of foxl2 was observed in many teleosts, its role and regulative mechanism in fish remained largely unclear. In this study, we first identified two transcript variants of foxl2a and its homologous gene foxl2b in zebrafish, and revealed their specific expression in follicular layer cells in a sequential and divergent fashion during ovary differentiation, maturation, and maintenance. Then, homozygous foxl2a mutants (foxl2a-/-) and foxl2b mutants (foxl2b-/-) were constructed and detailed comparisons, such as sex ratio, gonadal histological structure, transcriptome profiling, and dynamic expression of gonadal development-related genes, were carried out. Initial ovarian differentiation and oocyte development occur normally both in foxl2a-/- and foxl2b-/- mutants, but foxl2a and foxl2b disruptions result in premature ovarian failure and partial sex reversal, respectively, in adult females. In foxl2a-/- female mutants, sox9a-amh/cyp19a1a signaling was upregulated at 150 days postfertilization (dpf) and subsequently oocyte apoptosis was triggered after 180 dpf. In contrast, dmrt1 expression was greater at 105 dpf and increased several 100-fold in foxl2b-/- mutated ovaries at 270 dpf, along with other testis-related genes. Finally, homozygous foxl2a-/-/foxl2b-/- double mutants were constructed in which complete sex reversal occurs early and testis-differentiation genes robustly increase at 60 dpf. Given mutual compensation between foxl2a and foxl2b in foxl2b-/- and foxl2a-/- mutants, we proposed a model in which foxl2a and foxl2b cooperate to regulate zebrafish ovary development and maintenance, with foxl2b potentially having a dominant role in preventing the ovary from differentiating as testis, as compared to foxl2a.

Divergent Expression Patterns and Function Implications of Four nanos Genes in a Hermaphroditic Fish, Epinephelus coioides

Multiple nanos genes have been characterized in several fishes, but the functional implications of their various expression patterns remain unclear. In this study, we identified and characterized four nanos genes from a hermaphroditic fish orange-spotted grouper, Epinephelus coioides. Ecnanos1a and Ecnanos1b show divergent expression patterns, and the dynamic expression change of Ecnanos1a in pituitaries during sex change is associated with testis differentiation and spermatogenesis. Ecnanos2 and Ecnanos3 might be germline stem cells (GSCs) and primordial germ cells (PGCs)-specific markers, respectively. Significantly, Ecnanos3 3'-untranslated region (UTR) is necessary for PGC specific expression, where a non-canonical "GCACGTTT" sequence is required for miR-430-mediated repression of Ecnanos3 RNA. Furthermore, grouper Dead end (Dnd) can relieve miR-430 repression in PGCs by associating with a 23 bp U-rich region (URR) in Ecnanos3 3'-UTR. The current study revealed the functional association of multiple nanos genes with PGC formation and germ cell development in orange-spotted grouper, and opened up new possibilities for developing biotechnologies through utilizing the associations between Ecnanos3 and PGCs or between Ecnanos2 and GSCs in the hermaphroditic fish.

Identification and characterization of a Masculinizer (Masc) gene involved in sex differentiation in Artemia

The sex of relatively primitive animals such as invertebrates is mostly determined by environmental factors and chromosome ploidy. Heteromorphic chromosomes may also play an important role, as in the ZW system in lepidopterans. However, the mechanisms of these various sex determination systems are still largely undefined. In the present study, a Masculinizer gene (Ar-Masc) was identified in the crustacean Artemia franciscana Kellogg 1906. Sequence analysis revealed that the 1140-bp full-length open reading frame of Ar-Masc encodes a 380-aa protein containing two CCCH-type zinc finger domains having a high degree of shared identities with the MASC protein characterized in the silkworm Bombyx mori, which has been determined to participate in the production of male-specific splice variants. Furthermore, although Ar-Masc could be detected in almost all stages in both sexual and parthenogenetic Artemia, there were significant variations in expression between these two reproductive modes. Firstly, qRT-PCR and Western blot analysis showed that levels of both Ar-Masc mRNA and protein in sexual nauplii were much higher than in parthenogenetic nauplii throughout the hatching process. Secondly, both sexual and parthenogenetic Artemia had decreased levels of Ar-Masc along with the embryonic developmental stages, while the sexual ones had a relatively higher and more stable expression than those of parthenogenetic ones. Thirdly, immunofluorescence analysis determined that sexual individuals had higher levels of Ar-MASC protein than parthenogenetic individuals during embryonic development. Lastly, RNA interference with dsRNA showed that gene silencing of Ar-Masc in sexual A. franciscana caused the female-male ratio of progeny to be 2.19:1. These data suggest that Ar-Masc participates in the process of sex determination in A. franciscana, and provide insight into the evolution of sex determination in sexual organisms.

Loss of stat3 function leads to spine malformation and immune disorder in zebrafish

STAT (Signal Transducers and Activators of Transcription) gene family members have been revealed to be involved in cell growth and differentiation in vertebrates. Despite their physiological importance, their functions are poorly studied at organ and systemic levels. In this study, we performed a genome-wide analysis using data from invertebrates to vertebrates to identify STAT genes and analyze their evolutionary history. Interestingly, the STAT gene family undergoes genome duplications during the evolutionary history with STAT3 homologues firstly appearing in the basal extant vertebrate, sea lamprey, suggesting its possible roles in spine formation. To investigate the functions of stat3 in fish species, TALEN technology was performed to generate mutant zebrafish lines. Stat3 mutant zebrafish showed no obvious defects at early developmental stage but displayed severe lateral and vertical curvature of the spine (scoliosis), spine fracture and the incomplete bone joints with narrower junction between vertebrae at early juvenile stage, as indicated by Alizarin red and Alcian blue staining, radiography and micro-computed tomography (MicroCT) analysis. Transcriptome analysis reveals dramatic alterations in a number of genes involved in immune and infection response, skeletal development and somatic growth, especially downregulated expression of collagen gene family, in the juvenile stat3 mutant zebrafish. Moreover, most of the collagen genes were detected to have abnormal expression pattern during the formation of spine deformities in stat3 mutants. Our data reveal that stat3 is specially expressed in vertebrates and required for normal spine development and immune function in zebrafish.