Genomics of sablefish (Anoplopoma fimbria): expressed genes, mitochondrial phylogeny, linkage map and identification of a putative sex gene

Gonadal sex differentiation and effects of dietary methyltestosterone treatment in sablefish (Anoplopoma fimbria)

Methods for sex control are needed to establish monosex aquaculture of sablefish (Anoplopoma fimbria). Here we conducted the first characterization of sex differentiation by histology and hormonal sex reversal experiment in sablefish. Ovarian differentiation was first discernible at ~80 mm fork length (FL) and characterized by development of lamellar structures and onset of meiosis. Testes exhibited a dual-lobe appearance over much of their length and remained non-meiotic until males were ≥520 mm FL (2 years post-fertilization). Juveniles with undifferentiated gonads were provided diets containing 0 (control), 5 or 50 mg 17α-methyltestosterone (MT)/kg for 2 months. Following treatment, controls possessed either ovaries or non-meiotic testes, whereas MT-treated fish exhibited meiotic testes (60% of the fish), intersex gonads (~30%), or gonads that appeared sterile (~10%). A genetic sex marker revealed that all intersex fish were genetic females, although other females appeared to be completely sex reversed (i.e., neomales). One year after treatment, MT-treated fish possessed non-meiotic testes similar to control males or intersex gonads with reduced ovarian features, presumably due to atresia following MT withdrawal. Milt collected from neomales and genetic males 3 years post-treatment permitted sperm motility analyses; however, neomale sperm were virtually immotile. These results demonstrated that sablefish are differentiated gonochorists and that MT treatment from 76 to 196 mm FL induced permanent masculinization of a portion of the genetic females, but acquisition of sperm motility was impaired. Earlier administration of MT may be necessary to sex reverse a higher proportion of genetic females and reduce negative effects on fertility.

Molecular characterization of the gonadal kisspeptin system: Cloning, tissue distribution, gene expression analysis and localization in sablefish (Anoplopoma fimbria)

The kisspeptin system plays pivotal roles in the regulation of vertebrate reproduction. Classically, kisspeptin produced in the brain stimulates brain gonadotropin-releasing hormone signaling, which in turn activates the pituitary-gonad axis. Expression of the kisspeptin system has also been documented in peripheral tissues, including gonads of mammals and fishes. However, the fish gonadal kisspeptin system remained uncharacterized. Herein we report identification and characterization of four kisspeptin system mRNAs (kisspeptin 1 (kiss1), kiss2, and G protein-coupled receptor 54-1 (gpr54-1) and gpr54-2) in sablefish, Anoplopoma fimbria. Sablefish predicted protein sequences were highly similar to those of other marine teleosts, but less so to freshwater teleosts. Tissue distribution analyses revealed that all four kisspeptin-system transcripts were expressed in both brain and gonad. However, kiss2 was the predominant transcript in the gonads and the only transcript detected in ovulated eggs. Ontogenetic analysis of kiss2 expression in juvenile sablefish gonads demonstrated that levels were low during sex differentiation but increased with fish size and gonadal development. Dramatic increases in kiss2 mRNA occurred during primary oocyte growth, while levels remained relatively low in testes. In situ hybridization revealed that kiss2 mRNA was localized to cytoplasm of perinucleolus stage oocytes, suggesting it could play a local role in oogenesis or could be synthesized and stored within oocytes as maternal mRNA. This represents the first study to focus on the gonadal kisspeptin system in fishes and provides important tools for further investigation of both the gonadal and brain kisspeptin systems in sablefish.

Gene expression analysis at the onset of sex differentiation in turbot (Scophthalmus maximus)

BACKGROUND

