Bioremediation of Crude Oil by Rhizosphere Fungal Isolates in the Presence of Silver Nanoparticles

Background: This research work focuses on the utilization of indigenous fungi for in situ bioremediation of crude oil in the presence of silver nanoparticles. Methods: Two fungi belonging to two different genera showed promising crude oil-degrading abilities. Fungal isolates were identified based on internal transcribed spacer rDNA sequence analysis. Gas chromatography-mass spectrometry analysis of the crude oil remaining in the culture medium after seven days was performed. The response surface method (RSM) designed by Box-Behnken was used to establish a mathematical model. Inter-simple sequence repeat (ISSR) primers were used to examine the genetic variation of fungal isolates. Results: Gas chromatography-mass spectrometry (GC-MS) analysis after seven days showed that the optimum biodegradation of crude oil was 57.8%. The crude oil degradation rate was significantly affected by a temperature of 30 °C, pH value of 7, crude oil concentration of 4 g/L, a 1:1 ratio between A. flavus AF15 and T. harzianum TH07, and an silver nanoparticle (AgNP) concentration of 0.05 g. Molecular characterization in fungal isolates is extremely valuable when using ISSR markers. Conclusions: Two fungal isolates showed promising crude oil-degrading abilities with positive effect of low concentrations of AgNPs on biodegradation. RSM is an efficient mathematical method to optimize the microbial biodegradation of crude oil.


Introduction
Saudi Arabia, which has almost a quarter of the world's oil reserves, is the largest producer and exporter of oil. Petroleum is transported around the world via pipelines, roads, ships, and trains, posing a significant environmental threat in the event of spills. Crude oil is an extremely toxic carcinogenic, is strongly mutagenic, and consists of teratogenic complex compounds. It is recognized as a serious threat to ecosystems and it takes several years or decades to recover from many environmental problems after the event of a spill [1]. The presence of crude oil in the natural environment is a critical problem because it causes gradual soil degradation and occasionally leads to the permanent destruction of soil and loss of fertility [2]. This necessitates the utilization of methods that are environmentally friendly for cleaning up oil spills. One is the use of biological agents which, compared with physicochemical approaches, are better performing and more cost-effective [3]. Bioremediation is feasible, targeted, and capable of achieving high removal efficiency at a low cost.
The biodegradation of crude oil by microorganisms is a "green" alternative for treating hazardous contaminants due to their lack of environmentally degrading effects, as well as the simplicity of their low-energy design, construction, and operation [4,5]. No fungal species exist with the metabolic

Energy Dispersive Spectroscopy (EDS)
For EDS, samples were prepared on a copper substrate by the drop coating of AgNPs. Elemental analysis of single particles was carried out using a JEOL (JSM-6380 LA) equipped with a scanning electron microscope.

Crude Oil Degradation by Fungal Isolates
The ability of four fungal isolates to degrade crude oil was evaluated. First, 100 mL of mineral salt medium (MSM) placed in 250 mL flasks was supplemented with 1% crude oil, prepared according to method [20]. The pure fungal isolate (2 mL, 1 × 10 6 CFU/mL) was inoculated into the MSM. Then, the flasks were incubated at 30 • C and 150 rpm, pH 7.5, for 14 days. At the end of the incubation period, three and seven days, 100 mL of each medium was thoroughly shaken with carbon tetrachloride (50 mL, three times) in a separating funnel, and the three fractions were collected in the case of the crude oil sample. The collected organic layer was dried over anhydrous sodium sulphate. The solvent was removed using a rotary evaporator until a constant weight was achieved. The oil sample was accurately weighed, the percentage of the biodegraded oil was calculated. All treatments, except the sterile control, were performed in triplicate.

Emulsification Activity, E24
To determine the emulsification activity, an equal amount of hexadecane was added to the cell-free culture broth, obtained by centrifugation at 10,000 rpm for 5 min. Samples were incubated for 24 h, and the emulsion index (E24) was calculated according to the following equation [21]: (E24) = 100 × (height of the emulsion layer/total height) (1)

Surface Tension
Surface tension was measured using the Du Nouy platinum ring method with a Krüss K6 tensiometer (A. Krüss: Hamburg, Germany). The fungal supernatant solution (50 mL) was tested at 25 • C to evaluate the surface tension of bio and chemical surfactants [22]. The value of the surface tension was expressed as mN/m.

