Identification of conserved genes triggering puberty in European sea bass males (Dicentrarchus labrax) by microarray expression profiling
© The Author(s). 2017
Received: 7 November 2016
Accepted: 25 May 2017
Published: 5 June 2017
The Erratum to this article has been published in BMC Genomics 2017 18:513
Spermatogenesis is a complex process characterized by the activation and/or repression of a number of genes in a spatio-temporal manner. Pubertal development in males starts with the onset of the first spermatogenesis and implies the division of primary spermatogonia and their subsequent entry into meiosis. This study is aimed at the characterization of genes involved in the onset of puberty in European sea bass, and constitutes the first transcriptomic approach focused on meiosis in this species.
European sea bass testes collected at the onset of puberty (first successful reproduction) were grouped in stage I (resting stage), and stage II (proliferative stage). Transition from stage I to stage II was marked by an increase of 11ketotestosterone (11KT), the main fish androgen, whereas the transcriptomic study resulted in 315 genes differentially expressed between the two stages. The onset of puberty induced 1) an up-regulation of genes involved in cell proliferation, cell cycle and meiosis progression, 2) changes in genes related with reproduction and growth, and 3) a down-regulation of genes included in the retinoic acid (RA) signalling pathway. The analysis of GO-terms and biological pathways showed that cell cycle, cell division, cellular metabolic processes, and reproduction were affected, consistent with the early events that occur during the onset of puberty. Furthermore, changes in the expression of three RA nuclear receptors point at the importance of the RA-signalling pathway during this period, in agreement with its role in meiosis.
The results contribute to boost our knowledge of the early molecular and endocrine events that trigger pubertal development and the onset of spermatogenesis in fish. These include an increase in 11KT plasma levels and changes in the expression of several genes involved in cell proliferation, cell cycle progression, meiosis or RA-signalling pathway. Moreover, the results can be applied to study meiosis in this economically important fish species for Mediterranean countries, and may help to develop tools for its sustainable aquaculture.
KeywordsSpermatogenesis Cell cycle Meiosis Retinoic acid Teleosts
Puberty in fish, as in other vertebrates, comprises the developmental process during which an immature individual acquires for the first time the ability to undergo sexual reproduction [1, 2]. In teleost males, puberty is tightly regulated and implies the proliferation and division of spermatogonia (mitotic phase), their subsequent entry into meiosis with the appearance of spermatocytes (meiotic phase), and the final formation of the spermatids and the haploid mature spermatozoa [2, 3]. A species-specific number of genetically determined divisions characterize the mitotic phase , whereas the meiotic phase remains under the influence of the retinoic acid (RA) signalling pathway [4, 5]. Somatic Sertoli cells are important players during spermatogenesis, exhibiting a high mitotic activity, particularly at the beginning of each seasonal cycle , and are essential for the proliferation and differentiation of germ cells . Mitosis and meiosis reveal thus as key processes for the onset of puberty in vertebrates. Meiosis is of particular importance since it implies the recombination and reduction of the genetic material, essential to ensure the correct formation of gametes, and therefore guarantees the reproduction and maintenance of the species.
Spermatogenesis is marked by the functional stimulation of the brain-pituitary-gonad (BPG) axis, responsible for its neuroendocrine control . The brain is the central organ that integrates the circuits that sense the internal and external stimuli and secretes different neuropeptides that control the production of gonadotropins from the pituitary. It is generally accepted that gonadotropins (follicle stimulating hormone; Fsh and luteinizing hormone; Lh), and androgens are the main internal stimuli for vertebrate spermatogenesis. Both gonadotropins become activated by a number of factors among which the metabolic status of the individual, in terms of body size and visceral fat content, or the photoperiod, are worth mentioning [1, 2, 8, 9]. 11ketotestosterone (11KT) is the main androgen in fish and plays an important role in the progression of spermatogenesis . It is involved in the proliferation of spermatogonia towards meiosis  and mediates the action of several factors produced by Sertoli cells like antimüllerian hormone (Amh) and insulin-like growth factors (Igfs) at the start of the cycle . In addition, Fsh has been shown to stimulate 11KT production in several fish species [13–17]. In fact, Fsh receptor is present not only in Sertoli cells but also in Leydig cells, the somatic cells with steroidogenic capability, as shown in Senegalese sole  African catfish  and zebrafish [15–17].
The European sea bass (Dicentrarchus labrax) is an important fish species for marine aquaculture that after intense research on its reproductive function has become a model for both basic and applied research. A number of studies focused on the endocrine control of reproduction shed light on the process and aided to develop protocols for its control in captivity (reviewed by ). Increased growth rates under intensive culture resulted in precocious puberty in about 20–30% of males by the end of the first year of life , something that in normal conditions occurs during the second year . As in other farmed fish, precocious maturation is one of the main drawbacks for its culture , resulting in important economic losses since by the time of marketing during the second year, precocious males exhibit a smaller size than that of normal-maturing males [20, 21]. Several studies have shown that this problem can be partially solved by a well-planned strategy of photoperiod control [18, 21–23], although clearly, more work is needed to understand why males mature precociously. Brain factors including kisspeptins [24–27], gonadotropin-releasing hormones (Gnrhs) [28, 29], leptin and their receptors  have been characterized in this species. Regarding pituitary gonadotropins, both Fsh and Lh receptors have been cloned , and different assays are currently available to measure Fsh and Lh levels in plasma and pituitary [31–33]. Moreover, recombinant European sea bass gonadotropins have been used as a biotechnological approach in gene therapy for assisted reproduction [34, 35]. In addition, secretion patterns of sex steroids , sex steroid receptors, and several steroidogenic enzymes [30, 37] during the seasonal cycle have also been reported.
Despite all the previous knowledge, little is known about the molecular machinery that triggers puberty in European sea bass males, apart from a study reporting the possible involvement of several genes in the early events preceding gonadal maturation . The recent availability of several molecular resources, including a partially annotated European sea bass genome database  gave us the opportunity to study this process using a high throughput strategy. The aim of the present work is to compare the transcriptome from European sea bass testes before and right at the start of puberty, and to identify potential genes and pathways involved in the process. This will boost our knowledge of the onset of pubertal development from a basic perspective and will help to implement tools for the improvement of a sustainable aquaculture.
Animals, rearing conditions and samplings
European sea bass hatched in April 2008 at the Ecloserie Marine de Gravelines (EMG) in the North of France and were grown there until 100 days post hatching (dph) when they were transported to our aquaria facilities at the Institute of Aquaculture Torre la Sal (IATS), a research centre belonging to the Spanish National Research Council (CSIC) in the Spanish Mediterranean coast. The facilities were approved for animal experimentation by the Ministry of Agriculture and Fisheries and by the Department of Fisheries from the Generalitat Valenciana (certificate number ES120330001055). Fish were reared in 2000 l round fiberglass tanks under natural conditions of photoperiod and temperature until the end of the experiment. In December 2008, coinciding with the first breeding season in this species, fish were subjected to abdominal massage to check for the presence of sperm. This allowed us to eliminate males exhibiting precocious puberty that could interfere with our results. These procedures were repeated every 20 days for a period of about 3 months. Prior to the start of the second breeding season (August 2009), coinciding with European sea bass normal puberty, a sampling procedure was designed to obtain testes covering the first stages of spermatogenesis. Samplings (15 fish per sampling point) were performed every 10 days starting in mid-August and finishing by the end of October. At each sampling point, fish were anesthetized with 2-phenoxyethanol (0.2 ml·l−1) and blood was taken from the caudal vein for plasma sex-steroid measurements. Fish were subsequently sacrificed by quickly severing their spinal cord and gonads were dissected for histology (the central part of the gonad) whereas the rest was kept at −80 °C for further RNA extractions needed for the different analyses including microarray hybridizations, validations, and tissue specific expression studies. Fish were treated in agreement with the Spanish regulations (Royal Decree Act 53/2013) and the European legislation (2010/63 EU) concerning the protection of animals used for experimental and other scientific purposes. All steps were taken to reduce suffering of the animals.
After dissection, the central part of the testes was immediately fixed in 4% formaldehyde: 1% glutaraldehyde in phosphate buffered saline (PBS; ). Tissues were washed in PBS and dehydrated in an increasing series of ethanol 70–96%. Samples were embedded in glycol methacrilate resin (Technovit 7100; Heraeus, Kulzer, Germany), sectioned at 3–4 μm, and stained as in . The stages of testicular development and the type and abundance of germ cells in each stage were assessed according to  and , respectively.
Steroid analysis by enzyme immune assay (EIA)
Plasma levels of 11KT were determined by enzyme immune assay (EIA) in 20 fish selected from each developmental stage, using the protocol by . Briefly, antibodies were used at a final concentration of 1:200,000 and the tracer (Cayman chemicals, MI, USA) was diluted at 1:50 Ellman Units (UE)/ml (used at 0.1042 EU/ml). The sensitivity of the assay was around 0.003 ng/ml (Bi/B0 = 90%) and half displacement (Bi/B0 = 50%) occurred around 0.03 ng/ml (slope = −1.018). The inter-assay coefficient of variation (n = 2 plates) was 1.72%.
RNA isolation and cDNA synthesis
For hybridizations and real-time validations, testes (approx. 50–100 mg) were homogenized in Trizol (Invitrogen, Carlsbad, CA) using the FastPrep® Instrument (Qbiogene, Inc., Carlsbad, CA), a tissue homogenizer with ceramic spheres as a lysing matrix. Total RNA was extracted from the lysate with the PureLink™ RNA mini Kit (Invitrogen), following the manufacturer’s instructions. Briefly, RNA was phase separated, washed, and finally eluted in DEPC water. For the tissue-distribution study, tissues including telencephalon, hypothalamus, cerebellum, spleen, gills, head kidney, kidney, liver, testis, ovary, heart and gut were homogenized in a thioglycerol-based buffer included in the Maxwell® 16 LEV simplyRNA tissue kit (Promega, Madison, WI). The homogenates were used for RNA isolation with the Maxwell® 16 instrument (Promega) following the manufacturer’s instructions that include a DNase treatment. Nevertheless, an additional test was done on the RNAs to discard any possible DNA contamination. For microarray hybridizations RNA quality was assessed with a Bioanalyzer 2100 (RNA 6000 Nano LabChip kit Agilent, Spain) and only RNAs with RIN values higher than 8.5 were used. For other downstream applications such as quantitative real time PCR (qPCR) or conventional PCR (tissue expression study), RNA quantification was done with a Nanodrop 2000c (Thermo Scientific, Wilmington, DE) and stored at −80 °C until further cDNA synthesis. Total RNA (3 μg) was reverse transcribed to cDNA with Superscript III (Invitrogen) and random hexamers following the manufaturer’s instructions. Protection of RNA from ribonucleases during cDNA synthesis was done by including 40 units of RNAse inhibitors (RNasin, Promega). The reaction was inactivated at 70 °C for 15 min.
Microarray hybridization and analysis
RNA labelling, hybridizations, and scanning were performed at the Autonomous University of Barcelona (UAB). Total RNA (100 ng) was amplified and Cy3-labeled with One-Color Microarray Gene Expression Analysis (Low Input Quick Amp Labelling kit, Agilent) along with One-Color RNA SpikeIn Kit (Agilent) following the manufacturer’s instructions. The resulting cRNA was purified (RNeasy mini spin columns; Qiagen), quantified with a Nanodrop ND-1000 and checked with a Bioanalyzer 2100 as previously described. Amplified samples (1.65 μg per sample) were hybridized to a custom oligonucleotide high-density European sea bass microarray (Agilent 4 × 44 K design format; http://www.agilent.com/) containing 60-mer oligonucleotides with a linker directly spotted on glass slides using the Agilent’s SurePrint Tecnology. Three samples from each testicular stage, selected after histological examination and 11KT plasma levels, were used for microarray hybridizations, each of them consisting of a pool RNAs from six males. The pools were used as biological replicates and thus independent samples for microarray hybridizations. In addition, and since each microarray plate can hold up to four samples (4 × 44 design), one sample from each stage was randomly chosen and hybridized in both plates as a quality control to check for possible inter-plate hybridization differences. The probes contained in the microarray (GEO accession number GPL13443) cover 13,199 unique sequences of Dicentrarchus labrax that include 6275 annotated transcripts, each with 3 specific probes, and 6924 ESTs with 1 probe/target sequence. Assuming that a typical diploid teleost genome is expected to have 26–28 thousand protein coding loci, the microarray used for the study should cover about half of the genes of the species. Hybridizations were done at 65 °C for 17 h (GE Hybridization Kit; Agilent). Washes were conducted as recommended by the manufacturer using Agilent’s Gene Expression Wash Pack with stabilization and drying solution and arrays were scanned with a G2505B (Agilent). Several quality control features and spot intensities were extracted with Agilent’s Feature Extraction software v10.4. Finally, data were analyzed with GeneSpring software v10.1. Percentile shift normalization was used to adjust all spot intensities in the array (percentile target = 75). Principal Component Analysis (PCA) was used as a quality control on samples and allowed to decrease the number of false positives before the statistical analysis. Normalized data were filtered by comparison of the standard deviation expression among groups (filter by expression). Statistical analyses were performed on filtered data using a t-test. Significant differences in the transcriptomic profile between early stages of spermatogenesis (data filtered at a fold change (FC) expression of 2) were set at p < 0.01. The corresponding study was deposited at the Gene Expression Omnibus (GEO-NCBI) database under the accession number GSE47400.
Gene annotation and enrichment analysis
The web-based tools Genecards (http://www.genecards.org), Uniprot (http://www.uniprot.org) and AmiGO 2 (amigo.http://amigo.geneontology.org) were used to assign gene names, synonyms and functions to the differentially expressed genes (DEGs) found after microarray hybridizations. The annotation of the sequences was manually curated, improving the accuracy of the information obtained from the microarray used for this study. A further improvement was added implementing the Blast2Go software  that enriched the number of GO-term annotations. A list containing all genes included in our custom-made microarray was used as a reference set to evaluate the enrichment in GO-terms in the subset of DEGs. The resulting data were analysed with Fisher’s exact test with multiple testing correction of the false discovery rate. In addition, annotated DEGs were ascribed to functional biological pathways using the Kyoto Encyclopaedia of Genes and Genomes (http://www.genome.jp/kegg) and the possible altered metabolic pathways were assessed.
Array validation by quantitative real-time PCR (qPCR)
EST sequences of the DEGs were used as a query in Blast searches against the European sea bass genome and GeneBank databases in order to position the selected DEGs in their corresponding genes. Primers for the amplification of the DEGs were designed in areas covering intron-exon boundaries to check for genomic contamination using Primer 3 (http://primer3.ut.ee). Primers (Additional file 1) were checked by conventional PCR and the amplified fragments sequenced to verify their identities. qPCR analyses were performed with an iCycler iQ™ (BioRad Labs., Inc.) using SYBR® Green (PCR Master Mix; Applied Biosystems). PCR reactions were run in triplicate in optically clear 96-well plates in a final 20 μl volume containing 10 μl of 2× Sybr Green Master mix, 10 pmol of each primer and 5 μl of diluted cDNA (1:50 for the target genes or 1:500 for the reference gene). Cycling parameters included an initial denaturation at 95 °C for 3 min, followed by 40 cycles at 95 °C for 15 s and annealing-extension at 60–72 °C for 1 min ending with an extension at 72 °C for 1 min. A final temperature dissociation step was done to ensure the presence of just one product. qPCR data were collected with iCycler™ iQ optical system software (v. 3.0, BioRad). The cycle threshold (Ct) was calculated as the average of three replicates per sample. Gene expression analyses were conducted using the Q-Gene core module . Briefly, for each gene the amplification efficiency (E) was calculated from the slope of the linear correlation between Cts and the logarithm of the amount of serially diluted RNA, used as a standard, following the eq. E = 10(−1/slope). E values for the different genes were within the range of 93.5–101.8%. Values were normalized (normalized expression; NE) to the constitutively expressed reference gene 18S rRNA in each sample (n = 8 individual fish per stage and gene) according to the eq. NE = (Eref)Ctref/(Etarget)Cttarget. 18S rRNA was considered a good reference gene since it exhibited the best bestkeeper index when comparing different developmental stages . In addition, the expression of this gene remains constant in many physiological conditions such as differentiation and proliferation  making it a suitable reference gene for this study.
Tissue specific expression
The expression of the selected DEGs was assessed in different tissues including telencephalon, hypothalamus, cerebellum, spleen, gills, head kidney, kidney, liver, testis, ovary, heart and gut. PCR reactions were performed with an initial denaturation of 5 min at 94 °C, and then 34 cycles with the following characteristics: denaturation at 94 °C for 30 s, annealing at 60 °C for 30 s, and extension at 72 °C for 30 s. A final extension of 2 min at 72 °C was added at the end of the 34 cycles.
Sequencing, cloning, and phylogenetic studies of European sea bass cyp26a1
Based on the ESTs of the microarray, the full sequence of European sea bass cyp26a1 was localized along the genome. Specific primers were designed in 3′- and 5′-UTR flanking regions to amplify its full-coding sequence. The fragment was cloned into a bacterial vector using the pGEM T-easy cloning kit (Promega Corp., Madison, WI), and amplified in E. coli competent cells following the manufacturer’s instructions. Several colonies were selected, grown in liquid LB and finally sequenced with an automatic ABI 3100 Genetic Analyser (Applied Biosystems, Foster City, CA), using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems). The identity of the clones was confirmed after sequencing, multiple alignment comparisons, and phylogenetic analysis. An in silico study of the 5′ upstream 1500 bp of the flanking promoter sequence using MatInspector and (Promo) Transfac v.8.3 was used to identify the presence of putative binding sites for transcription factors that could be involved in the activation or repression of cyp26a1 transcription. An alignment of known Cyp26 proteins from vertebrates, either compiled from GenBank/EMBL or predicted in ENSEMBL, was made with clustalW2. For the phylogenetic tree, the distances were computed with the Poisson correction method  and the evolutionary history was inferred using the Neighbor-Joining method  after a bootstrap test using 1000 replicates. The phylogenetic analysis was carried out in MEGA v.4 .
Student’s t-test for hormonal analysis, microarray hybridizations and gene expression levels was used to reveal significant differences between stage I and stage II. In all cases, significant differences were accepted at p < 0.05 except for microarray hybridizations for which differences were accepted at p < 0.01.
Sample selection: histological and hormonal classification
Gene ontology analysis of annotated transcripts significantly affected during the onset of puberty in European sea bass testis (Fisher’s exact test with multiple corrections for FDR)
Sample frequency (N = 141)
Array frequency (N = 7681)
GO:0007049. Cell cycle
uhrf2, dmc1, nsl1, ccnd2, kif2c, cenpi, thbs1 , jmy , chaf1b, psmd3, ncapg, rbbp4, syce1, ttk, sycp2, pold3, camk2d , ndc80, aspm, rad9b, aurkb, psmb7, bub3, dtymk, ppef1, cenph, anln, slbp, spc25, nup37, mad2l1bp, ccne2, trip13, cdc28, nsmce2, cenpf
GO:0051301. Cell division
nup37, nsmce2, anl, aurkb, cenph, aspm, ndc80, mad2l1bp, bub3, ppef1, cdc28
GO:0007059. Chromosome segregation
nup37, nsmce2, aurkb, cenph, ndc80, kif2c, madl1bp, cenpf, ttk, bub3, syce1, ncapg
GO:0006259. DNA metabolic process
uhrf2, dmc1, cacybp, jmy , chaf1b, pcna, ncapg, rbbp4, pold3, rfc3, rad9b, cry2 , smarcc2, asf1b, fen1, trip13, ruvbl2, nt5e , mcm3, nsmce2
GO:0010467. Gene expression
uhrf2, thbs1 , eed, asf1b
tgfbr1 , adamts1 , cenpi, bub3, syce1, ihh , sycp2, dmc1, hsf2bp, amh , trip13, sdf1a, ttk, fosla
GO:0044237. Cellular metabolic process
uhrf2, dmc1, rab6a , hck, cacybp, thbs1 , jmy , chaf1b, pcna, psmd3, ncapg, cry2 , gpd1, cad, rbbp4, nme1, ttk, pold3, camk2d , rfc3, lox , srsf7, rad9b, eed, ifitm2 , fen1, smarcc2, asf1b, nme3, aurkb, shmt1, cpsf3, atp6v0e1, dtymk, nsmce2, ppef1, u2af35, slbp, trip13, lnx1 , fkbp8 , nnt, ruvbl2, tgfbr1 , nt5e , ubr7, mcm3
GO:0050794. Regulation of cellular process
bnc1 , ihh , agrp2 , rab6a , thbs1 , c-fosla , psmd3, q2laq1 , amh , ect2, cad, trip13, camk2d , ndc80, rad9b, stat3 , cry2 , sdf1a , ppef1, depdc1b, crabp1, isg20l2, lbr, ap2s1, cyp26a1 , fkbp8 , tfrc, adamts1 , tgfbr1
GO:0043169. Cation binding
ihh , mgp , thbs1 , cad, tppp3 , ppef1, slc25a25 , anxa2
GO:0043167. Ion binding
ihh , mgp , thbs1 , cad, tppp3 , ppef1, slc25a25 , anxa2
GO:0046872. Metal ion binding
ihh , mgp , thbs1 , cad, tppp3 , ppef1, slc25a25 , anxa2
bnc1 , apeh, uhrf2, dmc1, ihh , ccnd2, cacybp, cenpi, jmy , c-fosla , chaf1b, pcna, trhb , sh3bgrl3 , rbbp4, syce1, sycp2, pold3, camk2d , rfc3, ndc80, aspm, lox , srsf7, rad9b, eed, stat3 , cry2 , smarcc2, asf1b, smox , aurkb, shmt1, cpsf3, lbr, psmb7, b9d2, bub3, gfi1 , fen1, ass1, ppef1, u2af35, cenph, chrac1, slbp, nup37, cenpf, mad2l1bp, ccne2, trip13, isg20l2, ruvbl2, mcm3, cdc28
uhrf2, dmc1, nsl1, kif2c, cenpi, pcna, ncapg, rbbp4, syce1, pold3, rfc3, ndc80, rad9b, eed, asf1b, aurkb, bub3, u2af35, spc25, nup37, mad2l1bp, ruvbl2, mcm3, nsmce2, cenpf
GO:0000228. Nuclear chromosome
dmc1, pcna, rbbp4, syce1, pold3, ndc80, rad9b, eed, aurkb, ruvbl2, mcm3
bnc1 , ihh , thbs1 , dsg2 , cad, camk2d , slco2a1 , ppef1, rims1, ap2s1, tfrc, tgfbr1 , nt5e , synpo
Affected KEGG pathways at the onset of European sea bass puberty
- Nucleotide metabolism
nme3, pold3, cad, tyms, nt5e , dtymk
- Retinol metabolism
- Lipid metabolism
pnpla2 , gpd1, agrp2 , ptges
- Amino acid metabolism
odc1, pycr2, smox , ass1, cad, shmt1, pah
- Cell cycle
bub3, mcm3, pcna, ttk, ccnd2, mad2l1bp, cdc28, aurkb, ndc80, spc25
mcm3, cdc28, dmc1, mad2l2, sycp2, syce1, cenph, cenpi, cenpf, ndc80, spc25
- DNA replication and repair
mcm3, pold3, fen1, pcna, rfc3
- Focal adhesion
actb , ccnd2, thbs1
Genetic information processing
- FoxO signalling pathway
tgfbr1 , stat3 , plk4, ccnd2
- Hippo signalling pathway
amh , tgfbr1 , ccnd2, actb
- Jak-STAT signalling pathway
stat3 , socs3 , ccnd2
- TGF beta signalling pathway
amh , tgfbr1 , thbs1
- MAPK signalling pathway
gpd1, cdc28, tgfbr1
- TNF signalling pathway
- Toll-like receptor signalling pathway
- Wnt signalling pathway
- PI3K-Akt signalling pathway
- cAMP signalling pathway
- Hedgehog signalling pathway
- Rap1 signalling pathway
- Prolactin signalling pathway
stat3 , socs3 , ccnd2
- Renin secretion
- Thyroid hormone signalling pathway
Molecular cloning of European sea bass cyp26a1 and phylogenetic analysis
Tissue specific expression
The present study contributed to identify differences in gene expression during the early stages of pubertal development in European sea bass males using a custom-made microarray. However, and despite the fact that the differences observed at the transcriptome and the steroid level are related to puberty (first successful reproduction), it is possible that similar changes could be found in successive reproductive seasons since they mark the transition between quiescence and the start of gametogenesis. The use of long oligo-based microarrays has been shown to have a higher sensitivity for detection but a lower specificity  and therefore could have a reduced ability to discriminate between similar transcripts produced by the same locus, paralogs or similar members of large gene families. To circumvent this problem, we cloned the full-length sequence of several selected DEGs and then validated the results with qPCR obtaining a good correlation between both methods. In addition, the study revealed the complete sequences of some transcripts for the first time in the European sea bass, adding contrasted information to the microarray that was based on EST sequences. However, for other DEGs that were not annotated in the microarray, we could not find any match to reveal their identities in any of the databases searched including Genebank, Ensembl, and Uniprot. The different stages of spermatogenesis (stage I and stage II) used for the study were classified by histology and their corresponding 11KT levels were further confirmed by EIA, demonstrating that the increase of circulating 11KT marks the initiation of pubertal development. A similar result was found in other teleosts including eel , goldfish , zebrafish  or trout  and 11KT measured from the mucus of carps was suitable to differentiate males from females [59, 60]. However, no correlation between 11KT and gonad developmental stage in either male or female carp could be found . Our results open the possibility to explore the use of 11KT as a non-lethal marker for the onset of puberty in this species helping to manage the fish farms stocks to separate precocious from non-precocious European sea bass males.
The transcriptome response revealed that cell proliferation, cell cycle and meiosis progression were pathways preferentially affected during the onset of male puberty. The spindle assembly checkpoint (SAC) is a control mechanism of dividing cells that ensures the correct segregation of chromosomes by blocking cell cycle progression until kinetochores are properly connected to the spindle . In our study, several genes coding for SAC proteins (bub3) and SAC protein regulators (mad2l1bp and ttk), were upregulated in stage II testes as well as other important kinetochore-associated transcripts such as ndc80, spc25, aurkb and cdc28 (Additional file 2). ndc80 and spc25 code for essential proteins of the Ndc80 complex, needed for SAC activity  while aurkb (aurora kinase b) controls kinetochore orientation during meiosis . Defects in cdc28 function result in delays in the exit from mitosis and in meiosis impairment among others . Transcripts coding for centromere proteins like cenph, cenpi or cenpf, were also upregulated. In the case of cenpi, in addition to its role in centromere formation, it is involved in the response of gonadal tissues to Fsh . The upregulation of cenpi in stage II coincides with the initiation of the gradual increase in Fsh plasma levels in European sea bass during early spermatogenesis [1, 33]. This is in line with the role of Fsh inducing germ cell proliferation and marks the onset of spermatogenesis through the activation of spermatogenesis-related genes . Several transcripts like sycp2 and syce1, coding for proteins of the synaptonemal complex and trip13, required for the completion of meiosis , were upregulated in stage II testes. Moreover, pcna, essential for DNA replication and a molecular marker of dividing cells , also increased during the onset of spermatogenesis in agreement with the active mitosis of spermatogonia typical of this period . The role of pcna in the proliferation of germ cells has been described in several teleosts and is currently used in a number of fish species as a marker of spermatogenesis progression . All these results are supported by the analysis of affected biological pathways (Table 2) that showed higher expression of genes involved in cellular processes, particularly those involved in cell cycle, meiosis and DNA replication and repair. In addition, several signalling pathways involved in testicular development such as Wnt, MAPK, hedgehog and TGF beta signalling pathways  were altered during early puberty in European sea bass. Altogether, the upregulation of the above mentioned genes is indicative of an active period of mitosis, reflect the need for a tight control of the correct division of the cells, and constitutes an indicator for the progression of meiosis typical of this stage.
A second group of DEGs includes those implicated in reproduction and growth. Amh is involved in gonadal development and steroidogenesis in vertebrates and induces the regression of müllerian ducts in mammals during male embryogenesis . Although fish do not have müllerian ducts, amh homologues have been identified in several teleosts , suggesting evolutionary conserved functions for this gene. The role of Amh as a meiosis inhibiting factor was first shown in eels  and recently in zebrafish , induced by the increase of circulating 11KT that blocked amh expression facilitating spermatogenesis completion. In teleosts, amh has a key role in early testicular maturation with highest levels in pre-spermatogenic testis and lowest during spawning . In European sea bass, administration of recombinant Fsh induced spermatogonial proliferation and differentiation into spermatocytes, due to the increase of 11KT levels and the concomitant suppression of amh expression . Moreover, amh mRNA and protein expression was detected in Sertoli cells of prepubertal European sea bass, and the signal decreased during spermatogenesis . Our results showed a decrease of amh levels during early spermatogenesis in agreement with its role as an inhibitor of spermatogenesis progression. This is supported by the increase of pcna, cenpi and 11KT levels in the same testicular stage. Likewise, a downregulation of amh during the reproductive cycle has been found in the testicular transcriptome of rainbow trout  and in precocious Atlantic salmon , further demonstrating that the inhibitory effect of Amh of the onset of puberty can be extended to all fish species so far studied. In addition, the study of the biological pathways affected during the onset of puberty, also show the importance of amh in several signalling pathways including those of TGF beta, hippo and cAMP . Few studies are available for sgII in fish, apart from those in goldfish [75–77]. SgII is widely distributed in secretory granules of neurons and endocrine cells  and is the precursor of secretoneurin, a bioactive neuropeptide capable to induce Lh secretion [75, 77]. Our results show low sgII expression by the onset of spermatogenesis, in agreement with the low Lh levels found at this stage in the European sea bass [1, 30]. It would be very interesting to determine whether sgII, and therefore secretoneurin, also increase during later stages of spermatogenesis, coinciding with the surge of Lh, to experimentally test this hypothesis. Pioneering studies in sea bream point at the relevance of aquaporin 1 (aqp1) in fish reproduction due to its role in water intake during oocyte hydration prior to spawning  and in the activation of sperm motility during the last stages of spermatogenesis [80, 81]. Our results show low aqp1 levels during early spermatogenesis, in agreement with its prominent role in sperm maturation during the last stages of spermatogenesis. Agrp2 (agouti-related protein 2) is an orexigenic peptide with a key role in the regulation of energy balance in mammals  and fish . In this regard, our study shows that the lipid metabolism pathway where agrp2 was included was affected during the early stages of pubertal development further supporting its role in energy balance. A direct link between leptin, the most powerful orexigenic neuropeptide in fish , and the AGRP system has been suggested in European sea bass males . Moreover, abundant Agrp expression was found in mouse pachytene-spermatocytes and immunohistochemistry revealed that Agrp co-localized with Scp3, a meiotic-specific protein of the synaptonemal complex . Although agrp2 has been characterized in European sea bass testis , this is the first time its involvement in spermatogenesis is suggested, possibly due to the specific energy requirements during spermatogenesis and the decrease in food intake. This is in agreement with the downregulation of agrp2 in European sea bass brain after long-term fasting  and its decrease in testis during early spermatogenesis (present study), and link the appetite and growth system with reproduction . Moreover, a transcriptomic study of trout testis revealed that Fsh administration induced the increase of igfbp6 . This strong Fsh-induced upregulation was present during early spermatogenesis, including germ cell proliferation and meiosis, and was associated to the effect of the Igf-signalling pathway on spermatogenesis progression . Our results also show a clear upregulation of igfbp6, coinciding with the first stages of pubertal development and the increase of 11KT plasma levels.
The last group of DEGs is associated with the RA-signalling pathway. RA has been proposed as a meiosis inducing factor in tetrapods including mammals , birds , amphibians  and fish [4, 5, 54, 90]. Two transcripts coding for binding proteins, one in charge of retinol transport (rbp4) through the blood stream and another one (crabp1) in charge of the translocation of RA to the nucleus of the target cells  were differentially expressed in European sea bass transcriptome. In addition, cyp26a1, responsible for the degradation of intracellular RA and essential for the maintenance of RA homeostasis , was affected. The decrease of cyp26a1 in stage II is associated with a decrease in the translocation and transport of RA brought about by the downregulation of rbp4 and crabp1, in order to maintain the homeostasis of RA that otherwise, and at high levels can be toxic for the cell . In addition, retinol metabolism was one of the affected metabolic pathways found in the present study. It seems thus plausible that in European sea bass, the suppression of RA degradation and the concomitant increase in the availability of RA could be partially responsible for triggering the onset of meiosis. In zebrafish testes, cyp26a1 was expressed in germ cells entering meiosis, while in females, a downregulation was found in oocytes during meiosis resumption . Likewise, in medaka, RA was found to act directly on Sertoli cells, Leydig cells, and pre-meiotic germ cells with a decrease of cyp26a1 expression by the time of meiosis resumption, whereas in ovaries, RA-transcriptional activity is highest in meiotic oocytes . In addition, in vivo Fsh-injection to pre-spermatogenic zebrafish males induced the onset of spermatogenesis and resulted in changes of several enzymes involved in the RA-signalling pathway, including a decrease in cyp26a1 expression, although no effect was found after ex vivo culture of pre-spermatogenic testes with Fsh . Moreover, the administration of an inhibitor of RA synthesis in combination with a deficient diet of vitamin A (a precursor of RA) to adult zebrafish also induced a downregulation of cyp26a1, most likely to increase intracellular RA levels, although spermatogenesis was still disrupted, and fertility compromised . The above mentioned studies suggest that a decrease in cyp26 expression is associated with the onset of spermatogenesis and the initiation of meiosis and prompted us to clone and obtain the full length of its cDNA. The alignment of the deduced protein sequence with other Cyp26 proteins available from other vertebrates confirmed its identity as Cyp26a1 revealing slight differences in homology among teleosts, further supporting their high conservation due to its pivotal role controlling RA levels. The tree shows the presence of a common ancestor cyp26 protein, strengthening the hypothesis of an independent functionalization of its coding gene prior to the two rounds of genome duplication in vertebrates . Moreover, the study of the promoter showed the presence of binding sites for cell cycle regulators and for RA nuclear receptors including Ppar, Rxr and Rar, and also DR1 and DR5 sites (RA-responsive elements) indicating the role of RA, via the interaction with its nuclear receptors, in the regulation of cyp26a1 transcription in the European sea bass. A similar result has been shown for the cyp26a1 promoter in zebrafish [94, 95] and medaka , although the development of functional studies is clearly needed to confirm the capability of RA to induce the regulation of cyp26a1 in the European sea bass. To gain more insight on the importance of the RA-signalling pathway in meiosis we studied the expression of several nuclear receptors involved in RA binding including rarα, raxrα, and pparγ. The receptors appeared ubiquitously expressed and at similar levels in all the tissues studied, reflecting the general actions and the importance of RA in numerous biological processes throughout evolution  and its involvement in the proliferation and differentiation of many cell types .
To the best of our knowledge, this is the first transcriptomic study focussed on the early stages of puberty, and aimed at the identification of molecular and endocrine signals triggering the start of the initial spermatogenic wave in European sea bass. Increases in androgen plasma levels, particularly 11KT, mark the transition between testicular stage I and stage II. This opens the possibility to explore the use of 11KT in the management of European sea bass stocks in aquaculture farms to separate precocious from non-precocious males. The study improved the annotation of different genes of the microarray and helped to increase the knowledge of several mechanisms and biological pathways involved in early stages of puberty. Altogether, the study shows that the onset of spermatogenesis is characterized by the activation of genes involved in cell cycle progression and division including mitosis and meiosis. The differential expression of several components of the RA-signalling pathway suggests their important role in the onset of meiosis. This work lays the foundation for an in-depth study of the RA-signalling signalling pathway and its role in the onset of meiosis in fish. A future increase in the sequencing of the European sea bass gonad transcriptome and the use of RNA-seq technologies will help to shed light on the molecular pathways involved in relevant aspects of the reproductive process of this economically important species and will aid to develop comparative studies on gonadal differentiation and maturation in teleosts.
We thank Sebastián Boltaña, Nerea Roher and Benjamí Piña for technical support on microarray processing data.
This work was supported by Aquagenomics project (CSD2007–0002) from the Spanish Ministry of Science and Innovation to SZ and a MICINN project (AGL2011–28890) to AG. PM was supported by a grant from the University of Antofagasta MECE2 (ANT0806).
Availability of data and materials
The complete design was deposited at the Gene Expression Omnibus NCBI database GSE47400 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=exgnmkgcfbwnhkp&acc=GSE47400 and the platform that validates the microarray can be found at GPL13443.
MB, AG and SZ designed the study; PM, BC and MB conducted the experiments, performed the molecular biology studies and the statistical and microarray analyses. MB drafted the initial manuscript; MB, PM, BC, AG and SZ wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing financial and non-financial interests.
Consent for publication
Ethics approval and consent to participate
Fish were treated in agreement with the Spanish regulations (Royal Decree Act 53/2013) and the European legislation (2010/63 EU) concerning the protection of animals used for experimental and other scientific purposes. The experiments were performed following the approval from the institutional regulations on animal experimental handling review board from the IATS-CSIC and from the CSIC ethics committee (project number CSD2007–002). All steps were taken to reduce suffering of the animals.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Carrillo M, Espigares F, Felip A, Escobar S, Molés G, Rodríguez R, et al. Updating control of puberty in male European sea bass: a holistic approach. Gen Comp Endocrinol. 2015;221:42–53.Google Scholar
- Schulz RW, de França LR, Lareyre J-J, LeGac F, Chiarini-Garcia H, Nobrega RH, et al. Spermatogenesis in fish. Gen Comp Endocrinol. 2010;165(3):390–411.Google Scholar
- Nóbrega R, Batlouni S, França L. An overview of functional and stereological evaluation of spermatogenesis and germ cell transplantation in fish. Fish Physiol Biochem. 2009;35(1):197–206.View ArticlePubMedGoogle Scholar
- Feng R, Fang L, Cheng Y, He X, Jiang W, Dong R, et al. Retinoic acid homeostasis through aldh1a2 and cyp26a1 mediates meiotic entry in Nile tilapia (Oreochromis niloticus). Sci Rep. 2015;5:10131.View ArticlePubMedPubMed CentralGoogle Scholar
- Adolfi MC, Herpin A, Regensburger M, Sacquegno J, Waxman JS, Schartl M. Retinoic acid and meiosis induction in adult versus embryonic gonads of medaka. Sci rep. 2016;6:34281.View ArticlePubMedPubMed CentralGoogle Scholar
- Griswold MD. The central role of Sertoli cells in spermatogenesis. Semin Cell dev Biol. 1998;9:411–6.View ArticlePubMedGoogle Scholar
- Zohar Y, Muñoz-Cueto JA, Elizur A, Kah O. Neuroendocrinology of reproduction in teleost fish. Gen Comp Endocrinol. 2010;165(3):438–55.View ArticlePubMedGoogle Scholar
- Escobar S, Rocha A, Felip A, Carrillo M, Zanuy S, Kah O, et al. Leptin receptor gene in the European sea bass (Dicentrarchus labrax): cloning, phylogeny, tissue distribution and neuroanatomical organization. Gen Comp Endocrinol. 2016;229:100–11.Google Scholar
- Taranger GL, Carrillo M, Schulz RW, Fontaine P, Zanuy S, Felip A, et al. Control of puberty in farmed fish. Gen Comp Endocrinol. 2010;165(3):483–515.Google Scholar
- Schulz RW, Miura T. Spermatogenesis and its endocrine regulation. Fish Physiol Biochem. 2002;26(1):43–56.View ArticleGoogle Scholar
- Miura T, Miura CI. Molecular control mechanisms of fish spermatogenesis. Fish Physiol Biochem. 2003;28(1–4):181–6.View ArticleGoogle Scholar
- Rolland AD, Lareyre J-J, Goupil A-S, Montfort J, Ricordel M-J, Esquerré D, et al. Expression profiling of rainbow trout testis development identifies evolutionary conserved genes involved in spermatogenesis. BMC Genomics. 2009;10(1):546.Google Scholar
- Chauvigné F, Verdura S, Mazón MJ, Duncan N, Zanuy S, Gómez A, et al. Follicle-stimulating hormone and luteinizing hormone mediate the androgenic pathway in Leydig cells of an evolutionary advanced teleost. Biol Reprod. 2012;87(2):1–11.Google Scholar
- García-López A, Bogerd J, Granneman JC, van Dijk W, Trant JM, Taranger GL, et al. Leydig cells express follicle-stimulating hormone receptors in African catfish. Endocrinol. 2009;150(1):357–65.Google Scholar
- García-López A, de Jonge H, Nóbrega RH, de Waal PP, van Dijk W, Hemrika W, et al. Studies in zebrafish reveal unusual cellular expression patterns of gonadotropin receptor messenger ribonucleic acids in the testis and unexpected functional differentiation of the gonadotropins. Endocrinol. 2010;151(5):2349–60.Google Scholar
- Skaar K, Nobrega R, Magaraki A, Olsen L, Schulz R, Male R. Proteolytically activated, recombinant anti-müllerian hormone inhibits androgen secretion, proliferation, and differentiation of spermatogonia in adult zebrafish testis organ cultures. Endocrinol. 2011;152(9):3527–40.View ArticleGoogle Scholar
- Crespo D, Assis LH, Furmanek T, Bogerd J, Schulz RW. Expression profiling identifies Sertoli and Leydig cell genes as Fsh targets in adult zebrafish testis. Mol Cell Endocrinol. 2016;437:237–51.View ArticlePubMedGoogle Scholar
- Begtashi I, Rodríguez L, Moles G, Zanuy S, Carrillo M. Long-term exposure to continuous light inhibits precocity in juvenile male European sea bass (Dicentrarchus labrax, L.). I. Morphological aspects. Aquaculture. 2004;241(1):539–59.View ArticleGoogle Scholar
- Zanuy S, Carrillo M, Felip A, Rodríguez L, Blázquez M, Ramos J, et al. Genetic, hormonal and environmental approaches for the control of reproduction in the European sea bass (Dicentrarchus labrax L.). Aquaculture. 2001;202(3):187–203.Google Scholar
- Felip A, Zanuy S, Carrillo M. Comparative analysis of growth performance and sperm motility between precocious and non-precocious males in the European sea bass (Dicentrarchus labrax, L.). Aquaculture. 2006;256(1):570–8.View ArticleGoogle Scholar
- Felip A, Zanuy S, Muriach B, Cerdá-Reverter JM, Carrillo M. Reduction of sexual maturation in male Dicentrarchus labrax by continuous light both before and during gametogenesis. Aquaculture. 2008;275:347–55.View ArticleGoogle Scholar
- Rodríguez L, Begtashi I, Zanuy S, Carrillo M. Long-term exposure to continuous light inhibits precocity in European male sea bass (Dicentrarchus labrax, L.): hormonal aspects. Gen Comp Endocrinol. 2005;140(2):116–25.View ArticlePubMedGoogle Scholar
- Rodríguez L, Carrillo M, Sorbera LA, Zohar Y, Zanuy S. Effects of photoperiod on pituitary levels of three forms of GnRH and reproductive hormones in the male European sea bass (Dicentrarchus labrax, L.) during testicular differentiation and first testicular recrudescence. Gen Comp Endocrinol. 2004;136(1):37–48.View ArticlePubMedGoogle Scholar
- Alvarado M, Carrillo M, Felip A. Expression of kisspeptins and their receptors, gnrh-1/gnrhr-II-1a and gonadotropin genes in the brain of adult male and female European sea bass during different gonadal stages. Gen Comp Endocrinol. 2013;187:104–16.View ArticlePubMedGoogle Scholar
- Escobar S, Felip A, Gueguen M-M, Zanuy S, Carrillo M, Kah O, et al. Expression of kisspeptins in the brain and pituitary of the European sea bass (Dicentrarchus labrax). J Comp Neurol. 2013;521(4):933–48.Google Scholar
- Escobar S, Servili A, Espigares F, Gueguen M-M, Brocal I, Felip A, et al. Expression of kisspeptins and kiss receptors suggests a large range of functions for kisspeptin systems in the brain of the European sea bass. PLoS One. 2013;8(7):e70177.Google Scholar
- Felip A, Zanuy S, Pineda R, Pinilla L, Carrillo M, Tena-Sempere M, et al. Evidence for two distinct KiSS genes in non-placental vertebrates that encode kisspeptins with different gonadotropin-releasing activities in fish and mammals. Mol Cell Endocrinol. 2009;312(1):61–71.Google Scholar
- González-Martínez D, Zmora N, Mañanos E, Saligaut D, Zanuy S, Zohar Y, et al. Immunohistochemical localization of three different prepro-GnRHs in the brain and pituitary of the European sea bass (Dicentrarchus labrax) using antibodies to the corresponding GnRH-associated peptides. J Comp Neurol. 2002;446(2):95–113.Google Scholar
- González-Martínez D, Madigou T, Mañanos E, Cerdá-Reverter JM, Zanuy S, Kah O, et al. Cloning and expression of gonadotropin-releasing hormone receptor in the brain and pituitary of the European sea bass: an in situ hybridization study. Biol Reprod. 2004;70(5):1380–91.Google Scholar
- Rocha A, Zanuy S, Carrillo M, Gómez A. Seasonal changes in gonadal expression of gonadotropin receptors, steroidogenic acute regulatory protein and steroidogenic enzymes in the European sea bass. Gen Comp Endocrinol. 2009;162(3):265–75.View ArticlePubMedGoogle Scholar
- Mateos J, Mañanós EL, Swanson P, Carrillo M, Zanuy S. Purification of luteinizing hormone (LH) in the sea bass (Dicentrarchus labrax) and development of a specific immunoassay. Ciencias Marinas. 2006;32(2):271–83.Google Scholar
- Molés G, Gómez A, Carrillo M, Rocha A, Mylonas CC, Zanuy S. Determination of Fsh quantity and bioactivity during sex differentiation and oogenesis in European sea bass. Biol Reprod. 2011;85(4):848–57.View ArticlePubMedGoogle Scholar
- Molés G, Gómez A, Carrillo M, Zanuy S. Development of a homologous enzyme-linked immunosorbent assay for European sea bass FSH. Reproductive cycle plasma levels in both sexes and in yearling precocious and non-precocious males. Gen Comp Endocrinol. 2012;176(1):70–8.View ArticlePubMedGoogle Scholar
- Mazón MJ, Gómez A, Yilmaz O, Carrillo M, Zanuy S. Administration of follicle-stimulating hormone in vivo triggers testicular recrudescence of juvenile European sea bass (Dicentrarchus labrax). Biol Reprod. 2014;90(1):6. 1-10View ArticlePubMedGoogle Scholar
- Mazón MJ, Molés G, Rocha A, Crespo B, Lan-Chow-Wing O, Espigares F, et al. Gonadotropins in European sea bass: endocrine roles and biotechnological applications. Gen Comp Endocrinol. 2015;221:31–41.Google Scholar
- Prat F, Zanuy S, Carrillo M, De Mones A, Fostier A. Seasonal changes in plasma levels of gonadal steroids of sea bass, Dicentrarchus labrax L. Gen Comp Endocrinol. 1990;78(3):361–73.View ArticlePubMedGoogle Scholar
- Blázquez M, González A, Papadaki M, Mylonas C, Piferrer F. Sex-related changes in estrogen receptors and aromatase gene expression and enzymatic activity during early development and sex differentiation in the European sea bass (Dicentrarchus labrax). Gen Comp Endocrinol. 2008;158(1):95–101.View ArticlePubMedGoogle Scholar
- Crespo B, Gómez A, Mazón MJ, Carrillo M, Zanuy S. Isolation and characterization of Ff1 and Gsdf family genes in European sea bass and identification of early gonadal markers of precocious puberty in males. Gen Comp Endocrinol. 2013;191:155–67.View ArticlePubMedGoogle Scholar
- Tine M, Kuhl H, Gagnaire P-A, Louro B, Desmarais E, Martins RS, et al. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat Commun. 2014;5:5770.Google Scholar
- McDowell E, Trump B. Histologic fixatives suitable for diagnostic light and electron microscopy. Arch Pathol Lab Med. 1976;100(8):405–14.PubMedGoogle Scholar
- Bennett HS, Wyrick AD, Lee SW, McNeil JH. Science and art in preparing tissues embedded in plastic for light microscopy, with special reference to glycol methacrylate, glass knives and simple stains. Stain Tech. 1976;51(2):71–97.View ArticleGoogle Scholar
- Espigares F, Rocha A, Molés G, Gómez A, Carrillo M, Zanuy S. New insights into the factors mediating the onset of puberty in sea bass. Gen Comp Endocrinol. 2015;224:176–85.View ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.View ArticlePubMedGoogle Scholar
- Müller PY, Janovjak H, Miserez AR, Dobbie Z. Processing of gene expression data generated by quantitative real-time RT-PCR. BioTechniques. 2002;32:1372–9.PubMedGoogle Scholar
- Mitter K, Kotoulas G, Magoulas A, Mulero V, Sepulcre P, Figueras A, et al. Evaluation of candidate reference genes for QPCR during ontogenesis and of immune-relevant tissues of European seabass (Dicentrarchus labrax). Comp Biochem Physiol B: Biochem Mol Biol. 2009;153(4):340–7.Google Scholar
- Fuster G, Vicente R, Coma M, Grande M, Felipe A, Pinto GV. One-step reverse transcription polymerase chain reaction for semiquantitative analysis of mRNA expression. Methods Find exp Clin Pharmacol. 2002;24(5):253–9.View ArticlePubMedGoogle Scholar
- Zuckerkandl E, Pauling L. Evolutionary divergence and convergence in proteins. In: Bryson V, Vogel HJ, editors. Evolving genes and proteins. New York: Academic Press; 1965. p. 97–166.View ArticleGoogle Scholar
- Saitou N, Nei M. The Neighbor-Joining method - a new method for reconstructing Phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.View ArticlePubMedGoogle Scholar
- White JA, Beckett-Jones B, Guo Y-D, Dilworth FJ, Bonasoro J, Jones G, et al. cDNA cloning of human retinoic acid-metabolizing enzyme (hP450RAI) identifies a novel family of cytochromes P450 (CYP26). J Biol Chem. 1997;272(30):18538–41.Google Scholar
- Bowles J, Koopman P. Retinoic acid, meiosis and germ cell fate in mammals. Development. 2007;134(19):3401–11.View ArticlePubMedGoogle Scholar
- Smith CA, Roeszler KN, Bowles J, Koopman P, Sinclair AH. Onset of meiosis in the chicken embryo; evidence of a role for retinoic acid. BMC dev Biol. 2008;8:85.View ArticlePubMedPubMed CentralGoogle Scholar
- Wallacides A, Chesnel A, Chardard D. Flament S, Dumond H. Evidence for a conserved role of retinoic acid in urodele amphibian meiosis onset. Dev Dyn. 2009;238(6):1389–98.View ArticlePubMedGoogle Scholar
- Rodríguez-Marí A, Cañestro C, BreMiller RA, Catchen JM, Yan Y-L, Postlethwait JH. Retinoic acid metabolic genes, meiosis, and gonadal sex differentiation in zebrafish. PLoS One. 2013;8(9):e73951.View ArticlePubMedPubMed CentralGoogle Scholar
- Relógio A, Schwager C, Richter A, Ansorge W, Valcárcel J. Optimization of oligonucleotide-based DNA microarrays. Nucleic Acids res. 2002;30(11):e51.View ArticlePubMedPubMed CentralGoogle Scholar
- Miura T, Yamauchi K, Takahashi H, Nagahama Y. Hormonal induction of all stages of spermatogenesis in vitro in the male Japanese eel (Anguilla japonica). P Natl Acad Sci USA. 1991;88(13):5774–8.View ArticleGoogle Scholar
- Kobayashi M, Aida K, Stacey NE. Induction of testis development by implantation of 11-ketotestosterone in female goldfish. Zool Sci. 1991;8(2):389–93.Google Scholar
- Rolland AD, Lardenois A, Goupil A-S, Lareyre J-J, Houlgatte R, Chalmel F, et al. Profiling of androgen response in rainbow trout pubertal testis: relevance to male gonad development and spermatogenesis. PLoS One. 2013;8(1):e53302.Google Scholar
- Schultz D, Perez N, Tan CK, Mendez A, Capo T, Snodgrass D, et al. Concurrent levels of 11-ketotestosterone in fish surface mucus, muscle tissue and blood. J Appl Ichthyol. 2005;21(5):394–8.Google Scholar
- Schultz D, Perez N, Mendez A, Snodgrass D, Serafy J, Prince E, et al. Tracking gender factors in fish surface mucus: temporal patterns in individual koi (Cyprinus carpio). J Appl Ichthyol. 2007;23(2):184–8.Google Scholar
- Lara-Gonzalez P, Westhorpe FG, Taylor SS. The spindle assembly checkpoint. Curr Biol. 2012;22(22):R966–80.View ArticlePubMedGoogle Scholar
- Kline-Smith SL, Sandall S, Desai A. Kinetochore-spindle microtubule interactions during mitosis. Curr Opin Cell Biol. 2005;17(1):35–46.View ArticlePubMedGoogle Scholar
- Monje-Casas F, Prabhu VR, Lee BH, Boselli M, Amon A. Kinetochore orientation during meiosis is controlled by aurora B and the monopolin complex. Cell. 2007;128(3):477–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Mendenhall MD, Hodge AE. Regulation of Cdc28 cyclin-dependent protein kinase activity during the cell cycle of the yeast Saccharomyces cerevisiae. Microbiol Mol Biol R. 1998;62(4):1191–243.Google Scholar
- Slegtenhorst-Eegdeman K, Post M, Baarends W, Themmen A, Grootegoed J. Regulation of gene expression in Sertoli cells by follicle-stimulating hormone (FSH): cloning and characterization of LRPR1, a primary response gene encoding a leucine-rich protein. Mol Cell Endocrinol. 1995;108(1):115–24.View ArticlePubMedGoogle Scholar
- Li X, Schimenti JC. Mouse pachytene checkpoint 2 (trip13) is required for completing meiotic recombination but not synapsis. PLoS Genet. 2007;3(8):e130.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelman Z. PCNA: structure, functions and interactions. Oncogene. 1997;14(6):629–40.View ArticlePubMedGoogle Scholar
- Windley SP, Wilhelm D. Signaling pathways involved in mammalian sex determination and gonad development. Sex dev. 2015;9(6):297–315.View ArticlePubMedGoogle Scholar
- Josso N, Picard J, Rey R, Di Clemente N. Testicular anti-Müllerian hormone: history, genetics, regulation and clinical applications. Pediatric Endocrinol Rev. 2006;3(4):347–58.Google Scholar
- Pfennig F, Standke A, Gutzeit HO. The role of Amh signaling in teleost fish-multiple functions not restricted to the gonads. Gen Comp Endocrinol. 2015;223:87–107.View ArticlePubMedGoogle Scholar
- Miura T, Miura C, Konda Y, Yamauchi K. Spermatogenesis-preventing substance in Japanese eel. Development. 2002;129(11):2689–97.PubMedGoogle Scholar
- Wu G-C, Chiu P-C, Lyu Y-S, Chang C-F. The expression of amh and amhr2 is associated with the development of gonadal tissue and sex change in the protandrous black porgy, Acanthopagrus Schlegeli. Biol Reprod. 2010;83(3):443–53.View ArticlePubMedGoogle Scholar
- Rocha A, Zanuy S, Gómez A. Conserved anti-Müllerian hormone: anti-Müllerian hormone type-2 receptor specific interaction and intracellular signaling in Teleosts. Biol Reprod. 2016;87(2):35. 31-11Google Scholar
- Guiry A, Flynn D, Hubert S, O'Keeffe AM, LeProvost O, White SL, et al. Testes and brain gene expression in precocious male and adult maturing Atlantic salmon (Salmo salar). BMC Genomics. 2010;11(1):211.Google Scholar
- Blázquez M, Bosma PT, Chang JP, Docherty K, Trudeau VL. Gamma-Aminobutyric acid up-regulates the expression of a novel Secretogranin-II messenger ribonucleic acid in the goldfish pituitary. Endocrinol. 1998;139(12):4870–80.View ArticleGoogle Scholar
- Zhao E, Basak A, Wong AOL, Ko W, Chen A, López GC, et al. The Secretogranin II-derived peptide Secretoneurin stimulates luteinizing hormone secretion from Gonadotrophs. Endocrinol. 2009;150(5):2273–82.Google Scholar
- Zhao E, Hu H, Trudeau VL. Secretoneurin as a hormone regulator in the pituitary. Regul Peptides. 2010;165(1):117–22.View ArticleGoogle Scholar
- Blázquez M, Shennan KI. Basic mechanisms of secretion: sorting into the regulated secretory pathway. Biochem Cell Biol. 2000;78(3):181–91.View ArticlePubMedGoogle Scholar
- Fabra M, Raldúa D, Bozzo MG, Deen PM, Lubzens E, Cerdá J. Yolk proteolysis and aquaporin-1o play essential roles to regulate fish oocyte hydration during meiosis resumption. Dev Biol. 2006;295(1):250–62.View ArticlePubMedGoogle Scholar
- Zilli L, Schiavone R, Chauvigné F, Cerdá J, Storelli C, Vilella S. Evidence for the involvement of aquaporins in sperm motility activation of the teleost gilthead sea bream (Sparus aurata). Biol Reprod. 2009;81(5):880–8.View ArticlePubMedGoogle Scholar
- Boj M, Chauvigné F, Zapater C, Cerdá J. Gonadotropin-activated androgen-dependent and independent pathways regulate aquaporin expression during teleost (Sparus aurata) spermatogenesis. PLoS One. 2015;10(11):e0142512.View ArticlePubMedPubMed CentralGoogle Scholar
- Girardet C, Butler AA. Neural melanocortin receptors in obesity and related metabolic disorders. BBA- Mol Basis Dis. 2014;1842(3):482–94.View ArticleGoogle Scholar
- Guillot R, Cortés R, Navarro S, Mischitelli M, Garía-Herranz V, Sánchez E, et al. Behind melanocortin antagonist overexpression in the zebrafish brain: a behavioral and transcriptomic approach. Horm Behav. 2016;82:87–100.Google Scholar
- Volkoff H, Canosa L, Unniappan S, Cerda-Reverter J, Bernier N, Kelly S, et al. Neuropeptides and the control of food intake in fish. Gen Comp Endocrinol. 2005;142(1):3–19.Google Scholar
- Ilnytska O, Sözen MA, Dauterive R, Argyropoulos G. Control elements in the neighboring ATPase gene influence spatiotemporal expression of the human agouti-related protein. J Mol Biol. 2009;388(2):239–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Agulleiro MJ, Cortés R, Leal E, Ríos D, Sánchez E, Cerdá-Reverter JM. Characterization, tissue distribution and regulation by fasting of the agouti family of peptides in the sea bass (Dicentrarchus labrax). Gen Comp Endocrinol. 2014;205:251–9.View ArticlePubMedGoogle Scholar
- Reinecke M. Insulin-like growth factors and fish reproduction. Biol Reprod. 2010;82(4):656–61.View ArticlePubMedGoogle Scholar
- Sambroni E, Lareyre J-J, Le Gac F. Fsh controls gene expression in fish both independently of and through steroid mediation. PLoS One. 2013;8(10):e76684.View ArticlePubMedPubMed CentralGoogle Scholar
- Sambroni E, Rolland AD, Lareyre J-J, Le Gac F. Fsh and Lh have common and distinct effects on gene expression in rainbow trout testis. J Mol Endocrinol. 2013;50(1):1–18.View ArticlePubMedGoogle Scholar
- Lau E-L, Lee M-F, Chang C-F. Conserved sex-specific timing of meiotic initiation during sex differentiation in the protandrous black porgy Acanthopagrus schlegelii. Biol Reprod. 2013;88(6):150. 151-113View ArticlePubMedGoogle Scholar
- Schug TT, Berry DC, Shaw NS, Travis SN, Noy N. Dual transcriptional activities underlie opposing effects of retinoic acid on cell survival. Cell. 2007;129(4):723–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung S, Wolgemuth D. Role of retinoid signaling in the regulation of spermatogenesis. Cytogenet Genome Res. 2004;105(2–4):189–202.View ArticlePubMedPubMed CentralGoogle Scholar
- Pradhan A, Olsson PE. Inhibition of retinoic acid synthesis disrupts spermatogenesis and fecundity in zebrafish. Gen Comp Endocrinol. 2015;217-218:81–91.View ArticlePubMedGoogle Scholar
- Hu P, Tian M, Bao J, Xing G, Gu X, Gao X, et al. Retinoid regulation of the zebrafish cyp26a1 promoter. Dev Dyn. 2008;237(12):3798–808.Google Scholar
- Loudig O, Babichuk C, White J, Abu-Abed S, Mueller C, Petkovich M. Cytochrome P450RAI (CYP26) promoter: a distinct composite retinoic acid response element underlies the complex regulation of retinoic acid metabolism. Mol Endocrinol. 2000;14(9):1483–97.View ArticlePubMedGoogle Scholar
- Albalat R. The retinoic acid machinery in invertebrates: ancestral elements and vertebrate innovations. Mol Cell Endocrinol. 2009;313(1):23–35.View ArticlePubMedGoogle Scholar
- Mark M, Ghyselinck NB, Chambon P. Function of retinoid nuclear receptors: lessons from genetic and pharmacological dissections of the retinoic acid signaling pathway during mouse embryogenesis. Annu Rev Pharmacol Toxicol. 2006;46:451–80.View ArticlePubMedGoogle Scholar