Next Article in Journal
Comparative Analysis of Phytohormone Biosynthesis Genes Responses to Long-Term High Light in Tolerant and Sensitive Wheat Cultivars
Previous Article in Journal
Shading and Water Addition Alleviate the Elemental Limitations of the Early Restoration Community in a Stressful Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Structural and Functional Responses of Rhizosphere Bacteria to Biodegradable Microplastics in the Presence of Biofertilizers

1
Marine Ecology and Human Factors Assessment Technical Innovation Center of Natural Resources Ministry, Tsinghua Shenzhen International Graduate School, Shenzhen 518055, China
2
Shenzhen Public Platform for Screening and Application of Marine Microbial Resources, Institute for Ocean Engineering, Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
3
Shenzhen Key Laboratory of Advanced Technology for Marine Ecology, Institute for Ocean Engineering, Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
4
Ecological Fertilizer Research Institute, Shenzhen Batian Ecological Engineering Co., Ltd., Shenzhen 518057, China
*
Authors to whom correspondence should be addressed.
Plants 2024, 13(18), 2627; https://doi.org/10.3390/plants13182627
Submission received: 18 August 2024 / Revised: 9 September 2024 / Accepted: 12 September 2024 / Published: 20 September 2024
(This article belongs to the Section Plant–Soil Interactions)

Abstract

:
Biodegradable microplastics (Bio-MPs) are a hot topic in soil research due to their potential to replace conventional microplastics. Biofertilizers are viewed as an alternative to inorganic fertilizers in agriculture due to their potential to enhance crop yields and food safety. The use of both can have direct and indirect effects on rhizosphere microorganisms. However, the influence of the coexistence of “Bio-MPs and biofertilizers” on rhizosphere microbial characteristics remains unclear. We investigated the effects of coexisting biofertilizers and Bio-MPs on the structure, function, and especially the carbon metabolic properties of crop rhizosphere bacteria, using a pot experiment in which polyethylene microplastics (PE-MPs) were used as a reference. The results showed that the existence of both microplastics (MPs) changed the physicochemical properties of the rhizosphere soil. Exposure to MPs also remarkably changed the composition and diversity of rhizosphere bacteria. The network was more complex in the Bio-MPs group. Additionally, metagenomic analyses showed that PE-MPs mainly affected microbial vitamin metabolism. Bio-MPs primarily changed the pathways related to carbon metabolism, such as causing declined carbon fixation/degradation and inhibition of methanogenesis. After partial least squares path model (PLS-PM) analysis, we observed that both materials influenced the rhizosphere environment through the bacterial communities and functions. Despite the degradability of Bio-MPs, our findings confirmed that the coexistence of biofertilizers and Bio-MPs affected the fertility of the rhizosphere. Regardless of the type of plastic, its use in soil requires an objective and scientifically grounded approach.

1. Introduction

Traditional agricultural practices have long relied on inorganic fertilizers to manage soil and increase crop yields. However, the use of these fertilizers has resulted in various negative consequences, including increased soil salinity, reduced soil fertility, an exacerbation of the greenhouse effect, and contribution to water body eutrophication [1]. In response to these challenges, alternative methods that rely on organic inputs to enhance nutrient availability and safeguard the environment have been advocated in recent decades. Biofertilizers are viewed as an alternative to inorganic fertilizers in agriculture due to their potential to enhance crop yields and ensure food safety [2]. These organic products contain specific microorganisms [3] capable of injecting various nutrients into the soil by fixing nitrogen, dissolving or mineralizing phosphate and potassium, releasing plant growth regulators, and producing antibiotics [3]. Serving as microbial inoculants, biofertilizers enhance soil physicochemical properties, microbial community diversity, and plant growth [3]. Useful microbial populations in agriculture include plant growth-promoting rhizobacteria (PGPR), N2-fixing cyanobacteria, plant disease-suppressive bacteria, oil-toxicant degrading microbes, and actinomycetes [4,5].
Environmental contamination by microplastics (MPs) is now considered an emerging threat to soil ecosystems, and agricultural lands act as potential sinks for MPs [6]. Recent reports indicate that microplastic fragments can make up up to 7% of soil by dry weight [7]. Estimates suggest that 43–63 and 30–44 million tons of MPs are introduced annually to European and North American farmlands, respectively [8]. This microplastic contamination in agroecosystems often originates from farm practices, especially when using plastic mulches. Among these mulches, the main production materials are durable polymers, such as polyethylene (PE) and polypropylene (PP), which are not biodegradable and can persist in the environment for decades [9]. Some biodegradable mulches have been developed to achieve greener and more environmentally friendly plastic mulches [10,11]. As an alternative to traditional agricultural mulches, biodegradable polymers such as polylactic acid (PLA) and polybutylene terephthalate (PBAT) are increasingly being used [12]. However, the degradation of many biodegradable polymers under environmental conditions can be slow or incomplete, leading to the accumulation of Bio-MPs in soil [13].
The rhizosphere, which refers to the narrow zone surrounding and influenced by plant roots, is a hotspot for numerous organisms and is considered as one of the most complex ecosystems on Earth [14]. Each gram of root can contain up to 10 [11] rhizosphere microorganisms [15]. Rhizosphere microorganisms can promote plant growth through various mechanisms, including nitrogen fixation [16], promoting the uptake of micronutrients such as Fe [17], releasing cations through mineral weathering [18]. They can also produce antibiotics [14], and regulate the plant’s immune system to protect the plant from pathogen attacks [19]. The coexistence of biofertilizers and MPs is becoming increasingly common due to their parallel appearance in the environment [1]. And the use of both can have direct and indirect effects on rhizosphere microorganisms. Previous research has primarily focused on studying the effects of soil physicochemical properties, soil enzyme activities, and plant growth and development in the presence of biofertilizers or MPs alone. Although some studies have examined the effects on soil microorganisms, few have specifically investigated the rhizosphere microenvironment [20,21,22,23,24]. Therefore, further studies are necessary to investigate the responses of rhizosphere microbial communities in the coexistence of biofertilizers and MPs and enhance our understanding of the impact on the rhizosphere environment. For conventional MPs, MPs can change the features of soil (e.g., pH value, water evaporation, and aggregate stability) [5,25], microbial profiles (e.g., diversity, composition, and organic matter metabolic ability) [26], and the physicochemical properties of crops (e.g., biomass, growth, and root development) [2]. The interaction between Bio-MPs and biofertilizers and their impact on soil microecological characteristics are relatively unknown compared with conventional plastic particles. Little is known about the effects of Bio-MPs (as regards coexistence with biofertilizer) on soil microecological processes, especially the rhizospheric bacterial response.
In this study, we investigated the responses of rhizosphere microbial communities in the coexistence of biofertilizers and Bio-MPs and enhanced our understanding of the direct or indirect impact on the rhizosphere. In addition, biodegradable plastics exhibit a different picture of rhizosphere microorganisms than traditional plastics. We used Bacillus amyloliquefaciens biofertilizer as an example, and two different types of MPs, PE-MPs (the most typical staple material for agricultural films and packaging materials) and Bio-MPs (poly (butylene adipate-co-butylene terephthalate)-poly(lactic acid) (PBAT-PLA) blends, an increasingly popular biodegradable plastic) as the target materials to investigate. We investigated the effects of MPs on the rhizosphere microorganisms of Brassica parachinensis in the context of biofertilizer application. Specifically, we (1) analyzed the changes in the bacterial diversity, composition, and network complexity of the rhizosphere microbial community using 16S rRNA sequencing and (2) analyzed the changes in the functioning (especially the carbon cycle potential) of the rhizosphere microbial community using metagenomic tools. The ultimate goal was to evaluate the impact of degradable plastic particles on the soil microecological environment.

