Identification and biochemical characterization of four wood-associated glucuronoxylan methyltransferases in Populus

Flexure wood formation via growth reprogramming in hybrid aspen involves jasmonates and polyamines and transcriptional changes resembling tension wood development

Stem bending in trees induces flexure wood but its properties and development are poorly understood. Here, we investigated the effects of low-intensity multidirectional stem flexing on growth and wood properties of hybrid aspen, and on its transcriptomic and hormonal responses. Glasshouse-grown trees were either kept stationary or subjected to several daily shakes for 5 wk, after which the transcriptomes and hormones were analyzed in the cambial region and developing wood tissues, and the wood properties were analyzed by physical, chemical and microscopy techniques. Shaking increased primary and secondary growth and altered wood differentiation by stimulating gelatinous-fiber formation, reducing secondary wall thickness, changing matrix polysaccharides and increasing cellulose, G- and H-lignin contents, cell wall porosity and saccharification yields. Wood-forming tissues exhibited elevated jasmonate, polyamine, ethylene and brassinosteroids and reduced abscisic acid and gibberellin signaling. Transcriptional responses resembled those during tension wood formation but not opposite wood formation and revealed several thigmomorphogenesis-related genes as well as novel gene networks including FLA and XTH genes encoding plasma membrane-bound proteins. Low-intensity stem flexing stimulates growth and induces wood having improved biorefinery properties through molecular and hormonal pathways similar to thigmomorphogenesis in herbaceous plants and largely overlapping with the tension wood program of hardwoods.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

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

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

Analytical Techniques for Determining the Role of Domain of Unknown Function 579 Proteins in the Synthesis of O-Methylated Plant Polysaccharides

Matrix polysaccharides are a diverse group of structurally complex carbohydrates and account for a large portion of the biomass consumed as food or used to produce fuels and materials. Glucuronoxylan and arabinogalactan protein are matrix glycans that have sidechains decorated with 4-O-methyl glucuronosyl residues. Methylation is a key determinant of the physical properties of these wall glycopolymers and consequently affects both their biological function and ability to interact with other wall polymers. Indeed, there is increasing interest in determining the distribution and abundance of methyl-etherified polysaccharides in different plant species, tissues, and developmental stages. There is also a need to understand the mechanisms involved in their biosynthesis. Members of the Domain of Unknown Function (DUF) 579 family have been demonstrated to have a role in the biosynthesis of methyl-etherified glycans. Here we describe methods for the analysis of the 4-O-methyl glucuronic acid moieties that are present in sidechains of arabinogalactan proteins. These methods are then applied toward the analysis of loss-of-function mutants of two DUF579 family members that lack this modification in muro. We also present a procedure to assay DUF579 family members for enzymatic activity in vitro using acceptor oligosaccharides prepared from xylan of loss-of-function mutants. Our approach facilitates the characterization of enzymes that modify glycosyl residues during cell wall synthesis and the structures that they generate.

A systems genetics analysis in Eucalyptus reveals coordination of metabolic pathways associated with xylan modification in wood-forming tissues

Acetyl- and methylglucuronic acid decorations of xylan, the dominant hemicellulose in secondary cell walls (SCWs) of woody dicots, affect its interaction with cellulose and lignin to determine SCW structure and extractability. Genes and pathways involved in these modifications may be targets for genetic engineering; however, little is known about the regulation of xylan modifications in woody plants. To address this, we assessed genetic and gene expression variation associated with xylan modification in developing xylem of Eucalyptus grandis × Eucalyptus urophylla interspecific hybrids. Expression quantitative trait locus (eQTL) mapping identified potential regulatory polymorphisms affecting gene expression modules associated with xylan modification. We identified 14 putative xylan modification genes that are members of five expression modules sharing seven trans-eQTL hotspots. The xylan modification genes are prevalent in two expression modules. The first comprises nucleotide sugar interconversion pathways supplying the essential precursors for cellulose and xylan biosynthesis. The second contains genes responsible for phenylalanine biosynthesis and S-adenosylmethionine biosynthesis required for glucuronic acid and monolignol methylation. Co-expression and co-regulation analyses also identified four metabolic sources of acetyl coenxyme A that appear to be transcriptionally coordinated with xylan modification. Our systems genetics analysis may provide new avenues for metabolic engineering to alter wood SCW biology for enhanced biomass processability.

Poplar carbohydrate-active enzymes: whole-genome annotation and functional analyses based on RNA expression data

