Are scattered microsatellites weak chromosomal markers? Guided mapping reveals new insights into Trachelyopterus (Siluriformes: Auchenipteridae) diversity

Abstract

The scattered distribution pattern of microsatellites is a challenging problem in fish cytogenetics. This type of array hinders the identification of useful patterns and the comparison between species, often resulting in over-limited interpretations that only label it as "scattered" or "widely distributed". However, several studies have shown that the distribution pattern of microsatellites is non-random. Thus, here we tested whether a scattered microsatellite could have distinct distribution patterns on homeologous chromosomes of closely related species. The clustered sites of 18S and 5S rDNA, U2 snRNA and H3/H4 histone genes were used as a guide to compare the (GATA)_n microsatellite distribution pattern on the homeologous chromosomes of six Trachelyopterus species: T. coriaceus and Trachelyopterus aff. galeatus from the Araguaia River basin; T. striatulus, T. galeatus and T. porosus from the Amazonas River basin; and Trachelyopterus aff. coriaceus from the Paraguay River basin. Most species had similar patterns of the (GATA)_n microsatellite in the histone genes and 5S rDNA carriers. However, we have found a chromosomal polymorphism of the (GATA)_n sequence in the 18S rDNA carriers of Trachelyopterus galeatus, which is in Hard-Weinberg equilibrium and possibly originated through amplification events; and a chromosome polymorphism in Trachelyopterus aff. galeatus, which combined with an inversion polymorphism of the U2 snRNA in the same chromosome pair resulted in six possible cytotypes, which are in Hardy-Weinberg disequilibrium. Therefore, comparing the distribution pattern on homeologous chromosomes across the species, using gene clusters as a guide to identify it, seems to be an effective way to further the analysis of scattered microsatellites in fish cytogenetics.

Introduction

Microsatellites, also known as short tandem repeats (STRs, [1]) or simple sequence repeats (SSRs, [2]), are stretches of DNA that consist of tandemly repeating di-, tri-, tetra- or penta-nucleotide motifs arranged through eukaryotic genomes [3,4] in both coding and noncoding regions [4,5]. Most microsatellites are located in the nucleus (nuSSR), although they can also be found in mitochondria (mtSSR) and chloroplasts (cpSSR) [4]. They are one of the most abundant and variable types of DNA sequences in the genome [3,5,6] and occur primarily due to slipped-strand mispairing and subsequent error(s) during DNA replication, repair or recombination [4,7]. However, the activity of transposable elements, mainly non-LTR retrotransposons, has also been reported as a major source of new microsatellites and their movement throughout the genome [4,8–11].

SSRs are recognized as powerful informative markers of genetic diversity and variability in both animal and plants [12–14]. In the last decades, they have been associated to several evolutionary and diversification aspects (see [3,10,15], including as potential causes of major structural chromosomal rearrangements, which seems to be facilitated by their high flexibility and low stability that creates fragile chromosomal sites [16,17]. In chromosomal mapping, although studies are still scarce for most species, they have provided a useful tool for understand the genome and chromosomal evolution of many different taxa [18–20]. They can be cytogenetically identified in large accumulations on a few chromosomal pairs [20,21] or with scattered signals throughout the chromosomes [20,22,23]. Furthermore, SSRs are commonly described associated to heterochromatic regions (e.g., [24,25]), often participating in its origin and increase (cHC) (for reviews, see [26–28]), and with crucial roles in the origin and evolution of specific chromosomes, mainly B chromosomes [22,23,29] and sex chromosomes [21,30,31]. However, in complement A chromosomes and populations without sex or B chromosomes, the cytogenetic studies are still concentrated only in the type of array and presence or absence of the marker, considering the whole karyotype [32,33].

One of the most challenging problems in expanding the use of microsatellites in fish cytogenetics is the scattered distribution pattern (e.g., [34–36], Table 1), because this type of array exhibit apparently random signals on the chromosomes. As a consequence, in most species without B or sex chromosomes, it hampers the description of useful distribution patterns to discuss evolutionary aspects as well as the comparison between karyotypes, often leading to oversimplified interpretations that simply label it as widely distributed or scattered. However, different lines of evidence indicated that the distribution of microsatellites is nonrandom [10]. They can influence in several aspects of the genome, including the nucleosome packing [37], methylation [38], high order chromatin structure [39–41] and splicing [42–44]. Microsatellites can also have enhancer functions [45–47], modulate gene expression [48–50] and participate of gene activity, recombination, DNA replication, cell cycle and mismatch repair (MMR) system [10,51]. Therefore, it is unlikely to believe that the scattered distribution pattern of SSRs in cytogenetic mappings is random and that it could not be helpful as chromosomal markers. In this scenario, considering that the scattered pattern consists in signals throughout the chromosomes, the association or colocalization with clustered sites in cytogenetic mapping, as already reported with several other types of DNA sequences [11,13,52], may provide an additional method to compare species and karyotypes. The clusters of previously mapped repetitive elements can be used as guide to identify homeologous chromosomes among species or populations, enabling comparative analyses beyond those with B or sex chromosomes.

[Figure omitted. See PDF.]

Auchenipteridae, known as the driftwood catfishes, is one of the 43 Siluriformes families. They are restricted but widely distributed in South America and have 25 genera and 128 valid species [86]. Among them, Trachelyopterus (with 17 valid species) is one of the most studied genera owing to high morphological similarity and, hence, a controversial taxonomic history [87]. The physical localization of 18S rDNA, 5S rDNA, H3/H4 histone genes and snRNA U2 is well documented in six species: T. coriaceus, T. striatulus, Trachelyopterus aff. galeatus (cited as Parauchenipterus galeatus) [88–90], T. porosus [90,91], T. galeatus (cited as P. galeatus), Trachelyopterus aff. coriaceus (cited as Trachelyopterus sp.) [23,90]; whereas the SSR (GATA)_n sites were described in only three species: T. galeatus [22,23], T. porosus [22] and Trachelyopterus aff. coriaceus (cited as Trachelyopterus sp.) [23]. In B chromosomes, the (GATA)_n sequence has a particularly clustered pattern, and provided valuable data about the origin and evolution of these extra chromosomes in Auchenipteridae. On the other hand, in complement A chromosomes it presented a scattered distribution pattern and was not able to reveal diversity or specific traits [22,23]. Therefore, these Trachelyopterus species constitutes an excellent model to be investigated through the integrated mapping, since it has (1) scattered signal of a microsatellite in all complement A chromosomes; (2) well known clusters of 18S rDNA, 5S rDNA, H3/H4 histone genes and U2 snRNA and; (3) they are closely related species with high similar karyotypes, which could test the integrated mapping efficiency.

In this paper the microsatellite (GATA)_n distribution pattern was compared among the carriers of the 5S and 18S rDNA, H3/H4 histone genes and U2 snRNA sites, in six Trachelyopterus species. The clustered repetitive elements sites were used as a guide to specific chromosomal pairs. We intended to test whether a widely distributed microsatellite could have distinct distribution patterns on homeologous chromosomes among closely related species, revealing unknown chromosome arrangements as well as architecture interaction between them.

Materials and methods

Ethics statement

Fish collections were authorized by Instituto Chico Mendes de Conservação da Biodiversidade (ICMBio, Permit number 49379), and the experimental procedures were approved by the Ethics Committee on Animal Experimentation and Practical Classes at Unioeste (09/13-CEEAAP/Unioeste). The use of the genetic data in this study was authorized by the National System for the Management of Genetic Heritage and Associated Traditional Knowledge (SisGen n° AFFBCC7).

Specimens and DNA extraction

Four Trachelyopterus species, including two possible new ones, were analyzed: Trachelyopterus galeatus Linnaeus, 1766 (04 males and 07 females) from Catalão lake, Manaus, Amazon River basin, 03°09’47’’S; 59°54’2’’W (Instituto Nacional de Pesquisas na Amazônia—INPA 57939); Trachelyopterus coriaceus Valenciennes, 1840 (04 males and 03 females) from Araguaia River, São Miguel do Araguaia, 13°08’52,7’’S; 50°25’02,8’’W (Museu da Universidade de São Paulo—MZUSP 106766); Trachelyopterus porosus Eigenmann and Eigenmann, 1888 (04 males and 04 females) from Catalão lake, Manaus, Amazon River basin, 03°09’47’’S; 59°54’29"W (Instituto Nacional de Pesquisas da Amazônia—INPA 57940); Trachelyopterus aff. galeatus (06 males and 10 females) from Araguaia River (suggested as new species), São Miguel do Araguaia, 13°08’52,7’’S; 50°25’02,8’W (MZUSP 110803); Trachelyopterus striatulus Steindachner, 1877 (03 males and 03 females) from Verde lagoon, Doce River basin, Marliéria, 19°49’44,5’’S; 42°37’52.5’’W (MZUSP 109798); and Trachelyopterus aff. coriaceus (02 males and 01 female) from Arrombado lagoon (suggested as new species), Bento Gomes River basin, Poconé, 16°25’40,9’’S; 56°25’07,4’’W (MZUSP 110806). Genomic DNA was extracted from the liver of all species according to Sambrook et al. (2001).

Probes labeling and sequencing

The 18S rDNA and 5S rDNA probes were obtained through polymerase chain reaction—PCR [90] using the primers NS1 and NS8 [92], 5SA and 5SB [93]. The 18S rDNA probes were labeled with Biotin-16-dUTP (Bio-Nick-Translation, Roche) according to manufacturer’s instructions and detected with Antibiotin-Avidin-FITC / Antiavidin-Biotin (Roche). The 5S rDNA were labeled through PCR with Tetramethyl-Rhodamine-5-dUTP (Roche). The H3 and H4 histone genes fragments were amplified from Trachelyopterus galeatus genome using the primers H3F and H3R [94]; and H4F2s and H4F2er [95]. Both probes were labeled through PCR with Tetramethyl-Rhodamine-5-dUTP (Roche) [90]. The U2 snRNA fragments were amplified from Trachelyopterus galeatus genome with the primers U2F and U2R [96]. The U2 snRNA probes were labeled through PCR with Fluorescein-12-dUTP (Roche) [90]. The (GATA)_n probes were generated through PCR in the absence of a template, using only the primers [97] and labeled with Tetramethyl-Rhodamine-5-dUTP (Roche).

All PCR products were sequenced in both ways, forward and reverse, using the ABI 3730 DNA Analyzer with the BigDye Terminator v3.1 Cycle Sequencing Kit (code 4337456) and the Sequencing Analysis software 5.3.1. The consensus sequence was generated by the Bioedit Sequence Alignment Editor [98]. The fragments identity was confirmed through BLASTn 2.11.0 (National Center for Biotechnology Information) [99].

Cytogenetic analyses and Fluorescent in situ hybridization

The samples were treated with a 0.02% colchicine solution (1 ml/100g of body weight) for 30–40 min and sequentially euthanized by clove oil overdose [100] (according to the ethics committee on animal experimentation and practical classes at Unioeste: 09/13—CEEAAP / Unioeste). The mitotic chromosomes were obtained from anterior kidney cells [101]. The chromosome morphology was classified according to [102].

Fluorescent in situ hybridization (FISH) was carried out with 77% of stringency [103] with some suggested modifications [104]. Clusters of 18S rDNA, 5S rDNA, snRNA U2, H3 and H4 histone genes were used as guide to integrative (guided) mapping with the SSR (GATA)_n. The digital images were captured by the DP Controller 3.2.1.276 software using an Olympus DP71 digital camera connected to the BX61 epifluorescence microscope (Olympus America Inc., Center Valley, PA, United States of America). The Hardy–Weinberg equilibrium (HWE) and Chi-squared test were performed using the Hardy-Weinberg (HW) testing program [105].

Results

DNA sequencing analysis

The U2 snRNA and both H3 and H4 histone genes were sequenced with forward and reverse primers. The U2 snRNA had a 199-bp consensus sequence (Genbank ID: OK166555.1) with 96,4% similarity with Parabotia fasciatus U2 snRNA gene (Genbank ID: MG874999.1). The histone gene H3 had a 448-bp consensus sequence (Genbank ID: OK166556.1) with 91.57% similarity with the Pimelodus microstoma H3 histone gene (Genbank ID: MT094432.1). The histone gene H4 had a 269-bp consensus sequence (Genbank ID: OK166557.1) with 90.09% similarity with the Pimelodus microstoma H4 histone gene (Genbank ID: MT094433.1).

Cytogenetic analysis

Trachelyopterus aff. galeatus–Araguaia River basin.

This species had 2n = 58 chromosomes for both sexes (7 males and 9 females). The microsatellite (GATA)_n were found spread throughout the chromosomes. The chromosomal pair 24, which bears the 18S rDNA, H3 and H4 histone genes loci, exhibited two blocks of the SSR (GATA)_n in the terminal position of the long arm and none in the short arm (Figs 1 and 4). The chromosomal pair 25, which bears another H3 and H4 histone genes loci, also had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 4). The chromosomal pair 26, which bear the U2 snRNA locus, exhibited a chromosome polymorphism associated to the SSR (GATA)_n distribution pattern (Fig 2), which combined with the U2 snRNA distribution pattern evidenced three chromosomal forms, including a new one, here referred as C chromosomal form: (A) U2 snRNA sites in the short arm with SSR (GATA)_n blocks only in the long arm of the chromosome; (B) U2 snRNA sites in the long arm with SSR (GATA)_n blocks only in the long arm of the chromosome; and (C) U2 snRNA sites in the short arm with SSR (GATA)_n blocks in both the short and long arm of the chromosome. From all chromosomal forms, six combinations could be hypothesized and four of them were found in the sample: AA (1 female), AB (2 males and 3 females), AC, BB (1 male and 5 females), BC, CC (2 males and 2 females) (Figs 2 and 5, Table 2). The χ2 value for Hardy-Weinberg equilibrium (HWE) was 15.99 (Df = 03; p = <0.05). The chromosomal pair 3 (5S rDNA) had three blocks of the SSR (GATA)_n, two in the short arm (one in the terminal position and one in interstitial position) and one in the terminal position of the long arm (Fig 3).

