Construction of a chromosome-level genome and variation map for the Pacific oyster Crassostrea gigas

Transcriptome profiling of triploid Crassostrea gigas gills indicates the host immune mechanism against bacterial infection

As an important member of global aquaculture, oysters (Crassostrea gigas) have significant economic value. With the development of commercial aquaculture, the frequent occurrence of diseases caused by Vibrio alginolyticus has become a hindrance to high-density aquaculture. Gill tissue, as an important component of immune system of the oysters, plays the key point in the face of invasion by foreign substances. Compared to the diploid oyster, the triploid oyster presents a higher growth rate and lower growth investment, making it a more ideal model for studying oyster immune defense. In this study, triploid oysters were as the research subject, and gill tissues attacked by V. alginolyticus were sequenced. By analyzing samples from different time points, 1746 DEGs were obtained. The KEGG and GO functional enrichment analysis showed that gill tissues mainly participate in immune function through the PIK3-Akt signaling pathway and the MAPK signaling pathway. The protein interaction network revealed three genes (CASP8, CASP9 and PIK3CA) that play core roles in immune defense by analyzing the interaction relationship between genes. Finally, qRT-PCR verified the expression of key genes. This study provides a more effective scientific basis for disease prevention and control of oysters and other bivalve shellfish, and helps to promote the sustainable development of aquaculture.

Comparative transcriptomic analysis reveals a differential acid response mechanism between estuarine oyster (Crassostrea ariakensis) and Pacific oyster (Crassostrea gigas)

Ocean and coastal acidification (OCA) poses a significant and rapidly emerging threat to mollusks. The physiological resilience of mollusks to OCA varies considerably; however, the underlying molecular mechanisms remain poorly understood. Seawater in estuaries, being more susceptible to acidification than that in open coastal zones, may enhance the tolerance of resident mollusks to low pH levels. Here, we conducted a comparative analysis between estuarine oysters (Crassostrea ariakensis) and Pacific oysters (Crassostrea gigas) using physiological phenotype and transcriptomic analyses to reveal differential acid-tolerance mechanisms in response to constant pH of 7.8. Our findings indicated that survival and respiration rates of C. ariakensis, which inhabits estuaries with fluctuating pH levels, were higher than those of C. gigas, which inhabits open coastal zones with relative stable pH conditions. Acid-responsive genes identified in C. gigas, including molecular chaperones and immune-related genes, exhibited higher constitutive expression in C. ariakensis under control conditions. Co-expression analyses revealed that C. ariakensis mitigated the effects of low pH by expressing genes involved in ion transporter activity and translation control. C. gigas activated genes associated with glycolipid metabolism while inhibiting cell division during acid stress. These findings suggested that C. ariakensis has evolved into a more energy-efficient regulatory network than C. gigas, incorporating both front-loading and responsive mechanisms to maintain acidbase homeostasis. This study is the first to investigate acid-tolerance differences between mollusks inhabiting estuarine and open coastal environments and provides critical insights into the resilience of mollusks in increasingly acidified oceans.

Genome of Kumamoto Oyster Crassostrea sikamea Provides Insights Into Bivalve Evolution and Environmental Adaptation

The Kumamoto oyster, Crassostrea sikamea, is a marine bivalve naturally distributed along the coasts of southern China and southern Japan, with a hatchery population that has been under domestication in the United States since its introduction from Japan in the 1940s. To understand its evolutionary history and environmental adaptation, we produced a chromosome-level genome assembly of C. sikamea and conducted whole-genome resequencing of 141 individuals from the US hatchery population and six wild populations from China and Japan. The assembled genome of C. sikamea has a size of 616 Mb covering all 10 chromosomes with a contig N50 of 4.21 Mb and a scaffold N50 of 62.25 Mb. Phylogenetic analysis indicated that C. sikamea diverged from the Crassostrea angulata and Crassostrea gigas clade about 9.9 million years ago. Synteny analysis revealed significant chromosomal rearrangements during bivalve evolution leading to oysters, but remarkable conservation of all 10 oyster chromosomes over ~180 million years, a surprising disparity in chromosomal evolution. Phylogenetic analysis produced three distinct clusters for the US, Japanese, and Chinese populations, with the US population closer to the Japanese population, confirming its origin. No differentiation was detected among the five Chinese populations, indicating strong gene flow. Between the US and Japan populations, 402 genes exhibited selection signals, including three myosin heavy chain genes that were also differentiated in domesticated lines of the eastern oyster, suggesting changes in these genes may be important for domestic production. Among the 768 genes showing selection signals between natural populations of Japan and China, genes related to stress response are most enriched, suggesting responding to environmental stress is critical for local adaptation. These findings provide insights into bivalve evolution and environmental adaptation, as well as useful resources for comparative genomics and genetic improvement of cultured Kumamoto oyster stocks.

Optimization of Microsatellite Multiplex PCRs for Triploidy Verification and Genetic Diversity Assessment in the Pacific Oyster, Crassostrea gigas

The ploidy detection is crucial for the oyster industry. The objective of this study was to develop a method that verifies ploidy of the triploid Pacific oyster Crassostrea gigas by analyzing the diversity of triploid through microsatellite multiplex PCRs using fluorescent universal primers. We developed four information-rich multiplex PCR panels, comprising a total of 12 genomic microsatellites located in the genome of the C. gigas, distributed across seven chromosomes with an average of 14 alleles per locus. Each panel used M13(-21) primers labeled with specific fluorochrome dyes, and the forward primers for each locus were appended with M13(-21) sequences. We validated the approach to infer ploidy using flow cytometry as a reference, finding > 95% agreement between these methods, and demonstrated its potential utility to infer aneuploidy. Genotyping of 496 triploid samples from eight populations yielded 10 or more alleles per locus in 99.63% of samples in a single capillary electrophoresis. The correct assignment of triploidy depends on the number of markers with three unique allele fragments (MNM). Using semi-strict criteria of three unique alleles at one or more loci, the detection accuracy rate was 95.26% for triploids. Using the strict criteria of three unique alleles at two or more loci, the detection accuracy rate was 98.34%. Populations with reduced genetic diversity due to selective breeding were better suited for the semistrict criterion, maximizing triploid detection. And cultured populations were more suitable for evaluation using the strict criteria, which effectively reduced false-positive diploid assignment and increased triploid detection accuracy. The markers developed in this study were highly polymorphic and effective for assessing genetic diversity and distinguishing populations, providing a reliable tool for triploid detection and analysis in oyster breeding.

