Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii)
BMC Genomics volume 17, Article number: 217 (2016)
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.
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.
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.
Bluefin tunas (Thunnus spp.) provide some of the highest valued fish in the fresh and frozen international fish market and as such, they are highly susceptible to over-fishing, resulting in strict fishing regulations and quotas, which limit the available catch . In order to maintain a sustainable supply of bluefin tuna to meet the ever-growing demand, without seasonal or regional constrains, and relieve the fishing pressure from wild stocks, bluefin tuna supply should turn to aquaculture based production systems . To achieve this goal, bluefin tuna broodstock must be bred in captivity and therefore extensive research has been invested into facilitating the reproduction in captivity and broodstock management of three major bluefin tuna species: Pacific bluefin tuna (PBT, Thunnus orientalis), Atlantic bluefin tuna (BFT, Thunnus thynnus thynnus) and Southern bluefin tuna (SBT, Thunnus maccoyii).
Molecular techniques and tools are becoming increasingly helpful to identify key pathways and biological processes involved in bluefin tunas reproductive biology. Such molecular research includes the discovery of the Kisspeptin system in SBT , treatment of gonadotropin releasing hormone agonist (Gnrha) to induce spawning in BFT and SBT [4, 5], description of gene expression patterns in gonads of BFT [6, 7] and identification of germ cell specific molecular markers in PBT . Development of these molecular tools requires identification of genes and control elements coded in the DNA sequence and expressed as mRNA transcripts.
Germ cell transplantation (GCT) is an innovative molecular biotechnology for the production of inter-species surrogates, capable of facilitating easy handling, shorter generation time and a more economical management of large-bodied broodstock, such as the bluefin tuna. This technology involves the isolation of undifferentiated type-A spermatogonia (ASG) from a donor fish testis and transplanting these germ cell precursors into larvae of a host fish that will become a surrogate producing donor-derived gametes upon sexual maturation [9, 10]. The isolation of spermatogonia from the donor testis is achieved by mechanical and enzymatic dissociation of the testis, often followed by a second cell selection step using a discontinuous Percoll density gradient [11–13].
Implementation of GCT technology for a new species, however, requires development of molecular tools similar to those described above to identify transplantable germ cell populations in the donor gonads, to track the migration and settlement of the transplanted germ cells once in the surrogate host and to detect the presence of donor-derived DNA and RNA in the maturing hosts’ gonads and gametes [10, 14]. This approach was demonstrated previously in this study to distinguish between transplanted SBT germ cells and yellowtail kingfish (YTK, Seriola lalandi) hosts endogenous cells, using the 3’ untranslated region (UTR) of the germ cell-specific vasa gene. The transplanted SBT cells, however, did not proliferate and further differentiate in the YTK host, most likely due to molecular incompatibilities derived from the evolutionary distance between the two species . Moreover, because the host species should preferably be as phylogenetically close to the donor species, a higher level of homology is expected between the genes of the two species, therefore a wide range of molecular markers is needed to ensure that some would be divergent enough to allow for species-specific identification.
Isolation of genes in non-model organisms, such as the bluefin tuna species, has typically relied on degenerate-primer polymerase chain reaction (PCR) amplification of candidate genes, followed by sequencing. This method, however, is time consuming as it needs to be performed for each individual gene with substantial trial-and-error to clone the gene of interest. Furthermore, it requires prior knowledge of conserved regions of the candidate genes in other species, preferably as phylogenetically close as possible to the target species. This requirement presents a major bottleneck for gene discovery for the SBT, because there is a lack of gene sequences available in the public databases: currently less than 300 combined nucleotide and protein sequences and only 13 annotated genes can be found for SBT (taxonomy ID: 8240) in the National Center for Biotechnology Information (NCBI) databases . However, the recently published genome of the PBT , together with nucleotide and protein sequences from the entire Thunnus genus (taxonomy id: 8234) available on NCBI, provide good reference for comparative discovery of genes in the closely related SBT, to overcome the lack of publicly available sequence data.
Gene discovery methods are evolving rapidly from the traditional methods described above, to high throughput next generation sequencing (NGS), as a result of decreasing costs, fast processing times and a plethora of emerging analysis tools [18–21]. In a similar manner, gene expression data acquired with NGS RNA-Seq can cover the entire transcriptome of a sample in a single analysis, and serve as an alternative to individual gene real-time quantitative PCR and higher throughput microarrays, both of which require prior sequence knowledge .
This study aimed at identifying genes differentially expressed in male and female gonads of SBT. Special attention was given to genes known to be involved in germ cell differentiation and proliferation, to develop molecular tools for implementation of GCT for SBT. Specifically, markers for undifferentiated transplantable gonadal stem cells, ASG from testes and oogonia from ovaries, were sought after to enable their isolation and detection with molecular methods before and after transplantation. To achieve this, transcriptomes of SBT gonadal cells were assembled: crude cell extracts from ovary, testis, Percoll-enriched germ cells (as used in GCT experiments) and oogonia-enriched filtered cells. The transcriptomes were constructed using a combined approach of de novo and genome guided assembly of NGS RNA-Seq, with the PBT genome as a closely-related species reference. The transcriptomes were used to compare gene expression of SBT ovary and testis cells and evaluate currently practiced ASG and oogonia cell enrichment methods used to prepare germ cells for GCT. In addition, the discovery of reproduction-related genes and their expression profiles in male and female SBT gonads is described as a foundation for future use of molecular tools in SBT reproductive research towards successful broodstock management and spawning in captivity.
SBT gonads were collected during a commercial harvest at Cleanseas Tuna Ltd. sea-cages offshore Port Lincoln, South Australia. Twelve (12) fish were killed and immobilized by the harvest crew, then measured for weight (gilled and gutted), total length (TL) and examined for sex and gonad collection (details in Table 1). Similar sized fish were selected to try and minimize size and sexual development related variation between and within groups. The gonads were washed briefly with ice cold 10 mM phosphate buffered solution (PBS, P4117; Sigma-Aldrich) immediately after removal, transferred into 100 mL clean PBS and kept on ice during transport back to the laboratory, approximately 4 h after collection. Gonado-somatic index (GSI), presented as percentage, was calculated for each fish as the ratio between gonad weight to gilled and gutted weight (in grams).
All experimental procedures were performed in accordance with the University of the Sunshine Coast and Animal Ethics Committee guidelines (approval numbers AN/A/11/58), following the Australian Code for the Care and Use of Animals for Scientific Purposes .
The gonads were washed again in the laboratory in ice cold 10 mM PBS, weighed and sliced into 5 mm cross-sections. A whole cross-section slice from the middle part of each gonad was further cut into small pieces (5×5 mm) and stored in RNAlater (R0901; Sigma-Aldrich) at −20 °C until further processing. Testis sample from each male was further enzymatically dissociated, filtered and centrifuged through a discontinuous Percoll gradient (17-5445-01; GE Healthcare) to obtain ASG-enriched cell suspension, following the methods described in . Enriched ASG cell suspensions were stored in 1 mL of RNAlater at −20 °C until further processing. Ovary sample from each female was dissociated using the same methods described above, however, oogonia enrichment was performed without Percoll gradient, but rather by consecutive filtering of the dissociated cell suspension through sterile 150 and 30 μm filters (CellTrics 04-004-2326/9; Partec) respectively, as described in .
Cross-section slices from each gonad were fixed in 4 % paraformaldehyde for 24 h, then dehydrated in graduated ethanol washes (30, 50, 70, 90, 95 and 100 % for 2 h at 4 °C) and xylene (twice for 40 min at room temperature). The cross-sections were finally embedded in paraffin (x3 times for 40 min at 65 °C) and mounted in blocks. Tissues in the blocks were sliced to 7 μm sections using a microtome (UM-MS355; ProSciTech), and sections were then dried on superfrost slides (LBS4951; LabServ), stained with Harris hematoxylin and eosin (H&E) and permanently mounted using DePex (BDH Chemicals) .
Stained sections were visualized under a compound microscope (DM5500B; Leica Microsystems) and digital images were acquired using a digital camera (DFC550; Leica Microsystems) coupled with an image acquiring software (LAS v4.3; Leica Microsystems). Gonad maturation stage was assessed visually and fitted onto a categorical scale (ranging from 1–5 for testes and 1–9 for ovaries), based on the presence and morphology of germ cell populations in the gonads, as described in [26, 27].
RNA was extracted from five biological replicates (originating from different animals) of each group: enriched cells or crude extract from both testes and ovaries. Prior to RNA extraction, the samples were centrifuged at 4000 rpm for 3 min and the RNAlater supernatant was discarded by pipetting. Fifty milligram (50 mg) of crude gonad sample was placed in a 2 mL clean, RNase free micro-centrifuge tube, containing 1 mL RNAZol®; RT (RN 190; Molecular Research Center) and 5 ceramic (zirconium oxide) beads (2.8 mm, KT03961-1-002.2; Precellys). The tissue was homogenized using a Qiagen TissueLyser II (85300; Qiagen) at a frequency of 30 Hz, for 60 s, 2–6 times until no solid tissue was visible. The enriched cell suspensions were homogenized for 30 s using a motorized pellet pestle (Z359971; Sigma-Aldrich). Total RNA was then extracted according to the RNAZol® RT manufacturer instructions and stored at −80 °C. RNA quality and quantity were assessed visually using gel electrophoresis, then concentration was calculated from a spectrophotometer measurement (Nanodrop 2000; Thermo Scientific) and finally an RNA integrity number (RIN) was determined for each sample by a 2100 BioAnalyzer (G2947CA; Agilent).
Best quality RNA extracts were selected from for both crude and enriched extracts of 5 female and 5 male gonads (n = 20) with RIN >6, preferably using gonad sample from the same individual fish for both crude and enriched preparation when RNA quality allowed. The RNA samples were then sent to the Australian Genome Research Facility (AGRF; Brisbane, QLD) for library preparation and RNA-Seq on an Illumina HiSeq2500 platform, producing 100 base pair (bp) paired end reads. All five replicates from each treatment were multiplexed on each lane, to a total of 4 lanes, achieving coverage of (8–12) × 106 bp per sample.
Short read sequences obtained from RNA-Seq were downloaded and processed in-house at the University of the Sunshine Coast. Preliminary quality check of the reads was performed using FastQC (v0.11.2) and further 3’ and 5’ ends quality trimming and Illumina adapters removal were performed by Trimmomatic (v0.32) . The processed reads were mapped to a reference genome of the closely related PBT, T. orientalis , using BBMap spliced reads aligner (v33.65), with an adjusted seed length of k = 12 to account for cross-species mismatches. The resulting BAM files with the aligned reads were then used as an input for Trinity (r20140717) , to assemble the reads into a genome-guided de novo transcriptome of SBT gonad cells.
Gene and protein annotation
Open reading frames (ORFs) were predicted from the assembled transcripts using TransDecoder (r20140704); ORF was considered as complete by the presence of a starting methionine amino acid and an ending stop codon. Transcripts and predicted peptides were annotated by sequence alignment similarity search, BLAST (v2.2.30+) , to protein databases (NCBI nr, UniProt, Swiss-Prot and KOBAS, e-value <1×10−3) and by hidden Markov models protein domain identification using HMMER3.1  against the HMMER/Pfam protein database (v27.0) . Based on these annotations, Gene Ontology (GO)  and Kyoto Encyclopedia of Genes and Genomes (KEGG)  terms were assigned to each putative protein. Furthermore, prediction for presence of secretory signal peptide and trans-membrane topology were performed using SignalP (v4.1)  and TMHMM (v2.0) , respectively. The resulting annotation output files were further processed and cleaned to remove duplicates, select best matching annotation and identify errors or missing values.
The number of reads mapping to each transcript was estimated by RSEM . Preliminary exploratory data analysis (EDA) of the estimated counts of each sample was performed by computing and plotting a between-samples distance matrix and principal component analysis (PCA) to identify sample-related bias. The estimated counts were fed to edgeR  to identify statistically significant differentially expressed (DE) transcripts between treatment groups. Both tools were initiated automatically using scripts from the Trinity pipeline. Transcripts were considered as DE when the absolute value of the log base-2 of the fold change was greater than 2 (|log2FC>2|, positive or negative for either over- or under-expression respectively) with an adjusted p-value <1×10−3 to reduce the false discovery rate(FDR). When multiple transcripts were annotated to the same DE gene, log2FC was calculated as the mean with corresponding standard errors (SE).
Gene set enrichment analysis (GSEA) was performed for the GO and KEGG annotated DE genes in each comparison using the goseq R package . In brief, the total appearances of each GO and KEGG function was accumulated across DE genes, taking into account the length of each gene and determining enriched functions in each comparison group.
The annotation and expression data files were then combined and loaded onto a lightweight, standalone relational SQLite database (http://www.sqlite.org/), using the scripts provided in the Trinotate (v2.0.1, https://trinotate.github.io/) pipeline. This allowed fast and easy retrieval of sequences, annotation and expression data using any combination of conditional filtering and ordering. A complete bioinformatics data processing and analyses workflow is presented in Fig. 1.
Statistical analysis and further data summarizing was performed using the R (v3.1.2) statistical programming language . Specifically, relevant data was retrieved from the SQLite database with RSQLite and sqldf packages and then further processed using various tools from the Bioconductor R package repository (http://www.bioconductor.org/) .
Availability of supporting data
The data set containing the raw RNA-Seq reads generated in this study was deposited to the NCBI Sequence Read Archive database (study accession number: SRP059929), available at http://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA288188.
Results and discussion
Southern bluefin tuna gonad samples
Twelve SBT were sampled on May 27, 2014 from Cleanseas Tuna Ltd. sea-pens located in Port Lincoln bay, South Australia. Females mean body weight was 38.75 kg ± 340 kg with a mean GSI of 0.124 % ± 0.058 %; males mean body weight and GSI were 37.67 kg ±5.56 kg and 0.0378 % ± 0.0075 % respectively (n = 6 for each sex, Table 1). Histology of the sampled fish gonads revealed that the female ovaries show an asynchronous development pattern and were categorized at stages 1–3 (on a scale of 1–9 ) and comprised mainly of unyolked oocytes up to 50 μm in diameter, with a few early yolked oocytes up to 130 μm, in diameter in the more developed ovaries (Table 1 and Fig. 2 a–b). The males were at stages 2–3 (on a scale of 1–5 ), displaying all germ cell populations, from ASG to spermatozoa; however, the majority of the germ cells were at early stages of spermatogonia with increasing number of spermatocytes in stage 3 testes (Table 1 and Fig. 2 c–d). These body weights, GSI and gonad maturation stages are consistent with the previously reported range of values for early maturing male and sexually immature female SBT. For example, females have been reported as being sexually immature at GSI <0.22 % and gonad development stages below 4, while males start maturing at GSI <0.06 % and gonad development stage 3 .
Gonad sampling, preparation and RNA extraction workflow is presented in Fig. 3. Testis samples produced higher quality RNA following extraction, with average RIN of 7.8 for the crude extracts and 9.12 for the Percoll-enriched cell extracts, whereas ovary cells crude extracts RIN averaged at 7.1 (Table 2). A RIN for the filtered ovary cells extract could not be determined due to a low molecular weight peak present in the gel electrophoresis. The same artifact affected the crude ovary cells extract as well, resulting in the observed lower RIN numbers compared with testis cells extracts. All attempts to extract higher quality RNA without this artifact have failed.
Apparently, fish ovarian RNA contain a high peak of 5S ribosomal RNA, which changes the RNA profile on the gel and interferes with correct estimation of the 18S and 28S rRNA peaks, resulting in poor RIN numbers [42, 43]. Despite these abnormalities in the quality of the ovary cell RNA extracts, the RNA samples did pass the required quality checks at the RNA-Seq provider and were further processed and sequenced.
Southern bluefin tuna gonadal cells transcriptome assembly
A combined total of 1.55×109 short reads, 100 bp each, were sequenced by the Illumina Hiseq2500 RNA-Seq from all the samples and all treatments; of those reads, 95.4 % (1.48×109 reads) passed adapter removal, length and quality control (QC). The raw reads that passed QC were mapped to the genome of the closely related PBT (T. orientalis), giving an average mapping rate of 95.51 % and resulting in 1.34×109 mapped reads, an average of 6.72×107 reads/sample. The mapped reads were then assembled using a combined approach of genome-guided and de novo assembly, as detailed in Methods subsection ‘Assembly’, resulting in 128,065 contigs, with a median of N50 = 2094 bp and a longest contig of 20,234 bp long (Table 3 and Fig. 4).
When attempting to assemble the RNA-Seq transcriptome of a non-model species, whose genome is not available, two options for assembly can be considered: 1) A de novo approach, where transcripts are built solely from overlapping segments in the reads [29, 44–46]; or 2) A genome-based approach, where reads are aligned to a reference genome of a closely related species [47, 48]. A genome-based mapping can produce longer contigs and wider coverage thanks to the ability to close gaps of low read coverage. However, it uses a more lenient mapping strategy (smaller k-mer size compared with a de novo approach) to account for mismatches between the species, thus increasing the chance of introducing artifacts and sequencing errors . This trade-off must be considered when deciding which strategy to utilize, depending on the aim of the study. In this study, a combined approach was chosen, using the genome of the closely related PBT as a reference and mapping the reads to it followed by a de novo assembly while taking into account sequence dissimilarities and novel transcripts [29, 47, 49–51]. This strategy was selected to discover as many genes as possible that are involved in the reproduction and gonad development in the SBT, while using tools like homology-based annotation and coding region (CD) prediction to assist in quality filtering of the transcriptome. When higher fidelity is required, for example when the focus of the research is on discovery of isoforms and variant detection, a de novo approach with higher stringency setting might be more suitable .
Multiple ORFs were predicted from each transcript, to a total of 136,465 ORFs. The ORF list was further processed to remove redundancies (ORFs sharing exactly the same peptide sequence) and then annotated to known proteins. For each ORF, the best hit annotation was selected, discarding ORFs that had no hit within any of the protein databases used, to provide a final list of 84,204 annotated ORFs from 300 different species (Table 4). Only 145 out of the 84,204 ORFs (less than 0.2 %) were annotated back to species of the Thunnus genus (highlighted in Fig. 5). The low number of ORFs annotated to Thunnus genus emphasizes the lack of publicly available annotated tuna gene data. Currently, 6754 nucleotide and 2120 protein sequences from the entire Thunnus genus (taxonomy id: 8234) are available on NCBI databases. These include mainly mitochondrial DNA used in phylogenetic and energy metabolism studies (approximately 50 % of the sequences), as well as genes involved in the optic, locomotion and reproductive systems, the main functional systems studied in the tunas [17, 52–55]. These sequences contain a high level of redundancy, with multiple isoforms and partial sequences for each gene.
Differential expression exploratory data analysis
Differentially expressed (DE) gene analysis was performed on the estimated counts of the reads mapped to the CD of predicted ORFs of each transcript; the count estimation was performed by RSEM software.
Prior to the DE gene analysis, an EDA was performed to detect major factors or biases affecting the variability between the different analysis groups and samples. A distance matrix was computed between each sample/treatment combination and similar combinations were clustered together, as presented by the heatmap and dendrogram in Fig. 6 a. The samples from each replicate cluster together by the tissue type and treatment: testis cells, Percoll-enriched testis cells, ovary cells and enriched ovary cells. The greatest variation was found between testis cells and ovary cells, both crude and enriched (green squares in top left and bottom right of Fig. 6 a). Medium variation was found between the Percoll-enriched testis cells and the crude testis cells extracts (black squares in top right of Fig. 6 a). The least variation was found between enriched ovary cells and crude ovary cells extracts (dark red squares in bottom left of Fig. 6 a).
Principal component analysis (PCA) was performed on the regularized log-transformed counts, using tissue and treatment (crude/enriched) as factors (Fig. 6 b). The tissue factor was found to contribute to 90 % of the variance, whereas the treatment factor for only 2 %. The testis cell samples, both crude and enriched extracts grouped tightly together, with a small but distinct difference between the crude and the Percoll-enriched cells. The ovary cell samples, however, were widely scattered with less obvious separation between the 2 treatments.
A second PCA was performed, this time for each tissue separately and including the individual fish which the samples were taken from as a factor. The testis cells showed a similar pattern as observed in the first analysis, with clear separation between crude and enriched cell samples, meaning that the enrichment treatment contributed to 61 % of the variance, with no significant effect of the individual fish (Fig. 6 c). The ovary cells exhibited some separation by treatment (x-axis in Fig. 6 d); however, the individual fish that each sample originated from had a significant contribution to the variance, at the same level as the treatment effect (y-axis in Fig. 6 d).
The results of the EDA suggest that the sampled male fish were more uniform in their gonadal cell expression profile (Fig. 6 c) and shown a significant effect caused by the Percoll enrichment treatment, whereas the sampled females exhibited an individual-dependent expression profile (Fig. 6 d), with little effect caused by the filtration enrichment treatment. This suggests that the Percoll enrichment protocol used for the testis cells was more effective and produced a significant difference in the composition of the resultant cell suspension, whereas the female ovary cell filtration had little (if any) germ cell enrichment effect. An alternative explanation might be that the sampled female fish had inherently greater gonadal development variability between individuals than that of the sampled males. The GSI distribution of the sampled females supports the latter, with a standard error which accounts for 46 % of the mean (0.124 % ± 0.058 %), while in the males, the GSI standard error was less than 20 % of the mean (0.0378 % ± 0.0075 %), as detailed in Table 1. Both male and female GSI values were below the reported minimum for sexually mature SBT . High variability within groups in gene expression data of fish gonads was reported in similar RNA-Seq transcriptome analysis of the sharpsnout seabream (Diplodus puntazzo) .
The filtration method used to enrich undifferentiated oogonia cells was based on physical dimensions of the pre-vitellogenic oogonia, which are the smallest germ cells in the fish ovary [24, 56]. This method was used after a failed attempt to obtain sufficient numbers of oogonia when using the same Percoll enrichment method that was used to enrich spermatogonia from testis. However, the ovary cell count data EDA demonstrated high variability within treatment groups and non-significant effect of the filtration enrichment treatment, which contributed to just over 20 % of the count data variance (according to the PCA of the ovary cell samples, Fig. 6 d). Based on these results, it was concluded that any DE transcripts found in this comparison could not be accounted for the enrichment procedure and therefore this comparison was discarded from the subsequent DE transcript analysis.
These results, along with the lower RNA quality of the enriched ovary cell extractions (Table 2), suggest that an improved method for oogonia separation is needed, both for GCT and for comparative gene expression purposes.
DE genes in crude ovary vs. testis cells
DE transcript analysis was performed on a pairwise comparison (contrasts) of crude ovary cells vs. crude testis cells. The resulting DE transcripts of each analysis were filtered to include only annotated ORFs that passed the thresholds of |log2FC >2| (positive or negative for either over- or under-expression, respectively) and p-adjusted-value <1×103. The average rate of annotated ORFs out of the total DE ORFs in this comparison was 81.2 % and these will be referred to from now on as “genes”. The number of DE genes in each tissue/treatment combination within each analysis and the ratio between annotated and non-annotated ones are presented in Table 5 and Fig. 7 a.
Altogether, 17,849 annotated ORFs were found to be DE when comparing ovary to testis cells; 10,823 of them over-expressed in the testis and just over 7000 were over-expressed in the ovary cells. The total number of DE genes and the bias towards over-expression in testis match reported results of DE gene analysis in gonads of Nile tilapia (Oreochromis niloticus)  zebrafish (Danio rerio) [58, 59] and sharpsnout seabream (D. puntazzo) , where similar differences were found. According to these studies, this bias towards male transcripts had been suggested to occur as a result of evolutionary pressure, driven by male competitiveness.
Each gene was further annotated by assigning GO terms associated with its UniProt and Swiss-Prot annotation. Similar GO molecular function terms were then clustered together and a DE GSEA was performed to determine enriched terms in each analysis set (Fig. 8, full list available in Additional file 1). Enriched GO molecular functions with a q-value <0.01 in the ovary vs. testis cells DE analysis are presented in Fig. 8, separated by the tissue in which each term is over-expressed. Review of these DE enriched GO terms reveals that the majority are general cell maintenance molecular functions related to metabolism and protein translation and synthesis, such as nucleic acid and RNA binding and ATPase activity. In contrast to the ovary, testis cells contain an abundance of genes with motility-related molecular functions.
An in silico enrichment of specific germ cell and gonadal development-related genes and markers was used to refine the DE analysis results. An in-depth analysis of a subset of germ cell and gonadal development-related genes revealed additional expression patterns (selected genes are presented in Fig. 9 and a full list in Additional file 2) as detailed below.
Steroid synthesis related genes
Tissue-specific (and consequently sex-specific) DE steroid synthesis genes are notable: 3- β-hydroxysteroid dehydrogenase (3b-hsd), expressed in Leydig cells and responsible for biosynthesis of steroids, mainly progesterone and testosterone , was highly over-expressed in testis cells. On the other hand, 17- β-hydroxysteroid dehydrogenase 8 (17b-hsd8) and 17b-hsd12, that regulate the concentration of biologically active estrogens and androgens and are involved in synthesis of estradiol (E2) , were both over-expressed in ovary cells.
Corticosteroids and their receptors play an important role in regulating reproductive functions by acting as transcription factors in the gonad somatic cells . A member of the corticosteroids synthesis pathway, which was found to be over-expressed in high levels in testis cells, is 11- β-hydroxysteroid dehydrogenase 2 (11b2-hsd), an enzyme which converts gonad maturation inhibiting cortisol or corticosterone into inactive cortisone or 11-dihydrocorticosterone respectively [62, 63].
Other genes involved in the synthesis of steroids from cholesterol are the P450 cytochrome enzymes. The 17- α-hydroxylase, 17, 20 lyase (cyp17a1) was over-expressed in the testis cells, whereas aromatase (cyp19a) was highly over-expressed in ovarian cells. cyp17a1 and cyp19a were shown to express both in germ cells and Leydig cells in zebrafish (D. rerio) testis, with high expression levels of cyp17a1 in Leydig cells. However, exposure to exogenous steroids, such as E2, can dramatically lower the expression levels of cyp17a1 in Leydig cells and lead to a disruption of spermatogenesis .
A comparison between XX female to XY male Nile tilapia (O. niloticus) gonads found similar expression patterns to the ones detailed above for 3b-hsd, 17b-hsd8, 17b-hsd12a1, 11b2-hsd, cyp17a1 and cyp19a [65, 66]. The cyp2k1, which was highly over-expressed in the ovary cells, is known for its lauric acid (ω-1) and (ω-2), E2 and progesterone hydroxylation activities [67, 68]. It had also been suggested as a biomarker to determine prior exposure to produced formation water, that is, water that are discharged from oil-well platforms and produced as a by-product along with the oil . Interestingly, these oil-well platforms are located at the North-West Shelf of Australia, on the migration route of SBT, between the only reported spawning grounds in South-East of Java, Indonesia  and their site of capture for ranching in the Great Australian Bight in South Australia.
The above mentioned steroid synthesis genes demonstrate an expected expression pattern, with androgen synthesis enzymes over-expressed in the testis cells, while estrogen synthesis enzymes over-expressed in the ovary cells. The function of cyp2k1 has not been well defined in the context of sex steroids synthesis, but rather in the metabolism of carcinogenic substances, and therefore its high level of over-expression in the ovary cells requires further investigation.
A second set of genes that were DE between ovary and testis cells are genes associated with sex determination and sex-specific gonad development. One such gene that was highly expressed in the testis cells is the DNA binding motif domain Double-sex and mab-3 related transcription factor 1 (dmrt1). The dmrt1 is known to be involved in a male-specific vertebrate sex-determining pathway and is the key sex-determining factor in medaka (Oryzias latipes) [71–73], with similar function suggested in tilapia (O. niloticus)  and rainbow trout (Oncorhynchus mykiss) . In contrast, dmrt5, a less studied member of the dmrt family, was found to be over-expressed in the ovary cells. In zebrafish (D. rerio), dmrt5 was shown to be involved in the endocrine control of the reproductive system by controlling corticotrope and gonadotrope differentiation in the pituitary . In zebrafish gonads, dmrt5 expression was limited to developing germ cells, especially spermatogonia, spermatocytes, spermatids and sperm cells in the testis; and developing oocytes, including early perinucleolus stage oocyte, late yolk vesicle stage oocyte, and oil drop stage oocyte in the ovary . Its expression pattern in the SBT therefore could be indicative of the developmental stage of the gonad.
The SRY-related HMG-box 3 (sox3) transcription factor, is often mentioned along with dmrt1 as a male-specific gene, involved in early stages of testis development and growth, mainly in mammals [77, 78], but also in sex-changing protandrous (Acanthopagrus schlegeli) and protogynous (Epinephelus coioides) fish [79, 80]. In SBT however, sox3 was found highly over-expressed (log2FC=5.8) in the ovary cells, indicating it might have additional roles in ovarian differentiation not previously recognized. Another Sox transcription factor gene, sox30, had an opposing expression pattern and was moderately over-expressed (log2FC=4.6) in testis cells. Just recently identified in non-mammalian vertebrates, sox30 was found to be expressed exclusively in germ cells within Nile tilapia (O. niloticus) testis, indicating its possible involvement in spermatogonial differentiation and spermatogenesis, like in human and mouse [81, 82].
The wingless-type MMTV integration site family member 4 (wnt4) gene, was moderately over-expressed in the ovary cells (log2FC=5.6). Wnt4 is one of the few genes that were identified as ovarian-determining genes in mammals; it regulates the development of the female reproductive tract by suppressing male features development in XX gonads [83–85]. In fish, wnt4 has been reported to be important to oocyte maturation and is expressed in the developing ovary of both protandrous (A. schlegeli) and protogynous (E. coioides) fish [86, 87]. Similar to wnt4, the forkhead-box L2 (foxl2) transcription factor, which was highly over-expressed (log2FC=6.9) in the SBT ovary cells, has been associated with female gonad differentiation through its function as a transcriptional regulator of aromatase (cyp19a), a function that is conserved across a variety of vertebrates, including mammals and fish [84, 88]. In fish, foxl2 was shown to be expressed exclusively in somatic cells of developing (pre-vitellogenic) ovaries of gonochoristic fish such as tilapia [66, 89], medaka (O. latipes and Oryzias luzonensis) [90, 91], Japanese flounder (Paralichthys olivaceus)  and Rare minnow (Gobiocypris rarus)  to name a few. In sex-reversing hermaphrodite fish, however, foxl2 expression pattern was not as strictly defined: it was still highly expressed in early ovarian developmental stages gonads of protogynous species (E. coioides, Halichoeres trimaculatus and Monopterus albus) [79, 87, 94, 95] but its expression levels declined in transition to testis, and vice versa in protandrous species (A. schlegeli and D. puntazzo) [43, 86].
Endocrine regulation-related genes
Another set of genes that were found to be over-expressed in testis cells compared to ovary cells included various gonadotropin pathway receptors: luteinizing hormone (lh) receptor, follicle stimulating hormone (fsh) receptor and gonadotropin releasing hormone (gnrh) receptor. These genes were moderately over-expressed in the testis (log2FC up to 5), along with high over-expression (log2FC=6.87) of the kisspeptin1 receptor (kiss1r). This synchronous expression pattern was previously reported during testis germ cell maturation in the fathead minnow (Pimephales promelas)  and the Scombrid chub mackerel (Scomber japonicus) [97–100]. The platelet-derived growth factor (pdgf) and its receptor (pdgf receptor b) were both moderately over-expressed (log2FC between 3 and 4) in the testis cells. The Pdgf signaling system, through its several isoforms and receptors, has been shown to play a critical role in regulating the development and functional control of various tissues in humans and mice [101, 102]. In testicular development, pdgf was shown to be expressed by Sertoli cells to promote proliferation and/or differentiation of Leydig cells, and is required for intra-gonad cell migration, both in fetal and adult testis [83, 103].
In rainbow trout (O. mykiss) developing embryos, pdgf receptor expression was found to be in tight association with dmrt1 and therefore was suggested as a male sex-determining gene .
Another genes of interest are tachykinin3, which encodes a neuropeptide that plays an important role in neuroendocrine regulation of reproduction and its associate kisspeptin, known as the puberty onset gene in mammals, where it controls pulsatile gonadotropin-releasing hormone release [3, 104]. Both tachykinin3 and kisspeptin were over-expressed in the SBT ovary cells. While tachykinin3 displayed some of the highest over-expression (log2FC=11.88) among the gonad development-related genes in the ovary, kisspeptin was moderately over-expressed (log2FC=4.33) in an interesting contrast to the high over-expression of its receptor in the testis.
A member of the Notch protein family, notch2, was slightly over-expressed in the ovary cells. Involvement of Notch2 protein in ovarian development has not been thoroughly studied in fish so far; in mouse however, it was recently suggested that Notch2 cooperates with gonadotropins to regulate ovarian follicle formation and growth [105, 106]. This finding in SBT suggests it might have a similar role in fish.
Germ cell-specific genes
A final set of genes that were DE in this comparison are germ cell-specific genes. The lymphocyte antigen 75 (ly75), glial cell line-derived neurotrophic factor (GDNF) family receptor alpha 1 (gfra1), nanos2, notch1 and promyelocytic leukaemia zinc finger (plzf), which were previously reported as ASG-specific markers in GCT studies , were over-expressed in the testis at moderate levels (log2FC between 2.5 and 5).
The POU domain, class 5, transcription factor 1 (pou5f1), also known as octamer-binding transcription factor 4 (oct4) and nanos1 are key genes in stem cell maintenance and regulation [108–110]. Both were over-expressed in the ovary, with very high expression levels of pou5f1 (log2FC=9.15), indicating strong stem cell activity in the female ovaries.
The aquaporin10 gene, encoding for a water (and other solubles) channel, was found at very high levels at the ovary cells (log2FC=11), an approximate x2000 fold change compared with its expression in testis cells. Similar expression pattern (though with lower fold change) was reported in a comparison of gene expression levels between mature ovaries and testis of BFT (T. thynnus thynnus) . Aquaporin10 has been reported to facilitate oocyte hydration in the maturation of buoyant (pelagic) eggs of the gilthead sea bream (Sparus aurata); though oocyte hydration occurs at late stages of oocyte maturation, aquaporin10 was highly expressed in early vitellogenic ovaries [111, 112].
Crude ovary vs. testis cells DE summary
The gonad development gene expression patterns detailed above suggest that the ovaries of the sampled female fish were at an early stage of gonadal development, mainly expressing high levels of transcripts for stem cells maintenance (nanos1 and pou5f1/oct4) and early stages of oocyte development-related genes (aquaporin10, foxl2 and dmrt5). In contrast, genes identified as over-expressed in testis cells, encode for proteins whose function indicates a more advanced gonad developmental stage, with increased expression of receptors for reproductive hormones (fsh, lh, gnrh and kisspeptin receptors) as well as enzymes involved in sex steroids and hormone synthesis and signaling.
These findings are supported by the previously determined minimal weights for sexual maturation in SBT, in which males display sexual maturation earlier than females, at a minimal weight of 50 kg and GSI >0.065%, while females exhibit sexual maturation signs only from a weight of 100 kg and GSI >0.71 % . The sexual maturation process in SBT males progresses at a faster rate than in females, with all stages of germ cells (ASG, type-B spermatogonia - BSG, spermatocytes, spermatids and spermatozoa) observed already in immature males in the second stage of spermatogenesis (observed in males weighing 30–71 kg) . The fish that were sampled in this study, in comparison, weighed 38 kg on average, with the heaviest male weighing just under 50 kg, with a GSI =0.046 %, determined to be at stage 3 (out of 5) of spermatogenesis (early maturing testis) . The heaviest female weighed 42.2 kg, with a GSI =0.22 % and determined to be at stage 3 (out of 9) of gonad development (immature ovary, details in Table 1 and Fig. 2) . These results confirm that the males were approaching minimal size for sexual maturation, whereas the females were still in early stages of gonadal development and in preparation for puberty.
DE genes in crude vs. Percoll-enriched testis cells
DE gene analysis was performed on a pairwise comparison (contrasts) of crude testis cells vs. Percoll-enriched testis cells. The resulting DE transcripts of each analysis were filtered as described for the crude testis vs. crude ovary (subsection ‘DE genes in crude ovary vs. testis cells’). The average rate of annotated ORF out of the total DE genes in this comparison was 82 %. The number of DE genes in each tissue/treatment combination within each analysis and the ratio between annotated and non-annotated ones are presented in Table 5 and Fig. 7 b. The genes were further annotated by GO terms and then analyzed as previously described.
GSEA of GO molecular functions associated with DE genes in crude testis compared with Percoll-enriched testis cells was performed. The genes that were found over-expressed in Percoll-enriched cells are mainly involved in transcript translation (RNA polymerase and DNA binding) and cell motility (microtubule, ion channels and ATPase activities) (Fig. 7 b, full list available in Additional file 3). Similar enrichment of these functions was observed in microarray gene expression of rainbow trout (O. mykiss) testis cells; specifically, enrichment of genes related to RNA metabolism and DNA replication functions was associated with advanced-stage, proliferating spermatogonia, whereas cell motility functions were enriched in post-meiotic testicular germ cells . Within the molecular functions enriched in the crude testis cells, a few functions related to germ cell activity stand out: G-protein coupled receptors, growth factor, chemokine receptor, tumor necrosis factor receptor, platelet-derived growth factor receptor and interleukin 1 binding activities (highlighted in Fig. 7 b). These molecular functions were reported to be associated with somatic cells regulation of spermatogonia maintenance and differentiation [107, 113, 114].
An in silico enrichment of gonad development and germ cell-related genes was performed to reveal specific expression patterns out of the entire DE genes in the crude vs. Percoll-enriched testis cells (selected genes are presented in Fig. 8 b and a full list in Additional file 4).
Germ cell-specific genes
A set of germ cell-specific genes was found to be over-expressed at the crude testis cells when compared to the Percoll-enriched cells (Fig. 8 b). These include the genes ly75, plzf, notch1, nanos1, gfra1, chemokine receptor 4 (cxcr4), platelet-derived growth factor receptor (pdgf receptor) and Wilm’s tumor suppressor 1 (wt1), all previously reported to be expressed in ASG and involved in germ cells maintenance [107, 115]. In human, WT1 encodes a zinc finger transcription factor and RNA-binding protein that can activate or repress numerous target genes, including PDGF, thus playing a key role in testis development and maintenance [115, 116]. The G-protein-coupled chemokine receptor Cxcr4 is expressed in germ cells and directs their migration to the genital ridge in the developing embryo through detection and binding of its ligand, Sdf1 (Cxcl12) [117, 118]. The Cxcr4–Sdf1 signaling system has been found to be highly conserved across fish and other vertebrate species, enabling the migration of transplanted germ cells to the genital ridge of host larvae in cross-species GCT experiments .
In contrast, importins 4 and 5 (ipo4/5) were slightly over-expressed (log2FC of 3 and 2.4, respectively) in the Percoll-enriched testis cells. Importin proteins are nucleo-cytoplasmic transport proteins that play a central role in cargo-specific protein transport through the nuclear pores and into the nucleus, thus tightly regulating access of transcription factors to the genomic DNA. In mammalian testis, Ipo4 and Ipo5, both importin- β family members, were up-regulated in advanced stages of post-meiosis spermatogenesis, i.e. in spermatocytes and spermatids; and their function was suggested to be involved in chromatin remodeling [120–122].
Slightly over-expressed in the Percoll-enriched cells (log2FC=2.7), is the Ago3 Piwi domain, a protein of the Argonaute (Ago) protein family, which is required for the self-renewal of germ-line stem cells as it maintains transposon silencing in the germ-line genome through its interaction with small Piwi-associated RNA (piRNA) [123–125]. The proliferating cell nuclear antigen (pcna), which was also over-expressed in the Percoll-enriched testis cells, plays an important role in BSG proliferation initiation before entering meiosis . Pcna was used as a molecular marker for successful proliferation of transplanted yellowtail (Seriola quinqueradiata) spermatogonia after GCT into nibe croaker (Nibea mitsukurii) surrogate hosts .
Somatic cell-specific genes
Other genes that demonstrated over-expression in the crude testis cells are genes that are known to be expressed in somatic cells, supporting the developing germ cells: anti müllerian hormone and its receptor (amh/amhr2), lh and fsh receptors, gonadal soma-derived factor 1 (gsdf1), insulin growth factor 1 receptor (igf1r), sox9 and pdgf (Fig. 8 b). Gsdf1 was shown to localize predominately in somatic cells surrounding spermatogonia in the testis of medaka (O. latipes) and rainbow trout (O. mykiss) and plays an essential role in spermatogonia differentiation [128, 129]. Expressed by Sertoli cells, platelet-derived endothelial cell growth factor (PD-ECGF) and members of the transforming growth factor (TGF-b) superfamily of glycoproteins (activins, inhibins and amh), mediate spermatogonia proliferation and differentiation under control of the sex steroids E2 and 11-Ketotestosterone (11KT) (reviewed in [114, 130]). The insulin-like growth factor 1 (Igf1) plays a key role in regulation of bony fish growth, differentiation and reproduction and was shown to stimulate proliferation of testicular germ cells in vivo [131, 132]. It appears that Igf1 can directly stimulate germ cells, thanks to the localization of its receptor (Igf1r) both in Sertoli cells and germ cells (i.e., ASG) [133, 134]. The igf1r was moderately over-expressed (log2FC=4.1) in SBT crude testis cells.
Unlike the mammalian Sox9 gene, which is considered a male-determining gene , the medaka (O. latipes) and zebrafish (D. rerio) sox9b is not necessary for testis differentiation, but it is required for germ cell proliferation and survival. In medaka, somatic cells supporting undifferentiated germ cells, organized in units named ‘germinal cradles’, were shown to express sox9b predominately in males (Sertoli cells surrounding spermatogonia), but also in females (granulosa cells surrounding oogonia) .
Crude testis vs. Percoll-enriched testis cells DE summary
Both analyses, the GO molecular function GSEA and the specific gonad development DE genes, found consistent differential gene expression patterns between the Percoll-enriched and the crude testis cells. The genes found to be over-expressed at the Percoll-enriched testis cells, compared with crude testis cells, are known to be expressed exclusively in germ cells: ago3 piwi domain, pcna, ipo4, ipo5 and sox30. This expression pattern, together with over-expression of somatic cell specific genes (such as the Leydig cells-expressed 3b-hsd) in the crude testis cells, demonstrates that the Percoll enrichment method successfully enriched testicular germ cell populations over somatic cells. However, most well-known ASG-specific markers (ly75, plzf, notch1, nanos1 and gfra1) were found to be slightly over-expressed in the crude testis cells rather than in the Percoll-enriched cells, suggesting that the Percoll-enrichment protocol failed to enrich ASG. The over-expression of pcna, ipo4 and ipo5 in the Percoll-enriched cells indicates the enrichment of differentiated testicular germ cells that are towards, or post meiosis, such as proliferating BSG and post-meiosis spermatocytes and spermatids, in the cell suspension.
The Percoll-enrichment method used in this study did not appear to yield high proportions of enriched ASG compared with crude testis cells, yet this enrichment method is known to obtain transplantable ASG, capable of migrating and colonizing the host gonads, as demonstrated both by this group  and in previous studies [11, 127, 136]. A single pluripotent germ cell was demonstrated to be sufficient for successful GCT and therefore even a lowly ASG-enriched cell suspension is suitable for the purpose of GCT . The ASG enrichment protocol can be optimized by experimentally changing the Percoll concentration in each layer of the gradient, as previously suggested in [12, 138], as well as by dissecting just the peripheral, ASG-rich part of the testis, as demonstrated in . Such improvements could achieve better enrichment of transplantable ASG, which has the potential to increase the rate of successful migration and colonization of the transplanted ASG in the surrogate host gonads .
Germ cell separation methods for DE analysis
For the purpose of differential gene expression and cell type-specific gene discovery, which require better separation between germ cell populations in ovary and testis cells, more stringent cell selection methods might be required. Such methods include single-cell laser capture micro-dissection , or specific green fluorescence protein (GFP) labeling of the cells, followed by fluorescence activated cell sorting (FACS) [141–144]. It is worth noting though, that the use of GFP-labeled germ cells for FACS requires molecular engineering methods, which are not always applicable for non-model species lacking available gene data, such as the SBT.
Additional annotated germ cell-specific genes
Additional germ cell-specific genes that did not pass the threshold for significantly DE genes were identified from the annotated SBT gonadal cells transcriptome. These include vasa and dead-end, both used extensively as germ cell markers in fish gonad development and GCT studies [13, 15, 145–149]. These studies have shown that vasa is expressed in all germ cell populations, starting from the embryonic primordial germ cells (PGC), through all levels of gametogenesis [140, 150–152], while dead-end was found more specific to undifferentiated germ cell populations (PGC, ASG and oogonia) [15, 107, 139]. The presence of these genes in germ cell populations in both ovaries and testes explains why they have not been detected as significantly DE.
Discovering the full sequence of these genes in SBT enabled the development of SBT-specific germ cell markers that were used successfully for both in situ hybridization and PCR analyses of transplanted SBT spermatogonial cells in a GCT study performed by this group .
A high throughput analysis of gene expression in SBT gonad cells was performed in this study. The majority of genes described in this study were sequenced for the first time in SBT. The wealth of SBT gonadal and germ cell-related gene sequences identified (complete list in Additional files 2 and 4) provide a valuable resource for further reproductive molecular biology studies of the commercially important SBT and stands out in comparison to the low number of annotated tuna (Thunnus spp.) genes currently available. SBT ovaries and testes were found to have clearly distinct gene expression patterns, with thousands of DE genes involved in tissue structure, metabolism, cellular processes and gonad development. In silico enrichment of genes that are specific to, or play a major role in gonad development and germ cell molecular functions allowed for a more focused investigation of their DE patterns. The expression patterns of these gonadal genes supported the observed differences in sexual maturation stages between the sampled male and female fish and expand the existing molecular knowledge of sex differentiation and sexual maturation in SBT. Based on the differential gene expression between crude testis and Percoll gradient germ cell-enriched fraction, we suggest that the Percoll enrichment method, which is currently used in GCT experiments, should be adjusted and tested for the specific study species.
Furthermore, a shortlist of genes associated with specific germ cell activity were identified (Table 6). These may be used as SBT-specific germ cell molecular markers to distinguish between transplanted SBT germ cells and the surrogate host endogenous germ cells in future GCT studies. These genes could also be used in molecular studies of germ and soma cell interactions, in an attempt to unveil the molecular mechanisms governing germ cell proliferation and maturation to try and artificially recover these functions in naturally non-compatible GCT surrogate hosts.
FAO. The State of World Fisheries and Aquaculture 2014 (SOFIA). http://www.fao.org/publications/card/en/c/097d8007-49a4-4d65-88cd-fcaf6a969776. Accessed 2015-04-24.
Gándara F. The Challenge of Bluefin Tuna Domestication. The SELFDOTT project. Spain. 2012. http://cordis.europa.eu/documents/documentlibrary/117900771EN6.pdf. Accessed 6 Jul 2015.
Nocillado JN, Biran J, Lee YY, Levavi-Sivan B, Mechaly AS, Zohar Y, et al. The Kiss2 receptor (Kiss2r) gene in Southern Bluefin Tuna, Thunnus maccoyii and in Yellowtail Kingfish, Seriola lalandi - functional analysis and isolation of transcript variants. Mol Cell Endocrinol. 2012; 362(1-2):211–20. [doi:10.1016/j.mce.2012.06.024, 22824208.
Mylonas CC, Fostier A, Zanuy S. Broodstock management and hormonal manipulations of fish reproduction. Gen Comp Endocrinol; 165(3):516–34. [doi:10.1016/j.ygcen.2009.03.007].
Knibb W, Miller A, Deichmann M, Lamont R, Zohar Y, Foster C, et al. Statistical multivariate analysis to identify potential spawning cues for captive southern bluefin tuna, Thunnus maccoyii. World Aquaculture. 2014.
Chini V, Cattaneo AG, Rossi F, Bernardini G, Terova G, Saroglia M, et al. Genes expressed in Blue Fin Tuna (Thunnus thynnus) liver and gonads. Gene; 410(1):207–13. [doi:10.1016/j.gene.2007.12.012].
Gardner LD, Jayasundara N, Castilho PC, Block B. Microarray gene expression profiles from mature gonad tissues of Atlantic bluefin tuna, Thunnus thynnus in the Gulf of Mexico. BMC Genomics; 13(1):530. [doi:10.1186/1471-2164-13-530.23036107. Accessed 2015-04-11].
Nagasawa K, Miwa M, Yazawa R, Morita T, Takeuchi Y, Yoshizaki G. Characterization of lymphocyte antigen 75 (Ly75/CD205) as a potential cell-surface marker on spermatogonia in Pacific bluefin tuna thunnus orientalis. Fish Sci; 78(4):791–800. [doi:10.1007/s12562-012-0501-9].
Takeuchi Y, Yoshizaki G, Takeuchi T. Biotechnology: Surrogate broodstock produces salmonids. Nature; 430(7000):629–30.
Lacerda S, Costa G, Campos-Junior P, Segatelli T, Yazawa R, Takeuchi Y, et al. Germ cell transplantation as a potential biotechnological approach to fish reproduction. Fish Physiol Biochem.:1–9. [doi:10.1007/s10695-012-9606-4].
Takeuchi Y, Higuchi K, Yatabe T, Miwa M, Yoshizaki G. Development of Spermatogonial Cell Transplantation in Nibe Croaker, Nibea mitsukurii (Perciformes, Sciaenidae). Biol Reprod; 81(6):1055–63. [doi:10.1095/biolreprod.109.077701].
Heidari B, Gifani M, Shirazi A, Zarnani A-H, Baradaran B, Naderi MM, et al. Enrichment of Undifferentiated Type A Spermatogonia from Goat Testis Using Discontinuous Percoll Density Gradient and Differential Plating. Avicenna J Med Biotechnol; 6(2):94–103. 24834311. Accessed 2015-06-07.
Pšenička M, Saito T, Linhartová Z, Gazo I. Isolation and transplantation of sturgeon early-stage germ cells. Theriogenology; 83(6):1085–92. [doi:10.1016/j.theriogenology.2014.12.010, 25559841. Accessed 2015-06-10].
Yoshizaki G, Fujinuma K, Iwasaki Y, Okutsu T, Shikina S, Yazawa R, et al. Spermatogonial transplantation in fish: A novel method for the preservation of genetic resources. Comp Biochem Physiol Part D Genomics Proteomics; 6(1):55–61. [doi:10.1016/j.cbd.2010.05.003].
Bar I, Smith A, Bubner E, Yoshizaki G, Takeuchi Y, Yazawa R, et al. Assessment of yellowtail kingfish (Seriola lalandi) as a surrogate host for the production of southern bluefin tuna (Thunnus maccoyii) seed via spermatogonial germ cell transplantation. Reprod Fertil Dev. [doi:10.1071/RD15136. Accessed 2015-07-21].
Sayers EW, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin, V, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res; 37:5–15. [doi:10.1093/nar/gkn741. 18940862].
Nakamura Y, Mori K, Saitoh K, Oshima K, Mekuchi M, Sugaya T, et al. Evolutionary changes of multiple visual pigment genes in the complete genome of Pacific bluefin tuna. PNAS; 110(27):11061–6. [doi:10.1073/pnas.1302051110. 23781100].
Zhang W, Chen J, Yang Y, Tang Y, Shang J, Shen B. A Practical Comparison of De Novo Genome Assembly Software Tools for Next-Generation Sequencing Technologies. PLoS ONE; 6(3):17915. [doi:10.1371/journal.pone.0017915. Accessed 2015-04-15].
Bao S, Jiang R, Kwan W, Wang B, Ma X, Song YQ. Evaluation of next-generation sequencing software in mapping and assembly. J Hum Genet; 56(6):406–414. [doi:10.1038/jhg.2011.43. Accessed 2015-04-15].
Lee HC, Lai K, Lorenc MT, Imelfort M, Duran C, Edwards D. Bioinformatics tools and databases for analysis of next-generation sequence data. Brief Funct Genomics; 11(1):12–24. [doi:10.1093/bfgp/elr037. 22184335. Accessed 2015-04-15].
Nekrutenko A, Taylor J. Next-generation sequencing data interpretation: enhancing reproducibility and accessibility. Nat Rev Genet; 13(9):667–672. [doi:10.1038/nrg3305. Accessed 2015-04-15].
Morozova O, Hirst M, Marra MA. Applications of New Sequencing Technologies for Transcriptome Analysis. Annu. Rev. Genomics Hum. Genet; 10(1):135–151. [doi:10.1146/annurev-genom-082908-145957. 19715439. Accessed 2015-03-09].
NHMRC. Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 8th Edition: National Health and Medical Research Council. https://www.nhmrc.gov.au/_files_nhmrc/publications/attachments/ea28_code_care_use_animals_131209.pdf.
Yoshizaki G, Ichikawa M, Hayashi M, Iwasaki Y, Miwa M, Shikina S, et al. Sexual plasticity of ovarian germ cells in rainbow trout. Development; 137(8):1227–30. [doi:10.1242/dev.044982. 20223765. Accessed 2014-04-09].
Fischer AH, Jacobson KA, Rose J, Zeller R. Hematoxylin and eosin staining of tissue and cell sections. CSH Protoc. 2008;4986. 21356829.
Mackie M, Lewis P. Assessment of Gonad Staging Systems and Other Methods Used in the Study of the Reproductive Biology of Narrow-barred Spanish Mackerel, Scomberomorus Commerson, in Western Australia.
Bubner E, Farley J, Thomas P, Bolton T, Elizur A. Assessment of reproductive maturation of southern bluefin tuna (Thunnus maccoyii) in captivity. Aquaculture; 364-365:82–95. [doi:10.1016/j.aquaculture.2012.08.005].
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity. Nat Protoc; 8(8):1494–1512. [doi:10.1038/nprot.2013.084. 23845962].
Camacho C, Madden T, Ma N, Tao T, Agarwala R, Morgulis A. BLAST Command Line Applications User Manual. http://www.ncbi.nlm.nih.gov/books/NBK1763/. Accessed 2014-04-04.
Eddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol; 7(10):1002195. [doi:10.1371/journal.pcbi.1002195. Accessed 2015-04-11].
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy, SR, et al. Pfam: the protein families database. Nucleic Acids Res; 42:222–30. [doi:10.1093/nar/gkt1223.24288371. Accessed 2015-04-11].
GO Consortium. Gene ontology annotations and resources; 41:530–5. [doi:10.1093/nar/gks1050. Accessed 2015-04-21].
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Meth; 8(10):785–6. [doi:10.1038/nmeth.1701. Accessed 2015-04-21.]
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biolo; 305(3):567–80. [doi:10.1006/jmbi.2000.4315. Accessed 2015-04-21].
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics; 12(1):323. [doi:10.1186/1471-2105-12-323, 21816040. Accessed 2015-04-22].
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics; 26(1):139–40. [doi:10.1093/bioinformatics/btp616, 19910308. Accessed 2015-04-22].
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol; 11(2):14. [doi:10.1186/gb-2010-11-2-r14. Accessed 2015-10-21].
Ihaka R, Gentleman R. R: A Language for Data Analysis and Graphics. J Comput Graph Stat; 5(3):299–314. [doi:10.1080/10618600.1996.10474713. Accessed 2015-04-26.]
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit, S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol; 5(10):80. doi:10.1186/gb-2004-5-10-r80, 15461798. Accessed 2015-04-26.
Diaz de Cerio O, Rojo-Bartolomé I, Bizarro C, Ortiz-Zarragoitia M, Cancio I. 5S rRNA and Accompanying Proteins in Gonads: Powerful Markers to Identify Sex and Reproductive Endocrine Disruption in Fish. Environ Sci Technol; 46(14):7763–71. [doi:10.1021/es301132b. Accessed 2015-06-09].
Manousaki T, Tsakogiannis A, Lagnel J, Sarropoulou E, Xiang JZ, Papandroulakis N, et al. The sex-specific transcriptome of the hermaphrodite sparid sharpsnout seabream (Diplodus puntazzo). BMC Genomics; 15(1):655. [doi:10.1186/1471-2164-15-655, 25099474. Accessed 2015-06-09].
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotech; 29(7):644–52. [doi:10.1038/nbt.1883. Accessed 2014-02-25].
Li B, Fillmore N, Bai Y, Collins M, Thomson JA, Stewart R, et al. Evaluation of de novo transcriptome assemblies from RNA-Seq data. bioRxiv:006338. [doi:10.1101/006338. Accessed 2014-10-29].
Moreton J, Dunham SP, Emes RD. A consensus approach to vertebrate de novo transcriptome assembly from RNA-seq data: assembly of the duck (anas platyrhynchos) transcriptome. Front Genet. 5(190). doi:10.3389/fgene.2014.00190.
Lu B, Zeng Z, Shi T. Comparative study of de novo assembly and genome-guided assembly strategies for transcriptome reconstruction based on RNA-Seq. Sci China Life Sci; 56(2):143–55. [doi:10.1007/s11427-013-4442-z. Accessed 2014-03-22].
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14(4):36. doi:10.1186/gb-2013-14-4-r36. 23618408. Accessed 2014-03-23.
Bao E, Jiang T, Girke T. BRANCH: boosting RNA-Seq assemblies with partial or related genomic sequences10. Bioinformatics; 29:1250–9. [doi:10.1093/bioinformatics/btt127.23493323. Accessed 2014-04-09].
Bao E, Jiang T, Girke T. AlignGraph: algorithm for secondary de novo genome assembly guided by closely related references. Bioinformatics; 30(12):319–28. [doi:10.1093/bioinformatics/btu291, 24932000].
Takeyama H, Chow S, Tsuzuki H, Matunaga T. Mitochondrial DNA sequence variation within and between tuna Thunnus species and its application to species identification. J Fish Biol; 58(6):1646–57. [doi:10.1111/j.1095-8649.2001.tb02319.x. Accessed 2015-06-08].
Fitzgibbon QP, Seymour RS, Ellis D, Buchanan J. The energetic consequence of specific dynamic action in southern bluefin tuna Thunnus maccoyii. J Exp Biol; 210(2):290–8. [doi:10.1242/jeb.02641.17210965. Accessed 2015-06-08].
Shadwick RE, Syme DA. Thunniform swimming: muscle dynamics and mechanical power production of aerobic fibres in yellowfin tuna (Thunnus albacares). J Exp Biol; 211(10):1603–11. [doi:10.1242/jeb.013250.18456888. Accessed 2015-06-08].
Viñas J, Tudela S. A Validated Methodology for Genetic Identification of Tuna Species (Genus Thunnus). PLoS ONE; 4(10):7606. [doi:10.1371/journal.pone.0007606. Accessed 2015-06-08].
Lubzens E, Young G, Bobe J, Cerdà J. Oogenesis in teleosts: How fish eggs are formed. Gen Comp Endocrinol; 165(3):367–89. [doi:10.1016/j.ygcen.2009.05.022].
Böhne A, Sengstag T, Salzburger W. Comparative Transcriptomics in East African Cichlids Reveals Sex- and Species-Specific Expression and New Candidates for Sex Differentiation in Fishes. Genome Biol Evol; 6(9):2567–85. [doi:10.1093/gbe/evu200. Accessed 2015-05-18].
Sreenivasan R, Cai M, Bartfai R, Wang X, Christoffels A, Orban L. Transcriptomic Analyses Reveal Novel Gene s with Sexually Dimorphic Expression in the Zebrafish Gonad and Brain. PLoS ONE; 3(3):1791. [doi:10.1371/journal.pone.0001791, 18335061. Accessed 2015-08-23].
Small CM, Carney GE, Mo Q, Vannucci M, Jones AG. A microarray analysis of sex- and gonad-biased gene expression in the zebrafish: Evidence for masculinization of the transcriptome. BMC Genomics; 10(1):1. [doi:10.1186/1471-2164-10-579. Accessed 2015-08-23].
Kobayashi T, Nakamura M, Kajiura-Kobayashi H, Young G, Nagahama Y. Immunolocalization of steroidogenic enzymes (P450scc, P450c17, P450arom, and 3 β-HSD) in immature and mature testes of rainbow trout (Oncorhynchus mykiss). Cell Tissue Res; 292(3):573–7. Accessed 2015-05-19.
Mindnich R, Adamski J. Zebrafish 17beta-hydroxysteroid dehydrogenases: An evolutionary perspective. Mol Cell Endocrinol; 301(1-2):20–6. [doi:10.1016/j.mce.2008.12.002. Accessed 2015-05-19].
Milla S, Wang N, Mandiki SNM, Kestemont P. Corticosteroids: Friends or foes of teleost fish reproduction?. Comp Biochem Physiol A Mol Integr Physiol; 153(3):242–251. [doi:10.1016/j.cbpa.2009.02.027. Accessed 2015-05-17].
Ozaki Y, Higuchi M, Miura C, Yamaguchi S, Tozawa Y, Miura T. Roles of 11 β-Hydroxysteroid Dehydrogenase in Fish Spermatogenesis. Endocrinology; 147(11):5139–46. [doi:10.1210/en.2006-0391. Accessed 2015-05-19].
Hinfray N, Nóbrega RH, Caulier M, Baudiffier D, Maillot-Maréchal E, Chadili E, et al. Cyp17a1 and Cyp19a1 in the zebrafish testis are differentially affected by oestradiol. J Endocrinol; 216(3):375–88. [doi:10.1530/JOE-12-0509. 23283575. Accessed 2015-05-21].
Tao W, Yuan J, Zhou L, Sun L, Sun Y, Yang S, et al. Characterization of Gonadal Transcriptomes from Nile Tilapia (Oreochromis niloticus) Reveals Differentially Expressed Genes. PLoS ONE; 8(5):63604. [doi:10.1371/journal.pone.0063604. Accessed 2015-03-09].
Kobayashi T, Nagahama Y. Molecular Aspects of Gonadal Differentiation in a Teleost Fish, the Nile Tilapia; 3(2-3):108–17. [doi:10.1159/000223076].
Yang Y-H, Miranda CL, Henderson MC, Wang-Buhler JL, Buhler DR. Heterologous Expression of Cyp2k1 and Identification of the Expressed Protein (bv-Cyp2k1) As Lauric Acid (ω-1)-Hydroxylase and Aflatoxin B1exo-Epoxidase. Drug Metab Dispos; 28(11):1279–83. 11038153. Accessed 2015-06-07.
James MO. Steroid catabolism in marine and freshwater fish. J Steroid Biochem Mol Biol; 127(3-5):167–75. [doi:10.1016/j.jsbmb.2010.10.003. Accessed 2015-06-07].
Zhu S, King SC, Haasch ML. Biomarker induction in tropical fish species on the Northwest Shelf of Australia by produced formation water. Mar Environ Res; 65(4):315–24. [doi:10.1016/j.marenvres.2007.11.007. Accessed 2015-05-21].
Patterson TA, Evans K, Carter TI, Gunn JS. Movement and behaviour of large southern bluefin tuna (Thunnus maccoyii) in the Australian region determined using pop-up satellite archival tags. Fish Oceanogr; 17(5):352–67. [doi:10.1111/j.1365-2419.2008.00483.x. Accessed 2015-09-08].
Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. Nature; 417(6888):559–63. Accessed 2015-05-20.
Bratuś A, Słota E. DMRT1/Dmrt1, the Sex Determining or Sex Differentiating Gene in Vertebrata. Folia Biologica; 54(3-1):81–6. [doi:10.3409/173491606778557563].
Herpin A, Schartl M. Dmrt1 genes at the crossroads: a widespread and central class of sexual development factors in fish. FEBS J; 278(7):1010–9. [doi:10.1111/j.1742-4658.2011.08030.x. Accessed 2015-05-20].
Cavileer T, Hunter S, Okutsu T, Yoshizaki G, Nagler JJ. Identification of Novel Genes Associated with Molecular Sex Differentiation in the Embryonic Gonads of Rainbow Trout (Oncorhynchus mykiss); 3(4):214–24. [doi:10.1159/000228722. Accessed 2015-12-18].
Graf M, Teo Qi-Wen E-R, Sarusie MV, Rajaei F, Winkler C. Dmrt5 controls corticotrope and gonadotrope differentiation in the zebrafish pituitary. Mol Endocrinol; 29(2):187–99. [doi:10.1210/me.2014-1176, 25489906].
Guo Y, Li Q, Gao S, Zhou X, He Y, Shang X, et al. Molecular cloning, characterization, and expression in brain and gonad of Dmrt5 of zebrafish. Biochem Biophys Res Commun; 324(2):569–75. [doi:10.1016/j.bbrc.2004.09.085. Accessed 2015-08-24].
Raverot G, Weiss J, Park SY, Hurley L, Jameson JL. Sox3 expression in undifferentiated spermatogonia is required for the progression of spermatogenesis. Developmental Biology; 283(1):215–25. [doi:10.1016/j.ydbio.2005.04.013. Accessed 2015-05-31].
Sarkar A, Hochedlinger K. The Sox Family of Transcription Factors: Versatile Regulators of Stem and Progenitor Cell Fate. Cell Stem Cell; 12(1):15–30. [doi:10.1016/j.stem.2012.12.007. Accessed 2015-05-31].
Xia W, Zhou L, Yao B, Li C-J, Gui J-F. Differential and spermatogenic cell-specific expression of DMRT1 during sex reversal in protogynous hermaphroditic groupers. Mol Cell Endocrinol; 263(1-2):156–72. [doi:10.1016/j.mce.2006.09.014. Accessed 2015-05-31].
Shin HS, An KW, Park MS, Jeong MH, Choi CY. Quantitative mRNA expression of sox3 and DMRT1 during sex reversal, and expression profiles after GnRHa administration in black porgy, Acanthopagrus schlegeli. Comp Biochem Physiol B Biochem Mol Biol; 154(1):150–6. [doi:10.1016/j.cbpb.2009.05.013. Accessed 2015-05-31].
Han F, Wang Z, Wu F, Liu Z, Huang B, Wang D. Characterization, phylogeny, alternative splicing and expression of Sox30 gene. BMC Mol Biol; 11(1):98. [doi:10.1186/1471-2199-11-98, 21143990. Accessed 2015-05-31].
Han F, Dong Y, Liu W, Ma X, Shi R, Chen H, et al. Epigenetic Regulation of Sox30 Is Associated with Testis Development in Mice. PLoS ONE; 9(5):97203. [doi:10.1371/journal.pone.0097203. Accessed 2015-05-02].
Brennan J, Capel B. One tissue, two fates: molecular genetic events that underlie testis versus ovary development. Nat Rev Genet; 5(7):509–21. [doi:10.1038/nrg1381. Accessed 2015-05-30].
Ottolenghi C, Pelosi E, Tran J, Colombino M, Douglass E, Nedorezov T, et al. Loss of Wnt4 and Foxl2 leads to female-to-male sex reversal extending to germ cells. Hum Mol Genet; 16(23):2795–804. [doi:10.1093/hmg/ddm235, 17728319. Accessed 2015-05-31].
Wu GC, Chang C-F. wnt4 Is Associated with the Development of Ovarian Tissue in the Protandrous Black Porgy, Acanthopagrus schlegeli. Biol Reprod; 81(6):1073–82. [doi:10.1095/biolreprod.109.077362, 19605791. Accessed 2015-05-31].
Chen H, Li S, Xiao L, Zhang Y, Li G, Liu X, et al. Wnt4 in protogynous hermaphroditic orange-spotted grouper (Epinephelus coioides): Identification and expression. Comp Biochem Physiol B Biochem Mol Biol; 183:67–74. [doi:10.1016/j.cbpb.2015.01.002. Accessed 2015-05-31].
Caburet S, Georges A, L’Hôte D, Todeschini AL, Benayoun BA, Veitia RA. The transcription factor FOXL2: At the crossroads of ovarian physiology and pathology. Mol Cell Endocrinol; 356(1-2):55–64. [doi:10.1016/j.mce.2011.06.019. Accessed 2015-06-09].
shou Wang D, Kobayashi T, yan Zhou L, Nagahama Y. Molecular cloning and gene expression of Foxl2 in the Nile tilapia, Oreochromis niloticus. Biochem Biophys Res Commun; 320(1):83–9. [doi:10.1016/j.bbrc.2004.05.133. Accessed 2015-06-09].
Nakamoto M, Matsuda M, Wang DS, Nagahama Y, Shibata N. Molecular cloning and analysis of gonadal expression of Foxl2 in the medaka, Oryzias latipes. Biochem Biophys Res Commun; 344(1):353–61. [doi:10.1016/j.bbrc.2006.03.137. Accessed 2015-06-09].
Nakamoto M, Muramatsu S, Yoshida S, Matsuda M, Nagahama Y, Shibata N. Gonadal sex differentiation and expression of Sox9a2, Dmrt1, and Foxl2 in Oryzias luzonensis. Genesis; 47(5):289–99. [doi:10.1002/dvg.20498. Accessed 2015-06-09].
Yamaguchi T, Yamaguchi S, Hirai T, Kitano T. Follicle-stimulating hormone signaling and Foxl2 are involved in transcriptional regulation of aromatase gene during gonadal sex differentiation in Japanese flounder, Paralichthys olivaceus. Biochem Biophys Res Commun; 359(4):935–40. [doi:10.1016/j.bbrc.2007.05.208. Accessed 2015-06-09].
Jiang W, Yang Y, Zhao D, Liu X, Duan J, Xie S, Zhao H. Effects of sexual steroids on the expression of foxl2 in Gobiocypris rarus. Comp Biochem Physiol B Biochem Mol Biol; 160(4):187–93. [doi:10.1016/j.cbpb.2011.08.005. Accessed 2015-06-09].
Kobayashi Y, Horiguchi R, Nozu R, Nakamura M. Expression and localization of forkhead transcriptional factor 2 (Foxl2) in the gonads of protogynous wrasse, Halichoeres trimaculatus. Biol Sex Differ; 1(1):3. Accessed 2015-06-09.
Hu Q, Guo W, Gao Y, Tang R, Li D. Molecular cloning and analysis of gonadal expression of Foxl2 in the rice-field eel Monopterus albus. Sci. Rep; 4:6884. [doi:10.1038/srep06884. Accessed 2015-06-09].
Filby AL, van Aerle R, Duitman J, Tyler CR. The Kisspeptin/Gonadotropin-Releasing Hormone Pathway and Molecular Signaling of Puberty in Fish. Biol Reprod; 78(2):278–89. [doi:10.1095/biolreprod.107.063420, 17978278. Accessed 2015-04-26].
Zohar Y, Muñoz-Cueto JA, Elizur A, Kah O. Neuroendocrinology of reproduction in teleost fish. Gen Comp Endocrinol; 165(3):438–55. [doi:10.1016/j.ygcen.2009.04.017].
Selvaraj S, Ohga H, Kitano H, Nyuji M, Yamaguchi A, Matsuyama M. Peripheral Administration of Kiss1 Pentadecapeptide Induces Gonadal Development in Sexually Immature Adult Scombroid Fish. Zool Sci; 30(6):446–54. [doi:10.2108/zsj.30.446. Accessed 2015-05-21].
Ohga H, Fujinaga Y, Selvaraj S, Kitano H, Nyuji M, Yamaguchi A, et al. Identification, characterization, and expression profiles of two subtypes of kisspeptin receptors in a scombroid fish (chub mackerel). Gen Comp Endocrinol; 193:130–40. [doi:10.1016/j.ygcen.2013.07.016. Accessed 2015-05-21].
Selvaraj S, Kitano H, Ohga H, Yamaguchi A, Matsuyama M. Expression changes of mRNAs encoding kisspeptins and their receptors and gonadotropin-releasing hormones during early development and gonadal sex differentiation periods in the brain of chub mackerel (scomber japonicus). 20–32; 222. [doi:10.1016/j.ygcen.2014.09.019].
Betsholtz C, Karlsson L, Lindahl P. Developmental roles of platelet-derived growth factors. BioEssays; 23(6):494–507. [doi:10.1002/bies.1069. Accessed 2015-05-30].
Basciani S, Mariani S, Arizzi M, Ulisse S, Rucci N, Jannini EA, et al. Expression of Platelet-Derived Growth Factor-A (PDGF-A):PDGF-B, and PDGF Receptor- α and - β during Human Testicular Development and Disease. J Clin Endocrinol Metab; 87(5):2310–9. [doi:10.1210/jcem.87.5.8476. Accessed 2015-05-30].
Vanorny DA, Prasasya RD, Chalpe AJ, Kilen SM, Mayo KE. Notch Signaling Regulates Ovarian Follicle Formation and Coordinates Follicular Growth. Mol Endocrinol; 28(4):499–511. [doi:10.1210/me.2013-1288. Accessed 2015-05-22].
Prasasya RD, Mayo KE. Notch Pathway Genes and Activity Are Regulated By Gonadotropins in the Mouse Ovary. In: Female Reproduction: Regulation of Follicle Growth. Meeting Abstracts. p. 0051–0051. Endocrine Society. http://press.endocrine.org/doi/abs/10.1210/endo-meetings.2014.RE.5.MON-0051. Accessed 2015-05-21.
Lacerda SMdSN, Costa GMJ, de França LR. Biology and identity of fish spermatogonial stem cell. Gen Comp Endocrinol; 207:56–65. [doi:10.1016/j.ygcen.2014.06.018. Accessed 2014-12-01].
Sánchez-Sánchez AV, Camp E, García-España A, Leal-Tassias A, Mullor JL. Medaka Oct4 is expressed during early embryo development, and in primordial germ cells and adult gonads; 239(2):672–9. [doi:10.1002/dvdy.22198].
Wang D, Manali D, Wang T, Bhat N, Hong N, Li Z, et al. Identification of Pluripotency Genes in the Fish Medaka. Int J Biol Sci; 7(4):440–51. 21547061. Accessed 2015-05-20.
Fabra M, Raldú D, Bozzo MG, Deen PMT, Lubzens E, Cerdà J. Yolk proteolysis and aquaporin-1o play essential roles to regulate fish oocyte hydration during meiosis resumption. Devel Biol; 295(1):250–62. [doi:10.1016/j.ydbio.2006.03.034. Accessed 2015-06-01].
Rolland AD, Lareyre J-J, Goupil AS, Montfort J, Ricordel MJ, Esquerré D, et al. Expression profiling of rainbow trout testis development identifies evolutionary conserved genes involved in spermatogenesis. BMC Genomics; 10(1):546. [doi:10.1186/1471-2164-10-546, 19925684. Accessed 2015-06-08].
Miura C, Shimizu Y, Uehara M, Ozaki Y, Young G, Miura T. Gh is produced by the testis of Japanese eel and stimulates proliferation of spermatogonia. Reproduction; 142(6):869–77. [doi:10.1530/REP-11-0203. Accessed 2015-05-23].
Davies R, Moore A, Schedl A, Bratt E, Miyahawa K, Ladomery M, et al. Multiple Roles for the Wilms’ Tumor Suppressor, WT1. Cancer Res; 59(7 Supplement):1747–51. 10197591. Accessed 2015-05-23.
Gashler AL, Bonthron DT, Madden SL, Rauscher FJ, Collins T, Sukhatme VP. Human platelet-derived growth factor A chain is transcriptionally repressed by the Wilms tumor suppressor WT1. PNAS; 89(22):10984–8. 1332065. Accessed 2015-05-31.
Doitsidou M, Reichman-Fried M, Stebler J, Köprunner M, Dörries J, Meyer D, et al. Guidance of Primordial Germ Cell Migration by the Chemokine SDF-1. Cell; 111(5):647–59. [doi:10.1016/s0092-8674(02)01135-2].
Knaut H, Werz C, Geisler R, The Tubingen Screen C, Nusslein-Volhard C. A zebrafish homologue of the chemokine receptor Cxcr4 is a germ-cell guidance receptor. Nature; 421(6920):279–82.
Saito T, Goto-Kazeto R, Kawakami Y, Nomura K, Tanaka H, Adachi, S, et al. The Mechanism for Primordial Germ-Cell Migration Is Conserved between Japanese Eel and Zebrafish. PLoS ONE; 6(9):24460. [doi:10.1371/journal.pone.0024460].
Pradeepa MM, Manjunatha S, Sathish V, Agrawal S, Rao MRS. Involvement of Importin-4 in the Transport of Transition Protein 2 into the Spermatid Nucleus. Mol Cell Biol; 28(13):4331–41. [doi:10.1128/MCB.00519-07, 17682055. Accessed 2015-06-01].
Butler R, Jans DA, Loveland KL, Major A, Miyamoto Y, Young JC. Putting things in place for fertilization: Discovering roles for importin proteins in cell fate and spermatogenesis. Asian J Androl; 17(4):537–44. [doi:10.4103/1008-682X.154310].
Houwing S, Kamminga LM, Berezikov E, Cronembold D, Girard A, van den Elst H, et al. A Role for Piwi and piRNAs in Germ Cell Maintenance and Transposon Silencing in Zebrafish. Cell; 129(1):69–82. [doi:10.1016/j.cell.2007.03.026].
O’Donnell KA, Boeke JD. Mighty Piwis Defend the Germline against Genome Intruders. Cell; 129(1):37–44.
Juliano C, Wang J, Lin H. Uniting Germline and Stem Cells: the Function of Piwi Proteins and the piRNA Pathway in Diverse Organisms. Annu Rev Genet; 45:447–469. [doi:10.1146/annurev-genet-110410-132541, 21942366].
Miura C, Miura T, Yamashita M. PCNA Protein Expression during Spermatogenesis of the Japanese Eel (Anguilla japonica). Zool Sci; 19(1):87–91. [doi:10.2108/zsj.19.87. Accessed 2015-05-23].
Higuchi K, Takeuchi Y, Miwa M, Yamamoto Y, Tsunemoto K, Yoshizaki G. Colonization, proliferation, and survival of intraperitoneally transplanted yellowtail Seriola quinqueradiata spermatogonia in nibe croaker Nibea mitsukurii recipient. Fish Sci; 77(1):69–77. [doi:10.1007/s12562-010-0314-7].
Sawatari E, Shikina S, Takeuchi T, Yoshizaki G. A novel transforming growth factor-beta superfamily member expressed in gonadal somatic cells enhances primordial germ cell and spermatogonial proliferation in rainbow trout (Oncorhynchus mykiss). Dev Biol; 301(1):266–75. [doi:10.1016/j.ydbio.2006.10.001, 17109839].
Shibata Y, Paul-Prasanth B, Suzuki A, Usami T, Nakamoto M, Matsuda M, et al. Expression of gonadal soma derived factor (GSDF) is spatially and temporally correlated with early testicular differentiation in medaka. Gene Expr Patterns; 10(6):283–9. [doi:10.1016/j.gep.2010.06.005].
Lacerda S, Aponte PM, Costa GMJ, Campos-Junior PHA, Segatelli TM, Silva MA, et al. An overview of spermatogonial stem cell physiology, niche and transplantation in fish. Anim Reprod; 9(4):798–808. Accessed 2015-05-23.
Loir M. In vitro approach to the control of spermatogonia proliferation in the trout. Mol Cell Endocrinol; 102(1-2):141–150. [doi:10.1016/0303-7207(94)90107-4. Accessed 2015-05-31].
Nader MR, Miura T, Ando N, Miura C, Yamauchi K. Recombinant Human Insulin-Like Growth Factor I Stimulates All Stages of 11-Ketotestosterone-Induced Spermatogenesis in the Japanese Eel, Anguilla japonica, In Vitro. Biol Reprod; 61(4):944–7. [doi:10.1095/biolreprod61.4.944, 10491628. Accessed 2015-05-31].
Perrot V, Moiseeva EB, Gozes Y, Chan SJ, Funkenstein B. Insulin-Like Growth Factor Receptors and Their Ligands in Gonads of a Hermaphroditic Species, the Gilthead Seabream (Sparus aurata): Expression and Cellular Localization. Biol Reprod; 63(1):229–41. [doi:10.1095/biolreprod61.4.944, 10491628. Accessed 2015-05-31].
Nakamura S, Watakabe I, Nishimura T, Toyoda A, Taniguchi Y, Tanaka M. Analysis of Medaka sox9 Orthologue Reveals a Conserved Role in Germ Cell Maintenance. PLoS ONE; 7(1):29982. [doi:10.1371/journal.pone.0029982].
Yazawa R, Takeuchi Y, Higuchi K, Yatabe T, Kabeya N, Yoshizaki G. Chub Mackerel Gonads Support Colonization, Survival, and Proliferation of Intraperitoneally Transplanted Xenogenic Germ Cells. Biol Reprod; 82(5):896–904. doi:10.1095/biolreprod.109.081281.
Saito T, Goto-Kazeto R, Fujimoto T, Kawakami Y, Arai K, Yamaha E. Inter-species transplantation and migration of primordial germ cells in cyprinid fish. Int J Dev Biol; 54(10):1481.
Herrid M, Davey RJ, Hutton K, Colditz IG, Hill JR. A comparison of methods for preparing enriched populations of bovine spermatogonia. Reprod Fertil Dev; 21(3):393–9. [doi:10.1071/RD08129].
Yazawa R, Takeuchi Y, Morita T, Ishida M, Yoshizaki G. The Pacific bluefin tuna (Thunnus orientalis) dead end gene is suitable as a specific molecular marker of type A spermatogonia. Mol Reprod Dev; 80(10):871–80. [doi:10.1002/mrd.22224. Accessed 2013-10-28].
Viñas J, Piferrer F. Stage-Specific Gene Expression During Fish Spermatogenesis as Determined by Laser-Capture Microdissection and Quantitative-PCR in Sea Bass (Dicentrarchus labrax) Gonads. Biol Reprod; 79(4):738–47. [doi:10.1095/biolreprod.108.069708].
Kobayashi T, Yoshizaki G, Takeuchi Y, Takeuchi T. Isolation of highly pure and viable primordial germ cells from rainbow trout by GFP-dependent flow cytometry. Mol Reprod Dev; 67(1):91–100. [doi:10.1002/mrd.20003].
Yano A, Suzuki K, Yoshizaki G. Flow-Cytometric Isolation of Testicular Germ Cells from Rainbow Trout(Oncorhynchus mykiss) Carrying the Green Fluorescent Protein Gene Driven by Trout vasa Regulatory Regions. Biol Reprod; 78(1):151–8. [doi:10.1095/biolreprod.107.064667].
Yoshizaki G, Okutsu T, Ichikawa M, Hayashi M, Takeuchi Y. Sexual plasticity of rainbow trout germ cells. Anim Reprod; 7(3):187–96.
Kise K, Yoshikawa H, Sato M, Tashiro M, Yazawa R, Nagasaka Y, et al. Flow-Cytometric Isolation and Enrichment of Teleost Type A Spermatogonia Based on Light-Scattering Properties. Biol Reprod; 86(4):107. [doi:10.1095/biolreprod.111.093161, 22219211. Accessed 2015-06-07].
Tanaka M, Kinoshita M, Kobayashi D, Nagahama Y. Establishment of Medaka (Oryzias latipes) Transgenic Lines with the Expression of Green Fluorescent Protein Fluorescence Exclusively in Germ Cells: A Useful Model to Monitor Germ Cells in a Live Vertebrate. Proc Natl Acad Sci USA; 98(5):2544–9. 3055092.
Yoshizaki G, Takeuchi Y, Sakatani S, Takeuchi T. Germ cell-specific expression of green fluorescent protein in transgenic rainbow trout under control of the rainbow trout vasa-like gene promoter. Int J Dev Biol; 44(3):323–6.
Kobayashi T, Kajiura-Kobayashi H, Nagahama Y. Differential expression of vasa homologue gene in the germ cells during oogenesis and spermatogenesis in a teleost fish, tilapia, Oreochromis niloticus. Mech Dev; 99(1-2):139–42. [doi:10.1016/s0925-4773(00)00464-0].
Åkugor A, Tveiten H, Krasnov A, Andersen Ø. Knockdown of the germ cell factor Dead end induces multiple transcriptional changes in Atlantic cod (Gadus morhua) hatchlings. Anim Reprod Sci; 144(3-4):129–37. [doi:10.1016/j.anireprosci.2013.12.010. Accessed 2014-02-24].
Pacchiarini T, Sarasquete C, Cabrita E. Development of interspecies testicular germ-cell transplantation in flatfish. Reprod Fertil Dev; 26(5):690–702. Accessed 2014-11-26.
Yoon C, Kawakami K, Hopkins N. Zebrafish vasa homologue RNA is localized to the cleavage planes of 2- and 4-cell-stage embryos and is expressed in the primordial germ cells. Development; 124(16):3157–65.
Nagasawa K, Takeuchi Y, Miwa M, Higuchi K, Morita T, Mitsuboshi, T, et al. cDNA cloning and expression analysis of a Vasa-like gene in Pacific bluefin tuna Thunnus orientalis. Fish Sci; 75(1):71–9. [doi:10.1007/s12562-008-0021-9].
Olsen LC, Aasland R, Fjose A. A vasa-like gene in zebrafish identifies putative primordial germ cells. Mech Dev; 66(1-2):95–105. [doi:10.1016/s0925-4773(97)00099-3].
Dumstrei K, Mennecke R, Raz E. Signaling pathways controlling primordial germ cell migration in zebrafish. J Cell Sci; 117(20):4787–95. [doi:10.1242/jcs.01362].
Kurokawa H, Saito D, Nakamura S, Katoh-Fukui Y, Ohta K, Baba T, et al. Germ cells are essential for sexual dimorphism in the medaka gonad. Proc Natl Acad Sci; 104(43):16958–63. [doi:10.1073/pnas.0609932104].
Costa GMJ, Avelar GF, Rezende-Neto JV, Campos-Junior PHA, Lacerda SMSN, Andrade BSC, et al. Spermatogonial Stem Cell Markers and Niche in Equids. PLoS ONE; 7(8):44091. [doi:10.1371/journal.pone.0044091].
Santos Nassif Lacerda SM, Costa GMJ, da Silva MdA, Almeida Campos-Junior PH, Segatelli TM, Peixoto MTD, et al. Phenotypic characterization and in vitro propagation and transplantation of the Nile tilapia (Oreochromis niloticus) spermatogonial stem cells. Gen Comp Endocrinol; 192:95–106. [doi:10.1016/j.ygcen.2013.06.013, Accessed 2015-04-26].
Nakajima S, Hayashi M, Kouguchi T, Yamaguchi K, Miwa M, Yoshizaki G. Expression patterns of gdnf and gfr α1 in rainbow trout testis. Gene Expr Patterns; 14(2):111–20. [doi:10.1016/j.gep.2014.01.006 Accessed 2014-11-28].
Bellaïche J, Goupil AS, Sambroni E, Lareyre J-J, Gac FL. Gdnf-Gfra1 Pathway Is Expressed in a Spermatogenetic-Dependent Manner and Is Regulated by Fsh in a Fish Testis. Biol Reprod; 91(4):94. [doi:10.1095/biolreprod.114.119834 25165121].
Nagasawa K, Shikina S, Takeuchi Y, Yoshizaki G. Lymphocyte Antigen 75 (Ly75/CD205) Is a Surface Marker on Mitotic Germ Cells in Rainbow Trout. Biol Reprod; 83(4):597–606. [doi:10.1095/biolreprod.109.082081 20554922. Accessed 2015-04-26].
Köprunner M, Thisse C, Thisse B, Raz E. A zebrafish nanos-related gene is essential for the development of primordialgerm cells. Genes Dev; 15(21):2877–85. [doi:10.1101/gad.212401].
Hashimoto Y, Maegawa S, Nagai T, Yamaha E, Suzuki H, Yasuda K, Inoue K. Localized maternal factors are required for zebrafish germ cell formation. Dev Biol; 268(1):152–61. [doi:10.1016/j.ydbio.2003.12.013].
Saito T, Fujimoto T, Maegawa S, Inoue K, Tanaka M, Arai K, et al. Visualization of primordial germ cells in vivo using GFP-nos1 3’UTR mRNA. Int J Dev Biol; 50(8):691–9.
Bellaiche J, Lareyre J-J, Cauty C, Yano A, Allemand I, Gac FL. Spermatogonial Stem Cell Quest: nanos2, Marker of a Subpopulation of Undifferentiated A Spermatogonia in Trout Testis. Biol Reprod; 90(4):79. [doi:10.1095/biolreprod.113.116392 24554733].
Sada A, Hasegawa K, Pin PH, Saga Y. NANOS2 Acts Downstream of Glial Cell Line-Derived Neurotrophic Factor Signaling to Suppress Differentiation of Spermatogonial Stem Cells. Stem Cells; 30(2):280–91. [doi:10.1002/stem.790].
Yano A, von Schalburg K, Cooper G, Koop BF, Yoshizaki G. Identification of a molecular marker for type A spermatogonia by microarrayanalysis using gonadal cells from pvasa-GFP transgenic rainbow trout (oncorhynchus mykiss); 76(3):246–54. [doi:10.1002/mrd.20947].
Miura C, Miura T. Analysis of Spermatogenesis Using an Eel Model. Aqua-Biosci Monogr; 4(4):105–29. [doi:10.5047/absm.2011.00404.0105 Accessed 2015-05-27].
Froschauer A, Khatun MM, Sprott D, Franz A, Rieger C, Pfennig F, et al. oct4-EGFP reporter gene expression marks the stem cells in embryonic development and in adult gonads of transgenic medaka. Mol Reprod Dev; 80(1):48–58. [doi:10.1002/mrd.22135. Accessed 2013-03-30 ].
Onichtchouk D. Pou5f1/oct4 in pluripotency control: Insights from zebrafish; 50(2):75–85. [doi:10.1002/dvg.20800].
This project was supported by the Australian Seafood Cooperative Research Centre and Fisheries Research and Development Cooperation Innovation grant to AE (2011/730), a University of the Sunshine Coast Ph.D fellowship and HDR travel grant to IB, an Australian Academy of Sciences travel grant to Japan to AE and an ARC Future Fellowship to SC (FT110100990). We are grateful to Cleanseas Tuna LTD. for generously providing us with freshly harvested SBT and assisting with sampling and physical measurements. Final and special thanks to the Lincoln Marine Science Centre for the use of their lab facilities.
The authors declare that they have no competing interests.
IB conducted sample collection and preparation for histology and RNA-Seq, processed and analyzed the sequencing data, participated in the design of the study and drafted the manuscript. SC participated in the design of the study and helped to draft the manuscript. AE conceived of the study, participated in its design, assisted in data analysis and helped to draft the manuscript. All authors have read and approved the final manuscript.
Gene ontology (GO) terms enrichment in Southern bluefin tuna ovary vs. testis cells transcriptome. Enriched GO terms in ovary vs. testis cells, along with the corresponding ontology, term description, number of genes associated (total count and differentially expressed) and the p-value and p-value (false discovery rate adjusted p-value) statistics as determined by goseq. (XLSX 50.5 kb)
Differentially expressed gonad development and germ cell-related genes in ovary vs. testis cells. ORF – Open reading frame. Blast hit species, e-value, Bit Score – Blast best hit details. Log2FC – Fold change values (log2 transformed) as determined by edgeR, positive values represent over-expression in testis cells, negative values are over-expressed in the ovary cells. p-value, FDR (false discovery rate adjusted p-value) – edgeR differential expression analysis statistics. KEGG – Kyoto encyclopedia of genes and genomes function. (XLSX 35.7 kb)
Gene ontology (GO) terms enrichment in Southern bluefin tuna crude testis vs. Percoll-enriched testis cells transcriptome. Enriched GO terms in crude testis vs. Percoll-enriched testis cells, along with the corresponding ontology, term description, number of genes associated (total count and differentially expressed) and the p-value and p-value (false discovery rate adjusted p-value) statistics as determined by goseq. (XLSX 71.6 kb)
Differentially expressed gonad development and germ cell-related genes in crude vs. Percoll-enriched testis cells. ORF – Open reading frame. Blast hit species, e-value, Bit Score – Blast best hit details. Log2FC – Fold change values (log2 transformed) as determined by edgeR, positive values represent over-expression in testis cells, negative values are over-expressed in the ovary cells. p-value, FDR (false discovery rate adjusted p-value) – edgeR differential expression analysis statistics. KEGG – Kyoto encyclopedia of genes and genomes function. (XLSX 26.5 kb)
About this article
Cite this article
Bar, I., Cummins, S. & Elizur, A. Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii). BMC Genomics 17, 217 (2016). https://doi.org/10.1186/s12864-016-2397-8