Skip to main content

Cassava/peanut intercropping improves soil quality via rhizospheric microbes increased available nitrogen contents



Intercropping, an essential cultivation pattern in modern agricultural systems, increases crop yields and soil quality. Cassava and peanut intercropping systems exhibit advantages in solar utilization and cadmium absorption, etc. However, the inner mechanisms need to be elucidated. In this study, Illumina MiSeq platform was used to reveal the rhizospheric microbes and soil quality in cassava/peanut intercropping systems, and the results provided a reference for the application of this method in studying other intercropping systems.


Both intercropping cassava/peanut (IP) and intercropping peanut/cassava (IC) systems significantly increased available N, available K, pH value, and urease activity, comparing with that in monocropping cassava (MC) and monocropping peanut (MP) system. However, there were few effects on the total N, total P, total K, available P, organic matter, protease activity, catalase activity, sucrase activity, and acid phosphatase activity. Both IP and MP soils contained more bacteria and fungi than those in the IC and MC soils, which were mainly made of Proteobacteria and Actinobacteria. Intercropping remarkably increased the number of Nitrospirae in IP and IC soils comparing those in MC and MP soils. Redundancy analysis (RDA) revealed that the abundances of DA101, Pilimelia, and Ramlibacter were positively correlated to the soil quality. These results suggest that intercropping enhances the available nitrogen content of soil through increasing the quantity of rhizospheric microbes, especially that of DA101 and Pilimelia.


The cassava/peanut intercropping system improves soil quality through increasing the available nitrogen content and abundance of DA101, Pilimelia, and Ramlibacter in the soil.


Intercropping, an old and effective planting method of growing one crop/animals alongside another, can increase yields, reduce pests [1] and weeds [2], etc. There are many intercropping models, such as chickpea/maize [3], maize/lablab [4], rubber/plantain [5], even marine finfish in shrimp ponds [6]. Intercropping showed many advantages, like higher yield [7, 8], higher light interception and utilization rate [9], phytoremediation of heavy metal contaminated soils [10], enhancing iron nutrition [11] and phosphorus availability [12], etc.

Cassava (Manihot esculenta Crantz) is an important food crop in the world [13], especially in Africa [14]. Since its importance, breeding high yield and low cyanide cultivars [14], cultivation [15, 16], genome evolution [17] as well as diseases control were widely studied [18]. Peanut (Arachis hypogaea Linn.) is a legume crop that has edible seeds and oils. In cassava, cultivation and nitrogen fixation [19], abiotic stress-responsive ESTs [20], disease resistance genes [21] were well studied. Interestingly, continuous cropping of peanut changes the soil bacterial community [22]. Previous studies showed that peanut/maize intercropping changed rhizosphere and nutrient concentrations in shoots [23], but showed no effects on Aspergillus flavus in soil [24].

Recently, the inner mechanisms of intercropping are found related to microorganisms [25]. The progress in molecular [26] and microbiome techniques [27] provided new tools for the elucidation of the mechanisms of intercropping. For instance, molecular mechanisms of microbial disease control in intercropping were related to signals triggered by neighboring plants [28]. Sugarcane-soybean intercropping increased microbial diversity [29]. Microbiome-dependent immunity was related to soil organic matter content [30]. It was reported in previous studies that intercropping can change the soil microecology, as indicated by increasing farmland biodiversity [31]. Intercropping can effectively improve the mobilization and uptake of nitrogen (N), phosphorus (P), potassium (K), and micronutrients via interspecific interactions in the rhizosphere soil [23]. In addition, legume/cereal intercropping systems could improve the utilization of phosphorus (P) by root exudation of organic acids from legume crops which also improve legume N uptake by enhanced nodulation of legume crops [32]. Soil microbe and soil enzyme activities play essential roles in nutrient cycling, organic matter decomposition, and suppression of soil-borne pathogens in the rhizosphere [33, 34]. Plants can release root exudates, thereby affecting the rhizospheric microbial community [35]. Changes can influence the potential activities of soil enzymes in microbial community composition. Due to the quantitative and qualitative differences between the root exudates of intercropping and monocropping systems, differences in the microbial community can be observed [36]. Many studies have investigated the changes in the biochemical and microbial characteristics of rhizospheric soils caused by intercropping [37]. For the alfalfa/rye intercropping system, intercropping affected the soil microbial composition and soil enzymatic activities [38]. Through phospholipid fatty acid (PLFA) analysis, it was found that the soil urease activities, invertase activities, and the soil gram-negative (G-) bacterial abundance were significantly increased in the peanut/Atractylodes lancea system [39]. Many studies have shown that the maize/peanut intercropping system can facilitate the acquisition of Fe and Zn by peanut and improve the yields of both crops [23].

Cassava/peanut intercropping is a typical intercropping cultivation mode in southern China. Since the original spacing of the cassava crop remains unchanged when intercropped with peanut, thereby it provides a distinct yield advantage. Although researches were conducted regarding the selection of the cassava/peanut intercropping model and the associated yield benefits, study on the molecular mechanisms underlying this cultivation system remains insufficient. Previous studies on cassava/peanut intercropping mainly focused on the uptake and utilization of nutrients, photosynthesis, agronomic traits, yield, efficiency, and nutrient conversion efficiency in the soil [40]. However, the influence of the microecological soil environment in the cassava/peanut intercropping system remains unknown. In this study, the inner mechanisms of cassava/peanut intercropping were elucidated through the analysis of rhizospheric soil quality and microbial community using the Illumina MiSeq platform. We found that cassava/peanut intercropping enhances the quantity of DA101 and Pilimelia, then facilitate nitrogen use efficiency in plants. The results of this study provide a reference to applying this pattern in studying other intercropping systems.


