Navigating the transcriptional roadmap regulating plant secondary cell wall deposition

Identification of a novel viral factor inducing tumorous symptoms by disturbing vascular development in planta

Plant viruses induce various disease symptoms that substantially impact agriculture, but the underlying mechanisms of viral disease in plants are poorly understood. Kobu-sho is a disease in gentian that shows gall formation with ectopic development of lignified cells and vascular tissues such as xylem. Here, we show that a gene fragment of gentian Kobu-sho-associated virus, which is designated as Kobu-sho-inducing factor (KOBU), induces gall formation accompanied by ectopic development of lignified cells and xylem-like tissue in Nicotiana benthamiana. Transgenic gentian expressing KOBU exhibited tumorous symptoms, confirming the gall-forming activity of KOBU. Surprisingly, KOBU expression can also induce differentiation of an additional leaf-like tissue on the abaxial side of veins in normal N. benthamiana and gentian leaves. Transcriptome analysis with Arabidopsis thaliana expressing KOBU revealed that KOBU activates signaling pathways that regulate xylem development. KOBU protein forms granules and plate-like structures and co-localizes with mRNA splicing factors within the nucleus. Our findings suggest that KOBU is a novel pleiotropic virulence factor that stimulates vascular and leaf development. IMPORTANCE While various mechanisms determine disease symptoms in plants depending on virus-host combinations, the details of how plant viruses induce symptoms remain largely unknown in most plant species. Kobu-sho is a disease in gentian that shows gall formation with ectopic development of lignified cells and vascular tissues such as xylem. Our findings demonstrate that a gene fragment of gentian Kobu-sho-associated virus (GKaV), which is designated as Kobu-sho-inducing factor, induces the gall formation accompanied by the ectopic development of lignified cells and xylem-like tissue in Nicotiana benthamiana. The molecular mechanism by which gentian Kobu-sho-associated virus induces the Kobu-sho symptoms will provide new insight into not only plant-virus interactions but also the regulatory mechanisms underlying vascular and leaf development.

The ARF2-MYB6 module mediates auxin-regulated petal expansion in rose

In cut rose (Rosa hybrida), the flower-opening process is closely associated with vase life. Auxin induces the expression of transcription factor genes that function in petal growth via cell expansion. However, the molecular mechanisms underlying the auxin effect during flower opening are not well understood. Here, we identified the auxin-inducible transcription factor gene RhMYB6, whose expression level is high during the early stages of flower opening. Silencing of RhMYB6 delayed flower opening by controlling petal cell expansion through down-regulation of cell expansion-related genes. Furthermore, we demonstrated that the auxin response factor RhARF2 directly interacts with the promoter of RhMYB6 and represses its transcription. Silencing of RhARF2 resulted in larger petal size and delayed petal movement. We also showed that the expression of genes related to ethylene and petal movement showed substantial differences in RhARF2-silenced petals. Our results indicate that auxin-regulated RhARF2 is a critical player that controls flower opening by governing RhMYB6 expression and mediating the crosstalk between auxin and ethylene signaling.

The transcription factor PbrMYB24 regulates lignin and cellulose biosynthesis in stone cells of pear fruits

Lignified stone cell content is a key factor used to evaluate fruit quality, influencing the economic value of pear (Pyrus pyrifolia) fruits. However, our understanding of the regulatory networks of stone cell formation is limited due to the complex secondary metabolic pathway. In this study, we used a combination of co-expression network analysis, gene expression profiles, and transcriptome analysis in different pear cultivars with varied stone cell content to identify a hub MYB gene, PbrMYB24. The relative expression of PbrMYB24 in fruit flesh was significantly correlated with the contents of stone cells, lignin, and cellulose. We then verified the function of PbrMYB24 in regulating lignin and cellulose formation via genetic transformation in homologous and heterologous systems. We constructed a high-efficiency verification system for lignin and cellulose biosynthesis genes in pear callus. PbrMYB24 transcriptionally activated multiple target genes involved in stone cell formation. On the one hand, PbrMYB24 activated the transcription of lignin and cellulose biosynthesis genes by binding to different cis-elements [AC-I (ACCTACC) element, AC-II (ACCAACC) element and MYB-binding sites (MBS)]. On the other hand, PbrMYB24 bound directly to the promoters of PbrMYB169 and NAC STONE CELL PROMOTING FACTOR (PbrNSC), activating the gene expression. Moreover, both PbrMYB169 and PbrNSC activated the promoter of PbrMYB24, enhancing gene expression. This study improves our understanding of lignin and cellulose synthesis regulation in pear fruits through identifying a regulator and establishing a regulatory network. This knowledge will be useful for reducing the stone cell content in pears via molecular breeding.

Comparative functional analysis of PdeNAC2 and AtVND6 in the tracheary element formation

Tracheary elements (i.e., vessel elements and tracheids) are highly specialized, non-living cells present in the water-conducting xylem tissue. In angiosperms, proteins in the VASCULAR-RELATED NAC-DOMAIN (VND) subgroup of the NAC transcription factor family (e.g., AtVND6) are required for the differentiation of vessel elements through transcriptional regulation of genes responsible for secondary cell wall (SCW) formation and programmed cell death (PCD). Gymnosperms, however, produce only tracheids, the mechanism of which remains elusive. Here, we report functional characteristics of PdeNAC2, a VND homolog in Pinus densiflora, as a key regulator of tracheid formation. Interestingly, our molecular genetic analyses show that PdeNAC2 can induce the formation of vessel element-like cells in angiosperm plants, demonstrated by transgenic overexpression of either native or NAC domain-swapped synthetic genes of PdeNAC2 and AtVND6 in both Arabidopsis and hybrid poplar. Subsequently, genome-wide identification of direct target genes of PdeNAC2 and AtVND6 revealed 138 and 174 genes as putative direct targets, respectively, but only 17 genes were identified as common direct targets. Further analyses have found that PdeNAC2 does not control some AtVND6-dependent vessel differentiation genes in angiosperm plants, such as AtVRLK1, LBD15/30, and pit-forming ROP signaling genes. Collectively, our results suggest that different target gene repertoires of PdeNAC2 and AtVND6 may contribute to the evolution of tracheary elements.

PtrWOX13A Promotes Wood Formation and Bioactive Gibberellins Biosynthesis in Populus trichocarpa

WUSCHEL-related homeobox (WOX) genes are plant-specific transcription factors (TFs) involved in multiple processes of plant development. However, there have hitherto no studies on the WOX TFs involved in secondary cell wall (SCW) formation been reported. In this study, we identified a Populus trichocarpa WOX gene, PtrWOX13A, which was predominantly expressed in SCW, and then characterized its functions through generating PtrWOX13A overexpression poplar transgenic lines; these lines exhibited not only significantly enhanced growth potential, but also remarkably increased SCW thicknesses, fiber lengths, and lignin and hemicellulose contents. However, no obvious change in cellulose content was observed. We revealed that PtrWOX13A directly activated its target genes through binding to two cis-elements, ATTGATTG and TTAATSS, in their promoter regions. The fact that PtrWOX13A responded to the exogenous GAs implies that it is responsive to GA homeostasis caused by GA inactivation and activation genes (e.g., PtrGA20ox4, PtrGA2ox1, and PtrGA3ox1), which were regulated by PtrWOX13A directly or indirectly. Since the master switch gene of SCW formation, PtrWND6A, and lignin biosynthesis regulator, MYB28, significantly increased in PtrWOX13A transgenic lines, we proposed that PtrWOX13A, as a higher hierarchy TF, participated in SCW formation through controlling the genes that are components of the known hierarchical transcription regulation network of poplar SCW formation, and simultaneously triggering a gibberellin-mediated signaling cascade. The discovery of PtrWOX13A predominantly expressed in SCW and its regulatory functions in the poplar wood formation has important implications for improving the wood quality of trees via genetic engineering.

