Modeling cross-regulatory influences on monolignol transcripts and proteins under single and combinatorial gene knockdowns in Populus trichocarpa.

Accurate manipulation of metabolites in monolignol biosynthesis is a key step for controlling lignin content, structure, and other wood properties important to the bioenergy and biomaterial industries. A crucial component of this strategy is predicting how single and combinatorial knockdowns of monolignol specific gene transcripts influence the abundance of monolignol proteins, which are the driving mechanisms of monolignol biosynthesis. Computational models have been developed to estimate protein abundances from transcript perturbations of monolignol specific genes. The accuracy of these models, however, is hindered by their inability to capture indirect regulatory influences on other pathway genes. Here, we examine the manifestation of these indirect influences on transgenic transcript and protein abundances, identifying putative indirect regulatory influences that occur when one or more specific monolignol pathway genes are perturbed. We created a computational model using sparse maximum likelihood to estimate the resulting monolignol transcript and protein abundances in transgenic Populus trichocarpa based on targeted knockdowns of specific monolignol genes. Using in-silico simulations of this model and root mean square error, we showed that our model more accurately estimated transcript and protein abundances, in comparison to previous models, when individual and families of monolignol genes were perturbed. We leveraged insight from the inferred network structure obtained from our model to identify potential genes, including PtrHCT, PtrCAD, and Ptr4CL, involved in post-transcriptional and/or post-translational regulation. Our model provides a useful computational tool for exploring the cascaded impact of single and combinatorial modifications of monolignol specific genes on lignin and other wood properties.

Complete proteomic-based enzyme reaction and inhibition kinetics reveal how monolignol biosynthetic enzyme families affect metabolic flux and lignin in Populus trichocarpa.

We established a predictive kinetic metabolic-flux model for the 21 enzymes and 24 metabolites of the monolignol biosynthetic pathway using Populus trichocarpa secondary differentiating xylem. To establish this model, a comprehensive study was performed to obtain the reaction and inhibition kinetic parameters of all 21 enzymes based on functional recombinant proteins. A total of 104 Michaelis-Menten kinetic parameters and 85 inhibition kinetic parameters were derived from these enzymes. Through mass spectrometry, we obtained the absolute quantities of all 21 pathway enzymes in the secondary differentiating xylem. This extensive experimental data set, generated from a single tissue specialized in wood formation, was used to construct the predictive kinetic metabolic-flux model to provide a comprehensive mathematical description of the monolignol biosynthetic pathway. The model was validated using experimental data from transgenic P. trichocarpa plants. The model predicts how pathway enzymes affect lignin content and composition, explains a long-standing paradox regarding the regulation of monolignol subunit ratios in lignin, and reveals novel mechanisms involved in the regulation of lignin biosynthesis. This model provides an explanation of the effects of genetic and transgenic perturbations of the monolignol biosynthetic pathway in flowering plants.

Monolignol pathway 4-coumaric acid:coenzyme A ligases in Populus trichocarpa: novel specificity, metabolic regulation, and simulation of coenzyme A ligation fluxes.

4-Coumaric acid:coenzyme A ligase (4CL) is involved in monolignol biosynthesis for lignification in plant cell walls. It ligates coenzyme A (CoA) with hydroxycinnamic acids, such as 4-coumaric and caffeic acids, into hydroxycinnamoyl-CoA thioesters. The ligation ensures the activated state of the acid for reduction into monolignols. In Populus spp., it has long been thought that one monolignol-specific 4CL is involved. Here, we present evidence of two monolignol 4CLs, Ptr4CL3 and Ptr4CL5, in Populus trichocarpa. Ptr4CL3 is the ortholog of the monolignol 4CL reported for many other species. Ptr4CL5 is novel. The two Ptr4CLs exhibited distinct Michaelis-Menten kinetic properties. Inhibition kinetics demonstrated that hydroxycinnamic acid substrates are also inhibitors of 4CL and suggested that Ptr4CL5 is an allosteric enzyme. Experimentally validated flux simulation, incorporating reaction/inhibition kinetics, suggested two CoA ligation paths in vivo: one through 4-coumaric acid and the other through caffeic acid. We previously showed that a membrane protein complex mediated the 3-hydroxylation of 4-coumaric acid to caffeic acid. The demonstration here of two ligation paths requiring these acids supports this 3-hydroxylation function. Ptr4CL3 regulates both CoA ligation paths with similar efficiencies, whereas Ptr4CL5 regulates primarily the caffeic acid path. Both paths can be inhibited by caffeic acid. The Ptr4CL5-catalyzed caffeic acid metabolism, therefore, may also act to mitigate the inhibition by caffeic acid to maintain a proper ligation flux. A high level of caffeic acid was detected in stem-differentiating xylem of P. trichocarpa. Our results suggest that Ptr4CL5 and caffeic acid coordinately modulate the CoA ligation flux for monolignol biosynthesis.

