A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis.

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis.

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5-9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

Cytogenomics Unveil Possible Transposable Elements Driving Rearrangements in Chromosomes 2 and 4 of Solea senegalensis.

Cytogenomics, the integration of cytogenetic and genomic data, has been used here to reconstruct the evolution of chromosomes 2 and 4 of Solea senegalensis. S. senegalensis is a flat fish with a karyotype comprising 2n = 42 chromosomes: 6 metacentric + 4 submetacentric + 8 subtelocentric + 24 telocentric. The Fluorescence in situ Hybridization with Bacterial Artificial Chromosomes (FISH-BAC) technique was applied to locate BACs in these chromosomes (11 and 10 BACs in chromosomes 2 and 4, respectively) and to generate integrated maps. Synteny analysis, taking eight reference fish species (Cynoglossus semilaevis, Scophthalmus maximus, Sparus aurata, Gasterosteus aculeatus, Xiphophorus maculatus, Oryzias latipes, Danio rerio, and Lepisosteus oculatus) for comparison, showed that the BACs of these two chromosomes of S. senegalensis were mainly distributed in two principal chromosomes in the reference species. Transposable Elements (TE) analysis showed significant differences between the two chromosomes, in terms of number of loci per Mb and coverage, and the class of TE (I or II) present. Analysis of TE divergence in chromosomes 2 and 4 compared to their syntenic regions in four reference fish species (C. semilaevis, S. maximus, O. latipes, and D. rerio) revealed differences in their age of activity compared with those species but less notable differences between the two chromosomes. Differences were also observed in peaks of divergence and coverage of TE families for all reference species even in those close to S. senegalensis, like S. maximus and C. semilaevis. Considered together, chromosomes 2 and 4 have evolved by Robertsonian fusions, pericentric inversions, and other chromosomal rearrangements mediated by TEs.

High-Density Linkage Maps Based on Genotyping-by-Sequencing (GBS) Confirm a Chromosome-Level Genome Assembly and Reveal Variation in Recombination Rate for the Pacific Oyster Crassostrea gigas.

Studies of linkage and linkage mapping have advanced genetic and biological knowledge for over 100 years. In addition to their growing role, today, in mapping phenotypes to genotypes, dense linkage maps can help to validate genome assemblies. Previously, we showed that 40% of scaffolds in the first genome assembly for the Pacific oyster Crassostrea gigas were chimeric, containing single nucleotide polymorphisms (SNPs) mapping to different linkage groups. Here, we merge 14 linkage maps constructed of SNPs generated from genotyping-by-sequencing (GBS) methods with five, previously constructed linkage maps, to create a compendium of nearly 69 thousand SNPs mapped with high confidence. We use this compendium to assess a recently available, chromosome-level assembly of the C. gigas genome, mapping SNPs in 275 of 301 contigs and comparing the ordering of these contigs, by linkage, to their assembly by Hi-C sequencing methods. We find that, while 26% of contigs contain chimeric blocks of SNPs, i.e., adjacent SNPs mapping to different linkage groups than the majority of SNPs in their contig, these apparent misassemblies amount to only 0.08% of the genome sequence. Furthermore, nearly 90% of 275 contigs mapped by linkage and sequencing are assembled identically; inconsistencies between the two assemblies for the remaining 10% of contigs appear to result from insufficient linkage information. Thus, our compilation of linkage maps strongly supports this chromosome-level assembly of the oyster genome. Finally, we use this assembly to estimate, for the first time in a Lophotrochozoan, genome-wide recombination rates and causes of variation in this fundamental process.

The genomic structure of the highly-conserved dmrt1 gene in Solea senegalensis (Kaup, 1868) shows an unexpected intragenic duplication.