[Figure omitted. See PDF.]

The 18S rDNA probes are in green signal and the (GATA)n probes in red signal.

[Figure omitted. See PDF.]

U2 snRNA probes in green signal and (GATA)n probes in red signal. A1A1: Both chromosomes with U2 snRNA sites in the short arm and SSR (GATA)n only in the long arm; A1A2: One chromosome with U2 snRNA site in the long arm and another chromosome with the U2 snRNA site in the short arm, both with the SSR (GATA)n only in the long arm of the chromosome; A2A2: Both chromosomes with U2 snRNA sites in the short arm and the SSR (GATA)n only in the long arm; A3A3: Both chromosomes with U2 snRNA sites in the short arm of the chromosome and the SSR (GATA)n in the long and short arm of the chromosome.

[Figure omitted. See PDF.]

The 5S rDNA probes are in green signal and the (GATA)n probes in red signal.

[Figure omitted. See PDF.]

Trachelyopterus galeatus–Amazon River basin.

This species had 2n = 58 chromosomes for both sexes (7 males and 9 females). The chromosomal pair 20, which bears the 18S rDNA, H3 and H4 histone genes loci, has shown a chromosome polymorphism (Fig 1). It consists in two different organizations of the SSR (GATA)_n in the chromosome: two SSR (GATA)_n blocks in the terminal and interstitial position of long arm of the chromosome (A form); or two SSR (GATA)_n blocks in the terminal position and interstitial position of the long arm and one in the interstitial position of the short arm of the chromosome (B form). Three combinations of these polymorphic forms were found in this population: (AA) both chromosomes with two SSR (GATA)_n blocks in the terminal and proximal position of long arm of the chromosome (2 males and 1 female); (AB) one chromosome with two SSR (GATA)_n blocks in the terminal and proximal position of long arm of the chromosome and another chromosome with two SSR (GATA)_n blocks in the terminal and proximal position of the long arm of the chromosome and one in the proximal position of the short arm of the chromosome (2 males and 5 females); (BB) both chromosomes with two SSR (GATA)_n blocks in the terminal and proximal position of the long arm of the chromosome and one in the proximal position of the short arm of the chromosome (1 male). The chromosomal pair 28, bearing the U2 snRNA locus, had two blocks of the SSR (GATA)_n in the long arm and one in terminal position of the short arm. The chromosome polymorphism was found in Hardy-Weinberg equilibrium (HWE) (χ2 = 1.099; Df = 01; p = <0.05). The chromosomal pair 24, also bearing a H3 and H4 histone genes loci, had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 2). The chromosomal pair 14, bearing the first 5S rDNA locus, exhibited three blocks of the SSR (GATA)_n, two in the short arm (one in the terminal position and one in interstitial position) and one in the terminal position of the long arm (Fig 4). The chromosomal pair 16, bearing the second 5S rDNA locus, had two blocks of the SSR (GATA)_n, one in the terminal position of the short and another in the long arm (Fig 3).

[Figure omitted. See PDF.]

The H3/H4 histone genes probes are in red signal, The 18S rDNA and (GATA)n probes in green signal.

Trachelyopterus porosus–Amazon River basin. This species had 2n = 58 chromosomes for both sexes (4 males and 4 females). The microsatellite (GATA)_n were found spread throughout the chromosomes. The chromosomal pair 23, which bears the 18S rDNA, H3 and H4 histone genes loci, had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Figs 1 and 4). The chromosomal pair 24, also bearing a H3 and H4 histone genes loci, had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 4). The chromosomal pair 26, which bears the U2 snRNA locus, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 2). The chromosomal pair 3, which bears the first 5S rDNA locus, exhibited had two blocks of the SSR (GATA)_n, one in the terminal position of the short and another in the long arm. The chromosomal pair 4, bearing the second 5S rDNA locus, had three blocks of the SSR (GATA)_n, two in the short arm (one in the terminal position and one in interstitial position) and one in the terminal position of the long arm (Fig 3).

