Abstract Simple Summary With the acceleration of global industrialization, marine heavy metal pollution is becoming increasingly severe, posing significant threats to marine ecosystems and human health. To comprehensively understand the response mechanisms of mollusks to cadmium (Cd) stress, we selected the Pacific oyster (Crassostrea gigas) as the research species. We aim to delve into the molecular mechanisms underlying the response of oysters to cadmium exposure. Our results indicate that Cd exposure significantly increases Cd concentration in oyster gill tissues, eliciting strong physiological and metabolic reactions, including enhanced oxidative stress response and disrupted energy metabolism, signaling the oyster’s biological response to acute Cd stress. This study provides new insights into the signal transduction and metabolic pathways involved in the acute biological response of bivalves and other marine organisms to heavy metals. Abstract Marine mollusks, including oysters, are highly tolerant to high levels of cadmium (Cd), but the molecular mechanisms underlying their molecular response to acute Cd exposure remain unclear. In this study, the Pacific oyster Crassostrea gigas was used as a biological model, exposed to acute Cd stress for 96 h. Transcriptomic analyses of their gills were performed, and metabolomic analyses further validated these results. In our study, a total of 111 differentially expressed metabolites (DEMs) and 2108 differentially expressed genes (DEGs) were identified under acute Cd exposure. Further analyses revealed alterations in key genes and metabolic pathways associated with heavy metal stress response. Cd exposure triggered physiological and metabolic responses in oysters, including enhanced oxidative stress and disturbances in energy metabolism, and these changes revealed the biological response of oysters to acute Cd stress. Moreover, oysters could effectively enhance the tolerance and detoxification ability to acute Cd exposure through activating ABC transporters, enhancing glutathione metabolism and sulfur relay system in gill cells, and regulating energy metabolism. This study reveals the molecular mechanism of acute Cd stress in oysters and explores the molecular mechanism of high tolerance to Cd in oysters by using combined metabolomics and transcriptome analysis. Keywords: cadmium stress, Crassostrea gigas, metabolomic, transcriptomic, energy metabolism 1. Introduction Cadmium, a highly toxic heavy metal, tends to accumulate in the food chain. The growth of industries like mining and electroplating in recent years has led to an increase in the release of wastewater containing Cd. This ongoing discharge has resulted in the accumulation of cadmium elements in rivers and coastal waters [[50]1]. The Pacific oysters Crassostrea gigas primarily inhabit marine environments including estuaries, bays, and intertidal zones along the Pacific coast. Due to their high bioaccumulation capacity, they are capable of accumulating 104 times more Cd in their soft tissues compared to the concentration found in seawater [[51]2,[52]3]. Consequently, Cd pollution poses certain risks to human health and other predators’ food safety while also causing economic losses for the aquaculture industry [[53]4]. Currently, marine Cd pollution is recognized as a global issue with significant implications for marine organisms. Most bivalves are benthic filter feeders with a strong ability to accumulate heavy metals. Due to deficiencies in their metabolic mixed-oxidase systems, they release internal pollutants more slowly than fish and crustaceans, resulting in higher levels of accumulation in their bodies, which significantly affects growth, immunity, metabolism, and other processes [[54]5,[55]6]. Despite the elevated levels of heavy metals in the oysters’ bodies, they do not exhibit apparent tissue deformities, suggesting that they possess specialized mechanisms for accumulating, storing, and detoxifying heavy metals. Research on the heavy metal tolerance mechanisms of bivalves has been increasing annually, mainly focusing on aspects such as toxicology, immunology, antioxidant enzyme activity, functional genes, and DNA damage [[56]7,[57]8,[58]9,[59]10]. Marine species contaminated with heavy metals tend to accumulate substantial quantities of these metals, which markedly impact their growth, immunity, metabolism, and various other physiological processes. For instance, the Mediterranean mussel Mytilus galloprovincialis shows significant responses to Cd exposure by increasing metallothionein levels and antioxidant enzyme activity, which are biomarkers of metal exposure [[60]11]. In the clam Ruditapes philippinarum, Cd exposure increases the energy demand under cadmium stress by upregulating glycolysis and the TCA cycle pathways. [[61]12]. Additionally, the freshwater mussel Anodonta woodiana exhibits DNA damage and protein denaturation from cadmium exposure but mitigates these effects by activating repair genes and antioxidant defenses, demonstrating cadmium tolerance [[62]13]. Although different bivalves adopt various strategies to cope with cadmium pollution, they all exhibit adaptation and resistance mechanisms. However, there have been few studies on the response mechanisms of mollusks under acute Cd exposure, and none specifically on oysters. Metabolomics is a technique used to examine the qualitative and quantitative aspects of small molecule metabolites present in an organism after undergoing perturbations, such as alterations in genetics or surroundings. It aims to identify the types of metabolites and their patterns of change and to elucidate the relationships between metabolites and physiological changes [[63]14]. Transcriptomics refers to a methodology that studies gene transcription and transcriptional regulatory patterns at the global level within cells. It involves the investigation of gene expression at the RNA level [[64]15]. In bivalves, heavy metals tend to accumulate primarily in their gill tissues, which exhibit significant responses to metal toxicity [[65]16]. The gills of bivalve mollusks are primary target organs for metal accumulation due to their direct contact with heavy metals in water, often exhibiting the highest levels of pollutants [[66]17]. Studies have shown that among the gills, hepatopancreas, and mantle of Meretrix meretrix, the gills contain the highest levels of Cd. Furthermore, dissolved metals are primarily absorbed through direct adsorption onto the gills and mantle, while particulate metals tend to be absorbed via the digestive organs along with food [[67]18]. The gills also serve as a major barrier against damage and pathological factors mediated by environmental pollutants, making them the first site of toxic metal impact [[68]19]. Therefore, we chose to analyze the gill tissues of oysters for subsequent analysis. By employing a combined approach of transcriptomics and metabolomics, we can uncover the correlations between genes and metabolites, identify key regulatory nodes of genes and metabolites, and elucidate the regulatory mechanisms of metabolic pathways within organisms. Our study aims to provide further evidence of the toxicological mechanisms occurring in oysters under acute Cd exposure. In this study, we used the Pacific oyster C. gigas as a toxicological model to investigate the molecular mechanisms of molecular response and metal toxicity under acute Cd exposure, and we revealed the oxidative stress response induced by acute Cd exposure in the oyster. Through the enrichment analysis of transcriptome and metabolome data, we also studied the toxicity effects of acute Cd exposure on the oyster and the mechanisms regulating the strength of tolerance from the translational and metabolic levels. Through these analyses, our objectives were to (1) clarify the molecular mechanisms of oxidative stress and energy metabolism disorders in oysters under acute Cd exposure, and (2) explore the molecular mechanisms and potential molecular mechanisms of high Cd tolerance response in C. gigas. Our results have provided fundamental insights into the regulatory pathways of complex Cd responses. 2. Materials and Methods 2.1. Animal Materials and Cd Treatment C. gigas oysters were sourced from the Kongtong Island Breeding Base in Yantai, Shandong Province, China, and underwent a one-week adaptation period in a seawater tank. These one-year-old C. gigas had an average shell length of 69.07 mm (range: 61.92–76.23 mm) and an average weight of 45.71 g (range: 39.84–54.68 g) and were not in the breeding season. During both the adaptation phase and subsequent exposure to Cd, the oysters were kept in seawater filtered to maintain a temperature of 24 ± 2 °C. To minimize the introduction of metals through the food chain, they received a diet of freshly harvested microalgae (comprising both Diatom and Chlorella vulgaris) every alternate day. The oysters were allocated across six 60 L tanks, each containing seawater that had been filtered (three biological replicates for both the experimental and control groups), with 10 oysters per tank. A dilution of Cd was created by dissolving CdCl[2] in deionized water, achieving a concentration of 1 mg/L. The Cd concentrations were determined based on previous studies and our preliminary research. After a pre-experiment, the 96 h LC50 (semi-lethal concentration) was determined to be 65.75 mg/L. The stock solution was added to three tanks in the experimental group, with the Cd concentration in acute exposure set at one-twentieth of the LC50, resulting in a final seawater Cd concentration of 3.4 mg/L [[69]20]. Furthermore, concentrations higher than environmental standards were employed to provoke a significant stress response in the organism, thereby guaranteeing the validity of the results. In a study, the concentration of cadmium in the treatment group for greenfin horse-faced filefish (Thamnaconus septentrionalis) was set at 50 mg/L [[70]21]. Therefore, in research on the response of oysters to cadmium stress, a cadmium concentration of 3.4 mg/L in seawater is considered a reasonable level. The oysters were exposed to Cd for 96 h. Throughout the experiment, the seawater was changed daily, and the survival status of the oysters was monitored. Deceased individuals were promptly removed and recorded. For sample collection, gill tissues were sampled from three oysters in each of the six tanks (acute Cd exposure and control groups). The samples were promptly cryogenically preserved in liquid nitrogen and stored at a temperature of −80 °C. These samples were reserved for subsequent analyses, including transcriptomics, metabolomics, and the determination of physiological and biochemical indicators. 2.2. Cd Concentration and Biochemical Parameter Analysis Gill samples weighing 0.5 g were first homogenized, lyophilized, and then completely digested with a mixture of HCl (20%) and HNO[3] (80%). These samples were sealed and heated in an oven at 150 °C for a duration of 10 h. Upon cooling, the concentration of Cd in the oyster gill tissues was assessed using inductively coupled plasma mass spectrometry (ICP-MS) (Thermo Scientific iCAP Q ICP-MS, manufactured in Waltham, MA, USA). Additionally, 0.1 g of gill tissue from both groups (the Cd exposure group and the control group) were homogenized in a chilled phosphate buffer solution. The mixture was centrifuged at a speed of 3000× g for 15 min, resulting in the collection of the supernatant. Commercial assay kits were employed to quantify enzyme activity levels, including superoxide dismutase (SOD), catalase (CAT), malondialdehyde (MDA), and glutathione peroxidase (GSH-PX), within the collected supernatant [[71]22]. The overall quantification of soluble proteins was conducted using the Bradford technique, with bovine serum albumin (BSA) serving as the standard reference, and measurements were recorded at a wavelength of 595 nm. Enzyme assays were conducted in triplicate to ensure accuracy and reproducibility. 2.3. Metabolite Extraction, Detection, and Analysis Both control individuals and those exposed to acute Cd (9 samples/group) had their gill tissues cryopreserved in liquid nitrogen before being air-dried and pulverized into fine powder. This powder was then reconstituted in an 80% methanol solution and diluted with LC-MS-grade water to reach a final concentration of 53% methanol. After centrifugation (1500 g, 4 °C, 20 min), the resulting supernatant was analyzed using the LC-MS system. Metabolite annotation was carried out leveraging databases such as KEGG and LIPIDMaps. To analyze the data, we utilized MetaX 1.4.2 Software, which is specifically designed for interpreting metabolomics data and includes principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA). Statistical significance was determined using univariate analysis based on p-values. Metabolites showing variable importance in projection (VIP) > 1, p-value < 0.05, and fold change (FC) ≥ 2 or FC ≤ 0.5 were classified as DEMs. The roles and associated metabolic pathways of these metabolites were explored using the KEGG database. Pathway enrichment analysis of differential metabolites was conducted, deeming a pathway significantly enriched when the ratio x/n exceeded y/N, and the pathway’s p-value was <0.05. 2.4. Sampling, RNA Extraction, Library Construction, and RNA-Seq Analysis RNA was extracted from the gill tissues of nine oysters per group, with every three oysters pooled as one sample, for RNA sequencing of both the acute Cd exposure and control groups. Subsequent analyses utilized only high-quality, clean data. A reference genome index was built, and the clean, paired-end reads were mapped to this reference using Hisat2 v2.0.5 [[72]23]. The FPKM metric was used to quantify gene expression levels. To control the false discovery rate and adjust p-values, we employed the Benjamini and Hochberg method. In our study, a p-value ≤ 0.05 and |log[2] (fold change)| ≥ 1 were used as the thresholds to define DEGs. Enrichment analyses for the DEGs were executed using DAVID v6.8, leveraging the reference genome as background. The DEGs served as the validation set for GO and KEGG analyses, enabling the assessment of DEG distribution in various functions and pathways between the acute Cd-exposed and control groups. 2.5. Expression Validation Using Quantitative Real-Time PCR To validate the RNA-Seq findings, qRT-PCR was conducted on 12 selected DEGs ([73]Figure S1). Primers for these DEGs were designed using Primer Premier 6.0, and the 2^−ΔΔCt technique was employed to quantify their relative expression levels [[74]24]. 2.6. Integrative Analysis of Metabolome and Transcriptome DEGs (p-value < 0.05 and |log2FC| > 1) and DEMs (VIP > 1 and p < 0.05) were used for the integrative analysis of the acute Cd exposure and control groups. The Spearman method analyzed the correlation coefficients for integrating the metabolome and transcriptome data. A heat plot visually represented the relationships between genes and metabolites. 2.7. Statistical Analysis The comparison of Cd concentration and biochemical parameters between the acute Cd exposure and control groups was performed using a t-test, with a significance threshold set at p < 0.05. This analysis was conducted using SPSS 19.0 software. 3. Results 3.1. Cd Concentration and Biochemical Parameter Analysis [75]Figure 1a shows that the concentration of Cd in the gill tissues of oysters was 28-fold higher in the acute exposure group compared to the control group. [76]Figure 1b demonstrates that the survival rate of oysters after 96 h of acute Cd exposure was 83.3%. Our preliminary experiment revealed that adult oysters exhibit a remarkable tolerance to acute cadmium stress, surpassing that of other species. Among these, the activities of SOD, CAT, and GPx in the gill tissues of the Cd-treated group significantly increased ([77]Figure 1c,d,f), along with a significant elevation in MDA content ([78]Figure 1e). The results indicate a significant increase in these enzyme activity indices, preliminarily suggesting that acute Cd stress leads to an augmentation of oxidative stress and an enhanced molecular response within oyster gill cells. Figure 1. [79]Figure 1 [80]Open in a new tab Biochemical alterations of the acute Cd exposure C. gigas. (a) Increased intracellular Cd level in Cd-exposed oyster individuals. The bar chart depicts mean levels and standard deviation (SD) values, with n = 3. Significance differences were determined by t-test (p < 0.05), denoted by asterisks (*). (b) The survival rate of oysters after 96 h of acute Cd exposure. (c–f) After Cd exposure, oxidative stress markers such as SOD (c), CAT (d), MDA (e), and GPx (f) were measured. Data were presented as the mean ± SD (n = 3). Distinguishing letters were assigned to indicate significant differences (p < 0.05). Star plots depict the impact of Cd exposure on gill samples. 3.2. Transcriptomic Changes Induced by Acute Cd Exposure The volcano plot visually represents the distribution of DEGs between the groups treated with Cd and those under control conditions. Compared to the control group, the exposed oysters exhibited 1132 genes upregulated and 976 genes downregulated internally ([81]Figure 2a). A visual representation of the hierarchical clustering analysis conducted on all the DEGs revealed distinct segregation in response to Cd exposure ([82]Figure 2b). To investigate the relationship between the DEGs and cell function, we performed enrichment analysis of GO and KEGG pathways. As shown in [83]Figure 2c, significantly enriched GO terms include proteolysis, glutathione metabolic process; protein catabolic process; amino acid activation; organic substance catabolic process; tRNA metabolic process (ontology: biological process); proteasome complex; endopeptidase complex; proteasome core complex; peptidase complex (ontology: cellular component); catalytic activity, acting on a tRNA; ATP-dependent peptidase activity; cofactor binding; oxidoreductase activity; and metalloendopeptidase activity (ontology: molecular function). Figure 2. [84]Figure 2 [85]Open in a new tab Analysis of oyster transcriptome following a 96-hour exposure to Cd. (a) A volcano plot of DEGs is graphically represented with upregulated genes in red and downregulated genes in blue. (b) Hierarchical clustering based on the DEGs, where red signifies upregulation and blue signifies downregulation. (c) Enrichment of DEGs in GO terms categorized into cellular components, biological processes, and molecular functions. (d) KEGG pathway enrichment analysis of DEGs, with colors indicating p-value significance and bubble size reflecting the count of enriched genes. The enrichment of level 2 KEGG signaling pathways such as proteasome, protein processing in the endoplasmic reticulum, ABC transporters, sulfur relay system, amino sugar and nucleotide sugar metabolism, protein export, ECM–receptor interaction, glutathione metabolism, and the peroxisome pathway ([86]Figure 2d). However, this does not definitively establish that these physiological processes occurred within the organism. Therefore, we proceed to further examine the metabolome to ascertain whether metabolic pathways or metabolites relevant to the transcriptome results can be identified, thereby confirming the occurrence of physiological activities induced by acute Cd exposure within oyster gill cells. 3.3. Metabolomic Changes Induced by Acute Cd Exposure To examine metabolic alterations following acute Cd exposure, intracellular metabolites were scrutinized using an LC-MS system. Distinct metabolomic profiles emerged in the OPLS-DA model after Cd exposure ([87]Figure 3a,b). A total of 111 differentially expressed metabolites (DEMs) were identified, characterized by VIP values > 1.0 and p-values < 0.05. Among these, the Cd-exposed group exhibited elevated levels in 85 metabolites, whereas decreased levels were observed in 26 metabolites ([88]Figure 3d). The unsupervised hierarchical clustering analysis of the 110 differentially expressed metabolites further emphasized the distinct expression patterns between the two groups. ([89]Figure 3c). To delve deeper into the disrupted metabolic pathways, KEGG pathway analysis was conducted using MetaboAnalyst 5.0. The result showed the top 20 significantly enriched KEGG pathways ([90]Figure 4a). Of these, we identified metabolic pathways related to Cd exposure, including glutathione metabolism, arginine and proline metabolism, fatty acid biosynthesis, sulfur relay system, and ABC transporters, which were significantly altered ([91]Figure 4b). Figure 3. [92]Figure 3 [93]Open in a new tab Analysis of oyster metabolomes following a 96-hour exposure to Cd. (a) OPLS-DA score plot of the metabolomic data. (b) PLS-DA sorting test plot of the metabolomic data. (c) Hierarchical clustering was performed using 111 DEMs, with red indicating upregulation and blue representing downregulation, respectively. (d) Sample comparisons for the matchstick diagram. The top 20 metabolites of up and down are displayed in the matchstick diagram. The x-axis of the matchstick diagram represents log[2] (Fold Change) values, the y-axis represents metabolites, and the size of the points corresponds to VIP values. The metabolites that are upregulated and downregulated are represented by the red and blue points, respectively. Figure 4. [94]Figure 4 [95]Open in a new tab Pathway analysis of the DEMs. (a) Top 20 KEGG pathway results. p-values are represented by colors, while pathway impact is indicated by the size of the bubbles. (b) KEGG regulatory network diagram. Red circles represent individual metabolic pathways, yellow circles depict enzyme information related to specific substances, green circles indicate background substances for a metabolic pathway, purple circles represent information on molecular modules of a certain substance category, blue circles represent chemical interactions involving a specific substance, and green squares denote differentially expressed substances identified in this comparison. 3.4. Integrated Analysis of Metabolomic and Transcriptomic Data To conduct a comprehensive analysis of the impact of Cd on biochemical pathways using systems biology methods, we performed a joint-pathway analysis in MetaboAnalyst by integrating the DEMs and DEGs at the pathway level. [96]Figure 5 shows the results of the correlation analysis of the top 82 DEMs and the top 33 DEGs. The data shows that a total of 25 pathways were enriched, including arginine and proline metabolism, aminoacyl-tRNA biosynthesis, histidine metabolism, phenylalanine metabolism, beta-alanine metabolism, nitrogen metabolism, the biosynthesis of amino acids, thiamine metabolism, the sulfur relay system, porphyrin and chlorophyll metabolism, ABC transporters, glyoxylate and dicarboxylate metabolism, butanoate metabolism, purine metabolism, tyrosine metabolism, taurine and hypotaurine metabolism, alanine, aspartate, and glutamate metabolism, 2-oxocarboxylic acid metabolism, arginine biosynthesis, glutathione metabolism, cysteine and methionine metabolism, fatty acid biosynthesis, pantothenate and CoA biosynthesis, alpha-linolenic acid metabolism, and pyrimidine metabolism ([97]Figure 6a). Figure 5. [98]Figure 5 [99]Open in a new tab Correlation heatmap analysis of DEGs and DEMs. The DEMs are depicted on the y-axis, while the x-axis displays the DEGs. A correlation coefficient less than 0 is indicative of a negative correlation, and a coefficient greater than 0 indicates a positive correlation. Negative correlations are symbolized by the color blue, while positive correlations are represented by the color red. Figure 6. [100]Figure 6 [101]Open in a new tab Pathway enrichment analysis of DEGs and DEMs. (a) Bubble plot of KEGG enrichment for DEGs and DEMs. The x-axis represents the ratio of the number of enriched differential metabolites or genes annotated to metabolites or genes in the pathway to the total number in that pathway (ratio). The y-axis represents KEGG pathways jointly enriched in the metabolome and transcriptome. The count indicates the number of enriched metabolites or genes in the pathway. The colors represent p-values, with brighter colors indicating smaller p-values and more significant pathway enrichment. (b) iPath pathway map of shared enriched pathways. Colored boxes represent enriched pathways, nodes depict various biochemical molecules, lines represent biochemical reactions, and blue lines within the pathways indicate pathways jointly enriched with DEGs and DEMs. The construction and visualization of the gene-metabolism network were carried out utilizing Metscape. The differential genes and metabolites were mapped to the KEGG database for identification of common pathways, leading to the discovery of key biochemical and signal transduction pathways associated with these molecules. The results indicated energy metabolism, amino acid metabolism, carbohydrate metabolism, metabolism of terpenoids and polyketides, glutathione metabolism, glyoxylic acid and dicarboxylic acid metabolism, the sulfur relay system, and ABC transporters were significantly enriched ([102]Figure 6b). We further categorized these physiological processes in conjunction with the iPath diagram to elucidate the primary impact of Cd exposure. The bold blue lines in the iPath pathway map represent commonly enriched pathways identified by both omics analyses. Among them, the ABC transporter, glutathione metabolism, sulfur relay system, and energy metabolism were found to be significant. We selected these four pathways that exhibited the strongest association with Cd exposure for further investigation. Furthermore, we selected 12 core genes for qPCR validation. The results demonstrated that the changes in gene expression levels were consistent with those observed in the transcriptome analysis, ensuring the reliability of the findings ([103]Figure S1). By integrating DEMs and DEGs from these pathways, we attempted to identify the core genes and metabolites responsive to acute Cd exposure. These four pathways were significantly enriched in both the transcriptomic and metabolomic results. Based on our findings, we constructed a comprehensive toxicological mechanism map encompassing differential genes and metabolites involved in these pathways as well as their intricate interactions ([104]Figure 7). Figure 7. [105]Figure 7 [106]Open in a new tab Main biological pathway responses to acute Cd exposure in oysters. The blue boxes represent four core Cd exposure response pathways. The orange boxes represent DEMs associated with the main pathways. The green boxes represent downregulated DEGs related to the main pathways. The red elliptical frames represent upregulated DEGs associated with the main pathways. 4. Discussion 4.1. The Toxicological Response Mechanisms of Oysters under Cd Exposure Transcriptomics is a pivotal approach for investigating cellular functions and examining gene expression at the RNA level. Metabolomics unveils mechanisms underlying biological activities by elucidating patterns of changes in small-molecule metabolites. Through correlated analysis of transcriptomics and metabolomics, critical metabolic pathways and regulatory genes can be more accurately identified, thereby providing deeper insights into the response mechanisms within oysters under acute Cd stress [[107]25]. The oyster gill tissue exhibits a specific response to Cd exposure, acting as a primary barrier against damage and pathological factors [[108]26]. When gill cells face acute cadmium exposure, they undergo a series of physiological and biochemical changes. In our transcriptomic and metabolomic KEGG results, we observed significant changes in the following pathways. The enhancement of the ABC transporter pathway in the gill tissues indicates its critical role in cadmium efflux. Kim et al. demonstrated that ABC transporters are significantly activated under cadmium exposure, accelerating cadmium efflux and reducing intracellular toxic accumulation [[109]27]. Furthermore, we found significant changes in the glutathione metabolism within gill tissues. Under cadmium stress, gill tissues showed a marked increase in GSH levels, indicating a high capacity for synthesizing and utilizing GSH to counteract cadmium-induced oxidative stress, thus enhancing detoxification and protective mechanisms. This differential response highlights the unique role of gills in managing oxidative damage. Similarly, Espinoza et al. emphasized the importance of GSH in the detoxification of heavy metal Cd in the gill tissues of coho salmon Oncorhynchus kisutch [[110]28]. Additionally, the energy metabolism pathways in gill tissues also exhibited significant alterations. Our study detected notable changes in metabolic pathways related to ATP production and utilization, indicating that cadmium exposure disrupts the energy homeostasis in gill cells. This forces the cells to alter gene expression to support detoxification processes and damage repair mechanisms. Zhang et al.’s research supports this view, showing that gill cells in greenfin horse-faced filefish Thamnaconus septentrionalis adjust glycolysis, protein degradation, and fatty acid metabolism to meet the nutritional and energy demands under cadmium stress [[111]21,[112]29]. Our study uses an integrated approach combining transcriptomic and metabolomic analyses to gain a deeper and more comprehensive understanding of the acute toxicological mechanisms in Pacific oysters. It reveals the unique response mechanisms of oyster gill tissues to Cd exposure, including significant changes in ABC transporter proteins, glutathione metabolism, energy metabolism, and oxidative damage pathways. These findings not only confirm previous research results but also provide new insights into the molecular mechanisms underlying the acute stress responses of aquatic invertebrates to heavy metal Cd exposure. While it is well-known that oysters exhibit extreme tolerance to cadmium, our focus is to observe the significant molecular changes during acute cadmium stress. This approach aims to identify the most critical physiological changes and uncover the survival strategies oysters employ to cope with high cadmium concentrations. Understanding these immediate molecular responses is important for determining key biomarkers and physiological pathways involved in their remarkable tolerance. 4.1.1. The Role of ABC Transport Proteins in Cd Tolerance in Oysters During acute Cd stress, oyster gill cells exhibit an immediate response in both their cell membranes and organelles [[113]30]. Among these, ATP-binding cassette (ABC) transport proteins constitute a class of membrane proteins that play crucial roles in transport and efflux functions on the cell membrane [[114]31] ([115]Figure 7). The ABC transporters have a crucial function in cellular detoxification as they actively remove harmful substances and their byproducts from the cell through ATP-dependent binding and hydrolysis mechanisms [[116]32]. The ABC transport mechanism has been identified and characterized in various aquatic multicellular eukaryotes, including bivalve mussel Mytilus edulis [[117]33,[118]34], teleost fish [[119]35,[120]36,[121]37], and crustaceans [[122]38,[123]39]. Previous studies have indicated that the activation of ABC transport proteins in bacteria, plants, and animals is linked to protective mechanisms against the toxic metal cadmium. These proteins have shown their capability to alleviate the detrimental impacts of Cd exposure on individual development and viability [[124]40,[125]41,[126]42]. When the organism is subjected to cadmium stress, ABC transport proteins act as cadmium efflux pumps, enhancing Cd tolerance by reducing the intracellular Cd content [[127]43]. The activity and gene expression of ABC transport proteins have been identified in a range of aquatic organisms, serving to respond to changes in environmentally exogenous substances [[128]44,[129]45]. Studies employing gene knockout have demonstrated that ABC transporters play a pivotal role in Arabidopsis thaliana, yeast Saccharomyces cerevisiae, and marine ciliates Paramecium marinus, enhancing the organisms’ Cd efflux rate and tolerance [[130]46,[131]47,[132]48]. ABC transport proteins are considered key components of the multidrug resistance (MDR) system, capable of blocking the uptake of compounds by cells and controlling the efflux of metabolites. Abcb1 is a gene that synthesizes MDR-related ABC transport proteins. In aquatic invertebrates, exposure to 5 μg/L Cd upregulates the expression of ABCB1 in the gill tissue of the clam Ruditapes philippinarum, significantly activating this transporter and expediting the efflux of Cd^2+ [[133]49,[134]50]. In addition to the detoxification role of ABC transport proteins, acute Cd exposure is also associated with a significant molecular response in mollusks. For instance, the molecular response in oyster gill cells involves the upregulation of genes related to oxidative stress pathways. SOD is an antioxidative enzyme utilized for scavenging superoxide radicals, representing highly deleterious oxidative species. Elevated SOD activity signifies enhanced capacity within the oyster organism to combat oxidative stress [[135]51]. CAT primarily functions to degrade hydrogen peroxide, a toxic oxidative substance. Elevated CAT activity indicates an enhanced efficiency of oyster gill cells in clearing hydrogen peroxide [[136]52]. MDA serves as a lipid peroxidation byproduct, typically utilized as a marker for oxidative stress. Glutathione peroxidase (Gpx), with its principal role in degrading peroxides, including lipid peroxides, signifies an elevated defense against oxidative stress within oyster gill cells under acute Cd exposure when its activity increases [[137]53]. The elevation in MDA levels signifies an increased degree of membrane lipid oxidation in oyster gill cells under Cd treatment, indicating a certain level of oxidative damage to the cell membrane. Enzyme activities significantly increase, indicating an enhanced antioxidant defense mechanism. In this study, we observed a significant upregulation of genes such as ABCA1, ABCB1, and ABCC3, along with a notable downregulation of the ABCC4 gene under acute Cd exposure. Associated differential metabolites included glutamate, aspartate, branched-chain amino acids, neutral amino acids, and histidine. Our findings indicate that under acute Cd stress, genes associated with ABC transport proteins are expressed in oyster gill cells, leading to the generation of related metabolites. This suggests a crucial role of ABC transporters in alleviating Cd toxicity and maintaining the balance of cadmium concentrations within oyster gill cells. Moreover, the high expression of these ABC transporter-related genes and associated metabolites may contribute to the oyster’s high accumulation and tolerance to Cd, reflecting its adaptive response to Cd exposure. 4.1.2. Glutathione Metabolism Plays a Detoxifying Role Against Cd Upon the ingress of cadmium ions into the cellular environment, a series of physiological reactions is triggered. Oysters exhibit an exceptional tolerance and accumulation capacity for Cd, which indicates that Cd does not predominantly exist in a free state within the cells. It is inferred that oysters possess certain proteins capable of binding with Cd to form complexes, enabling Cd to remain in a stable state within the cell. This mechanism significantly mitigates the cytotoxic effects of Cd. Our findings reveal that the secondary pathway, glutathione metabolism, possesses this capability ([138]Figure 7). It executes vital antioxidative and detoxification duties within the cellular environment, especially in combating oxidative stress induced by heavy metals like Cd [[139]54]. Glutathione (GSH) represents a critically important internal antioxidant, with a plethora of essential physiological roles such as free radical elimination and detoxification, contributing significantly to the organism’s biochemical defense mechanism [[140]55]. Previous research has demonstrated the protective role of glutathione (GSH) in isolated fish gill cells against the cytotoxicity of metal ions. Cd can significantly increase the depletion of GSH in cells, with GSH acting as a direct chelator for intracellular metal ions. Metal ions, due to their strong affinity for sulfur, can bind to the thiol groups within GSH molecules. A portion of this remains inside the cell as stable metal chelate proteins, while another is expelled through transport proteins on the cell membrane, thus facilitating detoxification and reducing the loss inflicted by heavy metals on cells [[141]56,[142]57]. Furthermore, glutathione functions as an essential cofactor for glutathione peroxidase, facilitating the breakdown of hydrogen peroxide and lipid peroxides generated during redox processes [[143]58]. In this study, the increase in glutathione transferase levels may be attributed to the activation of the glutathione metabolic pathway. The most common oxidized form of glutathione in cells is glutathione disulfide (GSSG). Previous research has indicated that the GSH/GSSG redox buffer participates in many vital cellular functions, including cell differentiation, proliferation, apoptosis, ferroptosis, signal transduction, cytokine production, and molecular responses [[144]59]. The initial and pivotal step in glutathione biosynthesis involves the coupling of glutamate and cysteine, facilitated by glutamate–cysteine ligase. Subsequently, the produced γ-glutamylcysteine combines with glycine, culminating in the formation of GSH by the enzyme GSH synthetase [[145]60]. In this research, the expression of genes related to the glutathione metabolic pathway, such as PEPN, LAP3, and MGST3, was significantly increased, whereas the expression of TRYS was significantly decreased. Differential metabolites include Glutathione disulfide (GSSG), Glutathione, R-S-Cysteinyl-glycine, Mercapturic acid, and L-Ornithine. These genes may influence the pathways associated with glutathione metabolism, thus enhancing the survival rate and tolerance of oysters under acute Cd stress and leading to the generation of related metabolites. The sulfur relay system constitutes the third pathway, functioning in the thiolation of sulfur within tRNA and the pathway for cysteine production. A tight interplay exists between this system and glutathione metabolism, collaboratively ensuring the equilibrium of sulfur metabolism and antioxidative protection within cells. The sources of sulfur for the sulfur relay system comprise an array of enzymes and molecules connected to sulfur, like cysteine, methionine, and cystathionine, among others [[146]61]. These molecules are precursor substances for the synthesis of glutathione. Among them, cysteine is a critical component in the synthesis of glutathione because it contains sulfur, serving as the sulfur source for glutathione molecules. In European seabass (Dicentrarchus labrax), alterations in the gut microbiota impact immune status and digestive performance. The organism can adapt to these changes through the regulation of the sulfur relay system [[147]62]. Glutathione is synthesized from cysteine, cystathionine, and methionine, with both cysteine and cystathionine being directly related to the sulfur relay system. These precursor substances undergo a series of enzymatic reactions to ultimately produce glutathione. The sulfur relay system, by providing a source of sulfur, maintains the sulfur metabolic balance within oysters, sustains glutathione levels, and enhances antioxidative defense [[148]63,[149]64]. In this study, genes related to the sulfur relay system, such as NCS2, NCS8, MOAE, and MOAC, were significantly upregulated. Differential metabolites include GTP, SAM, and MOCO. This indicates that in oysters under acute Cd exposure, pathways related to the sulfur relay system are activated and work in conjunction with the glutathione metabolic pathway to maintain intracellular antioxidative capacity, playing a critical role in cellular homeostasis. Our results underscore the critical role of the glutathione metabolic pathway in this process, as it exerts important antioxidative and detoxifying effects within the cellular environment. It serves as a significant internal antioxidant, playing a central role in combating oxidative stress induced by heavy metals such as cadmium. The upregulation of genes associated with glutathione metabolism and the activation of the sulfur relay system enhance the survival and tolerance of oysters under acute cadmium exposure, thereby strengthening the molecular response of oysters to cadmium-induced oxidative damage. This interconnected network of pathways emphasizes the complex biological strategies adopted by oysters to counter cadmium stress, reflecting their ability to enhance immunity and maintain cellular homeostasis in the face of heavy metal toxicity. Understanding these interrelated mechanisms provides fundamental information for studying defense mechanisms and cadmium tolerance in mollusks. 4.1.3. Acute Cadmium Exposure Affects Gill Energy Metabolism The maintenance of physiological functions and the regulation of stress heavily rely on energy metabolism. Typically, organisms are adequately supplied with energy to meet their requirements, while any excess is stored [[150]65]. However, the energy balance of marine organisms can be disturbed by heavy metal stress. As a result, organisms may need to expend additional energy in order to activate stress responses, facilitate detoxification processes, and engage in damage repair mechanisms [[151]66]. In energy metabolism, glycolysis is a crucial pathway for glucose metabolism. Glycolysis is the process of breaking down a glucose molecule into two pyruvate molecules, resulting in a net production of two ATP molecules and two NADH molecules. In the presence of oxygen, pyruvate enters the mitochondria, where it undergoes further metabolism via the TCA cycle and oxidative phosphorylation, leading to a significant generation of ATP [[152]67]. We discovered that in the Cd-treated groups, the levels of metabolites involved in glycolysis significantly increased, such as dodecanoic acid, acetyl-CoA, and ACP. In addition, the downregulation of the FASN gene caught our attention, indicating that the energy metabolism process was promoted under Cd exposure ([153]Figure 7). Previous research has found that with the increase in Cd concentration in mice, the glycolysis pathway was significantly enhanced [[154]68]. Furthermore, Cd can trigger oxidative stress reactions within cells, resulting in an overproduction of oxygen free radicals. These radicals can harm biomacromolecules within cells, such as proteins, lipids, and DNA, consequently disrupting normal cellular metabolism [[155]69,[156]70,[157]71,[158]72]. In our results, Cd exposure induced changes in fatty acid metabolism, alanine, and arginine metabolism. Alanine is an energy substrate primarily involved in the tricarboxylic acid cycle, while arginine primarily maintains stable ADP concentrations and can act as an energy buffer in tissues requiring high energy. In the study of rainbow trout Oncorhynchus mykiss, it was observed that the production of fatty acids and arginine under Cd exposure affected cellular homeostasis and metabolism [[159]73,[160]74]. Our study also indicates that the inhibitory effect of Cd exposure may lead to disruptions in the energy metabolism pathways of oyster organisms and oxidative stress damage. The tricarboxylic acid cycle is a sequence of biochemical reactions employed by aerobic organisms to produce energy via the oxidation of acetyl-CoA obtained from carbohydrates, fatty acids, and proteins [[161]75]. Our results reveal that metabolomic analysis showed an accumulation of free fatty acids and some small peptide molecules, such as dodecanoic acid, acetyl-CoA, and ACP. It has been reported that Cd exposure alters the abundance of metabolites related to energy metabolism (ATP, ACP, phosphatidylcholine, lactate, and succinate) in the liver of the halibut Paralichthys olivaceus [[162]76]. In our results, we also detected relevant metabolites among which the increased levels of ACP may result in cellular damage. ACP is a lysosomal enzyme with activity in multiple biological processes, including cathepsin, pathological necrosis, autophagy, and phagocytosis. Metals can accumulate in lysosomes, leading to increased toxicity and alterations in lysosomal function or structure, consequently causing cellular damage [[163]77,[164]78]. In the transcriptome results, we found that the FASN gene was significantly downregulated. Previous studies have shown that the downregulation of FASN can alleviate energy metabolism disorders induced by potassium dichromate in LMH cells. This indicates that the TCA cycle pathway is activated and enhanced under Cd stress, leading to energy disruption and structural–functional changes in oyster gill cells. In summary, Cd exposure significantly activates the glycolytic pathway, increasing the levels of metabolites involved in glycolysis. Cd stress triggers oxidative stress reactions within cells, leading to the generation of oxygen radicals that damage cellular biomolecules, resulting in alterations in fatty acid and amino acid metabolism, ultimately causing disruptions in energy metabolism and cellular dysfunction. These findings reveal the potential impact of heavy metal pollution on the energy balance and physiological functions of marine organisms. 5. Conclusions This study, through a comprehensive analysis of the metabolomic and transcriptomic responses of Pacific oysters C. gigas to acute Cd exposure, delves into the molecular mechanisms of oysters’ response to heavy metal Cd pollution, especially in marine environments. Our results indicate that Cd exposure significantly increases Cd concentration in oyster gill tissues, eliciting strong physiological and metabolic reactions, including enhanced oxidative stress response and disrupted energy metabolism, signaling the oyster’s biological response to acute Cd stress. Furthermore, oysters can enhance their tolerance and detoxification capabilities under acute Cd exposure by activating key pathways such as ABC transporters, regulating glutathione metabolism, participating in the sulfur relay system, and adjusting energy metabolism, thereby bolstering their molecular response to Cd-induced cellular damage. Our study not only offers new perspectives on the acute biological responses of marine organisms like bivalves to heavy metals but also provides important scientific evidence for the study of the biological effects and mechanisms of heavy metal pollution in marine ecosystems. Supplementary Materials The following supporting information can be downloaded at: [165]https://www.mdpi.com/article/10.3390/ani14162318/s1. Figure S1. Comparison of qPCR and RNA-Seq results of 12 hub genes. Expression levels of key genes were standardized using the expression levels of β-actin genes. The x-axis represents the Cd-exposed group and the control group, while the y-axis represents the fold change in gene expression comparing the Cd-exposed group to the control group. The letters indicate the level of significant difference between groups (p < 0.05). [166]animals-14-02318-s001.zip^ (203.2KB, zip) Author Contributions Conceptualization, L.D. and Y.S.; Methodology, L.D., Y.S. and P.W.; Software, L.D., Y.S. and X.X.; Validation, Y.S. and Z.L.; Formal analysis, L.D. and Y.S.; Investigation, L.D., Y.S., B.L. and Z.L.; Resources, Y.F.; Data curation, M.C. and Y.X.; Writing—original draft, L.D.; Writing—review & editing, C.C. and W.W.; Visualization, B.L., C.C. and W.W.; Supervision, Y.F., G.S. and J.Y.; Funding acquisition, Z.W. and J.Y. All authors have read and agreed to the published version of the manuscript. Institutional Review Board Statement Not applicable. Informed Consent Statement Not applicable. Data Availability Statement The raw data supporting the conclusions of this article will be made available by the authors upon request. Conflicts of Interest Bin Li and Zhongping Wang was employed by the Yantai Kongtong Island Industrial Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Funding Statement This work was supported by the Earmarked Fund for Agriculture Seed Improvement Project of Shandong Province, China (No. 2020LZGC016), and the Fund of Central Government Guiding Local for Scientific and Technological Development (No. YDZX2021017). Footnotes Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. References