Arabidopsis VASCULAR-RELATED NAC-DOMAIN6 directly regulates the genes that govern programmed cell death and secondary wall formation during xylem differentiation

VND Genes Redundantly Regulate Cell Wall Thickening during Parasitic Nematode Infection

Plant parasitic root-knot nematodes are major agricultural pests worldwide, as they infect plant roots and cause substantial damages to crop plants. Root-knot nematodes induce specialized feeding cells known as giant cells (GCs) in the root vasculature, which serve as nutrient reservoirs for the infecting nematodes. Here, we show that the cell walls of GCs thicken to form pitted patterns that superficially resemble metaxylem cells. Interestingly, VASCULAR-RELATED NAC-DOMAIN1 (VND1) was found to be upregulated, while the xylem-type programmed cell death marker XYLEM CYSTEINE PEPTIDASE 1 was downregulated upon nematode infection. The vnd2 and vnd3 mutants showed reduced secondary cell wall pore size, while the vnd1 vnd2 vnd3 triple mutant produced significantly fewer nematode egg masses when compared with the wild type. These results suggest that the GC development pathway likely shares common signaling modules with the metaxylem differentiation pathway and VND1, VND2, and VND3 redundantly regulate plant-nematode interaction through secondary cell wall formation.

Aspartic proteases gene family: Identification and expression profiles during stem vascular development in tobacco

Aspartic proteases (APs) constitute a large family in plants and are widely involved in diverse biological processes, like chloroplast metabolism, biotic and abiotic stress responses, and reproductive development. In this study, we focused on overall analysis of the APs genes in tobacco. Our analysis included the phylogeny and cis-elements in the cell wall-associated promoters of these genes. To characterize the expression patterns of APs genes in stem vascular development. The tissue expression analysis showed that NtAED3-like was preferentially expressed in the differentiating xylem and phloem cells of the vascular system. Based on histochemical staining analysis showed that the NtAED3-like gene was specifically expressed in stem vascular tissue, root vascular tissue, and petiole vascular tissue. The TdT-mediated dUTP nick-end labeling (TUNEL) assay illustrated a delayed progression of programmed cell death (PCD) within the xylem of the ko-ntaed3a-like mutant, relative to the wild type. The mutant ko-ntaed3a-like exhibited a phenotype of thinning stem circumference and changed in xylem structure and lignin content. In addition, the two-dimension heteronuclear single quantum coherent nuclear magnetic resonance (2D-HSQC) analysis of three milled wood lignins (MWLs) showed that the content of β-O-4 connection in ko-ntaed3a-like decreased slightly compared with wild type. In conclusion, this study provides our understanding of the regulation of vascular tissue development by the NtAED3-like gene in tobacco and provides a better basis for determining the molecular mechanism of the aspartic protease in secondary cell wall (SCW) development.

LONESOME HIGHWAY-TARGET OF MONOPTEROS5 transcription factor complex promotes a predifferentiation state for xylem vessel differentiation in the root apical meristem by inducing the expression of VASCULAR-RELATED NAC-DOMAIN genes

In Arabidopsis thaliana, heterodimers comprising two bHLH family proteins, LONESOME HIGHWAY (LHW) and TARGET OF MONOPTEROS5 (TMO5) or its homolog TMO5-LIKE 1 (T5L1) control vascular development in the root apical meristem (RAM). The LHW-TMO5/T5L1 complex regulates vascular cell proliferation, vascular pattern organization, and xylem vessel differentiation; however, the mechanism of preparation for xylem vessel differentiation in the RAM remains elusive. We examined the relationship between LHW-T5L1 and VASCULAR-RELATED NAC-DOMAIN (VND) genes, which are key regulators of vessel differentiation, using reverse genetics approaches. LHW-T5L1 upregulated the expression of VND1, VND2, VND3, VND6, and VND7 but not that of other VNDs. The expression of VND1-VND3 in the RAM was decreased in lhw. In vnd1 vnd2 vnd3 triple loss-of-function mutant roots, metaxylem differentiation was delayed, and VND6 and VND7 expression was reduced. Furthermore, transcriptome analysis of VND1-overexpressing cells revealed that VND1 upregulates genes involved in the synthesis of secondary cell wall components. These results suggest that LHW-T5L1 upregulates VND1-VND3 at the early stages of vascular development in the RAM, and VNDs promote a predifferentiation state for xylem vessels by triggering low levels of VND6 and VND7 as well as genes for the synthesis of secondary cell wall materials.

Effect of potassium fertilization on storage root number, yield, and appearance quality of sweet potato (Ipomoea batatas L.)

Increasing storage root number is a pivotal approach to enhance both storage root (SR) yield and appearance quality of sweet potato. Here, 2-year field experiments were conducted to investigate the effect of 0 (K0), 120 (K1), 240 (K2), and 360 (K3) kg ha-1 potassium fertilizer (K2O) on lignin metabolism, root growth, storage root yield, and uniformity. The results demonstrated that potassium (K) application led to a decrease in the activities of key enzymes involved in lignin biosynthesis, including phenylalanine deaminase (PAL), 4-coumarate coenzyme A ligase (4-CL), cinnamic acid dehydrogenase (CAD), polyphenol oxidase (PPO), and peroxidase (POD). This resulted in a significant reduction in lignin and G-type lignin contents in potential SRs compared to K0 treatment within 10-30 days after planting (DAP). BJ553 exhibited a significant decrease in PAL activity, as well as lignin and G-type contents at 10 DAP, whereas YS25 showed delayed effects until 20 DAP. However, the number and distribution of secondary xylem conduits as well as the mid-column diameter area in roots were increased in K2 treatment. Interestingly, K2 treatment exhibited significantly larger potential SR diameter than other treatments at 15, 20, and 25 DAP. At harvest, K2 treatment increased the SR number, the single SR weight, and overall yield greatly compared with K0 treatment, with an average increase of 19.12%, 16.54%, and 16.92% respectively. The increase of SR number in BJ553 was higher than that of YS25. Furthermore, K2 treatment exhibited the lowest coefficient of variation for both SR length and diameter, indicating a higher yield of middle-sized SRs. In general, appropriate potassium application could effectively suppress lignin biosynthesis, leading to a reduction in the degree of pericycle lignification in potential SRs. This promotes an increase in the number of storage roots and ultimately enhances both yield and appearance quality of sweet potato. The effect of potassium fertilizer on lignin metabolism in BJ553 roots was earlier and resulted in a greater increase in the SR number compared to YS25.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

The Carbon Flow Shifts from Primary to Secondary Metabolism during Xylem Vessel Cell Differentiation in Arabidopsis thaliana

Xylem vessel cell differentiation is characterized by the deposition of a secondary cell wall (SCW) containing cellulose, hemicellulose and lignin. VASCULAR-RELATED NAC-DOMAIN7 (VND7), a plant-specific NAC (NAM, ATAF1/2, and CUC2) transcription factor, is a master regulator of xylem vessel cell differentiation in Arabidopsis (Arabidopsis thaliana). Previous metabolome analysis using the VND7-inducible system in tobacco BY-2 cells successfully revealed significant quantitative changes in primary metabolites during xylem vessel cell differentiation. However, the flow of primary metabolites is not yet well understood. Here, we performed a metabolomic analysis of VND7-inducible Arabidopsis T87 suspension cells. Capillary electrophoresis-time-of-flight mass spectrometry quantified 57 metabolites, and subsequent data analysis highlighted active changes in the levels of UDP-glucose and phenylalanine, which are building blocks of cellulose and lignin, respectively. In a metabolic flow analysis using stable carbon 13 (13C) isotope, the 13C-labeling ratio specifically increased in 3-phosphoglycerate after 12 h of VND7 induction, followed by an increase in shikimate after 24 h of induction, while the inflow of 13C into lactate from pyruvate was significantly inhibited, indicating an active shift of carbon flow from glycolysis to the shikimate pathway during xylem vessel cell differentiation. In support of this notion, most glycolytic genes involved in the downstream of glyceraldehyde 3-phosphate were downregulated following the induction of xylem vessel cell differentiation, whereas genes for the shikimate pathway and phenylalanine biosynthesis were upregulated. These findings provide evidence for the active shift of carbon flow from primary metabolic pathways to the SCW polymer biosynthetic pathway at specific points during xylem vessel cell differentiation.

