Transcriptomic analyses reveal novel genes with sexually dimorphic expression in the zebrafish gonad and brain

Zebrafish gender-specific anxiety-like behavioral and physiological reactions elicited by caffeine

Caffeine exerts a biphasic effect on zebrafish behavior. High doses of caffeine have been associated with increased stress and anxiety, whereas low doses have been found to enhance performance on tasks requiring focus and attention. However, the sex-specific nature of these biphasic effects on behavior and physiology remains unclear. This study assessed the behavioral responses and hormone levels in male and female zebrafish after acute exposure to caffeine ranging from 0.3 to 600mg/L. The results showed no significant difference in caffeine intake between males and females after acute exposure at each concentration. Caffeine-induced behavioral and physiological responses indicated a threshold dosage existed between 30 and 300mg/L. Female fish displayed increased anxiety-like behavioral phenotypes, i.e., latency to upper and freezing, whereas males exhibited more erratic movement following acute exposure to a high-dose treatment. In addition, females exhibited a significant increase in whole-body cortisol levels, while males experienced a testosterone elevation at 300mg/L of caffeine acute exposure. There was a significant decrease in the duration of erratic movements in males treated with the androgen receptor antagonist flutamide compared to the control group. The transcriptome analysis uncovered 511 and 592 up-regulated and 761 and 922 down-regulated differential expression genes in males and females, respectively, compared to the control. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analysis revealed that caffeine has the potential to impact various pathways in zebrafish, including phototransduction and steroid hormone biosynthesis. Our findings demonstrate that testosterone and cortisol play a combined role in regulating stress responses in both behavior and physiology. Furthermore, our study highlights the significance of encompassing both male and female zebrafish as a model system.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

Comparative transcriptome analysis of ovaries and testes reveals sex-biased genes and pathways in zebrafish

Zebrafish (Danio rerio) is a widely recognized and extensively studied model organism in scientific research. The regulatory mechanism of gonadal development and differentiation of this species has aroused considerable attention. Nonetheless, the major sex-biased genes and pathways associated with gonadal development remain elusive. Therefore, to comprehend this intricate process, gonadal transcriptome sequencing was carried out to identify differentially expressed genes (DEGs) between the testes and ovaries of adult zebrafish. The preliminary assessment yielded a total of 23,529,272 and 23,521,368 clean reads from the cDNA libraries of ovaries and testes. Afterward, a comparative analysis of the transcriptome revealed 3,604 upregulated and 11,371 downregulated DEGs in the ovaries compared to the testes. Of these genes, 428 were exclusively expressed in females, while 3,516 were exclusively expressed in males. Additionally, further assessments were conducted to explore the functions associated with these DEGs in various biological processes. The data revealed their involvement in sex-biased pathways, such as progesterone-mediated oocyte maturation, oocyte meiosis, cytokine-cytokine receptor interaction, and cardiac muscle contraction. Finally, the expression levels of 14 sex-biased DEGs (cdc20, ccnb1, ypel3, chn1, bmp15, rspo1, tnfsf10, egfra, acta2, cox8a, gsdf, dmrt1, star, and cyp17a1) associated with the enriched pathways were subjected to further validation through qRT-PCR. The data acquired from these investigations offer valuable resources to support further exploration of the mechanisms governing sexual dimorphism and gonadal development in zebrafish.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Development of epigenetic biomarkers for the identification of sex and thermal stress in fish using DNA methylation analysis and machine learning procedures

The sex ratio is a key ecological demographic parameter crucial for population viability. However, the epigenetic mechanisms operating during gonadal development regulating gene expression and the sex ratio remain poorly understood. Moreover, there is interest in the development of epigenetic markers associated with a particular phenotype or as sentinels of environmental effects. Here, we profiled DNA methylation and gene expression of 10 key genes related to sex development and stress, including steroidogenic enzymes, and growth and transcription factors. We provide novel information on the sex-related differences and on the influence of elevated temperature on these genes in zebrafish, a species with mixed genetic and environmental influences on sex ratios. We identified both positive (e.g., amh, cyp11c and hsd11b2) and negative (e.g., cyp11a1 and dmrt1) correlations in unexposed males, and negative correlation (amh) in exposed females between DNA methylation and gene expression levels. Further, we combined DNA methylation analysis with machine learning procedures and found a series of informative CpGs capable not only of correctly identifying sex (based on cyp19a1a DNA methylation levels) but also of identifying whether males and females had been exposed to abnormally elevated temperature when young (based on amh and foxl2a DNA methylation levels, respectively). This was achieved in the absence of conspicuous morphological alterations of the gonads. These DNA methylation-based epigenetic biomarkers represent molecular resources that can correctly recapitulate past thermal history and pave the way for similar findings in other species to assess potential ecological effects of environmental disturbances in the context of climate change.

Female zebrafish (Danio rerio) demonstrate stronger preference for established shoals over newly-formed shoals in the three-tank open-swim preference test

Zebrafish (Danio rerio) share a considerable amount of biological similarity with mammals, including identical or homologous gene expression pathways, neurotransmitters, hormones, and cellular receptors. Zebrafish also display complex social behaviors like shoaling and schooling, making them an attractive model for investigating normal social behavior as well as exploring impaired social function conditions such as autism spectrum disorders. Newly-formed and established shoals exhibit distinct behavior patterns and inter-member interactions that can convey the group’s social stability. We used a three-chamber open-swim preference test to determine whether individual zebrafish show a preference for an established shoal over a newly-formed shoal. Results indicated that both sexes maintained greater proximity to arena zones nearest to the established shoal stimulus. In addition, we report the novel application of Shannon entropy to discover sex differences in systematicity of responses not revealed by unit-based measurements; male subjects spent more time investigating between the two shoals than female subjects. This novel technique using established versus newly-formed shoals can be used in future studies testing transgenics and pharmacological treatments that mimic autism spectrum disorder and other disorders that affect social interaction.

Sex-specific DNA methylation: impact on human health and development

Human evolution has shaped gender differences between males and females. Over the years, scientific studies have proposed that epigenetic modifications significantly influence sex-specific differences. The evolution of sex chromosomes with epigenetics as the driving force may have led to one sex being more adaptable than the other when exposed to various factors over time. Identifying and understanding sex-specific differences, particularly in DNA methylation, will help determine how each gender responds to factors, such as disease susceptibility, environmental exposure, brain development and neurodegeneration. From a medicine and health standpoint, sex-specific methylation studies have shed light on human disease severity, progression, and response to therapeutic intervention. Interesting findings in gender incongruent individuals highlight the role of genetic makeup in influencing DNA methylation differences. Sex-specific DNA methylation studies will empower the biotechnology and pharmaceutical industry with more knowledge to identify biomarkers, design and develop sex bias drugs leading to better treatment in men and women based on their response to different diseases.

