Cinnamoyl CoA reductase, the first committed enzyme of the lignin branch biosynthetic pathway: cloning, expression and phylogenetic relationships

Plant terrestrialization: an environmental pull on the evolution of multi-sourced streptophyte phenolics

Phenolic compounds of land plants are varied: they are chemodiverse, are sourced from different biosynthetic routes and fulfil a broad spectrum of functions that range from signalling phytohormones, to protective shields against stressors, to structural compounds. Their action defines the biology of land plants as we know it. Often, their roles are tied to environmental responses that, however, impacted already the algal progenitors of land plants, streptophyte algae. Indeed, many streptophyte algae successfully dwell in terrestrial habitats and have homologues for enzymatic routes for the production of important phenolic compounds, such as the phenylpropanoid pathway. Here, we synthesize what is known about the production of specialized phenolic compounds across hundreds of millions of years of streptophyte evolution. We propose an evolutionary scenario in which selective pressures borne out of environmental cues shaped the chemodiversity of phenolics in streptophytes. This article is part of the theme issue 'The evolution of plant metabolism'.

Multi-omics analyses reveal the mechanisms underlying the responses of Casuarina equisetifolia ssp. incana to seawater atomization and encroachment stress

Casuarina equisetifolia trees are used as windbreaks in subtropical and tropical coastal zones, while C. equisetifolia windbreak forests can be degraded by seawater atomization (SA) and seawater encroachment (SE). To investigate the mechanisms underlying the response of C. equisetifolia to SA and SE stress, the transcriptome and metabolome of C. equisetifolia seedlings treated with control, SA, and SE treatments were analyzed. We identified 737, 3232, 3138, and 3899 differentially expressed genes (SA and SE for 2 and 24 h), and 46, 66, 62, and 65 differentially accumulated metabolites (SA and SE for 12 and 24 h). The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that SA and SE stress significantly altered the expression of genes related to plant hormone signal transduction, plant-pathogen interaction, and starch and sucrose metabolism pathways. The accumulation of metabolites associated with the biosynthetic pathways of phenylpropanoid and amino acids, as well as starch and sucrose metabolism, and glycolysis/gluconeogenesis were significantly altered in C. equisetifolia subjected to SA and SE stress. In conclusion, C. equisetifolia responds to SA and SE stress by regulating plant hormone signal transduction, plant-pathogen interaction, biosynthesis of phenylpropanoid and amino acids, starch and sucrose metabolism, and glycolysis/gluconeogenesis pathways. Compared with SA stress, C. equisetifolia had a stronger perception and response to SE stress, which required more genes and metabolites to be regulated. This study enhances our understandings of how C. equisetifolia responds to two types of seawater stresses at transcriptional and metabolic levels. It also offers a theoretical framework for effective coastal vegetation management in tropical and subtropical regions.

Functional impacts of PtrMYB203 on phenylpropanoid pathway regulation and wood properties in hybrid poplar

The phenylpropanoid pathway is vital for plant growth and development, producing lignin and flavonoids. This study investigates PtrMYB203, a homolog of MYB repressors of proanthocyanidin (PA) biosynthesis in Populus trichocarpa, as a transcriptional repressor in the phenylpropanoid pathway of hybrid poplar (Populus alba x P. glandulosa). Overexpression of PtrMYB203 (35S::PtrMYB203) in hybrid poplar detrimentally impacted plant growth and development. Histological analysis revealed irregular xylem vessel formation and decreased lignin content, corroborated by Klason lignin assays. Moreover, 35S::PtrMYB203 transgenic poplars exhibited significant decreases in anthocyanin and PA accumulations in callus tissues, even under high light conditions. Quantitative RT-PCR analysis and protoplast-based transcriptional activation assay confirmed the downregulation of lignin and flavonoid biosynthesis genes. This genetic modification also alters the expression of several MYB transcription factors, essential for phenylpropanoid pathway regulation. Remarkably, saccharification efficiency in the 35S::PtrMYB203 poplar was improved by over 34% following hot water treatment alone. These findings suggest PtrMYB203 as a potential genetic target for enhancing wood properties for bioenergy production, providing valuable insights into the manipulation of metabolite pathways in woody perennials to advance wood biotechnology.

Transcriptome analysis reveals activation of detoxification and defense mechanisms in smoke-exposed Merlot grape (Vitis vinifera) berries

A significant consequence of climate change is the rising incidence of wildfires. When wildfires occur close to wine grape (Vitis vinifera) production areas, smoke-derived volatile phenolic compounds can be taken up by the grape berries, negatively affecting the flavor and aroma profile of the resulting wine and compromising the production value of entire vineyards. Evidence for the permeation of smoke-associated compounds into grape berries has been provided through metabolomics; however, the basis for grapevines' response to smoke at the gene expression level has not been investigated in detail. To address this knowledge gap, we employed time-course RNA sequencing to observe gene expression-level changes in grape berries in response to smoke exposure. Significant increases in gene expression (and enrichment of gene ontologies) associated with detoxification of reactive compounds, maintenance of redox homeostasis, and cell wall fortification were observed in response to smoke. These findings suggest that the accumulation of volatile phenols from smoke exposure activates mechanisms that render smoke-derived compounds less reactive while simultaneously fortifying intracellular defense mechanisms. The results of this work lend a better understanding of the molecular basis for grapevines' response to smoke and provide insight into the origins of smoke-taint-associated flavor and aroma attributes in wine produced from smoke-exposed grapes.

Identification of key genes involved in lignin and flavonoid accumulation during Tilia tuan seed maturation

KEY MESSAGE

Transcriptomics and phenotypic data analysis identified 24 transcription factors (TFs) that play key roles in regulating the competitive accumulation of lignin and flavonoids. Tilia tuan Szyszyl. (T. tuan) is a timber tree species with important ecological and commercial value. However, its highly lignified pericarp results in a low seed germination rate and a long dormancy period. In addition, it is unknown whether there is an interaction between the biosynthesis of flavonoids and lignin as products of the phenylpropanoid pathway during seed development. To explore the molecular regulatory mechanism of lignin and flavonoid biosynthesis, T. tuan seeds were harvested at five stages (30, 60, 90, 120, and 150 days after pollination) for lignin and flavonoid analyses. The results showed that lignin accumulated rapidly in the early and middle stages (S1, S3, and S4), and rapid accumulation of flavonoids during the early and late stages (S1 and S5). High-throughput RNA sequencing analysis of developing seeds identified 50,553 transcripts, including 223 phenylpropanoid biosynthetic pathway genes involved in lignin accumulation grouped into 3 clusters, and 106 flavonoid biosynthetic pathway genes (FBPGs) grouped into 2 clusters. Subsequent WGCNA and time-ordered gene co-expression network (TO-GCN) analysis revealed that 24 TFs (e.g., TtARF2 and TtWRKY15) were involved in flavonoids and lignin biosynthesis regulation. The transcriptome data were validated by qRT-PCR to analyze the expression profiles of key enzyme-coding genes. This study revealed that there existed a competitive relationship between flavonoid and lignin biosynthesis pathway during the development of T. tuan seeds, that provide a foundation for the further exploration of molecular mechanisms underlying lignin and flavonoid accumulation in T. tuan seeds.

GhTKPR1_8 functions to inhibit anther dehiscence and reduce pollen viability in cotton

Sporopollenin, as the main component of the pollen exine, is a highly resistant polymer that provides structural integrity under unfavourable environmental conditions. Tetraketone α-pyrone reductase 1 (TKPR1) is essential for sporopollenin formation, catalyzing the reduction of tetraketone carbonyl to hydroxylated α-pyrone. The functional role of TKPR1 in male sterility has been reported in flowering plants such as maize, rice, and Arabidopsis. However, the molecular cloning and functional characterization of TKPR1 in cotton remain unaddressed. In this study, we identified 68 TKPR1s from four cotton species, categorized into three clades. Transcriptomics and RT-qPCR demonstrated that GhTKPR1_8 exhibited typical expression patterns in the tetrad stage of the anther. GhTKPR1_8 was localized to the endoplasmic reticulum. Moreover, ABORTED MICROSPORES (GhAMS) transcriptionally activated GhTKPR1_8 as indicated by luciferase complementation tests. GhTKPR1_8-knockdown inhibited anther dehiscence and reduced pollen viability in cotton. Additionally, overexpression of GhTKPR1_8 in the attkpr1 mutant restored its male sterile phenotype. This study offers novel insights into the investigation of TKPR1 in cotton while providing genetic resources for studying male sterility.

The key pathways for drought tolerance in Cerasus humilis were unveiled through transcriptome analysis

Drought stress exerts a significant impact on the growth, development, and yield of fruit trees. Cerasus humilis is an endemic drought-resistant fruit tree in northern China. To elucidate the underlying mechanism of drought resistance in C. humilis, comprehensive physiological measurements and transcriptome analysis were conducted on the leaves of C. humilis subjected to 15- or 22-days of drought stress. We identified multiple GO terms and KEGG pathways associated with the drought stress response by performing GO and KEGG analysis on DEGs. Furthermore, through the prediction of transcription factors (TFs) and analysis of their expression levels, we observed differential expression patterns among most members of stress-responsive TF families as the duration of drought stress increased. WGCNA analysis was performed on the transcriptome to identify gene cluster modules that exhibited a strong correlation with the durations of drought. Subsequently, these modules underwent GO and KEGG enrichment analyses. The study revealed that the TF-mediated lignin biosynthesis pathway, along with the plant hormone signal transduction pathway, played a prominent role in responding to drought stress of C. humilis. Gene profiling analysis, qRT-PCR, and determination of phytohormone and lignin contents further supported this hypothesis. The hierarchical gene regulatory network was finally constructed based on DEGs from the aforementioned key enriched pathways to predict the gene regulatory mechanisms in response to stress for C. humilis. The findings from this study provide valuable insights into how C. humilis copes with drought stress while analyzing crucial gene pathways associated with its resistance from a TF perspective. This research is significant for the genetic breeding of economic forests.

Identification of a cinnamoyl-CoA reductase from Cinnamomum cassia involved in trans-cinnamaldehyde biosynthesis

MAIN CONCLUSION

The identification of a functional cinnamoyl-CoA reductase enzyme from Cinnamomum cassia involved in trans-cinnamaldehyde biosynthesis offers the potential for enhancing trans-cinnamaldehyde production through genetic engineering. A significant accumulation of trans-cinnamaldehyde has been found in the bark tissues of C. cassia, used in traditional Chinese medicine. trans-Cinnamaldehyde exhibits various pharmacological properties such as anti-inflammatory, analgesic, and protection of the stomach and the digestive tract. However, further elucidation and characterization of the biosynthetic pathway for trans-cinnamaldehyde is required. In this study, we conducted an integrated analysis of trans-cinnamaldehyde accumulation profiles and transcriptomic data from five different C. cassia tissues to identify the genes involved in its biosynthesis. The transcriptome data we obtained included nearly all genes associated with the trans-cinnamaldehyde pathway, with the majority demonstrating high abundance in branch barks and trunk barks. We successfully cloned four C. cassia cinnamoyl-CoA reductases (CcCCRs), a key gene in trans-cinnamaldehyde biosynthesis. We found that the recombinant CcCCR1 protein was the only one that more efficiently converted cinnamoyl-CoA into trans-cinnamaldehyde. CcCCR1 exhibited approximately 14.7-fold higher catalytic efficiency (kcat/Km) compared to the Arabidopsis thaliana cinnamoyl-CoA reductase 1 (AtCCR1); therefore, it can be utilized for engineering higher trans-cinnamaldehyde production as previously reported. Molecular docking studies and mutagenesis experiments also validated the superior catalytic activity of CcCCR1 compared to AtCCR1. These findings provide valuable insights for the functional characterization of enzyme-coding genes and hold potential for future engineering of trans-cinnamaldehyde biosynthetic pathways.

Novel lignin-based extracellular barrier in glandular trichome

Successful biochemical reactions in organisms necessitate compartmentalization of the requisite components. Glandular trichomes (GTs) act as compartments for the synthesis and storage of specialized compounds. These compounds not only are crucial for the survival of plants under biotic and abiotic stresses but also have medical and commercial value for humans. However, the mechanisms underlying compartmentalization remain unclear. Here we identified a novel structure that is indispensable for the establishment of compartments in cucumber GTs. Silica, a specialized compound, is deposited on the GTs and is visible on the surface of the fruit as a white powder, known as bloom. This deposition provides resistance against pathogens and prevents water loss from the fruits1. Using the cucumber bloomless mutant2, we discovered that a lignin-based cell wall structure in GTs, named 'neck strip', achieves compartmentalization by acting as an extracellular barrier crucial for the silica polymerization. This structure is present in the GTs of diverse plant species. Our findings will enhance the understanding of the biosynthesis of unique compounds in trichomes and provide a basis for improving the production of compounds beneficial to humans.

Integrated Transcriptome and Metabolome Analyses Reveal Bamboo Culm Color Formation Mechanisms Involved in Anthocyanin Biosynthetic in Phyllostachys nigra

Phyllostachys nigra has green young culms (S1) and purple black mature culms (S4). Anthocyanins are the principal pigment responsible for color presentation in ornamental plants. We employ a multi-omics approach to investigate the regulatory mechanisms of anthocyanins in Ph. nigra. Firstly, we found that the pigments of the culm of Ph. nigra accumulated only in one to four layers of cells below the epidermis. The levels of total anthocyanins and total flavonoids gradually increased during the process of bamboo culm color formation. Metabolomics analysis indicated that the predominant pigment metabolites observed were petunidin 3-O-glucoside and malvidin O-hexoside, exhibiting a significant increase of up to 9.36-fold and 13.23-fold, respectively, during pigmentation of Ph. nigra culm. Transcriptomics sequencing has revealed that genes involved in flavonoid biosynthesis, phenylpropanoid biosynthesis, and starch and sucrose metabolism pathways were significantly enriched, leading to color formation. A total of 62 differentially expressed structural genes associated with anthocyanin synthesis were identified. Notably, PnANS2, PnUFGT2, PnCHI2, and PnCHS1 showed significant correlations with anthocyanin metabolites. Additionally, certain transcription factors such as PnMYB6 and PnMYB1 showed significant positive or negative correlations with anthocyanins. With the accumulation of sucrose, the expression of PnMYB6 is enhanced, which in turn triggers the expression of anthocyanin biosynthesis genes. Based on these findings, we propose that these key genes primarily regulate the anthocyanin synthesis pathway in the culm and contribute to the accumulation of anthocyanin, ultimately resulting in the purple-black coloration of Ph. nigra.

Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.)

BACKGROUND

Teak (Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak (Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2).

METHODS AND RESULTS

The CAD1, MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m-3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0-21% SNP pairs) at P < 0.05 and high LD decay (33-934 bp) at R2 = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly (P < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06-7.04 × 10-4) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait.

CONCLUSION

Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.

Reinterpretation of anthocyanins biosynthesis in developing black rice seeds through gene expression analysis

The biosynthesis of anthocyanins is still questionable in regulating the quantities of anthocyanins biosynthesized in rice seeds and the expression levels of transcription factors and the structural genes involved in the biosynthetic pathway of anthocyanins. We herein investigated the relationship between the accumulated anthocyanin contents and the expression levels of genes related to the biosynthesis of anthocyanins in rice seeds. Liquid chromatography/mass spectrometry-mass spectrometry analysis of cyanidin 3-glucoside (C3G) in rice seeds showed no accumulation of C3G in white and red rice cultivars, and the differential accumulation of C3G among black rice cultivars. RNA-seq analysis in rice seeds, including white, red, and black rice cultivars, at twenty days after heading (DAH) further exhibited that the genes involved in the biosynthesis of anthocyanins were differentially upregulated in developing seeds of black rice. We further verified these RNA-seq results through gene expression analysis by a quantitative real-time polymerase chain reaction in developing seeds of white, red, and black rice cultivars at 20 DAH. Of these genes related to the biosynthesis of anthocyanins, bHLHs, MYBs, and WD40, which are regulators, and the structural genes, including chalcone synthase (CHS), flavanone 3-hydroxylase (F3H), flavonoid 3´-hydroxylase (F3´H), dihydroflavonol 4-reductase (DFR), and anthocyanidin synthase (ANS), were differentially upregulated in black rice seeds. The correlation analysis revealed that the quantities of C3G biosynthesized in black rice seeds were positively correlated to the expression levels of bHLHs, MYBs and WD40, CHS, F3H, F3´H, DFR, and ANS. In addition, we present bHLH2 (LOC_Os04g47040) and MYBs (LOC_Os01g49160, LOC_Os01g74410, and LOC_Os03g29614) as new putative transcription factor genes for the biosynthesis of anthocyanins in black rice seeds. It is expected that this study will help to improve the understanding of the molecular levels involved in the biosynthesis of anthocyanins in black rice seeds.

