Open access peer-reviewed chapter

Composition and Nature of Heterochromatin in the Electrical Fish (Knifefishes) Gymnotus (Gymnotiformes: Gymnotidae)

Written By

Maelin da Silva, Daniele Aparecida Matoso, Vladimir Pavan Margarido, Eliana Feldberg and Roberto Ferreira Artoni

Submitted: 28 October 2020 Reviewed: 12 April 2021 Published: 07 May 2021

DOI: 10.5772/intechopen.97673

From the Edited Volume

Cytogenetics - Classical and Molecular Strategies for Analysing Heredity Material

Edited by Marcelo L. Larramendy and Sonia Soloneski

Chapter metrics overview

415 Chapter Downloads

View Full Metrics

Abstract

Fishes of the genus Gymnotus have been suggested as a good model for biogeographic studies in the South American continent. In relation to heterochromatin, species of this genus have blocks preferably distributed in the centromeric region. The content of these regions has been shown to be variable, with description of transposable elements, pseudogenes of 5S rDNA and satellite sequences. In G. carapo Clade, although geographically separated, species with 2n = 54 chromosomes share the distribution of many 5S rDNA sites, a unique case within the genus. Here, repetitive DNA sequences from G. sylvius (2n = 40) and G. paraguensis (2n = 54) were isolated and mapped to understand their constitution. The chromosome mapping by FISH showed an exclusive association in the centromeres of all chromosomes. However, the cross-FISH did not show positive signs of interspecific hybridization, indicating high levels of heterochromatic sequence specificity. In addition, COI-1 sequences were analyzed in some species of Gymnotus, which revealed a close relationship between species of clade 2n = 54, which have multiple 5S rDNA sites. Possibly, the insertion of retroelements or pseudogenization and dispersion of this sequence occurred before the geographic dispersion of the ancestor of this clade from the Amazon region to the hydrographic systems of Paraná-Paraguay, a synapomorphy for the group.

Keywords

  • FISH
  • Biogeography
  • Satellite DNA
  • rDNA 5S
  • C0t-1

1. Introduction

Repetitive DNA sequences are broadly distributed in eukaryotes genomes [1] and are classified into two categories: 1) repetitive sequences arranged in tandem as satellite, minisatellite, or microsatellite DNAs composed of hundreds of base pairs repeated a thousand times or more in each genome; and 2) moderately to highly repetitive sequences spread throughout the genome as retroelements or transposable mobile elements [2].

Copies of repetitive sequences are commonly associated with heterochromatin regions that can be visualized by C banding. These sequences are extremely important to the functional and structural organization of the eukaryote genome, composing, for example, the pericentromeric heterochromatin regions [3, 4]. The heterochromatin in fish chromosomes is largely located in pericentromeric regions and has structural functions [5, 6, 7].

The studies about the location of repetitive sequences on chromosomes has broadened the knowledge of the structural organization of chromatin in fish, revealing an association of ribosomal DNA, telomeric sequences, transposition elements and satellite sequences in chromosomal rearrangements and weak break points [8, 9], in the fixation of sex chromosomes [10, 11], in the expansion of heterochromatin [12, 13] and in gene regulation [14]. This advances in molecular cytogenetics have demonstrated that repetitive DNA sequences are useful as chromosomal markers in studies of species evolution and can provided valuable information about sex chromosome systems and chromosomal rearrangements [15]. The mapping of Non-long terminal repeat (non-LTR) retrotransposable elements, the Rex in the fish species, for example, demonstrated strong FISH signals in heterochromatin regions [16].

Neotropical electric fish species, order Gymnotiformes, have shown their heterochromatin to be preferentially distributed in the pericentromeric regions of their chromosomes [13, 17, 18].

The investigation of repetitive sequences in Gymnotus seems promising for understanding the chromosomal evolutionary dynamics in the genus. DNA probes and chromosomal painting were applied to investigate chromosomal rearrangements in two species of the G. carapo complex, and rearrangements were found between the two species, involving several pairs of chromosomes, corroborating the existence of cryptic species in this group, in addition to the recent speciation between them [19, 20]. Analysis of the satellitome of some species by Next Generation Sequencing (NGS) revealed sets of conserved satellite sequences, but the CA, GA and GAG motifs when mapped revealed a useful band for identification of homologous chromosomes [21].

The ribosomal DNA mapping can also provide new answers about chromosomal evolution in the genus, and even serve as a tool for understanding geographic distribution patterns [22]. In the species that had the 18S rDNA mapped, the proposal of only one pair carrying the conserved nucleolus organizing regions (NORs) is plausible [17, 18, 23, 24]. Regarding the 5S rDNA, the group behaves as an attractive model for evolutionary studies, showing a species-specific pattern. The evolution of this gene family receives special attention in the species that comprise the group G. carapo with a diploid number of 54 chromosomes, which have from 14 to 19 pairs identified with this sequence [21]. This situation is totally different from the pattern observed for other species within the group and shared only by species carrying 2n = 54 chromosomes.

Advertisement

2. Biogeography of electric fish and repetitive DNA sequences

The complex history of the formation of the South American rivers is fundamental to explain the diversity and distribution of aquatic biota in this region. Successive continental geomorphological changes, such as the one that resulted in the formation of the Andes, associated with historic and biological factors allow the identification of patterns that led to the formation of the largest and most diverse freshwater ichthyofauna on the planet [25, 26]. Such changes alter the drainage scenario forming lakes, capture headwaters and basins of varying sizes, or even isolate populations for certain periods, favoring the diversification of biota by vicariance and allopatric, in addition to promoting subsequent drainage coalescence, leading to enrichment and contact of organisms [26, 27].

The Gymnotiformes order comprises electric fish or knifefishes. Endemic to the Neotropical region [28], which are widely distributed, from the Pampas in Argentina to Chiapas, Mexico and reach their greatest diversity and abundance in the Amazon basin [29, 30, 31]. Members of the Gymnotiformes order, are unique in their ability to produce and recognize electrical signals, never left the South America plate, since their electrosensory system is not capable of functioning in brackish or salt water [21].

