Ectopic Expression of PtrLBD39 Retarded Primary and Secondary Growth in Populus trichocarpa.

Primary and secondary growth of trees are needed for increments in plant height and stem diameter, respectively, affecting the production of woody biomass for applications in timber, pulp/paper, and related biomaterials. These two types of growth are believed to be both regulated by distinct transcription factor (TF)-mediated regulatory pathways. Notably, we identified PtrLBD39, a highly stem phloem-specific TF in Populus trichocarpa and found that the ectopic expression of PtrLBD39 in P. trichocarpa markedly retarded both primary and secondary growth. In these overexpressing plants, the RNA-seq, ChIP-seq, and weighted gene co-expression network analysis (WGCNA) revealed that PtrLBD39 directly or indirectly regulates TFs governing vascular tissue development, wood formation, hormonal signaling pathways, and enzymes responsible for wood components. This regulation led to growth inhibition, decreased fibrocyte secondary cell wall thickness, and reduced wood production. Therefore, our study indicates that, following ectopic expression in P. trichocarpa, PtrLBD39 functions as a repressor influencing both primary and secondary growth.

PtrVCS2 Regulates Drought Resistance by Changing Vessel Morphology and Stomatal Closure in Populus trichocarpa.

Drought has severe effects on plant growth, forest productivity, and survival throughout the world. Understanding the molecular regulation of drought resistance in forest trees can enable effective strategic engineering of novel drought-resistant genotypes of tree species. In this study, we identified a gene, PtrVCS2, encoding a zinc finger (ZF) protein of the ZF-homeodomain transcription factor in Populus trichocarpa (Black Cottonwood) Torr. & A. Gray. ex Hook. Overexpression of PtrVCS2 (OE-PtrVCS2) in P. trichocarpa resulted in reduced growth, a higher proportion of smaller stem vessels, and strong drought-resistance phenotypes. Stomatal movement experiments revealed that the OE-PtrVCS2 transgenics showed lower stomata apertures than wild-type plants under drought conditions. RNA-seq analysis of the OE-PtrVCS2 transgenics showed that PtrVCS2 regulates the expression of multiple genes involved in regulation of stomatal opening and closing, particularly the PtrSULTR3;1-1 gene, and several genes related to cell wall biosynthesis, such as PtrFLA11-12 and PtrPR3-3. Moreover, we found that the water use efficiency of the OE-PtrVCS2 transgenic plants was consistently higher than that of wild type plants when subjected to chronic drought stress. Taken together, our results suggest that PtrVCS2 plays a positive role in improving drought adaptability and resistance in P. trichocarpa.

A novel synthetic-genetic-array-based yeast one-hybrid system for high discovery rate and short processing time.

Eukaryotic gene expression is often tightly regulated by interactions between transcription factors (TFs) and their DNA cis targets. Yeast one-hybrid (Y1H) is one of the most extensively used methods to discover these interactions. We developed a high-throughput meiosis-directed yeast one-hybrid system using the Magic Markers of the synthetic genetic array analysis. The system has a transcription factor-DNA interaction discovery rate twice as high as the conventional diploid-mating approach and a processing time nearly one-tenth of the haploid-transformation method. The system also offers the highest accuracy in identifying TF-DNA interactions that can be authenticated in vivo by chromatin immunoprecipitation. With these unique features, this meiosis-directed Y1H system is particularly suited for constructing novel and comprehensive genome-scale gene regulatory networks for various organisms.

Dimerization of PtrMYB074 and PtrWRKY19 mediates transcriptional activation of PtrbHLH186 for secondary xylem development in Populus trichocarpa.