Sucrose Signaling Contributes to the Maintenance of Vascular Cambium by Inhibiting Cell Differentiation

Plants produce sugars by photosynthesis and use them for growth and development. Sugars are transported from source-to-sink organs via the phloem in the vasculature. It is well known that vascular development is precisely controlled by plant hormones and peptide hormones. However, the role of sugars in the regulation of vascular development is poorly understood. In this study, we examined the effects of sugars on vascular cell differentiation using a vascular cell induction system named 'Vascular Cell Induction Culture System Using Arabidopsis Leaves' (VISUAL). We found that sucrose has the strongest inhibitory effect on xylem differentiation, among several types of sugars. Transcriptome analysis revealed that sucrose suppresses xylem and phloem differentiation in cambial cells. Physiological and genetic analyses suggested that sucrose might function through the BRI1-EMS-SUPPRESSOR1 transcription factor, which is the central regulator of vascular cell differentiation. Conditional overexpression of cytosolic invertase led to a decrease in the number of cambium layers due to an imbalance between cell division and differentiation. Taken together, our results suggest that sucrose potentially acts as a signal that integrates environmental conditions with the developmental program.

NACs, generalist in plant life

Plant-specific NAC proteins constitute a major transcription factor family that is well-known for its roles in plant growth, development, and responses to abiotic and biotic stresses. In recent years, there has been significant progress in understanding the functions of NAC proteins. NAC proteins have a highly conserved DNA-binding domain; however, their functions are diverse. Previous understanding of the structure of NAC transcription factors can be used as the basis for their functional diversity. NAC transcription factors consist of a target-binding domain at the N-terminus and a highly versatile C-terminal domain that interacts with other proteins. A growing body of research on NAC transcription factors helps us comprehend the intricate signalling network and transcriptional reprogramming facilitated by NAC-mediated complexes. However, most studies of NAC proteins have been limited to a single function. Here, we discuss the upstream regulators, regulatory components and targets of NAC in the context of their prospective roles in plant improvement strategies via biotechnology intervention, highlighting the importance of the NAC transcription factor family in plants and the need for further research.

UBC6, a ubiquitin-conjugating enzyme, participates in secondary cell wall thickening in the inflorescence stem of Arabidopsis

Secondary cell wall (SCW) thickening in plant inflorescence stems is a complicated cellular process that is essential for stem strength and biomass. Although Arabidopsis NAC transcription factor (TF) 1 (NST1) regulates the SCW thickening in anther walls, the single T-DNA-insertion mutant (nst1) does not show disrupted SCW thickening in anther endothecium, interfascicular fibers or xylem. To better understand the regulatory mechanism of this process, we generated an ethyl methanesulfonate (EMS)-mutagenized Arabidopsis population with the nst1 background. scd5 (SCW-defective mutant 5) was isolated in a forward genetic screen from the EMS mutant library, which displayed not only less lignin deposition in the interfascicular fiber and xylem than the wild type but also a pendent inflorescence stem. The EMS-induced mutation associated with the scd5 phenotype was found in the 5th exon of At2G46030 that encodes a ubiquitin-conjugating enzyme (UBC6), we thereby renamed the allele nst1 ubc6. Overexpressing UBC6 in nst1 ubc6 rescued the defective SCW, whereas disrupting UBC6 in nst1 by the CRISPR/Cas9 system caused a phenotype similar to that observed in nst1 ubc6. UBC6 was localized to the nucleus and plasma membrane, and possessed E2 ubiquitin-conjugating activity in vitro. MYB7 and MYB32 are considered as transcription repressors in the phenylpropanoid pathway and are involved in NAC TF-related transcriptional regulation in SCW thickening. UBC6 can interact with MYB7 and MYB32 and positively mediate the degradation of MYB7 and MYB32 by the 26S proteasome. Overall, these results indicated the contribution of UBC6 to SCW thickening in Arabidopsis inflorescence stems.

How to explore what is hidden? A review of techniques for vascular tissue expression profile analysis

The evolution of plants to efficiently transport water and assimilates over long distances is a major evolutionary success that facilitated their growth and colonization of land. Vascular tissues, namely xylem and phloem, are characterized by high specialization, cell heterogeneity, and diverse cell components. During differentiation and maturation, these tissues undergo an irreversible sequence of events, leading to complete protoplast degradation in xylem or partial degradation in phloem, enabling their undisturbed conductive function. Due to the unique nature of vascular tissue, and the poorly understood processes involved in xylem and phloem development, studying the molecular basis of tissue differentiation is challenging. In this review, we focus on methods crucial for gene expression research in conductive tissues, emphasizing the importance of initial anatomical analysis and appropriate material selection. We trace the expansion of molecular techniques in vascular gene expression studies and discuss the application of single-cell RNA sequencing, a high-throughput technique that has revolutionized transcriptomic analysis. We explore how single-cell RNA sequencing will enhance our knowledge of gene expression in conductive tissues.

Endosperm cell death promoted by NAC transcription factors facilitates embryo invasion in Arabidopsis

In flowering plants, two fertilization products develop within the limited space of the seed: the embryo and the surrounding nutritive endosperm. The final size of the endosperm is modulated by the degree of embryo growth. In Arabidopsis thaliana, the endosperm expands rapidly after fertilization, but later gets invaded by the embryo that occupies most of the seed volume at maturity, surrounded by a single remaining aleurone-like endosperm layer.1,2,3,4 Embryo invasion is facilitated by the endosperm-expressed bHLH-type transcription factor ZHOUPI, which promotes weakening of endosperm cell walls.5,6 Endosperm elimination in zou mutants is delayed, and embryo growth is severely affected; the endosperm finally collapses around the dwarf embryo, causing the shriveled appearance of mature zou seeds.5,6,7 However, whether ZHOUPI facilitates mechanical endosperm destruction by the invading embryo or whether an active programmed cell death (PCD) process causes endosperm elimination has been subject to debate.2,8 Here we show that developmental PCD controlled by multiple NAC transcription factors in the embryo-adjacent endosperm promotes gradual endosperm elimination. Misexpressing the NAC transcription factor KIRA1 in the entire endosperm caused total endosperm elimination, generating aleurone-less mature seeds. Conversely, dominant and recessive higher-order NAC mutants led to delayed endosperm elimination and impaired cell corpse clearance. Promoting PCD in the zhoupi mutant partially rescued its embryo growth defects, while the endosperm in a zhoupi nac higher-order mutant persisted until seed desiccation. These data suggest that a combination of cell wall weakening and PCD jointly facilitates embryo invasion by an active auto-elimination of endosperm cells.

