Molecular control of wood formation in trees

Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem

Silver birch (Betula pendula Roth) is an economically important species in Northern Europe. The current research focused on the molecular background of different xylogenesis scenarios in the birch trunks. The study objects were two forms of silver birch, silver birch trees, and Karelian birch trees; the latter form is characterized by the formation of two types of wood, non-figured (straight-grained) and figured, respectively, while it is currently not clear which factors cause this difference. We identified VND/NST/SND genes that regulate secondary cell wall biosynthesis in the birch genome and revealed differences in their expression in association with the formation of xylem with different ratios of structural elements. High expression levels of BpVND7 accompanied differentiation of the type of xylem which is characteristic of the species. At the same time, the appearance of figured wood was accompanied by the low expression levels of the VND genes and increased levels of expression of NST and SND genes. We identified BpARF5 as a crucial regulator of auxin-dependent vascular patterning and its direct target—BpHB8. A decrease in the BpARF5 level expression in differentiating xylem was a specific characteristic of both Karelian birch with figured and non-figured wood. Decreased BpARF5 level expression in non-figured trees accompanied by decreased BpHB8 and VND/NST/SND expression levels compared to figured Karelian birch trees. According to the results obtained, we suggested silver birch forms differing in wood anatomy as valuable objects in studying the regulation of xylogenesis.

Investigation Into Different Wood Formation Mechanisms Between Angiosperm and Gymnosperm Tree Species at the Transcriptional and Post-transcriptional Level

Enormous distinctions of the stem structure and cell types between gymnosperms and angiosperms tree species are expected to cause quite different wood physical and mechanical attributes, however, the molecular mechanisms underlying the differing wood morphology are still unclear. In this study, we compared the transcriptomes obtained by RNA-Seq between Populus alba × P. glandulosa clone 84K, and Larix kaempferi (Lamb.) Carr trees. Available genome resource served as reference for P. alba × P. glandulosa and the Iso-Seq results of a three-tissues mixture (xylem, phloem, and leaf) were used as the reference for L. kaempferi to compare the xylem-specifically expressed genes and their alternative splicing model. Through screening, we obtained 13,907 xylem-specifically expressed genes (5,954 up-regulated, 7,953 down-regulated) in the xylem of P. alba × P. glandulosa, and 2,596 xylem-specifically expressed genes (1,648 up-regulated, 948 down-regulated) in the xylem of L. kaempferi. From the GO and KEGG analyses, some genes associated with two wood formation-related pathways, namely those for phenylpropanoid biosynthesis, and starch and sucrose metabolism, were successfully screened. Then the distributions and gene expression models between P. alba × P. glandulosa and L. kaempferi in those pathways were compared, which suggested differential wood formation processes between the angiosperm and gymnosperm trees. Furthermore, a Weight Gene Co-expression Network Analysis (WGCNA) for total xylem-specifically expressed genes in two species was conducted, from which wood formation-related modules were selected to build a co-expression network for the two tree species. The genes within this co-expression network showed different co-expression relationships between the angiosperm and gymnosperm woody species. Comparing the alternative splicing events for wood formation-related genes suggests a different post-transcriptional regulation process exists between the angiosperm and gymnosperm trees. Our research thus provides the foundation for the in-depth investigation of different wood formation mechanisms of angiosperm and gymnosperm species.

BEL1-like Homeodomain Protein BLH6a Is a Negative Regulator of CAl5H2 in Sinapyl Alcohol Monolignol Biosynthesis in Poplar

Lignin is one of the major components of xylem cell walls in tree stems. The lignin in the wood of most flowering plants (dicotyledonous angiosperms) is typically polymerized from three monolignol precursors, coniferyl alcohol, sinapyl alcohol, and p-coumaroyl alcohol, resulting in guaiacyl (G), syringyl (S), and hydroxyphenyl (H) subunits, respectively. In this study, we focus on the transcriptional regulation of a coniferaldehyde 5-hydroxylase (CAld5H2) gene, which encodes a key enzyme for sinapyl alcohol biosynthesis. We carried out a yeast one-hybrid (Y1H) screen to identify candidate upstream transcription factors (TFs) regulating CAld5H2. We obtained 12 upstream TFs as potential regulators of CAld5H2. One of these TF genes, BLH6a, encodes a BEL1-like homeodomain (BLH) protein and negatively regulated the CAld5H2 promoter activity. The direct regulation of CAld5H2 promoter by BLH6a was supported by chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) and dominant repression of BLH6a in transgenic plants. Luciferase complementation imaging analyses showed extensive protein-protein interactions among these 12 TFs. We propose that BLH6a is a negative regulator of CAld5H2, which acts through combinatorial regulation of multiple TFs for sinapyl alcohol (S monolignol) biosynthesis in poplar.

Ubiquitylome analysis reveals the involvement of ubiquitination in the bast fiber growth of ramie

A total of 249 sites from 197 proteins showed a differential ubiquitination level in the fiber development of ramie barks. The function of two differentially ubiquitinated proteins for fiber growth was demonstrated. Ubiquitination is one of the most common post-translational modifications of proteins, and it plays essential roles in plant growth and development. However, the involvement of ubiquitination in the growth of plant fibers remains largely unknown. We compared the ubiquitylome of the top and middle stems of ramie bark, with different fiber growth stages. We identified 249 differentially ubiquitinated sites in 197 proteins in fiber-developing barks in the stems and found that seven were homologs of Arabidopsis proteins associated with fiber growth. Overexpression of the differentially ubiquitinated proteins, RWA3 homolog whole_GLEAN_10024150 and MYB protein whole_GLEAN_10015497, significantly promoted fiber growth in transgenic Arabidopsis, indicating their involvement in this process. We also found that the abundance of these proteins decreased when their ubiquitination levels increased and vice versa in the fiber-developing bark. These results indicated that the abundance of these two proteins was adjusted through ubiquitin-dependent degradation. Collectively, our findings provide important insights into the involvement of ubiquitination in the growth of ramie fibers.

Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa

Poplar wood is the main source of renewable biomass energy worldwide, and is also considered to be a model system for studying woody plants. The Full-length cDNA Over-eXpressing (FOX) gene hunting system is an effective method for generating gain-of-function mutants. Large numbers of novel genes have successfully been identified from many herbaceous plants according to the phenotype of gain-of-function mutants under normal or abiotic stress conditions using this system. However, the system has not been used for functional gene identification with high-throughput mutant screening in woody plants. In this study, we constructed a FOX library from the Chinese white poplar, Populus tomentosa. The poplar cDNA library was constructed into the plant expression vector pEarleyGate101 and further transformed into Arabidopsis thaliana (thale cress). We collected 1749 T1 transgenic plants identified by PCR. Of these, 593 single PCR bands from different transgenic lines were randomly selected for sequencing, and 402 diverse sequences of poplar genes were isolated. Most of these genes were involved in photosynthesis, environmental adaptation, and ribosome biogenesis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation. We characterized in detail two mutant lines carrying PtoCPCa or PtoWRKY13 cDNA insertions. Phenotypic characterization showed that overexpression of these genes in A. thaliana affected trichome development or secondary cell wall (SCW) deposition, respectively. Together, the Populus-FOX-Arabidopsis library generated in our experiments will be helpful for efficient discovery of novel genes in poplar.

Synergies and Entanglement in Secondary Cell Wall Development and Abiotic Stress Response in Trees

