Introduction

Nutrient mobilization is essential to life on earth (Read and Perez‐Moreno 2003; Treseder and Lennon 2015, Bardgett 2016). Decomposition recycles organic material to structure soils, sequester carbon, and shape trophic interactions (Cebrian 1999; Gessner et al. 2010; Matulich and Martiny 2015). Numerous microbial functional groups possess the enzymatic machinery to contribute to decomposition (Swift et al. 1979; Wall and Moore 1999; Waldrop and Firestone 2004; d’Annunzio et al. 2008; Wei et al. 2009; Shah et al. 2016), but the relationship between microbial composition and nutrient cycling in soils likely depends on the quality, identity, quantity, and frequency of deposited organic materials (Li et al. 2019; Glassman et al. 2018; Chapman et al. 2013; Zeng et al. 2017). In woodlands, dead plant material first encounter soil communities as litter. Litter introduction, including leaves, roots, stems or woody material, can promote specific microbial taxa, perhaps as a result of differences in polyphenol to nitrogen ratio (Diedhiou et al. 2009; Baldwin et al. 1983; Kuiters 1990; Hättenschwiler and Vitousek 2000). Numerous studies show that leaf litter input can influence the soil microbial profile, including microbial relative abundance, activity, and function (Diedhiou et al. 2009; Aneja et al. 2006; Mukhopadhyay and Joy 2010). Despite the large number of studies that have examined the relationship between leaf litter and soil microorganisms, little is known about bark litter in this context.

Bark is diverse in chemical composition, and its stem volume varies across species (Harkin and Rowe 1971; Gregoire et al. 1993). The outer bark is external to the vascular cambium and sieves elements that help translocate nutrients from leaves to roots (Rosell 2016; Notaguchi and Okamoto 2015). The targeted role of bark is to protect stems from fire, desiccation, herbivores, and pathogens (Dantas and Pausas 2013; Cernusak and Cheesman 2015; Zas et al. 2011; Pearce 1996; Mullick 1977). Because of a primary role in protection rather than absorption, bark, in contrast to leaves and fine roots, is highly recalcitrant and slow to degrade. However, several forces, including freeze–thaw events and disturbance from animals, can fragment bark. Bark fragmentation and displacement creates a litter source that is low in nutrient content (Dossa et al. 2016, 2018; Harmon 1989; Shorohova and Kapitsa 2014), but rich in antimicrobial tannins that can impact microbial profiles (Field and Lettinga 1992; Updegraff and Grant 1975). For example, bark contains antimicrobial tannins of a relatively high molecular weight, compared to wood (Feng et al. 2013). As an anti-microbial compound, tannins were reported to lower mycotoxin production and suppress nitrogen mineralization (Peng et al. 2018; Zhang and Laanbroek 2018).

Determining the relationship between bark litter and soil microorganisms will provide further insight into nutrient cycling and brown food-web processes. Soil communities that are regularly exposed to litter deposited from a nearby plant can be thought of as part of the ‘home soil’ for that plant, which is central to the home-field advantage (HFA) hypothesis (Diepen et al. 2017; Strickland et al. 2009; Ayres et al. 2009b). According to HFA, soil microorganisms that are most efficient at metabolizing litter occur in the ‘home soil’, as a result of a specialized decomposer-litter relationship (Palozzi and Lindo 2018). The HFA hypothesis has been supported in numerous studies that have examined leaf litter decomposition (Ayres et al. 2009b; Veen et al. 2015; Wang et al. 2013; Perez et al. 2013), but to our knowledge this is the first study to focus on bark in this context. As it pertains to bark decomposition, Shorohova et al. (2012), estimated bark decomposition by observing spruce, pine, and birch stumps; and Wu et al. (2019) examined HFA for intact wood (i.e. branches and bark covered wood).

