Molecular and cellular effects of temperature in oysters Crassostrea brasiliana exposed to phenanthrene

Emerging and legacy contaminants on the Brazilian southern coast (Santa Catarina): A multi-biomarker approach in oysters Crassostrea gasar (Adanson, 1757)

Coastal environments, such as those in the Santa Catarina State (SC, Brazil), are considered the primary receptors of anthropogenic pollutants. In this study, our objective was to evaluate the levels of emerging contaminants (ECs) and persistent organic pollutants (POPs) in indigenous Crassostrea gasar oysters from different regions of SC coast in the summer season (March 2022). Field collections were conducted in the São Francisco do Sul, Itajaí, Florianópolis and Laguna coastal zones. We analyzed the bioaccumulation levels of 75 compounds, including antibiotics (AB), endocrine disruptors (ED), non-steroidal anti-inflammatory drugs (NSAIDs), polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs) and pesticides. Furthermore, we assessed biomarker responses related to biotransformation, antioxidant defense, heat shock protection and oxidative damage in oysters' gills. Prevalence of ECs was observed in the central and southern regions, while the highest concentrations of POPs were detected in the central-northern regions of SC. Oysters exhibited an induction in biotransformation systems (cyp2au1 and cyp356a1, sult and GST activity) and antioxidant enzymes activities (SOD, CAT and GPx). Higher susceptibility to lipid peroxidation was observed in the animals from Florianópolis compared to other regions. Correlation analyses indicated possible associations between contaminants and environmental variables in the biomarker responses, serving as a warning related to climate change. Our results highlight the influence of anthropogenic activities on SC, serving as baseline of ECs and POPs levels in the coastal areas of Santa Catarina, indicating more critical zones for extensive monitoring, aiming to conserve coastal regions.

Cefotaxime Exposure-Caused Oxidative Stress, Intestinal Damage and Gut Microbial Disruption in Artemia sinica

Cefotaxime (CTX) is an easily detectable antibiotic pollutant in the water environment, but little is known about its toxic effects on aquatic invertebrates, especially on the intestine. Here, we determined the oxidative stress conditions of A. sinica under CTX exposure with five concentrations (0, 0.001, 0.01, 0.1 and 1 mg/L) for 14 days. After that, we focused on changes in intestinal tissue morphology and gut microbiota in A. sinica caused by CTX exposure at 0.01 mg/L. We found malondialdehyde (MDA) was elevated in CTX treatment groups, suggesting the obvious antibiotic-induced oxidative stress. We also found CTX exposure at 0.01 mg/L decreased the villus height and muscularis thickness in gut tissue. The 16S rRNA gene analysis indicated that CTX exposure reshaped the gut microbiota diversity and community composition. Proteobacteria, Actinobacteriota and Bacteroidota were the most widely represented phyla in A. sinica gut. The exposure to CTX led to the absence of Verrucomicrobia in dominant phyla and an increase in Bacteroidota abundance. At the genus level, eleven genera with an abundance greater than 0.1% exhibited statistically significant differences among groups. Furthermore, changes in gut microbiota composition were accompanied by modifications in gut microbiota functions, with an up-regulation in amino acid and drug metabolism functions and a down-regulation in xenobiotic biodegradation and lipid metabolism-related functions under CTX exposure. Overall, our study enhances our understanding of the intestinal damage and microbiota disorder caused by the cefotaxime pollutant in aquatic invertebrates, which would provide guidance for healthy aquaculture.

Effect of elevated temperature, sea water acidification, and phenanthrene on the expression of genes involved in the shell and pearl formation of economic pearl oyster (Pinctada radiata)

Our study aims to examine the effect of some stressors on the gene expression levels of shell matrix proteins in a pearl oyster. Oysters were exposed to the different combinations of the temperature, pH, and phenanthrene concentration is currently measured in the Persian Gulf and the predicted ocean warming and acidification for 28 days. The expression of all the studied genes was significantly downregulated. Time and temperature had the greatest effects on the decreases in n19 and n16 genes expression, respectively. Aspein and msi60 genes expression were highly influenced by pH. Pearlin was affected by double interaction temperature and phenanthrene. Moreover, a correlation was observed among the expression levels of studied genes. This study represents basic data on the relationship between mRNA transcription genes involved in the shell and pearl formation and climate changes in pollutant presence conditions and acclimatizing mechanism of the oyster to the future scenario as well.

