Effects of Soil Amendment with Chinese Herbal Medicine on Cassava Growth and Rhizosphere Microbial Diversity


 This study aimed to determine the effects of different dosages of Chinese herbal compound soil amendment on cassava growth and rhizosphere microbial diversity. In this study, Cassava SC8 growth was determined under four treatments. The growth was promoted by the application of the Chinese herbal compound soil amendment powder and agent. After 240 days of treatment, we measure its height, diameter of the stem, maximum perimeter of cassava, tuber number, weight of single tuber, yield per hectare, the mineral elements in cassava, the ash content in starch and the microbial diversity in the soil. Compared with the control group, the height, diameter of the stem, maximum perimeter of cassava, tuber number, the weight of a single tuber, yield per hectare in the treatment groups all show increase, and the application of 1200 kg of Chinese herbal compound soil amendment per hm2 (T2) had the strongest positive effect on the cassava growth. As for the water content of tubers, that of T2 decreased by 5.96% compared with T4, while no substantial difference among other groups. The Fe content in T2 decreased by 29.46%, and Zn in T3 increased by 39.48%, while the content of other metal ions in four treatments showed no substantial difference. At the same time, the experiment also found that, as bacterial abundance in soil, Streptophyta dominated in T2, and Fimbriimonas is a unique genera of T2. In addition, the abundances of Nocardioides and unclassfied_Spartobacteria in the T2 were greatly considerable than other treatment groups. As fungal abundance in soil, after treated with different concentrations, the composition of soil fungi differed greatly among the groups, and the samples at different concentrations contained more than 20 species of unique bacteria. At the genus level, compared with the control group, the relative abundance of unclassfied_Chaetomiaceae in the treatment group was considerably reduced, and the relative abundance of Psathyrella in the T1 increase substantially. We speculates that these advantages and the existence of bacteria and fungi can affect rhizosphere soil microorganisms, further improve soil effective nutrients and enzyme activities, thereby affecting the growth and physiological indicators of cassava plants.


Introduction
Cassava (Manihot esculenta Crantz) is one of the three key tuber crops in the world, known as the "king of starch," "special crop" or "underground granary." Moreover, it is the food crop in tropical and subtropical areas, mainly used for starch production, and also used in medicine, textile production and energy generation 1 . However, cassava has long been regarded as a crude crop. During the planting process, farmers often apply only a small amount of basal fertilizer or leave the soil unfertilized 2 , which negatively affecting the yield and quality of cassava.
Soil amendments always be treated as a valuable source of nutrients for impoverished agricultural soils, with the potential to improve crop quality. In the recent decades, studies have shown that soil amendment can effectively facilitate plant growth and quality, including the yield, emergence rate, aboveground biomass and other physiological indicators 3,4,5 and improve rhizosphere microbial diversity 6 . Wu 3 applied an environmentally friendly soil amendment (acid soil amendment and nutritive acid soil amendment) developed from agricultural waste resources and humic acid to cassava planting elds in Hainan Province and found that it could substantially improve the yield and other agronomic traits of cassava. In addition, other studies have found that the application of soil amendment applied in lettuce (Lactuca sativa L.), ginkgo (Ginkgo biloba L.) and walnut (Juglans regia) can promote their growth and improve various physiological indexes 7,8,9,10,11,12 . Microorganisms in the soil are highly diverse and abundant 13 . These microbes play a pivotal role in the establishment of microbial diversity. Recent studies has proved many different biotic or abiotic factors will cause the change of microbial community structure and composition which may in uence crop production in certain aspects 14,15 . In the last years, the application of organic matter (OM) in the agricultural production always in form of applying crop residues, manure, compost into soil, which has gained much attention on develop soil amendments as a sustainable approach for the reutilization of these byproducts 16,17,18,19 . For example, livestock manure-derived amendments can be very useful to improve the quality of agricultural soil, as they increase the content of soil OM and stimulate microbial activity and biomass, thus enhancing crop yield 20 .
Although there has been considerable research on the improvement of plant physiological indexes promoted by soil amendments, the effect of Chinese herbal medicine on the growth of cassava has not yet been reported. Taking cassava as the research object, which are the optimal combination selected by multi-site and multi-crop plot tests for many years on crops such as tomatoes, strawberries, tea and Radix pseudostellariae in Guizhou from 2018 to 2020. In this study, we use residue of Chinese herbal medicine as raw materials, making them into powder and water as organic soil amendment. Artemisia artemisiae and A. argyi are common herbs in powder and water. Artemisia carvifolia and Lonicera japonica were applied as the main raw material preparation for the production of crop root soil powder (a local substitute for conventional organic fertilizer or chemical fertilizer), intended to act as a fertilizer and also assist in the prevention and control of plant diseases and insect pests.
Artemisia annua is widely distributed in all parts of China, and its medicinal parts are the dry ground parts. This plant has a considerable growthpromoting effect, and the appropriate concentration of extract can regulate the growth of rice, wheat and other crops, and improve seed vigor and activities of various enzymes 21,22 . Chen also found that a variety of lipid-solvent and water-impregnated compounds in A. annua can regulate and control plant growth 23 . Artemisia annua mother liquor ethanol extract compound has allelopathic effect on wheat, and the appropriate concentration has a strong promoting effect on wheat seed germination 24 . Flavonoids are secondary metabolites of plants and have a variety of important biological functions 25 . According to previous studies by the project team, avonoids are included in the chemical components of bacteriostatic agents and antiviral water agents 12,26 , which may affect the microbial diversity in the soil.
In this study, the residue of A. artemisiae and A. argyi were used to prepare water spraying agent on cassava leaves, and the soil improvement powder was also applied to roots as organic soil amendment. We found that the combination of the water and powder was shown to improve the growth, yield and quality of cassava. This study also aimed to determine the effects of different dosages of soil amendment powder on cassava growth indexes and microbial diversity. The results of this study are intended to provide support for improving cassava cultivation technology and raising production and planting level.

