Generation of a chromosome‐level genome assembly for Pacific halibut ( Hippoglossus stenolepis ) and characterization of its sex‐determining genomic region

Advancing genetic improvement in the omics era: status and priorities for United States aquaculture

BACKGROUND

The innovations of the "Omics Era" have ushered in significant advancements in genetic improvement of agriculturally important animal species through transforming genetics, genomics and breeding strategies. These advancements were often coordinated, in part, by support provided over 30 years through the 1993-2023 National Research Support Project 8 (NRSP8, National Animal Genome Research Program, NAGRP) and affiliate projects focused on enabling genomic discoveries in livestock, poultry, and aquaculture species. These significant and parallel advances demand strategic planning of future research priorities. This paper, as an output from the May 2023 Aquaculture Genomics, Genetics, and Breeding Workshop, provides an updated status of genomic resources for United States aquaculture species, highlighting major achievements and emerging priorities.

MAIN TEXT

Finfish and shellfish genome and omics resources enhance our understanding of genetic architecture and heritability of performance and production traits. The 2023 Workshop identified present aims for aquaculture genomics/omics research to build on this progress: (1) advancing reference genome assembly quality; (2) integrating multi-omics data to enhance analysis of production and performance traits; (3) developing resources for the collection and integration of phenomics data; (4) creating pathways for applying and integrating genomics information across animal industries; and (5) providing training, extension, and outreach to support the application of genome to phenome. Research focuses should emphasize phenomics data collection, artificial intelligence, identifying causative relationships between genotypes and phenotypes, establishing pathways to apply genomic information and tools across aquaculture industries, and an expansion of training programs for the next-generation workforce to facilitate integration of genomic sciences into aquaculture operations to enhance productivity, competitiveness, and sustainability.

CONCLUSION

This collective vision of applying genomics to aquaculture breeding with focus on the highlighted priorities is intended to facilitate the continued advancement of the United States aquaculture genomics, genetics and breeding research community and industries. Critical challenges ahead include the practical application of genomic tools and analytical frameworks beyond academic and research communities that require collaborative partnerships between academia, government, and industry. The scope of this review encompasses the use of omics tools and applications in the study of aquatic animals cultivated for human consumption in aquaculture settings throughout their life-cycle.

Sex Influences the Genetic Structure of Greenland Halibut in the North Atlantic

Greenland halibut (Reinhardtius hippoglossoides) is a commercially important species in the North Atlantic whose spatial population structure has not yet been fully determined across its entire range. We genotyped individuals from across the North Atlantic using a subset of informative single nucleotide polymorphic (SNP) markers to assess their usability as a SNP panel. We assessed whether these purportedly structured SNPs had any association with sex. We found several of these loci to be in sex-determining chromosomes and that their inclusion generated genetic structure mainly in males. The population structure without the sex-associated SNPs was weak and followed an isolation-by-distance pattern, likely with a large regional population on each side of the North Atlantic. We discuss how different sex ratios in the samples and/or an evolving sex-determination system in this species likely caused the inclusion of sex-associated loci in the panel. We found suggestive evidence of polymorphisms at sex-determining chromosomes differentiating males on east and west locations, indicating evolution of the sex-determination system. These results highlight the importance of documenting sex-based differences in genetic studies and call for a better understanding of genomic architecture to understand sex-determination systems across the whole distribution of sexually dimorphic species.

Fully phased genome assemblies and graph-based genetic variants of the olive flounder, Paralichthys olivaceus

The olive flounder, Paralichthys olivaceus, also known as the Korean halibut, is an economically important flatfish in East Asian countries. Here, we provided four fully phased genome assemblies of two different olive flounder individuals using high-fidelity long-read sequencing and their parental short-read sequencing data. We obtained 42-44 Gb of ~15-kb and ~Q30 high-fidelity long reads, and their assembly quality values were ~53. We annotated ~30 K genes, ~170-Mb repetitive sequences, and ~3 M 5-methylcytosine positions for each genome assembly, and established a graph-based draft pan-genome of the olive flounder. We identified 5 M single-nucleotide variants and 100 K structural variants with their genotype information, where ~13% of the variants were possibly fixed in the two Korean individuals. Based on our chromosome-level genome assembly, we also explored chromosome evolution in the Pleuronectiformes family, as reported earlier. Our high-quality genomic resources will contribute to future genomic selection for accelerating the breeding process of the olive flounder.

Population genomic analysis reveals a polygenic sex determination system in Apostichopus japonicus

The sea cucumber Apostichopus japonicus, a key species in Chinese aquaculture, plays a significant evolutionary role within the Echinodermata phylum. However, the sex determination mechanism in this species remains poorly understood. Here, we conducted extensive sex surveys and sampling of eight wild populations, investigating the sex-related SNPs and insertion or deletions (indels) through bulk segregation analysis (BSA) and genome-wide association study (GWAS) analysis. Our findings suggest that A. japonicus employs a polygenic sex determination (PSD) system, with solute carrier family 8 (SLC8A) being the candidate gene for sex determination, encoding sodium-calcium exchanger (NCX1). The analysis of normalized sequencing depth reveals variations across chromosomes 6, 13, 14, 16, and 18, supporting the PSD system. We also identified 541.656 kb of male-specific sequences and screened five markers (C77185, C98086, C64977, C125, and C876) for molecular sex identification. Overall, this study provides new insights into A. japonicus sex determination, highlighting a complex multi-gene mechanism rather than a simple XX/XY system.

