Damien R Finn

Sungeun Lee

Anders Lanzén

Michel Bertrand

Graeme W Nicol

Christina Hazard
email: christina.hazard@ec-lyon.fr

Ikerbasque

Cropping systems

Keywords: Biodiversity, Community Dynamics, Conservation Agriculture, Conventional Agriculture, Arbuscular Mycorrhizal Fungi

Introduction

The soil microbial community is primarily controlled by both plant-dependent factors and soil physico-chemistry [START_REF] Philippot | Going back to the roots: the microbial ecology of the rhizosphere[END_REF], particularly physico-chemistry related to soil fertility management [START_REF] Guo | Fertilization regime has a greater effect on soil microbial community structure than crop rotation and growth stage in an agroecosystem[END_REF][START_REF] Semenov | Long-term fertilization rather than plant species shapes rhizosphere and bulk soil prokaryotic communities in agroecosystems[END_REF]. It is a complex assembly including Bacteria, Archaea, obligate plant-symbiont arbuscular mycorrhizal fungi (AMF) and saprotrophic, non-mycorrhizal fungi. Bacteria regulate many essential processes, including nitrogen fixation [START_REF] Rosswall | Microbiological regulation of the biogeochemical nitrogen cycle[END_REF], nitrification [START_REF] Kowalchuk | Ammonia-oxidizing bacteria: A model for molecular microbial ecology[END_REF], decomposition of complex organic matter (i.e. cellulose, hemicellulose, chitin) [START_REF] Schimel | Microbial control over carbon cycling in soil[END_REF] and emissions of greenhouse gases methane and nitrous oxide [START_REF] Hanson | Methanotrophic bacteria[END_REF][START_REF] Barnard | Global change, nitrification, and denitrification: A review[END_REF]. In the soil environment, methanogenic archaea produce methane and members of the phylum Thaumarchaeota oxidise ammonia [START_REF] Prosser | Relative contributions of archaea and bacteria to aerobic ammonia oxidation in the environment[END_REF] with their activity often distinct from bacterial ammonia oxidisers. AMF are vital for the growth of many plant species by facilitating nutrient uptake (e.g. phosphorus and nitrogen) from soil and transfer to roots [START_REF] Koide | Extraradical hyphae of the mycorrhizal fungus Glomus intraradices can hydrolyse organic phosphate[END_REF], and also improve soil aggregation and carbon sequestration via hyphal growth [START_REF] Borie | Effects of tillage systems on soil characteristics, glomalin and mycorrhizal propagules in a Chilean Ultisol[END_REF]. Finally, saprotrophic fungi also assist in the decomposition and recycling of complex organic matter [START_REF] Setala | Decomposition rate of organic substrates in relation to the species diversity of soil saprophytic fungi[END_REF]. The combined activity of these microorganisms benefit plant-growth, often in exchange for root exudates as an organic carbon source [START_REF] Cordovez | Ecology and Evolution of Plant Microbiomes[END_REF].

Relatively high alpha-diversity (e.g. richness, evenness) in soil microbial communities is associated with improved plant net primary production, nitrogen mineralisation and microbial biomass production, amongst other functions [START_REF] Delgado-Baquerizo | Microbial diversity drives multifunctionality in terrestrial ecosystems[END_REF]. Therefore, cropping systems that have the potential to improve alpha-diversity are desirable. Tillage, for example, has a particularly detrimental effect on AMF community richness [START_REF] Sale | Impact of conservation tillage and organic farming on the diversity of arbuscular mycorrhizal fungi[END_REF][START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF]. Mineral nitrogen fertilisation favours Actinobacteria, Proteobacteria and ammonia oxidising bacteria (AOB) over ammonia oxidising archaea (AOA) [START_REF] Fierer | Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients[END_REF][START_REF] Illescas | Effect of inorganic N top dressing and Trichoderma harzianum seed-inoculation on crop yield and the shaping of root microbial communities of wheat plants cultivated under high basal N fertilization[END_REF]. Organic fertilisation increases richness of Verrucomicrobia, Planctomycetes, Acidobacteria and if manure is applied, Firmicutes associated with the gut [START_REF] Wessen | Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management[END_REF]. Organic amendments can alleviate negative impacts of salt-stress on alphadiversity [START_REF] Szoboszlay | Contrasting microbial community responses to salinization and straw amendment in a semiarid bare soil and its wheat rhizosphere[END_REF]. Crop rotations have the potential to increase bacterial alpha-diversity, while simultaneously decreasing that of fungi [START_REF] Zhou | Soil microbial communities in cucumber monoculture and rotation systems and their feedback effects on cucumber seedling growth[END_REF]. Pesticides have the potential to decrease nitrogen-fixing and nitrifying bacteria, among others [START_REF] Johnsen | Pesticide effects on bacterial diversity in agricultural soils -a review[END_REF]. Permanent cover crops structure the community in many ways: a) their root systems change soil structure that stabilises soil moisture content and increases overall microbial richness [START_REF] Vukicevich | Cover crops to increase soil microbial diversity and mitigate decline in perennial agriculture. A review[END_REF]; b) AMFassociated with cover crops remain stable in the soil, and can even transfer nutrients to commercial crops when they are present [START_REF] Cheng | Effects of mycorrhizal roots and extraradical hyphae on N-15 uptake from vineyard cover crop litter and the soil microbial community[END_REF]; and c) the rhizosphere of cover crops can promote richness of specific root disease-suppressive pseudomonads and fungal Trichoderma spp. [START_REF] Vukicevich | Cover crops to increase soil microbial diversity and mitigate decline in perennial agriculture. A review[END_REF]. Employing cropping systems that utilise management practices aimed at maximising species richness can thus maximise the benefits of soil microbial function on plant growth.

A common practice in agricultural systems that cultivate annual crops sown and harvested within a year, e.g. wheat (Triticum aestivum L.) or rapeseed (Brassica napus), is to have a fallow period between harvest and sowing of the next season's crop. A crop is not grown within this fallow period. Due to the importance of plantmicrobe interactions in shaping the soil microbial community, the removal of an active root system at harvest and subsequent fallow period likely acts as a major disturbance event. Where cropping systems that foster higher alpha-diversity are employed, the negative effect of disturbance on community composition may be less pronounced, i.e. they would be relatively more stable. This is an example of the insurance hypothesis, whereby communities that display greater species richness should be more robust to disturbance [START_REF] Ives | Stability and species richness in complex communities[END_REF] and ultimately ensure more stable functionality. Similarly, low agriculture intensity would be expected to promote relatively higher alpha-diversity, and to yield more stable relationships between individual taxa (i.e. network connections; e.g. [START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF].

