Chromosomal Rearrangements and Origin of the Multiple XX/XY1Y2 Sex Chromosome System in Harttia Species (Siluriformes: Loricariidae)

The Neotropical genus Harttia comprises species with extensive chromosomal remodeling and distinct sex chromosome systems (SCSs). So far, three different SCSs with male heterogamety have been characterized in the group. In some species, the presence of the XX/XY₁Y₂ SCS is associated with a decrease in diploid numbers and several chromosomal rearrangements, although a direct relation to sex chromosome differentiation has not been shown yet. Here, we aimed to investigate the differentiation processes that have led to the establishment of the rare XX/XY₁Y₂ SCS and track its evolutionary history among other Harttia species. For that, four whole chromosome painting probes derived from chromosome 1 of H. torrenticola (HTO-1), chromosomes 9 and X of H. carvalhoi (HCA-9 and HCA-X), and chromosome X from H. intermontana (HIN-X) were applied in nine Harttia species. Homeologous chromosome blocks were located in Harttia species and demonstrated that Robertsonian (Rb) fusions originated HTO-1, HCA-9, and HCA-X chromosomes, while Rb fissions explain Y₁ and Y₂ sex chromosomes. Specifically, in H. intermontana, HCA-X, HCA-9, and the NOR-bearing chromosome demonstrated that homeologous blocks were used in the HIN-X and metacentric pair 2 origins. Consequently, diploid numbers changed between the studied species. Overall, the data also reinforce the existence of unstable genomic sites promoting chromosomal differentiation and remodeling within the genus Harttia.

Introduction

Although sex determination can be environmentally determined in some species, it is usually genetically regulated, and often associated with the presence of sex chromosomes (Furman et al., 2020). According to a widely accepted model, sex chromosomes arise from an autosomal pair due to the emergence of a sex-specific locus in one of the homologous (Bull, 1983; Charl-esworth, 2002). Over time, the ancestral homologous pair undergoes divergences in its genetic composition, recombination rate, and morphology, leading to sex chromosomes differentiation (Charlesworth et al., 2005; Bachtrog et al., 2014). Thus, sex chromosomes can be recognized according to their size and shape in a karyotype (Ghigliotti et al., 2016). However, sometimes sex chromosomes are indistinguishable concerning their gross morphology, and so defined as homomorphic ones (Chalopin et al., 2015). The most common examples of heteromorphic systems are the XX/XY, where the Y chromosome is restricted to males, and the ZZ/ZW, where the W chromosome is restricted to females (Graves, 2006). Being observed in most mammals and birds, respectively, both systems show different levels of genetic divergence (Graves, 2006; Chalopin et al., 2015). Although little is known about why and how different SCSs have evolved (Tree of Sex Consortium, 2014), the processes associated with their evolutionary origin and differentiation among vertebrates have awoken considerable interest (Devlin and Nagahama, 2002). In fishes, a high diversity of sex-determining mechanisms and SCSs with independent origins can be found (Devlin and Nagahama, 2002; Sember et al., 2021), thus making comparative evolutionary analyzes somewhat difficult. Fishes often present high plasticity concerning sex chromosomes, including none or only subtle changes between the sex pair to major chromosomal rearrangements and size differences (Devlin and Nagahama, 2002).

Among the multiple SCSs, the following types were already identified in fishes X₁X₁X₂X₂/X₁X₂Y, X₁X₁X₂X₂/X₁Y₁X₂Y₂, XX/XY₁Y₂, Z₁Z₁Z₂Z₂/Z₁Z₂W₁W₂, and ZZ/ZW₁W₂ (Kitano and Peichel, 2012). While the X₁X₁X₂X₂/X₁X₂Y system is well-represented among several fish families, the XX/XY₁Y₂ system is found only in a few (Kitano and Peichel, 2012; Sember et al., 2021). In contrast to simple SCSs, where repetitive DNAs play an essential role in sex chromosome differentiation (Yano et al., 2014; Schemberger et al., 2019), multiple SCSs appear forced by divergent evolutionary trends. It appears that chromosomal rearrangements are more relevant to the evolutionary process of multiple SCSs than the accumulation of repetitive sequences (Almeida et al., 2015). For this reason, molecular cytogenetic procedures based on fluorescence in situ hybridization (FISH), e.g., using whole chromosome painting (WCP) probes, has been successfully applied in different fish groups, providing new insights into the differentiation of sex chromosomes, especially for multiple ones (Cioffi et al., 2011; Blanco et al., 2014; Oliveira et al., 2018; Moraes et al., 2019).