A major challenge for sustainable food, fuel, and fiber production is simultaneous genetic improvement of yield, biomass quality, and resilience to episodic environmental stress and climate change. For Populus and other forest trees, quality traits involve alterations in the secondary cell wall (SCW) of wood for traditional uses, as well as for a growing diversity of biofuels and bioproducts. Alterations in wood properties that are desirable for specific end uses can have negative effects on growth and stress tolerance. Understanding of the diverse roles of SCW genes is necessary for the genetic improvement of fast-growing, short-rotation trees that face perennial challenges in their growth and development. Here, we review recent progress into the synergies and antagonisms of SCW development and abiotic stress responses, particularly, the roles of transcription factors, SCW biogenesis genes, and paralog evolution.

Elicitors Modulate Young Sandalwood (Santalum album L.) Growth, Heartwood Formation, and Concrete Oil Synthesis

Five chemical elicitors––6-benzyladenine (BA), ethephon (ETH), methyl jasmonate (MeJA), hydrogen peroxide (H₂O₂) and calcium chloride (CaCl₂)––were used to treat 1- and 5-year-old sandal trees (Santalum album L.) to assess their effects on growth, heartwood formation and concrete oil synthesis. The results showed that some newly formed branches in stems that were induced by BA and ETH displayed leaf senescence and developed new smaller and light-green leaves. The relative percentage of concrete oil from the heartwood of water-treated trees (0.65%) was significantly lower than that from trees treated with 4 mM H₂O₂ (2.85%) and 4 mM BA (2.75%) within one year. Four mM BA, H₂O₂ and CaCl₂ induced a significantly higher level of sesquiterpenoids than heartwood treated with 2 mM of these elicitors. Four mM MeJA induced significantly less sesquiterpenoids than heartwood treated with 2 mM MeJA. Morphological, physiological, and chromatographic–spectrometric technologies were integrated to trace the potential function of these exogenously applied chemical elicitors. The results may have important applications and provide a better understanding of the molecular mechanism of heartwood formation and hardening in young sandalwood trees.

DELLA-NAC Interactions Mediate GA Signaling to Promote Secondary Cell Wall Formation in Cotton Stem

Gibberellins (GAs) promote secondary cell wall (SCW) development in plants, but the underlying molecular mechanism is still to be elucidated. Here, we employed a new system, the first internode of cotton, and the virus-induced gene silencing method to address this problem. We found that knocking down major DELLA genes via VIGS phenocopied GA treatment and significantly enhanced SCW formation in the xylem and phloem of cotton stems. Cotton DELLA proteins were found to interact with a wide range of SCW-related NAC proteins, and virus-induced gene silencing of these NAC genes inhibited SCW development with downregulated biosynthesis and deposition of lignin. The findings indicated a framework for the GA regulation of SCW formation; that is, the interactions between DELLA and NAC proteins mediated GA signaling to regulate SCW formation in cotton stems.

Characterization of a High Hierarchical Regulator, PtrGATA12, Functioning in Differentially Regulating Secondary Wall Component Biosynthesis in Populus trichocarpa

In plants, GATA transcription factors (TFs) have been reported to play vital roles in to a wide range of biological processes. To date, there is still no report about the involvement and functions of woody plant GATA TFs in wood formation. In this study, we described the functional characterization of a Populus trichocarpa GATA TF, PtrGATA12, which encodes a nuclear-localized transcriptional activator predominantly expressing in developing xylem tissues. Overexpression of PtrGATA12 not only inhibited growths of most phenotypic traits and biomass accumulation, but also altered the expressions of some master TFs and pathway genes involved in secondary cell wall (SCW) and programmed cell death, leading to alternated SCW components and breaking forces of stems of transgenic lines. The significant changes occurred in the contents of hemicellulose and lignin and SCW thicknesses of fiber and vessel that increased by 13.5 and 10.8%, and 20.83 and 11.83%, respectively. Furthermore, PtrGATA12 bound directly to the promoters of a battery of TFs and pathway genes and activated them; the binding sites include two cis-acting elements that were specifically enriched in their promoter regions. Taken together, our results suggest PtrGATA12, as a higher hierarchical TF on the top of PtrWND6A, PtrWND6B, PtrMYB152, and PtrMYB21, exert a coordinated regulation of SCW components biosynthesis pathways through directly and indirectly controlling master TFs, middle-level TFs, and further downstream pathway genes of the currently known hierarchical transcription network that governs SCW formation.

Identifying transcription factors that reduce wood recalcitrance and improve enzymatic degradation of xylem cell wall in Populus

Developing an efficient deconstruction step of woody biomass for biorefinery has been drawing considerable attention since its xylem cell walls display highly recalcitrance nature. Here, we explored transcriptional factors (TFs) that reduce wood recalcitrance and improve saccharification efficiency in Populus species. First, 33 TF genes up-regulated during poplar wood formation were selected as potential regulators of xylem cell wall structure. The transgenic hybrid aspens (Populus tremula × Populus tremuloides) overexpressing each selected TF gene were screened for in vitro enzymatic saccharification. Of these, four transgenic seedlings overexpressing previously uncharacterized TF genes increased total glucan hydrolysis on average compared to control. The best performing lines overexpressing Pt × tERF123 and Pt × tZHD14 were further grown to form mature xylem in the greenhouse. Notably, the xylem cell walls exhibited significantly increased total xylan hydrolysis as well as initial hydrolysis rates of glucan. The increased saccharification of Pt × tERF123-overexpressing lines could reflect the improved balance of cell wall components, i.e., high cellulose and low xylan and lignin content, which could be caused by upregulation of cellulose synthase genes upon the expression of Pt × tERF123. Overall, we successfully identified Pt × tERF123 and Pt × tZHD14 as effective targets for reducing cell wall recalcitrance and improving the enzymatic degradation of woody plant biomass.

A metabolite roadmap of the wood-forming tissue in Populus tremula

Wood, or secondary xylem, is the product of xylogenesis, a developmental process that begins with the proliferation of cambial derivatives and ends with mature xylem fibers and vessels with lignified secondary cell walls. Fully mature xylem has undergone a series of cellular processes, including cell division, cell expansion, secondary wall formation, lignification and programmed cell death. A complex network of interactions between transcriptional regulators and signal transduction pathways controls wood formation. However, the role of metabolites during this developmental process has not been comprehensively characterized. To evaluate the role of metabolites during wood formation, we performed a high spatial resolution metabolomics study of the wood-forming zone of Populus tremula, including laser dissected aspen ray and fiber cells. We show that metabolites show specific patterns within the wood-forming zone, following the differentiation process from cell division to cell death. The data from profiled laser dissected aspen ray and fiber cells suggests that these two cell types host distinctly different metabolic processes. Furthermore, by integrating previously published transcriptomic and proteomic profiles generated from the same trees, we provide an integrative picture of molecular processes, for example, deamination of phenylalanine during lignification is of critical importance for nitrogen metabolism during wood formation.

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.

DNA methylation and its effects on gene expression during primary to secondary growth in poplar stems

Background

As an important epigenetic mark, 5-methylcytosine (5mC) methylation is involved in many DNA-dependent biological processes and plays a role during development and differentiation of multicellular organisms. However, there is still a lack of knowledge about the dynamic aspects and the roles of global 5mC methylation in wood formation in tree trunks. In this study, we not only scrutinized single-base resolution methylomes of primary stems (PS), transitional stems (TS), and secondary stems (SS) of Populus trichocarpa using a high-throughput bisulfite sequencing technique, but also analyzed the effects of 5mC methylation on the expression of genes involved in wood formation.