Integrated analysis of the miRNA and mRNA expression profiles in Leiocassis longirostris at gonadal maturation

Leiocassis longirostris is a commercially important fish species that shows a sexually dimorphic growth pattern. A lack of molecular data from the gonads of this species has hindered research and selective breeding efforts. In this study, we conducted a comprehensive analysis of the expression profile of miRNA and mRNA to explore their regulatory roles in the gonadal maturation stage of L. longirostris. We identified 60 differentially expressed miRNAs and 20,752 differentially expressed genes by sequencing. A total of 90 miRNAs and 21 target genes involved in gonad development and sex determination were identified. Overall, the results of this study enhance our understanding of the molecular mechanisms underlying sex determination and differentiation and provide valuable genomic information for the selective breeding of L. longirostris.

Genome biology of the Darkedged Splitfin, Girardinichthys multiradiatus, and the evolution of sex chromosomes and placentation

Viviparity evolved independently about 150 times in vertebrates and more than 20 times in fish. Several lineages added to the protection of the embryo inside the body of the mother, the provisioning of nutrients, and physiological exchange. This often led to the evolution of a placenta. Among fish, one of the most complex systems serving the function of the placenta is the embryonal trophotaenia/ovarian luminal epithelium of the goodeid fishes. For a better understanding of this feature and others of this group of fishes, high-quality genomic resources are essential. We have sequenced the genome of the darkedged splitfin, Girardinichthys multiradiatus. The assembly is chromosome level and includes the X and Y Chromosomes. A large male-specific region on the Y was identified covering 80% of Chromosome 20, allowing some first inferences on the recent origin and a candidate male sex determining gene. Genome-wide transcriptomics uncovered sex-specific differences in brain gene expression with an enrichment for neurosteroidogenesis and testis genes in males. The expression signatures of the splitfin embryonal and maternal placenta showed overlap with homologous tissues including human placenta, the ovarian follicle epithelium of matrotrophic poeciliid fish species and the brood pouch epithelium of the seahorse. Our comparative analyses on the evolution of embryonal and maternal placenta indicate that the evolutionary novelty of maternal provisioning development repeatedly made use of genes that already had the same function in other tissues. In this way, preexisting modules are assembled and repurposed to provide the molecular changes for this novel trait.

Profile of gene expression changes during estrodiol-17β-induced feminization in the Takifugu rubripes brain

Background

As the critical tissue of the central nervous system, the brain has been found to be involved in gonad development. Previous studies have suggested that gonadal fate may be affected by the brain. Identifying brain-specific molecular changes that occur during estrodiol-17β (E₂) -induced feminization is crucial to our understanding of the molecular control of sex differentiation by the brains of fish.

Results

In this study, the differential transcriptomic responses of the Takifugu rubripes larvae brain were compared after E₂ treatment for 55 days. Our results showed that 514 genes were differentially expressed between E₂-treated-XX (E-XX) and Control-XX (C-XX) T. rubripes, while 362 genes were differentially expressed between E₂-treated-XY (E-XY) and Control-XY (C-XY). For example, the expression of cyp19a1b, gnrh1 and pgr was significantly up-regulated, while st, sl, tshβ, prl and pit-1, which belong to the growth hormone/prolactin family, were significantly down-regulated after E₂ treatment, in both sexes. The arntl1, bhlbe, nr1d2, per1b, per3, cry1, cipc and ciart genes, which are involved in the circadian rhythm, were also found to be altered. Differentially expressed genes (DEGs), which were identified between E-XX and C-XX, were significantly enriched in neuroactive ligand-receptor interaction, arachidonic acid metabolism, cytokine-cytokine receptor interaction and the calcium signaling pathway. The DEGs that were identified between E-XY and C-XY were significantly enriched in tyrosine metabolism, phenylalanine metabolism, arachidonic acid metabolism and linoleic acid metabolism.

Conclusion

A number of genes and pathways were identified in the brain of E₂-treated T. rubripes larvae by RNA-seq. It provided the opportunity for further study on the possible involvement of networks in the brain-pituitary-gonadal axis in sex differentiation in T. rubripes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08158-0.

Advantages of omics technology for evaluating cadmium toxicity in zebrafish

In the last decade, several advancements have been made in omics technologies and they have been applied extensively in diverse research areas. Especially in toxicological research, omics technology can efficiently and accurately generate relevant data on the molecular dynamics associated with adverse outcomes. Toxicomics is defined as the combination of toxicology and omics technologies and encompasses toxicogenomics, toxicoproteomics, and toxicometabolomics. This paper reviews the trend of applying omics technologies to evaluate cadmium (Cd) toxicity in zebrafish (D. rerio). Cd is a toxic heavy metal posing several environmental concerns; however, it is being used widely in everyday life. Zebrafish embryos and larvae are employed as standard models for many toxicity tests because they share 71.4% genetic homology with humans. This study summarizes the toxicity of Cd on the nerves, liver, heart, skeleton, etc. of zebrafish and introduces detailed omics techniques to understand the results of the toxicomic studies. Finally, the trend of toxicity evaluation in the zebrafish model of Cd based on omics technology is presented.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

Analyses of mRNA-seq and miRNA-seq of the brain reveal the sex differences of gene expression and regulation before and during gonadal differentiation in 17β-estradiol or 17α-methyltestosterone-induced olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is a commercially important flatfish species cultured in East Asia. Female flounders generally grow more rapidly than males, therefore control of the sex ratio seems to be a proposed way to increase production. However, the sex determination gene and sex determination mechanism have yet been elucidated. The brain is an important organ that is involved in gonadal development. To explore the sex differences of gene expression in the brain before and during the flounder gonadal differentiation, we used messenger RNA (mRNA)-seq technology to investigate transcriptomes of male and female brains. Between female and male brains, 103 genes were differentially expressed before ovarian differentiation, 16 genes were differentially expressed before testicular differentiation, and 64 genes were differentially expressed during gonadal differentiation. According to annotation and Kyoto Encyclopedia of Genes and Genomes information, the differentially expressed genes (DEGs) were involved in circadian rhythm, circadian rhythm-fly, circadian entrainment, dopaminergic synapse, calcium signaling, glutamatergic synapse, taste transduction, herpes simplex infection, long-term depression, retrograde endocannabinoid signaling, and the synaptic vesicle cycle pathways. MicroRNA (miRNA)-seq was performed during the gonadal differentiation and the target genes of miRNAs were predicted. Integrated analysis of mRNA-seq and miRNA-seq showed that 29 of the 64 DEGs were regulated by the differentially expressed miRNAs during the gonadal differentiation. Our study provides a basis for further studies of brain sex differentiation and the molecular mechanism of sex determination in olive flounder.