PtrHAT22, as a higher hierarchy regulator, coordinately regulates secondary cell wall component biosynthesis in Populus trichocarpa

Homeodomain-leucine zipper (HD-Zip) II transcription factors (TFs) have been reported to play vital roles in diverse biological processes of plants. However, it remains unclear whether HD-Zip II TFs regulate secondary cell wall (SCW) in woody plants. In this study, we performed the functional characterization of a Populus trichocarpa HD-Zip II TF, PtrHAT22, which encodes a nuclear localized transcription repressor predominantly expressing in secondary developing tissues. Overexpression of PtrHAT22 showed arrested growths, including reduced heights and diameters above the ground, small leaves, and decreased biomass. Meanwhile, the contents of lignin, cellulose, and thickness of SCW significantly decreased, whilst the content of hemicellulose obviously increased in PtrHAT22 transgenic poplar. The expressions of some wood-associated TFs and structural genes significantly changed accordingly with the alternations of SCW characteristics in PtrHAT22 transgenic poplar. Furthermore, PtrHAT22 directly repressed the promoter activities of PtrMYB20, PtrMYB28, and PtrCOMT2, and bind two cis-acting elements that were specifically enriched in their promoter regions. Taken together, our results suggested that PtrHAT22, as a higher hierarchy TF like PtrWNDs, exerted coordination regulation of poplar SCW component biosynthesis through directly and indirectly regulating structural genes and different hierarchy TFs of SCW formation network.

A vacuolar hexose transport is required for xylem development in the inflorescence stem

In Angiosperms, the development of the vascular system is controlled by a complex network of transcription factors. However, how nutrient availability in the vascular cells affects their development remains to be addressed. At the cellular level, cytosolic sugar availability is regulated mainly by sugar exchanges at the tonoplast through active and/or facilitated transport. In Arabidopsis (Arabidopsis thaliana), among the genes encoding tonoplastic transporters, SUGAR WILL EVENTUALLY BE EXPORTED TRANSPORTER 16 (SWEET16) and SWEET17 expression has been previously detected in the vascular system. Here, using a reverse genetics approach, we propose that sugar exchanges at the tonoplast, regulated by SWEET16, are important for xylem cell division as revealed in particular by the decreased number of xylem cells in the swt16 mutant and the accumulation of SWEET16 at the procambium-xylem boundary. In addition, we demonstrate that transport of hexoses mediated by SWEET16 and/or SWEET17 is required to sustain the formation of the xylem secondary cell wall. This result is in line with a defect in the xylem cell wall composition as measured by Fourier-transformed infrared spectroscopy in the swt16swt17 double mutant and by upregulation of several genes involved in secondary cell wall synthesis. Our work therefore supports a model in which xylem development partially depends on the exchange of hexoses at the tonoplast of xylem-forming cells.

Histone Deacetylation Controls Xylem Vessel Cell Differentiation via Transcriptional Regulation of a Transcription Repressor Complex OFP1/4-MYB75-KNAT7-BLH6

Xylem vessels are indispensable tissues in vascular plants that transport water and minerals. The differentiation of xylem vessel cells is characterized by secondary cell wall deposition and programmed cell death. These processes are initiated by a specific set of transcription factors, called VASCULAR-RELATED NAC-DOMAIN (VND) family proteins, through the direct and/or indirectly induction of genes required for secondary cell wall deposition and programmed cell death. In this study, we explored novel regulatory factors for xylem vessel cell differentiation in Arabidopsis thaliana. We tested the effects of cellular stress inducers on VND7-induced differentiation of xylem vessel cells with the VND7-VP16-GR system, in which VND7 activity is post-translationally induced by dexamethasone application. We established that the histone deacetylase (HDAC) inhibitors trichostatin A (TSA) and sirtinol inhibited VND7-induced xylem vessel cell differentiation. The inhibitory effects of TSA and sirtinol treatment were detected only when they were added at the same time as the dexamethasone application, suggesting that TSA and sirtinol mainly influence the early stages of xylem vessel cell differentiation. Expression analysis revealed that these HDAC inhibitors downregulated VND7-downstream genes, including both direct and indirect targets of transcriptional activation. Notably, the HDAC inhibitors upregulated the transcript levels of negative regulators of xylem vessel cells, OVATE FAMILY PROTEIN1 (OFP1), OFP4, and MYB75, which are known to form a protein complex with BEL1-LIKE HOMEODOMAIN6 (BLH6) to repress gene transcription. The KDB system, another in vitro induction system of ectopic xylem vessel cells, demonstrated that TSA and sirtinol also inhibited ectopic formation of xylem vessel cells, and this inhibition was partially suppressed in knat7-1, bhl6-1, knat7-1 bhl6-1, and quintuple ofp1 ofp2 ofp3 ofp4 ofp5 mutants. Thus, the negative effects of HDAC inhibitors on xylem vessel cell differentiation are mediated, at least partly, by the abnormal upregulation of the transcriptional repressor complex OFP1/4-MYB75-KNAT7-BLH6. Collectively, our findings suggest that active regulation of histone deacetylation by HDACs is involved in xylem vessel cell differentiation via the OFP1/4-MYB75-KNAT7-BLH6 complex.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Xylem vessel cell differentiation: A best model for new integrative cell biology?

Xylem vessels transport water and essential low-molecular-weight compounds throughout vascular plants. To achieve maximum performance as conductive tissues, xylem vessel cells undergo secondary cell wall deposition and programmed cell death to produce a hollow tube-like structure with a rigid outer shell. This unique process has been explored in detail from a cell biology and molecular biology perspective, culminating in the identification of the master transcriptional switches of xylem vessel cell differentiation, the VASCULAR-RELATED NAC-DOMAIN (VND) proteins. High-resolution analyses of xylem vessel cell differentiation have since accelerated and are now moving toward single cell-level dissection from a variety of directions. In this review, we introduce the current model of xylem vessel cell differentiation and discuss possible future directions in this field.

Multi-omics network-based functional annotation of unknown Arabidopsis genes

Unraveling gene function is pivotal to understanding the signaling cascades that control plant development and stress responses. As experimental profiling is costly and labor intensive, there is a clear need for high-confidence computational annotation. In contrast to detailed gene-specific functional information, transcriptomics data are widely available for both model and crop species. Here, we describe a novel automated function prediction method, which leverages complementary information from multiple expression datasets by analyzing study-specific gene co-expression networks. First, we benchmarked the prediction performance on recently characterized Arabidopsis thaliana genes, and showed that our method outperforms state-of-the-art expression-based approaches. Next, we predicted biological process annotations for known (n = 15 790) and unknown (n = 11 865) genes in A. thaliana and validated our predictions using experimental protein-DNA and protein-protein interaction data (covering >220 000 interactions in total), obtaining a set of high-confidence functional annotations. Our method assigned at least one validated annotation to 5054 (42.6%) unknown genes, and at least one novel validated function to 3408 (53.0%) genes with computational annotations only. These omics-supported functional annotations shed light on a variety of developmental processes and molecular responses, such as flower and root development, defense responses to fungi and bacteria, and phytohormone signaling, and help fill the information gap on biological process annotations in Arabidopsis. An in-depth analysis of two context-specific networks, modeling seed development and response to water deprivation, shows how previously uncharacterized genes function within the respective networks. Moreover, our automated function prediction approach can be applied in future studies to facilitate gene discovery for crop improvement.

