DNA-binding specificities of plant transcription factors and their potential to define target genes

Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes

Time-series transcriptomes of a biological process obtained under different conditions are useful for identifying the regulators of the process and their regulatory networks. However, such data are 3D (gene expression, time, and condition), and there is currently no method that can deal with their full complexity. Here, we developed a method that avoids time-point alignment and normalization between conditions. We applied it to analyze time-series transcriptomes of developing maize leaves under light-dark cycles and under total darkness and obtained eight time-ordered gene coexpression networks (TO-GCNs), which can be used to predict upstream regulators of any genes in the GCNs. One of the eight TO-GCNs is light-independent and likely includes all genes involved in the development of Kranz anatomy, which is a structure crucial for the high efficiency of photosynthesis in C₄ plants. Using this TO-GCN, we predicted and experimentally validated a regulatory cascade upstream of SHORTROOT1, a key Kranz anatomy regulator. Moreover, we applied the method to compare transcriptomes from maize and rice leaf segments and identified regulators of maize C₄ enzyme genes and RUBISCO SMALL SUBUNIT2 Our study provides not only a powerful method but also novel insights into the regulatory networks underlying Kranz anatomy development and C₄ photosynthesis.

Epigenetic control of UV-B-induced flavonoid accumulation in Artemisia annua L

UV-B-induced flavonoid biosynthesis is epigenetically regulated by site-specific demethylation of AaMYB1, AaMYC, and AaWRKY TF-binding sites inAaPAL1promoter-causing overexpression ofAaPALgene inArtemisia annua. The present study was undertaken to understand the epigenetic regulation of flavonoid biosynthesis under the influence of ultraviolet-B radiation using Artemisia annua L. as an experimental model. In-vitro propagated and acclimatized plantlets were treated with UV-B radiation (2.8 W m^-2; 3 h), which resulted in enhanced accumulation of total flavonoid and phenolics content as well as eleven individual flavonoids measured through HPLC-DAC. Expression of eight genes (phenylanaline ammonia lyase, cinnamate-4-hydroxylase, 4-coumarate: CoA ligase; chalcone synthase, chalcone isomerase, cinnamoyl reductase, flavonoid-3'-hydroxylase, and flavones synthase) from upstream and downstream flavonoid biosynthetic pathways was measured through RT-PCR and RT-Q-PCR and all were variably induced under UV-B irradiation. Among them, AaPAL1 transcript and its protein were most significantly upregulated. Global DNA methylation analysis revealed hypomethylation of genomic DNA in A. annua. Further epigenetic characterization of promoter region of AaPAL1 revealed cytosine demethylation at five sites, which in turn caused epigenetic activation of six transcription factor-binding sites including QELEMENT, EBOXBNNAPA/MYCCONSENSUSAT, MYBCORE, MYBCOREATCYCB1, and GCCCORE. MYB transcription factors are positive regulators of flavonoid biosynthesis. Epigenetic activation of transcription-enhancing cis-regulatory elements in AaPAL1 promoter and subsequent overexpression of AaMYB1 and AaMYC and AaWRKY transcription factors under UV-B irradiation may probably be the reason for higher AaPAL1 expression and hence greater biosynthesis of flavonoids in A. annua L. The present study is the first report that provides mechanistic evidence of epigenetic regulation of flavonoid biosynthesis under UV-B radiation in A. annua L.

A Regulatory Module Controlling GA-Mediated Endosperm Cell Expansion Is Critical for Seed Germination in Arabidopsis

A key component of seed germination is the interplay of mechanical forces governing embryo growth and the surrounding restraining endosperm tissue. Endosperm cell separation is therefore thought to play a critical role in the control of this developmental transition. Here we demonstrate that in Arabidopsis thaliana seeds, endosperm cell expansion is a key component of germination. Endosperm cells expand to accommodate embryo growth prior to germination. We show that this is an actively regulated process supported by spatiotemporal control of the cell expansion gene EXPANSIN 2 (EXPA2). The NAC transcription factors NAC25 and NAC1L were identified as upstream regulators of EXPA2 expression, gibberellin-mediated endosperm expansion, and seed germination. The DELLA protein RGL2 repressed activation of the EXPA2 promoter by NAC25/NAC1L. Taken together, our findings uncover a key role of the GA/DELLA-NAC25/NAC1L-EXPA2 network in regulating endosperm cell expansion to control the seed-to-seedling transition.

Functional dissection of the ARGONAUTE7 promoter

ARGONAUTES are the central effector proteins of RNA silencing which bind target transcripts in a small RNA-guided manner. Arabidopsis thaliana has 10 ARGONAUTE (AGO) genes, with specialized roles in RNA-directed DNA methylation, post-transcriptional gene silencing, and antiviral defense. To better understand specialization among AGO genes at the level of transcriptional regulation we tested a library of 1497 transcription factors for binding to the promoters of AGO1,AGO10, and AGO7 using yeast 1-hybrid assays. A ranked list of candidate DNA-binding TFs revealed binding of the AGO7 promoter by a number of proteins in two families: the miR156-regulated SPL family and the miR319-regulated TCP family, both of which have roles in developmental timing and leaf morphology. Possible functions for SPL and TCP binding are unclear: we showed that these binding sites are not required for the polar expression pattern of AGO7, nor for the function of AGO7 in leaf shape. Normal AGO7 transcription levels and function appear to depend instead on an adjacent 124-bp region. Progress in understanding the structure of this promoter may aid efforts to understand how the conserved AGO7-triggered TAS3 pathway functions in timing and polarity.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

Enriched networks 'nucleoside/nucleotide and ribonucleoside/ribonucleotide metabolic processes' and 'response to stimulus' potentially conferred to drought adaptation of the epiphytic orchid Dendrobium wangliangii

Dendrobium wangliangii is an endangered and epiphytic orchid with tolerance to seasonally extreme arid conditions and occurs exclusively in the hot-dry valley area of southwestern China. To reveal its molecular basis responsible for ecological adaptation, large-scale transcriptome sequencing was performed using Illumina sequencing with pooled mRNA extracted from whole plants and pseudobulbs during drought and rainy seasons. Based on the target transcript selection, the differentially expressed genes were related to 8 well-known drought-tolerant categories, and to morphological traits in resistance to water stress including pseudobulbs and roots. Further gene ontology enrichment analysis revealed that 'nucleoside/nucleotide and ribonucleoside/ribonucleotide metabolic processes' and 'response to stimulus' were the two most important aspects in resistance to drought stress with respect to the whole plant. In addition, the difference in the number and category of differentially expressed genes in whole plant and stem suggested the involvement of genes specifically localized in the stem, such as GTP-binding protein, lipases, signaling related transcripts and those involved in the ATP metabolic process. The comprehensive analysis of the epiphytic orchid in response to water deprivation indicates that integral tactics lead to active adaptation as a basal defense response to drought stress by the endangered epiphyte, including the collaboration of metabolic processes, responses to a various stimulus and other candidate genes contribute to its extreme drought tolerance. Insights from this study can be further utilized to understand stress-responsive genes in other medicinally important species and to improve the drought tolerance of food crops.

Analysis of bHLH coding genes of Cicer arietinum during heavy metal stress using biological network

bHLH family of transcription factors play important role in regulating many cellular and physiological functions in plants. These proteins are also known to be involved in response to several abiotic stress types. Cicer arietinum is an important source of protein in food across the globe. Considerable differential expression in the bHLH family of proteins during heavy metal exposure in Cicer arietinum was observed by microarray data analysis. The study aimed to construct a Pearson coefficient correlation based network of bHLH coding genes in the plant. Microarray data of Cicer arietinum recorded under cadmium and chromium stress (GSE86807) from GEO at NCBI was used for analysis. The network constructed from expression data set of the 85 bHLH coding genes revealed 10 hub genes that are connected with topological genes. These hub genes are stress responsive genes that may also be regarded as the marker genes for heavy metal response. Our analysis reported a new set of reference genes (hub genes) that have potentially significant role in development of stress tolerant crops.

AP2/ERF Transcription Factor Regulatory Networks in Hormone and Abiotic Stress Responses in Arabidopsis

Dynamic environmental changes such as extreme temperature, water scarcity and high salinity affect plant growth, survival, and reproduction. Plants have evolved sophisticated regulatory mechanisms to adapt to these unfavorable conditions, many of which interface with plant hormone signaling pathways. Abiotic stresses alter the production and distribution of phytohormones that in turn mediate stress responses at least in part through hormone- and stress-responsive transcription factors. Among these, the APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) family transcription factors (AP2/ERFs) have emerged as key regulators of various stress responses, in which they also respond to hormones with improved plant survival during stress conditions. Apart from participation in specific stresses, AP2/ERFs are involved in a wide range of stress tolerance, enabling them to form an interconnected stress regulatory network. Additionally, many AP2/ERFs respond to the plant hormones abscisic acid (ABA) and ethylene (ET) to help activate ABA and ET dependent and independent stress-responsive genes. While some AP2/ERFs are implicated in growth and developmental processes mediated by gibberellins (GAs), cytokinins (CTK), and brassinosteroids (BRs). The involvement of AP2/ERFs in hormone signaling adds the complexity of stress regulatory network. In this review, we summarize recent studies on AP2/ERF transcription factors in hormonal and abiotic stress responses with an emphasis on selected family members in Arabidopsis. In addition, we leverage publically available Arabidopsis gene networks and transcriptome data to investigate AP2/ERF regulatory networks, providing context and important clues about the roles of diverse AP2/ERFs in controlling hormone and stress responses.

