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

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.

The microRNA7833-AUX6 module plays a critical role in wood development by modulating cellular auxin influx in Populus tomentosa

The critical role of auxin on secondary vascular development in woody plants has been demonstrated. The concentration gradient of endogenous indole-3-acetic acid and the cellular and molecular pathways contributing to the auxin-directed vascular organization and wood growth have been uncovered in recent decades. However, our understanding of the roles and regulations of auxin influx in wood formation in trees remains limited. Here, we reported that a microRNA, miR7833, participates in the negative regulation of stem cambial cell division and secondary xylem development in Populus tomentosa. The miR7833 is mainly expressed in the vascular cambium during stem radical growth and specifically targets and represses two AUX/LAX family auxin influx carriers, AUX5 and AUX6, in poplar. We further revealed that poplar AUX6, the most abundant miR7833 target in the stem, is preferentially enriched in the developing xylem and is a positive regulator for cell division and differentiation events during wood formation. Moreover, inhibition of auxin influx carriers by 1-naphthoxyacetic acids abolished the regulatory effects of miR7833 and AUX6 on secondary xylem formation in poplar. Our results revealed the essential roles of the miR7833-AUX6 module in regulating cellular events in secondary xylem development and demonstrated an auxin influx-dependent mechanism for wood formation in poplar.

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

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

PagUNE12 encodes a basic helix-loop-helix transcription factor that regulates the development of secondary vascular tissue in poplar

Secondary growth in woody plants generates new cells and tissues via the activity of the vascular cambium and drives the radial expansion of stems and roots. It is regulated by a series of endogenous factors, especially transcription factors. Here, we cloned the basic helix-loop-helix (bHLH) transcription factor gene UNFERTILIZED EMBRYO SAC12 (UNE12) from poplar (Populus alba × Populus glandulosa Uyeki) and used biochemical, molecular, and cytological assays to investigate the biological functions and regulatory mechanism of PagUNE12. PagUNE12 mainly localized in the nucleus and possessed transcriptional activation activity. It was widely expressed in vascular tissues, including primary phloem and xylem and secondary phloem and xylem. Poplar plants overexpressing PagUNE12 showed significantly reduced plant height, shorter internodes, and curled leaves compared with wild-type plants. Optical microscopy and transmission electron microscopy revealed that overexpressing PagUNE12 promoted secondary xylem development, with thicker secondary cell walls than wild-type poplar. Fourier transform infrared spectroscopy, confocal Raman microscopy, and 2D Heteronuclear Single Quantum Correlation analysis indicated that these plants also had increased lignin contents, with a lower relative abundance of syringyl lignin units and a higher relative abundance of guaiacyl lignin units. Therefore, overexpressing PagUNE12 promoted secondary xylem development and increased the lignin contents of secondary xylem in poplar, suggesting that this gene could be used to improve wood quality in the future.

Single-cell transcriptomics unveils xylem cell development and evolution

BACKGROUND

Xylem, the most abundant tissue on Earth, is responsible for lateral growth in plants. Typical xylem has a radial system composed of ray parenchyma cells and an axial system of fusiform cells. In most angiosperms, fusiform cells comprise vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed by tracheids in other vascular plants such as gymnosperms. Little is known about the developmental programs and evolutionary relationships of these xylem cell types.

RESULTS

Through both single-cell and laser capture microdissection transcriptomic profiling, we determine the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Based on cross-species analyses of single-cell clusters and overlapping trajectories, we reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages, whereas the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed trait, exhibit strong transcriptomic similarity to vessel elements rather than libriform fibers.

CONCLUSIONS

This evo-devo framework provides a comprehensive understanding of the formation of xylem cell lineages across multiple plant species spanning over a hundred million years of evolutionary history.

The Inducible Accumulation of Cell Wall-Bound p-Hydroxybenzoates Is Involved in the Regulation of Gravitropic Response of Poplar

Lignin in Populus species is acylated with p-hydroxybenzoate. Monolignol p-hydroxybenzoyltransferase 1 (PHBMT1) mediates p-hydroxybenzoylation of sinapyl alcohol, eventually leading to the modification of syringyl lignin subunits. Angiosperm trees upon gravistimulation undergo the re-orientation of their growth along with the production of specialized secondary xylem, i.e., tension wood (TW), that generates tensile force to pull the inclined stem or leaning branch upward. Sporadic evidence suggests that angiosperm TW contains relatively a high percentage of syringyl lignin and lignin-bound p-hydroxybenzoate. However, whether such lignin modification plays a role in gravitropic response remains unclear. By imposing mechanical bending and/or gravitropic stimuli to the hybrid aspens in the wild type (WT), lignin p-hydroxybenzoate deficient, and p-hydroxybenzoate overproduction plants, we examined the responses of plants to gravitropic/mechanical stress and their cell wall composition changes. We revealed that mechanical bending or gravitropic stimulation not only induced the overproduction of crystalline cellulose fibers and increased the relative abundance of syringyl lignin, but also significantly induced the expression of PHBMT1 and the increased accumulation of p-hydroxybenzoates in TW. Furthermore, we found that although disturbing lignin-bound p-hydroxybenzoate accumulation in the PHBMT1 knockout and overexpression (OE) poplars did not affect the major chemical composition shifts of the cell walls in their TW as occurred in the WT plants, depletion of p-hydroxybenzoates intensified the gravitropic curving of the plantlets in response to gravistimulation, evident with the enhanced stem secant bending angle. By contrast, hyperaccumulation of p-hydroxybenzoates mitigated gravitropic response. These data suggest that PHBMT1-mediated lignin modification is involved in the regulation of poplar gravitropic response and, likely by compromising gravitropism and/or enhancing autotropism, negatively coordinates the action of TW cellulose fibers to control the poplar wood deformation and plant growth.

