Bioremoval of heavy metals from aqueous solution using dead biomass of indigenous fungi derived from fertilizer industry effluents: isotherm models evaluation and batch optimization

The present work investigated the utilization of dead biomass of the highly multi-heavy metals tolerant indigenous fungal strain NRCA8 isolated from the mycobiome of fertilizer industry effluents that containing multiple heavy metal ions at high levels to remove Pb2+, Ni2+, Zn2+, and Mn2+ as multiple solutes from multi-metals aqueous solutions for the first time. Based on morphotype, lipotype and genotype characteristics, NRCA8 was identified as Cladosporium sp. NRCA8. The optimal conditions for the bioremoval procedure in the batch system were pH 5.5 for maximum removal (91.30%, 43.25%, and 41.50%) of Pb2+, Zn2+ and Mn2+ but pH 6.0 supported the maximum bioremoval and uptake of Ni2+ (51.60% and 2.42 mg/g) by NRCA8 dead biomass from the multi-metals aqueous solution, respectively. The 30 min run time supported the highest removal efficiency and uptake capacity of all heavy metals under study. Moreover, the equilibrium between the sorbent NRCA8 fungal biomass and sorbates Ni2+, Pb2+ and Zn2+ was attained after increasing the dead biomass dose to 5.0 g/L. Dead NRCA8 biomass was described by scanning electron microscopy, energy-dispersive X-ray spectroscopy and Fourier transform infrared spectrometer before and after biosorption of Pb2+, Ni2+, Zn2+ and Mn2+ under multiple metals system. The Langmuir, Freundlich and Dubinin-Kaganer-Radushkevich isotherms were applied to characterize the adsorption equilibrium between Pb2+, Ni2+, Mn2+ and Zn2+ and the adsorbent NRCA8. By comparing the obtained coefficient of regression (R2) by Freundlich (0.997, 0.723, 0.999, and 0.917), Langmiur (0.974, 0.999, 0.974, and 0.911) and Dubinin-Radushkevich (0.9995, 0.756, 0.9996 and 0.900) isotherms values for Pb2+, Zn2+, Ni2+ and Mn2+ adsorption, respectively, it was found that the isotherms are proper in their own merits in characterization the possible of NRCA8 for removal of Pb2+, Zn2+, Ni2+ and Mn2+. DKR isotherm is the best for Pb2+ and Ni2+ (0.9995 and 0.9996) while Langmiur isotherm giving a good fit to the Zn2+ sorption (0.9990) as well as Freundlich isotherm giving a good fit to the Mn2+ sorption (0.9170). The efficiencies of Cladosporium sp. NRCA8 dead biomass for bioremoval of heavy metals from real wastewater under the optimized conditions were Pb2+, Ag+, Mn2+, Zn2+ and Al3+ ˃ Ni2+ ˃ Cr6+ ˃ Co2+ ˃ Fe3+ ˃ Cu2+ ˃ Cd2+. Dead NRCA8 biomass showed efficient ability to adsorb and reduce harmful components in the industrial effluents to a level acceptable for discharge into the environment.

Several technologies have been developed to remove heavy metals from aqueous systems (Sumalatha et al. 2023a).Among the various physical, chemical and biological procedures used to treat industrial effluents loaded with metals them bioadsorption has been documented as a promising biotechnology (El-Gendy and El-Bondkly 2016).Biosorption can be defined as a fast process independent of energy on which biological materials or biopolymers acting as sorbents to remove pollutants, such as heavy metals from wastewater through metabolically mediated or physico-chemical pathways of uptake (Sumalatha et al. 2022a, b;Khan et al. 2022).Biosorption has many advantages, such as low cost removing contaminants even in dilute concentrations, using of biomass for removal of heavy metals, cheaper production of biomass (live and dead bacteria or fungi) can be used as biosorbents for the process of biosorption, multiple heavy metals uptake at a time, treatment of large volume of effluents, no necessity for chemical additions as extremely selective for uptake and removal of specific metals and it have great potential to be an economic method for heavy metals removal from industrial effluent (Khan et al. 2022;El-Gendy et al. 2011, 2017a, b).Bioadsorption of adsorbates from the aqueous solutions to the surface of the fungal biomass as green sorbent has advantages as small handling time, slight space requisite, low energy and chemical consumptions, low-cost and effective environmentally friendly technolog (Alzahrani et al. 2017;Alzahrani and El-Gendy 2019).
In dematiaceous hyphomycetes, Cladosporium is the largest of the genera, and its species have been described as being among the most common fungi in both indoor and outdoor environments (Becchimanzi et al. 2021;El-Gendy et al. 2017b).Previous works proved that the dead fungal biomass could be an efficient option to reduce various toxic compounds and heavy metals from industrial wastewater with great adsorption capacity because it's big specific external area and effective surface interaction along with its ability endure temperature difference (Paria et al. 2022;Ab Rhaman et al. 2022).Hence, the goal of this search was to investigation the potential applies of dead fungal biomass as adsorbent for Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ from aqueous solutions, optimization of the batch adsorption operating conditions under multi-metals system including pH, adsorption time, biomass dosage and initial metals ions concentration.Furthermore, the Langmuir, Freundlich and Dubinin-Kaganer-Radushkevich (DKR) isotherm systems were applied to appropriate the equilibrium isotherm of these heavy metals on the surface of the fungal strain NRCA8 under the optimized conditions.