2. Results

2.1. Changes in the Physicochemical Properties of the Rhizosphere

The inclusion of 0.2% and 2.0% Bio-MPs led to a significant increase in SOC (96.87% and 176.81%) and OM (90.60% and 158.03%), which also showed dose dependence (Figure 1a,b). By comparison, PE-MPs increased SOC (27.57%) and OM (20.15%) at only high concentrations (2.0%), and the effects were significantly weaker than those of Bio-MPs (Figure 1a,b). Additionally, the presence of Bio-MPs also caused a significant increase in the pH of the rhizosphere soil, changing it from weakly acidic (6.79) to weakly alkaline (7.92–8.11) compared to the BF group (Figure 1c). For other physicochemical parameters (TN, TP, TK, AHN, and AK), there were slight variations to Bio-MPs and exposed concentration (Figure S2). Both MP materials significantly decreased the AP level compared with the blank group. The decrease was approximately one of 3.5–4 times (Figure S2).

2.2. Changes in Bacterial Communities in the Rhizospheric Environment

The raw sequences were subjected to quality control and sequence screening before being optimized based on the 97% similarity threshold. This resulted in 18,485 bacterial operational taxonomic units (Figure S3a,b). The richness and diversity of the communities in the different treatment groups were assessed using the ACE and Shannon indexes. Compared to the BF group, the ACE index of bacteria of the PE-H group significantly reduced (Figure 2a). Both treatment groups of Bio-MPs also significantly reduced the alpha-diversity of bacteria and followed the concentration-dependent effect (Figure 2a,b). For β-diversity, PERMANOVA revealed significant differences in bacterial community composition structure among the six treatment groups (F = 4.584, p < 0.01) (Figure 2c).
As for the bacterial composition, the community in the rhizosphere varied significantly depending on the treatments. Among all the groups, Proteobacteria was the dominant taxa at the phyla level, followed by Actinobacteria, Firmicutes, and Bacteroidetes (Figure 2d).
A linear discriminant analysis effect size (LEfSe) was performed to clarify the differences among the treatment groups. A total of 83 biomarkers were identified (Figure S3c). The six groups of CK, BF, PE-L, PE-H, Bio-L, and Bio-H had 11, 25, 6, 8, 6, and 27 biomarkers, respectively. The BF group showed significant enrichment of Acidobacteria, Nitrospirae, and Firmicutes. After MP exposure, the indicators were changed. The indicators in the PE groups were Janthinobacterium, Bacillus, and Parasegitibacter. Similarly, the biomarkers in the Bio groups were Comamonas, Cupriavidus, and Burkholderia.
We found that adding plastic particles changed the sensitivity of rhizosphere microorganisms to environmental parameters after conducting a redundancy analysis (Figure S4a). The Bio groups were positively (significantly) correlated with pH, SOC, OM, and TK. However, these factors were negatively (significantly) correlated with the bacterial community in the CK, BF, and PE groups. Meanwhile, AP, AHN, and AK showed a significant negative correlation with Bio groups and a significant positive correlation with CK, BF, and PE groups (Figure S4a). More detailed results of the relative abundance of key bacterial taxa and soil environmental factors are displayed in the heat maps (Figure S4b).

2.3. The Network Profiles of the Rhizosphere Bacterial Community

The co-occurrence model of different treatments (100 high-ranking bacteria at the genus level) was analyzed to reveal rhizosphere bacterial community responses to biofertilizers and MPs. The main topological parameters include edges, nodes, the percentage of positive/negative edges, average path length, graphic density, network diameter, average degree, and modularity, as summarized in Figure 3. Overall, the percentage of positive edges (56.25–86.07%) outnumbered the negative edges in all test groups (13.93–43.75%) (Figure 3), indicating that bacterial connection types tended to be positively correlated. The graphic density and the number of edges in the BF group decreased by 38.2% and 35.9%, respectively, compared to group CK. The percentage of positive edges decreased by 19.4%, while the percentage of negative edges increased by 120%. These findings suggest that the addition of biofertilizers weakened the linkage of the rhizosphere bacterial community, resulting in low network complexity. Additionally, the average degree and clustering coefficient decreased by 37.2% and 2.0%, respectively, as the network aggregation decreased (Figure 3). PE-H and Bio-H had opposite effects compared to BF. The graph density, number of edges, average degree, and clustering coefficient of PE-H decreased by 12.5%, 14.4%, 13.5%, and 1.3%, respectively. In contrast, the parameters mentioned above increased in the Bio-H group by 106.8%, 112.8%, 110.7%, and 19.4%, respectively (Figure 3). These results suggest that the presence of 2.0% PE-MPs reduces the complexity and density of the network, whereas 2.0% Bio-MPs enhances it.

2.4. Rhizosphere Bacterial Function Based on Metagenomic Analysis

Comparative metagenomic analyses were conducted on the CK, BF, PE-H, and Bio-H groups to determine the effects of MPs present with biofertilizers on rhizosphere bacterial function. The gene sets from each treatment group were annotated using the KEGG primary and secondary functional classification pathways (Figure S5a,b). Figure 4a demonstrates a reduction in the relative abundance of most secondary functions in the Bio-H group, with a stronger reduction than in PE-H. These functions include “amino acid metabolism”, “carbohydrate metabolism”, “lipid metabolism”, and “energy metabolism”. In contrast, “cell motility”, “cell growth and death”, “cellular community and signal transduction”, and “membrane transport” were upregulated (Figure 4a). The clustering results indicated that the BF and PE-H groups had similar features. However, treated Bio-H and PE-H showed the greatest difference in functional composition (Figure 4a).
In addition to the common core functions, representative key functions (LDA score > 2.5) for each group were also determined using LEfSe analysis. Within the KEGG tertiary functional classification pathways, the BF group was enriched in “oxidative phosphorylation” and “metabolic pathways”, with only “fatty acid metabolism” enriched in PE-H. Similarly, the Bio-H group showed enrichment in four pathways: “fatty acid degradation”, “benzoate degradation”, “valine-leucine and isoleucine degradation”, and “two-component system” (Figure 4b).
Considering the metabolic potential of rhizosphere microorganisms under the influence of plastic particles (PE-MPs and Bio-MPs), these enriched functional genes are closely related to carbon metabolism. Therefore, we researched carbon cycle genes in the rhizosphere environment further. Fifteen pathways related to carbon cycle genes (e.g., carbon fixation, central metabolism, and methane metabolism) were identified in the rhizosphere across all treatment groups (Figure S5c). A total of 184 genes (KO numbers) are involved in the metabolic pathways of the carbon cycle, with some KO numbers being associated with multiple pathways. The tricarboxylic acid cycle was the predominant pathway, constituting 17.39% of the overall carbon metabolism function (Figure S5c).
Figure 4c shows that the distribution of genes related to the carbon cycle varies across different groups. Biofertilizers partially downregulated the expression of carbon cycle genes. The introduction of MPs also results in alterations to the carbon cycle. The expression of carbon cycle genes was upregulated in PE-H. At the same time, many genes were further downregulated in Bio-H (Figure 4c). The pathways of carbon fixation and central metabolism did not show significant differences among the four groups (Figure 4d). However, the presence of biofertilizers and MPs had a significant impact on the pathways of methane metabolism. Compared to the CK group, the BF and PE-H groups significantly increased methane production, whereas BF-Bio20 had a significant inhibitory effect on the conversion of methane to methanol (Figure 4d).
A partial (geographic distance-corrected) Mantel test was performed to explore the correlation between carbon cycle genes and environmental factors and to clarify the contribution of environmental factors to carbon metabolism (Figure 5). TK had a significantly strong correlation with functional genes related to methane metabolism in terms of taxonomic and functional composition (p < 0.05). Significant positive correlations were also observed between TK-TN, pH-SOC, pH-OM, SOC-OM, and AK-AHN. Bacterial species also contributed to environmental parameters related to carbon metabolism. Figure S6 shows that multiple species of Acidovorax, Nitrospira, Nocardioides, Noviherbaspirillum, Paraburkholderia, Phenylobacterium, Pseudolabrys, Trinickia, and Massilia contribute to at least one carbon cycle pathway. Devosia was found exclusively in the central metabolism of the CK and BF groups. Opitutus was identified only in the central metabolism of the BF group, while Methylobacillus was identified in the methane metabolism of the CK group.

