A de novo transcriptome assembly approach elucidates the dynamics of ovarian maturation in the swordfish (Xiphias gladius)

Determination of Hg and Se in swordfish (Xiphias gladius) from Mediterranean Sea: Implications for nutritional recommendations during pregnancy and childhood

Swordfish (Xiphias gladius) are important marine predators facing potential mercury (Hg) contamination. This study analyzed Hg levels in swordfish muscle tissue from two Mediterranean Sea stocks. Italian specimens showed significantly higher Hg concentrations (1.4 ± 0.8 mg kg-1 ww) than Spanish ones (0.8 ± 0.5 mg kg-1 ww), with many samples exceeding EU's safety limit (1 mg kg-1 ww). Selenium (Se) content and Se:Hg ratios were examined, along with the maximum safe fish meals per month for vulnerable consumers. These results highlight the urgent need for monitoring and mitigation strategies to reduce health risks, especially for susceptible populations.

Global expression pattern of genes containing positively selected sites in European anchovy (Engraulis encrasicolus L.) may shed light on teleost reproduction

European anchovy is a multiple-spawning and highly fecundate pelagic fish with high economic and ecological significance. Although fecundity is influenced by nutrition, temperature and weight of spawners, high reproductive capacity is related to molecular processes in the ovary. The ovary is an essential and complex reproductive organ composed of various somatic and germ cells, which interact to facilitate the development of the ovary and functional oocytes. Revealing the ovarian transcriptome profile of highly fecundate fishes provides insights into oocyte production in teleosts. Here we use a comprehensive tissue-specific RNA sequencing which yielded 102.3 billion clean bases to analyze the transcriptional profiles of the ovary compared with other organs (liver, kidney, ovary, testis, fin, cauda and gill) and juvenile tissues of European anchovy. We conducted a comparative transcriptome and positive selection analysis of seven teleost species with varying fecundity rates to identify genes potentially involved in oogenesis and oocyte development. Of the 2,272 single copies of orthologous genes found, up to 535 genes were under positive selection in European anchovy and these genes are associated with a wide spectrum of cellular and molecular functions, with enrichments such as RNA methylation and modification, ribosome biogenesis, DNA repair, cell cycle processing and peptide/amide biosynthesis. Of the 535 positively selected genes, 55 were upregulated, and 45 were downregulated in the ovary, most of which were related to RNA and DNA transferase, developmental transcription factors, protein kinases and replication factors. Overall, our analysis of the transcriptome level in the ovarian tissue of a teleost will provide further insights into molecular processes and deepen our genetic understanding of egg production in highly fecund fish.

A Window of Vulnerability: Chronic Environmental Stress Does Not Impair Reproduction in the Swordfish Xiphias gladius

Xiphias gladius is an important fishing resource. The Mediterranean stock is affected by overfishing and is declining. In this light, the aim of this study was to evaluate the cross-talk among metabolism, stress response, immune system and reproduction in immature and mature females, coupling histological and transcriptomic approaches. The transcriptome of livers from 3 immature and 3 mature females was analyzed using the Artificial Intelligence RNA-Seq. For the histological analysis, ovary and liver samples were collected from 50 specimens caught during the reproductive season in the Mediterranean Sea. A total of 750 genes were differentially expressed between the livers. The gene ontologtabey analysis showed 91 upregulated and 161 downregulated biological process GO terms. Instead, the KEGG enrichment analysis revealed 15 enriched pathways. Furthermore, the binding occurring between estrogen receptors and aryl hydrocarbon receptor nuclear translocator, upregulated in mature females, could be liable for the inhibition of detoxification pathway. Indeed, at the histological level, mature females showed a higher density and number of melanomacrophage centers, biomarkers of stress. The present findings reveal the cross-talk among response to environmental stressors, metabolism and reproduction, highlighting that mature females invest a lot of energy in reproduction instead of immune response and detoxification.

Transcriptome analysis of flathead grey mullet (Mugil cephalus) ovarian development induced by recombinant gonadotropin hormones

