
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39223262
70920
10.1038/s41598-024-70920-7
Article
Satellite DNAs and the evolution of the multiple X1X2Y sex chromosomes in the wolf fish Hoplias malabaricus (Teleostei; Characiformes)
http://orcid.org/0000-0002-8084-1787
Toma Gustavo Akira 1
http://orcid.org/0000-0003-4441-9615
Sember Alexandr 2
http://orcid.org/0000-0002-7595-0014
Goes Caio Augusto Gomes 3
http://orcid.org/0000-0002-6856-2152
Kretschmer Rafael 4
http://orcid.org/0000-0001-8845-3845
Porto-Foresti Fabio 3
http://orcid.org/0000-0001-6487-2834
Bertollo Luiz Antônio Carlos 1
http://orcid.org/0000-0003-1672-3054
Liehr Thomas thomas.liehr@med.uni-jena.de

5
http://orcid.org/0000-0003-4565-6065
Utsunomia Ricardo 3
http://orcid.org/0000-0003-4340-1464
de Bello Cioffi Marcelo 1
1 https://ror.org/00qdc6m37 grid.411247.5 0000 0001 2163 588X Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, SP 13565-905 Brazil
2 https://ror.org/053avzc18 grid.418095.1 0000 0001 1015 3316 Laboratory of Fish Genetics, Institute of Animal Physiology and Genetics, Czech Academy of Sciences, 277 21 Liběchov, Czech Republic
3 https://ror.org/00987cb86 grid.410543.7 0000 0001 2188 478X Faculdade de Ciências, UNESP, Bauru, SP 17033-360 Brazil
4 https://ror.org/05msy9z54 grid.411221.5 0000 0001 2134 6519 Departamento de Ecologia, Zoologia e Genética, Instituto de Biologia, Universidade Federal de Pelotas, Pelotas, RS 96010-610 Brazil
5 grid.9613.d 0000 0001 1939 2794 Jena University Hospital, Institute of Human Genetics, Friedrich Schiller University, 07747 Jena, Germany
2 9 2024
2 9 2024
2024
14 204025 2 2024
22 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Multiple sex chromosomes usually arise from chromosomal rearrangements which involve ancestral sex chromosomes. There is a fundamental condition to be met for their long-term fixation: the meiosis must function, leading to the stability of the emerged system, mainly concerning the segregation of the sex multivalent. Here, we sought to analyze the degree of differentiation and meiotic pairing properties in the selected fish multiple sex chromosome system present in the wolf-fish Hoplias malabaricus (HMA). This species complex encompasses seven known karyotype forms (karyomorphs) where the karyomorph C (HMA-C) exhibits a nascent XY sex chromosomes from which the multiple X1X2Y system evolved in karyomorph HMA-D via a Y-autosome fusion. We combined genomic and cytogenetic approaches to analyze the satellite DNA (satDNA) content in the genome of HMA-D karyomorph and to investigate its potential contribution to X1X2Y sex chromosome differentiation. We revealed 56 satDNA monomers of which the majority was AT-rich and with repeat units longer than 100 bp. Seven out of 18 satDNA families chosen for chromosomal mapping by fluorescence in situ hybridization (FISH) formed detectable accumulation in at least one of the three sex chromosomes (X1, X2 and neo-Y). Nine satDNA monomers showed only two hybridization signals limited to HMA-D autosomes, and the two remaining ones provided no visible FISH signals. Out of seven satDNAs located on the HMA-D sex chromosomes, five mapped also to XY chromosomes of HMA-C. We showed that after the autosome-Y fusion event, the neo-Y chromosome has not substantially accumulated or eliminated satDNA sequences except for minor changes in the centromere-proximal region. Finally, based on the obtained FISHpatterns, we speculate on the possible contribution of satDNA to sex trivalent pairing and segregation.

Keywords

Meiosis
Sex trivalent
FISH
Multiple sex chromosomes
Satellitome
Subject terms

Genetics
Cytogenetics
Evolutionary biology
Genomics
http://dx.doi.org/10.13039/501100001807 Fundação de Amparo à Pesquisa do Estado de São Paulo 2022/00427-3 2023/00955-2 Toma Gustavo Akira de Bello Cioffi Marcelo http://dx.doi.org/10.13039/501100003593 Conselho Nacional de Desenvolvimento Científico e Tecnológico 302928/2021-9 de Bello Cioffi Marcelo Friedrich-Schiller-Universität Jena (1010)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Multiple sex chromosomes arise, in most cases, from chromosomal rearrangements between autosomes and existing sex chromosomes1–3. This way, a new genomic portion becomes sex-linked and may both restart the sex chromosome differentiation process4,5 and provide genetic material for sex chromosome turnover and/or hybrid incompatibility between the populations differing by their sex chromosome systems6–9. The analysis of multiple sex chromosomes may hence provide important insights into early phases of restarted sex chromosome differentiation, spreading over the newly added autosomal part, and shed light on the drivers of chromosome fusion1,3,10–12.

The most prevalent configuration of multiple sex chromosomes recorded thus far in nature is ♀X1X1X2X2/♂X1X2Y created by centric or tandem fusion between a Y chromosome and an autosome1,2,13,14. The system is clearly distinguishable both by a consistent difference between male and female diploid chromosome number (2n), and a remarkably large biarmed (meta- or submetacentric) neo-Y chromosome present in one copy in males and entirely missing in females. While neo-Y hence achieves instant heteromorphy (i.e. it can be distinguished cytologically), the process is not analogous to heteromorphy achieved by a gradual sequence differentiation and gene content degradation15. Hence, in multiple sex chromosome systems, the heteromorphy does not necessarily correlate with the degree of sex chromosome differentiation3. To date, out of the 81 known cases of fish multiple sex chromosomes (2,16–19 and references therein), 65 are of the X1X2Y type and the overwhelming majority of them exhibits signs of low differentiation, inferred mostly from the degree of repetitive DNA and heterochromatin accumulation2,17,18.

Y-autosome fusion results in the formation of a sex trivalent in the male first meiotic division. Its correct pairing and segregation are critical for stability and fixation as a multiple sex chromosome system given that deviations from these processes may lead to meiotic arrest or reduced fertility and/or viability of the progeny3,5,20,21. Although the properties of sex trivalents have been investigated in diverse eukaryotic lineages22–27 (including fishes28,29), little is known about the mechanisms related to their correct alignment and segregation, and the role of repetitive DNAs in these processes.

Repetitive DNA, i.e. the highly variable fraction of the genome present in multiple dispersed or tandemly arrayed copies, has a multilayered influence on the evolution of sex chromosomes30–33. Firstly, repetitive DNAs usually populate the region of abolished recombination on sex chromosomes, thus contributing to gene content decay and further spread of the region of differentiation on the Y (or, by analogy, W) chromosome by repeat amplification and by triggering heterochromatinization30–34. Secondly, repetitive DNA may induce intrachromosomal rearrangements via non-allelic homologous recombination which may lead, e.g., to inversion or deletion and this way to a further spread of non-recombining region and/or—similarly to repeat amplification—potentially to Y/W sex chromosome heteromorphy4,35. Thirdly, repeat content may also promote sex chromosome-autosome fusions36,37 or shield the newly added autosomal part from the effect of dosage compensation if given enough time23,38,39. Last but not least, certain repeats have been found to directly rewire the expression patterns within the sex determination regulatory network40,41 and to facilitate sex chromosome pairing in meiosis21,42.

Satellite DNA (satDNA) is a highly abundant repetitive DNA class with an extreme diversity in monomer size and nucleotide sequence as well as distribution patterns (mostly tandemly organized but also dispersed, present both in heterochromatin and euchromatin)43–45. As these repeats often occupy key chromosomal regions such as centromeres and telomeres, satDNA sequences play important roles in chromosome architecture and function, including centromere organization and kinetochore assembly, which allow proper chromatid cohesion and segregation in cell division46–48. SatDNAs are among the fastest evolving sequences in the genome and therefore they generate highly species-specific landscapes in terms of quantity, distribution, and diversity of satDNA monomer collection (the so-called satellitome)43,44,49,50. During the past decade, the advancements in bioinformatic tools allowed a steep increase in studies aimed at satellitome characterization and its subsequent mapping on chromosomes by fluorescence in situ hybridization (FISH)43,49,51. Given the above-mentioned satDNA properties, this type of investigation has already provided important novel insights into satDNA evolution at inter-specific and inter-population levels52–55, chromosomal markers for genomic composition of hybrids and polyploids56,57 (peri)centromeric repeat turnover58–60 and the evolutionary dynamics of sex61,62 and B chromosomes63, with all of these aspects already being studied also in teleost fishes64–70.

In the present work, we aimed to analyze the differentiation and meiotic pairing properties of a selected X1X2Y multiple sex chromosome system through the lens of satellitome analysis. For this, we chose the wolf-fish Hoplias malabaricus (Characiformes, Erythrinidae) as a model. This species complex (hereafter abbreviated as HMA) comprises seven karyotype forms (karyomorphs) denoted A–G, which represent populations differing by their karyotype features, including the presence of four different karyomorph-specific sex chromosome systems (based on their configuration and degree of differentiation): (1) a well-differentiated, heteromorphic XX/XY (karyomorph HMA-B), non-homologous to any of the following sex chromosome systems), (2) nascent XX/XY systems (HMA-C and HMA-F; these are mutually non-homeologous), and two multiple sex chromosome systems: (3) X1X1X2X2/X1X2Y in HMA-D (derived from XX/XY of HMA-C), and (4) XX/XY1Y2 in HMA-G) (derived from XX/XY of HMA-F)71–74.

Molecular cytogenetic analyses showed that the X and Y chromosomes of HMA-C differ just slightly by degree of their (peri)centromeric repeats amplification, and that the X1X1Y system of HMA-D evolved from HMA-C XY system through a tandem (end-to-end) fusion between the short (p) arms of a Y chromosome and a p-arms of one homolog from autosome pair 2074. Meiotic analysis indicated that the neo-Y, the X1 and X2 chromosomes pair regularly by forming a trivalent without obvious irregularities in the first meiotic division28. This was further corroborated by repetitive DNA analysis74, whole chromosome painting (WCP)72 and comparative genomic hybridization (CGH)75. However, a detailed study on the sequence differentiation on these multiple sex chromosomes is still missing.

We herein analyzed the satellitome of HMA-D karyomorph by a combination of genomic and cytogenetic approaches. We focused specifically on the satDNA landscape on X1X2Y multiple sex chromosomes to trace the origin and pathways of differentiation of this system. To determine whether there was a significant change in the Y-linked satDNA patterns after the Y-autosome fusion, we also investigated the chromosomal distribution of selected satDNA monomers in the homeologous XX/XY sex chromosome system of HMA-C. The results showed that the neo-Y chromosome has not yet undergone major changes in its satDNA landscape after the fusion event which corroborates the previous assumptions that this multiple sex chromosomes is still poorly differentiated.

Results

Satellitome features of HMA-D