Harttia is a Neotropical fish group comprising species with distinct diploid numbers and karyotypic variations emerged by extensive evolutionary conserved chromosomal rearrangements (Blanco et al., 2013, 2014, 2017; Deon et al., 2020; Sassi et al., 2020; Sassi et al., 2021). The chromosomal number ranges from 2n = 52 to 62, including B chromosomes and different SCSs (Blanco et al., 2012, 2017; Deon et al., 2020; Sassi et al., 2020; Sassi et al., 2021). In phylogenetic reconstructions, three distinct clades were proposed for the genus, thus reinforcing the extensive diversification experienced by the lineage; also, it is grouping the species according to their South American distribution: (I) from the Guyana shield rivers; (II) from the northern Brazilian rivers; and (III) from the Brazilian south/southeast rivers (Londoño-Burbano and Reis, 2021). Three SCSs were detected so far: (1) the X₁X₁X₂X₂/X₁X₂Y system, present in H. punctata, H. duriventris, and H. villasboas, and (2) a proto/neo-XX/XY system in H. rondoni, both belonging to clade II, and (3) the XX/XY₁Y₂ system in H. carvalhoi, H. intermontana, and Harttia sp.1, species which belong to clade III (Centofante et al., 2006; Blanco et al., 2017; Sassi et al., 2020; Sassi et al., 2021; Deon et al., submitted). Chromosomal data compared to a phylogenetic framework indicate that ancestral karyotype with 2n = 58 chromosomes and without a differentiated SCS could represent a plesiomorphic condition for clade III (Deon et al., submitted). Belonging to the same clade III, the species H. torrenticola has a karyotype composed by 2n = 56, undifferentiated sex chromosomes (Blanco et al., 2013), and a large metacentric pair being morphologically similar to the X chromosome of H. carvalhoi. WCP-FISH experiments using X₁ and X₂ probes derived from H. punctata, confirmed that chromosomes that gave rise to the X₁X₂Y and the XY₁Y₂ systems are evolutionary independent (Deon et al., submitted).

Here, we aimed to investigate the differentiation processes that have led to the establishment of the rare XX/XY₁Y₂ SCS and to track its evolutionary history among other Harttia species. For that, we performed a WCP-FISH investigation using four distinct sex chromosome-specific probes hybridized in several species. The results allowed us to identify the main rearrangements involved in the origin of this unique SCS. Besides, the data provide new insights into the origin and evolution of such a rare XY-derived SCS, consequently increasing our knowledge about the evolution of vertebrate sex chromosomes.

Materials and Methods

Individuals and Chromosome Preparation

Representatives of Harttia species analyzed in this study are summarized in Table 1. Specimens were collected with the authorization of the Chico Mendes Institute for Biodiversity Conservation (ICMBIO), System of Authorization and Information about Biodiversity (SISBIO-Licenses No. 10538-3 and 15117-2), and National System of Genetic Resource Management and Associated Traditional Knowledge (SISGEN-A96FF09), Brazil. Species were identified based on their morphological features by Dr. Oswaldo Oyakawa curator of the fish collection of the Museu de Zoologia da Universidade de São Paulo (MZUSP), Brazil.

TABLE 1

TABLE 1. Collection sites of the studied species, diploid chromosome number (2n), and sample size (N).

Mitotic chromosomes were obtained from kidney cells, according to Bertollo et al. (2015). All procedures agreed with the Ethics Committee of Animal Usage of the Universidade Federal de São Carlos (Process number CEUA 1853260315), Brazil.

Chromosome Microdissection, Probe Preparation, and Labeling

Fifteen copies of each target chromosome were isolated by glass-needle-based microdissection, and obtained DNA was amplified by oligonucleotide primed-polymerase chain reaction (DOP-PCR) as described in Yang et al. (2009). Chromosomes were chosen based on their morphology - bi-armed chromosomes that were suspected to be originated from Robertsonian fusions were targeted: the largest metacentric (HCA-X), and the largest submetacentric (HCA-9) from H. carvalhoi; the largest metacentric (HIN-X) from H. intermontana, and the largest metacentric (HTO-1) from H. torrenticola (Figure 1). Probes were labeled with Spectrum Orange-dUTP or Spectrum Green-dUTP (Vysis, Downers Grove, United States) in a secondary DOP-PCR, using 1 μL of the primarily amplified product as a template DNA (Yang and Graphodatsky, 2009).

FIGURE 1

FIGURE 1. Schematic representation of the microdissected chromosomes used for probe construction in comparative WCP-FISH analyses. 1) the X chromosome of H. carvalhoi (HCA-X, in yellow); 2) chromosome 9 of H. carvalhoi (HCA-9, in light pink); 3) the X chromosome of H. intermontana (HIN-X, in purple); 4) chromosome 1 of H. torrenticola (HTO-1, in green).

Fluorescence in situ Hybridization

Slides were prepared and pre-treated according to Yang et al. (2009) and denatured in 70% formamide/2x SSC for 3 min at 72°C. For each slide, 12 μL of hybridization solution (containing 0.2 μg of each labeled probe, 50% formamide, 2×SSC, 10% dextran sulfate, and 5 μg of salmon sperm DNA) was denatured for 10 min at 75°C and incubated to pre-hybridize for 1 h at 37°C. To block the hybridization of high-copy repeat sequences, 20 μg of C₀t-1 DNA, directly prepared from H. carvalhoi, H. torrenticola, and H. intermontana male genomes were used, according to Zwick et al. (1997). Hybridization took place for 48 h at 37°C in a moist chamber. Post-hybridization washes were performed in 1×SSC for 5 min at 65°C, and 5 min in 4×SSC/Tween at room temperature. Finally, the slides were counterstained with 4’,6-diamidino-2-phenylindole (DAPI) in Vectashield mounting medium (Vector, Burlingame, CA, United States).

Image Analyses and Processing

Metaphase plates were captured using an Olympus BX50 light microscope (Olympus Corporation, Ishikawa, Japan) coupled with a CoolSNAP camera. The images were processed using Image-Pro Plus 4.1 software (Media Cybernetics, Silver Spring, MD, United States). The figures were edited and organized using Adobe Photoshop CC 2020 (San Jose, CA, United States) software.

Results