Heat Shock Factor 5 Is Essential for Spermatogenesis in Zebrafish

Heat shock factors (Hsfs) are transcription factors that regulate responses to heat shock and other environmental stimuli. Four heat shock factors (Hsf1-4) have been characterized from vertebrates to date. In addition to stress response, they also play important roles in development and gametogenesis. Here, we study the fifth member of heat shock factor family, Hsf5, using zebrafish as a model organism. Mutant hsf5^-/- males, generated by CRISPR/Cas9 technique, were infertile with drastically reduced sperm count, increased sperm head size, and abnormal tail architecture, whereas females remained fertile. We show that Hsf5 is required for progression through meiotic prophase 1 during spermatogenesis as suggested by the accumulation of cells in the leptotene and zygotene-pachytene stages and increased apoptosis in post-meiotic cells. hsf5^-/- mutants show gonadal misregulation of a substantial number of genes with roles in cell cycle, apoptosis, protein modifications, and signal transduction, indicating an important role of Hsf5 in early stages of spermatogenesis.

Transcriptome profiling analysis of sex-based differentially expressed mRNAs and lncRNAs in the brains of mature zebrafish (Danio rerio)

Background

Similar to humans, the zebrafish brain plays a central role in regulating sexual reproduction, maturation and sexual behavior. However, systematic studies of the dimorphic patterns of gene expression in the brain of male and female zebrafish are lacking.

Results

In this study, the mRNA and lncRNA expression profiles were obtained from the brain tissue samples of the three male and three female zebrafish by high-throughput transcriptome sequencing. We identified a total of 108 mRNAs and 50 lncRNAs with sex-based differential expression. We randomly selected four differentially expressed genes for RT-qPCR verification and the results certified that the expression pattern showed a similar trend between RNA-seq and RT-qPCR results. Protein-protein interaction network analysis, Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to obtain the biological significance of differentially expressed mRNA in the brain dimorphism of zebrafish. Finally, a Pearson correlation analysis was performed to construct the co-expression network of the mRNAs and lncRNAs.

Conclusions

We found that 12 new lncRNAs not only have significant gender specificity in the brain of zebrafish, and this finding may provide a clue to further study of the functional difference between male and female zebrafish brain.

Identification of sex differences in zebrafish (Danio rerio) brains during early sexual differentiation and masculinization using 17α-methyltestoterone

Sexual behavior in teleost fish is highly plastic. It can be attributed to the relatively few sex differences found in adult brain transcriptomes. Environmental and hormonal factors can influence sex-specific behavior. Androgen treatment stimulates behavioral masculinization. Sex dimorphic gene expression in developing teleost brains and the molecular basis for androgen-induced behavioral masculinization are poorly understood. In this study, juvenile zebrafish (Danio rerio) were treated with 100 ng/L of 17 alpha-methyltestosterone (MT) during sexual development from 20 days post fertilization to 40 days and 60 days post fertilization. We compared brain gene expression patterns in MT-treated zebrafish with control males and females using RNA-Seq to shed light on the dynamic changes in brain gene expression during sexual development and how androgens affect brain gene expression leading to behavior masculinization. We found modest differences in gene expression between juvenile male and female zebrafish brains. Brain aromatase (cyp19a1b), prostaglandin 3a synthase (ptges3a), and prostaglandin reductase 1 (ptgr1) were among the genes with sexually dimorphic expression patterns. MT treatment significantly altered gene expression relative to both male and female brains. Fewer differences were found among MT-treated brains and male brains compared to female brains, particularly at 60 dpf. MT treatment upregulated the expression of hydroxysteroid 11-beta dehydrogenase 2 (hsd11b2), deiodinase, iodothyronine, type II (dio2), and gonadotrophin releasing hormones (GnRH) 2 and 3 (gnrh2 and gnrh3) suggesting local synthesis of 11-ketotestosterone, triiodothyronine, and GnRHs in zebrafish brains which are influenced by androgens. Androgen, estrogen, prostaglandin, thyroid hormone, and GnRH signaling pathways likely interact to modulate teleost sexual behavior.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

Spatiotemporal expression profile of embryonic and adult ankyrin repeat and EF-hand domain containing protein 1-encoding genes ankef1a and ankef1b in zebrafish

Recent human next-generation sequencing (NGS) studies indicate a correlation between ANKEF1 (ankyrin repeat and EF-hand domain containing protein 1) expression and cilia formation or function. Additionally, a single study conducted in the African clawed frog (Xenopus laevis) showed ankef1 is down-regulated after pharmacological fibroblast growth factor (FGF) inhibition and plays a role in protocadherin-mediated cell protrusion and adhesion. That study also revealed a critical role for ankef1 in the embryonic development of the frog, with morphants exhibiting phenotypes including spina bifida and a shortened body axis. Interestingly, while little is known about ANKEF1 function in other vertebrate systems, recent proteomic data has shown ANKEF1 enriched in ciliated cells. Likewise, publicly available EST profile databases imply ANKEF1 expression in multiple human tissues, including high levels in the testes. Together, these previous studies suggest an important role for ANKEF1 in ciliated tissues and during embryonic development. Here, we report cloning of zebrafish (Danio rerio) ankef1a, as well as its paralog, ankef1b, and expression analyses by whole-mount in situ hybridization (WISH) and quantitative polymerase chain reaction (qPCR) during embryonic development and in adult tissues. WISH shows both forms are ubiquitously expressed early in development, with more discrete expression of both transcripts in embryonic tissues known to precede or possess motile cilia, including dorsal forerunner cells (DFC) and the otic vesicles, respectively. Additionally, both transcripts are enriched in the developing pharynx and swim bladder. Our qPCR results indicate enhanced expression in the testes, along with increased expression in brain. Certainly, our experiments in the zebrafish model system with ankef1a and ankef1b provide a solid foundation for future studies to uncover the molecular pathways through which Ankef1 acts in both healthy and disease states.

Transcriptome profiling reveals male- and female-specific gene expression pattern and novel gene candidates for the control of sex determination and gonad development in Xenopus laevis