The SPB-Box Transcription Factor AaSPL2 Positively Regulates Artemisinin Biosynthesis in Artemisia annua L

Artemisinin, an important compound produced by Artemisia annua, is the active ingredient in the treatment of malaria. Jasmonic acid, one of the phytohormones, is an important elicitor of artemisinin biosynthesis by enhancing transcription levels of transcription factors. SPL transcription factors are plant-specific transcription factors of plant growth, development, and secondary metabolism regulation. However, to date, the SPL transcription factors that regulate artemisinin biosynthesis is currently unclear. Here, we show that an SPL transcription factor can positively regulate artemisinin biosynthesis by binding to the promoter of artemisinin biosynthetic pathway genes. We screened AaSPL2 by gene expression profiles analysis in 14 SPL transcription factors. We demonstrated that AaSPL2 can activate the promoter of DBR2 by dual-LUC assy. Moreover, in the AaSPL2 overexpression plants, the artemisinin content was increased by 33-86%, and in the AaSPL2 -RNAi transgenic plants, artemisinin content was decreased by 33-65%. These data suggest that AaSPL2 and DBR2 interact with a "GTAC" cis-element in the DBR2 promoter, mediating the transcriptional activation of DBR2 in response to JA and resulting in the improvement on artemisinin content.

AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

