Skip to main content

Engineered reduction of S-adenosylmethionine alters lignin in sorghum

Abstract

Background

Lignin is an aromatic polymer deposited in secondary cell walls of higher plants to provide strength, rigidity, and hydrophobicity to vascular tissues. Due to its interconnections with cell wall polysaccharides, lignin plays important roles during plant growth and defense, but also has a negative impact on industrial processes aimed at obtaining monosaccharides from plant biomass. Engineering lignin offers a solution to this issue. For example, previous work showed that heterologous expression of a coliphage S-adenosylmethionine hydrolase (AdoMetase) was an effective approach to reduce lignin in the model plant Arabidopsis. The efficacy of this engineering strategy remains to be evaluated in bioenergy crops.

Results

We studied the impact of expressing AdoMetase on lignin synthesis in sorghum (Sorghum bicolor L. Moench). Lignin content, monomer composition, and size, as well as biomass saccharification efficiency were determined in transgenic sorghum lines. The transcriptome and metabolome were analyzed in stems at three developmental stages. Plant growth and biomass composition was further evaluated under field conditions. Results evidenced that lignin was reduced by 18% in the best transgenic line, presumably due to reduced activity of the S-adenosylmethionine-dependent O-methyltransferases involved in lignin synthesis. The modified sorghum features altered lignin monomer composition and increased lignin molecular weights. The degree of methylation of glucuronic acid on xylan was reduced. These changes enabled a ~20% increase in glucose yield after biomass pretreatment and saccharification compared to wild type. RNA-seq and untargeted metabolomic analyses evidenced some pleiotropic effects associated with AdoMetase expression. The transgenic sorghum showed developmental delay and reduced biomass yields at harvest, especially under field growing conditions.

Conclusions

The expression of AdoMetase represents an effective lignin engineering approach in sorghum. However, considering that this strategy potentially impacts multiple S-adenosylmethionine-dependent methyltransferases, adequate promoters for fine-tuning AdoMetase expression will be needed to mitigate yield penalty.

Background

Dedicated bioenergy crops and residues from agriculture and forestry represent a renewable source of fermentable sugars for the production of advanced bioproducts and biofuels [1, 2]. A large fraction of the sugars contained in lignocellulosic biomass derive from the cell wall polymers cellulose and hemicellulose, which are surrounded by a protective, hydrophobic and resistant polymer called lignin [3]. The main units in lignin arise from the oxidative polymerization of methoxylated hydroxycinnamyl alcohols synthesized from p-coumarate [4]. These include coniferyl and sinapyl alcohols that generate the guaiacyl (G) and syringyl (S) lignin units, respectively. The presence of small amounts of p-hydroxyphenyl alcohol (H units) in lignin has also been documented, in addition to the occurrence of the flavone tricin (T units) in certain plant species (Fig. 1a), as well as a range of other phenolic compounds that form rare and atypical lignin monomers [5,6,7]. Several approaches have been proposed to make the conversion of lignocellulosic feedstocks into bioproducts cheaper and more efficient [8]. For example, the manipulation of lignin content and monomer composition via genetic engineering represents an important milestone to reduce lignin recalcitrance and facilitate the enzymatic deconstruction of cellulose and hemicellulose into simple sugars [9, 10].

Fig. 1
figure 1

Expression of AdoMetase in sorghum. a Simplified lignin biosynthetic pathway showing the involvement of the two AdoMet-utilizing enzymes COMT and CCoAOMT. Dashed arrows denote multiple steps. AdoCys, S-adenosylhomocysteine. b Construct used for sorghum transformation. Genetic elements in the T-DNA include a selectable marker gene for kanamycin resistance (nptII) placed between a duplicated 35S promoter (p2 × 35S) and a terminator (T35S) from the cauliflower mosaic virus, as well as the AdoMetase synthetic gene located between the sorghum COMT gene promoter (pSbCOMT) and the Agrobacterium tumefaciens nopaline synthase gene terminator (TNOS). c AdoMetase expression in four independent transgenic lines. The level of AdoMetase expression compared to the PP2A gene is shown. First-strand cDNA obtained from wild type (WT) were used as negative control. Values are means ± SD of four biological replicates (n = 4 plants). ND not detected

S-adenosylmethionine (AdoMet) is a universal methyl-group donor produced in the cytosol and needed for the methylation of a large number of metabolites as well as a precursor for the synthesis of polyamines, nicotianamine, 5ʹ-deoxyadenosyl radicals, and ethylene [11, 12]. The lignin biosynthetic pathway contains two enzymes that utilize AdoMet for transmethylation reactions: caffeoyl-CoA O-methyltransferase (CCoAOMT) is involved in the 3-O-methylation step that produces feruloyl-CoA, and caffeic acid O-methyltransferase (COMT) methylates 5-hydroxyconiferaldehyde to produce sinapaldehyde. COMT is a promiscuous enzyme that has the ability to catalyze both 5-O- and 3-O-methylation of several lignin pathway intermediates including acids and alcohols, suggesting a functional redundancy between CCoAOMT and COMT in lignin synthesis [13]. In addition, a few studies have demonstrated the role of COMT in the synthesis of tricin [14,15,16] (Fig. 1a). As an important reaction in one-carbon metabolism, AdoMet is regenerated from ATP and methionine (Met) by AdoMet synthetase via the Met cycle. The characterization of several mutants in Arabidopsis, sorghum, and maize highlighted the importance of AdoMet as a methyl donor for lignin biosynthesis. For example, mutation in one of the AdoMet synthetase genes in Arabidopsis results in concomitant reductions of AdoMet and lignin [17]. Moreover, the alteration of genes responsible for the synthesis of 5-methyltetrahydrofolate—a cofactor used by Met synthase for the synthesis of Met and AdoMet regeneration—reduces the lignin content in maize, sorghum, switchgrass, and Arabidopsis [18,19,20,21,22]. In particular, these genetic modifications affect folylpolyglutamate synthase (FPGS) and 5-methyltetrahydrofolate reductase (MTHFR) in the folate cycle that mediates one-carbon metabolism. For example, both the brown midrib4 (bm4) gene in maize and the brown midrib19 (bmr19) gene in sorghum encode an FPGS. The corresponding bm4 and bmr19 mutants show an increased lignin S/G ratio, as also observed in the Arabidopsis fpgs mutant [18, 19, 21, 23]. Silencing of FPGS in switchgrass reduced the lignin content and resulted in a significant biomass yield penalty when grown in the field [20]. A mutation in the brown midrib2 (bm2) gene that encodes a MTHFR in maize also results in a higher lignin S/G ratio [22, 23].

We previously demonstrated in Arabidopsis that the heterologous expression of AdoMet hydrolase from E. coli phage (AdoMetase, EC 3.13.2.2) was an effective approach to reduce AdoMet and lignin [24]. This engineering strategy is potentially transferable to any crop amenable to genetic transformation and could enable the modification of lignin without any prior knowledge of the CCoAOMT/COMT genes involved in lignification or the genes that participate in AdoMet synthesis. By using tissue-specific promoters, this approach also offers an opportunity for reducing AdoMet preferentially in lignifying tissues while avoiding pleiotropic effects resulting from the modification of other metabolic pathways that consume AdoMet. The major demand for one-carbon units is allocated for lignin and polysaccharide biosynthesis, and reducing AdoMet specifically in the tissues that produce lignified secondary cell walls may have limited impact on cellular functions that have lower AdoMet demand [21, 25].

In the current study, we assessed the impact of expressing AdoMetase in the bioenergy crop sorghum (Sorghum bicolor L. Moench) using the promoter of a COMT gene (bmr12) involved in lignification. The transgenic lines have reductions in lignin content, show altered lignin monomer composition, and feature increased lignin molecular weights. These modifications result in improvements of biomass saccharification efficiency compared to the wild type, under both greenhouse and field growth conditions. The analysis of the transcriptome and metabolome in stems of the transgenic lines revealed some changes in metabolic processes other than lignin biosynthesis, indicating some pleiotropic effects associated with AdoMetase expression, which may explain the reduction in biomass yield and developmental delay observed in some lines.

Results

Generation of sorghum expressing AdoMetase

