Evolution of the Proto Sex-Chromosome in Solea senegalensis

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

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.

Chromosome anchoring in Senegalese sole (Solea senegalensis) reveals sex-associated markers and genome rearrangements in flatfish

The integration of physical and high-density genetic maps is a very useful approach to achieve chromosome-level genome assemblies. Here, the genome of a male Senegalese sole (Solea senegalensis) was de novo assembled and the contigs were anchored to a high-quality genetic map for chromosome-level scaffolding. Hybrid assembled genome was 609.3 Mb long and contained 3403 contigs with a N50 of 513 kb. The linkage map was constructed using 16,287 informative SNPs derived from ddRAD sequencing in 327 sole individuals from five families. Markers were assigned to 21 linkage groups with an average number of 21.9 markers per megabase. The anchoring of the physical to the genetic map positioned 1563 contigs into 21 pseudo-chromosomes covering 548.6 Mb. Comparison of genetic and physical distances indicated that the average genome-wide recombination rate was 0.23 cM/Mb and the female-to-male ratio 1.49 (female map length: 2,698.4 cM, male: 2,036.6 cM). Genomic recombination landscapes were different between sexes with crossovers mainly concentrated toward the telomeres in males while they were more uniformly distributed in females. A GWAS analysis using seven families identified 30 significant sex-associated SNP markers located in linkage group 18. The follicle-stimulating hormone receptor appeared as the most promising locus associated with sex within a region with very low recombination rates. An incomplete penetrance of sex markers with males as the heterogametic sex was determined. An interspecific comparison with other Pleuronectiformes genomes identified a high sequence similarity between homologous chromosomes, and several chromosomal rearrangements including a lineage-specific Robertsonian fusion in S. senegalensis.

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.

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.

Development of whole-genome multiplex assays and construction of an integrated genetic map using SSR markers in Senegalese sole

The Senegalese sole (Solea senegalensis) is an economically important flatfish species. In this study, a genome draft was analyzed to identify microsatellite (SSR) markers for whole-genome genotyping. A subset of 224 contigs containing SSRs were preselected and validated by using a de novo female hybrid assembly. Overall, the SSR density in the genome was 886.7 markers per megabase of genomic sequences and the dinucleotide motif was the most abundant (52.4%). In silico comparison identified a set of 108 SSRs (with di-, tetra- or pentanucleotide motifs) widely distributed in the genome and suitable for primer design. A total of 106 markers were structured in thirteen multiplex PCR assays (with up to 10-plex) and the amplification conditions were optimized with a high-quality score. Main genetic diversity statistics and genotyping reliability were assessed. A subset of 40 high polymorphic markers were selected to optimize four supermultiplex PCRs (with up to 11-plex) for pedigree analysis. Theoretical exclusion probabilities and real parentage allocation tests using parent–offspring information confirmed their robustness and effectiveness for parental assignment. These new SSR markers were combined with previously published SSRs (in total 229 makers) to construct a new and improved integrated genetic map containing 21 linkage groups that matched with the expected number of chromosomes. Synteny analysis with respect to C. semilaevis provided new clues on chromosome evolution in flatfish and the formation of metacentric and submetacentric chromosomes in Senegalese sole.

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.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

High dynamism for neo-sex chromosomes: satellite DNAs reveal complex evolution in a grasshopper

A common characteristic of sex chromosomes is the accumulation of repetitive DNA, which accounts for their diversification and degeneration. In grasshoppers, the X0 sex-determining system in males is considered ancestral. However, in some species, derived variants like neo-XY in males evolved several times independently by Robertsonian translocation. This is the case of Ronderosia bergii, in which further large pericentromeric inversion in the neo-Y also took place, making this species particularly interesting for investigating sex chromosome evolution. Here, we characterized the satellite DNAs (satDNAs) and transposable elements (TEs) of the species to investigate the quantitative differences in repeat composition between male and female genomes putatively associated with sex chromosomes. We found a total of 53 satDNA families and 56 families of TEs. The satDNAs were 13.5% more abundant in males than in females, while TEs were just 1.02% more abundant in females. These results imply differential amplification of satDNAs on neo-Y chromosome and a minor role of TEs in sex chromosome differentiation. We showed highly differentiated neo-XY sex chromosomes owing to major amplification of satDNAs in neo-Y. Furthermore, chromosomal mapping of satDNAs suggests high turnover of neo-sex chromosomes in R. bergii at the intrapopulation level, caused by multiple paracentric inversions, amplifications, and transpositions. Finally, the species is an example of the action of repetitive DNAs in the generation of variability for sex chromosomes after the suppression of recombination, and helps understand sex chromosome evolution at the intrapopulation level.

Gene clusters related to metamorphosis in Solea senegalensis are highly conserved

The flatfish, Solea senegalensis has considerable scientific interest and commercial value. The metamorphosis in this species occurs between 12 and 19 days after hatching and it takes about 1 week to complete. Eleven Bacterial Artificial Chromosomes (BAC) clones containing the various candidate genes involved in the process of metamorphosis: thyroxine 5 deiodinase 3 (dio3); forkhead box protein E4 (foxe4); melatonin receptor type 1C (mel1c); calsequestrin 1b (casq1b); thyrotropin subunit beta (tshβ); thyrotropin-releasing hormone receptor 1, 2, and 3 (trhr1, trhr2, trhr3); thyroid hormone receptor α a and b (thrαa, thrαb); and thyroid hormone receptor beta (thrβ) were analyzed by multiple Fluorescence in situ Hybridization (mFISH) and Next Generation Sequencing (NGS) techniques. The mFISH technique localized the 11 BAC clones on 12 different chromosome pairs because three of them, specifically the trhr1a, trhr2 and thrβ BAC clones, showed double signals. This signal duplication indicates a duplication of the genomic region inserted within the BAC clone, which provides evidence for the Teleost-Specific Whole Genome Duplication (TS-WGD). Micro-synteny and phylogenetic analysis showed that Cynoglossus semilaevis is the nearest species to S. senegalensis and that Danio rerio is the most distant one. The tshβ BAC clone was highly conserved as the genes belonging to this BAC were located on a single chromosome in all the species studied. These genes participate in proliferation, migration and cell-death, which are key processes during metamorphosis. Overall, micro-synteny analysis showed that most candidate genes are found in conserved genomic surroundings.