Overexpression of Plasma Membrane H+-ATPase in Guard Cells Enhances Light-Induced Stomatal Opening, Photosynthesis, and Plant Growth in Hybrid Aspen

Stomata in the plant epidermis open in response to light and regulate CO₂ uptake for photosynthesis and transpiration for uptake of water and nutrients from roots. Light-induced stomatal opening is mediated by activation of the plasma membrane (PM) H⁺-ATPase in guard cells. Overexpression of PM H⁺-ATPase in guard cells promotes light-induced stomatal opening, enhancing photosynthesis and growth in Arabidopsis thaliana. In this study, transgenic hybrid aspens overexpressing Arabidopsis PM H⁺-ATPase (AHA2) in guard cells under the strong guard cell promoter Arabidopsis GC1 (AtGC1) showed enhanced light-induced stomatal opening, photosynthesis, and growth. First, we confirmed that AtGC1 induces GUS expression specifically in guard cells in hybrid aspens. Thus, we produced AtGC1::AHA2 transgenic hybrid aspens and confirmed expression of AHA2 in AtGC1::AHA2 transgenic plants. In addition, AtGC1::AHA2 transgenic plants showed a higher PM H⁺-ATPase protein level in guard cells. Analysis using a gas exchange system revealed that transpiration and the photosynthetic rate were significantly increased in AtGC1::AHA2 transgenic aspen plants. AtGC1::AHA2 transgenic plants showed a>20% higher stem elongation rate than the wild type (WT). Therefore, overexpression of PM H⁺-ATPase in guard cells promotes the growth of perennial woody plants.

Transcriptional landscape of highly lignified poplar stems at single-cell resolution

BACKGROUND

Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood.

RESULTS

Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees.

CONCLUSIONS

These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees.

Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus

High-throughput single-cell RNA sequencing (scRNA-seq) has advantages over traditional RNA-seq to explore spatiotemporal information on gene dynamic expressions in heterogenous tissues. We performed Drop-seq, a method for the dropwise sequestration of single cells for sequencing, on protoplasts from the differentiating xylem of Populus alba × Populus glandulosa. The scRNA-seq profiled 9,798 cells, which were grouped into 12 clusters. Through characterization of differentially expressed genes in each cluster and RNA in situ hybridizations, we identified vessel cells, fiber cells, ray parenchyma cells and xylem precursor cells. Diffusion pseudotime analyses revealed the differentiating trajectory of vessels, fiber cells and ray parenchyma cells and indicated a different differentiation process between vessels and fiber cells, and a similar differentiation process between fiber cells and ray parenchyma cells. We identified marker genes for each cell type (cluster) and key candidate regulators during developmental stages of xylem cell differentiation. Our study generates a high-resolution expression atlas of wood formation at the single cell level and provides valuable information on wood formation.

Molecular self-organization of wood lignin-carbohydrate matrix

The analysis of the state of research on the chemical composition, functional nature and structure of the main components of the lignin-carbohydrate matrix allows considering the wood substance as a thermodynamically self-organizing nanobiocomposite system. Features of biosynthesis of the wood matrix main biopolymers, the formation of their functional nature and structure determine the complex hierarchical organization of cell walls. The supramolecular level of biosynthesis considers the interaction of cell wall components. On the one hand, these are questions of dynamics of cell walls synthesis and processes of self-organization that control the formation of chaotic objects of biological origin; on the other hand, it is the question of thermodynamic compatibility of plant tissue components. Various models of structural organization are currently being considered, focusing on various features (biological, chemical, structural) of wood substance. At the same time, the lignin-carbohydrate matrix is a three-component system of natural polymers: lignin-hemicelluloses-cellulose, the state of which is described by specific values of thermodynamic parameters that characterize the degree of its stability. The new approach proposed in this paper allows considering the plant lignin-carbohydrate matrix from the standpoint of physical chemistry of polymer as quasi-equilibrium, thermodynamically limited ordered system of biopolymers. Thus, the biochemical processes of synthesis and self-organization lead to the formation of a complex multicomponent system of wood substance, considered as a nanobiocomposite. This determines the need to study the applicability of the fundamental cycle "structure-functional nature-properties" from the standpoint of physical chemistry of biopolymers both for the investigation of plant objects and for the development of modern technologies for complex processing based on the principles of "green chemistry".

Disparate Dynamics of Gene Body and cis-Regulatory Element Evolution Illustrated for the Senescence-Associated Cysteine Protease Gene SAG12 of Plants

