Biogeography and organic matter removal shape long-term effects of timber harvesting on forest soil microbial communities

Towards repeated clear-cutting of boreal forests - a tipping point for biodiversity?

Boreal forests are important carbon sinks and host a diverse array of species that provide important ecosystem functions. Boreal forests have a long history of intensive forestry, in which even-aged management with clear-cutting has been the dominant harvesting practice for the past 50-80 years. As a second cycle of clear-cutting is emerging, there is an urgent need to examine the effects of repeated clear-cutting events on biodiversity. Clear-cutting has led to reduced numbers of old and large trees, decreased volumes of dead wood of varied decay stages and diameters, and altered physical and chemical compositions of soils. The old-growth boreal forest has been fragmented and considerably reduced. Here, we review short- and long-term (≥50 years) effects of clear-cutting on boreal forest biodiversity in four key substrates: living trees, dead wood, ground and soil. We then assess landscape-level changes (habitat fragmentation and edge effects) on this biodiversity. There is evidence for long-term community changes after clear-cutting for several taxa: epiphytic lichens; saproxylic fungi, bryophytes and insects; epigeic bryophytes; and soil snails, bacteria, and ectomycorrhizal fungi. Long-term declines in species richness were found for saproxylic fungi, bryophytes and true flies. However, for the majority of taxa, long-term effects of clear-cutting are not well understood. On the landscape level, reduced connectivity to old-growth forests has negative effects on several species of fungi, lichens, bryophytes and insects, notably among Red-Listed species. Furthermore, altered microclimate near clear-cut edges negatively affects epiphytic lichens and epigeic arthropods, implying complex effects of habitat fragmentation. Repeated cycles of clear-cutting might pose even stronger pressures on boreal forest biodiversity due to continued fragmentation of old-growth forests and accumulation of extinction debts. Examining the broad effects of forestry on biodiversity across the boreal biome is crucial: (i) to increase our knowledge of long-term and landscape-level effects of former clear-cutting; and (ii) to gain a better understanding of how forestry will affect biodiversity and, subsequently, ecosystem functioning, with repeated cycles of clear-cutting.

Preliminary construction of a microecological evaluation model for uranium-contaminated soil

With the extensive development of nuclear energy, soil uranium contamination has become an increasingly prominent problem. The development of evaluation systems for various uranium contamination levels and soil microhabitats is critical. In this study, the effects of uranium contamination on the carbon source metabolic capacity and microbial community structure of soil microbial communities were investigated using Biolog microplate technology and high-throughput sequencing, and the responses of soil biochemical properties to uranium were also analyzed. Then, ten key biological indicators as reliable input variables, including arylsulfatase, biomass nitrogen, metabolic entropy, microbial entropy, Simpson, Shannon, McIntosh, Nocardioides, Lysobacter, and Mycoleptodisus, were screened by random forest (RF), Boruta, and grey relational analysis (GRA). The optimal uranium-contaminated soil microbiological evaluation model was obtained by comparing the performance of three evaluation methods: partial least squares regression (PLS), support vector regression (SVR), and improved particle algorithm (IPSO-SVR). Consequently, partial least squares regression (PLS) has a higher R2 (0.932) and a lower RMSE value (0.214) compared to the other. This research provides a new evaluation method to describe the relationship between soil ecological effects and biological indicators under nuclear contamination.

Spatial distribution patterns across multiple microbial taxonomic groups

Even in the vertical dimension, soil bacterial communities are spatially distributed in a distance-decay relationship (DDR). However, whether this pattern is universal among all soil microbial taxonomic groups, and how body size influences this distribution, remains elusive. Our study consisted of obtaining 140 soil samples from two adjacent ecosystems in the Yellow River Delta (YRD), both nontidal and tidal, and measuring the DDR between topsoil and subsoil for bacteria, archaea, fungi and protists (rhizaria). Our results showed that the entire community generally fitted the DDR patterns (P < 0.001), this was also true at the kingdom level (P < 0.001, with the exception of the fungal community), and for most individual phyla (47/75) in both ecosystems and with soil depth. Meanwhile, these results presented a general trend that the community turnover rate of nontidal soils was higher than tidal soils (P < 0.05), and that the rate of topsoil was also higher than that of subsoil (P < 0.05). Additionally, microbial spatial turnover rates displayed a negative relationship with body sizes in nontidal topsoil (R² = 0.29, P = 0.009), suggesting that the smaller the body size of microorganisms, the stronger the spatial limitation was in this environment. However, in tidal soils, the body size effect was negligible, probably owing to the water's fluidity. Moreover, community assembly was judged to be deterministic, and heterogeneous selection played a dominant role in the different environments. Specifically, the spatial distance was much more influential, while the soil salinity in these ecosystems was the major environmental factor in selecting the distributions of microbial communities. Overall, this study revealed that microbial community compositions at different taxonomic levels followed relatively consistent distribution patterns and mechanisms in this coastal area.

