Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation

Conserved satellite DNA motif and lack of interstitial telomeric sites in highly rearranged African Nothobranchius killifish karyotypes

Using African annual killifishes of the genus Nothobranchius from temporary savannah pools with rapid karyotype and sex chromosome evolution, we analysed the chromosomal distribution of telomeric (TTAGGG)n repeat and Nfu-SatC satellite DNA (satDNA; isolated from Nothobranchius furzeri) in 15 species across the Nothobranchius killifish phylogeny, and with Fundulosoma thierryi as an out-group. Our fluorescence in situ hybridization experiments revealed that all analysed taxa share the presence of Nfu-SatC repeat but with diverse organization and distribution on chromosomes. Nfu-SatC landscape was similar in conspecific populations of Nothobranchius guentheri and Nothobranchius melanospilus but slightly-to-moderately differed between populations of Nothobranchius pienaari, and between closely related Nothobranchius kuhntae and Nothobranchius orthonotus. Inter-individual variability in Nfu-SatC patterns was found in N. orthonotus and Nothobranchius krysanovi. We revealed mostly no sex-linked patterns of studied repetitive DNA distribution. Only in Nothobranchius brieni, possessing multiple sex chromosomes, Nfu-SatC repeat occupied a substantial portion of the neo-Y chromosome, similarly as formerly found in the XY sex chromosome system of turquoise killifish N. furzeri and its sister species Nothobranchius kadleci-representatives not closely related to N. brieni. All studied species further shared patterns of expected telomeric repeats at the ends of all chromosomes and no additional interstitial telomeric sites. In summary, we revealed (i) the presence of conserved satDNA class in Nothobranchius clades (a rare pattern among ray-finned fishes); (ii) independent trajectories of Nothobranchius sex chromosome differentiation, with recurrent and convergent accumulation of Nfu-SatC on the Y chromosome in some species; and (iii) genus-wide shared tendency to loss of telomeric repeats during interchromosomal rearrangements. Collectively, our findings advance our understanding of genome structure, mechanisms of karyotype reshuffling, and sex chromosome differentiation in Nothobranchius killifishes from the genus-wide perspective.

Fast satellite DNA evolution in Nothobranchius annual killifishes

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius.

Revision of the Schistura cincticauda species group (Teleostei, Nemacheilidae) using molecular and morphological markers

To approach the taxonomy of large and complex animal groups it is of advantage to focus on species groups with shared derived character state. We investigate the composition, morphological characteristics and relationships of and within the Schistura cincticauda species group, whose members are small freshwater fishes that inhabit streams and rivers in eastern Myanmar and western and southern Thailand. A phylogenetic analysis using molecular genetic markers demonstrated the monophyly of this group; a combined genetic and morphological analysis revealed the inclusion of at least twelve species. They share the presence of a pair of black marks on the lower lip, one on each side of the median interruption (these marks may be reduced to few melanophores or even missing in some individuals). Additionally, all species share a small body size (max. 60 mm SL), an incomplete lateral line reaching at most to vertical through anal-fin base, and the absence of sexual dimorphism. Each of the 12 species is diagnosed by a unique combination of character states in fin ray numbers, anus position, presence/absence of an axillary pelvic lobe, and colour pattern. The distribution areas of several species overlap and five cases of syntopic occurrence are known. Five unnamed species are described herein.

Chromosomes of Asian Cyprinid Fishes: Genomic Differences in Conserved Karyotypes of 'Poropuntiinae' (Teleostei, Cyprinidae)

The representatives of cyprinid lineage 'Poropuntiinae' with 16 recognized genera and around 100 species form a significant part of Southeast Asian ichthyofauna. Cytogenetics are valuable when studying fish evolution, especially the dynamics of repetitive DNAs, such as ribosomal DNAs (5S and 18S) and microsatellites, that can vary between species. Here, karyotypes of seven 'poropuntiin' species, namely Cosmochilus harmandi, Cyclocheilichthys apogon, Hypsibarbus malcomi, H. wetmorei, Mystacoleucus chilopterus, M. ectypus, and Puntioplties proctozysron occurring in Thailand were examined using conventional and molecular cytogenetic protocols. Variable numbers of uni- and bi-armed chromosomes indicated widespread chromosome rearrangements with a stable diploid chromosome number (2n) of 50. Examination with fluorescence in situ hybridization using major and minor ribosomal probes showed that Cosmochilus harmandi, Cyclocheilichthys apogon, and Puntioplites proctozystron all had one chromosomal pair with 5S rDNA sites. However, more than two sites were found in Hypsibarbus malcolmi, H. wetmorei, Mystacoleucus chilopterus, and M. ectypus. The number of chromosomes with 18S rDNA sites varied amongst their karyotypes from one to three; additionally, comparative genomic hybridization and microsatellite patterns varied among species. Our results reinforce the trend of chromosomal evolution in cyprinifom fishes, with major chromosomal rearrangements, while conserving their 2n.

Robertsonian Fusion Site in Rineloricaria pentamaculata (Siluriformes: Loricariidae): Involvement of 5S Ribosomal DNA and Satellite Sequences

Cytogenetic studies demonstrated that unstable chromosomal sites in armored catfishes (Loricariidae) triggered intense karyotypic diversification, mainly derived from Robertsonian rearrangements. In Loricariinae, the presence of ribosomal DNA (rDNA) clusters and their flanking repeated regions (such as microsatellites or partial transposable element sequences) was proposed to facilitate chromosomal rearrangements. Hence, this study aimed to characterize the numerical chromosomal polymorphism observed in Rineloricaria pentamaculata and evaluate the chromosomal rearrangements which originated diploid chromosome number (2n) variation, from 56 to 54. Our data indicate a centric fusion event between acrocentric chromosomes of pairs 15 and 18, bearing 5S rDNA sites on their short (p) arms. This chromosome fusion established the numerical polymorphism, decreasing the 2n from original 56 (karyomorph A) to 55 in karyomorph B and 54 in karyomorph C. Although vestiges of telomeric sequences were evidenced at the fusion point, no 5S rDNA was detected in this region. The acrocentric chromosomes involved in the origin of the fusion were enriched with (CA)n and (GA)n microsatellites. Repetitive sequences in the acrocentric chromosomes subtelomeres have facilitated the rearrangement. Our study thus reinforces the view on the important role of particular repetitive DNA classes in promoting chromosome fusions which frequently drive Rineloricaria karyotype evolution.

Chromosomal rearrangements and the first indication of an ♀X1 X1 X2 X2 /♂X1 X2 Y sex chromosome system in Rineloricaria fishes (Teleostei: Siluriformes)