A construct containing the promoter region of the sorghum COMT gene (bmr12) followed by a codon-optimized AdoMetase open reading frame was built for Agrobacterium-mediated sorghum transformation (Fig. 1b). The pSbCOMT promoter was chosen in an attempt to express AdoMetase preferentially in lignifying tissues. Several primary T0 transformants were regenerated from calli and analyzed by Taqman real-time PCR to identify lines containing a single copy of the transgene. Four of these single-copy events were grown in the T1 generation and wild-type (WT) segregants and homozygous plants were identified. T2 seeds from one WT segregant plant from each line were pooled and used as controls for further experiments. The expression of AdoMetase in each transgenic line was confirmed by qPCR using first-strand cDNA synthesized from total RNAs extracted from stems of 3-week-old plants (Fig. 1c). Consistent with AdoMetase activity, the analysis of metabolites extracted from 3-week-old seedlings showed in transgenics some reductions of AdoMet (by 11–26%) and the presence of the two AdoMetase products 5ʹ-methylthioadenosine and homoserine (Fig. 2). Measurements of growth parameters including the number of days to panicle emergence, number of flowering tillers, height of the main tiller, and dry weight from stover at plant maturity revealed that transgenics were shorter and had delayed flowering compared to the WT. In addition, lines AdoMetase #3 and #4 had reduced biomass yield (Table 1).

Fig. 2
figure 2

Adometase enzymatic reaction (a) and quantification of AdoMet, 5ʹ-methylthioadenosine, and homoserine in 3-week-old seedlings from wild-type (WT) and AdoMetase lines (b). Values are means ± SD of five biological replicates (n = 5 plants). Asterisks indicate a significant difference from the WT using the unpaired Student’s t test (*P < 0.05). ND not detected

Table 1 Growth parameters of wild-type and AdoMetase lines

Cell wall composition analysis

The chemical composition of extractive-free cell wall residues obtained from leaf and stem biomass was determined using a two-step acid hydrolysis procedure. The content of acid-insoluble lignin was significantly reduced in the AdoMetase lines #1 (−18%), #2 (−11%), and #4 (−7%). Similar amount of glucose and xylose was obtained from the cell wall of the different lines, indicating that cellulose and xylan contents do not differ between WT and transgenics (Table 2). However, glucuronic acid was significantly increased in the AdoMetase lines #1 and #2, whereas 4-O-methylglucuronic acid was decreased in the AdoMetase lines #2 and #4. Glucuronic acid typically occurs as side chains on the xylan backbone and its degree of methylation (or 4-O-methylglucuronic/glucuronic acid ratio) was significantly reduced in the AdoMetase lines #1, #2, and #4 compared to WT (Table 2). Arabinose, which is also found as side chains on xylan, was not consistently changed among the four AdoMetase lines compared to WT (Table 2). The aromatics p-coumarate and ferulate were released from cell wall residues by mild alkaline hydrolysis and quantified. All transgenic lines showed reductions of ferulate (−8 to 15%) and the AdoMetase line #1 had 15% less p-coumarate, although the biosynthesis of the latter does not require AdoMet (Table 2).

Table 2 Chemical composition of cell walls isolated from wild-type and AdoMetase lines grown in the greenhouse until full maturity

Lignin degree of polymerization and monomeric composition

2D-HSQC NMR spectroscopy was used to analyze the monomeric composition of lignin (Fig. 3). The data showed a relative increase of S units and a decrease of G units in all four AdoMetase lines, which resulted in higher S/G ratios compared to the WT. A decrease of tricin was also observed in the transgenics. H units were not quantified due to the presence of overlapping signals that arise from protein contamination. The size distribution of lignin oligomers in cell walls was examined using gel permeation chromatography (Fig. 4). The distribution curves of the lignin molecular weights showed both decreased abundance in the lower molecular weight range (100–3000 Da) and increased abundance of higher molecular weight lignin (>3500 Da) in the case of the AdoMetase lines compared to WT. Accordingly, the estimated weight average molecular weight (MW) and number average molecular weight (MN) of lignin are both higher in transgenics, whereas lignin dispersities (Đ) remain similar to that of the WT.

Fig. 3
figure 3

Lignin monome composition in wild-type and AdoMetase lines. Extractive-free stem and leaf biomass from fully mature plants was analyzed by 2D-HSQC NMR spectroscopy. Regions of partial short-range 13C–1H HSQC spectra are shown. For each line, the relative amount of lignin subunits including syringyl (S), guaiacyl (G), and tricin (T), as well as S/G values are provided. pCA: p-coumarate, FA: ferulate, H: p-hydroxyphenyl

Fig. 4
figure 4

Gel-permeation chromatography (GPC) traces of solubilized lignin derived from wild-type and AdoMetase extractive-free biomass. Traces show the molecular weight distribution of lignin oligomers. For each line, lignin weight average molecular weight (MW), number average molecular weight (MN), and dispersity (Đ) are indicated

Saccharification potential

The recalcitrance of engineered sorghum towards enzymatic degradation was evaluated by measuring the amount of sugars released from the biomass after dilute acid and mild alkaline pretreatments followed by a 48-h enzymatic hydrolysis using a mixture of cell wall degrading enzymes (Novozymes Cellic® CTec3). As shown in Fig. 5, higher amount of glucose was released from all AdoMetase lines. After alkaline pretreatment, glucose titers were ~18% higher in hydrolysates from the AdoMetase lines #1 and #2, and ~12% higher for AdoMetase #4. A significant increase in xylose (+15%) was also observed in the case of the AdoMetase #1 (Fig. 5a). After acid pretreatment, 10–21% more glucose was released from the biomass of transgenics compared to the WT (Fig. 5b). Considering that cell walls from both the control and transgenics contain similar amounts of glucose and xylose (Table 2), we conclude that the engineered lines have reduced cell wall recalcitrance and improved saccharification efficiency.

Fig. 5
figure 5

Biomass saccharification of the AdoMetase sorghum lines grown in the greenhouse. The amounts of glucose and xylose released after dilute alkaline (a) and dilute acid (b) pretreatments followed by 48 h of enzymatic digestion with an enzyme cocktail are shown. Values are means ± SE of five biological replicates (n = 5 plants). Asterisks indicate significant differences from the wild type (WT) using the unpaired Student’s t test (*P < 0.05)

Transcriptome analysis

The impact of AdoMetase expression on the transcriptome was investigated by conducting RNA-seq analysis in stems of AdoMetase lines #1 and #3 at three developmental stages. In WT control plants, principal component analysis (PCA) demonstrated a clear distinction between stems analyzed at the three leaf, growing point differentiation (GPD), and boot stages (Additional File 1: Fig. S1a). A total of 28,734 transcripts were identified across all stages in WT, representing 89% of the 32,160 predicted protein-coding genes found in the sorghum genome. The expression level of these genes at the three stages is shown in Additional File 2: Dataset S1. Several transcripts were specific to the three leaf, GPD, and boot stage, while 25,334 genes were found to be expressed at all stages (Additional File 1: Fig. S1b). At each stage, principal component analysis of the transcripts allowed us to distinguish the different genotypes and showed a separation between AdoMetase #1 and #3, indicating alternate transcriptional changes between the two lines (Additional File 1: Fig. S2). The largest number of differentially expressed genes (DEGs) in transgenics compared to WT was observed at the boot stage, accounting for 958 upregulated and 847 downregulated genes in AdoMetase #1, and for 793 upregulated and 426 downregulated genes in AdoMetase #3 (Fig. 6a, Additional File 1: Fig. S3, Additional File 3: Dataset S2). Some overlap was found among DEGs in the two transgenic lines since 369 upregulated genes and 197 downregulated genes were common to both (Fig. 6b, c). However, only a limited number of genes were found consistently upregulated or downregulated at all three developmental stages in a given transgenic line (Additional File 1: Fig. S4). We performed GO and KEGG analyses to determine which biological processes, molecular functions, and metabolic pathways the DEGs identified in transgenics were significantly enriched in. The KEGG analysis revealed in the transgenics at the boot stage several DEGs associated with phenylpropanoid biosynthesis (i.e., 15 peroxidase genes), cytokinin synthesis, fatty acid elongation, cutin and suberin biosynthesis, and chromatin remodeling (Additional File 3: Dataset S2, Additional File 1: Fig. S5a). In particular, several genes associated with DNA metabolism encoding 19 different histone proteins, one histone demethylase, one histone methyltransferase, three DNA cytosine 5-methyltransferases, two DNA-3-methyladenine glycosylases, and one DNA repair protein RAD51 are overexpressed in AdoMetase #3 at the boot stage (yellow highlights in Additional File 3: Dataset S2). The GO analysis also indicated an enrichment of DEGs in biological processes related to cell wall organization, transmembrane transport, and stress and detoxification (Additional File 1: Fig. S5b).

