Comparative genomic analysis of the R2R3 MYB secondary cell wall regulators of Arabidopsis, poplar, rice, maize, and switchgrass

Grass lignin: biosynthesis, biological roles, and industrial applications

Lignin is a phenolic heteropolymer found in most terrestrial plants that contributes an essential role in plant growth, abiotic stress tolerance, and biotic stress resistance. Recent research in grass lignin biosynthesis has found differences compared to dicots such as Arabidopsis thaliana. For example, the prolific incorporation of hydroxycinnamic acids into grass secondary cell walls improve the structural integrity of vascular and structural elements via covalent crosslinking. Conversely, fundamental monolignol chemistry conserves the mechanisms of monolignol translocation and polymerization across the plant phylum. Emerging evidence suggests grass lignin compositions contribute to abiotic stress tolerance, and periods of biotic stress often alter cereal lignin compositions to hinder pathogenesis. This same recalcitrance also inhibits industrial valorization of plant biomass, making lignin alterations and reductions a prolific field of research. This review presents an update of grass lignin biosynthesis, translocation, and polymerization, highlights how lignified grass cell walls contribute to plant development and stress responses, and briefly addresses genetic engineering strategies that may benefit industrial applications.

Interpreting omics data with pathway enrichment analysis

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

Genome-Wide Identification Analysis of the R2R3-MYB Transcription Factor Family in Cymbidium sinense for Insights into Drought Stress Responses

Cymbidium sinense represents a distinctive Orchidaceae plant that is more tolerant than other terrestrial orchids. Studies have shown that many members of the MYB transcription factor (TF) family, especially the R2R3-MYB subfamily, are responsive to drought stress. This study identified 103 CsMYBs; phylogenetic analysis classified these genes into 22 subgroups with Arabidopsis thaliana. Structural analysis showed that most CsMYB genes contained the same motifs, three exons and two introns, and showed a helix-turn-helix 3D structure in each R repeat. However, the members of subgroup 22 contained only one exon and no intron. Collinear analysis revealed that C. sinense had more orthologous R2R3-MYB genes with wheat than A. thaliana and rice. Ka/Ks ratios indicated that most CsMYB genes were under purifying negative selection pressure. Cis-acting elements analysis revealed that drought-related elements were mainly focused on subgroups 4, 8, 18, 20, 21, and 22, and Mol015419 (S20) contained the most. The transcriptome analysis results showed that expression patterns of most CsMYB genes were upregulated in leaves in response to slight drought stress and downregulated in roots. Among them, members in S8 and S20 significantly responded to drought stress in C. sinense. In addition, S14 and S17 also participated in these responses, and nine genes were selected for the real-time reverse transcription quantitative PCR (RT-qPCR) experiment. The results were roughly consistent with the transcriptome. Our results, thus, provide an important contribution to understanding the role of CsMYBs in stress-related metabolic processes.

Downregulation of Brassica napus MYB69 (BnMYB69) increases biomass growth and disease susceptibility via remodeling phytohormone, chlorophyll, shikimate and lignin levels

MYB transcription factors are major actors regulating plant development and adaptability. Brassica napus is a staple oil crop and is hampered by lodging and diseases. Here, four B. napus MYB69 (BnMYB69s) genes were cloned and functionally characterized. They were dominantly expressed in stems during lignification. BnMYB69 RNA interference (BnMYB69i) plants showed considerable changes in morphology, anatomy, metabolism and gene expression. Stem diameter, leaves, roots and total biomass were distinctly larger, but plant height was significantly reduced. Contents of lignin, cellulose and protopectin in stems were significantly reduced, accompanied with decrease in bending resistance and Sclerotinia sclerotiorum resistance. Anatomical detection observed perturbation in vascular and fiber differentiation in stems, but promotion in parenchyma growth, accompanied with changes in cell size and cell number. In shoots, contents of IAA, shikimates and proanthocyanidin were reduced, while contents of ABA, BL and leaf chlorophyll were increased. qRT-PCR revealed changes in multiple pathways of primary and secondary metabolisms. IAA treatment could recover many phenotypes and metabolisms of BnMYB69i plants. However, roots showed trends opposite to shoots in most cases, and BnMYB69i phenotypes were light-sensitive. Conclusively, BnMYB69s might be light-regulated positive regulators of shikimates-related metabolisms, and exert profound influences on various internal and external plant traits.

Transcriptome and metabolite profiling to identify genes associated with rhizome lignification and the function of ZoCSE in ginger (Zingiber officinale)

Ginger (Zingiber officinale Roscoe) is an important spice crop in China, and fresh ginger rhizomes are consumed as vegetable in Sichuan and Chongqing. However, tissue lignification accelerates with rhizome maturation, resulting in the loss of edible quality. To understand the molecular mechanisms of texture modification during rhizome development, we investigated lignin accumulation patterns and identified the key genes associated with lignin biosynthesis using gas chromatography-mass spectrometry (GC-MS), liquid chromatography-tandem mass spectrometry (LC-MS/MS) and RNA-sequencing (RNA-Seq). Results showed that the contents of total lignin and its precursors exhibited notable declines with tissue maturation. However, the lignin composition was remarkably modified and syringyl lignin was deposited in mature rhizomes, leading to ginger lignification. Transcriptome analysis displayed 32 lignin biosynthetic genes were dramatically downregulated with rhizome development, including caffeoylshikimate esterase (CSE ), 4-coumarate-CoA ligase , laccase , cinnamoyl-CoA reductase , cinnamyl-alcohol dehydrogenase , peroxidase and caffeic acid 3-O-methyltransferase , indicating that lignin reduction might be attributed to deficiency in intermediates or the downregulation of key biosynthetic enzymes. Furthermore, overexpressing ZoCSE in Nicotiana benthamiana L. enhanced the total lignin content, suggesting its fundamental role in lignin biosynthesis. RNA-Seq also identified candidate lignin production regulators, including hormone-related genes and NAC/MYB transcription factors (ZoNAC1 , ZoNAC4 , ZoMYB14 and ZoMYB17 ). This result provides a molecular basis for lignin accumulation in ginger.

Co-expression network analyses of anthocyanin biosynthesis genes in Ruellia (Wild Petunias; Acanthaceae)

Background

Anthocyanins are major pigments contributing to flower coloration and as such knowledge of molecular architecture underlying the anthocyanin biosynthetic pathway (ABP) is key to understanding flower color diversification. To identify ABP structural genes and associated regulatory networks, we sequenced 16 transcriptomes generated from 10 species of Ruellia and then conducted co-expression analyses among resulting data.

