1 Introduction
Arsenic (As) contamination is a growing global concern due to its severe ecological and health impacts (
Appelo et al., 2002). In China, groundwater arsenic levels can reach up to 12 mg/L (
Kim et al., 2012). The long-term irrigation with contaminated water has increased soil arsenic concentrations, in Bangladesh, where levels rose from 57 to 83 mg/kg (
Alam and Sattar, 2000;
Panda et al., 2010). Due to its stability under aerobic conditions, As(V) is more commonly found in agricultural environments where crops are cultivated (
Wang et al., 2021b;
Cai et al., 2022). This persistent use of contaminated water for crop irrigation has led to significant arsenic accumulation in grains and vegetables (
Bednar et al., 2002;
Tripathi et al., 2007;
Zhu et al., 2008). When plants absorb inorganic arsenic, reactive oxygen species (ROS) such as hydrogen peroxide (H
2O
2) and hydroxyl radicals (•OH) are generated, causing oxidative stress (
Hartley-Whitaker et al., 2001). Consuming arsenic-contaminated crops is linked to various health issues, including metabolic disorders, kidney disease, cardiovascular conditions, cancer, and maternal-fetal complications (
Khan et al., 2021).
Apart from arsenic, the presence of emerging phenolic contaminants in the consumable portions of vegetables has raised growing health concerns (
Matschullat, 2000;
Meharg and Hartley-Whitaker, 2002;
Huang et al., 2019;
Yao et al., 2020;
Kovačič et al., 2023). However, the effects of arsenic on these emerging pollutants, particularly under conditions of ubiquitous iodine ion (I
−) presence, remain poorly understood. Bisphenol pollutants, such as bisphenol F (BPF), have garnered significant attention due to their widespread use and high bioaccumulation potential in plants (
Dodgen et al., 2013), BPF has been identified in a range of crops and fruits, with concentrations ranging from 3.3 to 4.6 mg/kg in
C. viridevar. and up to 7.2 mg/kg in
Bracteatum G. faberi (
Lu et al., 2013;
Huang et al., 2019). Compared to bisphenol A (BPA), BPF exhibits more pronounced adverse effects on spring rapeseed yield, showing significant inhibition in soil contaminated with 5 mg/kg of BPF (
Zaborowska et al., 2023). Despite being a BPA substitute, BPF still poses significant reproductive and genetic toxicity risks, evidenced by a high detection rate of 66.5% in the urine of both adults and children (
Lehmler et al., 2018;
Ding et al., 2022).
Phenolic pollutants like BPF can react with reactive halogens in the presence of free radicals to form more toxic organic pollutants. Iodine, ubiquitous in the environment, primarily exists as an inorganic iodine species (I
− and IO
3−) (
Hetzel, 2009). In soils, iodine concentrations range from 0.1 to 660 mg/kg, with a median value of 3 mg/kg (
Johnson, 2003;
Smyth and Johnson, 2011). Natural surface water iodine concentrations vary from 0.004 to 0.787 µmol/L, with some drinking water sources containing up to 32.417 µmol/L (
Tang et al., 2013). Groundwater generally exhibits higher iodine concentrations, reaching up to 114.173 µmol/L in Denmark, 267.716 µmol/L in Japan, and 48.000 µmol/L in Chile. In China, groundwater iodine levels are notably high, with concentrations of 36.228 µmol/L detected in Shandong and peaking at 222.519 µmol/L in the Guanzhong Basin (
Duan et al., 2016;
Wang et al., 2021c).
Iodide (I
−) readily transforms into reactive iodine species (RIS) such as hypoiodous acid (HIO), iodine radicals (I
•), iodide ions (I
2−), and triiodide ions (I
3−) through various pathways involving reactive oxygen species (ROS) and enzymatic reactions (
Buxton et al., 1988;
Milenković and Stanisavljev, 2012;
Chen et al., 2021;
Mackeown et al., 2022;
Wang et al., 2022a). RIS exhibits selective reactivity, primarily targeting phenolic compounds like BPF and
α-methyl carbonyl compounds such as acetaldehyde and acetone (
Bichsel and Von Gunten, 2000;
Ding and Zhang, 2009;
Li et al., 2018). The iodine-containing byproducts formed in these reactions are generally more toxic than other halogenated byproducts. Specifically, iodinated disinfection byproducts (I-DBPs) demonstrate more potent cytotoxicity and genotoxicity compared to their brominated or chlorinated counterparts (
Wagner and Plewa, 2017;
Dong et al., 2019;
Gonsioroski et al., 2020). Furthermore, arsenic and iodide co-contamination of groundwater is common in many regions (
Pi et al., 2015;
Wang et al., 2018). Arsenic stress in plants induces a burst of ROS (
Begum et al., 2016;
Tiwari and Sarangi, 2017). And a significant increase in the activity of various antioxidant enzymes (
Mishra et al., 2011;
Vezza et al., 2022). These biochemical changes can influence the transformation of I
− into RIS within plants, thereby affecting the iodination of phenolic pollutants.
This study employs arsenate (As(V)) as a representative metalloid pollutant and BPF as a representative phenolic pollutant to investigate 1) the effects of As(V) on the iodination of BPF in Brassica chinensis L. with I− and 2) the iodination pathways, identification of iodinated products, and evaluation of their ecological risks.
2 Materials and methods
2.1 Chemicals
BPF (≥ 98%), salicylhydroxamic acid (SHAM, ≥ 98%), methyl alcohol (≥ 99%), diphenyliodonium chloride (DPI, ≥ 98%), 4-Hydroxy Tempo (≥ 98%), and chromatography-grade acetonitrile (≥99%) were obtained from Macklin (Shanghai, China). Na2HAsO4·7H2O (≥ 98%) and Dimethyl-p-phenylenediamine (DPD, ≥ 98%) were purchased from Aladdin (Shanghai, China). Potassium iodide (≥99%), 30% hydrogen peroxide (≥ 30%), disodium hydrogen phosphate (≥ 99%), and sodium dihydrogen phosphate (≥ 99%) were sourced from Sinopharm Chemical Reagent Co., Ltd., Shanghai, China. Sample details are available in Table S1 of the supplementary material.
2.2 Cultivation and management of Brassica chinensis L.
The seeds of
Brassica chinensis L. were sterilized by immersing them in 10% H
2O
2 for 10 min, followed by soaking at 30 °C for 12 h. They were then germinated in the dark at 26 °C for two days. After germination, the seedlings were transferred to a hydroponic system containing 100% Hoagland nutrient solution under sterile conditions and cultivated in a sterile growth chamber. The plants were maintained under an 18-h light/6-h dark photoperiod at 26 °C/24 °C (
Dong et al., 2023). The nutrient solution was renewed every two days to prevent bacterial growth, and after two weeks, it was replaced with a 25% diluted solution. All procedures, including nutrient solution preparation and seedling transfer, were conducted aseptically to ensure a sterile environment. The specific operation details and cultivation parameters are provided in the supplementary material Text S1.
2.3 Experimental treatment of Brassica chinensis L.
The experiment used Brassica chinensis L. with stable growth for 45 d, exposed to a mixed nutrient solution containing I−, BPF, and different As(V) levels. Three groups were established by adding a peroxidase inhibitor (SHAM), NADPH oxidase inhibitor (DPI), and hydrogen peroxide scavenger (4-Hydroxy-Tempo), with samples collected for analysis after 5 d. The specific details of the exposure experiment and the selection of inhibitor concentrations are referenced in the supplementary material Text S2.
2.4 Assay methods
2.4.1 Homogenization of Brassica chinensis L. tissue samples
Plant root or leaf tissues were quickly frozen in liquid nitrogen and homogenized in 0.1 mol/L phosphate-buffered saline (pH 7.2, 4 °C) using a SCIENTZ-48LT grinder at a ratio of 0.1 g tissue per 1 mL buffer. The homogenate was centrifuged at 11500 g for 15 min at 4 °C (Shanghai Flying Pigeon, TGL-16G, China), and the supernatant was collected for further analysis.
2.4.2 Determination of reactive iodine species
Brassica chinensis L
. root and leaf samples (0.5 g) were homogenized, followed by centrifugation to collect the supernatant. The RIS (I
2, I
3−, and HIO) in the tissue samples were measured in this study using the DPD (
N,
N-diethyl-
p-phenylenediamine) method. A 5 mL sample was sequentially mixed with 0.5 mL of 0.5 mol/L phosphate buffer (pH 6.0) and 0.5 mL of standard DPD solution. Optical measurements were performed using a UV spectrophotometer (SHIMADZU, UV-2600, Japan), completed rapidly within 20 s. This method has been widely used for detecting RIS (
Brion and Silverstein, 1999;
Chen et al., 2021). The total amount of RIS (I
2, I
3−, and HIO) was quantified using the DPD method by measuring the formation of DPD
• +, which exhibits strong absorbance at 551 nm (ε = 21000 M
−1·cm
−1). Due to the complexity of the tissue sample, the presence of H
2O
2 in it can be rapidly oxidized by peroxidases, leading to the accuracy of its absorbance at 551 nm, and the presence of other organic substances, such as chlorophyll, may also affect the accuracy of absorbance at 551 nm. It is, therefore necessary to pre-treat the grinding solution of the tissue samples before measurement to ensure that the solution being measured is inorganic. In Text S3 and Fig. S1, we give the specific operation of the pretreatment step. We conducted preliminary pre-tests on the grinding solution of untreated
Brassica chinensis L. roots and leaves post-pretreatment to validate the effectiveness of the pretreatment process. Its experimental results showed that the pretreatment was effective in eliminating the interference of organic substances (chlorophyll, etc.) and substances such as peroxidase, and that its color development and absorbance were almost the same as those of deionized water (Figs. S2 and S3), which also ensured the accuracy of the subsequent RIS measurement by the DPD method.
2.4.3 Determination of hydrogen peroxide content, POD activity, and SOD activity, as well as the ability to produce O2•− and •OH
Fresh Brassica chinensis L. roots and leaves (0.1 g) were ground into a homogenate, and the supernatant was collected after centrifugation. The grinding solution was prepared according to the instructions provided with the Nanjing jiancheng reagent kit. The levels of H2O2, hydroxyl radicals (•OH), superoxide anions (O2•−), and the activities of peroxidase (POD) and superoxide dismutase (SOD) were assessed using the kit’s protocol.
2.4.4 Analysis of transformed products and assessment of BPF levels
Brassica chinensis L
. root and leaf samples (0.5 g) were homogenized, and the supernatant was collected after centrifugation. The extract underwent solid-phase extraction (SPE), was eluted with methanol, concentrated to 1 mL, and filtered through a 0.22 μm membrane. BPF transformation products and pathways were analyzed via LC-MS/MS (TSQ Quantum Access MAX, Thermo Scientific, USA) with negative electrospray ionization (ESI). The separation utilized a reverse-phase C18 column (50 mm × 2 mm, 4 μm) with modified HPLC conditions as described by
Yu et al. (
2020) and
Guo et al. (
2021). The UV detection was conducted at a wavelength of 243 nm using a mobile phase composed of acetonitrile and water in a 60:40 ratio, with a flow rate of 0.5 mL/min and an injection volume of 20 μL. The mass spectrometry setup featured a 3000 V capillary voltage, 8.0 eV capillary energy, 550.0 Vpp collision cell RF, and an MS scan range of 60–500 m/z. The content of BPF in fresh root and leaf samples was quantified using HPLC under the same pretreatment and chromatographic conditions as LC-MS/MS.
2.5 Transcriptomics and q-PCR analysis
Brassica chinensis L. was subjected to various treatments for 5 d, after which 0.2 g samples of root or leaf tissues were collected, thoroughly washed with deionized water, flash-frozen in liquid nitrogen, and stored at −80 °C. RNA was isolated from tissue samples, and its concentration and purity were assessed with a spectrophotometer (NanoDrop 2000, Thermo Scientific, USA). RNA integrity was evaluated using agarose gel electrophoresis, and the RNA integrity number (RIN) was measured with a bioanalyzer (Agilent 2100, Agilent, USA). mRNA was extracted from total RNA using magnetic beads with Oligo (dT), which binds to the polyA tail of mRNA via A-T base pairing. Sequencing was performed on the Illumina platform (Shanghai Meiji Biotechnology, China) targeting short-read sequences. A fragmentation buffer was introduced to randomly cleave the mRNA into ~300 bp fragments, which were then isolated via magnetic bead screening. Double-stranded cDNA was synthesized by first generating first-strand cDNA from mRNA using reverse transcriptase and random hexamer primers, followed by second-strand synthesis. The double-stranded cDNA’s sticky ends were converted to blunt ends using End Repair Mix, followed by the addition of an ‘A’ base to the 3' ends to facilitate ligation with Y-shaped adapters. Sequencing was performed on the Illumina platform, and differentially expressed genes shared between roots and leaves were selected for q-PCR analysis. Real-time quantitative PCR was conducted using the Applied Biosystems StepOne™ Real-Time PCR system, with primer sequences detailed in Supplementary Table S2. Reference genomes were from Ensembl Genomes.
2.6 Toxicity assessment
The ecological toxicity variations of BPF and its iodinated products were comprehensively predicted using two software tools: Ecological Structure Activity Relationships (ECOSAR, version 2.2) and the Toxicity Estimation Software Tool (T·E·S·T, version 5.1.2) (
Chen et al., 2018;
Svedberg et al., 2023;
Xu et al., 2024). Both tools utilize advanced QSAR models to predict compound toxicity from molecular structures. These tools have been internationally recognized and validated for predictions (
Claeys et al., 2013;
Khan et al., 2019). The ECOSAR model assessed the acute and chronic toxicity of BPF and its iodinated derivatives on three aquatic species: fish, daphnia, and green algae. The T·E·S·T software evaluated the acute toxicity of BPF and its iodinated derivatives on four organisms: fathead minnow (Pimephales promelas), Daphnia magna, and Tetrahymena pyriformis. Toxicity parameters encompassed LC
50 (median lethal concentration), IGC
50 (50% inhibitory growth concentration), developmental toxicity, and mutagenicity.
2.7 Cytotoxicity analysis
HepG2 cells were obtained from the China Cell Bank. Detailed procedures for cell culture and passaging are outlined in Text S4 (
Hercog et al., 2019;
Wang et al., 2021a;
Yu and Liu, 2023). A slurry was prepared by homogenizing 0.1 g of fresh plant roots and leaves, followed by centrifugation. The supernatant was then collected and filtered through a 0.22 μm membrane to eliminate macromolecular proteins. To investigate the toxicity of organic compounds in the plant samples, the supernatant was passed through an SPE solid-phase extraction column, eluted with methanol, and then used to expose HepG2 cells. HepG2 cells were treated with a 0.5% concentration of both the organic extract solution and the self-prepared BPF solution for 48 h (Fig. S4). This concentration ensured uniformity and minimized the potential toxicity of methanol as a solvent (
Nguyen et al., 2020). The BPF concentrations in the prepared solutions matched those found in 0.1 g of homogenized leaf or root tissue per 1 mL of grinding buffer for each treatment group. The Nanjing Jiancheng Bioengineering Institute assay kits were used to measure cell viability (CCK-8 assay) and lactate dehydrogenase (LDH) activity in the cell culture medium.
2.8 Data analysis
Each experimental treatment was performed three times, and the outcomes were presented as the average ± standard deviation (SD). Statistical analysis was conducted using SPSS 25.0, with one-way ANOVA applied to assess significance at
p < 0.05. Post-hoc tests were conducted using LSD, Tukey’s b, and Waller-Duncan methods.
t-tests with Studentized range distributions were used to evaluate pairwise comparisons of group means, and multiple comparisons were also conducted using
t-statistics. Moreover, Bayesian techniques were utilized for further statistical analysis (
Gu et al., 2020;
Kundu et al., 2024). Gaussian 09W calculates 2FED
2HOMO per carbon before iodine electrophilic substitution.
3 Results and discussion
3.1 As(V) exposure enhances the iodination of BPF in Brassica chinensis L.
3.1.1 The iodinated products of BPF in Brassica chinensis L.
Fig.1(a)–Fig.1(c) and Fig.1(d)–Fig.1(f) display the full scan mass spectra of BPF transformation products in the roots and leaves of cabbage exposed to 0, 25, and 100 μmol/L As(V), with different iodinated products highlighted in different colored boxes (the corresponding EIC plots are shown in the supporting information, Fig. S5). The identification of these iodine-containing transformation products (namely TP1–TP7) is based on the characteristic ion m/z 127 found in the secondary mass spectrometry results (Figs. S6–S12), and the m/z values of iodine-containing ions detected in the samples are listed. The structures of the substances were inferred based on the secondary mass spectrometry data. For example, TP5 is inferred to be a substance containing two iodine because molecular ion peaks are observed at m/z 211 and 210 which differ from the precursor ion by 254 and 255 Da, indicating the loss of two iodine atoms or one iodine and one HI, and, molecular ion peaks are found at m/z 127 and 128, which further confirms that it contains two iodine. At the same time, it was also observed that the molecular ion peak containing m/z was 232, which was 233 different from the precursor ion, that is, C7H6OI, which was believed to be dropped by the breakdown of the carbon bridge bond formed by two carbons in TP5, and the molecular ion peak with a difference of 127 m/z was found to be 105 with the molecular ion peak of 232, which proved that the molecular ion peak of 232 also contained iodine. These data support our description of the structure of TP5 (Fig. S10). Due to the lack of standards, further structural confirmation could not be performed, so theoretical calculations were used to determine the specific positions of iodine substitution. The carbon atoms that are most likely to be electrophilically substituted by BPF, TP2-1, TP3-1, TP4-1, TP4, TP6-1, and TP7-1 are shown in Fig. S13. The corresponding values are detailed in Table S3.
3.1.2 As(V) increased the type and content of iodinated products in Brassica chinensis L.
As shown in Fig.2(a), the number of iodinated products in roots increased with rising As(V) concentrations. Notably, the peak areas of newly formed TP5 and TP6 in roots also increased with higher As(V) concentrations, and at 100 μmol/L As(V), a root-specific iodinated product, TP7, was observed. Similarly, in the leaves, the types of iodinated products increased after adding As(V), and the peak area of iodinated product TP2 increased with the addition of As(V) concentration, reaching its maximum at 100 μmol/L As(V) (Fig.2(b)).
Combined with the results of Fig.1 and Fig.2, the addition of As(V) increases the type and number of iodide products in the roots, while the diversity of iodide products in the leaves enriches the diversity of iodide products, although the iodine products formed in the roots and leaves are not the same. Thus, the addition of As(V) enhances the iodination of BPF.
3.1.3 The generation pathway of iodinated products during the BPF transformation in Brassica chinensis L.
Based on the m/z values in the primary mass spectra and a review of relevant literature, the iodinated products of BPF may be formed through the following three generation pathways: (I) BPF can be directly electrophilically replaced by RIS to form the iodinated product TP1. (II) BPF may react with •OH radicals to initially produce intermediates A1-1 or A1-2 (m/z 231), which then convert into TP2-1 (m/z 247) before undergoing substitution by RIS to form the iodinated product TP2. This hydroxylation process is well-documented and has been studied in previous research (
Morales-Roque et al., 2009;
Porcar-Santos et al., 2022). (III) Reactive species within plant tissues may directly oxidize to attack BPF, resulting in
β bond cleavage and the formation of
p-hydroxybenzyl alcohol (A2, m/z 123) and phenol (A3). As phenolic compounds, A2 and A3 may undergo chemical oxidation or react with peroxidase (in H
2O
2) to form phenoxy radicals. These phenoxy radicals can then be coupled at different locations to produce dimers or oligomers (
Henriksen et al., 1999;
Laurenti et al., 2003;
Yang et al., 2020;
Jia et al., 2022;
Liu et al., 2022). A2 may undergo self-coupling to form TP4-1 (m/z 213), which can then be substituted by reactive iodine species (RIS) to produce TP4. The diiodinated product TP5 may arise through two pathways: further electrophilic substitution of TP4 by RIS or direct electrophilic substitution of TP4-1 by additional RIS. The formation of TP6 is likely initiated by the self-coupling of A3 to form TP6-1, followed by electrophilic substitution with RIS. In addition, A3 may promote TP7 production through two possible pathways: direct coupling to TP1 or conjugation to BPF to form TP7-1, which is then substituted by RIS electrophiles, leading to the formation of TP7 (Fig.3).
3.2 Iodination mechanisms of BPF with As(V) exposure in Brassica chinensis L.
3.2.1 As(V) promotes the iodination of BPF by elevating POD activities in Brassica chinensis L.
Numerous studies have shown that As(V) contamination can induce changes in the activities of various antioxidant enzymes in plants, including POD, SOD, catalase (CAT), and glutathione S-transferase (GST) (
Mascher et al., 2002;
Bianucci et al., 2017;
Tzean et al., 2024). Without As(V), the POD activity in both roots and leaves did not show significant increases. However, upon the addition of As(V), POD activity in roots and leaves increased with rising As(V) concentrations, with the changes in root POD activity being more pronounced than those in leaves (Fig.4(a)). Previous studies have demonstrated that Class III peroxidases (POD), such as horseradish peroxidase (HRP), can catalyze the transformation of I
− into RIS (HIO, I
2, I
3−, etc.) in the presence of H
2O
2, which then react with organic pollutants to produce iodinated products (Tab.1, Eqs. (1)–(3)) (
Wang et al., 2022a).
The levels of RIS were measured, revealing that after the addition of As(V), RIS in both roots and leaves increased with rising As(V) concentrations, peaking at 100 μmol/L As(V) (Fig.4(b)). Since previous studies have shown that POD can oxidize I− to RIS, we infer that As(V) may promote the production of RIS in the roots and leaves of Brassica chinensis L. by stimulating the activity of POD. To further confirm this, SHAM, a potent POD inhibitor, was introduced to the group with the highest RIS production. The results showed a reduction of RIS by 36.12% in roots and 13.85% in leaves (Fig.4(c)). These findings indicate that As(V) likely promotes the conversion of I− to RIS by enhancing POD activity in the roots and leaves of Brassica chinensis L., thereby increasing the iodination of BPF within the plant.
3.2.2 As(V)-induced oxidative stress enhancing the iodination of BPF in Brassica chinensis L.
3.2.2.1 The role of hydrogen peroxide
After adding As(V), the H
2O
2 levels in both the roots and leaves increased to varying degrees, reaching a peak at 100 μmol/L As(V) (Fig.5(a)). Although exposure to BPF alone also caused a slight increase. H
2O
2 is typically generated through the conversion of O
2•− under the action of superoxide dismutase (SOD) (
Quan et al., 2008). After the addition of As(V), the SOD activity and the ability to produce O
2•− in both roots and leaves increased to varying degrees, which explains the elevated H
2O
2 levels in different tissues of
Brassica chinensis L. after As(V) exposure (Fig. S14). Notably, H
2O
2 levels in leaves were higher than those in roots. This could be attributed to the higher POD activity in roots compared to leaves, leading to the conversion and decomposition of more H
2O
2 in root tissues (Fig.4(a)). H
2O
2, as the longest-lived ROS (
Bhattacharjee, 2005;
Quan et al., 2008), possesses sufficient oxidative capacity to convert I
− into RIS (Tab.1, Eqs. (2)–(5)) (
Chen et al., 2021;
Wang et al., 2022a;
Xu et al., 2023;
Shao et al., 2024). To further verify its role, 4-Hydroxy Tempo, an effective H
2O
2 scavenger (
Wang et al., 2022b;
Deng et al., 2023), was added to the experimental group with the highest RIS production. This resulted in a 40.73% reduction in RIS levels in roots and a 19.32% reduction in leaves (Fig.5(b)). Compared to the use of peroxidase inhibitors, the more pronounced suppression of RIS levels may be attributed to 4-Hydroxy Tempo not only eliminating H
2O
2 but also indirectly limiting the POD-mediated pathway for RIS production. These findings suggest that As(V) likely enhances the iodination of BPF in
Brassica chinensis L. by increasing H
2O
2 levels in roots and leaves, thereby promoting RIS production.
3.2.2.2 The role of hydroxyl groups
As shown in the Fig.6(a), in the absence of As(V), there was no significant variation in the ability of leaves to produce •OH, while roots showed a slight increase, though not markedly different. After the addition of various concentrations of As(V), the ability to produce •OH increased to varying degrees in leaves, whereas in roots, a significant increase was observed only at a concentration of 100 μmol/L As(V) (Fig.6(a)).
O
2•−, as the initial reactive oxygen species, can be converted into H
2O
2, which can then generate hydroxyl radicals (•OH) through photolysis in the Haber-Weiss reaction and reactions with transition metals (Eqs. (6) and (7)) (
Chen and Schopfer, 1999;
Richards et al., 2015). Therefore, it plays a critical role in the formation of •OH. •OH, as the most oxidizing species among ROS, can easily convert I
− into various RIS (I
•, I
2•−, I
3•−, I
2, HIO, I
3−) (Eqs. (2), (3), (5) and (8)–(13)), however, O
2•− itself cannot directly oxidize I
−. Previous studies have shown that NADPH oxidase is involved in arsenic-induced ROS production in various plants, including peanut (
Bianucci et al., 2017), rice (
Kushwaha et al., 2019), and tomato (
Siddiqui et al., 2024). In plants, NADPH oxidase (respiratory burst oxidase homologs) serves as a major source of O
2•−, and its impairment can significantly reduce •OH generation (
Renew et al., 2005;
Glyan’ko and Ischenko, 2010;
Heyno et al., 2011).
Therefore, the effective NADPH oxidase inhibitor DPI was added in the experimental group with the highest RIS production. After the addition of DPI, RIS levels in roots decreased by 40.54%. In comparison, in leaves they decreased by 32.25% (Fig.6(b)). In leaves, the inhibition rate of RIS production with DPI was further improved compared to that with 4-Hydroxy Tempo, likely due to the additional suppression of •OH-mediated oxidation of I− to RIS, on top of inhibiting POD and H2O2 pathways. However, in roots, the inhibition rate of RIS production with DPI was not further enhanced compared to 4-Hydroxy Tempo. This phenomenon occurs because the increase in •OH production in roots after As(V) exposure is not as significant as in leaves, resulting in a less pronounced contribution of •OH to RIS formation (Fig.6(a)). These findings suggest that As(V) may enhance the iodination of BPF in Brassica chinensis L. by increasing the production of •OH in roots and leaves, which in turn boosts RIS formation.
3.3 Using PCA analysis to conduct a correlation analysis on the mechanism of enhancing Brassica chinensis L. BPF iodization by As(V)
Principal component analysis (PCA) has advantages in dimensionality reduction, multicollinearity, revealing global structure, visualization and synthesis of information, and is suitable for high-dimensional data analysis and pattern recognition. However, the Spearman and Pearson correlation is more suitable for analyzing the specific relationship between the two variables, so the PCA analysis was used to explore the relationship between the ability of roots and leaves to produce •OH, H2O2 content, POD activity, and RIS changes in the experimental group with different As(V) concentrations. The results indicate that in both roots and leaves, H2O2, •OH, and POD were positively correlated with RIS (Fig.7). In roots, the increase in RIS production with rising As(V) concentrations was more closely related to changes in H2O2 content and POD activity, while •OH also showed a positive correlation but appeared to play a less dominant role, aligning with our earlier hypothesis based on the results of adding DPI. H2O2 likely plays a critical role in root RIS production due to its dual capacity to directly oxidize I− into RIS and to serve as a substrate for POD to convert I− into RIS (Fig.7(a)). In leaves, POD did not appear to dominate RIS production, while •OH and H2O2 were more strongly associated with RIS. This may be due to the lower POD activity in leaves than roots, leaving more H2O2 unutilized, indirectly promoting •OH production. Both H2O2 and •OH can independently oxidize I− into RIS, explaining the formation of the polyhydroxylated iodinated product TP2 in leaves, which was not observed in roots. Compared to H2O2, •OH possesses stronger oxidative power and reacts with I− more rapidly, which may explain its closer correlation with RIS changes in response to As(V) concentrations (Fig.7(b)). Furthermore, the correlations among components in roots under different As(V) treatments were more pronounced than those in leaves, indicating a more distinct relationship in roots.
3.4 Transcriptomic analysis and validation of the enhancement of BPF iodination mechanism in Brassica chinensis L. by As(V)
A transcriptomic analysis was conducted to further explore the mechanism by which As(V) enhances BPF iodination in Brassica chinensis L. The results showed that as As(V) concentration increased, the number of differentially expressed genes (DEGs) in the roots and leaves reached 3824 and 1618, respectively, accounting for 20.43% and 15.55% of the total DEGs in the three groups of roots and leaves (Fig.8(a) and Fig.8(b)). These proportions are relatively high, indicating that adding As(V) can cause significant changes and differential impacts on gene expression in Brassica chinensis L.
Based on the principle of prioritizing genes commonly expressed in both roots and leaves (highlighted with boxes), followed by genes regulating the same proteins (annotated with identical text), and finally genes with similar regulatory functions (indicated by the same text color), potential genes enhancing the iodination of BPF in Brassica chinensis L. were selected for analysis. The results showed that genes in both roots and leaves responsive to H2O2, ROS, and those regulating POD and SOD were predominantly upregulated to varying degrees after As(V) treatment. This further supports the proposed three pathways by which As(V) enhances the iodination of BPF in Brassica chinensis L.
At the same time, it can be observed that genes regulating POD in roots exhibit a stronger upregulated response to increasing As(V) concentrations compared to those in leaves. The shared gene
Bra017830 in both roots and leaves highlights this phenomenon more clearly, aligning with our PCA results (Fig.8(c)). This further suggests that the POD pathway may be a crucial mechanism by which As(V) enhances BPF iodination in roots. Notably, since POD can couple phenoxy radicals to form coupled products (
Laurenti et al., 2003;
Liu et al., 2022), this further explains why at 100 μmol/L As(V), the root contains the coupled iodinated product TP7, which is absent in leaves.
Genes responsive to ROS were more strongly upregulated in leaves than in roots, as indicated by the upregulation of three shared genes
(Bra024875,
Bra005687,
Bra000964). Since •OH is a type of ROS, these gene expression results further support our PCA findings and help explain why the polyhydroxylated iodinated product TP2 was present in leaves but absent in roots after As(V) treatment. Notably, genes regulating monooxygenases (such as P450 enzymes) were downregulated in roots but strongly upregulated in leaves following As(V) exposure, as observed in the shared genes
Bra008510,
Bra020484, and
Bra015998. Monooxygenases, including P450 enzymes, are known for their ability to hydroxylate diphenol-like compounds efficiently (
Wang et al., 2020;
Zhou and Hong, 2021), and their upregulation may also contribute to the formation of TP2 in leaves but not in roots (Fig.8(d)).
To verify the reliability of the transcriptomic results, we selected four genes (Bra008510, Bra017830, Bra024875, and Bra038612) commonly expressed in roots and leaves, as shown in Fig.8(c) and Fig.8(d), for q-PCR validation. The results were almost identical to the transcriptomic findings (Fig. S15), further confirming the reliability of the transcriptomic analysis.
Based on the two analyses mentioned above, the mechanism by which As(V) enhances the iodination of BPF in Brassica chinensis L. likely involves three pathways: increasing POD activity, enhancing H2O2 content, and boosting the ability to produce •OH. POD activity and H2O2 content dominant roots, while the effects on •OH and H2O2 are more significant in leaves.
3.5 ECOSAR and T·E·S·T software toxicity evaluation of BPF iodinated products in Brassica chinensis L.
Seven iodinated products were analyzed for toxicity using ECOSAR and T·E·S·T software. The results revealed that TP1, TP4, TP5, TP6, and TP7 exhibited increased acute toxicity to most species. Notably, TP2 showed contradictory toxicity trends: an increase in T·E·S·T analysis and a decrease in ECOSAR analysis, likely due to differences in target species.
Whether using ECOSAR analysis or T·E·S·T analysis, the most toxic substance in the roots is TP5, which is generated after the addition of As(V) and increases with the concentration of As(V). ECOSAR analysis revealed that the toxicity of TP5 to different species increased by 2.57 to 96.30 times compared to BPF. In the T·E·S·T analysis, the toxicity also increased by 10.66 to 93.67 times. The second most toxic substance in the roots is TP7, which is produced at a concentration of 100 μmol/L. Toxicity analysis using ECOSAR and T·E·S·T found that its toxicity increased by 2.98 to 32.50 times and 7.89 to 26.30 times, respectively. In the leaves, TP1 and TP6, newly generated after the addition of As(V), showed an increase in toxicity during the ECOSAR and T·E·S·T evaluations by 1.53 to 4.85 times, 2.19 to 4.01 times, and 1.30 to 3.20 times, 1.68 to 3.01 times, respectively. These findings indicate that the iodinated products in Brassica chinensis L. influenced by As(V) could present greater environmental hazards and potentially impact human well-being (Tab.2 and Tab.3).
3.6 Toxicity of iodinated products extracted from different tissues of Brassica chinensis L. to HepG2 cells
Previous studies have indicated that bisphenol analogs exhibit cytotoxicity toward HepG2 cells (
Vidyashankar et al., 2014;
Geng et al., 2017;
Hercog et al., 2019;
Yue et al., 2019). Therefore, to further validate the toxicity of iodinated products in
Brassica chinensis L., HepG2 cells were used for toxicity assessment in this study.
It was found that BPF levels in the different As(V) treatment groups did not cause any significant change in the viability of HepG2 cells (Fig.9(a) and Fig.9(b)). Similarly, lactate dehydrogenase (LDH), is a stable cytosolic enzyme rapidly released into the culture medium upon cell membrane rupture due to cell death (
Horinouchi et al., 2015;
Zhang et al., 2016). The activity of LDH in the cell culture medium was also not significantly altered by the changes in BPF content in roots and leaves after exposure to word self-prepared BPF solution (Fig.9(c) and Fig.9(d)). When As(V) was added, organic extracts from the roots resulted in a further decrease in HepG2 cell viability and an increase in LDH levels in the cell culture medium, with a significant difference in effect at 100 μmol/L. This indicates higher levels of harmful organic substances produced in the roots, a phenomenon attributed to the increased content of the iodide products TP5 and TP6 in the roots, as well as the production of TP7 at 100 μmol/L. Similarly, the organic extract from the leaves also resulted in varying decreases in cell viability and increases in LDH levels after As(V) treatment. This indicates that more types of harmful organic substances are produced in the leaves, namely the toxic iodide products TP1 and TP6 formed in the leaves when As(V) is added.
4 Conclusions
Arsenate (As(V)) significantly enhances the iodination of bisphenol F (BPF) in Brassica chinensis L. by inducing bursts of ROS, specifically hydrogen peroxide (H2O2) and hydroxyl radicals (•OH), and by increasing POD activity. H2O2, due to its longer lifespan, plays a crucial role in As(V)-mediated BPF iodination in roots and leaves. However, the underlying mechanisms differ: in roots, iodination is predominantly driven by elevated POD activity, whereas in leaves, it relies more on •OH bursts. Furthermore, exposure to As(V) results in the formation and accumulation of toxic iodinated byproducts, including TP5 and TP6 in roots, with the highly toxic TP7 emerging at 100 μmol/L As(V). In leaves, As(V) exposure leads to the production of additional toxic iodinated products, TP1 and TP6. These toxic compounds may be released into the soil or pore water via the root system, posing significant ecological risks. Moreover, these iodinated byproducts can enter the human body through dietary intake, directly or indirectly, potentially leading to more severe health hazards.
This study highlights the dual threat of arsenic contamination and the enhanced formation of toxic iodinated pollutants, underscoring the need for comprehensive risk assessments and mitigation strategies to protect both environmental and public health.