# Effects of Sphagnum Leachate on Competitive Sphagnum Microbiome Depend on Species and Time

Samuel Hamard, Bjorn J M Robroek, Pierre-Marie Allard, Constant Signarbieux, Shuaizhen Zhou, Tongchai Saesong, Flore de Baaker, Alexandre Buttler, Geneviève Chiapusio, Jean-Luc Wolfender, et al.

► To cite this version:

Samuel Hamard, Bjorn J M Robroek, Pierre-Marie Allard, Constant Signarbieux, Shuaizhen Zhou, et al.. Effects of Sphagnum Leachate on Competitive Sphagnum Microbiome Depend on Species and Time. Frontiers in Microbiology, 2019, 10, pp.1-17. 10.3389/fmicb.2019.02042. hal-02377974

HAL Id: hal-02377974

https://hal.science/hal-02377974

Submitted on 24 Nov 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Check for
updates

# Effects of *Sphagnum* Leachate on
# Competitive *Sphagnum* Microbiome
# Depend on Species and Time

Samuel Hamard<fnref n="1" />,<fnref n="2" />,³,⁴, Bjorn J. M. Robroek<fnref n="2" />,<fnref n="5" />,⁶, Pierre-Marie Allard⁷,
Constant Signarbieux<fnref n="2" />,³, Shuaizhen Zhou⁷, Tongchai Saesong⁷,⁸, Flore de Baaker⁷†,
Alexandre Buttler<fnref n="2" />,³,⁹, Geneviève Chiapusio⁹,<fnref n="10" />, Jean-Luc Wolfender⁷, Luca Bragazza<fnref n="2" />,³,¹¹
and Vincent E. J. Jassey<fnref n="1" />,<fnref n="2" />,³*

OPEN ACCESS

**Edited by:**

Javier A. Ceja-Navarro,
Lawrence Berkeley National
Laboratory, United States

**Reviewed by:**

Tomasz Mieczan,
University of Life Sciences of Lublin,
Poland

Talia Jewell,
General Automation Lab
Technologies, Inc., United States

*Correspondence:

Vincent E. J. Jassey
vincent.jassey@univ-tlse3.fr

†Present address:

Flore de Baaker,
Pfizer Manufacturing, Puurs, Belgium

**Specialty section:**

This article was submitted to
Terrestrial Microbiology,
a section of the journal
Frontiers in Microbiology

Received: 31 August 2018

Accepted: 20 August 2019

Published: 06 September 2019

**Citation:**

Hamard S, Robroek BJM,
Allard P-M, Signarbieux C, Zhou S,
Saesong T, de Baaker F, Buttler A,
Chiapusio G, Wolfender J-L,
Bragazza L and Jassey VEJ (2019)
Effects of Sphagnum Leachate on
Competitive Sphagnum Microbiome
Depend on Species and Time.
Front. Microbiol. 10:2042.
doi: 10.3389/fmicb.2019.02042

Plant specialized metabolites play an important role in soil carbon (C) and nutrient
fluxes. Through anti-microbial effects, they can modulate microbial assemblages and
associated microbial-driven processes, such as nutrient cycling, so to positively or
negatively cascade on plant fitness. As such, plant specialized metabolites can be used
as a tool to supplant competitors. These compounds are little studied in bryophytes.
This is especially notable in peatlands where *Sphagnum* mosses can dominate the
vegetation and show strong interspecific competition. *Sphagnum* mosses form carpets
where diverse microbial communities live and play a crucial role in *Sphagnum* fitness
by regulating C and nutrient cycling. Here, by means of a microcosm experiment, we
assessed to what extent moss metabolites of two *Sphagnum* species (S. *fallax* and
S. *divinum*) modulate the competitive *Sphagnum* microbiome, with particular focus on
microbial respiration. Using a reciprocal leachate experiment, we found that interactions
between *Sphagnum* leachates and microbiome are species-specific. We show that
both *Sphagnum* leachates differed in compound richness and compound relative
abundance, especially sphagnum acid derivates, and that they include microbial-related
metabolites. The addition of S. *divinum* leachate on the S. *fallax* microbiome immediately
reduced microbial respiration (–95%). Prolonged exposition of S. *fallax* microbiome to
S. *divinum* leachate destabilized the food web structure due to a modulation of microbial
abundance. In particular, leachate addition decreased the biomass of testate amoebae
and rotifers but increased that of ciliates. These changes did not influence microbial
CO₂ respiration, suggesting that the structural plasticity of the food web leads to its
functional resistance through the replacement of species that are functionally redundant.
In contrast, S. *fallax* leachate neither affected S. *divinum* microbial respiration, nor

<footnote n="1">¹ ECOLAB, Laboratoire d'Ecologie Fonctionnelle et Environnement, Université de Toulouse, CNRS, Toulouse, France,</footnote>

<footnote n="2">² Laboratory of Ecological Systems (ECOS), Ecole Polytechnique Fédérale de Lausanne (EPFL), School of Architecture, Civil and Environmental Engineering, Lausanne, Switzerland, ³ Swiss Federal Institute for Forest, Snow and Landscape Research (WSL), Lausanne, Switzerland, ⁴ Laboratoire de Géologie, UMR 8538, CNRS-ENS, Ecole Normale Supérieure, Paris, France,</footnote>

<footnote n="5">⁵ School of Biological Sciences, University of Southampton, Southampton, United Kingdom, ⁶ Aquatic Ecology and Environmental Biology Group, Faculty of Science, Institute for Water and Wetland Research, Radboud University, Nijmegen, Netherlands, ⁷ School of Pharmaceutical Sciences, University of Geneva, University of Lausanne, Geneva, Switzerland, ⁸ Department of Pharmaceutical Chemistry and Pharmacognosy, Faculty of Pharmaceutical Sciences and Center of Excellence for Innovation in Chemistry, Naresuan University, Phitsanulok, Thailand, ⁹ Laboratoire Chrono-Environnement, Université Bourgogne Franche Comté, UMR CNRS 6249 USC INRA, Montbéliard, France,</footnote>

<footnote n="10">¹⁰ Laboratoire Carttel, Université Savoie Mont Blanc INRA 042, Domaine Universitaire Belledonne, Le Bourget-du-Lac, France, ¹¹ Department of Life Science and Biotechnologies, University of Ferrara, Ferrara, Italy</footnote>
microbial biomass. We, however, found that *S. fallax* leachate addition stabilized the food web structure associated to *S. divinum* by changing trophic interactions among species. The differences in allelopathic effects between both *Sphagnum* leachates might impact their competitiveness and affect species distribution at local scale. Our study further paves the way to better understand the role of moss and microbial specialized metabolites in peatland C dynamics.

**Keywords:** allelopathy and allelochemicals, metabolomics, microbial networks, microbial respiration and biomass, peatland, plant competition, plant-exudates, soil food-web

## INTRODUCTION

Plant species composition and diversity play a fundamental role in regulating ecological processes such as nutrient and carbon (C) fluxes through their linkages with belowground components. Notably, plants are known to put a selective pressure on soil microbes and their activities (Robroek et al., 2015, 2017a; Cúcio et al., 2016; Sánchez-Canñizares et al., 2017), and to drive microbial assemblages in soil (Berg and Smalla, 2009; Schlatter et al., 2015; Tkacz et al., 2015). This selective effect is performed by plant-derived chemicals, either through the amount and form of carbon and other nutrients that plants provide to the soil (Huang et al., 2014; Cline and Zak, 2015; Sasse et al., 2018), or through specialized metabolites - i.e., allelochemicals - that directly interact with microbes (Bertin et al., 2003; Musilova et al., 2016; Niro et al., 2016). These plant-derived chemicals allow the establishment of mutual, competitive and defensive relationships between specific plants and microorganisms (Latif et al., 2017), and may have indirect effects on competing plants. Such effects could arise from either direct phytotoxic or stimulatory effects on the microbial community on which competing plants rely for nutrients uptake. This possibility is exemplified by plants that inhibit mutualistic bacteria or fungi that competitive plants require for optimal growth, but that the donor plant does not need (Cipollini et al., 2012). Allelopathy can also modify plant-microbe interactions by favoring growth of pathogenic or parasitic microbes that harm competitive plants (Cipollini et al., 2012; Geisen et al., 2018). Understanding the effects of plant-derived metabolites on the specific plant microbiome as a competitive mechanism is therefore of key importance to understand the dynamics of ecological processes within ecosystems.

The majority of studies on plant-microbial interactions focus on vascular plant-dominated ecosystems, largely ignoring the importance of cryptogams. Yet, cryptogams occur in many terrestrial habitats and contribute significantly to global ecosystem functions such as nitrogen and C cycle (Elbert et al., 2012). The processes linked to cryptogam metabolites (i.e., regulation of microbial activity, indirect plant competition) in cryptogam-dominated ecosystems remain poorly explored, while they could be rather important for ecosystem functions (Asplund et al., 2013; Asplund and Wardle, 2013). This is especially the case in ombrotrophic peatlands, where cryptogams such as *Sphagnum* mosses can dominate the ecosystem (Yu et al., 2011). *Sphagnum* mosses form expansive carpets that provide a habitat for a large diversity of microbial communities (Gilbert et al., 1998;

Jassey et al., 2013, 2015; Bragina et al., 2014; Mieczan et al., 2015a,b). The *Sphagnum* microbiome is structured in a microbial food web constituted by bacteria, fungi (decomposers), protists (producers, predators, top-predators), and small-sized metazoan (predators, top-predators). The functioning of this *Sphagnum*-associated food web critically determines the cycling of C and nutrients (Jassey et al., 2015) through the microbial loop (Gilbert et al., 1998) –a trophic pathway through which C and nutrients from organic matter are returned to higher trophic levels through their incorporation into bacterial and fungal biomasses. Empirical observations suggest that the structure and activity of the *Sphagnum* microbiome determine *Sphagnum* fitness (Kostka et al., 2016). For instance, microbial functional guilds such as methanotrophic bacteria, photosynthetic protists and nitrogen-fixing cyanobacteria clearly benefit the host-species by providing a source of C and/or nitrogen that enhances moss productivity (Jassey et al., 2015; Kostka et al., 2016; Carrell et al., 2019). Hence, it can be argued that any shift in the composition of the *Sphagnum* microbiome can modulate *Sphagnum* productivity and, ultimately, ecosystem C cycling.