Fungal Adhesion to Hydrocarbons (FATH)
The conidial suspension of isolates (6.8 mL) was mixed with 3.2 mL of benzene (Merck) in acid-washed test tubes. The mixture was vortexed for 1 min and following the method of a previous study [23]. Hydrophobicity was calculated as follows: FATH = initial absorbance − absorbance after phase separation/initial absorbance × 100 (2)

Biodegradation of Crude Oil in the Presence of Fungal Isolates and/or AgNPs
This experiment was conducted to assess the impact of two available fungal isolates and/or different concentrations of AgNPs on crude oil biodegradation. The experiment included a fungal isolate (1 × 10 6 spores/mL) and an aliquot of 2 mL of inoculum, which was inoculated into MSM (100 mL) in a 250 mL Erlenmeyer flask. The cultures were incubated on a temperature-controlled shaker incubator at 150 rpm and 30 • C for seven days using 1 g crude petroleum oil as the sole carbon source. The second experiment consisted of the same previously stated steps, but with only three different concentrations of AgNPs (0.05, 0.1, and 0.2 g). All experiments were performed in triplicates. A control sample was used without the inoculum. Next, flasks were incubated at 30 • C and 150 rpm, pH 7.5, for seven days [24]. The percentage of biodegradation, FATH, surface tension, and emulsification power were determined at days three and seven, and the crude oil samples were extracted from all experiments.

Extraction and Analysis of Crude Oil Degradation by GC-MS
Extraction of crude oil carried out according to pervious method, and alterations in its chemical composition were studied by chromatographic analysis, such as gas chromatography-mass spectrometry (GC-MS). The oil extracted from the biodegraded crude oil was monitored using a Perkin-Elmer GC Clarus 500 system (PerkinElmer: Boston, MA, US) flame ionization detector. Crude oil was separated using a capillary column (30 m × 0.25 mm × 0.25 µm) under the conditions described in a previous study [25]. The inlet temperature program was 50 • C/min. The initial temperature of the oven was maintained at 60 • C for 2 min and increased linearly at a rate of 6 • C/min and was maintained at a final temperature of 300 • C. The operating temperatures of the injector and detector were 300 • C and 320 • C, respectively. The detector temperature was 300 • C. Then, 1 µL of sample was injected with a 1:50 split ratio, and the total run time was less than 35 min.

Identifying Significant Variables Using Box-Behnken Design (BBD)
The biodegradation process was carried out by conducting experiments to determine which variables significantly affected biodegradation efficiency. The BBD for RSM was used to study the combined effect of temperature, pH, and crude oil concentration on the biodegradation of crude oil over three levels. The independent parameters and coding levels of crude oil degradation are shown in Table 1. The independent variables were studied at three different levels: low (−1), medium (0), and high (+1). Crude oil-based MSM was incubated in a shaker at a rate of 160 r/min for seven days and contained as described previously [26]. The absorbance was measured at 225 nm with petroleum ether as a reference. Each sample was measured three times, and the average value was taken.

Genomic DNA Extraction
The fungal isolate spore suspensions (1 × 10 6 spore/mL) were inoculated into double-layer media in 50 mm Petri dishes, one with liquid peptone yeast glucose (1200 µL), and the other with solid potato dextrose agar as a film. DNA extraction was performed according to a previously described protocol [27].

Data Analysis
Statistical analysis was done using SPSS (version 15,IBM SPSS Statistics, International Business Machines Corporation, Armonk, NY, USA) was used to assess the ability to utilize petroleum hydrocarbons among the 4 fungal isolates based on their abilities. A randomized complete block design was used in the present study. The standard deviation is calculated by SPSS using the descriptive procedure.

Molecular Identification
The morphological and microscopic fungal isolate characteristics and molecular tools used for isolate identification were compared with those of reference isolates in the NCBI GenBank. These four isolates were classified as A. terreus KC462061, A. flavus AF15, T. harzianum TH07, and F. solani FS12, as shown in Table 2. In the BLAST-based analysis, the ITS region sequence of each identified fungal isolate showed that isolate A. terreus KC462061 was 99% similar to A. terreus GU 966497 ( Figure S1). Using the same method, A. flavus AF15 showed 98% similarity with A. flavus KY488467 ( Figure S2), and T. harzianum TH07 showed 99% similarity with T. harzianum KC569346 ( Figure S3), while F. solani FS12 showed 98% similarity with F. solani MF136402 ( Figure S4). The sequence of the ITS1-5.8S-ITS2 region appeared to be a credible molecular approach for fungal identification and discrimination [30]. The ITS region is able to differentiate and provide accurate and rapid identification of many fungi genera, such as Aspergillus, Fusarium, and Trichoderma at the species level [31].

