Genome-wide characterization of the lignification toolbox in Arabidopsis

Characterization and Expression Profiling of Camellia sinensis Cinnamate 4-hydroxylase Genes in Phenylpropanoid Pathways

Cinnamate 4-hydroxylase (C4H), a cytochrome P450-dependent monooxygenase, participates in the synthesis of numerous polyphenoid compounds, such as flavonoids and lignins. However, the C4H gene number and function in tea plants are not clear. We screened all available transcriptome and genome databases of tea plants and three C4H genes were identified and named CsC4Ha, CsC4Hb, and CsC4Hc, respectively. Both CsC4Ha and CsC4Hb have 1518-bp open reading frames that encode 505-amino acid proteins. CsC4Hc has a 1635-bp open reading frame that encodes a 544-amino acid protein. Enzymatic analysis of recombinant proteins expressed in yeast showed that the three enzymes catalyzed the formation of p-coumaric acid (4-hydroxy trans-cinnamic acid) from trans-cinnamic acid. Quantitative real-time PCR (qRT-PCR) analysis showed that CsC4Ha was highly expressed in the 4th leaf, CsC4Hb was highly expressed in tender leaves, while CsC4Hc was highly expressed in the young stems. The three CsC4Hs were induced with varying degrees by abiotic stress treatments. These results suggest they may have different subcellular localization and different physiological functions.

A gene encoding maize caffeoyl-CoA O-methyltransferase confers quantitative resistance to multiple pathogens

Alleles that confer multiple disease resistance (MDR) are valuable in crop improvement, although the molecular mechanisms underlying their functions remain largely unknown. A quantitative trait locus, qMdr_9.02, associated with resistance to three important foliar maize diseases-southern leaf blight, gray leaf spot and northern leaf blight-has been identified on maize chromosome 9. Through fine-mapping, association analysis, expression analysis, insertional mutagenesis and transgenic validation, we demonstrate that ZmCCoAOMT2, which encodes a caffeoyl-CoA O-methyltransferase associated with the phenylpropanoid pathway and lignin production, is the gene within qMdr_9.02 conferring quantitative resistance to both southern leaf blight and gray leaf spot. We suggest that resistance might be caused by allelic variation at the level of both gene expression and amino acid sequence, thus resulting in differences in levels of lignin and other metabolites of the phenylpropanoid pathway and regulation of programmed cell death.

Spatial regulation of monolignol biosynthesis and laccase genes control developmental and stress-related lignin in flax

BACKGROUND

Bast fibres are characterized by very thick secondary cell walls containing high amounts of cellulose and low lignin contents in contrast to the heavily lignified cell walls typically found in the xylem tissues. To improve the quality of the fiber-based products in the future, a thorough understanding of the main cell wall polymer biosynthetic pathways is required. In this study we have carried out a characterization of the genes involved in lignin biosynthesis in flax along with some of their regulation mechanisms.

RESULTS

We have first identified the members of the phenylpropanoid gene families through a combination of in silico approaches. The more specific lignin genes were further characterized by high throughput transcriptomic approaches in different organs and physiological conditions and their cell/tissue expression was localized in the stems, roots and leaves. Laccases play an important role in the polymerization of monolignols. This multigenic family was determined and a miRNA was identified to play a role in the posttranscriptional regulation by cleaving the transcripts of some specific genes shown to be expressed in lignified tissues. In situ hybridization also showed that the miRNA precursor was expressed in the young xylem cells located near the vascular cambium. The results obtained in this work also allowed us to determine that most of the genes involved in lignin biosynthesis are included in a unique co-expression cluster and that MYB transcription factors are potentially good candidates for regulating these genes.

CONCLUSIONS

Target engineering of cell walls to improve plant product quality requires good knowledge of the genes responsible for the production of the main polymers. For bast fiber plants such as flax, it is important to target the correct genes from the beginning since the difficulty to produce transgenic material does not make possible to test a large number of genes. Our work determined which of these genes could be potentially modified and showed that it was possible to target different regulatory pathways to modify lignification.

The effect of drought stress on the expression of key genes involved in the biosynthesis of phenylpropanoids and essential oil components in basil (Ocimum basilicum L.)

Basil (Ocimum basilicum L.), a medicinal plant of the Lamiaceae family, is used in traditional medicine; its essential oil is a rich source of phenylpropanoids. Methylchavicol and methyleugenol are the most important constituents of basil essential oil. Drought stress is proposed to enhance the essential oil composition and expression levels of the genes involved in its biosynthesis. In the current investigation, an experiment based on a completely randomized design (CRD) with three replications was conducted in the greenhouse to study the effect of drought stress on the expression level of four genes involved in the phenylpropanoid biosynthesis pathway in O. basilicum c.v. Keshkeni luvelou. The genes studied were chavicol O-methyl transferase (CVOMT), eugenol O-methyl transferase (EOMT), cinnamate 4-hydroxylase (C4H), 4-coumarate coA ligase (4CL), and cinnamyl alcohol dehydrogenase (CAD). The effect of drought stress on the essential oil compounds and their relationship with the expression levels of the studied genes were also investigated. Plants were subjected to levels of 100%, 75%, and 50% of field capacity (FC) at the 6-8 leaf stage. Essential oil compounds were identified by gas chromatography/mass spectrometry (GC-MS) at flowering stage and the levels of gene expression were determind by real time PCR in plant leaves at the same stage. Results showed that drought stress increased the amount of methylchavicol, methyleugenol, β-Myrcene and α-bergamotene. The maximum amount of these compounds was observed at 50% FC. Real-time PCR analysis revealed that severe drought stress (50% FC) increased the expression level of CVOMT and EOMT by about 6.46 and 46.33 times, respectively, whereas those of CAD relatively remained unchanged. The expression level of 4CL and C4H reduced under drought stress conditions. Our results also demonstrated that changes in the expression levels of CVOMT and EOMT are significantly correlated with methylchavicol (r = 0.94, P ≤ 0.05) and methyleugenol (r = 0.98, P ≤ 0.05) content. Thus, drought stress probably increases the methylchavicol and methyleugenol content, in part, through increasing the expression levels of CVOMT and EOMT.

Fine-tuning of the flavonoid and monolignol pathways during apple early fruit development

A coordinated regulation of different branches of the flavonoid pathway was highlighted that may contribute to elucidate the role of this important class of compounds during the early stages of apple fruit development. Apple (Malus × domestica Borkh.) is an economically important fruit appreciated for its organoleptic characteristics and its benefits for human health. The first stages after fruit set represent a very important and still poorly characterized developmental process. To enable the profiling of genes involved in apple early fruit development, we combined the suppression subtractive hybridization (SSH) protocol to next-generation sequencing. We identified and characterized genes induced and repressed during fruit development in the apple cultivar 'Golden Delicious'. Our results showed an opposite regulation of genes coding for enzymes belonging to flavonoid and monolignol pathways, with a strong induction of the former and a simultaneous repression of the latter. Two isoforms of phenylalanine ammonia-lyase and 4-coumarate:CoA ligase, key enzymes located at the branching point between flavonoid and monolignol pathways, showed opposite expression patterns during the period in analysis, suggesting a possible regulation mechanism. A targeted metabolomic analysis supported the SSH results and revealed an accumulation of the monomers catechin and epicatechin as well as several forms of procyanidin oligomers in apple fruitlets starting early after anthesis, together with a decreased production of other classes of flavonoids such as some flavonols and the dihydrochalcone phlorizin. Moreover, gene expression and metabolites accumulation of 'Golden Delicious' were compared to a wild apple genotype of Manchurian crabapple (Malus mandshurica (Maxim.) Kom.). Significant differences in both gene expression and metabolites accumulation were found between the two genotypes.

How Embryophytic is the Biosynthesis of Phenylpropanoids and their Derivatives in Streptophyte Algae?