Microbial gradual shifts during the process of species replacement in Taihang Mountain

Introduction

Understanding microbial gradual shifts along species replacement can help elucidate the mechanisms driving secondary succession, and predict microbial responses to changing environments. However, how climate-induced species replacement alters microbial processes, and whether microbial shifts follow predictable assembly trajectories remain unclear.

Methods

Using space-for-time substitution approach, we studied shifts in bacterial and fungal communities in the succession from Leptodermis oblonga to Vitex negundo var. heterophylla shrubland in Taihang Mountain.

Results and Discussion

Species replacement, induced by climate related environmental change, significantly increased the above-ground biomass of shrublands, and TP and TK contents in topsoil. The succession from L. oblonga to V. negundo var. heterophylla communities resulted in the gradually replacement of cold-tolerant microbes with warm-affinity ones, and alterations of microbial communities involved in soil biogeochemical processes. Soil and plant variables, such as above-ground biomass, soil pH, total phosphorus, and total potassium, well explained the variations in microbial communities, indicating that the coordinated changes in plant communities and soil properties during secondary succession caused accompanied shifts in microbial diversity and composition.

Ecological insights into soil health according to the genomic traits and environment-wide associations of bacteria in agricultural soils

Soil microbiomes are sensitive to current and previous soil conditions, and bacterial 'bioindicators' of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

Microbial community shifts correspond with suppression of decomposition 25 years after liming of acidic forest soils

Microbial community structure and function regularly covary with soil pH, yet effects of these interactions on soil carbon are rarely tested experimentally within natural ecosystems. We investigated the enduring (25 year) impacts of liming on microbial community structure and decomposition at an acidic northern hardwood forest, where experimental liming increased pH one unit and surprisingly doubled the organic carbon stocks of the forest floor. We show that this increase in carbon storage corresponded with restructuring of the bacterial and fungal communities that drive decomposition. In the Oe horizon, liming reduced the activities of five extracellular enzymes that mediate decomposition, while the Oa horizon showed an especially large (64%) reduction in the activity of a sixth, peroxidase, which is an oxidative enzyme central to lignocellulose degradation. Decreased enzyme activities corresponded with loss of microbial taxa important for lignocellulose decay, including large reductions in the dominant ectomycorrhizal genera Russula and Cenococcum, saprotrophic and wood decaying fungi, and Actinobacteria (Thermomonosporaceae). These results demonstrate the importance of pH as a dominant regulator of microbial community structure and illustrate how changes to this structure can produce large, otherwise unexpected increases in carbon storage in forest soils.

Tree Cover Species Modify the Diversity of Rhizosphere-Associated Microorganisms in Nothofagus obliqua (Mirb.) Oerst Temperate Forests in South-Central Chile

Chilean native forests have been subjected to several types of disturbances, with one of them being the replacement by exotic species. Pinus radiata D. Don is a widespread exotic tree that forms extensive plantations in southern Chile. It covers extended areas, affecting the landscape, biodiversity, and ecosystem services associated with native forest ecosystems. Although advances in assessing the impact of exotic plant species have been conducted, few studies have focused on the alteration of soil microorganisms. This study aimed to characterize the rhizosphere bacterial and fungal communities associated with the tree species Nothofagus obliqua inside a native forest stand and within a P. radiata plantation growing nearby. We used a 16S rRNA gene and ITS region metabarcoding approach. Using bioinformatics, diversity indices, relative abundance, preferential taxa, and predicted functions and guilds were estimated. The β-diversity analysis showed that both factors, the type of soil (rhizosphere or bulk soil) and the type of site (native forest or P. radiata plantation), were significant, with the site explaining most of the variation among bacterial and fungal communities. Proteobacteria and Basidiomycota were the most abundant bacterial and fungal phyla in both types of soil and sites. Similarly, bacteria showed similar abundant taxa at the family level, independent of the soil type or the site. The main fungal taxa associated with native forests were Tricholomataceae and Cantharellales, whereas in P. radiata plantations, Russulaceae and Hyaloscyphaceae were the most abundant families. The main bacteria functional groups were chemoheterotrophy and aerobic chemoheterotrophy, without significant differences between the type of soil or sites. Overall, these results demonstrate that the composition and diversity of bacterial and fungal communities associated with native N. obliqua forest are influenced by the surrounding forest, and mainly depend on the site’s characteristics, such as the lignin-rich wood source. These results improve our understanding of the impact of native forest replacement on soil microbial communities, which can alter microbial-related soil ecosystem services.

