5S rRNA and accompanying proteins in gonads: powerful markers to identify sex and reproductive endocrine disruption in fish

Method for the molecular and quantitative identification of oocytes and their developmental stage in teleost fish gonads

Fish display diverse reproductive strategies and their gametogenesis is influenced by numerous genetic, physiological and environmental factors. The analysis of 5S rRNA expression levels in gonads has been proposed as useful method for the molecular identification of the presence of oocytes in fish tissues. The present method provides an easy and unbiased approach to analyse the expression of tRNAs and 5S rRNA in teleost gonads and stablish the presence and developmental stage of oocytes. Total RNA extracted from gonads is analysed through capillary electrophoresis in a Bioanalyzer 2100 (Agilent Technologies) using Small RNA Assays. Electropherograms allow quantifying the concentrations of tRNAs, 5S rRNA and 5.8S rRNA per sample and calculate their tRNA/5.8S rRNA and 5S/5.8S rRNA indices. Both indices clearly differentiate ovaries from testes and can be used to identify testes that present oocytes due to exposure to environmental xenoestrogens. The tRNA/5.8S and 5S/5.8S indices show the highest values in ovaries in previtellogenic stage, values decreasing as they advance towards maturity.•Detailed molecular method to sex fish and quantitatively identify the maturity stage of females.•tRNA levels in gonads can help in the study of teleost reproduction (female fecundity assessment, molecular gonad sexing) and environmental health assessment.

Otolith shape and elemental signatures provide insights into the connectivity of euryhaline Chelon labrosus inhabiting two close estuaries with different burdens of xenoestrogens in the Southern Bay of Biscay

Intersex gonads have been observed in thicklip grey mullet Chelon labrosus, inhabiting estuaries with high burdens of xenoestrogens in the Southeast Bay of Biscay, but knowledge of population connectivity among estuaries is lacking for this euryhaline fish species. This study investigates the population structure of C. labrosus using otolith shape and elemental signatures of 60 adults (overall length ∼ 38 cm) from two estuaries 21 nautic miles apart, one with a high incidence of intersex condition (Gernika), and the other one pristine (Plentzia). Otolith shape analyses were performed using elliptical Fourier descriptors, while elemental signatures of whole sagittae were obtained by inductively coupled plasma mass spectrophotometry. Univariate and multivariate statistics were applied to determine if otolith signatures show patterns of homogeneity between estuaries. The data indicated significant differences in the otolith shape and elemental composition between mullets of Gernika and Plentzia. Elemental differences were mainly driven by Sr, Li (both higher in Plentzia) and Ba (higher in Gernika). The 98% re-classification success rate obtained from stepwise linear discriminant function analysis suggests that Gernika and Plentzia individuals form separated population units. The limited connectivity between these two close estuaries would indicate a different life history of exposure to chemicals, which might explain the high prevalence of intersex condition in Gernika and its absence in Plenztia.

Lack of genetic structure in euryhaline Chelon labrosus from the estuaries under anthropic pressure in the Southern Bay of Biscay to the coastal waters of the Mediterranean Sea

Over the last decade, xenoestrogenic effects have been reported in populations of thicklip grey mullet Chelon labrosus from contaminated estuaries in the Bay of Biscay, resulting in intersex condition. To understand the level of gene flow in individuals of different Basque estuaries microsatellite markers were used to evaluate the population structure and connectivity of C. labrosus from estuaries of the Basque coast. 46 microsatellites were tested and 10 validated for the analysis of 204 individuals collected from 5 selected Basque estuaries and 2 outgroups in the Bay of Cadiz and Thermaic Gulf. The polymorphic microsatellites revealed 74 total alleles, 2-19 alleles per locus. The mean observed heterozygosity (0.49 ± 0.02) was lower than the expected one (0.53 ± 0.01). There was no evidence of genetic differentiation (F_ST = 0.0098, P = 0.0000) among individuals or sites. Bayesian clustering analysis revealed a single population in all sampled locations. The results of this study indicate widespread genetic homogeneity and panmixia of C. labrosus across the current sampling areas spanning the Atlantic and Mediterranean basins. The hypothesis of panmixia could therefore be well supported so individuals inhabiting estuaries with high prevalence of intersex condition should be considered as members of the same single genetic group as those inhabiting adjacent estuaries without incidence of xenoestrogenicity.

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.

