Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-02T22:50:42.520Z Has data issue: false hasContentIssue false

Molecular organization of 5S rDNA in sharks of the genus Rhizoprionodon: insights into the evolutionary dynamics of 5S rDNA in vertebrate genomes

Published online by Cambridge University Press:  17 February 2009

DANILLO PINHAL
Affiliation:
Departamento de Morfologia, Instituto de Biociências, UNESP – Universidade Estadual Paulista, CEP 18618-000, Botucatu, SP, Brazil
CARLOS S. ARAKI
Affiliation:
Departamento de Morfologia, Instituto de Biociências, UNESP – Universidade Estadual Paulista, CEP 18618-000, Botucatu, SP, Brazil
OTTO B. F. GADIG
Affiliation:
Campus Litoral Paulista, UNESP – Universidade Estadual Paulista, São Vicente, SP, Brazil
CESAR MARTINS*
Affiliation:
Departamento de Morfologia, Instituto de Biociências, UNESP – Universidade Estadual Paulista, CEP 18618-000, Botucatu, SP, Brazil
*
*Corresponding author. Tel:/Fax: +55(14)38116264. e-mail: [email protected]

Summary

In this study, we attempted a molecular characterization of the 5S rDNA in two closely related species of carcharhiniform sharks, Rhizoprionodon lalandii and Rhizoprionodon porosus, as well as a further comparative analysis of available data on lampreys, several fish groups and other vertebrates. Our data show that Rhizoprionodon sharks carry two 5S rDNA classes in their genomes: a short repeat class (termed class I) composed of ~185 bp repeats, and a large repeat class (termed class II) arrayed in ~465 bp units. These classes were differentiated by several base substitutions in the 5S coding region and by completely distinct non-transcribed spacers (NTS). In class II, both species showed a similar composition for both the gene coding region and the NTS region. In contrast, class I varied extensively both within and between the two shark species. A comparative analysis of 5S rRNA gene sequences of elasmobranchs and other vertebrates showed that class I is closely related to the bony fishes, whereas the class II gene formed a separate cartilaginous clade. The presence of two variant classes of 5S rDNA in sharks likely maintains the tendency for dual ribosomal classes observed in other fish species. The present data regarding the 5S rDNA organization provide insights into the dynamics and evolution of this multigene family in the fish genome, and they may also be useful in clarifying aspects of vertebrate genome evolution.

Type
Paper
Copyright
Copyright © 2009 Cambridge University Press

1. Introduction

Tandemly arrayed copies of 5S rDNA consist of highly conserved 120 base pair (bp) transcribing sequence, separated from each other by a variable non-transcribed spacer (NTS) (see Long & David, Reference Long and David1980). The high level of conservation of the 5S rRNA genes arises from their essential molecular function of enhancing protein synthesis by stabilizing the ribosome structure (Barciszewska et al., Reference Barciszewska, Szymanski, Erdmann and Barciszewski2001). While the 5S rRNA genes are conserved even among unrelated taxa, their NTSs show extensive variation, and no relationships are usually detectable between the NTSs of closely related species. Thus, the NTSs seem to be subject to rapid evolution, which makes this region both an important resource for studies concerning the organization and evolution of multigene families and genomes and also a marker to trace recent events of evolution. The intense dynamism of 5S rDNA repeats generates variant classes of 5S rDNA, which have been reported in the genomes of several vertebrates, from lampreys to mammals (Komiya et al., Reference Komiya, Hasegawa and Takemura1986; Hallenberg et al., Reference Hallenberg, Nederby-Nielsen and Frederiksen1994; Frederiksen et al., Reference Frederiksen, Cao, Lomholt, Levan and Hallemberg1997). Particularly in marine and freshwater ray-finned fishes, vast numbers of structural and functional data have demonstrated the occurrence of a dual size-class pattern of organization of the 5S rDNA (Moran et al., Reference Móran, Martínez, Garcia-Vásquez and Pendás1996; Céspedes et al., Reference Céspedes, Garcia, Carrera, Gonzalez, Fernandez, Hernandez and Martin1999; Rocco et al., Reference Rocco, Russo, Stingo, Aprea and Odierna1999; Deiana et al., Reference Deiana, Cau, Salvadori, Coluccia, Cannas, Milia and Tagliavini2000; Martins & Galletti, Reference Martins and Galetti2001a; Wasko et al., Reference Wasko, Martins, Wright and Galetti2001; Martins et al., Reference Martins, Wasko, Oliveira, Porto-Foresti, Parise-Maltempi, Wright and Foresti2002; Messias et al., Reference Messias, Ferreira, Wasko, Oliveira, Foresti and Martins2003; Tigano et al., Reference Tigano, Rocco, Ferrito, Costagliola, Pappalardo and Stingo2004; Alves-Costa et al., Reference Alves-Costa, Wasko, Oliveira, Foresti and Martins2006).

The Chondrichthyes are possibly the least-studied fish group with respect to genome structure and evolution. Although genomic studies in fish have increased in the past decade, cartilaginous fishes remain among the least examined in this aspect. For example, studies focused on 5S rDNA have been carried out only on seven species, covering less than 1% of the nearly 1100 living species. These experiments include the characterization of the 5S rRNA gene in one carcharhiniform shark, Scyliorhinus caniculus (Wegnez et al., Reference Wegnez, Denis, Mazabraud and Clerot1978), and more comprehensive analyses of the entire 5S rDNA repeats in six members of Rajiformes: Raja polystigma and Taeniura lymma (Rocco et al., Reference Rocco, Costagliola, Fiorillo, Tinti and Stingo2005); and Raja asterias, Raja clavata, R. polystigma, Raja miraletus and Dipturus oxyrinchus (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006). In the last study, the authors reported the presence of two size classes of 5S rDNA in Rajidae, suggesting that a dual system, such as that detected in bony fishes, could be a common pattern for the 5S rDNA organization in elasmobranchs.

Rhizoprionodon is a genus of the family Carcharhinidae represented worldwide by seven species of small coastal sharks (Compagno, Reference Compagno1984). Since they are abundant and globally distributed from tropical to subtropical oceans, Rhizoprionodon sharks form the basis of important commercial and recreational fisheries as well as substantial artisanal fisheries (Motta et al., Reference Motta, Gadig, Namora and Braga2005). Rhizoprionodon sharks have accounted for about 60% of all sharks landed in Brazil (Motta et al., Reference Motta, Gadig, Namora and Braga2005), and intensive coastal fishing has led these species to be given ‘vulnerable’ status (Rosa et al., Reference Rosa, Gadig, Santos Motta and Namora2004). Several features allow for the discrimination of Rhizoprionodon lalandii from Rhizoprionodon porosus, such as head and nose curvature, dorsal fin length and position relative to pectoral fins, and body colour; however, it is difficult to distinguish these two species from one another once the head and fins are removed at the fisheries (Gadig, Reference Gadig2001). Currently, there are no molecular data on mitochondrial or nuclear markers concerning population or forensic analysis in these two sharpnose sharks. This prevents inspection programs from properly identifying unprocessed products and inhibits these programs' ability to enforce endangered species laws. Moreover, it makes it impossible to estimate catch and trade on a species-specific basis and thus to assess the sustainability and environmental impact of fisheries. To ameliorate these difficulties, increasing efforts have been directed to the discovery of molecular markers that help in the identification of shark species.

Here, in order to advance understanding of the dynamics and evolution of 5S rDNA arrays in Chondrichthyes, and also to contribute to the development of genetic markers that help in the management of sharks, an evaluation of the genetic structure of the 5S rDNA in the genome of R. lalandii and R. porosus was conducted. In addition, we carried out a comparative analysis of our data with the available data on lampreys, elasmobranchs and several model organisms such as humans, chickens, Xenopus and fugu, among other vertebrates, in an attempt to elucidate the evolutionary patterns that guide the 5S rDNA arrays.

