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.

No abstract

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.

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.

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.