Results

Complete coding sequences for 12 candidate structural loci representing eight genes plus nine candidate regulatory loci were assembled. Analysis of non-synonymous/synonymous (dn/ds) mutation rates indicated all identified loci are under purifying selection, suggesting overall selection to prevent the accumulation of deleterious mutations. Additionally, upstream enzymes have lower rates of molecular evolution compared to downstream enzymes. However, site-specific tests of selection yielded evidence for positive selection at several sites, including four in F3'H2 and five in DFR3, and these sites are located in protein binding regions. A species-level phylogenetic tree constructed using a newly implemented hybrid transcriptome–RADseq approach implicates several flower color transitions among the 10 species. We found evidence of both regulatory and structural mutations to F3′5'H in helping to explain the evolution of red flowers from purple-flowered ancestors.

Conclusions

Sequence comparisons and co-expression analyses of ABP loci revealed that mutations in regulatory loci are likely to play a greater role in flower color transitions in Ruellia compared to mutations in underlying structural genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01955-x.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Identification of candidate MYB transcription factors that influence CslF6 expression in barley grain

(1,3;1,4)-β-Glucan is a non-cellulosic polysaccharide required for correct barley grain fill and plant development, with industrial relevance in the brewing and the functional food sector. Barley grains contain higher levels of (1,3;1,4)-β-glucan compared to other small grain cereals and this influences their end use, having undesirable effects on brewing and distilling and beneficial effects linked to human health. HvCslF6 is the main gene contributing to (1,3;1,4)-β-glucan biosynthesis in the grain. Here, the transcriptional regulation of HvCslF6 was investigated using an in-silico analysis of transcription factor binding sites (TFBS) in its putative promoter, and functional characterization in a barley protoplast transient expression system. Based on TFBS predictions, TF classes AP2/ERF, MYB, and basic helix-loop-helix (bHLH) were over-represented within a 1,000 bp proximal HvCslF6 promoter region. Dual luciferase assays based on multiple HvCslF6 deletion constructs revealed the promoter fragment driving HvCslF6 expression. Highest HvCslF6 promoter activity was narrowed down to a 51 bp region located -331 bp to -382 bp upstream of the start codon. We combined this with TFBS predictions to identify two MYB TFs: HvMYB61 and HvMYB46/83 as putative activators of HvCslF6 expression. Gene network analyses assigned HvMYB61 to the same co-expression module as HvCslF6 and other primary cellulose synthases (HvCesA1, HvCesA2, and HvCesA6), whereas HvMYB46/83 was assigned to a different module. Based on RNA-seq expression during grain development, HvMYB61 was cloned and tested in the protoplast system. The transient over-expression of HvMYB61 in barley protoplasts suggested a positive regulatory effect on HvCslF6 expression.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Distinct and Overlapping Functions of Miscanthus sinensis MYB Transcription Factors SCM1 and MYB103 in Lignin Biosynthesis

Cell wall recalcitrance is a major constraint for the exploitation of lignocellulosic biomass as a renewable resource for energy and bio-based products. Transcriptional regulators of the lignin biosynthetic pathway represent promising targets for tailoring lignin content and composition in plant secondary cell walls. However, knowledge about the transcriptional regulation of lignin biosynthesis in lignocellulosic feedstocks, such as Miscanthus, is limited. In Miscanthus leaves, MsSCM1 and MsMYB103 are expressed at growth stages associated with lignification. The ectopic expression of MsSCM1 and MsMYB103 in N. benthamiana leaves was sufficient to trigger secondary cell wall deposition with distinct sugar and lignin compositions. Moreover, RNA-seq analysis revealed that the transcriptional responses to MsSCM1 and MsMYB103 overexpression showed an extensive overlap with the response to the NAC master transcription factor MsSND1, but were distinct from each other, underscoring the inherent complexity of secondary cell wall formation. Furthermore, conserved and previously described promoter elements as well as novel and specific motifs could be identified from the target genes of the three transcription factors. Together, MsSCM1 and MsMYB103 represent interesting targets for manipulations of lignin content and composition in Miscanthus towards a tailored biomass.

Phylogenetic Occurrence of the Phenylpropanoid Pathway and Lignin Biosynthesis in Plants

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

AP2/ERF and R2R3-MYB family transcription factors: potential associations between temperature stress and lipid metabolism in Auxenochlorella protothecoides

BACKGROUND

Both APETALA2/Ethylene Responsive Factor (AP2/ERF) superfamily and R2R3-MYB family were from one of the largest diverse families of transcription factors (TFs) in plants, and played important roles in plant development and responses to various stresses. However, no systematic analysis of these TFs had been conducted in the green algae A. protothecoides heretofore. Temperature was a critical factor affecting growth and lipid metabolism of A. protothecoides. It also remained largely unknown whether these TFs would respond to temperature stress and be involved in controlling lipid metabolism process.

RESULTS

Hereby, a total of six AP2 TFs, six ERF TFs and six R2R3-MYB TFs were identified and their expression profiles were also analyzed under low-temperature (LT) and high-temperature (HT) stresses. Meanwhile, differential adjustments of lipid pathways were triggered, with enhanced triacylglycerol accumulation. A co-expression network was built between these 18 TFs and 32 lipid-metabolism-related genes, suggesting intrinsic associations between TFs and the regulatory mechanism of lipid metabolism.

CONCLUSIONS

This study represented an important first step towards identifying functions and roles of AP2 superfamily and R2R3-MYB family in lipid adjustments and response to temperature stress. These findings would facilitate the biotechnological development in microalgae-based biofuel production and the better understanding of photosynthetic organisms' adaptive mechanism to temperature stress.

Overexpression of PtrMYB121 Positively Regulates the Formation of Secondary Cell Wall in Arabidopsis thaliana

MYB transcription factors have a wide range of functions in plant growth, hormone signaling, salt, and drought tolerances. In this study, two homologous transcription factors, PtrMYB55 and PtrMYB121, were isolated and their functions were elucidated. Tissue expression analysis revealed that PtrMYB55 and PtrMYB121 had a similar expression pattern, which had the highest expression in stems. Their expression continuously increased with the growth of poplar, and the expression of PtrMYB121 was significantly upregulated in the process. The full length of PtrMYB121 was 1395 bp, and encoded protein contained 464 amino acids including conserved R2 and R3 MYB domains. We overexpressed PtrMYB121 in Arabidopsis thaliana, and the transgenic lines had the wider xylem as compared with wild-type Arabidopsis. The contents of cellulose and lignin were obviously higher than those in wild-type materials, but there was no significant change in hemicellulose. Quantitative real-time PCR demonstrated that the key enzyme genes regulating the synthesis of lignin and cellulose were significantly upregulated in the transgenic lines. Furthermore, the effector-reporter experiment confirmed that PtrMYB121 bound directly to the promoters of genes relating to the synthesis of lignin and cellulose. These results suggest that PtrMYB121 may positively regulate the formation of secondary cell wall by promoting the synthesis of lignin and cellulose.