To date, many studies have expanded our understanding of ecological factors governing HFA, including edaphic environment, landscape position, global change, and litter traits (Hobbie et al. 2006; Veen et al. 2018; Diepen et al. 2017; Freschet et al. 2012; Fanin et al. 2016). However, most studies that attempt to test the HFA hypothesis hinge on labile litter models, and do not consider recalcitrant litter such as bark. To advance our understanding of this process, this study tests the HFA hypothesis for recalcitrant bark litter. While the definition of home-field advantage (HFA) does not include scale, previous studies suggest that HFA mainly occurs between ecosystems (Ayres et al. 2009a; Wang et al. 2013; Cizungu et al. 2014; Chomel et al. 2015; Sun and Zhao 2016; Asplund et al. 2018). Here we performed a study that tests HFA within an ecosystem, beneath the canopy of two temperate tree species. We address two specific questions: (1) Is the ‘home’ environment a key factor in determining rates of bark decomposition? (2) Do decomposition of different bark types impact microbial assemblages in the adherent soil? To address these questions, bark from a softwood tree species, eastern hemlock, and a hardwood tree species, white oak, were reciprocally transplanted into ‘home’ and ‘away’ environments, in a mixed hardwood forest in central Pennsylvania. We hypothesized that (1) the ‘home’ environment would increase decomposition rates, and (2) that decomposing bark can initiate change in soil microbial composition.

Methods

Study system

Mature white oak and eastern hemlock are widespread at the Susquehanna - Shale Hills Critical Zone Observatory (40° 39′ N, 77° 54′ W) in central Pennsylvania. The average temperature in this region is 10 °C and the mean annual precipitation is 1,006 mm. The underlying soil is derived from shale parent rock residuum (Hasenmueller et al. 2017). The common soil series at this site is Berks–Weikert (Order: Inceptisol, Great Group: Dystrudepts). The percentage of channery shale has been reported to increase with soil depth (White et al. 2015), and the pH in the top 10 cm of soil is ~ 4.0 (Malik 2019). Bark from eastern hemlock (Tsuga canadensis) and white oak (Quercus alba) were used for this study, as these two species have contrasting life histories and traits. Eastern hemlock is a softwood, late successional, gymnosperm; while white oak is a hardwood, mid-successional, angiosperm. In addition, these two species face different environmental challenges. Eastern hemlock is on the decline from an invasive sap-sucking Hemipteran, known as hemlock wooly adelgid (Mahan et al. 2004), while white oak is on the decline as a result of fire suppression (Abrams 2003). Importantly, these two species deposit distinct types of litter. Eastern hemlock deposits chemically recalcitrant needles, whereas white oak deposits more labile broad leaves. We would expect these differences in litter input to drive differentiation of soil microbial communities, making this an effective system for contrasting ‘home’ environment effects on bark decomposition.

Experimental design

This experiment tested for the interactive effect of bark type and soil type on bark decomposition and soil microbial composition. We used a full factorial design featuring two tree species × two soil types (‘home’ and ‘away’) with 25 replicates for each, for a total of 100 observations. Since bark volume is related to stem diameter (Paine et al. 2010), mid-late successional trees were used, so that stem diameter ratios were relatively consistent across individual trees. The ‘home’ soil for each bark sample was designated as the soil below the canopy of the tree from which it was collected, while the ‘away’ soil was designated as the soil below the canopy of a randomly selected tree of the other species (Figs. 1b, S1).

Fig. 1
figure 1

Bark introduction may impact soil microbes. To test this hypothesis, bark was introduced into ‘home’ and ‘away’ soils. It is expected that bark decomposition and solubilization can influence microbial assemblages in ‘home’ and ‘away’ soils and vice versa. Soil regimes can be influenced by land use history, legacy effects, and dominant plants (a). Despite this, bark introduction (b) may or may not have an effect on microbial assemblages (c, d). Note illustration of white oak and eastern hemlock are provided by World Almanac for Kids.com, an imprint of Infobase. Reprinted by permission of the Publisher

Bark excision

Bark was excised from 50 individual trees, with an excision no deeper between 1.0 and 1.5 cm. The depth of excision, which was consistent across tree individuals, included bark tissue from the dead cell exterior to the outer cambium. The excised bark was then split in half, with one portion designated for burial in the soil community beneath its parent tree (‘home’ environment), and the other designated for burial in the ‘away’ soil environment (Fig. 1). After the split, each portion was roughly 5 cm × 2.5 cm, while the bark depth was between 1.0 and 1.5 cm. The consistency of bark dimensions was important since decomposability can relate to surface area (Dossa et al. 2018). A small bark subsample (~ 1 g) was frozen at – 20 °C for subsequent analysis of microbial composition. Prior to burial, bark was returned to the lab for initial assessment of bark mass. Bark was oven dried for 20 h at 40 °C since water content can vary across tree species (Rosell et al. 2015), and presumably individual samples.

