- 1Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing, China
- 2Jinhejianong Biotechnology Co., Ltd., Beijing, China
Microbial management offers a sustainable pathway to enhance crop performance by optimizing plant-associated microbiomes. However, integrated strategies that concurrently target both the rhizosphere and phyllosphere to improve fruit tree productivity and quality remain underexplored. This study systematically evaluated the effects of combined soil and foliar microbial applications on the yield, fruit quality, and microbiome dynamics of ‘Yuluxiang’ pear. We compared conventional fertilization (CK) with two treatments: CK plus a soil-applied anti-replant disease agent (CF) and CK plus both the soil agent and a foliar growth-promoting inoculant (CFW). Microbial applications significantly increased the yield by up to 60.4% in CF treatment, and enhanced key fruit quality parameters, including soluble solids content (increased by 17.2% in CF and 16.7% in CFW) and fruit shape index. These agronomic improvements were closely associated with a targeted restructuring of bacterial communities in both the phyllosphere and rhizosphere. Specifically, beneficial phyla such as Actinomycetota were enriched in the phyllosphere under CFW treatment, while Bacillota increased in the rhizosphere under microbial amendments. Furthermore, co-occurrence network analysis revealed that microbial applications fostered more complex and cooperative microbial networks, with increased nodes and edges across both compartments. This work demonstrates that an integrated soil and foliar microbiome management strategy can mitigate replant disease constraints and elevate fruit quality, providing a practical approach for sustainable orchard production.
1 Introduction
Pears, as a key economic fruit within the Rosaceae family, are cultivated on a global scale (Kim et al., 2024). The ‘Yuluxiang’ pear, a prominent cultivar in Chinese orchards, is highly appreciated for its desirable sensory qualities, including a crisp bite, abundant juiciness, and high sweetness (Lozano et al., 2023). However, the reliable expression of these superior traits remains a challenge (Kumar et al., 2019). In many production areas, common issues include undersized fruit, poor skin coloration, and low soluble sugar content, which collectively diminish fruit marketability and grower profitability (Abeyrathna et al., 2023). Refining nutrient management practices is therefore essential to mitigate these constraints and achieve consistent fruit quality.
Conventional chemical fertilization, the traditional orchard management practice, focuses on supplementing macronutrients but fails to address the underlying microbial imbalance and pathogen accumulation in replanted soils. This long-term reliance can exacerbate soil degradation and negatively alter the structure and function of soil microbial communities (Jin et al., 2022). Meanwhile, the sustainable production is critically constrained by replant disease (RD), a pervasive problem in major fruit-producing regions. This issue is exemplified by apple replant disease (ARD), a specific manifestation of RD primarily driven by pathogen accumulation in the rhizosphere, notably fungi such as Fusarium, Rhizoctonia, and Cylindrocarpon species, and oomycetes like Pythium and Phytophthora (Geng et al., 2022). A key aggravating factor is the root-secreted toxic accumulation of phenolic acids (e.g., cinnamic acid, phlorizin), which impairs the root antioxidant system and further promotes pathogen proliferation (Duan et al., 2025). These pathological and biochemical stresses collectively result in reduced tree vigor, significant yield losses (20–50%), and deteriorated fruit quality, including lower soluble solid content and firmness (Duan et al., 2022). Conventional control has long relied on broad-spectrum soil chemical fumigants like methyl bromide. However, due to their high toxicity, ozone-depleting potential, and adverse environmental impacts, the search for sustainable alternatives is imperative (Geng et al., 2022). Among these alternatives, soil-applied microbial agents that suppress replant disease offer a promising strategy to mitigate disease pressure and improve plant health (Durán et al., 2021). Furthermore, combining disease-suppressing with growth-promoting microbial agents may yield synergistic benefits, enhancing plant stress resistance and fruit quality through the targeted modulation of soil microbial communities (Mishu et al., 2025). Complementing soil-based strategies, foliar fertilization serves as a targeted nutritional supplement that enhances nutrient use efficiency and directly influences fruit development and metabolism. It may also secondarily alter the composition and interactions of the rhizosphere bacterial community (Demir et al., 2024). Given the vital role of microbial assemblages in nutrient cycling and plant health, this study investigates the effects of microbial agents on the yield, fruit quality, and associated bacterial community dynamics of ‘Yuluxiang’ pears. The objective is to establish a foundation for a more sustainable and effective integrated nutrient management system.
Biological control agents (BCAs), particularly those based on Bacillus spp., represent promising eco-friendly alternatives (Duan et al., 2025; Nzimande et al., 2025). Their efficacy stems from multiple mechanisms, including: (1) the synthesis of antimicrobial compounds (e.g., lipopeptides, siderophores, lytic enzymes) and volatile organic compounds (VOCs) for pathogen inhibition; (2) the degradation of phytotoxic phenolic acids; and (3) the modulation of soil microbial communities and key enzyme activities (e.g., urease, phosphatase) to enhance nutrient cycling (Chen et al., 2020; Duan et al., 2022). Field trials with specific strains, such as Bacillus amyloliquefaciens QSB-6, have validated these mechanisms, demonstrating significant plant growth promotion alongside reduced pathogen abundance and phenolic acid content in replanted soils (Duan et al., 2022). However, despite this demonstrated potential, current research and practical application remain predominantly focused on the single soil application of these disease-suppressing agents.
This soil-centric focus overlooks the potential of complementary aboveground strategies. The phyllosphere (leaf surface) constitutes a critical, yet distinct, microbial habitat integral to plant health and productivity (Leveau, 2019). The foliar application of plant growth-promoting bacteria (PGPB) can directly colonize the phyllosphere, enhancing photosynthetic efficiency, boosting antioxidant enzyme activities, and improving nutrient absorption (Hirata et al., 2025). This approach has proven effective in various crops. For instance, in sweet cherry, foliar-applied PGPB significantly increased yield, shoot growth, and leaf mineral content (e.g., N, P, K, Fe, Zn) (Esitken et al., 2006). Similarly, in apple, foliar application of Pseudomonas and Bacillus strains at key phenological stages enhanced plant growth and yield (Kordatzaki et al., 2022). These benefits are mediated through the diverse and complementary mechanisms of action employed by PGPB, which collectively support plant development and stress resilience (Acuña et al., 2024). These can be categorized as direct or indirect. Direct mechanisms enhance plant nutrition and growth, including biological nitrogen fixation, phytohormone production, and mineral solubilization (e.g., phosphorus) (Agake et al., 2022). Indirect mechanisms primarily involve biocontrol, such as synthesizing antibiotics, lytic enzymes, and antifungal metabolites to suppress phytopathogens, and modulating plant defense signaling pathways (Guardiola-Márquez et al., 2022). Collectively, these interactions enhance the plant’s ability to mitigate abiotic and biotic stresses, thereby supporting sustained health and productivity (Efthimiadou et al., 2020). Furthermore, the efficacy of such treatments can be optimized by using indigenous bacterial strains, which are often more competitive and effective due to their adaptation to local conditions (Katsenios et al., 2022).
Microbial co-occurrence networks are a powerful framework for elucidating complex interactions within soil microbiomes (Lin et al., 2024). By identifying statistically robust co-occurrence patterns among microbial taxa, these networks capture both the structural organization and functional potential of microbial communities (Feng et al., 2025). Evidence suggests that more complex and modular networks support enhanced ecosystem functions, including nutrient cycling, organic matter turnover, pathogen suppression, and ultimately, plant productivity (Zhang et al., 2025). For example, introducing rhizobia has been shown to increase fungal network connectivity in the soybean rhizosphere (Xu et al., 2020), while inoculation with Rhizopus irregularis can enhance bacterial network complexity in maize roots (Darriaut et al., 2025). Integrating studies of plant phyllosphere and rhizosphere microbiomes with symbiotic network analysis can provide mechanistic insights into microbial community assembly and interactions. However, integrated strategies combining soil-applied disease suppressants with foliar growth promoters remain largely unexplored. This represents both a significant knowledge gap and a promising avenue for developing integrated orchard management systems.
Therefore, to address critical knowledge gaps in the literature, including the lack of comparative application studies, unclear effects on soil and leaf microbial networks, and limited data linking microbial modulation to fruit quality, this study evaluated three management regimes: (1) conventional fertilization (CK); (2) conventional fertilization plus a soil-applied anti-replant disease agent (CF); and (3) conventional fertilization plus the soil-applied agent and a foliar growth-promoting microbial agent (CFW). By analyzing fruit quality, yield, and bacterial community dynamics in the rhizosphere and phyllosphere, this research aims to establish a scientific foundation for an integrated microbiome management strategy to mitigate replant disease and enhance the sustainable production of ‘Yuluxiang’ pears.
2 Materials and methods
2.1 Experimental site and design
The study was conducted in a high-density orchard located in Xi County, Shanxi Province, China (110°55′E, 36°41′N), a major apple-producing region. The site has a temperate continental monsoon climate and a loess hilly landform. The climate is characterized by hot, rainy summers (average temperature 21 °C; mean rainfall 303 mm) and cold, dry winters (average temperature 3.5 °C; mean rainfall 75.6 mm). The average annual sunshine duration is 2,740.9 h, and the frost-free period is 150–160 days.
The orchard was established in 2015 on flat terrain with deep, fertile soil characterized by moderately brown topsoil. The experimental tree species is the 11-year-old ‘Yuluxiang’ pear (Pyrus bretschneideri), planted in an east–west orientation with a spacing of 1.0 meters × 3.5 meters. The experiment included three treatments: conventional fertilization control (CK), conventional fertilization plus anti-replant disease agent (CF), and conventional fertilization plus anti-replant disease agent combined with plant growth-promoting microbial inoculant (CFW). A randomized block design was employed, with five replicates per treatment, establishing a total of 15 experimental plots.
The control treatment (CK) followed conventional local fertilization practices, which consisted of applying 2 tons of organic fertilizer per mu (approximately 0.067 hectares) in autumn via 30 cm deep strip trenches placed on one side of each tree. For the CF treatment, 4 kg per mu of a commercial anti-replant disease agent was incorporated into this standard fertilizer regimen. This agent is a microbial consortium primarily composed of Bacillus amyloliquefaciens, Bacillus subtilis, Trichoderma harzianum, Bacillus pumilus, Paenibacillus mucilaginosus, and Clonostachys spp., with a guaranteed minimum viable count of ≥1.0 × 109 CFU mL−1. The CFW treatment received the same soil application of the anti-replant disease agent, supplemented with a foliar spray of a growth-promoting microbial inoculant. The foliar inoculant, with a viable count of ≥1.0 × 109 CFU mL−1, contained a consortium of Bacillus amyloliquefaciens, Bacillus subtilis, and Trichoderma harzianum. The foliar application was performed twice: during the fruit expansion stage (mid-July) and the fruit maturation stage (early August). All fertilizers and microbial amendments were supplied by Jinhua Jianong (Beijing) Biotechnology Co., Ltd.
2.2 Sample collection
During the pear tree’s late expansion and ripening stages, leaf samples and rhizosphere soil samples were collected one week after spraying the microbial foliar agent. Five pear trees were randomly selected from each experimental plot. Ten healthy, fully expanded leaves were randomly collected from each tree, covering the upper, middle, and lower canopy layers and distributed across the four cardinal directions (east, west, south, north). Collected leaves were immediately placed in sterile plastic bags. Beneath each pear tree, roots were gently excavated with a shovel. Loose soil adhering to the roots was shaken off, while tightly adhering rhizosphere soil was brushed off and collected. This rhizosphere soil was sieved through a 2 mm mesh and placed into sterile bags. Analyze each experimental plot as an independent biological replicate, consolidating multiple subsamples collected from the same plot. Leaf and soil samples were stored at −80 °C for subsequent DNA isolation and characterization of leaf surface and rhizosphere bacterial communities.
2.3 Measurement of yield and quality indicators
In each plot, one pear tree is randomly selected. Using an electronic balance with an accuracy of 0.1 kilograms, the fruit yield of each pear tree is measured to determine the total yield of the entire tree. The yield per unit area was extrapolated based on the orchard planting density using the formula:
The planting density was calculated from the specific row (3.5 m) and plant (1.0 m) spacing.
From each cluster, randomly select 10 fruits and measure the maximum vertical height (from the base of the receptacle to the apex of the calyx) and the maximum equatorial diameter using a DL91150 electronic digital caliper. These parameters are recorded as longitudinal diameter and transverse diameter, respectively, with an accuracy of 1 millimeter. The fruit shape index was calculated as the ratio of longitudinal diameter to transverse diameter. Subsequently, pulp tissue was uniformly collected along the fruit’s equatorial region, pooled, and ground in a mortar. The ground mixture was filtered through four layers of sterile gauze to obtain clear juice. Total water-soluble solids content in the juice was determined using a handheld refractometer.
From each plot, 10 additional randomly selected fruits exhibiting uniform size, minimal mechanical damage, and no insect infestation were equilibrated at room temperature (20 ± 1 °C) for at least 2 h. A texture analyzer was then used to perform puncture tests on two relatively smooth areas of the fruit’s equatorial region, measuring fruit firmness.
2.4 DNA extraction, PCR, Illumina sequencing, and bioinformatics analysis
DNA was extracted from treated leaf samples and rhizosphere soil samples using the E. Z. N. A.® Plant DNA Kit (Omega Bio-tek, Norcross, GA, U. S.) and the FastDNA® Soil Spin Kit (MP Biomedicals, U. S.), respectively, according to the instructions. The integrity and purity of extracted DNA were assessed by 1% agarose gel electrophoresis, and concentration as well as purity were measured with a NanoDrop® One spectrophotometer (Thermo Fisher Scientific, USA). RCR amplification was performed using the extracted DNA as a template, employing specific bacterial primers 338F and 806R (F: 5′-ACTCCTACGGGAGGCAGCA-3′; R: 5′- GGACTACHVGGGTWTCTAAT-3′). Amplifications were carried out using Premix Taq polymerase (Takara Biotechnology, Dalian, China). PCR products were quantified with GeneTools analysis software (v4.03.05.0, SynGene). Based on concentration, amplicons were pooled in equimolar ratios and purified with the E. Z. N. A.® Gel Extraction Kit (Omega Bio-tek, Norcross, GA, USA). The target DNA fragments were eluted in TE buffer, and the final DNA concentration was determined using a NanoDrop 2000C spectrophotometer (Thermo Fisher Scientific). After equalizing the samples, a paired-end (PE) library was constructed using a NEXTFLEX Rapid DNA-Seq Kit. Equal amounts of purified amplicon were pooled and subjected to high-throughput sequencing (Illumina Nextseq 2000 platform PE150 sequencing mode; Majorbio Bio-pharm Technology, Shanghai, China).
Sequencing data underwent denoising using QIIME 2 (Caporaso Lab, Northern Arizona University, Flagstaff, AZ, USA) combined with the DADA2 plugin (Benjamin Callahan, North Carolina State University, Raleigh, NC, USA) on optimized sequences following quality control assembly. Sequences processed by DADA2 are typically referred to as ASVs (Amplifier Sequence Variants). Sequences annotated to chloroplast and mitochondrial origins were removed from all samples. To minimize the impact of sequencing depth on subsequent alpha and beta diversity analyses, the number of sequences per sample was downsampled to 20,000. After downsampling, the average sequence coverage (Good’s coverage) per sample remained at 99.09%. Species taxonomic analysis of ASVs was performed using the Naive Bayes classifier in Qiime2, based on the Sliva 16S rRNA gene database (v 138.2). In this study, the microbial raw sequences were deposited in the SRA database short-read archive SRR36362282.
2.5 Statistical analysis
All statistical analyses were performed using IBM SPSS Statistics (v27.1) and R (v4.3.1). All statistical analyses were conducted with experimental plots as independent biological replicates (n = 5 per treatment). Microbial sequencing data were derived from pooled samples across all plots. For agronomic data (yield and fruit quality parameters), treatment effects were assessed by one-way analysis of variance (ANOVA). When ANOVA indicated a significant overall effect (p < 0.05), means were compared using Tukey’s honestly significant difference (HSD) post-hoc test. General data visualization was performed using GraphPad Prism (v9.5). Bacterial community analyses were conducted in R. Alpha diversity (Chao1 and Shannon) was analyzed via QIIME2, the ggplot2 package (3.5.1) in R software, and SPSS 27.1. Differences in bacterial community structure (beta-diversity) among treatments were tested using Permutational Multivariate Analysis of Variance (PERMANOVA) based on Bray–Curtis dissimilarities with 9,999 permutations, implemented in the vegan package (v2.6–4). The homogeneity of group dispersions was checked. Principal coordinates analysis (PCoA) based on Bray–Curtis distance was used to visualize community composition.
Microbial co-occurrence networks were constructed separately for each treatment and compartment. Pairwise Spearman correlations were calculated among operational taxonomic units (OTUs) with a prevalence > 10%. Only robust correlations with |r| > 0.7 and a p-value < 0.01 were retained. Networks were built and their topological properties (nodes, edges, average degree, modularity) were calculated using the igraph package (v1.5.1). Mantel tests based on Pearson’s correlation were performed between the Bray–Curtis distance matrix of bacterial communities and the Euclidean distance matrices of pear yield and quality traits, and the relative abundance of dominant bacterial phyla.
3 Results
3.1 Impact of soil and foliar microbial applications on yield and quality of pear
Microbial agent applications significantly enhanced the yield and quality of ‘Yuluxiang’ pears (Table 1). The CF treatment produced the highest yield (37,035 kg ha−1), which was significantly greater than both the CFW (33,660 kg ha−1) and CK (23,085 kg ha−1) treatments. This constitutes a 60.4% yield increase over the conventional fertilization (CK) treatment. Fruit quality was significantly influenced by the microbial treatments. Both microbial treatments significantly enhanced fruit sweetness. The soluble solids content (SSC) reached 15.67% in CF and 15.60% in CFW, representing increases of 17.2 and 16.7%, respectively, over the CK. The fruit shape index indicated that microbial applications resulted in a shape closer to spherical (1.0). The indices for CF (0.94) and CFW (0.96) were both significantly higher than that of CK (0.86). For firmness, the CFW treatment yielded the highest value (4.43), which was significantly greater than that of CF (4.33) and CK (4.10).
3.2 Dynamics of phyllosphere and rhizosphere bacterial communities
Alpha diversity of bacterial communities differed markedly between the phyllosphere and rhizosphere and was significantly affected by treatment (Table 2). In the phyllosphere, and most notably at the fruit maturation stage, treatment effects were pronounced. Regarding richness, the CFW treatment significantly reduced the Chao1 index compared to the CK and CF treatments (p < 0.05). Similarly, for community diversity, the Shannon index was lowest under the CFW treatment, showing a significant decrease relative to both CK and CF (p < 0.05). In contrast, the CF treatment alone significantly increased the Shannon index compared to the CK control (p < 0.05). In the rhizosphere, both microbial treatments (CF and CFW) resulted in moderate reductions in bacterial richness during the fruit expansion stage compared to CK (p < 0.05). However, during fruit maturation, the CF treatment substantially enhanced richness relative to CK (p < 0.05). In contrast to these shifts in richness, Shannon diversity indices remained relatively consistent across all treatments and both growth stages.
Table 2. Alpha diversity of phyllosphere and rhizosphere bacterial communities in ‘Yuluxiang’ pears under different treatments at two growth stages.
The application of microbial agents significantly altered bacterial community composition at the phylum level in both the phyllosphere and rhizosphere, with more pronounced shifts occurring in the phyllosphere (Figure 1). During the fruit expansion stage, these phyllosphere communities were restructured. The CF treatment dramatically reduced the abundance of Bacillota from 64.85% (CK) to 3.70%, while promoting Pseudomonadota, which increased from 4.21 to 50.62%. The relative abundance of Actinomycetota also increased from 29.50 to 43.96% under the CF treatment. In contrast, the CFW treatment predominantly enriched Actinomycetota (77.43%), markedly reducing both Bacillota (9.23%) and Pseudomonadota (7.19%) relative to CK. During the fruit maturation stage, Actinomycetota remained the dominant phylum across all treatments, reaching its highest relative abundance (88.19%) in the CFW treatment. Notably, the CF treatment sustained a higher proportion of Bacillota (12.13%) than both CK (4.22%) and CFW (4.80%).
Figure 1. Mean relative abundance of bacteria at the phylum level during the fruit expansion stage and maturation stage in the phyllosphere (a,b) and rhizosphere (c,d), respectively. CK, conventional fertilization, control; CF, conventional fertilization with the addition of an anti-replant disease agent; CFW, conventional fertilization with the addition of an anti-replant disease agent and foliar application of a growth-promoting microbial agent.
In contrast to the phyllosphere, the rhizosphere bacterial communities exhibited greater stability in response to the microbial treatments. During the fruit expansion stage, Pseudomonadota remained the dominant phylum across all treatments, with CK, CF, and CFW showing relative abundances of 28.26, 24.35, and 23.56%, respectively. Acidobacteriota and Actinomycetota abundances were stable, although Actinomycetota showed a slight increase from 14.51% (CK) to 16.30% (CFW). Similarly, Chloroflexota and Gemmatimonadota displayed minimal fluctuations among treatments. This pattern of stability persisted into the fruit maturation stage. Pseudomonadota maintained its dominance, and Actinomycetota abundance remained stable, while Acidobacteriota decreased by 7.2% in CFW compared to CK. Notably, Bacillota abundance increased significantly in both CF and CFW treatments compared to CK, representing an 89.1% increase. Principal coordinates analysis (PCoA) based on Bray-Curtis distances revealed distinct clustering of bacterial communities according to treatment (Figure 2). In the phyllosphere during the fruit maturation stage, communities from CFW treatments formed distinct clusters separate from the CK control (PERMANOVA, p < 0.05), with PC1 explaining 57.83% of the variance. In contrast, rhizosphere communities demonstrated a different response pattern. During the fruit expansion stage, communities from CFW treatments formed distinct clusters separate from the CK control (PERMANOVA, p < 0.05), with PC1 explaining 35.83% of the variance.
Figure 2. PCoA of bacterial communities based on Bray-Curtis distance in phyllosphere (a,b) and rhizosphere (c,d) under different treatments at fruit expansion stage and maturation stage. CK, conventional fertilization, control; CF, conventional fertilization with the addition of an anti-replant disease agent; CFW, conventional fertilization with the addition of an anti-replant disease agent and foliar application of a growth-promoting microbial agent.
3.3 Dynamic changes in the network structure of bacterial communities on leaf surfaces and in the rhizosphere
Fertilization regimes significantly influenced the topological structure of bacterial co-occurrence networks in both the rhizosphere and phyllosphere of ‘YuluXiang’ pear (Figure 3). Compared to the CK control, microbial amendments (CF and CFW) generally increased network complexity (i.e., the number of nodes and edges) across both sampling periods. During the fruit expansion stage, compared to CK, the CF treatment augmented the network by 28 nodes and 430 edges. The CFW treatment exerted a more substantial influence, increasing the network by 54 nodes and 1,554 edges relative to CK. During the fruit maturation stage, the effect of the CF treatment diminished relative to the expansion stage, with only 5 additional nodes and 344 edges. In contrast, CFW maintained a strong effect, increasing node and edge counts by 75 and 2,652, respectively, compared to CK.
Figure 3. Microbial networks in the phyllosphere (a–f) and rhizosphere (g–l) regions under different fertilization treatments during the fruit expansion stage (a–c, g–i) and maturation stage (d–f, j–l). CK represents conventional fertilization (a,d,g,j); CF denotes conventional fertilization with the addition of an anti-replant disease agent (b,e,h,k); CFW indicates conventional fertilization with the addition of an anti-replant disease agent and foliar application of a growth-promoting microbial agent (c,f,i,l).
3.4 Interrelationships between bacterial community composition and yield, quality, and dominant bacteria
Mantel tests assessed correlations of both pear yield/quality and dominant phyla abundance with the overall bacterial community structure. During the fruit expansion stage, a strong positive correlation was observed between overall phyllosphere community composition and the matrix of yield and quality traits. Separately, at the phylum level, Bacillota abundance showed a significant negative correlation with yield and quality (Figure 4a). During the fruit maturation, the correlation between overall community structure and yield/quality was not significant. At the phylum level, only Actinobacteria showed a weak positive correlation, whereas Pseudomonadota, Bacteroidota, and Myxococcota were significantly negatively correlated with yield and quality (Figure 4b). For the rhizosphere during the fruit expansion stage, the overall bacterial community structure showed no significant correlation with yield or quality metrics. In a separate analysis at the phylum level, the relative abundance of Pseudomonadota was significantly negatively correlated with yield and quality (Figure 4c). During the fruit maturation stage, a strong positive correlation was observed between the overall rhizosphere community composition and yield/quality. Furthermore, the composition was significantly correlated with the combined abundance of the major bacterial phyla (Figure 4d).
Figure 4. (a–d) Relationships between bacterial community composition (based on Bray–Curtis distance) and fruit phenotypic traits or dominant bacterial taxa in the phyllosphere [(a) fruit expansion stage; (b) maturation stage] and rhizosphere [(c) fruit expansion stage; (d) maturation stage]. The bacterial community composition was related to each fruit phenotypic trait or dominant bacterial taxon via a Mantel test. Line width corresponds to the Mantel’s r statistic (thick = strong correlation, medium = moderate correlation, thin = weak correlation), and line color denotes statistical significance based on permutations (blue = non-significant, yellow = p < 0.05, red = highly significant). Pairwise Pearson’s correlations among fruit phenotypic traits and dominant bacterial taxa are shown in the upper triangle, with a color gradient representing the correlation coefficient (solid line = positive correlation, dashed line = negative correlation).
4 Discussion
4.1 Effects of microbial agents on yield and fruit quality
The application of microbial agents significantly improved the yield and fruit quality of ‘Yuluxiang’ pears (Table 1). This aligns with established roles of plant growth-promoting bacteria (PGPB) in enhancing nutrient acquisition, stimulating plant growth, and modifying soil microbial community structure (Liu et al., 2025). The CF treatment yielded the highest production, likely because the soil-applied anti-replant agent suppressed soil-borne pathogens and improved nutrient availability, thereby promoting fruit set and development (Demir et al., 2024). This mechanism is supported by studies showing that microbial inoculants can induce beneficial shifts in resident microbial communities, thereby enhancing plant vigor and productivity (Joshi et al., 2019). Notably, both the CF and CFW treatments significantly increased the soluble solids content of the fruit, indicating enhanced sugar accumulation. This enhancement in the soluble solids content may be associated with the ability of PGPB to facilitate nutrient solubilization and modulate hormone-mediated growth processes, mechanisms documented in studies of plant growth regulation (Qin et al., 2021). Furthermore, microbial agents may influence photosynthate partitioning, thereby directing more carbohydrates to the fruit and increasing the soluble solids content (Lan et al., 2025). The improved fruit shape index under microbial treatments implies a role for PGPB in regulating fruit morphogenesis, likely mediated through hormonal signaling pathways (Zhang et al., 2022). Such effects have been documented in other horticultural crops, where microbial inoculants were found to enhance fruit size and uniformity by modulating endogenous hormone levels and improving photosynthetic efficiency (Wu et al., 2025). These observations underscore the potential of microbial agents as a sustainable tool for optimizing fruit appearance and yield.
4.2 Shifts in bacterial community structure in response to microbial treatments
Our results demonstrate that microbial applications significantly restructured bacterial community composition (Figure 1) and altered alpha diversity (Table 2) in both the phyllosphere and rhizosphere, with more pronounced effects in the phyllosphere. This is consistent with Qin et al. (2019), who reported that foliar application of BCAs could substantially modify the phyllosphere microbiota, leading to improved plant health. The reduction in phyllosphere bacterial richness under CFW treatment may reflect a selective enrichment of beneficial taxa such as Actinomycetota, which are known for their pathogen suppression and role in nutrient cycling (García et al., 2025). In the rhizosphere, microbial treatments significantly increased the abundance of Bacillota during the fruit maturation stage (Figure 1d). The relative stability of the rhizosphere community compared to the phyllosphere underscores its resilience but also highlights the potential for long-term modulation through repeated microbial inoculations (Díaz-Rodríguez et al., 2025). This is consistent with previous work showing that microbial inoculants can induce both transient and persistent changes in soil communities, depending on the inoculant and application strategy (Trabelsi and Mhamdi, 2013).
Beyond community composition, microbial applications altered the topological properties of bacterial co-occurrence networks (Figure 3). Both CF and CFW treatments increased the number of nodes and edges across sampling periods, indicating greater network complexity and enhanced microbial interactions. This is consistent with reports linking greater network complexity, particularly in the phyllosphere, to improved plant health and disease suppression (Hilber-Bodmer et al., 2017). The increased density and predominance of positive correlations in these networks suggest a shift toward stronger cooperative interactions and greater functional redundancy. This may enhance the community’s stability and resilience against pathogen invasion (Xing et al., 2024). Moreover, the enrichment of dominant taxa within these networks, particularly Bacillota and Actinomycetota, implies that microbial inoculants not only introduced beneficial strains but also reinforced the ecological roles of indigenous beneficial taxa (García et al., 2025). This principle extends to other beneficial taxa. For example, the genus Pseudomonas is well-known for promoting plant growth and suppressing pathogens via secondary metabolites, enzymes, and hormones (Soliman et al., 2023). The integration of such beneficial taxa into more complex and cooperative networks enhances functional redundancy and ecological stability. This restructured environment likely suppresses pathogens through intensified competition for resources and spatial niches (Cui et al., 2025). Principal coordinates analysis further confirmed that microbial treatments drove the formation of distinct bacterial clusters over time (Figure 2). The clear separation of communities in the PCoA plots indicates that inoculants induced significant and persistent shifts in microbiome structure, which aligns with the established capacity of microbial amendments to reshape resident communities (Berg et al., 2021). Collectively, these findings underscore the potential of tailored microbial applications to steer microbiome assembly.
4.3 Linking microbial community dynamics to yield and quality
Mantel test analyses revealed significant correlations between the structure of bacterial communities in both the phyllosphere and rhizosphere and the yield and quality of ‘Yuluxiang’ pears (Figure 4). The strong positive correlation between early-stage phyllosphere community structure and fruit yield/quality underscores the critical role of early microbial colonization in shaping development (Mehlferber et al., 2023). This early establishment of a beneficial microbiome may enhance the plant’s physiological state and direct resource allocation toward reproductive growth (Khoiri et al., 2021). This aligns with the understanding that a more complex phyllosphere network can suppress pathogens and improve plant health (Mehlferber et al., 2023). Such complexity fosters a competitive environment that reduces the available ecological niches for pathogens, a mechanism likely contributing to the improved outcomes observed here.
The relationship between the rhizosphere microbial community and plant performance also showed significant dynamics. The shift from a non-significant correlation at the fruit expansion stage to a strong positive correlation at maturation suggests that the influence of microbial inoculants on the rhizosphere microbiome intensifies during critical fruit development. This shift coincided with a significant increase in the abundance of Bacillota under the CF and CFW treatments (Figure 4d). Bacillota species (e.g., Bacillus spp.) are well-documented for their plant growth-promoting and biocontrol traits, including the production of antimicrobial compounds, lytic enzymes, and volatile organic compounds (Marco et al., 2022), as well as their ability to induce systemic resistance in plants (Balthazar et al., 2020). Their enrichment in the rhizosphere likely contributed to the observed improvements in fruit quality, particularly the increase in soluble solids and sugars, by enhancing nutrient availability and uptake (Lan et al., 2025). This aligns perfectly with the findings of, who reported that soil application of PGPB enhances plant growth and physiological performance through rhizosphere activity. This mechanism likely involves improved mineral nutrition and a more robust root system, thereby enhancing photoassimilate production and translocation to developing fruits (Sun and Shahrajabian, 2023). The distinct responses of phyllosphere and rhizosphere communities to microbial applications, evident in our PCoA and Mantel test results, indicate these compartments require targeted management to optimize plant outcomes. The phyllosphere, being an open and harsh environment, may require repeated applications or specific consortia to establish a stable, beneficial community (Hirata et al., 2025). In contrast, the rhizosphere community, while more buffered, exhibits a delayed but significant restructuring in response to inoculation, which has a profound impact on yield (Kui et al., 2021). This intricate interplay between the two compartments, mediated by the plant’s systemic physiology, ultimately determines the final economic output (Hacquard et al., 2022). Therefore, a precision management strategy that separately targets the phyllosphere and rhizosphere microbiomes at key growth stages presents a powerful approach to enhance fruit economic traits and advance sustainable orchard production.
5 Conclusion
This study shows that microbial agent application is an effective strategy for enhancing the yield and fruit quality of ‘Yuluxiang’ pears. Among the treatments, soil application of an anti-replant disease agent (CF) led to the greatest yield increase, while the combined soil and foliar microbial treatment (CFW) optimally enhanced fruit quality traits, including soluble solids content, shape, and firmness. These agronomic improvements were strongly associated with a targeted restructuring of the microbiome, characterized by the enrichment of beneficial taxa (e.g., Actinomycetota) and the development of more complex, cooperative microbial networks in both compartments. Notably, the response of the microbial community was compartment-specific: phyllosphere communities were rapidly reshaped, whereas rhizosphere communities exhibited a more gradual but ultimately significant restructuring correlated with final yield. Our results found that comprehensive microbial community management strategies targeting specific compartments can achieve sustainable pear tree production while enhancing fruit quality. Our findings demonstrate that an integrated, compartment-targeted microbiome management strategy can promote sustainable pear production: soil-applied microbial agents alleviate replant disease and enhance yield, while foliar-applied PGPB optimize fruit quality. Together, these applications provide a practical, science-based framework for building more resilient and productive pear cultivation systems.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, SRR36362282.
Author contributions
NL: Writing – original draft, Formal analysis, Data curation. YZ: Formal analysis, Data curation, Writing – original draft. ZR: Writing – review & editing, Methodology. XW: Writing – review & editing, Investigation, Data curation. HL: Formal analysis, Writing – review & editing, Supervision. AZ: Funding acquisition, Writing – review & editing, Supervision.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National Key Research and Development Program of China (2023YFD1900501).
Conflict of interest
ZR and XW were employed by Jinhejianong Biotechnology Co., Ltd.
The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Abeyrathna, R. M. R. D., Nakaguchi, V. M., Minn, A., and Ahamed, T. (2023). Recognition and counting of apples in a dynamic state using a 3D camera and deep learning algorithms for robotic harvesting systems. Sensors 23:3810. doi: 10.3390/s23083810,
Acuña, J. J., Rilling, J. I., Inostroza, N. G., Zhang, Q., Wick, L. Y., Sessitsch, A., et al. (2024). Variovorax sp. strain P1R9 applied individually or as part of bacterial consortia enhances wheat germination under salt stress conditions. Sci. Rep. 14:2070. doi: 10.1038/s41598-024-52535-0,
Agake, S., Plucani do Amaral, F., Yamada, T., Sekimoto, H., Stacey, G., Yokoyama, T., et al. (2022). Plant growth-promoting effects of viable and dead spores of Bacillus pumilus TUAT1 on Setaria viridis. Microbes Environ. 37:ME21060. doi: 10.1264/jsme2.ME21060,
Balthazar, C., Cantin, G., Novinscak, A., Joly, D. L., and Filion, M. (2020). Expression of putative defense responses in Cannabis primed by Pseudomonas and/or Bacillus strains and infected by Botrytis cinerea. Front. Plant Sci. 11:572112. doi: 10.3389/fpls.2020.572112,
Berg, G., Kusstatscher, P., Abdelfattah, A., Cernava, T., and Smalla, K. (2021). Microbiome modulation—toward a better understanding of plant microbiome response to microbial inoculants. Front. Microbiol. 12:650610. doi: 10.3389/fmicb.2021.650610,
Chen, K., Tian, Z., He, H., Long, C., and Jiang, F. (2020). Bacillus species as potential biocontrol agents against citrus diseases. Biol. Control 151:104419. doi: 10.1016/j.biocontrol.2020.104419
Cui, X., Xu, D., Zhang, Y., Huang, S., Wei, W., Ma, G., et al. (2025). Straw-enhanced soil bacterial robustness via resource-driven niche dynamics in tea plantations, South Henan, China. Microorganisms 13:832. doi: 10.3390/microorganisms13040832,
Darriaut, R., Lailheugue, V., Wastin, J., Tran, J., Martins, G., Ballestra, P., et al. (2025). Rhizobacteria from vineyard and commercial arbuscular mycorrhizal fungi induce synergistic microbiome shifts within grapevine root systems. Sci. Rep. 15:27884. doi: 10.1038/s41598-025-12673-5,
Demir, H., Saka, A. K., Uçan, U., Akgün, İ. H., and Yalçı, H. K. (2024). Impact of effective micro-organisms (EM) on the yield, growth and bio-chemical properties of lettuce when applied to soil and leaves. BMC Plant Biol. 24:1189. doi: 10.1186/s12870-024-05980-y,
Díaz-Rodríguez, A. M., Parra Cota, F. I., Cira Chávez, L. A., García Ortega, L. F., Estrada Alvarado, M. I., Santoyo, G., et al. (2025). Microbial inoculants in sustainable agriculture: advancements, challenges, and future directions. Plants 14:191. doi: 10.3390/plants14020191,
Duan, Y., Ma, Z., Jia, Y., Xing, H., Mao, Z., Mao, K., et al. (2025). Bacillus vallismortis LRB-5: a promising biocontrol agent for mitigating apple replant disease through pathogen suppression and growth promotion. Stress Biol. 5:51. doi: 10.1007/s44154-025-00246-5,
Duan, Y., Zhou, Y., Li, Z., Chen, X., Yin, C., and Mao, Z. (2022). Effects of Bacillus amyloliquefaciens QSB-6 on the growth of replanted apple trees and the soil microbial environment. Horticulturae 8:83. doi: 10.3390/horticulturae8010083
Durán, P., Tortella, G., Sadowsky, M. J., Viscardi, S., Barra, P. J., and Mora, M. d. l. L. (2021). Engineering multigenerational host-modulated microbiota against soilborne pathogens in response to global climate change. Biology 10:865. doi: 10.3390/biology10090865
Efthimiadou, A., Katsenios, N., Chanioti, S., Giannoglou, M., Djordjevic, N., and Katsaros, G. (2020). Effect of foliar and soil application of plant growth promoting bacteria on growth, physiology, yield and seed quality of maize under Mediterranean conditions. Sci. Rep. 10:21060. doi: 10.1038/s41598-020-78034-6,
Esitken, A., Pirlak, L., Turan, M., and Sahin, F. (2006). Effects of floral and foliar application of plant growth promoting rhizobacteria (PGPR) on yield, growth and nutrition of sweet cherry. Sci. Hortic. 110, 324–327. doi: 10.1016/j.scienta.2006.07.023
Feng, T., Li, J., Mao, X., Jin, X., Cheng, L., Xie, H., et al. (2025). A comparative analysis of the rhizosphere microbial communities among three species of the Salix genus. PeerJ 13:e19182. doi: 10.7717/peerj.19182,
García, G., Carlin, M., and de Jesus Cano, R. (2025). Holobiome harmony: linking environmental sustainability, agriculture, and human health for a thriving planet and one health. Microorganisms 13:514. doi: 10.3390/microorganisms13030514
Geng, W., Lv, Y., Duan, Y., Wang, H., Jiang, W., Zhang, R., et al. (2022). Preparation of composite microbial culture and its biocontrol effect on apple replant disease. Sci. Hortic. 303:111236. doi: 10.1016/j.scienta.2022.111236
Guardiola-Márquez, C. E., Santos-Ramírez, M. T., Segura-Jiménez, M. E., Figueroa-Montes, M. L., and Jacobo-Velázquez, D. A. (2022). Fighting obesity-related micronutrient deficiencies through biofortification of Agri-food crops with sustainable fertilization practices. Plants 11:3477. doi: 10.3390/plants11243477,
Hacquard, S., Wang, E., Slater, H., and Martin, F. (2022). Impact of global change on the plant microbiome. New Phytol. 234, 1907–1909. doi: 10.1111/nph.18187,
Hilber-Bodmer, M., Schmid, M., Ahrens, C. H., and Freimoser, F. M. (2017). Competition assays and physiological experiments of soil and phyllosphere yeasts identify Candida subhashii as a novel antagonist of filamentous fungi. BMC Microbiol. 17:4. doi: 10.1186/s12866-016-0908-z,
Hirata, R., Fujinawa, Y., and Mine, A. (2025). Toward a plant biotechnological application of phyllosphere bacteria. Plant Biotechnol. 42, 251–264. doi: 10.5511/plantbiotechnology.25.0411a,
Jin, L., Jin, N., Wang, S., Li, J., Meng, X., Xie, Y., et al. (2022). Changes in the microbial structure of the root soil and the yield of Chinese baby cabbage by chemical fertilizer reduction with bio-organic fertilizer application. Microbiol. Spectr. 10, e01215–e01222. doi: 10.1128/spectrum.01215-22,
Joshi, H., Bishnoi, S., Choudhary, P., and Mundra, S. (2019). Role of effective microorganisms (EM) in sustainable agriculture. Int. J. Curr. Microbiol. Appl. Sci. 8, 172–181. doi: 10.20546/ijcmas.2019.803.024
Katsenios, N., Andreou, V., Sparangis, P., Djordjevic, N., Giannoglou, M., Chanioti, S., et al. (2022). Assessment of plant growth promoting bacteria strains on growth, yield and quality of sweet corn. Sci. Rep. 12:11598. doi: 10.1038/s41598-022-16044-2,
Khoiri, A. N., Cheevadhanarak, S., Jirakkakul, J., Dulsawat, S., Prommeenate, P., Tachaleat, A., et al. (2021). Comparative metagenomics reveals microbial signatures of sugarcane phyllosphere in organic management. Front. Microbiol. 12:623799. doi: 10.3389/fmicb.2021.623799,
Kim, J. S., Chung, H., Park, B., Veerappan, K., and Kim, Y.-K. (2024). Chloroplast genomesequencing and divergence analysis of 18 Pyrus species: insights into intron length polymorphisms and evolutionary processes. Front. Genet. 15:1468596. doi: 10.3389/fgene.2024.1468596,
Kordatzaki, G., Katsenios, N., Giannoglou, M., Andreou, V., Chanioti, S., Katsaros, G., et al. (2022). Effect of foliar and soil application of plant growth promoting bacteria on kale production and quality characteristics. Scientia Hortic. 301:111094. doi: 10.1016/j.scienta.2022.111094
Kui, L., Xiang, G., Wang, Y., Wang, Z., Li, G., Li, D., et al. (2021). Large-scale characterization of the soil microbiome in ancient tea plantations using high-throughput 16S rRNA and internal transcribed spacer amplicon sequencing. Front. Microbiol. 12:745225. doi: 10.3389/fmicb.2021.745225,
Kumar, S., Kirk, C., Deng, C. H., Shirtliff, A., Wiedow, C., Qin, M., et al. (2019). Marker-trait associations and genomic predictions of interspecific pear (Pyrus) fruit characteristics. Sci. Rep. 9:9072. doi: 10.1038/s41598-019-45618-w,
Lan, J., Dong, X., He, R., Huang, Q., Liu, L., Liu, J., et al. (2025). Impact of rain-shelter cultivation on rhizosphere microecology and kiwifruit quality. Front. Microbiol. 16:1600236. doi: 10.3389/fmicb.2025.1600236,
Leveau, J. H. (2019). A brief from the leaf: latest research to inform our understanding of the phyllosphere microbiome. Curr. Opin. Microbiol. 49, 41–49. doi: 10.1016/j.mib.2019.10.002,
Lin, X., Han, X., Yang, J., Liu, F., Li, Y., and Chen, Z. (2024). Network of soil fungi and the microfauna community under diverse anthropic disturbances under Chrysopogon zizanioides planting in the reservoir. Plants 13:393. doi: 10.3390/plants13030393,
Liu, Z., Luo, C., Zheng, K., Sun, Y., Ru, J., Ma, Y., et al. (2025). Native mixed microbe inoculants (M1H) optimize soil health to promote Cajanus cajan growth: the soil fungi are more sensitive than bacteria. Front. Microbiol. 16:1521064. doi: 10.3389/fmicb.2025.1521064,
Lozano, L., Iglesias, I., Puy, J., and Echeverria, G. (2023). Performance of an expert sensory panel and instrumental measures for assessing eating fruit quality attributes in a pear breeding programme. Foods 12:1426. doi: 10.3390/foods12071426,
Marco, S., Loredana, M., Riccardo, V., Raffaella, B., Walter, C., and Luca, N. (2022). Microbe-assisted crop improvement: a sustainable weapon to restore holobiont functionality and resilience. Hortic. Res. 9:uhac160. doi: 10.1093/hr/uhac160,
Mehlferber, E. C., Debray, R., Conover, A. E., Sherman, J. K., Kaulbach, G., Reed, R., et al. (2023). Phyllosphere microbial associations improve plant reproductive success. Front. Plant Sci. 14:1273330. doi: 10.3389/fpls.2023.1273330,
Mishu, N. J., Hasan, M. R., Islam, S. M. N., Nayeema, J., and Hossain, M. M. (2025). Additive effects of Trichoderma isolates for enhancing growth, suppressing southern blight and modulating plant defense enzymes in tomato. PLoS One 20:e0329368. doi: 10.1371/journal.pone.0329368,
Nzimande, S. N., Sanusi, I. A., Yobo, K., Ramchuran, S. O., and Kana, G. E. B. (2025). Process development for antifungal production by Bacillus subtilis BS20: nanoparticle supplementation, process optimization and preliminary scale-up studies. Bioprocess Biosyst. Eng. 48, 1707–1718. doi: 10.1007/s00449-025-03205-6,
Qin, C., Tao, J., Liu, T., Liu, Y., Xiao, N., Li, T., et al. (2019). Responses of phyllosphere microbiota and plant health to application of two different biocontrol agents. AMB Express 9:42. doi: 10.1186/s13568-019-0765-x,
Qin, L., Tian, P., Cui, Q., Hu, S., Jian, W., Xie, C., et al. (2021). Bacillus circulans GN03 alters the microbiota, promotes cotton seedling growth and disease resistance, and increases the expression of phytohormone synthesis and disease resistance-related genes. Front. Plant Sci. 12:644597. doi: 10.3389/fpls.2021.644597,
Soliman, S. A., Abdelhameed, R. E., and Metwally, R. A. (2023). In vivo and in vitro evaluation of the antifungal activity of the PGPR Bacillus amyloliquefaciens RaSh1 (MZ945930) against Alternaria alternata with growth promotion influences on Capsicum annuum L. plants. Microb. Cell Factories 22:70. doi: 10.1186/s12934-023-02080-8,
Sun, W., and Shahrajabian, M. H. (2023). The application of arbuscular mycorrhizal fungi as microbial biostimulant, sustainable approaches in modern agriculture. Plants 12:3101. doi: 10.3390/plants12173101,
Trabelsi, D., and Mhamdi, R. (2013). Microbial inoculants and their impact on soil microbial communities: a review. Biomed. Res. Int. 2013, 1–11. doi: 10.1155/2013/863240,
Wu, R., Li, Y., Meng, J., and Han, J. (2025). Effects of the addition of microbial agents after Dazomet fumigation on the microbial community structure in soils with continuous cropping of strawberry (Fragaria × Ananassa Duch.). Microorganisms 13:1178. doi: 10.3390/microorganisms13061178,
Xing, W., Gai, X., Xue, L., Li, S., Zhang, X., Ju, F., et al. (2024). Enriched rhizospheric functional microbiome may enhance adaptability of Artemisia lavandulaefolia and Betula luminifera in antimony mining areas. Front. Microbiol. 15:1348054. doi: 10.3389/fmicb.2024.1348054,
Xu, H., Yang, Y., Tian, Y., Xu, R., Zhong, Y., and Liao, H. (2020). Rhizobium inoculation drives the shifting of rhizosphere fungal community in a host genotype dependent manner. Front. Microbiol. 10:3135. doi: 10.3389/fmicb.2019.03135,
Zhang, C., Cui, L., and Fang, J. (2022). Genome-wide association study of the candidate genes for grape berry shape-related traits. BMC Plant Biol. 22:42. doi: 10.1186/s12870-022-03434-x,
Keywords: bacterial communities, fruit quality, microbial agents, phyllosphere, rhizosphere
Citation: Luo N, Zhang Y, Ren Z, Wang X, Li H and Zhang A (2026) Microbial agents enhance the yield and quality of pears by regulating the composition and networks of microbial communities in the phyllosphere and rhizosphere. Front. Microbiol. 17:1763579. doi: 10.3389/fmicb.2026.1763579
Edited by:
Lei Jiang, Zhejiang Agriculture and Forestry University, ChinaReviewed by:
Yang Lyu, China Agricultural University, ChinaXinxin Wang, Hebei Agricultural University, China
Copyright © 2026 Luo, Zhang, Ren, Wang, Li and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhifeng Ren, YmpyemZAMTI2LmNvbQ==; Aiping Zhang, emhhbmdhaXBpbmdAY2Fhcy5jbg==
†These authors have contributed equally to this work
Yulou Zhang1†