Grass secondary cell walls, Brachypodium distachyon as a model for discovery

A key aspect of plant growth is the synthesis and deposition of cell walls. In specific tissues and cell types including xylem and fibre, a thick secondary wall comprised of cellulose, hemicellulose and lignin is deposited. Secondary cell walls provide a physical barrier that protects plants from pathogens, promotes tolerance to abiotic stresses and fortifies cells to withstand the forces associated with water transport and the physical weight of plant structures. Grasses have numerous cell wall features that are distinct from eudicots and other plants. Study of the model species Brachypodium distachyon as well as other grasses has revealed numerous features of the grass cell wall. These include the characterisation of xylosyl and arabinosyltransferases, a mixed-linkage glucan synthase and hydroxycinnamate acyltransferases. Perhaps the most fertile area for discovery has been the formation of lignins, including the identification of novel substrates and enzyme activities towards the synthesis of monolignols. Other enzymes function as polymerising agents or transferases that modify lignins and facilitate interactions with polysaccharides. The regulatory aspects of cell wall biosynthesis are largely overlapping with those of eudicots, but salient differences among species have been resolved that begin to identify the determinants that define grass cell walls.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

The MYB transcription factor CiMYB42 regulates limonoids biosynthesis in citrus

BACKGROUND

Limonoids are major bioactive compounds that are produced by the triterpenoid metabolic pathway. The detailed biochemical process of limonoid biosynthesis and the mechanism of its molecular regulation remain elusive. The identification of transcription factors that regulate limonoid biosynthetic pathways is very important for understanding the underlying regulatory mechanisms. This information could also provide tools for manipulating biosynthesis genes to modulate limonoid production.

RESULTS

In this study, the CiMYB42 transcription factor was isolated to identify its role in limonoid biosynthesis. Multiple alignment analysis and phylogenetic analysis demonstrated that CiMYB42 is a typical R2R3MYB transcription factor that shares high similarity of its amino acid sequence with AtMYB42. Limonoids contents were higher in Citrus sinensis and Citrus grandis than in other species. Limonoid accumulation during leaf development also showed diverse trends in different genotypes. The expression of CiMYB42 was significantly related to the limonoid content and the expression of CiOSC in some citrus accessions. The overexpression of CiMYB42 in sweet orange resulted in significant accumulation of limonin, whereas the downregulation of CiMYB42 by RNAi resulted in a dwarf phenotype and less nomilin accumulation. Furthermore, the results of a yeast one-hybrid assay and EMSA indicated that CiMYB42 binds exclusively to the TTGTTG sequence (type II MYB core) in the promoter of CiOSC. Together, these results suggest that CiMYB42 positively regulates limonoid biosynthesis by regulating the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core) of its promoter.

CONCLUSIONS

CiMYB42 is an important transcription activator involved in limonoid biosynthesis that regulates the expression of CiOSC by binding to the TTGTTG sequence (type II MYB core).

Overexpression of OsMYB305 in Rice Enhances the Nitrogen Uptake Under Low-Nitrogen Condition

Excessive nitrogen fertilizer application causes severe environmental degradation and drives up agricultural production costs. Thus, improving crop nitrogen use efficiency (NUE) is essential for the development of sustainable agriculture. Here, we characterized the roles of the MYB transcription factor OsMYB305 in nitrogen uptake and assimilation in rice. OsMYB305 encoded a transcriptional activator and its expression was induced by N deficiency in rice root. Under low-N condition, OsMYB305 overexpression significantly increased the tiller number, shoot dry weight and total N concentration. In the roots of OsMYB305-OE rice lines, the expression of OsNRT2.1, OsNRT2.2, OsNAR2.1, and OsNiR2 was up-regulated and 15NO3 - influx was significantly increased. In contrast, the expression of lignocellulose biosynthesis-related genes was repressed so that cellulose content decreased, and soluble sugar concentration increased. Certain intermediates in the glycolytic pathway and the tricarboxylic acid cycle were significantly altered and NADH-GOGAT, Pyr-K, and G6PDH were markedly elevated in the roots of OsMYB305-OE rice lines grown under low-N condition. Our results revealed that OsMYB305 overexpression suppressed cellulose biosynthesis under low-nitrogen condition, thereby freeing up carbohydrate for nitrate uptake and assimilation and enhancing rice growth. OsMYB305 is a potential molecular target for increasing NUE in rice.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

KNAT2/6b, a class I KNOX gene, impedes xylem differentiation by regulating NAC domain transcription factors in poplar

Wood formation is the terminal differentiation of xylem mother cells derived from cambial initials, and negative regulators play important roles in xylem differentiation. The molecular mechanism of the negative regulator of xylem differentiation PagKNAT2/6b was investigated. PagKNAT2/6b is an ortholog of Arabidopsis KNAT2 and KNAT6 that is highly expressed in phloem and xylem. Compared to nontransgenic control plants, transgenic poplar plants overexpressing PagKNAT2/6b present with altered vascular patterns, characterized by decreased secondary xylem with thin cell walls containing less cellulose, xylose and lignin. RNA sequencing analyses revealed that differentially expressed genes are enriched in xylem differentiation and secondary wall synthesis functions. Expression of NAM/ATAF/CUC (NAC) domain genes including PagSND1-A1, PagSND1-A2, PagSND1-B2 and PagVND6-C1 is downregulated by PagKNAT2/6b, while PagXND1a is directly upregulated. Accordingly, the dominant repression form of PagKNAT2/6b leads to increased xylem width per stem diameter through downregulation of PagXND1a. PagKNAT2/6b can inhibit cell differentiation and secondary wall deposition during wood formation in poplar by modulating the expression of NAC domain transcription factors. Direct activation of PagXND1a by PagKNAT2/6b is a key node in the negative regulatory network of xylem differentiation by KNOXs.

Silencing Folylpolyglutamate Synthetase1 (FPGS1) in Switchgrass (Panicum virgatum L.) Improves Lignocellulosic Biofuel Production

