3.1 Impact of hydroquinone on Daphnia magna
Figure 1a shows the
D. magna dose-response curve to HQ.
The calculated EC
50 for HQ is 0.142 (0.104-0.204) µg/mL, indicating high toxicity of this product on
D. magna. The toxic effects of HQ on this organism have been documented in previous studies [
43] and when calculating the EC
50, the values obtained after 24 hours of exposure, are very similar to ours, with EC
50=0.150 µg/mL [
44,
45]. After 48 hours of exposure, values are slightly higher at 0.25-0.28 µg/mL [
17]. Interestingly, another aquatic crustacean species, also belonging to Branchiopoda,
Ceriodaphnia dubia, shows a very similar sensitivity to HQ as
D. magna with EC
50 values of 0.15 µg/mL as well [
22].
D. magna is a good indicator of water quality since it is exposed to toxics through a dual pathway: surface exposure and also through its diet as it is a filter-feeding organism. HQ is a relatively small molecule (MW = 110.11 g/mol) and electrically neutral, with a pKa of approximately 9.9 and 11.6 [
35], which might facilitate its passage through cell membranes. Changes in membrane permeability can affect the integrity of the cell membranes of
D. magna, subsequently altering cellular homeostasis and leading to cell death. However, it is not very lipophilic compound (LogKow=0.59) [
35].
On the other hand, it is soluble in water (73 g/L at 25°C) [
33] which enhances its bioavailability. Therefore, the digestive tract may be the main route of exposure to these organisms, facilitating the entry of HQ into
D. magna, which could lead to cardiac [
46] and nervous [
47] disturbances. It could also act by inducing oxidative stress [
48] or affecting the protein content in the hemolymph, as observed in other invertebrates [
49]. Similar to benzene, HQ can inhibit the activity of certain enzymes such as topoisomerase II [
50], negatively impacting essential cellular processes for the survival of
D. magna. This, in conjunction, would explain the high toxicity of HQ observed on this organism.
3.2. Impact of hydroquinone on V. fisheri
Toxicity of HQ to the bacteria
V. fischeri is illustrated in
Figure 1b and the obtained EC
50 was 1.446 (1.155-1.796) µg/mL. Limited data exist on the toxicity of HQ to
V. fischeri, as studies typically focus on the toxicity of by-products, including HQ, generated during the decomposition of various products such as paracetamol [
51], benzidine [
52], benzoquinone [
53], sulfamethoxazole[
54], sulfanilamide [
55] or clofibric acid [
56] among others. It is noteworthy that almost all studies agree that HQ is one of the most toxic by-products, even more than the original product.
The EC
50 value for
V. fischeri exposed to HQ (as dimethomorph intermediate on TiO
2 suspension) in a 2% NaCl solution was measured at 0.08 mg/L[
57] but the exposure time was only 5 minutes. A. Santos et.al., [
10] reported an EC
50 of 0.041 mg/L (15 minutes) in
V. fischeri during the Catalytic Oxidation of phenol. These results are challenging to compare due to different experimental conditions, and in our case, the exposure was for 30 minutes. Nevertheless, all results suggest that HQ is highly toxic to this aquatic indicator.
The Gram-negative outer covering of V. fischeri may partially shield the bacterium from intracellular exposure to HQ, acting as a selective barrier. Due to its size, HQ may face challenges in traversing the porins of the outer membrane of the Gram-negative wall or interacting with its lipopolysaccharides. Alternatively, it could be expelled by efflux pumps. This may explain its somewhat lower toxicity compared to D. magna. However, once inside the prokaryotic cell, it is likely to have toxicity mechanisms similar to those observed in D. magna.
To the best of our knowledge, no information is available regarding the mechanism of action of HQ on
V. fisheri. However, documented inhibitory effects on the growth of pathogenic bacteria such as
Pseudomonas Aeruginosa,
Klebsiella pneumoniae, and
Escherichia coli have been reported [
58,
59,
60]. Additionally, studies on other bacteria within the genus
Vibrio [
61,
62] observe antimicrobial activity of HQ derivatives. Interestingly, these derivatives appear to downregulate genes of
Vibrio spp implicated in motility, protease synthesis, indol, and capsular polysaccharide production, suggesting a potential mechanism of action [
61].
The substantial impact of HQ on both V. fischeri and D. magna suggests potential significant effects on river ecosystems. However, assessing its effects on complete communities, such as microbial ones, is essential for a more realistic diagnosis.
In
Figure 2, the genetic sequencing of river microbial communities can be seen. River microorganism sequencing generated a total of 65615 reads, all of which passed quality filters with a 100% success rate. Sequencing comprehensively covered all taxonomic levels, achieving >95% for phylum, class, and order, >50% for family and genus, and 23.33% for species.
Figure 2a displays the most prevalent taxa (>2%) for river microorganisms at each taxonomic level. In
Figure 2b, A visual representation is provided, illustrating the most prominently observed phyla with pie chart slices indicating their respective percentages.
3.3. Impact on river microbial communities: growth and Community-level physiological profiling (CLPP)
Three predominant phyla were: Cyanobacteria (41.4% of the bacterial reads), Proteobacteria (29%), and Bacteroidetes (12.2%). Notably, 16.5% of bacterial reads remained unidentified, highlighting the presence of novel sequences in this study.
The Cyanobacteria phylogenetically belong to oxygenic phototrophic bacteria frequently found in rivers [
63,
64]. Almost all Cyanobacteria were classified within the class Oscillatoriophycidaeae (94.4%), with the majority falling under the order Chroococcales, a dominant group in freshwater biotopes[
65].
Within Proteobacteria, we encountered three predominant classes: Gammaproteobacteria and Betaproteobacteria, exhibiting similar abundances at 34.7% and 31.5%, respectively, and Bacteroidetes at 13.22%. Proteobacteria, a prolific phylum of Gram-negative bacteria in freshwater bacterial communities [
66] demonstrates rapid growth in response to organic nutrients (Madigan et al., 2015). Gammaproteobacteria, known for its high taxonomic diversity, featured Alteromonadales as the most prevalent order (31%), a representative of river microbial communities [
67,
68]. Notably, the order Pseudomonadales (8% of Gammaproteobacteria reads) includes the Pseudomonadaceae and Moraxellaceae families, some of whose members, such as Pseudomonas, play an active role in the degradation of phenolic compounds [
69,
70]. Betaproteobacteria were predominantly of the Burkholderiales order (74.7%), and among the Alphaproteobacteria, Rhodobacterales stood out (42.5%).
Within Bacteroidetes, Gram-negative anaerobic bacteria with significant involvement in the degradation of humic materials and polymers [
71], we found two dominant classes: Flavobacteria (54% of Bacteroidetes) and Sphingobacteriia (40.5%).
Freshwater microbial communities have been suggested as excellent bioindicators for assessing the impact of micropollutants in river ecosystems [
72] because disruptions at this level can have consequences throughout the trophic levels [
72,
73], leading to unpredictable effects on the ecological balance of the aquatic environment [
74]. These communities serve as the foundation of the aquatic food web, particularly among primary producers, and also play a significant role in organic matter decomposition, thereby contributing to nutrient cycling and energy exchange, as well as the degradation of pollutants [
75,
76].
While our results indicate high toxicity of HQ in various aquatic indicators, it is surprising how the impact on the growth and metabolic capacity of these microbial communities appears to be buffered, as if these communities could effectively withstand HQ's toxic effects.
In
Figure 3 the effect of HQ on river microbial communities, measured as AWCD, can be seen. Furthermore,
Figure 4 illustrates the impact of this product on the microbial profile of the community, compared to the control.
Furthermore, among the diversity of taxa, there may also exist varying metabolic capabilities, with certain bacteria potentially possessing mechanisms capable of degrading HQ. These microorganisms may derive greater advantages than others, potentially reshaping ecological interactions where the dominant flora that degrades HQ hydroquinone can be gradually formed [
77,
78]. As can be seen (
Figure 3), although at the beginning there were small differences, after 72 hours the growth of the community exposed to HQ practically matched those of the control, possibly due to these readjustments within the community.
Figure 3.
AWCD values of metabolized substrates in Biolog EcoPlates based on 168h incubation of river microorganisms exposed to Hydroquinone. Each point is the average value of three replicates.
Figure 3.
AWCD values of metabolized substrates in Biolog EcoPlates based on 168h incubation of river microorganisms exposed to Hydroquinone. Each point is the average value of three replicates.
Zhang et. al. [
77] observed that concentrations of HQ at 100 mg/L (the highest tested in this study) in wastewater treatment plants resulted in the establishment of a stable community dominated by the same taxa we have identified in our samples (Cyanobacteria, Proteobacteria, and Bacteroidetes). These taxa showed minimal variation in their relative abundance compared to the control [
27]. Specifically, the abundance of Cyanobacteria remained largely unaffected
, Bacteroidetes showed a slight increase, and
Proteobacteria exhibited a minor decrease in this study. The limited impact on Cyanobacteria, which constitute nearly half of our samples, may explain the minimal metabolic changes observed in our study, even at the highest concentration. Proteobacteria, as the largest group of gram-negative bacteria with a wide range of metabolic pathways and a major role in the degradation of phenolic compounds [
79], could withstand the HQ impact despite experiencing a modest decline (on the order of 10% at 100 mg/L, according to Zhang et. al. In fact, several members of this group present in our samples have been reported to be able to metabolize HQ.
Among the Gammaproteobacteria we found Pseudomonadales (specifically
Pseudomonas genera) and members of the Moraxellaceae family, both proficient in utilizing and degrading HQ [
69,
70,
80,
81]. Additionally, within the Betaproteobacteria, we observed the presence of Burkholderiales, also capable of following HQ degradation pathways[
5]
.
On the other hand,
Bacteroidetes are known for their capacity to degrade various complex carbon compounds, including HQ [
82], potentially increasing in number to compensate for the loss of Proteobacteria.
Beyond these changes in community structure reported, our results demonstrate that the final result of this taxonomic rearrangement within the community is that the metabolic capacity of the entire community is minimally affected by HQ (
Figure 4). Only a decrease in the ability to metabolize polymers at the highest concentration of 100 µg/mL (p = 0.02) appears to occur. All other changes in the metabolic profile of the microbial community are not significant at any of the concentrations tested. This would be consistent with studies showing that functional genes for carbohydrate metabolism and energy metabolism were maintained at a high level following HQ exposure [
27].
Figure 4.
Metabolic effect differentiation by carbon sources of the river microorganisms exposed in different concentrations to Hydroquinone respect to the control (Y axis). Each point is the average value of three replicates. The significance of differences from the control is indicated by p-values (t—Student), and the dispersion of values among the three replicates is represented by the coefficient of variation (CV).
Figure 4.
Metabolic effect differentiation by carbon sources of the river microorganisms exposed in different concentrations to Hydroquinone respect to the control (Y axis). Each point is the average value of three replicates. The significance of differences from the control is indicated by p-values (t—Student), and the dispersion of values among the three replicates is represented by the coefficient of variation (CV).
Therefore, although initially, the microbial flora was stressed by the influent HQ, which may even trigger the secretion of secondary metabolites that increase toxicity [
27,
28] the microbial community, after a succession of biological communities, gradually forms a dominant flora capable of degrading or tolerating HQ. As a result, the metabolic capacity of the microbial community remains stable, and it is foreseeable that the impact of HQ on rivers will be minimal.
In many countries, the implementation of maximum concentration limits for the industrial discharge of phenols has been established [
3,
83]. These limits typically range from low mg/L to μg/L, depending on the specific discharge location and the flow characteristics of the watercourse (EC, Commission Implementing Decision (EU), 2018) [
120]. While these levels may provide protection for microbial communities, it is not necessarily guaranteed for other aquatic organisms, such as
D. magna.
3.5. Impact of Hydroquinone on Eisenia fetida
Our results demonstrate that, despite
E. fetida being the most resilient bioindicator among the four tested, it still exhibits detectable toxicity. The dose-response of the earthworm exposed to HQ can be seen in
Figure 5b with LC
50 of 234.05 (184.13-281.18) mg/kg. When comparing the toxicity values of HQ on
E. fetida to other phenolic compounds of plant origin (non-quinones), such as tannic acid, the latter shows much higher values (LC
50 > 2000 µg/L) [
88]. However, to the best of our knowledge, the impact of HQ on earthworms, particularly
E. fetida, has not been previously investigated. While some evidence of toxicity can be found in the literature, it often pertains to compounds within the HQ family or chemically distinct derivatives, and it may involve different earthworm species. For instance, exposure studies involving various polyesters containing HQ, among other compounds, showed an
E. fetida survival rate exceeding 80% after 14 days, suggesting a moderate levels of toxicity to these bioindicators [
89].
Interestingly, Osman [
90] observed that additional earthworm species, including
L. rubellus and
A. chlorotica, seem to exhibit susceptibility to oxidative stress induced by quinones. This susceptibility may be attributed to their deficiency or notably low levels of DT-diaphorase, an enzyme recognized for its significant role in quinone detoxification.
Exposure of earthworms to HQ can occur through the ingestion of particles carrying the active product [
91] and through percutaneous means. Earthworms possess a highly water-absorbent and water-loss-tolerant cuticle, allowing for significant water exchange through the body wall [
92]. HQ's relatively low molecular weight and slight hydrophobic nature could enable its permeability in biological membranes [
93]. However, it is likely that ingestion, in this case, is what triggers the cytotoxic effects.
Earthworms play a crucial role in soil health and fertility as they decompose organic matter and mix the soil, improving its structure and enhancing its ability to retain water and nutrients, thereby allowing plants to access these nutrients. Therefore, their decline or reduction can have significant consequences for soil fertility [
94].
The activity of these organisms is intimately connected to that of soil microorganisms, as earthworms have an important role in promoting microbial activity, likely by feeding on microorganisms or by selecting and stimulating specific microbial groups [
95].
3.6 Impact on Soil microbial communities: growth and Community-level physiological profiling (CLPP)
Figure 6 show the great diversity of soil taxa. In this case the total reads were of 61347 and the 100% passing quality filtering. It was possible to identify >90% of taxa at the taxonomic level of Phylum, Class, Order and Family, 88.63% of Genus and only 24.23 % of species.
Figure 6a displays the relative abundance of the main taxons within each taxonomic level of the most prevalent taxa (>2%).
Figure 6b, a visual representation highlights the most prominently detected phyla.
In our samples, we observed a predominance of two bacterial phyla: Actinobacteria, which constituted 48.7% of the bacterial reads, and Proteobacteria, making up 34.6% of the composition. Additionally, we detected a smaller proportion of Firmicutes, accounting for 8.0% of the total reads. This taxonomic distribution aligns with the typical bacterial diversity encountered in uncontaminated edaphic ecosystems where Proteobacteria are usually very abundant [
96,
97], Actinobacteria phyla are well represented [
71] and Firmicutes are frequently detected [
98,
99,
100].
Among the Actinobacteria, the Class Actinobacteria predominates (68.0 %), practically all belonging to the order Actynomycetales, ubiquitous in different soil types [
98,
99,
101,
102].
More than half of Proteobacteria were Alphaproteobacteria (60.9%) followed by Deltaproteobacteria (17.9%) and Gammaproteobacteria (15.7%). Almost all Alphaproteobacteria are of the order Sphingomonadales with a small representation of the order Rhizobiales (8.11% of the Alphaproteobacteria).
Figure 6.
(a) Relative abundance of genetically sequenced microorganisms from river within their taxonomic classifications at each level. (b) Illustration of phyla that are most prominently observed in soil. The significance of differences from the control is indicated by p-values (t—Student), and the dispersion of values among the three replicates is represented by the coefficient of variation (CV).
Figure 6.
(a) Relative abundance of genetically sequenced microorganisms from river within their taxonomic classifications at each level. (b) Illustration of phyla that are most prominently observed in soil. The significance of differences from the control is indicated by p-values (t—Student), and the dispersion of values among the three replicates is represented by the coefficient of variation (CV).
Among the Deltaproteobacteria, Myxococcales predominate (52.9%), all of them belonging to the family Cystobacteraceae and the genus Cystobacter. In Gammaproteobacteria, all the Pseudomonadaceae family are Pseudomonas.
Among the Firmicutes, Bacilli (54.1%) and Clostridia (33.3%) are the predominant class.
In the
Figure 7, the effects of HQ on community growth measured as AWCD are depicted.
As can be observed, microbial communities also appear to withstand HQ exposure well, except at concentrations greater than 100 µg/mL (p=0.05). In this case, there is no initial growth decline followed by subsequent recovery, as seen in the case of river microorganisms. Instead, at 100 µg/mL, growth is partially inhibited right from the beginning of HQ exposure. This heightened sensitivity of soil microbial communities compared to aquatic ones is consistent with findings from other studies where soil or sediment microorganisms seem to be more vulnerable to potentially toxic compounds than aquatic microorganisms [
103,
104]. This observation has also been noted for products or extracts of plant origin [
41,
105].
Moreover, at the metabolic level (see
Figure 8), the concentration of 100 µg/mL induces a significant decrease in the ability to metabolize not only polymers (p=0.012), as observed in the case of river microbial communities, but also carbohydrates (p=0.006). Nevertheless, at lower concentrations, there are no significant changes in the metabolic profile (p>0.05) for any metabolite.
There are very few studies that have examined the effect of HQ on soil microbial communities. Nevertheless, there is evidence that HQ may indeed impact microbial growth. Chen [
106] observed that soils amended with HQ experience a decrease in the growth of cultivable microbial populations, with HQ being the most toxic dihydroxybenzene compared to other phenolic compounds such as resorcinol and catechol. It has also been reported that soil microorganisms' exposure can lead to minor changes, such as an increase in the relative abundance of groups involved in fermentation and cellulolysis [
31], which, in some way, may account for the slight variations in the metabolic profile we have detected.
The use of HQ as a urease inhibitor [
18] to prevent urease from breaking down into urea, thus increasing the availability of NH
3/NH
4+ for plant uptake [
107], has led to a limited number of studies examining the effect of HQ on soil microorganisms, especially in nitrification and denitrification processes, with varying results. On one hand, HQ, in line with our findings, appears to induce minimal changes in the community composition and functional profiles of the soil microbial community, with little impact on ureolysis groups [
31,
108]. However, other authors have reported that ammonia oxidation microbes were inhibited following HQ application [
32] or that HQ delays urea hydrolysis, subsequently affecting nitrification and denitrification [
109]. Nevertheless, there are limited reports on the effects of long-term HQ application on the soil nitrification and denitrification microbial community. Our results, however, do not indicate significant changes in the capacity to metabolize substrates potentially involved in nitrogen metabolism, such as carboxylic and ketonic acids, amino acids, or amines and amides.
The resilience exhibited by these soil microorganisms to HQ at concentrations below 100 µg/mL may stem from strategies akin to those described for aquatic microorganisms. In this scenario, we also encounter a significant diversity of taxonomic groups, making the replacement of sensitive species with more resistant ones, capable of degrading HQ, an expected occurrence.
According to genetic sequencing, we have identified several genera, including
Pseudomonas (3.31% of the total reads) and
Burkholderia (
Figure 6, within the section "other Proteobacteria") and members belonging to the order Rhizobiales, all of them able to metabolize HQ [
69,
110,
111]. Furthermore, as previously discussed, taxonomic groups within the Sphingomonadaceae family (constituting 19.6% of total soil reads) have been found to possess mechanisms for safeguarding against HQ exposure [
112].
Other mechanisms, such as the production of specific enzymes for phenolic compound detoxification, as described in Actinomycetales members (constituting 32.94% of total reads in our samples) [
113], and the formation of biofilms, as demonstrated by Corynebacteriaceae within the Actinomycetales order (representing 32.94% of total reads), able to metabolize HQ [
114], are also plausible. In fact, the microbial diversity, structure, and function of a biofilm imparts a high metabolic capacity. It has been reported that biofilms are capable of removing more than 95% of phenolic compounds, including HQ) [
30].
Therefore, our findings suggest that unless occurring at exceptionally high concentrations rarely encountered in the environment, the impact of HQ on soil microbial communities is likely to have minimal effects on microbial growth and will not significantly impair their metabolic capacity.