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

OysterDB: A Genome Database for Ostreidae

The molluscan family Ostreidae, commonly known as oysters, is an important molluscan group due to its economic and ecological importance. In recent years, an abundance of genomic data of Ostreidae species has been generated and available in public domain. However, there is still a lack of a high-efficiency database platform to store and distribute these data with comprehensive tools. In this study, we developed an oyster genome database (OysterDB) to consolidate oyster genomic data. This database includes eight oyster genomes and 208,923 protein-coding gene annotations. Bioinformatic tools, such as BLAST and JBrowse, are integrated into the database to provide a user-friendly platform for homologous sequence searching, visualization of genomes, and screen for candidate gene information. Moreover, OysterDB will be continuously updated with ever-growing oyster genomic resources and facilitate future studies for comparative and functional genomic analysis of oysters ( http://oysterdb.com.cn/ ).

Physiological and transcriptomic analysis provides new insights into osmoregulation mechanism of Ruditapes philippinarum under low and high salinity stress

The Manila clam (Ruditapes philippinarum) is a commercially important marine bivalve, which inhabits the estuarine and mudflat areas. The osmoregulation is of great significance for molluscs adaptation to salinity fluctuations. In this study, we investigated the effects of low salinity (10 psu) and high salinity (40 psu) stress on survival and osmoregulation of the R. philippinarum. The results of physiological parameters showed that the ion (Na+, K+, Cl-) concentrations and Na+/K+-ATPase (NKA) activity of R. philippinarum decreased significantly under low salinity stress, but increased significantly under high salinity stress, indicating that there are differences in physiological adaptation of osmoregulation of R. philippinarum. In addition, we conducted the transcriptome analysis in the gills of R. philippinarum exposed to low (10 psu) and high (40 psu) salinity challenge for 48 h using RNA-seq technology. A total of 153 and 640 differentially expressed genes (DEGs) were identified in the low salinity (LS) group and high salinity (HS) group, respectively. The immune (IAP, TLR6, C1QL4, Ank3), ion transport (Slc34a2, SLC39A14), energy metabolism (PCK1, LDLRA, ACOX1) and DNA damage repair-related genes (Gadd45g, HSP70B2, GATA4) as well as FoxO, protein processing in endoplasmic reticulum and endocytosis pathways were involved in osmoregulation under low salinity stress of R. philippinarum. Conversely, the ion transport (SLC6A7, SLC6A9, SLC6A14, TRPM2), amino acid metabolism (GS, TauD, ABAT, ALDH4A1) and immune-related genes (MAP2K6, BIRC7A, CTSK, GVIN1), and amino acid metabolism pathways (beta-Alanine, Alanine, aspartate and glutamate, Glutathione) were involved in the process of osmoregulation under high salinity stress. The results obtained here revealed the difference of osmoregulation mechanism of R. philippinarum under low and high salinity stress through physiological and molecular levels. This study contributes to the assessment of salinity adaptation of bivalves in the context of climate change and provides useful information for marine resource conservation and aquaculture.

Genome-wide identification of the HSP70 genes in Pacific oyster Magallana gigas and their response to heat stress

Heat shock protein 70 (HSP70), the most prominent and well-characterized stress protein in animals, plays an important role in assisting animals in responding to various adverse conditions. In the present study, a total of 113 HSP70 gene family members were identified in the updated genome of Magallana gigas (designated MgHSP70) (previously known as Crassostrea gigas). There were 75, 12, 11, and 8 HSP70s located in the cytoplasm, nucleus, mitochondria, and endoplasmic reticulum, respectively, and 7 HSP70s were located in both the nucleus and cytoplasm. Among 113 MgHSP70 genes, 107 were unevenly distributed in 8 chromosomes of M. gigas with the greatest number in chromosome 07 (61 genes, 57.01%). The MgHSP70 gene family members were mainly assigned into five clusters, among which the HSPa12 subfamily underwent lineage-specific expansion, consisting of 89 members. A total of 68 MgHSP70 genes (60.18%) were tandemly duplicated and formed 30 gene pairs, among which 14 gene pairs were under strong positive selection. In general, the expression of MgHSP70s was tissue-specific, with the highest expression in labial palp and gill and the lowest expression in adductor muscle and hemocytes. There were 35, 31, and 47 significantly upregulated genes at 6, 12, and 24 h after heat shock treatment (28 °C), respectively. The expression patterns of different tandemly duplicated genes exhibited distinct characteristics after shock treatment, indicating that these genes may have different functions. Nevertheless, genes within the same tandemly duplicated group exhibit similar expression patterns. Most of the tandemly duplicated HSP70 gene pairs showed the highest expression levels at 24 h. This study provides a comprehensive description of the MgHSP70 gene family in M. gigas and offers valuable insights into the functions of HSP70 in the mollusc adaptation of oysters to environmental stress.

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.

Morphological and Molecular Analysis Identified a Subspecies of Crassostrea ariakensis (Fujita, 1913) along the Coast of Asia

Crassostrea ariakensis (Fujita, 1913) is one of the most important economic and ecological oysters that is naturally distributed along the coast of Asia, separated by the Yangtze River estuary. They are usually compared as different populations, while there is no consensus on whether C. ariakensis in northern and southern areas should be considered as two species or subspecies. Here, we analyzed morphological characteristics, COI, 16s rRNA, mitogenome sequences, and species delimitation analysis (ASAP and PTP) to resolve the intraspecific taxonomic status of the C. ariakensis. Phylogenetic and ASAP analysis highlight that C. ariakensis was divided into N-type and S-type. PTP was unable to differentiate between the two types of C. ariakensis. The divergence time of N-type and S-type C. ariakinsis is estimated to be 1.6 Mya, using the relaxed uncorrelated lognormal clock method. Additionally, significant morphological differences exist between the two groups in terms of the adductor muscle scar color. Despite these differences, the COI (0.6%) and 16S rRNA (0.6%) genetic distance differences between N-type and S-type C. ariakensis has not yet reached the interspecific level. These results suggest that N-type and S-type C. ariakensis should be treated as different subspecies and renamed as C. ariakensis ariakensis subsp. nov and C. ariakensis meridioyangtzensis subsp. nov.

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.

The prokaryotic and eukaryotic microbiome of Pacific oyster spat is shaped by ocean warming but not acidification

Pacific oysters (Magallana gigas, a.k.a. Crassostrea gigas), the most widely farmed oysters, are under threat from climate change and emerging pathogens. In part, their resilience may be affected by their microbiome, which, in turn, may be influenced by ocean warming and acidification. To understand these impacts, we exposed early-development Pacific oyster spat to different temperatures (18°C and 24°C) and pCO2 levels (800, 1,600, and 2,800 µatm) in a fully crossed design for 3 weeks. Under all conditions, the microbiome changed over time, with a large decrease in the relative abundance of potentially pathogenic ciliates (Uronema marinum) in all treatments with time. The microbiome composition differed significantly with temperature, but not acidification, indicating that Pacific oyster spat microbiomes can be altered by ocean warming but is resilient to ocean acidification in our experiments. Microbial taxa differed in relative abundance with temperature, implying different adaptive strategies and ecological specializations among microorganisms. Additionally, a small proportion (~0.2% of the total taxa) of the relatively abundant microbial taxa were core constituents (>50% occurrence among samples) across different temperatures, pCO2 levels, or time. Some taxa, including A4b bacteria and members of the family Saprospiraceae in the phyla Chloroflexi (syn. Chloroflexota) and Bacteroidetes (syn. Bacteroidota), respectively, as well as protists in the genera Labyrinthula and Aplanochytrium in the class Labyrinthulomycetes, and Pseudoperkinsus tapetis in the class Ichthyosporea were core constituents across temperatures, pCO2 levels, and time, suggesting that they play an important, albeit unknown, role in maintaining the structural and functional stability of the Pacific oyster spat microbiome in response to ocean warming and acidification. These findings highlight the flexibility of the spat microbiome to environmental changes.IMPORTANCEPacific oysters are the most economically important and widely farmed species of oyster, and their production depends on healthy oyster spat. In turn, spat health and productivity are affected by the associated microbiota; yet, studies have not scrutinized the effects of temperature and pCO2 on the prokaryotic and eukaryotic microbiomes of spat. Here, we show that both the prokaryotic and, for the first time, eukaryotic microbiome of Pacific oyster spat are surprisingly resilient to changes in acidification, but sensitive to ocean warming. The findings have potential implications for oyster survival amid climate change and underscore the need to understand temperature and pCO2 effects on the microbiome and the cascading effects on oyster health and productivity.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis)

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

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.