Apoptosis and autophagy-related gene transcription during ovarian follicular atresia in European hake (Merluccius merluccius)

Follicular atresia is an energy-saving oocyte resorption process that can allow the survival of female fish when environmental conditions are unfavourable and at the expense of fecundity. This study investigated the transcription levels of apoptosis and autophagy-related genes during atresia in the European hake that can show episodes of increased follicular atresia throughout the reproductive cycle. 169 female individuals were collected from the Bay of Biscay, and the ovaries were analysed using histological and molecular methods. Different levels of atresia were histologically detected in 73.7% of the ovaries analysed and the TUNEL assay identified apoptotic nuclei in follicles from both previtellogenic and vitellogenic stages. Transcripts of beclin-1 and ptenb were up-regulated in the ovaries containing atretic follicles, whereas p53, caspase-3, cathepsin D and dapk1 were up-regulated only in ovaries presenting vitellogenic atretic follicles. Our results indicate different implications of apoptotic vs autophagic processes leading to atresia during oocyte development, vitellogenesis being the moment of maximal apoptotic and autophagic activity in atretic hakes. The analysed genes could provide early warning biomarkers to identify follicular atresia in fish and evaluate fecundity in fish stocks.

A case study: temporal trends of environmental stressors and reproductive health of smallmouth bass (Micropterus dolomieu) from a site in the Potomac River Watershed, Maryland, USA

Decades of poor reproductive success and young-of-the-year survival, combined with adult mortality events, have led to a decline in the smallmouth bass (SMB; Micropterus dolomieu) population in sections of the Potomac River. Previous studies have identified numerous biologic and environmental stressors associated with negative effects on SMB health. To better understand the impact of these stressors, this study was conducted at the confluence of Antietam Creek and the Potomac River from 2013 to 2019 to identify temporal changes associated with SMB reproductive health. Surface water samples were collected and analyzed for over 300 organic contaminants, including pesticides, phytoestrogens, pharmaceuticals, hormones and total estrogenicity (E2Eq). Adult SMB were collected and sampled for multiple endpoints, including gene transcripts associated with reproduction (molecular), histopathology (cellular), and organosomatic indices (tissue). In males, biomarkers of estrogenic endocrine disruption, including testicular oocytes (TO) and plasma vitellogenin (Vtg) were assessed. Numerous agriculture-related contaminants or land use patterns were associated with gene transcript abundance in both male and female SMB. Positive associations between pesticides in the immediate catchment with TO severity and E2Eq with plasma Vtg in males were identified. In males, the prevalence of TO and detectable levels of plasma Vtg, liver vitellogenin transcripts (vtg) and testis vtg were high throughout the study. Peaks of complex mixtures of numerous contaminants occurred during the spring/early summer when spawning and early development occurs and to a lesser extent in fall/winter during recrudescence. Management practices to reduce exposure during these critical and sensitive periods may enhance reproductive health of these economically important sportfishes.

First discovered, long out of sight, finally visible: ribosomal DNA

With the advent of long-read sequencing, previously unresolvable genomic elements are being revisited in an effort to generate fully complete reference genomes. One such element is ribosomal DNA (rDNA), the highly conserved genomic region that encodes rRNAs. Genomic structure and content of the rDNA are variable in both prokarya and eukarya, posing interesting questions about the biology of rDNA. Here, we consider the types of variation observed in rDNA - including locus structure and number, copy number, and sequence variation - and their known phenotypic consequences. With recent advances in long-read sequencing technology, incorporating the full rDNA sequence into reference genomes is within reach. This knowledge will have important implications for understanding rDNA biology within the context of cell physiology and whole-organism phenotypes.

A new experimental model for the investigation of sequential hermaphroditism

