A Conserved Core of Programmed Cell Death Indicator Genes Discriminates Developmentally and Environmentally Induced Programmed Cell Death in Plants

Graft incompatibility between pepper and tomato can be attributed to genetic incompatibility between diverged immune systems

Graft compatibility is the capacity of two plants to form cohesive vascular connections. Tomato and pepper are incompatible graft partners; however, the underlying cause of graft rejection between these two species remains unknown.We diagnosed graft incompatibility between tomato and diverse pepper varieties based on weakened biophysical stability, decreased growth, and persistent cell death using trypan blue and TUNEL assays. Transcriptomic analysis of cell death in the junction was performed using RNA-sequencing, and molecular signatures for incompatible graft response were characterized based on meta-transcriptomic comparisons with other biotic processes.We show that tomato is broadly incompatible with diverse pepper cultivars. These incompatible graft partners activate prolonged transcriptional changes that are highly enriched for defense processes. Amongst these processes was broad NLR upregulation and hypersensitive response. Using transcriptomic datasets for a variety of biotic stress treatments, we identified a significant overlap in the genetic profile of incompatible grafting and plant parasitism. In addition, we found over 1000 genes that are uniquely upregulated in incompatible grafts.Based on NLR overactivity, DNA damage, and prolonged cell death we have determined that tomato and pepper graft incompatibility is likely caused by a form of genetic incompatibility, which triggers a hyperimmune-response.

DNA methylation and lipid metabolism are involved in GA-induced maize aleurone layers PCD as revealed by transcriptome analysis

BACKGROUND

The aleurone layer is a part of many plant seeds, and during seed germination, aleurone cells undergo PCD, which is promoted by GA from the embryo. However, the numerous components of the GA signaling pathway that mediate PCD of the aleurone layers remain to be identified. Few genes and transcriptomes have been studied thus far in aleurone layers to improve our understanding of how PCD occurs and how the regulatory mechanism functions during PCD. Our previous studies have shown that histone deacetylases (HDACs) are required in GA-induced PCD of aleurone layer. To further explore the molecular mechanisms by which epigenetic modifications regulate aleurone PCD, we performed a global comparative transcriptome analysis of embryoless aleurones treated with GA or histone acetylase (HAT) inhibitors.

RESULTS

In this study, a total of 7,919 differentially expressed genes (DEGs) were analyzed, 2,554 DEGs of which were found to be common under two treatments. These identified DEGs were involved in various biological processes, including DNA methylation, lipid metabolism and ROS signaling. Further investigations revealed that inhibition of DNA methyltransferases prevented aleurone PCD, suggesting that active DNA methylation plays a role in regulating aleurone PCD. GA or HAT inhibitor induced lipoxygenase gene expression, leading to lipid degradation, but this process was not affected by DNA methylation. However, DNA methylation inhibitor could regulate ROS-related gene expression and inhibit GA-induced production of hydrogen peroxide (H2O2).

CONCLUSION

Overall, linking of lipoxygenase, DNA methylation, and H2O2 may indicate that GA-induced higher HDAC activity in aleurones causes breakdown of lipids via regulating lipoxygenase gene expression, and increased DNA methylation positively mediates H2O2 production; thus, DNA methylation and lipid metabolism pathways may represent an important and complex signaling network in maize aleurone PCD.

Genome-wide characterization of the soybean DOMAIN OF UNKNOWN FUNCTION 679 membrane protein gene family highlights their potential involvement in growth and stress response

The DMP (DUF679 membrane proteins) family is a plant-specific gene family that encodes membrane proteins. The DMP family genes are suggested to be involved in various programmed cell death processes and gamete fusion during double fertilization in Arabidopsis. However, their functional relevance in other crops remains unknown. This study identified 14 genes from the DMP family in soybean (Glycine max) and characterized their physiochemical properties, subcellular location, gene structure, and promoter regions using bioinformatics tools. Additionally, their tissue-specific and stress-responsive expressions were analyzed using publicly available transcriptome data. Phylogenetic analysis of 198 DMPs from monocots and dicots revealed six clades, with clade-I encoding senescence-related AtDMP1/2 orthologues and clade-II including pollen-specific AtDMP8/9 orthologues. The largest clade, clade-III, predominantly included monocot DMPs, while monocot- and dicot-specific DMPs were assembled in clade-IV and clade-VI, respectively. Evolutionary analysis suggests that soybean GmDMPs underwent purifying selection during evolution. Using 68 transcriptome datasets, expression profiling revealed expression in diverse tissues and distinct responses to abiotic and biotic stresses. The genes Glyma.09G237500 and Glyma.18G098300 showed pistil-abundant expression by qPCR, suggesting they could be potential targets for female organ-mediated haploid induction. Furthermore, cis-acting regulatory elements primarily related to stress-, hormone-, and light-induced pathways regulate GmDMPs, which is consistent with their divergent expression and suggests involvement in growth and stress responses. Overall, our study provides a comprehensive report on the soybean GmDMP family and a framework for further biological functional analysis of DMP genes in soybean or other crops.

RAB7 GTPases as coordinators of plant endomembrane traffic

The ras gene from rat brain (RAB) family of small GTPases is highly conserved among eukaryotes and regulates endomembrane trafficking pathways. RAB7, in particular, has been linked to various processes involved in regulating endocytic and autophagic pathways. Plants have several copies of RAB7 proteins that reflect the intricacy of their endomembrane transport systems. RAB7 activity regulates different pathways of endomembrane trafficking in plants: (1) endocytic traffic to the vacuole; (2) biosynthetic traffic to the vacuole; and (3) recycling from the late endosome to the secretory pathway. During certain developmental and stress related processes another pathway becomes activated (4) autophagic trafficking towards the vacuole that is also regulated by RAB7. RAB7s carry out these functions by interacting with various effector proteins. Current research reveals many unexplored RAB7 functions in connection with stress responses. Thus, this review describes a comprehensive summary of current knowledge of plant RAB7's functions, discusses unresolved challenges, and recommends prospective future research directions.

Industrial chicory genome gives insights into the molecular timetable of anther development and male sterility

Industrial chicory (Cichorium intybus var. sativum) is a biannual crop mostly cultivated for extraction of inulin, a fructose polymer used as a dietary fiber. F1 hybrid breeding is a promising breeding strategy in chicory but relies on stable male sterile lines to prevent self-pollination. Here, we report the assembly and annotation of a new industrial chicory reference genome. Additionally, we performed RNA-Seq on subsequent stages of flower bud development of a fertile line and two cytoplasmic male sterile (CMS) clones. Comparison of fertile and CMS flower bud transcriptomes combined with morphological microscopic analysis of anthers, provided a molecular understanding of anther development and identified key genes in a range of underlying processes, including tapetum development, sink establishment, pollen wall development and anther dehiscence. We also described the role of phytohormones in the regulation of these processes under normal fertile flower bud development. In parallel, we evaluated which processes are disturbed in CMS clones and could contribute to the male sterile phenotype. Taken together, this study provides a state-of-the-art industrial chicory reference genome, an annotated and curated candidate gene set related to anther development and male sterility as well as a detailed molecular timetable of flower bud development in fertile and CMS lines.

Functional specialization of chloroplast vesiculation (CV) duplicated genes from soybean shows partial overlapping roles during stress-induced or natural senescence

