The evolution of insular woodiness.

Insular woodiness (IW)-the evolutionary transition from herbaceousness toward woodiness on islands-is one of the most iconic features of island floras. Since pioneering work by Darwin and Wallace, a number of drivers of IW have been proposed, such as 1) competition for sunlight requiring plants with taller and stronger woody stems and 2) drought favoring woodiness to safeguard root-to-shoot water transport. Alternatively, IW may be the indirect result of increased lifespan related to 3) a favorable aseasonal climate and/or 4) a lack of large native herbivores. However, information on the occurrence of IW is fragmented, hampering tests of these potential drivers. Here, we identify 1,097 insular woody species on 375 islands and infer at least 175 evolutionary transitions on 31 archipelagos, concentrated in six angiosperm families. Structural equation models reveal that the insular woody species richness on oceanic islands correlates with a favorable aseasonal climate, followed by increased drought and island isolation (approximating competition). When continental islands are also included, reduced herbivory pressure by large native mammals, increased drought, and island isolation are most relevant. Our results illustrate different trajectories leading to rampant convergent evolution toward IW and further emphasize archipelagos as natural laboratories of evolution, where similar abiotic or biotic conditions replicated evolution of similar traits.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera.

BACKGROUND: Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera. RESULTS: We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis. CONCLUSIONS: This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis.

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

Identification and characterization of functionally relevant SSR markers in natural Dalbergia odorifera populations.

BACKGROUND: Dalbergia odorifera is a rare and precious rosewood specie, which is valued for its amber tones, abstract figural patterns, and impermeability to water and insects. However, the information on genetic diversity and marker-assisted selection breeding of D. odorifera is still limited. Simple sequence repeat (SSR) markers are an ideal tool for genetic diversity analysis and marker-assisted molecular breeding for complex traits. RESULTS: Here, we have developed SSR markers within candidate genes and used them to explore the genetic diversity among D. odorifera germplasm resources. A total of 635 SSR loci were identified. The proportions of mono-, di- and tri-nucleotide repeat motifs were 52.28%, 22.99% and 21.42%, respectively. From these, a total of 114 SSR primers were synthesized, of which 24 SSR markers displayed polymorphism (polymorphic information content (PIC) > 0.25). Subsequently, these polymorphic markers were used for the genetic diversity analysis of 106 D. odorifera individuals from 11 natural populations. According to the genetic diversity analysis of D. odorifera natural populations, the average observed heterozygosity (Ho) was 0.500, the average expected heterozygosity (He) was 0.524, and the average Shannon's information index (I) was 0.946. These indicated that the natural populations had moderate genetic diversity. AMOVA analysis showed that 5% of the total variation was within the individuals of a population, whereas 95% of the variation was among the individuals of the populations, indicating a high degree of genetic variation between populations. On the basis of their genetic structures, these populations could be divided into four groups. CONCLUSIONS: Our study provides important experimental resources for genetic studies and assists in the program of molecular breeding of D. odorifera wood formation.

PagARGOS promotes low-lignin wood formation in poplar.

Wood formation, which occurs mainly through secondary xylem development, is important not only for supplying raw material for the 'ligno-chemical' industry but also for driving the storage of carbon. However, the complex mechanisms underlying the promotion of xylem formation remain to be elucidated. Here, we found that overexpression of Auxin-Regulated Gene involved in Organ Size (ARGOS) in hybrid poplar 84 K (Populus alba × Populus tremula var. glandulosa) enlarged organ size. In particular, PagARGOS promoted secondary growth of stems with increased xylem formation. To gain further insight into how PagARGOS regulates xylem development, we further carried out yeast two-hybrid screening and identified that the auxin transporter WALLS ARE THIN1 (WAT1) interacts with PagARGOS. Overexpression of PagARGOS up-regulated WAT1, activating a downstream auxin response promoting cambial cell division and xylem differentiation for wood formation. Moreover, overexpressing PagARGOS caused not only higher wood yield but also lower lignin content compared with wild-type controls. PagARGOS is therefore a potential candidate gene for engineering fast-growing and low-lignin trees with improved biomass production.

