PagJAZ5 regulates cambium activity through coordinately modulating cytokinin concentration and signaling in poplar.

Wood is resulted from the radial growth paced by the division and differentiation of vascular cambium cells in woody plants, and phytohormones play important roles in cambium activity. Here, we identified that PagJAZ5, a key negative regulator of jasmonate (JA) signaling, plays important roles in enhancing cambium cell division and differentiation by mediating cytokinin signaling in poplar 84K (Populus alba × Populus glandulosa). PagJAZ5 is preferentially expressed in developing phloem and cambium, weakly in developing xylem cells. Overexpression (OE) of PagJAZ5m (insensitive to JA) increased cambium activity and xylem differentiation, while jaz mutants showed opposite results. Transcriptome analyses revealed that cytokinin oxidase/dehydrogenase (CKXs) and type-A response regulators (RRs) were downregulated in PagJAZ5m OE plants. The bioactive cytokinins were significantly increased in PagJAZ5m overexpressing plants and decreased in jaz5 mutants, compared with that in 84K plants. The PagJAZ5 directly interact with PagMYC2a/b and PagWOX4b. Further, we found that the PagRR5 is regulated by PagMYC2a and PagWOX4b and involved in the regulation of xylem development. Our results showed that PagJAZ5 can increase cambium activity and promote xylem differentiation through modulating cytokinin level and type-A RR during wood formation in poplar.

PagERF81 regulates lignin biosynthesis and xylem cell differentiation in poplar.

Lignin is a major component of plant cell walls and is essential for plant growth and development. Lignin biosynthesis is controlled by a hierarchical regulatory network involving multiple transcription factors. In this study, we showed that the gene encoding an APETALA 2/ethylene-responsive element binding factor (AP2/ERF) transcription factor, PagERF81, from poplar 84 K (Populus alba × P. glandulosa) is highly expressed in expanding secondary xylem cells. Two independent homozygous Pagerf81 mutant lines created by gene editing, produced significantly more but smaller vessel cells and longer fiber cells with more lignin in cell walls, while PagERF81 overexpression lines had less lignin, compared to non-transgenic controls. Transcriptome and reverse transcription quantitative PCR data revealed that multiple lignin biosynthesis genes including Cinnamoyl CoA reductase 1 (PagCCR1), Cinnamyl alcohol dehydrogenase 6 (PagCAD6), and 4-Coumarate-CoA ligase-like 9 (Pag4CLL9) were up-regulated in Pagerf81 mutants, but down-regulated in PagERF81 overexpression lines. In addition, a transient transactivation assay revealed that PagERF81 repressed the transcription of these three genes. Furthermore, yeast one hybrid and electrophoretic mobility shift assays showed that PagERF81 directly bound to a GCC sequence in the PagCCR1 promoter. No known vessel or fiber cell differentiation related genes were differentially expressed, so the smaller vessel cells and longer fiber cells observed in the Pagerf81 lines might be caused by abnormal lignin deposition in the secondary cell walls. This study provides insight into the regulation of lignin biosynthesis, and a molecular tool to engineer wood with high lignin content, which would contribute to the lignin-related chemical industry and carbon sequestration.

An Efficient Agrobacterium-Mediated Transformation Method for Hybrid Poplar 84K (Populus alba × P. glandulosa) Using Calli as Explants.

A highly efficient Agrobacterium-mediated transformation method is needed for the molecular study of model tree species such as hybrid poplar 84K (Populus alba × P. glandulosa cv. '84K'). In this study, we report a callus-based transformation method that exhibits high efficiency and reproducibility. The optimized callus induction medium (CIM1) induced the development of calli from leaves with high efficiency, and multiple shoots were induced from calli growing on the optimized shoot induction medium (SIM1). Factors affecting the transformation frequency of calli were optimized as follows: Agrobacterium concentration sets at an OD600 of 0.6, Agrobacterium infective suspension with an acetosyringone (AS) concentration of 100 µM, infection time of 15 min, cocultivation duration of 2 days and precultivation duration of 6 days. Using this method, transgenic plants are obtained within approximately 2 months with a transformation frequency greater than 50%. Polymerase chain reaction (PCR), reverse transcription-PCR (RT-PCR) and ß-galactosidase (GUS) histochemical staining analyses confirmed the successful generation of stable transformants. Additionally, the calli from leaves were subcultured and used to obtain new explants; the high transformation efficiency was still maintained in subcultured calli after 6 cycles. This method provides a reference for developing effective transformation protocols for other poplar species.

