Stage-dependent differences in RNA composition and content affect the outcome of expression profiling in roach (Rutilus rutilus) ovary

High production of transfer RNAs identifies the presence of developing oocytes in ovaries and intersex testes of teleost fish

5S rRNA is highly transcribed in fish oocytes and this transcription levels can be used to identify the presence of oocytes in the intersex testes of fish exposed to xenoestrogens. Similar to 5S rRNA, tRNAs are transcribed by RNA polymerase III (Pol-III) in eukaryotes, so this study focuses in the analysis of the levels of expression of tRNAs in the gonads (ovaries and testes) of eight teleost species as a possible new oocyte molecular marker. Total RNA extracted from gonads of six commercial teleost species in the Biscay Bay, from the pollution sentinel species thicklip grey mullet (Chelon labrosus) known present intersex testes in response to xenoestrogens in Gernika estuary and from the laboratory model species Danio rerio were analysed through capillary electrophoresis. Bioanalyzer electropherograms were used to quantify the concentrations of tRNAs, 5S and 5.8S rRNA. All studied ovaries expressed significantly higher levels of tRNAs and 5S rRNA than testes. A tRNA to 5.8S rRNA index was calculated which differentiates ovaries from testes, and identifies some intersex testes in between testes and ovaries in mullets. The tRNA/5.8S ratio was highest in ovaries in previtellogenic stage, decreasing towards maturity. Thus, strong oocyte expression of tRNAs is an additional proof of high activity levels of Pol-III during early stages of oocyte development in teleost ovaries. Incidentally, we observed that miRNA concentrations were always higher in testes than ovaries. The indexing approach developed in the present study could have multiple applications in teleost reproduction research and in the development of early molecular markers of intersex condition.

Lipoprotein receptors in ovary of eel, Anguilla australis: molecular characterisation of putative vitellogenin receptors

Lipoprotein receptors, including low-density lipoprotein receptor (LDLr) relatives (Lrs) and LDLr-related proteins (Lrps), belong to the LDLr supergene family and participate in diverse physiological functions. In this study, novel sequences of lr and lrp genes expressed in the ovary of the short-finned eel, Anguilla australis, during early gonadal development are presented. The genes encoding the LDLr-like, Lrp1-like, Lrp1b-like, Lrp3, Lrp4-like, Lrp5-like, Lrp6, Lrp10, Lrp11, Lrp12-like, and Lr11-like proteins were found and identified by sequence and structure analysis, in addition to phylogenetic analysis. Genes encoding proteins previously implicated in follicle development and vitellogenin (Vtg) uptake in oviparous vertebrates were also identified, i.e. lr8 (including lr8 + and lr8- variants) and lrp13; their identification was reinforced by conserved synteny with orthologues in other teleost fish. Compared to other lr/lrp genes, the genes encoding Lr8 + , Lr8-, and Lrp13 were highly expressed in ovary during early development, decreasing as oocyte development advanced when induced by hypophysation. Furthermore, lr8 + , lr8-, and lrp13 were dominantly expressed in the ovary when compared with 17 other tissues. Finally, this study successfully detected the expression of both lr8 variants, which showed different expression patterns to those reported in other oviparous vertebrates and provided the first characterisation of Lrp13 in Anguilla sp. We propose that lr8 + , lr8-, and lrp13 encode putative Vtg receptors in anguillid eels.

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.

Quantitative transcriptomics, and lipidomics in evaluating ovarian developmental effects in Atlantic cod (Gadus morhua) caged at a capped marine waste disposal site

