Effect of Biochar on Rhizosphere Soil Microbial Diversity and Metabolism in Tobacco-Growing Soil

: In this study, four different biochar application rates and a control were set up using indoor potted tobacco, to study the effects of biochar on the microbial diversity and metabolism of tobacco-growing soil. The ﬁve treatments were as follows: control—0% biochar ( w / w ) + 26 g fertilizer/pot; biochar treatments—1% biochar ( w / w ) + 26 g fertilizer/pot, 2% biochar ( w / w ) + 26 g fertilizer/pot, 3% biochar ( w / w ) + 26 g fertilizer/pot, and 4% biochar ( w / w ) + 26 g fertilizer/pot. We found that biochar increases the microbial diversity of soils and simultaneously changes the microbial community structure. Under the inﬂuence of biochar, soil urease activity increased by 18%, invertase activity increased by 23.40%, polyphenol oxidase activity increased by 59.50%, and catalase activity increased by 30.92%. Biochar also signiﬁcantly increased the microbial biomass carbon and nitrogen content of the soil. Soil microbial biomass nitrogen had a positive correlation on bacterial diversity, with the highest coefﬁcient, while soil microbial biomass carbon had a positive correlation on fungal diversity, with the highest coefﬁcient. The microbial diversity and metabolic capacity of soil are improved under the inﬂuence of biochar, and soil enzyme activity and microbial biomass carbon and nitrogen have positive impacts on soil microbial diversity.


Introduction
Tobacco is widely cultivated in China. However, the problems faced by tobacco planting are the same as those faced by China's other agricultural products, and the degradation of tobacco planting soil has become more and more serious, leading to unhealthy development of tobacco plants and serious diseases and pests. These factors have led to a great loss of tobacco farmers. In view of the serious problems existing in tobacco planting soil, an effective material and application method is needed to improve the soil, improve the growth of tobacco plants, and ensure the income of tobacco farmers.
Biochar is a material obtained by pyrolysis of crop straw, wood materials, livestock manure or other organic materials at a low temperature (300-700 • C) and low oxygen environment, and it is a carbon rich material with rich surface functional groups, high porosity, specific surface area, and strong adsorption [1]. As a new type of high carbon organic material, biochar had been favored by many research fields in recent years. Biochar, known as "black gold", is attracting increasing attention in upcoming fields, while the relevant scientific research reports have also been increasing year by year. This is mainly due to the huge potential of biochar in soil improvement, ecological restoration, and rational utilization of waste [2].
In the past two decades, an increasing number of studies have focused on the effects and mechanisms of biochar on soil rhizosphere microorganisms. The German microbiologist Lorenz Hiltner defined the rhizosphere in 1904 as the soil around the root system that is affected by the growth of the root system. Some rhizosphere microorganisms are

Materials
An indoor pot experiment was carried out in the Modern Tobacco Science and Technology Park of Henan Agricultural University, Xuchang City, Henan Province. In the large field of the Science and Technology Park, five bulk soil samples (0-20 cm soil layer) were randomly selected from an area of approximately 335 m 2 . The soil was mixed, and the samples were naturally air-dried in the shade at room temperature. Then, the soil was passed through a 2-mm sieve before use. The basic information of the 0-20 cm soil layer was as follows: soil type, cinnamon; organic matter content, 1.02%; pH, 6.4; EC, 2.6 mS/cm; total nitrogen, 0.713 g kg −1 ; total carbon, 7.772 g kg −1 ; total sulfur, 0.627 mg kg −1 ; pH, 6.5; available nitrogen, 0.098 g kg −1 ; available phosphorus, 0.026 g kg −1 ; and available potassium, 0.207 g kg −1 . The continuous flow carbonization method described by Nielsen et al. [22] was used to produce biochar from peanut shells (Wang et al. [23] thought that the biochar made from peanut shells had a better effect on soil improvement). Briefly, the pyrolysis temperature was increased to 450 • C at a rate of 26 • C min −1 and held for 30 min. The basic physical and chemical properties of the obtained biochar were as follows: specific surface area, 16.72 m 2 g −1 ; pH 8.6; conductivity, 2.3 ds m −1 ; TN, 7.9 g kg −1 ; total carbon, 773 g kg −1 ; total potassium, 4.1 g kg −1 ; total phosphorus, 2.3 g kg −1 ; and total sulfur, 0.39 g kg −1 . The tested tobacco variety was 'Zhongyan 100' (CF965, a flue-cured tobacco variety).

