ABSTRACT Studying the obstacles associated with continuous cropping is necessary for sustainable agricultural production. Phenolic acids play an important role in continuous cropping

systems, although their mechanism of action in these systems remains unclear. Using High-performance Liquid Chromatography, we characterized the changes in phenolic acid contents in soils

that had been continuously cropped with tobacco for different time periods and evaluated the interactions between soil physicochemical properties, bacterial community structure and

diversity, and phenolic acids. Prolonged continuous cropping was associated with a significant increase in the content of phenolic acids and a significant decrease in soil pH and bacterial

a potential factor leading to changes in bacterial community structure. Continuous cropping lowered the soil pH, which stimulated phenolic acid accumulation and consequently altered the

bacterial community structure and diversity, ultimately impacting tobacco plant growth. SIMILAR CONTENT BEING VIEWED BY OTHERS EFFECT OF STRAW BIOCHAR AMENDMENT ON TOBACCO GROWTH, SOIL

PROPERTIES, AND RHIZOSPHERE BACTERIAL COMMUNITIES Article Open access 20 October 2021 EFFECTS OF SOIL TEXTURE AND NITROGEN FERTILISATION ON SOIL BACTERIAL COMMUNITY STRUCTURE AND NITROGEN

UPTAKE IN FLUE-CURED TOBACCO Article Open access 22 November 2021 IMPACTS OF CONTINUOUS AND ROTATIONAL CROPPING PRACTICES ON SOIL CHEMICAL PROPERTIES AND MICROBIAL COMMUNITIES DURING PEANUT

CULTIVATION Article Open access 17 February 2022 INTRODUCTION The continuous planting of the same crop or related crops in the same plot every year results in poor plant growth and an

increase in pests and diseases1,2,3. Continuous cropping has become widespread in China, with more than 20% of the land displaying the severe negative consequences of continuous cropping,

such as reduced crop yields and poor soil health. This has become more severe as cultivated land resources have become scarcer and agriculture has become more profit-driven4. The problems

associated with continuous cropping have thus attracted increasing attention from scientists. Long-term continuous cropping results in a series of soil problems5,6,7 that lead to increased

crop pests and diseases, and reduced yield and quality8. Numerous studies have shown that the effect of continuous cropping on the soil environment mainly includes the following four

aspects: deterioration of soil physicochemical properties8,9, reduction in enzyme activity10,11,12, accumulation of toxic substances13,14,15, and changes in microbial diversity16,17.

Although previous studies have thoroughly assessed the above four aspects, research regarding their interactions, especially the relationship between toxic substances and other environmental

factors, is still unclear and requires further investigation. Numerous studies have shown that the deterioration of the soil micro-ecological environment is a primary obstacle in continuous

cropping18, and the role of toxic substance accumulation has garnered increasing research interest13,14,15. There is substantial evidence that the accumulation of phenolic acids is an

important contributor to the issues associated with continuous cropping19,20. Previous studies have shown that phenolic acids, including p-hydroxybenzoic acid, vanillic acid, p-coumaric

acid, and ferulic acid, are widely present in the soil and inhibit the growth of crops21,22. Phenolic acids are released into the soil by plant evaporation, leaching, root secretion, litter,

and the decomposition of their residues23,24. Their contents are affected by soil physicochemical properties, especially pH25,26, while microbes degrade phenolic acids and affect their

allelopathy27. Phenolic acids, as organic substances, will inevitably affect the community structure and diversity of microorganisms while simultaneously providing a C and N source for the

microorganisms. Although many studies have shown that the accumulation of phenolic acids alters the microbial community structure16, its specific mechanism of influence is unclear. Soil

microorganisms are important for maintaining soil quality. They regulate the decomposition of plant residues and soil organic matter, the biochemical cycle, and the formation of soil

structure28,29. Many studies have shown that long-term continuous cropping leads to the transformation of the soil microbial community structure from “bacterial type” to “fungal type”16,

which destroys the soil microbial community structure, leading to an increase in harmful pathogens and a decrease in beneficial flora, which ultimately leads to poor soil health17,30. It is

thus necessary that the relationships between soil environmental factors, particularly phenolic acids and soil microorganisms, are studied in continuous cropping systems. Previous studies

have shown that when phenolic acids enter the soil, they are not only able to inhibit seed germination and seedling growth31, but can also alter the microbial community structure and

diversity as a result of their particular properties or by acting as signaling molecules19. Phenolic acids can significantly affect the biomass, diversity, and community structure of soil

microbes, and selectively increase the specific microbial species in the soil32,33. For example, vanillic acid in soybean root exudates had a considerable impact on rhizosphere microbial

communities33. In addition, phenolic acids can also indirectly stimulate soil-borne pathogenic microorganisms and increase plan morbidity, which is a major factor contributing to their role

in continuous cropping34. For instance, coumarin in cucumber root exudates can significantly affect rhizosphere soil microbial communities and promote the growth of soil-borne pathogens35,

while vanillin can significantly affect the community structure and abundance of _Fusarium_ and _Trichoderma_ in cucumber rhizomes36. Based on the above studies, it is evident that phenolic

acids are a major factor contributing to the obstacles associated with continuous cropping and should be evaluated further. In the present study, the changes in physicochemical properties,

soil enzyme activities, bacterial diversity, and phenolic acids in soil that had been continuously cropped with tobacco for different numbers of years were assessed. The objectives of this

study were to: 1) explore the factors contributing to the variation in environmental characteristics, particularly the phenolic acid contents, in soil that had been continuously cropped for

different durations, 2) evaluate whether any correlation between soil physiochemical properties, bacterial abundance, bacterial community structure, and phenolic acids exists, 3) improve our

understanding of the contribution of phenolic acids to the soil micro-ecological environment under continuous cropping systems, and (4) assess the impacts of continuously cropped soil on

tobacco plant growth. RESULTS PHYSIOLOGICAL INDEXES OF THE TOBACCO PLANTS The difference in the maximum net photosynthetic rate (Pnmax) of the tobacco plants under continuous cropping for 4,

6, and 8 years was not significant; however, it was significantly reduced after continuous cropping for 14 years (Fig. 1a). The Pnmax of “T16” was 39.85%, 38.04%, 33.47%, and 16.87% lower

than that of the first four treatments, respectively, and the root/shoot ratio also showed the same decreasing trend (Fig. 1b). Although the difference in root activity [measured using

triphenyl tetrazolium chloride (TTC)] between “T6” and “T8” and “T14” and “T16” was not significant, the activity of “T4” was significantly higher than the other treatments (Fig. 1c).

