Arabidopsis BBX14 is involved in high light acclimation and seedling development.

The development of photosynthetically competent seedlings requires both light and retrograde biogenic signaling pathways. The transcription factor GLK1 functions at the interface between these pathways and receives input from the biogenic signal integrator GUN1. BBX14 was previously identified, together with GLK1, in a core module that mediates the response to high light (HL) levels and biogenic signals, which was studied by using inhibitors of chloroplast development. Our chromatin immunoprecipitation-Seq experiments revealed that BBX14 is a direct target of GLK1, and RNA-Seq analysis suggests that BBX14 may function as a regulator of the circadian clock. In addition, BBX14 plays a role in chlorophyll biosynthesis during early onset of light. Knockout of BBX14 results in a long hypocotyl phenotype dependent on a retrograde signal. Furthermore, the expression of BBX14 and BBX15 during biogenic signaling requires GUN1. Investigation of the role of BBX14 and BBX15 in GUN-type biogenic (gun) signaling showed that the overexpression of BBX14 or BBX15 caused de-repression of CA1 mRNA levels, when seedlings were grown on norflurazon. Notably, transcripts of the LHCB1.2 marker are not de-repressed. Furthermore, BBX14 is required to acclimate plants to HL stress. We propose that BBX14 is an integrator of biogenic signals and that BBX14 is a nuclear target of retrograde signals downstream of the GUN1/GLK1 module. However, we do not classify BBX14 or BBX15 overexpressors as gun mutants based on a critical evaluation of our results and those reported in the literature. Finally, we discuss a classification system necessary for the declaration of new gun mutants.

Natural genetic variation in GLK1-mediated photosynthetic acclimation in response to light.

BACKGROUND: GOLDEN-like (GLK) transcription factors are central regulators of chloroplast biogenesis in Arabidopsis and other species. Findings from Arabidopsis show that these factors also contribute to photosynthetic acclimation, e.g. to variation in light intensity, and are controlled by retrograde signals emanating from the chloroplast. However, the natural variation of GLK1-centered gene-regulatory networks in Arabidopsis is largely unexplored. RESULTS: By evaluating the activities of GLK1 target genes and GLK1 itself in vegetative leaves of natural Arabidopsis accessions grown under standard conditions, we uncovered variation in the activity of GLK1 centered regulatory networks. This is linked with the ecogeographic origin of the accessions, and can be associated with a complex genetic variation across loci acting in different functional pathways, including photosynthesis, ROS and brassinosteroid pathways. Our results identify candidate upstream regulators that contribute to a basal level of GLK1 activity in rosette leaves, which can then impact the capacity to acclimate to different environmental conditions. Indeed, accessions with higher GLK1 activity, arising from habitats with a high monthly variation in solar radiation levels, may show lower levels of photoinhibition at higher light intensities. CONCLUSIONS: Our results provide evidence for natural variation in GLK1 regulatory activities in vegetative leaves. This variation is associated with ecogeographic origin and can contribute to acclimation to high light conditions.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system.

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Cell identity specification in plants: lessons from flower development.

Multicellular organisms display a fascinating complexity of cellular identities and patterns of diversification. The concept of 'cell type' aims to describe and categorize this complexity. In this review, we discuss the traditional concept of cell types and highlight the impact of single-cell technologies and spatial omics on the understanding of cellular differentiation in plants. We summarize and compare position-based and lineage-based mechanisms of cell identity specification using flower development as a model system. More than understanding ontogenetic origins of differentiated cells, an important question in plant science is to understand their position- and developmental stage-specific heterogeneity. Combinatorial action and crosstalk of external and internal signals is the key to cellular heterogeneity, often converging on transcription factors that orchestrate gene expression programs.

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks.

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

PUMPKIN, the Sole Plastid UMP Kinase, Associates with Group II Introns and Alters Their Metabolism.

