Document Detail

Evidence for 5S rDNA horizontal transfer in the toadfish Halobatrachus didactylus (Schneider, 1801) based on the analysis of three multigene families.
Jump to Full Text
MedLine Citation:
PMID:  23039906     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
BACKGROUND: The Batrachoididae family is a group of marine teleosts that includes several species with more complicated physiological characteristics, such as their excretory, reproductive, cardiovascular and respiratory systems. Previous studies of the 5S rDNA gene family carried out in four species from the Western Atlantic showed two types of this gene in two species but only one in the other two, under processes of concerted evolution and birth-and-death evolution with purifying selection. Here we present results of the 5S rDNA and another two gene families in Halobatrachus didactylus, an Eastern Atlantic species, and draw evolutionary inferences regarding the gene families. In addition we have also mapped the genes on the chromosomes by two-colour fluorescence in situ hybridization (FISH).
RESULTS: Two types of 5S rDNA were observed, named type α and type β. Molecular analysis of the 5S rDNA indicates that H. didactylus does not share the non-transcribed spacer (NTS) sequences with four other species of the family; therefore, it must have evolved in isolation. Amplification with the type β specific primers amplified a specific band in 9 specimens of H. didactylus and two of Sparus aurata. Both types showed regulatory regions and a secondary structure which mark them as functional genes. However, the U2 snRNA gene and the ITS-1 sequence showed one electrophoretic band and with one type of sequence. The U2 snRNA sequence was the most variable of the three multigene families studied. Results from two-colour FISH showed no co-localization of the gene coding from three multigene families and provided the first map of the chromosomes of the species.
CONCLUSIONS: A highly significant finding was observed in the analysis of the 5S rDNA, since two such distant species as H. didactylus and Sparus aurata share a 5S rDNA type. This 5S rDNA type has been detected in other species belonging to the Batrachoidiformes and Perciformes orders, but not in the Pleuronectiformes and Clupeiformes orders. Two hypotheses have been outlined: one is the possible vertical permanence of the shared type in some fish lineages, and the other is the possibility of a horizontal transference event between ancient species of the Perciformes and Batrachoidiformes orders. This finding opens a new perspective in fish evolution and in the knowledge of the dynamism of the 5S rDNA. Cytogenetic analysis allowed some evolutionary trends to be roughed out, such as the progressive change in the U2 snDNA and the organization of (GATA)n repeats, from dispersed to localized in one locus. The accumulation of (GATA)n repeats in one chromosome pair could be implicated in the evolution of a pair of proto-sex chromosomes. This possibility could situate H. didactylus as the most highly evolved of the Batrachoididae family in terms of sex chromosome biology.
Authors:
Manuel A Merlo; Ismael Cross; José L Palazón; María Ubeda-Manzanaro; Carmen Sarasquete; Laureana Rebordinos
Related Documents :
23417006 - Rapid gene diversification of microcystis cyanophages revealed by long- and short-term ...
19844226 - Computational methods for discovering structural variation with next-generation sequenc...
23919726 - Description and molecular differentiation of a new staphylocystoides (cyclophyllidea: h...
23360976 - A new species and a new record of diatrypaceae from panama.
7828936 - Characterization of levj, a sucrase/fructanase-encoding gene from actinomyces naeslundi...
11978196 - Cloning and characterization of a novel isoform of calpastatin in human adult testis.
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2012-10-07
Journal Detail:
Title:  BMC evolutionary biology     Volume:  12     ISSN:  1471-2148     ISO Abbreviation:  BMC Evol. Biol.     Publication Date:  2012  
Date Detail:
Created Date:  2013-01-15     Completed Date:  2013-06-03     Revised Date:  2013-07-11    
Medline Journal Info:
Nlm Unique ID:  100966975     Medline TA:  BMC Evol Biol     Country:  England    
Other Details:
Languages:  eng     Pagination:  201     Citation Subset:  IM    
Affiliation:
Laboratorio Genética, Facultad de Ciencias del Mar y Ambientales, CACYTMAR, Universidad de Cádiz, Puerto Real (Cádiz), 11510, Spain.
Data Bank Information
Bank Name/Acc. No.:
GENBANK/JN406320;  JN406321;  JN406322;  JN406323;  JN406324;  JN406325;  JN406326;  JN406327;  JN406328;  JN406329;  JN406330;  JN406331;  JN406332;  JN406333;  JN406334;  JN406335;  JN406336;  JN406337;  JN406338;  JN406339;  JN406340;  JN406341;  JN406342;  JN406343;  JN406344;  JN406345;  JN406346;  JN406347;  JN406348;  JN406349;  JN406350;  JN406351;  JN406352;  JN406353;  JN406354;  JN406355;  JN406356;  JN406357;  JN406358;  JN406359;  JN406360;  JN406361;  JN406362;  JN406363;  JN406364;  JN406365;  JN406366;  JN406367;  JN406368;  JN406369;  JN406370
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Animals
Base Sequence
Batrachoidiformes / genetics*
DNA, Ribosomal / chemistry,  classification,  genetics*
DNA, Ribosomal Spacer / chemistry,  classification,  genetics
Electrophoresis, Agar Gel
Gene Transfer, Horizontal*
In Situ Hybridization, Fluorescence
Molecular Sequence Data
Multigene Family*
Perciformes / genetics
Phylogeny
RNA, Ribosomal, 5S / genetics*
RNA, Small Nuclear / genetics
Sequence Analysis, DNA
Sequence Homology, Nucleic Acid
Chemical
Reg. No./Substance:
0/DNA, Ribosomal; 0/DNA, Ribosomal Spacer; 0/RNA, Ribosomal, 5S; 0/RNA, Small Nuclear; 0/U2 small nuclear RNA
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): BMC Evol Biol
Journal ID (iso-abbrev): BMC Evol. Biol
ISSN: 1471-2148
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Merlo et al.; licensee BioMed Central Ltd.
open-access:
Received Day: 29 Month: 3 Year: 2012
Accepted Day: 3 Month: 10 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 7 Month: 10 Year: 2012
Volume: 12First Page: 201 Last Page: 201
PubMed Id: 23039906
ID: 3544641
Publisher Id: 1471-2148-12-201
DOI: 10.1186/1471-2148-12-201