The origin of land plants from algae is a long-standing question in evolutionary biology. It is becoming increasingly clear that many characters that were once assumed to be 'embryophyte specific' can in fact be found in their closest algal relatives, the streptophyte algae. One such case is the phenylpropanoid pathway. While biochemical data indicate that streptophyte algae harbor lignin-like components, the phenylpropanoid core pathway, which serves as the backbone of lignin biosynthesis, has been proposed to have arisen at the base of the land plants. Here we revisit this hypothesis using a wealth of new sequence data from streptophyte algae. Tracing the biochemical pathway towards lignin biogenesis, we show that most of the genes required for phenylpropanoid synthesis and the precursors for lignin production were already present in streptophyte algae. Nevertheless, phylogenetic analyses and protein structure predictions of one of the key enzyme classes in lignin production, cinnamyl alcohol dehydrogenase (CAD), suggest that CADs of streptophyte algae are more similar to sinapyl alcohol dehydrogenases (SADs). This suggests that the end-products of the pathway leading to lignin biosynthesis in streptophyte algae may facilitate the production of lignin-like compounds and defense molecules. We hypothesize that streptophyte algae already possessed the genetic toolkit from which the capacity to produce lignin later evolved in vascular plants.

Functional Characterization of Populus PsnSHN2 in Coordinated Regulation of Secondary Wall Components in Tobacco

Wood formation is a biological process during which the most abundant lignocellulosic biomass on earth is produced. Although a number of transcription factors have been linked to the regulation of wood formation process, none of them has been demonstrated to be a higher hierarchical regulator that coordinately regulates secondary wall biosynthesis genes. Here, we identified a Populus gene, PsnSHN2, a counterpart of the Arabidopsis AP2/ERF type transcription factor, SHINE2. PsnSHN2 is predominantly expressed in xylem tissues and acted evidently as a high hierarchical transcriptional activator. Overexpression of PsnSHN2 in tobacco significantly altered the expression of both transcription factors and biosynthesis genes involved in secondary wall formation, leading to the thickened secondary walls and the changed cell wall composition. The most significant changes occurred in the contents of cellulose and hemicellulose that increased 37% and 28%, respectively, whereas the content of lignin that decreased 34%. Furthermore, PsnSHN2 activated or repressed the promoter activities of transcription factors involved in secondary wall biosynthesis and bound to five cis-acting elements enriched in the promoter regions of these transcription factors. Taken together, our results suggest PsnSHN2 coordinately regulate secondary wall formation through selective up/down-regulation of its downstream transcription factors that control secondary wall formation.

Solanum tuberosum StCDPK1 is regulated by miR390 at the posttranscriptional level and phosphorylates the auxin efflux carrier StPIN4 in vitro, a potential downstream target in potato development

Among many factors that regulate potato tuberization, calcium and calcium-dependent protein kinases (CDPKs) play an important role. CDPK activity increases at the onset of tuber formation with StCDPK1 expression being strongly induced in swollen stolons. However, not much is known about the transcriptional and posttranscriptional regulation of StCDPK1 or its downstream targets in potato development. To elucidate further, we analyzed its expression in different tissues and stages of the life cycle. Histochemical analysis of StCDPK1::GUS (β-glucuronidase) plants demonstrated that StCDPK1 is strongly associated with the vascular system in stems, roots, during stolon to tuber transition, and in tuber sprouts. In agreement with the observed GUS profile, we found specific cis-acting elements in StCDPK1 promoter. In silico analysis predicted miR390 to be a putative posttranscriptional regulator of StCDPK1. Quantitative real time-polymerase chain reaction (qRT-PCR) analysis showed ubiquitous expression of StCDPK1 in different tissues which correlated well with Western blot data except in leaves. On the contrary, miR390 expression exhibited an inverse pattern in leaves and tuber eyes suggesting a possible regulation of StCDPK1 by miR390. This was further confirmed by Agrobacterium co-infiltration assays. In addition, in vitro assays showed that recombinant StCDPK1-6xHis was able to phosphorylate the hydrophilic loop of the auxin efflux carrier StPIN4. Altogether, these results indicate that StCDPK1 expression is varied in a tissue-specific manner having significant expression in vasculature and in tuber eyes; is regulated by miR390 at posttranscriptional level and suggest that StPIN4 could be one of its downstream targets revealing the overall role of this kinase in potato development.

Characterization of the cinnamoyl-CoA reductase (CCR) gene family in Populus tomentosa reveals the enzymatic active sites and evolution of CCR

Two distinct cinnamoyl-coenzyme A reductases (CCRs) from Populus tomentosa were cloned and studied and active sites in CCRs were further identified based on sequence divergence, molecular simulation, and site-directed mutants. Cinnamoyl-coenzyme A (CoA) reductase (CCR) is the first committed gene in the lignin-specific pathway and plays a role in the lignin biosynthesis pathway. In this study, we cloned 11 genes encoding CCR or CCR-like proteins in Populus tomentosa. An enzymatic assay of the purified recombinant P. tomentosa (Pto) CCR and PtoCCR-like proteins indicated that only PtoCCR1 and PtoCCR7 had detectable activities toward hydroxycinnamoyl-CoA esters. PtoCCR1 exhibited specificity for feruloyl-CoA, with no detectable activity for any other hydroxycinnamoyl-CoA esters. However, PtoCCR7 catalyzed p-coumaroyl-CoA, caffeoyl-CoA, feruloyl-CoA, and sinapoyl-CoA with a preference for feruloyl-CoA. Site-directed mutations of selected amino acids divergent between PtoCCR1 and 7, combined with modeling and docking, showed that A132 in CCR7 combined with the catalytic triad might comprise the catalytic center. In CCR7, L192, F155, and H208 were identified as the substrate-binding sites, and site-directed mutations of these amino acids showed obvious changes in catalytic efficiency with respect to both feruloyl-CoA and sinapoyl-CoA. Mutant F155Y exhibited greater catalytic efficiency for sinapoyl-CoA compared with that of wild-type PtoCCR7. Finally, recent genome duplication events provided the foundation for CCR divergence. This study further identified the active sites in CCRs and the evolutionary process of CCRs in terrestrial plants.

Impact of Postharvest Nitric Oxide Treatment on Lignin Biosynthesis-Related Genes in Wax Apple (Syzygium samarangense) Fruit

The role of nitric oxide (NO) during storage in wax apple through NO (10 μL/L) fumigate fruit was investigated. Wax apple fruit treated with NO had a significantly lower rate of weight loss, a softening index, and loss of firmness during storage. The transcriptional profile of 10 genes involved in lignin biosynthesis has been analyzed using quantitative real-time polymerase chain reaction (qRT-PCR). The qRT-PCR analysis showed nine genes regulated in the wax apple (p < 0.05) upon NO fumigation, which coincided with the enzyme activity results (NO group lower than control group in peroxidase, phenylalanine ammonia-lyase, and 4-coumarate-CoA ligase), whose total lignin content decreased upon treatment with NO. These results indicate that NO treatment can effectively delay the softening and senescence of wax apple fruit and play an important regulatory role in lignin biosynthesis.

Regulation of plant secondary metabolism and associated specialized cell development by MYBs and bHLHs

Plants are unrivaled in the natural world in both the number and complexity of secondary metabolites they produce, and the ubiquitous phenylpropanoids and the lineage-specific glucosinolates represent two such large and chemically diverse groups. Advances in genome-enabled biochemistry and metabolomic technologies have greatly increased the understanding of their metabolic networks in diverse plant species. There also has been some progress in elucidating the gene regulatory networks that are key to their synthesis, accumulation and function. This review highlights what is currently known about the gene regulatory networks and the stable sub-networks of transcription factors at their cores that regulate the production of these plant secondary metabolites and the differentiation of specialized cell types that are equally important to their defensive function. Remarkably, some of these core components are evolutionarily conserved between secondary metabolism and specialized cell development and across distantly related plant species. These findings suggest that the more ancient gene regulatory networks for the differentiation of fundamental cell types may have been recruited and remodeled for the generation of the vast majority of plant secondary metabolites and their specialized tissues.

Systematic Analysis of the 4-Coumarate:Coenzyme A Ligase (4CL) Related Genes and Expression Profiling during Fruit Development in the Chinese Pear

