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.

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.

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 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.

Transcriptomic Profiling of Cryptomeria fortunei Hooibrenk Vascular Cambium Identifies Candidate Genes Involved in Phenylpropanoid Metabolism

Cryptomeria fortunei Hooibrenk (Chinese cedar) is a coniferous tree from southern China that has an important function in landscaping and timber production. Lignin is one of the key components of secondary cell walls, which have a crucial role in conducting water and providing mechanical support for the upward growth of plants. It is mainly biosynthesized via the phenylpropanoid metabolic pathway, of which the molecular mechanism remains so far unresolved in C. fortunei. In order to obtain further insight into this pathway, we performed transcriptome sequencing of the C. fortunei cambial zone at 5 successive growth stages. We generated 78,673 unigenes from transcriptome data, of which 45,214 (57.47%) were successfully annotated in the non-redundant protein database (NR). A total of 8975 unigenes were identified to be significantly differentially expressed between Sample_B and Sample_A after analyzing their expression profiles. Of the differentially expressed genes (DEGs), 6817 (75.96%) and 2158 (24.04%) were up- and down-regulated, respectively. 83 DEGs were involved in phenylpropanoid metabolism, 37 DEGs that encoded v-Myb avian myeloblastosis viral oncogene homolog (MYB) transcription factor (TF), and many candidates that encoded lignin synthesizing enzymes. These findings contribute to understanding the expression pattern of C. fortunei cambial zone transcriptome. Furthermore, our results provide additional insight towards understanding the molecular mechanisms of wood formation in C. fortunei.

ETHYLENE RESPONSE FACTOR39-MYB8 complex regulates low-temperature-induced lignification of loquat fruit

Flesh lignification is a specific chilling response that causes deterioration in the quality of stored red-fleshed loquat fruit (Eribotrya japonica) and is one aspect of wider chilling injury. APETALA2/ETHLENE RESPONSIVE FACTOR (AP2/ERF) transcription factors are important regulators of plant low-temperature responses and lignin biosynthesis. In this study, the expression and action of 27 AP2/ERF genes from the red-fleshed loquat cultivar 'Luoyangqing' were investigated in order to identify transcription factors regulating low-temperature-induced lignification. EjERF27, EjERF30, EjERF36, and EjERF39 were significantly induced by storage at 0 °C but inhibited by a low-temperature conditioning treatment (pre-storage at 5 °C for 6 days before storage at 0 °C, which reduces low-temperature-induced lignification), and their transcript levels positively correlated with flesh lignification. A dual-luciferase assay indicated that EjERF39 could transactivate the promoter of the lignin biosynthetic gene Ej4CL1, and an electrophoretic mobility shift assay confirmed that EjERF39 recognizes the DRE element in the promoter region of Ej4CL1. Furthermore, the combination of EjERF39 and the previously characterized EjMYB8 synergistically transactivated the Ej4CL1 promoter, and both transcription factors showed expression patterns correlated with lignification in postharvest treatments and red-fleshed 'Luoyangqing' and white-fleshed 'Ninghaibai' cultivars with different lignification responses. Bimolecular fluorescence complementation and luciferase complementation imaging assays confirmed direct protein-protein interaction between EjERF39 and EjMYB8. These data indicate that EjERF39 is a novel cold-responsive transcriptional activator of Ej4CL1 that forms a synergistic activator complex with EjMYB8 and contributes to loquat fruit lignification at low temperatures.

Insights of Molecular Mechanism of Xylem Development in Five Black Poplar Cultivars

Black poplar (Populus deltoides, P. nigra, and their hybrids) is the main poplar cultivars in China. It offers interesting options of large-scale biomass production for bioenergy due to its rapid growth and high yield. Poplar wood properties were associated with chemical components and physical structures during wood formation. In this study, five poplar cultivars, P. euramericana 'Zhonglin46' (Pe1), P. euramericana 'Guariento' (Pe2), P. nigra 'N179' (Pn1), P. deltoides 'Danhong' (Pd1), and P. deltoides 'Nanyang' (Pd2), were used to explore the molecular mechanism of xylem development. We analyzed the structural differences of developing xylem in the five cultivars and profiled the transcriptome-wide gene expression patterns through RNA sequencing. The cross sections of the developing xylem showed that the cell wall thickness of developed fiber in Pd1 was thickest and the number of xylem vessels of Pn1 was the least. A total of 10,331 differentially expressed genes were identified among 10 pairwise comparisons of the five cultivars, most of them were related to programmed cell death and secondary cell wall thickening. K-means cluster analysis and Gene Ontology enrichment analysis showed that the genes highly expressed in Pd1 were related to nucleotide decomposition, metabolic process, transferase, and microtubule cytoskeleton; whereas the genes highly expressed in Pn1 were involved in cell wall macromolecule decomposition and polysaccharide binding processes. Based on a weighted gene co-expression network analysis, a large number of candidate regulators for xylem development were identified. And their potential regulatory roles to cell wall biosynthesis genes were validated by a transient overexpression system. This study provides a set of promising candidate regulators for genetic engineering to improve feedstock and enhance biofuel conversion in the bioenergy crop Populus.