Results

The overall average percentages of CG, CHG, and CHH methylation in poplar stems were ~ 53.6%, ~ 37.7%, and ~ 8.5%, respectively, and the differences of 5mC in genome-wide CG/CHG/CHH contexts among PS, TS, and SS were statistically significant (p < 0.05). The evident differences in CG, CHG, and CHH methylation contexts among 2 kb proximal promoters, gene bodies, and 2 kb downstream regions were observed among PS, TS, and SS. Further analysis revealed a perceptible global correlation between 5mC methylation levels of gene bodies and transcript levels but failed to reveal a correlation between 5mC methylation levels of proximal promoter regions and transcript levels. We identified 653 and 858 DMGs and 4978 and 4780 DEGs in PS vs TS and TS vs SS comparisons, respectively. Only 113 genes of 653 DMGs and 4978 DEGs, and 114 genes of 858 DMGs and 4780 DEG were common. Counterparts of some of these common genes in other species, including Arabidopsis thaliana, are known to be involved in secondary cell wall biosynthesis and hormone signaling. This indicates that methylation may directly modulate wood formation genes and indirectly attune hormone signaling genes, which in turn impact wood formation.

Conclusions

DNA methylation only marginally affects pathway genes or regulators involved in wood formation, suggesting that further studies of wood formation should lean towards the indirect effects of methylation. The information and data we provide here will be instrumental for understanding the roles of methylation in wood formation in tree species.

Regulation of Lignin Biosynthesis by Post-translational Protein Modifications

Post-translational modification of proteins exerts essential roles in many biological processes in plants. The function of these chemical modifications has been extensively characterized in many physiological processes, but how these modifications regulate lignin biosynthesis for wood formation remained largely unknown. Over the past decade, post-translational modification of several proteins has been associated with lignification. Phosphorylation, ubiquitination, glycosylation, and S-nitrosylation of transcription factors, monolignol enzymes, and peroxidases were shown to have primordial roles in the regulation of lignin biosynthesis. The main discoveries of post-translational modifications in lignin biosynthesis are discussed in this review.

Transcriptional analysis of arogenate dehydratase genes identifies a link between phenylalanine biosynthesis and lignin biosynthesis

Biogenesis of the secondary cell wall in trees involves the massive biosynthesis of the phenylalanine-derived polymer lignin. Arogenate dehydratase (ADT) catalyzes the last, and rate-limiting, step of the main pathway for phenylalanine biosynthesis. In this study, we found that transcript levels for several members of the large ADT gene family, including ADT-A and ADT-D, were enhanced in compression wood of maritime pine, a xylem tissue enriched in lignin. Transcriptomic analysis of maritime pine silenced for PpMYB8 revealed that this gene plays a critical role in coordinating the deposition of lignin with the biosynthesis of phenylalanine. Specifically, it was found that ADT-A and ADT-D were strongly down-regulated in PpMYB8-silenced plants and that they were transcriptionally regulated through direct interaction of this transcription factor with regulatory elements present in their promoters. Another transcription factor, PpHY5, exhibited an expression profile opposite to that of PpMYB8 and also interacted with specific regulatory elements of ADT-A and ADT-D genes, suggesting that it is involved in transcriptional regulation of phenylalanine biosynthesis. Taken together, our results reveal that PpMYB8 and PpHY5 are involved in the control of phenylalanine formation and its metabolic channeling for lignin biosynthesis and deposition during wood formation in maritime pine.

Wood Architecture and Composition Are Deeply Remodeled in Frost Sensitive Eucalyptus Overexpressing CBF/DREB1 Transcription Factors

Eucalypts are the most planted trees worldwide, but most of them are frost sensitive. Overexpressing transcription factors for CRT-repeat binding factors (CBFs) in transgenic Eucalyptus confer cold resistance both in leaves and stems. While wood plays crucial roles in trees and is affected by environmental cues, its potential role in adaptation to cold stress has been neglected. Here, we addressed this question by investigating the changes occurring in wood in response to the overexpression of two CBFs, taking advantage of available transgenic Eucalyptus lines. We performed histological, biochemical, and transcriptomic analyses on xylem samples. CBF ectopic expression led to a reduction of both primary and secondary growth, and triggered changes in xylem architecture with smaller and more frequent vessels and fibers exhibiting reduced lumens. In addition, lignin content and syringyl/guaiacyl (S/G) ratio increased. Consistently, many genes of the phenylpropanoid and lignin branch pathway were upregulated. Most of the features of xylem remodeling induced by CBF overexpression are reminiscent of those observed after long exposure of Eucalyptus trees to chilling temperatures. Altogether, these results suggest that CBF plays a central role in the cross-talk between response to cold and wood formation and that the remodeling of wood is part of the adaptive strategies to face cold stress.

PtiCYP85A3, a BR C-6 Oxidase Gene, Plays a Critical Role in Brassinosteroid-Mediated Tension Wood Formation in Poplar

In angiosperm trees, the gelatinous layer (G-layer) takes a great part of the fiber cell wall in the tension wood (TW). However, the mechanism underlying G-layer formation in poplar is largely unknown. In this work, we demonstrate that G-layer formation in poplar TW cells is regulated by brassinosteroid (BR) and its signaling. PtiCYP85A3, a key BR biosynthesis gene, was predominantly expressed in the xylem of TW, accompanied with a relatively higher castasterone (CS) accumulation, than in the xylem of opposite wood (OW). A wider expression zone of BZR1, a key transcriptional factor in BR singling pathway, was also observed in G-fiber cells on TW side than in wood fiber cells on the OW side, as indicated by immunohistochemistry assays. Transgenic poplar plants overexpressing PtiCYP85A3 produced thicker G-layer with higher cellulose proportion, and accumulated more BZR1 protein in the xylem of TW than did the wild type (WT) plants. Expression of most TW-associated CesAs, which were induced by 2, 4-epibrassinolide, an active BR, and inhibited by brassinazole, a BR biosynthesis inhibitor, were also up-regulated in the xylem of TW in transgenic plants compared to that in WT plants. Further studies with dual-luciferase assays demonstrated that the promoters of PtiCesAs were activated by PtiMYB128, a TW specific transcription factor, which was then regulated by BZR1. All these results indicate that BR plays a crucial role in the G-layer formation of TW fiber cells by regulating the expression of BZR1, PtiMYB128, and PtiCesAs in poplar.

Construction of a hierarchical gene regulatory network centered around a transcription factor

We have modified a multitude of transcription factors (TFs) in numerous plant species and some animal species, and obtained transgenic lines that exhibit phenotypic alterations. Whenever we observe phenotypic changes in a TF's transgenic lines, we are always eager to identify its target genes, collaborative regulators and even upstream high hierarchical regulators. This issue can be addressed by establishing a multilayered hierarchical gene regulatory network (ML-hGRN) centered around a given TF. In this article, a practical approach for constructing an ML-hGRN centered on a TF using a combined approach of top-down and bottom-up network construction methods is described. Strategies for constructing ML-hGRNs are vitally important, as these networks provide key information to advance our understanding of how biological processes are regulated.

Transcriptomics and Proteomics Reveal the Cellulose and Pectin Metabolic Processes in the Tension Wood (Non-G-Layer) of Catalpa bungei