As an economically important crop, apple is one of the most cultivated fruit trees in temperate regions worldwide. Recently, a large number of high-quality transcriptomic and epigenomic datasets for apple were made available to the public, which could be helpful in inferring gene regulatory relationships and thus predicting gene function at the genome level. Through integration of the available apple genomic, transcriptomic, and epigenomic datasets, we constructed co-expression networks, identified functional modules, and predicted chromatin states. A total of 112 RNA-seq datasets were integrated to construct a global network and a conditional network (tissue-preferential network). Furthermore, a total of 1,076 functional modules with closely related gene sets were identified to assess the modularity of biological networks and further subjected to functional enrichment analysis. The results showed that the function of many modules was related to development, secondary metabolism, hormone response, and transcriptional regulation. Transcriptional regulation is closely related to epigenetic marks on chromatin. A total of 20 epigenomic datasets, which included ChIP-seq, DNase-seq, and DNA methylation analysis datasets, were integrated and used to classify chromatin states. Based on the ChromHMM algorithm, the genome was divided into 620,122 fragments, which were classified into 24 states according to the combination of epigenetic marks and enriched-feature regions. Finally, through the collaborative analysis of different omics datasets, the online database AppleMDO (http://bioinformatics.cau.edu.cn/AppleMDO/) was established for cross-referencing and the exploration of possible novel functions of apple genes. In addition, gene annotation information and functional support toolkits were also provided. Our database might be convenient for researchers to develop insights into the function of genes related to important agronomic traits and might serve as a reference for other fruit trees.

Transcription factor DUO1 generated by neo-functionalization is associated with evolution of sperm differentiation in plants

Evolutionary mechanisms underlying innovation of cell types have remained largely unclear. In multicellular eukaryotes, the evolutionary molecular origin of sperm differentiation is unknown in most lineages. Here, we report that in algal ancestors of land plants, changes in the DNA-binding domain of the ancestor of the MYB transcription factor DUO1 enabled the recognition of a new cis-regulatory element. This event led to the differentiation of motile sperm. After neo-functionalization, DUO1 acquired sperm lineage-specific expression in the common ancestor of land plants. Subsequently the downstream network of DUO1 was rewired leading to sperm with distinct morphologies. Conjugating green algae, a sister group of land plants, accumulated mutations in the DNA-binding domain of DUO1 and lost sperm differentiation. Our findings suggest that the emergence of DUO1 was the defining event in the evolution of sperm differentiation and the varied modes of sexual reproduction in the land plant lineage.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean’s agronomic traits.

Transcription factors involved in basal immunity in mammals and plants interact with the same MAMP-responsive cis-sequence from Arabidopsis thaliana

WRKY and NF-κB transcription factors, involved in innate immunity in plants and mammals, interact with the same cis-sequence. Novel microbe-associated molecular pattern (MAMP)-responsive cis-sequences, designated type II WT-boxes, are required for flg22-responsive gene expression in Arabidopsis thaliana protoplasts. While type I WT-boxes like TGACTTTT and CGACTTTT interact with WRKY transcription factors (TFs), the question remained which TFs bind to the type II WT-boxes GGACTTTC, GGACTTTT, and GGACTTTG. Surprisingly, a bioinformatic analysis predicts mouse (Mus musculus) NF-κB p65 as a TF interacting with type II WT-boxes. NF-κB p65, like WRKY factors in plants, plays a role in innate immunity in mammals. Therefore, the interaction of NF-κB p65 with type II WT-boxes was tested experimentally. NF-κB p65 requires the WT-boxes GGACTTTC, GGACTTTT, and GGACTTTG for activating reporter gene expression in plant cells. NF-κB p65 directly binds to WT-box containing synthetic promoters in vitro and requires the WT-box for binding. Earlier studies indicate that the sequence GGACTTTC is also required for WRKY26 mediated reporter gene activation. Here it is shown that WRKY26, like NF-κB p65, binds to the sequence GGACTTTC. Consistent with other recent studies, type II WT boxes are WRKY binding sites and the distinction between type I and type II no longer applies.

Redundancy, Feedback, and Robustness in the Arabidopsis thaliana BZR/BEH Gene Family

Organismal development is remarkably robust, tolerating stochastic errors to produce consistent, so-called canalized adult phenotypes. The mechanistic underpinnings of developmental robustness are poorly understood, but recent studies implicate certain features of genetic networks such as functional redundancy, connectivity, and feedback. Here, we examine the BZR/BEH gene family, whose function contributes to embryonic stem development in the plant Arabidopsis thaliana, to test current assumptions on functional redundancy and trait robustness. Our analyses of BZR/BEH gene mutants and mutant combinations revealed that functional redundancy among these gene family members is not necessary for trait robustness. Connectivity is another commonly cited determinant of robustness; however, we found no correlation between connectivity among gene family members or their connectivity with other transcription factors and effects on developmental robustness. Instead, our data suggest that BEH4, the earliest diverged family member, modulates developmental robustness. We present evidence indicating that regulatory cross-talk among gene family members is integrated by BEH4 to promote wild-type levels of developmental robustness. Further, the chaperone HSP90, a known determinant of developmental robustness, appears to act via BEH4 in maintaining robustness of embryonic stem length. In summary, we demonstrate that even among closely related transcription factors, trait robustness can arise through the activity of a single gene family member, challenging common assumptions about the molecular underpinnings of robustness.

Mining Late Embryogenesis Abundant (LEA) Family Genes in Cleistogenes songorica, a Xerophyte Perennial Desert Plant

Plant growth and development depends on its ability to maintain optimal cellular homeostasis during abiotic and biotic stresses. Cleistogenes songorica, a xerophyte desert plant, is known to have novel drought stress adaptation strategies and contains rich pools of stress tolerance genes. Proteins encoded by Late Embryogenesis Abundant (LEA) family genes promote cellular activities by functioning as disordered molecules, or by limiting collisions between enzymes during stresses. To date, functions of the LEA family genes have been heavily investigated in many plant species except perennial monocotyledonous species. In this study, 44 putative LEA genes were identified in the C. songorica genome and were grouped into eight subfamilies, based on their conserved protein domains and domain organizations. Phylogenetic analyses indicated that C. songorica Dehydrin and LEA_2 subfamily proteins shared high sequence homology with stress responsive Dehydrin proteins from Arabidopsis. Additionally, promoter regions of CsLEA_2 or CsDehydrin subfamily genes were rich in G-box, drought responsive (MBS), and/or Abscisic acid responsive (ABRE) cis-regulatory elements. In addition, gene expression analyses indicated that genes from these two subfamilies were highly responsive to heat stress and ABA treatment, in both leaves and roots. In summary, the results from this study provided a comprehensive view of C. songoricaLEA genes and the potential applications of these genes for the improvement of crop tolerance to abiotic stresses.

Circadian control of abscisic acid biosynthesis and signalling pathways revealed by genome-wide analysis of LHY binding targets

The LATE ELONGATED HYPOCOTYL (LHY) transcription factor functions as part of the oscillatory mechanism of the Arabidopsis circadian clock. This paper reports the genome-wide analysis of its binding targets and reveals a role in the control of abscisic acid (ABA) biosynthesis and downstream responses. LHY directly repressed expression of 9-cis-epoxycarotenoid dioxygenase enzymes, which catalyse the rate-limiting step of ABA biosynthesis. This suggested a mechanism for the circadian control of ABA accumulation in wild-type plants. Consistent with this hypothesis, ABA accumulated rhythmically in wild-type plants, peaking in the evening. LHY-overexpressing plants had reduced levels of ABA under drought stress, whereas loss-of-function mutants exhibited an altered rhythm of ABA accumulation. LHY also bound the promoter of multiple components of ABA signalling pathways, suggesting that it may also act to regulate responses downstream of the hormone. LHY promoted expression of ABA-responsive genes responsible for increased tolerance to drought and osmotic stress but alleviated the inhibitory effect of ABA on seed germination and plant growth. This study reveals a complex interaction between the circadian clock and ABA pathways, which is likely to make an important contribution to plant performance under drought and osmotic stress conditions.

The DNA binding landscape of the maize AUXIN RESPONSE FACTOR family

AUXIN RESPONSE FACTORS (ARFs) are plant-specific transcription factors (TFs) that couple perception of the hormone auxin to gene expression programs essential to all land plants. As with many large TF families, a key question is whether individual members determine developmental specificity by binding distinct target genes. We use DAP-seq to generate genome-wide in vitro TF:DNA interaction maps for fourteen maize ARFs from the evolutionarily conserved A and B clades. Comparative analysis reveal a high degree of binding site overlap for ARFs of the same clade, but largely distinct clade A and B binding. Many sites are however co-occupied by ARFs from both clades, suggesting transcriptional coordination for many genes. Among these, we investigate known QTLs and use machine learning to predict the impact of cis-regulatory variation. Overall, large-scale comparative analysis of ARF binding suggests that auxin response specificity may be determined by factors other than individual ARF binding site selection.

LMI1 homeodomain protein regulates organ proportions by spatial modulation of endoreduplication

Here, Vuolo et al. investigated the mechanisms controlling the relative size of leaves compared with their lateral appendages (stipules). Using genetics, live imaging, and modeling, they demonstrate that the LATE MERISTEM IDENTITY1 (LMI1) homeodomain protein regulates stipule proportions via an endoreduplication-dependent trade-off that limits tissue size despite increasing cell growth.

How the interplay between cell- and tissue-level processes produces correctly proportioned organs is a key problem in biology. In plants, the relative size of leaves compared with their lateral appendages, called stipules, varies tremendously throughout development and evolution, yet relevant mechanisms remain unknown. Here we use genetics, live imaging, and modeling to show that in Arabidopsis leaves, the LATE MERISTEM IDENTITY1 (LMI1) homeodomain protein regulates stipule proportions via an endoreduplication-dependent trade-off that limits tissue size despite increasing cell growth. LM1 acts through directly activating the conserved mitosis blocker WEE1, which is sufficient to bypass the LMI1 requirement for leaf proportionality.

OsPHR3 affects the traits governing nitrogen homeostasis in rice

BACKGROUND

Phosphate (Pi) and Nitrogen (N) are essential macronutrients required for plant growth and development. In Arabidopsis thaliana (Arabidopsis), the transcription factor PHR1 acts as a Pi central regulator. PHL1 is a homolog of PHR1 and also plays a role in maintaining Pi homeostasis. In rice (Oryza sativa), OsPHR1-4 are the orthologs of PHR1 and have been implicated in regulating sensing and signaling cascades governing Pi homeostasis.

RESULTS

Here the role of OsPHR3 was examined in regulating the homeostasis of N under different Pi regimes. Deficiencies of different variants of N exerted attenuating effects on the relative expression levels of OsPHR3 in a tissue-specific manner. For the functional characterization of OsPHR3, its Tos17 insertion homozygous mutants i.e., osphr3-1, osphr3-2, and osphr3-3 were compared with the wild-type for various morphophysiological and molecular traits during vegetative (hydroponics with different regimes of N variants) and reproductive (pot soil) growth phases. During vegetative growth phase, compared with the wild-type, OsPHR3 mutants showed significant variations in the adventitious root development, influx rates of 15N-NO3- and 15N-NH4+, concentrations of total N, NO3- and NH4+ in different tissues, and the relative expression levels of OsNRT1.1a, OsNRT2.4, OsAMT1;1, OsNia1 and OsNia2. The effects of the mutation in OsPHR3 was also explicit on the seed-set and grain yield during growth in a pot soil. Although Pi deficiency affected total N and NO3- concentration, the lateral root development and the relative expression levels of some of the NO3- and NH4+ transporter genes, its availability did not exert any notable regulatory influences on the traits governing N homeostasis.

CONCLUSIONS

OsPHR3 plays a pivotal role in regulating the homeostasis of N independent of Pi availability.

ChIP-Seq reveals that QsMYB1 directly targets genes involved in lignin and suberin biosynthesis pathways in cork oak (Quercus suber)

BACKGROUND

Gene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy.

RESULTS

Results provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly.

CONCLUSION

To our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.

Coordination of Chloroplast Development through the Action of the GNC and GLK Transcription Factor Families

Fundamental questions regarding how chloroplasts develop from proplastids remain poorly understood despite their central importance to plant life. Two families of nuclear transcription factors, the GATA NITRATE-INDUCIBLE CARBON-METABOLISM-INVOLVED (GNC) and GOLDEN TWO-LIKE (GLK) families, have been implicated in directly and positively regulating chloroplast development. Here, we determined the degree of functional overlap between the two transcription factor families in Arabidopsis (Arabidopsis thaliana), characterizing their ability to regulate chloroplast biogenesis both alone and in concert. We determined the DNA-binding motifs for GNC and GLK2 using protein-binding microarrays; the enrichment of these motifs in transcriptome datasets indicates that GNC and GLK2 are repressors and activators of gene expression, respectively. ChIP-seq analysis of GNC identified PHYTOCHROME INTERACTING FACTOR and brassinosteroid activity genes as targets whose repression by GNC facilitates chloroplast biogenesis. In addition, GNC targets and represses genes involved in ERECTA signaling and thereby facilitates stomatal development. Our results define key regulatory features of the GNC and GLK transcription factor families that contribute to the control of chloroplast biogenesis and photosynthetic activity, including areas of independence and cross talk.

Genome-wide analysis of purple acid phosphatase structure and expression in ten vegetable species

Background

Acquisition of external phosphorus (P) and optimisation of internal P are essential for plant growth and development, and insufficient availability of P in soils is a major challenge in agriculture. Members of the purple acid phosphatase (PAP) family of enzymes are candidates for increasing P use efficiency. Herein, we identified PAP homologs in the genomes of 10 vegetable species, along with Arabidopsis thaliana and Amborella trichopoda as references, to provide fundamental knowledge for this family.

Results

Phylogenetic analysis of protein sequences revealed nine distinct clades, indicating that functional differentiation of extant PAPs was established prior to the emergence of early angiosperms, and conserved among homologs in each clade. Analysis of transcript abundance in different tissues (root, stem, leaf, flower, and fruit) and following phosphates (Pi) starvation treatments from published RNA-seq transcriptome datasets facilitated comprehensive evaluation of expression patterns, and some groups of tissue-specific and Pi starvation-induced PAPs were characterised. Conserved motifs identified from upstream sequences of homologs that are highly expressed in particular tissues or following starvation treatment suggests that divergence in PAP gene expression is associated with cis-acting elements in promoters.

Conclusions

The genome-wide analysis of PAP enzyme structure and transcriptional expression patterns advance our understanding of PAP family in vegetables genomes. Therefore, PAP homologs with known enzyme structures and expression profiles could serve as targets for plant breeding and/or genetic engineering programs to improve P acquisition and use.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5022-1) contains supplementary material, which is available to authorized users.