A regulatory network driving shoot lignification in rapidly growing bamboo

Woody bamboo is environmentally friendly, abundant, and an alternative to conventional timber. Degree of lignification and lignin content and deposition affect timber properties. However, the lignification regulatory network in monocots is poorly understood. To elucidate the regulatory mechanism of lignification in moso bamboo (Phyllostachys edulis), we conducted integrated analyses using transcriptome, small RNA, and degradome sequencing followed by experimental verification. The lignification degree and lignin content increased with increased bamboo shoot height, whereas phenylalanine ammonia-lyase and Laccase activities first increased and then decreased with shoot growth. Moreover, we identified 11,504 differentially expressed genes (DEGs) in different portions of the 13th internodes of different height shoots; most DEGs associated with cell wall and lignin biosynthesis were upregulated, whereas some DEGs related to cell growth were downregulated. We identified a total of 1,502 miRNAs, of which 687 were differentially expressed. Additionally, in silico and degradome analyses indicated that 5,756 genes were targeted by 691 miRNAs. We constructed a regulatory network of lignification, including 11 miRNAs, 22 transcription factors, and 36 enzyme genes, in moso bamboo. Furthermore, PeLAC20 overexpression increased lignin content in transgenic Arabidopsis (Arabidopsis thaliana) plants. Finally, we proposed a reliable miRNA-mediated "MYB-PeLAC20" module for lignin monomer polymerization. Our findings provide definite insights into the genetic regulation of bamboo lignification. In addition to providing a platform for understanding related mechanisms in other monocots, these insights could be used to develop strategies to improve bamboo timber properties.

Wood Formation under Severe Drought Invokes Adjustment of the Hormonal and Transcriptional Landscape in Poplar

Drought is a severe environmental stress that exerts negative effects on plant growth. In trees, drought leads to reduced secondary growth and altered wood anatomy. The mechanisms underlying wood stress adaptation are not well understood. Here, we investigated the physiological, anatomical, hormonal, and transcriptional responses of poplar to strong drought. Drought-stressed xylem was characterized by higher vessel frequencies, smaller vessel lumina, and thicker secondary fiber cell walls. These changes were accompanied by strong increases in abscisic acid (ABA) and antagonistic changes in salicylic acid in wood. Transcriptional evidence supported ABA biosynthesis and signaling in wood. Since ABA signaling activates the fiber-thickening factor NST1, we expected upregulation of the secondary cell wall (SCW) cascade under stress. By contrast, transcription factors and biosynthesis genes for SCW formation were down-regulated, whereas a small set of cellulose synthase-like genes and a huge array of genes involved in cell wall modification were up-regulated in drought-stressed wood. Therefore, we suggest that ABA signaling monitors normal SCW biosynthesis and that drought causes a switch from normal to "stress wood" formation recruiting a dedicated set of genes for cell wall biosynthesis and remodeling. This proposition implies that drought-induced changes in cell wall properties underlie regulatory mechanisms distinct from those of normal wood.

Posttranscriptional regulation of cellulose synthase genes by small RNAs derived from cellulose synthase antisense transcripts

Transcriptional regulatory mechanisms governing plant cell wall biosynthesis are incomplete. Expression programs that activate wall biosynthesis are well understood, but mechanisms that control the attenuation of gene expression networks remain elusive. Previous work has shown that small RNAs (sRNAs) derived from the HvCESA6 (Hordeum vulgare, Hv) antisense transcripts are naturally produced and are capable of regulating aspects of wall biosynthesis. Here, we further test the hypothesis that CESA-derived sRNAs generated from CESA antisense transcripts are involved in the regulation of cellulose and broader cell wall biosynthesis. Antisense transcripts were detected for some but not all members of the CESA gene family in both barley and Brachypodium distachyon. Phylogenetic analysis indicates that antisense transcripts are detected for most primary cell wall CESA genes, suggesting a possible role in the transition from primary to secondary cell wall biosynthesis. Focusing on one antisense transcript, HvCESA1 shows dynamic expression throughout development, is correlated with corresponding sRNAs over the same period and is anticorrelated with HvCESA1 mRNA expression. To assess the broader impacts of CESA-derived sRNAs on the regulation of cell wall biosynthesis, transcript profiling was performed on barley tissues overexpressing CESA-derived sRNAs. Together, the data support the hypothesis that CESA antisense transcripts function through an RNA-induced silencing mechanism, to degrade cis transcripts, and may also trigger trans-acting silencing on related genes to alter the expression of cell wall gene networks.

Wood transcriptome analysis of Pinus densiflora identifies genes critical for secondary cell wall formation and NAC transcription factors involved in tracheid formation

Although conifers have significant ecological and economic value, information on transcriptional regulation of wood formation in conifers is still limited. Here, to gain insight into secondary cell wall (SCW) biosynthesis and tracheid formation in conifers, we performed wood tissue-specific transcriptome analyses of Pinus densiflora (Korean red pine) using RNA sequencing. In addition, to obtain full-length transcriptome information, PacBio single molecule real-time iso-sequencing was carried out using RNAs from 28 tissues of P. densiflora. Subsequent comparative tissue-specific transcriptome analysis successfully pinpointed critical genes encoding key proteins involved in biosynthesis of the major secondary wall components (cellulose, galactoglucomannan, xylan and lignin). Furthermore, we predicted a total of 62 NAC (NAM, ATAF1/2 and CUC2) family transcription factor members and identified seven PdeNAC genes preferentially expressed in developing xylem tissues in P. densiflora. Protoplast-based transcriptional activation analysis found that four PdeNAC genes, homologous to VND, NST and SND/ANAC075, upregulated GUS activity driven by an SCW-specific cellulose synthase promoter. Consistently, transient overexpression of the four PdeNACs induced xylem vessel cell-like SCW deposition in both tobacco (Nicotiana benthamiana) and Arabidopsis leaves. Taken together, our data provide a foundation for further research to unravel transcriptional regulation of wood formation in conifers, especially SCW formation and tracheid differentiation.

De novo transcriptome analysis of white teak (Gmelina arborea Roxb) wood reveals critical genes involved in xylem development and secondary metabolism

BACKGROUND

Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees.

RESULTS

For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR.

CONCLUSIONS

We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

Targeted mutation of transcription factor genes alters metaxylem vessel size and number in rice roots