Switchgrass (Panicum virgatum L.) is a lignocellulosic perennial grass with great potential in bioenergy field. Lignocellulosic bioenergy crops are mostly resistant to cell wall deconstruction, and therefore yield suboptimal levels of biofuel. The one-carbon pathway (also known as C1 metabolism) is critical for polymer methylation, including that of lignin and hemicelluloses in cell walls. Folylpolyglutamate synthetase (FPGS) catalyzes a biochemical reaction that leads to the formation of folylpolyglutamate, an important cofactor for many enzymes in the C1 pathway. In this study, the putatively novel switchgrass PvFPGS1 gene was identified and its functional role in cell wall composition and biofuel production was examined by RNAi knockdown analysis. The PvFPGS1-downregulated plants were analyzed in the field over three growing seasons. Transgenic plants with the highest reduction in PvFPGS1 expression grew slower and produced lower end-of-season biomass. Transgenic plants with low-to-moderate reduction in PvFPGS1 transcript levels produced equivalent biomass as controls. There were no significant differences observed for lignin content and syringyl/guaiacyl lignin monomer ratio in the low-to-moderately reduced PvFPGS1 transgenic lines compared with the controls. Similarly, sugar release efficiency was also not significantly different in these transgenic lines compared with the control lines. However, transgenic plants produced up to 18% more ethanol while maintaining congruent growth and biomass as non-transgenic controls. Severity of rust disease among transgenic and control lines were not different during the time course of the field experiments. Altogether, the unchanged lignin content and composition in the low-to-moderate PvFPGS1-downregulated lines may suggest that partial downregulation of PvFPGS1 expression did not impact lignin biosynthesis in switchgrass. In conclusion, the manipulation of PvFPGS1 expression in bioenergy crops may be useful to increase biofuel potential with no growth penalty or increased susceptibility to rust in feedstock.

Induction of PrMADS10 on the lower side of bent pine tree stems: potential role in modifying plant cell wall properties and wood anatomy

The molecular mechanisms underlying inclination responses in trees are unclear. In this study, we identified a MADS-box transcription factor differentially expressed early after inclination in the stems of Pinus radiata D. Don. PrMADS10 has a CDS of 582 bp and encodes a group II MADS-box transcription factor. We measured highest accumulation of this transcript on the lower side of inclined pine stems. In an effort to identify putative targets, we stably transformed Arabidopsis thaliana with a 35S::PrMADS10 construct. Transcriptome analysis revealed 1,219 genes differentially-expressed, with 690 and 529 genes up- and down-regulated respectively, when comparing the transgenic and wild-type. Differentially-expressed genes belong to different biological processes, but were enriched in cell wall remodeling and phenylpropanoid metabolic functions. Interestingly, lignin content was 30% higher in transgenic as compared to wild-type plants consistent with observed changes in gene expression. Differentially expressed transcription factors and phenylpropanoid genes were analyzed using STRING. Several MYB and NAC transcription factors showed interactions with genes of the phenylpropanoid pathway. Together, these results implicate PrMADS10 as a regulatory factor, triggering the expression of other transcription factors and genes involved in the synthesis of lignin.

OsMYB108 loss-of-function enriches p-coumaroylated and tricin lignin units in rice cell walls

Breeding approaches to enrich lignins in biomass could be beneficial to improving the biorefinery process because lignins increase biomass heating value and represent a potent source of valuable aromatic chemicals. However, despite the fact that grasses are promising lignocellulose feedstocks, limited information is yet available for molecular-breeding approaches to upregulate lignin biosynthesis in grass species. In this study, we generated lignin-enriched transgenic rice (Oryza sativa), a model grass species, via targeted mutagenesis of the transcriptional repressor OsMYB108 using CRISPR/Cas9-mediated genome editing. The OsMYB108-knockout rice mutants displayed increased expressions of lignin biosynthetic genes and enhanced lignin deposition in culm cell walls. Chemical and two-dimensional nuclear magnetic resonance (NMR) analyses revealed that the mutant cell walls were preferentially enriched in γ-p-coumaroylated and tricin lignin units, both of which are typical and unique components in grass lignins. NMR analysis also showed that the relative abundances of major lignin linkage types were altered in the OsMYB108 mutants.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

Genome-Wide Identification of R2R3-MYB Transcription Factors Regulating Secondary Cell Wall Thickening in Cotton Fiber Development

MYB proteins represent one of the largest transcription factor (TF) families in plants, some of which act as key transcriptional regulators of secondary cell wall (SCW) biosynthesis. Cotton (Gossypium hirsutum) fiber is thought to be an ideal single-cell model to study cell elongation and SCW biosynthesis. However, little knowledge regarding the TFs controlling fiber SCW biosynthesis, particularly for R2R3-MYBs is known. By far, no comprehensive genome-wide analysis of the secondary wall-associated R2R3-MYBs has been reported in cultivated tetraploid upland cotton. In this study, we identified 419 R2R3-MYB genes by systematically examining the cotton genome. A combination of phylogenetic, RNA-seq and co-expression analyses indicated that 36 R2R3-MYBs were either preferentially or highly expressed in 20 day post anthesis (dpa) fibers and are putative SCW regulators. Among these MYB genes, 22 MYBs are homologs of known SCW MYB proteins and the other 14 MYBs are novel proteins without prior reported SCW biosynthesis-related functions. Finally, we highlighted on the roles of two MYBs named GhMYB46_D13 and GhMYB46_D9, both of which displayed the highest expression in 20 dpa fibers. Expression of GhMYB46_D13 or GhMYB46_D9 individually in Arabidopsis resulted in ectopic SCW deposition in transgenic plants. Furthermore, both GhMYB46_D13 and GhMYB46_D9 were able to activate the cotton fiber SCW cellulose synthase gene promoters. Thus, we have identified 36 R2R3-MYBs as potential SCW regulators in cotton fibers that represent strong candidates for further functional studies during fiber development and SCW thickening.

Rice Genome-Scale Network Integration Reveals Transcriptional Regulators of Grass Cell Wall Synthesis