2.5. Partial Least Squares Path Model (PLS-PM) Analysis

PLS-PM analysis was conducted to investigate the potential relationships between biofertilizers and MPs, soil physicochemical properties, bacterial diversity, and related carbon cycling potential. The results (Figure 6) indicate that biofertilizers and MPs had a direct and significant positive correlation with soil pH, SOC, and OM and a negative correlation with soil nutrient content and bacterial diversity, affecting rhizosphere carbon cycling. Simultaneously, the addition of biofertilizers and MPs has a direct impact on the carbon cycle in the rhizosphere. In the natural environment, the coexistence of biofertilizers and MPs will affect the physical and chemical properties of soil bacterial activities in the rhizosphere and drive carbon turnover in the rhizosphere environment.

3. Materials and Methods

3.1. Preparation of Biofertilizers and Characterization of MPs

The B. amyloliquefaciens biofertilizer was provided by the Shenzhen Batian Ecological Science and Technology Co., Ltd., Shenzhen, China [27]. The PE-MPs and PBAT-PLA Bio-MPs were purchased from Huachuang Chemicals Co., Ltd., Dongguan, China. Scanning electron microscopy (Phenom ProX, Eindhoven, The Netherlands), Fourier transform infrared spectroscopy (FT-IR) (Thermo Scientific Nicolet iS-50, Waltham, MA, USA), and laser particle size analysis (Malvern Mastersizer 2000, Malvern, PA, USA) were used to determine the morphology, composition, and particle size of the two types of MPs (Figure S1). The production process of the biofertilizer and information on the MPs can be found in Supplementary Information Texts S1 and S2.

3.2. Exposure Experiment of MPs

A 28-day greenhouse pot experiment was conducted from 28 November to 26 December 2022, at the Shenzhen Batian Ecological Fertilizer Research Center in Guangdong Province, China. The test soil came from farmland in Guangming District that had not been treated with plastic mulch. Table S1 shows the physicochemical properties of the soil at the site.
The biofertilizers used in our potting experiments were referenced in our previous work (Wang et al. [22]). Each pot contained three 14-day-old B. parachinensis plants with the same growth. We used two concentrations of MPs (0.2% and 2.0% w/w, dry weight ratio of MPs to soil), with references from previous studies [25], to simulate MP concentrations in the real environment. The two concentrations represented low- and high-strength mulch use, respectively.
Six treatment groups were designed, including group CK (an untreated control with no biofertilizer or MPs added), group BF (with biofertilizer only), groups PE-L/PE-H (with biofertilizer and 0.2%/2.0% of PE-MPs), and groups Bio-L/Bio-H (with biofertilizer and 0.2%/2.0% of Bio-MPs). Each treatment group comprised six biological replicates (pots) containing 4 kg of test soil and three B. parachinensis seedlings.

3.3. Physicochemical Properties of Rhizosphere Soil

Rhizosphere soil sampling is described in Supplementary Information Text S3. The physicochemical properties of the rhizosphere soil were determined as follows: pH was determined using a potentiometric method with a soil/water ratio of 1:2.5 w/v. Soil organic carbon (SOC) and organic matter (OM) contents were determined using the potassium dichromate oxidation–external heating method [28]. Total nitrogen (TN) content was determined using the Kjeldahl method. The study analyzed the contents of total phosphorus (TP), total potassium (TK), and alkaline hydrolysis nitrogen (AHN) using molybdenum blue colorimetry, flame photometry, and the alkali diffusion method, respectively. The available phosphorus (AP) content was determined using the molybdenum-antimony inverse colorimetric method after soaking the inter-root soil in ammonium hydrogen fluoride and sodium bicarbonate solutions [28]. The available potassium (AK) content was determined using ammonium acetate leaching and flame photometric methods [29].

3.4. DNA Extraction, PCR Amplification, and Metagenomic Analysis

The bacterial communities in the rhizosphere soil of the six treatment groups were characterized using high-throughput sequencing. DNA was extracted from the rhizosphere soil using cetyltrimethylammonium bromide (CTAB). DNA purity and concentration were assessed using 1% agarose gel electrophoresis. The gene amplification of different regions (V3–V4) of bacterial 16S rRNA was performed using primers 341F (CCTAYGGGRBGCASCAG) and 806R (GGACTACNNGGGGTATCTAAT). The DNA extraction process and assessment of microbial diversity are described in Text S4 of the Supplementary Information.
Metagenomic sequencing analysis was performed for bacterial function. The CK, BF, PE-H, and Bio-H groups were collected as targets. Metagenomic sequencing was performed using the NEBNext@ UltraTM DNA Library Prep Kit for Illumina (New England Biolabs, Massachusetts, USA) according to the manufacturer’s recommendations. Briefly, indexing codes were added to the sequences for each set of samples. DNA samples were sonicated to obtain 350 bp fragments, which were then end-polished, a-tailed, and ligated to full-length aptamers. This process was followed by PCR amplification and purification of the PCR products using the AMPure XP system. Library size distribution was analyzed on an Agilent 2100 Bioanalyzer and quantified using real-time PCR. Index-coded samples were clustered on the cBot clustering system according to the manufacturer’s instructions, and the DNA libraries were finally sequenced on the Illumina Novaseq 6000 platform, generating 150 bp paired reads. Detailed sequence quality control and functional assembly information are provided in Supplementary Information Text S5.

3.5. Statistical Analyses

Statistical analyses were performed using IBM SPSS Statistics 25.0, specifically one-way analysis of variance (ANOVA) with a p-value ≤ 0.05 for significance, and implemented in SPSS 25 software (IBM, New York, NY, USA). All experimental parameters were measured at least in triplicate, whose results are presented as the mean ± standard deviation. Data were analyzed and graphically presented using R Studio (4.1.0) and OriginPro 2019. Two alpha diversity indices, comprising ACE and Simpson, were analyzed using the “vegan” and “picante” packages in R (v4.1.0). To further characterize the microbial community, we used OriginPro 2019 (OriginLab, Northampton, USA) software to draw histograms of the relative abundance of bacteria at the phylum level. Comparisons between multiple groups were achieved using the linear discriminant analysis effect size (LEfSe) method to identify biomarkers that differed significantly in abundance between groups (Kruskal–Wallis rank sum test, linear discriminant analysis (LDA) score > 3.5). In addition, redundancy analysis (RDA) was used to assess the effect of soil physicochemical properties on the microbial community, and its significance was determined by Monte Carlo permutation (permu = 999). We conducted a correlation heatmap analysis in R software to investigate the relationships between environmental factors and species composition at the phylum level. Spearman’s coefficient (|ρ| > 0.6) was used to assess the degree of correlation between species at the genus level, these visualized using Gephi 0.10 software. The Circos diagram shows the abundance of KEGG primary and secondary functional classification pathways. Partial Least Squares Path Modeling (PLS-PM) with the “plspm” package in R was used to explore the direct, indirect, and interactive effects of MPs, soil physicochemical characteristics, microbial community diversity, and carbon cycle functional gene abundance. The“psych”and“ggplot2”packages in R software were used for correlation analyses among different treatments.