Metal ions solution and factory effluents preparation
Multi-metals ions solution stock composed of Pb 2+ [Pb (CH 3 COO) 2 ], Ni 2+ (NiCl 2 .6H 2 O), Mn 2+ (MnSO 4 •H 2 O), and Zn 2+ [Zn(CH 3 CO 2 ) 2 ] was prepare by dissolving a proper amount of each metal in deionized water to prepare a concentration of 1000 mg/L.The preferred concentrations of C o = 10, 25, 50, and 100 mg/L of each ion in the multi-metals solution were prepared by dilution of the stock solution.pH was from 2.0 to 6.0 by 0.1 M HCl and 0.1 M NaOH solutions.Real wastewater samples belonging to fertilizer industry were collected from the drainage areas of fertilizer industry wastewater at the Manqabad industrial regions of Assiut, Egypt in poly ethylene bottles of 1 L. The wastewater samples were gathered, filtered and divided into three portions.The first was processed immediately for the isolation of their fungal mycobiome; the second was located in sterile flasks including 2.5 mL nitric acid and stay in at 4 °C until analyses for their characteristics before treatment within 24 h of collection by using Agilent 5100 Synchronous Vertical Dual View (SVDV) ICP-OES, with Agilent Vapor Generation Accessory VGA 77.The third kept in − 80 °C until treatment with the dead biomass of selected fungus based on the optimization experiments followed by analyses for its Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ contents and other properties.

Isolation of mycobiome from the fertilizer industrial effluents
The fungal biosorbents were isolated from fertilizer manufacturing industrial wastewater (drainage areas, Manqabad, Assiut, Egypt).Contaminated water drainage samples were filtered, serially diluted using the serial dilution technique and inoculated into potato dextrose agar (PDA) medium.The plates were incubated for 10 days at 28 °C.The hyper multi-metals tolerant isolate NRCA8 was selected, identified and analyzed for its biosorption efficiency and uptake capacity of different heavy metals from the multimetals aqueous solutions as previously described.More over the ability of this strain to enhance the industrial wastewater properties was evaluated under the conditions optimized earlier in the batch process including pH, contact time, biomass dosages and initial heavy metals concentrations.
Screening of fungal strains for the higher multi-metals tolerant strain under single and multi-metals types Different concentrations of Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ were prepared individually and in different combinations of binary, tertiary and quartet of these metals in a series increasing concentrations of these heavy metals from 0.1 to 6.0 g/L in equal mass ratios.These concentrations of each metal added to the Czapek yeast extract agar (CYEA) medium, individually to investigate the resistance of fungal mycobiome of fertilizer industrial wastewater and selected the hyper multi-metals tolerant strain.All experiments were conducted in duplicate.The inoculated plates were incubated at 30 °C for 7 days to confirm their growth at these varied metals concentrations individually and in mixtures.Growth test was conducted and calculated the higher tolerance concentration for each treatment.The strain which showed the higher tolerance concentrations under single, binary and multi-metals conditions was chosen for the further study.
Description and identification of the excessive tolerant isolate NRCA8 Description and identification of the strain under study NRCA8 was performed based on its morphotype, lipotype and genotype as earlier described (Crous et al. 2009;Domsch et al. 1980;Kujur and Patel 2014;Ogórek et al. 2012;Schubert et al. 2007;Stahl and Klug 1996;Weete 1980).The fungus morphology and colony characteristics on PDA, malt extract agar (MEA) and oat agar (OA) media (Difco Laboratories, USA) were recorded and photographed after 14 days of incubation at 25 °C.The surface and reverse colors were rated using the charts of Rayner (1970).For micro-morphological observation, preparations of colonies made on PDA were mounted onto Shear's solution and conidial development and branching patterns were studies (Crous et al. 2009;Schubert et al. 2007).Temperature range for the growth of fungal isolate was established on PDA cultures after 14 days at temperatures ranging from 5 to 40 °C at intervals of 5 °C.Also pH range for growth was determined by culturing NRCA8 at pH ranged from 3.0 to 11.0 on PDA.

Molecular identification of the hyper tolerant isolate NRCA8
DNA was extracted and purified from culture NRCA8 growing on PDA after 7 days of incubation at 28 °C, the DNA of isolate used as the template for the PCR, for isolate NRCA8, partial sequence of the rDNA with primers ITS1 and ITS4 was amplified for the internal transcribed spacer (ITS) region of rDNA, amplification, PCR product purification and sequencing in both directions were performed following the previous protocol (El-Bondkly 2012;El-Gendy et al. 2018;El-Bondkly and El-Gendy 2022;White et al. 1990).Finally, the software SeqMan (DNAStarLasergene) was applied.The sequence achieved was matched with other fungal sequences put in the NCBI database using the BLASTn.Phylogenetic reconstruction was made with the phylogenetic marker (ITS) approved for perfect identification at the species level using maximum likelihood (ML) analyses, with the MEGA11 software (Kumar et al. 2018;Tamura et al. 2011Tamura et al. , 2021)).

The dead fungal biomass (adsorbent) preparation
The dead biomass of the NRCA8 isolate was used as natural biosorbent for Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ from aqueous solutions under multi-metals system.Ten-days old culture spores (10 6 spore/mL) were transferred individually into 500 mL Erlenmeyer flasks, each including 100 mL potato dextrose broth medium and incubated at 30 °C, 150 rpm on a rotary shaker for 10 days.The biomass of isolate was pelletized by filtration using filter papers (Whatman No. 1), washed with 0.1 M NaCl followed by deionized water.Dead biomass was obtained by autoclaving, washed with 0.1 M NaCl, pre weighted and dried in an oven at 60 °C followed by crushing to a fine powder and kept in sterile polyethylene bottles at 4 °C until use.

Evaluation of adsorption performance of the selected fungal isolate NRCA8
The biosorption trials were conducted under multimetals system in quick-fit flasks containing biosorbent dosage 1 g/L of the dead biomass in working volume 50 mL aliquots of a mixture of Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ at a concentration of 100 mg/L (25 mg/L for each).Flasks were kept on rotary shakers (150 rpm) at 30 °C and pH 5.0 for 30 min.The samples were filtrated and the concentration of each metal ion in the multi-metals solution was estimated.All samples were digested using Anton-Paar microwave digestion method (APHA 2017).The supernatants were analyzed for residual heavy metals.Heavy metal solutions without biomass were assisted as control, trials were conducted in duplicate and average values were calculated.The data was recorded in percentage using the following equations: the metal removal efficiency percentage (R%) was given in the multi-metals solution of Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ or in real wastewater according to the equation (Li et al. 2009;Singh et al. 2007): where, R is the biosorption efficiencies percentage (%) of each metal separately; C o is the initial metal concentration (mg/L), and C f is the equilibrium or final concentration for each metal calculated separately.Moreover, the adsorption capacity of fungal biomass was assessed by the equation: wherever q t is adsorption capacity (mg/g) and C t is metal concentration (mg/L) at time = t (min), V is the solution volume (L), and m is adsorbent mass (g).

Optimization of the biosorption batch factors
Batch trials were carried out in 250 mL Erlenmeyer flasks containing multi-metals solution of Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ solutes, following optimization process.The impact of pH on metal biosorption by NRCA8 biomass was evaluated by varying the initial solution pH (2, 4, 5, 5.5, and 6) by dilute HCl or NaOH, the pH values of the experiment are not set higher than 6 for avoiding any metal precipitation as hydroxide as well as the effect of the adsorption time was evaluated at various contact times (10, 20, 30, 60, 90, and 180 min) at the optimum pH.Moreover, the effects of the dead fungal biomass dosage (1, 2, 5, and 10 mg/mL) was evaluated at the optimum pH and contact time.The initial concentration of Pb 2+ , Ni 2+ , Mn 2+ , and Zn 2+ metals (C o = 10, 25, 50, and 100 mg/L of each metal in the mixture) were evaluated at the proceed optimum conditions.In each experiment, flasks were permitted to equilibrate on a rotating shaker and samples were collected after the appropriate period, aqueous solutions were filtered and each filtrat analyzed for residual concentration of each metal as well as its biosorption efficiency (%) and uptake (mg/g) were quantified to determine the optimum process parameters for maximum metal ion biosorption.

Adsorption isotherm
The uptake of heavy metal ions by the fungal biomass was evaluated by the various adsorption isotherm systems containing Langmuir, Freundlich and Dubinin-Kaganer-Radushkevich (DKR) models.
The linear form of Langmuir equation is characterized by the following equation: where Q e is the equilibrium uptake capacity (mg/g), C e is the concentration of adsorbate molecule remaining in solution at equilibrium (mg/L), Q max is the maximum ions uptake per unit mass of fungi, (mg/g) related to adsorption capacity that represents monolayer coverage and K L is the Langmuir constant equivalent to the enthalpy of adsorption (L/mg).Therefore, the linear plot of Ce/qe versus Ce gives a straight line of slope 1/qmax and intercepts 1/ (qmaxK L ) (Dąbrowski 2001).The Langmuir system were applied to determine the separation factor RL, as stated by Eq. (4) according to (Fawzy et al. 2018).
The Freundlich isotherm system suggests the heterogeneous adsorption of the surface that has unequal available sites.The linear equation can be written as follows: Q e (mg/g) is the amount of metal ion adsorbed on adsorbent at equilibrium, C e (mg/L) the equilibrium concentration of metal ion in the solution, K f (mg 1−1/n L 1/n g −1 ) is a Freundlich isotherm constant describing the adsorption capacity and n is empirical parameter related with multiple layer coverage (Ayawei et al. 2017;Hamdaoui and Naffrechoux 2007).
The results were also fitted with DKR isotherm model to estimate the nature of sorption process as chemical or physical and estimate the mean energy of sorption.The linear equation of DKR isotherm is where q e is the number of metal ion adsorbed per unit weight of adsorbent (mol/g), q m is the maximum sorption capacity, b is the activity coefficient related to mean sorption energy, and Ɛ is the Polanyi potential, which is equal to: where R is the gas constant (kJ/kmol K) and T is the temperature (K).By plotting a relationship between lnq e and ε 2 , β and q DR can be obtained.(D-R) isotherm parameter β used to determine adsorption energy E (KJ/mol) as follows: Scanning electron microscope (SEM), energy-dispersive X-ray spectroscopy (EDX) and Fourier transform infrared spectroscopy (FTIR) analyses SEM-EDX method was applied to define the chemical description of fungal biomass before and after adsorption of heavy metal ions Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ .The fungal biomass amended with multi-metals solution of these ions at initial concentration of 100 mg/L of each was used for SEM analysis (SEM Quanta FEG 250 with field emission gun, FEI Company-Netherlands) at the Central Laboratory (National Research Centre, Egypt).The functional groups on the biomass surface were determined by a Fourier transform infrared spectrometer (Broker Vertex80v, Germany).Confirmation of presence of metal ions on the fungal biomass surface was tested using EDX analysis using an X-ray micro-analyzer connected to a scanning electron microscope.The individual ratios given represent the average of ten measurements.In FT-IR were used to measure the transmittance spectra recorded in the range of 4000-400 cm −1 with resolution 4 cm −1 at the Central Laboratory of National Research Centre, Egypt to define the vibration frequency groups in the biosorbent NRCA8 before and after biosorption of heavy metals under study from the multi-metals solution.

Wastewater analyses
For the analyses of wastewater, samples collected from the fertilizer industry effluents, Mankapad, Assuit, Egypt were exposed to centrifugation at 2000 rpm for 2 min, filtration by Whatman filter paper with 0.2 μm pore size, acid digestion according to APHA (2017) followed by determination the initial Vol:. ( 1234567890) and final concentrations of metal ions concentrations before and after treatment with NRCA8.Moreover, the other parameters including total suspended solids (TSS), total dissolved solids (TDS), oil and grease, chemical oxygen demand (COD), nitrogen (N), and phosphorus (P) were evaluated following the standard approaches for the analysis of water and wastewater (APHA 2017).

Statistical analysis
The data was statistically processed by analyzes of variance (ANOVA), followed by Tukey's tests when significant effects were detected (P ≤ 0.05).Data were expressed as means ± standard error.

