Genes involved in sexual development have been identified in many teleost fishes, including the turbot, a flatfish of high commercial value in aquaculture. In this species, a major sex determination (SD) Quantitative Trait Locus has been identified, and Sox2 has been proposed as the SD gene. Although RT-PCR, qPCR, and transcriptomic analyses have been performed on turbot gonads, histological studies remain limited. Here, we examined the tissue-specific expression patterns of several sex-related genes in turbot using in situ hybridization on paraffin sections (SISH) of differentiated ovaries and testes of juvenile and adult specimens. Vasa, Foxl2, Amh, Sox2 and Sox17 transcripts were detected in both male and female gonads, whereas Cyp19a1a and Sox19 expression was restricted to the ovary; these results support previous gene expression analyses and suggest a role for these genes in gonadal development and reproductive functions in this species. SISH provides complementary information to molecular analyses by identifying the specific cell types expressing the sex-related genes analyzed in gonadal tissue, thereby offering a more comprehensive understanding of gonadal differentiation in turbot. Comparison with results reported for the gonads of other teleost species revealed similar tissue-specific gene expression patterns during sexual development.

It is known that the insulin-like growth factor (Igf) system is involved in gonadal development and, in turn, reproduction in vertebrates. In this study, the circulating levels of Igf-1 and follicle-stimulating hormone (Fsh) were assessed as physiological signals related to gonadal maturation stage in juvenile male European sea bass at 1 year of age. The seasonal profiles of Igf-1 and Fsh were analyzed in non-precocious (NP) and precocious (PR) fish during the early onset of puberty. Subsequently, gonadal expression patterns of the three igf genes and their two receptors, igf-1r and igf-2r, were studied during their first year of life. The analysis was continued to evaluate the mRNA expression profiles of three genes (sox17, scf, nf2) mediating spermatogenesis progression. Collectively, these results demonstrated that PR males have higher circulating levels of Igf-1 as compared to NP fish before and during the putative gonadal development period (PGD). Plasma Fsh levels were significantly higher in PR males as compared to NP fish before the start of the PGD period, although this was followed by a decrease in Fsh levels and a subsequent recovery, coinciding with the reproductive season (February). Results showed that the total variance associated with the appearance of precocity may be explained by the contribution of factors such as body size and plasma levels of Igf-1, Fsh and sex steroids in male sea bass. Expression of igfs was high during stage I, and this was significantly correlated with the expression of genes related to cell cycle progression, including sox17, scf and nf2. Finally, the decreased expression of all these gonadal factors in precocious animals suggests their potential involvement in regulating the proliferative growth of spermatogonia and supporting cells during the early stages of spermatogenesis, thus evoking the onset of puberty in male sea bass.

Accumulating evidence suggests that the insulin-like growth factor (Igf) system plays a crucial role in gonadal development and, in turn, reproduction. The aim of this study was to assess the circulating levels of Igf-1 and follicle-stimulating hormone (Fsh) as indicators of the gonadal maturation stage in juvenile male European sea bass at 1 year of age. The seasonal profiles of Igf-1 and Fsh were analyzed in non-precocious (NP) and precocious (PR) fish during the early onset of puberty. Subsequently, gonadal expression patterns of igf-1, -2 and -3 and their receptors, igf-1r and igf-2r, were studied during the first reproductive cycle of male sea bass. The analysis was extended to determine the mRNA expression profiles of three genes (sox17, scf, nf2) mediating spermatogenesis progression. Collectively, these results demonstrated that PR males have higher circulating levels of Igf-1 as compared to NP fish before and during the putative gonadal development period (PGD). Circulating Fsh levels were significantly higher in PR males as compared to NP fish before the start of the PGD period, although this was followed by a decrease in Fsh levels and a subsequent recovery, coinciding with the reproductive season (February). Body size and circulating plasma levels (Igf-1, Fsh, sex steroids) significantly contributed to explaining the total variance associated with early onset of puberty in male sea bass. Expression of igfs was high during stage I, and this was significantly correlated with the expression of genes related to cell cycle progression, including sox17, scf and nf2. Finally, the decreased expression of all these gonadal factors in precocious animals suggests their potential involvement in regulating the proliferative growth of spermatogonia and supporting cells during the early stages of spermatogenesis, thus evoking the onset of puberty in male sea bass.

