Impact Factor 4.019

The world's most-cited Microbiology journal

Original Research ARTICLE

Front. Microbiol., 22 August 2016 | https://doi.org/10.3389/fmicb.2016.01278

Spatio-Temporal Variation of Core and Satellite Arbuscular Mycorrhizal Fungus Communities in Miscanthus giganteus

Christopher J. Barnes1*†, Caitlin A. Burns1, Christopher J. van der Gast2, Niall P. McNamara3 and Gary D. Bending1
  • 1School of Life Sciences, University of Warwick, Coventry, UK
  • 2NERC Centre for Ecology and Hydrology, Wallingford, UK
  • 3NERC Centre for Ecology and Hydrology – Lancaster Environment Centre, Lancaster, UK

Arbuscular mycorrhizal fungi (AMF) are a group of obligate plant symbionts which can promote plant nutrition. AMF communities are diverse, but the factors which control their assembly in space and time remain unclear. In this study, the contributions of geographical distance, environmental heterogeneity and time in shaping AMF communities associated with Miscanthus giganteus (a perennial grass originating from south-east Asia) were determined over a 13 months period. In particular, the community was partitioned into core (abundant and persistent taxa) and satellite (taxa with low abundance and persistence) constituents and the drivers of community assembly for each determined. β-diversity was exceptionally low across the 140 m line transects, and there was limited evidence of geographical scaling effects on the composition of the core, satellite or combined communities. However, AMF richness and community composition changed over time associated with fluctuation within both the core and satellite communities. The degree to which AMF community variation was explained by soil properties was consistently higher in the core community than the combined and satellite communities, suggesting that the satellite community had considerable stochasticity associated with it. We suggest that the partitioning of communities into their core and satellite constituents could be employed to enhance the variation explained within microbial community analyses.

Introduction

Rhizosphere fungi play a major role in terrestrial ecosystems, shaping plant community structure and ecosystem function though parasitism, mutualism, and saprophytism (Berendsen et al., 2012). Arbuscular mycorrhizal fungi (AMF) are one of the most studied microbial inhabitants of the rhizosphere (Verbruggen et al., 2012). These fungi, from the phylum Glomeromycota, are obligate symbionts, which can form mutualistic symbioses with most plant species (Gosling et al., 2013). AMF assimilate resources from the soil, particularly phosphorus and other nutrients, which are traded with the plant in return for carbon-based metabolites (Kiers et al., 2011; Orwin et al., 2011). Consequently AMF can be a significant sink of plant assimilates, and can have far reaching effects within ecosystems, influencing nutrient cycles, soil stability and water retention (Kucey, 1987; Hodge et al., 2001; Crow et al., 2009).

A wide range of factors are known to affect the specific composition of AMF communities which assemble in the root zone. AMF communities can vary greatly over time, particularly over seasonal time frames (Escudero and Mendoza, 2005; Dumbrell et al., 2009; Pereira e Silva et al., 2012; Bennett et al., 2013). This seasonality in AMF communities may be driven directly through changing abiotic factors associated with annual cycles, such as temperature and rainfall, or indirectly through changes in carbon availability from host plants (Dumbrell et al., 2011). However, more long-term changes in AMF community composition over years and decades may also occur in response to plant age (Daniell et al., 2001; Husband et al., 2002) and changes in aboveground biodiversity (van der Heijden et al., 1998; Wardle et al., 2004).

Arbuscular mycorrhizal fungi (AMF) communities also show significant spatial variability, responding to factors such as soil properties (Dumbrell et al., 2009; Hazard et al., 2013), vegetation and climate (Treseder et al., 2004; Tedersoo et al., 2012; Peay et al., 2013). Although some fungi have far-reaching dispersal patterns, many may face the same barriers to dispersal as larger organisms, including mountains and oceans, anthropogenic pressures, and ecological factors (Treseder and Cross, 2006; Tedersoo et al., 2014). These spatial scaling effects have been demonstrated for fungi and AMF over local, landscape and continental geographic distances (Tedersoo et al., 2003; Green et al., 2004; Lilleskov et al., 2004; Peay et al., 2007; van der Gast et al., 2011), although some studies have found limited evidence for dispersal limitation affecting AMF community assembly over regional geographical scales (An et al., 2008; Hazard et al., 2013).

The majority of microbial biogeographical studies analyze spatial or temporal scaling alone, with few considering both together (Matsuda and Hijii, 1998; Courty et al., 2008; Pereira e Silva et al., 2012). However, since microbial communities are known to display great temporal and spatial variation (Bell, 2010; Barnes et al., 2016), these variables should be considered together, and the relative effects of each elucidated.

Communities can be divided into ‘core’ species, which are locally abundant and regionally persistent, and ‘satellite’ species which are both regionally and locally rare (Hanski, 1982; Magurran, 2007). Magurran and Henderson (2003) hypothesized that in macro-organisms, core species are well-adapted to surroundings, whereas, satellite species are under limitations of dispersal. This approach has proved to be a useful tool to understand ecological principles shaping communities of macro-organisms (Pärtel et al., 2001; Unterseher et al., 2011b; Supp et al., 2015) but has only infrequently been implemented in analyses of microbial communities (Ulrich and Zalewski, 2006; Unterseher et al., 2011a; Rogers et al., 2013; Lindh et al., 2015).

Miscanthus giganteus is a perennial grass that has been widely grown as a bioenergy crop within Europe (Kahle et al., 2001). Plantations remain untilled throughout their lifetime, thus soil communities can develop over several decades, in a soil matrix which remains relatively undisturbed (Christian et al., 2008). Diverse AMF communities have been found to associate with Miscanthus grasses (An et al., 2008). Miscanthus plantations provide an excellent system to study the spatial and temporal scaling of AMF, the homogenous plantations limiting confounding effects of disturbance, host species and genotypic variation on AMF communities (Corredor et al., 2014).

The aim of this study was to investigate the relative importance and interactions of spatial and temporal drivers of community assembly in the AMF communities of M. giganteus, with a particular focus on the assembly of core and satellite communities.

Materials and Methods

Study Site and Root Sampling

The study site was a commercial perennial M. giganteus bioenergy plantation, which was established in 2006 near Lincoln, UK. Prior to growing M. giganteus, the field had been historically cropped with an oilseed rape-wheat rotation. The M. giganteus was planted at a density of 10,000 rhizomes ha-1. There was no nitrogen fertilization applied during or subsequent to establishment. The soil was a sandy loam with 53% sand, 32% silt, and 15% clay.