The stunning sexual transformation commonly triggered by age, size or social context in some fishes is one of the best examples of phenotypic plasticity thus far described. To date our understanding of this process is dominated by studies on a handful of subtropical and tropical teleosts, often in wild settings. Here we have established the protogynous New Zealand spotty wrasse, Notolabrus celidotus, as a temperate model for the experimental investigation of sex change. Captive fish were induced to change sex using aromatase inhibition or manipulation of social groups. Complete female-to-male transition occurred over 60 days in both cases and time-series sampling was used to quantify changes in hormone production, gene expression and gonadal cellular anatomy. Early-stage decreases in plasma 17β-estradiol (E2) concentrations or gonadal aromatase (cyp19a1a) expression were not detected in spotty wrasse, despite these being commonly associated with the onset of sex change in subtropical and tropical protogynous (female-to-male) hermaphrodites. In contrast, expression of the masculinising factor amh (anti-Müllerian hormone) increased during early sex change, implying a potential role as a proximate trigger for masculinisation. Collectively, these data provide a foundation for the spotty wrasse as a temperate teleost model to study sex change and cell fate in vertebrates.

Non-coding RNA Expression Patterns of Two Different Teleost Gonad Maturation Stages

Non-coding RNAs (ncRNAs) are involved in several different regulatory pathways including reproduction. In teleost fish, efficacious reproduction is heavily dependent on the completion of the reproductive cycle. The presence of ncRNA, however, and their expression dynamics and putative regulatory role in mature and immature gonads have not yet been extensively explored. Therefore, the abundance of ncRNAs in mature and immature female sharpsnout seabream (Diplodus puntazzo) was investigated. The sharpsnout seabream is a rudimentary hermaphrodite which, in captivity, displays dysfunctions in the gonad maturation process. Our analyses revealed a gonad specific read length distribution with two main peaks representing miRNAs (21-26 nt) and PIWI RNA (27-34 nt). Besides, distinct expression patterns for several ncRNA biotypes including microRNAs (miRNAs), PIWI RNAs (piRNAs), and ribosomal RNAs (rRNAs) were detected. Identified miRNA accounted to 938, corresponding to ~ 13% of obtained transcripts. Among the differential expressed ncRNAs, 10 (~ 7%) were annotated as miRNA, out of which 2 were found in higher abundance in immature gonads (miR-125c and miR-24) and 8 (miR-451, miR-7a, miR-122-1, miR190a, miR129, ENSGACT00000029608, ENSGACT00000029489, and ENSGACT00000029667) were found to be higher expressed in mature gonads. Putative miRNA targets, including long non-coding RNAs (lncRNAs) and genes, are proposed. Target genes are involved in several processes of fish oocyte development, such as steroidogenesis, proteolysis, and apoptosis, and may explain hormone regulation. This study demonstrates a gonad maturation biased ncRNA profile which in turn may support the role of ncRNAs in ovarian physiology and reproductive performance of fish, stressing the specific function of each RNA biotype in oocyte development.

The effects of 11-ketotestosterone implants on transcript levels of gonadotropin receptors, and foxl2 and dmrt1 genes in the Previtellogenic ovary of cultured beluga (Huso huso)

The in vivo effect of 11-ketotestosterone (11KT) on transcript levels of the gonadotropin receptors (fshr and lhr) and sex differentiation-related genes (dmrt1 and foxl2) was examined in the ovaries of immature female beluga. For this purpose, six fish were treated with implants containing 2.5 mg 11KT and a placebo group of six females of the same age and gametogenic stage were given a blank implant. The implants were intraperitoneally inserted into 4-year-old females at the previtellogenic stage (mean body weight 5580 ± 165 g) and maintained under culture conditions for 8 weeks. Ovary samples for gene expression analysis of lhr, fshr, dmrt1 and foxl2 were collected by biopsy at 3 and 8 weeks post implantation. Diameters of oocytes increased in response to 11KT treatment, both at 3 and at 8 weeks post implantation, but no obvious changes were evident in cytology. Three weeks of 11KT treatment did not affect target gene expression, but a tendency for a time-dependent decrease of lhr and dmrt1 mRNA levels was observed in both treatment and placebo groups. By 8 weeks of treatment, however, 11KT implants provoked the upregulation of fshr and foxl2 transcript levels. Furthermore, lhr and dmrt1 transcript abundances recovered by 8 weeks of exposure in both blank- and 11KT-implanted beluga. These results suggest that 11KT, either directly or indirectly, may affect gametogenesis and regulate some key components of the reproductive axis in female beluga.

