Thermotolerance Divergence Revealed by the Physiological and Molecular Responses in Two Oyster Subspecies of Crassostrea gigas in China

Bivalves under extreme weather events: A comparative study of five economically important species in the South China sea during marine heatwaves

Marine heatwaves (MHWs) are increasing in frequency and intensity, threatening marine organisms and ecosystems they support. Yet, little is known about impacts of intensifying MHWs on ecologically and economically important bivalves cultured in the South China Sea. Here, we compared survival and physiological responses of five bivalve species, Pinctada fucata, Crassostrea angulata, Perna viridis, Argopecten irradians and Paphia undulata, to two consecutive MHWs events (3 days of thermal exposure to + 4 °C or + 8 °C, following 3 days of recovery under ambient conditions). While P. fucata, P. viridis, and P. undulata are native to the South China Sea region, C. angulata and A. irradians are not. Individuals of P. fucata, C. angulata and P. viridis had higher stress tolerance to MHWs than A. irradians and P. undulata, the latter already experiencing 100% mortality under +8 °C conditions during the first event. With increasing intensity of MHWs, standard metabolic rates of all five species increased significantly, in line with significant depressions of function-related energy-metabolizing enzymes (CMA, NKA, and T-ATP). Likewise, activities of antioxidant enzymes (SOD, CAT, and MDA) and shell mineralization-related enzymes (AKP and ACP) responded significantly to MHWs, despite species-specific performances observed. These findings demonstrate that some bivalve species can likely fail to accommodate intensifying MHWs events in the South China Sea, but some may persist. If this is the case, then one would expect substantial loss of fitness in bivalve aquaculture in the South China Sea under intensifying MHWs conditions.

Transcriptomic analysis of mRNA and miRNA reveals new insights into the regulatory mechanisms of Anadara granosa responses to heat stress

Temperature fluctuations resulting from climate change and global warming pose significant threats to various species. The blood clam, Anadara granosa, a commercially important marine bivalve, predominantly inhabits intertidal mudflats that are especially susceptible to elevated temperatures. This vulnerability has led to noticeable declines in the survival rates of A. granosa larvae, accompanied by an increase in malformations. Despite these observable trends, there is a lack of comprehensive research on the regulatory mechanisms underlying A. granosa's responses to heat stress. In this study, we examined the survival rates of A. granosa under varying high temperature conditions, selecting 34 °C as heat stress temperature. Enzyme activity assays have shed light on A. granosa's adaptive response to heat stress, revealing its ability to maintain redox balance and transition from aerobic to anaerobic metabolic pathways. Subsequently, mRNA and miRNA transcriptome analyses were conducted, elucidating several key responses of A. granosa to heat stress. These responses include the upregulation of transcription and protein synthesis, downregulation of proteasome activity, and metabolic pattern adjustments. Furthermore, we identified miRNA-mRNA networks implicated in heat stress responses, potentially serving as valuable candidate markers for A. granosa's heat stress response. Notably, we validated the involvement of agr-miR-3199 in A. granosa's heat stress response through its regulation of the target gene Foxj1. These findings not only deepen our understanding of the molecular mechanisms underlying the blood clam's response to heat stress but also offer valuable insights for enhancing heat stress resilience in the blood clam aquaculture industry. Moreover, they contribute to improved cultivation strategies for molluscs in the face of global warming and have significant implications for the conservation of marine resources and the preservation of ecological balance.

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.

Temperature effects on plasmalogen profile and quality characteristics in Pacific oyster (Crassostrea gigas) during depuration

The quality of Pacific oyster (Crassostrea gigas) can be affected by many factors during depuration, in which temperature is the major element. In this study, we aim to determine the quality and plasmalogen changes in C. gigas depurated at different temperatures. The quality was significantly affected by temperature, represented by varying survival rate, glycogen content, total antioxidant capacity, alkaline phosphatase activity between control and stressed groups. Targeted MS analysis demonstrated that plasmalogen profile was significantly changed during depuration with PUFA-containing plasmalogen species being most affected by temperature. Proteomics analysis and gene expression assay further verified that plasmalogen metabolism is regulated by temperature, specifically, the plasmalogen synthesis enzyme EPT1 was significantly downregulated by high temperature and four plasmalogen-related genes (GPDH, PEDS, Pex11, and PLD1) were transcriptionally regulated. The positive correlations between the plasmalogen level and quality characteristics suggested plasmalogen could be regarded as a quality indicator of oysters during depuration.

