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

Systems genetic analysis of lignin biosynthesis in Populus tremula

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

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.

Whole-genome sequencing and antimicrobial potential of bacteria isolated from Polish honey

The aim of this study was the whole-genome analysis and assessment of the antimicrobial potential of bacterial isolates from honey harvested in one geographical location-the north of Poland. In total, 132 strains were derived from three honey samples, and the antimicrobial activity of CFAM (cell-free after-culture medium) was used as a criterion for strain selection and detailed genomic investigation. Two of the tested isolates (SZA14 and SZA16) were classified as Bacillus paralicheniformis, and one isolate (SZB3) as Bacillus subtilis based on their ANI and phylogenetic analysis relatedness. The isolates SZA14 and SZA16 were harvested from the same honey sample with a nucleotide identity of 98.96%. All three isolates have been found to be potential producers of different antimicrobial compounds. The secondary metabolite genome mining pipeline (antiSMASH) identified 14 gene cluster coding for non-ribosomal peptide synthetases (NRPs), polyketide synthases (PKSs), and ribosomally synthesized and post-translationally modified peptides (RiPPs) that are potential sources of novel antibacterials. The BAGEL4 analysis revealed the presence of nine putative gene clusters of interest in the isolates SZA14 and SZA16 (including the presence of six similar clusters present in both isolates, coding for the production of enterocin Nkr-5-3B, haloduracin-alpha, sonorensin, bottromycin, comX2, and lasso peptide), and four in B. subtilis isolate SZB3 (competence factor, sporulation-killing factor, subtilosin A, and sactipeptides). The outcomes of this study confirm that honey-derived Bacillus spp. strains can be considered potential producers of a broad spectrum of antimicrobial agents. KEY POINTS: • Bacteria of the genus Bacillus are an important component of honey microbiota. • Honey-derived Bacillus spp. strains are potential producers of new antimicrobials.

Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry

Production of biofuel from lignocellulosic biomass is relatively low due to the limited knowledge about natural cell wall loosening and cellulolytic processes in plants. Industrial separation of cellulose fiber mass from lignin, its saccharification and alcoholic fermentation is still cost-ineffective and environmentally unfriendly. Assuming that the green transformation is inevitable and that new sources of raw materials for biofuels are needed, we decided to study cell death-a natural process occurring in plants in the context of reducing the recalcitrance of lignocellulose for the production of second-generation bioethanol. "Members of the enzyme families responsible for lysigenous aerenchyma formation were identified during the root hypoxia stress in Arabidopsis thaliana cell death mutants. The cell death regulatory genes, LESION SIMULATING DISEASE 1 (LSD1), PHYTOALEXIN DEFICIENT 4 (PAD4) and ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) conditionally regulate the cell wall when suppressed in transgenic aspen. During four years of growth in the field, the following effects were observed: lignin content was reduced, the cellulose fiber polymerization degree increased and the growth itself was unaffected. The wood of transgenic trees was more efficient as a substrate for saccharification, alcoholic fermentation and bioethanol production. The presented results may trigger the development of novel biotechnologies in the biofuel industry.

Impact of xylan on field productivity and wood saccharification properties in aspen

Xylan that comprises roughly 25% of hardwood biomass is undesirable in biorefinery applications involving saccharification and fermentation. Efforts to reduce xylan levels have therefore been made in many species, usually resulting in improved saccharification. However, such modified plants have not yet been tested under field conditions. Here we evaluate the field performance of transgenic hybrid aspen lines with reduced xylan levels and assess their usefulness as short-rotation feedstocks for biorefineries. Three types of transgenic lines were tested in four-year field tests with RNAi constructs targeting either Populus GT43 clades B and C (GT43BC) corresponding to Arabidopsis clades IRX9 and IRX14, respectively, involved in xylan backbone biosynthesis, GATL1.1 corresponding to AtGALT1 involved in xylan reducing end sequence biosynthesis, or ASPR1 encoding an atypical aspartate protease. Their productivity, wood quality traits, and saccharification efficiency were analyzed. The only lines differing significantly from the wild type with respect to growth and biotic stress resistance were the ASPR1 lines, whose stems were roughly 10% shorter and narrower and leaves showed increased arthropod damage. GT43BC lines exhibited no growth advantage in the field despite their superior growth in greenhouse experiments. Wood from the ASPR1 and GT43BC lines had slightly reduced density due to thinner cell walls and, in the case of ASPR1, larger cell diameters. The xylan was less extractable by alkali but more hydrolysable by acid, had increased glucuronosylation, and its content was reduced in all three types of transgenic lines. The hemicellulose size distribution in the GALT1.1 and ASPR1 lines was skewed towards higher molecular mass compared to the wild type. These results provide experimental evidence that GATL1.1 functions in xylan biosynthesis and suggest that ASPR1 may regulate this process. In saccharification without pretreatment, lines of all three constructs provided 8-11% higher average glucose yields than wild-type plants. In saccharification with acid pretreatment, the GT43BC construct provided a 10% yield increase on average. The best transgenic lines of each construct are thus predicted to modestly outperform the wild type in terms of glucose yields per hectare. The field evaluation of transgenic xylan-reduced aspen represents an important step towards more productive feedstocks for biorefineries.