Gene regulation networks precisely orchestrate the expression of genes that are closely associated with defined physiological and developmental processes such as leaf senescence in plants. The Arabidopsis thaliana senescence-associated gene 12 (AtSAG12) encodes a cysteine protease that is (i) involved in the degradation of chloroplast proteins and (ii) almost exclusively expressed during senescence. Transcription factors, such as WRKY53 and WRKY45, bind to W-boxes in the promoter region of AtSAG12 and play key roles in its activation. Other transcription factors, such as bZIPs, might have accessory functions in their gene regulation, as several A-boxes have been identified and appear to be highly overrepresented in the promoter region compared to the whole genome distribution but are not localized within the regulatory regions driving senescence-associated expression. To address whether these two regulatory elements exhibiting these different properties are conserved in other closely related species, we constructed phylogenetic trees of the coding sequences of orthologs of AtSAG12 and screened their respective 2000 bp promoter regions for the presence of conserved cis-regulatory elements, such as bZIP and WRKY binding sites. Interestingly, the functional relevant upstream located W-boxes were absent in plant species as closely related as Arabidopsis lyrata, whereas an A-box cluster appeared to be conserved in the Arabidopsis species but disappeared in Brassica napus. Several orthologs were present in other species, possibly because of local or whole genome duplication events, but with distinct cis-regulatory sites in different locations. However, at least one gene copy in each family analyzed carried one W-box and one A-box in its promoter. These gene differences in SAG12 orthologs are discussed in the framework of cis- and trans-regulatory factors, of promoter and gene evolution, of genetic variation, and of the enhancement of the adaptability of plants to changing environmental conditions.

Wood transcriptome analysis of Pinus densiflora identifies genes critical for secondary cell wall formation and NAC transcription factors involved in tracheid formation

Although conifers have significant ecological and economic value, information on transcriptional regulation of wood formation in conifers is still limited. Here, to gain insight into secondary cell wall (SCW) biosynthesis and tracheid formation in conifers, we performed wood tissue-specific transcriptome analyses of Pinus densiflora (Korean red pine) using RNA sequencing. In addition, to obtain full-length transcriptome information, PacBio single molecule real-time iso-sequencing was carried out using RNAs from 28 tissues of P. densiflora. Subsequent comparative tissue-specific transcriptome analysis successfully pinpointed critical genes encoding key proteins involved in biosynthesis of the major secondary wall components (cellulose, galactoglucomannan, xylan and lignin). Furthermore, we predicted a total of 62 NAC (NAM, ATAF1/2 and CUC2) family transcription factor members and identified seven PdeNAC genes preferentially expressed in developing xylem tissues in P. densiflora. Protoplast-based transcriptional activation analysis found that four PdeNAC genes, homologous to VND, NST and SND/ANAC075, upregulated GUS activity driven by an SCW-specific cellulose synthase promoter. Consistently, transient overexpression of the four PdeNACs induced xylem vessel cell-like SCW deposition in both tobacco (Nicotiana benthamiana) and Arabidopsis leaves. Taken together, our data provide a foundation for further research to unravel transcriptional regulation of wood formation in conifers, especially SCW formation and tracheid differentiation.

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.

Thermophilic microbial deconstruction and conversion of natural and transgenic lignocellulose

The potential to convert renewable plant biomasses into fuels and chemicals by microbial processes presents an attractive, less environmentally intense alternative to conventional routes based on fossil fuels. This would best be done with microbes that natively deconstruct lignocellulose and concomitantly form industrially relevant products, but these two physiological and metabolic features are rarely and simultaneously observed in nature. Genetic modification of both plant feedstocks and microbes can be used to increase lignocellulose deconstruction capability and generate industrially relevant products. Separate efforts on plants and microbes are ongoing, but these studies lack a focus on optimal, complementary combinations of these disparate biological systems to obtain a convergent technology. Improving genetic tools for plants have given rise to the generation of low-lignin lines that are more readily solubilized by microorganisms. Most focus on the microbiological front has involved thermophilic bacteria from the genera Caldicellulosiruptor and Clostridium, given their capacity to degrade lignocellulose and to form bio-products through metabolic engineering strategies enabled by ever-improving molecular genetics tools. Bioengineering plant properties to better fit the deconstruction capabilities of candidate consolidated bioprocessing microorganisms has potential to achieve the efficient lignocellulose deconstruction needed for industrial relevance.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood

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

MYB Transcription Factors and Its Regulation in Secondary Cell Wall Formation and Lignin Biosynthesis during Xylem Development

The secondary wall is the main part of wood and is composed of cellulose, xylan, lignin, and small amounts of structural proteins and enzymes. Lignin molecules can interact directly or indirectly with cellulose, xylan and other polysaccharide molecules in the cell wall, increasing the mechanical strength and hydrophobicity of plant cells and tissues and facilitating the long-distance transportation of water in plants. MYBs (v-myb avian myeloblastosis viral oncogene homolog) belong to one of the largest superfamilies of transcription factors, the members of which regulate secondary cell-wall formation by promoting/inhibiting the biosynthesis of lignin, cellulose, and xylan. Among them, MYB46 and MYB83, which comprise the second layer of the main switch of secondary cell-wall biosynthesis, coordinate upstream and downstream secondary wall synthesis-related transcription factors. In addition, MYB transcription factors other than MYB46/83, as well as noncoding RNAs, hormones, and other factors, interact with one another to regulate the biosynthesis of the secondary wall. Here, we discuss the biosynthesis of secondary wall, classification and functions of MYB transcription factors and their regulation of lignin polymerization and secondary cell-wall formation during wood formation.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

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.