2. Material and methods

(i) Animal sampling and DNA isolation

Two shark species, R. lalandii and R. porosus, were collected from different areas of the Brazilian coast, in the western South Atlantic. All the specimens of R. lalandii were from the São Paulo coast in southern Brazil (n=30), whereas those of R. porosus were from distant sampling points, presumably from different populations (São Paulo in southern Brazil, and Sergipe and Ceará in northern Brazil; n=28). Fin clips and gill tissues were collected from fresh or recently frozen animals and immediately immersed in tubes containing pure ethanol. Individuals were identified and voucher specimens were preserved in the fish collection of the Laboratório de Biologia e Genética de Peixes–UNESP, Botucatu, SP, Brazil, under collection numbers LBP-3001 (R. lalandii) and LBP-3155 (R. porosus). Genomic DNA was successfully isolated from shark samples based on the protocol of Aljanabi & Martinez (Reference Aljanabi and Martinez1997), which is an alternative procedure without environmentally hazardous reagents such as phenol and chloroform.

(ii) PCR amplification, cloning and sequencing

The entire 5S rRNA gene and NTS were amplified from the total genomic DNA by PCR. The set of primers 5SA (5′-TAC GCC CGA TCT CGT CCG ATC-3′) and 5SB (5′-CAG GCT GGT ATG GCC GTA AGC-3′), based on the 5S gene sequence of Salmo gardnerii, described by Komiya & Takemura (Reference Komiya and Takemura1979) and applied successfully to other fish species (Martins & Galetti, Reference Martins and Galetti2001a; Wasko et al., Reference Wasko, Martins, Wright and Galetti2001), was used for PCR. The primers 5SA and 5SB were designed to amplify the entire NTS and 118 bp of the 5S rRNA gene. After cloning and sequencing of the PCR fragments obtained, a second set of primers Cart5S2F (5′-TGG GAG ACC GCC TGG GAA-3′) and Cart5S2R (5′-CCA AGT ACT AAC CAG GCC-3′) was designed to amplify the annealing region of the primers 5SA and 5SB, allowing the complete sequencing of the 5S rRNA gene. The PCR reactions had 25 μl of total volume, which contained 2 units of Taq polymerase, 1×Taq buffer, 1·5 mM of MgCl2, 200 μM of dNTPs, 100 pmol of each primer and ~30 ng of genomic DNA. The background reactions were 35 cycles of 1 min at 95°C, 30 s at 55°C and 45 s at 72°C, with a 5-min final extension at 72°C. A negative control was also included to test for any contamination. PCR products were assessed by electrophoresis in 1·25% agarose gels and visualized by ethidium bromide staining and ultraviolet illumination. The amplified DNA fragments were purified using the kit GFX PCR Purification (GE Healthcare) and inserted into pGEM-T plasmids (Promega) that were used to transform Escherichia coli DH5α competent cells (Invitrogen), according to Sambrook & Russel (Reference Sambrook and Russel2001). Positive clones were sequenced on the ABI Prism 3100 automatic DNA sequencer (Perkin-Elmer) with the kit BigDye Terminator Cycle Sequencing (Perkin-Elmer), following the manufacturer's instructions.

(iii) Sequence analysis

