Aquatic risk assessment of a novel strobilurin fungicide: A microcosm study compared with the species sensitivity distribution approach

Ecotoxicity and Risk to Zooplankton of a Novel Insecticide Used in Rice Paddies in a Simulated Ecosystem

Diamide insecticides are widely used in rice paddies and pose a potential threat to aquatic organisms. However, the risk research related to their application in major rice-producing areas is very limited, especially mesocosm research to simulate the impact on aquatic ecosystems of long-term exposure, as well as exposure analysis based on local models and local scenarios. To assess potential risks from a novel diamide insecticide (tetrachlorantraniliprole) to aquatic nontarget organisms in the field over long-term exposure, an outdoor mesocosm study was performed, and the environmental concentrations were predicted by the multimedia paddy-pond model (TOPRICE). The mesocosm experiment showed that tetrachlorantraniliprole mainly stayed in the aqueous phase after entering the water body. Although the chemical dissipated quickly in the aqueous phase (half-life of 0.79-1.5 days), it showed toxic effects on zooplankton communities. Cladocerans, represented by Simocephalus vetulus, were most sensitive to tetrachlorantraniliprole stress. Significant short-term toxicity to cladocerans occurred in all treatment groups, but all recovered within 8 weeks except for the highest concentration group (30.0 µg /L). Based on the ecological recovery results, 7.74 µg tetrachlorantraniliprole/L (nominal concentration, 10.0 µg /L) is suggested to be the no-observed-ecological-adverse-effect concentration (NOEAEC) for the zooplankton community. When this NOEAEC was compared with predicted environmental concentrations (PECs; the PECs in natural ponds simulated by the TOPRICE model for 148 application scheme combinations in major rice-producing areas), a relatively high risk of applying tetrachlorantraniliprole during the rice tillering stage was found. The present study makes a positive contribution to the hypothesis that the current Tier 1 approaches for global acute risk assessment have a sufficient protective effect for assessing the risk of tetrachlorantraniliprole to aquatic organisms. Also, the present results should help us to gain a fuller understanding of the ecological risk of diamide insecticides in aquatic ecosystems and their rational application schemes. Environ Toxicol Chem 2024;43:429-439. © 2023 SETAC.

Swimming behavior and energy metabolism of the calanoid copepod invader Sinodiaptomus sarsi

The appearance of invasive species threatens the integrity of aquatic ecosystems. Much is known about dispersal and introduction mechanisms while little is known on the biological properties of invasive species, such as behavior and energy efficiency, allowing them to successfully colonize new environments and compete with native species. This study examines the functional features of the Asian invasive copepod Sinodiaptomus sarsi (Rylov, 1923) that has invaded Europe since 2016. We focused on the energy metabolism and kinematic parameters of the main swimming types (i.e., gliding, hovering, small relocation jumps, and the escape reaction) of females and males of S. sarsi. Based on the above parameters, the mechanical energy for swimming and the respiration energy needed for movement were calculated. Females and males spend up to 95% of time hovering and slowly gliding at a speed of up to 0.5 cm s-1. During the remaining time, the average swimming speed was 8 cm s-1 by small jumps. In contrast, the average speed was 42 cm s-1 during escape swimming. Non-ovigerous females moved faster than ovigerous females during all relocation swimming types except for upward gliding. While performing small jumps with a frequency of 0.79 Hz, the respiration rate of active non-ovigerous females (0.32 ± 0.03 µg O2 ind-1 h-1) was 2.1 times higher than that of anesthetized individuals. The respiration energy associated with movement was 2.6 * 10-3 J h-1, while the total mechanical energy was only 4.2% of this value. The low energy cost of feeding along with the high speed of locomotion may explain the success of this Asian invader in European waters.

SYP-3343 drives abnormal vascularization in zebrafish through regulating endothelial cell behavior

