Analysis of three multigene families as useful tools in species characterization of two closely-related species, Dicentrarchus labrax, Dicentrarchus punctatus and their hybrids

U2 and U4 snDNA Comparative Chromosomal Mapping in the Neotropical Fish Genera Apareiodon and Parodon (Characiformes: Parodontidae)

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis, and Parodon hilarii. However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon, aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon. Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW1W2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Comparative Cytogenetic and Sequence Analysis of U Small Nuclear RNA Genes in Three Ancistrus Species (Siluriformes: Loricariidae)

Ancistrus presents a wide karyotypic diversity, resulting from numeric and structural chromosomal rearrangements. It has been proposed that some genome-specific regions containing repetitive units could organize prone-to-break DNA sites in Loricariidae, triggering chromosomal rearrangements such as Robertsonian fusions (Rb fusions), centric fissions, translocations, and inversions. The tandemly repeats of the small nuclear RNAs (snRNAs) gene families are considered good cytogenetic markers for understanding chromosomal remodeling events among closely related species, but these snRNAs have been scarcely analyzed in Ancistrus. This study presented the nucleotide sequencing and comparative in situ location of U snRNA sequences from Ancistrus aguaboensis, Ancistrus cf. multispinis, and Ancistrus sp. (2n = 50, 52, and 50, respectively), aiming to provide information about snRNA clusters in the genome and chromosome evolution in Ancistrus. U snRNA nucleotide sequences of Ancistrus presented identity to orthologous copies and folded their secondary structures correctly. In situ localization and karyotyping of the three Ancistrus species revealed clustered copies of U2 and U5 snRNA gene families to a single chromosome site, one chromosome pair bearing U1 snRNA sequence, and one main locus of U4 snRNA sequence, besides scattered signals along the chromosomes. Previous studies related the participation of the rRNA gene families in centric fusion events, contributing to chromosome rearrangements and karyotype plasticity present in Loricariidae. In this study, homeologies in U snRNA loci chromosomal locations were detected, indicating the occurrence of conserved sites of these gene families in these three Ancistrus species with 2n = 50 or 52 chromosomes.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

Karyotypic Diversification in Two Megaleporinus Species (Characiformes, Anostomidae) Inferred from In Situ Localization of Repetitive DNA Sequences

Anostomidae species have conserved diploid numbers (2n = 54), although comparative cytogenetic studies have demonstrated chromosomal rearrangements occurrence among them, especially in repetitive DNA rich regions. The location and distribution of ribosomal DNA (rDNA) and small nuclear RNAs (snRNAs) multigene families are highly dynamic in the genomes of several organisms. In this study, we in situ located the rDNA and snRNA sites in two populations of Megaleporinus obtusidens and a sample of Megaleporinus reinhardti to infer their chromosomal changes in the evolutionary lineages. Both species of Megaleporinus shared 2n = 54 chromosomes with the presence of ZZ/ZW sex chromosome system, but they diverged in relationship to the location of 5S and 45S rDNAs as well as the distribution of snRNAs sites. The characterization of the analyzed sequences revealed the presence of complete rDNA and snRNAs sequences as well as snRNAs containing transposable elements (TEs) and microsatellite repeats. After chromosomal mapping, the sequences encompassing TEs proved to be dispersed through autosomes and accumulated on sex chromosomes. The data demonstrate that intra- and interspecific chromosomal changes occurred involving the multigene family's sites in Megaleporinus karyotypes. Furthermore, we detected TE-like sequences in the differentiation of sex chromosome systems in M. obtusidens and M. reinhardti.

Contrasting Evolutionary Paths Among Indo-Pacific Pomacentrus Species Promoted by Extensive Pericentric Inversions and Genome Organization of Repetitive Sequences