Gender influences molecular and histological biomarkers in mature oysters Crassostrea gasar (Adanson, 1757) after pyrene exposure

Oysters are frequently used as sentinel organisms for monitoring effects of contaminants due to their sessile, filtering habits and bioaccumulation capacity. These animals can show elevated body burden of contaminants, such as pyrene (PYR). PYR can be toxic at a molecular level until the whole oyster, which can show reproductive and behavioral changes. Considering that biologic parameters, such as gender or reproductive stage can interfere in the toxic effects elicited by contaminants uptake, the aim of this study was to evaluate some molecular and histological responses in females and males of oyster Crassostrea gasar exposed to PYR (0.25 and 0.5 μM) for 24 h at the pre-spawning stage. PYR concentrations were analyzed in water and in tissues of female and male oysters. Gene transcripts related to biotransformation (CYP3475C, CYP2-like, CYP2AU1, CYP356A, GSTO-like, GSTM-like, SULT-like), stress (HSP70), and reproduction (Vitellogenin, Glycoprotein) were quantified in gills. In addition, histological analysis and histo-localization of CYP2AU1 mRNA transcripts in gills, mantle and digestive diverticulum were carried out. Females and males in pre-spawning stage bioconcentrated PYR in their tissues. Males were more sensitive to PYR exposure. CYP2AU1 transcripts were higher in males (p < 0.05), as well as tubular atrophy was observed only in males exposed to PYR (p < 0.05). As expected, vitellogenin transcripts were lower in males (p < 0.05). Given these results, it is suggested that levels of CYP2AU1 be a good biomarker of exposure to PYR in oyster C. gasar and that it is important to consider the gender for the interpretation of biomarker responses.

Molecular changes in oysters Crassostrea gigas (Thunberg, 1793) from aquaculture areas of Santa Catarina Island bays (Florianópolis, Brazil) reveal anthropogenic effects

Anthropogenic activities in coastal regions cause risks to the environmental and human health. Due to the carcinogenic and mutagenic potential, polycyclic aromatic hydrocarbons (PAH) are considered priority for monitoring. Most of the Brazilian production of Crassostrea gigas oysters are placed in the Bays of Santa Catarina Island. The aim of this study was to evaluate molecular responses (phase I and II of biotransformation and antioxidant defense) of C. gigas from six oyster farming areas potentially contaminated by sanitary sewage in Florianópolis Metropolitan (SC, Brazil): Santo Antônio de Lisboa, Sambaqui, Serraria, Caieira, Tapera, Imaruim. We evaluated the transcript levels of CYP1A1-like, CYP2-like, CYP2AU2-like, CYP356A1, GSTA1A-like, GSTO.4A-like, SULT-like, SOD-like and CAT-like by qRT-PCR. Only oysters from Caieira showed levels of thermotolerant coliforms allowed by the law. Chemicals analyses in soft tissues of oysters showed low to average levels of PAH in all monitored areas. Enhanced transcript levels of phase I (CYP1A1-like, CYP3564A1-like, CYP2-like and CYP2AU2-like) were observed in oysters from Serraria and Imaruí, suggesting higher biotransformation activity in these farming areas. Regarding phase II of biotransformation, GSTO.4A-like was up-regulated in oysters from Imaruí compared to Caieira and Santo Antônio de Lisboa. An upregulation of SOD-like and CAT-like were observed in oysters from Imaruí and Serraria, suggesting that oysters from these sites are facing higher prooxidant conditions compared to other areas. By integrating the biological and chemical data it is suggested that human-derived contaminants are affecting the oyster metabolism in some farming areas.

Analysis of Crassostrea gasar transcriptome reveals candidate genes involved in metal metabolism

Oysters have been extensively employed for monitoring of metal pollution in dynamic aquatic ecosystems. Therefore, the use of specific biomarkers can assist in discriminating the ecotoxicological implications of different elements in such complex environments. In this study, we revisited the sequencing data of gills and digestive glands transcripts in the mangrove oyster Crassostrea gasar and generated a reference transcriptome assembly from multiple assemblers, seven in total. Overall, we were able to identify a total of 11,917 transcripts, with 86.6% of them being functionally annotated and 1.4 times more than the first annotation. We screened the annotated transcripts to identify genes potentially involved in metals' transport, storage, and detoxification. Our findings included genes related to Zn distribution in cells (Zn transporters - ZIP, ZnT), metallothionein (MT-I and MT-IV), GSH biosynthesis, Ca⁺ transporter (NCX and ATP2B), and Cu distribution in cells (ATP7, ATOX1, CCS, and laccase-like). These results provided a reference transcriptome for additional insights into the transcriptional profile of C. gasar and other bivalves to better understand the molecular pathways underpinning metal tolerance and susceptibility. The study also provided an auxiliary tool for biomonitoring metal contamination in dynamic environments as estuaries.