Grasses have evolved distinct cell wall composition and patterning relative to dicotyledonous plants. However, despite the importance of this plant family, transcriptional regulation of its cell wall biosynthesis is poorly understood. To identify grass cell wall-associated transcription factors, we constructed the Rice Combined mutual Ranked Network (RCRN). The RCRN covers >90% of annotated rice (Oryza sativa) genes, is high quality, and includes most grass-specific cell wall genes, such as mixed-linkage glucan synthases and hydroxycinnamoyl acyltransferases. Comparing the RCRN and an equivalent Arabidopsis network suggests that grass orthologs of most genetically verified eudicot cell wall regulators also control this process in grasses, but some transcription factors vary significantly in network connectivity between these divergent species. Reverse genetics, yeast-one-hybrid, and protoplast-based assays reveal that OsMYB61a activates a grass-specific acyltransferase promoter, which confirms network predictions and supports grass-specific cell wall synthesis genes being incorporated into conserved regulatory circuits. In addition, 10 of 15 tested transcription factors, including six novel Wall-Associated regulators (WAP1, WACH1, WAHL1, WADH1, OsMYB13a, and OsMYB13b), alter abundance of cell wall-related transcripts when transiently expressed. The results highlight the quality of the RCRN for examining rice biology, provide insight into the evolution of cell wall regulation, and identify network nodes and edges that are possible leads for improving cell wall composition.

SECONDARY WALL ASSOCIATED MYB1 is a positive regulator of secondary cell wall thickening in Brachypodium distachyon and is not found in the Brassicaceae

Grass biomass is comprised chiefly of secondary walls that surround fiber and xylem cells. A regulatory network of interacting transcription factors in part regulates cell wall thickening. We identified Brachypodium distachyon SECONDARY WALL ASSOCIATED MYB1 (SWAM1) as a potential regulator of secondary cell wall biosynthesis based on gene expression, phylogeny, and transgenic plant phenotypes. SWAM1 interacts with cellulose and lignin gene promoters with preferential binding to AC-rich sequence motifs commonly found in the promoters of cell wall-related genes. SWAM1 overexpression (SWAM-OE) lines had greater above-ground biomass with only a slight change in flowering time while SWAM1 dominant repressor (SWAM1-DR) plants were severely dwarfed with a striking reduction in lignin of sclerenchyma fibers and stem epidermal cell length. Cellulose, hemicellulose, and lignin genes were significantly down-regulated in SWAM1-DR plants and up-regulated in SWAM1-OE plants. There was no reduction in bioconversion yield in SWAM1-OE lines; however, it was significantly increased for SWAM1-DR samples. Phylogenetic and syntenic analyses strongly suggest that the SWAM1 clade was present in the last common ancestor between eudicots and grasses, but is not in the Brassicaceae. Collectively, these data suggest that SWAM1 is a transcriptional activator of secondary cell wall thickening and biomass accumulation in B. distachyon.

A genome-wide analysis of the flax (Linum usitatissimum L.) dirigent protein family: from gene identification and evolution to differential regulation

Identification of DIR encoding genes in flax genome. Analysis of phylogeny, gene/protein structures and evolution. Identification of new conserved motifs linked to biochemical functions. Investigation of spatio-temporal gene expression and response to stress. Dirigent proteins (DIRs) were discovered during 8-8' lignan biosynthesis studies, through identification of stereoselective coupling to afford either (+)- or (-)-pinoresinols from E-coniferyl alcohol. DIRs are also involved or potentially involved in terpenoid, allyl/propenyl phenol lignan, pterocarpan and lignin biosynthesis. DIRs have very large multigene families in different vascular plants including flax, with most still of unknown function. DIR studies typically focus on a small subset of genes and identification of biochemical/physiological functions. Herein, a genome-wide analysis and characterization of the predicted flax DIR 44-membered multigene family was performed, this species being a rich natural grain source of 8-8' linked secoisolariciresinol-derived lignan oligomers. All predicted DIR sequences, including their promoters, were analyzed together with their public gene expression datasets. Expression patterns of selected DIRs were examined using qPCR, as well as through clustering analysis of DIR gene expression. These analyses further implicated roles for specific DIRs in (-)-pinoresinol formation in seed-coats, as well as (+)-pinoresinol in vegetative organs and/or specific responses to stress. Phylogeny and gene expression analysis segregated flax DIRs into six distinct clusters with new cluster-specific motifs identified. We propose that these findings can serve as a foundation to further systematically determine functions of DIRs, i.e. other than those already known in lignan biosynthesis in flax and other species. Given the differential expression profiles and inducibility of the flax DIR family, we provisionally propose that some DIR genes of unknown function could be involved in different aspects of secondary cell wall biosynthesis and plant defense.

Current Models for Transcriptional Regulation of Secondary Cell Wall Biosynthesis in Grasses

Secondary cell walls mediate many crucial biological processes in plants including mechanical support, water and nutrient transport and stress management. They also provide an abundant resource of renewable feed, fiber, and fuel. The grass family contains the most important food, forage, and biofuel crops. Understanding the regulatory mechanism of secondary wall formation in grasses is necessary for exploiting these plants for agriculture and industry. Previous research has established a detailed model of the secondary wall regulatory network in the dicot model species Arabidopsis thaliana. Grasses, branching off from the dicot ancestor 140-150 million years ago, display distinct cell wall morphology and composition, suggesting potential for a different secondary wall regulation program from that established for dicots. Recently, combined application of molecular, genetic and bioinformatics approaches have revealed more transcription factors involved in secondary cell wall biosynthesis in grasses. Compared with the dicots, grasses exhibit a relatively conserved but nevertheless divergent transcriptional regulatory program to activate their secondary cell wall development and to coordinate secondary wall biosynthesis with other physiological processes.

Kiwifruit R2R3-MYB transcription factors and contribution of the novel AcMYB75 to red kiwifruit anthocyanin biosynthesis

Red kiwifruit (Actinidia chinensis) is a popular fresh fruit with a high market value due to its unique color, caused by anthocyanin accumulation. The R2R3-MYB transcription factors (TFs) have important roles in plant development and anthocyanin metabolism. In this first comprehensive study of R2R3-MYBs in kiwifruit, a total of 93 R2R3-MYB genes, including five novel previously unannotated AcMYBs, were identified. Their phylogenic relationship, exon-intron structures, and conserved motifs were analyzed. Based on transcriptome data, 60 AcMYBs were expressed (FPKM > 1) across seven developmental stages of kiwifruit, revealing five expression patterns. One of the 5 newly identified R2R3 TFs, AcMYB75, showed an anthocyanin accumulation-linked expression pattern during fruit development. AcMYB75 localized to the nucleus and has an active transactivation domain, verifying it as a transcription factor. AcMYB75 protein specifically bound the promoter of the anthocyanin biosynthesis gene ANS in yeast one-hybrid system and in vivo. In 35 S:AcMYB75 Arabidopsis plants, anthocyanin significantly accumulated in leaves, and the expression of anthocyanin biosynthetic genes was greatly up-regulated. Together, these results suggest that AcMYB75 is involved in anthocyanin biosynthesis in kiwifruit. These findings will increase our understanding of AcMYBs involved in anthocyanin biosynthesis, and also benefit further functional characterization of R2R3-MYB genes in kiwifruit.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