SYP-3343 is a novel strobilurin fungicide with excellent and broad-spectrum antifungal activity, and its potential toxicity raises public health concerns. However, the vascular toxicity of SYP-3343 to zebrafish embryos is still not well understood. In the present study, we investigated the effects of SYP-3343 on vascular growth and its potential mechanism of action. SYP-3343 inhibited zebrafish endothelial cell (zEC) migration, altered nuclear morphology, and triggered abnormal vasculogenesis and zEC sprouting angiogenesis, resulting in angiodysplasia. RNA sequencing showed that SYP-3343 exposure altered the transcriptional levels of vascular development-related biological processes in zebrafish embryos including angiogenesis, sprouting angiogenesis, blood vessel morphogenesis, blood vessel development, and vasculature development. Whereas, the addition of NAC exerted an improvement effect on zebrafish vascular defects owing to SYP-3343 exposure. Additionally, SYP-3343 altered cell cytoskeleton and morphology, obstructed migration and viability, disrupted cell cycle progression, and depolarized mitochondrial membrane potential, as well as promoted apoptosis and reactive oxygen species (ROS) in HUVEC. SYP-3343 also caused an imbalance of the oxidation and antioxidant systems and irritated the alterations in the cell cycle- and apoptosis-related genes in HUVECs. Collectively, SYP-3343 has high cytotoxicity, possibly by up-regulating p53 and caspase3 expressions and bax/bcl-2 ratio via ROS, leading to malformed vascular development.

Adverse effects of SYP-3343 on zebrafish development via ROS-mediated mitochondrial dysfunction

As a newly-invented and highly-efficiency strobilurin fungicide, pyraoxystrobin (SYP-3343) has been recognized as a highly poisonous toxin for a variety of aquatic organisms. Nevertheless, the developmental toxicity and potential mechanism of SYP-3343 have not been well-documented. The results showed that SYP-3343 was relatively stable and maintained within the range of 20 % in 24 h, and the LC₅₀ value to embryos at 72 hpf was 17.13 μg/L. The zebrafish embryotoxicity induced by 1, 2, 4, and 8 μg/L SYP-3343 is demonstrated by repressive embryo incubation, enhancive mortality rate, abnormal heart rate, malformed morphological characteristic, and impaired spontaneous coiling, indicating SYP-3343 mostly exerted its toxicity in a dose- and time-dependent manner. Besides SYP-3343 was critically involved in regulating cell cycle, mitochondrial membrane potential, and reactive oxygen species production as well as zebrafish primary cells apoptosis, which can be mitigated using antioxidant N-acetyl-L-cysteine. A significant change occurred in total protein content, the biochemical indices, and antioxidant capacities owing to SYP-3343 exposure. Additionally, SYP-3343 altered the mRNA levels of heart development-, mitochondrial function-, and apoptosis-related genes in zebrafish embryos. These results indicated that SYP-3343 induced apoptosis accompanying reactive oxygen species-initiated mitochondrial dysfunction in zebrafish embryos.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

Bioaccumulation of Pyraoxystrobin and Its Predictive Evaluation in Zebrafish

This paper aims to understand the bioaccumulation of pyraoxystrobin in fish. Using a flow-through bioconcentration method, the bioconcentration factor (BCF) and clearance rate of pyraoxystrobin in zebrafish were measured. The measured BCF values were then compared to those estimated from three commonly used predication models. At the exposure concentrations of 0.1 μg/L and 1.0 μg/L, the maximum BCF values for pyraoxystrobin in fish were 820.8 and 265.9, and the absorption rate constants (K₁) were 391.0 d⁻¹ and 153.2 d⁻¹, respectively. The maximum enrichment occurred at 12 d of exposure. At the two test concentrations, the clearance rate constant (K₂) in zebrafish was 0.5795 and 0.4721, and the half-life (t_1/2) was 3.84 d and 3.33 d, respectively. The measured BCF values were close to those estimated from bioconcentration predication models.

