Genome Compositional Organization in Gars Shows More Similarities to Mammals than to Other Ray-Finned Fish

Cytogenetic Analysis of the Fish Genus Carassius Indicates Divergence, Fission, and Segmental Duplication as Drivers of Tandem Repeat and Microchromosome Evolution

Fishes of the genus Carassius are useful experimental vertebrate models for the study of evolutionary biology and cytogenetics. Carassius demonstrates diverse biological characteristics, such as variation in ploidy levels and chromosome numbers, and presence of microchromosomes. Those Carassius polyploids with ≥150 chromosomes have microchromosomes, but the origin of microchromosomes, especially in European populations, is unknown. We used cytogenetics to study evolution of tandem repeats (U1 and U2 small nuclear DNAs and H3 histone) and microchromosomes in Carassius from the Czech Republic. We tested the hypotheses whether the number of tandem repeats was affected by polyploidization or divergence between species and what mechanism drives evolution of microchromosomes. Tandem repeats were found in tetraploid and hexaploid Carassius gibelio, and tetraploid Carassius auratus and Carassius carassius in conserved numbers, with the exception of U1 small nuclear DNA in C. auratus. This conservation indicates reduction and/or loss in the number of copies per locus in hexaploids and may have occurred by divergence rather than polyploidization. To study the evolution of microchromosomes, we used the whole microchromosome painting probe from hexaploid C. gibelio and hybridized it to tetraploid and hexaploid C. gibelio, and tetraploid C. auratus and C. carassius. Our results revealed variation in the number of microchromosomes in hexaploids and indicated that the evolution of the Carassius karyotype is governed by macrochromosome fissions followed by segmental duplication in pericentromeric areas. These are potential mechanisms responsible for the presence of microchromosomes in Carassius hexaploids. Differential efficacy of one or both of these mechanisms in different tetraploids could ensure variability in chromosome number in polyploids in general.

Post-testicular sperm maturation in ancient holostean species

Fish speciation was accompanied by changes in the urogenital system anatomy. In evolutionarily modern Teleostei, male reproductive tracts are fully separated from the excretory system, while in evolutionarily ancient Chondrostei and Holostei, the excretory and reproductive tracts are not separated. Sturgeon post-testicular sperm maturation (PTSM) occurring as a result of sperm/urine mixing is phenomenologically well described, while, in holosteans, functional intimacy of seminal ducts with kidney ducts and the existence of PTSM still need to be addressed. In Lepisosteus platostomus (Holostei), sperm samples were collected from testes (TS), efferent ducts (EDS), and Wolffian ducts (WDS). While WDS was motile, no motility was found in TS and EDS. The existence of PTSM was checked by in vitro PTSM procedure. After TS and EDS incubation in seminal fluid from WDS, no more than 5% motile spermatozoa were observed in TS, whereas in EDS the motility percentage was up to 75%. Experimental dyeing of urogenital ducts in gars and sturgeons revealed some differences in the interconnection between sperm ducts and kidneys. It is concluded that post-testicular sperm maturation occurs in gars and suggests that infraclass Holostei occupies an intermediate evolutionary position between Teleostei and Chondrostei in the anatomical arrangement of the urogenital system.

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.

Integrating Genomic and Chromosomal Data: A Cytogenetic Study of Transancistrus santarosensis (Loricariidae: Hypostominae) with Characterization of a ZZ/ZW Sex Chromosome System

The plecos (Loricariidae) fish represent a great model for cytogenetic investigations due to their variety of karyotypes, including diploid and polyploid genomes, and different types of sex chromosomes. In this study we investigate Transancistrus santarosensis a rare loricariid endemic to Ecuador, integrating cytogenetic methods with specimens' molecular identification by mtDNA, to describe the the species karyotype. We aim to verify whether sex chromosomes are cytologically identifiable and if they are associated with the accumulation of repetitive sequences present in other species of the family. The analysis of the karyotype (2n = 54 chromosomes) excludes recent centric fusion and pericentromeric inversion and suggests the presence of a ZZ/ZW sex chromosome system at an early stage of differentiation: the W chromosome is degenerated but is not characterized by the presence of differential sex-specific repetitive DNAs. Data indicate that although T. santarosensis has retained the ancestral diploid number of Loricariidae, it accumulated heterochromatin and shows non-syntenic ribosomal genes localization, chromosomal traits considered apomorphic in the family.

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.