Diverse Origins of Near-Identical Antifreeze Proteins in Unrelated Fish Lineages Provide Insights Into Evolutionary Mechanisms of New Gene Birth and Protein Sequence Convergence

Determining the origins of novel genes and the mechanisms driving the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. Recently evolved fish antifreeze proteins (AFPs) offer a unique opportunity to explore these processes, particularly the near-identical type I AFP (AFPI) found in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages. Through comprehensive comparative analyses of newly sequenced genomes of winter flounder and grubby sculpin, along with available high-quality genomes of cunner and 14 other related species, the study revealed that near-identical AFPI proteins originated from distinct genetic precursors in each lineage. Each lineage independently evolved a de novo coding region for the novel ice-binding protein while repurposing fragments from their respective ancestors into potential regulatory regions, representing partial de novo origination-a process that bridges de novo gene formation and the neofunctionalization of duplicated genes. The study supports existing models of new gene origination and introduces new ones: the innovation-amplification-divergence model, where novel changes precede gene duplication; the newly proposed duplication-degeneration-divergence model, which describes new functions arising from degenerated pseudogenes; and the duplication-degeneration-divergence gene fission model, where each new sibling gene differentially degenerates and renovates distinct functional domains from their parental gene. These findings highlight the diverse evolutionary pathways through which a novel functional gene with convergent sequences at the protein level can evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

Transforming growth factor-β (TGF-β): A master signal pathway in teleost sex determination

Sex determination and differentiation in fish has always been a hot topic in genetic breeding of aquatic animals. With the advances in next-generation sequencing (NGS) in recent years, sex chromosomes and sex determining genes can be efficiently identified in teleosts. To date, master sex determination genes have been elucidated in 114 species, of which 72 species have sex determination genes belonging to TGF-β superfamily. TGF-β is the only signaling pathway that the largest proportion of components, which including ligands (amhy, gsdfy, gdf6), receptors (amhr, bmpr), and regulator (id2bby), have opportunity recognized as a sex determination gene. In this review, we focus on the recent studies about teleost sex-determination genes within TGF-β superfamily and propose several hypotheses on how these genes regulate sex determination process. Differing from other reviews, our review specifically devotes significant attention to all members of the TGF-β signal pathway, not solely the sex determination genes within the TGF-β superfamily. However, the functions of the paralogous genes of TGF superfamily are still needed ongoing research. Further studies are required to more accurately interpret the molecular mechanism of TGF-β superfamily sex determination genes.

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The convergently evolved fish antifreeze proteins provide excellent opportunities to investigate evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner and 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage. Each independently evolved a coding region for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms are consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Endocrine patterns associated with ovarian development in female Pacific halibut (Hippoglossus stenolepis)

The Pacific halibut (Hippoglossus stenolepis) is a large migratory demersal flatfish species that occupies a top trophic role in the North Pacific Ocean and Bering Sea ecosystems, where it also supports various fisheries. As a first attempt to characterize the endocrine mechanisms driving sexual maturation in this important species, we collected pituitary, ovarian and blood samples from Pacific halibut females captured in the wild that were classified histologically into various female developmental stages. We conducted gene expression analyses of gonadotropin beta subunits in the pituitary and observed that mRNA expression levels of fshb gradually increased throughout vitellogenesis, remained elevated until before ovulation and declined after spawning. In contrast, the mRNA expression levels of lhb markedly increased during oocyte maturation and remained elevated until after spawning. Ovarian mRNA expression levels of the gonadotropin receptor genes fshr and lhr peaked during oocyte maturation and before spawning, respectively, immediately following the developmental stage at which pituitary fshb and lhb mRNA expression first reached maximum levels. The ovarian gene expression patterns of steroidogenic enzyme genes cyp19a1 and hsd20b2 paralleled those of fshr and lhr, respectively. Testosterone and 17β-estradiol (E2) plasma levels increased concomitantly with fshr and cyp19a1 mRNA expression levels, and vitellogenin plasma levels increased throughout vitellogenesis and reached maximum levels prior to spawning. These results are consistent with the notion that in female Pacific halibut, as in other teleosts, vitellogenesis and oocyte maturation and ovulation are likely under the control of pituitary gonadotropic hormones Fsh and Lh, respectively.

Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

SNP Panel and Genomic Sex Identification in Atlantic Halibut (Hippoglossus hippoglossus)

The ability to identify sex is necessary in population biology for a proper understanding of the dynamics of a population. In Atlantic halibut, phenotypic sex identification is not possible due to the lack of significant external morphological differences. We developed an Illumina SNP panel for Atlantic halibut with 4000 SNPs spread evenly throughout the genome with a minor allele frequency MAF ≥ 0.4, except for N = 249 SNPs located in a sex-determining region on chromosome 12, N = 176 of these SNPs were selected to genetically identify male and female individuals using a DAPC analysis. The genomic identification of sex allows for non-lethal sex determination and validation of sex identification in the field. The SNP panel is a new genomic resource for Atlantic halibut that will make it possible to generate the genotypic data for the large number of individuals needed to estimate population abundance using genomics and the Close Kin Mark Recapture (CKMR) approach, an emerging component of fisheries management and stock monitoring.