Over-expression of transcription factor ARK1 gene leads to down-regulation of lignin synthesis related genes in hybrid poplar '717'

Improving wood growth rate and wood quality are worthy goals in forest genetics and breeding research. The ARK1 gene is one member of the ARBORKNOX family in all plants, which play an essential role in the process of plant growth and development, but the mechanism associated with its gene network regulation is poorly investigated. In order to generate over-expression transgenic hybrid poplar, the agrobacterium-mediated transformation was used to obtain transgenic hybrid poplar ‘717’ plants to provide insight into the function of the ARK1 gene in poplar. Moreover, the morphology of transgenic plants was observed, and transcriptome analysis was performed to explore the ARK1 gene function. The results showed that there were significant differences in pitch, stem diameter, petiole length, leaf width, leaf length and seedling height between ARK1 transgenic seedlings and non-transgenic seedlings. The transgenic seedlings usually had multiple branches and slender leaves, with some leaves not being fully developed. The results of transcriptome analysis showed that the differentially expressed genes were involved in the growth of poplars, including proteins, transcription factors and protein kinases. Genes related to the positive regulation in plant hormone signal transduction pathways were up-regulated, and the genes related to lignin synthesis were down-regulated. The RT-qPCR analysis confirmed the expression levels of the genes involved in the plant hormone signal transduction pathways and phenylpropanoid pathway. In conclusion, the ARK1 gene had a positive regulatory effect on plant growth, and the gene’s coding enzymes related to lignin synthesis were down-regulated.

MicroRNA6443-mediated regulation of FERULATE 5-HYDROXYLASE gene alters lignin composition and enhances saccharification in Populus tomentosa

Ferulate 5-hydroxylase (F5H) is a limiting enzyme involved in biosynthesizing sinapyl (S) monolignol in angiosperms. Genetic regulation of F5H can influence S monolignol synthesis and therefore improve saccharification efficiency and biofuel production. To date, little is known about whether F5H is post-transcriptionally regulated by endogenous microRNAs (miRNAs) in woody plants. Here, we report that a microRNA, miR6443, specifically regulates S lignin biosynthesis during stem development in Populus tomentosa. In situ hybridization showed that miR6443 is preferentially expressed in vascular tissues. We further identified that F5H2 is the direct target of miR6443. Overexpression of miR6443 decreased the transcript level of F5H2 in transgenic plants, resulting in a significant reduction in S lignin content. Conversely, reduced miR6443 expression by short tandem target mimics (STTM) elevated F5H2 transcripts, therefore increasing S lignin composition. Introduction of a miR6443-resistant form of F5H2 into miR6443-overexpression plants restored lignin ectopic composition, supporting that miR6443 specifically regulated S lignin biosynthesis by repressing F5H2 in P. tomentosa. Furthermore, saccharification assays revealed decreased hexose yields by 7.5-24.5% in miR6443-overexpression plants compared with the wild-type control, and increased hexoses yields by 13.2-14.6% in STTM6443-overexpression plants. Collectively, we demonstrate that miR6443 modulates S lignin biosynthesis by specially regulating F5H2 in P. tomentosa.

Characterization of the Gh4CL gene family reveals a role of Gh4CL7 in drought tolerance

BACKGROUND

The function of 4-coumarate-CoA ligases (4CL) under abiotic stresses has been studied in plants, however, limited is known about the 4CL genes in cotton (G. hirsutum L.) and their roles in response to drought stress.

RESULTS

We performed genome-wide identification of the 4CL genes in G. hirsutum and investigated the expression profiles of the identified genes in various cotton tissues and in response to stress conditions with an aim to identify 4CL gene(s) associated with drought tolerance. We identified 34 putative 4CL genes in G. hirsutum that were clustered into three classes. Genes of the same class usually share a similar gene structure and motif composition. Many cis-elements related to stress and phytohormone responses were found in the promoters of the Gh4CL genes. Of the 34 Gh4CL genes, 26 were induced by at least one abiotic stress and 10 (including Gh4CL7) were up-regulated under the polyethylene glycol (PEG) simulated drought stress conditions. Virus-induced gene silencing (VIGS) in cotton and overexpression (OE) in Arabidopsis thaliana were applied to investigate the biological function of Gh4CL7 in drought tolerance. The Gh4CL7-silencing cotton plants showed more sensitive to drought stress, probably due to decreased relative water content (RWC), chlorophyll content and antioxidative enzyme activity, increased stomatal aperture, and the contents of malondialdehyde (MDA) and hydrogen peroxide (H₂O₂). Arabidopsis lines overexpressing Gh4CL7, however, were more tolerant to drought treatment, which was associated with improved antioxidative enzyme activity, reduced accumulation of MDA and H₂O₂ and up-regulated stress-related genes under the drought stress conditions. In addition, compared to their respective controls, the Gh4CL7-silencing cotton plants and the Gh4CL7-overexpressing Arabidopsis lines had a ~ 20% reduction and a ~ 10% increase in lignin content, respectively. The expression levels of genes related to lignin biosynthesis, including PAL, CCoAOMT, COMT, CCR and CAD, were lower in Gh4CL7-silencing plants than in controls. Taken together, these results demonstrated that Gh4CL7 could positively respond to drought stress and therefore might be a candidate gene for improvement of drought tolerance in cotton.