Test materials
The experimental cassava SC8 was planted in the agricultural science base of Hainan (ridges 12 m in length, 1.2 m in width, and 0.6 m between ridges; 0.5 m between plants). When the plant height was about 60 cm, the plant was topdressed with Chinese herbal compound soil improvement powder and agent. The Chinese herbal compound powder was prepared with residue of multiple herbs such as A. carvifolia, Artemisia argyi and Sophora avescens, according to the method described by Wang 12 , and the agent was prepared by the method described by Zheng 26 . The powder and agent were made up according to the formula of Chinese medicine, which was reported by Wang 12 . The raw powder of Chinese medicinal materials residue prepared as an agent was crushed into crude medicinal material powder (≤40 mm), which was mixed with hulled wheat bran and roasted for 45 min. The mixture turned slightly yellow and was then applied to the soil powder. The organic matter content of the powder was 90%, total nitrogen, phosphorus and potassium content was 7.2%, and pH was 5.8, according to the test by Guizhou Provincial Quality Inspection Institute, which was the optimal combination selected based on multi-site and multi-crop plot tests for many years on crops such as tomatoes, strawberries, tea and R. pseudostellariae in Guizhou from 2018 to 2020 12,26 . All materials were provided by Guizhou Yichuan Biotechnology Co., Ltd.
Experimental research and eld studies on plants complies with relevant institutional, national, and international guidelines and legislation Experimental equipment Experimental equipment included an electronic balance, ruler, tape measure and caliper.
Chinese herbal compound soil amendment fertilizer was used as topdressing, and no fertilizer was applied to the control. The agent were mixed, diluted 15 times and then sprayed on cassava regularly (every 10-15 days).
The compound soil amendment powder was used as topdressing, and a distance of application was 15 cm between each seedling to avoid fertilizer burn. After fertilizer application, 2.0 kg of compound antibacterial agent and antiviral agent was taken and diluted 15 times (1 kg of water solution for each 15 L volume sprayer), and then lled with clear water. The spray was applied once every 10-15 days. During spray irrigation, the nozzle of the sprayer was put to the root and sprayed for about 3 seconds, and then applied to the leaf surface. The amount of reagent sprayed for each fertilizer application shall comply with the following standards: foliar spray with water droplets on the leaf surface but do not fall.
Determination of physiological indexes and bacterial diversity after treatment with chinese herbal compound soil conditioner After 240 days of plant treatment, measure its height from the base of the plant with a tape measure near the ground; use a caliper to measure the diameter of the stem at the base; The weight of a single tuber is measured using an electronic balance. Three replicates were set for each index of each treatment, and each replicate had 5 plants. Yield per hectare is calculated from the yield per plant.
Quality treatment of cassava tubers: Choose different treatments, moderately sized and mature tubers, rinsed with pure water, dry them, and cut them into small pieces, two samples were taken from the pieces, and repeat 3 pieces for samples each treatment. One sample was put in an oven at 60°C for 3 days, and then its water content was measured. The other was wrapped in tin foil and quickly placed in dry ice to measure the diversity of the bacteria. The rhizosphere soil was taken from the soil 0.5 cm near the tubers, put in a kraft paper envelope and put in dry ice for testing. Each treatment was repeated 3 times.
Root tuber starch content and nutrient composition determination: Use ICPOES to determine the mineral elements in cassava; quickly determine the ash content in cassava starch based on near-infrared spectroscopy.
The experimental data were statistically analyzed using WPS 2019, and the variance and correlation were analyzed using SPSS 20.0.
The diversity of cassava rhizosphere bacterial ora was determined by Shenggong Bioengineering (Shanghai) Co., Ltd.