The genome sequence of the Pacific oyster, Magallana gigas (Thunberg, 1793)

We present a genome assembly from an individual Magallana gigas (the Pacific oyster; Mollusca; Bivalvia; Ostreida; Ostreidae). The genome sequence is 564.0 megabases in span. Most of the assembly is scaffolded into 10 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 18.23 kilobases in length. Gene annotation of this assembly on Ensembl identified 19,775 protein-coding genes.

Genetic Variants Affecting FADS2 Enzyme Dynamics and Gene Expression in Cogenetic Oysters with Different PUFA Levels Provide New Tools to Improve Unsaturated Fatty Acids

Long-chain polyunsaturated fatty acids (LC-PUFAs) are crucial for human health and cannot be produced internally. Bivalves, such as oysters, serve as valuable sources of high-quality PUFAs. The enzyme fatty acid desaturase (FADS) plays a key role in the metabolism of LC-PUFAs. In this study, we conducted a thorough genome-wide analysis of the genes belong to the FADS family in Crassostrea gigas and Crassostrea angulata, with the objective of elucidating the function of the FADS2 and investigating the genetic variations that affect PUFA biosynthesis. We identified six FADS genes distributed across four chromosomes, categorized into three subfamilies. The coding region of FADS2 revealed five non-synonymous mutations that were shown to influence protein structure and stability through molecular dynamics simulations. The promoter region of FADS2 contains ten SNPs and three indels significantly correlated with PUFA content. These genetic variations may explain the differences in PUFA levels observed between the two oyster species and could have potential applications in enhancing PUFA content. This study improves the molecular understanding of PUFA metabolism in oysters and presents a potential strategy for selecting oysters with high PUFA levels.

Genome-wide analysis of the cytochrome P450 gene family in Pacific oyster Crassostrea gigas and their expression profiles during gonad development

The cytochrome P450 (CYP) gene superfamily plays a significant role in various physiological processes, producing different compounds such as hormones, fatty acids, and biomolecules. However, little information is known their roles during gonad development in Pacific oyster (Crassostrea gigas). In this study, total of 116 CgCYP (Crassostrea gigas cytochrome P450) genes were identified and their expression pattern was analyzed for the first time. The relative molecular weights of these CgCYP genes ranged from 63.52 to 113.41 kDa, and the length of encoded amino acids ranged from 103 to 993. And total 26 cis-acting elements of these CgCYP genes were identified. GO and KEGG enrichment analysis showed some CgCYP genes are essential for the metabolism of male and female sex hormones. Additionally, expression anslysis showed 69 CgCYP genes were over-expressed in early gonad development and triploid infertile individuals. More importantly, expression levels of CgCYP1, CgCYP15, CgCYP34, CgCYP46, CgCYP69, CgCYP87, CgCYP88, and CgCYP103, were found to be significantly higher in female gonad, suggesting their important roles in female gonad development. The results of this study will provide a better understanding of the CgCYP genes in the gonad development of Pacific oyster.

Chromosome-level genome assembly of the Suminoe oyster Crassostrea ariakensis in south China

The Suminoe oyster Crassostrea ariakensis (Fujita, 1913) is one of the most important ecological and fishery bivalve mollusks with a worldwide distribution. Here, we reported an improved high-quality chromosomal-level genome assembly of C. ariakensis inhabiting the South China Sea, using Nanopore technology, Illumina sequencing, and high-throughput chromosomal conformation capture analysis. The assembled genome size is 631.73 Mb, with contig N50 length of 5.36 Mb and scaffold N50 length of 61.15 Mb, and is assigned to 10 chromosomes. A total of 29,357 protein-coding genes are predicted, 96.68% of which are functionally annotated. The genome contains 347.11 Mb (54.94%) of repetitive elements and 1130 noncoding RNAs. This improved genome assembly of south C. ariakensis is an important resource for understanding oyster diversity and evolution, and provides insights into genetic improvement, protection and management of oyster resource.

SNP Fingerprinting for Germplasm Identification of the Fast-Growing Pacific oyster (Crassostrea gigas) "Haida No. 1" Variety

The Pacific oyster (Crassostrea gigas) is a global aquaculture species of economic significance. Selective breeding programs have been conducted to produce multiple strains with fast growth as well as other desirable traits. However, due to the phenotypic plasticity of oysters, challenges existed for precise germplasm identification among selectively bred strains. In this work, we identified selection signatures of three fast-growing Pacific oyster strains originated from wild populations collected from China, Japan, and Korea, respectively, which were used for development of SNP-based molecular fingerprinting for precise identification of germplasm. We performed whole-genome resequencing of 59 oysters from three selectively bred strains and a wild population for genome-wide SNP analyses. Population structure analysis with these SNPs revealed significant genetic differentiation among the selectively bred strains. Based on the FST index, we identified 41, 49, and 36 strain-specific SNPs from the three selectively bred strains. Taking into account the "hitch-hiking effect" that occurs in the genome during positive selection, we identified two, three, and two molecular fingerprints for the three strains, respectively. We validated the molecular fingerprints of the China selectively bred strain (i.e., "Haida No. 1" variety) with a separate population of 42 oysters with diverse genetic background, demonstrating the accuracy of germplasm identification of over 96%. This work provides a reliable tool for precise germplasm identification of the "Haida No. 1" variety as well as other two selectively bred strains, which is valuable in germplasm conservation and breeding design in the C. gigas.

Novel Ser74 of NF-κB/IκBα phosphorylated by MAPK/ERK regulates temperature adaptation in oysters