Results obtained by HCA-X and HCA-9 probes are summarized in Figure 2 and Table 2. In H. carvalhoi (52♀/53♂ - XX/XY₁Y₂), the HCA-X probe successfully identified their X chromosomes in females and the X, Y_1, and Y₂ chromosomes in males. Small centromeric signals in both acrocentric pairs 23 and 24 were also evidenced. In agreement, the HCA-9 probe correctly recognized the submetacentric pair 9 (Figures 2A,B). Similarly, in Harttia sp. 1 (56♀/57♂ - XX/XY₁Y₂), the HCA-X probe detected the X chromosome pair in females and the X, Y₁, and Y₂ chromosomes in males, besides small centromeric signals in both 20 and 24 acrocentric pairs (Figures 2C,D). The HCA-9 hybridized to 21 and 26 acrocentric pairs (Figures 2C,D). In females of H. intermontana (52♀/53♂ - XX/XY₁Y₂), the HCA-X probe stained the long (q) arms of the chromosomes X and 2 (Figure 2E). In males, this probe gave signals on Xq, the Y₂ chromosome, and the 2q (Figure 2F), as well as in the centromeric region of the pair 24 in both males and females (Figures 2E,F). The HCA-9 probe detected the short (p) arms of the X chromosome and the 20q distal region in females, and the Xp arms, the Y₁ chromosome, and the 20q distal region in males (Figures 2E,F). In H. punctata (58♀/57♂ - X₁X₁X₂X₂/X₁X₂Y), the HCA-X probe hybridized on the submetacentric pairs 9 and 11, while the HCA-9 probe showed signals on the metacentric pair 8 and subtelocentric pair 19 (Figures 2G,H).

FIGURE 2

FIGURE 2. Whole chromosome painting by FISH using HCA-X (green) and HCA-9 (red) probes among Harttia species that possess SCSs. The numbers of the labeled chromosome pairs are highlighted in the images. In (A, B) metaphases of H. carvalhoi female and male, respectively; (C, D) metaphases of Harttia sp. 1 female and male, respectively; (E, F) metaphases of H. intermontana female and male, respectively; and (G, H) metaphases of H. punctata female and male, respectively.Bar = 5 µm.

TABLE 2

TABLE 2. Main localization of the WCP probes in Harttia species. Some small signals were not considered. SCSs means sex chromosome systems.

The HTO-1 probe, derived from H. torrenticola, showed the same results obtained with the HCA-X probe when tested on those species with an identical large metacentric pair (Supplementary Figure S1). In contrast, the HIN-X probe, from H. intermontana, showed different results than those obtained applying HCA-X and HTO-1 probes (Figure 3). In H. carvalhoi, HIN-X hybridized on the Xq and 9q arms in females (Figure 3A), and on the Xq arms, Y₂ chromosome, and 9q arms in males (Figure 3B). In Harttia sp. 1, HIN-X labeled the Xq arms and the acrocentric 21 pair in females (Figure 3C), and these same chromosomes, as well as the Y₂ chromosome, in males (Figure 3D). In H. torrenticola, HIN-X stained the 1q arms and the acrocentric 23 pair (Figures 3E,F).

FIGURE 3

FIGURE 3. Whole chromosome painting by FISH using the HIN-X probe (red) among Harttia species that shared a large metacentric pair. The numbers of the labeled chromosome pairs are highlighted in the images. In (A, B) metaphases of H. carvalhoi female and male, respectively; (C, D) metaphases of Harttia sp. 1 female and male, respectively; and (E, F) metaphases of H. torrenticola female and male, respectively. Bar = 5 µm.

Harttia species without the heteromorphic sex chromosomes were also used as targets for the comparative WCP-FISH using the HCA-X and HCA-9 probes (Figure 4; Table 2). In H. kronei (58♀♂) HCA-X hybridized in the subtelocentric pairs 17 and 19, while the HCA-9 hybridized in chromosome pairs 8 and 13 (Figure 4A). In H. gracilis (58♀♂), HCA-X marked the submetacentric 11 and the subtelocentric 22, besides centromeric signals in the acrocentric pairs 26 and 28 (Figure 4B), and the HCA-9 probe was detected in chromosomes 10 and 21 (Figure 4B). In H. longipinna (58♀♂), HCA-X hybridized in the subtelocentric pairs 15 and 17, besides the centromeric region of the acrocentric pairs 23 and 25, and HCA-9 hybridized in the metacentric 7, and submetacentric 10 (Figure 4C). In H. loricariformis (56♀♂), HCA-X was detected in the submetacentric pair 9, subtelocentric 20, and the centromeric region of the chromosome 25 (Figure 4D), while the HCA-9 probe hybridized in chromosome pairs 8 and 11 (Figure 4D). Finally, in H. torrenticola (56♀♂), HCA-X hybridized in chromosome 1 and the centromeric region of the acrocentric pairs 22 and 25, while the HCA-9 probe presents signals of hybridization in pairs 8 and 23 (Figure 4E).

FIGURE 4

FIGURE 4. Whole chromosome paint by FISH using HCA-X (green) and HCA-9 probes (red) in Harttia species without heteromorphic SCS. The numbers of the labeled chromosome pairs are highlighted in the images. Metaphases of H. kronei (A), H. gracilis (B), H. longipinna (C), H. loricariformis (D), and H. torrenticola (E).Bar = 5 µm.

Discussion

