De novo transcriptome analysis to search for sex-differentiation genes in the Siberian sturgeon

hsd17b1 is a key gene for ovarian differentiation of the Siberian sturgeon

This is the first work using gonads from undifferentiated, genetically-sexed Siberian sturgeon describing expression changes in genes related to steroid synthesis and female and male sex differentiation. One factor identified as relevant for ovarian differentiation was the gene coding for the enzyme Hsd17b1, which converts estrone into estradiol-17β. hsd17b1 was highly activated in female gonads at 2.5 months of age, around the onset of sex differentiation, preceding activation of two other genes involved in estrogen production (cyp19a1 and foxl2). hsd17b1 was also strongly repressed in males. Two known foxl2 paralogs are found in Siberian sturgeon-foxl2 and foxl2l-but only foxl2 appeared to be associated with ovarian differentiation. With regard to the male pathway, neither 11-oxygenated androgens nor classic male genes (amh, dmrt1, sox9, and dhh) were found to be involved in male sex differentiation, leaving open the question of which genes participate in early male gonad development in this ancient fish. Taken together, these results indicate an estrogen-dependence of female sex differentiation and 11-oxygenated androgen-independence of male sex differentiation.

Genome-wide association study reveals genomic loci of sex differentiation and gonadal development in Plectropomus leopardus

Introduction: Plectropomus leopardus, a commercially significant marine fish, is primarily found in the Western Pacific regions and along the coast of Southeast Asia. A thorough analysis of the molecular mechanisms involved in sex differentiation is crucial for gaining a comprehensive understanding of gonadal development and improving sex control breeding. However, the relevant fundamental studies of P. leopardus are relatively lacking. Methods: In this study, a genome-wide association study (GWAS) was conducted to investigate the genetic basis mechanism of sex differentiation and gonadal developmental traits in P. leopardus utilizing about 6,850,000 high-quality single-nucleotide polymorphisms (SNPs) derived from 168 individuals (including 126 females and 42 males) by the genome-wide efficient mixed-model association (GEMMA) algorithm. Results: The results of these single-trait GWASs showed that 46 SNP loci (-log10 p > 7) significantly associated with sex differentiation, and gonadal development traits were distributed in multiple different chromosomes, which suggested the analyzed traits were all complex traits under multi-locus control. A total of 1,838 potential candidate genes were obtained by considering a less-stringent threshold (-log10 p > 6) and ±100 kb regions surrounding the significant genomic loci. Moreover, 31 candidate genes were identified through a comprehensive analysis of significant GWAS peaks, gene ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, including taf7, ddx6, apoeb, sgk1, a2m, usf1, hsd3b7, dll4, xbp1, tet3, esr1, and gli3. These trait-associated genes have been shown to be involved in germline development, male sex differentiation, gonad morphogenesis, hormone receptor binding, oocyte development, male gonad development, steroidogenesis, estrogen-synthetic pathway, etc. Discussion: In the present study, multiple genomic loci of P. leopardus associated with sex differentiation and gonadal development traits were identified for the first time by using GWAS, providing a valuable resource for further research on the molecular genetic mechanism and sex control in P. leopardus. Our results also can contribute to understanding the genetic basis of the sex differentiation mechanism and gonadal development process in grouper fish.