Trachelyopterus coriaceus–Araguaia River basin. This species had 2n = 58 chromosomes for both sexes (4 males and 3 females). The microsatellite (GATA)_n were found spread throughout the chromosomes. The chromosomal pair 23, which bears the 18S rDNA, H3 and H4 histone genes loci, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Figs 1 and 4). The chromosomal pair 28, bearing the U2 snRNA loci, had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 2). The chromosomal pair 3, which bears the first 5S rDNA locus, had three blocks of the SSR (GATA)_n, two in the short arm (one in the terminal position and one in interstitial position) and one in the terminal position of the long arm. The chromosomal pair 16, bearing the second 5S rDNA locus, had two blocks of the SSR (GATA)_n, one in the terminal position of the short and another in the long arm (Fig 3).

Trachelyopterus aff. coriaceus–Paraguay river basin.

This species had 2n = 58 chromosomes for both sexes (2 males and 1 female). The microsatellite (GATA)_n were found spread throughout the chromosomes. The chromosomal pair 22, bearing the 18S rDNA locus, had two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 1). The chromosomal pair 23, which bears the H3 and H4 histone genes loci, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 4). The chromosomal pair 27 bearing the U2 snRNA locus, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 2). The chromosomal pair 16, bearing the first 5S rDNA locus, had two blocks of the SSR (GATA)_n, one in the terminal position of the short and another in the long arm. The chromosomal pair 18, which bears the second 5S rDNA locus, had three blocks of the SSR (GATA)_n, two in the short arm (one in the terminal position and one in interstitial position) and one in the terminal position of the long arm (Fig 3).

Trachelyopterus striatulus–Doce River basin.

This species had 2n = 58 chromosomes for both sexes (3 males and 3 females). The microsatellite (GATA)_n were found spread throughout the chromosomes. The chromosomal pair 23, which bears the 18S rDNA, H3 and H4 histone genes loci, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Figs 1 and 4). The chromosomal pair 18, also bearing a H3 and H4 histone genes loci, had had two blocks of the SSR (GATA)_n; however, one in terminal position of the long arm and another in the interstitial position of the short arm (Fig 4). The chromosomal pair 28, which bears the U2 snRNA locus, exhibited two blocks of the SSR (GATA)_n in the long arm and none in the short arm (Fig 2). The chromosomal pair 10, bearing the first 5S rDNA locus, had three blocks of the SSR (GATA)_n, two in the long arm (one in the terminal position and one in interstitial position) and one in the terminal position of the short arm. The chromosomal pair 13 and 15, which bears another two 5S rDNA locus, exhibited two blocks of the SSR (GATA)_n, one in the terminal position of the short and another in the long arm (Fig 3).

Discussion

