- Research article
- Open Access
Molecular organization and chromosomal localization of 5S rDNA in Amazonian Engystomops(Anura, Leiuperidae)
BMC Genetics volume 13, Article number: 17 (2012)
For anurans, knowledge of 5S rDNA is scarce. For Engystomops species, chromosomal homeologies are difficult to recognize due to the high level of inter- and intraspecific cytogenetic variation. In an attempt to better compare the karyotypes of the Amazonian species Engystomops freibergi and Engystomops petersi, and to extend the knowledge of 5S rDNA organization in anurans, the 5S rDNA sequences of Amazonian Engystomops species were isolated, characterized, and mapped.
Two types of 5S rDNA, which were readily differentiated by their NTS (non-transcribed spacer) sizes and compositions, were isolated from specimens of E. freibergi from Brazil and E. petersi from two Ecuadorian localities (Puyo and Yasuní). In the E. freibergi karyotypes, the entire type I 5S rDNA repeating unit hybridized to the pericentromeric region of 3p, whereas the entire type II 5S rDNA repeating unit mapped to the distal region of 6q, suggesting a differential localization of these sequences. The type I NTS probe clearly detected the 3p pericentromeric region in the karyotypes of E. freibergi and E. petersi from Puyo and the 5p pericentromeric region in the karyotype of E. petersi from Yasuní, but no distal or interstitial signals were observed. Interestingly, this probe also detected many centromeric regions in the three karyotypes, suggesting the presence of a satellite DNA family derived from 5S rDNA. The type II NTS probe detected only distal 6q regions in the three karyotypes, corroborating the differential distribution of the two types of 5S rDNA.
Because the 5S rDNA types found in Engystomops are related to those of Physalaemus with respect to their nucleotide sequences and chromosomal locations, their origin likely preceded the evolutionary divergence of these genera. In addition, our data indicated homeology between Chromosome 5 in E. petersi from Yasuní and Chromosomes 3 in E. freibergi and E. petersi from Puyo. In addition, the chromosomal location of the type II 5S rDNA corroborates the hypothesis that the Chromosomes 6 of E. petersi and E. freibergi are homeologous despite the great differences observed between the karyotypes of the Yasuní specimens and the others.
Ribosomal RNAs (rRNA) are molecules that combine to form the basic structures of the small and large ribosomal subunits. rRNAs are transcribed by two distinct multigene families. The 45S rDNA family synthesizes the 18S, 5.8S and 28S rRNAs, and the 5S rDNA family transcribes the 5S rRNA (as reviewed in references [1–3]). The 5S ribosomal genes are found as conserved copies of 120-bp sequences arranged in tandem and interspersed with variable non-transcribed spacers (NTSs) that differ in length and nucleotide composition (see reviews in  and ). Within the conserved 120-bp sequence lies an internal control region (ICR) consisting of 3 characteristic regions: the A box, an intermediate element (IE), and the C box . These regions act as promoters for transcription; the A box is a general binding sequence for RNA polymerase III, and the intermediate element and the C box are interaction sites for the transcription factor TFIIIA . Another characteristic typically found in the presumed-functional 5S genes is a poly-T terminator region, as initially reported by Korn and Brown .
Despite the large variations observed for the NTSs found in 5S rDNA, several studies have reported that conserved elements are present in these regions and that they play an important role in the regulation of 5S rRNA gene expression, similar to TATA-like sequences , oligonucleotide sequences , and the D box found in mammals [7, 8]. It has also been reported that a C localized at position -1 of many previously described 5S rRNA gene sequences guarantees a correct and efficient transcription start .
The sequences of 5S rDNA have been used as genetic and cytogenetic markers in evolutionary studies and in the identification and comparison of species, hybrids and strains [9–13]. The 120-bp region is highly conserved even among unrelated species, making it possible to isolate the 5S rDNA repeats of several species based on the available sequence from another species that is not necessarily closely related. An interesting feature of the 5S rDNA NTS that is most likely widespread throughout all vertebrate groups is the presence of two differently sized sequences, a subject that has been extensively studied in fish (see reviews in [3, 14, 15]).
For anurans, knowledge combining the molecular organization of the 5S rDNA sequences and their chromosomal locations are known for Xenopus laevis, X. borealis, X. muelleri [16, 17], n families Ascaphidae, Discoglossidae e Ranidae [13, 17], Physalaemus ephippifer  and Physalaemus cuvieri . In Xenopus species  and P. cuvieri , two types of 5S rDNA sequences were found and mapped to distinct chromosomal sites. In Xenopus species, the distinct types of 5S rDNA, named oocyte- and somatic types, show a differential pattern of gene expression [6, 16, 17, 21–26]. In addition to these cases with combined molecular and chromosomal analyses, studies restricted to the nucleotide sequence of the transcribed region of the 5S rDNA have been reported for the anurans Gastrotheca riobambae , Bufo americanus, Rana pipiens and Rana catesbeiana .
Based on a genetic analysis of morphological data, the Engystomops genus was revalidated by Nascimento et al.  to allocate species previously grouped in the Physalaemus pustulosus group. Currently, nine species comprise the Engystomops genus (the source studies are cited in reference ), and some reproductive and genetic evidences suggest the occurrence of undescribed cryptic species among the Amazonian Engystomops [31–33]. The cytogenetic analyses of different populations of the Amazonian species Engystomops petersi and Engystomops freibergi revealed interesting inter- and intraspecific divergences, and the possible involvement of these cytogenetic variations in incipient speciation has been suggested [33, 34]. Nevertheless, the lack of available cytogenetic markers prevents the proper identification of chromosomal homeologies in these populations, and consequently, the hypothesized rearrangements responsible for the karyotypic divergences remain unknown. In the present study, we isolated and characterized the 5S rDNA sequences of E. freibergi and two populations of E. petersi and mapped their chromosomal locations, extending the knowledge of 5S rDNA organization at the genomic and chromosomal levels in these amphibians. These data facilitate the identification of possible homeologous chromosomes among the Engystomops, serving as an important contribution for further studies on chromosomal divergence in anurans.
5S ribosomal gene characterization and molecular analysis
PCR amplification of the segments containing the 5S rDNA from E. freibergi from Acre (ZUEC 9647), E. petersi from Yasuní (QCAZ 34948), and E. petersi from Puyo (QCAZ 34937) using the primers 5S-A and 5S-B (Figure 1) generated bands with fragments of approximately 750 and 200 bp in all experiments. The cloning of the E. freibergi sequences resulted in six recombinant colonies, three of which carried an insert of 201 bp, one with an insert of 764 bp, and two with inserts of 765 bp. All the cloned sequences contained a 118-bp region that corresponded to the coding sequence of the 5S ribosomal gene (Figure 2a and 2b) and NTS regions of either 84 bp (5SACR 201-1 to 3) or approximately 647 bp (5SACR 764 and 5SACR765-1 to 2). A comparison of the 84-bp fragments revealed four base substitutions among the sequences. A comparison of the 647-bp fragments showed 3 base substitutions and a deletion-insertion at position 292 (Table 1; Figure 2a and 2b).
Three recombinant clones containing the 5S rDNA of E. petersi from Yasuní were recovered. The insert 5SYAS201 was shorter than the others (5SYAS766 and 5SYAS774) due to variations in their NTSs, which were 84 bp long in the insert 5SYAS201 and approximately 650 bp long in the others (Table 1). When compared with the clone 5SYAS766, clone 5SYAS774 showed nine additional nucleotides in a microsatellite DNA region (indicated at positions 546 to 554 in Figure 2b).
Nineteen recombinant clones containing the 5S rDNA of E. petersi from Puyo were obtained. Eleven had larger NTSs (649 bp in clone 5SPUY766 and 652 bp in clones 5SPUY769-1 to 10), whereas the remaining inserts had small NTSs of 84 bp (Figure 2). The differences in the lengths of the larger NTS segments resulted from the variation of the number of repetitions of a microsatellite DNA sequence present in these regions (positions 552 to 554 in Figure 2b). Thus, based on the NTSs, the sequences obtained can be classified into two types: type I 5S rDNA (with the small, 84-bp NTS) and type II 5S rDNA (with the large, ~650-bp NTS) (Table 1). Interestingly, the presumed 5S rDNA coding region was also characteristic for each of these types (see Table 1), as could be inferred from the maximum likelihood analysis of this specific region, in which all the type I sequences were clustered together and apart from the type II sequences group (Figure 3).
As expected, comparison of the 120 bp of the coding regions of the type I and type II 5S rDNA of Engystomops with those of the other anuran species and selected fish species available in GenBank revealed great similarity (Figure 4). When the type I and type II sequences of P. cuvieri were excluded from the analysis, there was a higher similarity between the Engystomops type I 5S RNA gene and the remaining 5S sequences (88-90%) than between the latter and the Engystomops type II gene (82%). When the Engystomops and P. cuvieri 5S genes were compared, a higher similarity was observed among the sequences of the same type (Figure 4).
A higher similarity was also observed between the Engystomops and the P. cuvieri type I NTS (97% -100%) than between the type II NTS of the Engystomops species (~650 bp) and that of P. cuvieri (~580 bp) (90% - 92%). Because the type I NTS sequences were the same in all the Engystomops analyzed, and the type II NTS sequences were quite similar (average pairwise similarity: 98%; Figure 2), to better illustrate the comparison between the Engystomops and P. cuvieri NTSs, only the pairwise alignments between the sequences of single specimens of E. petersi from Puyo and of P. cuvieri are shown in Figure 5. The three elements (A box, intermediate element and C box) of the 5S gene internal control region (ICR) were identified in the presumed coding regions of all the sequences isolated from the three species of Engystomops (Figure 4). Interestingly, the control elements of the Engystomops type II sequences differed more from those available in GenBank than did the elements of the Engystomops type I sequences.
All sequences of the type I 5S rDNA contained the control element TATA-box in their NTS regions, which was located 25 bp upstream from the coding region. In addition, the type I sequences also had a TATA-like element located 13 nucleotides upstream from the TATA-box (Figure 2a). In the NTS of the type II sequences, an element similar to the TATA-box motif was also detected, but it was very distant from the +1 position of the presumed coding region of the 5S gene, occupying positions -345 to -350 (Figure 2b). Additionally, in the type II 5S rDNA sequences, were found two hexanucleotides (at positions -53 to -48 and -21 to -16; Figure 2b) resembling the regulatory hexanucleotides previously described for Xenopus .
Both the type I and II Engystomops 5S sequences showed a T-rich sequence containing 4 to 5T-residues starting at position -119, which may correspond to the poly-T termination region of the 5S gene described in the literature [6, 35]. Interestingly, in addition to the poly-T region, the type I sequences showed a GAACAAA sequence very similar to the sequence GAAACAA, which is found downstream from the 5S rRNA gene in fish and has been suggested to act as a terminal region .
The consensus secondary structures for the presumed type I and type II 5S rRNA is shown in Figure 6. The secondary structure of all analyzed 5S rRNAs consists of five helices (I-V), two hairpin loops (C and E), two internal loops (B and D), and a hinge region (A), arranged into the three-helix junction.
Chromosome mapping of 5S rDNA
As expected, the karyotypes of the specimens of E. freibergi and E. petersi analyzed here were the same as those described previously by Targueta et al. . Therefore, our study includes two of the three karyological groups recognized by Targueta et al.  among E. petersi specimens: the Puyo and Yasuní karyological groups. It is interesting to notice that the karyotypes of the specimens of E. petersi from Puyo are more similar to those of the specimens of E. freibergi than to the karyotype of E. petersi from Yasuní. To describe the mapping of the 5S rDNA sequences in these karyotypes, the previously proposed chromosome classification scheme  was used.
In the E. freibergi karyotype, a FISH probe containing the entire repeating unit of the type I 5S rDNA hybridized pericentromerically to the short arm of Chromosome 3 and also distally to the long arm of Chromosome 6 (Figure 7a). When using the probe containing the entire repeating unit of the type II 5S rDNA, only the distal region of the long arm of Chromosome 6 was detected (Figure 7b). Neither the site detected in Chromosome 3 nor the site on Chromosome 6 coincided with any nucleolar organizer regions (NORs) reported by Targueta et al.  (Figure 7).
Because these data suggested a differential localization of both types of 5S rDNA sequences in the E. freibergi karyotype, further analyses were performed using more specific probes that exclusively contained either the type I or type II NTS. In the karyotypes of E. freibergi and E. petersi from Puyo, which are morphologically very similar (for details, see ), the type I NTS probe detected the proximal region of the short arm of Chromosome 3 but did not detect a distal region of Chromosome 6. This probe also detected the centromeric region of several chromosomes (Figure 8b). In contrast, the type II NTS probe detected only a distal site in the long arm of Chromosome 6 in the karyotypes of both E. freibergi (Figure 8a) and E. petersi from Puyo (Figure 8c). In the karyotype of E. petersi from Yasuní, the type I NTS probe detected the pericentromeric region of the short arm of Chromosome 5 and the centromeric regions of several chromosomes (Figure 8d), whereas the type II NTS probe detected only a distal site on the long arm of Chromosome 6 (Figure 8e).
Molecular Organization of the 5S rDNA in Engystomops
This study verified the occurrence of two types of 5S rDNA in the genomes of the Amazonian species of Engystomops, a feature widely documented in many vertebrates, including fish [10, 11, 36–42], the anurans Xenopus laevis and Xenopus borealis [6, 16, 21–24, 43, 44] chickens [45–47], and mammals [12, 48, 49]. Similarly to the findings reported in these studies, the two types of 5S rDNA sequences found in E. freibergi and E. petersi varied slightly in their corresponding coding regions, with the main difference between them found in the NTS region, which varied in length (84 bp for the type I 5S rDNA and approximately 650 bp for the type II 5S rDNA) and nucleotide composition. These two types of 5S rDNA do not, however, appear to be related to the dual system observed in Xenopus (i.e., the oocyte and somatic types) [6, 13, 16, 20–22, 24] because we found no similarities between the NTSs of the oocyte- or somatic-type sequences to those in the present study.
Since Ohno's publication , the origin of genic variants has been attributed to events of sequence duplication followed by processes that result in the divergence of the duplicated sequences. This hypothesis has been corroborated by many studies (reviewed in references [51–55]) and may also explain the presence of two types of 5S rDNA. Gene duplication may result from unequal crossing over, retropositioning, or chromosomal (or genomic) duplication , and the outcomes of these events are quite different, including neofunctionalization, pseudogene origin or simple preservation of gene duplicates [50, 56, 57].
Although our data for the Amazonian Engystomops do not allow us to elucidate the events involved in the origin of either type of the 5S rDNA, those events may have preceded the divergence of Engystomops and Physalaemus. Such an inference follows from the observation of higher nucleotide divergence between the sequences of the type I and type II 5S rDNAs of each species in these genera than among the sequences of the same type found in distinct species.
With respect to the functionality of the 5S rDNA sequences found in Engystomops, the analysis of the secondary structure of the presumed rRNAs shows that both types of sequences are consistent with the general eukaryotic 5S rRNA structure [58, 59], suggesting that the type I and type II 5S rDNA sequences may have transcriptional potential. The functionality of the type I 5S rDNA sequences of Engystomops is corroborated by the recognition, in this type of sequences, of elements that are similar to those considered to be important for the transcriptional activity of 5S rRNA genes. Those type I 5S rDNA elements are: (i) sequences quite similar to the ICR elements; (ii) a T-rich region downstream from the presumed coding region; (iii) a TATA-box located 25 nucleotides upstream from the coding region; and (iv) the nucleotide C at position -1. In addition, in the type I 5S rDNA of Engystomops, the presence of a GAACAAA segment was noted, which is very similar to the sequence GAAACAA suggested to act as a terminal region of the 5S gene transcription in fish . A region tentatively named the TATA-like region was also found in the Engystomops type I 5S rDNA, located 12 nucleotides upstream from the TATA-box. Campo et al.  reported an additional TATA-like region in the NTS regions of the 5S rDNA sequences of the fishes Merluccius merluccius, Merluccius senegalensis, and Merluccius capensis, and suggested that this TATA-box may serve as a "backup". The same hypothesis may be considered for the similar sequence found in the NTS of the type I 5S rDNA of Engystomops.
In contrast, some doubts remain about the transcriptional potential of the type II 5S rDNA sequences of Engystomops. In the type II repeats, the nucleotide at position -1 is a T and not a C; the ICR segment differed more from the ICRs of the other vertebrates used for comparison than the ICR of the type I repeats; and no TATA-box was found in the NTS. The only segment that resembles a TATA-like motif in the type II repeats of Engystomops was observed very distant from the region considered to be the coding region, approximately at position -420. However, a T-rich region downstream from the presumed coding region is also present in the Engystomops type II 5S rDNA.
It is also interesting to note that the results of previous experiments with Xenopus suggest that the oligonucleotides AGAAGC and AAAAGT, located at positions -28 to -23 and -18 to -13, respectively, may be involved in the initiation of 5S rDNA transcription instead of a TATA-box . In the type II 5S rDNA sequences of Engystomops, the hexanucleotides AGAAGC and GCAAGT were found at positions -53 to -48 and -21 to -16, respectively. The similarity of these oligonucleotides to those described for Xenopus, despite the low similarity of the remaining NTS sequence, is an interesting issue to be considered in further analyses of the functionality of the Engystomops type II 5S rDNA.
Physical Mapping of the 5S rDNA in the Engystomopskaryotypes X evolutionary diversification in 5S rDNA
The FISH assays suggest the existence of two sites of sequence accumulation of the 5S rDNA in the karyotypes of E. freibergi and E. petersi, one on 3p (E. freibergi and E. petersi from Puyo) or 5p (E. petersi from Yasuní) and another on 6q. The results of these assays also suggest that the former site is exclusive to or preferentially constitutes type I sequences, whereas the latter, on 6q, is associated with type II sequences.
Ribosomal DNA repeating units are evolutionarily dynamic and appear to be able to spread throughout the genome, creating new rDNA loci [3, 60–62]. The presence of two distinct 5S rDNA sequence types organized in different chromosomal regions or even on different chromosomes has been described for several fish, e.g., Salmo solar , Oncorhynchus mykiss , Coregonus artedi, C. zenithicus , and Oreochromis niloticus .
In anurans, the first chromosome mapping experiments for the 5S rRNA genes were conducted in Xenopus laevis. Using specific probes, Harper et al.  revealed a differential distribution of the two types of 5S rDNA in the Xenopus karyotypes, mapping the somatic-type 5S rDNA to the distal end of the long arm of Chromosome 9 in Xenopus laevis and X. borealis, and the oocyte-type to the distal ends of the majority of Xenopus laevis chromosomes. The authors also mapped a trace oocyte-type 5S rDNA in X. laevis, which is a minor class of the oocyte type, to the distal end of the long arm of Chromosome 13.
In addition to these data for the Xenopus karyotypes, only the 5S rDNA chromosomal sites were detected in the karyotypes of Physalaemus ephippifer  and Physalaemus cuvieri . A probe containing the entire repeat of the type I 5S rDNA of P. cuvieri detected a pericentromeric region of the short arm of Chromosome 3 in both Physalaemus species karyotypes [18, 19] whereas a probe with the entire repeat of the type II 5S rDNA of P. cuvieri preferentially detected a distal region of the long arm of Chromosome 6 . Similarly to the above-mentioned cases, a differential localization of the two types of 5S rDNA found in the Engystomops species was observed in this study.
As mentioned above, taking into account the similarity of the sequences, probably the origin of the two types of 5S rDNA found in Engystomops and Physalaemus species preceded the divergence of these genera. Apparently, the origin of this dual-system involved translocation or transposition events that lead to the separation of two groups of 5S rDNA sequences, favoring the dominance of divergence forces over homogenization processes between these groups. On the other hand, the homogeneity of the 5S rDNA repeats clustered in the same chromosomal site was maintained, what may be explained by concerted evolution. As a result of these processes, two distinct types of 5S rDNA sequences, occupying different chromosomal sites, have arisen. It is worth mentioning that purifying selection may also have been involved in this scenario. In addition to concerted evolution, purifying selection has been invoked to justify the homogenization in a gene family [64, 65]. Since the comparison between both types of 5S rDNA sequences of Engystomops showed a higher variation between the presumed coding-regions than between their NTSs, it is likely that purifying selection has been acting over these coding-regions, avoiding high level of divergence.
Another intriguing finding of this study was the hybridization of the probe that corresponds to the type I NTS to the centromeric region of various chromosomes, sites which were not detected in this analysis by the probes that potentially contain the transcribed region of the 5S rDNA. A possible explanation for this result is that the sequences associated with the centromeric regions are segments of satellite DNA derived from the 5S rDNA, a phenomenon previously reported for the fish Hoplias malabaricus  and the frog Physalaemus cuvieri .
The 5S rDNA cluster has been reported to be linked to the major rDNA sequences [4, 36, 37, 67–69] but, in some cases, is localized to different chromosomes [18, 63, 70–73]. The differential chromosomal localization of the 5S and 45S rDNAs is the prevalent condition, not only in the cited examples but also in other groups, including plants [74–77].
Several authors have previously discussed the prevalence of different chromosomal sites for the 45S and 5S rDNAs over their linkage in other organisms, and a probable explanation is intrinsically related to the repetitive nature of these sequences. Others have suggested that because tandem repeated sequences are frequently involved in events of unequal crossing-over and gene conversion, the separation of the two great families of ribosomal DNA at different chromosomal sites would avoid disruptive interference in its organization such as undesired rearrangements between the 45S and 5S arrays [3, 11].
The 5S rDNA chromosomal sites unlinked to the NORs of the Engystomops species represent new cytogenetic markers to be considered for their karyotypic comparison. Based on classic cytogenetic techniques, CMA3 and DAPI staining, and in situ localization of nucleolar rDNA, Targueta et al.  described the three Engystomops karyotypes of the present study and noted that recognizing chromosomal homeology was difficult, especially between the Yasuní karyotype and the other two. In addition, difficulties in differentiating the three very morphologically similar chromosome pairs (pairs 3, 6, and 8) in the Yasuní karyotype have been reported. In the present study, we were able to suggest a homeology between Chromosome 6 of the E. petersi karyotype from Yasuní and Chromosome 6 of E. freibergi and E. petersi from Puyo based on the mapping of the type II 5S rDNA sequences. Additionally, this chromosome site also constitutes a distinctive marker for Chromosome 6 in the karyotype of E. petersi from Yasuní, distinguishing it from Chromosome 8 and the NOR-bearing Chromosome 3. Finally, the mapping of the type I 5S rDNA sequence to Chromosome 5 of the specimens from Yasuní and Chromosome 3 of E. freibergi and of E. petersi from Puyo, which are similar chromosomes in size and morphology, may suggests that these chromosomes are homeologous.
In addition to the recognition of chromosomal homeologies among the Engystomops species, the 5S rDNA mapping performed here allows for a better cytogenetic comparison of Engystomops with its sister genus, Physalaemus. Chromosomes 3 of P. cuvieri  and P. ephippifer , bearing the type I 5S rDNA sequence, are morphologically similar to Chromosome 5 of E. petersi from Yasuní and Chromosome 3 of E. freibergi and E. petersi from Puyo. Therefore, the homeology of all these chromosomes can be strongly inferred. Similarly, we can deduce homeology among the metacentric Chromosome 6 of the Engystomops species and Chromosome 5 of P. cuvieri, which all carry the type II 5S rDNA sequences. Chromosome 5 or 6 of P. ephippifer likely also bears this type of sequence; however, the presence of this sequence has not been verified with a specific probe for type II 5S rDNA sequences . These cytogenetic data suggest that the chromosomal sites of the 5S rDNA may be conserved in these leiuperid genera, which has not been observed for the NORs [18, 33, 78–84]. Therefore, the two 5S rDNA arrays appear to be independent units of evolution in Engystomops species, and further studies of their functionality and their relation to a possible centromeric DNA satellite sequence are necessary to provide a better understanding of the evolution of these sequences.
In Amazonian Engystomops, two types of 5S rDNA were found and mapped to distinct chromosomal sites. Because these rDNA types are related to those found in Physalaemus with respect to their nucleotide sequences and chromosomal locations, their origin likely preceded the evolutionary divergence of these genera. In addition, our data revealed chromosomal homeologies among the three karyotypes of the Amazonian Engystomops, representing an important contribution for further studies of karyotype evolution in this genus.
Specimens of E. petersi from Puyo (Provincia de Pastaza) and Estación Científica de Yasuní (Provincia de Olleana), located in the Ecuadorian Amazon, and specimens of E. freibergi from the Tejo estuary, state of Acre, Brazil were studied. The vouchers specimens are at the Museo de Zoología de la Pontificia Universidad Católica del Ecuador (QCAZ), Quito, Ecuador or at the Museu de Zoologia "Prof. Adão José Cardoso" (ZUEC) at the Universidade Estadual de Campinas (UNICAMP), Campinas, state of São Paulo, Brazil.
Isolation and cloning of the 5S gene
Genomic DNA was isolated from liver and muscle samples stored at -80°C in the tissue collection deposited at the Laboratório de Estudos Cromossômicos em Anuros at IB-UNICAMP, Brazil. The genomic DNA isolation was performed as described , with a TNES solution (250 mM Tris-HCl pH 7.5, 2 M NaCl, 100 mM EDTA pH 8.0, and 2.5% SDS). After electrophoresis in 1% agarose, the DNA quality was evaluated, and its quantity was estimated. The entire repeating unit of the 5S rDNA (which includes the presumed transcribed region and the NTS) was isolated by PCR using the primers 5S-A (5'-TACGCCCGATCTCGTCCGATC-3') and 5S-B (5'-CAGGCTGGTATGGCCGTAAGC-3') . The amplified fragments were purified using the GFX PCR and Gel Band DNA Purification kit (GE Healthcare - Little Chalfont, Buckinghamsire, UK) and cloned into the pGEM-T vector (pGEM-T easy Vector - Promega - Madison, WI, USA) according to the manufacturer's instructions. Recombinant vectors were used to transform competent Escherichia coli of the JM109 strain (Fermentas), and the cloned fragments were sequenced. Repeating units of two different sizes were obtained and named type I 5S rDNA and type II 5S rDNA. Based on their nucleotide sequences, specific primers were designed to exclusively isolate the presumed 120-bp transcribed region from the type I 5S rDNA (5S120T1-F: 5'-GCTTTCGCTTACGGCCATACC-3'; 5S120T1-R: 5'-AGCTTACAGCACCTGGTATTC-3') and the type II 5S rDNA (5S120T2-F: 5'-GTCTCTGCTTACGGCCATACC-3'; 5S120T2-R: 5'-AGTCTACAGCACCCGCGCTTC-3'). Specific primers were also designed to isolate the NTS region from the type I (5ST1-F: 5'-GCTGTAAGCTTTTTGTTTTTTGAA-3'; 5ST1-R: 5'-GAAAGCTCAGGGCCTGTACAG-3') and type II (5ST2-F: 5'-GCTGTAGACCGTTTATTTACCTT-3'; 5ST2-R: 5'-AGAGACAGGCCTCTCACTTGC-3') 5S rDNA repeats. The annealing sites of the primers used in this study are indicated in Figure 1. The resulting amplified fragments were all sequenced and analyzed.
DNA sequencing and analysis
The cloned fragments were amplified by PCR using the T7 and SP6 primers. The amplified products were purified using the GFX PCR and Gel Band DNA Purification kit (GE Healthcare - Little Chalfont, Buckinghamsire, UK) and used directly as templates in amplification reactions using the BigDye Terminator chemistry version 3.1 (Applied Biosystems - Austin, TX, USA) according to the manufacturer's recommendations. Each cloned fragment was bi-directionally sequenced in an automatic DNA sequencer. The sequences were edited using BioEdit version 7.0.1  and aligned with ClustalW software. The fragments obtained were compared with other sequences available in the GenBank database  and in the 5S Ribosomal RNA Database . All presumed coding regions of the generated 5S rDNA sequences were assembled into a data matrix, and a Maximum Likelihood (ML) analysis was conducted using the PAUP* 4.0β10 software  with the evolution model K80, which was selected by Modeltest 3.7  for this data set. Nodal support for the ML arrangement was assessed through a non-parametric bootstrap analysis , with a heuristic search based on 1,000 pseudoreplicates.
Chromosome preparations and Fluorescent in situHybridization (FISH)
Chromosome preparations were made from the intestine and testis cell suspensions obtained from the E. freibergi specimens ZUEC 14435, ZUEC 14439, ZUEC 14440, and ZUEC14458, and the E. petersi specimens QCAZ 34946, QCAZ 34947, QCAZ 34940, and QCAZ 34942. All of these cell suspensions were available at the collection of amphibian material deposited at the Laboratório de Estudos Cromossômicos em Anuros at IB-UNICAMP, Brazil, and were used previously by Targueta et al.  to describe the karyotypes of E. freibergi and E. petersi. The cell suspensions were spotted onto clean slides, and the chromosome preparations were hybridized with the 5S rDNA fragments isolated as described above and PCR-labeled with either biotin or digoxigenin. The hybridization was performed as described elsewhere . The biotin-labeled probes were detected with an anti-biotin antibody (goat anti-biotin - Vector - Burlingame CA, USA), which is recognized by a FITC-conjugated secondary antibody (anti-goat IgG-FITC - Vector - Burlingame CA, USA). The digoxigenin-labeled probes were detected with an anti-digoxigenin antibody conjugated with rhodamine. The chromosomes were counterstained with DAPI (0.5 μg/mL).
short arm of Chromosome 3
short arm of Chromosome 5
long arm of Chromosome 6.
Long EO, Dawid IB: Repeated genes in eukaryotes. Ann Rev Biochem. 1980, 49: 727-764. 10.1146/annurev.bi.49.070180.003455.
Drouin G, de Sá MM: The concerted evolution of 5S ribosomal genes linked to the repeat units of other multigene families. Mol Biol Evol. 1995, 12 (3): 481-493.
Martins C, Wasko AP: Organization and evolution of 5S ribosomal DNA in the fish genome. Focus on Genome Research Chapter X. Edited by: Clyde R. 2004, Williams: Nova Science Publishers, New York, USA, Inc, 335-363.
Gerbi SA: Evolution of Ribosomal DNA. Molecular Evolutionary Genetics. Edited by: MacIntyre RJ. 1985, New York: Plenum Publishing, 419-517.
Pieler T, Hamm J, Roeder RG: The 5S gene internal control region is composed of three distinct sequence elements, organized as two functional domains with variable spacing. Cell. 1987, 48 (1): 91-100. 10.1016/0092-8674(87)90359-X.
Korn LJ, Brown DD: Nucleotide sequence of Xenopus borealis oocyte 5S DNA: Comparison of sequences that flank several related eukaryotic genes. Cell. 1978, 15: 1145-1156. 10.1016/0092-8674(78)90042-9.
Hallenberg C, Frederiksen S: Effect of mutations in the upstream promoter on the transcription of human 5S rRNA genes. Biochim Biophys Acta. 2001, 1520 (2): 169-173.
Nederby-Nielsen J, Hallenberg C, Frederiksen S, Sorensen PD, Lomholt B: Transcription of human 5S rRNA genes is influenced by an upstream DNA sequence. Nucleic Acids Res. 1993, 21 (16): 3631-3636. 10.1093/nar/21.16.3631.
Fernández-Tajes J, Méndez J: Identification of the razor clam species Ensis arcuatus, E. siliqua, E. directus, E. macha, and Solen marginatus using PCR-RFLP analysis of the 5S rDNA region. J Agric Food Chem. 2007, 55: 7278-7282. 10.1021/jf0709855.
Ferreira IA, Oliveira C, Venere PC, Galetti PM, Martins C: 5S rDNA variation and its phylogenetic inference in the genus Leporinus (Characiformes: Anostomidae). Genetica. 2007, 129: 253-257. 10.1007/s10709-006-0005-6.
Martins C, Galetti PM: Chromosomal localization of 5S rDNA genes in Leporinus fish (Anostomidae, Characiformes). Chrom Res. 1999, 7: 363-367. 10.1023/A:1009216030316.
Suzuki H, Moriwaki K, Sakurai S: Sequences and evolutionary analysis of mouse 5S rDNAs. Mol Biol Evol. 1994, 11 (4): 704-710.
Vitelli L, Batistoni R, Andronico F, Nardi I, Barsacchipilone G: Chromosomal localization of 18S + 28S and 5S ribosomal-RNA genes in evolutionarily diverse anuran amphibians. Chromosoma. 1982, 84 (4): 475-491. 10.1007/BF00292849.
Pasolini P, Costagliola D, Rocco L, Tinti F: Molecular Organization of 5S rDNAs in Rajidae (Chondrichthyes): Structural Features and Evolution of Piscine 5S rRNA Genes and Nontranscribed Intergenic Spacers. J Mol Evol. 2006, 62: 564-574. 10.1007/s00239-005-0118-z.
Pinhal D1, Araki CS, Gadig OBF, Martins C: Molecular organization of 5S rDNA in sharks of the genus Rhizoprionodon: insights into the evolutionary dynamics of 5S rDNA in vertebrate genomes. Genet Res. 2009, 91: 61-72. 10.1017/S0016672308009993.
Harper ME, Price J, Korn LJ: Chromosonal mapping of Xenopus 5S genes: somatic-type versus oocyte-type. Nucleic Acids Res. 1983, 11 (8): 2313-2323. 10.1093/nar/11.8.2313.
Nietfeld W, Digweedt M, Mentzel H, Meyerhof W, Köster M, Knöchel W, Erdmannt VA, Pieler T: Oocyte and somatic 5S ribosomal RNA and 5S RNA encoding genes in Xenopus tropicalis. Nucleic Acids Res. 1988, 16: 8803-8815. 10.1093/nar/16.18.8803.
Nascimento J, Quinderé YRSD, Recco-Pimentel SM, Lima JRF, Lourenço LB: Heteromorphic Z and W sex chromosomes in Physalaemus ephippifer (Steindachner, 1864) (Anura, Leiuperidae). Gentica. 2010, 138: 1127-1132. 10.1007/s10709-010-9501-9.
Vittorazzi SE, Lourenço LB, Del-Grande ML, Recco-Pimentel SM: Satellite DNA Derived from 5S rDNA in Physalaemus cuvieri (Anura, Leiuperidae). Cytogenet Genome Res. 2011, 134 (2): 101-107. 10.1159/000325540.
Schmid M, Vitelli I, Batistoni R: Chromosome Banding in Amphibia XI. Constitutive heterochromatin, nucleolus organizers, 18S+28S and 5S ribossomal RNA genes in Ascaphidae, Pipidae, Discoglossidae and Pelobatidae. Chromosoma. 1987, 95: 271-284. 10.1007/BF00294784.
Carroll D, Brown DD: Repeating units of Xenopus laevis oocyte-type 5S DNA are heterogeneous in length. Cell. 1976, 7: 467-475. 10.1016/0092-8674(76)90198-7.
Fedoroff NV, Brown DD: The nucleotide sequence of oocyte 5S DNA in Xenopus laevis. I. The AT-rich spacer. Cell. 1978, 13: 701-716. 10.1016/0092-8674(78)90220-9.
Miller JR, Cartwright EM, Brownlee GG, Fedoroff NV, Brown DD: The nucleotide sequence of oocyte 5S DNA in Xenopus laevis. II. The GC-rich region. Cell. 1978, 13: 717-725. 10.1016/0092-8674(78)90221-0.
Peterson RC, Doering JL, Brown DD: Characterization of two Xenopus somatic 5S DNAs and one minor oocyte-specific 5S DNA. Cell. 1980, 20: 131-141. 10.1016/0092-8674(80)90241-X.
Wolffe AP, Brown DD: Developmental regulation of two 5S ribosomal RNA genes. Science. 1988, 241 (4873): 1626-1632. 10.1126/science.3420414.
Chipev CC, Wolffe AP: Chromosomal Organization of Xenopus laevis Oocyte and Somatic 5S rRNA Genes In Vivo. Mol Cell Biol. 1992, 12 (1): 45-55.
del Pino EM, Murphy C, Masson PH, Gall JG: 5S rRNA-encoding genes of the marsupial frog Gastrotheca riobambae. Gene. 1992, 111: 235-238. 10.1016/0378-1119(92)90692-I.
Gaskins CJ, Smith JF, Ogilvie MK, Hanas JS: Comparison of the sequence and structure of Transcription Factor-IIIA from Bufo americanus and Rana pipiens. Gene. 1992, 120 (2): 197-206. 10.1016/0378-1119(92)90094-6.
Nascimento LB, Caramaschi U, Cruz CAG: Taxonomic review of the species group of the genus Physalameus Fitzinger, 1826 with revalidation of the genera Engystomops Jimenez-de-la-Espada, 1872 and Eupemphix Steindachner, 1836 (Amphibia, Anura, Leptodactylidae). Arquivos do Museu Nacional do Rio de Janeiro. 2005, 63: 297-320.
Frost DR: Amphibian Species of the World: An online reference. [http://research.amnh.org/vz/herpetology/amphibia/]
Funk CW, Caldwell JP, Peden CE, Padial JM, la Riva ID, Cannatella DC: Tests of biogeographic hypotheses for diversiWcation in the Amazonian forest frog, Physalaemus petersi. Mol Phylogen Evol. 2007, 44 (2): 825-837. 10.1016/j.ympev.2007.01.012.
Ron SR, Guerra MA: Mate choice and courtship signal differentiation promotes speciation in an Amazonian frog. Behav Ecol. 2008, 19 (6): 1128-1135. 10.1093/beheco/arn098.
Targueta CP, Rivera M, Souza MB, Recco-Pimentel SM, Lourenço LB: Cytogenetic contributions for the study of the Amazonian Engystomops (Anura; Leiuperidae) assessed in the light of phylogenetic relationships. Mol Phylogenet Evol. 2010, 54 (3): 709-725. 10.1016/j.ympev.2009.10.018.
Lourenço LB, Recco-Pimentel SM, Cardoso AJ: Two karyotypes and heteromorphic sex chromosomes in Physalaemus petersi (Anura, Leptodactylidae). Can J Zool. 1999, 77 (4): 624-631.
Bogenhaggen DF, Brown DD: Nucleotide sequences in Xenopus 5S DNA required for transcription termination. Cell. 1981, 24 (1): 261-270. 10.1016/0092-8674(81)90522-5.
Pendás AM, Moran P, Freije JP, Garcia-Vazquez E: Chromosomal mapping and nucleotide sequence of two tandem repeats of Atlantic salmon 5S rDNA. Cytogenet Cell Genet. 1994, 67: 31-36. 10.1159/000133792.
Móran P, Martinez JL, Garcia-Vazquez E, Pendás AM: Sex chromosome linkage of 5S rDNA in rainbow trout (Oncorhynchus mykiss). Cytogenet Cell Genet. 1996, 75 (2-3): 145-150. 10.1159/000134466.
Martins C, Galetti PM: Organization of 5S rDNA in species of the fish Leporinus: two different genomic locations are characterized by distinct nontranscribed spacers. Genome. 2001, 44: 903-910.
Martins C, Galetti PM: Two 5S rDNA arrays in Neotropical fish species: is it a general rule for fishes?. Genetica. 2001, 111: 439-446. 10.1023/A:1013799516717.
Martins C, Wasko AP, Oliveira C, Porto-Foresti F, Parise-Maltempi PP, Wrightb JM, Foresti F: Dynamics of 5S rDNA in the tilapia (Oreochromis niloticus) genome: repeat units, inverted sequences, pseudogenes and chromosome loci. Cytogenet Genome Res. 2002, 98: 78-85. 10.1159/000068542.
Alves-Costa FA, Wasko AP, Oliveira C, Foresti F, Martins C: Genomic organization and evolution of the 5S ribosomal DNA in Tilapiini fishes. Genetica. 2006, 127: 243-252. 10.1007/s10709-005-4013-8.
Campo D, Machado-Schiaffino G, Horreo JL, Garcia-Vazquez E: Molecular organization and evolution of 5S rDNA in the genus Merluccius and their phylogenetic implications. J Mol Evol. 2009, 68: 208-216. 10.1007/s00239-009-9207-8.
Pardue ML, Brown DD, Birnstiel ML: Location of the genes for 5S ribosomal RNA in Xenopus laevis. Chromosoma. 1973, 42: 191-203. 10.1007/BF00320940.
Komiya H, Hasegawa M, Takemura S: Differentiation of oocyte- and somatic-type 5S rRNAs in animals. J Biochem. 1986, 100: 369-374.
Brownlee GG, Cartwright EM: The nucleotide sequence of the 5S RNA of chicken embryo fibroblasts. Nucleic Acids Res. 1975, 2: 2279-2288. 10.1093/nar/2.12.2279.
Lazar E, Haendler B, Jacob M: Two 5S genes are expressed in chicken somatk cells. Nucleic Acids Res. 1983, 11 (22): 7735-7741. 10.1093/nar/11.22.7735.
Daniels LM, Delany ME: Molecular and cytogenetic organization of the 5S ribosomal DNA array in chicken (Gallus gallus). Chrom Res. 2003, 11: 305-317. 10.1023/A:1024008522122.
Hart RP, Folk WR: Structure and Organization of a Mammalian 5 S Gene. Cluster J Biol Chem. 1982, 257 (19): 11706-11711.
Jensen LR, Frederiksen S: The 5S rRNA genes in Macaca fascicularis are organized in two large tandem repeats. Biochim Biophys Acta. 2000, 1492 (2-3): 537-542.
Ohno S: Evolution by Gene Duplication. 1970, New York, Heidelberg, Berlin: Springer-Verlag
Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18 (6): 292-298. 10.1016/S0169-5347(03)00033-8.
Lynch M: The Origins of Genome Architecture. 2007, Indiana University: Sinauer Associates Inc
Zhou Q, Wang W: On the origin and evolution of new genes - A genomic and experimental perspective. J Genet Genomics. 2008, 35: 639-648. 10.1016/S1673-8527(08)60085-5.
Van de Peer Y, Maere S, Meyer A: The evolutionary significance of ancient genome duplications. Nat Rev Genet. 2009, 10: 725-732. 10.1038/nrg2600.
Kaessmann H, Vinckenbosch N, Long M: RNA-based gene duplication: mechanistic and evolutionary insights. Nat Rev Genet. 2009, 10 (1): 19-31.
Ohno S: So much 'junk DNA' in our genome. Brookhaven Symp Biol. 1972, 23: 366-370.
Kaessmann H: Origins, evolution, and phenotypic impact of new genes. Genome Res. 2010, 20: 1313-1326. 10.1101/gr.101386.109.
Delihas N, Andersen J: Generalized structures of the 5S ribosomal RNAs. Nucl Acids Res. 1982, 10: 7323-7344. 10.1093/nar/10.22.7323.
Barciszewska MZ, Szymanski M, Erdmann VA, Barciszewski J: 5S ribosomal RNA. Biomacromolecules. 2000, 1: 297-302. 10.1021/bm000293o.
Iborra FJ, Cook PR: The interdependence of nuclear structure and function. Curr Opin Cell Biol. 2002, 14: 780-785. 10.1016/S0955-0674(02)00386-1.
Raskina O, Barber JC, Nevo E, Belyayev A: Repetitive DNA and chromosomal rearrangements: speciation-related events in plant genomes. Chrom Evol. 2008, 120: 351-357.
Zhang X, Eickbush MT, Eickbush TH: Role of recombination in the long-term retention of transposable elements in rRNA gene loci. Genetics. 2008, 180: 1617-1626. 10.1534/genetics.108.093716.
Sajdak SL, Reed KM, Phillips RB: Intraindividual and interspecies variation in the 5S rDNA of Coregonid fish. J Mol Evol. 1998, 46: 680-688. 10.1007/PL00006348.
Nei M, Rogozin IB, Piontkivska H: Purifying selection and birth-and-death evolution in the ubiquitin gene family. PNAS. 2000, 97 (20): 10866-10871.
Pinhal D, Yoshimura TS, Araki CS, Martins C: The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays. BMC Evol Biol. 2011, 11: 151-165. 10.1186/1471-2148-11-151.
Martins C, Ferreira IA, Oliveira C, Foresti F, Galetti PM: A tandemly repetitive centromeric DNA sequence of the fish Hoplias malabaricus (Characiformes: Erythrinidae) is derived from 5S rDNA. Genetica. 2006, 127: 133-141. 10.1007/s10709-005-2674-y.
Fujiwara A, Abe S, Yamaha E, Yamazaki F, Yoshida MC: Chromosomal localization and heterochromatin association of ribosomal RNA gene loci and silver-stained nucleolar organizer regions in salmonid fishes. Chrom Res. 1998, 6: 463-471. 10.1023/A:1009200428369.
Dobigny G, Ozouf-Costaz C, Bonillo C, Volobouev V: Evolution of rRNA gene clusters and telomeric repeats during explosive genome repatterning in Taterillus X (Rodentia, Gerbillinae). Cytogenet Genome Res. 2003, 103: 94-103. 10.1159/000076296.
Cabral-de-Mello DC, Oliveira SG, Moura RC, Martins C: Chromosomal organization of the 18S and 5S rRNAs and histone H3 genes in Scarabaeinae coleopterans: insights into the evolutionary dynamics of multigene families and heterochromatin. BMC Genet. 2011, 12: 88-100.
Martinez JL, Morán P, Garcia-Vazquez E, Pendás AM: Chromosomal localization of the major and 5S rRNA genes in the European eel (Anguilla anguilla). Cytogen Cell Genet. 1996, 73: 149-152. 10.1159/000134328.
Zijlstra C, Mellink CHM, deHaan NA, Bosma AA: Localization of the 18S, 5.8S and 28S rRNA genes and the 5S rRNA genes in the babirusa and the white-lipped peccary. Cytogenet Cell Genet. 1997, 77: 273-277. 10.1159/000134596.
Gornung E, Innocentiis SD, Annesi F, Sola L: Zebrafish 5S rRNA genes map to the long arms of chromosome 3. Chrom Res. 2008, 8: 362-
Srikulnath K, Matsubara K, Uno Y, Thongpan A, Suputtitada S, Apisitwanich S, Matsuda Y, Nishida C: Karyological characterization of the butterfly lizard (Leiolepis reevesii rubritaeniata, Agamidae, Squamata) by molecular cytogenetic approach. Cytogenet Gen Res. 2009, 125 (3): 213-223. 10.1159/000230005.
Li X, Guo R, Pedersen C, Hayman D, Langridge P: Physical localization of rRNA genes by two-colour fluorescent in-situ hybridization and sequence analysis of the 5s rRNA gene in Phalaris coerulescens. Hereditas. 1997, 126: 289-294.
Ansari HA, Ellison NW, Reader SM, Badaeva ED, Friebe B, Miller TE, Williams WE: Molecular cytogenetic organization of 5S and 18S-26S rDNA loci in white clover (Trifolium repens L.) and Related Species. Ann Bot. 1999, 83: 199-206. 10.1006/anbo.1998.0806.
Siljak-Yakovlev S, Cerbah M, Coulaud J, Stoian V, Brown SC, Zoldos V, Jelenic S, Papes D: Nuclear DNA content, base composition, heterochromatin and rDNA in Picea omorika and Picea abies. Theor Appl Genet. 2002, 104: 505-512. 10.1007/s001220100755.
Puizina J, Weiss-Schneeweiss H, Pedrosa-Harand A, Kamenjarin J, Trinajstic I, Riha K, Schweizer D: Karyotype analysis in Hyacinthella dalmatica (Hyacinthaceae) reveals vertebrate-type telomere repeats at the chromosome ends. Genome. 2003, 46 (6): 1070-1076. 10.1139/g03-078.
Lourenço LB, Recco-Pimentel SM, Cardoso AJ: Polymorphism of the nucleolus organizer regions (NORs) in Physalaemus petersi (Amphibia, Anura, Leptodactylidae) detected by silver staining and fluorescence in situ hybridization. Chrom Res. 1998, 6: 621-628. 10.1023/A:1009253410553.
Silva APZ, Haddad CFB, Kasahara S: Nucleolus organizer regions in Physalaemus cuvievi (Anura, Leptodactylidae), with evidence of a unique case of Ag-NOR variability. Hereditas. 1999, 131: 135-141.
Silva APZ, Baldissera FA, Haddad CFB, Kasahara S: Karyotypes and nucleolar organizer regions in four species of the genus Physalaemus (Anura, Leptodactylidae). IheringiaSér Zool. 2000, 88: 159-164.
Amaral MJLV, Cardoso AJ, Recco-Pimentel SM: Cytogenetic analysis of three Physalaemus species (Amphibia, Anura). Caryologia. 2000, 53 (2): 283-288.
Quinderé YRSD, Lourenço LB, Andrade GV, Tomatis C, Baldo D, Recco-Pimentel SM: Polytypic and polymorphic cytogenetic variations in the widespread anuran Physalaemus cuvieri (Anura, Leiuperidae) with emphasis on nucleolar organizing regions. Biol Res. 2009, 42: 79-92.
Tomatis C, Baldo D, Kolenc F, Borteiro C: Chromosomal Variation in the Species of the Physalaemus henselii Group (Anura: Leiuperidae). J Herpet. 2009, 43: 555-560. 10.1670/08-122R1.1.
Milani M, Cassini CS, Recco-Pimentel SM, Lourenço LB: Karyotypic data detect interpopulational variation in Physalaemus olfersiiand the first case of supernumerary chromosome in the genus. Animal Biol J. 2011, 2 (1): 21-28.
Hall T: Bioedit. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html]
Szymanski M, Barciszewska MZ, Erdmann VA, Barciszewski J: 5S Ribosomal RNA Database. [http://rose.man.poznan.pl/5SData]
Swofford D: PAUP 4.0β10. [http://paup.csit.fsu.edu/index.html]
Posada D: [http://darwin.uvigo.es/software/modeltest.html]
Viegas-Péquignot E: In situ Hybridization to chromosomes biotinylated probes. situ Hybridization: A practical approach. Edited by: Wilkinson DG. 1992, The Practical Approaches Series, 137-157.
Acknowledgements and funding
The authors gratefully acknowledge Ailín Blasco for helping in the field work and Shirlei M. Recco Pimentel for the constant encouragement and assistance in obtaining some tissue samples. This work was supported by the Brazilian agencies FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo), CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) and CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), and the Ecuadorian agency SENESCYT (Secretaría Nacional de Educación Superior, Ciencia, Tecnología y Innovación).
The authors declare that they have no competing interests.
DSR acquired the data and drafted the manuscript. MR helped to collect the specimens and obtain the chromosome preparations and revised the manuscript. LBL designed and coordinated the study and helped draft the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.