Distinct metabolic responses to thermal stress between invasive freshwater turtle Trachemys scripta elegans and native freshwater turtles in China

Different responses or tolerance to thermal stress between invasive and native species can affect the outcome of interactions between climate change and biological invasion. However, knowledge about the physiological mechanisms that modulate the interspecific differences in thermal tolerance is limited. The present study analyzes the metabolic responses to thermal stress by the globally invasive turtle, Trachemys scripta elegans, as compared with two co-occurring native turtle species in China, Pelodiscus sinensis and Mauremys reevesii. Changes in metabolite contents and the expression or enzyme activities of genes involved in energy sensing, glucose metabolism, lipid metabolism, and tricarboxylic acid (TCA) cycle after exposure to gradient temperatures were assessed in turtle juveniles. Invasive and native turtles showed distinct metabolic responses to thermal stress. T. scripta elegans showed greater transcriptional regulation of energy sensors than the native turtles. Enhanced anaerobic metabolism was needed by all three species under extreme heat conditions, but phosphoenolpyruvate carboxykinase and lactate dehydrogenase in the invader showed stronger upregulation or stable responses than the native species, which showed inhibition by high temperatures. These contrasts were pronounced in the muscles of the three species. Regulation of lipid metabolism was observed in both T. scripta elegans and P. sinensis but not in M. reevesii under thermal stress. Thermal stress did not inhibit the TCA cycle in turtles. Different metabolic responses to thermal stress may contribute to interspecific differences in thermal tolerance. Overall, our study further suggested the potential role of physiological differences in mediating interactions between climate change and biological invasion.

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.

Redox control of antioxidants, metabolism, immunity, and development at the core of stress adaptation of the oyster Crassostrea gigas to the dynamic intertidal environment

This review uses the marine bivalve Crassostrea gigas to highlight redox reactions and control systems in species living in dynamic intertidal environments. Intertidal species face daily and seasonal environmental variability, including temperature, oxygen, salinity, and nutritional changes. Increasing anthropogenic pressure can bring pollutants and pathogens as additional stressors. Surprisingly, C. gigas demonstrates impressive adaptability to most of these challenges. We explore how ROS production, antioxidant protection, redox signaling, and metabolic adjustments can shed light on how redox biology supports oyster survival in harsh conditions. The review provides (i) a brief summary of shared redox sensing processes in metazoan; (ii) an overview of unique characteristics of the C. gigas intertidal habitat and the suitability of this species as a model organism; (iii) insights into the redox biology of C. gigas, including ROS sources, signaling pathways, ROS-scavenging systems, and thiol-containing proteins; and examples of (iv) hot topics that are underdeveloped in bivalve research linking redox biology with immunometabolism, physioxia, and development. Given its plasticity to environmental changes, C. gigas is a valuable model for studying the role of redox biology in the adaptation to harsh habitats, potentially providing novel insights for basic and applied studies in marine and comparative biochemistry and physiology.

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.

Triploid Pacific oysters exhibit stress response dysregulation and elevated mortality following heatwaves