In plants, 4-coumarate:coenzyme A ligases (4CLs), comprising some of the adenylate-forming enzymes, are key enzymes involved in regulating lignin metabolism and the biosynthesis of flavonoids and other secondary metabolites. Although several 4CL-related proteins were shown to play roles in secondary metabolism, no comprehensive study on 4CL-related genes in the pear and other Rosaceae species has been reported. In this study, we identified 4CL-related genes in the apple, peach, yangmei, and pear genomes using DNATOOLS software and inferred their evolutionary relationships using phylogenetic analysis, collinearity analysis, conserved motif analysis, and structure analysis. A total of 149 4CL-related genes in four Rosaceous species (pear, apple, peach, and yangmei) were identified, with 30 members in the pear. We explored the functions of several 4CL and acyl-coenzyme A synthetase (ACS) genes during the development of pear fruit by quantitative real-time PCR (qRT-PCR). We found that duplication events had occurred in the 30 4CL-related genes in the pear. These duplicated 4CL-related genes are distributed unevenly across all pear chromosomes except chromosomes 4, 8, 11, and 12. The results of this study provide a basis for further investigation of both the functions and evolutionary history of 4CL-related genes.

Transcriptome Analysis of Stem and Globally Comparison with Other Tissues in Brassica napus

Brassica napus is one of the most important oilseed crops in the world. However, there is currently no enough stem transcriptome information and comparative transcriptome analysis of different tissues, which impedes further functional genomics research on B. napus. In this study, the stem transcriptome of B. napus was characterized by RNA-seq technology. Approximately 13.4 Gb high-quality clean reads with an average length of 100 bp were generated and used for comparative transcriptome analysis with the existing transcriptome sequencing data of roots, leaves, flower buds, and immature embryos of B. napus. All the transcripts were annotated against GO and KEGG databases. The common genes in five tissues, differentially expressed genes (DEGs) of the common genes between stems and other tissues, and tissue-specific genes were detected, and the main biochemical activities and pathways implying the common genes, DEGs and tissue-specific genes were investigated. Accordingly, the common transcription factors (TFs) in the five tissues and tissue-specific TFs were identified, and a TFs-based regulation network between TFs and the target genes involved in 'Phenylpropanoid biosynthesis' pathway were constructed to show several important TFs and key nodes in the regulation process. Collectively, this study not only provided an available stem transcriptome resource in B. napus, but also revealed valuable comparative transcriptome information of five tissues of B. napus for future investigation on specific processes, functions and pathways.

Transcriptome-wide identification and characterization of CAD isoforms specific for podophyllotoxin biosynthesis from Podophyllum hexandrum

Podophyllotoxin (ptox) is a therapeutically important lignan derived from Podophyllum hexandrum and is used as a precursor for the synthesis of anticancer drugs etoposide, teniposide and etopophose. In spite of its enormous economic significance, genomic information on this endangered medicinal herb is scarce. We have performed de novo transcriptome analysis of methyl jasmonate (MeJA)-treated P. hexandrum cell cultures exhibiting enhanced ptox accumulation. The results revealed the maximum up-regulation of several isoforms of cinnamyl alcohol dehydrogenase (CAD). CAD catalyzes the synthesis of coniferyl alcohol and sinapyl alcohol from coniferaldehyde (CAld) and sinapaldehyde respectively. Coniferyl alcohol can produce both lignin and lignan while sinapyl alcohol produces only lignin. To isolate the CAD isoforms favoring ptox, we deduced full length cDNA sequences of four CAD isoforms: PhCAD1, PhCAD2, PhCAD3 and PhCAD4 from the contigs of the transcriptome data. In vitro enzyme assays indicated a higher affinity for CAld over sinapaldehyde for each isoform. In silico molecular docking analyses also suggested that PhCAD3 has a higher binding preference with CAld over sinapaldehyde, followed by PhCAD4, PhCAD2, and PhCAD1, respectively. The transgenic cell cultures overexpressing these isoforms independently revealed that PhCAD3 favored the maximum accumulation of ptox as compared to lignin followed by PhCAD4 and PhCAD2, whereas, PhCAD1 favored both equally. Together, our study reveals transcriptome-wide identification and characterization of ptox specific CAD isoforms from P. hexandrum. It provides a useful resource for future research not only on the ptox biosynthetic pathway but on overall P. hexandrum, an endangered medicinal herb with immense therapeutic importance.

Enzymatic characterization of two acetyl-CoA synthetase genes from Populus trichocarpa

The acetyl-CoA synthetase (ACS) family is a subfamily of adenylate-forming enzymes, which has a close evolutionary relationship with the 4-coumarate:CoA ligase (4CL) family. In this study, two ACS genes were cloned from Populus trichocarpa and were named PtrACS1 and PtrACS2. Bioinformatics characterization of PtrACS1 and PtrACS2 showed that they contained the key ACS residues and a putative peroxisome targeting sequence 1 (PTS1) at the end of the C-terminal sequence. Real-time PCR results showed that PtrACS1 and PtrACS2 were expressed in the phloem, xylem, leaves, and roots of one-year-old P. trichocarpa, but were expressed primarily in the leaves. The ACS enzyme activity was higher in leaves than other tissues in P. trichocarpa. Two overexpressed recombinant proteins showed no catalytic activity toward the substrates of 4CL, but did have notable catalytic activity toward sodium acetate and substrates of ACS. The relative activities of PtrACS1 and PtrACS2 were 194.16 ± 11.23 and 422.25 ± 21.69 μM min(-1) mg(-1), respectively. The K m and V max of PtrACS1 were 0.25 mM and 698.85 μM min(-1) mg(-1), while those for PtrACS2 were 0.72 mM and 245.96 μM min(-1) mg(-1), respectively. Our results revealed that both proteins belong to the ACS family, and provide a theoretical foundation for the identification and functional analysis of members of the adenylate-forming enzyme superfamily.

Evaluation of a novel promoter from Populus trichocarpa for mature xylem tissue specific gene delivery

Wood (i.e., secondary xylem) is an important raw material for many industrial applications. Mature xylem (MX) tissue-specific genetic modification offers an effective means to improve the chemical and physical properties of the wood. Here, we describe a promoter that drives strong gene expression in a MX tissue-specific manner. Using whole-transcriptome genechip analyses of different tissue types of poplar, we identified five candidate genes that had strong expression in the MX tissue. The putative promoter sequences of the five MX-specific genes were evaluated for their promoter activity in both transgenic Arabidopsis and poplar. Among them, we found the promoter of Potri.013G007900.1 (called the PtrMX3 promoter) had the strongest activity in MX and thus was further characterized. In the stem and root tissues of transgenic Arabidopsis plants, the PtrMX3 promoter activity was found exclusively in MX tissue. MX-specific activity of the promoter was reproduced in the stem tissue of transgenic poplar plants. The PtrMX3 promoter activity was not influenced by abiotic stresses or exogenously applied growth regulators, indicating the PtrMX3 promoter is bona fide MX tissue-specific. Our study provides a strong MX-specific promoter for MX-specific modifications of woody biomass.

Isolation and characterization of the 4-coumarate:coenzyme A ligase (4CL1) promoter from Eucalyptus camaldulensis

The most important enzyme of the phenylpropanoid pathway, 4-coumarate:coenzyme A ligase (4CL), is encoded by several homologous genes including 4CL1. The 4CL1 promoter is a tissue-specific gene expression element, particularly active in the secondary xylem or older stems. In this study, the 1127 bp 5'- upstream region of the 4CL1 coding sequence from Eucalyptus camaldulensis, Euc4CL1, was isolated and characterized. Essential putative cis-elements in the Euc4CL1 promoter included: a TATA-box at -22/-28 position, two CCAAT-boxes at -256/-260 and -277/-281 positions, respectively, an AC-element at -328/-336 and A-boxes at -115/-120 and -990/-995 positions. To investigate the effect of the Euc4CL1 promoter on gene expression, a plant transformation vector, pEuc4CL1p, containing the reporter gene for β-glucuronidase (GUS) under the control of Euc4CL1 promoter was constructed based on the pBI101 backbone and introduced in tobacco plants. Stable expression of the GUS gene in transgenic lines was analysed by a histochemical GUS assay. The results indicated the specific expression of the GUS gene in the stem xylem cells of transgenic tobacco lines was controlled by the Euc4CL1 promoter. The observations suggest the isolated Euc4CL1 promoter is a potential candidate for driving the expression of a foreign gene in plant xylem tissues.

Quantitative proteomics reveals protein profiles underlying major transitions in aspen wood development

Background

