# Characterization of Bacterial and Fungal Community Dynamics by High-Throughput Sequencing (HTS) Metabarcoding during Flax Dew-Retting

Christophe Djemiel, Sebastien Grec, Simon Hawkins

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

## To cite this version:

Christophe Djemiel, Sebastien Grec, Simon Hawkins. Characterization of Bacterial and Fungal Community Dynamics by High-Throughput Sequencing (HTS) Metabarcoding during Flax Dew-Retting. Frontiers in Microbiology, 2017, 8, pp.2052. 10.3389/fmicb.2017.02052. hal-03285503

HAL Id: hal-03285503

https://hal.univ-lille.fr/hal-03285503

Submitted on 13 Jul 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Characterization of Bacterial and Fungal Community Dynamics by High-Throughput Sequencing (HTS) Metabarcoding during Flax Dew-Retting

Christophe Djemiel†, Sébastien Grec† and Simon Hawkins*

Univ. Lille, Centre National de la Recherche Scientifique, UMR 8576 - Unité de Glycobiologie Structurale et Fonctionnelle, Lille, France

## OPEN ACCESS

**Edited by:**
Diana Elizabeth Marco,
National Scientific Council (CONICET),
Argentina

**Reviewed by:**
Alessio Mengoni,
University of Florence, Italy
Ulas Karaoz,
University of California, Berkeley,
United States

*Correspondence:
Simon Hawkins
simon.hawkins@univ-lille1.fr

†Christophe Djemiel
orcid.org/0000-0002-5659-7876
†Sébastien Grec
orcid.org/0000-0003-4143-4035

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

**Received:** 10 July 2017
**Accepted:** 06 October 2017
**Published:** 20 October 2017

**Citation:**
Djemiel C, Grec S and Hawkins S (2017) Characterization of Bacterial and Fungal Community Dynamics by High-Throughput Sequencing (HTS) Metabarcoding during Flax Dew-Retting. Front. Microbiol. 8:2052. doi: 10.3389/fmicb.2017.02052

Flax dew-retting is a key step in the industrial extraction of fibers from flax stems and is dependent upon the production of a battery of hydrolytic enzymes produced by micro-organisms during this process. To explore the diversity and dynamics of bacterial and fungal communities involved in this process we applied a high-throughput sequencing (HTS) DNA metabarcoding approach (16S rRNA/ITS region, Illumina Miseq) on plant and soil samples obtained over a period of 7 weeks in July and August 2014. Twenty-three bacterial and six fungal phyla were identified in soil samples and 11 bacterial and four fungal phyla in plant samples. Dominant phyla were Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes (bacteria) and Ascomycota, Basidiomycota, and Zygomycota (fungi) all of which have been previously associated with flax dew-retting except for Bacteroidetes and Basidiomycota that were identified for the first time. Rare phyla also identified for the first time in this process included Acidobacteria, CKC4, Chlorobi, Fibrobacteres, Gemmatimonadetes, Nitrospirae and TM6 (bacteria), and Chytridiomycota (fungi). No differences in microbial communities and colonization dynamics were observed between early and standard flax harvests. In contrast, the common agricultural practice of swath turning affects both bacterial and fungal community membership and structure in straw samples and may contribute to a more uniform retting. Prediction of community function using PICRUSt indicated the presence of a large collection of potential bacterial enzymes capable of hydrolyzing backbones and side-chains of cell wall polysaccharides. Assignment of functional guild (functional group) using FUNGuild software highlighted a change from parasitic to saprophytic trophic modes in fungi during retting. This work provides the first exhaustive description of the microbial communities involved in flax dew-retting and will provide a valuable benchmark in future studies aiming to evaluate the effects of other parameters (e.g., year-to year and site variability etc.) on this complex process.

**Keywords:** flax dew-retting, bacterial and fungal microbiota dynamics, 16S rRNA and ITS amplicons, metabarcoding, high-throughput sequencing, HTS, CAZyme predictions, trophic modes
## BACKGROUND

Land plants fix ~123 billion tons of carbon per year (Beer et al., 2010) of which an important part becomes channeled into the production of lignocellulosic biomass in plant cell walls (Kuhad and Singh, 1993; Boerjan et al., 2003; Zhou et al., 2011). Soil microflora function as key decomposers in various ecosystems (Soliveres et al., 2016) and are able to degrade this biomass, consisting of lignin and polysaccharide polymers such as cellulose, hemicelluloses, and pectins, by producing a set of synergistically acting hydrolytic enzymes (Warren, 1996; Lynd et al., 2002; Kubicek et al., 2014; Cragg et al., 2015).

During this process, monosaccharides are released and used by microorganisms for energy production thereby contributing to maintenance of the carbon cycle. The microbial diversity associated with this biomass degradation can vary depending on plant cell wall structure and the stage of the decomposition (Akin, 2008; Ventorino et al., 2015; Montella et al., 2017). Microbial dynamics can also vary depending on site location, soil composition, plant species, and biomass architecture (Schneider et al., 2012; Voříšková and Baldrian, 2013; Cardenas et al., 2015; Ventorino et al., 2015). Several investigations have reported that fungal communities change during leaf (e.g., beech, oak, maize) litter decay with an initial predominance of species assigned to the Ascomycota phylum, replaced gradually by Basidyomycota (Schneider et al., 2012; Kuramae et al., 2013; Voříšková and Baldrian, 2013). Bacterial dynamics generally involve changes in the relative proportions of species assigned to Proteobacteria, Actinobacteria, and Bacteroidetes depending on sampling location and wood species to be degraded (Ventorino et al., 2015). The composition and succession of different microbial communities is presumably related to their capacity to degrade and utilize biomass present at a given moment. In this context it is interesting to note that particular microbial CAZymes such as endo- and exo-cellulases, xylanases, pectinases, and peroxidases have been associated with specific ecological groups during plant cell wall decomposition (Eastwood et al., 2011; Zhao et al., 2014; Ventorino et al., 2015).

In this study we investigate the microflora associated with a particularly interesting, and ancient example of human exploitation of microbial lignocellulose degradation known as retting that is believed to date back to the Upper Paleolithic and/or Neolithic (Gübitz and Cavaco-Paulo, 2001; Kvavadze et al., 2009). This process is still used today and constitutes the first step in the industrial separation of long bast fibers from the stems of different fiber species such as flax, hemp, jute, and kenaf (Md. Tahir et al., 2011) used for textiles and composites (Campilho, 2015; Pil et al., 2016). During retting, bast fiber bundles become progressively separated from the surrounding stem tissues, and inter-fiber cohesion is reduced via the action of hydrolytic enzymes produced by straw-colonizing microorganisms (Rosemberg, 1965; Zhang et al., 2005; Md. Tahir et al., 2011; Akin, 2013; Preisner et al., 2014). Retting is performed by either leaving plants on the soil (dew- or field-retting), or by placing them in ponds, rivers, or water tanks (water-retting). Although good quality fibers are produced by water retting this process is more labor intensive and associated with extensive

water pollution. Currently the majority of the world's flax fiber is produced by dew-retting (Akin, 2013; Preisner et al., 2014).

The main challenge during retting is to facilitate fiber cohesion without degrading cellulosic fibers by over-retting (Brown and Sharma, 1984; Akin et al., 1998; Henriksson et al., 1999). Since this process relies on enzymes produced by colonizing microorganisms, a better knowledge of the different groups/species involved should enable a greater understanding and control of this complex process. Although various bacteria and fungi have been identified in a number of different studies by using isolation and culturing approaches (Sharma, 1986a; Henriksson et al., 1997), such a strategy is not powerful enough to obtain a complete inventory of the microorganisms present as only a small percentage of taxa can be successfully cultured under laboratory conditions (Staley and Konopka, 1985; Amann et al., 1995). More recently, molecular tools such as 16S rRNA gene amplification were used to identify new bacteria during bamboo, hemp, and flax retting (Tamburini et al., 2003; Fu et al., 2011; Ribeiro et al., 2015) and 18S rRNA gene amplification was used to identify fungi during hemp retting (Ribeiro et al., 2015). Nevertheless, these approaches are unable to generate an exhaustive inventory of the retting microbiome.

Over the last decade, microbial ecology studies have greatly benefited from the use of high throughput sequencing (HTS) technologies that can produce an exhaustive inventory of bacteria and fungi from complex samples such as soil, litter compost, rumen, and the midgut of cellulosic-feeding insects via targeted-metagenomics (Hirsch et al., 2010; Ihrmark et al., 2012; Suenaga, 2012). These approaches were also used to study plant-microbe interactions (Knief, 2014; Peršoh, 2015) in rhizospheres or endospheres (Lundberg et al., 2012; Bodenhausen et al., 2013; Beckers et al., 2017). Only two studies have reported the use of HTS technologies (Ion Torrent PGM system) on kenaf retting (Visi et al., 2013) and more recently during water-retting of flax (Zhao et al., 2016). In both studies, the microbial community analysis was limited to bacterial domain, despite the importance of fungal taxa in the production of extracellular hydrolytic enzymes (Schneider et al., 2012). In this work, we report the first exhaustive HTS microbial inventory focusing on both bacterial and fungal communities using rRNA amplicon sequencing during dew-retting of flax.

## METHODS

### Experimental Design-Study Site-Sampling