Comparative Transcriptome and WGCNA Analysis Reveal Molecular Responses to Salinity Change in Larvae of the Iwagaki Oyster Crassostrea Nippona

The Iwagaki oyster Crassostrea nippona is an important aquaculture species with significant potential for large-scale oyster farming. It is susceptible to the fluctuated salinity in the coastal area. In this study, we compared the transcriptome of Crassostrea nippona larvae under variant conditions with low-salinity stress (28, 20, 15, 10, and 5 practical salinity units (psu)) for 24 h. KEGG enrichment analysis of differentially expressed genes (DEGs) from pairwise comparisons identified several free amino acid metabolism pathway (taurine and hypotaurine, arginine and proline, glycine, and beta-alanine) contributing to the salinity change adaptation and activated "lysosome" and "apoptosis" pathway in response to the low-salinity stress (10 and 5 psu). Trend analysis revealed sustained upregulation of transmembrane transport-related genes (such as SLC family) and downregulation of ribosomal protein synthesis genes faced with decreasing salinities. In addition, 9 biomarkers in response to low-salinity stress were identified through weighted gene co-expression network analysis (WGCNA) and validated by qRT-PCR. Our transcriptome analysis provides a comprehensive view of the molecular mechanisms and regulatory networks underlying the adaptive responses of oyster larvae to hypo-salinity conditions. These findings contribute to our understanding of the complex biological processes involved in oyster resilience and adaptation to changing environmental conditions.

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Genome architecture and selective signals compensatorily shape plastic response to a new environment

