Next Article in Journal
Phenolic Composition and Antioxidant Activity of Tanacetum parthenium Cultivated in Different Regions of Ukraine: Insights into the Flavonoids and Hydroxycinnamic Acids Profile
Next Article in Special Issue
Effects of Common Biochar and Acid-Modified Biochar on Growth and Quality of Spinach in Coastal Saline Soils
Previous Article in Journal
Plant Resistance Inducers Affect Multiple Epidemiological Components of Plasmopara viticola on Grapevine Leaves
Previous Article in Special Issue
Response of Soil Microenvironment and Crop Growth to Cyclic Irrigation Using Reclaimed Water and Brackish Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combination of Biochar and Trichoderma harzianum Can Improve the Phytoremediation Efficiency of Brassica juncea and the Rhizosphere Micro-Ecology in Cadmium and Arsenic Contaminated Soil

1
State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi’an University of Technology, Xi’an 710048, China
2
PowerChina Northwest Engineering Corporation Limited, Xi’an 710065, China
*
Author to whom correspondence should be addressed.
Submission received: 13 July 2023 / Revised: 1 August 2023 / Accepted: 10 August 2023 / Published: 14 August 2023

Abstract

:
Phytoremediation is an environment-friendly method for toxic elements remediation. The aim of this study was to improve the phytoremediation efficiency of Brassica juncea and the rhizosphere soil micro-ecology in cadmium (Cd) and arsenic (As) contaminated soil. A field experiment was conducted with six treatments, including a control treatment (CK), two treatments with two contents of Trichoderma harzianum (T1: 4.5 g m−2; T2: 9 g m−2), one biochar treatment (B: 750 g m−2), and two combined treatments of T1B and T2B. The results showed Trichoderma harzianum promoted the total chlorophyll and translocation factor of Brassica juncea, while biochar promoted plant biomass compared to CK. T2B treatment showed the best results, which significantly increased Cd accumulation by 187.49–308.92%, and As accumulation by 125.74–221.43%. As a result, the soil’s total Cd content was reduced by 19.04% to 49.64% and total As contents by 38.76% to 53.77%. The combined amendment increased the contents of soil available potassium, phosphorus, nitrogen, and organic matter. Meanwhile, both the activity of glutathione and peroxidase enzymes in plants, together with urease and sucrase enzymes in soil, were increased. Firmicutes (dominant bacterial phylum) and Ascomycota (dominant fungal phylum) showed positive and close correlation with soil nutrients and plant potentially toxic elements contents. This study demonstrated that phytoremediation assisted by biochar and Trichoderma harzianum is an effective method of soil remediation and provides a new strategy for enhancing plant remediation efficiency.

1. Introduction

Industrial activities and agricultural production lead to the accumulation of potentially toxic elements (PTEs) in the soil environment, which is a big danger to the environment and human health [1,2]. Cadmium (Cd) and arsenic (As) can be effectively integrated into food chains through soil–plant systems. These two elements are very poisonous. They can damage the skin and raise the chance of getting cancer, according to IARC (International Agency for Research on Cancer) [3,4]. Therefore, effective treatment of Cd- or As-contaminated soils is urgently needed.
Phytoremediation can reduce the level or toxicity of PTEs in the soil environment. Compared to physical (plowing), chemical (soil cleansing), and electrical (electro-extraction) soil remediation methods [5], phytoremediation is a beneficial remediation technique because it prevents secondary pollution and can be used for massive-scale applications [6,7]. Plants with high biomass and good tolerance to PTEs have been discovered for phytoremediation in recent years [8]. Brassica juncea is an energy plant belonging to the Brassica genus in the cruciferous family that is highly adaptable and can grow in a variety of harsh environments [9,10]. Brassica juncea is a perfect plant source for phytoremediation that has been proven because it has high productivity, a well-established root system, and good resistance to a range of PTEs [11,12,13].
Nonetheless, phytoremediation has a drawback of low effectiveness caused by insufficient mass and PTEs accumulation ability. Therefore, investigators in and abroad have committed to combining amendments on phytoremediation in rehabilitating soil environments contaminated with PTEs [14,15]. Biochar was completely verified to be capable of lowering the biotoxicity of PTEs due to its high porosity and ion exchange capacity [16,17]. By boosting microbial activity in PTEs-polluted soils, biochar can improve the extraction of PTEs from Brassica napus, as previous research demonstrated [18]. Biochar can also directly contribute to the availability of phosphorus [19] and alleviate plant potassium deficiency by increasing exchangeable potassium in soil [20], as well as improving soil quality [21,22].
Furthermore, plant growth-promoting fungi (PGPF) have been found to greatly benefit plant growth and their ability to resist invasive plant PTEs by their hosts [23,24]. Numerous studies have explored how PGPF, such as Aspergillus, Penicillium, and Trichoderma, associated with their host, are extremely beneficial to plant growth and stress resistance to invasive PTEs [23,25,26]. Trichoderma harzianum (TH), known as one of PGPF, has been shown to have high resistance to Cd and Pb in the experimental area with high PTEs contamination [27], with which endomycorrhizae are formed during mutualistic symbiosis with plant roots [28]. Then, endomycorrhizae can provide a variety of services to the host plant, including nutrients providing for plant growth, as well as enhancing abiotic and biotic stress tolerance [29]. Coincidentally, biochar can create appropriate habitats for beneficial PGPF, promote rapid colonization of crop rhizosphere by fungal inoculants, and mitigate biotic and abiotic plant stress [30]. An efficient and novel technique for the sustainable recovery of polluted soils, recognized by many researchers, is the fusion of biochar with functional fungal strains to improve phytoremediation [31,32].
Based on the studies above, we hypothesize that the effectiveness of phytoremediation for soils polluted with PTEs can be improved by applying biochar and TH together. However, there is a lack of research on how soil properties and functions are restored during phytoremediation with biochar and PGPF on soils contaminated with PTEs, and how these properties change dynamically. Therefore, a field trial was conducted by us to evaluate the hypothesis that cooperative interactions occurred between biochar and TH, which promote the phytoremediation of PTEs by Brassica juncea (i.e., plant stabilization and extraction). The primary goals of this research are as follows: (1) to examine the joint impacts of biochar and TH on the development of Brassica juncea, soil condition, and rhizosphere microbial communities; (2) to uncover the mechanisms behind the combined impact of biochar and TH on the efficiency of phytoremediation, from the viewpoint of plant–microorganism interactions; (3) to determine the coupling mechanisms between soil environmental factors and soil microbial factors; (4) to determine the crucial factors and mechanisms that influence the enhancement of phytoremediation in mustard plants by biochar and TH.

2. Results

2.1. Plant Growth and PTEs Content

After 26 weeks’ growth of Brassica juncea, the Cd contents both in the shoot and root were decreased with the application of biochar compared to that in the control treatment, while TH enhanced the Cd uptake both in the root and shoot, which is beneficial for Cd transferring from soil to aboveground parts (Figure 1). The highest Cd uptake by Brassica juncea was found in the T2B treatment, of which the Cd contents in the root and shoot were increased by 35.30% and 41.17% compared to that in the control (Figure 1A). It could also be found that all treatments showed a greater As uptake compared to that in control (Figure 1B). Combined biochar and TH increased PTEs accumulation and translocation factor in Brassica juncea (Figure 1C). The combined treatment produced the greatest biomass of Brassica juncea, in which the biomass of shoot and root were increased by 125.85% to 242.18% and 103.13% to 185.94%, respectively. The total amount of chlorophyll increased significantly under the treatment of biochar and TH (Table S3).

2.2. Soil PTEs Content and Availability

Figure 2A,B depict the analysis of samples for T-Cd and T-As from different soil layers. They were taken to assess the impact of six treatments on the PTEs content in the soil. In comparison to the control treatment, the elimination of PTEs in the root zone soil was improved by biochar and TH; furthermore, the combined treatment of biochar and TH showed a better remediation effect. For T2B treatment, the contents of Cd and As decreased by 28.54% and 38.76% in the rhizosphere soil, 27.73% and 40.74% in the 0–15 cm deep soil layer, 49.64% and 53.77% in the 15–30 cm deep soil layer, and 19.04% and 51.30% in the 30–45 cm deep soil layer. For T2B treatment, the available Cd content in the rhizosphere soil was significantly increased by 21.06% compared to the single application of biochar (Figure 2C). The morphological distribution of As in rhizosphere soil and a surface layer of 15 cm soil showed similar results compared to that in control, in which the Ca-As ratio increased significantly and As activity decreased (Figure 2D,E).

2.3. Enzyme Activity and Soil Physicochemical Properties