In Harttia species, diploid numbers range from 52 to 62 chromosomes (Deon et al., 2020; Sassi et al., 2020; Sassi et al., 2021). Data from phylogeny reconstructions indicate that 58 chromosomes and no large biarmed chromosomes could correspond to a plesiomorphic karyotype condition for species distributed on south and southeast Brazilian drainages—the clade III (Deon et al., submitted). These chromosomal features (Figure 5) include the absence of morphologically differentiated sex chromosomes and a single location of the 5S and 45S rDNA sites in medium-sized bi-armed chromosomes (Deon et al., 2020; Deon et al., submitted). Here, the ancestral reconstructions of the Harttia karyotype, using both HCA-X and HCA-9 probes, demonstrated that two chromosome pairs were probably related to the origin of the H. carvalhoi chromosomes X and 9 (Figure 6). Thus, in situ localizations also reaffirm the role of Robertsonian fusions as the main rearrangements responsible for reducing the diploid number in H. carvalhoi. These homeologous chromosome pairs (unfused chromosomes) are shared by H. kronei, H. loricariformis, H. longipinna, and H. gracilis (Figure 6). As common features, H. kronei, H. longipinna, and H. gracilis kept 2n = 58 chromosomes and the absence of morphologically differentiated SCSs (Blanco et al., 2017), with chromosomal diversification events mainly occurring by repositioning of the rDNA sites in their karyotypes (Deon et al., 2020; Deon et al., submitted). Although H. loricariformis decreased the diploid number to 2n = 56, this species shares the homologous chromosome pairs to HCA-X and HCA-9 as highlighted in H. kronei by WCP-FISH. The presence of interstitial telomeric sites in a subtelocentric chromosome of H. loricariformis karyotype suggests an origin by Robertsonian fusions (Blanco et al., 2017). The current data thus support the hypothesis on the occurrence of a chromosomal fusion event in H. loricariformis karyotype, and that this corresponds to an independent evolutionary event being not associated with the chromosomes X and 9 of H. carvalhoi.

FIGURE 5

FIGURE 5. Schematic representation of the phylogenetic relationships among Harttia species adapted from the Londoño-Burbano and Reis (2021). Harttia sp.1 and H. intermontana are not being represented since they were not included in such previous analysis.

FIGURE 6

FIGURE 6. Schematic representation summarizing the distribution of the WCP probes obtained in this study: HCA-X (yellow), HCA-9 (light pink), HIN-X (purple) and HTO-1 (green) in Harttia species (A) without the largest metacentric chromosome pair -H. kronei, H. gracilis, H. longipinna, H. loricariformis, and H. punctata- and (B) with the largest metacentric pair in the karyotype -H. carvalhoi, H. torrenticola, H. intermontana and Harttia sp. 1. The highlighted boxes show the male condition and the different composition of the Y₁ and Y₂ chromosomes in H. carvalhoi, H. intermontana and Harttia sp. 1. Note the overlapping of the HCA-9 and HIN-X probes (light pink and purple), and the HCA-X and HTO-1 (yellow and green) probes.

Data also showed that the chromosomal rearrangements that led to the XX/XY₁Y₂ SCS were triggered within the branch with H. torrenticola (Figure 5). The phylogenetic branch grouping H. carvalhoi and H. torrenticola (Covain et al., 2016; Londoño-Burbano and Reis, 2021; Figure 5) was diversified by Robertsonian fusions, initially giving rise to a large metacentric pair, like that found in the H. torrenticola karyotype. Indeed, the large homeologous chromosome regions shared between the chromosomes 1 of H. torrenticola (HTO-1) and the X chromosome of H. carvalhoi and Harttia sp. 1, corroborate that a single evolutionary event of chromosomal fusion would have generated the large metacentric pair in these species. Although H. torrenticola does not show sex chromosome heteromorphism related to the metacentric chromosome 1, in H. carvalhoi and Harttia sp. 1, this chromosome corresponds to the X sex chromosome, with additional rearrangements triggering the origin of the Y₁ and Y₂ chromosomes. According to former suggestions (Blanco et al., 2017; Deon et al., 2020) centric fission on the largest metacentric formed the Y₁ and Y₂ chromosomes in Harttia sp. 1, which are also shared by H. carvalhoi. The Y₁ and Y₂ positive hybridizations using the HCA-X probe reiterate that centric fission is the main rearrangement related to the origin of the multiple XX/XY₁Y₂ SCS of H. carvalhoi and Harttia sp. 1.

However, different from Harttia sp. 1 (2n = 56 in females and 57 in males), H. carvalhoi diversified its karyotype by other chromosomal fusions, reducing the diploid number to 2n = 52 in females and 2n = 53 in males. Based on the HCA-9 WCP-FISH experiments a chromosome fusion between the subtelocentric pairs 21 and 26, like those found in Harttia sp. 1, triggered the origin of pair 9 of H. carvalhoi (Figures 5, 6). In fact, despite some morphological alterations, this chromosome pair is represented by two other, homologous pairs in H. kronei, H. loricariformis, H. longipinna, H. gracilis, and H. torrenticola. Based on conventional cytogenetic studies, Deon et al. (2020) proposed the same origin of the XX/XY₁Y₂ system in H. carvalhoi, Harttia sp. 1, and H. intermontana. However, the use of the HCA-X, HCA-9, and HIN-X probes enabled now to evidence that additional rearrangements are associated with the XX/XY₁Y₂ system of H. intermontana. The X chromosome of this species comprises the 9q arms and one arm of the X chromosome of H. carvalhoi, indicating a reciprocal translocation between these two chromosome pairs in its origin (Figure 6). After that, centric fission in one of the X chromosomes, followed by a pericentric inversion in one of the resulted elements, generated the Y₁ and Y₂ chromosomes in males (Figure 7). It is relevant to notice that the Y₁ chromosome of H. intermontana is derived from chromosome 9 of H. carvalhoi, thus different from the Y₁ chromosome of H. carvalhoi and Harttia sp. 1. In the same way, the metacentric pair 2 of H. intermontana was originated from species-specific chromosomal rearrangements, implying a translocation between the acrocentric chromosome bearing the 45S site and one X chromosome arm. Indeed, the chromosome pair of H. intermontana bears the 45S rDNA locus (Deon et al., 2020), a site prone to breaks in Harttia karyotypes, leading to extensive chromosomal remodeling events (Deon et al., 2020; Deon et al., submitted).