Comparative Transcriptome Analysis Reveals the Effect of Lignin on Storage Roots Formation in Two Sweetpotato (Ipomoea batatas (L.) Lam.) Cultivars

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important crops with high storage roots yield. The formation and expansion rate of storage root (SR) plays a crucial role in the production of sweet potato. Lignin affects the SR formation; however, the molecular mechanisms of lignin in SR development have been lacking. To reveal the problem, we performed transcriptome sequencing of SR harvested at 32, 46, and 67 days after planting (DAP) to analyze two sweet potato lines, Jishu25 and Jishu29, in which SR expansion of Jishu29 was early and had a higher yield. A total of 52,137 transcripts and 21,148 unigenes were obtained after corrected with Hiseq2500 sequencing. Through the comparative analysis, 9577 unigenes were found to be differently expressed in the different stages in two cultivars. In addition, phenotypic analysis of two cultivars, combined with analysis of GO, KEGG, and WGCNA showed the regulation of lignin synthesis and related transcription factors play a crucial role in the early expansion of SR. The four key genes swbp1, swpa7, IbERF061, and IbERF109 were proved as potential candidates for regulating lignin synthesis and SR expansion in sweet potato. The data from this study provides new insights into the molecular mechanisms underlying the impact of lignin synthesis on the formation and expansion of SR in sweet potatoes and proposes several candidate genes that may affect sweet potato yield.

Nicotiana benthamiana XYLEM CYSTEINE PROTEASE genes facilitate tracheary element formation in interfamily grafting

Grafting is a plant propagation technique widely used in agriculture. A recent discovery of the capability of interfamily grafting in Nicotiana has expanded the potential combinations of grafting. In this study, we showed that xylem connection is essential for the achievement of interfamily grafting and investigated the molecular basis of xylem formation at the graft junction. Transcriptome and gene network analyses revealed gene modules for tracheary element (TE) formation during grafting that include genes associated with xylem cell differentiation and immune response. The reliability of the drawn network was validated by examining the role of the Nicotiana benthamiana XYLEM CYSTEINE PROTEASE (NbXCP) genes in TE formation during interfamily grafting. Promoter activities of NbXCP1 and NbXCP2 genes were found in differentiating TE cells in the stem and callus tissues at the graft junction. Analysis of a Nbxcp1;Nbxcp2 loss-of-function mutant indicated that NbXCPs control the timing of de novo TE formation at the graft junction. Moreover, grafts of the NbXCP1 overexpressor increased the scion growth rate as well as the fruit size. Thus, we identified gene modules for TE formation at the graft boundary and demonstrated potential ways to enhance Nicotiana interfamily grafting.

Genomic characterization of the NAC transcription factors, directed at understanding their functions involved in endocarp lignification of iron walnut (Juglans sigillata Dode)

The NAC (NAM, ATAF1/2, and CUC2) transcription factors (TF), one of the largest plant-specific gene families, play important roles in the regulation of plant growth and development, stress response and disease resistance. In particular, several NAC TFs have been identified as master regulators of secondary cell wall (SCW) biosynthesis. Iron walnut (Juglans sigillata Dode), an economically important nut and oilseed tree, has been widely planted in the southwest China. The thick and high lignified shell derived endocarp tissues, however, brings troubles in processing processes of products in industry. It is indispensable to dissect the molecular mechanism of thick endocarp formation for further genetic improvement of iron walnut. In the present study, based on genome reference of iron walnut, 117 NAC genes, in total, were identified and characterized in silico, which involves only computational analysis to provide insight into gene function and regulation. We found that the amino acids encoded by these NAC genes varied from 103 to 1,264 in length, and conserved motif numbers ranged from 2 to 10. The JsiNAC genes were unevenly distributed across the genome of 16 chromosomes, and 96 of these genes were identified as segmental duplication genes. Furthermore, 117 JsiNAC genes were divided into 14 subfamilies (A-N) according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana and common walnut (Juglans regia). Furthermore, tissue-specific expression pattern analysis demonstrated that a majority of NAC genes were constitutively expressed in five different tissues (bud, root, fruit, endocarp, and stem xylem), while a total of 19 genes were specifically expressed in endocarp, and most of them also showed high and specific expression levels in the middle and late stages during iron walnut endocarp development. Our result provided a new insight into the gene structure and function of JsiNACs in iron walnut, and identified key candidate JsiNAC genes involved in endocarp development, probably providing mechanistic insight into shell thickness formation across nut species.

Lysine 27 of histone H3.3 is a fine modulator of developmental gene expression and stands as an epigenetic checkpoint for lignin biosynthesis in Arabidopsis

Chromatin is a dynamic platform within which gene expression is controlled by epigenetic modifications, notably targeting amino acid residues of histone H3. Among them is lysine 27 of H3 (H3K27), the trimethylation of which by the Polycomb Repressive Complex 2 (PRC2) is instrumental in regulating spatiotemporal patterns of key developmental genes. H3K27 is also subjected to acetylation and is found at sites of active transcription. Most information on the function of histone residues and their associated modifications in plants was obtained from studies of loss-of-function mutants for the complexes that modify them. To decrypt the genuine function of H3K27, we expressed a non-modifiable variant of H3 at residue K27 (H3.3^K27A ) in Arabidopsis, and developed a multi-scale approach combining in-depth phenotypical and cytological analyses, with transcriptomics and metabolomics. We uncovered that the H3.3^K27A variant causes severe developmental defects, part of them are reminiscent of PRC2 mutants, part of them are new. They include early flowering, increased callus formation and short stems with thicker xylem cell layer. This latest phenotype correlates with mis-regulation of phenylpropanoid biosynthesis. Overall, our results reveal novel roles of H3K27 in plant cell fates and metabolic pathways, and highlight an epigenetic control point for elongation and lignin composition of the stem.

Ectopic Expression of FvVND4c Promotes Secondary Cell Wall Thickening and Flavonoid Accumulation in Fragaria vesca

Secondary cell wall (SCW) thickening has a significant effect on the growth and development of plants, as well as in the resistance to various biotic and abiotic stresses. Lignin accounts for the strength of SCW. It is synthesized through the phenylpropanoid pathway that also leads to flavonoid synthesis. The coupling strategies for lignin and flavonoid syntheses are diverse in plants. How their syntheses are balanced by transcriptional regulation in fleshy fruits is still unclear. The diploid strawberry (Fragaria vesca) is a model for fleshy fruits research due to its small genome and wide scope of genetic transformation. SCW thickening is regulated by a multilevel transcriptional regulatory network wherein vascular-related NAC domains (VNDs) act as key regulators. In this study, we systematically characterized VNDs in Fragaria vesca and explored their functions. The overexpression of FvVND4c in diploid strawberry fruits resulted in SCW thickening and fruit color changes accompanied with the accumulation of lignin and flavonoids. Genes related to these phenotypes were also induced upon FvVND4c overexpression. Among the induced genes, we found FvMYB46 to be a direct downstream regulator of FvVND4c. The overexpression of FvMYB46 resulted in similar phenotypes as FvVND4c, except for the color change. Transcriptomic analyses suggest that both FvVND4c and FvMYB46 act on phenylpropanoid and flavonoid biosynthesis pathways, and induce lignin synthesis for SCW. These results suggest that FvVND4c and FvMYB46 cooperatively regulate SCW thickening and flavonoid accumulation in Fragaria vesca.

The RLCK-VND6 module coordinates secondary cell wall formation and adaptive growth in rice