Wood formation is controlled by transcriptional regulatory networks (TRNs) involving regulatory homeostasis determined by combinations of transcription factor (TF)-DNA and TF-TF interactions. Functions of TF-TF interactions in wood formation are still in the early stages of identification. PtrMYB074 is a woody dicot-specific TF in a TRN for wood formation in Populus trichocarpa. Here, using yeast two-hybrid and bimolecular fluorescence complementation, we conducted a genome-wide screening for PtrMYB074 interactors and identified 54 PtrMYB074-TF pairs. Of these pairs, 53 are novel. We focused on the PtrMYB074-PtrWRKY19 pair, the most highly expressed and xylem-specific interactor, and its direct transregulatory target, PtrbHLH186, the xylem-specific one of the pair's only two direct TF target genes. Using transient and CRISPR-mediated transgenesis in P. trichocarpa coupled with chromatin immunoprecipitation and electrophoretic mobility shift assays, we demonstrated that PtrMYB074 is recruited by PtrWRKY19 and that the PtrMYB074-PtrWRKY19 dimers are required to transactive PtrbHLH186. Overexpressing PtrbHLH186 in P. trichocarpa resulted in retarded plant growth, increased guaiacyl lignin, a higher proportion of smaller stem vessels and strong drought-tolerant phenotypes. Knowledge of the PtrMYB074-PtrWRKY19-PtrbHLH186 regulation may help design genetic controls of optimal growth and wood formation to maximize beneficial wood properties while minimizing negative effects on growth.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood.

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

Hierarchical Transcription Factor and Chromatin Binding Network for Wood Formation in Black Cottonwood (Populus trichocarpa).

Wood remains the world's most abundant and renewable resource for timber and pulp and is an alternative to fossil fuels. Understanding the molecular regulation of wood formation can advance the engineering of wood for more efficient material and energy productions. We integrated a black cottonwood (Populus trichocarpa) wood-forming cell system with quantitative transcriptomics and chromatin binding assays to construct a transcriptional regulatory network (TRN) directed by a key transcription factor (TF), PtrSND1-B1 (secondary wall-associated NAC-domain protein). The network consists of four layers of TF-target gene interactions with quantitative regulatory effects, describing the specificity of how the regulation is transduced through these interactions to activate cell wall genes (effector genes) for wood formation. PtrSND1-B1 directs 57 TF-DNA interactions through 17 TFs transregulating 27 effector genes. Of the 57 interactions, 55 are novel. We tested 42 of these 57 interactions in 30 genotypes of transgenic P. trichocarpa and verified that ∼90% of the tested interactions function in vivo. The TRN reveals common transregulatory targets for distinct TFs, leading to the discovery of nine TF protein complexes (dimers and trimers) implicated in regulating the biosynthesis of specific types of lignin. Our work suggests that wood formation may involve regulatory homeostasis determined by combinations of TF-DNA and TF-TF (protein-protein) regulations.

The AREB1 Transcription Factor Influences Histone Acetylation to Regulate Drought Responses and Tolerance in Populus trichocarpa.

Plants develop tolerance to drought by activating genes with altered levels of epigenetic modifications. Specific transcription factors are involved in this activation, but the molecular connections within the regulatory system are unclear. Here, we analyzed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment and examined its association with transcriptomes in Populus trichocarpa under drought stress. We revealed that abscisic acid-Responsive Element (ABRE) motifs in promoters of the drought-responsive genes PtrNAC006, PtrNAC007, and PtrNAC120 are involved in H3K9ac enhancement and activation of these genes. Overexpressing these PtrNAC genes in P trichocarpa resulted in strong drought-tolerance phenotypes. We showed that the ABRE binding protein PtrAREB1-2 binds to ABRE motifs associated with these PtrNAC genes and recruits the histone acetyltransferase unit ADA2b-GCN5, forming AREB1-ADA2b-GCN5 ternary protein complexes. Moreover, this recruitment enables GCN5-mediated histone acetylation to enhance H3K9ac and enrich RNA polymerase II specifically at these PtrNAC genes for the development of drought tolerance. CRISPR editing or RNA interference-mediated downregulation of any of the ternary members results in highly drought-sensitive P trichocarpa Thus, the combinatorial function of the ternary proteins establishes a coordinated histone acetylation and transcription factor-mediated gene activation for drought response and tolerance in Populus species.

CRISPR-Cas9 editing of CAFFEOYL SHIKIMATE ESTERASE 1 and 2 shows their importance and partial redundancy in lignification in Populus tremula × P. alba.