Phosphorylation of Ser32 and Ser36 controls the degradation of IκBα is the conserved cascade mechanisms of immune core signaling pathway, NF-κB pathway in metazoans, but it's response to abiotic stress and the presence of novel phosphorylation mechanisms in other species remain unclear. Herein, we reported a novel heat-induced phosphorylation site (Ser74) at oysters' major IκBα, which independently regulated ubiquitination-proteasome degradation without the requirement of phosphorylation at S32 and S36. And this site was phosphorylated by ERK/MAPK pathway, which then promoted REL nuclear translocation to activate cell survival related genes to defend heat-stress. The MAPK-NF-κB cascade exhibited divergent thermal responses and adaptation patterns between two congeneric oyster species with differential habitat temperatures, indicating its involvement in shaping temperature adaptation. This study demonstrated that the existence of complex and unique phosphorylation-mediated signaling transduction mechanism in marine invertebrates, and expanded our understanding of the evolution and function of established classical pathway crosstalk mechanisms.

PI3K-AKT-mediated phosphorylation of Thr260 in CgCaspase-3/6/7 regulates heat-induced activation in oysters

Cysteine-aspartic proteases (caspases) are critical drivers of apoptosis, exhibiting expansion and domain shuffling in mollusks. However, the functions and regulatory mechanisms of these caspases remain unclear. In this study, we identified a group of Caspase-3/6/7 in Bivalvia and Gastropoda with a long inter-subunit linker (IL) that inhibits cleavage activation. Within this region, we found that conserved phosphorylation at Thr260 in oysters, mediated by the PI3K-AKT pathway, suppresses heat-induced activation. This mechanism is involved in divergent temperature adaptation between two allopatric congeneric oyster species, the relatively cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. Our study elucidates the role of these effector caspase members and their long IL in bivalves, revealing that the PI3K-AKT pathway phosphorylates Thr260 on CgCASP3/6/7's linker to inhibit heat-induced activation. These findings provide insights into the evolution and function of apoptotic regulatory mechanisms in bivalves.

An amplicon panel for high-throughput and low-cost genotyping of Pacific oyster

Maintaining genetic diversity in cultured shellfish can be challenging due to high variance in individual reproductive success, founder effects, and rapid genetic drift, but is important to retain adaptive potential and avoid inbreeding depression. To support broodstock management and selective breeding in cultured Pacific oysters (Crassostrea (Magallana) gigas), we developed an amplicon panel targeting 592 genomic regions and SNP variants with an average of 50 amplicons per chromosome. Target SNPs were selected based on elevated observed heterozygosity or differentiation in Pacific oyster populations in British Columbia, Canada. The use of the panel for parentage applications was evaluated using multiple generations of oysters from a breeding program on Vancouver Island, Canada (n = 181) and families selected for Ostreid herpesvirus-1 resistance from the Molluscan Broodstock Program in Oregon, USA (n = 136). Population characterization was evaluated using wild, naturalized, farmed, or hatchery oysters sampled throughout the Northern Hemisphere (n = 189). Technical replicates showed high genotype concordance (97.5%; n = 68 replicates). Parentage analysis found suspected pedigree and sample handling errors, demonstrating the panel's value for quality control in breeding programs. Suspected null alleles were identified and found to be largely population dependent, suggesting population-specific variation impacting target amplification. Null alleles were identified using existing data without the need for pedigree information, and once they were removed, assignment rates increased to 93.0 and 86.0% of possible assignments in the two breeding program datasets. A pipeline for analyzing the amplicon sequence data from sequencer output, amplitools, is also provided.

MAPK/ERK-PK(Ser11) pathway regulates divergent thermal metabolism of two congeneric oyster species

Pyruvate kinase (PK), as a key rate-limiting enzyme in glycolysis, has been widely used to assess the stress tolerance and sensitivity of organisms. However, its phosphorylation regulatory mechanisms mainly focused on human cancer research, with no reports in marine organisms. In this study, we firstly reported a conserved PK Ser11 phosphorylation site in mollusks, which enhanced enzyme activity by promoting substrate binding, thereby regulating divergent thermal metabolism of two allopatric congeneric oyster species with differential habitat temperature. It was phosphorylated by ERK kinase, and regulated by the classical MAPK pathway. The MAPK/ERK-PK signaling cascade responded to increased environmental temperature and exhibited stronger activation pattern in the relatively thermotolerant species (Crassostrea angulata), indicating its involvement in shaping temperature adaptation. These findings highlight the presence of complex and unique phosphorylation-mediated signaling transduction mechanisms in marine organisms, and provide new insights into the evolution and function of the crosstalk between classical pathways.

Exploration of Molecular Mechanisms of Immunity in the Pacific Oyster (Crassostrea gigas) in Response to Vibrio alginolyticus Invasion

Over the years, oysters have faced recurring mass mortality issues during the summer breeding season, with Vibrio infection emerging as a significant contributing factor. Tubules of gill filaments were confirmed to be in the hematopoietic position in Crassostrea gigas, which produce hemocytes with immune defense capabilities. Additionally, the epithelial cells of oyster gills produce immune effectors to defend against pathogens. In light of this, we performed a transcriptome analysis of gill tissues obtained from C. gigas infected with Vibrio alginolyticus for 12 h and 48 h. Through this analysis, we identified 1024 differentially expressed genes (DEGs) at 12 h post-injection and 1079 DEGs at 48 h post-injection. Enrichment analysis of these DEGs revealed a significant association with immune-related Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. To further investigate the immune response, we constructed a protein-protein interaction (PPI) network using the DEGs enriched in immune-associated KEGG pathways. This network provided insights into the interactions and relationships among these genes, shedding light on the underlying mechanisms of the innate immune defense mechanism in oyster gills. To ensure the accuracy of our findings, we validated 16 key genes using quantitative RT-PCR. Overall, this study represents the first exploration of the innate immune defense mechanism in oyster gills using a PPI network approach. The findings provide valuable insights for future research on oyster pathogen control and the development of oysters with enhanced antimicrobial resistance.

The role of the balance between energy production and ammonia detoxification mediated by key amino acids in divergent hypersaline adaptation among crassostrea oysters

