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

Climate change can impair bacterial pathogen defences in sablefish via hypoxia-mediated effects on adaptive immunity

Low-oxygen levels (hypoxia) in aquatic habitats are becoming more common because of global warming and eutrophication. However, the effects on the health/disease status of fishes, the world's largest group of vertebrates, are unclear. Therefore, we assessed how long-term hypoxia affected the immune function of sablefish, an ecologically and economically important North Pacific species, including the response to a formalin-killed Aeromonas salmonicida bacterin. Sablefish were held at normoxia or hypoxia (100% or 40% air saturated seawater, respectively) for 6-16 weeks, while we measured a diverse array of immunological traits. Given that the sablefish is a non-model organism, this involved the development of a species-specific methodological toolbox comprised of qPCR primers for 16 key immune genes, assays for blood antibacterial defences, the assessment of blood immunoglobulin (IgM) levels with ELISA, and flow cytometry and confocal microscopy techniques. We show that innate immune parameters were typically elevated in response to the bacterial antigens, but were not substantially affected by hypoxia. In contrast, hypoxia completely prevented the ∼1.5-fold increase in blood IgM level that was observed under normoxic conditions following bacterin exposure, implying a serious impairment of adaptive immunity. Since the sablefish is naturally hypoxia tolerant, our results demonstrate that climate change-related deoxygenation may be a serious threat to the immune competency of fishes.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo)

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Speciation genomics and the role of depth in the divergence of rockfishes ( Sebastes ) revealed through Pool‐seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth‐segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth‐segregated and have different degrees of divergence: S. chlorostictus–S. rosenblatti and S. crocotulus–S. miniatus. We were able to reliably identify “islands of divergence” in the species pair with more recent divergence (S. chlorostictus–S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Heterochiasmy and the establishment of gsdf as a novel sex determining gene in Atlantic halibut

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly from a male and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male-only and female-only meiotic recombination regions (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4.5 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. Our findings indicate that heterochiasmy may facilitate the evolution of genetic sex determination systems relying on linkage of sexually antagonistic loci to a sex-determining factor.

Even closely related fish species can have different sex chromosomes, but this turn-over of sex determination systems is poorly understood. Here, we used large-scale genome sequencing to determine the DNA sequence of the Atlantic halibut chromosomes and compared sequencing data from males and females to identify the sex chromosomes. We show that males have much higher gene activity of the gene gonadal somatic cell derived factor (gsdf), which is located on the sex chromosomes and has a role in testicular development. The genome contains many mobile DNA sequences, transposable elements (TEs), one placed in front of gsdf, enhancing its activity. This made gsdf the sex determining factor, thereby creating a new Y-chromosome. We further describe how all Atlantic halibut chromosomes behave similar to sex chromosomes in that most regions only recombine in one sex. This phenomenon may contribute to the rapid turn-over of genetic sex determination systems in fish. Our results highlight the molecular events creating a new Y-chromosome and show that the new Atlantic halibut Y was formed less than 4.5 million years ago. Future studies in Atlantic halibut and closely related species can shed light on mechanisms contributing to sex chromosome evolution in fish.

A Chromosome-Level Genome Assembly of Yellowtail Kingfish (Seriola lalandi)

Yellowtail kingfish (Seriola lalandi) is a pelagic marine piscivore with a circumglobal distribution. It is particularly suitable for open ocean aquaculture owing to its large body size, fast swimming, rapid growth, and high economic value. A high-precision genome is of great significance for future genetic breeding research and large-scale aquaculture in the open ocean. PacBio, Illumina, and Hi-C data were combined to assemble chromosome-level reference genome with the size of 648.34 Mb (contig N50: 28.52 Mb). 175 contigs was anchored onto 24 chromosomes with lengths ranging from 12.28 to 34.59 Mb, and 99.79% of the whole genome sequence was covered. The BUSCOs of genome and gene were 94.20 and 95.70%, respectively. Gene families associated with adaptive behaviors, such as olfactory receptors and HSP70 gene families, expanded in the genome of S. lalandi. An analysis of selection pressure revealed 652 fast-evolving genes, among which mkxb, popdc2, dlx6, and ifitm5 may be related to rapid growth traits. The data generated in this study provide a valuable resource for understanding the genetic basis of S. lalandi traits.

Reference genome of Lumpfish Cyclopterus lumpus Linnaeus provides evidence of male heterogametic sex determination through the AMH pathway

Teleosts exhibit extensive diversity of sex determination (SD) systems and mechanisms, providing the opportunity to study the evolution of sex determination and sex chromosomes. Here we sequenced the genome of the Common Lumpfish (Cyclopterus lumpus Linnaeus), a species of increasing importance to aquaculture, and identified the SD region and master SD locus using a 70K SNP array and tissue-specific expression data. The chromosome-level assembly identified 25 diploid chromosomes with a total size of 572.89 Mb, a scaffold N50 of 23.86 Mb, and genome annotation predicted 21,480 protein-coding genes. Genome wide association analysis located a highly sex-associated region on chromosome 13, suggesting that anti-Müllerian hormone (AMH) is the putative SD factor. Linkage disequilibrium and heterozygosity across chromosome 13 support a proto-XX/XY system, with an absence of widespread chromosome divergence between sexes. We identified three copies of AMH in the Lumpfish primary and alternate haplotype assemblies localized in the SD region. Comparison to sequences from other teleosts suggested a monophyletic relationship and conservation within the Cottioidei. One AMH copy showed similarity to AMH/AMHY in a related species and was also the only copy with expression in testis tissue, suggesting this copy may be the functional copy of AMH in Lumpfish. The two other copies arranged in tandem inverted duplication were highly similar, suggesting a recent duplication event. This study provides a resource for the study of early sex chromosome evolution and novel genomic resources that benefits Lumpfish conservation management and aquaculture.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

Lessons from an unusual vertebrate sex-determining gene

So far, very few sex-determining genes have been identified in vertebrates and most of them, the so-called 'usual suspects', evolved from genes which fulfil essential functions during sexual development and are thus already tightly linked to the process that they now govern. The single exception to this 'usual suspects' rule in vertebrates so far is the conserved salmonid sex-determining gene, sdY (sexually dimorphic on the Y chromosome), that evolved from a gene known to be involved in regulation of the immune response. It is contained in a jumping sex locus that has been transposed or translocated into different ancestral autosomes during the evolution of salmonids. This special feature of sdY, i.e. being inserted in a 'jumping sex locus', could explain how salmonid sex chromosomes remain young and undifferentiated to escape degeneration. Recent knowledge on the mechanism of action of sdY demonstrates that it triggers its sex-determining action by deregulating oestrogen synthesis that is a conserved and crucial pathway for ovarian differentiation in vertebrates. This result suggests that sdY has evolved to cope with a pre-existing sex differentiation regulatory network. Therefore, 'limited options' for the emergence of new master sex-determining genes could be more constrained by their need to tightly interact with a conserved sex differentiation regulatory network rather than by being themselves 'usual suspects', already inside this sex regulatory network. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Evolution of master sex determiners: TGF-β signalling pathways at regulatory crossroads

To date, more than 20 different vertebrate master sex-determining genes have been identified on different sex chromosomes of mammals, birds, frogs and fish. Interestingly, six of these genes are transcription factors (Dmrt1- or Sox3- related) and 13 others belong to the TGF-β signalling pathway (Amh, Amhr2, Bmpr1b, Gsdf and Gdf6). This pattern suggests that only a limited group of factors/signalling pathways are prone to become top regulators again and again. Although being clearly a subordinate member of the sex-regulatory network in mammals, the TGF-β signalling pathway made it to the top recurrently and independently. Facing this rolling wave of TGF-β signalling pathways, this review will decipher how the TGF-β signalling pathways cope with the canonical sex gene regulatory network and challenge the current evolutionary concepts accounting for the diversity of sex-determining mechanisms. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Evolution, Expression, and Function of Gonadal Somatic Cell-Derived Factor

Fish gonads develop in very diverse ways different from mammalian gonads. This diversity is contributed by species-specific factors. Gonadal somatic cell-derived factor (Gsdf) is one such factor. The gsdf gene exists mostly in teleosts and is absent in many tetrapods, probably as a result of two gene losses during evolution. The gsdf transcript is expressed mainly in gonadal somatic cells, including Sertoli cell in testis and granulosa cells in ovary; however, these gonadal somatic cells can surround many types of germ cells at different developmental stages depending on the fish species. The function of gsdf is also variable. It is involved in germ cell proliferation, testicular formation, ovarian development and even male sex determination. Here, we summarize the common and diverse expression, regulation and functions of gsdf among different fish species with aspect of evolution.

Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY. By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus

The genetic regulatory network governing vertebrate gonadal differentiation appears less conserved than previously thought. Here, we investigated the gonadal development of Astyanax mexicanus Pachón cavefish by looking at primordial germ cells (PGCs) migration and proliferation, gonad histology, and gene expression patterns. We showed that PGCs are first detected at the 80% epiboly stage and then reach the gonadal primordium at 1 day post-fertilization (dpf). However, in contrast to the generally described absence of PGCs proliferation during their migration phase, PGCs number in cavefish doubles between early neurula and 8-9 somites stages. Combining both gonadal histology and vasa (germ cell marker) expression patterns, we observed that ovarian and testicular differentiation occurs around 65 dpf in females and 90 dpf in males, respectively, with an important inter-individual variability. The expression patterns of dmrt1, gsdf, and amh revealed a conserved predominant male expression during cavefish gonadal development, but none of the ovarian differentiation genes, i. e., foxl2a, cyp19a1a, and wnt4b displayed an early sexually dimorphic expression, and surprisingly all these genes exhibited predominant expression in adult testes. Altogether, our results lay the foundation for further research on sex determination and differentiation in A. mexicanus and contribute to the emerging picture that the vertebrate sex differentiation downstream regulatory network is less conserved than previously thought, at least in teleost fishes.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with XX/XY sex determination

Phenotypic sex of an organism is determined by molecular changes in the gonads, so-called molecular sex differentiation, which should precede the rise of cellular or anatomical sex-distinguishing features. This study characterized molecular and morphological sex differentiation in sablefish (Anoplopoma fimbria), a marine teleost with established XX/XY genotypic sex determination. Next generation sequencing was conducted on sablefish ovarian and testicular mRNAs to obtain sequences for transcripts associated with vertebrate sex determination and differentiation and early reproductive development. Gene-specific PCRs were developed to determine the distribution and ontogenetic gonadal expression of transcription, growth, steroidogenic and germline factors, as well as gonadotropin and steroid receptors. Molecular changes associated with sex differentiation were first apparent in both XY- and XX-genotype sablefish at ~ 60 mm in body length and prior to histological signs of sex differentiation. The earliest and most robust markers of testicular differentiation were gsdf, amh, dmrt1, cyp11b, star, sox9a, and fshr. Markedly elevated mRNA levels of several steroidogenesis-related genes and ar2 in differentiating testes suggested that androgens play a role in sablefish testicular differentiation. The earliest markers of ovarian differentiation were cyp19a1a, lhcgr, foxl2, nr0b1, and igf3. Other transcripts such as figla, zp3, and pou5f3 were expressed predominantly in XX-genotype fish and significantly increased with the first appearance and subsequent development of primary oocytes. This study provides valuable insight to the developmental sequence of events associated with gonadal sex differentiation in marine teleosts with XX/XY sex determination. It also implicates particular genes in processes of male and female development and establishes robust molecular markers for phenotypic sex in sablefish, useful for ongoing work related to sex control and reproductive sterilization.

Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758)

BACKGROUND

Goldfish is an important model for various areas of research, including neural development and behavior and a species of significant importance in aquaculture, especially as an ornamental species. It has a male heterogametic (XX/XY) sex determination system that relies on both genetic and environmental factors, with high temperatures being able to produce female-to-male sex reversal. Little, however, is currently known on the molecular basis of genetic sex determination in this important cyprinid model. Here we used sequencing approaches to better characterize sex determination and sex-chromosomes in an experimental strain of goldfish.

RESULTS

Our results confirmed that sex determination in goldfish is a mix of environmental and genetic factors and that its sex determination system is male heterogametic (XX/XY). Using reduced representation (RAD-seq) and whole genome (pool-seq) approaches, we characterized sex-linked polymorphisms and developed male specific genetic markers. These male specific markers were used to distinguish sex-reversed XX neomales from XY males and to demonstrate that XX female-to-male sex reversal could even occur at a relatively low rearing temperature (18 °C), for which sex reversal has been previously shown to be close to zero. We also characterized a relatively large non-recombining region (~ 11.7 Mb) on goldfish linkage group 22 (LG22) that contained a high-density of male-biased genetic polymorphisms. This large LG22 region harbors 373 genes, including a single candidate as a potential master sex gene, i.e., the anti-Mullerian hormone gene (amh). However, no sex-linked polymorphisms were detected in the coding DNA sequence of the goldfish amh gene.

CONCLUSIONS

These results show that our goldfish strain has a relatively large sex locus on LG22, which is likely the Y chromosome of this experimental population. The presence of a few XX males even at low temperature also suggests that other environmental factors in addition to temperature could trigger female-to-male sex reversal. Finally, we also developed sex-linked genetic markers, which will be important tools for future research on sex determination in our experimental goldfish population. However, additional work would be needed to explore whether this sex locus is conserved in other populations of goldfish.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

The complete mitochondrial genome sequencing of Anisakis simplex isolated from Anoplopoma fimbria

In this study, zoonotic Anisakis simplex was isolated and identified from the outermost layer of the stomach of diseased Anoplopoma fimbria at an industrial farm in Liaoning, North China (122.1842 E, 39.2616 N). With the completion of A. simplex mitochondrial genome sequencing, the full-length mitochondrial genome of A. simplex was assembled and analyzed. All results indicate that the complete mitochondrial genome of A. simplex was 13,899 bp. There were 20 tRNAs and 12 protein-coding genes (PCGs), and two rRNA all located at the heavy (H) strand. Besides, the phylogenetic tree of 19 A. simplex isolated from different host species was constructed. The results showed that A. simplex isolated from A. fimbria was clustered with Oncorhynchus nerka isolate in a clade. To sum up, our research results would further provide essential data for systematics and evolution study of A. simplex.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

The Genetics and Epigenetics of Sex Change in Fish

Fish show extraordinary sexual plasticity, changing sex naturally as part of their life cycle or reversing sex because of environmental stressors. This plasticity shows that sexual fate is not an irreversible process but the result of an ongoing tug-of-war for supremacy between male and female signaling networks. The behavioral, gonadal, and morphological changes involved in this process are well described, yet the molecular events that underpin those changes remain poorly understood. Epigenetic modifications emerge as a critical link between environmental stimuli, the onset of sex change, and subsequent maintenance of sexual phenotype. Here we synthesize current knowledge of sex change, focusing on the genetic and epigenetic processes that are likely involved in the initiation and regulation of sex change. We anticipate that better understanding of sex change in fish will shed new light on sex determination and development in vertebrates and on how environmental perturbations affect sexual fate.

Inhibition of ovarian development and instances of sex reversal in genotypic female sablefish (Anoplopoma fimbria) exposed to elevated water temperature

This study determined high temperature effects on ovarian development in a marine groundfish species, sablefish (Anoplopoma fimbria), with potential application in sex reversal or sterilization for aquaculture. Monosex female (XX-genotype) sablefish larvae (∼30 mm) were randomly divided into three groups and exposed to control (15.6 °C ± 0.8 °C), moderate (20.4 °C ± 0.5 °C), or high (21.7 °C ± 0.5 °C) temperatures for 19 weeks. Treated fish were then tagged and transferred to ambient seawater (11.2 °C ± 2.3 °C) for one year to determine whether temperature effects on reproductive development were maintained post-treatment. Fish were periodically sampled for gonadal histology, gene expression and plasma 17β-estradiol (E2) analyses to assess gonadal development. Short-term (4-week) exposure to elevated temperatures had only minor effects, whereas longer exposure (12-19 weeks) markedly inhibited ovarian development. Fish from the moderate and high treatment groups had significantly less developed ovaries relative to controls, and mRNA levels for germ cell (vasa, zpc) and apoptosis-associated genes (p53, casp8) generally indicated gonadal degeneration. The high treatment group also had significantly reduced plasma E2 levels and elevated gonadal amh gene expression. After one year at ambient temperatures, however, ovaries of moderate and high treatment fish exhibited compensatory recovery and were indistinguishable from controls. Two genotypic females possessing immature testes (neomales) were observed in the high treatment group, indicating sex reversal had occurred (6% rate). These results demonstrate that extreme elevated temperatures may inhibit ovarian development or trigger sex reversal. High temperature treatment is likely not an effective sterilization method but may be preferable for sablefish neomale broodstock production.

Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio)