Transcriptomic Analysis of Heat Stress Response in Brassica rapa L. ssp. pekinensis with Improved Thermotolerance through Exogenous Glycine Betaine

Chinese cabbage (Brassica rapa L. ssp. pekinensis) is sensitive to high temperature, which will cause the B. rapa to remain in a semi-dormancy state. Foliar spray of GB prior to heat stress was proven to enhance B. rapa thermotolerance. In order to understand the molecular mechanisms of GB-primed resistance or adaptation towards heat stress, we investigated the transcriptomes of GB-primed and non-primed heat-sensitive B. rapa ‘Beijing No. 3’ variety by RNA-Seq analysis. A total of 582 differentially expressed genes (DEGs) were identified from GB-primed plants exposed to heat stress relative to non-primed plants under heat stress and were assigned to 350 gene ontology (GO) pathways and 69 KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. The analysis of the KEGG enrichment pathways revealed that the most abundantly up-regulated pathways were protein processing in endoplasmic reticulum (14 genes), followed by plant hormone signal transduction (12 genes), ribosome (8 genes), MAPK signaling pathway (8 genes), homologous recombination (7 genes), nucleotide excision repair metabolism (5 genes), glutathione metabolism (4 genes), and ascorbate and aldarate metabolism (4 genes). The most abundantly down-regulated pathways were plant-pathogen interaction (14 genes), followed by phenylpropanoid biosynthesis (7 genes); arginine and proline metabolism (6 genes); cutin, suberine, and wax biosynthesis (4 genes); and tryptophan metabolism (4 genes). Several calcium sensing/transducing proteins, as well as transcription factors associated with abscisic acid (ABA), salicylic acid (SA), auxin, and cytokinin hormones were either up- or down-regulated in GB-primed B. rapa plants under heat stress. In particular, expression of the genes for antioxidant defense, heat shock response, and DNA damage repair systems were highly increased by GB priming. On the other hand, many of the genes involved in the calcium sensors and cell surface receptors involved in plant innate immunity and the biosynthesis of secondary metabolites were down-regulated in the absence of pathogen elicitors in GB-primed B. rapa seedlings. Overall GB priming activated ABA and SA signaling pathways but deactivated auxin and cytokinin signaling pathways while suppressing the innate immunity in B. rapa seedlings exposed to heat stress. The present study provides a preliminary understanding of the thermotolerance mechanisms in GB-primed plants and is of great importance in developing thermotolerant B. rapa cultivars by using the identified DEGs through genetic modification.

Transcriptome analysis identifies differentially expressed genes involved in lignin biosynthesis in barley

Lignin is an essential metabolite for plant growth but negatively affects the quality of forage barley. Genetic modification of quality traits to improve the forage digestibility requires an understanding of the molecular mechanism of lignin biosynthesis. RNA-Seq was used to quantify transcripts differentially expressed among leaf, stem and spike tissues from two barley genotypes. A total of 13,172 differentially expressed genes (DEGs) were identified, of which much more up-regulated DEGs were detected from the contrasting groups of leaf vs spike (L-S) and stem vs spike (S-S), and down-regulated DEGs were dominant in the group of stem vs leaf (S-L). 47 DEGs were successfully annotated to the monolignol pathway and six of them were candidate genes regulating the lignin biosynthesis. The qRT-PCR assay verified the expression profiles of the six candidate genes. Among them, four genes might positively regulate the lignin biosynthesis during forage barley development in terms of the consistency of their expression levels and changes of lignin content among the tissues, while the other two genes may have the reverse effects. These findings provide target genes for further investigations on molecular regulatory mechanisms of lignin biosynthesis and genetic resources for improvement of forage quality in barley molecular breeding programme.

A new biotechnology for in-planta gene editing and its application in promoting flavonoid biosynthesis in bamboo leaves

BACKGROUND

Bamboo is a perennial and renewable biomass forest resource and its leaf flavonoid is an antioxidant for biological and pharmacological research. The established genetic transformation and gene editing systems in bamboo are significantly limited by the dependence on bamboo regeneration capability. The way to improve the flavonoid content in bamboo leaves through biotechnology is still not feasible.

RESULTS

Here, we developed an in-planta, Agrobacterium-mediated gene expression method for exogenous genes via wounding and vacuum in bamboo. We demonstrated that the RUBY served as a reporter efficiently expressed in bamboo leaves and shoots, albeit unable to integrate into the chromosome. We have also developed a gene editing system by creating an in situ mutant of the bamboo violaxanthin de-epoxidase (PeVDE) gene in bamboo leaves, with lower NPQ values under the fluorometer, which can serve as a native reporter for gene editing. Furthermore, the bamboo leaves with increased flavonoid content were achieved by knocking out the cinnamoyl-CoA reductase genes.

CONCLUSIONS

Our method can be applied for the functional characterization of novel genes in a short time and is helpful for bamboo leaf flavonoid biotechnology breeding in the future.

OsCCRL1 is Essential for Phenylpropanoid Metabolism in Rice Anthers

Phenylpropanoid metabolism and timely tapetal degradation are essential for anther and pollen development, but the underlying mechanisms are unclear. In the current study, to investigate this, we identified and analyzed the male-sterile mutant, osccrl1 (cinnamoyl coA reductase-like 1), which exhibited delayed tapetal programmed cell death (PCD) and defective mature pollen. Map-based cloning, genetic complementation, and gene knockout revealed that OsCCRL1 corresponds to the gene LOC_Os09g32020.2, a member of SDR (short-chain dehydrogenase/reductase) family enzyme. OsCCRL1 was preferentially expressed in the tapetal cells and microspores, and localized to the nucleus and cytoplasm in both rice protoplasts and Nicotiana benthamiana leaves. The osccrl1 mutant exhibited reduced CCRs enzyme activity, less lignin accumulation, delayed tapetum degradation, and disrupted phenylpropanoid metabolism. Furthermore, an R2R3 MYB transcription factor OsMYB103/OsMYB80/OsMS188/BM1, involved in tapetum and pollen development, regulates the expression of OsCCRL1. Finally, the osmyb103 osccrl1 double mutants, exhibited the same phenotype as the osmyb103 single mutant, further indicating that OsMYB103/OsMYB80/OsMS188/BM1 functions upstream of OsCCRL1. These findings help to clarify the role of phenylpropanoid metabolism in male sterility and the regulatory network underlying the tapetum degradation.

Molecular cloning and characterization of Cinnamoyl-CoA reductase promoter gene from Asarum sieboldii Miq

Asarum sieboldii Miq., a perennial herb of the family Aristolochiaceae, is widely used in China to treat cold, fever, aphthous stomatitis, toothache, gingivitis, and rheumatoid arthritis. Methyleugenol is the most representative pharmacological constituent of this medicinal herb. Cinnamoyl-CoA reductase (CCR), which has been well known for occupying a critical position in the lignin biosynthesis pathway, is also shared with the biosynthesis of methyleugenol. To better understand the regulatory mechanisms of methyleugenol biosynthesis, a 1530-bp long promoter region of the AsCCR1 gene was isolated. PLACE and PlantCARE analysis affirmed the existence of the core promoter elements such as TATA and CAAT boxes, abiotic stress-responsive cis-regulation elements like abscisic acid-responsive element, G-box, and MBS in the isolated sequence. The histochemical assay suggested that it was a constitutive promoter, highly expressed in the root tissue. Moreover, the region of -200 bp to ATG (start codon) was enough to drive the expression of It GUS gene. Treatments with low temperature and high concentration of gibberellin or abscisic acid demonstrated the abiotic stress-induced expression of the AsCCR1 promoter. Overall, this study revealed the isolation and characterization of the AsCCR1 promoter. Moreover, it also provided a candidate gene for molecular breeding in A. sieboldii to enhance its pharmacological potential.

Transcriptomics and metabolomics reveal tolerance new mechanism of rice roots to Al stress