The chloroplast hosts photosynthesis and a variety of metabolic pathways that are essential for plant viability and acclimation processes. In this study, we show that the sole plastid UMP kinase (PUMPKIN) in Arabidopsis (Arabidopsis thaliana) associates specifically with the introns of the plastid transcripts trnG-UCC, trnV-UAC, petB, petD, and ndhA in vivo, as revealed by RNA immunoprecipitation coupled with deep sequencing (RIP-Seq); and that PUMPKIN can bind RNA efficiently in vitro. Analyses of target transcripts showed that PUMPKIN affects their metabolism. Null alleles and knockdowns of pumpkin were viable but clearly affected in growth, plastid translation, and photosynthetic performance. In pumpkin mutants, the levels of many plastid transcripts were reduced, while the amounts of others were increased, as revealed by RNA-Seq analysis. PUMPKIN is a homomultimeric, plastid-localized protein that forms in vivo RNA-containing megadalton-sized complexes and catalyzes the ATP-dependent conversion of UMP to UDP in vitro with properties characteristic of known essential eubacterial UMP kinases. A moonlighting function of PUMPKIN combining RNA and pyrimidine metabolism is discussed.

Differences in DNA Binding Specificity of Floral Homeotic Protein Complexes Predict Organ-Specific Target Genes.

Floral organ identities in plants are specified by the combinatorial action of homeotic master regulatory transcription factors. However, how these factors achieve their regulatory specificities is still largely unclear. Genome-wide in vivo DNA binding data show that homeotic MADS domain proteins recognize partly distinct genomic regions, suggesting that DNA binding specificity contributes to functional differences of homeotic protein complexes. We used in vitro systematic evolution of ligands by exponential enrichment followed by high-throughput DNA sequencing (SELEX-seq) on several floral MADS domain protein homo- and heterodimers to measure their DNA binding specificities. We show that specification of reproductive organs is associated with distinct binding preferences of a complex formed by SEPALLATA3 and AGAMOUS. Binding specificity is further modulated by different binding site spacing preferences. Combination of SELEX-seq and genome-wide DNA binding data allows differentiation between targets in specification of reproductive versus perianth organs in the flower. We validate the importance of DNA binding specificity for organ-specific gene regulation by modulating promoter activity through targeted mutagenesis. Our study shows that intrafamily protein interactions affect DNA binding specificity of floral MADS domain proteins. Differential DNA binding of MADS domain protein complexes plays a role in the specificity of target gene regulation.

The Chloroplast RNA Binding Protein CP31A Has a Preference for mRNAs Encoding the Subunits of the Chloroplast NAD(P)H Dehydrogenase Complex and Is Required for Their Accumulation.

Chloroplast RNA processing requires a large number of nuclear-encoded RNA binding proteins (RBPs) that are imported post-translationally into the organelle. Most of these RBPs are highly specific for one or few target RNAs. By contrast, members of the chloroplast ribonucleoprotein family (cpRNPs) have a wider RNA target range. We here present a quantitative analysis of RNA targets of the cpRNP CP31A using digestion-optimized RNA co-immunoprecipitation with deep sequencing (DO-RIP-seq). This identifies the mRNAs coding for subunits of the chloroplast NAD(P)H dehydrogenase (NDH) complex as main targets for CP31A. We demonstrate using whole-genome gene expression analysis and targeted RNA gel blot hybridization that the ndh mRNAs are all down-regulated in cp31a mutants. This diminishes the activity of the NDH complex. Our findings demonstrate how a chloroplast RNA binding protein can combine functionally related RNAs into one post-transcriptional operon.

AIL and HDG proteins act antagonistically to control cell proliferation.

Aintegumenta-like (AIL) transcription factors are key regulators of cell proliferation and meristem identity. Although AIL functions have been well described, the direct signalling components of this pathway are largely unknown. We show that baby boom (BBM) and other AIL proteins physically interact with multiple members of the L1-expressed homeodomain glabrous (HDG) transcription factor family, including HDG1, HDG11 and HDG12. Overexpression of HDG1, HDG11 and HDG12 restricts growth due to root and shoot meristem arrest, which is associated with reduced expression of genes involved in meristem development and cell proliferation pathways, whereas downregulation of multiple HDG genes promotes cell overproliferation. These results suggest a role for HDG proteins in promoting cell differentiation. We also reveal a transcriptional network in which BBM and HDG1 regulate several common target genes, and where BBM/AIL and HDG regulate the expression of each other. Taken together, these results suggest opposite roles for AIL and HDG proteins, with AILs promoting cell proliferation and HDGs stimulating cell differentiation, and that these functions are mediated at both the protein-protein interaction and transcriptional level.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity.

BACKGROUND: Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system. RESULTS: The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae. CONCLUSIONS: We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

