1. Introduction
Dyeing clothes is a prehistoric process. This process involved the application of early natural dyes for furs and textiles of vegetable origin, though some dyes were of animal origin. More complex coloring materials were developed over thousands of years. The woad (natural indigo) for example was obtained from the plant
Indigofera tinctoria, and the Tyrian purple was extracted from the gland of a purple snail and developed by the Phoenicians, whereas the Alizarine was taken out from madder Campeachi wood [
1]. By the end of the 19th and early 20th centuries the synthetic dyes industry was established in many countries and thousands of dye molecules have been synthesized and produced at a large scale [
2].
Despite the chemical diversity of dyes, dye molecules share a common chemical structure. In fact, each dye molecule has four components namely a chromophore group, an auxochrome group, a solubilizing group, and a matrix. The chromophore groups are responsible for the absorption of the light energy and the creation of the dye's color through the excitation of electrons. The auxochrome groups help with the dye fixation into the support, while, the solubilizing ones ensure the solubility of the molecule in water or organic solvents. The remaining parts of the dye molecule form the matrix or the skeleton [
1]
.
Dyes are classified according to several parameters including: color, chemical structure, application, manufacturer, synthesis route, fastness, and date invented. However, based on their chemical structure and the chromophore groups, the following dye families were identified: azo, anthraquinone, nitroso, nitro, indigoid, cyanine, phthalocyanine, and triphenylmethane [
3].
Azo dyes are characterized by two aromatic groups linked to each other by an azo bond (-N=N-). They are classified based on the number of azo linkages mono azo dyes, diazo dyes, etc. the number of azo groups varies from 1 to 4. Other than textile industries, this family of dyes is used in various fields such as pharmaceuticals, cosmetics, food, paint, paperwork, etc. Their success is mostly due to the stability of coloring, the ease of a coupling reaction between the dyes and the support, the high molar extinction coefficient (capacity to absorb light), the flexibility of the coloring structure, and their adaptability to a variety of applications [
4,
5,
6].
Approximately 70% of the dyes, used in the textile industry, are of the azo type. However, during the coloring process, non-adsorbed dyes are estimated between 15 and 20% and are discharged into the wastewater [
7,
8,
9]. Due to their toxicities, industries using this type of dyes are currently attempting to minimize their negative impact on the environment. This includes improving their binding to the matrix or their degradation once discharged into industrial wastewater, using biological or physicochemical processes.
Many studies have demonstrated that the released sewage contains, other than dyes, toxic molecules like heavy metals. Once released in the environment, the wastewater may affect both the human health and ecosystem [
10,
11,
12,
13]. Many health issues including cancer, chronic diseases, and skin irritation have been associated with exposure to azo dyes [
10]. Besides, the death of aquatic organisms and the stunting of plant growth were mentioned as a consequence of the release of untreated textile wastewater [
14]. To treat sewage from the textile industry, many attempts have been made and many physicochemical methods were developed (e.g., filtration, adsorption, coagulation/flocculation). Those treatments were mostly used at the outset. Nevertheless, their unwanted outcomes like the formation of secondary mud, the limited efficacy, and the high cost, has prompted industries to look for alternative biological methods that are especially eco-friendly and low-cost and where plants or microorganisms and/or their enzymes can be used [
15]. Yeast [
16], bacteria [
15,
17,
18,
19], algae [
20], and fungi [
21] have widely been used for this purpose. Several studies demonstrated the efficiency of white-rot fungi such as
Trametes versicolor [
22],
T. trogii [
23,
24], and
C. gallica [
7,
25], and other fungi like
Aspergillus niger [
26] in the removal of textile dyes using their enzymes or biomasses.
White-rot fungi
secrete a number of oxidoreductases that are involved
in lignin depolymerization [
27,
28]
. These oxidoreductases encompass heme-containing peroxidases (manganese, lignin, and versatile peroxidases) and copper-dependent polyphenol oxidases named laccases (E.C. 1.10.3.2). Laccases from white-rot fungi exhibit a higher redox potential (0.720-0.790 V) compared to other fungal, bacterial, or plant laccases (0.400-0.700 V). High redox potential laccases do not oxidize lignin directly but through small aromatic compounds (laccase-mediators) that can attack lignin after their oxidation in the active site [
29,
30,
31]. Although
laccases are widely used, their efficiency in removing pollutants is sometimes limited; this has prompted the use of a laccase-mediator system to enhance the laccase activity. Mediators allowed the active center of the enzyme to interact with large molecules of substrates or substrates with a high redox potential. Several laccase mediators have been studied including mediators of natural molecules (e.g., 3-hydroxy-anthranilic acid, syringaldehyde, vanillin etc.) or synthetic molecules such as HBT, TEMPO, Violoric acid etc.) [
32].
In this paper, we aim to study the biodegradation of the four azo-bond dye Sirius grey by the laccase mediator system and to optimize its decolorization conditions using a response surface methodology approach.
2. Materials and Methods
2.1. Chemicals
2,6-Dimethoxyphenol (DMP) and 1-hydroxybenzotrizole (HBT) were obtained from Sigma-Aldrich. Sirius grey GB was obtained from a textile factory located in Ksar Helal (Tunisia). Its properties are summarized in
Table 1 and its chemical structure is depicted in
Figure 1.
2.2. Media and culture conditions
Potato-Dextrose-Agar (PDA) medium was used for a short-term conservation of the fungal strain. After growth at 30 °C, plates were stored at 4 °C, and sub-cultured monthly. Laccase production by
C. gallica was performed in liquid medium as described by Zouari-Mechichi et al. [
24]. Media were inoculated by glass-beads homogenized mycelium (1%) and laccase production was induced by CuSO
4 at a final concentration of 300 μM. Cultures were incubated at 30 °C and 160 rpm, and when maximum laccase production was reached (8–10 days), laccase-rich supernatant was separated from biomass by filtration on 3 M filter paper and stored at -20 °C until utilization.
2.3. Fungal strain isolation and identification
2.3.1. Isolation
The fungus used in this study was the newly isolated strain BS9. To isolate this strain, a piece of fungal crust growing on a decayed Eucalyptus (Eucalyptus globulus Labill.) wood was inoculated on a PDA medium and incubated at 28 °C, the growing mycelium was transferred several times on the same solid medium until a pure culture was obtained. The culture was tested for laccase activity production on a PDA medium supplemented with 0.01% of guaiacol. When oxidized the non-colored substrate turns orange, indicating the production of a phenol oxidase activity.
2.3.2. Fungal DNA extraction, amplification and sequencing
Total genomic DNA was extracted from fungal culture using DNA extraction kit EM13 (Blirt S.A., Gdańsk, Poland) according to the manufacturer's instructions. The internal transcribed spacer (ITS) was selected as markers for molecular identification purposes. ITS is widely considered as a universal fungal barcode of first choice [
33]. The following primers pair was used: ITS1F (5'-CTT GGT CAT TTA GAG GAA GTA A-3') and ITS4R (5'-TCC TCC GCT TAT TGA TAT GC-3') for ITS amplification [
34]. PCR mix consisted of: 10 µL of 2x TaqNova-RED PCR Master Mix (Blirt S.A., Gdańsk, Poland), 4 µL of water, 1.5 µL of each forward and reverse primers in 10 µM concentration, and 3 µL of template DNA for a final reaction volume of 20 µL. Presence of the expected product was checked on 1% agarose gel with Midori Green (Nippon Genetics Europe, Düren, Germany) PCR products were cleaned using DNA purification kit EM26 (Blirt S.A., Gdańsk, Poland) and bidirectionally sequenced using BigDye™ Terminator v3.1 Cycle Sequencing Kit (Thermofisher Scientific, Waltham, MA, USA). The product was cleaned using Zetadex-50 Fine (EMP BIOTECH GMBH, Berlin, Germany) and sent to the external company Genomed (Warsaw, Poland) for reading. The obtained sequences were deposited in the GenBank database under the accession number: OR234862.
2.3.3. Phylogenetic analysis
The obtained reads were assembled using the CAP3 algorithm implemented in UGENE v.37.0 [
35]. Resulting consensus sequences were compared with the NCBI nucleotide database using BLASTn algorithm for preliminary taxonomic placement. Additional 9 reference sequences were retrieved from GenBank to prepare a phylogenetic tree (
Table 2). All sequences were aligned using the muscle algorithm as implemented in the SeaView program [
36]. Subsequently, the alignments were trimmed in the trimAl program using the automated algorithm [
37]. The phylogenetic tree was calculated using the maximum likelihood approach in RAxML software [
38] using the GTR+G+I substitution model. The tree robustness was assessed by bootstrap analyses with 1000 replicates. The isolate was assigned to species based on its position on the phylogenetic tree.
2.4. Enzyme assay
Laccase assay was performed with 10 mM 2,6-dimethoxyphenol (DMP) as substrate in 100 mM tartrate buffer, pH 5 (ε
469 nm = 27,500 M
−1 cm
−1) [
39]. One unit of laccase activity was defined as the amount of enzyme oxidizing 1μmol of substrate per minute.
2.5. Dye decolorization experiments
All experiments were performed in 2 mL disposable cuvettes containing 1.5 mL final reaction volume. The reaction mixture contained 100 mM tartrate buffer pH 3 to 6, HBT, 50 mg/L of dye, and culture filtrate (0.5 U/mL laccase). The reaction was initiated by the addition of culture filtrate. The decolorization was followed by measuring the absorbance of the dye solution at the maximum wavelength as indicated in
Table 1. All experiments were performed in triplicates; controls did not contain laccase. The incubation was carried out for all experiments at 30 °C for 4 hours. pH of the tartrate buffer, dye concentration, HBT concentration, and enzyme concentration were the independent variables parameters that were optimized in this study. The decolorization was calculated as (Equation (1)):
where Absorbance t
0 is the absorbance of the reaction mixture at the maximum wavelength of the dye before incubation with the enzyme and Absorbance t
f is the absorbance of the reaction mixture after incubation 1.5 mL reaction mixture in 100 mM tartrate buffer. The decolorization rate was determined as followed (Equation (2)):
2.6. Box–Behnken design
As mentioned before, the Sirius grey decolorization yield in percentage (designated as y1), and its decolorization rate in percentage of color removal per minute (designated y2) were considered as the experimental studied responses (Equation (1) and (2), respectively).
The aim of this study was the determination of the diverse influences of the four reaction studied factors namely: HBT concentration (mM), pH, initial dye concentration (mg/L), and initial enzyme concentration (U/mL), on the studied responses: Decolorization yield (%) and rate (%/min).
Table 3 presents the values of coded and uncoded levels used in this work.
Since each response can be influenced by one or more factors, the multivariate study using the Response Surface Methodology (RSM), with Box–Behnken design, can be beneficial. In fact, the use of this methodology provides the possibility to determine the best polynomial multivariable model involving the coefficients calculation and statistical tests on one hand, and the accurate identification of the optimum responses and the carrying out of the relative conditions on the other hand.
27 tested runs, repeated in triplicates, were analyzed via the Box–Behnken design and presented in
Table 4 experimental responses values as function of experimental conditions.
The adopted model in this case had the following form with four studied factors (Equation (3)):