Global ocean salinity is changing under rapid climate change and intensified anthropogenic activity. Increased differences in salinity threaten marine biodiversity, organismal survival, and evolution, particularly sessile invertebrates dwelling in highly fluctuating intertidal and estuarine environments. Comparing the responses of closely related species to salinity changes can provide insights into the adaptive mechanisms underlying inter- and intraspecific divergence in salinity tolerance, but are poorly understood in marine bivalves. We collected wild individuals of four Crassostrea species, in addition to two populations of the same species from their native habitats and determined the dynamics of hydrolyzed amino acids (HAAs) and transcriptional responses to hypersaline stress. In response to hypersaline stress, species/populations inhabiting natural high-salinity sea environments showed higher survival and less decline in HAAs than that of congeners inhabiting low-salinity estuaries. Thus, native environmental salinity shapes oyster tolerance. Notably, a strong negative correlation between the decline in HAAs and survival indicated that the HAAs pool could predict tolerance to hypersaline challenge. Four HAAs, including glutamine (Glu), aspartic acid (Asp), alanine (Ala) and glycine (Gly), were identified as key amino acids that contributed substantially to the emergency response to hypersaline stress. High-salinity-adapted oyster species only induced substantial decreases in Glu and Asp, whereas low-salinity-adapted congeners further incresaed Ala and Gly metabolism under hypersaline stress. The dynamics of the content and gene expression responsible for key amino acids pathways revealed the importance of maintaining the balance between energy production and ammonia detoxification in divergent hypersaline responses among oyster species/populations. High constructive or plastic expression of evolutionarily expanded gene copies in high-salinity-adapted species may contribute to their greater hypersaline tolerance. Our findings reveal the adaptive mechanism of key amino acids in salinity adaptation in marine bivalves and provide new avenues for the prediction of adaptive potential and aquaculture with high-salinity tolerant germplasms.

Whole genome sequencing of Crassostrea ariakensis (Mollusca: Ostreidae) and C. hongkongensis expands understandings of stress resistance in sessile oysters

To understand the environmental adaptations among sessile bivalves lacking adaptive immunity, a series of analyses were conducted, with special emphasis on the widely distributed C. ariakensis. Employing Pacbio sequencing and Hi-C technologies, whole genome for each of a C. ariakensis (southern China) and C. hongkongensis individual was generated, with the contig N50 reaching 6.2 and 13.0 Mb, respectively. Each genome harbored over 30,000 protein-coding genes, with approximately half of each genome consisting of repeats. Genome alignment suggested possible introgression between C. gigas and C. ariakensis (northern China), and re-sequencing data corroborated this result and indicated significant gene flow between C. gigas and C. ariakensis. These introgressed candidates, well-represented by genes related to immunity and osmotic pressure, may be associated with environmental stresses. Gene family dynamics modeling suggested immune-related genes were well represented among the expanded genes in C. ariakensis. These outcomes could be attributed to the spread of C. ariakensis.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Genome-wide identification and expression analysis of Dmrt gene family and their role in gonad development of Pacific oyster (Crassostrea gigas)

Doublesex and Mab-3-related transcription factor (Dmrt) is a type of transcription factor with a zinc-finger DM structural domain, which plays a significant role in sex determination and differentiation in animals. Although Dmrt has been reported in many vertebrates and invertebrates, it has rarely been studied in bivalves. In this study, a total of three members of the Dmrt gene family were identified and characterized in Crassostrea gigas, and all these CgDmrt genes contained a conserved DM domain. Analysis of the phylogenetic tree and gene structure revealed that Dmrt genes clustered on one branch may have similar functions in bivalves. Expression profiling of CgDmrt mRNA in different tissues and stages of gonad development indicated that CgDmrt3 exhibited sexually dimorphic expression and played an important role in the development of the male gonad in C. gigas. Furthermore, analysis of CgDmrt mRNA expression between fertile triploids and sterile triploids showed that CgDmrt3 may be involved in sperm production. Collectively, the systematic analysis of the CgDmrt genes will provide potential insights into the function of these genes in gonadal development.

Genome-Wide Association Analysis of Heat Tolerance in F2 Progeny from the Hybridization between Two Congeneric Oyster Species

As the world's largest farmed marine animal, oysters have enormous economic and ecological value. However, mass summer mortality caused by high temperature poses a significant threat to the oyster industry. To investigate the molecular mechanisms underlying heat adaptation and improve the heat tolerance ability in the oyster, we conducted genome-wide association analysis (GWAS) analysis on the F2 generation derived from the hybridization of relatively heat-tolerant Crassostrea angulata ♀ and heat-sensitive Crassostrea gigas ♂, which are the dominant cultured species in southern and northern China, respectively. Acute heat stress experiment (semi-lethal temperature 42 °C) demonstrated that the F2 population showed differentiation in heat tolerance, leading to extremely differentiated individuals (approximately 20% of individuals die within the first four days with 10% survival after 14 days). Genome resequencing and GWAS of the two divergent groups had identified 18 significant SNPs associated with heat tolerance, with 26 candidate genes located near these SNPs. Eleven candidate genes that may associate with the thermal resistance were identified, which were classified into five categories: temperature sensor (Trpm2), transcriptional factor (Gata3), protein ubiquitination (Ube2h, Usp50, Uchl3), heat shock subfamily (Dnajc17, Dnaja1), and transporters (Slc16a9, Slc16a14, Slc16a9, Slc16a2). The expressional differentiation of the above genes between C. gigas and C. angulata under sublethal temperature (37 °C) further supports their crucial role in coping with high temperature. Our results will contribute to understanding the molecular mechanisms underlying heat tolerance, and provide genetic markers for heat-resistance breeding in the oyster industry.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

Transcriptome profiling explores the immune defence mechanism of triploid Pacific oyster (Crassostrea gigas) blood against Vibrio alginolyticus based on protein interaction networks