DNA Extraction and Illumina MiSeq Sequencing
Genomic DNA was extracted from freeze-dried soil sample with E.Z.N.ATM Mag-Bind Soil DNA kits. Quanti cations of bacterial 16S rDNA and fungal ITS rDNA were performed using the Qubit3.0 DNA detection kits. We selected the bacterial V3-V4 regions of the 16S rRNA gene using the individual barcoded primers 341F/805R and fungal 1-2 regions of ITS gene using the individual barcoded primers ITS1F/ITS2R for high-throughput sequencing analysis 27 . The Polymerase Chain Reaction ampli cation conditions were as follows, the reaction mixtures in each tube contained 1-µL of target DNA, 15-µL of 2 ×Taq master Mix (2 ×, TaKaRa, Japan), 1-µL of forward and reverse primers, and 30-µL of sterile distilled water, with a total volume of 48-µL. A two-step PCR approach was as follows: Thermal cycler conditions for primary PCRs consisted of 3 min at 93°C, followed by 5 cycles of 30 s at 94°C, 20 s at 45°C and 30 s at 65°C, followed by 20 cycles of 20 s at 94°C, 20 s at 55°C and 30 s at 72°C, and nal extension it for 5 min at 72°C. Secondary PCRs comprised of denaturation for 3 min at 95°C, followed by 5 cycles of 20 s at 94°C, 20 s at 55°C and 30 s at 72°C, and nal extension it for 5 min at 72°C. After PCR ampli cation, the amplicons were quanti ed. Next, samples were loaded onto an Illumina MiSeq highthroughput sequencing platform for paired-end sequencing 28 , and sequenced by Sangong Bioengineering (Shanghai, China). Illumina MiSeq raw sequences were processed and analysed using the Quantitative Insights into Microbial Ecology (QIIME) software package (version 1.8.0) 29 . After the PCR products of the target fragments of soil microbes were sequenced, the effective sequence was processed for quality control, splicing, and precise removal of impurities to obtain an optimized sequence. The quality-ltered sequences were subsequently clustered in Operational Taxonomic Units (OTUs) at 97% similarity and annotated using the Ribosomal Database Project (RDP) and the Silva bacteria database to determine the phylogeny and relative abundance of the OTUs 30 .

Bioinformatics and Statistical Analysis
The R software (version 4.0.2) was used for bioinformatics analysis, and some plots were generated using the "ggplot2" package. Use the "vegan" package to calculate the number of microorganisms, abundance and diversity index based on the 16S and ITS OTU table. The relative abundances of bacteria were expressed as percentages. T test or two-way analysis of variance (ANOVA) with the Duncan's multiple range test was performed for multiple comparisons to determine the substantial differences of the total number and α-diversity index of bacteria and fungi in the rhizosphere soil, and use Excel 2019 software is visualized. All statistical tests performed in this study were considered signi cant at p < 0.05 using SPSS version 20.0 software. Differences were considered signi cant at P-values less than 0.05. The effects of different dosages Chinese herbal medicine on microbial core and unique OTUs in each soil environment were analyzed according to the method provided by Zhao 31

and Shade and
Handelsman 32 . The OTUs that consistently appeared in the three biological replicates of all soil samples were regarded as core microbiomes, while the some OTUs that only presented in the three biological replicates of one soil samples were considered as unique microbiomes. These results were visualized using the Venn diagrams drawn. Principal coordinates analyses (PCoA) based on Bray-Curtis distance were used to evaluate the differences among microbial communities of different rhizosphere soil. Perform Hierarchical cluatering analysis based on the beta diversity distance matrix. Linear discriminant analysis (LDA) and effect size (LEfSe) analyses were performed using the LEfSe tool 33 . Differences in rhizosphere bacterial abundance were analyzed by LDA EffectSize (LEfSe). LEfSe analysis uses the Kruskal-Wallis rank sum test to detect considerably different abundances and generates LDA scores to estimate the effect size (threshold: ≥2).