Rineloricaria is the most diverse genus within the freshwater fish subfamily Loricariinae, and it is widely distributed in the Neotropical region. Despite limited cytogenetic data, records from southern and south-eastern Brazil suggest a high rate of chromosomal rearrangements in this genus, mirrored in remarkable inter- and intraspecific karyotype variability. In the present work, we investigated the karyotype features of Rineloricaria teffeana, an endemic representative from northern Brazil, using both conventional and molecular cytogenetic techniques. We revealed different diploid chromosome numbers (2n) between sexes (33♂/34♀), which suggests the presence of an ♀X₁ X₁ X₂ X₂ /♂X₁ X₂ Y multiple sex chromosome system. The male-limited Y chromosome was the largest and the only biarmed element in the karyotype, implying Y-autosome fusion as the most probable mechanism behind its origination. C-banding revealed low amounts of constitutive heterochromatin, mostly confined to the (peri)centromeric regions of most chromosomes (including the X₂ and the Y) but also occupying the distal regions of a few chromosomal pairs. The chromosomal localization of the 18S ribosomal DNA (rDNA) clusters revealed a single site on chromosome pair 4, which was adjacent to the 5S rDNA cluster. Additional 5S rDNA loci were present on the autosome pair 8, X₁ chromosome, and in the presumed fusion point on the Y chromosome. The probe for telomeric repeat motif (TTAGGG)_n revealed signals of variable intensities at the ends of all chromosomes except for the Y chromosome, where no detectable signals were evidenced. Male-to-female comparative genomic hybridization revealed no sex-specific or sex-biased repetitive DNA accumulations, suggesting a presumably low level of neo-Y chromosome differentiation. We provide evidence that rDNA sites might have played a role in the formation of this putative multiple sex chromosome system and that chromosome fusions originate through different mechanisms among different Rineloricaria species.

Small Body, Large Chromosomes: Centric Fusions Shaped the Karyotype of the Amazonian Miniature Fish Nannostomus anduzei (Characiformes, Lebiasinidae)

Miniature refers to species with extraordinarily small adult body size when adult and can be found within all major metazoan groups. It is considered that miniature species have experienced severe alteration of numerous morphological traits during evolution. For a variety of reasons, including severe labor concerns during collecting, chromosomal acquisition, and taxonomic issues, miniature fishes are neglected and understudied. Since some available studies indicate possible relationship between diploid chromosome number (2n) and body size in fishes, we aimed to study one of the smallest Neotropical fish Nannostomus anduzei (Teleostei, Characiformes, Lebiasinidae), using both conventional (Giemsa staining, C-banding) and molecular cytogenetic methods (FISH mapping of rDNAs, microsatellites, and telomeric sequences). Our research revealed that N. anduzei possesses one of the lowest diploid chromosome numbers (2n = 22) among teleost fishes, and its karyotype is entirely composed of large metacentric chromosomes. All chromosomes, except for pair number 11, showed an 18S rDNA signal in the pericentromeric region. 5S rDNA signals were detected in the pericentromeric regions of chromosome pair number 1 and 6, displaying synteny to 18S rDNA signals. Interstitial telomeric sites (ITS) were identified in the centromeric region of pairs 6 and 8, indicating that centric fusions played a significant role in karyotype evolution of studied species. Our study provides further evidence supporting the trend of diploid chromosome number reduction along with miniaturization of adult body size in fishes.

Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci

Homomorphic sex chromosomes and their turnover are common in teleosts. We investigated the evolution of nascent sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci, focusing on their under-studied repetitive landscape. We combined bioinformatic analyses of the repeatome with molecular cytogenetic techniques, including comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, ribosomal RNA genes (rDNA) and bacterial artificial chromosomes (BACs), and immunostaining of SYCP3 and MLH1 proteins to mark lateral elements of synaptonemal complexes and recombination sites, respectively. Both species share the same heteromorphic XY sex chromosome system, which thus evolved prior to their divergence. This was corroborated by sequence analysis of a putative master sex determining (MSD) gene gdf6Y in both species. Based on their divergence, differentiation of the XY sex chromosome pair started approximately 2 million years ago. In all populations, the gdf6Y gene mapped within a region rich in satellite DNA on the Y chromosome long arms. Despite their heteromorphism, X and Y chromosomes mostly pair regularly in meiosis, implying synaptic adjustment. In N. kadleci, Y-linked paracentric inversions like those previously reported in N. furzeri were detected. An inversion involving the MSD gene may suppress occasional recombination in the region, which we otherwise evidenced in the N. furzeri population MZCS-121 of the Limpopo clade lacking this inversion. Y chromosome centromeric repeats were reduced compared with the X chromosome and autosomes, which points to a role of relaxed meiotic drive in shaping the Y chromosome repeat landscape. We speculate that the recombination rate between sex chromosomes was reduced due to heterochiasmy. The observed differences between the repeat accumulations on the X and Y chromosomes probably result from high repeat turnover and may not relate closely to the divergence inferred from earlier SNP analyses.

Karyotype variability in six Amazonian species of the family Curimatidae (Characiformes) revealed by repetitive sequence mapping

Fishes of the Curimatidae family represent one of the most important freshwater ichthyofauna groups of Central and South America, with 117 recognized species distributed in eight genera. In this study, six species - Curimata inornata, Curimatella dorsalis, and Psectrogaster falcata collected from the Lower Araguaia River, Pará, Brazil; Curimata vittata, Curimatella meyeri, and Psectrogaster rutiloides collected from the Catalão Lake, Amazonas, Brazil - were cytogenetically analyzed, investigate the occurrence and distribution of repetitive DNA classes in the karyotypes. All species had 2n=54 metacentric/submetacentric chromosomes. Despite the conservative diploid number, we observed variations in the karyotypic structure among species. Ribosomal DNA (rDNA) 18S and 5S were found in single or multiple sites, with the first report of synteny in Curimatella dorsalis, and the occurrence of several interstitial telomeric sequences (ITSs) in species of the genera Curimatella and Psectrogaster. Interspecific karyotypic diversity both concerning structure and location/position of the nucleolar organizer regions (NOR) and ribosomal DNA, suggesting the occurrence of several non-Robertsonian rearrangements driving the evolution of this family.

High chromosomal mobility of rDNA clusters in holocentric chromosomes of Triatominae, vectors of Chagas disease (Hemiptera-Reduviidae)

The subfamily Triatominae (Hemiptera-Reduviidae) includes more than 150 blood-sucking species, potential vectors of the protozoan Trypanosoma cruzi, causative agent of Chagas disease. A distinctive cytogenetic characteristic of this group is the presence of extremely stable chromosome numbers. Unexpectedly, the analyses of the chromosomal location of ribosomal gene clusters and other repetitive sequences place Triatominae as a significantly diverse hemipteran subfamily. Here, we advance the understanding of Triatominae chromosomal evolution through the analysis of the 45S rDNA cluster chromosomal location in 92 Triatominae species. We found the 45S rDNA clusters in one to four loci per haploid genome with different chromosomal patterns: On one or two autosomes, on one, two or three sex chromosomes, on the X chromosome plus one to three autosomes. The movement of 45S rDNA clusters is discussed in an evolutionary context. Our results illustrate that rDNA mobility has been relatively common in the past and in recent evolutionary history of the group. The high frequency of rDNA patterns involving autosomes and sex chromosomes among closely related species could affect genetic recombination and the viability of hybrid populations, which suggests that the mobility of rDNA clusters could be a driver of species diversification.

Using Species Groups to Approach the Large and Taxonomically Unresolved Freshwater Fish Family Nemacheilidae (Teleostei: Cypriniformes)

Large animal families with unresolved taxonomy are notoriously difficult to handle with respect to their biodiversity, systematics, and evolutionary history. We approach a large and taxonomically unresolved family of freshwater fishes (Nemacheilidae, >600 species) by proposing, on the basis of morphologic data, a species group within the family and study its phylogeny with conclusions regarding its diversity, taxonomy, and biogeographic history. Phylogenetic analyses of two mitochondrial and three nuclear genes of 139 specimens, representing about 46 species (17 candidate species from the proposed species-group, plus 29 comparative species), revealed that the proposed species group does not form a distinct monophyletic lineage, but that the candidate and comparative species mixed in three different lineages. However, the results revealed more than 20% of undescribed species within the ingroup and showed that species do not cluster according to the presently recognised genera. At least one of the genetic clades shows signs of an eastward range expansion during the second half of Miocene from north India via Myanmar into Laos, western China, and western Thailand. We conclude that the approach of picking monophyletic lineages to study biodiversity, systematics, and evolutionary history helps to open the door to large animal families.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