Polyploidy has been suggested to negatively impact environmental stress tolerance, resulting in increased susceptibility to extreme climate events. In this study, we compared the genomic and physiological response of diploid (2n) and triploid (3n) Pacific oysters (Crassostrea gigas) to conditions present during an atmospheric heatwave that impacted the Pacific Northwestern region of the United States in the summer of 2021. Climate stressors were applied either singly (single stressor; elevated seawater temperature, 30°C) or in succession (multiple stressor; elevated seawater temperature followed by aerial emersion at 44°C), replicating conditions present within the intertidal over a tidal cycle during the event. Oyster mortality rate was elevated within stress treatments with respect to the control and was significantly higher in triploids than diploids following multiple stress exposure (36.4% vs. 14.8%). Triploids within the multiple stressor treatment exhibited signs of energetic limitation, including metabolic depression, a significant reduction in ctenidium Na+ /K+ ATPase activity, and the dysregulated expression of genes associated with stress response, innate immunity, glucose metabolism, and mitochondrial function. Functional enrichment analysis of ploidy-specific gene sets identified that biological processes associated with metabolism, stress tolerance, and immune function were overrepresented within triploids across stress treatments. Our results suggest that triploidy impacts the transcriptional regulation of key processes that underly the stress response of Pacific oysters, resulting in downstream shifts in physiological tolerance limits that may increase susceptibility to extreme climate events that present multiple environmental stressors. The impact of chromosome set manipulation on the climate resilience of marine organisms has important implications for domestic food security within future climate scenarios, especially as triploidy induction becomes an increasingly popular tool to elicit reproductive control across a wide range of species used within marine aquaculture.

Parallel evolution in Crassostrea oysters along the latitudinal gradient is associated with variation in multiple genes involved in adipogenesis

Parallel diversification provides a proper framework for studying the role of natural selection in evolution. Yet, empirical studies from ecological 'non-model' species of invertebrates are limited at the whole genome level. Here, we presented a chromosome-scale genome assembly for Crassostrea angulata and investigated the parallel genomic evolution in oysters. Specifically, we used population genomics approaches to compare two southern-northern oyster species pairs (C. angulata-C. gigas and southern-northern C. ariakensis) along the coast of China. The estimated divergence time of C. angulata and C. gigas is earlier than that of southern and northern C. ariakensis, which aligns with the overall elevated genome-wide divergence. However, the southern-northern C. ariakensis FST profile represented more extremely divergent "islands". Combined with recent reciprocal hybridization studies, we proposed that they are currently at an early stage of speciation. These two southern-northern oyster species pairs exhibited significant repeatability in patterns of genome-wide differentiation, especially in genomic regions with extremely high and low divergence. This suggested that divergent and purifying selection has contributed to the genomic parallelism between southern and northern latitudes. Top differentiated genomic regions shared in these two oyster species pairs contained candidate genes enriched for functions in energy metabolism, especially adipogenesis, which are closely related to reproductive behaviours. These genes might be good candidates for further investigation in vivo. In conclusion, our results suggest that similar divergent selection and shared genomic features could predictably transform standing genetic variation within one species pair into differences in another.

Differentiation and decreased genetic diversity in field contaminated oysters Crassostrea hongkongensis: Identification of selection signatures

The extent to which chemical contamination affects the population structure and genetic diversity of natural populations remains elusive. Here, we used the whole-genome resequencing and transcriptome to diagnose the effects of long-term exposure to multiple elevated chemical pollutants on the population differentiation and genetic diversity in oysters Crassostrea hongkongensis in a typically polluted Pearl River Estuary (PRE) of Southern China. Population structure revealed an obvious differentiation between the PRE oysters and those collected from a nearby clean Beihai (BH) individuals, while no significant differentiation was observed among individuals collected from the three pollution sites within PRE due to the high gene flow. The decreased genetic diversity in the PRE oysters reflected the long-term effects of chemical pollutants. Selective sweeps between BH and PRE oysters revealed that chemical defensome genes, including glutathione S-transferase, zinc transporter, were responsible for their differentiation, sharing common metabolic process of other pollutants. Combined with the genome-wide association analysis, 25 regions containing 77 genes were identified to be responsible for the direct selection regions of metals. Linkage disequilibrium blocks and haplotypes within these regions provided the biomarkers of permanent effects. Our results provide important insights to the genetic mechanisms underlying the rapid evolution under chemical contamination in marine bivalves.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

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.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

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.

Predicting selection-response gradients of heat tolerance in a widespread reef-building coral

