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.

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.

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.

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.

Stem vascular cambium cells in forest trees produce wood for materials and energy. WOX4 affects the proliferation of such cells in Populus. Here we show that PtrWOX4a is the most highly expressed stem vascular-cambium-specific (VCS) gene in P. trichocarpa, and its expression is controlled by the product of the second most highly expressed VCS gene, PtrVCS2, encoding a zinc finger protein. PtrVCS2 binds to the PtrWOX4a promoter as part of a PtrWOX13a–PtrVCS2–PtrGCN5-1–PtrADA2b-3 protein tetramer. PtrVCS2 prevented the interaction between PtrGCN5-1 and PtrADA2b-3, resulting in H3K9, H3K14 and H3K27 hypoacetylation at the PtrWOX4a promoter, which led to fewer cambium cell layers. These effects on cambium cell proliferation were consistent across more than 20 sets of transgenic lines overexpressing individual genes, gene-edited mutants and RNA interference lines in P. trichocarpa. We propose that the tetramer–PtrWOX4a system may coordinate genetic and epigenetic regulation to maintain normal vascular cambium development for wood formation.

A MYB transcription factor mediates feedback mechanisms in a hierarchical transcriptional regulatory network for wood formation in Populus species.