Individual sequences were subjected to BLASTn (Altschul et al., Reference Altschul, Gish, Miller, Myers and Lipman1990) searches on the National Center for Biotechnology Information (NCBI) website (http://www.ncbi.nlm.nih.gov/blast). All sequences were aligned using the software DAMBE (Xia & Xie, Reference Xia and Xie2001), using the options: gap open penalty=10 and gap extension penalty=0·1. The final alignment was checked visually and corrected when necessary. The 5S rDNA sequences retrieved from Rhizoprionodon sharks have been deposited in GenBank under the accession numbers FJ517166–FJ517257. Additional molecular data were obtained from different vertebrate taxa available at NCBI database and from the sharks Galeocerdo cuvier and Alopias superciliosus (present paper, GenBank accession numbers FJ539130–FJ539131) and were used in the comparative and evolutionary analysis.

Sequence composition and diversity among clones were examined by means of the program BioEdit 7.0 (Hall, Reference Hall1999). Nucleotide saturation, substitution patterns and genetic distance were examined in MEGA 4 (Tamura et al., Reference Tamura, Dudley, Nei and Kumar2007). Nucleotide saturation was observed by plotting the absolute number of transitions (Ti) and transversions (Tv) against genetic distance values. Maximum-likelihood (ML) phylogenetic relationships among 5S sequences were analysed using the website version of the program PhyML (Guindon & Gascuel, Reference Guindon and Gascuel2003; Guindon et al., Reference Guindon, Lethiec, Duroux and Gascuel2005). ML phylogenetic analyses were carried out using the Tamura–Nei nucleotide substitution model (Tamura & Nei, Reference Tamura and Nei1993), incorporating rate variation (G) and PINVAR with four G-distributed rate classes (Swofford et al., Reference Swofford, Olsen, Waddell and Hillis1996). This model was selected based on a hierarchical hypothesis test of alternative models implemented with Modeltest 3.6 (Posada & Crandall, Reference Posada and Crandall1998). The Ti/Tv ratio, gamma shape parameter and proportion of non-variant sites were estimated by ML from a neighbour-joining tree (BIONJ). Bootstrap resampling (Felsenstein, Reference Felsenstein1985) was applied to assess support for individual nodes, using 500 replicates with random additions and TBR branch swapping. The Bayesian-likelihood method of phylogenetic analysis (Huelsenbeck et al., Reference Huelsenbeck, Ronquist, Nielsen and Bollback2001) was used to evaluate tree topologies of vertebrates through the estimation of probabilities using MrBayes v.3.0 (Ronquist & Huelsenbeck, Reference Ronquist and Huelsenbeck2003). Four chains were run simultaneously for 1 000 000 generations using MrBayes analysis. Every 100th generation was sampled, and the asymptote of likelihood score was detected with the SUMP command. The above procedure was repeated twice. All the sampled topologies beneath the asymptote (12 000 generations in the present study) were discarded from the population of trees considered in the subsequent majority-rule consensus. Consensus trees were produced with the program TreeExplorer implemented in MEGA 4 (Tamura et al., Reference Tamura, Dudley, Nei and Kumar2007).

3. Results

(i) Sequence analysis of 5S rDNA arrays from Rhizoprionodon sharks

The pairs of primers 5SA and 5SB, and Cart5S2F and Cart5S2R allowed the amplification by PCR of the entire 5S rDNA repeats in Rhizoprionodon. Two fragments of roughly 185 and 465 bp (hereafter referred to as class I and class II, respectively) were obtained for both species (Fig. 1). Nucleotide sequencing and BLASTn searches confirmed that the obtained sequences correspond to 5S rDNA repeat units. These units consisted of a 120 bp coding region (5S gene) as well as a NTS of variable length and composition (Fig. 2).

Fig. 1. Agarose gel showing the two 5S rDNA types of R. lalandii (1) and R. porosus (2) amplified by PCR using the primers 5SA and 5SB. M refers to the 1 kb plus molecular weight marker in base pairs (bp).

Fig. 2. Representative sequences of 5S rDNA class I (5S I) and class II (5S II) from R. lalandii (RL) and R. porosus (RP). The start and end points of transcription are indicated by +1 and +120, respectively; the nucleotide substitutions of specific classes are indicated in grey shading, and the ICRs (A box, IE and C box) are indicated in black shading. TATA-like sequences are underlined. SSRs are in boldface. Hyphens represent gaps, and dots indicate identical nucleotides.

Sequencing of 46 recombinant clones showed that each class of 5S rDNA was characterized by distinct NTS and 5S rRNA gene sequences (Fig. 2), which showed variable levels of genetic divergence (Table 1). The GC contents in the 5S genes of class I and class II of R. lalandii were 55·2 and 53·1%, respectively, whereas R. porosus showed rates of 54·7% for class I and 52·5% for class II. The 5S genes of sharks contained some class-specific base substitutions. Regarding R. porosus, the class I and class II coding sequences contained six and two base substitutions, respectively, whereas in R. lalandii, the sequences contained four substitutions for each class. Common alterations detected in the 5S genes of class I in relation to those of class II are listed in Table 2. A transition at position 116 (A→G) caused the 5S class II gene of sharks to lose a HindIII restriction site, and it is extensively conserved among bony fishes and previously detected in the 5S rDNA of Rajidae (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006). In general, these base shifts were detected outside internal control regions (ICRs). When in ICRs, they are present in Box A and/or Box C, but never in the intermediate element (IE) (Fig. 2). Another clear disparity between the 5S rDNA classes was the number of thymidine residues located in the T-rich tail region of the 5S genes: five thymidines in class I and four in class II (Fig. 2). These repeated sequences correspond to transcriptional terminators (Korn & Brown, Reference Korn and Brown1978). TATA-like elements, believed to be regulatory regions for 5S gene transcription, also appear in the NTSs of sharks. These elements were found to be modified to TAAA in NTS I and AATT in NTS II (Fig. 2).

Table 1. Genetic distance of the 5S genes and NTS classes among Rhizoprionodon sharks. NC, number of clones; L, length in base pairs; GD, genetic distance; SE, standard error

Table 2. Diagnostic nucleotide sites present in 5S gene class I and class II from R. lalandii and R. porosus

Phylogenetic analyses of 5S rDNA sequences were carried out between and within each Rhizoprionodon species, analysing the entire 5S rDNA and the 5S gene and NTS sequences separately (Fig. 3ad). Both approaches supported the existence of a 5S rDNA dual pattern in all specimens checked. Additionally, no differences were found among paralogous copies of 5S genes or NTSs from the same 5S rDNA class of Rhizoprionodon sharks, even from different populations. This result was confirmed by the clustering of 5S rDNA classes into two distinct branches on all assembled ML trees (Fig. 3a). Our data also showed that it was not only variability in the NTS that contributed to the separation of 5S rDNA into two classes. In addition, 5S rRNA genes bearing consistent variability also increased the difference between the two 5S rDNA classes, as illustrated by the 96% divergence of relationship trees (Fig. 3d).

Fig. 3. Recovered ML trees for (a) the entire 5S rDNA sequences, (b) 5S gene only, (c) all NTS sequences and (d) NTS II only. Note the existence of two main branches in the trees discriminating the two 5S rDNA classes. In these trees, sequence names correspond to the first two letters of their species names (RL, R. lalandii and RP, R. porosus) and the number of the clones analysed. The scale bar indicates genetic distance based on the Tamura–Nei model analysis. Bootstrap values under 70 were omitted.

Comparing R. lalandii and R. porosus, it is feasible to assert that these sharks carry similar paralogous copies of 5S genes in their respective 5S rDNA classes, given that the sequences are randomly distributed in tree branches of the phylogenetic reconstructions (Fig. 3a, d). In addition, the base-to-base comparison of the 5S genes from both 5S rDNA classes did not permit us to recognize species-specific nucleotide sites.

Regarding NTSs, relationships among the different classes between or within species were only slight because these classes were separated into individual branches on ML trees by a 100% bootstrap value (Fig. 3c). NTS class I had 64–68 bp, and class II had 335–343 bp; mean values of GC content were as follows: class I=49·6% and class II=61·4% in R. porosus; class I=53·4% and class II=62·1% in R. lalandii. Considering all clones, the mean genetic distance of NTS sequences was very different between classes (NTS I: 0·389±0·096; NTS II: 0·025±0·007). While NTS class I was highly variable even within individuals from the same population in a species, large identical sequences were shared by NTS class II in the Rhizoprionodon species (e.g. positions 29–97). Besides this unexpected interspecific similarity, five species-specific sites could be detected (at positions 98, 211, 254, 337 and 338), enough to separate all R. lalandii sequences from those of R. porosus (Table 3). The length of NTS II paralogous sequences varied only slightly in R. porosus, ranging from 335 to 343 bp, whereas all sequences in R. lalandii were 343 bp. The variation of the former species was related to the presence of microsatellites [TCCC]n at positions 251–270 in this species.

Table 3. Diagnostic nucleotide sites present in NTS class I sequences of Rhizoprionodon sharks

(ii) Comparative analysis of 5S rDNA arrays in vertebrata

Coding sequences of 5S rDNA arrays from Rhizoprionodon sharks were aligned and compared with those of a large number of extant vertebrates, providing information about their relationships. Bayesian analysis was used to construct the phylogenetic relationships among paralogous and orthologous 5S rRNA gene sequences of eight elasmobranch species (two orders), several teleost orders, amphibians, reptiles, birds and mammals (see Supplementary Material at http://journals.cambridge.org/grh). In the Bayesian method, Rhizoprionodon sharks clustered according to the 5S rDNA classes on a consensus tree obtained from the analysis of 1000 bootstrap replicates (Fig. 4). Rhizoprionodon 5S rRNA genes of class I clustered close to bony fishes, to Raja type I and to Vertebrata 5S genes, showing it to be the most widespread variety in existence. On the other hand, Rhizoprionodon 5S genes of class II always grouped with other elasmobranch species in an isolated clade, bootstrapping by 91% similarity. Despite the high level of base substitutions in 5S genes of class II (Table 2), we found a higher level of transitions versus transversions (R=2.1) in all sequences, suggesting that homogenizing forces were acting across the transcribed region, thus preserving its function in the ribosome structure.

Fig. 4. A phylogenetic tree of 5S rRNA genes from Elasmobranchii, Teleostei and several other vertebrates using the Bayesian-likelihood method. Bootstrap values under 70 were omitted.

4. Discussion

(i) Double-class structure of 5S rDNA in Rhizoprionodon sharks: is this a rule in sharks?

Although extensive taxon sampling from shark lineages is indispensable to defining the number of 5S rDNA classes, our data provided new insights into the genomic organization of 5S rDNA in elasmobranchs. PCR and sequence analysis unequivocally demonstrated for the first time that sharks can bear two size classes of 5S rDNA arrays carrying distinct NTS sequences and, surprisingly, highly variable 5S genes. Different units of 5S rDNA, especially those characterizing a double-class pattern, have also been previously documented in several teleost species (Pendás et al., Reference Pendás, Móran, Martínez and Garcia-Vásquez1995; Sadjak et al., Reference Sajdak, Reed and Phillips1998; Martins & Galetti, Reference Martins and Galetti2001a; Alves-Costa et al., Reference Alves-Costa, Wasko, Oliveira, Foresti and Martins2006) and some elasmobranchs of the family Rajidae (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006).

Sharks exhibited a higher similarity in their 5S rDNA units within a specific class between two species than between the two classes in the same species. This probably occurred because 5S rDNA sequence tends to homogenize the different copies that are arrayed in the same cluster, which can differ extensively from the copies of a second 5S cluster (Martins & Galetti, Reference Martins and Galetti2001a; Martins & Wasko, Reference Martins and Wasko2004). This observation was clearly documented in the organization of the 5S rDNA in the Nile tilapia Oreochromis niloticus (Martins et al., Reference Martins, Wasko, Oliveira, Porto-Foresti, Parise-Maltempi, Wright and Foresti2002) and in the South American species of the genus Leporinus (Martins & Galetti, Reference Martins and Galetti2001a). It has been suggested that multigene families evolve according to homogenization processes governed by molecular drive and concerted evolution, resulting in a sequence similarity of the repeat units that is greater within than between species (Dover, Reference Dover1986; Elder & Turner, Reference Elder and Turner1995). Molecular pathways that lead to such concerted evolution are unequal crossing-over, gene conversion and replication slippage, but the rate at which the variant repeat types become homogenized depends on other aspects such as the number of repeats in an array, the strength of natural selection and the effective population size (Schlötterer & Tautz, Reference Schlötterer and Tautz1994). Thus, the discrepancy observed between the 5S rDNA classes obtained (i.e. variable length and base composition) must be the product of normal levels of divergence between orthologous sequences in the Rhizoprionodon genome.

Although the existence of variant copies of 5S rDNA has been reported in several mammals and amphibians (Komiya et al., Reference Komiya, Hasegawa and Takemura1986; Little & Braaten, Reference Little and Braaten1989; Leah et al., Reference Leah, Frederiksen, Engberg and Sorensen1990; Frederiksen et al., Reference Frederiksen, Cao, Lomholt, Levan and Hallemberg1997), these copies did not form distinct repetitive classes in the species investigated. In contrast, the presence of two 5S rDNA classes seems to be a general trend for teleost fishes (Martins & Galetti, Reference Martins and Galetti2001b; Martins & Wasko, Reference Martins and Wasko2004). Moreover, it is likely that the two 5S rDNA classes correspond to two distinct subfamilies of 5S genes. These findings uphold the complex organization of this repeated element in fishes. Furthermore, the distinct 5S rDNA classes detected in teleost fishes were well characterized by highly variable classes of NTSs and a single, highly conserved class of 5S gene.

Conversely, the present data, added to those of Rajidae from Pasolini et al. (Reference Pasolini, Costagliola, Rocco and Tinti2006), not only account for a double-class 5S rDNA pattern, but also for a double form of the 5S gene shared by elasmobranchs. This assumption is based on the consistent mutations observed in sequences responsible for the clustering of 5S genes of class II from sharks and rays into an ‘elasmobranch clade’ on relationship trees (Fig. 4) and is supported by a large dataset that we generated for sharks (46 complete 5S gene sequences of two species) plus other data from Rajidae acquired in GenBank (30 complete 5S gene sequences of six species). The presence of control elements in the NTS-I and NTS-II, the T-rich tail in the 3′ end of the gene, and conserved ICRs in 5S genes of both 5S rDNA classes of Rhizoprionodon species suggest that both classes of 5S rDNA might be functional in sharks. The occurrence of 5S genes that are differently regulated in somatic and oocyte cells has been previously described for several vertebrates, including lampreys and teleost fishes (Komiya et al., Reference Komiya, Hasegawa and Takemura1986), and it could explain the dual 5S rDNA pattern that is widely reported among teleosts. However, both somatic and oocyte types are similar to the genes found in the 5S rDNA class I of sharks. Considering that a dual 5S rDNA system, similar to that detected among vertebrates, was also described for plants of distant taxa such as Liliopsida and Coniferopsida (Shibata & Hizume, Reference Shibata and Hizume2002; Besendorfer et al., Reference Besendorfer, Krajačić-Sokol, Jeleni, Puizina, Mlinarec, Sviben and Papes2005), it seems probable that the presence of a dual system of repeat units of 5S rDNA has a major biological role in the cells. Although the information on Rhizoprionodon needs to be corroborated by further studies on a larger number of Elasmobranchii, including many sharks and other freshwater and marine rays, the present evidence suggests a group-specific evolutionary history of 5S rDNA in these animals. Additional studies will be helpful to clarify their expression and functionality.

(ii) NTS features

Numerous studies have shown that the 5S rDNA is a valuable molecular marker for fish species (Pendás et al., Reference Pendás, Móran, Martínez and Garcia-Vásquez1995; Sadjak et al., Reference Sajdak, Reed and Phillips1998; Carrera et al., Reference Carrera, Garcia, Céspedes, González, Fernández, Asensio, Hernández and Martin2000), essentially because of the low intra-specific and high inter-specific variability of NTS domains. Among mammals, sequence differences in the ribosomal NTS were sufficient to distinguish the mouse subspecies Mus musculus domesticus from M. m. musculus (Suzuki et al., Reference Suzuki, Moriwaki and Sakurai1994), suggesting that this region evolved at an appropriate rate for differentiating closely related species. The substitution rate of 0·58%/Myr within the NTS of mice (Suzuki et al., Reference Suzuki, Moriwaki and Sakurai1994) is as rapid as that of mammalian pseudogenes and non-functional sequences (Li et al., Reference Li, Luo and Wu1985).

Among the two Rhizoprionodon sharks, a variable degree of similarity was found within each NTS class. Whereas the NTS II showed a highly conserved pattern, the NTS I was diverse and appeared to have evolved faster than the other class. However, NTS II also separated R. lalandii from R. porosus in phylogenetic trees, due to five consistent base substitutions, showing that the NTS II can be useful for phylogenetic inference in shark species. This result is in agreement with those obtained for the neotropical characiform Leporinus (Ferreira et al., Reference Ferreira, Oliveira, Venere, Galetti and Martins2006) and Brycon (Wasko et al., Reference Wasko, Martins, Wright and Galetti2001), in which the NTS effectively separated closely related species, but it contrasts with the results obtained by Pasolini et al. (Reference Pasolini, Costagliola, Rocco and Tinti2006) for the equivalent NTS II existing in rays, as well as those by Robles et al. (Reference Robles, de la Herran, Ludwig, Rejon, Rejon and Garrido-Ramos2005) for sturgeons and by Sajdak et al. (Reference Sajdak, Reed and Phillips1998) for Coregonus. In view of these contradictory results, NTS sequences must be used with caution in evolutionary surveys.

Apart from these sequence variations, it is interesting that the remaining 340 nucleotides of NTS II of R. lalandii and R. porosus matched almost perfectly. Given the high mutation rates of NTS observed in every fish genome so far studied (for review see Martins & Wasko, Reference Martins and Wasko2004), as well as in mammals (Suzuki et al., Reference Suzuki, Moriwaki and Sakurai1994), the large sequence conservation of NTS II in the Rhizoprionodon genome indicates that these non-coding domains are not as free from selective pressure as other non-coding sequences. In this context, it has been recently demonstrated that the presence of conserved elements located within this NTS also plays an important role in the regulation of the 5S gene expression in mammals (Nederby-Nielsen et al., Reference Nederby-Nielsen, Hallenberg, Frederiksen, Sorensen and Lomholt1993; Hallenberg & Frederiksen, Reference Hallenberg and Frederiksen2001). In fishes, a conserved TATA-like sequence has been observed upstream of the 5S gene in several teleosts, including Salmo salar (Pendás et al., Reference Pendás, Moran, Freije and Garcia-Vasquez1994), Carassius auratus (Murakami & Fugitani, Reference Murakami and Fujitani1998), Coregonus (Sajdak et al., Reference Sajdak, Reed and Phillips1998), Gasterosteus aculeatus (Rocco et al., Reference Rocco, Russo, Stingo, Aprea and Odierna1999), Acheilognathus tabira and Cyprinus carpio (Inafuku et al., Reference Inafuku, Nabeyama, Kikuma, Saitoh, Kubota and Kohono2000), O. niloticus (Martins et al., Reference Martins, Wasko, Oliveira and Wright2000), Brycon (Wasko et al., Reference Wasko, Martins, Wright and Galetti2001) and Leporinus (Martins & Galetti, Reference Martins and Galetti2001a), suggesting a possible influence on the transcription level of this gene. In elasmobranchs, such sequences were also detected in Rajidae (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006), but no data regarding sharks have been reported previously. In the present study, the analysis of upstream sequences of 5S genes from Rhizoprionodon identified an upstream TATA-like control element at positions −25 in NTS II and −30 in NTS I; interestingly, these are the same positions recognized in Rajidae (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006). Although the role of these sequences for 5S gene transcription is still unknown, the high degree of conservation of TATA-like sequence positions in the two elasmobranch species currently studied is intriguing. This subject suggests some structural pattern for these elements shared by elasmobranchs, in opposition to ray-finned fishes, where the TATA-like sequences had an irregular position in different genera.

Microsatellite elements in the NTSs have also been found in some fish species (Deiana et al., Reference Deiana, Cau, Salvadori, Coluccia, Cannas, Milia and Tagliavini2000; Gornung et al., Reference Gornung, De Innocentiis, Annesi and Sola2000; Alves-Costa et al., Reference Alves-Costa, Wasko, Oliveira, Foresti and Martins2006), including elasmobranchs (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006). These elements could be linked to regulatory functions of the 5S rRNA gene, which is not yet fully elucidated (Ota et al., Reference Ota, Tateno and Gojobori2003). A similar example was found in Hox cluster sequences of sharks and mammals, which share simple repetitive elements in the intervening sequence between coding genes, implying conservation of large-scale structural properties of non-coding regions (Kim et al., Reference Kim, Amemiya, Bailey, Kawasaki, Mezey, Miller, Minoshima, Shimizu, Wagner and Ruddle2000). In Rhizoprionodon sharks, the single typical microsatellite (TCCC)n recovered in NTS class II was highly conserved, especially in the R. porosus genome, where the same number of repeats was recovered in all the clones studied (Fig. 2). However, in general, the scarcity of microsatellites in NTS of Rhizoprionodon sharks compared with their abundance in rajid rays (Pasolini et al., Reference Pasolini, Costagliola, Rocco and Tinti2006) eliminates the possibility of some global biological role of the 5S rDNA in elasmobranchs, but it does not negate its usefulness to exemplify the dynamic evolution experienced by NTS domains, even in related genomes.

(iii) Comparative analysis of 5S genes among several vertebrate lineages

According to many studies, modern elasmobranchs represent the vertebrate basal condition because of their long (about 400 mya) evolutionary history, although their relationships with subsequent vertebrate groups are unclear. Cartilaginous fishes have been believed to correspond to the sister group of Actinopterygii and Sarcopterygii (including Tetrapoda), and to have diverged from these vertebrate lineages about 500 mya (Kumar & Hedges, Reference Kumar and Hedges1998; Nelson, Reference Nelson2006). The ancient condition of elasmobranchs, together with their remarkable endurance, explains the interest in this group as pivotal in understanding early vertebrate evolution (Grogan & Lund, Reference Grongan and Lund2004).

Molecular studies have included sharks or rays when describing the structure and expression of particular genes and proteins in model organisms (e.g. human, mouse and chicken) to gain perspective on conservation of function (Martin, Reference Martin2001). A more complete characterization of the cartilaginous fish genome is, therefore, essential for gaining insights into the ancestral condition of other vertebrate genomes. Towards such a characterization, our analyses on sharks showed that the clusters bearing the 5S rDNA classes I and II of Rhizoprionodon seem to evolve rapidly and independently. This can be confirmed by examination of the variability found in NTS and, surprisingly, in 5S gene sequences. Both the 5S rDNA classes of teleosts as well as the different copies of 5S rDNA arrays, which comprise a unique class in the amphibians, reptiles, birds and mammals thus far studied, seem to have arisen from the same ancestral 5S rDNA class I found in Rhizoprionodon sharks, Rajidae rays and lampreys. Conversely, the 5S rDNA class II probably embodies derived sequences of Elasmobranchii, as confirmed by the Bayesian analysis (Fig. 4). Interestingly, these 5S rDNA class II sequences from sharks have no recognizable homologues in lampreys, which also bear the ancestral 5S rDNA class I. It remains unclear if lampreys may also have two classes of 5S rDNA, as detected in Elasmobranchii and Teleostei, but it is probable that they contain variant copies corresponding to a single class, as is well documented in tetrapod genomes (Rosenthal & Doering, Reference Rosenthal and Doering1983; Suzuki et al., Reference Suzuki, Moriwaki and Sakurai1994; Frederiksen et al., Reference Frederiksen, Cao, Lomholt, Levan and Hallemberg1997; Jensen & Frederiksen, Reference Jensen and Frederiksen2000). Regarding the number of 5S rDNA classes in the genome of sharks, agarose electrophoresis profiles of 5S rDNA-PCR products from Sphyrna lewini, G. cuvier, Carcharhinus obscurus, Carcharhinus leucas, Carcharhinus limbatus, Carcharhinus acronotus, A. superciliosus and Isurus oxyrinchus have indicated the extensive variability present in 5S rDNA unit size (Pinhal et al., Reference Pinhal, Gadig, Wasko, Oliveira, Foresti and Martins2008). Each species showed an exclusive pattern with one to four bands of different molecular weights, possibly referring to distinct paralogous copies of 5S rDNA arrays. We believe that these copies correspond to the two classes of 5S rDNA, whereas extra copies might represent variants and pseudogenes, similar to those previously detected in rat, mouse and human cells (Emerson & Roeder, Reference Emerson and Roeder1984; Little & Braaten, Reference Little and Braaten1989; Leah et al., Reference Leah, Frederiksen, Engberg and Sorensen1990; Nederby-Nielsen et al., Reference Nederby-Nielsen, Hallenberg, Frederiksen, Sorensen and Lomholt1993) and in fish genomes (Murakami & Fujitami, Reference Murakami and Fujitani1998; Roest-Crollius, Reference Roest-Crollius, Jaillon, Bernot, Dasilva, Bouneau, Fischer, Fizames, Wincker, Brottier, Quetier, Saurin and Weissenbach2000; Martins et al., Reference Martins, Wasko, Oliveira, Porto-Foresti, Parise-Maltempi, Wright and Foresti2002; Robles et al., Reference Robles, de la Herran, Ludwig, Rejon, Rejon and Garrido-Ramos2005). Species of sharks that showed merely one band on gels (e.g. G. cuvier and A. superciliosus) could either carry two classes of a similar mass or one single class, although more than PCR data are necessary to test these hypotheses. The loss of a cluster might have occurred secondarily in fish taxa that bear only one class of 5S rDNA, as reported for Synbranchus marmoratus (Messias et al., Reference Messias, Ferreira, Wasko, Oliveira, Foresti and Martins2003), Acipenser sturio (Tagliavini et al., Reference Tagliavini, Williot, Congiu, Chicca, Lanfredi, Rossi and Fontana1999) and G. aculeatus (Rocco et al., Reference Rocco, Russo, Stingo, Aprea and Odierna1999), perhaps caused by a strong reduction of the NTS (less than 60–80 bp length) below the minimum required to sustain the array and the dynamics of the 5S rDNA in the genome (Martins & Galetti, Reference Martins and Galetti2001b). To date, the short-class array length of 5S rDNA in the R. lalandii and R. porosus genomes has always been found to be greater than the apparently critical 60 bp length (Table 1), the same length as found for G. cuvier and A. superciliosus (data not shown but available upon request).

Considering the three round hypothesis (named ‘3R’) of genome duplications in the evolutionary history of chordates (McLysaght et al., Reference McLysaght, Hokamp and Wolfe2002; Donoghue & Purnell, Reference Donoghue and Purnell2005), it is likely that such events created new paralogous copies in 5S rDNA families. Indeed, recent molecular studies have suggested that the genome duplication event took place in the teleostean stem lineages, and many genes constitute more than two families in teleostean genomes (Hoegg et al., Reference Hoegg, Brinkmann, Taylor and Meyer2004; Crow et al., Reference Crow, Stadler, Lynch, Amemiya and Wagner2005), making it feasible for 5S rDNA arrays. Nevertheless, almost all studies on this subject have claimed that 5S rDNA arrays in the teleostean genome evolved into a double-class (or double-family) system (see Martins & Wasko, Reference Martins and Wasko2004). Therefore, we can conclude that the extra copies found dispersed throughout their genomes probably represent variant or even inactive forms and do not constitute a third class of 5S rDNA. The shark genome, in contrast, has been suggested to have undergone only two rounds of duplications (Froschauer et al., Reference Froschauer, Braasch and Volff2006), and duplicated segments encompassing genes and non-codifying sequences could evolve independently in individual species, being lost or maintained in subsequent generations, as described in other fish taxa (Crollius & Weissenbach, Reference Crollius and Weissenbach2005). These conclusions provide support for the different varieties of 5S rDNA classes currently detected in these animals.

Although extensive taxa sampling is needed to deduce the origin of 5S rDNA classes of sharks in vertebrate evolution, we used corresponding sequences available from model organisms such as human, Xenopus, chicken and fugu, added to those of basal taxa such as lampreys, to generate a consistent phylogenetic report. In such an evolutionary scenario, the 5S rDNA class I in the shark and lamprey genome seems to represent the widespread ancestral condition, retained in the course of evolution and still extant in tetrapod genomes (Fig. 4), whereas the 5S rDNA class II apparently evolved in parallel in Elasmobranchii, explaining the existence of the paralogous copies of 5S genes in these animals.

Additional studies including representatives from the family Potamotrygonidae of freshwater stingrays, additional lineages of Vertebrata, and shark species such as those of the orders Lamniformes and Squatiniformes, and including sequences of different genes will be very valuable for testing the relationship pattern found in the current study and for improving our knowledge of the evolution of 5S rDNA.

The authors are grateful to FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) and CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) for financial support.

References

Aljanabi, S. M. & Martinez, I. (1997). Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Research 25, 46924693.CrossRefGoogle ScholarPubMed
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). Basic local alignment search tool. Journal of Molecular Biology 215, 403410.Google Scholar
Alves-Costa, F. A., Wasko, A. P., Oliveira, C., Foresti, F. & Martins, C. (2006). Genomic organization and evolution of the 5S ribosomal DNA in Tilapiini fishes. Genetica 127, 243252.CrossRefGoogle ScholarPubMed
Barciszewska, M. Z., Szymanski, M., Erdmann, V. A. & Barciszewski, J. (2001). Structure and functions of 5S rRNA. Acta Biochimica Polonica 48, 191198.CrossRefGoogle Scholar
Besendorfer, V., Krajačić-Sokol, I., Jeleni, S., Puizina, J., Mlinarec, J., Sviben, T. & Papes, D. (2005). Two classes of 5S rDNA unit arrays of the silver fir, Abies alba Mill.: structure, localization and evolution. Theoretical and Applied Genetics 110, 730741.CrossRefGoogle ScholarPubMed
Carrera, E., Garcia, T., Céspedes, A., González, I., Fernández, A., Asensio, L. M., Hernández, P. E. & Martin, R. (2000). Differentiation of smoked Salmo solar, Onconhynchus mykiss and Brama raii using the nuclear marker 5S rDNA. International Journal of Food Science and Technology 35, 401406.CrossRefGoogle Scholar
Céspedes, A., Garcia, T., Carrera, E., Gonzalez, I., Fernandez, A., Hernandez, P. E. & Martin, R. (1999). Identification of sole (Solea solea) and Greenland halibut (Reinhardtius hippoglossoides) by PCR amplification of the 5S rDNA gene. Journal of Agricultural and Food Chemistry 47, 10461050.CrossRefGoogle ScholarPubMed
Compagno, L. J. V. (1984). FAO Species Catalogue, Volume 4: Sharks of the World. Part 2: Carcharhiniformes. FAO Fisheries. Synopsis 125, 251655.Google Scholar
Crollius, H. R. & Weissenbach, J. (2005). Fish genomics and biology. Genome Research 15, 16751682.CrossRefGoogle Scholar
Crow, K. D., Stadler, P. F., Lynch, V. J., Amemiya, C. & Wagner, G. P. (2005). The ‘fish specific’ Hox cluster duplication is coincident with the origin of teleosts. Molecular Biology and Evolution 23, 121136.Google Scholar
Deiana, A. M., Cau, A., Salvadori, S., Coluccia, E., Cannas, R., Milia, A. & Tagliavini, J. (2000). Major and 5S ribosomal sequences of the largemouth bass Micropterus salmoides (Perciformes, Centrarchidae) are localized in GC-rich regions of the genome. Chromosome Research 8, 213218.Google Scholar
Donoghue, P. C. J. & Purnell, M. A. (2005). Gene duplication, extinction, and vertebrate evolution. Trends in Ecology and Evolution 20, 312319.CrossRefGoogle ScholarPubMed
Dover, G. A. (1986). Linkage disequilibrium and molecular drive in the rDNA gene family. Genetics 122, 249252.CrossRefGoogle Scholar
Elder, J. F. Jr, & Turner, B. J. (1995). Concerted evolution of repetitive DNA sequences in eukaryotes. The Quarterly Review of Biology 70, 277320.CrossRefGoogle ScholarPubMed
Emerson, B. M. & Roeder, R. G. (1984). Isolation and genomic arrangement of active and inactive forms of mammalian 5S RNA genes. Journal of Biological Chemistry 259, 79167925.Google Scholar
Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783791.CrossRefGoogle ScholarPubMed
Ferreira, I. A., Oliveira, C., Venere, P. C., Galetti, P. M. Jr & Martins, C. (2006). 5S rDNA variation and its phylogenetic inference in the genus Leporinus (Characiformes: Anostomidae). Genetica 129, 253257.CrossRefGoogle ScholarPubMed
Frederiksen, S., Cao, H., Lomholt, B., Levan, G. & Hallemberg, C. (1997). The rat 5S rRNA bona fide gene repeat maps to chromosome 19q12→qter and the pseudogene repeat maps to 12q12. Cytogenetics and Cell Genetics 76, 101106.Google Scholar
Froschauer, F., Braasch, I. & Volff, J. N. (2006). Fish genomes, comparative genomics and vertebrate evolution. Current Genomics 7, 4357.CrossRefGoogle Scholar
Gadig, O. B. F. (2001). Tubarões da costa brasileira. PhD Thesis, São Paulo State University, 343 pp.Google Scholar
Gornung, E., De Innocentiis, S., Annesi, F. & Sola, L. (2000). Zebrafish 5S rRNA genes map to the long arms of chromosome 3. Chromosome Research 8, 362.CrossRefGoogle Scholar
Grongan, E. D. & Lund, R. (2004). The origin and relationships of early Chondrichthyes. In Biology of Sharks and Rays, pp. 331. New York: CRC Press.Google Scholar
Guindon, S. & Gascuel, O. (2003). A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology 52, 696704.CrossRefGoogle ScholarPubMed
Guindon, S., Lethiec, F., Duroux, P. & Gascuel, O. (2005). PHYML Online – a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Research 33, W557W559.CrossRefGoogle ScholarPubMed
Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 41, 9598.Google Scholar
Hallenberg, C., Nederby-Nielsen, J. & Frederiksen, S. (1994). Characterization of 5S rRNA genes from mouse. Gene 142, 291295.CrossRefGoogle ScholarPubMed
Hallenberg, C. & Frederiksen, S. (2001). Effect of mutations in the upstream promoter on the transcription of human 5S rRNA genes. Biochimica et Biophysica Acta 1520, 169173.CrossRefGoogle ScholarPubMed
Hoegg, S., Brinkmann, H., Taylor, J. S. & Meyer, A. (2004). Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish. Journal of Molecular Evolution 59, 190203.Google Scholar
Huelsenbeck, J. P., Ronquist, F., Nielsen, R. & Bollback, J. P. (2001). Bayesian inference of phylogeny and its impact of evolutionary biology. Science 294, 23102314.CrossRefGoogle Scholar
Inafuku, J., Nabeyama, M., Kikuma, Y., Saitoh, J., Kubota, S. & Kohono, S. (2000). Chromosomal location and nucleotide sequences of 5S ribosomal DNA of two cyprinid species (Osteichthyes, Pisces). Chromosome Research 8, 193199.Google Scholar
Jensen, L. R. & Frederiksen, S. (2000). The 5S rRNA genes in Macaca fascicularis are organized in two large tandem repeats. Biochimica et Biophysica Acta 1492, 537542.CrossRefGoogle ScholarPubMed
Kim, C. B., Amemiya, C., Bailey, W., Kawasaki, K., Mezey, J., Miller, W., Minoshima, S., Shimizu, N., Wagner, G. & Ruddle, F. (2000). Hox cluster genomics in the horn shark, Heterodontus francisci. Proceedings of the National Academy of Sciences of the USA 97, 16551660.Google Scholar
Komiya, H. & Takemura, S. (1979). Nucleotide sequence of 5S ribosomal RNA from rainbow trout (Salmo gairdnerii) liver. Journal of Biochemistry 86, 10671080.CrossRefGoogle ScholarPubMed
Komiya, H., Hasegawa, M. & Takemura, S. (1986). Differentiation of oocyte- and somatic-type 5S rRNAs in animals. Journal of Biochemistry 100, 369374.Google Scholar
Korn, L. J. & Brown, D. D. (1978). Nucleotide sequence of Xenopus borealis oocyte 5S DNA: comparison of sequences that flank several eukaryotic genes. Cell 15, 11451156.CrossRefGoogle Scholar
Kumar, S. & Hedges, S. B. (1998). A molecular timescale for vertebrate evolution. Nature 392, 917920.Google Scholar
Leah, R., Frederiksen, S., Engberg, J. & Sorensen, P. D. (1990). Nucleotide sequence of a mouse 5S rRNA variant gene. Nucleic Acids Research 18, 7441.Google Scholar
Li, W.-H., Luo, C.-C. & Wu, C.-I. (1985). A new method for estimating the synonymous and nonsynonymous rates of nucleotides substitution considering the relative likelihood of nucleotide and codon changes. Molecular Biology and Evolution 2, 150174.Google ScholarPubMed
Little, R. & Braaten, D. (1989). Genomic organization of human 5S rDNA and sequence of one tandem repeat. Genomics 4, 376383.Google Scholar
Long, E. O. & David, I. B. (1980). Repeated genes in Eukaryotes. Annual Reviews in Biochemistry 49, 727764.Google Scholar
Martin, A. (2001). The phylogenetic placement of Chondrichthyes: inferences from analysis of multiple genes and implications for comparative studies. Genetica 111, 349357.CrossRefGoogle ScholarPubMed
Martins, C., Wasko, A. P., Oliveira, C. & Wright, J. M. (2000). Nucleotide sequence of 5S rDNA and localization of the ribosomal RNA genes to metaphase chromosomes of the tilapiine cichlid fish, Oreochromis niloticus. Hereditas 133, 3946.CrossRefGoogle ScholarPubMed
Martins, C. & Galetti, P. M. Jr (2001 a). Organization of 5S rDNA in species of the fish Leporinus: two different genomic locations are characterized by distinct non-transcribed spacers. Genome 44, 903910.Google Scholar
Martins, C. & Galetti, P. M. Jr (2001 b). Two 5S rDNA arrays in neotropical fish species: is it a general rule for fishes? Genetica 111, 439446.CrossRefGoogle ScholarPubMed
Martins, C., Wasko, A. P., Oliveira, C., Porto-Foresti, F., Parise-Maltempi, P. P., Wright, J. M. & Foresti, F. (2002). Dynamics of 5S rDNA in the tilapia (Oreochromis niloticus) genome: repeat units, inverted sequences, pseudogenes and chromosome loci. Cytogenetics and Genome Research 98, 7885.Google Scholar
Martins, C. & Wasko, A. P. (2004). Organization and evolution of 5S ribosomal DNA in the fish genome. In Focus on Genome Research, pp. 335363. New York: Nova Science Publishers.Google Scholar
McLysaght, A., Hokamp, K. & Wolfe, K. H. (2002). Extensive genomic duplication during early chordate evolution. Nature 31, 200202.Google ScholarPubMed
Messias, L. H., Ferreira, D. C., Wasko, A. P., Oliveira, C., Foresti, F. & Martins, C. (2003). 5S rDNA organization in the fish Synbranchus marmoratus (Synbranchidae, Synbranchiformes). Hereditas 139, 228231.Google Scholar
Móran, P., Martínez, J. L., Garcia-Vásquez, E. & Pendás, A. M. (1996). Sex linkage of 5S rDNA in rainbow trout (Oncorhynchus mykiss). Cytogenetics and Cell Genetics 75, 145150.CrossRefGoogle ScholarPubMed
Motta, F. S., Gadig, O. B. F., Namora, R. C. & Braga, F. M. S. (2005). Size and sex compositions, length–weight relationship, and occurrence of the Brazilian sharpnose shark, Rhizoprionodon lalandii, caught by artisanal fishery from southeastern Brazil. Fisheries Research 74, 116126.Google Scholar
Murakami, M. & Fujitani, H. (1998). Characterization of repetitive DNA sequences carrying 5S rDNA of the triploid ginbuna (Japanese silver crucian carp, Carassius auratus langsdorfi). Genes and Genetic Systems 73, 920.Google Scholar
Nederby-Nielsen, J., Hallenberg, C., Frederiksen, S., Sorensen, P. D. & Lomholt, B. (1993). Transcription of human 5S rRNA genes is influenced by an upstream DNA sequence. Nucleic Acids Research 26, 36313636.CrossRefGoogle Scholar
Nelson, J. S. (2006). Fishes of the World. New York: John Wiley and Sons.Google Scholar
Ota, K., Tateno, Y. & Gojobori, T. (2003). Highly differentiated and conserved sex chromosome in fish species (Aulopus japonicus: Teleostei, Aulopidae). Gene 317, 187193.CrossRefGoogle ScholarPubMed
Pasolini, P., Costagliola, D., Rocco, L. & Tinti, F. (2006). Molecular organization of 5S rDNA in Rajidae (Chondrichthyes): structural features and evolution of piscine 5S rRNA genes and nontranscribed intergenic spacers. Journal of Molecular Evolution 62, 564574.Google Scholar
Pendás, A. M., Moran, P., Freije, J. P. & Garcia-Vasquez, E. (1994). Chromosomal mapping and nucleotide sequence of two tandem repeats of Atlantic salmon 5S rDNA. Cytogenetics and Cell Genetics 67, 3136.CrossRefGoogle Scholar
Pendás, A. M., Móran, P., Martínez, J. L. & Garcia-Vásquez, E. (1995). Applications of 5S rDNA in Atlantic salmon, brown trout, and in Atlantic salmon×brown trout hybrid identification. Molecular Ecology 4, 275276.CrossRefGoogle Scholar
Pinhal, D., Gadig, O. B. F., Wasko, A. P., Oliveira, C., Foresti, F. & Martins, C. (2008). Discrimination of shark species by simple PCR of 5S rDNA repeats. Genetics and Molecular Biology 31, 361365.CrossRefGoogle Scholar
Posada, D. & Crandall, K. A. (1998). Modeltest: testing the model of DNA substitution. Bioinformatics 14, 817818.Google Scholar
Robles, F., de la Herran, R., Ludwig, A., Rejon, C. R., Rejon, M. R. & Garrido-Ramos, M. A. (2005). Genomic organization and evolution of the 5S ribosomal DNA in the ancient fish sturgeon. Genome 48, 1828.CrossRefGoogle ScholarPubMed
Rocco, L., Russo, C., Stingo, V., Aprea, G. & Odierna, G. (1999). Characterization of 5S rDNA in Gasterosteus aculeatus (Teleostei, Gasterosteidae). Italian Journal of Zoology 66, 285289.CrossRefGoogle Scholar
Rocco, L., Costagliola, D., Fiorillo, M., Tinti, F. & Stingo, V. (2005). Molecular and chromosomal analysis of ribosomal cistrons in two cartilaginous fish, Taeniura lymma and Raja montagui (Chondrichthyes, Batoidea). Genetica 123, 245253.CrossRefGoogle ScholarPubMed
Roest-Crollius, H., Jaillon, O., Bernot, A., Dasilva, C., Bouneau, L., Fischer, C., Fizames, C., Wincker, P., Brottier, P., Quetier, F., Saurin, W. & Weissenbach, J. (2000). Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence. Nature Genetics 25, 235238.CrossRefGoogle ScholarPubMed
Ronquist, F. & Huelsenbeck, J. P. (2003). MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 15721574.Google Scholar
Rosa, R. S., Gadig, O. B. F., Santos Motta, F. & Namora, R. C. (2004). Rhizoprionodon lalandii. In: IUCN 2008. 2008 IUCN Red List of Threatened Species. Available at http://www.iucnredlist.org (downloaded on 15 November 2008).Google Scholar
Rosenthal, D. S. & Doering, J. L. (1983). The genomic organization of dispersed tRNA and 5 S RNA genes in Xenopus laevis. Journal of Biological Chemistry 258, 74027410.Google Scholar
Sajdak, S. L., Reed, K. M. & Phillips, R. B. (1998). Intraindividual and interspecies variation in the 5S rDNA of coregonid fish. Journal of Molecular Evolution 46, 680688.CrossRefGoogle ScholarPubMed
Sambrook, J. & Russel, D. W. (2001). Molecular Cloning: A Laboratory Manual. New York: Cold Spring Harbor Laboratory Press.Google Scholar
Schlötterer, C. & Tautz, D. (1994). Chromosomal homogeneity of Drosophila ribosomal DNA arrays suggests intrachromosomal exchanges drive concerted evolution. Current Biology 4, 777783.CrossRefGoogle ScholarPubMed
Shibata, F. & Hizume, M. (2002). Evolution of 5S rDNA units and their chromosomal localization in Allium cepa and Allium schoenoprasum revealed by microdissection and FISH. Theoretical and Applied Genetics 105, 167172.CrossRefGoogle ScholarPubMed
Suzuki, H., Moriwaki, K. & Sakurai, S. (1994). Sequences and evolutionary analysis of mouse 5S rDNAs. Molecular Biology and Evolution 11, 704710.Google Scholar
Swofford, D. L., Olsen, G. J., Waddell, P. J. & Hillis, D. M. (1996). Phylogenetic inference. In Molecular Systematics, pp. 407514. Sunderland, MA: Sinauer Associates.Google Scholar
Tagliavini, J., Williot, P., Congiu, L., Chicca, M., Lanfredi, M., Rossi, R. & Fontana, F. (1999). Molecular cytogenetic analysis of the karyotype of the European Atlantic sturgeon, Acipenser sturio. Heredity 83, 520525.Google Scholar
Tamura, K. & Nei, M. (1993). Estimation of the number of nucleotide substitution in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution 10, 512526.Google Scholar
Tamura, K., Dudley, J., Nei, M. & Kumar, S. (2007). MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24, 15961599.Google Scholar
Tigano, C., Rocco, L., Ferrito, V., Costagliola, D., Pappalardo, A. M. & Stingo, V. (2004). Chromosomal mapping and molecular characterization of ribosomal RNA genes in Lebias fasciata (Teleostei, Cyprinodontidae). Genetica 121, 95100.Google Scholar
Wasko, A. P., Martins, C., Wright, J. M. & Galetti, P. M. Jr (2001). Molecular organization of 5S rDNA in fishes of the genus Brycon. Genome 44, 893902.CrossRefGoogle ScholarPubMed
Wegnez, M., Denis, H., Mazabraud, A. & Clerot, J. C. (1978). RNA accumulation during oogenesis of the dogfish Scyliorhinus caniculus. Biochemical research on oogenesis. Developmental Biology 62, 99111.Google Scholar
Xia, X. & Xie, Z. (2001). DAMBE: data analysis in molecular biology and evolution. Journal of Heredity 92, 371373.Google Scholar
Figure 0

Fig. 1. Agarose gel showing the two 5S rDNA types of R. lalandii (1) and R. porosus (2) amplified by PCR using the primers 5SA and 5SB. M refers to the 1 kb plus molecular weight marker in base pairs (bp).

Figure 1

Fig. 2. Representative sequences of 5S rDNA class I (5S I) and class II (5S II) from R. lalandii (RL) and R. porosus (RP). The start and end points of transcription are indicated by +1 and +120, respectively; the nucleotide substitutions of specific classes are indicated in grey shading, and the ICRs (A box, IE and C box) are indicated in black shading. TATA-like sequences are underlined. SSRs are in boldface. Hyphens represent gaps, and dots indicate identical nucleotides.

Figure 2

Table 1. Genetic distance of the 5S genes and NTS classes among Rhizoprionodon sharks. NC, number of clones; L, length in base pairs; GD, genetic distance; SE, standard error

Figure 3

Table 2. Diagnostic nucleotide sites present in 5S gene class I and class II from R. lalandii and R. porosus

Figure 4

Fig. 3. Recovered ML trees for (a) the entire 5S rDNA sequences, (b) 5S gene only, (c) all NTS sequences and (d) NTS II only. Note the existence of two main branches in the trees discriminating the two 5S rDNA classes. In these trees, sequence names correspond to the first two letters of their species names (RL, R. lalandii and RP, R. porosus) and the number of the clones analysed. The scale bar indicates genetic distance based on the Tamura–Nei model analysis. Bootstrap values under 70 were omitted.

Figure 5

Table 3. Diagnostic nucleotide sites present in NTS class I sequences of Rhizoprionodon sharks

Figure 6

Fig. 4. A phylogenetic tree of 5S rRNA genes from Elasmobranchii, Teleostei and several other vertebrates using the Bayesian-likelihood method. Bootstrap values under 70 were omitted.