A supernumerary "B-sex" chromosome drives male sex determination in the Pachón cavefish, Astyanax mexicanus

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now.^1,2 B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species^3,4 that have often been considered as selfish genetic elements that behave as genome parasites.^5,6 The observation that in some species Bs can be either restricted or predominant in one sex^7-14 raised the interesting hypothesis that Bs could play a role in sex determination.¹⁵ The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene, which can promote male sex determination in this species.

Against the mainstream: exceptional evolutionary stability of ZW sex chromosomes across the fish families Triportheidae and Gasteropelecidae (Teleostei: Characiformes)

Teleost fishes exhibit a breath-taking diversity of sex determination and differentiation mechanisms. They encompass at least nine sex chromosome systems with often low degree of differentiation, high rate of inter- and intra-specific variability, and frequent turnovers. Nevertheless, several mainly female heterogametic systems at an advanced stage of genetic differentiation and high evolutionary stability have been also found across teleosts, especially among Neotropical characiforms. In this study, we aim to characterize the ZZ/ZW sex chromosome system in representatives of the Triportheidae family (Triportheus auritus, Agoniates halecinus, and the basal-most species Lignobrycon myersi) and its sister clade Gasteropelecidae (Carnegiella strigata, Gasteropelecus levis, and Thoracocharax stellatus). We applied both conventional and molecular cytogenetic approaches including chromosomal mapping of 5S and 18S ribosomal DNA clusters, cross-species chromosome painting (Zoo-FISH) with sex chromosome-derived probes and comparative genomic hybridization (CGH). We identified the ZW sex chromosome system for the first time in A. halecinus and G. levis and also in C. strigata formerly reported to lack sex chromosomes. We also brought evidence for possible mechanisms underlying the sex chromosome differentiation, including inversions, repetitive DNA accumulation, and exchange of genetic material. Our Zoo-FISH experiments further strongly indicated that the ZW sex chromosomes of Triportheidae and Gasteropelecidae are homeologous, suggesting their origin before the split of these lineages (approx. 40-70 million years ago). Such extent of sex chromosome stability is almost exceptional in teleosts, and hence, these lineages afford a special opportunity to scrutinize unique evolutionary forces and pressures shaping sex chromosome evolution in fishes and vertebrates in general.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Measurement of Chromosomal Arms and FISH Reveal Complex Genome Architecture and Standardized Karyotype of Model Fish, Genus Carassius

The widely distributed ray-finned fish genus Carassius is very well known due to its unique biological characteristics such as polyploidy, clonality, and/or interspecies hybridization. These biological characteristics have enabled Carassius species to be successfully widespread over relatively short period of evolutionary time. Therefore, this fish model deserves to be the center of attention in the research field. Some studies have already described the Carassius karyotype, but results are inconsistent in the number of morphological categories for individual chromosomes. We investigated three focal species: Carassius auratus, C. carassius and C. gibelio with the aim to describe their standardized diploid karyotypes, and to study their evolutionary relationships using cytogenetic tools. We measured length (q+plength) of each chromosome and calculated centromeric index (i value). We found: (i) The relationship between q+plength and i value showed higher similarity of C. auratus and C. carassius. (ii) The variability of i value within each chromosome expressed by means of the first quartile (Q1) up to the third quartile (Q3) showed higher similarity of C. carassius and C. gibelio. (iii) The fluorescent in situ hybridization (FISH) analysis revealed higher similarity of C. auratus and C. gibelio. (iv) Standardized karyotype formula described using median value (Q2) showed differentiation among all investigated species: C. auratus had 24 metacentric (m), 40 submetacentric (sm), 2 subtelocentric (st), 2 acrocentric (a) and 32 telocentric (T) chromosomes (24m+40sm+2st+2a+32T); C. carassius: 16m+34sm+8st+42T; and C. gibelio: 16m+22sm+10st+2a+50T. (v) We developed R scripts applicable for the description of standardized karyotype for any other species. The diverse results indicated unprecedented complex genomic and chromosomal architecture in the genus Carassius probably influenced by its unique biological characteristics which make the study of evolutionary relationships more difficult than it has been originally postulated.

One northward, one southward: Contrasting biogeographical history in two benthic freshwater fish genera across Southeast Asia (Teleostei: Cobitoidea: Nemacheilus, Pangio)

Southeast Asia is one of the world's biodiversity hotspots, and the high level of diversity and endemism was reached by colonisation events as well as internal diversification. We investigate the phylogenetic relationships and biogeographic history of the loach genus Nemacheilus, which is widely distributed and common across freshwaters of Southeast Asia. In addition we present the ancestral range reconstruction of the related loach genus Pangio that commonly occurs in the same region as Nemacheilus. Our results reveal that the species currently classified as Nemacheilus in fact are a polyphyletic assemblage; most species are now retaining in a monophyletic Nemacheilus sensu stricto and five species belong to different lineages. We further indicate the existence of hidden diversity within Nemacheilus in the form of several undescribed species. Three major clades (Selangoricus, Masyae and Ornatus) are found within the genus Nemacheilus sensu stricto. These clades generally correspond to the species groups formerly defined on the basis of their pigmentation pattern. The biogeographic analyses show that Nemacheilus most likely originated in mainland Southeast Asia and subsequently expanded in a southward direction to Borneo, Sumatra and Java and the southern Malay Peninsula. In contrast, the genus Pangio originated in Sundaland, from where it extended several times northwards into Indochina and to northern India. Our results demonstrate that small freshwater fishes with restricted dispersal ability are very helpful for the reconstruction of biogeographic history. The contrasting biogeographic history of these two groups of small, benthic and related fish show how complex and case-specific the processes that lead to the biodiversity richness of Southeast Asia are.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

Repetitive DNA Profiles Reveal Evidence of Rapid Genome Evolution and Reflect Species Boundaries in Ground Beetles

Genome architecture is a complex, multidimensional property of an organism defined by the content and spatial organization of the genome's component parts. Comparative study of entire genome architecture in model organisms is shedding light on mechanisms underlying genome regulation, evolution, and diversification, but such studies require costly analytical approaches which make extensive comparative study impractical for most groups. However, lower-cost methods that measure a single architectural component (e.g., distribution of one class of repeats) have potential as a new data source for evolutionary studies insofar as that measure correlates with more complex biological phenomena, and for which it could serve as part of an explanatory framework. We investigated copy number variation (CNV) profiles in ribosomal DNA (rDNA) as a simple measure reflecting the distribution of rDNA subcomponents across the genome. We find that signatures present in rDNA CNV profiles strongly correlate with species boundaries in the breve species group of Bembidion, and vary across broader taxonomic sampling in Bembidion subgenus Plataphus. Profiles of several species show evidence of re-patterning of rDNA-like sequences throughout the genome, revealing evidence of rapid genome evolution (including among sister pairs) not evident from analysis of traditional data sources such as multigene data sets. Major re-patterning of rDNA-like sequences has occurred frequently within the evolutionary history of Plataphus. We confirm that CNV profiles represent an aspect of genomic architecture (i.e., the linear distribution of rDNA components across the genome) via fluorescence in-situ hybridization. In at least one species, novel rDNA-like elements are spread throughout all chromosomes. We discuss the potential of copy number profiles of rDNA, or other repeats, as a low-cost tool for incorporating signal of genomic architecture variation in studies of species delimitation and genome evolution. [Bembidion; Carabidae; copy number variation profiles; rapid genome evolution; ribosomal DNA; species delimitation.].