Different species of *Sphagnum* coexist in peatlands. Each species has a specific productivity and performance (Gunnarsson, 2005; Robroek et al., 2007), and each species is associated to a specific microbiome (Opelt et al., 2007; Bragina et al., 2012). It has been observed that *Sphagnum* species often grow in spatially-structured population in response to interspecific competition (Ingerpuu and Vellak, 2013). Apart from environmental effects such as water-table level or niche separation (Bragazza, 1997; Robroek et al., 2007), the mechanisms that drive competition are not well determined. *Sphagnum* exudates can affect the growth of other *Sphagnum* species (Ingerpuu and Vellak, 2013), suggesting that *Sphagnum* metabolites may play a role in interspecific competition. This hypothesis is supported by the fact that *Sphagnum* mosses produce a variety of specialized metabolites (Rasmussen et al., 1995; Opelt et al., 2007; Chiapusio et al., 2018), such as phenolic acid derivatives (Rudolph and Samland, 1985; Rasmussen et al., 1995), with potential allelopathic effects (Verhoeven and Liefveld, 1997). Phenolic acids, particularly sphagnum acid, are highly water-soluble compounds (Rasmussen et al., 1995) and have been long suspected of antimicrobial effects in peatlands (Verhoeven and Liefveld, 1997; Binet et al., 2017). Recent studies further suggest that polyphenol compounds could be an important factor modulating the structure of microbial assemblages in *Sphagnum* peatlands (Jassey et al., 2011b,c, 2013). Altogether, these studies indicate an important possible role of *Sphagnum*
exudates in affecting peatland microbial activity and suggest their potential implication in interspecific *Sphagnum* competition by alteration of the microbiome, ultimately affecting *Sphagnum* fitness, peatland primary productivity and the peatland C cycle. Understanding the role of *Sphagnum* exudates in driving microbial communities is thus crucial to better understand peatland C dynamics.

Here we focus on two widespread *Sphagnum* species, i.e., *S. fallax* and *S. divinum*, that often co-occur together, to assess whether *Sphagnum* interspecific competition and ecosystem functioning can be indirectly mediated by a shift in their respective microbiome. We tested how each of these two *Sphagnum* species affected the structure and the functioning of the microbial community associated to the competitive species. Using a reciprocal leachate-transfer experiment, we assessed whether a brief exposition to allochthonous *Sphagnum* leachates (i.e., leachates from the competing species) altered the CO₂ respiration of the microbiome. Second, we tested if a prolonged exposition to allochthonous *Sphagnum* leachates altered the microbial food-web structure and functioning, emphasizing on microbial activity and C related processes. We hypothesized that: (1) a short-term exposure to allochthonous leachates will inhibit microbial respiration due to species-specific metabolites inhibitory effects, and that this effect is universal across the two species; (2) a prolonged exposure to allochthonous leachate will alter the structure of microbiome causing a decrease of microbial metabolism.

# MATERIALS AND METHODS

## Field Sampling and Leachate Collection

In April 2015, we collected 15 intact shallow cores (diam. 11 cm; depth 15 cm) of *Sphagnum fallax* and *S. divinum* (30 cores in total) in the Store Mosse National Park, Sweden (57°17′54 N, 14°00′39 E, permit 521-895-2011). Cores were extracted in a habitat that was low in vascular plant cover (<5%), but occasional specimens were carefully removed after which the cores were placed in PVC pipes (hereafter referred as microcosms) that were open from the top and closed at the bottom. Microcosms were then transported to the laboratory facilities at the Ecole Polytechnique Fédérale de Lausanne, Switzerland, where they were kept in a growth chamber (20°C, 70% RH, 14 h/10 h day/night photoperiod, PAR intensity 200 µmol m⁻² s⁻¹) for 2 weeks to acclimate. During the acclimation period, all microcosms were watered daily with 12 mL of artificial rainwater (Garrels and Christ, 1965), which sufficed to keep water levels at field conditions (−1 cm for *S. fallax* and −3 cm for *S. divinum*). During the acclimation period, we collected the leachates from all microcosms. To do so, we first drained the microcosms and added fresh artificial rainwater. This was repeated after 2 days, with the difference that this time all leachate was collected. Leachates from microcosms with the same species were combined in a bulk leachate, resulting in one leachate for *S. fallax* (*L*<sub>SF</sub>) and one for *S. divinum* (*L*<sub>SD</sub>). After collection, both leachates were filtered at 0.2 µm to remove microorganisms, and frozen (−20°C) until utilization.

## Leachate Chemical Characterization

After leachate collection, a set of chemical analyses were performed for both *L*<sub>SF</sub> and *L*<sub>SD</sub> leachates. The concentrations of dissolved organic carbon (DOC) and nitrogen (DON) were quantified by combustion using a Shimadzu analyzer (TOC-V CPH). The quality of DOC was determined by spectroscopy analysis by measuring the absorbance within the range of wavelengths 250–665 nm (Jaffrain et al., 2007). Spectral slopes (S<sub>250-665</sub>, nm⁻¹) were calculated using linear least squares regressions with Ln-transformed absorptions. High S<sub>250-665</sub> values indicate low molecular weight material and/or decreasing aromaticity of the leachate (Hansen et al., 2016). Leachate phenolic content was determined using the Folin and Ciocalteu's reagent with gallic acid as standard (Jassey et al., 2011b). Phenolic content was expressed as mg of equivalent gallic acid per volume of leachate (mg eq. gallic acid·L⁻¹).

A detailed chemical analysis of the metabolites present in *Sphagnum* leachates, *Sphagnum* mosses and bog water was performed using Ultra-High Pressure Liquid Chromatography coupled with High Resolution Mass Spectrometry (UHPLC-HRMS; see details below). To disentangle the provenance of the metabolites (*Sphagnum* versus peat), several comparisons were made. First, we compared leachate metabolic composition to water and methanolic extracts of lyophilized *S. fallax* and *S. divinum* (Jassey et al., 2011a). These extracts were prepared by grinding 0.05 g dry weight of lyophilized *Sphagnum* (0–3 cm; capitulum) in 5 mL of water (mili-Q) or methanol (80/20 methanol/distilled water v/v) using metal ball grinder. We also compared leachate metabolic profiles to the one of the bog water as comparison. *S. fallax*, *S. divinum* and bog water were sampled in the Store Mosse site during the core sampling campaign. All samples were filtered through SPE columns, and directly injected in the UHPLC-HRMS. In total, we analyzed seven metabolic profiles through UHPLC-HRMS: the bog water, two *Sphagnum* leachates, two *Sphagnum*-water extracts, and two *Sphagnum*-methanolic extracts.

## UHPLC-HRMS Analysis

We qualitatively assessed the metabolic composition of both leachates using High-resolution Mass Spectrometry (MS) and molecular networking. Briefly, chromatographic separation was performed on a Waters Acquity UPLC system interfaced to a Q-Exactive Focus mass spectrometer (Thermo Fisher Scientific, Bremen, Germany), using a heated electrospray ionization (HESI-II) source. Thermo Scientific Xcalibur 3.1 software was used for instrument control. The LC conditions were as follows: column, Waters BEH C18 50 × 2.1 mm, 1.7 µm; mobile phase, (A) water with 0.1% formic acid; (B) acetonitrile with 0.1% formic acid; flow rate, 600 µL·min⁻¹; injection volume, 1 µL; gradient, linear gradient of 2–100% B over 6 min and isocratic at 100% B for 0.6 min. An Acquity IClass UPLC photodiode array detector was used to acquire PDA spectra, which were collected from 210 to 450 nm. In positive ion mode, diisooctyl phthalate C₂₄H₃₈O₄ [M + H]⁺ ion (m/z 391.28429) was used as internal lock mass. The optimized HESI-II parameters were as follows: source voltage, 4.0 kV (pos); sheath gas flow rate (N2), 55 units; auxiliary gas flow rate, 15 units; spare gas
flow rate, 3.0; capillary temperature, 275.00°C (pos), S-Lens RF Level, 45. The mass analyzer was calibrated using a mixture of caffeine, methionine-arginine-phenylalanine-alanine-acetate (MRFA), sodium dodecyl sulfate, sodium taurocholate, and Ultramark 1621 in an acetonitrile/methanol/water solution containing 1% formic acid by direct injection. The data-dependent MS/MS events were performed on the three most intense ions detected in full scan MS (Top3 experiment). The MS/MS isolation window width was 1 Da, and the stepped normalized collision energy (NCE) was set to 15, 30, and 45 units. In data-dependent MS/MS experiments, full scans were acquired at a resolution of 35,000 FWHM (at *m/z* 200) and MS/MS scans at 17 500 FWHM both with an automatically determined maximum injection time. After being acquired in a MS/MS scan, parent ions were placed in a dynamic exclusion list for 2.0 s.

## LC-MS Data Processing

ThermoRAW MS data were converted to the open MS format (.mzXML) using the MSConvert software, part of the ProteoWizard package (Chambers et al., 2012). The converted files were then treated using the MzMine 2.36 software suite. The parameters were adjusted as follow: the centroid mass detector was used for mass detection with a noise level set to 1.0E6 for MS level set to 1, and to 0 for MS level set to 2, respectively. The ADAP chromatogram builder was used and set to a minimum group size of scans of 5, minimum group intensity threshold of 1.0E5, minimum highest intensity of 1.0E5 and *m/z* tolerance of 8.0 ppm. For chromatogram deconvolution, the algorithm used was the wavelets (ADAP). The intensity window S/N was used as S/N estimator with a signal to noise ratio set at 25, a minimum feature height at 10,000, a coefficient area threshold at 100, a peak duration ranges from 0.02 to 0.9 min and the RT wavelet range from 0.02 to 0.05 min. Isotopes were detected using the isotopes peaks grouper with a *m/z* tolerance of 5.0 ppm, a RT tolerance of 0.02 min (absolute), the maximum charge set at 2 and the representative isotope used was the most intense. An adduct (Na⁺, K⁺, NH₄⁺, ACN⁺, CH₃OH⁺, Isopropanol⁺) search was performed with the RT tolerance set at 0.1 min and the maximum relative peak height at 500%. A complex search was also performed using [M + H]⁺ for ESI positive mode, with the RT tolerance set at 0.1 min and the maximum relative peak height at 500%. A custom database search was finally performed using the Dictionary of Natural Products 2018 (v. 26.2) database<fnref n="1" />, restricted to moss metabolites. Peak alignment was performed using the join aligner method (*m/z* tolerance at 8 ppm), absolute RT tolerance 0.065 min, weight for *m/z* at 10 and weight for RT at 10. The peak list was gap-filled with the same RT and *m/z* range gap filler (*m/z* tolerance at 8 ppm). Eventually the resulting aligned peaklist was filtered using the peak-list rows filter option in order to keep only features associated with MS2 scans. Full parameters are available as .xml file as supporting information (Sphagnol_profiles_MzMineparams.mzmine).