The orderly deposition of secondary cell wall (SCW) in plants is implicated in various biological programs and is precisely controlled. Although many positive and negative regulators of SCW have been documented, the molecular mechanisms underlying SCW formation coordinated with distinct cellular physiological processes during plant adaptive growth remain largely unclear. Here, we report the identification of Cellulose Synthase co-expressed Kinase1 (CSK1), which encodes a receptor-like cytoplasmic kinase, as a negative regulator of SCW formation and its signaling cascade in rice. Transcriptome deep sequencing of developing internodes and genome-wide co-expression assays revealed that CSK1 is co-expressed with cellulose synthase genes and is responsive to various stress stimuli. The increased SCW thickness and vigorous vessel transport in csk1 indicate that CSK1 functions as a negative regulator of SCW biosynthesis. Through observation of green fluorescent protein-tagged CSK1 in rice protoplasts and stable transgenic plants, we found that CSK1 is localized in the nucleus and cytoplasm adjacent to the plasma membrane. Biochemical and molecular assays demonstrated that CSK1 phosphorylates VASCULAR-RELATED NAC-DOMAIN 6 (VND6), a master SCW-associated transcription factor, in the nucleus, which reduces the transcription of a suite of SCW-related genes, thereby attenuating SCW accumulation. Consistently, genetic analyses show that CSK1 functions upstream of VND6 in regulating SCW formation. Interestingly, our physiological analyses revealed that CSK1 and VND6 are involved in abscisic acid-mediated regulation of cell growth and SCW deposition. Taken together, these results indicate that the CSK1-VND6 module is an important component of the SCW biosynthesis machinery, which coordinates SCW accumulation and adaptive growth in rice. Our study not only identifies a new regulator of SCW biosynthesis but also reveals a fine-tuned mechanism for precise control of SCW deposition, offering tools for rationally tailoring agronomic traits.

HD-ZIP III-dependent local promotion of brassinosteroid synthesis suppresses vascular cell division in Arabidopsis root apical meristem

Spatiotemporal control of cell division in the meristem is vital for plant growth. In the stele of the root apical meristem (RAM), procambial cells divide periclinally to increase the number of vascular cell files. Class III homeodomain leucine zipper (HD-ZIP III) proteins are key transcriptional regulators of RAM development and suppress the periclinal division of vascular cells in the stele; however, the mechanism underlying the regulation of vascular cell division by HD-ZIP III transcription factors (TFs) remains largely unknown. Here, we performed transcriptome analysis to identify downstream genes of HD-ZIP III and found that HD-ZIP III TFs positively regulate brassinosteroid biosynthesis-related genes, such as CONSTITUTIVE PHOTOMORPHOGENIC DWARF (CPD), in vascular cells. Introduction of pREVOLUTA::CPD in a quadruple loss-of-function mutant of HD-ZIP III genes partly rescued the phenotype in terms of the vascular defect in the RAM. Treatment of a quadruple loss-of-function mutant, a gain-of-function mutant of HD-ZIP III, and the wild type with brassinosteroid and a brassinosteroid synthesis inhibitor also indicated that HD-ZIP III TFs act together to suppress vascular cell division by increasing brassinosteroid levels. Furthermore, brassinosteroid application suppressed the cytokinin response in vascular cells. Together, our findings suggest that the suppression of vascular cell division by HD-ZIP III TFs is caused, at least in part, by the increase in brassinosteroid levels through the transcriptional activation of brassinosteroid biosynthesis genes in the vascular cells of the RAM. This elevated brassinosteroid level suppresses cytokinin response in vascular cells, inhibiting vascular cell division in the RAM.

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.

Functional identification of moss PpGATA1 provides insights into the evolution of LLM-domain B-GATA transcription factors in plants

B-GATA transcription factors with the LLM domain (LLM-domain B-GATAs) play important roles in developmental processes and environmental responses in flowering plants. Their characterization can therefore provide insights into the structural and functional evolution of functional gene families. Phylogenetic and sequence analysis suggests that LLM-domain B-GATAs evolved from ancestral GATA transcription factors before the divergence of chlorophyte algae and Streptophyta. We compared the function of PpGATA1, a LLM-domain B-GATA gene in moss Physcomitrium patens, with Arabidopsis thaliana counterparts and showed that, in P. patens, PpGATA1 controls growth and greening in haploid gametophytes, while in transgenic Arabidopsis it affects germination, leaf development, flowering time, greening and light responses in diploid sporophytes. These PpGATA1 functions are similar to those of Arabidopsis counterparts, AtGNC, AtGNL and AtGATA17. PpGATA1 was able to complement the role of GNC and GNL in a gnc gnl double mutant, and the LLM domains of PpGATA1 and GNC behaved similarly. The functions of LLM-domain B-GATAs regulating hypocotyl elongation and cotyledon epinasty in flowering plants pre-exist before the divergence of mosses and the lineage leading to flowering plants. This study sheds light on adaption of PpGATA1 and its homologs to new developmental designs during the evolution.

Transcriptome analysis in contrasting maize inbred lines and functional analysis of five maize NAC genes under drought stress treatment

Drought substantially influences crop growth and development. NAC (NAM, ATAF1/2, and CUC2) transcription factors (TFs) have received much attention for their critical roles in drought stress responses. To explore the maize NAC genes in response to drought stress, the transcriptome sequencing data of NAC TFs in two maize inbred lines, the drought tolerance line H082183 and the sensitive line Lv28, were used to screen the differentially expressed genes (DEGs). There were 129 maize NAC protein-coding genes identified, of which 15 and 20 NAC genes were differentially expressed between the two genotypes under MD and SD treatments, respectively. Meanwhile, the phylogenetic relationship of 152 non-redundant NAC family TFs in maize was generated. The maize NAC family proteins were grouped into 13 distinct subfamilies. Five drought stress-responsive NAC family members, which were designed as ZmNAP, ZmNAC19, ZmNAC4, ZmJUB1(JUBGBRUNNEN1), and ZmNAC87, were selected for further study. The expression of ZmNAP, ZmNAC19, ZmNAC4, ZmJUB1, and ZmNAC87 were significantly induced by drought, dehydration, polyethylene glycol (PEG) stress, and abscisic acid (ABA) treatments. The overexpressing Arabidopsis of these five NAC genes was generated for functional characterization, respectively. Under different concentrations of NaCl, D-mannitol stress, and ABA treatments, the sensitivity of ZmNAP-, ZmNAC19-, ZmNAC4-, ZmJUB1-, and ZmNAC87-overexpressing lines was significantly increased at the germination stage compared to the wild-type lines. The overexpression of these five NAC members significantly improved the drought stress tolerance in transgenic Arabidopsis. Yeast two-hybrid screening analysis revealed that ZmNAP may cooperatively interact with 11 proteins including ZmNAC19 to activate the drought stress response. The above results inferred that ZmNAP, ZmNAC19, ZmNAC4, ZmJUB1, and ZmNAC87 may play important roles in the plant response to drought stress and may be useful in bioengineering breeding and drought tolerance improvement.

Modulation of lignin biosynthesis for drought tolerance in plants

Lignin is a complex polymer that is embedded in plant cell walls to provide physical support and water protection. For these reasons, the production of lignin is closely linked with plant adaptation to terrestrial regions. In response to developmental cues and external environmental conditions, plants use an elaborate regulatory network to determine the timing and location of lignin biosynthesis. In this review, we summarize the canonical lignin biosynthetic pathway and transcriptional regulatory network of lignin biosynthesis, consisting of NAC and MYB transcription factors, to explain how plants regulate lignin deposition under drought stress. Moreover, we discuss how the transcriptional network can be applied to the development of drought tolerant plants.