Transcriptomic basis for reinforcement of elm antiherbivore defence mediated by insect egg deposition.

Plant responses to insect egg depositions are known to shape subsequent defensive responses to larvae hatching from the eggs. Elm (Ulmus minor) leaves, on which elm leaf beetles laid their eggs, mount a more efficient defence against larvae hatching from the eggs. However, the molecular mechanisms of this egg-mediated, improved defence are insufficiently understood and have so far only been studied in annual plants. We analysed the dynamics of transcriptomic changes in larval feeding-damaged elm leaves with and without prior egg deposition using de novo assembled RNA-seq data. Compared to egg-free leaves, egg deposition-treated leaves showed earlier and/or faster transcriptional regulations, as well as slightly enhanced differential transcriptional regulation after the onset of larval feeding. These early responding transcripts were overrepresented in gene ontology terms associated with post-translational protein modification, signalling and stress (defence) responses. We found evidence of transcriptional memory in initially egg deposition-induced transcripts whose differential expression was reset prior to larval hatching, but was more rapidly induced again by subsequent larval feeding. This potential memory effect of prior egg deposition, as well as the earlier/faster and enhanced feeding-induced differential regulation of transcripts in egg deposition-treated leaves, may contribute to the egg-mediated reinforcing effect on the elm's defence against larvae. Hence, our study shows that a plant's experience of a stress-indicating environmental cue (here: insect eggs) can push the dynamics of the plant's transcriptomic response to subsequent stress (here: larval feeding). Such experience-mediated acceleration of a stress-induced plant response may result in improved stress resistance.

The BABY BOOM Transcription Factor Activates the LEC1-ABI3-FUS3-LEC2 Network to Induce Somatic Embryogenesis.

Somatic embryogenesis is an example of induced cellular totipotency, where embryos develop from vegetative cells rather than from gamete fusion. Somatic embryogenesis can be induced in vitro by exposing explants to growth regulators and/or stress treatments. The BABY BOOM (BBM) and LEAFY COTYLEDON1 (LEC1) and LEC2 transcription factors are key regulators of plant cell totipotency, as ectopic overexpression of either transcription factor induces somatic embryo formation from Arabidopsis (Arabidopsis thaliana) seedlings without exogenous growth regulators or stress treatments. Although LEC and BBM proteins regulate the same developmental process, it is not known whether they function in the same molecular pathway. We show that BBM transcriptionally regulates LEC1 and LEC2, as well as the two other LAFL genes, FUSCA3 (FUS3) and ABSCISIC ACIDINSENSITIVE3 (ABI3). LEC2 and ABI3 quantitatively regulate BBM-mediated somatic embryogenesis, while FUS3 and LEC1 are essential for this process. BBM-mediated somatic embryogenesis is dose and context dependent, and the context-dependent phenotypes are associated with differential LAFL expression. We also uncover functional redundancy for somatic embryogenesis among other Arabidopsis BBM-like proteins and show that one of these proteins, PLETHORA2, also regulates LAFL gene expression. Our data place BBM upstream of other major regulators of plant embryo identity and totipotency.

Evolution of DNA-Binding Sites of a Floral Master Regulatory Transcription Factor.

Flower development is controlled by the action of key regulatory transcription factors of the MADS-domain family. The function of these factors appears to be highly conserved among species based on mutant phenotypes. However, the conservation of their downstream processes is much less well understood, mostly because the evolutionary turnover and variation of their DNA-binding sites (BSs) among plant species have not yet been experimentally determined. Here, we performed comparative ChIP (chromatin immunoprecipitation)-seq experiments of the MADS-domain transcription factor SEPALLATA3 (SEP3) in two closely related Arabidopsis species: Arabidopsis thaliana and A. lyrata which have very similar floral organ morphology. We found that BS conservation is associated with DNA sequence conservation, the presence of the CArG-box BS motif and on the relative position of the BS to its potential target gene. Differences in genome size and structure can explain that SEP3 BSs in A. lyrata can be located more distantly to their potential target genes than their counterparts in A. thaliana. In A. lyrata, we identified transposition as a mechanism to generate novel SEP3 binding locations in the genome. Comparative gene expression analysis shows that the loss/gain of BSs is associated with a change in gene expression. In summary, this study investigates the evolutionary dynamics of DNA BSs of a floral key-regulatory transcription factor and explores factors affecting this phenomenon.