Denovo RNA-Seq analysis of ovary and testis reveals potential differentially expressed transcripts associated with gonadal unsynchronization development in Onychostoma macrolepis

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered wild species. Their endangered extinction might be due to their low fertility. To further illustrate the molecular mechanism of gonad development of the male and female O. macrolepis, the present study carried out de novo testicular and ovarian transcriptome sequencing. By comparing ovary and testis, 30,869 differentially expressed unigenes (9870 in female, 20999 in male) were identified. In addition, KEGG and GO analysis suggested that the Hedgehog signaling pathway have important roles in testis maintenance and spermatogenesis, whereas the Hippo signaling pathway and Wnt signaling pathway are likely to participate in ovary maintenance. RT-qPCR analysis results were consistent with transcriptome sequencing that all of gender differentiation-related genes (FOXL2, GDF9, WNT4, CYP19A1, SOX9 and GATA4), temperature-enriched genes (NOVA1, CTGF and NR4A1), clock-related genes (PER2, PER3, CRY1, CRY2, BMAL1 and CIPC) were significantly differential expression in testis compared with ovaries. Taken together, these results revealed a potential molecular mechanism that low fertility of the O. macrolepis might strong correlate with the gonadal dyssynchrony development of the male and female, which might provide theoretical basis and technical support for artificial reproduction and germplasm resource protection of the O. macrolepis.

Distinct patterns of gene expression during regeneration and asexual reproduction in the annelid Pristina leidyi

Regeneration, the ability to replace lost body parts, is a widespread phenomenon in the animal kingdom often connected to asexual reproduction or fission, since the only difference between the two appears to be the stimulus that triggers them. Both developmental processes have largely been characterized; however, the molecular toolkit and genetic mechanisms underlying these events remain poorly unexplored. Annelids, in particular the oligochaete Pristina leidyi, provide a good model system to investigate these processes as they show diverse ways to regenerate, and can reproduce asexually through fission under laboratory conditions. Here, we used a comparative transcriptomics approach based on RNA-sequencing and differential gene expression analyses to understand the molecular mechanisms involved in anterior regeneration and asexual reproduction. We found 291 genes upregulated during anterior regeneration, including several regeneration-related genes previously reported in other annelids such as frizzled, paics, and vdra. On the other hand, during asexual reproduction, 130 genes were found upregulated, and unexpectedly, many of them were related to germline development during sexual reproduction. We also found important differences between anterior regeneration and asexual reproduction, with the latter showing a gene expression profile more similar to that of control individuals. Nevertheless, we identified 35 genes that were upregulated in both conditions, many of them related to cell pluripotency, stem cells, and cell proliferation. Overall, our results shed light on the molecular mechanisms that control anterior regeneration and asexual reproduction in annelids and reveal similarities with other animals, suggesting that the genetic machinery controlling these processes is conserved across metazoans.

Sex-relevant genes in the embryo stage of Chinese soft-shelled turtles as revealed by RNA-Seq analysis

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an economically important and rare freshwater aquaculture species in China. The male turtles have better economic benefits than females due to their faster growth speed. In order to explore the mechanism of sex differentiation and determination of P. sinensis, transcriptome analysis was carried out using embryo samples from the sex determination, sex differentiation and shelling stages, respectively. The gender types of embryos were identified by genotyping and histological analysis. In all, 277,230 and 273,859 genes were identified from embryos of male and female turtles, respectively. Cluster analysis of gene expression patterns exhibited trends for all differentially expressed genes (DEGs), which can be classified into three subclusters. In the sex determination, sex differentiation, and shelling stages 975, 43, and 300 DEGs, respectively, were identified by comparing the male and female groups. KEGG analysis was used to explore the function of the DEGs. Fgf9, foxl2, mex3b, sox7,tgfβ-3, wnt4, wnt7a, and wt1 were upregulated in ovary development stages and chfr, ampk, aldh2, dao, glt1d1, hvcn1, psmd9, srsf9 and ubc were upregulated in testis development stages, indicating that these genes play important roles in the sex development of P. sinensis. The results of this study provide preliminary insights into the sex determination and differentiation of Chinese soft-shelled turtles.

Isolation and Characterization of Highly Pure Type A Spermatogonia From Sterlet (Acipenser ruthenus) Using Flow-Cytometric Cell Sorting