An inhibitor of apoptosis (SfIAP) interacts with SQUAMOSA promoter-binding protein (SBP) transcription factors that exhibit pro-cell death characteristics

Despite the importance of proper cell death regulation across broad evolutionary distances, an understanding of the molecular machinery underpinning this fundamental process in plants remains largely elusive. This is despite its critical importance to development, homeostasis, and proper responses to stress. The identification of endogenous plant regulators of cell death has been hindered by the fact that many core regulators of cell death in animals are absent in plant genomes. Remarkably, numerous studies have shown that the ectopic expression of animal prosurvival genes in plants can suppress cell death imposed by many stresses. In this study, we capitalize on the ectopic expression of one of these animal prosurvival genes, an inhibitor of apoptosis from Spodoptera frugiperda (SfIAP), to identify novel cell death regulators in plants. A yeast two-hybrid assay was conducted using SfIAP as bait to screen a tomato cDNA library. This screen identified several transcription factors of the SQUAMOSA promoter-binding protein (SBP) family as potential SfIAP binding partners. We confirmed this interaction in vivo for our top two interactors, SlySBP8b and SlySBP12a, using coimmunoprecipitation. Interestingly, overexpression of SlySBP8b and SlySBP12a induced cell death in Nicotiana benthamiana leaves. Overexpression of these two transcription factors also induced the accumulation of reactive oxygen species and enhanced the growth of the necrotrophic pathogen Alternaria alternata. Fluorescence microscopy confirmed the nuclear localization of both SlySBP8b and SlySBP12a, while SlySBP12a was also localized to the ER membrane. These results suggest a prodeath role for SlySBP8b and SlySBP12a and implicate ER membrane tethering as a means of regulating SlySBP12a activity.

Pioneer Factors in Animals and Plants-Colonizing Chromatin for Gene Regulation

Unlike most transcription factors (TF), pioneer TFs have a specialized role in binding closed regions of chromatin and initiating the subsequent opening of these regions. Thus, pioneer TFs are key factors in gene regulation with critical roles in developmental transitions, including organ biogenesis, tissue development, and cellular differentiation. These developmental events involve some major reprogramming of gene expression patterns, specifically the opening and closing of distinct chromatin regions. Here, we discuss how pioneer TFs are identified using biochemical and genome-wide techniques. What is known about pioneer TFs from animals and plants is reviewed, with a focus on the strategies used by pioneer factors in different organisms. Finally, the different molecular mechanisms pioneer factors used are discussed, highlighting the roles that tertiary and quaternary structures play in nucleosome-compatible DNA-binding.

Expression of the Intracellular COPT3-Mediated Cu Transport Is Temporally Regulated by the TCP16 Transcription Factor