Experimental Design
Approximately 20 kg of dry sieved soil (soil that passed through a 2-mm sieve) was placed into each pot (n = 15, diameter 0.45 m × height 0.40 m. Each pot holds 25 kg of soil). The experimental design included five experimental treatments: the control treatment [0% biochar (w/w) + 26 g fertilizer/pot (T0)] and four biochar treatments [1% biochar (w/w) + 26 g fertilizer/pot (T1); 2% biochar (w/w) + 26 g fertilizer/pot (T2); 3% biochar (w/w) + 26 g fertilizer/pot (T3); and 4% biochar (w/w) + 26 g fertilizer/pot (T4)]. The weight of biochar applied per plant of each treatment was as follows: 0.00 kg biochar/pot (T0); 0.25 kg biochar/pot (T1); 0.5 kg biochar/pot (T2); 0.75 kg biochar/pot (T3); 1.00 kg biochar/pot (T4). Under normal circumstances, 15,000 tobacco plants are planted per hectare of tobacco field, so the amount of biochar used in the pot experiment was converted into the amount used in tobacco farmland: 0.00 kg biochar/hm 2 (T0); 3750 kg biochar/hm 2 (T1); 7500 kg biochar/hm 2 (T2); 11,250 kg biochar/hm 2 (T3); 15,000 kg biochar/hm 2 (T4). Fertilizer ratio: m(K 2 SO 4 ): m[Ca(H 2 PO 4 ) 2 ]:m(NH 4 NO 3 ) =2:1:1, which is consistent with tobacco fertilization recommendations [24]. The biochar and fertilizer were fully mixed with the soil in the pot. Note that the volume differences caused by the addition of biochar and fertilizer to the soil were negligible. Each treatment was repeated three times, and each repetition included 15 pots. The biochar and fertilizer were applied 4 days before transplanting, and one tobacco seedling was transplanted into each pot. The distance between the pots in the greenhouse was 0.30 m, and the experimental layout was a random block design. During transplanting, 1000 mL of water was added per pot. The pots were then watered with 200 mL every 5 days. Forty days after transplanting, the watering frequency was increased to every 3 days. After sampling, the unsampled tobacco plants continued to grow until they died naturally. The dead tobacco plants (including the roots) were then removed from the soil. During this process, little soil was lost from the pot. The pots were watered with 50 mL every week until tobacco was planted again the following year. Sixty days after tobacco transplantation, three uniformly growing and representative tobacco plants were selected from each treatment. The tobacco roots were carefully removed from the pot, the soil attached to the root system gently shaken off, and a brush applied to the lateral roots. The soil around the root system was gently brushed off, to collect the soil within 2 mm of the root surface (rhizosphere soil). The rhizosphere soil sample was then passed through a 2-mm sieve. Ten grams of rhizosphere soil was taken from each plant, placed in a 10-mL centrifuge tube wrapped in aluminum foil, and stored in liquid nitrogen for later use.

Soil Physical and Chemical Indicators
Two soil samples, each 500 g, were collected from the pots after the rhizosphere soil samples had been collected. One sample was naturally air-dried in a cool environment and the other was sealed and stored at 4 • C. Soil enzyme activity was measured using an enzyme activity kit produced, and the soil MBC and MBN were measured using the chloroform fumigation-potassium sulfate method [25]. The Elementar Vario MAX CN analyzer was used to measure the total organic carbon (TOC) and TN content of the biochar, using the Dumas combustion method.

Measurement Methods Soil Bacterial DNA Extraction and PCR Amplification
When extracting the total DNA of the soil microorganisms following the instructions of the EZNA ® soil DNA kit (Omega Bio-tek, Norcross, GA, USA), 1% agarose gel electrophoresis was used to detect the DNA extraction quality. Purity and concentration were determined using a NanoDrop2000 system. Then, 806R (5 -GGACTACHVGGGTWTCTAAT-3 ) and 338F (5 -ACTCCTACGGGAGGCAGCAG-3 ) were used to perform PCR amplification on the 16SrRNA V3-V4 variable region, while SSU0817F 5 -TTAGCATGGAATAATRRAA TAGGA-3 and 1196R 5 -TCTGGACCTGGTGAGTTTCC-3 were used for amplification of the 18SrRNA V5-V7 variable region. The amplification procedure was as follows: 95 • C pre-denaturation for 3 min; 27 cycles of 95 • C denaturation for 30 s, annealing at 55 • C for 30 s, and extension at 72 • C for 30 s; 72 • C stable extension for 10 min, and finally storage at 4 • C (PCR instrument: ABI Gene Amp ® Type 9700). The PCR reaction system consisted of 5 × Trans Start Fast Pfu buffer 4 µL, 2.5 mM dNTPs 2 µL, downstream primer (5 µM) 0.8 µL, upstream primer (5 µM) 0.8 µL, Trans Start Fast Pfu DNA polymerase 0.4 µL, template DNA 10 ng, and supplement to 20 µL. Three replicates were used for each experiment.

