Arabidopsis AZG2 transports cytokinins in vivo and regulates lateral root emergence.

Cytokinin and auxin are key regulators of plant growth and development. During the last decade transport mechanisms have turned out to be the key for the control of local and long-distance hormone distributions. In contrast with auxin, cytokinin transport is poorly understood. Here, we show that Arabidopsis thaliana AZG2, a member of the AZG purine transporter family, acts as cytokinin transporter involved in root system architecture determination. Even though purines are substrates for both AZG1 and AZG2, we found distinct transport mechanisms. The expression of AZG2 is restricted to a small group of cells surrounding the lateral root (LR) primordia and induced by auxins. Compared to the wild-type (WT), mutants carrying loss-of-function alleles of AZG2 have higher LR density, suggesting that AZG2 is part of a regulatory pathway in LR emergence. Moreover, azg2 is partially insensitive to exogenous cytokinin, which is consistent with the observation that the cytokinin reporter TCSnpro :GFP showed lower fluorescence signal in the roots of azg2 compared to the WT. These results indicate a defective cytokinin signalling pathway in the region of LR primordia. The integration of AZG2 subcellular localization and cytokinin transport capacity data allowed us to propose a local cytokinin : auxin signalling model for the regulation of LR emergence.

Role of BASIC PENTACYSTEINE transcription factors in a subset of cytokinin signaling responses.

Cytokinin plays diverse roles in plant growth and development, generally acting by modulating gene transcription in target tissues. The type-B Arabidopsis response regulators (ARR) transcription factors have emerged as primary targets of cytokinin signaling and are required for essentially all cytokinin-mediated changes in gene expression. The diversity of cytokinin function is likely imparted by the activity of various transcription factors working with the type-B ARRs to alter specific sets of target genes. One potential set of co-regulators modulating the cytokinin response are the BARLEY B-RECOMBINANT/BASIC PENTACYSTEINE (BBR/BPC) family of plant-specific transcription factors. Here, we show that disruption of multiple BPCs results in reduced sensitivity to cytokinin. Further, the BPCs are necessary for the induction of a subset of genes in response to cytokinin. We identified direct in vivo targets of BPC6 using ChIP-Seq and found an enrichment of promoters of genes differentially expressed in response to cytokinin. Further, a significant number of BPC6 regulated genes are also direct targets of the type-B ARRs. Potential cis-binding elements for a number of other transcription factors linked to cytokinin action are enriched in the BPC binding fragments, including those for the cytokinin response factors (CRFs). In addition, several BPCs interact with a subset of type-A ARRs. Consistent with these results, a significant number of genes whose expression is altered in bpc mutant roots are also mis-expressed in crf1,3,5,6 and type-A arr3,4,5,6,7,8,9,15 mutant roots. These results suggest that the BPCs are part of a complex network of transcription factors that are involved in the response to cytokinin.

The Enzyme-Like Domain of Arabidopsis Nuclear ß-Amylases Is Critical for DNA Sequence Recognition and Transcriptional Activation.

Plant BZR1-BAM transcription factors contain a ß-amylase (BAM)-like domain, characteristic of proteins involved in starch breakdown. The enzyme-derived domains appear to be noncatalytic, but they determine the function of the two Arabidopsis thaliana BZR1-BAM isoforms (BAM7 and BAM8) during transcriptional initiation. Removal or swapping of the BAM domains demonstrates that the BAM7 BAM domain restricts DNA binding and transcriptional activation, while the BAM8 BAM domain allows both activities. Furthermore, we demonstrate that BAM7 and BAM8 interact on the protein level and cooperate during transcriptional regulation. Site-directed mutagenesis of residues in the BAM domain of BAM8 shows that its function as a transcriptional activator is independent of catalysis but requires an intact substrate binding site, suggesting it may bind a ligand. Microarray experiments with plants overexpressing truncated versions lacking the BAM domain indicate that the pseudo-enzymatic domain increases selectivity for the preferred cis-regulatory element BBRE (BZR1-BAM Responsive Element). Side specificity toward the G-box may allow crosstalk to other signaling networks. This work highlights the importance of the enzyme-derived domain of BZR1-BAMs, supporting their potential role as metabolic sensors.

