The two giant sister species of the Southern Ocean, Dissostichus eleginoides and Dissostichus mawsoni, differ in karyotype and chromosomal pattern of ribosomal RNA genes

Chromosome-level genome assembly and annotation of the Patagonian toothfish Dissostichus eleginoides

The Patagonian toothfish (Dissostichus eleginoides) belongs to the Actinopterygii class, and the suborder Notothenioidei, which lives in cold waters in the Southern Hemisphere. We performed assembly and annotation, and we integrated the Illumina short-read sequencing for polishinng, PacBio long-read sequencing for contig-level assembly, and Hi-C sequencing technology to obtain high-quality of chromosome-level genome assembly. The final assembly analysis resulted in a total of 495 scaffolds, a genome size of 844.7 Mbp and an N50 length of 36 Mbp. Among these data, we confirmed 24 scaffolds exceeded 10 Mbp and classified as chromosome-level. The completeness of BUSCO rate was over 97%. A total gene set of 32,224 was identified. Furthermore, we analyzed the presence of AFGP genes, classified into Antarctic and sub-Antarctic categories through phylogenetic analysis. This study provides a useful resource for the genomic analysis of Patagonian toothfish and genetic insights into the comparison with Antarctic fishes.

De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome

BACKGROUND

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited.

RESULTS

A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome.

CONCLUSIONS

We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Single Copies of the 5S rRNA Inserted into 45S rDNA Intergenic Spacers in the Genomes of Nototheniidae (Perciformes, Actinopterygii)

In the vast majority of Animalia genomes, the 5S rRNA gene repeats are located on chromosomes outside of the 45S rDNA arrays of the nucleolar organiser (NOR). We analysed the genomic databases available and found that a 5S rDNA sequence is inserted into the intergenic spacer (IGS) between the 45S rDNA repeats in ten species of the family Nototheniidae (Perciformes, Actinopterigii). We call this sequence the NOR-5S rRNA gene. Along with Testudines and Crocodilia, this is the second case of a close association between four rRNA genes within one repetitive unit in deuterostomes. In both cases, NOR-5S is oriented opposite the 45S rDNA. None of the three nucleotide substitutions compared to the canonical 5S rRNA gene influenced the 5S rRNA secondary structure. In transcriptomes of the Patagonian toothfish, we only found NOR-5S rRNA reads in ovaries and early embryos, but not in testis or somatic tissues of adults. Thus, we consider the NOR-5S gene to be a maternal-type 5S rRNA template. The colocalization of the 5S and 45S ribosomal genes appears to be essential for the equimolar production of all four rRNAs in the species that show rDNA amplification during oogenesis. Most likely, the integration of 5S and NOR rRNA genes occurred prior to Nototheniidae lineage diversification.

Physical Mapping of Repeated Sequences on Fish Chromosomes by Fluorescence In Situ Hybridization (FISH)

The opportunity to map genes and noncoding DNA sequences on the chromosomes by fluorescence in situ hybridization (FISH) has greatly enhanced the potential for fish karyotyping and comparative cytogenetics. The use of FISH allowed for significant advances in our understanding of the fish genome architecture, especially when applied to the study of the repetitive component of the genome, that is generally underestimated in the bioinformatic assembly. Here we provide a step-by-step protocol for FISH of repeated sequences onto chromosomes of fish species.

Comparative study of four Mystus species (Bagridae, Siluriformes) from Thailand: insights into their karyotypic diversity

Karyotypes of four catfishes of the genus Mystus Scopoli, 1777 (family Bagridae), M. atrifasciatus Fowler, 1937, M. mysticetus Roberts, 1992, M. singaringan (Bleeker, 1846) and M. wolffii (Bleeker, 1851), were analysed by conventional and Ag-NOR banding as well as fluorescence in situ hybridization (FISH) techniques. Microsatellite d(GC)₁₅, d(CAA)₁₀, d(CAT)₁₀ and d(GAA)₁₀ repeat probes were applied in FISH. The obtained data revealed that the four studied species have different chromosome complements. The diploid chromosome numbers (2n) and the fundamental numbers (NF) range between 52 and 102, 54 and 104, 56 and 98, or 58 and 108 in M. mysticetus, M. atrifasciatus, M. singaringan or M. wolffii, respectively. Karyotype formulae of M. mysticetus, M. atrifasciatus, M. singaringan and M. wolffii are 24m+26sm+4a, 26m+24sm+2a, 24m+18sm+14a and 30m+22sm+6a, respectively. A single pair of NORs was identified adjacent to the telomeres of the short arm of chromosome pairs 3 (metacentric) in M. atrifasciatus, 20 (submetacentric) in M. mysticetus, 15 (submetacentric) in M. singaringan, and 5 (metacentric) in M. wolffii. The d(GC)₁₅, d(CAA)₁₀, d(CAT)₁₀ and d(GAA)₁₀ repeats were abundantly distributed in species-specific patterns. Overall, we present a comparison of cytogenetic and molecular cytogenetic patterns of four species from genus Mystus providing insights into their karyotype diversity in the genus.