Fig. 6
figure 6

Analysis of DEGs in stems of the AdoMetase lines at the boot stage. a Volcano plots of transcripts identified in stems from WT and AdoMetase lines. The number of downregulated (in blue) and upregulated (in red) transcripts in AdoMetase #1 and AdoMetase #3 compared to WT control is indicated on each plot (log2-fold change +2/−2 and P value < 0.05). Grey dots represent transcripts that are not differentially expressed. Venn diagrams of upregulated (b) and downregulated (c) genes in AdoMetase #1 and #3

Metabolome analysis

Separate aliquots of tissue samples used for transcriptome analysis were used for nontargeted metabolomics. We extracted metabolites with methanol and conducted metabolite profiling with a liquid chromatography–mass spectrometry platform that utilized gradient-based hydrophilic interaction (HILIC) both in the positive and negative ionization modes. The positive ionization allowed the detection of the largest number of unique features (5449) in the samples from WT plants, including 345, 147, and 173 features that are specific to the three leaf, GPD, and boot stage, respectively, whereas 3659 features were detected in all stages (Additional File 1: Fig. S6a, b). Principal component analysis (PCA) demonstrated a clear distinction between the three growth stages and between the different genotypes (Additional File 1: Fig. S6c, d, Fig. S7). Several features had different abundance in AdoMetase #1 and #3 compared to the WT. The highest number of differentially abundant features was found at the boot stage, representing 14.6% and 11.5% of the detected features in AdoMetase #1 and #3, respectively (Fig. 7a, Additional File 1: Fig. S8, Additional File 4: Dataset S3). At the boot stage, 244 more abundant features and 60 less abundant features were common in both transgenic lines compared to the WT (Fig. 7b, c). We also observed 327 features that are absent from both transgenic lines and detected 16 features that are new in both lines (Fig. 7d, e). Moreover, a majority of the differentially abundant features evidenced in the transgenic lines was only observed at a specific growth stage (Additional File 1: Fig. S9). Next, using GNPS, a subset of differentially abundant features observed in transgenics were classified and putatively identified: A majority of the differentially abundant metabolites identified in transgenics belong to the class of phenylpropanoids, alkaloids, and amino acids (Additional File 1: Fig. S10). As examples, the methylated flavonoids selgin, tricin, isorhamnetin and homoeriodictyol were decreased in the AdoMetase lines, whereas the AdoMetase product homoserine and a few derivatives such as threonine, Met, and isoleucine were significantly increased in the transgenics compared to WT (Additional File 1: Fig. S11). The content of the methylated flavone and tricin precursor chrysoeriol was unchanged in transgenics (data not shown), indicating that AdoMet levels are not limiting for its synthesis from luteolin.

Fig. 7
figure 7

Analysis of differentially abundant features in stems of the AdoMetase lines at the boot stage. a Volcano plots of features detected in stems from WT and transgenic lines using HILIC chromatography (positive ionization mode). The number of decreased (in blue) and increased (in red) features in AdoMetase #1 and AdoMetase #3 compared to WT control is indicated on each plot (log2-fold change +2/−2 and P value < 0.05). Gray dots represent features that are not differentially abundant. Venn diagrams of features that are more abundant (b), less abundant (c), new (d), and absent (e) in AdoMetase #1 and #3

Field testing

The sorghum lines AdoMetase #1 and #3 were grown in the field until the soft dough stage in two different locations to evaluate their performance under natural environment. A composition analysis was performed on extractive-free cell wall residues obtained from stems and leaves. In AdoMetase #1, the content of acid-insoluble lignin was significantly reduced by 18% and 13% at the two different field sites, respectively, whereas line Adometase #3 showed no significant change (P < 0.05) (Table 3). For the Davis site, glucuronic acid content was increased by 43% and 16% in AdoMetase #1 and #3, respectively, whereas 4-O-methylglucuronic acid was reduced in AdoMetase #1 (8–27%) at both sites. As a result, the degree of glucuronic acid methylation was significantly reduced in both lines in the first site (P < 0.05). The amount of cell-wall-bound p-coumarate and ferulate in the transgenic lines was unchanged compared to the WT. Measurement of stover biomass showed yield reductions for both AdoMetase lines in the two locations (Additional File 1: Fig. S12). The yield penalty was more pronounced at the Davis site where the two lines were affected differently (69% reduction for AdoMetase #1 and 49% for AdoMetase #3) compared to the KARE site (~23% less biomass for both lines). In addition, the WT produced less biomass at the KARE site, possibly due to the sandy soil type and warmer temperatures at this location.

Table 3 Chemical composition of cell walls isolated from wild-type and AdoMetase lines grown in the field until the soft dough stage

Saccharification assays were conducted on the biomass from stems and leaves. As previously observed in the greenhouse experiment, the two engineered sorghum lines showed improvements of saccharification efficiency, releasing up to 18% and 24% more glucose after dilute alkaline and dilute acid pretreatments, respectively (Fig. 8).

Fig. 8
figure 8

Biomass saccharification of the sorghum lines AdoMetase #1 and AdoMetase #3 grown in two different field sites in California. The data shows the amounts of glucose and xylose released after dilute alkaline and dilute acid pretreatments followed by 48 h of enzymatic digestion with an enzyme cocktail. Values are means ± SE of four biological replicates (n = 4 plots). For each site and pretreatment method, asterisks indicate significant differences from the wild type (WT) using the unpaired Student’s t test (*P < 0.1, **P < 0.05)

Discussion

In sorghum, genetic mutations that introduce premature stop codons in a gene encoding a COMT (Sobic.007G047300) lead to brown midrib (bmr) phenotypes due to changes in lignin composition [26]. For example, biomass from the well-studied bmr12 COMT mutant contains ~6% less acid-insoluble lignin and shows a decrease in the lignin S/G ratio and tricin content, resulting in slight improvements of sugar yields (<10%) after dilute acid or mild alkali biomass pretreatments followed by enzymatic saccharification [14, 27, 28]. These observations also suggest that other COMTs besides Sobic.007G047300 are involved in lignin biosynthesis in sorghum since S and T units remain present in the bmr12 mutant. In the current study, the heterologous expression of AdoMetase in sorghum reduces AdoMet content and potentially affects the activity of any OMTs that participate in lignin synthesis, resulting in up to 18–19% reductions of acid-insoluble lignin in the greenhouse and the field. The AdoMetase sorghum lines show delayed flowering, as previously observed in bmr12 mutants, but they also show reductions of plant height and biomass yields, which is in contrast with bmr12 plants [29,30,31,32]. Such a short phenotype was previously observed in the sorghum bmr19 mutant that is presumably affected in AdoMet synthesis due to a mutation in a FPGS gene [33]. Thus, by using the pSbCOMT promoter to drive AdoMetase expression, our results suggest that not only COMT but also CCoAOMTs are impacted in their activity in the AdoMetase lines. First, unlike the bmr12 mutants that have reduced lignin S/G ratios, the AdoMetase lines at maturity show an increased S/G, as previously observed in the bmr19 mutant and maize plants silenced in a CCoAOMT gene or carrying mutations in FPGS (bm4) or MTHFR (bm2) [19, 22, 34, 35]. Arguably, it is not excluded that AdoMetase expression declines at later stages during the reproductive growth phase, which would marginally affect lignin biosynthesis in maturing stems that typically produce S-rich lignin, and result in lignin with higher S/G ratio as observed in fully mature transgenic plants [36]. Second, an increase of ferulate and a decrease of p-coumarate bound to the cell wall were measured in bmr12 mutants [28], which is in contrast with our AdoMetase lines that show no change in cell-wall-bound p-coumarate and a reduction of ferulate at maturity (Table 2). The lignin in the AdoMetase lines shows increased MW and MN, which is different from the data usually observed in engineered plants that have higher lignin S/G [37]. This difference is possibly due to the reduction of T units since MW and MN were shown to be negatively correlated with tricin content in isolated lignin fractions, consistent with the exclusive occurrence of tricin as pendant end-groups in lignin and its putative role as nucleation sites for lignin polymerization [5, 38].