Triploid oysters have provided the oyster industry with many benefits, such as fast growth rates, meat quality improvement, and increased oyster production and economic benefits, since the first report on triploid oysters was published. The development of polyploid technology has remarkably increased the output of triploid oysters to meet the increasing demand of consumers for Crassostrea gigas in the past decades. At present, research on triploid oyster has mainly focused on breeding and growth, but studies on the immunity of triploid oysters are limited. According to recent reports, Vibrio alginolyticus is a highly virulent strain that can cause disease and death in shellfish, shrimp, as well as serious economic losses. V. alginolyticus may be a reason why oysters die during summer. Therefore, using V. alginolyticus to explore the resistance and immune defense mechanisms of triploid oysters against pathogens presents practical significance. Transcriptome analysis of gene expression was performed in triploid C. gigas at 12 and 48 h after infection with V. alginolyticus, and the respective 2257 and 191 differentially expressed genes (DEGs) were identified. The results of GO and KEGG enrichment analyses showed that multiple significantly enriched GO terms and KEGG signaling pathways are associated with immunity. A protein-protein interaction network was constructed to investigate the interaction relationship of immune-related genes. Finally, we verified the expression situation of 16 key genes using quantitative RT-PCR. This study is the first to use the PPI network in exploring the immune defense mechanism of triploid C. gigas blood to fill the gap in the immune mechanism of triploid oysters and other mollusks, and provide valuable reference for future triploid farming and pathogen prevention and control.

Species-specific CgCaspase-Cg-5 in the pacific oyster induces haemocyte apoptosis by regulating the mRNA expression of apoptosis-related genes in the early stage of immune response

Caspases are cysteinyl aspartate-specific proteinases, playing critical roles in apoptotic pathway to induce apoptosis and inflammatory response. In this study, the expanded repertoire of Caspases was revealed in the Pacific oyster Crassostrea gigas, and a total of 30 Caspases were identified from the genomic and stress-induced transcriptomic databases of the Pacific oyster. They were clustered into CgCaspase-2/9, CgCaspase-8/10, CgCaspase-3/6/7, CgCaspase-Cg, and CgCaspase-L. CgCaspase-Cg subgroup was found to be specifically expanded after a positive selection in oyster with average Ka/Ks of 0.50. The mRNA expression of CgCaspase-Cg-5 was found to be obviously induced against various bacterial and viral stimulations or environmental stresses. The relative expression level of CgCaspase-Cg-5 in haemocytes increased and reached the peak at 6 h after Vibrio splendidus stimulation, which was 5.57-fold of that in the control group (p < 0.01). In the oysters whose CgCaspase-Cg-5 expression was knocked down, the mRNA expression of apoptosis-related genes including CgBcl2, CgBax, CgCaspase3 and CgCaspase9 changed significantly at 12 h after V. splendidus stimulation. The expression of CgBax, CgCaspase3 and CgCaspase9 decreased, which was 0.64-fold (p < 0.05), 0.53-fold (p < 0.05) and 0.62-fold (p < 0.01), while the expression of CgBcl2 increased, which was 2.81-fold (p < 0.01) of that in the EGFP-dsRNA group, respectively. Meanwhile, the apoptotic rate of haemocytes (1.90 ± 0.71%) significantly decreased compared to that in the EGFP-dsRNA group (5.40 ± 0.72%) (p < 0.05), and the histological damages of widened cell spacing, gill filament swelling and loose cytoplasm were observed in the CgCaspase-Cg-5-knockdown oysters after V. splendidus stimulation. Collectively, CgCaspase-Cg subgroup was specifically expanded in oyster and some bivalve species, and species-specific CgCaspase-Cg-5 regulated the mRNA expression of the apoptosis-related genes to induce haemocyte apoptosis in the early stage of immune response. This provided insight into the evolutionary and functional characteristics of Caspase repertoire in the Pacific oyster and highlighted the important role of CgCaspase-Cg-5 in the response to pathogen infection and environmental stresses.

Transgenerational effects of intertidal environment on physiological phenotypes and DNA methylation in Pacific oysters

Climate change and intensifying human activity are posing serious threats to marine organisms. The fluctuating intertidal zone forms a miniature ecosystem of a rapidly changing environment for studying biological adaptation. Transgenerational plasticity (TGP), an evolutionary phenomenon in which parental experience influences offspring phenotypes, provides an avenue for adaptation, but the molecular mechanism was poorly understood in marine molluscs. In this study, wild Pacific oysters (Crassostrea gigas), which were collected from intertidal zones, were used to conduct two-generation breeding in a subtidal area combined with a heat shock experiment in the laboratory to investigate the intertidal environment-induced TGP under temperate subtidal condition and thermally exposed condition, respectively. We showed that TGP could influence the physiological phenotypes related to the status of oxidation and energy in non-stress-exposed subtidal offspring for at least two generations. Genomic DNA methylation exhibited heritable divergence between intertidal and subtidal oysters, and 1655 (or 42.83 %) differentially methylated genes (DMGs) in F₀ were continuously reserved to F₂, which may mediate physiological TGP by participating in biological processes including macromolecule metabolism, cellular responses to stress, and the positive regulation of molecular function, especially fatty acid metabolism. The intertidal experience also influenced the thermal plasticity of physiological phenotypes within and across generations. Totally, 320 (or 14.74 %) specific thermal response DMGs in the intertidal F₀ generation were identified in F₁ and F₂, participating in pathways including carbohydrate, lipid, and energy metabolism, signal transduction, and the organismal immune system, which suggested transgenerational intertidal effect mediated by these genes could positively contribute to stress adaptation and had potential applications for aquaculture. This study demonstrates an epigenetic mechanism for TGP in stress adaptation in marine molluscs, and provides new avenues to improve the stress adaptation for marine resource conservation and aquaculture.

Experimental DNA Demethylation Reduces Expression Plasticity and Thermal Tolerance in Pacific Oysters

