Gonadal transcriptome analysis of the common carp, Cyprinus carpio: Identification of differentially expressed genes and SSRs

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

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

Genetic diversity of the common carp black strain population based on mtDNA (D-loop and cytb)

The common strain black carp (Cyprinus carpio var. baisenensis) is a culturally important carp strain that is raised and cultured in Guangxi Province, China. Its color reflects the interactions between the Burau people and their surrounding environment. The population of the common carp black strain was isolated and cultured in a rice-fish integration system. To explore the genetic diversity and protection of germplasm resources, we analyzed mitochondrial DNA (mtDNA) sequences, specifically the displacement loop (D-loop) and cytochrome b (Cytb), using single-nucleotide polymorphisms (SNP). We compared these sequences with those from four other local common carp populations. The study included a total of 136 adult common carps from five strain populations: the common black carp strain (HJ = 31), Jian (F = 30), Heilongjiang (H = 10), Songpu (S = 31), and Saijiang (SJ = 34). The results of the Cytb and D-loop analyses showed that the Heilongjiang carp (H) and Saijiang (SJ) populations had the highest levels of haplotype diversity (0.867 ± 0.034785) and nucleotide diversity (π = 0.0063 ± 0.000137 and 0.0093 ± 0.000411), respectively. On the other hand, the Common carp black strain population (HJ) exhibited the lowest haplotype diversity in both Cytb and D-loop, with haplotype 2 being the most commonly observed among the populations. Private haplotypes dominated the five common carp populations, which were significantly different at P<0.001. Furthermore, analyzing the coefficient of genetic differentiation (Fst), the highest genetic difference was observed between Saijiang (SJ) and Heilongjiang (H) (Fst = 0.963), whereas the lowest was observed between Songpu (S) and the Common carp black strain population (HJ) (Fst = 0.019) for the Cytb gene sequences. For the D-loop, the Common carp black strain population (HJ) and Songpu (S) (Fst = 0.7) had the highest values, and Heilongjiang (H) and Common black carp strain (HJ) had an Fst of 0.125. Additionally, the AMOVA analysis revealed a higher level of variance for the Cytb and D-loop genes, indicating lower genetic diversity within the local carp community. On the other hand, the phylogenetic tree analysis showed that the five carp populations were closely related and formed a distinct cluster. The distinct cluster of populations suggests a common ancestor or recent gene flow, possibly due to geographic proximity or migration history, and unique genetic characteristics, possibly due to adaptations or selective pressures. The results of this study provide valuable insights into the genetic diversity of the common strain black carp, which can have implications for conservation, breeding programs, evolutionary studies, and fisheries management.

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

Identification and Characterization of Sex-Biased miRNAs in the Golden Pompano (Trachinotus blochii)

The golden pompano (Trachinotus blochii) is a marine fish of considerable commercial importance in China. It shows notable sexual size dimorphism; the growth rate of females is faster than that of males. Therefore, sex-biased research is of great importance in T. blochii breeding. However, there have been few studies on sex differentiation and mechanisms underlying sex determination in T. blochii. MicroRNAs (miRNAs) play crucial roles in sex differentiation and determination in animals. However, limited miRNA data are available on fish. In this study, two small RNA libraries prepared from the gonads of T. blochii were constructed and sequenced. The RNA-seq analysis yielded 1366 known and 69 novel miRNAs with 289 significantly differentially expressed miRNAs (p < 0.05). Gene ontology (GO) analysis confirmed that the TFIIA transcription factor complex (GO: 0005672) was the most significantly enriched GO term. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the differentially expressed miRNAs and target genes were mainly related to sex determination and gonadal developmental signaling pathways, specifically the Wnt signaling pathway, MAPK signaling pathway, and steroid biosynthetic pathway. MiRNA-mRNA co-expression network analysis strongly suggested a role for sex-biased miRNAs in sex determination/differentiation and gonadal development. For example, gata4, foxo3, wt1, and sf1 genes were found to be regulated by bta-miR-2898; esr2 and foxo3 by novel_176, and ar by oar-let-7b. Quantitative real-time polymerase chain reaction analysis of selected mRNAs and miRNAs validated the integrated analysis. This study established a set of sex-biased miRNAs that are potential regulatory factors in gonadal development in T. blochii. These results provide new insight into the function of miRNAs in sex differentiation and determination in T. blochii and highlight some key miRNAs for future studies.

Exposure of carbon nanotubes affects testis and brain of common carp

Carbon nanotubes production has been rapidly increasing for many potential applications, however, the environmental impact of this nanomaterial needs to be comprehended. The present work focused on unraveling the effects of single-walled carbon nanotubes (SWCNT) in the common carp, Cyprinus carpio. The physicochemical properties of SWCNT were analyzed with X-ray diffraction, Fourier transforms infra-red, UV-Vis absorption, transmission electron microscopy (TEM), and Raman spectroscopy before testing for exposure impact. The effects of SWCNT, were investigated by exposing to two doses viz., 10 and 50 μg/L, for 7 days in adult common carp, in vivo. Expression of key steroidogenic and transcription factor genes related to testis and brain were downregulated after the treatment. The concomitant decreases in serum testosterone and 11-ketotestosterone levels revealed the impact of SWCNT after exposure. Further, SWCNT exposure induced antioxidant enzymes namely glutathione-S-transferases, superoxide dismutase, and catalase in both testis and brain. Concurrently, histological and TEM analysis of testis revealed structural disarray. In addition, SWCNT treatment, in testicular and brain primary cell cultures decreased cell viability with an increase of reactive oxygen species levels, leading to a significant elevation of apoptotic cells. In line with this, low mitochondrial membrane potential and DNA damage were also observed during post SWCNT treatment. Taken together, transient exposure of SWCNT causes toxic effects and alters testicular and brain function in the common carp. Thus, the discharge of carbon nanotubes poses a greater risk to the aquatic environment warranting regulatory measures.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

sox19 regulates ovarian steroidogenesis in common carp

In teleost, ovarian steroidogenesis governed by the neuroendocrine system is also regulated by several transcription factors of gonadal origin. Investigating the synchronized interactions between the transcriptional and the hormonal factors is vital to comprehend the mechanisms that lead to gonadal differentiation. This study signifies the role of sry-related box (sox) 19 in ovarian steroidogenesis regulation of the common carp, Cyprinus carpio. Analysis of tissue distribution displayed higher sox19 expression in brain and ovary, and gonadal ontogeny showed higher expression of sox19 at 80 days post hatch (dph). Higher sox19 mRNA expression during spawning and increase of sox19 post human chorionic gonadotropin induction substantiate gonadotropin dependency. Estradiol-17β treatment but not 17α-methyl-di-hydroxy-testosterone to 50 dph common carp for inducing mono-sex, elevated sox19 expression substantially. Sox19 protein was observed in granulosa cells of the follicular layer in common carp ovary. Higher sox19 expression was detected in isolated granulosa and theca cells, in vitro. Transient gene silencing with sox19-siRNA caused downregulation of various ovary-related genes including those specific to activator protein-1 factors, fibroblast growth factors, wnt-signaling, steroidogenic genes along with certain transcription factors. Serum 17α, 20β-dihydroxy-4-pregnen-3-one and estradiol-17β reduced significantly post sox19 silencing, in vivo. Concomitantly, a decrease in aromatase activity was detected post sox19-siRNA treatment, in vivo. This study demonstrates the impact of sox19 in the regulation of common carp ovarian growth and steroidogenesis.

Functional characterization of two 20β-hydroxysteroid dehydrogenase type 2 homeologs from Xenopus laevis reveals multispecificity