Based on our RNA-seq data, we did not find any significant changes in the expression of known genes (i.e. Bmr2, Bmr6, Bmr12, Bmr19, Bmr30, and SbMYB60) or putative genes involved in lignin biosynthesis in the stems of the AdoMetase lines at the three different growth stages (Additional File 3: Dataset S2, log2-fold change +2/−2 and P < 0.05) [18, 26, 39,40,41,42,43]. Moreover, none of the seven CCoAOMT genes found in the sorghum genome showed differential expression in transgenics compared to the WT [44]. We found that 29 out of the 40 OMT genes present in the sorghum genome are expressed in stems (Additional File 2: Dataset S1) [45]. Among those, only two genes are significantly downregulated in both AdoMetase lines: Sobic.007G170500 (three-leaf stage) and Sobic.002G079500 (boot stage) (Additional File 3: Dataset S2). Bmr12 (Sobic.007G047300) is the OMT gene with the highest expression level, showing DESeq2 normalized values of 23,525, 26,202.5, 33,304 at the three leaf, GPD, and boot stage, respectively, followed by Sobic.009G197000 (i.e. DESeq2 values of 8,296, 11,474, and 18,576 at the three above-mentioned stages). The 27 other OMT genes show much lower expression levels in the stem (i.e. DESeq2 values <645). Overall, this data indicates that the reduction of lignin measured in the AdoMetase lines is probably the result of lower cytosolic AdoMet pools rather than downregulation of the genes involved in lignin biosynthesis. Moreover, the reduction of the degree of glucuronic acid methylation on xylan indicates a diminution of the AdoMet pool inside the Golgi where the corresponding glucuronic acid methyltransferases are located [46].

The RNA-seq data indicates that biological processes other than lignification and xylan methylation may be impacted in the AdoMetase lines, which is perhaps not surprising considering the ubiquitous role of the AdoMet cofactor. Notably, the expression of AdoMetase modifies the expression of several genes encoding histone and DNA base methyltransferases, as wells as histone proteins. These changes in the expression of enzymes involved in DNA metabolism could result in epigenetic modifications that affect gene expression more broadly, and explain the delay in growth and development observed in the AdoMetase lines under field conditions [47]. Changes in the expression of genes involved in the metabolism of the cytokinin zeatin, which derives from the nucleobase adenine, also suggest some possible perturbations of hormone signaling in the AdoMetase lines. For example, one of the genes encoding 1-aminocyclopropane-1-carboxylate synthase (Sobic.001G105900), an ethylene biosynthetic enzyme that uses AdoMet as cofactor, is upregulated in both AdoMetase lines at the three-leaf stage, presumably as part of a compensation mechanism (Additional File 3: Data S2). In addition, the overexpression of multiple peroxidase genes in transgenics is indicative of a stress response due to higher levels of reactive oxygen species.

At the metabolite level, several compounds belonging to the class of phenylpropanoids and alkaloids are differentially abundant in the AdoMetase lines. O-methyltransferases are indeed involved in the methylation of a wide range of secondary metabolites such as alkaloids and aromatics, so that reductions in AdoMet pools are expected to impact their relative amounts [48]. For example, we showed that two flavones synthesized by SbCOMT/Bmr12 (selgin and tricin) are reduced in metabolite extracts from the AdoMetase lines compared to wild-type plants. Other unknown O-methyltransferases involved in the synthesis of methylated compounds such as isorhamnetin and homoeriodictyol are probably affected by a reduction of AdoMet since these two flavonoids are also reduced in transgenics (Additional File 1: Fig. S11a). The compound class ‘amino acids’ is also represented in the differentially abundant metabolites identified in the AdoMetase lines. In particular, the AdoMetase product homoserine appears to be partially recycled into threonine, isoleucine, and Met, which are all increased in transgenics (Additional File 1: Fig. S11b). The increase in Met could represent a cellular response to compensate for the reduction of AdoMet pools since Met is the AdoMet precursor. Met accumulation may also occur from the recycling of 5ʹ-methylthioadenosine (i.e. the other AdoMetase product) via the Yang cycle because the early intermediate 5-methylthioribose of this cycle is increased in the transgenic lines (Additional File 1: Fig. S11b). A significant increase of S-adenosylhomocysteine (AdoCys) was also observed (Additional File 1: Fig. S11b), which is somehow unexpected since AdoCys is a by-product of AdoMet-dependent methyltransferases that usually undergoes hydrolysis into adenosine and homocysteine by AdoCys hydrolase. Whether higher AdoCys content is the consequence of AdoCys hydrolase reverse reaction or another unknown enzymatic activity in transgenics remains to be determined [49, 50]. Interestingly, augmenting AdoCys levels was achieved in switchgrass by silencing AdoCys hydrolase, which successfully led to lignin reductions, supposedly via inhibition of COMTs/CCoAOMTs by AdoCys [51].

Conclusions

This study demonstrates the effectiveness of expressing AdoMetase in sorghum to reduce lignin and improve enzymatic conversion of biomass to simple sugars. Our work also highlights in transgenic lines some pleiotropic effects on metabolic pathways that require AdoMet for methylation reactions. Thus, fine-tuning AdoMetase expression by using promoters active at later developmental stages should be evaluated in order to mitigate adverse effects on plant growth. For example, a promoter that is induced in response to shorter day length could be appropriate in the case of late planting during the growing season. Moreover, introducing AdoMetase in energy sorghum varieties could represent a viable approach to reduce biomass recalcitrance and maintain high biomass yields [52].

Materials and methods

Plant growth conditions