Latitudinal Cline in Chromosome Numbers of Ice Cod A. glacialis (Gadidae) from Northeast Greenland

The ice cod Arctogadus glacialis (Peters, 1872) is one of the few fish species endemic to the Arctic. With a circumpolar distribution, the species is confined to the fjords and shelves of the Arctic seas. Biological information on A. glacialis is scarce, with genomic information restricted to microsatellites. Within the frame of the TUNU-Programme: Arctic Ocean Fishes-Diversity, Adaptation and Conservation, we studied A. glacialis at the chromosomal level to explore fish diversity and evolutionary aspects. The analysis of over 50 individuals from the Northeast Greenland fjords between latitudes 71°09' N and 76°42' N revealed a remarkable intraspecific diversity epitomized by chromosome numbers spanning from 28 to 33, the occurrence of putative B chromosomes, and diversified patterns of distribution of heterochromatin and rDNAs. The number of B chromosomes followed a latitudinal gradient from 0-2 in the north to 2-5 in the south. Considering the benthic and rather stationary life history of this species, the observed chromosomal differences might have arisen independently, possibly driven and/or fostered by the dynamics of repetitive sequences, and are being fixed in relatively isolated fjord populations. The resulting latitudinal cline we observe today might have repercussions on the fate of local populations facing the ongoing climate-driven environmental changes.

Cytogenetic characterization of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902) through analysis of mitotic chromosomes from early larvae

This paper describes the cytogenetic features of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902), a keystone species of the Antarctic coastal marine ecosystem. Conventional cytogenetic analyses and physical mapping of repetitive DNA sequences were performed on metaphase plates obtained through direct chromosome preparation from P. antarctica early larvae. The Antarctic silverfish have a diploid number (2n) = 48, and a karyotype made up of a majority of two-armed chromosomes (karyotype formula36m/sm + 10st + 2a, fundamental number = 94). Major ribosomal gene repeats were detected on three chromosome pairs (20, 21, and 23), in correspondence of dim DAPI stained regions. Long Interspersed Nuclear Elements (LINEs) were abundant and wide spread over all chromosomes. Overall, the cytogenetic data presented herein are consistent with a long independent cytogenetic and evolutionary history for the species. The large number of two-armed chromosomes, indicative of highly-rearranged karyotype, coupled with a diploid number of 48, a presumed primitive character for this fish group, and the spread of the major ribosomal genes on three chromosome pairs, make the Antarctic silverfish distinct from all other notothenioid species.

The genomic basis for colonizing the freezing Southern Ocean revealed by Antarctic toothfish and Patagonian robalo genomes

Background

The Southern Ocean is the coldest ocean on Earth but a hot spot of evolution. The bottom-dwelling Eocene ancestor of Antarctic notothenioid fishes survived polar marine glaciation and underwent adaptive radiation, forming >120 species that fill all water column niches today. Genome-wide changes enabling physiological adaptations and the rapid expansion of the Antarctic notothenioids remain poorly understood.

Results

We sequenced and compared 2 notothenioid genomes—the cold-adapted and neutrally buoyant Antarctic toothfish Dissostichus mawsoni and the basal Patagonian robalo Eleginops maclovinus, representing the temperate ancestor. We detected >200 protein gene families that had expanded and thousands of genes that had evolved faster in the toothfish, with diverse cold-relevant functions including stress response, lipid metabolism, protein homeostasis, and freeze resistance. Besides antifreeze glycoprotein, an eggshell protein had functionally diversified to aid in cellular freezing resistance. Genomic and transcriptomic comparisons revealed proliferation of selcys–transfer RNA genes and broad transcriptional upregulation across anti-oxidative selenoproteins, signifying their prominent role in mitigating oxidative stress in the oxygen-rich Southern Ocean. We found expansion of transposable elements, temporally correlated to Antarctic notothenioid diversification. Additionally, the toothfish exhibited remarkable shifts in genetic programs towards enhanced fat cell differentiation and lipid storage, and promotion of chondrogenesis while inhibiting osteogenesis in bone development, collectively contributing to the achievement of neutral buoyancy and pelagicism.