The role of senescence-associated gene101 (PagSAG101a) in the regulation of secondary xylem formation in poplar.

Wood is produced by the accumulation of secondary xylem via proliferation and differentiation of the cambium cells in woody plants. Identifying the regulators involved in this process remains a challenging task. In this study, we isolated PagSAG101a, the homolog of Arabidopsis thaliana SAG101, from a hybrid poplar (Populus alba × Populus glandulosa) clone 84K and investigated its role in secondary xylem development. PagSAG101a was expressed predominantly in lignified stems and localized in the nucleus. Compared with non-transgenic 84K plants, transgenic plants overexpressing PagSAG101a displayed increased plant height, internode number, stem diameter, xylem width, and secondary cell wall thickness, while opposite phenotypes were observed for PagSAG101a knock-out plants. Transcriptome analyses revealed that differentially expressed genes were enriched for those controlling cambium cell division activity and subsequent secondary cell wall deposition during xylem formation. In addition, the tandem CCCH zinc finger protein PagC3H17, which positively regulates secondary xylem width and secondary wall thickening in poplar, could bind to the promoter of PagSAG101a and mediate the regulation of xylem differentiation. Our results support that PagSAG101a, downstream of PagC3H17, functions in wood development.

PagGRF12a interacts with PagGIF1b to regulate secondary xylem development through modulating PagXND1a expression in Populus alba × P. glandulosa.

Growth-regulating factors (GRFs) are important regulators of plant development and growth, but their possible roles in xylem development in woody plants remain unclear. Here, we report that Populus alba × Papulus glandulosa PagGRF12a negatively regulates xylem development in poplar. PagGRF12a is expressed in vascular tissues. Compared to non-transgenic control plants, transgenic poplar plants overexpressing PagGRF12a exhibited reduced xylem width and plants with repressed expression of PagGRF12a exhibited increased xylem width. Xylem NAC domain 1 (XND1) encodes a NAC domain transcription factor that regulates xylem development and transcriptional analyses revealed that PagXND1a is highly upregulated in PagGRF12a-overexpressing plants and downregulated in PagGRF12a-suppressed plants, indicating that PagGRF12a may regulate xylem development through PagXND1a. Transient transcriptional assays and chromatin immunoprecipitation-polymerase chain reaction assays confirmed that PagGRF12a directly upregulates PagXND1a. In addition, PagGRF12a interacts with the GRF-Interacting Factor (GIF) PagGIF1b, and this interaction enhances the effects of PagGRF12a on PagXND1a. Our results indicate that PagGRF12a inhibits xylem development by upregulating the expression of PagXND1a.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees.

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

KNAT2/6b, a class I KNOX gene, impedes xylem differentiation by regulating NAC domain transcription factors in poplar.

Wood formation is the terminal differentiation of xylem mother cells derived from cambial initials, and negative regulators play important roles in xylem differentiation. The molecular mechanism of the negative regulator of xylem differentiation PagKNAT2/6b was investigated. PagKNAT2/6b is an ortholog of Arabidopsis KNAT2 and KNAT6 that is highly expressed in phloem and xylem. Compared to nontransgenic control plants, transgenic poplar plants overexpressing PagKNAT2/6b present with altered vascular patterns, characterized by decreased secondary xylem with thin cell walls containing less cellulose, xylose and lignin. RNA sequencing analyses revealed that differentially expressed genes are enriched in xylem differentiation and secondary wall synthesis functions. Expression of NAM/ATAF/CUC (NAC) domain genes including PagSND1-A1, PagSND1-A2, PagSND1-B2 and PagVND6-C1 is downregulated by PagKNAT2/6b, while PagXND1a is directly upregulated. Accordingly, the dominant repression form of PagKNAT2/6b leads to increased xylem width per stem diameter through downregulation of PagXND1a. PagKNAT2/6b can inhibit cell differentiation and secondary wall deposition during wood formation in poplar by modulating the expression of NAC domain transcription factors. Direct activation of PagXND1a by PagKNAT2/6b is a key node in the negative regulatory network of xylem differentiation by KNOXs.