Characterization of Fungal AgNPs
AgNPs were examined using TEM, SEM, and EDS. Figure 1a illustrates the TEM micrograph of the biosynthesized AgNPs. The TEM micrograph shows that the as-synthesized NPs had a spherical morphology, in which nanoparticles were not conglomerated. The particle size ranged from 35 to 60 nm. TEM imaging is an essential method for observing the formation and stabilization of AgNPs.
The size and shape of the biosynthesized nanoparticles appear to depend on the type of microorganisms and other factors, such as temperature and pH of the medium, thus, a considerable size variability in AgNPs is produced by different Aspergillus spp. [32].
The SEM image of the AgNPs is shown in Figure 1b. The SEM micrograph of the AgNPs did not show a uniform surface but had three dimensions. The nanoparticles tended to be oval to spherical in shape. Most of the nanoparticles were aggregated, and no individual particles were present. Several AgNPs clusters were present, which could be attributed to the aggregation of NPs formed throughout specimen preparation [33].
Furthermore, the elemental structure of the synthesized specimen was also measured by EDX analysis (Figure 1c), showing the formation of pure AgNPs. However, the mass present was approximately 85% of the sample. The EDS profile also showed a strong silver signal along with weak carbon peaks, strongly suggesting that Ag was the major element. Clear identification of the synthesized AgNPs of the elemental composition profile was achieved. The intense signal at 3 keV strongly suggested that Ag was the main element showing optical absorption in this range due to the resonance of the surface plasmon [34]. The size and shape of the biosynthesized nanoparticles appear to depend on the type of microorganisms and other factors, such as temperature and pH of the medium, thus, a considerable size variability in AgNPs is produced by different Aspergillus spp. [32].
The SEM image of the AgNPs is shown in Figure 1b. The SEM micrograph of the AgNPs did not show a uniform surface but had three dimensions. The nanoparticles tended to be oval to spherical in shape. Most of the nanoparticles were aggregated, and no individual particles were present. Several AgNPs clusters were present, which could be attributed to the aggregation of NPs formed throughout specimen preparation [33].
Furthermore, the elemental structure of the synthesized specimen was also measured by EDX analysis (Figure 1c), showing the formation of pure AgNPs. However, the mass present was approximately 85% of the sample. The EDS profile also showed a strong silver signal along with weak carbon peaks, strongly suggesting that Ag was the major element. Clear identification of the synthesized AgNPs of the elemental composition profile was achieved. The intense signal at 3 keV strongly suggested that Ag was the main element showing optical absorption in this range due to the resonance of the surface plasmon [34].

Crude Oil Degradation by Fungal Isolates
During the fungal degradation of crude oil, the investigated fungal isolates showed perfect levels of crude oil biodegradation ability (Table 3). A. flavus AF15 showed the maximum biodegradation ability of 41.94%, followed by T. harzianum TH07 (37.36%), while third-ranking was Aspergillus terreus KC462061 with moderate biodegradation ability whereas the minimum biodegradation ability (19.42%) was recorded for F. solani FS12. In addition to the surface tension, FATH, and emulsification activity, assays were also performed for the four fungal isolates (Table 3). A. flavus AF15 had the best physicochemical properties with high emulsification activity (E24; 51.57%) and low FATH (36.74%) and surface tension (39.80 N/m), followed by T. harzianum TH07.

Crude Oil Degradation by Fungal Isolates
During the fungal degradation of crude oil, the investigated fungal isolates showed perfect levels of crude oil biodegradation ability (Table 3). A. flavus AF15 showed the maximum biodegradation ability of 41.94%, followed by T. harzianum TH07 (37.36%), while third-ranking was Aspergillus terreus KC462061 with moderate biodegradation ability whereas the minimum biodegradation ability (19.42%) was recorded for F. solani FS12. In addition to the surface tension, FATH, and emulsification activity, assays were also performed for the four fungal isolates (Table 3). A. flavus AF15 had the best physicochemical properties with high emulsification activity (E24; 51.57%) and low FATH (36.74%) and surface tension (39.80 N/m), followed by T. harzianum TH07.
Aspergillus oryzae KR029081 decreased the surface tension to 23 N/m, had higher emulsification activity (E24; 75.6%), and showed a decrease in FATH (26.5%). This indicates a very close correlation between some physical properties and crude oil degradation by fungi [35]. Penicillium spp. KY883662 reduced the surface tension and exhibited a cell surface hydrophobicity of more than 70%. Furthermore, it showed effective crude oil-degrading activity and crude oil emulsification [36].