where ŷk are the modeled studied responses: decolorization yield of Sirius grey after 4 hours in % (k = 1), and decolorization rate in % of color removal/min (k = 2); β0, βi, βii and βij are the model’s intercepts, linear, quadratic and interactions coefficients, respectively; xi is the coded level of variable factors, n is the number of factors (n = 4)
2.7. Design of experiments and statistical analysis
The experimental Design, the model’s coefficients determinations, statistical analysis of the model quality and of the different factors’ influences, figures drawings, and the optimization protocol were carried out within Minitab® 19.2020.1 Statistical Software (64-bit) (© 2020 Minitab, LLC All rights reserved). The models’ coefficients were determined using the least-squares method. Analysis of variance (ANOVA) was used to identify the level of significance of the studied model and the tested factors and their interactions with a confidence level at 95% (p < 0.05). The coefficient of determination (R2), the adjusted coefficient of determination (R2adj), and root mean square error (RMSE) were chosen to quantify the model fitting quality.
2.8. Phytotoxicity assay
A phytotoxicity assay of the treated and untreated dye solutions was carried out using radish seeds (
Raphanus sativus). A Whatman filter paper was initially soaked with 2 mL of sterile distilled water then 5 mL of treated and untreated dye solution were poured. Ten seeds were distributed on the paper and the Petri dishes were incubated in the dark at 22 °C for 7 days. Germination index (GI) was calculated according to Equation (4):
4. Discussion
Strain BS9 shares more than 99% similarity of its ITS1-ITS4 region of rDNA with members of the species
C. gallica. This fungus is known for its capacity to grow on several woods to produce ligninolytic enzymes [
7,
40]. Based on this fact,
C. gallica was also shown to degrade many pollutants including dyes [
7,
25], hydrocarbons [
41], phenols [
42] and bisphenol A [
43]. Recently,
C. gallica was shown to be able to degrade antibiotics [
44]. In the same work a proteomic analysis showed the presence of one major secreted laccase, although the presence of several laccase genes in the genome.
Decolorization of a wide range of synthetic and textile dyes using laccases from basidiomycetes has been investigated in recent years [
45]. For this reason, we used
C. gallica for the decolorization of Sirius grey; the latter belongs to the azo compounds that contain one or more azo groups (N=N) and most of them are xenobiotics [
46]. By using the culture filtrate from
C. gallica, 48% decolorization of Sirius grey was achieved. This decolorization yield was improved by adding 1 mM of HBT to the reaction mixture. Under these conditions, the decolorization yield increased to 81%. Therefore, the use of a mediator (such as HBT) is necessary, especially for certain laccase with low redox potential or in case the substrate is highly recalcitrant. Indeed, laccase mediators are low molecular weight molecules with a significant redox potential, enabling it to act as an electron messenger between the substrate and the enzyme [
47,
48].
The experimental design that was performed, used four variable factors namely: initial enzyme concentration, initial dye concentration, initial HBT concentration and pH. Ben Ayed et al. [
7] found that these factors had significant effect on Reactive black 5 (RB5) decolorization using a laccase-like activity of cell free supernatant from
C. gallica. The optimized conditions obtained for laccase concentration, HBT concentration and pH were 1 U/ml, 50 mg/L, 1 mM and pH 5, respectively with a maximum decolorization yield of 87%.
In this study, the crude laccase of
C. gallica was used in the decolorization experiments, and the presence of 4 azo groups in the Sirius grey, makes its treatment more challenging. The decolorization rate obtained here (87% after 4 hours) is significant when compared to the results reported by Daassi et al. [
25]. In their study, they used partially purified
C. gallica laccase for the treatment of three different groups of dyes. However, the RB5 and Bismarck brown R (BBR), which are diazoic dyes, did not show significant decolorization rate. Concerning BBR, the rate was approximately 47.1% over 24 hours, whereas for RB5, this rate did not exceed 70%, even after 24 hours of incubation in the presence of 1 mM HBT.
To identify the influence of the interactions between studied factors on decolorization yield and rate, 3D-surface responses were designed (
Figure 4). Increasing the HBT concentration to its highest level (1 mM) followed by the increase of pH to 4-5 resulted in enhancing decolorization yield to a level of 80% (
Figure 4a). This effect was observed in the interactions between pH × dye concentration and pH × enzyme concentration (
Figure 4d,e, respectively). However, increasing the pH beyond 5 led to a reduction of decolorization yield and this aligns with Forootanfar et al. [
49] observations. These findings were explained by the fact that hydroxide anions could bind to the enzyme at acid pH and this affects negatively the electrons transfer. In contrast, Aksu and Tezer [
50] considered the possibility of basic azo dyes to become charged positively at higher pH and this affects their interactions with the mediator and enzyme. The effect of this interaction on the decolorization rate showed that the highest rate was achieved at pH 3, independently of the other factors concentration (
Figure 6a,d,e). This can be explained by the fact that pH 3 matches with the optimum pH for the used enzyme, allowing the decolorization to reach a high speed. This aligns with previous studies claiming that fungal laccases are active at 3-5 pH range [
7,
51]. In addition, the stability of the dye can be affected at pH 3 increasing even more the rate of decolorization as reported before by Yin et al. [
52] and Ben Ayed et al. [
7].
Figure 4b shows the interaction HBT concentration × Dye concentration at pH 4.5 and 0.6 U/mL of Enzyme. Increasing the HBT concentration to its maximum level at different dye concentrations resulted in an increase in the decolorization rate from values less than 1.80%/min (for 0.2 mM of HBT and 150 mg/L of Dye) to values more than 2.25% (for 1 mM of HBT and 50, 100 and 150 mg/L of Dye). This means that higher HBT concentrations improved the dye oxidation by laccase. However, at 150 mg/L dye concentration (and 1 mM of HBT), the decolorization yield decreased slightly and this was likely due to enzyme and HBT saturation. As a matter of fact, excessive dye concentration can lead to enzyme inhibition and/or unproductive reactions by intensifying competition for the enzyme’s active sites and substrate saturation. This hypothesis about enzyme-substrate concentration was discussed by Benzina et al. [
53]. The effect of saturation can also be seen in the Enzyme concentration × HBT concentration interaction (
Figure 4c). In fact, increasing HBT concentration up to 0.8 mM and that of the enzyme to 0.8 U/mL could led to an increase in the decolorization yield reaching 85%, but a higher concentration of enzyme and HBT led to a slight decrease attaining 80%. Accordingly, the presence of the enzyme and the mediator facilitates an efficient dye cleavage.
The effect of the interaction Dye concentration × Enzyme concentration on the decolorization yield is depicted in
Figure 4f. it can be see that treating high dye concentrations at low enzyme levels (and at HBT concentration of 0.6 mM) might cause a saturation effect of the enzyme and therefore of the decolororization yield. However, when the substrate concentration was decreased while increasing that of the enzyme, the decolorization yield was enhanced. This suggests that higher enzyme concentrations provide more active sites for the degradation of of the dye.
With regard to the decolorization rate,
Figure 6 represents the effects of factors interactions.
Figure 6b illustrates the interaction between HBT and dye concentrations. The lowest decolorization rate (1.8%/min) was observed for a minimal HBT concentration (0.2 mM) and a maximum concentration of Sirius grey concentration (150 mg/L). Conversely, a high response (2.2%/min) was obtained with a maximum mediator level (1 mM of HBT) and at different dye concentrations, or at low concentrations of both dye and HBT. These results were expected since they, aligning with the mediator role in improving electron transfer between the enzyme and substrate, and enhancing the treatment rate [
54,
55]. In fact, higher dye concentrations necessitate greater mediator concentrations to accelerate the decolorization rate. The need of mediator was also shown in the interaction HBT concentration × Enzyme concentration (
Figure 6c). Indeed, this interaction showed that a higher response of approximately 3%/min was reached at higher concentrations of both factors (1 mM of HBT and 1 U/mL of Enzyme). So, as the levels of the enzyme and mediator was increased, the rate of dye decolorization was accelerated. The Dye concentration × Enzyme concentration interaction presented in
Figure 6f exhibited linearity, showcasing a high response (2.8%/min) at 1 U/mL of laccase for various dye concentrations, maintaining a pH of 4.5 and 0.6 mM of HBT. Increasing enzyme concentration means increasing the active site number which boosts the decolorization rate [
56].
As textile industry effluents could be used for the irrigation of some crops [
57,
58], it is necessary to evaluate their phytotoxicity. In the present study, the toxicity of the treated and untreated Sirius grey solution was evaluated by measuring the germination index of radish seeds. It was found that the Germination Index (%IG) was significantly increased after treatment of Sirius grey by the supernatant of
C. gallica compared to the dye solution. This indicates that the treatment with laccase has effectively minimized the toxicity of the dye to lower levels compared to that of untreated dye.