Pomacentrus (damselfishes) is one of the most characteristic groups of fishes in the Indo-Pacific coral reef. Its 77 described species exhibit a complex taxonomy with cryptic lineages across their extensive distribution. Periods of evolutionary divergences between them are very variable, and the cytogenetic events that followed their evolutionary diversification are largely unknown. In this respect, analyses of chromosomal divergence, within a phylogenetic perspective, are particularly informative regarding karyoevolutionary trends. As such, we conducted conventional cytogenetic and cytogenomic analyses in four Pomacentrus species (Pomacentrus similis, Pomacentrus auriventris, Pomacentrus moluccensis, and Pomacentrus cuneatus), through the mapping of repetitive DNA classes and transposable elements, including 18S rDNA, 5S rDNA, (CA)₁₅, (GA)₁₅, (CAA)₁₀, Rex6, and U2 snDNA as markers. P. auriventris and P. similis, belonging to the Pomacentrus coelestis complex, have indistinguishable karyotypes (2n = 48; NF = 48), with a peculiar syntenic organization of ribosomal genes. On the other hand, P. moluccensis and P. cuneatus, belonging to another clade, exhibit very different karyotypes (2n = 48, NF = 86 and 92, respectively), with a large number of bi-armed chromosomes, where multiple pericentric inversions played a significant role in their karyotype organization. In this sense, different chromosomal pathways followed the phyletic diversification in the Pomacentrus genus, making possible the characterization of two well-contrasting species groups regarding their karyotype features. Despite this, pericentric inversions act as an effective postzygotic barrier in many organisms, which appear to be also the case for P. moluccensis and P. cuneatus; the extensive chromosomal similarities in the two species of P. coelestis complex suggest minor participation of chromosomal postzygotic barriers in the phyletic diversification of these species.

Evolutionary Dynamics of rDNAs and U2 Small Nuclear DNAs in Triportheus (Characiformes, Triportheidae): High Variability and Particular Syntenic Organization

Multigene families correspond to a group of genes tandemly repeated, showing enormous diversity in both number of units and genomic organization. In fishes, unlike rDNAs that have been well explored in cytogenetic studies, U2 small nuclear RNA (snRNA) genes are poorly investigated concerning their chromosomal localization. All Triportheus species (Characiformes, Triportheidae) studied so far carry a ZZ/ZW sex chromosomes system, where the W chromosome contains a huge 18S rDNA cistron. In some species the syntenic organization of rDNAs on autosomes was also verified. To explore this particular organization, we performed three-color-fluorescence in situ hybridization using 5S, 18S rDNA, and U2 snRNA genes as probes in eight Triportheus species. This work represents the first one analyzing the chromosomal distribution of U2 snRNA genes in genomes of Triportheidae. The variability in number of rDNA clusters, and the divergent syntenies for these three multigene families, put in evidence their evolutionary dynamism, revealing a much more complex organization of these genes than previously supposed for closely related species. Our study also provides additional data on the accumulation of repetitive sequences in the sex-specific chromosome. Besides, the chromosomal organization of U2 snDNAs among fish species is also reviewed.

Chromosome mapping of repetitive DNAs in sergeant major fishes (Abudefdufinae, Pomacentridae): a general view on the chromosomal conservatism of the genus

Species of the Abudefduf genus (sergeant-majors) are widely distributed in the Indian, Pacific and Atlantic oceans, with large schools inhabiting rocky coastal regions and coral reefs. This genus consists of twenty recognized species are of generalist habit, showing typical characteristics of colonizers. Some populations maintain gene flow between large oceanic areas, a condition that may influence their cytogenetic features. A number of species have been shown to be invaders and able to hybridize with local species. However, cytogenetic data in this genus are restricted to few species. In this way, the present study includes the chromosomal investigation, using conventional (Giemsa staining, Ag-NOR and C-banding) and molecular (in situ mapping of six different repetitive DNA classes) approaches in four Abudefduf species from different oceanic regions (A. bengalensis and A. sexfasciatus from the Indo-Pacific, A. vaigiensis from the Indian and A. saxatilis from the Atlantic oceans, respectively), to investigate the evolutionary events associated with the chromosomal diversification in this group. All species share a similar karyotype (2n = 48; NF = 52), except A. sexfasciatus (2n = 48; NF = 50), which possesses a characteristic pericentric inversion in the NOR-bearing chromosomal pair. Mapping of repetitive sequences suggests a chromosomal conservatism in this genus. The high karyotypic similarity between allopatric species of Abudefduf may be related to the success of natural viable hybrids among species with recent secondary contact.

Microsatellite Organization in the B Chromosome and A Chromosome Complement in Astyanax (Characiformes, Characidae) Species