Moderate acidification mitigates the toxic effects of phenanthrene on the mitten crab Eriocheir sinensis

Freshwater acidification and phenanthrene may result in complex adverse effects on aquatic animals. Juvenile Chinese mitten crabs (Eriocheir sinensis) were exposed to different pH levels (7.8, 6.5, and 5.5) under phenanthrene (PHE) (0 (control) and 50 μg/L) conditions for 14 days. Antioxidant and transcriptomic responses were determined under stress conditions to evaluate the physiological adaptation of crabs. Under the control pH 7.8, PHE led to significantly reduced activities of superoxide dismutase (SOD), catalase (CAT), glutathione reductase (GR) and glutathione S-transferase (GST), but increased glutathione peroxidase (GSH-Px), 7-ethoxyresorufin-o-deethylase (EROD) activities, and malondialdehyde (MDA) levels. However, moderate acidification (pH 6.5) changed PHE effects by increasing antioxidant enzymes. Acidification generally reduced SOD, GPx, GST and EROD activities, but increased CAT, GR, MDA. Compared with pH7.8 group, pH7.8 × PHE and pH6.5 × PHE groups had 1148 and 1498 differentially expressed genes, respectively, with "Biological process" being the main category in the two experimental groups. pH7.8 × PHE treatment caused significant enrichment of disease and immune-related pathways, while under pH6.5 × PHE, more pathways related to metabolism, detoxification, environmental information processing, and energy supply were significantly enriched. Thus, PHE had a significant inhibitory effect on antioxidant performance in crabs, while moderate acidification (pH6.5) mitigated the toxic effects of PHE. Overall, moderate acidification has a positive effect on the defense against the negative effects of PHE in Chinese mitten crabs, and this study provides insights into the defense mechanism of crustaceans in response to combined stress of acidification and PHE.

Impact of Polycyclic Aromatic Hydrocarbon Accumulation on Oyster Health

In the past decade, the Deepwater Horizon oil spill triggered a spike in investigatory effort on the effects of crude oil chemicals, most notably polycyclic aromatic hydrocarbons (PAHs), on marine organisms and ecosystems. Oysters, susceptible to both waterborne and sediment-bound contaminants due to their filter-feeding and sessile nature, have become of great interest among scientists as both a bioindicator and model organism for research on environmental stressors. It has been shown in many parts of the world that PAHs readily bioaccumulate in the soft tissues of oysters. Subsequent experiments have highlighted the negative effects associated with exposure to PAHs including the upregulation of antioxidant and detoxifying gene transcripts and enzyme activities such as Superoxide dismutase, Cytochrome P450 enzymes, and Glutathione S-transferase, reduction in DNA integrity, increased infection prevalence, and reduced and abnormal larval growth. Much of these effects could be attributed to either oxidative damage, or a reallocation of energy away from critical biological processes such as reproduction and calcification toward health maintenance. Additional abiotic stressors including increased temperature, reduced salinity, and reduced pH may change how the oyster responds to environmental contaminants and may compound the negative effects of PAH exposure. The negative effects of acidification and longer-term salinity changes appear to add onto that of PAH toxicity, while shorter-term salinity changes may induce mechanisms that reduce PAH exposure. Elevated temperatures, on the other hand, cause such large physiological effects on their own that additional PAH exposure either fails to cause any significant effects or that the effects have little discernable pattern. In this review, the oyster is recognized as a model organism for the study of negative anthropogenic impacts on the environment, and the effects of various environmental stressors on the oyster model are compared, while synergistic effects of these stressors to PAH exposure are considered. Lastly, the understudied effects of PAH photo-toxicity on oysters reveals drastic increases to the toxicity of PAHs via photooxidation and the formation of quinones. The consequences of the interaction between local and global environmental stressors thus provide a glimpse into the differential response to anthropogenic impacts across regions of the world.

Common and particular biochemical responses of Unio tumidus to herbicide, pharmaceuticals and their combined exposure with heating