Evidence for 5S rDNA Horizontal Transfer in the toadfish Halobatrachus didactylus (Schneider, 1801) based on the analysis of three multigene families
Manuel A Merlo1 Email: alejandro.merlo@uca.es
Ismael Cross1 Email: ismael.cross@uca.es
José L Palazón23 Email: jluis.palazon@icman.csic.es
María Úbeda-Manzanaro12 Email: maria.ubedamanzanaro@icman.csic.es
Carmen Sarasquete2 Email: carmen.sarasquete@icman.csic.es
Laureana Rebordinos1 Email: laureana.rebordinos@uca.es
1Laboratorio Genética, Facultad de Ciencias del Mar y Ambientales, CACYTMAR, Universidad de Cádiz, Puerto Real (Cádiz), 11510, Spain
2Instituto de Ciencias Marinas de Andalucía – CSIC, Polígono Río San Pedro, Puerto Real (Cádiz), 11510, Spain
3Instituto de Investigaciones Cientificas, Universidad de Oriente, Isla de Margarita, Venezuela

Background

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 [5-7]. 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.


Results
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].

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

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

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.

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.

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.


Discussion

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 [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 [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. 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 [43]. 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 [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].

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.


Conclusions

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.


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

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


Competing interests

The authors declare that they have no competing interests.


Authors’ contributions

MAM carried out the molecular genetic studies, the sequence alignment, the phylogenetic analysis, the cytogenetic techniques, and drafted the manuscript. IC participated in the sequence and phylogenetic analysis and helped to draft the manuscript. JLP, MUM and CS carried out the sampling, tissue extraction, and chromosome preparations. LR conceived and coordinated the study, participated in its design, discussed the results and corrected the manuscript. All authors read and approved the final manuscript.


Acknowledgements

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.


References
Greenfield DW,Allenbatrachus, a new genus of Indo-Pacific toadfish (Batrachoididae)Pac SciYear: 199751306313
Palazón JL,Nirchio M,Sarasquete C,Conventional karyotype and nucleolar organiser regions of the toadfish Halobatrachus didactylus (Schneider, 1801) (Pisces: Batrachoididae)SciYear: 200367445449
Roux C,Whitehead PJP, Bauchot ML, Hureau JC, Nielsen JL, Tortonese EBatrachoididaeFishes of North-eastern Atlantic and MediterraneanYear: 19863Paris: UNESCO13601361
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 PesqYear: 198246323330
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 epitheliumHistol HistopatholYear: 200722233517128408
Palazón JL,Arias AM,Sarasquete C,Aspects of the reproductive biology of the toadfish, Halobatrachus didactylus (Schneider, 1801) (Pisces: Batrachoididae)SciYear: 200165131138
Sarasquete C,Estudios de las series eritrocítica, leucocítica y trombocítica de la sangre periférica del pez sapo marino, Halobatrachus didactylusInv PesqYear: 198347403412
Ú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 familyGeneticaYear: 201013878779410.1007/s10709-010-9460-120440541
Ú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 BatrachoididaeGenomeYear: 20105372373010.1139/G10-04820924421
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)GeneticaYear: 200713119520010.1007/s10709-006-9131-417192824
Nei M,Rooney AP,Concerted and Birth-and-Death evolution of multigene familiesAnnu Rev GenetYear: 20053912115210.1146/annurev.genet.39.073003.11224016285855
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 genesGenomeYear: 200649798610.1139/g05-06816462904
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 genomesGenet Res, CambYear: 200991617210.1017/S001667230800999319220932
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 GenetYear: 20094763564410.1007/s10528-009-9255-119633948
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 vivoJ Biol ChemYear: 2001276317863179210.1074/jbc.M10127320011431466
Szymanski M,Barciszewska MZ,Erdmann VA,Barciszewski J,5S ribosomal RNA databaseNucleic Acids ResYear: 20023017617810.1093/nar/30.1.17611752286
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 redundantMol Cel BiolYear: 1989916821690
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 repeatsHereditasYear: 200313923223615061807
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 hybridJ GenetYear: 20108916317110.1007/s12041-010-0022-120861567
Martins C,Wasko AP,Williams CROrganization and evolution of 5S ribosomal DNA in the fish genomeFocus on genome researchYear: 2004New York: Nova Science Publishers335363
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 eperlanusGeneticaYear: 2007131293510.1007/s10709-006-9110-917061144
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 evolutionJ Mol EvolYear: 20107041342610.1007/s00239-010-9341-320386892
Fujiwara M,Inafuku J,Takeda A,Watanabe A,Fujiwara A,Kohno S,Kubota S,Molecular organization of 5S rDNA in bitterlings (Cyprinidae)GeneticaYear: 200913535536510.1007/s10709-008-9294-218648989
Finn RN,Kristoffersen BA,Vertebrate vitellogenin gene duplication in relation to the "3R Hypothesis": Correlation to the Pelagic egg and the Oceanic radiation of TeleostsPLoS OneYear: 20072e16910.1371/journal.pone.000016917245445
Frost LS,Leplae R,Summers AO,Toussaint A,Mobile genetic elements: the agents of open source evolutionNat Rev MicrobiolYear: 2005372273210.1038/nrmicro123516138100
Treangen TJ,Rocha EP,Horizontal Transfer, not duplication, drives the expansion of protein families in prokaryotesPLoS GenetYear: 20117e100128410.1371/journal.pgen.100128421298028
Keeling PJ,Palmer JD,Horizontal gene transfer in eukaryotic evolutionNat Rev GenetYear: 2008960561810.1038/nrg238618591983
Schaack S,Gilbert C,Feschotte C,Promiscuous DNA: horizontal transfer of transposable elements and why it matters for eukaryotic evolutionTrends Ecol EvolYear: 20102553754610.1016/j.tree.2010.06.00120591532
Kidwell MG,Horizontal transfer of P-elements and other short inverted repeat transposonsGeneticaYear: 19928627528610.1007/BF001337261334912
Deprá M,Panzera Y,Ludwig A,Valente VLS,Loreto ELS,hosimary: a new hAT transposon group involved in horizontal transferMol Genet GenomicsYear: 201028345145910.1007/s00438-010-0531-x20339871
Klasson L,Kambris Z,Cook PE,Walker T,Sinkins SP,Horizontal gene transfer between Wolbachia and the mosquito Aedes aegyptiBMC GenomicsYear: 200910334110.1186/1471-2164-10-3319154594
Pace JK,Gilbert C,Clark MS,Feschotte C,Repeated horizontal transfer of a DNA transposon in mammals and other tetrapodsProc Natl Acad Sci USAYear: 2008105170231702810.1073/pnas.080654810518936483
Roulin A,Piegu B,Wing RA,Panaud O,Evidence of multiple horizontal transfers of the long terminal repeat retrotransposon RIRE1 within the genus OryzaPlant JYear: 20085395095918088314
de Boer JG,Yazawa R,Davidson WS,Koop BF,Bursts and horizontal evolution of DNA transposons in the speciation of pseudotetraploid salmonidsBMC GenomicsYear: 2007842243110.1186/1471-2164-8-42218021408
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 maizeNat GenetYear: 200537997100210.1038/ng161516056225
Kapitonov VV,Jurka J,A novel class of SINE elements derived from 5S rRNAMol Biol EvolYear: 20032069470210.1093/molbev/msg07512679554
Rooney AP,Ward TJ,Evolution of a large ribosomal RNA multigene family in filamentous fungi: Birth and death of a concerted evolution paradigmPNASYear: 20051025084508910.1073/pnas.040968910215784739
Piskurek O,Okada N,Poxviruses as possible vectors for horizontal transfer of retrotransposons from reptiles to mammalsProc Natl Acad Sci USAYear: 2007104120461205110.1073/pnas.070053110417623783
Raychoudhury R,Baldo L,Oliveira DCSG,Werren JH,Modes of acquisition of Wolbachia: Horizontal transfer, hybrid introgression, and codivergence in the Nasonia species complexEvolutionYear: 20096316518310.1111/j.1558-5646.2008.00533.x18826448
Gilbert C,Schaack S,Pace JK,Brindley PJ,Feschotte C,A role for host-parasite interactions in the horizontal transfer of transposons across phylaNatureYear: 20104641347135010.1038/nature0893920428170
Jehle JA,Nickel A,Vlak JM,Backhaus H,Horizontal escape of the novel Tc1-like lepidopteran transposon TCp3.2 into Cydia pomonella granulovirusJ Mol EvolYear: 19984621522410.1007/PL000062969452523
Anderson MT,Seifert HS,Opportunity and means: horizontal gene transfer from the human host to a bacterial pathogenmBioYear: 20112e00005e0001121325040
Spadafora C,Mini-Review: developments in reproductive biology and medicineHum ReprodYear: 20082373574010.1093/humrep/dem42518270183
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 swineTheriogenologyYear: 2009721163117010.1016/j.theriogenology.2009.06.03319767086
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 uptakeMol Rep DevYear: 20036428429110.1002/mrd.10230
Khoo H-W,Sperm-mediated gene transfer studies on zebrafish in SingaporeMol Rep DevYear: 20005627828010.1002/(SICI)1098-2795(200006)56:2+<278::AID-MRD14>3.0.CO;2-A
Smith KR,The role of sperm-mediated gene transfer in genome mutation and evolutionMed HypothesesYear: 20025943343710.1016/S0306-9877(02)00150-012208184
Dell'Anno A,Corinaldesi C,Degradation and turnover of extracellular DNA in marine sediments: Ecological and methodological considerationsAppl Environ MicrobiolYear: 2004704384438610.1128/AEM.70.7.4384-4386.200415240325
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 conversionEMBO JYear: 19971658859810.1093/emboj/16.3.5889034341
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 sequencesEMBO JYear: 1999183783379210.1093/emboj/18.13.378310393193
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 hybridsGenes Genet SystYear: 20108534134910.1266/ggs.85.34121317546
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 familiesAquaculture ResYear: 201210.1111/j.1365-2109.2012.03103.x
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 ChemYear: 2006543765376710.1021/jf053059416719493
Cuello P,Boyd DB,Dye MJ,Proudfoot NJ,Murphy S,Transcription of the human U2 snRNA genes continues beyond the 3' box in vivoEMBO JYear: 1999182867287710.1093/emboj/18.10.286710329632
Parry HD,Tebb G,Mattaj IW,The Xenopus U2 gene PSE is a single, compact, element required for transcription initiation and 3' end formationNucleic Acids ResYear: 1989173633364410.1093/nar/17.10.36332734097
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 animalsMar BiotechnolYear: 20091130130610.1007/s10126-008-9153-218937008
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 DNAMol Phylogenet EvolYear: 1992122323010.1016/1055-7903(92)90019-D1342939
Gornung E,Colangelo P,Annesi F,5S ribosomal RNA genes in six species of Mediterranean grey mullets: genomic organization and phylogenetic inferenceGenomeYear: 20075078779510.1139/G07-05817893718
Nirchio M,Turner BJ,Pérez JE,Gaviria JI,Cequea H,Karyotypes of three species of toadfish (Batrachoididae: Teleostei) from VenezuelaSciYear: 20026614
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 chromosomesPNASYear: 2006103181901819510.1073/pnas.060527410317110446
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 functionBioinformaticsYear: 20031968168510.1093/bioinformatics/btg06712691979
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 mouseProc Natl Acad Sci USAYear: 1982793798380210.1073/pnas.79.12.37986954524
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 fishJ Mol EvolYear: 19903045646210.1007/BF02101117
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 chromosomeSex DevYear: 2010417618510.1159/00030972620502069
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 speciesCytogenet Genome ResYear: 201113429530210.1159/00032948121720155
Drouin G,Moniz de Sá M,The concerted evolution of 5S ribosomal genes linked to the repeat units of other multigene familiesMol Biol EvolYear: 1995124814937739390
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 rDNACytogenet Cell GenetYear: 199467313610.1159/0001337928187548
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 analysesMol Phylogenet EvolYear: 200745698010.1016/j.ympev.2007.04.01817590356
Cross I,Caracterización poblacional y citogenética de Crassostrea angulata. Análisis molecular del ADNr 5S en ostreidos (PhD Thesis).Year: 2005Cadiz University: Department of Biomedicine, Biotechnology and Public Health
Cross I,Díaz E,Sanchez I,Rebordinos L,Molecular and cytogenetic characterization of Crassostrea angulata chromosomesAquacultureYear: 200524713514410.1016/j.aquaculture.2005.02.039
Katoh K,Misawa K,Kuma K-I,Miyata T,MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transformNucleic Acids ResYear: 2002303059306610.1093/nar/gkf43612136088
Hall TA,BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NTNucleic Acids Symp SerYear: 1999419598
Altschul SF,Gish W,Miller W,Myers EW,Lipman DJ,Basic local alignment search toolJ Mol BiolYear: 19902154034102231712
Librado P,Rozas J,DnaSP v5: a software for comprehensive analysis of DNA polymorphism dataBioinformaticsYear: 2009251451145210.1093/bioinformatics/btp18719346325
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 MethodsMol Biol EvolYear: 2011282731273910.1093/molbev/msr12121546353
Kimura M,A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequencesJ Mol EvolYear: 19801611112010.1007/BF017315817463489
Reuter JS,Mathews DH,RNAstructure: software for RNA secondary structure prediction and analysisBMC BioinformaYear: 20101112913710.1186/1471-2105-11-129
Saitou N,Nei M,The neighbor-joining method: a new method for reconstructing phylogenetic treesMol Biol EvolYear: 198744064253447015
Felsenstein J,Confidence limits on phylogenies: an approach using the bootstrapEvolutionYear: 19853978379110.2307/2408678
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.0Syst BiolYear: 20105930732110.1093/sysbio/syq01020525638
Posada D,jModelTest: phylogenetic model averagingMol Biol EvolYear: 2008251253125610.1093/molbev/msn08318397919
Lanave C,Preparata G,Saccone C,Serio G,A new method for calculating evolutionary substitution ratesJ Mol EvolYear: 198420869310.1007/BF021019906429346
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 schemesSyst BiolYear: 20116068569910.1093/sysbio/syr04121540409
Weigant J,Ried T,Nederlof PM,van der Ploeg M,Tanke HJ,Raap AK,In situ hybridization with fluoresceinated DNANucleic Acids ResYear: 1991193237324110.1093/nar/19.12.32372062640
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 sequencesMol Phylogenet EvolYear: 200638657810.1016/j.ympev.2005.09.01116290035
Ijdo JW,Wells RA,Baldini A,Reeders ST,Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCRNucleic Acids ResYear: 199119478010.1093/nar/19.17.47801891373

Figures

[Figure ID: F1]
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).



[Figure ID: F2]
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.



[Figure ID: F3]
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 ID: F4]
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.



[Figure ID: F5]
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.



[Figure ID: F6]
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).



Tables
[TableWrap ID: T1] Table 1 

Polymorphism by gene region


Region Length (bp) %GC S h π
5S rRNA (α)
120
52.50-55.83
12
10
0.018±0.003
NTS (α)
289-304
36.75-39.79
41
14
0.023±0.006
5S rRNA (β)
120
56.67-57.50
3
4
0.004±0.002
NTS (β)
77-87
38.96-43.02
4
5
0.009±0.003
18S rRNA
182
48.90-49.45
3
4
0.004±0.001
ITS-1
388
62.89-63.40
9
5
0.008±0.002
5.8S rRNA
93-94
50.54-51.06
1
2
0.002±0.002
U2 snRNA
189
35.45-38.10
12
7
0.024±0.003
U2 spacer 623-631 36.61-37.88 31 8 0.017±0.002

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


[TableWrap ID: T2] Table 2 

Localization summary of the five probes used inHalobatrachus didactylusand other Batrachoididae species


Species 2n FN (GATA)n 5S rDNA 18S rRNA U2 snRNA Gene co-localization
Halobatrachus didactylus (present work)
46
72
1 SMT (qSc)
1 SMT (qSc)
1 SMT (pTl)
1 SMT (qI)
None
Batrachoides manglae[8]
46
48
Dispersed (++)
2 A (1 qStl and 1 qSc)
1 STL (pTl)
Dispersed (+++) 1 A (qStl)
None
Porichthys plectrodon[8]
44
62
Dispersed (+++)
1 MT (qI)
1 STL (pTl)
Dispersed (+)
None
Thalassophryne maculosa[8]
46
64
Dispersed (+++)
1 STL (qTl)
1 STL (pTl)
Dispersed (+) 1 STL (qI)
18S rRNA and U2 snRNA
Amphichthys cryptocentrus[8] 46 52 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.


[TableWrap ID: T3] Table 3 

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


Sequence Specimen N GenBank acc. no.
5S rDNA
♂0 (type α)
5
JN406335- JN406337; JN406339- JN406343
♂0 (type β)
3
♀0 (type α)
5
JN406320- JN406326; JN406338
♀0 (type β)
3
♂1 (type α)
6
JN406344- JN406352
♂1 (type β)
3
♀3 (type α)
3
JN406327- JN406334
♀3 (type β)
5
ITS-1
♂0
3
JN406365- JN406366; JN406370
♂2
3
JN406367- JN406369
♀4
3
JN406362- JN406364
U2 snRNA ♂0
3
JN406356- JN406358
♀0
3
JN406353- JN406355
♂1 3 JN406359- JN406361

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



Article Categories:
  • Research Article


Previous Document:  Women's views and adherence to low molecular weight heparin therapy during pregnancy and the puerper...
Next Document:  Modelling biological activities of nanoparticles.