FIGURE 7

FIGURE 7. Schematic model representing the rearrangements occurred from H. carvalhoi to H. intermontana, and the evolvement of HCA-X (yellow) and HCA-9 probes (light pink) on the origin of the XX/XY₁Y₂ SCS. Considering H. intermontana, a derived species from H. carvalhoi, a centric fusion between the long arms of the chromosome 1 (X) and the long arms of the chromosome 9 of H. carvalhoi, lead the origin of the X chromosome in H. intermontana. An additional fusion with part of the X chromosome and part of the chromosome 23 (bearing the 18S rDNA site) gave rise to the second-largest metacentric chromosome pair in H. intermontana (pair 2). It is worth mentioning that the different origin of the X chromosome directly reflects on the genomic composition of the Y chromosomes: while the Y₁ chromosome corresponds to the long arms of the chromosome 9, the Y2 chromosome corresponds to a part of the X chromosome of H. carvalhoi.

According to molecular-phylogenetic reconstructions, H. punctata - 2n = 58♀/57♂, X₁X₁/X₂X₂/X₁X₂Y (Blanco et al., 2014) - belongs to Harttia’s clade II (Covain et al., 2016; Londoño-Burbano and Reis, 2021), and the WCP results here obtained evidenced a similar hybridization condition to those found in H. kronei, i.e., the HTO-1, HCA-X, HCA-9, and HIN-X chromosomes were not related to the karyotype diversification of H. punctata, highlighting a probable plesiomorphic condition.

It was demonstrated that sex chromosomes could emerge independently and follow distinct differentiation patterns, even among closely related species (Cioffi et al., 2013). Our WCP-FISH data also indicated independent origins for the X₁X₂Y and XY₁Y₂ SCSs of Harttia lineage, as previously proposed (Deon et al., 2020; Sassi et al., 2020). The X chromosome of the XX/XY₁Y₂ system originated by fusion of two autosome pairs, leading to the largest metacentric in the karyotype. This fusion could set up a putative homomorphic XX/XY SCS, with subsequent centric fission originating the Y₁ and Y₂ chromosomes, as proposed by Blanco et al. (2013). Thus, a set of diverse chromosomal rearrangements probably triggered the differentiation of the same or different SCSs within the Harttia lineage, suggesting that sex chromosome turnover may play an important role in the speciation processes of this group.

Evolutionarily conserved breakpoint regions (ECBRs), inside or adjacent to rDNA clusters, were proposed to occur in some Loricariidae lineages, leading to extensive chromosomal remodeling (Barros et al., 2017; Glugoski et al., 2018; Deon et al., 2020; Deon et al., submitted). In the Harttia clade III from the south/southeast Brazilian region, several rearrangements adjacent to the rDNAs sites have been extensively reused in the chromosomal diversification (Deon et al., 2020; Deon et al., submitted), including the origin of the X₁X₁X₂X₂/X₁X₂Y SCS in Harttia clade II (Deon et al., submitted). In contrast, as rDNAs were not involved in the origin of the XX/XY₁Y₂ system, other unstable sites likely occur in the genomes of the species from Harttia clade III, as well.

In addition, some hybridization signals besides those indicating the discussed main rearrangements were also detected in small chromosomal regions. They correspond to a highly rearranged repetitive DNA unit shared among Harttia species. In H. carvalhoi, H. gracilis, H. longipinna, H. torrenticola, and Harttia sp. 1, they localize close to the nucleolar organizer region (NOR) and in a small acrocentric pair. In H. loricariformis and H. intermontana only the signal close to the NOR site and in the small acrocentric chromosome were detected, respectively. The mechanisms responsible for the instability of Harttia genome are not fully understood (Deon et al., 2020; Deon et al., submitted). However, repetitive DNA clusters scattered at some genome locations are likely candidates for chromosomal breaks and rearrangements. Cytogenetic data indicate that these sequences are reused in several chromosome rearrangements, including the Robertsonian ones responsible for the origin of the SCSs and a 2n decrease in Harttia.

Interstitial telomeric sites (ITS) are common features in some Harttia genomes (Blanco et al., 2017; Deon et al., 2020). H. carvalhoi and H. torrenticola, for example, present an ITS in the large metacentric chromosome (Blanco et al., 2017), indicating its origin by Robertsonian fusion. In contrast, this ITS was lost during the chromosomal evolution of Harttia sp. 1. It is known that ITS are hotspots for breakages (Slijepcevic et al., 1997) and that telomeric DNA damages can be irreparable, causing persistent DNA-damage-response activation (Fumagalli et al., 2012), or remaining as fragile sites (Sfeir et al., 2009). According to Slijepcevic (2016), both ITS and terminal telomeric sequences are naturally prone to breakage, leading to chromosome plasticity. Therefore, the rearrangements observed in the X and 2 chromosomes of H. intermontana may have been triggered by the instability generated by the ITS in the X chromosome of H. carvalhoi.