Background: Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches. Results: The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a similar to 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD. Conclusions: The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

ABSTRACTA major challenge in evolutionary biology is to find an explanation for the variation in sex-determining (SD) systems across taxa and to understand the mechanisms driving sex chromosome differentiation. We studied the turbot, holding a ZW/ZZ SD system and no sex chromosome heteromorphism, by combining classical genetics and genomics approaches to disentangle the genetic architecture of this trait. RAD-Seq was used to genotype 18,214 SNPs on 1,135 fish from 36 families and a genome wide association study (GWAS) identified a ∼ 6 Mb region on LG5 associated with sex (P < 0.05). The most significant associated markers were located close to sox2, dnajc19 and fxr1 genes. A segregation analysis enabled narrowing down the associated region and evidenced recombination suppression in a region overlapping the candidate genes. A Nanopore/Illumina assembly of the SD region using ZZ and WW individuals identified a single SNP fully associated with Z and W chromosomes. RNA-seq from 5-90 day-old fish detected the expression along the gonad differentiation period of a short non-coding splicing variant (ncRNA) included in a vertebrate-conserved long non-coding RNA overlapping sox2. qPCR showed that sox2 was the only differentially expressed gene between males and females at 50-55 days post fertilization, just prior the beginning of gonad differentiation. More refined information on the involvement of secondary genetic and environmental factors and their interactions on SD was gathered after the analysis of a broad sample of families. Our results confirm the complex nature of SD in turbot and support sox2 as its main driver.

A major challenge in evolutionary biology is to find an explanation for the variation in sex-determining (SD) systems across taxa and to understand the mechanisms driving sex chromosome differentiation. We studied the turbot, holding a ZW/ZZ SD system and no sex chromosome heteromorphism, by combining classical genetics and genomics approaches to disentangle the genetic architecture of this trait. RAD-Seq was used to genotype 18,214 SNPs on 1,135 fish from 36 families and a genome wide association study (GWAS) identified a ∼ 6 Mb region on LG5 associated with sex (P < 0.05). The most significant associated markers were located close to , and genes. A segregation analysis enabled narrowing down the associated region and evidenced recombination suppression in a region overlapping the candidate genes. A Nanopore/Illumina assembly of the SD region using ZZ and WW individuals identified a single SNP fully associated with Z and W chromosomes. RNA-seq from 5-90 day-old fish detected the expression along the gonad differentiation period of a short non-coding splicing variant (ncRNA) included in a vertebrate-conserved long non-coding RNA overlapping . qPCR showed that was the only differentially expressed gene between males and females at 50-55 days post fertilization, just prior the beginning of gonad differentiation. More refined information on the involvement of secondary genetic and environmental factors and their interactions on SD was gathered after the analysis of a broad sample of families. Our results confirm the complex nature of SD in turbot and support as its main driver.

Repetitive DNA plays a fundamental role in the organization, size and evolution of eukaryotic genomes. The sequencing of the turbot revealed a small and compact genome, as in all flatfish studied to date. The assembly of repetitive regions is still incomplete because it is difficult to correctly identify their position, number and array. The combination of classical cytogenetic techniques along with high quality sequencing is essential to increase the knowledge of the structure and composition of these sequences and, thus, of the structure and function of the whole genome. In this work, the in silico analysis of H1 histone, 5S rDNA, telomeric and Rex repetitive sequences, was compared to their chromosomal mapping by fluorescent in situ hybridization (FISH), providing a more comprehensive picture of these elements in the turbot genome. FISH assays confirmed the location of H1 in LG8; 5S rDNA in LG4 and LG6; telomeric sequences at the end of all chromosomes whereas Rex elements were dispersed along most chromosomes. The discrepancies found between both approaches could be related to the sequencing methodology applied in this species and also to the resolution limitations of the FISH technique. Turbot cytogenomic analyses have proven to add new chromosomal landmarks in the karyotype of this species, representing a powerful tool to investigate targeted genomic sequences or regions in the genetic and physical maps of this species.