Based on the five subsequent rounds of analysis via Tandem Repeat Analyzer (TAREAN) pipeline51 we revealed 56 satDNA monomers in HMA-D (Supplementary Table 1). Among them, HmaSat03 was the only one characterized previously76, corresponding to the 5SHindIII-DNA sequence. Forty HmaSatDNA monomers had the proportion of AT-pairs over 50%. The repeat unit lengths (RUL) varied from 6 to 2774 bp (average of 493 bp), with 39 HmaSatDNA monomers being shorter, and 17 ones longer than 100 bp (Supplementary Table 1). By calculating the sex abundance ratio (M/F) (see Material and Methods section for explanation) we found seven male (HmaSat07; HmaSat14; HmaSat18; HmaSat32; HmaSat34; HmaSat42; HmaSat56) and six female-biased HmaSatDNAs (HmaSat11; HmaSat17; HmaSat25; HmaSat36; HmaSat48; HmaSat51) (Supplementary Table 1). The homology search between HmaSatDNAs revealed the presence of five superfamilies (SF1-SF5) (Supplementary Table 1).

Hybridization patterns of selected HmaSatDNA monomers in the HMA-D karyomorph

Among the 18 mapped HmaSatDNA monomers (marked by asterisk in Supplementary Table 1), only two (HmaSat14 and HmaSat56) did not produce specific signals on male and female mitotic chromosomes after FISH (data not shown). This outcome may be attributed to their potentially low abundance in the genome or their arrangement into clusters that are too small for effective cytological detection. The remaining examined satDNA monomers provided varied hybridization patterns (Figs. 1 and 2 and Supplementary Figs. 1 and 2). Only two hybridization signals, restricted to a certain autosome pair, were recorded for nine HmaSatDNA monomers, namely HmaSat11, HmaSat17, HmaSat25, HmaSat32, HmaSat34, HmaSat36, HmaSat42, HmaSat48, and HmaSat51, with HmaSat25 and HmaSat32 showing adjacent co-localized arrangement (Fig. 1 and Supplementary Figs. 1 and 2). On the other hand, FISH with HmaSat07 probe yielded strong (peri)centromeric signals exclusively on sex chromosomes, extended towards the p- and q-arms of the X1 and neo-Y, respectively (Fig. 2 and Supplementary Figs. 1 and 2). HmaSat07 was hence the only sex chromosome-specific satDNA monomer revealed in the present study. Finally, the remaining six satDNA monomers displayed FISH signals on multiple autosomes, neo-Y, and either X1 or X2 sex chromosomes. HmaSat01 and HmaSat18 exhibited highly similar patterns, with strong clusters being detected in the telomeric regions of many autosomes along with X2 and neo-Y chromosomes. On the other hand, HmaSat03, HmaSat04 and HmaSat05 shared the pattern of (peri)centromeric signals, being present on many (HmaSat04) or almost all (HmaSat03 and HmaSat05) chromosomes of the complement, including neo-Y and either X2 (HmaSat03) or X1 (HmaSat04 and HmaSat05) chromosome. A specific pattern was produced by the probe for HmaSat02, with prominent clusters at many terminal and interstitial regions across the entire chromosome set (Fig. 2).Fig. 1 Male mitotic metaphases of Hoplias malabaricus, karyomorph D (HMA-D) after FISH with a set of HmaSatDNA monomers located exclusively on autosomes. SatDNA family names are indicated on the right bottom, in green (Atto-488-dUTP) or red (Atto-550-dUTP). The arrows indicate the X1, X2, and neo-Y sex chromosomes. Bar = 5 μm.

Fig. 2 Male mitotic metaphases of Hoplias malabaricus, karyomorph D (HMA-D) after FISH with a set of HmaSatDNA monomers located both on autosomes and sex chromosomes. SatDNA family names are indicated on the right bottom, in green (Atto-488-dUTP) or red (Atto-550-dUTP). The arrows indicate the X1, X2 and neo-Y sex chromosomes. Bar = 5 μm.

Patterns of HmaSatDNA monomers on the meiotic sex trivalent formed by X1X2Y multiple sex chromosomes of HMA-D

To better understand the organization of HmaSatDNA monomers on the X1X2Y multiple sex chromosomes of HMA-D during their meiotic pairing into the trivalent structure, we sought to map the seven HmaSatDNAs that display detectable clusters on these sex chromosomes (Fig. 2). A prominent cluster was found for the HmaSat01 and HmaSat18 monomers, corresponding to the segment formed by the terminal regions of the X2 and neo-Y chromosomes (Fig. 3). Four sites that resided on the X1, X2, and neo-Y chromosomes were displayed by the HmaSat02-specific probe. In contrast, the HmaSat03 probe displayed two signals, one of which populated the q-arms of the neo-Y chromosome and the second located in the centromeric region of the X2 chromosome (Fig. 3). All the three remaining HmaSatDNA monomers (HmaSat04, HmaSat05, and HmaSat07) presented two FISH signals in the (peri)centromeric region of the X1 and neo-Y chromosomes (Fig. 3).Fig. 3 Meiotic chromosomes of male Hoplias malabaricus karyomorph D (HMA-D) after FISH with a set of HmaSatDNA monomers located both on autosomes and sex chromosomes. Arrows indicate the sex trivalent. Bar = 5 μm.

Chromosomal location of HmaSatDNA monomers in HMA-C

We mapped the subset of seven HmaSatDNA monomers, which hybridized to HMA-D sex chromosomes, also to metaphases of H. malabaricus males of karyomorph C (HMA-C) to examine their possible presence on nascent XY sex chromosomes (Fig. 4 and Supplementary Fig. 3). Most of FISH signals resembled the patterns revealed on the HMA-D sex chromosomes, with the exception of HmaSat01 (Fig. 2). The same satDNA monomers found on the X1 and the (peri)centromeric region of the neo-Y, (HmaSat02; HmaSat04; HmaSat05; HmaSat07) were located in the (peri)centromeric regions of the corresponding X and Y chromosomes in the karyomorph HMA-C (Fig. 4 and Supplementary Fig. 3). In a similar way, the HmaSat01 and HmaSat02 monomers, with signals on the X2 and neo-Y (HmaSat01), or all multiple sex chromosome (HmaSat02) of the HMA-D, also formed telomeric signals in both nascent sex chromosomes of the HMA-C (Fig. 4 and Supplementary Fig. 3). In turn, the HmaSat03 and HmaSat18 monomers were the only ones without detectable hybridization signals on the X and Y chromosomes of HMA-C. In the case of HmaSat03, the situation particularly reflects the fact that clusters of this satDNA reside on the X2 sex chromosome of HMA-D, i.e. the newly added chromosomal part (Fig. 2), therefore this satDNA joined the sex chromosome system only after the Y-autosome fusion. On the other hand, HmaSat18 was entirely absent from the HMA-C chromosomes, at least under the FISH resolution (Fig. 4 and Supplementary Fig. 3). Lastly, it is worth noting, that the probe corresponding to HmaSat05 monomer marked the centromeric regions of all but one chromosome pair in the HMA-C complement, which, after comparison with the HMA-D pattern (particularly the lack of signal on X2), indicates that the sole chromosome pair without the signal represents highly probably the linkage group (pair 20) which has been incorporated into the X1X2Y system in karyomorph HMA-D.Fig. 4 Male mitotic metaphases of Hoplias malabaricus, karyomorph C (HMA-C) after FISH with a set of HmaSatDNA monomers that produced detectable signals on X1X2Y sex chromosomes of karyomorph D (HMA-D). SatDNA family names are indicated on the right bottom, in green (Atto-488-dUTP) or red (Atto-550-dUTP). The arrows point on the X and Y chromosomes. Bar = 5 μm.

Minimum spanning trees (MST)

We selected two satDNA monomers to generate minimum spanning trees: (1) the HmaSat02, as the satDNA with the most pronounced hybridization pattern, with telomeric and interstitial signals located on all chromosomes of the HMA-D complement, including multiple sex chromosome system, and (2) the HmaSat18, which displayed terminal signals on several autosomes and in the X2 and neo-Y chromosomes of this karyomorph. The MST of HmaSat02 revealed the predominance of three abundant haplotypes, besides the presence of many monomers with few mutations. All of them displayed a very similar abundance to all haplotypes identified in males and females (Fig. 5). The MST of HmaSat18 revealed the presence of many haplotypes. Although the predominant ones were shared between males and females, unique male-specific haplotypes were observed (Fig. 6).Fig. 5 Linear MST of HmaSat02 satDNA obtained from female (red) and male (green) reads. The diameter of the circle is proportional to the abundance of the haplotype. Black circles represent one base of divergence between haplotypes.

Fig. 6 Linear MST of HmaSat18 satDNA obtained from female (red) and male (green) reads. The diameter of the circle is proportional to the abundance of the haplotype. Black circles represent one base of divergence between haplotypes.

Discussion

In this work, we studied the satDNA content and its chromosomal distribution in one cytotype (karyomorph D) of the wolf-fish H. malabaricus with multiple sex chromosome system of the X1X2Y type. The main purposes of this study were to investigate, (1) whether certain satDNAs could be associated with the breakpoint regions involved in the Y-autosome fusion which gave rise to the mentioned multiple sex chromosome system, and (2) whether the neo-Y chromosome already underwent a significant degree of differentiation that could be reflected by the satDNA accumulation. The opportunity to compare the satDNA distribution also with the patterns on the XY sex chromosomes of closely related karyomorph (HMA-C), from which X1X2Y of HMA-D originated via Y-autosome fusion74, qualifies H. malabaricus as a useful study model for these research questions.

A recent burst in satellitome studies has led to significant advancements in our understanding of the teleost satDNA evolutionary dynamics and its possible impact on various genomic processes66,68–70,77. Thus far, a particular attention has been given to characiform fishes64,66–68,77. Among the characiform satellitomes, the diversity of satDNA ranges from 30 (Piaractus mesopotamicus) to an exceptional diversity of 164 monomers (Megaleporinus macrocephalus)68,78. In the present study, we revealed 56 satDNA monomers in the genome of H. malabaricus—karyomorph D (Supplementary Table 1), which falls within the above-mentioned range reported for characiforms. We further revealed an extraordinary wide range of satDNA monomer sizes, ranging from 6 to 2774 bp, with the predominance of those exhibiting higher proportion of AT-pairs (> 50%) (Supplementary Table 1), similarly to what has been found in some other vertebrates such as mammals and amphibians79–81, but differing from reptiles and birds with a predominantly CG-rich content82,83.

Sixteen out of 18 chromosomally mapped HmaSatDNA monomers showed detectable clusters pointing on their organization into high-copy number arrays. We observed the following three major distribution patterns: (1) only two hybridization signals restricted to a certain autosome pair (nine monomers; Fig. 1); (2) clusters located on the X1, X2, and/or neo-Y sex chromosomes in addition to several-to-all autosomes (six monomers; Figs. 2, 7), and (3) a single satDNA, HmaSat07, accumulated exclusively in the pericentromeric regions and surrounding area on the X1 and neo-Y sex chromosome (Figs. 2, 7).Fig. 7 Partial karyotypes of sex chromosomes from H. malabaricus karyomorphs HMA-D and HMA-C after FISH with selected satDNA monomers, including the patterns on X1X2Y sex chromosome meiotic trivalent (HMA-D). The sex chromosomes were retrieved from chromosome spreads presented in Figs. 2, 3 and 4. Bar = 5 μm.