Expression Pattern of nos1 in the Developing Nervous System of Ray-Finned Fish

Fish have colonized nearly all aquatic niches, making them an invaluable resource to understand vertebrate adaptation and gene family evolution, including the evolution of complex neural networks and modulatory neurotransmitter pathways. Among ancient regulatory molecules, the gaseous messenger nitric oxide (NO) is involved in a wide range of biological processes. Because of its short half-life, the modulatory capability of NO is strictly related to the local activity of nitric oxide synthases (Nos), enzymes that synthesize NO from L-arginine, making the localization of Nos mRNAs a reliable indirect proxy for the location of NO action domains, targets, and effectors. Within the diversified actinopterygian nos paralogs, nos1 (alias nnos) is ubiquitously present as a single copy gene across the gnathostome lineage, making it an ideal candidate for comparative studies. To investigate variations in the NO system across ray-finned fish phylogeny, we compared nos1 expression patterns during the development of two well-established experimental teleosts (zebrafish and medaka) with an early branching holostean (spotted gar), an important evolutionary bridge between teleosts and tetrapods. Data reported here highlight both conserved expression domains and species-specific nos1 territories, confirming the ancestry of this signaling system and expanding the number of biological processes implicated in NO activities.

On the Base Composition of Transposable Elements

Transposable elements exhibit a base composition that is often different from the genomic average and from hosts’ genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.

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.

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.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Nucleotide composition of transposable elements likely contributes to AT/GC compositional homogeneity of teleost fish genomes

BACKGROUND

Teleost fish genome size has been repeatedly demonstrated to positively correlate with the proportion of transposable elements (TEs). This finding might have far-reaching implications for our understanding of the evolution of nucleotide composition across vertebrates. Genomes of fish and amphibians are GC homogenous, with non-teleost gars being the single exception identified to date, whereas birds and mammals are AT/GC heterogeneous. The exact reason for this phenomenon remains controversial. Since TEs make up significant proportions of genomes and can quickly accumulate across genomes, they can potentially influence the host genome with their own GC content (GC%). However, the GC% of fish TEs has so far been neglected.

RESULTS

The genomic proportion of TEs indeed correlates with genome size, although not as linearly as previously shown with fewer genomes, and GC% negatively correlates with genome size in the 33 fish genome assemblies analysed here (excluding salmonids). GC% of fish TE consensus sequences positively correlates with the corresponding genomic GC% in 29 species tested. Likewise, the GC contents of the entire repetitive vs. non-repetitive genomic fractions correlate positively in 54 fish species in Ensembl. However, among these fish species, there is also a wide variation in GC% between the main groups of TEs. Class II DNA transposons, predominant TEs in fish genomes, are significantly GC-poorer than Class I retrotransposons. The AT/GC heterogeneous gar genome contains fewer Class II TEs, a situation similar to fugu with its extremely compact and also GC-enriched but AT/GC homogenous genome.

CONCLUSION

Our results reveal a previously overlooked correlation between GC% of fish genomes and their TEs. This applies to both TE consensus sequences as well as the entire repetitive genomic fraction. On the other hand, there is a wide variation in GC% across fish TE groups. These results raise the question whether GC% of TEs evolves independently of GC% of the host genome or whether it is driven by TE localization in the host genome. Answering these questions will help to understand how genomic GC% is shaped over time. Long-term accumulation of GC-poor(er) Class II DNA transposons might indeed have influenced AT/GC homogenization of fish genomes and requires further investigation.

Adaptive Radiation from a Chromosomal Perspective: Evidence of Chromosome Set Stability in Cichlid Fishes (Cichlidae: Teleostei) from the Barombi Mbo Lake, Cameroon