Xylan glucuronic acid side chains fix suberin-like aliphatic compounds to wood cell walls

Wood is the most important repository of assimilated carbon in the biosphere, in the form of large polymers (cellulose, hemicelluloses including glucuronoxylan, and lignin) that interactively form a composite, together with soluble extractives including phenolic and aliphatic compounds. Molecular interactions among these compounds are not fully understood. We have targeted the expression of a fungal α-glucuronidase to the wood cell wall of aspen (Populus tremula L. × tremuloides Michx.) and Arabidopsis (Arabidopsis thaliana (L.) Heynh), to decrease contents of the 4-O-methyl glucuronopyranose acid (mGlcA) substituent of xylan, to elucidate mGlcA's functions. The enzyme affected the content of aliphatic insoluble cell wall components having composition similar to suberin, which required mGlcA for binding to cell walls. Such suberin-like compounds have been previously identified in decayed wood, but here, we show their presence in healthy wood of both hardwood and softwood species. By contrast, γ-ester bonds between mGlcA and lignin were insensitive to cell wall-localized α-glucuronidase, supporting the intracellular formation of these bonds. These findings challenge the current view of the wood cell wall composition and reveal a novel function of mGlcA substituent of xylan in fastening of suberin-like compounds to cell wall. They also suggest an intracellular initiation of lignin-carbohydrate complex assembly.

The electron transport mechanism of downflow Leersia hexandra Swartz constructed wetland-microbial fuel cell when used to treat Cr(VI) and p-chlorophenol

Constructed wetland-microbial fuel cells are used to treat heavy metal and/or refractory organic wastewater. However, the electron transport mechanism of downflow Leersia hexandra constructed wetland-microbial fuel cells (DLCW-MFCs) is poorly understood when used to treat composite-polluted wastewater containing Cr(VI) and p-chlorophenol (4-CP) (C&P). In this study, metagenomics and in situ electrochemical techniques were used to investigate the electrochemical properties and the electricigens and their dominant gene functions. The DLCW-MFC was used to treat C&P and single-pollutant wastewater containing Cr(VI) (SC) and 4-CP (SP). The results showed that C&P had a higher current response and charge transfer capability and lower solution resistance plus charge transfer resistance. The anode bacteria solution of C&P contained more electron carriers (RF, FMN, FAD, CoQ10, and Cyt c). Metagenomic sequencing indicated that the total relative abundance of the microorganisms associated with electricity production (Desulfovibrio, Pseudomonas, Azospirillum, Nocardia, Microbacterium, Delftia, Geobacter, Acinetobacter, Bacillus, and Clostridium) was the highest in C&P (4.24%). However, Microbacterium was abundant in SP (0.12%), which exerted antagonistic effects on other electricigens. Among the 10 electricigens based on gene annotation, C&P had a higher overall relative abundance of the Unigene gene annotated to the KO pathway and CAZy level B compared with SC and SP, which were 1.31% and 0.582% respectively. Unigene153954 (ccmC), Unigene357497 (coxB), and Unigene1033667 (ubiG) were related to the electron carrier Cyt c, electron transfer, and CoQ biosynthesis, respectively. These were annotated to Desulfovibrio, Delftia, and Pseudomonas, respectively. Unigene161312 (AA1) used phenols and other substrates as electron donors and was annotated to Pseudomonas. Other functional carbohydrate enzyme genes (e.g., GT2, GT4, and GH31) used carbohydrates as donors and were annotated to other electricigens. This study provides a theoretical basis for electron transfer to promote the development of CW-MFCs.