Biodegradation of Crude Oil by Two Fungal Isolates and/or AgNPs
In the present study, two fungal isolates (A. flavus AF15 and T. harzianum TH07) were grown together or separately with/without three different concentrations of AgNPs, and three different concentrations of AgNPs were used separately. The biodegradation percentage of crude oil is listed in Table 4. Some factors related to biodegradation were assayed after three and seven days. The results showed that all treatments utilized for the degradation of crude oil were potential degraders, with a biodegradation percentage ranging from 12.4% to 57.8% after seven days. The highest extent of biodegradation was 57.8% with a combination of A. flavus AF15 and T. harzianum TH07, and 0.05 g AgNPs. The degradation level was very slight when three different concentrations of AgNPs were used separately, ranging from 12.4% to 15.1% after seven days.
The fungal biodegradation of aromatic hydrocarbons mechanisms is generally related to detoxification and bioremediation by many enzymes. Three very different-yet not mutually exclusive-enzymatic mechanisms have been proposed: (1) oxidation of the aromatic ring by intracellular Cytochromes P450 monooxygenases (CYP monooxygenases) for detoxification, (2) the nonspecific and coincidental oxidative action of extracellular aromatics by excreted lignin-degrading peroxidases, (3) the same kind of coincidental action involving extracellular laccases [37]. Our results were similar to [38] who isolated crude oil-degrading fungal strain, Aspergillus sp. RFC-1 (Rumaila field oil in Basra, Iraq); they reduced surface tension and also found an inverse relationship between potent the crude oil biodegradation of fungi and a decrease in surface tension. Correlation of Fungal Biodegradation and crude oil depends on fungi increasing emulsification activity and reduced surface tension, increasing the fungal growth rate on hydrocarbons. Fungi also play a critical role in modifying cell surface hydrophobicity. Fungal isolates can act as emulsifying agents by decreasing surface tension, improving crude oil biodegradation through two mechanisms. The first mechanism promotes the emulsification of hydrophobic compounds to mycelial pellet structures. The second mechanism induces high cell surface hydrophobicity, thus increasing the direct physical contact between cells and poor water-soluble substrates [38,39].
Furthermore, fungal consortia play a key role as highly efficient agents for the bioremediation of contaminated environments [40]. A consortium of an enzyme laccase, manganese peroxidase (MnP), Lignin peroxidase (LiP), and catalase showed a significant level of enzyme activity compared to a solo enzyme. A fungal consortium provides a greater spectrum of enzymes compared to an individual fungus, and remarkably higher enzyme activity due to the synergistic effect of the pooled genotypes. Moreover, it establishes very effective bioremediation agents [6]. Most Basidiomycete fungi are capable of producing extracellular ligninolytic enzymes (laccase, manganese peroxidase, lignin peroxidase, and versatile peroxidase) and accessory enzymes (H 2 O 2 -generating enzymes), both of which are responsible for the degradation of lignin. Because of the nonspecific nature of an extracellular fungal ligninolytic system, this could be a useful and powerful tool for bioremediation purposes [41]. The degradation percentage of crude oil by Bacillus licheniformis was up to 60% when the medium contained specific concentrations of Fe 2 O 3 and Zn 5 OH 8 Cl 2 nanoparticles and biosurfactant [42]. Genetic modification technologies (GM) and genome-wide (-omics) technologies are available that can act as enabling technologies, so-called 'ecogenomics', to develop bioremediation in the field. Ecogenomic approaches could be used to characterize contaminated sites and monitor the bioremediation process [37].

