Metabolic profiling of two white-rot fungi during 4-hydroxybenzoate conversion reveals biotechnologically relevant biosynthetic pathways

White-rot fungi are efficient organisms for the mineralization of lignin and polysaccharides into CO2 and H2O. Despite their biotechnological potential, WRF metabolism remains underexplored. Building on recent findings regarding the utilization of lignin-related aromatic compounds as carbon sources by WRF, we aimed to gain further insights into these catabolic processes. For this purpose, Trametes versicolor and Gelatoporia subvermispora were incubated in varying conditions - in static and agitation modes and different antioxidant levels - during the conversion of 4-hydroxybenzoic acid (a lignin-related compound) and cellobiose. Their metabolic responses were assessed via transcriptomics, proteomics, lipidomics, metabolomics, and microscopy analyses. These analyses reveal the significant impact of cultivation conditions on sugar and aromatic catabolic pathways, as well as lipid composition of the fungal mycelia. Additionally, this study identifies biosynthetic pathways for the production of extracellular fatty acids and phenylpropanoids - both products with relevance in biotechnological applications - and provides insights into carbon fate in nature.

Epigenetic regulation of lignin biosynthesis in wood formation

Lignin, a major wood component, is the key limiting factor for wood conversion efficiency. Its biosynthesis is controlled by transcriptional regulatory networks involving transcription factor (TF)-DNA interactions. However, the epigenetic mechanisms underlying these interactions in lignin biosynthesis remain largely unknown. Here, using yeast one-hybrid, chromatin immunoprecipitation, and electrophoretic mobility shift assays, we identified that PtrbZIP44-A1, a key wood-forming TF, directly interacts with the promoters of PtrCCoAOMT2 and PtrCCR2, genes involved in the monolignol biosynthetic pathway. We used yeast two-hybrid, bimolecular fluorescence complementation, biochemical analyses, transient and CRISPR-mediated transgenesis in Populus trichocarpa to demonstrate that PtrHDA15, a histone deacetylase, acts as an epigenetic inhibitor and is recruited by PtrbZIP44-A1 for chromatin histone modifications to repress PtrCCoAOMT2 and PtrCCR2, leading to reduced lignin deposition. In transgenic lines overexpressing PtrbZIP44-A1 or PtrHDA15, histone acetylation at the promoters of PtrCCoAOMT2 and PtrCCR2 decreased, reducing their expression and lignin content. Conversely, in loss-of-function ptrbzip44-a1 and ptrhda15 mutants, histone acetylation levels at PtrCCoAOMT2 and PtrCCR2 promoters increased, enhancing target gene expression and lignin content. Our study uncovered an epigenetic mechanism that suppresses lignin biosynthesis. This finding may help fill a knowledge gap between epigenetic regulation and lignin biosynthesis during wood formation in Populus.

Simultaneous Down-Regulation of Dominant Cinnamoyl CoA Reductase and Cinnamyl Alcohol Dehydrogenase Dramatically Altered Lignin Content in Mulberry

Mulberry (Morus alba L.) is a significant economic tree species in China. The lignin component serves as a critical limiting factor that impacts both the forage quality and the conversion efficiency of mulberry biomass into biofuel. Cinnamoyl CoA reductase (CCR; EC 1.21.1.44) and cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.95) are the key enzymes that catalyze the final two reductive steps in the biosynthesis of monolignols. In this study, we conducted a comprehensive functional analysis to validate the predominant CCR genes involved in monolignol biosynthesis. In this study, we initially validated the predominant CCR genes implicated in monolignol biosynthesis through an extensive functional analysis. Phylogenetic analysis, tissue-specific expression profiling and enzymatic assays indicated that MaCCR1 is the authentic CCR involved in lignin biosynthesis. Furthermore, the expression level of MaCCR1 exhibited a significant positive correlation with lignin content, and the down-regulation of MaCCR1 via virus-induced gene silencing resulted in altered lignin content in mulberry. The down-regulation of MaCCR1 and MaCAD3/4, both individually and concurrently, exhibited markedly different effects on lignin content and mulberry growth. Specifically, the simultaneous down-regulation of MaCCR1 and MaCAD3/4 significantly altered lignin content in mulberry, resulting in dwarfism of the plants. Conversely, the down-regulation of MaCAD3/4 alone not only decreased lignin content but also led to an increase in biomass. These findings offer compelling evidence elucidating the roles of MaCCRs in mulberry and identify specific target genes, thereby providing a crucial foundation for the genetic modification of lignin biosynthesis.

Comprehensive, Genome-Wide Identification and Expression Analyses of Phenylalanine Ammonia-Lyase Family under Abiotic Stresses in Brassica oleracea

Phenylalanine ammonia-lyase (PAL) acts as the rate-limiting enzyme for anthocyanin biosynthesis through the phenylpropanoid pathway, a crucial component of plant secondary metabolism. The PAL gene family plays a crucial role in plants' defense and stress responses, but its in silico identification and expression analyses in Brassica oleracea under different abiotic stresses remain unexplored. In this study, nine BolPAL, seven BrPAL, four AtPAL, and seventeen BnPAL genes were obtained from the genomes of B. oleracea, Brassica rapa, Arabidopsis thaliana, and Brassica napus, respectively. Segmental duplication and purifying selection are the causes of the BolPAL gene's amplification and evolution. The BolPAL genes with comparable intron-exon architectures and motifs were grouped together in the same clade. Three categories comprised the cis-regulatory elements: abiotic stressors, phytohormones, and light. According to the results of the qRT-PCR experiments, the majority of the BolPAL genes were expressed highly under MeJA, a low temperature, and a high temperature, and they were downregulated under ABA. Under white light (100 µmol m-2 s-1) with 50, 100, or 150 µmol m-2 s-1 far-red (FR), only a small number of the PAL genes were expressed at 50 and 100 µmol m-2 s-1 FR, while the majority of the PAL genes were slightly elevated at 150 µmol m-2 s-1 FR. This work offers a theoretical foundation for molecular breeding research to investigate the role of BolPAL genes and their role in anthocyanin biosynthesis.

Exploring Lignin Biosynthesis Genes in Rice: Evolution, Function, and Expression

Lignin is nature's second most abundant vascular plant biopolymer, playing significant roles in mechanical support, water transport, and stress responses. This study identified 90 lignin biosynthesis genes in rice based on phylogeny and motif constitution, and they belong to PAL, C4H, 4CL, HCT, C3H, CCoAOMT, CCR, F5H, COMT, and CAD families. Duplication events contributed largely to the expansion of these gene families, such as PAL, CCoAOMT, CCR, and CAD families, mainly attributed to tandem and segmental duplication. Microarray data of 33 tissue samples covering the entire life cycle of rice suggested fairly high PAL, HCT, C3H, CCoAOMT, CCR, COMT, and CAD gene expressions and rather variable C4H, 4CL, and F5H expressions. Some members of lignin-related genes (OsCCRL11, OsHCT1/2/5, OsCCoAOMT1/3/5, OsCOMT, OsC3H, OsCAD2, and OsPAL1/6) were expressed in all tissues examined. The expression patterns of lignin-related genes can be divided into two major groups with eight subgroups, each showing a distinct co-expression in tissues representing typically primary and secondary cell wall constitutions. Some lignin-related genes were strongly co-expressed in tissues typical of secondary cell walls. Combined HPLC analysis showed increased lignin monomer (H, G, and S) contents from young to old growth stages in five genotypes. Based on 90 genes' microarray data, 27 genes were selected for qRT-PCR gene expression analysis. Four genes (OsPAL9, OsCAD8C, OsCCR8, and OsCOMTL4) were significantly negatively correlated with lignin monomers. Furthermore, eleven genes were co-expressed in certain genotypes during secondary growth stages. Among them, six genes (OsC3H, OsCAD2, OsCCR2, OsCOMT, OsPAL2, and OsPAL8) were overlapped with microarray gene expressions, highlighting their importance in lignin biosynthesis.

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.