Lignins are cell wall-located aromatic polymers that provide strength and hydrophobicity to woody tissues. Lignin monomers are synthesized via the phenylpropanoid pathway, wherein CAFFEOYL SHIKIMATE ESTERASE (CSE) converts caffeoyl shikimate into caffeic acid. Here, we explored the role of the two CSE homologs in poplar (Populus tremula × P. alba). Reporter lines showed that the expression conferred by both CSE1 and CSE2 promoters is similar. CRISPR-Cas9-generated cse1 and cse2 single mutants had a wild-type lignin level. Nevertheless, CSE1 and CSE2 are not completely redundant, as both single mutants accumulated caffeoyl shikimate. In contrast, the cse1 cse2 double mutants had a 35% reduction in lignin and associated growth penalty. The reduced-lignin content translated into a fourfold increase in cellulose-to-glucose conversion upon limited saccharification. Phenolic profiling of the double mutants revealed large metabolic shifts, including an accumulation of p-coumaroyl, 5-hydroxyferuloyl, feruloyl and sinapoyl shikimate, in addition to caffeoyl shikimate. This indicates that the CSEs have a broad substrate specificity, which was confirmed by in vitro enzyme kinetics. Taken together, our results suggest an alternative path within the phenylpropanoid pathway at the level of the hydroxycinnamoyl-shikimates, and show that CSE is a promising target to improve plants for the biorefinery.

The microRNA476a-RFL module regulates adventitious root formation through a mitochondria-dependent pathway in Populus.

For woody plants, clonal propagation efficiency is largely determined by adventitious root (AR) formation at the bases of stem cuttings. However, our understanding of the molecular mechanisms contributing to AR morphogenesis in trees remains limited, despite the importance of vegetative propagation, currently the most common practice for tree breeding and commercialization. Here, we identified Populus-specific miR476a as a regulator of wound-induced adventitious rooting that acts by orchestrating mitochondrial homeostasis. MiR476a exhibited inducible expression during AR formation and directly targeted several Restorer of Fertility like (RFL) genes encoding mitochondrion-localized pentatricopeptide repeat proteins. Genetic modification of miR476a-RFL expression revealed that miR476a/RFL-mediated dynamic regulation of mitochondrial homeostasis influences AR formation in poplar. Mitochondrial perturbation via exogenous application of a chemical inhibitor indicated that miR476a/RFL-directed AR formation depends on mitochondrial regulation that acts via auxin signaling. Our results thus establish a microRNA-directed mitochondrion-auxin signaling cascade required for AR development, providing insights into the role of mitochondrial regulation in the developmental plasticity of plants.

MYB Transcription Factor161 Mediates Feedback Regulation of Secondary wall-associated NAC-Domain1 Family Genes for Wood Formation.

Wood formation is a complex process that involves cell differentiation, cell expansion, secondary wall deposition, and programmed cell death. We constructed a four-layer wood formation transcriptional regulatory network (TRN) in Populus trichocarpa (black cottonwood) that has four Secondary wall-associated NAC-Domain1 (PtrSND1) transcription factor (TF) family members as the top-layer regulators. We characterized the function of a MYB (PtrMYB161) TF in this PtrSND1-TRN, using transgenic P trichocarpa cells and whole plants. PtrMYB161 is a third-layer regulator that directly transactivates five wood formation genes. Overexpression of PtrMYB161 in P. trichocarpa (OE-PtrMYB161) led to reduced wood, altered cell type proportions, and inhibited growth. Integrative analysis of wood cell-based chromatin-binding assays with OE-PtrMYB161 transcriptomics revealed a feedback regulation system in the PtrSND1-TRN, where PtrMYB161 represses all four top-layer regulators and one second-layer regulator, PtrMYB021, possibly affecting many downstream TFs in, and likely beyond, the TRN, to generate the observed phenotypic changes. Our data also suggested that the PtrMYB161's repressor function operates through interaction of the base PtrMYB161 target-binding system with gene-silencing cofactors. PtrMYB161 protein does not contain any known negative regulatory domains. CRISPR-based mutants of PtrMYB161 in P. trichocarpa exhibited phenotypes similar to the wild type, suggesting that PtrMYB161's activator functions are redundant among many TFs. Our work demonstrated that PtrMYB161 binds to multiple sets of target genes, a feature that allows it to function as an activator as well as a repressor. The balance of the two functions may be important to the establishment of regulatory homeostasis for normal growth and development.

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.