Intriguingly, although the most abundant satDNA sequence in the genome very often occupies centromeric regions of all chromosomes, implying its functional significance for this crucial chromosomal region84, our data on distribution patterns imply this is not the case for HmaSat01, the most abundant satDNA monomer in HMA-D. Previously, the 5SHindIII-DNA satDNA, corresponding to the herein characterized HmaSat03 monomer, has been proposed to play a role in centromere function or structure in this fish76, however, our current analysis showed that other satDNA, namely HmaSat05, is a more likely candidate for this role based on the higher proportion of (peri)centromeric signals per chromosome set (Fig. 2). Yet, this monomer does not form detectable clusters on each centromeric region in the complement, with a single autosome pair (putatively the autosomal addition creating X1X2Y system in HMA-D) lacking the signal in HMA-C, and more chromosomes without the signal in HMA-D (Figs. 1, 2, 4). Our present findings therefore suggest a gradual change in the (peri)centromeric satDNA content in H. malabaricus, with analogous pattern being recently reported for African Nothobranchius killifishes85. In H. malabaricus, this process seems to be considerably fast given that HmaSat05 monomer occupies the centromeric regions only in HMA-C and HMA-D, while it is not present in these regions on chromosomes of closely related karyomorphs A, F, and G (Supplementary Fig. 4). It is generally thought that the mechanisms driving fast turnover of centromeric repeats and associated proteins may be responsible for a reproductive barrier between diverging populations of the same or closely related species, due to hybrid incompatibility related to chromosome segregation86–88. This, together with possible incompatibilities related to different sex chromosome systems89, might have contributed to the previously reported absence of hybrids between H. malabaricus karyomorphs90, thus supporting the proposal that the karyomorphs of H. malabaricus represent already distinct evolutionary units71. A functional role of HmaSat05 monomer in HMA-C and HMA-D centromeres will need to be further assessed by testing the interaction between this satDNA and key centromeric proteins such as CENP-A52,58.

From this point on, we will attempt to address the two main queries raised beforehand, starting with the inquiry of whether changes in satDNAs content accompanied the rearrangement(s) involved in the origin of the X1X2Y multiple sex chromosome system present in HMA-D.

The association between repetitive DNA accumulation and the occurrence of chromosomal rearrangements along with overall elevated genome instability has been extensively recorded in diverse lineages across the tree of life91,92, including humans93. In the case of X1X2Y sex chromosome system of HMA-D, most neo-Y-linked clusters of HmaSatDNA monomers are concentrated in the centromeric region and its surrounding area (Fig. 7). From the herein mapped monomers, only HmaSat02 and HmaSat03 repeats populate the putative breakpoint junction of the Y-autosome tandem fusion (Figs. 2 and 3), which has been inferred from previous WCP experiments72, and comparisons of the position of repetitive DNA clusters between the X2 and neo-Y (HMA-D) and the ancestral XY sex chromosomes in HMA-C74. Whether or not these sequences directly contributed to the rearrangement remains an open question. The other monomers, HmaSat04 and HmaSat05, forming strong (peri)centromeric accumulations on X1 and neo-Y, are likely absent from the ancestor of the X2 chromosome as demonstrated by the lack of FISH signals in Figs. 2 and 3. Contrary to other studies61,62, which evidenced a high rate of elimination and turnover of satDNAs populating neo-sex chromosomes, our findings indicate that both HmaSat04 and HmaSat05 satDNAs were not excised during the Y-autosome fusion event or lost through the subsequent process of centromere inactivation/elimination (reviewed by Barra and Fachinetti94, and McKinley and Cheeseman95).

The second question was whether or not the neo-Y chromosome underwent a significant degree of differentiation, which could be inferred from the level of satDNA accumulation. For this, we took the opportunity of direct comparison of satDNA landscapes between the multiple X11X2Y sex chromosomes of HMA-D and standard XY sex chromosome system of HMA-C from which the system of HMA-D evolved. Species with heteromorphic sex chromosomes usually exhibit sex-linked variations in the accumulation of certain repetitive DNA sequences; they can be either substantially enriched or even limited to sex chromosome(s)31,32,43,61,62, which has been reported also in teleosts64,67,68. When narrowed down to multiple sex chromosomes, those so far revealed in teleost fishes show, in the overwhelming majority, low degree of repeat accumulation2. In the case of X1X2Y system of HMA-D, our results corroborate the previous findings from other FISH-based studies73–75 by showing non-elevated degree of satDNA accumulation, with a limited difference between the synteny blocks of X1/X2 versus neo-Y chromosome (Supplementary Table 1; Fig. 7). Minor changes linked to neo-Y satDNA content may be related to the loss of the HmaSat01 accumulation from the terminal portion of q-arms (on both neo-Y and X1 chromosome), based on the comparison with patterns found in the XY sex chromosome system of HMA-C (see Fig. 7). Nevertheless, it is not possible to rule out the alternative hypothesis that HmaSat01 repeats populated/amplified to cytologically detectable extent on the XY sex chromosomes of HMA-C more recently, i.e. the sex chromosome-autosome fusion in HMA-D might predate this accumulation.

The causal relationship between repetitive DNA accumulation and recombination suppression on sex chromosomes has been a matter of a long-standing debate. It seems more likely that repetitive DNA accumulates in the region of suppressed recombination and then it contributes to further expansion of this region, to genetic decay, and to sex chromosome heteromorphy4,30–35,96. Important question is then: what happens during the first meiotic division when sex chromosomes need to pair efficiently despite possible differences in size, morphology and/or repetitive DNA/heterochromatin content? Accumulation of repeated DNAs and heterochromatin could hamper the establishment and stability of the trivalent, preventing thereby its proper segregation, leading to high frequency of unbalanced gametes and thus being avoided by natural selection20–22,28,97,98. Bearing this in mind, we speculate that the lack of major differences in satDNA accumulation between the members of X1X2Y sex chromosome system of HMA-D, along with generally low degree of repeat accumulation might have contributed to faithful pairing and segregation of these sex chromosomes (Fig. 7). The presence of sex chromosome-specific repeat segments in the area of assumed sequence differentiation, both in the XY sex chromosome system of HMA-C and X1X2Y system of HMA-D, might facilitate smooth sex chromosome pairing within this region, without asynapsis, as has been already proposed for other studied teleost fishes99,100.

Conclusions

Taken together, our present study provides novel insight into the aspects linked to emergence and differentiation of representative fish multiple sex chromosome system—X1X2Y present in H. malabaricus, karyomorph D. We showed that the satDNA content remained rather unchanged during the evolution of this system after Y-autosome fusion. Last but not least, by unravelling the satellitome features of this H. malabaricus karyomorph, we provided new tools for further studies on satDNA evolution in characiform fishes.

Material and methods

Sampling and chromosome preparation

We analyzed 20 individuals of H. malabaricus karyomorph HMA-C (11 males; 9 females) and 17 individuals of karyomorph HMA-D (10 males; 7 females). HMA-C individuals were collected in Bento Gomes River (Poconé, Mato Grosso State, Brazil) and HMA-D in the Monjolinho stream (São Carlos, São Paulo State, Brazil). The sex of sampled individuals was identified by morphological analysis of gonadal tissues, and each individual was assigned to specific karyomorph based on the characteristic cytogenetic profiles (e.g., 2n, karyotype structure) uncovered in previous studies71,73,74. Sampling was authorized by the Brazilian Environmental Agency ICMBIO/SISBIO (License 48628-14) and SISGEN (A96FF09). Mitotic chromosomes were obtained from kidney biopsies following the protocol by Bertollo et al.101, while meiotic spreads were prepared from testes of HMA-D males according to Kligerman and Bloom102. Briefly, we treated the animals with benzocaine and applied a colchicine solution (0.025%) for 40 min. Subsequently, we excised kidney samples and testicular tissue, which were fragmented and treated with a hypotonic solution (0.075 M KCl) for 20 min at 37 °C. Both mitotic cell suspensions and fragmented testes were preserved in Carnoy I fixative (3:1 methanol/glacial acetic acid). Chromosome slides were stained with a 10% Giemsa solution (pH 6.8). All experiments followed the ethical guidelines sanctioned by the Ethics Committee on Animal Experimentation of the Universidade Federal de São Carlos (Process Number 7994170423). The authors complied with ARRIVE guidelines.

DNA extraction and next generation sequencing

Genomic DNA (gDNA) of one male and one female of HMA-D was extracted using liver tissue by the standard phenol–chloroform-isoamylalcohol procedure103. Illumina paired-end libraries with 150 bp read length were prepared from the isolated HMA-D gDNAs and sequenced on the BGISEQ-500 platform at BGI (BGI Shenzhen Corporation, Shenzhen, China). The sequencing yielded a total of 1.55 Gb and 1.4 Gb of raw reads for males and females, respectively. The genomic reads obtained were deposited in the Sequence Read Archive (SRA) under accession numbers SRR27786732 (male) and SRR27786731 (female).

Analysis of H. malabaricus satDNAs

