Effect of Humic Acid on Soil Physical and Chemical Properties, Microbial Community Structure, and Metabolites of Decline Diseased Bayberry

In recent years, bayberry decline disease has caused significant damage to the bayberry industry. In order to evaluate whether humic acid can be used to effectively control the disease, this research examined the nutritional growth and fruit quality of bayberry, soil physical and chemical properties, soil microbial community structure, and metabolites. Results indicated that the application of humic acid not only improved the vigor and fruit quality of diseased trees, but also increased the diversity of microbial communities in the rhizosphere soil. A great increase was observed in the relative abundance of bacterial genus Mycobacterium and Crossiella; fungal genus Fusarium and Coniosporium. In contrast, a significant decrease was observed in the relative abundance of bacterial genus Acidothermus, Bryobacter, Acidibacter, fungal genus of Geminibasidium and Mycena. Analysis of redundancies (RDA) for microbial communities and soil characteristics showed that the main four variables, including available nitrogen, phosphorus, potassium, and calcium, had a great effect on the composition of bacterial and fungal communities in bayberry rhizosphere soil at the genus level. The main four variables had a greater effect on bacterial communities than on fungal communities. In addition, ABC transporter, arginine and proline metabolism, galactose metabolism, and glutathione metabolism were significantly affected by humic acid, which changed the content of 81 metabolites including 58 significantly down-regulated metabolites such as isohexonic acid and carinitine, and 23 significantly up-regulated metabolites such as acidic acid, guaninosuccinate, lyxose, 2-monoolein, epicatechin, and pentonolactone. These metabolites also significantly correlated with rhizosphere soil microbiota at the phylum, order, and genus levels. In conclusion, the results demonstrated the role of humic acid on plant growth and fruit quality, as well as rhizosphere soil characteristics, microbiota, and secondary metabolites, which provides novel insights into the control of bayberry decline disease.


Introduction
Myrica rubra is an important fruit tree in southern China, with a planting area of about 334,000 hectares and an annual output of about 950,000 tons. It is also distributed in Japan, France, Kenya, the United States, and Brazil [1,2]. Bayberry fruit is rich in juice, and it is sour, sweet, and delicious. It not only contains carbohydrates, organic acids, proteins, and other substances, but is also rich in vitamin C and anthocyanin

Effect of Humic Acid on Vegetative Growth and Fruit Quality
The findings of this study demonstrated that the vitality of weakened trees was considerably enhanced compared to the control after one year of application of humic acid. In fact, the treatment by humic acid resulted in increases of 28.62%, 11.03%, 12.13%, 14.86%, 10.68%, and 17.75% in the branch length, branch diameter, leaf length, leaf width, leaf thickness, and chlorophyll content, respectively, as compared to the control (Table 1). Similarly, the treatment by humic acid significantly increased fruit weight, soluble solids, total sugar, and vitamin C content by 36.30%, 16.19%, 24.72%, and 173.15%, but significantly decreased titratable acid levels by 50% ( Table 2). This data revealed that humic acid had a mitigating effect on bayberry decline disease.  35 ± 0.14 * LD and HA represent decline diseased bayberry in the absence and presence of humic acid, respectively. "*" shows significant increases compared to the control of the decline diseased trees (p < 0.05). Table 2. Effects of humic acid on fruit quality of decline diseased bayberry.
In agreement with the result of this study, our previous studies have shown that HA had the same effect on bayberry as fertilizer and fungicide on the growth of the plant and the quality of the fruit [3,22]. Furthermore, these data are also comparable to the findings of Nardi et al. [11], who found that humic acid can directly impact the photosynthesis of plants. In addition, the application of humic acid can significantly increase the total root length, root surface area, root tip number, and root volume of cucumber seedlings under low nitrogen conditions and promote plant height, stem diameter, and leaf area [23].

Effect of Humic Acid in Microbial Community Diversity
There was a difference in the number of operational taxonomical units (OTUs) between the presence and the absence of humic acid in the bacterial V3 + V4 region and fungal ITS region of the decline diseased bayberry trees. In fact, the average number of bacterial OTUs in the absence of and presence of humic acid are 1352.33 (range: 1202 to 1497) and 1564.83 (range: 1315 to 1764) (Figure 1), while the number of bacterial OTUs in the rhizo-sphere soil of decline diseased trees in the presence of humic acid is 15.17% higher than that of the absence of humic acid.
tween the presence and the absence of humic acid in the bacterial V3 + V4 fungal ITS region of the decline diseased bayberry trees. In fact, the average bacterial OTUs in the absence of and presence of humic acid are 1352.33 (ran 1497) and 1564.83 (range: 1315 to 1764) (Figure 1), while the number of bacter the rhizosphere soil of decline diseased trees in the presence of humic acid higher than that of the absence of humic acid. The average number of fungal OTUs in decline diseased bayberry trees (varying from 594 to 886) and 790.33 (varying from 761 to 845), respectively, in and the presence of humic acid ( Figure 1). In contrast, the number of fungi O zosphere soil in decline diseased bayberry trees was 5.28% higher in the presen acid than in the absence of humic acid. Therefore, it suggests that humic acid h effect on the bacterial operational taxonomic units than that of fungi of rhizosp decline diseased bayberry trees.
Compared with the decline diseased trees, the application of humic aci 11.03% and 3.09%, respectively, increase in the Chao1 index and Shannon ind ria; a 2.76% and 4.74%, respectively, increase in the Chao1 index and Shann fungi in the rhizosphere of bayberry ( Figure 1). The result of this study sugge mic acid had a greater effect on the bacterial richness of rhizosphere soil of de eased trees than on fungi. In agreement with the result of this study, our previ have shown that HA had the same effect on the variety and richness of the r soil of decline diseased bayberry as a compound fertilizer and fungicide [22,2

Effect of Humic Acid in Soil Microbial Community Structure
According to the principal coordinate's analysis (PCoA) of the bacterial structure, the six replicates of the diseased control and humic acid treatment into two groups. The diseased control was well separated from humic acid The average number of fungal OTUs in decline diseased bayberry trees was 750.67 (varying from 594 to 886) and 790.33 (varying from 761 to 845), respectively, in the absence and the presence of humic acid ( Figure 1). In contrast, the number of fungi OTUs of rhizosphere soil in decline diseased bayberry trees was 5.28% higher in the presence of humic acid than in the absence of humic acid. Therefore, it suggests that humic acid had a greater effect on the bacterial operational taxonomic units than that of fungi of rhizosphere soil in decline diseased bayberry trees.
Compared with the decline diseased trees, the application of humic acid caused a 11.03% and 3.09%, respectively, increase in the Chao1 index and Shannon index of bacteria; a 2.76% and 4.74%, respectively, increase in the Chao1 index and Shannon index of fungi in the rhizosphere of bayberry ( Figure 1). The result of this study suggests that humic acid had a greater effect on the bacterial richness of rhizosphere soil of declining diseased trees than on fungi. In agreement with the result of this study, our previous studies have shown that HA had the same effect on the variety and richness of the rhizosphere soil of decline diseased bayberry as a compound fertilizer and fungicide [22,24].

Effect of Humic Acid in Soil Microbial Community Structure
According to the principal coordinate's analysis (PCoA) of the bacterial community structure, the six replicates of the diseased control and humic acid treatments were split into two groups. The diseased control was well separated from humic acid treatments, demonstrating that humic acid considerably affected the bacterial community structure of rhizosphere soil (Figure 2A). Furthermore, PCoA analysis of fungal community structure suggested that six replicates of the diseased control and humic acid treatment were divided into two distinct groups ( Figure 2B). In addition, six replicates of the humic acid treatments revealed a better bacterial diversity compared to the diseased control than those of the fungal community structure (Figure 2). In accordance with the findings of this study, the structure and percentage contribution of soil microorganisms differed between fertilizer or fungicide-treated and untreated soil [22,24].
x FOR PEER REVIEW  Figure 3A). Compared with the decline diseas humic acid treatment significantly increased the relative abundance of Mycobacte Crossiella by 16.06% and 25.32%, respectively; while it significantly reduced th abundance of Acidothermus, Bryobacter, and Acidibacter by 23.74%, 36.38% and 30 spectively ( Figure 3A). Interestingly, the beneficial effect of Mycobacterium and on plant growth has been documented in previous studies [25,26]. For example, exhibits bacteriostatic function [25], whereas Mycobacterium is helpful to plan Mycobacterium Mya-zh01 can efficiently synthesize and secrete the plant growth indole-3-acetic acid, resulting in the increase of the number and length of ro height, number of leaves, and leaf length of Doritaenopsis [26]. This finding demonstrated that the application of humic acid resulted in a significant change in the bacterial and fungal community structure at the phylum ( Figure S1), order ( Figure S2), and genus levels ( Figure 3) compared to the decline disease control. Among the top 15 bacterial genera, Acidothermus, Mycobacterium, Bryobacter, Acidibacter, and Crossiella were the main five genera ( Figure 3A). Compared with the decline disease control, humic acid treatment significantly increased the relative abundance of Mycobacterium and Crossiella by 16.06% and 25.32%, respectively; while it significantly reduced the relative abundance of Acidothermus, Bryobacter, and Acidibacter by 23.74%, 36.38% and 30.50%, respectively ( Figure 3A). Interestingly, the beneficial effect of Mycobacterium and Crossiella on plant growth has been documented in previous studies [25,26]. For example, Crossiella exhibits bacteriostatic function [25], whereas Mycobacterium is helpful to plant growth. Mycobacterium Mya-zh01 can efficiently synthesize and secrete the plant growth hormone indole-3-acetic acid, resulting in the increase of the number and length of roots, plant height, number of leaves, and leaf length of Doritaenopsis [26].
Additionally, putting Mya-zh01 on Doritaenopsis seeds made them grow faster [26]. In vitro and in silico studies revealed that the strains prevented pathogenic Gram-positive and Gram-negative bacteria and fungi from growing. These results indicated that Mycobacterium and Crossiella may be advantageous for bayberry growth. The Acidibacter reduction is similar to the declined diseased bayberry treated by compound fertilizer and bio-organic fertilizer [3]. However, in contrast with the result of this study, Bryobacter has been found to contribute to high maize yield by increasing phosphorus flow and inhibiting toxic aluminum and manganese flow from soils to plants [27]. Among the top 15 fungi, Cladophialophora, Geminibasidium, Mycena, Fusarium, and Coniosporium were the main five genera ( Figure 3B).
The relative abundance of Fusarium and Coniosporium increased significantly by 80.74% and 29.87%, respectively, compared with the decline diseased trees; the relative abundances of Geminibasidium and Mycena were significantly reduced by 30.71% and 86.27%, respectively, and no significant change was observed on Cladophialophora ( Figure 3B). Fusarium spp. is a well-known soil-borne plant pathogen that causes severe vascular wilt in economically important crops worldwide [28]. The function of the Coniosporium has not been reported. Geminibasidium reduction in decline disease bayberry rhizosphere soil is comparable to a previous study in the rhizosphere soil [3] but different from compound fertilizer and bio-organic fertilizer treated soil [24]. By increasing nitrogen absorption and NH4+ assimilation, the plant growth-promoting fungus MF23 (Mycena sp.) increases the yield of Dendrobium officinale (Orchidaceae) [29]. Compared to the previous reports, spectively ( Figure 3A). Interestingly, the beneficial effect of Mycobacterium and on plant growth has been documented in previous studies [25,26]. For example exhibits bacteriostatic function [25], whereas Mycobacterium is helpful to pla Mycobacterium Mya-zh01 can efficiently synthesize and secrete the plant growth indole-3-acetic acid, resulting in the increase of the number and length of r height, number of leaves, and leaf length of Doritaenopsis [26]. Additionally, putting Mya-zh01 on Doritaenopsis seeds made them grow In vitro and in silico studies revealed that the strains prevented pathogenic Gra and Gram-negative bacteria and fungi from growing. These results indicated bacterium and Crossiella may be advantageous for bayberry growth. The Acidiba tion is similar to the declined diseased bayberry treated by compound fertilize organic fertilizer [3]. However, in contrast with the result of this study, Bryobact found to contribute to high maize yield by increasing phosphorus flow and toxic aluminum and manganese flow from soils to plants [27]. Among the top Cladophialophora, Geminibasidium, Mycena, Fusarium, and Coniosporium were the genera ( Figure 3B). In conclusion, beneficial bacteria increased while harmful bacteria decreased in bayberry root soil. Beneficial bacteria can utilize organic matter, degrade pollutants, improve stress tolerance, and enhance the growth of diseased bayberry trees to make them develop towards the level of healthy trees, reduce bayberry disease, and promote healthy bayberry tree growth.