Transcriptional plasticity interacts with natural selection in complex ways and is crucial for the survival of species under rapid climate change. How 3D genome architecture affects transcriptional plasticity and its interaction with genetic adaptation are unclear. We transplanted estuarine oysters to a new environment and found that genes located in active chromatin regions exhibited greater transcriptional plasticity, and changes in these regions were negatively correlated with selective signals. This indicates a trade-off between 3D active regions and selective signals in shaping plastic responses to a new environment. Specifically, a mutation, lincRNA, and changes in the accessibility of a distal enhancer potentially affect its interaction with the ManⅡa gene, which regulates the muscle function and survival of oysters. Our findings reveal that 3D genome architecture compensates for the role of genetic adaptation in environmental response to new environments and provide insights into synergetic genetic and epigenetic interactions critical for fitness-related trait and survival in a model marine species.

Energetic budget of diploid and triploid eastern oysters during a summer die-off

Triploid oysters are widely used in off-bottom aquaculture of eastern oysters, Crassostrea virginica. However, farmers of the Gulf of Mexico (GoM) and Atlantic coast estuaries have observed unresolved, late-spring die-offs of triploid oysters, threatening the sustainability of triploid aquaculture. To investigate this, the physiological processes underlying oyster growth (e.g., feeding, respiration) and mortality of one-year-old diploid and triploid oysters were compared in early summer following an uptick in mortality. It was predicted that higher triploid mortality was the result of energetic imbalances (increased metabolic demands and decreased feeding behavior). Oyster clearance rates, percentage of time valves were open, absorption efficiency, oxygen consumption rates (basal and routine), ammonia excretion rate were measured in the laboratory and scope for growth was calculated. In addition, their condition index, gametogenic stage, Perkinsus marinus infection level, and mortality were measured. Mortality of triploids in the laboratory was greater than for diploids, mirroring mortality observed in a related field study. The physiological parameters measured, however, could not explain triploid mortality. Scope for growth, condition index, and clearance rates of triploids were greater than for diploids, suggesting sufficient energy reserves, while all other measurements where similar between the ploidies. It remains to be determined whether mortality could be caused from disruption of energy homeostasis at the cellular level.

Cardiac performance and heart gene network provide dynamic responses of bay scallop Argopecten irradians irradians exposure to marine heatwaves