Soybean is a globally important legume crop which is highly sensitive to drought. The identification of genes of particular relevance for drought responses provides an important basis to improve tolerance to environmental stress. Chloroplast Vesiculation (CV) genes have been characterized in Arabidopsis and rice as proteins participating in a specific chloroplast-degradation vesicular pathway (CVV) during natural or stress-induced leaf senescence. Soybean genome contains two paralogous genes encoding highly similar CV proteins, CV1 and CV2. In this study, we found that expression of CV1 was differentially upregulated by drought stress in soybean contrasting genotypes exhibiting slow-wilting (tolerant) or fast-wilting (sensitive) phenotypes. CV1 reached higher induction levels in fast-wilting plants, suggesting a negative correlation between CV1 gene expression and drought tolerance. In contrast, autophagy (ATG8) and ATI-PS (ATI1) genes were induced to higher levels in slow-wilting plants, supporting a pro-survival role for these genes in soybean drought tolerance responses. The biological function of soybean CVs in chloroplast degradation was confirmed by analyzing the effect of conditional overexpression of CV2-FLAG fusions on the accumulation of specific chloroplast proteins. Functional specificity of CV1 and CV2 genes was assessed by analyzing their specific promoter activities in transgenic Arabidopsis expressing GUS reporter gene driven by CV1 or CV2 promoters. CV1 promoter responded primarily to abiotic stimuli (hyperosmolarity, salinity and oxidative stress), while the promoter of CV2 was predominantly active during natural senescence. Both promoters were highly responsive to auxin but only CV1 responded to other stress-related hormones, such as ABA, salicylic acid and methyl jasmonate. Moreover, the dark-induced expression of CV2, but not of CV1, was strongly inhibited by cytokinin, indicating similarities in the regulation of CV2 to the reported expression of Arabidopsis and rice CV genes. Finally, we report the expression of both CV1 and CV2 genes in roots of soybean and transgenic Arabidopsis, suggesting a role for the encoded proteins in root plastids. Together, the results indicate differential roles for CV1 and CV2 in development and in responses to environmental stress, and point to CV1 as a potential target for gene editing to improve crop performance under stress without compromising natural development.

Repressive ZINC FINGER OF ARABIDOPSIS THALIANA proteins promote programmed cell death in the Arabidopsis columella root cap

Developmental programmed cell death (dPCD) controls a plethora of functions in plant growth and reproduction. In the root cap of Arabidopsis (Arabidopsis thaliana), dPCD functions to control organ size in balance with the continuous stem cell activity in the root meristem. Key regulators of root cap dPCD including SOMBRERO/ANAC033 (SMB) belong to the NAC family of transcription factors. Here, we identify the C2H2 zinc finger protein ZINC FINGER OF ARABIDOPSIS THALIANA 14 ZAT14 as part of the gene regulatory network of root cap dPCD acting downstream of SMB. Similar to SMB, ZAT14-inducible misexpression leads to extensive ectopic cell death. Both the canonical EAR motif and a conserved L-box motif of ZAT14 act as transcriptional repression motifs and are required to trigger cell death. While a single zat14 mutant does not show a cell death-related phenotype, a quintuple mutant knocking out 5 related ZAT paralogs shows a delayed onset of dPCD execution in the columella and the adjacent lateral root cap. While ZAT14 is co-expressed with established dPCD-associated genes, it does not activate their expression. Our results suggest that ZAT14 acts as a transcriptional repressor controlling a so far uncharacterized subsection of the dPCD gene regulatory network active in specific root cap tissues.

The large GTPase AtGBPL3 links nuclear envelope formation and morphogenesis to transcriptional repression

Guanylate binding proteins (GBPs) are prominent regulators of immunity not known to be required for nuclear envelope formation and morphogenesis. Here we identify the Arabidopsis GBP orthologue AtGBPL3 as a lamina component with essential functions in mitotic nuclear envelope reformation, nuclear morphogenesis and transcriptional repression during interphase. AtGBPL3 is preferentially expressed in mitotically active root tips, accumulates at the nuclear envelope and interacts with centromeric chromatin as well as with lamina components transcriptionally repressing pericentromeric chromatin. Reduced expression of AtGBPL3 or associated lamina components similarly altered nuclear morphology and caused overlapping transcriptional deregulation. Investigating the dynamics of AtGBPL3-GFP and other nuclear markers during mitosis (1) revealed that AtGBPL3 accumulation on the surface of daughter nuclei precedes nuclear envelope reformation and (2) uncovered defects in this process in roots of AtGBPL3 mutants, which cause programmed cell death and impair growth. AtGBPL3 functions established by these observations are unique among dynamin-family large GTPases.

A molecular framework for grain number determination in barley

Flowering plants with indeterminate inflorescences often produce more floral structures than they require. We found that floral primordia initiations in barley (Hordeum vulgare L.) are molecularly decoupled from their maturation into grains. While initiation is dominated by flowering-time genes, floral growth is specified by light signaling, chloroplast, and vascular developmental programs orchestrated by barley CCT MOTIF FAMILY 4 (HvCMF4), which is expressed in the inflorescence vasculature. Consequently, mutations in HvCMF4 increase primordia death and pollination failure, mainly through reducing rachis greening and limiting plastidial energy supply to developing heterotrophic floral tissues. We propose that HvCMF4 is a sensory factor for light that acts in connection with the vascular-localized circadian clock to coordinate floral initiation and survival. Notably, stacking beneficial alleles for both primordia number and survival provides positive implications on grain production. Our findings provide insights into the molecular underpinnings of grain number determination in cereal crops.

Spatial and temporal regulation of parent-of-origin allelic expression in the endosperm

Genomic imprinting promotes differential expression of parental alleles in the endosperm of flowering plants and is regulated by epigenetic modification such as DNA methylation and histone tail modifications in chromatin. After fertilization, the endosperm develops through a syncytial stage before it cellularizes and becomes a nutrient source for the growing embryo. Regional compartmentalization has been shown both in early and late endosperm development, and different transcriptional domains suggest divergent spatial and temporal regional functions. The analysis of the role of parent-of-origin allelic expression in the endosperm as a whole and the investigation of domain-specific functions have been hampered by the inaccessibility of the tissue for high-throughput transcriptome analyses and contamination from surrounding tissue. Here, we used fluorescence-activated nuclear sorting (FANS) of nuclear targeted GFP fluorescent genetic markers to capture parental-specific allelic expression from different developmental stages and specific endosperm domains. This approach allowed us to successfully identify differential genomic imprinting with temporal and spatial resolution. We used a systematic approach to report temporal regulation of imprinted genes in the endosperm, as well as region-specific imprinting in endosperm domains. Analysis of our data identified loci that are spatially differentially imprinted in one domain of the endosperm, while biparentally expressed in other domains. These findings suggest that the regulation of genomic imprinting is dynamic and challenge the canonical mechanisms for genomic imprinting.

Cork Development: What Lies Within

The cork layer present in all dicotyledonous plant species with radial growth is the result of the phellogen activity, a secondary meristem that produces phellem (cork) to the outside and phelloderm inwards. These three different tissues form the periderm, an efficient protective tissue working as a barrier against external factors such as environmental aggressions and pathogen attacks. The protective function offered by cork cells is mainly due to the abundance of suberin in their cell walls. Chemically, suberin is a complex aliphatic network of long chain fatty acids and alcohols with glycerol together with aromatic units. In most woody species growing in temperate climates, the first periderm is replaced by a new functional periderm upon a few years after being formed. One exception to this bark development can be found in cork oak (Quercus suber) which display a single periderm that grows continuously. Quercus suber stands by its thick cork layer development with continuous seasonal growth. Cork raw material has been exploited by man for centuries, especially in Portugal and Spain. Nowadays, its applications have widened vastly, from the most known product, stoppers, to purses or insulating materials used in so many industries, such as construction and car production. Research on how cork develops, and the effect environmental factors on cork oak trees is extremely important to maintain production of good-quality cork, and, by maintaining cork oak stands wealthy, we are preserving a very important ecosystem both by its biodiversity and its vital social and economic role in areas already showing a population declination.

Killing me softly - Programmed cell death in plant reproduction from sporogenesis to fertilization

Regulated or programmed cell death (RCD or PCD) is a fundamental biological principle integral to a considerable variety of functions in multicellular organisms. In plants, different PCD processes are part of biotic and abiotic stress responses, but also occur as an essential aspect of unperturbed plant development. PCD is particularly abundant during plant reproduction, eliminating unwanted or no longer needed cells, tissues, or organs in a precisely controlled manner. Failure in reproductive PCD can have detrimental consequences for plant reproduction. Here we shed a light on the latest research into PCD mechanisms in plant reproduction from sex determination over sporogenesis to pollination and fertilization.

Sphingolipids are involved in insect egg-induced cell death in Arabidopsis

