Evolutionary Dynamics of 5S rDNA and Recurrent Association of Transposable Elements in Electric Fish of the Family Gymnotidae (Gymnotiformes): The Case of Gymnotus mamiraua

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

Karyotype description and comparative chromosomal mapping of rDNA and U2 snDNA sequences in Eigenmannialimbata and E.microstoma (Teleostei, Gymnotiformes, Sternopygidae)

The genus Eigenmannia Jordan et Evermann,1896 includes electric fishes endemic to the Neotropical region with extensive karyotype variability and occurrence of different sex chromosome systems, however, cytogenetic studies within this group are restricted to few species. Here, we describe the karyotypes of Eigenmannialimbata (Schreiner et Miranda Ribeiro, 1903) and E.microstoma (Reinhardt, 1852) and the chromosomal locations of 5S and 18S rDNAs (ribosomal RNA genes) and U2 snDNA (small nuclear RNA gene). Among them, 18S rDNA sites were situated in only one chromosomal pair in both species, and co-localized with 5S rDNA in E.microstoma. On the other hand, 5S rDNA and U2 snRNA sites were observed on several chromosomes, with variation in the number of sites between species under study. These two repetitive DNAs were observed co-localized in one chromosomal pair in E.limbata and in four pairs in E.microstoma. Our study shows a new case of association of these two types of repetitive DNA in the genome of Gymnotiformes.

Chromosome Painting in Gymnotus carapo "Catalão" (Gymnotiformes, Teleostei): Dynamics of Chromosomal Rearrangements in Cryptic Species

The genus Gymnotus is a large monophyletic group of freshwater weakly-electric fishes, with wide distribution in Central and South America. It has 46 valid species divided into six subgenera (Gymnotus, Tijax, Tigre, Lamontianus, Tigrinus and Pantherus) with large chromosome plasticity and diploid numbers (2n) ranging from 34 to 54. Within this rich diversity, there is controversy about whether Gymnotus (Gymnotus) carapo species is a single widespread species or a complex of cryptic species. Cytogenetic studies show different diploid numbers for G. carapo species, ranging from 40 to 54 chromosomes with varied karyotypes found even between populations sharing the same 2n. Whole chromosome painting has been used in studies on fish species and recently has been used for tracking the chromosomal evolution of Gymnotus and assisting in its cytotaxonomy. Comparative genomic mapping using chromosome painting has shown more complex rearrangements in Gymnotus carapo than shown in previous studies by classical cytogenetics. These studies demonstrate that multiple chromosome pairs are involved in its chromosomal reorganization, suggesting the presence of a complex of cryptic species due to a post zygotic barrier. In the present study, metaphase chromosomes of G. carapo occidentalis "catalão" (GCC, 2n = 40, 30m/sm+10st/a) from the Catalão Lake, Amazonas, Brazil, were hybridized with whole chromosome probes derived from the chromosomes of G. carapo (GCA, 2n = 42, 30m/sm+12st/a). The results reveal chromosome rearrangements and a high number of repetitive DNA sites. Of the 12 pairs of G. carapo chromosomes that could be individually identified (GCA 1-3, 6, 7, 9, 14, 16 and 18-21), 8 pairs (GCA 1, 2, 6, 7, 9, 14, 20, 21) had homeology conserved in GCC. Of the GCA pairs that are grouped (GCA [4, 8], [5, 17], [10, 11] and [12, 13, 15]), most kept the number of signals in GCC (GCA [5, 17], [10, 11] and [12, 13, 15]). The remaining chromosomes are rearranged in the GCC karyotype. Analysis of both populations of the G. carapo cytotypes shows extensive karyotype reorganization. Along with previous studies, this suggests that the different cytotypes analyzed here may represent different species and supports the hypothesis that G. carapo is not a single widespread species, but a group of cryptic species.

Archolaemus janeae (Gymnotiformes, Teleostei): First insights into karyotype and repetitive DNA distribution in two populations of the Amazon