Selection of the high metals tolerant fungal strain
The mycobiome derived from fertilizer industrial effluents were evaluated and screened for their tolerance of Ni 2+ , Pb 2+ , Mn 2+ and Zn 2+ in single and multimetals systems, among them the strain under the isolation code NRCA8 showed the highest metals tolerant behavior against these heavy metals (Table 1).Data in Table 1 indicated that the highest concentrations of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ that the strain was able to grow were 5.0, 2.5, 3.2, and 2.6 g/L in the presence of single metal in growth medium while it was estimated at 5.4, 6.0, 5.5, 2.7, 2.5, and 2.7 g/L in the presence of the binary combinations (Pb 2+ + Ni 2+ ), (Pb 2+ + Zn 2+ ), (Pb 2+ + Mn 2+ ), (Ni 2+ + Zn 2+ ), (Ni 2+ + Mn 2+ ) and (Zn 2+ + Mn 2+ ), respectively (Table 1).Moreover, in the multi-metals system in the growth medium composed of (Pb 2+ + Ni 2+ + Zn 2+ ), (Pb 2+ + Ni 2+ + Mn 2+ ), (Pb 2+ + Zn 2+ + Mn 2+ ), (Ni 2+ + Mn 2+ + Zn 2+ ) and (Pb 2+ + Ni 2+ + Zn 2+ + Mn 2+ ), individually the highest initial multimetals concentration of each ion the strain was able to grow at 4.1, 4.0, 5.0, 3.45, and 4.0 g/L, respectively (Table 1).Then the fungal strain that showed high adaptive tolerance to Ni 2+ , Pb 2+ , Mn 2+ and Zn 2+ under single, binary, ternary and quaternary metals systems was then selected for further studies.Strain NRCA8 has created exceptional performance and can reduce concentrations of Ni 2+ , Pb 2+ , Mn 2+ and Zn 2+ well below the irrigation threshold level set by the Food and Agriculture Organization (FAO).In agreement with our results Madhuri et al. (2022) who stated that fungi are help maintain tolerance to heavy metals in various contaminated sites by developing different methods of resistance against different heavy metals, and have the potential to survive through adapting or mutating at high concentrations of heavy metals as well as can also decrease heavy metals from environment to some extent, and thus these approaches would facilitate the development of enhanced methods for the bioremediation of heavy metals in the environment.Furthermore, Dey et al. (2016) reported that the fungal isolates A. terreus AML02, Paecilomyces fumosoroseus 4099, Beauveria bassiana 4580, A. terreus PD-17, and A. fumigatus PD-18 exposed to a mixture of multiple metals (Cd, Cr, Cu, Ni, Pb and Zn) of different concentrations (6, 12, 18, and 30 mg/L) showed high metals tolerance index for each metal in the mixture suggesting their better adaptability to multi-metals stress.
Characterization and identification of highest tolerant isolate NRCA8 The reproductive structures of the selected strain NRCA8 after 14 days of growth on different media are shown in Table 2 and Fig. 1a-e.The NRCA8 isolate on PDA formed distinct microscopic features including a short cylindrical mycelium, 1.5-4.5 μm wide, olivaceous-brown, smooth or verruculose concerning the base of conidiophores, thick-walls.
On the other hand, isolate NRCA8 included fatty acids of carbon chain lengths extending from 14 to 22.The predominant fatty acids detected in strain NRCA8 were 16:0 (24.61 ± 1.39%), 18:2 (24.20 ± 1.55%), 18:1 (22.37 ± 1.45%) and 22:6 (10.00 ± 0.84%), respectively followed by 18:1 ω9c (6.24 ± 0.51%), 16:1 (4.10 ± 0.36%), 20:4 (3.40 ± 0.39%) and 14:1 (3.19 ± 0.28%) (Table 2), which propose their prospective structural role in cell membranes as earlier described (El-Gendy et al. 2017b).However, fatty acids that characterized less than 1% of the total fatty acid content were 14:0 (1.22 ± 0.16%) and 18:0 (0.67 ± 0.04%) (Table 2).Depending on the morphological, microscopic, physiological and biochemical features, strain NRCA8 was defined as a strain belonging to Ascomycota, family; Davidiellaceae and genus; Cladosporium.Earlier investigates described fatty acid types as a reliable method to illustrate fungal species and strains (Kujur and Patel 2014;Stahl and Klug 1996).The amplification of the PCR product of the ITS regions of the rDNA of isolated strain NRCA8, generated fragments of 530 bp by ITS1 and ITS4 primers as well as presented to GenBank (accession no.ON667856).The results of the comparisons with Blastn they correspond to a local alignment leading to putative identifications.Therefore, to have greater reliability in the identification, phylogenetic analyzes were performed.A relative analysis by Blastn established that ITS area sequence from isolate NRCA8 had an important identity to the genus Cladosporium (Fig. 2).Comparison of isolate NRCA8 with the sequences of reference species in the bank database displayed that isolate NRCA8 exhibited a relationship of 97.90% with Cladosporium sp.LF183, and 97.72% with Cladosporium sp.NP23-10-5, and Cladosporium sp.BB23-3-1 (Fig. 2).Phylogenetic analyzes based on 530 bp were performed by MEGA11, the phylogenetic tree achieved by using the maximum likelihood system is explained in Fig. 2. According to the sequence analysis of the ITS area, along with its phenotypic and chemotypic features, isolate NRCA8 was recognized as Cladosporium species and labelled as Cladosporium sp.NRCA8.The most prominent markers of fungal phylogenetic are the ITS areas of the rDNA sequences.Previously, the fungal strains NRCF5, Gen 9, Gen 20, and ALAA-20 were identified as Aspergillus sp.NRCF5, Trichoderma sp.Gen 9, Cladosporium sp.Gen 20, and Fusarium sp.ALAA-20, respectively based their phenotypic and chemotypic features together with ITS sequences analysis (El-Bondkly 2012;El-Bondkly et al. 2021;El-Gendy et al. 2017b).