MYB-Mediated Regulation of Anthocyanin Biosynthesis.

Anthocyanins are natural water-soluble pigments that are important in plants because they endow a variety of colors to vegetative tissues and reproductive plant organs, mainly ranging from red to purple and blue. The colors regulated by anthocyanins give plants different visual effects through different biosynthetic pathways that provide pigmentation for flowers, fruits and seeds to attract pollinators and seed dispersers. The biosynthesis of anthocyanins is genetically determined by structural and regulatory genes. MYB (v-myb avian myeloblastosis viral oncogene homolog) proteins are important transcriptional regulators that play important roles in the regulation of plant secondary metabolism. MYB transcription factors (TFs) occupy a dominant position in the regulatory network of anthocyanin biosynthesis. The TF conserved binding motifs can be combined with other TFs to regulate the enrichment and sedimentation of anthocyanins. In this study, the regulation of anthocyanin biosynthetic mechanisms of MYB-TFs are discussed. The role of the environment in the control of the anthocyanin biosynthesis network is summarized, the complex formation of anthocyanins and the mechanism of environment-induced anthocyanin synthesis are analyzed. Some prospects for MYB-TF to modulate the comprehensive regulation of anthocyanins are put forward, to provide a more relevant basis for further research in this field, and to guide the directed genetic modification of anthocyanins for the improvement of crops for food quality, nutrition and human health.

Molecular and Metabolic Insights into Anthocyanin Biosynthesis for Leaf Color Change in Chokecherry (Padus virginiana).

Chokecherry (Padus virginiana L.) is an important landscaping tree with high ornamental value because of its colorful purplish-red leaves (PRL). The quantifications of anthocyanins and the mechanisms of leaf color change in this species remain unknown. The potential biosynthetic and regulatory mechanisms and the accumulation patterns of anthocyanins in P. virginiana that determine three leaf colors were investigated by combined analysis of the transcriptome and the metabolome. The difference of chlorophyll, carotenoid and anthocyanin content correlated with the formation of P. virginiana leaf color. Using enrichment and correlation network analysis, we found that anthocyanin accumulation differed in different colored leaves and that the accumulation of malvidin 3-O-glucoside (violet) and pelargonidin 3-O-glucoside (orange-red) significantly correlated with the leaf color change from green to purple-red. The flavonoid biosynthesis genes (PAL, CHS and CHI) and their transcriptional regulators (MYB, HD-Zip and bHLH) exhibited specific increased expression during the purple-red periods. Two genes encoding enzymes in the anthocyanin biosynthetic pathway, UDP glucose-flavonoid 3-O-glucosyl-transferase (UFGT) and anthocyanidin 3-O-glucosyltransferase (BZ1), seem to be critical for suppressing the formation of the aforesaid anthocyanins. In PRL, the expression of the genes encoding for UGFT and BZ1 enzymes was substantially higher than in leaves of other colors and may be related with the purple-red color change. These results may facilitate genetic modification or selection for further improvement in ornamental qualities of P. virginiana.

Reciprocal cross-regulation of VND and SND multigene TF families for wood formation in Populus trichocarpa.