Background: Treatment with recombinant gonadotropin hormones (rGths), follicle-stimulating hormone (rFsh) and luteinizing hormone (rLh), was shown to induce and complete vitellogenesis to finally obtain viable eggs and larvae in the flathead grey mullet (Mugil cephalus), a teleost arrested at early stages of gametogenesis in intensive captivity conditions. This study aimed to investigate the transcriptomic changes that occur in the ovary of females during the rGths-induced vitellogenesis. Methods: Ovarian samples were collected through biopsies from the same five females at four stages of ovarian development. RNASeq libraries were constructed for all stages studied, sequenced on an Illumina HiSeq4000, and a de novo transcriptome was constructed. Differentially expressed genes (DEGs) were identified between stages and the functional properties of DEGs were characterized by comparison with the gene ontology and Kyoto Encyclopedia. An enrichment analysis of molecular pathways was performed. Results: The de novo transcriptome comprised 287,089 transcripts after filtering. As vitellogenesis progressed, more genes were significantly upregulated than downregulated. The rFsh application induced ovarian development from previtellogenesis to early-to-mid-vitellogenesis with associated pathways enriched from upregulated DEGs related to ovarian steroidogenesis and reproductive development, cholesterol metabolism, ovarian growth and differentiation, lipid accumulation, and cell-to-cell adhesion pathways. The application of rFsh and rLh at early-to-mid-vitellogenesis induced the growth of oocytes to late-vitellogenesis and, with it, the enrichment of pathways from upregulated DEGs related to the production of energy, such as the lysosomes activity. The application of rLh at late-vitellogenesis induced the completion of vitellogenesis with the enrichment of pathways linked with the switch from vitellogenesis to oocyte maturation. Conclusion: The DEGs and enriched molecular pathways described during the induced vitellogenesis of flathead grey mullet with rGths were typical of natural oogenesis reported for other fish species. Present results add new knowledge to the rGths action to further raise the possibility of using rGths in species that present similar reproductive disorders in aquaculture, the aquarium industry as well as the conservation of endangered species.

Are Cell Junctions Implicated in the Regulation of Vitellogenin Uptake? Insights from an RNAseq-Based Study in Eel, Anguilla australis

At the onset of puberty, ovarian follicles become competent to incorporate large amounts of vitellogenin (Vtg). Using an RNAseq-based approach, transcriptomes from pre-vitellogenic (PV) and early vitellogenic (EV) ovaries from wild-caught eel, Anguilla australis, were compared to investigate the expression of specific genes encoding cell junction proteins that could be involved in regulating Vtg uptake. Partial support was found for the mechanical barrier hypothesis proposing that the access of Vtg to the oolemma is restricted by a tight junction (TJ) network within the granulosa cell layer, which changes between the PV and EV stage. Among 25 genes encoding TJ-constituting proteins, five were down-regulated and two were up-regulated. A chemical barrier hypothesis stating that gap junctions (GJs) are involved in modulating Vtg uptake was not supported, as only five GJs were found to be expressed in the ovary with no significant changes in expression between stages. Furthermore, the endocytic pathway was found to be up-regulated during the PV-EV transition. Finally, the study showed that gene expression patterns may help identify suitable candidates involved in the regulation of Vtg uptake, and provided novel sequence data for A. australis, including putative Vtg receptors corresponding to Lr8 and Lrp13 members of the low-density lipoprotein receptor family.

Probiotic Administration Mitigates Bisphenol A Reproductive Toxicity in Zebrafish

Although the use of bisphenol A (BPA) has been banned in a number of countries, its presence in the environment still creates health issues both for humans and wildlife. So far, BPA toxicity has been largely investigated on different biological processes, from reproduction to development, immune system, and metabolism. In zebrafish, Danio rerio, previous studies revealed the ability of environmentally relevant concentrations of this contaminant to significantly impair fertility via epigenetic modification. In addition, several studies demonstrated the ability of different probiotic strains to improve organism health. This study provides information on the role of the probiotic mixture SLAb51 to counteract adverse BPA effects on reproduction. A 28-day trial was set up with different experimental groups: BPA, exposed to 10 µg/L BPA; P, receiving a dietary supplementation of SLAb51 at a final concentration of 10⁹ CFU/g; BPA+P exposed to 10 µg/L BPA and receiving SLAb51 at a final concentration of 10⁹ CFU/g and a C group. Since oocyte growth and maturation represent key aspects for fertility in females, studies were performed on isolated class III (vitellogenic) and IV (in maturation) follicles and liver, with emphasis on the modulation of the different vitellogenin isoforms. In males, key signals regulating spermatogenesis were investigated. Results demonstrated that in fish exposed to the combination of BPA and probiotic, most of the transcripts were closer to C or P levels, supporting the hypothesis of SLAb51 to antagonize BPA toxicity. This study represents the first evidence related to the use of SLAb51 to improve reproduction and open new fields of investigation regarding its use to reduce endocrine disrupting compound impacts on health.