Optimization of operating conditions for different heavy metals removal in multi-metals aqueous solution
The pH is significant factor that acting an important role in the removal of heavy metals by fungi since it impacts the the metals speciation in the solution and the surface characteristics of the fungi (El-Gendy et al. 2017a).As revealed in Fig. 3, an increase in the pH from 2.0 to 5.5 significantly improved the removal efficiency of Pb 2+ , Zn 2+ and Mn 2+ from 61.6, 21.6%, and 17.6 to 91.30%, 43.25, and 41.50%, respectively as well as the uptake of Pb 2+ , Zn 2+ and Mn 2+ (5.86, 3.78, and 1.7 mg/g) was attained at pH 5.5, respectively (Fig. 3).This may be attributed to that in a highly acidic environment (pH 2.0), the adsorption locates of NRCA8 became saturated with a positively charged hydrogen ion (H + ), which could compete with Pb 2+ , Ni 2+ , Mn 2+ and Zn 2+ for the active sites that resulted in decrease the heavy metals under study in their aqueous mixture (Fig. 3).With increasing pH, negatively charged OH − tends to dominate the adsorption sites (El-Gendy and El-Bondkly 2016; Fei and Hu 2022).Furthermore, the results in Fig. 3 show that the multi-metals sorption efficiency of Pb 2+ , Zn 2+ and Mn 2+ onto NRCA8 doesn't improved by rising pH value to pH 6.0 and a slightly increase in Ni 2+ removal (51.60%).Zhang et al. (2020) reported that low pH resulted a negative effect on the adsorption performance because it can accelerate the dissolution/ oxidation of functional groups on the surface of the adsorbent, which leads to the release of ions and desorption of the entrapped heavy metals (Table 3).
Contact time is just one the most important influences affecting biosorption efficiency.The rapid uptake of Mn 2+ , Ni 2+ , Pb 2+ and Zn 2+ on NRCA8 (1.70, 2.15, 5.86 and 3.78 mg/g) along with the highest removal efficiency (28.81%, 40.86%, 91.56% and 32.17), respectively were achieved and equilibrium was practically reached after 30 min contact time.The rapid adsorption capability through the initial phase is likely because of the abundance of vacant active locations on the NRCA8 dead biomass and the high concentration gradient of solutes Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ in the multi-metals solution as previously reported for the nickel, zinc and mercury biosorption by dead biomass of metal tolerant fungi (Alzahrani et al. 2017;Alzahrani and El-Gendy 2019).There was no significant change in equilibrium concentration for Pb 2+ , Ni 2+ and Zn 2+ ions after 30 min but only a slight decrease for Mn 2+ ion was observed, these decreases may be attributed to the agglomeration of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ onto the NRCA8 active sites, the difficulty of occupying the remaining binding locations because of forces concerning the solute molecules that solid and bulk stages and  permanent interaction can occur.Then Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ ions successfully diffuse from the boundary layer surrounding the NRCA8 particles to the bulk solution as previously reported (Alothman et al. 2020;Chen et al. 2019;Kumar et al. 2022).
In this research, it was detected that the equilibrium between the sorbent NRCA8 biomass and sorbates Ni 2+ , and Pb 2+ in the multimetals solution was achieved after increasing the biomass dose to 5.0 g/L, where their bioremoval efficiency reached their maximum values (51.6%, and 91.72%, respectively); afterwards a non-significant decrease in the removal capacity was observed (Table 4).Additionally, the maximum adsorption of Mn 2+ and Zn 2+ (47.46% and 43.25%, respectively) were accomplished by further increasing the biosorbent does to 10.0 g/L from the multi-metals aqueous solution.The increased bioremoval efficiency with increasing doses of NRCA8 could be attributed to the larger surface zone of the biosorbent, which in turn increased the availability of active locations for metal ions (Ayele et al. 2021;El-Gendy et al. 2011).
The bioremoval (%) of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ were evaluated under different initial multi-metals concentrations extending from 10 to 100 mg/L at the optimal operating circumstances for each metal in the multi-metals aqueous solution as shown in Table 5.The improved biosorption at first phases can be accredited to the larger driving force of metal ions into the fungal surface and abundance of vacant binding locations on the biosorbent surface (Alzahrani and El-Gendy 2019;Kumar et al. 2019).The multi-sorption of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ by NRCA8 dead biomass from the solution was significantly decreased with increasing the early metal concentrations of each metal in the mixture (Table 5), the reduction in heavy metals bioremoval can be associated to the exaggerated quantity of ions about saturation of all available sorption locations on the NRCA8 fungal biomass surface and accomplishment of equilibrium among the sorbent and sorbate, thus averting further adsorption of heavy metal ions (Tu et al. 2018).