The Arabidopsis GAGA-Binding Factor BASIC PENTACYSTEINE6 Recruits the POLYCOMB-REPRESSIVE COMPLEX1 Component LIKE HETEROCHROMATIN PROTEIN1 to GAGA DNA Motifs.

Polycomb-repressive complexes (PRCs) play key roles in development by repressing a large number of genes involved in various functions. Much, however, remains to be discovered about PRC-silencing mechanisms as well as their targeting to specific genomic regions. Besides other mechanisms, GAGA-binding factors in animals can guide PRC members in a sequence-specific manner to Polycomb-responsive DNA elements. Here, we show that the Arabidopsis (Arabidopsis thaliana) GAGA-motif binding factor protein basic pentacysteine6 (BPC6) interacts with like heterochromatin protein1 (LHP1), a PRC1 component, and associates with vernalization2 (VRN2), a PRC2 component, in vivo. By using a modified DNA-protein interaction enzyme-linked immunosorbant assay, we could show that BPC6 was required and sufficient to recruit LHP1 to GAGA motif-containing DNA probes in vitro. We also found that LHP1 interacts with VRN2 and, therefore, can function as a possible scaffold between BPC6 and VRN2. The lhp1-4 bpc4 bpc6 triple mutant displayed a pleiotropic phenotype, extreme dwarfism and early flowering, which disclosed synergistic functions of LHP1 and group II plant BPC members. Transcriptome analyses supported this synergy and suggested a possible function in the concerted repression of homeotic genes, probably through histone H3 lysine-27 trimethylation. Hence, our findings suggest striking similarities between animal and plant GAGA-binding factors in the recruitment of PRC1 and PRC2 components to Polycomb-responsive DNA element-like GAGA motifs, which must have evolved through convergent evolution.

Elucidating the evolutionary conserved DNA-binding specificities of WRKY transcription factors by molecular dynamics and in vitro binding assays.

WRKY transcription factors constitute a large protein family in plants that is involved in the regulation of developmental processes and responses to biotic or abiotic stimuli. The question arises how stimulus-specific responses are mediated given that the highly conserved WRKY DNA-binding domain (DBD) exclusively recognizes the 'TTGACY' W-box consensus. We speculated that the W-box consensus might be more degenerate and yet undetected differences in the W-box consensus of WRKYs of different evolutionary descent exist. The phylogenetic analysis of WRKY DBDs suggests that they evolved from an ancestral group IIc-like WRKY early in the eukaryote lineage. A direct descent of group IIc WRKYs supports a monophyletic origin of all other group II and III WRKYs from group I by loss of an N-terminal DBD. Group I WRKYs are of paraphyletic descent and evolved multiple times independently. By homology modeling, molecular dynamics simulations and in vitro DNA-protein interaction-enzyme-linked immunosorbent assay with AtWRKY50 (IIc), AtWRKY33 (I) and AtWRKY11 (IId) DBDs, we revealed differences in DNA-binding specificities. Our data imply that other components are essentially required besides the W-box-specific binding to DNA to facilitate a stimulus-specific WRKY function.

Volatiles of two growth-inhibiting rhizobacteria commonly engage AtWRKY18 function.

Interactions with the (a)biotic environment play key roles in a plant's fitness and vitality. In addition to direct surface-to-surface contact, volatile chemicals can also affect the physiology of organism. Volatiles of Serratia plymuthica and Stenotrophomonas maltophilia significantly inhibited growth and induced H(2) O(2) production in Arabidopsis in dual culture. Within 1 day, transcriptional changes were observed by promoter-GUS assays using a stress-inducible W-box-containing 4xGST1 construct. Expression studies performed at 6, 12 and 24 h revealed altered transcript levels for 889 genes and 655 genes in response to Se. plymuthica or St. maltophilia volatiles, respectively. Expression of 162 genes was altered in both treatments. Meta-analysis revealed that specifically volatile-responsive genes were significantly overlapping with those affected by abiotic stress. We use the term mVAMP (microbial volatile-associated molecular pattern) to describe these volatile-specific responses. Genes responsive to both treatments were enriched for W-box motifs in their promoters, and were significantly enriched for transcription factors (ERF2, ZAT10, MYB73 and WRKY18). The susceptibility of wrky18 mutant lines to volatiles was significantly delayed, suggesting an indispensable role for WRKY18 in bacterial volatile responses.