Controlling sex ratios is essential for the aquaculture industry, especially in those species with sex dimorphism for relevant productive traits, hence the importance of knowing how the sexual phenotype is established in fish. Turbot, a very important fish for the aquaculture industry in Europe, shows one of the largest sexual growth dimorphisms amongst marine cultured species, being all-female stocks a desirable goal for the industry. Although important knowledge has been achieved on the genetic basis of sex determination (SD) in this species, the master SD gene remains unknown and precise information on gene expression at the critical stage of sex differentiation is lacking. In the present work, we examined the expression profiles of 29 relevant genes related to sex differentiation, from 60 up to 135 days post fertilization (dpf), when gonads are differentiating. We also considered the influence of three temperature regimes on sex differentiation.

RESULTS

The first sex-related differences in molecular markers could be observed at 90 days post fertilization (dpf) and so we have called that time the onset of sex differentiation. Three genes were the first to show differential expression between males and females and also allowed us to sex turbot accurately at the onset of sex differentiation (90 dpf): cyp19a1a, amh and vasa. The expression of genes related to primordial germ cells (vasa, gsdf, tdrd1) started to increase between 75-90 dpf and vasa and tdrd1 later presented higher expression in females (90-105 dpf). Two genes placed on the SD region of turbot (sox2, fxr1) did not show any expression pattern suggestive of a sex determining function. We also detected changes in the expression levels of several genes (ctnnb1, cyp11a, dmrt2 or sox6) depending on culture temperature.

CONCLUSION

Our results enabled us to identify the first sex-associated genetic cues (cyp19a1a, vasa and amh) at the initial stages of gonad development in turbot (90 dpf) and to accurately sex turbot at this age, establishing the correspondence between gene expression profiles and histological sex. Furthermore, we profiled several genes involved in sex differentiation and found specific temperature effects on their expression.

Development of approaches to induce puberty in cultured female sablefish (Anoplopoma fimbria)

Efforts to establish sustainable and efficient aquaculture production of sablefish (Anoplopoma fimbria) have been constrained by delayed puberty in cultured females. This study integrates a series of experiments aimed at gaining an understanding of the reproductive physiology of puberty in female sablefish. We detected transcripts for the dopamine D2 receptor (drd2) in brain, pituitary and ovary of sablefish, and prepubertal females exhibited significantly elevated brain and pituitary drd2 expression relative to wild maturing females. Treatments with sustained-release cholesterol pellets containing testosterone (T) and the dopamine D2 receptor antagonist, metoclopramide (Met), stimulated expression of pituitary luteinizing hormone beta subunit (lhb) and follicle-stimulating hormone beta subunit (fshb), respectively, in prepubertal females, whereas a combination of T and gonadotropin-releasing hormone agonist (GnRHa) had a strong synergistic effect on lhb expression (2000-fold higher than control). Although T induced a significant increase in the maximum ovarian follicle volume, none of the treatments tested stimulated onset of vitellogenesis. Using liquid chromatography/tandem mass spectrometry, we demonstrated that Met stimulated production of T by previtellogenic ovarian follicles in vitro, whereas gonadotropin preparations enhanced 17α-hydroxyprogesterone, androstenedione (A4), T and 17β-estradiol (E2) production. Treatment with T increased production of A4, 11β-hydroxyandrostenedione, 11β-hydroxytestosterone, E2, 11-ketotestosterone, and 5α-dihydrotestosterone (DHT). Interestingly, in the presence of high doses of T the previtellogenic ovary preferentially produced A4 and DHT over any other metabolite. Our data suggest the existence of dopamine inhibition of the reproductive axis in female sablefish. Treatments with Met and T elevated gonadotropin mRNAs in prepubertal females but failed to stimulate the transition into vitellogenic growth, suggesting a possible failure in pituitary gonadotropin protein synthesis/release. Previtellogenic ovarian follicles of sablefish are equipped to synthesize steroids, including those required for vitellogenic growth, and DHT, a steroid hormone whose role in reproduction of fishes remains unknown.

Plasticity of gene-regulatory networks controlling sex determination: of masters, slaves, usual suspects, newcomers, and usurpators