Over the last several decades zebrafish (Danio rerio) has become a major model organism for the study of vertebrate development and physiology. Given this, it may be surprising how little is known about the mechanism that zebrafish use to determine sex. While zebrafish are a gonochoristic species (having two sexes) that do not switch sex as adults, it was appreciated early on that sex ratios obtained from breeding lab domesticated lines were not typically a 1:1 ratio of male and female, suggesting that sex was not determined by a strict chromosomal mechanism. Here we will review the recent progress toward defining the genetic mechanism for sex determination in both wild and domesticated zebrafish.

Repeated Evolution Versus Common Ancestry: Sex Chromosome Evolution in the Haplochromine Cichlid Pseudocrenilabrus philander

Why sex chromosomes turn over and remain undifferentiated in some taxa, whereas they degenerate in others, is still an area of ongoing research. The recurrent occurrence of homologous and homomorphic sex chromosomes in distantly related taxa suggests their independent evolution or continued recombination since their first emergence. Fishes display a great diversity of sex-determining systems. Here, we focus on sex chromosome evolution in haplochromines, the most species-rich lineage of cichlid fishes. We investigate sex-specific signatures in the Pseudocrenilabrus philander species complex, which belongs to a haplochromine genus found in many river systems and ichthyogeographic regions in northern, eastern, central, and southern Africa. Using whole-genome sequencing and population genetic, phylogenetic, and read-coverage analyses, we show that one population of P. philander has an XX-XY sex-determining system on LG7 with a large region of suppressed recombination. However, in a second bottlenecked population, we did not find any sign of a sex chromosome. Interestingly, LG7 also carries an XX-XY system in the phylogenetically more derived Lake Malawi haplochromine cichlids. Although the genomic regions determining sex are the same in Lake Malawi cichlids and P. philander, we did not find evidence for shared ancestry, suggesting that LG7 evolved as sex chromosome at least twice in haplochromine cichlids. Hence, our work provides further evidence for the labile nature of sex determination in fishes and supports the hypothesis that the same genomic regions can repeatedly and rapidly be recruited as sex chromosomes in more distantly related lineages.

Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod

The genetic mechanisms determining sex in teleost fishes are highly variable and the master sex determining gene has only been identified in few species. Here we characterize a male-specific region of 9 kb on linkage group 11 in Atlantic cod (Gadus morhua) harboring a single gene named zkY for zinc knuckle on the Y chromosome. Diagnostic PCR test of phenotypically sexed males and females confirm the sex-specific nature of the Y-sequence. We identified twelve highly similar autosomal gene copies of zkY, of which eight code for proteins containing the zinc knuckle motif. 3D modeling suggests that the amino acid changes observed in six copies might influence the putative RNA-binding specificity. Cod zkY and the autosomal proteins zk1 and zk2 possess an identical zinc knuckle structure, but only the Y-specific gene zkY was expressed at high levels in the developing larvae before the onset of sex differentiation. Collectively these data suggest zkY as a candidate master masculinization gene in Atlantic cod. PCR amplification of Y-sequences in Arctic cod (Arctogadus glacialis) and Greenland cod (Gadus macrocephalus ogac) suggests that the male-specific region emerged in codfishes more than 7.5 million years ago.

Genetic regulation of sex determination and maintenance in zebrafish (Danio rerio)

Over the last several decades zebrafish (Danio rerio) has become a major model organism for the study of vertebrate development and physiology. Given this, it may be surprising how little is known about the mechanism that zebrafish use to determine sex. While zebrafish are a gonochoristic species (having two sexes) that do not switch sex as adults, it was appreciated early on that sex ratios obtained from breeding lab domesticated lines were not typically a 1:1 ratio of male and female, suggesting that sex was not determined by a strict chromosomal mechanism. Here we will review the recent progress toward defining the genetic mechanism for sex determination in both wild and domesticated zebrafish.

