Molecular organization and chromosomal localization of 5S rDNA in Amazonian Engystomops (Anura, Leiuperidae)

Evidence for the Transcription of a Satellite DNA Widely Found in Frogs

BACKGROUND

The satellite DNA (satDNA) PcP190 has been identified in multiple frog species from seven phylogenetically distant families within Hyloidea, indicating its broad distribution. This satDNA consists of repeats of approximately 190 bp and exhibits a highly conserved region (CR) of 120 bp, which is similar to the transcribed region of 5S ribosomal DNA (rDNA), and a hypervariable region (HR) that varies in size and nucleotide composition among and within species. Here, to improve our understanding of PcP190 satDNA, we searched for evidence of its transcription in the available transcriptomes of Rhinella marina (Bufonidae) and Engystomops pustulosus (Leptodactylidae), two phylogenetically distantly related species.

METHODS

We first characterized the 5S rDNA and PcP190 sequences in these species by searching for them in available genome assemblies. Next, we used the PcP190 (CR and HR) and 5S rDNA sequences of each species as queries to search for these sequences in RNA-seq libraries.

RESULTS

We identified two types of 5S rDNA in each analyzed species, with a new type found in E. pustulosus. Our results also revealed a novel type of PcP190 sequence in R. marina and a new subtype of PcP-1 in E. pustulosus. Transcriptome analyses confirmed the expected transcription of the 5S rRNA gene and showed transcription of both the CR and HR of the PcP190 satDNA in both species and in different tissues.

CONCLUSIONS

As the entire repeat of this satDNA is susceptible to transcription, the high variability observed in the HR cannot be attributed to transcriptional activity confined to the CR.

Comparative cytogenetics among Boana species (Anura, Hylidae): focus on evolutionary variability of repetitive DNA

Boana comprises a diverse genus of Neotropical treefrogs, currently rearranged into seven taxonomic species groups. Although cytogenetic studies have demonstrated diversity in its representatives, the chromosomal mapping of repetitive DNA sequences is still scarce. In this study, Boana albopunctata, Boana faber, and Boana prasina were subjected to in situ localization of different repetitive DNA units to evaluate trends of chromosomal evolution in this genus. Boana faber and B. prasina had 2n=24 chromosomes, while B. albopunctata has 2n=22 and an intra-individual variation related to the presence/absence of one B chromosome. The location of 45S rDNA sites was different in the analyzed karyotypes, corroborating with what was found in the distinct phylogenetic groups of Boana. We presented the first description of 5S rDNA in a Boana species, which showed markings resulting from transposition/translocation mechanisms. In situ localization of microsatellite loci proved to be a helpful marker for karyotype comparison in Boana, commonly with cis accumulation in the heterochromatin. On the other hand, genomic dispersion of microsatellites may be associated with hitchhiking effects during the spreading of transposable elements. The obtained results corroborated the independent diversification of these lineages of species from three distinct phylogenetic groups of Boana.

Cytogenetic characterization and mapping of the repetitive DNAs in Cycloramphus bolitoglossus (Werner, 1897): More clues for the chromosome evolution in the genus Cycloramphus (Anura, Cycloramphidae)