The African clawed frog, Xenopus laevis, is a versatile model for biomedical research and is largely similar to mammals in terms of organ development, anatomy, physiology, and hormonal signaling mechanisms. Steroid hormones control a variety of processes and their levels are regulated by hydroxysteroid dehydrogenases (HSDs). The subfamily of 20β-HSD type 2 enzymes currently comprises eight members from teleost fish and mammals. Here, we report the identification of three 20β-HSD type 2 genes in X. tropicalis and X. laevis and the functional characterization of the two homeologs from X. laevis. X. laevis Hsd20b2.L and Hsd20b2.S showed high sequence identity with known 20β-HSD type 2 enzymes and mapped to the two subgenomes of the allotetraploid frog genome. Both homeologs are expressed during embryonic development and in adult tissues, with strongest signals in liver, kidney, intestine, and skin. After recombinant expression in human cell lines, both enzymes co-localized with the endoplasmic reticulum and catalyzed the conversion of cortisone to 20β-dihydrocortisone. Both Hsd20b2.L and Hsd20b2.S catalyzed the 20β-reduction of further C₂₁ steroids (17α-hydroxyprogesterone, progesterone, 11-deoxycortisol, 11-deoxycorticosterone), while only Hsd20b2.S was able to convert corticosterone and cortisol to their 20β-reduced metabolites. Estrone was only a poor and androstenedione no substrate for both enzymes. Our results demonstrate multispecificity of 20β-HSD type 2 enzymes from X. laevis similar to other teleost 20β-HSD type 2 enzymes. X. laevis 20β-HSD type 2 enzymes are probably involved in steroid catabolism and in the generation of pheromones for intraspecies communication. A role in oocyte maturation is unlikely.

Interactive role of Wnt signaling and Zn in regulating testicular function of the common carp, Cyprinus carpio

Wnt signaling is conserved among all species and plays a significant role in various cellular processes including reproduction. The present study identified significant involvement of wnt4a, wnt5b, and wnt8a signaling in the testicular growth of common carp,Cyprinus carpio. Predominant expression of wnt4a, wnt5b, and wnt8a was found in the gonads and Wnt4a was localized in spermatocytes and interstitial cells. Ontogeny and testicular phase-wise analysis signified the importance of wnt isofoms analyzed in this study. Specific pathway activation of Wnt signaling revealed that Wnt4a and Wnt5b act through non-canonical while Wnt8a prefers the canonical pathway. The Wnt signaling regulates several steroidogenic enzyme and testis-related genes which was confirmed by the Wnt blockade experiments. Incidentally, zinc (Zn) is an essential trace element involved in the progression of spermatogenesis in teleosts. In adult male carp, a single administration of Zn at different doses elevated the expression of Wnt and Zn transporter genes and a single-dose (30 μg/g body weight [BW]) of Zn treatment elevated steroidogenic enzyme and testis-related genes which coincided with elevated androgens. Conversely, single-dose administration of Zn chelator to the Zn administered (30 μg/g BW) fish reversed the effects emphasizing a prominent role of Zn in the testicular function perhaps through Wnt signaling. Similar effects were observed in the in vitro experiments using the Zn chelator. Bioaccumulation of Zn and histological analysis revealed the importance of Zn in progression of spermatogenesis and sperm motility. Various assays related to cell viability and proliferation exhibited the role of Zn in promoting spermatogenic cell progression. Flow cytometric analysis confirmed Zn-induced elevation of Wnt and Zn transporter genes in germ and supporting cells. Furthermore, the effects of Zn are dose-related in carp. Taken together, it seems wnt4a, wnt5b, and wnt8a play an important role in testis and exposure of Wnt inhibitor, canonical as well as non-canonical activators, and Zn confirmed that Zn regulates Wnt signaling vis-à-vis promoting spermatogenesis in the common carp.

Role of sox30 in regulating testicular steroidogenesis of common carp

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

Transcriptome characterization of BPG axis and expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine

Background