Ocean temperatures continue to rise owing to climate change, but it is unclear whether heat tolerance of marine organisms will keep pace with warming. Understanding how tolerance scales from individuals to species and quantifying adaptive potentials is essential to forecasting responses to warming. We reproductively crossed corals from a globally distributed species (Acropora tenuis) on the Great Barrier Reef (Australia) from three thermally distinct reefs to create 85 offspring lineages. Individuals were experimentally exposed to temperatures (27.5, 31 and 35.5°C) in adult and two critical early life stages (larval and settlement) to assess acquired heat tolerance via outcrossing of offspring phenotypes by comparing five physiological responses (photosynthetic yields, bleaching, necrosis, settlement and survival). Adaptive potentials and physiological reaction norms were calculated across three stages to integrate heat tolerance at different biological scales. Selective breeding improved larval survival to heat by 1.5-2.5× but did not result in substantial enhancement of settlement, although population crosses were significantly different. Under heat stress, adults were less variable compared with larval responses in warmer reefs than in the cooler reef. Adults and offspring also differed in their mean population responses, likely underpinned by heat stress imposing strong divergent selection on adults. These results have implications for downstream selection during reproduction, evidenced by variability in a conserved heat tolerance response across offspring lineages. These results inform our ability to forecast the impacts of climate change on wild populations of corals and will aid in developing novel conservation tools such as the assisted evolution of at-risk species.

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.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6–26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

Transcriptomics Analysis and Re-sequencing Reveal the Mechanism Underlying the Thermotolerance of an Artificial Selection Population of the Pacific Oyster

The Pacific oyster is a globally important aquaculture species inhabiting the intertidal environment, which experiences great temperature variation. Mass deaths in the summer pose a major challenge for the oyster industry. We initiated an artificial selection breeding program in 2017 using acute heat shock treatments of the parents to select for thermotolerance in oysters. In this study, we compared the respiration rate, summer survival rate, gene expression, and gene structure of F₂ selected oysters and non-selected wild oysters. A transcriptional analysis revealed global divergence between the selected and control groups at the larval stage, including 4764 differentially expressed genes, among which 79 genes were heat-responsive genes. Five heat shock proteins were enriched, and four of the six genes (five heat stock genes in the enriched GO terms and KEGG pathways and BAG4) were differentially expressed in 1-year-old oysters. Integration of the transcriptomic and re-sequencing data of the selected and the control groups revealed 1090 genes that differentiated in both gene structure and expression. Two SNPs (single nucleotide polymorphism) that may mediate the expression of CGI_10022585 and CGI_10024709 were validated. In addition, the respiration rate of 1-year-old oysters varied significantly between the selected group and the control group at room temperature (20°C). And the summer survival rate of the selected population was significantly improved. This study not only shows that artificial selection has a significant effect on the gene structure and expression of oysters, but it also helps reveal the mechanism underlying their tolerance of high temperature as well as the ability of oysters to adapt to climate change.

Transient Receptor Potential (TRP) Channels in the Pacific Oyster (Crassostrea gigas): Genome-Wide Identification and Expression Profiling after Heat Stress between C. gigas and C. angulata

Transmembrane proteins are involved in an array of stress responses, particularly in thermo-sensation and thermo-regulation. In this study, we performed a genome-wide identification and characterization of the Transient Receptor Potential (TRP) genes in the Pacific oyster (Crassostrea gigas) and investigated their expression profiles after heat stress to identify critical TRPs potentially associated with thermal regulation. A total of 66 TRP genes were identified in the C. gigas, which showed significant gene expansion and tandem duplication. Meta-analysis of the available RNA-Seq data generated from samples after acute heat stress revealed a set of heat-inducible TRPs. Further examination of their expression profiles under chronic heat stress, and comparison between C. gigas and C. angulata, two oyster species with different tolerance levels to heat stress, led to the identification of TRPC3.6, TRPC3.7, and TRPV4.7 as important TRPs involved in thermal regulation in oysters. This work provided valuable information for future studies on the molecular mechanism of TRP mediated thermal tolerance, and identification of diagnostic biomarker for thermal stress in the oysters.

Temperature adaptations of the thermophilic snail Echinolittorina malaccana: insights from metabolomic analysis

