Quanzi Li

Towards a Systems Approach for Lignin Biosynthesis in Populus trichocarpa: Transcript Abundance and Specificity of the Monolignol Biosynthetic Genes

As a step toward a comprehensive description of lignin biosynthesis in Populus trichocarpa, we identified from the genome sequence 95 phenylpropanoid gene models in 10 protein families encoding enzymes for monolignol biosynthesis. Transcript abundance was determined for all 95 genes in xylem, leaf, shoot and phloem using quantitative real-time PCR (qRT-PCR). We identified 23 genes that most probably encode monolignol biosynthesis enzymes during wood formation. Transcripts for 18 of the 23 are abundant and specific to differentiating xylem. We found evidence suggesting functional redundancy at the transcript level for phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate:CoA ligase (4CL), p-hydroxycinnamoyl-CoA:quinate shikimate p-hydroxycinnamoyltransferase (HCT), caffeoyl-CoA O-methyltransferase (CCoAOMT) and coniferyl aldehyde 5-hydroxylase (CAld5H). We carried out an enumeration-based motif identification and discriminant analysis on the promoters of all 95 genes. Five core motifs correctly discriminate the 18 xylem-specific genes from the 77 non-xylem genes. These motifs are similar to promoter elements known to regulate phenylpropanoid gene expression. This work suggests that genes in monolignol biosynthesis are regulated by multiple motifs, often related in sequence.

Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa

Laccases, as early as 1959, were proposed to catalyze the oxidative polymerization of monolignols. Genetic evidence in support of this hypothesis has been elusive due to functional redundancy of laccase genes. An Arabidopsis double mutant demonstrated the involvement of laccases in lignin biosynthesis. We previously identified a subset of laccase genes to be targets of a microRNA (miRNA) ptr-miR397a in Populus trichocarpa. To elucidate the roles of ptr-miR397a and its targets, we characterized the laccase gene family and identified 49 laccase gene models, of which 29 were predicted to be targets of ptr-miR397a. We overexpressed Ptr-MIR397a in transgenic P. trichocarpa. In each of all nine transgenic lines tested, 17 PtrLACs were down-regulated as analyzed by RNA-seq. Transgenic lines with severe reduction in the expression of these laccase genes resulted in an ∼40% decrease in the total laccase activity. Overexpression of Ptr-MIR397a in these transgenic lines also reduced lignin content, whereas levels of all monolignol biosynthetic gene transcripts remained unchanged. A hierarchical genetic regulatory network (GRN) built by a bottom-up graphic Gaussian model algorithm provides additional support for a role of ptr-miR397a as a negative regulator of laccases for lignin biosynthesis. Full transcriptome–based differential gene expression in the overexpressed transgenics and protein domain analyses implicate previously unidentified transcription factors and their targets in an extended hierarchical GRN including ptr-miR397a and laccases that coregulate lignin biosynthesis in wood formation. Ptr-miR397a, laccases, and other regulatory components of this network may provide additional strategies for genetic manipulation of lignin content.

Membrane protein complexes catalyze both 4- and 3-hydroxylation of cinnamic acid derivatives in monolignol biosynthesis

The hydroxylation of 4- and 3-ring carbons of cinnamic acid derivatives during monolignol biosynthesis are key steps that determine the structure and properties of lignin. Individual enzymes have been thought to catalyze these reactions. In stem differentiating xylem (SDX) of Populus trichocarpa, two cinnamic acid 4-hydroxylases (PtrC4H1 and PtrC4H2) and a p-coumaroyl ester 3-hydroxylase (PtrC3H3) are the enzymes involved in these reactions. Here we present evidence that these hydroxylases interact, forming heterodimeric (PtrC4H1/C4H2, PtrC4H1/C3H3, and PtrC4H2/C3H3) and heterotrimeric (PtrC4H1/C4H2/C3H3) membrane protein complexes. Enzyme kinetics using yeast recombinant proteins demonstrated that the enzymatic efficiency (Vmax/km) for any of the complexes is 70–6,500 times greater than that of the individual proteins. The highest increase in efficiency was found for the PtrC4H1/C4H2/C3H3-mediated p-coumaroyl ester 3-hydroxylation. Affinity purification-quantitative mass spectrometry, bimolecular fluorescence complementation, chemical cross-linking, and reciprocal coimmunoprecipitation provide further evidence for these multiprotein complexes. The activities of the recombinant and SDX plant proteins demonstrate two protein-complex–mediated 3-hydroxylation paths in monolignol biosynthesis in P. trichocarpa SDX; one converts p-coumaric acid to caffeic acid and the other converts p-coumaroyl shikimic acid to caffeoyl shikimic acid. Cinnamic acid 4-hydroxylation is also mediated by the same protein complexes. These results provide direct evidence for functional involvement of membrane protein complexes in monolignol biosynthesis.

Complete Proteomic-Based Enzyme Reaction and Inhibition Kinetics Reveal How Monolignol Biosynthetic Enzyme Families Affect Metabolic Flux and Lignin in Populus trichocarpa