: Catalpa bungei is an economically important tree with high-quality wood and highly valuable to the study of wood formation. In this work, the xylem microstructure of C. bungei tension wood (TW) was observed, and we performed transcriptomics, proteomics and Raman spectroscopy of TW, opposite wood (OW) and normal wood (NW). The results showed that there was no obvious gelatinous layer (G-layer) in the TW of C. bungei and that the secondary wall deposition in the TW was reduced compared with that in the OW and NW. We found that most of the differentially expressed mRNAs and proteins were involved in carbohydrate polysaccharide synthesis. Raman spectroscopy results indicated that the cellulose and pectin content and pectin methylation in the TW were lower than those in the OW and NW, and many genes and proteins involved in the metabolic pathways of cellulose and pectin, such as galacturonosyltransferase (GAUT), polygalacturonase (PG), endoglucanase (CLE) and β-glucosidase (BGLU) genes, were significantly upregulated in TW. In addition, we found that the MYB2 transcription factor may regulate the pectin degradation genes PG1 and PG3, and ARF, ERF, SBP and MYB1 may be the key transcription factors regulating the synthesis and decomposition of cellulose. In contrast to previous studies on TW with a G-layer, our results revealed a change in metabolism in TW without a G-layer, and we inferred that the change in the pectin type, esterification and cellulose characteristics in the TW of C. bungei may contribute to high tensile stress. These results will enrich the understanding of the mechanism of TW formation.

Parallel evolution of arborescent carrots (Daucus) in Macaronesia

PREMISE

Despite intensive research, the pathways and driving forces behind the evolution of derived woodiness on oceanic islands remain obscure. The genus Daucus comprises mostly herbs (therophytes, hemicryptophytes) with few rosette treelets (chamaephytes) endemic to various Macaronesian archipelagos, suggesting their independent evolution. To elucidate the evolutionary pathways to derived woodiness, we examined phylogenetic relationships and the habit and secondary xylem evolution in Daucus and related taxa.

METHODS

Sixty taxa were surveyed for molecular markers, life history, and habit traits. Twenty-one species were considered for wood anatomical characters. A dated phylogeny was estimated using Bayesian methods. The evolution of selected traits was reconstructed using parsimony and maximum likelihood.

RESULTS

Daucus dispersed independently to the Canary Islands (and subsequently to Madeira), Cape Verde, and the Azores in the late Miocene and Pleistocene. Life span, reproductive strategy, and life form were highly homoplastic; the ancestor of Daucus was probably a monocarpic, biennial hemicryptophyte. Rosette treelets evolved independently in the Canarian-Madeiran lineage and in Cape Verde, the latter within the last 0.13 Myr. Treelets and hemicryptophytes did not differ in wood anatomy. Pervasive axial parenchyma in wood occurred more often in polycarpic rather than monocarpic species.

CONCLUSIONS

Life span and life form in Daucus are evolutionarily labile and may change independently of wood anatomy, which is related to plant reproductive strategy rather than to life form. Insular woodiness may evolve rapidly (as demonstrated in D. bischoffii), and in Daucus, it does not seem to be an adaptation to lower the risk of xylem embolism.

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.

Overexpression of PdC3H17 Confers Tolerance to Drought Stress Depending on Its CCCH Domain in Populus

Plant CCCH zinc finger proteins control growth, development, and stress responses mainly at the post-transcriptional level. Currently, limited reports are available about the roles of plant CCCH proteins in drought tolerance. In this study, we provided evidence showing that PdC3H17 from Populus deltoides × P. euramericana involves drought tolerance and response. Overexpression of PdC3H17 in poplar caused dwarf, resulted in higher stem water potential, and showed increased photosynthetic and ROS-scavenging abilities, thereby enhancing tolerance to drought stress, compared to controls. Accordingly, after drought treatment the stem elongation and thickening rates of these overexpression lines were higher than those of the controls. However, overexpression of the coding region excluding the CCCH domain of PdC3H17 roughly exhibited WT-like physiological and drought-resistant phenotypes, indicating the requirement of the CCCH domain for PdC3H17 controlling these processes. In addition, N-terminal sequence of PdC3H17 was found to possess transcriptional activity ability in yeast cells. Together, our results suggest that PdC3H17 may depend on its CCCH domain to control drought tolerance in Populus.

Identification and characterization of circular RNAs during wood formation of poplars in acclimation to low nitrogen availability

Circular RNA (circRNA) identification and expression profiles, and construction of circRNAs-miRNAs-mRNAs networks indicates that circRNAs are involved in wood formation of poplars in acclimation to low nitrogen availability. Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in wood formation of poplars in acclimation to low nitrogen (N) availability are currently unknown. Here, we undertook a systematic identification and characterization of circRNAs in the wood of Populus × canescens exposed to either 50 (low N) or 500 (normal N) µM NH₄NO₃ using rRNA-depleted RNA-sequencing. A total of 2,509 unique circRNAs were identified, and 163 (ca. 6.5%) circRNAs were significantly differentially expressed (DE) under low N condition. We observed a positive correlation between the expression patterns of DE circRNAs and their hosting protein-coding genes. Moreover, circRNAs-miRNAs-mRNAs' networks were identified in the wood of poplars under low N availability. For instance, upregulated several circRNAs, such as circRNA1226, circRNA 1732, and circRNA392 induced increases in nuclear factor Y, subunit A1-A (NFYA1-A), NFYA1-B, and NFYA10 transcript levels via the mediation of miR169b members, which is in line with reduced xylem width and cell layers of the xylem in the wood of low N-supplied poplars. Upregulation of circRNA1006, circRNA1344, circRNA1941, circRNA901, and circRNA146 caused increased transcript level of MYB61 via the mediation of a miR5021 member, corresponding well to the higher lignin concentration in the wood of low N-treated poplars. Overall, these results indicated that DE circRNAs play an essential role in regulating gene expression via circRNAs-miRNAs-mRNAs' networks to modulate wood anatomical and chemical properties of poplars in acclimation to low N availability.

Over-expression of poplar NAC15 gene enhances wood formation in transgenic tobacco

BACKGROUND

NAC (NAM/ATAF/CUC) is one of the largest plant-specific transcription factor (TF) families known to play significant roles in wood formation. Acting as master gene regulators, a few NAC genes can activate secondary wall biosynthesis during wood formation in woody plants.

RESULTS

In the present study, firstly, we screened 110 differentially expressed NAC genes in the leaves, stems, and roots of di-haploid Populus simonii×P. nigra by RNA-Seq. Then we identified a nucleus-targeted gene, NAC15 gene, which was one of the highly expressed genes in the stem among 110 NAC family members. Thirdly, we conducted expression pattern analysis of NAC15 gene, and observed NAC15 gene was most highly expressed in the xylem by RT-qPCR. Moreover, we transferred NAC15 gene into tobacco and obtained 12 transgenic lines overexpressing NAC15 gene (TLs). And the relative higher content of hemicellulose, cellulose and lignin was observed in the TLs compared to the control lines containing empty vector (CLs). It also showed darker staining in the culms of the TLs with phloroglucinol staining, compared to the CLs. Furthermore, the relative expression level of a few lignin- and cellulose-related genes was significantly higher in the TLs than that in the CLs.

CONCLUSIONS

The overall results indicated that NAC15 gene is highly expressed in the xylem of poplar and may be a potential candidate gene playing an important role in wood formation in transgenic tobacco.

PdWND3A, a wood-associated NAC domain-containing protein, affects lignin biosynthesis and composition in Populus

BACKGROUND

