Who is who in litter decomposition? Metaproteomics reveals major microbial players and their biogeochemical functions.

Leaf-litter decomposition is a central process in carbon cycling; however, our knowledge about the microbial regulation of this process is still scarce. Metaproteomics allows us to link the abundance and activity of enzymes during nutrient cycling to their phylogenetic origin based on proteins, the 'active building blocks' in the system. Moreover, we employed metaproteomics to investigate the influence of environmental factors and nutrients on the decomposer structure and function during beech litter decomposition. Litter was collected at forest sites in Austria with different litter nutrient content. Proteins were analyzed by 1-D-SDS-PAGE followed by liquid-chromatography and tandem mass-spectrometry. Mass spectra were assigned to phylogenetic and functional groups by a newly developed bioinformatics workflow, assignments being validated by complementary approaches. We provide evidence that the litter nutrient content and the stoichiometry of C:N:P affect the decomposer community structure and activity. Fungi were found to be the main producers of extracellular hydrolytic enzymes, with no bacterial hydrolases being detected by our metaproteomics approach. Detailed investigation of microbial succession suggests that it is influenced by litter nutrient content. Microbial activity was stimulated at higher litter nutrient contents via a higher abundance and activity of extracellular enzymes.

Introduction

Leaf-litter decomposition is the physical and chemical breakdown of dead plant material, a key ecosystem process that provides energy for microbial growth, releases nutrients for plant growth, influences ecosystem class="Chemical">carbon (C) storage, aclass="Chemical">nd thus may have a loclass="Chemical">ng-term efclass="Chemical">n class="Chemical">fect on the climate (Chapin ). Investigations of litter decomposition at different geographic locations have shed light on the factors that influence decomposition rates, for example, that litter decomposition occurs above a certain threshold in mean annual temperature and depends on litter quality (Zhang ; Prescott, 2010). In spite of the vast amount of litter decomposition studies, little is known about specific functions contributed by the different microbial groups involved in litter decomposition, mainly due to the lack of suitable methods available for in-depth studies of ecosystem functioning.

Leaf litter consists mainly of class="Chemical">cellulose, but also class="Chemical">n class="Chemical">lignin, hemicellulose, pectin and proteins (Yadav and Malanson, 2007), requiring numerous enzymes for degradation. The presence and activity of extracellular enzymes link environmental nutrient availability to microbial nutrient demand (Sinsabaugh ), which in turn is determined by the elemental stoichiometry (that is, the C:N:P ratio) of the decomposer's biomass (Sterner and Elser, 2002). Within the microbial biomass, phylogenetic groups such as fungi and bacteria exhibit different nutrient demands and thereby influence decomposition and the carbon use efficiency within the ecosystem (Sinsabaugh ; Keiblinger ). Studies on extracellular enzyme activities may reflect the constraints on microbial biomass stoichiometry and enzyme relationships to litter decomposition, that is, for hydrolytic enzymes ratios of C, N and P acquisition activities have been assessed and they converged on 1:1:1 (Sinsabaugh ). However, so far no attempt has been made to link enzyme activities to the respective protein abundances, their microbial origin, or the biogeochemistry of leaf litter.

In the last decade, metaproteomics has greatly advanced our understanding of microorganisms and their geochemical envclass="Chemical">ironmeclass="Chemical">nt (reviewed iclass="Chemical">n Newmaclass="Chemical">n aclass="Chemical">nd Baclass="Chemical">nfield, 2002; Schclass="Chemical">neider aclass="Chemical">nd Riedel, 2010; Wilmes aclass="Chemical">nd Boclass="Chemical">nd, 2009). Iclass="Chemical">n coclass="Chemical">ntrast to DNA aclass="Chemical">nd RNA, most proteiclass="Chemical">ns have aclass="Chemical">n iclass="Chemical">ntriclass="Chemical">nsic metabolic fuclass="Chemical">nctioclass="Chemical">n aclass="Chemical">nd caclass="Chemical">n thus be used to relate specific microbial activities to deficlass="Chemical">ned orgaclass="Chemical">nisms iclass="Chemical">n multispecies commuclass="Chemical">nities. Therefore, the ideclass="Chemical">ntificatioclass="Chemical">n of the microbial proteiclass="Chemical">ns of a giveclass="Chemical">n habitat together with the aclass="Chemical">nalysis of their phylogeclass="Chemical">netic origiclass="Chemical">n aclass="Chemical">nd their temporal distributioclass="Chemical">n is expected to provide fuclass="Chemical">ndameclass="Chemical">ntally class="Chemical">new iclass="Chemical">nsights iclass="Chemical">nto the role of microbial diversity iclass="Chemical">n biogeochemical processes. Receclass="Chemical">nt large-scale characterisatioclass="Chemical">ns of the eclass="Chemical">ntire proteiclass="Chemical">n complemeclass="Chemical">nt of diverse microbial commuclass="Chemical">nities have proveclass="Chemical">n useful to describe microbial fuclass="Chemical">nctioclass="Chemical">ns iclass="Chemical">n soil (Schulze ; Waclass="Chemical">ng ), lake aclass="Chemical">nd grouclass="Chemical">nd class="Chemical">n class="Chemical">water (Schulze ; Benndorf , 2009), and leaf phyllosphere (Delmotte ). We have recently employed a proteomics approach to investigate a model system, consisting of the litter degrading fungus Aspergillus nidulans and the bacterium Pectobacterium carotovorum grown on sterilized leaf litter (Schneider ). In the current study, we expand our proteomics analyses to environmental litter samples.

The major goals of the present study were: (I) to link structure and function of microbial communities involved in leaf-litter decomposition and (II) to investigate the impact of leaf-litter nutrient content (stoichiometry) and season on the abundance and activity of the leaf-litter decomposers. We analyzed the metaproteome from beech (class="Species">Fagus sylvatica L.) litter collected iclass="Chemical">n class="Chemical">n class="Chemical">February and May 2009 at four Austrian sampling-sites that differed in their nutrient content. The newly developed bioinformatics pipeline ‘PROteomics result Pruning & Homology group ANotation Engine' (PROPHANE, Schneider ) was employed to assign proteins to their phylogenetic and functional origin. Phospholipid fatty acid (PLFA) analyses, enzymatic measurements and qualitative assessment of fungal communities were used for data validation.

Materials and methods

Litter sampling sites

Beech leaf litter was collected at four forest sites in Austria: (I) Achenkirch (AK), Tyrol; (II) Klausenleopoldsdorf (KL), Lower Austria; (III) Ort (OR), Gmunden, Upper Austria; and (IV) Schottenwald (SW), Vienna that difclass="Chemical">fered iclass="Chemical">n their class="Chemical">nutrieclass="Chemical">nt coclass="Chemical">nteclass="Chemical">nt. AK is located iclass="Chemical">n the Acheclass="Chemical">ntal of the North Tyroleaclass="Chemical">n limestoclass="Chemical">ne Alps. It is characterised by its locatioclass="Chemical">n oclass="Chemical">n calcareous bedrock that results iclass="Chemical">n a high soil pH (Kitzler ). OR is situated at aclass="Chemical">n elevatioclass="Chemical">n of about 700 m above sea level aclass="Chemical">nd class="Chemical">northeast orieclass="Chemical">ntatioclass="Chemical">n. The forest coclass="Chemical">nsists mostly of beech, but also of various class="Chemical">n class="Species">conifers, and has an average slope angle of 34%. The soil is a brown loam on carbonatic rock. SW is situated in direct vicinity of Vienna (Kitzler ) similar to KL that is located about 40 km south-west of Vienna (Kitzler ). The predominant vegetation, stand age, exposition and elevation, as well as soil type and texture of the different sampling sites are described in Table 1.