Cassava/peanut intercropping enhanced soil physicochemical properties

As shown in Table 1, after three continuous years planting cassava and peanut between March and July of each year, both in the monocropping and intercropping system, the soil physicochemical properties were significantly changed. The available N, P, K increased nearly 20-folds, and the organic matter increased by almost 40% compared to the control soils (P < 0.01). IP (IP, i.e., planting in peanut former cassava field) and IC (IC, i.e., planting cassava in former peanut field) cultivation patterns significantly promoted the soil nutrient contents, especially in available N (P < 0.01) and pH value.

Table 1 Basic soil physicochemical properties in the rhizospheric soils of MP, MC, IP and IC cultivation patterns

To figure out which enzymes play a crucial role in the intercropping system, five major soil enzyme activities were measured (Table 2). The catalase, sucrase, protease, and acid phosphatase activities in the rhizospheric soil of IC and IP cultivation patterns showed no difference with those in the monocropping cassava (MC) and monocropping peanut (MP) plants. However, urease activity was highest in IC, which was 78.5% more than that in MC. These results suggested that intercropping cassava in the peanut field significantly enhances urease activity and available N contents.

Table 2 The activities of five major enzymes in the rhizospheric soils of MP, MC, IP and IC cultivation patterns

Cassava/peanut intercropping increased the quantity of culturable microbial in the rhizospheric soil

Since the rhizospheric soil physicochemical properties changed after IC and IP cultivation, we tested the microbial quantity to characterize the relationship between microbes and soil physicochemical properties (Table 3). IP cultivation pattern induced a significant increase in the bacterial abundance, fungal abundance, and total microbial amount comparing with that in MC. In contrast, microbial Shannon-Wiener diversity index in the rhizospheric soil of IP was less than that of MC. These results indicated that intercropping peanut in cassava fields increased microbial quantity but decreased microbial diversity. IC significantly reduced the bacterial abundance, fungal abundance, and total microbial amount comparing with that of MP. However, the microbial Shannon-Wiener diversity index in the rhizospheric soil of IC was increased. These results indicated that intercropping cassava in peanut fields decreased microbial amount but increased microbial diversity.

Table 3 Microbial quantity in the rhizospheric soils of MP, MC, IP and IC cultivation patterns

Cassava/peanut intercropping changed the microbial community

To identify the change of microbial community in the intercropping systems, a total of 43,519 valid sequences were yielded by 16S rRNA gene sequencing using the Illumina MiSeq platform, which represented the vast diversity of the bacterial community (Figs. 1, 2). The taxonomic distribution at the phylum level is shown in Fig. 1. Proteobacteria was the most abundant phylum, accounting for 28.24 to 37.45% of the total valid reads in all the samples, with an average relative abundance of 34.45%. Actinobacteria was the second most abundant phylum, with an average relative abundance of 20.70%. The other dominant phyla were Acidobacteria (12.55–18.32%, with an average value of 14.79%), Chloroflexi (7.19–8.40%, with an average value of 7.76%), Gemmatimonadetes (3.80–5.83%, with an average value of 4.63%), Nitrospirae (3.27–5.81%, with an average value of 4.16%), Planctomycetes (1.60–4.36%, with an average value of 2.96%), Verrucomicrobia (1.46–4.51%, with an average value of 2.71%), and Bacteroidetes (1.03–3.75%, with an average value of 1.91%). Importantly, the percentages of Nitrospirae, Verrucomicrobia and Gemmatimonadetes in the rhizospheric soils of the IP and IC intercropping systems were more than those in the monocropping systems. Bacteroidetes and Planctomycetes were also more abundant in the rhizospheric soil of the intercropping system IP than those in the monocropping system MP. These phyla were also less abundant in MC than that in IC. Other phyla, such as Proteobacteria, Actinobacteria, Acidobacteria and Chloroflexi did not exhibit a significant difference between the monoculture and intercropping systems. On the genera levels (Fig. 2), four kinds of soils showed different dominant genera, for instance, Aquicella, Chthonomonas, Kribbella, DA101 and Nitrospira were more abundant in the rhizospheric soil of the cassava plants in the IP and IC. Actinoallomurus and Streptomyces were still more abundant in the MC. Comparing with the rhizospheric soil of cassava, the rhizospheric soil of peanut exhibited higher diversity of dominant genera. There were 15 dominant genera in the rhizospheric soil of the intercropped peanut plants, including Optitutus, A4, Chthoniobacter, Flavisolibater, Dokdonella, and Pilimelia. However, these genera were not highly abundant in the peanut plants of monoculture system, which had different dominant genera, such as Chryseobacterium, Alicyclobacillus, Escherichia, Ralstonia, and Hypomicrobium. Thus, these results clearly demonstrated that the cassava/peanut intercropping changed the microbial community, which became distinct from that of the monocropping system. Furthermore, the microbial community in the rhizospheric soils of peanut and cassava plants is different due to different planting patterns, i.e. monocropping and intercropping systems.

Fig. 1

Taxonomic classification of bacterial reads in the rhizospheric soils of different cultivation patterns. MC, monocropping cassava; MP, monocropping peanut; IC, intercropping peanut/cassava; IP, intercropping cassava/peanut

Fig. 2