Secondary cell wall (SCW) biosynthesis is the biological process that generates wood, an important renewable feedstock for materials and energy. NAC domain transcription factors, particularly Vascular-Related NAC-Domain (VND) and Secondary Wall-Associated NAC Domain (SND) proteins, are known to regulate SCW differentiation. The regulation of VND and SND is important to maintain homeostasis for plants to avoid abnormal growth and development. We previously identified a splice variant, PtrSND1-A2IR , derived from PtrSND1-A2 as a dominant-negative regulator, which suppresses the transactivation of all PtrSND1 family members. PtrSND1-A2IR also suppresses the self-activation of the PtrSND1 family members except for its cognate transcription factor, PtrSND1-A2, suggesting the existence of an unknown factor needed to regulate PtrSND1-A2 Here, a splice variant, PtrVND6-C1IR , derived from PtrVND6-C1 was discovered that suppresses the protein functions of all PtrVND6 family members. PtrVND6-C1IR also suppresses the expression of all PtrSND1 members, including PtrSND1-A2, demonstrating that PtrVND6-C1IR is the previously unidentified regulator of PtrSND1-A2 We also found that PtrVND6-C1IR cannot suppress the expression of its cognate transcription factor, PtrVND6-C1PtrVND6-C1 is suppressed by PtrSND1-A2IR Both PtrVND6-C1IR and PtrSND1-A2IR cannot suppress their cognate transcription factors but can suppress all members of the other family. The results indicate that the splice variants from the PtrVND6 and PtrSND1 family may exert reciprocal cross-regulation for complete transcriptional regulation of these two families in wood formation. This reciprocal cross-regulation between families suggests a general mechanism among NAC domain proteins and likely other transcription factors, where intron-retained splice variants provide an additional level of regulation.

CAD1 and CCR2 protein complex formation in monolignol biosynthesis in Populus trichocarpa.

Lignin is the major phenolic polymer in plant secondary cell walls and is polymerized from monomeric subunits, the monolignols. Eleven enzyme families are implicated in monolignol biosynthesis. Here, we studied the functions of members of the cinnamyl alcohol dehydrogenase (CAD) and cinnamoyl-CoA reductase (CCR) families in wood formation in Populus trichocarpa, including the regulatory effects of their transcripts and protein activities on monolignol biosynthesis. Enzyme activity assays from stem-differentiating xylem (SDX) proteins showed that RNAi suppression of PtrCAD1 in P. trichocarpa transgenics caused a reduction in SDX CCR activity. RNAi suppression of PtrCCR2, the only CCR member highly expressed in SDX, caused a reciprocal reduction in SDX protein CAD activities. The enzyme assays of mixed and coexpressed recombinant proteins supported physical interactions between PtrCAD1 and PtrCCR2. Biomolecular fluorescence complementation and pull-down/co-immunoprecipitation experiments supported a hypothesis of PtrCAD1/PtrCCR2 heterodimer formation. These results provide evidence for the formation of PtrCAD1/PtrCCR2 protein complexes in monolignol biosynthesis in planta.

Phosphorylation is an on/off switch for 5-hydroxyconiferaldehyde O-methyltransferase activity in poplar monolignol biosynthesis.

Although phosphorylation has long been known to be an important regulatory modification of proteins, no unequivocal evidence has been presented to show functional control by phosphorylation for the plant monolignol biosynthetic pathway. Here, we present the discovery of phosphorylation-mediated on/off regulation of enzyme activity for 5-hydroxyconiferaldehyde O-methyltransferase 2 (PtrAldOMT2), an enzyme central to monolignol biosynthesis for lignification in stem-differentiating xylem (SDX) of Populus trichocarpa. Phosphorylation turned off the PtrAldOMT2 activity, as demonstrated in vitro by using purified phosphorylated and unphosphorylated recombinant PtrAldOMT2. Protein extracts of P. trichocarpa SDX, which contains endogenous kinases, also phosphorylated recombinant PtrAldOMT2 and turned off the recombinant protein activity. Similarly, ATP/Mn(2+)-activated phosphorylation of SDX protein extracts reduced the endogenous SDX PtrAldOMT2 activity by ∼ 60%, and dephosphorylation fully restored the activity. Global shotgun proteomic analysis of phosphopeptide-enriched P. trichocarpa SDX protein fractions identified PtrAldOMT2 monophosphorylation at Ser(123) or Ser(125) in vivo. Phosphorylation-site mutagenesis verified the PtrAldOMT2 phosphorylation at Ser(123) or Ser(125) and confirmed the functional importance of these phosphorylation sites for O-methyltransferase activity. The PtrAldOMT2 Ser(123) phosphorylation site is conserved across 93% of AldOMTs from 46 diverse plant species, and 98% of the AldOMTs have either Ser(123) or Ser(125). PtrAldOMT2 is a homodimeric cytosolic enzyme expressed more abundantly in syringyl lignin-rich fiber cells than in guaiacyl lignin-rich vessel cells. The reversible phosphorylation of PtrAldOMT2 is likely to have an important role in regulating syringyl monolignol biosynthesis of P. trichocarpa.