Karyotype evolution and preliminary molecular assessment of genera in the family Scorpiopidae (Arachnida: Scorpiones)

The scorpions represent an ancient and morphologically conserved order of arachnids. Despite that, their karyotypes may differ considerably even among closely related species. In this study, we identify the trends of the karyotype evolution in the family Scorpiopidae based on integrating cytogenetic data and multi-locus molecular phylogenetic approaches. We detected considerable variability in diploid numbers of chromosomes (from 48 to 147), 18S rRNA gene cluster positions (from terminal to pericentromeric) at the interspecific level. Moreover, we identified independent fusions, fissions and inversions in the evolution of the family Scorpiopidae, leading to a remarkable diversification of the karyotypes. The dynamic system of the karyotype changes in this group is further documented by the presence of interstitial telomeric sequences (ITS) in two species. The cytogenetic differences observed among the analyzed species highlight the potential of this type of data for species-level taxonomy in scorpion lineages with monocentric chromosomes. Additionally, the results of our phylogenetic analyses support the monophyly of the family Scorpiopidae, but rendered several genera para- or polyphyletic.

Present and Future Salmonid Cytogenetics

Salmonids are extremely important economically and scientifically; therefore, dynamic developments in their research have occurred and will continue occurring in the future. At the same time, their complex phylogeny and taxonomy are challenging for traditional approaches in research. Here, we first provide discoveries regarding the hitherto completely unknown cytogenetic characteristics of the Anatolian endemic flathead trout, Salmo platycephalus, and summarize the presently known, albeit highly complicated, situation in the genus Salmo. Secondly, by outlining future directions of salmonid cytogenomics, we have produced a prototypical virtual karyotype of Salmo trutta, the closest relative of S. platycephalus. This production is now possible thanks to the high-quality genome assembled to the chromosome level in S. trutta via soft-masking, including a direct labelling of repetitive sequences along the chromosome sequence. Repetitive sequences were crucial for traditional fish cytogenetics and hence should also be utilized in fish cytogenomics. As such virtual karyotypes become increasingly available in the very near future, it is necessary to integrate both present and future approaches to maximize their respective benefits. Finally, we show how the presumably repetitive sequences in salmonids can change the understanding of the overall relationship between genome size and G+C content, creating another outstanding question in salmonid cytogenomics waiting to be resolved.

Highly Rearranged Karyotypes and Multiple Sex Chromosome Systems in Armored Catfishes from the Genus Harttia (Teleostei, Siluriformes)

Harttia comprises an armored catfish genus endemic to the Neotropical region, including 27 valid species with low dispersion rates that are restricted to small distribution areas. Cytogenetics data point to a wide chromosomal diversity in this genus due to changes that occurred in isolated populations, with chromosomal fusions and fissions explaining the 2n number variation. In addition, different multiple sex chromosome systems and rDNA loci location are also found in some species. However, several Harttia species and populations remain to be investigated. In this study, Harttia intermontana and two still undescribed species, morphologically identified as Harttia sp. 1 and Harttia sp. 2, were cytogenetically analyzed. Harttia intermontana has 2n = 52 and 2n = 53 chromosomes, while Harttia sp. 1 has 2n = 56 and 2n = 57 chromosomes in females and males, respectively, thus highlighting the occurrence of an XX/XY₁Y₂ multiple sex chromosome system in both species. Harttia sp. 2 presents 2n = 62 chromosomes for both females and males, with fission events explaining its karyotype diversification. Chromosomal locations of the rDNA sites were also quite different among species, reinforcing that extensive rearrangements had occurred in their karyotype evolution. Comparative genomic hybridization (CGH) experiments among some Harttia species evidenced a shared content of the XY₁Y₂ sex chromosomes in three of them, thus pointing towards their common origin. Therefore, the comparative analysis among all Harttia species cytogenetically studied thus far allowed us to provide an evolutionary scenario related to the speciation process of this fish group.

Chromosomes of Asian cyprinid fishes: Variable karyotype patterns and evolutionary trends in the genus Osteochilus (Cyprinidae, Labeoninae, "Osteochilini")

The Cyprinidae family is a highly diversified but demonstrably monophyletic lineage of cypriniform fishes. Among them, the genus Osteochilus contains 35 recognized valid species distributed from India, throughout Myanmar, Laos, Thailand, Malaysia, Indonesian archipelago to southern China. In this study, karyotypes and other chromosomal characteristics of five Osteochilus species occurring in Thailand, namely O. lini, O. melanopleura, O. microcephalus, O. vittatus and O. waandersii were examined using conventional and molecular cytogenetic protocols. Our results showed they possessed diploid chromosome number (2n) invariably 2n = 50, but the ratio of uni- and bi-armed chromosomes was highly variable among their karyotypes, indicating extensive chromosomal rearrangements. Only one chromosome pair bearing 5S rDNA sites occurred in most species, except O. melanopleura, where two sites were detected. In contrast, only one chromosomal pair bearing 18S rDNA sites were observed among their karyotypes, but in different positions. These cytogenetic patterns indicated that the cytogenomic divergence patterns of these Osteochilus species were largely corresponding to the inferred phylogenetic tree. Similarly, different patterns of the distributions of rDNAs and microsatellites across genomes of examined species as well as their different karyotype structures indicated significant evolutionary differentiation of Osteochilus genomes.

Patterns of Sex Chromosome Differentiation in Spiders: Insights from Comparative Genomic Hybridisation

Spiders are an intriguing model to analyse sex chromosome evolution because of their peculiar multiple X chromosome systems. Y chromosomes were considered rare in this group, arising after neo-sex chromosome formation by X chromosome-autosome rearrangements. However, recent findings suggest that Y chromosomes are more common in spiders than previously thought. Besides neo-sex chromosomes, they are also involved in the ancient X1X2Y system of haplogyne spiders, whose origin is unknown. Furthermore, spiders seem to exhibit obligatorily one or two pairs of cryptic homomorphic XY chromosomes (further cryptic sex chromosome pairs, CSCPs), which could represent the ancestral spider sex chromosomes. Here, we analyse the molecular differentiation of particular types of spider Y chromosomes in a representative set of ten species by comparative genomic hybridisation (CGH). We found a high Y chromosome differentiation in haplogyne species with X1X2Y system except for Loxosceles spp. CSCP chromosomes exhibited generally low differentiation. Possible mechanisms and factors behind the observed patterns are discussed. The presence of autosomal regions marked predominantly or exclusively with the male or female probe was also recorded. We attribute this pattern to intraspecific variability in the copy number and distribution of certain repetitive DNAs in spider genomes, pointing thus to the limits of CGH in this arachnid group. In addition, we confirmed nonrandom association of chromosomes belonging to particular CSCPs at spermatogonial mitosis and spermatocyte meiosis and their association with multiple Xs throughout meiosis. Taken together, our data suggest diverse evolutionary pathways of molecular differentiation in different types of spider Y chromosomes.

Resolving an unnoticed diversity within the Schistura robertsi species complex (Teleostei: Nemacheilidae) using molecules and morphology

The Schistura robertsi species complex is a group of freshwater fish inhabiting streams in southeast Myanmar as well as in western and southern Thailand. In southern Thailand, the distribution exceeds the biogeographically important 'Surat Thani - Krabi line'. The complex is believed to include five described and one undescribed species, but monophyly and systematics of the group have never been studied explicitly. The present study aims to resolve the number of species within the Schistura robertsi group as well as their distribution areas and phylogenetic relations. We analysed mitochondrial and nuclear sequence data of 86 specimens from 47 localities and 18 morphological characters of 193 specimens. The phylogenetic analyses revealed the S. robertsi complex to be monophyletic and to be composed of ten major lineages. Six of them correspond to the known described or undescribed species, but another four newly identified clades reveal the existence of an overlooked diversity within the group. All genetic lineages are statistically highly supported and all are morphologically diagnosable, suggesting that they represent distinct species. The distribution areas of several clades overlap, the cases of direct co-occurrence show no sign of hybridisation.

Taxonomic Diversity Not Associated with Gross Karyotype Differentiation: The Case of Bighead Carps, Genus Hypophthalmichthys (Teleostei, Cypriniformes, Xenocyprididae)

The bighead carps of the genus Hypophthalmichthys (H. molitrix and H. nobilis) are important aquaculture species. They were subjected to extensive multidisciplinary research, but with cytogenetics confined to conventional protocols only. Here, we employed Giemsa-/C-/CMA3- stainings and chromosomal mapping of multigene families and telomeric repeats. Both species shared (i) a diploid chromosome number 2n = 48 and the karyotype structure, (ii) low amount of constitutive heterochromatin, (iii) the absence of interstitial telomeric sites (ITSs), (iv) a single pair of 5S rDNA loci adjacent to one major rDNA cluster, and (v) a single pair of co-localized U1/U2 snDNA tandem repeats. Both species, on the other hand, differed in (i) the presence/absence of remarkable interstitial block of constitutive heterochromatin on the largest acrocentric pair 11 and (ii) the number of major (CMA3-positive) rDNA sites. Additionally, we applied here, for the first time, the conventional cytogenetics in H. harmandi, a species considered extinct in the wild and/or extensively cross-hybridized with H. molitrix. Its 2n and karyotype description match those found in the previous two species, while silver staining showed differences in distribution of major rDNA. The bighead carps thus represent another case of taxonomic diversity not associated with gross karyotype differentiation, where 2n and karyotype structure cannot help in distinguishing between genomes of closely related species. On the other hand, we demonstrated that two cytogenetic characters (distribution of constitutive heterochromatin and major rDNA) may be useful for diagnosis of pure species. The universality of these markers must be further verified by analyzing other pure populations of bighead carps.

Chromosome-Level Genome Reveals the Origin of Neo-Y Chromosome in the Male Barred Knifejaw Oplegnathus fasciatus

The barred knifejaw, Oplegnathus fasciatus, is characterized by an X₁X₂Y system with a neo-Y chromosome for males. Here, a chromosome-level genome was assembled to investigate the origin of neo-Y chromosome to the male O. fasciatus. Twenty-three chromosomes corresponding to the male karyotypes were scaffolded to 762-Mb genome with a contig N50 length of 2.18 Mb. A large neo-Y chromosome (Ch9) in the male O. fasciatus genome was also assembled and exhibited high identity to those of the female chromosomes Ch8 and Ch10. Chromosome rearrangements events were detected in the neo-chromosome Ch9. Our results suggested that a centric fusion of acrocentric chromosomes Ch8 and Ch10 should be responsible for the formation of the X₁X₂Y system. The high-quality genome will not only provide a solid foundation for further sex-determining mechanism research in the X₁X₂Y system but also facilitate the artificial breeding aiming to improve the yield and disease resistance for Oplegnathus.

•

Construction of a chromosome-level reference genome for the male O. fasciatus

•

Identification of the origin of neo-Y chromosome to the X₁X₂Y system

•

Accurate comparisons of sequences and genes between female X₁X₁X₂X₂ and male X₁X₂Y

An Insight into the Chromosomal Evolution of Lebiasinidae (Teleostei, Characiformes)

Lebiasinidae fishes have been historically neglected by cytogenetical studies. Here we present a genomic comparison in eleven Lebiasinidae species, in addition to a review of the ribosomal DNA sequences distribution in this family. With that, we develop ten sets of experiments in order to hybridize the genomic DNA of representative species from the genus Copeina, Copella, Nannostomus, and Pyrrhulina in metaphase plates of Lebiasina melanoguttata. Two major pathways on the chromosomal evolution of these species can be recognized: (i) conservation of 2n = 36 bi-armed chromosomes in Lebiasininae, as a basal condition, and (ii) high numeric and structural chromosomal rearrangements in Pyrrhulininae, with a notable tendency towards acrocentrization. The ribosomal DNA (rDNA) distribution also revealed a marked differentiation during the chromosomal evolution of Lebiasinidae, since both single and multiple sites, in addition to a wide range of chromosomal locations can be found. With some few exceptions, the terminal position of 18S rDNA appears as a common feature in Lebiasinidae-analyzed species. Altogether with Ctenoluciidae, this pattern can be considered a symplesiomorphism for both families. In addition to the specific repetitive DNA content that characterizes the genome of each particular species, Lebiasina also keeps inter-specific repetitive sequences, thus reinforcing its proposed basal condition in Lebiasinidae.

First Record on Sex Chromosomes in a Species of the Family Cynodontidae: Cynodon gibbus (Agassiz, 1829)

The fish family Cynodontidae belongs to the superfamily Curimatoidea, together with the Hemiodontidae, Serrasalmidae, Parodontidae, Prochilodontidae, Chilodontidae, Curimatidae, and Anostomidae. The majority of the species of this superfamily that have been analyzed to date have a diploid chromosome number of 2n = 54. Differentiated sex chromosomes (with female heterogamety) have been observed only in the Prochilodontidae, Parodontidae, and Anostomidae. The present study provides the first description of differentiated sex chromosomes in the cynodontid species Cynodon gibbus, which has a ZZ/ZW system, and shows that repetitive DNA has played a fundamental role in the differentiation of these sex chromosomes.

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

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

Centric Fusions behind the Karyotype Evolution of Neotropical Nannostomus Pencilfishes (Characiforme, Lebiasinidae): First Insights from a Molecular Cytogenetic Perspective

Lebiasinidae is a Neotropical freshwater family widely distributed throughout South and Central America. Due to their often very small body size, Lebiasinidae species are cytogenetically challenging and hence largely underexplored. However, the available but limited karyotype data already suggested a high interspecific variability in the diploid chromosome number (2n), which is pronounced in the speciose genus Nannostomus, a popular taxon in ornamental fish trade due to its remarkable body coloration. Aiming to more deeply examine the karyotype diversification in Nannostomus, we combined conventional cytogenetics (Giemsa-staining and C-banding) with the chromosomal mapping of tandemly repeated 5S and 18S rDNA clusters and with interspecific comparative genomic hybridization (CGH) to investigate genomes of four representative Nannostomus species: N. beckfordi, N. eques, N. marginatus, and N. unifasciatus. Our data showed a remarkable variability in 2n, ranging from 2n = 22 in N. unifasciatus (karyotype composed exclusively of metacentrics/submetacentrics) to 2n = 44 in N. beckfordi (karyotype composed entirely of acrocentrics). On the other hand, patterns of 18S and 5S rDNA distribution in the analyzed karyotypes remained rather conservative, with only two 18S and two to four 5S rDNA sites. In view of the mostly unchanged number of chromosome arms (FN = 44) in all but one species (N. eques; FN = 36), and with respect to the current phylogenetic hypothesis, we propose Robertsonian translocations to be a significant contributor to the karyotype differentiation in (at least herein studied) Nannostomus species. Interspecific comparative genome hybridization (CGH) using whole genomic DNAs mapped against the chromosome background of N. beckfordi found a moderate divergence in the repetitive DNA content among the species’ genomes. Collectively, our data suggest that the karyotype differentiation in Nannostomus has been largely driven by major structural rearrangements, accompanied by only low to moderate dynamics of repetitive DNA at the sub-chromosomal level. Possible mechanisms and factors behind the elevated tolerance to such a rate of karyotype change in Nannostomus are discussed.