Copper is an essential element in plants. When scarce, copper is acquired from extracellular environment or remobilized from intracellular sites, through members of the high affinity copper transporters family COPT located at the plasma membrane and internal membrane, respectively. Here, we show that COPT3 is an intracellular copper transporter, located at a compartment of the secretory pathway, that is mainly expressed in pollen grains and vascular bundles. Contrary to the COPT1 plasma membrane member, the expression of the internal COPT3 membrane transporter was higher at 12 h than at 0 h of a neutral photoperiod day under copper deficiency. The screening of a library of conditionally overexpressed transcription factors implicated members of the TCP family in the COPT3 differential temporal expression pattern. Particularly, in vitro, TCP16 was found to bind to the COPT3 promoter and down-regulated its expression. Accordingly, TCP16 was mainly expressed at 0 h under copper deficiency and induced at 12 h by copper excess. Moreover, TCP16 overexpression resulted in increased sensitivity to copper deficiency, whereas the tcp16 mutant was sensitive to copper excess. Both copper content and the expression of particular copper status markers were altered in plants with modified levels of TCP16. Consistent with TCP16 affecting pollen development, the lack of COPT3 function led to altered pollen morphology. Furthermore, analysis of copt3 and COPT3 overexpressing plants revealed that COPT3 function exerted a negative effect on TCP16 expression. Taken together, these results suggest a differential daily regulation of copper uptake depending on the external and internal copper pools, in which TCP16 inhibits copper remobilization at dawn through repression of intracellular transporters.

Regulatory Divergence in Wound-Responsive Gene Expression between Domesticated and Wild Tomato

The evolution of transcriptional regulatory mechanisms is central to how stress response and tolerance differ between species. However, it remains largely unknown how divergence in cis-regulatory sites and, subsequently, transcription factor (TF) binding specificity contribute to stress-responsive expression divergence, particularly between wild and domesticated species. By profiling wound-responsive gene transcriptomes in wild Solanum pennellii and domesticated S. lycopersicum, we found extensive wound response divergence and identified 493 S. lycopersicum and 278 S. pennellii putative cis-regulatory elements (pCREs) that were predictive of wound-responsive gene expression. Only 24-52% of these wound response pCREs (depending on wound response patterns) were consistently enriched in the putative promoter regions of wound-responsive genes across species. In addition, between these two species, their differences in pCRE site sequences were significantly and positively correlated with differences in wound-responsive gene expression. Furthermore, ∼11-39% of pCREs were specific to only one of the species and likely bound by TFs from different families. These findings indicate substantial regulatory divergence in these two plant species that diverged ∼3-7 million years ago. Our study provides insights into the mechanistic basis of how the transcriptional response to wounding is regulated and, importantly, the contribution of cis-regulatory components to variation in wound-responsive gene expression between a wild and a domesticated plant species.

Evolutionary history of HOMEODOMAIN LEUCINE ZIPPER transcription factors during plant transition to land

Plant transition to land required several regulatory adaptations. The mechanisms behind these changes remain unknown. Since the evolution of transcription factors (TFs) families accompanied this transition, we studied the HOMEODOMAIN LEUCINE ZIPPER (HDZ) TF family known to control key developmental and environmental responses. We performed a phylogenetic and bioinformatics analysis of HDZ genes using transcriptomic and genomic datasets from a wide range of Viridiplantae species. We found evidence for the existence of HDZ genes in chlorophytes and early-divergent charophytes identifying several HDZ members belonging to the four known classes (I-IV). Furthermore, we inferred a progressive incorporation of auxiliary motifs. Interestingly, most of the structural features were already present in ancient lineages. Our phylogenetic analysis inferred that the origin of classes I, III, and IV is monophyletic in land plants in respect to charophytes. However, class IIHDZ genes have two conserved lineages in charophytes and mosses that differ in the CPSCE motif. Our results indicate that the HDZ family was already present in green algae. Later, the HDZ family expanded accompanying critical plant traits. Once on land, the HDZ family experienced multiple duplication events that promoted fundamental neo- and subfunctionalizations for terrestrial life.

Comparative analysis of binding patterns of MADS-domain proteins in Arabidopsis thaliana

BACKGROUND

Correct flower formation requires highly specific temporal and spatial regulation of gene expression. In Arabidopsis thaliana the majority of the master regulators that determine flower organ identity belong to the MADS-domain transcription factor family. The canonical DNA binding motif for this transcription factor family is the CArG-box, which has the consensus CC(A/T)₆GG. However, so far, a comprehensive analysis of MADS-domain binding patterns has not yet been performed.

RESULTS

Eight publicly available ChIP-seq datasets of MADS-domain proteins that regulate the floral transition and flower formation were analyzed. Surprisingly, the preferred DNA binding motif of each protein was a CArG-box with an NAA extension. Furthermore, motifs of other transcription factors were found in the vicinity of binding sites of MADS-domain transcription factors, suggesting that interaction of MADS-domain proteins with other transcription factors is important for target gene regulation. Finally, conservation of CArG-boxes between Arabidopsis ecotypes was assessed to obtain information about their evolutionary importance. CArG-boxes that fully matched the consensus were more conserved than other CArG-boxes, suggesting that the perfect CArG-box is evolutionary more important than other CArG-box variants.

CONCLUSION

Our analysis provides detailed insight into MADS-domain protein binding patterns. The results underline the importance of an extended version of the CArG-box and provide a first view on evolutionary conservation of MADS-domain protein binding sites in Arabidopsis ecotypes.

Distinct tissue-specific transcriptional regulation revealed by gene regulatory networks in maize

BACKGROUND

Transcription factors (TFs) are proteins that can bind to DNA sequences and regulate gene expression. Many TFs are master regulators in cells that contribute to tissue-specific and cell-type-specific gene expression patterns in eukaryotes. Maize has been a model organism for over one hundred years, but little is known about its tissue-specific gene regulation through TFs. In this study, we used a network approach to elucidate gene regulatory networks (GRNs) in four tissues (leaf, root, SAM and seed) in maize. We utilized GENIE3, a machine-learning algorithm combined with large quantity of RNA-Seq expression data to construct four tissue-specific GRNs. Unlike some other techniques, this approach is not limited by high-quality Position Weighed Matrix (PWM), and can therefore predict GRNs for over 2000 TFs in maize.

RESULTS

Although many TFs were expressed across multiple tissues, a multi-tiered analysis predicted tissue-specific regulatory functions for many transcription factors. Some well-studied TFs emerged within the four tissue-specific GRNs, and the GRN predictions matched expectations based upon published results for many of these examples. Our GRNs were also validated by ChIP-Seq datasets (KN1, FEA4 and O2). Key TFs were identified for each tissue and matched expectations for key regulators in each tissue, including GO enrichment and identity with known regulatory factors for that tissue. We also found functional modules in each network by clustering analysis with the MCL algorithm.

CONCLUSIONS