Evaluation of adsorption isotherm models
Data analysis by isothermal equilibrium adsorption is important for understanding the mechanism of metal ions adsorption into NRCA8.In the current work, the adsorption data between the adsorbate ions (Pb 2+ , Ni 2+ , Mn 2+ and Zn 2+ ) and the adsorbent (dead fungal biomass) at equilibrium condition were evaluated by different models including Langmuir, Freundlich and DKR adsorption isotherms.
The Langmuir isothermal model has been extensively used in adsorption studies.It describes the deposition of adsorbates on the free surface of the adsorbent and the formation of a monolayer adsorbate on the outer surface of the adsorbent.The parameters of the Langmuir isotherm system for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ adsorption on NRCA8 showed that Q max = 18.05, 1.77, 16.89, and 4.12 mg/g as well as K = 21.54,18.36, 0.75, and 1.68 L/mg, respectively (Table 6, Fig. 4).Moreover, the determination of the coefficient of the linear equation for Langmuir was R 2 = 0.974, 0.999, 0.974, and 0.911 for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ respectively.This noticing indicated that the trial inputs pH, temperature, biomass dosage, contact time and initial metal ions concentration showed positive and linear impacts on the model output beside no over-fitting problem occurred during prediction.Therefore, the Langmuir system can adequately define the adsorption mechanism of Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ onto fungal dead biomass at equilibrium.In the current study the separation factor (RL) for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ was estimated to be 0.0019, 0.0054, 0.1177, and 0.0562, then the RL values lies in the range of 0 < RL < 1 indicating that adsorption of these heavy metals on surface NRCA8 is preferred at both low and high initial concentration as the RL values are very close to zero.The adsorption is satisfactory at 0 < RL < 1, irreversible at RL = 0, linear at RL = 1, and unfavorable at RL > 1 (Khayyun et al. 2019).Then the Langmiur system provides a good fit to the sorption procedure, which is confirmed by the positive values that gotten for the Langmiur constants presented in Table 6.Similarly, Sumalatha (2023b) proved that Langmuir model and Pseudo second order rate models best suited the adsorption process of heavy metal Cr(VI) from synthetic medium using biodegradable natural polymeric biosorbent, suggesting that ions were adsorbed in monolayer due to their chemical affinity and the thermodynamic characteristics demonstrated the process possibility, spontaneity, and exothermic nature of adsorption.
From the Freundlich model parameters listed in Table 6 and Fig. 5, n = 1.825, 5.740, 1.304, and 1.617 as well as K f = 2.235, 0.830, 0.326, and 0.207 for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ , respectively.While the regression coefficients for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ were R 2 = 0.997, 0.723, 0.999, and 0.917, respectively (Table 6).El-Gendy et al. (2017a) reported that the values of K f and n define the steepness, the isotherm curvature and the adsorption capability of the adsorbents increases at a higher K f value.Khayyun et al. (2019) and El-Gendy et al. (2011, 2017a) described that once 1/n values are in range of 0.1 < 1/n < 1, the adsorption method is desirable, irreversible at 1/n = 0, and unfavorable at 1/n > 1. Hence the 1/n determinants of 0.548, 0.174, 0.767, and 0.618 suggested that the isothermal type desirable, favorable and refer to the strong interaction between fungal biomass and Pb 2+ , Ni 2+ and Mn 2+ under multi-metals sorption condition (Table 6, Fig. 5) but it was not much desirable for Zn 2+ adsorption.In line with our results John et al. (2022) reported that the experimental findings of Pb(II) biosorption from aqueous solutions were fitted with Langmuir, Freundlich and Temkin isotherms but Freundlich model gave good fit with R 2 = 0.99 and a maximum adsorption capacity of 36.49mg/g using dried turmeric leaves powder as biosorbent.Vol:.( 1234567890) The DKR isotherm model in Table 6 and Fig. 6 was selected to assessment the characteristics porosity of the dead biomass and the apparent energy of adsorption as well as describe the equilibrium between the adsorbates and the adsorbent.The maximum sorption capacity (saturation capacity in mol/g), X m , and values describing the total specific micropore volume of the sorbent were estimated for Pb 2+ , Zn 2+ , Mn 2+ and Ni 2+ (5.14 × 10 −4 , 4.83 × 10 −5 , 1.91 × 10 −4 , and 1.19 × 10 −3 mol/g), respectively (Table 6).The parameters of the DKR model recorded in Table 6 indicated that the determination of the linear equation coefficient for DKR was R 2 = 0.9995, 0.7560, 0.9996, and 0.9000 for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ , respectively.The type of adsorption can be predicted by calculating the adsorption energy E, and free energy E for Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ (10.86, 15.83, 8.02, and  10.01 kJ/mol) were positive values, which showed the endothermic nature of the heavy metals sorption method by NRCA8 dead biomass.In the present work, the energy values for heavy metals sorption on NRCA8 biomass were between 8 and 16 kJ/ mol which indicated that the sorption process can  be construed by ion exchange as reported before (Embaby et al. 2021;Igwe et al. 2007).
Characterization of dead NRCA8 biomass before and after bioadsorption of Pb 2+ , Ni 2+ , Mn 2+ and Zn 2+ by SEM-EDX and FTIR spectroscopy The SEM-EDX and IR analyzes are appeared in Figs.7a, b and 8a, b; respectively.The SEM analysis in Fig. 7a and b showed the surface morphology of untreated and treated dead NRCA8 biomass.The unloaded NRCA8 dead biomass displayed the incidence of distinct, regular, intact long rod, cylindrical sheets or even ribbon shaped mycelial fibers which were extremely branched and tangled as well as some few uniform conical spiral like structures were observed (Fig. 7a).Furthermore, unloaded NRCA8 dead biomass included huge numbers of unoccupied pores with heterogeneous aggregation of single particles (Fig. 7a).However, the SEM image upon treatment the biomass of Cladosporium sp.NRCA8 biomass with the heavy metals Pb 2+ , Ni 2+ , Mn 2+ and Zn 2+ , from the multi-metals aqueous solution under the optimized conditions (Fig. 7b) showed blocking of the vacant sites, representing the connection of metal ions onto the surface of NRCA8 dead biomass.Conversely, significant morphological variations were detected in the hyphal shape including distorted, shrunken and irregular expansion of thick fungal mycelial mass with some broad conical deformations in the existing of various heavy metals (Fig. 7b).These changes after treatment might be because of the uptake and accumulation of metal ions lead to variations in physiological, morphological, cellular and molecular levels as well as stress conditions that cause an increase in the area of interaction of metal ions with the fungal biomass.Morphological variations in the fungal mycelia below heavy metals stress have been stated previously (Liaquat et al. 2020), which might be because of oxidation of protein and DNA molecules, variations in ultrastructure, or inhibition of antioxidant defense system in cells (Chen et al. 2014).Bankar et al. (2018) reported that a study of the dimorphism on the marine yeast Yarrowia lipolytica revealed a change in morphology as elongated, oval or rounded in reaction to various heavy metals stress containing Pb 2+ and Cd 2+ .Furthermore, the EDX was applied to know the chemical composition of NRCA8 dead biomass before and following the uptake of Pb 2+ , Ni 2+ , Mn 2+ and Zn 2+ (Fig. 7a, b).The elemental profile of EDX indicated that the C and O contents were 53.3% and 46.7% in terms of weight (Fig. 7a).It was observed that the detectable amounts of N 3+ , C 4+ , O 2+ , Pb 2+ , Mn 2+ , Ni 2+ and Zn 2+ were detected to be adsorbed on the fungal mycelia after the handling with multi-metals solution and metals blend (Fig. 7b).A noticeable peak of Pb 2+ , Mn 2+ and Zn 2+ adsorbed by the cell wall fungus was detected in the EDX.Conversely, there were no important peaks for the Ni 2+ and Mn 2+ were existing viewing an intracellular accumulation of these metals rather than binding on the cell surface (Fig. 7b).Also, laser scanning microscopy images proposed that an intracellular localization primarily within vacuoles and vesicles (Sharma et al. 2022;Traxler et al. 2022).Interestingly, after treatment with Pb 2+ , Zn 2+ , Ni 2+ and Mn 2+ ; it is evident from Fig. 7b that the amount of these ions on the surface increases significantly from 0.0% in the untreated sample to 2.44, 2.23, 0.92 and 0.60, respectively in the multimetals-treated sample.The amount of O (II) decreases from 46.7 to 26.82% while the amount of C and N was increased from 53.3 and 0.0% in untreated sample to 61.57% and 5.41%, respectively in loaded NRCA8 sample (Fig. 7b).