Increasing seawater temperatures pose a great threat to marine organisms, especially those settled in fluctuating intertidal areas. DNA methylation, which can be induced by environmental variation, can influence gene expression and mediate phenotypic plasticity. However, the regulatory mechanisms of DNA methylation in gene expression-mediated adaptation to environmental stress have rarely been elucidated. In this study, DNA demethylation experiments were conducted on a typical intertidal species, the Pacific oyster (Crassostrea gigas), to determine the direct role of DNA methylation in regulating gene expression and adaptability under thermal stress. The global methylation level and the expression level of DNA methyltransferases (DNMT1, DNMT3a) showed an accordant variation trend under high temperatures, supporting that the genomic methylation status was catalyzed by DNMTs. DNA methylation inhibitor 5-Azacytidine (5-Aza) effectively inhibited DNA methylation level and decreased methylation plasticity at the 6th hour in thermal conditions. In total, 88 genes were identified as candidate DNA methylation-regulated thermal response genes; they exhibited reduced expression plasticity in response to heat stress, possibly caused by the decreased methylation plasticity. Post-heat shock, the thermal tolerance indicated by the survival curve was reduced when oysters were pretreated with 5-Aza, meaning that DNA demethylation negatively affected thermal adaptation in oysters. This study provides direct evidence for the crucial role of DNA methylation in mediating stress adaptation in marine invertebrates and contributes to the theoretical foundations underlying marine resource conservation and aquaculture.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

A chromosome-level genome assembly of Ostrea denselamellosa provides initial insights into its evolution

The oyster Ostrea denselamellosa is a live-bearing species with a sharp decline in the natural population. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in O. denselamellosa. Here, we carried out the first whole genome sequencing at the chromosome-level in O. denselamellosa. Our studies yielded a 636 Mb assembly with scaffold N50 around 71.80 Mb. 608.3 Mb (95.6% of the assembly) were anchored to 10 chromosomes. A total of 26,412 protein-coding genes were predicted, of which 22,636 (85.7%) were functionally annotated. By comparative genomics, we found that long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) made up a larger proportion in O. denselamellosa genome than in other oysters'. Moreover, gene family analysis showed some initial insight into its evolution. This high-quality genome of O. denselamellosa provides a valuable genomic resource for studies of evolution, adaption and conservation in oysters.

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.

A high-quality chromosome-level genome assembly of the bivalve mollusk Mactra veneriformis

Mactra veneriformis (Bivalvia: Mactridae) is a bivalve mollusk of major economic importance in China. Decreased natural yields of M. veneriformis have led to an urgent need for genomic resources. To address this problem and the currently limited knowledge of molecular evolution in this genus, we here report a high-quality chromosome-level genome assembly of M. veneriformis. Our approach yielded a 939.32 Mb assembled genome with an N50 contig length of 7,977.84 kb. Hi-C scaffolding of the genome resulted in assembly of 19 pseudochromosomes. Repetitive elements made up ∼51.79% of the genome assembly. A total of 29,315 protein-coding genes (PCGs) were predicted in M. veneriformis. Construction of a genome-level phylogenetic tree demonstrated that M. veneriformis and Ruditapes philippinarum diverged around 231 million years ago (MYA). Inter-species comparisons revealed that 493 gene families have undergone expansion and 449 have undergone contraction in the M. veneriformis genome. Chromosome-based macrosynteny analysis revealed a high degree of synteny between the 19 chromosomes of M. veneriformis and those of Patinopecten yessoensis. These results suggested that M. veneriformis has a similar karyotype to that of P. yessoensis, and that a highly conserved 19-chromosome karyotype was formed in the early differentiation stages of bivalves. In summary, the genomic resources generated in this work serve as a valuable reference for investigating the molecular mechanisms underlying biological functions in M. veneriformis and will facilitate future genetic improvement and disease treatment in this economically important species. Furthermore, the assembled genome greatly improves our understanding of early genomic evolution of the Bivalvia.

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.

Massive expansion of P-selectin genes in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria: evidence from comparative genomics of Bivalvia

Background

P-selectin is a molecule participating in the inflammatory response through mediating cellular adhesion and essential for wound repair. However, studies regarding P-selectin in Bivalvia are rare. This study identified 90 P-selectin genes among nine bivalve genomes and classified them into 4 subfamilies according to phylogenetic analysis.

Results

Notable P-selectin gene expansion was observed in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria. The synteny analysis revealed that P-selectin gene expansion was mostly caused by tandem duplication. In addition, the expression profiles of P-selectin genes in S. constricta showed that many P-selectins were specifically highly expressed in the gills, and the P-selectin expression patterns changed dramatically under low salt stress and ammonia nitrogen stress.

Conclusions

The massive expansion of P-selectins may facilitate the tolerance to environmental stresses. This study sheds light on the characterizations and expression profiles of P-selectin genes in Bivalvia and provides an integrated framework for further investigation of the role of P-selectins in the environmental tolerance of bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08861-6.

Pooled resequencing of larvae and adults reveals genomic variations associated with Ostreid herpesvirus 1 resistance in the Pacific oyster Crassostrea gigas

The Ostreid herpesvirus 1 (OsHV-1) is a lethal pathogen of the Pacific oyster (Crassostrea gigas), an important aquaculture species. To understand the genetic architecture of the defense against the pathogen, we studied genomic variations associated with herpesvirus-caused mortalities by pooled whole-genome resequencing of before and after-mortality larval samples as well as dead and surviving adults from a viral challenge. Analysis of the resequencing data identified 5,271 SNPs and 1,883 genomic regions covering 3,111 genes in larvae, and 18,692 SNPs and 28,314 regions covering 4,863 genes in adults that were significantly associated with herpesvirus-caused mortalities. Only 1,653 of the implicated genes were shared by larvae and adults, suggesting that the antiviral response or resistance in larvae and adults involves different sets of genes or differentiated members of expanded gene families. Combined analyses with previous transcriptomic data from challenge experiments revealed that transcription of many mortality-associated genes was also significantly upregulated by herpesvirus infection confirming their importance in antiviral response. Key immune response genes especially those encoding antiviral receptors such as TLRs and RLRs displayed strong association between variation in regulatory region and herpesvirus-caused mortality, suggesting they may confer resistance through transcriptional modulation. These results point to previously undescribed genetic mechanisms for disease resistance at different developmental stages and provide candidate polymorphisms and genes that are valuable for understanding antiviral immune responses and breeding for herpesvirus resistance.

