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

The 'dark matter' in the plant genomes: non-coding and unannotated DNA sequences associated with open chromatin

Sequencing of complete plant genomes has become increasingly more routine since the advent of the next-generation sequencing technology. Identification and annotation of large amounts of noncoding but functional DNA sequences, including cis-regulatory DNA elements (CREs), have become a new frontier in plant genome research. Genomic regions containing active CREs bound to regulatory proteins are hypersensitive to DNase I digestion and are called DNase I hypersensitive sites (DHSs). Several recent DHS studies in plants illustrate that DHS datasets produced by DNase I digestion followed by next-generation sequencing (DNase-seq) are highly valuable for the identification and characterization of CREs associated with plant development and responses to environmental cues. DHS-based genomic profiling has opened a door to identify and annotate the 'dark matter' in sequenced plant genomes.

A modular analysis of the auxin signalling network

Auxin is essential for plant development from embryogenesis onwards. Auxin acts in large part through regulation of transcription. The proteins acting in the signalling pathway regulating transcription downstream of auxin have been identified as well as the interactions between these proteins, thus identifying the topology of this network implicating 54 Auxin Response Factor (ARF) and Aux/IAA (IAA) transcriptional regulators. Here, we study the auxin signalling pathway by means of mathematical modeling at the single cell level. We proceed analytically, by considering the role played by five functional modules into which the auxin pathway can be decomposed: the sequestration of ARF by IAA, the transcriptional repression by IAA, the dimer formation amongst ARFs and IAAs, the feedback loop on IAA and the auxin induced degradation of IAA proteins. Focusing on these modules allows assessing their function within the dynamics of auxin signalling. One key outcome of this analysis is that there are both specific and overlapping functions between all the major modules of the signaling pathway. This suggests a combinatorial function of the modules in optimizing the speed and amplitude of auxin-induced transcription. Our work allows identifying potential functions for homo- and hetero-dimerization of transcriptional regulators, with ARF:IAA, IAA:IAA and ARF:ARF dimerization respectively controlling the amplitude, speed and sensitivity of the response and a synergistic effect of the interaction of IAA with transcriptional repressors on these characteristics of the signaling pathway. Finally, we also suggest experiments which might allow disentangling the structure of the auxin signaling pathway and analysing further its function in plants.

AtNIGT1/HRS1 integrates nitrate and phosphate signals at the Arabidopsis root tip

Nitrogen and phosphorus are among the most widely used fertilizers worldwide. Nitrate (NO3(-)) and phosphate (PO4(3-)) are also signalling molecules whose respective transduction pathways are being intensively studied. However, plants are continuously challenged with combined nutritional deficiencies, yet very little is known about how these signalling pathways are integrated. Here we report the identification of a highly NO3(-)-inducible NRT1.1-controlled GARP transcription factor, HRS1, document its genome-wide transcriptional targets, and validate its cis-regulatory elements. We demonstrate that this transcription factor and a close homologue repress the primary root growth in response to P deficiency conditions, but only when NO3(-) is present. This system defines a molecular logic gate integrating P and N signals. We propose that NO3(-) and P signalling converge via double transcriptional and post-transcriptional control of the same protein, HRS1.

AtNIGT1/HRS1 integrates nitrate and phosphate signals at the Arabidopsis root tip

Nitrogen and phosphorus are among the most widely used fertilizers worldwide. Nitrate (NO₃⁻) and phosphate (PO₄³⁻) are also signaling molecules whose respective transduction pathways are being intensively studied. However, plants are continuously challenged with combined nutritional deficiencies, yet very little is known about how these signaling pathways are integrated. Here we report the identification of a highly NO₃⁻-inducible NRT1.1-controlled GARP transcription factor, HRS1, document its genome-wide transcriptional targets, and validate its cis-regulatory-elements. We demonstrate that this transcription factor and a close homolog repress primary root growth in response to P deficiency conditions, but only when NO₃⁻ is present. This system defines a molecular logic gate integrating P and N signals. We propose that NO₃⁻ and P signaling converge via double transcriptional and post-transcriptional control of the same protein, HRS1