The (GATA)_n repeats, molecular components of the Bkm satellite DNA (Banded krait minor-satellite DNA), are widely distributed in the genome of higher organisms [106]. As expected, in all Trachelyopterus species of this study it was found scattered throughout the chromosomes (summarized in Fig 5). This type of array and has been described to different SSRs and taxa, such as plants [107,108], amphibians [20], fungus [109] and fish [22,23,33]. In Auchenipteridae, it was already reported in T. porosus and T. galeatus [22], Trachelyopterus aff. coriaceus (cited as Trachelyopterus sp. [23], and now it can also be seen in T. coriaceus, T. striatulus and Trachelyopterus aff. galeatus. This type of scattered array could be explained by the activity of transposable elements, which are a major source of new microsatellites and can also drive it throughout the genome [8–11] and/or to chromosomal rearrangements, as already proposed for a close genus, Hypostomus [33]. Though the (GATA)_n repeats have been associated to sex chromosome in different organisms [21,110–112], no differences between males and females have been found in these populations.

[Figure omitted. See PDF.]

In the integrated mapping with the 18S rDNA, the chromosomal pair of all species had almost the same SSR (GATA)_n distribution pattern, characterized by the absence of the microsatellite only in the short arm (Fig 1). It was also reported in a close species, Glanidium ribeiroi [53], but it was unexplored yet. Usually, SSRs can constitute a larger fraction of noncoding DNA, but are rare in protein-coding regions [10] mainly due to negative selection against frameshift mutations in coding regions, as evidenced in plants, primates, and microorganisms [113]. Comparative studies showed that only repeats in multiples of three may develop evenly in both regions [5,113], since RNA bases are read as triplets and other types could result in frameshift mutations [113,114]. In this way, the absence of the (GATA)_n sequence in the short arm of the 18S rDNA chromosomes, might suggest a negative impact of the microsatellite near the coding areas of these species.

However, this hypothesis is contradicted by some evidences: (a) T. galeatus from the Amazon River basin had a (GATA)_n block in the short arm of the 18S rDNA carrier; (b) there is overlaid signal between the (GATA)_n sequence and the 5S rDNA and U2 snRNA sites; and (c) the (GATA)_n sequence is distributed throughout the chromosomes, indicating that it could also be near other unmapped gene sequences. Therefore, the non-existence of (GATA)_n signal in the short arm might only be related to a spatial issue, since besides the 18S rDNA, the short arm of most species also carries the H3 and H4 histone genes and, consequently, there could not be enough space for large amounts of (GATA)_n sequence detectable through Fluorescent in situ Hybridization, which needs targets of at least 1kb to express significative results [115].

In contrast to the other species of our sample, T. galeatus from the Amazon River basin was the only one to possess the (GATA)_n sequence in the short arm of the 18S rDNA carrier. However, it was identified in a polymorphic state, which appears to be neutral, since the heterozygous has no different adaptive value compared to other forms (Table 2). Although other mechanisms have been proposed to explain how microsatellites can arise and expand over time, e.g., double-stranded DNA recombination (unequal crossing over and gene conversion), mismatch/double-strand break repair, and retrotransposition; replication slippage is still considered one of the main mechanisms thus far [3,4,116]. It also seems to be the most parsimonious to explain the origin of the microsatellite (GATA)_n in the short arm of the 18S rDNA carrier in T. galeatus from the Amazon River basin. In this case, a small number of repeats (proto-microsatellite) is required before DNA polymerase slippage can extend the number of repeats, originating a new microsatellite [114], which seems not to be a problem in the genome of these Trachelyopterus species, as it is highly enriched with the (GATA)_n repeats. Once generated, the new microsatellite can undergo several reduction or expansion events over time (see [3]), which can lead to the formation of large microsatellite blocks, such as the one visualized in T. galeatus from the Amazon River basin. On the other hand, the polymorphic state observed in this population may have arisen from a crossing between individuals possessing the new trait and those with the original condition.

The (GATA)_n distribution pattern in the U2 snRNA chromosomes was conserved among most species. However, Trachelyopterus aff. galeatus from Araguaia River basin presented a chromosomal polymorphism in the same chromosomal pair (26) that was also reported to have a U2 snRNA chromosomal polymorphism [90]. Combining both polymorphic markers, U2 snRNA and SSR (GATA)_n resulted in three chromosomal forms. The B chromosomal form is exclusive of this population; whereas the C chromosomal form was also found in T. galeatus from Amazon River basin and the and A form is present in T. striatulus, T. coriaceus, T. porosus, and Trachelyopterus aff. coriaceus. Although both polymorphic states interact to compose the chromosome arrangement in Trachelyopterus aff. galeatus, they seem to be originated in different evolutive events. The U2 snRNA polymorphism is an exclusive trait of Trachelyopterus aff. galeatus and product of a pericentric inversion [90]. On the other hand, the additional (GATA)_n block that characterizes the polymorphism in Trachelyopterus aff. galeatus is not an exclusive trait and can also be seen in T. galeatus from Amazon River. Thus, the most parsimonious hypothesis is that the (GATA)_n polymorphism in Trachelyopterus aff. galeatus might be originated through hybridization in secondary contact zones [117,118], which could be facilitated by the historical and geomorphological aspects of the Araguaia River floodplains, known by constant ichthyofaunistic exchange across surrounding hydrographic systems during the neotectonics reactivations in the Transbrasiliano Lineament during the formation of the Araguaia depression [119–121].

Interestingly, the U2 snRNA inversion polymorphism in Trachelyopterus aff. galeatus was reported in Hardy-Weinberg equilibrium [90], in which, the spread of the polymorphism could be associated to the neutrality of the rearrangement, since it suggests that there is no change in adaptative value among the genotypes or in the host fitness [122], as reported for water beetles [123] and blackflies [124]. In this state, the polymorphism is essentially influenced by genetic drift and migration [122]. However, the polymorphic state of U2 snRNA with the (GATA)_n sequence, resulted in Hardy-Weinberg disequilibrium (Table 2), suggesting that the combined arrangement may be under the effect of different forces beyond just genetic drift and gene flow.

Furthermore, all genotypes were found in a similar proportion compared to the expected by the Hardy-Weinberg equilibrium test (Table 2), except the ones with the C chromosomal form, in which, none of the heterozygous were found in the sample (AC and BC) and the homozygous (CC) presented a three times higher frequency than expected. In some cases, heterozygous originated from chromosome rearrangements can suffer severe reductions in fitness [125], zygotic lethality [126] or hybrid [127,128], especially when it involves change in gene order within a chromosome (inversions) [125] or when the hybrids carry multiple rearrangements [125,126]. In this scenario, the presence of multiple chromosomal polymorphism in the same chromosomal pair, which origin of both could be related to major chromosome rearrangements (inversions), associated to the absence of the C heterozygous and higher frequency of the C homozygous, may suggest the existence of distinct evolutive pressures over it compared to other genotypes. However, analyzes with a larger sample size, since the C form in heterozygosis may just not have been collected, are still needed to clarify it.

In contrast to the previously discussed markers, the (GATA)_n chromosomal mapping on the 5S rDNA and H3/H4 histone gene carriers did not reveal new information about the structure of these chromosomal pairs. However, using the (GATA)_n distribution pattern on the 5S rDNA carriers, the possible chromosome homeologies between the species could be inferred (Fig 3). Since the 5S rDNA is a marker usually found on multiple chromosomal pairs in Trachelyopterus, the chromosomal correspondence is difficult to suggest without additional information about the organization of each chromosome, a gap that could be partially filled with the distribution of the (GATA)_n sequence on these chromosomes. Through the (GATA)_n mapping on the 5S rDNA chromosome pairs, three main chromosomal arrangements could be evidenced among the species (Fig 3), in which, the chromosomal arrangement (1) is present in all species, whereas the chromosomal arrangement (2) is present in T. striatulus, T. galeatus, Trachelyopterus aff. coriaceus and T. coriaceus, and the chromosomal arrangement (3) is present only in T. striatulus and T. porosus.

The microsatellite mapping combined with the H3 and H4 histone genes confirmed the distribution of SSR (GATA)_n already pointed out through the integration with 18S rDNA for the species that have this synteny (T. striatulus, T. galeatus, Trachelyopterus aff. galeatus, T. porosus and T. coriaceus) as well as the polymorphism in the chromosomal pair 20 of T. galeatus from the Amazon River basin. No new arrangement could be observed, and even for species that have multiple sites of H3 and H4 histone genes (T. striatulus, T. galeatus, Trachelyopterus aff. galeatus and T. porosus), both chromosomal pairs showed the same microsatellite distribution pattern.

Integrated mapping perspectives to scattered microsatellites in Neotropical fishes

To date, there is no similar approach in cytogenetic mapping of microsatellites comparing Neotropical fish species. Most studies have focused in the origin and evolution of specific chromosomes, mainly B chromosomes (13 out of 81 studied species—16,04%) and sex chromosomes (24 out of 81 studied species– 29,62%); whereas others focused only on the type of array and presence or absence of the marker (Table 1). Nonetheless, of all cytogenetically analyzed species through physical mapping of microsatellites, 76.14% had at least one scattered microsatellite (67 out of 88 analyzed species), and in most studies, they could not be used to differ the species or populations.

Although the integrated mapping did not reveal large chromosomal rearrangements and most species had a similar distribution pattern of the (GATA)_n sequence in the analyzed chromosomes, it proved that the scattered microsatellite (GATA)_n has a non-random distribution, reiterating the existence of organization even in scattered microsatellites, which can be better described through a smaller scale analysis, comparing specific chromosomes between species. Through the integrated mapping a more accurate (GATA)_n pattern were described to the chromosome carriers of the 18S and 5S rDNA, H3 and H4 histone genes and U2 snRNA, and even though three species/populations used in this study were already mapped with (GATA)_n, none of them were able to detect the (GATA)_n polymorphism in T. galeatus from Amazonas River basin and the presence of three chromosomal arrangements of the (GATA)_n sequence in the 5S rDNA carrier. Likewise, the (GATA)_n polymorphism in the U2 snRNA chromosomes of Trachelyopterus aff. galeatus from the Araguaia River basin would possibly go unnoticed without using the U2 snRNA site as a guide. Therefore, the integrated mapping (guided) proved to be an efficient methodology to reveal cryptic chromosomal arrangements of scattered microsatellites.

Cytotaxonomically, the integrated mapping showed the divergence in one more marker between T. galeatus from the Amazon River basin and Trachelyopterus aff. galeatus from the Araguaia River, which reiterate the existence of a possible new species, as already proposed through ribosomal markers [89] and integrated mapping with H3/H4 histone genes and U2 snRNA [90]. Methodologically, the integrated mapping turned the widely distributed SSR (GATA)_n, with contributions in Trachelyopterus only to B chromosomes origin and evolution, into promising marker to distinguish other Auchenipteridae species. Thus, with the advancement of repetitive elements cytogenetic, we expect that the integrated mapping can further add to Trachelyopterus and also to other species with scattered distributed microsatellites.

Acknowledgments

We thank to Dr. Heraldo Antonio Britski and Dr. José Luís Olivan Birindelli for the identification of the specimens; Pedro Luis Gallo and Luis Henrique da Silva for their assistance in field activities; to Chico Mendes Institute for Biodiversity Conservation for the authorization to collect the animals; to Minas Gerais State Forestry Institute for the logistical support during activities in the Parque Estadual do Rio Doce; and Sr. Osmar Pereira de Barros, owner of the Ana Maria farm in São Miguel do Araguaia, for the support and authorization to collect the animals in his property.

Citation: Haerter CAG, Blanco DR, Traldi JB, Feldberg E, Margarido VP, Lui RL (2023) Are scattered microsatellites weak chromosomal markers? Guided mapping reveals new insights into Trachelyopterus (Siluriformes: Auchenipteridae) diversity. PLoS ONE 18(6): e0285388. https://doi.org/10.1371/journal.pone.0285388

About the Authors:

Chrystian Aparecido Grillo Haerter

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Instituto Nacional de Pesquisas da Amazônia, Manaus, Brasil

Daniel Rodrigues Blanco

Roles: Conceptualization, Data curation, Supervision, Validation, Visualization, Writing – review & editing

Affiliation: Universidade Tecnológica Federal do Paraná, Campus Santa Helena, Paraná, Brasil

Josiane Baccarin Traldi

Roles: Conceptualization, Data curation, Investigation, Validation, Visualization, Writing – review & editing

Affiliation: Departamento de Genética, Instituto de Ciências Biológicas, Universidade Federal do Amazonas, Manaus, Brasil

ORICD: https://orcid.org/0000-0003-3019-7441

Eliana Feldberg

Roles: Data curation, Validation, Visualization, Writing – review & editing

Affiliation: Instituto Nacional de Pesquisas da Amazônia, Manaus, Brasil

Vladimir Pavan Margarido

Roles: Conceptualization, Data curation, Funding acquisition, Resources, Supervision, Validation, Visualization, Writing – review & editing

Affiliation: Universidade Estadual do Oeste do Paraná, Centro de Ciências Biológicas e da Saúde, Cascavel, Paraná, Brasil

Roberto Laridondo Lui

Roles: Conceptualization, Data curation, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

E-mail: [email protected]

Affiliation: Universidade Estadual do Oeste do Paraná, Centro de Ciências Biológicas e da Saúde, Cascavel, Paraná, Brasil

ORICD: https://orcid.org/0000-0003-4310-4865

References

1. Edwards A, Civitello A, Hammond HA, Caskey CT. DNA typing and genetic mapping with trimeric and tetrameric tandem repeats. Am J Hum Genet. 1991;49(4):746–56. Epub 1991/10/01. pmid:1897522; PubMed Central PMCID: PMC1683171.

2. Jacob HJ, Lindpaintner K, Lincoln SE, Kusumi K, Bunker RK, Mao YP, et al. Genetic mapping of a gene causing hypertension in the stroke-prone spontaneously hypertensive rat. Cell. 1991;67(1):213–24. Epub 1991/10/04. pmid:1655275.

3. Ellegren H. Microsatellites: simple sequences with complex evolution. Nat Rev Genet. 2004;5(6):435–45. Epub 2004/05/22. pmid:15153996.

4. Kalia RK, Rai MK, Kalia S, Singh R, Dhawan AK. Microsatellite markers: an overview of the recent progress in plants. Euphytica. 2011;177:309–24. https://doi.org/10.1007/s10681-010-0286-9.

5. Toth G, Gaspari Z, Jurka J. Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000;10(7):967–81. Epub 2000/07/19. pmid:10899146; PubMed Central PMCID: PMC310925.

6. Weber JL. Informativeness of human (dC-dA)n.(dG-dT)n polymorphisms. Genomics. 1990;7(4):524–30. Epub 1990/08/01. pmid:1974878.

7. Levinson G, Gutman GA. Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol. 1987;4(3):203–21. Epub 1987/05/01. pmid:3328815.

8. Messier W, Li SH, Stewart CB. The birth of microsatellites. Nature. 1996;381(6582):483. Epub 1996/06/06. pmid:8632820.

9. Nadir E, Margalit H, Gallily T, Ben-Sasson SA. Microsatellite spreading in the human genome: evolutionary mechanisms and structural implications. Proc Natl Acad Sci U S A. 1996;93(13):6470–5. Epub 1996/06/25. pmid:8692839; PubMed Central PMCID: PMC39047.

10. Li YC, Korol AB, Fahima T, Beiles A, Nevo E. Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review. Mol Ecol. 2002;11(12):2453–65. Epub 2002/11/28. pmid:12453231.

11. Grandi FC, An W. Non-LTR retrotransposons and microsatellites: Partners in genomic variation. Mob Genet Elements. 2013;3(4):e25674. Epub 2013/11/07. pmid:24195012; PubMed Central PMCID: PMC3812793.

12. Peakall R, Gilmore S, Keys W, Morgante M, Rafalski A. Cross-species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: implications for the transferability of SSRs in plants. Mol Biol Evol. 1998;15(10):1275–87. Epub 1998/10/27. pmid:9787434.

13. Cuadrado A, Jouve N. Similarities in the chromosomal distribution of AG and AC repeats within and between Drosophila, human and barley chromosomes. Cytogenet Genome Res. 2007;119(1–2):91–9. Epub 2007/12/28. pmid:18160787.

14. Ernetti JR, Gazolla CB, Recco-Pimentel SM, Luca EM, Bruschi DP. Non-random distribution of microsatellite motifs and (TTAGGG)n repeats in the monkey frog Pithecopus rusticus (Anura, Phyllomedusidae) karyotype. Genet Mol Biol. 2020;42(4):e20190151. Epub 2020/01/23. pmid:31968045; PubMed Central PMCID: PMC7198017.

15. Li YC, Korol AB, Fahima T, Nevo E. Microsatellites within genes: structure, function, and evolution. Mol Biol Evol. 2004;21(6):991–1007. Epub 2004/02/14. pmid:14963101.

16. Raskina O, Barber JC, Nevo E, Belyayev A. Repetitive DNA and chromosomal rearrangements: speciation-related events in plant genomes. Cytogenet Genome Res. 2008;120(3–4):351–7. Epub 2008/05/28. pmid:18504364.

17. Barros AV, Wolski MA, Nogaroto V, Almeida MC, Moreira-Filho O, Vicari MR. Fragile sites, dysfunctional telomere and chromosome fusions: What is 5S rDNA role? Gene. 2017;608:20–7. Epub 2017/01/24. pmid:28111257.

18. Piscor D, Parise-Maltempi PP. Microsatellite Organization in the B Chromosome and A Chromosome Complement in Astyanax (Characiformes, Characidae) Species. Cytogenet Genome Res. 2016;148(1):44–51. Epub 2016/03/19. pmid:26992246.

19. Travenzoli NM, Lima BA, Cardoso DC, Dergam JA, Fernandes-Salomao TM, Lopes DM. Cytogenetic Analysis and Chromosomal Mapping of Repetitive DNA in Melipona Species (Hymenoptera, Meliponini). Cytogenet Genome Res. 2019;158(4):213–24. Epub 2019/07/29. pmid:31352441.

20. Bueno GP, Gatto KP, Gazolla CB, Leivas PT, Struett MM, Moura M, et al. 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). PLoS One. 2021;16(1):e0245128. Epub 2021/01/14. pmid:33439901; PubMed Central PMCID: PMC7806164.