In Brassicaceae, hypersensitive-like programmed cell death (HR-like) is a central component of direct defenses triggered against eggs of the large white butterfly (Pieris brassicae). The signaling pathway leading to HR-like in Arabidopsis (Arabidopsis thaliana) is mainly dependent on salicylic acid (SA) accumulation, but downstream components are unclear. Here, we found that treatment with P. brassicae egg extract (EE) triggered changes in expression of sphingolipid metabolism genes in Arabidopsis and black mustard (Brassica nigra). Disruption of ceramide (Cer) synthase activity led to a significant decrease of EE-induced HR-like whereas SA signaling and reactive oxygen species levels were unchanged, suggesting that Cer are downstream activators of HR-like. Sphingolipid quantifications showed that Cer with C16:0 side chains accumulated in both plant species and this response was largely unchanged in the SA-induction deficient2 (sid2-1) mutant. Finally, we provide genetic evidence that the modification of fatty acyl chains of sphingolipids modulates HR-like. Altogether, these results show that sphingolipids play a key and specific role during insect egg-triggered HR-like.

KIL1 terminates fertility in maize by controlling silk senescence

Plant flowers have a functional life span during which pollination and fertilization occur to ensure seed and fruit development. Once flower senescence is initiated, the potential to set seed or fruit is irrevocably lost. In maize, silk strands are the elongated floral stigmas that emerge from the husk-enveloped inflorescence to intercept airborne pollen. Here we show that KIRA1-LIKE1 (KIL1), an ortholog of the Arabidopsis NAC (NAM (NO APICAL MERISTEM), ATAF1/2 (Arabidopsis thaliana Activation Factor1 and 2) and CUC (CUP-SHAPED COTYLEDON 2)) transcription factor KIRA1, promotes senescence and programmed cell death (PCD) in the silk strand base, ending the window of accessibility for fertilization of the ovary. Loss of KIL1 function extends silk receptivity and thus strongly increases kernel yield following late pollination. This phenotype offers new opportunities for possibly improving yield stability in cereal crops. Moreover, despite diverging flower morphologies and the substantial evolutionary distance between Arabidopsis and maize, our data indicate remarkably similar principles in terminating floral receptivity by PCD, whose modulation offers the potential to be widely used in agriculture.

TALE-induced cell death executors: an origin outside immunity?

Phytopathogenic bacteria inject effector proteins into plant host cells to promote disease. Plant resistance (R) genes encoding nucleotide-binding leucine-rich repeat (NLR) proteins mediate the recognition of functionally and structurally diverse microbial effectors, including transcription-activator like effectors (TALEs) from the bacterial genus Xanthomonas. TALEs bind to plant promoters and transcriptionally activate either disease-promoting host susceptibility (S) genes or cell death-inducing executor-type R genes. It is perplexing that plants contain TALE-perceiving executor-type R genes in addition to NLRs that also mediate the recognition of TALE-containing xanthomonads. We present recent findings on the evolvability of TALEs, which suggest that the native function of executors is not in plant immunity, but possibly in the regulation of developmentally controlled programmed cell death (PCD) processes.

Abiotic Stresses in Plants and Their Markers: A Practice View of Plant Stress Responses and Programmed Cell Death Mechanisms

Understanding how plants cope with stress and the intricate mechanisms thereby used to adapt and survive environmental imbalances comprise one of the most powerful tools for modern agriculture. Interdisciplinary studies suggest that knowledge in how plants perceive, transduce and respond to abiotic stresses are a meaningful way to design engineered crops since the manipulation of basic characteristics leads to physiological remodeling for plant adaption to different environments. Herein, we discussed the main pathways involved in stress-sensing, signal transduction and plant adaption, highlighting biochemical, physiological and genetic events involved in abiotic stress responses. Finally, we have proposed a list of practice markers for studying plant responses to multiple stresses, highlighting how plant molecular biology, phenotyping and genetic engineering interconnect for creating superior crops.

An ARF1-binding factor triggering programmed cell death and periderm development in pear russet fruit skin

Plants have a cuticular membrane (CM) and periderm membrane (PM), which act as barriers to terrestrial stresses. The CM covers primary organs with a continuous hydrophobic layer of waxes embedded in cutin, while the PM stacks with suberized cells outermost to the secondary tissues. The formation of native periderm is regulated by a postembryonic meristem phellogen that produces suberized phellem (cork) outwardly. However, the mechanism controlling phellogen differentiation to phellem remains to be clarified. Here, map-based cloning in a pear F1 population with segregation for periderm development in fruit skin facilitated the identification of an aspartic acid repeat deletion in Pyrus Periderm Programmed Cell Death 1.1 (PyPPCD1.1) that triggers phellogen activity for cork formation in pear russet fruit skin. PyPPCD1.1 showed preferential expression in pear fruit skin, and the encoded protein shares a structural similarity to that of the viral capsid proteins. Asp deletion in PyPPCD1.1 weakened its nuclear localization but increased its accumulation in the chloroplast. Both PyPPCD1.1 and its recessive allele directly interact with ADP-ribosylation factor 1 (ARF1). PyPPCD1.1 triggered PCD in an ARF1-dependent manner. Thus, this study identified the switch gene for PCD and periderm development and provided a new molecular regulatory mechanism underlying the development of this trait.

Oxidative stress response and programmed cell death guided by NAC013 modulate pithiness in radish taproots

Radish, Raphanus sativus L., is an important root crop that is cultivated worldwide. Owing to its evolutionary proximity to Arabidopsis thaliana, radish can be used as a model root crop in research on the molecular basis of agronomic traits. Pithiness is a significant defect that reduces the production of radish with commercial value; however, traditional breeding to eliminate this trait has thus far been unsuccessful. Here, we performed transcriptomics and genotype-by-sequencing (GBS)-based quantitative trait locus (QTL) analyses of radish inbred lines to understand the molecular basis of pithiness in radish roots. The transcriptome data indicated that pithiness likely stems from the response to oxidative stress, leading to cell death of the xylem parenchyma during the root-thickening process. Subsequently, we narrowed down a list of candidates responsible for pithiness near a major QTL and found polymorphisms in a radish homologue of Arabidopsis ANAC013 (RsNAC013), an endoplasmic reticulum bound NAC transcription factor that is targeted to the nucleus to mediate the mitochondrial retrograde signal. We analysed the effects of polymorphisms in RsNAC013 using Arabidopsis transgenic lines overexpressing RsNAC013 alleles as well as in radish inbred lines bearing these alleles. This analysis indicated that non-synonymous variations within the coding sequence result in different levels of RsNAC013 activities, thereby providing a genetic condition for root pithiness. The elevated oxidative stress or hypoxia that activates RsNAC013 for mitochondrial signalling enhances this process. Collectively, this study serves as an exemplary case of translational research taking advantage of the extensive information available from a model organism.

Initiation and Execution of Programmed Cell Death and Regulation of Reactive Oxygen Species in Plants

Programmed cell death (PCD) plays crucial roles in plant development and defence response. Reactive oxygen species (ROS) are produced during normal plant growth, and high ROS concentrations can change the antioxidant status of cells, leading to spontaneous cell death. In addition, ROS function as signalling molecules to improve plant stress tolerance, and they induce PCD under different conditions. This review describes the mechanisms underlying plant PCD, the key functions of mitochondria and chloroplasts in PCD, and the relationship between mitochondria and chloroplasts during PCD. Additionally, the review discusses the factors that regulate PCD. Most importantly, in this review, we summarise the sites of production of ROS and discuss the roles of ROS that not only trigger multiple signalling pathways leading to PCD but also participate in the execution of PCD, highlighting the importance of ROS in PCD.

Senescence-Associated Glycine max (Gm)NAC Genes: Integration of Natural and Stress-Induced Leaf Senescence

Leaf senescence is a genetically regulated developmental process that can be triggered by a variety of internal and external signals, including hormones and environmental stimuli. Among the senescence-associated genes controlling leaf senescence, the transcriptional factors (TFs) comprise a functional class that is highly active at the onset and during the progression of leaf senescence. The plant-specific NAC (NAM, ATAF, and CUC) TFs are essential for controlling leaf senescence. Several members of Arabidopsis AtNAC-SAGs are well characterized as players in elucidated regulatory networks. However, only a few soybean members of this class display well-known functions; knowledge about their regulatory circuits is still rudimentary. Here, we describe the expression profile of soybean GmNAC-SAGs upregulated by natural senescence and their functional correlation with putative AtNAC-SAGs orthologs. The mechanisms and the regulatory gene networks underlying GmNAC081- and GmNAC030-positive regulation in leaf senescence are discussed. Furthermore, new insights into the role of GmNAC065 as a negative senescence regulator are presented, demonstrating extraordinary functional conservation with the Arabidopsis counterpart. Finally, we describe a regulatory circuit which integrates a stress-induced cell death program with developmental leaf senescence via the NRP-NAC-VPE signaling module.