Sturgeons are among the most ancient linages of actinopterygians. At present, many sturgeon species are critically endangered. Surrogate production could be used as an affordable and a time-efficient method for endangered sturgeons. Our study established a method for identifying and isolating type A spermatogonia from different developmental stages of testes using flow cytometric cell sorting (FCM). Flow cytometric analysis of a whole testicular cell suspension showed several well-distinguished cell populations formed according to different values of light scatter parameters. FCM of these different cell populations was performed directly on glass slides for further immunocytochemistry to identify germ cells. Results showed that the cell population in gate P1 on a flow cytometry plot (with high forward scatter and high side scatter parameter values) contains the highest amount of type A spermatogonia. The sorted cell populations were characterized by expression profiles of 10 germ cell specific genes. The result confirmed that setting up for the P1 gate could precisely sort type A spermatogonia in all tested testicular developmental stages. The P2 gate, which was with lower forward scatter and side scatter values mostly, contained type B spermatogonia at a later maturing stage. Moreover, expressions of plzf, dnd, boule, and kitr were significantly higher in type A spermatogonia than in later developed germ cells. In addition, plzf was firstly found as a reliable marker to identify type A spermatogonia, which filled the gap of identification of spermatogonial stem cells in sterlet. It is expected to increase the efficiency of germ stem cell culture and transplantation with plzf identification. Our study thus first addressed a phenotypic characterization of a pure type A spermatogonia population in sterlet. FCM strategy can improve the production of sturgeons with surrogate broodstock and further the analysis of the cellular and molecular mechanisms of sturgeon germ cell development.

Serum amyloid A is a positive acute phase protein in Russian sturgeon challenged with Aeromonas hydrophila

The immune system of sturgeons, one of the most ancient and economically valuable fish worldwide, is poorly understood. The lack of molecular tools and data about infection biomarkers hinders the possibility to monitor sturgeon health during farming and detect infection outbreaks. To tackle this issue, we mined publicly available transcriptomic datasets and identified putative positive acute-phase proteins (APPs) of Russian sturgeons that could be induced by a bacterial infection and monitored using non-invasive methods. Teleost literature compelled us to focus on five promising candidates: hepcidin, a warm acclimation associated hemopexin, intelectin, serum amyloid A protein (SAA) and serotransferrin. Among them, SAA was the most upregulated protein at the mRNA level in the liver of sturgeons challenged with heat-inactivated or live Aeromonas hydrophila. To assess whether this upregulation yielded increasing SAA levels in circulation, we developed an in-house ELISA to quantify SAA levels in sturgeon serum. Circulating SAA rose upon bacterial challenge and positively correlated with hepatic saa expression. This is the first time serum SAA has been quantified in an Actinopterygii fish. Since APPs vary across different fish species, our work sheds light on sturgeon acute-phase response, revealing that SAA is a positive APP with potential value as infection biomarker.

Siberian sturgeon multi-tissue reference transcriptome database

Motivation: Siberian sturgeon is a long lived and late maturing fish farmed for caviar production in 50 countries. Functional genomics enable to find genes of interest for fish farming. In the absence of a reference genome, a reference transcriptome is very useful for sequencing based functional studies.

Results: We present here a high-quality transcriptome assembly database built using RNA-seq reads coming from brain, pituitary, gonadal, liver, stomach, kidney, anterior kidney, heart, embryonic and pre-larval tissues. It will facilitate crucial research on topics such as puberty, reproduction, growth, food intake and immunology. This database represents a major contribution to the publicly available sturgeon transcriptome reference datasets.

Availability: The database is publicly available at http://siberiansturgeontissuedb.sigenae.org

Supplementary information:  Supplementary data are available at Database online.

RNA-Seq transcriptome data of undifferentiated and differentiated gonads of Siberian sturgeon

RNA-Seq transcriptome data from twenty Siberian sturgeon gonads at different developmental stages is described: ten undifferentiated gonads, six gonads of immature males and four gonads from immature females. Siberian sturgeon, Acipenser baerii, is long-lived, late-maturing fish farmed in 50 countries but its production remains on a craftsman scale when compared to industrial species. Sturgeon genetic and physiological studies are less developed than for industrial fish. The data presented hereafter enables fundamental studies on the regulatory mechanisms of sturgeon gonad development, which can further be applied both in aquaculture and in fundamental research.