Within the UK, M. giganteus growth is most rapid in spring and early summer, although maximum height is not reached until September or October (Christian et al., 2008). Soil samples were collected on 6th October 2010, and 21st June, 9th August and 10th October 2011, in order to capture AMF community development over a complete growing season, and to allow comparison of late growth stage communities across 2 years. Sampling was performed along line transects (Figure 1), starting at the western edge of the field site and heading due east. An initial sampling site was located 25 m into the field to avoid edge effects, and each of the seven subsequent sampling sites was located at 20 m intervals along the transect. At each sampling site, four subsamples were taken 1 m from the central position in north, south, east, and west directions. In order to avoid resampling previously disturbed sites, entire transects were shifted 3 m east at each time point.

FIGURE 1
www.frontiersin.org

FIGURE 1. Overview of a line transect. Sampling started 25 m within the field site to avoid edge effects, whilst sampling locations were 20 m apart. At each sampling location 4 cm × 15 cm deep soil cores were taken (subsamples) and analyzed independently.

At each subsampling position, a soil core (0–15 cm) was collected and stored at 4°C for a maximum of 1 week prior to processing. Miscanthus produces a highly branched, dense network of roots growing from a rhizome. We collected roots under approximately 2 mm from the soil using tweezers. Healthy white roots were separated from darker senescent Miscanthus roots, and adhering soil was washed off using deionized water. The roots were cut to 1 cm lengths, thoroughly mixed, and stored at either -20°C, or in 20% (v/v) ethanol at 4°C, for molecular and staining analyses respectively.

Soil Analysis

Soil was oven dried at 60°C prior to undergoing analysis, as previously described (Barnes et al., 2016). Briefly, soil pH was measured with a Russell model RL 150 pH meter, in a solution ratio of 1:5 with deionized water. NO3- and NH4+ were extracted from soil with 0.5 M K2SO4 and quantified using a Foss FIAStar 5000 flow injection analyser. % C and % N were measured by combustion using a Leco CN2000 analyser. Available P was extracted with 0.5 M NaHCO3, while available K and Mg were extracted with 1 M NH4NO3. Available P, K, and Mg were subsequently quantified using a Jobin Yvon Ultima 2 inductively coupled plasma optical emission spectrometer.

Mycorrhizal Staining and Colonization

Roots were stained to visualize AMF and to calculate root length colonization (RLC) by AMF (Grace and Stribley, 1991). RLC was estimated in each sample using 10 fields of view from 10 roots, using methods previously described (Brundrett et al., 1996).

DNA Extraction and PCR

Roots were washed, blotted dry, and 200 mg added to lysis tubes provided with the PowerSoil DNA isolation kit (MP Biomedicals, Cambridge, UK), and mechanically lysed in a TissueLyser (QAIGEN, UK) using three separate 30 s pulses at 30 Hz, before undergoing extraction as per manufacturers instructions. 18S rDNA fragments were amplified using the AMF specific primers AML1 (5′-ATCAACTTTCGATGGTAGGTAGGATAGA-3′) and AML2 (5′-GAACCCAAACACTTTGGTTTCC-3′) fluorescently labeled with 6-FAM (Lee et al., 2008). AMF PCR cycle conditions were as follows: 94°C 5 min, 40 cycles of 94°C 1 min, 57°C 1 min, 72°C 1 min, then 72°C 10 min, 12°C. PCR products were purified using a QIA-quick PCR purification kit.

Clone Library and Processing

To determine the most efficient restriction enzymes for terminal restriction fragment length polymorphism (TRFLP) analysis, and confirm specificity of the AMF PCR, AMF community 18S rRNA genes were sequenced. Two clone libraries established, one for the October-10 transect and another for June-11. For each library, 1 μl of equilibrated (25 ng/μl) DNA from each of the 32 subsamples was combined to produce a pooled sample. 18S rDNA was amplified from each pooled sample using the method described above, with unlabeled AML1 and AML2 primers. PCR products were purified using QIAquick Gel Extraction Kit (QIAGEN House, Sussex, UK), following the manufacturer’s protocol. Cloning was carried out using the QIAGEN PCR cloning plus kit (QIAGEN House, Sussex, UK), and TempliPHI Amplification kit (GE Healthcare Ltd, UK), following manufacturer’s protocols. Products were sequenced on an Applied Biosystems 3130X1 automated capillary sequencer. Sequences were quality checked, leaving 84 sequences in October-10 and 87 sequences in June-11, which were then aligned, and used for hierarchical clustering at 99% similarity using DNASTAR 9 (Lasenger, Madison, WI, USA). A total of seven operational taxonomic units (OTUs) were identified and taxonomy was assigned by BLAST searching against the nucleotide collection (nr/nt) database of the NCBI. Sequence accession numbers (KU937164, KU937170, KU937161, KU937154, KU937168, KU937155, KU937150), relative abundances and taxonomy are given in Supplementary Table S1.

Restriction enzymes were chosen by fragmenting clone library sequences in silico with all restriction enzymes listed on the bioperl Restriction Enzyme Mapping Application (REMA) http://bioperl.macaulay.ac.uk. The restriction enzyme which generated the highest terminal restriction fragment (TRF) richness from the clone library, Ase1 (New England Biolabs, UK), was chosen for AMF 18S rDNA. However these fragments could not be subsequently linked to specific TRFs, as found in other studies (Gosling et al., 2013).

Fragment Length Analysis

Terminal restriction fragment length polymorphism was performed on the 6-FAM labeled AMF community PCR products to assess diversity of AMF. TRFLP reactions contained 300 ng of PCR product, 10% buffer and two units of restriction enzyme, made up to 30 μl with pure water. Samples were incubated for 4 h at 37°C, followed by 20 min at 65°C to denature enzymes. Mixtures of 1 μl digestion product, 9.85 μl formamide, and GeneScanTM Size Standard (LIZ® 1200) were denatured for 5 min at 95°C, cooled on ice, and analyzed on an Applied Biosystems 3130X1 automated capillary sequencer. TRF profiles were analyzed using GeneMarker v1.8 (SoftGenetics, Stat College, PA, USA) and characterized by peak area. TRFs above 1% relative abundance in a given sample, and between 100 and 550 bp, were included in the analysis.

Partitioning of Core and Satellite Communities

‘Core’ TRFs, were differentiated from ‘satellite’ TRFs using a persistence threshold (Figure 2; Unterseher et al., 2011b; Siqueira et al., 2012). At every time point, persistence of each TRF (i.e., the number of samples it was found in) was plotted against relative abundance. TRFs found in >75% samples were classified as core taxa, since this approximately the point of inflection of the exponential curve in which TRF abundance showed substantially higher persistence than TRF found in <75% of samples. Those TRF found in <75% of samples, were classified as satellite (low abundance and persistence).

FIGURE 2
www.frontiersin.org