MYB-mediated upregulation of lignin biosynthesis in Oryza sativa towards biomass refinery

Lignin encrusts lignocellulose polysaccharides, and has long been considered an obstacle for the efficient use of polysaccharides during processes such as pulping and bioethanol fermentation. However, lignin is also a potential feedstock for aromatic products and is an important by-product of polysaccharide utilization. Therefore, producing biomass plant species exhibiting enhanced lignin production is an important breeding objective. Herein, we describe the development of transgenic rice plants with increased lignin content. Five Arabidopsis thaliana (Arabidopsis) and one Oryza sativa (rice) MYB transcription factor genes that were implicated to be involved in lignin biosynthesis were transformed into rice (O. sativa L. ssp. japonica cv. Nipponbare). Among them, three Arabidopsis MYBs (AtMYB55, AtMYB61, and AtMYB63) in transgenic rice T₁ lines resulted in culms with lignin content about 1.5-fold higher than that of control plants. Furthermore, lignin structures in AtMYB61-overexpressing rice plants were investigated by wet-chemistry and two-dimensional nuclear magnetic resonance spectroscopy approaches. Our data suggested that heterologous expression of AtMYB61 in rice increased lignin content mainly by enriching syringyl units as well as p-coumarate and tricin moieties in the lignin polymers. We contemplate that this strategy is also applicable to lignin upregulation in large-sized grass biomass plants, such as Sorghum, switchgrass, Miscanthus and Erianthus.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

Co-expression network analysis reveals transcription factors associated to cell wall biosynthesis in sugarcane

Sugarcane is a hybrid of Saccharum officinarum and Saccharum spontaneum, with minor contributions from other species in Saccharum and other genera. Understanding the molecular basis of cell wall metabolism in sugarcane may allow for rational changes in fiber quality and content when designing new energy crops. This work describes a comparative expression profiling of sugarcane ancestral genotypes: S. officinarum, S. spontaneum and S. robustum and a commercial hybrid: RB867515, linking gene expression to phenotypes to identify genes for sugarcane improvement. Oligoarray experiments of leaves, immature and intermediate internodes, detected 12,621 sense and 995 antisense transcripts. Amino acid metabolism was particularly evident among pathways showing natural antisense transcripts expression. For all tissues sampled, expression analysis revealed 831, 674 and 648 differentially expressed genes in S. officinarum, S. robustum and S. spontaneum, respectively, using RB867515 as reference. Expression of sugar transporters might explain sucrose differences among genotypes, but an unexpected differential expression of histones were also identified between high and low Brix° genotypes. Lignin biosynthetic genes and bioenergetics-related genes were up-regulated in the high lignin genotype, suggesting that these genes are important for S. spontaneum to allocate carbon to lignin, while S. officinarum allocates it to sucrose storage. Co-expression network analysis identified 18 transcription factors possibly related to cell wall biosynthesis while in silico analysis detected cis-elements involved in cell wall biosynthesis in their promoters. Our results provide information to elucidate regulatory networks underlying traits of interest that will allow the improvement of sugarcane for biofuel and chemicals production.

Genome-Wide Identification of R2R3-MYB Genes and Expression Analyses During Abiotic Stress in Gossypium raimondii

The R2R3-MYB is one of the largest families of transcription factors, which have been implicated in multiple biological processes. There is great diversity in the number of R2R3-MYB genes in different plants. However, there is no report on genome-wide characterization of this gene family in cotton. In the present study, a total of 205 putative R2R3-MYB genes were identified in cotton D genome (Gossypium raimondii), that are much larger than that found in other cash crops with fully sequenced genomes. These GrMYBs were classified into 13 groups with the R2R3-MYB genes from Arabidopsis and rice. The amino acid motifs and phylogenetic tree were predicted and analyzed. The sequences of GrMYBs were distributed across 13 chromosomes at various densities. The results showed that the expansion of the G. Raimondii R2R3-MYB family was mainly attributable to whole genome duplication and segmental duplication. Moreover, the expression pattern of 52 selected GrMYBs and 46 GaMYBs were tested in roots and leaves under different abiotic stress conditions. The results revealed that the MYB genes in cotton were differentially expressed under salt and drought stress treatment. Our results will be useful for determining the precise role of the MYB genes during stress responses with crop improvement.

A Novel R2R3-MYB Transcription Factor BpMYB106 of Birch (Betula platyphylla) Confers Increased Photosynthesis and Growth Rate through Up-regulating Photosynthetic Gene Expression

We isolated a R2R3-MYB transcription factor BpMYB106, which regulates photosynthesis in birch (Betula platyphylla Suk.). BpMYB106 mainly expresses in the leaf and shoot tip of birch, and its protein is localized in the nucleus. We further fused isolated a 1588 bp promoter of BpMYB106 and analyzed it by PLACE, which showed some cis-acting elements related to photosynthesis. BpMYB106 promoter β-glucuronidase (GUS) reporter fusion studies gene, the result, showed the GUS reporter gene in transgenic birch with BpMYB106 promoter showed strong activities in shoot tip, cotyledon margins, and mature leaf trichomes. The overexpression of BpMYB106 in birch resulted in significantly increased trichome density, net photosynthetic rate, and growth rate as compared with the wild-type birch. RNA-Seq profiling revealed the upregulation of several photosynthesis-related genes in the photosynthesis and oxidative phosphorylation pathways in the leaves of transgenic plants. Yeast one-hybrid analysis, coupled with transient assay in tobacco, revealed that BpMYB106 binds a MYB binding site MYB2 in differentially expressed gene promoters. Thus, BpMYB106 may directly activate the expression of a range of photosynthesis related genes through interacting with the MYB2 element in their promoters. Our study demonstrating the overexpression of BpMYB106-a R2R3-MYB transcription factor-upregulates the genes of the photosynthesis and oxidative phosphorylation pathways to improve photosynthesis.

Cotton GhMYB7 is predominantly expressed in developing fibers and regulates secondary cell wall biosynthesis in transgenic Arabidopsis