Aspartic proteases modulate programmed cell death and secondary cell wall synthesis during wood formation in poplar

Programmed cell death (PCD) is essential for wood development in trees. However, the determination of crucial factors involved in xylem PCD of wood development is still lacking. Here, two Populus trichocarpa typical aspartic protease (AP) genes, AP17 and AP45, modulate xylem maturation, especially fibre PCD, during wood formation. AP17 and AP45 were dominantly expressed in the fibres of secondary xylem, as suggested by GUS expression in APpro::GUS transgenic plants. Cas9/gRNA-induced AP17 or AP45 mutants delayed secondary xylem fibre PCD, and ap17ap45 double mutants showed more serious defects. Conversely, AP17 overexpression caused premature PCD in secondary xylem fibres, indicating a positive modulation in wood fibre PCD. Loss of AP17 and AP45 did not alter wood fibre wall thickness, whereas the ap17ap45 mutants showed a low lignin content in wood. However, AP17 overexpression led to a significant decrease in wood fibre wall thickness and lignin content, revealing the involvement in secondary cell wall synthesis during wood formation. In addition, the ap17ap45 mutant and AP17 overexpression plants resulted in a significant increase in saccharification yield in wood. Overall, AP17 and AP45 are crucial modulators in xylem maturation during wood development, providing potential candidate genes for engineering lignocellulosic wood for biofuel utilization.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages)

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

MdFLP enhances drought tolerance by regulating MdNAC019 in self-rooted apple stocks

Self-rooted apple stocks are widely used for the production of apples worldwide. However, self-rooted apple stocks are weak due to shallow roots and poor grounding, resulting in poor drought resistance. Therefore, it is essential to understand the molecular mechanisms to develop self-rooted apple stock cultivars with drought resistance. We reported that MdFLP, an R2R3-MYB transcription factor, directly binds to the promoter of MdNAC019, activating its transcription and consequently enhancing drought tolerance in self-rooted apple stocks. In addition, MdFLP indirectly activates the transcriptional expression of abiotic stress-related genes, namely, MdERF6 and MdZAT10. The plants overexpressing MdFLP displayed stronger drought tolerance, whereas MdFLP-RNAi plants showed weak drought tolerance compared with non-transgenic "Gala" plants, indicating that MdFLP regulates drought tolerance in self-rooted apple stocks. Altogether, we believe that our findings provide novel insights into the functions of MdFLP in the regulation of drought tolerance in self-rooted apple stocks.

Current Understanding of the Genetics and Molecular Mechanisms Regulating Wood Formation in Plants

Unlike herbaceous plants, woody plants undergo volumetric growth (a.k.a. secondary growth) through wood formation, during which the secondary xylem (i.e., wood) differentiates from the vascular cambium. Wood is the most abundant biomass on Earth and, by absorbing atmospheric carbon dioxide, functions as one of the largest carbon sinks. As a sustainable and eco-friendly energy source, lignocellulosic biomass can help address environmental pollution and the global climate crisis. Studies of Arabidopsis and poplar as model plants using various emerging research tools show that the formation and proliferation of the vascular cambium and the differentiation of xylem cells require the modulation of multiple signals, including plant hormones, transcription factors, and signaling peptides. In this review, we summarize the latest knowledge on the molecular mechanism of wood formation, one of the most important biological processes on Earth.

Secondary cell wall patterning-connecting the dots, pits and helices

All plant cells are encased in primary cell walls that determine plant morphology, but also protect the cells against the environment. Certain cells also produce a secondary wall that supports mechanically demanding processes, such as maintaining plant body stature and water transport inside plants. Both these walls are primarily composed of polysaccharides that are arranged in certain patterns to support cell functions. A key requisite for patterned cell walls is the arrangement of cortical microtubules that may direct the delivery of wall polymers and/or cell wall producing enzymes to certain plasma membrane locations. Microtubules also steer the synthesis of cellulose-the load-bearing structure in cell walls-at the plasma membrane. The organization and behaviour of the microtubule array are thus of fundamental importance to cell wall patterns. These aspects are controlled by the coordinated effort of small GTPases that probably coordinate a Turing's reaction-diffusion mechanism to drive microtubule patterns. Here, we give an overview on how wall patterns form in the water-transporting xylem vessels of plants. We discuss systems that have been used to dissect mechanisms that underpin the xylem wall patterns, emphasizing the VND6 and VND7 inducible systems, and outline challenges that lay ahead in this field.

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.

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.

Enhancement of Secondary Cell Wall Formation in Poplar Xylem Using a Self-Reinforced System of Secondary Cell Wall-Related Transcription Factors

The secondary cell wall (SCW) in the xylem is one of the largest sink organs of carbon in woody plants, and is considered a promising sustainable bioresource for biofuels and biomaterials. To enhance SCW formation in poplar (Populus sp.) xylem, we developed a self-reinforced system of SCW-related transcription factors from Arabidopsis thaliana, involving VASCULAR-RELATED NAC-DOMAIN7 (VND7), SECONDARY WALL-ASSOCIATED NAC-DOMAIN PROTEIN 1/NAC SECONDARY WALL THICKENING-PROMOTING FACTOR3 (SND1/NST3), and MYB46. In this system, these transcription factors were fused with the transactivation domain VP16 and expressed under the control of the Populus trichocarpa CesA18 (PtCesA18) gene promoter, creating the chimeric genes PtCesA18pro::AtVND7:VP16, PtCesA18pro::AtSND1:VP16, and PtCesA18pro::AtMYB46:VP16. The PtCesA18 promoter is active in tissues generating SCWs, and can be regulated by AtVND7, AtSND1, and AtMYB46; thus, the expression levels of PtCesA18pro::AtVND7:VP16, PtCesA18pro::AtSND1:VP16, and PtCesA18pro::AtMYB46:VP16 are expected to be boosted in SCW-generating tissues. In the transgenic hybrid aspens (Populus tremula × tremuloides T89) expressing PtCesA18pro::AtSND1:VP16 or PtCesA18pro::AtMYB46:VP16 grown in sterile half-strength Murashige and Skoog growth medium, SCW thickening was significantly enhanced in the secondary xylem cells, while the PtCesA18pro::AtVND7:VP16 plants showed stunted xylem formation, possibly because of the enhanced programmed cell death (PCD) in the xylem regions. After acclimation, the transgenic plants were transferred from the sterile growth medium to pots of soil in the greenhouse, where only the PtCesA18pro::AtMYB46:VP16 aspens survived. A nuclear magnetic resonance footprinting cell wall analysis and enzymatic saccharification analysis demonstrated that PtCesA18pro::AtMYB46:VP16 influences cell wall properties such as the ratio of syringyl (S) and guaiacyl (G) units of lignin, the abundance of the lignin β-aryl ether and resinol bonds, and hemicellulose acetylation levels. Together, these data indicate that we have created a self-reinforced system using SCW-related transcription factors to enhance SCW accumulation.

Functional Analysis of Poplar Sombrero-Type NAC Transcription Factors Yields a Strategy to Modify Woody Cell Wall Properties