The function of the plant cell wall in plant-microbe interactions

The plant cell wall is an interface of plant-microbe interactions. The ability of microbes to decompose cell wall polysaccharides contributes to microbial pathogenicity. Plants have evolved mechanisms to prevent cell wall degradation. However, the role of the cell wall in plant-microbe interactions is not well understood. Here, we discuss four functions of the plant cell wall-physical defence, storage of antimicrobial compounds, production of cell wall-derived elicitors, and provision of carbon sources-in the context of plant-microbe interactions. In addition, we discuss the four families of cell surface receptors associated with plant cell walls (malectin-like receptor kinase family, wall-associated kinase family, leucine-rich repeat receptor-like kinase family, and lysin motif receptor-like kinase family) that have been the subject of several important studies in recent years. This review summarises the findings on both plant cell wall and plant immunity, improving our understanding and may provide impetus to various researchers.

Integrative functional analysis uncovers metabolic differences between Candida species

Candida species are a dominant constituent of the human mycobiome and associated with the development of several diseases. Understanding the Candida species metabolism could provide key insights into their ability to cause pathogenesis. Here, we have developed the BioFung database, providing an efficient annotation of protein-encoding genes. Along, with BioFung, using carbohydrate-active enzyme (CAZymes) analysis, we have uncovered core and accessory features across Candida species demonstrating plasticity, adaption to the environment and acquired features. We show a greater importance of amino acid metabolism, as functional analysis revealed that all Candida species can employ amino acid metabolism. However, metabolomics revealed that only a specific cluster of species (AGAu species—C. albicans, C. glabrata and C. auris) utilised amino acid metabolism including arginine, cysteine, and methionine metabolism potentially improving their competitive fitness in pathogenesis. We further identified critical metabolic pathways in the AGAu cluster with biomarkers and anti-fungal target potential in the CAZyme profile, polyamine, choline and fatty acid biosynthesis pathways. This study, combining genomic analysis, and validation with gene expression and metabolomics, highlights the metabolic diversity with AGAu species that underlies their remarkable ability to dominate they mycobiome and cause disease.

Metabolic differences between Candida species are uncovered using the BioFung database alongside genomic and metabolic analysis.

The Land-Sea Connection: Insights Into the Plant Lineage from a Green Algal Perspective

The colonization of land by plants generated opportunities for the rise of new heterotrophic life forms, including humankind. A unique event underpinned this massive change to earth ecosystems-the advent of eukaryotic green algae. Today, an abundant marine green algal group, the prasinophytes, alongside prasinodermophytes and nonmarine chlorophyte algae, is facilitating insights into plant developments. Genome-level data allow identification of conserved proteins and protein families with extensive modifications, losses, or gains and expansion patterns that connect to niche specialization and diversification. Here, we contextualize attributes according to Viridiplantae evolutionary relationships, starting with orthologous protein families, and then focusing on key elements with marked differentiation, resulting in patchy distributions across green algae and plants. We place attention on peptidoglycan biosynthesis, important for plastid division and walls; phytochrome photosensors that are master regulators in plants; and carbohydrate-active enzymes, essential to all manner of carbohydratebiotransformations. Together with advances in algal model systems, these areas are ripe for discovering molecular roles and innovations within and across plant and algal lineages.

Leaf necrosis resulting from downregulation of poplar glycosyltransferase UGT72A2