Wood development is of outstanding interest both to basic research and industry due to the associated cellulose and lignin biomass production. Efforts to elucidate wood formation (which is essential for numerous aspects of both pure and applied plant science) have been made using transcriptomic analyses and/or low-resolution sampling. However, transcriptomic data do not correlate perfectly with levels of expressed proteins due to effects of post-translational modifications and variations in turnover rates. In addition, high-resolution analysis is needed to characterize key transitions. In order to identify protein profiles across the developmental region of wood formation, an in-depth and tissue specific sampling was performed.

Results

We examined protein profiles, using an ultra-performance liquid chromatography/quadrupole time of flight mass spectrometry system, in high-resolution tangential sections spanning all wood development zones in Populus tremula from undifferentiated cambium to mature phloem and xylem, including cell expansion and cell death zones. In total, we analyzed 482 sections, 20–160 μm thick, from four 47-year-old trees growing wild in Sweden. We obtained high quality expression profiles for 3,082 proteins exhibiting consistency across the replicates, considering that the trees were growing in an uncontrolled environment. A combination of Principal Component Analysis (PCA), Orthogonal Projections to Latent Structures (OPLS) modeling and an enhanced stepwise linear modeling approach identified several major transitions in global protein expression profiles, pinpointing (for example) locations of the cambial division leading to phloem and xylem cells, and secondary cell wall formation zones. We also identified key proteins and associated pathways underlying these developmental landmarks. For example, many of the lignocellulosic related proteins were upregulated in the expansion to the early developmental xylem zone, and for laccases with a rapid decrease in early xylem zones. We observed upregulation of two forms of xylem cysteine protease (Potri.002G005700.1 and Potri.005G256000.2; Pt-XCP2.1) in early xylem and their downregulation in late maturing xylem. Our data also show that Pt-KOR1.3 (Potri.003G151700.2) exhibits an expression pattern that supports the hypothesis put forward in previous studies that this is a key xyloglucanase involved in cellulose biosynthesis in primary cell walls and reduction of cellulose crystallinity in secondary walls.

Conclusion

Our novel multivariate approach highlights important processes and provides confirmatory insights into the molecular foundations of wood development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2458-z) contains supplementary material, which is available to authorized users.

CCoAOMT Down-Regulation Activates Anthocyanin Biosynthesis in Petunia

Anthocyanins and volatile phenylpropenes (isoeugenol and eugenol) in petunia (Petunia hybrida) flowers have the precursor 4-coumaryl coenzyme A (CoA) in common. These phenolics are produced at different stages during flower development. Anthocyanins are synthesized during early stages of flower development and sequestered in vacuoles during the lifespan of the flowers. The production of isoeugenol and eugenol starts when flowers open and peaks after anthesis. To elucidate additional biochemical steps toward (iso)eugenol production, we cloned and characterized a caffeoyl-coenzyme A O-methyltransferase (PhCCoAOMT1) from the petals of the fragrant petunia 'Mitchell'. Recombinant PhCCoAOMT1 indeed catalyzed the methylation of caffeoyl-CoA to produce feruloyl CoA. Silencing of PhCCoAOMT1 resulted in a reduction of eugenol production but not of isoeugenol. Unexpectedly, the transgenic plants had purple-colored leaves and pink flowers, despite the fact that cv Mitchell lacks the functional R2R3-MYB master regulator ANTHOCYANIN2 and has normally white flowers. Our results indicate that down-regulation of PhCCoAOMT1 activated the anthocyanin pathway through the R2R3-MYBs PURPLE HAZE (PHZ) and DEEP PURPLE, with predominantly petunidin accumulating. Feeding cv Mitchell flowers with caffeic acid induced PHZ expression, suggesting that the metabolic perturbation of the phenylpropanoid pathway underlies the activation of the anthocyanin pathway. Our results demonstrate a role for PhCCoAOMT1 in phenylpropene production and reveal a link between PhCCoAOMT1 and anthocyanin production.

Lignin biosynthesis in wheat (Triticum aestivum L.): its response to waterlogging and association with hormonal levels

BACKGROUND

Lignin is an important structural component of plant cell wall that confers mechanical strength and tolerance against biotic and abiotic stressors; however it affects the use of biomass such as wheat straw for some industrial applications such as biofuel production. Genetic alteration of lignin quantity and quality has been considered as a viable option to overcome this problem. However, the molecular mechanisms underlying lignin formation in wheat biomass has not been studied. Combining molecular and biochemical approaches, the present study investigated the transcriptional regulation of lignin biosynthesis in two wheat cultivars with varying lodging characteristics and also in response to waterlogging. It also examined the association of lignin level in tissues with that of plant hormones implicated in the control of lignin biosynthesis.

RESULTS

Analysis of lignin biosynthesis in the two wheat cultivars revealed a close association of lodging resistance with internode lignin content and expression of 4-coumarate:CoA ligase1 (4CL1), p-coumarate 3-hydroxylase1 (C3H1), cinnamoyl-CoA reductase2 (CCR2), ferulate 5-hydroxylase2 (F5H2) and caffeic acid O-methyltransferase2 (COMT2), which are among the genes highly expressed in wheat tissues, implying the importance of these genes in mediating lignin deposition in wheat stem. Waterlogging of wheat plants reduced internode lignin content, and this effect is accompanied by transcriptional repression of three of the genes characterized as highly expressed in wheat internode including phenylalanine ammonia-lyase6 (PAL6), CCR2 and F5H2, and decreased activity of PAL. Expression of the other genes was, however, induced by waterlogging, suggesting their role in the synthesis of other phenylpropanoid-derived molecules with roles in stress responses. Moreover, difference in internode lignin content between cultivars or change in its level due to waterlogging is associated with the level of cytokinin.

CONCLUSION

Lodging resistance, tolerance against biotic and abiotic stresses and feedstock quality of wheat biomass are closely associated with its lignin content. Therefore, the findings of this study provide important insights into the molecular mechanisms underlying lignin formation in wheat, an important step towards the development of molecular tools that can facilitate the breeding of wheat cultivars for optimized lignin content and enhanced feedstock quality without affecting other lignin-related agronomic benefits.

Transcriptional Profiles of Hybrid Eucalyptus Genotypes with Contrasting Lignin Content Reveal That Monolignol Biosynthesis-related Genes Regulate Wood Composition

Eucalyptus species constitutes the most widely planted hardwood trees in temperate and subtropical regions. In this study, we compared the transcript levels of genes involved in lignocellulose formation such as cellulose, hemicellulose and lignin biosynthesis in two selected 3-year old hybrid Eucalyptus (Eucalyptus urophylla × Eucalyptus grandis) genotypes (AM063 and AM380) that have different lignin content. AM063 and AM380 had 20.2 and 35.5% of Klason lignin content and 59.0 and 48.2%, α-cellulose contents, respectively. We investigated the correlation between wood properties and transcript levels of wood formation-related genes using RNA-seq with total RNAs extracted from developing xylem tissues at a breast height. Transcript levels of cell wall construction genes such as cellulose synthase (CesA) and sucrose synthase (SUSY) were almost the same in both genotypes. However, AM063 exhibited higher transcript levels of UDP-glucose pyrophosphorylase and xyloglucan endotransglucoxylase than those in AM380. Most monolignol biosynthesis-related isozyme genes showed higher transcript levels in AM380. These results indicate monolignol biosynthesis-related genes may regulate wood composition in Eucalyptus. Flavonoids contents were also observed at much higher levels in AM380 as a result of the elevated transcript levels of common phenylpropanoid pathway genes, phenylalanine ammonium lyase, cinnamate-4-hydroxylase (C4H) and 4-coumarate-CoA ligase (4CL). Secondary plant cell wall formation is regulated by many transcription factors. We analyzed genes encoding NAC, WRKY, AP2/ERF, and KNOX transcription factors and found higher transcript levels of these genes in AM380. We also observed increased transcription of some MYB and LIM domain transcription factors in AM380 compared to AM063. All these results show that genes related to monolignol biosynthesis may regulate the wood composition and help maintain the ratio of cellulose and lignin contents in Eucalyptus plants.

Integrative analysis and expression profiling of secondary cell wall genes in C4 biofuel model Setaria italica reveals targets for lignocellulose bioengineering

