Skip to main content

Inactivation of LACCASE8 and LACCASE5 genes in Brachypodium distachyon leads to severe decrease in lignin content and high increase in saccharification yield without impacting plant integrity



Dedicated lignocellulosic feedstock from grass crops for biofuel production is extensively increasing. However, the access to fermentable cell wall sugars by carbohydrate degrading enzymes is impeded by lignins. These complex polymers are made from reactive oxidized monolignols in the cell wall. Little is known about the laccase-mediated oxidation of monolignols in grasses, and inactivation of the monolignol polymerization mechanism might be a strategy to increase the yield of fermentable sugars.


LACCASE5 and LACCASE8 are inactivated in a Brachypodium double mutant. Relative to the wild type, the lignin content of extract-free mature culms is decreased by 20–30% and the saccharification yield is increased by 140%. Release of ferulic acid by mild alkaline hydrolysis is also 2.5-fold higher. Interfascicular fibers are mainly affected while integrity of vascular bundles is not impaired. Interestingly, there is no drastic impact of the double mutation on plant growth.


This work shows that two Brachypodium laccases with clearly identified orthologs in crops are involved in lignification of this model plant. Lignification in interfascicular fibers and metaxylem cells is partly uncoupled in Brachypodium. Orthologs of these laccases are promising targets for improving grass feedstock for cellulosic biofuel production.


The optimum utilization of plant biomass as an alternative to fossil carbon in order to produce biomolecules, bioenergy and biomaterials is a major challenge for the next decade [1]. Lignocellulose biomass is mainly composed of hemicellulose and cellulose polysaccharides, and phenolic polymers called lignins. It is well established that lignins negatively affect the cellulose-to-ethanol conversion process [2], but lignins also represent potential sources of high-value products [1]. As grasses, such as Miscanthus spp. and switchgrass (Panicum virgatum) are considered promising bioenergy crops, there is increasing interest in their lignin composition and deposition. Lignins are polymerized by the oxidative coupling of p-coumaryl, coniferyl and sinapyl alcohols (monolignols), which give rise to p-hydroxyphenyl (H), guaiacyl (G) and syringyl (S) units within the lignin polymer, respectively. In addition, lignins formed in grasses (Poaceae) have additional features, the main ones being their association to ferulic acid (FA) and p-coumaric acid (CA) [3,4,5,6,7,8].

The polymerization of monolignols is facilitated by peroxidase and laccase proteins [9]. In planta, the role of laccases in lignification was first revealed in Arabidopsis and poplar plants [10,11,12]. Less is known about laccases in monocotyledonous plants [13, 14], although recent mutant analyses in Brachypodium distachyon (Brachypodium) have proposed roles in lignification for LACCASE5 (LAC5, Bradi1g66720) and potentially LACCASE6 (LAC6, Bradi1g74320) [15]. While the lignin contents of the mature culms from wild type (WT) and lac6 mutant were similar, the lac5 mutant had a lower lignin content (~ 10% lower) [15]. Considering the limited reductions in the lac5 mutant, we continued to search for additional components that facilitate lignin formation in the model Brachypodium plant. The present study identifies and characterizes the function of LACCASE8 (LAC8, Bradi2g23370), which has an important function in the lignification of cell walls in Brachypodium distachyon stems.

Results and discussion

LAC5 and LAC8 are close paralogs with similar expression pattern