Xenopus laevis is an amphibian (frog) species widely used in developmental biology and genetics. To unravel the molecular machinery regulating sex differentiation of Xenopus gonads, we analyzed for the first time the transcriptome of developing amphibian gonads covering sex determination period. We applied microarray at four developmental stages: (i) NF50 (undifferentiated gonad during sex determination), (ii) NF53 (the onset of sexual differentiation of the gonads), (iii) NF56 (sexual differentiation of the gonads), and (iv) NF62 (developmental progression of differentiated gonads). Our analysis showed that during the NF50, the genetic female (ZW) gonads expressed more sex-specific genes than genetic male (ZZ) gonads, which suggests that a robust genetic program is realized during female sex determination in Xenopus. However, a contrasting expression pattern was observed at later stages (NF56 and NF62), when the ZW gonads expressed less sex-specific genes than ZZ gonads, i.e., more genes may be involved in further development of the male gonads (ZZ). We identified sexual dimorphism in the expression of several functional groups of genes, including signaling factors, proteases, protease inhibitors, transcription factors, extracellular matrix components, extracellular matrix enzymes, cell adhesion molecules, and epithelium-specific intermediate filaments. In addition, our analysis detected a sexually dimorphic expression of many uncharacterized genes of unknown function, which should be studied further to reveal their identity and if/how they regulate gonad development, sex determination, and sexual differentiation. Comparison between genes sex-specifically expressed in developing gonads of Xenopus and available transcriptome data from zebrafish, two reptile species, chicken, and mouse revealed significant differences in the genetic control of sex determination and gonad development. This shows that the genetic control of gonad development is evolutionarily malleable.

Changes in the morphology and gene expression of developing zebrafish gonads

Zebrafish gonadal sexual differentiation is an important but poorly understood subject. The difficulty in investigating zebrafish sexual development lies in its sex determination plasticity, the lack of morphological tools to distinguish juvenile females from males, and the lack of sex chromosomes in laboratory strains. Zebrafish sexual differentiation starts at around 8 days post-fertilization when germ cells start to proliferate. The number of germ cells determines the future sex of the gonad. Gonads with more germ cells differentiate into ovaries, whereas a reduced germ cell number leads to male-biased sexual differentiation. Genes controlling sexual differentiation in pre-meiotic gonads encode proteins such as transcription factors, the transforming growth factor (TGF)-β family of signaling proteins, and RNA-binding proteins. These proteins coordinately control germ cell proliferation/meiosis/maintenance and gonadal somatic cell differentiation, leading to stepwise differentiation of gonads. Morphological changes in differentiating gonads are characterized by the appearance of oocytes containing condensed chromatin, followed by incorporation of vitellogenin and oocyte maturation. Marker genes and morphological characteristics help distinguish the steps in zebrafish gonadal differentiation during this important sex-determining stage.

Concentration and timing of application reveal strong fungistatic effect of tebuconazole in a Daphnia-microparasitic yeast model

Given the importance of pollutant effects on host-parasite relationships and disease spread, the main goal of this study was to assess the influence of different exposure scenarios for the fungicide tebuconazole (concentration×timing of application) on a Daphnia-microparasitic yeast experimental system. Previous results had demonstrated that tebuconazole is able to suppress Metschnikowia bicuspidata infection at ecologically-relevant concentrations; here, we aimed to obtain an understanding of the mechanism underlying the anti-parasitic (fungicidal or fungistatic) action of tebuconazole. We exposed the Daphnia-yeast system to four nominal tebuconazole concentrations at four timings of application (according to the predicted stage of parasite development), replicated on two Daphnia genotypes, in a fully crossed experiment. An "all-or-nothing" effect was observed, with tebuconazole completely suppressing infection from 13.5μgl^-1 upwards, independent of the timing of tebuconazole application. A follow-up experiment confirmed that the suppression of infection occurred within a narrow range of tebuconazole concentrations (3.65-13.5μgl^-1), although a later application of the fungicide had to be compensated for by a slight increase in concentration to elicit the same anti-parasitic effect. The mechanism behind this anti-parasitic effect seems to be the inhibition of M. bicuspidata sporulation, since tebuconazole was effective in preventing ascospore production even when applied at a later time. However, this fungicide also seemed to affect the vegetative growth of the yeast, as demonstrated by the enhanced negative effect of the parasite (increasing mortality in one of the host genotypes) at a later time of application of tebuconazole, when no signs of infection were observed. Fungicide contamination can thus affect the severity and spread of disease in natural populations, as well as the inherent co-evolutionary dynamics in host-parasite systems.

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.

Heat-induced masculinization in domesticated zebrafish is family-specific and yields a set of different gonadal transcriptomes

Understanding environmental influences on sex ratios is important for the study of the evolution of sex-determining mechanisms and for evaluating the effects of global warming and chemical pollution. Fishes exhibit sexual plasticity, but the underlying mechanisms of environmental effects on their reproduction are unclear even in the well-established teleost research model, the zebrafish. Here we established the conditions to study the effects of elevated temperature on zebrafish sex. We showed that sex ratio response to elevated temperature is family-specific and typically leads to masculinization (female-to-male sex reversal), resulting in neomales. These results uncovered genotype-by-environment interactions that support a polygenic sex determination system in domesticated (laboratory) zebrafish. We found that some heat-treated fish had gene expression profiles similar to untreated controls of the same sex, indicating that they were resistant to thermal effects. Further, most neomales had gonadal transcriptomes similar to that of regular males. Strikingly, we discovered heat-treated females that displayed a normal ovarian phenotype but with a "male-like" gonadal transcriptome. Such major transcriptomic reprogramming with preserved organ structure has never been reported. Juveniles were also found to have a male-like transcriptome shortly after exposure to heat. These findings were validated by analyzing the expression of genes and signaling pathways associated with sex differentiation. Our results revealed a lasting thermal effect on zebrafish gonads, suggesting new avenues for detection of functional consequences of elevated temperature in natural fish populations in a global warming scenario.

Finding clues to the riddle of sex determination in zebrafish

How sex is determined has been one of the most intriguing puzzles in biology since antiquity. Although a fundamental process in most metazoans, there seems to be myriad of ways in which sex can be determined - from genetic to environmental sex determination. This variation is limited mainly to upstream triggers with the core of sex determination pathway being conserved. Zebrafish has gained prominence as a vertebrate model system to study development and disease. However, very little is known about its primary sex determination mechanism. Here we review our current understanding of the sex determination in zebrafish. Zebrafish lack identifiable heteromorphic sex chromosomes and sex is determined by multiple genes, with some influence from the environment. Recently, chromosome 4 has been identified as sex chromosome along with few sex-linked loci on chromosomes 5 and 16. The identities of candidate sex-linked genes, however, have remained elusive. Sex in zebrafish is also influenced by the number of meiotic oocytes in the juvenile ovary, which appear to instruct retention of the ovarian fate. The mechanism and identity of this instructive signal remain unknown. We hypothesize that sex in zebrafish is a culmination of combinatorial effects of the genome, germ cells and the environment with inputs from epigenetic factors translating the biological meaning of this interaction.