Woody cells generate lignocellulosic biomass, which is a promising sustainable bioresource for wide industrial applications. Woody cell differentiation in vascular plants, including the model plant poplar (Populus trichocarpa), is regulated by a set of NAC family transcription factors, the VASCULAR-RELATED NAC-DOMAIN (VND), NAC SECONDARY CELL WALL THICKENING PROMOTING FACTOR (NST)/SND, and SOMBRERO (SMB) (VNS)-related proteins, but the precise contributions of each VNS protein to wood quality are unknown. Here, we performed a detailed functional analysis of the poplar SMB-type VNS proteins PtVNS13-PtVNS16. PtVNS13-PtVNS16 were preferentially expressed in the roots of young poplar plantlets, similar to the Arabidopsis thalianaSMB gene. PtVNS13 and PtVNS14, as well as the NST-type PtVNS11, suppressed the abnormal root cap phenotype of the Arabidopsis sombrero-3 mutant, whereas the VND-type PtVNS07 gene did not, suggesting a functional gap between SMB- or NST-type VNS proteins and VND-type VNS proteins. Overexpressing PtVNS13-PtVNS16 in Arabidopsis seedlings and poplar leaves induced ectopic xylem-vessel-like cells with secondary wall deposition, and a transient expression assay showed that PtVNS13-16 transactivated woody-cell-related genes. Interestingly, although any VNS protein rescued the pendant stem phenotype of the Arabidopsis nst1-1 nst3-1 mutant, the resulting inflorescence stems exhibited distinct cell wall properties: poplar VNS genes generated woody cell walls with higher enzymatic saccharification efficiencies compared with Arabidopsis VNS genes. Together, our data reveal clear functional diversity among VNS proteins in woody cell differentiation and demonstrate a novel VNS-based strategy for modifying woody cell wall properties toward enhanced utilization of woody biomass.

Plant secondary cell wall proteome analysis with an inducible system for xylem vessel cell differentiation

Plant cell walls are typically composed of polysaccharide polymers and cell wall proteins (CWPs). CWPs account for approximately 10% of the plant cell wall structure and perform a wide range of functions. Previous studies have identified approximately 1000 CWPs in the model plant Arabidopsis thaliana; however, the analyses mainly targeted primary cell walls, which are generated at cell division. In contrast, little is known about CWPs in secondary cell walls (SCWs), which are rigid and contain the phenolic polymer lignin. Here, we performed a cell wall proteome analysis to obtain novel insights into CWPs in SCWs. To this end, we tested multiple methods for cell wall extraction with cultured Arabidopsis cells carrying the VND7-VP16-GR system, with which cells can be transdifferentiated into xylem-vessel-like cells with lignified SCWs by dexamethasone treatment. We then subjected the protein samples to in-gel trypsin digestion followed by LC-MS/MS analysis. The different extraction methods resulted in the detection of different cell wall fraction proteins (CWFPs). In particular, centrifugation conditions had a strong impact on the extracted CWFP species, resulting in the increased number of identified CWFPs. We successfully identified 896 proteins as CWFPs in total, including proteases, expansins, purple phosphatase, well-known lignin-related enzymes (laccase and peroxidase), and 683 of 896 proteins were newly identified CWFPs. These results demonstrate the usefulness of our CWP analysis method. Further analyses of SCW-related CWPs could be expected to produce information useful for understanding the roles of CWPs in plant cell functions.

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.

Bph30 confers resistance to brown planthopper by fortifying sclerenchyma in rice leaf sheaths

Phloem-feeding insects cause massive losses in agriculture and horticulture. Host plant resistance to phloem-feeding insects is often mediated by changes in phloem composition, which deter insect settling and feeding and decrease viability. Here, we report that rice plant resistance to the phloem-feeding brown planthopper (BPH) is associated with fortification of the sclerenchyma tissue, which is located just beneath the epidermis and a cell layer or two away from the vascular bundle in the rice leaf sheath. We found that BPHs prefer to feed on the smooth and soft region on the surface of rice leaf sheaths called the long-cell block. We identified Bph30 as a rice BPH resistance gene that prevents BPH stylets from reaching the phloem due to the fortified sclerenchyma. Bph30 is strongly expressed in sclerenchyma cells and enhances cellulose and hemicellulose synthesis, making the cell walls stiffer and sclerenchyma thicker. The structurally fortified sclerenchyma is a formidable barrier preventing BPH stylets from penetrating the leaf sheath tissues and arriving at the phloem to feed. Bph30 belongs to a novel gene family, encoding a protein with two leucine-rich domains. Another member of the family, Bph40, also conferred resistance to BPH. Collectively, the fortified sclerenchyma-mediated resistance mechanism revealed in this study expands our understanding of plant-insect interactions and opens a new path for controlling planthoppers in rice.

Calcium signaling contributes to xylem vessel cell differentiation via post-transcriptional regulation of VND7 downstream events

Secondary cell walls (SCWs) accumulate in specific cell types of vascular plants, notably xylem vessel cells. Previous work has shown that calcium ions (Ca²⁺) participate in xylem vessel cell differentiation, but whether they function in SCW deposition remains unclear. In this study, we examined the role of Ca²⁺ in SCW deposition during xylem vessel cell differentiation using Arabidopsis thaliana suspension-cultured cells carrying the VND7-inducible system, in which VND7 activity can be post-translationally upregulated to induce transdifferentiation into protoxylem-type vessel cells. We observed that extracellular Ca²⁺ concentration was a crucial determinant of differentiation, although it did not have consistent effects on the transcription of VND7-downstream genes as a whole. Increasing the Ca²⁺ concentration reduced differentiation but the cells could generate the spiral patterning of SCWs. Exposure to a calcium-channel inhibitor partly restored differentiation but resulted in abnormal branched and net-like SCW patterning. These data suggest that Ca²⁺ signaling participates in xylem vessel cell differentiation via post-transcriptional regulation of VND7-downstream events, such as patterning of SCW deposition.

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.

Gene co-expression network analysis identifies BEH3 as a stabilizer of secondary vascular development in Arabidopsis

In plants, vascular stem cells located in the cambium continuously undergo self-renewal and differentiation during secondary growth. Recent advancements in cell sorting techniques have enabled access to the transcriptional regulatory framework of cambial cells. However, mechanisms underlying the robust control of vascular stem cells remain unclear. Here, we identified a new cambium-related regulatory module through co-expression network analysis using multiple transcriptome datasets obtained from an ectopic vascular cell transdifferentiation system using Arabidopsis cotyledons, Vascular cell Induction culture System Using Arabidopsis Leaves (VISUAL). The cambium gene list included a gene encoding the transcription factor BES1/BZR1 Homolog 3 (BEH3), whose homolog BES1 negatively affects vascular stem cell maintenance. Interestingly, null beh3 mutant alleles showed a large variation in their vascular size, indicating that BEH3 functions as a stabilizer of vascular stem cells. Genetic analysis revealed that BEH3 and BES1 perform opposite functions in the regulation of vascular stem cells and the differentiation of vascular cells in the context of the VISUAL system. At the biochemical level, BEH3 showed weak transcriptional repressor activity and functioned antagonistically to other BES/BZR members by competing for binding to the brassinosteroid response element. Furthermore, mathematical modeling suggested that the competitive relationship between BES/BZR homologs leads to the robust regulation of vascular stem cells.

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.

RNA-Seq analysis reveals potential regulators of programmed cell death and leaf remodelling in lace plant (Aponogeton madagascariensis)

BACKGROUND