## LC-MS Data Analysis: Molecular Networks Generation

In order to further identify compounds in MS dataset, we used the molecular networking (MN) approach that group metabolites by structural similarity (Wang et al., 2016). The MN approach is based on the organization and visualization of tandem MS data through a spectral similarity map, revealing the presence of similar MS fragmentations patterns. As structurally related compounds tend to share similar fragmentation spectra, nodes gathered together create clusters of structural analogs. The generated molecular networks were annotated using experimental spectral libraries (GNPS libraries<fnref n="2" />) and an *in silico* fragmented database of natural products using a previously detailed dereplication strategy (Allard et al., 2016). In the network created, edges were filtered to have a cosine score above 0.65 and more than 6 matched peaks. Further edges between two nodes were kept in the network if and only if each of the nodes appeared in each other's respective top 10 most similar nodes. The spectra in the network were then searched against GNPS' spectral libraries. The library spectra were filtered in the same manner as the input data. All matches kept between network spectra and library spectra were required to have a score above 0.7 and at least 6 matched peaks. The output was visualized using Cytoscape 3.6 software<fnref n="3" />. The nodes of the generated networks were colored according to the number of MS/MS triggered in a given sample, thus offering a semi-quantitative information. The size of the nodes was mapped according to the total sum of precursor ions intensities. Molecular networks are available on the GNPS servers at the following addresses: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=00de48872c8b4d76b108a40fdfb7ea0a and https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=60f9f10962ae4aa0952a54c9eb9fc21d.

## Experimental Setup and Measurements

To assess the allelopathic effects of *Sphagnum* leachate addition on potential neighboring competing species, we set up a targeted reciprocal leachate experiment. We first divided the species-specific mesocosms into two experimental groups. One group, consisting of 10 microcosms (five per species), was used to test the instantaneous effect of allochthonous leachate addition on the *Sphagnum* microbiome. The second group, consisting of 20 microcosms (10 per species), was used to test the effect of prolonged exposure to allochthonous leachates.

### Instantaneous Leachate Effects

We sampled 3 g of *Sphagnum* shoots (0–4 cm depth) from all microcosms and placed them in falcon tubes (n = 5). From these shoots, we extracted the microbiome by shaking at 40 rpm for 1.5 h in 30 mL of Mili-Q water®. The extracted microbiome was recovered by filtration at 1 mm. For each microcosm, 0.8 mL of microbial extract was transferred to a Microresp™ (Campbell et al., 2003) 96-deep-well microplate, allocating 8 wells to each microcosm. Subsequently, we added 0.2 mL

<footnote n="1">¹http://dnp.chemnetbase.com</footnote>

<footnote n="2">²http://gnps.ucsd.edu</footnote>

<footnote n="3">³www.cytoscape.org</footnote>
artificial rainwater solution (+ water) to four of these technically replicated microbiomes while the other four microbiomes received 0.2 mL of allochthonous leachates (+ leachate). On the whole, this resulted in four incubation treatments, i.e., two for *S. fallax* microbiome (SF microbiome + water and SF microbiome + SD leachate), and two for *S. divinum* microbiome (SD microbiome + water and SD microbiome + SF leachate). After substrate addition (rainwater or allochthonous leachates), the 96-deep-well microplate was sealed with a 96-well detection microplate containing agar gel and cresol red as indicator dye (Campbell et al., 2003), and incubated at 20°C in the dark. Discoloration of the indicator gel was measured using spectroscopy at 570 nm at irregular intervals (BioTek SynergyMX). The absorbance values were normalized at a given time by the initial absorbance values, after which the percentage of CO₂ released from each well was calculated (Campbell et al., 2003). Values for the four technical replicates were averaged after data-quality check.

In order to understand the effects of leachate composition on microbial respiration, we focused on the effects of *S. divinum* compounds on *S. fallax* microbial CO₂ respiration. We used a fractionation of methanolic extracts of *S. divinum*. Methanolic extracts presented the advantage to be more concentrated than leachates which made the fractionation possible. 3.27 g of the methanolic extract of *S. divinum* were separated by Flash chromatography. Chromatographic conditions: PuriFlash® C18 HQ column (15 µm particle size, 120G), gradient H₂O (+0.1% FA)/MeOH (+0.1% FA) (25:75-100:0 in 109 min, 100:0 isocratic from 109 to 164 min), 2 µL (5 mg/mL), flow rate 30 mL/min, UV detection 200, 254, 280, and 366 nm, 18 mL fractions. The 136 resulting fractions were pooled in 23 fractions according to their thin layer chromatography profiles. The 23 fractions were dissolved in dimethyl sulfoxide (DMSO) at a concentration of 10 µg.L⁻¹. As a next step, we incubated 0.8 mL of *S. fallax* microbial extract with the 0.2 mL of all obtained *S. divinum* metabolite fractions using an identical approach as aforementioned (Microresp™). Microbial incubation with DMSO was used as a control. Following the same protocol as previously described, we assessed microbial respiration at irregular intervals.

## Prolonged Effect of Leachates

Ten microcosms from each species were randomly divided into two groups. The first group (n = 5) was watered daily with 12 mL artificial rainwater, and hence served as a control (C). The second group received daily 12 mL of allochthonous leachate; *S. fallax* was watered with 12 mL of L~SD~ and *S. divinum* was watered with 12 mL of L~SF~. Essentially this resulted in four experimental treatments: SF-C, SF-L~SD~, SD-C, and SD-L~SF~. The experiment lasted 3 weeks. Due to experimental constraints (i.e., limited leachate availability) we were not able to apply autochthonous leachate addition. The experiment was achieved over a 3 weeks period, which guarantees the microbiome to have turned over multiple times (Schönborn, 1965; Schmidt et al., 2007). Throughout the experiment, the position of the cores in the growth chamber was spatially randomized and we kept the same conditions applied during the acclimation period.

## Microbial biomass and community structure

At the end of the 3 weeks leaching experiment, phospholipids fatty acids (PLFA) biomarkers were used to estimate the biomass of fungi, gram-negative, gram-positive, and actinomycete bacteria (Denef et al., 2009). PLFA were extracted from 0.25 g lyophilized *Sphagnum* shoot over a 12 h period in a solvent phase comprising 3.0 mL 50 mmol.L⁻¹ phosphate buffer (pH 7.0), 3.8 mL chloroform and 7.6 mL methanol (Börjesson et al., 1998). PLFA 19:0 (Larodan, Malmö, Sweden) was added as an internal standard to the phospholipid fraction. PLFA were methylated to form fatty acid methyl esters using 1 mL of 0.2 mol.L⁻¹ methanolic KOH (Sundh et al., 1997; Chowdhury and Dick, 2012) and analyzed on a gas chromatograph coupled to a mass spectrometer. Results were expressed as micrograms of PLFA per gram of *Sphagnum* dry mass (µg PLFA.g⁻¹ DM).

The biomass of predators (bacterivores, fungivores, and omnivores) such as ciliates, testate amoebae, rotifers and nematodes was estimated using inverted microscopy (Olympus IX71, × 400, Utermöhl method). To this aim, we collected 3 g of fresh *Sphagnum* shoots (0-4 cm depth), fixed them in 20 mL of glutaraldehyde (2% final concentration) and stored at 4°C in the dark before analyses. Testate amoebae, ciliates, rotifers, and nematodes were extracted from mosses following the method described in Jassey et al. (2011b). From a 3 mL subsample, we identified microbial predators with the appropriate taxonomic literature (Lynn, 2006). The 3 mL aliquots were settled for 1 h in the counting chamber, which was more than enough according to ciliate sinking rates (Claessens and Prast, 2008). Cells were identified and enumerated across 50% of the total chamber area using transects. The abundance of each species was converted into biovolume (µm³), based on geometrical shapes and dimensions measured under the microscope (length or diameter; width, and height) and then into biomass using conversion factors as given in Gilbert et al. (1998). The biomass data were expressed as micrograms of C per gram of *Sphagnum* dry mass (µgC g⁻¹ DM).

## Microbial respiration and enzyme activity

We extracted the microbial communities from each of the 20 microcosms as described above. We incubated the microbial extracts in MicroResp™ with rain water for 8 h, following the previous protocol. Microbial enzyme activity was quantified in microplates following Jassey et al. (2016). Fluorescence of fluorescein diacetate (FDA), a proxy for total enzyme activity (Green et al., 2006), was monitored spectrophotometrically with an excitation wavelength of 365 nm and emission detection at 450 nm (BioTek, SynergyMX). Spectrophotometric measurements were made every hour during incubation until a plateau was reached. We based the calculation of enzyme activities on incubation times of 3 h for FDA hydrolysis. Methodological controls containing boiled enzyme extracts were further assayed. Overall enzyme activity was then calculated by subtracting the mean fluorescence of boiled controls from the mean fluorescence of extracts wells. FDA enzyme activity was converted into µmoles per gram dry weight per min (µmol min⁻¹ g⁻¹ DM).
## Statistical Analyses

All statistical analyses were performed in R version 3.5.0 (R Core Team, 2018). Analysis of variance (ANOVA) was applied to test the effects of *Sphagnum* species, leachate treatments and their interaction (fixed effects) on CO₂ release, microbial enzyme activity and microbial biomasses. Prior analyses, we checked the normality and the homoscedasticity of the data; the data were log-transformed when necessary. For time-related measurements (kinetics of microbial respiration) a linear mixed model was used with time, species and treatment as fixed variables and specific microcosm in function of time as a random effect (Pinheiro and Bates, 2000). Similar models were used to assess the effect of each *S. divinum* fraction on microbial respiration. Differences among metabolite relative abundance in both leachates were tested using chisquared-test.

## Food Web Numerical Analyses

### Food Web Constructions