Archolaemus, one of the five genera of Neotropical freshwater fish of the family Sternopygidae (Gymnotiformes), was long considered a monotypic genus represented by Archolaemus blax. Currently, it consists of six species, most of them occurring in the Amazon region. There are no cytogenetic data for species of this genus. In the present study, we used classical cytogenetics (conventional staining and C-banding) and molecular cytogenetics (probes of telomeric sequences and multigenic families 18S rDNA, 5S rDNA, and U2 snDNA) to study the karyotype of Archolaemus janeae from Xingu and Tapajós rivers in the state of Pará (Brazil). The results showed that the two populations have identical karyotypes with 46 chromosomes: four submetacentric and 42 acrocentric (2n = 46; 4m/sm + 42a). Constitutive heterochromatin occurs in the centromeric region of all chromosomes, in addition to small bands in the interstitial and distal regions of some pairs. The 18S rDNA occurs in the distal region of the short arm of pair 2; the 5S rDNA occurs in five chromosome pairs; and the U2 snDNA sequence occurs in chromosome pairs 3, 6, and 13. No interstitial telomeric sequence was observed. These results show karyotypic similarity between the studied populations suggesting the existence of a single species and are of great importance as a reference for future cytotaxonomic studies of the genus.

Chromosomal Diversity in Two Allopatric Populations of Farlowella hahni Meinken 1937 (Teleostei: Siluriformes): Cytogenetics and Cytochrome b Analyses

Farlowella is the second richest genus in Loricariinae, broadly distributed in freshwater streams and rivers of South America. In this article, we aimed to expand on the cytogenetic and molecular data available for two allopatric populations of Farlowella hahni. Both populations had diploid chromosome number 58, but with karyotype differences, indicative of chromosomal rearrangements. C-banding showed large heterochromatic blocks at telomeric regions in acrocentric chromosomes in both populations. Fluorescence in situ hybridization (FISH) revealed a single 18S rDNA site in both populations and a single 5S rDNA site for individuals from lower Paraná River basin (native region) and multiple 5S rDNA sites for individuals from upper Paraná River basin (non-native region). Mitochondrial sequence analyses did not separate the two F. hahni populations. The cytogenetic and molecular data obtained are relevant in a preliminary study and suggested the existence of cryptic diversity and the hypothesis that at least two Farlowella lineages may coexist in the Paraná basin.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Chromosomal mapping of repetitive sequences in Hyphessobrycon eques (Characiformes, Characidae): a special case of the spreading of 5S rDNA clusters in a genome

Cytogenetic data showed a variation in diploid chromosome number in the genus Hyphessobrycon ranging from 2n = 46 to 52, and studies involving repetitive DNA sequences are scarce in representatives of this genus. The purpose of this paper was the chromosomal mapping of repetitive sequences (rDNA, histone genes, U snDNA and microsatellites) and investigation of the amplification of 5S rDNA clusters in the Hyphessobrycon eques genome. Two H. eques populations displayed 2n = 52 chromosomes, with the acrocentric pair No. 24 bearing Ag-NORs corresponding with CMA₃+/DAPI-. FISH with a 18S rDNA probe identified the NORs on the short (p) arms of the acrocentric pairs Nos. 22 and 24. The 5S rDNA probe visualized signals on almost all chromosomes in genomes of individuals from both populations (40 signals); FISH with H3 histone probe identified two chromosome pairs, with the pericentromeric location of signals; FISH with a U2 snDNA probe identified one chromosome pair bearing signals, on the interstitial chromosomal region. The mononucleotide (A), dinucleotide (CA) and tetranucleotide (GATA) repeats were observed on the centromeric/pericentromeric and/or terminal positions of all chromosomes, while the trinucleotide (CAG) repeat showed signals on few chromosomes. Molecular analysis of 5S rDNA and non-transcribed spacers (NTS) showed microsatellites (GATA and A repeats) and a fragment of retrotransposon (SINE3/5S-Sauria) inside the sequences. This study expanded the available cytogenetic data for H. eques and demonstrated to the dispersion of the 5S rDNA sequences on almost all chromosomes.

Effects of environmental pollution on the rDNAomics of Amazonian fish