Effect of Humic Acid on Soil Nutrient Status
Compared with the decline in diseased control, the content of pH, organic matter, alkali hydrolyzed nitrogen, exchangeable calcium, and exchangeable magnesium increased by 7.74%, 174.82%, 231.91%, 335.93%, and 316.10%, respectively, after applying humic acid treatment. In comparison, the content of available phosphorus and available potassium decreased by 86.48% and 72.96%, respectively (Table 3). Fourier-transform infrared (FTIR) spectroscopy revealed the complexation of Zn sulfate by HA through its S and C functional groups. In both oxisols, solution Zn increased due to the combined use of Zn and HA [30].
In the oxisol with lower organic matter content, HA can assure an adequate supply of residual Zn while increasing the growth of brachiaria cultivated in sequence with maize. Calcium phosphate particles coated with humic substances are a potential plant biostimulant from the circular economy. The obtained material showed promising results in its potential to elicit phosphorous uptake and foliar translocation by plants [31]. Humic acid reduces harmful aluminum ions in the soil, lowers soil acidity, promotes the release of alkali hydrolyzable nitrogen, exchangeable calcium, and exchangeable magnesium, and improves plant uptake, utilization, and transport of phosphorus and potassium. Improved bayberry soil physical, chemical, and nutritional balance. LD and HA represent the absence and presence of humic acid in diseased bayberry trees, respectively. "*" and "#" represents significant (p < 0.05) increase and decrease compared to the control of the decline diseased trees, respectively.