Characterizing Transcriptional Networks in Male Rainbow Darter (Etheostoma caeruleum) that Regulate Testis Development over a Complete Reproductive Cycle

Intersex is a condition that has been associated with exposure to sewage effluents in male rainbow darter (Etheostoma caeruleum). To better understand changes in the transcriptome that are associated with intersex, we characterized annual changes in the testis transcriptome in wild, unexposed fish. Rainbow darter males were collected from the Grand River (Ontario, Canada) in May (spawning), August (post-spawning), October (recrudescence), January (developing) and March (pre-spawning). Histology was used to determine the proportion of spermatogenic cell types that were present during each period of testicular maturation. Regression analysis determined that the proportion of spermatozoa versus spermatocytes in all stages of development (R² ≥ 0.58) were inversely related; however this was not the case when males were in the post-spawning period. Gene networks that were specific to the transition from developing to pre-spawning stages included nitric oxide biosynthesis, response to wounding, sperm cell function, and stem cell maintenance. The pre-spawning to spawning transition included gene networks related to amino acid import, glycogenesis, Sertoli cell proliferation, sperm capacitation, and sperm motility. The spawning to post-spawning transition included unique gene networks associated with chromosome condensation, ribosome biogenesis and assembly, and mitotic spindle assembly. Lastly, the transition from post-spawning to recrudescence included gene networks associated with egg activation, epithelial to mesenchymal transition, membrane fluidity, and sperm cell adhesion. Noteworthy was that there were a significant number of gene networks related to immune system function that were differentially expressed throughout reproduction, suggesting that immune network signalling has a prominent role in the male testis. Transcripts in the testis of post-spawning individuals showed patterns of expression that were most different for the majority of transcripts investigated when compared to the other stages. Interestingly, many transcripts associated with female sex differentiation (i.e. esr1, sox9, cdca8 and survivin) were significantly higher in the testis during the post-spawning season compared to other testis stages. At post-spawning, there were higher levels of estrogen and androgen receptors (esr1, esr2, ar) in the testis, while there was a decrease in the levels of sperm associated antigen 1 (spag1) and spermatogenesis associated 4 (spata4) mRNA. Cyp17a was more abundant in the testis of fish in the pre-spawning, spawning, and post-spawning seasons compared to those individuals that were recrudescent while aromatase (cyp19a) did not vary in expression over the year. This study identifies cell process related to testis development in a seasonally spawning species and improves our understanding regarding the molecular signaling events that underlie testicular growth. This is significant because, while there are a number of studies characterizing molecular pathways in the ovary, there are comparatively less describing transcriptomic patterns in the testis in wild fish.

Sex Change in Clownfish: Molecular Insights from Transcriptome Analysis

Sequential hermaphroditism is a unique reproductive strategy among teleosts that is displayed mainly in fish species living in the coral reef environment. The reproductive biology of hermaphrodites has long been intriguing; however, very little is known about the molecular pathways underlying their sex change. Here, we provide the first de novo transcriptome analyses of a hermaphrodite teleost´s undergoing sex change in its natural environment. Our study has examined relative gene expression across multiple groups-rather than just two contrasting conditions- and has allowed us to explore the differential expression patterns throughout the whole process. Our analysis has highlighted the rapid and complex genomic response of the brain associated with sex change, which is subsequently transmitted to the gonads, identifying a large number of candidate genes, some well-known and some novel, involved in the process. The present study provides strong evidence of the importance of the sex steroidogenic machinery during sex change in clownfish, with the aromatase gene playing a central role, both in the brain and the gonad. This work constitutes the first genome-wide study in a social sex-changing species and provides insights into the genetic mechanism governing social sex change and gonadal restructuring in protandrous hermaphrodites.

Maternal stress has divergent effects on gene expression patterns in the brains of male and female threespine stickleback

Maternal stress can have long-term effects on neurodevelopment that can influence offspring performance and population evolutionary trajectories. To examine the mechanistic basis for these neurodevelopmental effects of maternal stress, we used RNA-seq to assess differential gene expression across the brain transcriptome of adult male and female threespine stickleback (Gasterosteus aculeatus) from stressed and unstressed mothers. We identified sexually divergent effects of maternal stress on the brain transcriptome. In males, genes that were upregulated by maternal stress were enriched for processes involved in synaptic function and organization and steroid hormone-mediated signalling pathways, whereas in females genes that were upregulated by maternal stress were enriched for processes involved in protein translation and metabolic functions. The expression of several genes involved in the hypothalamic-pituitary-interrenal response to stress and epigenetic processes such as the regulation of DNA methylation patterns and miRNAs increased in males and not in females. These data suggest that maternal stress has markedly different effects on cellular pathways in the brains of male and female offspring of mothers that are exposed to stress, which could have important implications when assessing the long-term ecological and evolutionary impacts of stress across generations.

Investigation of septin biology in vivo using zebrafish

The zebrafish (Danio rerio) is an important animal model to study cell biology in vivo. Benefits of the zebrafish include a fully annotated reference genome, an easily manipulable genome (for example, by morpholino oligonucleotide or CRISPR-Cas9), and transparent embryos for noninvasive, real-time microscopy using fluorescent transgenic lines. Zebrafish have orthologues of most human septins, and studies using larvae were used to investigate the role of septins in vertebrate development. The zebrafish larva is also an established model to study the cell biology of infection and has recently been used to visualize septin assembly during bacterial infection in vivo. Here, we describe protocols for the study of septins in zebrafish, with emphasis on techniques used to investigate the role of septins in host defense against bacterial infection.

Sex differences in DNA methylation and expression in zebrafish brain: a test of an extended 'male sex drive' hypothesis

