Abstract
Modern agricultural practices have vastly increased crop production but negatively affected soil health. As such, there is a call to develop sustainable, ecologically-viable approaches to food production. Mixed-cropping of plant varieties can increase yields, although impacts on plant-associated microbial communities are unclear, despite their critical role in plant health and broader ecosystem function. We investigated how mixed-cropping between two field pea (Pisum sativum L.) varieties (Winfreda and Ambassador) influenced root-associated microbial communities and yield. The two varieties supported significantly different fungal and bacterial communities when grown as mono-crops. Mixed-cropping caused changes in microbial communities but with differences between varieties. Root bacterial communities of Winfreda remained stable in response to mixed-cropping, whereas those of Ambassador became more similar to Winfreda. Conversely, root fungal communities of Ambassador remained stable under mixed-cropping, and those of Winfreda shifted towards the composition of Ambassador. Microbial co-occurrence networks of both varieties were stronger and larger under mixed-cropping, which may improve stability and resilience in agricultural soils. Both varieties produced slightly higher yields under mixed-cropping, although overall Ambassador plants produced higher yields than Winfreda plants. Our results suggest that variety diversification may increase yield and promote microbial interactions.
Similar content being viewed by others
Introduction
Agriculture is a key area for the security of human health and well-being, and productivity may need to increase two-fold by 2050 to provide sufficient food for our rapidly growing human population1. Modern farming practices such as intensive tillage and large-scale input of inorganic fertilisers have helped drive high agricultural outputs, but at the cost of widespread environmental degradation and ecosystem instability2,3,4,5. Furthermore, historic and continued domestication and selection of crop species have driven a severe reduction in genetic diversity in modern cultivars6,7. This intensive approach to agriculture has negatively affected biodiversity and associated ecosystem service provision at multiple biological levels, resulting in reduced stability and resilience in agricultural systems8,9,10,11,12,13. There is now a concerted call to focus efforts on developing sustainable, integrated management processes that minimise intensive practices and instead focus on ecological solutions to increasing food production13,14,15,16.
Plant-associated microbial communities are critical to host function and fitness17,18. Root and rhizosphere communities are particularly vital, and the composition of these has significant implications for nutrient dynamics, immune function, pathogen susceptibility, and stress tolerance, among other traits19,20,21,22,23. For example, mycorrhizal fungi and nitrogen fixing bacteria are responsible for around 80% of nitrogen and 75% of phosphorus assimilation by plants in temperate forests24. The vast array of benefits conferred by plant-associated microbial communities has made harnessing the plant microbiome a key area of agricultural research in recent years, as its potential to revolutionise sustainable agriculture is recognised25,26,27. Practices such as the addition of organic amendments to soil, the use of cover crops, increasing crop diversity in fields, and reducing soil disturbance have all shown to be effective at improving microbial diversity in soil, potentially reducing the need for resource inputs whilst maintaining high yields13,14,16,28. As such, a deeper understanding of the ecological functions of the plant-associated microbial communities may help in the development of new farming practices with reduced ecological damage.
There is considerable variation in rhizosphere composition between plant species, and concurrently, mixed-cropping between crop species can cause significant changes in microbial community composition and diversity, with implications for productivity29,30,31,32. Greater plant diversity also provides increased ecosystem service provision and associated biodiversity33,34,35,36 and thus, such as an approach may have benefits at multiple levels. Furthermore, different varieties or cultivars of the same crop also support different root-associated microbial communities37,38,39,40,41. A meta-analysis of 91 studies indicated that growing different varieties of the same crop together can increase yield and stability, particularly under biotic and abiotic stressors such as disease or weather variability42. Intra-cropping different varieties can also increase below- and above-ground invertebrate diversity39 [but also see40]. As such, mixed-cropping between different varieties of the same crop can improve crop performance and ecosystem service provision16. However, the effects of combining different crop varieties on microbial community composition are poorly understood.
Here we use a glasshouse pot experiment to identify how mixed-cropping between two field pea (Pisum sativum L.) varieties influenced yield and root bacterial and fungal communities. We hypothesised that mixed-cropping will alter bacterial and fungal communities in comparison to mono-cropped varieties, with concurrent increases in yield, microbial diversity, and species interactions.
Methods
Glasshouse experiment and yield data
The study was conducted in a greenhouse using two field pea varieties; Ambassador (Tamar Organics, UK) and Winfreda (Heritage Seed Library, UK), which were grown in 3 L pots containing potting soil (Moorland Gold, UK). We collected three samples of potting soil at the start of the experiment and froze these at −80 °C for subsequent analysis of microbial communities. We sowed two seeds per pot, equidistant from one another and the sides of the pots. Each pot contained either two Ambassador seeds, two Winfreda seeds, or one of each, with 8 pots per mono-crop group and 11 pots of the mixed-crop. Pots were randomly distributed across the growing space in the glasshouse, and plants were watered daily with tap water. The study ran for six weeks, during which time all plants flowered and produced pods. We then removed all pods and weighed these collectively to gain yield data for each plant. These data were not normally distributed and so to test for differences in yield between the four treatment groups (mono-cropped Ambassador, mixed-cropped Ambassador, mixed-cropped Winfreda, and mono-cropped Winfreda), we used a generalised linear model with quasipoisson distribution in the car package43, with a Tukey’s post hoc analysis conducted using the multcomp package44 in RStudio45,46.
DNA extraction, 16S rRNA and ITS rRNA sequencing
At the end of the experiment, we isolated 250 mg of the root and associated soil close to the base of the above-ground part of the plant, and froze these samples at −80 °C until subsequent DNA extraction. We extracted DNA from root and soil samples using the Qiagen PowerSoil Pro kit (Qiagen, UK) according to the manufacturer’s instructions and quantified it using the high sensitivity assay kit on a QubitTM 3.0 Fluorometer (ThermoFisher Scientific, UK). We conducted 16S rRNA gene amplicon sequencing of the V4 region using dual indexed forward and reverse primers according to Griffiths & Harrison et al.47 and Kozich et al.48, and ITS rRNA gene sequencing according to Nguyen et al.49. For 16S rRNA amplification, we ran PCRs in duplicate using Solis BioDyne 5x HOT FIREPol® Blend Master Mix, 2 μM primers and 1 μl of sample DNA, with thermocycling conditions of: 95 °C for 15 minutes; 28 cycles of (95 °C for 20 s, 50 °C for 60 s, 72 °C for 60 s) and a final extension at 72 °C for 10 minutes. For ITS rRNA gene amplification, we used 4ul of DNA with thermocycling conditions of 95° for 10 minutes; 35 cycles of (95° for 30 s, 54° for 45 s, 72° for 60 sec) and a final extension at 72° for 10 minutes. We combined PCR replicates and cleaned these using HighPrepTM PCR clean up beads (MagBio, USA) according to the manufacturers’ instructions. We quality checked products using the Agilent 2200 TapeStation and quantified these on a QubitTM 3.0 Fluorometer. Samples were then pooled to equimolar concentrations. ITS and 16S rRNA gene amplicon sequencing was conducted on two separate runs using paired-end reads (2 × 250 bp) with v2 chemistry on the Illumina MiSeq platform at the University of Salford. Negative (extraction blanks) and positive (mock community) controls were included in the sequencing run.
Pre-processing of amplicon sequencing data
We conducted all analyses in RStudio (v1.0.153)46 for R (v3.4.1)45. We processed 16S rRNA and ITS rRNA gene amplicon sequences in DADA2 v1.5.050 using the default pipelines. A total of 4,057,971 raw sequence reads from 64 samples were generated from 16S rRNA gene sequencing, and 881,602 from ITS rRNA gene sequencing. Once paired end reads were merged, modal contig length was 253 bp (range of 240–260 bp) for 16S rRNA data and 219 bp (range of 210–434) for ITS rRNA data. For bacteria, we removed sequence variants (SVs) with length > 260 bp (46 out of 8516 SVs; 0.005% of total sequences) along with chimeras and two SVs found in the negative controls. As the ITS rRNA gene region is highly variable (Lindahl et al. 2013), we did not filter sequence variants based on length but we did remove chimeras and one SV identified in the negative controls. DADA2 identified all expected unique SVs in both the bacterial and fungal mock communities.
SVs with fewer than 100 reads across all samples were removed, leaving a median of 40870 SVs per sample (range of 14020–95216) for 16S rRNA data and 5684 SVs per sample (range of 2653–13954) for ITS rRNA data. We assigned 16S taxonomy using the SILVA v132 database51,52 and ITS taxonomy using the UNITE v7.2 database53 and exported the final SV tables, taxonomy tables and sample metadata to the phyloseq package54 for further analysis. To provide greater taxonomic detail about unidentified SVs and to stop the removal of these during analyses that agglomerate to a given taxonomic level, we fully annotated the taxonomy tables to species level using higher levels assignments (e.g. SV3 was named “Family_Prevotellaceae” at the genus and species levels).
Root community analyses
We constructed stacked plots to visualise the relative composition of root and soil bacterial and fungal communities at the class level according to treatment group (mono-cropped Ambassador, mixed-cropped Ambassador, mixed-cropped Winfreda, and mono-cropped Winfreda). We visualised the variation in the relative abundance of root communities of bacterial and fungal communities using NMDS plots (with Bray-Curtis distance) in phyloseq. We used a permutational ANOVA (PERMANOVA; adonis) in the vegan package55 to identify variation in the bacterial and fungal root communities according to treatment group. We determined the effect of variety, cropping type and their interaction on the relative abundance of the top ten bacterial and fungal families using one-way ANOVAs, and visualised this variation using box plots. We further visualised variation between treatment groups by plotting the relative abundance of genera as a heatmap (only bacterial genera with a combined relative abundance of > 1% across all four treatment groups were plotted, whereas all fungal genera were included). We calculated alpha-diversity (species richness) across the two microbial kingdoms by subsampling the raw SV count table to a standardised number of reads (equal to the sample with the lowest number of reads) using an iterative approach (100 times), averaging the diversity estimates from each trial. We tested for differences according to treatment group using a general linear model with quasipoisson distribution.
To identify relationships between root fungal and bacterial communities, we extracted Jensen-Shannon divergence and Jaccard distance matrices between all samples for both fungal and bacterial communities in the phyloseq.54 and vegan55 packages. We used mantel tests to correlate fungal and bacterial community distances for each measure separately. We rarefied the bacterial and fungal data sets to the same size as the sample with the lowest number of reads (13332 reads for bacteria and 1824 reads for fungi). We agglomerated these to genus level and then merged the phyloseq objects for bacterial and fungal communities for each plant. Using these cross-kingdom data, we calculated the co-occurrence between each pair of microbial genera by constructing a Spearman’s correlation coefficient matrix in the bioDist package56,57. We calculated the number of interactions with p < 0.05 for each treatment, in addition to those with −0.50 > rho > 0.50, and −0.75 > rho > 0.75 as indicators of network stability. We visualised those with rho > 0.75 (strong positive interactions) and rho < −0.75 (strong negative interactions) for the four treatment groups separately using network plots produced in igraph58.
Results
Yield data
There were significant differences in pod yield between the four treatment groups (X2 = 17.436, d.f. = 3, p < 0.001). Tukey’s post hoc analysis showed the yield for the mono-cropped Winfreda plants was significantly lower than that of mono-cropped Ambassador plants (p = 0.002) and mixed-cropped Ambassador plants (p = 0.009; Fig. 1). In general, yields were higher for the Ambassador variety than Winfreda, regardless of whether these were grown as mixed- or mono-crop, although yields were much more variable for Ambassador (Fig. 1). The average yield of each variety was higher when grown as a mixed-cropped than a mono-crop (Fig. 1), although these were not statistically significant (all p > 0.05).
Bacterial community composition
The taxonomic composition of bacterial communities in the initial experimental soil was different to those of pea plant roots, although it had similar classes represented but at different proportions (Fig. S1). Across both soil and roots, the main bacterial classes present were Actinobacteria, Alphaproteobacteria, Anaerolineae, Bacilli, Bacteroidia, Clostridia, Deltaproteobacteria, Gammaproteobacteria, and in plant roots only, Verrucomicrobiae (Fig. S1).
The taxonomic composition of the roots of pea plants were significantly different according to treatment group (F3,56 = 2.423, R2 = 0.078, p = 0.001; Fig. 2a). The root bacterial communities of the two varieties grown as mono-crops were most different to one another (Fig. 2a). The roots of the Winfreda variety had similar bacterial community compositions regardless of whether they were mono-cropped or mixed-cropped (Fig. 2a). The roots of the mixed-cropped Ambassador variety had a root bacterial community composition that was intermediate between that of mono-cropped Ambassador plants and Winfreda plants (Fig. 2a).
The ten most abundant bacterial families in the roots of the two pea varieties are shown in Fig. 3a. Variety had a significant effect on the relative abundance of five of the ten most abundant root bacterial families (Rhodanobacteraceae, Bacillaceae, Microscillaceae, Sphingomonadaceae, Chitinophagaceae), and the interaction between variety and cropping strategy had a significant effect on one (Rhizobiaceae) (all p < 0.05). In the case of Rhizobiaceae, mixed-cropping led to an increase in abundance of this family in the roots of Winfreda plants compared with its mono-cropped counterpart. Similar variation between treatment groups in the relative abundance of taxa is evident at the genus level (Fig. S2).
Fungal community composition
As with bacterial communities, the fungal community of soil at the start of the experiment was considerably different to that of pea plant root communities, with a relatively even distribution of fungal classes (Fig. S3). Four of the most dominant fungal classes in soil were also the most dominant in root communities; Agaricomycetes, Orbiliomycetes, Sordariomycetes and a group of unidentified fungi (Fig. S3). Treatment group had a significant effect on root fungal community composition (F3,56 = 1.618, R2 = 0.078, p = 0.024; Fig. 2b). The average root fungal communities of mono-cropped Winfreda pea plants were different to those of the other three treatment groups, which clustered together closely (Fig. 2b). Thus, converse to bacterial communities, root fungal communities of the Ambassador variety remained stable despite mixed-cropping, whereas those of the Winfreda variety plants became more similar to Ambassador under mixed-cropping (Fig. 2b).
The ten most abundant fungal families in the roots of the two pea varieties are shown in Fig. 3b. There was a significant effect of variety on the relative abundance of two of the ten most abundant root fungal families (Pleurotaceae and unidentified Coniochaetales), and cropping strategy had a significant effect on one (Agaricaceae) (all p < 0.05). Some amount of variation in the relative abundance of taxa between treatment groups is also evident at the genus level (Fig. S4). In contrast to the large number of bacterial genera identified through amplicon sequencing, relatively few fungal genera were identified (Fig. S4).
Cross-kingdom interactions
There were no significant differences in total microbial SV richness between treatment groups (F3,56 = 0.841, p = 0.477; Fig. 5). There was a significant relationship between fungal and bacterial communities of roots for both Jensen-Shannon divergence (JSD; r = 0.198, p = 0.003) and Jaccard’s distance measures (r = 0.144, p = 0.001). As fungal community distance increased, so did bacterial community distance, or in other words, roots with more similar fungal communities also had more similar bacterial communities (Fig. S6).
Co-occurrence analysis revealed that the roots of both pea varieties had a considerably greater number of statistically significant (p < 0.05) microbial interactions when mono-cropped than mixed-cropped (Table 1). Both varieties had larger networks (i.e. more microbial interactions) under a mono-cropping strategy than a mixed-cropping strategy, however the majority of these were relatively weak (i.e. −0.75 < rho < 0.75). When considering only strong (i.e. −0.75 > rho > 0.75) or strong positive (rho > 0.75) interactions, both varieties formed larger, stronger networks under a mixed-cropping strategy (Table 1). Approximately 70% of these interactions were positive (Table 1; Fig. 4a–d). The majority of microbial interactions were between bacteria, although fungal genera were also involved in both positive and negative interactions (Table 1; Fig. 4a–d).
Discussion
Here we demonstrate that mixed-cropping between two varieties of field peas alters root microbial community composition in comparison to cropping with single varieties (mono-cropping). The two varieties supported significantly different root bacterial and fungal communities to one another when mono-cropped. Under mixed-cropping, both fungal and bacterial communities tended to shift towards a composition that was intermediate to those of mono-cropped plants. Both bacterial and fungal growth and reproduction are dependent on nutrient type, quality and quantity59,60, and nutrient uptake and exudate deposition often varies between crop varieties61,62,63,64. Thus, changes in microbial communities may reflect differences in plant exudates produced under mixed-cropping in comparison to mono-cropping. There is growing evidence that plants, including crops, can recognise the relatedness of neighbours, possibly through root exudates or rhizosphere communities, and respond with changes in phenotype, resource use, exudate production, and gene expression65,66,67,68,69. Here, we provide novel evidence that plant root microbial communities also respond to the relatedness of their neighbours, although the mechanisms for this are not clear.
There were differences between the two varieties in the responses of bacterial and fungal communities to mixed-cropping. Under mixed-cropping, root bacterial communities of Ambassador shifted whilst those of Winfreda remained stable. Conversely, root fungal communities of Winfreda shifted under mixed-cropping, whilst those of Ambassador remained stable. That root fungal and bacterial communities of the two varieties responded differently to mixed-cropping suggests they are using plant- or microbe-derived carbon differently, or that the two kingdoms are undergoing niche partitioning. For example, there is evidence of variation in resource use for different taxonomic groups of bacteria and fungi, particularly in relation to nitrogen and carbon70,71,72,73. Although the two kingdoms do not necessarily exploit entirely different energy channels60, broadly speaking, bacteria are primarily considered r-strategists; fast-growing organisms that thrive in nutrient-rich, unstable environments with a reliance on labile carbon59,60,74. On the other hand, fungi tend to adopt both an r-strategy as well as a more a slow-growing K-strategy, making use of stable, low-nutrient environments and consuming both labile and recalcitrant carbon59,60,74. Thus, differences in microbial community structure are likely to reflect variation in exudate-derived carbon and/or nitrogen sources arising from the interaction between crop variety and cropping strategy. Despite varietal differences in the responses of bacterial and fungal communities to mixed-cropping, we found a broad, significant, correlation between the composition of the two microbial kingdoms in pea roots. Although these cross-kingdom interactions remain relatively under-studied, their role is of growing interest75. For example, cross-kingdom interactions may be important for biofilm production76,77, and fungal communities can influence bacterial colonisation through modulation of carbon, nitrogen, and pH78,79. Indeed, community composition may also be determined, in-part, through microbially-derived carbon production and other forms of interaction59,66,80,81. Soil pH also has a considerable effect on bacterial community composition, and a lesser effect on fungal communities82,83. Therefore, a greater understanding of how cropping strategies influence exudate production, nutrient cycling, and soil parameters may provide insight into the mechanisms driving community changes.
The co-occurrence results show that root communities form weaker networks under mono-cropping in comparison to mixed-cropping. This suggests plants provide a less stable environment for microbial community development when planted as a genetic monoculture. One might expect more unstable environments to arise from competition between crop varieties84, however, our results suggest mixed-cropping may promote larger, more stable microbial networks. This may have implications for a range of host traits and environmental aspects, such as yield output, pathogen resistance, carbon deposition, and soil fertility, to name a few16,60. However, in our study, mixed-cropping did not lead to increases in microbial diversity, which would likely be of benefit to agricultural soils and associated biodiversity5,16. Our results also show that in some cases, mixed-cropping introduced microbes to the roots of the other variety, and in other cases, some microbes appear in mixed-cropped plants without necessarily being introduced by the other variety. This further supports the hypothesis that mixed-cropping leads to novel exudates that alter microbial community composition. Further work is required to understand the function of the different members of the root microbial community, and how different varieties can be combined to facilitate the development of communities that maximise improvements to food production and carbon storage, among other desirable attributes for agricultural systems16.
Root microbial communities at the end of the study were considerably different to the microbial communities associated with the soil in which they were grown, albeit with similar groups, indicating that pea plants select their microbial communities from their environment. This is a well-known phenomenon in plants known as the ‘rhizosphere effect’ (reviewed in74). Previous studies have also shown that the microbial communities that inhabit the rhizosphere are significantly different in the early development stages compared with the latter vegetative, flowering, bolting, and senescence stages, during which times, microbial communities tend to converge85,86,87,88. This is likely linked to changes in root secretions, which increase throughout development with higher root secretion of sugars in early growth stages, which are then replaced by phenolics and specific amino acids in late growth stages, with a reduction in secretions post-flowering64,88,89,90. We collected root samples after flowering so that we could also obtain yield data, which suggests we sampled the end-state microbial communities rather the earlier, dynamic communities, which may have shown more pronounced differences according to variety and cropping strategy.
We found much larger bacterial than fungal communities associated with the roots of field peas. As a legume and a member of the Fabaceae family, peas have a strong association with rhizobia bacteria, which play an important role in nitrogen fixation, particularly in nutrient-poor soil19. Furthermore, bacteria are expected to dominate over fungi in high-nutrient environments59,91,92. Thus, the use of high-quality compost rather than nutrient-poor soil in this study is likely to have affected the community composition and may further explain the relatively low fungal diversity identified in our study. In addition, soil composition and structure also influence root architecture, exudate production and plant signalling, with implications for microbial communities19,93,94,95. Therefore, the response of microbial communities of peas (and other crops) to mixed-cropping may be even more pronounced in nutrient-depleted or highly-disturbed soils, as are commonly associated with agricultural fields. Indeed, in such soils, the use of different crop genotypes may be a critical tool for maximising the availability of nutrients for crop growth95. The effects of mixed-cropping on nutrient quality of the crop would also be of considerable interest for human health96,97.
Although Ambassador plants had significantly higher yields than Winfreda, we found that mixed-cropping increased the average yield of both varieties relative to its mono-cropped counterpart. Plants may be benefitting from mixed-cropping via release from negative plant-soil feedbacks conferred by mono-cropped rhizospheres66,98,99, or through altered carbon translocation100. Thus, mixed-cropping between varieties may have the potential to increase yields through selection of complimentary varieties that promote positive plant-soil feedback relationships or carbon assimilation. A meta-analysis of 91 intercropping studies found a negative correlation between nutrient content of soil and the net yield benefit from the use of intra-cropping42, compared with soils high in nutrients, which resulted in low to no yield increases. This suggests that when nutrient content of the soil is low, a higher yield benefit from intercropping is likely. This may be attributed to a greater level of organic matter mineralisation by soil microbes when plants are mixed-cropped, increasing nutrient availably to plants42,101. Thus, greater improvements to yields arising from mixed-cropping strategies may also be more evident in nutrient-poor soils. Further trials with more genotypes, particularly under field settings, may identify complimentary varieties that maximise yields and restore microbial diversity to agricultural landscapes16. Given that there is a heritable basis to root exudates93 and microbial communities38,102, such studies would ideally be conducted in combination with host genotype analyses across the spectrum of host diversity. Furthermore, the implications of such cropping strategies for associated biodiversity and ecosystem service provision are a particularly important aspect to consider.
In conclusion, we show that mixed-cropping between two field pea varieties leads to differential changes in bacterial and fungal community composition of roots. This may be related to changes in soil chemistry or microbe-microbe interactions, or result from differences in plant exudate production arising from inter-variety competition. Mixed-cropped plants supported larger, more stable microbial co-occurrence networks, suggesting this approach may improve resilience and stability of root-associated communities and subsequently, the crop. More work is required to identify mechanistic drivers of variation in microbial communities arising from mixed-cropping, in addition to combinations of genotypes that can maximise yield production, disease resistance, and wider ecosystem provision in agricultural settings. Such an approach may also provide an opportunity to restore microbial diversity to degraded agricultural ecosystems103,104.
Data availability
All sequencing data produced during this study are available on the NCBI SRA database under project numbers PRJNA505960 and PRJNA526546. All analysis code is available on request from Rachael Antwis.
References
Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl. Acad. Sci. 108, 20260–20264 (2011).
Godfray, H. C. J. et al. Food Security: The Challenge of Feeding 9 Billion People. Science (80-.). 327, 812–819 (2010).
Nkonya, E., Mirzabaev, A. & von Braun, J. Economics of Land Degradation and Improvement – A Global Assessment for Sustainable Development (2016).
Chen, Y. H., Gols, R. & Benrey, B. Crop Domestication and Its Impact on Naturally Selected Trophic Interactions. Annu. Rev. Entomol. 60, 35–58 (2014).
Tsiafouli, M. A. et al. Intensive agriculture reduces soil biodiversity across Europe. Glob. Chang. Biol. 21, 973–985 (2014).
Doebley, J. F., Gaut, B. S. & Smith, B. D. Review The Molecular Genetics of Crop Domestication. 1309–1321 (2006).
Shi, J. & Lai, J. Patterns of genomic changes with crop domestication and breeding. Curr. Opin. Plant Biol. 24, 47–53 (2015).
Flohre, A., Rudnick, M., Traser, G., Tscharntke, T. & Eggers, T. Does soil biota benefit from organic farming in complex vs. simple landscapes? Agric. Ecosyst. Environ. 141, 210–214 (2011).
Chaudhry, V., Rehman, A., Mishra, A., Chauhan, P. S. & Nautiyal, C. S. Changes in Bacterial Community Structure of Agricultural Land Due to Long-Term Organic and Chemical Amendments. Microb. Ecol. 64, 450–460 (2012).
Hartmann, M., Frey, B., Mayer, J., Mäder, P. & Widmer, F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 9, 1177–1194 (2015).
Leff, J. W. et al. Predicting the structure of soil communities from plant community taxonomy, phylogeny, and traits. ISME J. 12, 1794–1805 (2018).
Pershina, E. et al. Comparative analysis of prokaryotic communities associated with organic and conventional farming systems. PLoS One 10, 1–16 (2015).
Bommarco, R., Kleijn, D. & Potts, S. G. Ecological intensification: Harnessing ecosystem services for food security. Trends Ecol. Evol. 28, 230–238 (2013).
Godfray, H. C. J. & Garnett, T. Food security and the role of sustainable intensification. Ecos 6–11, https://doi.org/10.1071/ec14167 (2014).
Gurr, G. M. et al. Multi-country evidence that crop diversification promotes ecological intensification of agriculture. Nat. Plants 2, 22–25 (2016).
Barot, S. et al. Designing mixtures of varieties for multifunctional agriculture with the help of ecology. A review. Agron. Sustain. Dev. 37, (2017).
Turner. The Plant Microbiome. Adv. Bot. Res. 69, 279–309 (2014).
Lareen, A., Burton, F. & Schäfer, P. Plant root-microbe communication in shaping root microbiomes. Plant Mol. Biol. 90, 575–587 (2016).
Franche, C., Lindström, K. & Elmerich, C. Nitrogen-fixing bacteria associated with leguminous and non-leguminous plants. Plant Soil 321, 35–59 (2009).
Marcel, B. Functional biology of plant phosphate uptake at root and mycorrhiza interfaces. New Phytol. 173, 11–26 (2007).
Doornbos, R. F., Van Loon, L. C. & Bakker, P. A. H. M. Impact of root exudates and plant defense signaling on bacterial communities in the rhizosphere. A review. Agron. Sustain. Dev. 32, 227–243 (2012).
Timmusk, S. et al. Drought-tolerance of wheat improved by rhizosphere bacteria from harsh environments: Enhanced biomass production and reduced emissions of stress volatiles. PLoS One 9 (2014).
Thomas, R. Turner et al. The Plant Microbiome. Genome Biol. 14, 1–10 (2013).
Van Der Heijden, M. G. A., Bardgett, R. D. & Van Straalen, N. M. The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310 (2008).
Bender, S. F., Wagg, C. & van der Heijden, M. G. A. An Underground Revolution: Biodiversity and Soil Ecological Engineering for Agricultural Sustainability. Trends Ecol. Evol. 31, 440–452 (2016).
Quiza, L., St-arnaud, M., Yergeau, E. & Rey, T. R. Harnessing phytomicrobiome signaling for rhizosphere microbiome engineering. 6, 1–11 (2015).
Toju, H. et al. Core microbiomes for sustainable agroecosystems. Nat. Plants 4, 247–257 (2018).
Wezel, A. et al. Agroecological practices for sustainable agriculture. A review. Agron. Sustain. Dev. 34, 1–20 (2014).
Granzow, S. et al. The Effects of Cropping Regimes on Fungal and Bacterial Communities of Wheat and Faba Bean in a Greenhouse Pot Experiment Differ between Plant Species and Compartment. Front. Microbiol. 8, 1–22 (2017).
Lian, T. et al. Impact of intercropping on the coupling between soil microbial community structure, activity, and nutrient-use efficiencies. PeerJ 7, e6412 (2019).
Li, S. & Wu, F. Diversity and co-occurrence patterns of soil bacterial and fungal communities in seven intercropping systems. Front. Microbiol. 9, 1–13 (2018).
Yang, Z. et al. Variation of bacterial community diversity in rhizosphere soil of sole-cropped versus intercropped Wheat field after harvest. PLoS One 11, 1–18 (2016).
Thakur, M. P. et al. Plant diversity drives soil microbial biomass carbon in grasslands irrespective of global environmental change factors. Glob. Chang. Biol. 21, 4076–4085 (2015).
Lange, M. et al. Plant diversity increases soil microbial activity and soil carbon storage. Nat. Commun. 6 (2015).
De Deyn, G. et al. Additional carbon sequestration benefits of grassland diversity restoration. J. Appl. Ecol. 48, 600–608 (2011).
Scherber, C. et al. Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. Nature 468, 553–556 (2010).
Inceoǧlu, Ö., Salles, J. F., Van Overbeek, L. & Van Elsas, J. D. Effects of plant genotype and growth stage on the betaproteobacterial communities associated with different potato cultivars in two fields. Appl. Environ. Microbiol. 76, 3675–3684 (2010).
Peiffer, J. A. et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. 110, 6548–6553 (2013).
Aira, M., Gómez-Brandón, M., Lazcano, C., Bååth, E. & Domínguez, J. Plant genotype strongly modifies the structure and growth of maize rhizosphere microbial communities. Soil Biol. Biochem. 42, 2276–2281 (2010).
Mahoney, A. K., Yin, C. & Hulbert, S. H. Community Structure, Species Variation, and Potential Functions of Rhizosphere-Associated Bacteria of Different Winter Wheat (Triticum aestivum) Cultivars. Front. Plant Sci. 8, 1–14 (2017).
Bouffaud, M.-L., Poirier, M.-A., Muller, D. & Moënne-Loccoz, Y. Root microbiome relates to plant host evolution in maize and other Poaceae. Environ. Microbiol. 16, 2804–2814 (2014).
Reiss, E. R. & Drinkwater, L. E. Cultivar mixtures: A meta-analysis of the effect of intraspecific diversity on crop yield: A. Ecol. Appl. 28, 62–77 (2018).
Fox, J. & Weisberg, S. An R Companion to Applied Regression, Second Edition (Thousand Oaks CA: Sage. 2011).
Hothorn, T., Bretz, F., Westfall, P., Heiberger, R. M. & Schuetzenmeister, A. multcomp: Simultaneous Inference in General Parametric Models. 50, 346–363 (2013).
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, https://www.R-project.org/ (2017).
RStudio Team. RStudio: Integrated Development for R. RStudio, Inc., Boston, MA, http://www.rstudio.c (2016).
Griffiths, S. M. et al. Genetic variability and ontogeny predict microbiome structure in a disease-challenged montane amphibian. ISME J. 1, https://doi.org/10.1038/s41396-018-0167-0 (2018).
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120 (2013).
Nguyen, N. H., Smith, D., Peay, K. & Kennedy, P. Parsing ecological signal from noise in next generation amplicon sequencing (2014).
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).
Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596 (2013).
Yilmaz, P. et al. The SILVA and ‘all-species Living Tree Project (LTP)’ taxonomic frameworks. Nucleic Acids Res. 42, 643–648 (2014).
UNITE. UNITE general FASTA release. Version 01.12.2017 (2017).
McMurdie, P. J. & Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS One 8 (2013).
Oksanen, J. et al. vegan: Community Ecology Package (2018).
Williams, R. J., Howe, A. & Hofmockel, K. S. Demonstrating microbial co-occurrence pattern analyses within and between ecosystems. Front. Microbiol. 5, 1–10 (2014).
Ding, B., Gentleman, R. & Carey, V. bioDist: Different distance measures. R package version 1.54.0 (2018).
Csárdi, G. & Nepusz, T. The igraph software package for complex network research. InterJournal Complex Syst. 1695, 1695 (2006).
Rousk, J. & Frey, S. D. Revisiting the hypothesis that fungal-to-bacterial dominance characterizes turnover of soil organic matter and nutrients. Ecol. Monogr. 85, 457–472 (2015).
de Vries, F. T. & Caruso, T. Eating from the same plate? Revisiting the role of labile carbon inputs in the soil food web. Soil Biol. Biochem. 102, 4–9 (2016).
Barraclough, P. B. et al. Nitrogen efficiency of wheat: Genotypic and environmental variation and prospects for improvement. Eur. J. Agron. 33, 1–11 (2010).
Schmidt, J. E., Bowles, T. M. & Gaudin, A. C. M. Using Ancient Traits to Convert Soil Health into Crop Yield: Impact of Selection on Maize Root and Rhizosphere Function. Front. Plant Sci. 7, 1–11 (2016).
Pérez-Jaramillo, J. E., Mendes, R. & Raaijmakers, J. M. Impact of plant domestication on rhizosphere microbiome assembly and functions. Plant Mol. Biol. 90, 635–644 (2016).
Aulakh, M. S., Wassmann, R., Bueno, C., Kreuzwieser, J. & Rennenberg, H. Characterization of Root Exudates at Different Growth Stages of Ten Rice (Oryza sativa L.) Cultivars. Plant Biol. 3, 298–298 (2002).
Semchenko, M., Saar, S. & Lepik, A. Plant root exudates mediate neighbour recognition and trigger complex behavioural changes. New Phytol. 204, 631–637 (2014).
Hodge, A. & Fitter, A. H. Microbial mediation of plant competition and community structure. Funct. Ecol. 27, 865–875 (2013).
Murphy, G. P., Swanton, C. J., Van Acker, R. C. & Dudley, S. A. Kin recognition, multilevel selection and altruism in crop sustainability. J. Ecol. 105, 930–934 (2017).
Schofield, E. J. et al. Cultivar Differences and Impact of Plant-Plant Competition on Temporal Patterns of Nitrogen and Biomass Accumulation. 10, 1–9 (2019).
Biedrzycki, M. L., Venkatachalam, L. & Bais, H. P. Transcriptome analysis of Arabidopsis thaliana plants in response to kin and stranger recognition. Plant Signal. Behav. 6, 1515–1524 (2011).
Kramer, S. et al. Resource partitioning between bacteria, fungi, and protists in the detritusphere of an agricultural soil. Front. Microbiol. 7, 1–12 (2016).
Lennon, J. T. Mapping the Niche Space of Soil Microorganisms. Bull. Ecol. Soc. Am. 93, 351–351 (2012).
Baran, R. et al. Exometabolite niche partitioning among sympatric soil bacteria. Nat. Commun. 6 (2015).
Smith, M. E. et al. Investigating niche partitioning of ectomycorrhizal fungi in specialized rooting zones of the monodominant leguminous tree Dicymbe corymbosa. New Phytol. 215, 443–453 (2017).
Philippot, L., Raaijmakers, J. M., Lemanceau, P. & Van Der Putten, W. H. Going back to the roots: The microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799 (2013).
Menezes, A. B. D., Richardson, A. E. & Thrall, P. H. Linking fungal – bacterial co-occurrences to soil ecosystem function ScienceDirect Linking fungal – bacterial co-occurrences to soil ecosystem function. Curr. Opin. Microbiol. 37, 135–141 (2017).
Frey-Klett, P. et al. Bacterial-Fungal Interactions: Hyphens between Agricultural, Clinical, Environmental, and Food Microbiologists. Microbiol. Mol. Biol. Rev. 75, 583–609 (2011).
van Overbeek, L. S. & Saikkonen, K. Impact of Bacterial-Fungal Interactions on the Colonization of the Endosphere. Trends Plant Sci. 21, 230–242 (2016).
Johnston, S. R., Hiscox, J., Savoury, M., Boddy, L. & Weightman, A. J. Highly competitive fungi manipulate bacterial communities in decomposing beech wood (Fagus sylvativa). FEMS Microbiol. Ecol., https://doi.org/10.1093/femsec/fiy225 (2018).
Hassani, M. A., Durán, P. & Hacquard, S. Microbial interactions within the plant holobiont. Microbiome 6, 58 (2018).
Sasse, J., Martinoia, E. & Northen, T. Feed Your Friends: Do Plant Exudates Shape the Root Microbiome? Trends Plant Sci. 23, 25–41 (2018).
Venturi, V. & Keel, C. Signaling in the Rhizosphere. Trends Plant Sci. 21, 187–198 (2016).
Fierer, N. Embracing the unknown: Disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590 (2017).
Rousk, J. et al. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351 (2010).
Simonsen, A. K., Chow, T. & Stinchcombe, J. R. Reduced plant competition among kin can be explained by Jensen’s inequality. Ecol. Evol. 4, 4454–4466 (2014).
Chaparro, J. M., Badri, D. V. & Vivanco, J. M. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 8, 790–803 (2014).
Pfeiffer, S. et al. Rhizosphere microbiomes of potato cultivated in the High Andes show stable and dynamic core microbiomes with different responses to plant development. FEMS Microbiol. Ecol. 93, fiw242 (2017).
Micallef, S. A., Channer, S., Shiaris, M. P. & Colón-Carmona, A. Plant age and genotype impact the progression of bacterial community succession in the Arabidopsis rhizosphere. Plant Signaling and Behavior 4, 777–780 (2009).
Zhalnina, K. et al. Rhizosphere Microbial Community Assembly. Nat. Microbiol. 3, 470–480 (2018).
Chaparro, J. M. et al. Root Exudation of Phytochemicals in Arabidopsis Follows Specific Patterns That Are Developmentally Programmed and Correlate with Soil Microbial Functions. PLoS One 8, e55731 (2013).
Badri, D. V. & Vivanco, J. M. Regulation and function of root exudates. Plant, Cell Environ. 32, 666–681 (2009).
Strickland, M. S. & Rousk, J. Considering fungal: Bacterial dominance in soils - Methods, controls, and ecosystem implications. Soil Biol. Biochem. 42, 1385–1395 (2010).
Fierer, N., Strickland, M. S., Liptzin, D., Bradford, M. A. & Cleveland, C. C. Global patterns in belowground communities. Ecol. Lett. 12, 1238–1249 (2009).
Iannucci, A., Fragasso, M., Beleggia, R., Nigro, F. & Papa, R. Evolution of the Crop Rhizosphere: Impact of Domestication on Root Exudates in Tetraploid Wheat (Triticum turgidum L.). Front. Plant Sci. 8 (2017).
Neumann, G. et al. Root exudation and root development of lettuce (Lactuca sativa L. cv. Tizian) as affected by different soils. Front. Microbiol. 5, 1–6 (2014).
Rengel, Z. & Marschner, P. Nutrient availability and management in the rhizosphere: Exploiting genotypic differences. New Phytol. 168, 305–312 (2005).
Cooper, R. Re-discovering ancient wheat varieties as functional foods. 5, 138–143 (2015).
Goicoechea, N. & Antolín, M. C. Increased nutritional value in food crops. Microb. Biotechnol. 10, 1004–1007 (2017).
Kulmatiski, A., Beard, K. H., Stevens, J. R. & Cobbold, S. M. Plant-soil feedbacks: A meta-analytical review. Ecol. Lett. 11, 980–992 (2008).
Eck, J. L., Stump, S. M., Delavaux, C. S., Mangan, S. A. & Comita, L. S. Evidence of within-species specialization by soil microbes and the implications for plant community diversity, https://doi.org/10.1073/pnas.1810767116 (2019).
Deyn, G. B. D., Quirk, H., Oakley, S., Ostle, N. J. & Bardgett, R. D. Increased Plant Carbon Translocation Linked to Overyielding in Grassland Species Mixtures. PLoS One 7, e45926 (2012).
Hinsinger, P. et al. P for Two, Sharing a Scarce Resource: Soil Phosphorus Acquisition in the Rhizosphere of Intercropped Species. Plant Physiol. 156, 1078–1086 (2011).
Schweitzer, J. A. et al. Plant-Soil-Microorganism Interactions: Heritable Relationship between Plant Genotype and Associated Soil Microorganisms. Ecology 89, 773–781 (2008).
Chateil, C. et al. Crop genetic diversity benefits farmland biodiversity in cultivated fields. Agric. Ecosyst. Environ. 171, 7 (2013).
Dubs, F. et al. Positive effects of wheat variety mixtures on aboveground arthropods are weak and variable. 33, 66–78 (2018).
Acknowledgements
We are sincerely grateful to the Heritage Seed Library (part of the Henry Doubleday Research Association) for providing the Winfreda seeds used in this study. Our thanks also to BEI Resources, the US Forest Service, and the University of Wisconsin-Madison for provision of bacterial and fungal mock communities for amplicon sequencing.
Author information
Authors and Affiliations
Contributions
R.E.A. conceived the idea and designed the methodology; R.E.A. and A.H. collected the data; R.E.A. and S.S.B. analysed the data; R.E.A. and A.H. wrote the manuscript. All authors contributed to the drafts and gave final approval for publication.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Horner, A., Browett, S.S. & Antwis, R.E. Mixed-Cropping Between Field Pea Varieties Alters Root Bacterial and Fungal Communities. Sci Rep 9, 16953 (2019). https://doi.org/10.1038/s41598-019-53342-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-53342-8