To elucidate the effect of prolonged leachate addition on microbial interactions within the *Sphagnum* microbiome, we used a network approach based on the force of trophic interactions among microbial species. For each treatment (i.e., SF-C, SF-L<sub>SD</sub>, SD-C, SD-L<sub>SF</sub>), we built an average food web based on the observations from the five replicates. These networks were built in four steps. Firstly, we assigned every species or generic group assessed with PLFA (i.e., gram-positive and gram-negative bacteria, actinomycete, other bacteria, saprophytic fungi, arbuscular mycorrhiza, and diatoms) to a trophic group: decomposers, primary producer, consumer and top-predator. We then identified the feeding habit of consumers and top-predators (**Supplementary Table S1**) using microscope observation and literature (Gilbert et al., 2000, 2003; Mieczan, 2009; Wilkinson and Mitchell, 2010; Jassey et al., 2012; personam observations). Based on these feeding habits, we constructed a general table of hypothetical trophic links between microbial species and/or groups. We obtained two tables of hypothetical links, one associated to *S. fallax* microbiome, the other associated to *S. divinum* microbiome. Secondly, we used basic rules to transform the general tables of hypothetical links into microcosm-specific tables of effective trophic links. A link between a predator and a prey was considered to be effective in a microcosm (i) if the link existed in the general table of hypothetical links, (ii) if the predator and the prey coexisted in microcosm, and (iii) if the predator was less abundant than the prey. Based on previous observations in peatlands (Jassey et al., 2013, 2015; Reczuga et al., 2018), we assumed that predators were always less abundant than bacteria, fungi and algae quantified by PLFA. Thirdly, we weighted the effective trophic links between a predator and a prey by the relative abundance of the predator amongst microcosms of a *Sphagnum* species to take into account the intensity of interactions. For instance, an effective link between a predator and a prey was considered more intense in microcosm A than in microcosm B if the predator was more abundant in microcosm A compared to microcosm B, both microcosms belonging to the same *Sphagnum* species. Finally, we built an average table of links for each treatment by averaging

the microcosm-specific tables of weighted links according to their treatments. We only kept links with a weight higher than 0.2 in order to remove weaker interactions (**Supplementary Figure S1**), assuming that they were unlikely. Such threshold resulted in the removal of 20% of weaker links in *S. fallax* and *S. divinum* microcosms.

### Food Web Structure Analyses

Networks were then produced using the *igraph* R package (Csárdi and Nepusz, 2006). We extracted the core properties of the webs to evaluate whether substructures responded to leachate addition according to Ma et al. (2019). Each network was then analyzed in terms of connectivity (C), density of links within the network (edge D), core size and density of links within the core (Phir). Connectivity was calculated as the number of links divided by the square of the number of species in the network. Density of links within the network was calculated as the number of effective links divided by the number of total possible links. Core size and density of links within the core were calculated as in Ma et al. (2019). Indices of network beta diversity were also calculated between networks (Poisot et al., 2012) using betalink R package (Poisot et al., 2016). Especially, we calculated network beta diversity between average and hypothetical networks (constructed with hypothetical trophic links as explained above) to see if treatments altered the distance with hypothetical networks.

### Food Web Robustness Analyses

We tested the robustness and specificity of each network using a series of null models (Robroek et al., 2017b; Ma et al., 2019). We used three scenarios to simulate new networks. The first scenario (hereafter R<sub>I</sub>) comprised 1000 randomizations of trophic links among microbes whilst keeping the total number of links and species within network intact. The second scenario (hereafter R<sub>R</sub>) comprised the removal of trophic links of two species taken randomly, and is based on the observed effect of leachate addition on food web structures (see Results for details). Species removal was repeated until the whole combinations of species removal was reached (i.e., 231 combinations). The last scenario (hereafter R<sub>IR</sub>) was a combination between R<sub>I</sub> and R<sub>R</sub>: 30 combinations of two species were randomly chosen based on our observations. For each combination, the links of the two species were removed before proceeding to 100 randomizations of the trophic links within the networks as described above (3000 reshuffled networks in total). On each simulated networks, we calculated structural indices (beta diversity, C, edgeD, Coresize, and Phir) and compared them with the structural indices found in original networks from leachate addition treatments using standardized effect size as described in Robroek et al. (2017b). Finally, we proceeded to a targeted species removal within networks based on species that dramatically changed their network position between control and leachate addition treatments (*Hyalosphenia papillo* and *Hyalosphenia elegans* in *S. fallax* networks and *Assulina musorum* and *Hyalosphenia papillo* in *S. divinum* networks). Starting from the control networks, we proceeded to a target removal of these species when they were brought
to network periphery in leachate addition treatments. On the contrary, we proceeded to a replacement of their effective links by their hypothetic links when they were brought to the heart of networks in leachate addition treatments. Again, we measured structural network indices and compared them with the observed structural indices.

# RESULTS

## Global Chemical Characterization of Sphagnum Leachate

Leachate composition was species-specific with global compounds being more concentrated in L\_SF than in L\_SD: water-soluble phenolic concentration was higher in L\_SF (4.87 mg L⁻¹) than in L\_SD (3.92 mg L⁻¹) as well as DOC (L\_SF = 39 mg L⁻¹; L\_SD = 11.55 mg L⁻¹ in) and DON (L\_SF = 1.95 mg L⁻¹; L\_SD = 0.89 mg L⁻¹). The spectral slope (S250-665) of DOC was lower in L\_SF (0.067) than in L\_SD (0.144), which indicates that the aromaticity and/or molecular weight of L\_SF was higher than L\_SD.

According to the metabolite composition, we found that L\_SF and L\_SD contained almost the same chemical compounds but in different relative proportions (Figure 1A). L\_SD (524 metabolites) was slightly richer in metabolites than L\_SF (516 metabolites); 2% of these leachate-metabolites being species-specific, some of them being present only in L\_SD (Supplementary Figure S2). Despite similar metabolite richness (Figure 1B), nearly half of the more concentrated compounds (relative abundance > 2.5%) differed between L\_SF and L\_SD (P < 0.01, Chi-test; Figure 1A). Furthermore, only 25% of leachate-metabolites were common with bog water metabolites. This comparison indicates that bog water and Sphagnum-leachates composition is different in terms of molecular compounds. The relative abundance of shared metabolites between Sphagnum leachates and bog water also differed (P < 0.01, Chi-test) (Figure 1A). Finally, we found according to the molecular network (Supplementary Figure S3) that Sphagnum leachates were closely related to water and methanolic Sphagnum extracts, which indicates that the metabolites from Sphagnum leachates are mostly released by Sphagnum mosses and do not come from the peat or vascular plants.

Focusing on metabolites common to methanolic/water extractions and leachates profiles, several annotated metabolites corresponded to aminoacids derivatives. Phenylalanine for example was found to be present in both Sphagnum species extract (water, methanolic) and leachates but not in bog water. Typical Sphagnum metabolites such as sphagnum acid were found to be present in water and methanolic extracts but not in leachates (Figure 2A). We however detected a derivate of sphagnum acid, sphagnum acid methyl ester, in L\_SD but not L\_SF. Finally, we also detected bacteria-related metabolites such as aminobacteriohopane and bacteriohopanetetrol derivatives (Figure 2B). 35-aminobacteriohopane-32,33,34-triol was found to be present in methanolic, water extracts of both species, in the bog water, but only in L\_SD. Another one [32,33,34,35-Bacteriohopanetetrol (21βH,32R,33R,34S)-form 35-O-(6-Amino-6-deoxy-β-D-glucopyranoside)], was detected

in leachates, and water/methanolic extracts of S. divinum but only observed in the S. fallax methanolic extract (Figure 2B).

## Instantaneous Effect of an Allochthonous Leachate Addition on Microbial Respiration

Allochthonous leachates addition effects on Sphagnum microbiome CO₂ respiration are time and species-specific [F(1, 173) = 3.73, P = 0.05]. The respiration of S. fallax-associated microbiome was strongly inhibited by L\_SD addition and reached only 5% of the control rates within the first 30 h [F(1, 8) = 10.38, P = 0.018; Figure 3A]. This inhibition was, however, time-limited so that after 46 h of incubation microbial respiration recovered and even released 78% more CO₂ than controls [F(1, 8) = 3.8, P = 0.10]. We did not find an effect of L\_SF on the microbial respiration of S. divinum microcosms [F(1, 81) = 0.56, P = 0.48; Figure 3B].

The incubation of S. fallax associated microbiome with S. divinum metabolic fractions as substrate confirmed the inhibition of microbial respiration on short-term. After 6 h of incubation, S. divinum metabolic fractions reduced S. fallax microbial respiration by 12% on average compared to the DMSO control (Figure 3C). We found that 19 fractions out of 23 induced an inhibition effect on microbial respiration, although only the fraction F11 was significant (-36%, t = 2.81, P = 0.047). Furthermore, the effect of S. divinum fractions varied over time [F(23, 240) = 1.77, P = 0.019]. Indeed, the fractions that had minimal effect after 6 h had a positive effect on microbial respiration after 21 h: F2 + 15% (t = 2.53, P = 0.045), F16 + 37% (t = 7.48, P < 0.001), F17 + 52% (t = 4.00, P = 0.019), F18 + 17% (t = 3.14, P = 0.021), F19 + 20% (t = 4.11, P = 0.008), F20 + 24% (t = 5.138, P = 0.004). Overall, 11 fractions out of 23 showed an inhibitory effect after 21 h of incubation, in particular the fractions F11 (-25%, t = 2.71, P < 0.05) and F9 (-6%, t = 2.81, P < 0.05).

## Effect of a Prolonged Allochthonous Leachate Addition on Microbial Food-Web Structure and Function

By the end of the 3 weeks leachate experiment, we did not find differences in CO₂ respiration between the microbial communities that had been exposed every day to allochthonous leachate (L\_SD and L\_SF treatment) and those that had been watered only with rain water (C treatment). This result was apparent in the S. divinum microbiome [F(1, 8) = 1.22, P = 0.3] and the S. fallax microbiome [F(1, 8) = 2.68, P = 0.14; Supplementary Figure S4]. Similarly, we did not find any significant effect of treatment on overall enzyme activity (FDA) for both Sphagnum species (Supplementary Figure S5).

The analysis of the Sphagnum microbiome community composition showed a clear difference in terms of species diversity and species biomass between S. fallax and S. divinum microcosms. The first two axis of the PCA (Supplementary Figure S6) showed that Hyalosphenia papilio (mixotrophic testate amoebae), Amphileptus sp. (ciliates), Lecane quadridentata and Polyarthra sp. (rotifers) were dominating in S. fallax
<figure><img src="image_3.png" /></figure>

microbiome, while *Hyalosphenia elegans* (heterotrophic testate amoebae), *Colurella obtusa* (rotifers) and nematodes dominated the *S. divinum* microbiome. We found that 3 weeks of allochthonous leachate addition tended to alter the composition of microbial communities in *S. fallax* microcosms, but not in *S. divinum* microcosms (**Figure 4**). In the *S. fallax* microbiome, we observed a decrease in biomass of some predators, so that the biomass of testate amoebae decreased by 45% [*F*(1, 8) = 4.44, *P* = 0.07] mainly due to a decrease of mixotrophic testate amoebae [−50%, *F*(1, 8) = 3.63, *P* = 0.09]. Additionally, the biomass of rotifers decreased by 42% [*F*(1, 8) = 2.47, *P* = 0.16], although these trends were not significant.

The addition of allochthonous leachates caused changes in the food web structure, both in *S. fallax* and *S. divinum* microbiome