Effect of Humic Acid on RDA of Soil Properties and Microbial Communities
The composition of bacterial and fungal communities in bayberry rhizosphere soil was affected significantly by different soil physical and chemical properties at the genus level ( Figure 4) ( Table 4). In the RDA result, the first and second axes explained 71.13% and 75.78% of the cumulative variance of the rhizosphere microbial community-factor correction, respectively, at the bacterial genus level ( Figure 4A), 31.22% and 39.28% of the cumulative variance of the rhizosphere microbial community-factor correction, respectively, at the fungal genus level ( Figure 4B). For the bacterial communities at the genus level, the contributions of four main variables were shown as below: 28.8% for exchangeable calcium, 25.0% for available phosphorus, 16.5% for available potassium, 13.0% for available nitrogen, and 10.0% for pH. The contribution of organic matter to bacteria (4.2%) was smaller than to fungi (32.6%), which may be because the decline disease was caused by a fungal pathogen, resulting in a greater number of fungi than bacteria in the rhizosphere of diseased bayberry trees [32].
The contributions of four main variables, including organic matter, available phosphorus, pH and available potassium, explained 32.6%, 18.5%, 14.7% and 9.4% of the fungal community, respectively, at the genus level. Results showed that the four main variables, including available nitrogen, phosphorus, potassium, and calcium, significantly affected the composition of bacterial and fungal communities in bayberry rhizosphere soil at the genus level. In agreement with the result of this study, a previous study has also shown that the growth of soil microbe is affected by a series of environmental factors such as soil pH, organic matter content, magnesium content, available nitrogen, phosphorus, and calcium [33]. However, this study revealed that the four main variables had a greater effect on bacterial communities than on fungal communities.
The result of this study revealed the complexity of the relationship between microbial growth and soil nutrient elements. Furthermore, the result obtained demonstrated that the order of contribution rates of soil physical and chemical indexes and nutrient elements to bacteria and fungi is different from that of untreated healthy trees, diseased trees, fertilizertreated diseased trees and fungicide-treated diseased trees [3,22,24]. This data suggests that a balance of nutrients may be the best for microbial growth. cium, 25.0% for available phosphorus, 16.5% for available potassium, 13.0% for available nitrogen, and 10.0% for pH. The contribution of organic matter to bacteria (4.2%) was smaller than to fungi (32.6%), which may be because the decline disease was caused by a fungal pathogen, resulting in a greater number of fungi than bacteria in the rhizosphere of diseased bayberry trees [32].

Change in Rhizosphere Soil Metabolomics
Using GC-MS analysis, 223 rhizosphere soil metabolites were identified in this study. Furthermore, a score map of metabolites was made using orthogonal partial least squaresdiscriminant analysis (OPLS-DA, a supervised statistical method of discriminant analysis), which can reduce intragroup differences, increase intergroup differences, and remove the impact of irrelevant factors on the experimental data to make accurate predictions of different samples. According to these findings, humic acid exhibited a great influence in the metabolomics distribution of rhizosphere soils from decline diseased bayberry. Indeed, in the absence of HA, the decline diseased samples were all distributed in the negative area of t [1] (principal component 1).
In contrast, in the presence of HA, the decline disease samples were all distributed in the positive area of t [1] ( Figure 5). Furthermore, the parameters of the OPLS-DA are R 2 X(cum) = 0.574, R 2 Y(cum) = 0.997, Q 2 (cum) = 0.881. The Q 2 values were greater than 0.5, which showed that the model had good interpretation and prediction ability. The cluster separation effect of LD and HA indicated that the application of humic acid significantly changed the metabolic structure of the rhizosphere soil of the decline-diseased bayberry.

Analysis of Differential Metabolites
The cluster analysis showed the contents and chang soil under LD and HA, while the major metabolites with and HA were sugar, amino, organic, and secondary. Am metabolites, 23 were up-regulated by 11.30%-3547.82%, w and 58 were down-regulated by 15.55%-79.49%, with d ( Figure 6; Table 5). The 19 metabolites with more than 50 were palatinitol, 1-monoolein, zymosterol, deoxycho melezitose, isohexonic acid, carinitine, proline, 1-kestose the exogenous application of 1 mM carnitine mitigates th increasing mitosis and decreasing DNA damage cause