The lace plant (Aponogeton madagascariensis) is an aquatic monocot that develops leaves with uniquely formed perforations through the use of a developmentally regulated process called programmed cell death (PCD). The process of perforation formation in lace plant leaves is subdivided into several developmental stages: pre-perforation, window, perforation formation, perforation expansion and mature. The first three emerging "imperforate leaves" do not form perforations, while all subsequent leaves form perforations via developmentally regulated PCD. PCD is active in cells called "PCD cells" that do not retain the antioxidant anthocyanin in spaces called areoles framed by the leaf veins of window stage leaves. Cells near the veins called "NPCD cells" retain a red pigmentation from anthocyanin and do not undergo PCD. While the cellular changes that occur during PCD are well studied, the gene expression patterns underlying these changes and driving PCD during leaf morphogenesis are mostly unknown. We sought to characterize differentially expressed genes (DEGs) that mediate lace plant leaf remodelling and PCD. This was achieved performing gene expression analysis using transcriptomics and comparing DEGs among different stages of leaf development, and between NPCD and PCD cells isolated by laser capture microdissection.

RESULTS

Transcriptomes were sequenced from imperforate, pre-perforation, window, and mature leaf stages, as well as PCD and NPCD cells isolated from window stage leaves. Differential expression analysis of the data revealed distinct gene expression profiles: pre-perforation and window stage leaves were characterized by higher expression of genes involved in anthocyanin biosynthesis, plant proteases, expansins, and autophagy-related genes. Mature and imperforate leaves upregulated genes associated with chlorophyll development, photosynthesis, and negative regulators of PCD. PCD cells were found to have a higher expression of genes involved with ethylene biosynthesis, brassinosteroid biosynthesis, and hydrolase activity whereas NPCD cells possessed higher expression of auxin transport, auxin signalling, aspartyl proteases, cysteine protease, Bag5, and anthocyanin biosynthesis enzymes.

CONCLUSIONS

RNA sequencing was used to generate a de novo transcriptome for A. madagascariensis leaves and revealed numerous DEGs potentially involved in PCD and leaf remodelling. The data generated from this investigation will be useful for future experiments on lace plant leaf development and PCD in planta.

Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem

Silver birch (Betula pendula Roth) is an economically important species in Northern Europe. The current research focused on the molecular background of different xylogenesis scenarios in the birch trunks. The study objects were two forms of silver birch, silver birch trees, and Karelian birch trees; the latter form is characterized by the formation of two types of wood, non-figured (straight-grained) and figured, respectively, while it is currently not clear which factors cause this difference. We identified VND/NST/SND genes that regulate secondary cell wall biosynthesis in the birch genome and revealed differences in their expression in association with the formation of xylem with different ratios of structural elements. High expression levels of BpVND7 accompanied differentiation of the type of xylem which is characteristic of the species. At the same time, the appearance of figured wood was accompanied by the low expression levels of the VND genes and increased levels of expression of NST and SND genes. We identified BpARF5 as a crucial regulator of auxin-dependent vascular patterning and its direct target—BpHB8. A decrease in the BpARF5 level expression in differentiating xylem was a specific characteristic of both Karelian birch with figured and non-figured wood. Decreased BpARF5 level expression in non-figured trees accompanied by decreased BpHB8 and VND/NST/SND expression levels compared to figured Karelian birch trees. According to the results obtained, we suggested silver birch forms differing in wood anatomy as valuable objects in studying the regulation of xylogenesis.

Transgene CpNAC68 from Wintersweet (Chimonanthus praecox) Improves Arabidopsis Survival of Multiple Abiotic Stresses

The NAC (NAM, ATAFs, CUC) family of transcription factors (TFs) play a pivotal role in regulating all processes of the growth and development of plants, as well as responses to biotic and abiotic stresses. Yet, the functions of NACs from non-model plant species remains largely uncharacterized. Here, we characterized the stress-responsive effects of a NAC gene isolated from wintersweet, an ornamental woody plant that blooms in winter when temperatures are low. CpNAC68 is clustered in the NAM subfamily. Subcellular localization and transcriptional activity assays demonstrated a nuclear protein that has transcription activator activities. qRT-PCR analyses revealed that CpNAC68 was ubiquitously expressed in old flowers and leaves. Additionally, the expression of CpNAC68 is induced by disparate abiotic stresses and hormone treatments, including drought, heat, cold, salinity, GA, JA, and SA. Ectopic overexpression of CpNAC68 in Arabidopsis thaliana enhanced the tolerance of transgenic plants to cold, heat, salinity, and osmotic stress, yet had no effect on growth and development. The survival rate and chlorophyll amounts following stress treatments were significantly higher than wild type Arabidopsis, and were accompanied by lower electrolyte leakage and malondialdehyde (MDA) amounts. In conclusion, our study demonstrates that CpNAC68 can be used as a tool to enhance plant tolerance to multiple stresses, suggesting a role in abiotic stress tolerance in wintersweet.

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.

BEL1-like Homeodomain Protein BLH6a Is a Negative Regulator of CAl5H2 in Sinapyl Alcohol Monolignol Biosynthesis in Poplar

Lignin is one of the major components of xylem cell walls in tree stems. The lignin in the wood of most flowering plants (dicotyledonous angiosperms) is typically polymerized from three monolignol precursors, coniferyl alcohol, sinapyl alcohol, and p-coumaroyl alcohol, resulting in guaiacyl (G), syringyl (S), and hydroxyphenyl (H) subunits, respectively. In this study, we focus on the transcriptional regulation of a coniferaldehyde 5-hydroxylase (CAld5H2) gene, which encodes a key enzyme for sinapyl alcohol biosynthesis. We carried out a yeast one-hybrid (Y1H) screen to identify candidate upstream transcription factors (TFs) regulating CAld5H2. We obtained 12 upstream TFs as potential regulators of CAld5H2. One of these TF genes, BLH6a, encodes a BEL1-like homeodomain (BLH) protein and negatively regulated the CAld5H2 promoter activity. The direct regulation of CAld5H2 promoter by BLH6a was supported by chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) and dominant repression of BLH6a in transgenic plants. Luciferase complementation imaging analyses showed extensive protein-protein interactions among these 12 TFs. We propose that BLH6a is a negative regulator of CAld5H2, which acts through combinatorial regulation of multiple TFs for sinapyl alcohol (S monolignol) biosynthesis in poplar.

The BOP-type co-transcriptional regulator NODULE ROOT1 promotes stem secondary growth of the tropical Cannabaceae tree Parasponia andersonii

Tree stems undergo a massive secondary growth in which secondary xylem and phloem tissues arise from the vascular cambium. Vascular cambium activity is driven by endogenous developmental signalling cues and environmental stimuli. Current knowledge regarding the genetic regulation of cambium activity and secondary growth is still far from complete. The tropical Cannabaceae tree Parasponia andersonii is a non-legume research model of nitrogen-fixing root nodulation. Parasponia andersonii can be transformed efficiently, making it amenable for CRISPR-Cas9-mediated reverse genetics. We considered whether P. andersonii also could be used as a complementary research system to investigate tree-related traits, including secondary growth. We established a developmental map of stem secondary growth in P. andersonii plantlets. Subsequently, we showed that the expression of the co-transcriptional regulator PanNODULE ROOT1 (PanNOOT1) is essential for controlling this process. PanNOOT1 is orthologous to Arabidopsis thaliana BLADE-ON-PETIOLE1 (AtBOP1) and AtBOP2, which are involved in the meristem-to-organ-boundary maintenance. Moreover, in species forming nitrogen-fixing root nodules, NOOT1 is known to function as a key nodule identity gene. Parasponia andersonii CRISPR-Cas9 loss-of-function Pannoot1 mutants are altered in the development of the xylem and phloem tissues without apparent disturbance of the cambium organization and size. Transcriptomic analysis showed that the expression of key secondary growth-related genes is significantly down-regulated in Pannoot1 mutants. This allows us to conclude that PanNOOT1 positively contributes to the regulation of stem secondary growth. Our work also demonstrates that P. andersonii can serve as a tree research system.

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.

