Cropping practices manipulate abundance patterns of root and soil microbiome members paving the way to smart farming
© The Author(s). 2018
Received: 2 August 2017
Accepted: 17 December 2017
Published: 16 January 2018
Harnessing beneficial microbes presents a promising strategy to optimize plant growth and agricultural sustainability. Little is known to which extent and how specifically soil and plant microbiomes can be manipulated through different cropping practices. Here, we investigated soil and wheat root microbial communities in a cropping system experiment consisting of conventional and organic managements, both with different tillage intensities.
While microbial richness was marginally affected, we found pronounced cropping effects on community composition, which were specific for the respective microbiomes. Soil bacterial communities were primarily structured by tillage, whereas soil fungal communities responded mainly to management type with additional effects by tillage. In roots, management type was also the driving factor for bacteria but not for fungi, which were generally determined by changes in tillage intensity. To quantify an “effect size” for microbiota manipulation, we found that about 10% of variation in microbial communities was explained by the tested cropping practices. Cropping sensitive microbes were taxonomically diverse, and they responded in guilds of taxa to the specific practices. These microbes also included frequent community members or members co-occurring with many other microbes in the community, suggesting that cropping practices may allow manipulation of influential community members.
Understanding the abundance patterns of cropping sensitive microbes presents the basis towards developing microbiota management strategies for smart farming. For future targeted microbiota management—e.g., to foster certain microbes with specific agricultural practices—a next step will be to identify the functional traits of the cropping sensitive microbes.
Agricultural intensification has resulted in an increased production of staple crops such as wheat, rice, and maize and lead to greater food security for a continuously growing world population [1, 2]. Despite these benefits, there is increasing awareness about the adverse environmental impacts arising from the intensive practices of modern agriculture. These include increased greenhouse gas emissions and nutrient leaching as a result of intensive fertilizer application , increased soil erosion , and detrimental effects on biodiversity [5, 6]. To alleviate such deleterious effects, an ecological intensification has been proposed that focuses on meeting standards of environmental quality while promoting and maintaining organisms that provide beneficial ecosystem services [7, 8]. A number of practices improve the sustainability of agriculture, including organic farming  and reduced or no-tillage . These practices aim to enhance soil fertility while maintaining crop yields through supporting a diverse and active soil biota . Soil biota includes microbes such as bacteria and fungi that collectively function as a microbiome. Bacteria and fungi regulate many ecosystem processes and play key roles in nutrient cycling through decomposition of organic matter, and transformation and fixation of important soil nutrients like nitrogen and phosphorus .
Aside from the environmental benefits of organic agriculture  and less intensive tillage regimes , there is still debate about the effects of these cropping practices on belowground microbial communities. In general, arable management affects community composition and diversity; although such effects may depend on the microbial kingdom being studied and the different farming systems being compared [14, 15]. However, there are few agricultural experiments comparing conventional and organic farming practices  and fewer that compare different management types and tillage intensities . Therefore, an agricultural experiment combining these two aspects at a single site allows to separate the effects of management type and tillage on microbial communities and minimize variation caused by soil spatial heterogeneity. The Farming System and Tillage experiment (FAST) was established in 2009 near Zürich to address this for the main arable cropping systems in Switzerland (Additional file 1: Fig. S1). These cropping systems are, namely, conventional (C) and organic (O) management types, with different tillage intensities (no-tillage (NT), reduced-tillage (RT), and intensive tillage (IT)). The FAST experiment compares the 4 main cropping systems C-IT, C-NT, O-IT, O-RT. Wittwer et al.  described the design of the FAST experiment and provided agronomic insights into the effects of different farming practices on winter wheat and maize yields. However, the effects of management types and different tillage intensities on belowground microbial communities at the site have not been explored.
Soil microbial communities and their associated functions largely determine the productivity of agroecosystems . The composition of the soil microbiome presents the major driver in shaping the bacterial and fungal communities associated with plant roots [20, 21]. Root microbiome is an important determinant for plant growth and health by assisting in nutrient uptake, supporting abiotic stress tolerance and protecting the host from pathogens [20, 22]. A number of recent microbial community surveys have described the root microbiomes of Arabidopsis thaliana [23, 24], clover , maize , rice , sugarcane , and grapevine  and reported significant effects of soil type on root microbiome composition. If soil and root microbial communities are closely linked, root microbial communities may also be affected by agronomic practices . To date, the effects of agricultural practices on root microbial communities remain still poorly understood, owing to contrasting reports and the use of low-resolution fingerprinting methods . Using high-throughput sequencing, we aimed to unravel how root microbial communities respond to conventional and organic agriculture and various tillage regimes.
Members of the soil and root microbiome interact directly and indirectly with each other, and a tool for better understanding of these potential interactions is co-occurrence network analysis [32, 33]. Long used in the social sciences to analyze relationships between humans , network analyses have recently been applied in soil microbial ecology to explore patterns of community assembly , visualize response patterns of different taxonomic groups to agronomic practices , and to identify individual microbiome members that significantly influence community composition . It was recently shown that soils under conventional and organic management harbor distinct microbial networks in each farming system . To date, the effects of different cropping practices on co-occurrence patterns in the root microbiome remain unexplored.
From the perspective of microbiome management, it is important to understand which microbes are sensitive to cropping practices and whether they possess specific network properties. Microbes that frequently co-occur with many others are referred to as keystone taxa because they may play an ecologically important role by determining community dynamics and microbiome functioning [37–39]. It is unclear whether keystone taxa in soil and root microbiomes are responsive to cropping practices. More importantly, are cropping sensitive microbes solitary community members, or do they belong to guilds of simultaneously responding taxa? Are they frequent or not? Such information is relevant for implementing agricultural management strategies to promote specific microbes that contribute to soil fertility and plant health.
With these ideas in mind, we investigated the impact of cropping practices at the FAST experimental site on soil and root bacterial and fungal communities in winter wheat using amplicon sequencing and network analysis. We specifically asked: (1) Do soil and root microbial communities differ in their responses to management type and tillage intensities? (2) Which microbes are the indicator taxa for particular cropping practices (conventional vs. organic; reduced vs. intensive tillage)? (3) How do cropping practices impact co-occurrence patterns of soil- and root-associated microbes? (4) What are the network characteristics (abundance, degree of co-occurrence, and keystoneness) of cropping sensitive microbes?
Soil and root microbiota
We conducted separate bacterial and fungal community profiling of 16 soil and 16 winter wheat root samples from of the FAST experiment (Additional file 1: Fig. S1) to investigate the effects of management type and tillage intensity on microbial communities. The bacterial community profiling yielded a total of 639,440 high-quality sequences (range 11,192–37,255; median 18,122; Additional file 2). Fungal profiling yielded 962,619 sequences, ranging between 9138 and 48,750 sequences per sample (median 30,284). We identified 2972 bacterial, 3 archaeal, and 1975 fungal operational taxonomic units (OTUs) across all samples (Additional file 1: Fig. S2).
For α-diversity analyses, we rarified the communities to 11,000 (bacteria) and 9000 (fungi) sequences per sample, which captured most of the observed OTU richness (Additional file 1: Fig. S5). Soils supported higher species richness than roots with bacterial communities being greater in richness than fungi (Additional file 1: Fig. S5 and Table S3). In both soil and root communities, bacteria and fungi richness was highest in O-IT samples with significant effect for bacterial communities in root samples. To conclude, plant root and soil microbiota differ markedly in richness, composition, and taxonomy.
Cropping system effects on soil and root microbial communities
Different patterns were observed for the root microbiota. Root bacteria formed four distinct clusters in the ordination with axis 1 again separated the samples by management type and axis 2 separated the samples by tillage intensity. Pairwise PERMANOVA comparisons detected significant differences between the two conventional treatments and O-RT but not O-IT samples. For root fungal communities, CAP separated the O-RT samples along axis 1 and the C-NT samples from the other treatments on axis 2. PERMANOVA also confirmed a general effect of cropping system, but no pairwise differences on community dissimilarity were found.
Since β-diversity can be driven by true biological differences, differences in group dispersion (variance), or both , we tested for differences in dispersion for both soil and root microbiota using BETADISP. The lack of significance in these dispersion tests suggested that differences between cropping systems were driven primarily by true biological differences and not an artifact of differences of within-group dispersion (Additional file 1: Table S4). In summary, while tillage-driven differences were seen in the soil bacterial community, the management type appeared to be the main driving factor in root bacteria. Conversely, root fungal communities did not strongly respond to management type induced changes in soil and instead were determined by changes in the tillage intensity.
Identifying cropping sensitive OTUs
As indicator OTUs were solely identified based on correlation, we validated them using likelihood ratio tests implemented in edgeR (; Additional file 3). Finally, we defined the OTUs that were supported by both methods as cropping sensitive OTUs (hereafter: csOTUs). In soil, we found a total of 53 and 70 bacteria and fungi csOTUs, respectively (Additional file 1: Fig. S6). As approximation for an “effect size” of cropping practices on microbial communities, we calculated these bacteria and fungi csOTUs to account for 8.3 and 9.9% of the total soil community sequences. Similarly, we identified 62 and 36 csOTUs for root bacteria and fungi, corresponding to an effects size of 14.2 and 5.0%, respectively. Consistent with the previous conclusion that cropping practices affected soil and root communities differently, we saw little overlap between bacteria and fungi csOTUs comparing root and soil samples. While the identified csOTUs responded to specific cropping systems, they did not exhibit a particular taxonomic pattern with cropping system (Additional file 1: Figs. S7-S12 and Supplementary Results). Taken together, each cropping system supports a specialized subset of soil and root bacteria and fungi, while the majority of the communities are shared between management types and tillage regimes.
Cropping effects on microbial co-occurrence patterns
Lastly, we explored the extent to which management types and tillage regimes impacted co-occurrence patterns in microbial communities. We first constructed separate co-occurrence networks for soil and root bacterial and fungal communities and determined their properties (see the “Methods” section). Consistent with the α-diversity analyses (Additional file 1: Fig. S5), the soil bacteria network comprised the highest number of significantly co-occurring OTUs, followed by intermediate and similar numbers in the soil fungi and root bacteria networks (Additional file 1: Fig. S13). Consequently, network connectivity (measured by average number of connections per OTU) was higher in the soil bacteria and soil fungi networks than the root bacteria network. The root fungi network comprised the fewest OTUs and was the least complex. We also mapped the csOTUs (as defined in Additional file 1: Fig. S6) into the microbial networks, and we found them agglomerating according to management type and/or tillage intensity (Additional file 1: Fig. S13).
Properties of soil and root meta co-occurrence networks
Taken together, we show that the differing cropping practices markedly alter co-occurrence patterns of numerous bacteria and fungi, and these impacts affected microbiome members largely independent of their abundance and connectivity.
While considering the effects of agricultural practices on soil microbial communities, previous studies have often been limited to the examination of single factors like management type [42–44], tillage [45–47], soil amendments [48–50], or focused on either the bacterial or fungal kingdom. Hartmann et al.  investigated effects of management type on soil bacteria and fungi in a multi-factor long-term agricultural experiment. Our motivation was to expand on these efforts by examining cropping system effects on root microbes and by also including the effect of tillage into the experimental design. To our knowledge, this is the first study to investigate how conventional and organic managements, both with different tillage regimes, influence bacterial and fungal communities in both soil and root compartments. The strength of the FAST experiment is that management and tillage effects can be studied factorially and independent of possibly confounding variables such as soil types. However, a broader generalization of the presented findings requires further studies with greater sample numbers, from multiple sites, across different climatic conditions and accounting for temporal and special variation in root and soil microbiomes.
Differential responses of soil and root microbiomes
The specific sets of microbes in roots and soil (Fig. 1) explain the clear separation by compartment in the ordination analyses (Additional file 1: Fig. S4). This distinctiveness of the plant root microbiome was also found in previous studies of other plant species [23–26, 28, 29, 51]. Plants recruit a root microbiome in their early life stages from a larger pool of soil microbes [27, 52]. The initial composition of this soil microbial pool is the most influential factor determining the composition of root microbial communities [23, 24, 53, 54]. Therefore, we also investigated if root bacterial and fungal communities reflected cropping system-driven differences in soil microbial communities (see Additional file 1: Supplementary Discussion).
We observed compartment-specific responses of bacteria and fungi to the tested cropping practices (Fig. 2, Additional file 1: Table S4). Dissimilarities in soil bacteria were generally driven by differences in tillage regimes, whereas management type, together with tillage, was influential for the soil fungi. Notably, the most influential cropping factors driving differences in soil communities were not necessarily the most influential in root communities. In the root microbiome, we found that the management type was the most influential factor for the root bacteria, while tillage intensity explained most of the variation in the root fungi. Taken together, our results demonstrate that agricultural management affects soil and root microbial communities differently.
We hypothesize that a combination of timing and nutrient characteristics of root and soil compartments could explain the differential responses of soil and root microbiomes to the cropping practices. By timing, we refer to the different time points between the effective cropping practices (mostly before seeding) and the harvest or sampling of the crop. We assume that the soil microbiome would exhibit the most pronounced differences in response to tillage or manure fertilization shortly after application and that such effects would gradually decline over the length of the growing season until a soil-type specific equilibrium is reached again after disturbance. For the root microbiome, however, crops are sown soon after tillage or manure application and the roots recruit microbes from the most divergent conditions so that the pronounced differences between cropping practices may be “fixed” for a longer time. Hence, primary colonizers that coin the root microbiome assembly at early stages would explain the preservation of precedent management differences. In addition, the nutritional characteristics of the root compartment may contribute to preserving specific management differences. Compared to an oligotrophic soil environment, we consider plant roots a copiotrophic compartment due to the continuous secretion of root exudates. Nutrient-rich organic fertilizers mainly contain copiotrophic gut bacteria from cattle that may also find favorable conditions in the root compartment. We see support for this idea as there is a marked impact of organic management on the root bacteria and not on soil (Figs. 2 and 3) and because the only bacteria with high degrees of co-occurrence were exclusively found in the root microbiome networks and were specific to organic farming (Fig. 5). Future experiments are needed to test these hypotheses. Such experiments would include, e.g., the quantification of cropping-induced microbiome differences of soil and root samples in time series throughout the growing season or manure application tests that uncouple nutritional from microbial components (e.g., applications of nutrient-free microbial extracts from slurry or extracts with inactivated slurry microbiota).
We confirmed the soil to be more diverse than the root microbial communities [20, 55]; however, we only found marginal impacts of the cropping systems on bacteria and fungi α-diversity (Additional file 1: Fig. S5 and Table S3, see Supplementary Discussion). Hence, the different cropping practices affected species richness to smaller degree than community composition. This is consistent with previous observations that species richness was less variable in their responses to environmental factors (i.e., different cropping systems) than species composition [56, 57]. Changes in microbial community composition may not necessarily lead to altered diversity or richness because changes of some taxonomic groups may be compensated by changes in others  and because univariate measures of diversity and richness mask relationships between individual and groups of taxa .
Cropping sensitive microbes
We identified cropping sensitive csOTUs in both soil and root microbial communities (Additional file 1: Fig. S6), and they function as indicator taxa to explain the β-diversity patterns by cropping practices (Fig. 2). For example, the higher relative abundance of bacteria csOTUs from the Firmicutes in organically managed plots (Additional file 1: Fig. S7) was congruent with a separation by management type in CAP analysis of soil and root communities. The association of Firmicutes OTUs to organic plots that receive manure fertilizer was found earlier . In our case, we noted that OTUs representing four families within the Firmicutes, Peptostreptococcaceae (bOTU23, bOTU119, genus not assigned), Clostridiaceae (bOTU341, genus Clostridium), Erysipelotrichaceae (bOTU36, genus Turicibacter), and Lachnospiraceae (bOTU1403, genus Butyrivibrio) had higher abundances in soil and root samples from organically managed plots (Additional file 1: Figs. S9 and S11). It is possible that the higher abundance of these OTUs is a direct result of manure application, as bacteria from these families have previously been isolated from cattle manure  or reported in such community surveys  and are also common in waste products of other livestock .
Although inference of ecological function from OTU data must be interpreted cautiously, we inspected the csOTUs for taxa with known functions of potential importance in agriculture. Notably in soil fungi, we found two OTUs from the genus Gibberella (fOTU57, fOTU32, family Nectriaceae) that were responsive to tillage intensities and had higher abundances in no and reduced tillage samples (Additional file 1: Fig. S10). Gibberella, specifically Gibberella zeae (fOTU57), is a teleomorph of Fusarium graminearum. This pathogen of wheat causes Fusarium head blight disease, which is responsible for wheat yield losses worldwide . Similarly, in root fungi, we noted an Alternaria OTU (fOTU63 family Pleosporaceae) with a higher abundance in C-NT samples (Additional file 1: Fig. S12). Species of this genus are also known pathogens of wheat and cause leaf blight disease . These examples could suggest that less intensive tillage systems may favor potentially pathogenic taxa. In a study examining the functional role of plant-beneficial Pseudomonads and soil suppressiveness at the FAST experiment, the soil from the O-RT plots tended to be more suppressive to the soil-borne pathogen Pythium ultimum to than the soil from C-NT plots (personal communication, Dr. M. Maurhofer, ETH Zurich). It is generally difficult to infer ecological function of a microbe solely based upon a taxonomy assignment . Thus, hypotheses about microbial functions of csOTUs need to be tested using other methods such as (meta-)genome or (meta-)transcriptome sequencing or by functional assays with isolated strains to experimentally test how the cropping sensitive microbes affect plant performance .
Cropping system effects on microbial co-occurrence
In both soil and root meta-networks, we identified modules containing high proportions of OTUs responding similarly to different cropping practices (Fig. 4, Additional file 1: Fig. S14). We observed that csOTUs grouped in distinct modules that reflected the different cropping systems. We concluded that larger groups of microbes responded in a similar manner to the specific cropping practices and therefore, clustered together in the soil and root microbial networks. The soil csOTUs exhibited low to medium degrees of co-occurrence in the soil network (Fig. 5), revealing that cropping practices did not affect the highly co-occurring soil microbes, which possibly belong to “core microbiome” members . This observation suggests that only the “accessory soil microbiome” could be manipulated through cropping practices. In contrast, the csOTUs in the root microbiome—in particular the ones that were sensitive to organic farming—included members with high degrees of co-occurrence (see keystones below). This possibly means that influential community members can also be manipulated with organic cropping practices in the root microbiome. We see additional support for this hypothesis in the observation that csOTUs also included abundant microbiome members.
Keystone taxa are thought to frequently interact with many other taxa, thereby playing an important role in the overall community [67, 68]. We found the effects of cropping system were mostly limited to non-keystone taxa despite significant effects of cropping system on β-diversity and network patterns (Figs. 2 and 4). Nevertheless, we found five keystone OTUs to be cropping sensitive in the root bacteria (Additional file 1: Table S5). Three of these—bOTU23 and bOTU119 (both Peptostreptococcaceae) and bOTU36 (Erysipelotrichaceae)—are common bacteria in cattle manure or livestock waste samples [59–61], and they had higher abundances in organically managed plots (Additional file 1: Fig. S11). This finding suggests the hypothesis that manure application to soil may introduce taxa to the root microbiome with keystone function. Hence, the possible introduction of microbes from manure and their particular influence on root microbiome functioning presents a high research priority.
It is important to stress that co-occurrence networks visualize correlative relationships between taxa that include true ecological interactions (e.g., mutualism), but also non-random processes (e.g., niche-overlap), and therefore, do not necessarily reflect direct interactions between taxa [33, 69]. Future experiments will assess whether the identified keystone or cropping sensitive species directly influence other members of the microbiome or indirectly influence host performance and fitness, thereby affecting other community members . Nevertheless, co-occurrence networks are a useful tool for exploring abundance patterns in complex microbial communities and could be useful in designing future experiments. For example, in combination with reference stocks of microbial isolates, plant growth experiments can be conducted in which the presence/absence or relative abundance of keystone taxa identified by network analysis can be manipulated and the effects on plant growth and development can be scored .
The concept of “smart farming” postulates the use of state-of-the-art (originally sensing) technology to improve the quality, quantity, and sustainability of agricultural production . Its central promises are targeted and site-specific interventions with “intelligent” agricultural management. Here, we propose that agricultural microbiota manipulations and management strategies shall also be considered as “smart farming.” The goal is to integrate beneficial plant microbiome traits (e.g., those improving plant growth, nutrient use efficiency, abiotic stress tolerance, and disease resistance) into sustainable agricultural production .
As a basis for implementing microbiota management strategies into smart cropping systems, we showed here to which extent and how the different cropping practices permit the manipulation of soil and root microbiota. The types of land management and tillage intensities had marked influence on dominant or well-connected bacteria and fungi in both soil and roots. Follow-up studies now need to identify the beneficial traits of cropping sensitive microbes in order to define the microbiome functions that can be manipulated through cropping practices.
The FAST experiment
All samples in this study were collected from the Farming Systems and Tillage (FAST) experiment near Zürich, Switzerland (47° 26′ 20″ N 8° 31′ 40″ E). For a detailed description of the FAST experiment see Wittwer et al. . Briefly, the FAST experiment comprises two replicates established beside each other on the same field. The first replicate started in summer 2009 (FAST I) and the second in summer 2010 (FAST II), following a staggered start design. The FAST experiment was designed to compare conventional (C) and organic (O) managements coupled with different tillage regimes. The FAST experiment compares the four main cropping systems C-IT, C-NT, O-IT, and O-RT. Conventional plots receive synthetic mineral fertilizers, post-emergence herbicides and pesticides and are subjected to either intensive tillage (IT) or no-tillage (NT, with additional use of glyphosate). The corresponding cropping systems are referred to as conventional with intensive tillage (C-IT) or conventional without tillage (C-NT). Organically managed plots are fertilized with cattle slurry, did not receive synthetic herbicides or pesticides, and are subjected to either intensive tillage (IT) or reduced tillage (RT). The cropping systems are referred to as organic with intensive tillage (O-IT) and organic with reduced tillage (O-RT). A full-factorial design would formally require an “O-NT” treatment instead of an “O-RT” treatment. While scientifically sound, a no-till regime under organic management is not agronomically practical because of insufficient weed control without reduced tillage. Additionally, the FAST experiments comprises four cover crop treatments that are applied at the subplot level; however, for this study we only collected root and soil samples from the cover crop treatment planted with a legume species (e.g., Vicia sp.).
Sample collection and DNA extraction
Soil and root samples from Triticum aestivum were collected at flowering stage in June 2014 from the second experimental replication (FAST II; Additional file 1: Fig. S1). The FAST experiment was cropped with the same winter wheat variety (cv. Titlis) but differed in seed coating between organic (untreated) and conventional (against seed-borne pathogens) systems (details: ). In total, 32 samples were collected (4 cropping systems (C-IT, C-NT, O-IT, O-RT) * 4 replicates * 2 sample types (soil and root)). Five soil cores (at 10–20 cm depth) were collected in each plot between wheat rows, pooled and immediately frozen at − 80 °C until DNA extraction. Additional bulk soil was collected for chemical analysis (see Additional file 1: Supplementary Methods). In each sampled subplot, whole root systems corresponding to a rooting depth of ~ 10 cm were collected from five plants and pooled. The roots were rinsed with tap water to remove soil debris, dried by blotting with sterile paper, and stored at − 80 °C until DNA extraction. Our sampling method does not discriminate between microbes inhabiting the inner root tissue and the root surface and for simplicity; we refer to these combined habitats of root-associated microbes as “root” samples.
The entire root systems were first lyophilized for 48 h and then ground to a fine powder in a ball mill. DNA was extracted from a 300-mg soil or root (dry weight) subsample using the NucleoSpin Soil DNA extraction kit (Machery-Nagel GmbH & Co. KG, Düren, Germany) according to the manufacturer’s instructions, except each sample was extracted twice and the supernatants pooled to maximize DNA yield. Extracted DNA was quantified using a Quant-iT Picogreen dsDNA Assay Kit (Invitrogen, Eugene, OR, USA) on a Varian Cary Eclipse fluorescence spectrometer (Agilent Technologies, Santa Clara, CA, USA).
PCR, library preparation, and sequencing
The 16S rRNA gene amplicon library was generated using the PCR primers 799F  and 1193R . The ITS amplicon library was generated using the PCR primers fITS7  and ITS4 . The primers were extended at the 5′end with an error-tolerant barcode for multiplexed library sequencing (Additional file 2). We refer to Additional file 1: Supplementary Methods for details in PCR setup, cycling conditions (Additional file 1: Table S1) and the protocol for library preparation. The libraries were sequenced on the MiSeq Instrument (Illumina, San Diego, USA) using a 600-cycle v3 Sequencing kit, paired-end 2 × 300 cycle sequencing mode at the Functional Genomics Center Zurich (www.fgcz.ch).
Raw reads were processed using a custom-developed bioinformatics pipeline whose command-line based script is provided as Additional file 4. Reads were pre-quality filtered and trimmed at the 3′-end to 280 bp using PRINSEQ  and then merged with FLASH . Sequences were de-multiplexed using Cutadapt  and were quality-filtered with PRINSEQ. For operational taxonomic unit (OTU) delineation the 16S rRNA gene sequences were trimmed to the fixed length of 360 bp, sorted by abundance, de-replicated, and clustered to OTUs (≥ 97%, singletons removed) with UPARSE . Chimeric sequences were screened using UCHIME  against the GOLD database  and removed. Taxonomy assignment was performed using the SILVA database (v119;  with the RDP classifier as implemented in QIIME . ITS sequences were processed similarly, except they were trimmed to 220 bp and chimeric sequences were screened against the UNITE database . Taxonomy was assigned using the UNITE database (v7.0) with the RDP classifier in QIIME.
Data analysis in R
All statistical analyses were conducted in R v3.3.0 . The R script and all necessary input files are provided as Additional file 5. Additionally, a workflow of the data analysis steps presented below and the figures generated from each step is given in Additional file 1: Fig. S2. Briefly, the bacteria OTUs (bOTUs) and taxonomy tables were filtered to exclude OTUs classified as chloroplasts and mitochondria. Similarly, fungi OTUs (fOTUs) classified as plant, protist, or whose kingdom or phylum was unassigned were removed.
Rarefaction analysis was performed in QIIME on the filtered OTU tables that were exported from R. The OTU tables were rarefied from 1000 to 37,000 (bacteria) or 1000 to 48,000 (fungi) sequences per sample with a step size of 1000 and 100 iterations at each step. Estimates of α-diversity (observed OTU richness) were calculated at each rarefaction level in QIIME (Additional file 1: Fig. S3a). We tested the effects of sample type and cropping system on observed species richness for each kingdom individually. For this, we randomly selected one file containing α-diversity estimations at 11,000 (bacteria) and 9000 (fungi) sequences per sample from QIIME. We tested for differences between soil and root sample using Student’s t test. We then assessed the effects of experimental block and cropping system on observed species richness using two-way ANOVA within each kingdom and sample type separately. Because cropping system was confounded within experimental block, we did not test for the Block*Cropping System interaction. Significant differences between cropping systems were assessed using Tukey’s Honest Significant Differences test using the R package TukeyC .
We conducted a general analysis of β-diversity on the bacterial and fungal communities comparing soil and root samples together (Additional file 1: Fig. S3a) and then subsequently we performed more specific hypothesis testing on the soil and root communities individually (Additional file 1: Fig. S3b). For the general analysis, we normalized the filtered OTU sequence counts for each microbial kingdom separately using the “trimmed means of M” (TMM) method with the BioConductor package edgeR (10) and expressed the normalized counts as relative abundance counts per million (CPM). We then performed unconstrained principle coordinates analysis (PCoA) on Bray-Curtis dissimilarities to quantify the major variance components of β-diversity in each kingdom. Ordination analyses were performed using the R package phyloseq . We tested for sample type and cropping system effects on community dissimilarity with permutational analysis of variance (PERMANOVA) using the functions adonis in the vegan package with 104 permutations .
For the separate in-depth analyses of each microbial kingdom and for each sample type (soil and root), we additionally applied the following sequence count threshold to the OTU tables: we selected OTUs with at least two sequences (avoiding single-count OTUs) in at least four samples (the number of replicates per treatment). We considered OTUs remaining after this thresholding step to be the soil and root communities. We normalized the communities using the TMM method and expressed the values as relative abundance CPM. We then performed multivariate analysis of microbial diversity based on the steps outlined by Anderson and Willis . This included a constrained analysis of principal coordinates (CAP) testing the effect of the cropping systems, statistical testing of the cropping system hypothesis, and identification of the OTUs responsible for the observed effects (see below). All ordination analyses were performed using the R package phyloseq . Statistical significance of the CAP was assessed using the permutest function in the vegan package  with 104 permutations. We tested for cropping system effects on community dissimilarity with permutational analysis of variance (PERMANOVA) and permutational analysis of multivariate dispersions (BETADISP) using the functions adonis and betadisp, respectively, in the vegan package with 104 permutations. Where applicable, pairwise differences between the cropping systems were assessed with the function pairwise.perm.manova from the package RVAideMemoire .
Identification of cropping sensitive OTUs (csOTUs)
We employed complementary approaches to identify the OTUs responsible for the observed effects. We used correlation based indicator species analysis with the R package indicspecies  to calculate the point-biserial correlation coefficient (r) of an OTU’s positive association to one or a combination of cropping systems. The analysis was conducted with 104 permutations and considered significant at p < 0.05. Additionally, we tested for differential OTU abundance between one or more of the cropping systems of soil and root communities (same thresholded OTU tables) of both kingdoms using likelihood ratio tests (LRT) with the R package edgeR . OTUs whose abundances were identified as differing between one or more of the cropping systems at a false discovery rate (FDR) corrected value of p < 0.05 were considered to be cropping system responsive. We then defined OTUs that were confirmed by both indicator species analysis and LRT as cropping sensitive OTUs (csOTUs).
We visualized the significant (p < 0.05) OTU associations to one or more of the different cropping system from the indicator species analysis using bipartite networks. The networks were constructed using the Fruchterman-Reingold layout with 104 permutations as implemented in the R package igraph .
We constructed two types of co-occurrence networks. For all networks, we utilized the TMM normalized CPM counts and conducted Spearman rank correlations between OTUs and visualized the positive, significant correlations (ρ > 0.7 and p < 0.001). All networks were visualized with the Fruchterman-Reingold layout with 104 permutations in igraph.
For the in-depth assessment of soil and root bacterial and fungal communities, we performed Spearman rank correlations between all pairs of bacteria and all pairs of fungi OTUs within the soil and root communities separately. We calculated the descriptive and topological network properties with igraph. These included: the total number of network nodes (representing OTUs), total number of edges (connections between nodes representing positive, significant correlations between OTUs), and degrees of co-occurrence (number of direct correlations to a node).
We then constructed meta-networks to visualize correlations between bacteria and fungi in the soil and root communities. For this, we combined the TMM normalized CPM counts of bacteria and fungi into separate OTU tables for the soil and root communities. We performed Spearman rank correlations between all pairs of bOTUs and fOTUs. We calculated the network properties mentioned above, and additionally, to explore community structure within the soil and root meta-networks, we identified network modules. These are substructures of nodes with a higher density of edges within groups than between them. For this we utilized the greedy optimization of modularity algorithm  as implemented in igraph.
Microbial taxa that frequently co-occur with other taxa in microbial co-occurrence networks are thought to be ecologically important and potentially play a key role within the microbiome [37, 38]. We identified keystone OTUs separately for the soil and root meta-networks and defined them as those nodes within the top 1% of node degree values of each network. We prioritized this simple definition over a more complex method (e.g., based on high degree and low betweenness centrality) because both definitions uncovered largely the same sets of keystone OTUs (data not shown).
We thank Dr. Lucy Poveda from the Functional Genomics Centre Zurich for technical support in MiSeq sequencing and Dr. Diane Bürge from Agroscope for performing the physical and chemical analyses of the soil samples.
This work was supported by a grant from the Swiss National Science Foundation SNSF (grant PDFMP3_137136) awarded to MVDH and Bernhard Schmid (University of Zurich, Switzerland). The SNSF also supports KS’s research (grant 31003A_165891).
Availability of data and materials
The raw sequencing data is available from the European Nucleotide Archive. The 16S rRNA gene library was run together with data presented in Hartman et al.  and is available under the study accession PRJEB15152 (Sample: SAMEA54297418). The MiSeq run containing the ITS reads is available under the study accession PRJEB21595. We provide the command-line script for the bioinformatic sequence processing as Additional file 4. Similarly, we provide the R script for data analysis and all necessary input files as Additional file 5.
KH, MVDH, and KS conceived the study, participated in its design, and wrote the manuscript. RAW maintains the Farming System and Tillage field experiment and helped with sampling. KH, SB, and KS conducted the experiments and analyzed the data. JCW performed the bioinformatic analysis. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Tilman D. Global environmental impacts of agricultural expansion: the need for sustainable and efficient practices. Proc Natl Acad Sci U S A. 1999;96:5995–6000.View ArticlePubMedPubMed CentralGoogle Scholar
- Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, et al. The challenge of food security. Science. 2012;327:812–8.View ArticleGoogle Scholar
- Matson P, Parton W, Power A, Swift M. Agricultural intensification and ecosystem properties. Science. 1997;277:504–9.View ArticlePubMedGoogle Scholar
- Pimentel D, Harvey C, Resosudarmo P, Sinclair K, Kurz D, McNair M, et al. Environmental and economic costs of soil erosion and conservation benefits. Science. 1995;267:1117–23.View ArticlePubMedGoogle Scholar
- Vitousek PM, Mooney HA, Lubchenco J, Melillo JM. Human domination of earth’ s ecosystems. Science. 1997;277:494–9.View ArticleGoogle Scholar
- Stoate C, Boatman ND, Borralho RJ, Carvalho CR, de Snoo GR, Eden P. Ecological impacts of arable intensification in Europe. J Environ Manag. 2001;63:337–65.View ArticleGoogle Scholar
- Bommarco R, Kleijn D, Potts SG. Ecological intensification: harnessing ecosystem services for food security. Trends Ecol Evol. 2013;28:230–8.View ArticlePubMedGoogle Scholar
- Cassman KG. Ecological intensification of cereal production systems: yield potential, soil quality, and precision agriculture. Proc Natl Acad Sci U S A. 1999;96:5952–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Reganold JP, Wachter JM. Organic agriculture in the twenty-first century. Nat Plants. 2016;2:15221.View ArticlePubMedGoogle Scholar
- Hobbs PR, Sayre K, Gupta R. The role of conservation agriculture in sustainable agriculture. Philos Trans R Soc London B Biol Sci. 2008;363:543–55.View ArticlePubMedGoogle Scholar
- Mäder P, Fliessbach A, Dubois D, Gunst L, Fried P, Niggli U. Soil fertility and biodiversity in organic farming. Science. 2002;296:1694–7.View ArticlePubMedGoogle Scholar
- Bender SF, Wagg C, van der Heijden MGA. An underground revolution: biodiversity and soil ecological engineering for agricultural sustainability. Trends Ecol Evol. 2016;31:440–52.View ArticlePubMedGoogle Scholar
- Gomiero T, Pimentel D, Paoletti MG. Environmental impact of different agricultural management practices: conventional vs. organic agriculture. CRC Crit Rev Plant Sci. 2011;30:95–124.View ArticleGoogle Scholar
- Bünemann EK, Schwenke GD, Zwieten L. Van. Impact of agricultural inputs on soil organisms—a review. Aust J Soil Res. 2006;44:379.View ArticleGoogle Scholar
- Postma-Blaauw M, de Goede R, Ecology S, February N, Jack H. Soil biota community structure and abundance under agricultural intensification and extensification. Ecology. 2010;91:460–73.View ArticlePubMedGoogle Scholar
- Esperschütz J, Gattinger A, Mäder P, Schloter M, Fließbach A. Response of soil microbial biomass and community structures to conventional and organic farming systems under identical crop rotations. FEMS Microbiol Ecol. 2007;61:26–37.View ArticlePubMedGoogle Scholar
- Treonis AM, Austin EE, Buyer JS, Maul JE, Spicer L, Zasada IA. Effects of organic amendment and tillage on soil microorganisms and microfauna. Appl Soil Ecol. 2010;46:103–10.View ArticleGoogle Scholar
- Wittwer RA, Dorn B, Jossi W, van der Heijden MGA. Cover crops support ecological intensification of arable cropping systems. Sci Rep. 2017;7 https://doi.org/10.1038/srep41911.
- van der Heijden MGA, Bardgett RD, Van Straalen NM. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett. 2008;11:296–310.View ArticlePubMedGoogle Scholar
- Bulgarelli D, Schlaeppi K, Spaepen S, Ver Loren van Themaat E, Schulze-Lefert P. Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol. 2013;64:807–38.View ArticlePubMedGoogle Scholar
- Vandenkoornhuyse P, Quaiser A, Duhamel M, Le Van A, Dufrense A. The importance of the microbiome of the plant holobiont. New Phytol. 2015;206:1196-206.Google Scholar
- Berendsen RL, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.View ArticlePubMedGoogle Scholar
- Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, et al. Defining the core Arabidopsis thaliana root microbiome. Nature. 2012;488:86–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Bulgarelli D, Rott M, Schlaeppi K, Ver Loren van Themaat E, Ahmadinejad N, Assenza F, et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature. 2012;488:91–5.View ArticlePubMedGoogle Scholar
- Hartman K, van der Heijden MG, Roussely-Provent V, Walser J-C, Schlaeppi K. Deciphering composition and function of the root microbiome of a legume plant. Microbiome. 2017;5 https://doi.org/10.1186/s40168-016-0220-z.
- Peiffer JA, Spor A, Koren O, Jin Z, Tringe SG, Dangl JL, et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc Natl Acad Sci U S A. 2013;110:6548–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Edwards J, Johnson C, Santos-Medellín C, Lurie E, Podishetty NK, Bhatnagar S, et al. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc Natl Acad Sci. 2015;112:E911–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Yeoh YK, Paungfoo-Lonhienne C, Dennis PG, Robinson N, Ragan MA, Schmidt S, et al. The core root microbiome of sugarcanes cultivated under varying nitrogen fertiliser application. Environ Microbiol. 2015;18:1338–51.View ArticlePubMedGoogle Scholar
- Zarraonaindia I, Owens SM, Weisenhorn P, West K, Hampton-Marcell J, Lax S, et al. The soil microbiome influences grapevine-associated microbiota. MBio. 2015;6:e02527–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Seghers D, Wittebolle L, Top EM, Verstraete W, Siciliano S. Impact of agricultural practices on the Zea mays L. endophytic community. Appl Environ Microbiol. 2004;70:1475–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Nannipieri P, Ascher J, Ceccherini MT, Landi L, Pietramellara G, Renella G. Microbial diversity and soil functions. Eur J Soil Sci. 2003;54:655–70.View ArticleGoogle Scholar
- Fuhrman JA. Microbial community structure and its functional implications. Nature. 2009;459:193–9.View ArticlePubMedGoogle Scholar
- Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2012;10:538–50.View ArticlePubMedGoogle Scholar
- Newman MEJ. The structure and function of complex networks. SIAM Rev. 2003;45:167–256.View ArticleGoogle Scholar
- Barberan A, Bates ST, Casamayor EO, Fierer N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012;6:343–51.View ArticlePubMedGoogle Scholar
- Hartmann M, Frey B, Mayer J, Mäder P, Widmer F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 2014;9:1177–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Agler MT, Ruhe J, Kroll S, Morhenn C, Kim ST, Weigel D, et al. Microbial hub taxa link host and abiotic factors to plant microbiome variation. PLoS Biol. 2016;14:e1002352.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Heijden MGA, Hartmann M. Networking in the plant microbiome. PLoS Biol. 2016;14:e1002378.View ArticlePubMedPubMed CentralGoogle Scholar
- Banerjee S, Kirkby CA, Schmutter D, Bissett A, Kirkegaard JA, Richardson AE. Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biol Biochem. 2016;97:188–98.View ArticleGoogle Scholar
- Anderson MJ. A new method for non parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.Google Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.View ArticlePubMedGoogle Scholar
- Sugiyama A, Vivanco JM, Jayanty SS, Manter DK. Pyrosequencing assessment of soil microbial communities in organic and conventional potato farms. Plant Dis. 2010;94:1329–35.View ArticleGoogle Scholar
- Chaudhry V, Rehman A, Mishra A, Chauhan PS, Nautiyal CS. Changes in bacterial community structure of agricultural land due to long-term organic and chemical amendments. Microb Ecol. 2012;64:450–60.View ArticlePubMedGoogle Scholar
- Schneider S, Hartmann M, Enkerli J, Widmer F. Fungal community structure in soils of conventional and organic farming systems. Fungal Ecol. 2010;3:215–24.View ArticleGoogle Scholar
- Upchurch R, Chiu CY, Everett K, Dyszynski G, Coleman DC, Whitman WB. Differences in the composition and diversity of bacterial communities from agricultural and forest soils. Soil Biol Biochem. 2008;40:1294–305.View ArticleGoogle Scholar
- Drijber RA, Doran JW, Parkhurst AM, Lyon DJ. Changes in soil microbial community structure with tillage under long-term wheat- fallow management. Soil Biol Biochem. 2000;32:1419–30.View ArticleGoogle Scholar
- Wang Z, Chen Q, Liu L, Wen X, Liao Y. Responses of soil fungi to 5-year conservation tillage treatments in the drylands of northern China. Appl Soil Ecol. 2016;101:132–40.View ArticleGoogle Scholar
- Francioli D, Schulz E, Lentendu G, Wubet T, Buscot F, Reitz T. Mineral vs. organic amendments: microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front Microbiol. 2016;7 https://doi.org/10.3389/fmicb.2016.01446.
- Zhen Z, Liu H, Wang N, Guo L, Meng J, Ding N, et al. Effects of manure compost application on soil microbial community diversity and soil microenvironments in a temperate cropland in China. PLoS One. 2014;9:e108555.View ArticlePubMedPubMed CentralGoogle Scholar
- Lentendu G, Wubet T, Chatzinotas A, Wilhelm C, Buscot F, Schlegel M. Effects of long-term differential fertilization on eukaryotic microbial communities in an arable soil: a multiple barcoding approach. Mol Ecol. 2014;23:3341–55.View ArticlePubMedGoogle Scholar
- Bulgarelli D, Garrido-Oter R, Münch PC, Weiman A, Dröge J, Pan Y, et al. Structure and function of the bacterial root microbiota in wild and domesticated barley. Cell Host Microbe. 2015;17:392–403.View ArticlePubMedPubMed CentralGoogle Scholar
- Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014;8:790–803.View ArticlePubMedGoogle Scholar
- Haney CH, Samuel BS, Bush J, Ausubel FM. Associations with rhizosphere bacteria can confer an adaptive advantage to plants. Nat Plants. 2015;1:1–9.Google Scholar
- Bonito G, Reynolds H, Robeson MS, Nelson J, Hodkinson BP, Tuskan G, et al. Plant host and soil origin influence fungal and bacterial assemblages in the roots of woody plants. Mol Ecol. 2014;23:3356–70.View ArticlePubMedGoogle Scholar
- Jumpponen A, Johnson LC. Can rDNA analyses of diverse fungal communities in soil and roots detect effects of environmental manipulations--a case study from tallgrass prairie. Mycologia. 2005;97:1177–94.PubMedGoogle Scholar
- Ernest SKM, Brown JH. Homeostasis and compensation: the role of species and resources in ecosystem stability. Ecology. 2001;82:2118–32.View ArticleGoogle Scholar
- Hartmann M, Widmer F. Community structure analyses are more sensitive to differences in soil bacterial communities than anonymous diversity indices. Appl Environ Microbiol. 2006;72:7804–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu T, Chellemi DO, Martin KJ, Graham JH, Rosskopf EN. Discriminating the effects of agricultural land management practices on soil fungal communities. Soil Biol Biochem. 2007;39:1139–55.View ArticleGoogle Scholar
- Ziemer CJ. Newly cultured bacteria with broad diversity isolated from eight-week continuous culture enrichments of cow feces on complex polysaccharides. Appl Environ Microbiol. 2014;80:574–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Shanks OC, Kelty CA, Archibeque S, Jenkins M, Newton RJ, McLellan SL, et al. Community structures of fecal bacteria in cattle from different animal feeding operations. Appl Environ Microbiol. 2011;77:2992–3001.View ArticlePubMedPubMed CentralGoogle Scholar
- Verbarg S, Göker M, Scheuner C, Schumann P, Stackebrandt E. The families Erysipelotrichaceae emend., Coprobacillaceae fam. Nov., and Turicibacteraceae fam. Nov. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes: Firmicutes and Tenericutes. 4th ed. Berlin-Heidelberg: Springer-Verlag; 2014. p. 79–105.Google Scholar
- Goswami RS, Kistler HC. Heading for disaster: Fusarium graminearum on cereal crops. Mol Plant Pathol. 2004;5:515–25.View ArticlePubMedGoogle Scholar
- Mercado Vergnes D, Renard ME, Duveiller E, Maraite H. Identification of Alternaria spp. on wheat by pathogenicity assays and sequencing. Plant Pathol. 2006;55:585–93.View ArticleGoogle Scholar
- Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Schlaeppi K, Bulgarelli D. The plant microbiome at work. Mol Plant-Microbe Interact. 2015;28:212–7.View ArticlePubMedGoogle Scholar
- Shade A, Peter H, Allison SD, Baho DL, Berga M, Bürgmann H, et al. Fundamentals of microbial community resistance and resilience. Front Microbiol. 2012;3:1–19.View ArticleGoogle Scholar
- Berry D, Widder S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front Microbiol. 2014;5:1–14.View ArticleGoogle Scholar
- Ma B, Wang H, Dsouza M, Lou J, He Y, Dai Z, et al. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 2016;10:1891–901.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiss S, Van Treuren W, Lozupone C, Faust K, Friedman J, Deng Y, et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 2016. https://doi.org/10.1669-81.
- Wolfert S, Ge L, Verdouw C, Bogaardt MJ. Big data in smart farming - a review. Agric Syst. 2017;153:69–80.View ArticleGoogle Scholar
- Busby PE, Soman C, Wagner MR, Friesen ML, Kremer J, Bennett A, et al. Research priorities for harnessing plant microbiomes in sustainable agriculture. PLoS Biol. 2017;15:e2001793.View ArticlePubMedPubMed CentralGoogle Scholar
- Chelius MK, Triplett EW. The diversity of archaea and bacteria in association with the roots of Zea Mays L. Microb Ecol. 2001;41:252–63.View ArticlePubMedGoogle Scholar
- Bodenhausen N, Horton MW, Bergelson J. Bacterial communities associated with the leaves and the roots of Arabidopsis thaliana. PLoS One. 2013;8:e56329.View ArticlePubMedPubMed CentralGoogle Scholar
- Ihrmark K, Bödeker ITM, Cruz-Martinez K, Friberg H, Kubartova A, Schenck J, et al. New primers to amplify the fungal ITS2 region--evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol Ecol. 2012;82:666–77.View ArticlePubMedGoogle Scholar
- White T, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis M, Gelfland D, Sninsky J, White T, editors. PCR Protoc. A Guid. To Methods Appl. San Diego: Academic Press; 1990. p. 315–22.Google Scholar
- Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Magoč T, Salzberg SL. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;21:2957–63.Google Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.View ArticleGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.View ArticlePubMedGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Reddy TBK, Thomas AD, Stamatis D, Bertsch J, Isbandi M, Jansson J, et al. The Genomes OnLine Database (GOLD) v.5: a metadata management system based on a four level (meta)genome project classification. Nucleic Acids Res. 2015;43:D1099–106.View ArticlePubMedGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Abarenkov K, Henrik Nilsson R, Larsson K-H, Alexander I, Eberhardt U, Erland S, et al. The UNITE database for molecular identification of fungi – recent updates and future perspectives. New Phytol. 2010;186:281–5.View ArticlePubMedGoogle Scholar
- R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015.Google Scholar
- Faria JC, Jelihovschi EG, Allaman Bezerra I. TukeyC: conventional Tukey test. 2016.Google Scholar
- McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.View ArticlePubMedPubMed CentralGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, et al. Vegan: community ecology package. 2015.Google Scholar
- Anderson MJ, Willis TJ. Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology. 2003;84:511–25.View ArticleGoogle Scholar
- Hervé M. RVAideMemoire: diverse basic statistical and graphical functions. 2016.Google Scholar
- De Cáceres M, Legendre P, Moretti M. Improving indicator species analysis by combining groups of sites. Oikos. 2010;119:1674–84.View ArticleGoogle Scholar
- Csardi G, Nepusz T: The igraph software package for complex network research, InterJournal, Complex Systems. 2006;1695. http://igraph.org.
- Clauset A, Newman MEJ, Moore C. Finding community structure in very large networks. Phys Rev E. 2004;70:66111.View ArticleGoogle Scholar