The organization of microsatellites in B and sex chromosomes has been linked to chromosomal evolution in a number of animal groups. Here, the chromosomal organizations of (CA)15, (GA)15, (CG)15, (GACA)4, and (GATA)8 microsatellites were examined in several Astyanax species with different diploid numbers: Astyanax mexicanus (2n = 50 + 1 B chromosome), A. altiparanae (2n = 50), A. marionae (2n = 48), A. fasciatus (2n = 46), and A. schubarti (2n = 36). The (CA)15 and (GA)15 microsatellites were dispersed across the chromosomes of A. altiparanae and A. fasciatus but were also observed as clusters (CA and GA for A. altiparanae, and CA for A. fasciatus). In A. marionae and A. schubarti, the (CA)15 and (GA)15 microsatellites were dispersed but were also observed as clustered signals and coincident with heterochromatic regions. In all 4 of these species, the (CG)15 and (GACA)4 microsatellites were dispersed across chromosomes, and the (GATA)8 microsatellite was co-localized with 5S rDNA. In A. mexicanus, the (CA)15, (GA)15, (CG)15, (GATA)8, and (GACA)4 microsatellites were weakly detected and dispersed across the chromosomes of the A complement. On the B chromosome, signals for the different microsatellites were weak, strong, absent, weak, and absent, respectively. The distribution of microsatellites and the locational relationship between microsatellites and 5S rDNA are discussed, and a possible evolutionary pathway is proposed for microsatellites in Astyanax.

Characterization of three different clusters of 18S-26S ribosomal DNA genes in the sea urchin P. lividus: Genetic and epigenetic regulation synchronous to 5S rDNA

We previously reported the characterization 5S ribosomal DNA (rDNA) clusters in the common sea urchin Paracentrotus lividus and demonstrated the presence of DNA methylation-dependent silencing of embryo specific 5S rDNA cluster in adult tissue. In this work, we show genetic and epigenetic characterization of 18S-26S rDNA clusters in this specie. The results indicate the presence of three different 18S-26S rDNA clusters with different Non-Transcribed Spacer (NTS) regions that have different chromosomal localizations. Moreover, we show that the two largest clusters are hyper-methylated in the promoter-containing NTS regions in adult tissues, as in the 5S rDNA. These findings demonstrate an analogous epigenetic regulation in small and large rDNA clusters and support the logical synchronism in building ribosomes. In fact, all the ribosomal RNA genes must be synchronously and equally transcribed to perform their unique final product.

Molecular Cytogenetic Analysis of the European Hake Merluccius merluccius (Merlucciidae, Gadiformes): U1 and U2 snRNA Gene Clusters Map to the Same Location

The European hake (Merluccius merluccius) is a highly valuable and intensely fished species in which a long-term alive stock has been established in captivity for aquaculture purposes. Due to their huge economic importance, genetic studies on hakes were mostly focused on phylogenetic and phylogeographic aspects; however chromosome numbers are still not described for any of the fifteen species in the genus Merluccius. In this work we report a chromosome number of 2n = 42 and a karyotype composed of three meta/submetacentric and 18 subtelo/telocentric chromosome pairs. Telomeric sequences appear exclusively at both ends of every single chromosome. Concerning rRNA genes, this species show a single 45S rDNA cluster at an intercalary location on the long arm of subtelocentric chromosome pair 12; the single 5S rDNA cluster is also intercalary to the long arm of chromosome pair 4. While U2 snRNA gene clusters map to a single subcentromeric position on chromosome pair 13, U1 snRNA gene clusters seem to appear on almost all chromosome pairs, but showing bigger clusters on pairs 5, 13, 16, 17 and 19. The brightest signals on pair 13 are coincident with the single U2 snRNA gene cluster signals. Therefore, the use of these probes allows the unequivocal identification of at least 7 of the chromosome pairs that compose the karyotype of Merluccius merluccius thus opening the way to integrate molecular genetics and cytological data on the study of the genome of this important species.

Repetitive DNA Sequences and Evolution of ZZ/ZW Sex Chromosomes in Characidium (Teleostei: Characiformes)

Characidium constitutes an interesting model for cytogenetic studies, since a large degree of karyotype variation has been detected in this group, like the presence/absence of sex and supernumerary chromosomes and variable distribution of repetitive sequences in different species/populations. In this study, we performed a comparative cytogenetic analysis in 13 Characidium species collected at different South American river basins in order to investigate the karyotype diversification in this group. Chromosome analyses involved the karyotype characterization, cytogenetic mapping of repetitive DNA sequences and cross-species chromosome painting using a W-specific probe obtained in a previous study from Characidium gomesi. Our results evidenced a conserved diploid chromosome number of 2n = 50, and almost all the species exhibited homeologous ZZ/ZW sex chromosomes in different stages of differentiation, except C. cf. zebra, C. tenue, C. xavante and C. stigmosum. Notably, some ZZ/ZW sex chromosomes showed 5S and/or 18S rDNA clusters, while no U2 snDNA sites could be detected in the sex chromosomes, being restricted to a single chromosome pair in almost all the analyzed species. In addition, the species Characidium sp. aff. C. vidali showed B chromosomes with an inter-individual variation of 1 to 4 supernumerary chromosomes per cell. Notably, these B chromosomes share sequences with the W-specific probe, providing insights about their origin. Results presented here further confirm the extensive karyotype diversity within Characidium in contrast with a conserved diploid chromosome number. Such chromosome differences seem to constitute a significant reproductive barrier, since several sympatric Characidium species had been described during the last few years and no interespecific hybrids were found.