Reactive species (RS) causing oxidative stress are unavoidable by-products of various plant metabolic processes, such as photosynthesis, respiration or photorespiration. In leaves, flavonoids scavenge RS produced during photosynthesis and protect plant cells against deleterious oxidative damages. Their biosynthesis and accumulation are therefore under tight regulation at the cellular level. Glycosylation has emerged as an essential biochemical reaction in the homeostasis of various specialized metabolites such as flavonoids. This article provides a functional characterization of the Populus tremula x P. alba (poplar) UGT72A2 coding for a UDP-glycosyltransferase that is localized in the chloroplasts. Compared with the wild type, transgenic poplar lines with decreased expression of UGT72A2 are characterized by reduced growth and oxidative damages in leaves, as evidenced by necrosis, higher content of glutathione and lipid peroxidation products as well as diminished soluble peroxidase activity and NADPH to NADP+ ratio under standard growing conditions. They furthermore display lower pools of phenolics, anthocyanins and total flavonoids but higher proanthocyanidins content. Promoter analysis revealed the presence of cis-elements involved in photomorphogenesis, chloroplast biogenesis and flavonoid biosynthesis. The UGT72A2 is regulated by the poplar MYB119, a transcription factor known to regulate the flavonoid biosynthesis pathway. Phylogenetic analysis and molecular docking suggest that UGT72A2 could glycosylate flavonoids; however, the actual substrate(s) was not consistently evidenced with either in vitro assays nor analyses of glycosylated products in leaves of transgenic poplar overexpressing or downregulated for UGT72A2. This article provides elements highlighting the importance of flavonoid glycosylation regarding protection against oxidative stress in poplar leaves and raises new questions about the link between this biochemical reaction and regulation of the redox homeostasis system.

Combination of system biology and classical approaches for developing biorefinery relevant lignocellulolytic Rasamsonia emersonii strain

The objective of this study was to develop thermophilic fungus Rasamsonia emersonii using integrated system biology tools (genomics, proteomics and transcriptional analysis) in combination with classical strain breeding approaches. Developed hyper cellulolytic mutant strain M36 showed endoglucanase (476.35 U/ml), β-glucosidase (70.54 U/ml), cellobiohydrolase (15.17 U/ml), FPase (4.89 U/ml) and xylanase (485.21 U/ml) on cellulose/gram flour based production medium. Comparison of the expression profile at proteome and transcriptional level of the developed strain and wild type parent gave detailed insight into the up-regulation of different CAZymes including glycosyl hydrolases (GH5, GH6, GH7, GH3, GH10) and auxiliary enzymes (lytic polysaccharide monooxygenase, swollenin) at system level. Furthermore, the potential of lignocellulolytic enzyme produced by the developed strain and custom designed cocktail spiked with heterologously expressed lytic polysaccharide monooxygenase from Mycothermus thermophiloides were analyzed for the hydrolysis of biorefinery relevant unwashed pretreated rice straw slurry (PRAJ and IOCL) @17% substrate loading rate.

Plant-microbe interactions in the apoplast: Communication at the plant cell wall

The apoplast is a continuous plant compartment that connects cells between tissues and organs and is one of the first sites of interaction between plants and microbes. The plant cell wall occupies most of the apoplast and is composed of polysaccharides and associated proteins and ions. This dynamic part of the cell constitutes an essential physical barrier and a source of nutrients for the microbe. At the same time, the plant cell wall serves important functions in the interkingdom detection, recognition, and response to other organisms. Thus, both plant and microbe modify the plant cell wall and its environment in versatile ways to benefit from the interaction. We discuss here crucial processes occurring at the plant cell wall during the contact and communication between microbe and plant. Finally, we argue that these local and dynamic changes need to be considered to fully understand plant-microbe interactions.

Metabolic Pathways Involved in the Drought Stress Response of Nitraria tangutorum as Revealed by Transcriptome Analysis

Drought resistance in plants is controlled by multiple genes. To identify the genes that mediate drought stress responses and to assess the associated metabolic pathways in the desert shrub Nitraria tangutorum, we conducted a transcriptome analysis of plants under control (maximum field capacity) and drought (20% of the maximum field capacity) conditions. We analyzed differentially expressed genes (DEGs) of N. tangutorum and their enrichment in the KEGG metabolic pathways database, and explored the molecular biological mechanisms underlying the answer to its drought tolerance. Between the control and drought groups, 119 classified metabolic pathways annotated 3047 DEGs in the KEGG database. For drought tolerance, nitrate reductase (NR) gene expression was downregulated, indicating that NR activity was decreased to improve drought tolerance. In ammonium assimilation, drought stress inhibited glutamine formation. Protochlorophyllide reductase (1.3.1.33) expression was upregulated to promote chlorophyll a synthesis, whereas divinyl reductase (1.3.1.75) expression was downregulated to inhibit chlorophyll-ester a synthesis. The expression of the chlorophyll synthase (2.5.1.62) gene was downregulated, which affected the synthesis of chlorophyll a and b. Overall, drought stress appeared to improve the ability to convert chlorophyll b into chlorophyll a. Our data serve as a theoretical foundation for further elucidating the growth regulatory mechanism of desert xerophytes, thereby facilitating the development and cultivation of new, drought-resistant genotypes for the purpose of improving desert ecosystems.