Deciphering the Evolutionary History of Arowana Fishes (Teleostei, Osteoglossiformes, Osteoglossidae): Insight from Comparative Cytogenomics

Arowanas (Osteoglossinae) are charismatic freshwater fishes with six species and two genera (Osteoglossum and Scleropages) distributed in South America, Asia, and Australia. In an attempt to provide a better assessment of the processes shaping their evolution, we employed a set of cytogenetic and genomic approaches, including i) molecular cytogenetic analyses using C- and CMA₃/DAPI staining, repetitive DNA mapping, comparative genomic hybridization (CGH), and Zoo-FISH, along with ii) the genotypic analyses of single nucleotide polymorphisms (SNPs) generated by diversity array technology sequencing (DArTseq). We observed diploid chromosome numbers of 2n = 56 and 54 in O. bicirrhosum and O. ferreirai, respectively, and 2n = 50 in S. formosus, while S. jardinii and S. leichardti presented 2n = 48 and 44, respectively. A time-calibrated phylogenetic tree revealed that Osteoglossum and Scleropages divergence occurred approximately 50 million years ago (MYA), at the time of the final separation of Australia and South America (with Antarctica). Asian S. formosus and Australian Scleropages diverged about 35.5 MYA, substantially after the latest terrestrial connection between Australia and Southeast Asia through the Indian plate movement. Our combined data provided a comprehensive perspective of the cytogenomic diversity and evolution of arowana species on a timescale.

Extensive Chromosomal Reorganization in Apistogramma Fishes (Cichlidae, Cichlinae) Fits the Complex Evolutionary Diversification of the Genus

Neotropical cichlid fishes are one of the most diversified and evolutionarily successful species assemblages. Extremely similar forms and intraspecific polychromatism present challenges for the taxonomy of some of these groups. Several species complexes have a largely unknown origin and unresolved evolutionary processes. Dwarf cichlids of the genus Apistogramma, comprising more than a hundred species, exhibit intricate taxonomic and biogeographic patterns, with both allopatric and sympatric distributions. However, karyotype evolution and the role of chromosomal changes in Apistogramma are still unknown. In the present study, nine South American Apistogramma species were analyzed using conventional cytogenetic methods and the mapping of repetitive DNA sequences [18S rDNA, 5S rDNA, and (TTAGGG)n] by fluorescence in situ hybridization (FISH). Our results showed that Apistogramma has unique cytogenetic characteristics in relation to closely related groups, such as a reduced 2n and a large number of bi-armed chromosomes. Interspecific patterns revealed a scenario of remarkable karyotypic changes, including a reduction of 2n, the occurrence of B-chromosomes and evolutionary dynamic of rDNA tandem repeats. In addition to the well-known pre-zygotic reproductive isolation, the karyotype reorganization in the genus suggests that chromosomal changes could act as postzygotic barriers in areas where Apistogramma congeners overlap.

Comparative Cytogenetics and Neo-Y Formation in Small-Sized Fish Species of the Genus Pyrrhulina (Characiformes, Lebiasinidae)

Although fishes have traditionally been the subject of comparative evolutionary studies, few reports have concentrated on the application of multipronged modern molecular cytogenetic techniques (such as comparative genomic hybridization = CGH and whole chromosome painting = WCP) to analyze deeper the karyotype evolution of specific groups, especially the historically neglected small-sized ones. Representatives of the family Lebiasinidae (Characiformes) are a notable example, where only a few cytogenetic investigations have been conducted thus far. Here, we aim to elucidate the evolutionary processes behind the karyotype differentiation of Pyrrhulina species on a finer-scale cytogenetic level. To achieve this, we applied C-banding, repetitive DNA mapping, CGH and WCP in Pyrrhulina semifasciata and P. brevis. Our results showed 2n = 42 in both sexes of P. brevis, while the difference in 2n between male and female in P. semifasciata (♂41/♀42) stands out due to the presence of a multiple X₁X₂Y sex chromosome system, until now undetected in this family. As a remarkable common feature, multiple 18S and 5S rDNA sites are present, with an occasional synteny or tandem-repeat amplification. Male-vs.-female CGH experiments in P. semifasciata highlighted the accumulation of male-enriched repetitive sequences in the pericentromeric region of the Y chromosome. Inter-specific CGH experiments evidenced a divergence between both species’ genomes based on the presence of several species-specific signals, highlighting their inner genomic diversity. WCP with the P. semifasciata-derived Y (PSEMI-Y) probe painted not only the entire metacentric Y chromosome in males but also the X₁ and X₂ chromosomes in both male and female chromosomes of P. semifasciata. In the cross-species experiments, the PSEMI-Y probe painted four acrocentric chromosomes in both males and females of the other tested Pyrrhulina species. In summary, our results show that both intra- and interchromosomal rearrangements together with the dynamics of repetitive DNA significantly contributed to the karyotype divergence among Pyrrhulina species, possibly promoted by specific populational and ecological traits and accompanied in one species by the origin of neo-sex chromosomes. The present results suggest how particular evolutionary scenarios found in fish species can help to clarify several issues related to genome organization and the karyotype evolution of vertebrates in general.

Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes)

Oplegnathus fasciatus and O. punctatus (Teleostei: Centrarchiformes: Oplegnathidae), are commercially important rocky reef fishes, endemic to East Asia. Both species present an X1X2Y sex chromosome system. Here, we investigated the evolutionary forces behind the origin and differentiation of these sex chromosomes, with the aim to elucidate whether they had a single or convergent origin. To achieve this, conventional and molecular cytogenetic protocols, involving the mapping of repetitive DNA markers, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) were applied. Both species presented similar 2n, karyotype structure and hybridization patterns of repetitive DNA classes. 5S rDNA loci, besides being placed on the autosomal pair 22, resided in the terminal region of the long arms of both X1 chromosomes in females, and on the X1 and Y chromosomes in males. Furthermore, WCP experiments with a probe derived from the Y chromosome of O. fasciatus (OFAS-Y) entirely painted the X1 and X2 chromosomes in females and the X1, X2, and Y chromosomes in males of both species. CGH failed to reveal any sign of sequence differentiation on the Y chromosome in both species, thereby suggesting the shared early stage of neo-Y chromosome differentiation. Altogether, the present findings confirmed the origin of the X1X2Y sex chromosomes via Y-autosome centric fusion and strongly suggested their common origin.

Genomic Organization of Repetitive DNA Elements and Extensive Karyotype Diversity of Silurid Catfishes (Teleostei: Siluriformes): A Comparative Cytogenetic Approach