The autosomal Gsdf gene plays a role in male gonad development in Chinese tongue sole (Cynoglossus semilaevis)

Gsdf is a key gene for testicular differentiation in teleost. However, little is known about the function of Gsdf in Chinese tongue sole (Cynoglossus semilaevis). In this study, we obtained the full-length Gsdf gene (CS-Gsdf), and functional characterization revealed its potential participation during germ cell differentiation in testes. CS-Gsdf transcription was predominantly detected in gonads, while the levels in testes were significantly higher than those in ovaries. During the different developmental stages in male gonads, the mRNA level was significantly upregulated at 86 dph, and a peak appeared at 120 dph; then, the level decreased at 1 and 2 yph. In situ hybridization revealed that CS-Gsdf mRNA was mainly localized in the Sertoli cells, spermatogonia, and spermatids in mature testes. After CS-Gsdf knockdown in the male testes cell line by RNA interference, a series of sex-related genes was influenced, including several sex differentiation genes, CS-Wnt4a, CS-Cyp19a1a and CS-Star. Based on these data, we speculated that CS-Gsdf may play a positive role in germ differentiation and proliferation via influencing genes related to sex differentiation.

Diverse and variable sex determination mechanisms in vertebrates

Sex is prevalent in nature and sex determination is one of the most fundamental biological processes, while the way of initiating female and male development exhibits remarkable diversity and variability across vertebrates. The knowledge on why and how sex determination mechanisms evolve unusual plasticity remains limited. Here, we summarize sex determination systems, master sex-determining genes and gene-regulatory networks among vertebrates. Recent research advancements on sex determination system transition are also introduced and discussed in some non-model animals with multiple sex determination mechanisms. This review will provide insights into the origin, transition and evolutionary adaption of different sex determination strategies in vertebrates, as well as clues for future perspectives in this field.