Conversely, the malondialdehyde (MDA) content of “T16” increased significantly by 69.04%, 48.80%, 44.25%, and 13.42%, respectively, compared with the first four treatments (Fig. 1d). SOIL

PHYSICOCHEMICAL CHARACTERISTICS AND ENZYME ACTIVITIES The soil pH, available P, and organic matter (OM) content decreased with the long-term continuous cropping of tobacco (Fig. 2a,b,d),

while the available K content increased over time (Fig. 2e). Soil catalase (CAT), urease, acid phosphatase, and invertase activities were assessed to measure the potential turnover rates of

N or C in the soils from the five sites. As shown in Fig. 3, the activity of these four enzymes showed a trend of increasing first and then decreasing with the increase in continuous

cropping duration. In addition, urease and sucrase activity peaked in the sixth year, while CAT and acid phosphatase activities peaked in the eighth year. PHENOLIC ACID CONTENTS The contents

of the phenolic acids in the five soils were measured by High-performance Liquid Chromatography (HPLC). As shown in Fig. 4, the contents of phloroglucinol, coumalic acid, _P_-hydroxybenzoic

acid, vanillic acid, ferulic acid, and the total amount increased with the long-term continuous cropping of tobacco. The content of the above phenolic acids in “T16” increased by 61.41%,

168.00%, 267.72%, 280.82%, 223.33%, 893.75%, and 109.92%, respectively, compared with “T4”. In addition, the contents of benzoic acid and cinnamic acid appeared to increase first and then

decrease with increasing cropping duration, peaking in the eighth year. BACTERIAL COMMUNITY COMPOSITION All the sample sequences were classified into 47 phyla using Mothur. The top 10 phyla

are shown in Fig. 5. The overall bacterial composition of the different treatments was similar, but the distribution of each phylum varied. _Proteobacteria_ (50.27% of all sequence reads),

_Actinobacteria_ (13.16%), _Bacteroidetes_ (14.47%), _Acidobacteria_ (7.23%), and _Gemmatimonadetes_ (3.64%) were the five most dominant phyla, accounting for 88.76% of the reads.

Additionally, the relative abundance of _Proteobacteria_ decreased with the long-term continuous cropping of tobacco, while that of _Acidobacteria_ and _Actinobacteria_ appeared to increase

initially and then decrease as the cropping duration increased. The abundance of _Bacteroidetes_ first decreased but then increased substantially by T16. Furthermore, comparison of the

relative abundances of the top 35 classified bacterial genera showed significant variations among the five sites (Fig. 6), particularly with regards to T16 and the other treatments. The

relative abundances of _Chlorobium_, _Deltaproteobacteria_, and _Betaproteobacteria_ significantly decreased with increased cropping duration. BACTERIAL ALPHA DIVERSITY The Good’s coverage

values of the five samples were 97.0%, 96.8%, 97.0%, 98.4%, and 97.5%, indicating that sequencing reads were sufficient for this analysis (Fig. 7a). The Shannon, Simpson, ACE, and Chao1

indexes decreased with the duration of cropping, although not all interactions were significant. However, all of the above indexes decreased significantly after 16 years of monoculture (Fig.

 7b–e). Furthermore, the number of observed species in the bacterial communities significantly decreased after 14 and 16 years of continuous cropping (Fig. 7f). BACTERIAL BETA DIVERSITY

Bray-Curtis dissimilarity hierarchical clustering analysis showed that the samples from the same sites grouped together and were more similar than the samples between different sites (Fig. 

8). The “T4”, “T6”, “T8”, and “T14” sites grouped together and were separated from “T16”; the “T4”, “T6”, and “T8” sites grouped together and were separated from “T14”; and “T4” and “T6”

grouped together and were separated from “T8”, which suggested that the soil bacterial community structure might be affected by monoculture duration. UniFrac-weighted PCoA based on the

operational taxonomic unit (OTU) composition also clearly demonstrated variations among these different soil samples, with the first two axes explaining 52.26% and 24.38% of the total

variation for the bacterial data. In addition, the bacterial communities in the four-year soil sample (T4), six-year soil sample (T6), and eight-year soil sample (T8) were obviously

separated from the 14-year soil sample (T14) and the 16-year soil sample (T16). Furthermore, the four-year and six-year soil samples had the most similar bacterial community memberships

(Fig. 9). The unweighted UniFrac algorithm displayed similar results; for clarity, only the weighted UniFrac PCoA plot is shown here. CORRELATIONS BETWEEN SOIL PHYSICOCHEMICAL PROPERTIES AND

PHENOLIC ACIDS The Mantel test revealed a significant (r = 0.700, _P_ < 0.01) correlation between the soil properties and phenolic acid contents. Spearman’s correlation coefficient was

used to assess the relationships between the soil properties and phenolic acids (Table 1). Phloroglucinol, _P_-hydroxybenzoic acid, vanillin, ferulic acid, and total amount were negatively

correlated with soil pH. No other significant correlations were observed, with the exception of P-hydroxybenzoic acid and available K. The above results indicate that soil pH and phenolic

acid content are closely associated. Furthermore, stepwise regression analysis also demonstrated that pH was the only physicochemical property that was significantly associated with phenolic

acid content (Table 2). CORRELATIONS BETWEEN SOIL PHYSICOCHEMICAL PROPERTIES AND ABUNDANT PHYLA The Mantel test analysis revealed a significant (r = 0.440, _P_ < 0.01) correlation

between the soil properties and the abundance of the bacterial phyla. The first two components of the RDA explained 56.2% and 18.7% of the total variation in bacterial phyla data (Fig. 10).

The first component (RDA1) separated “T4” and “T6” from the “T8”, “T14”, and “T16” sites, and the second component also exhibited the same trend (RDA 2). Spearman’s correlation was used to

assess the relationships between the soil properties and abundance of the bacterial phyla (Table 3). The abundance of _Proteobacteria_ was positively associated with soil OM. The relative

abundances of _Bacteroidetes_ and _Firmicutes_ were positively correlated with soil AN, while the abundances of _Verrucomicrobia_, _Acidobacteria_, and _Planctomycetes_ were negatively

correlated with soil AN. CORRELATIONS BETWEEN PHENOLIC ACIDS AND ABUNDANT PHYLA The Mantel test analysis revealed a significant (r = 0.480, _P_ < 0.01) correlation between the soil

phenolic acids and the abundance of the bacterial phyla. The first two components of the RDA explained 54.8% and 26.2% of the total variance in the bacterial phyla data (Fig. 11). The first

component (RDA1) separated “T4”, “T6”, “T8”, and “T14” from “T16”, and the second component (RDA 2) separated “T4”, “T6”, and “T8” from “T14” and “T16”. Spearman’s correlation was used to

assess the relationships between the soil phenolic acids and bacterial phylum abundance (Table 4). The abundances of _Bacteroidetes_ and _Firmicutes_ were positively correlated with coumalic

acid, vanillic acid, vanillin, and total amount. The abundance of _Actinobacteria_ was positively associated with cinnamic acid. DISCUSSION Previous field experiment studies have shown that

continuous cropping leads to a decrease in tobacco plant height, biomass, and root activity37. In this study, in addition to the decrease in growth indicators (Fig. 1a–c), we also found

that continuous cropping caused a significant increase in MDA content in the tobacco plants (Fig. 1d). MDA content is representative of the degree of membrane peroxidation in plant cells38,

and the significant increase in MDA content in this study indicates that cells experience oxidative damage under monoculture stress, which may result in an inhibition of tobacco plant

growth. Based on the above analysis, we speculate that continuous cropping leads to the accumulation of toxic substances in the soil (such as phenolic acids), which may cause cell membrane

peroxidation via allelopathy39,40,41. Additionally, toxic phenolic acids can also indirectly suppress tobacco plant growth by altering the soil micro-ecological environment42,43. Continuous

cropping was thus associated with a stressful growth environment for the tobacco plants. Our results showed that continuous cropping leads to soil acidification and nutrient imbalances as

well as changes in enzyme activity, characterized by an initial increase followed by a decrease. Specifically, our findings indicate that continuous cropping results in the acidification of

the soil, and the acidified soil provides suitable breeding conditions for pathogens44, which is one of the major obstacles associated with continuous cropping. The increase in soil

available K content may be due to the excessive quantities of K2SO4 that are added during cultivation, as tobacco is a K-favoring crop. This excessive K+ is exchanged with cations, such as

Ca2+ and Mg2+, which further leads to soil compaction and lower OM contents, which in turn aggravates the issues associated with continuous cropping. Previous studies have reported that

continuous cropping may result in decreased enzyme activity10, but our study found that as the duration of continuous cropping increased, the activity of the four enzymes (CAT, urease,

invertase, acid phosphatase) increased first and then decreased (Fig. 3). One explanation for this trend is that prior to the peak period of enzyme activity, the tobacco plants and root

residues in the surface soil accumulate yearly, providing sufficient substrates for the enzymes and thus resulting in increased enzyme activity. Once the continuous cropping exceeds the peak

period of enzyme activity, the accumulation of pathogens and toxic substances in the soil eventually leads to a decrease in enzyme activity. In recent years, increasing numbers of studies

have demonstrated that phenolic acids have strong allelopathic activity45,46. Although the accumulation of phenolic acids in continuously cropped soil has also been widely confirmed, only

3–4 species have been studied in many reports20,47,48. In this study, the contents of phloroglucinol, coumalic acid, p-hydroxybenzoic acid, vanillic acid, vanillin, and ferulic acid

accumulated significantly with increased continuous cropping duration. A sharp increase in the contents of these phenolic acids was observed after 14 years of continuous cropping, which may

be due to the deterioration of the microbial community structure, which resulted in the slow degradation of phenolic acids, eventually leading to their accumulation. Interestingly, we found

that benzoic acid and cinnamic acid first increased and then decreased with the increase in continuous cropping duration, peaking after eight years. It may be that benzoic acid and cinnamic

acid accumulate at the initial stage of continuous cropping, and as the duration of the continuous cropping increases, the microorganisms that degrade these two phenolic acids become

dominant bacteria. This finding has theoretical implications for the screening of novel strains that degrade phenolic acids. The observed changes in the total phenolic acid content in the

tobacco-planted soils confirm that long-term continuous cropping leads to the accumulation of phenolic acids that have negative consequences for soil health, further verifying that phenolic

acids constitute an important obstacle in continuous cropping. Our results indicate that continuous cropping leads to changes in soil bacterial community structure and reduced bacterial

diversity. Sequence analysis revealed that _Proteobacteria_, _Actinobacteria_, _Bacteroidetes_, _Acidobacteria_, and _Gemmatimonadetes_ were the most abundant phyla in the five soil samples

(Fig. 5), and this result is generally consistent with the results of previous studies49,50,51,52. However, some differences in the bacterial species and abundance exist between this study

and previous studies, which may be due to differences in climate, vegetation, and soil type. The overall trend of our data indicated that the relative abundances of _Proteobacteria_,

_Chlorobia_, _Deltaproteobacteria_, and _Betaproteobacteria_ decreased significantly with increasing cropping duration. Furthermore, there was a great difference between T16 and the other

treatments (Fig. 6), indicating that long-term continuous cropping most significantly contributes to changes in bacterial community structure. As the period of continuous cropping increased,

bacterial diversity showed a significant downward trend (Fig. 7), indicating that reduced bacterial diversity is another major contributor to the issues associated with continuous

cropping53,54. In addition, both hierarchical clustering and PCoA analyses indicated that continuous cropping of tobacco strongly influenced the variation in bacterial community structure.

Particularly, the 14-year and 16-year soils differed from the four-year, six-year, and eight-year soils (Figs 8 and 9), indicating that the long-term continuous planting of tobacco appears

to result in a shift in the original state of the soil microbial communities. The Mantel test indicated a certain correlation between soil physicochemical properties and phenolic acids.

Specifically, the correlation analysis demonstrated that the contents of phloroglucinol, _P_-hydroxybenzoic acid, vanillin, ferulic acid, and total amount were negatively associated with

soil pH (Table 1). Although the study by Whitehead _et al_.25 also showed a negative correlation between pH and phenolic acid content, no significant correlation was observed between the pH

values of 6–8. In order to further clarify the relationship between the indicators of the physicochemical properties and phenolic acids, we used stepwise regression analysis to eliminate the

variables that were unrelated or less influential. The results showed that pH has a certain inhibitory effect on the content of phloroglucinol, _P_-hydroxybenzoic acid, vanillin, ferulic

acid, and total amount (Table 2); that is, the lower the pH value, the more obvious the accumulation of phenolic acids. Phenolic acids are introduced into the soil through several pathways,

such as root exudation, pollen transmission, and spoilage degradation of the residue55. After entering the soil, they are consumed and decomposed as substrates56. However, under continuous

cropping conditions, soil acidification is caused by the excessive application of acidic chemical fertilizers57. The acidic environment provides suitable conditions for the accumulation of

phenolic acids25, and the accumulation of phenolic acids exacerbates the acidification of the soil, eventually resulting in an undesirable feedback loop. Based on the above analysis, we

speculate that continuous cropping leads to soil acidification, which is turn leads to the accumulation of phenolic acids. We thus recommend that the soil pH is adjusted in agricultural

production to reduce the accumulation of toxic phenolic acids. Previous studies showed that the soil microbial community distribution is associated with these types and concentrations of

soil phenolic acids and inorganic nutrients33,58,59. In this study, the Mantel test showed that phenolic acids and soil physicochemical properties were significantly correlated with

bacterial abundance (Tables 3 and 4). Notably, the correlation between phenolic acids (r = 0.480**) and bacterial abundance was higher than the soil physicochemical properties (r = 0.440**).

In terms of the correlation between physicochemical properties and bacterial abundance, the abundance of _Proteobacteria_ was positively associated with soil OM, while that of

_Bacteroidetes_ and _Firmicutes_ was positively correlated with soil AN. In contrast, the abundances of _Verrucomicrobia_, _Acidobacteria_, and _Planctomycetes_ were negatively correlated

with soil AN. However, pH was not significantly correlated with bacterial abundance (Fig. 10 and Table 3), suggesting that pH may indirectly affect bacterial community structure by affecting

other environmental factors (phenolic acids, for example) in the soil. Although the total phenolic acid content of the soil was found to increase significantly following two years of

continuous _Rehmannia glutinosa_ cropping, and expression analysis of different proteins in the rhizosphere soil of a continuous cropping system indicated that the accumulation of phenolic

acids caused changes in microbial community structure47,48, the above studies did not report on the types of phenolic acids detected. In this study, we further clarified the relationship

between different phenolic acids and bacterial abundance. Specifically, coumalic acid, vanillic acid, vanillin, and total amount were significantly associated with _Bacteroidetes_ and

_Firmicutes_ (Fig. 11 and Table 4). Phenolic acids can provide microorganisms with a carbon source and energy for growth and reproduction56, which leads to changes in the microbial community

structure60,61. Previous studies have shown that phenol 2,4-di-tert-butylphenol and vanillic acid can selectively enhance the types of certain microorganisms in the soil33. In this study,

in addition to vanillic acid, we also found that an increase in coumalic acid and vanillin levels could increase the abundance of _Bacteroidetes_ and _Firmicutes_. This might be because

these two bacteria can selectively use coumalic acid, vanillic acid, and vanillin as substrates for growth and reproduction. Based on the above results, we speculate that the accumulation of

phenolic acids alters the community structure of the soil bacteria by altering the relative abundance of _Bacteroidetes_ and _Firmicutes_. This leads to a deterioration in the soil

micro-ecological environment may due to the proliferation of pathogenic bacteria and the reduction of beneficial bacteria, etc.33,60,61, which ultimately indirectly inhibits tobacco plant

growth. In conclusion, continuous cropping resulted in a lowered soil pH, which stimulated the accumulation of phenolic acids and consequently altered the bacterial community structure and

diversity, ultimately leading to poor tobacco plant growth. The above results provide a theoretical basis for further understanding the mechanisms associated with the obstacles of continuous

cropping. On this basis, in agricultural production, soil pH, OM, and soil microbiota can be improved through various methods, such as the use of quicklime, green manure, and organic

fertilizers, which would provide a better micro-ecological environment for crop growth. METHODS SITE DESCRIPTION AND SAMPLING The experimental site was located at a contiguous flue-cured

tobacco cultivation demonstration area in Shidian County, Yunnan Province, China (99°14′E, 24°36′N), which has a tropical subtropical monsoon climate. The mean annual temperature and

precipitation in the area are 17 °C and 1120 mm, respectively. Before the tobacco was planted, the area was a paddy field with the same physical and chemical properties, geographical

factors, and ecological environment. Once the tobacco was planted, the annual planting varieties, fertilization regime, and management methods remained the same. No major environmental

events occurred during this period. The land in this area has been cultivated by the same farmer, and thus the history of the tobacco planting, including the continuous cropping period,

could be determined. At all sites, the soil texture is sandy loam and the flue-cured tobacco variety that had been planted was ‘K326’ (_Nicotiana tabacum_). The annual nitrogen (N)

application rate (pure N) was 120 kg∙hm−2, the phosphorus (P) application rate (P2O5) was 90 kg∙hm−2, and the potassium (K) application rate (K2O) was 321 kg∙hm−2. The fertilizer was

provided by a local tobacco company, and the cultivation management practices were carried out in accordance with local high-quality tobacco production protocols. The soil samples were

collected on September 20, 2017, from five tobacco fields with continuous cropping histories of 4, 6, 8, 14, and 16 years, and were marked as “T4”, “T6”, “T8”, “T14”, and “T16,”

respectively. The five tobacco fields were each approximately 2000 m2 in size. For each field, three biological samples were randomly collected, and each biological sample consisted of a

composite of five soil cores randomly collected from the tobacco-planted soil (3.0 cm in diameter, 0–20 cm in depth). After passing through a 2-mm sieve, each soil sample was divided into

three subsamples: one portion for the pot experiment, one portion was air-dried for soil characteristics and soil enzyme activity analysis, while the remainder was snap-frozen in liquid

nitrogen and stored at −80 °C for the determination of microbes and phenolic acids. POT TRIALS A pot experiment was performed to evaluate how the soils from the five different cropping

treatments impacted the growth and development of tobacco at the seedling stage from 2017 to 2018. The flue-cured tobacco seedling variety used was K326 and had been cultivated by floating

seedlings. Each pot (inner diameter 15 cm, height 15 cm) contained 1.0 kg of soil and one flue-cured tobacco seedling. Each treatment was in triplicate (three blocks), and each block had 10

pots. All the pots were randomly placed in a light incubator with an average temperature of 23 °C and an average humidity of 55%. After 30 d, root/shoot ratio, malondialdehyde content (MDA),

root activity, and the maximum net photosynthetic rate (Pnmax) of the flue-cured tobacco plants were measured. Determination of MDA content was by the thiobarbituric acid method. The root

activity was measured according to the triphenyl tetrazolium chloride (TTC) method. The net photosynthetic rate was measured using a LI-6400 portable photosynthetic apparatus (Li-COR Inc.,

USA) with a blue-red light source probe (Li-6400-02B). The photosynthetically active radiant (PAR) was set to 1800 μmol/(m2•s), the CO2 injection system was set to 400 μmol/mol, the gas flow

rate was 500 mmol/s, and the leaf temperature was set to 25 °C. For the above indicators, five seedlings exhibiting uniform growth in each block were selected, and the average value was