4. Discussion

4.1. Bacterial Composition Responses to Biofertilizers and MPs

As expected, the addition of biofertilizer resulted in significantly high rhizosphere SOC and OM content, attributed to its 29.27% exogenous organic matter (Text S1). After MP exposure, especially to Bio-MPs, the SOC and OM contents were further increased, likely due to their ability to release endogenous organic matter during degradation. Biodegradable plastics comprising 60–80% of carbon undergo breakdown in the natural environment [30]. Bio-MPs and their degradation products act as supplementary biologically effective carbon, directly affecting soil OM content. In addition, we found that both MP additions significantly increased the pH of the rhizosphere, changing it from weakly acidic to weakly alkaline. This increase in pH is consistent with the findings of Liu et al. [31]. The potential reason was that the addition of MPs led to the competitive adsorption of hydrophobic organic compounds with OM, which in turn altered the acid–base balance of the soil [30].
In addition to the nonbiological parameters (physicochemical factors) of the rhizosphere environment, exposure to MPs also affects the biological parameters of the rhizosphere. The inclusion of high concentrations of MPs, particularly in the Bio group, led to a significant decrease in the richness and diversity of the bacterial community compared to the sole BF group (Figure 2). We found that the negative effect of MPs was significantly greater than the promoting effect of BF when compared with the CK group. These findings are supported by previous research that showed the negative effects of MPs on rhizosphere bacterial diversity and richness [32]. Possible reasons include, as a hydrophobic polymer with a large specific surface area, how MPs can absorb multiple substances in the rhizosphere, which contributes to the enrichment of specific bacteria [30,32]. Additionally, some antioxidants, colorants, and plasticizers added to MPs are released into the environment and act as nutrient sources for certain bacteria [33]. Third, MPs may facilitate colonization by bacteria capable of degrading polymers as carbon sources. Although PE-MPs are less biodegradable and are unlikely to affect soil microorganisms directly, they can drive the selection of specific microbial taxa with potential degradation capabilities [34]. For instance, certain members of Actinobacteria can break down PE-MPs by producing hydrolases [35,36]. When compared to conventional PE-MPs, Bio-MPs decompose faster. In this case, Bio-MPs with rough surfaces release more carbon sources, which is conducive to bacterial colonization [30]. As a result, Bio-MPs have a greater impact on bacterial alpha-diversity. Meanwhile, the beta-diversity was also influenced by exposure to MPs and displayed concentration dependence (Figure 2b). Combined with Figure S4, we speculated that the dose and type of MP are co-drivers of bacterial community composition [37].
Further, significant differences in bacterial composition were found under different treatments, indicating that biofertilizers and MPs co-influence the structure of bacterial communities. Several studies have shown that Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes play critical roles in organic mineralization and C, N, and S cycling [38,39]. The content of Proteobacteria increased slightly after biofertilizer application, with the highest content found in the Bio-H group. Previous research has identified Proteobacteria as the primary decomposer of OM [40]. As the main taxonomic unit in soil samples, Proteobacteria play an integral role in the carbon cycle and tend to colonize nutrient-rich environments [41]. Acidobacteria plays a crucial role in the carbon cycle and can process a variety of carbon sources [41]. However, the abundance of Acidobacteria was significantly reduced in the Bio treatment group. Acidobacteria is a typical acidophilic bacterium that thrives in acidic and oligotrophic ecological niches rather than nutrient-rich alkaline environments [42]. In this study, adding MPs increased the rhizosphere pH and nutrient content, especially for Bio-MPs. This finding provides supplementary evidence for the reduction in Acidobacteria. The LEfSe results also indicate that the exposure of MPs could screen their biomarker (Figure S3c). Considering Bio groups as an example, Comamonas and Burkholderia were the main indicators in the Bio-L and Bio-H groups. Comamonas is a multifunctional bacterium capable of degrading phenols, PAHs, and heterocyclic aromatic hydrocarbons [43]. The enrichment of Comamonas may be related to the biodegradation of PBAT, which is a fully biodegradable aliphatic-aromatic co-polyester in Bio-MPs [44].
Network analysis has been widely used to explore co-occurrence patterns in soil microbial communities [45,46]. This study found that the addition of biofertilizers resulted in a decrease in network complexity and network aggregation, weakening the linkages of rhizosphere bacterial communities. These findings are inconsistent with previous research [47]. Furthermore, the percentage of negative edges increased while positive edges decreased, indicating that competition, rather than cooperation, was the prevailing factor in the bacterial community. The use of biofertilizers was inferred to introduce many beneficial bacteria into the rhizosphere. This condition led to overlapping dominant bacterial ecological niches [48], increasing competition, and decreasing linkages between rhizosphere bacterial communities. MP exposure further reduced the complexity and density of rhizosphere bacteria networks, which is consistent with the study of Wu et al. [49]. As organic pollutants, PE-MPs can impact the symbiotic pattern of rhizosphere bacterial communities. This situation can lead to the creation of synergistic bacterial symbiotic modules that promote the degradation of PE-MPs. Bio-MPs significantly impacted the total number of links in the rhizosphere bacterial network, particularly the percentage of positive edges, facilitating cooperation among bacteria. This outcome may be attributed to the degradable nature of Bio-MPs, resulting in the rhizosphere becoming a resource-rich environment. As a result, bacteria exhibit swift responses to the nutrients in the rhizosphere [35]. Based on these observations, we propose that Bio-MPs can promote the formation of complex bacterial community networks in the rhizosphere.

4.2. Rhizosphere Bacterial Function Responses to Biofertilizers and MPs