Prerequisites, performance and profits of transcriptional profiling the abiotic stress response.

During the last decade, microarrays became a routine tool for the analysis of transcripts in the model plant Arabidopsis thaliana and the crop plant species rice, poplar or barley. The overwhelming amount of data generated by gene expression studies is a valuable resource for every scientist. Here, we summarize the most important findings about the abiotic stress responses in plants. Interestingly, conserved patterns of gene expression responses have been found that are common between different abiotic stresses or that are conserved between different plant species. However, the individual histories of each plant affect the inter-comparability between experiments already before the onset of the actual stress treatment. This review outlines multiple aspects of microarray technology and highlights some of the benefits, limitations and also pitfalls of the technique. This article is part of a Special Issue entitled: Plant gene regulation in response to abiotic stress.

The SCO2 protein disulphide isomerase is required for thylakoid biogenesis and interacts with LHCB1 chlorophyll a/b binding proteins which affects chlorophyll biosynthesis in Arabidopsis seedlings.

The process of chloroplast biogenesis requires a multitude of pathways and processes to establish chloroplast function. In cotyledons of seedlings, chloroplasts develop either directly from proplastids (also named eoplasts) or, if germinated in the dark, via etioplasts, whereas in leaves chloroplasts derive from proplastids in the apical meristem and are then multiplied by division. The snowy cotyledon 2, sco2, mutations specifically disrupt chloroplast biogenesis in cotyledons. SCO2 encodes a chloroplast-localized protein disulphide isomerase, hypothesized to be involved in protein folding. Analysis of co-expressed genes with SCO2 revealed that genes with similar expression patterns encode chloroplast proteins involved in protein translation and in chlorophyll biosynthesis. Indeed, sco2-1 accumulates increased levels of the chlorophyll precursor, protochlorophyllide, in both dark grown cotyledons and leaves. Yeast two-hybrid analyses demonstrated that SCO2 directly interacts with the chlorophyll-binding LHCB1 proteins, being confirmed in planta using bimolecular fluorescence complementation (BIFC). Furthermore, ultrastructural analysis of sco2-1 chloroplasts revealed that formation and movement of transport vesicles from the inner envelope to the thylakoids is perturbed. SCO2 does not interact with the signal recognition particle proteins SRP54 and FtsY, which were shown to be involved in targeting of LHCB1 to the thylakoids. We hypothesize that SCO2 provides an alternative targeting pathway for light-harvesting chlorophyll binding (LHCB) proteins to the thylakoids via transport vesicles predominantly in cotyledons, with the signal recognition particle (SRP) pathway predominant in rosette leaves. Therefore, we propose that SCO2 is involved in the integration of LHCB1 proteins into the thylakoids that feeds back on the regulation of the tetrapyrrole biosynthetic pathway and nuclear gene expression.

Plant core environmental stress response genes are systemically coordinated during abiotic stresses.

Studying plant stress responses is an important issue in a world threatened by global warming. Unfortunately, comparative analyses are hampered by varying experimental setups. In contrast, the AtGenExpress abiotic stress experiment displays intercomparability. Importantly, six of the nine stresses (wounding, genotoxic, oxidative, UV-B light, osmotic and salt) can be examined for their capacity to generate systemic signals between the shoot and root, which might be essential to regain homeostasis in Arabidopsis thaliana. We classified the systemic responses into two groups: genes that are regulated in the non-treated tissue only are defined as type I responsive and, accordingly, genes that react in both tissues are termed type II responsive. Analysis of type I and II systemic responses suggest distinct functionalities, but also significant overlap between different stresses. Comparison with salicylic acid (SA) and methyl-jasmonate (MeJA) responsive genes implies that MeJA is involved in the systemic stress response. Certain genes are predominantly responding in only one of the categories, e.g., WRKY genes respond mainly non-systemically. Instead, genes of the plant core environmental stress response (PCESR), e.g., ZAT10, ZAT12, ERD9 or MES9, are part of different response types. Moreover, several PCESR genes switch between the categories in a stress-specific manner.