Flax plants (*Linum usitatissimum* L., Cultivar Lorea) were sown on 14 March 2014 near Martainneville (F-27210 Region Hauts-de-France) in the north of France (50°00'03"N and 1°42'27"E). Plants were cultivated and retted on a typical silt loam soil with a neutral/slightly acid pH (INRA Soil Analysis Laboratory, LAS, Arras, France, http://www.lille.inra.fr/las) (Supplementary Table 1). Climatic data during the retting period was obtained from the Abbeville meteorological station at 10 km from Martainneville (infoclimat: http://www.infoclimat.fr/observations-meteo/temps-reel/abbeville/07005.html) (Supplementary Figure 1). "Early" and "standard" flax cultures were pulled (up-rooted) on the 16.07.2014 and 24.07.2014, respectively and dew-retted
in the field until the 25.08.2014 (early cultures) and 05.09.2014 (standard cultures). Replicate straw (plant) and soil samples were collected at regular intervals (R0–R6) during retting from five different locations in the retting field chosen according to a non-systematic W pattern as previously described (Plassart et al., 2012) and shown in **Supplementary Figure 2**. For straw samples, the middle region (30 cm long × total swath height) of the swath was collected; for soil samples, cores (20 cm deep × 8 cm diameter) were used. Stem samples were directly stored at –20°C and soil samples were sieved (pore size <2.0 mm), homogenized and freeze-dried before storage at –80°C.

## DNA Extraction

DNA was extracted from 1 g sample using the GnS-GII (Plassart et al., 2012; Terrat et al., 2012, 2015). Briefly, samples were ground in 15 ml Falcon tubes containing a bead mix (ceramic, silica, and glass) and lysis buffer (100 mM Tris-HCl, pH 8; 100 mM EDTA, pH 8; 100 mM NaCl, 2% w/v and sodium dodecyl sulfate, 2% w/v) in a FastPrep®-24 (MP-Biomedicals, NY, USA) (3 × 30 s at 4, 000 s⁻¹ shaking). Proteins were precipitated by adding 100 µl of KAc (3 M) and nucleic acids recovered by isopropanol precipitation and washed with 70% ethanol, before drying and re-suspension in 100 µl water.

## DNA Purification, Quantification, and Normalization

DNA extracts were filtered through PVPP (PolyVinylPolyPyrrolidone) Micro Bio-Spin® Columns with Bio-Gel® P-6 (Bio-Rad) by a 4 min at 1,000 g, 10°C centrifugation. Collected samples were then purified using the Geneclean Turbo kit (MP-Biomedicals, NY, USA) following the manufacturer's instructions. DNA was quantified on a LightCycler 480 System (Roche) using the Quant-iT™ PicoGreen® dsDNA Assay kit (Invitrogen). Samples were normalized to a concentration of 5 ng/µl and the DNA from the five replicates pooled using the epMotion® 5075 TMX (eppendorf). Altogether, 16 soil samples and 14 stem samples were recovered for further analysis.

## Primers, PCR Amplification, and Sequencing

Bacterial 16S rDNA were amplified using the forward primer S-D-Bact-0341-a-S-17 described by Klindworth (Klindworth et al., 2013) coupled with a customized reverse primer S-D-Bact-0787-a-A-19, based on the 786r primer (Gołębiewski et al., 2014). Fungal ITS regions were amplified using the fITS7 forward primer 5.8S (Ihrmark et al., 2012) and the reverse primer ITS4_KYO1 (Toju et al., 2012; Bokulich and Mills, 2013). All primer sequences are given in Supplementary Table 2.

Amplifications were carried out in a total volume of 40 µl using 5 ng of DNA, 4 µl of 5x HOT FIREPol® Blend Master Mix with 7.5 mM MgCl₂ (Solis Biodyne, Tartu, Estonia), 0.8 µl (0.2 µM) of each primer. PCR1 conditions were: 15 min at 95°C, followed by 30 cycles of 20 s at 95°C, 30 s at 53°C, and 20 s at 72°C, and final elongation for 5 min at 72°C. Single multiplexing was performed using home-made 6 bp indexes that were added to reverse primer during a second PCR2 of 12 cycles using

indexed primers. The resulting PCR2 products were purified by HighPrep™ PCR (Magbio) clean-up system as described by the manufacturer, pooled and loaded onto the Illumina MiSeq cartridge according to the manufacturer instructions for a 2 × 250 bp paired-end sequencing on the GeT-PlaGe Genotoul Platform (INRA Castanet Tolosan, France). The quality of the run was checked internally using PhiX, and then each pair-end sequence was assigned to its sample with the help of the previously integrated index.

## Sequence Processing

A bioinformatic pipeline based on mothur v.1.37.4 (https://github.com/mothur/mothur/releases) (Schloss et al., 2009) was configured to process the bacterial 16S rRNA gene sequences. This pipeline uses the standard Schloss lab operating procedure (http://www.mothur.org/wiki/MiSeq_SOP). Pair-End (PE) FASTQ files were overlapped to form contiguous reads in a single FASTA file with zero differences to the primer sequence and a quality score threshold of 30. Sequences with the following characteristics were removed: ambiguous bases and mismatches, <300/> 500 bp, homopolymers >8 bp, overlap <30 bp. Bacterial sequences were aligned against both SILVA (SSU SILVA 123) and Greengenes (August 2013 release, for input PICRUSt) reference databases. A pre-clustering was done to reduce noise as recommended (Pruesse et al., 2007; Huse et al., 2010) allowing for up to four differences between sequences. Chimeras were detected and removed de novo with the UCHIME (version 4.2) algorithm (Edgar et al., 2011). The clustering of the non-chimeric sequences to Operational Taxonomic Units (OTUs) was done by de novo clustering at 0.03 cut-off of dissimilarity using neighbor based on genomic distance matrix. Finally, a general count sequence table for each OTU of all samples was generated to obtain the consensus taxonomy based on the Ribosomal Database Project's naïve bayesian classifier method (Wang et al., 2007) and for the future OTU-based analysis.

For processing ITS2 from fungal ribosomal ITS sequences, the recently described PIPITS v.1.3.3 pipeline (https://github.com/hsgweon/pipits/releases) was used (Gweon et al., 2015). Raw reads were prepared for ITS extraction and the chosen sub-regions extracted with the ITSx software tool (Bengtsson-Palme et al., 2013) before clustering and taxonomic assignation using the UNITE database (version 31.01.2016) (Abarenkov et al., 2010).

All parameters, algorithms and tools for the bioinformatic steps used in the two pipelines are given in Supplementary Table 3.

The microbial DNA sequencing data sets supporting the results in this article are available at the EBI ENA with accession number PRJEB20299.

## Statistical Analysis

All estimators used to measure the α-diversity and β-diversity were calculated applying mothur procedures following recommendations and parameters suggested by tutorials (Kozich et al., 2013).

Alpha-diversity was estimated with the chao1 non-parametric estimator (Chao, 1984) and evenness was measured with Heip's
estimator [$E_{\text{heip}} = (e^{H'} - 1)/(S-1)$ with $H'$ being Shannon's diversity index and $S$ the number of species] (Heip, 1974). Community diversity was estimated with Shannon's diversity index (Ludwig and Reynolds, 1988) and the inverse Simpson's index (Simpson, 1949). Microbial community coverage was tested by calculating the Good's non-parametric coverage estimator (Good, 1953; Esty, 1986) and verified by rarefaction curves. Differences in alpha diversities were evaluated using the Mann-Whitney-Wilcoxon test.

Beta-diversity was assessed using the Yue and Clayton theta similarity coefficient for community structure and the Jaccard index for community membership (Yue and Clayton, 2005; Barwell et al., 2015).

The non-parametric analysis of molecular variance (AMOVA) (Excoffier et al., 1992) was used to examine the significance of differences between and within different groups (Early vs. Standard and Before vs. After turning swaths) with a *p*-value ≤0.05 being considered as statistically significant.

The diversity indices are computed from a standardized file containing the count of OTUs for each sample.

Spearman rank correlation coefficients were calculated from generated dissimilarity matrices to look for any significant correlations between climatic conditions (temperature and rainfall) and bacterial and fungal community structure.

For population level analyses, several tools were used: Metastats (White et al., 2009) (White et al., 2009), LEfSe (Linear discriminant analysis Effect Size (Segata et al., 2011), and Indicator (from Mothur software).

The PICRUSt v.1.0.0 (https://github.com/picrust/picrust/releases) pipeline (http://picrust.github.io/picrust/) (Langille et al., 2013) was used to predict the functional composition of bacterial enzymatic activity abundance using 16S rDNA datasets. An OTUs table (input file) in BIOM format was generated using Mothur and then reference picked against the Greengenes database. Accuracy of metagenome predictions was controlled by measuring the weighted Nearest Sequenced Taxon Index (NSTI) scores that reflect the availability of reference genomes closely related to the most abundant microorganisms for each sample. To analyze the Carbohydrate Active enZymes (CAZymes) prediction, a pre-calculated table was used (https://sourceforge.net/projects/picrust/files/precalculated_files/).

The FUNGuild v1.0 database (https://github.com/UMNFuN/FUNGuild) was used to assign ecological functions (trophic modes) to each OTUs (Nguyen et al., 2016).

Graphic representations were produced using handmade scripts and based on Highcharts facilities (http://www.highcharts.com/) and jvenn plug-in (Bardou et al., 2014).

# RESULTS

## Metabarcoding and Sequencing

16S rDNA (bacterial) and ribosomal ITS (fungal) amplicons were sequenced using the Illumina MiSeq system. Redesigned primers (Supplementary Table 2) were used in order to avoid potential amplification of plant chloroplastic/mitochondrial DNA. Sequencing generated a very large data set ranging between 103,859 and 279,553 (average = 162,390 ± 33,164) bacterial

raw sequences and between 187,055 and 483,703 (average = 285,070 ± 62,420) fungal raw sequences (Supplementary Tables 4A,B). OTU tables listing all OTUs detected and their abundance normalized by a subsampling are given in Supplementary Table 5 (bacteria) and Supplementary Table 6 (fungi).

## Community Coverage and Diversity

To estimate how representative our samples were of the bacterial and fungal communities Good's coverage estimator was calculated for all samples (Supplementary Tables 6, 7). For bacterial samples Good's coverage values were greater than 99% for all straw samples and between 91 and 92% for soil samples indicating: (i) the high coverage of the sampling community and (ii) that the redesigned reverse primer did not significantly affect V3-V4 bacterial amplification (Supplementary Table 6). For fungal samples, Good's coverage estimators were above 99% for all samples confirming that the population is well-sampled (Supplementary Table 7). These results indicate that the sequencing depth used provides an accurate view of microbial community diversity and were also confirmed by the rarefaction curves (Supplementary Figures 3,4).

To analyze community diversity (alpha diversity) within our microbial samples we calculated Chao1 (species richness), Heip's (species evenness), and Inverse Simpson index metric estimators (Supplementary Figures 5, Supplementary Tables 7–9). For all estimators and all conditions [soil vs. plant (straw), early vs. standard cultures] community diversity was always higher in bacterial samples when compared to fungal samples. For both bacterial and fungal samples all indicators indicated that community diversity was higher in soil samples when compared to plant samples. In contrast, the same indicators revealed no difference in bacterial community diversity (both soil and plant samples) between early vs. standard cultures. For fungal communities, the situation was more complex. While no differences in species richness (Chao1) were observed between early vs. standard plant samples, Heip's estimator values suggested differences in sample evenness. The Inverse Simpson index values also suggested differences in community diversity between these two samples.

Examination of indicator values during the retting period (R0–R6) revealed a range of different profiles suggesting that sample community diversity evolves during this process (Supplementary Figure 5). When only the Inverse Simpson index is taken into account as an overall measure of community diversity (Supplementary Tables 7, 8) all profile types show an overall bimodal form with a peak/trough mainly occurring at R2 (5/8 profiles), but also at R3 (2/8 profiles) and R1 (1/8 profiles). Taken together, these results would suggest that sample community diversity changes (increases or decreases) at some point after R2 (and/or R3). The R2/R3 points are close to the moment when the stem swathes were turned and the observed change in community diversity values might be related to this process. When “early” and “standard” culture sample values are pooled (to provide sufficient data points) analyses shows that there is a significant effect (*p*-value <0.05) of swath turning on bacterial, but not fungal community diversity.
Comparison of soil Chao1, Heip's and Inverse Simpson estimators for the first retting point (R0) with those obtained at the sowing stage (R-1) (Supplementary Figure 5) show that both bacterial and fungal community diversity are always lower at R-1. Such an observation suggests that flax plants modify soil diversity through either a rhizosphere effect and/or the input of other organic material (e.g., leaves). Nevertheless, other abiotic effects (e.g., temperature, soil moisture content) may also have an effect and should not be neglected.

## Community Membership and Structure

To obtain an idea of the beta diversity between our samples we analyzed bacterial/fungal community membership and structure. Principal Coordinate analysis (PCoA) using Jaccard distances (Figure 1) clearly revealed that membership between soil samples and plant samples differed for both bacterial (Figure 1A) and fungal (Figure 1B) communities. While no differences in

community membership could be observed between early vs. standard culture plant samples (AMOVA centroid with *p*-value ≥ 0.05) the presence of two distinct clusters indicated that community membership clearly differed between early vs. standard soil samples. Statistical analyses also indicated that swath turning had a significant effect (AMOVA centroid with *p*-value < 0.05) on both bacterial and fungal community membership of plant samples, but not soil samples.

To analyze community structure we then used the Non-parametric MultiDimensional Scaling (NMDS) ordination of Yue & Clayton dissimilarities to determine distance matrices (Theta YC distances) between all samples (Figure 2). The results show that stress values for both bacterial (0.089) and fungal (0.085) communities are inferior to 0.1 as recommended by Mothur SOP (Standard Operating Procedure, https://www.mothur.org). Overall, and as observed for community membership data, clear differences in community structure occur

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Principal coordinate analysis (PCoA) of bacterial (A) and fungal (B) community membership based on Jaccard distances. Different colors indicate the source samples (light brown, soil early harvest; dark brown, soil standard harvest; light green, plant early harvest; dark green, plant standard harvest). Triangles, samples before swath turning; circles, samples after swath turning; lozenge, soil sample during sowing. (C,D) statistically significant clustering based on AMOVA.</figcaption></figure>
between soil and plant samples (early and standard cultures) for both bacterial (Figure 2A) and fungal (Figure 2B) samples. However, the fungal R0 (early and standard cultures) plant samples form a separate cluster from the other plant samples whereas bacterial R0 plant samples do not. The community structure of early and standard bacterial/fungal soil samples, but not plant samples, is also significantly different (*p*-value < 0.001). As observed for community membership, swath turning also appeared to modify community structure, but not necessarily in the same samples. For bacteria, swath turning had a significant effect on community structure in both soil and plant samples (cf. community membership, significant effect only in plant, but not soil, samples). In contrast for fungi, swath turning only had a significant effect on the community structure of plant, but not soil samples. Calculation of Spearman rank

correlation coefficients indicated that there was no significant correlation between climatic conditions (temperature, rainfall) and community structure (Supplementary Table 10).

## Taxonomic Distribution of Identified Bacteria and Fungi

To evaluate taxonomic distribution of identified bacteria and fungi, OTUs were analyzed to determine consensus taxonomy (Figure 3). Overall more phyla (bacteria and fungi) were present in soil samples when compared to plant samples with 23 (excluding unclassified) bacterial and six fungal phyla in soil samples and 11 bacterial and four fungal phyla in plant samples. Of these phyla, 8 (bacteria) and 2 (fungi) were not previously associated with flax dew-retting in the literature thereby underlining the interest of a metabarcoding approach

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Non-metric multidimensional scaling (NMDS) of bacterial (A) and fungal (B) community structure based on Yue & Clayton distance matrix. Different colors indicate the source samples (light brown, soil early harvest; dark brown, soil standard harvest; light green, plant early harvest; dark green, plant standard harvest). Triangles, samples before swath turning; circles, samples after swath turning; lozenge, soil sample during sowing. (C,D) statistically significant clustering based on AMOVA.</figcaption></figure>
for the identification of new microorganisms. Although the number of phyla identified in soil samples was higher than in plant samples, the most abundant taxa were the same in both cases as might be expected in an analysis at this level (Phyla: Bacteria—Proteobacteria (x̄ 60.64% ± 6.64); Fungi—Ascomycota (x̄ 76.29% ± 4.047) (Supplementary Tables 11, 12). For both bacteria and fungi, the type of culture (early vs. standard) appeared to have little effect on phyla relative abundance, neither in soil nor in plant samples. The relative abundance in both bacterial and fungal soil samples appeared to remain fairly

constant throughout the retting period. In contrast, relative abundance in bacterial plant samples was more dynamic being characterized by a relative increase and/or decrease in percentage relative abundances of Proteobacteria and Bacteroidetes at R2 (Figure 3A). The relative abundance in fungal plant samples appeared to be more stable throughout retting.

Subsequent analyses of plant samples at class level (Figure 4) indicated that the observed increase (Figure 3) in the % relative abundance of the Proteobacteria at R2 was mainly related to a substantial increase (>100%) in the relative

<figure><img src="image_6.png" /></figure>
abundance of the Gammaproteobacteria class in early samples (Figure 4A) correlated with a smaller reduction in relative abundances of the Flavobacteria and Sphingobacteria classes. Similarly, in standard samples the previously observed increase in Proteobacteria (Figure 4) could be related to the increase in Gammaproteobacteria and Betaproteobacteria, coupled with a decrease in relative abundance of Sphingobacteria at R2 (Figure 4B). Additional Proteobacteria peaks were also observed at R3 (Alphaproteobacteria) and R5 (Betaproteobacteria) but had

less overall impact on the Proteobacteria/Bacteroidetes ratio in standard cultures due to an increase in relative abundance of Bacteroidetes classes in latter stages of retting.

Examination of relative abundances of fungal classes revealed a different pattern. Although more classes were identified (8/9 classes in the Ascomycota and Basidiomycota, respectively), the class Dothidiomycetes was by far the dominant class in both early (Figure 4C) and standard (Figure 4D) samples with a relative abundance ranging from a “low” of 40% (R0) and arriving at a

<figure><img src="image_7.png" /></figure>
maximum of 60+ % (R2). The Sordariomycetes were the next most abundant class with a value of between 10 and 15% relative abundance.

The generation of community distance heatmaps (Figure 5) for the top ten bacterial and fungal OTUs provided more detailed information on the different taxonomic groups represented in classes identified in plant samples. For bacteria (Figure 5A), results underlined the abundance of Sphingomonas and Pseudomonas genera. For fungi (Figure 5B), the order Capnodiales represented by Cladosporium herbarum was clearly the most abundant group. When analyzed globally, three different profiles could be identified: (i) the OTU is present throughout the retting period (e.g., OTU00001, Sphingomonas sp. and OTU2685, C. herbarum); (ii) the OTU is present at the beginning and then decreases (e.g., OTU00002, Pseudomonas rhizosphaerae and OTU00006, Pantoea vagans); and (iii) the OTU is absent at the beginning and then increases during retting (e.g., OTU00003, Rhizobium genus, OTU00004, Massilia sp. and OTU1918, Altenaria sp.).

## Community-Level Analysis

Our results suggested that swath turning during retting had a significant effect on the microbial communities. To identify those OTUs most likely to explain differences highlighted by the diversity analyses between before- and after-turning - and that

could therefore represent potential biomarkers of this process, we used three different tests (Metastas, LEfSE, and Indicator). Our results (Supplementary Figure 6, Supplementary Tables 13–15) identified 7/8 bacterial and fungal OTUs, respectively “before,” and 14/6 bacterial and fungal OTUs “after,” swath turning in all three tests. Four and six of these OTUs are present in the top 10 bacterial/fungal OTUs, respectively (Figure 5).

## Hydrolytic Enzyme Potential and Trophic Mode Prediction

Dew-retting of flax straw occurs via the action of hydrolytic enzymes produced by microorganisms and we therefore used PICRUSt software followed by expert curation to predict the bacterial Carbohydrate Active enZyme (CAZy) families potentially present during dew-retting and playing a role in the degradation of cell wall polymers. For all stem samples the NSTI scores were around the 0.15 level considered as acceptable according to PICRUSt instructions. Our results (Figure 6, Supplementary Figure 7) show that a wide range of different enzymes targeting both the backbones and side chains of the major polysaccharide cell wall polymers (cellulose, hemicelluloses, pectins) are present. Altogether, 22, 32, and 6 CAZy families targeting pectin, hemicelluloses and cellulose polymers were identified (Supplementary Figure 7). Generally,

### A

<table><thead><tr><th>Phyla</th><th>Classes</th><th>Genus</th><th>Species (NCBI)</th></tr></thead><tbody><tr><td>Proteobacteria</td><td>Alphaproteobacteria *</td><td>Sphingomonas</td><td>Sphingomonas sp.</td></tr><tr><td>Proteobacteria</td><td>Gammaproteobacteria †*</td><td>Pseudomonas †*</td><td>Pseudomonas rhizosphaerae</td></tr><tr><td>Proteobacteria</td><td>Alphaproteobacteria *</td><td>Rhizobium</td><td>Rhizobium sp.</td></tr><tr><td>Proteobacteria</td><td>Betaproteobacteria *</td><td>Massilia</td><td>Massilia sp.</td></tr><tr><td>Proteobacteria</td><td>Gammaproteobacteria †*</td><td>Enterobacteriaceae (Family)</td><td>Pantoea vagans</td></tr><tr><td>Bacteroidetes</td><td>Sphingobacteria *</td><td>Pedobacter</td><td>Pedobacter sp.</td></tr><tr><td>Proteobacteria</td><td>Alphaproteobacteria *</td><td>Sphingomonas</td><td>Sphingomonas sp.</td></tr><tr><td>Proteobacteria</td><td>Gammaproteobacteria †*</td><td>Pseudomonas †*</td><td>Pseudomonas sp.</td></tr><tr><td>Bacteroidetes</td><td>Flavobacteria *</td><td>Epilithonimonas</td><td>Epilithonimonas sp.</td></tr><tr><td>Bacteroidetes</td><td>Flavobacteria *</td><td>Flavobacterium</td><td>Flavobacterium sp.</td></tr></tbody></table>

<table><thead><tr><th>OTUs (total count)</th></tr></thead><tbody><tr><td>OTU00001 (18312/22584)</td></tr><tr><td>OTU00002 (13903/15333)</td></tr><tr><td>OTU00003 (10603/10696)</td></tr><tr><td>OTU00004 (8978/11071)</td></tr><tr><td>OTU00006 (9294/6173)</td></tr><tr><td>OTU00005 (7096/7329)</td></tr><tr><td>OTU00007 (4120/6479)</td></tr><tr><td>OTU00008 (6557/3714)</td></tr><tr><td>OTU00009 (6115/3757)</td></tr><tr><td>OTU00010 (5983/2540)</td></tr></tbody></table>

Time point
Day of retting
Turning

Harvest (Maturity of fibers)

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

### B

<table><thead><tr><th>Phyla</th><th>Classes</th><th>Lowest level of taxonomic assignment by PIPITS (UNITE)</th><th>Species (NCBI)</th></tr></thead><tbody><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Capnodiales (Order) *</td><td>Cladosporium herbarum †</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Pleosporaceae (Family) *</td><td>Alternaria sp. †</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Fungi (Kingdom)</td><td>Epicoccum nigrum †</td></tr><tr><td>Basidiomycota</td><td>Tremellomycetes</td><td>Itersonilia perplexans (SH199072.07FU)</td><td>Itersonilia perplexans</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Capnodiales (Order) *</td><td>Cladosporium sp. †</td></tr><tr><td>Basidiomycota</td><td>Agaricostilbomycetes</td><td>Kondoa sp (SH476649.07FU)</td><td>Kondoa sp.</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Alternaria (Genus) *</td><td>Alternaria alternata †</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Pleosporales sp (SH215818.07FU)</td><td>Neoscochyta sp.</td></tr><tr><td>Ascomycota</td><td>Leotiomycetes *</td><td>Sclerotiniaceae (Family) *</td><td>Botrytis cinerea †</td></tr><tr><td>Ascomycota</td><td>Dothideomycetes *</td><td>Aureobasidium (Genus) *</td><td>Aureobasidium pullulans</td></tr></tbody></table>

<table><thead><tr><th>OTUs (total count)</th></tr></thead><tbody><tr><td>OTU2685 (133157/131260)</td></tr><tr><td>OTU1918 (39112/52323)</td></tr><tr><td>OTU2206 (26845/20359)</td></tr><tr><td>OTU448 (26450/18308)</td></tr><tr><td>OTU2906 (9114/5067)</td></tr><tr><td>OTU1100 (5504/8671)</td></tr><tr><td>OTU1949 (3406/6509)</td></tr><tr><td>OTU2268 (6886/2886)</td></tr><tr><td>OTU2877 (5679/2789)</td></tr><tr><td>OTU2280 (2759/5356)</td></tr></tbody></table>

Time point
Day of retting
Turning

Harvest (Maturity of fibers)

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

Type of samples

Plant

<caption>**FIGURE 5** | Heatmaps of the relative abundance of the top 10 OTUs of bacteria (A) and fungi (B) in plant samples (n = 14). OTUs previously associated with flax dew-retting† and water-retting*.</caption>
<figure><img src="image_10.png" /></figure>

the hydrolytic enzyme potential (all polymers) was greater during the first stages of retting (R0-R2/R3) compared to latter stages (R3/R4-R6) for both “early” and “standard” cultures. The drop in hydrolytic potential observed for R1 and R5 stages in “early” cultures is most likely related to the corresponding decrease in the most abundant bacterial OTU (e.g., Sphingomonas OTU00001, **Figure 5**). PICRUSt prediction does not exist for fungal OTUs and so hydrolytic enzyme potential cannot be directly predicted. Nevertheless, we were able to gain a relative idea of the overall hydrolytic enzyme potential by using the FUNGuild software that describes fungal trophic mode. Our results (**Figure 7**) show a progressive decrease in relative abundance of pathotrophs associated with a steady increase in saprotrophs and saprotrophs-pathotrophs as retting progresses. Pathogenic fungi generally produce a wider range of cell wall degrading enzymes than rot fungi and observed change in trophic mode during retting could suggest a decrease in hydrolytic enzyme diversity (Choi et al., 2013).

# DISCUSSION

## Microbial Identification and Retting Parameters

Previous studies using culture-based approaches and non-HTS metabarcoding have identified different bacteria and fungi phyla

present during retting including Actinobacteria, Firmicutes, Proteobacteria (bacteria), and Ascomycota and Zygomycota (fungi) (Lanigan, 1950; Rosemberg, 1965; Brown, 1984; Sharma, 1986a,b; Donaghy et al., 1990; Henriksson et al., 1997). Our results obtained using metabarcoding coupled with HTS not only identified these phyla, but also allowed the identification of new phyla not previously associated with dew-retting. Overall we identified 95 bacteria and 215 fungi species in dew-retted flax straw (plant) samples. HTS metabarcoding has been recently used to investigate bacterial (but not fungal) population dynamics in water-retted flax (Zhao et al., 2016). A comparison of relative abundances of the major bacterial phyla identified indicates that water retting is very different from dew-retting, despite the fact that the same lignocellulosic material is being degraded. Major phyla identified during water-retting were Firmicutes (genus Clostridium) and Proteobacteria (genera Azotobacter and Enterobacter). In contrast, Firmicutes were only present in low abundance during dew-retting and Azotobacter were absent. These differences can be most likely related to the anaerobic environment of water-retting compared to the more aerobic environment of dew-retting. Indeed, Clostridium is an obligate anaerobe and is known to be an agent of water-retting (Donaghy et al., 1990; Tamburini et al., 2003).

Phyla, identified in our study and not previously associated with flax dew-retting, included, for the bacteria,
<figure><img src="image_11.png" /><figcaption>**FIGURE 7** | Relative abundance of fungal OTUs classified by trophic mode. The trophic mode was assigned using the FUNGuild database and based on 35% of OTUs (i.e., total = 1048).</figcaption></figure>

Acidobacteria, Bacteroidetes, CKC4, Chlorobi, Fibrobacteres, Gemmatimonadetes, Nitrospirae and TM6; and for the fungi, Basidiomycota and Chytridiomycota. The Bacteroidetes phylum has been associated with cellulose degradation in agricultural soils (Schellenberger et al., 2010) and was previously detected in hemp dew-retting (Ribeiro et al., 2015) and flax water-retting (Zhao et al., 2016). Our observation of this phylum could indicate that it is also involved in flax dew-retting. Basidiomycota are linked to plant cell wall degradation in different ecosystems (Baldrian et al., 2008; Schneider et al., 2012; Kuramae et al., 2013; Voříšková and Baldrian, 2013; Rytioja et al., 2014) and were also detected in hemp dew-retting (Ribeiro et al., 2015).

Although the observation that new bacterial phyla (except for the Bacteroidetes) and fungal phyla represent less than 2% of the whole microbiota might suggest that they are not involved in the retting process, some of these phyla are related to microorganisms characterized as biomass degraders in previous studies (Zhao et al., 2014). This observation, together with the fact that low abundance OTUs can still contribute to the decomposition of plant matter (Baldrian et al., 2012) indicates that these phyla should not be ignored during the study of dew-retting.

A number of parameters potentially affecting microbial population structure during retting were examined. It is commonly admitted by farmers that the maturity of flax plants has a direct impact on the retting time and influences the choice for the pulling (up-rooting) date. Generally, straw from younger plants (flowering/green capsule stage) rets more quickly than that of more mature plants (yellow/brown capsule stage). This is thought to be related to differences in cell wall composition (e.g.,

pectin/lignin modifications and/or deposition) and water content (Meijer et al., 1995; Day et al., 2005; Akin, 2013). Our results showing that there was no significant difference in microbial communities and colonization dynamics between the early vs. standard cultures would suggest that differences in retting time may indeed be related to differences in cell wall structure and not to population differences.

Compared to litter decay that normally proceeds undisturbed, dew-retting is a semi-controlled process during which the straw swaths are turned by farmers to obtain a more uniform fiber separation. Our analyses revealed that this practice had a significant effect on both bacterial and fungal community membership and structure of the flax straw microbiome confirming a real microbiological effect of swath turning that probably contributes to a more uniform retting.

Although our results indicated no significant correlation between measured climatic conditions (temperature and rainfall) and community structures during the retting period it is important to remember that our study was conducted within a single year. It is possible that significant variations in community structures may occur between different seasons and further work is necessary to clarify this point.

## Microbial Dynamics

During dew-retting the relative abundance of the Bacteroidetes phylum increases while that of the Protobacteria decreases. A similar dynamic also occurs during biodegradation of field biomass from different angiosperm species (e.g., *Arundo donax*, *Eucalyptus camaldulensis*, and *Populus nigra*) suggesting, as might be expected, that similarities exist between the temporary
dew-retting ecosystem and degradation of lignocellulose in the field (Ventorino et al., 2015). Interestingly, the bacterial dynamics of flax dew-retting appear to be closer to that of field lignocellulose degradation than to that observed during flax water retting where Protobacteria increased during retting (Zhao et al., 2016). In this latter case, the phylum Proteobacteria was mainly represented by the genera Azotobacter that increased during retting and (to a much lesser extent) Enterobacter that remained constant. For fungal phyla we observed an increase in the relative abundance of Ascomycota at the expense of Basidiomycota in contrast to the situation generally observed during both field and forest litter decomposition (Schneider et al., 2012; Kuramae et al., 2013; Voříšková and Baldrian, 2013). The observed increase of Ascomycota was due to the saprophytic Altenaria species (Dang et al., 2015) that has previously been linked to later stages of dew-retting (Brown et al., 1986). In contrast, Altenaria species are more abundant during initial stages of litter decay (Snajdr et al., 2011). Our results also indicated that *C. herbarum* and *Epicoccum nigrum* contributed to the increase in Ascomycota relative abundance. During this stage less recalcitrant components of the biomass (pectins, and hemicelluloses) are progressively degraded (Dilly et al., 2001). Contrary to litter decay, dew-retting is a semi-controlled process and the challenge is to limit degradation of major quality related polymers such as crystalline cellulose. In this context, changes in the relative abundance of Ascomycota vs. Basidiomycota could represent an interesting bioindicator of retting progress.

More detailed information on population dynamics at different time points during retting was provided by analyzing the relative abundance of OTUs at different taxonomic rank (e.g., phyla, classes, or genus/species level). The most abundant bacterial OTU corresponded to Sphingomonas sp. that was present throughout most of the retting period in both early and standard cultures. Although Sphingomonas species have been previously identified during bamboo and hemp retting, as well as in forest litter microbiome, this is the first time they have been found in flax retting (Fu et al., 2011; Ribeiro et al., 2015) (Urbanová et al., 2015). These species are able to hydrolyze terminal non-reducing alpha-L-rhamnose residues in alpha-L-rhamnosides giving them the ability to degrade pectin (rhamnogalacturonan I and rhamnogalacturonan II) in the middle lamella (Hashimoto and Murata, 1998). Another *Sphingomonas* species, *S. paucimobilis* is also able to degrade lignin (Masai et al., 1999; de Gonzalo et al., 2016). The second most abundant OTU corresponded to *P. rhizosphaerae*, present during the early and medium retting stages but decreasing in latter stages. A number of Pseudomonas species have previously been associated with retting of different fiber plants (e.g., flax, hemp, jute, ramie) (Rosemberg, 1965; Munshi and Chattoo, 2008; Duan et al., 2012; Ribeiro et al., 2015). Pseudomonas sp. is considered as one of the most efficient lignin degradation bacterium (Shui Yang et al., 2007) and the genomes of both *Pseudomonas putida* and *Pseudomonas aeruginosa* contain genes encoding endoglucanases (Talia et al., 2012). Other abundant OTUs corresponded to Rhizobium, Pedobacter, and Flavobacterium that are known to show pectinase, cellulose, and hemicellulose activities (Mateos et al., 1992; McBride et al., 2009;

López-Mondéjar et al., 2016). In addition, Pedobacter has also been identified during bamboo and hemp retting (Fu et al., 2011; Ribeiro et al., 2015) or forest litter degradation (Urbanová et al., 2015). In contrast to Sphingomonas and Pseudomonas, these organisms become more abundant toward the end of the retting period and could be associated with "over-retting" when the structural integrity of the fiber starts to be degraded.

In contrast to the more evenly distributed abundance of the bacterial OTUs, fungal OTUs were dominated by one major species—*C. herbarum*—that rapidly increased during early retting. This species, as well as the third most abundant OTU (*E. nigrum*) are known to be common dew-retting agents and are believed to degrade cellulose (Brown, 1984). Of the other fungal OTUs, all have previously been associated with dew-/water-retting except for *Itersonilia perplexans*. Interestingly, our results also indicated that *Alternaria alternata* is present at the start of retting. Traditionally, the appearance of this species is used as a signal that retting is starting to go too far and that the swaths should be collected (Brown et al., 1986).

## Hydrolytic Enzyme Potential

Prediction of hydrolytic enzymes potentially present during retting was performed by using PICRUST (Langille et al., 2013). This software successfully predicts bacterial enzymatic activities represented in different databases (e.g., KEGG Ortholog, COGs, or CAZy). Overall, a large collection of enzyme activities targeting both the main backbones and side chains of the major polysaccharide polymers were identified. Based on OTU counts, ~38, 43, and 19 percent of the total hydrolytic enzyme potential targeted pectins, hemicelluloses, and cellulose, respectively. Despite the clear dynamics and significant changes in the straw microbiome these values remained constant throughout the retting period. Similar software does not exist for predicting fungal enzyme potential. This represents an important hurdle for obtaining a complete overview of the dew-retting process as fungi are major producers of extracellular hydrolytic enzymes (Schneider et al., 2012). Nevertheless, FUNGuild analysis showed that pathogenic taxa, present at the beginning of retting are progressively replaced by saprophytic fungi, more able to degrade lignocellulose. This change is most likely related to the fact that flax plants are still living when up-rooted.

In conclusion, we have shown that HTS metabarcoding is a powerful technique for analyzing complex bacterial and fungal community dynamics during flax dew-retting that can be used to identify different factors affecting the microbiota and—potentially—fiber isolation and quality. However, these results were obtained on samples retted in 1 year and it will be necessary to validate these data over several seasons. The use of PICRUST data allows a predictive study of potential bacterial hydrolytic activity but should be coupled in future studies with alternative meta-omics methods such as metatranscriptomic or metaproteomic coupled with metagenomics to facilitate the assembling with appropriate reference genomes (Schneider et al., 2012; Dai et al., 2015; Hesse et al., 2015; Kuske et al., 2015; Wu et al., 2015). Such an approach would not only allow confirmation of bacterial enzyme dynamics but would also enable identification of fungal enzymes involved in this process.
# AVAILABILITY OF DATA AND MATERIALS

The microbial DNA sequencing data sets supporting the results in this article are available at the EBI ENA with accession number PRJEB20299.

# AUTHOR CONTRIBUTIONS

Conceptualization: CD, SG, and SH; Methodology: CD, SG, and SH; Experimentation: CD; Bioinformatic and statistical analysis: CD; Writing—original draft: CD; Writing—Review and editing: SG and SH; Funding Acquisition: SG and SH.

# FUNDING

This work was funded within the framework of the collaborative French "Future project" SINFONI. CD thanks the region of Hauts-de-France and Bpifrance for their financial support.

# ACKNOWLEDGMENTS

The authors would like to thank the following people/organizations: The flax farmer C.A.L.I.R.A (Coopérative Agricole Linière de la Région d'Abbeville) for growing and retting flax plants, in particular Vincent DELAPORTE for his valuable assistance. The Genomic and Transcriptomic platform of Genopole Occitanie-Toulouse (INRA, GeT-Plage, http://get.genotoul.fr/) where sequencing was performed and more particularly Catherine ZHANCHETTA and Olivier BOUCHEZ. The Plateau d'Ecologie moléculaire et biochimie Evo-Eco-Paléo (CNRS-UMR 8198, Evolution, Ecologie et Paléontologie) for use of the LightCycler 480 and epMotion and more particularly Cécile GODÉ and Anne-Catherine HOLL for their technical support. Antoine PORTELETTE, Julien LE

ROY, Sandrine ARRIBAT, and Brigitte CHABBERT for their valuable assistance during the sampling. Alexandrine THORE for her invaluable help in the Supplementary Figure 2 graphical design.

# SUPPLEMENTARY MATERIAL

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

**Supplementary Figure 1** | Rainfall, daily amplitudes of temperature and moisture, and retting time points at Martainneville during the 2014 dew-retting period. Raw data from the Abbeville station available on the infoclimat website (https://www.infoclimat.fr).

**Supplementary Figure 2** | Schematic representation of the experimental site at Martainneville during the 2014 dew-retting campaign.

**Supplementary Figure 3** | Rarefaction curves of observed OTUs in bacterial (16S) soil (A) and stem (B) samples, defined at a 97% sequence similarity cut-off.

**Supplementary Figure 4** | Rarefaction curves of observed OTUs in fungal (ITS) soil (A) and plant (B) samples, defined at a 97% sequence similarity cut-off.

**Supplementary Figure 5** | Alpha diversity of bacterial (A–C) and fungal (D–F) communities in soil and plant samples during flax dew-retting. (A,D) Richness (Chao1), (B,E) Evenness (Heip's), and (C,F) Diversity (Inverse Simpson). Different colors indicate the source samples (light brown, soil early harvest; dark brown, soil standard harvest; light green, plant early harvest; dark green, plant standard harvest). Plots show the distribution of results of 20,548 bacterial and 42,436 fungal sequences subsampled from each sample 1,000 times and calculated for the average (error bars for Chao1 and Inverse Simpson correspond to lower and higher bound 95% confidence intervals and for Heip's represent standard deviation).

**Supplementary Figure 6** | Before- (A,C) and after- (B,D) swath turning biomarkers for bacteria (A,B) and fungi (C,D). Biomarkers were identified by analyzing differential OTUs abundance using Metastats, Indicator and LEfSe.

**Supplementary Figure 7** | Heatmaps showing the total counts of all CAZyme families in early and standard harvests obtained from the OTUs table (Greengenes Database used for consensus taxonomy) and generated by PICRUSt software.

# REFERENCES

Abarenkov, K., Nilsson, R. H., Larsson, K. H., Alexander, I. J., Eberhardt, U., Erland, S., et al. (2010). The UNITE database for molecular identification of fungi - recent updates and future perspectives. *New Phytol.* 186, 281–285. doi: 10.1111/j.1469-8137.2009.03160.x

Akin, D. E. (2008). Plant cell wall aromatics: influence on degradation of biomass. *Biofuels Bioprod.* *Biorefining* 2, 288–303. doi: 10.1002/bbb.76

Akin, D. E. (2013). Linen most useful: perspectives on structure, chemistry, and enzymes for retting flax. *ISRN Biotechnol.* 2013, 23. doi: 10.5402/2013/186534

Akin, D. E., Rigsby, L. L., Henriksson, G., and Eriksson, K.-E. L. (1998). Structural effects on flax stems of three potential retting fungi. *Text. Res. J.* 68, 515–519. doi: 10.1177/004051759806800708

Amann, R. I., Ludwig, W., and Schleifer, K.-H. (1995). Phylogenetic identification and *in situ* detection of individual microbial cells without cultivation. *Microbiol. Rev.* 59, 143–169.

Baldrian, P., Kolařík, M., Stursová, M., Kopecký, J., Valášková, V., Vítrovský, T., et al. (2012). Active and total microbial communities in forest soil are largely different and highly stratified during decomposition. *ISME J.* 6, 248–258. doi: 10.1038/ismej.2011.95

Baldrian, P., Valášková, V., and Valášková, V. (2008). Degradation of cellulose by basidiomycetous fungi. *FEMS Microbiol. Rev.* 32, 501–521. doi: 10.1111/j.1574-6976.2008.00106.x

Bardou, P., Mariette, J., Escudié, F., Djemiel, C., and Klopp, C. (2014). jvenn: an interactive Venn diagram viewer. *BMC Bioinformatics* 15:293. doi: 10.1186/1471-2105-15-293

Barwell, L. J., Isaac, N. J., and Kunin, W. E. (2015). Measuring β -diversity with species abundance data. *J. Anim. Ecol.* 84, 1112–1122. doi: 10.1111/1365-2656.12362

Beckers, B., Op De Beeck, M., Weyens, N., Boerjan, W., and Vangronsveld, J. (2017). Structural variability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of field-grown poplar trees. *Microbiome* 5, 25. doi: 10.1186/s40168-017-0241-2

Beer, C., Reichstein, M., Tomelleri, E., Ciais, P., Jung, M., Carvalhais, N., et al. (2010). Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. *Science* 329, 834–838. doi: 10.1126/science.1184984

Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., et al. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. *Methods Ecol. Evol.* 4, 914–919. doi: 10.1111/2041-210X.12073

Bodenhausen, N., Horton, M. W., and Bergelson, J. (2013). Bacterial communities associated with the leaves and the roots of *Arabidopsis thaliana*. *PLoS ONE* 8:e56329. doi: 10.1371/journal.pone.0056329

Boerjan, W., Ralph, J., and Baucher, M. (2003). Lignin biosynthesis. *Annu. Rev. Plant Biol.* 54, 519–546. doi: 10.1146/annurev.arplant.54.031902.134938
Bokulich, N. A., and Mills, D. A. (2013). Improved selection of internal transcribed spacer-specific primers enables quantitative, ultra-high-throughput profiling of fungal communities. *Appl. Environ. Microbiol.* 79, 2519–2526. doi: 10.1128/AEM.03870-12

Brown, A. E. (1984). *Epicoccum nigrum*, a primary saprophyte involved in the retting of flax. *Trans. Br. Mycol. Soc.* 83, 29–35. doi: 10.1016/S0007-1536(84)80242-9

Brown, A. E., and Sharma, H. S. S. (1984). Production of polysaccharide-degrading enzymes by saprophytic fungi from glyphosate-treated flax and their involvement in retting. *Ann. Appl. Biol.* 105, 65–74. doi: 10.1111/j.1744-7348.1984.tb02803.x

Brown, A. E., Sharma, H. S. S., and Black, D. L. R. (1986). Relationship between pectin content of stems of flax cultivars, fungal cell wall-degrading enzymes and pre-harvest retting. *Ann. Appl. Biol.* 109, 345–351. doi: 10.1111/j.1744-7348.1986.tb05326.x

Campilho, R. D. S. G. (2015). *Natural Fiber Composites*. Boca Raton, FL: CRC Press.

Cardenas, E., Kranabetter, J. M., Hope, G., Maas, K. R., Hallam, S., and Mohn, W. W. (2015). Forest harvesting reduces the soil metagenomic potential for biomass decomposition. *ISME J.* 9, 1–12. doi: 10.1038/ismej.2015.57

Chao, A. (1984). Nonparametric estimation of the number of classes in a population. *Scand. J. Stat.* 11, 265–270. doi: 10.1214/aoms/1177729949

Choi, J., Kim, K.-T., Jeon, J., and Lee, Y.-H. (2013). Fungal plant cell wall-degrading enzyme database: a platform for comparative and evolutionary genomics in fungi and Oomycetes. *BMC Genomics* 14(Suppl. 5):S7. doi: 10.1186/1471-2164-14-S5-S7

Cragg, S. M., Beckham, G. T., Bruce, N. C., Bugg, T. D., Distel, D. L., Dupree, P., et al. (2015). Lignocellulose degradation mechanisms across the Tree of Life. *Curr. Opin. Chem. Biol.* 29, 108–119. doi: 10.1016/j.cbpa.2015.10.018

Dai, X., Tian, Y., Li, J., Luo, Y., Liu, D., Zheng, H., et al. (2015). Metatranscriptomic analyses of plant cell wall polysaccharide degradation by microorganisms in the cow rumen. *Appl. Environ. Microbiol.* 81, 1375–1386. doi: 10.1128/AEM.03682-14

Dang, H. X., Pryor, B., Peever, T., and Lawrence, C. B. (2015). The Alternaria genomes database: a comprehensive resource for a fungal genus comprised of saprophytes, plant pathogens, and allergenic species. *BMC Genomics* 16:239. doi: 10.1186/s12864-015-1430-7

Day, A., Ruel, K., Neutelings, G., Crónier, D., David, H., Hawkins, S., et al. (2005). Lignification in the flax stem: evidence for an unusual lignin in bast fibers. *Planta* 222, 234–245. doi: 10.1007/s00425-005-1537-1

de Gonzalo, G., Colpa, D. I., Habib, M. H., and Fraaije, M. W. (2016). Bacterial enzymes involved in lignin degradation. *J. Biotechnol.* 236, 110–119. doi: 10.1016/j.jbiotec.2016.08.011

Dilly, O., Bartsch, S., Rosenbrock, P., Buscot, F., and Munch, J. C. (2001). Shifts in physiological capabilities of the microbiota during the decomposition of leaf litter in a black alder (*Alnus glutinosa* (Gaertn.) L.) forest. *Soil Biol. Biochem.* 33, 921–930. doi: 10.1016/S0038-0717(00)00239-X

Donaghy, J. A., Levett, P. N., and Haylock, R. W. (1990). Changes in microbial populations during anaerobic flax retting. *J. Appl. Bacteriol.* 69, 634–641. doi: 10.1111/j.1365-2672.1990.tb01556.x

Duan, S., Liu, Z., Feng, X., and Zheng, K. (2012). Diversity and characterization of ramie-degumming strains. *Sci. Agric.* 69, 119–125. doi: 10.1590/S0103-90162012000200006

Eastwood, D. C., Floudas, D., Binder, M., Majcherczyk, A., Schneider, P., Aerts, A., et al. (2011). The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. *Science* 333, 762–765. doi: 10.1126/science.1205411

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

Esty, W. W. (1986). The efficiency of good's nonparametric coverage estimator. *Ann. Stat.* 14, 1257–1260. doi: 10.1214/aos/1176350066

Excoffier, L., Smouse, P. E., and Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. *Genetics* 131, 479–491. doi: 10.1007/s00424-009-0730-7

Fu, J., Mueller, H., de Castro, J. V., Yu, C., Cavaco-Paulo, A., Guebitz, G. M., et al. (2011). Changes in the bacterial community structure and diversity during bamboo retting. *Biotechnol. J.* 6, 1262–1271. doi: 10.1002/biot.201100105

Gołebiewski, M., Deja-Sikora, E., Cichosz, M., Tretyn, A., and Wróbel, B. (2014). 16S rDNA pyrosequencing analysis of bacterial community in heavy metals polluted soils. *Microb. Ecol.* 67, 635–647. doi: 10.1007/s00248-013-0344-7

Good, I. J. (1953). The population frequencies of species and the estimation of population parameters. *Biometrika* 40, 237–264. doi: 10.1093/biomet/40.3.4-237

Gübitz, G. M., and Cavaco-Paulo, A. (2001). Biotechnology in the textile industry–perspectives for the new millennium. *J. Biotechnol.* 89, 89–90. doi: 10.1016/S0168-1656(01)00317-0

Gweon, H. S., Oliver, A., Taylor, J., Booth, T., Gibbs, M., Read, D. S., et al. (2015). PIPITS: an automated pipeline for analyses of fungal internal transcribed spacer sequences from the Illumina sequencing platform. *Methods Ecol. Evol.* 6, 973–980. doi: 10.1111/2041-210X.12399

Hashimoto, W., and Murata, K. (1998). alpha-L-rhamnosidase of Sphingomonas sp. R1 producing an unusual exopolysaccharide of sphingan. *Biosci. Biotechnol. Biochem.* 62, 1068–1074. doi: 10.1271/bbb.62.1068

Heip, C. (1974). A new index measuring evenness. *J. Mar. Biol. Assoc. U. K.* 54, 555–557. doi: 10.1017/S0025315400022736

Henriksson, G., Akin, D. E., Rigsby, L. L., Patel, N., and Eriksson, K.-E. L. (1997). Influence of chelating agents and mechanical pretreatment on enzymatic retting of flax. *Text. Res. J.* 67, 829–836. doi: 10.1177/004051759706 701107

Henriksson, G., Akin, D. E., Slomczynski, D., and Eriksson, K. E. L. (1999). Production of highly efficient enzymes for flax retting by *Rhizomucor pusillus*. *J. Biotechnol.* 68, 115–123. doi: 10.1016/S0168-1656(98)00192-8

Hesse, C. N., Mueller, R. C., Vuyisich, M., Gallegos-Graves, L. V., Gleasner, C. D., Zak, D. R., et al. (2015). Forest floor community metatranscriptomes identify fungal and bacterial responses to N deposition in two maple forests. *Front. Microbiol.* 6:337. doi: 10.3389/fmicb.2015.00337

Hirsch, P. R., Mauchline, T. H., and Clark, I. M. (2010). Culture-independent molecular techniques for soil microbial ecology. *Soil Biol. Biochem.* 42, 878–887. doi: 10.1016/j.soilbio.2010.02.019

Huse, S. M., Welch, D. M., Morrison, H. G., and Sogin, M. L. (2010). Ironing out the wrinkles in the rare biosphere through improved OTU clustering. *Environ. Microbiol.* 12, 1889–1898. doi: 10.1111/j.1462-2920.2010.02193.x

Ihrmark, K., Bödeker, I. T., Cruz-Martinez, K., Friberg, H., Kubartova, A., Schenck, J., et al. (2012). New primers to amplify the fungal ITS2 region–evaluation by 454-sequencing of artificial and natural communities. *FEMS Microbiol. Ecol.* 82, 666–677. doi: 10.1111/j.1574-6941.2012.01437.x

Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. *Nucleic Acids Res.* 41, e1. doi: 10.1093/nar/gks808

Knief, C. (2014). Analysis of plant microbe interactions in the era of next generation sequencing technologies. *Front. Plant Sci.* 5:216. doi: 10.3389/fplrs.2014.00216

Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. *Appl. Environ. Microbiol.* 79, 5112–5120. doi: 10.1128/AEM.01043-13

Kubicek, C. P., Starr, T. L., and Glass, N. L. (2014). Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. *Annu. Rev. Phytopathol.* 52, 427–451. doi: 10.1146/annurev-phyto-102313-045831

Kuhad, R. C., and Singh, A. (1993). Lignocellulose biotechnology: current and future prospects. *Crit. Rev. Biotechnol.* 13, 151–172. doi: 10.3109/07388559309040630

Kuramae, E. E., Hillekens, R. H., De Hollander, M., Van Der Heijden, M. G., Van Den Berg, M., Van Straalen, N. M., et al. (2013). Structural and functional variation in soil fungal communities associated with lit ter bags containing maize leaf. *FEMS Microbiol. Ecol.* 84, 519–531. doi: 10.1111/1574-6941.12080

Kuske, C. R., Hesse, C. N., Challacombe, J. F., Cullen, D., Herr, J. R., Mueller, R. C., et al. (2015). Prospects and challenges for fungal metatranscriptomics of complex communities. *Fungal Ecol.* 14, 133–137. doi: 10.1016/j.funeco.2014.12.005

Kvavadze, E., Bar-Yosef, O., Belfer-Cohen, A., Boaretto, E., Jakeli, N., Matskevich, Z., et al. (2009). 30,000-year-old wild flax fibers. *Science* 325, 1359–1359. doi: 10.1126/science.1175404
Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat. Biotechnol.* 31, 814–821. doi: 10.1038/nbt.2676

Lanigan, G., w. (1950). Relationship of fiber quality of bacterial flora in flax retting. *Nature* 165, 516–517. doi: 10.1038/165516a0

López-Mondéjar, R., Zühlke, D., Becher, D., Riedel, K., and Baldrian, P. (2016). Cellulose and hemicellulose decomposition by forest soil bacteria proceeds by the action of structurally variable enzymatic systems. *Sci. Rep.* 6:25279. doi: 10.1038/srep25279

Ludwig, J. A., and Reynolds, J. F. (1988). *Statistical Ecology: a Primer on Methods and Computing*. Available online at: http://dl.acm.org/citation.cfm?id=61827 (Accessed May 8, 2016).

Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core *Arabidopsis thaliana* root microbiome. *Nature* 488, 86–90. doi: 10.1038/nature11237

Lynd, L. R., Weimer, P. J., van Zyl, W. H., and Pretorius, I. S. (2002). Microbial cellulose utilization: fundamentals and biotechnology. *Microbiol. Mol. Biol. Rev.* 66, 506–577. doi: 10.1128/MMBR.66.3.506-577.2002

Masai, E., Katayama, Y., Nishikawa, S., and Fukuda, M. (1999). Characterization of *Sphingomonas paucimobilis* SYK-6 genes involved in degradation of lignin-related compounds. *J. Ind. Microbiol. Biotechnol.* 23, 364–373. doi: 10.1038/sj.jim.2900747

Mateos, P. F., Jimenez-Zurdo, J. I., Chen, J., Squartini, A. S., Haack, S. K., Martinez-Molina, E., et al. (1992). Cell-associated pectinolytic and cellulolytic enzymes in *Rhizobium leguminosarum* biovar trifolii. *Appl. Environ. Microbiol.* 58, 1816–1822.

McBride, M. J., Xie, G., Martens, E. C., Lapidus, A., Henrissat, B., Rhodes, R. G., et al. (2009). Novel features of the polysaccharide-digesting gliding bacterium *Flavobacterium johnsoniae* as revealed by genome sequence analysis. *Appl. Environ. Microbiol.* 75, 6864–6875. doi: 10.1128/AEM.01495-09

Md. Tahir, P., Ahmed, A. B., SaifulAzry, S. O. A., and Ahmed, Z. (2011). Retting process of some bast plant fibres and its effect on fibre quality: a review. *BioResources* 6, 5260–5281. doi: 10.15376/biores.6.4.5260-5281

Meijer, W. J. M., Vertregt, N., Rutgers, B., and van de Waart, M. (1995). The pectin content as a measure of the retting and rettability of flax. *Ind. Crops Prod.* 4, 273–284. doi: 10.1016/0926-6690(95)00041-0

Montella, S., Ventorino, V., Lombard, V., Henrissat, B., Pepe, O., and Faraco, V. (2017). Discovery of genes coding for carbohydrate-active enzyme by metagenomic analysis of lignocellulosic biomasses. *Sci. Rep.* 7:42623. doi: 10.1038/srep42623

Munshi, T. K., and Chattoo, B. B. (2008). Bacterial population structure of the jute-retting environment. *Microb. Ecol.* 56, 270–282. doi: 10.1007/s00248-007-9345-8

Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. *Fungal Ecol.* 20, 241–248. doi: 10.1016/j.funeco.2015.06.006

Peršoh, D. (2015). Plant-associated fungal communities in the light of meta'omics. *Fungal Divers.* 75, 1–25. doi: 10.1007/s13225-015-0334-9

Pil, L., Bensadoun, F., Pariset, J., and Verpoest, I. (2016). Why are designers fascinated by flax and hemp fibre composites? *Compos. Part A* *Appl. Sci. Manuf.* 83, 193–205. doi: 10.1016/j.compositesa.2015.11.004

Plassart, P., Terrat, S., Thomson, B., Griffiths, R., Dequiedt, S., Lelievre, M., et al. (2012). Evaluation of the ISO standard 11063 DNA extraction procedure for assessing soil microbial abundance and community structure. *PLoS ONE* 7:e44279. doi: 10.1371/journal.pone.0044279

Preisner, M., Wojtasik, W., Kulma, A., Zuk, M., and Szopa, J. (2014). "Flax Fiber," in *Kirk-Othmer Encyclopedia of Chemical Technology* (Hoboken, NJ: John Wiley & Sons, Inc.), 1–32. doi: 10.1002/0471238961.0612012401110914.a01.pub2

Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. *Nucleic Acids Res.* 35, 7188–7196. doi: 10.1093/nar/gkm864

Ribeiro, A., Pochart, P., Day, A., Mennuni, S., Bono, P., Baret, J.-L., et al. (2015). Microbial diversity observed during hemp retting. *Appl. Microbiol. Biotechnol.* 99, 4471–4484. doi: 10.1007/s00253-014-6356-5

Rosemberg, J. A. (1965). Bacteria responsible for the retting of Brazilian flax. *Appl. Microbiol.* 13, 991–992.

Rytioja, J., Hildén, K., Yuzon, J., Hatakka, A., de Vries, R. P., and Mäkelä, M. R. (2014). Plant-polysaccharide-degrading enzymes from basidiomycetes. *Microbiol. Mol. Biol. Rev.* 78, 614–649. doi: 10.1128/MMBR.00035-14

Schellenberger, S., Kolb, S., and Drake, H. L. (2010). Metabolic responses of novel cellulolytic and saccharolytic agricultural soil Bacteria to oxygen. *Environ. Microbiol.* 12, 845–861. doi: 10.1111/j.1462-2920.2009.02128.x

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microbiol.* 75, 7537–7541. doi: 10.1128/AEM.01541-09

Schneider, T., Keiblinger, K. M., Schmid, E., Sterflinger-Gleixner, K., Ellersdorfer, G., Roschitzki, B., et al. (2012). Who is who in litter decomposition? Metaproteomics reveals major microbial players and their biogeochemical functions. *ISME J.* 6, 1749–1762. doi: 10.1038/ismej.2012.11

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. *Genome Biol.* 12:R60. doi: 10.1186/gb-2011-12-6-r60

Sharma, H. (1986a). Enzymatic degradation of residual non-cellulosic polysaccharides present on dew-retted flax fibres. *Appl. Microbiol. Biotechnol.* 197, 358–362.

Sharma, H. (1986b). The role of bacteria in retting of desiccated flax during damp weather. *Appl. Microbiol. Biotechnol.* 24, 463–467.

Shui Yang, J., Ren Ni, J., Li Yuan, H., and Wang, E. (2007). Biodegradation of three different wood chips by Pseudomonas sp. PKE117. *Int. Biodeterior. Biodegrad.* 60, 90–95. doi: 10.1016/j.ibiod.2006.12.006

Simpson, E. H. (1949). Measurement of diversity. *Nature* 163, 688–688. doi: 10.1038/163688a0

Snajdr, J., Cajthaml, T., Valášková, V., Merhautová, V., Petránková, M., Spetz, P., et al. (2011). Transformation of *Quercus petraea* litter: successive changes in litter chemistry are reflected in differential enzyme activity and changes in the microbial community composition. *FEMS Microbiol. Ecol.* 75, 291–303. doi: 10.1111/j.1574-6941.2010.00999.x

Soliveres, S., van der Plas, F., Manning, P., Prati, D., Gossner, M. M., Renner, S. C., et al. (2016). Biodiversity at multiple trophic levels is needed for ecosystem multifunctionality. *Nature* 536, 456–459. doi: 10.1038/nature19092

Staley, J. T., and Konopka, A. (1985). Measurement of *in situ* activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. *Annu. Rev. Microbiol.* 39, 321–346. doi: 10.1146/annurev.mi.39.100185.001541

Suenaga, H. (2012). Targeted metagenomics: a high-resolution metagenomics approach for specific gene clusters in complex microbial communities. *Environ. Microbiol.* 14, 13–22. doi: 10.1111/j.1462-2920.2011.02438.x

Talia, P., Sede, S. M., Campos, E., Rorig, M., Principi, D., Tosto, D., et al. (2012). Biodiversity characterization of cellulolytic bacteria present on native Chaco soil by comparison of ribosomal RNA genes. *Res. Microbiol.* 163, 221–232. doi: 10.1016/j.resmic.2011.12.001

Tamburini, E., León, A. G., Perito, B., and Mastromei, G. (2003). Characterization of bacterial pectinolytic strains involved in the water retting process. *Environ. Microbiol.* 5, 730–736. doi: 10.1046/j.1462-2920.2003.00462.x

Terrat, S., Christen, R., Dequiedt, S., Lelièvre, M., Nowak, V., Regnier, T., et al. (2012). Molecular biomass and MetaTaxogenomic assessment of soil microbial communities as influenced by soil DNA extraction procedure. *Microb. Biotechnol.* 5, 135–141. doi: 10.1111/j.1751-7915.2011.00307.x

Terrat, S., Plassart, P., Bourgeois, E., Ferreira, S., Dequiedt, S., Adele-Dit-DeRenseville, N., et al. (2015). Meta-barcoded evaluation of the ISO standard 11063 DNA extraction procedure to characterize soil bacterial and fungal community diversity and composition. *Microb. Biotechnol.* 8, 131–142. doi: 10.1111/1751-7915.12162

Toju, H., Tanabe, A. S., Yamamoto, S., and Sato, H. (2012). High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. *PLoS ONE* 7:e40863. doi: 10.1371/journal.pone.0040863

Urbanová, M., Šnajdr, J., and Baldrian, P. (2015). Composition of fungal and bacterial communities in forest litter and soil is largely determined by dominant trees. *Soil Biol. Biochem.* 84, 53–64. doi: 10.1016/j.soilbio.2015.02.011
Ventorino, V., Aliberti, A., Faraco, V., Robertiello, A., Giacobbe, S., Ercolini, D., et al. (2015). Exploring the microbiota dynamics related to vegetable biomasses degradation and study of lignocellulose-degrading bacteria for industrial biotechnological application. *Sci. Rep.* 5:8161. doi: 10.1038/srep08161

Visi, D. K., D'Souza, N., Ayre, B. G., Webber Iii, C. L., and Allen, M. S. (2013). Investigation of the bacterial retting community of kenaf (*Hibiscus cannabinus*) under different conditions using next-generation semiconductor sequencing. *J. Ind. Microbiol. Biotechnol.* 40, 465–475. doi: 10.1007/s10295-013-1242-1

Voříšková, J., and Baldrian, P. (2013). Fungal community on decomposing leaf litter undergoes rapid successional changes. *ISME J.* 7, 477–486. doi: 10.1038/ismej.2012.116

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. *Appl. Environ. Microbiol.* 73, 5261–5267. doi: 10.1128/AEM.00062-07

Warren, R. A. (1996). Microbial hydrolysis of polysaccharides. *Annu. Rev. Microbiol.* 50, 183–212.

White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. *PLoS Comput. Biol.* 5:e1000352. doi: 10.1371/journal.pcbi.1000352

Wu, S., Ren, Y., Peng, C., Hao, Y., Xiong, F., Wang, G., et al. (2015). Metatranscriptomic discovery of plant biomass-degrading capacity from grass carp intestinal microbiomes. *FEMS Microbiol. Ecol.* 91:fiv107. doi: 10.1093/femsec/fiv107

Yue, J. C., and Clayton, M. K. (2005). A similarity measure based on species proportions. *Commun. Stat. Theor. Methods* 34, 2123–2131. doi: 10.1080/STA-200066418

Zhang, J., Henriksson, H., Szabo, I. J., Henriksson, G., and Johansson, G. (2005). The active component in the flax-retting system of the zygomycete *Rhizopus oryzae* sb is a family 28 polygalacturonase. *J. Ind. Microbiol. Biotechnol.* 32, 431–438. doi: 10.1007/s10295-005-0014-y

Zhao, D., Liu, P., Pan, C., Du, R., Ping, W., and Ge, J. (2016). Bacterial succession and metabolite changes during flax (*Linum usitatissimum* L.) retting with *Bacillus cereus* HDYM-02. *Sci. Rep.* 6:31812. doi: 10.1038/srep31812

Zhao, Z., Liu, H., Wang, C., and Xu, J. (2014). Correction: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. *BMC Genomics* 15:6. doi: 10.1186/1471-2164-15-6

Zhou, C.-H., Xia, X., Lin, C.-X., Tong, D.-S., and Beltramini, J. (2011). Catalytic conversion of lignocellulosic biomass to fine chemicals and fuels. *Chem. Soc. Rev.* 40, 5588–5617. doi: 10.1039/c1cs15124j

**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 © 2017 Djemiel, Grec and Hawkins. 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) or licensor 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.