A Nanopore Based Chromosome-Level Assembly Representing Atlantic Cod from the Celtic Sea

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

Whole genome assembly and annotation of the King Angelfish (Holacanthus passer) gives insight into the evolution of marine fishes of the Tropical Eastern Pacific

Holacanthus angelfishes are some of the most iconic marine fishes of the Tropical Eastern Pacific (TEP). However, very limited genomic resources currently exist for the genus. In this study we: (i) assembled and annotated the nuclear genome of the King Angelfish (Holacanthus passer), and (ii) examined the demographic history of H. passer in the TEP. We generated 43.8 Gb of ONT and 97.3 Gb Illumina reads representing 75× and 167× coverage, respectively. The final genome assembly size was 583 Mb with a contig N50 of 5.7 Mb, which captured 97.5% of the complete Actinoterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs). Repetitive elements accounted for 5.09% of the genome, and 33,889 protein-coding genes were predicted, of which 22,984 were functionally annotated. Our demographic analysis suggests that population expansions of H. passer occurred prior to the last glacial maximum (LGM) and were more likely shaped by events associated with the closure of the Isthmus of Panama. This result is surprising, given that most rapid population expansions in both freshwater and marine organisms have been reported to occur globally after the LGM. Overall, this annotated genome assembly provides a novel molecular resource to study the evolution of Holacanthus angelfishes, while facilitating research into local adaptation, speciation, and introgression in marine fishes.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

Improved high-quality reference genome of red drum facilitates the processes of resistance-related gene exploration

Sciaenops ocellatus is among the most important artificially introduced farmed fish across 11 countries and regions. However, the frequent occurrence of extreme weather events and breeding escapes have placed great pressure on local marine biodiversity and ecosystems. We reported the de novo assembly and annotation with a contig N50 of 28.30 Mb using PacBio HiFi sequencing and Hi-C technologies, which resulted in a 283-fold increase in contig N50 length and improvement in continuity and quality in complex repetitive region for S. ocellatus compared to the previous version. In total, 257.36 Mb of repetitive sequences accounted for 35.48% of the genome, and 22,845 protein-coding genes associated with a BUSCO value of 98.32%, were identified by genome annotation. Moreover, 54 hub genes rapidly responding to hypoosmotic stress were identified by WGCNA. The high-quality chromosome-scale S. ocellatus genome and candidate resistance-related gene sets will not only provide a genomic basis for genetic improvement via molecular breeding, but will also lay an important foundation for investigating the molecular regulation of rapid responses to stress.

Atlantic cod (Gadus morhua) MHC I localizes to endolysosomal compartments independently of cytosolic sorting signals

Major histocompatibility complex (MHC) class I and II are crucial for the adaptive immune system because they are involved in peptide presentation to T cells. Until recently, it was believed that MHC genes and their associated immune components had been conserved since their evolutionary emergence in jawed fish. However, sequencing of the Atlantic cod (Gadus morhua) genome revealed a loss of MHC class II genes, and an extreme expansion of MHC class I genes. These findings lead to the hypothesis that a loss of the MHC class II pathway coincided with a more versatile use of MHC class I, but so far there is no direct experimental evidence in support of this. To gain a deeper understanding of the function of the expanded MHC class I, we selected five MHC class I gene variants representing five of the six clades identified in previous studies and investigated their intracellular localization in human and Atlantic cod larval cells. Intriguingly, we uncovered that all selected MHC class I variants localize to endolysosomal compartments in Atlantic cod cells. Additionally, by introducing point mutations or deletions in the cytosolic tail, we found that hypothetical sorting signals in the MHC class I cytosolic tail do not influence MHC class I trafficking. Moreover, we demonstrated that in Atlantic cod, tapasin and MHC class I colocalize on endolysosomes suggesting that peptide-loading assistance and stabilization of MHC class I occurs outside the endoplasmic reticulum. Altogether, our results demonstrate that MHC class I from Atlantic cod is sorted to the endolysosomal system, which may indicate that it interacts with exogenous peptides for potential cross presentation.