The secondary cell wall in mature cotton fibers contains over 90% cellulose with low quantities of xylan and lignin. However, little is known regarding the regulation of secondary cell wall biosynthesis in cotton fibers. In this study, we characterized an R2R3-MYB transcription factor, GhMYB7, in cotton. GhMYB7 is expressed at a high level in developing fibers and encodes a MYB protein that is targeted to the cell nucleus and has transcriptional activation activity. Ectopic expression of GhMYB7 in Arabidopsis resulted in small, curled, dark green leaves and also led to shorter inflorescence stems. A cross-sectional assay of basal stems revealed that cell wall thickness of vessels and interfascicular fibers was higher in transgenic lines overexpressing GhMYB7 than in the wild type. Constitutive expression of GhMYB7 in Arabidopsis activated the expression of a suite of secondary cell wall biosynthesis-related genes (including some secondary cell wall-associated transcription factors), leading to the ectopic deposition of cellulose and lignin. The ectopic deposition of secondary cell walls may have been initiated before the cessation of cell expansion. Moreover, GhMYB7 was capable of binding to the promoter regions of AtSND1 and AtCesA4, suggesting that GhMYB7 may function upstream of NAC transcription factors. Collectively, these findings suggest that GhMYB7 is a potential transcriptional activator, which may participate in regulating secondary cell wall biosynthesis of cotton fibers.

Breeding maize for silage and biofuel production, an illustration of a step forward with the genome sequence

The knowledge of the gene families mostly impacting cell wall digestibility variations would significantly increase the efficiency of marker-assisted selection when breeding maize and grass varieties with improved silage feeding value and/or with better straw fermentability into alcohol or methane. The maize genome sequence of the B73 inbred line was released at the end of 2009, opening up new avenues to identify the genetic determinants of quantitative traits. Colocalizations between a large set of candidate genes putatively involved in secondary cell wall assembly and QTLs for cell wall digestibility (IVNDFD) were then investigated, considering physical positions of both genes and QTLs. Based on available data from six RIL progenies, 59 QTLs corresponding to 38 non-overlapping positions were matched up with a list of 442 genes distributed all over the genome. Altogether, 176 genes colocalized with IVNDFD QTLs and most often, several candidate genes colocalized at each QTL position. Frequent QTL colocalizations were found firstly with genes encoding ZmMYB and ZmNAC transcription factors, and secondly with genes encoding zinc finger, bHLH, and xylogen regulation factors. In contrast, close colocalizations were less frequent with genes involved in monolignol biosynthesis, and found only with the C4H2, CCoAOMT5, and CCR1 genes. Close colocalizations were also infrequent with genes involved in cell wall feruloylation and cross-linkages. Altogether, investigated colocalizations between candidate genes and cell wall digestibility QTLs suggested a prevalent role of regulation factors over constitutive cell wall genes on digestibility variations.

Exogenous GA₃ Application Enhances Xylem Development and Induces the Expression of Secondary Wall Biosynthesis Related Genes in Betula platyphylla

Gibberellin (GA) is a key signal molecule inducing differentiation of tracheary elements, fibers, and xylogenesis. However the molecular mechanisms underlying the effect of GA on xylem elongation and secondary wall development in tree species remain to be determined. In this study, Betula platyphylla (birch) seeds were treated with 300 ppm GA₃ and/or 300 ppm paclobutrazol (PAC), seed germination was recorded, and transverse sections of hypocotyls were stained with toluidine blue; the two-month-old seedlings were treated with 50 μM GA₃ and/or 50 μM PAC, transverse sections of seedling stems were stained using phloroglucinol–HCl, and secondary wall biosynthesis related genes expression was analyzed by real-time quantitative PCR. Results indicated that germination percentage, energy and time of seeds, hypocotyl height and seedling fresh weight were enhanced by GA₃, and reduced by PAC; the xylem development was wider in GA₃-treated plants than in the control; the expression of NAC and MYB transcription factors, CESA, PAL, and GA oxidase was up-regulated during GA₃ treatment, suggesting their role in GA₃-induced xylem development in the birch. Our results suggest that GA₃ induces the expression of secondary wall biosynthesis related genes to trigger xylogenesis in the birch plants.

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.

Functional Characterization of NAC and MYB Transcription Factors Involved in Regulation of Biomass Production in Switchgrass (Panicum virgatum)

Switchgrass is a promising biofuel feedstock due to its high biomass production and low agronomic input requirements. Because the bulk of switchgrass biomass used for biofuel production is lignocellulosic secondary walls, studies on secondary wall biosynthesis and its transcriptional regulation are imperative for designing strategies for genetic improvement of biomass production in switchgrass. Here, we report the identification and functional characterization of a group of switchgrass transcription factors, including several NACs (PvSWNs) and a MYB (PvMYB46A), for their involvement in regulating secondary wall biosynthesis. PvSWNs and PvMYB46A were found to be highly expressed in stems and their expression was closely associated with sclerenchyma cells. Overexpression of PvSWNs and PvMYB46A in Arabidopsis was shown to result in activation of the biosynthetic genes for cellulose, xylan and lignin and ectopic deposition of secondary walls in normally parenchymatous cells. Transactivation and complementation studies demonstrated that PvSWNs were able to activate the SNBE-driven GUS reporter gene and effectively rescue the secondary wall defects in the Arabidopsis snd1 nst1 double mutant, indicating that they are functional orthologs of Arabidopsis SWNs. Furthermore, we showed that PvMYB46A could activate the SMRE-driven GUS reporter gene and complement the Arabidopsis myb46 myb83 double mutant, suggesting that it is a functional ortholog of Arabidopsis MYB46/MYB83. Together, these results indicate that PvSWNs and PvMYB46A are transcriptional switches involved in regulating secondary wall biosynthesis, which provides molecular tools for genetic manipulation of biomass production in switchgrass.

The Eucalyptus grandis R2R3-MYB transcription factor family: evidence for woody growth-related evolution and function

The R2R3-MYB family, one of the largest transcription factor families in higher plants, controls a wide variety of plant-specific processes including, notably, phenylpropanoid metabolism and secondary cell wall formation. We performed a genome-wide analysis of this superfamily in Eucalyptus, one of the most planted hardwood trees world-wide. A total of 141 predicted R2R3-MYB sequences identified in the Eucalyptus grandis genome sequence were subjected to comparative phylogenetic analyses with Arabidopsis thaliana, Oryza sativa, Populus trichocarpa and Vitis vinifera. We analysed features such as gene structure, conserved motifs and genome location. Transcript abundance patterns were assessed by RNAseq and validated by high-throughput quantitative PCR. We found some R2R3-MYB subgroups with expanded membership in E. grandis, V. vinifera and P. trichocarpa, and others preferentially found in woody species, suggesting diversification of specific functions in woody plants. By contrast, subgroups containing key genes regulating lignin biosynthesis and secondary cell wall formation are more conserved across all of the species analysed. In Eucalyptus, R2R3-MYB tandem gene duplications seem to disproportionately affect woody-preferential and woody-expanded subgroups. Interestingly, some of the genes belonging to woody-preferential subgroups show higher expression in the cambial region, suggesting a putative role in the regulation of secondary growth.