Effects of different amendment concentrations on basic physiological indexes
Effects of different amendment concentrations on cassava plant height Tall cassava plants can use light energy more e ciently than shorter plants, thus improving the quality and yield of tuber blocks. Table 1 describes the effects of different concentrations on growth indexes of four groups. As shown in Table 1, there was a substantial difference in plant height of cassava plants treated with different dosage of modi ed preparation. All three treatments showed a considerable increase, among which T2 was highest, with an increase of 24.22% compared with T4. Relative to T4, T1 and T3 increased by 15.70% and 12.56%, respectively; however, there was no signi cant difference between the two groups, as shown in Fig. 1 (P < 0.05). Note: The values in the table are mean ± standard deviation; Different lowercase letters in the same column represent signi cant differences at the P < 0.05 level.

Effects of different amendment concentrations on plant base stem diameter
The stem diameter of cassava can directly re ect the growth potential of cassava. Generally speaking, the taller plants with strong stems have a strong growth potential, and the yield of tuber blocks is also high. As shown in Table 1 and Fig. 1, the diameter of the base stem of cassava was substantially different among plants treated with different dosages of modi ed preparation. All three treatments showed thickening, with T1 and T2 increased signi cantly relative to T4, by 9.68% and 18.82%, respectively (P 0.05). The increase of T3 was 4.61%, and it was not signi cantly different from that of T4.

Effects of different amendment concentrations on yield parameters
Tubers are the main product of cassava, and thus the main way to obtain economic bene ts; we found that soil amendments can substantially improve the quality and yield of cassava tubers. Table 1 and Fig. 2  The maximum tuber block circumference was signi cantly increased in T2, T3 and T4, but T1 was not signi cantly increased (P 0.05). In terms of the number of tuber, the increase of T1 and T3 was not substantial except for T2. The yield per plant in all three groups was signi cantly increased, but there was no signi cant difference between T1 and T3(P 0.05). The yield per mu increased signi cantly, but T1 and T3 did not increase signi cantly.

Effects of different concentration treatments on the quality of cassava tuber blocks
In uence of different concentration treatments on water content of root tubers Starch content is an important indicator for quality evaluation of cassava. Cassava with a high starch content is highly suitable for deep processing industries such as starch processing. Our study showed that the use of soil amendments changed the starch content, and the starch content gradually decreased with the increase of the application concentration, but there was no substantial difference between groups on the whole. The water content of treated cassava root tuber varied by treatment (Table 2). Relative to T4, a decrease in water content was shown in T2 and T3, of which T2 was greatest, with a decrease of 5.96%. In T1, the water content increased by 1.1%. Moreover, the dry matter content of T2 also increased substantially (an increase by 5.95% relative to T4).

Effects of different concentration treatments on the nutrient components of root tubers
The mineral content of cassava directly indicates the tapioca edible nutritional value, which is crucial to inform the breeding of good varieties in the future. Under the action of soil amendment, the mineral content of cassava showed a rising trend overall, including K as one of the highest levels of mineral elements in cassava showed increased substantially.   (Table 5), with a total of 10 distinct fungal phyla and 528 genus were detected across all samples.  The Shannon index curves indicated that a su cient amount and quantity of sequencing data was obtained ( Supplementary Fig. S1 and S2).

As shown in
Bacterial richness (abundance-based coverage estimator, ACE) and Chao index were used to estimate the abundance and diversity of microorganisms. Through the analyses of alpha diversity (Chao, Ace) ( Supplementary Fig. S3, S4, S5 and S6), we found that the diversity of bacteria in the soil samples was as following: T2 > T1 > T3 > T4, while the fungal diversity was dominated by T1, and followed by T2, T3, T4.
Compared with the control group, the experimental group had higher microbial diversity, and the soil samples treated with T1 concentration had the highest fungi richness. At the same time, the diversity of bacteria in T2 was the most abundant.