Knowing the factors responsible for sex determination in a species has significant theoretical and practical implications; the dmrt1 gene (Doublesex and Mab-3 (DM)-related Transcription factor 1) plays this role in diverse animal species. Solea senegalensis is a commercially important flat fish in which females grow 30% faster than males. It has 2n = 42 chromosomes and an XX / XY chromosome system for sex determination, without heteromorph chromosomes but with sex proto-chromosome. In the present study, we are providing the genomic structure and nucleotide sequence of dmrt1 gene obtained from cDNA from male and female adult gonads. A cDNA of 2027 containing an open-reading frame (ORF) of 1206 bp and encoding a 402 aa protein it is described for dmrt1 gene of S. senegalensis. Multiple mRNA isoforms indicating a high variable system of alternative splicing in the expression of dmrt1 of the sole in gonads were studied. None isoforms could be related to sex of individuals. The genomic structure of the dmrt1 of S. senegalensis showed a gene of 31400 bp composed of 7 exons and 6 introns. It contains an unexpected duplication of more than 10399 bp, involving part of the exon I, exons II and III and a SINE element found in the sequence that it is proposed as responsible for the duplication. A mature miRNA of 21 bp in length was localized at 336 bp from exon V. Protein-protein interacting networks of the dmrt1 gene showed matches with dmrt1 protein from Cynoglossus semilaevis and a protein interaction network with 11 nodes (dmrt1 plus 10 other proteins). The phylogenetic relationship of the dmrt1 gene in S. senegalensis is consistent with the evolutionary position of its species. The molecular characterization of this gene will enhance its functional analysis and the understanding of sex differentiation in Solea senegalensis and other flatfish.

A Comprehensive Integrated Genetic Map of the Complete Karyotype of Solea senegalensis (Kaup 1858).

Solea senegalensis aquaculture production has experienced a great increase in the last decade and, consequently, the genome knowledge of the species is gaining attention. In this sense, obtaining a high-density genome mapping of the species could offer clues to the aquaculture improvement in those aspects not resolved so far. In the present article, a review and new processed data have allowed to obtain a high-density BAC-based cytogenetic map of S. senegalensis beside the analysis of the sequences of such BAC clones to achieve integrative data. A total of 93 BAC clones were used to localize the chromosome complement of the species and 588 genes were annotated, thus almost reaching the 2.5% of the S. senegalensis genome sequences. As a result, important data about its genome organization and evolution were obtained, such as the lesser gene density of the large metacentric pair compared with the other metacentric chromosomes, which supports the theory of a sex proto-chromosome pair. In addition, chromosomes with a high number of linked genes that are conserved, even in distant species, were detected. This kind of result widens the knowledge of this species' chromosome dynamics and evolution.

Evolution of the Proto Sex-Chromosome in Solea senegalensis.

Solea senegalensis is a flatfish belonging to the Soleidae family within the Pleuronectiformes order. It has a karyotype of 2n = 42 (FN = 60; 6M + 4 SM + 8 St + 24 T) and a XX/XY system. The first pair of metacentric chromosomes has been proposed as a proto sex-chromosome originated by a Robertsonian fusion between acrocentric chromosomes. In order to elucidate a possible evolutionary origin of this chromosome 1, studies of genomic synteny were carried out with eight fish species. A total of 88 genes annotated within of 14 BACs located in the chromosome 1 of S. senegalensis were used to elaborate syntenic maps. Six BACs (BAC5K5, BAC52C17, BAC53B20, BAC84K7, BAC56H24, and BAC48P7) were distributed in, at least, 5 chromosomes in the species studied, and a group of four genes from BAC53B20 (grsf1, rufy3, slc4a4 and npffr2) and genes from BAC48K7 (dmrt2, dmrt3, dmrt1, c9orf117, kank1 and fbp1) formed a conserved cluster in all species. The analysis of repetitive sequences showed that the number of retroelements and simple repeat per BAC showed its highest value in the subcentromeric region where 53B20, 16E16 and 48K7 BACs were localized. This region contains all the dmrt genes, which are associated with sex determination in some species. In addition, the presence of a satellite "chromosome Y" (motif length: 860 bp) was detected in this region. These findings allowed to trace an evolutionary trend for the large metacentric chromosome of S. senegalensis, throughout different rearrangements, which could be at an initial phase of differentiation as sex chromosome.

Assessing the geographic scale of genetic population management with microsatellites and introns in the clam Ruditapes decussatus.

The clam Ruditapes decussatus is commercially important in southwestern Europe, suffering from population decline and hybridization with exotic Manila clam (R. philippinarum). Previous studies with intronic markers showed a genetic subdivision of the species in three races (Atlantic, West Mediterranean, and Adriatic-Aegean). However, detailed population genetic studies to help management of the main production areas in the southwest of Europe are missing. We have analyzed eight Atlantic and two Mediterranean populations from the Spanish coasts using 14 microsatellites and six intronic markers. Microsatellites confirmed the Atlantic and West Mediterranean races detected with introns and showed that genetic variability was higher in Mediterranean than in Atlantic populations. Both marker types showed that genetic differentiation of Atlantic populations was low and indicated that populations could be managed at the regional level in the case of Cantabrian and Gulf of Cadiz areas, but not in the case of Rias Baixas and the Mediterranean. This study shows the interest of including different types of markers in studies of genetic population structure of marine organisms.