Expression Analysis of Key Auxin Biosynthesis, Transport, and Metabolism Genes of Betula pendula with Special Emphasis on Figured Wood Formation in Karelian Birch

Auxin status in woody plants is believed to be a critical factor for the quantity and quality of the wood formed. It has been previously demonstrated that figured wood formation in Karelian birch (Betula pendula Roth var. carelica (Merckl.) Hämet-Ahti) is associated with a reduced auxin level and elevated sugar content in the differentiating xylem, but the molecular mechanisms of the abnormal xylogenesis remained largely unclear. We have identified genes involved in auxin biosynthesis (Yucca), polar auxin transport (PIN) and the conjugation of auxin with amino acids (GH3) and UDP-glucose (UGT84B1) in the B. pendula genome, and analysed their expression in trunk tissues of trees differing in wood structure. Almost all the investigated genes were overexpressed in Karelian birch trunks. Although Yucca genes were overexpressed, trunk tissues in areas developing figured grain had traits of an auxin-deficient phenotype. Overexpression of GH3s and UGT84B1 appears to have a greater effect on figured wood formation. Analysis of promoters of the differentially expressed genes revealed a large number of binding sites with various transcription factors associated with auxin and sugar signalling. These data agree with the hypothesis that anomalous figured wood formation in Karelian birch may be associated with the sugar induction of auxin conjugation.

Characterization of the UDP-glycosyltransferase UGT72 Family in Poplar and Identification of Genes Involved in the Glycosylation of Monolignols

Monolignols are the building blocks for lignin polymerization in the apoplastic domain. Monolignol biosynthesis, transport, storage, glycosylation, and deglycosylation are the main biological processes partaking in their homeostasis. In Arabidopsis thaliana, members of the uridine diphosphate-dependent glucosyltransferases UGT72E and UGT72B subfamilies have been demonstrated to glycosylate monolignols. Here, the poplar UGT72 family, which is clustered into four groups, was characterized: Group 1 UGT72AZ1 and UGT72AZ2, homologs of Arabidopsis UGT72E1-3, as well as group 4 UGT72B37 and UGT72B39, homologs of Arabidopsis UGT72B1-3, glycosylate monolignols. In addition, promoter-GUS analyses indicated that poplar UGT72 members are expressed within vascular tissues. At the subcellular level, poplar UGT72s belonging to group 1 and group 4 were found to be associated with the nucleus and the endoplasmic reticulum. However, UGT72A2, belonging to group 2, was localized in bodies associated with chloroplasts, as well as possibly in chloroplasts. These results show a partial conservation of substrate recognition between Arabidopsis and poplar homologs, as well as divergent functions between different groups of the UGT72 family, for which the substrates remain unknown.

Overexpression of a developing xylem cDNA library in transgenic poplar generates high mutation rate specific to wood formation

We investigated feasibility of the Full-length complementary DNA OvereXpression (FOX) system as a mutagenesis approach in poplar, using developing xylem tissue. The main goal was to assess the overall mutation rate and if the system will increase instances of mutants affected in traits linked to the xylem tissue. Indeed, we found a high mutation rate of 17.7%, whereas 80% of all mutants were significantly affected in cellulose, lignin and/or hemicellulose. Cell wall biosynthesis is a major process occurring during xylem development. Enrichment of mutants affected in cell wall composition suggests that the tissue source for the FOX library influenced the occurrence of mutants affected in a trait linked to this tissue. Additionally, we found that FLcDNAs from mutants affected in cell wall composition were homologous to genes known to be involved in cell wall biosynthesis and most recovered FLcDNAs corresponded to genes whose native expression was highest in xylem. We characterized in detail a mutant line with increased diameter. The phenotype was caused by a poplar homolog of LONELY GUY 1 (LOG1), which encodes an enzyme in cytokinin biosynthesis and significantly increased xylem proliferation. The causative role of LOG1 in the observed phenotype was further reaffirmed by elevated cytokinin concentration in the mutant and recapitulation overexpression experiment wherein multiple independent lines phenocopied the original FOX mutant. Our experiments show that the FOX approach can be efficiently used for gene discovery and molecular interrogation of traits specific to woody perennial growth and development.

Molecular Changes Concomitant with Vascular System Development in Mature Galls Induced by Root-Knot Nematodes in the Model Tree Host Populus tremula × P. alba

One of the most striking features occurring in the root-knot nematode Meloidogyne incognita induced galls is the reorganization of the vascular tissues. During the interaction of the model tree species Populus and M. incognita, a pronounced xylem proliferation was previously described in mature galls. To better characterise changes in expression of genes possibly involved in the induction and the formation of the de novo developed vascular tissues occurring in poplar galls, a comparative transcript profiling of 21-day-old galls versus uninfected root of poplar was performed. Genes coding for transcription factors associated with procambium maintenance and vascular differentiation were shown to be differentially regulated, together with genes partaking in phytohormones biosynthesis and signalling. Specific signatures of transcripts associated to primary cell wall biosynthesis and remodelling, as well as secondary cell wall formation (cellulose, xylan and lignin) were revealed in the galls. Ultimately, we show that molecules derived from the monolignol and salicylic acid pathways and related to secondary cell wall deposition accumulate in mature galls.