Towards the biogeography of prokaryotic genes

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats^1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analysed these traits from publicly available metagenomes derived from marine, soil, host-associated and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size—reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments—a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

Disturbance alters the forest soil microbiome

Billions of microorganisms perform critical below-ground functions in all terrestrial ecosystems. While largely invisible to the naked eye, they support all higher lifeforms, form symbiotic relationships with ~90% of terrestrial plant species, stabilize soils, and facilitate biogeochemical cycles. Global increases in the frequency of disturbances are driving major changes in the structure and function of forests. However, despite their functional significance, the disturbance responses of forest microbial communities are poorly understood. Here, we explore the influence of disturbance on the soil microbiome (archaea, fungi and bacteria) of some of the world's tallest and most carbon-dense forests, the Mountain Ash forests of south-eastern Australia. From 80 sites, we identified 23,277 and 19,056 microbial operational taxonomic units from the 0-10 cm and 20-30 cm depths of soil respectively. From this extensive data set, we found the diversity and composition of these often cryptic communities has been altered by human and natural disturbance events. For instance, the diversity of ectomycorrhizal fungi declined with clearcut logging, the diversity of archaea declined with salvage logging, and bacterial diversity and overall microbial diversity declined with the number of fires. Moreover, we identified key associations between edaphic (soil properties), environmental (slope, elevation) and spatial variables and the composition of all microbial communities. Specifically, we found that soil pH, manganese, magnesium, phosphorus, iron and nitrate were associated with the composition of all microbial communities. In a period of widespread degradation of global forest ecosystems, our findings provide an important and timely insight into the disturbance responses of soil microbial communities, which may influence key ecological functions.

Dissolved organic matter (DOM) quality drives biogeographic patterns of soil bacterial communities and their association networks in semi-arid regions

It is of great interest to elucidate the biogeographic patterns of soil microorganisms and their driving forces, which is fundamental to predicting alterations in microbial-mediated functions arising from environment changes. Although dissolved organic matter (DOM) represents an important resource for soil microorganisms, knowledge of how its quality affects microbial biogeography is limited. Here, we characterized soil bacterial communities and DOM quality in 45 soil samples collected from a 1500-km sampling transect through semi-arid regions in northern China which are currently suffering great pressure from climate change, using Illumina Miseq sequencing and fluorescence spectroscopy, respectively. We found that DOM quality (i.e. the source of DOM and the humification degree of DOM) had profound shaping influence on the biogeographic patterns exhibited by bacterial diversity, community composition and association networks. Specifically, the composition of bacteria community closely associated with DOM quality. Plant-derived DOM sustained higher bacterial diversity relative to microbial-derived DOM. Meanwhile, bacterial diversity linearly increased with increasing humification degree of DOM. Additionally, plant-derived DOM was observed to foster more complex bacterial association networks with less competition. Together, our work contributes to the factors underlying biogeographic patterns not only of bacterial diversity, community composition but also of their association networks and reports previously undocumented important role of DOM quality in shaping these patterns.

Changes in soil fungal community composition depend on functional group and forest disturbance type