The clupeoid fishes are ecologically and commercially important fish species worldwide that exhibit a high level of population fluctuation, accompanied by alteration of reproductive traits. However, knowledge about their reproductive physiology in order to understand mechanisms underlying such population dynamics is limited. The endocrine system along with the brain–pituitary–gonadal (BPG) axis is critical for regulating reproduction. The aims of this study were to provide transcript data and genes related to the BPG axis, and to characterize the expression profiles of ovarian steroidogenesis-related genes in the Japanese sardine (Sardinops melanostictus, Clupeidae).

Results

RNA sequencing was performed using the sardine brain, pituitary, and gonad in both sexes. A total of 290,119 contigs were obtained and 115,173 non-redundant ORFs were annotated. The genes differentially expressed between ovary and testis were strongly associated with GO terms related to gamete production. The tissue-specific profile of the abundance of transcripts was characterized for the major regulators in the BPG axis, such as gonadotropin-releasing hormone, gonadotropin, and steroidogenic enzyme. By comparing between ovary and testis, out of eight different 17β-hydroxysteroid dehydrogenase (Hsd17b) genes identified, higher hsd17b7 expression was found in testis, whereas higher expression of hsd17b8, hsd17b10, hsd17b12a, and hsd17b12b was found in ovary. The cDNAs encoding key endocrine factors in the ovarian steroidogenic pathway were cloned, sequenced, and quantitatively assayed. In the pituitary, follicle-stimulating hormone beta peaked during vitellogenesis, while luteinizing hormone beta peaked at the completion of vitellogenesis. In the ovary, follicle-stimulating hormone receptor and luteinizing hormone receptor were upregulated from mid- to late phase of vitellogenesis. Furthermore, three steroidogenic enzyme genes (cyp11a1, cyp17a1, and cyp19a1a) gradually increased their expression during ovarian development, accompanying a rise in serum estradiol-17β, while 3β-hydroxysteroid dehydrogenase and steroidogenic acute regulatory protein did not change significantly.

Conclusions

This is the first report of deep RNA sequencing analysis of Japanese sardine, in which many key genes involved in the BPG axis were identified. Expression profiles of ovarian steroidogenesis-related genes provide a molecular basis of the physiological processes underlying ovarian development in the sardine. Our study will be a valuable resource for clarifying the molecular biology of clupeoid fishes.

Common carp pentraxin gene: Evidence for its role in ovarian differentiation and growth

Pentraxins (PTX), belong to an evolutionarily conserved family, containing a PTX protein domain, having role in acute immunological responses and fertility in higher vertebrates. However, information regarding the action of ptx on reproduction is extremely limited in fish. To study this, ptx cDNA was cloned for downstream analysis. Tissue distribution and ontogeny expression analysis indicated the prevalence of ptx in ovary. Varied phase-wise expression during carp ovarian cycle and elevated ptx expression after human chorionic gonadotropin induction, in vitro and in vivo, indicated probable regulation of gonadotropin. In situ hybridization and immunohistochemistry revealed the presence of ptx transcript and protein in the follicular layer of stage-III/IV oocytes indicating a role in ovarian growth. To assess the functional significance of ptx, transient silencing was performed using follicular primary cell culture, in vitro and in common carp, in vivo, through ovary-targeted injection of PEI-siRNA. Transient silencing of ptx-siRNA reduced the expression of various genes/factors related to oogenesis such as transcription factors, several steroidogenic enzymes, and esrs genes. These alterations in expression suggested a plausible role for ptx in ovarian steroidogenesis either, directly or indirectly, which is evident from the changes in the serum estradiol-17β (E₂) and 17α,20β-dihydroxyprogesterone levels. Furthermore, downregulation of aromatase activity was also noticed after transient silencing. Increased ptx expression after E₂ induced sex reversal to juvenile carp showed the correlative role of ptx during ovarian differentiation and development. Taken together, these findings suggest that ptx exerts an important role during ovarian growth, maturation and/or recrudescence of common carp.

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.