Targeted Metagenomics of Retting in Flax: The Beginning of the Quest to Harness the Secret Powers of the Microbiota
- 1Univ. Lille, CNRS, UMR 8576 – UGSF – Unité de Glycobiologie Structurale et Fonctionnelle, Lille, France
- 2Université de Reims Champagne Ardenne, INRAE, UMR FARE A 614, Reims, France
The mechanical and chemical properties of natural plant fibers are determined by many different factors, both intrinsic and extrinsic to the plant, during growth but also after harvest. A better understanding of how all these factors exert their effect and how they interact is necessary to be able to optimize fiber quality for use in different industries. One important factor is the post-harvest process known as retting, representing the first step in the extraction of bast fibers from the stem of species such as flax and hemp. During this process microorganisms colonize the stem and produce hydrolytic enzymes that target cell wall polymers thereby facilitating the progressive destruction of the stem and fiber bundles. Recent advances in sequencing technology have allowed researchers to implement targeted metagenomics leading to a much better characterization of the microbial communities involved in retting, as well as an improved understanding of microbial dynamics. In this paper we review how our current knowledge of the microbiology of retting has been improved by targeted metagenomics and discuss how related ‘-omics’ approaches might be used to fully characterize the functional capability of the retting microbiome.
Natural fibers from different plant species have long been used by man to make textiles and are now being increasingly exploited as a viable replacement for synthetic fibers in composite materials. From a biological point of view the fibers used in textiles and composites are single cells, much longer than they are wide and characterized by the presence of a thick secondary cell wall. Such structures can be found in the xylem (wood fibers) and/or associated with the phloem in the outer tissues of the stem of non-woody plants (e.g., the bast fibers of flax, hemp, ramie), as well as in the walls of the fruits or the leaves of some species (e.g., cotton, kapok, sisal). In bast fiber species, the individual fiber cells (elementary fibers) are grouped together to form so-called fiber bundles.
The ‘quality’ of the fabricated textiles and natural fiber composites (NFCs) is determined by both the mechanical and chemical properties of the isolated fibers/fiber bundles and the subsequent industrial transformation process (Müssig, 2010; Summerscales et al., 2010; Bourmaud et al., 2018), and researchers have therefore explored the different factors influencing the quality of the finished products.
It is generally accepted that the mechanical and chemical properties of plant fibers result from fiber morphology (e.g., length, diameter), cell wall composition (i.e., what polymers are present and in what quantities) and organization (i.e., where in the cell wall the polymers are located and how they are organized and interact with one another). At the biological level, these aspects are determined by the combined and coordinated spatio-temporal expression of several hundreds of genes, and a major challenge is therefore to identify which ones are involved in this process and to understand their role(s). As a result, genomic and genetic studies aimed at improving fiber plants have so far mainly focused on genome sequencing and the identification of genes/traits associated with various agronomic traits (e.g., plant height, fiber yield, fiber diameter) or cell wall construction (e.g., cellulose/pectin/lignin biosynthesis). Flax was the second fiber plant to have its genome sequenced (Wang et al., 2012) after hemp (van Bakel et al., 2011) and can be accessed on the phytozome public database1. The sequences of other fiber plants followed a few years later, most likely as a result of the size and/or complexity of the genome with cotton, jute and ramie being sequenced, respectively in 2015 (Li et al., 2015), 2016 (Yuan et al., 2016), 2017 (Sarkar et al., 2017), and 2018 (Luan et al., 2018). In parallel, a number of whole genome transcriptomics studies by microarrays and more recently RNAseq, together with targeted qRT-PCR and in situ hybridization, have enabled the identification of genes likely to be involved in cell wall polymer biosynthesis in flax (Roach and Deyholos, 2007; Fenart et al., 2010; Huis et al., 2012; Chantreau et al., 2015; Zhang and Deyholos, 2016; Gorshkov et al., 2017; le Roy et al., 2017). Despite the identification of numerous candidates, examples of functional validation of cell wall genes in flax, as in other fiber species – are relatively limited despite the fact that this species can be genetically transformed and mutant populations are available (Wróbel-Kwiatkowska et al., 2007; Day et al., 2009; Chantreau et al., 2013, 2014).
Linking fiber phenotype (morphology, composition and organization) to gene expression profiles in field-grown flax is an extremely challenging task as these profiles are continually modified during plant growth according to developmental and environmental cues (Figure 1). Both microarrays and RNAseq have been used to investigate modifications in gene expression in response to drought stress (Dash et al., 2014), saline and alkaline stress (Yu et al., 2014), and infection by the fungus Fusarium oxysporum (Galindo-González and Deyholos, 2016). In both cases, comparison of the transcriptomes from control plants and stressed plants revealed significant changes in the expression of a number of genes likely to have a direct/indirect impact on fiber phenotype.
Figure 1. Holistic overview of intrinsic and extrinsic factors contributing to flax fiber phenotype and different scales of (meta-)omic investigation. Numbers 1–4 represent different stages in the development of the flax plant and number 5 represents the retting stage. Number 6 represents the soil microbial community that will interact with plant roots during plant growth (number 7) and with the plant stem during retting (number 8). Biotic effects (left side) represent how plant organic material (e.g., leaves) affect soil microbiota. Biotic effects (right side) represent how microbiota affect the flax plant during retting (progressive degradation of cell walls). Abiotic effects (soil) represent how soil parameters (e.g., soil type, pH, mineral disponibility, water content) affect the plant during growth. Abiotic effects (climate) represent how meteorological conditions affect plant growth, soil parameters, and microbiome during plant growth (numbers 1–4) and retting (number 5). Omic approaches indicate different scales of investigation that can be used to investigate intrinsic (plant) biology during growth (numbers 1–4). Meta-omic approaches indicate different scales of investigation that can be used to characterize the microbial community during plant growth and retting (numbers 1–5). This diagram underlines the importance of considering the plant-microbiome-environment as a ‘holobiont’ for a complete understanding of how numerous different factors contribute to shaping fiber phenotype.
The situation is further complicated by the fact that flax fibers are extracted from the stem by mechanical defibering that may alter their mechanical and chemical properties especially at the level of so-called ‘technical fibers’, corresponding to whole or fragmented fiber bundles. Furthermore, mechanical extraction itself is generally preceded by a retting step that can also modify fiber properties if not managed properly (Md.Tahir et al., 2011; Liu et al., 2017).
In field-/dew-retting, plants are pulled (uprooted) and left on the ground for several weeks during which they are colonized by microorganisms (e.g., fungi and bacteria) that produce cell wall degrading enzymes contributing to the progressive de-solidification of fiber bundles from the other stem tissues and partial dissociation of the fiber bundles (Figure 2A). Fibers can also be de-solidified by water-retting, a process in which flax plants are placed into water tanks and colonized mainly by anaerobic bacteria (Tamburini et al., 2003; Zhao et al., 2016). In both cases these natural biological processes affect fiber homogeneity and mechanical performance.
Figure 2. Multiscale presentation of retting process in flax. (A) Microphotographs illustrating progressive changes in flax stem morphology during retting. Appearance of a stem cross-section from a freshly up-rooted plant at the beginning of retting (left), at the optimal retting (middle), and over-retted (right). At optimal retting cortical parenchyma cells are no longer visible and fiber bundle dissociation has started. Bast fibers (f) and xylem (x). (B) Schematic representation of the different layers of a plant cell wall: middle lamella (ML), primary cell wall (PW), secondary cell wall (SW), S1 and S2 layers of the secondary cell wall (S1, S2). All cell wall layers are intact at the start of retting (left), at optimal retting stage (middle) the middle lamella and primary cell walls of cortical parenchyma cells are degraded (green/orange) and the middle lamella between neighboring fiber cells is starting to become degraded (green/orange). At over-retting, the secondary cell wall of fibers starts to become degraded (right). (C) Examples of different cell wall polymers: Rhamnogalactan I (RGI) pectin (left), xyloglucan hemicellulose (middle), and cellulose (right). Keys indicate sugar monomers making up polymers (top frame) and cell wall degrading enzyme activities (bottom frame). This part was adapted from Polizeli et al. (2016).
Until recently, our knowledge of the different microorganisms present during retting, the dynamics of colonization, and the complexity of the cell wall degrading enzyme arsenal has been severely limited by our incapacity to cultivate and study the great majority of these species. Furthermore, we know hardly anything about how assembly, activities (i.e., cell-wall degrading capacity) of the retting microbial community and interactions among members are affected by environmental conditions or by plant genotype (Figure 1). The fact that such considerations are important is demonstrated by the results of several studies that have shown that both abiotic factors (e.g., climate, soil properties) and biotic factors (microorganisms) continue to affect the phenotype of plant fibers after harvest during the process of retting (Martin et al., 2013; Liu et al., 2015; Mazian et al., 2018; Bleuze et al., 2020; Chabbert et al., 2020). As a result, retting of fiber species is largely based upon the farmer’s experience and remains a biological black box preventing the development of more objective approaches.
Currently, rapid advances in high-throughput sequencing (HTS) technologies are allowing us to exploit the rapidly expanding field of metagenomics to improve our understanding of retting in flax and other fiber plants (Djemiel et al., 2017; Chabbert et al., 2020). In this paper we review how the application of this approach promises to make a major contribution to the understanding of this complex process and we discuss the perspectives of how related strategies such as metatranscriptomics and metaproteomics could also lead to further advances.
The Retting of Flax, a Key Step for the Obtention of Quality Fibers
Retting – Definition and Types
One of the first definitions of retting is as follows: ‘Retting refers to the operation that textile plants are subjected to in order to free the fibers of the liber from the gum-resinous substance which binds them and keeps them attached to the woody stem of the plant’ (Renouard, 1890). The ‘art’ of retting is all about compromise: under-retting will give coarse technical fibers, consisting of several associated individual fibers and contaminated by wood debris and external tissues (cortical parenchyma and cutinized epidermis; Goodman et al., 2002), while over-retting will reduce the resistance of the fibers by affecting their integral structure (Rosemberg and de França, 1967; Akin et al., 1998).
Generally, two main types of retting are used: “water retting” and “field-/dew-retting.” For water retting, flax stems are harvested, grouped in bundles and submerged in water. Historically, this was done in natural basins (e.g., lakes, rivers or dams) for 5–7 days, followed by a period where stems were dried directly on the ground for one to 2 weeks (Akin, 2010). More recently, artificial pits or tanks have been used to wash the stems with clean water to remove residues. This technique makes it possible to control parameters influencing bacterial development including aeration and temperature, and also allows the inoculation of selected bacterial strains. Water-retting depends mainly on the action of anaerobic bacteria that colonize the stem and is related to fermentation (Donaghy et al., 1990). This method generates high quality flax fibers but is also a source of environmental pollution (Md.Tahir et al., 2011). In field-/dew-retting, flax plants are mechanically up-rooted (pulled) by specialized harvesters and the stems placed directly on the field in piles to form swathes (Figure 1; Meijer et al., 1995). Subsequently, the humidity brought by the morning dew and the alternation of rain and heat, favors the development of microorganisms (mainly bacteria and fungi) initially present on the stems and the colonization by the microflora from the soil on which the swathes are placed (Sharma and Faughey, 1999; Akin, 2010; Djemiel et al., 2017). During retting, the intervention of the farmer involves turning the swathes halfway through the process in order to ensure uniform retting over the entire swathe height, and surveying the progression of the process (Franck, 2005; Akin, 2013).
More recently, various alternative methods using commercial enzyme cocktails mainly containing pectinases and xylanases have been investigated in an attempt to provide better control of the process and therefore improve fiber yield and quality (Henriksson et al., 1997; Zhang et al., 2000; Akin et al., 2001). However, their use remains limited mainly because of the high cost of producing these enzymes in bioreactors.
Retting and Fiber Cell Walls
Individual flax fiber cells (elementary fibers) are grouped together to form fiber bundles that are found in the stem outer tissues surrounding the xylem and pith that constitute the inner tissues (Figure 2A). The bundles are located beneath the epidermis and between the cortical parenchyma and the phloem. During retting, a number of morphological changes occur, mainly in the outer tissues of the stem (Figure 2A). Parenchyma cells around and between fiber bundles that can be easily observed during the early stages of retting almost disappear at later stages (Chabbert et al., 2020). Fiber bundles themselves also become progressively destructured mainly due to the dissociation of individual fiber cells (Akin, 2003; Figure 2A). Individual fiber cells are integrated into stem tissues via their compound middle lamella that not only links them to other fibers, but also links the peripheral bundle cells to the cortical cells in the external stem tissues (Meijer et al., 1995; Akin, 2013; Figure 2B). Retting therefore involves the progressive enzymatic degradation of both cortical parenchyma cell walls and the fiber compound middle lamella. Chemically these structures are composed of different polymers. Cellulose, hemicellulose (mainly xyloglucan), pectic polysaccharides (homogalacturonan, HG), rhamnogalacturonan I (RG-I), and rhamnogalacturonan II (RG-II), as well as some structural proteins are found in the primary cell wall. In contrast, the middle lamella is mainly composed of pectins forming what is commonly called the cellular cement. Fibers cells, in addition to a middle lamella and primary cell wall, have thick secondary walls that are rich in cellulose (up to 80%) and contain non-cellulosic polysaccharides such as galactans and glucomannans (Morvan et al., 2003; Rihouey et al., 2017). Low levels of lignin (2–5%) can also be sometimes found in the middle lamella, primary cell wall, and S1 layer of the secondary cell wall (Day et al., 2005). The technical properties of fibers are due to a high content of crystalline cellulose, with microfibrils almost parallel to the main axis of the cells (Müller et al., 1998). The Figure 2B shows the general organization of the plant cell wall.
The degradation of cell wall polymers necessitates the intervention of enzymes belonging to a wide range of different enzyme families (Figure 2C). In recent studies on hemp and flax retting, exopolygalacturonase and α-L-arabinosidase, β-D-xylosidase and β-D-galactosidase, β-D-glucosidase and cellobiohydrolase activities were followed during the retting process (Bleuze et al., 2018; Chabbert et al., 2020). Pectin-degrading enzymes are distributed across nine families of carbohydrate-active enzymes (CAZymes), CE8, PL1, PL2, PL3, PL9, PL10, GH28, GH78, and GH88. GH28 polygalacturonases (PGAs) play an important role in the breakdown of pectins in pathogenic ‘fungi.’ Pectin esterases (family CE8) catalyze the de-esterification of pectins into pectate and methanol (Zhao et al., 2013, 2014). Several strains of filamentous ‘fungi,’ isolated from flax stems retted in the field, such as Rhizomucor pusillus (Mucoromycota formerly Zygomycota) and Fusarium lateritium (Ascomycota) have shown a capacity for decohesion of fibers, in particular thanks to a high level of pectinolytic activity (Henriksson et al., 1997). Pectinolytic activities (e.g., polygalacturonases or pectin lyases) have also been observed in various microorganisms during water-retting (Zhao et al., 2016). For example, Clostridium felsineum shows a strong pectinolytic activity and a good capacity for retting. Studies using enzymatic retting have shown that the use of polygalacturonases alone is sufficient for the decohesion of fibers (Evans et al., 2002; Akin et al., 2004). Some microorganisms used during retting such as Bacillus subtilis and Erwinia carotovora (syn. Pectobacterium carotovorum) are known to produce pectin lyases (Sharma, 1986). Saprophytic or parasitic microorganisms have CAZymes showing activities of hemicellulolytic types (vanden Wymelenberg et al., 2010). These activities, in particular those of xylanases, have been observed during water-retting (Donaghy et al., 1990), or in bacterial cultures resulting from the cultivation of dew retted flax stems (Sharma, 1986). Very recently, a study on retting with water in a fermenter, showed a cyclic evolution of mannanases, another hemicellulolytic enzyme (Zhao et al., 2016).
During retting the enzymatic degradation of cell wall polymers can be followed by monitoring the sugars released. For example, the release of fibers was correlated with a decrease in the concentration of galacturonic acid resulting from the breakdown of pectic material in the middle lamella at the end of retting (Rosemberg and de França, 1967). A recent study of retting dynamics by both scanning electron microscope (SEM) imaging and the analysis of certain wall polymers also confirmed the link between the degradation of primary cell walls and fiber middle lamella and the dissociation of fiber bundles (Chabbert et al., 2020). In contrast glucose and mannose/galactose levels remain stable suggesting that secondary cell wall polysaccharides are not degraded (Chabbert et al., 2020). The release of other cell wall polymers such as cutin, wax, and aromatic compounds (sinapyl alcohol and ferulic acid) has also been investigated by chemical analyzes and mass spectrometry in a comparative study of water- and field-retting (Morrison et al., 2000).
Factors Affecting Retting
At the plant level, the amount of pectins varies depending on the cultivar and will influence the duration of retting (Brown et al., 1986; Haag et al., 2017). Increasing maturity of the stems will also favor lignin deposition in the fiber middle lamella associated with difficulties in retting (Meijer et al., 1995; Pallesen, 1996). However, the influence of exogenous factors, either abiotic or biotic, also has an important impact on retting. In the first category, can be classified all the pedoclimatic factors such as the soil structure or the mineral elements (for example nitrogen, phosphorus and potassium), the physico-chemistry of the soil, the rotation of the cultures, the thickness of the swathe, as well as climatic conditions and seasonal variations (Sharma and Faughey, 1999). In the second category, the most influential factor will undoubtedly be the microbiota present in the soil and in/on the stem at the beginning of retting (Djemiel et al., 2017; Chabbert et al., 2020). Intriguingly, a recent study has revealed the existence of a very close link between the microbiome and the spatial and temporal development of the plant (composition, physiology), or even cultural conditions (Comeau et al., 2020). In the light of such observations it seems likely that plant growth and fiber development in flax will also be affected by the microbiome present. Furthermore, the microbiome itself will be influenced by plant growth and retting conditions. A more integrated approach is clearly necessary if we are to fully understand all of the factors affecting fiber quality (Figure 1).
Microbiology of Retting, What Have We Learned From Classical Approaches?
Since it is the microorganisms that produce the hydrolytic enzymes responsible for retting, an important step in our understanding of this process is to identify the organisms responsible. Historically, various bacteria and fungi were identified in a number of different studies using isolation and culture-based approaches. However, such strategies are 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. Moreover, these approaches are generally inappropriate for dynamic studies looking at how microbial communities evolve during retting. More recently, microbial retting studies have greatly benefited from the use of HTS technologies that have produced exhaustive inventories of bacteria and fungi linked with this process. A list of microorganisms identified by classical methods is given in Tables 1, 2 and those identified by metabarcoding approaches in Table 3.
Table 1. List of different bacterial species involved in flax retting and identified by classical approaches.
Table 2. List of different fungal species involved in flax retting and identified by classical approaches.
Table 3. Non-exhaustive list of studies on the microbiology of retting performed using a metabarcoding approach.
The first reported study focusing on the retting of flax with water was carried out by the French biologist and botanist Philippe Édouard Léon Van Tieghem (1839–1914), who concluded that the bacteria Bacillus amylobacter was probably responsible for the decomposition of pectins (Figure 3A; van Tieghem, 1879). The second genus to be associated with water-retting was the genus Clostridium and more particularly Clostridium felsineum (Figure 3A; Lanigan, 1950). A few years later, Rosemberg (1965) isolated 22 individuals including species of the Clostridium and Pseudomonas genera. Pseudomonas aeruginosa was recognized as the fastest species in fiber decohesion (Rosemberg, 1965). This study also identified, for the first time, Achromobacter parvulus which appears to be involved in the latter stages of retting.
Figure 3. Timelines showing the increasing number of bacterial and fungal species identified over time in flax water-retting (A) and dew-retting (B). Figure shows authors and year of published studies (top line), microbial organization level (individual, community, ecosystem) and name (middle line), culture method and type of analysis (bottom line). The color of phyla (HTS approach) and species names indicates the identification method type (classical vs. HTS). When phyla names are in black it means that they have been identified for the first time by metabarcoding approaches.
Later, Sharma (1986) isolated and identified six bacterial species from the Actinobacteria, Firmicutes, and Proteobacteria phyla as well as five fungal species belonging to the Ascomycota phylum (Figure 3B).
In 1990, a study at the University of Ulster (Northern Ireland) followed the evolution of anaerobic bacteria during water-retting, both at an industrial and laboratory scale. Bacillus licheniformis and Bacillus subtilis were the most dominant during the first and second phase, corresponding to the growth phase and the rapid pectinolytic phase (between 10–40 h) of water-retting. Clostridium acetobutylicurn and Clostridium felsineum appear during the last phase known as the slow pectinolytic phase (Figure 3A; Donaghy et al., 1990). All the bacteria considered belong to the Firmicutes phylum.
In 2003, another study provided further information about pectinolytic bacteria involved in water-retting. Although this study still involved a culture step to isolate bacteria, it was the first one to use a molecular marker approach, in this case by amplifying a partial region of the 16S rRNA gene, using the (amplified ribosomal DNA restriction analysis) ARDRA technique. All of the anaerobic strains were assigned to the Clostridium genus and the aerobic strains to the Bacillus or Paenibacillus genera. Anaerobic colonies with significant polygalacturonase activity belonged to two phylogenetic clusters assigned to the Clostridium acetobutylicum/Clostridium felsineum and Clostridium saccharobutylicum species. For aerobic bacteria, colonies with significant polygalacturonase activity belonged to two phylogenetic clusters assigned to the Bacillus subtilis species (Figure 3A; Tamburini et al., 2003). All these strains also belong to the Firmicutes phylum.
One of the first studies on the microbiology of dew-retting identified several genera and species from the Fungi kingdom belonging to the Ascomycota and Mucoromycota (formerly Zygomycota) phyla (Figure 3B; Brown, 1984). In the late 1990s, an American study looked at fungi involved in flax dew-retting in the United States, Netherlands, and France (Figure 3B; Henriksson et al., 1999). Seven strains of filamentous fungi (including six still not described in the literature to this day) and a yeast, could be identified following isolation, cultivation on synthetic media, and purification in order to compare their activities and their efficiency in fiber release. Another study involving in vitro retting tests was carried out by Fila and coworkers but no new fungal genus was identified (Figure 3B; Fila et al., 2001).
Some fungi have been characterized as over-retting actors, for example Fusarium lateritium and more particularly Epicoccum nigrum, recognized as being a primary saprophyte in retting (Brown, 1984; Henriksson et al., 1997; Akin et al., 1998). Another species, Rhizomucor pusillus has been observed to degrade part of the surface of the cuticle, thus probably allowing the entry of microorganisms (Henriksson et al., 1997).
The overall conclusion that can be drawn from these classical, culture-based studies of retting is that the type of microorganism identified depends heavily upon the retting type. Field retting is an aerobic environment and although both bacteria and fungi are identified it is generally the latter that predominate with species belonging mostly to the Ascomycota phylum and to a lesser extent the Mucoromycota phylum (formerly Zygomycota; Table 2). In contrast, water-retting mainly involves anaerobic microorganisms (bacteria). Nevertheless, several fungi have also been identified in various anaerobic environments and classified within two very close phyla, Neocallimastigomycota and Chytridiomycota (Griffith et al., 2010; Gruninger et al., 2014). Interestingly, some of these anaerobic fungi, are also found in the rumen, and have been shown to possess a number of GH genes probably resulting from horizontal transfers of bacterial genes (Garcia-Vallvé et al., 2000; Steenbakkers et al., 2001).
Metabarcoding and the Breakthrough in Knowledge on Retting Microbial Diversity
Metabarcoding, also known as targeted metagenomics, is an approach that has revolutionized the study of microbial communities. In recent years the decrease in sequencing costs has allowed it to become more democratic, leading to the use of this approach in a larger number of laboratories and studies. The large and growing number of publications that use metabarcoding reflects this evolution regardless of the habitat explored, including studies on the microbiota of plants (Figure 4). Metabarcoding consists in sequencing an amplicon (most often corresponding to one or more regions of a phylogenetic marker) on a high−throughput platform to characterize microbial diversity from complex/environmental samples, thereby enabling studies of alpha-diversity, community structure or taxonomy (Deiner et al., 2017; Piper et al., 2019; Tedersoo et al., 2019).
Figure 4. Graph indicating the number of publications containing the keywords ‘metabarcoding’ (pink) or ‘metabarcoding AND plant’ (green) published each year since 2011.
There are currently very few studies on the holistic interactions between microorganisms and the flax plant during its growth and their implication on the quality of fibers. The only interactions studied for the moment are those taking place during retting, and that concern the capacity of microorganisms to degrade the plant cell wall and therefore in fine to affect fiber quality for their future use (Table 3). The first investigation using the metabarcoding approach investigated bacterial communities during kenaf water-retting under different controlled conditions (Visi et al., 2013). The authors used the 16S rRNA gene as the molecular marker and performed sequencing with an Ion Torrent Personal Genome Machine (PGM) system. The operational taxonomic unit (OTU) richness values from Visi et al. (2013) seem to be extremely high compared to other water-retting studies (Table 3). More recently, the evolution of both bacterial and fungal communities during water-retting of kenaf was studied using 18S rRNA gene (Duan et al., 2020). However, the results obtained concerning the Fungi kingdom must be taken with care because precisions concerning the bioinformatics analysis are not available (e.g., sub-sampling size/rarefaction) and certain inaccuracies on the taxonomic affiliation occured (e.g., confusion between the different taxonomic levels; Opisthokonta has been classified as a subkingdom and a phylum). To date, only one study has been conducted on hemp retting by metabarcoding (Law et al., 2020). This study indicated an average OTU richness mean of one hundred. In comparison, the study by Ribeiro et al. (2015) identified 66 different bacterial phylotypes. Despite the fact that these two studies involved different environmental conditions –Law et al. (2020) conducted hemp retting in a greenhouse under controlled conditions, while Ribeiro et al. (2015) collected samples in the field – the OTU or phylotype distributions were similar for the Proteobacteria (85 vs. 77%), Bacteroidetes (8 vs. 11%) and Firmicute (1% or less) but differed for Actinobacteria (2.3 vs. 11%).
Until now, it is flax retting that has been the most studied by metabarcoding, with two studies conducted on field retting (Djemiel et al., 2017) and one on water-retting in tanks (Zhao et al., 2016; Figure 3A). Zhao et al. (2016) focused on the evolution of anaerobic bacteria using a partial 16S rRNA gene sequence region following the addition of a strain of Bacillus cereus HDYM-02, previously characterized in their laboratory as being very effective in the decohesion of fibers. Interestingly, this species had been identified previously as an efficient retting agent that did adversely affect fiber mechanical properties (Donaghy et al., 1990). The study by Zhao et al. (2016) demonstrated the power of HTS since they were able to identify five phyla (Actinobacteria, Cyanobacteria, Bacteroidetes, Firmicutes, and Proteobacteria) and eight classes (Flavobacteriia, Sphingobacteriia, Bacilli, Clostridia, Negativicutes, α-proteobacteria, β-proteobacteria, and γ-proteobacteria), some of which were not previously known to be involved in water-retting (e.g., Cyanobacteria and Bacteroidetes for phyla, and Flavobacteriia, Sphingobacteriia, Negativicutes, α-proteobacteria, and β-proteobacteria for classes; Zhao et al., 2016).
The interest of the metabarcoding approach for flax field-retting was also demonstrated by two recent studies (Djemiel et al., 2017; Chabbert et al., 2020). The first study explored the diversity and dynamics of bacterial and fungal communities involved in this process by using an HTS DNA metabarcoding approach (16S rRNA/Internal Transcribed Spacer (ITS) region, Illumina Miseq) on plant and soil samples collected over a period of 7 weeks in July and August 2014 (Djemiel et al., 2017). Twenty-three bacterial and six fungal phyla were identified in soil samples and 11 bacterial and four fungal phyla in plant samples. Dominant phyla identified 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. The use of this powerful method also allowed the identification of rare phyla never previously associated with retting: Acidobacteria, CKC4, Chlorobi, Fibrobacteres, Gemmatimonadetes, Nitrospirae, and TM6 (bacteria), and Chytridiomycota (fungi). These results perfectly illustrate the efficiency of metabarcoding for detecting phyla of low abundance. In addition, this study also revealed that the agricultural practice of swathe turning affected bacterial and fungal community structure, and most likely contributed to a more uniform retting. The use of the prediction tool Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) (Langille et al., 2013; Douglas et al., 2020); was also used on OTU tables to decipher the functional potential of bacterial communities. This analysis predicted a large collection of potential bacterial enzymes capable of hydrolyzing the backbones and side-chains of cell wall polysaccharides. These results demonstrate the interest of a combined metabarcoding and PICRUSt analysis to predict functional composition of bacterial communities. A similar approach reported by Chabbert et al. (2020) allowed a comparison of the evolution of bacterial and fungal retting communities under different environmental conditions. Eight out of the ten most abundant bacterial OTUs identified in this study were also among the top ten OTUs identified in the 2014 study by Djemiel et al. (2017) in the same geographical area, but in a different field and under different climatic conditions. In contrast, only half of the top ten fungal OTUs were retrieved. In this last study, enzyme activities were measured in addition to targeted metagenomics showing the presence of a high microbial diversity with potential enzymatic functions during retting. The comparison of these two studies could suggest that a year-to-year variability exists, perhaps more so in fungal communities compared to bacteria. Nevertheless, the identification of comparable enzymatic functions in the two studies would suggest a functional redundancy within the communities. In conclusion, the contribution of metabarcoding to the knowledge about the microbiology of retting is undeniable. In addition, the technique can also be adapted to ask questions about the living world in the widest possible way. For example, by refining the choice of primers metabarcoding can also be used to inventory often neglected organisms such as Archaea or microeukaryotes (Bahram et al., 2019).
Integrative Meta-Omic Approaches to Decipher the Metaphenome and Global Ecology of Dew-Retting
Our understanding of the microbiology of retting has clearly benefited from the use of targeted metagenomics/metabarcoding. In particular, this strategy has made it possible to make a breakthrough in our knowledge about microbial community diversity and structure during this process. However, while metabarcoding allows identification of microorganisms, the functional gene content of the community can only be predicted based on the available knowledge of the whole genome sequences of individual organisms (or phylogenetically related organisms) making up the community. Furthermore, metabarcoding cannot determine which functional genes are expressed. Nevertheless, recent advances in sequencing technologies and computational capabilities offer a unique way for different meta-omic approaches (e.g., metagenomics, metatranscriptomics, and metaproteomics) to fill the gap from taxonomic information to functions by generating information on the functional gene composition of the retting microbial community.
To our knowledge, while there are currently no publications reporting the use of these meta-omics approaches to study retting, recent studies have been performed on biologically similar (but not identical) processes such as litter degradation. This latter phenomenon is indeed commonly described as a breakdown of dead plant material, which is composed of cellulose, hemicellulose, lignin, pectin but also proteins and can therefore be compared to the retting process. In the following sections we present some examples of how different meta-omic approaches have been used to look at plant biomass breakdown and discuss how these approaches could be exploited to provide functional information on retting.
Shotgun metagenomics is a high-throughput, culture-independent genomic analysis of all microbes present in a sample. It starts with extraction of total DNA from a sample followed by library preparation and sequencing. Shotgun metagenomics is thus an untargeted global sequencing approach compared to the 16S/ITS/18S rRNA sequencing used for metabarcoding (targeted metagenomics) and which is sometimes erroneously referred to as metagenomics. In the former approach, the aim is to sequence all the genomes present, whereas, in the metabarcoding strategy only the rRNA amplicons are sequenced. Coupled with functional profiling, shotgun metagenomics enables the description of the total genetic content of the microbial community. The functional annotation of the obtained sequences requires the use of several databases, including well-known broad databases such as NCBI nt, ortholog databases such as Clusters of Orthologous Groups (COG) and Evolutionary genealogy of genes: Non-supervised Orthologous Groups (EggNOG) for inference of functional categories, pathway databases such as Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2014) for metabolic pathway reconstruction and specific databases such as CAZy (Lombard et al., 2014) for CAZymes prediction. In a plant biomass context, shotgun metagenomics is principally applied to the mining for novel genes encoding potentially new/interesting CAZymes (Montella et al., 2017) driven by the need to find new biocatalysts for biofuel production. Shotgun metagenomics is also used for whole-genome reconstruction of individual species, the so-called metagenome-assembled genomes (MAGs). This approach allows researchers to simultaneously uncover a microorganism’s functional potential and its identity, with a more resolutive, unbiased phylogenetic information than with metabarcoding. However, the sequencing depth obtained in a sample is critical for successful microbial genome reconstruction and the complexity of the microbial community usually hinders such an approach (Alteio et al., 2020).
With regard to retting, no shotgun metagenomic studies have been conducted to date, and insights into genes involved in fiber degradation are still only obtained via culture-dependent genomic studies (Datta et al., 2020). Nevertheless, other environments associated with cell wall breakdown such as crop soils, and forest soils and litters have been studied (Table 4). Albeit the environmental conditions, cell wall structure and microbial community composition differ from those encountered during fiber retting, such studies show how the genetic potential of cell-wall decomposing microorganisms and their associated functions can be identified using shotgun metagenomics and illustrate how such approaches could be used in future to provide more in-depth information about retting.
Table 4. Selected metagenomic studies indicating number of reads using different sequencing technologies for a variety of materials from biological situations that present similarities to retting.
In a study of Amazonian rainforest and rainforest-to-pasture converted soils, Kroeger et al. (2018) showed (using a classical shotgun metagenomic strategy) that the microbial community undergoes drastic taxonomic and functional shifts with land-use type. Dominant phyla were common for both soils while less abundant phyla were strongly affected by the land-use type. In the pasture soil, gene functions were predominantly associated with carbohydrate metabolism, dormancy and sporulation, and regulation and cell signaling whereas in the rainforest soil, they were mostly associated with transcription and vitamin production. Although these findings do not constitute a breakthrough on the effect of different land-type use on microbial community composition, the metagenomics analysis was undoubtedly improved by the reconstruction of 28 MAGs from both dominant and rare phyla. Interestingly, a Melainabacteria MAG from a rare lineage was only detected in the pasture soil and shown to contain 72 protein encoding genes not found in other Melainabacteria genomes available. Among the Melainabacteria MAG genes annotated some were potentially involved in carbohydrate utilization. By gaining insights into the Melainabacteria MAG this study provides an example of how shotgun metagenomics can be used to explore the functional diversity of rare biospheres.
A MAG approach can also be applied to explore the functional diversity among members of the same phylum (Alteio et al., 2020). In this study, the authors reconstructed 67 putatively novel Bacteroidetes MAGs from warmed experimental forest soils and showed clade-specific diversity and abundance of CAZymes, suggestive of an extensive potential for polysaccharide degradation among members of the clade. In contrast, 17 unclassified Bacteroidetes MAGs lack major CAZymes required to efficiently deconstruct polysaccharides, suggestive of a rather limited role in polysaccharide degradation.
The study of the effect of agricultural management practices on crop soil microbial communities is also a growing field for the use of shotgun metagenomics and aims to correlate these practices to microbial community taxonomic diversity and the pool of functions (Souza et al., 2018; Nelkner et al., 2019). In their study, Souza et al. (2018) analyzed soil samples that were subjected to four different soil and crop management practices (conventional tillage and no tillage, each with crop succession or crop rotation). Although their results revealed a limited impact of management practice on the abundance of each functional category from the SEED database, the use of KEGG pathways identified some practice-related differences in key enzymes involved in the metabolism of fructose and mannose.
Another interesting example of the use of shotgun metagenomics to understand microbial degradation of plant cell wall polymers concerns the recent work by Kundu et al. (2019). The authors firstly obtained sequences from public shotgun metagenomics datasets of DNA extracted from different gut compartments of the termite Nasutitermes corniger. They then combined these data with metabolic data collected in the literature, resulting in a set of 2,988 metabolic reactions, that they used to construct a species-wide lignocellulose network. They identified 15 key bacterial species, all exclusively found in one gut compartment, producing a variety of enzymes that can breakdown lignocellulose including endoglucanases, exoglucanases, endo-1,4-β-xylanase, β-glucosidase, β-xylosidase, and cellobiohydrolase. In turn, the degradation products fed two other key bacterial species involved in fermentation.
Shotgun metagenomics offers a unique view of the functional potential (i.e., the genes present) of the microbial community and how it is affected (or not) by its environment. It does not, however, provide any information about the functions that are actually realized (i.e., the genes expressed and translated into functional proteins) by each member of the community. A first step in obtaining this information can be gained from gene expression studies through RNA sequencing. When applied to complex microbial communities this is called metatranscriptomics. Metatranscriptomics provides a snapshot of gene expression and constitutes a first step in going from ‘what the microbial community can do’ to ‘what the microbial community does.’ By analyzing gene expression at the community level, metatranscriptomics allows scientists to determine which metabolic pathways are active as well as their relative importance, thereby enabling predictions to be made about the functional role of the community. In addition, by accumulating metatranscriptomic data over the time course of a biological process it becomes possible to obtain a more accurate picture of the functional dynamics. The principle of metatranscriptomics is similar to the one described above for shotgun metagenomics except that total RNA (and not DNA) from a sample is extracted. The RNA then has to be processed for mRNA enrichment and/or rRNA depletion since rRNA are highly abundant and will compete for sequencing, resulting in poor representation of mRNA. As in all expression studies, particular precautions have to be taken to avoid RNA degradation during sample transport and preparation which can be challenging for field studies. A functional annotation of quality also relies on the use of several databases (including annotated metagenomes) to produce the most comprehensive transcript function prediction possible (Shakya et al., 2019).
As for shotgun metagenomics, metatranscriptomics has not yet been used in the context of retting. Nevertheless, some studies on biologically related environments have been done (Table 5). In their study, Hesse et al. (2015) looked at how CAZymes expression profiles were affected by long-term N deposition from leaf litter in forest soils from two geographically distant sites. They showed that the expression of CAZymes from bacterial or fungal origin responds to N amendment, with both common- and site-specific responses. Nitrogen input tended to favor CAZymes from bacterial origin. The authors also focused on the dynamic balance between Ascomycota and Basidiomycota fungi and showed a decrease of CAZymes from Basidiomycota in N-amended soils compared to ambient soils, that could be possibly linked to an overall decrease in litter decomposition under such conditions. Another valuable information that can be gained from metatranscriptomics concerns the role of other – generally ignored or at best, paid scant attention - members of the microbial community such as protists and viruses. Both of these groups of organisms have the ability to shape the microbial community by interacting with other members, and can thus indirectly modify lignocellulose decomposition activity. Current studies in soils have so far only focused on the taxonomic diversity of protists and viruses, revealing their extreme yet still unexplored diversity (Geisen et al., 2015; Starr et al., 2019). Further studies will be needed to investigate their role and interplay with the whole microbial community.
Table 5. Selected metatranscriptomic studies indicating number of reads using different sequencing technologies for a variety of materials from biological situations that present similarities to retting.
Thanks to advances in sequencing technologies and analysis tools in the last decade, more and more shotgun metagenomics and metatranscriptomics data can be obtained and analyzed thereby enabling in depth descriptions of microbial community structure and functions. The complexity of some environments such as soil that was previously a major drawback for analysis can now be alleviated by using other techniques to experimentally reduce the complexity of the community (Alteio et al., 2020). Despite current challenges, these approaches can be viewed, more than ever, as essential strategies on the path to decipher the role of microbial communities, and constitute a basis for other non-genomic-based omics approaches.
Like metabarcoding and metagenomics, metaproteomics can be used to provide information about microbial community structure and the different factors that can modify it. For example, it was recently shown that the long-term application of fertilization led to significant changes in microbial community structure and function associated with increased microbial biomass (Wu et al., 2015). In another example, the use of a stable isotope-probing metaproteomics approach using 15N-labeled plant-derived organic matter in a soil microbial community provided novel information on community structure and microorganism feeding strategies (Starke et al., 2016). This study revealed that Proteobacteria was the most abundant phylum followed by Actinobacteria and Ascomycota. It also demonstrated copiotrophic behavior for Rhizobiales belonging to Proteobacteria, Actinomycetales belonging to Actinobacteria and Chroococcales belonging to Cyanobacteria as these phylotypes immediately incorporated 15N from the added plant tissue. Conversely, the fungal Saccharomycetales and the bacterial Enterobacteriales, Pseudomonadales, Sphingomonadales, and Xanthomonadales displayed slower 15N-assimilation. Such a strategy could also be expected to generate information about the retting process in flax and other fiber species.
While metatranscriptomics provides information about the genes expressed in a community only a metaproteomics approach will provide precise information about the different proteins responsible for the observed function (in this case cell wall degrading enzymatic activity) in a given habitat (Herbst et al., 2016). Recent biochemical studies (Bleuze et al., 2018; Chabbert et al., 2020) have highlighted the role of different enzyme families during hemp and flax retting, but they did not identify the proteins responsible for these activities. The astounding advances of the past decade that have led to fast and affordable protein sequencing technologies associated with the increase in database information have led to a rapid increase in the number of protein identifications (Table 6). Protein data will not only give more direct information about microbial activities compared to metagenomics and metatranscriptomics, but can also fill taxonomical gaps of nucleotide-based methods (Wilmes and Bond, 2004, 2006).
Table 6. Selected studies performed on biologically-similar to retting systems reporting detection or identification of proteins using different proteo- and metaproteomic methods.
As is the case for metagenomics and metatranscriptomics, there are actually very few descriptions of a metaproteomics approach being applied to flax retting. Nevertheless, some general conclusions and useful information can be obtained from the studies on biologically similar systems. In a recent metaproteomics study of lodgepole pine decay, a diverse array of carbohydrate-active enzymes (CAZymes) were identified, representing a total of 132 families or subfamilies among which 672 glycoside hydrolases (GHs) including highly expressed cellulases or hemicellulases (Hori et al., 2018). The authors suggested that the observed enzymatic diversity and the coexistence of brown and white rot fungi indicated the existence of complex interactions between fungal species and degradation strategies. Another very recent metaproteomics study on the degumming of kenaf bast (outer-stem tissue containing fibers) probably represents the biological situation most closely related to retting (Duan et al., 2020). Here the authors analyzed kenaf bast fragments immersed into microbial fermentation liquid collected from different sites. Microbial secretomics analysis identified 197 proteins, including 67 differentially expressed proteins (DEPs) including Rds1, pyruvate kinase I and aconitate hydratase peptides. However, no DEPs associated with the degradation of cell wall polymers were observed underlying the difficulty of identifying non-core metabolism proteins.
Another important biological interest of using metaproteomics to study retting concerns the identification of protein termini for maturation or proteolytic processing (Westermann et al., 2017). After translation, proteins undergo a maturation step to obtain their final, active conformation. Apart from chemical modification of amino acid side chains, this maturation often involves cleavage of the polypeptide resulting in new N- or C-termini. The best-known example of this is the removal of the signal peptide from the N-terminus of the majority of secreted proteins, including most proteins involved in lignocellulose degradation. Of interest is the fact that non-signal peptide proteolytic processing has been described for plant cell wall interacting enzymes from the fungal genera Cladosporium and Cryptococcus (Ritch et al., 1991; Amoresano et al., 2000), identified as dominant during hemp retting (Ribeiro et al., 2015). Studies on proteins from Trichoderma reesei (Ståhlberg et al., 1988) the best studied cellulolytic fungi, indicate that proteolytic events are important in the control of the activity of cell wall lytic proteins (Coller et al., 1998). Similar observations have been made for bacterial enzymes with potential involvement in the partial degradation of plant cell walls, for instance, Erwinia pectate lyases are activated by N-terminal processing (Shevchik et al., 1998), as are Pseudomonas endoglucanases (Huang et al., 1989). For the study of proteolytic processing, different targeted approaches have been developed, the flagship being COFRADIC: Combined Fractional Diagonal Chromatography (Marino et al., 2015; Tanco et al., 2015) although more accessible alternatives have also been developed (Schilling et al., 2010).
Meta-Omics – Limitations and Best Practices
Before undertaking any experimental steps of any meta-omics approaches, it is firstly necessary to clearly define the biological question and research objectives (Matallana-Surget et al., 2018) in order to determine which meta-omics approach is best adapted.
While metabarcoding is an extremely interesting tool for identifying bacterial, fungal or other communities within a given sample, one of the main limitations of this approach is that it only delivers relative population data and is not quantitative (Lamb et al., 2019) regardless of the methodology, the primers, or the sequencing platform used. Inclusion of mock communities during metabarcoding is one way of getting around this problem. Metabarcoding is based on PCR amplification that can generate chimeras or other bias linked to this method and denoising data is therefore an important step in analysis that should be done with care. Recent pipelines (Qiime2 AND DADA2; Hall and Beiko, 2018) based on the notion of ASV (Amplicon Sequence Variant) are less likely to induce false discovery and seem more appropriate to identify diversity, replacing the notion of OTUs built by clustering. Generally, the errors and uncertainties associated with metabarcoding can often be tempered by careful study design, appropriate primer choice and robust sampling and replication (Murray et al., 2015). For example, the mix of plant and microorganism material collected during retting requires a careful choice of primer design to prevent the amplification of chloroplast/mitochondrial DNA (Djemiel et al., 2017). Another limitation of metabarcoding is its inability to describe the functional gene content of the community despite advanced functional inference tools such as PICRUSt2 (Douglas et al., 2020). Such information can only be predicted from whole genome sequence data of individual phylogenetically related organisms. However, when considering the cost/quantity of information ratio, metabarcoding is undoubtedly the most affordable way to approach the microbiology of retting.
In contrast to metabarcoding, shotgun metagenomics provides information about the functional potential (i.e., the genes present) in the microbial community. This approach is of course much more informative than metabarcoding but requires a much greater sequencing depth and therefore a much higher cost. Likewise, the need to assemble large-scale sequencing data requires skills and adequate computer resources. Extensive reviews of best-practices can be found in the literature and provide useful guidelines, from study design to data analysis, and address common pitfalls of such an approach (Knight et al., 2018; Bharti and Grimm, 2019).
If shotgun metagenomics can answer the questions “who is there” and “what enzymatic potential” is present, it cannot, however, provide any information about the functions that are really present (i.e., the genes expressed and translated into functional proteins). Such information can be obtained by metatranscriptomics. Despite the fact that the concept of metatranscriptomics has been around for several years (Gorni, 2013), there is currently no established protocol for this technique. Nevertheless, a number of recommendations can be made of which the most important is probably the necessity for ribodepletion in order to remove rRNAs from samples subjected to sequencing. This step remains critical and must be carefully controlled as the approaches used are not always successful (Barua, 2017). As with shotgun metagenomics, computer processing - in particular the assembly of the metatranscriptome, as well as the cost associated with sequencing remains an obstacle to the deployment of this approach on a large scale. Best practices for metatranscriptomics, especially regarding data analysis, are now starting to be published (Shakya et al., 2019), as the evolution of analysis tools grows.
If metatranscriptomics provides data on gene expression in microbial communities only metaproteomics can provide information about the different proteins responsible for the observed biological activity of the community (Herbst et al., 2016). Whatever the origin of the microbiome (human, plant, soil…), successful shotgun metaproteomics still faces a number of conceptual and technological hurdles that need to be overcome (Blackburn and Martens, 2016; Heyer et al., 2017; Knight et al., 2018; Matallana-Surget et al., 2018; Abiraami et al., 2020). Currently, a major limitation for metaproteomics in such systems is the lack of effective and reproducible protein extraction protocols and standardized data analyses resulting from i) the tremendous heterogeneity of samples (i.e., plant, soil and litter), ii) the low protein yield that can be obtained from soil and litter matrices and iii) the wide range of protein abundance levels, Becher et al. (2013), Keiblinger et al. (2016), Keiblinger and Riedel (2018).
Another issue concerns the size and the complexity of multi-organism protein sequence databases that are likely to contain many highly similar orthologs (Blackburn and Martens, 2016). There is therefore a global need for the improvement of database quality, including grouping of redundant proteins as well as taxonomic and functional annotation. A possible solution might be the generation of non-redundant fusion metagenomes for each type of microbial community (May et al., 2016). The use itself of protein databases could also be standardized, since some researchers use comprehensive protein databases and others use diverse metagenomes, which differ in the processing state and origin (Heyer et al., 2017). Currently, state-of-the- art metaproteomics studies only achieve identification of 5–30% spectra (Heyer et al., 2017) and it can be expected that a closer cooperation between bioinformaticians and biologists will improve bioinformatic strategies and increase the number of identified spectra. In the future, metaproteome studies will not only enable researchers to precisely characterize members of microbial and fungal communities but will also allow them to identify the retting hydrolytic enzymatic potential under given conditions. Approaches linking phylogenetics and functionality, could also be expected to help gain deeper insights into terrestrial microbial ecology.
In conclusion, by using integrative meta-omic approaches such as meta-genomics, metatranscriptomics and metaproteomics, a much better characterization and understanding of the dynamics of the retting microbial community and its molecular interactions can be obtained. However, for all methods, the same ‘bottlenecks’ remain: standardization of appropriate experimental designs to deal with the tremendous heterogeneity of samples, the use of state-of-the-art methodologies coupled with optimal genomic sequence information, pertinent database selection, and appropriate bioinformatic tools to analyze, integrate and visualize comprehensive global data sets (Becher et al., 2013; Tanca et al., 2013; Allen White et al., 2017). As indicated above, each technique focuses on a subset of the biological interaction network, therefore combining such contemporary molecular tools with eco-physiological studies in a trait-based framework (Krause et al., 2014) will allow scientists to more precisely decipher the global ecology of dew-retting. In a similar approach to that used for soil or rhizosphere microbial communities (El Amrani et al., 2015; Rochfort et al., 2015) the further inclusion of meta-metabolomic techniques would ultimately allow access to the signaling network thereby leading to an even deeper characterization and associated prediction of the retting metaphenome. Nevertheless, the interpretation of such large sets of data is notoriously difficult, mainly because of the difficulty of successfully integrating huge amounts of data. This difficulty underlines the absolute necessity to develop new methods for investigating dynamic aspects of large-scale models (Zivy et al., 2015; do Amaral and Souza, 2017; Wang et al., 2018) in a more integrative way (Lê Cao et al., 2008; Rohart et al., 2017).
Concluding Remarks and Future Perspectives
A multitude of intrinsic and extrinsic factors, both during plant growth and after harvest, all contribute to establishing the final phenotype of plant fibers reaching the factory (Figure 1). Since these factors affect fiber morphology, cell wall composition and organization they will also impact on the behavior of the separated fibers during subsequent industrial processing. In this paper we have focused on the different meta-omic technologies that have been, or could be, exploited to improve our understanding of one of these many factors - the field retting step of flax. This step involves a complex interaction between the harvested plant material and microorganisms, both in the soil and on the plant at the time of harvest. However, the complexity of this interaction does not stop there since the nature of the plant material itself (e.g., fiber morphology and cell wall composition), as well as the composition and functional capacity of the microbiome are also impacted by genetics and environmental conditions. In the light of such observations it is becoming clear that not only do we need to implement a multi-omics systems biology approach, but that we also need to take into account all of the actors that contribute to producing the final phenotype. This idea is embodied in the concept of the ‘holobiont’ which states that we should no longer consider the plant as an isolated ‘stand-alone’ organism, but rather as an individual with its associated microbial communities and in which their overall interactions are modulated by the pedo-climatic environment (Simon et al., 2019). Their association plays a decisive role in terms of the biodiversity and functionality of this ecosystem in which all the partners influence each other. Moreover, plant/microorganism (symbiotic, pathogenic, saprophytic) and microorganism/microorganism interactions will be affected by different biotic factors (e.g., species, cultivar, age, health, and stages of plant development) and abiotic factors (e.g., soil physicochemical composition, climatic conditions (Comeau et al., 2020; Fitzpatrick et al., 2020). As demonstrated in this review the combination of related ‘-omics’ such as metatranscriptomics and metaproteomics, together with biochemical studies, offers the unique opportunity to investigate the particular holobiont constituted by the overall fiber plant growing cycle. There is no doubt that the comprehension of this complex ecological machinery promises to make a major contribution to the control of plant natural fiber quality. For this purpose, it will be necessary to integrate together large multiscale data sets consisting of fiber parameters, biological and environmental information with the development of predictive models to be processed by AI analyses. The collection of massive data would rely on the development of connected microsensors reporting different parameters in real time directly from the field. Flax farmers will then have at their disposal a ‘retting toolbox’ making use of molecular markers, metabolite, protein and gene expression, biochemistry and phenotyping (morphological and agronomic) data for targeting selected traits. The cultivation of flax will then enter the era of smart connected agriculture.
SG and SH revised the manuscript. All authors were involved in the conceptualization of the manuscript and wrote the manuscript. All authors approved the final version of the manuscript.
This work was funded within the frameworks of the collaborative French ‘Future project’ SINFONI and CPER of region Hauts-de-France ALIBIOTECH. CD thanks the region of Hauts-de-France and Bpifrance and NB the region of Hauts de France for their financial support. BC thanks the region Grand-Est and Feder (Europe) for financial support. Publication fees were supported by University of Lille, CNRS, UMR 8576 -UGSF.
Conflict of Interest
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.
Abiraami, T. V., Singh, S., and Nain, L. (2020). Soil metaproteomics as a tool for monitoring functional microbial communities: promises and challenges. Rev. Environ. Sci. Technol. 19, 73–102. doi: 10.1007/s11157-019-09519-8
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
Allen, L. A. (1946a). Recent developments in the microbiology of retting. Part II.—The part played by bacteria in retting the flax and in purifying the liquor. J. Soc. Chem. Ind. 63:5. doi: 10.1002/jctb.5000630201
Allen White, R., Borkum, M. I., Rivas-Ubach, A., Bilbao, A., Wendler, J. P., Colby, S. M., et al. (2017). From data to knowledge: the future of multi-omics data analysis for the rhizosphere. Rhizosphere 3, 222–229. doi: 10.1016/j.rhisph.2017.05.001
Alteio, L. V., Schulz, F., Seshadri, R., Varghese, N., Rodriguez-Reillo, W., Ryan, E., et al. (2020). Complementary metagenomic approaches improve reconstruction of microbial diversity in a forest soil. mSystems 5:e00768–19. doi: 10.1128/mSystems.00768-19
Amoresano, A., Andolfo, A., Corsaro, M. M., Zocchi, I., Petrescu, I., Gerday, C., et al. (2000). Structural characterization of a xylanase from psychrophilic yeast by mass spectrometry. Glycobiology 10, 451–458. doi: 10.1093/glycob/10.5.451
Badalato, N., Guillot, A., Sabarly, V., Dubois, M., Pourette, N., Pontoire, B., et al. (2017). Whole proteome analyses on ruminiclostridium cellulolyticum show a modulation of the cellulolysis machinery in response to cellulosic materials with subtle differences in chemical and structural properties. PLoS One 12:e0170524. doi: 10.1371/journal.pone.0170524
Bahram, M., Anslan, S., Hildebrand, F., Bork, P., and Tedersoo, L. (2019). Newly designed 16S rRNA metabarcoding primers amplify diverse and novel archaeal taxa from the environment. Environ. Microbiol. Rep. 11, 487–494. doi: 10.1111/1758-2229.12684
Bailly, J., Fraissinet-Tachet, L., Verner, M. C., Debaud, J. C., Lemaire, M., Wésolowski-Louvel, M., et al. (2007). Soil eukaryotic functional diversity, a metatranscriptomic approach. ISME J. 1, 632–642. doi: 10.1038/ismej.2007.68
Bastida, F., Algora, C., Hernández, T., and García, C. (2012). Feasibility of a cell separation-proteomic based method for soils with different edaphic properties and microbial biomass. Soil Biol. Biochem. 45, 136–138. doi: 10.1016/j.soilbio.2011.10.017
Bastida, F., Jehmlich, N., Torres, I. F., and García, C. (2018). The extracellular metaproteome of soils under semiarid climate: a methodological comparison of extraction buffers. Sci. Total Environ. 61, 707–711. doi: 10.1016/j.scitotenv.2017.11.134
Becher, D., Bernhardt, J., Fuchs, S., and Riedel, K. (2013). Metaproteomics to unravel major microbial players in leaf litter and soil environments: challenges and perspectives. Proteomics 13, 2895–2909. doi: 10.1002/pmic.201300095
Benndorf, D., Balcke, G. U., Harms, H., and von Bergen, M. (2007). Functional metaproteome analysis of protein extracts from contaminated soil and groundwater. ISME J. 1, 224–234. doi: 10.1038/ismej.2007.39
Benndorf, D., Vogt, C., Jehmlich, N., Schmidt, Y., Thomas, H., Woffendin, G., et al. (2009). Improving protein extraction and separation methods for investigating the metaproteome of anaerobic benzene communities within sediments. Biodegradation 20, 737–750. doi: 10.1007/s10532-009-9261-3
Berlemont, R., Allison, S. D., Weihe, C., Lu, Y., Brodie, E. L., Martiny, J. B. H., et al. (2014). Cellulolytic potential under environmental changes in microbial communities from grassland litter. Front. Microbiol. 5:639. doi: 10.3389/fmicb.2014.00639
Bleuze, L., Chabbert, B., Lashermes, G., and Recous, S. (2020). Hemp harvest time impacts on the dynamics of microbial colonization and hemp stems degradation during dew retting. Ind. Crops Prod. 145:112122. doi: 10.1016/j.indcrop.2020.112122
Bleuze, L., Lashermes, G., Alavoine, G., Recous, S., and Chabbert, B. (2018). Tracking the dynamics of hemp dew retting under controlled environmental conditions. Ind. Crops Prod. 123, 55–63. doi: 10.1016/j.indcrop.2018.06.054
Bourmaud, A., Beaugrand, J., Shah, D. U., Placet, V., and Baley, C. (2018). Towards the design of high-performance plant fibre composites. Prog. Mater. Sci. 97, 347–408. doi: 10.1016/j.pmatsci.2018.05.005
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
Callister, S. J., Fillmore, T. L., Nicora, C. D., Shaw, J. B., Purvine, S. O., Orton, D. J., et al. (2018). Addressing the challenge of soil metaproteome complexity by improving metaproteome depth of coverage through two-dimensional liquid chromatography. Soil Biol. Biochem. 125, 290–299. doi: 10.1016/j.soilbio.2018.07.018
Chabbert, B., Padovani, J., Djemiel, C., Ossemond, J., Lemaître, A., Yoshinaga, A., et al. (2020). Multimodal assessment of flax dew retting and its functional impact on fibers and natural fiber composites. Ind. Crops Prod. 148:112255. doi: 10.1016/j.indcrop.2020.112255
Chantreau, M., Chabbert, B., Billiard, S., Hawkins, S., and Neutelings, G. (2015). Functional analyses of cellulose synthase genes in flax (Linum usitatissimum) by virus-induced gene silencing. Plant Biotechnol. J. 13, 1312–1324. doi: 10.1111/pbi.12350
Chantreau, M., Grec, S., Gutierrez, L., Dalmais, M., Pineau, C., Demailly, H., et al. (2013). PT-Flax (phenotyping and TILLinG of flax): development of a flax (Linum usitatissimum L.) mutant population and TILLinG platform for forward and reverse genetics. BMC Plant Biol. 13:159. doi: 10.1186/1471-2229-13-159
Chantreau, M., Portelette, A., Dauwe, R., Kiyoto, S., Crônier, D., Morreel, K., et al. (2014). Ectopic lignification in the Flax lignified bast fiber1 mutant stem is associated with tissue-specific modifications in gene expression and cell wall composition. Plant Cell 26, 4462–4482. doi: 10.1105/tpc.114.130443
Chen, S., Rillig, M. C., and Wang, W. (2009). Improving soil protein extraction for metaproteome analysis and glomalin-related soil protein detection. Proteomics 9, 4970–4973. doi: 10.1002/pmic.200900251
Chourey, K., Jansson, J., Verberkmoes, N., Shah, M., Chavarria, K. L., Tom, L. M., et al. (2010). Direct cellular lysis/protein extraction protocol for soil metaproteomics. J. Prot. Res. 9, 6615–6622. doi: 10.1021/pr100787q
Coller, S. P., Schoisswohl, D., Baron, M., Parriche, M., and Kubicek, C. P. (1998). Role of endoproteolytic dibasic proprotein processing in maturation of secretory proteins in Trichoderma reesei. Appl. Environ. Microbiol. 64, 3202–3208. doi: 10.1128/aem.64.9.3202-3208.1998
Cui, J., Mai, G., Wang, Z., Liu, Q., Zhou, Y., Ma, Y., et al. (2019). Metagenomic insights into a cellulose-rich niche reveal microbial cooperation in cellulose degradation. Front. Microbiol. 10:618. doi: 10.3389/fmicb.2019.00618
Damon, C., Lehembre, F., Oger-Desfeux, C., Luis, P., Ranger, J., Fraissinet-Tachet, L., et al. (2012). Metatranscriptomics reveals the diversity of genes expressed by eukaryotes in forest soils. PLoS One 7:e28967. doi: 10.1371/journal.pone.0028967
Dash, P. K., Cao, Y., Jailani, A. K., Gupta, P., Venglat, P., Xiang, D., et al. (2014). Genome-wide analysis of drought induced gene expression changes in flax (Linum usitatissimum). GM Crops Food 5, 106–119. doi: 10.4161/gmcr.29742
Datta, S., Saha, D., Chattopadhyay, L., and Majumdar, B. (2020). Genome comparison identifies different bacillus species in a bast fibre-retting bacterial consortium and provides insights into pectin degrading genes. Sci. Rep. 10:8169. doi: 10.1038/s41598-020-65228-1
Day, A., Neutelings, G., Nolin, F., Grec, S., Habrant, A., Crônier, D., et al. (2009). Caffeoyl coenzyme A O-methyltransferase down-regulation is associated with modifications in lignin and cell-wall architecture in flax secondary xylem. Plant Physiol. Biochem. 47, 9–19. doi: 10.1016/j.plaphy.2008.09.011
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 Vries, M., Schöler, A., Ertl, J., Xu, Z., and Schloter, M. (2015). Metagenomic analyses reveal no differences in genes involved in cellulose degradation under different tillage treatments. FEMS Microbiol. Ecol. 91:fiv069. doi: 10.1093/femsec/fiv069
Deiner, K., Bik, H. M., Mächler, E., Seymour, M., Lacoursière-Roussel, A., Altermatt, F., et al. (2017). Environmental DNA metabarcoding: transforming how we survey animal and plant communities. Mol. Ecol. 26, 5872–5895. doi: 10.1111/mec.14350
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
Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688. doi: 10.1038/s41587-020-0548-6
Duan, S. W., Cheng, L. F., Feng, X. Y., Yang, Q., Liu, Z. Y., Zheng, K., et al. (2020). Insights on bio-degumming of kenaf bast based on metagenomic and proteomics. BMC Genom. 21:121. doi: 10.1186/s12864-020-6531-2
El Amrani, A., Dumas, A. S., Wick, L. Y., Yergeau, E., and Berthomé, R. (2015). “Omics” insights into PAH degradation toward improved green remediation biotechnologies. Environ. Sci. Technol. 49, 11281–11291. doi: 10.1021/acs.est.5b01740
Evans, J. D., Akin, D. E., and Foulk, J. A. (2002). Flax-retting by polygalacturonase-containing enzyme mixtures and effects on fiber properties. J. Biotechnol. 97, 223–231. doi: 10.1016/S0168-1656(02)00066-4
Fenart, S., Ndong, Y.-P. A., Duarte, J., Rivière, N., Wilmer, J., van Wuytswinkel, O., et al. (2010). Development and validation of a flax (Linum usitatissimum L.) gene expression oligo microarray. BMC Genom. 11:592. doi: 10.1186/1471-2164-11-592
Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., and Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007–1017. doi: 10.1038/ismej.2011.159
Fitzpatrick, C. R., Salas-González, I., Conway, J. M., Finkel, O. M., Gilbert, S., Russ, D., et al. (2020). The plant microbiome: from ecology to reductionism and beyond. Annu. Rev. Microbiol. 74:annurev–micro–022620–014327. doi: 10.1146/annurev-micro-022620-014327
Freedman, Z. B., Upchurch, R. A., Zak, D. R., and Cline, L. C. (2016). Anthropogenic N deposition slows decay by favoring bacterial metabolism: insights from metagenomic analyses. Front. Microbiol. 7:259. doi: 10.3389/fmicb.2016.00259
Galindo-González, L., and Deyholos, M. K. (2016). RNA-seq Transcriptome Response of Flax (Linum usitatissimum L.) to the pathogenic fungus fusarium oxysporum f. sp. lini. Front. Plant Sci. 7:1766. doi: 10.3389/fpls.2016.01766
Geisen, S., Tveit, A. T., Clark, I. M., Richter, A., Svenning, M. M., Bonkowski, M., et al. (2015). Metatranscriptomic census of active protists in soils. ISME J. 9, 2178–2190. doi: 10.1038/ismej.2015.30
Geng, A., Cheng, Y., Wang, Y., Zhu, D., Le, Y., Wu, J., et al. (2018). Transcriptome analysis of the digestive system of a wood-feeding termite (Coptotermes formosanus) revealed a unique mechanism for effective biomass degradation. Biotechnol. Biofuels 11, 1–14. doi: 10.1186/s13068-018-1015-1
Giagnoni, L., Magherini, F., Landi, L., Taghavi, S., Modesti, A., Bini, L., et al. (2011). Extraction of microbial proteome from soil: potential and limitations assessed through a model study. Eur. J. Soil Sci. 62, 74–81. doi: 10.1111/j.1365-2389.2010.01322.x
Goodman, A. M., Ennos, A. R., and Booth, I. (2002). A mechanical study of retting in glyphosate treated flax stems (Linum usitatissimum). Ind. Crops Prod. 15, 169–177. doi: 10.1016/S0926-6690(01)00108-X
Gorshkov, O., Mokshina, N., Gorshkov, V., Chemikosova, S., Gogolev, Y., and Gorshkova, T. (2017). Transcriptome portrait of cellulose-enriched flax fibres at advanced stage of specialization. Plant Mol. Biol. 93, 1–19. doi: 10.1007/s11103-016-0571-7
Goss-Souza, D., Mendes, L. W., Borges, C. D., Rodrigues, J. L. M., and Tsai, S. M. (2019). Amazon forest-to-agriculture conversion alters rhizosphere microbiome composition while functions are kept. FEMS Microbiol. Ecol. 95, 1–13. doi: 10.1093/femsec/fiz009
Griffith, G. W., Baker, S., Fliegerova, K., Liggenstoffer, A., van der Giezen, M., Voigt, K., et al. (2010). Anaerobic fungi: neocallimastigomycota. IMA Fung. 1, 181–185. doi: 10.5598/imafungus.2010.01.02.11
Gruninger, R. J., Puniya, A. K., Callaghan, T. M., Edwards, J. E., Youssef, N., Dagar, S. S., et al. (2014). Anaerobic fungi (phylum Neocallimastigomycota): advances in understanding their taxonomy, life cycle, ecology, role and biotechnological potential. FEMS Microbiol. Ecol. 90, 1–17. doi: 10.1111/1574-6941.12383
Haag, K., Padovani, J., Fita, S., Trouvé, J.-P., Pineau, C., Hawkins, S., et al. (2017). Influence of flax fibre variety and year-to-year variability on composite properties. Ind. Crops Prod. 98, 1–9. doi: 10.1016/j.indcrop.2016.12.028
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. Textile Res. J. 67, 829–836. doi: 10.1177/004051759706701107
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
Herbst, F.-A., Lünsmann, V., Kjeldal, H., Jehmlich, N., Tholey, A., von Bergen, M., et al. (2016). Enhancing metaproteomics-The value of models and defined environmental microbial systems. Proteomics 16, 783–798. doi: 10.1002/pmic.201500305
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
Heyer, R., Schallert, K., Zoun, R., Becher, B., Saake, G., and Benndorf, D. (2017). Challenges and perspectives of metaproteomic data analysis. J. Biotechnol. 261, 24–36. doi: 10.1016/j.jbiotec.2017.06.1201
Hori, C., Gaskell, J., Cullen, D., Sabat, G., Stewart, P. E., Lail, K., et al. (2018). Multi-omic analyses of extensively decayed pinus contorta reveal expression of a diverse array of lignocellulose-degrading enzymes. Appl. Environ. Microbiol. 84:e01133–18. doi: 10.1128/AEM.01133-18
Huis, R., Morreel, K., Fliniaux, O., Lucau-Danila, A., Fenart, S., Grec, S., et al. (2012). Natural hypolignification is associated with extensive oligolignol accumulation in flax stems. Plant Physiol. 158, 1893–1915. doi: 10.1104/pp.111.192328
Kanehisa, M., Goto, S., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2014). Data, information, knowledge and principle: back to metabolism in KEGG. Nucl. Acids Res. 42, D199–D205. doi: 10.1093/nar/gkt1076
Keiblinger, K. M., Fuchs, S., Zechmeister-Boltenstern, S., and Riedel, K. (2016). Soil and leaf litter metaproteomics—a brief guideline from sampling to understanding. FEMS Microbiol. Ecol. 92:fiw180. doi: 10.1093/femsec/fiw180
Keiblinger, K. M., Wilhartitz, I. C., Schneider, T., Roschitzki, B., Schmid, E., Eberl, L., et al. (2012). Soil metaproteomics – Comparative evaluation of protein extraction protocols. Soil Biol. Biochem. 54, 14–24. doi: 10.1016/j.soilbio.2012.05.014
Kipping, L., Holzscheck, N., Maurer, F., Muszynski, S., Noll, M., and Jehmlich, N. (2020). Microbial metaproteome data from decayed beech dead wood. Data Brief 29:105285. doi: 10.1016/j.dib.2020.105285
Knight, R., Vrbanac, A., Taylor, B. C., Aksenov, A., Callewaert, C., Debelius, J., et al. (2018). Best practices for analysing microbiomes. Nat. Rev. Microbiol. 16, 410–422. doi: 10.1038/s41579-018-0029-9
Krause, S., le Roux, X., Niklaus, P. A., van Bodegom, P. M., Lennon, J. T., Bertilsson, S., et al. (2014). Trait-based approaches for understanding microbial biodiversity and ecosystem functioning. Front. Microbiol. 5:251. doi: 10.3389/fmicb.2014.00251
Kroeger, M. E., Delmont, T. O., Eren, A. M., Meyer, K. M., Guo, J., Khan, K., et al. (2018). New biological insights into how deforestation in amazonia affects soil microbial communities using metagenomics and metagenome-assembled genomes. Front. Microbiol. 9:1635. doi: 10.3389/fmicb.2018.01635
Kundu, P., Manna, B., Majumder, S., and Ghosh, A. (2019). Species-wide metabolic interaction network for understanding natural lignocellulose digestion in termite gut microbiota. Sci. Rep. 9:16329. doi: 10.1038/s41598-019-52843-w
Kuramae, E. E., Hillekens, R. H. E., de Hollander, M., van der Heijden, M. G. A., van den Berg, M., van Straalen, N. M., et al. (2013). Structural and functional variation in soil fungal communities associated with litter bags containing maize leaf. FEMS Microbiol. Ecol. 84, 519–531. doi: 10.1111/1574-6941.12080
Lamb, P. D., Hunter, E., Pinnegar, J. K., Creer, S., Davies, R. G., and Taylor, M. I. (2019). How quantitative is metabarcoding: a meta-analytical approach. Mol. Ecol. 28, 420–430. doi: 10.1111/mec.14920
Langille, M. G. I., 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
Law, A. D., McNees, C. R., and Moe, L. A. (2020). The microbiology of hemp retting in a controlled environment: steering the hemp microbiome towards more consistent fiber production. Agronomy 10:492. doi: 10.3390/agronomy10040492
Lê Cao, K. A., Rossouw, D., Robert-Granié, C., and Besse, P. (2008). A sparse PLS for variable selection when integrating omics data. Statist. Appl. Genet. Mol. Biol. 7:35. doi: 10.2202/1544-6115.1390
le Roy, J., Blervacq, A.-S., Créach, A., Huss, B., Hawkins, S., and Neutelings, G. (2017). Spatial regulation of monolignol biosynthesis and laccase genes control developmental and stress-related lignin in flax. BMC Plant Biol. 17:124. doi: 10.1186/s12870-017-1072-9
Li, F., Fan, G., Lu, C., Xiao, G., Zou, C., Kohel, R. J., et al. (2015). Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat. Biotechnol. 33, 524–530. doi: 10.1038/nbt.3208
Liu, M., Fernando, D., Daniel, G., Madsen, B., Meyer, A. S., Ale, M. T., et al. (2015). Effect of harvest time and field retting duration on the chemical composition, morphology and mechanical properties of hemp fibers. Ind. Crops Prod. 69, 29–39. doi: 10.1016/j.indcrop.2015.02.010
Liu, M., Thygesen, A., Summerscales, J., and Meyer, A. S. (2017). Targeted pre-treatment of hemp bast fibres for optimal performance in biocomposite materials: a review. Ind. Crops Prod. 108, 660–683. doi: 10.1016/j.indcrop.2017.07.027
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucl. Acids Res. 42, D490–D495. doi: 10.1093/nar/gkt1178
Luan, M. B., Jian, J. B., Chen, P., Chen, J. H., Chen, J. H., Gao, Q., et al. (2018). Draft genome sequence of ramie, Boehmeria nivea (L.) Gaudich. Mol. Ecol. Resour. 18, 639–645. doi: 10.1111/1755-0998.12766
Martin, N., Mouret, N., Davies, P., and Baley, C. (2013). Influence of the degree of retting of flax fibers on the tensile properties of single fibers and short fiber/polypropylene composites. Ind. Crops Prod. 49, 755–767. doi: 10.1016/j.indcrop.2013.06.012
Marynowska, M., Goux, X., Sillam-Dussès, D., Rouland-Lefèvre, C., Roisin, Y., Delfosse, P., et al. (2017). Optimization of a metatranscriptomic approach to study the lignocellulolytic potential of the higher termite gut microbiome. BMC Genom. 18:681. doi: 10.1186/s12864-017-4076-9
Matallana-Surget, S., Jagtap, P. D., Griffin, T. J., Beraud, M., and Wattiez, R. (2018). Comparative metaproteomics to study environmental changes. Perspect. Methods Appl. 2018, 327–363. doi: 10.1016/B978-0-08-102268-9.00017-3
May, D. H., Timmins-Schiffman, E., Mikan, M. P., Harvey, H. R., Borenstein, E., Nunn, B. L., et al. (2016). An alignment-free “metapeptide” strategy for metaproteomic characterization of microbiome samples using shotgun metagenomic sequencing. J. Prot. Res. 15, 2697–2705. doi: 10.1021/acs.jproteome.6b00239
Mazian, B., Bergeret, A., Benezet, J.-C., and Malhautier, L. (2018). Influence of field retting duration on the biochemical, microstructural, thermal and mechanical properties of hemp fibres harvested at the beginning of flowering. Ind. Crops Prod. 116, 170–181. doi: 10.1016/j.indcrop.2018.02.062
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
Morrison, I. I. I. W. H., Archibald, D. D., Sharma, H. S. S., and Akin, D. E. (2000). Chemical and physical characterization of water- and dew-retted flax fibers. Ind. Crops Prod. 12, 39–46. doi: 10.1016/S0926-6690(99)00044-8
Morvan, C., Andème-Onzighi, C., Girault, R., Himmelsbach, D. S., Driouich, A., and Akin, D. E. (2003). Building flax fibres: more than one brick in the walls. Plant Physiol. Biochem. 41, 935–944. doi: 10.1016/j.plaphy.2003.07.001
Müller, M., Czihak, C., Vogl, G., Fratzl, P., Schober, H., and Riekel, C. (1998). Direct observation of microfibril arrangement in a single native cellulose fiber by microbeam small-angle X-ray Scattering. Macromolecules 31, 3953–3957. doi: 10.1021/ma980004c
Murray, D. C., Coghlan, M. L., and Bunce, M. (2015). From benchtop to desktop: important considerations when designing amplicon sequencing workflows. PLoS One 10:e0124671. doi: 10.1371/journal.pone.0124671
Nelkner, J., Henke, C., Lin, T. W., Pätzold, W., Hassa, J., Jaenicke, S., et al. (2019). Effect of long-term farming practices on agricultural soil microbiome members represented by metagenomically assembled genomes (mags) and their predicted plant-beneficial genes. Genes 10:424. doi: 10.3390/genes10060424
Piper, A. M., Batovska, J., Cogan, N. O. I., Weiss, J., Cunningham, J. P., Rodoni, B. C., et al. (2019). Prospects and challenges of implementing DNA metabarcoding for high-throughput insect surveillance. GigaScience 8:giz092. doi: 10.1093/gigascience/giz092
Polizeli, M. L. T. M., Vici, A. C., Scarcella, A. S. A., Cereia, M., and Pereira, M. G. (2016). “Enzyme System from Aspergillus in Current Industrial Uses and Future Applications in the Production of Second-Generation Ethanol,” in New and Future Developments in Microbial Biotechnology and Bioengineering, ed. V. K. Gupta (Amsterdam: Elsevier), 127–140. doi: 10.1016/b978-0-444-63505-1.00009-9
Renouard, A. (1890). Etudes sur la culture, le rouissage et le teillage du lin dans tous les pays du monde: France, Belgique, Angleterre, Autriche, Allemagne, Etats-Unis. Paris: Librairie Polytechnique de A. Baudry, 32.
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
Rihouey, C., Paynel, F., Gorshkova, T., and Morvan, C. (2017). Flax fibers: assessing the non-cellulosic polysaccharides and an approach to supramolecular design of the cell wall. Cellulose 24, 1985–2001. doi: 10.1007/s10570-017-1246-5
Ritch, T. G., Nipper, V. J., Akileswaran, L., Smith, A. J., Pribnow, D. G., and Gold, M. H. (1991). Lignin peroxidase from the basidiomycete Phanerochaete chrysosporium is synthesized as a preproenzyme. Gene 107, 119–126. doi: 10.1016/0378-1119(91)90304-T
Roach, M. J., and Deyholos, M. K. (2007). Microarray analysis of flax (Linum usitatissimum L.) stems identifies transcripts enriched in fibre-bearing phloem tissues. Mol. Genet. Genom. 278, 149–165. doi: 10.1007/s00438-007-0241-1
Rochfort, S., Ezernieks, V., Mele, P., and Kitching, M. (2015). NMR metabolomics for soil analysis provide complementary, orthogonal data to MIR and traditional soil chemistry approaches-A land use study. Mag. Res. Chem. 53, 719–725. doi: 10.1002/mrc.4187
Rohart, F., Gautier, B., Singh, A., and Lê Cao, K. A. (2017). mixOmics: an R package for ‘omics feature selection and multiple data integration. PLoS Comput. Biol. 13:e1005752. doi: 10.1371/journal.pcbi.1005752
Sarkar, D., Mahato, A. K., Satya, P., Kundu, A., Singh, S., Jayaswal, P. K., et al. (2017). The draft genome of Corchorus olitorius cv. JRO-524 (Navin). Genom. Data 12, 151–154. doi: 10.1016/j.gdata.2017.05.007
Schneider, T., Gerrits, B., Gassmann, R., Schmid, E., Gessner, M. O., Richter, A., et al. (2010). Proteome analysis of fungal and bacterial involvement in leaf litter decomposition. Proteomics 10, 1819–1830. doi: 10.1002/pmic.200900691
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
Sharma, H. S. S., and Faughey, G. J. (1999). Comparison of subjective and objective methods to assess flax straw cultivars and fibre quality after dew-retting. Ann. Appl. Biol. 135, 495–501. doi: 10.1111/j.1744-7348.1999.tb00879.x
Shevchik, V. E., Boccara, M., Vedel, R., and Hugouvieux-Cotte-Pattat, N. (1998). Processing of the pectate lyase Pell by extracellular proteases of Erwinia chrysanthemi 3937. Mol. Microbiol. 29, 1459–1469. doi: 10.1046/j.1365-2958.1998.01028.x
Souza, R. C., Cantão, M. E., Nogueira, M. A., Vasconcelos, A. T. R., and Hungria, M. (2018). Outstanding impact of soil tillage on the abundance of soil hydrolases revealed by a metagenomic approach. Br. J. Microbiol. 49, 723–730. doi: 10.1016/j.bjm.2018.03.001
Souza, R. C., Cantão, M. E., Vasconcelos, A. T. R., Nogueira, M. A., and Hungria, M. (2013). Soil metagenomics reveals differences under conventional and no-tillage with crop rotation or succession. Appl. Soil Ecol. 72, 49–61. doi: 10.1016/j.apsoil.2013.05.021
Souza, R. C., Hungria, M., Cantão, M. E., Vasconcelos, A. T. R., Nogueira, M. A., and Vicente, V. A. (2015). Metagenomic analysis reveals microbial functional redundancies and specificities in a soil under different tillage and crop-management regimes. Appl. Soil Ecol. 86, 106–112. doi: 10.1016/j.apsoil.2014.10.010
Ståhlberg, J., Johansson, G., and Pettersson, G. (1988). A binding-site-deficient, catalytically active, core protein of endoglucanase III from the culture filtrate of Trichoderma reesei. Eur. J. Biochem. 173, 179–183. doi: 10.1111/j.1432-1033.1988.tb13982.x
Starke, R., Kermer, R., Ullmann-Zeunert, L., Baldwin, I. T., Seifert, J., Bastida, F., et al. (2016). Bacteria dominate the short-term assimilation of plant-derived N in soil. Soil Biol. Biochem. 96, 30–38. doi: 10.1016/j.soilbio.2016.01.009
Starr, E. P., Nuccio, E. E., Pett-Ridge, J., Banfield, J. F., and Firestone, M. K. (2019). Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil. Proc. Natl. Acad. Sci. U.S.A. 116, 25900–25908. doi: 10.1073/pnas.1908291116
Steenbakkers, P. J. M., Li, X.-L., Ximenes, E. A., Arts, J. G., Chen, H., Ljungdahl, L. G., et al. (2001). Noncatalytic docking domains of cellulosomes of anaerobic fungi. J. Bacteriol. 183, 5325–5333. doi: 10.1128/JB.183.18.5325-5333.2001
Summerscales, J., Dissanayake, N. P. J., Virk, A. S., and Hall, W. (2010). A review of bast fibres and their composites. Part 1 - Fibres as reinforcements. Appl. Sci. Manufact. 41, 1329–1335. doi: 10.1016/j.compositesa.2010.06.001
Tamburini, E., Leon, A. G., Perito, B., and Mastromei, G. (2003). Characterization of bacterial pectinolytic strains involved in the water retting process. Environmental Microbiology 5, 730–736. doi: 10.1046/j.1462-2920.2003.00462.x
Tanca, A., Palomba, A., Deligios, M., Cubeddu, T., Fraumene, C., Biosa, G., et al. (2013). Evaluating the impact of different sequence databases on metaproteome analysis: insights from a lab-assembled microbial mixture. PLoS One 8:e82981. doi: 10.1371/journal.pone.0082981
Tanco, S., Gevaert, K., and van Damme, P. (2015). C-terminomics: targeted analysis of natural and posttranslationally modified protein and peptide C-termini. Proteomics 15, 903–914. doi: 10.1002/pmic.201400301
Tedersoo, L., Drenkhan, R., Anslan, S., Morales-Rodriguez, C., and Cleary, M. (2019). High-throughput identification and diagnostics of pathogens and pests: overview and practical recommendations. Mol. Ecol. Resour. 19, 47–76. doi: 10.1111/1755-0998.12959
van Bakel, H., Stout, J. M., Cote, A. G., Tallon, C. M., Sharpe, A. G., Hughes, T. R., et al. (2011). The draft genome and transcriptome of Cannabis sativa. Genome Biol. 12:R102. doi: 10.1186/gb-2011-12-10-r102
vanden Wymelenberg, A., Gaskell, J., Mozuch, M., Sabat, G., Ralph, J., Skyba, O., et al. (2010). Comparative transcriptome and secretome analysis of wood decay fungi Postia placenta and Phanerochaete chrysosporium. Appl. Environ. Microbiol. 76, 3599–3610. doi: 10.1128/AEM.00058-10
Visi, D. K., D’Souza, N., Ayre, B. G., Webber, I. I. I. 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. Indust. Microbiol. Biotechnol. 40, 465–475. doi: 10.1007/s10295-013-1242-1
Wang, J. P., Matthews, M. L., Williams, C. M., Shi, R., Yang, C., Tunlaya-Anukit, S., et al. (2018). Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis. Nat. Commun. 9:1579. doi: 10.1038/s41467-018-03863-z
Wang, W., Chen, D., Sun, X., Zhang, Q., Koide, R. T., Insam, H., et al. (2019). Impacts of mixed litter on the structure and functional pathway of microbial community in litter decomposition. Appl. Soil Ecol. 144, 72–82. doi: 10.1016/j.apsoil.2019.07.006
Wang, Z., Hobson, N., Galindo, L., Zhu, S., Shi, D., McDill, J., et al. (2012). The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 72, 461–473. doi: 10.1111/j.1365-313X.2012.05093.x
Westermann, B., Jacome, A. S. V., Rompais, M., Carapito, C., and Schaeffer-Reiss, C. (2017). Doublet N-terminal oriented proteomics for n-terminomics and proteolytic processing identification. Methods Mol. Biol. 1574, 77–90. doi: 10.1007/978-1-4939-6850-3_6
Williams, M. A., Taylor, E. B., and Mula, H. P. (2010). Metaproteomic characterization of a soil microbial community following carbon amendment. Soil Biol. Biochem. 42, 1148–1156. doi: 10.1016/j.soilbio.2010.03.021
Wilmes, P., and Bond, P. L. (2004). The application of two-dimensional polyacrylamide gel electrophoresis and downstream analyses to a mixed community of prokaryotic microorganisms. Environ. Microbiol. 6, 911–920. doi: 10.1111/j.1462-2920.2004.00687.x
Wróbel-Kwiatkowska, M., Starzycki, M., Zebrowski, J., Oszmiański, J., and Szopa, J. (2007). Lignin deficiency in transgenic flax resulted in plants with improved mechanical properties. J. Biotechnol. 128, 919–934. doi: 10.1016/j.jbiotec.2006.12.030
Wu, L., Wang, H., Zhang, Z., Lin, R., Zhang, Z., and Lin, W. (2011). Comparative metaproteomic analysis on consecutively Rehmannia glutinosa-monocultured rhizosphere soil. PLoS One 6:e20611. doi: 10.1371/journal.pone.0020611
Wu, Y., Li, C. H., Zhao, J., Xiao, Y. L., and Cao, H. (2015). Metaproteome of the microbial community in paddy soil after long-term treatment with mineral and organic fertilizers. Isr. J. Ecol. Evol. 61, 146–156. doi: 10.1080/15659801.2015.1137438
Yu, Y., Huang, W., Chen, H., Wu, G., Yuan, H., Song, X., et al. (2014). Identification of differentially expressed genes in flax (Linum usitatissimum L.) under saline-alkaline stress by digital gene expression. Gene 549, 113–122. doi: 10.1016/j.gene.2014.07.053
Yuan, D., Tang, Z., Wang, M., Gao, W., Tu, L., Jin, X., et al. (2016). The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres. Sci. Rep. 5:17662. doi: 10.1038/srep17662
Zhang, N., and Deyholos, M. K. (2016). RNASeq analysis of the shoot apex of flax (Linum usitatissimum) to identify phloem fiber specification genes. Front. Plant Sci. 7:950. doi: 10.3389/fpls.2016.00950
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.-R. (2013). Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 14:274. doi: 10.1186/1471-2164-14-274
Zhao, Z., Liu, H., Wang, C., and Xu, J.-R. (2014). Correction: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 15:6. doi: 10.1186/1471-2164-15-6
Keywords: cell wall, flax, soil, meta-omics, microbiota, holobiont, natural fibers, retting
Citation: Djemiel C, Goulas E, Badalato N, Chabbert B, Hawkins S and Grec S (2020) Targeted Metagenomics of Retting in Flax: The Beginning of the Quest to Harness the Secret Powers of the Microbiota. Front. Genet. 11:581664. doi: 10.3389/fgene.2020.581664
Received: 09 July 2020; Accepted: 21 September 2020;
Published: 27 October 2020.
Edited by:Gea Guerriero, Luxembourg Institute of Science and Technology, Luxembourg
Reviewed by:Michael Deyholos, University of British Columbia, Canada
Wenxiang Ping, Heilongjiang University, China
Copyright © 2020 Djemiel, Goulas, Badalato, Chabbert, Hawkins and Grec. 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.
*Correspondence: Sébastien Grec, firstname.lastname@example.org