Symptoms and yield loss caused by rice stripe mosaic virus

Background

Rice stripe mosaic virus (RSMV) is a tentative new Cytorhabdovirus species in family Rhabdoviridae transmitted by the leafhopper Recilia dorsalis. Although the virus was first detected in southern China in 2015, few studies have investigated rice symptoms and yield losses caused by RSMV infection.

Methods

In this study, we observed and systematically compared symptoms of three virally infected, representative varieties of indica, hybrid and japonica rice and determined the yield parameters of the artificially inoculated plants.

Results

The three RSMV-infected cultivated rice varieties exhibited slight dwarfing, striped mosaicism, stiff, crinkled or even twisted leaves, an increased number of tillers, delayed heading, cluster-shaped shortening of panicles and mostly unfilled grains. Slight differences in symptom occurrence time were observed under different environmental conditions. For example, mosaic symptoms appeared earlier and crinkling symptoms appeared later, with both symptoms later receding in some infected plants. Yield losses due to RSMV also differed among varieties. The most serious yield reduction was experienced by indica rice (cv. Meixiangzhan), followed by hybrid indica rice (cv. Wuyou 1179) and then japonica (cv. Nipponbare). Single panicle weight, seed setting rate and 1000-kernel weight were reduced in the three infected varieties compared with healthy plants—by 85.42, 94.85 and 31.56% in Meixiangzhan; 52.43, 53.06 and 25.65% in Wuyou 1179 and 25.53, 49.32 and 23.86% in Nipponbare, respectively.

Conclusions

Our findings contribute basic data for field investigations, formulation of prevention and control strategies and further study of the pathogenesis of RSMV.

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.

Wood Transcriptome Profiling Identifies Critical Pathway Genes of Secondary Wall Biosynthesis and Novel Regulators for Vascular Cambium Development in Populus

Wood, the most abundant biomass on Earth, is composed of secondary xylem differentiated from vascular cambium. However, the underlying molecular mechanisms of wood formation remain largely unclear. To gain insight into wood formation, we performed a series of wood-forming tissue-specific transcriptome analyses from a hybrid poplar (Populus alba × P. glandulosa, clone BH) using RNA-seq. Together with shoot apex and leaf tissue, cambium and xylem tissues were isolated from vertical stem segments representing a gradient of secondary growth developmental stages (i.e., immature, intermediate, and mature stem). In a comparative transcriptome analysis of the 'developing xylem' and 'leaf' tissue, we could identify critical players catalyzing each biosynthetic step of secondary wall components (e.g., cellulose, xylan, and lignin). Several candidate genes involved in the initiation of vascular cambium formation were found via a co-expression network analysis using abundantly expressed genes in the 'intermediate stem-derived cambium' tissue. We found that transgenic Arabidopsis plants overexpressing the PtrHAM4-1, a GRAS family transcription factor, resulted in a significant increase of vascular cambium development. This phenotype was successfully reproduced in the transgenic poplars overexpressing the PtrHAM4-1. Taken together, our results may serve as a springboard for further research to unravel the molecular mechanism of wood formation, one of the most important biological processes on this planet.

Co-expression networks for plant biology: why and how

Co-expression network analysis is one of the most powerful approaches for interpretation of large transcriptomic datasets. It enables characterization of modules of co-expressed genes that may share biological functional linkages. Such networks provide an initial way to explore functional associations from gene expression profiling and can be applied to various aspects of plant biology. This review presents the applications of co-expression network analysis in plant biology and addresses optimized strategies from the recent literature for performing co-expression analysis on plant biological systems. Additionally, we describe the combined interpretation of co-expression analysis with other genomic data to enhance the generation of biologically relevant information.

A systems biology view of wood formation in Eucalyptus grandis trees submitted to different potassium and water regimes

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

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

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

Genome-Wide Characterization of AspATs in Populus: Gene Expression Variation and Enzyme Activities in Response to Nitrogen Perturbations

Aspartate aminotransferase (AspAT) catalyzes a reversible transamination reaction between glutamate and oxaloacetate to yield aspartate and 2-oxoglutarate, exerting a primary role in amino acid biosynthesis and homeostasis of nitrogen (N) and carbon metabolism within all cellular organisms. While progress in biochemical characterization of AspAT has been made for decades, the molecular and physiological characteristics of different members of the AspAT gene family remain poorly known particularly in forest trees. Here, extensive genome-wide survey of AspAT encoding genes was implemented in black cottonwood (Populus trichocarpa Torr. & A. Gray), a model species of woody plants. Thorough inspection of the phylogenies, gene structures, chromosomal distribution, cis-elements, conserved motifs, and subcellular targeting resulted in the identification of 10 AspAT isogenes (PtAspAT1-10) in the Populus genome. RNA-seq along with quantitative real-time polymerase chain reaction (qRT-PCR) validation revealed that PtAspATs displayed diverse patterns of tissue-specific expression. Spatiotemporal expressions of homologous AspATs in the poplar hybrid clone ‘Nanlin895’ were further evaluated, showing that gene expressions varied depending on source-sink dynamics. The impact on AspAT transcripts upon N starvation and seasonal senescence showed the upregulation of five AspAT in leaves concurrent with drastic downregulation of six or more AspATs in roots. Additionally, marked reductions of many more AspATs transcripts were observed in roots upon N excess. Accordingly, AspAT activities were significantly suppressed upon N starvation by an in-gel assay, prompting the argument that enzyme activity was a more direct indicator of the growth morphology under a N stress regime. Taken together, the expression profiling and enzyme activities upon stress cues provide a theoretical basis for unraveling the physiological significance of specific gene(s) in regulation of N acquisition and remobilization in woody plants.