Different soil physicochemical properties after various remediation treatments were observed in Figure 3 and Table S4. All soil layers had a weakly alkaline pH, with values ranging from 8.04 to 8.89, as shown in Table S4. Both TH and biochar promoted the release of AN, AP, and AK in the root zone soil and the 15 cm soil layer. Meanwhile, biochar positively promoted OM content in all soil layers. The soil enzyme activities, such as urease, ALP, and sucrase in the rhizosphere, varied greatly. Compared to the control, TH significantly increased the enzyme activities in the rhizosphere soil (Figure 3A–C). Combined treatment increased the activities of urease and sucrase by 87.25% to 97.79% and 101.60% to 111.20%, respectively (Figure 3A,B). Biochar had a promoting effect on the activities of urease and sucrase but reduced ALP activity by 13.75% to 25.62% (Figure 3C).
During the phytoremediation of PTE-contaminated soils, changes in plant enzyme activity were also observed (Figure 3D,E). Both biochar and TH enhanced the POD activity and GSH content in plant shoots and roots compared to that in the control (Figure 3D,E). Concurrently, the root GSH content reached its maximum in the T2B treatment (Figure 3E). Both biochar and TH reduced the MDA content in plant shoots and roots, and the lowest MDA content in shoots and roots was observed in the T2B treatment, with reductions of 75.12% and 7.69%, respectively (Figure 3F).

2.4. Diversity and Community Structure of Rhizosphere Soil Microbial Community

Through high-throughput sequencing, all soil samples gave a total number of 341,197 effective bacterial gene sequences and 435,504 effective fungal gene sequences. The bacterial sequences were grouped into 12,816 Operational Taxonomic Units (OTUs), while the fungal sequences were grouped into 4628 OTUs. The sequencing data reflect the bulk of microbial diversity information in the remediation treatment sufficiently as the rarefaction curve nears saturation (Figure S1). According to the Shannon and Simpson indices, the alpha diversity indices, as well as the bacterial and fungal populations, increased in the B and T1 treatments. However, As the application rate of TH increased, soil bacterial and fungal diversity decreased, with bacterial diversity being greater than fungal diversity, as shown in Table S5.
The taxonomic classification of the entire qualified OTUs has been determined to include 80 different bacterial phyla and 16 fungal phyla. Specific physicochemical parameters and PTEs concentrations from different remediation treatments can affect microbial populations. In all treatments, Proteobacteria, Acidobacteriota, and Firmicutes were the prevalent bacterial phyla, constituting 76.65–87.10% of the total phylum level. The relative proportion of Proteobacteria varied between 21.29% and 57.06% (Figure 4A). The ternary plot analysis showed that Proteobacteria contributed the largest proportion to the ratio in the combined amendment, which exceeded 60% (Figure 4C). For the fungal community, Ascomycota is the dominant fungal phylum, with a relative abundance of 33.85–80.24% (Figure 4B). The proportion of Ascomycota in the combined treatment exceeded 50% (Figure 4D). The Proteobacteria phylum had 11 genera classified as belonging to it. The Firmicutes phylum had eight genera classified as belonging to it. The clustering analysis showed that there was good consistency between the B and T2B treatments in the bacterial community. The dominant bacterial genera in the soil under the T2B treatment are Bacillus, Lechevalieria, and Pseudomonas (Figure S2A). The Ascomycota phylum had 22 genera classified as belonging to it. In the treatment of TH and combined treatment, the predominant fungal genera were Trichoderma, Penicillium, and Gibberella (Figure S2B).

2.5. Responses of Phytoremediation to Microhabitat Properties of Rhizosphere Soil

In order to assess the impact of soil characteristics and PTEs levels on microbial populations in polluted farmland, a redundancy analysis was conducted (Figure 5). According to RDA, the structure of bacterial and fungal communities was mainly influenced by factors such as T-Cd, T-As, CEC, OM, AP, and AK (p-value < 0.05). The structure of fungal communities was significantly impacted by A-Cd, with a p-value of less than 0.001 (Figure 5A,B). Procrustes analysis was further performed to investigate the correlations between different soil environmental factors and soil microbes (Figure 5C,D). The results demonstrated that the overall variation in environmental factors was considerably uniform among bacterial communities.
In order to investigate the associations between soil properties and microbial communities, a Spearman correlation heatmap combined with UPGMA cluster analysis was utilized (Figure 6). Soil pH did not have a significant effect on the abundance of the bacterial and fungal phyla. However, the soil properties of CEC, OM, AN, AP, and AK had a significant positive association with Firmicutes (a dominant bacterial phylum) and Ascomycota (a dominant fungal phylum). The soil Cd and As exhibited a substantial inverse correlation with Acidobacteriota and Firmicutes (dominant bacterial phyla), as well as Ascomycota (a dominant fungal phylum).
The association between plant characteristics and rhizosphere soil ecological factors was measured by Pearson correlation analysis (Figure S3). The various plant and soil parameters were divided into two groups: Group A (pH, CEC, AN, AP, AK, OM, Shoot Cd, Root Cd, Shoot As, Root As, Shoot POD, Root POD, Shoot GSH, Root GSH, Urease, and Sucrase) and Group B (Shoot MDA, Root MDA, Cd, As, and ALP). Based on the Pearson correlation analysis, the factors of Group A and Group B demonstrated a pronounced negative link with each other, while the environmental factors within each group revealed a considerable positive connection with each other.
The levels of Cd and As in plants were inversely related to PTEs levels in various soil layers. In contrast, the levels of Cd and As in rhizosphere soil had a direct relationship with A-Cd and A-As (Figure 7A). Soil physical and chemical characteristics positively influenced the functions of soil urease and sucrase, whereas soil PTEs adversely affected them, with A-Cd and A-As being the most important factors (Figure 7B). The accumulation of Cd and As in Brassica juncea was found to be inversely related to soil PTEs levels, with T-Cd and T-As being the most crucial factors (Figure 7C). The correlation of MDA activity in plants with soil environment was opposite to that of POD and GSH, with T-As, A-As, and A-Cd being the most critical factors (Figure 7D).

3. Discussion

3.1. Soil PTEs Content and Plant PTEs Accumulation

Phytoremediation was an effective method in promoting the transport of PTEs in soil layers at different depths, transforming availability to organisms of PTEs in the rhizosphere soil [33]. PTEs transformation can be facilitated by biochar into less available forms in soil without reducing the whole quantity of PTEs [34], and the poisonousness of Cd and As to soil microbes and plants was mitigated by biochar [35]. The combination of biochar and TH significantly reduced T-Cd and T-As content in soils at different depths during phytoremediation, with the greatest reduction observed in soil layers at 30 cm and 40 cm depths (Figure 2). Additionally, greatly improved performance of plant-assisted cleanup for PTEs-contaminated soil was achieved through the synthesis of biochar and TH treatment. The formation of overall chlorophyll and the accumulation of biomass in the plants was enhanced by TH (Table S3) [36,37]. The symbiotic relationship of endophytic fungi had various benefits for plants, such as the elevated production of phytohormones and the maintenance of the internal hormone equilibrium [38]. TH enhanced antioxidant enzyme activities and improved glutathione content in mustard. The plant’s capacity to transport and accumulate PTEs can be strengthened by altering its cellular mechanisms and biological processes [29,39,40]. The concentration of PTEs in plants demonstrated a significant negative link with the PTEs content in the soil. However, a remarkable positive association was discovered between the PTEs content in the inter-root area and the PTEs content in each soil layer (Figure 7A). The above findings mean that the combined amendment method can encourage the movement of cadmium from deeper soil layers to the surface. Additionally, it can significantly increase the PTEs accumulation of Brassica juncea, and it promotes the migration of cadmium and arsenic from plant roots to shoots.

3.2. Plant Stress Tolerance and Soil Physicochemical Properties

The implementation of biochar had a beneficial impact on AN, AP, and AK in rhizosphere soil and the top 15 cm soil layer, as well as on OM content in each soil layer (Table S4). It is typically believed that soil pH has a substantial impact on how accessible soil PTEs are [41]. However, since the pH was not improved in this study, the availability of soil PTEs was affected by other factors. OM was shown to be adversely associated with soil Cd and As (Figure S3), which means the increased soil OM content inhibited soil PTEs activity and reduced PTEs toxicity [42,43]. In this study, the combined treatment was found to heighten the performance of urease and sucrase enzymes while simultaneously curbing ALP activity. Based on the previous study, ALP activity has a significant relationship with soil AP. In this study, AP contents were improved greatly, especially in the treatment with biochar. ALP gene expression is induced in phosphorus-poor conditions, and ALP expression is inhibited by increases in AP content [44]. The interaction of PTEs ions with enzyme thiol groups may also play a part in the suppression of soil enzyme activity caused by PTEs [41]. The activity of the enzyme’s urease and sucrase in the soil correlated favorably with soil pH, CEC, OM, AN, AP, and AK, where AN and AP were the most significant factors. However, soil urease and sucrase enzyme activity exhibited a negative correlation with soil PTEs, where A-Cd and A-As were the most significant factors. Enhancing the activity of urease and sucrase enzymes, biochar, and TH increased soil nutrients levels and minimized the effects of PTEs [45]. Reactive oxygen species (ROS) can be greatly increased once PTEs have entered the plant tissues, producing cellular membrane damage and antioxidant enzyme deactivation, resulting in plant tissue damage and even death [46]. Biochar and TH reduced MDA activity in plants and significantly increased POD and GSH activities (Figure 3), which mitigated Cd and As-induced oxidative damage. The increase in POD and GSH activities was more pronounced in the roots than in the shoots, and this is related to TH’s capacity to encourage Brassica juncea’s uptake of PTEs [47]. GSH plays a defensive role by chelating PTEs and detoxifying foreign pollutants, which can reduce the damage of total PTEs to plants [48]. The increase in GSH content helps plants to accumulate PTEs, thereby improving the efficiency of remediation. A favorable relationship was found between soil nutrient levels and POD activity (Figure 7D). The improvement of soil nutrient status helps to improve the growth and enzyme activity of plants [49].