The priority list of freshwater pollutants is increasingly amended by pharmaceuticals. Their impact on the aquatic biota can be modulated by the presence of typical pollutants, like pesticides, and/or abnormal heating. The aim of this study was to elucidate potentially hazardous impact of combined environmental factors on the freshwater mussels by analyzing various sets of biochemical markers. We treated the bivalve molluscs of Unio tumidus with non-steroidal anti-inflammatory drug diclofenac (Dc, 2 nM), calcium antagonist and antihypertensive drug nifedipine (Nf, 2 nM) or organophosphonate glyphosate-based herbicide Roundup MAX (Rn, 79 nM of glyphosate) at 18 °C as well as with the mixture of these substances at 18 °C (Mix) or 25 °C (MixT) during 14 days. The concentrations used were correspondent to the environmentally relevant levels. The biomarkers of stress and toxicity were evaluated in digestive gland, except the lysosomal membrane stability measured in hemocytes. Exposures caused an oxidative stress due to the decreased SOD and GST activities and GSH/GSSG ratio, increased levels of thiobarbituric acid-reactive substances and protein carbonyls (with some exceptions). Dc increased cathepsin D activity in lysosomes. Nf increased lysosomal membrane stability and caspase-3 activity. Rn caused a dramatic distortion of metallo-thiolome due to increased levels of GSH and metallothionein-related thiols (MTSH) as well as depletion of Zn, Cu and Cd in the composition of metallothioneins, and decreased Zn/Cu molar ratio in the tissue. The particular toxicity of Rn was also attested by decreased lysosomal membrane stability and cholinesterase activity. Canonical discriminant analysis separated Rn-, Mix- and MixT-groups from the joint set of C-, Dc- and Nf-groups. Generally, compound-specific effects were expressed in U. tumidus responses to the mixtures, but in MixT-group some effects were particular or extremely strong. Multi-marker approach and integrative analysis proved to be a useful tool for understanding possible future risks to freshwater mussels under a combination of xenobiotics and warming climate.

Biotransformation and oxidative stress responses in relation to tissue contaminant burden in Clarias gariepinus exposed to simulated leachate from a solid waste dumpsite in Calabar, Nigeria

In this study, we have investigated biotransformation and oxidative stress responses in relation to tissue contaminant burden in the African sharptooth catfish (Clarias gariepinus) exposed to simulated leachate from a solid waste dumpsite in Calabar, Nigeria. Fish were exposed to simulated leachate, diluted to 0:0 (negative control), 1:10, 1:50, 1:100 and phenanthrene (a PAH: 50 μg/L used as a positive control) for 3, 7 and 14 days. Hepatic transcripts for cat, sod1, gpx1, gr, gst, cyp1a, cyp2d3, and cyp27 were analyzed by real-time PCR, while enzymatic assays for ethoxyresorufin O-deethylase (EROD), buthoxyresorufin O-deethylase (BROD), methoxyresorufin O-deethylase (MROD), pentoxyresorufin O-deethylase (PROD), catalase (CAT), glutathione peroxidase (GPx), glutathione reductase (GR), glutathione S-transferase (GST), uridine diphospho-glucuronosyltransferase (UDPGT) and lipid peroxidase (LPO) were measured using standard methods. In addition, protein expression for CYP1A, CYP3A and metallotheionin (MT) were measured by immunoblotting. Fish muscle samples were analyzed for selected group of contaminants after 14 days exposure showing significantly high uptake of heavy metals (Cd, Hg and Pb), polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), chlorophenols, organochlorine (OC) and organophosphate pesticides in exposed fish. We observed significant concentration- and time-specific increases in biotransformation and oxidative stress responses at transcript and functional (enzyme and protein) levels, that paralleled tissue contaminants bioaccumulation patterns, after exposure to the simulated leachates. Our results highlighted the potential environmental, wildlife and public health consequences from improper solid waste disposal. In addition, it also provides a scientific basis for local sensitization and inform legislative decisions and policy formulation towards sustainable environmental management of solid wastes in Nigeria and other developing countries.

Short-term spatiotemporal biomarker changes in oysters transplanted to an anthropized estuary in Southern Brazil