We characterized the satellite DNA content of HMA-D male and female independently. First, we examined the read quality (Q > 20 for all nucleotides), by using the Trimmomatic software version 3.0104. Then, a subset of 2 × 500,000 reads has been randomly subsampled and usedas input for the TAREAN analysis51. All recovered satDNAs were removed from the original library by using the DeconSeq software version 0.4.3105, and from the remaining dataset, a new randomly selected subset of 2 × 500,000 reads was created to pass through a second round of the TAREAN analysis. These steps (iterations) were repeated until no novel satDNAs were found. In the next step, we performed a screening and the removal of the non-satDNA repetitive DNA classes (e.g., multigene families, transposable elements etc.) from our dataset. Finally, we performed a homology search with RM_homology (https://github.com/fjruizruano/satminer/blob/master/rm_homology.py) to group the sequences into variants (> 95% similar), families (80–95% similar), and superfamilies (50–80% similar)49. The catalog of satellite DNAs was deposited on the GenBank with accession numbers PP262651-PP262706.

Estimating the abundance and genetic distances of satDNAs

The abundance of each satDNA (proportion of reads containing the same repeat vs. the entire pool of reads) was estimated by selecting 2 × 500,000 reads from each genomic library (HMA-D female and HMA-D male). Then, using the “cross-match” tool in the RepeatMasker version 4.1.6 106 and a custom python script (https://github.com/fjruizruano/ngsprotocols/blob/master/repeat_masker_run_big.py, accessed on 20/01/2023), the reads were aligned against each respective catalog. We used the script calcDivergenceFromAlign.py in the RepeatMasker106 to calculate the genetic distance of each satDNA (Kimura-2-parameter). All satDNAs were named “HmaSat” and numbered according to their decreasing abundance in the genome. Finally, to check for the possible presence of male-biased satDNAs, we calculated the male/female (M/F) abundance ratio of each HmaSat monomer given by RepeatMasker and selected those with an M/F ratio greater than 1.0.

Minimum spanning trees (MST)

To achieve the accurate scores of monomer diversity, we used PHYLOVIZ version 2.0107 to construct the minimum spanning trees (MSTs) for two satDNA monomers which were found on the sex chromosomes of HMA-D and, at the same time, their monomer sizes fell into the required range for the analysis, i.e. below 150 pb (reflecting the size of the sequencing raw reads): HmaSat02 and HmaSat18. First, we subsampled 2 × 5,000,000 reads from each genomic library (males and females HMA-D) and aligned the reads using Bowtie2 version 1.3.1108 against the two aforementioned satDNA monomers. Then, we removed the central area that corresponded to one monomer and, in order to avoid sequencing error bias by CD-HIT version 4.8.1109, we also removed singletons (monomers detected only once) and aligned each satDNA independently using MUSCLE version 5.0110.

Primer design and DNA amplification via polymerase chain reaction (PCR)

Out of the 56 HmaSatDNA monomers found by our analysis (see Supplementary Table 1), we selected 18 representative monomers (five most abundant ones and those with some differences in the abundance between sexes; indicated by asterisk in Supplementary Table 1) for physical mapping. For this purpose, we manually designed primers for each consensus monomer, using both multiple primer analyzer (Thermofisher) and OligoCalc (Biotools) to check for the hairpin formation and self-annealing. The PCR conditions followed Kretschmer et al.67 and it was performed in a T100 Bio-Rad Thermal Cycler (Bio-Rad Laboratories). Briefly, the amplification reaction consisted of an initial denaturation at 95 °C for 7 min, followed by 34 cycles with denaturation at 95 °C (45 s), annealing temperatures varying from 50 to 62 °C (60 s), extension at 72 °C (60 s), and the final extension at 72 °C (7 min). All PCR products were quantified using a ThermoFisher NanoDrop spectrophotometer (ThermoFisher Scientific) and submitted to electrophoresis with a 1% or 2% agarose gel (Supplementary Fig. 5).

Fluorescence in situ hybridization (FISH)

To examine the chromosomal location of HmaSatDNA monomers, with the emphasis on their landscapes on the nascent XX/XY (HMA-C) and multiple X1X1X2X2/X1X2Y sex chromosomes (HMA-D), we performed FISH on both male and female metaphase spreads of the two karyomorphs. To unambiguously identify the sex chromosomes in both karyomorphs, we took the advantage of the previously presence of 18S ribosomal DNA (rDNA) on the neo-Y and X1 chromosome of HMA-D, and XY chromosomes of HMA-C74 and the characteristic distribution patterns of one satDNA monomer76. In addition, the X2 chromosome was identified according to its morphology and size, as it is the smallest submetacentric chromosome in the karyotype71,74. For FISH, we used the 18S rDNA probe (1400 pb-long segment of the 18S rRNA gene from H. malabaricus) was previously isolated by Cioffi et al.111. Finally, we also used the sex chromosome-specific marker, HmaSat07 satDNA described in the present study, in complementary experiments. Results from the sequential experiments (i.e. on the same metaphases) with 18S rDNA or HmaSat07, and each individual satDNA are provided in Supplementary Figs. 6 and 7. The HmaSat DNAs that displayed signals on the mitotic X1, X2, or neo-Y chromosomes of HMA-D, were mapped also onto male meiotic spreads in prophase I of this karyomorph and further onto HMA-C metaphases. All the 18 tested satDNA monomers along with the 18S rDNA fragment were labeled with Atto-488-dUTP; Atto-425-dUTP (green fluorescence), or Atto-550-dUTP (red fluorescence), using a nick-translation kit (Jena Bioscience, Jena, Germany) according to the manufacturer's instructions. The hybridization procedures were conducted under high-stringency conditions. After initial slide pretreatment steps with an RNAse (10 µg/ml in 2 × SSC solution) for 1 h 30 min at 37 °C and pepsin (50 µg/mL in 10 mM HCl) for 10 min, followed by denaturation in 70% formamide at 72 °C for 3 min, 15 s, and hybridization in a 37 °C moist chamber for at least 14 h112. Finally, the chromosome slides were stained with 4′,6-diamidino-2-phenylindole (DAPI) solution.

Microscopy and image processing

To verify the FISH patterns, 2n and karyotype structure of each karyomorph, we analyzed at least 30 metaphases per individual. The images were captured using an Olympus BX50 microscope (Olympus Corporation, Ishikawa, Japan) with CoolSNAP, and processed using the Image-Pro Plus 4.1 software (Media Cybernetics, Silver Spring, USA). We classified chromosomes according to their arm ratios113 as metacentric (m), submetacentric (sm), subtelocentric (st) or acrocentric (a).

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70920-7.

Acknowledgements

We are grateful to the staff of the Laboratory of “Molekulare Zytogenetik” (Institut für Humangenetik; Jena, Germany), specially: Stefanie Kankel; Luiza Person and Niklas Padutsch for assisting with the experiment's workflow. We also thank Luiz Henrique da Silva for helping with sample collection.

Author contributions

G.A.T., R.U. and M.B.C., conceived and designed research. G.A.T., C.A.G.G. and R.K., conducted experiments. G.A.T., A.S., R.U., M.B.C., C.A.G.G., and R.K., analyzed the data. G.A.T., R.U. and C.A.G.G., contributed with new methods. G.A.T., A.S., R.U., M.B.C., C.A.G.G., R.K., F.P.F., L.A.C.B. and T.L., wrote the paper.

Funding

Open Access funding enabled and organized by Projekt DEAL. This work was supported by São Paulo Research Foundation (FAPESP) grants 2022/00427-3 (G.A.T); 2023/00955-2 (M.B.C) and Brazilian National Council for Scientific and Technological Development (CNPq) grant number 302928/2021-9 (M.B.C). We acknowledge support by the German Research Foundation Projekt-Nr. 512648189 (T.L.), RVO: 67985904 of IAPG CAS, Liběchov (Czech Academy of Sciences) (A.S.) and the Open Access Publication Fund of the Thueringer Universitaetsund Landesbibliothek Jena. All authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. The catalog of satellite DNAs was deposited on the GenBank with Accession Numbers PP262651-PP262706, and raw reads are available in Sequence Read Archive (SRA-NCBI) under Accession Numbers SRR27786732 (male) and SRR27786731 (female). All other relevant data are within the paper and its supplementary material.

Competing interests

The authors declare no competing interests.

Ethical approval

Our sampling was authorized by the Brazilian Environmental Agency ICMBIO/SISBIO (License 48628-14) and SISGEN (A96FF09). All experiments followed the ethical guidelines sanctioned by the Ethics Committee on Animal Experimentation of the Universidade Federal de São Carlos (Process number 7994170423). The authors complied with ARRIVE guidelines.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Pennell MW Y Fuse? Sex chromosome fusions in fishes and reptiles PLoS Genet. 2015 11 e1005237 10.1111/j.1601-5223.1964.tb01953.x 25993542
Pennell, M. W. et al. Y Fuse? Sex chromosome fusions in fishes and reptiles. PLoS Genet. 11, e1005237. 10.1111/j.1601-5223.1964.tb01953.x (2015).25993542 10.1111/j.1601-5223.1964.tb01953.x
2. Sember A Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges Phil. Trans. R. Soc. B Biol. Sci. 2021 376 20200098 10.1098/rstb.2020.0098
Sember, A. et al. Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges. Phil. Trans. R. Soc. B Biol. Sci. 376, 20200098. 10.1098/rstb.2020.0098 (2021).10.1098/rstb.2020.0098
3. Charlesworth D When and how do sex-linked regions become sex chromosomes? Evolution 2021 75 569 581 10.1111/evo.14196 33592115
Charlesworth, D. When and how do sex-linked regions become sex chromosomes? Evolution 75, 569–581. 10.1111/evo.14196 (2021).33592115 10.1111/evo.14196
4. Bergero R Charlesworth D The evolution of restricted recombination in sex chromosomes Trends Ecol. Evol. 2009 24 94 102 10.1016/j.tree.2008.09.010 19100654
Bergero, R. & Charlesworth, D. The evolution of restricted recombination in sex chromosomes. Trends Ecol. Evol. 24, 94–102. 10.1016/j.tree.2008.09.010 (2009).19100654 10.1016/j.tree.2008.09.010
5. Blackmon H Demuth JP The fragile Y hypothesis: Y chromosome aneuploidy as a selective pressure in sex chromosome and meiotic mechanism evolution BioEssays 2015 37 942 950 10.1002/bies.201500040 26200104
Blackmon, H. & Demuth, J. P. The fragile Y hypothesis: Y chromosome aneuploidy as a selective pressure in sex chromosome and meiotic mechanism evolution. BioEssays 37, 942–950. 10.1002/bies.201500040 (2015).26200104 10.1002/bies.201500040
6. Kitano J A role for a neo-sex chromosome in stickleback speciation Nature 2009 461 1079 1083 10.1038/nature08441 19783981
Kitano, J. et al. A role for a neo-sex chromosome in stickleback speciation. Nature 461, 1079–1083. 10.1038/nature08441 (2009).19783981 10.1038/nature08441
7. Beaudry FEG Barrett SCH Wright SI Ancestral and neo-sex chromosomes contribute to population divergence in a dioecious plant Evolution 2020 74 256 269 10.1111/evo.13892 31808547
Beaudry, F. E. G., Barrett, S. C. H. & Wright, S. I. Ancestral and neo-sex chromosomes contribute to population divergence in a dioecious plant. Evolution 74, 256–269. 10.1111/evo.13892 (2020).31808547 10.1111/evo.13892
8. Wang S Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier Proc. Natl. Acad. U S A 2022 119 e2119382119 10.1073/pnas.2119382119
Wang, S. et al. Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier. Proc. Natl. Acad. U S A 119, e2119382119. 10.1073/pnas.2119382119 (2022).10.1073/pnas.2119382119
9. Smith DAS Neo Sex Chromosomes, colour polymorphism and male-killing in the African queen butterfly, Danaus chrysippus (L.) Insects 2019 10 291 10.3390/insects10090291 31505824
Smith, D. A. S. et al. Neo sex chromosomes, colour polymorphism and male-killing in the African queen butterfly, Danaus chrysippus (L.). Insects 10, 291. 10.3390/insects10090291 (2019).31505824 10.3390/insects10090291
10. Charlesworth D Charlesworth B Sex differences in fitness and selection for centric fusions between sex-chromosomes and autosomes Genet. Res. 1980 35 205 214 10.1017/S0016672300014051 6930353
Charlesworth, D. & Charlesworth, B. Sex differences in fitness and selection for centric fusions between sex-chromosomes and autosomes. Genet. Res. 35, 205–214. 10.1017/S0016672300014051 (1980).6930353 10.1017/S0016672300014051
11. Matsumoto T Kitano J The intricate relationship between sexually antagonistic selection and the evolution of sex chromosome fusions J. Theor. Biol. 2016 404 97 108 10.1016/j.jtbi.2016.05.036 27259387
Matsumoto, T. & Kitano, J. The intricate relationship between sexually antagonistic selection and the evolution of sex chromosome fusions. J. Theor. Biol. 404, 97–108. 10.1016/j.jtbi.2016.05.036 (2016).27259387 10.1016/j.jtbi.2016.05.036
12. Mora P Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies Mol. Ecol. 2024 5 e17256 10.1111/mec.17256
Mora, P. et al. Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies. Mol. Ecol. 5, e17256 (2024).10.1111/mec.17256
13. Pokorná M Altmanová M Kratochvíl L Multiple sex chromosomes in the light of female meiotic drive in amniote vertebrates Chromosome Res. 2014 22 35 44 10.1007/s10577-014-9403-2 24590843
Pokorná, M., Altmanová, M. & Kratochvíl, L. Multiple sex chromosomes in the light of female meiotic drive in amniote vertebrates. Chromosome Res. 22, 35–44. 10.1007/s10577-014-9403-2 (2014).24590843 10.1007/s10577-014-9403-2
14. Ma W-J Veltsos P The diversity and evolution of sex chromosomes in frogs Genes 2021 12 483 10.3390/genes12040483 33810524
Ma, W.-J. & Veltsos, P. The diversity and evolution of sex chromosomes in frogs. Genes 12, 483. 10.3390/genes12040483 (2021).33810524 10.3390/genes12040483
15. Charlesworth D Some thoughts about the words we use for thinking about sex chromosome evolution Phil. Trans. R. Soc. B Biol. Sci. 2022 377 20210314 10.1098/rstb.2021.0314
Charlesworth, D. Some thoughts about the words we use for thinking about sex chromosome evolution. Phil. Trans. R. Soc. B Biol. Sci. 377, 20210314. 10.1098/rstb.2021.0314 (2022).10.1098/rstb.2021.0314
16. de Araújo L Cytogenetic and molecular characterization of Eigenmannia aff . desantanai (Gymnotiformes: Sternopygidae): A first report of system of sex chromosomes ZW1W2/ZZ in Gymnotiformes Zebrafish 2023 20 77 85 10.1089/zeb.2022.0059 36930749
de Araújo, L. et al. Cytogenetic and molecular characterization of Eigenmannia aff. desantanai (Gymnotiformes: Sternopygidae): A first report of system of sex chromosomes ZW1W2/ZZ in Gymnotiformes. Zebrafish 20, 77–85. 10.1089/zeb.2022.0059 (2023).36930749 10.1089/zeb.2022.0059
17. Marajó L Chromosomal rearrangements and the first indication of an ♀X1X1X2X2/♂X1X2Y sex chromosome system in Rineloricaria fishes (Teleostei: Siluriformes) J. Fish Biol. 2023 102 443 454 10.1111/jfb.15275 36427042
Marajó, L. et al. Chromosomal rearrangements and the first indication of an ♀X1X1X2X2/♂X1X2Y sex chromosome system in Rineloricaria fishes (Teleostei: Siluriformes). J. Fish Biol. 102, 443–454. 10.1111/jfb.15275 (2023).36427042 10.1111/jfb.15275
18. Ferchaud A-L Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut ( Reinhardtius hippoglossoides ) G3 (Bethesda) 2022 12 376 10.1093/g3journal/jkab376
Ferchaud, A.-L. et al. Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides). G3 (Bethesda) 12, 376. 10.1093/g3journal/jkab376 (2022).10.1093/g3journal/jkab376
19. Sassi F Homeology of sex chromosomes in Amazonian Harttia armored catfishes supports the X-fission hypothesis for the X1X2Y sex chromosome system origin Sci. Rep. 2023 13 15756 10.1038/s41598-023-42617-w 37735233
Sassi, F. et al. Homeology of sex chromosomes in Amazonian Harttia armored catfishes supports the X-fission hypothesis for the X1X2Y sex chromosome system origin. Sci. Rep. 13, 15756. 10.1038/s41598-023-42617-w (2023).37735233 10.1038/s41598-023-42617-w
20. Sciurano RB Rahn MI Rey-Valzacchi G Coco R Solari AJ The role of asynapsis in human spermatocyte failure Int. J. Androl. 2012 35 541 549 10.1111/j.1365-2605.2011.01221.x 21977946
Sciurano, R. B., Rahn, M. I., Rey-Valzacchi, G., Coco, R. & Solari, A. J. The role of asynapsis in human spermatocyte failure. Int. J. Androl. 35, 541–549. 10.1111/j.1365-2605.2011.01221.x (2012).21977946 10.1111/j.1365-2605.2011.01221.x
21. Wolf KW How meiotic cells deal with non-exchange chromosomes BioEssays 1994 16 107 114 10.1002/bies.950160207 8147841
Wolf, K. W. How meiotic cells deal with non-exchange chromosomes. BioEssays 16, 107–114. 10.1002/bies.950160207 (1994).8147841 10.1002/bies.950160207
22. Nokkala S Kuznetsova VG Maryanska-Nadachowska A Nokkala C Holocentric chromosomes in meiosis. II. The modes of orientation and segregation of a trivalent Chromosome Res. 2006 14 559 565 10.1007/s10577-006-1053-6 16823618
Nokkala, S., Kuznetsova, V. G., Maryanska-Nadachowska, A. & Nokkala, C. Holocentric chromosomes in meiosis. II. The modes of orientation and segregation of a trivalent. Chromosome Res. 14, 559–565. 10.1007/s10577-006-1053-6 (2006).16823618 10.1007/s10577-006-1053-6
23. Noronha RCR Nagamachi CY O’Brien PCM Ferguson-Smith MA Pieczarka JC Neo-XY body: An analysis of XY1Y2 meiotic behavior in Carollia (Chiroptera, Phyllostomidae) by chromosome painting Cytogenet. Genome Res. 2009 124 37 43 10.1159/000200086 19372667
Noronha, R. C. R., Nagamachi, C. Y., O’Brien, P. C. M., Ferguson-Smith, M. A. & Pieczarka, J. C. Neo-XY body: An analysis of XY1Y2 meiotic behavior in Carollia (Chiroptera, Phyllostomidae) by chromosome painting. Cytogenet. Genome Res. 124, 37–43. 10.1159/000200086 (2009).19372667 10.1159/000200086
24. Castillo ERD Bidau CJ Martí DA Neo-sex chromosome diversity in Neotropical melanopline grasshoppers (Melanoplinae, Acrididae) Genetica 2010 138 775 786 10.1007/s10709-010-9458-8 20411310
Castillo, E. R. D., Bidau, C. J. & Martí, D. A. Neo-sex chromosome diversity in Neotropical melanopline grasshoppers (Melanoplinae, Acrididae). Genetica 138, 775–786. 10.1007/s10709-010-9458-8 (2010).20411310 10.1007/s10709-010-9458-8
25. Poggio MG Gaspe MS Papeschi AG Bressa MJ Cytogenetic study in a mutant of Triatoma infestans (Hemiptera: Reduviidae) carrying a spontaneous autosomal fusion and an extra chromosome Cytogenet. Genome Res. 2012 139 44 51 10.1159/000342875 23037941
Poggio, M. G., Gaspe, M. S., Papeschi, A. G. & Bressa, M. J. Cytogenetic study in a mutant of Triatoma infestans (Hemiptera: Reduviidae) carrying a spontaneous autosomal fusion and an extra chromosome. Cytogenet. Genome Res. 139, 44–51. 10.1159/000342875 (2012).23037941 10.1159/000342875
26. Farooq U Lovleen Saggoo MIS Male meiosis and behaviour of sex chromosomes in different populations of Rumex acetosa L. from the Western Himalayas, India Plant Syst. Evol. 2014 300 287 294 10.1007/s00606-013-0881-z
Farooq, U., Lovleen, & Saggoo, M. I. S. Male meiosis and behaviour of sex chromosomes in different populations of Rumex acetosa L. from the Western Himalayas, India. Plant Syst. Evol. 300, 287–294. 10.1007/s00606-013-0881-z (2014).10.1007/s00606-013-0881-z
27. Sember A Patterns of sex chromosome differentiation in spiders: Insights from comparative genomic hybridisation Genes 2020 11 849 10.3390/genes11080849 32722348
Sember, A. et al. Patterns of sex chromosome differentiation in spiders: Insights from comparative genomic hybridisation. Genes 11, 849. 10.3390/genes11080849 (2020).32722348 10.3390/genes11080849
28. Bertollo LAC Mestriner CA The X1X2Y sex chromosome system in the fish Hoplias malabaricus II. Meiotic analyses Chromosome Res. 1998 6 141 147 10.1023/A:1009243114124 9543017
Bertollo, L. A. C. & Mestriner, C. A. The X1X2Y sex chromosome system in the fish Hoplias malabaricus II. Meiotic analyses. Chromosome Res. 6, 141–147. 10.1023/A:1009243114124 (1998).9543017 10.1023/A:1009243114124
29. da Silva M Repetitive DNA and meiotic behavior of sex chromosomes in Gymnotus pantanal (Gymnotiformes, Gymnotidae) Cytogenet. Genome Res. 2011 135 143 149 10.1159/000330777 21860225
da Silva, M. et al. Repetitive DNA and meiotic behavior of sex chromosomes in Gymnotus pantanal (Gymnotiformes, Gymnotidae). Cytogenet. Genome Res. 135, 143–149. 10.1159/000330777 (2011).21860225 10.1159/000330777
30. Charlesworth D Charlesworth B Marais G Steps in the evolution of heteromorphic sex chromosomes Heredity 2005 95 118 128 10.1038/sj.hdy.6800697 15931241
Charlesworth, D., Charlesworth, B. & Marais, G. Steps in the evolution of heteromorphic sex chromosomes. Heredity 95, 118–128. 10.1038/sj.hdy.6800697 (2005).15931241 10.1038/sj.hdy.6800697
31. Ezaz T Deakin JE Repetitive sequence and sex chromosome evolution in vertebrates Adv. Evol. Biol. 2014 2014 1 9 10.1155/2014/104683
Ezaz, T. & Deakin, J. E. Repetitive sequence and sex chromosome evolution in vertebrates. Adv. Evol. Biol. 2014, 1–9. 10.1155/2014/104683 (2014).10.1155/2014/104683
32. Hobza R Impact of repetitive elements on the Y chromosome formation in plants Genes 2017 8 302 10.3390/genes8110302 29104214
Hobza, R. et al. Impact of repetitive elements on the Y chromosome formation in plants. Genes 8, 302. 10.3390/genes8110302 (2017).29104214 10.3390/genes8110302
33. Steinemann S Steinemann M Y chromosomes: Born to be destroyed BioEssays 2005 27 1076 1083 10.1002/bies.20288 16163733
Steinemann, S. & Steinemann, M. Y chromosomes: Born to be destroyed. BioEssays 27, 1076–1083. 10.1002/bies.20288 (2005).16163733 10.1002/bies.20288
34. Kent TV Uzunović J Wright SI Coevolution between transposable elements and recombination Phil. Trans. R. Soc. B Biol. Sci. 2017 372 20160458 10.1098/rstb.2016.0458
Kent, T. V., Uzunović, J. & Wright, S. I. Coevolution between transposable elements and recombination. Phil. Trans. R. Soc. B Biol. Sci. 372, 20160458. 10.1098/rstb.2016.0458 (2017).10.1098/rstb.2016.0458
35. Kratochvíl L Expanding the classical paradigm: what we have learnt from vertebrates about sex chromosome evolution. Phil. Trans. R. Soc. B. Biol. Sci. 2021 376 20200097 10.1098/rstb.2020.0097
Kratochvil, L. et al. Expanding the classical paradigm: what we have learnt from vertebrates about sex chromosome evolution. Phil. Trans. R. Soc. B. Biol. Sci. 376, 20200097 (2021).10.1098/rstb.2020.0097
36. Huang Z Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots Nat. Commun. 2022 13 944 10.1038/s41467-022-28585-1 35177601
Huang, Z. et al. Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots. Nat. Commun. 13, 944. 10.1038/s41467-022-28585-1 (2022).35177601 10.1038/s41467-022-28585-1
37. Flynn JM Hu KB Clark AG Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content Genetics 2023 224 062 10.1093/genetics/iyad062
Flynn, J. M., Hu, K. B. & Clark, A. G. Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content. Genetics 224, 062. 10.1093/genetics/iyad062 (2023).10.1093/genetics/iyad062
38. Dobigny G Ozouf-Costaz C Bonillo C Volobouev V Viability of X-autosome translocations in mammals: An epigenomic hypothesis from a rodent case-study Chromosoma 2004 113 34 41 10.1007/s00412-004-0292-6 15243753
Dobigny, G., Ozouf-Costaz, C., Bonillo, C. & Volobouev, V. Viability of X-autosome translocations in mammals: An epigenomic hypothesis from a rodent case-study. Chromosoma 113, 34–41. 10.1007/s00412-004-0292-6 (2004).15243753 10.1007/s00412-004-0292-6
39. Oliveira da Silva W The emergence of a new sex-system (XX/XY1Y2) suggests a species complex in the “monotypic” rodent Oecomys auyantepui (Rodentia, Sigmodontinae) Sci. Rep. 2022 12 8690 10.1038/s41598-022-12706-3 35610291
Oliveira da Silva, W. et al. The emergence of a new sex-system (XX/XY1Y2) suggests a species complex in the “monotypic” rodent Oecomys auyantepui (Rodentia, Sigmodontinae). Sci. Rep. 12, 8690. 10.1038/s41598-022-12706-3 (2022).35610291 10.1038/s41598-022-12706-3
40. Herpin A Transcriptional rewiring of the sex determining dmrt1 gene duplicate by transposable elements PLoS Genet. 2010 6 e1000844 10.1371/journal.pgen.1000844 20169179
Herpin, A. et al. Transcriptional rewiring of the sex determining dmrt1 gene duplicate by transposable elements. PLoS Genet. 6, e1000844. 10.1371/journal.pgen.1000844 (2010).20169179 10.1371/journal.pgen.1000844
41. Wang L Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish BMC Biol. 2022 20 5 10.1186/s12915-021-01205-y 34996452
Wang, L. et al. Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish. BMC Biol. 20, 5. 10.1186/s12915-021-01205-y (2022).34996452 10.1186/s12915-021-01205-y
42. Lohe AR Roberts PA An unusual Y chromosome of Drosophila simulans carrying amplified rDNA spacer without rRNA genes Genetics 1990 125 399 406 10.1093/genetics/125.2.399 2379820
Lohe, A. R. & Roberts, P. A. An unusual Y chromosome of Drosophila simulans carrying amplified rDNA spacer without rRNA genes. Genetics 125, 399–406. 10.1093/genetics/125.2.399 (1990).2379820 10.1093/genetics/125.2.399
43. Garrido-Ramos MA Satellite DNA: An evolving topic Genes 2017 8 230 10.3390/genes8090230 28926993
Garrido-Ramos, M. A. Satellite DNA: An evolving topic. Genes 8, 230. 10.3390/genes8090230 (2017).28926993 10.3390/genes8090230
44. Šatović-Vukšić E Plohl M Satellite DNAs—From localized to highly dispersed genome components Genes 2023 14 742 10.3390/genes14030742 36981013
Šatović-Vukšić, E. & Plohl, M. Satellite DNAs—From localized to highly dispersed genome components. Genes 14, 742. 10.3390/genes14030742 (2023).36981013 10.3390/genes14030742
45. López-Flores I Garrido-Ramos MA Garrido-Ramos MA The repetitive DNA content of eukaryotic genomes Repetitive DNA 2012 Karger Publishers 1 28
López-Flores, I. & Garrido-Ramos, M. A. The repetitive DNA content of eukaryotic genomes. In Repetitive DNA Vol. 7 (ed. Garrido-Ramos, M. A.) 1–28 (Karger Publishers, 2012).
46. Shapiro JA von Sternberg R Why repetitive DNA is essential to genome function Biol. Rev. 2005 80 227 250 10.1017/s1464793104006657 15921050
Shapiro, J. A. & von Sternberg, R. Why repetitive DNA is essential to genome function. Biol. Rev. 80, 227–250. 10.1017/s1464793104006657 (2005).15921050 10.1017/s1464793104006657
47. Louzada S Decoding the role of satellite DNA in genome architecture and plasticity—An evolutionary and clinical affair Genes 2020 11 72 10.3390/genes11010072 31936645
Louzada, S. et al. Decoding the role of satellite DNA in genome architecture and plasticity—An evolutionary and clinical affair. Genes 11, 72. 10.3390/genes11010072 (2020).31936645 10.3390/genes11010072
48. Talbert PB Henikoff S The genetics and epigenetics of satellite centromeres Genome Res. 2022 32 608 615 10.1101/gr.275351.121 35361623
Talbert, P. B. & Henikoff, S. The genetics and epigenetics of satellite centromeres. Genome Res. 32, 608–615. 10.1101/gr.275351.121 (2022).35361623 10.1101/gr.275351.121
49. Ruiz-Ruano FJ López-León MD Cabrero J Camacho JPM High-throughput analysis of the satellitome illuminates satellite DNA evolution Sci, Rep 2016 6 28333 10.1038/srep28333 27385065
Ruiz-Ruano, F. J., López-León, M. D., Cabrero, J. & Camacho, J. P. M. High-throughput analysis of the satellitome illuminates satellite DNA evolution. Sci, Rep 6, 28333. 10.1038/srep28333 (2016).27385065 10.1038/srep28333
50. Camacho JPM Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution BMC Biol. 2022 20 36 10.1186/s12915-021-01216-9 35130900
Camacho, J. P. M. et al. Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution. BMC Biol. 20, 36. 10.1186/s12915-021-01216-9 (2022).35130900 10.1186/s12915-021-01216-9
51. Novák P Neumann P Macas J Global analysis of repetitive DNA from unassembled sequence reads using RepeatExplorer2 Nat. Protoc. 2020 15 3745 3776 10.1038/s41596-020-0400-y 33097925
Novák, P., Neumann, P. & Macas, J. Global analysis of repetitive DNA from unassembled sequence reads using RepeatExplorer2. Nat. Protoc. 15, 3745–3776. 10.1038/s41596-020-0400-y (2020).33097925 10.1038/s41596-020-0400-y
52. Despot-Slade E Satellitome analyses in nematodes illuminate complex species history and show conserved features in satellite DNAs BMC Biol. 2022 20 259 10.1186/s12915-022-01460-7 36397071
Despot-Slade, E. et al. Satellitome analyses in nematodes illuminate complex species history and show conserved features in satellite DNAs. BMC Biol. 20, 259. 10.1186/s12915-022-01460-7 (2022).36397071 10.1186/s12915-022-01460-7
53. Palacios-Gimenez OM Eight million years of satellite DNA evolution in grasshoppers of the genus Schistocerca illuminate the ins and outs of the library hypothesis Genome Biol. Evol. 2020 12 88 102 10.1093/gbe/evaa018 32211863
Palacios-Gimenez, O. M. et al. Eight million years of satellite DNA evolution in grasshoppers of the genus Schistocerca illuminate the ins and outs of the library hypothesis. Genome Biol. Evol. 12, 88–102. 10.1093/gbe/evaa018 (2020).32211863 10.1093/gbe/evaa018
54. de Lima LG Ruiz-Ruano FJ In-depth satellitome analyses of 37 Drosophila species illuminate repetitive DNA evolution in the Drosophila genus Genome Biol. Evol. 2022 14 evac064 10.1093/gbe/evac064 35511582
de Lima, L. G. & Ruiz-Ruano, F. J. In-depth satellitome analyses of 37 Drosophila species illuminate repetitive DNA evolution in the Drosophila genus. Genome Biol. Evol. 14, evac064. 10.1093/gbe/evac064 (2022).35511582 10.1093/gbe/evac064
55. Schmidt N Repeat turnover meets stable chromosomes: Repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets Plant J. 2023 118 171 190 10.1111/tpj.16599 38128038
Schmidt, N. et al. Repeat turnover meets stable chromosomes: Repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets. Plant J. 118, 171–190. 10.1111/tpj.16599 (2023).38128038 10.1111/tpj.16599
56. Belyayev A The major satellite DNA families of the diploid Chenopodium album aggregate species: Arguments for and against the “library hypothesis” PLoS One 2020 15 e0241206 10.1371/journal.pone.0241206 33108401
Belyayev, A. et al. The major satellite DNA families of the diploid Chenopodium album aggregate species: Arguments for and against the “library hypothesis”. PLoS One 15, e0241206. 10.1371/journal.pone.0241206 (2020).33108401 10.1371/journal.pone.0241206
57. Heitkam T Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes Plant J. 2020 103 32 52 10.1111/tpj.14705 31981259
Heitkam, T. et al. Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes. Plant J. 103, 32–52. 10.1111/tpj.14705 (2020).31981259 10.1111/tpj.14705
58. Ávila Robledillo L Extraordinary sequence diversity and promiscuity of centromeric satellites in the legume tribe Fabeae Mol. Biol. Evol. 2020 37 2341 2356 10.1093/molbev/msaa090 32259249
Ávila Robledillo, L. et al. Extraordinary sequence diversity and promiscuity of centromeric satellites in the legume tribe Fabeae. Mol. Biol. Evol. 37, 2341–2356. 10.1093/molbev/msaa090 (2020).32259249 10.1093/molbev/msaa090
59. Bracewell R Chatla K Nalley MJ Bachtrog D Dynamic turnover of centromeres drives karyotype evolution in Drosophila Elife 2019 8 e49002 10.7554/eLife.49002 31524597
Bracewell, R., Chatla, K., Nalley, M. J. & Bachtrog, D. Dynamic turnover of centromeres drives karyotype evolution in Drosophila. Elife 8, e49002. 10.7554/eLife.49002 (2019).31524597 10.7554/eLife.49002
60. Nishihara H Stanyon R Tanabe H Koga A Replacement of owl monkey centromere satellite by a newly evolved variant was a recent and rapid process Genes Cells 2021 26 979 986 10.1111/gtc.12898 34570411
Nishihara, H., Stanyon, R., Tanabe, H. & Koga, A. Replacement of owl monkey centromere satellite by a newly evolved variant was a recent and rapid process. Genes Cells 26, 979–986. 10.1111/gtc.12898 (2021).34570411 10.1111/gtc.12898
61. Cabral-de-Mello DC The spread of satellite DNAs in euchromatin and insights into the multiple sex chromosome evolution in Hemiptera revealed by repeatome analysis of the bug Oxycarenus hyalinipennis Insect Mol. Biol. 2023 32 725 737 10.1111/imb.12868 37615351
Cabral-de-Mello, D. C. et al. The spread of satellite DNAs in euchromatin and insights into the multiple sex chromosome evolution in Hemiptera revealed by repeatome analysis of the bug Oxycarenus hyalinipennis. Insect Mol. Biol. 32, 725–737. 10.1111/imb.12868 (2023).37615351 10.1111/imb.12868
62. Ferretti ABSM Milani D Palacios-Gimenez OM Ruiz-Ruano FJ Cabral-de-Mello DC High dynamism for neo-sex chromosomes: Satellite DNAs reveal complex evolution in a grasshopper Heredity 2020 125 124 137 10.1038/s41437-020-0327-7 32499661
Ferretti, A. B. S. M., Milani, D., Palacios-Gimenez, O. M., Ruiz-Ruano, F. J. & Cabral-de-Mello, D. C. High dynamism for neo-sex chromosomes: Satellite DNAs reveal complex evolution in a grasshopper. Heredity 125, 124–137. 10.1038/s41437-020-0327-7 (2020).32499661 10.1038/s41437-020-0327-7
63. Ruban A Schmutzer T Scholz U Houben A How next-generation sequencing has aided our understanding of the sequence composition and origin of B chromosomes Genes 2017 8 294 10.3390/genes8110294 29068386
Ruban, A., Schmutzer, T., Scholz, U. & Houben, A. How next-generation sequencing has aided our understanding of the sequence composition and origin of B chromosomes. Genes 8, 294. 10.3390/genes8110294 (2017).29068386 10.3390/genes8110294
64. Crepaldi C Martí E Gonçalves ÉM Martí DA Parise-Maltempi PP Genomic differences between the sexes in a fish species seen through satellite DNAs Front. Genet. 2021 12 1885 10.3389/fgene.2021.728670
Crepaldi, C., Martí, E., Gonçalves, É. M., Martí, D. A. & Parise-Maltempi, P. P. Genomic differences between the sexes in a fish species seen through satellite DNAs. Front. Genet. 12, 1885. 10.3389/fgene.2021.728670 (2021).10.3389/fgene.2021.728670
65. Silva DMZDA High-throughput analysis unveils a highly shared satellite DNA library among three species of fish genus Astyanax Sci. Rep. 2017 7 12726 10.1038/s41598-017-12939-7 29018237
Silva, D. M. Z. D. A. et al. High-throughput analysis unveils a highly shared satellite DNA library among three species of fish genus Astyanax. Sci. Rep. 7, 12726. 10.1038/s41598-017-12939-7 (2017).29018237 10.1038/s41598-017-12939-7
66. Goes CAG Revealing the satellite DNA history in Psalidodon and Astyanax characid fish by comparative satellitomics Front. Genet. 2022 13 884072 10.3389/fgene.2022.884072 35801083
Goes, C. A. G. et al. Revealing the satellite DNA history in Psalidodon and Astyanax characid fish by comparative satellitomics. Front. Genet. 13, 884072. 10.3389/fgene.2022.884072 (2022).35801083 10.3389/fgene.2022.884072
67. Kretschmer R Satellitome analysis illuminates the evolution of ZW sex chromosomes of Triportheidae fishes (Teleostei: Characiformes) Chromosoma 2022 131 29 45 10.1007/s00412-022-00768-1 35099570
Kretschmer, R. et al. Satellitome analysis illuminates the evolution of ZW sex chromosomes of Triportheidae fishes (Teleostei: Characiformes). Chromosoma 131, 29–45. 10.1007/s00412-022-00768-1 (2022).35099570 10.1007/s00412-022-00768-1
68. Utsunomia R Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome Sci. Rep. 2019 9 5856 10.1038/s41598-019-42383-8 30971780
Utsunomia, R. et al. Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome. Sci. Rep. 9, 5856. 10.1038/s41598-019-42383-8 (2019).30971780 10.1038/s41598-019-42383-8
69. Toma GA Cytogenetics meets genomics: Cytotaxonomy and genomic relationships among color variants of the Asian Arowana ( Scleropages formosus ) Int. J. Mol. Sci. 2023 24 9005 10.3390/ijms24109005 37240350
Toma, G. A. et al. Cytogenetics meets genomics: Cytotaxonomy and genomic relationships among color variants of the Asian Arowana (Scleropages formosus). Int. J. Mol. Sci. 24, 9005. 10.3390/ijms24109005 (2023).37240350 10.3390/ijms24109005
70. Marta A Dedukh D Bartoš O Majtánová Z Janko K Cytogenetic characterization of seven novel satDNA markers in two species of spined loaches ( Cobitis ) and their clonal hybrids Genes 2020 11 617 10.3390/genes11060617 32512717
Marta, A., Dedukh, D., Bartoš, O., Majtánová, Z. & Janko, K. Cytogenetic characterization of seven novel satDNA markers in two species of spined loaches (Cobitis) and their clonal hybrids. Genes 11, 617. 10.3390/genes11060617 (2020).32512717 10.3390/genes11060617
71. Bertollo LAC Born GG Dergam JA Fenocchio AS Moreira-Filho O A biodiversity approach in the neotropical Erythrinidae fish, Hoplias malabaricus . Karyotypic survey, geographic distribution of cytotypes and cytotaxonomic considerations Chromosome Res. 2000 8 603 613 10.1023/A:1009233907558 11117356
Bertollo, L. A. C., Born, G. G., Dergam, J. A., Fenocchio, A. S. & Moreira-Filho, O. A biodiversity approach in the neotropical Erythrinidae fish, Hoplias malabaricus. Karyotypic survey, geographic distribution of cytotypes and cytotaxonomic considerations. Chromosome Res. 8, 603–613. 10.1023/A:1009233907558 (2000).11117356 10.1023/A:1009233907558
72. B Cioffi M Chromosomal evolution in lower vertebrates: sex chromosomes in neotropical fishes. Genes 8, 258 (2017). Genes 2017 8 258 10.3390/genes8100258 28981468
Cioffi M. B., Yano, C. F., Sember, A. & Bertollo, L. A. C. Chromosomal evolution in lower vertebrates: sex chromosomes in neotropical fishes. Genes 8, 258 (2017).28981468 10.3390/genes8100258
73. Cioffi MB Franco W Ferreira R Bertollo LAC Tirunilai P Chromosomes as tools for discovering Biodiversity—The case of Erythrinidae fish family Recent Trends Cytogenet Studies Methodol Appl 2012 InTech 125 146
Cioffi, M. B., Franco, W., Ferreira, R. & Bertollo, L. A. C. Chromosomes as tools for discovering Biodiversity—The case of Erythrinidae fish family. In Recent Trends Cytogenet Studies Methodol Appl (ed. Tirunilai, P.) 125–146 (InTech, 2012).
74. 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 10.1038/hdy.2010.18 20216570
Cioffi, M. B. & Bertollo, L. A. C. Initial steps in XY chromosome differentiation in Hoplias malabaricus and the origin of an X1X2Y sex chromosome system in this fish group. Heredity 105, 554–561. 10.1038/hdy.2010.18 (2010).20216570 10.1038/hdy.2010.18
75. Sember A Sex chromosome evolution and genomic divergence in the fish Hoplias malabaricus (Characiformes, Erythrinidae) Front. Genet. 2018 9 71 10.3389/fgene.2018.00071 29556249
Sember, A. et al. Sex chromosome evolution and genomic divergence in the fish Hoplias malabaricus (Characiformes, Erythrinidae). Front. Genet. 9, 71 (2018).29556249 10.3389/fgene.2018.00071
76. Martins C Ferreira IA Oliveira C Foresti F Galetti PM A tandemly repetitive centromeric DNA sequence of the fish Hoplias malabaricus (Characiformes: Erythrinidae) is derived from 5S rDNA Genetica 2006 127 133 10.1007/s10709-005-2674-y 16850219
Martins, C., Ferreira, I. A., Oliveira, C., Foresti, F. & Galetti, P. M. A tandemly repetitive centromeric DNA sequence of the fish Hoplias malabaricus (Characiformes: Erythrinidae) is derived from 5S rDNA. Genetica 127, 133. 10.1007/s10709-005-2674-y (2006).16850219 10.1007/s10709-005-2674-y
77. dos Santos RZ A Long-term conserved satellite DNA that remains unexpanded in several genomes of Characiformes fish is actively transcribed Genome Biol. Evol. 2021 13 evab002 10.1093/gbe/evab002 33502491
dos Santos, R. Z. et al. A Long-term conserved satellite DNA that remains unexpanded in several genomes of Characiformes fish is actively transcribed. Genome Biol. Evol. 13, evab002. 10.1093/gbe/evab002 (2021).33502491 10.1093/gbe/evab002
78. Goes CAG The satellite DNA catalogues of two Serrasalmidae (Teleostei, Characiformes): Conservation of general satDNA features over 30 million years Genes 2022 14 91 10.3390/genes14010091 36672835
Goes, C. A. G. et al. The satellite DNA catalogues of two Serrasalmidae (Teleostei, Characiformes): Conservation of general satDNA features over 30 million years. Genes 14, 91. 10.3390/genes14010091 (2022).36672835 10.3390/genes14010091
79. da Silva MJ Gazoni T Haddad CFB Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest Front. Genet. 2023 14 1101397 10.3389/fgene.2023.1101397 37065500
da Silva, M. J., Gazoni, T. & Haddad, C. F. B. Analysis in Proceratophrys boiei genome illuminates the satellite DNA content in a frog from the Brazilian Atlantic forest. Front. Genet. 14, 1101397. 10.3389/fgene.2023.1101397 (2023).37065500 10.3389/fgene.2023.1101397
80. Sena RS Identification and characterization of satellite DNAs in two-toed sloths of the genus Choloepus (Megalonychidae, Xenarthra) Sci. Rep. 2020 10 19202 10.1038/s41598-020-76199-8 33154538
Sena, R. S. et al. Identification and characterization of satellite DNAs in two-toed sloths of the genus Choloepus (Megalonychidae, Xenarthra). Sci. Rep. 10, 19202. 10.1038/s41598-020-76199-8 (2020).33154538 10.1038/s41598-020-76199-8
81. Vozdova M Satellite DNA in Neotropical deer species Genes 2021 12 123 10.3390/genes12010123 33478071
Vozdova, M. et al. Satellite DNA in Neotropical deer species. Genes 12, 123. 10.3390/genes12010123 (2021).33478071 10.3390/genes12010123
82. Peona V Kutschera VE Blom MPK Irestedt M Suh A Satellite DNA evolution in Corvoidea inferred from short and long reads Mol. Ecol. 2023 32 1288 1305 10.1111/mec.16484 35488497
Peona, V., Kutschera, V. E., Blom, M. P. K., Irestedt, M. & Suh, A. Satellite DNA evolution in Corvoidea inferred from short and long reads. Mol. Ecol. 32, 1288–1305. 10.1111/mec.16484 (2023).35488497 10.1111/mec.16484
83. Lisachov A Rumyantsev A Prokopov D Ferguson-Smith M Trifonov V Conservation of major satellite DNAs in snake heterochromatin Animals 2023 13 334 10.3390/ani13030334 36766223
Lisachov, A., Rumyantsev, A., Prokopov, D., Ferguson-Smith, M. & Trifonov, V. Conservation of major satellite DNAs in snake heterochromatin. Animals 13, 334. 10.3390/ani13030334 (2023).36766223 10.3390/ani13030334
84. Melters DP Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution Genome Biol. 2013 14 R10 10.1186/gb-2013-14-1-r10 23363705
Melters, D. P. et al. Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution. Genome Biol. 14, R10. 10.1186/gb-2013-14-1-r10 (2013).23363705 10.1186/gb-2013-14-1-r10
85. Voleníková A Fast satellite DNA evolution in Nothobranchius annual killifishes Chromosome Res. 2023 31 33 10.1007/s10577-023-09742-8 37985497
Voleníková, A. et al. Fast satellite DNA evolution in Nothobranchius annual killifishes. Chromosome Res. 31, 33. 10.1007/s10577-023-09742-8 (2023).37985497 10.1007/s10577-023-09742-8
86. Henikoff S Ahmad K Malik HS The centromere paradox: Stable inheritance with rapidly evolving DNA Science 2001 293 1098 1102 10.1126/science.1062939 11498581
Henikoff, S., Ahmad, K. & Malik, H. S. The centromere paradox: Stable inheritance with rapidly evolving DNA. Science 293, 1098–1102. 10.1126/science.1062939 (2001).11498581 10.1126/science.1062939
87. Kitaoka M Smith OK Straight AF Heald R Molecular conflicts disrupting centromere maintenance contribute to Xenopus hybrid inviability Curr. Biol. 2022 32 3939 3951.e6 10.1016/j.cub.2022.07.037 35973429
Kitaoka, M., Smith, O. K., Straight, A. F. & Heald, R. Molecular conflicts disrupting centromere maintenance contribute to Xenopus hybrid inviability. Curr. Biol. 32, 3939-3951.e6. 10.1016/j.cub.2022.07.037 (2022).35973429 10.1016/j.cub.2022.07.037
88. Ferree PM Barbash DA Species-specific heterochromatin prevents mitotic chromosome segregation to cause hybrid lethality in Drosophila PLoS Biol. 2009 7 e1000234 10.1371/journal.pbio.1000234 19859525
Ferree, P. M. & Barbash, D. A. Species-specific heterochromatin prevents mitotic chromosome segregation to cause hybrid lethality in Drosophila. PLoS Biol. 7, e1000234. 10.1371/journal.pbio.1000234 (2009).19859525 10.1371/journal.pbio.1000234
89. O’Neill MJ O’Neill RJ Sex chromosome repeats tip the balance towards speciation Mol. Ecol. 2018 27 3783 3798 10.1111/mec.14577 29624756
O’Neill, M. J. & O’Neill, R. J. Sex chromosome repeats tip the balance towards speciation. Mol. Ecol. 27, 3783–3798. 10.1111/mec.14577 (2018).29624756 10.1111/mec.14577
90. Utsunomia R Genetic differentiation among distinct karyomorphs of the wolf fish Hoplias malabaricus species complex (Characiformes, Erythrinidae) and report of unusual hybridization with natural triploidy J. Fish Biol. 2014 85 1682 1692 10.1111/jfb.12526 25263542
Utsunomia, R. et al. Genetic differentiation among distinct karyomorphs of the wolf fish Hoplias malabaricus species complex (Characiformes, Erythrinidae) and report of unusual hybridization with natural triploidy. J. Fish Biol. 85, 1682–1692. 10.1111/jfb.12526 (2014).25263542 10.1111/jfb.12526
91. Raskina O Barber JC Nevo E Belyayev A Repetitive DNA and chromosomal rearrangements: Speciation-related events in plant genomes Cytogenet. Genome Res. 2008 120 351 357 10.1159/000121084 18504364
Raskina, O., Barber, J. C., Nevo, E. & Belyayev, A. Repetitive DNA and chromosomal rearrangements: Speciation-related events in plant genomes. Cytogenet. Genome Res. 120, 351–357. 10.1159/000121084 (2008).18504364 10.1159/000121084
92. George CM Alani E Multiple cellular mechanisms prevent chromosomal rearrangements involving repetitive DNA Crit. Rev. Biochem. Mol. Biol. 2012 47 297 313 10.3109/10409238.2012.675644 22494239
George, C. M. & Alani, E. Multiple cellular mechanisms prevent chromosomal rearrangements involving repetitive DNA. Crit. Rev. Biochem. Mol. Biol. 47, 297–313. 10.3109/10409238.2012.675644 (2012).22494239 10.3109/10409238.2012.675644
93. Chen J-M Cooper DN Férec C Kehrer-Sawatzki H Patrinos GP Genomic rearrangements in inherited disease and cancer Semin. Cancer Biol. 2010 20 222 233 10.1016/j.semcancer.2010.05.007 20541013
Chen, J.-M., Cooper, D. N., Férec, C., Kehrer-Sawatzki, H. & Patrinos, G. P. Genomic rearrangements in inherited disease and cancer. Semin. Cancer Biol. 20, 222–233. 10.1016/j.semcancer.2010.05.007 (2010).20541013 10.1016/j.semcancer.2010.05.007
94. Barra V Fachinetti D The dark side of centromeres: Types, causes and consequences of structural abnormalities implicating centromeric DNA Nat. Commun. 2018 9 4340 10.1038/s41467-018-06545-y 30337534
Barra, V. & Fachinetti, D. The dark side of centromeres: Types, causes and consequences of structural abnormalities implicating centromeric DNA. Nat. Commun. 9, 4340. 10.1038/s41467-018-06545-y (2018).30337534 10.1038/s41467-018-06545-y
95. McKinley KL Cheeseman IM The molecular basis for centromere identity and function Nat. Rev. Mol. Cell Biol. 2016 17 16 29 10.1038/nrm.2015.5 26601620
McKinley, K. L. & Cheeseman, I. M. The molecular basis for centromere identity and function. Nat. Rev. Mol. Cell Biol. 17, 16–29. 10.1038/nrm.2015.5 (2016).26601620 10.1038/nrm.2015.5
96. Furman B Sex chromosome evolution: so many exceptions to the rules Genome Biol. Evol. 2020 12 750 10.1093/gbe/evaa081 32315410
Furman, B. L. et al. Sex chromosome evolution: so many exceptions to the rules. Genome Biol. Evol. 12, 750–763 (2020).32315410 10.1093/gbe/evaa081
97. Moreira-Filho O Bertollo LAC Galetti PM Distribution of sex chromosome mechanisms in neotropical fish and description of a ZZ/ZW system in Parodon hilarii (Parodontidae) Caryologia 1993 46 115 125 10.1080/00087114.1993.10797253
Moreira-Filho, O., Bertollo, L. A. C. & Galetti, P. M. Distribution of sex chromosome mechanisms in neotropical fish and description of a ZZ/ZW system in Parodon hilarii (Parodontidae). Caryologia 46, 115–125. 10.1080/00087114.1993.10797253 (1993).10.1080/00087114.1993.10797253
98. Östergren G The mechanism of co-orientation in bivalents and multivalents: The theory of orientation by pulling Hereditas 1951 37 85 156 10.1111/j.1601-5223.1951.tb02891.x
Östergren, G. The mechanism of co-orientation in bivalents and multivalents: The theory of orientation by pulling. Hereditas 37, 85–156. 10.1111/j.1601-5223.1951.tb02891.x (1951).10.1111/j.1601-5223.1951.tb02891.x
99. Štundlová J Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci Chromosome Res. 2022 30 309 333 10.1007/s10577-022-09707-3 36208359
Štundlová, J. et al. Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci. Chromosome Res. 30, 309–333. 10.1007/s10577-022-09707-3 (2022).36208359 10.1007/s10577-022-09707-3
100. Nanda I Evolution of the degenerated Y-chromosome of the swamp guppy, Micropoecilia picta Cells 2022 11 1118 10.3390/cells11071118 35406682
Nanda, I. et al. Evolution of the degenerated Y-chromosome of the swamp guppy, Micropoecilia picta. Cells 11, 1118. 10.3390/cells11071118 (2022).35406682 10.3390/cells11071118
101. Bertollo LAC Cioffi MB Moreira-Filho O Ozouf-Costaz C Pisano E Foresti F Almeida Toledo LF Direct chromosome preparation from Freshwater Teleost Fishes Fish cytogenetic techniques (Chondrichthyans and Teleosts) 2015 CRC Press 21 26
Bertollo, L. A. C., Cioffi, M. B. & Moreira-Filho, O. Direct chromosome preparation from Freshwater Teleost Fishes. In Fish cytogenetic techniques (Chondrichthyans and Teleosts) (eds Ozouf-Costaz, C. et al.) 21–26 (CRC Press, 2015).
102. Kligerman AD Bloom SE Rapid chromosome preparations from solid tissues of fishes Can. J. Fish Aquat. Sci. 1977 34 266 269 10.1139/f77-039
Kligerman, A. D. & Bloom, S. E. Rapid chromosome preparations from solid tissues of fishes. Can. J. Fish Aquat. Sci. 34, 266–269. 10.1139/f77-039 (1977).10.1139/f77-039
103. Sambrook J Russell DW Molecular Cloning: A Laboratory Manual 2001 Cold Spring Harbor Laboratory Press
Sambrook, J. & Russell, D. W. Molecular Cloning: A Laboratory Manual (Cold Spring Harbor Laboratory Press, 2001).
104. Bolger AM Lohse M Usadel B Trimmomatic: a flexible trimmer for Illumina sequence data Bioinformatics 2014 30 2114 2120 10.1093/bioinformatics/btu170 24695404
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. 10.1093/bioinformatics/btu170 (2014).24695404 10.1093/bioinformatics/btu170
105. Schmieder R Edwards R Quality control and preprocessing of metagenomic datasets Bioinformatics 2011 27 863 864 10.1093/bioinformatics/btr026 21278185
Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. 10.1093/bioinformatics/btr026 (2011).21278185 10.1093/bioinformatics/btr026
106. Smit, A. F. A., Hubley, R. & Green, P. RepeatMasker Open (1996–2010).
107. Nascimento M PHYLOViZ 2.0: Providing scalable data integration and visualization for multiple phylogenetic inference methods Bioinformatics 2016 33 128 129 10.1093/bioinformatics/btw582 27605102
Nascimento, M. et al. PHYLOViZ 2.0: Providing scalable data integration and visualization for multiple phylogenetic inference methods. Bioinformatics 33, 128–129. 10.1093/bioinformatics/btw582 (2016).27605102 10.1093/bioinformatics/btw582
108. Langmead B Salzberg SL Fast gapped-read alignment with Bowtie 2 Nat. Methods 2012 9 357 359 10.1038/nmeth.1923 22388286
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. 10.1038/nmeth.1923 (2012).22388286 10.1038/nmeth.1923
109. Fu L Niu B Zhu Z Wu S Li W CD-HIT: Accelerated for clustering the next-generation sequencing data Bioinformatics 2012 28 3150 3152 10.1093/bioinformatics/bts565 23060610
Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. 10.1093/bioinformatics/bts565 (2012).23060610 10.1093/bioinformatics/bts565
110. Edgar RC MUSCLE: Multiple sequence alignment with high accuracy and high throughput Nucleic Acids Res. 2004 32 1792 1797 10.1093/nar/gkh340 15034147
Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. 10.1093/nar/gkh340 (2004).15034147 10.1093/nar/gkh340
111. Cioffi MB Martins C Bertollo LAC Comparative chromosome mapping of repetitive sequences. Implications for genomic evolution in the fish Hoplias malabaricus BMC Genet. 2009 10 34 10.1186/1471-2156-10-34 19583858
Cioffi, M. B., Martins, C. & Bertollo, L. A. C. Comparative chromosome mapping of repetitive sequences. Implications for genomic evolution in the fish Hoplias malabaricus. BMC Genet. 10, 34. 10.1186/1471-2156-10-34 (2009).19583858 10.1186/1471-2156-10-34
112. Yano CF Bertollo LAC de Cioffi M Liehr T B. Fish-FISH: Molecular cytogenetics in fish species Fluorescence in situ Hybridization 2017 2 Springer 429 443
Yano, C. F., Bertollo, L. A. C. & de Cioffi, M. B. Fish-FISH: Molecular cytogenetics in fish species. In Fluorescence in situ Hybridization 2nd edn (ed. Liehr, T.) 429–443 (Springer, 2017).
113. Levan A Fredga K Sandberg AA Nomenclature for centromeric position on chromosomes Hereditas 1964 52 201 220 10.1111/j.1601-5223.1964.tb01953.x
Levan, A., Fredga, K. & Sandberg, A. A. Nomenclature for centromeric position on chromosomes. Hereditas 52, 201–220. 10.1111/j.1601-5223.1964.tb01953.x (1964).10.1111/j.1601-5223.1964.tb01953.x