Pollution is a growing environmental problem throughout the world, and the impact of human activities on biodiversity and the genetic variability of natural populations is increasingly preoccupying, given that adaptive processes depend on this variability, in particular that found in the repetitive DNA. In the present study, the mitochondrial DNA (COI) and the distribution of repetitive DNA sequences (18S and 5S rDNA) in the fish genome were analysed in fish populations inhabiting both polluted and unpolluted waters in the northern Amazon basin. The results indicate highly complex ribosomal sequences in the fish genome from the polluted environment because these sequences are involved primarily in the maintenance of genome integrity, mediated by a systematic increase in the number of copies of the ribosomal DNA in response to changes in environmental conditions.

Cytogenetic markers as tools in delimiting species of the highly diverse Neotropical fish Bryconamericus (Characiformes: Characidae)

ABSTRACT Bryconamericus is a highly diverse group of characid fishes, being cytogenetic a valuable tool for the delimitation of species. Bryconamericus aff. iheringii (Upper Uruguay/Lower Paraná), B. coeruleus (Upper Paraná), B. cf. ecai e B. cf. eigenmanni (Upper Uruguay) were studied cytogenetically, and presented 2n=52 chromosomes, with interpopulational/interspecific variation of karyotype and fundamental number. Heterochromatin was evidenced in pericentromeric, telomeric and interstitial regions, and it was shown to be an important cytogenetic marker. Single nucleolar organizing regions (NORs) were found in B. cf. eigenmanni, B. cf. ecai and B. aff. iheringii (Lower Paraná), and multiple in B. aff. iheringii (Upper Uruguay) and B. coeruleus, with occurrence of two patterns for the first species, and three for the second. The 5S/18S rDNA-FISH confirmed the location of the NORs and showed single 5S rDNA cistrons only in B. aff. iheringii (Lower Paraná), evidencing the dispersion of both genes, often co-located, in the karyotype of the others species. The data of this work contribute for the delimitation of the species of the genus. Co-localization of ribosomal genes may represent a plesiomorphic condition for the group, and their dispersion suggest the occurrence of duplication, pseudogeneization and transposition events mediated by mobile genetic elements.

Particular Chromosomal Distribution of Microsatellites in Five Species of the Genus Gymnotus (Teleostei, Gymnotiformes)

Microsatellites show great abundance in eukaryotic genomes, although distinct chromosomal distribution patterns might be observed, from small dispersed signals to strong clustered motifs. In Neotropical fishes, the chromosome mapping of distinct microsatellites was employed several times to uncover the origin and evolution of sex and supernumerary chromosomes, whereas a detailed comparative analysis considering different motifs at the chromosomal level is scarce. Here, we report the chromosomal location of several simple sequence repeats (SSRs) in distinct electric knife fishes showing variable diploid chromosome numbers to unveil the structural organization of several microsatellite motifs in distinct Gymnotus species. Our results showed that some SSRs are scattered throughout the genomes, whereas others are particularly clustered displaying intense genomic compartmentalization. Interestingly, the motifs CA, GA, and GAG exhibited a band-like pattern of hybridization, useful for the identification of homologous chromosomes. Finally, the colocalization of SSRs with multigene families is probably related to the association of microsatellites with gene spacers in this case.

Dynamics of tandemly repeated DNA sequences during evolution of diploid and tetraploid botiid loaches (Teleostei: Cobitoidea: Botiidae)

Polyploidization has played an important role in the evolution of vertebrates, particularly at the base of Teleostei-an enormously successful ray-finned fish group with additional genome doublings on lower taxonomic levels. The investigation of post-polyploid genome dynamics might provide important clues about the evolution and ecology of respective species and can help to decipher the role of polyploidy per se on speciation. Few studies have attempted to investigate the dynamics of repetitive DNA sequences in the post-polyploid genome using molecular cytogenetic tools in fishes, though recent efforts demonstrated their usefulness. The demonstrably monophyletic freshwater loach family Botiidae, branching to evolutionary diploid and tetraploid lineages separated >25 Mya, offers a suited model group for comparing the long-term repetitive DNA evolution. For this, we integrated phylogenetic analyses with cytogenetical survey involving Giemsa- and Chromomycin A3 (CMA3)/DAPI stainings and fluorescence in situ hybridization with 5S/45S rDNA, U2 snDNA and telomeric probes in representative sample of 12 botiid species. The karyotypes of all diploids were composed of 2n = 50 chromosomes, while majority of tetraploids had 2n = 4x = 100, with only subtle interspecific karyotype differences. The exceptional karyotype of Botia dario (2n = 4x = 96) suggested centric fusions behind the 2n reduction. Variable patterns of FISH signals revealed cases of intraspecific polymorphisms, rDNA amplification, variable degree of correspondence with CMA3+ sites and almost no phylogenetic signal. In tetraploids, either additivity or loci gain/loss was recorded. Despite absence of classical interstitial telomeric sites, large blocks of interspersed rDNA/telomeric regions were found in diploids only. We uncovered different molecular drives of studied repetitive DNA classes within botiid genomes as well as the advanced stage of the re-diploidization process in tetraploids. Our results may contribute to link genomic approach with molecular cytogenetic analyses in addressing the origin and mechanism of this polyploidization event.