The periwinkle snail Echinolittorina malaccana, for which the upper lethal temperature is near 55°C, is one of the most heat-tolerant eukaryotes known. We conducted a multi-level investigation - including cardiac physiology, enzyme activity, and targeted and untargeted metabolomic analyses - that elucidated a spectrum of adaptations to extreme heat in this organism. All systems examined showed heat intensity-dependent responses. Under moderate heat stress (37-45°C), the snail depressed cardiac activity and entered a state of metabolic depression. The global metabolomic and enzymatic analyses revealed production of metabolites characteristic of oxygen-independent pathways of ATP generation (lactate and succinate) in the depressed metabolic state, which suggests that anaerobic metabolism was the main energy supply pathway under heat stress (37-52°C). The metabolomic analyses also revealed alterations in glycerophospholipid metabolism under extreme heat stress (52°C), which likely reflected adaptive changes to maintain membrane structure. Small-molecular-mass organic osmolytes (glycine betaine, choline and carnitine) showed complex changes in concentration that were consistent with a role of these protein-stabilizing solutes in protection of the proteome under heat stress. This thermophilic species can thus deploy a wide array of adaptive strategies to acclimatize to extremely high temperatures.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)—allopatric subspecies with distinct thermal tolerances—make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies’ responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

Temporal proteomic profiling reveals insight into critical developmental processes and temperature-influenced physiological response differences in a bivalve mollusc

Background

Protein expression patterns underlie physiological processes and phenotypic differences including those occurring during early development. The Pacific oyster (Crassostrea gigas) undergoes a major phenotypic change in early development from free-swimming larval form to sessile benthic dweller while proliferating in environments with broad temperature ranges. Despite the economic and ecological importance of the species, physiological processes occurring throughout metamorphosis and the impact of temperature on these processes have not yet been mapped out.

Results

Towards this, we comprehensively characterized protein abundance patterns for 7978 proteins throughout metamorphosis in the Pacific oyster at different temperature regimes. We used a multi-statistical approach including principal component analysis, ANOVA-simultaneous component analysis, and hierarchical clustering coupled with functional enrichment analysis to characterize these data. We identified distinct sets of proteins with time-dependent abundances generally not affected by temperature. Over 12 days, adhesion and calcification related proteins acutely decreased, organogenesis and extracellular matrix related proteins gradually decreased, proteins related to signaling showed sinusoidal abundance patterns, and proteins related to metabolic and growth processes gradually increased. Contrastingly, different sets of proteins showed temperature-dependent abundance patterns with proteins related to immune response showing lower abundance and catabolic pro-growth processes showing higher abundance in animals reared at 29 °C relative to 23 °C.

Conclusion

Although time was a stronger driver than temperature of metamorphic proteome changes, temperature-induced proteome differences led to pro-growth physiology corresponding to larger oyster size at 29 °C, and to altered specific metamorphic processes and possible pathogen presence at 23 °C. These findings offer high resolution insight into why oysters may experience high mortality rates during this life transition in both field and culture settings. The proteome resource generated by this study provides data-driven guidance for future work on developmental changes in molluscs. Furthermore, the analytical approach taken here provides a foundation for effective shotgun proteomic analyses across a variety of taxa.

Genetic Characterization of Cupped Oyster Resources in Europe Using Informative Single Nucleotide Polymorphism (SNP) Panels

The Pacific oyster, Crassostrea gigas, was voluntarily introduced from Japan and British Columbia into Europe in the early 1970s, mainly to replace the Portuguese oyster, Crassostrea angulata, in the French shellfish industry, following a severe disease outbreak. Since then, the two species have been in contact in southern Europe and, therefore, have the potential to exchange genes. Recent evolutionary genomic works have provided empirical evidence that C. gigas and C. angulata exhibit partial reproductive isolation. Although hybridization occurs in nature, the rate of interspecific gene flow varies across the genome, resulting in highly heterogeneous genome divergence. Taking this biological property into account is important to characterize genetic ancestry and population structure in oysters. Here, we identified a subset of ancestry-informative makers from the most differentiated regions of the genome using existing genomic resources. We developed two different panels in order to (i) easily differentiate C. gigas and C. angulata, and (ii) describe the genetic diversity and structure of the cupped oyster with a particular focus on French Atlantic populations. Our results confirm high genetic homogeneity among Pacific cupped oyster populations in France and reveal several cases of introgressions between Portuguese and Japanese oysters in France and Portugal.