Plants were grown at the UC Berkeley greenhouse Oxford facility with minimum temperature set at 22 °C. Seeds were germinated directly on soil (Sunshine mix #4, Sun Gro) in one-gallon pots and plants were grown until seeds reached the black layer stage (i.e. full physiological maturity). One tablespoon of Osmocote Plus 15-9-12 was added to the soil biweekly until the flowering stage. Watering was stopped at the end of the growing period and pots containing plants were allowed to dry for another 3 weeks. Plants without their panicles were harvested, further dried in an oven at 50 °C for 5 days, and subsequently ground using a knife mill equipped with a 1-mm mesh (Model 4 Wiley Mill, Thomas Scientific). For cell wall composition analyses, samples were additionally ground into a fine powder using a ball mill (Mixer Mill MM 400, Retsch Inc.) with stainless-steel balls.

Design of the pSbCOMT:AdoMetase construct

The promoter sequence of the SbCOMT gene (1990 bp) was amplified by PCR using genomic DNA from Sorghum bicolor (variety BTx623) and the primers listed in Additional File 1: Table S1. The fragment was cloned into the pBca9145 vector [53] via In-Fusion cloning (Takara Bio USA) to generate a level-0 construct. A DNA sequence encoding AdoMetase from the enterobacteria phage T3 (GenBank: CAA28477.1) was codon-optimized for expression in sorghum and synthesized as a gene fragment by GenScript. The sequence contained flanking BsaI restriction sites plus extra homologous sequences for cloning into pBca9145. The pSbCOMT:AdoMetase construct was obtained using the jStack cloning method [53]. The corresponding level-0 and level-1 intermediate plasmids are listed in Additional File 1: Table S2. Plasmid sequences are available at the Inventory of Composable Elements source registry (http://public-registry.jbei.org).

Sorghum transformation and genotyping

The Agrobacterium tumefaciens EHA105 strain was used to transform sorghum (S. bicolor, variety Wheatland) using immature embryos as previously described [54]. TaqMan Real-time PCR assays (Thermo Fisher Scientific) were performed on gDNA isolated from primary transformants using primers specific to the nptII and AdoMetase gene sequences to identify transformants with a single-copy event. Plants homozygous for the transgene and WT segregants were identified using the same approach for four selected lines in the T1 generation. T2 seeds from one T1 WT segregant plant from each line were pooled and used as WT controls. Plants in the T3 generation were used for all greenhouse experiments.

Analysis of AdoMetase expression

AdoMetase expression was assessed using reverse transcription quantitative PCR (RT-qPCR). Total RNAs were extracted from stems of 3-week-old plants using the RNeasy Plant Mini Kit (Qiagen) and first-strand cDNA synthesis was conducted using the SuperScript IV First-Strand Synthesis kit (Thermo Fisher Scientific) as previously described [55]. qPCR on first-strand cDNA was performed using 35 cycles consisting of 5 s at 95 °C for denaturation and 15 s at 60 °C for annealing and amplification. The relative quantification of AdoMetase transcripts was calculated using the 2−ΔCT method and normalized to the reference gene PP2A [56]. The results are the average from four biological replicates. RT-qPCR primers are listed in Additional File 1: Table S1.

RNA-seq analysis

Four biological replicates from two AdoMetase lines (#1 and #3) and WT control were used for large-scale transcriptomic analysis. Entire stems from seedlings (three-leaf stage) and 3-cm stem sections collected 5 cm from the soil surface (growing point differentiation and boot stages) were flash-frozen in liquid nitrogen and stored at −80 °C until processing. Frozen samples were pulverized using a freezer mill (model 6875D, SPEX SamplePrep LLC). Each sample was ground twice for 1 min at a grinding rate of 15 counts per sec, including a 15 s break between the two cycles. Total RNA was extracted using the RNeasy Plant Mini Kit (Qiagen) and treated with RNAse-free DNAse. RNA sequencing, RNA-seq read filtering, and raw gene counts was done as described previously [57]. Filtered reads from each library were aligned to the Sorghum bicolor v5.1 reference genome (https://phytozome-next.jgi.doe.gov/info/Sbicolor_v5_1) using HISAT2 version 2.2.1 [58].

Metabolite analysis

To quantify AdoMet, 5ʹ-methylthioadenosine, homoserine, threonine, isoleucine, AdoCys, methionine, and 5-methylthioribose, metabolites were extracted from 3-week-old plants using 5% (w/v) trichloroacetic acid and analyzed using high-performance liquid chromatography (HPLC), electrospray ionization (ESI), and time-of-flight (TOF) mass spectrometry (MS) as previously described [24]. Flavonoids were extracted from stems and leaves of fully mature plants and quantified after acid-hydrolysis of the metabolite extracts using the LC–MS method described in [14]. Measurement of the metabolites was performed via the analysis of standard compounds using the same LC–MS methods. 5-Methylthioribose was identified based on accurate mass measurements of LC–MS peaks ([M + H]+ = 181.0529 at <3.5 ppm) and CID–MS/MS fragmentation patterns showing major [M + H]+ ions at 57.0340, 71.0491, 85.0284, 117.0546, 121.0318, 149.0267, and 163.0423 m/z using an collision energy of 10 eV.

Aliquots of the homogenized biomass samples used for RNA-seq analysis were used for untargeted metabolomics. To obtain metabolite extracts, frozen tissue powder was weighed and mixed with methanol at a ratio of 100 µL of solvent per mg of powder. Samples were vortexed for 1 min, then incubated at room temperature for 20 min with continuous mixing, centrifuged at 20,000×g for 5 min and the supernatant filtered through 0.45-µm PTFE filters. Three extraction control samples were obtained by following the same procedure, but without any tissue powder. Filtered metabolite extracts and controls were vacuum-dried and resuspended in 100% methanol containing an internal standard mix of isotopically labeled compounds for UHPLC–MS/MS normal phase chromatography using a HILIC column (InfinityLab Poroshell 120 HILIC-Z, 2.1 × 150 mm, 2.7 µm, Agilent) and following the method described in [57]. The raw data from the metabolite LC–MS runs is provided in Additional File 5: Dataset S4. For data analysis, features with 0.8 min < RT < 18.5 min and max peak height fold-change between sample and extraction control >3 were selected. LC–MS data was analyzed using custom Python code [59], with each detected peak assigned a level of confidence indicated by a score from 0 to 3 in the compound identification. We performed a Feature-Based Molecular Networking workflow using MZmine 2 [60] and Global Natural Products Social Molecular Networking (GNPS; http://gnps.ucsd.edu) [61]. The MZmine workflow was used to generate a list of features obtained from extracted ion chromatograms containing chromatographic peaks within a narrow m/z range and filtered to remove isotopes. For each feature, the most intense fragmentation spectrum was uploaded to GNPS for putative identification by comparison with mass spectra deposited in the database. Compound classes are attributed to identified compounds using ClassyFire [62] or NPClassifier [63].

Cell wall composition analysis

Ball-milled biomass (1 g) was sequentially extracted using a Dionex ASE 350 accelerated solvent extractor set to 7 min static extraction cycles (Thermo Fisher Scientific). Solvents (5 mL) were water (two cycles), 80% (v/v) ethanol–water (ten cycles), 50% (v/v) methanol–chloroform (one cycle), and acetone (one cycle). Acid-insoluble lignin, acid-soluble lignin (absorbance at 320 nm), cell wall monosaccharides (glucose, xylose, and arabinose), and cell-wall-bound aromatics (p-coumarate and ferulate) were measured as previously described [64, 65]. Glucuronic acid and 4-O-methylglucuronic acid were quantified using HPLC–ESI–qTOF–MS analysis and a BioRad Fermentation Monitoring HPX-87H 150 × 7.8 mm column. The mobile phase was 0.1% formic acid in water at a flow rate of 0.50 mL/min and a total run time of 13 min. Other parameters were as follows: Column temperature: 50 °C, maximum pressure: 100 bar, ionization mode: negative, gas temperature: 340 °C, drying gas flow: 11 L/min, nebulizer: 30 psig, fragmentor: 140 V, skimmer: 50 V, octopole radio frequency voltage: 195 V, VCap; 3500 V, acquisition range (m/z): 50–700, acquisition rate (spectra/s): 0.86.

2D-HSQC NMR

For each sorghum line, the sample consisted of a pool of equal amount of extractive-free cell wall residue from five plants. Preparation of the samples for 2D-HSQC NMR experiments was described previously [66, 67]. In a 5-mm NMR tube, ~ 100 mg of extracted biomass was added along with 1 ml of pre-mixed DMSO–d6/pyridine–d5 solvent (4:1, v/v) to form a gel. The NMR tubes were sealed and sonicated for 4 h with a 30 min interval every hour until the gel became apparently homogeneous. The 2D-HSQC NMR spectra were collected on a Bruker Avance I 800 MHz spectrometer equipped with a Bruker Triple Resonance Probe (TXI) at 310 K. A standard Bruker pulse sequence (hsqcetgpsisp2.2) was used with parameters that are typical for lignocellulose samples. Data were acquired using Bruker’s TopSpin software (version 4.1.0). The HSQC spectra were collected from 11 to −1 ppm in the F2 (1H) dimension with 1,024 data points and 53 ms acquisition time with an interscan pulse delay of 1 s, and from 165 to −10 ppm in the F1 (13C) dimension with 256 data points and 3.5 ms acquisition time. For each evolution period (t1) increment, 256 scans were recorded. The central DMSO solvent peak was used as a reference for the chemical shift calibration for all samples (δC 39.5 ppm, δH 2.5 ppm). All HSQC spectra were processed using typical cosine-squared apodization in both F2 and F1 dimensions, and the contours were integrated using MestreNOVA (version 14, Mestrelab Research). Peaks were assigned according to published data [66,67,68,69].

Gel permeation chromatography (GPC)

Pools of extracted biomass from five plants were used for each line. Samples were incubated in 2.5 mL of acetic acid and acetyl bromide (92:8) and stirred at 50 °C for 2 h to dissolve the lignin [70]. Acetylated lignin was dissolved in tetrahydrofuran and filtered through 0.2-µm PTFE filters. GPC analysis of biomass-derived lignin fragments was conducted using an Ecosec HLC-8320GPC (Tosoh Bioscience) equipped with Agilent Technologies PLgel 5 μm Mixed-D column and a diode array detector. Tetrahydrofuran, spiked with 250 ppm of butylated hydroxytoluene, was used as the mobile phase with a flow rate of 1 ml min−1 and a column temperature of 40 °C. Polystyrene oligomers ranging from 162 to 29,150 g mol−1 (Agilent Technologies #PL2013) were dissolved in tetrahydrofuran and used as standards for calibration using the conditions described above.

Saccharification

The biomass from leaves and stems (panicles removed) was used for saccharification assays. Dilute acid (1.2% w/v sulfuric acid) and alkali (0.25% w/v sodium hydroxide) were used for biomass pretreatments. The enzyme mixture Ctec3 (Novozymes) was used for polysaccharide hydrolysis. The procedure and measurement of released monosaccharides were performed as previously described [71, 72].

Field trial

Seeds (T4 generation) from AdoMetase #1, AdoMetase #3, and pooled wild-type segregants were planted in two field trials conducted at the University of California Davis Plant Sciences Research Farm (‘Davis site’, Davis, CA) and the Kearney Agricultural Research and Extension Centre (‘KARE site’, Parlier, CA) in 2023 under a USDA–APHIS Biotechnology Regulatory Services permit for regulated sorghum (BRS# 124-Q89IS3R, AUTH# 241864) as previously described [73, 74]. The Davis site has a Yolo clay loam soil (fine-silty, mixed thermic Fluventic Haploxerept) and the KARE site has a Hanford sandy loam soil. Seeds were planted on June 22, 2023. Percent germination was taken for each line, and adjusted for a target seeding rate of 197,600 seeds ha−1. The experiment was a randomized complete block design with four replicates. Row spacing was 0.76 m. Plots measured 3 × 6 m and comprised two rows. A commercial sorghum variety (NK8828 from S&W Seed Company) was planted as border rows to limit edge effects. Before planting, the amount of N–P–K fertilizer applied at the Davis site was 224, 91, and 22 kg ha−1, respectively. At the KARE site, the amount of N–P–K was 115, 58, and 233 kg ha−1, respectively, with an additional 77 kg ha−1 of in-season nitrogen from a urea source injected into the drip system. Sivanto was applied at a rate of 481 g ha−1 on August 7 in the KARE site for sugarcane aphid control. Irrigation water was applied utilizing surface furrow methods at the Davis site and drip irrigation at the KARE site to satisfy the fully watered evapotranspiration requirement for sorghum during the growing season. Evapotranspiration was estimated from nearby weather stations. All panicles on each plant were covered with pollinating bags prior to anthesis to prevent pollen flow according to APHIS permit regulations. Plot harvest took place on October 27, 2023. Panicles (heads with immature grain) were removed by hand. Remaining stover was harvested using a Wintersteiger Cibus forage chopper (Wintersteiger). Subsamples of approximately 500 g of fresh weight stover material were taken from each harvested plot, and dried at 55 °C for 1 week in a forced air oven to determine dry matter and calculate yield. For composition analyses and saccharification assays, stover samples were first milled using a Model 4 Wiley Mill equipped with a 1-mm mesh (Thomas Scientific) and further ground into powder as previously described.

Data analysis

Principal component analysis plots were obtained with MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/) and Venn diagrams were generated with the webtool available at https://bioinformatics.psb.ugent.be/webtools/Venn/. For the Gene Ontology (GO) enrichment analysis, the web-based g:Profiler tool was used to identify overrepresented biological processes in the AdoMetase lines [75]. The build-in algorithms of g:Profiler was used to correct for multiple testing and to identify the leading GO terms and remove the redundant terms in each function component. Leading GO terms that were significantly enriched (adjusted P value ≤0.05) in at least one AdoMetase line were visualized in a dot plot using the ggplot2 package v3.4.4 in R v4.3.1 [76, 77]. For the KEGG analysis, enrichKEGG function with the clusterProfiler package v4.8.3 [78] in R was applied to identified overrepresented KEGG categories in the AdoMetase lines based on the UniProt IDs of their differentially expressed genes. The P value was adjusted using FDR method with a cutoff of 0.1. The significantly enriched KEGG categories were visualized in a dot plot as described above.

Availability of data and materials

All the data supporting the conclusions of this article is included within the article and its additional files. The RNA-seq raw data sets in this study are available at the NCBI Sequence Read Archive under BioProject IDs PRJNA1098655–PRJNA1098666 and PRJNA1098671–PRJNA1098695. The raw metabolite data can be accessed at the MassIVE public depository under the data set identifier MSV000094762.

Abbreviations

AdoCys:

S-Adenosylhomocysteine

AdoMet:

S-Adenosylmethionine

AdoMetase:

S-adenosylmethionine hydrolase

bmr:

Brown midrib

CCoAOMT:

Caffeoyl-CoA O-methyltransferase

COMT:

Caffeic acid O-methyltransferase

FPGS:

Folylpolyglutamate synthase

GPC:

Gel-permeation chromatography

GPD:

Growing point differentiation

MTHFR:

5-Methyltetrahydrofolate reductase

PCA:

Principal component analysis

References

  1. Ashokkumar V, Venkatkarthick R, Jayashree S, Chuetor S, Dharmaraj S, Kumar G, et al. Recent advances in lignocellulosic biomass for biofuels and value-added bioproducts - A critical review. Bioresour Technol. 2022;344(Pt B):126195.

    Article  CAS  PubMed  Google Scholar 

  2. Langholtz M. 2023 Billion-Ton Report. Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States), Bioenergy Knowledge Discovery Framework (BEKDF)

  3. Meents MJ, Watanabe Y, Samuels AL. The cell biology of secondary cell wall biosynthesis. Ann Bot. 2018;121:1107–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Boerjan W, Ralph J, Baucher M. Lignin biosynthesis. Annu Rev Plant Biol. 2003;54:519–46.

    Article  CAS  PubMed  Google Scholar 

  5. Lan W, Rencoret J, Lu F, Karlen SD, Smith BG, Harris PJ, et al. Tricin-lignins: occurrence and quantitation of tricin in relation to phylogeny. Plant J. 2016;88:1046–57.

    Article  CAS  PubMed  Google Scholar 

  6. Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010;153:895–905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. del Rio JC, Rencoret J, Gutierrez A, Elder T, Kim H, Ralph J. Lignin monomers from beyond the canonical monolignol biosynthetic pathway: Another brick in the wall. ACS Sustain Chem Eng. 2020;8:4997–5012.

    Article  Google Scholar 

  8. Baral NR, Sundstrom ER, Das L, Gladden JM, Eudes A, Mortimer J, et al. Approaches for more efficient biological conversion of lignocellulosic feedstocks to biofuels and bioproducts. ACS Sustain Chem Eng. 2019;7:9062–79.

    Article  CAS  Google Scholar 

  9. Liu C-J, Eudes A. Lignin synthesis and bioengineering approaches toward lignin modification. In: Sibout R, editor. Advances in botanical research, vol. 104. Amsterdam: Elsevier; 2022. p. 41–96.

    Google Scholar 

  10. Yoshida K, Sakamoto S, Mitsuda N. In planta cell wall engineering: from mutants to artificial cell walls. Plant Cell Physiol. 2021;62:1813–27.

    Article  CAS  PubMed  Google Scholar 

  11. Hanson AD, Roje S. One-carbon metabolism in higher plants. Annu Rev Plant Physiol Plant Mol Biol. 2001;52:119–37.

    Article  CAS  PubMed  Google Scholar 

  12. Roje S. S-adenosyl-L-methionine: beyond the universal methyl group donor. Phytochemistry. 2006;67:1686–98.

    Article  CAS  PubMed  Google Scholar 

  13. Lam LPY, Lui ACW, Bartley LE, Mikami B, Umezawa T, Lo C. Multifunctional 5-hydroxyconiferaldehyde O-methyltransferases (CAldOMTs) in plant metabolism. J Exp Bot. 2024;75:1671–95.

    Article  CAS  PubMed  Google Scholar 

  14. Eudes A, Dutta T, Deng K, Jacquet N, Sinha A, Benites VT, et al. SbCOMT (Bmr12) is involved in the biosynthesis of tricin-lignin in sorghum. PLoS ONE. 2017;12: e0178160.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Fornalé S, Rencoret J, García-Calvo L, Encina A, Rigau J, Gutiérrez A, et al. Changes in cell wall polymers and degradability in maize mutants lacking 3’- and 5’-O-methyltransferases involved in lignin biosynthesis. Plant Cell Physiol. 2017;58:240–55.

    PubMed  Google Scholar 

  16. Lam PY, Tobimatsu Y, Matsumoto N, Suzuki S, Lan W, Takeda Y, et al. OsCAldOMT1 is a bifunctional O-methyltransferase involved in the biosynthesis of tricin-lignins in rice cell walls. Sci Rep. 2019;9:11597.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Shen B, Li C, Tarczynski MC. High free-methionine and decreased lignin content result from a mutation in the Arabidopsis S-adenosyl-L-methionine synthetase 3 gene. Plant J. 2002;29:371–80.

    Article  CAS  PubMed  Google Scholar 

  18. Adeyanju AO, Sattler SE, Rich PJ, Rivera-Burgos LA, Xu X, Ejeta G. Sorghum brown midrib19 (bmr19) gene links lignin biosynthesis to folate metabolism. Genes. 2021;12(5):660.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li L, Hill-Skinner S, Liu S, Beuchle D, Tang HM, Yeh C-T, et al. The maize brown midrib4 (bm4) gene encodes a functional folylpolyglutamate synthase. Plant J. 2015;81:493–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Mazarei M, Baxter HL, Srivastava A, Li G, Xie H, Dumitrache A, et al. Silencing folylpolyglutamate synthetase1 (FPGS1) in switchgrass (Panicum virgatum L.) improves lignocellulosic biofuel production. Front Plant Sci. 2020;11:843.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Srivastava AC, Chen F, Ray T, Pattathil S, Peña MJ, Avci U, et al. Loss of function of folylpolyglutamate synthetase 1 reduces lignin content and improves cell wall digestibility in Arabidopsis. Biotechnol Biofuels. 2015;8:224.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Tang HM, Liu S, Hill-Skinner S, Wu W, Reed D, Yeh C-T, et al. The maize brown midrib2 (bm2) gene encodes a methylenetetrahydrofolate reductase that contributes to lignin accumulation. Plant J. 2014;77:380–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Barrière Y, Ralph J, Méchin V, Guillaumie S, Grabber JH, Argillier O, et al. Genetic and molecular basis of grass cell wall biosynthesis and degradability. II. Lessons from brown-midrib mutants. C R Biol. 2004;327:847–60.

    Article  PubMed  Google Scholar 

  24. Eudes A, Zhao N, Sathitsuksanoh N, Baidoo EEK, Lao J, Wang G, et al. Expression of S-adenosylmethionine hydrolase in tissues synthesizing secondary cell walls alters specific methylated cell wall fractions and improves biomass digestibility. Front Bioeng Biotechnol. 2016;4:58.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Amthor JS. Efficiency of lignin biosynthesis: a quantitative analysis. Ann Bot. 2003;91:673–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bout S, Vermerris W. A candidate-gene approach to clone the sorghum Brown midrib gene encoding caffeic acid O-methyltransferase. Mol Genet Genomics. 2003;269:205–14.

    Article  CAS  PubMed  Google Scholar 

  27. Godin B, Nagle N, Sattler S, Agneessens R, Delcarte J, Wolfrum E. Improved sugar yields from biomass sorghum feedstocks: comparing low-lignin mutants and pretreatment chemistries. Biotechnol Biofuels. 2016;9:251.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Singh SK, Saulnier BK, Hodge DB. Lignin properties and cell wall response to deconstruction by alkaline pretreatment and enzymatic hydrolysis in brown midrib sorghums. Ind Crops Prod. 2022;178: 114566.

    Article  CAS  Google Scholar 

  29. Oliver AL, Pedersen JF, Grant RJ, Klopfenstein TJ, Jose HD. Comparative effects of the sorghum -6 and -12 genes. Crop Sci. 2005;45:2240.

    Article  CAS  Google Scholar 

  30. Pedersen JF, Toy JJ, Funnell DL, Sattler SE, Oliver AL, Grant RA. Registration of BN611, AN612, BN612, and RN613 sorghum genetic stocks with stacked bmr -6 and bmr -12 genes. J Plant Regist. 2008;2:258–62.

    Article  Google Scholar 

  31. Guragain YN, Srinivasa Rao P, Vara Prasad PV, Vadlani PV. Evaluation of brown midrib sorghum mutants as a potential biomass feedstock for 2,3-butanediol biosynthesis. Appl Biochem Biotechnol. 2017;183:1093–110.

    Article  CAS  PubMed  Google Scholar 

  32. Sattler SE, Funnell-Harris DL, Pedersen JF. Efficacy of singular and stacked brown midrib 6 and 12 in the modification of lignocellulose and grain chemistry. J Agric Food Chem. 2010;58:3611–6.

    Article  CAS  PubMed  Google Scholar 

  33. Funnell-Harris DL, Sattler SE, Toy JJ, O’Neill PM, Bernhardson LF. Response of sorghum lines carrying recently identified brown midrib (bmr) mutations to stalk rot pathogens and water deficit. Plant Pathol. 2023;72:752–64.

    Article  CAS  Google Scholar 

  34. Saballos A, Vermerris W, Rivera L, Ejeta G. Allelic association, chemical characterization and saccharification properties of brown midrib mutants of sorghum (Sorghum bicolor (L.) Moench). Bioenerg Res. 2008;1:193–204.

    Article  Google Scholar 

  35. Li X, Chen W, Zhao Y, Xiang Y, Jiang H, Zhu S, et al. Downregulation of caffeoyl-CoA O-methyltransferase (CCoAOMT) by RNA interference leads to reduced lignin production in maize straw. Genet Mol Biol. 2013;36:540–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Waliszewska H, Zborowska M, Stachowiak-Wencek A, Waliszewska B, Czekała W. Lignin transformation of one-year-old plants during anaerobic digestion (AD). Polymers. 2019;11(5):835.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Stewart JJ, Akiyama T, Chapple C, Ralph J, Mansfield SD. The effects on lignin structure of overexpression of ferulate 5-hydroxylase in hybrid poplar. Plant Physiol. 2009;150:621–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Zikeli F, Ters T, Fackler K, Srebotnik E, Li J. Fractionation of wheat straw dioxane lignin reveals molar mass dependent structural differences. Ind Crops Prod. 2016;91:186–93.

    Article  CAS  Google Scholar 

  39. Saballos A, Ejeta G, Sanchez E, Kang C, Vermerris W. A genome wide analysis of the cinnamyl alcohol dehydrogenase family in sorghum [Sorghum bicolor (L.) Moench] identifies SbCAD2 as the brown midrib6 gene. Genetics. 2009;181:783–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Saballos A, Sattler SE, Sanchez E, Foster TP, Xin Z, Kang C, et al. Brown midrib2 (Bmr2) encodes the major 4-coumarate:coenzyme A ligase involved in lignin biosynthesis in sorghum (Sorghum bicolor (L.) Moench). Plant J. 2012;70:818–30.

    Article  CAS  PubMed  Google Scholar 

  41. McKinley B, Rooney W, Wilkerson C, Mullet J. Dynamics of biomass partitioning, stem gene expression, cell wall biosynthesis, and sucrose accumulation during development of Sorghum bicolor. Plant J. 2016;88:662–80.

    Article  CAS  PubMed  Google Scholar 

  42. Scully ED, Gries T, Sarath G, Palmer NA, Baird L, Serapiglia MJ, et al. Overexpression of SbMyb60 impacts phenylpropanoid biosynthesis and alters secondary cell wall composition in Sorghum bicolor. Plant J. 2016;85:378–95.

    Article  CAS  PubMed  Google Scholar 

  43. Tetreault HM, Gries T, Liu S, Toy J, Xin Z, Vermerris W, et al. The sorghum (Sorghum bicolor) brown midrib 30 gene encodes a chalcone isomerase required for cell wall lignification. Front Plant Sci. 2021;12: 732307.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Rakoczy M, Femiak I, Alejska M, Figlerowicz M, Podkowinski J. Sorghum CCoAOMT and CCoAOMT-like gene evolution, structure, expression and the role of conserved amino acids in protein activity. Mol Genet Genomics. 2018;293:1077–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Liu Z, Fan M, Li C, Xu J-H. Dynamic gene amplification and function diversification of grass-specific O-methyltransferase gene family. Genomics. 2019;111:687–95.

    Article  CAS  PubMed  Google Scholar 

  46. Urbanowicz BR, Peña MJ, Ratnaparkhe S, Avci U, Backe J, Steet HF, et al. 4-O-methylation of glucuronic acid in Arabidopsis glucuronoxylan is catalyzed by a domain of unknown function family 579 protein. Proc Natl Acad Sci USA. 2012;109:14253–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Kumar S, Mohapatra T. Dynamics of DNA methylation and its functions in plant growth and development. Front Plant Sci. 2021;12: 596236.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Lashley A, Miller R, Provenzano S, Jarecki S-A, Erba P, Salim V. Functional diversification and structural origins of plant natural product methyltransferases. Molecules. 2022;28(1):43.

    Article  PubMed  PubMed Central  Google Scholar 

  49. De La Haba G, Cantoni GL. The enzymatic synthesis of S-adenosyl-L-homocysteine from adenosine and homocysteine. J Biol Chem. 1959;234:603–8.

    Article  PubMed  Google Scholar 

  50. Brzezinski K, Dauter Z, Jaskolski M. High-resolution structures of complexes of plant S-adenosyl-L-homocysteine hydrolase (Lupinus luteus). Acta Crystallogr D Biol Crystallogr. 2012;68(Pt 3):218–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Bai Z, Qi T, Liu Y, Wu Z, Ma L, Liu W, et al. Alteration of S-adenosylhomocysteine levels affects lignin biosynthesis in switchgrass. Plant Biotechnol J. 2018;16:2016–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Mullet J, Morishige D, McCormick R, Truong S, Hilley J, McKinley B, et al. Energy sorghum–a genetic model for the design of C4 grass bioenergy crops. J Exp Bot. 2014;65:3479–89.

    Article  PubMed  Google Scholar 

  53. Shih PM, Vuu K, Mansoori N, Ayad L, Louie KB, Bowen BP, et al. A robust gene-stacking method utilizing yeast assembly for plant synthetic biology. Nat Commun. 2016;7:13215.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Gurel S, Gurel E, Miller TI, Lemaux PG. Agrobacterium-mediated transformation of Sorghum bicolor using immature embryos. Methods Mol Biol. 2012;847:109–22.

    Article  CAS  PubMed  Google Scholar 

  55. Hao Z, Yogiswara S, Wei T, Benites VT, Sinha A, Wang G, et al. Expression of a bacterial 3-dehydroshikimate dehydratase (QsuB) reduces lignin and improves biomass saccharification efficiency in switchgrass (Panicum virgatum L.). BMC Plant Biol. 2021;21:56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Sudhakar Reddy P, Srinivas Reddy D, Sivasakthi K, Bhatnagar-Mathur P, Vadez V, Sharma KK. Evaluation of Sorghum [Sorghum bicolor (L.)] Reference genes in various tissues and under abiotic stress conditions for quantitative real-time PCR data normalization. Front Plant Sci. 2016;7:529–36.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Turumtay EA, Turumtay H, Tian Y, Lin C-Y, Chai YN, Louie K, et al. Expression in poplar of dehydroshikimate dehydratase induces transcriptional and metabolic changes in the phenylpropanoid pathway. J Exp Bot. 2024;75(16):4960–77. https://doi.org/10.1093/jxb/erae251.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Yao Y, Sun T, Wang T, Ruebel O, Northen T, Bowen BP. Analysis of metabolomics datasets with high-performance computing and metabolite atlases. Metabolites. 2015;5:431–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Pluskal T, Castillo S, Villar-Briones A, Oresic M. MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics. 2010;11:395.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Wang M, Carver JJ, Phelan VV, Sanchez LM, Garg N, Peng Y, et al. Sharing and community curation of mass spectrometry data with Global Natural Products Social Molecular Networking. Nat Biotechnol. 2016;34:828–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Djoumbou Feunang Y, Eisner R, Knox C, Chepelev L, Hastings J, Owen G, et al. ClassyFire: automated chemical classification with a comprehensive, computable taxonomy. J Cheminform. 2016;8:61.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Kim HW, Wang M, Leber CA, Nothias L-F, Reher R, Kang KB, et al. NPClassifier: A deep meural metwork-based structural classification tool for natural products. J Nat Prod. 2021;84:2795–807.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Tian Y, Lin C-Y, Park J-H, Wu C-Y, Kakumanu R, Pidatala VR, et al. Overexpression of the rice BAHD acyltransferase AT10 increases xylan-bound p-coumarate and reduces lignin in Sorghum bicolor. Biotechnol Biofuels. 2021;14:217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J. Determination of structural carbohydrates and lignin in biomass. In: Laboratory Analytical Procedure. Technical Report, 2008, NREL/TP-510-42618. National Renewable Energy Laboratory, Golden, CO.

  66. Kim H, Ralph J. Solution-state 2D NMR of ball-milled plant cell wall gels in DMSO-d(6)/pyridine-d(5). Org Biomol Chem. 2010;8:576–91.

    Article  CAS  PubMed  Google Scholar 

  67. Mansfield SD, Kim H, Lu F, Ralph J. Whole plant cell wall characterization using solution-state 2D NMR. Nat Protoc. 2012;7:1579–89.

    Article  CAS  PubMed  Google Scholar 

  68. Lan W, Yue F, Rencoret J, Del Río JC, Boerjan W, Lu F, et al. Elucidating tricin-lignin structures: Assigning correlations in HSQC Spectra of monocot lignins. Polymers (Basel). 2018;10(8):916.

    Article  PubMed  Google Scholar 

  69. del Río JC, Rencoret J, Prinsen P, Martínez ÁT, Ralph J, Gutiérrez A. Structural characterization of wheat straw lignin as revealed by analytical pyrolysis, 2D-NMR, and reductive cleavage methods. J Agric Food Chem. 2012;60:5922–35.

    Article  PubMed  Google Scholar 

  70. Choudhary H, Das L, Pelton JG, Sheps L, Simmons BA, Gladden JM, et al. Funneled depolymerization of ionic liquid-based biorefinery “heterogeneous” lignin into guaiacols over reusable palladium catalyst. Chem Eur J. 2023;29: e202300330.

    Article  CAS  PubMed  Google Scholar 

  71. Eudes A, Sathitsuksanoh N, Baidoo EEK, George A, Liang Y, Yang F, et al. Expression of a bacterial 3-dehydroshikimate dehydratase reduces lignin content and improves biomass saccharification efficiency. Plant Biotechnol J. 2015;13:1241–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Lin C-Y, Vuu KM, Amer B, Shih PM, Baidoo EEK, Scheller HV, et al. In-planta production of the biodegradable polyester precursor 2-pyrone-4,6-dicarboxylic acid (PDC): Stacking reduced biomass recalcitrance with value-added co-product. Metab Eng. 2021;66:148–56.

    Article  CAS  PubMed  Google Scholar 

  73. Lin C-Y, Tian Y, Nelson-Vasilchik K, Hague J, Kakumanu R, Lee MY, et al. Engineering sorghum for higher 4-hydroxybenzoic acid content. Metab Eng Commun. 2022;15: e00207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Tian Y, Yang M, Lin C-Y, Park J-H, Wu C-Y, Kakumanu R, et al. Expression of dehydroshikimate dehydratase in sorghum improves biomass yield, accumulation of protocatechuate, and biorefinery economics. ACS Sustain Chem Eng. 2022;10:12520–8.

    Article  CAS  Google Scholar 

  75. Kolberg L, Raudvere U, Kuzmin I, Adler P, Vilo J, Peterson H. g:Profiler-interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update). Nucleic Acids Res. 2023;51:W207–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Valero-Mora PM. ggplot2: elegant graphics for data analysis. J Stat Softw. 2010;35 Book Review 1.

  77. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2018. Available at https://www.R-project.org/.

  78. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 40: A universal enrichment tool for interpreting omics data. Innovation (Camb). 2021;2:100141–5.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are thankful to Stephen Temple for contributing to the sorghum transformation. We thank Novozymes for providing Cellic CTec3 cellulase.

Funding

This material is based upon work supported by the Joint BioEnergy Institute, U.S. Department of Energy, Office of Science, Biological and Environmental Research Program under Award Number DE-AC02-05CH11231 with Lawrence Berkeley National Laboratory. The work (proposal award DOI: https://doi.org/10.46936/10.25585/60008680) conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, was supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy′s National Nuclear Security Administration under contract DE-NA0003525. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript or allow others to do so, for the U.S. Government purpose. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. DOE or the U.S. Government.

Author information

Authors and Affiliations

Authors

Contributions

YT, YG, HC, EAT, YNC, HC, J-HP, C-YW, CDB, JD, KBL, TH, DC, KMV, EEKB, and AE performed experiments and analyzed the data. BPB, PMS, TRN, BAS, RH, JA, DHP, CDS, JCM, HVS, and AE supervised the research. AE designed the research and wrote the paper with contributions from all authors.

Corresponding author

Correspondence to Aymerick Eudes.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

AE is an inventor on a patent application related to this work (US20200354736A1). CDS has a financial interest in Cyklos Materials. Other authors have no competing interests to declare.

Additional information

Publisher's Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tian, Y., Gao, Y., Turumtay, H. et al. Engineered reduction of S-adenosylmethionine alters lignin in sorghum. Biotechnol Biofuels 17, 128 (2024). https://doi.org/10.1186/s13068-024-02572-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-024-02572-8

Keywords