Field burial and incubation

Dried bark samples were enclosed in cylindrical mesh tubes (10 cm long × 5 cm diameter “cores”) with 0.5 cm × 0.5 cm window openings. Two cores were physically paired at each burial site, with one containing bark from the adjacent tree (‘home’ treatment) and one containing bark from the other tree species (Fig. 1). Paired cores were filled with soil collected at the burial site, and buried horizontally ~ 37.5 cm from the trunk of the ‘home’ tree, similar to Malik (2019). Excavations for core burial were made to 10 cm, since decomposing woody debris frequently occurs at shallow depths (Posada et al. 2012). In addition to bark and decomposing wood being found at this depth, placing bark just few cm below ground enabled bark to stay at a constant depth across samples, as soil surface placement would have led to stochastic uneven burials. Field incubation occurred from July 2017 to June 2018, after which cores were removed from the field and brought to the laboratory for analysis.

Bark samples were removed from each core and loose soil was shaken off. Directly adherent soil was collected and transferred into 1.5-ml Eppendorf tubes and stored at − 20 °C before processing for analysis of microbial composition. Bark was then dried for 96 h at 90 °C. Bark dry mass was then weighed to assess decomposition based on proportionate mass loss.

High-throughput sequencing of bacterial and fungal composition

After a year of decomposition, microbial DNA was assessed in both soil and bark. Initial bark samples were placed in 15-ml sterile water, and agitated overnight at 300 rpm. The overnight agitation fragmented the bark which increased the surface area as well as access to internal microbes. The bark fragments were then centrifuged for 30 min at 21,000 rpm. The supernatant was removed and the pellet was used for DNA extraction. Bark fragments were pulverized with ceramic beads during the DNA extraction protocol, thus, enabling microbes living in the bark interior to also be lysed, sampled and extracted. This approach is common to root microbiome sampling (Lundberg et al. 2012). DNA extraction, which included both the bark pellet and the adherent soil colloids, underwent DNA extraction following the manufacturer’s recommendations (NucleoSpin 96 Soil, Macherey–Nagel, Bethlehem, PA, USA).

To examine the bacterial and fungal composition of the bark and adherent soil, a two-step amplification approach was used to produce Illumina-compatible amplicons. The bacterial primers 515F (5′-GTGYCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACNVGGGTWTCTAAT-3′) were used to target the v4 region of the 16S rRNA gene (Apprill et al. 2015; Parada et al. 2016), and the fungal primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and 58A2R (5′-CTGCGTTCTTCATCGAT-3′) were used to target the ITS1 region of the ribosomal RNA gene (Gardes and Bruns 1993; Martin and Rygiewicz 2005). Primers were designed with overhangs that allow the attachment of barcodes and standard Illumina overhang adaptors in a second PCR step.

Initial PCR was performed in mixtures of: 8 μl of 5Prime HotStart MasterMix (Quanta BioSciences Inc., Beverly, MA, USA), 0.2–4 μl of template DNA, 1 μl of each primer from 10 μM stocks, and the appropriate volume of molecular biology grade water to bring the reaction to 20 μl. The following protocol was performed for 16S rRNA gene amplicons: 94 °C for 3 min; 25 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 45 s; and then a final elongation for 10 min at 72 °C. Separately, ITS amplicons were produced using the following protocol: 94 °C for 3 min; 35 cycles of 94 °C for 20 s, 45 °C for 30 s, and 72 °C for 45 s; and then a final elongation for 5 min at 72 °C. Post-PCR, amplicons were purified using Mag-Bind TotalPure NGS (Omega Bio-Tek, Norcross, GA, USA) magnetic beads. Next, standard Illumina overhang adaptors and unique indexing barcodes were added to both 16S and ITS amplicons from each sample during a second PCR step where 5 μl of cleaned PCR product, 12.5 μl of 5Prime Hotstart Mastermix, 2.5 μl of water, and 2.5 μl of index primer (10 μM) were added to the reactions and amplified under the following conditions: 98 °C for 1 min; 8 cycles of 98 °C for 15 s, 55 °C for 30 s, and 72 °C for 20 s, followed by a final elongation for 5 min at 72 °C. All amplifications were made using Mastercycler Nexus Gradient PCR machines (Eppendorf, Hamburg, Germany). Amplicon yields were normalized using the SequalPrep Normalization Plate Kit (Invitrogen, Carlsbad, CA, USA) and pooled. The pool was concentrated using a Savant SpeedVac (Thermo Scientific, Waltham, MA, USA) for 3 h at 50 °C and subsequently run on an agarose gel (1.2%). The band of expected size was then excised and purified using the PureLink Quick Gel Extraction kit (Invitrogen, Carlsbad, CA, USA).