Transcriptional Regulation and Signaling of Developmental Programmed Cell Death in Plants

Developmental programmed cell death (dPCD) has multiple functions in plant growth and development, and is of great value for industrial production. Among them, wood formed by xylem dPCD is one of the most widely used natural materials. Therefore, it is crucial to explore the molecular mechanism of plant dPCD. The dPCD process is tightly regulated by genetic networks and is involved in the transduction of signaling molecules. Several key regulators have been identified in diverse organisms and individual PCD events. However, complex molecular networks controlling plant dPCD remain highly elusive, and the original triggers of this process are still unknown. This review summarizes the recent progress on the transcriptional regulation and signaling of dPCD during vegetative and reproductive development. It is hoped that this review will provide an overall view of the molecular regulation of dPCD in different developmental processes in plants and identify specific mechanisms for regulating these dPCD events. In addition, the application of plants in industrial production can be improved by manipulating dPCD in specific processes, such as xylogenesis.

Contrasting roles of GmNAC065 and GmNAC085 in natural senescence, plant development, multiple stresses and cell death responses

NACs are plant-specific transcription factors involved in controlling plant development, stress responses, and senescence. As senescence-associated genes (SAGs), NACs integrate age- and stress-dependent pathways that converge to programmed cell death (PCD). In Arabidopsis, NAC-SAGs belong to well-characterized regulatory networks, poorly understood in soybean. Here, we interrogated the soybean genome and provided a comprehensive analysis of senescence-associated Glycine max (Gm) NACs. To functionally examine GmNAC-SAGs, we selected GmNAC065, a putative ortholog of Arabidopsis ANAC083/VNI2 SAG, and the cell death-promoting GmNAC085, an ANAC072 SAG putative ortholog, for analyses. Expression analysis of GmNAC065 and GmNAC085 in soybean demonstrated (i) these cell death-promoting GmNACs display contrasting expression changes during age- and stress-induced senescence; (ii) they are co-expressed with functionally different gene sets involved in stress and PCD, and (iii) are differentially induced by PCD inducers. Furthermore, we demonstrated GmNAC065 expression delays senescence in Arabidopsis, a phenotype associated with enhanced oxidative performance under multiple stresses, higher chlorophyll, carotenoid and sugar contents, and lower stress-induced PCD compared to wild-type. In contrast, GmNAC085 accelerated stress-induced senescence, causing enhanced chlorophyll loss, ROS accumulation and cell death, decreased antioxidative system expression and activity. Accordingly, GmNAC065 and GmNAC085 targeted functionally contrasting sets of downstream AtSAGs, further indicating that GmNAC85 and GmNAC065 regulators function inversely in developmental and environmental PCD.

Local brassinosteroid biosynthesis enables optimal root growth

Brassinosteroid (BR) hormones are indispensable for root growth and control both cell division and cell elongation through the establishment of an increasing signalling gradient along the longitudinal root axis. Because of their limited mobility, the importance of BR distribution in achieving a signalling maximum is largely overlooked. Expression pattern analysis of all known BR biosynthetic enzymes revealed that not all cells in the Arabidopsis thaliana root possess full biosynthetic machinery, and that completion of biosynthesis relies on cell-to-cell movement of hormone precursors. We demonstrate that BR biosynthesis is largely restricted to the root elongation zone, where it overlaps with BR signalling maxima. Moreover, optimal root growth requires hormone concentrations to be low in the meristem and high in the root elongation zone, attributable to increased biosynthesis. Our finding that spatiotemporal regulation of hormone synthesis results in local hormone accumulation provides a paradigm for hormone-driven organ growth in the absence of long-distance hormone transport in plants.

Polycomb Repressive Complex 2 and KRYPTONITE regulate pathogen-induced programmed cell death in Arabidopsis

The Polycomb Repressive Complex 2 (PRC2) is well-known for its role in controlling developmental transitions by suppressing the premature expression of key developmental regulators. Previous work revealed that PRC2 also controls the onset of senescence, a form of developmental programmed cell death (PCD) in plants. Whether the induction of PCD in response to stress is similarly suppressed by the PRC2 remained largely unknown. In this study, we explored whether PCD triggered in response to immunity- and disease-promoting pathogen effectors is associated with changes in the distribution of the PRC2-mediated histone H3 lysine 27 trimethylation (H3K27me3) modification in Arabidopsis thaliana. We furthermore tested the distribution of the heterochromatic histone mark H3K9me2, which is established, to a large extent, by the H3K9 methyltransferase KRYPTONITE, and occupies chromatin regions generally not targeted by PRC2. We report that effector-induced PCD caused major changes in the distribution of both repressive epigenetic modifications and that both modifications have a regulatory role and impact on the onset of PCD during pathogen infection. Our work highlights that the transition to pathogen-induced PCD is epigenetically controlled, revealing striking similarities to developmental PCD.

Cellular oxidative stress in programmed cell death: focusing on chloroplastic 1O2 and mitochondrial cytochrome-c release

The programmed cell death (PCD) occurs when the targeted cells have fulfilled their task or under conditions as oxidative stress generated by ROS species. Thus, plants have to deal with the singlet oxygen ¹O₂ produced in chloroplasts. ¹O₂ is unlikely to act as a primary retrograde signal owing to its high reactivity and short half-life. In addition to its high toxicity, the ¹O₂ generated under an excess or low excitation energy might also act as a highly versatile signal triggering chloroplast-to-nucleus retrograde signaling (ChNRS) and nuclear reprogramming or cell death. Molecular and biochemical studies with the flu mutant, which accumulates protochlorophyllide in the dark, demonstrated that chloroplastic ¹O₂-driven EXECUTER-1 (EX1) and EX2 proteins are involved in the ¹O₂-dependent response. Both EX1 and EX2 are necessary for full suppression of ¹O₂-induced gene expression. That is, EXECUTER proteolysis via the ATP-dependent zinc protease (FtsH) is an integral part of ¹O₂-triggered retrograde signaling. The existence of at least two independent ChNRS involving EX1 and β-cyclocitral, and dihydroactinidiolide and OXI1, respectively, seem clear. Besides, this update also focuses on plant PCD and its relation with mitochondrial cytochrome-c (Cytc) release to cytosol. Changes in the dynamics and morphology of mitochondria were shown during the onset of cell death. The mitochondrial damage and translocation of Cytc may be one of the major causes of PCD triggering. Together, this current overview illustrates the complexity of the cellular response to oxidative stress development. A puzzle with the majority of its pieces still not placed.

Insights into Plant Programmed Cell Death Induced by Heavy Metals-Discovering a Terra Incognita

Programmed cell death (PCD) is a process that plays a fundamental role in plant development and responses to biotic and abiotic stresses. Knowledge of plant PCD mechanisms is still very scarce and is incomparable to the large number of studies on PCD mechanisms in animals. Quick and accurate assays, e.g., the TUNEL assay, comet assay, and analysis of caspase-like enzyme activity, enable the differentiation of PCD from necrosis. Two main types of plant PCD, developmental (dPCD) regulated by internal factors, and environmental (ePCD) induced by external stimuli, are distinguished based on the differences in the expression of the conserved PCD-inducing genes. Abiotic stress factors, including heavy metals, induce necrosis or ePCD. Heavy metals induce PCD by triggering oxidative stress via reactive oxygen species (ROS) overproduction. ROS that are mainly produced by mitochondria modulate phytotoxicity mechanisms induced by heavy metals. Complex crosstalk between ROS, hormones (ethylene), nitric oxide (NO), and calcium ions evokes PCD, with proteases with caspase-like activity executing PCD in plant cells exposed to heavy metals. This pathway leads to very similar cytological hallmarks of heavy metal induced PCD to PCD induced by other abiotic factors. The forms, hallmarks, mechanisms, and genetic regulation of plant ePCD induced by abiotic stress are reviewed here in detail, with an emphasis on plant cell culture as a suitable model for PCD studies. The similarities and differences between plant and animal PCD are also discussed.