21. Viana PF, Ezaz T, de Bello Cioffi M, Jackson Almeida B, Feldberg E. Evolutionary Insights of the ZW Sex Chromosomesin Snakes: A New Chapter Added by the AmazonianPuffing Snakes of the Genus Spilotes. Genes (Basel). 2019;10(4). Epub 2019/04/12. pmid:30970650; PubMed Central PMCID: PMC6523457.

22. Felicetti D, Haerter CAG, Baumgartner L, Paiz LM, Takagui FH, Margarido VP, et al. A New Variant B Chromosome in Auchenipteridae: The Role of (GATA)n and (TTAGGG)n Sequences in Understanding the Evolution of Supernumeraries in Trachelyopterus. Cytogenet Genome Res. 2021;161(1–2):70–81. Epub 2021/02/19. pmid:33601372.

23. Lui RL, Traldi JB, Blanco DR, Margarido VP, Mariotto S, Centofante L, et al. Possible Common Origin of B Chromosomes in Neotropical Fish (Siluriformes, Auchenipteridae) Reinforced by Repetitive DNA Mapping. Brazilian Archives of Biology and Technology. 2021;64:e21190494. https://doi.org/10.1590/1678-4324-2021190494.

24. Cuadrado A, Cardoso M, Jouve N. Physical organisation of simple sequence repeats (SSRs) in Triticeae: structural, functional and evolutionary implications. Cytogenet Genome Res. 2008;120(3–4):210–9. Epub 2008/05/28. pmid:18504349.

25. Viana PF, Feldberg E, Cioffi MB, de Carvalho VT, Menezes S, Vogt RC, et al. The Amazonian Red Side-Necked Turtle Rhinemys rufipes (Spix, 1824) (Testudines, Chelidae) Has a GSD Sex-Determining Mechanism with an Ancient XY Sex Microchromosome System. Cells. 2020;9(9). Epub 2020/09/17. pmid:32932633; PubMed Central PMCID: PMC7563702.

26. Trojer P, Reinberg D. Facultative heterochromatin: is there a distinctive molecular signature? Mol Cell. 2007;28(1):1–13. Epub 2007/10/16. pmid:17936700.

27. Allshire RC, Madhani HD. Ten principles of heterochromatin formation and function. Nat Rev Mol Cell Biol. 2018;19(4):229–44. Epub 2017/12/14. pmid:29235574; PubMed Central PMCID: PMC6822695.

28. Saksouk N, Simboeck E, Dejardin J. Constitutive heterochromatin formation and transcription in mammals. Epigenetics Chromatin. 2015;8:3. Epub 2015/03/20. pmid:25788984; PubMed Central PMCID: PMC4363358.

29. Milani D, Cabral-de-Mello DC. Microsatellite organization in the grasshopper Abracris flavolineata (Orthoptera: Acrididae) revealed by FISH mapping: remarkable spreading in the A and B chromosomes. PLoS One. 2014;9(5):e97956. Epub 2014/05/30. pmid:24871300; PubMed Central PMCID: PMC4037182.

30. Barbosa P, Leal EV, da Silva M, de Almeida MC, Moreira-Filho O, Artoni RF. Variability and evolutionary implications of repetitive DNA dynamics in genome of Astyanax scabripinnis (Teleostei, Characidae). Comp Cytogenet. 2017;11(1):143–62. Epub 2017/09/19. pmid:28919955; PubMed Central PMCID: PMC5599702.

31. Barcellos SA, Kretschmer R, de Souza MS, Costa AL, Degrandi TM, Dos Santos MS, et al. Karyotype Evolution and Distinct Evolutionary History of the W Chromosomes in Swallows (Aves, Passeriformes). Cytogenet Genome Res. 2019;158(2):98–105. Epub 2019/06/04. pmid:31158838.