FIGURE 2. Persistence and abundance of the twenty-one AMF TRFs, at each sampling time. Where the vertical line is the threshold for ‘core’ taxa at 75% persistence.

Statistical Analyses

Multiple pairwise t-test were performed to assess differences in α-diversity (i.e., TRF richness) of the combined, core and satellite AMF communities, and also to compare RLC and soil properties at each time point. β-diversity within the combined, core and satellite communities was assessed by constructing distance-decay relationships (DDRs) for each transect using the power law equation S = cDd, where S = pair-wise similarity for any two samples using Jaccard (SJ) binary (i.e., presence/absence) indices of similarity, c = constant, D = distance in meters between pair-wise samples (created via a Euclidean distance matrix), and d = rate of decay in similarity (Green et al., 2004). A Bray–Curtis similarity matrix was constructed for soil nutrients (N) before partial Mantel statistics were performed to further test for an independent geographical distance effect for each transect. In these analyses, the community matrix for each season was correlated against the distance matrix, whilst the soil compositional matrix was used to compensate for changing soil properties using the equation r(SD.N) for each season (Legendre and Fortin, 1989). Jaccard binary dissimilarity matrices, Euclidean distance matrices and Bray–Curtis matrices were constructed using the vegan package within R (Oksanen et al., 2016).

To determine whether communities were stable across time, the variability between the whole, core and satellite AMF composition (Jaccard presence/absence similarity matrixes of TRF data) at each sampling time was analyzed using an analysis of similarity (ANOSIM), and similarity percentages (SIMPER) calculated using Primer6 (version 1.0.2; Primer-E, Ivybridge, UK). Community matrices were performed using data on a presence/absence basis to limit the effects of primer bias on community analyses (Anderson et al., 2003). Further tests of spatial variation were performed via PERMANOVA using the adonis function within the vegan package of R (Oksanen et al., 2016), with community data from TRFLP analysis again transformed into Jaccard similarity indexes (using presence/absence data). These indices underwent analysis against soil properties (pH, C, N, NH4, NO3, P, K, Mg), using the sampling site (each consisting of four subsamples) as the strata in the individual transect analyses, and the date transects were performed in the combined analyses across all sampling times.

Results

Soil Properties

Soil properties were stable over time (Figure 3), with the exception of NH4+ which was significantly higher in June-11 compared to the other sampling times (P < 0.001). All parameters showed variability within the transects. Soil pH varied between 5.41 and 6.80, and had a mean of 6.11. Meanwhile, K varied between 44.9 (mg/kg-1) and 267.7 (mg/kg-1) and had a mean of 199.7 (mg/kg-1). Soil C varied between 1.29 and 2.90% and had a mean of 1.83%, whilst %N ranged between 0.06 and 0.22% with a mean of 0.122%. Average Soil NH4+ content was 3.13 (mg/kg-1) and ranged between 0.71 (mg/kg-1) and 9.70 (mg/kg-1). Mg ranged between 0.9 and 476.3 (mg/kg-1), with an average of 198.8 (mg/kg-1), and finally P varied from 14.4 (mg/kg-1) to 56.0 (mg/kg-1) and had an average content of 31.5 (mg/kg-1).

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) pH, (B) % C, (C) % N, (D) NH4+ (mg kg-1), (E) NO3- (mg kg-1), (F) P (mg kg-1), (G) K (mg kg-1), and (H) Mg (mg kg-1) at each sampling time point (the mean of four sub samples). Error bars represent 1 standard error of the mean.

Root Length Colonization by Arbuscular Fungi

Stained coils, arbuscules, hyphae, and spores were observed under a light microscope in all samples, confirming the presence of AMF in M. giganteus roots. RLC ranged from 17 to 80%, and varied across sampling times (Figure 4). Mean RLC was greatest in June-11 (68.1%) followed by August-11 (61.4%). Mean RLC of October 2010 and October 2011 was 34.9, and 37.4% respectively, which were both significantly lower than either of the summer time points (P < 0.001 for each comparison), but not significantly different to each other.

FIGURE 4
www.frontiersin.org

FIGURE 4. Boxplot of AMF root length colonization; October-10, June-11, August-2011, and October-11, where a, b, and c are significantly different means (t-test).

Sequencing of the AMF Community

After hierarchical clustering, seven AMF OTUs were identified, with no non-AMF OTUs found (Supplementary Table S1). Six OTUs were present in each library. Sequences assigned to two Glomus OTUs were the most abundant, with 36.0 and 15.0% of sequences assigned to Glomus sp. M09 in October-10 and June-11 respectively, whilst 27.9 and 36.3% of sequences were assigned to Glomus sp. NBR8. Further abundant sequences were related to Paraglomus laccatum, Diversispora sp. EE1 and a Glomeromycota sp. MIB8442 that were present in both sampling points, and accounted for an average of 23.3, 8.6, and 8.5% of total reads respectively over time. A further Paraglomus laccatum OTU was present in June-11, accounting for just 5% of reads, whilst an OTU related to Glomeromycota sp. WR864 was found in October-10 only, and accounted for just 1.2% of reads.

Core and Satellite Partitioning of AMF Communities

Twenty one TRFs were identified using TRFLP. Numbers of TRFs detected varied across time points, with averages of 11.0 in October-10 (ranging from 5 to 17 TRFs), 13.6 in June-11 (ranging from 8 to 18 TRFs), 13.5 in August-11 (ranging from 10 to 17 TRFs), and 13.7 in October-11 (ranging from 8 to 21 TRFs). October-10 had significantly lower TRF richness than June-11 (P < 0.001), August-11 (=P < 0.001), and October-11 (P < 0.001). There were no significant differences in TRF richness between the three sampling times in 2011. There was a positive log-normal relationship between persistence and abundance of TRFs at each sampling time (Figure 2). The number of core TRFs at each sampling time varied between 5 and 9, while there were between 8 and 16 satellite TRFs (Supplementary Figure S1). Five TRFs remained consistently at >75% persistence over time. A total of 15 of the 21 TRFs were considered persistent (core) in at least one time point. Across sampling times, the core AMF taxa represented 41.7% of TRF richness but 79.1% of TRF abundance.

Spatio-Temporal Variation of AMF Communities

The 369 bp TRF was most abundant AMF within the community, with an average relative abundance of 16.3%, whilst other dominant core TRFs included 536 and 528 bp TRFs with average abundances of 12.3 and 10.8% respectively (Table 1). The relative abundance of many of TRFs was highly variable over time. For example the mean relative abundance of the 372 and 537 bp TRFs peaked in summer months, with the 372 bp TRF increasing from 0.4% in October-10, to 1.2% in June-11, 1.6% in August-11 before decreasing to 0.3% in October-11, whilst the 537 bp TRF rose from 2.8% in October-10, to 5.8% in June-11, 5.2% in August-11, before falling to 3.4% in October-11. Some TRFs remained relatively stable throughout the sampling period, such as the 369 bp TRF, which varied between 16.9 and 23.2%. There were six TRFs that were never core at any sampling time and also occurred in low abundance, ranging from a relative abundance of 0.1% to a maximum of 1.4%.

TABLE 1
www.frontiersin.org

TABLE 1. Average relative abundance of each TRF at each sampling time.

Distance-decay relationships were calculated to quantify β-diversity of the whole AMF community at each time point. Due to limited β-diversity, and particularly, low numbers of satellite community members, DDRs could not be determined for the core and satellite communities. The turnover of taxa was found to be consistently low across the whole 13 months sampling period, with decay values of just -0.061, -0.022, -0.058, and -0.027 for October-10, June-11, August-11, and October-11, respectively (Figure 5). Partial mantel statistics correlating community similarity to the distance matrix, whilst compensating for variation in soil nutrients, demonstrated significant independent geographical distance effects for both June-11 and October-11 (R = 0.195, P = 0.009 and R = 0.262, P = 0.002 respectively, i.e., communities closer together were more similar, even when differences in soil properties were taken into account), but not in October-10 and August-11 (R = 0.064, P = 0.231 and R = 0.003, P = 0.453 respectively).

FIGURE 5
www.frontiersin.org

FIGURE 5. Distance-decay of the combined AMF community similarity for (A) Oct-10, (B) June-11, (C) Aug-11 and (D) Oct-11. Distance-decay values were calculated for the combined community using the formula S = cDddr.

Changes in the abundance and persistence of AMF TRFs were reflected in temporal shifts in whole, core and satellite community compositions (ANOSIM, P < 0.001 for each respectively). SIMPER analysis revealed that the whole community had an average similarity of 48.0% between time points, but varied between 40.9 and 56.4% (Table 2). The similarity between time points within the core community was considerably greater than either the whole or satellite communities, averaging 55.8% overall and ranging between 49.1 and 63.4%, whilst the satellite community had an average of just 28.6% community similarity between time points, ranging between 21.5 and 36.9%.

TABLE 2
www.frontiersin.org

TABLE 2. Percentage similarity of whole, core and satellite AMF communities between sampling times.

Regulation of Combined, Core, and Satellite Communities

The factors regulating the assembly of whole, core and satellite AMF communities at each time point were analyzed independently against soil properties (pH, C, N, NH4, NO3, P, K, Mg; Table 3). A sampling site effect (with each site consisting of the four subsamples) was also investigated using PERMANOVA. Due to strong autocorrelation, only cumulative variation was considered, as the community variation associated with individual soil properties could not be reliably disentangled. In October 2010, none of the variation in the whole and satellite communities could be explained, while 9.1% of variation in core communities could be attributed to soil properties. In contrast, in August-11, soil properties explained 54.5, 55.7, and 16.6% of variation in the whole, core and satellite communities respectively (Table 3). Meanwhile in June-11, soil properties explained 18.6, 20.7, and 9.9% of variation in the whole, core and satellite communities respectively. Similarly, variation in soil properties accounted for 12.1, 13.4, and 9.3% of variation in whole, core and satellite communities in October-11.

TABLE 3
www.frontiersin.org

TABLE 3. Adonis analyses of the variance in of whole, core and satellite AMF communities explained by soil properties at each sampling time, and all times combined.

Data from across the four sampling times were combined, and cumulative variation attributed to time and soil properties was determined. The core community had a higher percentage of variation explained than either the combined or satellite communities (27.3, 25.6, and 21.8% respectively). The core component had the least percentage of variation explained by temporal variation at 3.9%, compared to 5.0 and 6.5% for the whole and satellite communities respectively.

Discussion

In this work we show that within a monodominant M. giganteus plantation the AMF community was comprised of core and satellite taxa, the composition of which varied over time. However, composition of the core community showed less variation than the satellite community in both space and time over the 13 months sampling period. Furthermore a higher percentage of variation was explained by soil parameters at each time point for the core AMF community compared to that of either the satellite or whole communities. This greater spatio-temporal variation in the satellite community, coupled with less variation explained by measured parameters suggests a higher degree of stochasticity associated with satellite community assembly. However, it remains possible that the factors regulating the satellite community could reflect unmeasured parameters, such as climatic variables or interactions with other soil biota. In macro-organisms, Magurran and Henderson (2003) hypothesized that core species are well-adapted to their surroundings, whereas satellite species are under limitations of dispersal. We therefore suggest that the partitioning of AMF communities into core and satellite communities can be utilized to identify the community best suited to the local environment and which are also most likely to be regulated by it. Additionally, the ratio of relative abundance and richness of core and satellite community members may be a useful metric in describing community dynamics, providing insight into the degree of specialization and functional redundancy within microbial communities, yielding an approach that may ultimately improve the understanding of microbial community variation over space and time (van der Gast et al., 2011).

Whilst the core and satellite communities both showed differences in resilience to temporal variation, changes in the whole community were not purely driven by extinction and immigration of the satellite community, but largely matched dynamics of the core community. Similar results have been described before in the bacterioplankton within marine ecosystems (Lindh et al., 2016), in which only eight OTUs from a community of 1000s were persistent throughout an entire year sampling, with many OTUs core over only brief time periods. There was also considerable variation in the dynamics of individual taxa within the core and satellite communities over space and time. For example the core 372 bp TRF peaked in June-11 and August-11 whilst the core 401 bp TRF was lowest in the 2 October sampling points.

Declines in AMF RLC abundance matched the decline in TRF richness in the October-10 sampling points, but not in October-11, although it should, however, be noted that absent TRF are not necessarily extinct within the community but merely fall below the detection of the molecular analysis (Avis et al., 2006). Nonetheless There was a strong log-normal relationship between TRF persistence and abundance, as found in many population studies of macroorganisms and microorganisms (Hanski, 1982; Neubert et al., 2006; Dumbrell et al., 2011; Unterseher et al., 2011a,b). Partitioning of AMF communities into core and satellite communities provided similar persistence profiles to other studies of AMF fungi (Dumbrell et al., 2011). However, Unterseher et al. (2011a) found that AMF distribution did not fit a log-normal model, but did find that 40% of AMF taxa were found to be core, comprising of 95% of all OTUs. Although in our study AMF communities followed a log-normal core and satellite model, persistence profiles were comparable to Unterseher et al. (2011a), with 41.7% of taxa core, comprising 79.1% of AMF relative abundance.