Heatmap analysis of the dominant genera in the rhizospheric soils of different cultivation patterns. MC, monocropping cassava; MP, monocropping peanut; IC, intercropping peanut/cassava; IP, intercropping cassava/peanut

Principal component analysis (PCA) and redundancy analysis (RDA) of microbial communities and physicochemical properties in different cultivation patterns

Principal component analysis showed that samples from monocropping and intercropping systems were separated from each other. The MP and IP were distributed in quadrant 2 and quadrant 1, 4, respectively (Fig. 3). It suggested that the same crop in different planting patterns influences the microbial communities accordingly. This finding may be associated with the environmental parameters in the different planting patterns. Microbial community exhibited a high correlation with intrinsic ecological parameters. Relationships between the important ecological parameters and the microbial community were discerned by RDA (Fig. 4). The length of the arrow corresponding to an ecological parameter indicated the strength of the ecological parameter concerning the overall microbial community. The results of RDA suggested that there were significant differences in the bacterial communities in the four planting patterns. As shown in Fig. 4, the available N, catalase, organic matter, sucrase activity, acid phosphatase activity, urease activity, total N, pH, total K, and available K were positively correlated with the RDA axis 1, and they were strongly and significantly associated with the overall microbial communities in IC and IP. In contrast, the total P, available P, and protease activity were negatively correlated with the RDA axis 1. These results revealed that available N, pH, catalase activity, and sucrase activity had the most significant impact on the microbial communities. Additionally, the abundances of some microbial genera, such as DA101, Pilimelia, and Ramlibacter, were positively correlated with available N, which were also the dominant genera in the soil of intercropped peanut plants (Fig. 1).

Fig. 3

Principal component analyses of bacterial community in the rhizospheric soils of monocropping and intercropping systems based on Euclidean distance. MC, monocropping cassava; MP, monocropping peanut; IC, intercropping peanut/cassava; IP, intercropping cassava/peanut

Fig. 4

Redundancy analysis (RDA) of sequencing data of 16S rRNA gene and physicochemical properties in the rhizospheric soils of different cultivation patterns. MC, monocropping cassava; MP, monocropping peanut; IC, Intercropping peanut/cassava; IP, intercropping cassava/peanut


Cassava/peanut intercropping system changed physicochemical properties of rhizospheric soils

Intercropping systems show great importance in agronomy, plant physiology, and ecology [41]. In general, intercropping systems change the bacterial diversity of soils [42], decrease disease rates [43, 44], and increase yields and cadmium accumulation [45]. In this study, rhizospheric soils both in IP and IC systems showed higher available N and pH than those in MP and MC systems, respectively. These were consistent with previous studies. For instance, intercropping of green garlic (Allium sativum L.) with cucumber (Cucumis sativus L.) increased organic matter and available N, P and K contents in soils [46]. Eggplant/garlic relay intercropping system increased available N contents from 61.95 mg kg− 1 to 76.30 mg kg− 1 [47]. In sugarcane-soybean intercropping system, available N increased by 10% [29]. The increase of available N was related to the urease activities, which was the only significant change in IC and IP comparing with that in MC and MP (Table 2). The role of urease is to catalyze urea hydrolysis into ammonia and carbon dioxide [48, 49], which is common in higher plants, bacteria, fungi, and algae. Thus, we speculate that the changes of soil physicochemical properties may be related to urease produced by microbial community since both the dominant genera and quantity of the microbial community were changed in IC and IP comparing with that in MC and MP (Table 3).

Improvement of physicochemical properties of rhizospheric soils was related to the microbial community in cassava/peanut intercropping system

Through MiSeq sequencing analysis of the 16S rRNA gene, we found that the percentage of Nitrospirae, Verrucomicrobia, and Gemmatimonadetes in the rhizospheric soils of IP and IC intercropping systems were higher than those in the monocropping systems (Fig. 1). This finding is quite different from the result of a previous study using a different detection technique [50]. In the previous study, denaturing gradient gel electrophoresis (DGGE) was used to analyze microbial community structure and diversity, and only ten kinds of fungi and bacteria were identified [50]. In this study, the new MiSeq sequencing technique was used, and a total of 50 species were identified (Fig. 3). Proteobacteria and Actinobacteria were the most abundances bacteria in different intercropping systems [49], while other species were varying. For instance, Firmicutes was the dominant genus in mulberry/alfalfa intercropping soils [49], but it accounted for only 0.54–1.33% in this study (Fig. 1). These results suggested that each intercropping system had different mechanisms, especially in the composition of bacterial community. Based on genus analysis, PCA and RDA analysis, it was clear that the differences in the IC and IP systems comparing with those in MC and MP systems (Figs. 2, 3 and 4) were related to the abundance of DA101, Pilimeli, and Ramlibac. Then these bacteria increased enzyme activities and improved the physicochemical properties of rhizospheric soil (Fig. 5). Unclassified DA101 belongs to the family Chthoniobacteraceae, phylum Verrucomicrobia, which was found in the soils of pepper tree [51], and grassland [52]. The abundance of DA101 was increased in the soil of intercropping system, but it was negatively correlated with the total N content and pH at significant level [53]. In this study, we found that DA101 abundance was significantly and positively correlated with the available N content in cassava/peanut intercropping system, which was consistent with the change in the abundance of the phylum Verrucomicrobia in the intercropping system (Fig. 2-4). These results indicated that the available N content, and not the total N content was related to DA101 abundance. Recently, researchers assembled a draft genome of Candidatus Udaeobacter copiosus, which is a representative of the DA101 clade, and speculated that this organism is a soil oligotrophic bacterium which reduces its requirement for soil organic carbon by acquiring amino acids and vitamins from the environment [54]. The functions of DA101 need to be investigated more deeply. Ramlibacter tataouinensis TTB310T (strain TTB310) is a beta-proteobacterium isolated from sand particles [55]. Different species from this genus were isolated in tropical forest soil [56] and the rhizospheric soil of Mugunghwa [57]. Furthermore, studies revealed that the genus of Ramlibacter has phosphatase activity [58] and ginsenoside-converting activity [59]. The relationship between these bacteria and enzymes activities should be studied further to reveal their functions.

