"Holostei versus Halecostomi" Problem: Insight from Cytogenetics of Ancient Nonteleost Actinopterygian Fish, Bowfin Amia calva

Abandoning the Isochore Theory Can Help Explain Genome Compositional Organization in Fish

The organization of the genome nucleotide (AT/GC) composition in vertebrates remains poorly understood despite the numerous genome assemblies available. Particularly, the origin of the AT/GC heterogeneity in amniotes, in comparison to the homogeneity in anamniotes, is controversial. Recently, several exceptions to this dichotomy were confirmed in an ancient fish lineage with mammalian AT/GC heterogeneity. Hence, our current knowledge necessitates a reevaluation considering this fact and utilizing newly available data and tools. We analyzed fish genomes in silico with as low user input as possible to compare previous approaches to assessing genome composition. Our results revealed a disparity between previously used plots of GC% and histograms representing the authentic distribution of GC% values in genomes. Previous plots heavily reduced the range of GC% values in fish to comply with the alleged AT/GC homogeneity and AT-richness of their genomes. We illustrate how the selected sequence size influences the clustering of GC% values. Previous approaches that disregarded chromosome and genome sizes, which are about three times smaller in fish than in mammals, distorted their results and contributed to the persisting confusion about fish genome composition. Chromosome size and their transposons may drive the AT/GC heterogeneity apparent on mammalian chromosomes, whereas far less in fishes.

Cytogenetics Meets Genomics: Cytotaxonomy and Genomic Relationships among Color Variants of the Asian Arowana Scleropages formosus

Scleropages formosus (Osteoglossiformes, Teleostei) represents one of the most valued ornamental fishes, yet it is critically endangered due to overexploitation and habitat destruction. This species encompasses three major color groups that naturally occur in allopatric populations, but the evolutionary and taxonomic relationships of S. formosus color varieties remain uncertain. Here, we utilized a range of molecular cytogenetic techniques to characterize the karyotypes of five S. formosus color phenotypes, which correspond to naturally occurring variants: the red ones (Super Red); the golden ones (Golden Crossback and Highback Golden); the green ones (Asian Green and Yellow Tail Silver). Additionally, we describe the satellitome of S. formosus (Highback Golden) by applying a high-throughput sequencing technology. All color phenotypes possessed the same karyotype structure 2n = 50 (8m/sm + 42st/a) and distribution of SatDNAs, but different chromosomal locations of rDNAs, which were involved in a chromosome size polymorphism. Our results show indications of population genetic structure and microstructure differences in karyotypes of the color phenotypes. However, the findings do not clearly back up the hypothesis that there are discrete lineages or evolutionary units among the color phenotypes of S. formosus, but another case of interspecific chromosome stasis cannot be excluded.

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.

Advances in Vertebrate (Cyto)Genomics Shed New Light on Fish Compositional Genome Evolution

Cytogenetic and compositional studies considered fish genomes rather poor in guanine-cytosine content (GC%) because of a putative "sharp increase in genic GC% during the evolution of higher vertebrates". However, the available genomic data have not been exploited to confirm this viewpoint. In contrast, further misunderstandings in GC%, mostly of fish genomes, originated from a misapprehension of the current flood of data. Utilizing public databases, we calculated the GC% in animal genomes of three different, technically well-established fractions: DNA (entire genome), cDNA (complementary DNA), and cds (exons). Our results across chordates help set borders of GC% values that are still incorrect in literature and show: (i) fish in their immense diversity possess comparably GC-rich (or even GC-richer) genomes as higher vertebrates, and fish exons are GC-enriched among vertebrates; (ii) animal genomes generally show a GC-enrichment from the DNA, over cDNA, to the cds level (i.e., not only the higher vertebrates); (iii) fish and invertebrates show a broad(er) inter-quartile range in GC%, while avian and mammalian genomes are more constrained in their GC%. These results indicate no sharp increase in the GC% of genes during the transition to higher vertebrates, as stated and numerously repeated before. We present our results in 2D and 3D space to explore the compositional genome landscape and prepared an online platform to explore the AT/GC compositional genome evolution.

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. 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 I)'.

Quantitative Approach to Fish Cytogenetics in the Context of Vertebrate Genome Evolution