Analysis of Crude Oil Biodegradation by GC-MS
GC-MS of the crude oil treated by the fungal consortium, AgNPs, and control samples of crude oil are shown in Figure 2. From this figure, it can be observed that the fungal consortium and/or AgNPs behaved differently in the degradation of crude oil. Figure 2a shows the crude oil control with the main hydrocarbon compound peaks of C 14 , C 15 , C 16 , C 17 , C 24 , and C 27 , which were identified at retention times of 6.250, 6.945, 7.700, 8.350, 9.400, and 10.450 min by comparing the data with standard library compounds. All the hydrocarbons present in the crude oil identified using GC-MS analysis were effectively degraded by different treatments containing fungal consortia and AgNPs. Figure 2b displays a good biodegradation percentage (53.4%) when A. flavus AF15 combined with 0.05 g AgNPs. Figure 2c appears to show a weak biodegradation percentage when high concentrations of AgNPs (0.2 g) were used to separately express the negative effect of AgNPs. Figure 2d shows the highest biodegradation percentage (57.8%) when a fungal consortium consisting of A. flavus AF15 and T. Harzianum TH07 and AgNPs 0.05 g was used, and the positive effect of the AgNPs was apparent.

Analysis of Crude Oil Biodegradation by GC-MS
GC-MS of the crude oil treated by the fungal consortium, AgNPs, and control samples of crude oil are shown in Figure 2. From this figure, it can be observed that the fungal consortium and/or AgNPs behaved differently in the degradation of crude oil. Figure 2a shows the crude oil control with the main hydrocarbon compound peaks of C14, C15, C16, C17, C24, and C27, which were identified at retention times of 6.250, 6.945, 7.700, 8.350, 9.400, and 10.450 min by comparing the data with standard library compounds. All the hydrocarbons present in the crude oil identified using GC-MS analysis were effectively degraded by different treatments containing fungal consortia and AgNPs. Figure 2b displays a good biodegradation percentage (53.4%) when A. flavus AF15 combined with 0.05 g AgNPs. Figure 2c appears to show a weak biodegradation percentage when high concentrations of AgNPs (0.2 g) were used to separately express the negative effect of AgNPs. Figure 2d shows the highest biodegradation percentage (57.8%) when a fungal consortium consisting of A. flavusAF15 and T. Harzianum TH07 and AgNPs 0.05 g was used, and the positive effect of the AgNPs was apparent. GC-MS is a highly accurate, effective, and versatile analytical technique with numerous scientific applications. GC-MS has become a highly recommended tool for monitoring and tracking hydrocarbons and chemical pollutants in the environment. It can be used to screen the degradation products of hydrocarbons, even C 60 [43]. Using GC-MS data to analyze the biodegradation of crude oil by microbial consortium, it was found to be strongly degraded with a degradation ratio of up to 56% after four weeks of treatment [6]. GC-MS analysis showed that the A. oryzae KR029081 isolate was very unique and degraded with a hydrocarbon from C 14 to C 27 with biodegradation of 99% [35].

Optimization of Biodegradation Using RSM
Five significant variables were optimized using the BBD (Table 5). Table 5 shows good agreement between the predicted and experimental values. A total of 30 experiments were performed with a different combination of five variables. The maximum biodegradation of crude oil (57.21%) was achieved in run number 16 under the following conditions: temperature of 30 • C, pH 7, crude oil concentration of 4 g/L, 1:1 ratio between A. flavus AF15 and T. harzianum TH07, and AgNP concentration of 0.05 g. Analysis of regression coefficients and t-value of five variables are presented in Table 6. The results showed that temperature (confidence level, 99.81%), pH (99.66%), crude oil concentration (99.52%), ratio between A. flavus AF15 and T. harzianum TH07 (96.18%), and AgNP concentration (95.79%) were the most significant factors influencing biodegradation efficiency. Table 5. Full-factorial Box-Behnken design matrix of 5 independent variables and mean response for biodegradation percentage of crude oil.

Run
Temperature pH Crude Oil Conc.

Ratio between AF and TH Nano
Conc.

Biodegradation (%)
Actual Predicted  The three most important factors were temperature, pH, and AgNPs concentration. These factors were used to establish a mathematical model using RSM technology. The mathematical model has optimized conditions of biodegradation utilizing these factors, and the best degradation rate was 81.63% with a temperature of 30 • C, pH of 7.14, and TPH concentration of 4.83 g/L [26].   Table 6. Estimated regression coefficients for optimization of biodegradation efficiency using Box-Behnken design. The three most important factors were temperature, pH, and AgNPs concentration. These factors were used to establish a mathematical model using RSM technology. The mathematical model has optimized conditions of biodegradation utilizing these factors, and the best degradation rate was 81.63% with a temperature of 30 °C, pH of 7.14, and TPH concentration of 4.83 g/L [26].