Table 1

Site and soil characteristics of the beech litter-sampling sites

	Achenkirch (Wanek et al., 2010)	Klausenleopoldsdorf (Wanek et al., 2010)	Ort	Schottenwald (Wanek et al., 2010)
Location	47°35′N 11°39′E	48°07′N 16°03′E	47°51′N 13°42′E	48°14′N 16°15′E
Vegetation	Spruce-fir, beech	Beech	Mainly beech, conifers	Beech
Stand age (years)	127	62	90	142
Exposition	N	NNE	NE	SE
Elevation (m a.s.l.)	895	510	700	370
Soil type	Rendzic leptosols/chromic cambisol	Dystric cambisol over sandstone	Cambisol	Dystric cambisol over sandstone
Soil texture	Loam	Loam-loamy clay	Loamy silt	Silty loam
Soil pH 0–7 cm	6.4	4.6	3.7	4.4

Abbreviation: a.s.l., above sea level.

Litter sampling

At the respective sampling sites, 500 g of beech leaf litter was sampled in class="Chemical">February aclass="Chemical">nd May 2009 (exact sampliclass="Chemical">ng dates are listed iclass="Chemical">n Table 2). Iclass="Chemical">n class="Chemical">n class="Chemical">February, all sites were covered by snow, whereas in May all sites were free of snow. Three biological replicates were collected from three plots of 10 × 10 m within an area of 50 × 50 m. Only the L-horizon was sampled carefully, but excluding mineral soil material and O-horizon. Litter samples consisted predominantly of last autumn's foliage, which means that leaves were almost intact as beech litter at the investigated sites is completely decomposed during the course of 1 year. The litter was transported to the laboratory at ambient temperature, sorted, that is, leaves other than beech leaves were removed, cut into <0.5 cm2 pieces using a kitchen knife and mixed thoroughly. Per sampling site, 10–50 g aliquots of cut leaf material from the three biological replicates per sampling site were stored at −20 °C until litter nutrient content had been analyzed, and at −80 °C until the metaproteome was analyzed.

Table 2

Different sampling sites and their mean air-temperature (T) and precipitation (Prec)

Site	SamplingFeb	SamplingMay	TFeb (°C)	TMay (°C)	PrecFeb (mm)	PrecMay (mm)
Achenkirch	4th	18th	−1.66	5.88	55	65
Klausenleopoldsdorf	5th	27th	−0.11	9.64	37	46
Ort	27th	19th	0.88	11.0	51	57
Schottenwald	5th	27th	1.32	11.39	44	55

Abbreviation: Feb, February.

Samples were collected at the given dates in 2009. Climatic data were collected over a 3-month time period before litter sampling. Mean air temperature and precipitation were calculated by averaging the climatic data of this period.

Determination of environmental parameters and litter nutrient quality

Air temperature and precipitation at the respective sampling sites were monitored for 3 months before litter collection; mean air-temperature (T) and mean precipitation (Prec) were calculated by averaging these values.

class="Chemical">Water coclass="Chemical">nteclass="Chemical">nt (WC) was determiclass="Chemical">ned gravimetrically by dryiclass="Chemical">ng the litter at 105°C for 24 h. Dry samples were grouclass="Chemical">nd with a mill (Retsch MM2000, Haclass="Chemical">nau, Germaclass="Chemical">ny) to a ficlass="Chemical">ne homogeclass="Chemical">neous powder. The total class="Chemical">n class="Chemical">carbon (C) and nitrogen (N) contents of the litter were analyzed with an integrated oxidation and detection device (Leco CN2000, LECO corp. St Joseph, MI, USA). The ground samples were wet acidically oxidized (H2SO4+HNO3) in a microwave oven (CEM MARS Express), and elements (phosphorous (P), potassium (K), magnesium (Mg), manganese (Mn), calcium (Ca), iron (Fe)) were determined by inductively coupled plasma atomic emission spectrometry as described by Henschler (1988). Extractable nitrogen (N) in form of ammonium (NH4-N) and nitrate (NO3-N) and phosphorous (P) in form of phosphate (PO4-P) were determined from freshly shred litter samples, extracted with 1 M KCl at a ratio of 1:20, by shaking on a reciprocal shaker for 60 min at 200 oscillations. Extracts were stored at −20 °C until analysis. Ammonium was measured by the Berthelot reaction according to Schinner . Nitrate was determined with the vanadium chloride method as described by Hood-Nowotny . Phosphate was quantified photometrically based on the phosphomolybdate blue reaction (Schinner ) in a microplate reader (μQuant mQx200, Bio-Tek Instruments, Winooski, VT, USA).

PLFA analysis

class="Chemical">Lipids were extracted from 1 g sub-samples usiclass="Chemical">ng a modified Bligh aclass="Chemical">nd Dyer techclass="Chemical">nique (Hackl ). Extracted class="Chemical">n class="Chemical">PLFAs were analyzed with a HP 6890 Series GC-System and a 7683 series injector and auto sampler on a HP-5 capillary column and detected with a flame ionization detector. For identification of the fatty acid methyl esters, an external standard (bacterial acid methyl ester mix from SUPELCO) was used. For quantification of the peaks, methyl non-adecanoate fatty acid (19:0) was added. PLFA nomenclature is based upon Frostegard . The ratio fungal/bacterial PLFA was calculated by dividing the amount of 18:2ω6 through the amount of total bacterial PLFAs (Frostegard and Baath, 1996).

Cellulase and xylanase activity assay

Frozen cut leaf litter was ground in liquid class="Chemical">nitrogen; 600 mg of the grouclass="Chemical">nd material were extracted with 4 ml extractioclass="Chemical">n bufclass="Chemical">n class="Chemical">fer (0.1 M sodium phosphate pH 6.0, 0.1% Triton X-100, 0.1% polyvinylpolyrrolidone). After shaking at 4 °C for 1 h, extracts were centrifuged for 5 min at 14 000 g. The resulting supernatants were filter-sterilized by cellulose acetate filters and directly subjected to enzyme activity assays.

Cellulase and class="Chemical">xylanase activity were determiclass="Chemical">ned as described earlier (Riedel ; Köclass="Chemical">nig ). 250 μl filter-sterilized superclass="Chemical">nataclass="Chemical">nt were mixed with 250 μl 1% (w/v) class="Chemical">n class="Chemical">carboxymethyl cellulose (Sigma, Buchs, Switzerland) or 1.5% (w/v) xylan (from oat spelts, Boehringer Ingelheim, Lot Nr. 06995) in 0.1 M sodium phosphate buffer, pH 6.0, and incubated for 3.5 h (cellulase) or 2 h (xylanase) at 40 °C. The enzymatic release of reducing groups was determined with the dinitrosalicylic acid reagent (Wood and Bhat, 1988). One unit of enzyme activity was defined as the amount of enzyme needed to release 1 μmol of glucose-equivalent reducing groups per minute.

Metaproteome analysis

From the cut and stored (−80 °C) litter material, 5 g were ground in liquid class="Chemical">nitrogen aclass="Chemical">nd mixed with extractioclass="Chemical">n bufclass="Chemical">n class="Chemical">fer containing 1% SDS, 50 mM Tris/KOH, pH 7.0 in a 1:5 ratio (w/v). Samples were sonicated for 2 min, followed by boiling for 20 min and shaking at 4 °C for 1 h. To remove debris, extracts were centrifuged at 3000 g at 4 °C followed by 5 min centrifugation at 14 000 g and 4 °C. Supernatants were concentrated about 5-fold by vacuum-centrifugation (Eppendorf Vacuum Concentrator plus) at 30 °C. Concentrated supernatants of 25 μl were subjected to 1D SDS-PAGE (Laemmli, 1970) in a 12% polyacrylamide gel to clean samples from interfering substances (for example, humic acids) and to reduce sample complexity. Protein lanes were cut into seven slices, and gel slices were subjected to immediate in-gel tryptic digestion by employing sequencing grade-modified trypsin (Promega, Madison, WI, USA, reference V5111) (Shevchenko ). The resulting peptide mixtures were analyzed on a hybrid LTQ-Orbitrap mass spectrometer (ThermoFischer Scientific, Waltham, MA, USA) as described earlier (Schneider ).

Database searches

Using MASCOT (version no.2.2.04), MS and MS/MS data were searched against a database containing all proteins from UniRef100 (9808438 protein entries, downloaded from http://www.ebi.ac.uk/uniref/ at the 26 January 2010) and protein sequence information from a farm silage class="Species">soil metagenome (Triclass="Chemical">nge ) (184 374 eclass="Chemical">ntries, dowclass="Chemical">nloaded from http://img.jgi.doe.gov at 15 October 2009), as well as commoclass="Chemical">n coclass="Chemical">ntamiclass="Chemical">naclass="Chemical">nts such as keraticlass="Chemical">n aclass="Chemical">nd trypsiclass="Chemical">n (iclass="Chemical">n total, 9 993 117 proteiclass="Chemical">n eclass="Chemical">ntries). The followiclass="Chemical">ng search parameters were applied: (i) trypsiclass="Chemical">n was choseclass="Chemical">n as proteiclass="Chemical">n-digesticlass="Chemical">ng eclass="Chemical">nzyme aclass="Chemical">nd up to two missed cleavages were tolerated, (ii) carbamidomethylatioclass="Chemical">n of class="Chemical">n class="Chemical">cysteine was chosen as fixed modification, and (iii) oxidation of methionine was chosen as variable modification. Searches were performed with a parent-ion mass tolerance of ±5 ppm and a fragment-ion mass tolerance of ±0.8 Da. A second database search was performed with the X!Tandem (version 2007.01.01.1, http://www.thegpm.org/TANDEM/) search engine (Craig and Beavis, 2004) with similar settings.

Data processing

Scaffold (version Scaffold 3.0, Proteome Software, Portland, OR, USA) was used to validate and quantify MS/MS-based peptide and protein identifications from both search engines. Protein identifications were accepted if they were established at >95% peptide probability and 90% protein probability with at least one peptide uniquely assigned to a respective protein in one of our samples. Proteins that were identified with the same set of peptides were grouped to protein clusters to satisfy the principles of parsimony. False-discovery rate was determined according to Elias and Gygi (2007) by searching against a composite version of our protein reclass="Chemical">fereclass="Chemical">nce database, created by coclass="Chemical">ncateclass="Chemical">naticlass="Chemical">ng the target proteiclass="Chemical">n sequeclass="Chemical">nces with reversed sequeclass="Chemical">nces (19 986 234 eclass="Chemical">ntries). Starticlass="Chemical">ng from the Scaffold output files, all obtaiclass="Chemical">ned proteiclass="Chemical">n hits were assigclass="Chemical">ned to phylogeclass="Chemical">netic aclass="Chemical">nd fuclass="Chemical">nctioclass="Chemical">nal groups aclass="Chemical">nd assigclass="Chemical">nmeclass="Chemical">nts were validated by the class="Chemical">n class="Chemical">PROPHANE workflow (http://prophane.svn.sourceforge.net/viewvc/prophane/trunk/). Higher protein abundance is represented by a higher number of MS/MS spectra acquired from peptides of the respective protein. Thus, protein abundances were calculated based on the normalized spectral abundance factor (NSAF; Florens ; Zybailov ). This number allows relative comparison of protein abundances over different samples (Bantscheff ).

Statistical analyses

The whole dataset was tested for normal distribution using the Shapiro–Wilk test and homogeneity of variances using Levene's test. When the assumption of normal distribution was violated, the data were transformed as indicated. Difclass="Chemical">fereclass="Chemical">nces betweeclass="Chemical">n sampliclass="Chemical">ng time aclass="Chemical">nd site were either checked by t-test or ANOVA, Tukey class="Chemical">n class="Disease">HSD. We performed single linear regression (SLR) analysis to evaluate which factors predict changes in microbial community structure and function. To identify key-parameters affecting structure and function of the decomposing community, we performed a PCA. All statistical analyses were conducted using PASW (version no. 18.0.0) statistical software packages.

Results

Litter from the sampling sites varied in nutrient content and environmental parameters

The four sampling sites AK, KL, OR and SW showed difclass="Chemical">fereclass="Chemical">nces iclass="Chemical">n their eclass="Chemical">nvclass="Chemical">n class="Chemical">ironmental conditions and leaf-litter quality. At all sampling sites, mean air temperature and precipitation increased from Feb to May. SW was the warmest of the four sites, whereas AK was the coldest (Table 2). Precipitation was highest in AK and lowest in KL at both sampling time points (Table 2).

As litter quality parameters, pH, WC and litter nutrient content were determined. WC was similar in litter of all sites in class="Chemical">February (Table 3), which might be explaiclass="Chemical">ned by the sclass="Chemical">now coverage of the litter. Iclass="Chemical">n May, WCs were geclass="Chemical">nerally lower, but showed site-specific difclass="Chemical">n class="Chemical">ferences. Litter pH was lowest in OR and highest in SW and KL. With respect to nutrient content, litter showed significant site-specific differences in C, N and P (Table 3). Carbon content was highest in AK and lowest in SW, N content was highest in SW and lowest in OR, and P content was highest in SW and lowest in AK. Huge differences were also found for the micronutrients Fe and Mn with highest values in SW and lowest in AK. Litter stoichiometry (C:nutrient and N:P ratios) differed strongly between sites. Generally, SW was the nutrient-richest litter (lowest C:N, C:P and N:P ratios), whereas AK was the nutrient-poorest litter (highest C:P, highest N:P).

Table 3

Beech leaf-litter water content, pH and nutrient content at the four sampling sites in February and May

	Achenkirch		Klausenleopoldsdorf		Ort		Schottenwald
	Feb	May	Feb	May	Feb	May	Feb	May
Water content (%)	73.07±0.6ab	42.2±3.3c	79.39±2.0a	63.7±15.5ab	57.0±4.1bc	12.0±0.3d	77.5±1.6a	18.0±0.6d
pH	n.d.	4.9±0.1cd	4.7±0.1c	5.3±0.1e	4.5±0.1a	4.6±0.1b	5.2±0.1de	5.2±0.1de
C (%)	50.16±0.17d	50.49±0.08d	46.82±1.45bc	46.12±0.68bc	47.20±0.20bc	47.70±0.26c	45.47±0.39b	43.44±1.13a
N (%)	1.24±0.03d	1.27±0.02d	1.00±0.02b	1.14±0.02c	0.82±0.02a	1.07±0.01c	1.25±0.04d	1.45±0.08e
C:N ratio	40.31±0.76c	39.87±0.63c	46.76±2.32e	40.64±0.97c	57.75±1.25f	44.42±0.49d	36.43±1.19b	29.96±1.30a
P (%)	0.035±0.001a	0.039±0.006a	0.040±0.002a	0.040±0.013a	0.035±0.006a	0.056±0.003b	0.064±0.004b	0.104±0.002c
C:P ratio	1439.8±4.8d	1318.6±197.9d	1160.2±86.7cd	1237.4±361.5d	1361.2±200.7d	853.24±42.42bc	714.0±43.0ab	418.6±17.3a
N:P ratio	35.73±0.55e	33.03±4.48de	24.80±0.84bcd	30.33±8.32cde	23.62±3.86bc	19.22±1.16ab	19.59±0.55ab	13.99±0.82a
PO4 (ng g–1 DW)	11.51±6.15a	14.09±1.55a	14.87±4.97a	9.87±4.14a	22.73±5.39a	90.30±6.93b	9.38±2.48a	113.66±6.29c
NH4 (μg g–1 DW)	38.84±10.21b	21.81±3.06a	23.45±3.18a	40.05±8.08b	12.46±0.95a	46.02±4.96b	18.89±5.54a	70.21±5.27c
NO3 (μg g–1 DW)	4.39±1.84abc	1.00±0.48ab	7.31±1.87bc	0.40±0.63a	9.55±4.74c	1.75±0.16ab	5.09±2.14abc	17.37±3.57d
K (%)	0.102±0.005a	0.087±0.008a	0.185±0.018d	0.151±0.007bc	0.170±0.011cd	0.140±0.003b	0.163±0.006c	0.260±0.009e
Ca (%)	1.39±0.08ab	1.32±0.001a	1.42±0.03ab	1.51±0.40ab	1.35±0.08ab	1.47±0.03ab	1.76±0.06b	1.37±0.01ab
Mg (%)	0.209±0.013b	0.173±0.002ab	0.140±0.008ab	0.160±0.069ab	0.141±0.003ab	0.116±0.002a	0.182±0.012ab	0.163±0.002ab
Fe (ppm)	166.6±26.8a	343.8±22.9a	659.1±486.7a	818.5±195.3a	385.9±162.1a	573.6±17.3a	728.7±238.7a	1703.9±391.9b
Mn (ppm)	88.8±2.5a	109.6±2.0a	832.1±89.6b	1025.2±417.2b	1040.5±82.8b	1172.3±16.5b	1493.3±73.9c	2140.7±63.7d
Zn (ppm)	46.08±3.42b	160.00±6.06d	39.44±4.23b	40.83±8.25b	29.17±1.44a	41.17±2.38b	64.33±2.52c	67.19±3.76c

Abbreviations: Feb, February; n.d., not determined because of lack of material.

The different letters in a given line denote significant differences (P<0.05) between values as determined by ANOVA followed by Student–Newman–Keuls multiple range test.

Semi-quantitative metaproteome analysis

When litter samples prepared as described in the materials and methods section were analyzed by class="Chemical">SDS-PAGE, class="Chemical">no disticlass="Chemical">nct baclass="Chemical">nds were detected oclass="Chemical">n the 1D-gels, but proteiclass="Chemical">n separatioclass="Chemical">n resulted iclass="Chemical">n a smear covericlass="Chemical">ng the whole laclass="Chemical">ne (Supplemeclass="Chemical">ntary Figure 1). This might be a result of the high sample complexity aclass="Chemical">nd proteiclass="Chemical">n degradatioclass="Chemical">n iclass="Chemical">n the litter material. Furthermore, class="Chemical">n class="Chemical">humic substances, which interact with the Coomassie stain might be responsible for the intensive blue staining. Lanes were cut into seven equally sized gel slices and analyzed further by LC-MS/MS. Spectra were assigned to a total of 8895 proteins; protein hits that were identified based on the same set of peptides were subsequently merged to 1724 unique protein clusters. Decoy database searches resulted in a protein false-discovery rate below 1.5%. Protein assignments were validated using PROPHANE by testing cluster homology using multiple sequence alignment analyses. This resulted in 204 inconsistent clusters (amino-acid identity <50%), which were therefore excluded from further analyses.

Subsequently, class="Chemical">PROPHANE assigclass="Chemical">ned the remaiclass="Chemical">niclass="Chemical">ng proteiclass="Chemical">n clusters to phylogeclass="Chemical">netic aclass="Chemical">nd fuclass="Chemical">nctioclass="Chemical">nal groups aclass="Chemical">nd quaclass="Chemical">ntified cluster abuclass="Chemical">ndaclass="Chemical">nce by spectral couclass="Chemical">nticlass="Chemical">ng based oclass="Chemical">n the NSAF that is aclass="Chemical">n iclass="Chemical">ndicator of relative difclass="Chemical">n class="Chemical">ferences in protein cluster abundances. All functional or phylogenetic group abundances presented in the figures are based on NSAFs. Proteins belonging to the respective clusters and protein identification parameters are listed in Supplementary Table 1. Supplementary Table 2 shows NSAF-based cluster abundance and a list of representative proteins assigned to the different phyla and functional categories. Supplementary Table 3 shows the sequences and charge states of all peptides that were assigned to proteins or protein clusters. The complete MASCOT results dataset including MS/MS spectra information is provided on the PRIDE database (Vizcaino ) at http://www.ebi.ac.uk/pride/; accession number is 17171.

Litter microbial community differs between sampling sites and seasons

Independent of the sampling site, the majority of the spectra were assigned to three major phylogenetic groups: Viridiplantae (non-degraded leaf proteins, 45–64% of all assigned spectra), fungi (24–51%), and bacteria (5–24%). Moreover, spectra were assigned to Metazoa (3–9%), class="Species">Alveolata (0–1%) aclass="Chemical">nd Archaea (0–0.8%) depeclass="Chemical">ndiclass="Chemical">ng oclass="Chemical">n sampliclass="Chemical">ng-site aclass="Chemical">nd -time (Figure 1a). We used the abuclass="Chemical">ndaclass="Chemical">nce of placlass="Chemical">nt-derived spectra as a marker of the exteclass="Chemical">nd of litter decompositioclass="Chemical">n. The obtaiclass="Chemical">ned phylogeclass="Chemical">netic compositioclass="Chemical">n of the leaf-litter commuclass="Chemical">nity is difclass="Chemical">n class="Chemical">ferent from the proportions of database (DB) protein entries of the respective groups; in the DB, bacteria dominate the number of entries (62.2%) followed by Metazoa (16.7%), fungi (7.8%) and plants (6.3%).

Figure 1

Assignment of spectra to taxonomic groups of organisms at different sampling times and sites. Relative abundances were calculated from the sum of NSAFs found for each group at the respective sampling sites and time. (a) General taxonomy, (b) Bacterial phyla, (c) Proteobacteria, (d) Fungal phyla, (e) Fungal classes. Groups are only presented if their relative abundance is >0.5% in the respective sample. Sampling sites are presented according to increasing mean air-temperature (from left to right).

When comparing samples collected in n class="Chemical">February aclass="Chemical">nd May, a sigclass="Chemical">nificaclass="Chemical">nt decrease (t-test, t=2.64; P=0.019) of the spectra from placlass="Chemical">nts was observed at all sampliclass="Chemical">ng sites. Vice versa, the class="Chemical">number of fuclass="Chemical">ngal spectra iclass="Chemical">ncreased sigclass="Chemical">nificaclass="Chemical">ntly (t-test, t=−2.29; P=0.050). The class="Chemical">number bacterial spectra also iclass="Chemical">ncreased, but this was class="Chemical">not sigclass="Chemical">nificaclass="Chemical">nt.

A comparison of samples from the four sampling sites revealed site-specific difclass="Chemical">fereclass="Chemical">nces of the commuclass="Chemical">nity compositioclass="Chemical">n betweeclass="Chemical">n class="Chemical">n class="Chemical">February and May. The smallest decrease in the proportion of plant spectra was observed in AK (64 to 58%) followed by KL (60 to 45%), SW (45 to 25%) and OR (57 to 27% Figure 1a). Although the proportion of fungal spectra did not increase in AK, the proportion of fungi rose in KL from 25 to 31%, in SW from 32 to 46%, and in OR from 31 to 51%. A similar trend was observed for bacterial spectra (almost no changes in AK, 10 to 19% in KL, 7 to 11% in OR and 17 to 24% in SW). The ratio between the number of fungal and bacterial spectra (F/B-ratio) did not differ between February and May, but site-specific differences were observed. In February, the highest F/B-ratio (5.0) was found in AK and OR (4.5) followed by KL (2.6) and SW (1.9). In May the highest ratio was observed in OR (4.6) followed by AK (3.2), SW (2.4) and KL (1.6).

Fungi and bacteria are thought to be the main degraders of leaf litter and therefore their phylogeny was analyzed in more detail. Most of the bacterial spectra could be assigned to the Proteobacteria (between 59% and 90%) and the Actinobacteria (between 4% and 31% Figure 1b). In class="Chemical">February, γ-Proteobacteria domiclass="Chemical">nated iclass="Chemical">n AK (50%), iclass="Chemical">n KL (57%) aclass="Chemical">nd iclass="Chemical">n SW (37%), whereas iclass="Chemical">n OR, the majority of proteobacterial spectra assigclass="Chemical">nmeclass="Chemical">nts beloclass="Chemical">nged to the α-Proteobacteria (49%, Figure 1c). A sigclass="Chemical">nificaclass="Chemical">nt (t-test, t=−3.18; P=0.0066) absolute iclass="Chemical">ncrease iclass="Chemical">n α-proteobacterial spectra was observed from class="Chemical">n class="Chemical">February to May in AK and KL, whereas in SW and OR, γ- and β-proteobacterial spectra dominated, respectively.

Detailed analysis of the fungal phylogeny revealed only minor difclass="Chemical">fereclass="Chemical">nces iclass="Chemical">n their commuclass="Chemical">nity compositioclass="Chemical">n. A clear domiclass="Chemical">naclass="Chemical">nce (>80%) of ascomycotal spectra was observed at all sites aclass="Chemical">nd times. The proportioclass="Chemical">n of basidiomycotal spectra iclass="Chemical">ncreased slightly from class="Chemical">n class="Chemical">February to May at all sites (Figure 1d) (t-test, t=−2.99; P<0.009). Moreover, a significant decrease in the number of Mucoromycotina-spectra was observed from February to May (t-test, t=3.37; P=0.007) at all sites. No significant changes were observed for fungal classes when different sites and times were compared (Figure 1e); most spectra fit data for the classes Leotiomycetes, Sordariomycetes and Eurotiomycetes.

To validate data obtained from the metaproteomics analysis, class="Chemical">PLFA aclass="Chemical">nalysis was performed aclass="Chemical">nd the phylogeclass="Chemical">netic origiclass="Chemical">n of fuclass="Chemical">ngi preseclass="Chemical">nt iclass="Chemical">n the litter samples was iclass="Chemical">nvestigated by cultivatioclass="Chemical">n aclass="Chemical">nd subsequeclass="Chemical">nt ideclass="Chemical">ntificatioclass="Chemical">n (Supplemeclass="Chemical">ntary Iclass="Chemical">nformatioclass="Chemical">n). To eclass="Chemical">nable comparisoclass="Chemical">n of both datasets, spectra had to be assigclass="Chemical">ned to the major taxoclass="Chemical">nomic groups that caclass="Chemical">n be disticlass="Chemical">nguished by class="Chemical">n class="Chemical">PLFA (Figure 2). A significant correlation of proteome and PLFA data for fungi (SLR, r2=0.380; P=0.014), total bacteria (SLR, r2=0.394; P=0.009) and Gram-negative bacteria (SLR, r2=0.352; P=0.015) was observed; Gram-positive bacteria and Actinobacteria showed similar but insignificant trends (Figure 2).

Figure 2

Comparison of community structure based on NSAF and PLFA analyses. The upper part of the figure shows community structure based on NSAF. Data represent the mean of two biological replicates. The lower part of the figure depicts community structure based on PLFA analysis. Data represent the mean±SD of three biological replicates.

The assignment of spectra to difn class="Chemical">fereclass="Chemical">nt fuclass="Chemical">ngal groups was partly coclass="Chemical">nfirmed by cultivatioclass="Chemical">n experimeclass="Chemical">nts. The ficlass="Chemical">ndiclass="Chemical">ng that the Ascomycota phylum domiclass="Chemical">nates the cultivable fuclass="Chemical">ngi followed by Basidiomycota aclass="Chemical">nd Mucoromycoticlass="Chemical">na is iclass="Chemical">n good accordaclass="Chemical">nce with the metaproteome data (for details see Supplemeclass="Chemical">ntary Iclass="Chemical">nformatioclass="Chemical">n, Supplemeclass="Chemical">ntary Table 4). Neither the metaproteome aclass="Chemical">nalysis class="Chemical">nor the cultivatioclass="Chemical">n approach detected aclass="Chemical">ny members of the Glomeromycota. The abseclass="Chemical">nce of these arbuscular mycorrhizal fuclass="Chemical">ngi is iclass="Chemical">n agreemeclass="Chemical">nt with the assumptioclass="Chemical">n that they are oclass="Chemical">nly fouclass="Chemical">nd iclass="Chemical">n soil layers aclass="Chemical">nd iclass="Chemical">n close associatioclass="Chemical">n with placlass="Chemical">nt roots (Domsch ).

Fungi are the main producers of litter-degrading enzymes

Our metaproteomics analyses enabled us to link structure and function of the complex microbial community present in the leaf litter. The obtained spectra were classified into COG (prokaryotic proteins) and KOG (eukaryotic proteins) categories based on their respective protein assignments (Supplementary Figure 2). The majority of bacterial protein spectra were assigned to functional categories such as translation, cell wall/membrane/envelope biogenesis, post-translational modifications/protein turnover and conversion, as well as energy production and conversion. Dominant functional groups of fungal proteins were translation, chromatin structure and dynamics, post-translational modifications/protein turnover and conversion and energy production and conversion. The proportion of proteins related to transport functions was higher in fungi than in bacteria (Supplementary Figure 2). When comparing samples from n class="Chemical">February aclass="Chemical">nd May a sigclass="Chemical">nificaclass="Chemical">nt iclass="Chemical">ncrease of proteiclass="Chemical">ns assigclass="Chemical">ned to post-traclass="Chemical">nslatioclass="Chemical">nal modificatioclass="Chemical">n aclass="Chemical">nd proteiclass="Chemical">n turclass="Chemical">nover was observed for both bacteria (t-test, t=−3.79; P=0.002) aclass="Chemical">nd fuclass="Chemical">ngi (t-test, t=−2.51; P=0.025).

As we were mainly interested in litter-decomposing enzymes, the metaproteome was searched for proteins that were assigned as extracellular hydrolytic enzymes. Moreover, total cellulolytic and class="Chemical">xylanolytic activities were measured by photometric assays. Strikiclass="Chemical">ngly, spectra could oclass="Chemical">nly be assigclass="Chemical">ned to eclass="Chemical">nzymes of fuclass="Chemical">ngal origiclass="Chemical">n, iclass="Chemical">ncludiclass="Chemical">ng phosphatases, cellulases, pecticlass="Chemical">nases, class="Chemical">n class="Chemical">xylanases, lipases, amylases, chitinases and oxidoreductases (Figure 3). The most prominent hydrolases were cellulolytic enzymes, that is, exo- and endo-glucanases as well as β-glucosidases.

Figure 3

Relative abundance of hydrolytic enzymes at different sampling times and sites. (a) Abundance of different enzymatic classes at different sampling times. Data represent mean±SD of eight samples collected at the respective sampling time. No significant differences among sites were detected by t-test analysis. (b) Abundance of different enzymatic classes at different sampling sites. Data represent mean±SD of four samples collected at the respective sampling site. *Represent significant differences, P<0.05, ANOVA followed by LSD-test.

The abundance of litter-degrading extracellular enzymes is affected by the sampling site but not by season

No significant changes were observed when the abundance of difclass="Chemical">fereclass="Chemical">nt fuclass="Chemical">ngal litter-degradiclass="Chemical">ng eclass="Chemical">nzyme classes was compared betweeclass="Chemical">n class="Chemical">n class="Chemical">February and May (Figure 3a). However, significant differences in enzyme abundances based on the NSAF between the sampling sites were observed (ANOVA, F=10.32; P<0.01) (Figure 3b), which was particularly evident for cellulases (ANOVA, F=20.24; P<0.01). Cellulases were most abundant in SW followed by KL, OR and AK. Furthermore, obvious trends in the abundance of phosphatases were detected between the sampling sites; fungal phosphatases were most abundant in AK followed by KL, OR and SW.

We did not observe any change in cellulolytic activity between class="Chemical">February aclass="Chemical">nd May, which is iclass="Chemical">n good agreemeclass="Chemical">nt with cellulolytic eclass="Chemical">nzyme abuclass="Chemical">ndaclass="Chemical">nce based oclass="Chemical">n spectral couclass="Chemical">nticlass="Chemical">ng (Figure 4). The comparisoclass="Chemical">n of cellulase abuclass="Chemical">ndaclass="Chemical">nce aclass="Chemical">nd total cellulolytic activity showed a sigclass="Chemical">nificaclass="Chemical">nt but weak positive correlatioclass="Chemical">n (SLR, r2=0.285, P=0.033). Iclass="Chemical">n coclass="Chemical">ntrast, class="Chemical">no correlatioclass="Chemical">n was observed wheclass="Chemical">n abuclass="Chemical">ndaclass="Chemical">nce aclass="Chemical">nd activity of class="Chemical">n class="Chemical">xylanase(s) were compared (SLR, r2=0.009, P=0.279). Whereas xylanase abundance did not differ between seasons; xylanase activity was higher in May.

Figure 4

Comparison of cellulase and xylanase abundance and total cellulolytic and xylanolytic activities. The upper part of the figure shows enzyme abundances based on NSAFs. Data represent the mean of two biological replicates. The lower part of the figure shows enzyme activities based on photometric assays. Data represent the mean±SD of three biological replicates.

The phylogenetic origin of cellulase-producing fungi depends on the sampling time and site

Cellulases represented the most abundant enzyme class involved in litter degradation at the investigated sampling times; therefore, their phylogenetic origin was more precisely investigated (Figure 5). Although fungal community changed only slightly, dramatic changes in the cellulase-producing fungal taxa were observed when samples from class="Chemical">February aclass="Chemical">nd May were compared. Iclass="Chemical">n class="Chemical">n class="Chemical">February, Leotiomycetes were the main cellulase producers, whereas in May a general trend to a broader spectrum of cellulase producers was observed. Eurotiomycetes' cellulases were most abundant in AK, whereas Dothideomycetes' (KL), Leotiomycetes' (OR) and Sordariomycetes' (SW) cellulases dominated the other sampling sites.

Figure 5

Phylogenetic origin of spectra assigned to fungal cellulases at different sampling times and sites compared with the overall distribution of fungal classes. Data represent the average of two biological replicates.

Community structure and function is influenced by season and leaf-litter nutrient content

SLR (Supplementary Figure 3) and PCA were performed to elucidate the influence of stoichiometry, season and leaf-litter nutrient content on community structure and functionality, as well as to identify the set of variables with the strongest relation to the decomposer community.

In our PCA biplots (Figures 6a, b and Supplementary Figures 4a–d), each axis (=principal component) is calculated from a set of variables, which together account for a high percentage of the variability of the data and can therefore be regarded as important drivers. The position of the sampling sites at the respective season within the biplot displays their relation to the measured variables. During litter decomposition, plant proteins decrease whereas fungal and bacterial proteins increase. On the x-axis (class="CellLine">PC1) iclass="Chemical">n Figure 6a, remaiclass="Chemical">niclass="Chemical">ng placlass="Chemical">nt proteiclass="Chemical">ns (‘Viridiplaclass="Chemical">ntae'), which caclass="Chemical">n be coclass="Chemical">nsidered as markers for the amouclass="Chemical">nt of remaiclass="Chemical">niclass="Chemical">ng litter material, are allocated oclass="Chemical">n the left/class="Chemical">negative side together with the remaiclass="Chemical">niclass="Chemical">ng C aclass="Chemical">nd Ca coclass="Chemical">nteclass="Chemical">nt. Bacterial proteiclass="Chemical">ns ordiclass="Chemical">nate together with pH, K aclass="Chemical">nd the microclass="Chemical">nutrieclass="Chemical">nts class="Chemical">n class="Chemical">Fe and Mn. Fungal proteins, temperature (T) and N content ordinate together. These variables are located opposite to the remaining plant proteins, indicating that they are related to a reduction in litter plant proteins. P content is also related to fungal growth as referred from the negative ordination of log P in the PCA (inversed direction results from log transformation). The y-axis (PC2) is dominated by the antagonistic variables precipitation and archaeal proteins. The positioning of sampling sites within the biplot and their shift from February to May towards the upper right reveals the effect of time/season factors on litter proteins: the reduction of plant proteins and the accumulation of (mainly) fungal proteins is positively correlated to warmer conditions. In addition, the individual sampling sites are separated according to their nutrient composition.

Figure 6

PCA biplots of community structure and function obtained in the metaproteomics approach as well as enviromental data and leaf-litter nutrient content. (a) General taxonomy–community structure, (b) community function. Data points are from different sampling times (February—squares and May—diamonds) and four sampling sites. Stars indicate transformations of the variables to meet the assumption of normal distribution, P is log 10 transformed and Ca is (1/(1+log10(Ca)) transformed.; T, mean air-temperature; *indicates inversion of the vector orientation because of data transformation.

SLR (Supplementary Figure 3) and PCA were also performed to identify factors afclass="Chemical">fecticlass="Chemical">ng eclass="Chemical">nzyme abuclass="Chemical">ndaclass="Chemical">nce aclass="Chemical">nd activity at the difclass="Chemical">n class="Chemical">ferent sampling sites. PCA of community function (represented by enzyme abundance) and time/season factors and nutrient content (Figure 6b) revealed relationships which support our findings on the community structure: the x-axis (PC1) can be regarded as ‘decomposition' axis with remaining, not yet respired, litter C content on the left/negative side and total enzymes, pH, K and micronutrients on the right side. Xylanases and cellulases—the main C-degrading enzymes—are ordinated opposite to litter C content. Also, Ca was inversely related to total enzyme abundance. Similarly, phosphatases were inversely related to P content. The y-axis (PC2) is dominated by season factors, such as WC and T, as well as by N content. Sampling sites are clearly separated along the x-axis due to their different nutrient and enzyme concentrations. The relation between time/season and nutrient release is indicated by an allocation of the May samples towards the warmer, drier conditions in the direction of the lower right in this plot. As enzyme abundances did not change with season, this movement is independent of enzyme ordination.

Interestingly, bacterial abundance correlated positively with the total abundance of extracellular hydrolytic enzymes (Supplementary Figure 3; r2=0.33, P=0.025) and in particular with cellulases (r2=0.44, P=0.007). However, no significant correlation between fungal abundance and extracellular hydrolytic enzymes was found (Supplementary Figure 3; r2=0.05, P=0.419 for cellulases; r2=0.02, P=0.619 for total enzymes). The bacterial relationship to cellulase abundance resulted in a negative correlation of the F/B ratio and cellulase abundance (Supplementary Figure 3; r2=0.42, P=0.002).

Discussion

Community structure of litter decomposers as revealed by metaproteomics

Metaproteomics can provide detailed information on the succession of the active part of the leaf litter inhabiting community (Nocker and Camper, 2009) and it thus enabled us to follow succession on a higher taxonomic resolution compared to e. g. class="Chemical">PLFA (for a detailed review see Joergeclass="Chemical">nseclass="Chemical">n aclass="Chemical">nd Wicherclass="Chemical">n, 2008). Ascomycota aclass="Chemical">nd Basidiomycota are commoclass="Chemical">nly regarded as predomiclass="Chemical">naclass="Chemical">nt fuclass="Chemical">ngal phyla iclass="Chemical">n the soil-litter iclass="Chemical">nterface (Osoclass="Chemical">no aclass="Chemical">nd Takeda, 2006). The metaproteomics data showed that the fuclass="Chemical">ngal commuclass="Chemical">nity was domiclass="Chemical">nated by Ascomycota aclass="Chemical">nd coclass="Chemical">ntaiclass="Chemical">ned oclass="Chemical">nly a small proportioclass="Chemical">n of Basidiomycota aclass="Chemical">nd Mucoromycoticlass="Chemical">na at all sampliclass="Chemical">ng-sites aclass="Chemical">nd -times (Figure 1d). This iclass="Chemical">n agreemeclass="Chemical">nt with the results of our cultivatioclass="Chemical">n approach (Supplemeclass="Chemical">ntary Iclass="Chemical">nformatioclass="Chemical">n). For fuclass="Chemical">ngal successioclass="Chemical">n it has beeclass="Chemical">n observed that Mucoromycoticlass="Chemical">na beloclass="Chemical">ng to the first coloclass="Chemical">nizers followed by Ascomycota (Torres ). These fuclass="Chemical">ngi have limited ability to degrade class="Chemical">n class="Chemical">lignin and are mainly regarded as cellulose decomposers or sugar fungi (Osono, 2007). Basidiomycota, with their ability to degrade the recalcitrant lignin-containing litter material, appear only later in the decomposition process (Osono, 2007; Lundell ). Successional analysis of the fungal community showed a decrease of Mucoromycotina and Ascomycota and an increase of Basidiomycota from February to May (Figure 1d). This is in good agreement with the observed fungal succession described by Torres and Osono (2007).

For bacterial succession, an increase in the proportion of Proteobacteria from winter to spring was observed, whereas that of class="Disease">Actinobacteria and Verrucumicrobia decreased (Figure 1b). Chaclass="Chemical">nges iclass="Chemical">n the respective group abuclass="Chemical">ndaclass="Chemical">nces were validated by a class="Chemical">n class="Chemical">PLFA analysis, which showed similar trends (Figure 2). A reduction of Actinobacteria was unexpected, because they are known to be involved in decomposition of organic materials, and thus are important for organic matter turnover and C cycle (Kirby, 2006). In other studies, an increase in the abundance of Actinobacteria has been shown during later stages of litter decomposition (Torres ; Snajdr ). The same accounts for the absence of Acidobacteria; members of this bacterial phylum can degrade various polysaccharides including cellulose and xylan (Ward ). Based on RNA sequencing, Baldrian found Acidobacteria to be the main bacterial group that was enriched in an active litter inhabiting community. The low number of acidobacterial proteins in our metaproteomics dataset might be explained by the limited number of Acidobacteria protein entries in the reference database (0.003% of database entries), which might have led to an underestimation of the contribution of Acidobacteria to the bacterial population.

Microbial succession is mainly influenced by leaf-litter quality

Previously, it has been shown that litter quality afn class="Chemical">fects decomposiclass="Chemical">ng commuclass="Chemical">nity aclass="Chemical">nd ecosystem processes (Prescott, 2010; Stricklaclass="Chemical">nd aclass="Chemical">nd Rousk, 2010). A major factor iclass="Chemical">n terrestrial ecosystems is pH (Siclass="Chemical">nsabaugh ). Our results revealed that pH had aclass="Chemical">n impact oclass="Chemical">n the bacterial, but class="Chemical">not oclass="Chemical">n the fuclass="Chemical">ngal commuclass="Chemical">nity, as seeclass="Chemical">n from a stroclass="Chemical">ng class="Chemical">negative correlatioclass="Chemical">n betweeclass="Chemical">n pH aclass="Chemical">nd the F/B ratio. Bacterial abuclass="Chemical">ndaclass="Chemical">nce was higher at more basic pH values (Figure 6a, Supplemeclass="Chemical">ntary Figures 3 aclass="Chemical">nd 4) as it has also beeclass="Chemical">n observed iclass="Chemical">n other studies (Fierer aclass="Chemical">nd Jacksoclass="Chemical">n, 2006; Högberg ). The lower impact of pH oclass="Chemical">n fuclass="Chemical">ngal diversity might be explaiclass="Chemical">ned by the fact that fuclass="Chemical">ngi caclass="Chemical">n grow over a wide pH raclass="Chemical">nge (Peclass="Chemical">nalva ), whereas bacterial growth is restricted to a smaller pH raclass="Chemical">nge as has beeclass="Chemical">n showclass="Chemical">n iclass="Chemical">n pure culture studies with isolates from class="Chemical">natural soils (Baath, 1996).

Ecosystem processes are often controlled by the availability of N or P (Elser ). Although litter N varied only slightly in our approach, clear difclass="Chemical">fereclass="Chemical">nces betweeclass="Chemical">n litter P values were preseclass="Chemical">nt (Table 3). The SLR aclass="Chemical">nd PCA aclass="Chemical">nalysis revealed that P is a major factor iclass="Chemical">nflueclass="Chemical">nciclass="Chemical">ng the abuclass="Chemical">ndaclass="Chemical">nces of all major fuclass="Chemical">ngal aclass="Chemical">nd bacterial groups with higher abuclass="Chemical">ndaclass="Chemical">nces at higher P (Figure 6a, Supplemeclass="Chemical">ntary Figures 3 aclass="Chemical">nd 4). class="Chemical">n class="Chemical">Phosphorous is needed for DNA replication and transcription and, therefore, fast growing bacteria are particularly affected under P-limiting conditions (Elser , 2003). Our data showing higher F/B ratios at nutrient and especially P poor sites (Supplementary Figure 3) are in good agreement with several earlier studies (Hieber and Gessner, 2002; van der Wal ; Güsewell and Gessner, 2009).

The function of fungi and bacteria in the initial phase of litter decomposition

Fungi are thought to be key players during litter decomposition because of their ability to degrade recalcitrant compounds such as class="Chemical">lignin aclass="Chemical">nd their domiclass="Chemical">naclass="Chemical">nce iclass="Chemical">n the decompositioclass="Chemical">n of class="Chemical">n class="Chemical">cellulose and hemicellulose (de Boer ; Meidute ). However, our knowledge about the contribution of bacteria and fungi to this process is still scarce. Recently, the analysis of ligninolytic and chitinolytic enzymes by a targeted metatranscriptomics approach of a forest soil (Kellner ) was performed with a focus on fungal enzymes, though neglecting the bacterial part of the community. Moreover, the analysis of a forest-soil metatranscriptome indicated that the active microbial community identified from RNA sequencing is dominated by Ascomycota and it is significantly different from the total community as indicated by DNA-based sequencing (Baldrian ).

Our analysis of the entire community revealed that only fungi produced extracellular enzymes in the investigated envclass="Chemical">ironmeclass="Chemical">nts, corroboraticlass="Chemical">ng with earlier ficlass="Chemical">ndiclass="Chemical">ngs that dealt with a secretome study of two model orgaclass="Chemical">nisms growiclass="Chemical">ng oclass="Chemical">n beech litter as substrate (Schclass="Chemical">neider ). Cellulases were the most abuclass="Chemical">ndaclass="Chemical">nt eclass="Chemical">nzymatic class. Iclass="Chemical">nteresticlass="Chemical">ngly, the phylogeclass="Chemical">netic origiclass="Chemical">n of cellulases chaclass="Chemical">nged over time (Figure 5) with fuclass="Chemical">ngal classes (Leotiomycetes, Dothideomycetes aclass="Chemical">nd Sordariomycetes of the Ascomycota) growiclass="Chemical">ng iclass="Chemical">n associatioclass="Chemical">n with leaves before litter fall as pathogeclass="Chemical">ns or eclass="Chemical">ndophytes domiclass="Chemical">naticlass="Chemical">ng cellulase productioclass="Chemical">n iclass="Chemical">n class="Chemical">n class="Chemical">February. In May, the contribution of ascomycotal classes like the Eurotiomycetes and basidiomycotal classes derived from soil and known for their saprotrophic lifestyle (Domsch ) increased. The generally low proportion of basidiomycotal cellulases might result from a low number of respective database entries. Only 21% of fungal cellulase entries are from Basidiomycota, whereas 73% belong to the Ascomycota and 6% to other fungal groups. Baldrian were able to increase the number of Basidiomycota-derived transcripts of the exocellulase gene cbh1 by sequencing of Basidiomycota species with a previously unknown genome, supporting the hypothesis that basidiomycotal cellulase entries are underrepresented in public protein databases.

No enzymes of bacterial origin were detected; this is likely not a result of a low number of bacterial cellulases in the reclass="Chemical">fereclass="Chemical">nce database (almost 60% of cellulase DB eclass="Chemical">ntries are of bacterial origiclass="Chemical">n), but might rather be explaiclass="Chemical">ned by a low abuclass="Chemical">ndaclass="Chemical">nce of bacterial eclass="Chemical">nzymes iclass="Chemical">ndicaticlass="Chemical">ng the miclass="Chemical">nor importaclass="Chemical">nce of these eclass="Chemical">nzymes iclass="Chemical">n the overall decompositioclass="Chemical">n process. Duriclass="Chemical">ng the iclass="Chemical">nvestigated stages of litter decompositioclass="Chemical">n, bacteria might proliclass="Chemical">n class="Chemical">ferate on low molecular weight carbohydrates provided by fungal enzymes. This idea is supported by the observed strong positive correlation between bacterial abundance and the abundance of fungal enzymes (Supplementary Figure 3) and might be referred to as ‘cheating behaviour' (Velicer, 2003). Bacteria might contribute to the production of extracellular degrading enzymes later in the decomposition process as shown for Actinobacteria, a group known for its ability to utilize lignin-derived compounds (Kirby, 2006). Moreover, it was shown that bacteria dominate litter decomposition in particular micro niches where anaerobic or high temperature conditions prevail (Lynd ).

Extracellular enzymes and their influence on decomposition

Our data indicated that the production of extracellular enzymes was influenced by litter stoichiometry. Cellulases and total enzyme abundances were highest at lowest C:P and C:N ratios (Supplementary Figure 3). C:P ratios correlated positively with the remaining plant proteins (Viridiplantae) and negatively with proteins assigned to many microbial groups (that is, Bacteroidetes, Ascomycota, Basidiomycota (Supplementary Figure 3), which shows that protein degradation and microbial community were strongly impacted by P. Sterner and Elser (2002) observed that the availability of P might be a rate-limiting factor in the synthesis of cellulolytic enzymes. The present study confirms this finding. Phosphatases ordinated opposite to P in our PCA (Figure 6b) supporting the idea that they are specifically produced when microbes have a certain demand for P (Sinsabaugh ).

The same holds true for the release of C from litter material. Difclass="Chemical">fereclass="Chemical">nt models for C-acquisitioclass="Chemical">n propose a sequeclass="Chemical">ntial decompositioclass="Chemical">n of class="Chemical">n class="Chemical">polysaccharides, starting with hemicellulose and cellulose degradation followed by the removal of lignin (Berg and Mcclaugherty, 2008; Snajdr ). Recently, Klotzbücher suggested that lignin is degraded preferentially in early phases of litter decomposition. In our approach, the dominance of xylanases and in particular cellulases indicates a decomposition phase in which hemicellulose and mainly cellulose can be considered as the major C-source. No insights into the degradation of highly recalcitrant phenolic lignin compounds were obtained because no lignolytic enzymes were detected. However, the increase of Basidiomycota, (known as main lignin decomposers) in May indicates that lignin degradation might start later in the decomposition process. Furthermore, our data imply that in litter with higher abundance of cellulose-degrading enzymes the amount of remaining degradable C is lower (Figure 6b) and the abundance of microorganisms is higher (Figure 6a, Supplementary Figure 3). This finding indicates that enzyme production might be regarded as the ‘bottleneck' of C-degradation/mineralization from litter. Similar hypotheses were postulated for soil organic-matter decomposition by Kuzyakov and Schimel and Weintraub (2003).

Conclusion

After litter fall, microorganisms grow on the litter using easily accessible compounds. Microbes are afclass="Chemical">fected by litter quality aclass="Chemical">nd eclass="Chemical">nvclass="Chemical">n class="Chemical">ironmental factors such as temperature and pH. The investigation of microbial succession revealed that Ascomycota dominated the fungal decomposer community with respect to abundance (total protein abundance) and activity (production of extracellular enzymes) especially in early stages of litter decomposition. In May, contribution of Basidiomycota to community functionality had increased in comparison with February. This strongly supports the assumption that the active part of microbial community changed over time. Only fungal, but no bacterial, polymer-degrading enzymes were detected at all sampling sites and -times, a finding that strongly supports our conclusion that fungi are the main players in litter decomposition at the investigated sites. Furthermore, a strong relationship between leaf-litter nutrient content and microbial abundance as well as cellulase production was found, indicating that high N and P content stimulate microbial decomposer activity. Accordingly, litter decomposition was fastest at sites with highest abundances of extracellular enzymes, suggesting that the production of extracellular enzymes is a rate-limiting step in the decomposition process.

We are fully aware that the metaproteomics-based approach sufclass="Chemical">fers from certaiclass="Chemical">n limitatioclass="Chemical">ns; the commuclass="Chemical">nity structure aclass="Chemical">nd fuclass="Chemical">nctioclass="Chemical">nality assessed oclass="Chemical">n the proteiclass="Chemical">n level might be iclass="Chemical">nflueclass="Chemical">nced by the proteiclass="Chemical">n extractioclass="Chemical">n efficieclass="Chemical">ncy, aclass="Chemical">nd hampered by false or missiclass="Chemical">ng assigclass="Chemical">nmeclass="Chemical">nts of MS aclass="Chemical">nd MS/MS data to proteiclass="Chemical">n eclass="Chemical">ntries iclass="Chemical">n the reclass="Chemical">n class="Chemical">ference database. Nevertheless, we believe that the metaproteome analysis, together with the validation by well-established complementary approaches, provides a deep insight into molecular details of the leaf-litter decomposition process. It opens up a new level of information in which both microbial succession and the activity of certain phylogenetic groups can be analyzed on the basis of their proteomes. The latter represent the active ‘building-blocks' in the ecosystem. Thus, further proteome analyses, for example, the investigation of later phases during litter decomposition or other litter types, will help to answer open questions in the field of litter decomposition: how the decomposer community further develops, who contributes to the actual degradation and which factors influence the process.