Systems genetic analysis of lignin biosynthesis in Populus tremula.

The genetic control underlying natural variation in lignin content and composition in trees is not fully understood. We performed a systems genetic analysis to uncover the genetic regulation of lignin biosynthesis in a natural 'SwAsp' population of aspen (Populus tremula) trees. We analyzed gene expression by RNA sequencing (RNA-seq) in differentiating xylem tissues, and lignin content and composition using Pyrolysis-GC-MS in mature wood of 268 trees from 99 genotypes. Abundant variation was observed for lignin content and composition, and genome-wide association study identified proteins in the pentose phosphate pathway and arabinogalactan protein glycosylation among the top-ranked genes that are associated with these traits. Variation in gene expression and the associated genetic polymorphism was revealed through the identification of 312 705 local and 292 003 distant expression quantitative trait loci (eQTL). A co-expression network analysis suggested modularization of lignin biosynthesis and novel functions for the lignin-biosynthetic CINNAMYL ALCOHOL DEHYDROGENASE 2 and CAFFEOYL-CoA O-METHYLTRANSFERASE 3. PHENYLALANINE AMMONIA LYASE 3 was co-expressed with HOMEOBOX PROTEIN 5 (HB5), and the role of HB5 in stimulating lignification was demonstrated in transgenic trees. The systems genetic approach allowed linking natural variation in lignin biosynthesis to trees´ responses to external cues such as mechanical stimulus and nutrient availability.

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.

Analyses of high spatial resolution datasets identify genes associated with multi-layered secondary cell wall thickening in Pinus bungeana.

BACKGROUND AND AIMS: Secondary cell wall (SCW) thickening is a major cellular developmental stage determining wood structure and properties. Although the molecular regulation of cell wall deposition during tracheary element differentiation has been well established in primary growth systems, less is known about the gene regulatory processes involved in the multi-layered SCW thickening of mature trees. METHODS: Using third-generation [long-read single-molecule real-time (SMRT)] and second-generation [short-read sequencing by synthesis (SBS)] sequencing methods, we established a Pinus bungeana transcriptome resource with comprehensive functional and structural annotation for the first time. Using these approaches, we generated high spatial resolution datasets for the vascular cambium, xylem expansion regions, early SCW thickening, late SCW thickening and mature xylem tissues of 71-year-old Pinus bungeana trees. KEY RESULTS: A total of 79 390 non-redundant transcripts, 31 808 long non-coding RNAs and 5147 transcription factors were annotated and quantified in different xylem tissues at all growth and differentiation stages. Furthermore, using this high spatial resolution dataset, we established a comprehensive transcriptomic profile and found that members of the NAC, WRKY, SUS, CESA and LAC gene families are major players in early SCW formation in tracheids, whereas members of the MYB and LBD transcription factor families are highly expressed during late SCW thickening. CONCLUSIONS: Our results provide new molecular insights into the regulation of multi-layered SCW thickening in conifers. The high spatial resolution datasets provided can serve as important gene resources for improving softwoods.

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

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

Wood of trees: Cellular structure, molecular formation, and genetic engineering.

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

PagMYB128 regulates secondary cell wall formation by direct activation of cell wall biosynthetic genes during wood formation in poplar.