Our novel Python-based tool EVANGELIST allows the visualization of GC and repeats percentages along chromosomes in sequenced genomes and has enabled us to perform quantitative large-scale analyses on the chromosome level in fish and other vertebrates. This is a different approach from the prevailing analyses, i.e., analyses of GC% in the coding sequences that make up not more than 2% in human. We identified GC content (GC%) elevations in microchromosomes in ancient fish lineages similar to avian microchromosomes and a large variability in the relationship between the chromosome size and their GC% across fish lineages. This raises the question as to what extent does the chromosome size drive GC% as posited by the currently accepted explanation based on the recombination rate. We ascribe the differences found across fishes to varying GC% of repetitive sequences. Generally, our results suggest that the GC% of repeats and proportion of repeats are independent of the chromosome size. This leaves an open space for another mechanism driving the GC evolution in vertebrates.

GC and Repeats Profiling along Chromosomes-The Future of Fish Compositional Cytogenomics

The study of fish cytogenetics has been impeded by the inability to produce G-bands that could assign chromosomes to their homologous pairs. Thus, the majority of karyotypes published have been estimated based on morphological similarities of chromosomes. The reason why chromosome G-banding does not work in fish remains elusive. However, the recent increase in the number of fish genomes assembled to the chromosome level provides a way to analyse this issue. We have developed a Python tool to visualize and quantify GC percentage (GC%) of both repeats and unique DNA along chromosomes using a non-overlapping sliding window approach. Our tool profiles GC% and simultaneously plots the proportion of repeats (rep%) in a color scale (or vice versa). Hence, it is possible to assess the contribution of repeats to the total GC%. The main differences are the GC% of repeats homogenizing the overall GC% along fish chromosomes and a greater range of GC% scattered along fish chromosomes. This may explain the inability to produce G-banding in fish. We also show an occasional banding pattern along the chromosomes in some fish that probably cannot be detected with traditional qualitative cytogenetic methods.

Karyotype Characterisation of Two Australian Dragon Lizards (Squamata: Agamidae: Amphibolurinae) Reveals Subtle Chromosomal Rearrangements Between Related Species with Similar Karyotypes

Agamid lizards (Squamata: Agamidae) are karyotypically heterogeneous. Among the 101 species currently described from Australia, all are from the subfamily Amphibolurinae. This group is, with some exceptions, karyotypically conserved, and all species involving heterogametic sex show female heterogamety. Here, we describe the chromosomes of 2 additional Australian agamid lizards, Tympanocryptis lineata and Rankinia diemensis. These species are phylogenetically and cytogenetically sisters to the well-characterised Pogona vitticeps, but their sex chromosomes and other chromosomal characteristics are unknown. In this study, we applied advanced molecular cytogenetic techniques, such as fluorescence in situ hybridisation (FISH) and cross-species gene mapping, to characterise chromosomes and to identify sex chromosomes in these species. Our data suggest that both species have a conserved karyotype with P. vitticeps but with subtle rearrangements in the chromosomal landscapes. We could identify that T. lineata possesses a female heterogametic system (ZZ/ZW) with a pair of sex microchromosomes, while R. diemensis may have heterogametic sex chromosomes, but this requires further investigations. Our study shows the pattern of chromosomal rearrangements between closely related species, explaining the speciation within Australian agamid lizards of similar karyotypes.

Artificial whole genome duplication in paleopolyploid sturgeons yields highest documented chromosome number in vertebrates

Critically endangered sturgeons, having undergone three whole genome duplication events, represent an exceptional example of ploidy plasticity in vertebrates. Three extant ploidy groups, combined with autopolyploidization, interspecific hybridization and the fertility of hybrids are important issues in sturgeon conservation and aquaculture. Here we demonstrate that the sturgeon genome can undergo numerous alterations of ploidy without severe physiological consequences, producing progeny with a range of ploidy levels and extremely high chromosome numbers. Artificial suppression of the first mitotic division alone, or in combination with suppression of the second meiotic division of functionally tetraploid zygotes (4n, C-value = 4.15) of Siberian sturgeon Acipenser baerii and Russian sturgeon A. gueldenstaedtii resulted in progeny of various ploidy levels—diploid/hexaploid (2n/6n) mosaics, hexaploid, octoploid juveniles (8n), and dodecaploid (12n) larvae. Counts between 477 to 520 chromosomes in octoploid juveniles of both sturgeons confirmed the modal chromosome numbers of parental species had been doubled. This exceeds the highest previously documented chromosome count among vertebrates 2n ~ 446 in the cyprinid fish Ptychobarbus dipogon.