Root metaxylem vessels are responsible for axial water transport and contribute to hydraulic architecture. Variation in metaxylem vessel size and number can impact drought tolerance in crop plants, including rice, a crop that is particularly sensitive to drought. Identifying and validating candidate genes for metaxylem development would aid breeding efforts for improved varieties for drought tolerance. We identified three transcription factor candidate genes that potentially regulate metaxylem vessel size and number in rice based on orthologous annotations, published expression data, and available root and drought-related QTL data. Single gene knockout mutants were generated for each candidate using CRISPR-Cas9 genome editing. Root metaxylem vessel area and number were analyzed in 6-week-old knockout mutants and wild-type plants under well-watered and drought conditions in the greenhouse. Compared with wild type, LONESOME HIGHWAY (OsLHW) mutants had fewer, smaller metaxylem vessels in shallow roots and more, larger vessels in deep roots in drought conditions, indicating that OsLHW may be a repressor of drought-induced metaxylem plasticity. The AUXIN RESPONSE FACTOR 15 mutants showed fewer but larger metaxylem vessel area in well-watered conditions, but phenotypes were inconsistent under drought treatment. ORYZA SATIVA HOMEBOX 6 (OSH6) mutants had fewer, smaller metaxylem vessels in well-watered conditions with greater effects on xylem number than size. OSH6 mutants had larger shoots and more, deeper roots than the wild type in well-watered conditions, but there were no differences in performance under drought between mutants and wild type. Though these candidate gene mutants did not exhibit large phenotypic effects, the identification and investigation of candidate genes related to metaxylem traits in rice deepen our understanding of metaxylem development and are needed to facilitate incorporation of favorable alleles into breeding populations to improve drought stress tolerance.

A Comparative Analysis of Transcription Networks Active in Juvenile and Mature Wood in Populus

Juvenile wood (JW) and mature wood (MW) have distinct physical and chemical characters, resulting from wood formation at different development phases over tree lifespan. However, the regulatory mechanisms that distinguish or modulate the characteristics of JW and MW in relation to each other have not been mapped. In this study, by employing the Populus trees with an identical genetic background, we carried out RNA sequencing (RNA-seq) and whole genome bisulfite sequencing (WGBS) in JW and MW forming tissue and analyzed the transcriptional programs in association with the wood formation in different phrases. JW and MW of Populus displayed different wood properties, including higher content of cellulose and hemicelluloses, less lignin, and longer and larger fiber cells and vessel elements in MW as compared with JW. Significant differences in transcriptional programs and patterns of DNA methylation were detected between JW and MW. The differences were concentrated in gene networks involved in regulating hormonal signaling pathways responsible for auxin distribution and brassinosteroids biosynthesis as well as genes active in regulating cell expansion and secondary cell wall biosynthesis. An observed correlation between gene expression profiling and DNA methylation indicated that DNA methylation affected expression of the genes related to auxin distribution and brassinosteroids signal transduction, cell expansion in JW, and MW formation. The results suggest that auxin distribution, brassinosteroids biosynthesis, and signaling be the critical molecular modules in formation of JW and MW. DNA methylation plays a role in formatting the molecular modules which contribute to the transcriptional programs of wood formation in different development phases. The study sheds light into better understanding of the molecular networks underlying regulation of wood properties which would be informative for genetic manipulation for improvement of wood formation.

Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

Characterization of a High Hierarchical Regulator, PtrGATA12, Functioning in Differentially Regulating Secondary Wall Component Biosynthesis in Populus trichocarpa

In plants, GATA transcription factors (TFs) have been reported to play vital roles in to a wide range of biological processes. To date, there is still no report about the involvement and functions of woody plant GATA TFs in wood formation. In this study, we described the functional characterization of a Populus trichocarpa GATA TF, PtrGATA12, which encodes a nuclear-localized transcriptional activator predominantly expressing in developing xylem tissues. Overexpression of PtrGATA12 not only inhibited growths of most phenotypic traits and biomass accumulation, but also altered the expressions of some master TFs and pathway genes involved in secondary cell wall (SCW) and programmed cell death, leading to alternated SCW components and breaking forces of stems of transgenic lines. The significant changes occurred in the contents of hemicellulose and lignin and SCW thicknesses of fiber and vessel that increased by 13.5 and 10.8%, and 20.83 and 11.83%, respectively. Furthermore, PtrGATA12 bound directly to the promoters of a battery of TFs and pathway genes and activated them; the binding sites include two cis-acting elements that were specifically enriched in their promoter regions. Taken together, our results suggest PtrGATA12, as a higher hierarchical TF on the top of PtrWND6A, PtrWND6B, PtrMYB152, and PtrMYB21, exert a coordinated regulation of SCW components biosynthesis pathways through directly and indirectly controlling master TFs, middle-level TFs, and further downstream pathway genes of the currently known hierarchical transcription network that governs SCW formation.

Identification and Characterization of Secondary Wall-Associated NAC Genes and Their Involvement in Hormonal Responses in Tobacco (Nicotiana tabacum)

Secondary wall-associated NAC (SWN) genes are a subgroup of NAC (NAM, ATAF, and CUC) transcription factors (TF) that play a key role in regulating secondary cell wall biosynthesis in plants. However, this gene family has not been systematically characterized, and their potential roles in response to hormones are unknown in Nicotiana tabacum. In this study, a total of 40 SWN genes, of which 12 from Nicotiana tomentosiformis, 13 from Nicotiana sylvestris, and 15 from Nicotiana tabacum, were successfully identified. The 15 SWNs from Nicotiana tabacum were further classified into three groups, namely, vascular-related NAC domain genes (NtVNDs), NAC secondary wall thickening promoting factor genes (NtNSTs), and secondary wall-associated NAC domain genes (NtSNDs). The protein characteristic, gene structure, and chromosomal location of 15 NtSWNs (also named Nt1 to Nt15) were also analyzed. The NtVND and NtNST group genes had five conserved subdomains in their N-terminal regions and a motif (LP[Q/x] L[E/x] S[P/A]) in their diverged C- terminal regions. Some hormones, dark and low-temperature related cis-acting elements, were significantly enriched in the promoters of NtSWN genes. A comprehensive expression profile analysis revealed that Nt4 and Nt12 might play a role in vein development. Others might be important for stem development. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) revealed that in the NtNST group, genes such as Nt7, Nt8, and Nt13 were more sensitive than the genes in NtVND and NtSND groups under abiotic stress conditions. A transactivation assay further suggested that Nt7, Nt8, and Nt13 showed a significant transactivation activity. Overall, SWN genes were finally identified and characterized in diploid and tetraploid tobacco, revealing new insights into their evolution, variation, and homology relationships. Transcriptome, cis-acting element, qRT-PCR, and transactivation assay analysis indicated the roles in hormonal and stress responses, which provided further resources in molecular mechanism and genetic improvement.

NST- and SND-subgroup NAC proteins coordinately act to regulate secondary cell wall formation in cotton

Secondary cell wall (SCW) has a strong impact on plant growth and adaptation to the environments. Previous studies have shown that NAC (NAM, ATAF1/2, and CUC2) transcription factors act as key regulators of SCW biosynthesis. However, the regulatory network triggered by NAC proteins is largely unknown, especially in cotton, a model plant for SCW development studies. Here, we show that several cotton NAC transcription factors are clustered in the same group with Arabidopsis secondary wall NACs (SWNs), including secondary wall-associated NAC domain protein1 (SND1) and NAC secondary wall thickening promoting factor1/2 (NST1/2), so we name these cotton orthologs as SND1s and NST1s. We found that simultaneous silencing of SND1s and NST1s led to severe xylem and phloem developmental defect in cotton stems, however silencing either SND1s or NST1s alone had no visible phenotype. Silencing both SND1s and NST1s but not one subgroup caused decreased expression of a set of SCW-associated genes, while over-expression of cotton SWNs in tobacco leaves resulted in SCW deposition. SWNs could bind the promoter of MYB46 and MYB83, which are highly expressed in SCW-rich tissues of cotton. In total, our data provide evidence that cotton SWNs positively and coordinately regulate SCW formation.