Insights into teleost sex determination from the Seriola dorsalis genome assembly

Background

The assembly and annotation of a genome is a valuable resource for a species, with applications ranging from conservation genomics to gene discovery. Genomic resource development is especially important for species in culture, such as the California Yellowtail (Seriola dorsalis), the likely candidate for the establishment of commercial offshore aquaculture production in southern California. Genomic resource development for this species will improve the understanding of sex and other phenotypic traits, and allow for rapid increases in genetic improvement for and economic gain in culture production.

Results

We describe the assembly and annotation of the S. dorsalis genome, and present resequencing data from 45 male and 45 female wild-caught S. dorsalis used to identify a sex-determining region and marker in this species. The genome assembly captured approximately 93% of the total 685 MB genome with an average coverage depth of 180×. Using the assembled genome, resequencing data from the 90 fish were aligned to place boundaries on the sex-determining region. Sex-specific markers were developed based on a female-specific, 61 nucleotide deletion identified in that region. We hypothesize that Estradiol 17-beta-dehydrogenase is the putative sex-determining gene and propose a plausible genetic mechanism for ZW sex determination in S. dorsalis involving a female-specific deletion of a transcription factor binding motif that may be targeted by Sox3.

Conclusions

Understanding the mechanism of sex determination and development of assays to determine sex is critical both for management of wild fisheries and for development of efficient and sustainable aquaculture practices. In addition, this genome assembly for S. dorsalis will be a substantial resource for a variety of future research applications.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4403-1) contains supplementary material, which is available to authorized users.

Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish

BACKGROUND

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

RESULTS

Zebrafish expressed gsdf in somatic cells of bipotential gonads and expression continued in ovarian granulosa cells and testicular Sertoli cells. Homozygous gsdf knockout mutants delayed leaving the bipotential gonad state, but then became a male or a female. Mutant females ovulated a few oocytes, then became sterile, accumulating immature follicles. Female mutants stored excess lipid and down-regulated aromatase, gata4, insulin receptor, estrogen receptor, and genes for lipid metabolism, vitellogenin, and steroid biosynthesis. Mutant females contained less estrogen and more androgen than wild-types. Mutant males were fertile. Genomic analysis suggests that Gsdf, Bmp15, and Gdf9, originated as paralogs in vertebrate genome duplication events.

CONCLUSIONS

In zebrafish, gsdf regulates ovarian follicle maturation and expression of genes for steroid biosynthesis, obesity, diabetes, and female fertility, leading to ovarian and extra-ovarian phenotypes that mimic human polycystic ovarian syndrome (PCOS), suggesting a role for a related TGFB signaling molecule in the etiology of PCOS. Developmental Dynamics 246:925-945, 2017. © 2017 Wiley Periodicals, Inc.

Growth performance, oxidative stress, and non-specific immune responses in juvenile sablefish, Anoplopoma fimbria, by changes of water temperature and salinity

Juvenile sablefish, Anoplopoma fimbria (mean length 15.5 ± 1.9 cm, mean weight 68.5 ± 4.8 g), were used to evaluate the effects on growth, oxidative stress, and non-specific immune responses by changes of water temperature (8, 10, 12, 14, 16, 18, and 20 °C) and salinity (100 (35.0), 90 (31.5), 80 (28.0), 70 (24.5), 60 (21.0), 50 (17.5), and 40% (14.0) (‰)) for 4 months. The growth performance was significantly increased at the temperature of 12 and 14 °C, and the feed efficiency was notably decreased at the temperature of 18 °C. The growth performance and feed efficiency were also significantly decreased at low salinity. The antioxidant responses such as superoxide dismutase and catalase were significantly increased by the high temperature and decreased by the low salinity. The immune responses such as lysozyme and phagocytosis were elevated by the temperature of 18 °C and decreased by the salinity of 50%. The results of this study indicate that the growth performance of juvenile sablefish, A. fimbria, is influenced by the temperature and salinity, and the excessive temperature and salinity levels can affect the antioxidant and immune responses.

Sexually dimorphic expression in developing and adult gonads shows an important role of gonadal soma-derived factor during sex differentiation in olive flounder (Paralichthys olivaceus)

Gonadal soma-derived factor (gsdf) is a new member of transforming growth factor beta (TGF-β) superfamily. As a teleost- and gonad-specific growth factor, gsdf has been indicated to play an important role in early germ cell development. However, little is known about its significance in germ cell development of olive flounder (Paralichthys olivaceus). In the present study, a 1338 bp gsdf gene was isolated from P. olivaceus for the first time. Bioinformatic analysis revealed that the genomic structure and synteny relationship of gsdf in teleosts were conserved. Quantitative real-time PCR (qRT-PCR) showed that gsdf expressed before sex gonadal differentiation, and the expression level increased rapidly after initiation of sex differentiation in males. In adult individuals, the expression of gsdf was higher in testis than that in ovary (P<0.01). In situ hybridization (ISH) indicated that gsdf mRNA was detected in the somatic cells of both males and females, and also in the cytoplasm of oocytes. These results suggested that gsdf might play an important role as initial switches to promote testis differentiation and participate in early germ cell development, such as proliferation and differentiation of spermatogonia and oogonia in P. olivaceus.

Large-scale SNP discovery and construction of a high-density genetic map of Colossoma macropomum through genotyping-by-sequencing

Colossoma macropomum, or tambaqui, is the largest native Characiform species found in the Amazon and Orinoco river basins, yet few resources for genetic studies and the genetic improvement of tambaqui exist. In this study, we identified a large number of single-nucleotide polymorphisms (SNPs) for tambaqui and constructed a high-resolution genetic linkage map from a full-sib family of 124 individuals and their parents using the genotyping by sequencing method. In all, 68,584 SNPs were initially identified using minimum minor allele frequency (MAF) of 5%. Filtering parameters were used to select high-quality markers for linkage analysis. We selected 7,734 SNPs for linkage mapping, resulting in 27 linkage groups with a minimum logarithm of odds (LOD) of 8 and maximum recombination fraction of 0.35. The final genetic map contains 7,192 successfully mapped markers that span a total of 2,811 cM, with an average marker interval of 0.39 cM. Comparative genomic analysis between tambaqui and zebrafish revealed variable levels of genomic conservation across the 27 linkage groups which allowed for functional SNP annotations. The large-scale SNP discovery obtained here, allowed us to build a high-density linkage map in tambaqui, which will be useful to enhance genetic studies that can be applied in breeding programs.