Both soil properties and time played substantial roles in determining AMF community assembly. Variation in soil properties did have significant effects in regulating AMF communities, in agreement with an abundance of literature (Treseder et al., 2004; Cavagnaro et al., 2006; Gosling et al., 2013; Nouri et al., 2014). However, attempts to determine the edaphic properties regulating the AMF community were confounded by strong spatial autocorrelation between many of the properties. Whilst this is common in ecological studies, it is often not easily disentangled. Regardless, pH could be a significant factor in shaping the AMF community, having previously been linked to changes in AMF community composition of wild Miscanthus sinensis (An et al., 2008) and playing a near ubiquitous role in shaping microbial communities (Coughlan et al., 2000; Griffiths et al., 2011; Verbruggen et al., 2012). Geographical distance effects such as dispersal limitation have also been shown to influence fungal community assembly of other fungal groups at a local level (Green et al., 2004; Lilleskov et al., 2004; Peay et al., 2007; Barnes et al., 2016). Although, partial Mantel statistics suggest an independent distance effect within June-11 and October-11, β-diversity was exceptionally low across transects and studies of AMF biogeography at larger scales have found no evidence for geographical scaling (An et al., 2008; Hazard et al., 2013), suggesting very limited effects of dispersal limitation influencing the spatial scaling of AMF at the local level.

Both AMF community composition and the degree to which soil properties accounted for community variation were found to change substantially over the 13 months sampling period. RLC rates were highest in the summer months and lowest in the October sampling points, whilst October-10 also had significantly reduced TRF richness. Highest AMF colonization rates were found in the June-11 and August-11 samples, at peak growing times, before the onset of senescence for Miscanthus grown in temperate regions (Christian et al., 2008). This matches results from a number of previous studies in which seasonal variation in AMF communities have been suggested (Welsh et al., 2009; Tian et al., 2011; Dumbrell, 2013), although, further temporal replicates would be required to confirm this trend. The exact mechanism for temporal variation in AMF has yet to be fully determined. Importantly, soil properties remained relatively stable throughout the sampling period and could not explain temporal variation within the AMF community or changes in the influence of soil properties on AMF community composition between sampling times. Furthermore, the consistent sampling design of transects eliminated the possibility of geographical scaling effects driving the observed temporal variation within the AMF community. However, fluctuating carbon availability from host roots have been shown to influence the composition and abundance of root-associated microbes (de Neergaard and Gorissen, 2004; Jones et al., 2004), and can vary throughout the year (Saggar and Hedley, 2001) and with edaphic properties (Denef et al., 2009). Thus the observed spatial and temporal changes in AMF composition and RLC could reflect differences in exchange of resources between the symbionts.

Conclusion

By partitioning communities into their core and satellite constituents we show that both the core and satellite communities vary spatially and temporally, and that local temporal variation within fungal communities is not only driven by variation within the infrequent, satellite community members. We also demonstrate that the removal of the more stochastic satellite community improved our ability to identify the factors regulating community assembly. Due to selective PCR biases, there is a continuing limitation of using data produced via metabarcoding techniques quantitatively (Pawluczyk et al., 2015; Piñol et al., 2015; Blanckenhorn et al., 2016), this approach can, however, be implemented on a presence/absence basis, as in this study, to further the understanding of factors determining microbial community regulation.

Author Contributions