Comparative Analysis of the Lignification Process of Two Bamboo Shoots Stored at Room Temperature

Two types of bamboo shoots, high bamboo (Phyllostachys prominens) shoots (HBSes) and moso bamboo (Phyllostachys edulis) shoots (MBSes), underwent a fast post-harvest lignification process under room temperature storage. To explore the mechanism of lignification in two types of bamboo shoots after post-harvest during room temperature storage, the measurement of cell wall polymers (lignin and cellulose) and enzyme activities of phenylalanine ammonialyase (PAL) and peroxidase (POD), and relative expression of related transcription networks factors (TFs) were performed. The results suggested that the lignification process in HBSes is faster than that in MBSes because of incremental increase in lignin and cellulose contents within 6 days and the shorter shelf-life. Additionally, compared with the expression pattern of lignification-related TFs and correlation analysis of lignin and cellulose contents, MYB20, MYB43, MYB85 could function positively in the lignification process of two types of bamboo shoots. A negative regulator, KNAT7, could negatively regulate the lignin biosynthesis in two types of bamboo shoots. In addition, MYB63 could function positively in HBSes, and NST1 could function negatively in MBSes. Notably, MYB42 may function differently in the two types of bamboo shoots, that is, a positive regulator in HBSes, but a negative regulator in MBSes. Transcription networks provide a comprehensive analysis to explore the mechanism of lignification in two types of bamboo shoots after post-harvest during room temperature storage. These results suggest that the lignification of bamboo shoots was mainly due to the increased activity of POD, higher expression levels of MYB20, MYB43, MYB63, and MYB85 genes, and lower expression levels of KNAT7 and NST1 genes, and the lignification process of HBSes and MBSes had significant differences.

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

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

The chromatin-modifying protein HUB2 is involved in the regulation of lignin composition in xylem vessels

PIRIN2 (PRN2) was earlier reported to suppress syringyl (S)-type lignin accumulation of xylem vessels of Arabidopsis thaliana. In the present study, we report yeast two-hybrid results supporting the interaction of PRN2 with HISTONE MONOUBIQUITINATION2 (HUB2) in Arabidopsis. HUB2 has been previously implicated in several plant developmental processes, but not in lignification. Interaction between PRN2 and HUB2 was verified by β-galactosidase enzymatic and co-immunoprecipitation assays. HUB2 promoted the deposition of S-type lignin in the secondary cell walls of both stem and hypocotyl tissues, as analysed by pyrolysis-GC/MS. Chemical fingerprinting of individual xylem vessel cell walls by Raman and Fourier transform infrared microspectroscopy supported the function of HUB2 in lignin deposition. These results, together with a genetic analysis of the hub2 prn2 double mutant, support the antagonistic function of PRN2 and HUB2 in deposition of S-type lignin. Transcriptome analyses indicated the opposite regulation of the S-type lignin biosynthetic gene FERULATE-5-HYDROXYLASE1 by PRN2 and HUB2 as the underlying mechanism. PRN2 and HUB2 promoter activities co-localized in cells neighbouring the xylem vessel elements, suggesting that the S-type lignin-promoting function of HUB2 is antagonized by PRN2 for the benefit of the guaiacyl (G)-type lignin enrichment of the neighbouring xylem vessel elements.

Involvement of VNS NAC-domain transcription factors in tracheid formation in Pinus taeda

Vascular plants have two types of water-conducting cells, xylem vessel cells (in angiosperms) and tracheid cells (in ferns and gymnosperms). These cells are commonly characterized by secondary cell wall (SCW) formation and programmed cell death (PCD), which increase the efficiency of water conduction. The differentiation of xylem vessel cells is regulated by a set of NAC (NAM, ATAF1/2 and CUC2) transcription factors, called the VASCULAR-RELATED NAC-DOMAIN (VND) family, in Arabidopsis thaliana Linne. The VNDs regulate the transcriptional induction of genes required for SCW formation and PCD. However, information on the transcriptional regulation of tracheid cell differentiation is still limited. Here, we performed functional analysis of loblolly pine (Pinus taeda Linne) VND homologs (PtaVNS, for VND, NST/SND, SMB-related protein). We identified five PtaVNS genes in the loblolly pine genome, and four of these PtaVNS genes were highly expressed in tissues with tracheid cells, such as shoot apices and developing xylem. Transient overexpression of PtaVNS genes induced xylem vessel cell-like patterning of SCW deposition in tobacco (Nicotiana benthamiana Domin) leaves, and up-regulated the promoter activities of loblolly pine genes homologous to SCW-related MYB transcription factor genes and cellulose synthase genes, as well as to cysteine protease genes for PCD. Collectively, our data indicated that PtaVNS proteins possess transcriptional activity to induce the molecular programs required for tracheid formation, i.e., SCW formation and PCD. Moreover, these findings suggest that the VNS-MYB-based transcriptional network regulating water-conducting cell differentiation in angiosperm and moss plants is conserved in gymnosperms.

PIRIN2 suppresses S-type lignin accumulation in a noncell-autonomous manner in Arabidopsis xylem elements

PIRIN (PRN) genes encode cupin domain-containing proteins that function as transcriptional co-regulators in humans but that are poorly described in plants. A previous study in xylogenic cell cultures of Zinnia elegans suggested a role for a PRN protein in lignification. This study aimed to identify the function of Arabidopsis (Arabidopsis thaliana) PRN proteins in lignification of xylem tissues. Chemical composition of the secondary cell walls was analysed in Arabidopsis stems and/or hypocotyls by pyrolysis-gas chromatography/mass spectrometry, 2D-nuclear magnetic resonance and phenolic profiling. Secondary cell walls of individual xylem elements were chemotyped by Fourier transform infrared and Raman microspectroscopy. Arabidopsis PRN2 suppressed accumulation of S-type lignin in Arabidopsis stems and hypocotyls. PRN2 promoter activity and PRN2:GFP fusion protein were localised specifically in cells next to the vessel elements, suggesting a role for PRN2 in noncell-autonomous lignification of xylem vessels. Accordingly, PRN2 modulated lignin chemistry in the secondary cell walls of the neighbouring vessel elements. These results indicate that PRN2 suppresses S-type lignin accumulation in the neighbourhood of xylem vessels to bestow G-type enriched lignin composition on the secondary cell walls of the vessel elements. Gene expression analyses suggested that PRN2 function is mediated by regulation of the expression of the lignin-biosynthetic genes.

A NAC-type transcription factor confers aluminium resistance by regulating cell wall-associated receptor kinase 1 and cell wall pectin

Transcriptional regulation is important for plants to respond to toxic effects of aluminium (Al). However, our current knowledge to these events is confined to a few transcription factors. Here, we functionally characterized a rice bean (Vigna umbellata) NAC-type transcription factor, VuNAR1, in terms of Al stress response. We demonstrated that rice bean VuNAR1 is a nuclear-localized transcriptional activator, whose expression was specifically upregulated by Al in roots but not in shoot. VuNAR1 overexpressing Arabidopsis plants exhibit improved Al resistance via Al exclusion. However, VuNAR1-mediated Al exclusion is independent of the function of known Al-resistant genes. Comparative transcriptomic analysis revealed that VuNAR1 specifically regulates the expression of genes associated with protein phosphorylation and cell wall modification in Arabidopsis. Transient expression assay demonstrated the direct transcriptional activation of cell wall-associated receptor kinase 1 (WAK1) by VuNAR1. Moreover, yeast one-hybrid assays and MEME motif searches identified a new VuNAR1-specific binding motif in the promoter of WAK1. Compared with wild-type Arabidopsis plants, VuNAR1 overexpressing plants have higher WAK1 expression and less pectin content. Taken together, our results suggest that VuNAR1 regulates Al resistance by regulating cell wall pectin metabolism via directly binding to the promoter of WAK1 and induce its expression.