UGT72, a Major Glycosyltransferase Family for Flavonoid and Monolignol Homeostasis in Plants

Simple Summary

Phenylpropanoids are specialized metabolites playing crucial roles in plant developmental processes and in plant defense towards pathogens. The attachment of sugar moieties to these small hydrophobic molecules renders them more hydrophilic and increases their solubility. The UDP-glycosyltransferase 72 family (UGT72) of plants has been shown to glycosylate mainly two classes of phenylpropanoids, (i) the monolignols that are the building blocks of lignin, the second most abundant polymer after cellulose, and (ii) the flavonoids, which play determinant roles in plant interactions with other organisms and in response to stress. The purpose of this review is to bring an overview of the current knowledge of the UGT72 family and to highlight its role in the homeostasis of these molecules. Potential applications in pharmacology and in wood, paper pulp, and bioethanol production are given within the perspectives.

Abstract

Plants have developed the capacity to produce a diversified range of specialized metabolites. The glycosylation of those metabolites potentially decreases their toxicity while increasing their stability and their solubility, modifying their transport and their storage. The UGT, forming the largest glycosyltransferase superfamily in plants, combine enzymes that glycosylate mainly hormones and phenylpropanoids by using UDP-sugar as a sugar donor. Particularly, members of the UGT72 family have been shown to glycosylate the monolignols and the flavonoids, thereby being involved in their homeostasis. First, we explore primitive UGTs in algae and liverworts that are related to the angiosperm UGT72 family and their role in flavonoid homeostasis. Second, we describe the role of several UGT72s glycosylating monolignols, some of which have been associated with lignification. In addition, the role of other UGT72 members that glycosylate flavonoids and are involved in the development and/or stress response is depicted. Finally, the importance to explore the subcellular localization of UGTs to study their roles in planta is discussed.

Two conifer GUX clades are responsible for distinct glucuronic acid patterns on xylan

Wood of coniferous trees (softwood), is a globally significant carbon sink and an important source of biomass. Despite that, little is known about the genetic basis of softwood cell wall biosynthesis. Branching of xylan, one of the main hemicelluloses in softwood secondary cell walls, with glucuronic acid (GlcA) is critical for biomass recalcitrance. Here, we investigate the decoration patterns of xylan by conifer GlucUronic acid substitution of Xylan (GUX) enzymes. Through molecular phylogenetics we identify two distinct conifer GUX clades. Using transcriptional profiling we show that the genes are preferentially expressed in secondary cell wall forming tissues. With in vitro and in planta assays we demonstrate that conifer GUX enzymes from both clades are active glucuronyltransferases. Conifer GUX enzymes from each clade have different specific activities. While members of clade one add evenly spaced GlcA branches, the members of clade two are also capable of glucuronidating two consecutive xyloses. Importantly, these types of xylan patterning are present in softwood. As xylan patterning might modulate xylan-cellulose and xylan-lignin interactions, our results further the understanding of softwood cell wall biosynthesis and provide breeding or genetic engineering targets that can be used to modify softwood properties.

Pectin and Xylan Biosynthesis in Poplar: Implications and Opportunities for Biofuels Production

A potential method by which society's reliance on fossil fuels can be lessened is via the large-scale utilization of biofuels derived from the secondary cell walls of woody plants; however, there remain a number of technical challenges to the large-scale production of biofuels. Many of these challenges emerge from the underlying complexity of the secondary cell wall. The challenges associated with lignin have been well explored elsewhere, but the dicot cell wall components of hemicellulose and pectin also present a number of difficulties. Here, we provide an overview of the research wherein pectin and xylan biosynthesis has been altered, along with investigations on the function of irregular xylem 8 (IRX8) and glycosyltransferase 8D (GT8D), genes putatively involved in xylan and pectin synthesis. Additionally, we provide an analysis of the evidence in support of two hypotheses regarding GT8D and conclude that while there is evidence to lend credence to these hypotheses, there are still questions that require further research and examination.