While the higher sequencing depth of metagenomic sequencing can determine what the altered microbes are doing, 16S rRNA sequencing can explain how the microbial structure is altered in rooted environments. The breakdown of biofertilizers and MPs by microbes is a long-term biometabolic process involving OM. Therefore, rhizosphere bacteria were annotated with the largest “Metabolism” genes. Our work found that the genes identified were primarily associated with metabolic processes that support bacterial activity, which facilitates nutrient cycling by rhizosphere microbes. Using Bio-MPs resulted in a decrease in metabolic processes and an increase in “cellular processes” and “environmental information processing”. This result suggests that adding Bio-MPs leads to changes in bacterial rhizospheric function.
Considering that carbon is one of the most important elements for plant growth, and through differential metabolic gene analysis, we found that the functional pathway mainly involves carbon metabolism processes. Therefore, we focused on carbon and analyzed its metabolic differences in various groups (Figure 4a,b). The use of biofertilizers induced alterations in carbon cycle genes. After adding MPs, different types had varying effects on rhizosphere bacterial functions. Ma et al. found that a moderate amount of MPs promoted nutrient cycling, but exceeding the threshold concentration, which harmed soil nutrient cycling [50]. Our study found that adding 2.0% PE-MPs upregulated the expression of carbon cycle-related genes, promoting nutrient cycling (Figure 4b). These results are consistent with the study by Ma et al. [50]. However, 2.0% of Bio-MPs suppressed most functional genes related to carbon cycling. The threshold of Bio-MPs is low, meaning they can adversely affect carbon cycling at low concentrations. Upon further evaluation of the relative abundance of the various pathways of the carbon cycle, Bio-MPs had a significant inhibitory effect on methane oxidation (Figure 4d). This inhibition could partially explain the lack of activity observed in Methylobacillus, which utilizes methane oxidation-derived carbon as a nutrient in the Bio group (Figure S6c). The rhizosphere SOC and OM contents were significantly higher in the BF and PE-H groups. This increase may promote the growth of methanogenic bacteria [51], leading to increased methane production (Figure 4d). Ma et al. found a significant correlation between environmental factors, such as SOM and TN, and the abundance of functional genes related to methane metabolism in the rhizosphere after exposure to MPs [50]. This study conducted partial Mantel tests and found a significant and strong correlation between TK and methane metabolism-related functional genes (Figure 5), suggesting that methane metabolism is vulnerable to changes in environmental factors when biofertilizers and MPs coexist.
Changes in epigenetic function are responsible for the composition of microorganisms. We found that different bacterial species responded to different functions after calculating the relative abundance of bacteria involved in the three pathways (Figure S6). For instance, Nitrospira and Paraburkholderia were enriched in carbon fixation in the PE-H and central metabolism in the Bio-H group, respectively. Paraburkholderia is a bacterium capable of C–N fixation [52], while Nitrospira plays an important role in the N–C coupling cycling [53]. Trinickia, Mitsuaria, and Phenylobacterium were relatively enriched in carbon fixation and central metabolism, when treated with Bio-MPs. Trinickia is a PGPR [54], Mitsuaria can reduce the growth of pathogens and act as a biocontrol agent [55], and Phenylobacterium is associated with the degradation of PAHs [56]. These results indicated that the response of bacteria to MPs is typically nonlinear and microbes have structural variability and functional plasticity.

5. Conclusions

This study shows that the coexistence of biofertilizers and microplastics significantly impacts the rhizosphere environment and related bacteria, ultimately affecting the soil carbon cycle. Our work confirmed that MP exposure (i) altered the physicochemical parameters of rhizosphere soil and the bacterial composition, in particular, of the OM decomposing taxa; (ii) decreased the co-occurrence complex and the buffer ability to environmental stresses; and (iii) affected the relative abundance of carbon cycle bacteria, which in turn affects methane metabolism. PLS-PM analyses showed that adding biofertilizers and MPs improved SOC, OM, and the nutrients of soil but reduced its bacterial diversity and carbon cycling ability. Additionally, the linkage analysis shows that the ecological ramifications of exposure to Bio-MPs are bigger than conventional MPs for the rhizosphere of B. parachinensis. The results provide insight into the impact of biofertilizers and MPs on the rhizosphere microenvironment when they coexist. They warn us that although biodegradable plastics are a valuable alternative product, their negative effects in actual soil environments are worth paying attention to.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants13182627/s1, Figure S1: Scanning electron microscope, Fourier transform infrared spectrometer spectra, and laser particle size analyzer images for (a, c, e) PE MPs and (b, d, f) Bio MPs; Figure S2: The effect of MPs on soil property existence with biofertilizer: (a) TN; (b) TP; (c) TK; (d) AHN; (e) AP; and (f) AK. Data are presented as the mean ± SD (n = 6). Different letters (a, b, and c) were statistically different (p < 0.05) according to Duncan’s new multivariate range test; Figure S3: Effects of biofertilizers and MPs on rhizosphere bacterial community structure and composition. A venn diagram illustrating the number of (a) PE MPs and (b) Bio MPs OTUs in rhizosphere soil samples. (c) LEfSe analysis of bacteria in different treatment groups (LDA > 3.5). Starting from the center and moving outwards, each circle represents a different taxonomic level: domain, phylum, class, order, family, and genus. Taxa that exhibit significant differences are highlighted with colors; Figure S4: The correlation between rhizosphere soil bacterial communities and environmental factors. (a) RDA analysis between genus-level bacterial communities and environmental parameters. The red clippings indicate environmental parameters, and the length and direction of the arrows represent the degree of influence of soil physicochemical properties on the structure of bacterial communities and the positive and negative correlations between them. (b) Spearman’s correlation analysis was conducted to examine the relationship between major bacterial communities and environmental parameters at the phylum level. (Statistical significance was denoted by *, **, and *** for p-values less than 0.05, 0.01, and 0.001, respectively); Figure S5: The functional composition of KEGG at primary (a) and secondary (b) levels. (c) The abundance bubble map of the carbon cycle pathway; Figure S6: The abundance composition of genus level rhizosphere bacteria in (a) carbon fixation, (b) central metabolism, and (c) methane metabolism; Table S1: The chemical and physical properties of soil sample (0 d incubation).

Author Contributions

Conceptualization and investigation, X.C. and X.L.; writing—original draft preparation, X.C.; statistical analysis and figure drawing, X.C. and Z.C.; reading and revision of the manuscript, funding acquisition, and project administration, Z.W. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the NSFC [grant number 41976126]; the Shenzhen Science and Technology Program [grant numbers RCJC20200714114433069, KCXFZ20201221173211033, ZDSYS20230626091459009, and JCYJ20230807111759016], as well as the Innovation Team Project for Guangdong’s Universities [grant number 2023KCXTD052].

Data Availability Statement