Cichlid fishes are the subject of scientific interest because of their rapid adaptive radiation, resulting in extensive ecological and taxonomic diversity. In this study, we examined 11 morphologically distinct cichlid species endemic to Barombi Mbo, the largest crater lake in western Cameroon, namely Konia eisentrauti, Konia dikume, Myaka myaka, Pungu maclareni, Sarotherodon steinbachi, Sarotherodon lohbergeri, Sarotherodon linnellii, Sarotherodon caroli, Stomatepia mariae, Stomatepia pindu, and Stomatepia mongo. These species supposedly evolved via sympatric ecological speciation from a common ancestor, which colonized the lake no earlier than one million years ago. Here we present the first comparative cytogenetic analysis of cichlid species from Barombi Mbo Lake using both conventional (Giemsa staining, C-banding, and CMA₃/DAPI staining) and molecular (fluorescence in situ hybridization with telomeric, 5S, and 28S rDNA probes) methods. We observed stability on both macro and micro-chromosomal levels. The diploid chromosome number was 2n = 44, and the karyotype was invariably composed of three pairs of meta/submetacentric and 19 pairs of subtelo/acrocentric chromosomes in all analysed species, with the same numbers of rDNA clusters and distribution of heterochromatin. The results suggest the evolutionary stability of chromosomal set; therefore, the large-scale chromosomal rearrangements seem to be unlikely associated with the sympatric speciation in Barombi Mbo.

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.

Chromosomal Evolution and Evolutionary Relationships of Lebiasina Species (Characiformes, Lebiasinidae)

We present the first cytogenetic data for Lebiasina bimaculata and L. melanoguttata with the aim of (1) investigating evolutionary events within Lebiasina and their relationships with other Lebiasinidae genera and (2) checking the evolutionary relationships between Lebiasinidae and Ctenoluciidae. Both species have a diploid number 2n = 36 with similar karyotypes and microsatellite distribution patterns but present contrasting C-positive heterochromatin and CMA3+ banding patterns. The remarkable interstitial series of C-positive heterochromatin occurring in L. melanoguttata is absent in L. bimaculata. Accordingly, L. bimaculata shows the ribosomal DNA sites as the only GC-rich (CMA3+) regions, while L. melanoguttata shows evidence of a clear intercalated CMA3+ banding pattern. In addition, the multiple 5S and 18S rDNA sites in L. melanogutatta contrast with single sites present in L. bimaculata. Comparative genomic hybridization (CGH) experiments also revealed a high level of genomic differentiation between both species. A polymorphic state of a conspicuous C-positive, CMA3+, and (CGG)n band was found only to occur in L. bimaculata females, and its possible relationship with a nascent sex chromosome system is discussed. Whole chromosome painting (WCP) and CGH experiments indicate that the Lebiasina species examined and Boulengerella maculata share similar chromosomal sequences, thus supporting the relatedness between them and the evolutionary relationships between the Lebiasinidae and Ctenoluciidae families.

Integrative rDNAomics-Importance of the Oldest Repetitive Fraction of the Eukaryote Genome

Nuclear ribosomal RNA (rRNA) genes represent the oldest repetitive fraction universal to all eukaryotic genomes. Their deeply anchored universality and omnipresence during eukaryotic evolution reflects in multiple roles and functions reaching far beyond ribosomal synthesis. Merely the copy number of non-transcribed rRNA genes is involved in mechanisms governing e.g., maintenance of genome integrity and control of cellular aging. Their copy number can vary in response to environmental cues, in cellular stress sensing, in development of cancer and other diseases. While reaching hundreds of copies in humans, there are records of up to 20,000 copies in fish and frogs and even 400,000 copies in ciliates forming thus a literal subgenome or an rDNAome within the genome. From the compositional and evolutionary dynamics viewpoint, the precursor 45S rDNA represents universally GC-enriched, highly recombining and homogenized regions. Hence, it is not accidental that both rDNA sequence and the corresponding rRNA secondary structure belong to established phylogenetic markers broadly used to infer phylogeny on multiple taxonomical levels including species delimitation. However, these multiple roles of rDNAs have been treated and discussed as being separate and independent from each other. Here, I aim to address nuclear rDNAs in an integrative approach to better assess the complexity of rDNA importance in the evolutionary context.

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.

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.

The complete mitochondrial DNA of the Cuban gar (Atractosteus tristoechus)