Copy number variation of TdDof controls solid-stemmed architecture in wheat

Stem solidness is an important agronomic trait of durum (Triticum turgidum L. var. durum) and bread (Triticum aestivum L.) wheat that provides resistance to the wheat stem sawfly. This dominant trait is conferred by the SSt1 locus on chromosome 3B. However, the molecular identity and mechanisms underpinning stem solidness have not been identified. Here, we demonstrate that copy number variation of TdDof, a gene encoding a putative DNA binding with one finger protein, controls the stem solidness trait in wheat. Using map-based cloning, we localized TdDof to within a physical interval of 2.1 Mb inside the SSt1 locus. Molecular analysis revealed that hollow-stemmed wheat cultivars such as Kronos carry a single copy of TdDof, whereas solid-stemmed cultivars such as CDC Fortitude carry multiple identical copies of the gene. Deletion of all TdDof copies from CDC Fortitude resulted in the loss of stem solidness, whereas the transgenic overexpression of TdDof restored stem solidness in the TdDof deletion mutant pithless1 and conferred stem solidness in Kronos. In solid-stemmed cultivars, increased TdDof expression was correlated with the down-regulation of genes whose orthologs have been implicated in programmed cell death (PCD) in other species. Anatomical and histochemical analyses revealed that hollow-stemmed lines had stronger PCD-associated signals in the pith cells compared to solid-stemmed lines, which suggests copy number-dependent expression of TdDof could be directly or indirectly involved in the negative regulation of PCD. These findings provide opportunities to manipulate stem development in wheat and other monocots for agricultural or industrial purposes.

High-Throughput Single-Cell Transcriptome Profiling of Plant Cell Types

Single-cell transcriptome profiling of heterogeneous tissues can provide high-resolution windows into developmental dynamics and environmental responses, but its application to plants has been limited. Here, we used the high-throughput Drop-seq approach to profile >12,000 cells from Arabidopsis roots. This identified numerous distinct cell types, covering all major root tissues and developmental stages, and illuminated specific marker genes for these populations. In addition, we demonstrate the utility of this approach to study the impact of environmental conditions on developmental processes. Analysis of roots grown with or without sucrose supplementation uncovers changes in the relative frequencies of cell types in response to sucrose. Finally, we characterize the transcriptome changes that occur across endodermis development and identify nearly 800 genes with dynamic expression as this tissue differentiates. Collectively, we demonstrate that single-cell RNA-seq can be used to profile developmental processes in plants and show how they can be altered by external stimuli.

The application of single-cell transcriptome profiling to plants has been limited. Shulse et al. performed Drop-seq on Arabidopsis roots, generating a transcriptional resource for >12,000 cells across major populations. This revealed marker genes for distinct cell types, cell frequency changes resulting from sucrose addition, and genes dynamically regulated during development.

Tomato T2 ribonuclease LE is involved in the response to pathogens

T2 ribonucleases (RNases) are RNA-degrading enzymes that function in various cellular processes, mostly via RNA metabolism. T2 RNase-encoding genes have been identified in various organisms, from bacteria to mammals, and are most diverse in plants. The existence of T2 RNase genes in almost every organism suggests an important biological function that has been conserved through evolution. In plants, T2 RNases are suggested to be involved in phosphate scavenging and recycling, and are implicated in defence responses to pathogens. We investigated the function of the tomato T2 RNase LE, known to be induced by phosphate deficiency and wounding. The possible involvement of LE in pathogen responses was examined. Expression analysis showed LE induction during fungal infection and by stimuli known to be associated with pathogen inoculation, including oxalic acid and hydrogen peroxide. Analysis of LE-suppressed transgenic tomato lines revealed higher susceptibility to oxalic acid, a cell death-inducing factor, compared to the wild type. This elevated sensitivity of LE-suppressed lines was evidenced by visual signs of necrosis, and increased ion leakage and reactive oxygen species levels, indicating acceleration of cell death. Challenge of the LE-suppressed lines with the necrotrophic pathogen Botrytis cinerea resulted in accelerated development of disease symptoms compared to the wild type, associated with suppressed expression of pathogenesis-related marker genes. The results suggest a role for plant endogenous T2 RNases in antifungal activity.

Transcriptomics at Maize Embryo/Endosperm Interfaces Identifies a Transcriptionally Distinct Endosperm Subdomain Adjacent to the Embryo Scutellum

Seeds are complex biological systems comprising three genetically distinct tissues nested one inside another (embryo, endosperm, and maternal tissues). However, the complexity of the kernel makes it difficult to understand intercompartment interactions without access to spatially accurate information. Here, we took advantage of the large size of the maize (Zea mays) kernel to characterize genome-wide expression profiles of tissues at different embryo/endosperm interfaces. Our analysis identifies specific transcriptomic signatures in two interface tissues compared with whole seed compartments: the scutellar aleurone layer and the newly named endosperm adjacent to scutellum (EAS). The EAS, which appears around 9 d after pollination and persists for around 11 d, is confined to one to three endosperm cell layers adjacent to the embryonic scutellum. Its transcriptome is enriched in genes encoding transporters. The absence of the embryo in an embryo specific mutant can alter the expression pattern of EAS marker genes. The detection of cell death in some EAS cells together with an accumulation of crushed cell walls suggests that the EAS is a dynamic zone from which cell layers in contact with the embryo are regularly eliminated and to which additional endosperm cells are recruited as the embryo grows.

Overexpression of three related root-cap outermost-cell-specific C2H2-type zinc-finger protein genes suppresses the growth of Arabidopsis in an EAR-motif-dependent manner

The root meristem of Arabidopsis thaliana is protected by the root cap, the size of which is tightly regulated by the balance between the formative cell divisions and the dispersal of the outermost cells. We isolated an enhancer-tagged dominant mutant displaying the short and twisted root by the overexpression of ZINC-FINGER OF ARABIDOPSIS THALIANA1 (ZAT1) encoding an EAR motif-containing zinc-finger protein. The growth inhibition by ZAT1 was shared by ZAT4 and ZAT9, the ZAT1 homologues. The ZAT1 promoter was specifically active in the outermost cells of the root cap, in which ZAT1-GFP was localized when expressed by the ZAT1 promoter. The outermost cell-specific expression pattern of ZAT1 was not altered in the sombrero (smb) or smb bearskin1 (brn1) brn2 accumulating additional root-cap layers. In contrast, ZAT4-GFP and ZAT9- GFP fusion proteins were distributed to the inner root-cap cells in addition to the outermost cells where ZAT4 and ZAT9 promoters were active. Overexpression of ZAT1 induced the ectopic expression of PUTATIVE ASPARTIC PROTEASE3 involved in the programmed cell death. The EAR motif was essential for the growth inhibition by ZAT1. These results suggest that the three related ZATs might regulate the maturation of the outermost cells of the root cap.

Differential Gene Profiling of the Heartwood Formation Process in Taiwania cryptomerioides Hayata Xylem Tissues

Taiwania (Taiwania cryptomerioides) is an important tree species in Taiwan because of the excellent properties of its wood and fascinating color qualities of its heartwood (HW), as well as the bioactive compounds therein. However, limited information is available as to the HW formation of this species. The objective of this research is to analyze the differentially expressed genes (DEGs) during the HW formation process from specific Taiwania xylem tissues, and to obtain genes that might be closely associated with this process. The results indicated that our analyses have captured DEGs representative to the HW formation process of Taiwania. DEGs related to the terpenoid biosynthesis pathway were all up-regulated in the transition zone (TZ) to support the biosynthesis and accumulation of terpenoids. Many DEGs related to lignin biosynthesis, and two DEGs related to pinoresinol reductase (PrR)/pinoresinol lariciresinol reductase (PLR), were up-regulated in TZ. These DEGs together are likely involved in providing the precursors for the subsequent lignan biosynthesis. Several transcription factor-, nuclease-, and protease-encoding DEGs were also highly expressed in TZ, and these DEGs might be involved in the regulation of secondary metabolite biosynthesis and the autolysis of the cellular components of ray parenchyma cells in TZ. These results provide further insights into the process of HW formation in Taiwania.