Extensive Karyotype Reorganization in the Fish Gymnotus arapaima (Gymnotiformes, Gymnotidae) Highlighted by Zoo-FISH Analysis

The genus Gymnotus (Gymnotiformes) contains over 40 species of freshwater electric fishes exhibiting a wide distribution throughout Central and South America, and being particularly prevalent in the Amazon basin. Cytogenetics has been an important tool in the cytotaxonomy and elucidation of evolutionary processes in this genus, including the unraveling the variety of diploid chromosome number (2n = from 34 to 54), the high karyotype diversity among species with a shared diploid number, different sex chromosome systems, and variation in the distribution of several Repetitive DNAs and colocation and association between those sequences. Recently whole chromosome painting (WCP) has been used for tracking the chromosomal evolution of the genus, showing highly reorganized karyotypes and the conserved synteny of the NOR bearing par within the clade G. carapo. In this study, painting probes derived from the chromosomes of G. carapo (GCA, 2n = 42, 30 m/sm + 12 st/a) were hybridized to the mitotic metaphases of G. arapaima (GAR, 2n = 44, 24 m/sm + 20 st/a). Our results uncovered chromosomal rearrangements and a high number of repetitive DNA regions. From the 12 chromosome pairs of G. carapo that can be individually differentiated (GCA1-3, 6, 7, 9, 14, 16, and 18-21), six pairs (GCA 1, 9, 14, 18, 20, 21) show conserved homology with GAR, five pairs (GCA 1, 9, 14, 20, 21) are also shared with cryptic species G. carapo 2n = 40 (34 m/sm + 6 st/a) and only the NOR bearing pair (GCA 20) is shared with G. capanema (GCP 2n = 34, 20 m/sm + 14 st/a). The remaining chromosomes are reorganized in the karyotype of GAR. Despite the close phylogenetic relationships of these species, our chromosome painting studies demonstrate an extensive reorganization of their karyotypes.

Cytogenetics of two Farlowella species (Loricariidae: Loricariinae): implications on the taxonomic status of the species

ABSTRACT Farlowella is one of the most diverse genera of the Loricariinae, restricted to South America rivers. The taxonomic and phylogenetic relationships among its species are contentious and, while genetic studies would contribute to the understanding of their relationships, the only available datum refer to the karyotype description of only one species. In the present study two Amazonian species, Farlowella cf. amazonum and F. schreitmuelleri, were analyzed using conventional and molecular cytogenetic procedures. Both species had diploid chromosome number 58, but different fundamental numbers (NF) 116 and 112, respectively, indicative of chromosomal rearrangements. C-banding is almost poor, especially in F. cf. amazonum, and occurs predominantly in the centromeric and in some telomeric regions, although genome of F. schreitmuelleri possessed a much larger heterochromatin amount then those of F. cf. amazonum. The chromosomes bearing the NOR sites were likely the same for both species, corresponding to the 1st metacentric pair in F. cf. amazonum and to the 28th acrocentric in F. schreitmuelleri. The location of the 5S rDNA was species-specific marker. This study expanded the available cytogenetic data for Farlowella species and pointed the remarkable karyotype diversity among species/populations, indicating a possible species complex within genus.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.)

Background

Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes.

Results

The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30–38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10–30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides.

Conclusions

Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3774-7) contains supplementary material, which is available to authorized users.