Next Article in Journal
Continuous Monitoring of Advanced Hemodynamic Parameters during Hemodialysis Demonstrated Early Variations in Patients Experiencing Intradialytic Hypotension
Next Article in Special Issue
Embryonic Zebrafish as a Model for Investigating the Interaction between Environmental Pollutants and Neurodegenerative Disorders
Previous Article in Journal
Machine Learning Model in Obesity to Predict Weight Loss One Year after Bariatric Surgery: A Pilot Study
Previous Article in Special Issue
Essential Oils Produce Developmental Toxicity in Zebrafish Embryos and Cause Behavior Changes in Zebrafish Larvae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Behavioral Effects of the Mixture and the Single Compounds Carbendazim, Fipronil, and Sulfentrazone on Zebrafish (Danio rerio) Larvae

by
Samara da Silva Gomes
1,
Jadson Freitas da Silva
1,
Renata Meireles Oliveira Padilha
1,
João Victor Alves de Vasconcelos
2,
Luís Gomes de Negreiros Neto
2,
James A. Marrs
3 and
Pabyton Gonçalves Cadena
1,*
1
Department of Morphology and Animal Physiology, Universidade Federal Rural de Pernambuco, Av. Dom Manoel de Medeiros s/n, Dois Irmãos, Recife 52171-900, PE, Brazil
2
Department of Physics, Universidade Federal Rural de Pernambuco, Av. Dom Manoel de Medeiros s/n, Dois Irmãos, Recife 52171-900, PE, Brazil
3
Department of Biology, Indiana University Purdue University Indianapolis, 723 West Michigan, Indianapolis, IN 46202, USA
*
Author to whom correspondence should be addressed.
Biomedicines 2024, 12(6), 1176; https://doi.org/10.3390/biomedicines12061176
Submission received: 18 April 2024 / Revised: 13 May 2024 / Accepted: 22 May 2024 / Published: 25 May 2024
(This article belongs to the Special Issue Zebrafish Models for Development and Disease 4.0)

Abstract

:
Pesticides are often detected in freshwater, but their impact on the aquatic environment is commonly studied based on single compounds, underestimating the potential additive effects of these mixtures. Even at low concentrations, pesticides can negatively affect organisms, altering important behaviors that can have repercussions at the population level. This study used a multi-behavioral approach to evaluate the effects of zebrafish larvae exposure to carbendazim (C), fipronil (F), and sulfentrazone (S), individually and mixed. Five behavioral tests, thigmotaxis, touch sensitivity, optomotor response, bouncing ball test, and larval exploratory behavior, were performed to assess potential effects on anxiety, fear, and spatial and social interaction. Significant changes were observed in the performance of larvae exposed to all compounds and their mixtures. Among the single pesticides, exposure to S produced the most behavioral alterations, followed by F and C, respectively. A synergistic effect between the compounds was observed in the C + F group, which showed more behavioral effects than the groups exposed to pesticides individually. The use of behavioral tests to evaluate pesticide mixtures is important to standardize methods and associate behavioral changes with ecologically relevant events, thus creating a more realistic scenario for investigating the potential environmental impacts of these compounds.

1. Introduction

Pesticides are a recognized problem that goes beyond human health because they pose serious threats to other vertebrates, invertebrates, and ecosystems in general. The dispersal of these products can reach non-target organisms such as fish, which can suffer physiological changes and, in some cases, can even cause death. In addition, changes in the organism’s behavior, in turn, can reduce the fitness of an individual, leading to population decline and serious effects on the ecosystem [1,2]. The application of these products does not respect natural barriers such as relief and riparian forests, and it is almost impossible to control their course in the field [3].
Recently, a series of studies showed that pesticides including carbendazim, fipronil, and sulfentrazone were detected in aquatic ecosystems and even in food [4,5], indicating that assessing the risks of these pesticides and their mixtures is necessary and important because they are used in common crops [6,7,8]. It is known that the fungicide carbendazim (C9H9N3O2) can promote neural excitability (0.57–0.64 mg/L), promoting a hyperactive state in zebrafish larvae, leading to disorganized swimming patterns [9]. Just as the insecticide fipronil (C12H4Cl2F6N4OS) can cause hyperpolarization of neuronal membranes, causing convulsions and muscle spasms (0.50–2.00 mg/L). It can also damage the nervous system, including structures such as the optic thalamus (optic tectum) of aquatic vertebrates [10,11]. The neurotoxic effects of the herbicide sulfentrazone (C11H10Cl2F2N4O3S) are not well understood. However, it acts on mitochondrial complex IV and alters cellular energy supply, acting like an endocrine disruptor and affecting zebrafish heart development in the early stages (0.01–0.40 mg/L) [12].
The assessment of these chemicals in relation to the aquatic environment is mainly studied by determining lethal concentrations, which is essential from a regulatory point of view, but reflects only the worst-case scenario [2]. The subtle effects of sublethal exposure are poorly understood, but are important to understand environmental impact [13]. Among sublethal effects, behavior is increasingly being studied. Although the behavioral analysis was initially slow to be integrated into aquatic toxicology, its speed, sensitivity, and environmental relevancy compared to traditional lethality endpoints make it a promising source of information [2,14]. Several studies report behavioral toxic effects at concentrations that are orders of magnitude lower than lethal concentrations of the chemicals [2,15]. Pesticides are commonly found in the aquatic environment at low concentrations, which accumulate, and effects are amplified through the food chain [16].
Zebrafish (Danio rerio) is widely used as an animal model for behavioral studies [14,15], not only because of their morphological and genetic conservation with humans (60 to 80% homology) [17,18], but also behavioral similarities [2,19,20], exhibiting a wide range of complex behaviors including social interactions, anxiety, learning, memory, and avoidance behaviors that may be useful for modeling neurological and psychiatric diseases [15].
There is growing interest in the behavior of zebrafish larvae for large-scale testing [21]. A single reproductive cycle can produce hundreds of embryos that rapidly develop, with organs such as the brain, heart, liver, pancreas, kidneys, bones, muscles, and sensory systems maturing within 5 dpf (days post-fertilization) [22]. Behaviors, such as swimming, hunting, and fleeing, are exhibited during the first week of development [23,24]. In addition, neurotoxic agents can be easily studied by exposure to their water during the early developmental period [25].
Despite the obvious difference between zebrafish and humans, the behavioral paradigms of mammals and zebrafish are closely comparable, suggesting an evolutionarily conserved nature in many behaviors across species [19,23]. Even with neuroanatomical differences between mammals and teleost, research shows homologous functions in several areas of the zebrafish brain [26]. Zebrafish larvae exhibit a robust cognitive repertoire, such as memory, social and spatial learning, anxious behaviors, stress, and fear, just as in other vertebrates, and the same neurotransmitters and neuroendocrine systems are present [15,27].
Many tests are used to investigate toxicological effects on behavior, such as basic motor responses, sensorimotor responses, and/or learning and memory [13,28]. However, due to the high complexity of behavioral responses, there is currently no standardized method or assay, which leads to difficulties when comparing behavioral experiments and their results [29]. Therefore, the strategy with multi-behavioral testing is advisable to determine the interactions between nervous system structures [30], since chemicals do not always produce the same behavioral phenotype [31].
Therefore, the present study tests the pesticides carbendazim, fipronil, and sulfentrazone in zebrafish larvae, analyzing their effects through a multi-behavioral analysis, since this can reveal defects not only in neuronal cells, but also in neural functions [21]. Previous studies suggest that these pesticides alone and in combination cause neurotoxic effects in zebrafish [9,10,12,32,33].

2. Materials and Methods

2.1. Reagents and Solutions