Geographic variation in gene flow from a genetically distinct migratory ecotype drives population genetic structure of coastal Atlantic cod (Gadus morhua L.)

Identifying how physical and biotic factors shape genetic connectivity among populations in time and space is essential to our understanding of the evolutionary trajectory as well as the management of marine species. Atlantic cod is a widespread and commercially important marine species displaying several ecotypes with different life history strategies. Using three sets of SNPs: neutral, informative, and genome-inversion linked, we studied population genetic structure of ~2500 coastal Atlantic cod (CC) from 40 locations along Norway's 2500 km coastline, including nine fjords. We observed: (1) a genetic cline, suggesting a mechanism of isolation by distance, characterized by a declining F ST between CC and North East Arctic Cod (NEAC-genetically distinct migratory ecotype) with increasing latitude, (2) that in the north, samples of CC from outer-fjord areas were genetically more similar to NEAC than were samples of CC from their corresponding inner-fjord areas, (3) greater population genetic differentiation among CC sampled from outer-fjord areas along the coast, than among CC sampled from their corresponding inner-fjord areas, (4) genetic differentiation among samples of CC from both within and among fjords. Collectively, these results permit us to draw two main conclusions. First, that differences in the relative presence of the genetically highly distinct, migratory ecotype NEAC, declining from north to south and from outer to inner fjord, plays the major role in driving population genetic structure of the Norwegian CC. Second, that there is limited connectivity between CC from different fjords. These results suggest that the current management units implemented for this species in Norway should be divided into smaller entities. Furthermore, the situation where introgression from one ecotype drives population genetic structure of another, as is the case here, may exist in other species and geographical regions, thus creating additional challenges for sustainable fisheries management.

Supergene origin and maintenance in Atlantic cod

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Atlantic cod carries four supergenes linked to migratory lifestyle and environmental adaptations. Using whole-genome sequencing, the authors show that the genome inversions that underlie the supergenes originated at different times and show gene flux between supergene haplotypes.

A highly contiguous nuclear genome assembly of the mandarinfish Synchiropus splendidus (Syngnathiformes: Callionymidae)

The fish order Syngnathiformes has been referred to as a collection of misfit fishes, comprising commercially important fish such as red mullets as well as the highly diverse seahorses, pipefishes, and seadragons-the well-known family Syngnathidae, with their unique adaptations including male pregnancy. Another ornate member of this order is the species mandarinfish. No less than two types of chromatophores have been discovered in the spectacularly colored mandarinfish: the cyanophore (producing blue color) and the dichromatic cyano-erythrophore (producing blue and red). The phylogenetic position of mandarinfish in Syngnathiformes, and their promise of additional genetic discoveries beyond the chromatophores, made mandarinfish an appealing target for whole-genome sequencing. We used linked sequences to create synthetic long reads, producing a highly contiguous genome assembly for the mandarinfish. The genome assembly comprises 483 Mbp (longest scaffold 29 Mbp), has an N50 of 12 Mbp, and an L50 of 14 scaffolds. The assembly completeness is also high, with 92.6% complete, 4.4% fragmented, and 2.9% missing out of 4584 BUSCO genes found in ray-finned fishes. Outside the family Syngnathidae, the mandarinfish represents one of the most contiguous syngnathiform genome assemblies to date. The mandarinfish genomic resource will likely serve as a high-quality outgroup to syngnathid fish, and furthermore for research on the genomic underpinnings of the evolution of novel pigmentation.

Evidence for selection and spatially distinct patterns found in a putative zona pellucida gene in Pacific cod, and implications for management