Sequencing was performed at the Cornell University Biotechnology Resource Center Genomics Facility on an Illumina MiSeq (2 × 250 cycle, v2 kit). Bacterial 16S rRNA gene sequencing yielded a total of 4,541,319 contigs following merging of forward and reverse reads across 90 libraries. After quality filtering, 631,548 total reads were obtained across 90 samples. This provided an average of 7017 reads/sample. Fungal ITS sequencing yielded a total of 5,386,484 contigs after merging forward and reverse reads across 90 libraries. After quality filtering, we obtained 2,302,039 total reads across 90 samples, providing an average of 25,578 reads/sample. Raw sequence reads were submitted the NCBI SRA under the accession PRJNA587721.

Bioinformatic and data assessment

Read processing was performed mostly as described by Howard et al. (2017), based on a modified Brazilian Microbiome Project pipeline for both the 16S and ITS datasets (Pylro et al. 2014). Briefly, forward and reverse reads were merged (make.contigs) and primer sequences trimmed (trim.seqs, pdiffs = 2, maxambig = 0) in Mothur v1.39.5 (Schloss et al. 2009) as singletons were removed. Reads were then split by sample into individual files for downstream processing in MacQIIME v1.9.1 (Caporaso et al. 2010). Qiime-compatible labels were added to the sequences which were then dereplicated (-derep_fulllength) with VSEARCH v2.3.4 (Rognes et al. 2016), 97% OTUs picked (-cluster_otus) with USEARCH v7 (Edgar 2010). Chimeric reads were removed (-uchime_ref) in USEARCHv7 using the rdp_gold database. OTUs were then assigned taxonomy (classify.seqs) in Mothur v1.39.5 using the May 2013 greengenes release for 16S rRNA gene OTUs and the UNITEv7 database for fungal ITS OTUs. For the 16S rRNA gene dataset, OTUs classified as ‘Archaea’, ‘Eukaryota’, ‘chloroplast’, ‘mitochondria’, and ‘unknown’ were excluded from the dataset. OTUs classified as ‘Plantae’, ‘Protista’, ‘Animalia’, ‘Protozoa’, and ‘unknown’ were removed from the ITS dataset. The resulting OTU and taxonomy tables were analyzed using R v3.4.4.

Decomposition analysis and statistics

Samples were rarefied to obtain an even number of sequences per sample. The vegan v2.4.4 package (Oksanen et al. 2013) was used to calculate diversity indices and Bray–Curtis dissimilarities and perform principal-coordinates analyses (PCoA). The effect of initial bark and soil treatments on microbial community composition was assessed with global PERMANOVA tests (Anderson 2001) for both bacterial and fungal composition using the ‘adonis’ function in vegan. Group dispersion was then assessed using the ‘betadisper’ function in vegan and differences in distances to group centroids were tested using ANOVA. Here, group was set as the explanatory variable and distance to group centroid was set as the response variable. Tukey HSD was employed for post-hoc analysis.

Relative decomposition was evaluated as proportionate change in mass (Δmass/initial mass). Proportionate change in mass was used to test for HFA. Bark tissue decomposed beneath the canopy of the tree in which it originated (home), while a portion of the tissue decomposed beneath the canopy of a different species (‘away’). The difference in bark decomposition rate at ‘home’ versus ‘away’ was used to calculate the HFA index (Ayres et al. 2009b).

The effect of soil environment on decomposition was determined using analysis of variance (ANOVA). ANOVA was justified by assessing normality visually via kernel density plots in R version 3.3.3 and Levene’s test was employed to assess homogeneity of variance. For ANOVA, tree species (oak and hemlock) and soil environment (‘home’ and ‘away’) were set as explanatory variables, with decomposition set as the response variable. The explanatory variable consisted of four levels, hemlock bark in hemlock soil; hemlock bark in white oak soil; white oak bark in white oak soil; and white oak bark in hemlock soil (Fig. S1). Tukey's HSD was used for post-hoc analysis.