The biosynthesis of cellulose, lignin, and hemicelluloses in plant secondary cell walls (SCWs) is regulated by a hierarchical transcriptional regulatory network. This network features orthologous transcription factors shared between poplar and Arabidopsis, highlighting a foundational similarity in their genetic regulation. However, knowledge on the discrepant behavior of the transcriptional-level molecular regulatory mechanisms between poplar and Arabidopsis remains limited. In this study, we investigated the function of PagMYB128 during wood formation and found it had broader impacts on SCW formation compared to its Arabidopsis ortholog, AtMYB103. Transgenic poplar trees overexpressing PagMYB128 exhibited significantly enhanced xylem development, with fiber cells and vessels displaying thicker walls, and an increase in the levels of cellulose, lignin, and hemicelluloses in the wood. In contrast, plants with dominant repression of PagMYB128 demonstrated the opposite phenotypes. RNA sequencing and reverse transcription - quantitative polymerase chain reaction showed that PagMYB128 could activate SCW biosynthetic gene expression, and chromatin immunoprecipitation along with yeast one-hybrid, and effector-reporter assays showed this regulation was direct. Further analysis revealed that PagSND1 (SECONDARY WALL-ASSOCIATED NAC-DOMAIN PROTEIN1) directly regulates PagMYB128 but not cell wall metabolic genes, highlighting the pivotal role of PagMYB128 in the SND1-driven regulatory network for wood development, thereby creating a feedforward loop in SCW biosynthesis.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing.

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

Toward a better understanding of angiosperm xylogenesis: a new method for a cellular approach.

The kinetics of wood formation in angiosperms are largely unknown because their complex xylem anatomy precludes using the radial position of vessels and fibers to infer their time of differentiation. We analyzed xylogenesis in ring-porous ash (Fraxinus angustifolia) and diffuse-porous beech (Fagus sylvatica) over 1 yr and proposed a novel procedure to assess the period of vessel and fiber enlargement using a referential radial file (RRF). Our approach captured the dynamics of wood formation and provided a robust estimation of the kinetics of vessel and fiber enlargement. In beech, fibers and vessels had a similar duration of enlargement, decreasing from 14 to 5 d between April and July. In ash, wide vessels formed in April enlarged at a rate of 27 × 103 µm2 d-1 , requiring half the time of contemporary fibers (6 vs 12 d), and less time than the narrower vessels (14 d) formed in May. These findings reveal distinct cell-type-dependent mechanisms for differentiation in diffuse-porous and ring-porous trees, enhancing our understanding of angiosperm wood cell kinetics. Our approach presents an effective method for investigating angiosperm wood formation and provides a more accurate representation of vessel and fiber morphogenesis in wood formation models.

Photoperiod and temperature as dominant environmental drivers triggering secondary growth resumption in Northern Hemisphere conifers.

Wood formation consumes around 15% of the anthropogenic CO2 emissions per year and plays a critical role in long-term sequestration of carbon on Earth. However, the exogenous factors driving wood formation onset and the underlying cellular mechanisms are still poorly understood and quantified, and this hampers an effective assessment of terrestrial forest productivity and carbon budget under global warming. Here, we used an extensive collection of unique datasets of weekly xylem tissue formation (wood formation) from 21 coniferous species across the Northern Hemisphere (latitudes 23 to 67°N) to present a quantitative demonstration that the onset of wood formation in Northern Hemisphere conifers is primarily driven by photoperiod and mean annual temperature (MAT), and only secondarily by spring forcing, winter chilling, and moisture availability. Photoperiod interacts with MAT and plays the dominant role in regulating the onset of secondary meristem growth, contrary to its as-yet-unquantified role in affecting the springtime phenology of primary meristems. The unique relationships between exogenous factors and wood formation could help to predict how forest ecosystems respond and adapt to climate warming and could provide a better understanding of the feedback occurring between vegetation and climate that is mediated by phenology. Our study quantifies the role of major environmental drivers for incorporation into state-of-the-art Earth system models (ESMs), thereby providing an improved assessment of long-term and high-resolution observations of biogeochemical cycles across terrestrial biomes.

Exploring the Seasonal Dynamics and Molecular Mechanism of Wood Formation in Gymnosperm Trees.