Sexual dimorphism is one of the most pervasive and diverse features of animal morphology, physiology, and behavior. Despite the generality of the phenomenon itself, the mechanisms controlling how sex is determined differ considerably among various organismic groups, have evolved repeatedly and independently, and the underlying molecular pathways can change quickly during evolution. Even within closely related groups of organisms for which the development of gonads on the morphological, histological, and cell biological level is undistinguishable, the molecular control and the regulation of the factors involved in sex determination and gonad differentiation can be substantially different. The biological meaning of the high molecular plasticity of an otherwise common developmental program is unknown. While comparative studies suggest that the downstream effectors of sex-determining pathways tend to be more stable than the triggering mechanisms at the top, it is still unclear how conserved the downstream networks are and how all components work together. After many years of stasis, when the molecular basis of sex determination was amenable only in the few classical model organisms (fly, worm, mouse), recently, sex-determining genes from several animal species have been identified and new studies have elucidated some novel regulatory interactions and biological functions of the downstream network, particularly in vertebrates. These data have considerably changed our classical perception of a simple linear developmental cascade that makes the decision for the embryo to develop as male or female, and how it evolves.

Expression analysis of sex-determining pathway genes during development in male and female Atlantic salmon (Salmo salar)

We studied the expression of 28 genes that are involved in vertebrate sex-determination or sex-differentiation pathways, in male and female Atlantic salmon (Salmo salar) in 11 stages of development from fertilization to after first feeding. Gene expression was measured in half-sibs that shared the same dam. The sire of family 1 was a sex-reversed female (i.e., genetically female but phenotypically male), and so the progeny of this family are all female. The sire of family 2 was a true male, and so the offspring were 50% male and 50% female. Gene expression levels were compared among three groups: 20 female offspring of the cross between a regular female and the sex-reversed female (family 1, first group), ∼ 10 females from the cross between a regular female and a regular male (family 2, second group) and ∼ 10 males from this same family (family 2, third group). Statistically significant differences in expression levels between males and the two groups of females were observed for two genes, gsdf and amh/mis, in the last four developmental stages examined. SdY, the sex-determining gene in rainbow trout, appeared to be expressed in males from 58 days postfertilization (dpf). Starting at 83 dpf, ovarian aromatase, cyp19a, expression appeared to be greater in both groups of females compared with males, but this difference was not statistically significant. The time course of expression suggests that sdY may be involved in the upregulation of gsdf and amh/mis and the subsequent repression of cyp19a in males via the effect of amh/mis.

Comparative Analysis of the Shared Sex-Determination Region (SDR) among Salmonid Fishes

Salmonids present an excellent model for studying evolution of young sex-chromosomes. Within the genus, Oncorhynchus, at least six independent sex-chromosome pairs have evolved, many unique to individual species. This variation results from the movement of the sex-determining gene, sdY, throughout the salmonid genome. While sdY is known to define sexual differentiation in salmonids, the mechanism of its movement throughout the genome has remained elusive due to high frequencies of repetitive elements, rDNA sequences, and transposons surrounding the sex-determining regions (SDR). Despite these difficulties, bacterial artificial chromosome (BAC) library clones from both rainbow trout and Atlantic salmon containing the sdY region have been reported. Here, we report the sequences for these BACs as well as the extended sequence for the known SDR in Chinook gained through genome walking methods. Comparative analysis allowed us to study the overlapping SDRs from three unique salmonid Y chromosomes to define the specific content, size, and variation present between the species. We found approximately 4.1 kb of orthologous sequence common to all three species, which contains the genetic content necessary for masculinization. The regions contain transposable elements that may be responsible for the translocations of the SDR throughout salmonid genomes and we examine potential mechanistic roles of each one.

Genetic variability between complete mitochondrion genomes of the sablefish, Anoplopoma fimbria (Pallas, 1814)