Plant secondary cell wall is a renewable feedstock for biofuels and biomaterials production. Arabidopsis VASCULAR-RELATED NAC DOMAIN (VND) has been demonstrated to be a key transcription factor regulating secondary cell wall biosynthesis. However, less is known about its role in the woody species.

RESULTS

Here we report the functional characterization of Populus deltoides WOOD-ASSOCIATED NAC DOMAIN protein 3 (PdWND3A), a sequence homolog of Arabidopsis VND4 and VND5 that are members of transcription factor networks regulating secondary cell wall biosynthesis. PdWND3A was expressed at higher level in the xylem than in other tissues. The stem tissues of transgenic P. deltoides overexpressing PdWND3A (OXPdWND3A) contained more vessel cells than that of wild-type plants. Furthermore, lignin content and lignin monomer syringyl and guaiacyl (S/G) ratio were higher in OXPdWND3A transgenic plants than in wild-type plants. Consistent with these observations, the expression of FERULATE 5-HYDROXYLASE1 (F5H1), encoding an enzyme involved in the biosynthesis of sinapyl alcohol (S unit monolignol), was elevated in OXPdWND3A transgenic plants. Saccharification analysis indicated that the rate of sugar release was reduced in the transgenic plants. In addition, OXPdWND3A transgenic plants produced lower amounts of biomass than wild-type plants.

CONCLUSIONS

PdWND3A affects lignin biosynthesis and composition and negatively impacts sugar release and biomass production.

Competing Endogenous RNA Networks Underlying Anatomical and Physiological Characteristics of Poplar Wood in Acclimation to Low Nitrogen Availability

Although poplar plantations are often established on nitrogen (N)-poor soil, the physiological and molecular mechanisms underlying wood properties of poplars in acclimation to low N availability remain largely unknown. To investigate wood properties of poplars in acclimation to low N, Populus � canescens saplings were exposed to either 50 (low N) or 500 (normal N) �M NH4NO3 for 2 months. Low N resulted in decreased xylem width and cell layers of the xylem (the number of cells counted along the ray parenchyma on the stem cross section), narrower lumina of vessels and fibers, greater thickness of double fiber walls (the walls between two adjacent fiber cells), more hemicellulose and lignin deposition, and reduced cellulose accumulation in poplar wood. Consistently, concentrations of gibberellins involved in cell size determination and the abundance of various metabolites including amino acids, carbohydrates and precursors for cell wall biosynthesis were decreased in low N-supplied wood. In line with these anatomical and physiological changes, a number of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) were significantly differentially expressed. Competing endogenous RNA regulatory networks were identified in the wood of low N-treated poplars. Overall, these results indicate that miRNAs-lncRNAs-mRNAs networks are involved in regulating wood properties and physiological processes of poplars in acclimation to low N availability.

Integration of Quantitative Trait Loci Mapping and Expression Profiling Analysis to Identify Genes Potentially Involved in Ramie Fiber Lignin Biosynthesis

Ramie fibers, one of the most important natural fibers in China, are mainly composed of lignin, cellulose, and hemicellulose. As the high lignin content in the fibers results in a prickly texture, the lignin content is deemed to be an important trait of the fiber quality. In this study, the genetic basis of the fiber lignin content was evaluated, resulting in the identification of five quantitative trait loci (QTLs). Three genes, whole_GLEAN_10021050, whole_GLEAN_10026962, and whole_GLEAN_10009464 that were identified on the QTL regions of qLC7, qLC10, and qLC13, respectively, were found to be homologs of the Arabidopsis lignin biosynthetic genes. Moreover, all three genes displayed differential expression in the barks located in the top and middle parts of the stem, where lignin was not being synthesized and where it was being biosynthesized, respectively. Sequence comparison found that these three genes had wide variations in their coding sequences (CDSs) and putative promoter regions between the two parents, especially the MYB gene whole_GLEAN_10021050, whose protein had insertions/deletions of five amino acids and substitutions of two amino acids in the conserved domain. This evidence indicates that these three genes are potentially involved in lignin biosynthesis in ramie fibers. The QTLs identified from this study provide a basis for the improvement of lignin content and fiber quality in ramie breeding. The characterization of the three candidate genes here will be helpful for the future clarification of their functions in ramie.

Phosphorylation of LTF1, an MYB Transcription Factor in Populus, Acts as a Sensory Switch Regulating Lignin Biosynthesis in Wood Cells

Lignin is specifically deposited in plant secondary cell walls, and initiation of lignin biosynthesis is regulated by a variety of developmental and environmental signals. However, the mechanisms governing the regulation of lignin biosynthesis remain to be elucidated. In this study, we identified a lignin biosynthesis-associated transcription factor (LTF) from Populus, LTF1, which binds the promoter of a key lignin biosynthetic gene encoding 4-coumarate-CoA ligase (4CL). We showed that LTF1 in its unphosphorylated state functions as a regulator restraining lignin biosynthesis. When LTF1 becomes phosphorylated by PdMPK6 in response to external stimuli such as wounding, it undergoes degradation through a proteasome pathway, resulting in activation of lignification. Expression of a phosphorylation-null mutant version of LTF1 led to stable protein accumulation and persistent attenuation of lignification in wood cells. Taken together, our study reveals a mechanism whereby LTF1 phosphorylation acts as a sensory switch to regulate lignin biosynthesis in response to environmental stimuli. The discovery of novel modulators and mechanisms modifying lignin biosynthesis has important implications for improving the utilization of cell-wall biomass.

OsMYB108 loss-of-function enriches p-coumaroylated and tricin lignin units in rice cell walls

Breeding approaches to enrich lignins in biomass could be beneficial to improving the biorefinery process because lignins increase biomass heating value and represent a potent source of valuable aromatic chemicals. However, despite the fact that grasses are promising lignocellulose feedstocks, limited information is yet available for molecular-breeding approaches to upregulate lignin biosynthesis in grass species. In this study, we generated lignin-enriched transgenic rice (Oryza sativa), a model grass species, via targeted mutagenesis of the transcriptional repressor OsMYB108 using CRISPR/Cas9-mediated genome editing. The OsMYB108-knockout rice mutants displayed increased expressions of lignin biosynthetic genes and enhanced lignin deposition in culm cell walls. Chemical and two-dimensional nuclear magnetic resonance (NMR) analyses revealed that the mutant cell walls were preferentially enriched in γ-p-coumaroylated and tricin lignin units, both of which are typical and unique components in grass lignins. NMR analysis also showed that the relative abundances of major lignin linkage types were altered in the OsMYB108 mutants.

What Makes the Wood? Exploring the Molecular Mechanisms of Xylem Acclimation in Hardwoods to an Ever-Changing Environment

Wood, also designated as secondary xylem, is the major structure that gives trees and other woody plants stability for upright growth and maintains the water supply from the roots to all other plant tissues. Over recent decades, our understanding of the cellular processes of wood formation (xylogenesis) has substantially increased. Plants as sessile organisms face a multitude of abiotic stresses, e.g., heat, drought, salinity and limiting nutrient availability that require them to adjust their wood structure to maintain stability and water conductivity. Because of global climate change, more drastic and sudden changes in temperature and longer periods without precipitation are expected to impact tree productivity in the near future. Thus, it is essential to understand the process of wood formation in trees under stress. Many traits, such as vessel frequency and size, fiber thickness and density change in response to different environmental stimuli. Here, we provide an overview of our current understanding of how abiotic stress factors affect wood formation on the molecular level focussing on the genes that have been identified in these processes.