A second generation SNP and SSR integrated linkage map and QTL mapping for the Chinese mitten crab Eriocheir sinensis

The Chinese mitten crab Eriocheir sinensis is the most economically important cultivated crab species in China, and its genome has a high number of chromosomes (2n = 146). To obtain sufficient markers for construction of a dense genetic map for this species, we employed the recently developed specific-locus amplified fragment sequencing (SLAF-seq) method for large-scale SNPs screening and genotyping in a F1 full-sib family of 149 individuals. SLAF-seq generated 127,677 polymorphic SNP markers, of which 20,803 valid markers were assigned into five segregation types and were used together with previous SSR markers for linkage map construction. The final integrated genetic map included 17,680 SNP and 629 SSR markers on the 73 linkage groups (LG), and spanned 14,894.9 cM with an average marker interval of 0.81 cM. QTL mapping localized three significant growth-related QTL to a 1.2 cM region in LG53 as well as 146 sex-linked markers in LG48. Genome-wide QTL-association analysis further identified four growth-related QTL genes named LNX2, PAK2, FMRFamide and octopamine receptors. These genes are involved in a variety of different signaling pathways including cell proliferation and growth. The map and SNP markers described here will be a valuable resource for the E. sinensis genome project and selective breeding programs.

The Reversible Sex of Gonochoristic Fish: Insights and Consequences

Fish sex reversal is a means to understand sex determination and differentiation, but it is also used to control sex in aquaculture. This review discusses sex reversal in gonochoristic fish, with the coexistence of genetic and environmental influences. The different periods of fish sensitivity to sex reversal treatments are presented with the mechanisms implicated. The old players of sex differentiation are revisited with transcriptome data and loss of function studies following hormone- or temperature-induced sex reversal. We also discuss whether cortisol is the universal mediator of sex reversal in fish due to its implication in ovarian meiosis and 11KT increase. The large plasticity in fish for sex reversal is also evident in the brain, with a reversibility existing even in adulthood. Studies on epigenetics are presented, since it links the environment, gene expression, and sex reversal, notably the association of DNA methylation in sex reversal. Manipulations with exogenous factors reverse the primary sex in many fish species under controlled conditions, but several questions arise on whether this can occur under wild conditions and what is the ecological significance. Cases of sex reversal in wild fish populations are shown and their fitness and future perspectives are discussed.

Variations on a theme: Genomics of sex determination in the cichlid fish Astatotilapia burtoni

Background

Sex chromosomes change more frequently in fish than in mammals or birds. However, certain chromosomes or genes are repeatedly used as sex determinants in different members of the teleostean lineage. East African cichlids are an enigmatic model system in evolutionary biology representing some of the most diverse extant vertebrate adaptive radiations. How sex is determined and if different sex-determining mechanisms contribute to speciation is unknown for almost all of the over 1,500 cichlid species of the Great Lakes. Here, we investigated the genetic basis of sex determination in a cichlid from Lake Tanganyika, Astatotilapia burtoni, a member of the most species-rich cichlid lineage, the haplochromines.

Results

We used RAD-sequencing of crosses for two populations of A. burtoni, a lab strain and fish caught at the south of Lake Tanganyika. Using association mapping and comparative genomics, we confirmed male heterogamety in A. burtoni and identified different sex chromosomes (LG5 and LG18) in the two populations of the same species. LG5, the sex chromosome of the lab strain, is a fusion chromosome in A. burtoni. Wnt4 is located on this chromosome, representing the best candidate identified so far for the master sex-determining gene in our lab strain of A. burtoni.

Conclusions

Cichlids exemplify the high turnover rate of sex chromosomes in fish with two different chromosomes, LG5 and LG18, containing major sex-determining loci in the two populations of A. burtoni examined here. However, they also illustrate that particular chromosomes are more likely to be used as sex chromosomes. Chromosome 5 is such a chromosome, which has evolved several times as a sex chromosome, both in haplochromine cichlids from all Great Lakes and also in other teleost fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3178-0) contains supplementary material, which is available to authorized users.

Sexually dimorphic gene expressions in eels: useful markers for early sex assessment in a conservation context

Environmental sex determination (ESD) has been detected in a range of vertebrate reptile and fish species. Eels are characterized by an ESD that occurs relatively late, since sex cannot be histologically determined before individuals reach 28 cm. Because several eel species are at risk of extinction, assessing sex at the earliest stage is a crucial management issue. Based on preliminary results of RNA sequencing, we targeted genes susceptible to be differentially expressed between ovaries and testis at different stages of development. Using qPCR, we detected testis-specific expressions of dmrt1, amh, gsdf and pre-miR202 and ovary-specific expressions were obtained for zar1, zp3 and foxn5. We showed that gene expressions in the gonad of intersexual eels were quite similar to those of males, supporting the idea that intersexual eels represent a transitional stage towards testicular differentiation. To assess whether these genes would be effective early molecular markers, we sampled juvenile eels in two locations with highly skewed sex ratios. The combined expression of six of these genes allowed the discrimination of groups according to their potential future sex and thus this appears to be a useful tool to estimate sex ratios of undifferentiated juvenile eels.