Chromosomal Mapping of Repetitive DNA Sequences in Five Species of Astyanax (Characiformes, Characidae) Reveals Independent Location of U1 and U2 snRNA Sites and Association of U1 snRNA and 5S rDNA

Astyanax is a genus of Characidae fishes currently composed of 155 valid species. Previous cytogenetic studies revealed high chromosomal diversification among them, and several studies have been performed using traditional cytogenetic techniques to investigate karyotypes and chromosomal locations of 18S and 5S rDNA genes. However, only a few studies are currently available about other repetitive sequences. Here, the chromosomal location of small nuclear RNA genes, identified as U1 and U2 snRNA clusters, was established and compared to the distribution of 5S rDNA and histone clusters in 5 Astyanax species (A. paranae, A. fasciatus, A. bockmanni, A. altiparanae, and A. jordani) using FISH. The cytogenetic mapping of U1 and U2 snRNA demonstrated a conserved pattern in the number of sites per genome independent of the location in Astyanax species. The location of the U1 snRNA gene was frequently associated with 5S rDNA sequences, indicating a possible interaction between the distinct repetitive DNA families. Finally, comparisons involving the location of U1 and U2 snRNA clusters in the chromosomes of Astyanax species revealed a very diverse pattern, suggesting that many rearrangements have occurred during the diversification process of this group.

De novo assembly, characterization and functional annotation of Senegalese sole (Solea senegalensis) and common sole (Solea solea) transcriptomes: integration in a database and design of a microarray

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-952) contains supplementary material, which is available to authorized users.

Primary analysis of repeat elements of the Asian seabass (Lates calcarifer) transcriptome and genome

As part of our Asian seabass genome project, we are generating an inventory of repeat elements in the genome and transcriptome. The karyotype showed a diploid number of 2n = 24 chromosomes with a variable number of B-chromosomes. The transcriptome and genome of Asian seabass were searched for repetitive elements with experimental and bioinformatics tools. Six different types of repeats constituting 8–14% of the genome were characterized. Repetitive elements were clustered in the pericentromeric heterochromatin of all chromosomes, but some of them were preferentially accumulated in pretelomeric and pericentromeric regions of several chromosomes pairs and have chromosomes specific arrangement. From the dispersed class of fish-specific non-LTR retrotransposon elements Rex1 and MAUI-like repeats were analyzed. They were wide-spread both in the genome and transcriptome, accumulated on the pericentromeric and peritelomeric areas of all chromosomes. Every analyzed repeat was represented in the Asian seabass transcriptome, some showed differential expression between the gonads. The other group of repeats analyzed belongs to the rRNA multigene family. FISH signal for 5S rDNA was located on a single pair of chromosomes, whereas that for 18S rDNA was found on two pairs. A BAC-derived contig containing rDNA was sequenced and assembled into a scaffold containing incomplete fragments of 18S rDNA. Their assembly and chromosomal position revealed that this part of Asian seabass genome is extremely rich in repeats containing evolutionarily conserved and novel sequences. In summary, transcriptome assemblies and cDNA data are suitable for the identification of repetitive DNA from unknown genomes and for comparative investigation of conserved elements between teleosts and other vertebrates.

A rapid method for differentiating four species of the Engraulidae (anchovy) family

The four species of the Engraulidae family: European anchovy (Engraulis encrasicolus), Californian anchovy (Engraulis mordax), Peruvian anchoveta (Engraulis ringens), and Japanese anchovy (Engraulis japonicus) studied in this work are very similar morphologically, and it is very difficult to distinguish between them, especially when frozen or processed. We have used the 5S rDNA as a molecular marker to discriminate these four species and used specific primers designed for each species in the nontranscribed spacers (NTS) of these genes. Multiplex PCR was performed with three pairs of primers, and three different sizes were obtained: 597 bp E. encrasicolus, 598 bp E. japonicus, 380 bp E. mordax, and 250 bp E. ringens. For the species E. encrasicolus and E. japonicus, PCR-RFLP was used as an additional technique to distinguish between them because their NTS sequences showed considerable similarity.