The increased frequency of marine heat waves (MHWs) caused by global climate change is predicted to threaten the survival of economic bivalves, therefore having severely adverse effects on local ecological communities and aquaculture production. However, the study of scallops facing MHWs is still scarce, particularly in the scallop Argopecten irradians irradians, which has a significant share of "blue foods" in northern China. In the present study, bay scallop heart was selected to detect its cardiac performance, oxidative impairment and dynamic molecular responses, accompanied by assessing survival variations of individuals in the simulated scenario of MWHs (32 °C) with different time points (0 h, 6 h, 12 h, 24 h, 3 d, 6 d and 10 d). Notably, cardiac indices heart rate (HR), heart amplitude (HA), rate-amplitude product (RAP) and antioxidant enzyme activities superoxide dismutase (SOD) and catalase (CAT) all peaked at 24 h but sharply dropped on 3 d, coinciding with mortality. Transcriptome analysis revealed that the heart actively defended against heat stress at the acute stage (<24 h) via energy supply, misfolded proteins correction and enhanced signal transduction, whereas regulation of the defense response and apoptotic process combined with twice transcription initiation were the dominant responses at the chronic stage (3-10 d). In particular, HSP70 (heat shock protein 70), HSP90 and CALR (calreticulin) in the endoplasmic reticulum were identified as the hub genes (top 5 %) in the HR-associated module via WGCNA (weighted gene co-expression network analysis) trait-module analysis, followed by characterization of their family members and diverse expression patterns under heat exposure. Furthermore, RNAi-mediated knockdown of CALR expression (after 24 h) significantly weakened the thermotolerance of scallops, as evidenced by a drop of 1.31 °C in ABT (Arrhenius break temperature) between the siRNA-injected group and the control group. Our findings elucidated the dynamic molecular responses at the transcriptome level and verified the cardiac functions of CALR in bay scallops confronted with stimulated MHWs.

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.

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.

Adaptive divergence and underlying mechanisms in response to salinity gradients between two Crassostrea oysters revealed by phenotypic and transcriptomic analyses

Comparing the responses of closely related species to environmental changes is an efficient method to explore adaptive divergence, for a better understanding of the adaptive evolution of marine species under rapidly changing climates. Oysters are keystone species thrive in intertidal and estuarine areas where frequent environmental disturbance occurs including fluctuant salinity. The evolutionary divergence of two sister species of sympatric estuarine oysters, Crassostrea hongkongensis and Crassostrea ariakensis, in response to euryhaline habitats on phenotypes and gene expression, and the relative contribution of species effect, environment effect, and their interaction to the divergence were explored. After a 2-month outplanting at high- and low-salinity locations in the same estuary, the high growth rate, percent survival, and high tolerance indicated by physiological parameters suggested that the fitness of C. ariakensis was higher under high-salinity conditions and that of C. hongkongensis was higher under low-salinity conditions. Moreover, a transcriptomic analysis showed the two species exhibited differentiated transcriptional expression in high- and low-salinity habitats, largely caused by the species effect. Several of the important pathways enriched in divergent genes between species were also salinity-responsive pathways. Specifically, the pyruvate and taurine metabolism pathway and several solute carriers may contribute to the hyperosmotic adaptation of C. ariakensis, and some solute carriers may contribute to the hypoosmotic adaptation of C. hongkongensis. Our findings provide insights into the phenotypic and molecular mechanisms underlying salinity adaptation in marine mollusks, which will facilitate the assessment of the adaptive capacity of marine species in the context of climate change and will also provide practical information for marine resource conservation and aquaculture.

Population Genomics, Transcriptional Response to Heat Shock, and Gut Microbiota of the Hong Kong Oyster Magallana hongkongensis

The Hong Kong oyster Magallana hongkongensis, previously known as Crassostrea hongkongensis, is a true oyster species native to the estuarine-coast of the Pearl River Delta in southern China. The species—with scientific, ecological, cultural, and nutritional importance—has been farmed for hundreds of years. However, there is only limited information on its genetics, stress adaptation mechanisms, and gut microbiota, restricting the sustainable production and use of oyster resources. Here, we present population structure analysis on M. hongkongensis oysters collected from Deep Bay and Lantau Island in Hong Kong, as well as transcriptome analysis on heat shock responses and the gut microbiota profile of M. hongkongensis oysters collected from Deep Bay. Single nucleotide polymorphisms (SNPs), including those on the homeobox genes and heat shock protein genes, were revealed by the whole genome resequencing. Transcriptomes of oysters incubated at 25 °C and 32 °C for 24 h were sequenced which revealed the heat-induced regulation of heat shock protein pathway genes. Furthermore, the gut microbe community was detected by 16S rRNA sequencing which identified Cyanobacteria, Proteobacteria and Spirochaetes as the most abundant phyla. This study reveals the molecular basis for the adaptation of the oyster M. hongkongensis to environmental conditions.