Auxin-mediated Aux/IAA-ARF-HB signaling cascade regulates secondary xylem development in Populus

Wood development is strictly regulated by various phytohormones and auxin plays a central regulatory role in this process. However, how the auxin signaling is transducted in developing secondary xylem during wood formation in tree species remains unclear. Here, we identified an Aux/INDOLE-3-ACETIC ACID 9 (IAA9)-AUXIN RESPONSE FACTOR 5 (ARF5) module in Populus tomentosa as a key mediator of auxin signaling to control early developing xylem development. PtoIAA9, a canonical Aux/IAA gene, is predominantly expressed in vascular cambium and developing secondary xylem and induced by exogenous auxin. Overexpression of PtoIAA9m encoding a stabilized IAA9 protein significantly represses secondary xylem development in transgenic poplar. We further showed that PtoIAA9 interacts with PtoARF5 homologs via the C-terminal III/IV domains. The truncated PtoARF5.1 protein without the III/IV domains rescued defective phenotypes caused by PtoIAA9m. Expression analysis showed that the PtoIAA9-PtoARF5 module regulated the expression of genes associated with secondary vascular development in PtoIAA9m- and PtoARF5.1-overexpressing plants. Furthermore, PtoARF5.1 could bind to the promoters of two Class III homeodomain-leucine zipper (HD-ZIP III) genes, PtoHB7 and PtoHB8, to modulate secondary xylem formation. Taken together, our results suggest that the Aux/IAA9-ARF5 module is required for auxin signaling to regulate wood formation via orchestrating the expression of HD-ZIP III transcription factors in poplar.

Genome-Wide Association Studies to Improve Wood Properties: Challenges and Prospects

Wood formation is an excellent model system for quantitative trait analysis due to the strong associations between the transcriptional and metabolic traits that contribute to this complex process. Investigating the genetic architecture and regulatory mechanisms underlying wood formation will enhance our understanding of the quantitative genetics and genomics of complex phenotypic variation. Genome-wide association studies (GWASs) represent an ideal statistical strategy for dissecting the genetic basis of complex quantitative traits. However, elucidating the molecular mechanisms underlying many favorable loci that contribute to wood formation and optimizing GWAS design remain challenging in this omics era. In this review, we summarize the recent progress in GWAS-based functional genomics of wood property traits in major timber species such as Eucalyptus, Populus, and various coniferous species. We discuss several appropriate experimental designs for extensive GWAS in a given undomesticated tree population, such as omics-wide association studies and high-throughput phenotyping technologies. We also explain why more attention should be paid to rare allelic and major structural variation. Finally, we explore the potential use of GWAS for the molecular breeding of trees. Such studies will help provide an integrated understanding of complex quantitative traits and should enable the molecular design of new cultivars.

Transcriptome analysis of the genes regulating phytohormone and cellular patterning in Lagerstroemia plant architecture

Plant architecture is a popular research topic because plants with different growth habits that may generate economic or ornamental value are in great demand by orchards and nurseries. However, the molecular basis of the architecture of woody perennial plants is poorly understood due to the complexity of the phenotypic and regulatory relationships. Here, transcriptional profiling of dwarf and non-dwarf crapemyrtles was performed, and potential target genes were identified based on the phenotype, histology and phytohormone metabolite levels. An integrated analysis demonstrated that the internode length was explained mainly by cell number and secondarily by cell length and revealed important hormones in regulatory pathway of Lagerstroemia architecture. Differentially expressed genes (DEGs) involved in phytohormone pathways and cellular patterning regulation were analysed, and the regulatory relationships between these parameters were evaluated at the transcriptional level. Exogenous indole-3-acetic acid (IAA) and gibberellin A₄ (GA₄) treatments further indicated the pivotal role of auxin in cell division within the shoot apical meristem (SAM) and suggested an interaction between auxin and GA₄ in regulating the internode length of Lagerstroemia. These results provide insights for further functional genomic studies on the regulatory mechanisms underlying Lagerstroemia plant architecture and may improve the efficiency of woody plant molecular breeding.

Molecular cloning and characterization of a brassinosteriod biosynthesis-related gene PtoDWF4 from Populus tomentosa

Brassinosteroids (BRs) as steroid hormones play an important role in plant growth and development. However, little is known about how BRs affect secondary wall biosynthesis in woody plants. In this study, we cloned and characterized PtoDWF4, a homologus gene of Arabidopsis DWF4 encoding a cytochrome P450 protein, from Populus tomentosa. qRT-PCR analysis showed that PtoDWF4 was highly expressed in stems, especially in xylem. Overexpression of PtoDWF4 (PtoDWF4-OE) in poplar promoted growth rate and biomass yield, increased area and cell layers of xylem. Transgenic plants showed a significant increase in plant height and stem diameter compared with the wild type. In contrast, the CRISPR/Cas9-generated mutation of PtoDWF4 (PtoDWF4-KO) resulted in significantly decreased biomass production in transgenic plants. Further studies revealed that constitutive expression of PtoDWF4 up-regulated the expression of secondary cell wall (SCW) biosynthesis-related genes, whereas knock-out of PtoDWF4 down-regulated their expression. Quantitative analysis of cell wall components showed a significant increase in PtoDWF4-OE lines but a reduction in PtoDWF4-KO lines compared with wild-type plants. Taken together, our results indicate that PtoDWF4 plays a positive role in improving growth rate and elevating biomass production in poplar.

Identification of miRNAs Associated with Graft Union Development in Pecan [Carya illinoinensis (Wangenh.) K. Koch]

Pecan [Carya illinoinensis (Wangenh.) K. Koch] is a high-value fruit tree with a long juvenile period. The fruiting process of pecan seedlings can be largely accelerated through grafting. As non-coding small RNAs, plant miRNAs participate in various biological processes through negative regulation of gene expression. To reveal the roles of miRNAs in the graft union development of pecan, four small RNA libraries were constructed from the graft union at days 0, 8, 15, and 30 after grafting. A total of 47 conserved miRNAs belonging to 31 families and 39 novel miRNAs were identified. For identified miRNAs, 584 target genes were bioinformatically predicted, and 266 of them were annotated; 29 miRNAs (including 16 conserved and 13 novel miRNAs) were differentially expressed during the graft process. The expression profiles of 12 miRNA were further validated by quantitative reverse transcription PCR (qRT-PCR). In addition, qRT-PCR revealed that the expression levels of 3 target genes were negatively correlated with their corresponding miRNAs. We found that miRS26 might be involved in callus formation; miR156, miR160, miR164, miR166, and miRS10 might be associated with vascular bundle formation. These results indicate that the miRNA-mediated gene regulations play important roles in the graft union development of pecan.

Temporal analysis of Arabidopsis genes activated by Eucalyptus grandis NAC transcription factors associated with xylem fibre and vessel development