Our analysis was performed in two ways. First, the ANOVA model was decomposed into apriori contrasts using the multcomp package (Hothorn et al. 2013) followed by simultaneous tests of general linear hypotheses (Hothorn et al. 2008). While not directly related to the main goals of this study, contrasts also included the means of ‘eastern hemlock soil environment’ versus ‘white oak soil environment’; means of ‘eastern hemlock bark substrate’ versus ‘white oak bark substrate’. A two-way ANOVA was also performed, since it is a common statistical approach for this type of experimental design. A two-way ANOVA can reveal an interaction, such that the outcome of one factor (i.e. bark type) is dependent on the levels of another factor (i.e. soil type). To assess HFA, a one-sample t-test was performed on each species to determine whether HFA analysis was significant (i.e different from zero).

Hypothesis testing

As with leaf litter, we expected bark would decompose faster within soils near the tree from which it was collected, demonstrating a home-field advantage. If the first hypothesis is true, bark should decompose faster when it is in its “home” soil community, resulting in an HFA index above zero (Fig. 1). We also expect that different bark types will differentially shift microbial assemblages. If the second hypothesis is true, microbial assemblages, will sort according to the identity of the introduced bark, as opposed to the soil regime (Fig. 1).

For the two-way ANOVA, we would interpret a significant interaction between bark type and soil type to mean that bark decomposition is uniquely dependent on the combination of soil type and bark type, which may suggest HFA. A nonsignificant interaction may suggest that decomposition rate is influenced by species bark type, soil type, or both.

Results

Bark HFA

The calculated HFA index indicated a clear “home-field advantage” to decomposition for white oak bark (Fig. 2, white oak, t = 5.1446, df = 24, 95% CI [0.05,0.13], p = 2.883e−05) and a home-field disadvantage for eastern hemlock bark (Fig. 2, eastern Hemlock, t = − 2.928, df = 23, 95% CI [− 0.16, − 0.02], p = 0.007), meaning that bark decomposition was consistently faster in white oak soils. In addition, our soil treatments were significant predictors of bark decomposition (Fig. 3, ANOVA, F3,96 = 5.15, p = 0.002). Hence, when mean bark decomposition in white oak soils was compared to that in eastern hemlock soils, decomposition was greater in white oak soils (ANOVA, Simultaneous Tests for GLH, 95% CI [− 0.33, − 0.03], p = 0.01). Bark type was also a factor in decomposition rates; irrespective of soil type, eastern hemlock bark decomposed faster than white oak bark (S2, ANOVA, a priori contrasts, Simultaneous Tests for GLH, 95% CI [− 0.30, 0.00], p = 0.05). A two-way ANOVA showed that both bark source (F1,96 = 6.151, p = 0.01487) and soil type (F1,96 = 9.316, p = 0.00294) were both significant predictors of decomposition, but this test did not detect an interaction between these two factors (F1,96 = 0.011, p = 0.91614). This indicates decomposition rate is mainly influenced by bark type and soil type, while evidence for HFA is unclear.

Fig. 2
figure 2

Home Field Advantage (HFA) Index. To test if bark preferentially decomposed near to the tree from which it was collected, decomposition of eastern hemlock and white oak bark was observed in ‘home’ and ‘away’ soils. The difference in ‘home’ and ‘away’ decomposition rates is used to formulate the HFA index. Positive values indicate greater mass loss at home. A significant HFA was detected for white oak (One sample t-test, p = 2.883e−05), while a significant home-field disadvantage was detected for eastern hemlock (One sample t-test, p = 0.007)

Fig. 3
figure 3

Bark decomposition in soil communities. Depicted here are boxplots showing the effect of each treatment. Hashmarks within a boxplot represent bark in away soils. The notches in each boxplot represents 95% CI around the mean. The upper edge and lower edge of each boxplot is the 25% and 75% quartile, that make up the interquartile. The whiskers represent 1.5 × the interquartile. Points beyond the whiskers are outliers. Analysis of variance (ANOVA) suggest that these treatments significantly impacted bark decomposition (F3, 96 = 5.15, p = 0.002). Also, decomposition appeared to be lower in eastern hemlock soils. These differences were significant when contrasted (ANOVA, apriori contrast, white oak soils versus eastern hemlock soils 95% CI [− 0.33, − 0.03], p = 0.01). There were also significant differences in decomposition when comparing the two away treatments (Tukey HSD, p < 0.001)