Molecular Analysis and Chromosome Mapping of Repetitive DNAs in the Green Terror Andinoacara rivulatus (Cichlidae: Cichlasomatini)

Neotropical cichlids include hundreds of species whose taxonomy has benefited of molecular phylogeny and whose karyotype evolution has been related to the amount and distribution of different classes of repetitive sequences. This study provides the first integrative molecular (cytochrome c oxidase subunit 1 and 16S sequences) and cytogenetic analyses of wild samples of the green terror Andinoacara rivulatus, a cichlid naturally distributed in Ecuador and spread throughout the world as an aquarium pet. Molecular data revealed that sequences of green terror constitute a single monophyletic clade within the genus and allowed species attribution of uncertain samples previously cytogenetically analyzed. Chromosome number (2n = 48) conforms to the general trend observed within neotropical cichlids. However, mapping of different classes of repeated sequences (18S rDNA, 5S rDNA, U1 snDNA and telomeric) revealed the presence of features uncommon among representatives of these fishes, like multiple major rDNA sites, and suggested a recent occurrence of rearrangements (fusion/inversion) in two chromosome pairs.

Range distribution and contributions to taxonomy of Elops smithi (ELOPIFORMES: ELOPIDAE)

This study presents the first record of Elops smithi for northern Brazil. The evidence suggests this species is being misidentified incorrectly as Elops saurus in estuaries of the Western Atlantic Ocean. Here, morphological, molecular, and cytogenetic evidence identified all ladyfish specimens from one estuary in the region as E. smithi. Thus, at least Elops smithi occurs in the northern coast of Brazil and it is recommended that specimens from this region identified as E. saurus be further investigated with genetic and cytogenetic tools in order to assure a correct species identification.

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.

Cytogenetics, genomics and biodiversity of the South American and African Arapaimidae fish family (Teleostei, Osteoglossiformes)

Osteoglossiformes represents one of the most ancestral teleost lineages, currently widespread over almost all continents, except for Antarctica. However, data involving advanced molecular cytogenetics or comparative genomics are yet largely limited for this fish group. Therefore, the present investigations focus on the osteoglossiform family Arapaimidae, studying a unique fish model group with advanced molecular cytogenetic genomic tools. The aim is to better explore and clarify certain events and factors that had impact on evolutionary history of this fish group. For that, both South American and African representatives of Arapaimidae, namely Arapaima gigas and Heterotis niloticus, were examined. Both species differed markedly by diploid chromosome numbers, with 2n = 56 found in A. gigas and 2n = 40 exhibited by H. niloticus. Conventional cytogenetics along with fluorescence in situ hybridization revealed some general trends shared by most osteoglossiform species analyzed thus far, such as the presence of only one chromosome pair bearing 18S and 5S rDNA sites and karyotypes dominated by acrocentric chromosomes, resembling thus the patterns of hypothetical ancestral teleost karyotype. Furthermore, the genomes of A. gigas and H. niloticus display remarkable divergence in terms of repetitive DNA content and distribution, as revealed by comparative genomic hybridization (CGH). On the other hand, genomic diversity of single copy sequences studied through principal component analyses (PCA) based on SNP alleles genotyped by the DArT seq procedure demonstrated a very low genetic distance between the South American and African Arapaimidae species; this pattern contrasts sharply with the scenario found in other osteoglossiform species. Underlying evolutionary mechanisms potentially explaining the obtained data have been suggested and discussed.

Karyotype and chromosomal characteristics of rDNA of Cobitisstrumicae Karaman, 1955 (Teleostei, Cobitidae) from Lake Volvi, Greece

The karyotype of Greek cobitid fish Cobitisstrumicae Karaman, 1955, from Lake Volvi, Greece, a representative of one of its two major intraspecific phylogenetic lineages, was analysed by means of sequential Giemsa-staining, C-banding, silver-staining, CMA3 fluorescence banding and also by in situ hybridization (FISH) with rDNA probe. The diploid chromosome number was 2n = 50, karyotype composed of 10 pairs of metacentric to submetacentric and 15 pairs of subtelocentric to acrocentric chromosomes. The nucleolus organizer regions (NORs) as revealed by Ag- and CMA3 staining and FISH were situated in the telomeric region of the fourth submetacentric chromosome pair. The chromosomes contained very low content of C-positive heterochromatin. No heteromorphic sex chromosomes were detected. This first karyotype report for any species of lineage Bicanestrinia Băcescu, 1962 shows a simple karyotype dominated by acrocentric chromosomes and possessing single NOR-bearing chromosome pair. Cytotaxonomic implications of this finding for the taxonomy of the genus Cobitis Linnaeus, 1758 are further discussed.

