A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus).

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 ~ 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.

Comparative expression patterns of Sox2 and Sox19 genes in the forebrain of developing and adult turbot (Scophthalmus maximus).

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.

Cytogenomic analysis of several repetitive DNA elements in turbot (Scophthalmus maximus).

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.

Whole genome sequencing of turbot (Scophthalmus maximus; Pleuronectiformes): a fish adapted to demersal life.

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.

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

BACKGROUND: Controlling sex ratios is essential for the aquaculture industry, especially in those species with sex dimorphism for relevant productive traits, hence the importance of knowing how the sexual phenotype is established in fish. Turbot, a very important fish for the aquaculture industry in Europe, shows one of the largest sexual growth dimorphisms amongst marine cultured species, being all-female stocks a desirable goal for the industry. Although important knowledge has been achieved on the genetic basis of sex determination (SD) in this species, the master SD gene remains unknown and precise information on gene expression at the critical stage of sex differentiation is lacking. In the present work, we examined the expression profiles of 29 relevant genes related to sex differentiation, from 60 up to 135 days post fertilization (dpf), when gonads are differentiating. We also considered the influence of three temperature regimes on sex differentiation. RESULTS: The first sex-related differences in molecular markers could be observed at 90 days post fertilization (dpf) and so we have called that time the onset of sex differentiation. Three genes were the first to show differential expression between males and females and also allowed us to sex turbot accurately at the onset of sex differentiation (90 dpf): cyp19a1a, amh and vasa. The expression of genes related to primordial germ cells (vasa, gsdf, tdrd1) started to increase between 75-90 dpf and vasa and tdrd1 later presented higher expression in females (90-105 dpf). Two genes placed on the SD region of turbot (sox2, fxr1) did not show any expression pattern suggestive of a sex determining function. We also detected changes in the expression levels of several genes (ctnnb1, cyp11a, dmrt2 or sox6) depending on culture temperature. CONCLUSION: Our results enabled us to identify the first sex-associated genetic cues (cyp19a1a, vasa and amh) at the initial stages of gonad development in turbot (90 dpf) and to accurately sex turbot at this age, establishing the correspondence between gene expression profiles and histological sex. Furthermore, we profiled several genes involved in sex differentiation and found specific temperature effects on their expression.

Genetic architecture of sex determination in fish: applications to sex ratio control in aquaculture.

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.

Analysis of qPCR reference gene stability determination methods and a practical approach for efficiency calculation on a turbot (Scophthalmus maximus) gonad dataset.

BACKGROUND: Gene expression analysis by reverse transcription quantitative PCR (qPCR) is the most widely used method for analyzing the expression of a moderate number of genes and also for the validation of microarray results. Several issues are crucial for a successful qPCR study, particularly the selection of internal reference genes for normalization and efficiency determination. There is no agreement on which method is the best to detect the most stable genes neither on how to perform efficiency determination. In this study we offer a comprehensive evaluation of the characteristics of reference gene selection methods and how to decide which one is more reliable when they show discordant outcomes. Also, we analyze the current efficiency calculation controversy. Our dataset is composed by gonad samples of turbot at different development times reared at different temperatures. Turbot (Scophthalmus maximus) is a relevant marine aquaculture European species with increasing production in the incoming years. Since females largely outgrow males, identification of genes related to sex determination, gonad development and reproductive behavior, and analysis of their expression profiles are of primary importance for turbot industry. RESULTS: We analyzed gene stability of six reference genes: RPS4, RPL17, GAPDH, ACTB, UBQ and B2M using the comparative delta-CT method, Bestkeeper, NormFinder and GeNorm approaches in gonad samples of turbot. Supported by descriptive statistics, we found NormFinder to be the best method, while on the other side, GeNorm results proved to be unreliable. According to our analysis, UBQ and RPS4 were the most stable genes, while B2M was the least stable gene. We also analyzed the efficiency calculation softwares LinRegPCR, LREanalyzer, DART and PCR-Miner and we recommend LinRegPCR for research purposes since it does not systematically overestimate efficiency. CONCLUSION: Our results indicate that NormFinder and LinRegPCR are the best approaches for reference gene selection and efficiency determination, respectively. We also recommend the use of UBQ and RPS4 for normalization of gonad development samples in turbot.

Fine mapping and evolution of the major sex determining region in turbot (Scophthalmus maximus).

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.

A sex-associated sequence identified by RAPD screening in gynogenetic individuals of turbot (Scophthalmus maximus).

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.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones.

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.

Mapping of DNA sex-specific markers and genes related to sex differentiation in turbot (Scophthalmus maximus).

Production of all-female populations in turbot can increase farmer's benefits since sexual dimorphism in growth in this species is among the highest within marine fish, turbot females reaching commercial size 3-6 months earlier than males. Puberty in males occurs earlier than in females, which additionally slows their growth. Thus, elucidating the mechanisms of sex determination and gonad differentiation is a relevant goal for turbot production. A ZZ/ZW sex determination mechanism has been suggested for this species, and four sex-related quantitative trait loci (QTL) were detected, the major one located in linkage group (LG) 5 and the three minor ones in LG6, LG8, and LG21. In the present work, we carried out a linkage analysis for several sex-related markers: (1) three anonymous sex-associated RAPD and (2) several candidate genes related to sex determination and gonad differentiation in other species (Sox3, Sox6, Sox8, Sox9, Sox17, Sox19, Amh, Dmrta2, Cyp19a, Cyp19b). We focused our attention on their co-localization with the major and minor sex-related QTL trying to approach to the master sex-determining gene of this species. Previously described growth-related QTL were also considered since the association observed between growth and sex determination in fish. Amh, Dmrta2, and one RAPD were located in LG5, while Sox9 and Sox17 (LG21), Cyp19b (LG6), and a second RAPD (LG8) co-mapped with suggestive sex-related QTL, thus supporting further analyses on these genes to elucidate the genetic basis of this relevant trait for turbot farming.

