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.

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.

Estimating DNA-Binding Specificities of Transcription Factors Using SELEX-Seq.

Here we provide an updated protocol for the Systematic Evolution of Ligands followed by massively parallel sequencing (SELEX-seq) method to study protein-DNA interaction specificities. This in vitro method is used to characterize DNA-binding specificities of transcription factors (TFs). The procedure is based on cycles of immunoprecipitation of protein-DNA complexes, starting with a randomized DNA library of defined fragment length, followed by massively parallel sequencing. The updated protocol includes aspects of experimental design and procedure as well as basic instructions on data analysis.

A Hands-On Guide to Generate Spatial Gene Expression Profiles by Integrating scRNA-seq and 3D-Reconstructed Microscope-Based Plant Structures.

Transcriptome profiles of individual cells in the plant are strongly dependent on their relative position. Cell differentiation is associated with tissue-specific transcriptomic changes. For that reason, it is important to study gene expression changes in a spatial context, and therefore to link those to potential morphological changes over developmental time. Even though great experimental advances have been made in recording spatial gene expression profiles, those attempts are limited in the plant field. New computational approaches attempt to solve this problem by integrating spatial expression profiles of few marker genes with single-cell/single-nuclei RNA-seq (scRNA-seq) methodologies. In this chapter, we provide a practical guide on how to predict gene expression patterns in a 3D plant structure by combining scRNA-seq data and 3D microscope-based reconstructed expression profiles of a small set of reference genes. We also show how to visualize these results.

Dual specificity and target gene selection by the MADS-domain protein FRUITFULL.

How transcription factors attain their target gene specificity and how this specificity may be modulated, acquiring different regulatory functions through the development of plant tissues, is an open question. Here we characterized different regulatory roles of the MADS-domain transcription factor FRUITFULL (FUL) in flower development and mechanisms modulating its activity. We found that the dual role of FUL in regulating floral transition and pistil development is associated with its different in vivo patterns of DNA binding in both tissues. Characterization of FUL protein complexes by liquid chromatography-tandem mass spectrometry and SELEX-seq experiments shows that aspects of tissue-specific target site selection can be predicted by tissue-specific variation in the composition of FUL protein complexes with different DNA binding specificities, without considering the chromatin status of the target region. This suggests a role for dynamic changes in FUL TF complex composition in reshaping the regulatory functions of FUL during flower development.

A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data.

Cellular heterogeneity in growth and differentiation results in organ patterning. Single-cell transcriptomics allows characterization of gene expression heterogeneity in developing organs at unprecedented resolution. However, the original physical location of the cell is lost during this methodology. To recover the original location of cells in the developing organ is essential to link gene activity with cellular identity and function in plants. Here, we propose a method to reconstruct genome-wide gene expression patterns of individual cells in a 3D flower meristem by combining single-nuclei RNA-seq with microcopy-based 3D spatial reconstruction. By this, gene expression differences among meristematic domains giving rise to different tissue and organ types can be determined. As a proof of principle, the method is used to trace the initiation of vascular identity within the floral meristem. Our work demonstrates the power of spatially reconstructed single cell transcriptome atlases to understand plant morphogenesis. The floral meristem 3D gene expression atlas can be accessed at http://threed-flower-meristem.herokuapp.com .

The intervening domain is required for DNA-binding and functional identity of plant MADS transcription factors.

The MADS transcription factors (TF) are an ancient eukaryotic protein family. In plants, the family is divided into two main lineages. Here, we demonstrate that DNA binding in both lineages absolutely requires a short amino acid sequence C-terminal to the MADS domain (M domain) called the Intervening domain (I domain) that was previously defined only in type II lineage MADS. Structural elucidation of the MI domains from the floral regulator, SEPALLATA3 (SEP3), shows a conserved fold with the I domain acting to stabilise the M domain. Using the floral organ identity MADS TFs, SEP3, APETALA1 (AP1) and AGAMOUS (AG), domain swapping demonstrate that the I domain alters genome-wide DNA-binding specificity and dimerisation specificity. Introducing AG carrying the I domain of AP1 in the Arabidopsis ap1 mutant resulted in strong complementation and restoration of first and second whorl organs. Taken together, these data demonstrate that the I domain acts as an integral part of the DNA-binding domain and significantly contributes to the functional identity of the MADS TF.

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.

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.

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.

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.

Ectopic Transplastomic Expression of a Synthetic MatK Gene Leads to Cotyledon-Specific Leaf Variegation.

Chloroplasts (and other plastids) harbor their own genetic material, with a bacterial-like gene-expression systems. Chloroplast RNA metabolism is complex and is predominantly mediated by nuclear-encoded RNA-binding proteins. In addition to these nuclear factors, the chloroplast-encoded intron maturase MatK has been suggested to perform as a splicing factor for a subset of chloroplast introns. MatK is essential for plant cell survival in tobacco, and thus null mutants have not yet been isolated. We therefore attempted to over-express MatK from a neutral site in the chloroplast, placing it under the control of a theophylline-inducible riboswitch. This ectopic insertion of MatK lead to a variegated cotyledons phenotype. The addition of the inducer theophylline exacerbated the phenotype in a concentration-dependent manner. The extent of variegation was further modulated by light, sucrose and spectinomycin, suggesting that the function of MatK is intertwined with photosynthesis and plastid translation. Inhibiting translation in the transplastomic lines has a profound effect on the accumulation of several chloroplast mRNAs, including the accumulation of an RNA antisense to rpl33, a gene coding for an essential chloroplast ribosomal protein. Our study further supports the idea that MatK expression needs to be tightly regulated to prevent detrimental effects and establishes another link between leaf variegation and chloroplast translation.

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.

Long ncRNA A-ROD activates its target gene DKK1 at its release from chromatin.

Long ncRNAs are often enriched in the nucleus and at chromatin, but whether their dissociation from chromatin is important for their role in transcription regulation is unclear. Here, we group long ncRNAs using epigenetic marks, expression and strength of chromosomal interactions; we find that long ncRNAs transcribed from loci engaged in strong long-range chromosomal interactions are less abundant at chromatin, suggesting the release from chromatin as a crucial functional aspect of long ncRNAs in transcription regulation of their target genes. To gain mechanistic insight into this, we functionally validate the long ncRNA A-ROD, which enhances DKK1 transcription via its nascent spliced released form. Our data provide evidence that the regulatory interaction requires dissociation of A-ROD from chromatin, with target specificity ensured within the pre-established chromosomal proximity. We propose that the post-transcriptional release of a subset of long ncRNAs from the chromatin-associated template plays an important role in their function as transcription regulators.

Profiling Nucleosome Occupancy by MNase-seq: Experimental Protocol and Computational Analysis.

Nucleosomes are the basic repeating units of eukaryotic chromatin. They play important roles in chromatin compaction and gene regulation. Therefore, it is important to profile the in vivo locations of nucleosomes in the genome. Here we illustrate how to profile nucleosome occupancy at genome-wide scale using micrococcal nuclease (MNase) digestion combined with high throughput Illumina sequencing (MNase-seq). Nucleosome-associated DNA is relatively insensitive to digestion by micrococcal nuclease (MNase). Upon mild MNase treatment, the undigested nucleosomal DNA can be purified and sequenced allowing a precise localization of in vivo nucleosomes at a genome-wide level.

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.

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.

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.