Illumina Miseq Sequencing of Soil Microorganisms
After mixing the PCR products from the same samples, a 2% agarose gel was used to recover the PCR products. An Axy Prep DNAGel Extraction Kit (Axygen Biosciences, Union City, CA, USA) was used to purify the recovered products, which were then detected by 2% agarose gel electrophoresis. A Quantus™ Fluorometer (Promega, Madison, WI, USA) was used to detect and quantify the recovered products. We used the NEXTFLEX Rapid DNA-Seq Kit to build the library as follows: (1) the linker was added; (2) magnetic beads were used to screen and remove the self-linked fragments of the linker, (3) PCR amplification was used to enrich the library template; and (4) the magnetic beads were recovered to obtain the final library. Sequencing was performed using the Illumina MiSeq PE300 platform (Shanghai Meiji Biomedical Technology Co., Ltd., Shanghai, China).

Soil Microbial Data Processing
We used the original sequence of the Trimmomatic software for quality control and FLASH software for splicing.
(1) Bases with tail quality values less than 20 were filtered out and a 50-bp window set. If the average quality value in the window was less than 20, the bases were removed from the beginning of the window. Filtered reads of less than 50-bp after quality control were removed, as were reads containing N bases.
(2) The paired reads were merged into a sequence according to the overlap relationships between PE reads, while the minimum overlap length was 10 bp.
(3) The maximum allowable mismatch ratio in the overlap region of the spliced sequence was 0.2, and the unmatched sequence was screened.
(4) The samples were differentiated according to the barcodes and primers at the beginning and end of the sequence, and the sequence direction was adjusted. The allowed number of mismatches in the barcode was zero, and the maximum number of primer mismatches was 2.
UPARSE software (version 7.1, http://drive5.com/uparse/, accessed on 25 May 2021) was used to perform OTU clustering of sequences based on 97% similarity, and to eliminate chimeras. The RDP Classifier (http://rdp.cme.msu.edu/, accessed on 28 May 2021) was used to classify and annotate each sequence and compare it to the Silva database (SSU132) with a comparison threshold of 70%. The statistics of the bacterial and fungal samples are shown in Tables 1 and 2.

Data Analysis Methods
Chao index: This index uses the chao1 algorithm to estimate the number of OTUs contained in the sample. The calculation formula used in this analysis is as follows: where S chao1 = estimated number of OTUs; S obs = observed number of OTUs; n 1 = number of OTUs containing only one sequence ("singletons"); n 2 = number of OTUs containing only two sequences ("doubletons"). ACE index: an index used to estimate the number of OTUs in a community. The calculation formula used in this analysis is as follows: , 0 n i = the number of OTUs containing i sequences; S mall = the number of OTUs containing "abund" or less than "abund"; S abund = the number of OTUs with more than "abund" sequences; abund = the threshold of "advantage" OTUs. The default is 10.
Shannon index: One of the indices used to estimate the diversity of microorganisms in a sample. The calculation formula used in this analysis is as follows: where, S obs is the number of OTUs actually observed, n i is the number of sequences contained in the i-th OTU, and N is the total number of sequences. The coverage index refers to the coverage rate of each sample library. The higher the value, the higher the probability that each sample sequence will be detected, and the lower the probability of it not being detected. The calculation formula used in this analysis is as follows: C = 1 − n 1 N where n 2 = the number of OTUs containing only one sequence and N = the total number of sequences in the sample.
Histograms and line charts were drawn in Excel 2013. One-way ANOVA was used to analyze the different treatments using SPSS 20.0. One-way ANOVA and Duncan's method were used to compare the differences between the different treatments. Mothur (v.1.30.1, http://www.mothur.org/wiki/Schloss_SOP#Alpha_diversity, accessed on 15 June 2021) software was used for index analysis, and the OTU similarity level used for index evaluation was 97% (0.97).

Species Assessment of Bacteria and Fungi
The dilution curves of bacteria and fungi tended to be flat ( Figure 1); the coverage of bacteria reached 98.12% and the coverage of fungi reached 99.97% ( Figure 1). These results indicate that the detection rate of the soil sample microbial community was close to saturation, and the current sequencing volume can cover most of the species in the sample.

Alpha Diversity of Bacteria and Fungi
T1, T2, T3, and T4 were the treatment groups, and T0 was the control group (Table 3). The soil bacterial diversity Chao index values of the treatment groups were significantly higher than that of the control group. The T1, T2, T3, and T4 values were 69.61%, 63.42%, 58.18%, and 64.32% higher than the T0 value, respectively. The fungal diversity Chao index values of the treatment groups were also significantly higher than that of the control group. The T1, T2, T3, and T4 values were 13.51%, 15.37%, 8.29%, and 7.13% higher than the T0 value, respectively. The bacterial diversity ACE indices of T1, T2, T3, and T4 were all significantly higher than that of T0 (by 7.58%, 7.55%, 6.17%, and 7.33%, respectively). The fungal diversity ACE index of T2 was significantly higher than that of T0 (by 12.79%). There was no significant difference in bacterial Shannon for each treatment, and this index was the highest in T4; fungal Shannon values in T3 and T4 were higher than that of T0 (by 3.65%, 3.32%). The bacterial diversity coverage indices of T0 and T2 were significantly higher than those of T3 and T4, while the fungal diversity coverage indices of the different treatments were not significantly different (Table 3).

Species Assessment of Bacteria and Fungi
The dilution curves of bacteria and fungi tended to be flat ( Figure 1); the coverage of bacteria reached 98.12% and the coverage of fungi reached 99.97% ( Figure 1). These results indicate that the detection rate of the soil sample microbial community was close to saturation, and the current sequencing volume can cover most of the species in the sample.

Alpha Diversity of Bacteria and Fungi
T1, T2, T3, and T4 were the treatment groups, and T0 was the control group ( Table  3). The soil bacterial diversity Chao index values of the treatment groups were significantly higher than that of the control group. The T1, T2, T3, and T4 values were 69.61%, 63.42%, 58.18%, and 64.32% higher than the T0 value, respectively. The fungal diversity Chao index values of the treatment groups were also significantly higher than that of the control group. The T1, T2, T3, and T4 values were 13.51%, 15.37%, 8.29%, and 7.13% higher than the T0 value, respectively. The bacterial diversity ACE indices of T1, T2, T3, and T4 were all significantly higher than that of T0 (by 7.58%, 7.55%, 6.17%, and 7.33%, respectively). The fungal diversity ACE index of T2 was significantly higher than that of T0 (by 12.79%). There was no significant difference in bacterial Shannon for each treatment, and this index was the highest in T4; fungal Shannon values in T3 and T4 were higher than that of T0 (by 3.65%, 3.32%). The bacterial diversity coverage indices of T0 and T2 were significantly higher than those of T3 and T4, while the fungal diversity coverage indices of the different treatments were not significantly different (Table 3). Table 3. Analysis of diversity indices among different treatments. The data are means ± standard deviation. Different lowercase letters in the same row indicate significant differences between treatments (p < 0.05).

Beta Diversity of Bacteria and Fungi
Sample-level cluster analysis was performed on the community distance matrices of the bacteria (Figure 2A) and fungi ( Figure 2B). The results of sample hierarchical clustering show that the samples can be divided into four significantly different groups according to their bacteria, indicating that the bacterial community compositions of the treatment groups and the control were significantly different (Figure 2A). According to the fungi, the samples could be divided into three significantly different groups, indicating that the fungal community compositions of the treatment groups and the control group were significantly different ( Figure 2B).

Beta Diversity of Bacteria and Fungi
Sample-level cluster analysis was performed on the community distance matrices of the bacteria (Figure 2A) and fungi ( Figure 2B). The results of sample hierarchical clustering show that the samples can be divided into four significantly different groups according to their bacteria, indicating that the bacterial community compositions of the treatment groups and the control were significantly different (Figure 2A). According to the fungi, the samples could be divided into three significantly different groups, indicating that the fungal community compositions of the treatment groups and the control group were significantly different ( Figure 2B).

The Impact of Biochar on Soil Enzyme Activity
Biochar improved the activities of invertase, urease, polyphenol oxidase, and catalase. With increasing biochar, the soil urease activity showed a trend of first increasing and then decreasing. The treatment group values were significantly higher than that of the control group, and the differences were significant. T3 showed the highest urease activity (18.06% higher than that of T0) ( Figure 3A). The invertase activity values of the treatment groups were higher than that of the control group, and the differences were significant. As the amount of biochar increased, the soil invertase activity showed a trend of first increasing and then decreasing. The invertase activity in T4 was 23.40% higher than that in T0 ( Figure 3B). The change rule of soil polyphenol oxidase activity was similar to that of invertase, with T4 59.50% higher than T0 ( Figure 3C). The catalase activity values in the treatment groups were higher than that in the control group, and the difference was significant. The catalase activity of the T3 treatment was the highest, 30.92% higher than that of T0 ( Figure 3D).

The Impact of Biochar on Soil Microbial Biomass Carbon and Microbial Biomass Nitrogen
Biochar can significantly increase the content of MBC and MBN in the soil. Although the soil TOC and TN content also increased with biochar input, MBC/TOC and MBN/TN increased with biochar input ( Table 4). The soil MBC, MBN, MBC/TOC, and MBN/TN all showed trends of first increasing and then decreasing with increasing biochar amount.
The MBC values of T3 and T4 were significantly higher than those of the other treatments; they were 40.90% and 38.69% higher than those of T0, respectively. The MBN values of T2 and T3 were significantly higher than those of the other treatments; they were 42.66% and 32.14% higher than those of T0, respectively. The MBC/TOC of T3 was the highest, significantly higher than those of the other treatments, while the MBN/TN of T2 was the highest, also significantly higher than those of the other treatments (Table 4). ment groups were higher than that of the control group, and the differences were significant. As the amount of biochar increased, the soil invertase activity showed a trend of first increasing and then decreasing. The invertase activity in T4 was 23.40% higher than that in T0 ( Figure 3B). The change rule of soil polyphenol oxidase activity was similar to that of invertase, with T4 59.50% higher than T0( Figure 3C). The catalase activity values in the treatment groups were higher than that in the control group, and the difference was significant. The catalase activity of the T3 treatment was the highest, 30.92% higher than that of T0( Figure 3D).

The Impact of Biochar on Soil Microbial Biomass Carbon and Microbial Biomass Nitrogen
Biochar can significantly increase the content of MBC and MBN in the soil. Although the soil TOC and TN content also increased with biochar input, MBC/TOC and MBN/TN increased with biochar input ( Table 4). The soil MBC, MBN, MBC/TOC, and MBN/TN all showed trends of first increasing and then decreasing with increasing biochar amount. The MBC values of T3 and T4 were significantly higher than those of the other treatments;  Table 4. The impact of biochar on soil microbial biomass carbon (MBC), microbial biomass nitrogen (MBN), total organic carbon (TOC), and total nitrogen (TN). Different letters in the same column indicate significant differences between treatments (p < 0.05).

Correlation Analysis of Soil Enzyme Activity and Microbial Diversity
Soil urease, invertase, polyphenol oxidase, and catalase were positively correlated with soil microbial diversity (Figures 4 and 5). The correlation coefficients (R 2 ) between catalase and the diversities of bacteria and fungi were the highest, with scores of 0.2353 and 0.2296, respectively; the second highest R 2 values were for invertase, 0.1852 and 0.142, respectively. The correlation coefficients for the relationships between urease and the soil bacteria and fungi were 0.045 and 0.0796, respectively; while those for polyphenol oxidase were 0.1107 and 0.0852, respectively (Figures 4 and 5).
with soil microbial diversity (Figures 4 and 5). The correlation coefficients (R 2 ) between catalase and the diversities of bacteria and fungi were the highest, with scores of 0.2353 and 0.2296, respectively; the second highest R 2 values were for invertase, 0.1852 and 0.142, respectively. The correlation coefficients for the relationships between urease and the soil bacteria and fungi were 0.045 and 0.0796, respectively; while those for polyphenol oxidase were 0.1107 and 0.0852, respectively (Figures 4 and 5). . Regression analysis of bacterial ranking. The X-axis is the environmental factor, the Y-axis is the α-diversity ranking axis, and R 2 is the coefficient of determination, which represents the proportion of variation explained by the regression line. The larger the R 2 value, which indicates that the environmental factor shows differences in the community composition or α-diversity of the sample on the ranking axis, the higher is the degree of explanation of the difference in the diversity index. (A) The environmental factor is urease (UR); (B) the environmental factor is sucrase (SU); (C) the environmental factor is polyphenol oxidase (PO); (D) the environmental factor is catalase (CA).

Correlation Analysis of Soil Microbial Biomass Carbon and Nitrogen and Microbial Diversity
Soil MBC, MBN, TOC, TN, MBC/TOC, and MBN/TN had positive effects on soil microbial diversity (Figures 6 and 7). The coefficient of influence (R 2 ) of MBN on bacterial diversity was the highest, 0.3553; the second highest was that of TOC, 0.3553. The coefficients of influence of TN, MBC, MBC/TOC, and MBN/TN on soil bacteria were 0.2568, 0.2341, 0.0156, and 0.2132, respectively ( Figure 6). MBC, MBN, TOC, TN, MBC/TOC, and MBN/TN also had positive impacts on soil fungal diversity, but overall, the above indicators had a lower impact on soil fungal diversity than on soil bacterial diversity ( Figure 6). The coefficient of influence of MBC on fungal diversity was the highest, 0.18; the second highest was that of MBN/TN, 0.0958. The coefficients of influence of TOC, TN, MBN, and MBC/TOC on soil fungi (R 2 ) were 0.0775, 0.0824, 0.0285, and 0.0428, respectively (Figure 7). . Regression analysis of bacterial ranking. The X-axis is the environmental factor, the Y-axis is the α-diversity ranking axis, and R 2 is the coefficient of determination, which represents the proportion of variation explained by the regression line. The larger the R 2 value, which indicates that the environmental factor shows differences in the community composition or α-diversity of the sample on the ranking axis, the higher is the degree of explanation of the difference in the diversity index. (A) The environmental factor is urease (UR); (B) the environmental factor is sucrase (SU); (C) the environmental factor is polyphenol oxidase (PO); (D) the environmental factor is catalase (CA).

Correlation Analysis of Soil Microbial Biomass Carbon and Nitrogen and Microbial Diversity
Soil MBC, MBN, TOC, TN, MBC/TOC, and MBN/TN had positive effects on soil microbial diversity (Figures 6 and 7). The coefficient of influence (R 2 ) of MBN on bacterial diversity was the highest, 0.3553; the second highest was that of TOC, 0.3553. The coefficients of influence of TN, MBC, MBC/TOC, and MBN/TN on soil bacteria were 0.2568, 0.2341, 0.0156, and 0.2132, respectively ( Figure 6). MBC, MBN, TOC, TN, MBC/TOC, and MBN/TN also had positive impacts on soil fungal diversity, but overall, the above indicators had a lower impact on soil fungal diversity than on soil bacterial diversity ( Figure 6). The coefficient of influence of MBC on fungal diversity was the highest, 0.18; the second highest was that of MBN/TN, 0.0958. The coefficients of influence of TOC, TN, MBN, and MBC/TOC on soil fungi (R 2 ) were 0.0775, 0.0824, 0.0285, and 0.0428, respectively ( Figure  7).  is the α-diversity ranking axis, R is the coefficient of determination, which represents the proportion of variation explained by the regression line. The larger the R 2 , which indicates that the environmental factor shows differences in the community composition or α-diversity of the sample on the ranking axis, the higher is the degree of explanation of the difference in the diversity index. (A) The environmental factor is total organic carbon (TOC); (B) the environmental factor is microbial biomass carbon (MBC); (C) the environmental factor is TN (total nitrogen); (D) the environmental factor is microbial biomass nitrogen (MBN). (E) the environmental factor is MBC/TOC. (F) the environmental factor is MBN/TN. The α-diversity of soil bacteria refers to the abundance of bacteria in a local area or in the same habitat. It generally includes the Shannon index, Chao index, and ACE index. The results of this study show that biochar promotes an increase in soil bacterial α-diversity. The Shannon and Chao index values of the bacteria in T1 and T2 were significantly different from those in the control group, while the ACE index of T3 showed a highly significant difference from the control (Table 3). It is possible that, as biochar was imported into the soil, the physical and chemical properties of the soil changed, especially the soil carbon and nitrogen levels and the physical structure. This would lead to changes in soil bacterial abundance. An increase in the soil TOC content can increase the abundance of bacterial species in the soil, as organic carbon is the most important source of soil bacterial nutrition [26]. Therefore, with the increase in TOC content, the diversity index increased accordingly. This may also be related to the symbiotic relationship between soil rhizosphere microorganisms and plant roots [27]; biochar can stimulate plant roots to secrete more metabolites, and these root exudates can promote vigorous bacterial growth and reproduction, and enhance bacterial metabolism. In this study, the treatment with 1% (w/w) of biochar had a significant promotion effect on bacterial α-diversity. Although the amount of biochar used in other treatments was higher, bacterial α-diversity decreased. This may be because 1% (w/w) of biochar is more suitable for the improvement of bacterial α-diversity. Excessive biomass will cause more drastic changes in soil physical and chemical properties, which is not conducive to the improvement of bacteria α-diversity.

Bacterial Beta Diversity
This study showed that the bacterial species compositions of each treatment were significantly different, while the differences in the soil bacterial species within the biochar treatments were not obvious, and their similarity with the control was low (Figure 1). It is possible that biochar changes the soil environment and nutritional status and promotes the development of bacterial communities in a specific direction. Changes in the soil environment and its nutritional status are mainly caused by changes in organic carbon, moisture, and gas permeability [28,29]. Biochar can cause changes in all these indicators, subsequently changing the soil environment and guiding the bacterial population to develop in a certain direction. This was also confirmed by the fact that the different biochar treatments led to different changes in their soil bacterial population structures. However, different amounts of biochar have different effects on the composition of bacterial species, indicating that the appropriate amount of biochar has more obvious effects on the change of bacterial species composition. In this study, the treatment of 2% (w/w) biochar had a significant effect on improving bacterial species compositions. This may be because when the biochar reaches the most appropriate concentration of soil microorganisms, the excess biochar will affect the physical and chemical properties of the soil, especially the pH of the soil which will change greatly [28], resulting in a certain limit to the growth of soil microorganisms.

Fungal Alpha Diversity
In the study of fungal α-diversity in soil community ecology, the analysis of the microbial diversity in a single sample can reflect changes in the abundance of soil microbial communities, including the analysis of different types of diversity indices [29]. Soil microbial α-diversity is the microbial diversity in a uniform habitat in a local area; therefore, it is also called the microbial diversity within the habitat. The results of this study show that biochar improved the α-diversity index of soil fungi. The fungal α-diversity index first increased and then decreased with increasing biochar. Among the treatments, T2 and T3 had significant impacts on the soil fungal α-diversity (Table 3). This may be because the structural characteristics of the biochar changed the soil environment, providing a good habitat and nutrient source for the growth and reproduction of soil fungi [30]. Alternatively, increased levels of soil organic carbon can also change the physical structure of the soil, which can create good conditions for the improvement of soil fungal diversity. As the amount of biochar increased, soil fungal α-diversity did not show a continuous increasing trend but reached a maximum within a certain range. This may be because the biochar environment is more suitable for the growth and development of fungi at a certain concentration. The extent to which the physical and chemical properties of soil are affected by biochar is related to the amount of biochar applied. Soil affected by large doses of biochar is not suitable for the reproduction and growth of soil fungi [31], resulting in a decrease in the diversity of soil fungi. A certain amount of biochar can significantly improve the activity of fungi, but a large amount has an inhibitory effect. The reason for this is similar to bacteria. Excess biochar will aggravate change of the soil physical and chemical properties, and excessive change of soil physical and chemical properties is not conducive to the growth of soil fungi. Han's research results are consistent with this view [32].

Fungal Beta Diversity
The species composition of soil fungi is an ordered collection of different microorganisms that are simultaneously related to each other, under specific time and space conditions, in a habitat with obvious characteristics. Its basic characteristics include structure (such as predator relationships), community environment, boundary features, and distribution range [33]. In this study, soil fungal community structure changed significantly after biochar application. Compared with the control, the soil fungal communities of the biochar treatments changed in a specific direction. Soil fungal community structure also showed large differences between the treatment groups, depending on the amount of biochar applied ( Figure 2). This may be because the ecological conditions and soil types of the soils in the different treatments were similar, and these conditions promoted the same dominant fungal species in the different treatments. These changes affect the growth and development of soil fungi, resulting in changes in fungal community structure [34]. There were also significant differences in the soil fungal community structure between the different biochar treatments. This may be because different amounts of biochar affect soil physical and chemical properties differently. The growth and reproduction of fungi are greatly affected by these changes in the soil environment [1], leading to significant differences in the diversity and community structure of the soil fungi.

The Effect of Biochar on Soil Enzyme Activity
Soil enzyme activity is an important reference index for soil microbial activity, biochemical reaction capacity, soil nutrient cycling, and material metabolism. The activities of soil enzymes are extremely sensitive to environmental changes and can reflect changes in soil quality under various conditions [35]. Biochar can cause changes in soil physical and chemical properties through its own physical and chemical properties. It therefore has important impacts on soil enzyme activities [36]. Under the influence of biochar, enzyme activities related to the utilization of phosphorus and nitrogen in the soil can be significantly increased. Jones et al. [34] showed that the continuous application of 50 t ha −1 biochar to the soil significantly increased denitrification enzyme activity. Lehmann et al. believe that biochar can increase the activity of most soil enzymes involved in N mineralization [36]. Urease is a key regulator of the nitrogen cycle in the soil. It is mainly involved in promoting the hydrolysis of urea and is an important reference indicator that reflects soil nitrogen mineralization ability and nitrogen level. Catalase, an oxidoreductase, is an important participant in the decomposition of hydrogen peroxide in the soil. Sucrase is mainly used to hydrolyze sucrose in the soil and convert it into fructose or glucose to facilitate its absorption by plants and soil microorganisms [36]. Polyphenol oxidase is also an oxidoreductase. It is mainly derived from plant residues, root exudates, and soil microorganisms, and is an important participant in soil aromatic compound cycling [1]. Biochar increases soil organic carbon and active organic carbon, both of which are key indicators that affect soil enzyme activity. Therefore, the activities of soil invertase, soil urease, polyphenol oxidase, and catalase vary. These enzyme activities were all improved by biochar application, which is consistent with the results of Huang et al. [37]. This may be because soil enzyme activity is closely related to soil respiration intensity, microbial community structure characteristics, and organic matter content. After biochar application, the above indicators affect changes in soil enzyme activity. The level and coordination of the mineral element content in the soil are key factors that influence plant growth and development, especially plant leaf respiration and photosynthesis. It is the energy carrier and material basis for plant carbon assimilation and metabolism [38].

The Effect of Biochar on Soil Microbial Biomass Carbon and Microbial Biomass Nitrogen
Soil MBC and MBN are the most active components of soil organic carbon and nitrogen. Generally, MBC refers to a carbon component with a volume of 5 × 10 3 µo 3 [39]. MBC and MBN represent the nutrients in the soil that are easily used by plants and microorganisms and are the driving force for carbon and nitrogen mineralization. Moreover, MBC, MBN, and the absorption and circulation of nutrients such as nitrogen, phosphorus, potassium, sulfur, magnesium, and iron in the soil are closely related. Although MBC only accounts for approximately 1% of soil TOC and MBN about 2% of soil TN, they can reflect the activity, abundance, and metabolic intensity of soil microorganisms. In addition, MBC and MBN directly participate in soil metabolism. They represent essential turnover material for effective nutrients in the soil and occupy important positions in the cyclic transformation of soil materials. Therefore, MBC and MBN are of great significance in plant growth and metabolism, and soil fertility maintenance [40]. The results of this study show that biochar effectively increases the content of MBC in the soil, and this effect increases with increasing biochar input. Biochar is loose and porous, providing excellent living conditions for the growth of microorganisms [41]. In addition, the surface structure and carbon content of biochar promotes soil microorganism growth. Biochar can retain water and reduce nutrient leaching by adsorbing cations and anions [42], which indirectly improves microbial nutrient utilization, leading to an increase in MBC and MBN. Previous research [41], has shown that, owing to its chemical stability, biochar application can increase the level of soil MBC, and as the amount of biochar in the soil increases, soil MBC and MBN also increase. However, other studies have shown that biochar has an inhibitory effect on MBC and MBN after being applied to the soil, and it is hypothesized that this phenomenon may be caused by the adsorption of certain carbon and nitrogen compounds by the biochar, as well as its compounding with biochar. This limits the availability of substrates for microbial growth and ultimately leads to a decrease in soil MBC and MBN content [42]. This difference may exist because biochar stimulates the metabolism and reproduction of soil microorganisms [39]. When the microbial community abundance and metabolic intensity of the soil were strong, soil MBC and MBN increased accordingly. In a prosperous period of tobacco transplantation, most of the physiological indicators of tobacco plants reached their maximum values [43], and the temperature and rainfall also increased, promoting the vigorous metabolism of soil microorganisms. The material conversion rate accelerated, the net MBC and MBN in the soil increased, and the influence of biochar on soil MBC and MBN also reached a maximum.

The Relationship between Soil Microbial Diversity and Microbial Biomass Carbon and Nitrogen
"Soil enzymes" is the general term for specialized protein compounds with biocatalytic abilities that exist in soil. They are the products of microbial activity and are extremely small but effective soil components. Various biochemical processes of soil microorganisms are realized by the enzymes they produce. Therefore, the activity of soil enzymes can be used as an indicator of the intensity of soil biochemical processes and to evaluate soil fertility [44]. Soil MBC refers to the sum of the carbon in living and dead microorganisms with a volume of <5000 µm 3 in the soil. Soil microbial biomass refers to the total biomass in the soil with a volume of less than 5 × 10 3 µm 3 and is the most active component of soil organic matter. Between these two, MBC is the more important component. MBC is an easy-to-use nutrient pool in the soil, and the driving force for the decomposition of organic matter and the mineralization of N. It is closely related to the nutrient cycles of C, N, P, and S in the soil. As a component of soil activated carbon, although MBC only accounts for 1-4% of total soil organic carbon, it can not only reflect small changes in total soil carbon, but it also directly participates in the soil. Biochemical transformation is from a reservoir of effective plant nutrients in the soil and can promote the effectiveness of soil nutrients. Therefore, it plays an important role in soil fertility and plant nutrition. Soil enzyme activity and MBC and MBN are important indicators of the metabolic capacity of soil microorganisms. This study shows that soil enzyme activity is positively correlated with soil bacterial and fungal diversity, and that soil microbial carbon and nitrogen are positively correlated with soil microbial diversity. This may be because soil enzymes are mainly produced by soil microorganisms, and they have a strong correlation.; thus, soil enzyme activity and microbial diversity have a mutually promoting relationship [45]. Soil MBC and MBN are present in microorganisms and metabolites while soil MBC and MBN content is an important reference index for soil microbial metabolism. Moreover, when the intensity of soil microbial metabolism increases, the MBC and MBN content of the soil also increases.

Conclusions
This study aimed to study the effects of biochar on the metabolic capacity of tobacco rhizosphere microbial diversity. Biochar increased the diversity of soil bacteria and fungi, and the soil community structure also changed under the influence of biochar. Soil enzyme activity, MBC, and MBN also increased with increasing biochar input. In terms of the effect of biochar on the metabolism of tobacco-growing soil microbial diversity, treatment with 3% biochar was the most effective. Soil enzyme activity is positively related to microbial diversity, and the greatest correlation was found between catalase and invertase and bacterial diversity. There is a strong positive correlation between soil MBC and MBN and soil microbial diversity, and the correlation between soil MBN and bacterial diversity was the strongest; the correlation between soil MBC and fungal diversity is the largest. Biochar plays an important role in improving soil microbial diversity and facilitating microbial metabolism. Thus, our results provide guidance regarding the impressive benefits of biochar application to agricultural soils. Funding: This research was funded by the study on the mechanism and control technology of organic carbon based on carbon and nitrogen balance in tobacco field soil, China (110201902004), and study on improving tobacco leaf quality in Henan Production Area based on soil remediation technology, China (2021410000340260), and study on the improvement of ripening resistance and matching baking technology of upper leaves in Nanyang tobacco-growing area, China (2020411300270065).