Next Article in Journal
Understanding the Role of Nitronate Monooxygenases in Virulence of the Human Fungal Pathogen Aspergillus fumigatus
Next Article in Special Issue
Changes in Chemical Structure of Thermally Modified Spruce Wood Due to Decaying Fungi
Previous Article in Journal
Cryptococcus neoformans Genotypic Diversity and Disease Outcome among HIV Patients in Africa
Previous Article in Special Issue
Metabolic Diversity of Xylariaceous Fungi Associated with Leaf Litter Decomposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fungal Community Composition as Affected by Litter Chemistry and Weather during Four Years of Litter Decomposition in Rainshadow Coastal Douglas-Fir Forests

by
Philip-Edouard Shay
1,2,
Richard S. Winder
2,
C. Peter Constabel
1 and
J. A. (Tony) Trofymow
1,2,*
1
Centre for Forest Biology, Department of Biology, University of Victoria, Victoria, BC V8W 3N5, Canada
2
Pacific Forestry Centre, Canadian Forest Service, 506 West Burnside Road, Victoria, BC V8Z 1M5, Canada
*
Author to whom correspondence should be addressed.
Submission received: 21 May 2022 / Revised: 7 July 2022 / Accepted: 13 July 2022 / Published: 16 July 2022
(This article belongs to the Special Issue Fungi in Decomposition Processes)

Abstract

:
Climate and litter chemistry are major factors influencing litter decay, a process mediated by microbes, such as fungi, nitrogen-fixing bacteria and ammonia-oxidizing bacteria. Increasing atmospheric CO2 concentrations can decrease nitrogen (N) and increase condensed tannin (CT) content in foliar litter, reducing litter quality and slowing decomposition. We hypothesized that reduced litter quality inhibits microbes and is the mechanism causing decomposition to slow. Litterbags of Douglas-fir needles and poplar leaves with a range of N (0.61–1.57%) and CT (2.1–29.1%) treatment and natural acid unhydrolyzable residue (35.3–41.5%) concentrations were placed along climatic gradients in mature Douglas-fir stands of coastal British Columbia rainshadow forests. The structure (diversity, richness and evenness) and composition of microbial communities were analyzed using DGGE profiles of 18S, NifH-universal and AmoA PCR amplicons in foliar litter after 7, 12, 24 and 43 months of decay. High CT and low N concentrations in leaf litter were associated with changes in microbial community composition, especially fungi. Contrary to our hypothesis, high CT and low N treatments did not inhibit microbial colonization or diversity. The joint effects of air temperature and soil moisture on microbial community composition at our sites were more important than the effects of initial litter chemistry.

1. Introduction

Decaying litter is one of the largest terrestrial sources of CO2 released into the atmosphere [1], impacting carbon sequestration and nutrient-cycling regimes. Litter decomposition is a process mediated by microbial communities and primarily regulated by macro- and microclimate, substrate quality, litter size, microbiome diversity and exogenous nutrient availability [2,3,4]. Climate change resulting from ongoing anthropogenic CO2 emissions is predicted to lead to a generally warmer global climate and altered moisture regimes [5] increasing decay rates [6,7] via increased microbial metabolism [1,6] and altered microbial competitive landscapes [7]. The full extent of potential climate change impacts on soil microbial communities is undetermined [8]. Faster decomposition in northern forests could release vast quantities of carbon (C) now stored in soils and decaying plant matter. For example, up to 70.6Gt could be released from Canadian forests that would exacerbate climate change [9,10].
Greater atmospheric CO2 concentrations are predicted to increase condensed tannins (CTs) and acid unhydrolyzable residues (AURs) and decrease nitrogen (N) concentrations in foliage of woody plants [11,12,13,14]. Higher CT and AUR concentrations in litter have been shown to slow the decomposition of foliage [15,16,17,18,19] and N mineralization [3,20]. These cascading effects could impact associated microbial communities, but these have not been extensively investigated. The biosynthesis of CTs is under genetic as well as environmental control [21], leading to a wide range of CTs differing in polymer length, B-ring hydroxylation and stereochemistry [22]. Variation in CT structure is associated with altered protein-binding and oxidative properties [23], which can impact below-ground microbial communities in multiple ways. CTs can inhibit microbial growth via direct toxicity [24] or via complexation with substrates such as cell wall polysaccharides [25], nitrogenous compounds [16,26] and exoenzymes [27].
Three microbial functional groups play especially important roles in decomposition and nutrient cycling in temperate forests, viz. fungi, nitrogen-fixing bacteria and ammonia-oxidizing bacteria. Fungal communities are the most abundant microbes in soils of conifer forests and play an important role in litter decay and nutrient cycling [28]. A variety of saprophytic and mycorrhizal fungal species are able to decay litter using extracellular enzymes [29] and are biogeographically distributed according to climate, substrate quality, vegetation [30] and in some cases dispersal limitations [31,32]. Nitrogen is limited in temperate forests [33], but essential for meeting the stoichiometric requirements of microbes responsible for decay [34]. Free-living diazotrophs are important sources of N input into forest ecosystems [33,35] and can be responsible for exogenous N sources in litter with low N availability, thus influencing decay rates. By catalyzing the rate-limiting step for nitrification, ammonia-oxidizing bacteria act to increase NH3 availability [36]. Fungi are affected by CT amendments in diverse experimental systems [37]. Free-living diazotrophs are sensitive to environmental disturbance [38], and they also react to CT amendments in experimental microcosms [37].
In an effort to quantify the effects of litter chemistry and climate on litter decomposition and N mineralization, we designed a field study that used bags of leaf litter with different levels of CTs and N to assess the impact of these variables on in situ decay rates and associated microbes in the forest floor over a period of 43 months. This study spanned multiple field sites located along a temperature and moisture climatic gradient, each site composed of similar, typical over- and understory species dominated by mature Douglas-firs. Shifts in plant species along climate gradients can influence microbial composition in soils and confound field analyses of microbial responses to climate [39]. In this study, the similarity of vegetation along the climate gradient was ideal for assessing the impacts of climate on decay-associated microbes in the field. Furthermore, previously established biogeographical assessment of microbial communities in the forest floor and mineral soil at these sites [40] provided a unique baseline from which to compare microbial responses observed on decaying leaf litter.
Our particular objectives were to assess how the structure (diversity, richness, evenness) and composition of microbial functional groups associated with decaying litter might be affected by litter chemistry and climate gradients. This included assessment of temporal and spatial components, since microbial communities are known to shift with time [41] and space as litter decay proceeds [42]. Polymerase chain reaction coupled with denaturing gradient gel electrophoresis (PCR-DGGE) was used as a cost-effective technique for microbial community fingerprinting. PCR-DGGE is sufficient for assessment of community structure and diversity when functional groups are targeted, allowing for basic screening and supporting focused analysis of treatment effects on abundant (>0.5% relative abundance) functional groups [43,44,45].
In the present study, we tested the hypotheses that (1) slower decay rates associated with high CT levels and low N concentrations in leaf litter are caused by inhibition of microbial colonization and diversity and that (2) initial litter chemistry continued to be an important determinant of microbial community composition over 43 months of litter decay.

2. Methods

2.1. Litter Bags and Field Sites

To assess decay impacts in forest floor soils, leaf litter samples having a range of N (6–15 mg g−1) and CT (21–278 mg g−1; Table 1) levels were placed in fine, polypropylene mesh (0.5 mm opening) litter bags (25 × 25 cm) for decay on the forest floor, as described in [19,46]. Using litter-fall traps, naturally abscised leaf litter from Populus angustifolia (poplar) and Pseudotsuga menziesii (fir = Douglas-fir) was obtained in 2008 from a common garden established in 1991 at the Ogden Nature Center (Ogden, UT, USA) and in 1991 from the Shawnigan Research Forest [46], respectively, and oven-dried (50 °C) prior to sealed storage in darkness. Litter bags were filled with either 5 g of bulked fir needles or 4.6 g of poplar leaf litter (oven-dried weight, 12 h at 70 °C). High and low CT treatments for poplar litter were made by bulking leaves from clones with known CT level, while ensuring other nutrients (e.g., N) and phenols (e.g., salicortin) were kept constant between treatments [47,48]. High-N litter treatments were prepared by spraying an atomized glutamine solution onto individual leaf-litter samples, while control (low-N) treatments were sprayed with an equal amount of distilled H2O. Determination of C, N, CT and proximate chemistry (e.g., AUR) in leaf litter prior to decay was outlined in [19,49].
Four strings of litter bags, each with six litter types (Table 1), were laid on the forest floor at four replicated microplots (ca. 2.5 × 2 m) in close proximity (<23 m) to meteorological stations at nine sites on southeastern Vancouver Island, BC, Canada [40]. Sites were located along south, central and north transects spanning the coastal Douglas-fir biogeoclimatic (BGC) zone (DF), coastal western hemlock biogeoclimatic zone (WH) and the transition between the two (TR) [50]. Moisture probes and thermocouples connected to CR10 (Campbell Scientific, Logan, UT, USA), HOBO (models U30 and U12, Onset Computer Corporation, Bourne, MA, USA) and Decagon (model Em50 Digital, Decagon Devices, Pullman, WA, USA) data loggers and Ibuttons (Maxim Integrated, San Jose, CA, USA) recorded hourly measurements of air and soil temperature (30 cm above and below ground, respectively), along with soil moisture at each site.

2.2. Sample Collection and Processing