C2H2 zinc finger protein PagIDD15A regulates secondary wall thickening and lignin biosynthesis in poplar

Wood production is largely determined by the activity of cambial cell proliferation, and the secondary cell wall (SCW) thickening of xylem cells determines the wood property. In this study, we identified an INDETERMINATE DOMAIN (IDD) type C2H2 zinc finger transcription factor PagIDD15A as a regulator of wood formation in Populus alba × Populus glandulosa. Downregulation of PagIDD15A expression by RNA interference (RNAi) inhibited xylem development and xylem cell secondary wall thickening. RNA-seq analysis showed that PagPAL1, PagCCR2 and PagCCoAOMT1 were downregulated in the differentiating xylem of the PagIDD15A-RNAi transgenic plants, showing that PagIDD15A may regulate SCW biosynthesis through inhibiting lignin biosynthesis. The downregulation of PagVND6-B2, PagMYB10 and PagMYC4 and upregulation of PagWRKY12 in the differentiating xylem of RNAi transgenic plants suggest that PagIDD15A may also regulate these transcription factor (TF) genes to affect SCW thickening. RT-qPCR analysis in the phloem-cambium of RNAi transgenic demonstrates that PagIDD15A may regulate the expression of the genes associated with cell proliferation, including, PagSHR (SHORTROOT), PagSCR (SCARECROW), PagCYCD3;1 (CYCLIN D3;1) and PagSMR4 (SIAMESE-RELATED4), to affect the cambial activity. This study provides the knowledge of the IDD-type C2H2 zinc finger protein in regulating wood formation.

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

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

Blue light photoreceptor cryptochrome 1 promotes wood formation and anthocyanin biosynthesis in Populus

Blue light photoreceptor cryptochrome 1 (CRY1) in herbaceous plants plays crucial roles in various developmental processes, including cotyledon expansion, hypocotyl elongation and anthocyanin biosynthesis. However, the function of CRY1 in perennial trees is unclear. In this study, we identified two ortholog genes of CRY1 (PagCRY1a and PagCRY1b) from Populus, which displayed high sequence similarity to Arabidopsis CRY1. Overexpression of PagCRY1 substantially inhibited plant growth and promoted secondary xylem development in Populus, while CRISPR/Cas9-mediated knockout of PagCRY1 enhanced plant growth and delayed secondary xylem development. Moreover, overexpression of PagCRY1 dramatically increased anthocyanin accumulation. The further analysis supported that PagCRY1 functions specifically in response to blue light. Taken together, our results demonstrated that modulating the expression of blue light photoreceptor CRY1 ortholog gene in Populus could significantly influence plant biomass production and the process of wood formation, laying a foundation for further investigating the light-regulated tree growth.

Systematic characterization of gene families and functional analysis of PvRAS3 and PvRAS4 involved in rosmarinic acid biosynthesis in Prunella vulgaris

Prunella vulgaris is an important material for Chinese medicines with rosmarinic acid (RA) as its index component. Based on the chromosome-level genome assembly we obtained recently, 51 RA biosynthesis-related genes were identified. Sequence feature, gene expression pattern and phylogenetic relationship analyses showed that 17 of them could be involved in RA biosynthesis. In vitro enzymatic assay showed that PvRAS3 catalyzed the condensation of p-coumaroyl-CoA and caffeoyl-CoA with pHPL and DHPL. Its affinity toward p-coumaroyl-CoA was higher than caffeoyl-CoA. PvRAS4 catalyzed the condensation of p-coumaroyl-CoA with pHPL and DHPL. Its affinity toward p-coumaroyl-CoA was lower than PvRAS3. UPLC and LC-MS/MS analyses showed the existence of RA, 4-coumaroyl-3',4'-dihydroxyphenyllactic acid, 4-coumaroyl-4'-hydroxyphenyllactic acid and caffeoyl-4'-hydroxyphenyllactic acid in P. vulgaris. Generation and analysis of pvras3 homozygous mutants showed significant decrease of RA, 4-coumaroyl-3',4'-dihydroxyphenyllactic acid, 4-coumaroyl-4'-hydroxyphenyllactic acid and caffeoyl-4'-hydroxyphenyllactic acid and significant increase of DHPL and pHPL. It suggests that PvRAS3 is the main enzyme catalyzing the condensation of acyl donors and acceptors during RA biosynthesis. The role of PvRAS4 appears minor. The results provide significant information for quality control of P. vulgaris medicinal materials.

Genome-wide identification, expression profiling, and protein interaction analysis of the CCoAOMT gene family in the tea plant (Camellia sinensis)

BACKGROUND

The caffeoyl-CoA-O methyltransferase (CCoAOMT) family plays a crucial role in the oxidative methylation of phenolic substances and is involved in various plant processes, including growth, development, and stress response. However, there is a limited understanding of the interactions among CCoAOMT protein members in tea plants.

RESULTS

In this study, we identified 10 members of the CsCCoAOMT family in the genome of Camellia sinensis (cultivar 'HuangDan'), characterized by conserved gene structures and motifs. These CsCCoAOMT members were located on six different chromosomes (1, 2, 3, 4, 6, and 14). Based on phylogenetic analysis, CsCCoAOMT can be divided into two groups: I and II. Notably, the CsCCoAOMT members of group Ia are likely to be candidate genes involved in lignin biosynthesis. Moreover, through the yeast two-hybrid (Y2H) assay, we established protein interaction networks for the CsCCoAOMT family, revealing 9 pairs of members with interaction relationships.

CONCLUSIONS

We identified the CCoAOMT gene family in Camellia sinensis and conducted a comprehensive analysis of their classifications, phylogenetic and synteny relationships, gene structures, protein interactions, tissue-specific expression patterns, and responses to various stresses. Our findings shed light on the evolution and composition of CsCCoAOMT. Notably, the observed interaction among CCoAOMT proteins suggests the potential formation of the O-methyltransferase (OMT) complex during the methylation modification process, expanding our understanding of the functional roles of this gene family in diverse biological processes.

Genome-Wide Characterization of the Phenylalanine Ammonia-Lyase Gene Family and Their Potential Roles in Response to Aspergillus flavus L. Infection in Cultivated Peanut (Arachis hypogaea L.)