Disturbances have altered community dynamics in boreal forests with unknown consequences for belowground ecological processes. Soil fungi are particularly sensitive to such disturbances; however, the individual response of fungal guilds to different disturbance types is poorly understood. Here, we profiled soil fungal communities in lodgepole pine forests following a bark beetle outbreak, wildfire, clear-cut logging, and salvage-logging. Using Illumina MiSeq to sequence ITS1 and SSU rDNA, we characterized communities of ectomycorrhizal, arbuscular mycorrhizal, saprotrophic, and pathogenic fungi in sites representing each disturbance type paired with intact forests. We also quantified soil fungal biomass by measuring ergosterol. Abiotic disturbances changed the community composition of ectomycorrhizal fungi and shifted the dominance from ectomycorrhizal to saprotrophic fungi compared to intact forests. The disruption of the soil organic layer with disturbances correlated with the decline of ectomycorrhizal and the increase of arbuscular mycorrhizal fungi. Wildfire changed the community composition of pathogenic fungi but did not affect their proportion and diversity. Fungal biomass declined with disturbances that disrupted the forest floor. Our results suggest that the disruption of the forest floor with disturbances, and the changes in C and nutrient dynamics it may promote, structure the fungal community with implications for fungal biomass-C.

Distinct fungal successional trajectories following wildfire between soil horizons in a cold-temperate forest

Soil fungi represent a major component of below-ground biodiversity that determines the succession and recovery of forests after disturbance. However, their successional trajectories and driving mechanisms following wildfire remain unclear. We examined fungal biomass, richness, composition and enzymes across three soil horizons (Oe, A1 and A2) along a near-complete fire chronosequence (1, 2, 8, 14, 30, 49 and c. 260 yr) in cold-temperate forests of the Great Khingan Mountains, China. The importance of soil properties, spatial distance and tree composition were also tested. Ectomycorrhizal fungal richness and β-glucosidase activity were strongly reduced by burning and significantly increased with 'time since fire' in the Oe horizon but not in the mineral horizons. Time since fire and soil C : N ratio were the primary drivers of fungal composition in the Oe and A1/A2 horizons, respectively. Ectomycorrhizal fungal composition was remarkably sensitive to fire history in the Oe horizon, while saprotroph community was strongly affected by time since fire in the deeper soil horizon and this effect emerged 18 years after fire in the A2 horizon. Our study demonstrates pronounced horizon-dependent successional trajectories following wildfire and indicates interactive effects of time since fire, soil stoichiometry and spatial distance in the reassembly of below-ground fungal communities in a cold and fire-prone region.

Long-term harvesting of reeds affects greenhouse gas emissions and microbial functional genes in alkaline wetlands

Reed (Phragmites australis) is dominant vegetation in alkaline wetlands that is harvested annually due to its economic value. To reveal the effects of harvesting reeds on the emission of greenhouse gases (GHG), the annual soil physicochemical characteristics and flux of GHGs in a reed wetland without harvesting (NHRW) and with harvesting (HRW) were measured. The results showed that after the harvesting of reeds, the total organic carbon (TOC) and total nitrogen (TN) significantly decreased, and soil temperature significantly increased. The annual cumulative N₂O emissions decreased from 0.73 ± 0.20 kg ha^-1 to -0.57 ± 0.49 kg ha^-1 with the harvesting of reeds. The annual cumulative CH₄ emissions also decreased from 561.88 ± 18.61 kg ha^-1 to 183.13 ± 18.77 kg ha^-1 with the harvesting of reeds. However, harvesting of reeds had only a limited influence on the annual cumulative CO₂ emissions. A Pearson correlation analysis revealed that the CO₂ and N₂O emissions were more sensitive to temperature than the CH₄ emissions. Both structural equation modeling (SEM) analysis and slurry incubation confirmed that higher temperatures offset the reduction of CO₂ emissions after reed harvesting. Metagenomics showed that the abundance of functional genes involved in both GHG sink and source decreased with reed harvesting. This study presents a comprehensive view of reed harvesting on GHG emissions in alkaline wetlands, yielding new insight into the microbial response and offering a novel perspective on the potential impacts of wetland management.

Metaphenomic Responses of a Native Prairie Soil Microbiome to Moisture Perturbations