Several underutilized grasses have excellent potential for use as bioenergy feedstock due to their lignocellulosic biomass. Genomic tools have enabled identification of lignocellulose biosynthesis genes in several sequenced plants. However, the non-availability of whole genome sequence of bioenergy grasses hinders the study on bioenergy genomics and their genomics-assisted crop improvement. Foxtail millet (Setaria italica L.; Si) is a model crop for studying systems biology of bioenergy grasses. In the present study, a systematic approach has been used for identification of gene families involved in cellulose (CesA/Csl), callose (Gsl) and monolignol biosynthesis (PAL, C4H, 4CL, HCT, C3H, CCoAOMT, F5H, COMT, CCR, CAD) and construction of physical map of foxtail millet. Sequence alignment and phylogenetic analysis of identified proteins showed that monolignol biosynthesis proteins were highly diverse, whereas CesA/Csl and Gsl proteins were homologous to rice and Arabidopsis. Comparative mapping of foxtail millet lignocellulose biosynthesis genes with other C4 panicoid genomes revealed maximum homology with switchgrass, followed by sorghum and maize. Expression profiling of candidate lignocellulose genes in response to different abiotic stresses and hormone treatments showed their differential expression pattern, with significant higher expression of SiGsl12, SiPAL2, SiHCT1, SiF5H2, and SiCAD6 genes. Further, due to the evolutionary conservation of grass genomes, the insights gained from the present study could be extrapolated for identifying genes involved in lignocellulose biosynthesis in other biofuel species for further characterization.

Roles of lignin biosynthesis and regulatory genes in plant development

Lignin is an important factor affecting agricultural traits, biofuel production, and the pulping industry. Most lignin biosynthesis genes and their regulatory genes are expressed mainly in the vascular bundles of stems and leaves, preferentially in tissues undergoing lignification. Other genes are poorly expressed during normal stages of development, but are strongly induced by abiotic or biotic stresses. Some are expressed in non-lignifying tissues such as the shoot apical meristem. Alterations in lignin levels affect plant development. Suppression of lignin biosynthesis genes causes abnormal phenotypes such as collapsed xylem, bending stems, and growth retardation. The loss of expression by genes that function early in the lignin biosynthesis pathway results in more severe developmental phenotypes when compared with plants that have mutations in later genes. Defective lignin deposition is also associated with phenotypes of seed shattering or brittle culm. MYB and NAC transcriptional factors function as switches, and some homeobox proteins negatively control lignin biosynthesis genes. Ectopic deposition caused by overexpression of lignin biosynthesis genes or master switch genes induces curly leaf formation and dwarfism.

Suppression of Arabidopsis peroxidase 72 alters cell wall and phenylpropanoid metabolism

Class III peroxidases are glycoproteins with a major role in cell wall maturation such as lignin formation. Peroxidases are usually present in a high number of isoenzymes, which complicates to assign specific functions to individual peroxidase isoenzymes. Arabidopsis genome encodes for 73 peroxidases, among which AtPrx72 has been shown to participate in lignification. Here, we report by using knock out peroxidase mutants how the disruption of AtPrx72 causes thinner secondary walls in interfascicular fibres but not in the xylem of the stem. This effect is also age-dependent, and AtPrx72 function seems to be particularly important when lignification prevails over elongation processes. Finally, the suppression AtPrx72 leads to the down-regulation of lignin biosynthesis pathway, as well as genes and transcription factors involved in secondary wall thickening.

Functional characterization of a plastidal cation-dependent O-methyltransferase from the liverwort Plagiochasma appendiculatum

Caffeoyl CoA O-methyltransferases (CCoAOMTs), known to be involved in phenylpropanoid metabolism and lignin synthesis, have been characterized from several higher plant species, which also harbor CCoAOMT-like enzymes responsible for methylation of a variety of flavonoids, anthocyanins, coumarins and phenylpropanoids. Here, a gene encoding a CCoAOMT (PaOMT1) was isolated from a sequenced cDNA library of the liverwort species Plagiochasma appendiculatum, a species belonging to the Family Aytoniaceae. The full-length cDNA sequence of PaOMT1 contains 909 bp, and is predicted to encode a protein with 302 amino acids. The gene products were 40-50% identical to CCoAOMT sequences of other plants. Experiments based on recombinant PaOMT1 showed that the enzyme was able to methylate phenylpropanoids, flavonoids and coumarins, with a preference for the flavonoid quercetin (19). Although the substrate selectivity and biochemical feature of PaOMT1 is similar to CCoAOMT-like enzymes, the sequence alignment results indicated PaOMT1 is closer to true CCoAOMT enzymes. A phylogenetic analysis indicated that PaOMT1 is intermediate between true CCoAOMTs and CCoAOMT-like enzymes. The transient expression of a PaOMT1-GFP fusion in tobacco demonstrated that PaOMT1 is directed to the plastids. PaOMT1 may represent an ancestral form of higher plant true CCoAOMT and CCoAOMT-like enzymes. This is the first time an O-methyltransferase was characterized in liverworts.

The Arabidopsis RCC1 Family Protein TCF1 Regulates Freezing Tolerance and Cold Acclimation through Modulating Lignin Biosynthesis

Cell water permeability and cell wall properties are critical to survival of plant cells during freezing, however the underlying molecular mechanisms remain elusive. Here, we report that a specifically cold-induced nuclear protein, Tolerant to Chilling and Freezing 1 (TCF1), interacts with histones H3 and H4 and associates with chromatin containing a target gene, blue-copper-binding protein (BCB), encoding a glycosylphosphatidylinositol-anchored protein that regulates lignin biosynthesis. Loss of TCF1 function leads to reduced BCB transcription through affecting H3K4me2 and H3K27me3 levels within the BCB gene, resulting in reduced lignin content and enhanced freezing tolerance. Furthermore, plants with knocked-down BCB expression (amiRNA-BCB) under cold acclimation had reduced lignin accumulation and increased freezing tolerance. The pal1pal2 double mutant (lignin content reduced by 30% compared with WT) also showed the freezing tolerant phenotype, and TCF1 and BCB act upstream of PALs to regulate lignin content. In addition, TCF1 acts independently of the CBF (C-repeat binding factor) pathway. Our findings delineate a novel molecular pathway linking the TCF1-mediated cold-specific transcriptional program to lignin biosynthesis, thus achieving cell wall remodeling with increased freezing tolerance.

Identification of 4CL Genes in Desert Poplars and Their Changes in Expression in Response to Salt Stress

4-Coumarate:CoA ligase (4CL) genes are critical for the biosynthesis of plant phenylpropanoids. Here we identified 20 4CL genes in the genomes of two desert poplars (Populus euphratica and P. pruinosa) and salt-sensitive congener (P. trichocarpa), but 12 in Salix suchowensis (Salix willow). Phylogenetic analyses clustered all Salicaceae 4CL genes into two clades, and one of them (corresponding to the 4CL-like clade from Arabidopsis) showed signals of adaptive evolution, with more genes retained in Populus than Salix and Arabidopsis. We also found that 4CL12 (in 4CL-like clade) showed positive selection along the two desert poplar lineages. Transcriptional profiling analyses indicated that the expression of 4CL2, 4CL11, and 4CL12 changed significantly in one or both desert poplars in response to salt stress compared to that of in P. trichocarpa. Our results suggest that the evolution of the 4CL genes may have contributed to the development of salt tolerance in the two desert poplars.

Isolation and Functional Characterization of a Phenylalanine Ammonia-Lyase Gene (SsPAL1) from Coleus (Solenostemon scutellarioides (L.) Codd)