Regarding the karyotype, Gymnotus is also a diverse group, with a diploid number ranging from 34 chromosomes, verified in G. capanema [33] to 54, found in several species of the G. carapo group [18, 21, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34]. Sexual chromosomal heteromorphism has already been recorded in some, and all records follow a turnover model, with little heterochromatin and of recent origin [24, 35].

Among the Gymnotus, an intriguing fact has been reported for species that have 2n = 54 chromosomes, present in the carapo Clade: G. cf. carapo, G. inequilabiatus, G. paraguensis and G. mamiraua. These species share the diploid number and also a large number of 5S rDNA sites [13], but they are geographically separated in different hydrographic basins. Though, G. cf. carapo, G. inequilabiatus, G. paraguensis occur in the Paraná-Paraguay basin, and G. mamiraua is distributed in the Amazon basin (Figure 1). The number of 5S DNAr sites in G. paraguensis, G. inaequilabiatus and G. cf. carapo are 19, 17 and 15 bearing pairs, respectively [36, 37, 38], and 14 pairs in G. mamiraua [22].

Figure 1.

Location of species with 2n = 54 chromosomes in the Amazon and Paraná-Paraguay basins: a) G. mamiraua in the Amazon basin; b) G. carapo; c) G. paraguensis and d) G. inaequilabiatus. The pairs with 5S rDNA are represented schematically along with each species.

For G. paraguensis and G. mamiraua, the 5S rDNA are present in the pericentromeric region and when sequenced, a transposition element similar to Tc1 was identified in the NTS (non-transcribed spacer), suggested as one of the mechanisms used for the transposition of this rDNA in chromosomes [22, 36]. In G. inequilabiatus, which has 17 sites of 5S rDNA, the sequencing of some clones with classes I and II revealed sequences of the TATA type in NTS, which are associated with the regulation of the expression of the 5S rDNA gene [37].

As a rule, studies investigating this ribosomal gene associate an increase in the number of 5S rDNA sites with the prese-nce of pseudogenes, which would originate through duplication of copies, resulting from the transposition or duplication of genomic DNA, which would be facilitated by the organization tandem of this rDNA family. In addition, locating these sequences in the terminal portion of the chromosomes would facilitate the process of pseudogenization and the association with transposable elements [39].

Although the mechanism of origin and dispersion of these sequences in the 2n = 54 Gymnotus genome is still unclear, it is possible to assume that the increase in the number of 5S-like rDNA sites and fixation in the genome accompanies the migratory dispersion of these species. In the case of Gymnotus (2n = 54), infection by this transposable element may have probably occurred in a common ancestor for species with 2n = 54 chromosomes, before the final isolation between the Amazon and Paraná-Paraguay basins, since these species are currently found in different river basins. Gymnotus mamiraua was restricted to the Amazon basin show the smallest number of sites, 14 pairs, maybe the original satellite sequence like-5S rDNA; however, the isolated species in the Paraná-Paraguay basin experienced divergence from this satellite sequence, differing in number pairs with sequence, such as: G. inaequilabiatus (17 pairs), G. paraguensis (19 pairs), and G. carapo (15 pairs), (Figure 2) [13, 22].

Figure 2.

Topology obtained by the method of clustering neighbors (neighbor joining, NJ) for species of Gymnotus of the carapo clade, using the distance model Kimura-2-parameters (K2P) for the mitochondrial gene COI. The values present in the nodes are the bootstrap values (> 50%) calculated from 1000 replicates and the diploid number superimposed on the tree. The distribution of species is indicated on the map. * absence of coordinates, group not plotted on the map.