Sequence data used in this study were submitted to the Sequence Read Archive of the National Center for Biotechnology Information (accession number: 16S rRNA: RRJNA1039250; Macrogenome: PRJNA1080489).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bhardwaj, D.; Ansari, M.W.; Sahoo, R.K.; Tuteja, N. Biofertilizers function as key player in sustainable agriculture by improving soil fertility. plant tolerance and crop productivity. Microb. Cell Fact. 2014, 13, 66. [Google Scholar] [CrossRef] [PubMed]
  2. Mahanty, T.; Bhattacharjee, S.; Goswami, M.; Bhattacharyya, P.; Das, B.; Ghosh, A.; Tribedi, P. Biofertilizers: A potential approach for sustainable agriculture development. Environ. Sci. Pollut. Res. 2016, 24, 3315–3335. [Google Scholar] [CrossRef] [PubMed]
  3. Nosheen, S.; Ajmal, I.; Song, Y.D. Microbes as biofertilizers, a potential approach for sustainable crop production. Sustainability 2021, 13, 1868. [Google Scholar] [CrossRef]
  4. Singh, J.S.; Pandey, V.C.; Singh, D.P. Efficient soil microorganisms: A new dimension for sustainable agriculture and environmental development. Agric. Ecosyst. Environ. 2011, 140, 339–353. [Google Scholar] [CrossRef]
  5. Ibáñez, A.; Diez-Galán, A.; Cobos, R.; Calvo-Peña, C.; Barreiro, C.; Medina-Turienzo, J.; Sánchez-García, M.; Coque, J.J.R. Using Rhizosphere Phosphate Solubilizing Bacteria to Improve Barley (Hordeum vulgare) Plant Productivity. Microorganisms. 2021, 9, 1619. [Google Scholar] [CrossRef]
  6. Boots, B.; Russell, C.W.; Green, D.S. Effects of microplastics in soil ecosystems: Above and below ground. Environ. Sci. Technol. 2019, 53, 11496–11506. [Google Scholar] [CrossRef]
  7. Fuller, S.; Gautam, A. A procedure for measuring microplastics using pressurized fluid extraction. Environ. Sci. Technol. 2016, 50, 5774–5780. [Google Scholar] [CrossRef]
  8. Luca, N.; Martyn, F.; Sindre, L. Are agricultural soils dumps for microplastics of urban origin? Environ. Sci. Technol. 2016, 50, 10777–10779. [Google Scholar]
  9. Montazer, Z.; Habibi Najafi, M.B.; Levin, D.B. Challenges with verifying microbial degradation of polyethylene. Polymers 2020, 12, 123. [Google Scholar] [CrossRef]
  10. Silvia, C.; Riccardo, F.; Majid, G.; Jonan, M.A.; Stefano, L. Exploring sustainable alternatives: Wood distillate alleviates the impact of bioplastic in basil plants. Sci. Total Environ. 2023, 900, 166484. [Google Scholar]
  11. Azarnejad, N.; Celletti, S.; Ghorbani, M.; Fedeli, R.; Loppi, S. Dose-Dependent Effects of a Corn Starch-Based Bioplastic on Basil (Ocimum basilicum L.): Implications for Growth, Biochemical Parameters, and Nutrient Content. Toxics 2024, 12, 80. [Google Scholar] [CrossRef] [PubMed]
  12. Fan, P.; Yu, H.; Xi, B.D.; Tan, W.B. A review on the occurrence and influence of biodegradable microplastics in soil ecosystems: Are biodegradable plastics substitute or threat? Environ. Int. 2022, 163, 107244. [Google Scholar] [CrossRef] [PubMed]
  13. Qin, M.; Chen, C.Y.; Song, B.; Shen, M.C.; Cao, W.C.; Yang, H.L.; Zeng, G.M.; Gong, J.L. A review of biodegradable plastics to biodegradable microplastics: Another ecological threat to soil environments? J. Cleaner Prod. 2021, 312, 127816. [Google Scholar] [CrossRef]
  14. Rodrigo, M.; Paolina, G.; Jos, M. Raaijmakers The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 2013, 37, 634–663. [Google Scholar]
  15. Dilfuza, E.; Faina, K.; Shamil, V.; Laziza, G.; Zulfiya, K.; Ben, L. High incidence of plant growth-stimulating bacteria associated with the rhizosphere of wheat grown on salinated soil in Uzbekistan. Environ. Microbiol. 2007, 10, 1–9. [Google Scholar]
  16. Amanda, A.G.; Paula, M.D.J.; Rafaela, S.A.N.; Ligiane, A.F.; Karina, B.S.; Fatima, M.S.M. Genetic and symbiotic diversity of nitrogen-fixing bacteria isolated from agricultural soils in the western amazon by using cowpea as the trap plant. Appl. Environ. Microbiol. 2012, 78, 6726–6733. [Google Scholar]
  17. Petra, M.; David, C.; Zed, R. Rhizosphere interactions between microorganisms and plants govern iron and phosphorus acquisition along the root axis-model and research methods. Soil Biol. Biochem. 2011, 43, 883–894. [Google Scholar]
  18. Collignon, C.; Uroz, S.; Turpault, M.P.; Frey-Klett, P. Seasons differently impact the structure of mineral weathering bacterial communities in beech and spruce stands. Soil Biol. Biochem. 2011, 43, 2012–2022. [Google Scholar] [CrossRef]
  19. Christos, Z.; Corné, M.J.P. Modulation of host immunity by beneficial microbes. Mol. Plant-Microbe Interact. 2012, 25, 139–150. [Google Scholar]
  20. Li, C.D.; Cui, Q.; Li, Y.; Zhang, K.; Lu, X.; Zhang, Y. Effect of LDPE and biodegradable PBAT primary microplastics on bacterial community after four months of soil incubation. J. Hazard. Mater. 2022, 429, 128353. [Google Scholar] [CrossRef]
  21. Michaela, K.R.; Lucy, M.G.; Martine, G.; Charlotte, E.M.L.; Richard, P.E.; Dave, R.C.; Davey, L.J. LDPE and biodegradable PLA-PBAT plastics differentially affect plant-soil nitrogen partitioning and dynamics in a Hordeum vulgare mesocosm. J. Hazard. Mater. 2023, 447, 130825. [Google Scholar]
  22. Wang, T.; Cheng, K.K.; Huo, X.J.; Meng, P.P.; Cai, Z.H.; Wang, Z.K.; Zhou, J. Bioorganic fertilizer promotes pakchoi growth and shapes the soil microbial structure. Front. Plant Sci. 2022, 13, 1040437. [Google Scholar] [CrossRef]
  23. Sreejata, B.; Lluis, M.C.; Ana, M.P.; Jennifer, M.D. Biodegradable plastic mulch films: Impacts on soil microbial communities and ecosystem functions. Front. Microbiol. 2018, 9, 819. [Google Scholar]
  24. Dincă, L.C.; Grenni, P.; Onet, C.; Onet, A. Fertilization and soil microbial community: A review. Appl. Sci. 2022, 12, 1198. [Google Scholar] [CrossRef]
  25. De Souza Machado, A.A.; Lau, C.W.; Till, J.; Kloas, W.; Lehmann, A.; Becker, R.; Rillig, M.C. Impacts of microplastics on the soil biophysical environment. Environ. Sci. Technol. 2018, 52, 9656–9665. [Google Scholar] [CrossRef] [PubMed]
  26. Liu, H.F.; Yang, X.M.; Liu, G.B.; Liang, C.T.; Xue, S.; Chen, H.; Ritsema, C.J.; Geissen, V. Response of soil dissolved organic matter to microplastic addition in Chinese loess soil. Chemosphere 2017, 185, 907–917. [Google Scholar] [CrossRef]
  27. Xiao, H.; Liu, G.; Liu, P.L.; Zheng, F.L.; Zhang, J.Q.; Hu, F.N. Sediment transport capacity of concentrated flows on steep loessial slope with erodible beds. Sci. Rep. 2017, 7, 2350. [Google Scholar] [CrossRef]
  28. He, K.; Li, X.P.; Dong, L.L. The effects of flue gas desulfurization gypsum (FGD gypsum) on P fractions in a coastal plain soil. J. Soils Sediments. 2017, 18, 804–815. [Google Scholar] [CrossRef]
  29. Jiang, W.T.; Wang, X.K.; Xing, Y.Y.; Liu, X.H.; Cui, Z.G.; Yang, L.H. Enhancing rice production by potassium management: Recommended reasonable fertilization strategies in different inherent soil productivity levels for a sustainable rice production system. Sustainability 2019, 11, 6522. [Google Scholar] [CrossRef]
  30. Zhou, J.; Jia, R.; Brown, R.W.; Yang, Y.D.; Zeng, Z.H.; Jones, D.L.; Zang, H.D. The long-term uncertainty of biodegradable mulch film residues and associated microplastics pollution on plant-soil health. J. Hazard. Mater. 2023, 44, 130055. [Google Scholar] [CrossRef]
  31. Liu, D.; Iqbal, S.; Gui, H.; Xu, J.C.; An, S.; Xing, B.S. Nano-iron oxide (Fe3O4) mitigates the effects of microplastics on a ryegrass soil-microbe-plant system. ACS Nano. 2024, 17, 24867–24882. [Google Scholar] [CrossRef] [PubMed]
  32. Zhu, D.; Ma, J.; Li, G.; Rillig, M.C.; Zhu, Y.G. Soil plastispheres as hotpots of antibiotic resistance genes and potential pathogens. ISME J. 2021, 16, 521–532. [Google Scholar] [CrossRef] [PubMed]
  33. Hahladakis, J.N.; Velis, C.A.; Weber, R.; Iacovidou, E.; Purnell, P. An overview of chemical additives present in plastics: Migration, release, fate and environmental impact during their use, disposal and recycling. J. Hazard. Mater. 2018, 344, 179–199. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, M.J.; Zhao, Y.R.; Qin, X.; Jia, W.Q.; Chai, L.W.; Huang, M.K.; Huang, Y. Microplastics from mulching film is a distinct habitat for bacteria in farmland soil. Sci. Total Environ. 2019, 688, 470–478. [Google Scholar] [CrossRef] [PubMed]
  35. Abraham, J.; Ghosh, E.; Mukherjee, P.; Gajendiran, A. Microbial degradation of low density polyethylene. Environ. Prog. Sustainable Energy 2016, 36, 147–154. [Google Scholar] [CrossRef]
  36. Muhonja, C.N.; Magoma, G.; Imbuga, M.; Makonde, H.M. Molecular characterization of low-density polyethene (LDPE) degrading bacteria and fungi from Dandora Dumpsite, Nairobi, Kenya. Int. J. Microbiol. 2018, 2018, 4167845. [Google Scholar]
  37. Fei, Y.F.; Huang, S.Y.; Zhang, H.B.; Tong, Y.Z.; Wen, D.S.; Xia, X.Y.; Wang, H.; Luo, Y.M.; Barceló, D. Response of soil enzyme activities and bacterial communities to the accumulation of microplastics in an acid cropped soil. Sci. Total Environ. 2019, 707, 135634. [Google Scholar] [CrossRef]
  38. Ali, N.; Khan, S.; Li, Y.Y.; Zheng, N.G.; Yao, H.Y. Influence of biochars on the accessibility of organochlorine pesticides and microbial community in contaminated soils. Sci. Total Environ. 2019, 647, 551–560. [Google Scholar] [CrossRef]
  39. Awasthi, M.K.; Chen, H.Y.; Wang, Q.; Liu, T.; Duan, Y.M.; Awasthi, S.K.; Ren, X.N.; Tu, Z.N.; Li, J.; Zhao, J.; et al. Succession of bacteria diversity in the poultry manure composted mixed with clay: Studies upon its dynamics and associations with physicochemical and gaseous parameters. Bioresour. Technol. 2018, 267, 618–625. [Google Scholar] [CrossRef]
  40. Liang, J.; Tang, S.Q.; Gong, J.L.; Zeng, G.M.; Tang, W.W.; Song, B.; Zhang, P.; Yang, Z.X.; Luo, Y. Responses of enzymatic activity and microbial communities to biochar/compost amendment in sulfamethoxazole polluted wetland soil. J. Hazard. Mater. 2019, 385, 121533. [Google Scholar] [CrossRef]
  41. Xu, N.; Tan, G.C.; Wang, H.Y.; Gai, X.P. Effect of biochar additions to soil on nitrogen leaching, microbial biomass and bacterial community structure. Eur. J. Soil Biol. 2016, 74, 1–8. [Google Scholar] [CrossRef]
  42. Jaiswal, A.K.; Elad, Y.; Paudel, I.; Graber, E.R.; Cytryn, E.; Frenkel, O. Linking the Belowground Microbial Composition, Diversity and Activity to Soilborne Disease Suppression and Growth Promotion of Tomato Amended with Biochar. Sci. Rep. 2017, 7, 44382. [Google Scholar] [CrossRef] [PubMed]
  43. Ma, Q.; Qu, Y.Y.; Shen, W.L.; Zhang, Z.J.; Wang, J.W.; Liu, Z.Y.; Li, D.X.; Li, H.J.; Zhou, J.T. Bacterial community compositions of coking wastewater treatment plants in steel industry revealed by Illumina high-throughput sequencing. Bioresour. Technol. 2015, 179, 436–443. [Google Scholar] [CrossRef] [PubMed]
  44. Dammak, M.; Fourati, Y.; Tarrés, Q.; Delgado-Aguilar, M.; Mutjé, P.; Boufi, S. Blends of PBAT with plasticized starch for packaging applications: Mechanical properties, rheological behaviour and biodegradability. Ind. Crops Prod. 2020, 144, 112061. [Google Scholar] [CrossRef]
  45. She, J.Y.; Liu, J.; He, H.P.; Zhang, Q.; Lin, Y.Y.; Wang, J.; Yin, M.L.; Wang, L.L.; Wei, X.D.; Huang, Y.L.; et al. Microbial response and adaption to thallium contamination in soil profiles. J. Hazard. Mater. 2020, 423, 127080. [Google Scholar] [CrossRef]
  46. De Vries, F.T.; Griffiths, R.I.; Bailey, M.; Craig, H.; Girlanda, M.; Gweon, H.S.; Hallin, S.; Kaisermann, A.; Keith, A.; Kretzschmar, M.; et al. Soil bacterial networks are less stable under drought than fungal networks. Nat. Commun. 2018, 9, 3033. [Google Scholar] [CrossRef]
  47. Li, P.F.; Liu, J.; Saleem, M.; Li, G.L.; Luan, L.; Wu, M.; Li, Z.P. Reduced chemodiversity suppresses rhizosphere microbiome functioning in the mono-cropped agroecosystems. Microbiome. 2022, 10, 108. [Google Scholar] [CrossRef]
  48. Xiong, X.; Liao, H.; Xing, Y.F.; Han, X.K.; Wang, W.L.; Wan, W.J.; Huang, Q.Y.; Chen, W.L. Soil aggregates and fertilizer treatments drive bacterial interactions via interspecies niche overlap. Microbiol. Spectrum. 2022, 10, e02524-21. [Google Scholar] [CrossRef]
  49. Wu, C.C.; Ma, Y.J.; Wang, D.; Shan, Y.P.; Song, X.P.; Hu, H.Y.; Ren, X.L.; Ma, X.Y.; Cui, J.J.; Ma, Y. Integrated microbiology and metabolomics analysis reveal plastic mulch film residue affects soil microorganisms and their metabolic functions. J. Hazard. Mater. 2022, 423, 127258. [Google Scholar] [CrossRef]
  50. Ma, R.J.; Xu, Z.N.; Sun, J.Y.; Li, D.R.; Cheng, Z.; Niu, Y.L.; Guo, H.; Zhou, J.; Wang, T.C. Microplastics affect C, N, and P cycling in natural environments: Highlighting the driver of soil hydraulic properties. J. Hazard. Mater. 2023, 459, 132326. [Google Scholar] [CrossRef]
  51. Mai, Z.M.; Ye, M.; Wang, Y.S.; Foong, S.W.; Wang, L.; Sun, F.L.; Cheng, H. Characteristics of microbial community and function with the succession of mangroves. Front. Microbiol. 2021, 12, 764974. [Google Scholar] [CrossRef]
  52. Choi, G.M.; Im, W.T. Paraburkholderia azotifigens sp. nov., a nitrogen-fixing bacterium isolated from paddy soil. Int. J. Syst. Evol. Microbiol. 2018, 68, 310–316. [Google Scholar] [CrossRef] [PubMed]
  53. Yang, J.; He, J.H.; Jia, L.; Gu, H.Y. Integrating metagenomics and metabolomics to study the response of microbiota in black soil degradation. Sci. Total Environ. 2023, 899, 165486. [Google Scholar] [CrossRef] [PubMed]
  54. Fu, Y.S.; Yan, R.; Liu, D.L.; Jiang, S.W.; Cui, L.; Guo, X.W.; Wang, X.J.; Zhang, J.; Xiang, W.S. Trinickia diaoshuihuensis sp. nov., a plant growth promoting bacterium isolated from soil. Int. J. Syst. Evol. Microbiol. 2019, 69, 291–296. [Google Scholar] [CrossRef] [PubMed]
  55. Benitez, M.S.; Gardener, B.B.M. Linking Sequence to Function in Soil Bacteria: Sequence-Directed Isolation of Novel Bacteria Contributing to Soilborne Plant Disease Suppression. Appl. Environ. Microbiol. 2008, 75, 915–924. [Google Scholar] [CrossRef]
  56. Rodgers-Vieira, E.A.; Zhang, Z.F.; Adrion, A.C.; Gold, A.; Aitken, M.D. Identification of Anthraquinone-Degrading Bacteria in Soil Contaminated with Polycyclic Aromatic Hydrocarbons. Appl. Environ. Microbiol. 2015, 81, 3775–3781. [Google Scholar] [CrossRef]