Summarizing and exploring data of a decade of cytokinin-related transcriptomics

The genome-wide transcriptional response of the model organism Arabidopsis thaliana to cytokinin has been investigated by different research groups as soon as large-scale transcriptomic techniques became affordable. Over the last 10 years many transcriptomic datasets related to cytokinin have been generated using different technological platforms, some of which are published only in databases, culminating in an RNA sequencing experiment. Two approaches have been made to establish a core set of cytokinin-regulated transcripts by meta-analysis of these datasets using different preferences regarding their selection. Here we add another meta-analysis derived from an independent microarray platform (CATMA), combine all the meta-analyses available with RNAseq data in order to establish an advanced core set of cytokinin-regulated transcripts, and compare the results with the regulation of orthologous rice genes by cytokinin. We discuss the functions of some of the less known cytokinin-regulated genes indicating areas deserving further research to explore cytokinin function. Finally, we investigate the promoters of the core set of cytokinin-induced genes for the abundance and distribution of known cytokinin-responsive cis elements and identify a set of novel candidate motifs.

Evolutionary meandering of intermolecular interactions along the drift barrier

Many cellular functions depend on highly specific intermolecular interactions, for example transcription factors and their DNA binding sites, microRNAs and their RNA binding sites, the interfaces between heterodimeric protein molecules, the stems in RNA molecules, and kinases and their response regulators in signal-transduction systems. Despite the need for complementarity between interacting partners, such pairwise systems seem to be capable of high levels of evolutionary divergence, even when subject to strong selection. Such behavior is a consequence of the diminishing advantages of increasing binding affinity between partners, the multiplicity of evolutionary pathways between selectively equivalent alternatives, and the stochastic nature of evolutionary processes. Because mutation pressure toward reduced affinity conflicts with selective pressure for greater interaction, situations can arise in which the expected distribution of the degree of matching between interacting partners is bimodal, even in the face of constant selection. Although biomolecules with larger numbers of interacting partners are subject to increased levels of evolutionary conservation, their more numerous partners need not converge on a single sequence motif or be increasingly constrained in more complex systems. These results suggest that most phylogenetic differences in the sequences of binding interfaces are not the result of adaptive fine tuning but a simple consequence of random genetic drift.

Next-generation sequencing of genomic DNA fragments bound to a transcription factor in vitro reveals its regulatory potential

Several transcription factors (TFs) coordinate to regulate expression of specific genes at the transcriptional level. In Arabidopsis thaliana it is estimated that approximately 10% of all genes encode TFs or TF-like proteins. It is important to identify target genes that are directly regulated by TFs in order to understand the complete picture of a plant’s transcriptome profile. Here, we investigate the role of the LONG HYPOCOTYL5 (HY5) transcription factor that acts as a regulator of photomorphogenesis. We used an in vitro genomic DNA binding assay coupled with immunoprecipitation and next-generation sequencing (gDB-seq) instead of the in vivo chromatin immunoprecipitation (ChIP)-based methods. The results demonstrate that the HY5-binding motif predicted here was similar to the motif reported previously and that in vitro HY5-binding loci largely overlapped with the HY5-targeted candidate genes identified in previous ChIP-chip analysis. By combining these results with microarray analysis, we identified hundreds of HY5-binding genes that were differentially expressed in hy5. We also observed delayed induction of some transcripts of HY5-binding genes in hy5 mutants in response to blue-light exposure after dark treatment. Thus, an in vitro gDNA-binding assay coupled with sequencing is a convenient and powerful method to bridge the gap between identifying TF binding potential and establishing function.

Soybean miR172c targets the repressive AP2 transcription factor NNC1 to activate ENOD40 expression and regulate nodule initiation