Bioinformatic cis-element analyses performed in Arabidopsis and rice disclose bZIP- and MYB-related binding sites as potential AuxRE-coupling elements in auxin-mediated transcription.

BACKGROUND: In higher plants, a diverse array of developmental and growth-related processes is regulated by the plant hormone auxin. Recent publications have proposed that besides the well-characterized Auxin Response Factors (ARFs) that bind Auxin Response Elements (AuxREs), also members of the bZIP- and MYB-transcription factor (TF) families participate in transcriptional control of auxin-regulated genes via bZIP Response Elements (ZREs) or Myb Response Elements (MREs), respectively. RESULTS: Applying a novel bioinformatic algorithm, we demonstrate on a genome-wide scale that singular motifs or composite modules of AuxREs, ZREs, MREs but also of MYC2 related elements are significantly enriched in promoters of auxin-inducible genes. Despite considerable, species-specific differences in the genome structure in terms of the GC content, this enrichment is generally conserved in dicot (Arabidopsis thaliana) and monocot (Oryza sativa) model plants. Moreover, an enrichment of defined composite modules has been observed in selected auxin-related gene families. Consistently, a bipartite module, which encompasses a bZIP-associated G-box Related Element (GRE) and an AuxRE motif, has been found to be highly enriched. Making use of transient reporter studies in protoplasts, these findings were experimentally confirmed, demonstrating that GREs functionally interact with AuxREs in regulating auxin-mediated transcription. CONCLUSIONS: Using genome-wide bioinformatic analyses, evolutionary conserved motifs have been defined which potentially function as AuxRE-dependent coupling elements to establish auxin-specific expression patterns. Based on these findings, experimental approaches can be designed to broaden our understanding of combinatorial, auxin-controlled gene regulation.

Significance of light, sugar, and amino acid supply for diurnal gene regulation in developing barley caryopses.

The caryopses of barley (Hordeum vulgare), as of all cereals, are complex sink organs optimized for starch accumulation and embryo development. While their early to late development has been studied in great detail, processes underlying the caryopses' diurnal adaptation to changes in light, temperature, and the fluctuations in phloem-supplied carbon and nitrogen have remained unknown. In an attempt to identify diurnally affected processes in developing caryopses at the early maturation phase, we monitored global changes of both gene expression and metabolite levels. We applied the 22 K Barley1 GeneChip microarray and identified 2,091 differentially expressed (DE) genes that were assigned to six major diurnal expression clusters. Principal component analysis and other global analyses demonstrated that the variability within the data set relates to genes involved in circadian regulation, storage compound accumulation, embryo development, response to abiotic stress, and photosynthesis. The correlation of amino acid and sugar profiles with expression trajectories led to the identification of several hundred potentially metabolite-regulated DE genes. A comparative analysis of our data set and publicly available microarray data disclosed suborgan-specific expression of almost all diurnal DE genes, with more than 350 genes specifically expressed in the pericarp, endosperm, or embryo tissues. Our data reveal a tight linkage between day/night cycles, changes in light, and the supply of carbon and nitrogen. We present a model that suggests several phases of diurnal gene expression in developing barley caryopses, summarized as starvation and priming, energy collection and carbon fixation, light protection and chaperone activity, storage and growth, and embryo development.

The ABA-mediated switch between submersed and emersed life-styles in aquatic macrophytes.

