Source: https://www.nature.com/articles/s41396-018-0279-6?error=cookies_not_supported&code=4d4358f7-36c9-485c-87e6-81353f9966f1
Timestamp: 2019-04-23 08:58:33+00:00

Document:
Delignification, or lignin-modification, facilitates the decomposition of lignocellulose in woody plant biomass. The extant diversity of lignin-degrading bacteria and fungi is underestimated by culture-dependent methods, limiting our understanding of the functional and ecological traits of decomposers populations. Here, we describe the use of stable isotope probing (SIP) coupled with amplicon and shotgun metagenomics to identify and characterize the functional attributes of lignin, cellulose and hemicellulose-degrading fungi and bacteria in coniferous forest soils from across North America. We tested the extent to which catabolic genes partitioned among different decomposer taxa; the relative roles of bacteria and fungi, and whether taxa or catabolic genes correlated with variation in lignocellulolytic activity, measured as the total assimilation of 13C-label into DNA and phospholipid fatty acids. We found high overall bacterial degradation of our model lignin substrate, particularly by gram-negative bacteria (Comamonadaceae and Caulobacteraceae), while fungi were more prominent in cellulose-degradation. Very few taxa incorporated 13C-label from more than one lignocellulosic polymer, suggesting specialization among decomposers. Collectively, members of Caulobacteraceae could degrade all three lignocellulosic polymers, providing new evidence for their importance in lignocellulose degradation. Variation in lignin-degrading activity was better explained by microbial community properties, such as catabolic gene content and community structure, than cellulose-degrading activity. SIP significantly improved shotgun metagenome assembly resulting in the recovery of several high-quality draft metagenome-assembled genomes and over 7500 contigs containing unique clusters of carbohydrate-active genes. These results improve understanding of which organisms, conditions and corresponding functional genes contribute to lignocellulose decomposition.
The incomplete decomposition of woody biomass in coniferous forests is an important global carbon sink, with approximately one third of a gigaton of carbon accruing on an annual basis . Lignocellulose decomposition is influenced by the structure and function of microbial communities [2, 3] and terrestrial carbon cycling models increasingly parameterize these properties [4, 5]. However, efforts are constrained by our rudimentary knowledge of the composition and ecology of decomposer populations, stemming from limitations of culture-dependent methods and the complexity of soil communities. The best characterized decomposers inhabit forest soil litter, where conditions favour rapid lignocellulose degradation [6, 7], yet the decay of lignin-rich plant polymers in soil occurs in a continuum governed by conditions and substrate accessibility, resulting in diversified niches for decomposers [8,9,10]. To resolve the catabolic and ecological traits of decomposers, we must utilize culture-independent methods, like stable isotope probing (SIP), that better reflect in situ conditions.
The nature of microbial lignin-degradation is poorly described beyond the canonical breakdown of lignin in woody biomass by specialized, aerobic, litter-inhabiting wood-rot fungi. Yet, these fungi are largely absent in deeper mineral soil where conditions favour decomposition by bacteria , who represent the most likely lignin-degraders when oxygen is limited [12,13,14]. Several soil bacteria can degrade model lignin compounds in pure culture, suggesting a role in the catabolism of low-molecular weight, partially-degraded forms of lignin [15,16,17,18,19,20,21]. Knowledge of bacterial lignin-degradation in environmental contexts is limited, with all studies of soil degraders originating from tropical forests which describe active populations of predominantly Alpha and Gammaproteobacteria [13, 22]. The existing evidence for bacterial lignin-degradation demonstrates the need to characterize both bacterial and fungal populations and contrast their roles in different soil environments to better understand in situ the processes that govern the decomposition of lignocellulose.
Delignification rapidly increases the rate of lignocellulose decomposition and may have evolved primarily as a means of accessing more readily degradable plant carbohydrates, exemplified by the strategies of white and brown-rot Agaricomycota [23,24,25]. The capability to co-degrade lignin and other lignocellulosic polymers has not been thoroughly explored in bacteria. There is some evidence for the co-metabolism of lignocellulosic polymers by various Streptomyces spp. , while a comparative genomics study revealed that cellulose-degrading bacteria also possess higher numbers of hemicellulases . The extent to which catabolic traits are conserved within specific taxa versus communities will influence models of microbial decomposition. In the simplest case, the abundance of highly adapted, multi-substrate degrading taxa may predict rates of decomposition, which is supported by a small number of studies [2, 28]. However, forces of genomic streamlining in bacteria lead to functional diversification among closely related species, particularly in extra-cellular processes that produce common goods , evident in the species-level conservation of bacterial endoglucanases . A multi-substrate SIP experiment provides the means to identify whether forest soil decomposers can assimilate 13C from various lignocellulosic polymers and, with shotgun metagenomic sequencing, determine whether their genomes encode the necessary suite of catabolic genes.
To address the above knowledge gaps, we utilized SIP microcosm-based experiments to investigate the composition and degradative potential of hemicellulose, cellulose and lignin-degrading populations from organic and mineral layer soils in coniferous forests across North America. The identity of degraders and their genomic content were assessed using amplicon and shotgun sequencing of DNA enriched in 13C from labeled substrates. Lignocellulolytic activity was quantified according to the amount of 13C assimilated into total DNA and phospholipid fatty acids (PLFA). The in situ abundances of lignocellulolytic populations were determined on the basis of previously reported pyrotag libraries from the same field samples [31, 32]. This is currently the most comprehensive cultivation-independent study of lignocellulolytic populations in forest soils and yields new insights to the taxa responsible and the importance of certain catabolic gene families.
Soil samples were collected from five forest regions in North American chosen to encompass a range of climates, tree cover and soil conditions. Samples were collected from mature coniferous forest corresponding to ‘reference’ plots in the Long-Term Soil Productivity Study in northern Ontario (BSON and JPON), British Columbia (IDFBC), California (PPCA) and Texas (LPTX). Detailed descriptions of geographical locations, soil properties and sampling methods can be found in  as well as in Table S1. Samples were collected from three sites in each geographical region, except IDFBC where three plots were sampled from a single site. The organic layer (O-horizon) and the top 20 cm of mineral soil (A-horizon + occasionally upper B-horizon) were sampled separately (overview in Fig. 1a). Samples were kept at 4 °C during transport, sieved through 2-mm mesh and stored at −80 °C.
Microcosms were prepared with a minimum of three replicates by adding 13C-labeled hemicellulose, cellulose or lignin to soil from as many of the five regions as possible contingent on the availability of substrate (Fig. 1a). Microcosms were comprised of either 1 g (organic) or 2 g (mineral) dry wt soil in 30-mL serum vials adjusted to a moisture content of 60% (mineral) and 125% (organic) (w/v). Microcosms were wetted and pre-incubated in the dark at 20 °C for one-week prior to the addition of 1% (w/w) 13C-labeled cellulose or 0.8% (w/w) lignin. The hemicellulose-amended incubations were performed as part of a previous study examining the impacts of timber harvesting on decomposition using identical methods as described here for cellulose . Each microcosm with 13C-labeled substrate was paired with an identical ‘12C-control’ microcosm amended with the corresponding unlabeled substrate (~1.1% atom 13C). 12C-control microcosms were used to estimate the natural abundance of 13C in PLFAs and to control for the background presence of GC-rich DNA in higher density CsCl gradient fractions .
Bacterial cellulose was produced from Gluconacetobacter xylinus str. KCCM 10100 grown with either unlabeled or 13C-labeled glucose (99 atom % 13C, Cambridge Isotope Laboratories, MA, USA) in Yamanaka medium as described in the Supplementary Methods. DHP lignin was synthesized from ring-labeled coniferyl alcohol using horseradish peroxidase as described . Coniferyl alcohol was synthesized from ring-labeled vanillin (75 atom % 13C, Sigma Aldrich, CA) according to the reactions in Figure S1 with methods outlined in the Supplementary Methods. The DHP lignin had a weight average molecular weight (Mw) of 2,624 g·mol−1 which equates to ~14 polymerized units of coniferyl alcohol. Yields for cellulose and lignin were approximately 12% and 18% (w/w) and were labeled at ~99 and 60% atom % 13C, respectively.
Microcosms were incubated at 20 °C for 14 days (cellulose) or 60 days (lignin) or 2 days (hemicellulose; ). The incubation length for this study was optimized by preliminary time-course experiments described in . To overcome poor 13C-labeling of microbial biomass in organic soils from DHP-lignin, organic soils were mixed 1:3 with corresponding double-autoclaved mineral soil to dilute pre-existing organic matter. Following incubation, soil samples were lyophilized and stored at −80 °C until processing. All SIP-PLFA, SIP-pyrotag and SIP-metagenomic data were generated from the same microcosm unless otherwise stated. Due to limited quantities of 13C-substrate, SIP-lignin microcosms were prepared with soil from only BSON, PPCA and LPTX, while SIP-hemicellulose data was available for IDFBC and PPCA (see Table S2 for summary of data). To determine the lignolytic activity solely by bacteria, eight additional soil microcosms were prepared with 13C-lignin amended mineral soil and incubated with or without weekly additions of two fungicides: 1 mg·g−1 soil cylcoheximide plus 1.2 mg·g−1 soil fungizone (amphotericin B).
PLFAs were extracted from 0.75 g (organic) or 1.0 g (mineral) dry wt soil according to  and 13C-content was analyzed using ion ratio mass spectrometry (UBC Stable Isotope Facility) ported with gas chromatography as detailed in . Peak identification was based on retention time compared against two reference standards: bacterial acid methyl-ester standard (47080-0; Sigma–Aldrich, St. Louis) and a 37-Component fatty acid methyl-ester mix (47885-U; Sigma–Aldrich, St. Louis). For details on the treatment of unidentifiable, but 13C-labeled, peaks and taxonomic designations consult the Supplementary Methods.
DNA was extracted from soil (0.5 g) using the FastDNA™ Spin Kit for Soil (MPBio, Santa Ana, CA) following the manufacturer’s protocol. 13C-enriched DNA was separated using cesium chloride density centrifugation per methods outlined by  with minor modifications (see Supplementary Methods). DNA from heavy fractions (1.727–1.735 g·mL−1; typically F1-F7) was pooled, desalted and concentrated using Amicon Ultra-0.5 mL filters (EMD Millipore, MA, USA). DNA from 12C-controls was treated identically; however, the total DNA recovered from those heavy fractions of 12C-controls was insufficient for sequencing library preparation, resulting in the need to pool additional fractions (F1-F9). The % atom 13C of soil DNA extracts and DNA from heavy fractions was quantified according to . All other DNA quantitation was performed using Pico-Green fluorescent dye (ThermoFisher, MA, USA).
PCR amplification was performed on bacterial 16 S rRNA gene (V1–V3) using barcoded primers as described in . PCR reactions were performed in triplicate and pooled prior to purification. Samples were sequenced using the Roche 454 Titanium platform (GS FLX+) at the McGill University and Genome Québec Innovation Centre, yielding an average of 8400 bacterial reads per sample with a minimum length of 250-bp following quality processing. All 16 S rRNA gene libraries are stored at the European Nucleotide Archive (ENA) under the study accession PRJEB12502. Identically prepared 16 S rRNA gene amplicon libraries for SIP-hemicellulose  and for corresponding unincubated field soil samples  were re-analyzed in the present study, following retrieval from the ENA.
Whole shotgun metagenomes were prepared from 13C-enriched or 12C-control DNA using either the Nextera DNA Sample Preparation Kit (Illumina Inc., CA, USA) for SIP-cellulose microcosm (40–50 ng template DNA), or the Nextera XT DNA Sample Preparation Kit (Illumina Inc., CA, USA) for SIP-lignin (1 ng), due to lower recovery of 13C-enriched DNA. The two sequencing library preparations produced very similar metagenomic profiles in a direct comparison . Four libraries were multiplexed per lane of Illumina HiSeq 2500 (2 × 100-bp) and were grouped based on similar average fragment sizes based on fragment profiles from a 2100 Bioanalyzer (Agilent Technologies) All raw and assembled shotgun metagenomic data is available under ENA accession PRJEB12502, with an overview of accessioned data available in Table S3.
16S rRNA gene amplicon libraries for all three SIP-substrates were quality filtered and simultaneously processed using mothur according to the Schloss '454 SOP' (accessed November 2015; ) and clustered into operational taxonomic units (OTUs) at 0.01% dissimilarity. Taxonomic classification was performed using the RDP Classifier  with the GreenGenes database for 16S rRNA genes (database gg_13_8_99; August 2013). All OTU counts were normalized to total counts per thousand reads. All related count tables, classifications and sample data is provided as phyloseq objects in the Supplementary Data. Phylogenetic trees were prepared by placing sequences into the pre-aligned SILVA tree (; ‘SSURef_NR99_123_SILVA_12_07_15′) using maximum parsimony in ARB .
Shotgun metagenome libraries from SIP-cellulose and SIP-lignin were quality preprocessed with Trimmomatic (; v. 0.32) and FastX Toolkit (; v. 0.7). Individual metagenomes were assembled using Ray-meta . Metagenomes were also composited by location and soil layer and assembled with the low-memory assembler MEGAHIT (; v. 1.0.2). Contigs from composited metagenomes were binned into metagenome-assembled genomes (MAGs) using MetaBAT (; v. 0.18.6) based on read mapping with Bowtie2 . The quality of draft genome bins was assessed using CheckM  and the lowest common ancestor classification of raw, assembled and MAG data was performed with MEGAN (; v. 5.10.1) based on DIAMOND blastx (; v. 0.7.9) searches against a local version of the ‘nr’ NCBI database (downloaded October 2014). MAGs were classified to a specific taxon if >25% of contigs were uniformly classified at the genus level. ORF prediction was performed by Prodigal (; v. 2.6.2) and CAZy encoding ORFs were identified based on BLAST searches against a local CAZy database (downloaded 19 August 2015) and hmmscan searches (; HMMER v.3.1b1) using custom hmms for oxidative, lignolytic gene families (see Supplementary Methods) and hmms provided by dbCAN . Contigs containing clusters of three or more genes encoding CAZymes (‘CAZy gene clusters’) were recovered using a custom approach described in the Supplementary Methods. CAZyme families attributed xylanase, endoglucanase and ligninase activity can be found in the Supplementary Methods.
Statistics were performed using R v. 3.1.0  with general dependency on the following packages: reshape2, ggplot2, plyr [59,60,61], Hmisc  and phyloseq . Where necessary, P-values were adjusted according to Benjamini and Hochberg . The ‘vegan’ R-package  provided tools to calculate, non-parametric multidimensional scaling (metaMDS) on Bray-Curtis dissimilarities (‘vegdist). A combination of DESeq  and indicator species analysis (De Cáceres and Legendre 2009) were used to compare control and labeled microcosms to identify OTUs and taxonomic groups that incorporated 13C from the labeled substrate. An OTU or taxon had to be selected by at least one of these methods and be, on average, 3-fold more abundant in 13C libraries in at least one geographic region. Random forest classification, implemented in boruta , was used to identify taxa or CAZymes predictive of total 13C incorporation into DNA or PLFA. Correlation between features selected by boruta and total 13C were subsequently performed using ‘rcorr’ from the R package Hmisc. The relative importance of soil layer, geographic region, total carbon, CAZy composition and community structure in linear models for each substrate was assessed using the R package relaimpo  using the primary and secondary axes of NMDS as measures of community structure and CAZy composition.
We utilized SIP microcosm-based experiments to investigate the taxonomic identity and catabolic potential of lignocellulolytic populations from coniferous forest soils in five regions across North America. Microcosms contained either organic (top 2–5 cm) or mineral layer soil (5–20 cm) and were amended with either 13C-labeled hemicellulose, cellulose or lignin (Fig. 1a). Significant enrichment of 13C in soil DNA extracts (Fig. 1b) and PLFAs, and the recovery of heavy DNA from higher density gradient fractions (Fig. 1c) confirmed the catabolism and assimilation of 13C from substrates. The quantification of total 13C in PLFA and DNA was in close agreement (r = 0.74, p < 0.001). The successful targeting of lignocellulolytic populations was evident in the strong compositional differences in DNA and PLFA profiles from soils amended with 13C-labeled substrates compared to controls. A total of 503 OTUs were enriched in 13C-amplicon libraries, enabling the determination of which taxa had the capacity to degrade multiple lignocellulosic (Fig. 2; full list in Table S4). The enrichment of degradative populations was also apparent in shotgun metagenomes, as evidence by the increased proportion of reads comprising the assembly (Figure S2), the assembly quality (N5013C = 4000 vs. N5012C = 2100; t = 2.83, p = 0.005), and the recovery of a greater number of MAGs from metagenomes prepared from 13C-enriched DNA compared to controls (Fig. 3). The high quality of MAGs enabled the characterization of gene content to test our conclusions about individual capacities for catabolizing lignocellulosic polymers.
Members of the family Caulobacteraceae and order Burkholderiales dominated 13C-enriched DNA pools, comprising between 1–5% of 13C-metagenomes, and were highly represented in our collection of MAGs (Fig. 3). In Burkholderiales, three major families possessed taxa that could degrade one or more substrates: Burkholderiaceae, Comamonadaceae and Oxalobacteraceae (Fig. 2; Figure S3). In Caulobacteraceae, members of the genus Asticcacaulis assimilated 13C from both hemicellulose and cellulose, while members of Caulobacter demonstrated the capacity to assimilate all three substrates (Fig. 4a; full breakdown in Figure S4). Notably, several OTUs belonging to unclassified clades of Caulobacteraceae were highly enriched in 13C-DNA pools from lignin (Figure S5) and were significantly more abundant in mineral than organic soils (21.5 vs. 6.0 counts per thousand reads, respectively; t-test; p = 0.02). Differences in predisposition for cellulose- versus lignin-catabolism in Asticcacaulis versus Caulobacter, respectively, was evident in the catabolic gene content of metagenome assemblies (Figs. 4b, c).
We anticipated that lignocellulolytic species would commonly degrade more than one of the three lignocellulosic polymers tested, yet most incorporated 13C from only one (Fig. 2). No single OTU had increased relative abundance in 13C-DNA pools for all three substrates. Only five OTUs were enriched in 13C-DNA pools for two substrates and were, in all cases, enriched from 13C-lignin and 13C-cellulose. These OTUs were classified as Simplicispira (Burkholderiaceae), Aquincola (Comamonadaceae), unclassified Caulobacteraceae and two Sinobacteraceae (Figure S6). The only genera to possess OTUs that collectively incorporated 13C from all three substrates were Pelomonas and Caulobacter. More taxa possessed OTUs that collectively catabolized a combination of hemicellulose and cellulose (n = 7) than cellulose and lignin (n = 4), or hemicellulose and lignin (n = 3). Genera that utilized both hemicellulose and cellulose included Cellvibrio, Janthinobacterium, Cytophaga and Salinibacterium. All fungal taxa enriched in unassembled 13C-metagenomes were enriched from cellulose, with the sole exception of a putative member of Saccharomycetaceae from lignin. SIP tended to enrich relatively rare populations, evidenced by the small proportion of OTUs (41/503) shared between 13C-amplicon libraries and unincubated field samples representing in situ relative abundances.
Variation in the composition of 13C-assimilating bacterial populations was explained most by substrate use (PERMANOVA; R2 = 13%; p < 0.001), followed by the interaction of geographic region and substrate use (10%; p < 0.001), region (8%; p < 0.001) and the interaction of substrate and site (6%; p < 0.001; Table S6). The relative importance of the predictors of 13C-incorporation into DNA or PLFA differed for hemicellulose, cellulose and lignin-degrading populations (Fig. 5a). However, soil layer, forest region and CAZy content consistently ranked highly followed by taxonomic composition and total soil carbon content. Higher soil organic matter content slowed rates of 13C-assimilation, with total carbon and nitrogen negatively correlated with 13C enrichment of PLFAs in hemicellulose (ρ =−0.71; padj < 0.01) and lignin-amended mineral layer soils (ρ = −0.85; padj < 0.001; Fig. 5c), but not those with cellulose. The highest rates of 13C-assimilation were observed for hemicellulose and cellulose in the organic layer of IDFBC and for both cellulose and lignin in LPTX mineral soil (Fig. 6a).
Highly active cellulolytic populations in the IDFBC organic layer were uniquely dominated by bacteria in contrast to those from other regions and soil layers (Fig. 6b). These active cellulolytic populations were comprised largely of taxa not active elsewhere, namely: MIZ46 (Deltaproteobacteria), Cellvibrio (Gammaproteobacteria), and Planctomyces (vadinHA49) (Figure S7). In contrast, Ascomycota were the major degraders of cellulose in most other soils (i.e., BSON, PPCA and LPTX). The active lignolytic populations in BSON soils were also distinct from other forest soils, and were comprised of Cystobacteraceae (Deltaproteobacteria), FAC88 (Elusimicrobia, formerly Termite Group 1), Gaiellaceae and Solirubrobacterales (Actinobacteria). Twelve bacterial genera were identified by random forest classification whose relative abundance in 13C amplicon libraries was significantly positively correlated with 13C-enrichment of DNA (Table S7). The strongest correlations occurred for Caulobacter (cellulose; ρ = 0.6; padj < 0.001), unclassified Caulobacteraceae (lignin; Fig. 5b), and Janthinobacterium (cellulose; ρ = 0.59; padj < 0.001). No fungi were correlated with the 13C-enrichment of DNA or PLFA.
In most microcosms, bacteria were predominantly active in catabolizing and assimilating 13C from substrates and potentially other breakdown products (Fig. 6b; Figure S8). Gram-negative bacteria had the highest δ-13C enrichment of PLFAs in 72% of lignin-amended microcosms, while roughly equal numbers of cellulose- amended microcosms were dominated by gram-negative bacteria or fungi. The total quantity of 13C assimilated into PLFAs (µmol·g−1 dry wt soil) did not significantly differ between microcosms with greater bacterial versus fungal assimilation (Wilcoxon test; W = 194, p = 0.7). Gram-negative and gram-positive bacterial populations dominated the assimilation of hemicellulose in a comparable number of microcosms, 41% and 39%, respectively, though 13C-enrichment was 2-fold higher when gram-negative bacteria dominated (Wilcoxon test; W = 1261, p < 0.001).
Fungicide-treatment greatly reduced assimilation of 13C from lignin into fungal PLFAs, but also into gram-positive bacteria and, to a lesser extent, gram-negative bacteria (Fig. 7a). While overall bacterial activity was decreased by fungicide-treatment, the assimilation of 13C by gram-negative bacteria persisted at relatively high levels (20–50% of untreated incubations) when zero enrichment of fungal PLFAs had occurred (in BSON and PPCA). The major taxa incorporating 13C-lignin in fungicide-treated incubations were distinct from those in untreated ones and were predominantly gram-negative bacteria: Burkholderiaceae and Sphingobacteriaceae (Fig. 7b). Despite the reduction in 13C-enrichment of gram-positive PLFAs, the relative abundance of a major lignolytic Actinobacterial family (Nocardiaceae) and Actinobacteria overall (not shown) were increased in 13C-metagenomes after fungicide treatment. The impact of fungicide-treatment on the relative abundance of Caulobacteraceae differed by region, with only minor decreases observed in PPCA and LPTX.
Genes encoding lignin-modifying auxiliary activity (AA) enzymes were approximately three-fold more abundant in 13C-lignin versus 13C-cellulose metagenomes. Yet, many AA families were also enriched in 13C-cellulose versus control metagenomes, which corresponded to an overall increase in reads classified as fungal (Fig. 8). In contrast, the enriched AA gene families in 13C-lignin metagenomes (AA3, AA4 and AA6) were largely bacterial in origin. Overall, there was an 18-fold increase in fungal versus bacterial AA reads in 13C-cellulose metagenomes compared to a 4-fold increase in bacterial versus fungal AA reads in 13C-lignin metagenomes (Mann–Whitney, W = 92, p < 0.001).
Genes encoding several oxidative enzymes implicated in the depolymerization of lignin were enriched in 13C-lignin metagenomes and were classified to taxa previously identified as assimilating 13C from lignin in amplicon libraries. More specifically, genes encoding laccases (AA1), dye-decolouring peroxidases (DyP) and alcohol oxidases were the most enriched AA gene families. These were predominantly classified to Commamonadaceae, Sphingomonadaceae and Caulobacteraceae (Figure S9a). Nine MAGs contained genes predicted to encode the entire β-ketoadipate pathway, potentially involved in catabolizing aromatic lignin depolymerization products to central metabolites, and were classified to Burkholderia, Sphingomonas, Caulobacter and Sorangium. Aryl alcohol oxidase genes were the most significantly correlated (ρ > 0.3 and padj < 0.01) with total 13C-assimilation from lignin into PLFA.
13C-cellulose metagenomes were highly enriched in glycosyl hydrolase gene families (GHs). Approximately 37% of CAZy gene clusters (groups of 3 or more GHs) from 13C-cellulose metagenomes contained carbohydrate-binding module (CBM) genes and 4% of these also contained an endoglucanase gene. In comparison, only 22% of clusters in 13C-lignin metagenomes contained a CBM gene, and none of these contained endoglucanase genes. Thirty-three percent of endoglucanase gene-containing clusters contained a CBM, and were mostly classified to Asticcacaulis, Sorangium and Chthoniobacter (Figure S9b). Actinomycetales had the greatest diversity of endoglucanase gene families in CAZy gene clusters, while fungi predominantly had endoglucanase genes belonging to GH7 and GH9 families (Figure S9c). Few endoglucanase or GH gene families were correlated with cellulolytic activity (Table S7). Four MAGs contained CAZy gene clusters encoding enzymes putatively catabolizing all three substrates (i.e., xylanases, endoglucanases and ligninases) and were classified as Chaetomium (fungi), Caulobacter, Chthoniobacter and Sphingobium (Table S8).
We failed to find support for our initial hypothesis that individual species would degrade multiple components of lignocellulose. Instead, we identified a diverse array of taxa assimilating 13C from three major lignocellulosic polymers, demonstrating the degree of specialization in decomposer populations. While approximately one third of taxa enriched in 13C-DNA pools possessed members that could collectively degrade more than one substrate, and, though a handful of MAGs encoded a putative suite of catabolic enzymes, none of the individual bacterial OTUs (at a 99% similarity threshold) were definitively enriched from all three substrates. Even closely related taxa, like Caulobacter and Asticcacaulis, exhibited differences in their levels of enrichment on cellulose and lignin and in CAZy gene content. In contrast, the multi-substrate degradative capacity of fungi was evident in the increased relative abundances of genes encoding both delignifying (peroxidases) and cellulolytic enzymes in 13C-DNA pools. The recovery of a Myceliophthora MAG from 13C-metagenomes was consistent with its known capability for complete lignocellulose decomposition [69, 70]. Therefore, our findings suggest that complementation among functional guilds may be necessary for decomposition of lignocellulose by bacteria and not necessarily fungi. One caveat is the possibility that our method was biased towards identifying single substrate utilizing bacteria. The activity of multi-substrate users may have been masked by populations that grew faster on individual substrates, particularly the relatively labile hemicellulose, as was observed in other SIP experiments . At the very least, our study demonstrates the decomposition of lignocellulosic polymers can commonly occur via a division of labour among specialized taxa.
The substantial assimilation of carbon from our model lignin substrate by bacteria supports our hypothesis that bacteria contribute significantly to degradation of native forms of lignin in situ. Bacterial activity was particularly evident in deeper mineral layers of forest soil, suggesting that lignin decomposition can occur throughout the soil column. Interestingly, the most novel lignin-degraders from mineral soils belong to uncultured clades of Caulobacteraceae, Acidobacteria, Solirubrobacterales, Elusimicrobia, Nevskiales, and Cystobacteraceae. The low activity of fungal lignin degraders in our microcosms was unlikely due to unmet nutritional or environmental conditions, since metabolically competent taxa like Myceliophthora [69, 70] assimilated carbon from cellulose under the same conditions. One possible explanation is that fungi were involved in lignin-degradation but did not metabolize degradation products to become sufficiently labeled, perhaps because of bacterial cross-feeding, evident in the inhibitory effect of fungicide on lignin incorporation by bacteria. Indeed, the results from the fungicide-treatments indicated significant interactions occurred between bacteria and fungi in lignin degradation. However, fungi were not essential for bacterial lignin degradation, given the persistent and substantial incorporation of 13C-lignin by bacteria in fungicide-treated soils. This study provides strong evidence for the capacity of bacteria to degrade lignin, highlighting the need for further characterization of their activity and ecology in soil.
The substantial number of taxa identified by SIP with previously reported lignocellulolytic activity was testament to the success of previous culturing-based efforts to characterize soil decomposers and a validation of our SIP approach. Approximately 72% (14/19) of hemicellulolytic taxa and 74% (31/42) of cellulolytic taxa were previously reported to degrade the corresponding substrate in vitro. Far fewer lignolytic taxa (~28%, 8/29) had previously reported degradative activity, likely due to less frequent and formalized testing of lignin degradation. However, there was considerable agreement in the bacterial genera we identified and those identified in other culture-independent studies of lignin-degraders in forest soil, such as Caulobacter, Sphingomonas, Sphingobacterium, Nocardia, Telmatospirillum and Azospirillum [13, 22, 72]. A lower proportion of lignin-degrading taxa associated with mineral layer soil (5/16) had previously reported activity compared to those associated with the organic layer (8/12), suggesting the former may be more difficult to culture or were less frequently targeted for study.
Many of the novel cellulolytic groups we identified belong to cultivation-resistant phyla, Planctomycetes, Verrucomicrobia, Chloroflexi and Armatimonadetes (formerly OP10), which commonly predominate soil communities [73, 74]. Each of these phyla possesses at least one isolate capable of degrading cellulose [75,76,77,78] and have been designated cellulolytic in other SIP-cellulose experiments [28, 71, 79, 80]. Notably, all contigs that contained clusters of ten or more CAZymes (27 contigs), from 13C-cellulose metagenomes, were classified to taxa from the aforementioned groups and contained both CBMs (26/27) and endoglucanases (20/27). These findings supported our assertion that the richness of described lignocellulolytic taxa is underestimated in soils, which limits our understanding of the ecology of decomposition and may afford new types of biocatalysts for processing lignocellulosic biomass.
The widespread capacity among members of Caulobacteraceae to degrade the three components of lignocellulose was unexpected, yet consistent with reports of their enrichment on decaying wood [81, 82] and during early stages of litter decomposition . Caulobacter were first isolated from cellulose-amended lake-water , but are primarily known as oligotrophic, aquatic organisms . Their role in degradation of plant carbohydrates was first postulated based on the analysis of the C. crescentus genome  and subsequently demonstrated by growth on cellulose, xylose and vanillate [87,88,89,90]. Although we provide the first evidence for the role of Caulobacteraceae in degrading all three polymers of lignocellulose, several surveys of forest soils report enrichment of Caulobacter or Asticcacaulis in samples amended with cellulose [80, 91, 71] or lignin [13, 22]. The role of Caulobacteraceae in decomposition in forest soils may prove significant, given their relatively high abundance (0.5–2.5% of total libraries) and their capacity to adhere to insoluble polymers, like lignocellulose.
The taxonomy and catabolic capacity of lignocellulose degraders were moderately associated with variation in the total lignocellulolytic activity. The influence of compositional differences between soil layers was, in certain cases, overshadowed by the strong negative correlation between organic matter and 13C enrichment. Yet, in regression modeling, CAZy gene content and community composition had equivalent or greater explanatory power than other predictors. The abundances of several prominent lignocellulolytic taxa were significantly correlated with higher 13C assimilation, demonstrating that certain taxa play more significant roles, at least under conditions we tested. Strongly oxidative bacterial enzymes commonly studied in relation to lignin-degradation, such as DyP-type peroxidases and laccases [92,93,94], were less predictive of lignolytic activity than aryl alcohol oxidases, suggesting a greater role of the latter class of enzymes in bacterial lignin-degradation. Overall, fewer endoglucanase genes were predictive of 13C-assimilation compared to AA gene families. This difference suggests either the existence of a greater diversity of endoglucanases, minimizing the explanatory power of any single gene, or a narrower activity range of AA families with a correspondingly higher relevance to lignolytic activity. This result agrees with the finding that the composition of decomposer communities had increased explanatory power for the rate of litter decomposition during later stages when greater proportions of lignin remained (Cleveland et al. 2014).
The comprehensive study of all three major lignocellulosic polymers enabled us to examine the co-occurrence of lignocellulolytic traits, addressing several knowledge gaps in this area. The results indicated unexpected specialization among bacterial populations for degradation of individual lignocellulosic polymers and revealed several novel lignocellulolytic taxa, highlighting current limits to our knowledge of decomposition. This research supports the view that bacterial decomposition of oligomeric lignin is a ubiquitous soil process, with the potential to occur in deeper soil layers following early stages of litter decomposition. As hypothesized, variation in community composition was found to constrain lignocellulolytic activity across various forest types and soil layers in North America. The relationship between these communities and process rates should receive continuing study to refine our understanding of soil carbon stabilization and terrestrial carbon cycling models. Furthermore, the large number of degradative gene clusters from uncultivated lignocellulolytic taxa represent a trove of potentially novel enzymes for biotechnological applications.
Myneni RB, Dong J, Tucker CJ, Kaufmann RK, Liski J, Zhou L, et al. A large carbon sink in the woody biomass of Northern forests. PNAS. 2001;98:14784–9.
Strickland MS, Rousk J. Considering fungal: bacterial dominance in soils–methods, controls, and ecosystem implications. Soil Biol Biochem. 2010;42:1385–95.
Cleveland CC, Reed SC, Keller AB, Nemergut DR, O’Neill SP, Ostertag R, et al. Litter quality versus soil microbial community controls over decomposition: A quantitative analysis. Oecologia 2014;174:283–294.
Allison SD. A trait-based approach for modelling microbial litter decomposition. Ecol Lett. 2012;15:1058–70.
Wieder WR, Grandy AS, Kallenbach CM, Taylor PG, Bonan GB. Representing life in the Earth system with soil microbial functional traits in the MIMICS model. Geosci Model Dev. 2015;8:1789–808.
Melillo JM, Aber JD, Linkins AE, Ricca A, Fry B, Nadelhoffer KJ. Carbon and nitrogen dynamics along the decay continuum: Plant litter to soil organic matter. Ecol Arab L — Perspect Challenges 1989;198:53–62.
Cotrufo MF, Soong JL, Horton AJ, Campbell EE, Haddix ML, Wall DH, et al. Formation of soil organic matter via biochemical and physical pathways of litter mass loss. Nat Geosci 2015;8:776–9.
von Lützow M, Kögel-Knabner I, Ekschmitt K, Flessa H, Guggenberger G, Matzner E, et al. SOM fractionation methods: Relevance to functional pools and to stabilization mechanisms. Soil Biol Biochem 2007;39:2183–207.
Lehmann J, Kleber M. The contentious nature of soil organic matter. Nature 2015;528:60–8.
Klotzbücher T, Kaiser K, Guggenberger G, Gatzek C, Kalbitz K. A new conceptual model for the fate of lignin in decomposing plant litter. Ecology 2011;92:1052–62.
Ekschmitt K, Kandeler E, Poll C, Brune A, Buscot F, Friedrich M, et al. Soil-carbon preservation through habitat constraints and biological limitations on decomposer activity. J Plant Nutr Soil Sci 2008;171:27–35.
Benner R, Newell SY, Maccubbin AE, Hodson RE. Relative contributions of bacteria and fungi to rates of degradation of lignocellulosic detritus in salt-marsh sediments. Appl Environ Microbiol 1984;48:36–40.
Hall SJ, Silver WL, Timokhin VI, Hammel KE. Lignin decomposition is sustained under fluctuating redox conditions in humid tropical forest soils. Glob Change Biol. 2015;21:2818–28.
Brown ME, Barros T, Chang MC. Identification and characterization of a multifunctional dye peroxidase from a lignin-reactive bacterium. ACS Chem Biol. 2012;7:2074–81.
Crawford DL. Lignocellulose decomposition by selected Streptomyces strains. Appl Environ Microbiol. 1978;35:1041–5.
Davis JR, Goodwin L, Teshima H, Detter C, Tapia R, Han C, et al. Genome sequence of Streptomyces viridosporus strain T7A ATCC 39115, a lignin-degrading Actinomycete. Genome Announc. 2013;1:e00416–00413.
DeAngelisb KM, D’Haeseleer P, Chivian D, Fortney JL, Khudyakov J, Simmons B, et al. Complete genome sequence of Enterobacter lignolyticus SCF1. Stand Genom Sci. 2011;5:69–85.
Masai E, Katayama Y, Fukuda M. Genetic and biochemical investigations on bacterial catabolic pathways for lignin-derived aromatic compounds. Biosci Biotechnol Biochem. 2007;71:1–15.
Rashid GM, Taylor CR, Liu Y, Zhang X, Rea D, Fulop V, et al. Identification of manganese superoxide dismutase from Sphingobacterium sp. T2 as a novel bacterial enzyme for lignin oxidation. ACS Chem Biol. 2015;10:2286–94.
Taylor CR, Hardiman EM, Ahmad M, Sainsbury PD, Norris PR, Bugg TDH. Isolation of bacterial strains able to metabolize lignin from screening of environmental samples. J Appl Microbiol. 2012;113:521–30.
Pold G, Melillo JM, DeAngelis KM. Two decades of warming increases diversity of a potentially lignolytic bacterial community. Front Microbiol. 2015;6:480.
Eastwood DC, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A, et al. The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. Science (80-) 2011;333:762–5.
Ding SY, Liu YS, Zeng Y, Himmel ME, Baker JO, Bayer EA. How does plant cell wall nanoscale architecture correlate with enzymatic digestibility? Science. 2012;338:1055–60.
Zeng Y, Zhao S, Yang S, Ding SY. Lignin plays a negative role in the biochemical process for producing lignocellulosic biofuels. Curr Opin Biotech. 2014;27:38–45.
Větrovský T, Steffen KT, Baldrian P. Potential of cometabolic transformation of polysaccharides and lignin in lignocellulose by soil actinobacteria. PLoS ONE. 2014;9:e89108.
Medie FM, Davies GJ, Drancourt M, Henrissat B. Genome analyses highlight the different biological roles of cellulases. Nat Rev Microbiol. 2012;10:227–34.
Wilhelm RC, Cardenas E, Leung H, Szeitz A, Jensen LD, Mohn WW. Long-term enrichment of stress-tolerant cellulolytic soil populations following timber harvesting evidenced by multi-omic stable isotope probing. Front Microbiol. 2017;8:537.
Morris SA, Radajweski S, Willison TW, Murrel JC. Identification of the functionally active methanotroph population in a peat soil microcosm by stable-isotope probing. Appl Environ Microbiol. 2002;68:1446–53.
Berlemont R, Martiny AC. Phylogenetic distribution of potential cellulases in bacteria. Appl Environ Microbiol. 2013;79:1545–54.
Wilhelm RC, Cardenas E, Maas K, Leung H, McNeil L, Berch S, et al. Biogeography and organic matter removal shape long-term effects of timber harvesting on forest soil microbial communities. ISME J. 2017;11:2552–68.
Wilhelm RC, Cardenas E, Hahn AS, Leung H, Maas K, Hartmann M, et al. A multi-omic survey of forest soil microbial communities more than a decade after timber harvesting. Sci Data. 2017;4:170092.
Leung HT, Maas KR, Wilhelm RC, Mohn WW. Long-term effects of timber harvesting on hemicellulolytic microbial populations in coniferous forest soils. ISME J. 2016;10:363–75.
Youngblut N, Buckley DH. Intra-genomic variation in G + C content and its implications for DNA stable isotope probing. Env Microbiol. 2015;6:767–75.
Kirk TK, Brunow G. Synthetic 14C-labeled lignins. In: Methods in enzymology. (eds Colowick SP, Kaplan NO) Vol. 161;65–73 (Academic Press. Cambridge, MA, 1988).
Wilhelm RC, Szeitz S, Klassen TL, Mohn WW. Sensitive, efficient quantitation of 13C-enriched nucleic acids via ultrahigh-performance liquid chromatography–tandem mass spectrometry for applications in stable isotope probing. Appl Environ Microbiol. 2014;80:7206–11.
Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37:911–7.
Churchland C, Grayston SJ, Bengtson P. Spatial variability of soil fungal and bacterial abundance: consequences for carbon turnover along a transition from a forested to clear-cut site. Soil Biol Biochem. 2013;63:5–13.
Neufeld JD, Vohra J, Dumont MG, Leuders T, Manefield M, Friedrich, et al. DNA stable-isotope probing. Nat Protoc. 2007;2:860–6.
Hartmann M, Howes CG, VanInsberghe D, Yu H, Bachar D, Christen R, et al. Significant and persistent impact of timber harvesting on soil microbial communities in Northern coniferous forests. ISME J. 2012;6:2199–218.
Wilhelm RC. Deciphering decomposition and the effects of disturbance in forest soil microbial communities with metagenomics and stable isotope probing. Ph.D. thesis. University of British Columbia; 2016.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucl Acids Res. 2007;35:7188–96.
Westram R, Bader K, Prüsse E, Kumar Y, Meier H, Glöckner FO, et al. ARB: a software environment for sequence data. In: Handbook of molecular microbial ecology I: metagenomics and complementary approaches. 399–406 (John Wiley & Sons, Hoboken, NJ, 2011).
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. 2014;30:2114–20.
Gordon A, Hannon GJ. Fastx-toolkit. Computer program distributed by the author. 2010. http://hannonlab.cshl.edu/fastx_toolkit/index.html. Accessed Nov. 2015.
Boisvert S, Raymond F, Godzaridis É, Laviolette F, Corbeil J. Ray Meta: scalable de novo metagenome assembly and profiling. Genome Biol. 2012;13:R122.
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.
Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165.
Parks D, Imelfort M, Skennerton CT, Hugenholtz P, Tyson G. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.
Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–86.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.
Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 2010;11:119.
Finn RD, Clements J, Eddy S. HMMER web server: interactive sequence similarity searching. Nucl Acids Res. 2011;39:W29–37.
Yin Y, Mao X, Yan J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucl Acids Res. 2012;40:W445–51.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2015.
Wickham H. Reshaping data with the reshape package. J Stat Soft. 2007;21:1–20.
Wickham H. The split-apply-combine strategy for data analysis. J Stat Soft. 2011;40:1–29.
Wickham H. Elegant graphics for data analysis.. 2009. New York: Springer-Verlag; 2009.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2014;8:e61217.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B Met. 1995;Series B:289–300.
Kursa MB, Rudnicki WR. Feature selection with the boruta package. J Stat Softw. 2010;36:1–13.
Grömping U. Relative importance for linear regression in R: the package relaimpo. J Stat Softw. 2006;17:1–27.
Li K, Xu F, Eriksson KE. Comparison of fungal laccases and redox mediators in oxidation of a nonphenolic lignin model compound. Appl Environ Microbiol. 1999;65:2654–60.
Babot ED, Rico A, Rencoret J, Kalum L, Lund H, Romero J, et al. Towards industrially-feasible delignification and pitch removal by treating paper pulp with Myceliophthora thermophila laccase and a phenolic mediator. Biores Technol. 2011;102:6717–22.
Pepe-Ranney C, Campbell AN, Koechli CN, Berthrong S, Buckley DH. Unearthing the ecology of soil microorganisms using a high-resolution DNA-SIP approach to explore cellulose and xylose metabolism in soil. Front Microbiol. 2016;7:703.
Woo HL, Hazen TC, Simmons BA, DeAngelis KM. Enzyme activities of aerobic lignocellulolytic bacteria isolated from wet tropical forest soils. Syst Appl Microbiol. 2014;37:60–7.
Youssef NH, Elshahed MS. Diversity rankings among bacterial lineages in soil. ISME J 2009;3:305–13.
Bergmann GT, Bates ST, Eilers KG, Lauber CL, Caporaso JG, Walters WA, et al. The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol Biochem 2011;43:1450–5.
Sangwan P, Chen X, Hugenholtz P, Janssen PH. Chthoniobacter flavus gen. nov., sp. nov., the First Pure-Culture Representative of Subdivision Two, Spartobacteria classis nov., of the Phylum Verrucomicrobia. Appl Environ Microbiol 2004; 70: 5875–81.
Lladó S, Benada O, Cajthaml T, Baldrian P, Garcia-Fraile P. Silvibacterium bohemicum gen. nov. sp. nov., an acidobacterium isolated from coniferous soil in the Bohemian forest national park. Syst Appl Microbiol. 2015;39:14–9.
Dedysh SN, Kulichevskaya IS. Planctomycetes: Cell Structure, Origins and Biology. 2013;125–39.
Lee KCY, Morgan XC, Dunfield PF, Tamas I, McDonald IR, Stott MB. Genomic analysis of Chthonomonas calidirosea, the first sequenced isolate of the phylum Armatimonadetes. ISME J 2014;8:1522–33.
Eichorst SA, Kuske CR. Identification of cellulose-responsive bacterial and fungal communities in geographically and edaphically different soils by using stable isotope probing. Appl Environ Microbiol. 2012;78:2316–27.
Verastegui Y, Cheng J, Engel K, Kolczynski D, Mortimer S, Lavigne J, et al. Multisubstrate isotope labeling and metagenomic analysis of active soil bacterial communities. mBio. 2014;5:e01157–14.
Folman LB, Klein Gunnewiek PJA, Boddy L, de Boer W. Impact of white-rot fungi on numbers and community composition of bacteria colonizing beech wood from forest soil. FEMS Microbiol Ecol. 2008;63:181–91.
Valáková V, De Boer W, Klein Gunnewiek PJA, Pospíek M, Baldrian P. Phylogenetic composition and properties of bacteria coexisting with the fungus Hypholoma fasciculare in decaying wood. ISME J 2009;3:1218–21.
Smenderovac E. Microbial Community Structure and Function in Coarse Woody Debris and Boreal Forest Soil after Intensified Biomass Harvests. University of Toronto, PhD diss. 2014;140.
Henrici AT, Johnson DE. Studies of Freshwater Bacteria. J Bacteriol 1935;30:61–93.
Poindexter JS. The Caulobacters: ubiquitous unusual bacteria. Microbiol Rev. 1981;45:123–79.
Nierman WC, Feldblyum TV, Laub MT, Paulsen IT, Nelson KE, Eisen J, et al. Complete genome sequence of Caulobacter crescentus. PNAS. 2001;98:4136–41.
Hottes AK, Meewan M, Yang D, Arana N, Romero P, McAdams HH, et al. Transcriptional profiling of Caulobacter crescentus during growth on complex and minimal media. J Bacteriol 2004;186:1448–61.
Thanbichler M, Iniesta AA, Shapiro L. A comprehensive set of plasmids for vanillate - And xylose-inducible gene expression in Caulobacter crescentus. Nucleic Acids Res 2007;35.
Song Z, Vail A, Sadowsky MJ, Schilling JS. Influence of hyphal inoculum potential on the competitive success of fungi colonizing wood. Microb Ecol. 2015;69:758–67.
Presley GN, Payea MJ, Hurst LR, Egan AE, Martin BS, Periyannan GR. Extracellular gluco-oligosaccharide degradation by Caulobacter crescentus. Microbiol (United Kingdom) 2014;160:635–45.
Wang X, Sharp CE, Jones GM, Grasby SE, Brady AL, Dunfield PF. Stable-isotope probing identifies uncultured Planctomycetes as primary degraders of a complex heteropolysaccharide in soil. Appl Environ Microbiol 2015;81:4607–15.
Ausec L, van Elsas JD, Mandic-Mulec I. Two- and three-domain bacterial laccase-like genes are present in drained peat soils. Soil Biol Biochem 2011;43:975–83.
Colpa DI, Fraaije MW, Van Bloois E. DyP-type peroxidases: A promising and versatile class of enzymes. J Ind Microbiol Biotechnol 2014;41:1–7.
Singh R, Eltis LD. The multihued palette of dye-decolorizing peroxidases. Arch Biochem Biophys 2015;574:56–65.
This study was supported by a Large Scale Applied Research Project (162MIC) funded by Genome Canada and Genome BC. RCW was supported by a NSERC graduate scholarship and UBC Four Year Fellowship.
Correspondence to William W. Mohn.

References: v. 
 v. 
 v. 
 v. 
 v. 
 v. 
 v. 
 v.