Litter bags were collected after 7, 12, 24 and 43 months of decay in the field, and kept cool (4 °C) until processed in the laboratory within 36 h after removal. For microbial analysis, a representative subsample (0.25 g wet weight) of each litter sample was removed from litter bags using sterile forceps and gloves; the remaining litter fraction was used to assess % C mass remaining and N mineralization [19]. For each sampling period, four replicate subsamples of each litter treatment were pooled by site and homogenized manually using sterile technique. The pooled subsamples were stored at −20 °C prior to DNA extraction. For extraction, samples were first pulverized (1 min at 30 Hz) using a liquid N compatible bead beater (Retsch Mixer Mill MM400, Newtown, PA, USA). DNA was extracted from 0.25 g of pulverized tissue using PowerLyserTM PowerSoil DNA Isolation kits (according to manufacturer’s protocol and including a 10 min incubation period at 70 °C for optimal cell lysis). Extracted DNA was assessed for quantity (μg/mL) and quality (λ absorbance ratios 230/280 and 230/260) using an ND-1000 spectrophotometer (NanoDrop products, Thermo Fisher Scientific, Wilmington, DE, USA) to account for amplification bias that could occur during DNA extraction.
Primer sets for amplification of microbial functional groups each contained GC-clamps in order to prevent complete separation of complementary strands in PCR products during DGGE [40,51,52]. PCR-DGGE fingerprinting of fungi, ammonia-oxidizing β-proteobacteria and free-living diazotrophs was performed using primer sets 18S FF390/FR1, AmoA 1F’/2R, and NifH-universal ForA/ForB/Rev, respectively, according to modifications by [40] of existing protocols [53,54,55,56].
Pooled DNA from all samples in each sampling period was used as standard ladders to allow for band comparison within and between gels, and PCR runs (see Figure S1 from [40]). The migration distances of bands were determined using GeneTools software (Syngene, Frederick, MD, USA) and consolidated from at least two reference ladders per gel, followed by visual inspection. Operational taxonomic units (OTUs) obtained from PCR-DGGE gels were consolidated into a presence/absence matrix across time and space, while values for the sum of all peak heights for all bands in a lane (N), the total number of bands in a lane (S) and the peak height of the ith DGGE band in a lane (ni) were used for calculating diversity indices (Shannon’s diversity—Equation (1), richness—Equation (2) and Pielou’s evenness—Equation (3)) for each sample.
H = [ n i N ] l o g [ n i N ]
  d = [ S 1 ] l o g [ N ]
  J = H l o g [ S ]

2.3. Data Analysis

R Statistics software (version 3.3.1) was used for all statistical analyses [57]. We first modeled the multivariate response of DNA extraction efficiency (DNA concentration, 230/280 λ absorbance ratio and 230/260 λ absorbance ratio) to time, latitude, zone, litter type and their interactions, using constrained redundancy analyses (RDAs), followed by post hoc analysis of each trait using univariate regression. The effects of climate (Table 2), litter type (Table 1), time and space on microbial communities associated with decaying litter were determined with partially constrained RDA and correspondence analysis (CCA) [58,59] using the R package ‘vegan’ [60]. Climate, litter chemistry and time were each assessed as categorical (e.g., north DF, central DF and south DF) and continuous (e.g., mean, minimum and maximum annual air temperature) predictors (Table 1 and Table 2).
Categorical modeling of climatic effects was assessed by site (i.e., nine unique latitude and zone combinations), since latitude and zone did not always reflect clear climatic gradients. However, site categories represented both climatic and spatial effects, which were separated only in models using continuous predictors where spatial factors were determined by principal coordinates of neighbor matrices (PCNMs) [61]. Similar to methods first described to parse out spatial variation [58,59], microbial community responses were first assessed separately for time, climate, litter chemistry and PCNM effects, then jointly in partially constrained models, in order to determine the proportion of shared effects. Categorical models were also assessed after removing any significant effects of DNA extraction efficiency. The significance of models and individual predictors was assessed by permutation tests for RDA and CCA using reduced models.
Categorical models were reduced by step-wise elimination of least significant predictors. In continuous models, the number of climate and PCNM predictors (and litter chemistry predictors in the case of joint poplar and fir analyses) was reduced using step-wise model building according to AIC-like criteria [60]. This was followed by elimination of non-significant factors while assessing each separately. Candidate climate variables only included those with Pearson’s correlation r < 0.81 (Table 2). Elemental composition, proximate chemistry and CT variables were included as litter chemistry effects when assessing both poplar and fir data sets jointly, but only N and CT (poplar only) content were considered when assessing poplar and fir data sets separately. Post hoc analysis of litter chemistry effects in poplar litter was also performed on reduced data sets, removing the 43-month data. When considering all three targeted functional groups, best fit models were also confirmed using ‘leap’ model selection for multivariate responses [62]. In each case, models and predictors were deemed significant by permutation tests.
Responses in microbial community structure (Shannon’s diversity, richness and Pielou’s eveness) or composition (OTU presence or absence) of fungi, ammonia-oxidizing bacteria and nitrogen-fixing bacteria were assessed jointly and also separately for each functional group. Data sets for poplar and fir litter types were also treated jointly and separately. When modeling microbial structure, responses were assessed for multivariate normality, power transformed where appropriate and missing data (OTU absence causing null measurements) were imputed using the iterative PCA algorithm from the ‘missMDA’ R package in order to reduce the influence of such values in ordination analyses [63]. The effects of % C remaining in litter samples on overall microbial responses were also assessed with and without partial constraint by best-fit predictors. Results on the %C mass remaining in different litter treatments over 43 months have been previously published [19].

3. Results

Fingerprinting of microbial functional groups using PCR-DGGE identified 66 unique fungal OTUs, 23 NifH OTUs and 23 AmoA OTUs distributed among 192 samples.

3.1. DNA Quality

DNA quality from poplar litter was affected by time, site and litter type (p = 0.001, 0.013 and 0.009, respectively), specifically due to generally greater DNA concentration over time and in LTHN treatments, lower 260/280 λ ratios over time and in the SDF site and greater 260/230 ratios in 43-month samples (post hoc p < 0.05). DNA quality in fir litter was also affected by time (p = 0.001), due to a significant decrease in the 260/280 λ ratio over time of decay (p = 0.004). The 260/230 λ ratios in fir samples decaying for 12 and 24 months were on average half the values in samples at 7 and 43 months (post hoc p < 0.001). DNA concentrations increased in fir samples until 24 months of decay followed by a drop to 12-month levels (post hoc p < 0.019). On average, DNA concentration was 23.4% less in poplar than fir litter (p = 0.002), while the 260/280 λ ratio was 12.8% greater (p < 0.001). However, DNA quality generally did not confound significance of categorical models (Table 3). Only changes in fungal community composition associated with N treatments in fir litter, which were not significant when using continuous variables, could instead have been attributed to DNA concentration effects (Table 3B).

3.2. Temporal Effects

Time was significantly associated with up to 9% and 9.1% of the variation in overall microbial community structure and composition, respectively (Table 3—analyses of all taxa for both poplar and fir litter). No temporal effects were detected for the AmoA functional group, except for a 1.2% shift in AmoA community composition when considering time as a continuous variable (Figure 1A), despite lower explanatory power for overall microbial community structure and composition (4.3% and 5.7%, respectively; Figure 2) compared to categorical analyses. The structure of fungal communities in both poplar and fir litter changed over time (Table 3A and Table 4A). In poplar litter, fungal richness was greatest after 43 months of decay and was lowest during the 7-month sampling period (17 and 14 OTUs on average, respectively), while fungal community evenness generally decreased over time and diversity remained unchanged (p = 0.001). In fir litter, fungal richness decreased slightly until 24 months of decay (15–14 OTUs per fir sample, on average, in fir litter; p = 0.001), with concurrent increases in evenness and diversity, but was greatest in samples after 43 months (23 OTUs per fir sample, on average; p = 0.001).
The structure of NifH communities was only affected by time in poplar litter, not fir (Table 3A and Table 4A). The richness, evenness and diversity of NifH communities in poplar litter increased until 24 months of decay, on average being 17.6%, 22.3% and 33.5% greater after 24 months of decay compared to after 7 months of decay (p = 0.001). Composition of fungal and NifH communities also changed over time, with more pronounced differences found after the first year of decay (Figure 1A and Figure 3A). Overall microbial community composition was significantly associated with 3.1% of changes in the C remaining in litter (p = 0.001), which accounted for 1.3% of changes independent of time, litter chemistry, space or climate (p = 0.001). Overall community structure was not significantly associated with changes in C remaining in litter (p = 0.272), since the majority of associated effects (3.4 of 3.9%) were accounted for by other factors.

3.3. Litter Chemistry Effects

Litter type was significantly associated with microbial community structure (5.7%), but generally not when treating poplar and fir litter separately (Table 3A and Table 4A). Fungal richness, evenness and overall diversity were greater in fir than poplar litter (8.3%, 4.0% and 5.9%, respectively). Differences in fungal structure were primarily associated with greater diversity, evenness and, to a lesser extent, richness in litter with lower CT content, while greater AUR was associated with greater richness and less evenness (Table 4A). Greater ash content in fir litter, which was the chemical variable showing the least Z-score variance within litter species, was associated with the 18.5% greater NifH community evenness and 24.7% lower richness in this litter compared to poplar, as well as with different NifH community composition (Table 4A).
When assessing both poplar and fir litter, higher initial N content was associated with greater AmoA diversity and different AmoA community composition (Table 4A; Figure 1A,C). However, despite not detecting AmoA genes in any fir litter samples, N content in poplar litter was only deemed significant in determining AmoA community structure and composition when omitting samples collected after 43 months of decay (i.e., only samples decaying for 7, 12 and 24 months; p = 0.037 and 0.034, respectively; Table 4).
Overall, litter type was significantly associated with 4.1% of detected variance in microbial community composition (categorical analyses); 3.2% could be attributed to initial N, CT and AUR content in litter (Figure 1C and Figure 2B). Condensed tannin and AUR content accounted for 3% of differences in fungal composition among litter types. All but one of the OTUs were detected poplar samples, while only 60 out of 66 fungal OTUs and 18 out of 23 NifH OTUs were detected in fir litter. When treating litter species and functional groups separately, litter chemistry was only significant in describing fungal responses in poplar litter (Table 3B). CT and N treatments were associated with shifts in fungal community composition in poplar litter (Figure 1D), however, this was only significant when omitting samples collected after 43 months of decay (i.e., only samples decaying for 7, 12 and 24 months; p = 0.039 and 0.014, respectively; Table 4B). Litter chemistry treatments still accounted for 2.2% and 2.4% of overall changes in fungal OTUs when solely considering poplar (Figure 3D) or fir data sets, respectively (Table 3B); these effects were not distinctly clustered when models included all functional groups and litter species (Figure 3C).

3.4. Spatial and Climatic Effects

Site (latitude × zone) accounted for 11.2% and 7.3% of overall microbial community structure and composition, respectively (p = 0.001, Figure 3B, Table 3). The structure of each fungal and NifH community differed between sites in poplar litter (p ≤ 0.022), while only differences in fungal community structure were detected fir litter (p = 0.047, Table 3A). All fungal diversity indices were greater at northern sites, especially in the TR and WH zones, while lowest at southern sites. Diversity indices of NifH communities in poplar litter were greatest in central sites, especially in the DF zone. Incidence of AmoA OTUs was less frequent than for other functional groups, however, community structure was still associated with site in poplar litter (p = 0.001, Table 3A). Over time, the 23 AmoA OTUs we detected were largely found on poplar litter in NWH and DF sites. Reflecting structural indices, cluster analysis of microbial composition revealed generally different communities in litter decaying at NWH or CDF sites vs. the other remaining sites. Communities in the other sites divided into two groups dominated by northern or southern and central sites (Figure 3B).
Climate variables did not significantly describe microbial community structure in fir litter (p > 0.05; Table 4A). Despite accounting for 9.1% of microbial community structure in poplar litter, 6.8% of climatic effects on community structure were shared with spatial effects (Figure 2A), rendering independent climatic effects insignificant (Table 4A). Nonetheless, 3.5% lower MSmin, which was inversely correlated with PCNM 3 (r = −0.514), was associated with 13% greater diversity and 27% greater richness of the fungal community, while 4.9% lower MS, which correlated with PCNM 2 (r = 0.619), was associated with 5.8 times greater diversity and 4.6 times greater richness of the AmoA community (Table 4A). Increasing TA (1.6 °C) correlated with PCNM 3 (r = 0.575) and was associated with lower NifH community diversity, richness and evenness (Table 4A). The slightly lower AmoA community richness was also associated with greater TA (Table 4A). It was confounded by the joint effects of PCNMs 1, 2 and 3 (Table 4A), although spatio-climatic correlations were low (r < 0.316).
Climate accounted for 2.7% of microbial community composition independently, and shared a further 3.2% of effects with PCNM (Figure 2B). Independent climatic effects significantly accounted for community composition for each functional group and litter species subset, except for fungi in poplar litter, where the effects of MSmin and MSmax were confounded by negative correlations with PCNMs 3 (r = −0.514) and 1 (r = −0.666), respectively (Table 4B). MSmin, TA, TSrange and PET were all significantly associated with NifH community composition in poplar litter, while the effects of TAmin were confounded by correlations with PCNM 4 (r = 0.752; Table 4B). Nonetheless, TAmin was significantly associated with the composition of AmoA communities in poplar litter and NifH communities in fir litter (Table 4B). TSrange was significantly associated with NifH community composition in poplar litter and fungal community composition in fir litter (Table 4B). The significant effects of TAmax on AmoA community composition in poplar (p = 0.049) and composition of the fungal community in fir litter (p = 0.003) were only detected when sub-setting litter species and functional groups (Table 4B).
Spatial PCNM components 1, 3, 4 and 5 were significantly associated with overall microbial community structure (9.4%, p = 0.001) while PCNMs 1, 3, 5 and 6 were associated with composition (4.5%, p = 0.001; Figure 4). PCNM 3 was the only spatial component associated with microbial composition in fir litter, namely due to correlations with NifH community composition in each litter species (p < 0.008), but was confounded by climate effects (Table 4B and Figure 2B). Although PCNM 4 was significantly associated with overall microbial community structure in fir litter (p = 0.01), most likely due to fungi per joint analyses of poplar and fir samples (Table 4A), no spatial patterns were significantly attributed to the structure of individual functional groups. PCNM 4 was also associated with the composition of NifH communities in poplar litter (p = 0.045), but correlations with TAmin (r = 0.752) confounded such effects (Table 4B). In poplar litter, PCNM 1 was associated with the composition of each functional group (p < 0.002) and the structure of fungal and AmoA communities (p = 0.001 and 0.002, respectively).
Despite being correlated with less MSmax and greater TAmax (r = −0.666 and 0.661, respectively), PCNM 1 was only confounded by climate when assessing NifH community composition in poplar litter (p = 0.139). PCNM 3 was associated with the composition and structure of only fungal and NifH communities (p < 0.018), but responses were generally confounded by climate. PCNM 2 was uniquely associated with the structure of AmoA communities (p = 0.014), although confounded by correlations with MS (r = 0.619), while PCNM 5 was associated with both structure and composition (p = 0.001). PCNM 5 was also significant in describing overall NifH community composition (p = 0.039), but not when poplar and fir data were treated separately nor when accounting for climate. PCNM 6 was only associated with fungal community composition in poplar litter (p = 0.017).

4. Discussion

Litter chemistry and weather affected the composition of fungi and N-cycling communities in decaying litter, often reflecting trends in C decay and nutrient cycling. Spatial layout and temporal effects, which could not be completely isolated from climate effects, further explained significant portions of the structure and composition of microbial communities. Analyzing subsets of microbial communities in decaying leaf litter sometimes revealed more subtle (or incongruent) effects that were masked when all community data were used (e.g., the effects of TAmax on AmoA and fungal community structure), similar to observations made of baseline soil communities [40]. Carefully designed cluster analyses of constrained ordinations, accompanied by independent testing of explanatory variables (within potentially nested effects), revealed more than would otherwise be apparent (e.g., Figure 3C,D). We were able to observe macroecological trends impacting soil microbes, such as CT and N effects nested within time and litter species.

4.1. Ammonifying Communities Are Altered by Weather and Litter Nitrogen

Temporal shifts in the composition of the AmoA community progressed gradually over time, and hence could not be detected for this rarer functional group using only pairwise-like comparisons (Table 3). The increasing detection of AmoA functional genes in decaying litter with high N and at warmer sites with ample soil moisture is indicative of increasing ammonia availability in such rapidly decaying litter, albeit at earlier and later time frames, respectively [19]. Similarly, warmer temperatures have been associated with the biogeography of AmoA community composition in the forest floors at our sites [63], indicating that the influence of climate on AmoA communities in decaying litter is reflective of long-term trends at our sites. The absence of AmoA genes in fir litter decaying side-by-side with poplar litter (where they were detected) confirms that site-specific inoculation is not the reason for observed effects, nor is the endemicity of the litter.
Increased air temperature may alter the N forms available to trees, favoring species with affinities for NO3, unless the annual additive short-term effects of reduced litter N can offset the longer-term consequences of temperature. In the presence of a favorable future precipitation and humidity regime, changes in nutrient status in warmer climates could favor the growth of western hemlock, which preferentially takes up NO3 in contrast to fir [64,65]. This would be inconsistent with a hypothesized westward shift of the transition zone between coastal western hemlock (CWH) and coastal Douglas-fir (CDF) BGC zones [66].
The rapid decay rate of fir litter one year after the onset of decomposition [19] and concurrent increasing N mineralization in this litter type (unpublished) were never accompanied by the emergence of AmoA populations (i.e., never detected). We expect that mycorrhizae translocated surplus N from the in situ fir litter type to nearby trees or plants, preventing ammonia accumulation at the site of decay.

4.2. Weather Can Account for Early Succession in Fungal Composition

The chemical, climatic and PCNM variables selected in our final models thoroughly described the range of effects in our experimental design, since less than 1% of the effects of litter chemistry or site were not accounted for by continuous variables describing such treatments. The temporal effects on microbial community structure and composition not accounted for when using continuous time variables (4.7% and 3.4%, respectively) indicated non-linear trends in fungal and NifH community succession during decay. Successional shifts in fungal community composition were more pronounced after 24 months, and especially after 43 months of decay (Figure 3A), and did not affect fungal diversity in our decaying poplar and fir litter, similar to rapid successional changes observed in oak systems [41]. Increasing time periods between successive sampling could be responsible for the increasing divergence of community composition between sampling periods, and potentially indicate distinct annual shifts in community composition.
In a separate analysis, we were able to attribute a significant portion of the succession in fungal communities over the first 24 months of decay to weather regimes, when adjusting temperature and moisture variables over time in order to only reflect the weather litter bags were exposed to (see Supplementary Material for details). The ability of weather to explain successional changes in microbial community structure and composition runs counter to the classical view of microbial succession as a response to changing resource availability. Rapid turnover of microbial communities has been observed in response to climatic stress [67]. Our observations of microbial succession as a response to ongoing weather variation would support such findings, although confounded by exposure time and starting inoculum. Since climate impacts microbial biogeography [40], it is unclear if weather is directly accountable for observed fungal succession, or if instead climate is indirectly affecting succession via selection of available microbial inoculum in soils. Lin et al. [68] studied disturbance and fresh litter deposition in the forests of another conifer, Cryptomeria japonica, and showed that divergent fungal communities in different plots essentially converged and became identical after 21 months. Given the 43-month duration of our experiment, we therefore expect that the differences we observed in our communities were ultimately related to the impact of climate on community succession in the litter bags and/or its impact on the composition of the parent community. Experimental protocols using high-throughput sequencing for identifying microbial communities would allow direct comparison of communities in soils and decaying litter. Our results warrant further investigations into the extent of climate-induced microbial succession.

4.3. Short-Term Effects on Fungal Communities Driven by Litter Chemistry

Climate-induced shifts in foliar CT and N content [11,12,13,14] would affect litter fall chemistry and thus could alter below-ground fungal communities. The ranges of CT content (76–278 mg g−1) and N content (7.5–15 mg g−1) within poplar litter affected fungal composition up until 24 months of decay, and yet were not associated with the structure of fungal communities. These observations contradict our predictions that CTs inhibit overall microbial diversity. In this case, higher CT levels were simply associated with different fungal communities; we presume that a different cohort of decay fungi was selected from the more diverse background community. The effects of litter chemistry on community structure were mostly driven by litter species; parameters distinguishing poplar from fir litter (e.g., CT form and cutin concentration) were more important than intra-specific variation in CT and N concentrations in determining community diversity indexes. Intra-specific variation in CT and N concentrations may not be broad enough to elucidate relationships between litter chemistry and microbial community structure (i.e., diversity, richness and evenness) or alternatively the cohorts of decay fungi available to decompose litter are diverse enough that CT levels do not matter.
The changes in fungal composition associated with greater CT and lower N content were accompanied by reduced rates of decay [19]. The resolution of PCR-DGGE allowed linkage of changes in microbial community composition to C loss. The lack of association between our diversity indices (i.e., community structure) and litter-decay status supports previous findings using different microbial screening tools [69].
The different microbial communities detected on decaying poplar and fir litter mirror their differing rates of decay [19] and are likely indicative of distinct nutrient flows. It is possible that some of the fungal OTUs affected by our CT treatments were mycorrhizal species, as was the case in a study of the effects of CT from black spruce and sheep laurel [70]. Mycorrhizae strongly influence competition for nutrients between plants and decomposers, as well as rates of C sequestration [71]. Feedback loops between above- and below-ground communities were also demonstrated in grass/shrub systems, where plant community traits (e.g., C/N ratios) reflect their fundamental niche and generate positive feedback loops via microbial selection that can perpetuate the nutrient status of an ecosystem [72].

4.4. Diazotrophs Differ among Litter Types

Litter species had a strong influence on NifH communities and how they were affected by climate. Different NifH OTUs colonized decaying litter from poplar and fir, hence some of the detected diazotroph species may be species- or litter-specific. Winder et al. [37] identified several diazotrophic endophytes of poplar during microcosm incubations of poplar litter, a reminder that diazotrophs in our litter bags may not have been limited to free-living species. NifH community composition was further affected by the annual range in soil temperature in poplar litter, and by the minimum annual air temperature in fir litter (Table 4B). Apart from being a potential source of essential mineral elements (e.g., Mg, Fe, Mo), it is unclear how slight differences in initial ash content could be responsible for differing NifH community structure or composition between poplar and fir litter. Since ash was correlated with CT differences between litter species (r = −0.763), it is possible that the amount or type of CT, even in low-CT poplar treatments, may have selected for different diazotrophic species. Alternatively, ash content may correlate with other compounds or metabolites that we did not measure and vary in content between fir and poplar litters (e.g., waxes and cutins). Selective colonization may also be simply the result of significantly different C-decay trajectories and nutrient availability between litter species during decay, as suggested by Widmer et al. [54]. We generally detected more NifH genes (i.e., more richness) in poplar litter compared to fir litter. The slower decay rates for poplar litter from 7 to 43 months of decay compared to fir litter could suggest low N availability in this litter during this time period, however, the presence of genes from specific microbial communities (such as N-cycling soil bacteria) may not necessarily directly reflect ecological activity [73].

5. Conclusions

High CT and low N concentrations in leaf litter were associated with changes in microbial community composition, especially fungi, which may be directly responsible for slower early (<2 years) decay rates of poplar litter and lead to a net annual sequestration of litter carbon in forest soils. Contrary to our hypothesis, high CT and low N treatments did not inhibit microbial colonization or diversity. The effects of initial litter chemistry, AUR (35.3–41.5% dw), CT (2.1–29.1% dw) and N (0.61–1.57% dw) on microbial community composition in decaying litter were less than the joint effects of 1.6 and 3.4 °C differences in annual mean and minimum air temperature, respectively, and 4.1% and 26.9% differences in annual minimum and maximum soil moisture, respectively (Figure 2). Long-term microbial composition in forest soils may therefore be primarily affected by future climates, especially since the effects of litter chemistry diminish over time as the compositions of different types of litter converge. Furthermore, shared spatio-climatic and temporal–climatic effects made up a considerable portion of climate effects, and warrant further investigation in the face of climate change. The ability of available microbial inoculum to colonize both in situ and ex situ litter types shows promise for the capacity of below-ground systems to adapt to large-scale migration (natural or assisted) of trees in response to climate. However, differences in community assemblages were indicative of functional shifts that could compound over time and alter decay rates as climate change progresses. These differences should be carefully monitored where regions experience significant natural or deliberate shifts in forest cover. Baseline assessments of microbial biogeography in soils were helpful in interpreting our experimental results and will be of increasing importance in the analysis of future impacts of climate change on such a critical and lesser-known part of our ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/jof8070735/s1, Supplementary analysis of fungal communities associated with decay of Douglas-fir and poplar leaf litter over 24 months. Experimental design, analyses and model selection were performed as for data covering 43 months of data, however, continuous climate variables at each site were allowed to change over time in order to only reflect the time period litter was exposed to in the field, spatial effects (PCNM) were included in categorical modeling and additional predictors were included in modeling. Figure S1: Principal coordinates of neighbor matrices (PCNMs) of field-site geographic locations used in modeling. PCNMs 1, 3, 4 and 6 were associated with fungal community structure while only PCNMs 1 and 3 with composition; Figure S2: Fraction of (A) fungal community structure (diversity, richness, evenness) or (B) fungal community composition (OTU pres./abs.) explained by categorical or continuous predictors representing space, year, climate, litter chemistry or shared effects (shared), using RDA or CCA, respectively. Asterisks mark when conditioned models used to calculate a given fraction are significant (permutation test for RDA or CCA p < 0.05); Figure S3: CCA of fungal community composition from decaying leaf litter samples in response to litter type (categorical constraints) omitting the effects of years of decay, latitude, zone and space; Figure S4: CCA of fungal composition in response to: (A) litter chemistry alone (ratio of insoluble CT to N (insol.N), C/N ratio (C.N) and % C (C); 3.7% of total variance), (B) climate alone (maximum hourly and minimum mean daily soil moisture (MSmax, dMSmin), mean daily and minimum mean daily air temperature (dTA, dTAmin), and frost-free days (FFDs); 7.2% of total variance) and (C) both chemistry and climate (11% of total variance).

Author Contributions

Conceptualization, P.-E.S., J.A.T., R.S.W.; Methodology, P.-E.S., R.S.W., J.A.T.; Software, P.-E.S.; Formal Analysis, P.-E.S.; Investigation, P.-E.S., R.S.W., J.A.T.; Resources, R.S.W., J.A.T.; Data Curation, P.-E.S.; Writing—Original Draft Preparation, P.-E.S., J.A.T., R.S.W.; Writing—Review and Editing, P.-E.S., R.S.W., J.A.T., C.P.C.; Visualization, P.-E.S., J.A.T.; Supervision, C.P.C., J.A.T.; Project Administration, J.A.T.; Funding Acquisition, C.P.C., J.A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Sciences and Engineering Research Council (NSERC) of Canada—U. Victoria CREATE Program in Forests and Climate Change as a collaboration between the University of Victoria and the Pacific Forestry Centre (PFC) of Natural Resources Canada (NRCan). The APC was funded by the NRCan Canadian Forest Service—Climate Change Program and an NSERC Discovery Grant awarded to C. Peter Constabel.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw data are available on request to the lead author.

Acknowledgments

Special thanks to Nicholas von Wittgenstein (Centre for Forest Biology student intern) for putting together litter bags, Heather Klassen (BC Forest Service) for all the meteorological monitoring, Dave Dunn, Grace Ross and Rebecca Dixon (PFC Chemical Services) for proximate chemical analyses and help in the lab and Réal Roy (University of Victoria) and Doug Maynard (Canadian Forest Service) for their guidance and support. In particular, we acknowledge Thomas G. Whitham and the Cottonwood Research Group at Northern Arizona University for having provided us with poplar litter material that made this work possible, along with chemical and genetic information about the source plants.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schlesinger, W.H.; Andrews, J.A. Soil respiration and the global carbon cycle. Biogeochemistry 2000, 48, 7–20. [Google Scholar] [CrossRef]
  2. Palosuo, T.; Liski, J.; Trofymow, J.A.; Titus, B.D. Litter decomposition affected by climate and litter quality—Testing the Yasso model with litterbag data from the Canadian intersite decomposition experiment. Ecol. Mod. 2005, 189, 183–198. [Google Scholar] [CrossRef]
  3. Moore, T.R.; Trofymow, J.A.; Prescott, C.E.; Fyles, J.; Titus, B.D. Patterns of carbon, nitrogen and phosphorus dynamics in decomposing foliar litter in Canadian forests. Ecosystems 2006, 9, 46–62. [Google Scholar] [CrossRef]
  4. Prescott, C.E. Litter decomposition: What controls it and how can we alter it to sequester more carbon in forest soils? Biogeochemistry 2010, 101, 133–149. [Google Scholar] [CrossRef]
  5. Solomon, S.; Hegerl, G.; Heimann, M.; Hewitson, B.; Hoskins, B.; Joos, F.; Jouzel, J.; Kattsov, V.; Lohmann, U.; Matsuno, T.; et al. Contribution of Working Group 1, Fourth Assessment Report of the IPCC; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  6. Raich, J.W.; Potter, C.S. Global patterns of carbon-dioxide emissions from soils. Glob. Biogeochem. Cycles 1995, 9, 23–36. [Google Scholar] [CrossRef] [Green Version]
  7. Van Gestel, M.; Merckx, R.; Vlassak, K. Microbial biomass responses to soil drying and rewetting—The fate of fast-growing and slow-growing microorganisms in soils from different climates. Soil Biol. Biochem. 1993, 25, 109–123. [Google Scholar] [CrossRef]
  8. de Vries, F.T.; Shade, A. Controls on soil microbial community stability under climate change. Front. Microbiol. 2013, 4, 265. [Google Scholar] [CrossRef] [Green Version]
  9. Kurz, W.A.; Dymond, C.C.; White, T.M.; Stinson, G.; Shaw, C.H.; Rampley, G.J.; Smyth, C.; Simpson, B.N.; Neilson, E.T.; Trofymow, J.A.; et al. CBM-CFS3: A model of carbon-dynamics in forestry and land-use change implementing IPCC standards. Ecol. Mod. 2009, 220, 480–504. [Google Scholar] [CrossRef]
  10. Tarnocai, C.; Canadell, J.G.; Schuur, E.A.G.; Kuhry, P.; Mazhitova, G.; Zimov, S. Soil organic carbon pools in the northern circumpolar permafrost region. Glob. Biogeochem. Cycles 2009, 23, 1–11. [Google Scholar] [CrossRef]
  11. Mansfield, J.L.; Curtis, P.S.; Zak, D.R.; Pregitzer, K.S. Genotypic variation for condensed tannin production in trembling aspen (Populus tremuloides, Salicaceae) under elevated CO2 and in high- and low-fertility soil. Amer. J. Bot. 1999, 86, 1154–1159. [Google Scholar] [CrossRef] [Green Version]
  12. Moore, T.R.; Trofymow, J.A.; Taylor, B.; Prescott, C.; Camire, C.; Duschene, L.; Fyles, J.; Kozak, L.; Kranabetter, M.; Morrison, I.; et al. Rates of litter decomposition in Canadian forests. Glob. Chang. Biol. 1999, 5, 75–82. [Google Scholar] [CrossRef]
  13. Norby, R.J.; Cotrufo, M.F.; Ineson, P.; O’Neill, E.G.; Canadell, J.G. Elevated CO2, litter chemistry, and decomposition: A synthesis. Oecologia 2001, 127, 153–165. [Google Scholar] [CrossRef] [PubMed]
  14. Parsons, W.F.J.; Bockheim, J.G.; Lindroth, R.L. Independent, interactive, and species-specific responses of leaf litter decomposition to elevated CO2 and O3 in a northern hardwood forest. Ecosystems 2008, 11, 505–519. [Google Scholar] [CrossRef]
  15. Driebe, E.M.; Whitham, T.G. Cottonwood hybridization affects tannin and nitrogen content of leaf litter and alters decomposition. Oecologia 2000, 123, 99–107. [Google Scholar] [CrossRef]
  16. Hattenschwiler, S.; Vitousek, P.M. The role of polyphenols in terrestrial ecosystem nutrient cycling. Trends Ecol. Evol. 2000, 15, 238–243. [Google Scholar] [CrossRef]
  17. Hattenschwiler, S.; Hagerman, A.E.; Vitousek, P.M. Polyphenols in litter from tropical montane forests across a wide range in soil fertility. Biogeochemistry 2003, 64, 129–148. [Google Scholar] [CrossRef]
  18. Smyth, C.E.; Trofymow, J.A.; Kurz, W.A.; CIDET Working Group. Decreasing Uncertainty in CBM-CFS3 Estimates of Forest Soil C Sources and Sinks through Use of Long-Term Data from the Canadian Intersite Decomposition Experiment; Inf. Rep. BC-X-422; Natural Resources Canada, Canadian Forest Service: Victoria, BC, Canada, 2010. [Google Scholar]
  19. Shay, P.E.; Constabel, C.P.; Trofymow, J.A. Evidence for the role and fate of water-insoluble condensed tannins in the short-term reduction of carbon loss during litter decay. Biogeochemistry 2018, 137, 127–141. [Google Scholar] [CrossRef] [Green Version]
  20. Kraus, T.E.C.; Dahlgren, R.A.; Zasoski, R.J. Tannins in nutrient dynamics of forest ecosystems—A review. Plant Soil 2003, 256, 41–66. [Google Scholar] [CrossRef]
  21. Constabel, C.P. Molecular controls of proanthocyanidin synthesis and structure: Prospects for genetic engineering in crop plants. J. Agric. Food Chem. 2018, 66, 9882–9888. [Google Scholar] [CrossRef]
  22. Barbehenn, R.V.; Constabel, C.P. Tannins in plant-herbivore interactions. Phytochemistry 2011, 72, 1551–1565. [Google Scholar] [CrossRef]
  23. Hagerman, A.E.; Riedl, K.M.; Jones, G.A.; Sovik, K.N.; Ritchard, N.T.; Hartzfeld, P.W.; Riechel, T.L. High molecular weight plant polyphenolics (tannins) as biological antioxidants. J. Agric. Food Chem. 1998, 46, 1887–1892. [Google Scholar] [CrossRef] [PubMed]
  24. Scalbert, A. Antimicrobial properties of tannins. Phytochemistry 1991, 30, 3875–3883. [Google Scholar] [CrossRef]
  25. Tarascou, I.; Souquet, J.M.; Mazauric, J.P.; Carrillo, S.; Coq, S.; Canon, F.; Fulcrand, H.; Cheynier, V. The hidden face of food phenolic composition. Arch. Biochem. Biophys. 2010, 501, 16–22. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, L.H.; Lin, Y.M.; Ye, G.F.; Liu, X.W.; Lin, G.H. Changes in the N and P concentrations, N: P ratios, and tannin content in Casuarina equisetifolia branchlets during development and senescence. J. For. Res. 2008, 13, 302–311. [Google Scholar] [CrossRef] [Green Version]
  27. Joanisse, G.D.; Bradley, R.L.; Preston, C.M.; Munson, A.D. Soil enzyme inhibition by condensed litter tannins may drive ecosystem structure and processes: The case of Kalmia angustifolia. New Phytol. 2007, 175, 535–546. [Google Scholar] [CrossRef] [PubMed]
  28. Strickland, M.S.; Rousk, J. Considering fungal:bacterial dominance in soils—Methods, controls, and ecosystem implications. Soil Biol. Biochem. 2010, 42, 1385–1395. [Google Scholar] [CrossRef]
  29. Sinsabaugh, R.L.; Carreiro, M.M.; Repert, D.A. Allocation of extracellular enzymatic activity in relation to litter composition, N deposition, and mass loss. Biogeochemistry 2002, 60, 1–24. [Google Scholar] [CrossRef]
  30. Tedersoo, L.; Bahram, M.; Polme, S.; Koljalg, U.; Yorou, N.S.; Wijesundera, R.; Ruiz, L.V.; Vasco-Palacios, A.M.; Thu, P.Q.; Suija, A.; et al. Global diversity and geography of soil fungi. Science 2014, 346, 1078. [Google Scholar] [CrossRef] [Green Version]
  31. Peay, K.G.; Schubert, M.G.; Nguyen, N.H.; Bruns, T.D. Measuring ectomycorrhizal fungal dispersal: Macroecological patterns driven by microscopic propagules. Mol. Ecol. 2012, 21, 4122–4136. [Google Scholar] [CrossRef]
  32. Peay, K.G.; Bruns, T.D. Spore dispersal of basidiomycete fungi at the landscape scale is driven by stochastic and deterministic processes and generates variability in plant-fungal interactions. New Phytol. 2014, 204, 180–191. [Google Scholar] [CrossRef]
  33. Vitousek, P.M.; Howarth, R.W. Nitrogen limitation on land and in the sea—How can it occur. Biogeochemistry 1991, 13, 87–115. [Google Scholar] [CrossRef]
  34. Manzoni, S.; Trofymow, J.A.; Jackson, R.B.; Porporato, A. Stoichiometric controls on carbon, nitrogen, and phosphorus dynamics in decomposing litter. Ecol. Monogr. 2010, 80, 89–106. [Google Scholar] [CrossRef]
  35. Glick, B.R.; Patten, C.L.; Holguin, G.; Penrose, D.M. Biochemical and Genetic Mechanisms used by Plant Growth Promoting Bacteria; Imperial College Press: London, UK, 1991. [Google Scholar] [CrossRef]
  36. Barnard, R.; Leadley, P.W.; Hungate, B.A. Global change, nitrification, and denitrification: A review. Glob. Biogeochem. Cycles 2005, 19, 13. [Google Scholar] [CrossRef]
  37. Winder, R.S.; Lamarche, J.; Constabel, C.P.; Hamelin, R.C. The effects of high-tannin leaf litter from transgenic poplars on microbial communities in microcosm soils. Front. Microbiol. 2013, 4, 290. [Google Scholar] [CrossRef] [Green Version]
  38. Levy-Booth, D.J.; Winder, R.S. Quantification of nitrogen reductase and nitrate reductase gene abundance in soil of thinned and clear-cut Douglas-fir stands using real-time PCR. Appl. Environ. Microbiol. 2010, 76, 7116–7125. [Google Scholar] [CrossRef] [Green Version]
  39. Saetre, P.; Baath, E. Spatial variation and patterns of soil microbial community structure in a mixed spruce-birch stand. Soil Biol. Biochem. 2000, 32, 909–917. [Google Scholar] [CrossRef]
  40. Shay, P.E.; Winder, R.S.; Trofymow, J.A. Nutrient-cycling microbes in coastal Douglas-fir forests: Regional-scale correlation between communities, in situ climate, and other factors. Front. Microbiol. 2015, 6, 1097. [Google Scholar] [CrossRef] [Green Version]
  41. Voriskova, J.; Baldrian, P. Fungal community on decomposing leaf litter undergoes rapid successional changes. Int. Soc. Microbiol. Ecol. J. 2013, 7, 477–486. [Google Scholar]
  42. Ramette, A.; Tiedje, J.M. Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc. Nat. Acad. Sci. USA 2007, 104, 2761–2766. [Google Scholar] [CrossRef] [Green Version]
  43. Fromin, N.; Hamelin, J.; Tarnawski, S.; Roesti, D.; Jourdain-Miserez, K.; Forestier, N.; Teyssier-Cuvelle, S.; Gillet, F.; Aragno, M.; Rossi, P. Statistical analysis of denaturing gel electrophoresis (DGE) fingerprinting patterns. Environ. Microbiol. 2002, 204, 634–643. [Google Scholar]
  44. Fierer, N.; Jackson, R.B. The diversity and biogeography of soil bacterial communities. Proc. Nat. Acad. Sci. USA 2006, 103, 626–631. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Monroy, F.; van der Putten, W.H.; Yergeau, E.; Mortimer, S.R.; Duyts, H.; Bezemer, T.M. Community patterns of soil bacteria and nematodes in relation to geographic distance. Soil Biol. Biochem. 2012, 45, 1–7. [Google Scholar] [CrossRef]
  46. Trofymow, J.A.; CIDET Working Group. The Canadian Intersite Decomposition Experiment (CIDET): Project and Site Establishment Report; Inf. Rep. BC-X-378; Natural Resources Canada, Canadian Forest Service: Victoria, BC, Canada, 1998; 126p. [Google Scholar]
  47. Lamit, L.J.; Bowker, M.A.; Holeski, L.M.; Naesborg, R.R.; Wooley, S.C.; Zinkgraf, M.; Lindroth, R.L.; Whithama, T.G.; Gehring, C.A. Genetically-based trait variation within a foundation tree species influences a dominant bark lichen. Fungal Ecol. 2011, 4, 103–109. [Google Scholar] [CrossRef]
  48. Rehill, B.J.; Whitham, T.G.; Martinsen, G.D.; Schweitzer, J.A.; Bailey, J.K.; Lindroth, R.L. Developmental trajectories in cottonwood phytochemistry. J. Chem. Ecol. 2006, 32, 2269–2285. [Google Scholar] [CrossRef] [Green Version]
  49. Shay, P.E.; Trofymow, J.A.; Constabel, C.P. An improved butanol-HCl assay for quantification of water-soluble, acetone:methanol-soluble, and insoluble proanthocyanidins (condensed tannins). Plant Methods 2017, 13, 63. [Google Scholar] [CrossRef] [Green Version]
  50. Pojar, J.; Meidinger, D. British Columbia: The Environmental Settings. In Ecosystems of British Columbia; Meidinger, D.V., Pojar, J., Eds.; Special Report Series 6; BC Ministry of Forests: Victoria, BC, Canada, 2011; pp. 40–66. [Google Scholar]
  51. Sheffield, V.C.; Cox, D.R.; Lerman, L.S.; Myers, R.M. Attachment of a 40-base pair G +C-rich sequence (GC-clamp) to genomic DNA fragments by the polymerase chain reaction results in improved detection of single-base changes. Proc. Nat. Acad. Sci. USA 1987, 86, 232–236. [Google Scholar] [CrossRef] [Green Version]
  52. Kowalchuk, G.A.; Stephen, J.R.; DeBoer, W.; Prosser, J.I.; Embley, T.M.; Woldendorp, J.W. Analysis of ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria in coastal sand dunes by denaturing gradient gel electrophoresis and sequencing of PCR-amplified 16S ribosomal DNA fragments. Appl. Environ. Microbiol. 1997, 63, 1489–1497. [Google Scholar] [CrossRef] [Green Version]
  53. Vainio, E.J.; Hantula, J. Direct analysis of wood-inhabiting fungi using denaturing gradient gel electrophoresis of amplified ribosomal DNA. Mycol. Res. 2000, 104, 927–936. [Google Scholar] [CrossRef]
  54. Widmer, F.; Shaffer, B.T.; Porteous, L.A.; Seidler, R.J. Analysis of nifH gene pool complexity in soil and litter at a Douglas fir forest site in the Oregon Cascade Mountain Range. Appl. Environ. Microbiol. 1999, 65, 374–380. [Google Scholar] [CrossRef] [Green Version]
  55. Burgmann, H.; Meier, S.; Bunge, M.; Widmer, F.; Zeyer, J. Effects of model root exudates on structure and activity of a soil diazotroph community. Environ. Microbiol. 2005, 7, 1711–1724. [Google Scholar] [CrossRef]
  56. Nicolaisen, M.H.; Ramsing, N.B. Denaturing gradient gel electrophoresis (DGGE) approaches to study the diversity of ammonia-oxidizing bacteria. J. Microbiol. Methods 2002, 50, 189–203. [Google Scholar] [CrossRef]
  57. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  58. Borcard, D.; Legendre, P.; Drapeau, P. Partialling out the spatial component of ecological variation. Ecology 1992, 73, 1045–1055. [Google Scholar] [CrossRef] [Green Version]
  59. Legendre, P.; Borcard, D.; Peres-Neto, P.R. Analyzing beta diversity: Partitioning the spatial variation of community composition data. Ecol. Monogr. 2005, 75, 435–450. [Google Scholar] [CrossRef]
  60. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H. Vegan: Community Ecology Package; R Package Version 2.3-3.; 2016. Available online: http://CRAN.R-project.org/package=vegan (accessed on 15 December 2016).
  61. Borcard, D.; Legendre, P.; Avois-Jacquet, C.; Tuomisto, H. Dissecting the spatial structure of ecological data at multiple scales. Ecology 2004, 85, 1826–1832. [Google Scholar] [CrossRef] [Green Version]
  62. Cerdeira, J.O.; Silva, P.D.; Cadima, J.; Minhoto, M. Subselect: Selecting Variable Subsets; R Package Version 0.12-5; 2015. Available online: http://CRAN.R-project.org/package=subselect (accessed on 1 November 2015).
  63. Josse, J.; Husson, F. Handling missing values in exploratory multivariate data analysis methods. J. Soc. Fr. Stat. 2013, 153, 79–99. [Google Scholar]
  64. Gijsman, A.J. Soil-water content as a key factor determining the source of nitrogen (NH4+ or NO3) absorbed by Douglas fir (Pseudotsuga menziesii) and the pattern of rhizosphere pH along its roots. Can. J. For. Res. 1991, 21, 616–625. [Google Scholar] [CrossRef]
  65. Turner, D.P.; Sollins, P.; Leuking, M.; Rudd, N. Availability and uptake of inorganic nitrogen in a mixed old-growth coniferous forest. Plant Soil 1993, 148, 163–174. [Google Scholar] [CrossRef]
  66. Hamann, A.; Wang, T.L. Potential effects of climate change on ecosystem and tree species distribution in British Columbia. Ecology 2006, 87, 2773–2786. [Google Scholar] [CrossRef] [Green Version]
  67. Schimel, J.P.; Balser, T.C.; Wallenstein, M. Microbial stress-response physiology and its implications for ecosystem function. Ecology 2007, 88, 1386–1394. [Google Scholar] [CrossRef]
  68. Lin, W.-R.; Wang, P.-H.; Chen, W.-C.; Lai, C.-M.; Winder, R.S. Responses of soil fungal populations and communities to the thinning of Cryptomeria japonica forests. Microbes Environ. 2016, 31, 19–26. [Google Scholar] [CrossRef] [Green Version]
  69. Nannipieri, P.; Ascher, J.; Ceccherini, M.T.; Landi, L.; Pietramellara, G.; Renella, G. Microbial diversity and soil functions. Eur. J. Soil Sci. 2003, 54, 655–670. [Google Scholar] [CrossRef]
  70. Joanisse, G.D.; Bradley, R.L.; Preston, C.M.; Bending, G.D. Sequestration of soil nitrogen as tannin-protein complexes may improve the competitive ability of sheep laurel (Kalmia angustifolia) relative to black spruce (Picea mariana). New Phytol. 2009, 181, 187–198. [Google Scholar] [CrossRef] [PubMed]
  71. Averill, C.; Turner, B.L.; Finzi, A.C. Mycorrhiza-mediated competition between plants and decomposers drives soil carbon storage. Nature 2014, 505, 543. [Google Scholar] [CrossRef]
  72. Legay, N.; Baxendale, C.; Grigulis, K.; Krainer, U.; Kastl, E.; Schloter, M.; Bardgett, R.D.; Arnoldi, C.; Bahn, M.; Dumont, M.; et al. Contribution of above- and below-ground plant traits to the structure and function of grassland soil microbial communities. Ann. Bot. 2014, 114, 1011–1021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Graham, E.B.; Knelman, J.E.; Schindlbacher, A.; Siciliano, S.; Breulmann, M.; Yannarell, A.; Beman, J.M.; Abell, G.; Philippot, L.; Prosser, J.; et al. Microbes as engines of ecosystem function: When does community structure enhance predictions of ecosystem processes? Front. Microbiol. 2016, 7, 214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Constrained correspondence analysis showing the effects of (A) time, litter chemistry, climate and space (14.8% of total inertia) on microbial composition of fungi (*), NifH (+) and AmoA (X) communities. Partially constrained correspondence analysis (B) removing all effects except climate (2.75% of total inertia) or (C) litter chemistry (3.2% of total inertia) is also shown. The effects of both time and litter chemistry (D) on only fungal community composition in poplar litter (7.6% of total inertia). Point size decreases with time in panels A, C and D. Polygons represent k-means clustering based on CCAs 1 and 2.