T. harzianum Isolates
Of the 15 primers used, nine produced unambiguous fragments with repeatable patterns. 62 DNA bands were obtained with an average of 6.8 DNA bands per primer, and 40 DNA bands (64.5%) were polymorphic bands with an average of 4.4 polymorphic DNA bands per primer. The genetic similarity coefficient ranged from 0.82 to 0.95 Figure 5. The highest number of bands (10) was obtained with two primers, (TG) 8 G and (AG) 8 YT, while the lowest number five was obtained with primer (GA) 8 T (Figure 6a). Only one primer, (GA) 8 T, had good polymorphism (100%), and other polymorphic primers exhibited low values of polymorphism, ranging from 44.5% to 80% (Figure 6a). The genetic relationships among the nine T. harzianum isolates were also visualized by 2D PCA for nine ISSR primers (Figure 6b).
ISSR-PCR is a highly effective method for the characterization of Aspergillus and Trichoderma spp., and has also been shown to be useful for assessing the genetic diversity of two species. Moreover, it has been proven to be useful for DNA fingerprinting using ISSR genetic markers [28,44].

T. harzianum Isolates
Of the 15 primers used, nine produced unambiguous fragments with repeatable patterns. 62 DNA bands were obtained with an average of 6.8 DNA bands per primer, and 40 DNA bands (64.5%) were polymorphic bands with an average of 4.4 polymorphic DNA bands per primer. The genetic similarity coefficient ranged from 0.82 to 0.95 Figure 5. The highest number of bands (10) was obtained with two primers, (TG)8G and (AG)8YT, while the lowest number five was obtained with primer (GA)8T (Figure 6a). Only one primer, (GA)8T, had good polymorphism (100%), and other polymorphic primers exhibited low values of polymorphism, ranging from 44.5% to 80% (Figure 6a). The genetic relationships among the nine T. harzianum isolates were also visualized by 2D PCA for nine ISSR primers (Figure 6b).  ISSR-PCR is a highly effective method for the characterization of Aspergillus and Trichoderma spp., and has also been shown to be useful for assessing the genetic diversity of two species. Moreover, it has been proven to be useful for DNA fingerprinting using ISSR genetic markers [28,44].

Conclusions
The two fungal isolates used in this study demonstrated rapid crude oil biodegradation ability and when used together as a consortium displayed a cooperative effect that enhanced the biodegradation process. In addition, we have shown the positive effect of low concentrations of AgNPs on biodegradation. RSM is an efficient mathematical method to optimize the microbial biodegradation of crude oil. ISSR-PCR is a highly effective method for the characterization of fungi responsible for biodegrading crude oil.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Figure S1: Phylogenetic tree of A. terreus biocoAT (KC462061) based on ITS region and 5.8S sequences, Figure S2: Phylogenetic tree of A. flavus biocoAF15 based on ITS region and 5.8S sequences, Figure S3: Phylogenetic tree of T. harzianum biocoTH07 based on ITS region and 5.8S sequences, Figure S4: Phylogenetic tree of F. solani biocoFS12 based on ITS region and 5.8S sequences.

Conclusions
The two fungal isolates used in this study demonstrated rapid crude oil biodegradation ability and when used together as a consortium displayed a cooperative effect that enhanced the biodegradation process. In addition, we have shown the positive effect of low concentrations of AgNPs on biodegradation. RSM is an efficient mathematical method to optimize the microbial biodegradation of crude oil. ISSR-PCR is a highly effective method for the characterization of fungi responsible for biodegrading crude oil.
Supplementary Materials: The following are available online at http://www.mdpi.com/1660-4601/17/18/6564/s1, Figure S1: Phylogenetic tree of A. terreus biocoAT (KC462061) based on ITS region and 5.8S sequences, Figure S2: Phylogenetic tree of A. flavus biocoAF15 based on ITS region and 5.8S sequences, Figure S3: Phylogenetic tree of T. harzianum biocoTH07 based on ITS region and 5.8S sequences, Figure S4: Phylogenetic tree of F. solani biocoFS12 based on ITS region and 5.8S sequences.