A comprehensive review of strobilurin fungicide toxicity in aquatic species: Emphasis on mode of action from the zebrafish model

Strobilurins are popular fungicides used in agriculture on a global scale. Due to their widespread use as agrochemicals, they can enter aquatic environments at concentrations that can elicit adverse effects in organisms. This review synthesizes the current state of knowledge regarding the toxic effects of strobilurin fungicides on aquatic species, including algal species, Daphnia magna, and fish species, to determine risk to aquatic organisms and ecosystems. Data show that the toxicities of strobilurins vary widely across aquatic species. Strobilurins bind cytochrome bc1 in mitochondrial complex III in fungi, and as such, research in aquatic species has focused on mitochondria-related endpoints following exposures to strobilurins. In fish, studies into the activities of mitochondrial complexes and the expression of genes involved in the electron transfer chain have been conducted, converging on the theme that mitochondrial complexes and their enzymes are impaired by strobilurins. In general, the order of toxicity of strobilurins for fish species are pyraoxystrobin > pyraclostrobin ≈ trifloxystrobin > picoxystrobin > kresoxim-methyl > fluoxastrobin > azoxystrobin. In addition to mitochondrial toxicity, studies also report genotoxicity, immunotoxicity, cardiotoxicity, neurotoxicity, and endocrine disruption, and each of these events can potentially impact whole organism-level processes such as development, reproduction, and behavior. Screening data from the US Environmental Protection Agency ToxCast database supports the hypothesis that these fungicides may act as endocrine disruptors, and high throughput data suggest estrogen receptor alpha and thyroid hormone receptor beta can be activated by some strobilurins. It is recommended that studies investigate the potential for endocrine disruption by strobilurins more thoroughly in aquatic species. Based on molecular, physiological, and developmental outcomes, a proposed adverse outcome pathway is presented with complex III inhibition in the electron transfer chain as a molecular initiating event. This review comprehensively addresses sub-lethal toxicity mechanisms of strobilurin fungicides, important as the detection of strobilurins in aquatic environments suggests exposure risks in wildlife.

Periodic density as an endpoint of customized plankton community responses to petroleum hydrocarbons: A level of toxic effect should be matched with a suitable time scale

As an endpoint of community response to contaminants, average periodic density of populations (APDP) has been introduced to model species interactions in a community with 4 planktonic species. An ecological model for the community was developed by means of interspecific relationship including competition and predation to calculate the APDP. As a case study, we reported here the ecotoxicological effects of petroleum hydrocarbons (PHC) collected from Bohai oil field on densities of two algae, Platymonas subcordiformis and Isochrysis galbana, a rotifer, Brachionus plicatilis, and of a cladocera, Penilia avirostris, in single species and a microcosm experiment. Time scales expressing toxic effect increased with increasing levels of toxic effect from molecule to community. Remarkable periodic changes in densities were found during the tests in microcosm experiment, revealing a strong species reaction. The minimum time scale characterizing toxic effect at a community level should be the common cycle of population densities of the microcosm. In addition, the cycles of plankton densities shortened in general with increasing PHC, showing an evident toxic effect on the microcosm. Using APDP as the endpoint, a threshold concentration for the modeled microcosm was calculated to be 0.404 mg-PHC L^-1. The APDP was found to be more sensitive and reliable than the standing crops of populations as the endpoint. This indicated that the APDP, an endpoint at the community level, could be quantitatively related to the endpoints at the population level, and led to the quantitative concentration-toxic effect relationship at the community level.

Evaluation of pyraoxystrobin bioconcentration in zebrafish (Danio rerio) using modified QuEChERS extraction