The La Cage long-term experimental agricultural site was established by the Institut National de la Recherche pour l'Agriculture, l'Alimentation et l'Environnement (INRAE) in 1998 in Versailles, France. Its purpose is to compare the performance of four cropping systems; conventional (CON) versus those with a lesser intensity of disturbance -integrated (INT), conservation (CA) and organic (ORG). CA has been shown to increase soil organic carbon stocks, without increasing soil respiration rate, relative to other cropping systems [START_REF] Autret | Alternative arable cropping systems: A key to increase soil organic carbon storage? Results from a 16 year field experiment[END_REF]. Additionally, CA improved bacterial, fungal, earthworm and arthropod biomass relative to CON, while only ORG increased bacterial biomass relative to other cropping systems [START_REF] Henneron | Fourteen years of evidence for positive effects of conservation agriculture and organic farming on soil life[END_REF]. The benefits of CA are thought to be primarily derived from the incorporation of a permanent cover crop, which is unique to the CA cropping system and absent from CON, INT and ORG. The La Cage site therefore represents a valuable long-term model system to compare how various agricultural practices affect the alpha-diversity of microbial communities and how they change over a growing season. To address this, microbial community composition analyses were performed for prokaryotes, fungi and specifically for AMF on CON, CA, INT and ORG soils overwinter (December), part-way through the growing season (May) and immediately prior to harvest (July).

As the CON, CA, INT and ORG soil communities have been under the dual selection pressures of their respective cropping system and cyclical annual harvesting for 18 years, it was expected that both factors would play a role in determining composition, and the development of potential stable relationships between taxa.

We first sought to investigate whether cropping systems have differentially affected the alpha-diversity (as richness and Shannon index) of soil microbial communities in this long-term experimental system. Secondly, we sought to determine whether stable relationships between individual taxa have emerged under different cropping systems. Finally, we investigated whether cropping systems that promoted alpha-diversity would also promote relatively more stable microbial community compositions over time. This stability was measured as compositional similarity, whereby communities that remained similar over time were considered to be relatively stable, while those that demonstrated greater dissimilarity were considered less stable.

Materials and methods

La Cage experimental site

A detailed presentation of the La Cage long term trial has been reported previously [START_REF] Debaeke | Iterative design and evaluation of rule-based cropping systems: methodology and case studies. A review[END_REF]. A diagram of the La Cage experimental design is shown in Figure 1. The site is located 15 km southwest of Paris (48•48' N, 2•08' E). The soil is a deep loamy, well-drained Luvisol [START_REF] Wrb Iwg | World Reference Base for Soil Resources 2014, Update 2015[END_REF]. Texture includes: 15 -18.4% clay (< 2 µm particle size diameter), 16.5 -20.2% fine silt (2 -20 µm), 30.3 -43.2% coarse silt (20 -50 µm) and 18.4 -31.2% sand (50 -200 µm). Table 1 outlines specific practices, crop production (wheat, pea and rapeseed t ha -1 yr -1 ) and soil properties (total organic carbon TOC, total nitrogen TN, C:N ratio, pH and cation exchange capacity CEC) that differ between the four cropping systems. Briefly, CON is characterised by frequent tillage, pesticide application and high level of mineral N fertilisation. CA is no-till, infrequent pesticide application, has a permanent cover crop and receives N from mineral sources and a legume rotation. INT involves reduced tillage, infrequent pesticide application, and receives a reduced level of mineral N fertilisation. Finally, ORG involves frequent tillage, no pesticide application, and receives N from a legume rotation. Soil physico-chemical analyses were measured by INRAE in 2014, independent of this study, as described [START_REF] Debaeke | Iterative design and evaluation of rule-based cropping systems: methodology and case studies. A review[END_REF] from two pseudo-replications from subplots of each cropping system, from two independent blocks (n = 4 per cropping system). Annual crops wheat (Triticum aestivum L.) and rapeseed (Brassica napus) were grown on independent subplots (0.56 ha) within each block and undergo annual rotations between subplots at La Cage. Rapeseed was sown in August, and wheat in October. The four cropping systems (CON, CA, INT and ORG) were applied to each of the wheat/rapeseed subplots. The cover crops under CA were red fescue (Festuca rubra) and alfalfa (Medicago sativa). The cover crop is permanently maintained, but mechanically or chemically supressed prior to sowing wheat and rapeseed (i.e. a direct seeding living mulch-based system). For CON, INT and CA, mineral N application occurred yearly, with the exception of a rotating subplot (12 m x 20 m) that did not receive mineral N that year. Mineral fertiliser was applied in February. The ORG system does not involve fertilisation from mineral N sources or manure application, rather N comes from a legume rotation. During this study, ORG had an alfalfa crop rotation (sown in August) in place of rapeseed. Crops were harvested in July.

Soil sampling

From the CON, CA and INT cropping systems, five pseudo-replicate soil samples were taken from each subplot, respectively, from each combination of crop (wheat versus rapeseed (or alfalfa for ORG)) and mineral N rotation subplot (applied in 2015 versus applied in 2016). This yielded 20 pseudo-replications per CON, CA and INT cropping system. As ORG did not involve a mineral N rotating subplot, only 10 pseudo-replicate samples (five from wheat, five from alfalfa) were collected. Soil samples, at a diameter and depth of 4 and 20 cm, were randomly collected between densely distributed plants, with plot edges avoided and at least 5 meters between samples. Soil samples were collected in December 2015 (overwinter), May 2016 (during growing season)

and July 2016 (at harvest). Thus, a total of 70 samples were collected at each of the three time points, for a total of 210 samples. Finally, two spatially independent blocks were sampled. These independent blocks represented true replications of the cropping systems, as opposed to the pseudo-replicate samples taken within subplots, within each block. In total, 420 soil samples were collected for molecular analyses.

Soil samples were sieved (3.35 mm) and stored at -20°C.

Molecular analyses

DNA was extracted from 0.5 g soil using the MoBio PowerSoil DNA Kit (MoBio, Carlsbad, CA, USA) following the manufacturer's protocol. Prokaryote-and AMFspecific amplification of fragments of the small subunit rRNA gene (16S and 18S respectively) was carried out using the primer pairs 515F [START_REF] Parada | Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples[END_REF])/806R [START_REF] Apprill | Minor revisions to V4 region SSU rRNA gene primer greatly increases detection of SAR11 bacterioplankton[END_REF] and AMV4.5NF/AMDGR (Sato et al., 2005), respectively. Amplification of the fungal internal transcribed spacer 1 (ITS1) region was carried out with the primer pair ITS1F [START_REF] Gardes | ITS primers with enhanced specificity for basidiomycetes: application to the identification of mycorrhizae and rusts[END_REF])/ITS2 [START_REF] White | Amplification and direct sequencing of fungal ribosomal RNA genes for phylogentics[END_REF]. Primers had Illumina adaptor sequences attached. PCR was performed in 27 µl reactions with 22.5 µl of Invitrogen Platinum PCR SuperMix (Thermo Fisher, Carlsbad, CA, USA), 1.0 µl of forward and reverse primers (10 µM), and 2.5 µl of template DNA (5 ng/µl) on a Biometra T1 thermocycler (Biomentra GmbH, Göttingen, Germany). Thermocycling conditions for 16S rRNA gene and ITS assays were as follows: 94°C for 3 min; 35 cycles at 94°C for 45 s, 50°C for 1 min, 72°C for 90 s; and 72°C for 10 min. For 18S rRNA genes, 95°C for 5 min; 35 cycles at 95°C for 45 s, 56°C for 45 s, 72°C for 1 min; and 72°C for 7 min was used. Amplicons were bead purified using Agencourt AMPure XP (Beckman Coulter, Villepinte, France), followed by indexing PCR using the Nextera XT Index Kit (Illumina, San Diego, CA, USA) following manufacture's recommendations. Indexed amplicons were bead purified using Agencourt AMPure XP and quantified using a µDrop Plate (Thermo Fisher, Carlsbad, CA, USA).