(Figure 5). In the *S. fallax* microbiome, leachate addition led to a 30% decrease in connectivity and link density within the network (Figure 6). Similarly, the networks' core size was reduced from seven species in the SF-C networks to four species in the SF-L<sub>SD</sub> networks. Interestingly, two of the three species removed from the core were mixotrophs (*Hyalophenia papilio* and *Amphitrema wrightianum*). They were not only pulled away to the periphery of the cores, but they also lost the totality of their trophic links. We further found that SF-L<sub>SD</sub> network was more different with the hypothetical network than SF-C network as network beta diversity was higher in the first case than in the second one (Figure 6). In *S. divinum* microbiome, we found an opposite effect of allochthonous leachate addition. In particular, we found that core size was doubled in SD-L<sub>SF</sub> networks (10 species)
**A** Specialized metabolites of the bryosphere

<figure><img src="image_4.png" /></figure>

**B** Bacterial-related metabolites within the bryosphere

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Clusters individualized from the molecular network (see **Supplementary Figure S3**) and showing identified specialized metabolics found in *Sphagnum* leachates, extracts, and bog water. Metabolites are associated to the bryosphere (A) or to bacteria (B).</figcaption></figure>
<figure><img src="image_6.png" /></figure>

**FIGURE 3** | Response of microbial respiration to short-term allochthonous leachate addition. The CO₂ released in each well is calculated by normalizing with the initial value and combined according to the treatment. *S. fallax* microbiome (A) were either incubated with water (lightgreen) or with L~SD~ leachate (darkgreen). Similarly, *S. divinum* microbiome (B) was either incubated with water (lightred) or with L~SF~ (darkred). Error bars refer to standard error. (C) Shows the difference in respiration of *S. fallax* microbiome after incubation with different fractions of L~SD~ compared to DMSO only (control). Differences are expressed as percentage. Blue bars show differences of respiration after 6 h of incubation and red bars after 21 h of incubation. Error bars refer to standard error and asterisks indicate significant differences with control for each incubation time (*P* < 0.05).

compared to SD-C networks (five species). Two species, i.e., *Amphitrema musorum* and *Hyalophenia papilio*, that did not have any link in the SD-C network were further integrated to the core. The increasing of core size led to a decreasing of edge density within the core (0.11 and 0.40 in SD-L~SF~ and SD-C networks respectively). Allochthonous leachate addition did not affect neither connectivity and edge density of the networks, nor network beta diversity.

The networks obtained after randomization and/or random species loss strongly differed from SF-L~SD~ and SD-L~SF~ networks, both in term of network beta diversity and structural indices (SES were strongly different from 0, *P* < 0.001, **Supplementary Figure S7**). This result was found both in *S. fallax* and *S. divinum* microbiomes, irrespective of the null models tested (**Supplementary Figure S7**). This indicates that leachate-induced shifts in food web structures are due neither to random changes in interactions nor to random species loss. Finally, the target

removal of *H. papilio*, *H. elegans* or both species from SF-C network did not lead to an increased similarity with SF-L~SD~ network (**Supplementary Figure S8**). On the same way, the addition of hypothetical links of *A. musorum*, *H. papilio* or both species in SD-C network did not lead to an increased similarity with SD-L~SF~ network (**Supplementary Figure S8**).

DISCUSSION

*Sphagnum* specialized metabolites have long been suspected to drive microbial processes in peatlands (Verhoeven and Liefveld, 1997; Fenner and Freeman, 2011). In this study, we not only show that *Sphagnum* specialized metabolites comprise a mosaic of microbial and *Sphagnum* compounds, but also that these complex assemblages trigger a set of important changes in the *Sphagnum*-microbiome structure and its functioning. We found
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Structure of microbiome in *S. fallax* (top) and *S. divinum* (bottom) microcosms after 3 weeks of daily exposition to rain water (SF-C, SD-C) or allochthonous leachates (SF-LSD, SD-LSF). For each *Sphagnum* species, the two graphs on the left show the predator biomass with TP, total predators, i.e., the sum of HTA, heterotrophic testate amoebae; MTA, mixotrophic testate amoebae; C, ciliates; N, nematodes; R, rotifers. The two graphs on the right show the relative abundance of bacteria and fungi. No significant differences were detected. Error bars refer to standard error.</figcaption></figure>

that the effects of allochthonous *Sphagnum* leachate addition on competitive *Sphagnum*-microbiome were species-specific with antagonistic responses between *Sphagnum* species. While the addition of *S. divinum* leachates significantly reduced microbial CO₂ respiration and altered the structure of *S. fallax*-associated food web, *S. fallax* leachate addition had only slight effects on *S. divinum* microbiome with no effects on microbial respiration. These results are even more striking by considering that both *Sphagnum* leachates mostly contain similar chemical compounds but in different proportions, and that the few *Sphagnum*-specific metabolites were poorly concentrated. Furthermore, we found that *Sphagnum* microbiomes strongly differed between species, suggesting that *Sphagnum* leachates target specific microbial species. Species-specificity has been shown as an important factor in plant-plant allelopathic interactions (Inderjit Wardle et al., 2011). Indeed, plant species are known to produce different cocktail of specialized metabolites according to their

genus and environmental conditions, which can greatly influence their allelopathic interactions. Our results suggest that similar allelochemical mechanisms may exist between *Sphagnum* species and microorganisms. Although we cannot totally exclude that other factors than *Sphagnum* metabolites caused the observed differences in both *Sphagnum* microbiomes, the effect of *Sphagnum* metabolic fractions (**Figure 3C**) strongly suggests that *Sphagnum* metabolites can be a major driver of microbial CO₂ respiration in peatlands.

## Microbial CO₂ Respiration

Previous studies found that specialized compounds produced by *Sphagnum* mosses can inhibit bacterial growth (Mellegård et al., 2009) and microbial activity (Freeman et al., 2001; Fenner and Freeman, 2011). In this study, we show that such anti-microbial effects can be found in *Sphagnum* leachates but their magnitude is species-specific and their effects limited in time.
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Food web structure in SF-C, SF-L<sub>SD</sub>, SD-C, and SD-L<sub>SF</sub> microbiomes. Each point represent a microbial species or group and is colored according to its belonging to the core (blue) or to the periphery (yellow) of the network. Trophic links were weighted by the abundance of predators and the 20% weaker links were removed from the network. Abbreviations of names are specified in **Supplementary Table S1**.</figcaption></figure>

Allochthonous leachate addition caused immediate decrease of *S. fallax* microbial respiration since the first exposure, while *S. fallax* leachate addition did not influence *S. divinum* microbial respiration. This result indicates that *S. divinum* leachates exhibit anti-microbial properties that promptly reduce catabolic activities of microbes associated to *S. fallax*, whereas the opposite is not the case. The analysis of *S. divinum* metabolic fractions on the respiration of *S. fallax* microbiome further demonstrates

that not all *S. divinum* compounds inhibit microbial respiration. This result is supported by recent findings (Mellegård et al., 2009; Chiapusio et al., 2018) and suggests that anti-microbial effects of *Sphagnum* leachates depend both on metabolic cocktail and on the specific concentration of specialized metabolites (Chiapusio et al., 2018).

The general similarity in the composition of *S. divinum* and *S. fallax* leachates has, however, to be relativized by their low
<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Structural indices of SF-C (left, light green), SF-L<sub>SD</sub>, (left, dark green), SD-C (right, light red), and SD-L<sub>SF</sub> (right, dark red) food web networks. Structural indices comprised connectance (C), edge density within the network (EdgeD), core size, and edge density within the core (Phir). Beta diversity was also calculated with hypothetical networks. Each index was normalized by its overall range.</figcaption></figure>

concentration. We cannot exclude that some compounds were under the limit of detection of the UHPLC-HRMS. Nevertheless, slight differences with potential important consequences were found in the composition of both leachates. For instance, sphagnum acid methyl ester – a derivate from sphagnum acid known for its anti-microbial effect (Verhoeven and Liefveld, 1997) was only found in L<sub>SD</sub>. Bacteria-related metabolites belonging to the hopane group of pentacyclic triterpenoids and associated to *Rhodopseudomonas* sp. (van Winden et al., 2012) were also only identified in L<sub>SD</sub>. Molecules of that group are known to have cytotoxic and antibacterial properties (Nagumo et al., 1991; Rohmer et al., 1991). This result is important and indicates that the metabolites secreted by *S. divinum*-associated bacteria could have negative effects on *S. fallax* microbiome. Although further research is needed to specifically identify the compounds with anti-microbial effects in the *Sphagnum*-sphere, our findings bring evidences that *Sphagnum* leachates result from the combined exudation and secretion of metabolites from *Sphagnum* and associated microorganisms. The composition of the resulting chemical assemblage is likely to determine leachate allelopathic potential, and might explain why *S. divinum* leachate has a higher allelopathic effect than *S. fallax* leachate.

Our results indicate that the inhibitory effect of *Sphagnum*-sphere specialized metabolites on microbial respiration in *S. fallax* microbiome is time limited and disappears after 2 days. Following this inhibition period, microbial respiration recovered at an important rate and the amount of CO₂ respired overtook the one that was respired in control plots. This finding suggests that after the loss of inhibition, the nutrients and the DOC contained in *S. divinum* leachate might have stimulated microbial activity and respiration by supplying easily degradable C to microbes (Robroek et al., 2016). We found that prolonged addition of *S. divinum* leachates on *S. fallax* microbiome did not change microbial respiration

despite structural, but functionally redundant changes among the microbial network. This result provides nuance to the traditional assumption that the accumulation of *Sphagnum* specialized metabolites such as polyphenols interferes with microbial catabolic activity in peatlands (Freeman et al., 2001). In our study, microbial enzyme activity and respiration seems to be resilient to prolonged accumulation of *Sphagnum* specialized metabolites due to a selection among microbial species. This suggests that *Sphagnum* leachates may play a role in natural microbial selection (Callaway et al., 2005; Inderjit Wardle et al., 2011) and that microbial adaptation within the community can lead to an alleviation of allelopathic effects (Li et al., 2015). Such alleviation effects might have important consequences for the peatland C balance and suggest that a critical re-examination of the mechanisms driven by *Sphagnum* specialized metabolites in peatland C dynamics are urgently needed.

## *Sphagnum Microbiome and Microbial Networks*