In the present study, a previously capped waste disposal site at Kollevåg (Norway) was selected to study the effects of contaminant leakage on biomarkers associated with Atlantic cod (Gadus morhua) reproductive endocrinology and development. Immature cod were caged for 6 weeks at 3 locations, selected to achieve a spatial gradient of contamination, and compared to a reference station. Quantitative transcriptomic, and lipidomic analysis was used to evaluate the effects of the potential complex contaminant mixture on ovarian developmental and endocrine physiology. The number of expressed transcripts, with 0.75 log2-fold differential expression or more, varied among stations and paralleled the severity of contamination. Particularly, significant bioaccumulation of ∑PCB-7, ∑DDTs and ∑PBDEs were observed at station 1, compared to the other station, including the reference station. Respectively 1416, 698 and 719 differentially expressed genes (DEGs), were observed at stations 1, 2 and 3, compared to the reference station, with transcripts belonging to steroid hormone synthesis pathway being significantly upregulation. Transcription factors such as esr2 and ahr2 were increased at all three stations, with highest fold-change at Station 1. MetaCore pathway maps identified affected pathways that are involved in ovarian physiology, where some unique pathways were significantly affected at each station. For the lipidomics, sphingolipid metabolism was particularly affected at station 1, and these effects paralleled the high contaminant burden at this station. Overall, our findings showed a novel and direct association between contaminant burden and ovarian toxicological and endocrine physiological responses in cod caged at the capped Kollevåg waste disposal site.

Expression patterns of kiss2 and gpr54-2 in Monopterus albus suggest these genes may play a role in sex reversal in fish

Due to its exceptionally small genome size and protogynous hermaphroditism, Monopterus albus has been proposed as a model for vertebrate sexual development. The Kiss/GPR54 system is a central regulator of sexual development in most vertebrates, but its role in sex reversal remains hypothetical. In contrast to mammals, fishes often possess more than one copy of the kiss and gpr54 genes. Our objectives were to identify all kiss/gpr54 genes in the genome of M. albus and to assess their involvement in sex reversal via their expression patterns (qPCR) in females, males, and intersex specimens. We identified only two genes: kiss2 and gpr54‐2. kiss2 expression was extremely high in the gonads of males, intermediate in females, and low in intersex; and reduced in all tissues of intersex. gpr54 expression was also extremely high in the gonads of males, high in intersex, but low in females. gpr54 expression in brain was high in all three sexes. In conclusion, (a) kiss1 has been functionally replaced in M. albus; (b) the functions of gpr54‐2 in brain are not sex‐specific; (c) kiss2 appears to undergo a ‘reset’ in the expression during the sex change; and (d) sex‐specific expression patterns in the gonads indicate that these two genes may play a role in sex reversal in fish.

Trans-omics analyses revealed differences in hormonal and nutritional status between wild and cultured female Japanese eel (Anguilla japonica)

Long-term stock decline in the Japanese eel (Anguilla japonica) is a serious issue. To reduce natural resource utilization in Japan, artificial hormonal induction of maturation and fertilization in the Japanese eel has been intensively studied. Recent experiment on feminized (by feeding a commercial diet containing estradiol-17β for first half year) cultured female eels have shown ovulation problems, which is seldom observed in captured wild female eels. Therefore, the aim of this study is to try to investigate causes of ovulation problem frequently seen in cultured female eels by comparative trans-omics analyses. The omics data showed low growth hormone and luteinizing hormone transcription levels in the brain and low sex hormone–binding globulin transcription levels in the liver of the cultured female eels. In addition, it was found that high accumulation of glucose-6-phosphate and, maltose in the liver of the cultured female eel. It was also found that docosahexaenoic (DHA) acid, eicosapentaenoic acid (EPA) and arachidonic acid (ARA) ratios in cultured female eels were quite different from wild female eels. The data suggested that ovulation problem in cultured female eels was possibly resulted from prolonged intake of a high-carbohydrate diet and/or suboptimal DHA/EPA/ARA ratios in a diet.

Ribosome RNA Profiling to Quantify Ovarian Development and Identify Sex in Fish

Terminologies of ovary development, by somewhat subjective describing and naming main changes of oocytes, have been criticized for confusing and inconsistency of terms and classifications, and the incurred consequences impede communication among researchers. In the present work, we developed regression between ovary development and three ribosome RNA (rRNA) indexes, namely 5S rRNA percent, 18S rRNA percent, and 5S–18S rRNA ratio, using close relationship between volume percent of primary growth stage oocytes or gonadosomatic index and rRNA content, demonstrating species-specific quantification of ovary development can be established in species with either synchronous and asynchronous oogenesis. This approach may be extended to any species with primary growth oocytes, e.g. anurans and reptiles, to predict maturity stages in females. We further confirmed that 5S rRNA percent and 5S/18S rRNA ratio can serve as markers to distinguish sexes unambiguously. A micro-invasive sampling method may be invented for non-lethal prediction of ovary development and sex because only a small amount of ovary sample (<50 mg) is needed for the approach established in the current work. Researchers who work with ovary RNA-seq in these taxa should realize that insufficient depletion of rRNA will probably lead to incorrect quantification of gene expression and inaccurate conclusions.