3.3. Microbial Community Diversity and Structure Differences

The biological diversity and organization of soil organisms can be significantly impacted by PTEs contamination, which can have a dramatic impact on the physical and chemical characteristics of soil [50]. In the fungal community, the composition of microorganisms became more concentrated. TH inhibited the proliferation of other fungi in the soil by colonizing the plant rhizosphere [51]. But with increasing TH application rate, bacterial and fungal diversity within soil declined. The main cause is that TH made soil A-Cd and A-As concentrations higher; as a result, some of the microorganisms that could not tolerate PTEs died. The soil microbial quantity and diversity then decreased [52,53]. Proteobacteria, as the most abundant bacterial phylum in the soil environment, has been reported to be highly tolerant to PTEs in harsh conditions [54,55]. Under the combined treatment, the dominant bacterial genera were Bacillus, Lechevalieria, and Pseudomonas. These species are very competitive and essential for PTEs tolerance and the establishment of plants [56]. The fungal phylum Ascomycota showed adaptability to harsh environments and exhibited a wide distribution in areas with serene PTEs pollution [57,58]. Ascomycota is one of the largest phyla in the fungal kingdom [59] and possesses complex spores, which can mitigate PTEs stress by immobilizing or accumulating PTEs within their bodies [60]. Under the combined amendment, the predominant fungal genera were Trichoderma, Penicillium, and Gibberella. Trichoderma and Penicillium have strong abilities in accumulating and transporting PTEs, which helps to improve the adaptation to various harsh environments [57].

3.4. Modulation Mechanisms between Phytoremediation and Rhizosphere Microecology

Phytoremediation is closely related to the ecological characteristics of the rhizosphere, including soil physicochemical properties, PTEs contents, and microbiome [61]. Soil physicochemical properties (CEC, OM, AN, AP, AK) and T-Cd and T-As concentrations were the primary elements influencing the makeup of bacterial and fungal populations. Additionally, the key drivers of the composition of the soil microorganisms were soil nutrients, including AN and AP [62]. The predominant phylum of soil bacteria, Acidobacteriota, exhibits an extremely unfavorable association with plant Cd concentration (Figure 5). Previous studies have demonstrated a similar result, showing that enrichment of Acidobacteriota can decrease Cd concentration in plants and enhance the uptake of mineral elements by plants [63]. The number of PTEs in plants was positively correlated with the bacterial phylum Firmicutes, and a competitive relationship between Acidobacteriota and Firmicutes in the dominant bacterial species was found [64,65]. In addition, contents of both Cd and As in the rhizosphere were negatively correlated with the dominant bacterial phylum Firmicutes, indicating that Firmicutes are sensitive to the contents of Cd and As [55]. Additionally, it has been demonstrated that Firmicutes have several genes for PTEs susceptibility that help them to tolerate high amounts of PTEs pollution [66]. Cd and As contents of soil revealed both adversely associated with the dominant fungal phylum, Ascomycota. However, Cd contents both in root and shoot have a significantly positive correlation with Ascomycota. This fungal phylum acts as an endophytic fungus, which can reduce PTEs cytotoxicity in roots and enhance PTEs absorption [67]. The combined treatments showed significant changes in the diversity and bacterial and fungal community structure in the soil. Dominant bacteria and fungi can compete for available resources and also cooperate with each other to form stable coexisting communities. The complex relationship mentioned above is the main reason for improving the soil rhizosphere microecology [45] and promoting the enrichment efficiency of Brassica juncea for PTEs.

4. Materials and Methods

4.1. Materials

The experiment was conducted from March 2021 to November 2022 in Xi’an, Shaanxi, China (34°18′43″ N, 108°54′31″ E). A warm temperate monsoon climate prevails in the area. It has an annual rainfall of 609.6 mm, and the average temperature during the study period is 23.3°. The soil texture is silty clay loam (18.22% clay, 45.98% silt, and 35.80% sand). Large areas of soil pollution in the experimental area were caused by prolonged irrigation of sewage. The initial physicochemical characteristics of the contaminated soil, of which Cd and As were chosen as the remediation targets, are summarized in Table S1.
Seeds of Brassica juncea were procured from Gu Yue Flower and Plant Seeds Business Department of Quanzhou County, Shandong Province. All the seeds were disinfected using a 2% sodium hypochlorite solution for 30 min, and any remaining solvents on the seed surface were then washed away with three rinses in disinfected filtered water. The microbial agent TH was obtained from Shandong Lvlong Biological Technology Co., Ltd. (Jinan, China) with an effective bacterial count of ≥1 billion/g. The biochar used in this experiment was purchased from Zhenjiang Runwu Environmental Science and Technology Co., Ltd. (Wenzhou, China). The main properties of biochar are listed in Table S2.

4.2. Experimental Design and Sample Collection

Biochar and TH were chosen as the amendment to improve Brassica juncea. The application rate of biochar (B: 750 g m−2) and TH (T1: 4.5 g m−2 and T2: 9 g m−2) were determined based on previously reported studies [68]. Six treatments were established for the field experiments: (1) Brassica juncea treatment alone (CK), (2) Brassica juncea + biochar treatment (B), (3) Brassica juncea + 4.5 g m−2 TH treatment (T1), (4) Brassica juncea + 9 g m−2 TH treatment (T2), (5) Brassica juncea + biochar + 4.5 g m−2 TH treatment (T1B), and (6) Brassica juncea + biochar + 9 g m−2 TH treatment (T2B). Six treatments were applied in randomized plots with three replicates.
Each plot has an aggregate area of 4 m2 with dimensions of 2 m in length and 2 m in breadth. To minimize interference, a 50 cm wide isolation strip was placed among the plots. Before implementing the soil pollution treatment, a 30 cm depth of deep plowing was accomplished. After three days, biochar was spread on the designated treatment plots in the plowed layer soil of 0–15 cm. Then, seeds were sown after three days, and the microbial treatment was sprayed according to the experimental setup after Brassica juncea planting.
After 26 weeks of the growth period, Brassica juncea and soil were sampled. A portion of the plant samples was used to determine plant biomass and PTEs concentrations. Plant samples were oven-heated to 105° for 30 min to inactivate biological activity before drying to consistent weight at 80° for biomass measurements. Another portion of the samples was stored in liquid nitrogen for subsequent enzyme activity assays. Based on the growth characteristics of plant root systems, the rhizosphere soil and 0–15 cm, 15–30 cm, and 30–45 cm of soil were sampled for the study. The communities of soil microbes were examined using rhizosphere soil. The remaining samples were cleaned in the lab, allowed to air dry, and then stored in Zip-lock bags. The air-dried soil samples were sieved through a 1 mm nylon sieve for organic matter, available phosphorus, available potassium, and alkaline nitrogen and then sieved through a 0.15 mm nylon sieve for PTEs.

4.3. Soil and Plant Physicochemical Properties and Enzyme Activities

A pH-automated analyzer (Mettler Toledo, Greifensee, Switzerland) was used to measure the pH of the soil at a 1:2.5 (w/v) soil-to-water ratio [69]. The volumetric potassium dichromate method with external heating was implemented to measure the amount of soil organic matter (OM) [70]. Available phosphorus (AP) was measured using a Smart Chem automatic discontinuous chemical analyzer (SmartChem450, AMS Alliance, Courtaboeuf, Italy) [71]. Atomic absorption spectrophotometry was adopted to quantify the amount of available potassium (AK) [72]. The alkali diffusion method was implemented to measure the soil’s alkali hydrolyzable nitrogen (AN) [73]. Ammonium chloride–acetic acid exchange was implemented to determine the cation exchange capacity (CEC) [74].
Urease was determined utilizing the phenol–sodium hypochlorite colorimetric technique [75], alkaline phosphatase (ALP) was measured with the sodium phosphate colorimetric method [76], and sucrase was quantified utilizing the 3,5-dinitrosalicylic acid colorimetric technique [77]. Enzyme activities both in plant shoots and roots were also determined. Sedlak and Lindsay’s approach was adopted to determine the glutathione (GSH) content [78]. Enzyme solutions were prepared in pre-cooled mortar, and plant malondialdehyde (MDA) and peroxidase (POD) contents were determined using a UV-1800 spectrophotometer (Aoyi Instruments (Shanghai) Co., Ltd., Shanghai, China) [79].