Brassinosteroid regulation of wood formation in poplar.

Brassinosteroids have been implicated in the differentiation of vascular cell types in herbaceous plants, but their roles during secondary growth and wood formation are not well defined. Here we pharmacologically and genetically manipulated brassinosteroid levels in poplar trees and assayed the effects on secondary growth and wood formation, and on gene expression within stems. Elevated brassinosteroid levels resulted in increases in secondary growth and tension wood formation, while inhibition of brassinosteroid synthesis resulted in decreased growth and secondary vascular differentiation. Analysis of gene expression showed that brassinosteroid action is positively associated with genes involved in cell differentiation and cell-wall biosynthesis. The results presented here show that brassinosteroids play a foundational role in the regulation of secondary growth and wood formation, in part through the regulation of cell differentiation and secondary cell wall biosynthesis.

WUSCHEL-related homeobox gene PagWOX11/12a responds to drought stress by enhancing root elongation and biomass growth in poplar.

In plants, a large root system improves the uptake of water and nutrients, and is important for responding to drought stress. The poplar WUSCHEL-related homeobox (WOX) transcription factor promotes adventitious rooting, but its regulation of root growth in response to drought stress remains elusive. In this study, we found that PagWOX11/12a from hybrid poplar 84K (Populus alba×Populus glandulosa) is expressed predominantly in the roots and is strongly induced by drought stress. Compared with non-transgenic 84K plants, transgenic poplar plants overexpressing PagWOX11/12a displayed increased root biomass and enhanced drought tolerance, while opposite phenotypes were observed for PagWOX11/12a dominant repression plants. PagWOX11/12a functions as a nuclear transcriptional activator with a transactivation domain at the C-terminus. In addition, PagERF35 was found to specifically bind to a dehydration-responsive element (DRE) within the PagWOX11/12a promoter and activate PagWOX11/12a gene expression. These results indicate that PagERF35 may activate PagWOX11/12a expression in response to drought stress by promoting root elongation and biomass, thereby increasing drought tolerance of poplar.

Genome resequencing reveals demographic history and genetic architecture of seed salinity tolerance in Populus euphratica.

Populus euphratica is a dominant tree species in desert riparian forests and possesses extraordinary adaptation to salinity stress. Exploration of its genomic variation and molecular underpinning of salinity tolerance is important for elucidating population evolution and identifying stress-related genes. Here, we identify approximately 3.15 million single nucleotide polymorphisms using whole-genome resequencing. The natural populations of P. euphratica in northwest China are divided into four distinct clades that exhibit strong geographical distribution patterns. Pleistocene climatic fluctuations and tectonic deformation jointly shaped the extant genetic patterns. A seed germination rate-based salinity tolerance index was used to evaluate seed salinity tolerance of P. euphratica and a genome-wide association study was implemented. A total of 38 single nucleotide polymorphisms were associated with seed salinity tolerance and were located within or near 82 genes. Expression profiles showed that most of these genes were regulated under salt stress, revealing the genetic complexity of seed salinity tolerance. Furthermore, DEAD-box ATP-dependent RNA helicase 57 and one undescribed gene (CCG029559) were demonstrated to improve the seed salinity tolerance in transgenic Arabidopsis. These results provide new insights into the demographic history and genetic architecture of seed salinity tolerance in desert poplar.

An essential role for Arabidopsis Trs33 in cell growth and organization in plant apical meristems.