Phenylalanine ammonia-lyase (PAL) is an essential enzyme in the phenylpropanoid pathway, in which numerous aromatic intermediate metabolites play significant roles in plant growth, adaptation, and disease resistance. Cultivated peanuts are highly susceptible to Aspergillus flavus L. infection. Although PAL genes have been characterized in various major crops, no systematic studies have been conducted in cultivated peanuts, especially in response to A. flavus infection. In the present study, a systematic genome-wide analysis was conducted to identify PAL genes in the Arachis hypogaea L. genome. Ten AhPAL genes were distributed unevenly on nine A. hypogaea chromosomes. Based on phylogenetic analysis, the AhPAL proteins were classified into three groups. Structural and conserved motif analysis of PAL genes in A. hypogaea revealed that all peanut PAL genes contained one intron and ten motifs in the conserved domains. Furthermore, synteny analysis indicated that the ten AhPAL genes could be categorized into five pairs and that each AhPAL gene had a homologous gene in the wild-type peanut. Cis-element analysis revealed that the promoter region of the AhPAL gene family was rich in stress- and hormone-related elements. Expression analysis indicated that genes from Group I (AhPAL1 and AhPAL2), which had large number of ABRE, WUN, and ARE elements in the promoter, played a strong role in response to A. flavus stress.

Integrated physiological, metabolomic, and transcriptomic analyses elucidate the regulation mechanisms of lignin synthesis under osmotic stress in alfalfa leaf (Medicago sativa L.)

Alfalfa, an essential forage crop known for its high yield, nutritional value, and strong adaptability, has been widely cultivated worldwide. The yield and quality of alfalfa are frequently jeopardized due to environmental degradation. Lignin, a constituent of the cell wall, enhances plant resistance to abiotic stress, which often causes osmotic stress in plant cells. However, how lignin responds to osmotic stress in leaves remains unclear. This study explored the effects of osmotic stress on lignin accumulation and the contents of intermediate metabolites involved in lignin synthesis in alfalfa leaves. Osmotic stress caused an increase in lignin accumulation and the alteration of core enzyme activities and gene expression in the phenylpropanoid pathway. We identified five hub genes (CSE, CCR, CADa, CADb, and POD) and thirty edge genes (including WRKYs, MYBs, and UBPs) by integrating transcriptome and metabolome analyses. In addition, ABA and ethylene signaling induced by osmotic stress regulated lignin biosynthesis in a contradictory way. These findings contribute to a new theoretical foundation for the breeding of high-quality and resistant alfalfa varieties.

Deciphering the intricate hierarchical gene regulatory network: unraveling multi-level regulation and modifications driving secondary cell wall formation

Wood quality is predominantly determined by the amount and the composition of secondary cell walls (SCWs). Consequently, unraveling the molecular regulatory mechanisms governing SCW formation is of paramount importance for genetic engineering aimed at enhancing wood properties. Although SCW formation is known to be governed by a hierarchical gene regulatory network (HGRN), our understanding of how a HGRN operates and regulates the formation of heterogeneous SCWs for plant development and adaption to ever-changing environment remains limited. In this review, we examined the HGRNs governing SCW formation and highlighted the significant key differences between herbaceous Arabidopsis and woody plant poplar. We clarified many confusions in existing literatures regarding the HGRNs and their orthologous gene names and functions. Additionally, we revealed many network motifs including feed-forward loops, feed-back loops, and negative and positive autoregulation in the HGRNs. We also conducted a thorough review of post-transcriptional and post-translational aspects, protein-protein interactions, and epigenetic modifications of the HGRNs. Furthermore, we summarized how the HGRNs respond to environmental factors and cues, influencing SCW biosynthesis through regulatory cascades, including many regulatory chains, wiring regulations, and network motifs. Finally, we highlighted the future research directions for gaining a further understanding of molecular regulatory mechanisms underlying SCW formation.

The effect of nitrogen source and levels on hybrid aspen tree physiology and wood formation

Nitrogen can be taken up by trees in the form of nitrate, ammonium and amino acids, but the influence of the different forms on tree growth and development is poorly understood in angiosperm species like Populus. We studied the effects of both organic and inorganic forms of nitrogen on growth and wood formation of hybrid aspen trees in experimental conditions that allowed growth under four distinct steady-state nitrogen levels. Increased nitrogen availability had a positive influence on biomass accumulation and the radial dimensions of both xylem vessels and fibers, and a negative influence on wood density. An optimal level of nitrogen availability was identified where increases in biomass accumulation outweighed decreases in wood density. None of these responses depended on the source of nitrogen except for shoot biomass accumulation, which was stimulated more by treatments complemented with nitrate than by ammonium alone or the organic source arginine. The most striking difference between the nitrogen sources was the effect on lignin composition, whereby the abundance of H-type lignin increased only in the presence of nitrate. The differential effect of nitrate is possibly related to the well-known role of nitrate as a signaling compound. RNA-sequencing revealed that while the lignin-biosynthetic genes did not significantly (FDR <0.01) respond to added NO3 - , the expression of several laccases, catalysing lignin polymerization, was dependent on N-availability. These results reveal a unique role of nitrate in wood formation and contribute to the knowledge basis for decision-making in utilizing hybrid aspen as a bioresource.

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

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

Integrative analysis of metabolome, proteome, and transcriptome for identifying genes influencing total lignin content in Populus trichocarpa

Lignin, a component of plant cell walls, possesses significant research potential as a renewable energy source to replace carbon-based products and as a notable pollutant in papermaking processes. The monolignol biosynthetic pathway has been elucidated and it is known that not all monolignol genes influence the total lignin content. However, it remains unclear which monolignol genes are more closely related to the total lignin content and which potential genes influence the total lignin content. In this study, we present a combination of t-test, differential gene expression analysis, correlation analysis, and weighted gene co-expression network analysis to identify genes that regulate the total lignin content by utilizing multi-omics data from transgenic knockdowns of the monolignol genes that includes data related to the transcriptome, proteome, and total lignin content. Firstly, it was discovered that enzymes from the PtrPAL, Ptr4CL, PtrC3H, and PtrC4H gene families are more strongly correlated with the total lignin content. Additionally, the co-downregulation of three genes, PtrC3H3, PtrC4H1, and PtrC4H2, had the greatest impact on the total lignin content. Secondly, GO and KEGG analysis of lignin-related modules revealed that the total lignin content is not only influenced by monolignol genes, but also closely related to genes involved in the "glutathione metabolic process", "cellular modified amino acid metabolic process" and "carbohydrate catabolic process" pathways. Finally, the cinnamyl alcohol dehydrogenase genes CAD1, CADL3, and CADL8 emerged as potential contributors to total lignin content. The genes HYR1 (UDP-glycosyltransferase superfamily protein) and UGT71B1 (UDP-glucosyltransferase), exhibiting a close relationship with coumarin, have the potential to influence total lignin content by regulating coumarin metabolism. Additionally, the monolignol genes PtrC3H3, PtrC4H1, and PtrC4H2, which belong to the cytochrome P450 genes, may have a significant impact on the total lignin content. Overall, this study establishes connections between gene expression levels and total lignin content, effectively identifying genes that have a significant impact on total lignin content and offering novel perspectives for future lignin research endeavours.

Non-canonical two-step biosynthesis of anti-oomycete indole alkaloids in Kickxellales

BACKGROUND

Fungi are prolific producers of bioactive small molecules of pharmaceutical or agricultural interest. The secondary metabolism of higher fungi (Dikarya) has been well-investigated which led to > 39,000 described compounds. However, natural product researchers scarcely drew attention to early-diverging fungi (Mucoro- and Zoopagomycota) as they are considered to rarely produce secondary metabolites. Indeed, only 15 compounds have as yet been isolated from the entire phylum of the Zoopagomycota.

RESULTS