Hydrophytes comprise aquatic macrophytes from various taxa that are able to sustain and to complete their lifecycle in a flooded environment. Their ancestors, however, underwent adaptive processes to withstand drought on land and became partially or completely independent of water for sexual reproduction. Interestingly, the step backwards into the high-density aquatic medium happened independently several times in numerous plant taxa. For flowering plants, this submersed life-style is especially difficult as they need to erect their floral organs above the water surface to be pollinated. Moreover, fresh-water plants evolved the adaptive mechanism of heterophylly, which enabled them to switch between a submersed and an emersed leaf morphology. The plant hormone abscisic acid (ABA) is a key factor of heterophylly induction in aquatic plants and is a major switch between a submersed and emersed life. The mechanisms of ABA signal perception and transduction appear to be conserved throughout the evolution of basal plants to angiosperms and from terrestrial to aquatic plants. This review summarizes the interplay of environmental factors that act through ABA to orchestrate adaptation of plants to their aquatic environment.

Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system.

Sex determination is a central process for sexual reproduction and is often regulated by a sex determinant encoded on a sex chromosome. Rules that govern the evolution of sex chromosomes via specialization and degeneration following the evolution of a sex determinant have been well studied in diploid organisms. However, distinct predictions apply to sex chromosomes in organisms where sex is determined in the haploid phase of the life cycle: both sex chromosomes, female U and male V, are expected to maintain their gene functions, even though both are non-recombining. This is in contrast to the X-Y (or Z-W) asymmetry and Y (W) chromosome degeneration in XY (ZW) systems of diploids. Here, we provide evidence that sex chromosomes diverged early during the evolution of haploid liverworts and identify the sex determinant on the Marchantia polymorpha U chromosome. This gene, Feminizer, encodes a member of the plant-specific BASIC PENTACYSTEINE transcription factor family. It triggers female differentiation via regulation of the autosomal sex-determining locus of FEMALE GAMETOPHYTE MYB and SUPPRESSOR OF FEMINIZATION. Phylogenetic analyses of Feminizer and other sex chromosome genes indicate dimorphic sex chromosomes had already been established 430 mya in the ancestral liverwort. Feminizer also plays a role in reproductive induction that is shared with its gametolog on the V chromosome, suggesting an ancestral function, distinct from sex determination, was retained by the gametologs. This implies ancestral functions can be preserved after the acquisition of a sex determination mechanism during the evolution of a dominant haploid sex chromosome system.

Application of FLIM-FIDSAM for the in vivo analysis of hormone competence of different cell types.

Background fluorescence derived from subcellular compartments is a major drawback in high-resolution live imaging, especially of plant cells. A novel technique for contrast enhancement of fluorescence images of living cells expressing fluorescent fusion proteins termed fluorescence intensity decay shape analysis microscopy (FIDSAM) has been recently published and is applied here to plant cells expressing wild-type levels of a low-abundant membrane protein (BRI1-EGFP), demonstrating the applicability of FIDSAM to samples exhibiting about 80% autofluorescence. Furthermore, the combination of FIDSAM and fluorescence lifetime imaging microscopy enables the simultaneous determination and quantification of different ligand-specific responses in living cells with high spatial and temporal resolution even in samples with high autofluorescence background. Correlation of different responses can be used to determine the hormone ligand competence of different cell types as demonstrated here in BRI1-EGFP-expressing root and hypocotyl cells.

Analysis of plant regulatory DNA sequences by the yeast-one-hybrid assay.

Regulatory DNA sequences harbor the essential information to control specific gene expression changes and integrate information derived from upstream signaling cascades. This regulatory potential is mediated by direct binding of proteins, e.g., transcription factors, to defined stretches of DNA motifs in regulatory regions. The analysis of these DNA regions, at which several signaling pathways could merge to orchestrate gene expression, is still a challenging task. To date, the combination of functional approaches in the laboratory and computer aided sequence evaluation is frequently used for regulatory sequence analysis. The yeast-one-hybrid method is a possible approach to test for direct binding of plant proteins to DNA in a heterologous system. Moreover, it is the most frequently used method for the identification of DNA-binding proteins targeting a given DNA sequence by screening a cDNA library.

Analysis of plant regulatory DNA sequences by transient protoplast assays and computer aided sequence evaluation.