KEY MESSAGE: Trafficking protein particle (TRAPP) complexes subunit gene AtTrs33 plays an important role in keeping apical meristematic activity and dominance in Arabidopsis. TRAPP complexes, composed of multimeric subunits, are guanine-nucleotide exchange factors for certain Rab GTPases and are believed to be involved in the regulation of membrane trafficking, but the cases in Arabidopsis are largely unknown. Trs33, recently proposed to be a component of TRAPP IV, is non-essential in yeast cells. A single copy of Trs33 gene, AtTrs33, was identified in Arabidopsis. GUS activity assay indicated that AtTrs33 was ubiquitously expressed. Based on a T-DNA insertion line, we found that loss-of-function of AtTrs33 is lethal for apical growth. Knock-down or knock-in of AtTrs33 affects apical meristematic growth and fertility, which indicates that AtTrs33 plays an important role in keeping apical meristematic activity and dominance in Arabidopsis. Analysis of auxin responses and PIN1/2 localization indicate that impaired apical meristematic activity and dominance were caused by altered auxin responses through non-polarized PIN1 localization. The present study reported that AtTrs33 plays an essential role in Arabidopsis cell growth and organization, which is different with its homologue in yeast. These findings provide new insights into the functional divergence of TRAPP subunits.

The auxin receptor TIR1 homolog (PagFBL 1) regulates adventitious rooting through interactions with Aux/IAA28 in Populus.

Adventitious roots occur naturally in many species and can also be induced from explants of some tree species including Populus, providing an important means of clonal propagation. Auxin has been identified as playing a crucial role in adventitious root formation, but the associated molecular regulatory mechanisms need to be elucidated. In this study, we examined the role of PagFBL1, the hybrid poplar (Populus alba × P. glandulosa clone 84K) homolog of Arabidopsis auxin receptor TIR1, in adventitious root formation in poplar. Similar to the distribution pattern of auxin during initiation of adventitious roots, PagFBL1 expression was concentrated in the cambium and secondary phloem in stems during adventitious root induction and initiation phases, but decreased in emerging adventitious root primordia. Overexpressing PagFBL1 stimulated adventitious root formation and increased root biomass, while knock-down of PagFBL1 transcript levels delayed adventitious root formation and decreased root biomass. Transcriptome analyses of PagFBL1 overexpressing lines indicated that an extensive remodelling of gene expression was stimulated by auxin signalling pathway during early adventitious root formation. In addition, PagIAA28 was identified as downstream targets of PagFBL1. We propose that the PagFBL1-PagIAA28 module promotes adventitious rooting and could be targeted to improve Populus propagation by cuttings.

Retraction Note to: An assessment of transgenomics as a tool for gene discovery in Populus euphratica Oliv.

This article (Zhou et al. 2018) has been retracted by the authors because the sequence BIBAC 002A111F06 was incorrectly assigned to the wrong bacterial species. The BIBAC 002A111F06 sequence (GenBank Accession KC129717) reported in the paper was attributed to Populus euphratica Oliv. The BLAST search of this KC129717 sequence against the nr database at NCBI showed that it has very high similarity to a genomic sequence from the gram-negative bacteria Stenotrophomonas maltophilia. The bacterium associates with Populus euphratica Oliv. and DNA isolated from Populus euphratica Oliv. for the construction of the BIBAC clone library inlcuded DNA from Stenotrophomonas maltophilia. Therefore, the phenotype of the transgenic Arabidopsis line carrying the KC129717 sequence cannot be attributed to genes from Populus euphratica Oliv. The authors apologize for the confusion and misinterpretation of our data resulting from the incorrect sequence assignment. All authors agree to this retraction.

An assessment of transgenomics as a tool for gene discovery in Populus euphratica Oliv.

KEY MESSAGE: Transgenomics for gene discovery in Populus euphratica. Transgenomics, a member of the omics family of methodologies, is characterized as the introduction of DNA from one organism into another on a genome-wide scale followed by the identification of recipients with altered phenotypes. This strategy allows investigators to identify the gene(s) involved in these phenotypic changes. It is particularly promising for woody plants that have a long life cycle and for which molecular tools are limited. In this study, we constructed a large-insert binary bacterial artificial chromosome library of Populus euphratica, a stress-tolerant poplar species, which included 55,296 clones with average insert sizes of about 127 kb. To date, 1077 of the clones have been transformed into Arabidopsis thaliana via Agrobacterium by the floral dip method. Of these, 69 transgenic lines showed phenotypic changes represented by diverse aspects of plant form and development, 22 of which were reproducibly associated with the same phenotypic change. One of the clones conferring transgenic plants with increased salt tolerance, 002A1F06, was further analyzed and the 127,284 bp insert in this clone harbored eight genes that have been previously reported to be involved in stress resistance. This study demonstrates that transgenomics is useful in the study of functional genomics of woody plants and in the identification of novel gene(s) responsible for economically important traits. Thus, transgenomics can also be used for validation of quantitative trait loci mapped by molecular markers.