Carbohydrate-active enzymes (CAZymes) catalyze the formation and modification of glycoproteins, glycolipids, starch, secondary metabolites and cell wall biopolymers. They are key enzymes for the biosynthesis of food and renewable biomass. Woody biomass is particularly important for long-term carbon storage and as an abundant renewable natural resource for many industrial applications. This study presents a re-annotation of CAZyme genes in the current Populus trichocarpa genome assembly and in silico functional characterization, based on high-resolution RNA-Seq data sets. Altogether, 1914 CAZyme and expansin genes were annotated in 101 families. About 1797 of these genes were found expressed in at least one Populus organ. We identified genes involved in the biosynthesis of different cell wall polymers and their paralogs. Whereas similar families exist in poplar and Arabidopsis thaliana (with the exception of CBM13 found only in poplar), a few families had significantly different copy numbers between the two species. To identify the transcriptional coordination and functional relatedness within the CAZymes and other proteins, we performed co-expression network analysis of CAZymes in wood-forming tissues using the AspWood database (http://aspwood.popgenie.org/aspwood-v3.0/) for Populus tremula. This provided an overview of the transcriptional changes in CAZymes during the transition from primary to secondary wall formation, and the clustering of transcripts into potential regulons. Candidate enzymes involved in the biosynthesis of polysaccharides were identified along with many tissue-specific uncharacterized genes and transcription factors. These collections offer a rich source of targets for the modification of secondary cell wall biosynthesis and other developmental processes in woody plants.

Secondary cell wall biosynthesis

Contents Summary 1703 I. Introduction 1703 II. Cellulose biosynthesis 1705 III. Xylan biosynthesis 1709 IV. Glucomannan biosynthesis 1713 V. Lignin biosynthesis 1714 VI. Concluding remarks 1717 Acknowledgements 1717 References 1717 SUMMARY: Secondary walls are synthesized in specialized cells, such as tracheary elements and fibers, and their remarkable strength and rigidity provide strong mechanical support to the cells and the plant body. The main components of secondary walls are cellulose, xylan, glucomannan and lignin. Biochemical, molecular and genetic studies have led to the discovery of most of the genes involved in the biosynthesis of secondary wall components. Cellulose is synthesized by cellulose synthase complexes in the plasma membrane and the recent success of in vitro synthesis of cellulose microfibrils by a single recombinant cellulose synthase isoform reconstituted into proteoliposomes opens new doors to further investigate the structure and functions of cellulose synthase complexes. Most genes involved in the glycosyl backbone synthesis, glycosyl substitutions and acetylation of xylan and glucomannan have been genetically characterized and the biochemical properties of some of their encoded enzymes have been investigated. The genes and their encoded enzymes participating in monolignol biosynthesis and modification have been extensively studied both genetically and biochemically. A full understanding of how secondary wall components are synthesized will ultimately enable us to produce plants with custom-designed secondary wall composition tailored to diverse applications.

Analysis of carbohydrates and glycoconjugates by matrix-assisted laser desorption/ionization mass spectrometry: An update for 2013-2014

This review is the eighth update of the original article published in 1999 on the application of Matrix-assisted laser desorption/ionization mass spectrometry (MALDI) mass spectrometry to the analysis of carbohydrates and glycoconjugates and brings coverage of the literature to the end of 2014. Topics covered in the first part of the review include general aspects such as theory of the MALDI process, matrices, derivatization, MALDI imaging, fragmentation, and arrays. The second part of the review is devoted to applications to various structural types such as oligo- and poly- saccharides, glycoproteins, glycolipids, glycosides, and biopharmaceuticals. Much of this material is presented in tabular form. The third part of the review covers medical and industrial applications of the technique, studies of enzyme reactions, and applications to chemical synthesis. © 2018 Wiley Periodicals, Inc. Mass Spec Rev 37:353-491, 2018.

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

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

A Genome Wide Association Study of arabinoxylan content in 2-row spring barley grain

In barley endosperm arabinoxylan (AX) is the second most abundant cell wall polysaccharide and in wheat it is the most abundant polysaccharide in the starchy endosperm walls of the grain. AX is one of the main contributors to grain dietary fibre content providing several health benefits including cholesterol and glucose lowering effects, and antioxidant activities. Due to its complex structural features, AX might also affect the downstream applications of barley grain in malting and brewing. Using a high pressure liquid chromatography (HPLC) method we quantified AX amounts in mature grain in 128 spring 2-row barley accessions. Amounts ranged from ~ 5.2 μg/g to ~ 9 μg/g. We used this data for a Genome Wide Association Study (GWAS) that revealed three significant quantitative trait loci (QTL) associated with grain AX levels which passed a false discovery threshold (FDR) and are located on two of the seven barley chromosomes. Regions underlying the QTLs were scanned for genes likely to be involved in AX biosynthesis or turnover, and strong candidates, including glycosyltransferases from the GT43 and GT61 families and glycoside hydrolases from the GH10 family, were identified. Phylogenetic trees of selected gene families were built based on protein translations and were used to examine the relationship of the barley candidate genes to those in other species. Our data reaffirms the roles of existing genes thought to contribute to AX content, and identifies novel QTL (and candidate genes associated with them) potentially influencing the AX content of barley grain. One potential outcome of this work is the deployment of highly associated single nucleotide polymorphisms markers in breeding programs to guide the modification of AX abundance in barley grain.

Cytosol-Localized UDP-Xylose Synthases Provide the Major Source of UDP-Xylose for the Biosynthesis of Xylan and Xyloglucan

Xylan and xyloglucan are the two major cell wall hemicelluloses in plants, and their biosynthesis requires a steady supply of the sugar donor, UDP-xylose. UDP-xylose is synthesized through conversion of UDP-glucuronic acid (UDP-GlcA) by the activities of UDP-xylose synthase (UXS). There exist six UXS genes in the Arabidopsis thaliana genome; three of them (UXS1, UXS2 and UXS4) encode Golgi-localized enzymes and the other three (UXS3, UXS5 and UXS6) encode cytosol-localized enzymes. In this report, we investigated the contributions of these UXS genes in supplying UDP-xylose for the biosynthesis of xylan and xyloglucan. Expression analyses revealed that the six UXS genes exhibited distinct and overlapping expression patterns in different cell types of stems, root-hypocotyls and young seedlings, and that the relative enzymatic activity of UXS in the cytosol was 17 times higher than that in the Golgi. Among the six UXS genes, UXS3, UXS5 and UXS6 showed the highest expression in stems and were expressed predominantly in xylem cells and interfascicular fibers. Their predominant expression in secondary wall-forming cells was consistent with the finding that the expression of UXS3, UXS5 and UXS6 was directly activated by the secondary wall NAC master switches. Although simultaneous mutations of UXS1, UXS2 and UXS4 did not cause any apparent effects on plant growth and xylan biosynthesis, simultaneous down-regulation/mutations of UXS3, UXS5 and UXS6 led to a drastic reduction in secondary wall thickening, a severe deformation of xylem vessels, a significant decrease in xylan content without an apparent reduction in its chain length and an absence of GlcA side chains in xylan, which are reminiscent of the phenotypes of some known xylan-deficient mutants. Moreover, Immunolocalization with two xyloglucan monoclonal antibodies, LM15 and LM25, revealed a significant reduction in the amount of xylogulcan in the primary walls. These results demonstrate that the cytosol-localized UXS3, UXS5 and UXS6 play a predominant role in the supply of UDP-xylose for the biosynthesis of xylan and xyloglucan.

Evolutionary Conservation of Xylan Biosynthetic Genes in Selaginella moellendorffii and Physcomitrella patens

Xylan is a major cross-linking hemicellulose in secondary walls of vascular tissues, and the recruitment of xylan as a secondary wall component was suggested to be a pivotal event for the evolution of vascular tissues. To decipher the evolution of xylan structure and xylan biosynthetic genes, we analyzed xylan substitution patterns and characterized genes mediating methylation of glucuronic acid (GlcA) side chains in xylan of the model seedless vascular plant, Selaginella moellendorffii, and investigated GT43 genes from S. moellendorffii and the model non-vascular plant, Physcomitrella patens, for their roles in xylan biosynthesis. Using nuclear magentic resonance spectroscopy, we have demonstrated that S. moellendorffii xylan consists of β-1,4-linked xylosyl residues subsituted solely with methylated GlcA residues and that xylans from both S. moellendorffii and P. patens are acetylated at O-2 and O-3. To investigate genes responsible for GlcA methylation of xylan, we identified two DUF579 genes in the S. moellendorffii genome and showed that one of them, SmGXM, encodes a glucuronoxylan methyltransferase capable of adding the methyl group onto the GlcA side chain of xylooligomers. Furthermore, we revealed that the two GT43 genes in S. moellendorffii, SmGT43A and SmGT43B, are functional orthologs of the Arabidopsis xylan backbone biosynthetic genes IRX9 and IRX14, respectively, indicating the evolutionary conservation of the involvement of two functionally non-redundant groups of GT43 genes in xylan backbone biosynthesis between seedless and seed vascular plants. Among the five GT43 genes in P. patens, PpGT43A was found to be a functional ortholog of Arabidopsis IRX9, suggesting that the recruitment of GT43 genes in xylan backbone biosynthesis occurred when non-vascular plants appeared on land.

Suppression of PtrDUF579-3 Expression Causes Structural Changes of the Glucuronoxylan in Populus

DUF579 (domain unknown function 579) genes have been reported to play diverse roles in cell wall biosynthesis, such as in glucuronoxylan (GX) synthesis. As GX is a major type of hemicelluloses in hard wood species, how DUF579 genes function in wood formation remains to be demonstrated in planta. This study reports a Populus DUF579 gene, PtrDUF579-3, which is characterized for its function in wood cell wall formation. PtrDUF579-3 is localized in Golgi apparatus and catalyzes methylation of the glucuronic acid (GlcA) in GX biosynthesis. Suppression of PtrDUF579-3 expression in Populus caused a reduction in both the GlcA side chain number and GlcA side chain methylation on the GX backbone. The modified GX polymer through PtrDUF579-3 suppression was more susceptible to acid treatment and the PtrDUF579-3 suppressed plants displayed enhanced cellulose digestibility. These results suggest that PtrDUF579-3 is involved in GX biosynthesis and GX structure can be modified through PtrDUF579-3 suppression in Populus.

Molecular control of wood formation in trees

Wood (also termed secondary xylem) is the most abundant biomass produced by plants, and is one of the most important sinks for atmospheric carbon dioxide. The development of wood begins with the differentiation of the lateral meristem, vascular cambium, into secondary xylem mother cells followed by cell expansion, secondary wall deposition, programmed cell death, and finally heartwood formation. Significant progress has been made in the past decade in uncovering the molecular players involved in various developmental stages of wood formation in tree species. Hormonal signalling has been shown to play critical roles in vascular cambium cell proliferation and a peptide-receptor-transcription factor regulatory mechanism similar to that controlling the activity of apical meristems is proposed to be involved in the maintenance of vascular cambium activity. It has been demonstrated that the differentiation of vascular cambium into xylem mother cells is regulated by plant hormones and HD-ZIP III transcription factors, and the coordinated activation of secondary wall biosynthesis genes during wood formation is mediated by a transcription network encompassing secondary wall NAC and MYB master switches and their downstream transcription factors. Most genes encoding the biosynthesis enzymes for wood components (cellulose, xylan, glucomannan, and lignin) have been identified in poplar and a number of them have been functionally characterized. With the availability of genome sequences of tree species from both gymnosperms and angiosperms, and the identification of a suite of wood-associated genes, it is expected that our understanding of the molecular control of wood formation in trees will be greatly accelerated.

Diverse roles of PtrDUF579 proteins in Populus and PtrDUF579-1 function in vascular cambium proliferation during secondary growth

DUF579 (domain of unknown function 579) family proteins contain a DUF579 domain structure but vary greatly in their overall sequence similarity. Several DUF579 proteins have been found to play a role in cell wall biosynthesis in Arabidopsis, while DUF579 family genes have not yet been systematically investigated in Populus. In this study, the Populus DUF579 family proteins were found to be localized in different cell types and subcellular locations. The diverse expression patterns of the proteins indicate that they may perform different functions in Populus. Among the DUF579 family members, PtrDUF579-1 is found to be specifically expressed in vascular cambium zone cells where it is localized in the Golgi apparatus. Suppression of PtrDUF579-1 expression reduced plant height and stem diameter size. Cambium cell division and xylem tissue growth was inhibited while secondary cell wall formation was unchanged in PtrDUF579-1 suppressed plants. Cell walls analysis showed that the composition of the pectin fraction of the cambium cell wall was altered while other polysaccharides were not affected in PtrDUF579-1 suppressed plants. This observation suggest cambium expressed PtrDUF579-1 may affect cell wall biosynthesis and be involved in cambium cell proliferation in Populus. Overall, DUF579 family proteins play a diverse set of roles in Populus.