Tissue and cell-type co-expression networks of transcription factors and wood component genes in Populus trichocarpa.

MAIN CONCLUSION: Co-expression networks based on transcriptomes of Populus trichocarpa major tissues and specific cell types suggest redundant control of cell wall component biosynthetic genes by transcription factors in wood formation. We analyzed the transcriptomes of five tissues (xylem, phloem, shoot, leaf, and root) and two wood forming cell types (fiber and vessel) of Populus trichocarpa to assemble gene co-expression subnetworks associated with wood formation. We identified 165 transcription factors (TFs) that showed xylem-, fiber-, and vessel-specific expression. Of these 165 TFs, 101 co-expressed (correlation coefficient, r > 0.7) with the 45 secondary cell wall cellulose, hemicellulose, and lignin biosynthetic genes. Each cell wall component gene co-expressed on average with 34 TFs, suggesting redundant control of the cell wall component gene expression. Co-expression analysis showed that the 101 TFs and the 45 cell wall component genes each has two distinct groups (groups 1 and 2), based on their co-expression patterns. The group 1 TFs (44 members) are predominantly xylem and fiber specific, and are all highly positively co-expressed with the group 1 cell wall component genes (30 members), suggesting their roles as major wood formation regulators. Group 1 TFs include a lateral organ boundary domain gene (LBD) that has the highest number of positively correlated cell wall component genes (36) and TFs (47). The group 2 TFs have 57 members, including 14 vessel-specific TFs, and are generally less correlated with the cell wall component genes. An exception is a vessel-specific basic helix-loop-helix (bHLH) gene that negatively correlates with 20 cell wall component genes, and may function as a key transcriptional suppressor. The co-expression networks revealed here suggest a well-structured transcriptional homeostasis for cell wall component biosynthesis during wood formation.

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.

Systems biology of lignin biosynthesis in Populus trichocarpa: heteromeric 4-coumaric acid:coenzyme A ligase protein complex formation, regulation, and numerical modeling.

As a step toward predictive modeling of flux through the pathway of monolignol biosynthesis in stem differentiating xylem of Populus trichocarpa, we discovered that the two 4-coumaric acid:CoA ligase (4CL) isoforms, 4CL3 and 4CL5, interact in vivo and in vitro to form a heterotetrameric protein complex. This conclusion is based on laser microdissection, coimmunoprecipitation, chemical cross-linking, bimolecular fluorescence complementation, and mass spectrometry. The tetramer is composed of three subunits of 4CL3 and one of 4CL5. 4CL5 appears to have a regulatory role. This protein-protein interaction affects the direction and rate of metabolic flux for monolignol biosynthesis in P. trichocarpa. A mathematical model was developed for the behavior of 4CL3 and 4CL5 individually and in mixtures that form the enzyme complex. The model incorporates effects of mixtures of multiple hydroxycinnamic acid substrates, competitive inhibition, uncompetitive inhibition, and self-inhibition, along with characteristic of the substrates, the enzyme isoforms, and the tetrameric complex. Kinetic analysis of different ratios of the enzyme isoforms shows both inhibition and activation components, which are explained by the mathematical model and provide insight into the regulation of metabolic flux for monolignol biosynthesis by protein complex formation.

SND1 transcription factor-directed quantitative functional hierarchical genetic regulatory network in wood formation in Populus trichocarpa.

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.