Transcriptional and temporal response of Populus stems to gravi-stimulation.

Plants modify development in response to external stimuli, to produce new growth that is appropriate for environmental conditions. For example, gravi-stimulation of leaning branches in angiosperm trees results in modifications of wood development, to produce tension wood that pulls leaning stems upright. Here, we use gravi-stimulation and tension wood response to dissect the temporal changes in gene expression underlying wood formation in Populus stems. Using time-series analysis of seven time points over a 14-d experiment, we identified 8,919 genes that were differentially expressed between tension wood (upper) and opposite wood (lower) sides of leaning stems. Clustering of differentially expressed genes showed four major transcriptional responses, including gene clusters whose transcript levels were associated with two types of tissue-specific impulse responses that peaked at about 24-48 h, and gene clusters with sustained changes in transcript levels that persisted until the end of the 14-d experiment. Functional enrichment analysis of those clusters suggests they reflect temporal changes in pathways associated with hormone regulation, protein localization, cell wall biosynthesis and epigenetic processes. Time-series analysis of gene expression is an underutilized approach for dissecting complex developmental responses in plants, and can reveal gene clusters and mechanisms influencing development.

PagEXPA1 combines with PagCDKB2;1 to regulate plant growth and the elongation of fibers in Populus alba × Populus glandulosa.

Expansins are important plant cell wall proteins. They can loosen and soften the cell walls and lead to wall extension and cell expansion. To investigate their role in wood formation and fiber elongation, the PagEXPA1 that highly expressed in cell differentiation and expansion tissues was cloned from 84K poplar (Populus alba × P. glandulosa). The subcellular localization showed that PagEXPA1 located in the cell wall and it was highly expressed in primary stems and young leaves. Compared with non-transgenic 84K poplar, overexpression of PagEXPA1 can promote plant-growth, lignification, and fiber cell elongation, while PagEXPA1 Cas9-editing mutant lines exhibited the opposite phenotype. Transcriptome analysis revealed that DEGs were mainly enriched in some important processes, which are associated with cell wall formation and cellulose synthesis. The protein interaction prediction and expression analysis showed that PagCDKB2:1 and PagEXPA1 might have an interaction relationship. The luciferase complementary assay and bimolecular fluorescence complementary assay validated that PagEXPA1 can combined with PagCDKB2;1. So they promoted the expansion of xylem vascular tissues and the development of poplar though participating in the regulation of cell division and differentiation by programming the cell-cycle. It provides good foundation for molecular breeding of fast-growing and high-quality poplar varieties.

PagKNAT2/6b regulates tension wood formation and gravitropism by targeting cytokinin metabolism.

Tension wood is a specialized xylem tissue associated with gravitropism in angiosperm trees. However, few regulators of tension wood formation have been identified. The molecular mechanisms underpinning tension wood formation remain elusive. Here, we report that a Populus KNOTTED-like homeobox gene, PagKNAT2/6b, is involved in tension wood formation and gravity response. Transgenic poplar plants overexpressing PagKNAT2/6b displayed more sensitive gravitropism than controls, as indicated by increased stem curvature. Microscopic examination revealed greater abundance of fibre cells with a gelatinous cell wall layer (G-layer) and asymmetric growth of secondary xylem in PagKNAT2/6b overexpression lines. Conversely, PagKNAT2/6b dominant repression plants exhibited decreased tension wood formation and reduced response to gravity stimulation. Moreover, sensitivity to gravity stimulation showed a negative relationship with development stage. Expression of genes related to growth and senescence was affected in PagKNAT2/6b transgenic plants. More importantly, transcription activation and electrophoretic mobility shift assays suggested that PagKNAT2/6b promotes the expression of cytokinin metabolism genes. Consistently, cytokinin content was increased in PagKNAT2/6b overexpression plants. Therefore, PagKNAT2/6b is involved in gravitropism and tension wood formation, likely via modulation of cytokinin metabolism.