The most recent revision (v3.1) of the Brachypodium accession Bd21 genome [16, 17] induced changes in position of some LACCASE genes in the genome (Additional file 1) compared to earlier published data [15]. To update our data set and identify new genes involved in lignification, we conducted phylogenetic analysis of the entire Brachypodium laccase family. We used the Blastp tool from Phytozome ( to identify homologs of LAC5 (Bradi1g66720) in Zea mays, Oryza sativa and Setaria viridis (Additional file 2). We then selected proteins with the highest e-score (> 50%) and compared them with laccases published in [15]. We reconstructed a phylogeny tree with the maximum likelihood comparison method (Additional file 3). Despite some names changed for a few laccases (Additional file 1), no major alteration in clusters were found compared with our previous phylogeny analysis [15]. Consistently, LAC5 was most similar to and clustered together with LAC8 and LACCASE12 (LAC12, Bradi2g54740) protein sequences. Considering that the lignin content of the lac5 mutant was found to be only moderately reduced relative to the WT level [15], we hypothesized that another close paralog might be involved in lignification and that its activity might compensate for LAC5 deficiency. We thus analyzed in silico expression levels of the two closest paralogs LAC8 and LAC12 with the Gene Atlas tool from PlaNet ( and the BAR Brachypodium eFP Browser ( [18]. LAC8 but not LAC12 showed a similar expression pattern with LAC5 (, Fig. 1a). In line with this result, signal values relative to LAC8 and LAC5 transcript levels were both comparable and much higher than those of LAC12 in lignified tissues (Fig. 1a). We verified the expression pattern of LAC8 by qRT-PCR and confirmed data obtained in silico with BAR Brachypodium eFP Browser (Fig. 1b). Indeed, LAC8 is highly expressed in all lignified tissues of Brachypodium (internode, node, peduncle) but poorly detected in developing leaves and young spikelets (Fig. 1b). At last, it is interesting to note that LAC6, LAC8, LAC10 (Bradi2g54680) but not LAC12 belongs the co-expressed network of LAC5 (Additional file 4). Insofar as LAC8 is a close paralog of LAC5 with a highly similar expression pattern, we hypothesized that LAC8 was most probably involved in cell wall lignification.

Fig. 1
figure 1

Expression levels of LAC5, LAC8 and LAC12. a In silico expression levels of the three laccases obtained with the Gene Atlas tool from BAR Brachypodium eFP Browser ( [18]. b qPCR analysis shows the relative expression levels of LACCASE8 transcript in different tissues

Production and characterization of a lac5 lac8 double mutant

We used the TILLING (Targeting-Induced Local Lesions in Genomes) facility available for Brachypodium at University of Saclay, France (, to identify lac8 mutants from a sodium azide mutant collection [19]. Twenty independent lines with mutations in the screened genomic region were isolated using this approach (Fig. 2a). SIFT software analyses revealed which substitutions might partially or totally disrupt the protein activity (Fig. 2a) [19]. By so doing, we selected one homozygous line (line 5731) which harbored an induced stop codon mutation for further analysis. The premature stop codon in this line is predicted to produce a truncated LAC8 protein that lacks the last 29 amino acids of the highly conserved C-terminal region of the LACCASE protein family (Fig. 2b).

Fig. 2
figure 2

Allelic series of mutations in LAC8 gene (Bradi2g23370) identified by TILLING and peptide alignment of different Brachypodium laccase C-terminals. a Features of the mutagenized lines identified by TILLING. Lines harboring mutation in non-translated region (intron, 5′ or 3′UTR) were discarded from the table. b Alignment of amino acids present in the C-terminal region of all Brachypodium laccases present in the genome of Brachypodium distachyon Bd21-3 accession. The codon encoding tryptophan at position 551 is replaced by a stop codon in the mutant line Bd5731. Colors illustrate the conserved regions between proteins

To bypass a possible functional genetic redundancy with LAC5 as suggested by both expression patterns of LAC5 and LAC8, we crossed lac85731 homozygous line isolated in this study to the previously characterized lac54442 homozygous line [15]. We genotyped 136 F2 plants from the resulting F1 generation using high-resolution melting (HRM) qPCR. Frequency of both homozygous mutant alleles in F2 progeny revealed a typical 1:16 Mendelian-like segregation. In parallel with HRM qPCR analysis, we generated vibratome cross sections of dried stems from each genotype and performed phloroglucinol–HCl staining. This staining results in the formation of red stain in the presence of lignified cell walls [20]. Of all genotypes, it was rather obvious that the homozygous lac5 lac8 double mutant (LAC5LAC8, Fig. 3) showed the most substantial difference with WT (LAC5++LAC8++, Fig. 3) while this line did not show a severe growth phenotype when cultivated in the greenhouse or in the growth chamber (Fig. 4). UV excitation fluorescence microscopy has also been used as a sensitive method to document the spatial distribution of lignin in cell walls [21, 22]. Using two-photon UV fluorescence microscopy, we clearly observed that fluorescence at 420–460 nm was impacted in the secondary cell wall of interfascicular fibers (inFF) and mestome (mes) cells in 30-day-old plants (Fig. 5a–d). Surprisingly, fluorescence remained intense inside the fascicular bundle of the double mutant (Fig. 5b). This result suggests that lignins (which emit at 420–460 nm) from protoxylem (pXyl), metaxylem (mXyl) and intrafascicular fibers (FF) are not or moderately impacted. Dried senesced stems (90-day-old plants) treated with phloroglucinol–HCl confirmed these results (Fig. 5e, f). More importantly, while phloroglucinol–HCl staining was almost undetectable in the secondary cell wall of interfascicular fibers of lac5 lac8, it was still observed in primary cell wall (Fig. 5g). This result is consistent with the observation that some laccases direct lignification in the secondary cell wall [23] but were not localized to the cell corners/middle lamella in Arabidopsis [24]. In lac5 lac8 dried stems, the cell walls of interfascicular fibers were observed to be jagged and malformed probably because of the lack of cell strength and rigidity normally imparted by lignification in WT (Fig. 5h, i). We hypothesize that the lignin polymer functions as a ‘cement’ in the cell walls of WT plants and in the absence of this cement in the double mutant, the cell wall polysaccharides collapsed in the inner cell space when the plant is dried and cut.

Fig. 3
figure 3

Lignin patterning in different genotyped mutants. Cross sections of dried stems from F2 plants subjected to phloroglucinol–HCl are imaged. Each genotype respective to mutation in LACCASE5 or LACCASE8 is indicated. ++: WT homozygous, +−: WT heterozygous, −−: mutant homozygous

Fig. 4
figure 4

Growth phenotype of WT and of lac5 lac8 mutant

Fig. 5
figure 5

Two-photon fluorescence microscopy imaging and phloroglucinol–HCl staining of WT and lac5 lac8 lignified tissues. Cross sections of 30-day-old plants were imaged using two-photon fluorescence microscopy (green: 420–460-nm emission, purple: 495–540-nm emission ad) or stained with phloroglucinol–HCl prior imaging under visible microscopy (ei). WT: a, c, e, h; lac5 lac8: b, d, f, g, i. FF intrafascicular fibers, pXy protoxylem, mXy metaxylem, inFF interfascicular fibers, epi epidermis, mes mestome. Blue arrows show red staining of primary cell wall between interfascicular fiber cells in the double mutant

In addition to visible microscopy and two-photon UV fluorescence microscopy, we used Raman spectroscopy to characterize changes in the cell wall composition of interfascicular fibers and metaxylem cells (Fig. 6a). We first compared Raman spectra from mutant versus WT interfascicular fiber cells (Fig. 6b). After baseline correction and normalization on standard error, we noticed a substantial decrease in signal at the bands 1171, 1203, 1268 and 1337 cm−1 in the double mutant. Interestingly, these bands were attributed to lignin polymers in previous papers. Indeed, based on Raman measurements on standard lignins, 1268 cm−1 has been shown to be specific of G-unit lignin type and 1337 cm−1 has been assigned to both G and S units [25,26,27,28]. The 1203 cm−1 band was assigned to OCH3 present in both G and S units [26, 29]. The 1171 cm−1 band has been assigned to lignins and to esters (C–O–C) [30]. Bands at 1603 and 1630 cm−1 were also found slightly changed between mutant and WT. Overall, these signals are often attributed to aromatic compounds and thus to lignin polymers, ferulic and coumaric acids [31, 32].

Fig. 6
figure 6

Raman spectroscopy of interfascicular fiber cell wall and metaxylem cell wall from WT and lac5 lac8. a Scanning electron microscopy of lignified tissues in WT. Yellow arrows illustrate impacts of the laser on cell wall of interfascicular fiber cell and metaxylem. b Raman spectral acquisition of interfascicular fiber cell wall in mutant (red line) and WT (black line). c Raman spectral acquisition of metaxylem cell wall in mutant (red line) and WT (black line). inFF interfascicular fiber, mXy metaxylem. Wavelength numbers for major peaks are indicated

In the Raman spectroscopy, one new peak at 1370 cm−1 was detected for the lac5 lac8 mutant, and signals in the vicinity of 1424 and 1460 cm−1 were also found increased. These bands are attributable to C–C and C–H, and in consequence, it is difficult to assign them to a specific compound in the cell wall, even if they have been often attributed to lignin, cellulose and hemicellulose [30, 32, 33]. Interestingly, when we compared Raman spectra from metaxylem cell walls of WT versus those of double mutant (Fig. 6c), there was no detectable change in intensity of lignin-related peaks (1171, 1203, 1268 and 1337 cm−1) despite the fact that we found again a slight increase in signal in the region of 1330–1460 cm−1. These results confirm that lignification is specifically affected in interfascicular fiber cell wall, but it is not drastically affected in metaxylem cells as suggested by phloroglucinol–HCl staining and two-photon confocal analysis.

Consistent with the lignin being unaffected in the vessel secondary cell walls, we did not observe collapsed xylem vessels in the Brachypodium lac5 lac8 mutant (Fig. 5g), unlike the Arabidopsis lac4 lac17 double mutant [10] or in Arabidopsis lac4 lac17 lac11 triple mutants [11].

Lignin content and structure are strongly affected in the double lac5 lac8 mutant

Lignin analyses were performed on single and double mutants along with the WT. Both the gravimetric Klason method and the spectrometric acetyl bromide method confirmed that LAC5 deficiency induced a moderate decrease (about 10%) in the lignin level of Brachypodium mature culms, as compared to WT level (Table 1). A similar decrease (about 10%) was observed in the lac8 sample. By contrast, in the lac5 lac8 double mutant, the lignin level was found to be more severely reduced relative to the WT (20 to 30%, Table 1, [15]). This reduction provides the clue that, together with LAC5, the LAC8 gene is involved in the lignification of Brachypodium stem.

Table 1 Lignin content of extract-free mature culms from Brachypodium wild type and laccase mutants

Lignin structure was studied by thioacidolysis. This lignin-specific method provides p-hydroxyphenyl (H), guaiacyl (G) and syringyl (S) thioethylated monomers from H, G and S lignin units that are only bound together by β-O-4 bonds [34]. When expressed relative to the Klason lignin content, thioacidolysis yields were very similar in all samples, which indicates that the proportion of the target β-O-4 bonds is relatively consistent among the various lines (Table 2). Whatever the sample, H monomers were recovered as minor components. Relative to the control value, the molar frequency of G monomers was found to be decreased in all laccase mutants, whereas that of S monomers was concomitantly increased (Table 2). Not unexpectedly, such opposite variations induced dramatic changes in the S/G molar ratio from 1.85, for the control, up to 2.90 in lac5 lac8 mutant. These variations reveal that the deposition of G lignin units is more strongly affected in the lac5- and/or lac8-deficient mutants. This slight but systematic decrease in G lignin units, already observed in Brachypodium lac5 and Arabidopsis lac17 mutants [10, 15], supports the hypothesis of a spatial and/or temporal specificity of laccase for coniferyl alcohol. These laccases could be more active during the early lignification stages when G units are specifically deposited prior to S units [35].

Table 2 Determination of the main H, G and S monomers released by thioacidolysis of extract-free mature culms from Brachypodium wild type and laccase mutants

The amount of CA and FA ester-linked to Brachypodium cell walls was measured by mild alkaline hydrolysis. It is now well established that, in grass lignified tissues, most CA esters are linked to lignins while a lower amount is ester-linked to the arabinose substituents of arabinoxylans [36]. The amount of lignin-associated CA esters was reduced in the double mutant (Table 3), a phenomenon that could be related to its lower lignin content. Cell wall FA esters are predominantly linked to the arabinose substituents of arabinoxylans. These FA esters are proposed to act as initiation sites for lignification as first revealed by the identification of FA ether-linked to a G unit [8]. These oxidatively driven and resistant linkages between FA esters and lignins reduce the FA recovery after mild alkaline hydrolysis. On this basis, the dramatic increase in alkali-releasable FA from the double mutant may be directly related to its substantially lower lignin content.

Table 3 Determination of p-coumaric acid (CA) and ferulic acid (FA) released by mild alkaline hydrolysis of extract-free mature culms from Brachypodium wild type and laccase mutants

Because the drastic decrease in lignin content may deeply impact the cell wall’s susceptibility to enzymatic digestion, we performed saccharification assays without pre-treatment. Saccharification efficiency was evaluated both by the weight loss induced by the cellulose hydrolysis treatment and by the amount of glucose released in the reaction medium. As revealed in Table 4, saccharification efficiency was slightly increased in the single lac5 mutant, whereas it was not changed in the lac8 sample. By contrast and consistent with the severely reduced lignin content, this saccharification efficiency was substantially increased from the biomass of the lac5 lac8 double mutant (Table 4).

Table 4 Saccharification assays of extractive-free mature stem from Brachypodium wild type and laccase mutants


This work provides the first demonstration that LAC8 gene is involved in the lignification of Brachypodium culm. The impact of loss of function of LAC5 and LAC8 was specifically on the interfascicular fiber cells, which make up the bulk of the grass stem. This was reflected in the substantial increased saccharification yield of the double mutant lines. Cytology studies and Raman spectroscopy showed that lignification of metaxylem cells is not or only slightly impacted in lac5 lac8 double mutant, which permits these double mutants to continue to grow relatively normally. This result suggests that vascular tissues produce other enzymes that balance the loss of activity of LAC5 and LAC8. At last, close orthologs of the Brachypodium LAC5 and LAC8 genes are clearly identifiable in the genome of maize or rice, which confirms that these laccase genes are both promising targets to genetically redesign grass cell walls with improved susceptibility to enzymatic hydrolysis.

Materials and methods

Plant material and growth conditions

Single mutants were selected by TILLING in the collection of chemically induced Brachypodium mutants as described in [19]. lac5 (line 4442) and lac8 (line 5731) were generated by selfing the initial tilled lines. Brachypodium plants (accession Bd21-3) were grown in a greenhouse under long-day conditions (18-h light, 400-W sodium lamps). Day and night temperatures were 23 °C and 18 °C, respectively.

Histochemical staining

All the histochemical staining was performed on sections cut in the middle of the second internode from the top. Samples were hand-sectioned or embedded in 7% agarose before being transversely sectioned at a thickness of 50 µm using a vibratome (Leica VT1000S, Leica, Germany). Lignin deposition and composition were investigated histochemically by Wiesner staining (Phloroglucinol–HCl) as previously published [20]. All sections were observed under a Zeiss AxioPlan 2 microscope system with automatic exposure times.

TILLING and sequence analysis

PCR amplification and detection of mutations DNA amplification are based on nested PCR. The first PCR amplification is a standard PCR with target-specific primers and 10 ng of Brachypodium genomic DNA. Forward primer LCZ23_N1F1: CGCAGTCGCCAACCACACGCTGACGG and reverse primer LCZ23_N1R1: CGTGCAACACCATCGACCCATTCAGC were used. One milliliter of the first PCR product served as a template for the second nested PCR amplification, with a combination of specific primers carrying M13 tail and M13 universal primers, M13F700 (59-CACGACGTTGTAAAACGAC-39) and M13R800 (59-GGATAACATTTCACACAGG-39), labeled at the 5′ end with infrared dyes IRD700 and IRD800 (LI-CORH, Lincoln, Nebraska, USA), respectively. Mutation detection was carried out as described previously except for the second PCR. This PCR was carried out using 0.05 mM of specific primers carrying M13 tail and 0.1 mM of M13 universal primers. The identity of the mutations was determined by sequencing. Sequence analysis was performed with CODDLE software (Codons Optimized to Discover Deleterious Lesions) [37] and PARSESNP software [38]. Prediction of the impact of each mutation was made with SIFT software [39] as described in [19]. Multiple sequence alignment of full-length protein sequences was performed with ClustalW software (

Lignin content and structure determination

All main stems of each plant were collected and ground after removing spikelets and leaves. Ground samples were sequentially extracted at 60 °C with 50 mL of ethanol, water and ethanol. At each step, the samples were vortexed. These steps were repeated twice before sample drying. The extracted and dried samples, referred to as extract-free samples, were used for lignin analyses.

Lignin content was measured by the Klason method and the acetyl bromide method according to [40]. Lignin structure was studied by thioacidolysis, as previously described [41]. The lignin-derived thioacidolysis monomers were identified by gas chromatography–mass spectrometry as their trimethylsilylated derivatives. All the analyses were performed with at least three biological replicates. Significant differences were inferred by one-way ANOVA (Tukey’s HSD, P < 0.05).

PCR high-resolution melting (HRM) and quantitative PCR

Approximately 20 ng of leaf DNA extract was subjected to PCR amplification on Eppendorf Realplex2 Mastercycler using the SYBR Green kit (Bio-Rad) and using the following conditions: 95 °C for 5 min, followed by 45 cycles of 95 °C for 30 s, 62 °C for 30 s and 72 °C for 30 s. Amplifications were made with the following primers: LAC5fw-CCGGAGGTTGGGTCGCCATCAGGTT and LAC5rev-CATTTACGGTTAAGCAAGAACGTGTGCACG for LACCASE5 (Bradi1g66320); LAC8fw-GGGCATGCAGGTGTATGGAT and LAC8rev-CGGCAACTTCTGGTTCGGTA for LACCASE8 (Bradi2g23370). A high melting curve program was used to detect the different alleles with the following program: 95 °C for 10 s, 65 °C for 5 s and increase of 0.2 °C every 5 s until 95 °C. Fluorescence was quantified every 5 s. Relative expression of LAC8 was quantified as described in [15] using the following primers: qLAC8-fw-TACACGTTCAATGTGACAATGGCG and qLAC8-rev CTCACGCCGTGCCAGTGGAAG.

Phylogeny tree

An unrooted phylogenetic tree was reconstructed with PhyML in the platform ( [42]. Sequences were aligned with MUSCLE (V3.8.31) and then corrected with Gblocks (V0.91b). The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.1/3.0 aLRT) and graphically edited using TreeDyn (v198.3). Protein homologs from Brachypodium, rice, maize and Setaria were identified and downloaded from Phytozome V.12.1 (

Cell wall saccharification

Saccharification assays were performed as described by [10]. For each sample, 30 mg of extract-free samples was incubated with 4 mL of acetate buffer pH 4.5, containing 4 mg/mL commercial cellulase (cellulase Onozuka-R10; Serva) and 0.5 mg/mL NaN3 for 3 days at 45 °C on a carousel. After centrifugation, the glucose content of the supernatant was determined via enzymatic assay with the bioMérieux Kit (bioMérieux, Craponne, France). The pellet was washed twice with water, then freeze-dried and weighed to evaluate the weight loss.

Raman spectroscopy

Forty-micrometer sections of internodes were cut using a vibratome (HM 650 V, Microm Microtech), then placed with a drop of water on an aluminum slide and sealed with a coverslip and nail polish. Raman spectra were recorded using a microconfocal Raman inVia Reflex (Renishaw), equipped with a double-edge filter to eliminate the Rayleigh scattering with a cut at 100 cm−1, three wavelengths (532 nm, 633 nm and 785 nm) and two used different gratings (1800 and 1200 L/mm). The setup consisted of a confocal microscope aligned with an automated XYZ table, where the displacement motors generated 100 nm steps. Then, the focused power of the laser beam was checked for each wavelength to avoid any transformation or heating of the samples. The spectral resolution was smaller than 4 cm−1 with a precision better than 0.5 cm−1. Several Raman measurements were taken in secondary cell walls of metaxylem cell and sclerenchyma. Three to five biological repetitions per genotype were analyzed. A total of 20–30 measurements were taken for each cell type (metaxylem cell or sclerenchyma). Cosmic ray removal was applied using 4.2Wire software (Renishaw), and spectra were further analyzed using OriginLab Pro software. An average spectrum was obtained for each region of interest, which was baseline-corrected and normalized.

Scanning electron microscopy (SEM)

Two-millimeter hand-cross sections of 30-day-old WT stems were spray-coated with cobalt using a Cressington 208HR High Resolution Sputter Coater and observed with a Hitachi S-4700 field emission SEM.

Two-photon fluorescence microscopy

An Olympus FV1000 Multiphoton Laser Scanning Microscope with a tunable MaiTai BB DeepSee (710–990 nm) laser was used to stimulate lignin auto-fluorescence. Excitation (350 to 370 nm) of lignin auto-fluorescence is achieved by the absorption of two 730-nm photons [43]. During imaging, two emission channels (420–460 nm and 495–540 nm) are simultaneously collected. Samples were imaged using the Olympus XL Plan N25X objective.

Availability of data and materials

All data generated or analyzed during this study are included in the published article and its Additional files. Materials are available from the corresponding authors.


  1. Ragauskas AJ, Williams CK, Davison BH, Britovsek G, Cairney J, Eckert CA, et al. The path forward for biofuels and biomaterials. Science. 2006;311(5760):484–9.

    Article  CAS  Google Scholar 

  2. Wang Y, Chantreau M, Sibout R, Hawkins S. Plant cell wall lignification and monolignol metabolism. Front Plant Sci. 2013;4:220.

    PubMed  PubMed Central  Google Scholar 

  3. Ralph J, Hatfield RD, Quideau S, Helm RF, Grabber JH, Jung HJG. Pathway of p-coumaric acid incorporation into maize lignin as revealed by NMR. J Am Chem Soc. 1994;116(21):9448–56.

    Article  CAS  Google Scholar 

  4. Ralph JQS, Grabber JH, Hatfield RD. Identification and synthesis of new ferulic acid dehydrodimers present in grass cell walls. J Chem Soc Perkin Trans. 1994;1:3485–98.

    Article  Google Scholar 

  5. Grabber JH, Quideau S, Ralph J. p-coumaroylated syringyl units in maize lignin: implications for beta-ether cleavage by thioacidolysis. Phytochemistry. 1996;43(6):1189–94.

    Article  CAS  Google Scholar 

  6. Ralph J. Hydroxycinnamates in lignification. Phytochem Rev. 2010;9(1):65–83.

    Article  CAS  Google Scholar 

  7. Karlen SD, Zhang C, Peck ML, Smith RA, Padmakshan D, Helmich KE, et al. Monolignol ferulate conjugates are naturally incorporated into plant lignins. Sci Adv. 2016;2(10):e1600393.

    Article  Google Scholar 

  8. Jacquet G, Pollet B, Lapierre C, Mhamdi F, Rolando C. New ether-linked ferulic acid-coniferyl alcohol dimers identified in grass cell walls. J Agric Food Chem. 1995;43:2746–51.

    Article  CAS  Google Scholar 

  9. Vanholme R, Morreel K, Darrah C, Oyarce P, Grabber JH, Ralph J, et al. Metabolic engineering of novel lignin in biomass crops. New Phytol. 2012;196(4):978–1000.

    Article  CAS  Google Scholar 

  10. Berthet S, Demont-Caulet N, Pollet B, Bidzinski P, Cézard L, Le Bris P, et al. Disruption of LACCASE4 and 17 results in tissue-specific alterations to lignification of Arabidopsis thaliana stems. Plant Cell. 2011;23(3):1124–37.

    Article  CAS  Google Scholar 

  11. Zhao Q, Nakashima J, Chen F, Yin Y, Fu C, Yun J, et al. Laccase is necessary and nonredundant with peroxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell. 2013;25(10):3976–87.

    Article  CAS  Google Scholar 

  12. Lu S, Li Q, Wei H, Chang MJ, Tunlaya-Anukit S, Kim H, et al. Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Proc Natl Acad Sci USA. 2013;110(26):10848–53.

    Article  CAS  Google Scholar 

  13. Cesarino I, Araujo P, Sampaio Mayer JL, Vicentini R, Berthet S, Demedts B, et al. Expression of SofLAC, a new laccase in sugarcane, restores lignin content but not S:G ratio of Arabidopsis lac17 mutant. J Exp Bot. 2013;64(6):1769–81.

    Article  CAS  Google Scholar 

  14. Wang J, Feng J, Jia W, Fan P, Bao H, Li S, et al. Genome-wide identification of sorghum bicolor laccases reveals potential targets for lignin modification. Front Plant Sci. 2017;8:714.

    Article  Google Scholar 

  15. Wang Y, Bouchabke-Coussa O, Le Bris P, Antelme S, Soulhat C, Gineau E, et al. LACCASE 5 is required for lignification of the Brachypodium distachyon culm. Plant Physiol. 2015;168(1):192–204.

    Article  CAS  Google Scholar 

  16. Gordon SP, Contreras-Moreira B, Woods DP, Des Marais DL, Burgess D, Shu S, et al. Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure. Nat Commun. 2017;8(1):2184.

    Article  Google Scholar 

  17. Initiative TIB. Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010;463(7282):763–8.

    Article  Google Scholar 

  18. Sibout R, Proost S, Hansen BO, Vaid N, Giorgi FM, Ho-Yue-Kuang S, et al. Expression atlas and comparative coexpression network analyses reveal important genes involved in the formation of lignified cell wall in Brachypodium distachyon. New Phytol. 2017;215(3):1009–25.

    Article  CAS  Google Scholar 

  19. Dalmais M, Antelme S, Ho-Yue-Kuang S, Wang Y, Darracq O, Bouvier d’Yvoire M, et al. A TILLING platform for functional genomics in Brachypodium distachyon. PLoS ONE. 2013;8(6):e65503.

    Article  CAS  Google Scholar 

  20. Bouvier d’Yvoire M, Bouchabke-Coussa O, Voorend W, Antelme S, Cezard L, Legee F, et al. Disrupting the cinnamyl alcohol dehydrogenase 1 gene (BdCAD1) leads to altered lignification and improved saccharification in Brachypodium distachyon. Plant J. 2013;73(3):496–508.

    Article  Google Scholar 

  21. Donaldson LA, Knox JP. Localization of cell wall polysaccharides in normal and compression wood of radiata pine: relationships with lignification and microfibril orientation. Plant Physiol. 2012;158(2):642–53.

    Article  CAS  Google Scholar 

  22. Donaldson L, Vaidya A. Visualising recalcitrance by colocalisation of cellulase, lignin and cellulose in pretreated pine biomass using fluorescence microscopy. Sci Rep. 2017;7:44386.

    Article  Google Scholar 

  23. Schuetz M, Benske A, Smith RA, Watanabe Y, Tobimatsu Y, Ralph J, et al. Laccases direct lignification in the discrete secondary cell wall domains of protoxylem. Plant Physiol. 2014;166(2):798–807.

    Article  Google Scholar 

  24. Yi Chou E, Schuetz M, Hoffmann N, Watanabe Y, Sibout R, Samuels AL. Distribution, mobility, and anchoring of lignin-related oxidative enzymes in Arabidopsis secondary cell walls. J Exp Bot. 2018;69(8):1849–59.

    Article  Google Scholar 

  25. Agarwal U, Terashima N. FT-Raman study of dehydrogenation polymer (DHP) lignins. In:12th ISWPC international symposium on wood and pulping chemistry: Madison, Wisconsin, USA, June 9–12, 2003: proceedings, volume III : poster presentations. Madison, WI : University of Wisconsin-Madison, Dept. of Forest Ecology and Management, 2003. p. 123–6.

  26. Larsen KL, Barsberg S. Theoretical and Raman spectroscopic studies of phenolic lignin model monomers. J Phys Chem B. 2010;114(23):8009–21.

    Article  CAS  Google Scholar 

  27. Sun L, Varanasi P, Yang F, Loque D, Simmons BA, Singh S. Rapid determination of syringyl: guaiacyl ratios using FT-Raman spectroscopy. Biotechnol Bioeng. 2011;109(3):647–56.

    Article  Google Scholar 

  28. Agarwal UP, McSweeny JD, Ralph SA. FT–Raman investigation of milled-wood lignins: softwood, hardwood, and chemically modified black spruce lignins. J Wood Chem Technol. 2011;31(4):324–44.

    Article  CAS  Google Scholar 

  29. Agarwal UP, Kawai N. Self-absorption phenomenon in near-infrared Fourier transform Raman spectroscopy of cellulosic and lignocellulosic materials. Appl Spectrosc. 2005;59(3):385–8.

    Article  CAS  Google Scholar 

  30. Wiley J, Atalla R. Band assignments in the raman spectra of celluloses. Carbohydr Res. 1987;160:113–29.

    Article  CAS  Google Scholar 

  31. Ozparpucu M, Gierlinger N, Burgert I, Van Acker R, Vanholme R, Boerjan W, et al. The effect of altered lignin composition on mechanical properties of cinnamyl alcohol dehydrogenase (CAD) deficient poplars. Planta. 2018;247(4):887–97.

    Article  CAS  Google Scholar 

  32. Sene CFB, McCann MC, Wilson RH, Grinter R. Fourier-transform Raman and Fourier-transform infrared spectroscopy (an investigation of five higher plant cell walls and their components). Plant Physiol. 1994;106(4):1623–31.

    Article  CAS  Google Scholar 

  33. Kačuráková M, Wellner N, Ebringerová A, Hromádková Z, Wilson RH, Belton PS. Characterisation of xylan-type polysaccharides and associated cell wall components by FT-IR and FT-Raman spectroscopies. Food Hydrocolloids. 1999;13(1):35–41.

    Article  Google Scholar 

  34. Lapierre C. Determining lignin structure by chemical degradations. In: Heitner C, Dimmel D, Schmidt JA, editors. Lignin and lignans—advances in chemistry. Boca Raton: CRC Press, Taylor & Francis Group; 2010. p. 11–48.

    Chapter  Google Scholar 

  35. Terashima N, Fukushima K. Comprehensive model of the lignified plant cell wall. In: Jung HG, Buxton DR, Hatfield RD, Ralph R, editors. Forage cell wall structure and digestibility. Madison: American Society of Agronomy; 1993. p. 247–70.

    Google Scholar 

  36. Hatfield RD, Rancour DM, Marita JM. Grass cell walls: a story of cross-linking. Front Plant Sci. 2016;7:2056.

    PubMed  Google Scholar 

  37. Till BJ, Reynolds SH, Greene EA, Codomo CA, Enns LC, Johnson JE, et al. Large-scale discovery of induced point mutations with high-throughput TILLING. Genome Res. 2003;13(3):524–30.

    Article  CAS  Google Scholar 

  38. Taylor NE, Greene EA. PARSESNP: a tool for the analysis of nucleotide polymorphisms. Nucleic Acids Res. 2003;31(13):3808–11.

    Article  CAS  Google Scholar 

  39. Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11(5):863–74.

    Article  CAS  Google Scholar 

  40. Dence CW. The determination of lignin. In: Lin SY, Dence CW, editors. Methods in lignin chemistry. Springer: Berlin; 1992. p. 33–61.

    Chapter  Google Scholar 

  41. Lapierre C, Pollet B, Petit-Conil M, Toval G, Romero J, Pilate G, et al. Structural alterations of lignins in transgenic poplars with depressed cinnamyl alcohol dehydrogenase or caffeic acid O-methyltransferase activity have an opposite impact on the efficiency of industrial kraft pulping. Plant Physiol. 1999;119(1):153–63.

    Article  CAS  Google Scholar 

  42. Dereeper A, Audic S, Claverie JM, Blanc G. BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC Evol Biol. 2010;10:8.

    Article  Google Scholar 

  43. Zipfel WR, Williams RM, Webb WW. Nonlinear magic: multiphoton microscopy in the biosciences. Nat Biotechnol. 2003;21(11):1369–77.

    Article  CAS  Google Scholar 

Download references


We thank Anne-Laure Chateigner-Boutin for reading the manuscript. We also acknowledge Sylvie Durand for interesting discussions about Raman experiments. We thank Sebastien Santini for making available the site SEM picture of Brachypodium stem was made thanks to the UBC Bioimaging Facility


This work was funded by INRA, Agreenskills program (European Union’s 7th Framework Programme for research, technological development and demonstration under Grant Agreement Numbers FP7-267196 and FP7-609398, and INRA) and the Agence Nationale de la Recherche (ANR-14-CE19-0012-01) project. The IJPB benefits from the support of the LabEx Saclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS).

Author information

Authors and Affiliations



RS supervised the study, designed experiments, performed experiments, analyzed data, prepared figures and cowrote the manuscript. CL performed experiments, analyzed data, prepared tables and cowrote the manuscript. LS an MS performed experiments, analyzed data and edited the manuscript. MD and AB performed TILLING experiments. SA cultivated plants and selected genotypes. ADO analyzed Raman data and prepared figures. LC and FL performed lignin analyses. PLB, YW and CB genotyped and characterized mutants. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Catherine Lapierre or Richard Sibout.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1.

The laccase protein family in the sequenced Bd21 natural accession. Laccase proteins published in [15] are listed with their respective name regarding the first version (V1.0) or the most recent version (V3.1) of the Bd21 genome sequence.

Additional file 2.

Predicted LACCASE proteins in Brachypodium, maize, rice and Setaria. Putative laccases protein sequences of Zea mays, Oryza sativa, Setaria viridis and Brachypodium distachyon were obtained from proteomic databases available on Phytozome (

Additional file 3.

Phylogeny tree reconstructed with Brachypodium, maize, rice and Setaria LACCASES. The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program. The putative proteins sequences used to reconstruct the tree are available in Additional files 1 and 2. Branch length is proportional to the number of substitutions per site and represents evolutionary distance as indicated by the scale bar. The colored cluster highlights the positions of LACCASE 5 and 8 and their closest orthologs in other species.

Additional file 4.

Genes co-expressed with LAC5 (LAC8 is in bold red). The table containing genes found in the co-expression network was made with the search tool from PlaNet ( using LAC5 gene (Brad1g66720) as a bait.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Le Bris, P., Wang, Y., Barbereau, C. et al. Inactivation of LACCASE8 and LACCASE5 genes in Brachypodium distachyon leads to severe decrease in lignin content and high increase in saccharification yield without impacting plant integrity. Biotechnol Biofuels 12, 181 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: