Physiological and transcriptomic analyses to reveal underlying phenolic acid action in consecutive monoculture problem of Polygonatum odoratum
BMC Plant Biology volume 21, Article number: 362 (2021)
The root rot of fragrant solomonseal (Polygonatum odoratum) has occurred frequently in the traditional P. odoratum cultivating areas in recent years, causing a heavy loss in yield and quality. The phenolic acids in soil, which are the exudates from the P. odoratum root, act as allelochemicals that contribute to the consecutive monoculture problem (CMP) of the medicinal plant. The aim of this study was to get a better understanding of P. odoratum CMP.
The phenolic acid contents, the nutrient chemical contents, and the enzyme activities related to the soil nutrient metabolism in the first cropping (FC) soil and continuous cropping (CC) soil were determined, and the differentially expressed genes (DEGs) related to the regulation of the phenolic acids in roots were analyzed. The results showed that five low-molecule-weight phenolic acids were detected both in the CC soil and FC soil, but the phenolic acid contents in the CC soil were significantly higher than those in the FC soil except vanillic acid. The contents of the available nitrogen, available phosphorus, and available potassium in the CC soil were significantly decreased, and the activities of urease and sucrase in the CC soil were significantly decreased. The genomic analysis showed that the phenolic acid anabolism in P. odoratum in the CC soil was promoted. These results indicated that the phenolic acids were accumulated in the CC soil, the nutrient condition in the CC soil deteriorated, and the nitrogen metabolism and sugar catabolism of the CC soil were lowered. Meantime, the anabolism of phenolic acids was increased in the CC plant.
The CC system promoted the phenolic acid anabolism in P. odoratum and made phenolic acids accumulate in the soil.
Polygonatum odoratum (Mill.) Druce, popularly known as fragrant solomonseal, is a traditional Chinese perennial medicinal plant mainly cultivated in the southern parts of China and the other Southeast Asian countries such as Thailand and Vietnam. As a medicinal and edible plant, its rhizome has the functions of removing dryness, promoting secretion, and quenching thirst [1, 2], and it is widely welcomed especially in the Southeast Asian markets. However, the root rot of P. odoratum has occurred frequently in the planting areas in recent years due to the long-term continuous cropping (CC), causing the serious consecutive monoculture problem (CMP) with a sharp decline both in yield and quality. For example, in Shaodong County of Hunan Province, which was the largest traditional high-quality P. odoratum planting area in China, the medicinal peasants have almost given up the cultivation of P. odoratum due to this serious replant disease.
Researches have revealed that CMP is mainly caused by three factors covering the imbalance of soil nutrients, the shift in microorganisms towards pathogenic and allelopathic autotoxicity, and the root exudates of secondary metabolites, among which the root exudates are likely to be considered as the main factor to induce the root rot by influencing rhizosphere microbes and soil nutrients . The root exudates include low-molecule-weight substances such as polysaccharides, vitamins, nucleotides, and phenolic acids. Among them, phenolic acids are most likely to be considered as the allelochemicals that can change membrane permeability, inhibit nutrient uptake, and inactivate plant endogenous hormones to influence the normal physiological process .
Phenolic acids are a kind of substance containing active organic acids on their aromatic rings that are mainly synthesized by the shikimic acid pathway which provides phenylpropanoids or synthesized by the polyketide pathway in plants which provides simple phenolic acids directly. Shikimic acid is synthesized from 4-phosphate erythritose via the phosphoenolpyruvic acid and pentose phosphate pathway, and then the aromatic amino acids are subsequently used as precursors for the synthesis of phenolic acids . Phenolic acids in plants are products responding to environmental stresses, and they may be involved in the intracellular and intercellular signaling processes as signal transporting substances . When plant roots are infected, the phenolic acid levels in the root cells are promoted to act as antioxidants and inhibit infectious microorganisms directly .
Although phenolic acids are of vital importance for the plant’s interaction with the environment, few documents have been published on the relationship between the P. odoratum phenolic acids and the root rot disease. An early study showed that the low-molecule-weight substances extracted from the P. odoratum rhizosphere soil solutions exhibited an obvious self-poison effect on P. odoratum seedling growth and allelopathy effect on the rhizome spore germination rate . Our recent studies have revealed that the Polygonatum CMP could be alleviated or eradicated under the rotation regimes or with an aeroponic system [9, 10], and the miRNAs-regulated genes participating in the sugar and phenylpropanoid metabolism have been enriched through the mircoRNA sequencing on the CC and the first cropping (FC) P. odoratum . All these results prompted that the root rot in P. odoratum was likely related to the phenolic acid metabolism. Other researchers revealed that the root exudates including phenolic acids may not only involve in the signal transduction and hormone synthesis but also modify the rhizosphere soil environment to regulate the microbial community in the rhizosphere in some other medicinal plants such as Rehmannia glutinosa [5, 12]. Besides, recent researches have provided more information on the complex correlation among the root rot, phenolic acids, and CMP. Although phenolic acids in the CC soils stimulate the propagation of Fusarium which directly results in root rot, single phenolic acid also plays a role as a double-edged sword to decrease the Fusarium growth and production at a lower level, and the CC system can strengthen the phenolic acid accumulation in the rhizosphere soils of peanut, Panax notoginseng, and Chrysanthemum morifolium [13,14,15].
In this study, the low-molecule-weight phenolic acids from the rhizosphere soil of FC and CC P. odoratum were identified, and the relative chemical and physiological indices for evaluating the influences of the phenolic acids on the soil were also determined. Furthermore, the transcriptome profiles from the FC and CC rhizomes of P. odoratum were compared through high-throughput sequencing, and then the differentially expressed genes (DEGs) were analyzed. Then the pathway related to the phenolic acid metabolism was analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG). This research might provide more information to understand the regulatory function of the phenolic acids in P. odoratum under the CC system.
Rhizosphere soil phenolic acids
The rhizosphere soil phenolic acids were determined using the high-performance liquid chromatography (HPLC) method as compared to the standard samples of p-hydroxybenzoic acid, vanillic acid, syringic acid, cumaric acid, and ferulic acid with the retention time of 15.42, 19.29, 21.84, 27.16, and 28.94 min, respectively (Fig. S1A), and the above five phenolic acids were detected both in the FC and CC rhizosphere soil (Fig. S1B, Fig. S1C). R2 of the five regression curves was 0.999 (Table 1). The good linearity indicated that the determination method was reliable. The phenolic acid contents in the FC and CC rhizosphere soil were calculated according to the equations in Table 1. As shown in Table 2, p-hydroxybenzoic, syringic acid, cumaric acid, and ferulic acid contents in the CC soil were significantly higher than those in the FC soil, whereas vanillic acid in the CC and FC soil had no significant difference, although vanillic acid contents in the CC and FC soil exhibited much higher level than those of the other phenolic acids. The results showed that the phenolic acid contents in the CC soil were higher than those in the FC soil.
Rhizosphere soil chemical properties and enzyme activities
The rhizosphere soil elements of the CC and FC system reflecting the soil nutritional status were determined in this experiment. As shown in Table 3, the contents of total phosphorus (TP), total potassium (TK), available nitrogen (AN), available phosphorus (AP), available potassium (AK) in the CC soil were significantly lower than those in the FC soil, whereas the total nitrogen (TN) contents between the CC soil and FC soil had no significant change. Although TN in the CC soil showed no significant difference as compared to that in the FC soil, AN, which is the form of nitrogen that can be absorbed directly by plants, decreased significantly, and the other nutrient chemicals also decreased significantly. These results insisted that the CC soil nutritional status had deteriorated.
The rhizosphere soil enzymes related to the soil nutrient metabolism level were also determined in this experiment. As shown in Table 4, the activities of polyphenol oxidase (PPO), catalase (CAT), and acid phosphatase (ACP) had no significant differences between the CC soil and FC soil, whereas the activities of urease (UE) and sucrase (SC) in the CC soil were significantly lower than those in the FC soil. The results indicated that the soil oxidative metabolism and phosphorus utilization had not been changed between the CC soil and FC soil, but the levels of the nitrogen metabolism and sugar catabolism of the CC soil were lowered as compared to the FC soil.
RNA-seq analysis and qRT-PCR verification
After filtering out the low-quality sequences, 299,730,812 clean reads, which were accounted for 95.38% of the raw reads, were obtained from the samples of the CC and FC P. odoratum roots, and Q20 and Q30 values and GC contents were all in the reliable ranges (Table S1). Owing to the absence of the reference genome, the total clean reads were spliced to obtain 842,213 transcripts and 510,970 unigenes, and the N50 lengths of the transcripts and unigenes were 931 and 1167 bp, respectively, and 298,939 unigenes were annotated in the databases including NR, NT, KO, SwissProt, PFAM, GO, and KOG, whereas 176,156 unigenes were annotated in the GO database. Then a total number of 15,788 DEGs were screened out in the CC vs FC root tissues including 4843 up-regulated DEGs and 10,945 down-regulated DEGs (padj < 0.05).
These DEGs were significantly enriched in 307 GO functional items. The biological processes included “oxidation-reduction process”, “carbohydrate metabolic process”, “polysaccharide metabolic process”, etc.; the molecular functions included “catalytic activity”, “oxidoreductase activity”, “hydrolase activity”, etc.; the cell components included “cell wall”, “cell periphery”, “cytoskeletal part”, etc. (Fig. 1A). The up-regulated DEGs were mainly enriched in the metabolic process, the single-organism metabolic process, the catalytic activity, the oxidoreductase activity, etc. (Fig. 1B-D); the down-regulated DEGs were mainly enriched in the oxidation-reduction process, the carbohydrate metabolic process, the hydrolase activity (acting on glycosyl bonds), the hydrolase activity, the external encapsulating structure, the cell periphery, etc. (Fig. 1E-G). The results showed that the gene expressions of the CC plant were lower than those of the FC plant on the whole level, but the responses of P. odoratum to CMP were regulated by complicated networks, and lots of DEGs took part in the processes such as the secondary metabolism and the response reaction. It should be mentioned that even the up-regulated DEGs and the down-regulated DEGs co-existed in the same processes. For example, although 665 DEGs were significantly down-regulated in the carbohydrate metabolic process, the number of significantly up-regulated DEGs was also up to 177 (Table S2).
The KEGG enrichment analysis showed that DEGs were mainly enriched in the metabolic and biosynthetic pathways including phenylpropanoid biosynthesis, linoleic acid metabolism, phenylalanine metabolism, and starch and sucrose metabolism pathway, etc. The up-regulated DEGs were significantly enriched in the pathways such as “phenylalanine, tyrosine and tryptophan biosynthesis”, “phenylalanine metabolism”, “phenylpropanoid biosynthesis”, and “nitrogen metabolism” (Fig. 2A), whereas the down-regulated DEGs were significantly enriched in the pathways such as “plant hormone signal transduction”, “DNA replication”, “brassinosteroid biosynthesis”, and “sesquiterpenoid and triterpenoid biosynthesis” (Fig. 2B).
To eliminate those DEGs with relatively lower expressions, log2foldchange > 1.5 was selected as the screening condition, the pathways closely relative to the phenolic acid metabolism were enriched including “pentose phosphate pathway”, “pentose and glucoronate interconversions”, “starch and sucrose metabolism”, “tyrosine metabolism, and phenylalanine”, “tyrosine, and tryptophan biosynthesis” (Fig. S2), and then the significantly enriched DEGs were found including FBP, PGD, rpe, GPD, SORD, malZ, glgC, E22.214.171.124, aroK, aroC, TYDC, and E126.96.36.199 (Table 5). These DEGs were classified into two groups: down-regulated DEGs and up-regulated DEGs. Among them, FBP, PGD, rpe, SORD, malZ, aroK, aroC, and TYDC were up-regulated, whereas glgC, GPI, E188.8.131.52, and E184.108.40.206 were down-regulated (Fig. 3). Among these genes, it was found that the genes that encoded the enzymes in favor of the phenolic acid synthesis were up-regulated, whereas those unfavorable for the phenolic acid synthesis were down-regulated.
To verify the above results, a total of 12 up-regulated DEGs in the CC vs FC plants related to the synthesis of phenolic acids were selected from Table 5 and were tested at the mRNA expression level. The expression levels of DEGs after normalization were presented, and then their mRNA expressions of these DEGs were verified using the quantitative real-time PCR (qRT-PCR). As shown in Fig. 4, the mRNA levels of 11 DEGs were increased except that of 60,288.238790. These results proved that the enriched pathways related to the phenolic acid anabolism through RNA-seq analysis in this study were reliable.
The soil phenolic acids are usually released by plant leaching, root exudation, or decomposition of plant residues [16, 17]. For the low-molecule-weight phenolic acids, there are two groups including benzoic acid derivatives and cinnamic acid derivatives, which play an important role in plant growth and interactions with the environment. These molecules reduce water utilization, hydraulic conductivity, and nutrient uptake to influence the plant’s physiological state . In this study, five low-molecule-weight phenolic acids were identified from the P. odoratum rhizosphere soil, i.e., p-hydroxybenzoic acid, vanillic acid, syringic acid, cumaric acid, and ferulic acid, where p-hydroxybenzoic acid, vanillic acid, and syringic acid were from benzoic acid derivatives, and cumaric acid and ferulic acid were from cinnamic acid derivatives. Although their total contents were very subtle in the soil, the phenolic acid contents in the CC soil were significantly higher than those in the FC soil. These results supported that the CC system might result in an accumulation of phenolic acids secreted by the medicinal plant.
The content level and balance of nitrogen, phosphorus, and potassium in soil are essential for plant growth and development . It was observed that TP, TK, AN, AP, and AK were significantly reduced in the CC soil. Subsequently, this phenomenon might lead to nutrient deterioration in the CC soil for medicinal plant growth. Soil AN is the nitrogen form which is easy to be absorbed and used directly by plants, mainly including NH+4-N, NO-3-N, amino acid, amide, and hydrolysable protein nitrogen . The early research had proved that nitrogen in soil inhibits the formation of phenolic acids in plants, and there was a negative correlation between the total phenolic acids and the soil nitrogen content according to the “carbon/nutrient balance hypothesis” . It was also observed this correlation in our experiment. The result showed that the AN in the CC soil was significantly lower than that in the FC soil, whereas the phenolic acids in the CC soil were at a higher level.
The phenolic acids are detrimental to the rhizosphere soil environment and can change the soil microbial community . Phenolic acids have a hydroxyl group and a carboxyl group. For this reason, they can hydrogen-bond with the soil enzymes and change the soil pH value to worsen the soil’s physiological status. The activities of soil enzymes can reflect the soil health status . Soil SC hydrolyzes sucrose into monosaccharides, which are closely related to soil nutrition, whereas soil UE hydrolyzes urea to produce ammonia, which is the source of AN. In this experiment, it was found that the activities of UE and SC in the CC soil were decreased as compared to those in the FC soil. Under this situation, the utilization of sugar and the production of AN were weakened in the CC soil, and then the lack of nitrogen could lead to a comparatively higher content of phenolic acids in the CC plant as compared to that in the FC plant. Consequently, more phenolic acids were secreted and finally accumulated in the CC soil.
The high-throughput sequencing and relative qRT-PCR analysis also revealed that the gene expressions of the enzymes that catalyzed the key steps in the synthesis of the phenolic acids were regulated in the P. odoratum root tissues. The pentose pathway and relative pathways are the major source for metabolic intermediates in plants  and can provide erythrose-4-phosphate as a precursor for shikimic acid. In the study, 12 kinds of DEGs that were closely relative to the shikimic acid and chorismate metabolisms were enriched, where eight of them were up-regulated, and four of them were down-regulated (Fig. 5).
Among the up-regulated genes, PDG, rpe, and SORD encoded 6-phosphogluconate dehydrogenase, ribulose-phosphate 3-epimerase, and L-iditol 2-dehydrogenase, respectively, and they were the enzymes for synthesis of ribulose-5-phosphate, xylulose-5-phosphate, and xylusose, respectively [24,25,26]. The up-regulation of these genes was beneficial to the xylusoe-5-phosphate biosynthesis. FBP and malZ encoded fructose-1, 6-bisphosphatase and α-glucosidase, respectively. The up-regulation of them was favorable for the synthesis of fructose-6-phosphate, which was material for the synthesis of erthrose-4-phosphate, a precursor of shikmic acid [27, 28]. aroK and aroC encoded shikimate kinase and chorismate synthase, respectively [29, 30], and the up-regulation of them was beneficial to chorismate biosynthesis, which was the precursor of the low-molecule-weight phenolic acids. On the contrary, the down-regulation of glgC, E220.127.116.11, and GPI, which encoded glucose-1-phosphate adenylyltransferase, endoglucanase, and glucose-6-phosphate isomerase, respectively, could inhibit the cellodexitrin and amylose metabolism [31,32,33]. This process in turn could provide more materials for the erthrose-4-phosphate biosynthesis. Early researches have proved that the phenolic acid biosynthesis in plants is triggered by biotic and abiotic stresses [34, 35]. Likewise, the results in our study also supported that the biosynthesis of phenolic acids in the CC plant proceeded through the activation of the primary pathway (shikimic acid) and the suppression of the relative branch pathway (cellodexitrin and amylose).
Tyramine is the product of tyrosine decarboxylation, and partial phenolic acids in plants can be conjugated with tyramine especially when the plant is under stress [36, 37]. The tyramine-conjugated phenolic compounds act as phytoalexins to make plants adapt to the external stress, and the biosynthetic processes including the shikimate, phenylpropanoid, and arylmonoamine pathways are activated . The conjugation is also favorable for more phenolic acids to accumulate. In this experiment, it was also found that the expression of tyrosine decarboxylase in the CC plant was up-regulated and polyphenol oxidase was down-regulated. It was speculated that the change might help to promote the phenolic acid anabolism in P. odoratum in the CC soil and made them accumulate in the soil.
The inducement of the root rot of P. odoratum in the CC soils was very complicated. The multiple factors to cause CMP may mainly include soil element imbalance, soil microbial population transforming to harmful population, and autotoxicity caused by the secretions from plant roots, but the comprehensive mechanism for CMP remains unclear . For P. odoratum in the CC soils, the nutrient composition and the soil enzymes promoted the synthesis of phenolic acids. On the other hand, there might be multiple gene families activated together to accelerate the phenolic acid synthesis process coordinately to influence the CMP formation.
In this study, five low-molecule-weight phenolic acids were identified from the P. odoratum rhizosphere soil. Among them, the phenolic acid contents in the CC soil were significantly higher than those in the FC soil except vanillic acid. The contents of AN, AP, and AK in the CC soil were significantly higher than those in the FC soil, and the nitrogen metabolism and sugar catabolism of the CC soil were lowered. The genomic analysis showed that the phenolic acid anabolism in P. odoratum in the CC soil was promoted. This research presented us with a preliminary understanding of the role of phenolic acids in the CMP of P. odoratum. In perspective, proteomics and metabolomics are needed to look into more metabolic mechanism details of phenolic acid metabolism on the root rot of P. odoratum.
Experiment materials and chemicals
The P. odoratum cultivar “Zhushiwei” was selected as the experimental material and authenticated by associate Professor Zefa Liu, a horticulturist from Loudi Agricultural Institute, Hunan Province, China. The roots and rhizosphere soil of P. odoratum were collected in May, 2017 from the Gutang Town Experimental Station of Loudi City, Hunan Province, China with the permission of Loudi Agricultural Institute, Hunan Province, China. For the CC system, P. odoratum was cultivated in the land where the same plants had been harvested; for the FC system, P. odoratum was cultivated on the same date in the land near CC where the cabbages had been harvested. The collection of rhizosphere soil was according to Riley and Barber’s method [39, 40]. Briefly, the soil with a complete root system was excavated in the selected plot, the large soil without root was shaken off gently, and then the soil adhered to the root circumference (0–5 mm from the root circumference) was collected with a brush as the rhizosphere soil. The rhizosphere soil was air-dried under room temperature, passed through the 2 mm sieve, and stored at 4 °C until use.
Methanol, p-hydroxybenzoic acid, vanillic acid, syringic acid, cumaric acid, ferulic acid, and methanol were of chromatographic purity, and the other chemicals were of analytical purity.
Rhizosphere soil phenolic acid determination
The rhizosphere soil of 1.0 g was added with 3.0 mL 1 N NaOH, and then shaken at 4 °C overnight. After centrifuged with 8000×g at 4 °C for 10 min, the supernatant was adjusted to pH 2.5 with HCl. After the liquid was extracted two times with 10 mL ethyl acetate, it was concentrated to dryness under nitrogen, and then dissolved in methanol. The solution was filtered with a 0.45 μm membrane for the HPLC analysis.
The HPLC conditions were as follows: volume of 10 μL was loaded on a Kromasil C18 column (250 mm × 4.6 mm × 5 μm) (Akzo Nobel, Amsterdam, Netherlands) with a Rigol HPLC L3000 (Rigol Technologies, Beijing, China). The gradient elution was with eluent A 1% phosphoric acid and eluent B methanol. The elution gradients were 80% A + 20% B on 0–10 min, 70% A + 30% B on 10–20 min, 50% A + 50% B on 20–30 min, 50% A + 50% B on 30–40 min, 80% A + 20% B on 40–45 min, and 80% A + 20% B on 45–55 min, respectively, with the rate of 1 mL/min at 30 °C, and the signal was detected at 280 nm with an ultraviolet detector. For the standard curve, p-hydroxybenzoic acid, vanillic acid, syringic acid, cumaric acid, and ferulic acid were used as standard samples and detected as above. The external standard method was used for the quantitative analysis of the phenolic acids as per Wang et al. .
Determination of rhizosphere soil chemical properties and enzyme activity
The TN content of the rhizosphere soil was determined with the Kjeldahl method, the TP content of the rhizosphere soil was determined with the molybdenum-blue colorimetry after the soil was extracted with the mixture of concentrated sulfuric acid and perchloric acid, the TK content of the rhizosphere soil was determined with the flame spectrophotometry after the soil was extracted with the mixture of concentrated sulfuric acid and perchloric acid, the AN content of the rhizosphere soil was determined with the diffusion method after hydrolysis with the alkaline solution, the AP content of the rhizosphere soil determined with the molybdenum-blue colorimetry after the soil was extracted with sodium hydroxide, and the AK content of the rhizosphere soil was determined with the flame spectrophotometry after the soil was extracted with sodium acetate . The PPO activity of the rhizosphere soil was determined with the catechol method . The CAT activity of the rhizosphere soil was determined as per Kraus and Fletcher . The UE activity of the rhizosphere soil was determined with the indophenol blue colorimetric method . The ACP of the rhizosphere soil activity was determined as per Saa et al. . The SC activity of the rhizosphere soil was determined as per Li and Lu . All of these analyses were assigned using the kits of Suzhou Comin Biotechnology Co. Ltd. and were repeated three times at the biological sample level.
RNA-seq and bioinformatics analysis of root
The root samples of CC and FC plants were obtained randomly from the fields during the rhizome expansion stage, and each tissue had three biological repeats in total. The library preparation and sequencing were performed in the Novogene Technology Co., Ltd. (Beijing, China). The total RNA was extracted using the Invitrogen Kit (Thermo Fisher Scientific, Beijing, China) and detected with a nanophotometer (Implen Inc., CA, USA). The sequencing libraries were generated using NEBNext Ultra RNA Library Prep Kit for Illumina (New England Biolabs, MA, USA). The libraries were sequenced on Illumina Hiseq Platform and then 150 bp paired-end reads were generated.
The clean reads were obtained from the raw pair-end reads of mRNA sequencing data by removing adapters, low-quality reads, and the reads in which bases could not be determined, and then Trinity was employed to splice the clean-reads to obtain the de novo transcriptome assembly . For DEGs analysis, DESeq R package(1.10.1)was used to screen the differential expression (n = 3) with the gene symbol annotation (padj < 0.05) . The KEGG pathway enrichment analysis was performed using KOBAS (2.0) and the pathway with false discovery rate (FDR) ≤ 0.05 was defined as the significantly enriched pathway with DEGs .
The total RNA from the P. odoratum root was extracted with Trizol reagent according to the manufacture’s instruction. To obtain cDNA, the total RNA (0.5 μg) was reversely transcribed using a TUREscript cDNA Synthesis Kit (Aidlab Biotechnologies Co., Ltd., Beijing, China) with oligo dT as a primer. The candidate genes and the reference gene were analyzed with a qTOWER2.2 real-time PCR (Analytik Jena AG, Jena, Thuringia, Germany) with the primers in Table S3. The protocol for qRT-PCR was as follows: 3 min at 95 °C, followed by 10 s at 95 °C and 30 s at 58 °C with 39 cycles, and then was ended after a melt curve analysis from 60 °C to 95 °C with an increment of 1 °C for 4 s at each step. Three biological replicates were used for each analysis and the data were analyzed using the 2−ΔΔCt method.
All the experiments were repeated three times at the biological sample level, and the differences were analyzed using IBM SPSS 19.0 (International Business Machines Corporation, Armonk, NY, USA) with the T-student’s method (P < 0.05).
Availability of data and materials
All data supporting the findings were contained in the manuscript and its supplementary files except the RNA-seq raw data. And all the RNA-seq raw data were uploaded in the SRA of NCBI (PRJNA507291).
Consecutive monoculture problem
Differentially expressed genes
The first cropping
The Kyoto Encyclopedia of Genes and Genomes
High-performance liquid chromatography
quantitative real-time PCR
Chen Y, Yin L, Zhang X, Wang Y, Chen Q, Jin C, et al. Optimization of alkaline extraction and bioactivities of polysaccharides from rhizome of Polygonatum odoratum. Biomed Res Int. 2014;2014:504896–8. https://doi.org/10.1155/2014/504896.
Lan G, Chen H, Wang Z, Zhang W, Zhang L. Extraction of Polygonatum odoratum polysaccharides using response surface methodology and preparation of a compound beverage. Carbohydr Polym. 2011;86(3):1175–80. https://doi.org/10.1016/j.carbpol.2011.06.009.
Li M, Yang Y, Feng F, Zhang B, Yang C, Gu L, et al. Differential proteomic analysis of replanted Rehmannia glutinosa roots by iTRAQ reveals molecular mechanisms for formation of replant disease. BMC Plant Biol. 2017;17(1):116. https://doi.org/10.1186/s12870-017-1060-0.
Li Z, Wang Q, Ruan X, Pan C, Jiang D. Phenolics and plant allelopathy. Molecules. 2010;15(12):8933–52. https://doi.org/10.3390/molecules15128933.
Marchiosi R, dos Santos WD, Constantin RP, de Lima RB, Soares AR, Finger-Teixeira A, et al. Biosynthesis and metabolic actions of simple phenolic acids in plants. Phytochem Rev. 2020;11(4):865–906. https://doi.org/10.1007/s11101-020-09689-2.
Cheynier V, Comte G, Davies KM, Lattanzio V, Martens S. Plant phenolics: recent advances on their biosynthesis, genetics, and ecophysiology. Plant Physiol Biochem. 2013;72:1–20. https://doi.org/10.1016/j.plaphy.2013.05.009.
Weiß S, Winkelmann T. Transcriptome profiling in leaves representing aboveground parts of apple replant disease affected Malus domestica ‘M26’ plants. Sci Hortic. 2017;222:111–25. https://doi.org/10.1016/j.scienta.2017.05.012.
Xiao L, Liu Y, Xu B, Yu Y, Liu T. Study on allelopathy and identification of allelochemical in rhizosphere soil of Polygonatum. Chin Agric Sci Bulletin. 2015;31:163–7. https://doi.org/10.11924/j.issn.1000-6850.casb15070050.
Chen Y, Zhang X, Hu Y, Jin C. Effects of soil nutrients, enzyme activities and microbial composition in vegetable-Polygonatum Odoratum rotation regimes. Rev Chim. 2020;71(6):460–72. https://doi.org/10.37358/RC.20.6.8212.
Zou T, Jin C, Chen Y, Liu Z, Hu Y. Responses of Polygonatum odoratum seedlings in aeroponic culture to treatments of different ammonium: nitrate ratios. J Plant Nutr. 2017;40(20):2850–61. https://doi.org/10.1080/01904167.2017.1384011.
Chen Y, Liu Z, Tu N, Hu Y, Jin C, Luo Y, et al. Integrated transcriptome and microRNA profiles analysis reveals molecular mechanisms underlying the consecutive monoculture problem of Polygonatum odoratum. Cell Mol Biol. 2020;66(2):47–52. https://doi.org/10.14715/cmb/2020.66.2.7.
Yang Y, Li M, Li X, Chen X, Lin W, Zhang Z. Transcriptome-wide identification of the genes responding to replanting disease in Rehmannia glutinosa L. roots. Mol Biol Rep. 2015;42(5):881–92. https://doi.org/10.1007/s11033-014-3825-y.
Zhao Y, Cheng Y, Ma Y, Chen C, Xu F, Dong X. Role of phenolic acids from the rhizosphere soils of Panax notoginseng as a double-edge sword in the occurrence of root-rot disease. Molecules. 2018;23(4):819. https://doi.org/10.3390/molecules23040819.
Li P, Chen J, Li Y, Zhang K, Wang H. Possible mechanisms of control of Fusarium wilt of cut chrysanthemum by Phanerochaete chrysosporium in continuous cropping fields: a case study. Sci Rep. 2017;7(1):15994. https://doi.org/10.1038/s41598-017-16125-7.
Lin P, Zhao H, Tang Z, Zhang Y, Lin H, Shen Y, et al. Effects of continuous cropping on root exudates of different resistance peanut (Arachis hypogaea L.) varieties and allelochemicals content in soil. Chin J Oil Crop Sci. 2015;37:467–74. https://doi.org/10.7505/j.issn.1007-9084.2015.04.006.
Chen ZL, Krishnamurti GSR, Naidu R. Separation of phenolic acids in soil and plant tissue extracts by co-electroosmotic capillary electrophoresis with direct UV detection. Chromatographia. 2001;53(3-4):179–84. https://doi.org/10.1007/BF02491567.
Tian G, Bi Y, Sun Z, Zhang L. Phenolic acids in the plow layer soil of strawberry fields and their effects on the occurrence of strawberry anthracnose. Eur J Plant Pathol. 2015;143(3):581–94. https://doi.org/10.1007/s10658-015-0711-y.
Cecchi AM, Koskinen WC, Cheng HH, Haider K. Sorption–desorption of phenolic acids as affected by soil properties. Biol Fertil Soils. 2004;39(4):235–42. https://doi.org/10.1007/s00374-003-0710-6.
Lambers H, Chapin FS, Pons TL. Plant physiological ecology. New York: Springer; 1998. https://doi.org/10.1007/978-1-4757-2855-2.
Li W, Yang G, Chen H, Tian J, Zhang Y, Zhu Q, et al. Soil available nitrogen, dissolved organic carbon and microbial biomass content along altitudinal gradient of the eastern slope of Gongga Mountain. Acta Ecol Sin. 2013;33(5):266–71. https://doi.org/10.1016/j.chnaes.2013.07.006.
Muzika R. Terpenes and phenolics in response to nitrogen fertilization: a test of the carbon/nutrient balance hypothesis. Chemoecology. 1993;4(1):3–7. https://doi.org/10.1007/BF01245890.
Alkorta I, Aizpurua A, Riga P, Albizu I, Amézaga I, Garbisu C. Soil enzyme activities as biological indicators of soil health. Rev Environ Health. 2003;18(1):65–73. https://doi.org/10.1515/reveh.2003.18.1.65.
Zamboni N, Fischer E, Laudert D, Aymerich S, Hohmann H, Sauer U. The Bacillus subtilis yqjl gene encodes the NADP+-dependent 6-p-gluconate dehydrogenase in the pentose phosphate pathway. J Bacteriol. 2004;186(14):4528–34. https://doi.org/10.1128/JB.186.14.4528-4534.2004.
Pickl A, Schönheit P. The oxidative pentose phosphate pathway in the haloarchaeon Haloferax volcanii involves a novel type of glucose-6-phosphate dehydrogenase - the archaeal Zwischenferment. FEBS Lett. 2015;589(10):1105–11. https://doi.org/10.1016/j.febslet.2015.03.026.
Teige M, Melzer M, Süss K. Purification, properties and in situ localization of the amphibolic enzymes D-ribulose 5-phosphate 3-epimerase and transketolase from spinach chloroplasts. Eur J Biochem. 1998;252(2):237–44. https://doi.org/10.1046/j.1432-1327.1998.2520237.x.
Ye R, Wong S. Transcriptional regulation of the Bacillus subtilis glucitol dehydrogenase gene. J Bacteriol. 1994;176(3):3314–20. https://doi.org/10.1007/BF02181950.
Roy D, Ward P. Comparison of fructose-1, 6-bisphosphatase gene (fbp) sequences for the identification of lactobacillus rhamnosus. Curr Microbiol. 2004;49(5):313–20. https://doi.org/10.1007/s00284-004-4355-y.
Bauer S, Vasu P, Persson S, Mort AJ, Somerville C. Development and application of a suite of polysaccharide-degrading enzymes for analyzing plant cell walls. PNAS. 2006;103(30):11417–22. http://www.pnas.org/cgi/doi/10.1073/pnas.0604632103.
Oliveira JS, Pinto CA, Basso LA, Santos DS. Cloning and overexpression in soluble form of functional shikimate kinase and 5-enolpyruvylshikimate 3-phosphate synthase enzymes from mycobacterium tuberculosis. Protein Expr Purif. 2001;22(3):430–5. https://doi.org/10.1006/prep.2001.1457.
Charles IG, Lamb HK, Pickard D, Dougan G, Harwkins AR. Isolation, characterization and nucleotide sequences of the aroC genes encoding chorismate synthase from Salmonella typhi and Escherichia coli. J Gen Microbiol. 1990;136(2):353–8. https://doi.org/10.1099/00221287-136-2-353.
Leung PSC, Preiss J. Biosynthesis of bacterial glycogen: primary structure of Salmonella typhimurium ADPglucose synthetase as deduced from the nucleotide sequence of the glgC gene. J Bacteriol. 1987;169(1):4355–60. https://doi.org/10.1083/jcb.107.1.101.
Yoshida K, Komae K. A rice family 9 glycoside hydrolase isozyme with broad substrate specificity for hemicelluloses in yype II cell walls. Plant Cell Physiol. 2006;47(11):1541–54. https://doi.org/10.1093/pcp/pcl020.
Froman BE, Tait RC, Gottlieb LD. Isolation and characterization of the phosphoglucose isomerase gene from Escherichia coli. Mol Gen Genet. 1989;217(1):126–31. https://doi.org/10.1007/BF00330951.
Mandal ST, Chakraborty D, Dey S. Phenolic acids act as signaling molecules in plant-microbe symbioses. Plant Signal Behav. 2010;5(4):359–68. https://doi.org/10.4161/psb.5.4.10871.
Becerra-Moreno A, Redondo-Gil M, Benavides J, Nair V, Cisneros-Zevallos L, Jacobo-Velázquez. Combined effect of water loss and wounding stress on gene activation of metabolic pathways associated with phenolic biosynthesis in carrot. Front Plant Sci. 2015;6:837. https://doi.org/10.3389/fpls.2015.00837.
Li Y, Wang C, Wang F, Dong H, Guo S, Yang JX. Phenolic components and flavanones from Dendrobium andidum. Chin Pharm J. 2010;45:975–9.
Cho M, Lee S. Phenolic phytoalexins in rice: biological functions and biosynthesis. Int J Mol Sci. 2015;16(12):29120–33. https://doi.org/10.3390/ijms161226152.
Park LH, Lee S, Jung K, Hahn T, Cho M. Transcriptomic analysis of UV-treated rice leaves reveals UV-induced phytoalexin biosynthetic pathways and their regulatory networks in rice. Phytochem. 2013;96:57–71. https://doi.org/10.1016/j.phytochem.2013.08.012.
Riley D, Barber SA. Bicarbonate accumulation and pH changes at the soybean (Glycine max (L.) Merr.) root-soil interface. Soil Sci Soc Am J. 1969;33(6):905–8. https://doi.org/10.2136/sssaj1969.03615995003300060031x.
Riley D, Barber SA. Salt accumulation at the soybean (Glycine max ( L.) Merr.) root-soil interface. Soil Sci Soc Am J. 1970;34(1):154–5. https://doi.org/10.2136/sssaj1970.03615995003400010042x.
Wang J, Zhang L, Guo Q, Kou J, Yu B, Gu D. Quantitative analysis of seven phenolic acids in eight Yinqiao Jiedu serial preparations by quantitative analysis of multi-components with single-marker. Acta Pharm Sin. 2015;50(4):480–5. https://doi.org/10.16438/j.0513-4870.2015.04.001.
Zheng B. Technical guide for soil analysis. Beijing: China Agricultural Press; 2012. p. 154–212.
Arogba SS, Ajiboye OL, Ugboko LA, Essienette SY, Afolabi PO. Properties of polyphenol oxidase in mango (Mangifera indica) kernel. J Sci Food Agric. 1998;77(4):459–62. https://doi.org/10.1002/(SICI)1097-0010(199808)77:4<459::AID-JSFA61>3.0.CO;2-O.
Kraus TE, Fletcher RA. Paclobutrazol protect wheat seedlings from heat and paraquat injury. Is detoxification of active oxygen involved? Plant Cell Physiol. 1994;35(1):45–52. https://doi.org/10.1159/000163877.
Das S, Ganguly D, Mukherjee A, Chakraborty S, De TK. Soil urease activity of Sundarban mangrove ecosystem, India. Adv Microbi. 2017;7(8):617–32. https://doi.org/10.4236/aim.2017.78048.
Saa A, Trasar-Cepeda MC, Gil-Sotres F, Carballas T. Changes in soil phosphorus and acid phosphatase activity immediately following forest fires. Soil Biol Biochem. 1993;25(9):1223–30. https://doi.org/10.1016/0038-0717(93)90218-Z.
Li Z, Lu Z. Soil Sucrase: detection conditions based on DNS colorimetric. Chin Agric Sci Bulletin. 2016;32:171–6.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.
Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinf. 2005;21(2):3787–93. https://doi.org/10.2307/1592215.
We should like to thank Professor Yiwei Mo of College of Life Science, Shaoxing University for reading this manuscript.
This work was supported by the Scientific Research Project of Public Welfare Industry (Agriculture) (201503232–05) and the Scientific Research Fund of Hunan Provincial Education Department (20A277). The supporters had no role in study design, data collection, data analysis, data interpretation, the writing of the manuscript or decision to publish.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rhizosphere soil phenolic acid detection with HPLC. FC stands for first cropping, and CC stands for continuous cropping. A: standard samples, B: FC soil, C: CC soil. 1: p-hydroxybenzoic, 2: vanillic acid, 3: syringic acid, 4: cumaric acid, 5: ferulic acid.
Data output quality summary.
Up- and down-regulated DEGs in CC vs FC root tissues.
Enriched pathways relative to phenolic acid synthesis.
Primers designed for qRT-PCR.
Data of up-regulated DEGs expression level and qRT-PCR verification.
Data of soil enzyme activities in FC and CC rhizosphere soil.
Data of phenolic acid contents in FC and CC rhizosphere soil. Typical curve of standard samples. Typical curve of CC soil. Typical curve of FC soil.
Data of heatmap of DEGs relative to phonelic acid metabolism in cc vs fc root.
Data of chemical properties in FC and CC rhizosphere soil.
About this article
Cite this article
Ni, X., Jin, C., Liu, A. et al. Physiological and transcriptomic analyses to reveal underlying phenolic acid action in consecutive monoculture problem of Polygonatum odoratum. BMC Plant Biol 21, 362 (2021). https://doi.org/10.1186/s12870-021-03135-x