Evidence of the connection between South American watersheds is reflected in the evolutionary history of the fish. Other species corroborate the hypothesis of interconnection between these two systems. Migratory species of the genus Prochilodus are considered as indicators of this panorama, where the relationship of sister clades of that genus has been useful for understanding the separation of the basins [Magdalena (Orinoco (Amazonas-Paraná)]. The use of molecular clock methodologies, with species in the group, estimates that the isolation between the Amazonas-Paraná basins must have occurred between 2.3–4.1 million years, with an estimated coalescence between strains of 1–3-3 million years [40]. Morphological and molecular analyzes with rheophilic taxa increase the evidence of connection between these basins [41]. The final rupture of the connection between the Paraná-Paraguay and Amazonas-Orinoco basins is inferred as a recent event, with the final uplift of the Andes that changed the course of these basins and isolated them [25, 42], with final separation estimates still uncertain and very peculiar to each group of fish investigated [43].

In order to understand the relationships between Gymnotus, here, an analysis of a 556 bp (base pairs) COI mitochondrial gene fragment was conducted in 21 specimens from the carapo Clade from the central Amazon basin and four more specimens from a population in the Munim River, MA. Sequences deposited in the NCBI database were added to this analysis, with the access numbers described in Table 1. The analysis of 51 sequences resulted in eight clades (Figure 2). The data revealed a high percentage for the distance between the species identified as G. carapo (Table 1, Figure 2). The sequences of individuals identified as G. carapo do not form exclusive groups. Three clades with high rates of genetic distance were observed between G. carapo “Catalão” (2n = 40, XX/XY) and G. carapo Maranhão (2n = 42), but between G. carapo “Catalão” and G. carapo do the distance between Amazonas and Peru was 0.02% (Table 2). The species G. mamiraua forms an exclusive clade with an average intraspecific genetic distance of 0.12% and G. inaequilabiatus the average was 0.2% (Table 3). These two species showed an interspecific distance of only 0.6%.

SpeciesVoucherLocalityNCBI accession number
G. sylviusLBP7069Alto ParanáGU.064995.1
G. sylviusLBP31958Alto ParanáGU.701779.1
G. sylviusLBP7070Alto ParanáGU.702209.1
G. sylviusLBP8831Alto ParanáGU.7017821
G. sylviusLBP27380Alto ParanáGU.701778.1
G. sylviusLBP29096Alto ParanáGU.702207.1
G. sylviusLBP25852Alto ParanáGU.701758.1
G. sylviusLBP31933Alto ParanáGU.701767.1
G. sylviusLBP25853Alto ParanáGU.701762.1
G. sylviusLBPV27382Alto ParanáJN.988881.1
G. sylviusLBPV27381Alto ParanáJN.988880.1
G. sylviusLBP9823Alto ParanáGU.701780.1
G. pantanalMZUEL5644Alto ParanáKF.359492.1
G. pantanalLBP31929Alto ParanáGU.701776.1
G. pantanalLBP34742Alto ParanáGU.701775.1
G. pantanalLBP34743Alto ParanáGU.701774.1
G. pantanalLBP31927Alto ParanáGU.701773.1
G. pantanalLBP31932Alto ParanáGU.701763.1
G. carapo MaraGC7MBacia Costeira- MaranhãoXXXXXXXXX
G. carapo MaraGC9MBacia Costeira- MaranhãoXXXXXXXXX
G. carapo MaraGC11MBacia Costeira- MaranhãoXXXXXXXXX
G. carapo MaraGC18MBacia Costeira- MaranhãoXXXXXXXXX
G. carapo CatGCAT11029Amazônia CentralXXXXXXXXX
G. carapo CatGCAT 11028Amazônia CentralXXXXXXXXX
G. carapo CatGCAT 11066Amazônia CentralXXXXXXXXX
G. carapo CatGCAT 11477Amazonia CentralXXXXXXXXX
G. carapoGC2006Peru AmazonasKF533344
G. carapoGC2007Peru, AmazonasKF533345.1
G. mamirauaGM11078Amazônia CentralXXXXXXXXX
G. mamirauaGM11733Amazônia CentralXXXXXXXXX
G. mamirauaGM11731Amazônia CentralXXXXXXXXX
G. mamirauaGM10947Amazônia CentralXXXXXXXXX
G. mamirauaGM10948Amazônia CentralXXXXXXXXX
G. mamirauaGM10949Amazônia CentralXXXXXXXXX
G. mamirauaGM11730Amazônia CentralXXXXXXXXX
G. inaequilabiatusMZUEL5649Alto ParanáKF.359490.1
G. inaequilabiatusLBP26331Alto ParanáGU.701766.1
G. inaequilabiatusLBP31931Alto ParanáGU.701764.1
G. inaequilabiatusLBP25850Alto ParanáGU.701760.1
G. inaequilabiatusLBP7071Alto ParanáGU.702210.1
G. inaequilabiatusLBP29097Alto ParanáGU.702208.1
G. inaequilabiatusLBP34744Alto ParanáGU.701781.1
G. ucamaraGU11575Amazônia CentralXXXXXXXXX
G. ucamaraGU11802Amazônia CentralXXXXXXXXX
G. ucamaraGU11701Amazônia CentralXXXXXXXXX
G. ucamaraGU11698Amazônia CentralXXXXXXXXX
G. ucamaraGU11574Amazônia CentralXXXXXXXXX

Table 1.

Access numbers in the NCBI of the Gymnotus species of the carapo clade used for genetic distance verification through DNA barcode.

Intraspecific genetic distance
Gp 1 - G. sytlvius0.013526931
Gp 2 - G. Pantanal0
Gp 3 - G. carapo “Maranhão”0.000893724
Gp 4 - G. carapo “Catalão”0.000928937
Gp 5 - G carapo0.001858738
Gp 6 - G.mamiraua0.001241471
Gp 7 - G. inaequilabiatus0.002136221
Gp 8 - G. ucamara0.003103677

Table 2.

Intraspecific genetic distance based on mutations of the COI gene, using the K-2-P model.

Interspecific diversity
Gp- 1Gp-2Gp- 3Gp- 4Gp- 5Gp- 6Gp- 7
Gp 1 G. sylvius
Gp 2 G. Pantanal0.062
Gp 3 G. carapo “Maranhão”0.1680.156
Gp 4 G. carapo “Catalão”0.0300.0690.174
Gp 5 G. carapo0.0280.0710.1710.002
Gp 6 G. mamiraua0.0480.0770.1810.0550.053
Gp 7 G. inaequilabiatus0.0460.0740.1780.0530.0510.006
Gp 8 G. ucamara0.1790.1590.0510.1830.1810.1780.174

Table 3.

Matrix of the means of genetic distance (K2P) for the COI gene obtained among the different species of Gymnotus from carapo clade.

The results with the mitochondrial DNA COI obtained here validate the species of carapo Clade, with high rates of genetic distance between groups of species identified as G. carapo, confirming the probable existence of more than one taxonomic unit in this group. However, this did not occur for two species: G. mamiraua from the Amazon basin and G. inaequilabiatus from Alto Paraná, for which a low interspecific distance was detected and with a distance of less than 1% was detected, possibly causing recurrent speciation.

Cytogenetic data already pointed to the proximity between these two species, because in addition to sharing the same diploid number, both had many sites of 5S rDNA, a condition that is not common in fish. We suggest, by the results of diploid number and dispersion of 5S rDNAr, that the species G. paraguensis and G. carapo (2n = 54) from Alto Paraná are also related to the last two already mentioned (Figure 2). Thus, what has been verified is that the 5S rDNA has been a potential tool in helping to reconstruct the steps involved in the evolution and biogeographic history of the species of the genus Gymnotus, especially the carapo Clade, inferred both by the chromosomal mapping and by the molecular analysis of that gene.

2.1 New repetitive sequences studies in the genus Gymnotus

The eukaryote genome is characterized by presenting nucleotide sequences with varied arrangements, generally forming two large groups, gene regions and repetitive DNA sequences. In fish of the genus Gymnotus, the latter has been associated with heterochromatin. The location of the heterochromatin is reported to be preferentially organized in the centromeric and pericentromeric regions [18, 19, 36, 44].

The prospection of repetitive sequences by the technique of DNA reassociation kinetics (C0t-1) proves to be a safe and fast technique for obtaining copies of highly and moderately repetitive DNA sequences [45]. Thus, it is possible to build libraries and screening repetitive DNAs, and has been used to isolate the highly repeated fraction of the plant genome [46] and animals [47] to significantly expand our knowledge of the organization of their chromosomes.

In the present study, repetitive DNA sequences were isolated by C0t-1 (Figure 3) and mapped in two species of electric fish, G. paraguensis and G. sylvius. Our objective was identifying the heterochromatin compositions and to verify the presence of repetitive sequences originated from transposable elements. Twelve specimens (four females and eight males) of G. paraguensis and 21 specimens (seven females and 14 males) of G. sylvius collected at Piquiri River, Paraná – Brazil were analyzed.

Figure 3.

1% agarose gel electrophoresis of DNA treated for kinetic re-association method (C0t – 1). a) after one minute in autoclave the DNA appears as a trail; b) after treatment with S1 nuclease enzyme the DNA appears with defined length between 100 and 400 bp. (methodology described in complementary material).