4.4. Potentially Toxic Elements Analysis

The four acids, HNO3, HCl, HF, and HClO4 (5:4:3:2, v/v), were applied to dissolve the soil samples in digestion tubes on a graphite digestion instrument (EH-2542, Beijing Auwii Instrument Co. Ltd., Beijing, China). The total Cd (T-Cd) [80], available Cd (A-Cd) [81], and total As (T-As) [82] content were determined by atomic absorption spectrophotometry. Atomic absorption spectrometry (GF-AAS) was used to measure the quantity of Cd, whereas F-AAS was used to assess the quantity of As. Meanwhile, according to the extraction procedure described in the literature, the concentrations of A-As, Fe-As, and Ca-As in soil were extracted and determined by AFS [83]. The determination of PTEs content in plants is divided into two parts, including root Cd, shoot Cd, root As, and shoot As for both aboveground and underground parts of the plant. In a graphite digestion apparatus, the plant materials are digested using a 7:1 v/v solution of HNO3 and HClO4, and the PTEs content is measured using atomic absorption spectrophotometry [82,84]. Two licensed reference materials were employed for quality control and assurance analyses: Chinese soils (GBW07401, Institute of Geophysical and Geochemical Exploration (IGGE), Beijing, China) and laver (GBW08521, National Institute of Metrology of China (NIM), Beijing, China).

4.5. Potentially Toxic Elements Translocation and Accumulation

To assess the efficacy of phytoremediation, we estimated the accumulation of PTEs per plant (in μg plant−1) and the translocation factor (TF) [85]. In these equations, PTEsshoot and PTEsroot denote the amounts of Cd or As in the shoot and root, respectively, while Bshoot and Broot indicate the shoot and root mass, respectively.
TF = PTEsshoot/PTEsroot
PTEs accumulation (μg plant−1) = PTEsshoot × Bshoot + PTEsroot × Broot

4.6. DNA Isolation, Sequencing, and Microbial Diversity Analysis

DNA from rhizosphere soil was extracted. PCR was used to amplify and sequence the 16S rRNA and 18S rRNA genes of the bacteria and fungi, respectively. Following the procedure outlined by Gao et al. [86], primers ITS5-1737F (5′-GGAAGTAAAAGTCGTAACAAGG-3′) and ITS2-2043R (5′-GCTGCGTTCTTCATCGATGC-3′) were used to amplify fungal DNA targeting the fungal ITS region. PCR amplification of the V3–V4 variable region of the bacterial 16S rRNA gene was performed using the universal primers 515F (5′-CCTAYGGGRBGCASCAG-3′) and 806R (5′-GGACTACNNGGGTATCTAAT-3′). The TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, San Diego, CA, USA) was applied for creating sequencing libraries, and index codes were added to the libraries as instructed by the manufacturer. The 18S rRNA and 16S rRNA amplicon high-throughput paired-end sequencing was performed using Illumina (Kapa Biosciences, Woburn, MA, USA) NovaSeq6000 platform. The NCBI Sequence Read Archive (SRA) database has accepted the raw sequencing data from this study and assigned accession numbers for fungi and bacteria (PRJNA936863 and PRJNA936852, respectively).

4.7. Statistical Analysis