32. Lui RL, Blanco DR, Martinez JF, Margarido VP, Venere PC, Moreira-Filho O. The role of chromosomal fusion in the karyotypic evolution of the genus Ageneiosus (Siluriformes: Auchenipteridae. Neotropical Ichthyology. 2013;11:327–34. https://doi.org/10.1590/S1679-62252013005000004.

33. Traldi JB, Blanco DR, Vicari MR, Martinez Jde F, Lui RL, Artoni RF, et al. Physical mapping of (GATA)n and (TTAGGG)n sequences in species of Hypostomus (Siluriformes, Loricariidae). J Genet. 2013;92(1):127–30. Epub 2013/05/04. pmid:23640416.

34. da Silva MJ, de Araujo Vieira AP, Galvao Cipriano FM, Dos Santos Candido MR, de Oliveira EHC, Gimenez Pinheiro T, et al. The Karyotype of Salvator merianae (Squamata, Teiidae): Analyses by Classical and Molecular Cytogenetic Techniques. Cytogenet Genome Res. 2020;160(2):94–9. Epub 2020/02/18. pmid:32062647.

35. Yeesin P, Buasriyot P, Ditcharoen S, Chaiyasan P, Suwannapoom C, Juntaree S, et al. Comparative study of four Mystus species (Bagridae, Siluriformes) from Thailand: insights into their karyotypic diversity. Comp Cytogenet. 2021;15(2):119–36. Epub 2021/05/08. pmid:33959235; PubMed Central PMCID: PMC8093182.

36. Saenjundaeng P, Supiwong W, Sassi FMC, Bertollo LAC, Rab P, Kretschmer R, et al. Chromosomes of Asian cyprinid fishes: Variable karyotype patterns and evolutionary trends in the genus Osteochilus (Cyprinidae, Labeoninae, "Osteochilini"). Genet Mol Biol. 2020;43(4):e20200195. Epub 2020/11/07. pmid:33156892; PubMed Central PMCID: PMC7783954.

37. Iyer V, Struhl K. Poly(dA:dT), a ubiquitous promoter element that stimulates transcription via its intrinsic DNA structure. EMBO J. 1995;14(11):2570–9. Epub 1995/06/01. pmid:7781610; PubMed Central PMCID: PMC398371.

38. Fukuda K, Ichiyanagi K, Yamada Y, Go Y, Udono T, Wada S, et al. Regional DNA methylation differences between humans and chimpanzees are associated with genetic changes, transcriptional divergence and disease genes. J Hum Genet. 2013;58(7):446–54. Epub 2013/06/07. pmid:23739127.

39. Pathak RU, Mamillapalli A, Rangaraj N, Kumar RP, Vasanthi D, Mishra K, et al. AAGAG repeat RNA is an essential component of nuclear matrix in Drosophila. RNA Biol. 2013;10(4):564–71. Epub 2013/04/17. pmid:23588056; PubMed Central PMCID: PMC3710362.

40. McNeil JA, Smith KP, Hall LL, Lawrence JB. Word frequency analysis reveals enrichment of dinucleotide repeats on the human X chromosome and [GATA]n in the X escape region. Genome Res. 2006;16(4):477–84. Epub 2006/03/15. pmid:16533911; PubMed Central PMCID: PMC1457025.

41. Subramanian S, Mishra RK, Singh L. Genome-wide analysis of Bkm sequences (GATA repeats): predominant association with sex chromosomes and potential role in higher order chromatin organization and function. Bioinformatics. 2003;19(6):681–5. Epub 2003/04/15. pmid:12691979.

42. Cuppens H, Lin W, Jaspers M, Costes B, Teng H, Vankeerberghen A, et al. Polyvariant mutant cystic fibrosis transmembrane conductance regulator genes. The polymorphic (Tg)m locus explains the partial penetrance of the T5 polymorphism as a disease mutation. J Clin Invest. 1998;101(2):487–96. Epub 1998/02/07. pmid:9435322; PubMed Central PMCID: PMC508589.

43. Hefferon TW, Groman JD, Yurk CE, Cutting GR. A variable dinucleotide repeat in the CFTR gene contributes to phenotype diversity by forming RNA secondary structures that alter splicing. Proc Natl Acad Sci U S A. 2004;101(10):3504–9. Epub 2004/03/03. pmid:14993601; PubMed Central PMCID: PMC373492.

44. Hui J, Hung LH, Heiner M, Schreiner S, Neumuller N, Reither G, et al. Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing. EMBO J. 2005;24(11):1988–98. Epub 2005/05/13. pmid:15889141; PubMed Central PMCID: PMC1142610.

45. Kumar RP, Krishnan J, Pratap Singh N, Singh L, Mishra RK. GATA simple sequence repeats function as enhancer blocker boundaries. Nat Commun. 2013;4:1844. Epub 2013/05/16. pmid:23673629.

46. Gebhardt F, Zanker KS, Brandt B. Modulation of epidermal growth factor receptor gene transcription by a polymorphic dinucleotide repeat in intron 1. J Biol Chem. 1999;274(19):13176–80. Epub 1999/05/01. pmid:10224073.

47. Gymrek M, Willems T, Guilmatre A, Zeng H, Markus B, Georgiev S, et al. Abundant contribution of short tandem repeats to gene expression variation in humans. Nat Genet. 2016;48(1):22–9. Epub 2015/12/08. pmid:26642241; PubMed Central PMCID: PMC4909355.

48. Chen TM, Kuo PL, Hsu CH, Tsai SJ, Chen MJ, Lin CW, et al. Microsatellite in the 3’ untranslated region of human fibroblast growth factor 9 (FGF9) gene exhibits pleiotropic effect on modulating FGF9 protein expression. Hum Mutat. 2007;28(1):98. Epub 2006/12/13. pmid:17154280.

49. Kumar S, Bhatia S. A polymorphic (GA/CT)n- SSR influences promoter activity of Tryptophan decarboxylase gene in Catharanthus roseus L. Don. Sci Rep. 2016;6:33280. Epub 2016/09/14. pmid:27623355; PubMed Central PMCID: PMC5020687.

50. Joshi-Saha A, Reddy KS. Repeat length variation in the 5’UTR of myo-inositol monophosphatase gene is related to phytic acid content and contributes to drought tolerance in chickpea (Cicer arietinum L.). J Exp Bot. 2015;66(19):5683–90. Epub 2015/04/19. pmid:25888598.

51. Bagshaw ATM. Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes. Genome Biol Evol. 2017;9(9):2428–43. Epub 2017/09/29. pmid:28957459; PubMed Central PMCID: PMC5622345.

52. Ruiz-Ruano FJ, Cuadrado A, Montiel EE, Camacho JP, Lopez-Leon MD. Next generation sequencing and FISH reveal uneven and nonrandom microsatellite distribution in two grasshopper genomes. Chromosoma. 2015;124(2):221–34. Epub 2014/11/13. pmid:25387401.

53. Lui RL, Blanco DR, Traldi JB, Margarido VP, Moreira-Filho O. Karyotypic variation of Glanidium ribeiroi Haseman, 1911 (Siluriformes, Auchenipteridae) along the Iguazu river basin. Braz J Biol. 2015;75(4 Suppl 1):S215–21. Epub 2015/11/26. pmid:26602347.

54. Prizon AC, Bruschi DP, Gazolla CB, Borin-Carvalho LA, Portela-Castro ALB. Chromosome Spreading of the Retrotransposable Rex-3 Element and Microsatellite Repeats in Karyotypes of the Ancistrus Populations. Zebrafish. 2018;15(5):504–14. Epub 2018/07/27. pmid:30048232.

55. Conde-Saldana CC, Cunha MS, Albornoz-Garzon JG, Barreto CAV, Ibagon N, Villa-Navarro FA, et al. Karyotypic Divergence of Two Co-Occurring Species of Andean Climbing Catfishes (Siluriformes: Loricarioidei: Astroblepidae). Zebrafish. 2019;16(1):106–14. Epub 2018/11/21. pmid:30457940.

56. Blanco DR, Vicari MR, Lui RL, Artoni RF, de Almeida MC, Traldi JB, et al. Origin of the X1X1X2X2/X1X2Y sex chromosome system of Harttia punctata (Siluriformes, Loricariidae) inferred from chromosome painting and FISH with ribosomal DNA markers. Genetica. 2014;142(2):119–26. Epub 2014/03/01. pmid:24577679.

57. Sassi FMC, Deon GA, Moreira-Filho O, Vicari MR, Bertollo LAC, Liehr T, et al. Multiple Sex Chromosomes and Evolutionary Relationships in Amazonian Catfishes: The Outstanding Model of the Genus Harttia (Siluriformes: Loricariidae). Genes (Basel). 2020;11(10). Epub 2020/10/15. pmid:33050411; PubMed Central PMCID: PMC7600804.

58. Ferreira AMV, Viana PF, Zuanon J, Ezaz T, Cioffi MB, Takagui FH, et al. Cytogenetic Analysis of Panaqolus tankei Cramer & Sousa, 2016 (Siluriformes, Loricariidae), an Ornamental Fish Endemic to Xingu River, Brazil. Cytogenet Genome Res. 2021;161(3–4):187–94. Epub 2021/03/22. pmid:33744896.

59. Vanzela ALL, Swarca AC, Dias AL, Stolf R, Ruas PM, Ruas CF, et al. Differential Distribution of (GA)9+C Microsatellite on Chromosomes of Some Animal and Plant Species. Cytologia. 2002;67. https://10.1508/cytologia.67.9.

60. Gouveia JG, Wolf IR, Vilas-Boas LA, Heslop-Harrison JS, Schwarzacher T, Dias AL. Repetitive DNA in the Catfish Genome: rDNA, Microsatellites, and Tc1-Mariner Transposon Sequences in Imparfinis Species (Siluriformes, Heptapteridae). J Hered. 2017;108(6):650–7. Epub 2017/08/20. pmid:28821184.

61. Conde-Saldana CC, Barreto CAV, Villa-Navarro FA, Dergam JA. An Unusual Accumulation of Ribosomal Multigene Families and Microsatellite DNAs in the XX/XY Sex Chromosome System in the Trans-Andean Catfish Pimelodella cf. chagresi (Siluriformes:Heptapteridae). Zebrafish. 2018;15(1):55–62. Epub 2017/11/02. pmid:29090985.

62. Malimpensa GC, Traldi JB, Toyama D, Henrique-Silva F, Vicari MR, Moreira-Filho O. Chromosomal Mapping of Repeat DNA in Bergiaria westermanni (Pimelodidae, Siluriformes): Localization of 45S rDNA in B Chromosomes. Cytogenet Genome Res. 2018;154(2):99–106. Epub 2018/04/11. pmid:29635248.

63. Utsunomia R, Melo S, Scacchetti PC, Oliveira C, Machado MA, Pieczarka JC, et al. Particular Chromosomal Distribution of Microsatellites in Five Species of the Genus Gymnotus (Teleostei, Gymnotiformes). Zebrafish. 2018;15(4):398–403. Epub 2018/06/22. pmid:29927722.

64. Lima-Filho PA, Amorim KD, Cioffi MB, Bertollo LA, Molina WF. Chromosomal mapping of repetitive DNAs in Gobionellus oceanicus and G. stomatus (Gobiidae; Perciformes): A shared XX/XY system and an unusual distribution of 5S rDNA sites on the Y chromosome. Cytogenet Genome Res. 2014;144(4):333–40. Epub 2015/02/28. pmid:25720317.

65. Poltronieri J, Marquioni V, Bertollo LA, Kejnovsky E, Molina WF, Liehr T, et al. Comparative chromosomal mapping of microsatellites in Leporinus species (Characiformes, Anostomidae): unequal accumulation on the W chromosomes. Cytogenet Genome Res. 2014;142(1):40–5. Epub 2013/11/13. pmid:24217024.

66. de Moraes RLR, Sember A, Bertollo LAC, de Oliveira EA, Rab P, Hatanaka T, et al. Comparative Cytogenetics and Neo-Y Formation in Small-Sized Fish Species of the Genus Pyrrhulina (Characiformes, Lebiasinidae). Front Genet. 2019;10:678. Epub 2019/08/21. pmid:31428127; PubMed Central PMCID: PMC6689988.

67. de Moraes RLR, Bertollo LAC, Marinho MMF, Yano CF, Hatanaka T, Barby FF, et al. Evolutionary Relationships and Cytotaxonomy Considerations in the Genus Pyrrhulina (Characiformes, Lebiasinidae). Zebrafish. 2017;14(6):536–46. Epub 2017/08/03. pmid:28767325.

68. Sassi FMC, Oliveira EA, Bertollo LAC, Nirchio M, Hatanaka T, Marinho MMF, et al. Chromosomal Evolution and Evolutionary Relationships of Lebiasina Species (Characiformes, Lebiasinidae). Int J Mol Sci. 2019;20(12). Epub 2019/06/19. pmid:31208145; PubMed Central PMCID: PMC6628269.

69. Yano CF, Bertollo LA, Molina WF, Liehr T, Cioffi Mde B. Genomic organization of repetitive DNAs and its implications for male karyotype and the neo-Y chromosome differentiation in Erythrinus erythrinus (Characiformes, Erythrinidae). Comp Cytogenet. 2014;8(2):139–51. Epub 2014/08/26. pmid:25147625; PubMed Central PMCID: PMC4137284.

70. Terencio ML, Schneider CH, Gross MC, Vicari MR, Farias IP, Passos KB, et al. Evolutionary dynamics of repetitive DNA in semaprochilodus (characiformes, prochilodontidae): a fish model for sex chromosome differentiation. Sex Dev. 2013;7(6):325–33. Epub 2013/12/04. pmid:24296872.

71. de Oliveira EA, Bertollo LA, Yano CF, Liehr T, Cioffi Mde B. Comparative cytogenetics in the genus Hoplias (Characiformes, Erythrinidae) highlights contrasting karyotype evolution among congeneric species. Mol Cytogenet. 2015;8:56. Epub 2015/08/01. pmid:26225139; PubMed Central PMCID: PMC4518567.

72. Cioffi MB, Kejnovsky E, Bertollo LA. The chromosomal distribution of microsatellite repeats in the genome of the wolf fish Hoplias malabaricus, focusing on the sex chromosomes. Cytogenet Genome Res. 2011;132(4):289–96. Epub 2010/11/26. pmid:21099206.

73. de Oliveira EA, Sember A, Bertollo LAC, Yano CF, Ezaz T, Moreira-Filho O, et al. Tracking the evolutionary pathway of sex chromosomes among fishes: characterizing the unique XX/XY1Y2 system in Hoplias malabaricus (Teleostei, Characiformes). Chromosoma. 2018;127(1):115–28. Epub 2017/11/11. pmid:29124392.

74. de MCSF Perez MF, Oliveira VCS Deon GA, de Souza FHS Ferreira PHN, et al. High Genetic Diversity despite Conserved Karyotype Organization in the Giant Trahiras from Genus Hoplias (Characiformes, Erythrinidae). Genes (Basel). 2021;12(2). Epub 2021/02/14. pmid:33578790; PubMed Central PMCID: PMC7916553.

75. Ziemniczak K, Traldi JB, Nogaroto V, Almeida MC, Artoni RF, Moreira-Filho O, et al. In situ localization of (GATA)n and (TTAGGG)n repeated DNAs and W sex chromosome differentiation in Parodontidae (Actinopterygii: Characiformes). Cytogenet Genome Res. 2014;144(4):325–32. Epub 2015/02/11. pmid:25662193.

76. Pucci MB, Barbosa P, Nogaroto V, Almeida MC, Artoni RF, Scacchetti PC, et al. Chromosomal Spreading of Microsatellites and (TTAGGG)n Sequences in the Characidium zebra and C. gomesi Genomes (Characiformes: Crenuchidae). Cytogenet Genome Res. 2016;149(3):182–90. Epub 2016/11/03. pmid:27504623.

77. Scacchetti PC, Utsunomia R, Pansonato-Alves JC, Vicari MR, Artoni RF, Oliveira C, et al. Chromosomal Mapping of Repetitive DNAs in Characidium (Teleostei, Characiformes): Genomic Organization and Diversification of ZW Sex Chromosomes. Cytogenet Genome Res. 2015;146(2):136–43. Epub 2015/08/19. pmid:26277929.

78. Carvalho PC, de Oliveira EA, Bertollo LAC, Yano CF, Oliveira C, Decru E, et al. First Chromosomal Analysis in Hepsetidae (Actinopterygii, Characiformes): Insights into Relationship between African and Neotropical Fish Groups. Front Genet. 2017;8:203. Epub 2018/01/10. pmid:29312435; PubMed Central PMCID: PMC5733008.

79. Martins NF, Cioffi MB, Troy WP, Martinez JF, Moreira-Filho O, Bertollo LA. Differentiation and evolutionary relationships in Erythrinus erythrinus (Characiformes, Erythrinidae): occurrence and distribution of B chromosomes. Genet Mol Res. 2014;13(3):7094–101. Epub 2014/03/13. pmid:24615114.

80. Salgado FS, Cunha MS, Melo S, Dergam JA. Cytogenetic analysis of Hypomasticus copelandii and H. steindachneri: relevance of cytotaxonomic markers in the Anostomidae family (Characiformes). Comp Cytogenet. 2021;15(1):65–76. Epub 2021/03/30. pmid:33777329; PubMed Central PMCID: PMC7969579.

81. de ASDM Daniel SN, Camacho JP, Utsunomia R, Ruiz-Ruano FJ, Penitente M, et al. Origin of B chromosomes in the genus Astyanax (Characiformes, Characidae) and the limits of chromosome painting. Mol Genet Genomics. 2016;291(3):1407–18. Epub 2016/03/18. pmid:26984341.

82. Yano CF, Poltronieri J, Bertollo LA, Artoni RF, Liehr T, de Bello Cioffi M. Chromosomal mapping of repetitive DNAs in Triportheus trifurcatus (Characidae, Characiformes): insights into the differentiation of the Z and W chromosomes. PLoS One. 2014;9(3):e90946. Epub 2014/03/19. pmid:24632562; PubMed Central PMCID: PMC3954618.

83. Piscor D, Paiz LM, Baumgartner L, Cerqueira FJ, Fernandes CA, Lui RL, et al. Chromosomal mapping of repetitive sequences in Hyphessobrycon eques (Characiformes, Characidae): a special case of the spreading of 5S rDNA clusters in a genome. Genetica. 2020;148(1):25–32. Epub 2020/01/31. pmid:31997050.

84. Santos EOD, Deon GA, Almeida RB, Oliveira EA, Nogaroto V, Silva HPD, et al. Cytogenetics and DNA barcode reveal an undescribed Apareiodon species (Characiformes: Parodontidae). Genet Mol Biol. 2019;42(2):365–73. Epub 2019/07/02. pmid:31259363; PubMed Central PMCID: PMC6726146.

85. Traldi JB, Vicari MR, Martinez JF, Blanco DR, Lui RL, Azambuja M, et al. Recent Apareiodon species evolutionary divergence (Characiformes: Parodontidae) evidenced by chromosomal and molecular inference. Zoologischer Anzeiger. 2020:166–76. https://doi.org/10.1016/j.jcz.2020.10.010.

86. Eschmeyer’s Catalog of Fishes [Internet]. 2022.

87. Akama A. Sistemática dos gêneros Parauchenipterus Bleeker, 1862 e Trachelyopterus Valenciennes, 1840 (Siluriformes, Auchenipteridae): Universidade de São Paulo; 2004.

88. Lui RL, Blanco DR, Margarido VP, Moreira-Filho O. Chromosome characterization and biogeographic relations among three populations of the driftwood catfish Parauchenipterus galeatus (Linnaeus, 1766) (Siluriformes: Auchenipteridae) in Brazil. Biological Journal of the Linnean Society. 2010;99:648–56. https://doi.org/10.1111/j.1095-8312.2009.01389.x.

89. Santos DP, Felicetti D, Baumgartner L, Margarido VP, Blanco DR, Moreira-Filho O, et al. Contributions to the taxonomy of Trachelyopterus (Siluriformes): comparative cytogenetic analysis in three species of Auchenipteridae. Neotropical Ichthyology. 2021;19:e200115. https://doi.org/10.1590/1982-0224-2020-0115.

90. Haerter CAG, Margarido VP, Blanco DR, Traldi JB, Feldberg E, Lui RL. Contributions to Trachelyopterus (Siluriformes: Auchenipteridae) species diagnosis by cytotaxonomic autapomorphies: from U2 snRNA chromosome polymorphism to rDNA and histone gene synteny. Organisms Diversity & Evolution. 2022. https://doi.org/10.1007/s13127-022-00560-0.

91. Felicetti D. Análise genética de populações simpátricas de Trachelyopterus valenciennes, 1840 (auchenipteridae) das bacias Amazônica e Paraguai: Universidade Estadual do Oeste do Paraná; 2018.

92. White T, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis M, Gelfand D, Sninsky J, editors. PCR Protocols: A Guide to Methods and Applications: Academic Press; 1990. p. 315–22.

93. Martins C, Galetti PM Jr., Organization of 5S rDNA in species of the fish Leporinus: two different genomic locations are characterized by distinct nontranscribed spacers. Genome. 2001;44(5):903–10. Epub 2001/10/30. pmid:11681615.

94. Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas JM, et al. Histone H3 and U2 snRNA DNA sequences and arthropod evolution. Australian Journal of Zoology 1998;46:419–37.

95. Pineau P, Henry M, Suspene R, Marchio A, Dettai A, Debruyne R, et al. A universal primer set for PCR amplification of nuclear histone H4 genes from all animal species. Mol Biol Evol. 2005;22(3):582–8. Epub 2004/12/17. pmid:15601889.

96. Bueno D, Palacios-Gimenez OM, Cabral-de-Mello DC. Chromosomal Mapping of Repetitive DNAs in the Grasshopper Abracris flavolineata Reveal Possible Ancestry of the B Chromosome and H3 Histone Spreading. PLoS One. 2013;8(6):e66532. Epub 2013/07/05. pmid:23826099; PubMed Central PMCID: PMC3694960.

97. Ijdo JW, Wells RA, Baldini A, Reeders ST. Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCR. Nucleic Acids Res. 1991;19(17):4780. Epub 1991/09/11. pmid:1891373; PubMed Central PMCID: PMC328734.

98. Hall TA. BioEdit: A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999;41:95–8.

99. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997;25:3389–402. pmid:9254694

100. Griffiths SP. The use of clove oil as an anaesthetic and method or sampling intertidal rockpool fishes. Journal of Fish Biology. 2005;57:1453–64. https://doi.org/10.1111/j.1095-8649.2000.tb02224.x.

101. Bertollo LAC, Takahashi CS, Moreira-Filho O. Cytotaxonomic considerations on Hoplias lacerdae (Pisces, Erythrinidae). Revista Brasileira de Genética. 1978;1:103–20.

102. Levan A, Fredga K, Sandberg AA. Nomenclature for centromeric position on chromosomes. Hereditas. 1964;52:201–20. https://doi.org/10.1111/j.1601-5223.1964.tb01953.x.

103. Pinkel D, Straume T, Gray JW. Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization. Proc Natl Acad Sci U S A. 1986;83(9):2934–8. Epub 1986/05/01. pmid:3458254; PubMed Central PMCID: PMC323421.

104. Margarido VP, Moreira-Filho O. Karyotypic differentiation through chromosome fusion and number reduction in Imparfinis hollandi Genetics and Molecular Biology. 2008;31:235–8. https://doi.org/10.1590/S1415-47572008000200012.

105. Santos FAB, Lemes RB, Otto PA. HW_TEST, a program for comprehensive HARDY-WEINBERG equilibrium testing. Genet Mol Biol. 2020;43(2):e20190380. Epub 2020/05/19. pmid:32422648; PubMed Central PMCID: PMC7323896.

106. Morescalchi MA, Stingo V, Capriglione T. Cytogenetic analysis in Polypterus ornatipinnis (Actinopterygii, Cladistia, Polypteridae) and 5S rDNA. Mar Genomics. 2011;4(1):25–31. Epub 2011/03/25. pmid:21429462.

107. Chang SB, Yang TJ, Datema E, van Vugt J, Vosman B, Kuipers A, et al. FISH mapping and molecular organization of the major repetitive sequences of tomato. Chromosome Res. 2008;16(7):919–33. Epub 2008/08/09. pmid:18688733.

108. Kubat Z, Hobza R, Vyskot B, Kejnovsky E. Microsatellite accumulation on the Y chromosome in Silene latifolia. Genome. 2008;51(5):350–6. Epub 2008/04/29. pmid:18438438.

109. Micolino R, Cristiano MP, Cardoso DC. Population-Based Cytogenetic Banding Analysis and Phylogenetic Relationships of the Neotropical Fungus-Farming Ant Trachymyrmex holmgreni Wheeler, 1925. Cytogenetic and Genome Research. 2019;159:151–61. pmid:31683269

110. Epplen JT, McCarrey JR, Sutou S, Ohno S. Base sequence of a cloned snake W-chromosome DNA fragment and identification of a male-specific putative mRNA in the mouse. Proc Natl Acad Sci U S A. 1982;79(12):3798–802. Epub 1982/06/01. pmid:6954524; PubMed Central PMCID: PMC346515.

111. Epplen JT. On simple repeated GATCA sequences in animal genomes: a critical reappraisal. J Hered. 1988;79(6):409–17. Epub 1988/11/01. pmid:3062083.

112. Mazzoleni S, Augstenova B, Clemente L, Auer M, Fritz U, Praschag P, et al. Sex is determined by XX/XY sex chromosomes in Australasian side-necked turtles (Testudines: Chelidae). Sci Rep. 2020;10(1):4276. Epub 2020/03/11. pmid:32152354; PubMed Central PMCID: PMC7062838.

113. Metzgar D, Bytof J, Wills C. Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000;10(1):72–80. Epub 2000/01/25. pmid:10645952; PubMed Central PMCID: PMC310501.

114. Oliveira EJ, Pádua JG, Zucchi MI, Vencovsky R, Vieira MLC. Origin, evolution and genome distribution of microsatellites. Genetics and Molecular Biology. 2006;29:294–307. https://doi.org/10.1590/S1415-47572006000200018.

115. Schwarzacher T, Heslop-Harrison JS. Practical in Situ Hybridization. BIOS Scientific Publishers2000.

116. Richard GF, Kerrest A, Dujon B. Comparative genomics and molecular dynamics of DNA repeats in eukaryotes. Microbiol Mol Biol Rev. 2008;72(4):686–727. Epub 2008/12/05. pmid:19052325; PubMed Central PMCID: PMC2593564.

117. Abbott R, Albach D, Ansell S, Arntzen JW, Baird SJ, Bierne N, et al. Hybridization and speciation. J Evol Biol. 2013;26(2):229–46. Epub 2013/01/18. pmid:23323997.

118. Sedghifar A, Brandvain Y, Ralph P, Coop G. The Spatial Mixing of Genomes in Secondary Contact Zones. Genetics. 2015;201(1):243–61. Epub 2015/07/25. pmid:26205988; PubMed Central PMCID: PMC4566267.

119. Saadi A. Neotectônica da plataforma Brasileira: Esboço e Interpretação preliminares. Geonomos. 1993;1:1–15.

120. Saadi A, Bezerra FHR, Costa RD, Igreja HLS, Franzinelli E. Neotectonica da Plataforma Brasileira. In: Souza CLG, Suguio K, Oliveira AMS, Oliveira PE, editors. Quaternário do Brasil: Editora Holos; 2005. p. 211–34.

121. Lima FTC, Ribeiro AC. Continental-scale tectonic controls of biogeography and ecology. In: Albert JS, Reis RE, editors. Historical Biogeography of Neotropical Freshwater Fishes: University of California Press; 2011.

122. Hoffmann AA, Rieseberg LH. Revisiting the Impact of Inversions in Evolution: From Population Genetic Markers to Drivers of Adaptive Shifts and Speciation? Annu Rev Ecol Evol Syst. 2008;39:21–42. Epub 2008/12/01. pmid:20419035; PubMed Central PMCID: PMC2858385.

123. Aradottir GI, Angus RB. A chromosomal analysis of some water beetle species recently transferred from Agabus Leach to Ilybius Erichson, with particular reference to the variation in chromosome number shown by I. montanus Stephens (Coleoptera: Dytiscidae). Hereditas. 2004;140(3):185–92. Epub 2004/06/17. pmid:15198708; PubMed Central PMCID: PMC7202062.

124. Kuvangkadilok C, Boonkemtong C, Phayuhasena S, Baimai V. Larval polytene chromosomes of black flies (Simulium) from Thailand. I. Comparison among five species in the subgenus Gomphostilbia enderlein. Genetica. 2003;118(1):69–81. Epub 2003/05/08. pmid:12733666.

125. White S. Modes of Speciation. Paleobiology. 1978:373–97.

126. Searle JB. Chromosomal hybrid zones in eutherian mammals. In: Harrison RG, editor. Hybrid zones and the evolutionary process: Oxford University Press; 1993. p. 309–53.

127. Stebbins GL Jr. Variation and evolution in plants: Columbia University Press; 1950.

128. Levin DA. The long wait for hybrid sterility in flowering plants. New Phytol. 2012;196(3):666–70. Epub 2012/09/13. pmid:22966819.