FRUITFULL controls SAUR10 expression and regulates Arabidopsis growth and architecture.

MADS-domain transcription factors are well known for their roles in plant development and regulate sets of downstream genes that have been uncovered by high-throughput analyses. A considerable number of these targets are predicted to function in hormone responses or responses to environmental stimuli, suggesting that there is a close link between developmental and environmental regulators of plant growth and development. Here, we show that the Arabidopsis MADS-domain factor FRUITFULL (FUL) executes several functions in addition to its noted role in fruit development. Among the direct targets of FUL, we identified SMALL AUXIN UPREGULATED RNA 10 (SAUR10), a growth regulator that is highly induced by a combination of auxin and brassinosteroids and in response to reduced R:FR light. Interestingly, we discovered that SAUR10 is repressed by FUL in stems and inflorescence branches. SAUR10 is specifically expressed at the abaxial side of these branches and this localized activity is influenced by hormones, light conditions and by FUL, which has an effect on branch angle. Furthermore, we identified a number of other genes involved in hormone pathways and light signalling as direct targets of FUL in the stem, demonstrating a connection between developmentally and environmentally regulated growth programs.

The histone deacetylase inhibitor trichostatin a promotes totipotency in the male gametophyte.

The haploid male gametophyte, the pollen grain, is a terminally differentiated structure whose function ends at fertilization. Plant breeding and propagation widely use haploid embryo production from in vitro-cultured male gametophytes, but this technique remains poorly understood at the mechanistic level. Here, we show that histone deacetylases (HDACs) regulate the switch to haploid embryogenesis. Blocking HDAC activity with trichostatin A (TSA) in cultured male gametophytes of Brassica napus leads to a large increase in the proportion of cells that switch from pollen to embryogenic growth. Embryogenic growth is enhanced by, but not dependent on, the high-temperature stress that is normally used to induce haploid embryogenesis in B. napus. The male gametophyte of Arabidopsis thaliana, which is recalcitrant to haploid embryo development in culture, also forms embryogenic cell clusters after TSA treatment. Genetic analysis suggests that the HDAC protein HDA17 plays a role in this process. TSA treatment of male gametophytes is associated with the hyperacetylation of histones H3 and H4. We propose that the totipotency of the male gametophyte is kept in check by an HDAC-dependent mechanism and that the stress treatments used to induce haploid embryo development in culture impinge on this HDAC-dependent pathway.

Arabidopsis JAGGED links floral organ patterning to tissue growth by repressing Kip-related cell cycle inhibitors.

Plant morphogenesis requires coordinated cytoplasmic growth, oriented cell wall extension, and cell cycle progression, but it is debated which of these processes are primary drivers for tissue growth and directly targeted by developmental genes. Here, we used ChIP high-throughput sequencing combined with transcriptome analysis to identify global target genes of the Arabidopsis transcription factor JAGGED (JAG), which promotes growth of the distal region of floral organs. Consistent with the roles of JAG during organ initiation and subsequent distal organ growth, we found that JAG directly repressed genes involved in meristem development, such as CLAVATA1 and HANABA TARANU, and genes involved in the development of the basal region of shoot organs, such as BLADE ON PETIOLE 2 and the GROWTH REGULATORY FACTOR pathway. At the same time, JAG regulated genes involved in tissue polarity, cell wall modification, and cell cycle progression. In particular, JAG directly repressed KIP RELATED PROTEIN 4 (KRP4) and KRP2, which control the transition to the DNA synthesis phase (S-phase) of the cell cycle. The krp2 and krp4 mutations suppressed jag defects in organ growth and in the morphology of petal epidermal cells, showing that the interaction between JAG and KRP genes is functionally relevant. Our work reveals that JAG is a direct mediator between genetic pathways involved in organ patterning and cellular functions required for tissue growth, and it shows that a regulatory gene shapes plant organs by releasing a constraint on S-phase entry.

Structural determinants of DNA recognition by plant MADS-domain transcription factors.