Data statistics and analysis were performed using IBM SPSS Statistics 26. Duncan’s multiple range test was employed to calculate significance at a level of p < 0.05. SigmaPlot 14.0 was used for the quantification and plotting of soil and plant-related indicators. The R programming language (http://www.r-project.org, accessed on 1 May 2023) was used to investigate the probable connection between microbial populations and soil parameters of the environment. The correlation analysis approach of redundancy analysis (RDA) was utilized to determine the important environmental factors that influence microbial community structure.

5. Conclusions

Compared to the single amendment, the biochar and Trichoderma harzianum combined application considerably improved the phytoremediation efficiency of Brassica juncea and enhanced the rhizosphere microbial community and soil quality. The main results were as follows: Biochar significantly contributed to the accumulation of biomass in Brassica juncea, while the activity of antioxidant enzymes and PTEs transport was supported by TH. Biochar and TH synergistically enhance the biomass accumulation and PTEs transport of Brassica juncea. This results in increased accumulation of Cd and As in plants, thereby improving plant remediation efficiency. Using TH and biochar together decreased the MDA in plants and ALP in the soil while increasing the GSH, POD, soil urease, and sucrase. Additionally, the combined biochar and TH treatment increased AN, AP, and AK in rhizosphere soil and the top 15 cm soil layer, as well as on OM content in each soil layer by biochar. Meanwhile, this decreased the levels of Cd and As in all soil layers. Simultaneous application significantly enriched bacteria (mainly including Proteobacteria, Acidobacteriota, and Firmicutes) and fungi (Ascomycota) that promoted plant resistance to PTEs, and improved plant growth. RDA analysis and Spearman correlation analysis showed that soil physicochemical properties (CEC, OM) and PTEs (T-Cd, T-As) contents were the main factors affecting the structure of bacterial and fungal communities. Thus, the combined application of biochar and TH is highly recommended as an effective approach to improve the phytoremediation efficiency of Brassica juncea on Cd- or As-contaminated soil. This study provides new insights into the response mechanism of the combined application of biochar and Endophytic fungi (Trichoderma harzianum) in regulating the phytoremediation efficiency and rhizosphere microbial communities of soil. Further research is needed to explore the mechanisms of promoting plant remediation by combining biochar and PGPF.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/plants12162939/s1, Table S1: Soil physicochemical characteristics in the experimental field; Table S2: Biochar physicochemical characteristics in this experiment; Table S3: Plant shoots and roots biomass (mean ± S.E., n = 3); Table S4: Variation in soil physicochemical properties (pH, EC, CEC, AN, AP, AK, OM) in the inter-root: 15 cm, 30 cm, and 45 cm soil layers under different restoration treatments (mean ± S.E., n = 3); Table S5: Rhizosphere soil microbial community α-diversity index (mean ± S.E., n = 3); Figure S1: Rarefaction curves of bacteria (A) and fungi (B) at the 97% sequence similarity level; Figure S2: Heatmap of the 30 most abundant genera for bacteria (A) and fungi (B) across all 6 samples in rhizosphere soil; Figure S3: Pearson correlation between rhizosphere soil physicochemical properties, enzyme activity and plant potentially toxic elements accumulation.

Author Contributions

Conceptualization, S.Y. and B.Z.; data curation, S.Y. and B.Z.; formal analysis, S.Y., M.D. and H.Z.; funding acquisition, B.Z.; investigation, S.Y., T.C., Z.W., M.W., Q.Y. and Y.L.; methodology, S.Y., B.Z. and M.D.; project administration, X.C.; resources, B.Z., M.D., X.C. and W.C.; supervision, B.Z., M.D., X.C., H.W. and M.W.; visualization, S.Y., M.D., T.C. and Z.W.; writing—original draft, S.Y.; writing—review and editing, B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (41977007; 52222903), a research project of the State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China (QJNY-2019-01).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors gratefully acknowledge the Xi’an University of Technology for providing lab facilitation and experimental measurements.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Man, X.; Huang, H.; Chen, F.; Gu, Y.; Liang, R.; Wang, B.; Jordan, R.W.; Jiang, S. Anthropogenic impacts on the temporal variation of heavy metals in Daya Bay (South China). Mar. Pollut. Bull. 2022, 185, 114209. [Google Scholar] [CrossRef] [PubMed]
  2. Li, P.; Zhang, H.; Xu, L.; Li, H.; Hu, F.; Jiao, J. Investigation and evaluation of soil heavy metals in a wheat-maize cropping system in upland China. J. Agro-Environ. Sci. 2022, 41, 46–54. [Google Scholar]
  3. Kumar, A.; Subrahmanyam, G.; Mondal, R.; Cabral-Pinto, M.M.S.; Shabnam, A.A.; Jigyasu, D.K.; Malyan, S.K.; Fagodiya, R.K.; Khan, S.A.; Yu, Z.-G. Bio-remediation approaches for alleviation of cadmium contamination in natural resources. Chemosphere 2021, 268, 128855. [Google Scholar] [CrossRef]
  4. Mao, C.; Song, Y.; Chen, L.; Ji, J.; Li, J.; Yuan, X.; Yang, Z.; Ayoko, G.A.; Frost, R.L.; Theiss, F. Human health risks of heavy metals in paddy rice based on transfer characteristics of heavy metals from soil to rice. Catena 2019, 175, 339–348. [Google Scholar] [CrossRef]
  5. Saxena, G.; Purchase, D.; Mulla, S.I.; Saratale, G.D.; Bharagava, R.N. Phytoremediation of Heavy Metal-Contaminated Sites: Eco-environmental Concerns, Field Studies, Sustainability Issues, and Future Prospects. Rev. Environ. Contam. Toxicol. 2020, 249, 71–131. [Google Scholar] [PubMed] [Green Version]
  6. Liu, C.; Lin, H.; Dong, Y.; Li, B. Paper Increase of P and Cd bioavailability in the rhizosphere by endophytes promoted phytoremediation efficiency of Phytolacca acinosa. J. Hazard. Mater. 2022, 431, 128546. [Google Scholar] [CrossRef] [PubMed]
  7. Xing, J.; Li, C.; Li, Z.; Li, W.; Fang, A.; Li, A. Submerged macrophytes mediated remediation of molybdenum-contaminated sediments. Environ. Sci. Pollut. Res. Int. 2023, 30, 48962–48971. [Google Scholar] [CrossRef]
  8. Zadel, U.; Nesme, J.; Michalke, B.; Vestergaard, G.; Plaza, G.A.; Schroder, P.; Radl, V.; Schloter, M. Changes induced by heavy metals in the plant-associated microbiome of Miscanthus × giganteus. Sci. Total Environ. 2020, 711, 134433. [Google Scholar] [CrossRef]
  9. Tan, P.; Zeng, C.; Wan, C.; Liu, Z.; Dong, X.; Peng, J.; Lin, H.; Li, M.; Liu, Z.; Yan, M. Metabolic Profiles of Brassica juncea Roots in Response to Cadmium Stress. Metabolites 2021, 11, 383. [Google Scholar] [CrossRef]
  10. Dhaliwal, S.S.; Sharma, V.; Shukla, A.K.; Taneja, P.K.; Kaur, L.; Verma, V.; Kaur, M.; Kaur, J. Exploration of Cd transformations in Cd spiked and EDTA-chelated soil for phytoextraction by Brassica species. Environ. Geochem. Health 2022, 44, 1–13. [Google Scholar] [CrossRef]
  11. Zunaidi, A.A.; Lim, L.H.; Metali, F. Transfer of heavy metals from soils to curly mustard (Brassica juncea (L.) Czern.) grown in an agricultural farm in Brunei Darussalam. Heliyon 2021, 7, e07945. [Google Scholar] [CrossRef] [PubMed]
  12. Ying, R.; Xia, B.; Zeng, X.; Qiu, R.; Tang, Y.; Hu, Z. Adsorption of Cadmium by Brassica juncea (L.) Czern. and Brassica pekinensis (Lour.) Rupr in Pot Experiment. Sustainability 2022, 14, 429. [Google Scholar] [CrossRef]
  13. Ali, I.; Khan, M.J.; Shah, A.; Deeba, F.; Hussain, H.; Yazdan, F.; Khan, M.U.; Khan, M.D. Screening of various Brassica species for phytoremediation of heavy metals-contaminated soil of Lakki Marwat, Pakistan. Environ. Sci. Pollut. Res. 2022, 29, 37765–37776. [Google Scholar] [CrossRef] [PubMed]
  14. Bhatt, P.; Chaudhary, P.; Ahmad, S.; Bhatt, K.; Chandra, D.; Chen, S. Chapter 2—Recent advances in the application of microbial inoculants in the phytoremediation of xenobiotic compounds. In Unravelling Plant-Microbe Synergy; Chandra, D., Bhatt, P., Eds.; Academic Press: Cambridge, MA, USA, 2023; pp. 37–48. [Google Scholar]
  15. Srivastava, S.; Anand, V.; Singh, P.; Roy, A.; Pallavi, S.; Bist, V.; Kaur, J.; Srivastava, S.; Katiyar, R.; Srivastava, S. Chapter 16—Microbial systems as a source of novel genes for enhanced phytoremediation of contaminated soils. In Microbe Mediated Remediation of Environmental Contaminants; Kumar, A., Singh, V.K., Singh, P., Mishra, V.K., Eds.; Woodhead Publishing: Sawston, UK, 2021; pp. 177–198. [Google Scholar]
  16. Liu, J.; Sun, K.; Han, L. Effect of biochar on soil heavy metal speciation and bioavailability: A review. Environ. Chem. 2021, 40, 1643–1658. [Google Scholar]
  17. Li, F.; Shen, W.; Wu, X.; Hu, L.; Wang, W.; Wang, J. Remediation of Heavy Metal Contaminated Soil by Passivation of Biochar Complex Minerals. Chin. J. Soil Sci. 2020, 51, 195–200. [Google Scholar]
  18. Gasco, G.; Alvarez, M.L.; Paz-Ferreiro, J.; Mendez, A. Combining phytoextraction by Brassica napus and biochar amendment for the remediation of a mining soil in Riotinto (Spain). Chemosphere 2019, 231, 562–570. [Google Scholar] [CrossRef]
  19. Glaser, B.; Lehr, V.I. Biochar effects on phosphorus availability in agricultural soils: A meta-analysis. Sci. Rep. 2019, 9, 9338. [Google Scholar] [CrossRef] [Green Version]
  20. Wu, X.; Wang, D.; Riaz, M.; Zhang, L.; Jiang, C. Investigating the effect of biochar on the potential of increasing cotton yield, potassium efficiency and soil environment. Ecotoxicol. Environ. Saf. 2019, 182, 109451. [Google Scholar] [CrossRef]
  21. Liu, Y.; Zhou, J.; Sun, D.; Chen, H.; Qin, J.; Chen, G.; Qiu, R. Polyaspartic acid assisted-phytoremediation of cadmium-contaminated farmland: Phytoextraction efficiency, soil quality, and rhizosphere microbial community. Sci. Total Environ. 2023, 862, 160736. [Google Scholar] [CrossRef]
  22. Wu, R.; Long, M.; Tai, X.; Wang, J.; Lu, Y.; Sun, X.; Tang, D.; Sun, L. Microbiological inoculation with and without biochar reduces the bioavailability of heavy metals by microbial correlation in pig manure composting. Ecotoxicol. Environ. Saf. 2022, 248, 114294. [Google Scholar] [CrossRef]
  23. Geetha, N.; Sunilkumar, C.R.; Bhavya, G.; Nandini, B.; Abhijith, P.; Satapute, P.; Shetty, H.S.; Govarthanan, M.; Jogaiah, S. Warhorses in soil bioremediation: Seed biopriming with PGPF secretome to phytostimulate crop health under heavy metal stress. Environ. Res. 2023, 216, 114498. [Google Scholar] [CrossRef] [PubMed]
  24. Elsharkawy, M.M. Suppression of Potato virus Y Infection in Tobacco by Plant Growth Promoting Fungi. Egypt. J. Biol. Pest Control 2016, 26, 695–700. [Google Scholar]
  25. Wonglom, P.; Ito, S.-I.; Sunpapao, A. Volatile organic compounds emitted from endophytic fungus Trichoderma asperellum T1 mediate antifungal activity, defense response and promote plant growth in lettuce (Lactuca sativa). Fungal Ecol. 2020, 43, 100867. [Google Scholar] [CrossRef]
  26. Babu, A.G.; Shim, J.; Bang, K.-S.; Shea, P.J.; Oh, B.-T. Trichoderma virens PDR-28: A heavy metal-tolerant and plant growth-promoting fungus for remediation and bioenergy crop production on mine tailing soil. J. Environ. Manag. 2014, 132, 129–134. [Google Scholar] [CrossRef]
  27. Mohammadian, E.; Ahari, A.B.; Arzanlou, M.; Oustan, S.; Khazaei, S.H. Tolerance to heavy metals in filamentous fungi isolated from contaminated mining soils in the Zanjan Province, Iran. Chemosphere 2017, 185, 290–296. [Google Scholar] [CrossRef]
  28. Ghosh, S.K.; Panja, A. 13—Signatures of signaling pathways underlying plant-growth promotion by fungi. In Biocontrol Agents and Secondary Metabolites; Jogaiah, S., Ed.; Woodhead Publishing: Sawston, UK, 2021; pp. 321–346. [Google Scholar]
  29. Khalid, M.; Ur-Rahman, S.; Hassani, D.; Hayat, K.; Zhou, P.; Hui, N. Advances in fungal-assisted phytoremediation of heavy metals: A review. Pedosphere 2021, 31, 475–495. [Google Scholar] [CrossRef]
  30. Harindintwali, J.D.; Zhou, J.; Yang, W.; Gu, Q.; Yu, X. Biochar-bacteria-plant partnerships: Eco-solutions for tackling heavy metal pollution. Ecotoxicol. Environ. Saf. 2020, 204, 111020. [Google Scholar] [CrossRef]
  31. Tu, C.; Wei, J.; Guan, F.; Liu, Y.; Sun, Y.; Luo, Y. Biochar and bacteria inoculated biochar enhanced Cd and Cu immobilization and enzymatic activity in a polluted soil. Environ. Int. 2020, 137, 105576. [Google Scholar] [CrossRef]
  32. Chen, H.; Zhang, J.; Tang, L.; Su, M.; Tian, D.; Zhang, L.; Li, Z.; Hu, S. Enhanced Pb immobilization via the combination of biochar and phosphate solubilizing bacteria. Environ. Int. 2019, 127, 395–401. [Google Scholar] [CrossRef]
  33. Nedjimi, B. Phytoremediation: A sustainable environmental technology for heavy metals decontamination. SN Appl. Sci. 2021, 3, 286. [Google Scholar] [CrossRef]
  34. Guo, X.R.; Xu, J.Y.; He, D.Y.; Bu, D.R.; Lu, Y.Y.; Zhao, Y.X.; Chen, Y.R.; Tian, X.F. Combining paecilomyces variotii extracts and biochar for the remediation of alkaline Cd-contaminated soil. Biomass Convers. Biorefinery 2022, 12, 1–10. [Google Scholar] [CrossRef]
  35. Li, Z.; Wang, P.; Yue, X.; Wang, J.; Ren, B.; Qu, L.; Han, H. Effects of Bacillus thuringiensis HC-2 Combined with Biochar on the Growth and Cd and Pb Accumulation of Radish in a Heavy Metal-Contaminated Farmland under Field Conditions. Int. J. Environ. Res. Public Health 2019, 16, 3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Al-Askar, A.A.; Ezzat, A.S.; Ghoneem, K.M.; Saber, W.I.A. Trichoderma harzianum WKY5 and its Gibberellic Acid Control of Rhizoctonia solani, Improve Sprouting, Growth and Productivity of Potato. Egypt. J. Biol. Pest Control 2016, 26, 787–796. [Google Scholar]
  37. Mendes, J.B.S.; da Costa Neto, V.P.; de Sousa, C.D.A.; de Carvalho Filho, M.R.; Rodrigues, A.C.; Bonifacio, A. Trichoderma and bradyrhizobia act synergistically and enhance the growth rate, biomass and photosynthetic pigments of cowpea (Vigna unguiculata) grown in controlled conditions. Symbiosis 2020, 80, 133–143. [Google Scholar] [CrossRef]
  38. Zheng, J.; Xie, X.; Li, C.; Wang, H.; Yu, Y.; Huang, B. Regulation mechanism of plant response to heavy metal stress mediated by endophytic fungi. Int. J. Phytoremediation 2023, 25, 1–18. [Google Scholar] [CrossRef]
  39. Malik, A.; Butt, T.A.; Naqvi, S.T.A.; Yousaf, S.; Qureshi, M.K.; Zafar, M.I.; Farooq, G.; Nawaz, I.; Iqbal, M. Lead tolerant endophyte Trametes hirsuta improved the growth and lead accumulation in the vegetative parts of Triticum aestivum L. Heliyon 2020, 6, e04188. [Google Scholar] [CrossRef]
  40. El-Shahir, A.A.; El-Tayeh, N.A.; Ali, O.M.; Abdel Latef, A.A.H.; Loutfy, N. The Effect of Endophytic Talaromyces pinophilus on Growth, Absorption and Accumulation of Heavy Metals of Triticum aestivum Grown on Sandy Soil Amended by Sewage Sludge. Plants 2021, 10, 2659. [Google Scholar] [CrossRef]
  41. Tang, J.; Zhang, L.; Zhang, J.; Ren, L.; Zhou, Y.; Zheng, Y.; Luo, L.; Yang, Y.; Huang, H.; Chen, A. Physicochemical features, metal availability and enzyme activity in heavy metal-polluted soil remediated by biochar and compost. Sci. Total Environ. 2020, 701, 134751. [Google Scholar] [CrossRef]
  42. Fu, T.; Zhang, B.; Gao, X.; Cui, S.; Guan, C.-Y.; Zhang, Y.; Zhang, B.; Peng, Y. Recent progresses, challenges, and opportunities of carbon-based materials applied in heavy metal polluted soil remediation. Sci. Total Environ. 2023, 856, 158810. [Google Scholar] [CrossRef]
  43. Li, R.; Tan, W.; Wang, G.; Zhao, X.; Dang, Q.; Yu, H.; Xi, B. Nitrogen addition promotes the transformation of heavy metal speciation from bioavailable to organic bound by increasing the turnover time of organic matter: An analysis on soil aggregate level. Environ. Pollut. 2019, 255, 113170. [Google Scholar] [CrossRef]
  44. Bustamante, M.A.; Gomis, M.P.; Perez-Murcia, M.D.; Gangi, D.; Ceglie, F.G.; Paredes, C.; Perez-Espinosa, A.; Bernal, M.P.; Moral, R. Use of livestock waste composts as nursery growing media: Effect of a washing pre-treatment. Sci. Hortic. 2021, 281, 109954. [Google Scholar] [CrossRef]
  45. Liu, C.; Li, B.; Dong, Y.; Lin, H. Endophyte colonization enhanced cadmium phytoremediation by improving endosphere and rhizosphere microecology characteristics. J. Hazard. Mater. 2022, 434, 128829. [Google Scholar] [CrossRef] [PubMed]
  46. Hong, Y.K.; Kim, J.W.; Lee, S.P.; Yang, J.E.; Kim, S.C. Heavy metal remediation in soil with chemical amendments and its impact on activity of antioxidant enzymes in Lettuce (Lactuca sativa) and soil enzymes. Appl. Biol. Chem. 2020, 63, 42. [Google Scholar] [CrossRef]
  47. Sadeghi, F.; Samsampour, D.; Seyahooei, M.A.; Bagheri, A.; Soltani, J. Fungal endophytes alleviate drought-induced oxidative stress in mandarin (Citrus reticulata L.): Toward regulating the ascorbate-glutathione cycle. Sci. Hortic. 2020, 261, 108991. [Google Scholar] [CrossRef]
  48. Pal, M.; Janda, T.; Szalai, G. Interactions between plant hormones and thiol-related heavy metal chelators. Plant Growth Regul. 2018, 85, 173–185. [Google Scholar] [CrossRef] [Green Version]
  49. Jiao, P.-Y.; Guo, W.; Chen, Z.-L.; Liu, X.; Hu, Y.-L.; Wang, Y.-Z. Soil Enzyme Stoichiometric Characteristics of Pinus massoniana Plantations at Different Stand Ages in Mid-subtropical Areas. Huan Jing Ke Xue 2022, 43, 1059–1068. [Google Scholar] [CrossRef] [PubMed]
  50. Gutierrez, C.; Fernandez, C.; Escuer, M.; Campos-Herrera, R.; Beltran Rodriguez, M.E.; Carbonell, G.; Rodriguez Martin, J.A. Effect of soil properties, heavy metals and emerging contaminants in the soil nematodes diversity. Environ. Pollut. 2016, 213, 184–194. [Google Scholar] [CrossRef]
  51. Stella, T.; Covino, S.; Cvancarova, M.; Filipova, A.; Petruccioli, M.; D’Annibale, A.; Cajthaml, T. Bioremediation of long-term PCB-contaminated soil by white-rot fungi. J. Hazard. Mater. 2017, 324, 701–710. [Google Scholar] [CrossRef]
  52. Qi, Q.; Hu, C.; Lin, J.; Wang, X.; Tang, C.; Dai, Z.; Xu, J. Contamination with multiple heavy metals decreases microbial diversity and favors generalists as the keystones in microbial occurrence networks. Environ. Pollut. 2022, 306, 119406. [Google Scholar] [CrossRef]
  53. Ma, J.; Ullah, S.; Niu, A.; Liao, Z.; Qin, Q.; Xu, S.; Lin, C. Heavy metal pollution increases CH4 and decreases CO2 emissions due to soil microbial changes in a mangrove wetland: Microcosm experiment and field examination. Chemosphere 2021, 269, 128735. [Google Scholar] [CrossRef]
  54. Kou, B.; He, Y.; Wang, Y.; Qu, C.; Tang, J.; Wu, Y.; Tan, W.; Yuan, Y.; Yu, T. The relationships between heavy metals and bacterial communities in a coal gangue site. Environ. Pollut. 2023, 322, 121136. [Google Scholar] [CrossRef]
  55. Ma, S.; Qiao, L.; Liu, X.; Zhang, S.; Zhang, L.; Qiu, Z.; Yu, C. Microbial community succession in soils under long-term heavy metal stress from community diversity-structure to KEGG function pathways. Environ. Res. 2022, 214, 113822. [Google Scholar] [CrossRef]
  56. Al Disi, Z.; Al-Ghouti, M.A.; Zouari, N. Investigating the simultaneous removal of hydrocarbons and heavy metals by highly adapted Bacillus and Pseudomonas strains. Environ. Technol. Innov. 2022, 27, 102513. [Google Scholar] [CrossRef]
  57. Lin, Y.; Xiao, W.; Ye, Y.; Wu, C.; Hu, Y.; Shi, H. Adaptation of soil fungi to heavy metal contamination in paddy fields-a case study in eastern China. Environ. Sci. Pollut. Res. 2020, 27, 27819–27830. [Google Scholar] [CrossRef]
  58. De Agostini, A.; Caltagirone, C.; Caredda, A.; Cicatelli, A.; Cogoni, A.; Farci, D.; Guarino, F.; Garau, A.; Labra, M.; Lussu, M.; et al. Heavy metal tolerance of orchid populations growing on abandoned mine tailings: A case study in Sardinia Island (Italy). Ecotoxicol. Environ. Saf. 2020, 189, 110018. [Google Scholar] [CrossRef]
  59. Wijayawardene, N.N.; Hyde, K.D.; Dai, D.Q.; Tang, L.Z.; Aptroot, A.; Castaneda-Ruiz, R.F.; Druzhinina, I.S.; Cai, F.; Ekanayaka, A.H.; Erdogdu, M.; et al. A dynamic portal for a community-driven, continuously updated classification of Fungi and fungus-like organisms: Outlineoffungi.org. Mycosphere 2020, 11, 1514–1526. [Google Scholar] [CrossRef]
  60. Chen, Y.; Chen, Y.; Li, Y.; Wu, Y.; Zeng, Z.; Xu, R.; Wang, S.; Li, H.; Zhang, J. Changes of heavy metal fractions during co-composting of agricultural waste and river sediment with inoculation of Phanerochaete chrysosporium. J. Hazard. Mater. 2019, 378, 120757. [Google Scholar] [CrossRef]
  61. Wang, L.; Lin, H.; Dong, Y.; Li, B.; He, Y. Effects of endophytes inoculation on rhizosphere and endosphere microecology of Indian mustard (Brassica juncea) grown in vanadium-contaminated soil and its enhancement on phytoremediation. Chemosphere 2020, 240, 124891. [Google Scholar] [CrossRef]
  62. Li, Y.; Shi, X.; Ling, Q.; Li, S.; Wei, J.; Xin, M.; Xie, D.; Chen, X.; Liu, K.; Yu, F. Bacterial extracellular polymeric substances: Impact on soil microbial community composition and their potential role in heavy metal-contaminated soil. Ecotoxicol. Environ. Saf. 2022, 240, 113701. [Google Scholar] [CrossRef]
  63. Li, S.; Li, G.; Huang, X.; Chen, Y.; Lv, C.; Bai, L.; Zhang, K.; He, H.; Dai, J. Cultivar-specific response of rhizosphere bacterial community to uptake of cadmium and mineral elements in rice (Oryza sativa L.). Ecotoxicol. Environ. Saf. 2023, 249, 114403. [Google Scholar] [CrossRef]
  64. Lin, G.; Wang, K.; He, X.; Yang, Z.; Wang, L. Characterization of physicochemical parameters and bioavailable heavy metals and their interactions with microbial community in arsenic-contaminated soils and sediments. Environ. Sci. Pollut. Res. 2022, 29, 49672–49683. [Google Scholar] [CrossRef]
  65. Zhu, D.; Wang, Z.; Zhang, Z. Effects of heavy metal pollution and soil physicochemical properties on the Sphagnum farmland soil microbial community structure in Southern Guizhou, China. Int. J. Phytoremediation 2023, 25, 1–12. [Google Scholar] [CrossRef]
  66. Chen, Y.; Jiang, Y.; Huang, H.; Mou, L.; Ru, J.; Zhao, J.; Xiao, S. Long-term and high-concentration heavy-metal contamination strongly influences the microbiome and functional genes in Yellow River sediments. Sci. Total Environ. 2018, 637, 1400–1412. [Google Scholar] [CrossRef]
  67. Yadav, A.N.; Kour, D.; Kaur, T.; Devi, R.; Yadav, A. Endophytic fungal communities and their biotechnological implications for agro-environmental sustainability. Folia Microbiol. 2022, 67, 203–232. [Google Scholar] [CrossRef]
  68. Yuan, F.; Li, K.-Y.; Yang, H.; Deng, C.-J.; Liang, H.; Song, L.-H. Effects of Biochar Application on Yellow Soil Nutrients and Enzyme Activities. Huan Jing Ke Xue 2022, 43, 4655–4661. [Google Scholar] [CrossRef]
  69. Hong, S.B.; Piao, S.L.; Chen, A.P.; Liu, Y.W.; Liu, L.L.; Peng, S.S.; Sardans, J.; Sun, Y.; Penuelas, J.; Zeng, H. Afforestation neutralizes soil pH. Nat. Commun. 2018, 9, 7. [Google Scholar] [CrossRef] [Green Version]
  70. Wu, B.; Luo, S.; Luo, H.; Huang, H.; Xu, F.; Feng, S.; Xu, H. Improved phytoremediation of heavy metal contaminated soils by Miscanthus floridulus under a varied rhizosphere ecological characteristic. Sci. Total Environ. 2022, 808, 151995. [Google Scholar] [CrossRef]
  71. Page, A.L.; Page, A.L. Methods of Soil Analysis: Chemical and Microbiological Proerpteis; Amen Society of Agronomy: Madison, WI, USA, 1982. [Google Scholar]
  72. Liu, H.; Wang, C.; Xie, Y.; Luo, Y.; Sheng, M.; Xu, F.; Xu, H. Ecological responses of soil microbial abundance and diversity to cadmium and soil properties in farmland around an enterprise-intensive region. J. Hazard. Mater. 2020, 392, 122478. [Google Scholar] [CrossRef]
  73. Liu, H.; Xu, F.; Xie, Y.; Wang, C.; Zhang, A.; Li, L.; Xu, H. Effect of modified coconut shell biochar on availability of heavy metals and biochemical characteristics of soil in multiple heavy metals contaminated soil. Sci. Total Environ. 2018, 645, 702–709. [Google Scholar] [CrossRef]
  74. Wijayawardena, M.A.A.; Megharaj, M.; Naidu, R. Chapter Three—Exposure, Toxicity, Health Impacts, and Bioavailability of Heavy Metal Mixtures. In Advances in Agronomy; Sparks, D.L., Ed.; Academic Press: Cambridge, MA, USA, 2016; Volume 138, pp. 175–234. [Google Scholar]
  75. Bogati, K.A.; Golinska, P.; Sewerniak, P.; Burkowska-But, A.; Walczak, M. Deciphering the Impact of Induced Drought in Agriculture Soils: Changes in Microbial Community Structure, Enzymatic and Metabolic Diversity. Agronomy 2023, 13, 1417. [Google Scholar] [CrossRef]
  76. Tabatabai, M.A.; Bremner, J.M. Use of p-nitrophenyl phosphate for assay of soil phosphatase activity. Soil Biol. Biochem. 1969, 1, 301–307. [Google Scholar] [CrossRef]
  77. Wei, J.; Zhang, H.; Wang, Z.; Wang, H. Effect of Morchella Mycelium on Soil Enzymatic Activity. Biotechnol. Bull. 2020, 36, 211–217. [Google Scholar]
  78. Sedlak, J.; Lindsay, R.H. Estimation of total, protein-bound, and nonprotein sulfhydryl groups in tissue with Ellman’s reagent. Anal. Biochem. 1968, 25, 192–205. [Google Scholar] [CrossRef]
  79. Shi, Z.; Wen, M.; Zhang, J.; Tang, Z.; Wang, C. Effect of phenanthrene on the biological characteristics of earthworm casts and their relationships with digestive and anti-oxidative systems. Ecotoxicol. Environ. Saf. 2020, 193, 110359. [Google Scholar] [CrossRef]
  80. Wang, W.; Lu, T.; Liu, L.; Yang, X.; Li, X.; Qiu, G. Combined remediation effects of biochar, zeolite and humus on Cd-contaminated weakly alkaline soils in wheat farmland. Chemosphere 2022, 302, 134851. [Google Scholar] [CrossRef]
  81. Hu, J.; Wu, S.; Wu, F.; Leung, H.M.; Lin, X.; Wong, M.H. Arbuscular mycorrhizal fungi enhance both absorption and stabilization of Cd by Alfred stonecrop (Sedum alfredii Hance) and perennial ryegrass (Lolium perenne L.) in a Cd-contaminated acidic soil. Chemosphere 2013, 93, 1359–1365. [Google Scholar] [CrossRef]
  82. Sigrist, M.; Albertengo, A.; Beldoménico, H.; Tudino, M. Determination of As(III) and total inorganic As in water samples using an on-line solid phase extraction and flow injection hydride generation atomic absorption spectrometry. J. Hazard. Mater. 2011, 188, 311–318. [Google Scholar] [CrossRef]
  83. Shiowatana, J.; McLaren, R.G.; Chanmekha, N.; Samphao, A. Fractionation of Arsenic in Soil by a Continuous-Flow Sequential Extraction Method. J. Environ. Qual. 2001, 30, 1940–1949. [Google Scholar] [CrossRef] [Green Version]
  84. Tüzen, M. Determination of heavy metals in soil, mushroom and plant samples by atomic absorption spectrometry. Microchem. J. 2003, 74, 289–297. [Google Scholar] [CrossRef]
  85. Zhang, X.; Yu, J.; Huang, Z.; Li, H.; Liu, X.; Huang, J.; Zhuo, R.; Wu, Z.; Qin, X.; Gao, Y.; et al. Enhanced Cd phytostabilization and rhizosphere bacterial diversity of Robinia pseudoacacia L. by endophyte Enterobacter sp. YG-14 combined with sludge biochar. Sci. Total Environ. 2021, 787, 147660. [Google Scholar] [CrossRef]
  86. Gao, T.; Wang, X.; Liu, Y.; Wang, H.; Zuo, M.; He, Y.; Li, H.; Li, G.; Li, C.; Lis, X.; et al. Characteristics and diversity of microbial communities in lead-zinc tailings under heavy metal stress in north-west China. Lett. Appl. Microbiol. 2022, 74, 277–287. [Google Scholar] [CrossRef]
Figure 1. The concentration of Cd (A) and As (B) in the shoot and root of Brassica juncea. (C) Plant PTEs accumulation and translocation factor. Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Figure 1. The concentration of Cd (A) and As (B) in the shoot and root of Brassica juncea. (C) Plant PTEs accumulation and translocation factor. Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Plants 12 02939 g001
Figure 2. Cd concentration (A) and As concentration (B) in different soil layers. (C) Availability of Cd in the rhizosphere and 15 cm deep soil layer. Morphological classification of soil As in the rhizosphere soil (D) and the 15 cm depth soil layer (E). Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Figure 2. Cd concentration (A) and As concentration (B) in different soil layers. (C) Availability of Cd in the rhizosphere and 15 cm deep soil layer. Morphological classification of soil As in the rhizosphere soil (D) and the 15 cm depth soil layer (E). Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Plants 12 02939 g002
Figure 3. Changes in the activities of enzymes in the rhizosphere soil: urease (A), sucrase (B), ALP (C). Changes in antioxidant activities of plant shoot and root under different treatments: POD (D), GSH (E), MDA (F). Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Figure 3. Changes in the activities of enzymes in the rhizosphere soil: urease (A), sucrase (B), ALP (C). Changes in antioxidant activities of plant shoot and root under different treatments: POD (D), GSH (E), MDA (F). Different letters represent notable differences between the treatments at p < 0.05. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Plants 12 02939 g003
Figure 4. Relative percentage of the top 10 predominant microorganisms of bacteria (A) and fungi (B) at phylum level (%). Ternary plots of B (B), T (T1, T2), TB (T1B, T2B) samples of the top 10 species in terms of mean abundance at the phylum taxonomic level were selected for analysis: bacteria (C) and fungi (D). The relative abundance (weighted average) of an OTU is shown by the size of its circle. The specified compartments contribute to the total relative proportion, and this determines the circle’s position. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Figure 4. Relative percentage of the top 10 predominant microorganisms of bacteria (A) and fungi (B) at phylum level (%). Ternary plots of B (B), T (T1, T2), TB (T1B, T2B) samples of the top 10 species in terms of mean abundance at the phylum taxonomic level were selected for analysis: bacteria (C) and fungi (D). The relative abundance (weighted average) of an OTU is shown by the size of its circle. The specified compartments contribute to the total relative proportion, and this determines the circle’s position. CK, Brassica juncea alone; B, Brassica juncea + biochar; T1, Brassica juncea + 4.5 g m−2 TH; T2, Brassica juncea + 9 g m−2 TH; T1B, Brassica juncea + biochar + 4.5 g m−2 TH; T2B, Brassica juncea + biochar + 9 g m−2 TH.
Plants 12 02939 g004
Figure 5. Redundancy analysis of physicochemical properties in rhizosphere soil for bacteria (A) and fungi (B). Procrustes analysis was performed to investigate the association between environmental factors and all bacterial (C) and fungal (D) communities, following the NMDS outcomes (using the Bray–Curtis approach). The connection between the composition of the fungal community and environmental factors was demonstrated by the direction and size of the arrows. The ‘p-values’ reveal the consequences of significance tests for disparities among environmental factors and soil microbial communities. p ≤ 0.001 is represented by three symbols (***), 0.001 < p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Figure 5. Redundancy analysis of physicochemical properties in rhizosphere soil for bacteria (A) and fungi (B). Procrustes analysis was performed to investigate the association between environmental factors and all bacterial (C) and fungal (D) communities, following the NMDS outcomes (using the Bray–Curtis approach). The connection between the composition of the fungal community and environmental factors was demonstrated by the direction and size of the arrows. The ‘p-values’ reveal the consequences of significance tests for disparities among environmental factors and soil microbial communities. p ≤ 0.001 is represented by three symbols (***), 0.001 < p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Plants 12 02939 g005
Figure 6. The top 10 most prevalent bacterial (A) and fungal (B) groups were associated with soil physical and chemical characteristics and plant PTE buildup by the method of Spearman correlation heat map with UPGMA clustering. Blue and red illustrate negative and positive associations, respectively. p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Figure 6. The top 10 most prevalent bacterial (A) and fungal (B) groups were associated with soil physical and chemical characteristics and plant PTE buildup by the method of Spearman correlation heat map with UPGMA clustering. Blue and red illustrate negative and positive associations, respectively. p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Plants 12 02939 g006
Figure 7. (A) Pearson significance analysis of the PTEs in different soil horizons, plant, and the validity of PTEs in rhizosphere soil. (B) RDA of soil enzyme operation and soil physical and chemical attributes. (C) RDA of plant PTEs accumulation and soil physical and chemical attributes. (D) RDA of plant enzyme activity and soil physicochemical properties. p ≤ 0.001 is represented by three symbols (***), 0.001 < p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Figure 7. (A) Pearson significance analysis of the PTEs in different soil horizons, plant, and the validity of PTEs in rhizosphere soil. (B) RDA of soil enzyme operation and soil physical and chemical attributes. (C) RDA of plant PTEs accumulation and soil physical and chemical attributes. (D) RDA of plant enzyme activity and soil physicochemical properties. p ≤ 0.001 is represented by three symbols (***), 0.001 < p ≤ 0.01 is demonstrated by two asterisks (**), and 0.01 < p ≤ 0.05 is demonstrated by a single asterisk (*).
Plants 12 02939 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yao, S.; Zhou, B.; Duan, M.; Cao, T.; Wen, Z.; Chen, X.; Wang, H.; Wang, M.; Cheng, W.; Zhu, H.; et al. Combination of Biochar and Trichoderma harzianum Can Improve the Phytoremediation Efficiency of Brassica juncea and the Rhizosphere Micro-Ecology in Cadmium and Arsenic Contaminated Soil. Plants 2023, 12, 2939. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12162939

AMA Style

Yao S, Zhou B, Duan M, Cao T, Wen Z, Chen X, Wang H, Wang M, Cheng W, Zhu H, et al. Combination of Biochar and Trichoderma harzianum Can Improve the Phytoremediation Efficiency of Brassica juncea and the Rhizosphere Micro-Ecology in Cadmium and Arsenic Contaminated Soil. Plants. 2023; 12(16):2939. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12162939

Chicago/Turabian Style

Yao, Shaoxiong, Beibei Zhou, Manli Duan, Tao Cao, Zhaoquan Wen, Xiaopeng Chen, Hui Wang, Min Wang, Wen Cheng, Hongyan Zhu, and et al. 2023. "Combination of Biochar and Trichoderma harzianum Can Improve the Phytoremediation Efficiency of Brassica juncea and the Rhizosphere Micro-Ecology in Cadmium and Arsenic Contaminated Soil" Plants 12, no. 16: 2939. https://0-doi-org.brum.beds.ac.uk/10.3390/plants12162939

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

Article Metrics

Back to TopTop