Conclusions
The results of this work indicate toward the onsite use dead fungal biomass of the multi-metals hyper tolerant fungus Cladosporium sp.NRCA8 isolated from the mycobiome of the fertilizer industry effluents in the multi-adsorption of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ from multi-metals aqueous solutions and industrial wastewater.The highest multi-adsorption of these heavy metals was achieved by fungus NRCA8 after optimizing operating conditions including temperature; 30 °C, pH; 5.5-6.0,interaction period; 30 min, biosorbent dose; 5.0 g/L for uptake and removal of each ion for removal and uptake capacity and agitation speed; 150 rpm, respectively.Moreover, SEM showed the cellular changes of the fungus biomass resulting from the multi-adsorption of these heavy metals.EDX confirmed the involvement of biosorption and intracellular accumulation in the removal and uptake of these ions as well as FTIR showed that the bioremoval and uptake of Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ were attributed to the functional groups of Cladosporium sp.NRCA8 biomass included hydroxyl, carbonyl, ethers, amines/amides and amino.The adsorption data between the adsorbate ions (Pb 2+ , Ni 2+ , Zn 2+ and Mn 2+ ) and the adsorbent NRCA8 biomass at equilibrium condition were estimated by various adsorption isotherm systems containing Langmuir, Freundlich and DKR.Since there is Ethical approval This article does not contain studies with human participants or animals.

Fig. 2
Fig.2Phylogenetic tree generated by the maximum likelihood based system in the rDNA sequences of the ITS region with isolate NRCA8 belonging to the genus Cladosporium obtained through with 1000 repetitions.Saccharomyces cerevisiae was used as an outgroup

Fig. 3
Fig.3Influence of pH on the bioremoval efficiency and uptake capacity of different heavy metals by the dead biomass of NRCA8 strain

Fig. 7
Fig. 7 The SEM and EDX analysis of the surface morphology of dead NRCA8 biomass (a); before and (b); after adsorption

Fig. 8
Fig.8The FTIR analysis of the unloaded (a) and loaded (b) dead NRCA8 biomass before and after adsorption ◂

Table 1
Screening of fungal isolates obtained from the microbiome of the fertilizer industrial wastewater for the tolerance concentrations of Pb 2+

Table 3
Influence of contact time on the bioremoval efficiency (%) and uptake capacity (mg/g) of metal ions (Ni 2+ , Pb 2+ , Mn 2+ and Zn 2+ ) by the dead biomass of the fungus NRCA8

Table 5
Influence of initial metal concentrations on the bioremoval efficiency (%) and uptake capacity (mg/g) of metal ions(Ni 2+

Table 7
Assesment of fertilizer industrial effluent composition before and after adsorption by NRCA8 dead biomass *All determents in mg/L, pH in pH units, temperature in °C, and turbidity in NTU