CJB analyzed data and wrote the manuscript. CAB participated in writing the manuscript and performed laboratory work. CG provided statistical and analytical support. NM contributed to the writing of the manuscript and provided key materials. GB is the senior author and designed experiments and provided funding. They also supervised research activities and participated in the writing of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank Philip Law (The Institute of Cancer Research, London) for statistical support and were assisted in the field by Lawrence Davies and Cahyo Prayogo. The Lincolnshire bioenergy field site is supported by CEH project NEC03487, Jon Finch and the landowner Jonathan Wright. This project was funded by the Natural Environment Research Council as part of the Carbo-Biocrop project (http://www.carbo-biocrop.ac.uk/).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01278

References

An, G.-H., Miyakawa, S., Kawahara, A., Osaki, M., and Ezawa, T. (2008). Community structure of arbuscular mycorrhizal fungi associated with pioneer grass species Miscanthus sinensis in acid sulfate soils: habitat segregation along pH gradients. Soil Sci. Plant Nutr. 54, 517–528. doi: 10.1111/j.1747-0765.2008.00267.x

CrossRef Full Text | Google Scholar

Anderson, I. C., Campbell, C. D., and Prosser, J. I. (2003). Potential bias of fungal 18S rDNA and internal transcribed spacer polymerase chain reaction primers for estimating fungal biodiversity in soil. Environ. Microbiol. 5, 36–47. doi: 10.1046/j.1462-2920.2003.00383.x

CrossRef Full Text | Google Scholar

Avis, P. G., Dickie, I. A., and Mueller, G. M. (2006). A ‘dirty’ business: testing the limitations of terminal restriction fragment length polymorphism (TRFLP) analysis of soil fungi. Mol. Ecol. 15, 873–882. doi: 10.1111/j.1365-294X.2005.02842.x

CrossRef Full Text | Google Scholar

Barnes, C. J., van der Gast, C. J., Burns, C. A., McNamara, N. P., and Bending, G. D. (2016). Temporally variable geographical distance effects contribute to the assembly of root-associated fungal communities. Front. Microbiol. 7:195. doi: 10.3389/fmicb.2016.00195

CrossRef Full Text | Google Scholar

Bell, T. (2010). Experimental tests of the bacterial distance–decay relationship. ISME J. 4, 1357–1365. doi: 10.1038/ismej.2010.77

CrossRef Full Text | Google Scholar

Bennett, A. E., Daniell, T. J., Öpik, M., Davison, J., Moora, M., Zobel, M., et al. (2013). Arbuscular mycorrhizal fungal networks vary throughout the growing season and between successional stages. PLoS ONE 8:e83241. doi: 10.1371/journal.pone.0083241

CrossRef Full Text | Google Scholar

Berendsen, R. L., Pieterse, C. M. J., and Bakker, P. A. H. M. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. doi: 10.1016/j.tplants.2012.04.001

CrossRef Full Text | Google Scholar

Blanckenhorn, W. U., Rohner, P. T., Bernasconi, M. V., Haugstetter, J., and Buser, A. (2016). Is qualitative and quantitative metabarcoding of dung fauna biodiversity feasible? Environ. Toxicol Chem. 35, 1970–1977. doi: 10.1002/etc.3275

CrossRef Full Text | Google Scholar

Brundrett, M., Bougher, N., Dell, B., Grove, T., and Malajczuk, N. (1996). Working with Mycorrhizas in Forestry and Agriculture. Canberra: Australian Centre for International Agricultural Research.

Google Scholar

Cavagnaro, T. R., Jackson, L. E., Six, J., Ferris, H., Goyal, S., Asami, D., et al. (2006). Arbuscular mycorrhizas, microbial mommunities, nutrient availability, and soil aggregates in organic tomato production. Plant Soil 282, 209–225. doi: 10.1007/s11104-005-5847-7

CrossRef Full Text | Google Scholar

Christian, D. G., Riche, A. B., and Yates, N. E. (2008). Growth, yield and mineral content of Miscanthus ×giganteus grown as a biofuel for 14 successive harvests. Ind. Crops Prod. 28, 320–327. doi: 10.1016/j.indcrop.2008.02.009

CrossRef Full Text | Google Scholar

Corredor, A. H., Rees, K. V., and Vujanovic, V. (2014). Host genotype and health status influence on the composition of the arbuscular mycorrhizal fungi in Salix bioenergy plantations. For. Ecol. Manag. 314, 112–119. doi: 10.1016/j.foreco.2013.12.002

CrossRef Full Text | Google Scholar

Coughlan, A. P., Dalpé, Y., Lapointe, L., and Piché, Y. (2000). Soil pH-induced changes in root colonization, diversity, and reproduction of symbiotic arbuscular mycorrhizal fungi from healthy and declining maple forests. Can. J. For. Res. 30, 1543–1554. doi: 10.1139/x00-090

CrossRef Full Text | Google Scholar

Courty, P.-E., Franc, A., Pierrat, J.-C., and Garbaye, J. (2008). Temporal changes in the ectomycorrhizal community in two soil horizons of a temperate oak forest. Appl. Environ. Microbiol. 74, 5792–5801. doi: 10.1128/AEM.01592-08

CrossRef Full Text | Google Scholar

Crow, S. E., Lajtha, K., Filley, T. R., Swanston, C. W., Bowden, R. D., and Caldwell, B. A. (2009). Sources of plant-derived carbon and stability of organic matter in soil: implications for global change. Glob. Change Biol. 15, 2003–2019. doi: 10.1111/j.1365-2486.2009.01850.x

CrossRef Full Text | Google Scholar

Daniell, T. J., Husband, R., Fitter, A. H., and Young, J. P. W. (2001). Molecular diversity of arbuscular mycorrhizal fungi colonising arable crops. FEMS Microbiol. Ecol. 36, 203–209. doi: 10.1111/j.1574-6941.2001.tb00841.x

CrossRef Full Text | Google Scholar

de Neergaard, A., and Gorissen, A. (2004). Carbon allocation to roots, rhizodeposits and soil after pulse labelling: a comparison of white clover (Trifolium repens L.) and perennial ryegrass (Lolium perenne L.). Biol. Fertil. Soils 39, 228–234. doi: 10.1007/s00374-003-0699-x

CrossRef Full Text | Google Scholar

Denef, K., Roobroeck, D., Manimel Wadu, M. C. W., Lootens, P., and Boeckx, P. (2009). Microbial community composition and rhizodeposit-carbon assimilation in differently managed temperate grassland soils. Soil Biol. Biochem. 41, 144–153. doi: 10.1016/j.soilbio.2008.10.008

CrossRef Full Text | Google Scholar

Dumbrell, A. J. (2013). “Arbuscular mycorrhizal fungi throughout the year: using massively parallel pyrosequencing to quantify spatiotemporal seasonal dynamics,” in Molecular Microbial Ecology of the Rhizosphere, ed. F. J. de Bruijn (Hoboken, NJ: John Wiley & Sons, Inc.), 1113–1122.

Google Scholar

Dumbrell, A. J., Ashton, P. D., Aziz, N., Feng, G., Nelson, M., Dytham, C., et al. (2011). Distinct seasonal assemblages of arbuscular mycorrhizal fungi revealed by massively parallel pyrosequencing. New Phytol. 190, 794–804. doi: 10.1111/j.1469-8137.2010.03636.x

CrossRef Full Text | Google Scholar

Dumbrell, A. J., Nelson, M., Helgason, T., Dytham, C., and Fitter, A. H. (2009). Relative roles of niche and neutral processes in structuring a soil microbial community. ISME J. 4, 337–345. doi: 10.1038/ismej.2009.122

CrossRef Full Text | Google Scholar

Escudero, V., and Mendoza, R. (2005). Seasonal variation of arbuscular mycorrhizal fungi in temperate grasslands along a wide hydrologic gradient. Mycorrhiza 15, 291–299. doi: 10.1007/s00572-004-0332-3

CrossRef Full Text | Google Scholar

Gosling, P., Mead, A., Proctor, M., Hammond, J. P., and Bending, G. D. (2013). Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient. New Phytol. 198, 546–556. doi: 10.1111/nph.12169

CrossRef Full Text | Google Scholar

Grace, C., and Stribley, D. P. (1991). A safer procedure for routine staining of vesicular-arbuscular mycorrhizal fungi. Mycol. Res. 95, 1160–1162. doi: 10.1016/S0953-7562(09)80005-1

CrossRef Full Text | Google Scholar

Green, J. L., Holmes, A. J., Westoby, M., Oliver, I., Briscoe, D., Dangerfield, M., et al. (2004). Spatial scaling of microbial eukaryote diversity. Nature 432, 747–750. doi: 10.1038/nature03034

CrossRef Full Text | Google Scholar

Griffiths, R. I., Thomson, B. C., James, P., Bell, T., Bailey, M., and Whiteley, A. S. (2011). The bacterial biogeography of British soils. Environ. Microbiol. 13, 1642–1654. doi: 10.1111/j.1462-2920.2011.02480.x

CrossRef Full Text | Google Scholar

Hanski, I. (1982). Dynamics of regional distribution: the core and satellite species hypothesis. Oikos 38, 210–221. doi: 10.2307/3544021

CrossRef Full Text | Google Scholar

Hazard, C., Gosling, P., van der Gast, C. J., Mitchell, D. T., Doohan, F. M., and Bending, G. D. (2013). The role of local environment and geographical distance in determining community composition of arbuscular mycorrhizal fungi at the landscape scale. ISME J. 7, 498–508. doi: 10.1038/ismej.2012.127

CrossRef Full Text | Google Scholar

Hodge, A., Campbell, C. D., and Fitter, A. H. (2001). An arbuscular mycorrhizal fungus accelerates decomposition and acquires nitrogen directly from organic material. Nature 413, 297–299. doi: 10.1038/35095041

CrossRef Full Text | Google Scholar

Husband, R., Herre, E. A., Turner, S. L., Gallery, R., and Young, J. P. W. (2002). Molecular diversity of arbuscular mycorrhizal fungi and patterns of host association over time and space in a tropical forest. Mol. Ecol. 11, 2669–2678. doi: 10.1046/j.1365-294X.2002.01647.x

CrossRef Full Text | Google Scholar

Jones, D. L., Hodge, A., and Kuzyakov, Y. (2004). Plant and mycorrhizal regulation of rhizodeposition. New Phytol. 163, 459–480. doi: 10.1111/j.1469-8137.2004.01130.x

CrossRef Full Text | Google Scholar

Kahle, P., Beuch, S., Boelcke, B., Leinweber, P., and Schulten, H.-R. (2001). Cropping of Miscanthus in Central Europe: biomass production and influence on nutrients and soil organic matter. Eur. J. Agron. 15, 171–184. doi: 10.1016/S1161-0301(01)00102-2

CrossRef Full Text | Google Scholar

Kiers, E. T., Duhamel, M., Beesetty, Y., Mensah, J. A., Franken, O., Verbruggen, E., et al. (2011). Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis. Science 333, 880–882. doi: 10.1126/science.1208473

CrossRef Full Text | Google Scholar

Kucey, R. M. N. (1987). Increased phosphorus uptake by wheat and field beans inoculated with a phosphorus-solubilizing Penicillium bilaji strain and with vesicular-arbuscular mycorrhizal fungi. Appl. Environ. Microbiol. 53, 2699–2703.

Google Scholar

Lee, J., Lee, S., and Young, J. P. W. (2008). Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. FEMS Microbiol. Ecol. 65, 339–349. doi: 10.1111/j.1574-6941.2008.00531.x

CrossRef Full Text | Google Scholar

Legendre, P., and Fortin, M.-J. (1989). Spatial pattern and ecological analysis. Vegetatio 80, 107–138. doi: 10.1007/BF00048036

CrossRef Full Text | Google Scholar

Lilleskov, E. A., Bruns, T. D., Horton, T. R., Taylor, D., and Grogan, P. (2004). Detection of forest stand-level spatial structure in ectomycorrhizal fungal communities. FEMS Microbiol. Ecol. 49, 319–332. doi: 10.1016/j.femsec.2004.04.004

CrossRef Full Text | Google Scholar

Lindh, M. V., Sjöstedt, J., Andersson, A. F., Baltar, F., Hugerth, L. W., Lundin, D., et al. (2015). Disentangling seasonal bacterioplankton population dynamics by high-frequency sampling. Environ. Microbiol. 17, 2459–2476. doi: 10.1111/1462-2920.12720

CrossRef Full Text | Google Scholar

Lindh, M. V., Sjostedt, J., Ekstam, B., Casini, M., Lundin, D., Hugerth, L. W., et al. (2016). Bimodal occupancy-frequency distributions uncover the importance of regional dynamics in shaping marine microbial biogeography. BioRxiv doi: 10.1101/039883

CrossRef Full Text | Google Scholar

Magurran, A. E. (2007). Species abundance distributions over time. Ecol. Lett. 10, 347–354. doi: 10.1111/j.1461-0248.2007.01024.x

CrossRef Full Text | Google Scholar

Magurran, A. E., and Henderson, P. A. (2003). Explaining the excess of rare species in natural species abundance distributions. Nature 422, 714–716. doi: 10.1038/nature01547

CrossRef Full Text | Google Scholar

Matsuda, Y., and Hijii, N. (1998). Spatiotemporal distribution of fruitbodies of ectomycorrhizal fungi in an Abies firma forest. Mycorrhiza 8, 131–138. doi: 10.1007/s005720050225

CrossRef Full Text | Google Scholar

Neubert, K., Mendgen, K., Brinkmann, H., and Wirsel, S. G. R. (2006). Only a few fungal species dominate highly diverse mycofloras associated with the common reed. Appl. Environ. Microbiol. 72, 1118–1128. doi: 10.1128/AEM.72.2.1118-1128.2006

CrossRef Full Text | Google Scholar

Nouri, E., Breuillin-Sessoms, F., Feller, U., and Reinhardt, D. (2014). Phosphorus and nitrogen regulate arbuscular mycorrhizal symbiosis in Petunia hybrida. PLoS ONE 9:e90841. doi: 10.1371/journal.pone.0090841

CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., et al. (2016). Vegan: Community Ecology Package.

Google Scholar

Orwin, K. H., Kirschbaum, M. U. F., St John, M. G., and Dickie, I. A. (2011). Organic nutrient uptake by mycorrhizal fungi enhances ecosystem carbon storage: a model-based assessment. Ecol. Lett. 14, 493–502. doi: 10.1111/j.1461-0248.2011.01611.x

CrossRef Full Text | Google Scholar

Pärtel, M., Moora, M., and Zobel, M. (2001). Variation in species richness within and between calcareous (alvar) grassland stands: the role of core and satellite species. Plant Ecol. 157, 205–213. doi: 10.1023/A:1013938624716

CrossRef Full Text | Google Scholar

Pawluczyk, M., Weiss, J., Links, M. G., Aranguren, M. E., Wilkinson, M. D., and Egea-Cortines, M. (2015). Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples. Anal. Bioanal. Chem. 407, 1841–1848. doi: 10.1007/s00216-014-8435-y

CrossRef Full Text | Google Scholar

Peay, K. G., Baraloto, C., and Fine, P. V. (2013). Strong coupling of plant and fungal community structure across western Amazonian rainforests. ISME J. 7, 1852–1861. doi: 10.1038/ismej.2013.66

CrossRef Full Text | Google Scholar

Peay, K. G., Bruns, T. D., Kennedy, P. G., Bergemann, S. E., and Garbelotto, M. (2007). A strong species-area relationship for eukaryotic soil microbes: island size matters for ectomycorrhizal fungi. Ecol. Lett. 10, 470–480. doi: 10.1111/j.1461-0248.2007.01035.x

CrossRef Full Text | Google Scholar

Pereira e Silva, M. C., Dias, A. C., van Elsas, J. D., and Salles, J. F. (2012). Spatial and temporal variation of archaeal, bacterial and fungal communities in agricultural soils. PLoS ONE 7:e51554. doi: 10.1371/journal.pone.0051554

CrossRef Full Text | Google Scholar

Piñol, J., Mir, G., Gomez-Polo, P., and Agustí, N. (2015). Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods. Mol. Ecol. Resour. 15, 819–830. doi: 10.1111/1755-0998.12355

CrossRef Full Text | Google Scholar

Rogers, G. B., Gast, C. J., van der Cuthbertson, L., Thomson, S. K., Bruce, K. D., Martin, M. L., et al. (2013). Clinical measures of disease in adult non-CF bronchiectasis correlate with airway microbiota composition. Thorax 68, 731–737. doi: 10.1136/thoraxjnl-2012-203105

CrossRef Full Text | Google Scholar

Saggar, S., and Hedley, C. B. (2001). Estimating seasonal and annual carbon inputs, and root decomposition rates in a temperate pasture following field 14C pulse-labelling. Plant Soil 236, 91–103. doi: 10.1023/A:1011942619252

CrossRef Full Text | Google Scholar

Siqueira, T., Bini, L. M., Roque, F. O., Marques Couceiro, S. R., Trivinho-Strixino, S., and Cottenie, K. (2012). Common and rare species respond to similar niche processes in macroinvertebrate metacommunities. Ecography 35, 183–192. doi: 10.1007/s00442-013-2849-7

CrossRef Full Text | Google Scholar

Supp, S. R., Koons, D. N., and Ernest, S. K. M. (2015). Using life history trade-offs to understand core-transient structuring of a small mammal community. Ecosphere 6, 1–15. doi: 10.1890/ES15-00239.1

CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Põlme, S., Kõljalg, U., Yorou, N. S., Wijesundera, R., et al. (2014). Global diversity and geography of soil fungi. Science 346, 1256688. doi: 10.1126/science.1256688

CrossRef Full Text | Google Scholar

Tedersoo, L., Bahram, M., Toots, M., Diédhiou, A. G., Henkel, T. W., Kjøller, R., et al. (2012). Towards global patterns in the diversity and community structure of ectomycorrhizal fungi. Mol. Ecol. 21, 4160–4170. doi: 10.1111/j.1365-294X.2012.05602.x

CrossRef Full Text | Google Scholar

Tedersoo, L., Kõljalg, U., Hallenberg, N., and Larsson, K.-H. (2003). Fine scale distribution of ectomycorrhizal fungi and roots across substrate layers including coarse woody debris in a mixed forest. New Phytol. 159, 153–165. doi: 10.1046/j.1469-8137.2003.00792.x

CrossRef Full Text | Google Scholar

Tian, H., Drijber, R. A., Niu, X. S., Zhang, J. L., and Li, X. L. (2011). Spatio-temporal dynamics of an indigenous arbuscular mycorrhizal fungal community in an intensively managed maize agroecosystem in North China. Appl. Soil Ecol. 47, 141–152. doi: 10.1016/j.apsoil.2011.01.002

CrossRef Full Text | Google Scholar

Treseder, K. K., and Cross, A. (2006). Global distributions of arbuscular mycorrhizal fungi. Ecosystems 9, 305–316. doi: 10.1007/s10021-005-0110-x

CrossRef Full Text | Google Scholar

Treseder, K. K., Masiello, C. A., Lansing, J. L., and Allen, M. F. (2004). Species-specific measurements of ectomycorrhizal turnover under N-fertilization: combining isotopic and genetic approaches. Oecologia 138, 419–425. doi: 10.1007/s00442-003-1441-y

CrossRef Full Text | Google Scholar

Ulrich, W., and Zalewski, M. (2006). Abundance and co-occurrence patterns of core and satellite species of ground beetles on small lake islands. Oikos 114, 338–348. doi: 10.1111/j.2006.0030-1299.14773.x

CrossRef Full Text | Google Scholar

Unterseher, M., Jumpponen, A., Öpik, M., Tedersoo, L., Moora, M., Dormann, C. F., et al. (2011a). Species abundance distributions and richness estimations in fungal metagenomics – lessons learned from community ecology. Mol. Ecol. 20, 275–285. doi: 10.1111/j.1365-294X.2010.04948.x

CrossRef Full Text | Google Scholar

Unterseher, M., Westphal, B., Amelang, N., and Jansen, F. (2011b). 3,000 species and no end – species richness and community pattern of woodland macrofungi in Mecklenburg-Western Pomerania, Germany. Mycol. Prog. 11, 543–554. doi: 10.1007/s11557-011-0769-7

CrossRef Full Text | Google Scholar

van der Gast, C. J., Gosling, P., Tiwari, B., and Bending, G. D. (2011). Spatial scaling of arbuscular mycorrhizal fungal diversity is affected by farming practice. Environ. Microbiol. 13, 241–249. doi: 10.1111/j.1462-2920.2010.02326.x

CrossRef Full Text | Google Scholar

van der Heijden, M. G. A., Klironomos, J. N., Ursic, M., Moutoglis, P., Streitwolf-Engel, R., Boller, T., et al. (1998). Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72. doi: 10.1038/23932

CrossRef Full Text | Google Scholar

Verbruggen, E., Van Der Heijden, M. G. A., Weedon, J. T., Kowalchuk, G. A., and Röling, W. F. M. (2012). Community assembly, species richness and nestedness of arbuscular mycorrhizal fungi in agricultural soils. Mol. Ecol. 21, 2341–2353. doi: 10.1111/j.1365-294X.2012.05534.x

CrossRef Full Text | Google Scholar

Wardle, D. A., Bardgett, R. D., Klironomos, J. N., Setälä, H., Putten, W. H., and van der Wall, D. H. (2004). Ecological linkages between aboveground and belowground biota. Science 304, 1629–1633. doi: 10.1126/science.1094875

CrossRef Full Text | Google Scholar

Welsh, A. K., Burke, D. J., Hamerlynck, E. P., and Hahn, D. (2009). Seasonal analyses of arbuscular mycorrhizae, nitrogen-fixing bacteria and growth performance of the salt marsh grass Spartina patens. Plant Soil 330, 251–266. doi: 10.1007/s11104-009-0197-5

CrossRef Full Text | Google Scholar

Keywords: Glomeromycota, fungi, Miscanthus, core, satellite, arbuscular mycorrhizal fungi, spatio-temporal variation

Citation: Barnes CJ, Burns CA, van der Gast CJ, McNamara NP and Bending GD (2016) Spatio-Temporal Variation of Core and Satellite Arbuscular Mycorrhizal Fungus Communities in Miscanthus giganteus. Front. Microbiol. 7:1278. doi: 10.3389/fmicb.2016.01278

Received: 28 April 2016; Accepted: 02 August 2016;
Published: 22 August 2016.

Edited by:

Xavier Perret, University of Geneva, Switzerland

Reviewed by:

Vera Tai, University of British Columbia, Canada
Frank O’Neill Aylward, University of Hawai’i at Manoa, USA

Copyright © 2016 Barnes, Burns, van der Gast, McNamara and Bending. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Christopher J. Barnes, c.barnes@snm.ku.dk

Present address: Christopher J. Barnes, National History Museum of Denmark, University of Copenhagen, 83 Soelvgade, Copenhagen, Denmark