Equimolar concentrations of samples were pooled and sequenced on an Illumina MiSeq sequencer with V2 2x150 bp paired-end chemistry. The raw sequences, with associated metadata, are available at NCBI BioProject under accession PRJN609408.

Bioinformatic analyses

Sequencing data was manually inspected using FastQC [START_REF] Andrews | FastQC: a quality control tool for high throughput sequence data[END_REF], and processed using a workflow previously described in [START_REF] Martínez-Santos | Treated and untreated wastewater effluents alter river sediment bacterial communities involved in nitrogen and sulphur cycling[END_REF], modified for fungal ITS and AMF 18S rRNA gene analyses as described in detail below. Sequence read pairs were first merged (95% of read pairs merged), using vsearch allowing for a maximum of 20 differences and "staggered" read overlaps (-fastq_maxdiffs 20 -fastq_allowmergestagger) [START_REF] Rognes | Vsearch: a versatile open source tool for metagenomics[END_REF]. Thereafter, paired reads were trimmed to remove forward and reverse primers using cutadapt [START_REF] Martin | CUTADAPT removes adapter sequences from high-throughput sequencing reads[END_REF]. Sequences lacking the correct primer sequences were discarded (using thediscard-untrimmed argument of cutadapt). Paired 16S rRNA gene sequences were then trimmed to a length of 252 and shorter sequences were discarded using vsearch (-fastq_filter -fastq_trunclen 252). For fungal ITS and AMF 18S rRNA gene sequences, with a more variable amplicon length, no trimming was carried out but sequences shorter than 100 or 200 bp, respectively, were discarded (0.43% of merged reads). Further, for all amplicons, sequences with more than one expected error were discarded, using vsearch (-fastq_maxee 1). Sequence filtering as described above resulted in 87% (10,201,515 sequences), 90% (8,788,547 sequences) and 92%

(5,463,549 sequences) of the raw 16S rRNA, ITS and 18S rRNA read pairs, respectively. Clustering and denoising was performed using Swarm v2 with default parameters and fastidious mode (-f) [START_REF] Mahe | Swarm v2: highlyscalable and high-resolution amplicon clustering[END_REF], followed by reference based and de novo chimera checking using vsearch (UCHIME algorithm; [START_REF] Edgar | UCHIME improves sensitivity and speed of chimera detection[END_REF], with the RDP Gold reference sequences for 16S rRNA, SilvaMod v106 for 18S rRNA [START_REF] Lanzen | CREST -Classification Resources for Environmental Sequence Tags[END_REF] and UNITE [START_REF] Koljalg | UNITE: a database providing webbased methods for the molecular identification of ectomycorrhizal fungi[END_REF] for ITS. Of the filtered reads, 7.5%, 0.3% and 0.3% for 16S rRNA, ITS and 18S rRNA, respectively, were chimeric. SWARM OTUs were then further subjected to clustering using vsearch with a cut-off of 97% minimum similarity, and singletons retained (0.13%, 0.01% and 0.01% of 16S rRNA, ITS and 18S rRNA reads, respectively). Taxonomic classification of representative OTU sequences was then carried out using CREST [START_REF] Lanzen | CREST -Classification Resources for Environmental Sequence Tags[END_REF] with the Silva v123 reference database [START_REF] Quast | The SILVA ribosomal RNA gene database project: improved data processing and web-based tools[END_REF], except for ITS for which the UNITE database was used [START_REF] Koljalg | UNITE: a database providing webbased methods for the molecular identification of ectomycorrhizal fungi[END_REF]. OTUs not identified as bacteria or archaea, fungi and Glomeromycota were removed from the Prokaryote 16S rRNA (0.2% of reads), fungal ITS (10% of reads) and AMF 18S rRNA (70% of reads) datasets, respectively. The MaarjAM database was used for AMF taxa identification of the OTUs to genera (similarity ≥ 97%, query coverage ≥ 90 %, E-value ≤ 1e -100 ) [START_REF] Opik | The online database MaarjAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota)[END_REF].

Statistical analyses

All statistics were performed in R v4.0.0 [START_REF] Team | R: A language and environment for statistical computing[END_REF]. Initially Tukey's Honest Significant Difference (HSD) test was performed to test for differences in soil TOC, TN, C:N ratio, pH and CEC across cropping systems. Two pseudo-replicates from each subplot, per block, for each cropping system (n = 4 each) were tested.

Samples with less than 4000 prokaryote or fungal sequences, and less than 200 AMF sequences, were removed from the datasets prior to the following analyses (1, 4 and 19 samples were removed, respectively). Heatmaps showing the relative abundance (%) of phyla within prokaryotes, fungi and AMF, for each sample across time, cropping system and blocks (n = 420) were generated with the package 'gplots' [START_REF] Warnes | Package 'gplots[END_REF]. Due to their extensive taxonomic diversity, Proteobacteria were displayed at the class level. Only taxonomic groups that comprised more than at least 5% of the total community are shown. Effects of time, cropping system, crop (wheat/rapeseed (alfalfa for ORG)), year of N application (2015/2016) in subplots and blocks (1/2) on compositional abundances at the Phylum level were performed with the ANOVA-like differential expression (ALDEx2) [START_REF] Fernandes | ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-seq[END_REF][START_REF] Fernandes | Unifying the analysis of high-throughput sequencing datasets: characterizing RNAseq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis[END_REF][START_REF] Gloor | Microbiome datasets are compositional: and this is not optional[END_REF] in R. Phyla abundances were centred log-ratio (CLR) transformed with 128 Monte-Carlo samplings and, when found to have non-normal CLR-based univariate distributions [START_REF] Korkmaz | MVN: an R package for assessing multivariate normality[END_REF], effects were tested with Kruskal-Wallis in ALDEx2 with the 'aldex.kw' function. Probability values were Benjamini-Hochberg adjusted to account for false discovery rate.