Secondary cell wall (SCW) deposition in Arabidopsis is regulated among others by NAC transcription factors, where SND1 chiefly initiates xylem fibre differentiation while VND6 controls metaxylem vessel SCW development, especially programmed cell death and wall patterning. The translational relevance of Arabidopsis SCW regulation theory and the utility of characterized transcription factors as modular synthetic biology tools for improving commercial fibre crops is unclear. We investigated inter-lineage gene activation dynamics for potential fibre and vessel differentiation regulators from the widely grown hardwood Eucalyptus grandis (Myrtales). EgrNAC26, a VND6 homolog, and EgrNAC61, an SND1 homolog, were transiently expressed in Arabidopsis mesophyll protoplasts in parallel to determine early and late (i.e. 7 and 14 hours post-transfection) gene targets. Surprisingly, across the time series EgrNAC26 activated only a subset of SCW-related transcription factors and biosynthetic genes activated by EgrNAC61, specializing instead in targeting vessel-specific wall pit and programmed cell death markers. Promoters of EgrNAC26 and EgrNAC61 both induced reporter gene expression in vessels of young Arabidopsis plants, with EgrNAC61 also conferring xylem- and cork cambium-preferential expression in Populus. Our results demonstrate partial conservation, with notable exceptions, of SND1 and VND6 homologs in Eucalyptus and a first report of cork cambium expression for EgrNAC61.

Co-expression network of transcription factors reveal ethylene-responsive element-binding factor as key regulator of wood phenotype in Eucalyptus tereticornis

Suitability of wood biomass for pulp production is dependent on the cellular architecture and composition of secondary cell wall. Presently, systems genetics approach is being employed to understand the molecular basis of trait variation and co-expression network analysis has enabled holistic understanding of complex trait such as secondary development. Transcription factors (TFs) are reported as key regulators of meristematic growth and wood formation. The hierarchical TF network is a multi-layered system which interacts with downstream structural genes involved in biosynthesis of cellulose, hemicelluloses and lignin. Several TFs have been associated with wood formation in tree species such as Populus, Eucalyptus, Picea and Pinus. However, TF-specific co-expression networks to understand the interaction between these regulators are not reported. In the present study, co-expression network was developed for TFs expressed during wood formation in Eucalyptus tereticornis and ethylene-responsive element-binding factor, EtERF2, was identified as the major hub transcript which co-expressed with other secondary cell wall biogenesis-specific TFs such as EtSND2, EtVND1, EtVND4, EtVND6, EtMYB70, EtGRAS and EtSCL8. This study reveals a probable role of ethylene in determining natural variation in wood properties in Eucalyptus species. Understanding this transcriptional regulation underpinning the complex bio-processing trait of wood biomass will complement the Eucalyptus breeding program through selection of industrially suitable phenotypes by marker-assisted selection.

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.

A group of Populus trichocarpa DUF231 proteins exhibit differential O-acetyltransferase activities toward xylan

Wood represents the most abundant biomass produced by plants and one of its major components is acetyl xylan. Acetylation in xylan can occur at O-2 or O-3 of a xylosyl residue, at both O-2 and O-3 of a xylosyl residue, and at O-3 of a xylosyl residue substituted at O-2 with glucuronic acid. Acetyltransferases responsible for the regiospecific acetylation of xylan in tree species have not yet been characterized. Here we report the biochemical characterization of twelve Populus trichocarpa DUF231-containing proteins, named PtrXOATs, for their roles in the regiospecific acetylation of xylan. The PtrXOAT genes were found to be differentially expressed in Populus organs and among them, PtrXOAT1, PtrXOAT2, PtrXOAT9 and PtrXOAT10 exhibited the highest level of expression in stems undergoing wood formation. Activity assays of recombinant proteins demonstrated that all twelve PtrXOAT proteins were able to transfer acetyl groups from acetyl CoA onto a xylohexaose acceptor with PtrXOAT1, PtrXOAT2, PtrXOAT3, PtrXOAT11 and PtrXOAT12 having the highest activity. Structural analysis of the PtrXOAT-catalyzed reaction products using ¹H NMR spectroscopy revealed that PtrXOAT1, PtrXAOT2 and PtrXOAT3 mediated 2-O- and 3-O-monoacetylation and 2,3-di-O-acetylation of xylosyl residues and PtrXOAT11 and PtrXOAT12 only catalyzed 2-O- and 3-O-monoacetylation of xylosyl residues. Of the twelve PtrXOATs, only PtrXOAT9 and PtrXOAT10 were capable of transferring acetyl groups onto the O-3 position of 2-O-glucuronic acid-substituted xylosyl residues. Furthermore, when expressed in the Arabidopsis eskimo1 mutant, PtrXOAT1, PtrXAOT2 and PtrXOAT3 were able to rescue the defects in xylan acetylation. Together, these results demonstrate that the twelve PtrXOATs are acetyltransferases with different roles in xylan acetylation in P. trichocarpa.

Genes and gene clusters related to genotype and drought-induced variation in saccharification potential, lignin content and wood anatomical traits in Populus nigra

Wood is a renewable resource that can be employed for the production of second generation biofuels by enzymatic saccharification and subsequent fermentation. Knowledge on how the saccharification potential is affected by genotype-related variation of wood traits and drought is scarce. Here, we used three Populus nigra L. genotypes from habitats differing in water availability to (i) investigate the relationships between wood anatomy, lignin content and saccharification and (ii) identify genes and co-expressed gene clusters related to genotype and drought-induced variation in wood traits and saccharification potential. The three poplar genotypes differed in wood anatomy, lignin content and saccharification potential. Drought resulted in reduced cambial activity, decreased vessel and fiber lumina, and increased the saccharification potential. The saccharification potential was unrelated to lignin content as well as to most wood anatomical traits. RNA sequencing of the developing xylem revealed that 1.5% of the analyzed genes were differentially expressed in response to drought, while 67% differed among the genotypes. Weighted gene correlation network analysis identified modules of co-expressed genes correlated with saccharification potential. These modules were enriched in gene ontology terms related to cell wall polysaccharide biosynthesis and modification and vesicle transport, but not to lignin biosynthesis. Among the most strongly saccharification-correlated genes, those with regulatory functions, especially kinases, were prominent. We further identified transcription factors whose transcript abundances differed among genotypes, and which were co-regulated with genes for biosynthesis and modifications of hemicelluloses and pectin. Overall, our study suggests that the regulation of pectin and hemicellulose metabolism is a promising target for improving wood quality of second generation bioenergy crops. The causal relationship of the identified genes and pathways with saccharification potential needs to be validated in further experiments.

Transcriptomic Analysis Provides Insights into Grafting Union Development in Pecan (Carya illinoinensis)

Pecan (Carya illinoinensis), as a popular nut tree, has been widely planted in China in recent years. Grafting is an important technique for its cultivation. For a successful grafting, graft union development generally involves the formation of callus and vascular bundles at the graft union. To explore the molecular mechanism of graft union development, we applied high throughput RNA sequencing to investigate the transcriptomic profiles of graft union at four timepoints (0 days, 8 days, 15 days, and 30 days) during the pecan grafting process. After de novo assembly, 83,693 unigenes were obtained, and 40,069 of them were annotated. A total of 12,180 differentially expressed genes were identified between by grafting. Genes involved in hormone signaling, cell proliferation, xylem differentiation, cell elongation, secondary cell wall deposition, programmed cell death, and reactive oxygen species (ROS) scavenging showed significant differential expression during the graft union developmental process. In addition, we found that the content of auxin, cytokinin, and gibberellin were accumulated at the graft unions during the grafting process. These results will aid in our understanding of successful grafting in the future.

PtoMYB170 positively regulates lignin deposition during wood formation in poplar and confers drought tolerance in transgenic Arabidopsis