Plant MADS-domain transcription factors act as key regulators of many developmental processes. Despite the wealth of information that exists about these factors, the mechanisms by which they recognize their cognate DNA-binding site, called CArG-box (consensus CCW6GG), and how different MADS-domain proteins achieve DNA-binding specificity, are still largely unknown. We used information from in vivo ChIP-seq experiments, in vitro DNA-binding data and evolutionary conservation to address these important questions. We found that structural characteristics of the DNA play an important role in the DNA binding of plant MADS-domain proteins. The central region of the CArG-box largely resembles a structural motif called 'A-tract', which is characterized by a narrow minor groove and may assist bending of the DNA by MADS-domain proteins. Periodically spaced A-tracts outside the CArG-box suggest additional roles for this structure in the process of DNA binding of these transcription factors. Structural characteristics of the CArG-box not only play an important role in DNA-binding site recognition of MADS-domain proteins, but also partly explain differences in DNA-binding specificity of different members of this transcription factor family and their heteromeric complexes.

Characterization of MADS-domain transcription factor complexes in Arabidopsis flower development.

Floral organs are specified by the combinatorial action of MADS-domain transcription factors, yet the mechanisms by which MADS-domain proteins activate or repress the expression of their target genes and the nature of their cofactors are still largely unknown. Here, we show using affinity purification and mass spectrometry that five major floral homeotic MADS-domain proteins (AP1, AP3, PI, AG, and SEP3) interact in floral tissues as proposed in the "floral quartet" model. In vitro studies confirmed a flexible composition of MADS-domain protein complexes depending on relative protein concentrations and DNA sequence. In situ bimolecular fluorescent complementation assays demonstrate that MADS-domain proteins interact during meristematic stages of flower development. By applying a targeted proteomics approach we were able to establish a MADS-domain protein interactome that strongly supports a mechanistic link between MADS-domain proteins and chromatin remodeling factors. Furthermore, members of other transcription factor families were identified as interaction partners of floral MADS-domain proteins suggesting various specific combinatorial modes of action.

Arabidopsis class I and class II TCP transcription factors regulate jasmonic acid metabolism and leaf development antagonistically.

TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR1 (TCP) transcription factors control developmental processes in plants. The 24 TCP transcription factors encoded in the Arabidopsis (Arabidopsis thaliana) genome are divided into two classes, class I and class II TCPs, which are proposed to act antagonistically. We performed a detailed phenotypic analysis of the class I tcp20 mutant, showing an increase in leaf pavement cell sizes in 10-d-old seedlings. Subsequently, a glucocorticoid receptor induction assay was performed, aiming to identify potential target genes of the TCP20 protein during leaf development. The LIPOXYGENASE2 (LOX2) and class I TCP9 genes were identified as TCP20 targets, and binding of TCP20 to their regulatory sequences could be confirmed by chromatin immunoprecipitation analyses. LOX2 encodes for a jasmonate biosynthesis gene, which is also targeted by class II TCP proteins that are under the control of the microRNA JAGGED AND WAVY (JAW), although in an antagonistic manner. Mutation of TCP9, the second identified TCP20 target, resulted in increased pavement cell sizes during early leaf developmental stages. Analysis of senescence in the single tcp9 and tcp20 mutants and the tcp9tcp20 double mutants showed an earlier onset of this process in comparison with wild-type control plants in the double mutant only. Both the cell size and senescence phenotypes are opposite to the known class II TCP mutant phenotype in JAW plants. Altogether, these results point to an antagonistic function of class I and class II TCP proteins in the control of leaf development via the jasmonate signaling pathway.

PRI-CAT: a web-tool for the analysis, storage and visualization of plant ChIP-seq experiments.

Although several tools for the analysis of ChIP-seq data have been published recently, there is a growing demand, in particular in the plant research community, for computational resources with which such data can be processed, analyzed, stored, visualized and integrated within a single, user-friendly environment. To accommodate this demand, we have developed PRI-CAT (Plant Research International ChIP-seq analysis tool), a web-based workflow tool for the management and analysis of ChIP-seq experiments. PRI-CAT is currently focused on Arabidopsis, but will be extended with other plant species in the near future. Users can directly submit their sequencing data to PRI-CAT for automated analysis. A QuickLoad server compatible with genome browsers is implemented for the storage and visualization of DNA-binding maps. Submitted datasets and results can be made publicly available through PRI-CAT, a feature that will enable community-based integrative analysis and visualization of ChIP-seq experiments. Secondary analysis of data can be performed with the aid of GALAXY, an external framework for tool and data integration. PRI-CAT is freely available at http://www.ab.wur.nl/pricat. No login is required.