Transcriptomic analysis of female and male gonads in juvenile snakeskin gourami (Trichopodus pectoralis)

The snakeskin gourami (Trichopodus pectoralis) exhibits sexual dimorphism, particularly in body size. Since the snakeskin gourami is usually marketed during sexual maturation, the sexual size dimorphism has become an economically important trait. Sex-biased gene expression plays a key role in phenotypic sexual dimorphism. Therefore, using high-throughput RNA sequencing (RNA-seq) technology, we aimed to explore the differentially expressed genes (DEGs) in ovary and testis during sex differentiation in juvenile snakeskin gourami. Our results revealed a number of DEGs were demonstrated to be overexpressed in ovary (11,625 unigenes) and testis (16,120 unigenes), and the top 10 female-biased (rdh7, dnajc25, ap1s3, zp4, polb, parp12, trim39, gucy2g, rtbs, and fdxr) and male-biased (vamp3, nbl1, dnah2, ccdc11, nr2e3, spats1, pih1d2, tekt3, fbxo36, and mybl2) DEGs were suggested to be mainly associated with ovary and testis differentiation, respectively. Additionally, using real-time reverse transcription polymerase chain reaction (qRT-PCR), validation of the differential expression of 21 genes that were previously shown to be related to gonad development was performed (ar, bHLH, cyp19a1, daz, dead-end, esrb, esrrg, gnrhr, gpa, gsg1l, hsd17B, mospd1, nanos-1, nanos-2, p53, piwi-1, piwi-2, rerg, rps6ka, tgf-beta, and VgR). The results showed a significantly positive correlation (0.84; P < 0.001) between the results of RNA-seq and qRT-PCR. Therefore, RNA-seq analysis in our study identified global genes that were associated with ovary and testis differentiation in the juvenile phase of the snakeskin gourami. Our findings provide valuable transcriptomic bioinformation for further investigation of reproductive biology and applications of sex manipulation.

Transcription pattern of reproduction relevant genes along the brain-pituitary-gonad axis of female, male and intersex thicklip grey mullets, Chelon labrosus, from a polluted harbor

The reproductive cycle of teleost fishes is regulated by the brain-pituitary-gonad (BPG) axis. The transcription profile of genes involved in the reproduction signalling in the BPG-axis differs in females and males during the gametogenic cycle. Impacts of endocrine disrupting chemicals on these signalling pathways in fish are known, but the participation of the BPG-axis in the development of the intersex condition is not well understood. Intersex thicklip grey mullets (Chelon labrosus) have been identified in several estuaries from the SE Bay of Biscay, revealing the presence of feminizing contaminants in the area. In previous studies, transcription patterns of genes related with steroidogenesis and gamete growth have been shown to differ among female, male and intersex mullets. However, many components of the reproduction control have not been studied yet. The aim of this study was to assess the transcription levels of target BPG-axis genes in female, male and intersex mullets captured in the polluted harbour of Pasaia, during their gametogenic cycle. After histologically examining the gonads, the transcription levels of previously sequenced target genes were measured by qPCR: kiss2, gpr54 and gnrh1 in brain, fshβ and lhβ in pituitary and fshr and lhr in gonads. In both females and males, brain genes were most transcribed in early gametogenesis, proving their relevance in the onset of both oogenesis and spermatogenesis. Pituitary gonadotropins in females showed upregulation as oogenesis progressed, reaching the highest transcription levels at vitellogenic stage, while in males transcript levels were constant during spermatogenesis. Transcription levels of gonadotropin receptors showed different patterns in ovaries and testes, suggesting differing function in relation to gametogenesis and maturation. Intersex mullets showed transcription levels of brain target genes similar to those observed in females at cortical alveoli stage and to those in mid spermatogenic males. In intersex testes the transcription pattern of gonadotropin receptor fshr was downregulated in comparison to non-intersex testes.

Duplication and subfunctionalisation of the general transcription factor IIIA (gtf3a) gene in teleost genomes, with ovarian specific transcription of gtf3ab