Conclusion

Data obtained by WCP-FISH allowed to highlight small pieces of the complex chromosomal evolution that has taken place in Harttia species, with a particular emphasis on the origin of a rare multiple SCS and diploid number decrease. We demonstrated the existence of unstable genomic sites promoting chromosomal differentiation and remodeling, where homeologous chromosome blocks were identified after WCP experiments. Besides, we highlighted the distinct Robertsonian fusions and fissions that were involved in the origin the sex chromosomes. In this context, the genus Harttia has proved to be an excellent model for the study of evolution of sexual chromosome systems among Neotropical fish species. Next steps now will include a fine-scale analysis of the genetic content of the sex chromosomes in this group aiming to discover novel sex-determining genes, which is an inevitable next step towards fully understating this puzzling scenario.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by All procedures agreed with the Ethics Committee of Animal Usage of the Universidade Federal de São Carlos (Process number CEUA 1853260315), Brazil.

Author Contributions

GD, LG, FS, MC, and MV carried out the analysis and drafted the manuscript. TH, VN and AA-R helped in the analysis, drafted, and revised the manuscript. TL, LB, OM-F, MC and MV drafted and revised the manuscript. All authors read and approved the final version of the manuscript.

Funding

MC was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Proc. No. 302449/2018-3) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) (Proc. No. 2020/11772-8). FS was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) (Proc. No. 2020/002681-9). This work was financed in part by CAPES Finance code 001.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer MN declared a past co-authorship with the authors TH, LACB, TL, MDBC to the handling editor.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors are grateful to Dr. Osvaldo Oyakawa for the contributions in the animal´s identification.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.877522/full#supplementary-material

Supplementary Figure S1 | Hybridizations using the HTO-1 probe (red) in species that present a large metacentric chromosome pair in the karyotype: (A) H. carvalhoi female, (B) H. carvalhoi male, (C) Harttia sp.1 female, (D) Harttia sp. 1 male, (E) H. intermontana female and (F) H. intermontana male. Bar = 5 µm.

References

Almeida, J. S., Migues, V. H., Diniz, D., and Affonso, P. R. A. M. (2015). A Unique Sex Chromosome System in the Knifefish Gymnotus Bahianus with Inferences about Chromosomal Evolution of Gymnotidae. J. Hered. 106, 177–183. doi:10.1093/jhered/esu087

CrossRef Full Text | Google Scholar

Bachtrog, D., Mank, J. E., Peichel, C. L., Kirkpatrick, M., Otto, S. P., Ashman, T.-L., et al. (2014). Sex Determination: Why So many Ways of Doing it? Plos Biol. 12 (7), e1001899. doi:10.1371/journal.pbio.1001899

PubMed Abstract | CrossRef Full Text | Google Scholar

Barros, A. V., Wolski, M. A. V., Nogaroto, V., Almeida, M. C., Moreira-Filho, O., and Vicari, M. R. (2017). Fragile Sites, Dysfunctional Telomere and Chromosome Fusions: what Is 5S rDNA Role? Gene 608, 20–27. doi:10.1016/j.gene.2017.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertollo, L. A. C., Cioffi, M., and Moreira-Filho, O. (2015). “Direct Chromosome Preparation from Freshwater Teleost Fishes,” in Fish Cytogenetic Techniques (Chondrichthyans and Teleosts). Editors C. Ozouf-Costaz, E. Pisano, F. Foresti, and L. F. Almeida Toledo (CRC Press), 21–26.

Google Scholar

Blanco, D. R., Vicari, M. R., Artoni, R. F., Traldi, J. B., and Moreira-Filho, O. (2012). Chromosomal Characterization of Armored CatfishHarttialongipinna(Siluriformes, Loricariidae): First Report of B Chromosomes in the Genus. Zoolog. Sci. 29 (9), 604–609. doi:10.2108/zsj.29.604

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanco, D. R., Vicari, M. R., Lui, R. L., Artoni, R. F., de Almeida, M. C., Traldi, J. B., et al. (2014). Origin of the X1X1X2X2/X1X2Y Sex Chromosome System of Harttia Punctata (Siluriformes, Loricariidae) Inferred from Chromosome Painting and FISH with Ribosomal DNA Markers. Genetica 142, 119–126. doi:10.1007/s10709-014-9759-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanco, D. R., Vicari, M. R., Lui, R. L., Bertollo, L. A. C., Traldi, J. B., and Moreira-Filho, O. (2013). The Role of the Robertsonian Rearrangements in the Origin of the XX/XY1Y2 Sex Chromosome System and in the Chromosomal Differentiation in Harttia Species (Siluriformes, Loricariidae). Rev. Fish. Biol. Fish. 23 (1), 127–134. doi:10.1007/s11160-012-9283-5

CrossRef Full Text | Google Scholar

Blanco, D. R., Vicari, M. R., Lui, R. L., Traldi, J. B., Bueno, V., Martinez, J. d. F., et al. (2017). Karyotype Diversity and Evolutionary Trends in Armored Catfish Species of the Genus Harttia (Siluriformes: Loricariidae). Zebrafish 14, 169–176. doi:10.1089/zeb.2016.1377

PubMed Abstract | CrossRef Full Text | Google Scholar

Bull, J. J. (1983). Evolution of Sex Determining Mechanisms (Vol. 1). Benjamin-Cummings Publishing Company.

Google Scholar