MicroRNAs are noncoding RNAs that act as master regulators to modulate various biological processes by posttranscriptionally repressing their target genes. Repression of their target mRNA(s) can modulate signaling cascades and subsequent cellular events. Recently, a role for miR172 in soybean (Glycine max) nodulation has been described; however, the molecular mechanism through which miR172 acts to regulate nodulation has yet to be explored. Here, we demonstrate that soybean miR172c modulates both rhizobium infection and nodule organogenesis. miR172c was induced in soybean roots inoculated with either compatible Bradyrhizobium japonicum or lipooligosaccharide Nod factor and was highly upregulated during nodule development. Reduced activity and overexpression of miR172c caused dramatic changes in nodule initiation and nodule number. We show that soybean miR172c regulates nodule formation by repressing its target gene, Nodule Number Control1, which encodes a protein that directly targets the promoter of the early nodulin gene, ENOD40. Interestingly, transcriptional levels of miR172c were regulated by both Nod Factor Receptor1α/5α-mediated activation and by autoregulation of nodulation-mediated inhibition. Thus, we established a direct link between miR172c and the Nod factor signaling pathway in addition to adding a new layer to the precise nodulation regulation mechanism of soybean.

Using protein-binding microarrays to study transcription factor specificity: homologs, isoforms and complexes

Protein-DNA binding is central to specificity in gene regulation, and methods for characterizing transcription factor (TF)-DNA binding remain crucial to studies of regulatory specificity. High-throughput (HT) technologies have revolutionized our ability to characterize protein-DNA binding by significantly increasing the number of binding measurements that can be performed. Protein-binding microarrays (PBMs) are a robust and powerful HT platform for studying DNA-binding specificity of TFs. Analysis of PBM-determined DNA-binding profiles has provided new insight into the scope and mechanisms of TF binding diversity. In this review, we focus specifically on the PBM technique and discuss its application to the study of TF specificity, in particular, the binding diversity of TF homologs and multi-protein complexes.

Current trends and future directions in flower development research

Flowers, the reproductive structures of the approximately 400 000 extant species of flowering plants, exist in a tremendous range of forms and sizes, mainly due to developmental differences involving the number, arrangement, size and form of the floral organs of which they consist. However, this tremendous diversity is underpinned by a surprisingly robust basic floral structure in which a central group of carpels forms on an axis of determinate growth, almost invariably surrounded by two successive zones containing stamens and perianth organs, respectively. Over the last 25 years, remarkable progress has been achieved in describing the molecular mechanisms that control almost all aspects of flower development, from the phase change that initiates flowering to the final production of fruits and seeds. However, this work has been performed almost exclusively in a small number of eudicot model species, chief among which is Arabidopsis thaliana. Studies of flower development must now be extended to a much wider phylogenetic range of flowering plants and, indeed, to their closest living relatives, the gymnosperms. Studies of further, more wide-ranging models should provide insights that, for various reasons, cannot be obtained by studying the major existing models alone. The use of further models should also help to explain how the first flowering plants evolved from an unknown, although presumably gymnosperm-like ancestor, and rapidly diversified to become the largest major plant group and to dominate the terrestrial flora. The benefits for society of a thorough understanding of flower development are self-evident, as human life depends to a large extent on flowering plants and on the fruits and seeds they produce. In this preface to the Special Issue, we introduce eleven articles on flower development, representing work in both established and further models, including gymnosperms. We also present some of our own views on current trends and future directions of the flower development field.

Dehydration-induced endodormancy in crown buds of leafy spurge highlights involvement of MAF3- and RVE1-like homologs, and hormone signaling cross-talk