Climate change is causing shifts in precipitation patterns in the central grasslands of the United States, with largely unknown consequences on the collective physiological responses of the soil microbial community, i.e., the metaphenome. Here, we used an untargeted omics approach to determine the soil microbial community's metaphenomic response to soil moisture and to define specific metabolic signatures of the response. Specifically, we aimed to develop the technical approaches and metabolic mapping framework necessary for future systematic ecological studies. We collected soil from three locations at the Konza Long-Term Ecological Research (LTER) field station in Kansas, and the soils were incubated for 15 days under dry or wet conditions and compared to field-moist controls. The microbiome response to wetting or drying was determined by 16S rRNA amplicon sequencing, metatranscriptomics, and metabolomics, and the resulting shifts in taxa, gene expression, and metabolites were assessed. Soil drying resulted in significant shifts in both the composition and function of the soil microbiome. In contrast, there were few changes following wetting. The combined metabolic and metatranscriptomic data were used to generate reaction networks to determine the metaphenomic response to soil moisture transitions. Site location was a strong determinant of the response of the soil microbiome to moisture perturbations. However, some specific metabolic pathways changed consistently across sites, including an increase in pathways and metabolites for production of sugars and other osmolytes as a response to drying. Using this approach, we demonstrate that despite the high complexity of the soil habitat, it is possible to generate insight into the effect of environmental change on the soil microbiome and its physiology and functions, thus laying the groundwork for future, targeted studies.IMPORTANCE Climate change is predicted to result in increased drought extent and intensity in the highly productive, former tallgrass prairie region of the continental United States. These soils store large reserves of carbon. The decrease in soil moisture due to drought has largely unknown consequences on soil carbon cycling and other key biogeochemical cycles carried out by soil microbiomes. In this study, we found that soil drying had a significant impact on the structure and function of soil microbial communities, including shifts in expression of specific metabolic pathways, such as those leading toward production of osmoprotectant compounds. This study demonstrates the application of an untargeted multi-omics approach to decipher details of the soil microbial community's metaphenotypic response to environmental perturbations and should be applicable to studies of other complex microbial systems as well.

Bacterial contributions to delignification and lignocellulose degradation in forest soils with metagenomic and quantitative stable isotope probing

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 ¹³C-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 ¹³C-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.

Following the terrestrial tracks of Caulobacter - redefining the ecology of a reputed aquatic oligotroph

For the past 60 years Caulobacter spp. have been commonly attributed an aquatic and oligotrophic lifestyle yet are not uncommon in nutrient-rich or soil environments. This study evaluates the environmental and ecological associations of Caulobacter to reconcile past evidence, largely limited to culturing and microscopy, with currently available metagenomic and genomic data. The distribution of Caulobacter species and their characteristic adhesion-conferring genes, holdfast (hfaAB), were determined using collections of 10,641 16S rRNA gene libraries (196 studies) and 2625 shotgun metagenomes (190 studies) from a range of terrestrial and aquatic environments. Evidence for ecotypic variation was tested in 26 genomes sourced from soil, rhizosphere, plant, groundwater, and water. Caulobacter were, on average, fourfold more relatively abundant in soil than in aquatic environments, and abundant in decomposing wood, compost, and particulate matter (in air and water). Caulobacter holdfast genes were 35-fold more abundant in soils than aquatic environments. Ecotypic differences between soil and aquatic Caulobacter were evident in the environmental associations of several species and differences in genome size and content among isolates. However, most abundant species were common to both environments, suggesting populations exist in a continuum that was evident in the re-analysis of studies on the temporal dynamics of, and sources of bacterioplankton to, lakes and rivers. This study provides a new perspective on the ecological profile of Caulobacter, demonstrating that members of this genus are predominantly soil-borne, possess an overlooked role in plant matter decomposition and a dependency on water-mediated dispersal.

A metagenomic survey of forest soil microbial communities more than a decade after timber harvesting

The scarcity of long-term data on soil microbial communities in the decades following timber harvesting limits current understanding of the ecological problems associated with maintaining the productivity of managed forests. The high complexity of soil communities and the heterogeneity of forest and soil necessitates a comprehensive approach to understand the role of microbial processes in managed forest ecosystems. Here, we describe a curated collection of well replicated, multi-faceted data from eighteen reforested sites in six different North American ecozones within the Long-term Soil Productivity (LTSP) Study, without detailed analysis of results or discussion. The experiments were designed to contrast microbial community composition and function among forest soils from harvested treatment plots with varying intensities of organic matter removal. The collection includes 724 bacterial (16S) and 658 fungal (ITS2) amplicon libraries, 133 shotgun metagenomic libraries as well as stable isotope probing amplicon libraries capturing the effects of harvesting on hemicellulolytic and cellulolytic populations. This collection serves as a foundation for the LTSP Study and other studies of the ecology of forest soil and forest disturbance.