Twelfth International Symposium on Genetics in Aquaculture 2015 (ISGA XII), 21-27 June 2015, Santiago de Compostela, Spain.-- 1 page

The turbot, Scophthalmus maximus, belongs to the flatfishes (order Pleuronectiformes), which display substantial asymmetry of the olfactory organs and forebrain. Sox genes code for SRY-related HMG domain-bearing transcription factors involved in various developmental processes. Group B1 Sox genes as Sox2 and Sox19 appear to play major roles in neural development. Here, we characterized by in situ hybridization the developmental expression of Sox2 and Sox19 genes in metamorphic and postmetamorphic specimens and young adults of both sexes. Expression of S. maximus Sox2 (Sm-Sox2) and Sm-Sox19 mRNAs was detected in ependymal cells of different regions of the telencephalon, preoptic region, hypothalamus, and thalamus at all stages investigated. Sm-Sox2 expression but not Sm-Sox19 occurred in neurons located in particular regions such as the dorsal nucleus of the ventral telencephalon, the medial zone of the dorsal telencephalon, preoptic area and hypothalamus. Although Sm-Sox2 and Sm-Sox19 are expressed differentially in gonads, no sex differences in their expression were observed between male and female forebrains. We also investigated the topographical relation between Sox expression and cell proliferation using series double immunostained for a radial glial marker (BLBP) and cell proliferation marker (PCNA). Sm-Sox2 and Sm-Sox19 were strongly expressed in ependymal cells located in neurogenic niches revealed by the BLBP and PCNA immunostaining. Comparison with other teleosts indicates similar expression of Sox2 and Sox19 in the telencephalon, supporting conserved roles for both genes in teleost brains.

The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot.

International Symposium on Genetics in Aquaculture XII (ISGA XII), 21-27 June 2015, Santiago de Compostela, Spain.-- 1 page

The turbot is a flatfish with a ZW/ZZ sex determination system but with a still unknown sex determining gene(s), and with a marked sexual growth dimorphism in favor of females. To better understand sexual development in turbot we sampled young turbot encompassing the whole process of gonadal differentiation and conducted a comprehensive transcriptomic study on its sex differentiation using a validated custom oligomicroarray. Also, the expression profiles of 18 canonical reproduction-related genes were studied along gonad development. The expression levels of gonadal aromatase cyp19a1a alone at three months of age allowed the accurate and early identification of sex before the first signs of histological differentiation. A total of 56 differentially expressed genes (DEG) that had not previously been related to sex differentiation in fish were identified within the first three months of age, of which 44 were associated with ovarian differentiation (e.g., cd98, gpd1 and cry2), and 12 with testicular differentiation (e.g., ace, capn8 and nxph1). To identify putative sex determining genes, ∼4.000 DEG in juvenile gonads were mapped and their positions compared with that of previously identified sex- and growth-related quantitative trait loci (QTL). Although no genes mapped to the previously identified sex-related QTLs, two genes (foxl2 and 17βhsd) of the canonical reproduction-related genes mapped to growth-QTLs in linkage group (LG) 15 and LG6, respectively, suggesting that these genes are related to the growth dimorphism in this species.

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.

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.