Vegetative shoot growth from underground adventitious buds of leafy spurge is critical for survival of this invasive perennial weed after episodes of severe abiotic stress. To determine the impact that dehydration-stress has on molecular mechanisms associated with vegetative reproduction of leafy spurge, greenhouse plants were exposed to mild- (3-day), intermediate- (7-day), severe- (14-day) and extended- (21-day) dehydration treatments. Aerial tissues of treated plants were then decapitated and soil was rehydrated to determine the growth potential of underground adventitious buds. Compared to well-watered plants, mild-dehydration accelerated new vegetative shoot growth, whereas intermediate- through extended-dehydration treatments both delayed and reduced shoot growth. Results of vegetative regrowth further confirmed that 14 days of dehydration induced a full-state of endodormancy in crown buds, which was correlated with a significant (P < 0.05) change in abundance of 2,124 transcripts. Sub-network enrichment analyses of transcriptome data obtained from the various levels of dehydration treatment also identified central hubs of over-represented genes involved in processes such as hormone signaling (i.e., ABA, auxin, ethylene, GA, and JA), response to abiotic stress (DREB1A/2A, RD22) and light (PIF3), phosphorylation (MPK4/6), circadian regulation (CRY2, PHYA), and flowering (AGL20, AP2, FLC). Further, results from this and previous studies highlight homologs most similar to Arabidopsis HY5, MAF3, RVE1 and RD22 as potential molecular markers for endodormancy in crown buds of leafy spurge. Early response to mild dehydration also highlighted involvement of upstream ethylene and JA-signaling, whereas severe dehydration impacted ABA-signaling. The identification of conserved ABRE- and MYC-consensus, cis-acting elements in the promoter of leafy spurge genomic clones similar to Arabidopsis RVE1 (AT5G17300) implicates a potential role for ABA-signaling in its dehydration-induced expression. Response of these molecular mechanisms to dehydration-stress provides insights on the ability of invasive perennial weeds to adapt and survive under harsh environments, which will be beneficial for addressing future management practices.

HsfB2b-mediated repression of PRR7 directs abiotic stress responses of the circadian clock

The circadian clock perceives environmental signals to reset to local time, but the underlying molecular mechanisms are not well understood. Here we present data revealing that a member of the heat shock factor (Hsf) family is involved in the input pathway to the plant circadian clock. Using the yeast one-hybrid approach, we isolated several Hsfs, including Heat Shock Factor B2b (HsfB2b), a transcriptional repressor that binds the promoter of Pseudo Response Regulator 7 (PRR7) at a conserved binding site. The constitutive expression of HsfB2b leads to severely reduced levels of the PRR7 transcript and late flowering and elongated hypocotyls. HsfB2b function is important during heat and salt stress because HsfB2b overexpression sustains circadian rhythms, and the hsfB2b mutant has a short circadian period under these conditions. HsfB2b is also involved in the regulation of hypocotyl growth under warm, short days. Our findings highlight the role of the circadian clock as an integrator of ambient abiotic stress signals important for the growth and fitness of plants.

Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants

BACKGROUND

Members of the ancient land-plant-specific transcription factor AT-Hook Motif Nuclear Localized (AHL) gene family regulate various biological processes. However, the relationships among the AHL genes, as well as their evolutionary history, still remain unexplored.

RESULTS

We analyzed over 500 AHL genes from 19 land plant species, ranging from the early diverging Physcomitrella patens and Selaginella to a variety of monocot and dicot flowering plants. We classified the AHL proteins into three types (Type-I/-II/-III) based on the number and composition of their functional domains, the AT-hook motif(s) and PPC domain. We further inferred their phylogenies via Bayesian inference analysis and predicted gene gain/loss events throughout their diversification. Our analyses suggested that the AHL gene family emerged in embryophytes and further evolved into two distinct clades, with Type-I AHLs forming one clade (Clade-A), and the other two types together diversifying in another (Clade-B). The two AHL clades likely diverged before the separation of Physcomitrella patens from the vascular plant lineage. In angiosperms, Clade-A AHLs expanded into 5 subfamilies; while, the ones in Clade-B expanded into 4 subfamilies. Examination of their expression patterns suggests that the AHLs within each clade share similar expression patterns with each other; however, AHLs in one monophyletic clade exhibit distinct expression patterns from the ones in the other clade. Over-expression of a Glycine max AHL PPC domain in Arabidopsis thaliana recapitulates the phenotype observed when over-expressing its Arabidopsis thaliana counterpart. This result suggests that the AHL genes from different land plant species may share conserved functions in regulating plant growth and development. Our study further suggests that such functional conservation may be due to conserved physical interactions among the PPC domains of AHL proteins.

CONCLUSIONS

Our analyses reveal a possible evolutionary scenario for the AHL gene family in land plants, which will facilitate the design of new studies probing their biological functions. Manipulating the AHL genes has been suggested to have tremendous effects in agriculture through increased seedling establishment, enhanced plant biomass and improved plant immunity. The information gleaned from this study, in turn, has the potential to be utilized to further improve crop production.