Fish oogenesis is characterised by a massive growth of oocytes each reproductive season. This growth requires the stockpiling of certain molecules, such as ribosomal RNAs to assist the rapid ribosomal assembly and protein synthesis required to allow developmental processes in the newly formed embryo. Massive 5S rRNA expression in oocytes, facilitated by transcription factor 3A (Gtf3a), serves as marker of intersex condition in fish exposed to xenoestrogens. Our present work on Gtf3a gene evolution has been analysed in silico in teleost genomes and functionally in the case of the zebrafish Danio rerio. Synteny-analysis of fish genomes has allowed the identification of two gtf3a paralog genes, probably emerged from the teleost specific genome duplication event. Functional analyses demonstrated that gtf3ab has evolved as a gene specially transcribed in oocytes as observed in Danio rerio, and also in Oreochromis niloticus. Instead, gtf3aa was observed to be ubiquitously expressed. In addition, in zebrafish embryos gtf3aa transcription began with the activation of the zygotic genome (~8 hpf), while gtf3ab transcription began only at the onset of oogenesis. Under exposure to 100 ng/L 17β-estradiol, fully feminised 61 dpf zebrafish showed transcription of ovarian gtf3ab, while masculinised (100 ng/L 17α-methyltestosterone treated) zebrafish only transcribed gtf3aa. Sex related transcription of gtf3ab coincided with that of cyp19a1a being opposite to that of amh and dmrt1. Such sex dimorphic pattern of gtf3ab transcription was not observed earlier in larvae that had not yet shown any signs of gonad formation after 26 days of oestradiol exposure. Thus, gtf3ab transcription is a consequence of oocyte differentiation and not a direct result of estrogen exposure, and could constitute a useful marker of gonad feminisation and intersex condition.

POLR1B and neural crest cell anomalies in Treacher Collins syndrome type 4

Purpose

Treacher Collins syndrome (TCS) is a rare autosomal dominant mandibulofacial dysostosis, with a prevalence of 0.2–1/10,000. Features include bilateral and symmetrical malar and mandibular hypoplasia and facial abnormalities due to abnormal neural crest cell (NCC) migration and differentiation. To date, three genes have been identified: TCOF1, POLR1C, and POLR1D. Despite a large number of patients with a molecular diagnosis, some remain without a known genetic anomaly.

Methods

We performed exome sequencing for four individuals with TCS but who were negative for pathogenic variants in the known causative genes. The effect of the pathogenic variants was investigated in zebrafish.

Results

We identified three novel pathogenic variants in POLR1B. Knockdown of polr1b in zebrafish induced an abnormal craniofacial phenotype mimicking TCS that was associated with altered ribosomal gene expression, massive p53-associated cellular apoptosis in the neuroepithelium, and reduced number of NCC derivatives.

Conclusion

Pathogenic variants in the RNA polymerase I subunit POLR1B might induce massive p53-dependent apoptosis in a restricted neuroepithelium area, altering NCC migration and causing cranioskeletal malformations. We identify POLR1B as a new causative gene responsible for a novel TCS syndrome (TCS4) and establish a novel experimental model in zebrafish to study POLR1B-related TCS.

Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change

Bluehead wrasses undergo dramatic, socially cued female-to-male sex change. We apply transcriptomic and methylome approaches in this wild coral reef fish to identify the primary trigger and subsequent molecular cascade of gonadal metamorphosis. Our data suggest that the environmental stimulus is exerted via the stress axis and that repression of the aromatase gene (encoding the enzyme converting androgens to estrogens) triggers a cascaded collapse of feminizing gene expression and identifies notable sex-specific gene neofunctionalization. Furthermore, sex change involves distinct epigenetic reprogramming and an intermediate state with altered epigenetic machinery expression akin to the early developmental cells of mammals. These findings reveal at a molecular level how a normally committed developmental process remains plastic and is reversed to completely alter organ structures.

Toxicology tailored low density oligonucleotide microarray for the thicklip grey mullets (Chelon labrosus): Biomarker gene transcription profile after caging in a polluted harbour