High evolutionary dynamism in 5S rDNA of fish: state of the art

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.

Evidence for 5S rDNA horizontal transfer in the toadfish Halobatrachus didactylus (Schneider, 1801) based on the analysis of three multigene families

BACKGROUND

The Batrachoididae family is a group of marine teleosts that includes several species with more complicated physiological characteristics, such as their excretory, reproductive, cardiovascular and respiratory systems. Previous studies of the 5S rDNA gene family carried out in four species from the Western Atlantic showed two types of this gene in two species but only one in the other two, under processes of concerted evolution and birth-and-death evolution with purifying selection. Here we present results of the 5S rDNA and another two gene families in Halobatrachus didactylus, an Eastern Atlantic species, and draw evolutionary inferences regarding the gene families. In addition we have also mapped the genes on the chromosomes by two-colour fluorescence in situ hybridization (FISH).

RESULTS

Two types of 5S rDNA were observed, named type α and type β. Molecular analysis of the 5S rDNA indicates that H. didactylus does not share the non-transcribed spacer (NTS) sequences with four other species of the family; therefore, it must have evolved in isolation. Amplification with the type β specific primers amplified a specific band in 9 specimens of H. didactylus and two of Sparus aurata. Both types showed regulatory regions and a secondary structure which mark them as functional genes. However, the U2 snRNA gene and the ITS-1 sequence showed one electrophoretic band and with one type of sequence. The U2 snRNA sequence was the most variable of the three multigene families studied. Results from two-colour FISH showed no co-localization of the gene coding from three multigene families and provided the first map of the chromosomes of the species.

CONCLUSIONS

A highly significant finding was observed in the analysis of the 5S rDNA, since two such distant species as H. didactylus and Sparus aurata share a 5S rDNA type. This 5S rDNA type has been detected in other species belonging to the Batrachoidiformes and Perciformes orders, but not in the Pleuronectiformes and Clupeiformes orders. Two hypotheses have been outlined: one is the possible vertical permanence of the shared type in some fish lineages, and the other is the possibility of a horizontal transference event between ancient species of the Perciformes and Batrachoidiformes orders. This finding opens a new perspective in fish evolution and in the knowledge of the dynamism of the 5S rDNA. Cytogenetic analysis allowed some evolutionary trends to be roughed out, such as the progressive change in the U2 snDNA and the organization of (GATA)n repeats, from dispersed to localized in one locus. The accumulation of (GATA)n repeats in one chromosome pair could be implicated in the evolution of a pair of proto-sex chromosomes. This possibility could situate H. didactylus as the most highly evolved of the Batrachoididae family in terms of sex chromosome biology.

Genetic characterization of Plectorhinchus mediterraneus yields important clues about genome organization and evolution of multigene families

Background

Molecular and cytogenetic markers are of great use for to fish characterization, identification, phylogenetics and evolution. Multigene families have proven to be good markers for a better understanding of the variability, organization and evolution of fish species. Three different tandemly-repeated gene families (45S rDNA, 5S rDNA and U2 snDNA) have been studied in Plectorhinchus mediterraneus (Teleostei: Haemulidae), at both molecular and cytogenetic level, to elucidate the taxonomy and evolution of these multigene families, as well as for comparative purposes with other species of the family.

Results

Four different types of 5S rDNA were obtained; two of them showed a high homology with that of Raja asterias, and the putative implication of a horizontal transfer event and its consequences for the organization and evolution of the 5S rDNA have been discussed. The other two types do not resemble any other species, but in one of them a putative tRNA-derived SINE was observed for the first time, which could have implications in the evolution of the 5S rDNA. The ITS-1 sequence was more related to a species of another different genus than to that of the same genus, therefore a revision of the Hamulidae family systematic has been proposed. In the analysis of the U2 snDNA, we were able to corroborate that U2 snDNA and U5 snDNA were linked in the same tandem array, and this has interest for tracing evolutionary lines. The karyotype of the species was composed of 2n = 48 acrocentric chromosomes, and each of the three multigene families were located in different chromosome pairs, thus providing three different chromosomal markers.

Conclusions

Novel data can be extracted from the results: a putative event of horizontal transfer, a possible tRNA-derived SINE linked to one of the four 5S rDNA types characterized, and a linkage between U2 and U5 snDNA. In addition, a revision of the taxonomy of the Haemulidae family has been suggested, and three cytogenetic markers have been obtained. Some of these results have not been described before in any other fish species. New clues about the genome organization and evolution of the multigene families are offered in this study.