Fig. 5

Models of IP and IC enhanced soil physicochemical properties. MC, monocropping cassava; MP, monocropping peanut; IC, peanut/cassava intercropping; IP, cassava/peanut intercropping


The inner mechanism of cassava/peanut intercropping system was elucidated through analyzing the physicochemical properties of rhizospheric soils and microbial community. We found that the cassava/peanut intercropping system increased the quantities of DA101, Pilimelia, and Pilimelia, and thereby increased the available N content in the soil and improved soil quality.


Cassava and peanut plants

Cassava variety Huanan 205 and shade-tolerant peanut variety Guihua 836 were provided by the Cash Crops Research Institute of the Guangxi Academy of Agricultural Sciences. All cassava and peanut experiments were performed in the Lijian Scientific Base, which is certified for the field cultivation experiment by local government.

Experimental site and soil

The experiments were performed in the Lijian Scientific Base (23°14′25′′N, 108°03′42′′E) of the Guangxi Academy of Agricultural Sciences, Nanning City, Guangxi province, China. The field site was previously used for monocropping cassava and peanut. The tested soil was acid red loam, with total nitrogen content (total N), total phosphorus content (total P), total potassium content (total K), available nitrogen content (available N), available phosphorus content (available P), available potassium content (available K), the organic matter content and pH value of 1.34 g kg− 1, 0.53 g kg − 1, 12.6 g kg − 1,0.0705 g kg − 1, 0.0139 g kg − 1, 0.098 g kg − 1, 16.2 g kg − 1, and 5.8, respectively.

Experimental design and management

In March 2016–2018, cassava and peanut were planted simultaneously in the field. Monocropping cassava (MC), monocropping peanut (MP) crops were compared with peanut/cassava intercropping (IC, i.e., planting cassava in former peanut field) and cassava/peanut intercropping (IP, i.e., planting peanut in former cassava field) (Fig. 6). For MC, cassava planted with a row spacing of 1.1 m × 0.8 m and with equivalent line spacing. For MP, the peanut planted in a narrow-wide row spacing pattern. The line spacing for peanut in the wide line was 0.5 m, and the row spacing in the narrow row was 0.3 m × 0.16 m. For IC and IP, two lines of peanut were planted alongside one line of cassava. The line spacing between cassava and peanut was 0.4 m. The row spacing for cassava and peanut intercropping was 1.1 m × 0.8 m and 0.3 m × 0.16 m, respectively. The experiment was arranged in plots (6 m × 8 m) in a randomized design with three replicates in each treatment.

Fig. 6

Diagram of experimental design

Peanut was supplied with 450 kg ha− 1 compound granulated NPK fertilizers (N-P2O5-K2O = 15–15-15) and 750 kg. ha− 1 fused calcium-magnesium phosphate fertilizer (available P2O5 18%). Cassava was supplied with 750 kg Ha− 1 compound granulated NPK fertilizers. The crops were irrigated two times during crop growth period based on crop water requirement and soil water content. Pesticide and herbicide were applied about two months after sowing.

Soil sampling

In July 2016–2018, the time of harvesting mature peanut, ten plants of cassava and peanut per treatments were uprooted. The rhizospheric soil from both loose soil and cohesive soil from the plant roots were collected, mixed, and separated into three sealed virus-free bags for the following assays. One sample was maintained in the refrigerator at 4 °C and used for determining the culturable soil microbe. One sample was maintained in the refrigerator at − 80 °C and used for extracting soil DNA and high-throughput sequencing. The last sample was dried naturally, grounded, and sieved for determining the soil physicochemical properties.

Soil physicochemical property analysis

The physicochemical properties were measured according to the previous reports [49]. The available N contents were measured by the alkaline hydrolysis diffusion method. The available P was measured by sodium bicarbonate extraction/Mo-Sb colorimetry method. The available K was measured by ammonium acetate extraction/flame photometry method. Organic matter content was measured by the potassium bichromate titrimetric method. Soil enzyme activities were measured according to the previous reports. Catalase activity [60], sucrase activity, proteinase activity, urease activity [48], and acid phosphatase activity [3] were measured by permanganate titration, sodium thiosulfate titration, ninhydrin colorimetry, indophenol blue colorimetry, and the disodium phosphate benzene colorimetric method, respectively.

Determination of soil microbial quantity

Soil microbial quantity was measured by the conventional microculture method. Bacteria, fungi, and actinomycetes were cultured in beef extract-peptone medium, Martin medium, and Gao 1 medium, respectively. The Shannon-Wiener index was used to calculate the biodiversity index (H): H = −Σ(ni/N) × ln(ni/N). In this formula, ni is the microbial quantity of species i, and N is the total microbial quantity [61].

Soil DNA extraction, PCR amplification, high-throughput sequencing and analysis