The catfish family Siluridae contains 107 described species distributed in Asia, but with some distributed in Europe. In this study, karyotypes and other chromosomal characteristics of 15 species from eight genera were examined using conventional and molecular cytogenetic protocols. Our results showed the diploid number (2n) to be highly divergent among species, ranging from 2n = 40 to 92, with the modal frequency comprising 56 to 64 chromosomes. Accordingly, the ratio of uni- and bi-armed chromosomes is also highly variable, thus suggesting extensive chromosomal rearrangements. Only one chromosome pair bearing major rDNA sites occurs in most species, except for Wallago micropogon, Ompok siluroides, and Kryptoterus giminus with two; and Silurichthys phaiosoma with five such pairs. In contrast, chromosomes bearing 5S rDNA sites range from one to as high as nine pairs among the species. Comparative genomic hybridization (CGH) experiments evidenced large genomic divergence, even between congeneric species. As a whole, we conclude that karyotype features and chromosomal diversity of the silurid catfishes are unusually extensive, but parallel some other catfish lineages and primary freshwater fish groups, thus making silurids an important model for investigating the evolutionary dynamics of fish chromosomes.

Cryptic diversity and dynamic chromosome evolution in Alpine scorpions (Euscorpiidae: Euscorpius)

Over time, mountain biota has undergone complex evolutionary histories that have left imprints on its genomic arrangement, geographical distribution and diversity of contemporary lineages. Knowledge on these biogeographical aspects still lags behind for invertebrates inhabiting the Alpine region. In the present study, we examined three scorpion species of the subgenus Euscorpius (Alpiscorpius) from the European Alps using cytogenetic and molecular phylogenetic approaches to determine the variation and population structure of extant lineages at both chromosome and genetic level, and to provide an insight into the species diversification histories. We detected considerable intraspecific variability in chromosome complements and localization of the 18S rDNA loci in all studied species. Such chromosome differences were noticeable as the existence of three [in E. (A.) alpha and E. (A.) germanus] or four [in E. (A.) gamma] range-restricted karyotypic races. These races differed from one another either by 2n [in E. (A.) alpha 2n = 54, 60, 90; in E. (A.) gamma 2n = 58, 60, 88, 86-92], or by the karyotypic formula [in E. (A.) germanus 2n = 34m + 12sm; 36m + 10sm; 42m + 4sm]. Using mitochondrial (16S rRNA, COI) and nuclear (28S rDNA) genetic markers, we examined genetic variation and reconstructed phylogenetic relationships among the karyotypic races. Both approaches provided evidence for the existence of ten deeply divergent lineages exhibiting the features of local endemics and indicating the presence of cryptic species. Molecular dating analyses suggest that these lineages diversified during the Plio-Pleistocene and this process was presumably accompanied by dynamic structural changes in the genome organization.

Chromosomes of Asian cyprinid fishes: cytogenetic analysis of two representatives of small paleotetraploid tribe Probarbini

Background

Polyploidy, although still poorly explored, represents an important evolutionary event in several cyprinid clades. Herein, Catlocarpio siamensis and Probarbus jullieni - representatives of the paleotetraploid tribe Probarbini, were characterized both by conventional and molecular cytogenetic methods.

Results

Alike most other paleotetraploid cyprinids (with 2n = 100), both species studied here shared 2n = 98 but differed in karyotypes: C. siamensis displayed 18m + 34sm + 46st/a; NF = 150, while P. jullieni exhibited 26m + 14sm + 58st/a; NF = 138. Fluorescence in situ hybridization (FISH) with rDNA probes revealed two (5S) and eight (18S) signals in C. siamensis, respectively, and six signals for both probes in P. jullieni. FISH with microsatellite motifs evidenced substantial genomic divergence between both species. The almost doubled size of the chromosome pairs #1 in C. siamensis and #14 in P. jullieni compared to the rest of corresponding karyotypes indicated chromosomal fusions.

Conclusion

Based on our findings, together with likely the same reduced 2n = 98 karyotypes in the remainder Probarbini species, we hypothesize that the karyotype 2n = 98 might represent a derived character, shared by all members of the Probarbini clade. Besides, we also witnessed considerable changes in the amount and distribution of certain repetitive DNA classes, suggesting complex post-polyploidization processes in this small paleotetraploid tribe.

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

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

Evolutionary trends in animal ribosomal DNA loci: introduction to a new online database

Ribosomal DNA (rDNA) loci encoding 5S and 45S (18S-5.8S-28S) rRNAs are important components of eukaryotic chromosomes. Here, we set up the animal rDNA database containing cytogenetic information about these loci in 1343 animal species (264 families) collected from 542 publications. The data are based on in situ hybridisation studies (both radioactive and fluorescent) carried out in major groups of vertebrates (fish, reptiles, amphibians, birds, and mammals) and invertebrates (mostly insects and mollusks). The database is accessible online at www.animalrdnadatabase.com . The median number of 45S and 5S sites was close to two per diploid chromosome set for both rDNAs despite large variation (1-74 for 5S and 1-54 for 45S sites). No significant correlation between the number of 5S and 45S rDNA loci was observed, suggesting that their distribution and amplification across the chromosomes follow independent evolutionary trajectories. Each group, irrespective of taxonomic classification, contained rDNA sites at any chromosome location. However, the distal and pericentromeric positions were the most prevalent (> 75% karyotypes) for 45S loci, while the position of 5S loci was more variable. We also examined potential relationships between molecular attributes of rDNA (homogenisation and expression) and cytogenetic parameters such as rDNA positions, chromosome number, and morphology.

Molecular technique reveals high variability of 18S rDNA distribution in harvestmen (Opiliones, Phalangiidae) from South Africa

The knowledge of cytogenetics in the harvestmen family Phalangiidae has been based on taxa from the Northern Hemisphere. We performed cytogenetic analysis on Guruia africana (Karsch, 1878) (2n=24) and four species of the genus Rhampsinitus Simon, 1879 (2n=24, 26, 34) from South Africa. Fluorescence in situ hybridization with an 18S rDNA probe was used to analyze the number and the distribution of this cluster in the family Phalangiidae for the first time. The results support the cytogenetic characteristics typical for the majority of harvestmen taxa, i.e. the predominance of small biarmed chromosomes and the absence of morphologically well-differentiated sex chromosomes as an ancestral state. We identified the number of 18S rDNA sites ranging from two in R. qachasneki Kauri, 1962 to seven in one population of R. leighi Pocock, 1903. Moreover, we found differences in the number and localization of 18S rDNA sites in R. leighi between populations from two localities and between sexes of R. capensis (Loman, 1898). The heterozygous states of the 18S rDNA sites in these species may indicate the presence of XX/XY and ZZ/ZW sex chromosomes, and the possible existence of these systems in harvestmen is discussed. The variability of the 18S rDNA sites indicates intensive chromosomal changes during the differentiation of the karyotypes, which is in contrast to the usual uniformity in chromosomal morphology known from harvestmen so far.

Evolutionary insight on localization of 18S, 28S rDNA genes on homologous chromosomes in Primates genomes

We explored the topology of 18S and 28S rDNA units by fluorescence in situ hybridization (FISH) in the karyotypes of thirteen species representatives from major groups of Primates and Tupaia minor (Günther, 1876) (Scandentia), in order to expand our knowledge of Primate genome reshuffling and to identify the possible dispersion mechanisms of rDNA sequences. We documented that rDNA probe signals were identified on one to six pairs of chromosomes, both acrocentric and metacentric ones. In addition, we examined the potential homology of chromosomes bearing rDNA genes across different species and in a wide phylogenetic perspective, based on the DAPI-inverted pattern and their synteny to human. Our analysis revealed an extensive variability in the topology of the rDNA signals across studied species. In some cases, closely related species show signals on homologous chromosomes, thus representing synapomorphies, while in other cases, signal was detected on distinct chromosomes, leading to species specific patterns. These results led us to support the hypothesis that different mechanisms are responsible for the distribution of the ribosomal DNA cluster in Primates.