Our results show that *S. fallax* and *S. divinum* microbiomes differ and respond differently to allochthonous leachate addition to such an extent that a much stronger impact on *S. fallax* microbiome was observed compared to *S. divinum* microbiome. In both species, allochthonous leachate addition altered the food web structures, but in a different way. In *S. fallax* microbiome, L<sub>SD</sub> altered the structure of the network through a decreasing of connectivity, edge density and core size. These alterations indicate a destabilization of the food web, which is reinforced by the fact that L<sub>SD</sub> addition leads to a disparity with the hypothetical network. In fact, L<sub>SD</sub> disrupts the intensity of trophic interactions between microbial species through a modulation of microbial abundance and, especially, a decrease
of testate amoebae and rotifer abundance. More precisely, we found that mixotrophic testate amoebae, which combine both photosynthesis and predation, were dramatically isolated from the network in SF-L<sub>SD</sub> treatments. This suggests that the duality in their nutrition does not provide them an advantage to resist to anti-microbial metabolites. It further indicates that *S. divinum* leachates may influence the survival and/or reproduction of mixotrophic testate amoebae, explaining why these species are poorly represented in *S. divinum* microbiome. García-Palacios et al. (2016) reported similar inhibitory effect of vascular plant polyphenols on specific microbial functional groups such as nematodes in soils. Our findings further point to a novel mechanism by which *Sphagnum* metabolites can have anti-predator effects, in addition to the already recognized anti-bacterial and anti-fungal effects (Verhoeven and Liefveld, 1997; Mellegård et al., 2009; Binet et al., 2017; Chiapusio et al., 2018). In the *S. divinum* microbiome, we found that *S. fallax* leachate addition also altered the food web structure, but with opposite effects compared to *S. fallax*. Particularly, we found that connectivity and core size increased in SD-L<sub>SF</sub> treatment, which expresses a greater stability in the network (Csete and Doyle, 2004; Liu et al., 2011; Csermely et al., 2013). Therefore, it suggests that L<sub>SF</sub> did not stress *S. divinum* microbiome but rather benefited and stabilized *S. divinum* associated food web. Those changes in the organization and the structure of *Sphagnum* food webs show that they are structurally unstable and that a perturbation such as allochthonous leachate addition can quickly change them toward a novel state with novel interactions and species. Simulated targeted removal or addition of the more impacted species in term of abundance or connectivity from control networks gave very dissimilar networks than those found with leachate addition. These findings suggest that, in addition to the direct loss or gain of trophic links, allochthonous leachates strongly modulate the trophic interactions among non-affected species through cascading effects. However, as mentioned above, the structural changes do not lead to important functional alterations of the microbial food web due to functional redundancies among microbes. These findings suggest that *Sphagnum* associated microbial food web functioning is somehow resistant to allochthonous *Sphagnum* metabolites thanks to structural plasticity. This further indicates that structural plasticity of microbial food webs tends to alleviate the effects of allochthonous leachate on its functioning on long-term. The underlying biochemical mechanisms still have to be elucidated further, but our results underline that notions of perturbation durations are primordial in the comprehension of their effect to peatland C dynamics.

## Implications for Competition

Plant competitiveness is intimately linked with plant fitness and allocation of resources to different metabolic processes (growth, reproduction, defense, etc.) (Herms and Mattson, 1992). In *Sphagnum*, plant fitness strongly depends on the associated microbiome since microbial activity drives nutrient cycling and production of plant-assimilable nutrients (Weston et al., 2015; Kostka et al., 2016). Hence, by their effects on microbial network structure and functions, allochthonous

leachates probably alter *Sphagnum* fitness. For instance, the diminution of mixotrophic testate amoebae in *S. fallax* microbiome following *S. divinum* leachate addition could have decreased *Sphagnum* C uptake (Jassey et al., 2015). Further, the inhibition of microbial catabolic activity immediately after *S. divinum* leachate addition might have temporarily reduced nutrient availability within the *Sphagnum*-sphere. Over a longer time, the functional redundancy of microorganisms suggested that the microbial catabolic activities linked with C cycle were not affected by leachate addition. However, other microbial functions could have been impacted, especially those related to nitrogen or phosphate cycling (Carrell et al., 2019). This hypothesis is supported by Conkle and White (2012) who also found an immediate time-limited inhibition of microbial CO₂ respiration by antibiotics in wetland soils with a durable alteration of microbial N₂O respiration and phosphatase activity on longer term. As a consequence, an alteration of *Sphagnum* microbiome could impact *Sphagnum* competitiveness by affecting *Sphagnum* fitness. The increased leachate-resistance of *S. divinum* associated microbiome to allochthonous leachate suggests a better competitiveness than *S. fallax*. This result could have important consequences as it has been showed that hollow species (such as *S. fallax*) lose competitive strength under global warming (Breeuwer et al., 2008) so that the expected changes in *Sphagnum* distribution under global changes might be intensified by allelopathic mechanisms.

## CONCLUSION

Our findings have important implications for understanding how soil microbial communities respond to plant-associated chemical exudates. Although the exact nature of the metabolites interfering with microbial species still remains unknown, we refined the composition of *Sphagnum* leachates and extracts and showed that they can include microbial products. Further, our data show that addition of allochthonous leachates destabilize microbial network structure, primarily through anti-predators effects, with cascade consequences on microbial functioning and respiration. More research is needed to understand the exact genetic and physiologic mechanisms that define the observed response of microbial species to specialized metabolites. Our data reveal that microbiome alteration can be an important process in interspecific *Sphagnum* competition and indubitably in the C dynamics of peatlands. These findings suggest that facing the vegetation changes in peatlands due to climate warming and drought, specialized metabolites could significantly intervene in defense or offensive mechanisms through their effects on microbial communities.

## AUTHOR CONTRIBUTIONS

VJ, BR, and SH designed and implemented the microcosm study with the help of CS. SH, BR, and VJ sampled the
*Sphagnum* cores in the field and collected the leachates. SH and VJ collected the samples at the end of the experiment with the assistance of BR. SH performed MicroResp, characterized general characteristics of leachates and quantified enzyme activities with the assistance of VJ, BR, LB, and AB. VJ analyzed the microbial samples under the microscope. P-MA, SZ, and TS performed LC-MS analyses, peak data processing and molecular networks on *Sphagnum* leachates with assistance of GC and J-LW. SH and VJ performed statistical analyses and interpreted the data with assistance of BR, P-MA, LB, and AB. SH and VJ wrote the manuscript to which all authors contributed with discussions and text.

## FUNDING

BR was supported by a NWO-VENI grant awarded by The Netherlands Organization for Scientific Research (NWO-ALW: 863.10.014), by the British Ecological Society (small grant SR17/1427), and by the Dutch Foundation for the Conservation of Irish bogs. This research has been funded by the Swiss National

Science Foundation through the project SPHAGNOL (Grant No. 315260 149807 to LB) and supported by MIXOPEAT project (Grant No. ANR-17-CE01-0007 to VJ) from the French National Research Agency.

## ACKNOWLEDGMENTS

We are indebted to Länsstyrelsen i Jönköpings län and the staff of the Store Mosse National Park, particularly Arne Andersson, Dan Käll, Carina Härlin, and Johan Rova for granting site access (permission 521-895-2011) and make use of the infrastructure of the park. We are thankful to the three reviewers for their valuable comments on our manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02042/full#supplementary-material

## REFERENCES

Allard, P.-M., Péresse, T., Bisson, J., Gindro, K., Marcourt, L., Pham, V. C., et al. (2016). Integration of Molecular Networking And In-Silico MS/MS fragmentation for natural products dereplication. *Anal. Chem.* 88, 3317–3323. doi: 10.1021/acs.analchem.5b04804

Asplund, J., Bokhorst, S., and Wardle, D. A. (2013). Secondary compounds can reduce the soil micro-arthropod effect on lichen decomposition. *Soil Biol. Biochem.* 66, 10–16. doi: 10.1016/j.soilbio.2013.06.013

Asplund, J., and Wardle, D. A. (2013). The impact of secondary compounds and functional characteristics on lichen palatability and decomposition. *J. Ecol.* 101, 689–700. doi: 10.1111/1365-2745.12075

Berg, G., and Smalla, K. (2009). Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. *FEMS Microbiol. Ecol.* 68, 1–13. doi: 10.1111/j.1574-6941.2009.00654.x

Bertin, C., Yang, X., and Weston, L. A. (2003). The role of root exudates and allelochemicals in the rhizosphere. *Plant Soil* 256, 67–83. doi: 10.1023/A:1026290508166

Binet, P., Rouifed, S., Jassey, V. E. J., Toussaint, M.-L., and Chiapusio, G. (2017). Experimental climate warming alters the relationship between fungal root symbiosis and *Sphagnum* litter phenolics in two peatland microhabitats. *Soil Biol. Biochem.* 105, 153–161. doi: 10.1016/j.soilbio.2016.11.020

Börjesson, G., Sundh, I., Tunlid, A., and Svensson, B. H. (1998). Methane oxidation in landfill cover soils, as revealed by potential oxidation measurements and phospholipid fatty acid analyses. *Soil Biol. Biochem.* 30, 1423–1433. doi: 10.1016/s0038-0717(97)00257-5

Bragazza, L. (1997). *Sphagnum* niche diversification in two oligotrophic mires in the southern alps of italy. *Bryologist* 100, 507–515. doi: 10.2307/324413

Bragina, A., Berg, C., Cardinale, M., Shcherbakov, A., Chebotar, V., and Berg, G. (2012). *Sphagnum* mosses harbour highly specific bacterial diversity during their whole lifecycle. *ISME J.* 6, 802–813. doi: 10.1038/ismej.2011.151

Bragina, A., Oberauner-Wappis, L., Zachow, C., Halwachs, B., Thallinger, G. G., Müller, H., et al. (2014). The *Sphagnum* microbiome supports bog ecosystem functioning under extreme conditions. *Mol. Ecol.* 23, 4498–4510. doi: 10.1111/mec.12885

Breeuwer, A., Heijmans, M. M. P. D., Robroek, B. J. M., and Berendse, F. (2008). The effect of temperature on growth and competition between *Sphagnum* species. *Oecologia* 156, 155–167. doi: 10.1007/s00442-008-0963-968

Callaway, R. M., Ridenour, W. M., Laboski, T., Weir, T., and Vivanco, J. M. (2005). Natural selection for resistance to the allelopathic effects of invasive plants. *J. Ecol.* 93, 576–583. doi: 10.1111/j.1365-2745.2005.00994.x

Campbell, C. D., Chapman, S. J., Cameron, C. M., Davidson, M. S., and Potts, J. M. (2003). A rapid microtiter plate method to measure carbon dioxide evolved from carbon substrate amendments so as to determine the physiological profiles of soil microbial communities by using whole soil. *Appl. Environ. Microbiol.* 69, 3593–3599. doi: 10.1128/AEM.69.6.3593-3599.2003