The composition of bacterial and fungal communities
The top ten bacterial phyla in the tested samples were Cyanobacteria_Chloroplast Proteobacteria Actinobacteria Acidobacteria Chloro exi  (Fig. 3). And the top ten fungi phyla in the tested samples were Ascomycota,Basidiomycota Glomeromycota Mortierellomycota Mucoromycota Chytridiomycota Zoopagomycota Kickxellomycota Rozellomycota and Blastocladiomycota (Fig. 4). Among them, the dominant fungi phyla in all soil samples were Ascomycota and Basidiomycota, accounting for 80.3% and 12.9%, respectively. (Relative abundance > 10%).
α Under different concentration treatment, at the genus level, the heatmap indicated that Gp1 was the most abundant bacteria in the soil, followed by unclass ed_Ktedonobacterales and Gp3. Among them, the relative abundance of Streptophyta was substantially higher in T2 than that of other concentration treatments, accounting for about 7.3% of the identi ed sequences. As for fungi, the composition differed greatly between the groups after treatment with different concentrations, and the samples under different concentration treatments contained more than 20 species of peculiar fungi. In addition, at the fungal genus level, the relative abundance of unclass ed_Chaetomiaceae in the treatment group was considerably lower than that of T4, and the relative abundance of Psathyrella in the T1 increase substantially compared to other treatment groups (Fig. 5,6).
Venns show us the unique bacterial genus of T2 was Fimbriimonas, the unique bacterial genus of T3 was Candidatus_Procabacter and Propionibacterium, and T1 had no unique bacterial genus. The composition of soil fungal genus differed greatly between the groups after treatment with different concentrations, and the samples under different concentration treatments contained more than 20 species of peculiar fungi. (Fig. 7,8).
At the genus level, we performed a statistical analysis on the soil bacteria and fungi communities of different treatments, using the default logarithm (LDA) value of 2 to identify a total of 11 different bacterial groups and 16 unique fungal species (Supplementary Fig. S7, S8). The The community structure of bacterial and fungal communities The PCoA map based on the Bray-Curtis distance ( Supplementary Fig. S9, S10) showed that bacteria and fungi communities from treated soil samples (T1, T2, T3) are clearly separated from the control group (T4), respectively, explaining the 47.31% and 36.95% of the total variation.

Discussion
Growth-promoting effect of the amendment on cassava Cassava has been recognized as the best biomass energy crop in tropical and subtropical regions of China, and is one of the key crops for the development of biomass energy strategy in the "11th Five-Year Plan." It is an important goal of cassava cultivation to improve plant growth potential, root tuber yield and starch content 34,35 . In this paper, we found that amendment with Chinese herbal compound amendment powder, bacteriostatic agent and compound soil antivirals have a positive effect on the cassava plant growth. Various physiological indexes increased after using amendment, relative to T4, particularly in plant height, stem diameter, maximum circumference, plant tuber number and area yield. The physiological indexes increased may because the soil amendment powder can effectively improve soil structure, enrich soil nutrients and increase soil microbial activity, thus promoting the growth and yield of cassava plants 36 . In addition, some plant rhizosphere growth-promoting microorganisms in the soil also release plant hormones, such as indole acetic acid (IAA), abscisic acid (ABA), cytokinin (CTK) and gibberellin (GA), and some strains also promote plant growth by regulating ethylene levels 37,38,39,40,41 .
The root moisture content of cassava also substantially changed after the application of amendment powder and water. T2 was the lowest among the three treatments, while T1 showed a relative increase. In terms of root starch content, there was a decrease in all three treatments with the increase of the concentration of the improver, and the lowest value was reached when the application concentration was 1500 kg/hm 2 . The decrease of water content in root tubers may have occurred because soil amendment powder and water agent changed the water and fertilizer utilization rate of the plant, improved the absorption and conversion e ciency of soil nutrients, and thus increased the transformation and accumulation of dry matter in the plant 10 .

Effects of the amendment on nutrient indexes of cassava root tubers
The contents of most mineral elements, such as K, Ca, Na, Mg and Mn, changed following treatment. According to the above experimental results, T3 had the best effect in terms of the content of root tuber nutrients. This was probably because soil amendment powder contains a high ionic charge and can thus carry out ion exchange while increasing soil fertility, improve soil permeability to a certain extent and nally enhance nutrient accumulation. This supports the results of Yue 42 .
Moreover, spraying Chinese herbal compound bacteriostatic agent and antiviral water agent may change the structure of soil microbial community around the plant, such as bene cial bacteria and harmful bacteria. Studies have also shown that microorganisms determine soil nutrients, and increasing numbers of scientists have begun to apply microbial activity as an index to measure soil fertility 43 . For example, lactic acid bacteria and yeast can promote the transformation of nutrients in soil, enabling them to be absorbed by plants and participate in the nutrient cycle 44 .