Bark and soil specific clustering of microbial communities

Bark decomposition impacted microbial composition in nearby soil. A PCoA ordination of 16S rRNA gene OTUs based on Bray–Curtis distances showed that bark type influences bacterial and fungal composition, and these differences persisted in the bark-adherent soil, irrespective of soil type (Fig. 4a, PERMANOVA, F5, 74 = 16.01, p < 0.001). Similar patterns were observed for fungi (Fig. 4b, PERMANOVA, F5, 78 = 5.58, p < 0.001). Beta-dispersion analysis demonstrated that only white oak bark in the 16S dataset and hemlock bark in the ITS dataset differed in group dispersion compared to other treatment types (Supplemental Fig. S3; 16S: ANOVA F5-74 = 8.52, p = 2.07e−6; ITS: ANOVA, F5-78 = 10.72, p = 7.31e−8).

Fig. 4
figure 4

Similarities and dissimilarities of microbial communities across treatments. a Bacteria: depicted here are community similarities and dissimilarities among bacterial assemblages in the vicinity of decomposing bark (PERMANOVA, F5, 74 = 16.01, p < 0.001). b Fungi: depicted here are similarities and dissimilarities among fungal assemblages in the vicinity of decomposing bark (PERMANOVA, F5, 78 = 5.58, p < 0.001)

Discussion

In the context of this study, decomposition proceeded more rapidly in white oak soils than in eastern hemlock soils, while hemlock bark decomposed more rapidly than white oak bark overall. With respect to decomposition, the home-field advantage (HFA) hypothesis suggests that in particular points in space, the accumulation of specific litter types will shape the composition and activity of soil communities to promote those that are most efficient at metabolizing the available resources. Other hypotheses, including substrate quality, functional breadth, rhizosphere priming, and Gadgil effect, provide alternative views of the controls on decomposition pathways (Chapman et al. 2003; Hättenschwiler et al. 2011; Sinsabaugh and Shah 2011; Freschet et al. 2012; Sulman et al. 2017; Malik 2019; Fanin et al. 2016). Here we aimed to determine whether HFA could be detected for bark. Although the HFA index suggested preferential decomposition for white oak bark in “home” soils, it was clear that white oak soils also promoted the decomposition of hemlock bark, and we did not observe a significant bark type x soil type interaction based on a two-way ANOVA. Bark type also had a strong influence on the composition of adjacent soil microorganisms, regardless of soil type. Taken together, these findings suggest that the solubilization of bark and/or the introduction of bark-borne microorganisms may influence short-term local succession in microbial assemblages, but without lasting impacts, these trends are less likely to influence the suitability of recalcitrant litter decomposability.

Plant traits can affect bark decomposition

In this study, we used two trees with contrasting traits: softwood, coniferous eastern hemlock, versus hardwood deciduous, white oak. The deposition of litter with varied traits can influence soil microbial abundance, colonization, and enzymatic activity in different ways, depending on the traits of the plants producing litter, which can be shaped by how particular plants interact with the environment (Diedhiou et al. 2009; Aneja et al. 2006; Mukhopadhyay and Joy 2010). For example, plant populations responding to stress from insect pests (e.g. eastern hemlock) have been shown to increase production in polyphenolics (Constabel and Ryan 1998; Bhonwong et al. 2009; Dar et al. 2017), which can impact decomposition. Meanwhile adaptations to frequent fire disturbance (e.g. white oak) may result in thicker outer bark (Rosell 2016), which could reduce wood access to decomposers, and decrease the rate of carbon cycling. Our results suggest that bark decomposition was slower in soils beneath the canopy of eastern hemlock than those of white oak (Fig. 3), essentially creating a HFA for white oak and not eastern hemlock (Fig. 2); however, whether this effect is due to characteristics that are inherent to soils in which these trees grow or to soil conditioning by white oaks cannot be determined here. Bark traits are also likely to have influenced rates of decomposition both directly (hemlock decomposed more quickly than white oak; Figure S2) and through modification of soil microbial composition (Fig. 4a, b).