Controlling the sex ratio is essential in finfish farming. A balanced sex ratio is usually good for broodstock management, since it enables to develop appropriate breeding schemes. However, in some species the production of monosex populations is desirable because the existence of sexual dimorphism, primarily in growth or first time of sexual maturation, but also in color or shape, can render one sex more valuable. The knowledge of the genetic architecture of sex determination (SD) is convenient for controlling sex ratio and for the implementation of breeding programs. Unlike mammals and birds, which show highly conserved master genes that control a conserved genetic network responsible for gonad differentiation (GD), a huge diversity of SD mechanisms has been reported in fish. Despite theory predictions, more than one gene is in many cases involved in fish SD and genetic differences have been observed in the GD network. Environmental factors also play a relevant role and epigenetic mechanisms are becoming increasingly recognized for the establishment and maintenance of the GD pathways. Although major genetic factors are frequently involved in fish SD, these observations strongly suggest that SD in this group resembles a complex trait. Accordingly, the application of quantitative genetics combined with genomic tools is desirable to address its study and in fact, when applied, it has frequently demonstrated a multigene trait interacting with environmental factors in model and cultured fish species. This scenario has notable implications for aquaculture and, depending upon the species, from chromosome manipulation or environmental control techniques up to classical selection or marker assisted selection programs, are being applied. In this review, we selected four relevant species or fish groups to illustrate this diversity and hence the technologies that can be used by the industry for the control of sex ratio: turbot and European sea bass, two reference species of the European aquaculture, and salmonids and tilapia, representing the fish for which there are well established breeding programs.

Fish sex determination (SD) systems are varied, suggesting evolutionary changes including either multiple evolution origins of genetic SD from nongenetic systems (such as environmental SD) and/or turnover events replacing one genetic system by another. When genetic SD is found, cytological differentiation between the two members of the sex chromosome pair is often minor or undetectable. The turbot (Scophthalmus maximus), a valuable commercial flatfish, has a ZZ/ZW system and a major SD region on linkage group 5 (LG5), but there are also other minor genetic and environmental influences. We here report refined mapping of the turbot SD region, supported by comparative mapping with model fish species, to identify the turbot master SD gene. Six genes were located to the SD region, two of them associated with gonad development (sox2 and dnajc19). All showed a high association with sex within families (P = 0), but not at the population level, so they are probably partially sex-linked genes, but not SD gene itself. Analysis of crossovers in LG5 using two families confirmed a ZZ/ZW system in turbot and suggested a revised map position for the master gene. Genetic diversity and differentiation for 25 LG5 genetic markers showed no differences between males and females sampled from a wild population, suggesting a recent origin of the SD region in turbot. We also analyzed associations with markers of the most relevant sex-related linkage groups in brill (S. rhombus), a closely related species to turbot; the data suggest that an ancient XX/XY system in brill changed to a ZZ/ZW mechanism in turbot.

Understanding the genetic basis of sex determination mechanisms is essential for improving the productivity of farmed aquaculture fish species like turbot (Scophthalmus maximus). In culture conditions turbot males grow slower than females starting from eight months post-hatch, and this differential growth rate is maintained until sexual maturation is reached, being mature females almost twice as big as males of the same age. The goal of this study was to identify sex-specific DNA markers in turbot using comparative random amplified polymorphism DNA (RAPD) profiles in males and females to get new insights of the genetic architecture related to sex determination. In order to do this, we analyzed 540 commercial 10-mer RAPD primers in male and female pools of a gynogenetic family because of its higher inbreeding, which facilitates the detection of associations across the genome. Two sex-linked RAPD markers were identified in the female pool and one in the male pool. After the analysis of the three markers on individual samples of each pool and also in unrelated individuals, only one RAPD showed significant association with females. This marker was isolated, cloned and sequenced, containing two sequences, a microsatellite (SEX01) and a minisatellite (SEX02), which were mapped in the turbot reference map. From this map position, through a comparative mapping approach, we identified Foxl2, a relevant gene related to initial steps of sex differentiation, and Wnt4, a gene related with ovarian development, close to the microsatellite and minisatellite markers, respectively. The position of Foxl2 and Wnt4 was confirmed by linkage mapping in the reference turbot map.