Species richness and Shannon indices presented were calculated using nonrarefied data [START_REF] Mcmurdie | Waste not, want not: why rarefying microbiome data is inadmissible[END_REF]. Rarefying had no effect on resulting trends. Indices were visualised with box and whisker plots across time and between cropping systems. Linear mixed effects models were performed to test for differences in richness and Shannon index across time and between cropping system. Random effects were considered as crop, year of N application and block. Models were applied with the 'lme4' package as described [START_REF] Bates | Fitting linear mixed-effects models using lme4[END_REF]. Principle components analysis (PCA) was performed on CLR-transformed OTU matrices based on Euclidean distance [START_REF] Vermeesch | An R package for statistical provenance analysis[END_REF] via the 'vegan' package [START_REF] Oksanen | Vegan: Community Ecology Package. R package version 2[END_REF]. Zero handling of sparse matrices was corrected by adding 1 to all cells. Sample sizes were as above, with time, cropping system, crop, year of N application and block tested on CLR-transformed OTUs with PERMANOVA separately, with 999 permutations and based on Euclidean distance. Multivariate homogeneity of variance was tested as described [START_REF] Korkmaz | MVN: an R package for assessing multivariate normality[END_REF].

For network analyses, prokaryote, fungi and AMF OTUs were subsetted based on individual management practice to compare the occurrence and stability of significant taxon-taxon relationships over time (n = 120, 120, 120, 60 for CON, CA, INT and ORG, respectively). That is, each network based on a single cropping system included each of the three time points. Weighted networks were visualised from Spearman covariance matrices, with non-significant edge weights below the 99 th quantile trimmed, performed with the 'igraph' package [START_REF] Csardi | The igraph software package for complex network research[END_REF].

Spearman correlation has been shown to yield robust networks, outperforming compositional-based correlation methods, particularly in sparsely connected networks [START_REF] Hirano | Difficulty in inferring microbial community structure based on co-occurrence network approaches[END_REF]. Due to their relatively low number of nodes and possible edges, AMF were trimmed at the 95 th quantile. Nodes that lacked at least one edge were filtered out. Weighted networks were visualised with the force-directed Kamada-Kawai algorithm. Total nodes, edges (i.e. significant relationships), edge betweenness scores, clustering coefficients and diameter of each network was calculated as described [START_REF] Csardi | The igraph software package for complex network research[END_REF]. Jaccard Similarity (JS) was calculated for each sample over time, as a pairwise measure of community similarity between harvest and overwinter (July versus December), and harvest and during growing season (July versus May) communities.

The JS was chosen as it is a robust and simple measure of pairwise similarity commonly employed in the comparison of ecological communities [START_REF] Scheiner | Measuring pattern diversity[END_REF].

By measuring the JS after disturbance (December) and part-way through the growing season (May) we sought whether communities under certain practices, such as CA, would be more similar to the harvest community than other practices, such as CON.

Where communities displayed consistently high similarity, these were considered relatively stable, while those with consistently high dissimilarity were considered unstable.

Jaccard Similarity (JS) was calculated as:

JS = ∑ " !" # ! $ # "
Eq. 1 where s was the sum of shared OTUs in samples i and j, while ti and tj were the total OTUs in samples i and j, respectively. JS was calculated at the genus level for prokaryotes and fungi, and at the individual OTU level for AMF due to their relatively low richness. The July (harvest) time point for each sample was considered i, while j iterated through the December (overwinter) and May (part-way through cropping) time points. The JS of wheat versus rapeseed/alfalfa crop or year of N application were not compared, as each sample was only compared to its corresponding cropping system/crop/year of N application sample over time. Linear mixed effects models were performed to confirm whether cropping system affected JS at December and May time points, with crop, year of N application and block considered as random effects.

Results

Soil properties and microbial community composition

Table 1 lists soil properties that differed among the four cropping systems. TOC was greatest under CA and lowest under ORG. Only CA had significantly greater TN than the other schemes. The C:N ratio, soil pH and CEC did not differ based on management. Properties did not differ between blocks (data not shown).

Time was considered as a categorical variable with three levels: December sampling during winter with sown crops, May sampling part-way through the growing season and July sampling just prior to crop harvest. Relative abundances of dominant taxonomic groups are shown in Figure 2. The prokaryotic groups with relative abundance greater than 5% were Acidobacteria, Actinobacteria, Alphaproteobacteria, Bacteroidetes, Betaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Firmicutes, Planctomycetes, Thaumarchaeota and Verrucomicrobia. These relatively abundant groups were not affected by crop, N application or block (Kruskal-Wallis p > 0.05). Both time and cropping system had effects on the relative abundance of most prokaryote phyla, with time generally conveying a stronger effect. Cropping system affected the relative abundances of the dominant fungal Ascomycota, while most other phyla were primarily affected by time. Fungi were not affected by crop, N application or block (Kruskal-Wallis p > 0.05). Glomerales and Diversisporales were the most dominant AMF. Both were affected by time, while cropping system had no effect on AMF at the Order level. The Diversisporales and Glomerales were affected by crop (Kruskal-Wallis p < 0.05) and Diversisporales by block (Kruskal-Wallis p = 0.001).

Prokaryote richness and Shannon indices followed the trend Dec < May < July regardless of cropping system (Figure 3a). Both measurements increased in May and July relative to December (p < 0.05, Table 2) while cropping system had no effect on prokaryote alpha-diversity. Only year of N application had an appreciable random effect on prokaryote alpha-diversity (8.4 and 6.3% for richness and Shannon, respectively). Fungal richness was significantly highest in July, and interestingly, decreased in May relative to December (Table 2). The CA and INT systems had positive effects on fungal richness. The fungal Shannon index indicated negative effects of ORG and time between December and May. Random effects were very marginal for fungal diversity (< 2.4% for each). The AMF richness and Shannon were most sensitive to time and cropping system (Table 2). Relative to CON, the CA, INT and ORG systems all had strong positive effects on AMF alpha-diversity which continuously increased over time. However, crop (wheat versus rapeseed (alfalfa for ORG)) played an important background role for AMF, accounting for 15.4 -22.8% of observed variation in diversity as a random effect. (circles, Figure 4a). Pairwise PERMANOVA comparisons found prokaryote July and May communities to be more similar to each other than July and December (R 2 of 0.04 and 0.14, respectively, p < 0.001). Cropping system, crop, N application and block all had minor yet significant effects (R 2 = 0.06, 0.01, 0.05 and 0.02, respectively). Both time and cropping system affected fungal communities (R 2 = 0.10 and 0.12, respectively, p < 0.001) with crop, N application and block also having relatively minor yet significant effects (R 2 = 0.04, 0.01 and 0.01, respectively). The best separation could be observed as the CA system (blue) clustering in the top quadrats and the December communities (circles) clustering in the left-hand quadrats. Pairwise PERMANOVA also found July and May fungal communities to be more similar to each other than July and December (R 2 of 0.04 and 0.11, respectively, p < 0.001). Finally, AMF were primarily affected by cropping system (R 2 = 0.10, p < 0.001) more so than time (R 2 = 0.02, p < 0.001). This was evident by the CON communities (black) separating from communities under the other three cropping systems. Crop, N application and block had little effect (R 2 = 0.02, 0.01 and 0.02, respectively). Pairwise PERMANOVA comparisons of AMF communities also found July and May communities to be more similar than July and December (R 2 of 0.01 and 0.02, respectively, p < 0.001). However, it should be noted that CLR-transformed data displayed non-homogenous variance in prokaryote, fungal and AMF multivariate matrices [START_REF] Korkmaz | MVN: an R package for assessing multivariate normality[END_REF] and consequently the reported p values, as with all p values, should be interpreted with caution [START_REF] Colquhoun | An investigation of the false discovery rate and the misinterpretation of p-values[END_REF].