calculated following the measurements. DETERMINATION OF PHENOLIC ACID CONTENTS IN THE FIELD-COLLECTED SOIL The soil phenolic compounds were identified by High-performance Liquid

Chromatography (HPLC). One hundred grams of fresh sample was placed into a centrifuge tube (three biological replicates), to which 100 mL of 1 mol∙L−1 NaOH was added, allowing for extraction

to occur over 24 h. The samples were then shaken for 60 min. Following centrifugation, the solution was acidified to pH 2.5 using 12 mol∙L−1 HCl, and after 120 min, the humic acid was

removed by centrifugation, following which the supernatant was passed through a fiber membrane of 0.22 μm. The filtrate was then measured by HPLC, and the results were converted according to

the weight of the dried soil. The instrument used was an Agilent 1200 (Agilent Technologies, Santa Clara, CA, USA) with a SunFireTM C18 column (4.6 mm × 250 mm, 5 μm) with a flow rate of 1 

mL·min−1, a column temperature of 25 °C, and a UV detection wavelength of 280 nm. The specific gradients for mobile phase A (methanol) and mobile phase B (aqueous acetic acid solution with

pH = 2.5) were: 0 min, mobile phase A 30%, B 70%; 15 min, mobile phase A 50%, B 50%; 16 min, mobile phase A 70%, B is 30%; 30 min, mobile phase A 0%, B is 100%. According to the results of

previous studies, 10 phenolic acids were selected as standard samples for determination21,22. The standard samples (Sigma) included valproic acid, p-hydroxybenzoic acid, phloroglucinol,

vanillic acid, vanillin, ferulic acid, phthalic acid, benzoic acid, cinnamic acid, and salicylic acid, and the injection volume was 10 μL. There was a 10 min delay at the end of each cycle

to remove the effects of interfering components and ensure the stability and repeatability of the results. The samples were identified based on the retention time and peak area of the

analytical standards. The chromatogram of each standard sample is shown in Figs 12–17. DETERMINATION OF SOIL PHYSICOCHEMICAL PROPERTIES AND SOIL ENZYME ACTIVITIES The soil pH was measured by

a potentiometer using a water to soil ratio of 2.5:1 (NY/T 1377-2007). The organic matter (OM) was determined using the potassium dichromate-sulfuric heating method (NY/T 87-1988);

available N (AN) was determined by the alkali diffusion method (LY/T 1229-1999); available P (AP) was determined by NaHCO3 leaching and the molybdenum antimony colorimetric method (NY/T

1121.7-2006); available K (AK) was determined by ammonium acetate extraction and flame photometer (NY/T 889–2004). The main instruments used for the measurements included a flame photometer,

analytical balance, Kay-type nitrogen analyzer, electrothermal constant temperature blast drying oven, and an ultraviolet spectrophotometer. Catalase (CAT) activity was determined by KMnO4

titration; urease activity was determined using the indophenol blue colorimetric method; invertase activity was assessed using the 3,5-dinitrosalicylic acid colorimetric method; and acid

phosphatase activity was evaluated using the phenylphosphonium phosphate colorimetric method. All the above soil samples were measured in three replicates. DNA EXTRACTION, 16S RRNA GENE

AMPLIFICATION, AND SEQUENCING Total genomic DNA from 15 samples (three replicates from each group) was extracted using the cetyltrimethylammonium bromide/sodium dodecyl sulfate (CTAB/SDS)

method62. DNA concentration and purity were assessed on 1% agarose gels. Based on the concentration, the DNA was diluted to 1 ng/μL using sterile water. The 16S rDNA V4 region was amplified

using primers 515 F (5′-GTGCCAGCMGCCGCGG-3′) and 806 R (5′-GGACTACHVGGGTWTCTAAT-3′). All PCR reactions were carried out using Phusion® High-Fidelity PCR Master Mix (New England Biolabs,

Ipswich, MA, USA). The same volume of 1 × loading buffer (containing SYB green) was combined with the PCR products and subjected to 2% agarose gel electrophoresis for detection. PCR products

was mixed in equimolar ratios. The PCR products were then purified using a GeneJETTM Gel Extraction Kit (Thermo Scientific, Waltham, MA, USA). Sequencing libraries were generated using the

Ion Plus Fragment Library Kit 48 rxns (Thermo Scientific) following the manufacturer’s recommendations. The library quality was assessed using a Qubit@ 2.0 fluorometer (Thermo Scientific).

Ultimately, the library was sequenced on an Illumina HiSeq2500 platform and 250-bp paired-end reads were generated. SEQUENCING DATA ANALYSIS PAIRED-END READS QUALITY CONTROL Paired-end reads

were assigned to samples using Cutadapt63 based on their unique barcode and were truncated by cutting off the barcode and primer sequence. Quality filtering on the raw reads was performed

under specific filtering conditions to obtain high-quality clean reads according to the Cutadapt quality-controlled process. The reads were compared with the reference database UCHIME

algorithm64 to detect chimeric sequences, following which the chimeric sequences were removed65, resulting in clean reads. OTU CLUSTER AND SPECIES ANNOTATION Sequence analysis was performed

using Uparse software (Uparse version 7.0.1001)66. Sequences with ≥97% similarity were assigned to the same OTUs. Representative sequences for each OTU were then screened for further

annotation. For each representative sequence, the Silva Database67 was used for taxonomic information based on the Mothur algorithm. To evaluate the phylogenetic relationships of the

different OTUs and the differences in the dominant species in the different groups, multiple sequence alignment was conducted using MUSCLE (version 3.8.31)68. OTU abundance information was

normalized using the standard of the sequence number corresponding to the per sample with the least sequences. Subsequent analysis of alpha diversity and beta diversity were all performed

basing on this normalized data. ALPHA DIVERSITY Alpha diversity is used for analyzing the complexity of species diversity for a sample via six indices, including Observed-species, Chao1,

Shannon, Simpson, ACE, and Good’s coverage. All these indices were calculated for our samples with QIIME (version 1.7.0) and visualized with R software (version 2.15.3). BETA DIVERSITY Beta

diversity analysis was used to evaluate the differences in species complexity between the samples. Beta diversity on both weighted and unweighted UniFrac values were calculated using QIIME

software (version 1.7.0). Principal component analysis (PCA) was used to reduce the dimensionality of the original variables using the FactoMineR package and ggplot2 package in R software.

Principal Coordinate Analysis (PCoA) was performed to obtain the principal coordinates and visualize from complex multidimensional data. A distance matrix of the weighted or unweighted