Conclusions

Our study revealed a comprehensive landscape of evolutionary changes essential for Antarctic notothenioid cold adaptation and ecological expansion. The 2 genomes are valuable resources for further exploration of mechanisms underlying the spectacular notothenioid radiation in the coldest marine environment.

Cold Fusion: Massive Karyotype Evolution in the Antarctic Bullhead Notothen Notothenia coriiceps

Half of all vertebrate species share a series of chromosome fusions that preceded the teleost genome duplication (TGD), but we do not understand the causative evolutionary mechanisms. The "Robertsonian-translocation hypothesis" suggests a regular fusion of each ancestral acro- or telocentric chromosome to just one other by centromere fusions, thus halving the karyotype. An alternative "genome-stirring hypothesis" posits haphazard and repeated fusions, inversions, and reciprocal and nonreciprocal translocations. To study large-scale karyotype reduction, we investigated the decrease of chromosome numbers in Antarctic notothenioid fish. Most notothenioids have 24 haploid chromosomes, but bullhead notothen (Notothenia coriiceps) has 11. To understand mechanisms, we made a RAD-tag meiotic map with ∼10,000 polymorphic markers. Comparative genomics aligned about a thousand orthologs of platyfish and stickleback genes along bullhead chromosomes. Results revealed that 9 of 11 bullhead chromosomes arose by fusion of just two ancestral chromosomes and two others by fusion of three ancestral chromosomes. All markers from each ancestral chromosome remained contiguous, implying no inversions across fusion borders. Karyotype comparisons support a history of: (1) Robertsonian fusions of 22 ancestral chromosomes in pairs to yield 11 fused plus two small unfused chromosomes, like N. angustata; (2) fusion of one of the remaining two ancestral chromosomes to a preexisting fused pair, giving 12 chromosomes like N. rossii; and (3) fusion of the remaining ancestral chromosome to another fused pair, giving 11 chromosomes in N. coriiceps These results raise the question of what selective forces promoted the systematic fusion of chromosomes in pairs and the suppression of pericentric inversions in this lineage, and provide a model for chromosome fusions in stem teleosts.

Co-located 18S/5S rDNA arrays: an ancient and unusual chromosomal trait in Julidini species (Labridae, Perciformes)

Wrasses (Labridae) are extremely diversified marine fishes, whose species exhibit complex interactions with the reef environment. They are widely distributed in the Indian, Pacific and Atlantic oceans. Their species have displayed a number of karyotypic divergent processes, including chromosomal regions with complex structural organization. Current cytogenetic information for this family is phylogenetically and geographically limited and mainly based on conventional cytogenetic techniques. Here, the distribution patterns of heterochromatin, GC-specific chromosome regions and Ag-NORs, and the organization of 18S and 5S rDNA sites of the Atlantic species Thalassoma noronhanum (Boulenger, 1890), Halichoeres poeyi (Steindachner, 1867), Halichoeres radiatus (Linnaeus, 1758), Halichoeres brasiliensis (Bloch, 1791) and Halichoeres penrosei Starks, 1913, belonging to the tribe Julidini were analyzed. All the species exhibited 2n=48 chromosomes with variation in the number of chromosome arms among genera. Thalassoma noronhanum has 2m+46a, while species of the genus Halichoeres Rüppell, 1835 share karyotypes with 48 acrocentric chromosomes. The Halichoeres species exhibit differences in the heterochromatin distribution patterns and in the number and distribution of 18S and 5S rDNA sites. The occurrence of 18S/5S rDNA syntenic arrangements in all the species indicates a functionally stable and adaptive genomic organization. The phylogenetic sharing of this rDNA organization highlights a marked and unusual chromosomal singularity inside the family Labridae.

Molecular cytogenetic analyses of Epinephelus bruneus and Epinephelus moara (Perciformes, Epinephelidae)

Genus Epinephelus (Perciformes, Epinephelidae), commonly known as groupers, are usually difficult in species identification for the lack and/or change of morphological specialization. In this study, molecular cytogenetic analyses were firstly performed to identify the closely related species Epinephelus bruneus and E. moara in this genus. The species-specific differences of both fish species showed in karyotype, chromosomal distribution of nucleolar organizer regions (NORs) and localization of 18S rDNA. The heterochromatin (interstitial C-bands) and distribution pattern of telomere (TTAGGG)_n in E. bruneus revealed the chromosomal rearrangements and different karyotypic evolutionary characteristics compared to those in E. moara. The cytogenetic data suggested that the lineages of E. bruneus and E. moara were recently derived within the genus Epinephelus, and E. moara exhibited more plesiomorphic features than E. bruneus. All results confirmed that E. moara, which has long been considered a synonym of E. bruneus, is a distinct species in the family Epinephelidae. In addition, molecular cytogenetic analyses are useful in species differentiation and phylogenetic reconstruction in groupers.