Transcriptome analysis of ovarian maturation in a chondrostei Chinese sturgeon Acipenser sinensis

Chinese sturgeon (Acipenser sinensis) with an evolutionary history of over 200 million years, has a long lifespan, and an extremely late and asynchronous sexual maturation (8-18 years for males and 14-26 years for females), resulting in the difficulty of mature adult culture. However, little is known about the regulatory mechanisms of the transition among ovarian maturation stages in the Chinese sturgeon. We performed de novo transcriptome sequencing of the Chinese sturgeon at different ovarian maturation stages (FII, FIII, and FIV). The number of differentially expressed genes (DEGs) between FII and FIII/FIV (33,517/34,217) was more than that between FIII and FIV (22,123), suggesting that the transition from FII to FIII/FIV is more important than that from FIII to FIV for ovarian maturation. The number of upregulated genes was more than that of the downregulated genes, suggesting that increased gene expression was involved in the transition from FII to FIII/FIV. The representative pathways of DEGs were steroid biosynthesis, fatty acid biosynthesis, fatty acid elongation, glycerolipid metabolism, biosynthesis of unsaturated fatty acid, and α-linolenic acid metabolism. The differential expressions from the transcriptome sequencing were validated with real-time reverse-transcription polymerase chain reaction. We also found 13 genes in sexual development, female sex determination, gonadal development, ovarian maturation, ovarian follicle development, and oocyte development pathways, which were differently expressed among fish at FII, FIII, and FIV. We suggest that enhanced synthesis of steroid, unsaturated fatty acid, and α-linolenic acid may contribute to ovarian maturation of the Chinese sturgeon. These potential determinants provide a glimpse of the molecular architecture of ovary development in sturgeons.

2,2',4,4'-tetrabromodipheny ether (BDE-47) disrupts gonadal development of the Africa clawed frog (Xenopus laevis)

Previous studies have shown that BDE-47, one of the most abundant polybrominated diphenyl ethers (PBDEs) congeners, has a weak estrogenic activity, but it has remained unclear whether BDE-47 disrupts gonadal development and causes male-to-female sex reversal in lower vertebrates, with limited and controversial data. The present study aimed to determine the effects of BDE-47 on gonadal development in Xenopus laevis, a model amphibian species for studying adverse effects of estrogenic chemicals on reproductive development. X. laevis at stage 45/46 were exposed to BDE-47 (0.5, 5, 50 nM) in semi-static system, with 1 nM 17β-estradiol (E2) as the positive control. When reaching stage 53, tadpoles were examined for gonadal morphology, histology and sex-dimorphic gene expression. The phenotypic sex (gonadal morphology and histology) of each BDE-47-treated tadpole matched its genetic sex, showing no sex-reversal, whereas one half of genetic males treated with E2 displayed ovarian-like features. However, some genetic males (26%) in the 50 nM BDE-47 treatment group were found to contain more germ cells clumping together in the medulla, along with an increasing tendency of the gonad length/kidney length ratio in males, resembling feminizing outcomes of E2. These observations seem to suggest that BDE-47 exerted weak feminizing effects. However, BDE-47 induced increases in expression of both female-biased genes and male-biased genes in two sexes, which disagrees with feminizing outcomes, suggesting complicated effects of BDE-47 on gonadal development. Taken together, all results demonstrate that nanomolar BDE-47 disrupted gonadal development and exerted weak feminizing effects, but not resulted in male-to-female sex reversal in X. laevis.

Transcriptome sequencing of hybrid bester sturgeon: Responses to poly (I:C) in the context of comparative immunogenomics