Here, we showcase eight species of the order Kickxellales (phylum Zoopagomycota) as potent producers of the indole-3-acetic acid (IAA)-derived compounds lindolins A and B. The compounds are produced both under laboratory conditions and in the natural soil habitat suggesting a specialized ecological function. Indeed, lindolin A is a selective agent against plant-pathogenic oomycetes such as Phytophthora sp. Lindolin biosynthesis was reconstituted in vitro and relies on the activity of two enzymes of dissimilar evolutionary origin: Whilst the IAA-CoA ligase LinA has evolved from fungal 4-coumaryl-CoA synthetases, the subsequently acting IAA-CoA:anthranilate N-indole-3-acetyltransferase LinB is a unique enzyme across all kingdoms of life.

CONCLUSIONS

This is the first report on bioactive secondary metabolites in the subphylum Kickxellomycotina and the first evidence for a non-clustered, two-step biosynthetic route of secondary metabolites in early-diverging fungi. Thus, the generally accepted "gene cluster hypothesis" for natural products needs to be reconsidered for early diverging fungi.

Brassinosteroids affect wood development and properties of Fraxinus mandshurica

Introduction

Xylem development plays a crucial role in wood formation in woody plants. In recent years, there has been growing attention towards the impact of brassinosteroids (BRs) on this xylem development. In the present study, we evaluated the dynamic variation of xylem development in Fraxinus mandshurica (female parent, M8) and a novel interspecific hybrid F. mandshurica × Fraxinus sogdiana (1601) from May to August 2020.

Methods

We obtained RNA-Seq transcriptomes of three tissue types (xylem, phloem, and leaf) to identify the differences in xylem-differentially expressed genes (X-DEGs) and xylem-specifically expressed genes (X-SEGs) in M8 and 1601 variants. We then further evaluated these genes via weighted gene co-expression network analysis (WGCNA) alongside overexpressing FmCPD, a BR biosynthesis enzyme gene, in transient transgenic F. mandshurica.

Results

Our results indicated that the xylem development cycle of 1601 was extended by 2 weeks compared to that of M8. In addition, during the later wood development stages (secondary wall thickening) of 1601, an increased cellulose content (14%) and a reduced lignin content (11%) was observed. Furthermore, vessel length and width increased by 67% and 37%, respectively, in 1601 compared with those of M8. A total of 4589 X-DEGs were identified, including enzymes related to phenylpropane metabolism, galactose metabolism, BR synthesis, and signal transduction pathways. WGCNA identified hub X-SEGs involved in cellulose synthesis and BR signaling in the 1601 wood formation-related module (CESA8, COR1, C3H14, and C3H15); in contrast, genes involved in phenylpropane metabolism were significantly enriched in the M8 wood formation-related module (CCoAOMT and CCR). Moreover, overexpression of FmCPD in transient transgenic F. mandshurica affected the expression of genes associated with lignin and cellulose biosynthesis signal transduction. Finally, BR content was determined to be approximately 20% lower in the M8 xylem than in the 1601 xylem, and the exogenous application of BRs (24-epi brassinolide) significantly increased the number of xylem cell layers and altered the composition of the secondary cell walls in F. mandshurica.

Discussion

Our findings suggest that BR biosynthesis and signaling play a critical role in the differing wood development and properties observed between M8 and 1601 F. mandshurica.

Multiplex CRISPR editing of wood for sustainable fiber production

The domestication of forest trees for a more sustainable fiber bioeconomy has long been hindered by the complexity and plasticity of lignin, a biopolymer in wood that is recalcitrant to chemical and enzymatic degradation. Here, we show that multiplex CRISPR editing enables precise woody feedstock design for combinatorial improvement of lignin composition and wood properties. By assessing every possible combination of 69,123 multigenic editing strategies for 21 lignin biosynthesis genes, we deduced seven different genome editing strategies targeting the concurrent alteration of up to six genes and produced 174 edited poplar variants. CRISPR editing increased the wood carbohydrate-to-lignin ratio up to 228% that of wild type, leading to more-efficient fiber pulping. The edited wood alleviates a major fiber-production bottleneck regardless of changes in tree growth rate and could bring unprecedented operational efficiencies, bioeconomic opportunities, and environmental benefits.