In aquatic organisms inhabiting polluted waters genes are activated to build an adaptive/compensatory defence against the possible effects of pollutants. Such responses can be used as biomarkers of exposure to chemical compounds, outlining the molecular mechanisms activated under specific pollution scenarios. With the aim of exploiting such approach in environmental health assessment, toxicologically relevant gene fragments were sequenced in the thicklip grey mullet (Chelon labrosus) and a toxicologically tailored low-density (160 genes) oligonucleotide microarray was customised. The tool was validated comparing organ/sex specific gene expression profiles and characterising responses under laboratory exposure to model chemicals. Finally, juvenile mullets were caged in a polluted harbour and hepatic gene expression profiles analysed after 5 and 21 days of deployment. Cages were deployed in the inner (IH) and outer (OH) Pasaia harbour, Bay of Biscay. Mussels (Mytilus galloprovincialis) were also caged as biological matrix for chemical bioaccumulation analysis and stress biomarkers measurements. Slightly higher concentrations of chemicals (metals, tributyltin, PAHs, phthalates) were quantified in IH than in OH, fish bile metabolites also revealing higher availability of PAHs in IH. Lysosome membrane stability in mussels was reduced, indicating stress condition in both sites. The developed microarray discriminated mullets showing distinctive expression profiles depending on site and deployment time. Genes related to immune and hypoxia responses were regulated comparing IH and OH at day 5. Phase I and II biotransformation genes, such as cyp2, cyp3 and ugt, were up-regulated in IH, together with the aryl hydrocarbon receptor 2 (ahr2) and the ahr repressor. Similarly, TBT-binding proteins and genes involved in lipid metabolism (pparγ, cyp7) were up-regulated with deployment time. Even if nowadays higher throughput approaches for gene expression analyses are available, the developed mullet tool constitutes a comprehensive tool to assess molecular responses of mullets exposed to pollutants, although it remains to be explored whether it can be applied to assess pollutant exposure in active pollution monitorings and in environmental health assessment.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

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.

Alteration in molecular markers of oocyte development and intersex condition in mullets impacted by wastewater treatment plant effluents

Wastewater Treatment Plant (WWTP) discharges are an important source of endocrine disrupting chemicals (EDCs) into the aquatic environment. Fish populations inhabiting downstream of WWTP effluents show alterations in gonad and gamete development such as intersex condition, together with xenoestrogenic effects such as vitellogenin up-regulation. However, the molecular mechanisms participating in the development of intersex condition in fish are not elucidated. The aim of this study was to assess the impact of two WWTPs effluents (Gernika and Bilbao-Galindo situated in the South East Bay of Biscay) with different contaminant loads, in thicklip grey mullet (Chelon labrosus) populations inhabiting downstream, examining the presence and severity of intersex condition, during two seasons. Molecular markers of xenoestrogenicity and oocyte differentiation and development (vtgAa, cyp19a1a, cyp19a1b, cyp11b, foxl2, dmrt1 and gtf3a) were also studied. Intersex mullets were identified downstream of both WWTPs and vtgAa was upregulated in intersex and non intersex males. Sex dependent differential transcription levels of target genes were detected in mullets from Galindo. However, no such pattern was observed in mullets from Gernika, suggesting an attenuating effect over studied genes caused by a higher presence of EDCs in this site, as indicated by the elevated prevalence of intersex mullets in this population. In conclusion, no direct association between xenoestrogenic responses and intersex condition was established. Mullets from Gernika showed signs of severe EDC exposure compared to those from Galindo, as demonstrated by the higher prevalence of intersex males and the reduction in transcription profile differences between sexes of gametogenic gene markers.

De novo transcriptome sequencing and analysis of male, pseudo-male and female yellow perch, Perca flavescens

Transcriptome sequencing could facilitate discovery of sex-biased genes, biological pathways and molecular markers, which could help clarify the molecular mechanism of sex determination and sexual dimorphism, and assist with selective breeding in aquaculture. Yellow perch has unique gonad system and sexual dimorphism and is an alternative model to study mechanism of sex determination, sexual dimorphism and sexual selection. In this study, we performed the de novo assembly of yellow perch gonads and muscle transcriptomes by high throughput Illumina sequencing. A total of 212,180 contigs were obtained, ranging from 127 to 64,876 bp, and N50 of 1,066 bp. The assembly RNA-Seq contigs (≥200bp) were then used for subsequent analyses, including annotation, pathway analysis, and microsatellites discovery. No female- and pseudo-male-biased genes were involved in any pathways while male-biased genes were involved in 29 pathways, and neuroactive ligand receptor interaction and enzyme of trypsin (enzyme code, EC: 3.4.21.4) was highly involved. Pyruvate kinase (enzyme code, EC: 2.7.1.40), which plays important roles in cell proliferation, was highly expressed in muscles. In addition, a total of 183,939 SNPs, 11,286 InDels and 41,479 microsatellites were identified. This study is the first report on transcriptome information in Percids, and provides rich resources for conducting further studies on understanding the molecular basis of sex determinations, sexual dimorphism, and sexual selection in fish, and for population studies and marker-assisted selection in Percids.