Phenylalanine ammonia-lyase (PAL) is the first enzyme involved in the phenylpropanoid pathway and plays important roles in the secondary metabolisms, development and defense of plants. To study the molecular function of PAL in anthocyanin synthesis of Coleus (Solenostemon scutellarioides (L.) Codd), a Coleus PAL gene designated as SsPAL1 was cloned and characterized using a degenerate oligonucleotide primer PCR and RACE method. The full-length SsPAL1 was 2450 bp in size and consisted of one intron and two exons encoding a polypeptide of 711 amino acids. The deduced SsPAL1 protein showed high identities and structural similarities with other functional plant PAL proteins. A series of putative cis-acting elements involved in transcriptional regulation, light and stress responsiveness were found in the upstream regulatory sequence of SsPAL1. Transcription pattern analysis indicated that SsPAL1 was constitutively expressed in all tissues examined and was enhanced by light and different abiotic factors. The recombinant SsPAL1 protein exhibited high PAL activity, at optimal conditions of 60 °C and pH 8.2. Although the levels of total PAL activity and total anthocyanin concentration have a similar variation trend in different Coleus cultivars, there was no significant correlation between them (r = 0.7529, p > 0.1), suggesting that PAL was not the rate-limiting enzyme for the downstream anthocyanin biosynthetic branch in Coleus. This study enables us to further understand the role of SsPAL1 in the phenylpropanoid (flavonoids, anthocyanins) biosynthesis in Coleus at the molecular level.

Exogenous adenosine 5'-phosphoramidate behaves as a signal molecule in plants; it augments metabolism of phenylpropanoids and salicylic acid in Arabidopsis thaliana seedlings

Cells contain various congeners of the canonical nucleotides. Some of these accumulate in cells under stress and may function as signal molecules. Their cellular levels are enzymatically controlled. Previously, we demonstrated a signaling function for diadenosine polyphosphates and cyclic nucleotides in Arabidopsis thaliana and grape, Vitis vinifera. These compounds increased the expression of genes for and the specific activity of enzymes of phenylpropanoid pathways resulting in the accumulation of certain products of these pathways. Here, we show that adenosine 5'-phosphoramidate, whose level can be controlled by HIT-family proteins, induced similar effects. This natural nucleotide, when added to A. thaliana seedlings, activated the genes for phenylalanine:ammonia lyase, 4-coumarate:coenzyme A ligase, cinnamate-4-hydroxylase, chalcone synthase, cinnamoyl-coenzyme A:NADP oxidoreductase and isochorismate synthase, which encode proteins catalyzing key reactions of phenylpropanoid pathways, and caused accumulation of lignins, anthocyanins and salicylic acid. Adenosine 5'-phosphofluoridate, a synthetic congener of adenosine 5'-phosphoramidate, behaved similarly. The results allow us to postulate that adenosine 5'-phosphoramidate should be considered as a novel signaling molecule.

The expression of a rice secondary wall-specific cellulose synthase gene, OsCesA7, is directly regulated by a rice transcription factor, OsMYB58/63

A rice MYB transcription factor, OsMYB58/63, was found to directly upregulate the expression of a rice secondary wall-specific cellulose synthase gene, cellulose synthase A7 ( OsCesA7 ); in contrast, the Arabidopsis putative orthologs AtMYB58 and AtMYB63 have been shown to specifically activate lignin biosynthesis. Although indirect evidence has shown that grass plants are similar to but partially different from dicotyledonous ones in transcriptional regulation of lignocellulose biosynthesis, little is known about the differences. This study showed that a rice MYB transcription factor, OsMYB58/63, directly upregulated the expression of a rice secondary wall-specific cellulose synthase gene, cellulose synthase A7 (OsCesA7). Gene co-expression analysis showed that, in rice, OsMYB58/63 and several rice MYB genes were co-expressed with genes encoding lignocellulose biosynthetic enzymes. The expression levels of OsMYB55/61, OsMYB55/61-L, OsMYB58/63, and OsMYB42/85 were commonly found to be high in culm internodes and nodes. All four MYB transcription factors functioned as transcriptional activators in yeast cells. OsMYB58/63 most strongly transactivated the expression of OsCesA7 in rice protoplasts. Moreover, recombinant OsMYB58/63 protein was bound to two distinct cis-regulatory elements, AC-II and SMRE3, in the OsCesA7 promoter. This is in sharp contrast to the role of Arabidopsis orthologs, AtMYB58 and AtMYB63, which had been reported to specifically activate lignin biosynthesis. The promoter analysis revealed that AC elements, which are the binding sites for MYB58 and MYB63, were lacking in cellulose and xylan biosynthetic genes in Arabidopsis, but present in cellulose, xylan, and lignin biosynthetic genes in rice, implying that the difference of transcriptional regulation between rice and Arabidopsis is due to the distinct composition of promoters. Our results provide a new insight into transcriptional regulation in grass lignocellulose biosynthesis.

Large-Scale Transcriptome Analysis of Two Sugarcane Genotypes Contrasting for Lignin Content

Sugarcane is an important crop worldwide for sugar and first generation ethanol production. Recently, the residue of sugarcane mills, named bagasse, has been considered a promising lignocellulosic biomass to produce the second-generation ethanol. Lignin is a major factor limiting the use of bagasse and other plant lignocellulosic materials to produce second-generation ethanol. Lignin biosynthesis pathway is a complex network and changes in the expression of genes of this pathway have in general led to diverse and undesirable impacts on plant structure and physiology. Despite its economic importance, sugarcane genome was still not sequenced. In this study a high-throughput transcriptome evaluation of two sugarcane genotypes contrasting for lignin content was carried out. We generated a set of 85,151 transcripts of sugarcane using RNA-seq and de novo assembling. More than 2,000 transcripts showed differential expression between the genotypes, including several genes involved in the lignin biosynthetic pathway. This information can give valuable knowledge on the lignin biosynthesis and its interactions with other metabolic pathways in the complex sugarcane genome.

Poplar PdMYB221 is involved in the direct and indirect regulation of secondary wall biosynthesis during wood formation

Wood is formed by the successive addition of secondary xylem, which consists of cells with a conspicuously thickened secondary wall composed mainly of cellulose, xylan and lignin. Currently, few transcription factors involved in the direct regulation of secondary wall biosynthesis have been characterized in tree species. Here, we show that PdMYB221, a poplar ortholog of the Arabidopsis R2R3-MYB transcription factor AtMYB4, directly regulates secondary wall biosynthesis during wood formation. PdMYB221 is predominantly expressed in cells of developing wood, and the protein it encodes localizes to the nucleus and acts as a transcriptional repressor. Ectopic expression of PdMYB221 resulted in reduced cell wall thicknesses of fibers and vessels in Arabidopsis inflorescence stems. The amounts of cellulose, xylose, and lignin were decreased and the expression of key genes synthesizing the three components was suppressed in PdMYB221 overexpression plants. Transcriptional activation assays showed that PdMYB221 repressed the promoters of poplar PdCESA7/8, PdGT47C, PdCOMT2 and PdCCR1. Electrophoretic mobility shift assays revealed that PdMYB221 bound directly to the PdCESA8, PdGT47C, and PdCOMT2 promoters. Together, our results suggest that PdMYB221 may be involved in the negative regulation of secondary wall formation through the direct and indirect suppression of the gene expression of secondary wall biosynthesis.

Biochemical characterization of caffeoyl coenzyme A 3-O-methyltransferase from wheat

TaCCoAOMT1 is located in wheat chromosome 7A and highly expressed in stem and root. It is important for lignin biosynthesis, and associated with stem maturity but not lodging resistance. Caffeoyl coenzyme A 3-O-methyltransferases (CCoAOMTs) are one important class of enzymes to carry out the transfer of the methyl group from S-adenosylmethionine to the hydroxyl group, and play important roles in lignin and flavonoids biosynthesis. In the present study, sequences for CCoAOMT from the wheat genome were analyzed. One wheat CCoAOMT that belonged to bona fide subclade involved in lignin biosynthesis, namely TaCCoAOMT1, was obtained by the prokaryotic expression in E. coli. The three-dimensional structure prediction showed a highly similar structure of TaCCoAOMT1 with MsCCoAOMT. Recombinant TaCCoAOMT1 protein could only use caffeoyl CoA and 5-hydroxyferuloyl CoA as effective substrates and caffeoyl CoA as the best substrate. TaCCoAOMT1 had a narrow optimal pH and thermal stability. The TaCCoAOMT1 gene was highly expressed in wheat stem and root tissues, paralleled CCoAOMT enzyme activity. TaCCoAOMT1 mRNA abundance and enzyme activity increased linearly with stem maturity, but showed little difference between wheat lodging-resistant (H4546) and lodging-sensitive (C6001) cultivars in elongation, heading and milky stages. These data suggest that TaCCoAOMT1 is an important CCoAOMT for lignin biosynthesis that is critical for stem development, but not directly associated with lodging-resistant trait in wheat.

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome-wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.