Centofante, L., Bertollo, L. A., and Moreira-Filho, O. (2006). Cytogenetic Characterization and Description of an XX/XY1Y2 Sex Chromosome System in Catfish Harttia carvalhoi (Siluriformes, Loricariidae). Cytogenet. Genome Res. 112 (3-4), 320–324. doi:10.1159/000089887

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalopin, D., Volff, J.-N., Galiana, D., Anderson, J. L., and Schartl, M. (2015). Transposable Elements and Early Evolution of Sex Chromosomes in Fish. Chromosome Res. 23 (3), 545–560. doi:10.1007/s10577-015-9490-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Charl-esworth, B. (2002). “The Evolution of Chromosomal Sex Determination,” in The Genetics and Biology of Sex Determination. Editor R. Short, 578–585.

Google Scholar

Charlesworth, D., Charlesworth, B., and Marais, G. (2005). Steps in the Evolution of Heteromorphic Sex Chromosomes. Heredity 95 (2), 118–128. doi:10.1038/sj.hdy.6800697

PubMed Abstract | CrossRef Full Text | Google Scholar

Cioffi, M. B., Liehr, T., Trifonov, V., Molina, W. F., and Bertollo, L. A. (2013). Independent Sex Chromosome Evolution in Lower Vertebrates: a Molecular Cytogenetic Overview in the Erythrinidae Fish Family. Cytogenet. Genome Res. 141 (2-3), 186–194. doi:10.1159/000354039

PubMed Abstract | CrossRef Full Text | Google Scholar

Cioffi, M. B., Kejnovsky, E., and Bertollo, L. A. C. (2011). The Chromosomal Distribution of Microsatellite Repeats in the Genome of the Wolf Fish Hoplias malabaricus, Focusing on the Sex Chromosomes. Cytogenet. Genome Res. 132 (4), 289–296. doi:10.1159/000322058

PubMed Abstract | CrossRef Full Text | Google Scholar

Covain, R., Fisch-Muller, S., Oliveira, C., Mol, J. H., Montoya-Burgos, J. I., and Dray, S. (2016). Molecular Phylogeny of the Highly Diversified Catfish Subfamily Loricariinae (Siluriformes, Loricariidae) Reveals Incongruences with Morphological Classification. Mol. Phylogenet. Evol. 94, 492–517. doi:10.1016/j.ympev.2015.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

de Moraes, R. L. R., Sember, A., Bertollo, L. A. C., De Oliveira, E. A., Ráb, P., Hatanaka, T., et al. (2019). Comparative Cytogenetics and Neo-Y Formation in Small-Sized Fish Species of the Genus Pyrrhulina (Characiformes, Lebiasinidae). Front. Genet. 10, 678. doi:10.3389/fgene.2019.00678

PubMed Abstract | CrossRef Full Text | Google Scholar

de Oliveira, E. A., Sember, A., Bertollo, L. A. C., Yano, C. F., Ezaz, T., Moreira-Filho, O., et al. (2018). Tracking the Evolutionary Pathway of Sex Chromosomes Among Fishes: Characterizing the Unique XX/XY1Y2 System in Hoplias malabaricus (Teleostei, Characiformes). Chromosoma 127 (1), 115–128. doi:10.1007/s00412-017-0648-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Deon, G. A., Glugoski, L., Vicari, M. R., Nogaroto, V., Sassi, F. D. M. C., Cioffi, M. D. B., et al. (2020). Highly Rearranged Karyotypes and Multiple Sex Chromosome Systems in Armored Catfishes from the Genus Harttia (Teleostei, Siluriformes). Genes 11, 1366. doi:10.3390/genes11111366

PubMed Abstract | CrossRef Full Text | Google Scholar

Deon, G. A., Glugoski, L., Hatanaka, T., Sassi, F. M. C., Nogaroto, V., Bertollo, L. A. C., et al. (Submitted). Evolutionary Conserved Breakpoints Close to rDNA Regions Promote Remodeling of Homologous Chromosome Blocks Among Harttia Species (Siluriformes: Loricariidae). Genetics and Molecular Biology.

Google Scholar

Devlin, R. H., and Nagahama, Y. (2002). Sex Determination and Sex Differentiation in Fish: an Overview of Genetic, Physiological, and Environmental Influences. Aquaculture 208 (3-4), 191–364. doi:10.1016/S0044-8486(02)00057-1

CrossRef Full Text | Google Scholar

Fumagalli, M., Rossiello, F., Clerici, M., Barozzi, S., Cittaro, D., Kaplunov, J. M., et al. (2012). Telomeric DNA Damage Is Irreparable and Causes Persistent DNA-Damage-Response Activation. Nat. Cel Biol 14 (4), 355–365. doi:10.1038/ncb2466

CrossRef Full Text | Google Scholar

Furman, B. L. S., Metzger, D. C. H., Darolti, I., Wright, A. E., Sandkam, B. A., Almeida, P., et al. (2020). Sex Chromosome Evolution: So many Exceptions to the Rules. Genome Biol. Evol. 12 (6), 750–763. doi:10.1093/gbe/evaa081

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghigliotti, L., Cheng, C.-H. C., and Pisano, E. (2016). Sex Determination in Antarctic Notothenioid Fish: Chromosomal Clues and Evolutionary Hypotheses. Polar Biol. 39 (1), 11–22. doi:10.1007/s00300-014-1601-z

CrossRef Full Text | Google Scholar