Genetic differentiation has been observed in marine species even when no obvious barriers to gene flow exist, and understanding such differentiation is essential for effective fisheries management. Highly differentiated outlier loci can provide information on how genetic variation might not only contribute to local adaptation but may also be affected by historical demographic events. A locus which aligned to a predicted zona pellucida sperm-binding protein 3 gene (ZP3) in Atlantic cod (Gadus morhua) was previously identified as the highest outlier based on F ST in a RADseq study of Pacific cod (Gadus macrocephalus) across the West Coast of North America. However, because of the limited length of the RAD sequence and restricted geographic area of sampling, no conclusion on the functional significance of the observed variation was possible. In other marine species, ZP3 is involved in reproductive isolation, local adaptation, and has neofunctionalized as an antifreeze gene, and so it may provide important insights in functional population structure of Pacific cod. Here, we sequenced a 544-bp region of ZP3 in 230 Pacific cod collected from throughout their geographic range. We observed striking patterns of spatial structuring of ZP3 haplotypes, with a sharp break near Kodiak, Alaska, USA where populations within ~200 km of each other are nearly fixed for different haplotypes, contrasting a pattern of isolation by distance at other genetic markers in this region (F ST = 0.003). Phylogenetic analysis of ZP3 haplotypes revealed that the more southern haplotypes appear to be ancestral, with the northern haplotype evolving more recently, potentially in response to a novel selective pressure as Pacific cod recolonized northern latitudes after glaciation. The sharp break in haplotype frequencies suggests strong selective pressures are operating on small spatial scales and illustrates that selection can create high divergence even in marine species with ample opportunities for gene flow.

A long reads-based de-novo assembly of the genome of the Arlee homozygous line reveals chromosomal rearrangements in rainbow trout

Currently, there is still a need to improve the contiguity of the rainbow trout reference genome and to use multiple genetic backgrounds that will represent the genetic diversity of this species. The Arlee doubled haploid line was originated from a domesticated hatchery strain that was originally collected from the northern California coast. The Canu pipeline was used to generate the Arlee line genome de-novo assembly from high coverage PacBio long-reads sequence data. The assembly was further improved with Bionano optical maps and Hi-C proximity ligation sequence data to generate 32 major scaffolds corresponding to the karyotype of the Arlee line (2 N = 64). It is composed of 938 scaffolds with N50 of 39.16 Mb and a total length of 2.33 Gb, of which ∼95% was in 32 chromosome sequences with only 438 gaps between contigs and scaffolds. In rainbow trout the haploid chromosome number can vary from 29 to 32. In the Arlee karyotype the haploid chromosome number is 32 because chromosomes Omy04, 14 and 25 are divided into six acrocentric chromosomes. Additional structural variations that were identified in the Arlee genome included the major inversions on chromosomes Omy05 and Omy20 and additional 15 smaller inversions that will require further validation. This is also the first rainbow trout genome assembly that includes a scaffold with the sex-determination gene (sdY) in the chromosome Y sequence. The utility of this genome assembly is shown through the improved annotation of the duplicated genome loci that harbor the IGH genes on chromosomes Omy12 and Omy13.

Chromosome-level genome assembly of burbot (Lota lota) provides insights into the evolutionary adaptations in freshwater

The burbot (Lota lota) is the only member of the order Gadiformes adapted solely to freshwater. This species has the widest longitudinal range among freshwater fish worldwide. Burbot serves as a good model for studies on adaptive genome evolution from marine to freshwater environments. However, a high-quality reference genome of burbot has not yet been released. Here, the first chromosome-level genome of burbot was constructed using PacBio long sequencing and Hi-C technology. A total of 95.24 Gb polished PacBio sequences were generated, and the preliminary genome assembly was 575.83 Mb in size with a contig N50 size of 2.15 Mb. The assembled sequences were anchored to 22 pseudochromosomes by using Hi-C data. The final assembled genome after Hi-C correction was 575.92 Mb, with a contig N50 of 2.01 Mb and a scaffold N50 of 22.10 Mb. A total of 22,067 protein-coding genes were predicted, 94.82% of which were functionally annotated. Phylogenetic analyses indicated that burbot diverged with the Atlantic cod approximately 43.8 million years ago. In addition, 377 putative genes that appear to be under positive selection in burbot were identified. These positively selected genes might be involved in the adaptation to the freshwater environment. These genome data provide an invaluable resource for the ecological and evolutionary study of the order Gadiformes.