The pesticides used to prepare the test solutions were carbendazim (C) (lot # 002-18-54392, CAS: 10605-21-7, 50% (w/v)) formulated by Adama Brasil S/A (Londrina, Brazil), fipronil (F) (lot # 001/19, CAS: 120068-37-3, 2.5 (w/w)) formulated by Rogama Indústria e Comércio LTDA (Pindamonhangaba, Brazil) and sulfentrazone (S) (lot #1041-19-13767, CAS: 122836-35-5, 50% (w/v)) formulated by FMC Química do Brasil Ltd.a (Paulínia, Brazil), were purchased from commercial suppliers. The three pesticides tested were initially diluted in dimethyl sulfoxide (DMSO) (lot # 85713, CAS: 67-68-5, ≥99.9% purity, Dinâmica Química Contemporânea Ltd.a., Indaiatuba, Brazil), and these solutions were subsequently diluted in embryo medium (8.00 mM NaCl, 0.40 mM KCl, 0.60 mM KH2PO4, 0.35 mM Na2HPO4, 0.72 mM CaCl2, 1.23 mM MgSO4, and 0.35 mM NaHCO3 at pH 7.2) [34] to obtain the final concentrations (Table 1). The final nominal concentration of DMSO was less than 0.01% (v/v), and this is 50 times lower than that reported that affected behavioral parameters in zebrafish [35]. Pesticide concentrations were based on values equal to or lower than the maximum residue limits (MRL) of these pesticides in foods allowed by the Brazilian Health Regulatory Agency—ANVISA [6,7,8]. The pesticide residues were treated by advanced oxidation process (AOP) in a reactor using hydrogen peroxide/ultraviolet irradiation (16 w) before final discard.

2.2. Zebrafish Culture and Embryo Production

All protocols involving the animals were approved by the Ethics Committee for the Use of Animals (License No. 7373131021). Adult wild-type fish (1 year) were reared and housed at the Laboratório de Ecofisiologia e Comportamento Animal—LECA vivarium, Universidade Federal Rural de Pernambuco—UFRPE. Zebrafish were quarantined in 80 L aquariums to detect or confirm the absence of pathogens or diseases. They were housed under the following laboratory conditions, artificial aeration of 11 mg/L DO, a temperature of 25 ± 1 °C, pH 7.5 ± 0.5, and a 14/10 h (light/dark) cycle. The water was partially renewed once a week. Abiotic parameters such as dissolved oxygen, ammonia, nitrite, and nitrate were also measured and maintained within ideal ranges [36]. The animals were fed three times a day with 2× Fort Color® fish feed (30% crude protein), and 1× with live brine shrimp nauplii (Artemia ssp). To obtain embryos, zebrafish males and females, in the ratio of 2:1 [34], were placed in the spawning tanks (Alesco® Zebclean, Monte Mor, Brazil) for reproduction. After 30 min from the start of spawning the eggs were collected, the unfertilized eggs were removed, and the fertilized eggs (normal blastula development) [36] were washed and transferred to the exposure chambers.

2.3. Chemical Exposure

Toxicity tests were carried out following OECD 210 guidelines [36] with modifications. Briefly, embryos at 2 hpf (hours post-fertilization) were randomly added to sterile polystyrene chambers and exposed to 50 mL of the test solutions (Table 1). Animals were exposed to 144 hpf (or 6 dpf). The solutions were renewed daily, and an embryo medium without pesticides was used as a control. We chose not to feed the larvae, as they live up to 7 dpf exclusively on the nutrients provided by the yolk. According to Clift et al. [37], larvae fed at 6 and 7 dpf showed changes in swimming and resting speed compared to their non-fed congeners.

2.4. Behavioral Tests

Behavioral tests were conducted between 11 am–12 noon [38] using a filming system equipped with a camera (Canon 6D Mark and lens EF 75–300 mn f/4–5.6 III), LED board, and/or LCD screen to playback and record the videos of larvae. In total, 15 larvae were used for each experimental group, and all larvae were tested simultaneously and randomly and went through the same sequence of behavioral tests: thigmotaxis, touch sensitivity, optomotor response, exploratory behavior of larvae, and bouncing balls, respectively. Each test was repeated seven times, one replicate per week (n = 105 per group with 8 experimental groups ≈ 840 larvae). Behavioral tests were conducted only with larvae without visible morphological teratogenic effects [39,40]. Based on the previous screening test, pesticide concentrations used in this study did not produce mortality or significant morphological teratogenic effects. Initially, the larvae were distributed in 48-well plates, and subsequently the same larvae were transferred to 6-well plates, where they were organized into groups of 5 larvae per well for conducting a bouncing ball test. The behavioral evaluation of the larvae was performed through subsequent analysis of the videos recorded during the tests described below and performed in the sequence on the same day, respectively (Figure 1).

2.4.1. Thigmotaxis

The thigmotaxis test was used to analyze the propensity of the larvae to approach the walls of the plates to evaluate behavior similar to anxiety in zebrafish [23,39] because when the larva enters an unfamiliar environment, they tend to swim against the walls of the well to explore the new environment [40]. Before behavioral tests, the larvae underwent 15 min of acclimatization in 48-well plates in an environment with natural lighting, and each well was filled with 1 mL of the same test solution to which the animals were previously exposed. The response to the thigmotaxis test was recorded by photography so that all larvae could be evaluated in the same time interval. This response was measured according to the location of the larvae in the well, whether they were at the edge (positive) or in the center of the well (negative).

2.4.2. Touch Sensitivity

The touch sensitivity test was used to evaluate the larval response to mechanical stimuli. The larvae were gently touched in the posterior portion of the tail, and the response was recorded if the larva exhibited escape behavior. An escape behavior response was determined when larvae displayed swimming behaviors (positive) [39], and no response when larvae remained motionless after the touching (negative).

2.4.3. Exploratory Behavior of Larvae

The exploratory behavior of the larvae was evaluated during spontaneous swimming (without stimuli), where the larvae tend to swim at the limits of the wells (close to the wall), a tendency that is explained by the initial exploratory behavior and the avoidance of the central area of the well, where there is less protection (thigmotaxis). This assay was based on the protocol by Altenhofen et al. [41] and Pérez-Escudero et al. [42], with modifications. After 15 min of acclimatization, the exploratory activity of the larvae was observed in the 48-well plates for 5 min with artificial lighting (LED plate) without stimuli. However, the first and last minute of recording were excluded from the analysis, as it was observed that in the first minute, the larvae did not show swimming behaviors due to exposure to artificial light from the LED board, and in the last minute, the animals reduced their display of swimming behaviors, as was also observed by [42]. During the analysis, the animals that remained motionless during the recording were counted, and those that showed some swimming movement had their average speed (cm/s) and total trajectory covered (cm), as measured using the IdTracker software v 8.3 [43]. These were the parameters considered for exploring the new environment in this assay.

2.4.4. Optomotor Response

The optomotor response is used to measure the innate visual responses of larvae in which the individual swims in the direction of the optical flow and mimics naturalistic behavior, aiming to stabilize itself according to the flow of the video lines. The methodology was based on Cadena et al. [39] and Brastrom et al. [44], with modifications. The plates with animals were placed on a 21.5 inch LCD screen (Dell E2211H) in a room without external sounds and lights. The fish were acclimated for 15 min on a dark screen, and then the animation was started (video as black and white lines moving from right to left for 30 s and then left to right for another 30 s). Subsequently, the video recording was analyzed, and larval positioning was classified into two types with respect to the flow of the lines: low alignment (70% of the well), and high alignment (30% of the well), according to Cadena et al. [39].

2.4.5. Bouncing Ball Test

The bouncing ball test is used to observe the spatial interaction and social and escape responses of larvae to visual stimuli [21,45,46]. Larvae were transferred to 6-well plates (five larvae per well) and placed on an LCD screen (Dell E2211H), as described above. An animation with visual stimuli created in Microsoft PowerPoint v 2404 (Office 365) was played. The animation had red balls (diameter 1.35 cm) that moved from left to right on a square with a 2 cm trajectory in the bottom half of the well area (stimulus area) for 5 min after a 15 min acclimation on a dark screen. The percentage of larvae in the non-stimulated area (upper part of the well) during the 5 min session was taken to be indicative of the cognitive ability of the group to exhibit escape behavior against the red ball. The number of larvae in the same quadrant of the well was used to measure the percentage of grouping (positive, if three or more larvae were grouped) and social interaction among the larvae in the well (positive, if larvae show escape behavior from red ball). A time-point sampling method was used to collect the data (one analysis every 30 s) [45,46].

2.5. Statistical Analysis

Statistical analyses were performed using Origin Pro Academic 2015 (Origin Lab. Northampton, MA, USA). The normality was determined with the Shapiro–Wilk test with p < 0.05, and Levene’s test for homogeneity of variance was carried out on each variable. All data were presented as mean ± SD. For behavioral tests of thigmotaxis, sensitivity to touch, exploratory behavior, and bouncing balls, we used pooled samples from several animals. A one-way ANOVA was used for data that presented a normal distribution, followed by the Tukey test with p < 0.05 to determine the significant difference between groups. To evaluate the interactions of pesticides in the trajectory followed during the exploratory activity test, the simplex centroid design method was carried out using Statistica 14 (TIBCO, Palo Alto, CA, USA) according to Cadena et al. [40]. The results of the optomotor response test for each animal were analyzed using Bowker’s non-parametric test with p < 0.05 [44].

2.6. Cluster Analysis

The tests were evaluated using one-way ANOVA for the seven behaviors observed. For the results obtained, a number between 1 and 3 was assigned based on their statistical p-values, according to Richendrfer and Creton [47], with some modifications. Any behavior that was statistically different from the control was given the numbers 1 or 2 (1 for p < 0.01 or 2 for p < 0.05). Behaviors that were not statistically different from the control were given the number 3 (p > 0.05). The numbers assigned to the behavior scores were imported into Origin Pro Academic 2015 software (Origin Lab. Northampton, MA, USA). Colors were assigned to the numbers to illustrate; the green color corresponds to numbers 1 and 2 (p < 0.01 and p < 0.05) and black corresponds to 3 (p > 0.05). The results of the Optomotor Activity Test were not included in this analysis because the data were analyzed using a non-parametric test (Bowker’s test), and thus there was no p-value against the control to be compared.

3. Results

3.1. Thigmotaxis and Touch Sensitivity Tests

Larvae that were exposed to single pesticides and their mixtures. Thigmotaxis (Figure 2a) and touch sensitivity (Figure 2b) tests showed no significant difference from the control group.

3.2. Exploratory Behavior of Larvae

In this test, first, the animals that remained immobile (Figure 3a) during all the time intervals of the test were separated and used to determine the percentage of the immobile animals. Group F was the one that showed a significant difference, with a higher percentage of immobile animals. Later, with the animals that showed some swimming movement, the average speed (Figure 3b) and the total distance traveled (Figure 3c) were evaluated. Groups S, C + F, and F + S showed a significant difference in mean velocity, and these same groups showed a significant difference in the total distance traveled compared to the control. Results showed that groups with the lowest velocity (Figure 3b) were also groups of larvae with the shortest trajectory (Figure 3c). It was also observed that the groups C + F and F + S showed a synergistic effect (Figure 3d), as they showed a reduction compared to groups exposed to C and F individually. It was also possible to qualitatively evaluate the trajectory (Figure 4) of the larvae in each group, and we saw that even the groups that did not show a significant reduction compared to the control in their trajectory showed an abnormal trajectory during the tests.

3.3. Optomotor Response Test

All larvae exhibit optomotor responses in the first period of analysis with right-to-left line movement (Figure 5). However, it was observed in the optomotor response test that the larvae exposed to the pesticides F and S alone, the mixtures C + F and C + S showed no significant difference in the second period of analysis with left to right line movement compared to the control. The results demonstrated that these larvae exhibited low alignment with the direction of movement of the stripes and did not respond normally to the animation, indicating that these pesticides caused a deficiency in the optomotor system of zebrafish larvae. Interestingly, it was also observed that the groups exposed to the mixtures C + F and C + S showed no significant difference, suggesting that when the pesticides F and S are mixed with C, they have a synergistic effect, resulting in a higher percentage of larvae with low alignment compared to when the larvae were exposed to the three pesticides individually. However, when F and S were mixed (F + S), an antagonistic effect was observed, reducing the percentage of larvae with low alignment compared to that observed with the two single pesticides.

3.4. Bouncing Balls Test

In the bouncing balls test, two behaviors were observed: escape from the visual stimulus (Figure 6a) and the grouping (Figure 6b) of the larvae. All exposed groups showed a significant difference compared to the control in the two analyzed behaviors. About escape from the visual stimulus, larvae exposed to the single pesticides and their mixtures stayed significantly longer in the lower part of the well (area with stimulus) and did not show evasive behavior relative to the red ball, as the control did.

3.5. Cluster Analysis

The analysis of the behavioral results clustering by exposure group is described in Figure 7. Cluster analysis identifies the groups that presented the greatest number of different behaviors from the control. Based on this, we observed that among the isolated compounds, sulfentrazone was the group that presented the most differences from the control, followed by fipronil, and then carbendazim. We also saw that the binary mixture C + F showed more differences than the isolated compounds, even though it was 50% of the isolated concentration, suggesting that these mixed compounds potentiated the toxic effects in zebrafish larvae. This was similar to the group that was exposed to the binary mixture of F + S.

4. Discussion

Our study used a multi-behavioral approach to assess how single and mixed pesticides affected zebrafish larvae. Interestingly, this approach was very sensitive as a tool, and we recommend its inclusion in the toxicological tests. In addition, the analysis of different behavioral parameters can help to estimate more accurately the impact of pesticides on zebrafish behavior. The changes in the behavior result from one or a combination of molecular, biochemical, and physiological changes [48].
The group exposed to sulfentrazone was the one that showed the highest number of behavioral changes (Figure 7). Showed significant differences in mean speed (Figure 3b), distance traveled (Figure 3c), and presented an abnormal trajectory (Figure 4). According to Jiang et al. [12], exposure to sulfentrazone (0.01–3.60 mg/L for 30 days) produced adverse effects on calcium (Ca2+) channel regulation in zebrafish larvae, affecting the activities of Ca2+-ATPase, total Na+K+-ATPase, and Ca2+Mg2+-ATPase, as well as induced the transcription of many key calcium manipulation proteins, such as ATPase, Ca2+, transporter subunits Atp2b2, Atp2b3a, and Atp2b3b, and ATPase Na+/K+ transporter subunit Atp1b4. This biochemical process participates in muscle mechanics, and its alteration can compromise the locomotor performance of zebrafish [48]. Because the interaction of actin and myosin depends on the promotion of Ca2+, where Ca2+ is exported with high specificity by Ca2+-ATPases in the plasma membrane and the deregulation of Ca2+-ATPase activity, it can cause defects in the locomotor system of zebrafish [49]. We suggest that the alterations described above induced by sulfentrazone in zebrafish may contribute to the locomotor deficit observed in our study. Changes in swimming behavior can have important consequences for feeding behaviors (prey capture) and susceptibility to predation (predator escape), among others [48].
Larvae exposed to fipronil showed a significant increase in immobile behavior during the exploratory activity test (Figure 3a) and showed an abnormal trajectory in the test (Figure 4). Immobility of zebrafish was also observed by Wu et al. [10] when analyzing adult zebrafish exposed to fipronil (0.50–2.00 mg/L). They saw that after 24 h of exposure, there was a significant reduction in swimming and the traveled distance by the animals. The authors attributed these locomotion abnormalities to the neurotoxicity of fipronil in zebrafish brain tissue, possibly due to oxidative neural stress, inflammation, and apoptosis of neurons. This was confirmed with analysis of markers of oxidative stress (SOD2), inflammation (TNF-α), and apoptosis (caspase-3), showing a significant decrease in SOD2 and TNF-α levels and a significant increase in caspase-3 in zebrafish brain tissue. Park et al. [32] observed that fipronil caused developmental delays and motor neurons.
Fish generally adopt self-defense behaviors, like remaining motionless near the bottom or walls [50], which is an effective strategy to avoid capture by predators. However, fipronil resulted in the zebrafish larvae losing their ability to move, making larvae more vulnerable to predators. This group also showed significant differences in the optomotor response tests (Figure 5) and bouncing balls (Figure 6). Fipronil acts directly on γ-aminobutyric acid, preventing the closure of chloride channels (GABA) [51]. Reduced GABA receptor signaling may lead to loss of retinal ganglion cells [52]. Fipronil may impair the visual system.
Among the groups exposed to a single pesticide, carbendazim was the least toxic (Figure 7). The group exposed to carbendazim did not show reductions in average speed and traveled distance, but an abnormal trajectory was observed in the exploratory activity test (Figure 4). This change in swimming pattern was also observed by Zhang et al. [9]. The researchers associated changes in swimming behavior with changes in gene expression, namely the ctsbl gene, which plays a crucial role in antigen processing, regulation of cell death (apoptosis), autophagy, and metabolism, contributing to increased cell death mainly in the nervous system [9].
The carbendazim group also showed significant differences in the optomotor response test (Figure 5) and bouncing ball (Figure 6a,b). This result may be related to the impairment of the visual system in this exposed group. Xenopus laevis (clawed frog) larvae exposed to carbendazim when analyzed histologically by Yoon et al. [53] showed dysplasia and optic edema as the most common malformations. Furthermore, it was also observed that the retinal layers were poorly differentiated and thinner and the neural cells in the larval brain were disaggregated and invaded the cerebral ventricle, demonstrating that carbendazim strongly inhibited the differentiation of neural tissue. The authors suggest that this inhibition may be related to the interruption of tubulin synthesis, weakening cell junctions. Carbendazim interferes with the assembly of tubulin and the formation of microtubules, resulting in various malformations, genetic toxicity, and reproductive problems [54].
Our results showed that larvae exposed to binary mixtures C + F and F + S exhibited synergistic responses in zebrafish (Figure 3d) in the behavioral parameters affected in our work. Although the concentration of the binary mixture was 50% lower than the single pesticide groups, we observed reduced speed (Figure 3b), distance traveled (Figure 3c), and optomotor response (Figure 5). Synergistic toxicity is often found in organisms exposed to pesticide mixtures [55]. Larvae exposed to mixtures may experience oxidative stress, leading to the production of reactive oxygen species (ROS) [56], which could help explain our results. Several studies of pesticide mixtures with different modes of action have pointed to increased ROS production in zebrafish because of the synergy between pesticides [55,56,57]. Excessive ROS production was associated with apoptotic cell death in zebrafish embryos and larvae [58].
The binary and ternary mixtures, C + S, F + S, and C + F + S, also caused significant behavioral changes. We observed different toxicity patterns for mixtures depending on the parameters analyzed, where synergistic effects were observed. This type of response was also observed in other studies, where the effects on the same organism varied depending on the chosen endpoint [59,60]. In this way, we can also see that behavioral tests are useful tools and have the potential to reveal environmental stress and how combinations of effects produce heterogeneous behavioral changes in the environment. Future studies could elucidate the additive and synergistic mechanisms produced by mixed pesticides.

5. Conclusions

This study addressed the impact of pesticides, both individually and in mixtures, on zebrafish larvae using a multi-behavioral approach. This approach proved to be a sensitive and comprehensive tool for analyzing the systemic effects of one or more substances in zebrafish larvae. The results showed that among the individual pesticides, sulfentrazone caused the greatest number of behavioral changes in exposed larvae. Fipronil, on the other hand, caused immobile behavior. Carbendazim had less pronounced effects, but still affected the animals’ escape and grouping behavior. When we analyzed pesticide mixtures, we observed synergistic responses that affected speed, distance covered, and optomotor response. These results highlight the complexity of pesticide mixture toxicity and the sensitivity of behavioral tests, which can be used as initial indicators of environmental stress. Experiments testing pesticide mixtures will help us understand their additive and synergistic environmental exposure effects.

Author Contributions

Conceptualization, S.d.S.G. and P.G.C.; methodology, S.d.S.G., J.A.M. and P.G.C.; software, J.V.A.d.V. and L.G.d.N.N.; validation, S.d.S.G. and P.G.C.; formal analysis, S.d.S.G., J.F.d.S., R.M.O.P. and P.G.C.; investigation, S.d.S.G., J.F.d.S., R.M.O.P. and P.G.C.; resources, S.d.S.G., J.A.M. and P.G.C.; data curation, S.d.S.G., J.A.M. and P.G.C.; writing—original draft preparation, S.d.S.G. and P.G.C.; writing—review and editing, S.d.S.G., J.A.M. and P.G.C.; visualization, S.d.S.G. and P.G.C.; supervision, P.G.C. and J.A.M.; project administration, P.G.C.; acquisition financing, P.G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CNPq and CAPES. The authors would like to thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—CAPES for the scholarship to S.S.G. (88887.528986/2020-00), Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq for the productivity grant to P.G.C. (#306947/2020-0, and #423038/2018-4).

Institutional Review Board Statement

Methods and procedures used in this work are in accordance with guidelines and supervision of the Ethics Committee for the Use of Animals (License nº 7373131021) of Universidade Federal Rural de Pernambuco.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data and material used in this study are available from the corresponding author on request. The data are not publicly available due to complexity and lack of a standard format.

Acknowledgments

The authors would like to thank James Marrs (IUPUI) for helpful discussion, and teaching in the zebrafish behavior, and we thank students of Laboratório de Ecofisiologia e Comportamento Animal—LECA (Brazil) for inspiration.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Peterson, E.K.; Buchwalter, D.B.; Kerby, J.L.; LeFauve, M.K.; Varian-Ramos, C.W.; Swaddle, J.P. Integrative behavioral ecotoxicology: Bringing together fields to establish new insight to behavioral ecology, toxicology, and conservation. Curr. Zool. 2017, 63, 185–194. [Google Scholar] [CrossRef] [PubMed]
  2. Fitzgerald, J.A.; Konemann, S.; Krumpelmann, L.; Zupanic, A.; Vom Berg, C. Approaches to Test the Neurotoxicity of Environmental Contaminants in the Zebrafish Model: From Behavior to Molecular Mechanisms. Environ. Toxicol. Chem. 2021, 40, 989–1006. [Google Scholar] [CrossRef] [PubMed]
  3. Severo, E.S.; Marins, A.T.; Cerezer, C.; Costa, D.; Nunes, M.; Prestes, O.D.; Zanella, R.; Loro, V.L. Ecological risk of pesticide contamination in a Brazilian river located near a rural area: A study of biomarkers using zebrafish embryos. Ecotoxicol. Environ. Saf. 2020, 190, 110071. [Google Scholar] [CrossRef] [PubMed]
  4. Li, M.; Ma, X.; Saleem, M.; Wang, X.; Sun, L.; Yang, Y.; Zhang, Q. Biochemical response, histopathological change and DNA damage in earthworm (Eisenia fetida) exposed to sulfentrazone herbicide. Ecol. Indic. 2020, 115, 106465. [Google Scholar] [CrossRef]
  5. Kronbauer, E.A.; Biondo, E.; Zanetti, C.; Kolchinski, E.M. Agrotóxicos em água do rio e água tratada no Município de Encantado/RS. Ambiente Gestão Desenvolv. 2021, 14, 6–17. [Google Scholar] [CrossRef]
  6. ANVISA. Brazilian Health Regulatory Agency (ANVISA), Resolution N 294, 2019. Available online: https://www.gov.br/anvisa/pt-br/setorregulado/regularizacao/agrotoxicos/monografias/monografias-autorizadas/f/4351json-file-1 (accessed on 20 December 2023).
  7. ANVISA. Brazilian Health Regulatory Agency (ANVISA), Resolution N 294, 2019. Available online: https://www.gov.br/anvisa/pt-br/setorregulado/regularizacao/agrotoxicos/monografias/monografias-autorizadas/q-r-s/4513json-file-1 (accessed on 20 December 2023).
  8. ANVISA. Brazilian Health Regulatory Agency (ANVISA), IN 141, 2022. Available online: https://www.gov.br/anvisa/pt-br/setorregulado/regularizacao/agrotoxicos/monografias/monografias-autorizadas/c/4227json-file-1 (accessed on 20 December 2023).
  9. Zhang, W.; Fan, R.; Luo, S.; Ying, L.; Jin, Y.; Li, Y.; Xiong, M.; Yuan, X.; Jia, L.; Chen, Y. Antagonistic effects and mechanisms of carbendazim and chlorpyrifos on the neurobehavior of larval zebrafish. Chemosphere 2022, 293, 133522. [Google Scholar] [CrossRef]
  10. Wu, C.H.; Lu, C.W.; Hsu, T.H.; Wu, W.J.; Wang, S.E. Neurotoxicity of fipronil affects sensory and motor systems in zebrafish. Pestic. Biochem. Physiol. 2021, 177, 104896. [Google Scholar] [CrossRef] [PubMed]
  11. Monesson-Olson, B.; McClain, J.J.; Case, A.E.; Dorman, H.E.; Turkewitz, D.R.; Steiner, A.B.; Downes, G.B. Expression of the eight GABAA receptor alpha subunits in the developing zebrafish central nervous system. PLoS ONE 2018, 13, e0196083. [Google Scholar] [CrossRef] [PubMed]
  12. Jiang, J.; Wang, L.; Zhang, C.; Zhao, X. Health risks of sulfentrazone exposure during zebrafish embryo-larvae development at environmental concentration. Chemosphere 2022, 288, 132632. [Google Scholar] [CrossRef]
  13. Legradi, J.B.; Di Paolo, C.; Kraak, M.H.S.; van der Geest, H.G.; Schymanski, E.L.; Williams, A.J.; Dingemans, M.M.L.; Massei, R.; Brack, W.; Cousin, X.; et al. An ecotoxicological view on neurotoxicity assessment. Environ. Sci. Eur. 2018, 30, 46. [Google Scholar] [CrossRef]
  14. Kluver, N.; Konig, M.; Ortmann, J.; Massei, R.; Paschke, A.; Kuhne, R.; Scholz, S. Fish embryo toxicity test: Identification of compounds with weak toxicity and analysis of behavioral effects to improve prediction of acute toxicity for neurotoxic compounds. Environ. Sci. Technol. 2015, 49, 7002–7011. [Google Scholar] [CrossRef] [PubMed]
  15. Rosa, J.G.S.; Lima, C.; Lopes-Ferreira, M. Zebrafish Larvae Behavior Models as a Tool for Drug Screenings and Pre-Clinical Trials: A Review. Int. J. Mol. Sci. 2022, 23, 6647. [Google Scholar] [CrossRef] [PubMed]
  16. Syafrudin, M.; Kristanti, R.A.; Yuniarto, A.; Hadibarata, T.; Rhee, J.; Al-onazi, W.A.; Algarni, T.S.; Almarri, A.H.; Al-Mohaimeed, A.M. Pesticides in Drinking Water—A Review. Int. J. Environ. Res. Public Health 2021, 18, 468. [Google Scholar] [CrossRef] [PubMed]
  17. Howe, K.; Clark, M.D.; Torroja, C.F.; Torrance, J.; Berthelot, C.; Muffato, M.; Collins, J.E.; Humphray, S.; McLaren, K.; Matthews, L.; et al. The zebrafish reference genome sequence and its relationship to the human genome. Nature 2013, 496, 498–503. [Google Scholar] [CrossRef]
  18. Kalueff, A.V.; Gebhardt, M.; Stewart, A.M.; Cachat, J.M.; Brimmer, M.; Chawla, J.S.; Craddock, C.; Kyzar, E.J.; Roth, A.; Landsman, S.; et al. Towards a comprehensive catalog of zebrafish behavior 1.0 and beyond. Zebrafish 2013, 10, 70–86. [Google Scholar] [CrossRef] [PubMed]
  19. Gerlai, R. Evolutionary conservation, translational relevance and cognitive function: The future of zebrafish in behavioral neuroscience. Neurosci. Biobehav. Rev. 2020, 116, 426–435. [Google Scholar] [CrossRef] [PubMed]
  20. Gerlai, R. Zebrafish (Danio rerio): A newcomer with great promise in behavioral neuroscience. Neurosci. Biobehav. Rev. 2023, 144, 104978. [Google Scholar] [CrossRef] [PubMed]
  21. Pelkowski, S.D.; Kapoor, M.; Richendrfer, H.A.; Wang, X.; Colwill, R.M.; Creton, R. A novel high-throughput imaging system for automated analyses of avoidance behavior in zebrafish larvae. Behav. Brain Res. 2011, 223, 135–144. [Google Scholar] [CrossRef] [PubMed]
  22. Kimmel, C.B.; Ballard, W.W.; Kimmel, S.R.; Ullmann, B.; Schilling, T.F. Stages of embryonic development of the zebrafish. Dev. Dyn. Off. Publ. Am. Assoc. Anat. 1995, 203, 253–310. [Google Scholar] [CrossRef]
  23. Basnet, R.M.; Zizioli, D.; Taweedet, S.; Finazzi, D.; Memo, M. Zebrafish Larvae as a Behavioral Model in Neuropharmacology. Biomedicines 2019, 7, 23. [Google Scholar] [CrossRef]
  24. Nabinger, D.D.; Altenhofen, S.; Peixoto, J.V.; da Silva, J.M.K.; Gerlai, R.; Bonan, C.D. Feeding status alters exploratory and anxiety-like behaviors in zebrafish larvae exposed to quinpirole. Prog. Neuropsychopharmacol. Biol. Psychiatry 2021, 108, 110179. [Google Scholar] [CrossRef] [PubMed]
  25. Grandjean, P.; Landrigan, P.J. Developmental neurotoxicity of industrial chemicals. Lancet 2006, 368, 2167–2178. [Google Scholar] [CrossRef] [PubMed]
  26. Randlett, O.; Wee, C.L.; Naumann, E.A.; Nnaemeka, O.; Schoppik, D.; Fitzgerald, J.E.; Portugues, R.; Lacoste, A.M.; Riegler, C.; Engert, F.; et al. Whole-brain activity mapping onto a zebrafish brain atlas. Nat. Methods 2015, 12, 1039–1046. [Google Scholar] [CrossRef] [PubMed]
  27. Panula, P.; Sallinen, V.; Sundvik, M.; Kolehmainen, J.; Torkko, V.; Tiittula, A.; Moshnyakov, M.; Podlasz, P. Modulatory neurotransmitter systems and behavior: Towards zebrafish models of neurodegenerative diseases. Zebrafish 2006, 3, 235–247. [Google Scholar] [CrossRef] [PubMed]
  28. d’Amora, M.; Giordani, S. The Utility of Zebrafish as a Model for Screening Developmental Neurotoxicity. Front. Neurosci. 2018, 12, 976. [Google Scholar] [CrossRef] [PubMed]
  29. Ogungbemi, A.; Leuthold, D.; Scholz, S.; Küster, E. Hypo- or hyperactivity of zebrafish embryos provoked by neuroactive substances: A review on how experimental parameters impact the predictability of behavior changes. Environ. Sci. Eur. 2019, 31, 88. [Google Scholar] [CrossRef]
  30. Jordi, J.; Guggiana-Nilo, D.; Bolton, A.D.; Prabha, S.; Ballotti, K.; Herrera, K.; Rennekamp, A.J.; Peterson, R.T.; Lutz, T.A.; Engert, F. High-throughput screening for selective appetite modulators: A multibehavioral and translational drug discovery strategy. Sci. Adv. 2018, 4, eaav1966. [Google Scholar] [CrossRef]
  31. Leuthold, D.; Kluver, N.; Altenburger, R.; Busch, W. Can Environmentally Relevant Neuroactive Chemicals Specifically Be Detected with the Locomotor Response Test in Zebrafish Embryos? Environ. Sci. Technol. 2019, 53, 482–493. [Google Scholar] [CrossRef]
  32. Park, H.; Lee, J.Y.; Park, S.; Song, G.; Lim, W. Developmental toxicity of fipronil in early development of zebrafish (Danio rerio) larvae: Disrupted vascular formation with angiogenic failure and inhibited neurogenesis. J. Hazard. Mater. 2020, 385, 121531. [Google Scholar] [CrossRef]
  33. Fan, R.; Zhang, W.; Li, L.; Jia, L.; Zhao, J.; Zhao, Z.; Peng, S.; Yuan, X.; Chen, Y. Individual and synergistic toxic effects of carbendazim and chlorpyrifos on zebrafish embryonic development. Chemosphere 2021, 280, 130769. [Google Scholar] [CrossRef]
  34. Westerfield, M. The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Danio rerio), 4th ed.; University of Oregon Press: Eugene, OR, USA, 2000. [Google Scholar]
  35. Hedge, J.M.; Hunter, D.L.; Sanders, E.; Jarema, K.A.; Olin, J.K.; Britton, K.N.; Lowery, M. Influence of Methylene Blue or Dimethyl Sulfoxide on Larval Zebrafish Development and Behavior. Zebrafish 2023, 20, 132–145. [Google Scholar] [CrossRef]
  36. OECD 2013a. Fish, Early-Life Stage Toxicity Test. OECD: TG 210. Available online: https://doi.org/10.1787/9789264203785-en (accessed on 11 May 2024).
  37. Clift, D.; Richendrfer, H.; Thorn, R.J.; Colwill, R.M.; Creton, R. High-throughput analysis of behavior in zebrafish larvae: Effects of feeding. Zebrafish 2014, 11, 455–461. [Google Scholar] [CrossRef]
  38. Oliveri, A.N.; Bailey, J.M.; Levin, E.D. Developmental exposure to organophosphate flame retardants causes behavioral effects in larval and adult zebrafish. Neurotoxicol Teratol. 2015, 52, 220–227. [Google Scholar] [CrossRef] [PubMed]
  39. Cadena, P.G.; Cadena, M.R.S.; Sarmah, S.; Marrs, J.A. Folic acid reduces the ethanol-induced morphological and behavioral defects in embryonic and larval zebrafish (Danio rerio) as a model for fetal alcohol spectrum disorder (FASD). Reprod. Toxicol. 2020, 96, 249–257. [Google Scholar] [CrossRef]
  40. Cadena, P.G.; Sales Cadena, M.R.; Sarmah, S.; Marrs, J.A. Protective effects of quercetin, polydatin, and folic acid and their mixtures in a zebrafish (Danio rerio) fetal alcohol spectrum disorder model. Neurotoxicol. Teratol. 2020, 82, 106928. [Google Scholar] [CrossRef] [PubMed]
  41. Altenhofen, S.; Nabinger, D.D.; Wiprich, M.T.; Pereira, T.C.B.; Bogo, M.R.; Bonan, C.D. Tebuconazole alters morphological, behavioral and neurochemical parameters in larvae and adult zebrafish (Danio rerio). Chemosphere 2017, 180, 483–490. [Google Scholar] [CrossRef]
  42. Pérez-Escudero, A.; Vicente-Page, J.; Hinz, R.C.; Arganda, S.; De Polavieja, G.G. IdTracker: Tracking individuals in a group by automatic identification of unmarked animals. Nat. Methods 2014, 11, 743–748. [Google Scholar] [CrossRef] [PubMed]
  43. Romero-Ferrero, F.; Bergomi, M.G.; Hinz, R.C.; Heras, F.J.H.; de Polavieja, G.G. idtracker.ai: Tracking all individuals in small or large collectives of unmarked animals. Nat. Methods 2019, 16, 179–182. [Google Scholar] [CrossRef]
  44. Brastrom, L.K.; Scott, C.A.; Dawson, D.V.; Slusarski, D.C. A High-Throughput Assay for Congenital and Age-Related Eye Diseases in Zebrafish. Biomedicines 2019, 7, 28. [Google Scholar] [CrossRef]
  45. Nery, L.R.; Eltz, N.S.; Hackman, C.; Fonseca, R.; Altenhofen, S.; Guerra, H.N.; Freitas, V.M.; Bonan, C.D.; Vianna, M.R. Brain intraventricular injection of amyloid-beta in zebrafish embryo impairs cognition and increases tau phosphorylation, effects reversed by lithium. PLoS ONE 2014, 9, e105862. [Google Scholar] [CrossRef]
  46. Bridi, D.; Altenhofen, S.; Gonzalez, J.B.; Reolon, G.K.; Bonan, C.D. Glyphosate and Roundup((R)) alter morphology and behavior in zebrafish. Toxicology 2017, 392, 32–39. [Google Scholar] [CrossRef] [PubMed]
  47. Richendrfer, H.; Creton, R. Cluster analysis profiling of behaviors in zebrafish larvae treated with antidepressants and pesticides. Neurotoxicol Teratol. 2018, 69, 54–62. [Google Scholar] [CrossRef] [PubMed]
  48. Lopes-Ferreira, M.; Falcão, M.A.P.; Disner, G.R.; Lima, C. O Modelo Zebrafish e sua Contribuição ao Meio Ambiente. In Recurso água: Tecnologias e Pesquisas para o uso e a Conservação de Ecossistemas Aquáticos; Editora Cubo: São Carlos, Brazil, 2023. [Google Scholar]
  49. Yan, J.; Chen, S.; Zuo, Z.; He, C.; Yi, M. Graphene oxide quantum dot exposure induces abnormalities in locomotor activities and mechanisms in zebrafish (Danio rerio). J. Appl. Toxicol. 2020, 40, 794–803. [Google Scholar] [CrossRef] [PubMed]
  50. Colwill, R.M.; Creton, R. Locomotor behaviors in zebrafish (Danio rerio) larvae. Behav. Process. 2011, 86, 222–229. [Google Scholar] [CrossRef] [PubMed]
  51. Huang, Q.T.; Sheng, C.W.; Jiang, J.; Tang, T.; Jia, Z.Q.; Han, Z.J.; Zhao, C.Q. Interaction of insecticides with heteromeric GABA-gated chloride channels from zebrafish Danio rerio (Hamilton). J. Hazard. Mater. 2019, 366, 643–650. [Google Scholar] [CrossRef] [PubMed]
  52. Gramlich, O.W.; Godwin, C.R.; Wadkins, D.; Elwood, B.W.; Kuehn, M.H. Early Functional Impairment in Experimental Glaucoma Is Accompanied by Disruption of the GABAergic System and Inceptive Neuroinflammation. Int. J. Mol. Sci. 2021, 22, 7581. [Google Scholar] [CrossRef] [PubMed]
  53. Yoon, C.S.; Jin, J.H.; Park, J.H.; Yeo, C.Y.; Kim, S.J.; Hwang, Y.G.; Hong, S.J.; Cheong, S.W. Toxic effects of carbendazim and n-butyl isocyanate, metabolites of the fungicide benomyl, on early development in the African clawed frog, Xenopus laevis. Environ. Toxicol. 2008, 23, 131–144. [Google Scholar] [CrossRef]
  54. Lim, J.; Miller, M.G. The role of the benomyl metabolite carbendazim in benomyl-induced testicular toxicity. Toxicol. Appl. Pharmacol. 1997, 142, 401–410. [Google Scholar] [CrossRef]
  55. Wang, Y.; Lv, L.; Xu, C.; Wang, D.; Yang, G.; Wang, X.; Weng, H.; Wang, Q. Mixture toxicity of thiophanate-methyl and fenvalerate to embryonic zebrafish (Danio rerio) and its underlying mechanism. Sci. Total Environ. 2021, 756, 143754. [Google Scholar] [CrossRef]
  56. Shukla, S.; Jhamtani, R.C.; Dahiya, M.S.; Agarwal, R. Oxidative injury caused by individual and combined exposure of neonicotinoid, organophosphate and herbicide in zebrafish. Toxicol. Rep. 2017, 4, 240–244. [Google Scholar] [CrossRef]
  57. Shen, Q.; Truong, L.; Simonich, M.T.; Huang, C.; Tanguay, R.L.; Dong, Q. Rapid well-plate assays for motor and social behaviors in larval zebrafish. Behav. Brain Res. 2020, 391, 112625. [Google Scholar] [CrossRef] [PubMed]
  58. Wu, Y.; Huang, J.; Deng, M.; Jin, Y.; Yang, H.; Liu, Y.; Cao, Q.; Mennigen, J.A.; Tu, W. Acute exposure to environmentally relevant concentrations of Chinese PFOS alternative F-53B induces oxidative stress in early developing zebrafish. Chemosphere 2019, 235, 945–951. [Google Scholar] [CrossRef] [PubMed]
  59. Santos, M.J.G.; Soares, A.M.V.M.; Loureiro, S. Joint toxicity of three plant protection products to Triticum aestivum (L.) and Brassica rapa (L.). J. Soils Sediments 2011, 11, 990–999. [Google Scholar] [CrossRef]
  60. Loureiro, S.; Svendsen, C.; Ferreira, A.L.; Pinheiro, C.; Ribeiro, F.; Soares, A.M. Toxicity of three binary mixtures to Daphnia magna: Comparing chemical modes of action and deviations from conceptual models. Environ. Toxicol. Chem. 2010, 29, 1716–1726. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Experimental design. (a) After breeding and collection, the embryos were randomly allocated into experimental groups, with 15 animals per test solution. (b) After 144 h post-fertilization (hpf), larvae were randomly removed from polystyrene containers and placed in 48-well plates—one larva per well, with five larvae from each group per plate. After a period of acclimatization, the plates (one at a time) were positioned on an LED plate for video recording and photography. These videos and photographs were later analyzed for thigmotaxis, sensitivity to touch, and exploratory behavior. The LED board was replaced with an LCD screen, and after the acclimatization period, each of the boards was positioned on the screen to reproduce the animation of the optomotor activity test and record the videos for later analysis. (c) Larvae from each group were removed from the wells of 48-well plates and transferred to 6-well plates, with five larvae per well. After the acclimation period, the 6-well plates were placed on top of the LCD screen to play the animation of the bouncing ball test and record videos for later analysis. During all tests, the animals remained in the same test solution. The colors of the wells correspond to the test solutions to which the larvae are exposed.
Figure 1. Experimental design. (a) After breeding and collection, the embryos were randomly allocated into experimental groups, with 15 animals per test solution. (b) After 144 h post-fertilization (hpf), larvae were randomly removed from polystyrene containers and placed in 48-well plates—one larva per well, with five larvae from each group per plate. After a period of acclimatization, the plates (one at a time) were positioned on an LED plate for video recording and photography. These videos and photographs were later analyzed for thigmotaxis, sensitivity to touch, and exploratory behavior. The LED board was replaced with an LCD screen, and after the acclimatization period, each of the boards was positioned on the screen to reproduce the animation of the optomotor activity test and record the videos for later analysis. (c) Larvae from each group were removed from the wells of 48-well plates and transferred to 6-well plates, with five larvae per well. After the acclimation period, the 6-well plates were placed on top of the LCD screen to play the animation of the bouncing ball test and record videos for later analysis. During all tests, the animals remained in the same test solution. The colors of the wells correspond to the test solutions to which the larvae are exposed.
Biomedicines 12 01176 g001
Figure 2. Mean total percentage and standard deviation (SD) of responses of larvae exposed to pesticides and their mixtures in the thigmotaxis test (a) and responses of each group compared to the control group using a one-way ANOVA (F (7,55) = 2.06 p < 0.07) followed by Tukey’s test (C p = 0.82; F p = 0.24; S p = 0.11; C + F p = 0.08; C + S p = 0.06; F + S p = 0.43; C + F + S p = 0.43). Touch sensitivity response test (b) and responses of each group compared to the control group using a one-way ANOVA (F (7,55) = 2.11863 p < 0.05935) followed by Tukey’s test (C p = 0.92; F p = 0.09; S p = 0.92; C + F p = 0.97; C + S p = 1; F + S p = 1; C + F + S p = 1). Legend: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Figure 2. Mean total percentage and standard deviation (SD) of responses of larvae exposed to pesticides and their mixtures in the thigmotaxis test (a) and responses of each group compared to the control group using a one-way ANOVA (F (7,55) = 2.06 p < 0.07) followed by Tukey’s test (C p = 0.82; F p = 0.24; S p = 0.11; C + F p = 0.08; C + S p = 0.06; F + S p = 0.43; C + F + S p = 0.43). Touch sensitivity response test (b) and responses of each group compared to the control group using a one-way ANOVA (F (7,55) = 2.11863 p < 0.05935) followed by Tukey’s test (C p = 0.92; F p = 0.09; S p = 0.92; C + F p = 0.97; C + S p = 1; F + S p = 1; C + F + S p = 1). Legend: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Biomedicines 12 01176 g002
Figure 3. Mean total percentage and SD of larvae that remained immobile during the Exploratory Activity test (a). Groups were compared to the control group using a one-way ANOVA (F (7,55) = 2. 71277, p < 0.01) followed by Tukey’s test (C p = 0.84; F p = 0.01; S p = 0.88; C + F p = 0.70; C + S p = 0.98; F + S p = 0.13; C + F + S p = 0.88), p < 0.05 (*). Mean Speed of larvae swimming during the test (b). Groups were compared to the control group using a one-way ANOVA (F (7,435) = 4.0053, p < 0.01) followed by Tukey’s test (C p = 0.44; F p = 0.37; S p = 0.002; C + F p < 0.05; C + S p = 0.52; F + S p = 0.03; C + F + S p = 0.67), p < 0.05 (*). Mean Total Distance traveled by larvae that swam during the test (c). Groups were compared to the control group using a one-way ANOVA (F (7,435) = 4. 00565, p < 0.01) followed by Tukey’s test (C p = 0.44; F p = 0.37; S p = 0.003; C + F p < 0.05; C + S p = 0.52; F + S p = 0.03; C + F + S p = 0.68), p < 0.05 (*). (d) Ternary contour plot of the variables Carbendazim, Fipronil, and Sulfentrazone for a special cubic model (F (3,380) = 3.004, p = 0.03), evaluating the individual effects and interactions of these pesticides in their mixtures. Carbendazim and Sulfentrazone were more toxic by reducing the total distance traveled (green areas). Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Figure 3. Mean total percentage and SD of larvae that remained immobile during the Exploratory Activity test (a). Groups were compared to the control group using a one-way ANOVA (F (7,55) = 2. 71277, p < 0.01) followed by Tukey’s test (C p = 0.84; F p = 0.01; S p = 0.88; C + F p = 0.70; C + S p = 0.98; F + S p = 0.13; C + F + S p = 0.88), p < 0.05 (*). Mean Speed of larvae swimming during the test (b). Groups were compared to the control group using a one-way ANOVA (F (7,435) = 4.0053, p < 0.01) followed by Tukey’s test (C p = 0.44; F p = 0.37; S p = 0.002; C + F p < 0.05; C + S p = 0.52; F + S p = 0.03; C + F + S p = 0.67), p < 0.05 (*). Mean Total Distance traveled by larvae that swam during the test (c). Groups were compared to the control group using a one-way ANOVA (F (7,435) = 4. 00565, p < 0.01) followed by Tukey’s test (C p = 0.44; F p = 0.37; S p = 0.003; C + F p < 0.05; C + S p = 0.52; F + S p = 0.03; C + F + S p = 0.68), p < 0.05 (*). (d) Ternary contour plot of the variables Carbendazim, Fipronil, and Sulfentrazone for a special cubic model (F (3,380) = 3.004, p = 0.03), evaluating the individual effects and interactions of these pesticides in their mixtures. Carbendazim and Sulfentrazone were more toxic by reducing the total distance traveled (green areas). Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Biomedicines 12 01176 g003
Figure 4. Total trajectory traveled by a zebrafish larva from each experimental and control group in one well of a 48-well plate during the 3 min recording of the Exploratory Activity test. Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Figure 4. Total trajectory traveled by a zebrafish larva from each experimental and control group in one well of a 48-well plate during the 3 min recording of the Exploratory Activity test. Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Biomedicines 12 01176 g004
Figure 5. Mean total percentage of zebrafish larvae alignment to the well during the Optomotor Response test. If the larva is inside the red demarcated area, it is considered high alignment, and outside the area is low alignment. It was analyzed by Bowker’s symmetry test, p < 0.05 (*), to the same group among initial, middle, and final test positions. (Contr p < 0.05; C p < 0.05; F p = 0.16; S p = 1; C + F p = 0.16; C + S p = 0.88; F + S p < 0.05; C + F + S p < 0.05). Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone; I—Initial; M—Middle; F—Final.
Figure 5. Mean total percentage of zebrafish larvae alignment to the well during the Optomotor Response test. If the larva is inside the red demarcated area, it is considered high alignment, and outside the area is low alignment. It was analyzed by Bowker’s symmetry test, p < 0.05 (*), to the same group among initial, middle, and final test positions. (Contr p < 0.05; C p < 0.05; F p = 0.16; S p = 1; C + F p = 0.16; C + S p = 0.88; F + S p < 0.05; C + F + S p < 0.05). Key: Contr—Control; C—Carbendazim; F—Fipronil; S—Sulfentrazone; I—Initial; M—Middle; F—Final.
Biomedicines 12 01176 g005
Figure 6. Mean total percentage and SD of larvae that exhibited escape behavior about the visual stimulus and stayed longer in the upper part of the well, area without stimulus, during the Bouncing Ball test (a). The response of each group was compared to the control group by one-way ANOVA (F (7,55) = 17.11, p < 0.01) followed by Tukey’s test (C p < 0.05; F p < 0.05; S p < 0.05; C + F p = 0; C + S p < 0.05; F + S p = 0.03; C + F + S p < 0.05), p < 0.05 (*). Percent grouping of larvae during the Bouncing Ball test (b), groups were compared to the control group using a one-way ANOVA (F (7,55) = 6.20, p < 0.01) followed by Tukey’s test (C p < 0.05; F p < 0.05; S p < 0.05; C + F p < 0.05; C + S p < 0.05; F + S p < 0.05; C + F + S p = 0.01), p < 0.05 (*). Key: C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Figure 6. Mean total percentage and SD of larvae that exhibited escape behavior about the visual stimulus and stayed longer in the upper part of the well, area without stimulus, during the Bouncing Ball test (a). The response of each group was compared to the control group by one-way ANOVA (F (7,55) = 17.11, p < 0.01) followed by Tukey’s test (C p < 0.05; F p < 0.05; S p < 0.05; C + F p = 0; C + S p < 0.05; F + S p = 0.03; C + F + S p < 0.05), p < 0.05 (*). Percent grouping of larvae during the Bouncing Ball test (b), groups were compared to the control group using a one-way ANOVA (F (7,55) = 6.20, p < 0.01) followed by Tukey’s test (C p < 0.05; F p < 0.05; S p < 0.05; C + F p < 0.05; C + S p < 0.05; F + S p < 0.05; C + F + S p = 0.01), p < 0.05 (*). Key: C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Biomedicines 12 01176 g006
Figure 7. Cluster analysis of behaviors of larvae exposed to single and mixed pesticides at 6 dpf. p value compared to the control group as a variable. The green bars indicate values of 1–2, and the black bars indicate a value of 3. Key: C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Figure 7. Cluster analysis of behaviors of larvae exposed to single and mixed pesticides at 6 dpf. p value compared to the control group as a variable. The green bars indicate values of 1–2, and the black bars indicate a value of 3. Key: C—Carbendazim; F—Fipronil; S—Sulfentrazone.
Biomedicines 12 01176 g007
Table 1. Concentrations of the single compounds, binary, and ternary mixtures used to measure behavioral changes in zebrafish larvae exposed to the pesticides carbendazim (C), fipronil (F), and sulfentrazone (S). The layout of the simplex centroid design method for evaluating the interactions between the chemical compounds.
Table 1. Concentrations of the single compounds, binary, and ternary mixtures used to measure behavioral changes in zebrafish larvae exposed to the pesticides carbendazim (C), fipronil (F), and sulfentrazone (S). The layout of the simplex centroid design method for evaluating the interactions between the chemical compounds.
Simplex Centroid Design for Behavioral Assessment in 6 dpf (Days Post-Fertilization)
ComponentsCarbendazim (mg/L)Fipronil (mg/L)Sulfentrazone (mg/L)
C0.2000.00.0
F0.00.0500.0
S0.00.00.200
C + F0.1000.0250.0
C + S0.1000.00.100
F + S0.00.0250.100
C + F + S0.0600.0170.025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gomes, S.d.S.; da Silva, J.F.; Padilha, R.M.O.; de Vasconcelos, J.V.A.; Negreiros Neto, L.G.d.; Marrs, J.A.; Cadena, P.G. Behavioral Effects of the Mixture and the Single Compounds Carbendazim, Fipronil, and Sulfentrazone on Zebrafish (Danio rerio) Larvae. Biomedicines 2024, 12, 1176. https://doi.org/10.3390/biomedicines12061176

AMA Style

Gomes SdS, da Silva JF, Padilha RMO, de Vasconcelos JVA, Negreiros Neto LGd, Marrs JA, Cadena PG. Behavioral Effects of the Mixture and the Single Compounds Carbendazim, Fipronil, and Sulfentrazone on Zebrafish (Danio rerio) Larvae. Biomedicines. 2024; 12(6):1176. https://doi.org/10.3390/biomedicines12061176

Chicago/Turabian Style

Gomes, Samara da Silva, Jadson Freitas da Silva, Renata Meireles Oliveira Padilha, João Victor Alves de Vasconcelos, Luís Gomes de Negreiros Neto, James A. Marrs, and Pabyton Gonçalves Cadena. 2024. "Behavioral Effects of the Mixture and the Single Compounds Carbendazim, Fipronil, and Sulfentrazone on Zebrafish (Danio rerio) Larvae" Biomedicines 12, no. 6: 1176. https://doi.org/10.3390/biomedicines12061176

APA Style

Gomes, S. d. S., da Silva, J. F., Padilha, R. M. O., de Vasconcelos, J. V. A., Negreiros Neto, L. G. d., Marrs, J. A., & Cadena, P. G. (2024). Behavioral Effects of the Mixture and the Single Compounds Carbendazim, Fipronil, and Sulfentrazone on Zebrafish (Danio rerio) Larvae. Biomedicines, 12(6), 1176. https://doi.org/10.3390/biomedicines12061176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop