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.

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.

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.

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.

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.

Characterization of in vivo DNA-binding events of plant transcription factors by ChIP-seq: experimental protocol and computational analysis

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is a powerful technique for genome-wide identification of in vivo binding sites of DNA-binding proteins. The technique had been used to study many DNA-binding proteins in a broad variety of species. The basis of the ChIP-seq technique is the ability to covalently cross-link DNA and proteins that are located in very close proximity. This allows the use of an antibody against the (tagged) protein of interest to specifically enrich DNA-fragments bound by this protein. ChIP-seq can be performed using antibodies against the native protein or against tagged proteins. Using a specific antibody against a tag to immunoprecipitate tagged proteins eliminates the need for a specific antibody against the native protein and allows more experimental flexibility. In this chapter we present a complete workflow for experimental procedure and bioinformatic analysis that allows wet-lab biologists to perform and analyze ChIP-seq experiments.

Dynamics of chromatin accessibility and gene regulation by MADS-domain transcription factors in flower development

Background

Development of eukaryotic organisms is controlled by transcription factors that trigger specific and global changes in gene expression programs. In plants, MADS-domain transcription factors act as master regulators of developmental switches and organ specification. However, the mechanisms by which these factors dynamically regulate the expression of their target genes at different developmental stages are still poorly understood.

Results

We characterized the relationship of chromatin accessibility, gene expression, and DNA binding of two MADS-domain proteins at different stages of Arabidopsis flower development. Dynamic changes in APETALA1 and SEPALLATA3 DNA binding correlated with changes in gene expression, and many of the target genes could be associated with the developmental stage in which they are transcriptionally controlled. We also observe dynamic changes in chromatin accessibility during flower development. Remarkably, DNA binding of APETALA1 and SEPALLATA3 is largely independent of the accessibility status of their binding regions and it can precede increases in DNA accessibility. These results suggest that APETALA1 and SEPALLATA3 may modulate chromatin accessibility, thereby facilitating access of other transcriptional regulators to their target genes.

Conclusions

Our findings indicate that different homeotic factors regulate partly overlapping, yet also distinctive sets of target genes in a partly stage-specific fashion. By combining the information from DNA-binding and gene expression data, we are able to propose models of stage-specific regulatory interactions, thereby addressing dynamics of regulatory networks throughout flower development. Furthermore, MADS-domain TFs may regulate gene expression by alternative strategies, one of which is modulation of chromatin accessibility.

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.

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 CCW₆GG), 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.

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.

ChIP-seq Analysis in R (CSAR): An R package for the statistical detection of protein-bound genomic regions

BACKGROUND

In vivo detection of protein-bound genomic regions can be achieved by combining chromatin-immunoprecipitation with next-generation sequencing technology (ChIP-seq). The large amount of sequence data produced by this method needs to be analyzed in a statistically proper and computationally efficient manner. The generation of high copy numbers of DNA fragments as an artifact of the PCR step in ChIP-seq is an important source of bias of this methodology.

RESULTS

We present here an R package for the statistical analysis of ChIP-seq experiments. Taking the average size of DNA fragments subjected to sequencing into account, the software calculates single-nucleotide read-enrichment values. After normalization, sample and control are compared using a test based on the ratio test or the Poisson distribution. Test statistic thresholds to control the false discovery rate are obtained through random permutations. Computational efficiency is achieved by implementing the most time-consuming functions in C++ and integrating these in the R package. An analysis of simulated and experimental ChIP-seq data is presented to demonstrate the robustness of our method against PCR-artefacts and its adequate control of the error rate.

CONCLUSIONS

The software ChIP-seq Analysis in R (CSAR) enables fast and accurate detection of protein-bound genomic regions through the analysis of ChIP-seq experiments. Compared to existing methods, we found that our package shows greater robustness against PCR-artefacts and better control of the error rate.

Orchestration of floral initiation by APETALA1

The MADS-domain transcription factor APETALA1 (AP1) is a key regulator of Arabidopsis flower development. To understand the molecular mechanisms underlying AP1 function, we identified its target genes during floral initiation using a combination of gene expression profiling and genome-wide binding studies. Many of its targets encode transcriptional regulators, including known floral repressors. The latter genes are down-regulated by AP1, suggesting that it initiates floral development by abrogating the inhibitory effects of these genes. Although AP1 acts predominantly as a transcriptional repressor during the earliest stages of flower development, at more advanced stages it also activates regulatory genes required for floral organ formation, indicating a dynamic mode of action. Our results further imply that AP1 orchestrates floral initiation by integrating growth, patterning, and hormonal pathways.

Chromatin immunoprecipitation (ChIP) of plant transcription factors followed by sequencing (ChIP-SEQ) or hybridization to whole genome arrays (ChIP-CHIP)

Chromatin immunoprecipitation (ChIP) is a powerful technique to study interactions between transcription factors (TFs) and DNA in vivo. For genome-wide de novo discovery of TF-binding sites, the DNA that is obtained in ChIP experiments needs to be processed for sequence identification. The sequences can be identified by direct sequencing (ChIP-SEQ) or hybridization to microarrays (ChIP-CHIP). Given the small amounts of DNA that are usually obtained in ChIP experiments, successful and reproducible sample processing is challenging. Here we provide a detailed procedure for ChIP of plant TFs, as well as protocols for sample preparation for ChIP-SEQ and for ChIP-CHIP. Our ChIP procedure is optimized for high signal-to-noise ratio starting with tissue fixation, followed by nuclei isolation, immunoprecipitation, DNA amplification and purification. We also provide a guide for primary data analysis of ChIP-SEQ data. The complete protocol for ChIP-SEQ/ChIP-CHIP sample preparation starting from plant harvest takes approximately 7 d.

Target genes of the MADS transcription factor SEPALLATA3: integration of developmental and hormonal pathways in the Arabidopsis flower

The molecular mechanisms by which floral homeotic genes act as major developmental switches to specify the identity of floral organs are still largely unknown. Floral homeotic genes encode transcription factors of the MADS-box family, which are supposed to assemble in a combinatorial fashion into organ-specific multimeric protein complexes. Major mediators of protein interactions are MADS-domain proteins of the SEPALLATA subfamily, which play a crucial role in the development of all types of floral organs. In order to characterize the roles of the SEPALLATA3 transcription factor complexes at the molecular level, we analyzed genome-wide the direct targets of SEPALLATA3. We used chromatin immunoprecipitation followed by ultrahigh-throughput sequencing or hybridization to whole-genome tiling arrays to obtain genome-wide DNA-binding patterns of SEPALLATA3. The results demonstrate that SEPALLATA3 binds to thousands of sites in the genome. Most potential target sites that were strongly bound in wild-type inflorescences are also bound in the floral homeotic agamous mutant, which displays only the perianth organs, sepals, and petals. Characterization of the target genes shows that SEPALLATA3 integrates and modulates different growth-related and hormonal pathways in a combinatorial fashion with other MADS-box proteins and possibly with non-MADS transcription factors. In particular, the results suggest multiple links between SEPALLATA3 and auxin signaling pathways. Our gene expression analyses link the genomic binding site data with the phenotype of plants expressing a dominant repressor version of SEPALLATA3, suggesting that it modulates auxin response to facilitate floral organ outgrowth and morphogenesis. Furthermore, the binding of the SEPALLATA3 protein to cis-regulatory elements of other MADS-box genes and expression analyses reveal that this protein is a key component in the regulatory transcriptional network underlying the formation of floral organs.

Most regulatory genes encode transcription factors, which modulate gene expression by binding to regulatory sequences of their target genes. In plants in particular, which genes are directly controlled by these transcription factors, and the molecular mechanisms of target gene recognition in vivo, are still largely unexplored. One of the best-understood developmental processes in plants is flower development. In different combinations, transcription factors of the MADS-box family control the identities of the different types of floral organs: sepals, petals, stamens, and carpels. Here, we present the first genome-wide analysis of binding sites of a MADS-box transcription factor in plants. We show that the MADS-domain protein SEPALLATA3 (SEP3) binds to the regulatory regions of thousands of potential target genes, many of which are also transcription factors. We provide insight into mechanisms of DNA recognition by SEP3, and suggest roles for other transcription factor families in SEP3 target gene regulation. In addition to effects on genes involved in floral organ identity, our data suggest that SEP3 binds to, and modulates, the transcription of target genes involved in hormonal signaling pathways.

The key floral regulator SEPALLATA3 binds to the promoters of a large number of potential direct target genes to integrate different growth-related and hormonal pathways in flower development.