Forests, comprising 31% of the Earth's surface, play pivotal roles in regulating the carbon, water, and energy cycles. Despite being far less diverse than angiosperms, gymnosperms account for over 50% of the global woody biomass production. To sustain growth and development, gymnosperms have evolved the capacity to sense and respond to cyclical environmental signals, such as changes in photoperiod and seasonal temperature, which initiate growth (spring and summer) and dormancy (fall and winter). Cambium, the lateral meristem responsible for wood formation, is reactivated through a complex interplay among hormonal, genetic, and epigenetic factors. Temperature signals perceived in early spring induce the synthesis of several phytohormones, including auxins, cytokinins, and gibberellins, which in turn reactivate cambium cells. Additionally, microRNA-mediated genetic and epigenetic pathways modulate cambial function. As a result, the cambium becomes active during the summer, resulting in active secondary xylem (i.e., wood) production, and starts to become inactive in autumn. This review summarizes and discusses recent findings regarding the climatic, hormonal, genetic, and epigenetic regulation of wood formation in gymnosperm trees (i.e., conifers) in response to seasonal changes.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis.

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa.

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

The chromosome-scale genomes of Dipterocarpus turbinatus and Hopea hainanensis (Dipterocarpaceae) provide insights into fragrant oleoresin biosynthesis and hardwood formation.

Dipterocarpaceae are typical tropical plants (dipterocarp forests) that are famous for their high economic value because of their production of fragrant oleoresins, top-quality timber and usage in traditional Chinese medicine. Currently, the lack of Dipterocarpaceae genomes has been a limiting factor to decipher the fragrant oleoresin biosynthesis and gain evolutionary insights into high-quality wood formation in Dipterocarpaceae. We generated chromosome-level genome assemblies for two representative Dipterocarpaceae species viz. Dipterocarpus turbinatus Gaertn. f. and Hopea hainanensis Merr. et Chun. Our whole-genome duplication (WGD) analysis revealed that Dipterocarpaceae underwent a shared WGD event, which showed significant impacts on increased copy numbers of genes related to the biosynthesis of terpene, BAHD acyltransferases, fatty acid and benzenoid/phenylpropanoid, which probably confer to the formation of their characteristic fragrant oleoresin. Additionally, compared with common soft wood plants, the expansion of gene families was also found to be associated with wood formation, such as in CESA (cellulose synthase), CSLE (cellulose synthase-like protein E), laccase and peroxidase in Dipterocarpaceae genomes, which might also contribute to the formation of harder, stronger and high-density timbers. Finally, an integrative analysis on a combination of genomic, transcriptomic and metabolic data from different tissues provided further insights into the molecular basis of fragrant oleoresins biosynthesis and high-quality wood formation of Dipterocarpaceae. Our study contributes the first two representative genomes for Dipterocarpaceae, which are valuable genetic resources for further researches on the fragrant oleoresins and superior-quality timber, genome-assisted breeding and improvement, and conservation biology of this family.

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

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

Aspartic proteases modulate programmed cell death and secondary cell wall synthesis during wood formation in poplar.

Programmed cell death (PCD) is essential for wood development in trees. However, the determination of crucial factors involved in xylem PCD of wood development is still lacking. Here, two Populus trichocarpa typical aspartic protease (AP) genes, AP17 and AP45, modulate xylem maturation, especially fibre PCD, during wood formation. AP17 and AP45 were dominantly expressed in the fibres of secondary xylem, as suggested by GUS expression in APpro::GUS transgenic plants. Cas9/gRNA-induced AP17 or AP45 mutants delayed secondary xylem fibre PCD, and ap17ap45 double mutants showed more serious defects. Conversely, AP17 overexpression caused premature PCD in secondary xylem fibres, indicating a positive modulation in wood fibre PCD. Loss of AP17 and AP45 did not alter wood fibre wall thickness, whereas the ap17ap45 mutants showed a low lignin content in wood. However, AP17 overexpression led to a significant decrease in wood fibre wall thickness and lignin content, revealing the involvement in secondary cell wall synthesis during wood formation. In addition, the ap17ap45 mutant and AP17 overexpression plants resulted in a significant increase in saccharification yield in wood. Overall, AP17 and AP45 are crucial modulators in xylem maturation during wood development, providing potential candidate genes for engineering lignocellulosic wood for biofuel utilization.