The orchestrated regulation of hundreds of genes responding in a temporal, spatial, and conditional expression is in part mediated by the transient binding of transcription factors to their specific DNA motifs. The analysis of these cis-regulatory DNA sequences is still a challenging task. Therefore, the combination of the transient protoplast expression assay with computer aided sequence analysis is a preferred method for regulatory sequence analysis. The protocols given here describe the use of a reporter gene plasmid to investigate the effects of a transiently co-expressed transcription factor gene in planta. As a number of bioinformatic analysis tools for cis-elements are publicly available, we suggest a workflow for the electronic analysis of promoter sequences. This is in particular difficult, as most analysis programs have not been developed for the investigation of single sequences. We provide insight into bioinformatics tools that investigate cis-element presence, their distribution, and implications for further functional analyses.

Phylogenetic Analyses and GAGA-Motif Binding Studies of BBR/BPC Proteins Lend to Clues in GAGA-Motif Recognition and a Regulatory Role in Brassinosteroid Signaling.

Plant GAGA-motif binding factors are encoded by the BARLEY B RECOMBINANT / BASIC PENTACYSTEINE (BBR/BPC) family, which fulfill indispensable functions in growth and development. BBR/BPC proteins control flower development, size of the stem cell niche and seed development through transcriptional regulation of homeotic transcription factor genes. They are responsible for the context dependent recruitment of Polycomb repressive complexes (PRC) or other repressive proteins to GAGA-motifs, which are contained in Polycomb repressive DNA-elements (PREs). Hallmark of the protein family is the highly conserved BPC domain, which is required for DNA binding. Here we study the evolution and diversification of the BBR/BPC family and its DNA-binding domain. Our analyses supports a further division of the family into four main groups (I-IV) and several subgroups, to resolve a strict monophyletic descent of the BPC domain. We prove a polyphyletic origin for group III proteins, which evolved from group I and II members through extensive loss of domains in the N-terminus. Conserved motif searches lend to the identification of a WAR/KHGTN consensus and a TIR/K motif at the very C-terminus of the BPC-domain. We could show by DPI-ELISA that this signature is required for DNA-binding in AtBPC1. Additional binding studies with AtBPC1, AtBPC6 and mutated oligonucleotides consolidated the binding to GAGA tetramers. To validate these findings, we used previously published ChIP-seq data from GFP-BPC6. We uncovered that many genes of the brassinosteroid signaling pathway are targeted by AtBPC6. Consistently, bpc6, bpc4 bpc6, and lhp1 bpc4 bpc4 mutants display brassinosteroid-dependent root growth phenotypes. Both, a function in brassinosteroid signaling and our phylogenetic data supports a link between BBR/BPC diversification in the land plant lineage and the complexity of flower and seed plant evolution.

Phylogenetic and comparative gene expression analysis of barley (Hordeum vulgare) WRKY transcription factor family reveals putatively retained functions between monocots and dicots.

BACKGROUND: WRKY proteins belong to the WRKY-GCM1 superfamily of zinc finger transcription factors that have been subject to a large plant-specific diversification. For the cereal crop barley (Hordeum vulgare), three different WRKY proteins have been characterized so far as regulators in sucrose signaling, pathogen defense, and in response to cold and drought. However, their phylogenetic relationship remained unresolved. RESULTS: In this study, we used available sequence information to identify a minimum number of 45 barley WRKY transcription factor (HvWRKY) genes. According to their structural features, the HvWRKY factors were classified into the previously defined polyphyletic WRKY subgroups 1 to 3. Furthermore, we could assign putative orthologs of the HvWRKY proteins in Arabidopsis and rice. While in most cases clades of orthologous proteins were formed within each group or subgroup, other clades were composed of paralogous proteins for the grasses and Arabidopsis only, which is indicative of specific gene radiation events. To gain insight into their putative functions, we examined expression profiles of WRKY genes from publicly available microarray data resources and found group specific expression patterns. While putative orthologs of the HvWRKY transcription factors have been inferred from phylogenetic sequence analysis, we performed a comparative expression analysis of WRKY genes in Arabidopsis and barley. Indeed, highly correlative expression profiles were found between some of the putative orthologs. CONCLUSION: HvWRKY genes have not only undergone radiation in monocot or dicot species, but exhibit evolutionary traits specific to grasses. HvWRKY proteins exhibited not only sequence similarities between orthologs with Arabidopsis, but also relatedness in their expression patterns. This correlative expression is indicative for a putative conserved function of related WRKY proteins in monocot and dicot species.