Environmental factors in decomposition

Various soil traits can impact decomposition. Decomposition was consistently slower in eastern hemlock soils than that of white oak soils, suggesting a systematic difference between the tree-adjacent soils. Soil traits that covary with the preferred habitats of each tree species could overwhelm any potential HFA for bark decomposition. For instance, fine roots that play a role in nutrient, carbon and water cycling; can also impact decomposition (Cheng and Kuzyakov 2005; McCormack et al. 2015), while topographic control of soil moisture distribution (e.g. subsurface flow; slope position) can alter substrate decomposability (Burt and Butcher 1985; Nippgen et al. 2011). Differences in canopy cover can also impact soil moisture content (Duff et al. 1997). It is also possible that the legacy of eastern hemlock needles influenced soil pH, which may have impacted invertebrate activity (Beyer et al. 1987), and rates of decomposition (Ayanaba and Jenkinson 1990).

White oak soils and the functional breadth hypothesis

While this study was designed to test for HFA, there are competing hypotheses, including the functional breadth hypothesis (Keiser et al. 2014), which may provide an alternative explanation. Overall, bark decomposed more readily in white oak soils (Fig. 3), which could indicate that white oak soils have a broader range of functions related to decomposition than hemlock soils. A broader functional breadth of microbial activity is achieved when there is broad spectrum of labile and recalcitrant litter. This may explain why forests are expected to have a broader functional breadth than either grasslands or agroecosystems (Fanin et al. 2016). As it relates to litter deposition with a greater range along the labile–recalcitrant spectrum, hardwood soils may exhibit a greater range than hemlock soil due to dissimilarity in microbial machinery required to degrade the hemicellulose backbone. This explanation is partly based on the higher xylan and lower mannan content of hardwood hemicellulose than that of softwoods (Zhu and Pan 2010). Thus, the main hemicellulases involved in the degradation of hardwood hemicellulose are xylanases, beta‐xylosidases and xyloglucanases (Álvarez et al. 2016). In contrast, mannanases and beta‐mannosidases are preferentially used to degrade the hemicellulose backbone of softwoods (Álvarez et al. 2016). This contrast in enzymatic machinery suggests different modes of decomposition, which likely leads to differences in microbial assemblages and function, which is relevant to this study. At the labile end of the labile-recalcitrant spectrum, hardwoods may also broaden soil community function due to the frequency of broad leaf deposition that is often of higher N content and less acid-forming than hemlock leaf litter (Yates et al. submitted).

Decomposition and litter traits may shape community structure

Bark decomposition may impact local soil environments, at least in the short term. Variation in local microbial composition was explained more by bark type than soil environment (Fig. 4a, b). One explanation is that bark brings a resident microbiome to soil, which may persist through time. Alternatively, soil microbial composition may have been impacted during decomposition by the liberation of bark specific compounds. Dossa et al. (2018) showed that bark decomposition can influence soil faunal assemblages, indicating that environmental changes can be mediated at the local scale by recalcitrant litter at multiple trophic levels. As decomposition progresses, there is a change in a substrate’s rigidity, density and mass (Pyle and Brown 1998), which can reflect nitrogen mineralization and lignocellulose solubilization (Melillo et al. 1982; Pointing et al. 2003). Progressive change in bark decomposition was shown to impact soil microbial composition, which may in turn provide insight into nutrient cycling and carbon storage.

Conclusion

Common history among microbes and fallen litter is foundational to the “home-field advantage” hypothesis (Strickland et al. 2009; Ayres et al. 2009b). Mixed results were observed in this decomposition study. Although we observed higher decomposition of white oak bark in ‘home’ environments, decomposition of hemlock bark was also higher in white oak soils. This type of soil effect could be due either to environmental conditioning by white oaks, or by differences that are inherent to the locations that white oak and hemlock preferentially grow within a forest. Bark type rather than soil environment, was a significant predictor of microbial composition in the adherent soil community. Here we show that the decomposition of recalcitrant litter can influence microbial composition in the short-term, possibly due to an influx of bark-borne microorganisms or the liberation of polyphenolics, anti-microbial tannins and other chemical compounds, but that this may not necessarily lead to sustained advantages in decomposition for ‘home’ bark.