Cycloramphus bolitoglossus (Werner, 1897) is a rare species with a low population density in the Serra do Mar region of Paraná and Santa Catarina, in southern Brazil. Currently, it has been assigned to the Near Threatened (NT) category in the Brazilian List of Endangered Animal Species. Here, we described the karyotype of this species for the first time and investigated the patterns of some repetitive DNA classes in the chromosomes using molecular cytogenetic approaches. We isolated, sequenced and mapped the 5S rDNA and the satellite DNA PcP190 of C. bolitoglossus, as well as mapped the telomeric sequences and seven microsatellites motifies [(GA)15, (CA)15, (GACA)4, (GATA)8, (CAG)10, (CGC)10, and (GAA)]10. Cycloramphus bolitoglossus has 2n = 26 chromosomes and a fundamental number (FN) equal to 52, with a highly conserved karyotype compared to other genus members. Comparative cytogenetic under the phylogenetic context of genus allowed evolutionary interpretations of the morphological changes in the homologs of pairs 1, 3, and 6 along with the evolutionary history of Cycloramphus. Two subtypes of 5S rDNA type II were isolated in C. bolitoglossus genome, and several comparative analysis suggests mixed effects of concerted and birth-and-death evolution acting in this repetitive DNA. The 5S rDNA II subtype "a" and "b" was mapped on chromosome 1. However, their different position along chromosome 1 provide an excellent chromosome marker for future studies. PcP190 satellite DNA, already reported for species of the families Hylidae, Hylodidae, Leptodactylidae, and Odontophrynidae, is scattered throughout the C. bolitoglossus genome, and even non-heterochromatic regions showed hybridization signals using the PcP190 probe. Molecular analysis suggests that PcP190 satellite DNA exhibit a high-level of homogenization of this sequence in the genome of C. bolitoglossus. The PcP190 satDNA from C. bolitoglossus represents a novel sequence group, compared to other anurans, based on its hypervariable region. Overall, the present data on repetitive DNA sequences showed pseudogenization evidence and corroborated the hypothesis of the emergence of satDNA from rDNA 5S clusters. These two arguments that reinforced the importance of the birth-and-death evolutionary model to explain 5S rDNA patterns found in anuran genomes.

A Comparative Study of 5S rDNA Non-Transcribed Spacers in Elaeagnaceae Species

5S rDNA is organized as a cluster of tandemly repeated monomers that consist of the conservative 120 bp coding part and non-transcribed spacers (NTSs) with different lengths and sequences among different species. The polymorphism in the 5S rDNA NTSs of closely related species is interesting for phylogenetic and evolutional investigations, as well as for the development of molecular markers. In this study, the 5S rDNA NTSs were amplified with universal 5S1/5S2 primers in some species of the Elaeagnaceae Adans. family. The polymerase chain reaction (PCR) products of five Elaeagnus species had similar lengths near 310 bp and were different from Shepherdia canadensis (L.) Nutt. and Sh. argentea (Pusch.) Nutt. samples (260 bp and 215 bp, respectively). The PCR products were cloned and sequenced. An analysis of the sequences revealed that intraspecific levels of NTS identity are high (approximately 95–96%) and similar in the Elaeagnus L. species. In Sh. argentea, this level was slightly lower due to the differences in the poly-T region. Moreover, the intergeneric and intervarietal NTS identity levels were studied and compared. Significant differences between species (except E. multiflora Thunb. and E. umbellata Thunb.) and genera were found. Herein, a range of the NTS features is discussed. This study is another step in the investigation of the molecular evolution of Elaeagnaceae and may be useful for the development of species-specific DNA markers in this family.

Cytogenetic analysis of the genus Thoropa Cope, 1865 (Anura-Cycloramphidae) with evolutionary inferences based on repetitive sequences

Cytogenetics can be a useful tool to assist in taxonomic problems by adding information to the widely used morphological and molecular approaches. These taxonomic problems are especially common in anurans, once they are very diverse, highly polymorphic, and present many cryptic species. The genus Thoropa Cope, 1865 is composed of six specialist species that reproduce in rocky outcrops and are distributed throughout the Atlantic Forest and Cerrado ecotones. Phylogenetic studies point to possible cryptic species within the T. miliaris group. To assist in the evolutionary and taxonomic understanding of this group, classical cytogenetic techniques were used to find possible molecular markers for the genus through rDNA5S, rDNA18S, and U2snDNA probes and analyze their chromosome distribution in the group of T. miliaris. Despite the well conserved karyotype under conventional staining and classical techniques, such as Ag-NOR, our C-banding results showed differences in the centromeric heterochromatin concentration between two populations of T. miliaris. Furthermore, some differences among the populations and species were found for rDNA5S and U2snDNA. This study contributes to a better understanding of the evolutionary relationships within the genus; however, the use of different probe sequences, such as satDNA, is essential for a more robust cytogenetic analysis.

Molecular Organization and Chromosomal Localization Analysis of 5S rDNA Clusters in Autotetraploids Derived From Carassius auratus Red Var. (♀) × Megalobrama amblycephala (♂)

The autotetraploid fish (4n = 200, RRRR) (abbreviated as 4nRR) resulted from the whole genome duplication of red crucian carp (Carassius auratus red var., 2n = 100, RR) (abbreviated as RCC). During investigation of the influence of polyploidization on organization and evolution of the multigene family of 5S rDNA, molecular organization and chromosomal localization of the 5S rDNA were characterized in autotetraploid fish. By sequence analysis of the coding region (5S) and adjacent non-transcribed spacer (NTS), three distinct 5S rDNA units (type I: 203 bp; type II: 340 bp; and type III: 477bp) were identified and characterized in 4nRR. These 5S rDNA units were inherited from their female parent (RCC), in which obvious base variations in NTS and array recombination of repeat units were found. Using fluorescence in situ hybridization employing different 5S rDNA units as probes, these 5S rDNA clusters were localized in chromosomes of 4nRR, respectively, and showed obvious loss of chromosomal loci (type I and type II). Our data revealed genetic variation of the 5S rDNA multigene family in the genome of autopolyploid fish. Furthermore, results provided new insights into the evolutionary patterns of this vertebrate multigene family.

Isolation and characterization of 5S rDNA sequences in catfishes genome (Heptapteridae and Pseudopimelodidae): perspectives for rDNA studies in fish by C0t method

Sequences of 5S ribosomal RNA (rRNA) are extensively used in fish cytogenomic studies, once they have a flexible organization at the chromosomal level, showing inter- and intra-specific variation in number and position in karyotypes. Sequences from the genome of Imparfinis schubarti (Heptapteridae) were isolated, aiming to understand the organization of 5S rDNA families in the fish genome. The isolation of 5S rDNA from the genome of I. schubarti was carried out by reassociation kinetics (C0t) and PCR amplification. The obtained sequences were cloned for the construction of a micro-library. The obtained clones were sequenced and hybridized in I. schubarti and Microglanis cottoides (Pseudopimelodidae) for chromosome mapping. An analysis of the sequence alignments with other fish groups was accomplished. Both methods were effective when using 5S rDNA for hybridization in I. schubarti genome. However, the C0t method enabled the use of a complete 5S rRNA gene, which was also successful in the hybridization of M. cottoides. Nevertheless, this gene was obtained only partially by PCR. The hybridization results and sequence analyses showed that intact 5S regions are more appropriate for the probe operation, due to conserved structure and motifs. This study contributes to a better understanding of the organization of multigene families in catfish's genomes.

Evidence of birth-and-death evolution of 5S rRNA gene in Channa species (Teleostei, Perciformes)

In higher eukaryotes, minor rDNA family codes for 5S rRNA that is arranged in tandem arrays and comprises of a highly conserved 120 bp long coding sequence with a variable non-transcribed spacer (NTS). Initially the 5S rDNA repeats are considered to be evolved by the process of concerted evolution. But some recent reports, including teleost fishes suggested that evolution of 5S rDNA repeat does not fit into the concerted evolution model and evolution of 5S rDNA family may be explained by a birth-and-death evolution model. In order to study the mode of evolution of 5S rDNA repeats in Perciformes fish species, nucleotide sequence and molecular organization of five species of genus Channa were analyzed in the present study. Molecular analyses revealed several variants of 5S rDNA repeats (four types of NTS) and networks created by a neighbor net algorithm for each type of sequences (I, II, III and IV) did not show a clear clustering in species specific manner. The stable secondary structure is predicted and upstream and downstream conserved regulatory elements were characterized. Sequence analyses also shown the presence of two putative pseudogenes in Channa marulius. Present study supported that 5S rDNA repeats in genus Channa were evolved under the process of birth-and-death.

Unraveling the Sex Chromosome Heteromorphism of the Paradoxical Frog Pseudis tocantins

