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

  • Manuel A Merlo1,

    Affiliated with

    • Ismael Cross1,

      Affiliated with

      • José L Palazón2, 3,

        Affiliated with

        • María Úbeda-Manzanaro1, 2,

          Affiliated with

          • Carmen Sarasquete2 and

            Affiliated with

            • Laureana Rebordinos1Email author

              Affiliated with

              BMC Evolutionary Biology201212:201

              DOI: 10.1186/1471-2148-12-201

              Received: 29 March 2012

              Accepted: 3 October 2012

              Published: 7 October 2012



              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).


              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.


              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.


              The Lusitanian toadfish (Halobatrachus didactylus), a member of the Batrachoididae family, is found in natural habitats of soft sand and rocky sea beds, and reefs. The Batrachoididae family comprises 71 species, belonging to three subfamilies [1] that are widely distributed across the Atlantic, Pacific, and Indian Oceans. H. didactylus is the only species of this family found in the Iberian Peninsula [2] and populations are distributed along the coasts from the Bay of Biscay (Spain) to Ghana as well as the western Mediterranean [3].

              On the coasts of Portugal and southwestern Spain, this species is of moderate commercial importance. However, the real importance of this species is for research purposes rather than its commercial value. It has traditionally been used as a model animal in toxicology experiments [4], and in hematology, reproduction and histophysiology studies [57]. Genetic studies in the Batrachoididae family have been focused mainly on cytogenetics (reviewed by [8]); there is only one published paper in which this fish family has been studied from a molecular perspective [9]. For H. didactylus in particular, the karyotype and the localization of various repetitive sequences have been described [10].

              Multigene families are comprised of multiple genes that all descend from a common ancestral gene, have similar sequences and are functionally related [11]. Traditionally, it has been presumed that the evolution of multigene families is in accordance with a model known as the “Concerted Evolution Model”, in which the members of the family are homogenized by several different mechanisms such as unequal crossing-over and gene conversion [12]. This feature makes the units show more similarity within than between related species [13]. In the last two decades, a new model has been proposed to explain controversial situations found in some multigene families: the “Birth-and-Death Evolution Model”. Under this model new genes arise by duplications during the evolution of a group of organisms, and these new genes are either fixed as functional genes or become pseudogenes [14]. The histone and the immunoglobulin gene families are typical examples of birth-and-death evolution [11].

              The ribosomal genes are divided in two clusters: the major (45S rDNA) and the minor (5S rDNA) clusters. The 45S rDNA has a transcribed unit which comprises two external transcribed spacers (5’ ETS and 3’ ETS), the coding for 18S, 5.8S and 28S rRNAs, and two internal transcribed spacers, which separate the 18S rRNA from the 5.8S rRNA (ITS-1) and the 5.8S rRNA from the 28S rRNA (ITS-2); the transcribed units are separated by intergenic spacers (IGS). Meanwhile, the 5S rDNA comprises a conserved coding region of 120 bp and a non-transcribed spacer (NTS) which is variable among species in length and sequence. The two ribosomal clusters represent well-known examples of concerted evolution.

              The small nuclear RNAs (snRNA) are components of the small nuclear ribonucleoprotein particles (snRNP), a complex which participate in the splicing process of the mRNA precursors [12]. Most of the snRNA are transcribed by RNA polymerase II, except for U6 snRNA which is transcribed by RNA polymerase III [15]. Within the group of fishes, there have been very few studies made of the snRNA gene family at molecular and cytogenetic level. It has been presumed that the snRNA gene families undergo a concerted evolution [11].

              In this work a molecular characterization has been made of three multigene families: 5S rDNA, ITS-1 (from 45S rDNA) and the U2 snRNA gene of H. didactylus. In addition, double-FISH technique was used to ascertain the U2 snRNA gene localization and the possible co-localization of this probe with the 5S rDNA, 18S rRNA gene and (GATA)n probes. The co-localization of (GATA)n and the 18S rRNA gene was also tested. This technique not only provides various different cytogenetic markers, but also represents a starting point for the production of a genetic map in H. didactylus. The nucleotide variation has been estimated to assess the homogenizing forces which support the concerted evolution in the three multigene families. Finally, the results have been compared with those obtained from other members of the Batrachoididae family, and the evolutionary difference between those and H. didactylus has been established.


              Sequence analysis

              After PCR amplification, two electrophoretic bands of 5S rDNA sequence could be observed, named type α and type β, and whose sizes were 420 and 200 bp respectively. The amplification with the type β specific primer gave a 190 bp band in 9 specimens of H. didactylus (data not shown). The sequencing results demonstrate that the 5S rDNA sizes varied from 413 to 424 bp in type α, and from 197 to 207 bp in type β 5S rDNA sequences. In the former the variability was due a variable poly-A region at position −71, and in the latter by a 10 bp insertion-deletion (indel) localized immediately after the poly-T termination signal. In two clones an electrophoretic band of 840 bp was obtained; this corresponded to a dimeric form of 5S rDNA type α. From these findings, it would therefore be possible to ascertain the complete sequence for 5S rDNA of type α.

              Both types of 5S rDNA sequence conserved all promoter and regulatory regions necessary for transcription by RNA polymerase III (pol III), such as the TATA-like promoter element around the position −30 upstream from the transcription starting point, the poly-T terminator region situated at the 3’ end of the coding region and the Internal Control Regions (Box A, Intermediate Element and Box C) located inside coding regions between the positions +50 to +65, +68 to +73, and +81 to +98 respectively. The RNAstructure 5.2 program showed a stable secondary structure of 5S rRNA in both cases (type α and type β) (Figure 1), with similar free energy values, and were in accordance with secondary structures predicted in previously described models [16].
              Figure 1

              Secondary structure prediction of the two types of 5S rRNA. The free energy (ΔG) of the structure is shown, and was calculated considering the mean seawater temperature (15°C).

              Meanwhile, the ITS-1 and U2 snDNA sequences showed one electrophoretic band of 650 and 800 bp respectively. The sequencing results demonstrated that the ITS-1 sequence presented an exact size of 664 bp, and 820 bp in the U2 snDNA sequence. Furthermore, three putative promoter and regulator sequences in the U2 snRNA gene have been found (Figure 2). These include the Proximal Sequence Element (PSE), the Distal Sequence Element (DSE) and the 3’ box within the spacer region. The PSE and 3’ box were located by alignment with consensus sequences obtained from the Moronidae species, Dicentrarchus labrax and Dicentrarchus punctatus. From the alignment of 126 bp upstream from the coding region two conserved regions were detected (Figure 2a). The more proximal of these is between position −60 and −23, while the more distal region comprises a 17 bp box, the sequence of which has not previously been described as a promoter or regulator region. The 100 bp downstream from the coding region was also aligned and a more conserved box was found at 8 nucleotides downstream from the 3’ end of the coding region (Figure 2b). A putative DSE octameric sequence was localized at −311 bp from the transcription starting point (Figure 2d), and it was located in direct orientation with respect to the consensus octameric sequence [17].
              Figure 2

              U2 snRNA gene organization in Halobatrachus didactylus. (a) Different regions of the U2 snRNA represented schematically. (b) 3’ end alignment of the U2 spacer; 3’b: 3’ Box; D: Putative DSE element; RE: Regulator Element; P: PSE element. (c) 5' end alignment of the U2 spacer; Hdi: Halobatrachus didactylus; Dla: Dicentrarchus labrax; Dpu: Dicentrarchus punctatus. (d) Alignment of the DSE element with those obtained from other species.

              Variability analysis

              The three multigene families studied showed low values of nucleotide variability, in both coding regions and spacers. As expected, the spacer regions had higher values than coding regions, except in the U2 snRNA gene (Table 1). On the whole, the most variable multigene family was the U2 snDNA. The 5S rDNA type α also showed high nucleotide variability. The less variable genes were the 5S rDNA type β and the 18S rRNA-ITS-1-5.8S rRNA sequence. Furthermore, divergence values between the two types of 5S rDNA were different depending on the fragment considered in the study. Thus, the coding region showed a low value of divergence (d= 0.087 ± 0.028), whereas divergence in the NTS was higher (d= 0.751 ± 0.218).
              Table 1

              Polymorphism by gene region


              Length (bp)





              5S rRNA (α)






              NTS (α)






              5S rRNA (β)






              NTS (β)






              18S rRNA












              5.8S rRNA






              U2 snRNA






              U2 spacer






              %GC: GC content; S: number of polymorphic sites; h: number of haplotypes; π: nucleotide variability.

              The gene families under study were subjected to BLASTN search at the NCBI database; all coding regions (5S rRNA, 18S rRNA, 5.8S rRNA and U2 snRNA) showed high degrees of identity with others obtained from a wide variety of fish species. Spacer regions did not present any significant homologies, except for the 5S rRNA gene spacer (NTS) type β which, surprisingly, presented a high degree of homology with the NTS from Sparus aurata described by [18] (Figure 3). The analysis with the type β demonstrated that this type was present in all specimens of H. didactylus studied, as well as in all the species included for this experiment (see Methods section), except in Solea senegalensis and Engraulis encrasicholus (Figure 4). However, these two species showed positive bands in the parallel experiment with ITS-1 primers.
              Figure 3

              5S rDNA type β, common to Sparus aurata and Halobatrachus didactylus. Alignment of the two clones of Sparus aurata (GenBank Acc. No. AY330701.1 and AY330702.1) and the type β clones of Halobatrachus didactylus. A schematic representation of 5S rDNA is included, highlighting the A box, the Intermediate Element (IE), the C box, the poly-T terminator region, and the TATA-like box.

              Figure 4

              Agarose gel electrophoresis of 5S rDNA type β species-specific PCR. Positive products of 190 bp are marked by an arrow. Hdi: Halobatrachus didactylus; Ppl: Porichthys plectrodon; Tma: Thalassophryne maculosa; Acr: Amphichthys cryptocentrus; Bma: Batrachoides manglae; Sau: Sparus aurata (positive control); Ppa: Pagrus pagrus; Dsa: Diplodus sargus; Dla: Dicentrarchus labrax; Are: Argyrosomus regius; Pme: Plectorhinchus mediterraneus; Sse: Solea senegalensis; Een: Engraulis encrasicholus; C-: negative control.

              Phylogenetic analysis

              The Maximum Likelihood (ML) tree obtained (Figure 5) was compared with the Neighbor-Joining (NJ) tree previously described by [9] for other species of the Batrachoididae family (Western batrachoids). As in the NJ tree, the ML tree divided the Western batrachoids species in two clades, characterized by different NTS types. In the clade 2 (following the nomenclature devised by [9]) of these species, the types α and β of H. didactylus were grouped on the same branch, but within that branch, the type β clade formed a monophyletic group with the clade 2 of the Western batrachoids. The two sequences of S. aurata mentioned earlier were clustered within the H. didactylus type β clade with a high branch support. Therefore, S. aurata did not group with the two other sparid species, which formed a separate cluster from that of the Batrachoididae. Within clade 1 of the Western batrachoids, the clones of the species Amphichthys cryptocentrus were clustered in the same group, whereas in the NJ tree previously described by [9], these clones were divided into two groups.
              Figure 5

              Phylogenetic Maximum Likelihood tree of NTS sequences based on GTR substitution model. Type β clones of Halobatrachus didactylus are shown. Western batrachoid clones have been collapsed. Dicentrarchus labrax clones have been rooted as an outgroup. Sparus aurata clones are also shown, preceded by a rhombus.

              Cytogenetic analysis

              Four double-FISH experiments were carried out with the following: the U2 snRNA gene with 5S rDNA (Figure 6a); the 18S rRNA gene with U2 snRNA gene (Figure 6b); the (GATA)n sequence with U2 snRNA gene (Figure 6c); and the (GATA)n sequence with 18S rRNA gene (Figure 6d). The 5S rDNA probe only includes the type α. The double-FISH treatments with 5S-18S and 5S-GATA probes had been performed in previous work [10]. The U2 snDNA probe was localized in the middle of the long arm of a medium-size submeta/subtelocentric chromosome pair. No co-localization between the U2 snRNA gene probe and the remaining probes was detected. Double-FISH treatment with (GATA)n and 18S rRNA gene probes showed a non co-localized pattern.
              Figure 6

              Metaphase of Halobatrachus didactylus after FISH treatment. The following probes were used: (a) U2 snRNA gene (green signals) versus 5S rDNA (red signals); (b) 18S rRNA gene (green signals) versus U2 snRNA gene (red signals); (c) (GATA)n sequence (green signals) versus U2 snRNA gene (red signals); (d) (GATA)n sequence (green signals) versus 18S rRNA (red signals).


              Two types of 5S rDNA have been found in H. didactylus, named here as type α and type β. The NTS sequences between the two types do not share significant homologies. Furthermore, the degree of divergence of the coding region between the two types is more than 5 times higher than the nucleotide variability within the coding region of type α, and more than 28 times higher than such variability in type β (Table 1). Both types presented all the regulatory elements which make them functional genes. The existence of two types of 5S rDNA appears to be common in fish species [9, 19, 20], and this situation has been commonly explained as a dual expression system, where one type is expressed in both somatic and oocyte cells and the other type only in oocyte cells [21].

              In a previous study [9] with four species of the Batrachoididae family from the Western Atlantic, it has been found that two species also shared two types of 5S rDNA, and the two other species have only one type. Therefore, it cannot be concluded that a dual system of 5S rDNA is generally established in the Batrachoididae family. In the four Western Atlantic species the NTS sequences are each similar to the others (see Figure 3 of [9]), whereas in H. didactylus these sequences are distantly related to those of the four species. This feature could indicate that H. didactylus has evolved independently from the Western Atlantic species, due to the geographic isolation between them, which could have contributed to the greater divergence found between Eastern and Western Atlantic species.

              The sequences of both coding regions and the NTSs are quite conserved in H. didactylus, and the secondary structure agrees with previous described models [16] (Figure 1). Therefore, concerted evolution is the more feasible model which is acting in each locus of this multigene family. Under this model, the different units of the gene family tend to homogenize their sequences by means of mechanisms such as unequal cross-over and gene conversion [12], thus spreading mutation events along the multigene family. Concerted evolution makes the sequence more similar within the species than it is among related species [13]. Traditionally it has been established that the 5S rDNA gene family is a paradigm for the concerted evolution model; however, although this is the most common model observed for this gene, in more and more studies that are published, it is concluded that the 5S rDNA gene family can follow a different evolution model. Some authors have described 5S rDNA gene families which fit better with a birth-and-death evolution model than with the concerted model [9, 14, 22], or with a mixed model with elements of both birth-and-death and concerted evolution [23]. In the H. didactylus case, although the concerted model is acting at locus level, at genome level, the presence of two types and the between-species clustering of the multigene family indicate a birth-and-death evolution model.

              As mentioned above, surprisingly the 5S rDNA type β presented homology in the NTS sequence with those obtained in S. aurata[18] (Figure 3). This species and H. didactylus belong to different superorders (Acanthopterygii and Paracanthopterygii respectively), which diverged 55 My ago [24]. Both species could have inherited the type β from an ancestral species common to the two supeorders, and could have maintained in some lineages and lost in others along the fish evolution. Using the type β specific primer, we were able to conclude that this sequence is found in all specimens studied of H. didactylus, as well as in S. aurata and the remaining species belonging to Batrachoidiformes (Porichthys plectrodon, Thalassophryne maculosa, Amphichthys cryptocentrus and Batrachoides manglae) and Perciformes order (Pagrus pagrus, Diplodus sargus, Dicentrarchus labrax, Argyrosomus regius and Plectorhinchus mediterraneus) (Figure 4). However, the β type is not present in either the Pleuronectiformes (S. senegalensis) or Clupeiformes orders (E. encrasicholus). Nevertheless, taking into account the highly dynamism to which the NTS are subjected, and the period elapsed since the two superorders began to diverge, another theory could be applied: a Horizontal Transference (HT) between ancient species of Perciformes and Batrachoidiformes groups. Therefore, the new type found in H. didactylus was not generated by a duplication within the own genome, but rather by an “extragenomic duplication”.

              The HT phenomenon has been documented extensively in prokaryotes, and is a well-known mechanism of gene exchange [25]; it has been assumed that between 88% and 98% of the expansion of protein families in prokaryotes is due to HT [26]. However, the occurrence of HT between eukaryote organisms remains obscure [27]. It has been postulated that HT is a very important mechanism, the source for much biological innovation; indeed, some of the most significant events in cell evolution have been associated with HT, such as the origin of the primordial eukaryote cell, and HT has contributed to the characteristics of the multicellularity [28]. All documented cases of HT between eukaryotes involve transposable elements (TE), since their ability to mobilize and integrate within a genome makes them susceptible to horizontal transposable transference (HTT) [29]. Cases have been reported of HTT between a wide range of eukaryotes [3034]. To date, the transfer of host genes by TEs between different eukaryote species has not yet been observed, but it is known that TEs are able to capture and transduce sequences with considerable frequency within a species [35]. Thus, HTT could be the mechanism responsible for a lateral movement of genes between eukaryotes.

              It has been demonstrated that 5S rRNA genes and retro-transposons can interact with one another [36], and this interaction might be the cause of the pattern of evolution and the dispersed arrangement of some fungal organisms [37]. Similarly, this interaction could “open the door” for a 5S rRNA gene lateral transfer. The precise mechanism whereby a sequence is transferred horizontally from one species to another remains elusive, although various different mechanisms have been proposed in which some particles or organisms could take part as an intermediate, such as a virus [38], a bacteria [39] or a parasitic organism [40]. However, it is difficult to demonstrate the involvement of such a microorganism and virus, because of their capacity to erase any trace of non-essential DNA [28]. Few cases of HT from eukaryotes to virus [41] or bacteria [42] have been described. Typically, there are different ways to detect HT: (i) the high degree of similarity in sequences between two species phylogenetically distant; (ii) the rate of nucleotide variability of a transferred gene might be lower than that of the other orthologous genes; and (iii) in a phylogenetic analysis there should be incongruence between the phylogenies of the gene and of the organism [27]. The first way is the one used most often to ascertain the HT [28], and similarity can clearly be seen in the H. didactylus S. aurata case. Secondly, the nucleotide diversity of the putative transferred genes (type β) is lower than that of the orthologous genes (type α). Lastly, the phylogenetic analysis supports with high branch confidence an event of HT, since the S. aurata clones and type β are strongly grouped (Figure 5), thus leading to incongruence between the taxonomic situation and the NTS-based phylogeny.

              Another alternative way of horizontal transfer should not be discounted: the sperm-mediated gene transfer (SMGT). It is well-known that sperm cells are able to capture exogenous DNA and to transfer it to the oocyte at fertilization [43]. This capacity has led to sperm cells being used as vectors for transferring DNA in transgenic experiments [4446]. The frequency of sperm-exogene DNA uptake depends on several factors, such as the presence of inhibitory glycoproteins and the maintenance of optimum conditions for the DNA uptake in terms of quantity, length and primary structure of the exogene DNA [47]. It has been established that the greater the quantity of exogenous DNA, the higher the possibility of take-up by sperm cells [47, 48] reported that the DNA content in sediments is 3 to 4 orders of magnitude higher than in the water column. Taking this into account and the benthic behaviour of the majority of species of the Batrachoidiformes order, the possibility of a SMGT is plausible.

              It is commonly considered that the U2 snRNA gene also evolves in a concerted fashion [11]. However, there are few studies with this gene to conclude unequivocally that the concerted evolution is the general model followed by the U2 snRNA gene. The organization of U2 snRNA has only been investigated in humans [49], primates [50] and moronid fishes [51], and the same characteristics are shared, i.e., U2 snRNA are clustered in tandem arrays and undergo a concerted evolution model. The U2 sequence homology observed in many different specimens and units has led us to confirm a similar situation with the U2 snRNA of H. didactylus. In all the studies mentioned above (including the present work) the U2 snRNA gene is localized in one chromosome pair, although the position inside the chromosome is not conserved. A dispersed gene arrangement facilitates a birth-and-death evolution model [37]; therefore a clustered arrangement could be propitious for a concerted evolution model. Moreover, taking into account the results obtained by [8] with other species from the Batrachoididae family, it is possible that there has been a transition from a birth-and-death to a concerted evolution model within this fish family, since three distinct patterns of hybridization have been seen in only five species analyzed: dispersed, clustered, and both patterns together.

              Despite the scarcity of studies describing the molecular characterization of the U2 snRNA in fish species, the only four species characterized to date, Solea senegalensis[12], Dicentrarchus labrax and D. punctatus[51], and A. regius[52] all showed a close linkage with another member of the U snRNA family. This linkage has been applied as a useful tool for the identification of different species of sole [53]. However, in H. didactylus, the U2 was not found in linkage with any other U snRNA; this absence is, so far, the only instance found within a fish group. There are insufficient data to conclude whether the U2 linkage in fishes is a general rule or is an exceptional case of only some fish families. In either case, the analysis of this gene family could be an interesting tool for evolutionary and phylogenetic purposes.

              The search for promoter elements in the U2 spacer regions yielded four putative regions (Figure 2). The nearest to the 5’ coding region was located at −60 bp, and, as described by other authors [15, 54], the PSE element is located between positions −60 and −50; therefore, the box found at position −60 could be the PSE regulatory element. A putative 3’ box was also localized at 8 nt downstream from the coding region, which is close to the localization described by [54] for this box (9 nt downstream). The DSE enhancer element is independent of the position and orientation, so these characteristics make it difficult to find a possible DSE by aligning with other sequences. Despite this, it has been reported that the DSE element is usually found around 200–250 bp upstream from the PSE [55]. A putative DSE element has been observed in H. didactylus at −244 bp upstream from the PSE, so that could be the enhancer element. The presence of all these regulator elements favours the supposition of a transcriptionally active gene.

              The ITS-1 region without 18S and 5.8S coding regions comprises 388 bp, and the GC content is between 62.89% and 63.40%. Although this GC content is similar to the average for the Osteichthyes group (68.0% ± 4.2), the ITS-1 size is one of the smallest within this group, which varies between 318 and 1518 bp [56]. When the ITS-1 sequence was subjected to BLASTN search, an ITS-1 fragment aligned with another from many fish species, reflecting the moderate nucleotide conservation of this sequence compared with other spacer regions (IGS, NTS, etc.). This region undergoes some selective pressure, because their secondary structures are necessary for the rRNA maturation process [57]. These features make this sequence suitable for phylogenetic analysis at the species level.

              The relatively low value of the nucleotide variability in the spacer regions of the three sequences treated here (Table 1) leads us to accept that their evolution has been according to the concerted model. Unequal cross-over and gene conversion are the common mechanisms that the multigene families undergo with a concerted evolution model. However, coding regions undergo another additional mechanism, such as the purifying selection, which makes them less variable. The differences in the variability found in the coding region in each gene could be due to opposing differences between homogenizing and diversifying forces. Therefore, the dynamic of the concerted evolution process depends on two conflicting forces: the homogenizing forces and the forces which generate variability [58]. It has been postulated that the unequal cross-over and gene conversion mechanisms occur more frequently in chromosome regions adjacent to the telomeres [11], thus increasing the effectiveness of concerted evolution. On this point, cytogenetics data support the conclusions already obtained with nucleotide diversity, i.e., the concerted mechanisms are more frequent in the multigene family which is situated in the telomeric position, as is the 18S rRNA, resulting in less nucleotide variation. Conversely, the two most variable multigene families (5S rDNA and U2 snRNA gene) are localized in internal positions.

              None of the four probes used to hybridize in H. didactylus chromosomes were co-localized (Figure 6; Table 2); therefore they represent four different chromosomal markers, which is useful in a preliminary genome mapping. Comparisons with similar results obtained from other species of Batrachoididae family offer some evolutionary trends (Table 2; see also [8] and [10]). The five implicated species show a variable fundamental number (FN), which implies that large chromosomal re-arrangement events have occurred within this fish family, the most recent species being that with the highest FN [59]. Taking this into account, the most recent species is H. didactylus (FN= 72; [10]). Comparing the FISH results obtained in the present study with those from the other Batrachoididae species [8], some further relevant comments can be made. Firstly, U2 snRNA and (GATA)n sequences have been gradually changing from a dispersed organization throughout the genome, to become localized in a single locus as in H. didactylus. Second, the ancestral (plesiomorphic) localization for 5S rDNA in the Batrachoididae family is in only one chromosome pair and in an internal position of the long arm (q arm), because this organization has been observed in 3 out of the 5 species. Finally, the major ribosomal probe shows the most conservative pattern, which also matches with the plesiomorphic condition in the fish group, i.e., in an internal position and near the centromere [10].
              Table 2

              Localization summary of the five probes used in Halobatrachus didactylus and other Batrachoididae species





              5S rDNA

              18S rRNA

              U2 snRNA

              Gene co-localization

              Halobatrachus didactylus (present work)



              1 SMT (qSc)

              1 SMT (qSc)

              1 SMT (pTl)

              1 SMT (qI)


              Batrachoides manglae[8]



              Dispersed (++)

              2 A (1 qStl and 1 qSc)

              1 STL (pTl)

              Dispersed (+++) 1 A (qStl)


              Porichthys plectrodon[8]



              Dispersed (+++)

              1 MT (qI)

              1 STL (pTl)

              Dispersed (+)


              Thalassophryne maculosa[8]



              Dispersed (+++)

              1 STL (qTl)

              1 STL (pTl)

              Dispersed (+) 1 STL (qI)

              18S rRNA and U2 snRNA

              Amphichthys cryptocentrus[8]



              Dispersed (++)

              1 STL (qSc)

              1 STL (pTl)

              Dispersed (+)

              18S rRNA and 5S rDNA

              +: low abundance; ++: medium abundance; +++: high abundance; A: acrocentric; STL: subtelocentric; SMT: submetacentric; MT: metacentric; q: long arm; p: short arm; Sc: subcentromeric; I: internal; Stl; subtelomeric; Tl: telomeric.

              Repetitive DNAs, like GATA repeats, are sequences with an important role in the structural and functional evolution of genomes [60]. In particular, the (GATA)n sequence has been implicated in the sex chromosome development of several organisms, including humans [61], snakes [62], guppy fish [63] and wolf fish [64]. It is probable that the sex-determining genes arose from a mutation of a duplicated gene, resulting in a novel sex-determining gene and chromosome [65]. Thus, it cannot be discounted that the GATA clustering in a chromosome pair observed in H. didactylus, could be associated with a proto-sex chromosome pair.


              For 5S rDNA, a mixed evolutionary model has been proposed, in which the presence of two different types that group in a between-species pattern reveals a birth-and-death model, but the high homogeneity of the sequences in each type indicates a homogeneization derived by the concerted model. Moreover, the existence of the birth-and-death process could not be originated by a duplication of a pre-existing 5S rDNA; instead, an ancient HT event between Batrachoidiformes and Perciformes could have introduced the new type. The 5S rDNA provides the requirements for a HT because of its recognized ability to transpose to other loci [66]. Therefore future research on this type of transfer should also be directed towards ascertaining the specific mechanism by which the transfer took place.

              The U2 snRNA gene was not co-located with the other repetitive sequence, so four different chromosomal markers have been obtained in H. didactylus, which is important for producing a preliminary physical map of the species. Moreover, cytogenetic data have shown that the U2 snRNA gene have experienced a transition from a dispersed organization to a clustered organization during the Batrachoididae evolution. A similar case has been detected for (GATA)n repeats. All these characteristics make the Batrachoididae family an interesting group for studying the multigene families from the evolutionary perspective.


              Sampling, PCR amplification, cloning and sequencing

              Up to six specimens (3 males and 3 females) of H. didactylus were collected from natural populations in the Bay of Cadiz (SW Spain) for molecular purposes. Meanwhile, cytogenetic analysis was made from three specimens. Sampling procedures comply with the ARRIVE guidelines and with the UCACG052009 Regulation of the Bioethical Committee from the University of Cádiz.

              Genomic DNA was isolated from muscle tissue of the six specimens using FastDNA kit® (Q-Biogene). Extraction quality was validated by electrophoresis in agarose gel (1.5%) containing 0.5 μg ml-1 ethidium bromide. Three individuals were used to amplify the ITS-1 and U2 snRNA multigene families, while four were amplified for 5S rDNA (Table 3). These PCR amplifications were made using the following primers: A and B described by [67] for 5S rDNA amplification; Sp18 and Sp58 described by [68] for 18S rRNA-ITS-1-5.8S rRNA amplification; and U2ang-Fwd and U2ang-Rev described by [69] for U2 snDNA amplification. An internal reverse primer specific to the NTS of the 5S rDNA type β (5SBETAINT: 5’-CGAGGTGCAGGAGACAAACAC-3’) was designed to corroborate the amplification in the following fish species: H. didactylus, S. aurata (positive controls), four species of the Batrachoididae family (Porichthys plectrodon, Thalassophryne maculosa, Amphichthys cryptocentrus and Batrachoides manglae), two species of the Sparidae family (Pagrus pagrus and Diplodus sargus), Dicentrarchus labrax (Moronidae family), Argyrosomus regius (Sciaenidae family), Plectorhinchus mediterraneus (Haemulidae family), Solea senegalensis (Soleidae family) and Engraulis encrasicholus (Engraulidae family). Parallel to this experiment another PCR was carried out with the ITS-1 primers in the same species, as positive control for genomic DNA. Reactions were carried out in a final volume of 50 μl containing 60–80 ng of genomic DNA, 3 mM Cl2Mg, 300 μM dNTP, 0.2 pmol of the forward and reverse primers and 3 U of Taq polymerase (Euroclone). The PCR conditions were according to [70] and were performed in a Gene Amp® PCR System 2700 (Applied Biosystems) thermal cycler.
              Table 3

              Specimens used, number of clones sequenced (N), and GenBank accession number




              GenBank acc. no.

              5S rDNA

              ♂0 (type α)


              JN406335- JN406337; JN406339- JN406343

              ♂0 (type β)


              ♀0 (type α)


              JN406320- JN406326; JN406338

              ♀0 (type β)


              ♂1 (type α)


              JN406344- JN406352

              ♂1 (type β)


              ♀3 (type α)


              JN406327- JN406334

              ♀3 (type β)





              JN406365- JN406366; JN406370



              JN406367- JN406369



              JN406362- JN406364

              U2 snRNA



              JN406356- JN406358



              JN406353- JN406355



              JN406359- JN406361

              The specimens are denoted by a female or male symbol followed by a number.

              The PCR products were purified using the NucleoSpin® Extract II kit (Macherey-Nagel), and cloned into pGEM®-T Easy Vector (Promega). Plasmid DNA was extracted using NucleoSpin® Plasmid (Macherey-Nagel). DNA sequencing was performed with fluorescence-labeled terminator (BigDye Terminator 3.1 Cycle Sequencing Kit; Applied Biosystems) in an ABI3100 Genetic Analyzer.

              Sequence analysis

              From 7 to 9 clones per specimen were sequenced in 5S rDNA, and 3 clones per specimen in both ITS-1 and U2 snRNA gene (Table 3). Sequence data have been deposited with the GenBank Data Library under Accession Numbers JN406320 to JN406370.

              The sequences were aligned in MAFFT [71] using the slow and iterative refinement method strategy (FFT-NS-i). Consensus sequences were obtained with the Bioedit software [72]. A sequence-similarity search of coding and spacers sequences from the different multigene families was performed in BLASTN [73] to determine the similarities of the sequences obtained with other sequences from the GenBank database. Promoter and regulator regions of the 5S gene were detected by comparison with models previously described [20]. Regarding the U2 snRNA gene, the PSE and 3’ box elements were detected by aligning, respectively, the 5’ end and the 3’ end (respect to the coding region) with those obtained in the Moronidae fish family [51]. The putative DSE region was detected also comparing with those of the Moronidae family and with the consensus octameric sequence described for this region (ATGCAAAT) [17]. DnaSP version 5 program [74] was used to obtain the nucleotide variability (π), the number of polymorphic sites (s) and the number of haplotypes (h) in both coding and spacer regions of the three gene families. On the other hand, evolutionary distance between 5S rDNA types was calculated, using MEGA5 program [75], by the number of base substitutions per site from averaging over all sequence pairs, using the Kimura-2-parameter [76] in relation to evolutionary distances (d). Standard error estimates were obtained by a bootstrap procedure (1000 replicates), and complete deletion option was also applied.

              The dimers obtained from 5S rDNA type α were used to extract the complete sequence of the coding region, and this sequence was subjected to the RNAstructure 5.2 program [77]. A consensus sequence from 5S rDNA type β was also obtained to be folded in RNAstructure 5.2 program. The secondary structures obtained were compared with those predicted in previous described models [16].

              Phylogenetic analysis

              A phylogenetic analysis was carried out with the NTS sequences of H. didactylus and those of other four species from the Batrachoididae family, which were extracted from GenBank database (from GU645582 to GU645701). Six sequences of two different species of the Sparidae family were also included in the analysis, two from Pagrus pagrus (Acc. No. HM004380 and HM004384), two from Pagrus auriga (Acc. No. HM004389 and HM004390) and two from Sparus aurata (Acc. No. AY330701 and AY330702). In addition, two sequences of Dicentrarchus labrax were rooted as outgroup (Acc. No. HM014364 and HM014367). A Neighbor-Joining (NJ) [77] tree was firstly inferred in MEGA5 [75]. The Branch support was tested by the bootstrap method [78]. The evolutionary distances were computed using the p-distance method [79]. This NJ tree was used thereafter as the starting tree for the inference of a Maximum-Likelihood tree in the PhyML 3.0 program [80]. Previously, a statistical selection of best-fit models of nucleotide substitution had been performed by using the jModelTest program, version 0.1.1 [81], using the Akaike Information Criterion (AIC) for model selection. The GTR model [82] was chosen as the best for the PhyML 3.0 program (−lnL= 12156.5732, AIC= 25252.6211). The branch support of the ML tree was measured by applying the non-parametric version of the approximate likelihood ratio test (SH-aLRT) [83], and the option SPR was also chosen for tree improvement.

              Cytogenetic techniques

              Chromosome preparations were made from cephalic kidney as described by [2]. Slides were pretreated with RNase, pepsin and formaldehyde according to [84]. Finally, the samples were dehydrated in successive steps using ethanol and stored at −80°C up to the moment of hybridization.

              5S rDNA and U2 snRNA probes were labeled by the PCR-based method using PCR primers as previously mentioned. Similarly, 18S rDNA probe was labeled by PCR-based method using the primers 18S F and 28S R described by [85]. The (GATA)n probe was obtained by PCR as was described by [86], using the primers (GATA)7 and (TATC)7. The labelling of GATA probe was performed by Nick Translation procedure according to the manufacturer’s instructions.

              Double-FISH techniques were performed according to [10]. Finally, better images were obtained with an epifluorescence microscope (Axioskop 2 Plus, Zeiss), equipped with a cooled camera (CoolSnap, Photometrics© Inc.).



              This work was supported by grants from the Junta de Andalucía (Spain) to the PAI BIO-219 group (LR), CACYTMAR and the “Proyecto de excelencia” RNM-03074 (CS). MUM holds a CSIC JAE fellowship.

              Authors’ Affiliations

              Laboratorio Genética, Facultad de Ciencias del Mar y Ambientales, CACYTMAR, Universidad de Cádiz
              Instituto de Ciencias Marinas de Andalucía – CSIC, Polígono Río San Pedro
              Instituto de Investigaciones Cientificas, Universidad de Oriente


              1. Greenfield DW: Allenbatrachus, a new genus of Indo-Pacific toadfish (Batrachoididae). Pac Sci 1997, 51:306–313.
              2. Palazón JL, Nirchio M, Sarasquete C: Conventional karyotype and nucleolar organiser regions of the toadfish Halobatrachus didactylus (Schneider, 1801) (Pisces: Batrachoididae). Sci 2003, 67:445–449.
              3. Roux C: Batrachoididae. In Fishes of North-eastern Atlantic and Mediterranean. 3rd edition. Edited by: Whitehead PJP, Bauchot ML, Hureau JC, Nielsen JL, Tortonese E. Paris: UNESCO; 1986:1360–1361.
              4. Sarasquete C, Gutiérrez M, Establier R: Efecto del mercurio inorgánico (HgCl2) sobre la sangre y riñón del pez sapo marino, Halobatrachus didactylus (Schneider, 1801). Inv Pesq 1982, 46:323–330.
              5. Desantis S, Cirillo F, Deflorio M, Megalofonou P, Palazón JL, Sarasquete C, de Metrio G: Histochemical study of glycoconjugates in the toadfish, Halobatrachus didactylus oesophagus epithelium. Histol Histopathol 2007, 22:23–35.PubMed
              6. Palazón JL, Arias AM, Sarasquete C: Aspects of the reproductive biology of the toadfish, Halobatrachus didactylus (Schneider, 1801) (Pisces: Batrachoididae). Sci 2001, 65:131–138.
              7. Sarasquete C: Estudios de las series eritrocítica, leucocítica y trombocítica de la sangre periférica del pez sapo marino, Halobatrachus didactylus. Inv Pesq 1983, 47:403–412.
              8. Úbeda-Manzanaro M, Merlo MA, Palazón JL, Cross I, Sarasquete C, Rebordinos L: Chromosomal mapping of the major and minor ribosomal genes, (GATA)n and U2 snRNA gene by double-colour FISH in species of the Batrachoididae family. Genetica 2010, 138:787–794.PubMedView Article
              9. Úbeda-Manzanaro M, Merlo MA, Palazón JL, Sarasquete C, Rebordinos L: Sequence characterization and phylogenetic analysis of the 5S ribosomal DNA in species of the family Batrachoididae. Genome 2010, 53:723–730.PubMedView Article
              10. Merlo A, Cross I, Palazón JL, Sarasquete C, Rebordinos L: Chromosomal mapping of the major and minor ribosomal genes, (GATA)n and (TTAGGG)n by one-color and double-color FISH reveals three chromosomal markers in the toadfish Halobatrachus didactylus (Teleostei: Batrachoididae). Genetica 2007, 131:195–200.PubMedView Article
              11. Nei M, Rooney AP: Concerted and Birth-and-Death evolution of multigene families. Annu Rev Genet 2005, 39:121–152.PubMedView Article
              12. Manchado M, Zuasti E, Cross I, Merlo A, Infante C, Rebordinos L: Molecular characterization and chromosomal mapping of the 5S rRNA gene in Solea senegalensis: a new linkage to the U1, U2, and U5 small nuclear RNA genes. Genome 2006, 49:79–86.PubMedView Article
              13. Pinhal D, Araki CS, Gadig OBF, Martins C: Molecular organization of 5S rDNA in sharks of the genus Rhizoprionodon: insights into the evolutionary dynamics of 5S rDNA in vertebrate genomes. Genet Res, Camb 2009, 91:61–72.View Article
              14. Vierna J, González-Tizón AM, Martínez-Lage A: Long-term evolution of 5S ribosomal DNA seems to be driven by Birth-and-Death processes and selection in Ensis razor shells (Mollusca: Bivalvia). Biochem Genet 2009, 47:635–644.PubMedView Article
              15. McNamara-Schroeder KJ, Hennessey RF, Harding GA, Jensen RC, Stumph WE: The Drosophila U1 and U6 gene proximal sequence elements act as important determinants of the RNA polymerase specificity of small nuclear RNA gene promoters in vitro and in vivo. J Biol Chem 2001, 276:31786–31792.PubMedView Article
              16. Szymanski M, Barciszewska MZ, Erdmann VA, Barciszewski J: 5S ribosomal RNA database. Nucleic Acids Res 2002, 30:176–178.PubMedView Article
              17. Tebb G, Mattaj IW: The Xenopus laevis U2 gene Distal Sequence Element (enhancer) is composed of four subdomains that can act independently and are partly functionally redundant. Mol Cel Biol 1989, 9:1682–1690.
              18. Sola L, Rossi AR, Annesi F, Gornung E: Cytogenetic studies in Sparus auratus (Pisces, Perciformes): molecular organization of 5S rDNA and chromosomal mapping of 5S and 45S ribosomal genes and of telomeric repeats. Hereditas 2003, 139:232–236.PubMedView Article
              19. Zhu HP, Lu MX, Gao FY, Huang ZH, Yang LP, Gui JF: Chromosomal localization of rDNA genes and genomic organization of 5S rDNA in Oreochromis mossambicus, O. urolepis hornorum and their hybrid. J Genet 2010, 89:163–171.PubMedView Article
              20. Martins C, Wasko AP: Organization and evolution of 5S ribosomal DNA in the fish genome. In Focus on genome research. Edited by: Williams CR. New York: Nova Science Publishers; 2004:335–363.
              21. Ocalewicz K, Hliwa P, Krol J, Rábová M, Stabinski R, Ráb P: Karyotype and chromosomal characteristics of Ag-NOR sites and 5S rDNA in European smelt, Osmerus eperlanus. Genetica 2007, 131:29–35.PubMedView Article
              22. Freire R, Arias A, Insua A, Méndez J, Eirín-López JM: Evolutionary dynamics of the 5S rDNA gene family in the mussel Mytilus: mixed effects of birth-and-death and concerted evolution. J Mol Evol 2010, 70:413–426.PubMedView Article
              23. Fujiwara M, Inafuku J, Takeda A, Watanabe A, Fujiwara A, Kohno S, Kubota S: Molecular organization of 5S rDNA in bitterlings (Cyprinidae). Genetica 2009, 135:355–365.PubMedView Article
              24. Finn RN, Kristoffersen BA: Vertebrate vitellogenin gene duplication in relation to the "3R Hypothesis": Correlation to the Pelagic egg and the Oceanic radiation of Teleosts. PLoS One 2007, 2:e169.PubMedView Article
              25. Frost LS, Leplae R, Summers AO, Toussaint A: Mobile genetic elements: the agents of open source evolution. Nat Rev Microbiol 2005, 3:722–732.PubMedView Article
              26. Treangen TJ, Rocha EP: Horizontal Transfer, not duplication, drives the expansion of protein families in prokaryotes. PLoS Genet 2011, 7:e1001284.PubMedView Article
              27. Keeling PJ, Palmer JD: Horizontal gene transfer in eukaryotic evolution. Nat Rev Genet 2008, 9:605–618.PubMedView Article
              28. Schaack S, Gilbert C, Feschotte C: Promiscuous DNA: horizontal transfer of transposable elements and why it matters for eukaryotic evolution. Trends Ecol Evol 2010, 25:537–546.PubMedView Article
              29. Kidwell MG: Horizontal transfer of P-elements and other short inverted repeat transposons. Genetica 1992, 86:275–286.PubMedView Article
              30. Deprá M, Panzera Y, Ludwig A, Valente VLS, Loreto ELS: hosimary: a new hAT transposon group involved in horizontal transfer. Mol Genet Genomics 2010, 283:451–459.PubMedView Article
              31. Klasson L, Kambris Z, Cook PE, Walker T, Sinkins SP: Horizontal gene transfer between Wolbachia and the mosquito Aedes aegypti. BMC Genomics 2009, 10:33–41.PubMedView Article
              32. Pace JK, Gilbert C, Clark MS, Feschotte C: Repeated horizontal transfer of a DNA transposon in mammals and other tetrapods. Proc Natl Acad Sci USA 2008, 105:17023–17028.PubMedView Article
              33. Roulin A, Piegu B, Wing RA, Panaud O: Evidence of multiple horizontal transfers of the long terminal repeat retrotransposon RIRE1 within the genus Oryza. Plant J 2008, 53:950–959.PubMedView Article
              34. de Boer JG, Yazawa R, Davidson WS, Koop BF: Bursts and horizontal evolution of DNA transposons in the speciation of pseudotetraploid salmonids. BMC Genomics 2007, 8:422–431.PubMedView Article
              35. Morgante M, Brunner S, Pea G, Fengler K, Zuccolo A, Rafalski A: Gene duplication and exon shuffling by helitron-like transposons generate intraspecies diversity in maize. Nat Genet 2005, 37:997–1002.PubMedView Article
              36. Kapitonov VV, Jurka J: A novel class of SINE elements derived from 5S rRNA. Mol Biol Evol 2003, 20:694–702.PubMedView Article
              37. Rooney AP, Ward TJ: Evolution of a large ribosomal RNA multigene family in filamentous fungi: Birth and death of a concerted evolution paradigm. PNAS 2005, 102:5084–5089.PubMedView Article
              38. Piskurek O, Okada N: Poxviruses as possible vectors for horizontal transfer of retrotransposons from reptiles to mammals. Proc Natl Acad Sci USA 2007, 104:12046–12051.PubMedView Article
              39. Raychoudhury R, Baldo L, Oliveira DCSG, Werren JH: Modes of acquisition of Wolbachia: Horizontal transfer, hybrid introgression, and codivergence in the Nasonia species complex. Evolution 2009, 63:165–183.PubMedView Article
              40. Gilbert C, Schaack S, Pace JK, Brindley PJ, Feschotte C: A role for host-parasite interactions in the horizontal transfer of transposons across phyla. Nature 2010, 464:1347–1350.PubMedView Article
              41. Jehle JA, Nickel A, Vlak JM, Backhaus H: Horizontal escape of the novel Tc1-like lepidopteran transposon TCp3.2 into Cydia pomonella granulovirus. J Mol Evol 1998, 46:215–224.PubMedView Article
              42. Anderson MT, Seifert HS: Opportunity and means: horizontal gene transfer from the human host to a bacterial pathogen. mBio 2011, 2:e00005-e00011.PubMedView Article
              43. Spadafora C: Mini-Review: developments in reproductive biology and medicine. Hum Reprod 2008, 23:735–740.PubMedView Article
              44. Bacci ML, Zannoni A, de Cecco M, Fantinati P, Bernardini C, Galeati G, Spinaci M, Giovannoni R, Lavitrano M, Seren E, et al.: Sperm-mediated gene transfer-treated spermatozoa maintain good quality parameters and in vitro fertilization ability in swine. Theriogenology 2009, 72:1163–1170.PubMedView Article
              45. Lavitrano M, Forni M, Bacci ML, di Stefano C, Varzi V, Wang H, Seren E: Sperm mediated gene transfer in pig: selection of donor boars and optimization of DNA uptake. Mol Rep Dev 2003, 64:284–291.View Article
              46. Khoo H-W: Sperm-mediated gene transfer studies on zebrafish in Singapore. Mol Rep Dev 2000, 56:278–280.View Article
              47. Smith KR: The role of sperm-mediated gene transfer in genome mutation and evolution. Med Hypotheses 2002, 59:433–437.PubMedView Article
              48. Dell'Anno A, Corinaldesi C: Degradation and turnover of extracellular DNA in marine sediments: Ecological and methodological considerations. Appl Environ Microbiol 2004, 70:4384–4386.PubMedView Article
              49. Liao D, Pavelitz T, Kidd JR, Kidd KK, Weiner AM: Concerted evolution of the tandemly repeated genes encoding human U2 snRNA (The RNU2 locus) involves rapid intrachromosomal homogenization and rare interchromosomal gene conversion. EMBO J 1997, 16:588–598.PubMedView Article
              50. Pavelitz T, Liao D, Weiner AM: Concerted evolution of the tandem array enconding primate U2 snRNA (the RNU2 locus) is accompained by dramatic remodeling of the junctions with flanking chromosomal sequences. EMBO J 1999, 18:3783–3792.PubMedView Article
              51. Merlo MA, Cross I, Chairi H, Manchado M, Rebordinos L: Analysis of three multigene families as useful tools in species characterization of two closely-related species, Dicentrarchus labrax, Dicentrarchus punctatus and their hybrids. Genes Genet Syst 2010, 85:341–349.PubMedView Article
              52. Merlo MA, Cross I, Rodríguez-Rúa A, Manchado M, Rebordinos L: First approach to studying the genetics of the meagre (Argyrosomus regius; Asso, 1801) using three multigene families. Aquaculture Res 2012.
              53. Manchado M, Rebordinos L, Infante C: U1 and U2 Small Nuclear RNA Genetic Linkage: A Novel Molecular Tool for Identification of Six Sole Species (Soleidae, Pleuronectiformes). J Agric Food Chem 2006, 54:3765–3767.PubMedView Article
              54. Cuello P, Boyd DB, Dye MJ, Proudfoot NJ, Murphy S: Transcription of the human U2 snRNA genes continues beyond the 3' box in vivo. EMBO J 1999, 18:2867–2877.PubMedView Article
              55. Parry HD, Tebb G, Mattaj IW: The Xenopus U2 gene PSE is a single, compact, element required for transcription initiation and 3' end formation. Nucleic Acids Res 1989, 17:3633–3644.PubMedView Article
              56. Chow S, Ueno Y, Toyokawa M, Oohara I, Takeyama H: Preliminary analysis of length and GC content variation in the ribosomal first Internal Transcribed Spacer (ITS1) of marine animals. Mar Biotechnol 2009, 11:301–306.PubMedView Article
              57. Pleyte KA, Duncan SD, Phillips RB: Evolutionary Relationships of the salmonid fish genus Salvelinus inferred fron DNA sequences of the first Internal Transcribed Spacer (ITS 1) of ribosomal DNA. Mol Phylogenet Evol 1992, 1:223–230.PubMedView Article
              58. Gornung E, Colangelo P, Annesi F: 5S ribosomal RNA genes in six species of Mediterranean grey mullets: genomic organization and phylogenetic inference. Genome 2007, 50:787–795.PubMedView Article
              59. Nirchio M, Turner BJ, Pérez JE, Gaviria JI, Cequea H: Karyotypes of three species of toadfish (Batrachoididae: Teleostei) from Venezuela. Sci 2002, 66:1–4.
              60. Matsubara K, Tarui H, Toriba M, Yamada K, Nishida-Umehara C, Agata K, Matsuda Y: Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes. PNAS 2006, 103:18190–18195.PubMedView Article
              61. Subramanian S, Mishra RK, Singh L: Genome-wide analysis of Bkm sequences (GATA repeats): predominant association with sex chromosomes and potential role in higher order chromatin organization and function. Bioinformatics 2003, 19:681–685.PubMedView Article
              62. Epplen JT, McCarrey JR, Sutou S, Olmo S: Base sequence of a cloned snake W-chromosome DNA fragment and identification of a male-specific putative mRNA in the mouse. Proc Natl Acad Sci USA 1982, 79:3798–3802.PubMedView Article
              63. Nanda I, Feichtinger W, Schmid M, Schroeder JH, Zischler H, Epplen JC: Simple repetitive sequences are associated with differentiation of the sex chromosomes in the guppy fish. J Mol Evol 1990, 30:456–462.View Article
              64. Cioffi MB, Martins C, Vicari MR, Rebordinos L, Bertollo LAC: Differentiation of the XY sex chromosomes in the fish Hoplias malabaricus (Characiformes, Erythrinidae): Unusual accumulation of repetitive sequences on the X chromosome. Sex Dev 2010, 4:176–185.PubMedView Article
              65. Cioffi MB, Molina WF, Moreira-Filho O, Bertollo LAC: Chromosomal distribution of repetitive DNA sequences highlights the independent differentiation of multiple sex chromosomes in two closely related fish species. Cytogenet Genome Res 2011, 134:295–302.PubMedView Article
              66. Drouin G, Moniz de Sá M: The concerted evolution of 5S ribosomal genes linked to the repeat units of other multigene families. Mol Biol Evol 1995, 12:481–493.PubMed
              67. Pendás AM, Morán P, Freije JP, García-Vazquez E: Chromosomal mapping and nucleotide sequence of two tandem repeats of Atlantic salmon 5S rDNA. Cytogenet Cell Genet 1994, 67:31–36.PubMedView Article
              68. Kuriiwa K, Hanzawa N, Yoshino T, Kimura S, Nishida M: Phylogenetic relationships and natural hybridization in rabbitfishes (Teleostei: Siganidae) inferred from mitochondrial and nuclear DNA analyses. Mol Phylogenet Evol 2007, 45:69–80.PubMedView Article
              69. Cross I: Caracterización poblacional y citogenética de Crassostrea angulata. Análisis molecular del ADNr 5S en ostreidos. Cadiz University: Department of Biomedicine, Biotechnology and Public Health; 2005. [PhD Thesis]
              70. Cross I, Díaz E, Sanchez I, Rebordinos L: Molecular and cytogenetic characterization of Crassostrea angulata chromosomes. Aquaculture 2005, 247:135–144.View Article
              71. Katoh K, Misawa K, Kuma K-I, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res 2002, 30:3059–3066.PubMedView Article
              72. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 1999, 41:95–98.
              73. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215:403–410.PubMed
              74. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25:1451–1452.PubMedView Article
              75. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol 2011, 28:2731–2739.PubMedView Article
              76. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980, 16:111–120.PubMedView Article
              77. Reuter JS, Mathews DH: RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinforma 2010, 11:129–137.View Article
              78. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987, 4:406–425.PubMed
              79. Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution 1985, 39:783–791.View Article
              80. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol 2010, 59:307–321.PubMedView Article
              81. Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol 2008, 25:1253–1256.PubMedView Article
              82. Lanave C, Preparata G, Saccone C, Serio G: A new method for calculating evolutionary substitution rates. J Mol Evol 1984, 20:86–93.PubMedView Article
              83. Anisimova M, Gil M, Dufayard J-F, Dessimoz C, Gascuel O: Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol 2011, 60:685–699.PubMedView Article
              84. Weigant J, Ried T, Nederlof PM, van der Ploeg M, Tanke HJ, Raap AK: In situ hybridization with fluoresceinated DNA. Nucleic Acids Res 1991, 19:3237–3241.View Article
              85. Jansen G, Devaere S, Weekers PHH, Adriaens D: Phylogenetic relationships and divergence time estimate of African anguilliform catfish (Siluriformes: Clariidae) inferred from ribosomal gene and spacer sequences. Mol Phylogenet Evol 2006, 38:65–78.PubMedView Article
              86. Ijdo JW, Wells RA, Baldini A, Reeders ST: Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCR. Nucleic Acids Res 1991, 19:4780.PubMedView Article


              © Merlo et al.; licensee BioMed Central Ltd. 2012

              This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.