The Progression of Xylem Vessel Cell Differentiation is Dependent on the Activity Level of VND7 in Arabidopsis thaliana

Xylem vessels are important for water conduction in vascular plants. The VASCULAR-RELATED NAC-DOMAIN (VND) family proteins, master regulators of xylem vessel cell differentiation in Arabidopsis thaliana, can upregulate a set of genes required for xylem vessel cell differentiation, including those involved in secondary cell wall (SCW) formation and programmed cell death (PCD); however, it is not fully understood how VND activity levels influence these processes. Here, we examined the Arabidopsis VND7-VP16-GR line, in which VND7 activity is post-translationally activated by treatments with different concentrations of dexamethasone (DEX), a synthetic glucocorticoid. Our observations showed that 1 nM DEX induced weak SCW deposition, but not PCD, whereas 10 or 100 nM DEX induced both SCW deposition and PCD. The decreased chlorophyll contents and SCW deposition were apparent after 24 h of 100 nM DEX treatment, but became evident only after 48 h of 10 nM DEX treatment. Moreover, the lower DEX concentrations delayed the upregulation of VND7 downstream genes, and decreased their induction levels. They collectively suggest that the regulation of VND activity is important not only to initiate xylem vessel cell differentiation, but also regulate the quality of the xylem vessels through VND-activity-dependent upregulation of the PCD- and SCW-related genes.

MdMYB46 could enhance salt and osmotic stress tolerance in apple by directly activating stress-responsive signals

To expand the cultivation area of apple (Malus×domestica Borkh.) and select resistant varieties by genetic engineering, it is necessary to clarify the mechanism of salt and osmotic stress tolerance in apple. The MdMYB46 transcription factor was identified, and the stress treatment test of MdMYB46-overexpressing and MdMYB46-RNAi apple lines indicated that MdMYB46 could enhance the salt and osmotic stress tolerance in apple. In transgenic Arabidopsis and apple, MdMYB46 promoted the biosynthesis of secondary cell wall and deposition of lignin by directly binding to the promoter of lignin biosynthesis-related genes. To explore whether MdMYB46 could coordinate stress signal transduction pathways to cooperate with the formation of secondary walls to enhance the stress tolerance of plants, MdABRE1A, MdDREB2A and dehydration-responsive genes MdRD22 and MdRD29A were screened out for their positive correlation with osmotic stress, salt stress and the transcriptional level of MdMYB46. The further verification test demonstrated that MdMYB46 could activate their transcription by directly binding to the promoters of these genes. The above results indicate that MdMYB46 could enhance the salt and osmotic stress tolerance in apple not only by activating secondary cell wall biosynthesis pathways, but also by directly activating stress-responsive signals.

A Molecular Blueprint of Lignin Repression

Although lignin is essential to ensure the correct growth and development of land plants, it may be an obstacle to the production of lignocellulosics-based biofuels, and reduces the nutritional quality of crops used for human consumption or livestock feed. The need to tailor the lignocellulosic biomass for more efficient biofuel production or for improved plant digestibility has fostered considerable advances in our understanding of the lignin biosynthetic pathway and its regulation. Most of the described regulators are transcriptional activators of lignin biosynthesis, but considerably less attention has been devoted to the repressors of this pathway. We provide a comprehensive overview of the molecular factors that negatively impact on the lignification process at both the transcriptional and post-transcriptional levels.

Genome-wide identification and characterization of TALE superfamily genes in cotton reveals their functions in regulating secondary cell wall biosynthesis

BACKGROUND

Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.).

RESULTS

In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers.

CONCLUSION

We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

Genome-Wide Investigation of the NAC Gene Family and Its Potential Association with the Secondary Cell Wall in Moso Bamboo

NAC (NAM, ATAF, and CUC) transcription factors (TFs) are implicated in the transcriptional regulation of diverse processes and have been characterized in a number of plant species. However, NAC TFs are still not well understood in bamboo, especially their potential association with the secondary cell wall (SCW). Here, 94 PeNACs were identified and characterized in moso bamboo (Phyllostachys edulis). Based on their gene structures and conserved motifs, the PeNACs were divided into 11 groups according to their homologs in Arabidopsis. PeNACs were expressed variously in different tissues of moso bamboo, suggesting their functional diversity. Fifteen PeNACs associated with the SCW were selected for co-expression analysis and validation. It was predicted that 396 genes were co-expressed with the 15 PeNACs, in which 16 and 55 genes were involved in the lignin catabolic process and cellulose biosynthetic process respectively. As the degree of lignification in the growing bamboo shoots increased, all 15 PeNACs were upregulated with a trend of rising first and then decreasing except PeNAC37, which increased continuously. These results indicated that these PeNACs might play important roles in SCW biosynthesis and lignification in bamboo shoots. Seven of 15 PeNACs had been found positively co-expressed with seven PeMYBs, and they had similar expression patterns with those of the PeMYBs in bamboo shoots. The targeted sites of miR164 were found in 16 PeNACs, of which three PeNACs associated with SCW were validated to have an opposite expression trend to that of miR164 in growing bamboo shoots. In addition, three PeNACs were selected and verified to have self-activation activities. These results provide comprehensive information of the NAC gene family in moso bamboo, which will be helpful for further functional studies of PeNACs to reveal the molecular regulatory mechanisms of bamboo wood property.

Affinity-based high-resolution analysis of DNA binding by VASCULAR-RELATED NAC-DOMAIN7 via fluorescence correlation spectroscopy

VASCULAR-RELATED NAC-DOMAIN7 (VND7) is the master transcription factor for vessel element differentiation in Arabidopsis thaliana. To identify the cis-acting sequence(s) bound by VND7, we employed fluorescence correlation spectroscopy (FCS) to find VND7-DNA interactions quantitatively. This identified an 18-bp sequence from the promoter of XYLEM CYSTEINE PEPTIDASE1 (XCP1), a direct target of VND7. A quantitative assay for binding affinity between VND7 and the 18-bp sequence revealed the core nucleotides contributing to specific binding between VND7 and the 18-bp sequence. Moreover, by combining the systematic evolution of ligands by exponential enrichment (SELEX) technique with known consensus sequences, we defined a motif termed the Ideal Core Structure for binding by VND7 (ICSV). We also used FCS to search for VND7 binding sequences in the promoter regions of other direct targets. Taking these data together, we proposed that VND7 preferentially binds to the ICSV sequence. Additionally, we found that substitutions among the core nucleotides affected transcriptional regulation by VND7 in vivo, indicating that the core nucleotides contribute to vessel-element-specific gene expression. Furthermore, our results demonstrate that FCS is a powerful tool for unveiling the DNA-binding properties of transcription factors.

Analysis of Orthologous SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) Promotor Activity in Herbaceous and Woody Angiosperms

SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) is a master regulator of fibre secondary wall deposition in Arabidopsis thaliana (Arabidopsis), with homologs in other angiosperms and gymnosperms. However, it is poorly understood to what extent the fibre-specific regulation of the SND1 promoter, and that of its orthologs, is conserved between diverged herbaceous and woody lineages. We performed a reciprocal reporter gene analysis of orthologous SND1 promoters from Arabidopsis (AthSND1), Eucalyptus grandis (EgrNAC61) and Populus alba × P. grandidentata (PagWND1A) relative to secondary cell wall-specific Cellulose Synthase4 (CesA4) and CesA7 promoters, in both a non-woody (Arabidopsis) and a woody (poplar) system. β-glucuronidase (GUS) reporter analysis in Arabidopsis showed that the SND1 promoter was active in vascular tissues as previously reported and showed interfascicular and xylary fibre-specific expression in inflorescence stems, while reporter constructs of the woody plant-derived promoters were partial to the (pro)cambium-phloem and protoxylem. In transgenic P. tremula × P. alba plants, all three orthologous SND1 promoters expressed the GUS reporter similarly and preferentially in developing secondary xylem, ray parenchyma and cork cambium. Ours is the first study to reciprocally test orthologous SND1 promoter specificity in herbaceous and woody species, revealing diverged regulatory functions in the herbaceous system.

Understanding the molecular mechanisms underlying graft success in grapevine

BACKGROUND

Grafting is an intensive commercial practice required to protect the European grapevine against the Phylloxera pest. Rootstocks resistant to this pest are hybrids of American vine species with different levels of compatibility with European Vitis vinifera varieties. Aiming to understand what drives grafting compatibility in grapevine, a transcriptomic approach was used to search for master regulators of graft success. Two scion/rootstock combinations, with different levels of compatibility, were compared in a nursery-grafting context at two stages, at 21 and 80 days after grafting.

RESULTS

In the most compatible combination, an earlier and higher expression of genes signaling the metabolic and hormonal pathways as well as a reduced expression of genes of the phenolic metabolism and of the oxidative stress response was observed. At 80 days after grafting a higher expression of transcription factors regulating vascular maintenance, differentiation and proliferation was obtained in the most compatible combination. Moreover, lower expression levels of microRNAs potentially targeting important transcription factors related to plant development was observed in the more compatible combination when compared to the less compatible one.

CONCLUSION

In this context, a set of regulators was selected as potential expression markers for early prediction of a compatible grafting.

Understanding the molecular mechanisms underlying graft success in grapevine

BACKGROUND

Grafting is an intensive commercial practice required to protect the European grapevine against the Phylloxera pest. Rootstocks resistant to this pest are hybrids of American vine species with different levels of compatibility with European Vitis vinifera varieties. Aiming to understand what drives grafting compatibility in grapevine, a transcriptomic approach was used to search for master regulators of graft success. Two scion/rootstock combinations, with different levels of compatibility, were compared in a nursery-grafting context at two stages, at 21 and 80 days after grafting.

RESULTS

In the most compatible combination, an earlier and higher expression of genes signaling the metabolic and hormonal pathways as well as a reduced expression of genes of the phenolic metabolism and of the oxidative stress response was observed. At 80 days after grafting a higher expression of transcription factors regulating vascular maintenance, differentiation and proliferation was obtained in the most compatible combination. Moreover, lower expression levels of microRNAs potentially targeting important transcription factors related to plant development was observed in the more compatible combination when compared to the less compatible one.

CONCLUSION

In this context, a set of regulators was selected as potential expression markers for early prediction of a compatible grafting.

Identification and functional evaluation of accessible chromatin associated with wood formation in Eucalyptus grandis

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue.

The flowering hormone florigen accelerates secondary cell wall biogenesis to harmonize vascular maturation with reproductive development

Florigen, a proteinaceous hormone, functions as a universal long-range promoter of flowering and concurrently as a generic growth-attenuating hormone across leaf and stem meristems. In flowering plants, the transition from the vegetative phase to the reproductive phase entails the orchestration of new growth coordinates and a global redistribution of resources, signals, and mechanical loads among organs. However, the ultimate cellular processes governing the adaptation of the shoot system to reproduction remain unknown. We hypothesized that if the mechanism for floral induction is universal, then the cellular metabolic mechanisms underlying the conditioning of the shoot system for reproduction would also be universal and may be best regulated by florigen itself. To understand the cellular basis for the vegetative functions of florigen, we explored the radial expansion of tomato stems. RNA-Seq and complementary genetic and histological studies revealed that florigen of endogenous, mobile, or induced origins accelerates the transcription network navigating secondary cell wall biogenesis as a unit, promoting vascular maturation and thereby adapting the shoot system to the developmental needs of the ensuing reproductive phase it had originally set into motion. We then demonstrated that a remarkably stable and broadly distributed florigen promotes MADS and MIF genes, which in turn regulate the rate of vascular maturation and radial expansion of stems irrespective of flowering or florigen level. The dual acceleration of flowering and vascular maturation by florigen provides a paradigm for coordinated regulation of independent global developmental programs.

A systems biology view of wood formation in Eucalyptus grandis trees submitted to different potassium and water regimes

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

Modified expression of ZmMYB167 in Brachypodium distachyon and Zea mays leads to increased cell wall lignin and phenolic content

One of the challenges to enable targeted modification of lignocellulosic biomass from grasses for improved biofuel and biochemical production lies within our limited understanding of the transcriptional control of secondary cell wall biosynthesis. Here, we investigated the role of the maize MYB transcription factor ZmMYB167 in secondary cell wall biosynthesis and how modified ZmMYB167 expression in two distinct grass model species affects plant biomass and growth phenotypes. Heterologous expression of ZmMYB167 in the C3 model system Brachypodium led to mild dwarf phenotypes, increased lignin (~7% to 13%) and S-lignin monomer (~11% to 16%) content, elevated concentrations of cell wall-bound p-coumaric acid (~15% to 24%) and reduced biomass sugar release (~20%) compared to controls. Overexpression of ZmMYB167 in the C4 model system Zea mays increased lignin (~4% to 13%), p-coumaric acid (~8% to 52%) and ferulic acid (~13% to 38%) content but did not affect plant growth and development nor biomass recalcitrance. Taken together, modifying ZmMYB167 expression represents a target to alter lignin and phenolic content in grasses. The ZmMYB167 expression-induced discrepancies in plant phenotypic and biomass properties between the two grass model systems highlight the challenges and opportunities for MYB transcription factor-based genetic engineering approaches of grass biomass.

Organellar carbon metabolism is coordinated with distinct developmental phases of secondary xylem

Subcellular compartmentation of plant biosynthetic pathways in the mitochondria and plastids requires coordinated regulation of nuclear encoded genes, and the role of these genes has been largely ignored by wood researchers. In this study, we constructed a targeted systems genetics coexpression network of xylogenesis in Eucalyptus using plastid and mitochondrial carbon metabolic genes and compared the resulting clusters to the aspen xylem developmental series. The constructed network clusters reveal the organization of transcriptional modules regulating subcellular metabolic functions in plastids and mitochondria. Overlapping genes between the plastid and mitochondrial networks implicate the common transcriptional regulation of carbon metabolism during xylem secondary growth. We show that the central processes of organellar carbon metabolism are distinctly coordinated across the developmental stages of wood formation and are specifically associated with primary growth and secondary cell wall deposition. We also demonstrate that, during xylogenesis, plastid-targeted carbon metabolism is partially regulated by the central clock for carbon allocation towards primary and secondary xylem growth, and we discuss these networks in the context of previously established associations with wood-related complex traits. This study provides a new resolution into the integration and transcriptional regulation of plastid- and mitochondrial-localized carbon metabolism during xylogenesis.