Photoperiod modulates the reproductive axis of European sea bass through regulation of kiss1 and gnrh2 neuronal expression

The onset of puberty is characterized by activation of the brain-pituitary-gonad axis. However, the molecular and endocrine mechanism involved in the process of puberty and the influence of environmental conditions, such as photoperiod signalling, are not well understood in fish. In this study, 1-year-old male European sea bass (Dicentrarchus labrax) were exposed to photoperiod manipulation in combination with size-sorting. Two treatment groups, a puberty accelerating photoperiod (AP) group and a continuous light (LL) group, were studied from August to February. Our results indicate that AP and LL are able to entrain the rhythms of both kiss1 and gnrh2 mRNA levels in the brain, while kiss2 and gnrh1 mRNA expression does not seem to be directly affected by the photoperiod, at least during testicular growth. It is likely that AP and LL photoperiod regimes affected both plasma Fsh and 11-KT profiles, which might explain, respectively, the phase shift and reduction of testes maturation seen under these conditions. We therefore hypothesize that the unbalanced production of this androgen regulated by circulating Fsh might be limiting the stimulation of germ cell proliferation in European sea bass males. In summary, our study establishes that photoperiod modulates the expression of kiss1 and gnrh2 in the forebrain-midbrain, which may be involved in the translation of the light stimulus to activate the reproductive axis.

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.

New insights into the factors mediating the onset of puberty in sea bass

In populations of 1-year-old male European sea bass (Dicentrarchus labrax), only large males are able to acquire for the first time a functional competence of their reproductive axis; in other words, to attain puberty. To examine the causes and mechanisms involved in the onset of puberty in this species, a size sorting sampling was carried out to obtain two experimental groups of small and large male fish exhibiting different growth rates. As expected, only large fish reached full spermiogenesis (stage V of testicular development) by the end of the experiment. Our study suggests that fish size is a permissive condition to ensure full effectiveness of the hormonal (Gnrh1, gonadotropins and sexual steroids) actions. Thus, though small fish had endocrine profiles similar to those of large fish, their amplitude was much lower, and was most likely the reason why functional competence of the reproductive axis was not achieved. Moreover, this work provides evidence of the involvement of kisspeptin and Gnrh1 systems in the onset of puberty in a marine teleost fish. It also indicates that very likely kisspeptin and Gnrh1 may regulate gonadotropins and sex steroids at specific stages of testicular development.

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.

Development and validation of a mixed-tissue oligonucleotide DNA microarray for Atlantic bluefin tuna, Thunnus thynnus (Linnaeus, 1758)

BACKGROUND

The largest of the tuna species, Atlantic bluefin tuna (Thunnus thynnus), inhabits the North Atlantic Ocean and the Mediterranean Sea and is considered to be an endangered species, largely a consequence of overfishing. T. thynnus aquaculture, referred to as fattening or farming, is a capture based activity dependent on yearly renewal from the wild. Thus, the development of aquaculture practices independent of wild resources can provide an important contribution towards ensuring security and sustainability of this species in the longer-term. The development of such practices is today greatly assisted by large scale transcriptomic studies.

RESULTS

We have used pyrosequencing technology to sequence a mixed-tissue normalised cDNA library, derived from adult T. thynnus. A total of 976,904 raw sequence reads were assembled into 33,105 unique transcripts having a mean length of 893 bases and an N50 of 870. Of these, 33.4% showed similarity to known proteins or gene transcripts and 86.6% of them were matched to the congeneric Pacific bluefin tuna (Thunnus orientalis) genome, compared to 70.3% for the more distantly related Nile tilapia (Oreochromis niloticus) genome. Transcript sequences were used to develop a novel 15 K Agilent oligonucleotide DNA microarray for T. thynnus and comparative tissue gene expression profiles were inferred for gill, heart, liver, ovaries and testes. Functional contrasts were strongest between gills and ovaries. Gills were particularly associated with immune system, signal transduction and cell communication, while ovaries displayed signatures of glycan biosynthesis, nucleotide metabolism, transcription, translation, replication and repair.