Molecular markers of oocyte differentiation in European eel during hormonally induced oogenesis

Reproduction in captivity is a key study issue in Anguilla anguilla as a possible solution for its dwindling population. Understanding the mechanisms controlling the production of ribosomal building blocks during artificially induced oocyte maturation could be particularly interesting. Transcription levels of ribosomal biogenesis associated genes could be used as markers to monitor oogenesis. Eels from the Albufera Lagoon were injected with carp pituitary extract for 15weeks and ovaries in previtellogenic (PV) stage (non-injected), in early-, mid-, late-vitellogenesis (EV, MV, LV), as well as in migratory nucleus stage (MN) were analysed. 5S rRNA and related genes were highly transcribed in ovaries with PV oocytes. As oocytes developed, transcriptional levels of genes related to 5S rRNA production (gtf3a), accumulation (gtf3a, 42sp43) and nucleocytoplasmic transport (rpl5, rpl11) and the 5S/18S rRNA index decreased (PV>EV>MV>LV>MN). On the contrary, 18S rRNA was at its highest at MN stage while ubtf1 in charge of activating RNA-polymerase I and synthesising 18S rRNA behaved as 5S related genes. Individuals that did not respond (NR) to the treatment showed 5S/18S index values similar to PV females, while studied genes showed EV/LV-like transcription levels. Therefore, NR females fail to express the largest rRNAs, which could thus be taken as markers of successful vitellogenesis progression. In conclusion, we have proved that the transcriptional dynamics of ribosomal genes provides useful tools to characterize induced ovarian development in European eels. In the future, such markers should be studied as putative indicators of response to hormonal treatments and of the quality of obtained eel oocytes.

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.

Abnormal ovarian DNA methylation programming during gonad maturation in wild contaminated fish

There is increasing evidence that pollutants may cause diseases via epigenetic modifications. Epigenetic mechanisms such as DNA methylation participate in the regulation of gene transcription. Surprisingly, epigenetics research is still limited in ecotoxicology. In this study, we investigated whether chronic exposure to contaminants experienced by wild female fish (Anguilla anguilla) throughout their juvenile phase can affect the DNA methylation status of their oocytes during gonad maturation. Thus, fish were sampled in two locations presenting a low or a high contamination level. Then, fish were transferred to the laboratory and artificially matured. Before hormonal treatment, the DNA methylation levels of the genes encoding for the aromatase and the receptor of the follicle stimulating hormone were higher in contaminated fish than in fish from the clean site. For the hormone receptor, this hypermethylation was positively correlated with the contamination level of fish and was associated with a decrease in its transcription level. In addition, whereas gonad growth was associated with an increase in DNA methylation in fish from the clean site, no changes were observed in contaminated fish in response to hormonal treatment. Finally, a higher gonad growth was observed in fish from the reference site in comparison to contaminated fish.

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.

Up-regulation of gonadotropin mRNA-expression at the onset of gametogenesis in the roach (Rutilus rutilus): evidence for an important role of brain-type aromatase (cyp19a1b) in the pituitary