Glugoski, L., Giuliano-Caetano, L., Moreira-Filho, O., Vicari, M. R., and Nogaroto, V. (2018). Co-located hAT Transposable Element and 5S rDNA in an Interstitial Telomeric Sequence Suggest the Formation of Robertsonian Fusion in Armored Catfish. Gene 650, 49–54. doi:10.1016/j.gene.2018.01.099

PubMed Abstract | CrossRef Full Text | Google Scholar

Graves, J. A. M. (2006). Sex Chromosome Specialization and Degeneration in Mammals. Cell 124 (5), 901–914. doi:10.1016/j.cell.2006.02.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kitano, J., and Peichel, C. L. (2012). Turnover of Sex Chromosomes and Speciation in Fishes. Environ. Biol. Fish. 94 (3), 549–558. doi:10.1007/s10641-011-9853-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Londoño-Burbano, A., and Reis, R. E. (2021). A Combined Molecular and Morphological Phylogeny of the Loricariinae (Siluriformes: Loricariidae), with Emphasis on the Harttiini and Farlowellini. PloS one 16 (3), e0247747. doi:10.1371/journal.pone.0247747

PubMed Abstract | CrossRef Full Text | Google Scholar

Sassi, F. d. M. C., Deon, G. A., Moreira-Filho, O., Vicari, M. R., Bertollo, L. A. C., Liehr, T., et al. (2020). Multiple Sex Chromosomes and Evolutionary Relationships in Amazonian Catfishes: the Outstanding Model of the Genus Harttia (Siluriformes: Loricariidae). Genes 11 (10), 1179. doi:10.3390/genes11101179

PubMed Abstract | CrossRef Full Text | Google Scholar

Sassi, F. d. M. C., Moreira-Filho, O., Deon, G. A., Sember, A., Bertollo, L. A. C., Liehr, T., et al. (2021). Adding New Pieces to the Puzzle of Karyotype Evolution in Harttia (Siluriformes, Loricariidae): Investigation of Amazonian Species. Biology 10 (9), 922. doi:10.3390/biology10090922

PubMed Abstract | CrossRef Full Text | Google Scholar

Schemberger, M. O., Nascimento, V. D., Coan, R., Ramos, É., Nogaroto, V., Ziemniczak, K., et al. (2019). DNA Transposon Invasion and Microsatellite Accumulation Guide W Chromosome Differentiation in a Neotropical Fish Genome. Chromosoma 128, 547–560. doi:10.1007/s00412-019-00721-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sember, A., Nguyen, P., Perez, M. F., Altmanová, M., Ráb, P., and Cioffi, M. D. B. (2021). Multiple Sex Chromosomes in Teleost Fishes from a Cytogenetic Perspective: State of the Art and Future Challenges. Phil. Trans. R. Soc. B 376 (1833), 20200098. doi:10.1098/rstb.2020.0098

PubMed Abstract | CrossRef Full Text | Google Scholar

Sfeir, A., Kosiyatrakul, S. T., Hockemeyer, D., MacRae, S. L., Karlseder, J., Schildkraut, C. L., et al. (2009). Mammalian Telomeres Resemble Fragile Sites and Require TRF1 for Efficient Replication. Cell 138 (1), 90–103. doi:10.1016/j.cell.2009.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Slijepcevic, P., Hande, M. P., Bouffler, S. D., Lansdorp, P., and Bryant, P. E. (1997). Telomere Length, Chromatin Structure and Chromosome Fusigenic Potential. Chromosoma 106 (7), 413–421. doi:10.1007/s004120050263

PubMed Abstract | CrossRef Full Text | Google Scholar

Slijepcevic, P. (2016). Mechanisms of the Evolutionary Chromosome Plasticity: Integrating the 'Centromere-From-Telomere' Hypothesis with Telomere Length Regulation. Cytogenet. Genome Res. 148 (4), 268–278. doi:10.1159/000447415

PubMed Abstract | CrossRef Full Text | Google Scholar

Tree of Sex Consortium (2014). Tree of Sex: a Database of Sexual Systems. Sci. Data 1, 140015. doi:10.1038/sdata.2014.15

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, F., and Graphodatsky, A. S. (2009). “Animal Probes and ZOO-FISH,” in Fluorescence in Situ Hybridization (FISH)—application Guide (Berlin, Heidelberg: Springer), 323–346. doi:10.1007/978-3-540-70581-9_29

CrossRef Full Text | Google Scholar

Yang, F., Trifonov, V., Ng, B. L., Kosyakova, N., and Carter, N. P. (2009). “Generation of Paint Probes by Flow-Sorted and Microdissected Chromosomes,” in Fluorescence in Situ Hybridization (FISH) - Application Guide (Berlin, Heidelberg: Springer), 35–52. doi:10.1007/978-3-540-70581-9_3

CrossRef Full Text | Google Scholar

Yano, C. F., Poltronieri, J., Bertollo, L. A. C., Artoni, R. F., Liehr, T., and de Bello Cioffi, M. (2014). Chromosomal Mapping of Repetitive DNAs in Triportheus Trifurcatus (Characidae, Characiformes): Insights into the Differentiation of the Z and W Chromosomes. PLoS One 9 (3), e90946. doi:10.1371/journal.pone.0090946

PubMed Abstract | CrossRef Full Text | Google Scholar

Zwick, M. S., Hanson, R. E., Islam-Faridi, M. N., Stelly, D. M., Wing, R. A., Price, H. J., et al. (1997). A Rapid Procedure for the Isolation of C0t-1 DNA from Plants. Genome 40 (1), 138–142. doi:10.1139/g97-020

PubMed Abstract | CrossRef Full Text | Google Scholar