Figure 1. Constrained correspondence analysis showing the effects of (A) time, litter chemistry, climate and space (14.8% of total inertia) on microbial composition of fungi (*), NifH (+) and AmoA (X) communities. Partially constrained correspondence analysis (B) removing all effects except climate (2.75% of total inertia) or (C) litter chemistry (3.2% of total inertia) is also shown. The effects of both time and litter chemistry (D) on only fungal community composition in poplar litter (7.6% of total inertia). Point size decreases with time in panels A, C and D. Polygons represent k-means clustering based on CCAs 1 and 2.
Jof 08 00735 g001
Figure 2. Fraction of microbial community (A) structure and (B) composition explained by litter chemistry, time, climate, spatial PCNM and shared spatio-climatic effects, using partially constrained ordination. Data used in analyses either included all functional groups (all sp.) and both litter species (both), or subsets of data including only a specific functional group (fungi, NifH or AmoA) and only poplar or fir litter samples.
Figure 2. Fraction of microbial community (A) structure and (B) composition explained by litter chemistry, time, climate, spatial PCNM and shared spatio-climatic effects, using partially constrained ordination. Data used in analyses either included all functional groups (all sp.) and both litter species (both), or subsets of data including only a specific functional group (fungi, NifH or AmoA) and only poplar or fir litter samples.
Jof 08 00735 g002
Figure 3. Partially constrained CCA showing the effects of (A) time, (B) site and (C) CT level on microbial composition of fungi (*), NifH (+) and AmoA (X) communities associated with poplar and Douglas-fir (Fd) leaf litter with low and high nitrogen (LN and HN) and condensed tannin (LT and HT, poplar only) content and after 7, 12, 24 and 43 months of decay. Litter chemistry effects on (D) fungi communities associated with poplar litter only, where point size decreases with time. Polygons represent k-means clustering based on CCAs 1 and 2.
Figure 3. Partially constrained CCA showing the effects of (A) time, (B) site and (C) CT level on microbial composition of fungi (*), NifH (+) and AmoA (X) communities associated with poplar and Douglas-fir (Fd) leaf litter with low and high nitrogen (LN and HN) and condensed tannin (LT and HT, poplar only) content and after 7, 12, 24 and 43 months of decay. Litter chemistry effects on (D) fungi communities associated with poplar litter only, where point size decreases with time. Polygons represent k-means clustering based on CCAs 1 and 2.
Jof 08 00735 g003
Figure 4. Principal components of neighbor matrices (PCNMs) of spatial layout of field sites. PCNMs 1, 3, 4, 5 and 6 were significantly associated with microbial composition, while PCNMs 1, 2, 3, 4 and 5 were significantly associated with community structure. Point size is proportional to PCNM values for each sampling location.
Figure 4. Principal components of neighbor matrices (PCNMs) of spatial layout of field sites. PCNMs 1, 3, 4, 5 and 6 were significantly associated with microbial composition, while PCNMs 1, 2, 3, 4 and 5 were significantly associated with community structure. Point size is proportional to PCNM values for each sampling location.
Jof 08 00735 g004
Table 1. Initial litter chemistry of Pseudotsuga menziesii (fir = Fd) and Populus angustifolia (poplar) leaf litter with low and high nitrogen (LN and HN) and condensed tannins (LT and HT, poplar only), including carbon (C), nitrogen (N) and their ratio (C/N), phenols, total condensed tannins (CT) and proximate chemistry. Mean values for each litter type were used for modeling microbial community responses during decay. Values (±SE) are based on dry weight (70 °C overnight) and calculated from 6 and 4 replicates (for poplar and Fd litter, respectively) or 3 replicates (for proximate chemical fractions).
Table 1. Initial litter chemistry of Pseudotsuga menziesii (fir = Fd) and Populus angustifolia (poplar) leaf litter with low and high nitrogen (LN and HN) and condensed tannins (LT and HT, poplar only), including carbon (C), nitrogen (N) and their ratio (C/N), phenols, total condensed tannins (CT) and proximate chemistry. Mean values for each litter type were used for modeling microbial community responses during decay. Values (±SE) are based on dry weight (70 °C overnight) and calculated from 6 and 4 replicates (for poplar and Fd litter, respectively) or 3 replicates (for proximate chemical fractions).
SpeciesLitter TypeC
(mg g−1)
N
(mg g−1)
C/NPhenols
(mg g−1)
CT
(mg g−1)
Populus angustifoliaLT LN457.0 ± 1.18.1 ± 0.256.39 ±1.4518.3 ± 3.466.4 ± 8.1
LT HN461.4 ± 0.715.7 ± 0.329.42 ±0.5519.2 ± 4.485.8 ± 6.4
HT LN481.8 ± 1.26.9 ± 0.270.36 ±1.7260.8 ± 6.8264.9 ± 7.9
HT HN477.5 ± 1.314.3 ± 0.133.32 ±0.4054.3 ± 4.8291.5 ± 5.3
Pseudotsuga menziesiiFd LN512.6 ± 1.56.1 ± 0.183.52 ±1.3719.0 ± 3.120.8 ± 0.4
Fd HN511.0 ± 0.511.2 ± 0.245.46 ±0.7618.4 ± 1.520.9 ±0.3
SpeciesLitter TypeNPE a
(mg g−1)
WSE b
(mg g−1)
AHF c
(mg g−1)
AUR d
(mg g−1)
Ash
(mg g−1)
Populus angustifoliaLT LN137.6 ± 15.2130.7 ± 0.1304.2 ± 5.0414.8 ± 6.612.7 ± 0.6
LT HN152.4 ± 5.6128.4 ± 7.5307.8 ± 4.4398.1 ± 2.713.3 ± 1.6
HT LN177.7 ± 15.2154.0 ± 11.1283.8 ± 10.9373.1 ± 5.411.4 ± 0.8
HT HN230.2 ± 6.5132.3 ± 12.4274.3 ± 3.0353.0 ± 5.510.1 ± 0.5
Pseudotsuga menziesiiFd LN99.5 ± 3.5119.5 ± 1.0362.6 ± 0.6394.0 ± 3.424.4 ± 0.6
Fd HN121.6 ± 3.2117.2 ± 1.5355.1 ± 1.2382.6 ± 0.523.4 ± 0.4
a Non-polar extractables, b water-soluble extractables, c acid-hydrolyzable forage, d acid-unhydrolyzable residue.
Table 2. In situ climatic conditions measured at nine sites spanning three latitudes (north, central, south) and two biogeoclimatic (BGC) zones (coastal western hemlock (WH), coastal Douglas-fir (DF) or transition between the two (TR)) over 48 months (1 May 2011 to 30 April 2015).
Table 2. In situ climatic conditions measured at nine sites spanning three latitudes (north, central, south) and two biogeoclimatic (BGC) zones (coastal western hemlock (WH), coastal Douglas-fir (DF) or transition between the two (TR)) over 48 months (1 May 2011 to 30 April 2015).
LatitudeBEC ZoneSoil Moisture (%) aSoil Temp. (°C) aAir Temp. (°C) bPET c
NorthDF13.93
[3.85, 40.97]
9.55
[3.73, 15.37]
8.99
[−4.86, 23.90]
88.2
TR12.97
[2.77, 45.56]
9.64
[2.84, 15.85]
9.60
[−4.30, 27.18]
92.5
WH9.05
[3.03, 20.62]
8.97
[2.96, 14.71]
8.97
[−5.90, 26.32]
90.5
CentralDF9.83
[2.85, 23.07]
9.87
[3.37, 16.92]
10.09
[−3.52, 26.62]
93.6
TR9.2
[1.92, 19.52]
9.68
[2.47, 17.09]
9.44
[−6.67, 29.03]
92.7
WH13.86
[3.50, 26.84]
8.58
[2.14, 15.29]
8.72
[−5.56, 31.26]
90.3
SouthDF11.93
[4.03, 24.72]
9.94
[3.91, 16.18]
9.59
[−5.18, 30.00]
89.3
TR13.1
[3.47, 25.85]
8.96
[2.76, 14.48]
8.44
[−6.97, 28.16]
87.5
WH11.25
[6.00, 18.68]
8.40
[2.31, 13.55]
8.61
[−5.97, 28.50]
87.8
a Mean annual [mean of annual minimums, maximums] soil moisture or temperature measured 30 cm below forest floor (abbreviated as MS, MSmin and MSmax for soil moisture and TS, TSmin and TSmax for soil temperature). b Mean annual [mean of annual minimums, maximums] air temperature measured 30 cm above forest floor (abbreviated TA, TAmin and TAmax). c Mean potential evapotranspiration from May to September, using the Thornthwaite method.
Table 3. Significance (p-values) of categorical models and predictors in predicting microbial community structure (A; Shannon’s diversity, richness and Pielou’s evenness) or composition (B; presence/absence). Models constrained by categorical time, litter type and site predictors (constrained by treatment) or by the DNA concentration and λ absorbance (constrained by DNA) were reduced until only significant predictors remained (α-level = 0.05), with non-significant variables labelled ‘ns’. Best fit categorical models were then assessed after removing the effects of DNA quality on microbial responses (partially constrained by DNA), with ‘NA’ labeling cases where DNA quality is not significant. Data from poplar and Douglas-fir litter were treated separately (poplar only and fir only, respectively) and jointly (poplar and fir). Fungi, AmoA and NifH functional groups were also treated separately and jointly (all taxa).
Table 3. Significance (p-values) of categorical models and predictors in predicting microbial community structure (A; Shannon’s diversity, richness and Pielou’s evenness) or composition (B; presence/absence). Models constrained by categorical time, litter type and site predictors (constrained by treatment) or by the DNA concentration and λ absorbance (constrained by DNA) were reduced until only significant predictors remained (α-level = 0.05), with non-significant variables labelled ‘ns’. Best fit categorical models were then assessed after removing the effects of DNA quality on microbial responses (partially constrained by DNA), with ‘NA’ labeling cases where DNA quality is not significant. Data from poplar and Douglas-fir litter were treated separately (poplar only and fir only, respectively) and jointly (poplar and fir). Fungi, AmoA and NifH functional groups were also treated separately and jointly (all taxa).
(A) Structure Poplar and Fir Poplar Only Fir Only a
All taxaFungiNifHAmoAAll taxaFungiNifHAmoAAll taxaFungiNifH
Constrained by treatmentstime<0.001<0.001<0.001ns<0.001<0.001<0.001ns<0.001<0.001ns
litter type<0.0010.0030.005nsnsnsnsnsnsnsns
site<0.001<0.0010.0080.003<0.001<0.0010.031<0.0010.0130.026ns
model<0.001<0.001<0.001<0.001<0.001<0.0010.001<0.001<0.001<0.001ns
Constrained by DNAconcentrationnsnsns0.020.027nsns0.0080.043ns0.018
260/280 b0.0030.0330.005nsnsns0.036nsnsnsns
260/230 bnsnsnsnsnsnsnsns0.0260.002ns
model0.0020.0290.0050.0230.021ns0.0350.0090.0160.0020.018
Partially constrained by DNAtime<0.001<0.001<0.001ns<0.001NA0.001ns<0.001<0.001ns
litter type0.0030.0340.048nsnsNAnsnsnsnsns
site<0.001<0.0010.008<0.001<0.001NA0.022<0.0010.0680.047ns
model<0.001<0.001<0.001<0.001<0.001NA0.001<0.001<0.001<0.001ns
(B) Composition Poplar and Fir Poplar Only Fir Only a
All taxaFungiNifHAmoAAll taxaFungiNifHAmoAAll taxaFungiNifH
Constrained by treatmentstime0.005<0.001<0.001ns<0.001<0.001<0.001ns<0.001<0.001<0.001
litter type0.005<0.001<0.001ns0.0370.035nsns0.040.032ns
site0.005<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.0010.002
model<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001
Constrained by DNAconcentration0.005<0.0010.0340.0090.0020.008ns<0.001ns0.029ns
260/280 b0.0070.009<0.001nsnsns0.004nsnsns<0.007
260/230 bnsnsnsns0.03ns0.007nsnsnsns
model<0.001<0.001<0.0010.0060.0120.006<0.001<0.001ns0.0430.008
Partially constrained by DNAtime<0.001<0.001<0.001ns<0.001<0.001<0.001nsNA<0.001<0.001
litter type<0.001<0.001<0.001ns0.0450.041nsnsNA0.076ns
site<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001NA<0.001<0.001
model<0.001<0.001<0.001<0.001<0.001<0.001<0.001<0.001NA<0.001<0.001
aAmoA genes were not detected in fir litter, therefore the column was omitted. b Wavelength absorption ratios of extracted DNA solutions.
Table 4. Significance (p-values) of continuous models and predictors in predicting microbial community structure (A; Shannon’s diversity, richness and Pielou’s evenness) or composition (B; presence/absence). Partially constrained ordination analyses were used to test the effects of continuous time, initial litter chemistry, climate or principal coordinates of neighbor matrix (PCNM) predictors, while removing other respective effects. Climate and PCNM variables (and litter chemistry in the case of joint poplar and Fd analyses) were first reduced until only significant predictors remained (α-level = 0.05). Cases where no variables were deemed significant and when all of the variance accorded to a variable is confounded by partial constraints are labeled by ‘ns’. Data from poplar and Douglas-fir litter were treated separately (poplar only and fir only, respectively) and jointly (poplar and fir). Fungi, AmoA and NifH functional groups were also treated separately and jointly (all taxa). The column for AmoA is not included when using only Fd data since this gene was never detected in this litter type. Mean annual range in soil temperature is designated ‘TSrange’. Refer to Table 1 and Table 2 for other abbreviation definitions.
Table 4. Significance (p-values) of continuous models and predictors in predicting microbial community structure (A; Shannon’s diversity, richness and Pielou’s evenness) or composition (B; presence/absence). Partially constrained ordination analyses were used to test the effects of continuous time, initial litter chemistry, climate or principal coordinates of neighbor matrix (PCNM) predictors, while removing other respective effects. Climate and PCNM variables (and litter chemistry in the case of joint poplar and Fd analyses) were first reduced until only significant predictors remained (α-level = 0.05). Cases where no variables were deemed significant and when all of the variance accorded to a variable is confounded by partial constraints are labeled by ‘ns’. Data from poplar and Douglas-fir litter were treated separately (poplar only and fir only, respectively) and jointly (poplar and fir). Fungi, AmoA and NifH functional groups were also treated separately and jointly (all taxa). The column for AmoA is not included when using only Fd data since this gene was never detected in this litter type. Mean annual range in soil temperature is designated ‘TSrange’. Refer to Table 1 and Table 2 for other abbreviation definitions.
(A) Structure Poplar and Fir Poplar Only Fir Only a
All taxaFungiNifHAmoAAll taxaFungiNifHAmoAAll taxaFungiNifH
Partially constrained model significancetime<0.001<0.001<0.0010.366<0.001<0.001<0.0010.346<0.001<0.0010.1
litter chem.<0.0010.004<0.0010.0040.5080.3760.6110.1830.4610.4510.734
climate0.1920.5460.1210.6310.1880.0640.1620.614nsnsns
PCNM0.120.0110.0760.0080.012<0.0010.1570.0020.012nsns
Predictor significancelitter chem.N (0.026)
CT (0.008)
AUR (<0.001)
CT (0.027)
AUR (0.029)
Ash (<0.001)
N (0.007)CT (0.564)
N (0.297)
CT (0.173)
N (0.85)
CT (0.766)
N (0.341)
CT (0.711)
N b (0.078)
N (0.484)N (0.424)N (0.717)
climateMS (0.068)
MSmin (0.075)
MSmax (0.72)
TA (0.735)
TAmin (ns)
MSmin (0.239)
TAmin (0.959)
TA (0.121)MS (0.771)
TA (0.378)
MS (0.131)
MSmax (0.433)
TA
(0.3)
MSmin (0.071)TA (0.195)MS (0.796)
TA (0.332)
nsnsns
PCNMPCNM 1
(0.007)
PCNM 3
(0.826)
PCNM 4
(0.682)
PCNM 5
(ns)
PCNM 1
(0.003)
PCNM 3
(0.672)
PCNM 4
(0.462)
PCNM 3
(0.078)
PCNM 1
(0.003)
PCNM 2
(0.583)
PCNM 5
(0.122)
PCNM 1
(<0.001)
PCNM 3
(0.219)
PCNM 5
(0.832)
PCNM 1
(0.001)
PCNM 3
(0.429)
PCNM 3
(0.151)
PCNM 1
(0.002)
PCNM 2
(0.512)
PCNM 5
(0.044)
PCNM 4
(0.01)
nsns
(B) Composition Poplar and Fir Poplar Only Fir Only a
All taxaFungiNifHAmoAAll taxaFungiNifHAmoAAll taxaFungiNifH
Model significancetime<0.001<0.001<0.0010.033<0.001<0.001<0.0010.031<0.001<0.001<0.001
litter chem.<0.001<0.001<0.0010.0020.0230.0510.7180.2280.080.0740.268
climate<0.0010.03<0.0010.008<0.0010.153<0.0010.0020.003<0.0010.003
PCNM0.004<0.001<0.001<0.0010.002<0.0010.002<0.0010.296ns0.131
Predictor significancelitter chem.N (<0.001)
CT (<0.001)
AUR (<0.001)
CT (0.003)
AUR (<0.001)
Ash (<0.001)
Ash (<0.001)N (0.002)CT d (0.162)
N (0.024)
CT b,c (0.059)
N b (0.141)
CT (0.577)
N (0.692)
CT (0.533)
N b (0.132)
N (0.076)N (0.077)N (0.271)
climateMSmin (0.009)
MSmax (0.005)
TA (0.03)
TAmin (0.683)
TSrange d (ns)
PET
(ns)
MSmin (0.012)
TSrange (0.334)
MSmin (0.034)
TAmin (0.026)
TSmin (<0.001)
TSrange (0.016)
PET
(ns)
TAmin (0.013)MSmin (0.016)
MSmax (0.006)
TA (0.014)
TAmin (0.628)
TSrange (ns)
PET
(ns)
MSmin (0.474)
MSmax (0.078)
MSmin (0.005)
TA
(0.01)
TAmin (0.058)
TSrange (0.001)
PET (0.006)
TAmin (0.001)
TAmax (0.049)
TAmin (0.076)
TSrange (0.005)
TAmax (0.003)
TSrange (0.01)
MS (0.095)
TAmin (0.002)
PCNMPCNM 1 (0.002)
PCNM 3 (0.215)
PCNM 5
(ns)
PCNM 6
(ns)
PCNM 1 (<0.001)
PCNM 6
(0.01)
PCNM 1 (0.023)
PCNM 3 (0.002)
PCNM 4 (0.093)
PCNM 5
(ns)
PCNM 1 (<0.001)
PCNM 5 (0.001)
PCNM 1 (0.004)
PCNM 3 (0.074)
PCNM 5
(ns)
PCNM 6
(ns)
PCNM 1 (<0.001)
PCNM 3 (0.251)
PCNM 6 (0.007)
PCNM 1
(0.139)
PCNM 3
(0.001)
PCNM 4
(0.078)
PCNM 1 (0.002)
PCNM 5 (<0.001)
PCNM 3 (0.304)nsPCNM 3 (0.131)
aAmoA genes were not detected in fir litter, therefore the column was omitted. b Significant when only considering litter decaying for 7, 12 and 24 months. c Significant when using categorical time predictors. d Significant when using categorical time predictors while only considering litter decaying for 7, 12 and 24 months.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shay, P.-E.; Winder, R.S.; Constabel, C.P.; Trofymow, J.A. Fungal Community Composition as Affected by Litter Chemistry and Weather during Four Years of Litter Decomposition in Rainshadow Coastal Douglas-Fir Forests. J. Fungi 2022, 8, 735. https://0-doi-org.brum.beds.ac.uk/10.3390/jof8070735

AMA Style

Shay P-E, Winder RS, Constabel CP, Trofymow JA. Fungal Community Composition as Affected by Litter Chemistry and Weather during Four Years of Litter Decomposition in Rainshadow Coastal Douglas-Fir Forests. Journal of Fungi. 2022; 8(7):735. https://0-doi-org.brum.beds.ac.uk/10.3390/jof8070735

Chicago/Turabian Style

Shay, Philip-Edouard, Richard S. Winder, C. Peter Constabel, and J. A. (Tony) Trofymow. 2022. "Fungal Community Composition as Affected by Litter Chemistry and Weather during Four Years of Litter Decomposition in Rainshadow Coastal Douglas-Fir Forests" Journal of Fungi 8, no. 7: 735. https://0-doi-org.brum.beds.ac.uk/10.3390/jof8070735

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop