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  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  and populations are distributed along the coasts from the Bay of Biscay (Spain) to Ghana as well as the western Mediterranean .
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 , and in hematology, reproduction and histophysiology studies [5–7]. Genetic studies in the Batrachoididae family have been focused mainly on cytogenetics (reviewed by ); there is only one published paper in which this fish family has been studied from a molecular perspective . For H. didactylus in particular, the karyotype and the localization of various repetitive sequences have been described .
Multigene families are comprised of multiple genes that all descend from a common ancestral gene, have similar sequences and are functionally related . 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 . This feature makes the units show more similarity within than between related species . 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 . The histone and the immunoglobulin gene families are typical examples of birth-and-death evolution .
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 . Most of the snRNA are transcribed by RNA polymerase II, except for U6 snRNA which is transcribed by RNA polymerase III . 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 .
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.
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 .
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 .
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).
Polymorphism by gene region
5S rRNA (α)
5S rRNA (β)
%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  (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.
The Maximum Likelihood (ML) tree obtained (Figure 5) was compared with the Neighbor-Joining (NJ) tree previously described by  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 ) 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 , these clones were divided into two groups.
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 . 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.
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 .
In a previous study  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 ), 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  (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 , thus spreading mutation events along the multigene family. Concerted evolution makes the sequence more similar within the species than it is among related species . 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 . 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 (Figure 3). This species and H. didactylus belong to different superorders (Acanthopterygii and Paracanthopterygii respectively), which diverged 55 My ago . 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 ; it has been assumed that between 88% and 98% of the expansion of protein families in prokaryotes is due to HT . However, the occurrence of HT between eukaryote organisms remains obscure . 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 . 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) . Cases have been reported of HTT between a wide range of eukaryotes [30–34]. 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 . 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 , and this interaction might be the cause of the pattern of evolution and the dispersed arrangement of some fungal organisms . 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 , a bacteria  or a parasitic organism . 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 . Few cases of HT from eukaryotes to virus  or bacteria  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 . The first way is the one used most often to ascertain the HT , and similarity can clearly be seen in the H. didactylusS. 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 . This capacity has led to sperm cells being used as vectors for transferring DNA in transgenic experiments [44–46]. 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 . 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 . 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 , primates  and moronid fishes , 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 ; therefore a clustered arrangement could be propitious for a concerted evolution model. Moreover, taking into account the results obtained by  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, Dicentrarchus labrax and D. punctatus, and A. regius 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 . 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  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 . 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 . 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 . 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 . It has been postulated that the unequal cross-over and gene conversion mechanisms occur more frequently in chromosome regions adjacent to the telomeres , 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  and ). 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 . Taking this into account, the most recent species is H. didactylus (FN= 72; ). Comparing the FISH results obtained in the present study with those from the other Batrachoididae species , 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 .
Localization summary of the five probes used inHalobatrachus didactylusand other Batrachoididae species
+: 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 . In particular, the (GATA)n sequence has been implicated in the sex chromosome development of several organisms, including humans , snakes , guppy fish  and wolf fish . 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 . 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 . 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  for 5S rDNA amplification; Sp18 and Sp58 described by  for 18S rRNA-ITS-1-5.8S rRNA amplification; and U2ang-Fwd and U2ang-Rev described by  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  and were performed in a Gene Amp® PCR System 2700 (Applied Biosystems) thermal cycler.
Specimens used, number of clones sequenced (N), and GenBank accession number
GenBank acc. no.
♂0 (type α)
JN406335- JN406337; JN406339- JN406343
♂0 (type β)
♀0 (type α)
JN406320- JN406326; JN406338
♀0 (type β)
♂1 (type α)
♂1 (type β)
♀3 (type α)
♀3 (type β)
JN406365- JN406366; JN406370
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.
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  using the slow and iterative refinement method strategy (FFT-NS-i). Consensus sequences were obtained with the Bioedit software . A sequence-similarity search of coding and spacers sequences from the different multigene families was performed in BLASTN  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 . 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 . 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) . DnaSP version 5 program  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 , by the number of base substitutions per site from averaging over all sequence pairs, using the Kimura-2-parameter  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 . 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 .
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)  tree was firstly inferred in MEGA5 . The Branch support was tested by the bootstrap method . The evolutionary distances were computed using the p-distance method . This NJ tree was used thereafter as the starting tree for the inference of a Maximum-Likelihood tree in the PhyML 3.0 program . Previously, a statistical selection of best-fit models of nucleotide substitution had been performed by using the jModelTest program, version 0.1.1 , using the Akaike Information Criterion (AIC) for model selection. The GTR model  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) , and the option SPR was also chosen for tree improvement.
Chromosome preparations were made from cephalic kidney as described by . Slides were pretreated with RNase, pepsin and formaldehyde according to . 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 . The (GATA)n probe was obtained by PCR as was described by , 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.
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.
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
Greenfield DW: Allenbatrachus, a new genus of Indo-Pacific toadfish (Batrachoididae).Pac Sci 1997, 51:306–313.
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.
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.
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.
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
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.
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.
Ú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
Ú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
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
Nei M, Rooney AP: Concerted and Birth-and-Death evolution of multigene families.Annu Rev Genet 2005, 39:121–152.PubMedView Article
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
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
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
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
Szymanski M, Barciszewska MZ, Erdmann VA, Barciszewski J: 5S ribosomal RNA database.Nucleic Acids Res 2002, 30:176–178.PubMedView Article
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.
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
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
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.
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
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
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
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
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
Treangen TJ, Rocha EP: Horizontal Transfer, not duplication, drives the expansion of protein families in prokaryotes.PLoS Genet 2011, 7:e1001284.PubMedView Article
Keeling PJ, Palmer JD: Horizontal gene transfer in eukaryotic evolution.Nat Rev Genet 2008, 9:605–618.PubMedView Article
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
Kidwell MG: Horizontal transfer of P-elements and other short inverted repeat transposons.Genetica 1992, 86:275–286.PubMedView Article
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
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
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
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
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
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
Kapitonov VV, Jurka J: A novel class of SINE elements derived from 5S rRNA.Mol Biol Evol 2003, 20:694–702.PubMedView Article
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
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
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
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
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
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
Spadafora C: Mini-Review: developments in reproductive biology and medicine.Hum Reprod 2008, 23:735–740.PubMedView Article
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
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
Khoo H-W: Sperm-mediated gene transfer studies on zebrafish in Singapore.Mol Rep Dev 2000, 56:278–280.View Article
Smith KR: The role of sperm-mediated gene transfer in genome mutation and evolution.Med Hypotheses 2002, 59:433–437.PubMedView Article
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
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
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
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
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.
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
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
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
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
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
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
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.
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
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
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
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
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
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
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
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
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
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.
Cross I, Díaz E, Sanchez I, Rebordinos L: Molecular and cytogenetic characterization of Crassostrea angulata chromosomes.Aquaculture 2005, 247:135–144.View Article
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
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.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.J Mol Biol 1990, 215:403–410.PubMed
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data.Bioinformatics 2009, 25:1451–1452.PubMedView Article
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
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
Reuter JS, Mathews DH: RNAstructure: software for RNA secondary structure prediction and analysis.BMC Bioinforma 2010, 11:129–137.View Article
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees.Mol Biol Evol 1987, 4:406–425.PubMed
Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap.Evolution 1985, 39:783–791.View Article
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
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
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
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
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
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
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.