Figure 1. The effect of MPs on soil property existence with biofertilizer: (a) SOC; (b) OM; (c) pH. Data are presented as the mean ± SD (n = 6). Different letters (a, b, and c) were statistically different (p < 0.05) according to Duncan’s new multivariate range test.
Figure 1. The effect of MPs on soil property existence with biofertilizer: (a) SOC; (b) OM; (c) pH. Data are presented as the mean ± SD (n = 6). Different letters (a, b, and c) were statistically different (p < 0.05) according to Duncan’s new multivariate range test.
Plants 13 02627 g001
Figure 2. The effects of biofertilizers and MPs on rhizosphere bacterial community structure and composition. (a) ACE and (b) Shannon index. Statistical significance was indicated by *, **, and *** for p < 0.05, p < 0.01, and p < 0.001, respectively. (c) Non-metric multidimensional scaling (NMDS) plots with corresponding pressure values for different bacterial groups at the OTU level. (d) Relative abundance maps of bacterial distribution at the phylum levels.
Figure 2. The effects of biofertilizers and MPs on rhizosphere bacterial community structure and composition. (a) ACE and (b) Shannon index. Statistical significance was indicated by *, **, and *** for p < 0.05, p < 0.01, and p < 0.001, respectively. (c) Non-metric multidimensional scaling (NMDS) plots with corresponding pressure values for different bacterial groups at the OTU level. (d) Relative abundance maps of bacterial distribution at the phylum levels.
Plants 13 02627 g002
Figure 3. The effect of biofertilizers and MPs on the network of rhizosphere bacterial communities at the genus level. Nodes are colored according to module, with pink edges indicating a positive correlation and blue indicating a negative correlation. Connections are defined as strongly correlated (Spearman’s |ρ| > 0.6) and significantly correlated (p < 0.05). The topological parameters of microbial networks obtained in all treatments are shown in the table.
Figure 3. The effect of biofertilizers and MPs on the network of rhizosphere bacterial communities at the genus level. Nodes are colored according to module, with pink edges indicating a positive correlation and blue indicating a negative correlation. Connections are defined as strongly correlated (Spearman’s |ρ| > 0.6) and significantly correlated (p < 0.05). The topological parameters of microbial networks obtained in all treatments are shown in the table.
Plants 13 02627 g003
Figure 4. The effects of biofertilizer and MPs on the function and pathways related to carbon cycling in rhizosphere bacteria. (a) A heatmap of KEGG secondary functional clustering for the top 20 abundances in different treatment groups. Hierarchical clustering (Bray–Curtis) analysis of all samples based on the similarity of KEGG pathways; (b) LEfSe analysis of different treatment groups under KEGG tertiary levels of functional annotation. (c) A heat map of carbon cycle genes, and (d) the relative abundance of the carbon cycle pathway. The pie chart shows the relative abundance of each pathway for each sample clock, and the size of the pie chart shows the total relative abundance of that pathway. The * sign in the upper-right corner of the pie chart indicates the abundance of the pathway differs significantly between different groups of samples. The * and ** indicate significance at p < 0.05 and p < 0.01.
Figure 4. The effects of biofertilizer and MPs on the function and pathways related to carbon cycling in rhizosphere bacteria. (a) A heatmap of KEGG secondary functional clustering for the top 20 abundances in different treatment groups. Hierarchical clustering (Bray–Curtis) analysis of all samples based on the similarity of KEGG pathways; (b) LEfSe analysis of different treatment groups under KEGG tertiary levels of functional annotation. (c) A heat map of carbon cycle genes, and (d) the relative abundance of the carbon cycle pathway. The pie chart shows the relative abundance of each pathway for each sample clock, and the size of the pie chart shows the total relative abundance of that pathway. The * sign in the upper-right corner of the pie chart indicates the abundance of the pathway differs significantly between different groups of samples. The * and ** indicate significance at p < 0.05 and p < 0.01.
Plants 13 02627 g004
Figure 5. A partial Mantel test between functional genes and environmental factors under biofertilizer and MP exposure. The statistical significance based on 9999 permutations is indicated by the edge color, while the Mantel r statistic correlation for the corresponding distance is represented by the edge width. The Spearman correlation coefficient is shown by the color gradient. *** denotes that the significance level for the Spearman correlation coefficient is 1.
Figure 5. A partial Mantel test between functional genes and environmental factors under biofertilizer and MP exposure. The statistical significance based on 9999 permutations is indicated by the edge color, while the Mantel r statistic correlation for the corresponding distance is represented by the edge width. The Spearman correlation coefficient is shown by the color gradient. *** denotes that the significance level for the Spearman correlation coefficient is 1.
Plants 13 02627 g005
Figure 6. PLS-PM analysis showed the relationship between MP type, soil physicochemical properties, bacterial diversity, and rhizosphere carbon cycling (goodness of fit = 0.713). The width of the arrow indicates the strength of the path factor. The red and blue lines represent positive and negative effects, respectively. Solid and dashed lines indicate direct and indirect effects, respectively. In the PLS-PM model, the number on the line represents the total effect value (statistical significance was denoted by ** and *** for p-values less than 0.01 and 0.001).
Figure 6. PLS-PM analysis showed the relationship between MP type, soil physicochemical properties, bacterial diversity, and rhizosphere carbon cycling (goodness of fit = 0.713). The width of the arrow indicates the strength of the path factor. The red and blue lines represent positive and negative effects, respectively. Solid and dashed lines indicate direct and indirect effects, respectively. In the PLS-PM model, the number on the line represents the total effect value (statistical significance was denoted by ** and *** for p-values less than 0.01 and 0.001).
Plants 13 02627 g006
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

Cheng, X.; Li, X.; Cai, Z.; Wang, Z.; Zhou, J. The Structural and Functional Responses of Rhizosphere Bacteria to Biodegradable Microplastics in the Presence of Biofertilizers. Plants 2024, 13, 2627. https://doi.org/10.3390/plants13182627

AMA Style

Cheng X, Li X, Cai Z, Wang Z, Zhou J. The Structural and Functional Responses of Rhizosphere Bacteria to Biodegradable Microplastics in the Presence of Biofertilizers. Plants. 2024; 13(18):2627. https://doi.org/10.3390/plants13182627

Chicago/Turabian Style

Cheng, Xueyu, Xinyang Li, Zhonghua Cai, Zongkang Wang, and Jin Zhou. 2024. "The Structural and Functional Responses of Rhizosphere Bacteria to Biodegradable Microplastics in the Presence of Biofertilizers" Plants 13, no. 18: 2627. https://doi.org/10.3390/plants13182627

APA Style

Cheng, X., Li, X., Cai, Z., Wang, Z., & Zhou, J. (2024). The Structural and Functional Responses of Rhizosphere Bacteria to Biodegradable Microplastics in the Presence of Biofertilizers. Plants, 13(18), 2627. https://doi.org/10.3390/plants13182627

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