Comparative expression analysis in mature gonads, liver and brain of turbot (Scophthalmus maximus) by cDNA-AFLPS.

Turbot is one of the most important farmed fish in Europe. This species exhibits a considerable sexual dimorphism in growth and sexual maturity that makes the all-female production recommended for turbot farming. Our knowledge about the genetic basis of sex determination and the molecular regulation of gonad differentiation in this species is still limited. Our goal was to identify and compare gene expression and functions between testes and ovaries in adults in order to ascertain the relationship between the genes that could be involved in the gonad differentiation or related to the sex determination system. The identification of differentially expressed sex related genes is an initial step towards understanding the molecular mechanisms of gonad differentiation. For this, we carried out a transcriptome analysis based on cDNA-AFLP technique which allowed us to obtain an initial frame on sex-specific gene expression that will facilitate further analysis especially along the critical gonad differentiating period. With the aim of widening the study on sex-biased gene expression we reproduced the same experiments in two somatic tissues: liver and brain. We have selected the liver because it is the most analyzed one regarding sexual dimorphic gene expression and due to its importance in steroid hormones metabolism and the brain because the functional relationship between brain and gonad is documented. We found slight but important differences between sexes which deserve further investigation.

Identification of the major sex-determining region of turbot (Scophthalmus maximus).

Sex determination in fish is a labile character in evolutionary terms. The sex-determining (SD) master gene can differ even between closely related fish species. This group is an interesting model for studying the evolution of the SD region and the gonadal differentiation pathway. The turbot (Scophthalmus maximus) is a flatfish of great commercial value, where a strong sexual dimorphism exists for growth rate. Following a QTL and marker association approach in five families and a natural population, we identified the main SD region of turbot at the proximal end of linkage group (LG) 5, close to the SmaUSC-E30 marker. The refined map of this region suggested that this marker would be 2.6 cM and 1.4 Mb from the putative SD gene. This region appeared mostly undifferentiated between males and females, and no relevant recombination frequency differences were detected between sexes. Comparative genomics of LG5 marker sequences against five model species showed no similarity of this chromosome to the sex chromosomes of medaka, stickleback, and fugu, but suggested a similarity to a sex-associated QTL from Oreochromis spp. The segregation analysis of the closest markers to the SD region demonstrated a ZW/ZZ model of sex determination in turbot. A small proportion of families did not fit perfectly with this model, which suggests that other minor genetic and/or environmental factors are involved in sex determination in this species.

Characterization of a new HpaI centromeric satellite DNA in Salmo salar.

A highly repeated HpaI DNA family was revealed in Atlantic salmon (Salmo salar) and analyzed by Southern blotting and fluorescence in situ hybridization (FISH). In this report, we describe the nucleotide sequence, genomic structure and chromosomal localization of this HpaI repeat. This novel satellite appeared tandemly arrayed and located at centromeric areas of three acrocentric chromosome pairs as evidenced by FISH. The sequence was characterized by a high AT content (63%), a short consensus motif (A/T)(G/C)AAA(T/C) similar to other centromeric satellites motifs, and by short AT enriched stretches. The presence of this sequence in other salmonid species was also tested by Southern blot hybridization and used to analyze its evolution within this group.

Recombination analysis of the human minisatellite MsH42 suggests the existence of two distinct pathways for initiation and resolution of recombination at MsH42 in rat testes nuclear extracts.

We have previously described a GC-rich human minisatellite, termed MsH42, which exists in two allelic forms, long and short. Here, we have identified a third allele of medium length and localized the MsH42 locus in the chromosome 15q25.1 inside an intron belonging to a gene of unknown function. The recombinogenic potential of the three alleles was assayed in vitro incubating pBR322-based constructs containing two copies of the minisatellite MsH42 with its flanking sequences, in the presence of rat testes nuclear extracts. This assay system was configured to monitor only reciprocal exchange type events and not gene conversion. All MsH42 allelic sequences enhanced intramolecular homologous recombination promoting high rates (approximately 76%) of equal crossover, the long allele showing the highest recombinogenic activity. Removal of the MsH42 long allele flanking sequences, which are identical in the three alleles, provoked a decrease in the enhancement of recombination and in the frequency of equal crossovers, suggesting that these sequences are important for the recombinogenic activity and for the correct pairing between homologous sequences. The occurrence of some complex recombination events within the minisatellite MsH42 suggests the existence of processes related to polymerase slippage and unwinding with reinvasion during the repair synthesis. Our findings point toward the existence of two distinct biochemical pathways for initiation and resolution of recombination at the minisatellite MsH42. Finally, the in vitro recombination system employed in this study could provide an approach to dissect processes of repetitive DNA instability and recombination.