The sex drive hypothesis predicts that stronger selection on male traits has resulted in masculinization of the genome. Here we test whether such masculinizing effects can be detected at the level of the transcriptome and methylome in the adult zebrafish brain. Although methylation is globally similar, we identified 914 specific differentially methylated CpGs (DMCs) between males and females (435 were hypermethylated and 479 were hypomethylated in males compared to females). These DMCs were prevalent in gene body, intergenic regions and CpG island shores. We also discovered 15 distinct CpG clusters with striking sex-specific DNA methylation differences. In contrast, at transcriptome level, more female-biased genes than male-biased genes were expressed, giving little support for the male sex drive hypothesis. Our study provides genome-wide methylome and transcriptome assessment and sheds light on sex-specific epigenetic patterns and in zebrafish for the first time.

Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii)

BACKGROUND

Controlling and managing the breeding of bluefin tuna (Thunnus spp.) in captivity is an imperative step towards obtaining a sustainable supply of these fish in aquaculture production systems. Germ cell transplantation (GCT) is an innovative technology for the production of inter-species surrogates, by transplanting undifferentiated germ cells derived from a donor species into larvae of a host species. The transplanted surrogates will then grow and mature to produce donor-derived seed, thus providing a simpler alternative to maintaining large-bodied broodstock such as the bluefin tuna. Implementation of GCT for new species requires the development of molecular tools to follow the fate of the transplanted germ cells. These tools are based on key reproductive and germ cell-specific genes. RNA-Sequencing (RNA-Seq) provides a rapid, cost-effective method for high throughput gene identification in non-model species. This study utilized RNA-Seq to identify key genes expressed in the gonads of Southern bluefin tuna (Thunnus maccoyii, SBT) and their specific expression patterns in male and female gonad cells.

RESULTS

Key genes involved in the reproductive molecular pathway and specifically, germ cell development in gonads, were identified using analysis of RNA-Seq transcriptomes of male and female SBT gonad cells. Expression profiles of transcripts from ovary and testis cells were compared, as well as testis germ cell-enriched fraction prepared with Percoll gradient, as used in GCT studies. Ovary cells demonstrated over-expression of genes related to stem cell maintenance, while in testis cells, transcripts encoding for reproduction-associated receptors, sex steroids and hormone synthesis and signaling genes were over-expressed. Within the testis cells, the Percoll-enriched fraction showed over-expression of genes that are related to post-meiosis germ cell populations.

CONCLUSIONS

Gonad development and germ cell related genes were identified from SBT gonads and their expression patterns in ovary and testis cells were determined. These expression patterns correlate with the reproductive developmental stage of the sampled fish. The majority of the genes described in this study were sequenced for the first time in T. maccoyii. The wealth of SBT gonadal and germ cell-related gene sequences made publicly available by this study provides an extensive resource for further GCT and reproductive molecular biology studies of this commercially valuable fish.

Identification of Sex and Female's Reproductive Stage in Commercial Fish Species through the Quantification of Ribosomal Transcripts in Gonads

The estimation of maturity and sex of fish stocks in European waters is a requirement of the EU Data Collection Framework as part of the policy to improve fisheries management. On the other hand, research on fish biology is increasingly focused in molecular approaches, researchers needing correct identification of fish sex and reproductive stage without necessarily having in house the histological know-how necessary for the task. Taking advantage of the differential gene transcription occurring during fish sex differentiation and gametogenesis, the utility of 5S ribosomal RNA (5S rRNA) and General transcription factor IIIA (gtf3a) in the molecular identification of sex and gametogenic stage was tested in different economically-relevant fish species from the Bay of Biscay. Gonads of 9 fish species (, Atlantic, Atlantic-chub and horse mackerel, blue whiting, bogue, European anchovy, hake and pilchard and megrim), collected from local commercial fishing vessels were histologically sexed and 5S and 18S rRNA concentrations were quantified by capillary electrophoresis to calculate a 5S/18S rRNA index. Degenerate primers permitted cloning and sequencing of gtf3a fragments in 7 of the studied species. 5S rRNA and gtf3a transcript levels, together with 5S/18S rRNA index, distinguished clearly ovaries from testis in all of the studied species. The values were always higher in females than in males. 5S/18S rRNA index values in females were always highest when fish were captured in early phases of ovary development whilst, in later vitellogenic stages, the values decreased significantly. In megrim and European anchovy, where gonads in different oogenesis stages were obtained, the 5S/18S rRNA index identified clearly gametogenic stage. This approach, to the sexing and the quantitative non-subjective identification of the maturity stage of female fish, could have multiple applications in the study of fish stock dynamics, fish reproduction and fecundity and fish biology in general.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Transcriptomic Analyses Reveal Novel Genes with Sexually Dimorphic Expression in Yellow Catfish (Pelteobagrus fulvidraco) Brain

Yellow catfish (Pelteobagrus fulvidraco) is a pivotal freshwater aquaculture species in China. It shows sexual size dimorphism favoring male in growth. Whole transcriptome approach is required to get the overview of genetic toolkit for understanding the sex determination mechanism aiming at devising its monosex production. Beside gonads, the brain is also considered as a major organ for vertebrate reproduction. Transcriptomic analyses on the brain and of different developmental stages will provide the dynamic view necessary for better understanding its sex determination. In this regard, we have performed a de novo assembly of yellow catfish brain transcriptome by high throughput Illumina sequencing. A total number of 154,507 contigs were obtained with the lengths ranging from 201 to 27,822 bp and N50 of 2,101 bp, as well as 20,699 unigenes were identified. Of these unigenes, 13 and 54 unigenes were detected to be XY-specifically expressed genes (SEGs) for one and 2-year-old yellow catfish, while the corresponding numbers of XX-SEGs for those two stages were 19 and 13, respectively. Our work identifies a set of annotated genes that are candidate factors affecting sexual dimorphism as well as simple sequence repeat (SSR) and single nucleotide variation (SNV) in yellow catfish. To validate the expression patterns of the sex-related genes, we performed quantitative real-time PCR (qRT-PCR) indicating the reliability and accuracy of our analysis. The results in our study may enhance our understanding of yellow catfish sex determination and potentially help to improve the production of all-male yellow catfish for aquaculture.

Lasting effects of early exposure to temperature on the gonadal transcriptome at the time of sex differentiation in the European sea bass, a fish with mixed genetic and environmental sex determination

BACKGROUND

Sex in fish is plastic and in several species can be influenced by environmental factors. In sensitive species, elevated temperatures have a masculinizing effect. Previous studies on the effects of temperature on gene expression have been restricted to a few cognate genes, mostly related to testis or ovarian development, and analyzed in gonads once they had completed the process of sex differentiation. However, studies on the effect of temperature at the whole gonadal transcriptomic level are scarce in fish and, in addition, temperature effects at the time of sex differentiation at the transcriptomic level are also unknown. Here, we used the European sea bass, a gonochoristic teleost with a polygenic sex determination system influenced by temperature, and exposed larvae to elevated temperature during the period of early gonad formation. Transcriptomic analysis of the gonads was carried out about three months after the end of temperature exposure, shortly after the beginning of the process of sex differentiation.

RESULTS

Elevated temperature doubled the number of males with respect to untreated controls. Transcriptomic analysis of early differentiating female gonads showed how heat caused: 1) an up-regulation of genes related to cholesterol transport (star), the stress response (nr3c1) and testis differentiation (amh, dmrt, etc.), 2) a decrease in the expression of genes related to ovarian differentiation such as cyp19a1a, and 3) an increase in the expression of several genes related to epigenetic regulatory mechanisms (hdac11, dicer1, ehmt2, jarid2a, pcgf2, suz12, mettl22).

CONCLUSIONS

Taken together, the results of this study contribute to the understanding of how the early environment sets permanent changes that result in long-lasting consequences, in this case in the sexual phenotype. Results also show the usefulness of comparing the effects of heat on the behavior of cognate genes related to sex differentiation as well as that of genes involved in establishing and maintaining cell identity through epigenetic mechanisms.

Zebrafish sexual behavior: role of sex steroid hormones and prostaglandins

Background

Mating behavior differ between sexes and involves gonadal hormones and possibly sexually dimorphic gene expression in the brain. Sex steroids and prostaglandin E₂ (PGE₂) have been shown to regulate mammalian sexual behavior. The present study was aimed at determining whether exposure to sex steroids and prostaglandins could alter zebrafish sexual mating behavior.

Methods

Mating behavior and successful spawning was recorded following exposure to 17β-estradiol (E2), 11-ketotestosterone (11-KT), prostaglandin D₂ (PGD₂) and PGE₂ via the water. qRT-PCR was used to analyze transcript levels in the forebrain, midbrain, and hindbrain of male and female zebrafish and compared to animals exposed to E2 via the water.

Results

Exposure of zebrafish to sex hormones resulted in alterations in behavior and spawning when male fish were exposed to E2 and female fish were exposed to 11-KT. Exposure to PGD_2, and PGE₂ did not alter mating behavior or spawning success. Determination of gene expression patterns of selected genes from three brain regions using qRT-PCR analysis demonstrated that the three brain regions differed in gene expression pattern and that there were differences between the sexes. In addition, E2 exposure also resulted in altered gene transcription profiles of several genes.

Conclusions

Exposure to sex hormones, but not prostaglandins altered mating behavior in zebrafish. The expression patterns of the studied genes indicate that there are large regional and gender-based differences in gene expression and that E2 treatment alter the gene expression pattern in all regions of the brain.

Electronic supplementary material

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

Co-localization of growth QTL with differentially expressed candidate genes in rainbow trout

We tested whether genes differentially expressed between large and small rainbow trout co-localized with familial QTL regions for body size. Eleven chromosomes, known from previous work to house QTL for weight and length in rainbow trout, were examined for QTL in half-sibling families produced in September (1 XY male and 1 XX neomale) and December (1 XY male). In previous studies, we identified 108 candidate genes for growth expressed in the liver and white muscle in a subset of the fish used in this study. These gene sequences were BLASTN aligned against the rainbow trout and stickleback genomes to determine their location (rainbow trout) and inferred location based on synteny with the stickleback genome. Across the progeny of all three males used in the study, 63.9% of the genes with differential expression appear to co-localize with the QTL regions on 6 of the 11 chromosomes tested in these males. Genes that co-localized with QTL in the mixed-sex offspring of the two XY males primarily showed up-regulation in the muscle of large fish and were related to muscle growth, metabolism, and the stress response.

Molecular characterization and expression profile of three Fem-1 genes in Eriocheir sinensis provide a new insight into crab sex-determining mechanism

The FEM-1 protein of Caenorhabditis elegans plays a crucial role in the nematode sex-determination pathway. Here, we reported the characterization of three members of Fem-1 gene family in Eriocheir sinensis (designated EsFem-1a, EsFem-1b, and EsFem-1c), which were homologs of the nematode FEM-1 protein. The amino acid sequences of EsFem-1a, EsFem-1b, and EsFem-1c contained eight, nine, and eight ankyrin repeats, respectively. None of the ankyrin repeats had its own specific signature, and the evolution of ankyrin repeat was not completely independent. The predicted three-dimensional structure of EsFem-1 proteins exhibited highly similar superhelical conformation, especially the N-terminal six contiguous ankyrin repeats, which provided a binding surface for the protein-protein interaction. Phylogenetic tree based on the amino acid sequences revealed that EsFem-1a, EsFem-1b, and EsFem-1c were divided into three obvious separated clades. EsFem-1 genes were highly expressed in fertilized egg, 2-4 cell and blastula stage comparing with larval stage (P<0.01), which suggested they might be maternal genes. They also showed a certain degree of sexually dimorphic expression in some tissues. Notably, the highest expression of EsFem-1a was in the hepatopancreas, with EsFem-1b in testes and EsFem-1c in muscle (P<0.05), which indicated their potential role in a broad array of tissues. In addition, the genes initially involved in sex differentiation were not limited to those specifically expressed in the developing gonad. Taken together, these results suggested that EsFem-1 might function in crab early sex determination and late gonad development. The identification of Fem-1 gene family in E. sinensis provides a new insight into crab sex-determination mechanism.

Construction, De-Novo Assembly and Analysis of Transcriptome for Identification of Reproduction-Related Genes and Pathways from Rohu, Labeo rohita (Hamilton)

Rohu is a leading candidate species for freshwater aquaculture in South-East Asia. Unlike common carp the monsoon breeding habit of rohu restricts its seed production beyond season indicating strong genetic control over spawning. Genetic information is limited in this regard. The problem is exacerbated by the lack of genomic-resources. We identified 182 reproduction-related genes previously by Sanger-sequencing which were less to address the issue of seasonal spawning behaviour of this important carp. Therefore, the present work was taken up to generate transcriptome profile by mRNAseq. 16 GB, 72 bp paired end (PE) data was generated from the pooled-RNA of twelve-tissues from pre-spawning rohu using IlluminaGA-II-platform. There were 64.97 million high-quality reads producing 62,283 contigs and 88,612 numbers of transcripts using velvet and oases programs, respectively. Gene ontology annotation identified 940 reproduction-related genes consisting of 184 mainly associated with reproduction, 223 related to hormone-activity and receptor-binding, 178 receptor-activity and 355 embryonic-development related-proteins. The important reproduction-relevant pathways found in KEGG analysis were GnRH-signaling, oocyte-meiosis, steroid-biosynthesis, steroid-hormone biosynthesis, progesterone-mediated oocyte-maturation, retinol-metabolism, neuroactive-ligand-receptor interaction, neurotrophin-signaling and photo-transduction. Twenty nine simple sequence repeat containing sequences were also found out of which 12 repeat loci were polymorphic with mean expected-&-observed heterozygosity of 0.471 and 0.983 respectively. Quantitative RT-PCR analyses of 13-known and 6-unknown transcripts revealed differences in expression level between preparatory and post-spawning phase. These transcriptomic sequences have significantly increased the genetic-&-genomic resources for reproduction-research in Labeo rohita.