Wood formation is a complex developmental process under multi-level transcriptional control executed by a large set of transcription factors. However, only limited members have been characterized to be key regulators of lignin biosynthesis in poplar. Here we report the conserved and unique functions of PtoMYB170, a transcription factor identified from Populus tomentosa (Chinese white poplar), in lignin deposition and drought tolerance in comparison with its duplicate paralog PtoMYB216. PtoMYB170 is preferentially expressed in young leaves and xylem tissues. Overexpression of PtoMYB170 in transgenic poplar plants resulted in stronger lignification and more thickened secondary wall in xylem compared with wild-type plants, whereas the CRISPR/Cas9-generated mutation of PtoMYB170 weakened lignin deposition, thereby leading to a more flexible and collapsed xylem phenotype. Transient expression experiments demonstrated that PtoMYB170 specifically activated the expression of lignin biosynthetic genes, consistent with the function of PtoMYB216. However, GUS staining assays revealed that PtoMYB170 was specifically expressed in guard cells of transgenic Arabidopsis while PtoMYB216 was not. Heterologous expression of PtoMYB170 in Arabidopsis enhanced stomatal closure in the dark and resulted in drought tolerance of the transgenic plants through reduced water loss, indicating a diversified role from PtoMYB216. These results revealed the PtoMYB170-dependent positive transcriptional regulation on lignin deposition in poplar and its coordinated function in enhancing drought tolerance by promoting dark-induced stomatal closure.

Radiocarbon Dating of an Olive Tree Cross-Section: New Insights on Growth Patterns and Implications for Age Estimation of Olive Trees

The age of living massive olive trees is often assumed to be between hundreds and even thousands of years. These estimations are usually based on the girth of the trunk and an extrapolation based on a theoretical annual growth rate. It is difficult to objectively verify these claims, as a monumental tree may not be cut down for analysis of its cross-section. In addition, the inner and oldest part of the trunk in olive trees usually rots, precluding the possibility of carting out radiocarbon analysis of material from the first years of life of the tree. In this work we present a cross-section of an olive tree, previously estimated to be hundreds of years old, which was cut down post-mortem in 2013. The cross-section was radiocarbon dated at numerous points following the natural growth pattern, which was made possible to observe by viewing the entire cross-section. Annual growth rate values were calculated and compared between different radii. The cross-section also revealed a nearly independent segment of growth, which would clearly offset any estimations based solely on girth calculations. Multiple piths were identified, indicating the beginning of branching within the trunk. Different radii were found to have comparable growth rates, resulting in similar estimates dating the piths to the 19th century. The estimated age of the piths represent a terminus ante quem for the age of the tree, as these are piths of separate branches. However, the tree is likely not many years older than the dated piths, and certainly not centuries older. The oldest radiocarbon-datable material in this cross-section was less than 200 years old, which is in agreement with most other radiocarbon dates of internal wood from living olive trees, rarely older than 300 years.

Characterization and Functional Analysis of the Poplar Pectate Lyase-Like Gene PtPL1-18 Reveal Its Role in the Development of Vascular Tissues

Pectin is a major component of plant cell walls, and the structure of pectin impacts on the properties of wood. Although we know that pectate lyase (PL, EC 4.2.2.2) has a major influence on the structure of pectin, our knowledge of Pectate lyase-like genes (PLL) in tree species remains limited. To better understand the characteristics of PLL genes in trees and to identify novel PLL genes that are potentially involved in the development of wood, we performed comprehensive analyses of gene structures, phylogenetic relationships, chromosomal locations, gene duplication events, conserved protein motifs, and gene expression patterns of 30 PLLs in Populus trichocarpa (PtPL1s). We performed an in silico gene expression profiling and quantitative real-time PCR analysis and found that most of the PtPL1 genes from subgroups Ia and Ib were highly expressed in xylem. PtPL1-18 from subgroup Ia was preferentially expressed in developing primary xylem and in xylem cells that were developing secondary walls. Overexpression of PtPL1-18 in poplar reduced plant growth and xylem development. Reduced secondary cell wall thickening and irregular xylem cells were observed in the transgenic trees, probably due to their lower pectin content. Although pectin is not a major component of plant secondary cell walls, our results are consistent with the PtPL1 genes performing important functions during wood formation.

Repression of BLADE-ON-PETIOLE genes by KNOX homeodomain protein BREVIPEDICELLUS is essential for differentiation of secondary xylem in Arabidopsis root

Repression of boundary genes by KNOTTED1-like homeodomain transcription factor BREVIPEDICELLUS promotes the differentiation of phase II secondary xylem in Arabidopsis roots. Plant growth and development relies on the activity of meristems. Boundaries are domains of restricted growth that separate forming organs and the meristem. Class I KNOX homeodomain transcription factors are important regulators of meristem maintenance. Members of this class including BREVIDICELLUS also called KNOTTED-LIKE FROM ARABIDOPSIS THALIANA1 (BP/KNAT1) fulfill this function in part by spatially regulating boundary genes. The vascular cambium is a lateral meristem that allows for radial expansion of organs during secondary growth. We show here that BP/KNAT1 repression of boundary genes plays a crucial role in root secondary growth. In particular, exclusion of BLADE-ON-PETIOLE1/2 (BOP1/2) and other members of this module from xylem is required for the differentiation of lignified fibers and vessels during the xylem expansion phase of root thickening. These data reveal a previously undiscovered role for boundary genes in the root and shed light on mechanisms controlling wood development in trees.

Proteomics analysis reveals the molecular mechanism underlying the transition from primary to secondary growth of poplar

Wood is the most important natural source of energy and also provides fuel and fiber. Considering the significant role of wood, it is critical to understand how wood is formed. Integration of knowledge about wood development at the cellular and molecular levels will allow more comprehensive understanding of this complex process. In the present study, we used a comparative proteomic approach to investigate the differences in protein profiles between primary and secondary growth in young poplar stems using tandem mass tag (TMT)-labeling. More than 10,816 proteins were identified, and, among these, 3106 proteins were differentially expressed during primary to secondary growth. Proteomic data were validated using a combination of histochemical staining, enzyme activity assays, and quantitative real-time PCR. Bioinformatics analysis revealed that these differentially expressed proteins are related to various metabolic pathways, mainly including signaling, phytohormones, cell cycle, cell wall, secondary metabolism, carbohydrate and energy metabolism, and protein metabolism as well as redox and stress pathways. This large proteomics dataset will be valuable for uncovering the molecular changes occurring during the transition from primary to secondary growth. Further, it provides new and accurate information for tree breeding to modify wood properties.

Identifying gene coexpression networks underlying the dynamic regulation of wood-forming tissues in Populus under diverse environmental conditions

Trees modify wood formation through integration of environmental and developmental signals in complex but poorly defined transcriptional networks, allowing trees to produce woody tissues appropriate to diverse environmental conditions. In order to identify relationships among genes expressed during wood formation, we integrated data from new and publically available datasets in Populus. These datasets were generated from woody tissue and include transcriptome profiling, transcription factor binding, DNA accessibility and genome-wide association mapping experiments. Coexpression modules were calculated, each of which contains genes showing similar expression patterns across experimental conditions, genotypes and treatments. Conserved gene coexpression modules (four modules totaling 8398 genes) were identified that were highly preserved across diverse environmental conditions and genetic backgrounds. Functional annotations as well as correlations with specific experimental treatments associated individual conserved modules with distinct biological processes underlying wood formation, such as cell-wall biosynthesis, meristem development and epigenetic pathways. Module genes were also enriched for DNase I hypersensitivity footprints and binding from four transcription factors associated with wood formation. The conserved modules are excellent candidates for modeling core developmental pathways common to wood formation in diverse environments and genotypes, and serve as testbeds for hypothesis generation and testing for future studies.