Large-scale identification of gibberellin-related transcription factors defines group VII ETHYLENE RESPONSE FACTORS as functional DELLA partners

DELLA proteins are the master negative regulators in gibberellin (GA) signaling acting in the nucleus as transcriptional regulators. The current view of DELLA action indicates that their activity relies on the physical interaction with transcription factors (TFs). Therefore, the identification of TFs through which DELLAs regulate GA responses is key to understanding these responses from a mechanistic point of view. Here, we have determined the TF interactome of the Arabidopsis (Arabidopsis thaliana) DELLA protein GIBBERELLIN INSENSITIVE and screened a collection of conditional TF overexpressors in search of those that alter GA sensitivity. As a result, we have found RELATED TO APETALA2.3, an ethylene-induced TF belonging to the group VII ETHYLENE RESPONSE FACTOR of the APETALA2/ethylene responsive element binding protein superfamily, as a DELLA interactor with physiological relevance in the context of apical hook development. The combination of transactivation assays and chromatin immunoprecipitation indicates that the interaction with GIBBERELLIN INSENSITIVE impairs the activity of RELATED TO APETALA2.3 on the target promoters. This mechanism represents a unique node in the cross regulation between the GA and ethylene signaling pathways controlling differential growth during apical hook development.

A functional and evolutionary perspective on transcription factor binding in Arabidopsis thaliana

Understanding the mechanisms underlying gene regulation is paramount to comprehend the translation from genotype to phenotype. The two are connected by gene expression, and it is generally thought that variation in transcription factor (TF) function is an important determinant of phenotypic evolution. We analyzed publicly available genome-wide chromatin immunoprecipitation experiments for 27 TFs in Arabidopsis thaliana and constructed an experimental network containing 46,619 regulatory interactions and 15,188 target genes. We identified hub targets and highly occupied target (HOT) regions, which are enriched for genes involved in development, stimulus responses, signaling, and gene regulatory processes in the currently profiled network. We provide several lines of evidence that TF binding at plant HOT regions is functional, in contrast to that in animals, and not merely the result of accessible chromatin. HOT regions harbor specific DNA motifs, are enriched for differentially expressed genes, and are often conserved across crucifers and dicots, even though they are not under higher levels of purifying selection than non-HOT regions. Distal bound regions are under purifying selection as well and are enriched for a chromatin state showing regulation by the Polycomb repressive complex. Gene expression complexity is positively correlated with the total number of bound TFs, revealing insights in the regulatory code for genes with different expression breadths. The integration of noncanonical and canonical DNA motif information yields new hypotheses on cobinding and tethering between specific TFs involved in flowering and light regulation.

Maize and millet transcription factors annotated using comparative genomic and transcriptomic data

Background

Transcription factors (TFs) contain DNA-binding domains (DBDs) and regulate gene expression by binding to specific DNA sequences. In addition, there are proteins, called transcription coregulators (TCs), which lack DBDs but can alter gene expression through interaction with TFs or RNA Polymerase II. Therefore, it is interesting to identify and classify the TFs and TCs in a genome. In this study, maize (Zea mays) and foxtail millet (Setaria italica), two important species for the study of C4 photosynthesis and kranz anatomy, were selected.

Result

We conducted a comprehensive genome-wide annotation of TFs and TCs in maize B73 and in two strains of foxtail millet, Zhang gu and Yugu1, and classified them into families. To gain additional support for our predictions, we searched for their homologous genes in Arabidopsis or rice and studied their gene expression level using RNA-seq and microarray data. We identified many new TF and TC families in these two species, and described some evolutionary and functional aspects of the 9 new maize TF families. Moreover, we detected many pseudogenes and transposable elements in current databases. In addition, we examined tissue expression preferences of TF and TC families and identified tissue/condition-specific TFs and TCs in maize and millet. Finally, we identified potential C4-related TF and TC genes in maize and millet.

Conclusions

Our results significantly expand current TF and TC annotations in maize and millet. We provided supporting evidence for our annotation from genomic and gene expression data and identified TF and TC genes with tissue preference in expression. Our study may facilitate the study of regulation of gene expression, tissue morphogenesis, and C4 photosynthesis in maize and millet. The data we generated in this study are available at http://sites.google.com/site/jjlmmtf.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-818) contains supplementary material, which is available to authorized users.

