Methodology used is given in detail in Supplementary materials.
2.1. Chemicals
All reagents were of the Reagent grade or higher (S1 Appendix). They were obtained from Sigma-Aldrich (USA) or the Synbias (Ukraine). Roundup (Rn) formulation was Roundup MAX, Monsanto, USA, and chlorpromazine (Cpz) was of pharmaceutical grade, AMINAZINUM, Pat "Halychfarm", ATХ N05A A01.
2.2. Experimental groups
Adult bivalve molluscs Unio tumidus Philipson, 1788 (Unionidae) (~6 years old, ~8.5 cm length, and 60–70 g weight) were collected in a river site assumed to be a reference (Gnatyshyna et al. 2020). Specimens were transported to the laboratory and preacclimated to the laboratory conditions for up to seven days after the capture in the aerated, dechlorinated, softened tap water and fed 500 mg of Tropical SuperVit Basic contained beta-1.3/1.6-glucan twice a week. After that, molluscs were distributed randomly to five groups. The first group was exposed to the aquarium water only and was considered control (C). The Rn- and RnT-groups were exposed to organophosphonate pesticide Roundup MAX (17 µg L−1, corresponding to 6.1 µg L−1 or 36 nM of glyphosate) at the temperatures 18 oC and 25 oC, respectively. The temperature was increased gradually in the RnT‑group during 24 h. The Cpz-group was exposed to 18.0 µg L- 1 or 56 nM of Cpz, and the RnCpz-group - to a mixture of Rn and Cpz at 18 oC. The duration of exposure was 14 days. Water was changed and chemicals replenished every two days. Throughout the experiment, molluscs were fed with the same regularity. No mortality was detected during the experimental exposures.
After exposures, molluscs were dissected on ice. For lisosomal membrane stability and cholinesterase activity, the samples were utilised immediately. For all analyses except metallothioneins and caspase-3, tissues were homogenized (10% w/v) in 0.1 M phosphate buffer, pH 7.4, containing 100 mM KCl and 1 mM EDTA, as well as 0.1 mM phenylmethylsulfonyl fluoride (PMSF) for proteolysis inhibition. Homogenates were centrifuged at 6000g for 10 min, and the resulting supernatant was kept at -40 oC. For the analysis of cholinesterase, the 10% w/v homogenate was prepared in the same buffer without PMSF. The protein concentration was analysed in the 6000g supernatant according to the method of Lowry et al. (1951), using bovine serum albumin as the protein standard.
2.3 Assays for metallothionein and glutathione
Concentration of metallothioneins protein (MT) was assessed using 5,5’-dithio-bis(2-nitrobenzoic acid) (DTNB) reduction method after the ethanol/chloroform extraction from tissue 1:3 w/v homogenate in 20 mM Tris-sucrose buffer with 0.1% ß-mercaptoethanol, 0.5 mM PMSF, 6 µM leupeptine (Viarengo et al. 1997). The concentration of MT was expressed in µg g-1 tissue FW.
For the evaluation of Zn concentration in metallothioneins (Zn-MT), 70 mg of tissue per individual (to a total of 350 mg from five specimens in the group) were combined. Two replices for each group were accomplished. The samples were homogenized in 10 mM Tris-HCl buffer, pH 8.0 containing 10 mM β-mercaptoethanol and 0.1 mM PMSF and subjected to the isolation of thermostable supernatant. The MT‑contained fractions were isolated from the supernatant by size-exclusion chromatography on Sephadex G-50 with necessary adjustments needed to avoid their oxidation (Roesijadi and Fowler 1991). Low weight (approximately 7 kDa) fractions with high absorbance at 254 nm and high D254/D280 density ratio were identified as putative MTs-containing peak and pooled (to the total of 10 mL) for the Zn determination.
Total glutathione and oxidized glutathione (GSSG) concentrations were quantified by the glutathione reductase recycling assay (Griffith 1980) in the protein-free extract of homogenate using DTNB. Concentration was expressed as nmol g-1 FW. The concentration of the reduced glutathione (GSH) was calculated as the difference between the total glutathione and GSSG concentrations. The redox-index of glutathione (RI GSH) as the ratio of concentrations GSH/GSSG was calculated.
2.4. Oxidative stress and toxicity assays
Total antioxidant capacity (TAC) was determined as ABTS radical scavenging activity (Re et al. 1999). ABTS•+ radical cations (ABTS*) were generated from 2,2’-Azino-di-[3-ethylbenzthiazoline sulphonate] (ABTS) by potassium persulfate. The ascorbic acid was used as the reference compound. The reduction in absorbance of ABTS* solution was recorded at 734 nm. The result was compared with control (only ABTS* solution).
Protein carbonyls (PC) were determined as an index of protein oxidation in the sediment of 10% w/v homogenate in sulfosalicylic acid after its centrifugation by reaction with 2,4-dinitrophenylhydrazine (DNPH) (Reznick and Packer 1994). The concentration of carbonyls was expressed in nmol PC per g FW.
Cholinesterase (ChE, EC 3.1.1.7) activity was determined in the homogenate according to the colorimetric method of Ellman et al. (1961) at 25 oC. Acetylcholine iodide (ATCh) was used with DTNB as the thiol indicator. Enzyme activity was referred to the protein content.
Lysosomal membrane integrity was determined by Neutral Red retention (NRR) test based on the lysosomes ability to concentrate the dye as it was described in El Haj et al. (2019). The tissue samples (30 mg) were incubated for 2 h with a saline solution containing NR, washed and fixed in formaldehyde (0.5% in 1% CaCl2). After fixation, the formaldehyde impregnated tissue fragments were removed and frozen (up to one week). Dye was extracted in acid alcohol (1% acetic acid in 50% ethyl alcohol) and analyzed spectrophotometrically at 550-nm.
The ability of exposures to induce cytochrome P450 (CYP450) activity was quantified as 7-ethoxyresorufin O-deethylase (EROD) activity in the supernatant of 10% w/v homogenate by measuring the formation of resorufin at 572 nm (Klotz et al. 1984). The reaction was initiated by the addition of 0.5 mM NADPH. EROD activity was calculated using a molar extinction coefficient of 73.2 103 М–1·сm–1 and referred to the soluble protein concentration.
Glutathione S-transferase (GST, EC 2.5.1.18) activity was assayed using GSH and 1-chloro-2,4-dinitrobenzene (CDBN) as the substrate (Habig et al. 1974). The GST activity was expressed in nmol min-1·mg-1 soluble protein.
Caspase-3 (EC:3.4.22.56) activity as the marker of apoptosys was assayed colorimetrically in 25% w/v homogenate of digestive gland tissue based on the hydrolysis of peptide acetyl-Asp-Glu-Val-Asp p-nitroanilide (Ac-DEVD-pNA) by caspase-3 that produces a colored product p-nitroaniline (pNA). The pNA was detected at 405 nm (εmM = 10.5 mM‑1⋅cm‑1) (Du et al. 1997). The activity of caspase-3 was expressed in nmol pNA min-1 mg‑1 of soluble protein.
Citrate synthase (CS, EC 2.3.3.1) activity was measured in the 10% w/v homogenate of the digestive gland according to Flynn et al. (2015) as the maximum rate of increase in absorbance at 412 nm, caused by the production of a coenzyme A-SH and monitored by DNTB. CS enzyme activity was calculated by subtracting the background activity (negative control) from the CS enzyme activity (positive reaction) for each sample and quantified using the molar extinction coefficient of DTNB (14,150·M-1 cm-1).
2.5. Zinc concentration in the tissue and metallothioneins
The concentration of Zn was measured in the digestive gland tissue (Zn-t) and the pooled MT-containing eluate (Zn-MT) received by chromatography utilizing the reaction of the complexation of Zn(II) with 2-(5-bromo-2-pyridylazo)-5-[N-propyl-N-(3-sulfopropyl) amino]phenol disodium salt dihydrate (5-Br-PAPS) (Wang et al. 2018). The samples were dried at 105 °C for 24 h and then digested with HNO3. The received ash was dissolved in 1% trichloroacetic acid, then 0.6 M 5-Br-PAPS in carbonate buffer (pH 8.9) was added. The mixture was incubated for 30 min at 200 C. Zn concentration was evaluated from the absorbance of the metal-5-Br-PAPS complex at 560 nm. Fe ions were masked with citric acid, and Cu ions were masked with salicylaldoxime, deferoxamine, and sodium citrate. The detection limit was 0.1 µg·g–1 FW. Quality control was performed by the method of Standard Addition. Metal concentration in the tissue and MTs was expressed as nmol or μg·g-1 FW.
2.6. Statistical analysis
Results were expressed as mean ± standard deviation. Metal in MTs analysis was repeated in four samples for each of two independent combined from five individuals replices in a group, resulting in n = 8 for each group. For all other traits, the sample size was eight from eight individuals. Shapiro-Wilk test was used for the assessment of normality. Data were analyzed with parametric Student`s t-test significant at p < 0.05. Principal component analysis (PCA) was performed to assess the relations between measured parameters utilizing the rotation method Varimax. The adequacy of data was evaluated based on the value of the KMO and Bartlett’s test of sphericity. Canonical discriminant analysis was utilised for the separation of the exposed groups. The IBM SPSS Statistics version 26 software for Windows was used for calculations. Correlation was significant at p < 0.05 level (r > 0.304) and p < 0.01 (r > 0.393) (2-tailed), n = 40.
Integrated Biomarker Index (IBR) elaborated by Beliaeff and Burgeot (2002) was calculated for biomarkers (totally 14). The indices of CAP, EROD, GST, protein carbonyls (PC), GSH, GSSG, RI GSH, metallothionein concentration (MTSH), lysosomal integrity (NRR), caspase-3 activity (Cas-3), cholinesterase activity (ChE), citrate synthase (CS), Zn-metallothionein (Zn-MT) and Zn total (Zn-t) were used for the computation of data. The standardization of data was achieved by Xi calculation: Xi = (Mi-Mt)/SDi. For Ai calculation, the computation of data was made as to the Xe-Xc, assuming that the changes of the value in each exposed group (Xe) in relation to control (Xc) is corresponding to the stress or toxicity responses, and the value of standardized marker Xc was adjusted to zero. IBR value for each group was calculated as [(A1 × A2) + (A2 × A3) + … (An ×·A1)]/2. Since the value of IBR is dependent on the number of markers, the termed IBR value was given as IBR/n with n = 14.