Carrell, A. A., Kolton, M., Glass, J. B., Pelletier, D. A., Warren, M. J., Kostka, J. E., et al. (2019). Experimental warming alters the community composition, diversity, and N2 fixation activity of peat moss (*Sphagnum fallax*) microbiomes. *Glob. Change. Biol.* 25, 2993–3004. doi: 10.1111/gcb.14715

Chambers, M. C., Maclean, B., Burke, R., Amodei, D., Ruderman, D. L., Neumann, S., et al. (2012). A cross-platform toolkit for mass spectrometry and proteomics. *Nat. Biotechnol.* 30, 918–920. doi: 10.1038/nbt.2377

Chiapusio, G., Jassey, V. E. J., Bellvert, F., Comte, G., Weston, L. A., Delarue, F., et al. (2018). *Sphagnum* species modulate their phenolic profiles and mycorrhizal colonization of surrounding andromeda polifolia along peatland microhabitats. *J. Chem. Ecol.* 44, 1146–1157. doi: 10.1007/s10886-018-1023-1024

Chowdhury, T. R., and Dick, R. P. (2012). Standardizing methylation method during phospholipid fatty acid analysis to profile soil microbial communities. *J. Microbiol. Methods* 88, 285–291. doi: 10.1016/j.mimet.2011.12.008

Cipollini, D., Rigsby, C. M., and Barto, E. K. (2012). Microbes as targets and mediators of allelopathy in plants. *J. Chem. Ecol.* 38, 714–727. doi: 10.1007/s10886-012-0133-137

Claessens, M., and Prast, M. (2008). Concentration of fixed plankton samples via settling: how long is long enough? *J. Plankton Res.* 30, 57–64. doi: 10.1093/plankt/fbm095

Cline, L. C., and Zak, D. R. (2015). Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession. *Ecology* 96, 3374–3385. doi: 10.1890/15-0184.1

Conkle, J. L., and White, J. R. (2012). An initial screening of antibiotic effects on microbial respiration in wetland soils. *J. Environ. Scie. Health Part A* 47, 1381–1390. doi: 10.1080/10934529.2012.672315

Csárdi, G., and Nepusz, T. (2006). The igraph software package for complex network research. *Inter. J. Comp. Syst.* 1695, 1–9.

Csermely, P., London, A., Wu, L.-Y., and Uzzi, B. (2013). Structure and dynamics of core/periphery networks. *J. Complex Netw.* 1, 93–123. doi: 10.1093/comnet/cnt016
Csete, M., and Doyle, J. (2004). Bow ties, metabolism and disease. *Trends Biotechnol.* 22, 446-450. doi: 10.1016/j.tibtech.2004.07.007

Cúcio, C., Engelen, A. H., Costa, R., and Muyzer, G. (2016). Rhizosphere microbiomes of european seagrasses are selected by the plant, but are not species specific. *Front. Microbiol.* 7:440. doi: 10.3389/fmicb.2016.00440

Denef, K., Roobroeck, D., Manimel Wadu, M. C. W., Lootens, P., and Boeckx, P. (2009). Microbial community composition and rhizodeposit-carbon assimilation in differently managed temperate grassland soils. *Soil Biol. Biochem.* 41, 144-153. doi: 10.1016/j.soilbio.2008.10.008

Elbert, W., Weber, B., Burrows, S., Steinkamp, J., Büdel, B., Andreae, M. O., et al. (2012). Contribution of cryptogamic covers to the global cycles of carbon and nitrogen. *Nat. Geosci.* 5, 459-462. doi: 10.1038/ngeo1486

Fenner, N., and Freeman, C. (2011). Drought-induced carbon loss in peatlands. *Nat. Geosci.* 4, 895-900. doi: 10.1038/ngeo1323

Freeman, C., Ostle, N., and Kang, H. (2001). An enzymic "latch" on a global carbon store. *Nature* 409:149. doi: 10.1038/35051650

Garcia-Palacios, P., Shaw, E. A., Wall, D. H., and Hättenschwiler, S. (2016). Temporal dynamics of biotic and abiotic drivers of litter decomposition. *Ecol. Lett.* 19, 554-563. doi: 10.1111/ele.12590

Garrels, R. M., and Christ, C. L. (1965). *Minerals, Solutions, and Equilibria*. Manhattan, NY: Harper & Row.

Geisen, S., Mitchell, E. A. D., Adl, S., Bonkowski, M., Dunthorn, M., Ekelund, F., et al. (2018). Soil protists: a fertile frontier in soil biology research. *FEMS Microbiol. Rev.* 42, 293-323. doi: 10.1093/femsre/fuy006

Gilbert, D., Amblard, C., Bourdier, G., André-Jean, F., and Mitchell, E. A. D. (2000). Le régime alimentaire des thécamoebiens (Protista. Sarcodina). *L'Année Biol.* 39, 57-68. doi: 10.1016/S0003-5017(00)80001-X

Gilbert, D., Amblard, C., Bourdier, G., and Francez, A.-J. (1998). The microbial loop at the surface of a peatland:structure, function, and impact of nutrient input. *Microb. Ecol.* 35, 83-93. doi: 10.1007/s0024899 00062

Gilbert, D., Mitchell, E. A. D., Bourdier, G., and Francez, A.-J. (2003). Population dynamics and food preferences of the testate amoeba *Nebela tincta major*-bohemica-collaris Complex (Protozoa) in a Sphagnum Peatland. *Acta Protozool.* 42, 99-104.

Green, V. S., Stott, D. E., and Diack, M. (2006). Assay for fluorescein diacetate hydrolytic activity: optimization for soil samples. *Soil Biol. Biochem.* 38, 693-701. doi: 10.1016/j.soilbio.2005.06.020

Gunnarsson, U. (2005). Global Patterns of Sphagnum productivity. *J. Bryol.* 27, 269-279. doi: 10.1179/174328205X70029

Hansen, A. M., Kraus, T. E. C., Pellerin, B. A., Fleck, J. A., Downing, B. D., and Bergamaschi, B. A. (2016). Optical properties of dissolved organic matter (DOM): effects of biological and photolytic degradation. *Limnol. Oceanogr.* 61, 1015-1032. doi: 10.1002/lno.10270

Herms, D. A., and Mattson, W. J. (1992). The dilemma of plants: to grow or defend. *Q. Rev. Biol.* 67, 283-335. doi: 10.1086/417659

Huang, X.-F., Chaparro, J. M., Reardon, K. F., Zhang, R., Shen, Q., and Vivanco, J. M. (2014). Rhizosphere interactions: root exudates, microbes, and microbial communities. *Botany* 92, 267-275. doi: 10.1139/cjb-2013-2225

Inderjit Wardle, D. A., Karban, R., and Callaway, R. M. (2011). The ecosystem and evolutionary contexts of allelopathy. *Trends Ecol. Evol.* 26, 655-662. doi: 10.1016/j.tree.2011.08.003

Ingerpuu, N., and Vellak, K. (2013). Growth depends on neighbours: experiments with three Sphagnum L. species. *J. Bryol.* 35, 27-32. doi: 10.1179/1743282012Y. 0000000034

Jaffrain, J., Gérard, F., Meyer, M., and Ranger, J. (2007). Assessing the quality of dissolved organic matter in forest soils using ultraviolet absorption spectrophotometry. *Soil Sci. Soc. Am. J.* 71, 1851-1858. doi: 10.2136/sssj006. 0202

Jassey, V. E. J., Chiapusio, G., Binet, P., Buttler, A., Laggoun-Défarge, F., Delarue, F., et al. (2013). Above- and belowground linkages in *Sphagnum peatland*: climate warming affects plant-microbial interactions. *Global Chang. Biol.* 19, 811-823. doi: 10.1111/gcb.12075

Jassey, V. E. J., Chiapusio, G., Gilbert, D., Buttler, A., Toussaint, M.-L., and Binet, P. (2011a). Experimental climate effect on seasonal variability of polyphenol/phenoloxidase interplay along a narrow fen-bog ecological gradient in Sphagnum fallax. *Global Chang. Biol.* 17, 2945-2957. doi: 10.1111/j. 1365-2486.2011.02437.x

Jassey, V. E. J., Chiapusio, G., Mitchell, E. A. D., Binet, P., Toussaint, M.-L., and Gilbert, D. (2011b). Fine-scale horizontal and vertical micro-distribution patterns of testate amoebae along a narrow fen/bog gradient. *Microb. Ecol.* 61, 374-385. doi: 10.1007/s00248-010-9756-9759

Jassey, V. E. J., Gilbert, D., Binet, P., Toussaint, M.-L., and Chiapusio, G. (2011c). Effect of a temperature gradient on Sphagnum fallax and its associated living microbial communities: a study under controlled conditions. *Can. J. Microbiol.* 57, 226-235. doi: 10.1139/W10-116

Jassey, V. E. J., Chiapusio, G., Gilbert, D., Toussaint, M.-L., and Binet, P. (2012). Phenoloxidase and peroxidase activities in Sphagnum-dominated peatland in a warming climate. *Soil Biol. Biochem.* 46, 49-52. doi: 10.1016/j.soilbio.2011.11. 011

Jassey, V. E. J., Lamentowicz, M., Bragazza, L., Hofsommer, M. L., Mills, R. T. E., Buttler, A., et al. (2016). Loss of testate amoeba functional diversity with increasing frost intensity across a continental gradient reduces microbial activity in peatlands. *Eur. J. of Protistol.* 55, 190-202. doi: 10.1016/j.ejop.2016. 04.007

Jassey, V. E. J., Signarbieux, C., Hättenschwiler, S., Bragazza, L., Buttler, A., Delarue, F., et al. (2015). An unexpected role for mixotrophs in the response of peatland carbon cycling to climate warming. *Sci. Rep.* 5:6931. doi: 10.1038/srep16931

Kostka, J. E., Weston, D. J., Glass, J. B., Lilleskov, E. A., Shaw, A. J., and Turetsky, M. R. (2016). The Sphagnum microbiome: new insights from an ancient plant lineage. *New Phytol.* 211, 57-64. doi: 10.1111/nph.13993

Latif, S., Chiapusio, G., and Weston, L. A. (2017). "Chapter Two - Allelopathy and the Role of Allelochemicals in Plant Defence," in *Advances in Botanical Research How Plants Communicate with their Biotic Environment*, ed. G. Becard (Cambridge MA: Academic Press), 19-54. doi: 10.1016/bs.abr.2016.12. 001

Li, Y.-P., Feng, Y.-L., Chen, Y.-J., and Tian, Y.-H. (2015). Soil microbes alleviate allelopathy of invasive plants. *Sci. Bull.* 60, 1083-1091. doi: 10.1007/s11434- 015-0819-817