The complete mitochondrial genome of the sablefish, Anoplopoma fimbria (Genbank accession KP777542) is 16,507 bp in size and contains the typical 37 genes (13 protein-coding, 2 ribosomal RNA, and 22 transfer RNA) found in teleosts mitogenomes. The genome varies in 118 positions with respect to another mitogenome sablefish specimen.

Comparative transcriptomics in East African cichlids reveals sex- and species-specific expression and new candidates for sex differentiation in fishes

Males and females of the same species differ largely in gene expression, which accounts for most of the morphological and physiological differences and sex-specific phenotypes. Here, we analyzed sex-specific gene expression in the brain and the gonads of cichlid fishes from Lake Tanganyika belonging to four different lineages, so-called tribes (Eretmodini, Ectodini, Haplochromini, and Lamprologini), using the outgroup Nile tilapia (Oreochromis niloticus) as reference. The comparison between male and female brains revealed few differences between the sexes, consistent in all investigated species. The gonads, on the other hand, showed a large fraction of differentially expressed transcripts with the majority of them showing the same direction of expression in all four species. All here-studied cichlids, especially the three investigated mouth-breeding species, showed a trend toward more male- than female biased transcripts. Transcripts, which were female-biased in expression in all four species, were overrepresented on linkage group (LG)1 in the reference genome and common male-biased transcripts showed accumulation on LG23, the presumable sex chromosomes of the Nile tilapia. Sex-specific transcripts contained candidate genes for sex determination and differentiation in fishes,especially members of the transforming growth factor-b-superfamily and the Wnt-pathway and also prominent members of the sox-, dm-domain-, and high mobility group-box families. We further confirmed our previous finding on species/lineage-specific gene expression shifts in the sex steroid pathway, including synthesizing enzymes as the aromatase cyp19a1 and estrogen and androgen receptors.

A transcriptome derived female-specific marker from the invasive Western mosquitofish (Gambusia affinis)

Sex-specific markers are a prerequisite for understanding reproductive biology, genetic factors involved in sex differences, mechanisms of sex determination, and ultimately the evolution of sex chromosomes. The Western mosquitofish, Gambusia affinis, may be considered a model species for sex-chromosome evolution, as it displays female heterogamety (ZW/ZZ), and is also ecologically interesting as a worldwide invasive species. Here, de novo RNA-sequencing on the gonads of sexually mature G. affinis was used to identify contigs that were highly transcribed in females but not in males (i.e., transcripts with ovary-specific expression). Subsequently, 129 primer pairs spanning 79 contigs were tested by PCR to identify sex-specific transcripts. Of those primer pairs, one female-specific DNA marker was identified, Sanger sequenced and subsequently validated in 115 fish. Sequence analyses revealed a high similarity between the identified sex-specific marker and the 3´ UTR of the aminomethyl transferase (amt) gene of the closely related platyfish (Xiphophorus maculatus). This is the first time that RNA-seq has been used to successfully characterize a sex-specific marker in a fish species in the absence of a genome map. Additionally, the identified sex-specific marker represents one of only a handful of such markers in fishes.

Genetics and timing of sex determination in the East African cichlid fish Astatotilapia burtoni

Background

The factors determining sex are diverse in vertebrates and especially so in teleost fishes. Only a handful of master sex-determining genes have been identified, however great efforts have been undertaken to characterize the subsequent genetic network of sex differentiation in various organisms. East African cichlids offer an ideal model system to study the complexity of sexual development, since many different sex-determining mechanisms occur in closely related species of this fish family. Here, we investigated the sex-determining system and gene expression profiles during male development of Astatotilapia burtoni, a member of the rapidly radiating and exceptionally species-rich haplochromine lineage.

Results

Crossing experiments with hormonally sex-reversed fish provided evidence for an XX-XY sex determination system in A. burtoni. Resultant all-male broods were used to assess gene expression patterns throughout development of a set of candidate genes, previously characterized in adult cichlids only.

Conclusions