Effects Of The Amendment On Microbial Diversity Of Cassava
Through experiments, we found that the cassava plants in the T2 treatment group had great improvements in both phenotype and physiological indicators compared with other treatments. We further analyzed the bacterial abundance of the soil samples. It was found that in each treatment group Streptophyta is the dominant bacterial genera in T2 soil samples. Streptophyta was greatly important in forming biological soil crust,. Biological crusts played an essential components in forming arid and semi-arid ecosystems as they x nitrogen and carbon, and enhance soil fertility 45 .
Fimbriimonas is a unique genera of T2, and the abundances of Nocardioides and unclass ed_Spartobacteria in the T2 group were greatly different from those of other treatment groups. Shen et al. carried out redundancy analysis and Spearman correlation analysis on bananas and showed that Spartobacteria and other bacteria are dominant in bananas treated with high-concentration biological fertilizers, and it was positively correlated with soil pH, total nitrogen, total carbon and available phosphorus content. It was negatively correlated with the incidence of banana Fusarium wilt 46 . Chen found that the abundance of Nocardioides was positively correlated with the biomass accumulation and available phosphorus content of winter wheat during soil experiments on winter wheat 47 . At the same time, Fimbriimonas also had the characteristics of nitrogen xation and denitri cation, and played an important role in maintaining nitrogen balance 48 . All this showsed that the dominant or unique bacteria in T2 can affect the content of nutrient elements in the soil to regulate plant growth. In addition, when Liu studied the soil microbial mechanism of resistance to root rot in sugar beet varieties, they also found that the relative abundance of Nocardioides in samples with less disease was substantially increased compared with samples with more severe disease 49 . This may be because the presence of Nocardioides can enhance the ability of plants to resist diseases.
The author speculates that these advantages and the existence of unique bacteria can affect rhizosphere soil microorganisms, further improve soil effective nutrients and enzyme activities, and increase plant stress resistance, thereby affecting the growth and physiological indicators of cassava plants.

Necessity of cassava topdressing
Cassava has a strong ability to absorb water and fertilizer, and scienti cally informed fertilization in continuous cropping will lead to a decrease in soil fertility and the eld production of cassava 50 , which is often considered to be a strong fertilizer consumption 51 . The decrease in soil fertility leads to the unsustainable productivity of the cassava agro-ecosystem 36 . However, cultivation of casava doesn't need to apply base fertilizer.
Putthacharoen suggested that fertilization does not need to be applied to cassava in the rst two years of continuous cropping in newly reclaimed land 52 . However, fertilization should be applied to medium fertility soil or two years after continuous cropping. In addition to the application of chemical fertilizer as base fertilizer, Luo showed that the application of 600-750 kg/hm 2 bio-organic fertilizer increased the yield of fresh tubers by 57.3-74.8% compared with an unfertilized treatment, and promoted the growth of cassava stem and leaves, the transport and accumulation of assimilates to the underground part, and improved the yield of root tuber and the starch content 53,54 . However, excessive fertilization of cassava should be avoided, as the fertilizer utilization rate is low and excessive fertilization will lead to enhanced growth of the aboveground part and is not conducive to the formation of root tuber yield 55 . Studies have shown that the tuber starch content presents a parabolic change with the increase of nitrogen application rate 56 . Therefore, topdressing in the cultivation should be paid more attention to meet the nutrient requirements of cassava root tuber in the formation and expansion stages. Therefore, the use of Chinese herbal compound soil amendment powder as topdressing fertilizer could help to meet the demand of cassava growth economically and e ciently.
Our ndings indicate that 1200 kg/hm 2 treatment at the intermediate level performs best, and staged fertilization is conducive to reducing fertilizer waste and improving recovery e ciency, which is consistent with the research results of Luo 55 and Sangakkara 57 . Zhang showed that the application amount of N, P and K had no signi cant effect on the starch content of the tubers 58 . In the current study, the traditional Chinese herbal soil amendment powder used included rice bran as the main source of organic nutrients, and Artemisia annua, A. argyi and other medicinal materials, which are safe and non-toxic. It can also be used in combination with reduced chemical fertilizer or organic fertilizer, which can signi cantly improve crop yield quality and increase the cassava Se content. This is consistent with the Guizhou research team's results for tomato, strawberry, tea, R. pseudostellariae and Pinellia ternata crops from 2018 to 2020. Under the condition that most chemical agents were replaced by Chinese herbal compound bacteriostasis agent and antiviral water agent by spraying or irrigation root, 900-1200 kg/hm 2 Chinese herbal compound soil amendment powder per mu could signi cantly improve crop yield and quality.   Histogram of the relative abundances at the phyla level of the TOP 10 fungal communities

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download. SupplementaryFig.docx