Total genomic DNA was extracted from the samples using the FastDNA SPIN Kit (MP Biomedicals, Santa Ana, USA) according to the manufacturer’s instructions. The DNA concentration and purity were monitored on 1% agarose gels, and the DNA was diluted to 1 ng/μL using sterile water. The DNA was stored at − 80 °C until subsequent PCR amplification. The total genomic DNA was subjected to PCR amplification using the primer pair 515f/806r, which amplified the V4 region of the 16S rDNA gene [62], following a previously described protocol [63]. PCR was carried out in 30-μL reactions with 15 μL of Phusion® high-fidelity PCR master mix, 0.2 μM of forward and reverse primers, and approximately 10 ng of template DNA. Thermal cycling consisted of an initial denaturation at 98 °C for 1 min; followed by 30 cycles of denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, and elongation at 72 °C for 30 s; and a final step at 72 °C for 5 min. The PCR products were mixed with equal volumes of 1× loading buffer (containing SYBR Green), and electrophoresis conducted on a 2% agarose gel for detection. Samples with a bright band at 400–450 bp were chosen for further experiments. The PCR products were mixed at equal concentrations. Then, the mixture of PCR products were purified with the GeneJET Gel Extraction Kit (Thermo Scientific). Sequencing libraries were generated using the NEB Next® Ultra™ DNA Library Prep Kit for Illumina (NEB, USA) following the manufacturer’s recommendations, and the index codes were added. The library quality was assessed on a Qubit 2.0 fluorometer (Thermo Scientific) and an Agilent 2100 bioanalyzer system. Finally, the library was sequenced on an Illumina MiSeq platform by the Novogene Corporation (Beijing, China), and 250/300-bp paired-end reads were generated.

Illumina MiSeq sequencing analysis

Paired-reads from the original DNA fragments were merged based on a previously described method [64]. Sequencing reads were assigned to each sample according to the individual unique barcodes, then analyzed with the QIIME (Quantitative Insights Into Microbial Ecology) software package and the UPARSE pipeline [65]. The reads were first filtered by QIIME quality filters. Default settings for Illumina processing in QIIME were used. Then, the UPARSE pipeline was used to select operational taxonomic units (OTUs) with 97% similarity. For each OTU, a representative sequence was selected and used to assign taxonomic composition by the RDP classifier. Then, the estimated species richness was determined by rarefaction analysis [66]. Redundancy analysis (RDA) was performed to analyze the correlation between environmental factors and microbial community.

Data analysis

Principal component analysis (PCA) was performed with SIMCA software v.13.0 (Umetrics, Sweden) [67]. The means and standard errors of three repeats in MC, MP, IC, and IP planting patterns were analyzed by One-way variance analysis with SPSS24.0, and Duncan’s test of the homogeneity of variance was performed with the confidence level of 0.01 [68].

Availability of data and materials

Raw data of 16S rRNA gene obtained from all samples are accessible via NCBI SRA database under accession number PRJNA606845.



Intercropping peanut/cassava


Intercropping cassava/peanut


Monocropping cassava


Monocropping peanut


Operational taxonomic units


Principal component analysis