We could identify the onset of gonad sexual differentiation at 11–12 dpf. The expression profiles identified wnt4B and wt1A as the earliest gonad markers in A. burtoni. Furthermore we identified late testis genes (cyp19a1A, gsdf, dmrt1 and gata4), and brain markers (ctnnb1A, ctnnb1B, dax1A, foxl2, foxl3, nanos1A, nanos1B, rspo1, sf-1, sox9A and sox9B).

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0140-5) contains supplementary material, which is available to authorized users.

Genetics of sexual development: an evolutionary playground for fish

Teleost fishes are the most species-rich clade of vertebrates and feature an overwhelming diversity of sex-determining mechanisms, classically grouped into environmental and genetic systems. Here, we review the recent findings in the field of sex determination in fish. In the past few years, several new master regulators of sex determination and other factors involved in sexual development have been discovered in teleosts. These data point toward a greater genetic plasticity in generating the male and female sex than previously appreciated and implicate novel gene pathways in the initial regulation of the sexual fate. Overall, it seems that sex determination in fish does not resort to a single genetic cascade but is rather regulated along a continuum of environmental and heritable factors.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.

Transcriptome analysis of the mud crab (Scylla paramamosain) by 454 deep sequencing: assembly, annotation, and marker discovery

In this study, we reported the characterization of the first transcriptome of the mud crab (Scylla paramamosain). Pooled cDNAs of four tissue types from twelve wild individuals were sequenced using the Roche 454 FLX platform. Analysis performed included de novo assembly of transcriptome sequences, functional annotation, and molecular marker discovery. A total of 1,314,101 high quality reads with an average length of 411 bp were generated by 454 sequencing on a mixed cDNA library. De novo assembly of these 1,314,101 reads produced 76,778 contigs (consisting of 818,154 reads) with 5.4-fold average sequencing coverage. The remaining 495,947 reads were singletons. A total of 78,268 unigenes were identified based on sequence similarity with known proteins (E≤0.00001) in UniProt and non-redundant protein databases. Meanwhile, 44,433 sequences were identified (E≤0.00001) using a BLASTN search against the NCBI nucleotide database. Gene Ontology (GO) analysis indicated that biosynthetic process, cell part, and ion binding were the most abundant terms in biological process, cellular component, and molecular function categories, respectively. Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analysis revealed that 4,878 unigenes distributed in 281 different pathways. In addition, 19,011 microsatellites and 37,063 potential single nucleotide polymorphisms were detected from the transcriptome of S. paramamosain. Finally, thirty polymorphic microsatellite markers were developed and used to assess genetic diversity of a wild population of S. paramamosain. So far, existing sequence resources for S. paramamosain are extremely limited. The present study provides a characterization of transcriptome from multiple tissues and individuals, as well as an assessment of genetic diversity of a wild population. These sequence resources will facilitate the investigation of population genetic diversity, the development of genetic maps, and the conduct of molecular marker-assisted breeding in S. paramamosain and related crab species.

Co-option of Sox3 as the male-determining factor on the Y chromosome in the fish Oryzias dancena

Sex chromosomes harbour a primary sex-determining signal that triggers sexual development of the organism. However, diverse sex chromosome systems have been evolved in vertebrates. Here we use positional cloning to identify the sex-determining locus of a medaka-related fish, Oryzias dancena, and find that the locus on the Y chromosome contains a cis-regulatory element that upregulates neighbouring Sox3 expression in developing gonad. Sex-reversed phenotypes in Sox3(Y) transgenic fish, and Sox3(Y) loss-of-function mutants all point to its critical role in sex determination. Furthermore, we demonstrate that Sox3 initiates testicular differentiation by upregulating expression of downstream Gsdf, which is highly conserved in fish sex differentiation pathways. Our results not only provide strong evidence for the independent recruitment of Sox3 to male determination in distantly related vertebrates, but also provide direct evidence that a novel sex determination pathway has evolved through co-option of a transcriptional regulator potentially interacted with a conserved downstream component.