Analysis of Differential Metabolites
The cluster analysis showed the contents and changes of metabolites in rhizosphere soil under LD and HA, while the major metabolites with significant changes between LD and HA were sugar, amino, organic, and secondary. Among the 81 significantly changed metabolites, 23 were up-regulated by 11.30%-3547.82%, with aconitic acid being the most, and 58 were down-regulated by 15.55%-79.49%, with deoxycholic acid being the least ( Figure 6 and Table 5). The 19 metabolites with more than 50% reduction in relative contents were palatinitol, 1-monoolein, zymosterol, deoxycholic acid, hexitol, maltotriose, melezitose, isohexonic acid, carinitine, proline, 1-kestose, Tagatose and so on. However, the exogenous application of 1 mM carnitine mitigates the harmful effects of salt stress by increasing mitosis and decreasing DNA damage caused by oxidative stress on barley seedlings [34].
Sci. 2022, 23, x FOR PEER REVIEW Figure 6. Thermogram analysis of different metabolites in bayberry HA1-6 represent the absence and presence of humic acid in diseased b Table 5. The relative contents of the metabolites in rhizosphere soil we the decline diseased bayberry trees.

Metabolite Name
Relative Content Metabolite Name Figure 6. Thermogram analysis of different metabolites in bayberry rhizosphere soil. LD1-6 and HA1-6 represent the absence and presence of humic acid in diseased bayberry trees, respectively.   LD and HA represent the absence and presence of humic acid in diseased bayberry trees, respectively. "*" and "#" represents significant (p < 0.05) increase and decrease compared to the control of the decline diseased trees, respectively.
Proline is a multifunctional amino acid involved in plant adaptation to environmental constraints [35]. As the main fructooligosaccharide, 1-kestose shares similar physiological effects with other fructooligosaccharides. They have recently been determined to show more notable effects in promoting the growth of probiotics, including Faecalibacterium prausnitzii and Bifidobacterium, than those of other fructooligosaccharides [36]. Tagatose is a rare sugar decomposed by a few microorganisms and inhibits a broad range of phytopathogens [37]. In this study, the application of humic acid may increase the absorption and utilization of these substances by plants. Hence, the content of these secondary metabolites in the soil is lower than that of the control. The 23 metabolites with increased relative content may be beneficial to plant growth. The 12 metabolites with >50% increase of relative levels increased were 5-methoxytryptamine, diclofenac, pentonolactone, epicatechin, 2-monoolein, lyxose, 1-hexadecanol, guanidinosuccinate, beta-sitosterol, coniferin, aconitic acid, and dehydroabietic acid. The six relative contents of aconitic acid, guanidinosuccinate, lyxose, 2-monoolein, epicatechin, and pentonolactone were increased by 1101.31%-3547.28%. Interestingly, diclofenac could potentially act as an environmental contaminant disturbing the natural developmental processes of plants [38]. Furthermore, beta-sitosterol has been regarded as an important plant growth regulator [39]. Therefore, the effect of humic acid on the growth of bayberry trees may be, at least partially, attributed to the increased release of the above metabolites in rhizosphere soil.
Compared to control cells, salt-adapted cells accumulated higher sugars, amino acids, and intermediate metabolites in the shikimate pathway, such as coniferin [40]. Aconitic acid plays various biological roles within cells as an intermediate in the tricarboxylic acid cycle. It confers unique survival advantages to some plants as an antifeedant, antifungal, and means of storing fixed pools of carbon [41]. Aconitic acid has also been reported as an inhibitor of fermentation, an anti-inflammatory, and a possible nematicide [41]. The time-kill assay with 6.25 g·ml −1 of dehydroabietic acid showed that it killed Staphylococcus epidermidis (American Type Culture Collection 14,990) bacteria within 24 h [42]. The application of humic acid may increase the release of beneficial substances from the soil, reduce the absorption of toxic substances by plants, and increase the absorption of beneficial substances.