Pyraoxystrobin is a novel strobilurin fungicide that is widely used on many crops. The high log K_ow of pyraoxystrobin implies that it tends to accumulate in aquatic organisms. This study optimized the sorbents of QuEChERS (quick, easy, cheap, effective, rugged, and safe) using ¹³C-labelled pyraoxystrobin as the internal standard (IS). It has been established a QuEChERS-LC-MS/MS IS method to study the bioconcentration and elimination of pyraoxystrobin in zebrafish (Danio rerio). The results indicated that the method had satisfactory linearity between 0.234 and 15 μg L^-1 (R² = 0.9996). The limits of detection (LOD) and quantification (LOQ) for pyraoxystrobin were 0.01 and 0.03 μg L^-1, respectively. The LOQs of the method for water and zebrafish were 0.05 μg L^-1 and 0.01 mg/kg, respectively. The mean recovery of pyraoxystrobin in zebrafish and water at fortification levels of 0.01-0.3 mg kg^-1 and 0.05-1.5 μg L^-1 ranged from 98.31 to 105.61% and 101.87 to 108.48%, respectively, with a % RSD (relative standard deviation) of 0.94-3.57%. The bioconcentration has been evaluated. The bioconcentration factors for pyraoxystrobin in zebrafish were 1,792 and 3,505 after exposure to 0.5 μg L^-1 for 168 h and 0.05 μg L^-1 for 216 h, respectively. The half-life of pyraoxystrobin in zebrafish was 9.0-9.5 d.

Ecological risk assessment of petroleum hydrocarbons on aquatic organisms based on multisource data

From the perspective of ecological risk, this study uses a multisource data method to search for global data, uses the acute and chronic ratio method to process the data, uses the species-sensitive distribution method to evaluate the ecological risk that petroleum hydrocarbons pose to aquatic organisms, and evaluates the ecological risk of the water environment in five Chinese water bodies. The results are as follows. First, in an aquatic ecosystem, the toxicological effects of petroleum hydrocarbons were found to be more obvious on consumers, and the sensitivity of fish was found to be higher than that of crustaceans. Second, the acutely lethal effects of petroleum hydrocarbons, fluorene, and benzo [a] pyrene on aquatic ecosystems were fitted by using the documentary method of multisource data collection and a Log-logistic curve. Third, in the case study evaluation of five Chinese water bodies, the ecological risks of polycyclic aromatic hydrocarbons were ranked (from low to high) as fluorene < benzo [a] pyrene. The ecological risk values of benzo [a] pyrene were all greater than 1. These risks should not be underestimated, and prevention and control work should be performed.

Fungicides: An Overlooked Pesticide Class?

Fungicides are indispensable to global food security and their use is forecasted to intensify. Fungicides can reach aquatic ecosystems and occur in surface water bodies in agricultural catchments throughout the entire growing season due to their frequent, prophylactic application. However, in comparison to herbicides and insecticides, the exposure to and effects of fungicides have received less attention. We provide an overview of the risk of fungicides to aquatic ecosystems covering fungicide exposure (i.e., environmental fate, exposure modeling, and mitigation measures) as well as direct and indirect effects of fungicides on microorganisms, macrophytes, invertebrates, and vertebrates. We show that fungicides occur widely in aquatic systems, that the accuracy of predicted environmental concentrations is debatable, and that fungicide exposure can be effectively mitigated. We additionally demonstrate that fungicides can be highly toxic to a broad range of organisms and can pose a risk to aquatic biota. Finally, we outline central research gaps that currently challenge our ability to predict fungicide exposure and effects, promising research avenues, and shortcomings of the current environmental risk assessment for fungicides.

Acute toxicity and associated mechanisms of four strobilurins in algae

Strobilurins have been reported highly toxic to non-target aquatic organisms but few illustrated how they cause toxic effects on algae. This study investigated the acute toxicity of Kresoxim-methy (KRE), Pyraclostrobin (PYR), Trifloxystrobin (TRI) and Picoxystrobin (PIC) on two algae and their toxicity mechanisms. Four strobilurins showed lower toxic effects on Chlorella pyrenoidsa but higher on Chlorella vulgaris. bc1 complex activities in C. vulgaris were significantly inhibited by all strobilurins, suggesting bc 1 complex might be the target of strobilurin toxicity in algae. Moreover, SOD, CAT and POD activities were significantly up-regulated by all doses of KRE, PYR and PIC. In contrast, low concentrations of TRI stimulated SOD and POD activities but highest concentration significantly inhibited those activities. Comet assays showed damaged DNA in C. vulgaris by four strobulirins, suggesting their potential genotoxic threats to algae. The results illustrated acute toxicity by strobulirins on algae and their possible toxicity mechanisms.