Network analyses

Weighted networks of prokaryote, fungal and AMF communities were performed to visualise the stability of relationships between individual taxa under the four cropping systems, over time. Each node represents an OTU and each edge represents significant covariance. Table 3 lists node and edge numbers, mean and standard deviation of edge betweenness scores, clustering coefficients and the diameter of each network. Figure 5a are prokaryote networks for CON, CA, INT and ORG, respectively. While CON had the highest nodes, CA and INT had greater edges. Edge betweenness varied between the four systems. In the prokaryote networks, the clustering coefficients were relatively high under all cropping systems (0.51 -0.57) due to the formation of two distinct, strongly covarying clusters. Notably, one cluster tended to be dominated by Actinobacteria and Alphaproteobacteria (brown and blue nodes, respectively), whilst the other cluster tended to host a diversity of Acidobacteria, Bacteroidetes, Verrucomicrobia and other phyla. In fungal networks, CON had lower nodes than CA and INT yet higher edges. The edge betweenness was consistently low in CON, however, which indicated edge lengths between CA, INT and ORG nodes tended to be shorter, and therefore that taxa covaried more strongly.

Fungal clustering coefficients were all low, relative to prokaryotes, due to their networks tending toward a singular Ascomycota-dominated cluster with all other nodes sparsely connected and radiating outward (Figure 5b). Finally, the AMF CON network was most distinct from CA, INT and ORG with very low edges, high edge betweenness and similar clustering coefficient, despite having similar numbers of nodes. Visually this resulted in the AMF CON network having several densely packed, unconnected clusters including six nodes separate from the primary network. This differed markedly from the highly connected AMF networks under the other three systems, particularly the dense CA network. 4 shows mixed effects model results of JS differences between cropping systems at December and May timepoints. The intercepts for each mixed effects model followed the trend prokaryotes > fungi > AMF. Within each group, July (harvest) was more similar to May (during growing season) than December (overwinter). For prokaryotes, only CA had a weakly positive effect on community similarity (0.57 and 0.62, respectively) relative to CON (0.55 and 0.59, respectively). However, community similarity was strongly improved in fungal communities under CA (0.49) and INT (0.47) relative to CON (0.44) during overwinter, with this positive effect lasting under CA (0.51) relative to CON (0.49) during the growing season. The most marked differences between cropping systems were related to the AMF communities, despite having the lowest JS values and therefore having the most dissimilar compositions at overwinter and during the growing season. Similarity was greatest under CA (0.22) and ORG (0.26) relative to CON (0.15) during overwinter, which further improved in CA, INT and ORG (0.36, 0.35, 0.32 respectively) relative to CON (0.24) during the growing season.

Stability of communities over time

Random effects only explained minor variation in JS (< 10%) with the exception of crop (wheat versus rapeseed (alfalfa for ORG)) playing a larger role in background variation for AMF resistance (22.1%).

Discussion

General trends in high taxonomic ranks across time and cropping system

The removal of a crop at harvest is likely to act as a disturbance event for soil communities in agricultural systems. We sought to investigate whether cropping systems that increased microbial alpha-diversity would also improve the stability of community composition over a growing season. The four cropping systems compared here were diverse: CON with frequent tillage, frequent pesticide and mineral N application; CA with no tillage, infrequent pesticide application, an alfalfa and red fescue cover crop, and both mineral N and legume rotation for fertilisation; INT with reduced tillage, infrequent pesticide and reduced mineral N application; and ORG with frequent tillage, no pesticides and a legume rotation for fertilisation.

In general, the dominant prokaryotic phyla in these soils reflected what is known of abundant phyla in agricultural systems. Individual phyla contain organisms representing considerable metabolic and physiological diversity, but conserved traits are observed [START_REF] Finn | Functional trait relationships demonstrate life strategies in terrestrial prokaryotes[END_REF] with trends in changing abundance and diversity detected at this high taxonomic rank in soil systems. Nevertheless, it should be recognized that the behaviour of specific taxa may deviate from these general trends.

Bacteroidetes and Proteobacteria are dominant phyla in the wheat and rapeseed rhizosphere [START_REF] Ai | Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils[END_REF][START_REF] Gkarmiri | Identifying the active microbiome associated with roots and rhizosphere soil of oilseed rape[END_REF]. Unsurprisingly, Bacteroidetes, and many Proteobacteria decreased markedly in December in the absence of a mature plant host. These phyla include copiotrophic taxa that are enriched by amino acids and monosaccharides [START_REF] Fierer | Toward an ecological classification of soil bacteria[END_REF][START_REF] Goldfarb | Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance[END_REF][START_REF] Ho | Revisiting life strategy concepts in environmental microbial ecology[END_REF] which are common root exudates that tend to be less available in bulk soils [START_REF] Dennis | Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities?[END_REF][START_REF] Kuzyakov | Microbial hotspots and hot moments in soil: Concept & review[END_REF]. Conversely, the Actinobacteria thrived in December (Figure 2). While this phylum is also present in both the wheat and rapeseed rhizosphere [START_REF] Ai | Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils[END_REF][START_REF] Gkarmiri | Identifying the active microbiome associated with roots and rhizosphere soil of oilseed rape[END_REF] it includes many saprotrophs that excel at degrading and consuming dead microbial and plant biomass [START_REF] Barka | Taxonomy, physiology, and natural products of Actinobacteria[END_REF]. Thus, while the removal of the mature plant host appeared detrimental to phyla that included many plant-microbe symbionts, it was of benefit to saprotrophs. Furthermore, many Actinobacteria are well characterised psychrotolerant members of agricultural soils [START_REF] Subramanian | Cold stress tolerance in psychrotolerant soil bacteria and their conferred chilling resistance in tomato (Solanum lycoperscium Mill.) under low temperatures[END_REF] and cold deserts [START_REF] Siddiqui | Psychrophiles[END_REF] and the traits that allow these organisms to tolerate low temperatures may have contributed to their predominance in over-wintering December soils.