Genome-wide transcriptome analysis reveals the molecular mechanism of high temperature-induced floral abortion in Litchi chinensis

Background

Warm winter and hot spring attributed to global warming affected floral development and may induce floral abortion, resulted in poor flowering in litchi. To identify genes potentially involved in litchi floral abortion, six RNA-sequencing (RNA-Seq) libraries of the developing panicles (DPs) under low temperature (LT) conditions and the shrinking panicles (SPs) under high temperature (HT) conditions were constructed.

Results

3.07–8.97 × 10⁶ clean reads were generated. Digital expression of the DPs with that of the SPs was compared. As a result, 1320 up-regulated and 981 down-regulated differentially expressed genes (DEGs) were identified. From the enriched GO-term, 54 temperature responsive DEGs, 23 hormone homeostasis- or biosynthesis-related DEGs, 137 hormone signal transduction or responsive DEGs, and 18 flowering-related DEGs were identified. Partial Least Squares Structural Equation Modeling (PLS-SEM) analysis indicated that the effects of hormone-related DEGs on NACs, MYBs, WRKYs were stronger than that on flowering-related DEGs. Expression pattern analysis of the inflorescence in ‘Nuomici’ and ‘Huaizhi’ under LT and HT conditions showed that genes homologous to AIL6 (LcAIL6), LHY (LcLHY), MED16 (LcMED16), SKIP20 (LcSKIP20), POD20 (LcPOD20) in the two cultivars had similar expression trends.

Conclusion

This study elucidated the transcriptome in the HT-induced floral abortion and identified key genes involved in the process. NACs, MYBs, WRKYs may act as central players involved in the HT-induced floral abortion underlying hormonal control. Increased transcript in LcLHY, LcMED16, LcSKIP20, LcPOD20 and decreased transcript in LcAIL6 might be related to the inhibition of floral development. Our studies provided potential genes for the future molecular breeding of new cultivars that can reduce floral abortion under warm climates, and a novel clue to reveal the relationship of biological processes based on the RNA-seq data using PLS-SEM.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5493-8) contains supplementary material, which is available to authorized users.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

A New Calmodulin-Binding Protein Expresses in the Context of Secondary Cell Wall Biosynthesis and Impacts Biomass Properties in Populus

A greater understanding of biosynthesis, signaling and regulatory pathways involved in determining stem growth and secondary cell wall chemistry is important for enabling pathway engineering and genetic optimization of biomass properties. The present study describes a new functional role of PdIQD10, a Populus gene belonging to the IQ67-Domain1 family of IQD genes, in impacting biomass formation and chemistry. Expression studies showed that PdIQD10 has enhanced expression in developing xylem and tension-stressed tissues in Populus deltoides. Molecular dynamics simulation and yeast two-hybrid interaction experiments suggest interactions with two calmodulin proteins, CaM247 and CaM014, supporting the sequence-predicted functional role of the PdIQD10 as a calmodulin-binding protein. PdIQD10 was found to interact with specific Populus isoforms of the Kinesin Light Chain protein family, shown previously to function as microtubule-guided, cargo binding and delivery proteins in Arabidopsis. Subcellular localization studies showed that PdIQD10 localizes in the nucleus and plasma membrane regions. Promoter-binding assays suggest that a known master transcriptional regulator of secondary cell wall biosynthesis (PdWND1B) may be upstream of an HD-ZIP III gene that is in turn upstream of PdIQD10 gene in the transcriptional network. RNAi-mediated downregulation of PdIQD10 expression resulted in plants with altered biomass properties including higher cellulose, wall glucose content and greater biomass quantity. These results present evidence in support of a new functional role for an IQD gene family member, PdIQD10, in secondary cell wall biosynthesis and biomass formation in Populus.

The regulatory role of shikimate in plant phenylalanine metabolism