Cell wall modifications triggered by the down-regulation of Coumarate 3-hydroxylase-1 in maize

Coumarate 3-hydroxylase (C3H) catalyzes a key step of the synthesis of the two main lignin subunits, guaiacyl (G) and syringyl (S) in dicotyledonous species. As no functional data are available in regards to this enzyme in monocotyledonous species, we generated C3H1 knock-down maize plants. The results obtained indicate that C3H1 participates in lignin biosynthesis as its down-regulation redirects the phenylpropanoid flux: as a result, increased amounts of p-hydroxyphenyl (H) units, lignin-associated ferulates and the flavone tricin were detected in transgenic stems cell walls. Altogether, these changes make stem cell walls more degradable in the most C3H1-repressed plants, despite their unaltered polysaccharide content. The increase in H monomers is moderate compared to C3H deficient Arabidopsis and alfalfa plants. This could be due to the existence of a second maize C3H protein (C3H2) that can compensate the reduced levels of C3H1 in these C3H1-RNAi maize plants. The reduced expression of C3H1 alters the macroscopic phenotype of the plants, whose growth is inhibited proportionally to the extent of C3H1 repression. Finally, the down-regulation of C3H1 also increases the synthesis of flavonoids, leading to the accumulation of anthocyanins in transgenic leaves.

Genome-wide analysis of the lignin toolbox of Eucalyptus grandis

Lignin, a major component of secondary cell walls, hinders the optimal processing of wood for industrial uses. The recent availability of the Eucalyptus grandis genome sequence allows comprehensive analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway and identification of those mainly involved in xylem developmental lignification. We performed genome-wide identification of putative members of the lignin gene families, followed by comparative phylogenetic studies focusing on bona fide clades inferred from genes functionally characterized in other species. RNA-seq and microfluid real-time quantitative PCR (RT-qPCR) expression data were used to investigate the developmental and environmental responsive expression patterns of the genes. The phylogenetic analysis revealed that 38 E. grandis genes are located in bona fide lignification clades. Four multigene families (shikimate O-hydroxycinnamoyltransferase (HCT), p-coumarate 3-hydroxylase (C3H), caffeate/5-hydroxyferulate O-methyltransferase (COMT) and phenylalanine ammonia-lyase (PAL)) are expanded by tandem gene duplication compared with other plant species. Seventeen of the 38 genes exhibited strong, preferential expression in highly lignified tissues, probably representing the E. grandis core lignification toolbox. The identification of major genes involved in lignin biosynthesis in E. grandis, the most widely planted hardwood crop world-wide, provides the foundation for the development of biotechnology approaches to develop tree varieties with enhanced processing qualities.

Functional characterization of CCR in birch (Betula platyphylla × Betula pendula) through overexpression and suppression analysis

We cloned a Cinnamoyl-CoA Reductase gene (BpCCR1) from an apical meristem and first internode of Betula platyphylla and characterized its functions in lignin biosynthesis, wood formation and tree growth through transgenic approaches. We generated overexpression and suppression transgenic lines and analyzed them in comparison with the wild-type in terms of lignin content, anatomical characteristics, height and biomass. We found that BpCCR1 overexpression could increase lignin content up to 14.6%, and its underexpression decreased lignin content by 6.3%. Surprisingly, modification of BpCCR1 expression led to conspicuous changes in wood characteristics, including xylem vessel number and arrangement, and secondary wall thickness. The growth of transgenic trees in terms of height was also significantly influenced by the modification of BpCCR1 genes. We discuss the functions of BpCCR1 in the context of a phylogenetic tree built with CCR genes from multiple species.

Ancient origin of the biosynthesis of lignin precursors

BACKGROUND

Lignin plays an important role in plant structural support and water transport, and is considered one of the hallmarks of land plants. The recent discovery of lignin or its precursors in various algae has raised questions on the evolution of its biosynthetic pathway, which could be much more ancient than previously thought. To determine the taxonomic distribution of the lignin biosynthesis genes, we screened all publicly available genomes of algae and their closest non-photosynthetic relatives, as well as representative land plants. We also performed phylogenetic analysis of these genes to decipher the evolution and origin(s) of lignin biosynthesis.

RESULTS

Enzymes involved in making p-coumaryl alcohol, the simplest lignin monomer, are found in a variety of photosynthetic eukaryotes, including diatoms, dinoflagellates, haptophytes, cryptophytes as well as green and red algae. Phylogenetic analysis of these enzymes suggests that they are ancient and spread to some secondarily photosynthetic lineages when they acquired red and/or green algal endosymbionts. In some cases, one or more of these enzymes was likely acquired through lateral gene transfer (LGT) from bacteria.

CONCLUSIONS

Genes associated with p-coumaryl alcohol biosynthesis are likely to have evolved long before the transition of photosynthetic eukaryotes to land. The original function of this lignin precursor is therefore unlikely to have been related to water transport. We suggest that it participates in the biological defense of some unicellular and multicellular algae.

RNA-Seq Based De Novo Transcriptome Assembly and Gene Discovery of Cistanche deserticola Fleshy Stem

Backgrounds

Cistanche deserticola is a completely non-photosynthetic parasitic plant with great medicinal value and mainly distributed in desert of Northwest China. Its dried fleshy stem is a crucial tonic in traditional Chinese medicine with roles of mainly improving male sexual function and strengthening immunity, but few mechanistic studies have been conducted partly due to the lack of genomic and transcriptomic resources.

Results

In this study, we performed deep transcriptome sequencing in fleshy stem of C. deserticola, and about 80 million reads were generated using Illumina pair-end sequencing on HiSeq2000 platform. Using trinity assembler, we obtained 95,787 transcript sequences with transcript lengths ranging from 200bp to 15,698bp, having an average length of 950 bases and the N50 length of 1,519 bases. 63,957 transcripts were identified actively expressed with FPKM ≥ 0.5, in which 30,098 transcripts were annotated with gene descriptions or gene ontology terms by sequence similarity analyses against several public databases (Uniprot, NR and Nt at NCBI, and KEGG). Furthermore, we identified key enzyme genes involved in biosynthesis of lignin and phenylethanoid glycosides (PhGs) which are known to be the primary active ingredients. Four phenylalanine ammonia-lyase (PAL) genes, the first key enzyme in lignin and PhG biosynthesis, were identified based on sequences comparison and phylogenetic analysis. Two biosynthesis pathways of PhGs were also proposed for the first time.

Conclusions

In all, we completed a global analysis of the C. deserticola fleshy stem transcriptome using RNA-seq technology. A collection of enzyme genes related to biosynthesis of lignin and phenylethanoid glysides were identified from the assembled and annotated transcripts, and the gene family of PAL was also predicted. The sequence data from this study will provide a valuable resource for conducting future phenylethanoid glysides biosynthesis researches and functional genomic studies in this important medicinal plant.

High-level hemicellulosic arabinose predominately affects lignocellulose crystallinity for genetically enhancing both plant lodging resistance and biomass enzymatic digestibility in rice mutants

Rice is a major food crop with enormous biomass residue for biofuels. As plant cell wall recalcitrance basically decides a costly biomass process, genetic modification of plant cell walls has been regarded as a promising solution. However, due to structural complexity and functional diversity of plant cell walls, it becomes essential to identify the key factors of cell wall modifications that could not much alter plant growth, but cause an enhancement in biomass enzymatic digestibility. To address this issue, we performed systems biology analyses of a total of 36 distinct cell wall mutants of rice. As a result, cellulose crystallinity (CrI) was examined to be the key factor that negatively determines either the biomass enzymatic saccharification upon various chemical pretreatments or the plant lodging resistance, an integrated agronomic trait in plant growth and grain production. Notably, hemicellulosic arabinose (Ara) was detected to be the major factor that negatively affects cellulose CrI probably through its interlinking with β-1,4-glucans. In addition, lignin and G monomer also exhibited the positive impact on biomass digestion and lodging resistance. Further characterization of two elite mutants, Osfc17 and Osfc30, showing normal plant growth and high biomass enzymatic digestion in situ and in vitro, revealed the multiple GH9B candidate genes for reducing cellulose CrI and XAT genes for increasing hemicellulosic Ara level. Hence, the results have suggested the potential cell wall modifications for enhancing both biomass enzymatic digestibility and plant lodging resistance by synchronically overexpressing GH9B and XAT genes in rice.

Indole Glucosinolate Biosynthesis Limits Phenylpropanoid Accumulation in Arabidopsis thaliana

Plants produce an array of metabolites (including lignin monomers and soluble UV-protective metabolites) from phenylalanine through the phenylpropanoid biosynthetic pathway. A subset of plants, including many related to Arabidopsis thaliana, synthesizes glucosinolates, nitrogen- and sulfur-containing secondary metabolites that serve as components of a plant defense system that deters herbivores and pathogens. Here, we report that the Arabidopsis thaliana reduced epidermal fluorescence5 (ref5-1) mutant, identified in a screen for plants with defects in soluble phenylpropanoid accumulation, has a missense mutation in CYP83B1 and displays defects in glucosinolate biosynthesis and in phenylpropanoid accumulation. CYP79B2 and CYP79B3 are responsible for the production of the CYP83B1 substrate indole-3-acetaldoxime (IAOx), and we found that the phenylpropanoid content of cyp79b2 cyp79b3 and ref5-1 cyp79b2 cyp79b3 plants is increased compared with the wild type. These data suggest that levels of IAOx or a subsequent metabolite negatively influence phenylpropanoid accumulation in ref5 and more importantly that this crosstalk is relevant in the wild type. Additional biochemical and genetic evidence indicates that this inhibition impacts the early steps of the phenylpropanoid biosynthetic pathway and restoration of phenylpropanoid accumulation in a ref5-1 med5a/b triple mutant suggests that the function of the Mediator complex is required for the crosstalk.

Cloning and in silico analysis of a cinnamyl alcohol dehydrogenase gene in Pennisetum purpureum

Lignin is a major constituent of plant cell walls and indispensable to the normal growth of a plant. However, the presence of lignin complicates the structure of the plant cell walls and negatively influences pulping industry, lignocellulose utilization as well as forage properties. Cinnamyl alcohol dehydrogenase (CAD), a key enzyme involved in lignin biosynthesis, catalyses the last step in monolignol synthesis and has a major role in genetic regulation of lignin production. In the present study, a 1 342-bp cDNA fragment of CAD gene, named PpCAD, was isolated from Pennisetum purpureum using strategies of homologous clone and rapid amplification of cDNA end. It was translated into an intact protein sequence including 366 amino acid residues by ORF Finder. The genomic full-length DNA of PpCAD was a 3 738-bp sequence containing four exons and three introns, among which the 114-bp exon was considered to be a conserved region compared with other CADs. Basic bioinformatic analysis presumed that the PpCAD was a nonsecretory and hydrophobic protein with five possible transmembrane helices. The phylogenetic analysis indicated that the PpCAD belonged to the class of bona fide CADs involved in lignin synthesis and it showed a high similarity (nearly 90%) with CAD protein sequences of Sorghum bicolor, Panicum virgatum and Zea mays in Gramineae. Furthere, PpCAD amino acid sequence was demonstrated to have some conserved motifs such as Zn-binding site, Zn-catalytic centre and NADP(H) binding domain after aligning with other bona fide CADs. Three-dimensional homology modelling of PpCAD showed that the protein had some exclusive features of bona fide CADs.

Sclareol induces plant resistance to root-knot nematode partially through ethylene-dependent enhancement of lignin accumulation

The root-knot nematode (RKN) is one of the most devastating parasitic nematodes of plants. Although some secondary metabolites released by the host plant play roles as defense substances against parasitic nematodes, the mechanism underlying the induction of such defense responses is not fully understood. We found that sclareol, a natural diterpene known as an antimicrobial and defense-related molecule, inhibited RKN penetration of tomato and Arabidopsis roots. Sclareol induced genes related to ethylene (ET) biosynthesis and signaling and phenylpropanoid metabolism in Arabidopsis roots. In roots of ein2-1, an ET-insensitive mutant line, both sclareol-induced inhibition of RKN penetration and sclareol-induced enhancement of lignin accumulation were abolished. A mutant defective in lignin accumulation did not exhibit such inhibition. Sclareol also activated MPK3 and MPK6, Arabidopsis mitogen-activated protein kinases whose activation is required for triggering ET biosynthesis. Sclareol-induced inhibition of RKN penetration was exhibited by mutants of neither MPK3 nor MPK6. Treatment with a biosynthetic precursor of ET was insufficient compared with sclareol treatment to inhibit RKN penetration, suggesting the existence of an ET-independent signaling pathway leading to RKN resistance. These results suggested that sclareol induced resistance to RKN penetration partially through ET-dependent accumulation of lignin in roots.

Genome-wide identification of phenolic acid biosynthetic genes in Salvia miltiorrhiza

MAIN CONCLUSION

Twenty-nine genes related to phenolic acid biosynthesis were identified in the Salvia miltiorrhiza genome. Nineteen of these are described for the first time, with ten genes experimentally correlating to phenolic acid biosynthesis. Vast stores of secondary metabolites exist in plants, many of which possess biological activities related to human health. Phenolic acid derivatives are a class of valuable bioactive pharmaceuticals abundant in the widely used Chinese medicinal herb, Salvia miltiorrhiza. The biosynthetic pathway for phenolic acids differs in this species from that of other investigated plants. However, the molecular basis for this is unknown, with systematic analysis of the genes involved not yet performed. As the first step towards unraveling this complex biosynthetic pathway in S. miltiorrhiza, the current genome assembly was searched for putatively involved genes. Twenty-nine genes were revealed, 19 of which are described here for the first time. These include 15 genes predicted in the phenylpropanoid pathway; seven genes in the tyrosine-derived pathway; six genes encoding putative hydroxycinnamoyltransferases, and one CYP98A, namely CYP98A78. The promoter regions, gene structures and expression patterns of these genes were examined. Furthermore, conserved domains and phylogenetic relationships with homologous proteins in other species were revealed. Most of the key enzymes, including 4-coumarate: CoA ligase, 4-hydroxyphenylpyruvate reductase and hydroxycinnamoyltransferase, were found in multiple copies, each exhibiting different characteristics. Ten genes putatively involved in rosmarinic acid biosynthesis are also described. These findings provide a foundation for further analysis of this complex and diverse pathway, with potential to enhance the synthesis of water-soluble medicinal compounds in S. miltiorrhiza.

Pathways associated with lignin biosynthesis in lignomaniac jute fibres

We generated the bast transcriptomes of a deficient lignified phloem fibre mutant and its wild-type jute (Corchorus capsularis) using Illumina paired-end sequencing. A total of 34,163 wild-type and 29,463 mutant unigenes, with average lengths of 1442 and 1136 bp, respectively, were assembled de novo, ~77-79 % of which were functionally annotated. These annotated unigenes were assigned to COG (~37-40 %) and GO (~22-28 %) classifications and mapped to 189 KEGG pathways (~19-21 %). We discovered 38 and 43 isoforms of 16 and 10 genes of the upstream shikimate-aromatic amino acid and downstream monolignol biosynthetic pathways, respectively, rendered their sequence similarities, confirmed the identities of 22 of these candidate gene families by phylogenetic analyses and reconstructed the pathway leading to lignin biosynthesis in jute fibres. We also identified major genes and bast-related transcription factors involved in secondary cell wall (SCW) formation. The quantitative RT-PCRs revealed that phenylalanine ammonia-lyase 1 (CcPAL1) was co-down-regulated with several genes of the upstream shikimate pathway in mutant bast tissues at an early growth stage, although its expression relapsed to the normal level at the later growth stage. However, cinnamyl alcohol dehydrogenase 7 (CcCAD7) was strongly down-regulated in mutant bast tissues irrespective of growth stages. CcCAD7 disruption at an early growth stage was accompanied by co-up-regulation of SCW-specific genes cellulose synthase A7 (CcCesA7) and fasciclin-like arabinogalactan 6 (CcFLA6), which was predicted to be involved in coordinating the S-layers' deposition in the xylan-type jute fibres. Our results identified CAD as a promising target for developing low-lignin jute fibres using genomics-assisted molecular approaches.