A proteomic-based predictive kinetic metabolic-flux model was developed for monolignol biosynthesis in Populus trichocarpa. Absolute quantities of all monolignol pathway proteins and 189 kinetic parameters were generated to construct the model, which was experimentally validated in transgenic P. trichocarpa and provides a comprehensive description of the monolignol biosynthetic pathway. 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.

Splice variant of the SND1 transcription factor is a dominant negative of SND1 members and their regulation in Populus trichocarpa

Secondary Wall-Associated NAC Domain 1s (SND1s) are transcription factors (TFs) known to activate a cascade of TF and pathway genes affecting secondary cell wall biosynthesis (xylogenesis) in Arabidopsis and poplars. Elevated SND1 transcriptional activation leads to ectopic xylogenesis and stunted growth. Nothing is known about the upstream regulators of SND1. Here we report the discovery of a stem-differentiating xylem (SDX)-specific alternative SND1 splice variant, PtrSND1-A2IR, that acts as a dominant negative of SND1 transcriptional network genes in Populus trichocarpa. PtrSND1-A2IR derives from PtrSND1-A2, one of the four fully spliced PtrSND1 gene family members (PtrSND1-A1, -A2, -B1, and -B2). Each full-size PtrSND1 activates its own gene, and all four full-size members activate a common MYB gene (PtrMYB021). PtrSND1-A2IR represses the expression of its PtrSND1 member genes and PtrMYB021. Repression of the autoregulation of a TF family by its only splice variant has not been previously reported in plants. PtrSND1-A2IR lacks DNA binding and transactivation abilities but retains dimerization capability. PtrSND1-A2IR is localized exclusively in cytoplasmic foci. In the presence of any full-size PtrSND1 member, PtrSND1-A2IR is translocated into the nucleus exclusively as a heterodimeric partner with full-size PtrSND1s. Our findings are consistent with a model in which the translocated PtrSND1-A2IR lacking DNA-binding and transactivating abilities can disrupt the function of full-size PtrSND1s, making them nonproductive through heterodimerization, and thereby modulating the SND1 transcriptional network. PtrSND1-A2IR may contribute to transcriptional homeostasis to avoid deleterious effects on xylogenesis and plant growth.

SND1 Transcription Factor–Directed Quantitative Functional Hierarchical Genetic Regulatory Network in Wood Formation in Populus trichocarpa

Novel methods were developed and demonstrated for the discovery of genetic regulatory networks in wood-forming tissues. Transfection of protoplasts from differentiating xylem with the transcription factor gene Ptr-SND1-B1 and novel computational analysis revealed a three-level hierarchical genetic regulatory network that was verified by ChIP and Ptr-SND1-B1 overexpression in transgenic plants. Wood is an essential renewable raw material for industrial products and energy. However, knowledge of the genetic regulation of wood formation is limited. We developed a genome-wide high-throughput system for the discovery and validation of specific transcription factor (TF)–directed hierarchical gene regulatory networks (hGRNs) in wood formation. This system depends on a new robust procedure for isolation and transfection of Populus trichocarpa stem differentiating xylem protoplasts. We overexpressed Secondary Wall-Associated NAC Domain 1s (Ptr-SND1-B1), a TF gene affecting wood formation, in these protoplasts and identified differentially expressed genes by RNA sequencing. Direct Ptr-SND1-B1–DNA interactions were then inferred by integration of time-course RNA sequencing data and top-down Graphical Gaussian Modeling–based algorithms. These Ptr-SND1-B1-DNA interactions were verified to function in differentiating xylem by anti-PtrSND1-B1 antibody-based chromatin immunoprecipitation (97% accuracy) and in stable transgenic P. trichocarpa (90% accuracy). In this way, we established a Ptr-SND1-B1–directed quantitative hGRN involving 76 direct targets, including eight TF and 61 enzyme-coding genes previously unidentified as targets. The network can be extended to the third layer from the second-layer TFs by computation or by overexpression of a second-layer TF to identify a new group of direct targets (third layer). This approach would allow the sequential establishment, one two-layered hGRN at a time, of all layers involved in a more comprehensive hGRN. Our approach may be particularly useful to study hGRNs in complex processes in plant species resistant to stable genetic transformation and where mutants are unavailable.

Plant biotechnology for lignocellulosic biofuel production

Lignocelluloses from plant cell walls are attractive resources for sustainable biofuel production. However, conversion of lignocellulose to biofuel is more expensive than other current technologies, due to the costs of chemical pretreatment and enzyme hydrolysis for cell wall deconstruction. Recalcitrance of cell walls to deconstruction has been reduced in many plant species by modifying plant cell walls through biotechnology. These results have been achieved by reducing lignin content and altering its composition and structure. Reduction of recalcitrance has also been achieved by manipulating hemicellulose biosynthesis and by overexpression of bacterial enzymes in plants to disrupt linkages in the lignin-carbohydrate complexes. These modified plants often have improved saccharification yield and higher ethanol production. Cell wall-degrading (CWD) enzymes from bacteria and fungi have been expressed at high levels in plants to increase the efficiency of saccharification compared with exogenous addition of cellulolytic enzymes. In planta expression of heat-stable CWD enzymes from bacterial thermophiles has made autohydrolysis possible. Transgenic plants can be engineered to reduce recalcitrance without any yield penalty, indicating that successful cell wall modification can be achieved without impacting cell wall integrity or plant development. A more complete understanding of cell wall formation and structure should greatly improve lignocellulosic feedstocks and reduce the cost of biofuel production.