The isolated probes from G. paraguensis had lengths of 473 and 206 bp (Figure 4a) and when submitted to BLAST (http://www.ncbi.nlm.nih.gov/blast), clone 2 was found to have 95% identity with microsatellites from Salmo salar. Dinucleotide repetitions were observed in this clone. The clones of G. sylvius had lengths of 124, 202 and 123 bp (Figure 4b). The results of hybridization with total C0t-1 in both species were coincident with heterochromatic sites, according to the description by C banding [18].

Figure 4.

a) Alignment of two isolated clones of G. paraguensis, b) alignment of three clones of G. sylvius.

The heterochromatin of Gymnotus has been reported to be preferentially organized in the centromeric regions [18, 19], as detected in the present study (Figure 5a and b), and inmost species of the Gymnotidae family [36, 44]. Furthermore, heteromorphisms in length between homologous chromosomes of the NOR regions were observed in G. paraguensis, indicating differential accumulations of heterochromatin regions (Figure 5b). However, when the probes of one species were hybridized with the other (cross-FISH) no positive hybridization signals were observed (data not shown).

Figure 5.

Fluorescent in situ hybridization (FISH) performed with probes isolated for kinetic re-association of DNA (C0t–1). a) G. sylvius karyotype; b) G. paraguensis karyotype. Observe in b the heteromorphism of nucleolus organizer regions in pair 25. Bar = 10 μm.

According to Charlesworth et al. [48] and Topp and Dave [49], the regions located nearest to the centromere show fast evolutionary rates due to low recombination, initiating the accumulation of repetitive DNA sequences, which explains its specificity. This association between heterochromatin and repetitive sequences is fundamental to the organization of important chromosomal structures such as the centromere. In a study with Oreochromis niloticus using the GISH (Genomic in situ Hybridization) methodology, the heterochromatin present in pericentromeric regions was found to be species-specific [50] and composed of repetitive and transposable elements [51, 52, 53].

The repetitive elements isolated from G. sylvius and G. paraguensis in the present study were located in the pericentromeric region, coincident with the heterochromatin observed by C banding. High levels of specificity of the isolated probes and of the species’ genomes were assumed because no signals were observed in crossed FISH analyses - confirming the heterogeneous composition of the heterochromatin of these species.

Sequencing analyses showed exclusive sequences for both species, and although repetitive elements in the heterochromatin regions are present in distant eukaryotes groups such as Drosophila and plants [54, 55] as important structural regions of the genome, the structural functions of these sequences in G. sylvius and G. paraguensis are not yet known.

The Y chromosome of Eigenmannia virescens isolated by Henning et al. [56] had large amounts of heterochromatin and physical mapping with Y chromosome probe was performed with closely related species without differentiated sex chromosome systems, and the probes hybridized only in the centromeric and telomeric regions.

In addition to the two species analyzed in the present study, two other species of Gymnotiformes have had their repetitive DNA sequences analyzed and described. Claro [57] isolated the repetitive sequences of G. sylvius and G. carapo by enzymatic digestion with AluI and HaeIII. The isolated fragments with 300 bp showed dispersed distributions in both the species and similar locations using both enzymes. Furthermore, a transposable SINE element that labeled different regions in G. carapo was identified; all the sequences showed disperse labeling that was not coincident with heterochromatin regions, suggesting an important function in the evolution and organization of non-codifying DNA regions [58].

More recently, the publication of Satellitome results has been awaited, a global study by NGS and bioinformatics of all satellite sequences of Gymnotus species from the main clades. According to previous data released by the authors, the massive characterization of satellite DNA families processes of genetic differentiation and the dynamics of these sequences among the representatives stand out [59], which corroborates our findings in the present work.

Advertisement

3. Conclusions

Molecular studies with the multigene family 5S rDNA in electric fish (Gymnotidae) have advanced a lot in recent years. The chromosomal location and distribution have been particularly interesting, since all species of the genus Gymnotus, analyzed so far, show great variability. Its distribution on the chromosomes of species with 2n = 54 shows to be a biogeographic marker, suggesting that the speciation of this group was recent due to migratory expansion. On the other hand, the species analyzed in the present study showed repetitive sequences composed of microsatellite replications species-specific, present around the centromeres. Although the location of the heterochromatin in Gymnotus was conserved, they had different constitutions and could represent important evolutionary markers for cytotaxonomic studies of this group.

Advertisement

Acknowledgments

The current study was supported in part by INCT ADAPTA II funded by CNPq - Brazilian National Research Council (465540/2014-7). The authors are grateful to Miguel Airton Carvalho for his field assistance. The author M.D.S. received a scholarship from the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (grant 160155/2018-5).

Advertisement

Conflict of interest

The authors declare no conflict of interest.

Advertisement

Advertisement

A.1 Methodology for obtaining cytochrome oxidase sequences

For the Cytochrome oxidase I gene, five specimens of G. carapo “Catalão”, four specimens of G. carapo “Maranhão”, four specimens of G. ucamara and seven specimens of G. mamiraua were sequenced, using the primers BOL-COIfishF1 (5’TCAACYAATCAYAAAGATATYGGCAC´3′) and BOL-COI-fishR1 (5´-ACTTCYGGGTGRCCRAARAATCA- 3′) [60]. Genomic DNA extraction was performed using the chloroform-phenol method [61]. The PCR reactions were performed in a final volume of 25 μL, containing genomic DNA (100 ng), 10x buffer with 1.5 mM MgCl2, Taq DNA polymerase (5 U/μL), dNTPs (1 mM), pair of primers (5 pM) and Milli-Q water. The conditions for amplification were: 1 min 95° C, followed by 30 cycles of 1 min at 94° C, 1 min at 59° C, 1 min 30 sec at 72° C and final extension of 5 min at 72° C.

Sequences of Gymnotus species for the Clado carapo, available in the NCBI database, from the Paraná-Paraguay basin were added to our analysis: G. sylvius, G. inaequilabiatus, G. pantanal and G. carapo “Pantanal” (Table 1).

The calculation of intra and interspecific distance was performed using the Mega 5 software [62] with the Kimura-2-parameters evolutionary model. The identification was carried out according to the protocol established by DNA barcoding through the Neighbor Joining (NJ) analysis [63], which consists of looking for the tree with the lowest total sum of branches, using Kimura-2-parameters as an evolutionary model (K2P) [64]. The topology confidence test was performed with bootstrap analysis, containing 1000 replicates. Such analyzes were performed using the Mega 5 software [62].