Redundancy analysis


  1. 1.

    Singh A, Weisser WW, Hanna R, Houmgny R, Zytynska SE. Reduce pests, enhance production: benefits of intercropping at high densities for okra farmers in Cameroon. Pest Manag Sci. 2017;73:2017–27.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Liebman M, Dyck E. Crop rotation and intercropping strategies for weed management. Ecol Appl. 1993;3:92–122.

    PubMed  Article  Google Scholar 

  3. 3.

    Li SM, Li L, Zhang FS, Tang C. Acid phosphatase role in chickpea/maize intercropping. Ann Bot. 2004;94:297–303.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Maluleke MH, Addo-Bediako A, Ayisi KK. Influence of maize/lablab intercropping on lepidopterous stem borer infestation in maize. J Econom Entomol. 2005;98:384–8.

    Article  Google Scholar 

  5. 5.

    Tetteh EN, Abunyewa AA, Tuffour HO, Berchie JN, Acheampong PP, Twum-Ampofo K, Dawoe E, Logah V, Agbenyega O, Ennin SA, Nunoo I, Melenya C, Danquah EO, Barnes VR, Partey ST. Rubber and plantain intercropping: effects of different planting densities on soil characteristics. PLoS One. 2019;14:e0209260.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Damodaran D, Mojjada SK, Vase VK, Sukhdhane K, Abdul Azeez P, Kumar R. Intercropping of marine finfish in shrimp ponds: A maiden feasibility study. PLoS One. 2019;14:e0216648.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Chen P, Song C, Liu XM, Zhou L, Yang H, Zhang X, Zhou Y, Du Q, Pang T, Fu ZD, Wang XC, Liu WG, Yang F, Shu K, Du JB, Liu J, Yang WY, Yong TW. Yield advantage and nitrogen fate in an additive maize-soybean relay intercropping system. Sci Total Environ. 2019;657:987–99.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Dong N, Tang MM, Zhang WP, Bao XG, Wang Y, Christie P, Li L. Temporal differentiation of crop growth as one of the drivers of intercropping yield advantage. Sci Rep. 2018;8:3110.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  9. 9.

    Zhi XY, Han YC, Xing FF, Lei YP, Wang GP, Feng L, Yang BF, Wang ZB, Li XL, Xiong SW, Fan ZY, Li YB. How do cotton light interception and carbohydrate partitioning respond to cropping systems including monoculture, intercropping with wheat, and direct-seeding after wheat? PLoS One. 2019;14:e0217243.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    De Conti L, Ceretta CA, Melo GWB, Tiecher TL, Silva LOS, Garlet LP, Mimmo T, Cesco S, Brunetto G. Intercropping of young grapevines with native grasses for phytoremediation of cu-contaminated soils. Chemosphere. 2019;216:147–56.

    PubMed  Article  CAS  Google Scholar 

  11. 11.

    Dai J, Qiu W, Wang N, Nakanishi H, Zuo Y. Comparative transcriptomic analysis of the roots of intercropped peanut and maize reveals novel insights into peanut iron nutrition. Plant Physiol Biochem. 2018;127:516–24.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Darch T, Giles CD, Blackwell MSA, George TS, Brown LK, Menezes-Blackburn D, Shand CA, Stutter MI, Lumsdon DG, Mezeli MM, Wendler R, Zhang H, Wearing C, Cooper P, Haygarth PM. Inter- and intra-species intercropping of barley cultivars and legume species, as affected by soil phosphorus availability. Plant Soil. 2018;427:125–38.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Okezie BO, Kosikowski FV. Cassava as a food. Crit Rev Food Sci Nutr. 1982;17:259–75.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Onabolu A, Bokanga M, Tylleskar T, Rosling H. High cassava production and low dietary cyanide exposure in mid-West Nigeria. Public Health Nutr. 2001;4:3–9.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Schulthess F, Chabi-Olaye A, Gounou S. Multi-trophic level interactions in a cassava-maize mixed cropping system in the humid tropics of West Africa. Bull Entomol Res. 2004;94:261–72.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Aina OO, Dixon AG, Akinrinde EA. Effect of soil moisture stress on growth and yield of cassava in Nigeria. Pak J Biol Sci. 2007;10:3085–90.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Sardos J, McKey D, Duval MF, Malapa R, Noyer JL, Lebot V. Evolution of cassava (Manihot esculenta Crantz) after recent introduction into a South Pacific Island system: the contribution of sex to the diversification of a clonally propagated crop. Genome. 2008;51:912–21.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Patil BL, Fauquet CM. Cassava mosaic geminiviruses: actual knowledge and perspectives. Mol Plant Pathol. 2009;10:685–701.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    El-Akhal MR. Rincon a, Coba de la Pena T, Lucas MM, El Mourabit N, Barrijal S, Pueyo JJ. Effects of salt stress and rhizobial inoculation on growth and nitrogen fixation of three peanut cultivars. Plant Biol (Stuttg). 2013;15:415–21.

    CAS  Article  Google Scholar 

  20. 20.

    Kumari A, Kumar A, Wany A, Prajapati GK, Pandey DM. Identification and annotation of abiotic stress responsive candidate genes in peanut ESTs. Bioinformation. 2012;8:1211–9.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Kumar D, Kirti PB. Pathogen-induced SGT1 of Arachis diogoi induces cell death and enhanced disease resistance in tobacco and peanut. Plant Biotechnol J. 2015;13:73–84.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Chen M, Li X, Yang Q, Chi X, Pan L, Chen N, Yang Z, Wang T, Wang M, Yu S. Dynamic succession of soil bacterial community during continuous cropping of peanut (Arachis hypogaea L.). PLoS One. 2014;9:e101355.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  23. 23.

    Inal A, Gunes A, Zhang F, Cakmak I. Peanut/maize intercropping induced changes in rhizosphere and nutrient concentrations in shoots. Plant Physiol Biochem. 2007;45:350–6.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Horn BW, Greene RL, Dorner JW. Effect of corn and peanut cultivation on soil populations of Aspergillus flavus and A. parasiticus in southwestern Georgia. Appl Environ Microbiol. 1995;61:2472–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Chen M, Li X, Yang Q, Chi X, Pan L, Chen N, Yang Z, Wang T, Wang M, Yu S. Soil eukaryotic microorganism succession as affected by continuous cropping of peanut--pathogenic and beneficial fungi were selected. PLoS One. 2012;7:e40659.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Nelkner J, Henke C, Lin TW, Patzold W, Hassa J, Jaenicke S, Grosch R, Puhler A, Sczyrba A, Schluter A. Effect of long-term farming practices on agricultural soil microbiome members represented by metagenomically assembled genomes (MAGs) and their predicted plant-beneficial genes. Genes (Basel). 2019:10, 424.

  27. 27.

    Wang G, Schultz P, Tipton A, Zhang J, Zhang F, Bever JD. Soil microbiome mediates positive plant diversity-productivity relationships in late successional grassland species. Ecol Lett. 2019;22:1221–32.

    PubMed  Google Scholar 

  28. 28.

    Zhu S, Morel JB. Molecular mechanisms underlying microbial disease control in intercropping. Mol Plant-Microbe Interactions. 2019;32:20–4.

    CAS  Article  Google Scholar 

  29. 29.

    Lian T, Mu Y, Jin J, Ma Q, Cheng Y, Cai Z, Nian H. Impact of intercropping on the coupling between soil microbial community structure, activity, and nutrient-use efficiencies. Peer J. 2019;7:e6412.

    PubMed  Article  CAS  Google Scholar 

  30. 30.

    Haron MH, Tyler HL, Chandra S, Moraes RM, Jackson CR, Pugh ND, Pasco DS. Plant microbiome-dependent immune enhancing action of Echinacea purpurea is enhanced by soil organic matter content. Sci Rep. 2019;9:136.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  31. 31.

    Zhang F, Shen J, Li L, Liu X. An overview of rhizosphere processes related with plant nutrition in major cropping systems in China. Plant Soil. 2004;260:89–99.

    CAS  Article  Google Scholar 

  32. 32.

    Fang Q, Yu Q, Wang E, Chen Y, Zhang G, Wang J, Li L. Soil nitrate accumulation, leaching and crop nitrogen use as influenced by fertilization and irrigation in an intensive wheat–maize double cropping system in the North China plain. Plant Soil. 2006;284:335–50.

    CAS  Article  Google Scholar 

  33. 33.

    Burns RG, DeForest JL, Marxsen J, Sinsabaugh RL, Stromberger ME, Wallenstein MD, Weintraub MN, Zoppini A. Soil enzymes in a changing environment: current knowledge and future directions. Soil Biol Biochem. 2013;58:216–34.

    CAS  Article  Google Scholar 

  34. 34.

    Pignataro A, Moscatelli MC, Mocali S, Grego S, Benedetti A. Assessment of soil microbial functional diversity in a coppiced forest system. App Soil Ecol. 2012;62:115–23.

    Article  Google Scholar 

  35. 35.

    Kourtev P, Ehrenfeld J, Häggblom M. Experimental analysis of the effect of exotic and native plant species on the structure and function of soil microbial communities. Soil Biol Biochem. 2003;35:895–905.

    CAS  Article  Google Scholar 

  36. 36.

    Baudoin E, Benizri E, Guckert A. Impact of artificial root exudates on the bacterial community structure in bulk soil and maize rhizosphere. J Soil Biol Biochem. 2003;35:1183–92.

    CAS  Article  Google Scholar 

  37. 37.

    Bainard LD, Koch AM, Gordon AM, Klironomos JN. Growth response of crops to soil microbial communities from conventional monocropping and tree-based intercropping systems. Plant Soil. 2012;363:345–56.

    Article  CAS  Google Scholar 

  38. 38.

    Sun YM, Zhang NN, Wang ET, Yuan HL, Yang JS, Chen WX. Influence of intercropping and intercropping plus rhizobial inoculation on microbial activity and community composition in rhizosphere of alfalfa (Medicago sativa L.) and Siberian wild rye (Elymus sibiricus L.). FEMS Microbiol Ecol. 2009;70:62–70.

    PubMed  Article  CAS  Google Scholar 

  39. 39.

    Dai CC, Chen Y, Wang XX, Li PD. Effects of intercropping of peanut with the medicinal plant Atractylodes lancea on soil microecology and peanut yield in subtropical China. Agrofor Syst. 2012;87:417–26.

    Article  Google Scholar 

  40. 40.

    Polthanee A, Wanapat S, Mangprom P. Row arrangement of peanut in cassava-peanut intercropping: II. Kaen Kaset: Nutrient removal and nutrient balance in soil; 1998.

    Google Scholar 

  41. 41.

    Brooker RW, Bennett AE, Cong WF, Daniell TJ, George TS, Hallett PD, Hawes C, Iannetta PP, Jones HG, Karley AJ, Li L, McKenzie BM, Pakeman RJ, Paterson E, Schob C, Shen J, Squire G, Watson CA, Zhang C, Zhang F, Zhang J, White PJ. Improving intercropping: a synthesis of research in agronomy, plant physiology and ecology. New Phytol. 2015;206:107–17.

    PubMed  Article  Google Scholar 

  42. 42.

    Cao X, Liu S, Wang J, Wang H, Chen L, Tian X, Zhang L, Chang J, Wang L, Mu Z, Qiao Z. Soil bacterial diversity changes in different broomcorn millet intercropping systems. J Basic Microbiol. 2017;57:989–97.

    PubMed  Article  Google Scholar 

  43. 43.

    Boudreau MA. Diseases in intercropping systems. Annu Rev Phytopathol. 2013;51:499–519.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Damicone JP, Edelson JV, Sherwood JL, Myers LD, Motes JE. Effects of border crops and intercrops on control of cucurbit virus diseases. Plant Dis. 2007;91:509–16.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Deng Q, Deng Q, Wang Y, Li L, Long X, Ren S, Fan Y, Lin L, Xia H, Liang D, Wang J, Zhang H, Lv X, Wang Y. Effects of intercropping with Bidens species plants on the growth and cadmium accumulation of Ziziphus acidojujuba seedlings. Environ Monit Assess. 2019;191:342.

    PubMed  Article  CAS  Google Scholar 

  46. 46.

    Xiao X, Cheng Z, Meng H, Liu L, Li H, Dong Y. Intercropping of green garlic (Allium sativum L.) induces nutrient concentration changes in the soil and plants in continuously cropped cucumber (Cucumis sativus L.) in a plastic tunnel. PLoS One. 2013;8:e62173.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Wang M, Wu C, Cheng Z, Meng H, Zhang M, Zhang H. Soil chemical property changes in eggplant/garlic relay intercropping systems under continuous cropping. PLoS One. 2014;9:e111040.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Cordero I, Snell H, Bardgett RD. High throughput method for measuring urease activity in soil. Soil Biol Biochem. 2019;134:72–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Zhang MM, Wang N, Hu YB, Sun GY. Changes in soil physicochemical properties and soil bacterial community in mulberry (Morus alba L.)/alfalfa (Medicago sativa L.) intercropping system. Microbiologyopen. 2018;7:e00555.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 50.

    Xu HQ, Huang J, Liu ZF, Wei YX, Su BM, Li T. Effects of cassava-peanut intercropping on microbial amount, community structure and diversity in rhizosphere soils. [in Chinese]. J Southern Agricul. 2016;47:185–90.

    Google Scholar 

  51. 51.

    Dawkins K, Esiobu N. The invasive brazilian pepper tree (Schinus terebinthifolius) is colonized by a root microbiome enriched with Alphaproteobacteria and unclassified Spartobacteria. Front Microbiol. 2018;9:876.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Felske A, Akkermans AD. Prominent occurrence of ribosomes from an uncultured bacterium of the Verrucomicrobiales cluster in grassland soils. Lett App Microbiol. 1998;26:219–23.

    CAS  Article  Google Scholar 

  53. 53.

    Shen C, Ge Y, Yang T, Chu H. Verrucomicrobial elevational distribution was strongly influenced by soil pH and carbon/nitrogen ratio. J Soils Sediments. 2017;17:2449–56.

    CAS  Article  Google Scholar 

  54. 54.

    Brewer TE, Handley KM, Carini P, Gilbert JA, Fierer N. Genome reduction in an abundant and ubiquitous soil bacterium ‘Candidatus Udaeobacter copiosus’. Nat Microbiol. 2016;2:1–7.

    Google Scholar 

  55. 55.

    De Luca G, Barakat M, Ortet P, Fochesato S, Jourlin-Castelli C, Ansaldi M, Py B, Fichant G, Coutinho PM, Voulhoux R, Bastien O, Marechal E, Henrissat B, Quentin Y, Noirot P, Filloux A, Mejean V, DuBow MS, Barras F, Barbe V, Weissenbach J, Mihalcescu I, Vermeglio A, Achouak W, Heulin T. The cyst-dividing bacterium Ramlibacter tataouinensis TTB310 genome reveals a well-stocked toolbox for adaptation to a desert environment. PLoS One. 2011;6:e23784.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  56. 56.

    Zhang XJ, Feng GD, Yao Q, Wang YH, Yang SZ, Zhu HH. Ramlibacter humi sp. nov., isolated from tropical forest soil. Int J Syst Evol Microbiol. 2019;69:3460–4.

    PubMed  Article  Google Scholar 

  57. 57.

    Yan ZF, Trinh H, Moya G, Lin P, Li CT, Kook M, Yi TH. Ramlibacter rhizophilus sp. nov., isolated from rhizosphere soil of national flower Mugunghwa from South Korea. Int J Syst Evol Microbiol. 2017;67:3773–7.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Skouri-Panet F, Benzerara K, Cosmidis J, Ferard C, Caumes G, De Luca G, Heulin T, Duprat E. In vitro and in silico evidence of phosphatase diversity in the biomineralizing bacterium Ramlibacter tataouinensis. Front Microbiol. 2017;8:2592.

    PubMed  Article  Google Scholar 

  59. 59.

    Wang L, An DS, Kim SG, Jin FX, Kim SC, Lee ST, Im WT. Ramlibacter ginsenosidimutans sp. nov., with ginsenoside-converting activity. J Microbiol Biotechnol. 2012;22:311–5.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Mueller S, Riedel H-D, Stremmel W. Determination of catalase activity at physiological hydrogen peroxide concentrations. Anal Biochem. 1997;245:55–60.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Weaver MA, Krutz LJ, Zablotowicz RM, Reddy KN. Effects of glyphosate on soil microbial communities and its mineralization in a Mississippi soil. Pest Manag Sci. 2007;63:388–93.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Bergmann GT, Bates ST, Eilers KG, Lauber CL, Caporaso JG, Walters WA, Knight R, Fierer N. The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol Biochem. 2011;43:1450–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, Gormley N, Gilbert JA, Smith G, Knight R. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Magoc T, Salzberg SL. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Huang L, Fuchino H, Kawahara N, Narukawa Y, Hada N, Kiuchi F. Application of a new method, orthogonal projection to latent structure (OPLS) combined with principal component analysis (PCA), to screening of prostaglandin E2 production inhibitory flavonoids in Scutellaria root. J Nat Med. 2016;70:731–9.

    CAS  PubMed  Article  Google Scholar 

  68. 68.

    Guo M, Rupe MA, Wei J, Winkler C, Goncalves-Butruille M, Weers BP, Cerwick SF, Dieter JA, Duncan KE, Howard RJ, Hou Z, Loffler CM, Cooper M, Simmons CR. Maize ARGOS1 (ZAR1) transgenic alleles increase hybrid maize yield. J Exp Bot. 2013;65:249–60.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references


We thank for Wilking Biotechnologe Co., Ltd., (Nanning, Guangxi, China) for data collection and analysis, and thanks to Dr. Ning Yu for helps in the graph drawing.


These works were supported by the National Natural Science Foundation Project (31401318 & 31660371), Technology System for Modern Agriculture (CARS-13-Southern China), Guangxi Natural Science Foundation Project (2017GXNSFAA198144) and Scientific Project from the Guangxi Academy of Agricultural Sciences (2015JZ12 &2017JZ12), which had no role in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.

Author information




TX conceived and designed the experiments and wrote the manuscript. TR and HL revised the manuscript and provided critical advice. TX, ZR, JJ, HL, HZ, WH, XF, and HZ performed the experiments. SG and LJ analyzed the data. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Xiumei Tang or Ronghua Tang or Longfei He.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tang, X., Zhong, R., Jiang, J. et al. Cassava/peanut intercropping improves soil quality via rhizospheric microbes increased available nitrogen contents. BMC Biotechnol 20, 13 (2020).

Download citation


  • Cassava
  • Intercropping
  • Microbial community
  • Peanut
  • Rhizospheric soil