Morphological and Physiological Framework Underlying Plant Longevity in Arabidopsis thaliana

Monocarpic plants have a single reproductive phase, in which their longevity is developmentally programmed by molecular networks. In the reproductive phase of Arabidopsis thaliana, the inflorescence meristem (IM) maintains a central pool of stem cells and produces a limited number of flower primordia, which result in seed formation and the death of the whole plant. In this study, we observed morphological changes in the IM at cellular and intracellular resolutions until the end of the plant life cycle. We observed four biological events during the periods from 1 week after bolting (WAB) till the death of stem cells: (1) the gradual reduction in the size of the IM, (2) the dynamic vacuolation of IM cells, (3) the loss of the expression of the stem cell determinant WUSCHEL (WUS), and (4) the upregulation of the programmed cell death marker BIFUNCTIONAL NUCLEASE1 (BFN1) in association with the death of stem cells. These results indicate that the stem cell population gradually decreases in IM during plant aging and eventually is fully terminated. We further show that the expression of WUS became undetectable in IM at 3 WAB prior to the loss of CLAVATA3 (CLV3) expression at 5 WAB; CLV3 is a negative regulator of WUS. Moreover, clv3 plants showed delayed loss of WUS and lived 6 weeks longer compared with wild-type plants. These results indicated that the prolonged expression of CLV3 at 4-5 WAB may be a safeguard that inhibits the reactivation of WUS and promotes plant death. Finally, through transcriptome analysis, we determined that reactive oxygen species (ROS) are involved in the control of plant longevity. Our work presents a morphological and physiological framework for the regulation of plant longevity in Arabidopsis.

Genome-Wide Identification and Expression Analysis of the Metacaspase Gene Family in Gossypium Species

Metacaspases (MCs) are cysteine proteases that are important for programmed cell death (PCD) in plants. In this study, we identified 89 MC genes in the genomes of four Gossypium species (Gossypium raimondii, Gossypium barbadense, Gossypium hirsutum, and Gossypium arboreum), and classified them as type-I or type-II genes. All of the type-I and type-II MC genes contain a sequence encoding the peptidase C14 domain. During developmentally regulated PCD, type-II MC genes may play an important role related to fiber elongation, while type-I genes may affect the thickening of the secondary wall. Additionally, 13 genes were observed to be differentially expressed between two cotton lines with differing fiber strengths, and four genes (GhMC02, GhMC04, GhMC07, and GhMC08) were predominantly expressed in cotton fibers at 5–30 days post-anthesis (DPA). During environmentally induced PCD, the expression levels of four genes were affected in the root, stem, and leaf tissues within 6 h of an abiotic stress treatment. In general, the MC gene family affects the development of cotton fibers, including fiber elongation and fiber thickening while four prominent fiber- expressed genes were identified. The effects of the abiotic stress and hormone treatments imply that the cotton MC gene family may be important for fiber development. The data presented herein may form the foundation for future investigations of the MC gene family in Gossypium species.

Extracellular peptide Kratos restricts cell death during vascular development and stress in Arabidopsis

During plant vascular development, xylem tracheary elements (TEs) form water-conducting, empty pipes by genetically regulated cell death. Cell death is prevented from spreading to non-TEs by unidentified intercellular mechanisms, downstream of METACASPASE9 (MC9)-mediated regulation of autophagy in TEs. Here, we identified differentially abundant extracellular peptides in vascular-differentiating wild-type and MC9-down-regulated Arabidopsis cell suspensions. A peptide named Kratos rescued the abnormally high ectopic non-TE death resulting from either MC9 knockout or TE-specific overexpression of the ATG5 autophagy protein during experimentally induced vascular differentiation in Arabidopsis cotyledons. Kratos also reduced cell death following mechanical damage and extracellular ROS production in Arabidopsis leaves. Stress-induced but not vascular non-TE cell death was enhanced by another identified peptide, named Bia. Bia is therefore reminiscent of several known plant cell death-inducing peptides acting as damage-associated molecular patterns. In contrast, Kratos plays a novel extracellular cell survival role in the context of development and during stress response.

Plant proteases during developmental programmed cell death

Proteases are among the key regulators of most forms of programmed cell death (PCD) in animals. Many PCD processes have also been associated with protease expression or activation in plants, However, functional evidence for the roles and actual modes of action of plant proteases in PCD remains surprisingly limited. In this review, we provide an update on protease involvement in the context of developmentally regulated plant PCD. To illustrate the diversity of protease functions, we focus on several prominent developmental PCD processes, including xylem and tapetum maturation, suspensor elimination, endosperm degradation, and seed coat formation, as well as plant senescence processes. Despite the substantial advances in the field, protease functions are often only correlatively linked to developmental PCD, and the specific molecular roles of proteases in many developmental PCD processes remain to be elucidated.

Metacaspase gene family in Rosaceae genomes: Comparative genomic analysis and their expression during pear pollen tube and fruit development

Metacaspase (MC), which is discovered gene family with distant caspase homologs in plants, fungi, and protozoa, may be involved in programmed cell death (PCD) processes during plant development and respond abiotic and biotic stresses. To reveal the evolutionary relationship of MC gene family in Rosaceae genomes, we identified 8, 7, 8, 12, 12, and 23 MC genes in the genomes of Fragaria vesca, Prunus mume, Prunus persica, Pyrus communis, Pyrus bretschneideri and Malus domestica, respectively. Phylogenetic analysis suggested that the MC genes could be grouped into three clades: Type I*, Type I and Type II, which was supported by gene structure and conserved motif analysis. Microsynteny analysis revealed that MC genes present in the corresponding syntenic blocks of P. communis, P. bretschneideri and M. domestica, and further suggested that large-scale duplication events play an important role in the expansion of MC gene family members in these three genomes than other Rosaceae plants (F. vesca, P. mume and P. persica). RNA-seq data showed the specific expression patterns of PbMC genes in response to drought stress. The expression analysis of MC genes demonstrated that PbMC01 and PbMC03 were able to be detected in all four pear pollen tubes and seven fruit development stages. The current study highlighted the evolutionary relationship and duplication of the MC gene family in these six Rosaceae genomes and provided appropriate candidate genes for further studies in P. bretschneideri.

ABA-Induced Vegetative Diaspore Formation in Physcomitrella patens

The phytohormone abscisic acid (ABA) is a pivotal regulator of gene expression in response to various environmental stresses such as desiccation, salt and cold causing major changes in plant development and physiology. Here we show that in the moss Physcomitrella patens exogenous application of ABA triggers the formation of vegetative diaspores (brachycytes or brood cells) that enable plant survival in unfavorable environmental conditions. Such diaspores are round-shaped cells characterized by the loss of the central vacuole, due to an increased starch and lipid storage preparing these cells for growth upon suitable environmental conditions. To gain insights into the gene regulation underlying these developmental and physiological changes, we analyzed early transcriptome changes after 30, 60, and 180 min of ABA application and identified 1,030 differentially expressed genes. Among these, several groups can be linked to specific morphological and physiological changes during diaspore formation, such as genes involved in cell wall modifications. Furthermore, almost all members of ABA-dependent signaling and regulation were transcriptionally induced. Network analysis of transcription-associated genes revealed a large overlap of our study with ABA-dependent regulation in response to dehydration, cold stress, and UV-B light, indicating a fundamental function of ABA in diverse stress responses in moss. We also studied the evolutionary conservation of ABA-dependent regulation between moss and the seed plant Arabidopsis thaliana pointing to an early evolution of ABA-mediated stress adaptation during the conquest of the terrestrial habitat by plants.

The role of KDEL-tailed cysteine endopeptidases of Arabidopsis (AtCEP2 and AtCEP1) in root development