Estuarine ecosystems are increasingly being affected by pollution caused by anthropogenic activities. In this study, Crassostrea gasar oysters were transplanted and maintained for seven days at three sites (S1, S2, and S3) in the Laguna Estuarine System (LES)-situated in southern Brazil-that has been exposed to multiple anthropic stresses. On the basis of the concentrations of metal and organic pollutants in oysters, we identified marked spatial variations in pollutant levels, with S3 showing the highest concentration of Ag, Fe, Ni, Zn, and total polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), and linear alkylbenzenes (LABs), followed by S2 and S1. Along with the concentrations of pollutants, a set of biomarkers was analyzed. Oysters maintained at S3 showed enhanced protective defenses in gills, as observed by the increased levels of superoxide dismutase (SOD-like) and heat shock protein 90 (HSP90-like) transcripts and catalase (CAT) activity, concomitant with reduced lipid peroxidation (MDA) levels. Decreased antioxidant activities together with increased MDA levels are indicative of the digestive gland being more susceptible to pollutant-induced oxidative damage. Oysters transplanted into LES showed lower levels of cytochrome P450 transcripts (CYP356A1-like and CYP2AU1), and decreased glutathione S-transferase (GST) enzyme activity, suggesting lower biotransformation capacity. By integrating information regarding the concentration of metal and organic pollutants with that of molecular as well as biochemical biomarkers, our study provides novel insights into pollutant exposure and the potential biological impacts of such exposure on estuarine organisms in southern Brazil.

Decreased malondialdehyde levels in fish (Astyanax altiparanae) exposed to diesel: Evidence of metabolism by aldehyde dehydrogenase in the liver and excretion in water

Increased malondialdehyde (MDA) levels are commonly considered an indicator of lipid peroxidation derived from oxidative stress insults promoted by exposure of fish to pollutants. However, a decrease in MDA levels after xenobiotic exposure has been also reported, an effect that is mostly attributed to enhanced antioxidant defenses. In this study, we assessed whether pollutant-mediated MDA decrease would be associated with antioxidant enhancement or with its metabolism by aldehyde dehydrogenase (ALDH) in the liver and gills of lambari (Astyanax altiparanae) exposed to diesel oil (0.001, 0.01, and 0.1 mL/L). MDA levels were decreased in the liver of lambari exposed to diesel. The activities of the antioxidant enzymes, catalase (CAT) and glutathione peroxidase (GPx), were unchanged in the liver, while that of glucose-6-phosphate dehydrogenase (G6PDH) was decreased. In contrast, levels of total glutathione (tGSH) and the activity of glutathione S-transferase (GST) were increased in the liver, which partly support antioxidant protection against lipid peroxidation. More importantly, ALDH activity increased in a concentration-dependent manner, being negatively correlated with MDA levels, indicating MDA metabolism by ALDH. In the gills, diesel exposure increased MDA and lipid hydroperoxide levels, and promoted increases in antioxidant defenses, indicating oxidative stress. Curiously, ALDH activity was undetectable in the gills, supporting the possibility of direct MDA excretion in the water by the gills. Analyses of MDA in the water revealed increased levels of MDA in the aquaria in which the fish were exposed to diesel, compared to control aquaria. A second experiment was carried out in which the fish were intraperitoneally injected with MDA (10 mg/kg) and analyzed after 1, 6, and 12 h. MDA injection caused a time-dependent decrease in hepatic MDA levels, did not alter ALDH, CAT, GPx, and GST activities, and decreased G6PDH activity and tGSH levels. In the gills, MDA injection caused a slight increase in MDA levels after 1 h, but did not alter GPx, G6PDH, and GST activities. MDA injection also enhanced CAT activity and tGSH levels in the gills. MDA concentration in water increased progressively after 1, 6, and 12 h, supporting the hypothesis of direct MDA excretion as an alternative route for MDA elimination in fish. Our results suggest that the decreased MDA levels after exposure of lambari to diesel oil pollutant probably reflects an association between enhanced antioxidant protection, MDA metabolism, and MDA excretion in water.

Metal bioaccumulation, oxidative stress and antioxidant responses in oysters Crassostrea gasar transplanted to an estuary in southern Brazil