The Cuban gar (Atractosteus tristoechus) is an endemic lepisosteid living in Cuba. Among gars, this species is one of the most threatened and has the smallest natural distribution range. Lepisosteids are air-breathing fishes belonging to the Holostean, a basal non-teleost clade of actinopterygians. Recent studies have indicated that these fishes could be a 'bridge between tetrapods and teleost biomedical models'. Herein, we sequenced and assembled the first complete mitochondrial genome of A. tristoechus. The total length of the mt genome is 16,290 bp, containing the typical 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a 537 bp length control region.

Sex-specific differences in transcriptome profiles of brain and muscle tissue of the tropical gar

Background

The tropical gar (Atractosteus tropicus) is the southernmost species of the seven extant species of gar fishes in the world. In Mexico and Central America, the species is an important food source due to its nutritional quality and low price. Despite its regional importance and increasing concerns about overexploitation and habitat degradation, basic genetic information on the tropical gar is lacking. Determining genetic information on the tropical gar is important for the sustainable management of wild populations, implementation of best practices in aquaculture settings, evolutionary studies of ancient lineages, and an understanding of sex-specific gene expression. In this study, the transcriptome of the tropical gar was sequenced and assembled de novo using tissues from three males and three females using Illumina sequencing technology. Sex-specific and highly differentially expressed transcripts in brain and muscle tissues between adult males and females were subsequently identified.

Results

The transcriptome was assembled de novo resulting in 80,611 transcripts with a contig N50 of 3,355 base pairs and over 168 kilobases in total length. Male muscle, brain, and gonad as well as female muscle and brain were included in the assembly. The assembled transcriptome was annotated to identify the putative function of expressed transcripts using Trinotate and SwissProt, a database of well-annotated proteins. The brain and muscle datasets were then aligned to the assembled transcriptome to identify transcripts that were differentially expressed between males and females. The contrast between male and female brain identified 109 transcripts from 106 genes that were significantly differentially expressed. In the muscle comparison, 82 transcripts from 80 genes were identified with evidence for significant differential expression. Almost all genes identified as differentially expressed were sex-specific. The differentially expressed transcripts were enriched for genes involved in cellular functioning, signaling, immune response, and tissue-specific functions.

Conclusions

This study identified differentially expressed transcripts between male and female gar in muscle and brain tissue. The majority of differentially expressed transcripts had sex-specific expression. Expanding on these findings to other developmental stages, populations, and species may lead to the identification of genetic factors contributing to the skewed sex ratio seen in the tropical gar and of sex-specific differences in expression in other species. Finally, the transcriptome assembly will open future research avenues on tropical gar development, cell function, environmental resistance, and evolution in the context of other early vertebrates.

Electronic supplementary material

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

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

Bowfin belongs to an ancient lineage of nonteleost ray-finned fishes (actinopterygians) and is the only extant survivor of a once diverged group, the Halecomorphi or Amiiformes. Owing to the scarcity of extant nonteleost ray-finned lineages, also referred as "living fossils," their phylogenetic interrelationships have been the target of multiple hypotheses concerning their sister group relationships. Molecular and morphological data sets have produced controversial results; bowfin is considered as either the sister group to genome-duplicated teleosts (together forming the group of Halecostomi) or to gars (Lepisosteiformes; together forming the group of Holostei). However, any detailed cytogenetic analysis of bowfin chromosomes has never been performed to address this issue. Here we examined bowfin chromosomes by conventional (Giemsa-staining, C-banding, base-specific fluorescence and silver staining) and molecular (FISH with rDNA probes) cytogenetic protocols. We identified diploid chromosome number 2n = 46 with a middle-sized submetacentric chromosome pair as the major ribosomal DNA-bearing (45S rDNA), GC-positive and silver-positive element. The minor rDNA (5S rDNA) sites were localized in the pericentromeric region of one middle-sized acrocentric chromosome pair. Comparison with available cytogenetic data of other nonteleost actinopterygians (bichirs, sturgeons, gars) and teleost species including representative of basally branching lineages showed bowfin chromosomal characteristics more similar to the teleost type than to any other nonteleosts. Particularly striking differences were identified between bowfin and gars, the latter of which were found to mimic mammalian AT/GC genomic organisation. Such conclusion however contradicts the most recent phylogenomic results and raises the question what states are ancestral and what are derived.