By combining publicly available genome-wide expression data and network analysis, we can uncover GRNs at tissue-level resolution in maize. Since ChIP-Seq and PWMs are still limited in several model organisms, our study provides a uniform platform that can be adapted to any species with genome-wide expression data to construct GRNs. We also present a publicly available database, maize tissue-specific GRN (mGRN, https://www.bio.fsu.edu/mcginnislab/mgrn/ ), for easy querying. All source code and data are available at Github ( https://github.com/timedreamer/maize_tissue-specific_GRN ).

Arabidopsis SWC4 Binds DNA and Recruits the SWR1 Complex to Modulate Histone H2A.Z Deposition at Key Regulatory Genes

Deposition of the H2A.Z histone variant by the SWR1 complex (SWR1-C) in regulatory regions of specific loci modulates transcription. Characterization of mutations in Arabidopsis thaliana homologs of yeast SWR1-C has revealed a role for H2A.Z exchange in a variety of developmental processes. Nevertheless, the exact composition of plant SWR1-C and how it is recruited to target genes remains to be established. Here we show that SWC4, the Arabidopsis homolog of yeast SANT domain protein Swc4/Eaf2, is a DNA-binding protein that interacts with SWR1-C subunits. We demonstrate that the swc4-1 knockout mutant is embryo-lethal, while SWC4 RNAi knockdown lines display pleiotropic phenotypic alterations in vegetative and reproductive traits, including acceleration of flowering time, indicating that SWC4 controls post-embryonic processes. Transcriptomic analyses and genome-wide profiling of H2A.Z indicate that SWC4 represses transcription of a number of genes, including the floral integrator FT and key transcription factors, mainly by modulating H2A.Z deposition. Interestingly, SWC4 silencing does not affect H2A.Z deposition at the FLC locus nor expression of this gene, a master regulator of flowering previously shown to be controlled by SWR1-C. Importantly, we find that SWC4 recognizes specific AT-rich DNA elements in the chromatin regions of target genes and that SWC4 silencing impairs SWR1-C binding at FT. Collectively, our data suggest that SWC4 regulates plant growth and development by aiding SWR1-C recruitment and modulating H2A.Z deposition.

The Defense Phytohormone Signaling Network Enables Rapid, High-Amplitude Transcriptional Reprogramming during Effector-Triggered Immunity

The phytohormone network consisting of jasmonate, ethylene, PHYTOALEXIN-DEFICIENT4, and salicylic acid signaling is required for the two modes of plant immunity, pattern-triggered immunity (PTI), and effector-triggered immunity (ETI). A previous study showed that during PTI, the transcriptional responses of over 5000 genes qualitatively depend on complex interactions between the network components. However, the role of the network in transcriptional reprogramming during ETI and whether it differs between PTI and ETI remain elusive. Here, we generated time-series RNA-sequencing data of Arabidopsis thaliana wild-type and combinatorial mutant plants deficient in components of the network upon challenge with virulent or ETI-triggering avirulent strains of the foliar bacterial pathogen Pseudomonas syringae Resistant plants such as the wild type achieved high-amplitude transcriptional reprogramming 4 h after challenge with avirulent strains and sustained this transcriptome response. Strikingly, susceptible plants including the quadruple network mutant showed almost identical transcriptome responses to resistant plants but with several hours delay. Furthermore, gene coexpression network structure was highly conserved between the wild type and quadruple mutant. Thus, in contrast to PTI, the phytohormone network is required only for achieving high-amplitude transcriptional reprogramming within the early time window of ETI against this bacterial pathogen.

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.

Genome-Wide Screening and Characterization of the Dof Gene Family in Physic Nut (Jatropha curcas L.)

Physic nut (Jatropha curcas L.) is a species of flowering plant with great potential for biofuel production and as an emerging model organism for functional genomic analysis, particularly in the Euphorbiaceae family. DNA binding with one finger (Dof) transcription factors play critical roles in numerous biological processes in plants. Nevertheless, the knowledge about members, and the evolutionary and functional characteristics of the Dof gene family in physic nut is insufficient. Therefore, we performed a genome-wide screening and characterization of the Dof gene family within the physic nut draft genome. In total, 24 JcDof genes (encoding 33 JcDof proteins) were identified. All the JcDof genes were divided into three major groups based on phylogenetic inference, which was further validated by the subsequent gene structure and motif analysis. Genome comparison revealed that segmental duplication may have played crucial roles in the expansion of the JcDof gene family, and gene expansion was mainly subjected to positive selection. The expression profile demonstrated the broad involvement of JcDof genes in response to various abiotic stresses, hormonal treatments and functional divergence. This study provides valuable information for better understanding the evolution of JcDof genes, and lays a foundation for future functional exploration of JcDof genes.

Delineation of condition specific Cis- and Trans-acting elements in plant promoters under various Endo- and exogenous stimuli

BACKGROUND

Transcription factors (TFs) play essential roles during plant development and response to environmental stresses. However, the relationships among transcription factors, cis-acting elements and target gene expression under endo- and exogenous stimuli have not been systematically characterized.

RESULTS

Here, we developed a series of bioinformatics analysis methods to infer transcriptional regulation by using numerous gene expression data from abiotic stresses and hormones treatments. After filtering the expression profiles of TF-encoding genes, 291 condition specific transcription factors (CsTFs) were obtained. Differentially expressed genes were then classified into various co-expressed gene groups based on each CsTFs. In the case studies of heat stress and ABA treatment, several known and novel cis-acting elements were identified following our bioinformatics approach. Significantly, a palindromic sequence of heat-responsive elements is recognized, and also obtained from a 3D protein structure of heat-shock protein-DNA complex. Notably, overrepresented 3- and 4-mer motifs in an enriched 8-mer motif could be a core cis-element for a CsTF. In addition, the results suggest DNA binding preferences of the same CsTFs are different according to various conditions.

CONCLUSIONS

The overall results illustrate this study may be useful in identifying condition specific cis- and trans- regulatory elements and facilitate our understanding of the relationships among TFs, cis-acting elements and target gene expression.

Transcriptome analysis of pecan seeds at different developing stages and identification of key genes involved in lipid metabolism

Pecan is an economically important nut crop tree due to its unique texture and flavor properties. The pecan seed is rich of unsaturated fatty acid and protein. However, little is known about the molecular mechanisms of the biosynthesis of fatty acids in the developing seeds. In this study, transcriptome sequencing of the developing seeds was performed using Illumina sequencing technology. Pecan seed embryos at different developmental stages were collected and sequenced. The transcriptomes of pecan seeds at two key developing stages (PA, the initial stage and PS, the fast oil accumulation stage) were also compared. A total of 82,155 unigenes, with an average length of 1,198 bp from seven independent libraries were generated. After functional annotations, we detected approximately 55,854 CDS, among which, 2,807 were Transcription Factor (TF) coding unigenes. Further, there were 13,325 unigenes that showed a 2-fold or greater expression difference between the two groups of libraries (two developmental stages). After transcriptome analysis, we identified abundant unigenes that could be involved in fatty acid biosynthesis, degradation and some other aspects of seed development in pecan. This study presents a comprehensive dataset of transcriptomic changes during the seed development of pecan. It provides insights in understanding the molecular mechanisms responsible for fatty acid biosynthesis in the seed development. The identification of functional genes will also be useful for the molecular breeding work of pecan.

MYB52 Negatively Regulates Pectin Demethylesterification in Seed Coat Mucilage

Pectin, which is a major component of the plant primary cell walls, is synthesized and methyl-esterified in the Golgi apparatus and then demethylesterified by pectin methylesterases (PMEs) located in the cell wall. The degree of methylesterification affects the functional properties of pectin, and thereby influences plant growth, development and defense. However, little is known about the mechanisms that regulate pectin demethylesterification. Here, we show that in Arabidopsis (Arabidopsis thaliana) seed coat mucilage, the absence of the MYB52 transcription factor is correlated with an increase in PME activity and a decrease in the degree of pectin methylesterification. Decreased methylesterification in the myb52 mutant is also correlated with an increase in the calcium content of the seed mucilage. Chromatin immunoprecipitation analysis and molecular genetic studies suggest that MYB52 transcriptionally activates PECTIN METHYLESTERASE INHIBITOR6 (PMEI6), PMEI14, and SUBTILISIN-LIKE SER PROTEASE1.7 (SBT1.7) by binding to their promoters. PMEI6 and SBT1.7 have previously been shown to be involved in seed coat mucilage demethylesterification. Our characterization of two PMEI14 mutants suggests that PMEI14 has a role in seed coat mucilage demethylesterification, although its activity may be confined to the seed coat in contrast to PMEI6, which functions in the whole seed. Our demonstration that MYB52 negatively regulates pectin demethylesterification in seed coat mucilage, and the identification of components of the molecular network involved, provides new insight into the regulatory mechanism controlling pectin demethylesterification and increases our understanding of the transcriptional regulation network involved in seed coat mucilage formation.

Repression of Nitrogen Starvation Responses by Members of the Arabidopsis GARP-Type Transcription Factor NIGT1/HRS1 Subfamily

Nitrogen (N) is often a limiting nutrient whose availability determines plant growth and productivity. Because its availability is often low and/or not uniform over time and space in nature, plants respond to variations in N availability by altering uptake and recycling mechanisms, but the molecular mechanisms underlying how these responses are regulated are poorly understood. Here, we show that a group of GARP G2-like transcription factors, Arabidopsis thaliana NITRATE-INDUCIBLE, GARP-TYPE TRANSCRIPTIONAL REPRESSOR1/HYPERSENSITIVE TO LOW Pi-ELICITED PRIMARY ROOT SHORTENING1 proteins (NIGT1/HRS1s), are factors that bind to the promoter of the N starvation marker NRT2.4 and repress an array of N starvation-responsive genes under conditions of high N availability. Transient assays and expression analysis demonstrated that NIGT1/HRS1s are transcriptional repressors whose expression is regulated by N availability. We identified target genes of the NIGT1/HRS1s by genome-wide transcriptome analyses and found that they are significantly enriched in N starvation response-related genes, including N acquisition, recycling, remobilization, and signaling genes. Loss of NIGT1/HRS1s resulted in deregulation of N acquisition and accumulation. We propose that NIGT1/HRS1s are major regulators of N starvation responses that play an important role in optimizing N acquisition and utilization under fluctuating N conditions.

Predicting Transcription Factor Binding Sites and Their Cognate Transcription Factors Using Gene Expression Data

A transcription factor (TF) is a DNA binding protein that targets specific binding-sites (TFBSs) to regulate the transcript levels of its downstream genes. Thus, identifying the TF-TFBS pairs is a pivotal step in understanding the function of TFs and the regulatory network in an organism. Here, we describe two methods for predicting the TFBS of a given TF and for predicting the cognate TF of a given TFBS from a set of strongly co-expressed genes, using time-course transcriptome data of maize developing leaves.

The intracellular immune receptor Rx1 regulates the DNA-binding activity of a Golden2-like transcription factor

Plant nucleotide-binding leucine-rich repeat (NLR) proteins enable the immune system to recognize and respond to pathogen attack. An early consequence of immune activation is transcriptional reprogramming, and some NLRs have been shown to act in the nucleus and interact with transcription factors. The Rx1 NLR protein of potato is further able to bind and distort double-stranded DNA. However, Rx1 host targets that support a role for Rx1 in transcriptional reprogramming at DNA are unknown. Here, we report a functional interaction between Rx1 and NbGlk1, a Golden2-like transcription factor. Rx1 binds to NbGlk1 in vitro and in planta. NbGlk1 binds to known Golden2-like consensus DNA sequences. Rx1 reduces the binding affinity of NbGlk1 for DNA in vitro. NbGlk1 activates cellular responses to potato virus X, whereas Rx1 associates with NbGlk1 and prevents its assembly on DNA in planta unless activated by PVX. This study provides new mechanistic insight into how an NLR can coordinate an immune signaling response at DNA following pathogen perceptions.

A molecular framework controlling style morphology in Brassicaceae

Organ formation in multicellular organisms depends on the coordinated activities of regulatory components that integrate developmental and hormonal cues to control gene expression and mediate cell-type specification. For example, development of the Arabidopsis gynoecium is tightly controlled by distribution and synthesis of the plant hormone auxin. The functions of several transcription factors (TFs) have been linked with auxin dynamics during gynoecium development; yet how their activities are coordinated is not known. Here, we show that five such TFs function together to ensure polarity establishment at the gynoecium apex. The auxin response factor ETTIN (ARF3; herein, ETT) is a central component of this framework. Interaction of ETT with TF partners is sensitive to the presence of auxin and our results suggest that ETT forms part of a repressive gene-regulatory complex. We show that this function is conserved between members of the Brassicaceae family and that variation in an ETT subdomain affects interaction strengths and gynoecium morphology. These results suggest that variation in affinities between conserved TFs can lead to morphological differences and thus contribute to the evolution of diverse organ shapes.

Summary: Variation in interaction affinity between transcription factors of an ETTIN-containing complex underlies diversity of gynoecium style structure among members of the Brassicacea family.

Loss of ZNF32 augments the regeneration of nervous lateral line system through negative regulation of SOX2 transcription

Human zinc finger protein 32 (ZNF32) is a Cys2-His2 zinc-finger transcription factor that plays an important role in cell fate, yet much of its function remains unknown. Here, we reveal that the zebrafish ZNF32 homologue zfZNF32 is expressed in the nervous system, particularly in the lateral line system. ZfZNF32 knock-out zebrafish (zfZNF^−/−) were generated using the CRISPR-associated protein 9 system. We found that the regenerative capacity of the lateral line system was increased in zfZNF^−/− upon hair cell damage compared with the wild type. Moreover, SOX2 was essential for the zfZNF32-dependent modulation of lateral line system regeneration. Mechanistic studies showed that ZNF32 suppressed SOX2 transcription by directly binding to a consensus sequence (5′-gcattt-32) in the SOX2 promoter. In addition, ZNF32 localizes to the nucleus, and we have identified that amino acids 1-169 (Aa 1-169) and each of three independent nuclear localization signals (NLSs) in ZNF32 are indispensable for ZNF32 nuclear trafficking. Mutating the NLSs disrupted the inhibitory effect of ZNF32 in SOX2 expression, highlighting the critical role of the NLSs in ZNF32 function. Our findings reveal a pivotal role for ZNF32 function in SOX2 expression and regeneration regulation.

Ancient duons may underpin spatial patterning of gene expression in C4 leaves

If the highly efficient C4 photosynthesis pathway could be transferred to crops with the C3 pathway there could be yield gains of up to 50%. It has been proposed that the multiple metabolic and developmental modifications associated with C4 photosynthesis are underpinned by relatively few master regulators that have allowed the evolution of C4 photosynthesis more than 60 times in flowering plants. Here we identify a component of one such regulator that consists of a pair of cis-elements located in coding sequence of multiple genes that are preferentially expressed in bundle sheath cells of C4 leaves. These motifs represent duons as they play a dual role in coding for amino acids as well as controlling the spatial patterning of gene expression associated with the C4 leaf. They act to repress transcription of C4 photosynthesis genes in mesophyll cells. These duons are also present in the C3 model Arabidopsis thaliana, and, in fact, are conserved in all land plants and even some algae that use C3 photosynthesis. C4 photosynthesis therefore appears to have coopted an ancient regulatory code to generate the spatial patterning of gene expression that is a hallmark of C4 photosynthesis. This intragenic transcriptional regulatory sequence could be exploited in the engineering of efficient photosynthesis of crops.

Co-expression network analysis and cis-regulatory element enrichment determine putative functions and regulatory mechanisms of grapevine ATL E3 ubiquitin ligases

Arabidopsis thaliana Toxicos en Levadura (ATL) proteins are a subclass of the RING-H2 zinc finger binding E3 ubiquitin ligases. The grapevine (Vitis vinifera) ATL family was recently characterized, revealing 96 members that are likely to be involved in several physiological processes through protein ubiquitination. However, the final targets and biological functions of most ATL E3 ligases are still unknown. We analyzed the co-expression networks among grapevine ATL genes across a set of transcriptomic data related to defense and abiotic stress, combined with a condition-independent dataset. This revealed strong correlations between ATL proteins and diverse signal transduction components and transcriptional regulators, in particular those involved in immunity. An enrichment analysis of cis-regulatory elements in ATL gene promoters and related co-expressed genes highlighted the importance of hormones in the regulation of ATL gene expression. Our work identified several ATL proteins as candidates for further studies aiming to decipher specific grapevine resistance mechanisms activated in response to pathogens.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Genome-wide analyses of direct target genes of four rice NAC-domain transcription factors involved in drought tolerance

BACKGROUND

Plant stress responses and mechanisms determining tolerance are controlled by diverse sets of genes. Transcription factors (TFs) have been implicated in conferring drought tolerance under drought stress conditions, and the identification of their target genes can elucidate molecular regulatory networks that orchestrate tolerance mechanisms.

RESULTS

We generated transgenic rice plants overexpressing the 4 rice TFs, OsNAC5, 6, 9, and 10, under the control of the root-specific RCc3 promoter. We showed that they were tolerant to drought stress with reduced loss of grain yield under drought conditions compared with wild type plants. To understand the molecular mechanisms underlying this tolerance, we here performed chromatin immunoprecipitation (ChIP)-Seq and RNA-Seq analyses to identify the direct target genes of the OsNAC proteins using the RCc3:6MYC-OsNAC expressing roots. A total of 475 binding loci for the 4 OsNAC proteins were identified by cross-referencing their binding to promoter regions and the expression levels of the corresponding genes. The binding loci were distributed among the promoter regions of 391 target genes that were directly up-regulated by one of the OsNAC proteins in four RCc3:6MYC-OsNAC transgenic lines. Based on gene ontology (GO) analysis, the direct target genes were related to transmembrane/transporter activity, vesicle, plant hormones, carbohydrate metabolism, and TFs. The direct targets of each OsNAC range from 4.0-8.7% of the total number of up-regulated genes found in the RNA-Seq data sets. Thus, each OsNAC up-regulates a set of direct target genes that alter root system architecture in the RCc3:OsNAC plants to confer drought tolerance. Our results provide a valuable resource for functional dissection of the molecular mechanisms of drought tolerance.

CONCLUSIONS

Many of the target genes, including transmembrane/transporter, vesicle related, auxin/hormone related, carbohydrate metabolic processes, and transcription factor genes, that are up-regulated by OsNACs act as the cellular components which would alter the root architectures of RCc3:OsNACs for drought tolerance.

Auxin Response Factors: output control in auxin biology

The phytohormone auxin is involved in almost all developmental processes in land plants. Most, if not all, of these processes are mediated by changes in gene expression. Auxin acts on gene expression through a short nuclear pathway that converges upon the activation of a family of DNA-binding transcription factors. These AUXIN RESPONSE FACTORS (ARFs) are thus the effector of auxin response and translate the chemical signal into the regulation of a defined set of genes. Given the limited number of dedicated components in auxin signaling, distinct properties among the ARF family probably contribute to the establishment of multiple unique auxin responses in plant development. In the two decades following the identification of the first ARF in Arabidopsis, much has been learnt about how these transcription factors act, and how they generate unique auxin responses. Progress in genetics, biochemistry, genomics, and structural biology has helped to develop mechanistic models for ARF action. However, despite intensive efforts, many central questions are yet to be addressed. In this review, we highlight what has been learnt about ARF transcription factors, and identify outstanding questions and challenges for the near future.

Diversity of cis-regulatory elements associated with auxin response in Arabidopsis thaliana

The phytohormone auxin regulates virtually every developmental process in land plants. This regulation is mediated via de-repression of DNA-binding auxin response factors (ARFs). ARFs bind TGTC-containing auxin response cis-elements (AuxREs), but there is growing evidence that additional cis-elements occur in auxin-responsive regulatory regions. The repertoire of auxin-related cis-elements and their involvement in different modes of auxin response are not yet known. Here we analyze the enrichment of nucleotide hexamers in upstream regions of auxin-responsive genes associated with auxin up- or down-regulation, with early or late response, ARF-binding domains, and with different chromatin states. Intriguingly, hexamers potentially bound by basic helix-loop-helix (bHLH) and basic leucine zipper (bZIP) factors as well as a family of A/T-rich hexamers are more highly enriched in auxin-responsive regions than canonical TGTC-containing AuxREs. We classify and annotate the whole spectrum of enriched hexamers and discuss their patterns of enrichment related to different modes of auxin response.

The Role of Tomato WRKY Genes in Plant Responses to Combined Abiotic and Biotic Stresses

In the field, plants constantly face a plethora of abiotic and biotic stresses that can impart detrimental effects on plants. In response to multiple stresses, plants can rapidly reprogram their transcriptome through a tightly regulated and highly dynamic regulatory network where WRKY transcription factors can act as activators or repressors. WRKY transcription factors have diverse biological functions in plants, but most notably are key players in plant responses to biotic and abiotic stresses. In tomato there are 83 WRKY genes identified. Here we review recent progress on functions of these tomato WRKY genes and their homologs in other plant species, such as Arabidopsis and rice, with a special focus on their involvement in responses to abiotic and biotic stresses. In particular, we highlight WRKY genes that play a role in plant responses to a combination of abiotic and biotic stresses.

Chromatin Immunoprecipitation Sequencing (ChIP-Seq) for Transcription Factors and Chromatin Factors in Arabidopsis thaliana Roots: From Material Collection to Data Analysis

Chromatin immunoprecipitation combined with next-generation sequencing (ChIP-seq) is a powerful technique to investigate in vivo transcription factor (TF) binding to DNA, as well as chromatin marks. Here we provide a detailed protocol for all the key steps to perform ChIP-seq in Arabidopsis thaliana roots, also working on other A. thaliana tissues and in most non-ligneous plants. We detail all steps from material collection, fixation, chromatin preparation, immunoprecipitation, library preparation, and finally computational analysis based on a combination of publicly available tools.

Using CRISPR/Cas9 System to Introduce Targeted Mutation in Arabidopsis

Gene editing enables the modification of genetic information in a targeted manner. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a gene editing approach that has been widely used in both the animal and plant studies. Plant transcription factors play critical roles in plant development. Recent studies have used the CRISPER/Cas9 system to disrupt transcription factors and cis-elements in plants. In this chapter, we describe a protocol of using the CRISPR/Cas9 system to generate transcription factor mutants in the model plant species Arabidopsis thaliana. Our protocol contains five different steps: (1) adaptor primers design, (2) cloning of sgRNA expression cassettes, (3) Golden Gate cloning, (4) E. coli and Agrobacterium transformation, and (5) Arabidopsis transformation.