Developmental toxicity and potential mechanisms of pyraoxystrobin to zebrafish (Danio rerio)

As a newly developed, highly efficient strobilurin fungicide, pyraoxystrobin has been reported to be highly toxic to some aquatic organisms. However, the toxicity of pyraoxystrobin to different life stages of fish and the potential underlying mechanisms are still unknown. Hence, in the present study, the acute toxicity of pyraoxystrobin to different life stages of zebrafish (embryo, larva, and adult) was assessed. The developmental toxicity of pyraoxystrobin to zebrafish embryos and its effects on gene transcription in the embryo were also investigated. The results showed that the 96-h LC₅₀ values of pyraoxystrobin to embryos [2h post-fertilization (hpf)], 12h post-hatching (hph) larvae (84 hpf), 72 hph larvae (144 hpf), and adult zebrafish were 4.099, 1.069, 3.236, and 5.970µg/L, respectively. This suggests that pyraoxystrobin has very high toxicity to different life stages of zebrafish, while the newly hatched larvae constitute the most sensitive period of zebrafish to pyraoxystrobin. Decreased heart rate, hatching inhibition, growth regression, and morphological deformities were observed in zebrafish embryos after acute exposure to different concentrations of pyraoxystrobin. The rate of malformation increased in a time- and concentration-dependent manner in embryos, and the most pronounced abnormality was pericardial edema and yolk sac edema. Pyraoxystrobin (2 and 4μg/L) significantly altered the mRNA levels of genes related to mitochondrial respiratory chain and ATP synthesis (NDI, uqcrc, and ATPo6), oxidative stress (Mn-Sod, Cat, and Gpx), apoptosis (p53, Bcl2, Bax, and Cas3), and immune system (TNFα, IFN, and IL-1b) in zebrafish embryos. This result indicates that the alteration of these genes is a potential mechanism underlying the toxic effects of pyraoxystrobin on zebrafish.

Fate Characterization of Benzene Kresoxim-Methyl (a Strobilurin Fungicide) in Different Aerobic Soils

Benzene kresoxim-methyl (BKM) is a promising broad-spectrum strobilurin fungicide widely used to control fungal pathogens in crops. However, information on its environmental fate is limited. To broaden our understanding of this fungicide's kinetic fate in aerobic soils, we labeled BKM with C on its benzoate ring and used ultralow-level liquid scintillation counting coupled with high-performance liquid chromatography analysis. Results show that degradation, mineralization, and bound residue (BR) formation of BKM was controlled by soil type and microbial community composition. Degradation of BKM followed first-order dynamics, and the half-lives () were 51.7, 30.8, and 26.8 d for clay, loamy, and saline soils, respectively. After 100 d, about 0.13, 4.35, and 5.94% of the initial C-BKM was mineralized, and 14.43, 19.90, and 28.81% was formed as BRs in the clay, loamy, and saline soils, respectively. About 60 to 85% of the C-BKM residue in soil was extractable; of this fraction, 30 to 50% was composed of incomplete degradation intermediates. Up to 40% of extractable C-BKM in soil was readily available. Our results suggest that BKM and its incomplete intermediates had a relatively long persistence in soil, which may lead to exposure for nontarget organisms. Soil microbes may play a dominant role in controlling the fate of BKM in soil as sterilization sharply decreased its mineralization rate from 4.35 to 0.03%, increased from 30.8 to 85.6 d, and decreased the BR fraction from 19.90 to 3.25%.