Liu, Y.-Y., Slotine, J.-J., and Barabási, A.-L. (2011). Controllability of complex networks. *Nature* 473, 167-173. doi: 10.1038/nature10011

Lynn, D. H. (2006). The Sphagnum Ponds of Simmelried in Germany: a Biodiversity hot-spot for microscopic organisms. *Protozoological Monographs*, by M. Kreutz & W. Foissner. *J. Eukaryot. Microbiol.* 53, 543-543. doi: 10.1111/j.1550-7408.2006.00137.x

Ma, A., Lu, X., Gray, C., Raybould, A., Tamaddoni-Nezhad, A., Woodward, G., et al. (2019). Ecological networks reveal resilience of agro-ecosystems to changes in farming management | nature ecology & evolution. *Nat. Ecol. Evol.* 3, 260-264. doi: 10.1038/s41559-018-0757-752

Mellegård, H., Stalheim, T., Hormazabal, V., Granum, P. E., and Hardy, S. P. (2009). Antibacterial activity of Sphagnum acid and other phenolic compounds found in Sphagnum papillosum against food-borne bacteria. *Lett. Appl. Microbiol.* 49, 85-90. doi: 10.1111/j.1472-765X.2009.02622.x

Mieczan, T. (2009). Ciliates in Sphagnum peatlands: vertical micro-distribution, and relationships of species assemblages with environmental parameters. *Zool. Studies* 48, 33-48.

Mieczan, T., Michał, N., Adamczuk, M., and Bielańska-Grajner, I. (2015a). Stable isotope analyses revealed high seasonal dynamics in the food web structure of a peatbog. *Int. Rev. Hydrobiol.* 100, 141-150. doi: 10.1002/iroh.20150 1788

Mieczan, T., Niedźwiecki, M., and Tarkowska-Kukuryk, M. (2015b). Effects of rotifers, copepods and chironomid larvae on microbial communities in peatlands. *Eur. J. Protistol.* 51, 386-400. doi: 10.1016/j.ejop.2015.06.010

Musilova, L., Ridl, J., Polivkova, M., Macek, T., and Uhlik, O. (2016). Effects of Secondary plant metabolites on microbial populations: changes in community structure and metabolic activity in contaminated environments. *Int. J. Mol. Sci.* 17:1205. doi: 10.3390/ijms17081205

Nagumo, A., Takanashi, K., Hojo, H., and Suzuki, Y. (1991). Cytotoxicity of bacteriohopane-32-ol against mouse leukemia L1210 and P388 cells in vitro. *Toxicol. Lett.* 58, 309-313. doi: 10.1016/0378-4274(91)90042-90045

Niro, E., Marziaioli, R., De Crescenzo, S., D'Abrosca, B., Castaldi, S., Esposito, A., et al. (2016). Effects of the allelochemical coumarin on plants and soil microbial community. *Soil Biol. Biochem.* 95, 30-39. doi: 10.1016/j.soilbio.2015. 11.028

Opelt, K., Berg, C., Schönmann, S., Eberl, L., and Berg, G. (2007). High specificity but contrasting biodiversity of Sphagnum-associated bacterial and
plant communities in bog ecosystems independent of the geographical region. *ISME J.* 1, 502–516. doi: 10.1038/ismej.2007.58

Pinheiro, J. C., and Bates, D. M. (2000). "Linear mixed-effects models: basic concepts and examples," in *Mixed-Effects Models in Sand S-Plus. Statistics and Computing*, eds J. Chamber, W. Eddy, W. Hardle, S. Sheather, and L. Tierney (New York, NY: Springer).3–51.

Poisot, T., Baiser, B., Dunne, J. A., Kéfi, S., Massol, F., Mouquet, N., et al. (2016). mangal - making ecological network analysis simple. *Ecography* 39, 384–390. doi: 10.1111/ecog.00976

Poisot, T., Canard, E., Mouillot, D., Mouquet, N., and Gravel, D. (2012). The dissimilarity of species interaction networks. *Ecol. Lett.* 15, 1353–1361. doi: 10.1111/ele.12002

R Core Team, (2018). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing. Available at: http://www.R-project.org/

Rasmussen, S., Wolff, C., and Rudolph, H. (1995). Compartmentalization of phenolic constituents in *Sphagnum. Phytochemistry* 38, 35–39. doi: 10.1016/
0031-9422(94)00650-I

Reczuga, M. K., Lamentowicz, M., Mulot, M., Mitchell, E. A. D., Buttler, A., Chojnicki, B., et al. (2018). Predator–prey mass ratio drives microbial activity under dry conditions in *Sphagnum peatlands*. *Ecol. Evol.* 8, 5752–5764. doi: 10.1002/ece3.4114

Robroek, B. J. M., Albrecht, R. J. H., Hamard, S., Pulgarin, A., Bragazza, L., Buttler, A., et al. (2016). Peatland vascular plant functional types affect dissolved organic matter chemistry. *Plant Soil* 407, 135–143. doi: 10.1007/s11014-015-2710-2713

Robroek, B. J. M., Jassey, V. E. J., Beltman, B., and Hefting, M. M. (2017a). Diverse fen plant communities enhance carbon-related multifunctionality, but do not mitigate negative effects of drought. *R. Soc. open sci.* 4:170449. doi: 10.1098/
rsos.170449

Robroek, B. J. M., Jassey, V. E. J., Payne, R. J., Martí, M., Bragazza, L., Bleeker, A., et al. (2017b). Taxonomic and functional turnover are decoupled in European peat bogs. *Nat. Commun.* 1161, 1–9. doi: 10.1038/s41467-017-01350-5

Robroek, B. J. M., Jassey, V. E. J., Kox, M. A. R., Berendsen, R. L., Mills, R. T. E., Cécillon, L., et al. (2015). Peatland vascular plant functional types affect methane dynamics by altering microbial community structure. *J. Ecol.* 103, 925–934. doi: 10.1111/1365-2745.12413

Robroek, B. J. M., Limpens, J., Breeuwer, A., and Schouten, M. G. C. (2007). Effects of water level and temperature on performance of four *Sphagnum* mosses. *Plant Ecol.* 190, 97–107. doi: 10.1007/s11258-006-9193-9195

Rohmer, M., Bisseret, P., and Sutter, B. (1991). "The hopanoids, bacterial triterpenoids, and the biosynthesis of isoprenic units in prokaryotes," in *Progress in Drug Research/Fortschritte der Arzneimittelforschung/Progrès des Recherches Pharmaceutiques Progress in Drug Research/Fortschritte der Arzneimittelforschung/Progrès des recherches pharmaceutique*, eds J. A. Salmon, L. G. Garland, B. D. Hoyle, J. W. Costerton, N. Seiler, D. Raeburn, et al. (Basel: Birkhäuser Basel), 271–285. doi: 10.1007/978-3-0348-7139-6_6

Rudolph, H., and Samland, J. (1985). Occurrence and metabolism of *Sphagnum* acid in the cell walls of bryophytes. *Phytochemistry* 24, 745–749. doi: 10.1016/
S0031-9422(00)84888-84888

Sánchez-Cañizares, C., Jorrín, B., Poole, P. S., and Tkacz, A. (2017). Understanding the holobiont: the interdependence of plants and their microbiome. *Curr. Opin. Microbiol.* 38, 188–196. doi: 10.1016/j.mib.2017.07.001

Sasse, J., Martinoia, E., and Northen, T. (2018). Feed Your Friends: do plant exudates shape the root microbiome? *Trends Plant Sci.* 23, 25–41. doi: 10.1016/
j.tplants.2017.09.003

Schlatter, D. C., Bakker, M. G., Bradeen, J. M., and Kinkel, L. L. (2015). Plant community richness and microbial interactions structure bacterial communities in soil. *Ecology* 96, 134–142. doi: 10.1890/13-1648.1

Schmidt, S. K., Costello, E. K., Nemergut, D. R., Cleveland, C. C., Reed, S. C., Weintraub, M. N., et al. (2007). Biogeochemical consequences of rapid microbial turnover and seasonal succession in soil. *Ecology* 88, 1379–1385. doi: 10.1890/06-0164

Schönborn, W. (1965). Untersuchungen über die zoolchorellen-symbiose der hochmoor-testaceen. *Limnologica* 3, 173–176.

Sundh, I., Nilsson, M., and Borga, P. (1997). Variation in microbial community structure in two boreal peatlands as determined by analysis of phospholipid fatty acid profiles. *Appl. Environ. Microbiol.* 63, 1476–1482.

Tkacz, A., Cheema, J., Chandra, G., Grant, A., and Poole, P. S. (2015). Stability and succession of the rhizosphere microbiota depends upon plant type and soil composition. *ISME J.* 9, 2349–2359. doi: 10.1038/ismej.2015.41

van Winden, J. F., Talbot, H. M., Kip, N., Reichart, G.-J., Pol, A., McNamara, N. P., et al. (2012). Bacteriohopanepolyl signatures as markers for methanotrophic bacteria in peat moss. *Geochim. Cosmochim. Acta* 77, 52–61. doi: 10.1016/j.gca.
2011.10.026

Verhoeven, J. T. A., and Liefveld, W. M. (1997). The ecological significance of organochemical compounds in *Sphagnum. Acta Bot. Neerl.* 46, 117–130. doi: 10.1111/plb.1997.46.2.117

Wang, M., Carver, J. J., Phelan, V. V., Sanchez, L. M., Garg, N., Peng, Y., et al. (2016). Sharing and community curation of mass spectrometry data with global natural products social molecular networking. *Nat. Biotechnol.* 34, 828–837. doi: 10.1038/nbt.3597

Weston, D. J., Timm, C. M., Walker, A. P., Gu, L., Muchero, W., Schmutz, J., et al. (2015). *Sphagnum* physiology in the context of changing climate: emergent influences of genomics, modelling and host–microbiome interactions on understanding ecosystem function. *Plant, Cell Environ.* 38, 1737–1751. doi: 10.1111>pce.12458

Wilkinson, D. M., and Mitchell, E. A. D. (2010). Testate amoebae and nutrient cycling with particular reference to soils. *Geomicrobiol. J.* 27, 520–533. doi: 10.1080/01490451003702925

Yu, Z., Beilman, D. W., Frolking, S., MacDonald, G. M., Roulet, N. T., Camill, P., et al. (2011). Peatlands and their role in the global carbon cycle. *Eos Trans. AGU* 92, 97–98. doi: 10.1029/2011EO120001

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Hamard, Robroek, Allard, Signarbieux, Zhou, Saesong, de Baaker, Buttler, Chiapusio, Wolfender, Bragazza and Jassey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.