Sturgeons represent a substantial scientific interest due to their high economic value, endangered status and also as the most primitive group of ray-finned fishes. Rapid progress in knowledge of sturgeon immunity was achieved recently with use of RNA sequencing. We report transcriptome sequencing of gill, head kidney, and spleen of bester sturgeon (a hybrid of beluga Huso huso and sterlet Acipenser ruthenus) injected with synthetic double-stranded RNA (polyI:C). The composition of transcriptome and responses to treatment were examined in the context of comparative genomics with focus on immune genes. Sturgeon transcripts matched to 21.5 k different proteins (blastx). With reference to Atlantic salmon, the functional groups and pathways of the immune system were uniformly represented: at average 36.5 ± 0.8% genes were found. Immune genes comprise a significant fraction of transcriptome. Among twenty genes with highest transcription levels, five are specialized immune genes and two encode heme and iron binding proteins (serotransferrin and hemopexin) also known as acute phase proteins. Challenge induced multiple functional groups including apoptosis, cell cycle and a number of metabolic pathways. Treatment stimulated innate antiviral immunity, which is well conserved between sturgeon and salmon, the most responsive genes were mx, rsad2 (viperin), interferon induced protein 44 and protein with tetratricopeptide repeats 5, cd87 and receptor transporting protein 3. Results added to knowledge of immune phylogeny. Gain and loss of genes was assessed by comparison with genomes from different phylogenetic groups. Among differentially expressed genes, percentage of acquired and lost genes was much lower in comparison with genes present in all vertebrates. Innate antiviral immunity was subject to the greatest changes in evolution of jawed vertebrates. A significant fraction of genes (15%) was lost in mammals and only half of genes is annotated in public databases as involved in antiviral responses. Change of function may have an important role in evolution of immunity together with gain and loss of genes.

FoxL2 combined with Cyp19a1a regulate the spawning upstream migration in Coilia nasus

FoxL2 is a member of the forkhead/HNF-3-related family of transcription factors which provides tissue specific gene regulation. It is known to regulate ovarian aromatase, which plays a crucial role in ovarian development and mature. To understand the role of FoxL2/ovarian aromatase encoded gene Cyp19a1a during ovarian development and recrudescence, we identified cDNA characteristics of FoxL2 and Cyp19a1a, analyzed its temporal expression both at transcript and protein levels in the anadromous fish, Coilia nasus. Tissue distribution pattern revealed that FoxL2 mRNA expression level was highest in ovary, while Cyp19a1a mRNA was highest in brain. During the upstream migration cycle, in ovary, the FoxL2 mRNA temporal expression peaked at the multiplication stage (stage III in May), the Cyp19a1a mRNA expression peaked at the onset stage (stage I in March). It was found that their mRNA transcripts were maintained at high level during the migration stage (from stage I in March to stage VI in July). Additionally, the strongest immunolabeling positive signals of Cyp19a1a and FoxL2 proteins were mainly found in the cytoplasm of olfactory bulb cell, stratum granulare and neurogliocyte cells and development stage oocytes. Data indicated that FoxL2 and Cyp19a1a were inducible and functional in the C. nasus ovary development and migration process. Therefore, the present results can be regarded as evidence for indispensable roles of FoxL2 and Cyp19a1a in the ovary development and migratory behavior at gene expression patterns and encoded protein distribution level.

Green light irradiation during sex differentiation induces female-to-male sex reversal in the medaka Oryzias latipes

This study investigated whether irradiation of a specific light wavelength could affect the sex differentiation of fish. We first found that the photoreceptor genes responsible for receiving red, green, and ultraviolet light were expressed in the eyes of medaka during the sex differentiation period. Second, we revealed that testes developed in 15.9% of genotypic females reared under green light irradiation. These female-to-male sex-reversed fish (i.e. neo-males) showed male-specific secondary sexual characteristics and produced motile sperm. Finally, progeny tests using the sperm of neo-males (XX) and eggs of normal females (XX) revealed that all F1 offspring were female, indicating for the first time in animals that irradiation with light of a specific wavelength can trigger sex reversal.