Effect of Humic Acid on Metabolic Pathways
Using the KEGG (Kyoto Encyclopedia of Genes and Genomes) database, a pathway enrichment analysis of various metabolites in LD and HA was performed. There were significant differences between LD and HA in 10 metabolic pathways, including ABC transporters (p < 0.01), aminoacyl-tRNA biosynthesis (p < 0.01), galactose metabolism (p < 0.01), glutathione metabolism (p < 0.01), butanoate metabolism (p < 0.01), cyanoamino acid metabolism (p < 0.01), arginine and proline metabolism (p < 0.01), phosphotransferase system (PTS) (p < 0.05), amino sugar and nucleotide sugar metabolism (p < 0.05), and glycine, serine and threonine metabolism (p < 0.05) (Figure 7). The results demonstrated that the application of humic acid increased the concentration of metabolically active metabolites. The amino acid metabolism, sugar metabolism, and nucleotide metabolism are the principal metabolic pathways of the differential metabolites in rhizosphere soil from diseased bayberry trees between the absence and presence of humic acid. Interestingly, ABC transporters not only have a detoxification mechanism, and play an important role in the response of plant pathogenic microorganisms, regulation of heavy metals, and secondary metabolite transport, but also affect the absorption of nutrients by bacteria [43]. Furthermore, glutathione metabolism has a detoxifying effect, and the increase of its metabolic level is conducive to plant growth and microbial reproduction.
Arginine and proline are important amino acids in plants, while arginine has the functions of ammonia detoxification, hormone secretion, immune regulation, and nitro-gen storage. At the same time, proline is not only the energy which restores growth, but also has the function of maintaining protoplasm, environmental osmotic balance, and preventing water loss, and its metabolism also plays a good role in cold resistance, drought resistance and salt tolerance [44,45]. Therefore, the application of humic acid may improve the metabolic level of debilitated bayberry by promoting the transport of substances in the body and the absorption of nutrients by bacteria in the soil through ABC transporters.
Through galactose metabolism and glutathione metabolism, humic acid can inhibit harmful microorganisms and promote the growth of beneficial microorganisms, so as to improve the growth status of weak and diseased bayberry trees. Through the metabolism of arginine and proline, the resistance to cold stresses will be improved, which was able to promote the recovery and healthy growth of the root system of bayberry. The significant change in galactose metabolism after humic acid treatment is the same as those of diseased trees treated with fungicides and fertilizers, but the number of pathways was greater than that of fungicides and fertilizers [22,24].

Correlation of Soil Microorganisms with Metabolites
The relative content of bacteria at the phylum ( Figure S3), order ( Figure S4), and genus ( Figure 8) levels are correlated with the relative content of the main secondary metabolites in the soil. The bacterial genera closely related to the content of secondary metabolites were Jatrophihabitans, Bryobacter, Candidatus Soliactor, and Rhiodiomicrobium, which were significantly related to 33, 23, 21, and 12 metabolites, respectively. In detail, Jatrophihabitans and Rhizomacrobium were negatively correlated with 19 and 12 metabolites, respectively, with reduced relative content and positively correlated with 14 and 0 metabolites, respectively, with increased relative content.
Furthermore, Bryobacter and Candidatus Soliactor were positively correlated with 17 and 13 metabolites, respectively, with reduced relative content and negatively correlated with 6 and 8 metabolites, respectively, with increased relative content. Crossella and Mycobacterium with significantly increased relative content were correlated with 1 and 0 of 58 metabolites, respectively. Acidibacter and Acidothiermus with significantly reduced relative content were significantly correlated with the contents of 5 and 4 secondary metabolites, respectively. This data shows that some functional bacteria play an important role in the change of the content of secondary metabolites.
The relative content of fungi at the phylum ( Figure S5), order ( Figure S6), and genus ( Figure 9) levels are correlated with the relative content of the main secondary metabolites in the soil. The fungi closely related to the content of secondary metabolites are Mortierella, Umbelopsis, Rickenella, Penicillium, Fusarium and Geminibasidium, which are significantly related to 39,38,27,15,13, and 10 metabolites, respectively. There were 17, 7, and 3 species of Mortierella, Penicillium, and Fusarium that had a significant positive correlation with metabolites with increased relative content, and 22, 8, and 10 species of Mortierella, Penicillium, and Fusarium that had a significant negative correlation with metabolites with decreased relative content, respectively.
Umbelopsis, Rickenella, and Geminibasidium have 6, 4, and 0 species with significant negative correlation with metabolites with increased relative content, and 32, 23, and 10 species with significant positive correlation with metabolites with decreased relative content, respectively. This result shows that some functional fungi play an important role in the change of the content of secondary metabolites. Figure 8. Correlation analysis between the microorganism relative abundances at the ba nus level and the relative metabolite contents of the humic acid treatment. "*" and "**" re significant correlation at p < 0.05 and p < 0.01, respectively. The magnitude of the correla cient was indicated by the depth of the orange and blue scale, while the orange darker c greater positive correlation, and the blue darker color is the greater negative correlation.
Furthermore, Bryobacter and Candidatus Soliactor were positively correlate Figure 8. Correlation analysis between the microorganism relative abundances at the bacterial genus level and the relative metabolite contents of the humic acid treatment. "*" and "**" represents a significant correlation at p < 0.05 and p < 0.01, respectively. The magnitude of the correlation coefficient was indicated by the depth of the orange and blue scale, while the orange darker color is the greater positive correlation, and the blue darker color is the greater negative correlation. related to 39,38,27,15,13, and 10 metabolites, respectively. There were 17, 7, a of Mortierella, Penicillium, and Fusarium that had a significant positive corre metabolites with increased relative content, and 22, 8, and 10 species of Mort cillium, and Fusarium that had a significant negative correlation with metaboli creased relative content, respectively. Figure 9. Correlation analysis between the microorganism relative abundances at the level and the relative metabolite contents of the humic acid treatment. "*" and "**" significant correlation at p < 0.05 and p < 0.01, respectively. The magnitude of the corre cient was indicated by the depth of the orange and blue scale, while the orange darke greater positive correlation, and the blue darker color is the greater negative correlatio Figure 9. Correlation analysis between the microorganism relative abundances at the fungal genus level and the relative metabolite contents of the humic acid treatment. "*" and "**" represents a significant correlation at p < 0.05 and p < 0.01, respectively. The magnitude of the correlation coefficient was indicated by the depth of the orange and blue scale, while the orange darker color is the greater positive correlation, and the blue darker color is the greater negative correlation.