The Thaumarchaeota, Planctomycetes and Acidobacteria were all affected by cropping system. These three phyla are notably enriched under low-intensity practices [START_REF] Wessen | Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management[END_REF][START_REF] Ding | Effect of 35 years inorganic fertilizer and manure amendment on structure of bacterial and archaeal communities in black soil of northeast China[END_REF][START_REF] Li | Bacterial community structure after long-term organic and inorganic fertilization reveals important associations between soil nutrients and specific taxa involved in nutrient transformations[END_REF] while simultaneously decreasing under practices that involve adding mineral nitrogen [START_REF] Fierer | Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients[END_REF], Ramirez et al., 2012). Here, the relative abundance of Acidobacteria was greatest under ORG (25.9%) relative to the other systems (22.7,21.4 and 22% for CON,CA and INT,respectively). The Acidobacteria are also present in the wheat and rapeseed rhizosphere but may not necessarily grow directly on plant root exudates, as shown through isotope tracing [START_REF] Ai | Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils[END_REF][START_REF] Gkarmiri | Identifying the active microbiome associated with roots and rhizosphere soil of oilseed rape[END_REF]. It was hypothesised that the Acidobacteria are specialised to utilise recalcitrant plant organic matter in the rhizosphere, such as sloughed root cells. The capacity to utilise relatively recalcitrant plant matter would allow this taxonomic group to subsist through December until a new rhizosphere begins to take shape during the next growing season, and may explain why Acidobacteria were unaffected over the time of the growing season.

The most abundant fungal phylum was the Ascomycota (Figure 2). This phylum frequently dominates cropping and pasture systems [START_REF] De Castro | Diversity of soil fungal communities of Cerrado and its closely surrounding agriculture fields[END_REF][START_REF] Klaubauf | Molecular diversity of fungal communities in agricultural soils from Lower Austria[END_REF], Ma et al., 2013). They are a highly diverse phylum [START_REF] Lutzoni | Assembling the fungal tree of life: Progress, classification and evolution of subcellular traits[END_REF] that includes saprotrophs [START_REF] Voriskova | Fungal community on decomposing leaf litter undergoes rapid successional changes[END_REF], symbiotic endophytes of wheat and rapeseed [START_REF] Abdellatif | Endophytic hyphal compartmentalization is required for successful symbiotic Ascomycota association with root cells[END_REF][START_REF] Behie | Nutrient transfer in plant-fungal symbioses[END_REF][START_REF] Gkarmiri | Identifying the active microbiome associated with roots and rhizosphere soil of oilseed rape[END_REF], and a number of pathogens of wheat and other crops, including the genera Phaesophaeria, Alternaria and Fusarium [START_REF] Berbee | The phylogeny of plant and animal pathogens in the Ascomycota[END_REF]. Here they were particularly enriched in the CA, INT and ORG systems (58.5,58.8 and 61.6%,respectively) relative to CON (54.5%). The second most abundant fungal phyla here were the Zygomycota. This differs from other agricultural systems that show Basidiomycota as the second most dominant group [START_REF] De Castro | Diversity of soil fungal communities of Cerrado and its closely surrounding agriculture fields[END_REF][START_REF] Klaubauf | Molecular diversity of fungal communities in agricultural soils from Lower Austria[END_REF]. The Zygomycota have been noted as dominant saprotrophs when plant organic matter is depleted of lignin relative to cellulose, toward the latter stages of decomposition [START_REF] Osono | Organic chemical and nutrient dynamics in decomposing beech leaf litter in relation to fungal ingrowth and succession during 3-year decomposition processes in. a cool temperate deciduous forest in Japan[END_REF][START_REF] Vivelo | An evolutionary signal to fungal succession during plant litter decay[END_REF]. Here, the effect of time on Zygomycota favoured an increase in their abundance in May and July where a growing plant may have supplied them with lignin-depleted plant material.

Finally, analysis of AMF Glomeromycota identified Glomerales and Diversisporales to be the dominant orders. The overwhelming dominance of Glomerales in agricultural soils has been noted previously where three geographically separate soils under three different crops were dominated by > 90% Glomerales [START_REF] Helgason | Ploughing up the wood-wide web?[END_REF]. In arid agricultural soils the community can reflect roughly 50% Glomerales / 50% Diversisporales [START_REF] Li | Differences of arbuscular mycorrhizal fungal diversity and community between a cultivated land, an old field, and a never-cultivated field in a hot and arid ecosystem of southwest China[END_REF] and abiotic stress (as fungicides, pollutants) also drive communities to be dominated by Glomerales and Diversisporales [START_REF] Lenoir | Arbuscular mycorrhizal fungal responses to abiotic stresses: A review[END_REF]. However, it should be noted that certain studies have found Glomerales and Paraglomerales to dominate agricultural systems [START_REF] Sale | Impact of conservation tillage and organic farming on the diversity of arbuscular mycorrhizal fungi[END_REF][START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF]. Interestingly, brief 'blooms' of Diversisporales were apparent in December and several ORG May and July samples (Figure 2). These shifts may be a technical artefact arising from slight variation in relative abundance (%) being more obvious within the species poor (~ 250 OTUs) AMF community.

Time and cropping system effects on alpha-diversity and community composition

We first sought to investigate whether certain cropping systems would be beneficial for microbial alpha-diversity. Only time, and not cropping system, affected prokaryote alpha-diversity, with low biodiversity in December during the combined stressors of an absence of a mature plant host and low ambient temperature. Conversely, both richness and Shannon were highest in July near the culmination of harvesting with a mature plant host and mesophilic temperatures (Figure 3a). Although this study did not separately sequence the rhizosphere and bulk soil microenvironments, the presence of a mature plant likely increased the alpha-diversity, which further increased from May to July prior to harvest. Time was also the greatest determiner of prokaryote community composition, with December communities distinct from May and July (Figure 4a). It is unclear whether prokaryote biodiversity is generally greater in the rhizosphere than bulk soil [START_REF] Berg | Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere[END_REF], Philippot et al., 2013) with some studies suggesting selection for specific rhizosphere taxa reduces overall biodiversity [START_REF] Ai | Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils[END_REF] and others showing enriched biodiversity relative to bulk soil [START_REF] Mcpherson | Isolation and analysis of microbial communities in soil, rhizosphere, and roots in perennial grass experiments[END_REF]. Here, the presence of a mature crop was likely beneficial to prokaryote alpha-diversity and in shaping the community composition.

While fungal alpha-diversity was high in both December and July, cropping system played a greater role in fungal richness and Shannon index than observed for prokaryotes (Figure 3b). Specifically, CA and INT were beneficial for richness, while ORG was detrimental for diversity (Table 2). The ORG practice had significantly lower TOC than the other practices and produced the lowest crop biomass t ha -1 year -1 .

Conversely, the TOC-rich CA practice had the strongest positive effect on fungal richness and resulted in a unique community composition (Figure 4b). As soil TOC is predominantly plant C [START_REF] Dungait | Soil organic matter turnover is governed by accessibility not recalcitrance[END_REF] and assuming the majority of the Ascomycota, Basidiomycota and Zygomycota are saprotrophs reliant on plant debris from crops, the low TOC may be responsible for the decreased Shannon observed under ORG. At the global scale, climate and plant diversity are the most important drivers of fungal biodiversity [START_REF] Tedersoo | Global diversity and geography of soil fungi[END_REF]. However, at more local scales, for example within northern China and within the Mediterranean, TOC is the most important driver of fungal biodiversity [START_REF] Persiani | Mediterranean grassland soil fungi: Patterns of biodiversity, functional redundancy and soil carbon storage[END_REF][START_REF] Liu | Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China[END_REF]. The abundance of Ascomycota and Zygomycota phyla were also shown to be particularly dependent on TOC [START_REF] Liu | Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China[END_REF]. Thus, the TOC promoted by the differential cropping systems likely influenced overall fungal biodiversity.

The alpha-diversity of AMF followed clear trends of increasing from December, through May to July prior to harvest, with the CON system demonstrating consistently poor richness and Shannon indices. These effects were much stronger on the AMF than the other taxonomic groups (Table 2). Furthermore, community composition under the CON practice was most dissimilar to the other practices (Figure 4c). Tillage, and frequent pesticide application in the case of CON, have the potential to severely limit the alpha-diversity of AMF [START_REF] Sale | Impact of conservation tillage and organic farming on the diversity of arbuscular mycorrhizal fungi[END_REF][START_REF] Gottshall | Activity, diversity and function of arbuscular mycorrhizae vary with changes in agricultural management intensity[END_REF][START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF]. Conversely, the CA and ORG systems promoted AMF alpha-diversity the most, likely as a consequence of incorporating a cover crop [START_REF] Vukicevich | Cover crops to increase soil microbial diversity and mitigate decline in perennial agriculture. A review[END_REF][START_REF] Hontoria | The cover crop determines the AMF community composition in soil and in roots of maize after a ten-year continuous crop rotation[END_REF] and reduced management intensity of mineral N fertilisation and pesticide application [START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF]. The results here are in strong agreement with previous research on positive and negative drivers of AMF communities. Therefore, cropping system positively affected fungal richness, with CON being consistently relatively low and CA relatively high, while fungal Shannon diversity was low under ORG relative to CON. The AMF alpha-diversity increased under all systems relative to CON.

Cropping system effects on relationships between taxa

Another community aspect investigated was whether cropping systems could promote more stable relationships between individual taxa, which can be visualised through network analyses. Interestingly, all prokaryote networks showed the formation of two distinct clusters of significantly co-varying taxa (Figure 5a) accompanied by high clustering coefficients. Regardless of cropping system, one cluster tended toward low biodiversity dominated by potentially saprotrophic Actinobacteria and Alphaproteobacteria (the December overwinter cluster), while the other cluster tended toward higher biodiversity with all Proteobacteria classes, Bacteroidetes and Acidobacteria (the crop-associated July and May cluster). Taken together, these analyses suggest that while the absence of a mature plant and psychrotolerance potentially decreased the overall alpha-diversity in December, the prokaryote communities were flexible in that a novel community was established. Then, in the presence of a mature plant, a distinct crop-associated community emerged with increased alpha-diversity until harvest in July. It is likely that other factors associated with seasonal changes also contributed as community drivers (e.g. soil moisture).

Unlike the prokaryotes, the fungal community did not form distinct clusters of significantly covarying taxa -rather a central core of Ascomycota taxa was consistently present with other less-densely connected taxa radiating outward toward the periphery (Figure 5b). Despite being the second most dominant phylum, only few individual, scattered Zygomycota nodes were in the networks, indicating that only few of these taxa covaried significantly with others across time. As discussed above, the low TOC under ORG, which likely resulted in low biodiversity, was also likely responsible for sparse network connectivity between taxa in the ORG network, and therefore yielded fewer potential stable relationships between taxa.

Recently, [START_REF] Banerjee | Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots[END_REF] found a negative relationship between increasingly intensive agricultural management practices and node, edge and edge betweenness values within AMF community networks. Here the high-intensity CON system yielded the most distinct AMF network notable for its sparse connectivity. In contrast, the relatively low-intensity, permanent cover crop CA resulted in the most connected CA networks, indicative of strong and stable relationships between AMF over time.

Therefore, the strongest relationships between taxa were only observed for prokaryotes, and this was driven by temporal changes between overwinter versus crop-associated communities. An effect of cropping system on network structure was primarily demonstrated by the AMF CON and CA networks, which were either sparsely or highly connected, respectively.

Stability of soil microbial communities over time

The factors that confer stability to soil microbial communities have received a great deal of attention due to the role of microbes as drivers of ecosystem processes (Griffiths and Philippot, 2012, and references therein). Disturbance can come in many forms -drought [START_REF] Manzoni | Responses of soil microbial communities to water stress: results from a meta-analysis[END_REF], fire [START_REF] Banning | Effect of heat-induced disturbance on microbial biomass and activity in forest soil and the relationship between disturbance effects and microbial community structure[END_REF], metal [START_REF] Philippot | Effect of primary mild stresses on resilience and resistance of the nitrate reducer community to a subsequent severe stress[END_REF] and organic pollutants [START_REF] Girvan | Bacterial diversity promotes community stability and functional resilience after perturbation[END_REF]. Increasing biodiversity of a community can strengthen its capacity to respond to a disturbance event, termed the 'insurance hypothesis' [START_REF] Ives | Stability and species richness in complex communities[END_REF]. Here we considered crop harvesting followed by a fallow period as a disturbance, and aimed to determine whether cropping systems that promoted alpha-diversity would also demonstrate higher similarity in composition over a growing season, i.e. they would be more stable. This effect on community composition was measured with JS as a robust and simple pairwise metric of community similarity [START_REF] Scheiner | Measuring pattern diversity[END_REF] although it should be noted that it is unweighted by relative abundance. As one would expect, in all three groups, the similarity to the harvest community increased over time as communities transitioned from overwinter. Another general trend was that the composition of prokaryotes was most similar (i.e. stable) over time, fungi intermediate and AMF relatively poor (Figure 6, Intercepts in Table 4). This agreed with the insurance hypothesis at a broad taxonomic level, with the species-rich prokaryotes demonstrating greater stability than species-poor AMF, regardless of cropping system.