CONCLUSION

We characterized the 4CL gene family in upland cotton and revealed a role of Gh4CL7 in lignin biosynthesis and drought tolerance.

Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes

The phenylpropanoid pathway is an important route of secondary metabolism involved in the synthesis of different phenolic compounds such as phenylpropenes, anthocyanins, stilbenoids, flavonoids, and monolignols. The flux toward monolignol biosynthesis through the phenylpropanoid pathway is controlled by specific genes from at least ten families. Lignin polymer is one of the major components of the plant cell wall and is mainly responsible for recalcitrance to saccharification in ethanol production from lignocellulosic biomass. Here, we identified and characterized sugarcane candidate genes from the general phenylpropanoid and monolignol-specific metabolism through a search of the sugarcane EST databases, phylogenetic analysis, a search for conserved amino acid residues important for enzymatic function, and analysis of expression patterns during culm development in two lignin-contrasting genotypes. Of these genes, 15 were cloned and, when available, their loci were identified using the recently released sugarcane genomes from Saccharum hybrid R570 and Saccharum spontaneum cultivars. Our analysis points out that ShPAL1, ShPAL2, ShC4H4, Sh4CL1, ShHCT1, ShC3H1, ShC3H2, ShCCoAOMT1, ShCOMT1, ShF5H1, ShCCR1, ShCAD2, and ShCAD7 are strong candidates to be bona fide lignin biosynthesis genes. Together, the results provide information about the candidate genes involved in monolignol biosynthesis in sugarcane and may provide useful information for further molecular genetic studies in sugarcane.

Biochemical characterization and expression analysis of lignification in two pear (Pyrus ussuriensis Maxim.) varieties with contrasting stone cell content

As lignified stone cells reduce fruit quality, we investigated lignin deposition, phenolic metabolites, and expression of lignin biosynthetic genes during fruit development to elucidate the molecular mechanism of stone cell lignification using histological, biochemical, and transcriptional data from two Ussurian pear varieties (Jianba and Nanguo) with contrasting stone cell content. Lignin content and distribution coincided with stone cell accumulation. As per LC-MS analysis, Jianba exhibited higher levels of lignin monomers and hydroxycinnamates than Nanguo, consistently with lignin amount in each case. However, flavonoid content was much higher in Nanguo. Transcriptional data showed that most monolignol biosynthesis-related genes were particularly upregulated in Jianba during lignin accumulation; especially CCR and LAC, two monolignol biosynthesis-specific genes, were substantially upregulated in Jianba fruits at critical stages. Therefore, differences in stone cell content between "Jianba" and "Nanguo" may result from differential expression of lignin synthase genes located downstream of the lignin biosynthesis pathway. Taken together, our data may provide a deeper understanding of the molecular mechanism for stone cell lignification in pear fruit.

Functional Characteristics of Caffeoyl Shikimate Esterase in Larix Kaempferi and Monolignol Biosynthesis in Gymnosperms

Caffeoyl shikimate esterase (CSE) has been reported to be involved in lignin biosynthesis; however, studies of CSE in gymnosperms are lacking. In this study, CSE was successfully cloned from Larix kaempferi (LkCSE) based on Larix laricina transcriptome screening. LkCSE was likely to have catalytic activity based on homologous sequence alignment and phylogenetic analyses of CSEs from different species. In vitro assays with the recombinant enzyme validated the catalytic activity of LkCSE, indicating its function in converting caffeoyl shikimate into caffeate and shikimate. Additionally, the optimum reaction pH and temperature of LkCSE were determined to be 6.0 and 30 °C, respectively. The values of K_m and V_max of CSE for caffeoyl shikimate were 98.11 μM and 14.44 nM min^-1, respectively. Moreover, LkCSE was observed to have tissue expression specificity and was abundantly expressed in stems and leaves, especially stems, which was 50 times higher than the expression levels of roots. Lastly, translational fusion assays using LkCSE fused with green fluorescent proteins (GFP) in tobacco leaves indicated that LkCSE was localized in the plasma membrane and endoplasmic reticulum (ER). These results revealed that CSE clearly functions in gymnosperms and it is possible for LkCSE to interact with other ER-resident proteins and regulate mass flux in the monolignol biosynthesis pathway.