Plants encode a unique group of papain-type cysteine endopeptidases (CysEP) characterized by a C-terminal KDEL endoplasmic reticulum retention signal (KDEL-CysEP) and an unusually broad substrate specificity. The three Arabidopsis KDEL-CysEPs (AtCEP1, AtCEP2, and AtCEP3) are differentially expressed in vegetative and generative tissues undergoing programmed cell death (PCD). While KDEL-CysEPs have been shown to be implicated in the collapse of tissues during PCD, roles of these peptidases in processes other than PCD are unknown. Using mCherry-AtCEP2 and EGFP-AtCEP1 reporter proteins in wild type versus atcep2 or atcep1 mutant plants, we explored the participation of AtCEP in young root development. Loss of AtCEP2, but not AtCEP1 resulted in shorter primary roots due to a decrease in cell length in the lateral root (LR) cap, and impairs extension of primary root epidermis cells such as trichoblasts in the elongation zone. AtCEP2 was localized to root cap corpses adherent to epidermal cells in the rapid elongation zone. AtCEP1 and AtCEP2 are expressed in root epidermis cells that are separated for LR emergence. Loss of AtCEP1 or AtCEP2 caused delayed emergence of LR primordia. KDEL-CysEPs might be involved in developmental tissue remodeling by supporting cell wall elongation and cell separation.

Transcriptional networks orchestrating programmed cell death during plant development

Transcriptional gene regulation is a fundamental biological principle in the development of eukaryotes. It does control not only cell proliferation, specification, and differentiation, but also cell death processes as an integral feature of an organism's developmental program. As in animals, developmentally regulated cell death in plants occurs in numerous contexts and is of vital importance for plant vegetative and reproductive development. In comparison with the information available on the molecular regulation of programmed cell death (PCD) in animals, however, our knowledge on plant PCD still remains scarce. Here, we discuss the functions of different classes of transcription factors that have been implicated in the control of developmentally regulated cell death. Though doubtlessly representing but a first layer of PCD regulation, information on PCD-regulating transcription factors and their targets represents a promising strategy to understand the complex machinery that ensures the precise and failsafe execution of PCD processes in plant development.

Sorghum stem aerenchyma formation is regulated by SbNAC_D during internode development

Sorghum bicolor is a drought-resilient C4 grass used for production of grain, forage, sugar, and biomass. Sorghum genotypes capable of accumulating high levels of stem sucrose have solid stems that contain low levels of aerenchyma. The D-locus on SBI06 modulates the extent of aerenchyma formation in sorghum stems and leaf midribs. A QTL aligned with this locus was identified and fine-mapped in populations derived from BTx623*IS320c, BTx623*R07007, and BTx623*Standard broomcorn. Analysis of coding polymorphisms in the fine-mapped D-locus showed that genotypes that accumulate low levels of aerenchyma encode a truncated NAC transcription factor (Sobic.006G147400, SbNAC_d1), whereas parental lines that accumulate higher levels of stem aerenchyma encode full-length NAC TFs (SbNAC-D). During vegetative stem development, aerenchyma levels are low in nonelongated stem internodes, internode growing zones, and nodes. Aerenchyma levels increase in recently elongated internodes starting at the top of the internode near the center of the stem. SbNAC_D was expressed at low levels in nonelongated internodes and internode growing zones and at higher levels in regions of stem internodes that form aerenchyma. SbXCP1, a gene encoding a cysteine protease involved in programmed cell death, was induced in SbNAC_D genotypes in parallel with aerenchyma formation in sorghum stems but not in SbNAC_d1 genotypes. Several sweet sorghum genotypes encode the recessive SbNAC_d1 allele and have low levels of stem aerenchyma. Based on these results, we propose that SbNAC_D is the D-gene identified by Hilton (1916) and that allelic variation in SbNAC_D modulates the extent of aerenchyma formation in sorghum stems.

Arabidopsis ETHYLENE RESPONSE FACTOR 8 (ERF8) has dual functions in ABA signaling and immunity

BACKGROUND

ETHYLENE RESPONSE FACTOR (ERF) 8 is a member of one of the largest transcription factor families in plants, the APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) superfamily. Members of this superfamily have been implicated in a wide variety of processes such as development and environmental stress responses.

RESULTS

In this study we demonstrated that ERF8 is involved in both ABA and immune signaling. ERF8 overexpression induced programmed cell death (PCD) in Arabidopsis and Nicotiana benthamiana. This PCD was salicylic acid (SA)-independent, suggesting that ERF8 acts downstream or independent of SA. ERF8-induced PCD was abolished by mutations within the ERF-associated amphiphilic repression (EAR) motif, indicating ERF8 induces cell death through its transcriptional repression activity. Two immunity-related mitogen-activated protein kinases, MITOGEN-ACTIVATED PROTEIN KINASE 4 (MPK4) and MPK11, were identified as ERF8-interacting proteins and directly phosphorylated ERF8 in vitro. Four putative MPK phosphorylation sites were identified in ERF8, one of which (Ser103) was determined to be the predominantly phosphorylated residue in vitro, while mutation of all four putative phosphorylation sites partially suppressed ERF8-induced cell death in N. benthamiana. Genome-wide transcriptomic analysis and pathogen growth assays confirmed a positive role of ERF8 in mediating immunity, as ERF8 knockdown or overexpression lines conferred compromised or enhanced resistance against the hemibiotrophic bacterial pathogen Pseudomonas syringae, respectively.

CONCLUSIONS

Together these data reveal that the ABA-inducible transcriptional repressor ERF8 has dual roles in ABA signaling and pathogen defense, and further highlight the complex influence of ABA on plant-microbe interactions.

Transcriptional switch for programmed cell death in pith parenchyma of sorghum stems

Pith parenchyma cells store water in various plant organs. These cells are especially important for producing sugar and ethanol from the sugar juice of grass stems. In many plants, the death of pith parenchyma cells reduces their stem water content. Previous studies proposed that a hypothetical D gene might be responsible for the death of stem pith parenchyma cells in Sorghum bicolor, a promising energy grass, although its identity and molecular function are unknown. Here, we identify the D gene and note that it is located on chromosome 6 in agreement with previous predictions. Sorghum varieties with a functional D allele had stems enriched with dry, dead pith parenchyma cells, whereas those with each of six independent nonfunctional D alleles had stems enriched with juicy, living pith parenchyma cells. D expression was spatiotemporally coupled with the appearance of dead, air-filled pith parenchyma cells in sorghum stems. Among D homologs that are present in flowering plants, Arabidopsis ANAC074 also is required for the death of stem pith parenchyma cells. D and ANAC074 encode previously uncharacterized NAC transcription factors and are sufficient to ectopically induce programmed death of Arabidopsis culture cells via the activation of autolytic enzymes. Taken together, these results indicate that D and its Arabidopsis ortholog, ANAC074, are master transcriptional switches that induce programmed death of stem pith parenchyma cells. Thus, targeting the D gene will provide an approach to breeding crops for sugar and ethanol production.

NAC Transcription Factors ANAC087 and ANAC046 Control Distinct Aspects of Programmed Cell Death in the Arabidopsis Columella and Lateral Root Cap

Programmed cell death in plants occurs both during stress responses and as an integral part of regular plant development. Despite the undisputed importance of developmentally controlled cell death processes for plant growth and reproduction, we are only beginning to understand the underlying molecular genetic regulation. Exploiting the Arabidopsis thaliana root cap as a cell death model system, we identified two NAC transcription factors, the little-characterized ANAC087 and the leaf-senescence regulator ANAC046, as being sufficient to activate the expression of cell death-associated genes and to induce ectopic programmed cell death. In the root cap, these transcription factors are involved in the regulation of distinct aspects of programmed cell death. ANAC087 orchestrates postmortem chromatin degradation in the lateral root cap via the nuclease BFN1. In addition, both ANAC087 and ANAC046 redundantly control the onset of cell death execution in the columella root cap during and after its shedding from the root tip. Besides identifying two regulators of developmental programmed cell death, our analyses reveal the existence of an actively controlled cell death program in Arabidopsis columella root cap cells.

The tomato subtilase family includes several cell death-related proteinases with caspase specificity