Down-regulation of glycosyltransferase 8D genes in Populus trichocarpa caused reduced mechanical strength and xylan content in wood

Members of glycosyltransferase protein families GT8, GT43 and GT47 are implicated in the biosynthesis of xylan in the secondary cell walls of Arabidopsis. The Arabidopsis mutant irx8 has a 60% reduction in xylan. However, over-expression of an ortholog of Arabidopsis IRX8, poplar PoGT8D, in Arabidopsis irx8 mutant could not restore xylan synthesis. The functions of tree GT8D genes remain unclear. We identified two GT8 gene homologs, PtrGT8D1 and PtrGT8D2, in Populus trichocarpa. They are the only two GT8D members and are abundantly and specifically expressed in the differentiating xylem of P. trichocarpa. PtrGT8D1 transcript abundance was >7 times that of PtrGT8D2. To elucidate the genetic function of GT8D in P. trichocarpa, the expression of PtrGT8D1 and PtrGT8D2 was simultaneously knocked down through RNAi. Four transgenic lines had 85-94% reduction in transcripts of PtrGT8D1 and PtrGT8D2, resulting in 29-36% reduction in stem wood xylan content. Xylan reduction had essentially no effect on cellulose quantity but caused an 11-25% increase in lignin. These transgenics exhibit a brittle wood phenotype, accompanied by increased vessel diameter and thinner fiber cell walls in stem xylem. Stem modulus of elasticity and modulus of rupture were reduced by 17-29% and 16-23%, respectively, and were positively correlated with xylan content but negatively correlated with lignin quantity. These results suggest that PtrGT8Ds play key roles in xylan biosynthesis in wood. Xylan may be a more important factor than lignin affecting the stiffness and fracture strength of wood.

Chemical and spatial differentiation of syringyl and guaiacyl lignins in poplar wood via time-of-flight secondary ion mass spectrometry.

As a major component in plant cell walls, lignin is an important factor in numerous industrial processes, especially in wood saccharification and fermentation to biofuels. The ability to chemically differentiate and spatially locate lignins in wood cell structures provides an important contribution to the effort to improve these processes. The spatial distribution of the syringyl (S) and guaiacyl (G) lignins, both over larger regions and within a single cell wall, on poplar ( Populus trichocarpa ) wood cross-sections was determined via time-of-flight secondary ion mass spectrometry (ToF-SIMS). This is the first time that direct chemically specific mass spectrometric mapping has been employed to elucidate the spatial distribution of S and G lignins. In agreement with results obtained by UV microscopy, ToF-SIMS images clearly show that the guaiacyl lignin is predominantly located in the vessel cell walls of poplar wood while syringyl lignin is mainly located in the fiber cell walls. The G/S ratio in vessel cell walls was determined to be approximately twice that found in fiber cell walls. A combination of Bi ToF-SIMS spectral image acquisition and C(60) sputtering provided the ability to attain the combination of spatial resolution and signal-to-noise necessary to determine the distribution of S and G lignins in a single cell wall. By this technique, it was possible to demonstrate that more guaiacyl lignin is located in the middle lamella layer and more syringyl lignin is located in the inner cell wall area.

Comprehensive Quantification of Monolignol-Pathway Enzymes in Populus trichocarpa by Protein Cleavage Isotope Dilution Mass Spectrometry

The economic value of wood/pulp from many tree species is largely dictated by the quantity and chemical properties of lignin, which is directly related to the composition and linkages of monolignols comprising the polymer. Although much is known regarding the monolignol biosynthetic pathway, our understanding is still deficient due to the lack of quantitative information at the proteomic level. We developed an assay based on protein cleavage isotope dilution mass spectrometry (PC-IDMS) for the determination of all potential, primary enzymes involved in the biosynthesis of monolignols and the peroxidases responsible for their polymerization to form lignin in the model tree species, Populus trichocarpa. Described is the identification of quantitative surrogate peptides through shotgun analysis of native and recombinant proteins, optimization of trypsin proteolysis using fractional factorial design of experiments, and development of a liquid chromatography-selected reaction monitoring method for specific detection of all targeted peptides. Of the 25 targeted enzymes, three were undetected in the normal xylem tissues, and all but two of the detectable species showed good day-to-day precision (CV < 10%). This represents the most comprehensive assay for quantification of proteins regulating monolignol biosynthesis and will lead to a better understanding of lignin formation at a systems level.