The paradoxical frog Pseudis tocantins is the only species in the Hylidae family with known heteromorphic Z and W sex chromosomes. The Z chromosome is metacentric and presents an interstitial nucleolar organizer region (NOR) on the long arm that is adjacent to a pericentromeric heterochromatic band. In contrast, the submetacentric W chromosome carries a pericentromeric NOR on the long arm, which is adjacent to a clearly evident heterochromatic band that is larger than the band found on the Z chromosome and justify the size difference observed between these chromosomes. Here, we provide evidence that the non-centromeric heterochromatic bands in Zq and Wq differ not only in size and location but also in composition, based on comparative genomic hybridization (CGH) and an analysis of the anuran PcP190 satellite DNA. The finding of PcP190 sequences in P. tocantins extends the presence of this satellite DNA, which was previously detected among Leptodactylidae and Hylodidae, suggesting that this family of repetitive DNA is even older than it was formerly considered. Seven groups of PcP190 sequences were recognized in the genome of P. tocantins. PcP190 probes mapped to the heterochromatic band in Wq, and a Southern blot analysis indicated the accumulation of PcP190 in the female genome of P. tocantins, which suggests the involvement of this satellite DNA in the evolution of the sex chromosomes of this species.

The 5S rDNA in two Abracris grasshoppers (Ommatolampidinae: Acrididae): molecular and chromosomal organization

The 5S ribosomal DNA (rDNA) sequences are subject of dynamic evolution at chromosomal and molecular levels, evolving through concerted and/or birth-and-death fashion. Among grasshoppers, the chromosomal location for this sequence was established for some species, but little molecular information was obtained to infer evolutionary patterns. Here, we integrated data from chromosomal and nucleotide sequence analysis for 5S rDNA in two Abracris species aiming to identify evolutionary dynamics. For both species, two arrays were identified, a larger sequence (named type-I) that consisted of the entire 5S rDNA gene plus NTS (non-transcribed spacer) and a smaller (named type-II) with truncated 5S rDNA gene plus short NTS that was considered a pseudogene. For type-I sequences, the gene corresponding region contained the internal control region and poly-T motif and the NTS presented partial transposable elements. Between the species, nucleotide differences for type-I were noticed, while type-II was identical, suggesting pseudogenization in a common ancestor. At chromosomal point to view, the type-II was placed in one bivalent, while type-I occurred in multiple copies in distinct chromosomes. In Abracris, the evolution of 5S rDNA was apparently influenced by the chromosomal distribution of clusters (single or multiple location), resulting in a mixed mechanism integrating concerted and birth-and-death evolution depending on the unit.

First insights on the retroelement Rex1 in the cytogenetics of frogs

Background

While some transposable elements (TEs) have been found in the sequenced genomes of frog species, detailed studies of these elements have been lacking. In this work, we investigated the occurrence of the Rex1 element, which is widespread in fish, in anurans of the genus Physalaemus. We isolated and characterized the reverse transcriptase (RT)-coding sequences of Rex1 elements of five species of this genus.

Results

The amino acid sequences deduced from the nucleotide sequences of the isolated fragments allowed us to unambiguously identify regions corresponding to domains 3–7 of RT. Some of the nucleotide sequences isolated from Physlaemus ephippifer and P. albonotatus had internal deletions, suggesting that these fragments are likely not active TEs, despite being derived from a Rex1 element. When hybridized with metaphase chromosomes, Rex1 probes were revealed at the pericentromeric heterochromatic region of the short arm of chromosome 3 of the P. ephippifer karyotype. Neither other heterochromatin sites of the P. ephippifer karyotype nor any chromosomal regions of the karyotypes of P. albonotatus, P. spiniger and P. albifrons were detected with these probes.

Conclusions

Rex1 elements were found in the genomes of five species of Physalaemus but clustered in only the P. ephippifer karyotype, in contrast to observations in some species of fish, where large chromosomal sites with Rex1 elements are typically present.

Electronic supplementary material

The online version of this article (doi:10.1186/s13039-015-0189-5) contains supplementary material, which is available to authorized users.