Twenty years of physical mapping of major ribosomal RNA genes across the teleosts: A review of research

Molecular cytogenetic data on the number and position of 45S ribosomal RNA genes (rDNA; located in nucleolus organizing regions, NORs) detected by FISH in 330 species of 77 families and 22 orders of bony fishes (Teleostei) and, additionally, 11 species of basal ray-finned fishes are compiled and analyzed. The portion of species with single rDNA sites in the sample amounts to 72%. The percentage of species with multiple NORs decreases with increasing numbers of rDNA loci per genome, i.e. scarcely 3% of species carry 4 or more rDNA-bearing chromosome pairs. 43% of all rDNA sites analyzed occur terminally on the short arms of chromosomes or constitute them. In general, terminal rDNA sites account for 87% of all examined cases. Interspecific variation in the location of single rDNA sites among related taxa, polymorphisms of multiple NORs in some groups of teleosts and analytical outcomes on the subject are reviewed.

Karyotyping and cytogenetic mapping of Atlantic cod (Gadus morhua Linnaeus, 1758)

The Atlantic cod (Gadus morhua) is an important natural resource for northern societies and is now also considered to be a promising candidate for aquaculture. In recent years, much effort has been directed towards the development of genomic tools, and genome initiatives for Atlantic cod have been established. Despite the growing attention devoted to the Atlantic cod genomics, basic aspects of its genome structure and organization remain unknown. Thus, the present work aims to study cytogenetic features of the Atlantic cod as a contribution to the knowledge of this species' genome. The Atlantic cod displays a diploid number of 46 chromosomes, with a karyotypic formula 16 m/sm + 30 st/t. Conventional karyotyping was improved by chromosomal mapping of two classes of repetitive sequences. 18S rDNA clusters were assigned to pairs 2 and 4; small amounts of 18S rDNA clusters were occasionally detected on pair 5. These findings could not be related to the geographical origin of the specimens, but were consistent with the variability of these repeated genes in fish in general. 5S ribosomal gene clusters, apparently corresponding to a single 5S rDNA class, were detected on twelve chromosomes (pairs 11, 12, 14, 17, 20 and 21). The present update of the existing but meagre information on the karyotype of Atlantic cod, plus the first physical mapping of repetitive genes in this species herein, opens the way for an integrated approach that combines genetic and physical mapping with the assembly of the genome of this commercially important species.

The Patagonian toothfish: biology, ecology and fishery

Patagonian toothfish (Dissostichus eleginoides) is a large notothenioid fish that supports valuable fisheries throughout the Southern Ocean. D. eleginoides are found on the southern shelves and slopes of South America and around the sub-Antarctic islands of the Southern Ocean. Patagonian toothfish are a long-lived species (>50 years), which initially grow rapidly on the shallow shelf areas, before undertaking an ontogenetic migration into deeper water. Although they are active predators and scavengers, there is no evidence of large-scale geographic migrations, and studies using genetics, biochemistry, parasite fauna and tagging indicate a high degree of isolation between populations in the Indian Ocean, South Georgia and the Patagonian Shelf. Patagonian toothfish spawn in deep water (ca. 1000 m) during the austral winter, producing pelagic eggs and larvae. Larvae switch to a demersal habitat at around 100 mm (1-year-old) and inhabit relatively shallow water (<300 m) until 6-7 years of age, when they begin a gradual migration into deeper water. As juveniles in shallow water, toothfish are primarily piscivorous, consuming the most abundant suitably sized local prey. With increasing size and habitat depth, the diet diversifies and includes more scavenging. Toothfish have weakly mineralised skeletons and a high fat content in muscle, which helps neutral buoyancy, but limits swimming capacity. Toothfish generally swim with labriform motion, but are capable of more rapid sub-carangiform swimming when startled. Toothfish were first caught as a by-catch (as juveniles) in shallow trawl fisheries, but following the development of deep water longlining, fisheries rapidly developed throughout the Southern Ocean. The initial rapid expansion of the fishery, which led to a peak of over 40,000 tonnes in reported landings in 1995, was accompanied by problems of bird by-catch and overexploitation as a consequence of illegal, unreported and unregulated fishing (IUU). These problems have now largely been addressed, but continued vigilance is required to ensure that the species is sustainably exploited and the ecosystem effects of the fisheries are minimised.