In higher plants, the amino acid phenylalanine is a substrate of both primary and secondary metabolic pathways. The primary pathway that consumes phenylalanine, protein biosynthesis, is essential for the viability of all cells. Meanwhile, the secondary pathways are not necessary for the survival of individual cells, but benefit of the plant as a whole. Here we focus on the monolignol pathway, a secondary metabolic pathway in the cytosol that rapidly consumes phenylalanine to produce the precursors of lignin during wood formation. In planta monolignol biosynthesis involves a series of seemingly redundant steps wherein shikimate, a precursor of phenylalanine synthesized in the plastid, is transiently ligated to the main substrate of the pathway. However, shikimate is not catalytically involved in the reactions of the monolignol pathway, and is only needed for pathway enzymes to recognize their main substrates. After some steps the shikimate moiety is removed unaltered, and the main substrate continues along the pathway. It has been suggested that this portion of the monolignol pathway fulfills a regulatory role in the following way. Low phenylalanine concentrations (viz. availability) correlate with low shikimate concentrations. When shikimate concentratios are low, flux into the monolignol pathway will be limited by means of the steps requiring shikimate. Thus, when the concentration of phenylalanine is low it will be reserved for protein biosynthesis. Here we employ a theoretical approach to test this hypothesis. Simplified versions of plant phenylalanine metabolism are modelled as systems of ordinary differential equations. Our analysis shows that the seemingly redundant steps can be sufficient for the prioritization of protein biosynthesis over the monolignol pathway when the availability of phenylalanine is low, depending on system parameters. Thus, the phenylalanine precursor shikimate may signal low phenylalanine availability to secondary pathways. Because our models have been abstracted from plant phenylalanine metabolism, this mechanism of metabolic signalling, which we call the Precursor Shutoff Valve (PSV), may also be present in other biochemical networks comprised of two pathways that share a common substrate.

Global transcriptome and weighted gene co-expression network analyses reveal hybrid-specific modules and candidate genes related to plant height development in maize

Weighted gene co-expression network analysis was explored to find key hub genes involved in plant height regulation. Plant height, an important trait for maize breeding because of its close relatedness to lodging resistance and yield, has been reported to be determined by multiple qualitative and quantitative genes. However, few genes related to plant height have been characterized in maize. Herein, three different maize hybrids, with extremely distinct plant height, which were further classified into low (L), middle (M) and high (H) group, were selected for RNA sequencing at three key developmental stages, namely, jointing stage (I), big flare period (II) and tasseling stage (III). Intriguingly, transcriptome profiles for hybrids ranging from low to high group exhibited significantly similarity in both jointing stage and big flare period. However, remarkably larger differentially expressed genes could be detected between hybrid from low to either middle or high group in tasseling stage. These results were repeatedly observed in both phenotyping and gene ontology enrichment analysis, indicating that transition from big flare period to tasseling stage plays a critical role in determination of plant height. Furthermore, weighted gene co-expression network analysis was explored to find key hub genes involved in plant height regulation. Hundreds of candidate genes, encoding various transcription factors, and regulators involved in internode cell regulation and cell wall synthesis were identified in our network. More importantly, great majority of candidates were correlated to either metabolism or signaling pathway of several plant phytohormones. Particularly, numerous functionally characterized genes in gibberellic acid as well as brassinosteroids signaling transduction pathways were also discovered, suggesting their critical roles in plant height regulation. The present study could provide a modestly comprehensive insight into networks for regulation of plant height in maize.

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.

Ethylene signaling induces gelatinous layers with typical features of tension wood in hybrid aspen

The phytohormone ethylene impacts secondary stem growth in plants by stimulating cambial activity, xylem development and fiber over vessel formation. We report the effect of ethylene on secondary cell wall formation and the molecular connection between ethylene signaling and wood formation. We applied exogenous ethylene or its precursor 1-aminocyclopropane-1-carboxylic acid (ACC) to wild-type and ethylene-insensitive hybrid aspen trees (Populus tremula × tremuloides) and studied secondary cell wall anatomy, chemistry and ultrastructure. We furthermore analyzed the transcriptome (RNA Seq) after ACC application to wild-type and ethylene-insensitive trees. We demonstrate that ACC and ethylene induce gelatinous layers (G-layers) and alter the fiber cell wall cellulose microfibril angle. G-layers are tertiary wall layers rich in cellulose, typically found in tension wood of aspen trees. A vast majority of transcripts affected by ACC are downstream of ethylene perception and include a large number of transcription factors (TFs). Motif-analyses reveal potential connections between ethylene TFs (Ethylene Response Factors (ERFs), ETHYLENE INSENSITIVE 3/ETHYLENE INSENSITIVE3-LIKE1 (EIN3/EIL1)) and wood formation. G-layer formation upon ethylene application suggests that the increase in ethylene biosynthesis observed during tension wood formation is important for its formation. Ethylene-regulated TFs of the ERF and EIN3/EIL1 type could transmit the ethylene signal.

Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis

A multi-omics quantitative integrative analysis of lignin biosynthesis can advance the strategic engineering of wood for timber, pulp, and biofuels. Lignin is polymerized from three monomers (monolignols) produced by a grid-like pathway. The pathway in wood formation of Populus trichocarpa has at least 21 genes, encoding enzymes that mediate 37 reactions on 24 metabolites, leading to lignin and affecting wood properties. We perturb these 21 pathway genes and integrate transcriptomic, proteomic, fluxomic and phenomic data from 221 lines selected from ~2000 transgenics (6-month-old). The integrative analysis estimates how changing expression of pathway gene or gene combination affects protein abundance, metabolic-flux, metabolite concentrations, and 25 wood traits, including lignin, tree-growth, density, strength, and saccharification. The analysis then predicts improvements in any of these 25 traits individually or in combinations, through engineering expression of specific monolignol genes. The analysis may lead to greater understanding of other pathways for improved growth and adaptation.