Alternative Splicing Dynamics of the Hypothalamus-Pituitary-Ovary Axis During Pubertal Transition in Gilts

The timing of puberty in mammals marks the point at which reproduction becomes possible. Abnormalities in the timing of puberty may exert a series of negative effects on subsequent health outcomes. Alternative splicing (AS) has not only emerged as a significant factor in the transcription of genes but it is also reported to play a role in the timing of puberty. However, to date, the changes and dynamics of AS during the onset of puberty is extremely seldom explored. In the present study, we used gilts as a research model to investigated the dynamics of AS and differentially expressed AS (DEAS) events within the hypothalamus-pituitary-ovary (HPO) axis across pre-, in-, and post-puberty. We detected 3,390, 6,098, and 9,085 DEAS events in the hypothalamus, pituitary, and ovary when compared across pre-, in-, and post-pubertal stages, respectively. Within the entire HPO axis, we also identified 22,889, 22,857, and 21,055 DEAS events in the pre-, in-, and post-pubertal stages, respectively. Further analysis revealed that the differentially spliced genes (DSGs) associated with staged DEAS events were likely to be enriched in the oxytocin signaling pathway, thyroid hormone signaling pathway, GnRH signaling pathway, and oocyte meiosis signaling pathway. The DSGs associated with DEAS events across the entire HPO axis were enriched in endocytosis signaling pathway, the MAPK signaling pathway, and the Rap1 signaling pathway. Moreover. the ASs of TAC1, TACR3, CYP19A1, ESR1, ESRRA, and FSHR were likely to regulate the functions of the certain HPO tissues during the onset of puberty. Collectively, the AS dynamics and DEAS events were comprehensively profiled in hypothalamus, pituitary, and ovary across the pre-, in-, and post-pubertal stages in pigs. These findings may enhance our knowledge of how puberty is regulated by AS and shed new light on the molecular mechanisms underlying the timing of puberty in mammals.

Opsins and gonadal circadian rhythm in the swordfish (Xiphias gladius) ovary: Their potential roles in puberty and reproductive seasonality

Circadian rhythm is well known to play a pivotal role in reproduction but the presence of a gonadal circadian rhythm is opening a lot of questions about a local regulation of reproduction. In the present study, we first set to identify the key genes driving circadian rhythmicity from the gonadal transcriptome of the swordfish, a commercially relevant species undergoing overfishing, and we then investigated whether their transcriptional activity was influenced by the maturation stage. Finally, we explored whether seasonality had the ability to modulate the expression of these genes. Interestingly, we identified a wide number of circadian rhythm related genes in the transcriptome of the swordfish gonad including, among the others, clock, bmal1, cry2 and per2, which have been found to be differentially expressed between sexually immature and mature individuals sampled during the breeding season. This differential modulation was also found for melatonin biosynthesis genes (mel1b, asmt and tph1) as well as opsin genes (opsin4, tmt opsin, parapinopsin, VA opsin, rho and sws), known to be the primary receptors of light stimuli. These differences were not observed between mature and immature individuals when considering the non-breeding season suggesting that, despite the molecular machinery of mature gonads is able to respond to signals driving ovarian maturation, these signals are not present hence highlighting the potential role of seasonality in modulating the gonadal circadian rhythm. These results confirm the presence of a gonadal circadian rhythm in the swordfish and open new interesting questions about its role in driving puberty onset.

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.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.