Physiological Characteristic Changes and Full-Length Transcriptome of Rose (Rosa chinensis) Roots and Leaves in Response to Drought Stress

Rose (Rosa chinensis) is the most important ornamental crops worldwide. However, the physiological and molecular mechanism of rose under drought stress remains elusive. In this study, we analyzed the changes of photosynthetic and phytohormone levels in the leaves and roots of rose seedlings grown under control (no drought), mild drought (MD) and severe drought stress. The total chlorophyll content and water use efficiency were significantly enhanced under MD in rose leaves. In addition, the concentration of ABA was higher in the leaves compared to the roots, whereas the roots accumulated more IAA, methylindole-3-acetic acid and indole-3-propionic acid. We also constructed the first full-length transcriptome for rose, and identified 96,201,862 full-length reads of average length 1,149 bp that included 65,789 novel transcripts. A total of 3,657 and 4,341 differentially expressed genes (DEGs) were identified in rose leaves and roots respectively. KEGG pathway analysis showed enrichment of plant hormone, signal transduction and photosynthesis are among the DEGs. 42,544 alternatively spliced isoforms were also identified, and alternative 3' splice site was the major alternative splicing (AS) event among the DEGs. Variations in the AS patterns of three genes between leaves and roots indicated the possibility of tissue-specific posttranscriptional regulation in response to drought stress. Furthermore, 2,410 novel long non-coding RNAs were detected that may participate in regulating the drought-induced DEGs. Our findings identified previously unknown splice sites and new genes in the rose transcriptome, and elucidated the drought stress-responsive genes as well as their intricate regulatory networks.

Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development

Malectin domain (MD) is a ligand-binding protein motif of pro- and eukaryotes. It is particularly abundant in Viridiplantae, where it occurs as either a single (MD, PF11721) or tandemly duplicated domain (PF12819) called malectin-like domain (MLD). In herbaceous plants, MD- or MLD-containing proteins (MD proteins) are known to regulate development, reproduction, and resistance to various stresses. However, their functions in woody plants have not yet been studied. To unravel their potential role in wood development, we carried out genome-wide identification of MD proteins in the model tree species black cottonwood (Populus trichocarpa), and analyzed their expression and co-expression networks. P. trichocarpa had 146 MD genes assigned to 14 different clades, two of which were specific to the genus Populus. 87% of these genes were located on chromosomes, the rest being associated with scaffolds. Based on their protein domain organization, and in agreement with the exon-intron structures, the MD genes identified here could be classified into five superclades having the following domains: leucine-rich repeat (LRR)-MD-protein kinase (PK), MLD-LRR-PK, MLD-PK (CrRLK1L), MLD-LRR, and MD-Kinesin. Whereas the majority of MD genes were highly expressed in leaves, particularly under stress conditions, eighteen showed a peak of expression during secondary wall formation in the xylem and their co-expression networks suggested signaling functions in cell wall integrity, pathogen-associated molecular patterns, calcium, ROS, and hormone pathways. Thus, P. trichocarpa MD genes having different domain organizations comprise many genes with putative foliar defense functions, some of which could be specific to Populus and related species, as well as genes with potential involvement in signaling pathways in other tissues including developing wood.

The Spatio-Temporal Distribution of Cell Wall-Associated Glycoproteins During Wood Formation in Populus

Plant cell wall associated hydroxyproline-rich glycoproteins (HRGPs) are involved in several aspects of plant growth and development, including wood formation in trees. HRGPs such as arabinogalactan-proteins (AGPs), extensins (EXTs), and proline rich proteins (PRPs) are important for the development and architecture of plant cell walls. Analysis of publicly available gene expression data revealed that many HRGP encoding genes show tight spatio-temporal expression patterns in the developing wood of Populus that are indicative of specific functions during wood formation. Similar results were obtained for the expression of glycosyl transferases putatively involved in HRGP glycosylation. In situ immunolabelling of transverse wood sections using AGP and EXT antibodies revealed the cell type specificity of different epitopes. In mature wood AGP epitopes were located in xylem ray cell walls, whereas EXT epitopes were specifically observed between neighboring xylem vessels, and on the ray cell side of the vessel walls, likely in association with pits. Molecular mass and glycan analysis of AGPs and EXTs in phloem/cambium, developing xylem, and mature xylem revealed clear differences in glycan structures and size between the tissues. Separation of AGPs by agarose gel electrophoresis and staining with β-D-glucosyl Yariv confirmed the presence of different AGP populations in phloem/cambium and xylem. These results reveal the diverse changes in HRGP-related processes that occur during wood formation at the gene expression and HRGP glycan biosynthesis levels, and relate HRGPs and glycosylation processes to the developmental processes of wood formation.

