- 1College of Agronomy, State Key Laboratory of Crop Stress Biology in Arid Areas/Northwest A & F University, Yangling, China
- 2Belarusian State Agricultural Academy, Gorki, Belarus
Intercropping of cereals and legumes has been used in modern agricultural systems, and the soil microorganisms associated with legumes play a vital role in organic matter decomposition and nitrogen (N) fixation. This study investigated the effect of intercropping on the rhizosphere soil microbial composition and structure and how this interaction affects N absorption and utilization by plants to improve crop productivity. Experiments were conducted to analyze the rhizosphere soil microbial diversity and the relationship between microbial composition and N assimilation by proso millet (Panicum miliaceum L.) and mung bean (Vigna radiata L.) from 2017 to 2019. Four different intercropping row arrangements were evaluated, and individual plantings of proso millet and mung bean were used as controls. Microbial diversity and community composition were determined through Illumina sequencing of 16S rRNA and internal transcribed spacer (ITS) genes. The results indicated that intercropping increased N levels in the soil–plant system and this alteration was strongly dependent on changes in the microbial (bacterial and fungal) diversities and communities. The increase in bacterial alpha diversity and changes in unique operational taxonomic unit (OTU) numbers increased the soil N availability and plant N accumulation. Certain bacterial taxa (such as Proteobacteria) and fungal taxa (such as Ascomycota) were significantly altered under intercropping and showed positive responses to increased N assimilation. The average grain yield of intercropped proso millet increased by 13.9–50.1% compared to that of monoculture proso millet. Our data clearly showed that intercropping proso millet with mung bean altered the rhizosphere soil microbial diversity and community composition; thus, this intercropping system represents a potential mechanism for promoting N assimilation and increasing grain yield.
Introduction
With the increasing global demand for food, the relationship between crop production and food security should be determined and natural resources must be preserved (Banik and Sharma, 2009). Although industrial agriculture is directly beneficial to improving labor efficiency and crop production, intense fertilizer use has led to a series of ecological environmental problems, such as loss of diversity in ecosystems, decreases in soil fertility and aggravation of environmental pollution (Boardman et al., 2003; Jacobsen et al., 2013). Intercropping is a useful agricultural practice that permits the simultaneous growth of two or more crops in the same field, thereby improving the land use efficiency (Yu et al., 2017). This agricultural practice is a technological method based on the ecological principles of facilitation and complementarity (Duchene et al., 2017). Thus, intercropping has not only been adopted by developing countries but also European countries (Martin-Guay et al., 2017). Among different intercropping combinations, cereal–legume intercropping systems have become sustainable farming models because these types of crops are not competing for the same niche (Li et al., 2001) and have different nitrogen (N) use abilities or obtain N via different pathways (e.g., by mineral or organic fertilizer in cereals and N fixation in legumes) (Ghosh et al., 2009). Legumes reduce N input requirements by biological N fixation, which meets 50–60% of the N demand (Salvagiotti et al., 2008). Several studies have shown that intercropped crops use soil nutrients more efficiently than monocultured crops because of the higher recovery of N, increased yields of dry matter, and lack of negative impacts on the environment (Inal et al., 2007; Luo et al., 2016). Efficient utilization of N in belowground plant parts can promote the optimal growth of aboveground plant parts. Thus, the roles of effective planting patterns as a means of maintaining N supplies are valuable in modern agriculture.
Soil microorganisms account for a large part of the earth’s biodiversity and play a significant role in soil ecosystem biochemical processes, such as nutrient cycling and soil–borne pathogen suppression (Gomes et al., 2003; Zhang et al., 2018). Microbial diversity and composition are key determinants of their ecological functions. Most of these N-fixing microorganisms (diazotrophs) exist in free-living conditions and provide an average of 110 million tons of N per year to terrestrial ecosystems (Xu et al., 2019). In cereal–legume intercropping systems, N fixation by microorganisms associated with legumes increases the availability of N for plants, with soil N uptake by intercropped cereal greatly inducing N fixation by the root nodules in the legume rhizosphere, which stimulates N utilization by coordinating interspecific interactions (Latati et al., 2016). This practical technique influences N transformation processes by regulating microbial activities, diversities and community compositions (Lian et al., 2019; Vidal et al., 2019). Actinobacteria and Proteobacteria, two common and ubiquitous bacterial groups, are markedly affected by intercropping systems due to their biology and habitats (Stach and Bull, 2005; Gong et al., 2019b). Most N-fixing microorganisms belong to Proteobacteria and are aerobic organisms. Li et al. (2016) proposed that below-surface interactions in an maize/peanut intercropping system considerably altered the microbial structure of the soil and the dominant microbial species that are strongly linked to soil N availability. Chen et al. (2018) revealed that peanut/maize intercropping promotes plant N uptake and increases the abundance of N-cycling organisms and other beneficial rhizosphere bacteria. Lian et al. (2019) demonstrated that sugarcane/soybean intercropping in acidic soil increases microbial diversity and shifts soil microbial communities, which may stimulate N accumulation by the legume. Therefore, reasonable cereal–legume intercropping boosts soil N conversion by improving soil microbial activities and increases soil available N contents (Zuo and Zhang, 2009).
Proso millet (Panicum miliaceum L.) is one of the earliest cultivated crops in China, and it has a short growing season and is primarily planted as an important food source in the semiarid regions of China because of its high nutritional value and strong resistance to drought, saline-alkaline conditions and poor soil conditions (Seghatoleslami et al., 2008; Yang et al., 2018). Mung bean (Vigna radiata L.) is regarded as a functional food (Yao et al., 2019), and its functional components have been extracted and identified using analytical techniques. Mung bean also has a separate benefit as it engages in symbiosis with N–fixing bacteria (Choudhary and Agrawal, 2014). The intercropping combination of proso millet and mung bean has been regarded as a technically sound method based on the ecological principles of facilitation and complementarity. Thus, this practice has been rapidly applied in Northwest China (Gong et al., 2019a, 2020a, 2020b). To date, although field experiments have been conducted to evaluate the impacts of cereal–legume intercropping in recent years (Latati et al., 2016; Cao et al., 2017; Qian et al., 2018), information about crop productivity, soil microbial–mediated processes and N assimilation under proso millet–mung bean intercropping systems is still limited. Thus, we hypothesized that proso millet intercropped with mung bean will increase the soil microbial diversity, regulate the microbial community composition and efficiently contribute to N uptake and assimilation to obtain high grain yields.
The main purposes of this research were to (i) investigate the impacts of intercropping on the microbial diversity and community composition of proso millet rhizosphere soil and compare bacterial and fungal reactions to intercropping systems; (ii) explore N accumulation in soil and different plant organs and the changes in grain yield under intercropping systems; and (iii) analyze the potential microbial groups that contribute to N retention and increase grain yield. This study can provide insights into the fundamental processes of biodiversity enhancement in ecosystems for developing sustainable agriculture.
Materials and Methods
Experimental Locations
Field experiments were conducted from 2017 to 2019 at Northwest A&F University’s experimental site (37°56′26″N, 109°21′46″E) in Yulin City, Shaanxi Province, China (Supplementary Figure 1a). This area is characterized by a semiarid continental monsoon climate, and the precipitation and annual average temperature in the area are 400 mm and 8.3°C, respectively. The daily air temperature and precipitation during the crop growing seasons (May-September) in 2017, 2018, and 2019 are shown in Supplementary Figure 2. The soil pH was 8.6, the organic matter content was 7.34 g kg–1, the total N (TN) was 0.36 g kg–1, the total phosphorus (P) was 0.75 g kg–1, and the total potassium (K) was 17.88 g kg–1 of dry soil in the 0–20 cm soil layer before sowing.
Experimental Design and Treatments
The experimental design of this study was the same as that described in Gong et al. (2020a, b), and it consisted of a randomized complete block with four replicates. Four different intercropping patterns were designed as follows: 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M). Monoculture proso millet (MP) and mung bean (MM) were set as controls (Supplementary Figure 1b). All in–row distances were 0.33 m, and each experimental plot had an area of 30 m2 (6 × 5 m) and included at least three strips of proso millet and mung bean. Border rows were not used for sampling.
The cultivars ‘Shanmi–1’ (proso millet) and ‘Zhonglv–8’ (mung bean) were selected for use in this study. Proso millet was sown on 12 June and harvested on 23 September 2017, sown on 10 June and harvested on 25 September 2018, and sown on 10 June and harvested on 30 September 2019. Mung bean was sown on 28 May and harvested on 24 August 2017, sown on 18 May and harvested on 20 August 2018, and sown on 25 May and harvested on 30 August 2019. The previous season’s crops were yam, maize and potato for 2017, 2018, and 2019, respectively. The plant densities of proso millet and mung bean were 50 plants m–2 and 20 plants m–2, respectively, whereas the densities of both the monoculture crops and intercrops were the same. Basal fertilizers [120 kg (N) ha–1, 100 kg (P2O5) ha–1, and 75 kg (K2O) ha–1] were applied to the soils before planting (proso millet and mung bean) each year. During the growth period, fertilizers were not applied.
Plant Sampling and Grain Yield
Plants were sampled at the flowering stage of proso millet and at the filling stage of mung bean from 2017 to 2019. Three plants were randomly selected from the center of each plot (for a total of 12 plants) and split into stem, leaf, sheath and ear samples for proso millet and into stem, leaf, petiole and pod samples for mung bean. All aboveground samples were dried at 75°C until reaching a constant weight. The N content of different organs was measured by the Kjeldahl method after digestion with H2SO4–H2O2. At harvest, twenty proso millet plants were selected randomly in each plot (5 plants per treatment) to measure the ear length, ear number per plant, grain weight per plant and 1000-grain weight for proso millet. Similarly, the branch number per plant, pod number per plant, grain weight per plant and 100-grain weight for mung bean were investigated following the same sampling method. All the plants were harvested in each plot, and the grain yield was determined by weighing after air-drying for 2–3 weeks.
Soil Sampling and Analysis
Each soil sample was obtained from the rhizosphere of each intercropping pattern at the flowering stage of proso millet (60 days after sowing) and at the filling stage of mung bean (75 days after sowing) in August 2018 and 2019, respectively. Plants in each plot were selected using an ‘S’–shaped pattern and then homogenized to provide one composite sample per replicated site. During this sampling process, sterile paper was used to wipe the remains that were attached to the spade and sanitize the spade before collecting the next soil sample to avoid contamination between treatments and keep samples fresh. After the roots were gently shaken, the rhizosphere soils tightly attached to roots were collected and sieved through 2-mm mesh to remove stones and other residues. A portion of each soil sample was air-dried and used in the TN analysis, and another portion was stored at 4°C to determine the other N fractions. Subsamples for the molecular analyses were immediately homogenized in liquid nitrogen and stored at −80°C. Soil TN was measured as described by Li et al. (2016). Nitrate (NO3––N) and ammonium (NH4+–N) levels were analyzed based on standard methods using a continuous flow analyzer (Yu et al., 2019). Soil microbial biomass N (MBN) was determined using a chloroform fumigation–extraction method (Zhao et al., 2018). Four soil samples were analyzed per treatment.
Soil DNA Extraction, PCR Amplification and Sequencing Analysis
Microbial DNA of fresh rhizosphere soil (0.5 g) of proso millet was extracted four times (total of 2.0 g soil). The genomic DNA concentration and quality were estimated by 1.0% Sepharose gels. The bacterial 16S rRNA gene V3–V4 hypervariable region was amplified with primers 338F (ACTCCTACGGGAGGCAGCAG) and 806R (GGACTACHVGGGTWTCTAAT) (Zhang et al., 2015). A 10–digit barcode sequence was attached to the 5′ end of the forward and reverse primers in every soil sample (provided by Auwigene Company, Beijing). The PCR mixture included 4 μL of 12.5 Mm dNTP Mix, 5 μL of 10 × Ex Taq Buffer (Mg2+ plus), 200 nm of barcoded primers 16S–F and 16S–R, 1.25 U Ex Taq DNA polymerase, 2 μL of template DNA and 36.75 μL of ddH2O. The following procedure was used for PCR: an initial denaturation step for 2 min at 94°C, followed by 30 cycles of 94°C for 30 s, 57°C for 30 s and 72°C for 30 s, and a final extraction at 72°C for 5 min. The fungal internal transcribed spacer (ITS) region was amplified on an Eppendorf Mastercycler Gradient Thermocycler (Germany) with the primers ITS1F (5-CTTGGTCATTTAGAGGAAGTAA-3) and ITS2 (5-TGCGTTCTTCATCGATGC-3) (Wang et al., 2016). The 5′ ends of the two primers were tagged. Ultra–PAGE purified primers were obtained from Majorbio, China. The PCR mixture included 4 μl of 5 × FastPfu Buffer, 2.5 Mm dNTP mixture, 5 μM each primer, 2 μl of template DNA and 10 μl of H2O. The following procedure was used for PCR: an initial denaturation step for 2 min at 95°C, followed by 30 cycles of 95°C for 30 s, 55°C for 30 s and 72°C for 30 s, and a final extraction at 72°C for 5 min. Deep sequencing of bacterial and fungal samples was performed on the MiSeq platform at Allwegene Company in Beijing, China. Sequence data associated with this project have been deposited into the National Center for Biotechnology Information (NCBI) (accession numbers of PRJNA669229 for bacteria and PRJNA669216 for fungi).
Analysis of 16S rRNA and ITS Gene Data
For both bacterial and fungal reads, raw sequences were first trimmed and the reads were quality filtered, demultiplexed and processed on QIIME (Ren et al., 2018; Mcknight et al., 2019). Sequences were retained according to three criteria: (1) the barcodes and primers were explicit; (2) the length was greater than 200 bp; and (3) the quality score was higher than 30. All sequences were classified into different taxonomic groups by the Ribosomal Database Project classifier (Wang et al., 2007). The sequences, based on 97% similarity, were clustered into operational taxonomic units (OTUs) to produce rarefaction curves (Colwell and Coddington, 1994) and calculate the diversity and richness indices (Cole et al., 2008).
Statistical Analyses
The taxonomic alpha diversity, which represents the community diversity, was calculated by the Shannon index using Mothur software (v.1.30.1). The taxonomic beta diversity (the weighted UniFrac distances) illustrates the clustering of different samples and reflecting the microbial community structure, and it was determined through principal coordinate analysis (PCoA). The correlations among plant properties, soil N and soil microbial compositions were determined by a redundancy analysis (RDA) using the CANOCO 5.0 software package. We implemented a forward selection procedure according to the method described by Blanchet et al. (2008) to select a subset of plant properties, soil N and soil microbial compositions. The data were analyzed via a one-way ANOVA (SPSS 19.0, SPSS Inc., Chicago, IL, United States) under different intercropping patterns. The relationships among microbial diversity, plant properties and soil N were determined via Spearman’s correlation analysis. The differences between mean values were determined using the least significant difference (LSD) test (P < 0.05), as indicated by different letters. Moreover, the ANOVA was conducted with the standard design analysis method to determine the significance of year and treatment effects and their interactions.
Results
Grain Yield
The grain yield of proso millet under the intercropping patterns significantly increased over the three years studied (Table 1). Across all treatments and years, intercropping significantly increased the ear number per plant, ear length, grain weight per plant and 1000-grain weight compared with monoculture proso millet and resulted in grain yield improvements of 5.6–20.7% in 2017, 7.9–53.9% in 2018, and 28.3–75.4% in 2019. Among the different intercropping systems, the 2P4M treatment achieved the greatest grain yield. Compared with the monoculture mung bean, the grain yield of intercropped mung bean was 34.8–55.8% lower in 2017, 34.4–38.8% lower in 2018, and 26.5–46.5% lower in 2019 (Table 2). The maximum and minimum reductions occurred in the 4P2M and 2P4M treatments, respectively.
Table 1. Effect of intercropping on grain yield and yield components of proso millet in 2017, 2018, and 2019.
Table 2. Effect of intercropping on grain yield and yield components of mung bean in 2017, 2018, and 2019.
N Levels in the Soil-Plant System
The differences in the soil-plant system between the proso millet and mung bean strips are shown in Figures 1, 2 and Table 3. Compared with the monoculture proso millet, the average N content under the intercropping patterns in the stem, leaf, sheath, and ear tissues increased by 28.6, 15.2, 16.7, and 12.0%, respectively (Figure 1). Among the intercropping systems, the 2P4M treatment resulted in the greatest plant N improvement. Similarly, the intercropping systems achieved greater plant N contents in mung bean, which were 19.3, 13.0, 18.2, and 15.4% higher than that of the monoculture bean across all treatments and years (Figure 2).
Figure 1. Effect of intercropping patterns on soil TN (A), NO3–-N (B) NH4+–N (C), MBN (D) of proso millet. Different letters indicate significant differences (P < 0.05) among different intercropping patterns. Total nitrogen (TN), nitrate (NO3–-N), ammonium (NH4+-N), microbial biomass nitrogen (MBN). MP, 2P2M, 4P2M, 4P4M, and 2P4M represent the monoculture proso millet, 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M), respectively. * and ** Significant at the 0.05 and 0.01 probability levels, respectively. ns, no significant difference.
Figure 2. Effect of intercropping patterns on soil TN (A), NO3–-N (B) NH4+–N (C), MBN (D) of mung bean. Different letters indicate significant differences (P < 0.05) among different intercropping patterns. Total nitrogen (TN), nitrate (NO3–-N), ammonium (NH4+-N), microbial biomass nitrogen (MBN). MM, 2P2M, 4P2M, 4P4M, and 2P4M represent the monoculture mung bean, 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M), respectively. * and ** Significant at the 0.05 and 0.01 probability levels, respectively. ns, no significant difference.
Table 3. Effect of intercropping on plant N content of proso millet and mung bean in 2017, 2018, and 2019.
Intercropping greatly affected the soil N contents of these two species over the two years but showed differential effects (Table 3). For proso millet, the average TN, NO3––N, NH4+–N, and MBN were 67.7, 96.1, 23.5, and 28.1% higher under intercropping than monoculture, respectively. Among the intercropping patterns, 2P2M resulted in the maximum TN and NO3––N, which increased by 101.6 and 138.3%, respectively, while 4P4M resulted in the highest improvement in NH4+–N and MBN, which increased by 32.3% and 38.7%, respectively, compared with the MP. For mung bean, TN, NO3––N and NH4+–N were higher under intercropping than monoculture and the maximum values were achieved under the 4P4M treatment over the two years. Similarly, the large increase in MBN also corresponded to the intercropping system, and the 2P4M treatment achieved the maximum increase (55.5% higher than that under the MM treatment).
Soil Microbial Diversity
Venn graphs were constructed to evaluate the number and identity of the shared OTUs for proso millet from the five soil treatments (Figure 3). For bacteria, 1012 OTUs were jointly shared among the five treatments and 2P4M had the most OTUs that were specific to other areas (Figure 3A). For fungi, only 434 OTUs were common to the different treatments and specific OTUs increased under 4P2M and 2P4M and decreased under 2P2M and 4P4M compared with those under MP (Figure 3A).
Figure 3. Venn diagram displaying the OTU richness distribution of bacteria (A) and fungi (B) in five soil treatments. MP, 2P2M, 4P2M, 4P4M, and 2P4M represent the monoculture proso millet, 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M), respectively.
Microbial diversity under the different planting patterns was represented using an OTU-level approach. Figure 4A shows that intercropping clearly altered the soil bacterial alpha diversity (P < 0.01). The 4P4M treatment showed the highest diversity, and the MP treatment showed the lowest, with values ranging from 8.17 to 9.30. The Shannon index results displayed different changes for fungi, and the value only increased under 4P2M and decreased under 2P2M, 4P4M and 2P4M. To visualize and determine the similarities in the species composition data, the effects of different intercropping patterns on microbial community beta diversity were determined via a PCoA (Figures 4B,C). The soil bacterial communities in the 4P4M and 2P4M treatments (4 replicates of each treatment) were different from those in the 2P2M, 4P2M and MP treatments (Figure 4B). In contrast, the fungal community composition of the intercropping soil was very different from that of the MP soil, and a similar change trend was not observed (Figure 4C).
Figure 4. Effect of intercropping patterns on the soil microbial alpha diversity (Shannon index) (A) and beta diversity of proso millet (principle coordinates analysis, PCoA) (bacterial, B and fungal, C). Different letters indicate significant differences (P < 0.05) among different intercropping patterns. MP, 2P2M, 4P2M, 4P4M, and 2P4M represent the monoculture proso millet, 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M), respectively.
Soil Microbial Composition
The relative abundance of the soil bacterial communities showed seven predominant phyla (> 1%), namely, Actinobacteria (35.1%), Proteobacteria (26.5%), Chloroflexi (13.7%), Gemmatimonadetes (7.8%), Acidobacteria (8.5%), Firmicutes (2.2%) and Nitrospirae (1.4%) (Figure 5A and Supplementary Table 1). Among the microorganisms, the abundance of Proteobacteria was significantly increased under the 4P4M and 2P4M treatments. However, the average abundance of Actinobacteria was lower under the 4P4M and 2P4M treatments than the MP. Other phyla, such as Acidobacteria, Gemmatimonadetes, Chloroflexi, Nitrospirae and Firmicutes, were also altered by the cropping systems. Moreover, taxonomic classification showed that Actinobacteria was the most critical class and exhibited significant reductions under the 2P4M intercropping treatment (P < 0.05) (Supplementary Figure 3a and Supplementary Table 1). However, in Subgroup_6, the abundance of Acidobacteria was greater under the different intercropping patterns than the MP treatment. Further taxonomic classification showed that all principal bacterial groups at the order level (Supplementary Figure 4a and Supplementary Table 1) as well as Sphingomonadales, Nitrosomonadales, Myxococcales and Xanthomonadales were significantly reduced under the different intercropping patterns (P < 0.05).
Figure 5. Distribution of bacterial (A) and fungal (B) phylum communities under the intercropping patterns. MP, 2P2M, 4P2M, 4P4M, and 2P4M represent the monoculture proso millet, 2 rows of proso millet intercropped with 2 rows of mung bean (2P2M), 4 rows of proso millet intercropped with 2 rows of mung bean (4P2M), 4 rows of proso millet intercropped with 4 rows of mung bean (4P4M), and 2 rows of proso millet intercropped with 4 rows of mung bean (2P4M), respectively.
In terms of the fungal community composition in different intercropping systems across all samples (Figure 5B and Supplementary Table 2), the phyla Ascomycota, Mortierellomycota and Basidiomycota had average contributions of 64.5, 11.7, and 3.1%, respectively. In particular, the relative abundance of Basidiomycota was significantly higher in the intercropping patterns than the monoculture proso millet (P < 0.05) and could be ranked as 4P2M > 2P2M > 2P4M > 4P4M > MP. Conversely, Ascomycota and Mortierellomycota showed different change trends and were significantly affected by intercropping (P < 0.01). Within Ascomycota, the classes Sordariomycetes and Eurotiomycetes were the most abundant and had mean relative abundances of 32.06 and 15.63%, respectively (Supplementary Figure 3b and Supplementary Table 2). At the order level, the abundance of Microascales was significantly higher than that of Xylariales and Onygenales, and the abundance of Pezizales was significantly lower under intercropping than MP (P < 0.01) (Supplementary Figure 4b and Supplementary Table 2).
Relationship Among Plant Properties, Soil N and Soil Microbial Structure
Through the Spearman test method, the top 20 most abundant OTUs of the bacterial and fungal phyla were selected for correlation analysis. As shown in Figure 6, the size of the points represents the magnitude of phyla abundance while the thickness of the line represents the correlation size. The red line indicated a positive correlation, and the blue line showed a negative correlation. The figure shows that there was a strong positive correlation among the top four dominant groups of bacteria, namely, included Proteobacteria, Gemmatimonadetes, Chloroflexi, and Actinobacteria. Moreover, the three dominant fungal communities Mortierellomycota, Basidiomycota, and Ascomycota also showed strong correlations. These organisms play a vital role in N cycling in soil-plant systems.
Figure 6. Network interaction diagram of dominant bacterial (A) and fungal (B) groups. Red lines indicate a positive correlation while the blue lines show a negative correlation. The size of the points represents the magnitude of phyla abundance while the thickness of the line represents the correlation size.
The statistical analysis indicated that the bacterial alpha and beta diversities were markedly and positively related to the soil NO3––N, grain weight per plant, and grain yield, whereas only the fungal alpha diversity showed similar positive correlations (Table 4). Moreover, a RDA was conducted to quantify the relative influence of environmental parameters (four soil N properties, namely, TN, NO3––N, NH4+–N and MBN; and five plant properties, namely, stem N, leaf N, sheath N, ear N, and GY) on the microbial community composition at distinct taxonomic levels (phylum, class and order) (Figure 6 and Supplementary Figures 5, 6). The results showed that the environmental variables significantly affected the microbial community. In particular, TN and NO3––N were significantly associated with the changes in Proteobacteria and Firmicutes while NH4+–N and MBN were significantly associated with the changes in Chloroflexi and Acidobacteria, respectively. Except for the stem N and sheath N, other plant properties were significantly and negatively associated with the abundance of Actinobacteria, Gemmatimonadetes and Nitrospirae at the bacterial phylum level (Figure 7A). For Proteobacteria, the soil properties and plant properties were correlated with the changes in the abundance of Sphingomonadales, Nitrosomonadales, and Xanthomonadales (Supplementary Figure 6a), which belong to Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria, respectively (Supplementary Table 1). Moreover, Figure 7B shows that TN, NO3––N and MBN among the soil properties and stem N, and GY among the plant properties influenced the abundance of Ascomycota and Mortierellomycota. For Ascomycota, the abundance of the class Pezizomycetes and the orders Mortierellales, Glomerellales and Sordariales orders were sensitive to changes in soil and plant properties (Supplementary Figures 5b, 6b).
Table 4. Spearman’s rank correlation coefficients (R) between microbial diversity (i.e., alpha diversity: Shannon index; beta diversity: PCoA) and plant properties and soil N.
Figure 7. Ordination plots of the results from the redundancy analysis (RDA) to identify the relationships among the bacterial (A) and fungal (B) taxa (blue arrows) and the plant properties and soil N (red arrows) at the phylum level. Bacterial taxa: Actinobacteria (Acti), Proteobacteria (Prot), Chloroflexi (Chlo), Gemmatimonadetes (Gemm), Acidobacteria (Acid), Nitrospirae (Nitr), Firmicutes (Firm). Fungal taxa: Ascomycota (Asco), Mortierellomycota (Mort), Basidiomycota (Basi). Plant properties: Grain yield (GY), stem nitrogen (Stem N), leaf nitrogen (Leaf N), sheath nitrogen (Sheath N), ear nitrogen (Ear N). Soil N: Total nitrogen (TN), nitrate (NO3–-N), ammonium (NH4+–N), microbial biomass nitrogen (MBN).
Discussion
Effect of Intercropping on Crop Grain Yield
Positive interactions occur in complex symbiotic systems that enhance the growth of crops, especially in cereal–legume intercropping systems, and these interactions are beneficial for improving soil fertility and yield (Hauggaard-Nielsen et al., 2001; Qian et al., 2018). In this study, intercropping led to a significant improvement in proso millet grain yield and variations in the ear number per plant, ear length, grain weight per plant and 1000-grain weight, and the grain yield increasing by 5.6–20.7% in 2017, 7.9–53.9% in 2018, and 28.3–75.4% in 2019 (Table 1). In contrast, the grain yield of intercropped mung bean were lower than those under the MM treatment at 34.8–55.8% in 2017, 34.4–38.8% in 2018, and 26.5–46.5% in 2019 (Table 2). These results could be explained by the spatial canopy structure, which resulted in greater light capture for the taller proso millet and less incident radiation for the shorter mung bean (Gong et al., 2020b). In intercropping systems, two species compete for light and their competitive abilities shift with their niches. For the production performance of intercropped proso millet, the greater light area and intensity from the margin effects were directly increased under intercropping by the combination of tall and short species, resulting in an advantage of light energy utilization (Wang et al., 2015); moreover, these helpful effects might facilitate interactions among the soil environment (such as the nutrient content), microbes and plants (Zhao et al., 2017). This theory is supported by the close relationship among the soil microbial composition, soil biochemical parameters, plant N and grain yield. Such improvements are caused by the N2 fixation of legumes, which maximizes and sustains plant growth for the intercropped proso millet.
Effect of Intercropping on N Assimilation in the Soil–Plant System
Agricultural practices (i.e., tillage, cultivation, and fertilization) modify soil physicochemical properties and consequently alter soil nutrient conditions (Zhou et al., 2011; Jia et al., 2020). The directions and magnitudes of such differential responses are related to soil N because N is essential for plant growth and greatly affected by land use changes (Yang et al., 2017). The present results found that the TN, NO3––N, NH4+–N and MBN contents were significantly increased under the intercropping systems compared with the monoculture systems, including those of proso millet (Figure 1) and mung bean (Figure 2). This finding may be attributed to the complex biological diversity under intercropping systems that results in the transfer of N to soil via ions and root exudates and further facilitates the accumulation and decomposition of soil N fractions. In legume–mixed intercropping, legumes increase N2 fixation to provide higher N levels for the N utilization of the adjacent crop, thus yielding a growth advantage for the intercropped plant (Subedi and Ma, 2005; Zhao et al., 2017). Interestingly, N is an important resource for cereals and the N fixed by legumes can be used by intercropped cereals during their growing period (Shen and Chu, 2004). Soil N fractions are presented as N-containing compounds, such as nucleic acids, amino sugars, and amino acids, which vary with the root structure, nutritional quality, and litter deposition. Therefore, based on the significant effect of interspecific interactions on N uptake under intercropping, proso millet–mung bean intercrops increase the soil N supply following mung bean N2 fixation in the rhizosphere. Moreover, soil microbes may change N recycling and utilization processes and could transform N to many different N forms, including mobile forms, to promote plant N assimilation (Vidal et al., 2019). Duchene et al. (2017) observed that rhizosphere N availability was influenced by intricate microbial activities under intercropping treatments. Increases in soil N forms, such as TN, NO3––N, NH4+–N and MBN, could be produced via the high activity of N–fixing bacteria and ammonia–oxidizing bacteria under intercropping systems (Ramirez et al., 2012). Furthermore, belowground interspecific interactions may influence aboveground nutrient accumulation in cereal–legume intercropping systems. Compared with the N contents of monoculture crops, under the intercropping patterns, the average N contents in the stem, leaf, sheath, and ear tissues of proso millet increased by 28.6, 15.2, 16.7, and 12.0%, respectively, while the average N contents in the stem, leaf, petiole, and pod tissues of mung bean increased by 19.3, 13.0, 18.2, and 15.4%, respectively. This phenomenon indicated that higher soil nutrients may be expected to result in higher plant nutrients, which in turn results in the exudation of a more diverse range of organic compounds into the soil, thereby facilitating N assimilation. Additionally, it is generally accepted that belowground microbial communities can affect aboveground plant nutrient conditions by carrying out a wide spectrum of decomposition and metabolic processes. Thus, these findings indicated that intercropping could produce a large amount of N in rhizosphere soil and increase nutrient accumulation. Moreover, the habitats are fragile and the soil is barren and N limited on the Loess Plateau. Thus, intercropping systems alleviated N limitation and improved soil fertility in the present study.
Intercropping Altered the Soil Microbial Diversity and Community to Facilitate N Assimilation
The valuable influences of intercropping in ameliorating biodiversity, heightening farmland ecosystem balance, and diminishing the occurrence of harmful organisms have been well documented (LaMondia et al., 2002; Ning et al., 2017). The soil microbial community is a fundamental component of soil quality and is imperative for many ecological processes, such as energy flow, nutrient cycling, and organic matter turnover (Acosta-Martínez et al., 2008). Our sequencing results showed that intercropping increased the soil bacterial alpha diversity (Shannon index), especially under the 4P4M and 2P4M treatments (Figure 4A), which pointed to the invulnerable soil microbes in this agricultural ecosystem. The Shannon index accounts for species richness abundance and evenness of the species present in the sample (Shannon, 1948). Thus, the greater increase in alpha diversity under intercropping suggested that rare species might have been more abundant. The changes in the bacterial community could affect soil N availability, and this conclusion was verified by the significant positive effect of NO3––N on bacterial alpha diversity (Table 4). Soil microorganisms are pivotal players in driving biogeochemical cycles (Palomo et al., 2016); thus, the increased microbial diversity and N under intercropping would influence ecosystem processes and functions. Microbial biomass residues have been identified as significant sources of soil organic matter (Simpson et al., 2007; Miltner et al., 2012). The increase in soil microbial diversity and N may influence both the pool size and the chemical composition of soil carbon (Liang et al., 2011). In addition, links between microbial diversity and plant productivity, plant diversity, or nutrient acquisition have also been observed (Zak et al., 2003). Soil fungal alpha diversity was less affected by intercropping patterns compared with the bacteria, which may have been caused by the number of copiotrophic characteristics of bacteria (Ramirez et al., 2012). The nutrient supply under intercropping systems was more abundant than that under traditional monocultures and could provide more sufficient nutrients for bacterial growth. Although Ren et al. (2017b) found that the soil bacterial and fungal community characteristics were affected by altitude, they showed that soil bacterial alpha diversity fluctuated more along altitudinal gradients than fungal alpha diversity, which was likely because of soil N. N metabolism in soil is overwhelmingly associated with bacterial activity. Thus, a greater abundance of bacteria in the intercropping system illustrates that bacteria may play a crucial role in facilitating N metabolism in the rhizosphere. Compared with monocots, dicots, especially legumes, generate and secrete more organic compounds into the crop rhizosphere (Raghothama, 1999). Importantly, increases in the organic compound content in root exudates can provide more carbon sources for rhizosphere microorganism growth, which may be another reason for the change in the bacterial community in the intercropping treatments. These results illustrated that the intercropping treatment had some valuable effects on soil microbial activity, which in turn could facilitate the nutrient supply for plants as previously reported (Hao et al., 2003).
Microbial communities play an active role in the maintenance of soil functions in ecosystem processes (Van Der Heijden et al., 2008), such as decomposing organic matter and transforming inorganic N to organic N (Leiros et al., 2000; Prosser and Nicol, 2008). The soil bacterial community compositions revealed that the 4P4M and 2P4M treatments significantly increased the abundance of the phylum Proteobacteria. However, although Actinobacteria was the most abundant bacterial phylum, the average abundance of Actinobacteria was lower with intercropping under the same conditions compared with that under MP (Figure 5A). The change in the abundance of bacteria may be in connection with the habitat for low-nutrient consumption groups that are well adapted to unstable N levels in the soil (Stach and Bull, 2005). Chen et al. (2003) found that N fixation by Betaproteobacteria in the legume symbiotic system is widespread in nature. Our study revealed that TN and NO3––N were significantly positively related to Proteobacteria (Figure 7A). Consequently, N assimilation was enhanced by Proteobacteria. Furthermore, Sphingomonadales and Xanthomonadales are members of the Alphaproteobacteria and Gammaproteobacteria classes (Supplementary Table 1) and heterotrophic and N–fixing organisms. Sphingomonadales are also considered rhizosphere plant–promoting bacteria and can fix atmospheric N2 in symbiosis with plants (Regupathy et al., 2012). Hence, the increase in the abundance of Proteobacteria was beneficial for N accumulation and cycling. Moreover, other bacterial taxa, such as Gemmatimonadetes and Nitrospirae, exhibited strong responses to changes in plant N and soil N, and the changes in these two phyla were likely due to their own ecological characteristics (Ivanova et al., 2016). Zhang et al. (2015) observed that the abundance of ammonia oxidizers decreases in the crop rhizospheres under maize–faba bean intercropping, and this finding indirectly supports our explanation. The positive correlations among the dominant bacterial communities are shown in Figure 6A and indicate that under suitable conditions, these groups grow together to ensure their individual and collective survival and dominance over most other organisms in the soil. Altogether, proso millet intercropped with mung bean greatly affected the dominant bacterial communities, and the functions of these communities can provide feedback to facilitate N assimilation.
For the soil fungal community composition, the abundances of both Ascomycota and Mortierellomycota were 75% (Figure 5B), which is consistent with reports of these taxa on a global scale for arid farmlands (Li and Wu, 2018). However, these two dominant phyla were not responsive to environmentally induced changes in intercropping conditions and only Ascomycota was affected by TN and NO3––N (Figure 7B). This phenomenon could be explained by the functions of Ascomycota. The abundance of Ascomycota was increased by residue degradation, and these fungi can rapidly metabolize organic substrates of rhizodeposition in rhizosphere soil (Bastida et al., 2013). Hence, proso millet–mung bean intercropping may create suitable circumstances for Ascomycota and allow them to better exploit the easily degradable fraction of plant residues and facilitate N accumulation. Moreover, the growth rates of Ascomycota were increased by N availability, thereby expediting the decomposition of plant residues (Fontaine et al., 2011). In addition, according to the growth rate of the microorganisms (GRH; Elser et al., 1996; Ren et al., 2017a), the growth of most species is associated with a greater demand for N for the synthesis of ribosomal DNA and protein. Therefore, the relative abundance of this dominant microbial phylum was very sensitive to changes in soil N dynamics and availability. However, at the order level, Hypocreales (a member of Sordariomycetes in Ascomycota), which is related to the soil N fractions and plant N accumulation (Supplementary Figure 6), showed significant decreasing trends under the intercropping pattern (Supplementary Table 2) because Hypocreales is a rapidly growing plant tissue decomposer (Hannula et al., 2012) and the lower amount of litter at the proso millet flowering stage under intercropping led to a decrease in the abundance of this community. Moreover, studies have also shown that many species belonging to the order Hypocreales are insect pathogens (Castrillo et al., 2008; Riosvelasco et al., 2014) and increasing their abundance would contribute to the control of pests and diseases under field conditions. Furthermore, combined with the composition of the bacterial community, these discoveries illustrated that the microbial diversity and community composition (bacterial and fungal) in rhizosphere soil were altered under proso millet/mung bean intercropping, which promoted plant N assimilation and increased grain yield.
Conclusion
Our data clearly indicated that proso millet/mung bean intercropping altered the abundance of soil bacteria and fungi. The soil bacterial alpha diversity was higher under intercropping patterns and varied more than the fungal alpha diversity. With regard to the bacterial composition, the phylum Proteobacteria, which is positively related to soil N accumulation, significantly increased. The changes in the composition of Proteobacteria, such as in the orders Sphingomonadales and Xanthomonadales, influenced N assimilation. For the fungal community composition, Ascomycota was the principal phylum and was affected by TN and NO3––N. Microbial changes enhanced N uptake in the soil and promoted N accumulation in different plant organs, resulting in an improvement in grain yield. The average grain yield under the 2P4M intercropping treatment increased by 50.1% compared with that of the monoculture proso millet, showing the highest productivity. Our study shows that the intercropping of proso millet with mung bean influenced the microbial (bacterial and fungal) communities and provides insight into the roles of microbial biodiversity and ecological performance in improving crop production.
Data Availability Statement
Sequence data associated with this project have been deposited into the National Center for Biotechnology Information (NCBI) (accession numbers of PRJNA669229 for bacteria and PRJNA670554 for fungi).
Author Contributions
KD and XG performed the experiments and drafted the manuscript. AI and BF designed the project and guided the experiments. GZ and HW investigated the material characteristics. KD and XG analyzed the data and planted the material. BF organized and coordinated the whole project. All authors contributed to the article and approved the submitted version.
Funding
The research was supported by the National Natural Science Foundation of China (31371529), National Millet Crops Research and Development System (CARS–06–13.5–A26), National Science and Technology Supporting Plan (2014BAD07B03), and Minor Grain Crops Research and Development System of Shaanxi Province (NYKJ–2018–YL19). We are grateful to the reviewers for their constructive suggestions.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.601054/full#supplementary-material
References
Acosta-Martínez, V., Acosta-Mercado, D., Sotomayor-Ramírez, D., and Cruz-Rodríguez, L. (2008). Microbial communities and enzymatic activities under different management in semiarid soils. Appl. Soil Ecol. 38, 249–260. doi: 10.1016/j.apsoil.2007.10.012
Banik, P., and Sharma, R. C. (2009). Yield and resource utilization efficiency in baby corn—legume–intercropping system in the Eastern Plateau of India. J. Sustain. Agr. 33, 379–395. doi: 10.1080/10440040902834970
Bastida, F., Hernández, T., Albaladejo, J., and García, C. (2013). Phylogenetic and functional changes in the microbial community of long–term restored soils under semiarid climate. Soil Biol. Biochem. 65, 12–21. doi: 10.1016/j.soilbio.2013.04.022
Blanchet, F. G., Legendre, P., and Borcard, D. (2008). Forward selection of explanatory variables. Ecology 89, 2623–2632. doi: 10.1890/07-0986.1
Boardman, J., Poesen, J., and Evans, R. (2003). Socio–economic factors in soil erosion and conservation. Environ. Sci. Policy 6, 1–6. doi: 10.1016/S1462-9011(02)00120-X
Cao, X., Liu, S., Wang, J., Wang, H., Chen, L., Tian, X., et al. (2017). Soil bacterial diversity changes in different broomcorn millet intercropping systems. J. Basic Microb. 57, 989–997. doi: 10.1002/jobm.201700133
Castrillo, L. A., Griggs, M. H., and Vandenberg, J. D. (2008). Quantitative detection of Beauveria bassinan GHA (Ascomycota: Hypocreales), a potential microbial control agent of the emerald ash borer, by use of real-time PCR. Biol. Control 45, 163–169. doi: 10.1016/j.biocontrol.2008.01.005
Chen, J., Arafat, Y., Wu, L., Xiao, Z., Li, Q., Khan, M. A., et al. (2018). Shifts in soil microbial community, soil enzymes and crop yield under peanut/maize intercropping with reduced nitrogen levels. Appl. Soil Ecol. 124, 327–334. doi: 10.1016/j.apsoil.2017.11.010
Chen, W. M., Moulin, L., Bontemps, C., Vandamme, P., and Boivinmasson, C. (2003). Legume symbiotic nitrogen fixation by β–proteobacteria is widespread in nature. J. Bacteriol. 185, 7266–7272. doi: 10.1128/JB.185.24.7266-7272.2003
Choudhary, K. K., and Agrawal, S. (2014). Cultivar specificity of tropical mung bean (Vigna radiata L.) to elevated ultraviolet–B: changes in antioxidative defense system, nitrogen metabolism and accumulation of jasmonic and salicylic acids. Environ. Exp. Bot. 99, 122–132. doi: 10.1016/j.envexpbot.2013.11.006
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2008). The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, 141–145. doi: 10.1093/nar/gkn879
Colwell, R. K., and Coddington, J. A. (1994). Estimating terrestrial biodiversity through extrapolation. Philos. Trans. R. Soc. B 345, 101–118. doi: 10.2307/56143
Duchene, O., Vian, J. F., and Celette, F. (2017). Intercropping with legume for agroecological cropping systems: complementarity and facilitation processes and the importance of soil microorganisms. A review. Agr. Ecosyst. Environ. 240, 148–161. doi: 10.1016/j.agee.2017.02.019
Elser, J. J., Dobberfuhl, D. R., MacKay, N. A., and Schampel, J. H. (1996). Organism size, life history, and N:P stoichiometry toward a unified view of cellular and ecosystem processes. Bioscience 46, 674–684. doi: 10.2307/1312897
Fontaine, S., Henault, C., Aamor, A., Bdioui, N., Bloor, J., Maire, V., et al. (2011). Fungi mediate long term sequestration of carbon and nitrogen in soil through their priming effect. Soil Biol. Biochem. 43, 86–96. doi: 10.1016/j.soilbio.2010.09.017
Ghosh, P. K., Tripathi, A. K., Bandyopadhyay, K. K., and Manna, M. C. (2009). Assessment of nutrient competition and nutrient requirement in soybean/sorghum intercropping system. Eur. J. Agron. 31, 43–50. doi: 10.1016/j.eja.2009.03.002
Gomes, N. C. M., Fagbola, O., Costa, R., Rumjanek, N. G., Buchner, A., Mendona-Hagler, L., et al. (2003). Dynamics of fungal communities in bulk and maize rhizosphere soil in the tropics. Appl. Environ. Microbiol. 69, 3758–3766. doi: 10.1128/AEM.69.7.3758-3766.2003
Gong, X., Dang, K., Lv, S., Zhao, G., Tian, L., Luo, Y., et al. (2020a). Interspecific root interactions and water-use efficiency of intercropped proso millet and mung bean. Eur. J. Agron. 115:126034. doi: 10.1016/j.eja.2020.126034
Gong, X., Ferdinand, U., Dang, K., Li, J., Chen, G., Luo, Y., et al. (2020b). Boosting proso millet yield by altering canopy light distribution in proso millet/mung bean intercropping systems. Crop J. 8, 365–377. doi: 10.1016/j.cj.2019.09.009
Gong, X., Liu, C., Ferdinand, U., Dang, K., Zhao, G., Yang, P., et al. (2019a). Effect of intercropping on leaf senescence related to physiological metabolism in proso millet (Panicum miliaceum L.). Photosynthetica 57, 993–1006. doi: 10.32615/ps.2019.112
Gong, X., Liu, C., Li, J., Luo, Y., Yang, Q., Zhang, W., et al. (2019b). Responses of rhizosphere soil properties, enzyme activities and microbial diversity to intercropping patterns on the Loess Plateau of China. Soil Till. Res. 195:104355. doi: 10.1016/j.still.2019.104355
Hannula, S. E., Boschker, H. T., Boer, W., and Veen, J. A. (2012). 13C pulse-labeling assessment of the community structure of active fungi in the rhizosphere of a genetically starch-modified potato (Solanum tuberosum) cultivar and its parental isoline. New Phytol. 194, 784–799. doi: 10.1111/j.1469-8137.2012.04089.x
Hao, Y. R., Lao, X. R., Sun, W. H., and Peng, S. L. (2003). Interaction of roots and rhizosphere in the wheat/maize intercropping system. Rural Eco-Environ. 19, 18–22.
Hauggaard-Nielsen, H., Ambus, P., and Jensen, E. S. (2001). Interspecific competition, N use and interference with weeds in pea–barley intercropping. Field Crops Res. 70, 101–109. doi: 10.1016/S0378-4290(01)00126-5
Inal, A., Gunes, A., Zhang, F., and Cakmak, I. (2007). Peanut/maize intercropping induced changes in rhizosphere and nutrient concentrations in shoots. Plant Physiol. Biochem. 45, 350–356. doi: 10.1016/j.plaphy.2007.03.016
Ivanova, A. A., Kulichevskaya, I. S., Merkel, A. Y., Toshchakov, S. V., and Dedysh, S. N. (2016). High diversity of Planctomycetes in soils of two lichen-dominated sub-Arctic ecosystems of northwestern Siberia. Front. Microbiol. 7:2065. doi: 10.3389/fmicb.2016.02065
Jacobsen, S. E., Sørensen, M., Pedersen, S. M., and Weiner, J. (2013). Feeding the world: genetically modified crops versus agricultural biodiversity. Agron. Sustain. Dev. 33, 651–662. doi: 10.1007/s13593-013-0138-9
Jia, X., Zhong, Y., Liu, J., Zhu, G., Shangguan, Z., and Yan, W. (2020). Effects of nitrogen enrichment on soil microbial characteristics: from biomass to enzyme activities. Geoderma 366:114256. doi: 10.1016/j.geoderma.2020.114256
LaMondia, J. A., Elmer, W. H., Mervosh, T. L., and Cowles, R. S. (2002). Integrated management of strawberry pests by rotation and intercropping. Crop Prot. 21, 837–846. doi: 10.1016/S0261-2194(02)00050-9
Latati, M., Bargaz, A., Belarbi, B., Lazali, M., Benlahrech, S., Tellah, S., et al. (2016). The intercropping common bean with maize improves the rhizobial efficiency, resource use and grain yield under low phosphorus availability. Eur. J. Agron. 72, 80–90. doi: 10.1016/j.eja.2015.09.015
Leiros, M. C., Trasarcepeda, C., Seoane, S., and Gilsotres, F. (2000). Biochemical properties of acid soils under climax vegetation (Atlantic oakwood) in an area of the European temperate–humid zone (Galicia, NW Spain): general parameters. Soil Biol. Biochem. 32, 733–745. doi: 10.1016/S0038-0717(99)00195-9
Li, L., Sun, J., Zhang, F., Li, X., Rengel, Z., and Yang, S. (2001). Wheat/maize or wheat/soybean strip intercropping: II. Recovery or compensation of maize and soybean after wheat harvesting. Field Crops Res. 71, 173–181. doi: 10.1016/S0378-4290(01)00157-5
Li, Q., Wu, L., Jun, C., Khan, M. A., Luo, X., and Lin, W. (2016). Biochemical and microbial properties of rhizospheres under maize/peanut intercropping. J. Integr. Agr. 15, 101–110. doi: 10.1016/S2095-3119(15)61089-9
Li, S., and Wu, F. (2018). Diversity and cooccurrence patterns of soil bacterial and fungal communities in seven intercropping systems. Front. Microbiol. 9:1521. doi: 10.3389/fmicb.2018.01521
Lian, T., Mu, Y., Jin, J., Ma, Q., Cheng, Y., Cai, Z., et al. (2019). Impact of intercropping on the coupling between soil microbial community structure, activity, and nutrient–use efficiencies. Peer J. 7:e6412. doi: 10.7717/peerj.6412
Liang, C., Cheng, G., Wixon, D. L., and Balser, T. C. (2011). An Absorbing Markov chain approach to understanding the microbial role in soil carbon stabilization. Biogeochemistry 106, 303–309. doi: 10.1007/s10533-010-9525-3
Luo, S., Yu, L., Liu, Y., Zhang, Y., Yang, W., Li, Z., et al. (2016). Effects of reduced nitrogen input on productivity and N2O emissions in a sugarcane/soybean intercropping system. Eur. J. Agron. 81, 78–85. doi: 10.1016/j.eja.2016.09.002
Martin-Guay, M. O., Paquette, A., Dupras, A., and Rivest, D. (2017). The new green revolution: sustainable intensification of agriculture by intercropping. Sci. Total Environ. 615, 767–772. doi: 10.1016/j.scitotenv.2017.10.024
Mcknight, D. T., Huerlimann, R., Bower, D. S., Schwarzkopf, L., Alford, R. A., and Zenger, K. R. (2019). Methods for normalizing microbiome data: an ecological perspective. Methods Ecol. Evol. 10, 389–400. doi: 10.1111/2041-210X.13115
Miltner, A., Bombach, P., Schmidt-Brucken, B., and Kastner, M. (2012). SOM genesis: microbial biomass as a significant source. Biogeochemistry 111, 41–55. doi: 10.1007/s10533-011-9658-z
Ning, C., Qu, J., He, L., Yang, R., Chen, Q., Luo, S., et al. (2017). Improvement of yield, pest control and Si nutrition of rice by rice–water spinach intercropping. Field Crops Res. 208, 34–43. doi: 10.1016/j.fcr.2017.04.005
Palomo, A., Jane Fowler, S., Gulay, A., Rasmussen, S., Sicheritz-Ponten, T., and Smets, B. F. (2016). Metagenomic analysis of rapid gravity sand filter microbial communities suggests novel physiology of Nitrospira spp. ISME J. 10, 2569–2581. doi: 10.1038/ismej.2016.63
Prosser, J. I., and Nicol, G. W. (2008). Relative contributions of archaea and bacteria to aerobic ammonia oxidation in the environment. Environ. Microbiol. 10, 2931–2941. doi: 10.1111/j.1462-2920.2008.01775.x
Qian, X., Zang, H., Xu, H., Hu, Y., Ren, C., Guo, L., et al. (2018). Relay strip intercropping of oat with maize, sunflower and mung bean in semi–arid regions of Northeast China: yield advantages and economic benefits. Field Crops Res. 223, 33–40. doi: 10.1016/j.fcr.2018.04.004
Ramirez, K. S., Craine, J. M., and Fierer, N. (2012). Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. Global Change Biol. 18, 1918–1927. doi: 10.1111/j.1365-2486.2012.02639.x
Regupathy, T. V., Sun, H. L., Young, J. Y., and Young, H. R. (2012). Analysis of bacterial community in the ginseng soil using denaturing gradient gel electrophoresis (DGGE). Indian J. Microbiol. 52, 286–288. doi: 10.1007/s12088-011-0193-3
Ren, C., Chen, J., Deng, J., Zhao, F., Han, X., Yang, G., et al. (2017a). Response of microbial diversity to C:N:P stoichiometry in fine root and microbial biomass following afforestation. Biol. Fertil. Soils 53, 457–468. doi: 10.1007/s00374-017-1197-x
Ren, C., Wang, T., Xu, Y., Deng, J., Zhao, F., Yang, G., et al. (2018). Differential soil microbial community responses to the linkage of soil organic carbon fractions with respiration across land–use changes. Forest Ecol. Manag. 409, 170–178. doi: 10.1016/j.foreco.2017.11.011
Ren, C., Zhang, W., Zhong, Z. K., Han, X., Yang, G., Feng, Y., et al. (2017b). Differential responses of soil microbial biomass, diversity, and compositions to altitudinal gradients depend on plant and soil characteristics. Sci. Total Environ. 750, 610–611. doi: 10.1016/j.scitotenv.2017.08.110
Riosvelasco, C., Pérezcorral, D. A., Salasmarina, M., Berlanga-Reyes, D., Ornelas-Paz, J., Acosta Muñiz, C., et al. (2014). Pathogenicity of the Hypocreales Fungi Beauveria bassiana and Metarhizium anisopliae against insect pests of tomato. Southwestern Entomol. 39, 739–750. doi: 10.3958/059.039.0405
Salvagiotti, F., Cassman, K. G., Specht, J. E., Walters, D. T., Weiss, A., and Dobermann, A. (2008). Nitrogen uptake, fixation and response to fertilizer N in soybeans: a review. Field Crops Res. 108, 1–13. doi: 10.1016/j.fcr.2008.03.001
Seghatoleslami, M., Kafi, M., and Majidi, E. (2008). Effect of drought stress at different growth stages on yield and water use efficiency of five proso millet (Panicum miliaceum L.) genotypes. Pak. J. Bot. 40, 1427–1432. doi: 10.1127/0029-5035/2008/0087-0261
Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423. doi: 10.1109/9780470544242.ch1
Shen, Q., and Chu, G. (2004). Bi–directional nitrogen transfer in an intercropping system of peanut with rice cultivated in aerobic soil. Biol. Fertil. Soils 40, 81–87. doi: 10.1007/s00374-004-0737-3
Simpson, A. J., Simpson, M. J., Smith, E., and Kelleher, B. P. (2007). Microbially derived inputs to soil organic matter: are current estimates too low? Environ. Sci. Technol. 41, 8070–8076. doi: 10.1021/es071217x
Stach, E., and Bull, A. T. (2005). Estimating and comparing the diversity of marine actinobacteria. Antonie Van Leeuwenhoek. 87, 3–9. doi: 10.1007/s10482-004-6524-1
Subedi, K., and Ma, B. (2005). Effects of N–deficiency and timing of N supply on the recovery and distribution of labeled 15N in contrasting maize hybrids. Plant Soil 273, 189–202. doi: 10.1007/s11104-004-7540-7
Van Der Heijden, M. G., Bardgett, R. D., and Van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
Vidal, D. F., Trichet, P., Puzos, L., Bakker, M. R., Delerue, F., and Augusto, L. (2019). Intercropping N–fixing shrubs in pine plantation forestry as an ecologically sustainable management option. Forest Ecol. Manag. 437, 175–187. doi: 10.1016/j.foreco.2019.01.023
Wang, J., Du, X., Zhang, Y., Li, T., and Liao, X. (2016). Effect of substrate on identification of microbial communities in poultry carcass composting and microorganisms associated with poultry carcass decomposition. J. Agr. Food Chem. 64, 6838–6847. doi: 10.1021/acs.jafc.6b02442
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, Q., Sun, D., Hao, H., Zhao, X., Hao, W., and Liu, Q. (2015). Photosynthetically active radiation determining yields for an intercrop of maize with cabbage. Eur. J. Agron. 69, 32–40. doi: 10.1016/j.eja.2015.05.004
Xu, Y., Wang, T., Li, H., Ren, C., Chen, J., Yang, G., et al. (2019). Variations of soil nitrogen-fixing microorganism communities and nitrogen fractions in a Robinia pseudoacacia chronosequence on the Loess Plateau of China. Catena 174, 316–323. doi: 10.1016/j.catena.2018.11.009
Yang, F., Liao, D., Wu, X., Gao, R., Fan, Y., Raza, M. A., et al. (2017). Effect of aboveground and belowground interactions on the intercrop yields in maize–soybean relay intercropping systems. Field Crops Res. 203, 16–23. doi: 10.1016/j.fcr.2016.12.007
Yang, Q., Zhang, P., Qu, Y., Gao, X., Liang, J., Yang, P., et al. (2018). Comparison of physicochemical properties and cooking edibility of waxy and non–waxy proso millet (Panicum miliaceum L.). Food Chem. 257, 271–278. doi: 10.1016/j.foodchem.2018.03.009
Yao, M., Tian, Y., Yang, W., Huang, M., Zhou, S., and Liu, X. (2019). The multi–scale structure, thermal and digestion properties of mung bean starch. Int. J. Biol. Macromol. 131, 871–878. doi: 10.1016/j.ijbiomac.2019.03.102
Yu, H., Heerink, N., Jin, S., Berentsen, P., Zhang, L., and Werf, W. V. D. (2017). Intercropping and agroforestry in China–current state and trends. Agr. Ecosyst. Environ. 244, 52–61. doi: 10.1016/j.agee.2017.04.019
Yu, N., Zhao, C., Ma, B., Li, S., She, Z., Guo, L., et al. (2019). Impact of ampicillin on the nitrogen removal, microbial community and enzymatic activity of activated sludge. Bioresource Technol. 272, 337–345. doi: 10.1016/j.biortech.2018.10.048
Zak, D. R., Holmes, W. E., White, D. C., Peacock, A. D., and Tilman, D. (2003). Plant diversity, soil microbial communities, and ecosystem function, are there any links? Ecology 84, 2042–2050. doi: 10.1890/02-0433
Zhang, K., Cheng, X., Xiao, S., Yi, L., and Zhang, Q. (2018). Linking soil bacterial and fungal communities to vegetation succession following agricultural abandonment. Plant Soil 431, 19–36. doi: 10.1007/s11104-018-3743-1
Zhang, S., Jin, Y., Zeng, Z., Liu, Z., and Fu, Z. (2015). Subchronic exposure of mice to cadmium perturbs their hepatic energy metabolism and gut microbiome. Chem. Res. Toxicol. 28, 2000–2009. doi: 10.1021/acs.chemrestox.5b00237
Zhao, F., Ren, C., Han, X., Yang, G., Wang, J., and Doughty, R. (2018). Changes of soil microbial and enzyme activities are linked to soil C, N and P stoichiometry in afforested ecosystems. Forest Ecol. Manag. 427, 289–295. doi: 10.1016/j.foreco.2018.06.011
Zhao, M., Jones, C. M., Meijer, J., Lundquist, P. O., Fransson, P., Carlsson, G., et al. (2017). Intercropping affects genetic potential for inorganic nitrogen cycling by root–associated microorganisms in Medicago sativa and Dactylis glomerata. Appl. Soil Ecol. 119, 260–266. doi: 10.1016/j.apsoil.2017.06.040
Zhou, X., Yu, G., and Wu, F. (2011). Effects of intercropping cucumber with onion or garlic on soil enzyme activities, microbial communities and cucumber yield. Eur. J. Soil Biol. 47, 279–287. doi: 10.1016/j.ejsobi.2011.07.001
Keywords: cereal–legume intercropping, soil microbial community, soil–plant system, N availability, grain yield
Citation: Dang K, Gong X, Zhao G, Wang H, Ivanistau A and Feng B (2020) Intercropping Alters the Soil Microbial Diversity and Community to Facilitate Nitrogen Assimilation: A Potential Mechanism for Increasing Proso Millet Grain Yield. Front. Microbiol. 11:601054. doi: 10.3389/fmicb.2020.601054
Received: 31 August 2020; Accepted: 02 November 2020;
Published: 24 November 2020.
Edited by:
Markus Puschenreiter, University of Natural Resources and Life Sciences Vienna, AustriaReviewed by:
Muhammad Umar Khan, Fujian Agriculture and Forestry University, ChinaLinkun Wu, Fujian Agriculture and Forestry University, China
Copyright © 2020 Dang, Gong, Zhao, Wang, Ivanistau and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Baili Feng, fengbaili@nwsuaf.edu.cn
†These authors have contributed equally to this work