EDISA: extracting biclusters from multiple time-series of gene expression profiles.

BACKGROUND: Cells dynamically adapt their gene expression patterns in response to various stimuli. This response is orchestrated into a number of gene expression modules consisting of co-regulated genes. A growing pool of publicly available microarray datasets allows the identification of modules by monitoring expression changes over time. These time-series datasets can be searched for gene expression modules by one of the many clustering methods published to date. For an integrative analysis, several time-series datasets can be joined into a three-dimensional gene-condition-time dataset, to which standard clustering or biclustering methods are, however, not applicable. We thus devise a probabilistic clustering algorithm for gene-condition-time datasets. RESULTS: In this work, we present the EDISA (Extended Dimension Iterative Signature Algorithm), a novel probabilistic clustering approach for 3D gene-condition-time datasets. Based on mathematical definitions of gene expression modules, the EDISA samples initial modules from the dataset which are then refined by removing genes and conditions until they comply with the module definition. A subsequent extension step ensures gene and condition maximality. We applied the algorithm to a synthetic dataset and were able to successfully recover the implanted modules over a range of background noise intensities. Analysis of microarray datasets has lead us to define three biologically relevant module types: 1) We found modules with independent response profiles to be the most prevalent ones. These modules comprise genes which are co-regulated under several conditions, yet with a different response pattern under each condition. 2) Coherent modules with similar responses under all conditions occurred frequently, too, and were often contained within these modules. 3) A third module type, which covers a response specific to a single condition was also detected, but rarely. All of these modules are essentially different types of biclusters. CONCLUSION: We successfully applied the EDISA to different 3D datasets. While previous studies were mostly aimed at detecting coherent modules only, our results show that coherent responses are often part of a more general module type with independent response profiles under different conditions. Our approach thus allows for a more comprehensive view of the gene expression response. After subsequent analysis of the resulting modules, the EDISA helped to shed light on the global organization of transcriptional control. An implementation of the algorithm is available at http://www-ra.informatik.uni-tuebingen.de/software/IAGEN/.

Cis-motifs upstream of the transcription and translation initiation sites are effectively revealed by their positional disequilibrium in eukaryote genomes using frequency distribution curves.

BACKGROUND: The discovery of cis-regulatory motifs still remains a challenging task even though the number of sequenced genomes is constantly growing. Computational analyses using pattern search algorithms have been valuable in phylogenetic footprinting approaches as have expression profile experiments to predict co-occurring motifs. Surprisingly little is known about the nature of cis-regulatory element (CRE) distribution in promoters. RESULTS: In this paper we used the Motif Mapper open-source collection of visual basic scripts for the analysis of motifs in any aligned set of DNA sequences. We focused on promoter motif distribution curves to identify positional over-representation of DNA motifs. Using differentially aligned datasets from the model species Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster and Saccharomyces cerevisiae, we convincingly demonstrated the importance of the position and orientation for motif discovery. Analysis with known CREs and all possible hexanucleotides showed that some functional elements gather close to the transcription and translation initiation sites and that elements other than the TATA-box motif are conserved between eukaryote promoters. While a high background frequency usually decreases the effectiveness of such an enumerative investigation, we improved our analysis by conducting motif distribution maps using large datasets. CONCLUSION: This is the first study to reveal positional over-representation of CREs and promoter motifs in a cross-species approach. CREs and motifs shared between eukaryotic promoters support the observation that an eukaryotic promoter structure has been conserved throughout evolutionary time. Furthermore, with the information on positional enrichment of a motif or a known functional CRE, it is possible to get a more detailed insight into where an element appears to function. This in turn might accelerate the in depth examination of known and yet unknown cis-regulatory sequences in the laboratory.