Morphologically indistinguishable hybrid Carassius female with 156 chromosomes: A threat for the threatened crucian carp, C. carassius, L

The crucian carp Carassius carassius (Linnaeus, 1758), is native to many European freshwaters. Despite its wide distribution, the crucian carp is declining in both the number and sizes of populations across much of its range. Here we studied 30 individuals of a putative pure population from Helsinki, Finland. Despite clear external morphological features of C. carassius, an individual was of a higher ploidy level than the others. We therefore applied a set of molecular genetic (S7 nuclear and cytochrome b mitochondrial genes) and cytogenetic tools (sequential fluorescent 4’, 6-diamidino-2-phenylindole [DAPI], Chromomycin A₃ [CMA₃], C-banding and in situ hybridization [FISH] with both 5S and 28S ribosomal DNA probes) to determine its origin. While all examined characteristics of a diploid representative male (CCAHe2Fi) clearly corresponded to those of C. carassius, a triploid individual (CCAHe1Fi) was more complex. Phylogenetic analysis revealed that the nuclear genome of CCAHe1Fi contained three haploid sets: two C. gibelio and one C. carassius. However the mitochondrial DNA was that of C. gibelio, demonstrating its hybrid origin. The FISH revealed three strong (more intensive) 5S rDNA loci, confirming the triploid status, and an additional 24 weak (less intensive) signals were observed in the chromosome complement of CCAHe1Fi. On the other hand, only two strong and 16 weak 5S rDNA signals were visible on the chromosomes of the CCAHe2Fi male. 28S rDNA FISH revealed four strong signals in both CCAHe1Fi and CCAHe2Fi individuals. CMA₃ staining revealed four to six CMA₃-positive bands of CCAHe1Fi, while that of diploids contained only two to four. The fact that a polyploid hybrid Carassius female with a strong invasive potential may share morphological characters typical for endangered C. carassius highlights a need to combine genetic investigations of Carassius cryptic diversity with conservation measures of C. carassius in Europe.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Molecular cytogenetic differentiation of paralogs of Hox paralogs in duplicated and re-diploidized genome of the North American paddlefish (Polyodon spathula)

Background

Acipenseriformes is a basal lineage of ray-finned fishes and comprise 27 extant species of sturgeons and paddlefishes. They are characterized by several specific genomic features as broad ploidy variation, high chromosome numbers, presence of numerous microchromosomes and propensity to interspecific hybridization. The presumed palaeotetraploidy of the American paddlefish was recently validated by molecular phylogeny and Hox genes analyses. A whole genome duplication in the paddlefish lineage was estimated at approximately 42 Mya and was found to be independent from several genome duplications evidenced in its sister lineage, i.e. sturgeons. We tested the ploidy status of available chromosomal markers after the expected rediploidization. Further we tested, whether paralogs of Hox gene clusters originated from this paddlefish specific genome duplication are cytogenetically distinguishable.

Results

We found that both paralogs HoxA alpha and beta were distinguishable without any overlapping of the hybridization signal - each on one pair of large metacentric chromosomes. Of the HoxD, only the beta paralog was unequivocally identified, whereas the alpha paralog did not work and yielded only an inconclusive diffuse signal. Chromosomal markers on three diverse ploidy levels reflecting different stages of rediploidization were identified: quadruplets retaining their ancestral tetraploid condition, semi-quadruplets still reflecting the ancestral tetraploidy with clear signs of advanced rediploidization, doublets were diploidized with ancestral tetraploidy already blurred. Also some of the available microsatellite data exhibited diploid allelic band patterns at their loci whereas another locus showed more than two alleles.

Conclusions

Our exhaustive staining of paddlefish chromosomes combined with cytogenetic mapping of ribosomal genes and Hox paralogs and with microsatellite data, brings a closer look at results of the process of rediploidization in the course of paddlefish genome evolution. We show a partial rediploidization represented by a complex mosaic structure comparable with segmental paleotetraploidy revealed in sturgeons (Acipenseridae). Sturgeons and paddlefishes with their high propensity for whole genome duplication thus offer suitable animal model systems to further explore evolutionary processes that were shaping the early evolution of all vertebrates.

Electronic supplementary material

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

Karyotype description of the African weakly electric fish Campylomormyrus compressirostris in the context of chromosome evolution in Osteoglossiformes

Karyotyping is a basic method to investigate chromosomal evolution and genomic rearrangements. Sixteen genera within the basal teleost order Osteoglossiformes are currently described cytogenetically. Our study adds information to this chromosomal dataset by determining the karyotype of Campylomormyrus compressirostris, a genus of African weakly electric fish that has not been previously examined. Our results indicate a diploid chromosome number of 2n=48 (4sm+26m+18a) with a fundamental number of FN=72. This chromosome number is identical to the number documented for the sister taxon of the genus Campylomormyrus, i.e., Gnathonemus petersii (2n=48). These results support the close relationship of Campylomormyrus and Gnathonemus. However, the karyotype formula of C. compressirostris is different from Gnathonemus petersii, thereby confirming the high variability of karyotype formulae within the Mormyridae. We infer that the differences in chromosome number and formula of Campylomormyrus relative to other mormyrids may be caused by Robertsonian fusion and pericentric inversion. In addition to the karyotype description and classification of Campylomormyrus, a ChromEvol analysis was used to determine the ancestral haploid chromosome number of osteoglossiform taxa. Our results indicate a relatively conservative haploid chromosome number of n=24 for the most recent common ancestor of Osteoglossiformes and for most of the internal nodes of osteoglossiform phylogeny. Hence, we presume that the high chromosome variability evolved recently on multiple independent occasions. Furthermore, we suggest that the most likely ancestral chromosome number of Mormyridae is either n=24 or n=25. To the best of our knowledge this is the first attempt to determine and classify the karyotype of the weakly electric fish genus Campylomormyrus and to analyze chromosomal evolution within the Osteoglossiformes based on Maximum Likelihood and Bayesian Inference analyses.

Hypervariability of Nucleolus Organizer Regions in Bengal Slow Lorises, Nycticebus bengalensis (Primates, Lorisidae)

Slow lorises are a cryptic species complex, and thus genetic markers are needed to identify distinct evolutionary lineages or species. We examined the nucleolus organizer regions (NORs) of Bengal slow lorises (Nycticebus bengalensis) using FISH with 18S rDNA (rDNA-FISH) and silver nitrate staining (Ag-NOR stain). Ten individuals of the putatively single species N. bengalensis showed higher variability in localization than 3 other congeners, though their overall karyotypes were similar. The rDNA-FISH analysis detected a total of 18 loci, in contrast to previous studies of other slow loris species that revealed far fewer (6-10) loci. Eight of the 18 loci detected in the present analysis were found to be semi-stable localizations at 4 different chromosomes, while 10 were found to be unstable localizations at 5 other chromosomes. The semi-stable locations showed occasional presence/absence of variations for rDNA-FISH, and unstable locations were polymorphic among individuals, contributing to the higher variability of NORs in this taxon. We hypothesize that the larger numbers of rDNA loci found in N. bengalensis were introduced by genomic dispersion through ectopic recombination in association with terminal regions including rDNA. Such differences are potentially very powerful chromosomal markers to be used in species identification and conservation.