Mapping and dynamics of regulatory DNA and transcription factor networks in A. thaliana

Our understanding of gene regulation in plants is constrained by our limited knowledge of plant cis-regulatory DNA and its dynamics. We mapped DNase I hypersensitive sites (DHSs) in A. thaliana seedlings and used genomic footprinting to delineate ∼ 700,000 sites of in vivo transcription factor (TF) occupancy at nucleotide resolution. We show that variation associated with 72 diverse quantitative phenotypes localizes within DHSs. TF footprints encode an extensive cis-regulatory lexicon subject to recent evolutionary pressures, and widespread TF binding within exons may have shaped codon usage patterns. The architecture of A. thaliana TF regulatory networks is strikingly similar to that of animals in spite of diverged regulatory repertoires. We analyzed regulatory landscape dynamics during heat shock and photomorphogenesis, disclosing thousands of environmentally sensitive elements and enabling mapping of key TF regulatory circuits underlying these fundamental responses. Our results provide an extensive resource for the study of A. thaliana gene regulation and functional biology.

Class I TCP-DELLA interactions in inflorescence shoot apex determine plant height

Regulation of plant height, one of the most important agronomic traits, is the focus of intensive research for improving crop performance. Stem elongation takes place as a result of repeated cell divisions and subsequent elongation of cells produced by apical and intercalary meristems. The gibberellin (GA) phytohormones have long been known to control stem and internodal elongation by stimulating the degradation of nuclear growth-repressing DELLA proteins; however, the mechanism allowing GA-responsive growth is only slowly emerging. Here, we show that DELLAs directly regulate the activity of the plant-specific class I TCP transcription factor family, key regulators of cell proliferation. Our results demonstrate that class I TCP factors directly bind the promoters of core cell-cycle genes in Arabidopsis inflorescence shoot apices while DELLAs block TCP function by binding to their DNA-recognition domain. GAs antagonize such repression by promoting DELLA destruction and therefore cause a concomitant accumulation of TCP factors on promoters of cell-cycle genes. Consistent with this model, the quadruple mutant tcp8 tcp14 tcp15 tcp22 exhibits severe dwarfism and reduced responsiveness to GA action. Altogether, we conclude that GA-regulated DELLA-TCP interactions in inflorescence shoot apex provide a novel mechanism to control plant height.

Inference of transcriptional networks in Arabidopsis through conserved noncoding sequence analysis

Transcriptional regulation plays an important role in establishing gene expression profiles during development or in response to (a)biotic stimuli. Transcription factor binding sites (TFBSs) are the functional elements that determine transcriptional activity, and the identification of individual TFBS in genome sequences is a major goal to inferring regulatory networks. We have developed a phylogenetic footprinting approach for the identification of conserved noncoding sequences (CNSs) across 12 dicot plants. Whereas both alignment and non-alignment-based techniques were applied to identify functional motifs in a multispecies context, our method accounts for incomplete motif conservation as well as high sequence divergence between related species. We identified 69,361 footprints associated with 17,895 genes. Through the integration of known TFBS obtained from the literature and experimental studies, we used the CNSs to compile a gene regulatory network in Arabidopsis thaliana containing 40,758 interactions, of which two-thirds act through binding events located in DNase I hypersensitive sites. This network shows significant enrichment toward in vivo targets of known regulators, and its overall quality was confirmed using five different biological validation metrics. Finally, through the integration of detailed expression and function information, we demonstrate how static CNSs can be converted into condition-dependent regulatory networks, offering opportunities for regulatory gene annotation.

Expression-based network biology identifies immune-related functional modules involved in plant defense

Background

Plants respond to diverse environmental cues including microbial perturbations by coordinated regulation of thousands of genes. These intricate transcriptional regulatory interactions depend on the recognition of specific promoter sequences by regulatory transcription factors. The combinatorial and cooperative action of multiple transcription factors defines a regulatory network that enables plant cells to respond to distinct biological signals. The identification of immune-related modules in large-scale transcriptional regulatory networks can reveal the mechanisms by which exposure to a pathogen elicits a precise phenotypic immune response.