Fiber Cell-Specific Expression of the VP16-Fused Ethylene Response Factor 41 Protein Increases Biomass Yield and Alters Lignin Composition

Arabidopsis thaliana transcription factors belonging to the ERFIIId and ERFIIIe subclade (ERFIIId/e) of the APETALA 2/ethylene response factor (AP2/ERF) family enhance primary cell wall (PCW) formation. These transcription factors activate expression of genes encoding PCW-type cellulose synthase (CESA) subunits and other genes for PCW biosynthesis. In this study, we show that fiber-specific expression of ERF035-VP16 and ERF041-VP16, which are VP16-fused proteins of ERFIIId/e members, promote cell wall thickening in a wild-type background with a concomitant increase of alcohol insoluble residues (cell wall content) per fresh weight (FW) and monosaccharides related to the PCW without affecting plant growth. Furthermore, in the ERF041-VP16 lines, the total amount of lignin and the syringyl (S)/guaiacyl (G) ratio decreased, and the enzymatic saccharification yield of glucose from cellulose per fresh weight improved. In these lines, PCW-type CESA genes were upregulated and ferulate 5-hydropxylase1 (F5H1), which is necessary for production of the S unit lignin, was downregulated. In addition, various changes in the expression levels of transcription factors regulating secondary cell wall (SCW) formation were observed. In conclusion, fiber cell-specific ERF041-VP16 improves biomass yield, increases PCW components, and alters lignin composition and deposition and may be suitable for use in future molecular breeding programs of biomass crops.

XND1 Regulates Secondary Wall Deposition in Xylem Vessels through the Inhibition of VND Functions

Secondary wall deposition in xylem vessels is activated by Vascular-Related NAC Domain proteins (VNDs) that belong to a group of secondary wall NAC (SWN) transcription factors. By contrast, Xylem NAC Domain1 (XND1) negatively regulates secondary wall deposition in xylem vessels when overexpressed. The mechanism by which XND1 exerts its functions remains elusive. We employed the promoter of the fiber-specific Secondary Wall-Associated NAC Domain1 (SND1) gene to ectopically express XND1 in fiber cells to investigate its mechanism of action on secondary wall deposition. Ectopic expression of XND1 in fiber cells severely diminished their secondary wall deposition and drastically reduced the expression of SWN-regulated downstream transcription factors and secondary wall biosynthetic genes but not that of the SWN genes themselves. Transactivation analyses revealed that XND1 specifically inhibited SWN-activated expression of these downstream genes but not their MYB46-activated expression. Both the NAC domain and the C-terminus of XND1 were required for its inhibitory function and its NAC domain interacted with the DNA-binding domains of SWNs. XND1 was shown to be localized in the cytoplasm and the nucleus and its co-expression with VND6 resulted in the cytoplasmic sequestration of VND6. Furthermore, the C-terminus of XND1 was indispensable for the XND1-mediated cytoplasmic retention of VND6 and its fusion to VND6 was able to direct VND6 to the cytoplasm and render it unable to activate the gene expression. Since the XND1 gene is specifically expressed in xylem cells, these results indicate that XND1 acts through inhibiting VND functions to negatively regulate secondary wall deposition in xylem vessels.

NtMYB12a acts downstream of sucrose to inhibit fatty acid accumulation by targeting lipoxygenase and SFAR genes in tobacco

MYB12 promotes flavonol biosynthesis in plants by targeting several early biosynthesis genes (EBGs) of this pathway. The transcriptions of these EBGs are also induced by sucrose signal. However, whether MYB12 is activated by sucrose signal and what the other roles MYB12 has in regulating plant metabolism are poorly understood. In this study, two NtMYB12 genes were cloned from Nicotiana tabacum. Both NtMYB12a and NtMYB12b are involved in regulating flavonoids biosynthesis in tobacco. NtMYB12a is further shown to inhibit the accumulation of fatty acid (FA) in tobacco leaves and seeds. Post-translational activation and chromatin immunoprecipitation assays demonstrate that NtMYB12a directly promotes the transcriptions of NtLOX6, NtLOX5, NtSFAR4 and NtGDSL2, which encode lipoxygenase (LOX) or SFAR enzymes catalyzing the degradation of FA. NtLOX6 and NtLOX5 are shown to prevent the accumulation of FA in the mature seeds and significantly reduced the percentage of polyunsaturated fatty acids (PUFAs) in tobacco. Sucrose stimulates the transcription of NtMYB12a, and loss function of NtMYB12a partially suppresses the decrease of FA content in tobacco seedlings caused by sucrose treatment. The regulation of sucrose on the expression of NtLOX6 and NtGDSL2 genes is mediated by NtMYB12a, whereas those of NtLOX5 and NtSFAR4 genes are independent of sucrose.

HOMEODOMAIN GLABROUS2 regulates cellulose biosynthesis in seed coat mucilage by activating CELLULOSE SYNTHASE5

Numerous proteins involved in cellulose biosynthesis and assembly have been functionally characterized. Nevertheless, we have a limited understanding of the mechanisms underlying the transcriptional regulation of the genes that encode these proteins. Here, we report that HOMEODOMAIN GLABROUS2 (HDG2), a Homeobox-Leucine Zipper IV transcription factor, regulates cellulose biosynthesis in Arabidopsis (Arabidopsis thaliana) seed coat mucilage. HDG2 is a transcriptional activator with the transactivation domain located within its Leucine-Zipper domain. Transcripts of HDG2 were detected specifically in seed coat epidermal cells with peak expression at 10 d postanthesis. Disruptions of HDG2 led to seed coat mucilage with aberrant morphology due to a reduction in its crystalline cellulose content. Electrophoretic mobility shift and yeast one-hybrid assays, together with chromatin immunoprecipitation and quantitative PCR, provided evidence that HDG2 directly activates CELLULOSE SYNTHASE5 (CESA5) expression by binding to the L1-box cis-acting element in its promoter. Overexpression of CESA5 partially rescued the mucilage defects of hdg2-3. Together, our data suggest that HDG2 directly activates CESA5 expression and thus is a positive regulator of cellulose biosynthesis in seed coat mucilage.

NAC and MYB Families and Lignin Biosynthesis-Related Members Identification and Expression Analysis in Melilotus albus

Melilotus albus is an annual or biennial legume species that adapts to extreme environments via its high stress tolerance. NAC and MYB transcription factors (TFs) are involved in the regulation of lignin biosynthesis, which has not been studied in M. albus. A total of 101 MaNAC and 299 MaMYB members were identified based on M. albus genome. Chromosome distribution and synteny analysis indicated that some genes underwent tandem duplication. Ka/Ks analysis suggested that MaNACs and MaMYBs underwent strong purifying selection. Stress-, hormone- and development-related cis-elements and MYB-binding sites were identified in the promoter regions of MaNACs and MaMYBs. Five MaNACs, two MaMYBs and ten lignin biosynthesis genes were identified as presenting coexpression relationships according to weighted gene coexpression network analysis (WGCNA). Eleven and thirteen candidate MaNAC and MaMYB genes related to lignin biosynthesis were identified, respectively, and a network comprising these genes was constructed which further confirmed the MaNAC and MaMYB relationship. These candidate genes had conserved gene structures and motifs and were highly expressed in the stems and roots, and qRT-PCR further verified the expression patterns. Overall, our results provide a reference for determining the precise role of NAC and MYB genes in M. albus and may facilitate efforts to breed low-lignin-content forage cultivars in the future.