Genome-Wide Association and Genomic Prediction of Growth Traits in the European Flat Oyster (Ostrea edulis)

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes.

Genome-wide identification and characterization of superoxide dismutases in four oyster species reveals functional differentiation in response to biotic and abiotic stress

Background

Oysters inhabit in the intertidal zone and may be suffered from environmental stresses, which can increase the production of reactive oxygen species (ROS), resulting in mass mortality. Superoxide dismutases (SODs) protect oysters from ROS damage through different mechanisms compared with vertebrates. However, the molecular and functional differentiation in oyster SODs were rarely analyzed.

Result

In this study, a total of 13, 13, 10, and 8 candidate SODs were identified in the genome of Crassostrea gigas, Crassostrea virginica, Crassostrea hongkongensis, and Saccostrea glomerata respectively. The domain composition, gene structure, subcellular locations, conserved ligands, and cis-elements elucidated the SODs into five groups (Mn-SODs, Cu-only-SODs, Cu/Zn ion ligand Cu/Zn-SOD with enzyme activity, Zn-only-SODs, and no ligand metal ions Cu/Zn-SODs). For single domain Cu/Zn-SODs, only one cytosolic Cu/Zn-SOD (cg_XM_034479061.1) may conserve enzymatic activity while most extracellular Cu/Zn-SOD proteins appeared to lose SOD enzyme activity according to conserved ligand amino acid analysis and expression pattern under biotic and abiotic stress in C. gigas. Further, multi-domain-SODs were identified and some of them were expressed in response to biotic and abiotic stressors in C. gigas. Moreover, the expression patterns of these genes varied in response to different stressors, which may be due to the cis-elements in the gene promoter.

Conclusion

These findings revealed the most extracellular Cu/Zn-SOD proteins appeared to lose SOD enzyme activity in oysters. Further, our study revealed that only one cytosolic Cu/Zn-SOD (cg_XM_034479061.1) may conserve enzymatic activity of SOD. Moreover, the expression patterns of these genes varied in response to different stressors, which may be due to the cis-elements in the promoter. This study provides important insights into the mechanisms through which oysters adapt to harsh intertidal conditions, as well as potential biomarkers of stress response in related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08610-9.

DNA methylation differences between male and female gonads of the oyster reveal the role of epigenetics in sex determination

DNA methylation involved in sex determination mechanism by regulating gene expression related to sex determination networks are common in vertebrates. However, the mechanism linking epigenetics in invertebrates and sex determination has remained elusive. Here, methylome of the male and female gonads in the oyster Crassostrea gigas were conducted to explore the role of epigenetics in invertebrate sex determination. Comparative analysis of gonadal DNA methylation of females and males revealed that male gonads displayed a higher level of DNA methylation and a greater number of hypermethylated genes. Luxury genes presented hypomethylation, while housekeeping genes got hypermethylation. Genes in the conserved signaling pathways, rather than the key master genes in the sex determination pathway, were the major targets of substantial DNA methylation modification. The negative correlation of expression and promoter methylation in the diacylglycerol kinase delta gene (Dgkd) - a ubiquitously expressed gene - indicated DNA methylation may fine turn the expression of Dgkd and be involved in the process of sex determination. Dgkd can be used as an epigenetic marker to distinguish male C. gigas based on the different methylation regions in the promoter region. The results suggest that DNA methylation mechanisms played potential functional impacts in the sex determination in oysters, which is helpful to deepen the understanding of sex determination in invertebrate.

Genome-wide analysis of the hard clam mitogen-activated protein kinase kinase gene family and their transcriptional profiles under abiotic stress

Mitogen-activated protein kinase kinase (MAPKK) was the hub component of the Mitogen-activated protein kinase (MAPK) signaling pathway and played an important role in the cellular response to environmental stress. In this study, we identified five MmMAPKK genes in hard clam Mercenaria mercenaria and found that all MmMAPKK genes contain a conserved protein kinase domain. The MmMAPKK genes derived from dispersed duplication were unevenly distributed in three chromosomes. Although the genome size was highly variable among different bivalve mollusks, the number of MAPKK genes was relatively stable. Phylogenetic analysis showed that bivalve MAPKK was divided into five clades, and amino acid sequences of MAPKK from the same clade consisted of similar conserved motifs. The syntenic analysis demonstrated that MmMAPKKs had the highest number of homologous gene pairs with Cyclina sinensis. MmMAPKKs were ubiquitously expressed in all examined tissues, and all MmMAPKK genes were highly expressed in the ovary. MmMAPKK genes showed stress-specific expression under envirionmental stress. MmMAPKK7 showed an upregulated in heat and heat plus hypoxia stress while MmMAPKK1 showed an upregulated in hypoxic stress groups. Dynamic changes of MmMAPKK7, MmMAPKK6 and MmMAPKK1 in hemocytes were observed in response to air exposure. MmMAPKK4 significantly downregulated after air exposure for five days. MmMAPKK7 and MmMAPKK6 might participate in adaptation to low salinity stress. Our results provided useful information about MAPKK and laid a foundation for further studies on MAPKK evolution in the bivalve.

Chromosome-level genome and population genomic analysis provide insights into the evolution and environmental adaptation of Jinjiang oyster Crassostrea ariakensis

The Jinjiang oyster Crassostrea ariakensis, naturally distributing in estuarine regions with low salinity, is an important economic and ecological species in China. However, studies on its genomics and population genetics remain lacking. Here, we assembled the chromosome-level genome of a female C. ariakensis and re-sequenced 261 individuals from five locations in China representing three typical habitats. The C. ariakensis genome was 662.9 Mb with contig N50 length of 5.9 Mb using PacBio HiFi-CCS long reads, and 99.83% sequences were anchored onto 10 pseudochromosomes using Hi-C data. A total of 26,354 protein-coding genes were predicted. We identified three significantly expanded gene families which are closely associated with osmotic pressure regulation, including CDO, SLC13 and SDR. Population structure analysis revealed that the C. ariakensis from five locations were clustered into three typical groups (northern, southern and Shanghai) (K = 3) and their phylogenetic relationship was consistently correlated to their geographical distribution. Furtherly, the differentiation between northern and southern groups was clearly demonstrated by estimated population differentiation coefficient (F_ST  = 0.1154), and the PSMC distribution showed the two groups of effective population size separated at 0.1 Ma. Meanwhile gene flow from southern to Shanghai was detected. Selective sweep analysis between northern and southern group detected genes associated with heat response and salinity adaptation. This study could provide valuable genomic resources and information for further research on the molecular evolution, genetic breeding, biological function and evolutionary adaptation of C. ariakensis.