The sequences were aligned in the program Clustal W [65], using the editor BioEdit 7.0 [66], were submitted to BLAST in the NCBI database (http://www.ncbi.nlm.nih.gov).

Advertisement

A.2 Methodology for obtaining cytogenetic material

Mitotic chromosomes were obtained according to the protocol described by Bertollo et al. [67]. Genomic DNA extraction was performed using the chloroform-phenol method [61]. Repetitive DNA probes were obtained using the Cot-1 DNA technique as described by Zwick et al. [45] and adapted by Vicari et al. [68]. This technique is based on re-association kinetics and enzymatic digestion by S1 nuclease. The probes were labeled by nick translation with digoxigenin 11 dUTP (Roche®) and the signals were recognized by anti-digoxigenin-rhodamin (Roche®). The hybridization techniques followed the Pinkel et al. [69] protocol, with 77% of stringency (2.5 ng/μL of each probe, 50% of deionized formamide, 10% dextran sulphate, 2X SSC at 37°C for 18 hours). The chromosomes were counterstained with DAPI in Vectashield medium (Vector®). The chromosomal preparations were analyzed by epifluorescence microscopy using an Olympus BX41® microscope fitted with a CCD Olympus DP-71® digital camera. Image capture was performed using DP controler® software (Olympus). The probes with positive hybridization signals were purified and cloned using the pMOSBlue blunt ended RPN5110 cloning kit (Amersham Biociences®) for subsequent sequencing. The samples were sequenced in an ABI-PRISM 3100 Genetic Analyzer automatic sequencer at the ACTGene laboratory (Centro de Biotecnologia, UFRGS, Porto Alegre, RS).

The sequences were aligned and edited using the CLUSTAL W program [65] using the following parameters: weights 6.66 and 10.0 for opening and extension of gaps, respectively, for pairwise alignments, and weights 10.0 and 15.0 for opening and extension of gaps in the multiple alignments respectively. Clustering analysis was performed using the parsimony method on the PAUP program v. 4.0b10 [70]. The isolated fragments of repetitive DNA sequences (Cot-1) from G. sylvius and G. paraguensis had lengths between 100 and 400 bp, respectively (Figure 5). These fragments were used as probes for fluorescent in situ hybridization (FISH) and were found to label most of the pericentromeric regions of the chromosomes of G. sylvius (Figure 4a) and all of the chromosomes of G. paraguensis (Figure 4b).

References

  1. 1. Heslop-Harrison JSP, Schwarzacher T. Organisation of the plant genome in chromosomes. Plant J. 2011;66:18-33. DOI:10.1111/j.1365-313X.2011.04544.x
  2. 2. Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Gen. Res. 2005; 110:462-467. DOI: 10.1159/000084979.
  3. 3. Grady DL, Ratliff RL, Robinson DL, McCanlies EC, Meyne J, Moyzis RK. Highly conserved repetitive DNA sequences are present at human centromeres. Proc Natl. Acad. Sci. USA. 1992; 89:1695-1699.
  4. 4. Schueler MG, Higgins AW, Rudd MK, Gustashaw K, Willard H. Genomic and genetic definition of a functional human centromere. Science. 2001;294:109-115. DOI: 10.1126/science.1065042.
  5. 5. Garrido-Ramos MA, Jamilena M, Lozano R, Ruiz Rejón C, Ruiz Rejón M. Cloning and characterization of a fish centromeric satellite DNA. Cytogenetic and Genome Research. 1994;65:233-237. DOI:10.1159/000133637.
  6. 6. Martins C, Ferreira IA, Oliveira C, Foresti F, Galetti Jr. PM. A tandemly repetitive centromeric DNA sequence of the fish Hoplias malabaricus (Characiformes: Erythrinidae) is derived from 5S rDNA. Genetica. 2006;127: 133-141. DOI: 10.1007/s10709-005-2674-y.
  7. 7. Mazzuchelli J, Martins C, Genomic organization of repetitive DNAs in the cichlid fish Astronotus ocellatus. Genetica. 2009;136: 461-469. DOI: 10.1007/s10709-008-9346-7.
  8. 8. CioffiMB, Martins C, Bertollo LAC. Evidence spreading of associated transposable elements and ribosomal DNA in the fish Erythrinus erythrinus. Implications for genome change and karyoevolution in fish. BMC Evolutionary Biology. 2010;10: 271-280. DOI: 10.1186/1471-2148-10-271.
  9. 9. Silva M, Ribeiro ED, Matoso DA, Sousa LM, Hrbek T, Rapp Py-Daniel L. Feldberg, E. Chromosomal polymorphism in two species of Hypancistrus (Siluriformes: Loricariidae): an integrative approach forunderstanding their biodiversity. Genetica. 2014;142(2):127-139. DOI: 10.1007/s10709-014-9760-y.
  10. 10. Cioffi MB, Bertollo LAC. Initial steps in XY chromosome differentiation in Hoplias malabaricus and the origin of an X1X2Y sex chromosome system in this fish group. Heredity. 2010;105:554-561. DOI:10.1038/hdy.2010.18.
  11. 11. Borba RS, Silva EL, Parise-Maltempi PP. Chromosome mapping of retrotransposable elements Rex1 and Rex3 in Leporinus Spix, 1829 species (Characiformes: Anostomidae) and its relationships among heterochromatic segments and W sex chromosome. Mobile Genet Elements. 2013;3(6):e27460. DOI: 10.4161/mge.27460.
  12. 12. Mariotto S, Centofante L, Artoni RF, Vicari MR, Moreira-Filho O. Chromosomal diversification in ribosomal DNA sites in Ancistrus species (Loricariidae: Ancistrini) from three hydrographic basins of Mato Grosso, Brazil. Comparative Cytogenetics. 2011;5:289-300. DOI:10.1007/s1116001192159.
  13. 13. Silva M, Matoso DA, Artoni RF, Feldberg E. Karyotypic diversity and evolutionary trends in Neotropical electric fish of the genus Gymnotus (Gymnotiformes: Gymnotidae). Zebrafish. 2019;16:308-320. DOI: 10.1089/zeb.2018.1716.
  14. 14. Vicari MR, Artoni RF, Moreira-Filho O, Bertollo LAC. Colocalization of repetitive DNAs and silencing of major rRNA genes. A case report of the fish Astyanax janeiroensis. Cytogenetics Genome Research. 2008;122:67-72. DOI: 10.1159/000151318.
  15. 15. Ferreira IA, Martins C. Physical chromosome mapping of repetitive DNA sequences in Nile tilapia Oreochromis niloticus: evidences for a differential distribution of repetitive elements in the sex chromosomes. Micron. 2008;39:411-418. DOI:10.1016/j.micron. 2007.02.010
  16. 16. Silva FA, Guimarães EC, Carvalho NDM, Ferreira AMV, Schneider CH, Carvalho-Zilse, GA, Feldberg E, Gross MC. Transposable DNA elements in amazonian fish: from genome enlargement to genetic adaptation to stressful Environments. Cytogenetics Genome Research. 2020;10. DOI: 10.1159/000507104
  17. 17. Almeida-Toledo LF, Foresti F, Pe’Quignot EV, Daniel-Silva MFZ. XX/XY sex chromosome system with X heterochromatinization: an early stage of sex chromosome differentiation in the Neotropic electric eel Eigenmannia virescens. Cytogenet. Cell Genet. 2001;95:73-78. DOI: 10.1159/000057020.
  18. 18. Margarido VP, Bellafronte E, Moreira-Filho O. Cytogenetic analysis of three sympatric Gymnotus (Gymnotiformes, Gymnotidae) species verifies invasive species in the Upper Paraná River basin, Brazil. Journal of Fish Biology. 2007;70:155-164. DOI: 10.1111/j.1095-8649.2007.01365.x.
  19. 19. Milhomem SSR, Pieczarka JC, Crampton WGR, Silva DS, Souza ACP, Carvalho JR, Nagamachi CY. Chromosomal evidence for a putative cryptic species in the Gymnotus carapo species-complex (Gymnotiformes, Gymnotidae). BMC Genetics. 2008; 9:75. DOI: 10.1186/1471-2156-9-75.
  20. 20. Nagamachi CY, Pieczarka JC, Milhomem SSR, O’Brien PCM, Souza ACP, Ferguson-Smith MA. Multiple rearrangements in cryptic species of electric knifefish, Gymnotus carapo (Gymnotidae, Gymnotiformes) revealed by chromosome painting. BMC Genetics. 2010;11:2-9. DOI: 10.1186/1471-2156-11-28.
  21. 21. Utsunomia R, Melo S, Scacchetti PC, Oliveira C, Machado MA, Pieczarka JC, Nagamachi CY, Foresti F. Particular chromosomal distribution of microsatellites in five species of the genus Gymnotus (Teleostei, Gymnotiformes). Zebrafish. 2018;15:398-403. DOI:10.1089/zeb.2018.1570.
  22. 22. Silva M, Barbosa P, Artoni R F, Feldberg E. Evolutionarydynamics of 5S rDNA and recurrent association of transposable elementsinelectric fish of the family Gymnotidae (Gymnotiformes): the case of Gymnotus mamiraua. Cytogenet. Genome Res. 2016;149:297-303. DOI: 10.1159/000449431.
  23. 23. Fernades-Matioli FMC, Almeida-Toledo LF. A molecular phylogenetic analysis in Gymnotus species (Pisces: Gymnotiformes) with inferences on chromosome evolution. Caryologia 2001;54(1):23-30.
  24. 24. Silva EB, Margarido VP. An X1X1X2X2/X1X2Y multiple sex chromosome system in a new species of the genus Gymnotus (Pisces, Gymnotiformes). Environmental Biology of Fishes. 2005;73:293-297. DOI: 10.1007/s10641-005-2144-5.
  25. 25. Lundberg JG, Marsall LG, Horton B, Malabarba CSL, Wesselingh F. The stage for Neotropical fish diversification: A history of tropical South American rivers. In: Malabarba LR, Reis RE, Vari RP, Lucena CAS, Lucena ZMS. (eds) Phylogeny and classification of neotropical fishes. EDIPURCS, Porto Alegre, RS, 1998. p. 13-48.
  26. 26. Nelson, J. Fishes of the World. Hoboken, NJ, John Wiley and Sons. 2006.
  27. 27. Crampton WGR, Thorsen DH, Albert JS. Three new species from a diverse and sympatric assemblage of the electric fish Gymnotus (Ostariophysi: Gymnotidae) in the lowland Amazon Basin, with notes on ecology. Copeia. 2005:82-99.
  28. 28. Campos-da-Paz R. Family Gymnotidae. In: Reis RE, Kullander SO, Ferraris Jr. CJ. (eds.). Check List of the Freshwater Fishes of South and Central. EDIPURCS, Porto Alegre, 2003. p. 483-486.
  29. 29. Albert JS, Miller RR. Gymnotus maculosus: a new species of electric fish from Middle America (Teleostei: Gymnotoidei), with a key to the species of Gymnotus. Proceedings Biological Society Washington. 1995;108:662-678.
  30. 30. Albert JS, Campos-da-Paz R. Phylogenetic systematics of America knifefishs: a review of the available data. In: Malabarba L, Reis RE, Vari RP, Lucena CAS, Lucena ZMS. (ed) Phylogeny and Classification of Neotropical Fishes. EDIPUCRS, Porto Alegre, 1998. p. 419-446.
  31. 31. Oyakawa OT, Akama A, Mautari KC, Nolasco JC. Peixes de riachos da Mata Atlântica nas Unidades de Conservação do Vale do Rio Ribeira de Iguape no Estado de São Paulo. São Paulo. Ed. Neotrópica. 2006.
  32. 32. Alves-Gomes JA. The Mitochondrial Phylogeny of the South American Electric Fish (Teleostei: Gymnotiformes) and an Alternative Hypothesis for the Otophysan Historical Biogeography. In: Poyato-Ariza FJ, Grande T, Diogo R (ed.) Teleostean Fish Biology: Gonorynchiformes and Ostariophysan Relationships: A Comprehensive Review. Enfield, Science Publishers, 2009. p. 517-565.
  33. 33. Milhomem SSR, Crampton WGR, Pieczarka JC, Shetka GH, Silva DS Nagamachi, CY. Gymnotus capanema, a new species of electric knife fish (Gymnotiformes, Gymnotidae) from Eastern Amazonia, with comments on an unusual karyotype. Journal Fish Biology. 2012;80:802-815. DOI:10.1111/j.1095-8649.2012.03219.x.
  34. 34. Milhomem SSR, Pieczarka JC, Crampton WR, Souza ACP, Carvalho Jr. JR, Nagamachi CY. Differences in karyotype between two sympatric species of Gymnotus (Gymnotiformes:Gymnotidae) from the eastern Amazon of Brazil. Zootaxa. 2007;1397:55-62. DOI:10.1186/1471-2156-9-75.
  35. 35. Silva M, Matoso DA, Vicari MR, Almeida MC, Margarido VP, Artoni RF. Repetitive DNA and meiotic behavior of sex chromosomes in Gymnotus pantanal (Gymnotiformes, Gmynotidae). Cytogenetic and Genome Research. 2011;135:143-149. DOI:10.1159/000330777.
  36. 36. Silva M, Matoso DA, Vicari MR, Almeida MC, Margarido VP, Artoni RF. 2011. Physical mapping of 5S DNAr in two species of Knifefishes: Gymnotus pantanal and Gymnotus paraguensis (Gymnotiformes). Cytogenetic Genome Research. 2011;134:303-307. DOI:10.1159/000328998.
  37. 37. Scacchetti PC, Alves JC, Utsunomia R, Claro FL, Almeida-Toledo LF, Oliveira C, Foresti F. Molecular characterization and physical mapping of two classes of 5S rDNA in genomes of Gymnotus sylvius and G. inaequilabiatus (Gymnotiformes, Gymnotidae). Cytogenetics Genome Research. 2011;136:1-7. DOI:10.1159/000335658.
  38. 38. Scacchetti PC, Pansonato-Alves JC, Utsunomia R, Oliveira C, Foresti F. Karyotypic diversity in four species of the genus Gymnotus Linnaeus, 1758 (Teleostei, Gymnotiformes, Gymnotidae): physical mapping of ribosomal genes and telomeric sequences. Comparative Cytogenetics. 2012;5:223-235. DOI:10.3897/CompCytogen.v5i3.1375.
  39. 39. Sochorová J, Garcia S, Gálvez F. et al. Evolutionary trends in animal ribosomal DNA loci: introduction to a new online database. Chromosoma. 2018;127:141-150. DOI:10.1007/s00412-017-0651-8.
  40. 40. Sivasundar A, Eldregde B, Orti G. Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers. Molecular Ecology. 2001;10:407-417. DOI:10.1046/j.1365-294x.2001.01194.x.
  41. 41. Ribeiro AC, Jacob RM, Silva RRRS, Lima FCT, Ferreira DC, Ferreira KM, Mariguela TC, Pereira LHG, Oliveira C. Distribuitions and phylogeographic data of rheophilic freshwater fishes provide evidence on the geographic extension of a Central – Brazilian Amazonian palaeo plateau in the area of the present day Pantanal wetland. Neotropical Ichthyology.2013;11:319-326. DOI:10.1590/S1679-62252013000200010.
  42. 42. Hoorn, C.; Wessilingh, F.P. Introduction: Amazonia, landscape and species evolution. In: Hoorn, C.; Wessilingh, F.P. (Eds). Amazonia: landscape and species evolution. A look into the past. Wiley-Blackwell, Oxford, 2010. p.1-9.
  43. 43. Montoya-Burgos JI. Historical biogeography of the catfish genus Hypostomus (Siluriformes: Loricariidae), with implications on the diversification of Neotropical ichthyofauna. Molecular Ecology. 2003;12:1855-1867. DOI:10.1046/j.1365-294X.2003.01857.x.
  44. 44. Silva M, Matoso DA, Vicari MR, Almeida MC, Margarido VP, Artoni RF. Physical mapping of 5S rDNA in two species of knifefishes: Gymnotus pantanal and Gymnotus paraguensis (Gymnotiformes). Cytogenetic and Genome Research. 2011;134,303-307. DOI:10.1159/000328998.
  45. 45. Zwick MS, Hanson RE, McKnight TD, Islam-Faridi HM, Stelly DM, Wing RA, Price JH. A rapid procedure for the isolation of C0t-1 DNA from plants. Genome. 1997;40:138-142.
  46. 46. Hřibová E, Doleželová M, Town CD, Macas J, Doležel J. Isolation and characterization of the highly repeated fraction of the banana genome. Cytogenet. Genome Res. 2007;119:268-274. DOI:10.1159/000112073.
  47. 47. Lorscheider CA, Oliveira JIN, Thais AD, Nogaroto V, Martins-Santos IC, Vicari MR. Comparative cytogenetics among three sympatric Hypostomus species (Siluriformes: Loricariidae): An evolutionary analysis in a high endemic region. Braz. Arch. Biol. Technol. 2018;61:14. DOI:10.1590/1678-4324-2018180417
  48. 48. Chalesworth B, Sniegowski P, Stephan W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994;371:215-220.
  49. 49. Topp CN, Dawe RK. Reinterpreting pericentromeric heterochromatin. Curr. Opin. Plant. Biol. 2006;9:647-653. DOI:10.1016/j.pbi.2006.09.008.
  50. 50. Valente GT, Schneider CH, Gross MC, Feldeberg E, Martins C. Comparative cytogenetics of cichlid fishes through genomic in-situ hybridization (GISH) with emphasis on Oreochromis niloticus. Chromosome Resarch. 2009. DOI:10.1007/s10577-009-9067-5.
  51. 51. Oliveira C, Wright JM. Molecular cytogenetic analysis of heterochromatin in the chromosomes of tilapia, Oreochromis niloticus (Teleostei: Cichlidae). Chromosome Research. 1998;6:205-211.
  52. 52. Oliveira C, Joyce S, Chew K, Porto-Foresti F, Dobson MJ, Wright JM. 1999. A LINE2 repetitive DNA sequence from the cichlid fish, Oreochromis niloticus: sequence analysis and chromosomal distribution. Chromosoma 1999;108:457-468.
  53. 53. Grewal SIS, Jia S. Heterochromatin revisited. Nat. Rev. Genet. 2007;8:35-46. DOI:10.1038/nrg2008.
  54. 54. Pimpinelli S, Berloco M, Fanti L, Dimitri P, Bonaccors S, Marchettib E, Caizzi R. Transposable elements are stable structural components of Drosophila melanogaster heterochromatin. Proc. Natl. Acad. Sci. USA. 1995;92:3804-3808.
  55. 55. Presting GG, Malysheva L, Fuchs J, Schubert I. A Ty3/gypsy retrotransposon-like sequence localizes to the centromeric regions of cereal chromosomes. Plant. J. 1998;6:721-728. DOI:10.1046/ j.1365-313x.1998.00341.x
  56. 56. Henning F, Trifonov V, Almeida-Toledo LF. Use of chromosome microdissection in fish molecular cytogenetic. Genetics and Molecular Biology. 2008:31:279-281. DOI:10.1590/S1415-47572008000200022.
  57. 57. Claro FL. Gymnotus carapo e Gymnotus sylvius (Teleostei: Gymnotidae): Uma abordagem Citogenética-Molecular. [Dissertation]. São Paulo - Brasil: Universidade de São Paulo; 2008.
  58. 58. Buldyrev SV, Goldberger AC, Harlin S, Mategna RN, Matsa ME, Peng CK, Simmons M, Stanley HE. Long-range correlation proprieties of coding and noncoding DNA sequences. Gene Bank. Analysis. Phys. Rev. E. 1995;51:5084-5091.
  59. 59. Melo S. Análises cromossômicas e genômicas aplicadas ao estudo evolutivo e estrutural dos satelitomas em espécies do gênero Gymnotus (Teleostei: Gymnotiformes). [Thesis]. Botucatu –SP: Universidade Estadual Paulista "Júlio de Mesquita Filho", Instituto de Biociências de Botucatu; 2020.
  60. 60. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN. DNA barcoding Australia's fish species. Philosophical Transactions of the Royal Society Biological Sciences. 2005;360:1847-1857. DOI:
  61. 61. Sambrook J, Russell DW. Molecular cloning: a laboratory manual. Cold Spring Harbor Laboratory Press, New York. 2001. p. 2-5.
  62. 62. Tamura K, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Molecular Biology and Evolution. 2007; 24, 1596-1599. DOI:10.1093/molbev/msm092.
  63. 63. Saitou N, Nei M. The Neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Evolution. 1987;4:406-425.
  64. 64. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution. 1980;16:111-120.
  65. 65. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research. 1994;22:4673-4680.
  66. 66. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/96/NT. Nucleic Acids Symposium Series. 1999;41:95-98.
  67. 67. Bertollo LAC, Takahash CS, Moreira-Filho O. Cytotaxonomic considerations on Hoplias lacerdae (Pisces, Erythrinidae). Rev. Bras. Genet. 1978;1(2):103-120.
  68. 68. Vicari MR, Nogaroto V, Noleto RB, Cestari MM, Cioffi MB, Almeida MC, Moreira-Filho O, Bertollo LAC, Artoni RF. Satellite DNA and chromosomes in Neotropical fishes: methods, applications and perspectives. Journal of Fish Biology. 2010;76:1094-1116. DOI: 10.1111/j.1095-8649.2010.02564.x.
  69. 69. Pinkel D, Straume T, Gray JW. Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization. Proc. Natl. Acad. Sci. USA. 1986;83:2934-2938.
  70. 70. Swofford DL. PAUP Phylogenetic analysis using parsimony (and other methods). Version 4. Sinauer Associates, Sunderland, Massachusetts. 2002. DOI:10.1111/j.0014-3820.2002.tb00191.x

Written By

Maelin da Silva, Daniele Aparecida Matoso, Vladimir Pavan Margarido, Eliana Feldberg and Roberto Ferreira Artoni

Submitted: 28 October 2020 Reviewed: 12 April 2021 Published: 07 May 2021