Experimental Design
This study was performed in July 2019 on 15-year-old Dongkui bayberry trees from Majian Town, Lanxi City. The tested bayberry trees with the similar load, crown size, leaf loss (10-25% of the whole plant's leaves) and disease index (grade 3) [3]. The distance between planting rows was about 4 m × 5 m, while the orchard was managed routinely in the traditional manner. After harvesting the fruit, each decline diseased bayberry tree was treated by using 0 and 1400 g/plant of humic acid (90% effective content, Zhejiang Fengyu Ecological Technology Co., Ltd., Pujiang, China), which was applied in the drip line ditch of the roots. Each replicate consisted of a single tree, and each treatment contained 6 trees with similar vegetative and reproductive growth. The effect of humic acid on decline diseased bayberry trees was evaluated in June 2020 by measuring substantial variations in bayberry's spring shoots, leaves, fruits, and the microbial community structure of rhizosphere soils.

Measurement of Vegetative Growth Parameters
The collection and measurement of branch and leaf samples were carried out, as described by Ren et al. [3]. In detail, the diameters of 20 randomly selected ripe spring branches and 30 randomly selected leaves from 6 trees in each treatment were measured with a ruler and a digital Vernier caliper (Shanghai Daoju, Shanghai, China). Leaves were sampled from the fourth to eighth leaves down from the top of the living branches in the middle of the tree's circumference. The width and length (from tip to petiole) of the leaves were measured with a ruler, while the average thickness of leaves was determined with a digital Vernier caliper (Beaverton, OR, USA). The photosynthesis rate was calculated using a LI-6400 hand-held photosynthesis analyzer (LI-COR Company of America, Lincoln, NE, USA).

Measurement of Fruit Economic Characters
The collection and economic characters measurement of fruit samples were carried out, as described by Ren et al. [3]. In brief, 15 mature fruits were randomly selected from 6 trees in each treatment during fruit maturation phase. Fruit samples were routinely stored at −20 • C. After harvest, the weight and soluble solids content of single fruit were determined immediately. The fruit fresh weight was measured using an electronic balance (Shanghai Precision Instrument, Shanghai, China), while the contents of total soluble solids (TSS) was determined using an ATAGOPR-101a hand-held digital glucometer (Tokyo, Japan). In addition, the concentration of titratable acid in fruit samples was measured by acid-base titration, as described by [18], while the vitamin C concentration in fruit samples was quantified by using 2-6 dichloroindophenol titration, as described by Ghimire et al. [19].

Soil Sample Collection and Physical and Chemical Property Measurements
Soil sample collection and physical and chemical property measurements were performed, as described in our previous study [3]. After fruit harvest, rhizosphere soil samples were collected from the six trees in each treatment. A total of 2 kg of mixed soil samples (0-20 cm) were collected using the quartering method from the drip line surrounding the crown of the bayberry plant. After passing through a 0.45 mm filter, one-half sample was kept in a refrigerator at −80 • C for the extraction of genomic DNAs, and the other sample was air dried at room temperature to measure the soil physical and chemical properties. The content of organic matter was determined by using the K2CrO 7 oxidation method, while the pH was measured using a pH meter with a soil-to-water ratio of 1:2.5 [21]. The available N and available P were evaluated with the alkaline hydrolysis diffusion method and the anti-molybdenum antimony colorimetry method, respectively [46]. The available K was determined with the ammonium acetate extraction flame photometer method [47]. After extraction with ammonium acetate and iron, the exchangeable calcium and magnesium were measured with an ice3500 atomic absorption spectrophotometer [20].

Soil Genome Sequencing
Sample processing and microbiome sequencing were performed according to the method of [3]. In brief, the DNAs were extracted from the rhizosphere soil samples by using the E.Z.N.ATM Mag-Bind Soil DNA Kit (Biodiz, Lince, Peru) following the manufacturer's instructions (OMEGA, Norcross, GA, USA). The purity of the extracted DNAs was evaluated using a Nano-Drop (ND-1000) spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Bacterial diversity was determined by amplifying the V3-V4 region of the 16S rRNA gene using the primers 341F and 805R [48], while fungal diversity was determined by amplifying the ITS1 and ITS2 region using the primers ITS1F and ITS2 [3,15], respectively.

Gas Chromatography-Mass Spectrometry (GC-MS) Metabolomics Analysis
For detailed methods of soil sample processing and secondary metabolite detection, we refer to our previously published paper [3]. In order to increase the accuracy of the entire analysis, we also included one quality control (QC) sample for every 10 samples. The QC sample was prepared by mixing the extracts of all samples in the same amount, and each QC sample was the same amount as the tested sample. Analysis of GC-MS metabolomics of soil samples was conducted on Agilent Technologies Inc. 7890B-5977A GC/MSD GC-MS (Santa Clara, CA, USA), which was carried out, as described in our recent publication [3]. Furthermore, the metabolite information was obtained by analyzing the KEGG database after comparing the data obtained in this study with the standard spectrum library.

Statistical Analysis
Preliminary data were processed using Excel 2016 statistical software with a significance level of 0.05 for t-test, while the collected data were compared with those in the standard spectrum library, which was established by the National Institute of Standards and Technology. Furthermore, PCoA, community histograms and redundancy discriminant analysis (RDA) were performed by using R 3.5.1 (The R Foundation for Statistical Computing, Vienna, Austria). In addition, the heat map was drawn by using R package, Pheatmap version 1.0.12 (https://CRAN.R-project.org/package=pheatmap, accessed on 15 October 2022). Orthogonal partial least squares discriminant analysis (OPLS-DA) was made by using R 3.6.2 ropls. The α-diversity metrics and carry out the significance test (p < 0.05) was calculated using the t-test and SPSS 17.0 software (IBM, Chicago, IL, USA).

Conclusions
The application of humic acid could improve the vigor of diseased trees, the fruit quality, and the diversity and species number of bacterial communities in the rhizosphere soil. A great increase was observed in the relative abundance of bacteria Mycobacterium and Crossiella, and fungi Fusarium and Coniosporium. In contrast, a great decrease was also found in the relative abundance of bacteria Acidothermus, Bryobacter, Acidibacter, and fungi Geminibasidium and Mycena. RDA of microbial communities and soil characteristics revealed that the composition of bacterial and fungal communities in bayberry rhizosphere soil was significantly affected by the main four variables, including available nitrogen, phosphorus, potassium, and calcium, at genus levels.
The four main variables exhibited a greater influence in bacterial communities than in fungal communities. In addition, the significant correlation between microorganisms and secondary metabolites at the genus level confirmed the alteration in the organization of the microbial community. According to the findings, there was a great difference between the absence and presence of humic acid in the metabolomics distribution of rhizosphere soils from decline disease bayberry. In fact, GC-MS metabolomics analysis showed that 81 metabolites in the rhizosphere soil of diseased trees changed significantly after applying humic acid, 23 metabolites such as acidic acid, guanidinosuccinate, lyxose, 2-monoolein, epicatechin, and pentonolactone were significantly up-regulated, and 58 metabolites such as 1-monoolein, zymosterol, deoxycholic acid, hexitol, maltotriose, isohexonic acid and carinitine were significantly down-regulated.
There was a great difference between the absence and presence of humic acid in 10 metabolic pathways, including ABC transporters, galactose metabolism, glutathione metabolism, aminoacyl-tRNA biosynthesis, butanoate metabolism, cyanoamino acid metabolism, arginine and proline metabolism, phosphotransferase system, amino sugar and nucleotide sugar metabolism, and glycine, serine, and threonine metabolism. Overall, this study highlighted the great influence of humic acid in decline disease by regulating soil microbial community, physical and chemical properties, and secondary metabolites in the bayberry rhizosphere, which provided a novel strategy for the management of bayberry decline disease.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ijms232314707/s1, Figure S1: Relative abundance of bacteria (A) and fungi (B) at the phylum level; Figure S2: Relative abundance of bacteria (A) and fungi (B) at the order level; Figure S3: Heatmap of correlation analysis between the microorganism relative abundances at phylum level of bacteria and the metabolite relative contents of the humic acid in diseased bayberry trees treatment; Figure S4: Heatmap of correlation analysis between the microorganism relative abundances at order level of bacteria and the metabolite relative contents of the humic acid in diseased bayberry trees treatment; Figure S5: Heatmap of correlation analysis between the microorganism relative abundances at fungal phylum level of and the metabolite relative contents of the humic acid in diseased bayberry trees treatment; Figure S6: Correlation analysis between the microorganism relative abundances at the fungal order level and the relative metabolite contents of the humic acid treatment. "*" and "**" represents a significant correlation at p < 0.05 and p < 0.01, respectively. The magnitude of the correlation coefficient was indicated by the depth of the orange and blue scale, while the darker orange color is the greater positive correlation, and the darker blue color is the greater negative correlation.