A systematic analysis of lignin biosynthetic genes to quantitatively understand their effect on wood properties is still lacking. Here, the authors integrate transcriptomic, proteomic, fluxomic and phenomic data to quantify the impact of perturbations of transcript abundance on lignin biosynthesis and wood properties.

Identification of new protein-protein and protein-DNA interactions linked with wood formation in Populus trichocarpa

Cellular processes, such as signal transduction and cell wall deposition, are organized by macromolecule interactions. Experimentally determined protein-protein interactions (PPIs) and protein-DNA interactions (PDIs) relevant to woody plant development are sparse. To begin to develop a Populus trichocarpa Torr. & A. Gray wood interactome, we applied the yeast-two-hybrid (Y2H) assay in different ways to enable the discovery of novel PPIs and connected networks. We first cloned open reading frames (ORFs) for 361 genes markedly upregulated in secondary xylem compared with secondary phloem and performed a binary Y2H screen with these proteins. By screening a xylem cDNA library for interactors of a subset of these proteins and then recapitulating the process by using a subset of the interactors as baits, we ultimately identified 165 PPIs involving 162 different ORFs. Thirty-eight transcription factors (TFs) included in our collection of P. trichocarpa wood ORFs were used in a Y1H screen for binding to promoter regions of three genes involved in lignin biosynthesis resulting in 40 PDIs involving 20 different TFs. The network incorporating both the PPIs and PDIs included 14 connected subnetworks, with the largest having 132 members. Protein-protein interactions and PDIs validated previous reports and also identified new candidate wood formation proteins and modules through their interactions with proteins and promoters known to be involved in secondary cell wall synthesis. Selected examples are discussed including a PPI between Mps one binder (MOB1) and a mitogen-activated protein kinase kinase kinase kinase (M4K) that was further characterized by assays confirming the PPI as well as its effect on subcellular localization. Mapping of published transcriptomic data showing developmentally detailed expression patterns across a secondary stem onto the network supported that the PPIs and PDIs are relevant to wood formation, and also illustrated that wood-associated interactions involve gene products that are not upregulated in secondary xylem.

Enzyme-Enzyme Interactions in Monolignol Biosynthesis

The enzymes that comprise the monolignol biosynthetic pathway have been studied intensively for more than half a century. A major interest has been the role of pathway in the biosynthesis of lignin and the role of lignin in the formation of wood. The pathway has been typically conceived as linear steps that convert phenylalanine into three major monolignols or as a network of enzymes in a metabolic grid. Potential interactions of enzymes have been investigated to test models of metabolic channeling or for higher order interactions. Evidence for enzymatic or physical interactions has been fragmentary and limited to a few enzymes studied in different species. Only recently the entire pathway has been studied comprehensively in any single plant species. Support for interactions comes from new studies of enzyme activity, co-immunoprecipitation, chemical crosslinking, bimolecular fluorescence complementation, yeast 2-hybrid functional screening, and cell type-specific gene expression based on light amplification by stimulated emission of radiation capture microdissection. The most extensive experiments have been done on differentiating xylem of Populus trichocarpa, where genomic, biochemical, chemical, and cellular experiments have been carried out. Interactions affect the rate, direction, and specificity of both 3 and 4-hydroxylation in the monolignol biosynthetic pathway. Three monolignol P450 mono-oxygenases form heterodimeric and heterotetrameric protein complexes that activate specific hydroxylation of cinnamic acid derivatives. Other interactions include regulatory kinetic control of 4-coumarate CoA ligases through subunit specificity and interactions between a cinnamyl alcohol dehydrogenase and a cinnamoyl-CoA reductase. Monolignol enzyme interactions with other pathway proteins have been associated with biotic and abiotic stress response. Evidence challenging or supporting metabolic channeling in this pathway will be discussed.

MorphDB: Prioritizing Genes for Specialized Metabolism Pathways and Gene Ontology Categories in Plants

Recent times have seen an enormous growth of "omics" data, of which high-throughput gene expression data are arguably the most important from a functional perspective. Despite huge improvements in computational techniques for the functional classification of gene sequences, common similarity-based methods often fall short of providing full and reliable functional information. Recently, the combination of comparative genomics with approaches in functional genomics has received considerable interest for gene function analysis, leveraging both gene expression based guilt-by-association methods and annotation efforts in closely related model organisms. Besides the identification of missing genes in pathways, these methods also typically enable the discovery of biological regulators (i.e., transcription factors or signaling genes). A previously built guilt-by-association method is MORPH, which was proven to be an efficient algorithm that performs particularly well in identifying and prioritizing missing genes in plant metabolic pathways. Here, we present MorphDB, a resource where MORPH-based candidate genes for large-scale functional annotations (Gene Ontology, MapMan bins) are integrated across multiple plant species. Besides a gene centric query utility, we present a comparative network approach that enables researchers to efficiently browse MORPH predictions across functional gene sets and species, facilitating efficient gene discovery and candidate gene prioritization. MorphDB is available at http://bioinformatics.psb.ugent.be/webtools/morphdb/morphDB/index/. We also provide a toolkit, named "MORPH bulk" (https://github.com/arzwa/morph-bulk), for running MORPH in bulk mode on novel data sets, enabling researchers to apply MORPH to their own species of interest.

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

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