Novel candidate genes for lignin structure identified through genome-wide association study of naturally varying Populus trichocarpa

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (-log₁₀(p-value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

Cotton 4-coumarate-CoA ligase 3 enhanced plant resistance to Verticillium dahliae by promoting jasmonic acid signaling-mediated vascular lignification and metabolic flux

Lignins and their antimicrobial-related polymers cooperatively enhance plant resistance to pathogens. Several isoforms of 4-coumarate-coenzyme A ligases (4CLs) have been identified as indispensable enzymes involved in lignin and flavonoid biosynthetic pathways. However, their roles in plant-pathogen interaction are still poorly understood. This study uncovers the role of Gh4CL3 in cotton resistance to the vascular pathogen Verticillium dahliae. The cotton 4CL3-CRISPR/Cas9 mutant (CR4cl) exhibited high susceptibility to V. dahliae. This susceptibility was most probably due to the reduction in the total lignin content and the biosynthesis of several phenolic metabolites, e.g., rutin, catechin, scopoletin glucoside, and chlorogenic acid, along with jasmonic acid (JA) attenuation. These changes were coupled with a significant reduction in 4CL activity toward p-coumaric acid substrate, and it is likely that recombinant Gh4CL3 could specifically catalyze p-coumaric acid to form p-coumaroyl-coenzyme A. Thus, overexpression of Gh4CL3 (OE4CL) showed increasing 4CL activity that augmented phenolic precursors, cinnamic, p-coumaric, and sinapic acids, channeling into lignin and flavonoid biosyntheses and enhanced resistance to V. dahliae. Besides, Gh4CL3 overexpression activated JA signaling that instantly stimulated lignin deposition and metabolic flux in response to pathogen, which all established an efficient plant defense response system, and inhibited V. dahliae mycelium growth. Our results propose that Gh4CL3 acts as a positive regulator for cotton resistance against V. dahliae by promoting JA signaling-mediated enhanced cell wall rigidity and metabolic flux.

Identification and analysis of lignin biosynthesis genes related to fruit ripening and stress response in banana (Musa acuminata L. AAA group, cv. Cavendish)

BACKGROUND

Lignin is a key component of the secondary cell wall of plants, providing mechanical support and facilitating water transport as well as having important impact effects in response to a variety of biological and abiotic stresses.

RESULTS

In this study, we identified 104 genes from ten enzyme gene families related to lignin biosynthesis in Musa acuminata genome and found the number of MaCOMT gene family was the largest, while MaC3Hs had only two members. MaPALs retained the original members, and the number of Ma4CLs in lignin biosynthesis was significantly less than that of flavonoids. Segmental duplication existed in most gene families, except for MaC3Hs, and tandem duplication was the main way to expand the number of MaCOMTs. Moreover, the expression profiles of lignin biosynthesis genes during fruit development, postharvest ripening stages and under various abiotic and biological stresses were investigated using available RNA-sequencing data to obtain fruit ripening and stress response candidate genes. Finally, a co-expression network of lignin biosynthesis genes was constructed by weighted gene co-expression network analysis to elucidate the lignin biosynthesis genes that might participate in lignin biosynthesis in banana during development and in response to stresses.

CONCLUSION

This study systematically identified the lignin biosynthesis genes in the Musa acuminata genome, providing important candidate genes for further functional analysis. The identification of the major genes involved in lignin biosynthesis in banana provides the basis for the development of strategies to improve new banana varieties tolerant to biological and abiotic stresses with high yield and high quality.

A draft genome of the medicinal plant Cremastra appendiculata (D. Don) provides insights into the colchicine biosynthetic pathway

Cremastra appendiculata (D. Don) Makino is a rare terrestrial orchid with a high market value as an ornamental and Chinese traditional medicinal herb with a wide range of pharmacological properties. The pseudobulbs of C. appendiculata are one of the primary sources of the famous traditional Chinese medicine "Shancigu", which has been clinically used for treating many diseases, especially, as the main component to treat gout. The lack of genetic research and genome data restricts the modern development and clinical use of C. appendiculata. Here, we report a 2.3 Gb chromosome-level genome of C. appendiculata. We identify a series of candidates of 35 candidate genes responsible for colchicine biosynthesis, among which O-methyltransferase (OMT) gene exhibits an important role in colchicine biosynthesis. Co-expression analysis reveal purple and green-yellow module have close relationships with pseudobulb parts and comprise most of the colchicine pathway genes. Overall, our genome data and the candidate genes reported here set the foundation to decipher the colchicine biosynthesis pathways in medicinal plants.

The physiological and transcriptomic study of secondary growth in Neolamarckia cadamba stimulated by the ethylene precursor ACC

Though many biological roles of ethylene have been investigated intensively, the molecular mechanism of ethylene's action in woody plants remains unclear. In this study, we investigated the effects of exogenous 1-aminocyclopropane-1-carboxylic acid (ACC), the precursor of ethylene, on the growth of Neolamarckia cadamba seedlings, a fast-growing tropical tree. After 14 days of ACC treatment, the plants showed a reduced physiological morphology while stem diameter increased; however, this did not occur after the addition of 1-MCP. Meanwhile, the lignin content of N. cadamba also increased. Transcriptome analysis revealed that the expression of the ethylene biosynthesis and signaling genes ACC oxidase (ACO) and ethylene insensitive 3 (EIN3) were up-regulated mainly at the 6th hour and the 3rd day of the ACC treatment, respectively. The transcription levels of transcription factors, mainly in the basic helix-loop-helix (bHLH), ethylene response factor (ERF), WRKY and v-myb avian myeloblastosis viral oncogene homolog (MYB) families, involved in the ethylene signaling and secondary growth also increased significantly. Furthermore, in accordance to the increased lignification of the stem, the transcriptional level of key enzymes in the phenylalanine pathway were elevated after the ACC treatment. Our results revealed the physiological and molecular mechanisms underlying the secondary growth stimulated by exogenous ACC treatment on N. cadamba seedlings.

Biomass formation and sugar release efficiency of Populus modified by altered expression of a NAC transcription factor

Woody biomass is an important feedstock for biofuel production. Manipulation of wood properties that enable efficient conversion of biomass to biofuel reduces cost of biofuel production. Wood cell wall composition is regulated at several levels that involve expression of transcription factors such as wood-/secondary cell wall-associated NAC domains (WND or SND). In Arabidopsis thaliana, SND1 regulates cell wall composition through activation of its down-stream targets such as MYBs. The functional aspects of SND1 homologs in the woody Populus have been studied through transgenic manipulation. In this study, we investigated the role of PdWND1B, Populus SND1 sequence ortholog, in wood formation using transgenic manipulation through over-expression or silencing under the control of a vascular-specific 4-coumarate-CoA ligase (4CL) promoter. As compared with control plants, PdWND1B-RNAi plants were shorter in height, with significantly reduced stem diameter and dry biomass, whereas there were no significant differences in growth and productivity of PdWND1B over-expression plants. Conversely, PdWND1B over-expression lines showed a significant reduction in cellulose and increase in lignin content, whereas there was no significant impact on lignin content of downregulated lines. Stem carbohydrate composition analysis revealed a decrease in glucose, mannose, arabinose, and galactose, but an increase in xylose in the over-expression lines. Transcriptome analysis revealed upregulation of several downstream transcription factors and secondary cell wall related structural genes in the PdWND1B over-expression lines, partly explaining the observed phenotypic changes in cell wall chemistry. Relative to the control, glucose release efficiency and ethanol production from stem biomass was significantly reduced in over-expression lines. Our results show that PdWND1B is an important factor determining biomass productivity, cell wall chemistry and its conversion to biofuels in Populus.

Characterization, Expression Profiling, and Biochemical Analyses of the Cinnamoyl-CoA Reductase Gene Family for Lignin Synthesis in Alfalfa Plants

Cinnamoyl-CoA reductase (CCR) is a pivotal enzyme in plant lignin synthesis, which has a role in plant secondary cell wall development and environmental stress defense. Alfalfa is a predominant legume forage with excellent quality, but the lignin content negatively affects fodder digestibility. Currently, there is limited information on CCR characteristics, gene expression, and its role in lignin metabolism in alfalfa. In this study, we identified 30 members in the CCR gene family of Medicago sativa. In addition, gene structure, conserved motif, and evolution analysis suggested MsCCR1–7 presumably functioned as CCR, while the 23 MsCCR-likes fell into three categories. The expression patterns of MsCCRs/MsCCR-likes suggested their role in plant development, response to environmental stresses, and phytohormone treatment. These results were consistent with the cis-elements in their promoters. Histochemical staining showed that lignin accumulation gradually deepened with the development, which was consistent with gene expression results. Furthermore, recombinant MsCCR1 and MsCCR-like1 were purified and the kinetic parameters were tested under four substrates. In addition, three-dimensional structure models of MsCCR1 and MsCCR-like1 proteins showed the difference in the substrate-binding motif H212(X)2K215R263. These results will be useful for further application for legume forage quality modification and biofuels industry engineering in the future.

Molecular and biochemical characterization of two 4-coumarate: CoA ligase genes in tea plant (Camellia sinensis)

Two 4-coumarate: CoA ligase genes in tea plant involved in phenylpropanoids biosynthesis and response to environmental stresses. Tea plant is rich in flavonoids benefiting human health. Lignin is essential for tea plant growth. Both flavonoids and lignin defend plants from stresses. The biosynthesis of lignin and flavonoids shares a key intermediate, 4-coumaroyl-CoA, which is formed from 4-coumaric acid catalyzed by 4-coumaric acid: CoA ligase (4CL). Herein, we report two 4CL paralogs from tea plant, Cs4CL1 and Cs4CL2, which are a member of class I and II of this gene family, respectively. Cs4CL1 was mainly expressed in roots and stems, while Cs4CL2 was mainly expressed in leaves. The promoter of Cs4CL1 had AC, nine types of light sensitive (LSE), four types of stress-inducible (SIE), and two types of meristem-specific elements (MSE). The promoter of Cs4CL2 also had AC and nine types of LSEs, but only had two types of SIEs and did not have MSEs. In addition, the LSEs varied in the two promoters. Based on the different features of regulatory elements, three stress treatments were tested to understand their expression responses to different conditions. The resulting data indicated that the expression of Cs4CL1 was sensitive to mechanical wounding, while the expression of Cs4CL2 was UV-B-inducible. Enzymatic assays showed that both recombinant Cs4CL1 and Cs4CL2 transformed 4-coumaric acid (CM), ferulic acid (FR), and caffeic acid (CF) to their corresponding CoA ethers. Kinetic analysis indicated that the recombinant Cs4CL1 preferred to catalyze CF, while the recombinant Cs4CL2 favored to catalyze CM. The overexpression of both Cs4CL1 and Cs4CL2 increased the levels of chlorogenic acid and total lignin in transgenic tobacco seedlings. In addition, the overexpression of Cs4CL2 consistently increased the levels of three flavonoid compounds. These findings indicate the differences of Cs4CL1 and Cs4CL2 in the phenylpropanoid metabolism.

Molecular understanding of wood formation in trees

Trees convert and store the majority of their photosynthetic products in wood which is an essential renewable resource much in demand by human society. Formation of wood follows a process of consecutive cell developmental stages, from vascular cambium proliferation, cell expansion and differentiation, secondary cell wall deposition to programmed cell death, which is controlled by the functionality of complex molecular networks. What are the molecular networks involved in wood formation? How do the molecular networks act in a way to generate wood tissue during tree growth? What are the regulatory modules that lead to the formation of various wood characteristics? The answers to these questions are fundamental to understanding how trees grow, as well as how we can genetically engineer trees with desired properties of wood for human needs. In recent years, a great deal of interest has been invested in the elucidation of wood formation at the molecular level. This review summarizes the current state of understanding of the molecular process that guides wood formation in trees.

Transcriptional activation of rice CINNAMOYL-CoA REDUCTASE 10 by OsNAC5, contributes to drought tolerance by modulating lignin accumulation in roots

Drought is a common abiotic stress for terrestrial plants and often affects crop development and yield. Recent studies have suggested that lignin plays a crucial role in plant drought tolerance; however, the underlying molecular mechanisms are still largely unknown. Here, we report that the rice (Oryza sativa) gene CINNAMOYL-CoA REDUCTASE 10 (OsCCR10) is directly activated by the OsNAC5 transcription factor, which mediates drought tolerance through regulating lignin accumulation. CCR is the first committed enzyme in the monolignol synthesis pathway, and the expression of 26 CCR genes was observed to be induced in rice roots under drought. Subcellular localisation assays revealed that OsCCR10 is a catalytically active enzyme that is localised in the cytoplasm. The OsCCR10 transcript levels were found to increase in response to abiotic stresses, such as drought, high salinity, and abscisic acid (ABA), and transcripts were detected in roots at all developmental stages. In vitro enzyme activity and in vivo lignin composition assay suggested that OsCCR10 is involved in H- and G-lignin biosynthesis. Transgenic rice plants overexpressing OsCCR10 showed improved drought tolerance at the vegetative stages of growth, as well as higher photosynthetic efficiency, lower water loss rates, and higher lignin content in roots compared to non-transgenic (NT) controls. In contrast, CRISPR/Cas9-mediated OsCCR10 knock-out mutants exhibited reduced lignin accumulation in roots and less drought tolerance. Notably, transgenic rice plants with root-preferential overexpression of OsCCR10 exhibited higher grain yield than NT controls plants under field drought conditions, indicating that lignin biosynthesis mediated by OsCCR10 contributes to drought tolerance.

Seasonal changes in cambium activity from active to dormant stage affect the formation of secondary xylem in Pinus tabulaeformis Carr

Understanding the changing patterns of vascular cambium during seasonal cycles is crucial to reveal the mechanisms that control cambium activity and wood formation, but this area has been underexplored, especially in conifers. Here, we quantified the changing cellular morphology patterns of cambial zones during the active, transition and dormant stages. With the help of toluidine blue and periodic acid-Schiff staining to visualize cell walls and identify their constituents, we observed decreasing cambial cell layers, thickening of newly formed xylem cell walls and increased polysaccharide granules in phloem from June to the following March over the course of our collecting period. Pectin immunofluorescence showed that dormant-stage cambium can produce highly abundant de-esterified homogalacturonan and (1-4)-β-d-galactan epitopes, whereas active cambium can strong accumulate high methylesterified homogalacturonan. Calcofluor white staining and confocal Raman spectroscopy analysis revealed regular changes in the chemical composition of cell walls, such as relative lower cellulose deposition in transition stage in vascular cambium, and higher lignin accumulation was found in dormant stage in secondary xylem. Moreover, real-time quantitative polymerase chain reaction analysis suggested that various IAA (Aux/IAA protein), CesA, CslA and HDZ genes, as well as NAC, PME3 and PME4, may be involved in cambium activities and secondary xylem formation. Taken together, these findings provide new information about cambium activity and cell differentiation in the formation, structure and chemistry in conifers during the active-dormant transition.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

Identification and Characterization of Cinnamyl Alcohol Dehydrogenase Encoding Genes Involved in Lignin Biosynthesis and Resistance to Verticillium dahliae in Upland Cotton (Gossypium hirsutum L.)

Verticillium wilt, caused by the soil-borne fungus Verticillium dahliae, is one of the most devastating diseases in cotton (Gossypium spp.). Lignin in the cell wall forms a physical barrier to inhibit pathogen invasion, and defense-induced lignification reinforces secondary cell wall to prevent pathogens from further spreading. Cinnamyl alcohol dehydrogenases (CADs) catalyze the production of three main monolignols, p-coumaryl- (H), coniferyl- (G), and sinapyl-alcohols (S), which are the fundamental blocks of lignin. Here, we identified CAD genes in G. hirsutum, analyzed their expression profiles in cotton leaf, stem, and root from different developmental stages, and selected GhCAD35, GhCAD45, and GhCAD43, which were consistently induced by V. dahliae inoculation in G. hirsutum cultivars resistant or susceptible to V. dahliae. On the basis of confirmation of the in vitro enzymatic activity of the three proteins in generation of the three monolignols, we used virus-induced gene silencing (VIGS) to investigate the effects of silencing of GhCAD35, GhCAD45, or GhCAD43 on resistance to V. dahliae as well as on deposition and the composition of lignin. Silencing each of the three CADs impaired the defense-induced lignification and salicylic acid biosynthesis in stem, and compromised resistance to V. dahliae. Moreover, our study showed that silencing the three GhCADs severely affected the biosynthesis of S-lignin, leading to a decrease of the syringyl/guaiacyl (S/G) ratio. Heterogeneous overexpression of GhCAD35, GhCAD45, or GhCAD43 in Arabidopsis enhanced disease resistance. Taken together, our study demonstrates a role of the three GhCADs in defense-induced lignin biosynthesis and resistance to V. dahliae in G. hirsutum.

The role of senescence-associated gene101 (PagSAG101a) in the regulation of secondary xylem formation in poplar

Wood is produced by the accumulation of secondary xylem via proliferation and differentiation of the cambium cells in woody plants. Identifying the regulators involved in this process remains a challenging task. In this study, we isolated PagSAG101a, the homolog of Arabidopsis thaliana SAG101, from a hybrid poplar (Populus alba × Populus glandulosa) clone 84K and investigated its role in secondary xylem development. PagSAG101a was expressed predominantly in lignified stems and localized in the nucleus. Compared with non-transgenic 84K plants, transgenic plants overexpressing PagSAG101a displayed increased plant height, internode number, stem diameter, xylem width, and secondary cell wall thickness, while opposite phenotypes were observed for PagSAG101a knock-out plants. Transcriptome analyses revealed that differentially expressed genes were enriched for those controlling cambium cell division activity and subsequent secondary cell wall deposition during xylem formation. In addition, the tandem CCCH zinc finger protein PagC3H17, which positively regulates secondary xylem width and secondary wall thickening in poplar, could bind to the promoter of PagSAG101a and mediate the regulation of xylem differentiation. Our results support that PagSAG101a, downstream of PagC3H17, functions in wood development.

Genome-Wide Identification and Characterization of Caffeic Acid O-Methyltransferase Gene Family in Soybean

Soybean is one of the most important legumes, providing high-quality protein for humans. The caffeic acid O-methyltransferase (COMT) gene has previously been demonstrated to be a critical gene that regulates lignin production in plant cell walls and plays an important function in plant growth and development. However, the COMT gene family has not been studied in soybeans. In this study, 55 COMT family genes in soybean were identified by phylogenetic analysis and divided into two groups, I and II. The analysis of conserved domains showed that all GmCOMTs genes contained Methyltransferase-2 domains. Further prediction of cis-acting elements showed that GmCOMTs genes were associated with growth, light, stress, and hormonal responses. Eventually, based on the genomic data of soybean under different stresses, the results showed that the expression of GmCOMTs genes was different under different stresses, such as salt and drought stress. This study has identified and characterized the COMT gene family in soybean, which provides an important theoretical basis for further research on the biological functions of COMT genes and promotes revealing the role of GmCOMTs genes under stress resistance.

Differential transcriptome analysis and identification of genes related to resistance to blight in three varieties of Bambusa pervariabilis × Dendrocalamopsis grandis

Background

Bambusa pervariabilis × Dendrocalamopsis grandis is a fast-growing bamboo that is widely introduced in southern China and has great economic and ecological benefits. In recent years, a blight of B. pervariabilis × D. grandis caused by Arthrinium phaeospermum has led to much branch damage and even death of entire bamboo forests.

Methods

To screen for resistance genes in B. pervariabilis × D. grandis, transcriptome sequencing technology was used to compare the gene expression profiles of different varieties of B. pervariabilis × D. grandis with variable resistance and the same varieties under different treatments. The Clusters of Orthologous Groups of Proteins (COG) database; the Gene Ontology (GO) database; and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were used to annotate and analyse the differentially expressed genes.

Results

A total of 26,157 and 11,648 differentially expressed genes were obtained in the different varieties after inoculation with A. phaeospermum and the same varieties after inoculation A. phaeospermum or sterile water, respectively. There were 23 co-upregulated DGEs and 143 co-downregulated DEGs in #3 and #8, #6 and #8, #6 and #3. There were 50 co-upregulated DGEs and 24 co-downregulated DEGs in the same varieties after inoculation A. phaeospermum or sterile water. The results showed that many genes involved in cell wall composition synthesis, redox reactions and signal transduction were significantly different after pathogen infection. Twenty-one candidate genes for blight resistance, such as pme53, cad5, pod, gdsl-ll and Myb4l, were found. The qRT-PCR results were consistent with the sequencing results, verifying their authenticity. These results provide a foundation for the further exploration of resistance genes and their functions.

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.

Enzyme Complexes of Ptr4CL and PtrHCT Modulate Co-enzyme A Ligation of Hydroxycinnamic Acids for Monolignol Biosynthesis in Populus trichocarpa

Co-enzyme A (CoA) ligation of hydroxycinnamic acids by 4-coumaric acid:CoA ligase (4CL) is a critical step in the biosynthesis of monolignols. Perturbation of 4CL activity significantly impacts the lignin content of diverse plant species. In Populus trichocarpa, two well-studied xylem-specific Ptr4CLs (Ptr4CL3 and Ptr4CL5) catalyze the CoA ligation of 4-coumaric acid to 4-coumaroyl-CoA and caffeic acid to caffeoyl-CoA. Subsequently, two 4-hydroxycinnamoyl-CoA:shikimic acid hydroxycinnamoyl transferases (PtrHCT1 and PtrHCT6) mediate the conversion of 4-coumaroyl-CoA to caffeoyl-CoA. Here, we show that the CoA ligation of 4-coumaric and caffeic acids is modulated by Ptr4CL/PtrHCT protein complexes. Downregulation of PtrHCTs reduced Ptr4CL activities in the stem-differentiating xylem (SDX) of transgenic P. trichocarpa. The Ptr4CL/PtrHCT interactions were then validated in vivo using biomolecular fluorescence complementation (BiFC) and protein pull-down assays in P. trichocarpa SDX extracts. Enzyme activity assays using recombinant proteins of Ptr4CL and PtrHCT showed elevated CoA ligation activity for Ptr4CL when supplemented with PtrHCT. Numerical analyses based on an evolutionary computation of the CoA ligation activity estimated the stoichiometry of the protein complex to consist of one Ptr4CL and two PtrHCTs, which was experimentally confirmed by chemical cross-linking using SDX plant protein extracts and recombinant proteins. Based on these results, we propose that Ptr4CL/PtrHCT complexes modulate the metabolic flux of CoA ligation for monolignol biosynthesis during wood formation in P. trichocarpa.

Genome-wide analysis of the lignin toolbox for morus and the roles of lignin related genes in response to zinc stress

Mulberry (Morus, Moraceae) is an important economic plant with nutritional, medicinal, and ecological values. Lignin in mulberry can affect the quality of forage and the saccharification efficiency of mulberry twigs. The availability of the Morus notabilis genome makes it possible to perform a systematic analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway, providing the core genes for the lignin toolbox in mulberry. We performed genome-wide screening, which was combined with de novo transcriptome data for Morus notabilis and Morus alba variety Fengchi, to identify putative members of the lignin gene families followed by phylogenetic and expression profile analyses. We focused on bona fide clade genes and their response to zinc stress were further distinguished based on expression profiles using RNA-seq and RT-qPCR. We finally identified 31 bona fide genes in Morus notabilis and 25 bona fide genes in Fengchi. The putative function of these bona fide genes was proposed, and a lignin toolbox that comprised 19 genes in mulberry was provided, which will be convenient for researchers to explore and modify the monolignol biosynthesis pathway in mulberry. We also observed changes in the expression of some of these lignin biosynthetic genes in response to stress caused by excess zinc in Fengchi and proposed that the enhanced lignin biosynthesis in lignified organs and inhibition of lignin biosynthesis in leaf is an important response to zinc stress in mulberry.

Phylogenetic Occurrence of the Phenylpropanoid Pathway and Lignin Biosynthesis in Plants

The phenylpropanoid pathway serves as a rich source of metabolites in plants and provides precursors for lignin biosynthesis. Lignin first appeared in tracheophytes and has been hypothesized to have played pivotal roles in land plant colonization. In this review, we summarize recent progress in defining the lignin biosynthetic pathway in lycophytes, monilophytes, gymnosperms, and angiosperms. In particular, we review the key structural genes involved in p-hydroxyphenyl-, guaiacyl-, and syringyl-lignin biosynthesis across plant taxa and consider and integrate new insights on major transcription factors, such as NACs and MYBs. We also review insight regarding a new transcriptional regulator, 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase, canonically identified as a key enzyme in the shikimate pathway. We use several case studies, including EPSP synthase, to illustrate the evolution processes of gene duplication and neo-functionalization in lignin biosynthesis. This review provides new insights into the genetic engineering of the lignin biosynthetic pathway to overcome biomass recalcitrance in bioenergy crops.

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.

A peroxisomal β-oxidative pathway contributes to the formation of C6-C1 aromatic volatiles in poplar

Benzenoids (C6-C1 aromatic compounds) play important roles in plant defense and are often produced upon herbivory. Black cottonwood (Populus trichocarpa) produces a variety of volatile and nonvolatile benzenoids involved in various defense responses. However, their biosynthesis in poplar is mainly unresolved. We showed feeding of the poplar leaf beetle (Chrysomela populi) on P. trichocarpa leaves led to increased emission of the benzenoid volatiles benzaldehyde, benzylalcohol, and benzyl benzoate. The accumulation of salicinoids, a group of nonvolatile phenolic defense glycosides composed in part of benzenoid units, was hardly affected by beetle herbivory. In planta labeling experiments revealed that volatile and nonvolatile poplar benzenoids are produced from cinnamic acid (C6-C3). The biosynthesis of C6-C1 aromatic compounds from cinnamic acid has been described in petunia (Petunia hybrida) flowers where the pathway includes a peroxisomal-localized chain shortening sequence, involving cinnamate-CoA ligase (CNL), cinnamoyl-CoA hydratase/dehydrogenase (CHD), and 3-ketoacyl-CoA thiolase (KAT). Sequence and phylogenetic analysis enabled the identification of small CNL, CHD, and KAT gene families in P. trichocarpa. Heterologous expression of the candidate genes in Escherichia coli and characterization of purified proteins in vitro revealed enzymatic activities similar to those described in petunia flowers. RNA interference-mediated knockdown of the CNL subfamily in gray poplar (Populus x canescens) resulted in decreased emission of C6-C1 aromatic volatiles upon herbivory, while constitutively accumulating salicinoids were not affected. This indicates the peroxisomal β-oxidative pathway participates in the formation of volatile benzenoids. The chain shortening steps for salicinoids, however, likely employ an alternative pathway.

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.

Genome-wide identification and analysis of monolignol biosynthesis genes in Salix matsudana Koidz and their relationship to accelerated growth

Lignin plays an important role in plant growth and development. It serves as a raw material for the manufacture of paper, animal feed, and chemical fertilizers. However, the regulation of lignin biosynthesis genes and the composition of the relevant gene families remain unclear in many plant species. Here, we identified and characterized 11 families of monolignol biosynthesis genes in Salix matsudana Koidz. Based on phylogenetic analysis of lignin biosynthesis genes from nine angiosperm species (Arabidopsis thaliana, Oryza sativa, Zea mays, Solanum lycopersicum, S. suchowensis, S. purpurea, Populus euphratica, P. trichocarpa, and S. matsudana), the 11 gene families could be divided into two classes that differed in their apparent evolutionary history. We compared the distribution of lignin biosynthesis genes between the two sub-genomes (At and Bt) of S. matsudana and found that more duplicated genes were present in the Bt sub-genome. We analyzed RNA sequencing data from two parents of contrasting height and two of their F1 progeny, and detected 23 differentially expressed genes (DEGs) that may regulate accelerated growth. We analyzed the promoter regions of the lignin-related DEGs and identified several hormone-related (auxin, ethylene, and cytokinin) transcription factor binding sites. These results provide an important foundation for future studies on the molecular mechanisms and genetic regulation of lignin biosynthesis and its relationship to accelerated growth in forest trees.

Genome sequence and evolution of Betula platyphylla

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1-MKK2-MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

Genome-wide analysis of general phenylpropanoid and monolignol-specific metabolism genes in sugarcane

Lignin is the main component of secondary cell walls and is essential for plant development and defense. However, lignin is recognized as a major recalcitrant factor for efficiency of industrial biomass processing. Genes involved in general phenylpropanoid and monolignol-specific metabolism in sugarcane have been previously analyzed at the transcriptomic level. Nevertheless, the number of genes identified in this species is still very low. The recently released sugarcane genome sequence has allowed the genome-wide characterization of the 11 gene families involved in the monolignol biosynthesis branch of the phenylpropanoid pathway. After an exhaustive analysis of sugarcane genomes, 438 haplotypes derived from 175 candidate genes from Saccharum spontaneum and 144 from Saccharum hybrid R570 were identified as associated with this biosynthetic route. The phylogenetic analyses, combined with the search for protein conserved residues involved in the catalytic activity of the encoded enzymes, were employed to identify the family members potentially involved in developmental lignification. Accordingly, 15 candidates were identified as bona fide lignin biosynthesis genes: PTAL1, PAL2, C4H4, 4CL1, HCT1, HCT2, C3'H1, C3'H2, CCoAOMT1, COMT1, F5H1, CCR1, CCR2, CAD2, and CAD7. For this core set of lignin biosynthetic genes, we searched for the chromosomal location, the gene expression pattern, the promoter cis-acting elements, and microRNA targets. Altogether, our results present a comprehensive characterization of sugarcane general phenylpropanoid and monolignol-specific genes, providing the basis for further functional studies focusing on lignin biosynthesis manipulation and biotechnological strategies to improve sugarcane biomass utilization.

Impact of abiotic stress on the regulation of cell wall biosynthesis in Populus trichocarpa

Growth of biomass for lignocellulosic biofuels and biomaterials may take place on land unsuitable for foods, meaning the biomass plants are exposed to increased abiotic stresses. Thus, the understanding how this affects biomass composition and quality is important for downstream bioprocessing. Here, we analyzed the effect of drought and salt stress on cell wall biosynthesis in young shoots and xylem tissues of Populus trichocarpa using transcriptomic and biochemical methods. Following exposure to abiotic stress, stem tissues reduced vessel sizes, and young shoots increased xylem formation. Compositional analyses revealed a reduction in the total amount of cell wall polysaccharides. In contrast, the total lignin amount was unchanged, while the ratio of S/G lignin was significantly decreased in young shoots. Consistent with these observations, transcriptome analyses show that the expression of a subset of cell wall-related genes is tightly regulated by drought and salt stresses. In particular, the expression of a part of genes encoding key enzymes for S-lignin biosynthesis, caffeic acid O-methyltransferase and ferulate 5-hydroxylase, was decreased, suggesting the lower S/G ratio could be partly attributed to the down-regulation of these genes. Together, our data identifies a transcriptional abiotic stress response strategy in poplar, which results in adaptive changes to the plant cell wall.

Transgenic Poplar Designed for Biofuels

Members of the genus Populus (i.e., cottonwood, hybrid poplar) represent a promising source of lignocellulosic biomass for biofuels. However, one of the major factors negatively affecting poplar's efficient conversion to biofuel is the inherent recalcitrance to enzymatic saccharification due to cell wall components such as lignin. To this effect, there have been efforts to modify gene expression to reduce biomass recalcitrance by changing cell wall properties. Here, we review recent genetic modifications of poplar that led to change cell wall properties and the resulting effects on subsequent pretreatment efficacy and saccharification. Although genetic engineering's impacts on cell wall properties are not fully predictable, recent studies have shown promising improvement in the biological conversion of transgenic poplar to biofuels.

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

Background

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

Results

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

Conclusions

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