UniFrac values among samples was transformed into a new set of orthogonal axes, by which the maximum variation factor is demonstrated by the first principal coordinate, and the second

maximum by the second principal coordinate, and so on. PCoA analysis was displayed using the Weighted Correlation Network Analysis (WGCNA) package, stat packages, and ggplot2 package in R.

Unweighted Pair-group Method with Arithmetic Means (UPGMA) hierarchical clustering was used to interpret the distance matrix using average linkage and was conducted in QIIME software.

STATISTICAL ANALYSES Using SPSS version 22.0 (SPSS Inc., Armonk, NY, USA) and one-way analysis of variance (ANOVA) with Tukey’s Honest Significant Difference (HSD) tests were conducted for

multiple comparisons for all parameters, as were Spearman’s rank correlations and stepwise regression analysis. The correlations between the abundant bacterial phyla and soil characteristics

were determined by the Mantel test, and redundancy analysis (RDA) was carried out using the vegan package in R. Heatmap figures were generated using custom R scripts. DATA AVAILABILITY

Sequence data were deposited in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under the accession number SRP186892. REFERENCES * Wei, Z. &

Yu, D. Analysis of the succession of structure of the bacteria community in soil from long-term continuous cotton cropping in Xinjiang using high-throughput sequencing. _Archives of

Microbiology_ 200, 653–662 (2018). Article  CAS  Google Scholar  * Machado, S., Petrie, S., Rhinhart, K. & Qu, A. Long-term continuous cropping in the pacific northwest: tillage and

fertilizer effects on winter wheat, spring wheat, and spring barley production. _Soil & Tillage Research_ 94, 473–481 (2007). Article  Google Scholar  * Zhang, X. H., Lang, D. Y., Zhang,

E. H., Bai, C. C. & Wang, H. Z. Diurnal changes in photosynthesis and antioxidants ofangelica sinensisas influenced by cropping systems. _Photosynthetica (Prague)_ 51, 252–258 (2013).

Article  CAS  Google Scholar  * Yuan, L. _et al_. Responses of rice production, milled rice quality and soil properties to various nitrogen inputs and rice straw incorporation under

continuous plastic film mulching cultivation. _Field Crops Research_ 155, 164–171 (2014). Article  Google Scholar  * Wang, L. _et al_. Continuous plastic-film mulching increases soil

aggregation but decreases soil pH in semiarid areas of China. _Soil and Tillage Research_ 167, 46–53 (2017). Article  Google Scholar  * Njoroge, R., Otinga, A. N., Okalebo, J. R., Pepela, M.

& Merckx, R. Occurrence of poorly responsive soils in western Kenya and associated nutrient imbalances in maize (Zea mays L.). _Field crops research_ 210, 162–174 (2017). Article 

Google Scholar  * Xia, Z. C., Kong, C. H., Chen, L. C. & Wang, S. L. Allelochemical-mediated soil microbial community in long-term monospecific Chinese fir forest plantations. _Applied

Soil Ecology_ 96, 52–59 (2015). Article  Google Scholar  * Jones, A. R., Orton, T. G. & Dalal, R. C. The legacy of cropping history reduces the recovery of soil carbon and nitrogen after

conversion from continuous cropping to permanent pasture. _Agriculture, Ecosystems & Environment_ 216, 166–176 (2016). Article  CAS  Google Scholar  * Zhang, Y. _et al_. No-tillage with

continuous maize cropping enhances soil aggregation and organic carbon storage in Northeast China. _Geoderma_ 330, 204–211 (2018). Article  ADS  CAS  Google Scholar  * Xiong, W. _et al_.

The effect of long-term continuous cropping of black pepper on soil bacterial communities as determined by 454 pyrosequencing. _PloS one_ 10, e0136946 (2015). Article  Google Scholar  * Ma,

Y. H., Fu, S. L., Zhang, X. P., Zhao, K. & Chen, H. Y. Intercropping improves soil nutrient availability, soil enzyme activity and tea quantity and quality. _Applied soil ecology_ 119,

171–178 (2017). Article  Google Scholar  * Kotroczó, Z. _et al_. Soil enzyme activity in response to long-term organic matter manipulation. _Soil Biology and Biochemistry_ 70, 237–243

(2014). Article  Google Scholar  * Huang, L. F. _et al_. Plant-soil feedbacks and soil sickness: from mechanisms to application in agriculture. _Journal of chemical ecology_ 39, 232–242

(2013). Article  ADS  CAS  Google Scholar  * Nie, L. _et al_. Alleviating soil sickness caused by aerobic monocropping: responses of aerobic rice to soil oven-heating. _Plant and soil_ 300,

185–195 (2007). Article  CAS  Google Scholar  * Xiong, W. _et al_. Different continuous cropping spans significantly affect microbial community membership and structure in a vanilla-grown

soil as revealed by deep pyrosequencing. _Microbial ecology_ 70, 209–218 (2015). Article  Google Scholar  * She, S. _et al_. Significant relationship between soil bacterial community

structure and incidence of bacterial wilt disease under continuous cropping system. _Archives of microbiology_ 199, 267–275 (2017). Article  CAS  Google Scholar  * Caputo, F., Nicoletti, F.,

Picione, F. D. L. & Manici, L. M. Rhizospheric changes of fungal and bacterial communities in relation to soil health of multi-generation apple orchards. _Biological Control_ 88, 8–17

(2015). Article  Google Scholar  * Yang, R. _et al_. The effects of cattle manure and garlic rotation on soil under continuous cropping of watermelon (Citrullus lanatus L.). _PloS one_ 11,

e0156515 (2016). Article  Google Scholar  * Yin, C. _et al_. How to plant apple trees to reduce replant disease in apple orchard: a study on the phenolic acid of the replanted apple orchard.

_PloS one_ 11, e0167347 (2016). Article  Google Scholar  * Xie, X. G. _et al_. Reduction of soil-borne pathogen Fusarium solani reproduction in soil enriched with phenolic acids by

inoculation of endophytic fungus Phomopsis liquidambari. _BioControl_ 62, 111–123 (2017). Article  CAS  Google Scholar  * Whitehead, D. C. Identification of p-hydroxybenzoic, vanillic,

p-coumaric and ferulic acids in soils. _Nature_ 202, 417 (1964). Article  ADS  CAS  Google Scholar  * Flaig, W. Organic compounds in soil. _Soil Science_ 111, 19–33 (1971). Article  ADS  CAS

  Google Scholar  * Xiao, C. L. _et al_. Autotoxic effects of root exudates of soybean. _Allelopathy Journal_ 18, 121–127 (2006). Google Scholar  * Chen, L. _et al_. Trichoderma harzianum