Genome of the estuarine oyster provides insights into climate impact and adaptive plasticity

Understanding the roles of genetic divergence and phenotypic plasticity in adaptation is central to evolutionary biology and important for assessing adaptive potential of species under climate change. Analysis of a chromosome-level assembly and resequencing of individuals across wide latitude distribution in the estuarine oyster (Crassostrea ariakensis) revealed unexpectedly low genomic diversity and population structures shaped by historical glaciation, geological events and oceanographic forces. Strong selection signals were detected in genes responding to temperature and salinity stress, especially of the expanded solute carrier families, highlighting the importance of gene expansion in environmental adaptation. Genes exhibiting high plasticity showed strong selection in upstream regulatory regions that modulate transcription, indicating selection favoring plasticity. Our findings suggest that genomic variation and population structure in marine bivalves are heavily influenced by climate history and physical forces, and gene expansion and selection may enhance phenotypic plasticity that is critical for the adaptation to rapidly changing environments.

Overt and concealed genetic loads revealed by QTL mapping of genotype-dependent viability in the Pacific oyster Crassostrea gigas

Understanding the genetic bases of inbreeding depression, heterosis, and genetic load is integral to understanding how genetic diversity is maintained in natural populations. The Pacific oyster Crassostrea gigas, like many long-lived plants, has high fecundity and high early mortality (type-III survivorship), manifesting a large, overt, genetic load; the oyster harbors an even greater concealed genetic load revealed by inbreeding. Here, we map viability QTL (vQTL) in six interrelated F2 oyster families, using high-density linkage maps of single nucleotide polymorphisms generated by genotyping-by-sequencing (GBS) methods. Altogether, we detect 70 vQTL and provisionally infer 89 causal mutations, 11 to 20 per family. Genetic mortality caused by independent (unlinked) vQTL ranges from 94.2% to 97.8% across families, consistent with previous reports. High-density maps provide better resolution of genetic mechanisms, however. Models of one causal mutation present in both identical-by-descent (IBD) homozygotes and heterozygotes fit genotype frequencies at 37 vQTL; consistent with the mutation-selection balance theory of genetic load, 20 are highly deleterious, completely recessive mutations and 17 are less deleterious, partially dominant mutations. Another 22 vQTL require pairs of recessive or partially dominant causal mutations, half showing selection against recessive mutations linked in repulsion, producing pseudo-overdominance. Only eight vQTL appear to support the overdominance theory of genetic load, with deficiencies of both IBD homozygotes, but at least four of these are likely caused by pseudo-overdominance. Evidence for epistasis is absent. A high mutation rate, random genetic drift, and pseudo-overdominance may explain both the oyster's extremely high genetic diversity and a high genetic load maintained primarily by mutation-selection balance.

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

The Pacific oyster (Crassostrea gigas) genome is highly polymorphic and affluent in structural variations (SVs), a significant source of genetic variation underlying inter-individual differences. Here, we used two genome assemblies and 535 individuals of genome re-sequencing data to construct a comprehensive landscape of structural variations in the Pacific oyster. Through whole-genome alignment, 11,087 short SVs and 11,561 copy number variations (CNVs) were identified. While analysis of re-sequencing data revealed 511,170 short SVs and 979,486 CNVs, a total of 63,100 short SVs and 58,182 CNVs were identified in at least 20 samples and regarded as common variations. Based on the common short SVs, both Fst and Pi ratio statistical methods were employed to detect the selective sweeps between 20 oyster individuals from the fast-growing strain and 20 individuals from their corresponding wild population. A total of 514 overlapped regions (8.76 Mb), containing 746 candidate genes, were identified by both approaches, in addition with 103 genes within 61 common CNVs only detected in the fast-growing strains. The GO enrichment and KEGG pathway analysis indicated that the identified candidate genes were mostly associated with apical part of cell and were significantly enriched in several metabolism-related pathways, including tryptophan metabolism and histidine metabolism. This work provided a comprehensive landscape of SVs and revealed their responses to selection, which will be valuable for further investigations on genome evolution under selection in the oysters.

Noncoding variation and transcriptional plasticity promote thermal adaptation in oysters by altering energy metabolism

Genetic variation and phenotypic plasticity are both important to adaptive evolution. However, how they act together on particular traits remains poorly understood. Here, we integrated phenotypic, genomic, and transcriptomic data from two allopatric but closely related congeneric oyster species, Crassostrea angulata from southern/warm environments and Crassostrea gigas from northern/cold environments, to investigate the roles of genetic divergence and plasticity in thermal adaptation. Reciprocal transplantation experiments showed that both species had higher fitness in their native habitats than in nonnative environments, indicating strong adaptive divergence. The southern species evolved higher transcriptional plasticity, and the plasticity was adaptive, suggesting that increased plasticity is important for thermal adaptation to warm climates. Genome-wide comparisons between the two species revealed that genes under selection tended to respond to environmental changes and showed higher sequence divergence in noncoding regions. All genes under selection and related to energy metabolism exhibited habitat-specific expression with genes involved in ATP production and lipid catabolism highly expressed in warm/southern habitats, and genes involved in ATP consumption and lipid synthesis were highly expressed in cold/northern habitats. The gene for acyl-CoA desaturase, a key enzyme for lipid synthesis, showed strong selective sweep in the upstream noncoding region and lower transcription in the southern species. These results were further supported by the lower free fatty acid (FFA) but higher ATP content in southern species and habitat, pointing to significance of ATP/FFA trade-off. Our findings provide evidence that noncoding variation and transcriptional plasticity play important roles in shaping energy metabolism for thermal adaptation in oysters.