A Diverse Family of Host-Defense Peptides (Piscidins) Exhibit Specialized Anti-Bacterial and Anti-Protozoal Activities in Fishes

Conventional antibiotics and other chemical-based drugs are currently one of the most common methods used to control disease-related mortality in animal agriculture. Use of the innate immune system to decrease disease related mortalities is a novel alternative to conventional drugs. One component of the innate immune system is the host-defense peptides, also known as antimicrobial peptides. Host-defense peptides are typically small, amphipathic, α-helical peptides with a broad-spectrum of action against viral, bacterial, fungal, and/or protozoal pathogens. Piscidins are host-defense peptides first discovered in the hybrid striped bass (white bass, Morone chrysops, x striped bass, M. saxatilis). In this paper we identify four new piscidin isoforms in the hybrid striped bass and describe their tissue distributions. We also determine the progenitor species of origin of each piscidin (orthology) and propose a revised nomenclature for this newly described piscidin family based on a three class system. The Class I piscidins (22 amino acids in length; striped bass and white bass piscidin 1 and piscidin 3) show broad-spectrum activity against bacteria and ciliated protozoans, while the Class III piscidins (55 amino acids in length; striped bass and white bass piscidin 6 and striped bass piscidin 7) primarily show anti-protozoal activity. The Class II piscidins (44-46 amino acids in length; striped bass and white bass piscidin 4 and white bass piscidin 5) have a level of activity against bacteria and protozoans intermediate to Classes I and III. Knowledge of piscidin function and activity may help in the future development of disease-resistant lines of striped bass and white bass that could be used to produce superior hybrids for aquaculture.

Vertebrate sex-determining genes play musical chairs

Sexual reproduction is one of the most highly conserved processes in evolution. However, the genetic and cellular mechanisms making the decision of whether the undifferentiated gonad of animal embryos develops either towards male or female are manifold and quite diverse. In vertebrates, sex-determining mechanisms range from environmental to simple or complex genetic mechanisms and different mechanisms have evolved repeatedly and independently. In species with simple genetic sex-determination, master sex-determining genes lying on sex chromosomes drive the gonadal differentiation process by switching on a developmental program, which ultimately leads to testicular or ovarian differentiation. So far, very few sex-determining genes have been identified in vertebrates and apart from mammals and birds, these genes are apparently not conserved over a larger number of related orders, families, genera, or even species. To fill this knowledge gap and to better explore genetic sex-determination, we propose a strategy (RAD-Sex) that makes use of next-generation sequencing technology to identify genetic markers that define sex-specific segments of the male or female genome.

A PCR assay detects a male-specific duplicated copy of Anti-Müllerian hormone (amh) in the lingcod (Ophiodon elongatus)

BACKGROUND

Anti-Müllerian hormone (amh) or Müllerian-inhibiting substance (mis) is a member of the transforming growth factor-β family of hormones. This gene plays a key role in vertebrate male sex-determination by inhibiting the development of the Müllerian ducts, and has been shown to be the master sex-determinant in the Patagonian pejerrey.

RESULTS

In the lingcod, Ophiodon elongatus, both males and females share one copy of amh, however we have identified a second duplicate copy that appears solely in the male individuals. We have developed a PCR-based assay targeting the TGF-β domain of amh that provides a simple method with which to sex lingcod from a small amount of tissue. An analysis across 57 individuals gave a 100% success rate in identifying the phenotypic sex.

CONCLUSIONS

We present a simple method to sex lingcod through non-lethal tissue sampling. A third, independent, male-specific duplication of amh in a teleost fish has been identified in the lingcod.

Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish

The killifish Nothobranchius furzeri is the shortest-lived vertebrate that can be bred in the laboratory. Its rapid growth, early sexual maturation, fast aging, and arrested embryonic development (diapause) make it an attractive model organism in biomedical research. Here, we report a draft sequence of its genome that allowed us to uncover an intra-species Y chromosome polymorphism representing-in real time-different stages of sex chromosome formation that display features of early mammalian XY evolution "in action." Our data suggest that gdf6Y, encoding a TGF-β family growth factor, is the master sex-determining gene in N. furzeri. Moreover, we observed genomic clustering of aging-related genes, identified genes under positive selection, and revealed significant similarities of gene expression profiles between diapause and aging, particularly for genes controlling cell cycle and translation. The annotated genome sequence is provided as an online resource (http://www.nothobranchius.info/NFINgb).