Results

We have generated a large-scale immune co-expression network using a comprehensive set of Arabidopsis thaliana (hereafter Arabidopsis) transcriptomic data, which consists of a wide spectrum of immune responses to pathogens or pathogen-mimicking stimuli treatments. We employed both linear and non-linear models to generate Arabidopsis immune co-expression regulatory (AICR) network. We computed network topological properties and ascertained that this newly constructed immune network is densely connected, possesses hubs, exhibits high modularity, and displays hallmarks of a “real” biological network. We partitioned the network and identified 156 novel modules related to immune functions. Gene Ontology (GO) enrichment analyses provided insight into the key biological processes involved in determining finely tuned immune responses. We also developed novel software called OCCEAN (One Click Cis-regulatory Elements ANalysis) to discover statistically enriched promoter elements in the upstream regulatory regions of Arabidopsis at a whole genome level. We demonstrated that OCCEAN exhibits higher precision than the existing promoter element discovery tools. In light of known and newly discovered cis-regulatory elements, we evaluated biological significance of two key immune-related functional modules and proposed mechanism(s) to explain how large sets of diverse GO genes coherently function to mount effective immune responses.

Conclusions

We used a network-based, top-down approach to discover immune-related modules from transcriptomic data in Arabidopsis. Detailed analyses of these functional modules reveal new insight into the topological properties of immune co-expression networks and a comprehensive understanding of multifaceted plant defense responses. We present evidence that our newly developed software, OCCEAN, could become a popular tool for the Arabidopsis research community as well as potentially expand to analyze other eukaryotic genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-421) contains supplementary material, which is available to authorized users.

Down-regulation of AUXIN RESPONSE FACTORS 6 and 8 by microRNA 167 leads to floral development defects and female sterility in tomato

Auxin regulates the expression of diverse genes that affect plant growth and development. This regulation requires AUXIN RESPONSE FACTORS (ARFs) that bind to the promoter regions of these genes. ARF6 and ARF8 in Arabidopsis thaliana are required to promote inflorescence stem elongation and late stages of petal, stamen, and gynoecium development. All seed plants studied thus far have ARF6 and ARF8 orthologues as well as the microRNA miR167, which targets ARF6 and ARF8. Whether these genes have broadly conserved roles in flower development is not known. To address this question, the effects of down-regulation of ARF6 and ARF8 were investigated through transgenic expression of Arabidopsis MIR167a in tomato, which diverged from Arabidopsis before the radiation of dicotyledonous plants approximately 90-112 million years ago. The transgenic tomato plants overexpressing MIR167a exhibited reductions in leaf size and internode length as well as shortened petals, stamens, and styles. More significantly, the transgenic plants were female-sterile as a result of failure of wild-type pollen to germinate on the stigma surface and/or to grow through the style. RNA-Seq analysis identified many genes with significantly altered expression patterns, including those encoding products with functions in 'transcription regulation', 'cell wall' and 'lipid metabolism' categories. Putative orthologues of a subset of these genes were also differentially expressed in Arabidopsis arf6 arf8 mutant flowers. These results thus suggest that ARF6 and ARF8 have conserved roles in controlling growth and development of vegetative and flower organs in dicots.

Motif enrichment tool

The Motif Enrichment Tool (MET) provides an online interface that enables users to find major transcriptional regulators of their gene sets of interest. MET searches the appropriate regulatory region around each gene and identifies which transcription factor DNA-binding specificities (motifs) are statistically overrepresented. Motif enrichment analysis is currently available for many metazoan species including human, mouse, fruit fly, planaria and flowering plants. MET also leverages high-throughput experimental data such as ChIP-seq and DNase-seq from ENCODE and ModENCODE to identify the regulatory targets of a transcription factor with greater precision. The results from MET are produced in real time and are linked to a genome browser for easy follow-up analysis. Use of the web tool is free and open to all, and there is no login requirement. Address: http://veda.cs.uiuc.edu/MET/.