SQR-T037 rapidly degrades allelochemicals in rhizospheres of continuously cropped cucumbers. _Applied microbiology and biotechnology_ 89, 1653–1663 (2011). Article  CAS  Google Scholar  *

Whitehead, D. C., Dibb, H. & Hartley, R. D. Extractant pH and the release of phenolic compounds from soils, plant roots and leaf litter. _Soil Biology and Biochemistry_ 13, 343–348

(1981). Article  CAS  Google Scholar  * Johnson, R. M. & Pregitzer, K. S. Concentration of sugars, phenolic acids, and amino acids in forest soils exposed to elevated atmospheric CO2 and

O3. _Soil Biology and Biochemistry_ 39, 3159–316 (2007). Article  CAS  Google Scholar  * Inderjit. Soil microorganisms: an important determinant of allelopathic activity. _Plant and Soil_

274, 227–236 (2005). Article  CAS  Google Scholar  * Miki, T. & Doi, H. Leaf phenological shifts and plant–microbe–soil interactions can determine forest productivity and nutrient

cycling under climate change in an ecosystem model. _Ecological research_ 31, 263–274 (2016). Article  Google Scholar  * Macdonald, C. A. _et al_. Physiological, biochemical and molecular

responses of the soil microbial community after afforestation of pastures with Pinus radiata. _Soil Biology and Biochemistry_ 41, 1642–1651 (2009). Article  CAS  Google Scholar  *

Tewoldemedhin, Y. T., Mazzola, M., Labuschagne, I. & McLeod, A. A multi-phasic approach reveals that apple replant disease is caused by multiple biological agents, with some agents

acting synergistically. _Soil Biology and Biochemistry_ 43, 1917–1927 (2011). Article  CAS  Google Scholar  * Salvador, V. H. _et al_. Cinnamic acid increases lignin production and inhibits

soybean root growth. _PLoS One_ 8, e69105 (2013). Article  ADS  CAS  Google Scholar  * Trivedi, P. _et al_. Huanglongbing alters the structure and functional diversity of microbial

communities associated with citrus rhizosphere. _The ISME journal_ 6, 363–383 (2012). Article  CAS  Google Scholar  * Qu, X. H. & Wang, J. G. Effect of amendments with different phenolic

acids on soil microbial biomass, activity, and community diversity. _Applied Soil Ecology_ 39, 172–179 (2008). Article  Google Scholar  * Zhang, S. _et al_. Baicalin released from

Scutellaria baicalensis induces autotoxicity and promotes soilborn pathogens. _Journal of chemical ecology_ 36, 329–338 (2010). Article  CAS  Google Scholar  * Zhou, X. & Wu, F.

p-Coumaric acid influenced cucumber rhizosphere soil microbial communities and the growth of Fusarium oxysporum f. sp. cucumerinum Owen. _PloS one_ 7, e48288 (2012). Article  ADS  CAS 

Google Scholar  * Zhou, X., Jia, H., Ge, X. & Wu, F. Effects of vanillin on the community structures and abundances of Fusarium and Trichoderma spp. in cucumber seedling rhizosphere.

_Journal of Plant Interactions_ 13, 45–50 (2018). Article  CAS  Google Scholar  * Jia, Z. H., Yi, J. H., Su, Y. R. & Shen, H. Autotoxic substances in the root exudates from continuous

tobacco cropping. _Allelopathy Journal_ 27, 87–96 (2011). Google Scholar  * Ali, M. B., Hahn, E. J. & Paek, K. Y. Effects of temperature on oxidative stress defense systems, lipid

peroxidation and lipoxygenase activity in Phalaenopsis. _Plant physiology and Biochemistry_ 43, 213–223 (2005). Article  CAS  Google Scholar  * Li, Z. H., Wang, Q., Ruan, X., Pan, C. D.

& Jiang, D. A. Phenolics and plant allelopathy. _Molecules_ 15, 8933–8952 (2010). Article  CAS  Google Scholar  * Baziramakenga, R., Leroux, G. D. & Simard, R. R. Effects of benzoic

and cinnamic acids on membrane permeability of soybean roots. _Journal of Chemical Ecology_ 21, 1271–1285 (1995). Article  CAS  Google Scholar  * Blum, U. & Gerig, T. M. Relationships

between phenolic acid concentrations, transpiration, water utilization, leaf area expansion, and uptake of phenolic acids: nutrient culture studies. _Journal of Chemical Ecology_ 31,

1907–1932 (2005). Article  CAS  Google Scholar  * Gauri, S. S., Mandal, S. M., Dey, S. & Pati, B. R. Biotransformation of p-coumaric acid and 2,4-dichlorophenoxy acetic acid by

Azotobacter sp. strain SSB81. _Bioresource technology_ 126, 350–353 (2012). Article  CAS  Google Scholar  * Wang, Y. _et al_. Environmental behaviors of phenolic acids dominated their

rhizodeposition in boreal poplar plantation forest soils. _Journal of soils and sediments_ 16, 1858–1870 (2016). Article  CAS  Google Scholar  * Tagawa, M., Tamaki, H., Manome, A., Koyama,

O. & Kamagata, Y. Isolation and characterization of antagonistic fungi against potato scab pathogens from potato field soils. _FEMS microbiology letters_ 305, 136–142 (2010). Article 

CAS  Google Scholar  * Carlsen, S. C. _et al_. Allelochemicals in rye (Secale cereale L.): cultivar and tissue differences in the production of benzoxazinoids and phenolic acids. _Natural

product communications_ 4, 199–208 (2009). Article  ADS  CAS  Google Scholar  * Gulzar, A. & Siddiqui, M. B. Root-mediated allelopathic interference of bhringraj (Eclipta alba L.) Hassk.

on peanut (Arachis hypogaea) and mung bean (Vigna radiata). _Applied Soil Ecology_ 87, 72–80 (2015). Article  Google Scholar  * Wang, H. B. _et al_. Characterization of metaproteomics in

crop rhizospheric soil. _Journal of proteome research_ 10, 932–940 (2010). Article  Google Scholar  * Wu, L. _et al_. Comparative metaproteomic analysis on consecutively Rehmannia

glutinosa-monocultured rhizosphere soil. _Plos one_ 6, e20611 (2011). Article  ADS  CAS  Google Scholar  * Li, Y. C. _et al_. Variations of rhizosphere bacterial communities in tea (Camellia

sinensis L.) continuous cropping soil by high‐throughput pyrosequencing approach. _Journal of applied microbiology_ 121, 787–799 (2016). Article  CAS  Google Scholar  * Yang, L. _et al_.