A method for analyzing programmed cell death in xylem development by flow cytometry.

Programmed cell death (PCD) is a genetically regulated developmental process leading to the death of specific types of plant cells, which plays important roles in plant development and growth such as wood formation. However, an efficient method needs to be established to study PCD in woody plants. Flow cytometry is widely utilized to evaluate apoptosis in mammalian cells, while it is rarely used to detect PCD in plants, especially in woody plants. Here, we reported that the xylem cell protoplasts from poplar stem were stained with a combination of fluorescein annexin V-FITC and propidium iodide (PI) and then sorted by flow cytometry. As expected, living cells (annexin V-FITC negative/PI negative), early PCD cells (annexin V-FITC positive/PI negative), and late PCD cells (annexin V-FITC positive/PI positive) could be finely distinguished through this method and then subjected for quantitative analysis. The expression of cell-type- and developmental stages-specific marker genes was consistent with the cell morphological observation. Therefore, the newly developed fluorescence-activated cell sorting (FACS) method can be used to study PCD in woody plants, which will be beneficial for studying the molecular mechanisms of wood formation.

Bacillus thuringiensis Cry3Aa fused to a cellulase-binding peptide shows increased toxicity against the longhorned beetle.

Cry3 class toxins are used extensively for biological control of coleopteran larvae. We previously identified a peptide (PCx) from a phage display library that specifically binds Cx-cellulase from the midgut of Anoplophora glabripennis Motschulsky (Asian longhorn beetle) larvae. Here, we added a DNA fragment that encodes the peptide onto either end of the cry3Aa gene and tested the expressed PCx-Cry3Aa and Cry3Aa-PCx proteins for insecticidal activity in the longhorned beetle. An insect bioassay revealed that, compared with native Cry3Aa, the two modified Cry3Aa proteins had significantly higher lethality, with PCx-Cry3Aa exhibiting a mortality rate almost three times that of Cry3Aa. We also proposed that the increased lethality in larvae fed with PCx-Cry3Aa or Cry3Aa-PCx would be attributable to the binding of the toxin with Cx-cellulase, thereby increasing toxin retention in the midgut. The significantly enhanced insecticidal activity of Cry3Aa fused with the Cx-cellulase binding peptide provides a new strategy for increasing toxin efficacy against the longhorned beetle. These uniquely modified Cry3Aa proteins have potential use for pest control.

Integrated transcriptome and proteome analysis reveals brassinosteroid-mediated regulation of cambium initiation and patterning in woody stem.

Wood formation involves sequential developmental events requiring the coordination of multiple hormones. Brassinosteroids (BRs) play a key role in wood development, but little is known about the cellular and molecular processes that underlie wood formation in tree species. Here, we generated transgenic poplar lines with edited PdBRI1 genes, which are orthologs of Arabidopsis vascular-enriched BR receptors, and showed how inhibition of BR signaling influences wood development at the mRNA and/or proteome level. Six Populus PdBRI1 genes formed three gene pairs, each of which was highly expressed in basal stems. Simultaneous mutation of PdBRI1-1, -2, -3 and - 6, which are orthologs of the Arabidopsis vascular-enriched BR receptors BRI1, BRL1 and BRL3, resulted in severe growth defects. In particular, the stems of these mutant lines displayed a discontinuous cambial ring and patterning defects in derived secondary vascular tissues. Abnormal cambial formation within the cortical parenchyma was also observed in the stems of pdbri1-1;2;3;6. Transgenic poplar plants expressing edited versions of PdBRI1-1 or PdBRI1-1;2;6 exhibited phenotypic alterations in stem development at 4.5 months of growth, indicating that there is functional redundancy among these PdBRI1 genes. Integrated analysis of the transcriptome and proteome of pdbri1-1;2;3;6 stems revealed differential expression of a number of genes/proteins associated with wood development and hormones. Concordant (16%) and discordant (84%) regulation of mRNA and protein expression, including wood-associated mRNA/protein expression, was found in pdbri1-1;2;3;6 stems. This study found a dual role of BRs in procambial cell division and xylem differentiation and provides insights into the multiple layers of gene regulation that contribute to wood formation in Populus.