Phytaspases are Asp-specific subtilisin-like plant proteases that have been likened to animal caspases with respect to their regulatory function in programmed cell death (PCD). We identified twelve putative phytaspase genes in tomato that differed widely in expression level and tissue-specific expression patterns. Most phytaspase genes are tandemly arranged on tomato chromosomes one, four, and eight, and many belong to taxon-specific clades, e.g. the P69 clade in the nightshade family, suggesting that these genes evolved by gene duplication after speciation. Five tomato phytaspases (SlPhyts) were expressed in N. benthamiana and purified to homogeneity. Substrate specificity was analyzed in a proteomics assay and with a panel of fluorogenic peptide substrates. Similar to animal caspases, SlPhyts recognized an extended sequence motif including Asp at the cleavage site. Clear differences in cleavage site preference were observed implying different substrates in vivo and, consequently, different physiological functions. A caspase-like function in PCD was confirmed for five of the seven tested phytaspases. Cell death was triggered by ectopic expression of SlPhyts 2, 3, 4, 5, 6 in tomato leaves by agro-infiltration, as well as in stably transformed transgenic tomato plants. SlPhyts 3, 4, and 5 were found to contribute to cell death under oxidative stress conditions.

The pattern of 1-aminocyclopropane-1-carboxylate oxidase induction in the tomato leaf petiole abscission zone is independent of expression of the ribonuclease-LX-encoding LeLX gene

The abscission of tomato leaves occurs in the petiole abscission zone, and its late stage includes two spatially divided processes: cell separation and programmed cell death (PCD). Both of these processes are regulated by ethylene. The last step in ethylene biosynthesis is conversion of 1-aminocyclopropane-1-carboxylic acid to ethylene, which is catalysed by the enzyme 1-aminocyclopropane-1-carboxylate oxidase (ACO); however, the location of ACO in the leaf petiole abscission zone is not known. The tomato gene LeLX encodes ribonuclease LX, which is a marker for PCD and is induced by ethylene during abscission, but its association with ACO has not been explored. In a tomato transgenic line 1-7 with inhibited expression of LeLX showing delayed leaf abscission, the morphology and ultrastructure of the leaf petiole abscission zone was examined. In this zone of the cv.'VF36' and of a transgenic line 1-7, spatiotemporal differences in expression of LeACO1 and LeACO4 were analysed and ACO protein was detected immunohistochemically. In comparison to wild-type plants, there were no obvious morphological and ultrastructural features in the abscission zone of plants of a transgenic line 1-7 before and after abscission induction. LeACO1 expression was low before abscission induction, and increased 24 h after induction, although with no apparent spatial pattern. In contrast, LeACO4 was expressed before abscission induction, and its transcript level declined 24 h after induction on the distal side of the abscission zone fracture. In the LeLX-inhibited transgenic line, there were no significant differences in LeACO1 and LeACO4 expression in the petiole abscission zone, in comparison to wild-type plants. In addition, the ACO protein was immunolocalised to the vascular tissues that traverse the petiole abscission zone in plants of wild type and of a transgenic line 1-7; and additionally in the plane of future abscission zone fracture of transgenic-line plants. The results suggest temporal differential expression of the LeACO genes in tomato leaf petioles and vascular localisation of ACO1 protein. Additionally, the results indicate that expression of LeACO genes is not affected by suppression of the LeLX expression.

KIRA1 and ORESARA1 terminate flower receptivity by promoting cell death in the stigma of Arabidopsis

Flowers have a species-specific functional life span that determines the time window in which pollination, fertilization and seed set can occur. The stigma tissue plays a key role in flower receptivity by intercepting pollen and initiating pollen tube growth toward the ovary. In this article, we show that a developmentally controlled cell death programme terminates the functional life span of stigma cells in Arabidopsis. We identified the leaf senescence regulator ORESARA1 (also known as ANAC092) and the previously uncharacterized KIRA1 (also known as ANAC074) as partially redundant transcription factors that modulate stigma longevity by controlling the expression of programmed cell death-associated genes. KIRA1 expression is sufficient to induce cell death and terminate floral receptivity, whereas lack of both KIRA1 and ORESARA1 substantially increases stigma life span. Surprisingly, the extension of stigma longevity is accompanied by only a moderate extension of flower receptivity, suggesting that additional processes participate in the control of the flower's receptive life span.

Transcriptional Roadmap to Seasonal Variation in Wood Formation of Norway Spruce

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

The Role of Programmed Cell Death Regulator LSD1 in Nematode-Induced Syncytium Formation

Cyst-forming plant-parasitic nematodes are common pests of many crops. They inject secretions into host cells to induce the developmental and metabolic reprogramming that leads to the formation of a syncytium, which is the sole food source for growing nematodes. As in other host-parasite models, avirulence leads to rapid and local programmed cell death (PCD) known as the hypersensitive response (HR), whereas in the case of virulence, PCD is still observed but is limited to only some cells. Several regulators of PCD were analyzed to understand the role of PCD in compatible plant-nematode interactions. Thus, Arabidopsis plants carrying recessive mutations in LESION SIMULATING DISEASE1 (LSD1) family genes were subjected to nematode infection assays with juveniles of Heterodera schachtii. LSD1 is a negative and conditional regulator of PCD, and fewer and smaller syncytia were induced in the roots of lsd1 mutants than in wild-type Col-0 plants. Mutation in LSD ONE LIKE2 (LOL2) revealed a pattern of susceptibility to H. schachtii antagonistic to lsd1. Syncytia induced on lsd1 roots compared to Col0 showed significantly retarded growth, modified cell wall structure, increased vesiculation, and some myelin-like bodies present at 7 and 12 days post-infection. To place these data in a wider context, RNA-sequencing analysis of infected and uninfected roots was conducted. During nematode infection, the number of transcripts with changed expression in lsd1 was approximately three times smaller than in wild-type plants (1440 vs. 4206 differentially expressed genes, respectively). LSD1-dependent PCD in roots is thus a highly regulated process in compatible plant-nematode interactions. Two genes identified in this analysis, coding for AUTOPHAGY-RELATED PROTEIN 8F and 8H were down-regulated in syncytia in the presence of LSD1 and showed an increased susceptibility to nematode infection contrasting with lsd1 phenotype. Our data indicate that molecular regulators belonging to the LSD1 family play an important role in precise balancing of diverse PCD players during syncytium development required for successful nematode parasitism.

Transcriptome analysis of embryonic domains in Norway spruce reveals potential regulators of suspensor cell death

The terminal differentiation and elimination of the embryo-suspensor is the earliest manifestation of programmed cell death (PCD) during plant ontogenesis. Molecular regulation of suspensor PCD remains poorly understood. Norway spruce (Picea abies) embryos provide a powerful model for studying embryo development because of their large size, sequenced genome, and the possibility to obtain a large number of embryos at a specific developmental stage through somatic embryogenesis. Here, we have carried out global gene expression analysis of the Norway spruce embryo-suspensor versus embryonal mass (a gymnosperm analogue of embryo proper) using RNA sequencing. We have identified that suspensors have enhanced expression of the NAC domain-containing transcription factors, XND1 and ANAC075, previously shown to be involved in the initiation of developmental PCD in Arabidiopsis. The analysis has also revealed enhanced expression of Norway spruce homologues of the known executioners of both developmental and stress-induced cell deaths, such as metacaspase 9 (MC9), cysteine endopeptidase-1 (CEP1) and ribonuclease 3 (RNS3). Interestingly, a spruce homologue of bax inhibitor-1 (PaBI-1, for Picea abies BI-1), an evolutionarily conserved cell death suppressor, was likewise up-regulated in the embryo-suspensor. Since Arabidopsis BI-1 so far has been implicated only in the endoplasmic reticulum (ER)-stress induced cell death, we investigated its role in embryogenesis and suspensor PCD using RNA interference (RNAi). We have found that PaBI-1-deficient lines formed a large number of abnormal embryos with suppressed suspensor elongation and disturbed polarity. Cytochemical staining of suspensor cells has revealed that PaBI-1 deficiency suppresses vacuolar cell death and induces necrotic type of cell death previously shown to compromise embryo development. This study demonstrates that a large number of cell-death components are conserved between angiosperms and gymnosperms and establishes a new role for BI-1 in the progression of vacuolar cell death.