NAC-MYB-based transcriptional regulation of secondary cell wall biosynthesis in land plants

Plant cells biosynthesize primary cell walls (PCW) in all cells and produce secondary cell walls (SCWs) in specific cell types that conduct water and/or provide mechanical support, such as xylem vessels and fibers. The characteristic mechanical stiffness, chemical recalcitrance, and hydrophobic nature of SCWs result from the organization of SCW-specific biopolymers, i.e., highly ordered cellulose, hemicellulose, and lignin. Synthesis of these SCW-specific biopolymers requires SCW-specific enzymes that are regulated by SCW-specific transcription factors. In this review, we summarize our current knowledge of the transcriptional regulation of SCW formation in plant cells. Advances in research on SCW biosynthesis during the past decade have expanded our understanding of the transcriptional regulation of SCW formation, particularly the functions of the NAC and MYB transcription factors. Focusing on the NAC-MYB-based transcriptional network, we discuss the regulatory systems that evolved in land plants to modify the cell wall to serve as a key component of structures that conduct water and provide mechanical support.

WD40-Repeat Proteins in Plant Cell Wall Formation: Current Evidence and Research Prospects

The metabolic complexity of living organisms relies on supramolecular protein structures which ensure vital processes, such as signal transduction, transcription, translation and cell wall synthesis. In eukaryotes WD40-repeat (WDR) proteins often function as molecular "hubs" mediating supramolecular interactions. WDR proteins may display a variety of interacting partners and participate in the assembly of complexes involved in distinct cellular functions. In plants, the formation of lignocellulosic biomass involves extensive synthesis of cell wall polysaccharides, a process that requires the assembly of large transmembrane enzyme complexes, intensive vesicle trafficking, interactions with the cytoskeleton, and coordinated gene expression. Because of their function as supramolecular hubs, WDR proteins could participate in each or any of these steps, although to date only few WDR proteins have been linked to the cell wall by experimental evidence. Nevertheless, several potential cell wall-related WDR proteins were recently identified using in silico approaches, such as analyses of co-expression, interactome and conserved gene neighborhood. Notably, some WDR genes are frequently genomic neighbors of genes coding for GT2-family polysaccharide synthases in eukaryotes, and this WDR-GT2 collinear microsynteny is detected in diverse taxa. In angiosperms, two WDR genes are collinear to cellulose synthase genes, CesAs, whereas in ascomycetous fungi several WDR genes are adjacent to chitin synthase genes, chs. In this Perspective we summarize and discuss experimental and in silico studies on the possible involvement of WDR proteins in plant cell wall formation. The prospects of biotechnological engineering for enhanced biomass production are discussed.

Evaluation and integration of functional annotation pipelines for newly sequenced organisms: the potato genome as a test case

BACKGROUND

For most organisms, even if their genome sequence is available, little functional information about individual genes or proteins exists. Several annotation pipelines have been developed for functional analysis based on sequence, 'omics', and literature data. However, researchers encounter little guidance on how well they perform. Here, we used the recently sequenced potato genome as a case study. The potato genome was selected since its genome is newly sequenced and it is a non-model plant even if there is relatively ample information on individual potato genes, and multiple gene expression profiles are available.

RESULTS

We show that the automatic gene annotations of potato have low accuracy when compared to a "gold standard" based on experimentally validated potato genes. Furthermore, we evaluate six state-of-the-art annotation pipelines and show that their predictions are markedly dissimilar (Jaccard similarity coefficient of 0.27 between pipelines on average). To overcome this discrepancy, we introduce a simple GO structure-based algorithm that reconciles the predictions of the different pipelines. We show that the integrated annotation covers more genes, increases by over 50% the number of highly co-expressed GO processes, and obtains much higher agreement with the gold standard.

CONCLUSIONS

We find that different annotation pipelines produce different results, and show how to integrate them into a unified annotation that is of higher quality than each single pipeline. We offer an improved functional annotation of both PGSC and ITAG potato gene models, as well as tools that can be applied to additional pipelines and improve annotation in other organisms. This will greatly aid future functional analysis of '-omics' datasets from potato and other organisms with newly sequenced genomes. The new potato annotations are available with this paper.

Contrasting metabolism in perenniating structures of upland and lowland switchgrass plants late in the growing season

Background

Switchgrass (Panicum virgatum L.) is being developed as a bioenergy crop for many temperate regions of the world. One way to increase biomass yields is to move southern adapted lowland cultivars to more northern latitudes. However, many southerly adapted switchgrass germplasm can suffer significant winter kill in northerly climes.

Materials and Methods

Here, we have applied next-generation sequencing in combination with biochemical analyses to query the metabolism of crowns and rhizomes obtained from two contrasting switchgrass cultivars. Crowns and rhizomes from field-grown lowland (cv Kanlow) and upland (cv Summer) switchgrass cultivars were collected from three randomly selected post-flowering plants. Summer plants were senescing, whereas Kanlow plants were not at this harvest date.

Results

Principal component analysis (PCA) differentiated between both the Summer and Kanlow transcriptomes and metabolomes. Significant differences in transcript abundances were detected for 8,050 genes, including transcription factors such as WRKYs and those associated with phenylpropanoid biosynthesis. Gene-set enrichment analyses showed that a number of pathways were differentially up-regulated in the two populations. For both populations, protein levels and enzyme activities agreed well with transcript abundances for genes involved in the phenylpropanoid pathway that were up-regulated in Kanlow crowns and rhizomes. The combination of these datasets suggests that dormancy-related mechanisms had been triggered in the crowns and rhizomes of the Summer plants, whereas the crowns and rhizomes of Kanlow plants had yet to enter dormancy.

Conclusions

Delayed establishment of dormancy at more northerly latitudes could be one factor that reduces winter-survival in the high-yielding Kanlow plants. Understanding the cellular signatures that accompany the transition to dormancy can be used in the future to select plants with improved winter hardiness.