The prevalence of soluble aluminum (Al) ions is one of the major limitations to crop production worldwide on acid soils. Therefore, understanding the Al tolerance mechanism of rice and applying Al tolerance functional genes in sensitive plants can significantly improve Al stress resistance. In this study, transcriptomics and metabolomics analyses were performed to reveal the mechanism of Al tolerance differences between two rice landraces (Al-tolerant genotype Shibanzhan (KR) and Al-sensitive genotype Hekedanuo (MR) with different Al tolerance. The results showed that DEG related to phenylpropanoid biosynthesis was highly enriched in KR and MR after Al stress, indicating that phenylpropanoid biosynthesis may be closely related to Al tolerance. E1.11.1.7 (peroxidase) was the most significant enzyme of phenylpropanoid biosynthesis in KR and MR under Al stress and is regulated by multiple genes. We further identified that two candidate genes Os02g0770800 and Os06g0521900 may be involved in the regulation of Al tolerance in rice. Our results not only reveal the resistance mechanism of rice to Al stress to some extent, but also provide a useful reference for the molecular mechanism of different effects of Al poisoning on plants.

Genome-wide identification, characterization, and genetic diversity of CCR gene family in Dalbergia odorifera

INTRODUCTION

Lignin is a complex aromatic polymer plays major biological roles in maintaining the structure of plants and in defending them against biotic and abiotic stresses. Cinnamoyl-CoA reductase (CCR) is the first enzyme in the lignin-specific biosynthetic pathway, catalyzing the conversion of hydroxycinnamoyl-CoA into hydroxy cinnamaldehyde. Dalbergia odorifera T. Chen is a rare rosewood species for furniture, crafts and medicine. However, the CCR family genes in D. odorifera have not been identified, and their function in lignin biosynthesis remain uncertain.

METHODS AND RESULTS

Here, a total of 24 genes, with their complete domains were identified. Detailed sequence characterization and multiple sequence alignment revealed that the DoCCR protein sequences were relatively conserved. They were divided into three subfamilies and were unevenly distributed on 10 chromosomes. Phylogenetic analysis showed that seven DoCCRs were grouped together with functionally characterized CCRs of dicotyledons involved in developmental lignification. Synteny analysis showed that segmental and tandem duplications were crucial in the expansion of CCR family in D. odorifera, and purifying selection emerged as the main force driving these genes evolution. Cis-acting elements in the putative promoter regions of DoCCRs were mainly associated with stress, light, hormones, and growth/development. Further, analysis of expression profiles from the RNA-seq data showed distinct expression patterns of DoCCRs among different tissues and organs, as well as in response to stem wounding. Additionally, 74 simple sequence repeats (SSRs) were identified within 19 DoCCRs, located in the intron or untranslated regions (UTRs), and mononucleotide predominated. A pair of primers with high polymorphism and good interspecific generality was successfully developed from these SSRs, and 7 alleles were amplified in 105 wild D. odorifera trees from 17 areas covering its whole native distribution.

DISCUSSION

Overall, this study provides a basis for further functional dissection of CCR gene families, as well as breeding improvement for wood properties and stress resistance in D. odorifera.

Regional association and transcriptome analysis revealed candidate genes controlling plant height in Brassica napus

Plant height is a key morphological trait in rapeseed, which not only plays an important role in determining plant architecture, but is also an important characteristic related to yield. Presently, the improvement of plant architecture is a major challenge in rapeseed breeding. This work was carried out to identify genetic loci related to plant height in rapeseed. In this study, a genome-wide association study (GWAS) of plant height was performed using a Brassica 60 K Illumina Infinium SNP array and 203 Brassica napus accessions. Eleven haplotypes containing important candidate genes were detected and significantly associated with plant height on chromosomes A02, A03, A05, A07, A08, C03, C06, and C09. Moreover, regional association analysis of 50 resequenced rapeseed inbred lines was used to further analyze these eleven haplotypes and revealed nucleotide variation in the BnFBR12-A08 and BnCCR1-C03 gene regions related to the phenotypic variation in plant height. Furthermore, coexpression network analysis showed that BnFBR12-A08 and BnCCR1-C03 were directly connected with hormone genes and transcription factors and formed a potential network regulating the plant height of rapeseed. Our results will aid in the development of haplotype functional markers to further improve plant height in rapeseed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01337-1.

Silicon-Induced Tolerance against Arsenic Toxicity by Activating Physiological, Anatomical and Biochemical Regulation in Phoenix dactylifera (Date Palm)

Arsenic is a toxic metal abundantly present in agricultural, industrial, and pesticide effluents. To overcome arsenic toxicity and ensure safety for plant growth, silicon (Si) can play a significant role in its mitigation. Here, we aim to investigate the influence of silicon on date palm under arsenic toxicity by screening antioxidants accumulation, hormonal modulation, and the expression profile of abiotic stress-related genes. The results showed that arsenic exposure (As: 1.0 mM) significantly retarded growth attributes (shoot length, root length, fresh weight), reduced photosynthetic pigments, and raised reactive species levels. Contrarily, exogenous application of Si (Na₂SiO₃) to date palm roots strongly influenced stress mitigation by limiting the translocation of arsenic into roots and shoots as compared with the arsenic sole application. Furthermore, an enhanced accumulation of polyphenols (48%) and increased antioxidant activities (POD: 50%, PPO: 75%, GSH: 26.1%, CAT: 51%) resulted in a significant decrease in superoxide anion (O₂^•−: 58%) and lipid peroxidation (MDA: 1.7-fold), in silicon-treated plants, compared with control and arsenic-treated plants. The Si application also reduced the endogenous abscisic acid (ABA: 38%) under normal conditions, and salicylic acid (SA: 52%) and jasmonic acid levels (JA: 62%) under stress conditions as compared with control and arsenic. Interestingly, the genes; zeaxanthin epoxidase (ZEP) and 9-cis-epoxycarotenoid dioxygenase (NCED-1) involved in ABA biosynthesis were upregulated by silicon under arsenic stress. Likewise, Si application also upregulated gene expression of plant plasma membrane ATPase (PMMA-4), aluminum-activated malate transporter (ALMT) responsible for maintaining cellular physiology, stomatal conductance, and short-chain dehydrogenases/reductases (SDR) involved in nutrients translocation. Hence, the study demonstrates the remarkable role of silicon in supporting growth and inducing arsenic tolerance by increasing antioxidant activities and endogenous hormones in date palm. The outcomes of our study can be employed in further studies to better understand arsenic tolerance and decode mechanism.

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.

Combined nature and human selections reshaped peach fruit metabolome

Background

Plant metabolites reshaped by nature and human beings are crucial for both their lives and human health. However, which metabolites respond most strongly to selection pressure at different evolutionary stages and what roles they undertake on perennial fruit crops such as peach remain unclear.

Results

Here, we report 18,052 significant locus-trait associations, 12,691 expression-metabolite correlations, and 294,676 expression quantitative trait loci (eQTLs) for peach. Our results indicate that amino acids accumulated in landraces may be involved in the environmental adaptation of peaches by responding to low temperature and drought. Moreover, the contents of flavonoids, the major nutrients in fruits, have kept decreasing accompanied by the reduced bitter flavor during both domestication and improvement stages. However, citric acid, under the selection of breeders’ and consumers’ preference for flavor, shows significantly different levels between eastern and western varieties. This correlates with differences in activity against cancer cells in vitro in fruit from these two regions. Based on the identified key genes regulating flavonoid and acid contents, we propose that more precise and targeted breeding technologies should be designed to improve peach varieties with rich functional contents because of the linkage of genes related to bitterness and acid taste, antioxidant and potential anti-cancer activity that are all located at the top of chromosome 5.

Conclusions

This study provides powerful data for future improvement of peach flavor, nutrition, and resistance in future and expands our understanding of the effects of natural and artificial selection on metabolites.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02719-6.

Transcriptome Mining Provides Insights into Cell Wall Metabolism and Fiber Lignification in Agave tequilana Weber

Resilience of growing in arid and semiarid regions and a high capacity of accumulating sugar-rich biomass with low lignin percentages have placed Agave species as an emerging bioenergy crop. Although transcriptome sequencing of fiber-producing agave species has been explored, molecular bases that control wall cell biogenesis and metabolism in agave species are still poorly understood. Here, through RNAseq data mining, we reconstructed the cellulose biosynthesis pathway and the phenylpropanoid route producing lignin monomers in A. tequilana, and evaluated their expression patterns in silico and experimentally. Most of the orthologs retrieved showed differential expression levels when they were analyzed in different tissues with contrasting cellulose and lignin accumulation. Phylogenetic and structural motif analyses of putative CESA and CAD proteins allowed to identify those potentially involved with secondary cell wall formation. RT-qPCR assays revealed enhanced expression levels of AtqCAD5 and AtqCESA7 in parenchyma cells associated with extraxylary fibers, suggesting a mechanism of formation of sclerenchyma fibers in Agave similar to that reported for xylem cells in model eudicots. Overall, our results provide a framework for understanding molecular bases underlying cell wall biogenesis in Agave species studying mechanisms involving in leaf fiber development in monocots.

Comparative transcriptome analysis of resistant and susceptible wheat in response to Rhizoctonia cerealis

BACKGROUND

Sheath blight is an important disease caused by Rhizoctonia cerealis that affects wheat yields worldwide. No wheat varieties have been identified with high resistance or immunity to sheath blight. Understanding the sheath blight resistance mechanism is essential for controlling this disease. In this study, we investigated the response of wheat to Rhizoctonia cerealis infection by analyzing the cytological changes and transcriptomes of common wheat 7182 with moderate sensitivity to sheath blight and H83 with moderate resistance.

RESULTS

The cytological observation showed that the growth of Rhizoctonia cerealis on the surface and its expansion inside the leaf sheath tissue were more rapid in the susceptible material. According to the transcriptome sequencing results, a total of 88685 genes were identified in both materials, including 20156 differentially expressed genes (DEGs) of which 12087 was upregulated genes and 8069 was downregulated genes. At 36 h post-inoculation, compared with the uninfected control, 11498 DEGs were identified in resistant materials, with 5064 downregulated genes and 6434 upregulated genes, and 13058 genes were detected in susceptible materials, with 6759 downregulated genes and 6299 upregulated genes. At 72 h post-inoculation, compared with the uninfected control, 6578 DEGs were detected in resistant materials, with 2991 downregulated genes and 3587 upregulated genes, and 7324 genes were detected in susceptible materials, with 4119 downregulated genes and 3205 upregulated genes. Functional annotation and enrichment analysis showed that the main pathways enriched for the DEGs included biosynthesis of secondary metabolites, carbon metabolism, plant hormone signal transduction, and plant-pathogen interaction. In particular, phenylpropane biosynthesis pathway is specifically activated in resistant variety H83 after infection. Many DEGs also belonged to the MYB, AP2, NAC, and WRKY transcription factor families.

CONCLUSIONS

Thus, we suggest that the normal functioning of plant signaling pathways and differences in the expression of key genes and transcription factors in some important metabolic pathways may be important for defending wheat against sheath blight. These findings may facilitate further exploration of the sheath blight resistance mechanism in wheat and the cloning of related genes.

Characterization of the Transcriptome and Proteome of Brassica napus Reveals the Close Relation between DW871 Dwarfing Phenotype and Stalk Tissue

Rapeseed is a significant oil-bearing cash crop. As a hybrid crop, Brassica napus L. produces a high yield, but it also has drawbacks such as a tall stalk, easy lodging, and is not suitable for mechanized production. To address these concerns, we created the DW871 rapeseed dwarf variety, which has a high yield, high oil content, and is suitable for mechanized production. To fully comprehend the dwarfing mechanism of DW871 and provide a theoretical foundation for future applications of the variety, we used transcriptome and proteome sequencing to identify genes and proteins associated with the dwarfing phenotype, using homologous high-stalk material HW871 as a control. By RNA-seq and iTRAQ, we discovered 8665 DEGs and 50 DAPs. Comprehensive transcription and translation level analysis revealed 25 correlations, 23 of which have the same expression trend, involving monolignin synthesis, pectin-lignin assembly, lignification, glucose modification, cell wall composition and architecture, cell morphology, vascular bundle development, and stalk tissue composition and architecture. As a result of these results, we can formulate a hypothesis about the DW871 dwarfing phenotype: plant hormone signal transduction, such as IAA and BRs, is linked to the formation of dwarf phenotypes, and metabolic pathways related to lignin synthesis, such as phenylpropane biosynthesis, also play a role. Our works will contribute to a better understanding of the genes and proteins involved in the rapeseed dwarf phenotype, and we will propose new insights into the dwarfing mechanism of Brassica napus L.

PtrMYB120 functions as a positive regulator of both anthocyanin and lignin biosynthetic pathway in a hybrid poplar

Both anthocyanins and lignins are essential secondary metabolites in plant growth and development. Their biosynthesis is metabolically interconnected and diverges in the central metabolite 4-coumaroyl CoA of the phenylpropanoid pathway. Considerable progress has been made in understanding transcriptional regulation of genes involved in lignin and anthocyanin synthesis pathways, but the concerted regulation of these pathways is not yet fully understood. Here, we functionally characterized PtrMYB120, a R2R3-MYB transcription factor from Populus trichocarpa. Overexpression of PtrMYB120 in a hybrid poplar (i.e., 35S::PtrMYB120) was associated with increased anthocyanin (i.e., cyanidin 3-O-glucoside) accumulation and upregulation of anthocyanin biosynthetic genes. However, transgenic poplars with dominant suppression of PtrMYB120 function achieved by fusing the ERF-associated amphiphilic repression motif to PtrMYB120 (i.e., 35S::PtrMYB120-SRDX) had a dramatic decrease in not only anthocyanin but also Klason lignin content with downregulation of both anthocyanin and lignin biosynthetic genes. Indeed, 35S::PtrMYB120-SRDX poplars had irregularly shaped xylem vessels with reduced S-lignin content in stems, which was proportionally related to the level of the introduced PtrMYB120-SRDX gene. Furthermore, protoplast-based transcriptional activation assay using the PtrMYB120-GR system suggested that PtrMYB120 directly regulates genes involved in both anthocyanin and lignin biosynthesis, including chalcone synthase and ferulate-5 hydroxylase. Interestingly, the saccharification efficiency of line #6 of 35S::PtrMYB120-SRDX poplars, which had slightly reduced lignin content with a normal growth phenotype, was dramatically enhanced (>45%) by NaOH treatment. Taken together, our results suggest that PtrMYB120 functions as a positive regulator of both anthocyanin and lignin biosynthetic pathways and can be targeted to enhance saccharification efficiency in woody perennials.

Integrating genomics and multiplatform metabolomics enables metabolite quantitative trait loci detection in breeding-relevant apple germplasm

Apple (Malus × domestica) has commercial and nutritional value, but breeding constraints of tree crops limit varietal improvement. Marker-assisted selection minimises these drawbacks, but breeders lack applications for targeting fruit phytochemicals. To understand genotype-phytochemical associations in apples, we have developed a high-throughput integration strategy for genomic and multiplatform metabolomics data. Here, 124 apple genotypes, including members of three pedigree-connected breeding families alongside diverse cultivars and wild selections, were genotyped and phenotyped. Metabolite genome-wide association studies (mGWAS) were conducted with c. 10 000 single nucleotide polymorphisms and phenotypic data acquired via LC-MS and ¹ H NMR untargeted metabolomics. Putative metabolite quantitative trait loci (mQTL) were then validated via pedigree-based analyses (PBA). Using our developed method, 519, 726 and 177 putative mQTL were detected in LC-MS positive and negative ionisation modes, and NMR, respectively. mQTL were indicated on each chromosome, with hotspots on linkage groups 16 and 17. A chlorogenic acid mQTL was discovered on chromosome 17 via mGWAS and validated with a two-step PBA, enabling discovery of novel candidate gene-metabolite relationships. Complementary data from three metabolomics approaches and dual genomics analyses increased confidence in validity of compound annotation and mQTL detection. Our platform demonstrates the utility of multiomic integration to advance data-driven, phytochemical-based plant breeding.

Beyond Purple Tomatoes: Combined Strategies Targeting Anthocyanins to Generate Crimson, Magenta, and Indigo Fruit

The range of colours of many flowers and fruits is largely due to variations in the types of anthocyanins produced. The degree of hydroxylation on the B-ring affects the hue of these pigments, causing a shift from the orange end of the visible spectrum to the blue end. Besides colour, this modification can also affect other properties of anthocyanins, including the ability to protect the plant against different stresses or, when included in the human diet, to provide benefits for disease prevention. The level of hydroxylation of the B-ring is determined by the activity of two key hydroxylases, F3′H and F3′5′H, and by the substrate preference of DFR, an enzyme acting downstream in the biosynthetic pathway. We show that, in tomato, a strategy based on fruit-specific engineering of three regulatory genes (AmDel, AmRos1, AtMYB12) and a single biosynthetic gene (AmDFR), together with the availability of a specific mutation (f3′5′h), results in the generation of three different varieties producing high levels of anthocyanins with different levels of hydroxylation. These tomatoes show distinctive colours and mimic the classes of anthocyanins found in natural berries, thus providing unique near-isogenic material for different studies.

Genome-Wide Identification and Characterization of Short-Chain Dehydrogenase/Reductase (SDR) Gene Family in Medicago truncatula

Short-chain dehydrogenase/reductase (SDR) belongs to the NAD(P)(H)-dependent oxidoreductase superfamily. Limited investigations reveal that SDRs participate in diverse metabolisms. A genome-wide identification of the SDR gene family in M. truncatula was conducted. A total of 213 MtSDR genes were identified, and they were distributed on all chromosomes unevenly. MtSDR proteins were categorized into seven subgroups based on phylogenetic analysis and three types including 'classic', 'extended', and 'atypical', depending on the cofactor-binding site and active site. Analysis of the data from M. truncatula Gene Expression Atlas (MtGEA) showed that above half of MtSDRs were expressed in at least one organ, and lots of MtSDRs had a preference in a tissue-specific expression. The cis-acting element responsive to plant hormones (salicylic acid, ABA, auxin, MeJA, and gibberellin) and stresses were found in the promoter of some MtSDRs. Many genes of MtSDR7C,MtSDR65C, MtSDR110C, MtSDR114C, and MtSDR108E families were responsive to drought, salt, and cold. The study provides useful information for further investigation on biological functions of MtSDRs, especially in abiotic stress adaptation, in the future.

The evolution of the phenylpropanoid pathway entailed pronounced radiations and divergences of enzyme families

Land plants constantly respond to fluctuations in their environment. Part of their response is the production of a diverse repertoire of specialized metabolites. One of the foremost sources for metabolites relevant to environmental responses is the phenylpropanoid pathway, which was long thought to be a land-plant-specific adaptation shaped by selective forces in the terrestrial habitat. Recent data have, however, revealed that streptophyte algae, the algal relatives of land plants, have candidates for the genetic toolkit for phenylpropanoid biosynthesis and produce phenylpropanoid-derived metabolites. Using phylogenetic and sequence analyses, we here show that the enzyme families that orchestrate pivotal steps in phenylpropanoid biosynthesis have independently undergone pronounced radiations and divergence in multiple lineages of major groups of land plants; sister to many of these radiated gene families are streptophyte algal candidates for these enzymes. These radiations suggest a high evolutionary versatility in the enzyme families involved in the phenylpropanoid-derived metabolism across embryophytes. We suggest that this versatility likely translates into functional divergence, and may explain the key to one of the defining traits of embryophytes: a rich specialized metabolism.

Association of gene expression with syringyl to guaiacyl ratio in sugarcane lignin

A transcriptome analysis reveals the transcripts and alleles differentially expressed in sugarcane genotypes with contrasting lignin composition. Sugarcane bagasse is a highly abundant resource that may be used as a feedstock for the production of biofuels and bioproducts in order to meet increasing demands for renewable replacements for fossil carbon. However, lignin imparts rigidity to the cell wall that impedes the efficient breakdown of the biomass into fermentable sugars. Altering the ratio of the lignin units, syringyl (S) and guaiacyl (G), which comprise the native lignin polymer in sugarcane, may facilitate the processing of bagasse. This study aimed to identify genes and markers associated with S/G ratio in order to accelerate the development of sugarcane bioenergy varieties with modified lignin composition. The transcriptome sequences of 12 sugarcane genotypes that contrasted for S/G ratio were compared and there were 2019 transcripts identified as differentially expressed (DE) between the high and low S/G ratio groups. These included transcripts encoding possible monolignol biosynthetic pathway enzymes, transporters, dirigent proteins and transcriptional and post-translational regulators. Furthermore, the frequencies of single nucleotide polymorphisms (SNPs) were compared between the low and high S/G ratio groups to identify specific alleles expressed with the phenotype. There were 2063 SNP loci across 787 unique transcripts that showed group-specific expression. Overall, the DE transcripts and SNP alleles identified in this study may be valuable for breeding sugarcane varieties with altered S/G ratio that may provide desirable bioenergy traits.

Ethylene application at the immature stage of Fragaria chiloensis fruit represses the anthocyanin biosynthesis with a concomitant accumulation of lignin

Ethylene seems to play a secondary role in non-climacteric strawberry ripening compared to abscisic acid. However, this does not exclude that ethylene can regulate some specific events related to the ripening process. Preliminary experiments of applications of ethylene or its inhibitor 1-MCP to strawberry fruits have reinforced this hypothesis. Here, we reveal some previously non-covered physiological effects of ethylene using an in vitro strawberry ripening system. Fruits of Fragaria chiloensis treated with ethephon at the large green developmental stage showed inhibition of anthocyanin biosynthesis and downregulation of essential anthocyanin biosynthesis genes during the ripening. At the same time, ethylene stimulated lignin biosynthesis and remarkably upregulated the expression of FcPOD27. Since contrasting results have been reported when ethylene was applied at late ripening developmental stages, our findings support the hypothesis of a temporal-specific ethylene role in the ripening of strawberry fruits.

Using Gene Expression to Study Specialized Metabolism-A Practical Guide

Plants produce a vast array of chemical compounds that we use as medicines and flavors, but these compounds' biosynthetic pathways are still poorly understood. This paucity precludes us from modifying, improving, and mass-producing these specialized metabolites in suitable bioreactors. Many of the specialized metabolites are expressed in a narrow range of organs, tissues, and cell types, suggesting a tight regulation of the responsible biosynthetic pathways. Fortunately, with unprecedented ease of generating gene expression data and with >200,000 publicly available RNA sequencing samples, we are now able to study the expression of genes from hundreds of plant species. This review demonstrates how gene expression can elucidate the biosynthetic pathways by mining organ-specific genes, gene expression clusters, and applying various types of co-expression analyses. To empower biologists to perform these analyses, we showcase these analyses using recently published, user-friendly tools. Finally, we analyze the performance of co-expression networks and show that they are a valuable addition to elucidating multiple the biosynthetic pathways of specialized metabolism.

Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions

Phenylpropanoid metabolism is one of the most important metabolisms in plants, yielding more than 8,000 metabolites contributing to plant development and plant-environment interplay. Phenylpropanoid metabolism materialized during the evolution of early freshwater algae that were initiating terrestrialization and land plants have evolved multiple branches of this pathway, which give rise to metabolites including lignin, flavonoids, lignans, phenylpropanoid esters, hydroxycinnamic acid amides, and sporopollenin. Recent studies have revealed that many factors participate in the regulation of phenylpropanoid metabolism, and modulate phenylpropanoid homeostasis when plants undergo successive developmental processes and are subjected to stressful environments. In this review, we summarize recent progress on elucidating the contribution of phenylpropanoid metabolism to the coordination of plant development and plant-environment interaction, and metabolic flux redirection among diverse metabolic routes. In addition, our review focuses on the regulation of phenylpropanoid metabolism at the transcriptional, post-transcriptional, post-translational, and epigenetic levels, and in response to phytohormones and biotic and abiotic stresses.

Drought stress triggers proteomic changes involving lignin, flavonoids and fatty acids in tea plants

Drought stress triggers a series of physiological and biochemical changes in tea plants. It is well known that flavonoids, lignin and long-chain fatty acids play important roles in drought resistance. However, changes in proteins related to these three metabolic pathways in tea plants under drought stress have not been reported. We analysed the proteomic profiles of tea plants by tandem mass tag and liquid chromatography-tandem mass spectrometry. A total of 4789 proteins were identified, of which 11 and 100 showed up- and downregulation, respectively. The proteins related to the biosynthesis of lignin, flavonoids and long-chain fatty acids, including phenylalanine ammonia lyase, cinnamoyl-CoA reductase, peroxidase, chalcone synthase, flavanone 3-hydroxylase, flavonol synthase, acetyl-CoA carboxylase 1,3-ketoacyl-CoA synthase 6 and 3-ketoacyl-CoA reductase 1, were downregulated. However, the contents of soluble proteins, malondialdehyde, total phenols, lignin and flavonoids in the tea plants increased. These results showed that tea plants might improve drought resistance by inhibiting the accumulation of synthases related to lignin, flavonoids and long-chain fatty acids. The proteomic spectrum of tea plants provides a scientific basis for studying the pathways related to lignin, flavonoid and long-chain fatty acid metabolism in response to drought stress.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model species.

Identification and Functional Characterization of PtoMYB055 Involved in the Regulation of the Lignin Biosynthesis Pathway in Populus tomentosa

Wood, which is mainly composed of lignified secondary cell wall, is the most abundant biomass in woody plants. Previous studies have revealed that R2R3-type MYB transcription factors are important regulators of the formation of the secondary cell wall in vascular plants. In this study, we isolated the R2R3-type MYB transcription factor gene PtoMYB055, which is mainly expressed in xylem and phloem tissue, from Populus tomentosa and demonstrate that PtoMYB055 is a key regulator of lignin biosynthesis. PtoMYB055 as a transcriptional activator is localized to the nucleus. Overexpression of PtoMYB055 upregulates expression of lignin biosynthetic genes in transgenic poplar plants, resulting in ectopic deposition of lignin in phloem tissue and an increase in thickness of the secondary cell wall. In sum, PtoMYB055 is a transcriptional activator that is involved in regulating lignin biosynthesis during the formation of the secondary cell wall in poplar.

Implementing the CRISPR/Cas9 Technology in Eucalyptus Hairy Roots Using Wood-Related Genes

Eucalypts are the most planted hardwoods worldwide. The availability of the Eucalyptus grandis genome highlighted many genes awaiting functional characterization, lagging behind because of the lack of efficient genetic transformation protocols. In order to efficiently generate knock-out mutants to study the function of eucalypts genes, we implemented the powerful CRISPR/Cas9 gene editing technology with the hairy roots transformation system. As proofs-of-concept, we targeted two wood-related genes: Cinnamoyl-CoA Reductase1 (CCR1), a key lignin biosynthetic gene and IAA9A an auxin dependent transcription factor of Aux/IAA family. Almost all transgenic hairy roots were edited but the allele-editing rates and spectra varied greatly depending on the gene targeted. Most edition events generated truncated proteins, the prevalent edition types were small deletions but large deletions were also quite frequent. By using a combination of FT-IR spectroscopy and multivariate analysis (partial least square analysis (PLS-DA)), we showed that the CCR1-edited lines, which were clearly separated from the controls. The most discriminant wave-numbers were attributed to lignin. Histochemical analyses further confirmed the decreased lignification and the presence of collapsed vessels in CCR1-edited lines, which are characteristics of CCR1 deficiency. Although the efficiency of editing could be improved, the method described here is already a powerful tool to functionally characterize eucalypts genes for both basic research and industry purposes.

Large-scale Identification and Time-course Quantification of Ubiquitylation Events During Maize Seedling De-etiolation

The ubiquitin system is crucial for the development and fitness of higher plants. De-etiolation, during which green plants initiate photomorphogenesis and establish autotrophy, is a dramatic and complicated process that is tightly regulated by a massive number of ubiquitylation/de-ubiquitylation events. Here we present site-specific quantitative proteomic data for the ubiquitylomes of de-etiolating seedling leaves of Zea mays L. (exposed to light for 1, 6, or 12 h) achieved through immunoprecipitation-based high-resolution mass spectrometry (MS). Through the integrated analysis of multiple ubiquitylomes, we identified and quantified 1926 unique ubiquitylation sites corresponding to 1053 proteins. We analyzed these sites and found five potential ubiquitylation motifs, KA, AXK, KXG, AK, and TK. Time-course studies revealed that the ubiquitylation levels of 214 sites corresponding to 173 proteins were highly correlated across two replicate MS experiments, and significant alterations in the ubiquitylation levels of 78 sites (fold change >1.5) were detected after de-etiolation for 12 h. The majority of the ubiquitylated sites we identified corresponded to substrates involved in protein and DNA metabolism, such as ribosomes and histones. Meanwhile, multiple ubiquitylation sites were detected in proteins whose functions reflect the major physiological changes that occur during plant de-etiolation, such as hormone synthesis/signaling proteins, key C4 photosynthetic enzymes, and light signaling proteins. This study on the ubiquitylome of the maize seedling leaf is the first attempt ever to study the ubiquitylome of a C4 plant and provides the proteomic basis for elucidating the role of ubiquitylation during plant de-etiolation.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.

Novel motif is capable of determining CCR and CCR-like proteins based on the divergence of CCRs in plants

Cinnamoyl-coenzyme A reductases (CCRs) have been reported as key enzymes involved in monolignol biosynthesis. In this study, a motif-aware workflow based on a new signature motif effectively distinguished CCRs from CCR-like proteins. The divergence of CCRs and CCR-like sequences in Populus tomentosa Carr, Panicum virgatum L, Oryza sativa L and Selaginella moellendorffii Hieron suggests that NWYCY is not efficient for CCR recognition. The novel motif H202(X)2K205 (CCR-SBM or CCR substrate binding motif) was introduced to distinguish between CCRs and CCR-like proteins. The site-directed mutant R205K in Os(I)CCR-like and H202 in PtoCCR7 resulted in the rescue and loss of activity, respectively, further validating the fact that CCR-SBM is critical for maintaining CCR activity. The molecular docking using feruloyl-cinnamoyl-coenzyme A (CoA) as the ligand and binary PhCCR-NADP structures as receptors indicated an interaction between H202 and K205 with CoA moiety. The genuine CCRs and CCR-like proteins from several angiosperms and gymnosperms were screened using a motif-aware workflow and were validated using a biochemical assay. Our results suggest that the motif-aware workflow is efficient and effective for the identification of CCRs and CCR-like proteins in land plants and can be used as a more accurate way of identifying genuine CCRs among land plants.

An Integrated Analysis of the Rice Transcriptome and Metabolome Reveals Root Growth Regulation Mechanisms in Response to Nitrogen Availability

Nitrogen is an essential nutrient for plant growth and basic metabolic processes. Root systems play an important role in the ability of plants to obtain nutrients from the soil, and are closely related to the growth and development of above-ground plants. Root morphology analysis showed that root growth was induced under low-nitrogen conditions and inhibited under high-nitrogen conditions. To better understand the molecular mechanisms and metabolic basis underlying the rice root response to nitrogen availability, an integrated analysis of the rice root transcriptome and metabolome under three environmental conditions (low-, control, and high-nitrogen conditions) was conducted. A total of 262 and 262 differentially level metabolites were identified under low- and high-nitrogen conditions, respectively. A total of 696 and 808 differentially expressed genes were identified under low- and high-nitrogen conditions, respectively. For both the differentially expressed genes and metabolites, KEGG pathway analysis indicated that amino acid metabolism, carbon and nitrogen metabolism, phenylpropanoid metabolism, and phytohormones’ signal transduction were significantly affected by nitrogen availability. Additionally, variable levels of 65 transcription factors (TFs) were identified in rice leaves exposed to high and low nitrogen, covering 22 TF families. These results also indicate that there is a significant difference in the transcriptional regulation mechanisms of rice roots between low and high nitrogen. In summary, our study provides new information for a further understanding of the response of rice roots to low-nitrogen and high-nitrogen conditions.

PdWND3A, a wood-associated NAC domain-containing protein, affects lignin biosynthesis and composition in Populus

BACKGROUND

Plant secondary cell wall is a renewable feedstock for biofuels and biomaterials production. Arabidopsis VASCULAR-RELATED NAC DOMAIN (VND) has been demonstrated to be a key transcription factor regulating secondary cell wall biosynthesis. However, less is known about its role in the woody species.

RESULTS

Here we report the functional characterization of Populus deltoides WOOD-ASSOCIATED NAC DOMAIN protein 3 (PdWND3A), a sequence homolog of Arabidopsis VND4 and VND5 that are members of transcription factor networks regulating secondary cell wall biosynthesis. PdWND3A was expressed at higher level in the xylem than in other tissues. The stem tissues of transgenic P. deltoides overexpressing PdWND3A (OXPdWND3A) contained more vessel cells than that of wild-type plants. Furthermore, lignin content and lignin monomer syringyl and guaiacyl (S/G) ratio were higher in OXPdWND3A transgenic plants than in wild-type plants. Consistent with these observations, the expression of FERULATE 5-HYDROXYLASE1 (F5H1), encoding an enzyme involved in the biosynthesis of sinapyl alcohol (S unit monolignol), was elevated in OXPdWND3A transgenic plants. Saccharification analysis indicated that the rate of sugar release was reduced in the transgenic plants. In addition, OXPdWND3A transgenic plants produced lower amounts of biomass than wild-type plants.

CONCLUSIONS

PdWND3A affects lignin biosynthesis and composition and negatively impacts sugar release and biomass production.

Phenotypic and molecular characterization of Hessian fly resistance in diploid wheat, Aegilops tauschii

BACKGROUND

The Hessian fly (Mayetiola destructor), belonging to the gall midge family (Cecidomyiidae), is a devastating pest of wheat (Triticum aestivum) causing significant yield losses. Despite identification and characterization of numerous Hessian fly-responsive genes and associated biological pathways involved in wheat defense against this dipteran pest, their functional validation has been challenging. This is largely attributed to the large genome, polyploidy, repetitive DNA, and limited genetic resources in hexaploid wheat. The diploid progenitor Aegilops tauschii, D-genome donor of modern-day hexaploid wheat, offers an ideal surrogate eliminating the need to target all three homeologous chromosomes (A, B and D) individually, and thereby making the functional validation of candidate Hessian fly-responsive genes plausible. Furthermore, the well-annotated sequence of Ae. tauschii genome and availability of genetic resources amenable to manipulations makes the functional assays less tedious and time-consuming. However, prior to utilization of this diploid genome for downstream studies, it is imperative to characterize its physical and molecular responses to Hessian fly.

RESULTS

In this study we screened five Ae. tauschii accessions for their response to the Hessian fly biotypes L and vH13. Two lines were identified that exhibited a homozygous resistance response to feeding by both Hessian fly biotypes. Studies using physical measurements and neutral red staining showed that the resistant Ae. tauschii accessions resembled hexaploid wheat in their phenotypic responses to Hessian fly, that included similarities in larval developmental stages, leaf and plant growth, and cell wall permeability. Furthermore, molecular responses, characterized by gene expression profiling using quantitative real-time PCR, in select resistant Ae. tauschii lines also revealed similarities with resistant hexaploid wheat.

CONCLUSIONS

Phenotypic and molecular characterization of Ae. tauschii to Hessian fly infestation revealed resistant accessions that shared similarities to hexaploid wheat. Resembling the resistant hexaploid wheat, the Ae. tauschii accessions mount an early defense strategy involving defense proteins including lectins, secondary metabolites and reactive oxygen species (ROS) radicals. Our results reveal the suitability of the diploid progenitor for use as an ideal tool for functional genomics research in deciphering the wheat-Hessian fly molecular interactions.

Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis

BACKGROUND

The content of stone cells and lignin is one of the key factors affecting the quality of pear fruit. In a previous study, we determined the developmental regularity of stone cells and lignin in 'Dangshan Su' pear fruit 15-145 days after pollination (DAP). However, the development of fruit stone cells and lignin before 15 DAP has not been heavily researched.

RESULTS

In this study, we found that primordial stone cells began to appear at 7 DAP and that the fruit had formed a large number of stone cells at 15 DAP. Subsequently, transcriptome sequencing was performed on fruits at 0, 7, and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. During the 7-15 DAP period, a large number of key enzyme genes essential for lignin biosynthesis are gradually up-regulated, and their expression pattern is consistent with the accumulation of lignin in this period. Further analysis found that the biosynthesis of S-type lignin in 'Dangshan Su' pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H, and CAD. We cloned PbCCR1, 2 and analysed their functions in Chinese white pear lignin biosynthesis. PbCCR1 and 2 have a degree of functional redundancy; both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

CONCLUSIONS

Our results revealed that 'Dangshan Su' pear began to form a large number of stone cells and produce lignin after 7 DAP and mainly accumulated materials from 0 to 7 DAP. PbCCR1 is mainly involved in the biosynthesis of lignin in 'Dangshan Su' pear and plays a positive role in lignin biosynthesis.

Integrative Analysis of the Core Fruit Lignification Toolbox in Pear Reveals Targets for Fruit Quality Bioengineering

Stone cell content is an important factor affecting pear fruit flavor. Lignin, a major component of pear stone cells, hinders the quality and value of commercial fruit. The completion of the Chinese white pear (Pyrus bretschneideri) genome sequence provides an opportunity to perform integrative analysis of the genes encoding the eleven protein families (i.e., PAL, C4H, 4CL, HCT, C3H, CSE, CCoAOMT, CCR, F5H, COMT, and CAD) in the phenylpropanoid pathway. Here, a systematic study based on expression patterns and phylogenetic analyses was performed to identify the members of each gene family potentially involved in the lignification in the Chinese white pear. The phylogenetic analysis suggested that 35 P. bretschneideri genes belong to bona fide lignification clade members. Compared to other plants, some multigene families are expanded by tandem gene duplication, such as HCT, C3H, COMT, and CCR. RNA sequencing was used to study the expression patterns of the genes in different tissues, including leaf, petal, bud, sepal, ovary, stem, and fruit. Eighteen genes presented a high expression in fruit, indicating that these genes may be involved in the biosynthesis of lignin in pear fruit. Similarly to what has been observed for Populus trichocarpa, a bimolecular fluorescence complementation (BiFC) experiment indicated that P. bretschneideri C3H and C4H might also interact with each other to regulate monolignol biosynthesis in P. bretschneideri, ultimately affecting the stone cell content in pear fruits. The identification of the major genes involved in lignin biosynthesis in pear fruits provides the basis for the development of strategies to improve fruit quality.

Identification and expression profiling of genes governing lignin biosynthesis in Casuarina equisetifolia L

Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.

Biochemical Characterization of the Rice Cinnamyl Alcohol Dehydrogenase Gene Family

Cinnamyl alcohol dehydrogenase (CAD) is involved in the final step of the phenylpropanod pathway, catalyzing the NADPH-dependent reduction of hydroxy-cinnamaldehydes into the corresponding alcohols. The rice genome contains twelve CAD and CAD-like genes, collectively called OsCADs. To elucidate the biochemical function of the OsCADs, OsCAD1, 2, 6, and 7, which are highly expressed in rice, were cloned from rice tissues. The cloned OsCADs were heterologously expressed in Escherichia coli as His-tag fusion proteins. The activity assay of the recombinant OsCADs showed that OsCAD2, 6, and 7 have CAD activity toward hydroxycinnamaldehydes, but OsCAD1 has no detectable catalytic activity. The kinetic parameters of the enzyme reactions demonstrated that OsCAD2 has the highest catalytic activity among the examined enzymes. This result agrees well with the finding that the Zn binding and NADPH binding motifs and the residues constituting the substrate binding pocket in bona fide plant CADs were fully conserved in OsCAD2. Although they have large variations in the residue for the substrate binding pocket, OsCAD6 and 7 catalyzed the reduction of hydroxycinnamaldehydes with a similar efficiency. Alignment of amino acid sequences showed that OsCAD1 lacks the GxxxxP motif for NADPH binding and has mismatches in residues important in the reduction process, which could be responsible for the loss of catalytic activity. OsCAD2 belongs to CAD Class I with bona fide CADs from other plant species and is constitutively expressed throughout the developmental stages of rice, with preferential expression in actively lignifying tissues such as the root, stem, and panicle, suggesting that it is mainly involved in developmental lignification in rice. The expression of OsCAD2 was also induced by biotic and abiotic stresses such as Xanthomonas oryzae pv. oryzae (Xoo) infection and UV-irradiation, suggesting that it plays a role in the defense response of rice, in addition to a bona fide role in developmental lignification. OsCAD6 and 7 belong in CAD Class II. Their expression is relatively lower than that of OsCAD2 and is confined to certain tissues, such as the leaf sheath, stem, and panicle. The expression of OsCAD6 was stimulated by Xoo infection and UV-irradiation. Thus OsCAD6 appears to be an inducible OsCAD that is likely involved in the defense response of rice against biotic and abiotic stresses.

Cell wall composition and lignin biosynthetic gene expression along a developmental gradient in an Australian sugarcane cultivar

Sugarcane bagasse is an abundant source of lignocellulosic material for bioethanol production. Utilisation of bagasse for biofuel production would be environmentally and economically beneficial, but the recalcitrance of lignin continues to provide a challenge. Further understanding of lignin production in specific cultivars will provide a basis for modification of genomes for the production of phenotypes with improved processing characteristics. Here we evaluated the expression profile of lignin biosynthetic genes and the cell wall composition along a developmental gradient in KQ228 sugarcane. The expression levels of nine lignin biosynthesis genes were quantified in five stem sections of increasing maturity and in root tissue. Two distinct expression patterns were seen. The first saw highest gene expression in the youngest tissue, with expression decreasing as tissue matured. The second pattern saw little to no change in transcription levels across the developmental gradient. Cell wall compositional analysis of the stem sections showed total lignin content to be significantly higher in more mature tissue than in the youngest section assessed. There were no changes in structural carbohydrates across developmental sections. These gene expression and cell wall compositional patterns can be used, along with other work in grasses, to inform biotechnological approaches to crop improvement for lignocellulosic biofuel production.