There was a positive benefit to stability to all groups under CA, with this effect relatively more pronounced in the fungi and AMF. Soil aggregation and TOC can provide isolated microhabitats and/or support growth of microorganisms post-disturbance, and these edaphic properties have been strongly linked to stability [START_REF] Griffiths | Biological and physical resilience of soil amended with heavy metal-contaminated sewage sludge[END_REF][START_REF] Zhang | Does microbial habitat or community structure drive the functional stability of microbes to stresses following re-vegetation of a severely degraded soil?[END_REF]. As soils under CA were no-till and had a permanent cover crop, and likely had better soil aggregation than annually tilled CON and ORG [START_REF] West | Soil organic carbon sequestration rates by tillage and crop rotation: A global data analysis[END_REF], this may explain the generally higher stability observed under this cropping system. In contrast, the combination of annual tillage, mineral N fertilisation and frequent pesticide application was likely responsible for the poor stability of AMF observed under CON. Finally, it is possible that the stability of all communities measured here are high relative to other forms of disturbance in soil systems if communities have adapted to annual crop harvesting and fallow periods [START_REF] Griffiths | Insights into the resistance and resilience of the soil microbial community[END_REF]. With the exception of the cover crop in CA soils, each practice endures an annual fallow period which has been ongoing for 18 years. This may act as a form of pressure to strictly select for microbial species that can both survive during this fallow period in the absence of an actively growing crop and respond rapidly to a new crop during the growing season. Thus, one could reasonably suppose that the communities here are somewhat adapted to recovering after crop removal. Even so, crop removal and the ensuing overwinter period had dramatic effects on the composition of all communities, relative to the harvest period, particularly for the AMF.

In summary, relative to CON, CA had a consistently positive effect on prokaryote and fungal community stability over time. The AMF demonstrated greater stability under all cropping systems relative to CON. Thus, practices that increased alpha-diversity in fungi (and AMF) also appeared to convey more stable community compositions. The increased stability in prokaryotes under CA, without increased alpha-diversity, suggests that increasing alpha-diversity is not, however, absolutely necessary for a more stable community.

Conclusions

Presented here is an investigation of prokaryote, fungal and AMF soil communities that demonstrated markedly different changes over a growing season, under four agricultural cropping systems. These cropping systems had been implemented for 18 years, giving ample time to differentially shape these communities. Prokaryotes were flexible and formed a distinct saprotroph-dominated community in the absence of a mature crop overwinter, and a crop-associated community in the presence of a mature crop. This relatively species rich group demonstrated the highest stability over time.

The CA system was associated with the most diverse fungal community, which also increased its stability over the course of the growing season. Finally, the CON practice had a particularly detrimental effect on AMF alpha-diversity and stability relative to all other cropping systems. In conclusion, cropping systems that promoted alphadiversity led to increased stability of fungal (and AMF) communities over the growing season, with prokaryotes being relatively insensitive to the effect of soil management.

Funding

This work was supported by the AXA Research Fund. Rapeseed (t ha -1 yr -1 ) 4.5 -3.8 0.8 TOC 11.1 ± 0.9 a,b,c,d 13.2 ± 2.8 a,b,c 11 ± 1.6 a,b,c,d 10.4 ± 0.5 a,b,d TN 0.9 ± 0.1 a 1.1 ± 0.2 b 0.9 ± 0.1 a 0.9 ± 0.1 a C:N ratio 12.4 ± 1 a 11.5 ± 0.1 a 11.8 ± 0.2 a 11.6 ± 0.3 a pH 7.3 ± 0.3 a 6.9 ± 0.2 a 7.1 ± 0.4 a 7 ± 0.2 a CEC 9.9 ± 1 a 9.6 ± 0.3 a 9.5 ± 1 a 9.2 ± 0.3 a Soil properties were measured in all subplots in 2014 at INRAE. 1 Tillage occurs soon after crop harvest.

Figure 4

 4 Figure 4 displays the PCA of prokaryote (a), fungal (b) and AMF (c) community compositions. Prokaryotes were primarily affected by time (PERMANOVA, R 2 = 0.13, p < 0.001) as demonstrated by distinct separation of all December communities

Figure

  Figure 6a, b and c show JS indices for prokaryote, fungal and AMF communities over time. Table4shows mixed effects model results of JS differences between cropping

Figure 1 :

 1 Figure 1: Diagram of the La Cage experimental design. Two independent blocks (1 and 2) were sampled. Each block consists of four plots (1.12 ha each) in which the four cropping systems (Conventional, CON, Conservation, CA, Integrated, INT, and Organic, ORG) were randomly applied. Each plot is split into two subplots (0.56 ha) with annual rotations of wheat and rapeseed. CA and ORG also have legume rotations. During this study, ORG had an alfalfa crop rotation. For CON, INT and CA, mineral N application occurred yearly, with the exception of a rotating subplot (12 m x 20 m) that does not receive mineral N (-N) that year (i.e. applied in 2015 versus applied in 2016). Five pseudo-replicate samples were taken from each crop (wheat/rapeseed) and N (2015/2016) subplot, for CON, CA and INT (n = 20 each) from each block. Ten pseudo-replicate samples were taken from each crop (wheat/alfalfa) subplot for ORG (n = 10) from each block.

Figure 2 :

 2 Figure 2: Heatmap comparing relative abundances of prokaryote and fungal phyla, and AMF orders, across time and between cropping system. A total of 420 samples are shown. Units are relative abundance (%) per sample. Relative abundance was calculated separately for prokaryote, fungal and AMF communities. Where taxa differed over time (T) or by cropping system (CS) Kruskal-Wallis tests on centred logratio transformed relative abundances are shown as: p < 0.05 (*); p = 0.001 (**); p < 0.001 (***).

Figure 3 :

 3 Figure 3: Observed OTU richness and Shannon indices for: a) prokaryotes, b) fungi and c) AMF.

Figure 4 :

 4 Figure 4: Ordinations of community composition for: a) prokaryotes, b) fungi and c) AMF.

Figure 5 :

 5 Figure 5: Weighted networks under CON, CA, INT and ORG for a) prokaryotes, b) fungi and c) AMF. Nodes of major taxonomic groups have been coloured and enlarged.

Figure 6 :

 6 Figure 6: Observed JS over time and between cropping system for: a) prokaryotes, b) fungi and c) AMF.

  Figure 1

Table 1 :

 1 Management practices, crop production and soil properties of the La Cage cropping systems. Superscript letters indicate Tukey's HSD outcomes between cropping systems at the p < 0.05 level.

		CON	CA	INT	ORG
	Tillage 1	Annual	None	Biannual	Annual
	N fertilisation	Mineral-high	Mineral-	Mineral-	None
			reduced	reduced
	Legume rotation	No	Yes	No	Yes
	Wheat (t ha -1 yr -1 )	9.7	6.7	8.9	5.4
	Pea (t ha -1 yr -1 )	4.2	3.7	4.5	2.6

Table 4 :

 4 Mixed effects models of cropping system effects on JS at time point, with statistical significance noted as: (***) p < 0.001, (**) p = 0.001, (*) p < 0.05. Variance explained by random effects crop, year of N application and block are also shown.

		Variance explained by random
	Cropping System	effects (%)

Acknowledgements

We would like to thank INRAE, Versailles-Grignon for use of the La Cage field trial, and Gilles Grandeau for assistance with soil sampling.

Conflict of interest statement

The authors declare no conflict of interest in regard to the outcomes of this study.