Multiplex CRISPR editing of wood for sustainable fiber production.

The domestication of forest trees for a more sustainable fiber bioeconomy has long been hindered by the complexity and plasticity of lignin, a biopolymer in wood that is recalcitrant to chemical and enzymatic degradation. Here, we show that multiplex CRISPR editing enables precise woody feedstock design for combinatorial improvement of lignin composition and wood properties. By assessing every possible combination of 69,123 multigenic editing strategies for 21 lignin biosynthesis genes, we deduced seven different genome editing strategies targeting the concurrent alteration of up to six genes and produced 174 edited poplar variants. CRISPR editing increased the wood carbohydrate-to-lignin ratio up to 228% that of wild type, leading to more-efficient fiber pulping. The edited wood alleviates a major fiber-production bottleneck regardless of changes in tree growth rate and could bring unprecedented operational efficiencies, bioeconomic opportunities, and environmental benefits.

Functional redundancy of the two 5-hydroxylases in monolignol biosynthesis of Populus trichocarpa: LC-MS/MS based protein quantification and metabolic flux analysis.

Flowering plants have syringyl and guaiacyl subunits in lignin in contrast to the guaiacyl lignin in gymnosperms. The biosynthesis of syringyl subunits is initiated by coniferaldehyde 5-hydroxylase (CAld5H). In Populus trichocarpa there are two closely related CAld5H enzymes (PtrCAld5H1 and PtrCAld5H2) associated with lignin biosynthesis during wood formation. We used yeast recombinant PtrCAld5H1 and PtrCAld5H2 proteins to carry out Michaelis-Menten and inhibition kinetics with LC-MS/MS based absolute protein quantification. CAld5H, a monooxygenase, requires a cytochrome P450 reductase (CPR) as an electron donor. We cloned and expressed three P. trichocarpa CPRs in yeast and show that all are active with both CAld5Hs. The kinetic analysis shows both CAld5Hs have essentially the same biochemical functions. When both CAld5Hs are coexpressed in the same yeast membranes, the resulting enzyme activities are additive, suggesting functional redundancy and independence of these two enzymes. Simulated reaction flux based on Michaelis-Menten kinetics and inhibition kinetics confirmed the redundancy and independence. Subcellular localization of both CAld5Hs as sGFP fusion proteins expressed in P. trichocarpa differentiating xylem protoplasts indicate that they are endoplasmic reticulum resident proteins. These results imply that during wood formation, 5-hydroxylation in monolignol biosynthesis of P. trichocarpa requires the combined metabolic flux of these two CAld5Hs to maintain adequate biosynthesis of syringyl lignin. The combination of genetic analysis, absolute protein quantitation-based enzyme kinetics, homologous CPR specificity, SNP characterization, and ER localization provides a more rigorous basis for a comprehensive systems understanding of 5-hydroxylation in lignin biosynthesis.

Flux modeling for monolignol biosynthesis.

The pathway of monolignol biosynthesis involves many components interacting in a metabolic grid to regulate the supply and ratios of monolignols for lignification. The complexity of the pathway challenges any intuitive prediction of the output without mathematical modeling. Several models have been presented to quantify the metabolic flux for monolignol biosynthesis and the regulation of lignin content, composition, and structure in plant cell walls. Constraint-based models using data from transgenic plants were formulated to describe steady-state flux distribution in the pathway. Kinetic-based models using enzyme reaction and inhibition constants were developed to predict flux dynamics for monolignol biosynthesis in wood-forming cells. This review summarizes the recent progress in flux modeling and its application to lignin engineering for improved plant development and utilization.

Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis.

A multi-omics quantitative integrative analysis of lignin biosynthesis can advance the strategic engineering of wood for timber, pulp, and biofuels. Lignin is polymerized from three monomers (monolignols) produced by a grid-like pathway. The pathway in wood formation of Populus trichocarpa has at least 21 genes, encoding enzymes that mediate 37 reactions on 24 metabolites, leading to lignin and affecting wood properties. We perturb these 21 pathway genes and integrate transcriptomic, proteomic, fluxomic and phenomic data from 221 lines selected from ~2000 transgenics (6-month-old). The integrative analysis estimates how changing expression of pathway gene or gene combination affects protein abundance, metabolic-flux, metabolite concentrations, and 25 wood traits, including lignin, tree-growth, density, strength, and saccharification. The analysis then predicts improvements in any of these 25 traits individually or in combinations, through engineering expression of specific monolignol genes. The analysis may lead to greater understanding of other pathways for improved growth and adaptation.