Senegalese sole (Solea senegalensis) and common sole (S. solea) are two economically and evolutionary important flatfish species both in fisheries and aquaculture. Although some genomic resources and tools were recently described in these species, further sequencing efforts are required to establish a complete transcriptome, and to identify new molecular markers. Moreover, the comparative analysis of transcriptomes will be useful to understand flatfish evolution. A comprehensive characterization of the transcriptome for each species was carried out using a large set of Illumina data (more than 1,800 millions reads for each sole species) and 454 reads (more than 5 millions reads only in S. senegalensis), providing coverages ranging from 1,384x to 2,543x. After a de novo assembly, 45,063 and 38,402 different transcripts were obtained, comprising 18,738 and 22,683 full-length cDNAs in S. senegalensis and S. solea, respectively. A reference transcriptome with the longest unique transcripts and putative non-redundant new transcripts was established for each species. A subset of 11,953 reference transcripts was qualified as highly reliable orthologs (>97% identity) between both species. A small subset of putative species-specific, lineage-specific and flatfish-specific transcripts were also identified. Furthermore, transcriptome data permitted the identification of single nucleotide polymorphisms and simple-sequence repeats confirmed by FISH to be used in further genetic and expression studies. Moreover, evidences on the retention of crystallins crybb1, crybb1-like and crybb3 in the two species of soles are also presented. Transcriptome information was applied to the design of a microarray tool in S. senegalensis that was successfully tested and validated by qPCR. Finally, transcriptomic data were hosted and structured at SoleaDB. Transcriptomes and molecular markers identified in this study represent a valuable source for future genomic studies in these economically important species. Orthology analysis provided new clues regarding sole genome evolution indicating a divergent evolution of crystallins in flatfish. The design of a microarray and establishment of a reference transcriptome will be useful for large-scale gene expression studies. Moreover, the integration of transcriptomic data in the SoleaDB will facilitate the management of genomic information in these important species.

The 5S ribosomal DNA (rDNA) consists of one transcriptional unit of about 120 base pairs, which is separated from the next unit by a non-transcribed spacer (NTS). The coding sequence and the NTS together form a repeat unit which can be found in hundreds to thousands of copies tandemly repeated in the genomes. The NTS regions seem to be subject to rapid evolution. The first general model of evolution of these multigene families was referred to as divergent evolution, based on studies using hemoglobin and myoglobin as model systems. Later studies showed that nucleotide sequences of different multigene family members are more closely related within species than between species. This observation led to a new model of multigene family evolution, termed concerted evolution. Another model of evolution, named the birth-and-death model, has been found to be more suitable to explain the long-term evolution of these multigene families. According to this model, new genes originate by successive duplications, and these new genes are either maintained for a long time or are lost, or else degenerate into pseudogenes. In this review we describe different sources of variability in the 5S rDNA genes observed in several distinct fish species. This variability is mainly referred to NTSs and includes the presence of other multigene families (mainly LINEs, SINEs, non-LTR retrotransposons, and U snRNA families). Different types of microsatellites have also been found to contribute to the increase of variability in this region. Our recent results suggest that horizontal transfer contributes to the increase of diversity in the NTSs of some species. Variability in the 5S rDNA coding region affecting the stability of the structure, but without effects on the function of the 5S rRNA, is also described. Retrotransposons seem to be responsible for the high dynamism of 5S rDNA, while microsatellites acting as recombination hot spots could stabilize a wide variety of unusual DNA structures, affecting DNA replication and enhancing or decreasing promoter activity in gene expression. The relationship between the high variability found at molecular level and the low variability found at chromosomal level is also discussed.

Some units of the 5S rDNA of Solea senegalensis were amplified by PCR and sequenced. Three main PCR products (227, 441, and 2166 bp) were identified. The 227- and 441-bp fragments were characterized by highly divergent nontranscribed spacer sequences (referred to as NTS-I and NTS-II) that were 109 and 324 bp long, respectively, yet their coding sequences were nearly identical. The 2166-bp 5S rDNA unit was composed of two 5S rRNA genes separated by NTS-I and followed by a 1721-bp spacer containing the U2, U5, and U1 small nuclear RNA genes (snRNAs). They were inverted and arranged in the transcriptional direction opposite that of the 5S rRNA gene. This simultaneous linkage of 3 different snRNAs had never been observed before. The PCR products were used as probes in fluorescence in situ hybridization experiments to locate the corresponding loci on the chromosomes of S. senegalensis. A major 5S rDNA chromosomal site was located along most of the short arm of a submetacentric pair, while a minor site was detected near the centromeric region of an acrocentric pair.Key words: soleidae, pleuronectiformes, 5S rDNA, Solea, snRNAs linkage.