The present study assessed the spatial and temporal variations on metal bioaccumulation and biochemical biomarker responses in oysters Crassostrea gasar transplanted to two different sites (S1 and S2) at the Laguna Estuarine System (LES), southern Brazil, over a 45-days period. A multi-biomarker approach was used, including the evaluation of lipid peroxidation (MDA) levels, and antioxidant defense enzymes (CAT, GPx, GR and G6PDH) and phase II biotransformation enzyme (GST) in the gills and digestive gland of oysters in combination with the quantification of Al, Cd, Cu, Pb, Fe, Ni and Zn in both tissues. The exposed oysters bioaccumulated metals, especially Al, Cd and Zn in gills and digestive gland, with most prominent biomarker responses in the gills. Results showed that GPx, GR and G6PDH enzymes offered an increased and coordinated response possibly against metal (Zn, Ni, Cd and Cu) contamination in gills. GST was inversely correlated to Cd levels, being its activity significantly lowered over the 45-d exposure periods at S2. On contrary, in digestive gland GST was slightly positively correlated to Cd, revealing a compensatory mechanism between tissues to protect oysters' cells against oxidative damages, since MDA levels also decreased. CAT also appeared to be involved in the cellular protection against oxidative stress, being increased in gills. However, CAT was negatively correlated to Al levels, which might suggest a possible inhibitory effect of this metal in the gills of C. gasar. Differences between tissues were evident by the Integrative Biomarker Responses version 2 (IBRv2) indexes, which showed different pattern between tissues when studying the sites and exposure periods separately. This study provided evidence for the effectiveness of using a multi-biomarker approach in oyster C. gasar to monitor estuarine metal pollution.

Stress responses in Crassostrea gasar exposed to combined effects of acute pH changes and phenanthrene

Ocean acidification is a result of the decrease in the pH of marine water, caused mainly by the increase in CO₂ released in the atmosphere and its consequent dissolution in seawater. These changes can be dramatic for marine organisms especially for oysters Crassostrea gasar if other stressors such as xenobiotics are present. The effect of pH changes (6.5, 7.0 and 8.2) was assessed on the transcript levels of biotransformation [cytochromes P450 (CYP2AU1, CYP2-like2) and glutathione S-transferase (GSTΩ-like)] and antioxidant [superoxide dismutase (SOD-like), catalase (CAT-like) and glutathione peroxidase (GPx-like)] genes, as well as enzyme activities [superoxide dismutase, (SOD), catalase (CAT), glutathione reductase (GR), glutathione-S-transferases transferase (GST) and glucose-6-phosphate dehydrogenase (G6PDH)] and lipid peroxidation (MDA) in the gills of Crassostrea gasar exposed to 100 μg·L^-1 of phenanthrene (PHE) for 24 and 96 h. Likewise, the PHE burdens was evaluated in whole soft tissues of exposed oysters. The accumulation of PHE in oysters was independent of pH. However, acidification promoted a significant decrease in the transcript levels of some protective genes (24 h exposure: CYP2AU1 and GSTΩ-like; 96 h exposure: CAT-like and GPx-like), which was not observed in the presence of PHE. Activities of GST, CAT and SOD enzymes increased in the oysters exposed to PHE at the control pH (8.2), but at a lower pH values, this activation was suppressed, and no changes were observed in the G6PDH activity and MDA levels. Biotransformation genes showed better responses after 24 h, and antioxidant-coding genes after 96 h, along with the activities of antioxidant enzymes (SOD, CAT), probably because biotransformation of PHE increases the generation of reactive oxygen species. The lack of change in MDA levels suggests that antioxidant modulation efficiently prevented oxidative stress. The effect of pH on the responses to PHE exposure should be taken into account before using these and any other genes as potential molecular biomarkers for PHE exposure.

Genome-wide identification and expression of the entire 52 glutathione S-transferase (GST) subfamily genes in the Cu2+-exposed marine copepods Tigriopus japonicus and Paracyclopina nana

In this study, the entire glutathione S-transferases (GSTs), the major phase II detoxification enzyme, were identified in two marine copepod species Tigriopus japonicus and Paracyclopina nana. The genome-wide identification of GSTs in T. japonicus and P. nana resulted in 32 and 20 GSTs in total, respectively. Among the identified GSTs, two specific classes of GSTs, specifically sigma and delta/epsilon GSTs were the dominant form of cytosolic GSTs in T. japonicus, while delta/epsilon and mu classes were dominant cytosolic GSTs in P. nana. In addition, Membrane-Associated Proteins in Eicosanoid and Glutathione metabolism (MAPEG) family were found in relatively higher proportion compared to other classes. Moreover, sigma, delta/epsilon, and microsomal GSTs have shown to expand through tandem duplication. To validate the detoxification function of the identified GSTs, both copepods were exposed to copper (Cu²⁺) and the reactive oxygen species (ROS) level and GST activity were measured. With integration of phylogenetic analysis and xenobiotic-mediated GST mRNA expression patterns along with previous enzymatic activities, the functional divergence among species-specific GST genes was clearly observed. This study covers full identification of GST classes in two marine copepod species and their important role in marine environmental ecotoxicology.