Genome Structural Variation Landscape and Its Selection Signatures in the Fast-growing Strains of the Pacific Oyster, Crassostrea gigas

Genomic selection for resistance to one pathogenic strain of Vibrio splendidus in blue mussel Mytilus edulis

Introduction

The blue mussel is one of the major aquaculture species worldwide. In France, this species faces a significant threat from infectious disease outbreaks in both mussel farms and the natural environment over the past decade. Diseases caused by various pathogens, particularly Vibrio spp., have posed a significant challenge to the mussel industry. Genetic improvement of disease resistance can be an effective approach to overcoming this issue.

Methods

In this work, we tested genomic selection in the blue mussel (Mytilus edulis) to understand the genetic basis of resistance to one pathogenic strain of Vibrio splendidus (strain 14/053 2T1) and to predict the accuracy of selection using both pedigree and genomic information. Additionally, we performed a genome-wide association study (GWAS) to identify putative QTLs underlying disease resistance. We conducted an experimental infection involving 2,280 mussels sampled from 24 half-sib families containing each two full-sib families which were injected with V. splendidus. Dead and survivor mussels were all sampled, and among them, 348 dead and 348 surviving mussels were genotyped using a recently published multi-species medium-density 60K SNP array.

Results

From potentially 23.5K SNPs for M. edulis present on the array, we identified 3,406 high-quality SNPs, out of which 2,204 SNPs were successfully mapped onto the recently published reference genome. Heritability for resistance to V. splendidus was moderate ranging from 0.22 to 0.31 for a pedigree-based model and from 0.28 to 0.36 for a genomic-based model.

Discussion

GWAS revealed the polygenic architecture of the resistance trait in the blue mussel. The genomic selection models studied showed overall better performance than the pedigree-based model in terms of accuracy of breeding values prediction. This work provides insights into the genetic basis of resistance to V. splendidus and exemplifies the potential of genomic selection in family-based breeding programs in M. edulis.

Construction and analysis of the chromosome-level haplotype-resolved genomes of two Crassostrea oyster congeners: Crassostrea angulata and Crassostrea gigas

BACKGROUND

The Portuguese oyster Crassostrea angulata and the Pacific oyster C. gigas are two major Crassostrea species that are naturally distributed along the Northwest Pacific coast and possess great ecological and economic value. Here, we report the construction and comparative analysis of the chromosome-level haplotype-resolved genomes of the two oyster congeners.

FINDINGS

Based on a trio-binning strategy, the PacBio high-fidelity and Illumina Hi-C reads of the offspring of the hybrid cross C. angulata (♂) × C. gigas (♀) were partitioned and independently assembled to construct two chromosome-level fully phased genomes. The assembly size (contig N50 size, BUSCO completeness) of the two genomes were 582.4 M (12.8 M, 99.1%) and 606.4 M (5.46 M, 98.9%) for C. angulata and C. gigas, respectively, ranking at the top of mollusk genomes with high contiguity and integrity. The general features of the two genomes were highly similar, and 15,475 highly conserved ortholog gene pairs shared identical gene structures and similar genomic locations. Highly similar sequences can be primarily identified in the coding regions, whereas most noncoding regions and introns of genes in the same ortholog group contain substantial small genomic and/or structural variations. Based on population resequencing analysis, a total of 2,756 species-specific single-nucleotide polymorphisms and 1,088 genes possibly under selection were identified.

CONCLUSIONS

This is the first report of trio-binned fully phased chromosome-level genomes in marine invertebrates. The study provides fundamental resources for the research on mollusk genetics, comparative genomics, and molecular evolution.

Two parallel chromosome-level reference genomes to support restoration and aquaculture of European flat oyster Ostrea edulis

This volume of Evolutionary Applications sees the publication of two genomes for the European native flat oyster Ostrea edulis, a species of significant evolutionary, ecological and commercial value. Each is a highly contiguous chromosome-level assembly from individuals of different genetic backgrounds, which have been benchmarked against one another. This situation has resulted from the serendipitous discovery that two independent research groups were both deep into the process of building, annotating and investigating separately produced assemblies. Due to constraints with funder requirements and the need to recognize early career researchers for their work, alongside the technical challenge of integrating assemblies from two very different genomes, there was limited capacity to merge the sequences into one publication at the stage of discovery. This issue is likely to become very common over the next few years until the technologies for working with multiple genomes at once, for example, graph genomes, become commonplace in nonmodel species. Consequently, both of our teams have decided to collaborate rather than compete, recognizing the benefit to copublishing two separate genome resources for the research community, each with distinct scientific investigations, and working collaboratively to benchmark the assemblies.

Comparative Methylome Analysis Reveals Epigenetic Signatures Associated with Growth and Shell Color in the Pacific Oyster, Crassostrea gigas

Fast growth is one of the most important breeding goals for all economic species such as the Pacific oyster (Crassostrea gigas), an aquaculture mollusk with top global production. Although the genetic basis and molecular mechanisms of growth-related traits have been widely investigated in the oyster, the role of DNA methylation involved in growth regulation remains largely unexplored. In this study, we performed a comparative DNA methylome analysis of two selectively bred C. gigas strains with contrasted phenotypes in growth and shell color based on whole-genome bisulfite sequencing (WGBS). Genome-wide profiling of DNA methylation at the single-base resolution revealed that DNA methylations were widely spread across the genome with obvious hotspots, coinciding with the distribution of genes and repetitive elements. Higher methylation levels were observed within genic regions compared with intergenic and promoter regions. Comparative analysis of DNA methylation allowed the identification of 339,604 differentially methylated CpG sites (DMCs) clustering in 27,600 differentially methylated regions (DMRs). Functional annotation analysis identified 11,033 genes from DMRs which were enriched in biological processes including cytoskeleton system, cell cycle, signal transduction, and protein biosynthesis. Integrative analysis of methylome and transcriptome profiles revealed a positive correlation between gene expression and DNA methylation within gene-body regions. Protein-protein interaction (PPI) analysis of differentially expressed and methylated genes allowed for the identification of integrin beta-6 (homolog of human ITGB3) as a hub modulator of the PI3K/Akt signaling pathway that was involved in various growth-related processes. This work provided insights into epigenetic regulation of growth in oysters and will be valuable resources for studying DNA methylation in invertebrates.