Comparison of growth-related traits and gene expression profiles between the offspring of neomale (XX) and normal male (XY) rainbow trout

All-female lines of fish are created by crossing sex reversed (XX genotype) males with normal females. All-female lines avoid the deleterious phenotypic effects that are typical of precocious maturation in males. To determine whether all-female and mixed sex populations of rainbow trout (Oncorhynchus mykiss) differ in performance, we compared the growth and gene expression profiles in progeny groups produced by crossing a XX male and a XY male to the same five females. Body weight and length were measured in the resulting all-female (XX) and mixed sex (XX/XY) offspring groups. Microarray experiments with liver and white muscle were used to determine if the gene expression profiles of large and small XX offspring differ from those in large and small XX/XY offspring. We detected no significant differences in body length and weight between offspring groups but XX offspring were significantly less variable in the value of these traits. A large number of upregulated genes were shared between the large XX and large XX/XY offspring; the small XX and small XX/XY offspring also shared similar expression profiles. No GO category differences were seen in the liver or between the large XX and large XX/XY offspring in the muscle. The greatest differences between the small XX and small XX/XY offspring were in the genes assigned to the "small molecule metabolic process" and "cellular metabolic process" GO level 3 categories. Similarly, genes within these categories as well as the category "macromolecule metabolic process" were more highly expressed in small compared to large XX fish.

Limited sex-biased neural gene expression patterns across strains in Zebrafish (Danio rerio)

Background

Male and female vertebrates typically differ in a range of characteristics, from morphology to physiology to behavior, which are influenced by factors such as the social environment and the internal hormonal and genetic milieu. However, sex differences in gene expression profiles in the brains of vertebrates are only beginning to be understood. Fishes provide a unique complement to studies of sex differences in mammals and birds given that fish show extreme plasticity and lability of sexually dimorphic characters and behaviors during development and even adulthood. Hence, teleost models can give additional insight into sexual differentiation. The goal of this study is to identify neurotranscriptomic mechanisms for sex differences in the brain.

Results

In this study we examined whole-brain sex-biased gene expression through RNA-sequencing across four strains of zebrafish. We subsequently conducted systems level analyses by examining gene network dynamics between the sexes using weighted gene coexpression network analysis. Surprisingly, only 61 genes (approximately 0.4% of genes analyzed) showed a significant sex effect across all four strains, and 48 of these differences were male-biased. Several of these genes are associated with steroid hormone biosynthesis. Despite sex differences in a display of stress-related behaviors, basal transcript levels did not predict the intensity of the behavioral display. WGCNA revealed only one module that was significantly associated with sex. Intriguingly, comparing intermodule dynamics between the sexes revealed only moderate preservation. Further we identify sex-specific gene modules.

Conclusions

Despite differences in morphology, physiology, and behavior, there is limited sex-biased neural gene expression in zebrafish. Further, genes found to be sex-biased are associated with hormone biosynthesis, suggesting that sex steroid hormones may be key contributors to sexual behavioral plasticity seen in teleosts. A possible mechanism is through regulating specific brain gene networks.

Electronic supplementary material

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

Mugilid fish are sentinels of exposure to endocrine disrupting compounds in coastal and estuarine environments

Effects on fish reproduction can result from a variety of toxicity mechanisms first operating at the molecular level. Notably, the presence in the environment of some compounds termed endocrine disrupting chemicals (EDCs) can cause adverse effects on reproduction by interfering with the endocrine system. In some cases, exposure to EDCs leads to the animal feminization and male fish may develop oocytes in testis (intersex condition). Mugilid fish are well suited sentinel organisms to study the effects of reproductive EDCs in the monitoring of estuarine/marine environments. Up-regulation of aromatases and vitellogenins in males and juveniles and the presence of intersex individuals have been described in a wide array of mullet species worldwide. There is a need to develop new molecular markers to identify early feminization responses and intersex condition in fish populations, studying mechanisms that regulate gonad differentiation under exposure to xenoestrogens. Interestingly, an electrophoresis of gonad RNA, shows a strong expression of 5S rRNA in oocytes, indicating the potential of 5S rRNA and its regulating proteins to become useful molecular makers of oocyte presence in testis. Therefore, the use of these oocyte markers to sex and identify intersex mullets could constitute powerful molecular biomarkers to assess xenoestrogenicity in field conditions.

The sex-specific transcriptome of the hermaphrodite sparid sharpsnout seabream (Diplodus puntazzo)

Background

Teleosts are characterized by a remarkable breadth of sexual mechanisms including various forms of hermaphroditism. Sparidae is a fish family exhibiting gonochorism or hermaphroditism even in closely related species. The sparid Diplodus puntazzo (sharpsnout seabream), exhibits rudimentary hermaphroditism characterized by intersexual immature gonads but single-sex mature ones. Apart from the intriguing reproductive biology, it is economically important with a continuously growing aquaculture in the Mediterranean Sea, but limited available genetic resources. Our aim was to characterize the expressed transcriptome of gonads and brains through RNA-Sequencing and explore the properties of genes that exhibit sex-biased expression profiles.

Results

Through RNA-Sequencing we obtained an assembled transcriptome of 82,331 loci. The expression analysis uncovered remarkable differences between male and female gonads, while male and female brains were almost identical. Focused search for known targets of sex determination and differentiation in vertebrates built the sex-specific expression profile of sharpsnout seabream. Finally, a thorough genetic marker discovery pipeline led to the retrieval of 85,189 SNPs and 29,076 microsatellites enriching the available genetic markers for this species.

Conclusions

We obtained a nearly complete source of transcriptomic sequence as well as marker information for sharpsnout seabream, laying the ground for understanding the complex process of sex differentiation of this economically valuable species. The genes involved include known candidates from other vertebrate species, suggesting a conservation of the toolkit between gonochorists and hermaphrodites.

Electronic supplementary material

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