Characterization of nineteen microsatellite markers and development of multiplex PCRs for the wedge clam Donax trunculus (Mollusca: Bivalvia).

The wedge clam Donax trunculus is an Atlantic-Mediterranean warm-temperate species found from Senegal to the northern coast of France, including the Mediterranean and Black Sea. It is commercially exploited in several European countries and constitutes an important fishing resource due to its high economical value. To contribute to its conservation and management, nineteen microsatellite markers were isolated from two enriched genomic libraries. These loci were characterized in 30 clams from a single population from northwest Spain. The number of alleles per locus ranged from 2 to 17 and observed and expected heterozygosity varied from zero to 0.714 and from 0.078 to 0.950, respectively. Linkage disequilibrium was not detected and nine loci were in agreement with Hardy-Weinberg equilibrium. Fifteen polymorphic markers were arranged into three multiplex PCR sets to reduce both time and cost of microsatellite genotyping. This is the first time that polymorphic microsatellite markers have been reported for D. trunculus. These new markers provide a valuable resource for future population genetics studies and management and culture of this species.

Development and multiplex PCR amplification of microsatellite markers in the commercial clam Venerupis rhomboides (Mollusca: Bivalvia).

Venerupis rhomboides is a commercial clam whose production could be enhanced through effective management of natural and hatchery stocks. This study provides the first panel of microsatellite markers for the exploitation of this species according to genetic criteria. A total of 22 polymorphic microsatellite loci were isolated and characterized from two genomic libraries enriched for different motifs. The number of alleles per locus ranged from 2 to 14 in a sample of 20 clams from Spain, and the observed and expected heterozygosity from 0 to 0.95 and 0.05-0.901, respectively. Sixteen loci were in agreement with Hardy-Weinberg equilibrium after sequential Bonferroni correction and linkage disequilibrium between loci pairs was not detected. To reduce the cost of the genotyping process, tri- and pentaplex PCRs, amplifying a total of 13 microsatellites loci were optimized. The microsatellites developed here represent the first nuclear markers described in V. rhomboides and will be useful tools for genetic studies involving assessment of genetic variation and population structure of natural and cultivated populations, assignment testing, construction of genetic linkage maps and dissection of production traits.

Identification of Ensis siliqua samples and establishment of the catch area using a species-specific microsatellite marker.

European Council Regulation 104/2000 states that fishery products must be labeled to indicate commercial designation of species, the production method, and the catch area. Therefore, traceability of seafood implies knowledge of the species offered to retail and their origin. Ensis siliqua is a bivalve intensively fished in Europe and sold in fresh and canned forms. Although several published methods clearly differentiate Ensis genus species, none of those assess the origin of the commercial samples. In the present study, a microsatellite marker (Esi-UDC3055F) was developed to establish the catch area of E. siliqua samples. Amplification yielded a fragment of 275 or 302 base pairs, depending on whether they were Iberian or Irish populations. The usefulness of this method was also assessed in commercial samples. The results of this study provide a reliable methodology for the identification of catch area in European E. siliqua commercial samples. The coupling of this methodology with existing techniques for razor clam species identification provides a powerful tool for traceability and labeling enforcement.

Isolation of microsatellite markers and analysis of genetic diversity among east Atlantic populations of the sword razor shell Ensis siliqua: a tool for population management.

The sword razor shell Ensis siliqua (Linnaeus, 1758) is found mainly from Norway to the Atlantic coast of the Iberian Peninsula. It is intensively caught in Europe, being highly appreciated as seafood. To help in its conservation and management, five microsatellite markers were isolated and genetic variation was analyzed in samples from Ireland, Spain, and Portugal. The highly significant differentiation (θ = 0.287, P < 0.001) observed was mainly due to differences between samples from Irish and Iberian Peninsula localities, except Aveiro (its sample resembled the Irish samples, and it may be predominantly self-recruiting). These groups of samples showed significant differences in allelic richness that could be related to harvesting intensity. Moreover, microsatellites detected low but significant differentiation between Iberian localities (Celeiro and Olhão), and Aveiro differed significantly from Strangford Lough. Overall, results suggest that two independently evolving regions exist and that management strategies should be designed for each region.