A Standardized Synthetic Eucalyptus Transcription Factor and Promoter Panel for Re-engineering Secondary Cell Wall Regulation in Biomass and Bioenergy Crops

Re-engineering of transcriptional networks regulating secondary cell wall formation may allow the improvement of plant biomass in widely grown plantation crops such as Eucalyptus. However, there is currently a scarcity of freely available standardized biological parts (e.g., Phytobricks) compatible with Type IIS assembly approaches from forest trees, and there is a need to accelerate transcriptional network inference in nonmodel biomass crops. Here we describe the design and synthesis of a versatile three-panel biological parts collection of 221 secondary cell wall-related Eucalyptus grandis transcription factor coding sequences and 65 promoters that are compatible with GATEWAY, Golden Gate, MoClo, and GoldenBraid DNA assembly methods and generally conform to accepted Phytobrick syntaxes. This freely available resource is intended to accelerate synthetic biology applications in multiple plant biomass crops and enable reconstruction of secondary cell wall transcriptional networks using high-throughput assays such as DNA affinity purification sequencing (DAP-seq) and enhanced yeast one-hybrid (eY1H) screening.

Genome-wide identification and expression analysis of the MYB transcription factor in moso bamboo (Phyllostachys edulis)

The MYB family, one of the largest transcription factor (TF) families in the plant kingdom, plays vital roles in cell formation, morphogenesis and signal transduction, as well as responses to biotic and abiotic stresses. However, the underlying function of bamboo MYB TFs remains unclear. To gain insight into the status of these proteins, a total of 85 PeMYBs, which were further divided into 11 subgroups, were identified in moso bamboo (Phyllostachys edulis) by using a genome-wide search strategy. Gene structure analysis showed that PeMYBs were significantly different, with exon numbers varying from 4 to 13. Phylogenetic analysis indicated that PeMYBs clustered into 27 clades, of which the function of 18 clades has been predicted. In addition, almost all of the PeMYBs were differently expressed in leaves, panicles, rhizomes and shoots based on RNA-seq data. Furthermore, qRT-PCR analysis showed that 12 PeMYBs related to the biosynthesis and deposition of the secondary cell wall (SCW) were constitutively expressed, and their transcript abundance levels have changed significantly with increasing height of the bamboo shoots, for which the degree of lignification continuously increased. This result indicated that these PeMYBs might play fundamental roles in SCW thickening and bamboo shoot lignification. The present comprehensive and systematic study on the members of the MYB family provided a reference and solid foundation for further functional analysis of MYB TFs in moso bamboo.

Transcriptome Comparison Reveals Distinct Selection Patterns in Domesticated and Wild Agave Species, the Important CAM Plants

Agave species are an important family of crassulacean acid metabolism (CAM) plants with remarkable tolerance to heat and drought stresses (Agave deserti) in arid regions and multiple agricultural applications, such as spirit (Agave tequilana) and fiber (Agave sisalana) production. The agave genomes are commonly too large to sequence, which has significantly restricted our understanding to the molecular basis of stress tolerance and economic traits in agaves. In this study, we collected three transcriptome databases for comparison to reveal the phylogenic relationships and evolution patterns of the three agave species. The results indicated the close but distinctly domesticated relations between A. tequilana and A. sisalana. Natural abiotic and biotic selections are very important factors that have contributed to distinct economic traits in agave domestication together with artificial selection. Besides, a series of candidate unigenes regulating fructan, fiber, and stress response-related traits were identified in A. tequilana, A. sisalana, and A. deserti, respectively. This study represents the first transcriptome comparison within domesticated and wild agaves, which would serve as a guidance for further studies on agave evolution, environmental adaptation, and improvement of economically important traits.

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.

The quest for transcriptional hubs of lignin biosynthesis: beyond the NAC-MYB-gene regulatory network model

Lignin is an important secondary metabolite in plants. The biosynthesis of lignin is initiated by the transcriptional upregulation of genes encoding enzymes involved in monolignol biosynthesis and lignin polymerization. Based on studies of xylem differentiation over the last two decades, the NAC-MYB-based gene regulatory network (NAC-MYB-GRN) model is widely considered to underpin developmental lignin biosynthesis. We are now standing on the threshold of a new direction in transcriptional regulation research; the search for novel molecular hubs connecting developmental/environmental signals in lignin biosynthesis. Emerging genome-wide 'omics' technologies are a promising approach for understanding such hubs. Elucidating these molecular hubs may enable us to control lignification in harmony with plant development and environmental adaptation.

Blue Light Regulates Secondary Cell Wall Thickening via MYC2/MYC4 Activation of the NST1-Directed Transcriptional Network in Arabidopsis

Secondary cell walls (SCWs) are formed in some specific types of plant cells, providing plants with mechanical strength. During plant growth and development, formation of secondary cell walls is regulated by various developmental and environmental signals. The underlying molecular mechanisms are poorly understood. In this study, we analyzed the blue light receptor cryptochrome1 (cry1) mutant of Arabidopsis thaliana for its SCW phenotypes. During inflorescence stem growth, SCW thickening in the vasculature was significantly affected by blue light. cry1 plants displayed a decline of SCW thickening in fiber cells, while CRY1 overexpression led to enhanced SCW formation. Transcriptome analysis indicated that the reduced SCW thickening was associated with repression of the NST1-directed transcription regulatory networks. Further analyses revealed that the expression of MYC2/MYC4 that is induced by blue light activates the transcriptional network underlying SCW thickening. The activation is caused by direct binding of MYC2/MYC4 to the NST1 promoter. This study demonstrates that SCW thickening in fiber cells is regulated by a blue light signal that is mediated through MYC2/MYC4 activation of NST1-directed SCW formation in Arabidopsis.

Temporal analysis of Arabidopsis genes activated by Eucalyptus grandis NAC transcription factors associated with xylem fibre and vessel development

Secondary cell wall (SCW) deposition in Arabidopsis is regulated among others by NAC transcription factors, where SND1 chiefly initiates xylem fibre differentiation while VND6 controls metaxylem vessel SCW development, especially programmed cell death and wall patterning. The translational relevance of Arabidopsis SCW regulation theory and the utility of characterized transcription factors as modular synthetic biology tools for improving commercial fibre crops is unclear. We investigated inter-lineage gene activation dynamics for potential fibre and vessel differentiation regulators from the widely grown hardwood Eucalyptus grandis (Myrtales). EgrNAC26, a VND6 homolog, and EgrNAC61, an SND1 homolog, were transiently expressed in Arabidopsis mesophyll protoplasts in parallel to determine early and late (i.e. 7 and 14 hours post-transfection) gene targets. Surprisingly, across the time series EgrNAC26 activated only a subset of SCW-related transcription factors and biosynthetic genes activated by EgrNAC61, specializing instead in targeting vessel-specific wall pit and programmed cell death markers. Promoters of EgrNAC26 and EgrNAC61 both induced reporter gene expression in vessels of young Arabidopsis plants, with EgrNAC61 also conferring xylem- and cork cambium-preferential expression in Populus. Our results demonstrate partial conservation, with notable exceptions, of SND1 and VND6 homologs in Eucalyptus and a first report of cork cambium expression for EgrNAC61.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.