Duration of continuous cropping with straw return affects the composition and structure of soil bacterial communities in cotton fields. _Canadian journal of microbiology_ 64, 167–181 (2017).

Article  Google Scholar  * Tan, Y. _et al_. Diversity and composition of rhizospheric soil and root endogenous bacteria in Panax notoginseng during continuous cropping practices. _Journal

of basic microbiology_ 57, 337–344 (2017). Article  CAS  Google Scholar  * Chen, W. _et al_. Mechanisms by which organic fertilizer and effective microbes mitigate peanut continuous cropping

yield constraints in a red soil of south China. _Applied Soil Ecology_ 128, 23–34 (2018). Article  Google Scholar  * Chen, L. H. _et al_. Application of Trichoderma harzianum SQR‐T037

bio‐organic fertiliser significantly controls Fusarium wilt and affects the microbial communities of continuously cropped soil of cucumber. _Journal of the Science of Food and Agriculture_

92, 2465–2470 (2012). Article  CAS  Google Scholar  * Liu, X. _et al_. Microbial community diversities and taxa abundances in soils along a seven-year gradient of potato monoculture using

high throughput pyrosequencing approach. _PLoS One_ 9, e86610 (2014). Article  ADS  Google Scholar  * Chomel, M. _et al_. Secondary metabolites of Pinus halepensis alter decomposer organisms

and litter decomposition during afforestation of abandoned agricultural zones. _Journal of Ecology_ 102, 411–424 (2014). Article  Google Scholar  * Kefeli, V. I., Kalevitch, M. V. &

Borsari, B. Phenolic cycle in plants and environment. _J. Cell Mol. Biol_ 2, 13–18 (2003). Google Scholar  * Ma, B., Lv, X., Cai, Y., Chang, S. X. & Dyck, M. F. Liming does not

counteract the influence of long-term fertilization on soil bacterial community structure and its co-occurrence pattern. _Soil Biology & Biochemistry_ 123, 45–53 (2018). Article  CAS 

Google Scholar  * Blum, U. & Shafer, S. R. Microbial populations and phenolic acids in soil. _Soil Biology and Biochemistry_ 20, 793–800 (1988). Article  CAS  Google Scholar  * Siles, J.

A. & Margesin, R. Abundance and diversity of bacterial, archaeal, and fungal communities along an altitudinal gradient in alpine forest soils: what are the driving factors? _Microbial

ecology_ 72, 207–220 (2016). Article  Google Scholar  * Liu, J., Wang, X., Zhang, T. & Li, X. Assessment of active bacteria metabolizing phenolic acids in the peanut (Arachis hypogaea

L.) rhizosphere. _Microbiological research_ 205, 118–124 (2017). Article  CAS  Google Scholar  * Wu, H. _et al_. Mixed phenolic acids mediated proliferation of pathogens Talaromyces helicus

and Kosakonia sacchari in continuously monocultured Radix pseudostellariae rhizosphere soil. _Frontiers in Microbiology_ 7, 335 (2016). PubMed  PubMed Central  Google Scholar  * Han, Z.,

Sun, J., Lv, A. & Wang, A. Biases from different DNA extraction methods in intestine microbiome research based on 16S rDNA sequencing: a case in the koi carp, Cyprinus carpio var. Koi.

_MicrobiologyOpen_ 8, e00626 (2019). Article  Google Scholar  * Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. _EMBnet. journal_ 17, 10–12 (2011).

Article  Google Scholar  * Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. _Bioinformatics_ 27, 2194–2200

(2011). Article  CAS  Google Scholar  * Haas, B. J. _et al_. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. _Genome research_ 21, 494–504

(2011). Article  CAS  Google Scholar  * Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. _Nature methods_ 10, 996 (2013). Article  CAS  Google Scholar  *

Quast, C. _et al_. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. _Nucleic acids research_ 41, D590–D596 (2012). Article  Google Scholar  *

Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. _Nucleic acids research_ 32, 1792–1797 (2004). Article  CAS  Google Scholar  Download references

ACKNOWLEDGEMENTS We are grateful to Beijing Novogene Bioinformatics Technology Co., Ltd., for the 16S rRNA deep sequencing. This research was supported by the National Science Foundation of

China (31860357), Yunnan Graduate Scholar Newcomer Award (A2008057), and the China Tobacco Yunnan Industrial Co., Ltd (2017YN07). We would like to thank LetPub (www.letpub.com) for providing

linguistic assistance during the preparation of this manuscript. AUTHOR INFORMATION Author notes * Yuxiang Bai and Ge Wang contributed equally. AUTHORS AND AFFILIATIONS * Yunnan

Agricultural University, Kunming, 650201, China Yuxiang Bai, Ge Wang, Yadong Cheng, Puyou Shi, Chengcui Yang & Huanwen Yang * Yunnan Academy of Tobacco Agricultural Science, Kunming,

650021, China Zhaoli Xu Authors * Yuxiang Bai View author publications You can also search for this author inPubMed Google Scholar * Ge Wang View author publications You can also search for

this author inPubMed Google Scholar * Yadong Cheng View author publications You can also search for this author inPubMed Google Scholar * Puyou Shi View author publications You can also

search for this author inPubMed Google Scholar * Chengcui Yang View author publications You can also search for this author inPubMed Google Scholar * Huanwen Yang View author publications

You can also search for this author inPubMed Google Scholar * Zhaoli Xu View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS H.Y. and Z.X.

conceived the study and contributed to the design and interpretation of the research. Y.B. and G.W. carried out the experiments, analyzed the data, and wrote the manuscript. Y.B., P.S., C.Y.

and Y.C. contributed to the field experiments and collected samples. G.W. modified the manuscript. All authors read and approved the final manuscript. CORRESPONDING AUTHORS Correspondence

to Huanwen Yang or Zhaoli Xu. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER’S NOTE: Springer Nature remains neutral

with regard to jurisdictional claims in published maps and institutional affiliations. RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0

International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the

source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative

Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by

statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Bai, Y., Wang, G., Cheng, Y. _et al._ Soil acidification in continuously cropped

tobacco alters bacterial community structure and diversity via the accumulation of phenolic acids. _Sci Rep_ 9, 12499 (2019). https://doi.org/10.1038/s41598-019-48611-5 Download citation *

Received: 03 April 2019 * Accepted: 08 August 2019 * Published: 29 August 2019 * DOI: https://doi.org/10.1038/s41598-019-48611-5 SHARE THIS ARTICLE Anyone you share the following link with

will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt

content-sharing initiative