The present study characterized changes in key parameters of reproduction in adult roach (Rutilus rutilus) from Lake Grosser Mueggelsee (Berlin, Germany) during natural gametogenesis. Fish of both sexes were sampled in monthly intervals between April and August in order to cover the onset of gametogenesis. Investigated parameters included gonad histology, plasma levels of 17β-oestradiol (E2), testosterone (T), 11-ketotestosterone (11-KT), and 17,20β-dihydroxy-4-pregnen-3-one (17,20β-P) as well as the expression of gonadotropin subunits in the pituitary. Furthermore, the mRNA-expression of brain-type aromatase (cyp19a1b), androgen receptor (ar), and estrogen receptor isoforms was studied at the pituitary level. The onset of gametogenesis - as indicated by follicles with cortical alveoli in females and first spermatogonia B in males - was observed in July, accompanied by a significant up-regulation of follicle-stimulating hormone β (fshβ) mRNA in the pituitary in both sexes. On the other hand, luteinizing hormone β (lhβ) mRNA increased later on in August. In males, the increase of fshβ mRNA in July coincided with a rise in plasma 11-KT concentrations. In females, E2 in plasma increased later, not until August, shortly before true vitellogenesis (late cortical alveoli stage). Expression of sex steroid receptors in the pituitary revealed only minor seasonal fluctuations. Most pronounced, ar mRNA displayed the highest level pre-spawning in both sexes. Interestingly, cyp19a1b mRNA-expression in the pituitary increased in parallel with fshβ already before any changes in plasma E2 or T occurred. These data suggest an important role of pituitary FSH and aromatase at the onset of gametogenesis in the roach.

Reference genes to quantify gene expression during oogenesis in a teleost fish

Understanding the molecular events involved in the acquisition of competence during oogenesis is a key step to determine the secret of 'high quality' eggs for aquaculture. Quantitative real time polymerase chain reaction (qPCR) is the technique of election to determine changes in transcript abundance in such studies, but choosing reference genes for normalization, in particular during oogenesis, remains a challenge. In the present study, transcription of 6 functionally distinct genes, β actin (ACTB), cathepsin D (CTSD), cathepsin Z (CTSZ), elongation factor 1 α (EEF1A), TATA binding protein (TBP) and tubulin A (TUBA1A) was assessed as normalizers of bone morphogenetic protein (BMP) and activin membrane-bound inhibitor (BAMBI) gene expression in mRNA from Mozambique tilapia oocytes during oogenesis. Reverse transcription was equally efficient and varies little in all samples. Most of the genes considered for reference were stable during early stages of oogenesis but variations were observed during vitellogenesis. A single gene and up to 3 genes were shown to be insufficient for reliable normalization throughout the whole oogenesis. The combination of the genes ACTB, CTSD, EEF1A and CTSZ as reference was found to minimize variation and has the most stable expression pattern between maturation stages.

Risks of seawater ozonation in recirculation aquaculture--effects of oxidative stress on animal welfare of juvenile turbot (Psetta maxima, L.)

Ozone is frequently used for water treatment and disinfection in recirculating aquaculture systems. However, due to the fragmentary data on chronic toxicity of ozone produced oxidants (OPO) and its safe concentrations, the daily application of ozone in aquaculture is challenging. To evaluate the chronic effects of sublethal OPO concentrations, juvenile turbot (Psetta maxima, L.) were exposed to OPO concentrations of 0.06, 0.10 and 0.15 mg/l for 21 days. Gills were analysed for histopathological alterations and mRNA expression of heat shock protein 70 (hsp70), hsp90 as well as glutathione-S-transferase (gst) were determined in the gills and the liver after 1d, 7d and 21 d. Histopathologic findings confirmed adverse effects at 0.10-0.15 mg/l, but these (necrosis, lamellar clubbing, hypertrophy, hyperplasia) could only be observed after an extended exposure (mostly 21 d), and were considered as irreversible tissue damage. Hsp70 expression in the gills was only significantly increased at the highest OPO concentration (0.15 mg/l) on 1d and 7d, and returned to basic levels until day 21. Hsp90 mRNA was already increased at 0.10mg/l after 1 and 7 days of exposure, and again was comparable to the control group on day 21. In contrast, elevated gst mRNA expression was only observed on day 7 at 0.10mg and 0.15 mg/l. Although similar trends were observed in the liver for all markers, differences were only significant in exceptional cases due to the high individual variation observed. Thus, mRNA expression in the gills rather than in the liver is recommended as a marker to characterize OPO-induced oxidative stress in turbot. It has to be noted that mRNA expression returned to basic levels on day 21 regardless the actual OPO concentration, suggesting a collapse of adaptive mechanisms as a possible explanation for the observed tissue damage.