CONCLUSIONS

Sequence data generated from a novel mixed-tissue T. thynnus cDNA library provide an important transcriptomic resource that can be further employed for study of various aspects of T. thynnus ecology and genomics, with strong applications in aquaculture. Tissue-specific gene expression profiles inferred through the use of novel oligo-microarray can serve in the design of new and more focused transcriptomic studies for future research of tuna physiology and assessment of the welfare in a production environment.

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.

Gonadal transcriptome analysis of male and female olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is an important commercially cultured marine flatfish in China, Korea, and Japan, of which female grows faster than male. In order to explore the molecular mechanism of flounder sex determination and development, we used RNA-seq technology to investigate transcriptomes of flounder gonads. This produced 22,253,217 and 19,777,841 qualified reads from ovary and testes, which were jointly assembled into 97,233 contigs. Among them, 23,223 contigs were mapped to known genes, of which 2,193 were predicted to be differentially expressed in ovary and 887 in testes. According to annotation information, several sex-related biological pathways including ovarian steroidogenesis and estrogen signaling pathways were firstly found in flounder. The dimorphic expression of overall sex-related genes provides further insights into sex determination and gonadal development. Our study also provides an archive for further studies of molecular mechanism of fish sex determination.

Intersex condition and molecular markers of endocrine disruption in relation with burdens of emerging pollutants in thicklip grey mullets (Chelon labrosus) from Basque estuaries (South-East Bay of Biscay)

Endocrine disrupting chemicals (EDCs) interfere with the functioning of the endocrine system, causing reproductive and developmental disturbances in aquatic wildlife. Appearance of intersex gonads and elevated plasma levels of vitellogenin in male fish are well known biomarkers of exposure to xenoestrogenic EDCs. In the present study, intersex condition and transcription levels of vtg and cyp19a1b were assessed in five thicklip grey mullet populations from the Basque coast (Bay of Biscay). Levels of EDCs (estrogenic hormones, polycyclic musks, bisphenol-A, phthalates, alkylphenols and pesticides) were determined in water and fish bile. Intersex gonads were observed in three out of five mullet populations. Vtg and cyp19a1b were up-regulated in mullet populations with relatively higher EDCs load. Phthalates and pesticides were the most abundant EDCs in bile, followed by alkylphenols, musks, bisphenol-A and estrogenic hormones. Statistically significant correlations were found between concentrations of individual and total EDCs in bile and water samples and transcription levels of vtg and cyp19a1b.

Intersex in teleost fish: are we distinguishing endocrine disruption from natural phenomena?

Intersex is defined as the simultaneous presence of male and female gonadal tissue in a gonochoristic (fixed-sex) species. The intersex condition has been documented in both wild and laboratory animals, including fish, amphibians, and reptiles. In aquatic animals, intersex is often viewed as a signature effect of exposure to endocrine disrupting compounds. At least 37 fish species from 17 families have been identified with intersex gonads in 54 field survey studies. However, reports of the occurrence of intersex at reference sites have led to speculation that a baseline level of intersex is "normal". The objective of this critical review was to assess factors potentially associated with baseline levels of intersex in fish and to examine the mechanisms involved in the intersex condition in order to identify priority research areas. Based on current literature, the relationship between intersex and physiological parameters such as plasma sex steroids and vitellogenin is not well characterized or conclusive. Moreover, the literature is not definitive on whether field studies are distinguishing between natural intersex and intersex due to stressors. High throughput transcriptomics will improve understanding of how intersex condition manifests after exposure to aquatic pollution and it is recommended that studies consider both males with and without intersex that inhabit the same polluted site in order to differentiate pathways associated with xenobiotic responses versus molecular pathways associated with intersex. Other experimental design considerations for field studies examining intersex include data collection on life history (e.g. migratory patterns) and improved reference site characterization.