From Chromosomes to Genome: Insights into the Evolutionary Relationships and Biogeography of Old World Knifefishes (Notopteridae; Osteoglossiformes)

In addition to its wide geographical distribution, osteoglossiform fishes represent one of the most ancient freshwater teleost lineages; making it an important group for systematic and evolutionary studies. These fishes had a Gondwanan origin and their past distribution may have contributed to the diversity present in this group. However, cytogenetic and genomic data are still scarce, making it difficult to track evolutionary trajectories within this order. In addition, their wide distribution, with groups endemic to different continents, hinders an integrative study that allows a globalized view of its evolutionary process. Here, we performed a detailed chromosomal analysis in Notopteridae fishes, using conventional and advanced molecular cytogenetic methods. Moreover, the genetic distances of examined species were assessed by genotyping using diversity arrays technology sequencing (DArTseq). These data provided a clear picture of the genetic diversity between African and Asian Notopteridae species, and were highly consistent with the chromosomal, geographical, and historical data, enlightening their evolutionary diversification. Here, we discuss the impact of continental drift and split of Pangea on their recent diversity, as well as the contribution to biogeographical models that explain their distribution, highlighting the role of the Indian subcontinent in the evolutionary process within the family.

Neopterygian phylogeny: the merger assay

The phylogenetic relationships of the recently described genus †Ticinolepis from the Middle Triassic of the Monte San Giorgio are explored through cladistic analyses of the so far largest morphological dataset for fossil actinopterygians, including representatives of the crown-neopterygian clades Halecomorphi, Ginglymodi and Teleostei, and merging the characters from previously published systematic studies together with newly proposed characters. †Ticinolepis is retrieved as the most basal Ginglymodi and our results support the monophyly of Teleostei and Holostei, as well as Halecomorphi and Ginglymodi within the latter clade. The patterns of relationships within these clades mostly agree with those of previous studies, although a few important differences require future research. According to our results, ionoscopiforms are not monophyletic, caturids are not amiiforms and leptolepids and luisiellids form a monophyletic clade. Our phylogenetic hypothesis confirms the rapid radiation of the holostean clades Halecomorphi and Ginglymodi during the Early and Middle Triassic and the radiation of pholidophoriform teleosts during the Late Triassic. Crown-group Halecomorphi have an enormous ghost lineage throughout half of the Mesozoic, but ginglymodians and teleosts show a second radiation during the Early Jurassic. The crown-groups of Halecomorphi, Ginglymodi and Teleostei originated within parallel events of radiation during the Late Jurassic.

Vertebrate Genome Evolution in the Light of Fish Cytogenomics and rDNAomics

To understand the cytogenomic evolution of vertebrates, we must first unravel the complex genomes of fishes, which were the first vertebrates to evolve and were ancestors to all other vertebrates. We must not forget the immense time span during which the fish genomes had to evolve. Fish cytogenomics is endowed with unique features which offer irreplaceable insights into the evolution of the vertebrate genome. Due to the general DNA base compositional homogeneity of fish genomes, fish cytogenomics is largely based on mapping DNA repeats that still represent serious obstacles in genome sequencing and assembling, even in model species. Localization of repeats on chromosomes of hundreds of fish species and populations originating from diversified environments have revealed the biological importance of this genomic fraction. Ribosomal genes (rDNA) belong to the most informative repeats and in fish, they are subject to a more relaxed regulation than in higher vertebrates. This can result in formation of a literal 'rDNAome' consisting of more than 20,000 copies with their high proportion employed in extra-coding functions. Because rDNA has high rates of transcription and recombination, it contributes to genome diversification and can form reproductive barrier. Our overall knowledge of fish cytogenomics grows rapidly by a continuously increasing number of fish genomes sequenced and by use of novel sequencing methods improving genome assembly. The recently revealed exceptional compositional heterogeneity in an ancient fish lineage (gars) sheds new light on the compositional genome evolution in vertebrates generally. We highlight the power of synergy of cytogenetics and genomics in fish cytogenomics, its potential to understand the complexity of genome evolution in vertebrates, which is also linked to clinical applications and the chromosomal backgrounds of speciation. We also summarize the current knowledge on fish cytogenomics and outline its main future avenues.