MYB Transcription Factor161 Mediates Feedback Regulation of Secondary wall-associated NAC-Domain1 Family Genes for Wood Formation

Wood formation is a complex process that involves cell differentiation, cell expansion, secondary wall deposition, and programmed cell death. We constructed a four-layer wood formation transcriptional regulatory network (TRN) in Populus trichocarpa (black cottonwood) that has four Secondary wall-associated NAC-Domain1 (PtrSND1) transcription factor (TF) family members as the top-layer regulators. We characterized the function of a MYB (PtrMYB161) TF in this PtrSND1-TRN, using transgenic P trichocarpa cells and whole plants. PtrMYB161 is a third-layer regulator that directly transactivates five wood formation genes. Overexpression of PtrMYB161 in P. trichocarpa (OE-PtrMYB161) led to reduced wood, altered cell type proportions, and inhibited growth. Integrative analysis of wood cell-based chromatin-binding assays with OE-PtrMYB161 transcriptomics revealed a feedback regulation system in the PtrSND1-TRN, where PtrMYB161 represses all four top-layer regulators and one second-layer regulator, PtrMYB021, possibly affecting many downstream TFs in, and likely beyond, the TRN, to generate the observed phenotypic changes. Our data also suggested that the PtrMYB161's repressor function operates through interaction of the base PtrMYB161 target-binding system with gene-silencing cofactors. PtrMYB161 protein does not contain any known negative regulatory domains. CRISPR-based mutants of PtrMYB161 in P. trichocarpa exhibited phenotypes similar to the wild type, suggesting that PtrMYB161's activator functions are redundant among many TFs. Our work demonstrated that PtrMYB161 binds to multiple sets of target genes, a feature that allows it to function as an activator as well as a repressor. The balance of the two functions may be important to the establishment of regulatory homeostasis for normal growth and development.

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.

Arabidopsis XTH4 and XTH9 Contribute to Wood Cell Expansion and Secondary Wall Formation

Xyloglucan is the major hemicellulose of dicotyledon primary cell walls, affecting the load-bearing framework with the participation of xyloglucan endo-transglycosylase/hydrolases (XTHs). We used loss- and gain-of function approaches to study functions of XTH4 and XTH9 abundantly expressed in cambial regions during secondary growth of Arabidopsis (Arabidopsis thaliana). In secondarily thickened hypocotyls, these enzymes had positive effects on vessel element expansion and fiber intrusive growth. They also stimulated secondary wall thickening but reduced secondary xylem production. Cell wall analyses of inflorescence stems revealed changes in lignin, cellulose, and matrix sugar composition indicating an overall increase in secondary versus primary walls in mutants, indicative of higher xylem production compared with the wild type (since secondary walls were thinner). Intriguingly, the number of secondary cell wall layers compared with the wild type was increased in xth9 and reduced in xth4, whereas the double mutant xth4x9 displayed an intermediate number of layers. These changes correlated with specific Raman signals from the walls, indicating changes in lignin and cellulose. Secondary walls were affected also in the interfascicular fibers, where neither XTH4 nor XTH9 was expressed, indicating that these effects were indirect. Transcripts involved in secondary wall biosynthesis and cell wall integrity sensing, including THESEUS1 and WALL ASSOCIATED KINASE2, were highly induced in the mutants, indicating that deficiency in XTH4 and XTH9 triggers cell wall integrity signaling, which, we propose, stimulates xylem cell production and modulates secondary wall thickening. Prominent effects of XTH4 and XTH9 on secondary xylem support the hypothesis that altered xyloglucan affects wood properties both directly and via cell wall integrity sensing.