ClusterDraw web server: a tool to identify and visualize clusters of binding motifs for transcription factors

The dorsoventral patterning of Musca domestica embryos: insights into BMP/Dpp evolution from the base of the lower cyclorraphan flies

Background

In the last few years, accumulated information has indicated that the evolution of an extra-embryonic membrane in dipterans was accompanied by changes in the gene regulatory network controlled by the BMP/Dpp pathway, which is responsible for dorsal patterning in these insects. However, only comparative analysis of gene expression levels between distant species with two extra-embryonic membranes, like A. gambiae or C. albipunctata, and D. melanogaster, has been conducted. Analysis of gene expression in ancestral species, which evolved closer to the amnioserosa origin, could provide new insights into the evolution of dorsoventral patterning in dipterans.

Results

Here we describe the spatial expression of several key and downstream elements of the Dpp pathway and show the compared patterns of expression between Musca and Drosophila embryos, both dipterans with amnioserosa. Most of the analyzed gene showed a high degree of expression conservation, however, we found several differences in the gene expression pattern of M. domestica orthologs for sog and tolloid. Bioinformatics analysis of the promoter of both genes indicated that the variations could be related to the gain of several binding sites for the transcriptional factor Dorsal in the Md.tld promoter and Snail in the Md.sog enhancer. These altered expressions could explain the unclear formation of the pMad gradient in the M. domestica embryo, compared to the formation of the gradient in D. melanogaster.

Conclusion

Gene expression changes during the dorsal–ventral patterning in insects contribute to the differentiation of extra-embryonic tissues as a consequence of changes in the gene regulatory network controlled by BMP/Dpp. In this work, in early M. domestica embryos, we identified the expression pattern of several genes members involved in the dorsoventral specification of the embryo. We believe that these data can contribute to understanding the evolution of the BMP/Dpp pathway, the regulation of BMP ligands, and the formation of a Dpp gradient in higher cyclorraphan flies.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0102-5) contains supplementary material, which is available to authorized users.

Transcriptional activity of the short gastrulation primary enhancer in the ventral midline requires its early activity in the presumptive neurogenic ectoderm

The short gastrulation (sog) shadow enhancer directs early and late sog expression in the neurogenic ectoderm and the ventral midline of the developing Drosophila embryo, respectively. Here, evidence is presented that the sog primary enhancer also has both activities, with the late enhancer activity dependent on the early activity. Computational analyses showed that the sog primary enhancer contains five Dorsal (Dl)-, four Zelda (Zld)-, three Bicoid (Bcd)-, and no Single-minded (Sim)-binding sites. In contrast to many ventral midline enhancers, the primary enhancer can direct lacZ expression in the ventral midline as well as in the neurogenic ectoderm without a canonical Simbinding site. Intriguingly, the impaired transcriptional synergy between Dl and either Zld or Bcd led to aberrant and abolished lacZ expression in the neurogenic ectoderm and in the ventral midline, respectively. These findings suggest that the two enhancer activities of the sog primary enhancer are functionally consolidated and geographically inseparable. [BMB Reports 2016; 49(10): 572-577]

Midline enhancer activity of the short gastrulation shadow enhancer is characterized by three unusual features for cis-regulatory DNA

The shadow enhancer of the short gastrulation (sog) gene directs its sequential expression in the neurogenic ectoderm and the ventral midline of the developing Drosophila embryo. Here, we characterize three unusual features of the shadow enhancer midline activity. First, the minimal regions for the two different enhancer activities exhibit high overlap within the shadow enhancer, meaning that one developmental enhancer possesses dual enhancer activities. Second, the midline enhancer activity relies on five Single-minded (Sim)-binding sites, two of which have not been found in any Sim target enhancers. Finally, two linked Dorsal (Dl)- and Zelda (Zld)-binding sites, critical for the neurogenic ectoderm enhancer activity, are also required for the midline enhancer activity. These results suggest that early activation by Dl and Zld may facilitate late activation via the noncanonical sites occupied by Sim. We discuss a model for Zld as a pioneer factor and speculate its role in midline enhancer activity. [BMB Reports 2015; 48(10): 589-594]

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

DOI:http://dx.doi.org/10.7554/eLife.08445.001

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

DOI:http://dx.doi.org/10.7554/eLife.08445.002

Capicua is involved in Dorsal-mediated repression of zerknüllt expression in Drosophila embryo

The maternal transcription factor Dorsal (Dl) functions as both an activator and a repressor in a context-dependent manner to control dorsal-ventral patterning in the Drosophila embryo. Previous studies have suggested that Dl is an intrinsic activator and its repressive activity requires additional corepressors that bind corepressor-binding sites near Dl-binding sites. However, the molecular identities of the corepressors have yet to be identified. Here, we present evidence that Capicua (Cic) is involved in Dl-mediated repression in the zerknüllt (zen) ventral repression element (VRE). Computational and genetic analyses indicate that a DNA-binding consensus sequence of Cic is highly analogous with previously identified corepressorbinding sequences and that Dl failed to repress zen expression in lateral regions of cic mutant embryos. Furthermore, electrophoretic mobility shift assay (EMSA) shows that Cic directly interacts with several corepressor-binding sites in the zen VRE. These results suggest that Cic may function as a corepressor by binding the VRE. [BMB Reports 2014; 47(9):518-523]

Co-activation of microRNAs by Zelda is essential for early Drosophila development

Transcription factors and microRNAs (miRNAs) are two important classes of trans-regulators in differential gene expression. Transcription factors occupy cis-regulatory motifs in DNA to activate or repress gene transcription, whereas miRNAs specifically pair with seed sites in target mRNAs to trigger mRNA decay or inhibit translation. Dynamic spatiotemporal expression patterns of transcription factors and miRNAs during development point to their stage- and tissue-specific functions. Recent studies have focused on miRNA functions during development; however, much remains to explore regarding how the expression of miRNAs is initiated and how dynamic miRNA expression patterns are achieved by transcriptional regulatory networks at different developmental stages. Here, we focused on the identification, regulation and function of miRNAs during the earliest stage of Drosophila development, when the maternal-to-zygotic transition (MZT) takes place. Eleven miRNA clusters comprise the first set of miRNAs activated in the blastoderm embryo. The transcriptional activator Zelda is required for their proper activation and regulation, and Zelda binding observed in genome-wide binding profiles is predictive of enhancer activity. In addition, other blastoderm transcription factors, comprising both activators and repressors, the activities of which are potentiated and coordinated by Zelda, contribute to the accurate temporal and spatial expression of these miRNAs, which are known to function in diverse developmental processes. Although previous genetic studies showed no early phenotypes upon loss of individual miRNAs, our analysis of the miR-1; miR-9a double mutant revealed defects in gastrulation, demonstrating the importance of co-activation of miRNAs by Zelda during the MZT.

Context-dependent transcriptional interpretation of mitogen activated protein kinase signaling in the Drosophila embryo

Terminal regions of the Drosophila embryo are patterned by the localized activation of Mitogen Activated Protein Kinase (MAPK), which induces zygotic genes through relief of their repression by transcriptional repressor Capicua. The levels of MAPK activation at the anterior and posterior termini are close to each other, but the expression patterns of MAPK-target genes, such as zerknüllt (zen) and tailless (tll), display strong anterior-posterior (AP) asymmetry. This region-specific response to MAPK activation provides a clear example of context-dependent interpretation of inductive signaling, a common developmental effect that remains poorly understood. In the past, the AP asymmetry of zen expression was attributed to a mechanism that depends on MAPK substrate competition. We present data suggesting that the asymmetric expression of tll is generated by a different mechanism, based on feedforward control and multiple enhancers of the tll gene. A simple mathematical model of this mechanism correctly predicts how the wild-type expression pattern of tll changes in mutants affecting the anterior, dorsoventral, and terminal patterning systems and some of their direct targets.

Twine: display and analysis of cis-regulatory modules

Many algorithms analyze enhancers for overrepresentation of known and novel motifs, with the goal of identifying binding sites for direct regulators of gene expression. Twine is a Java GUI with multiple graphical representations ('Views') of enhancer alignments that displays motifs, as IUPAC consensus sequences or position frequency matrices, in the context of phylogenetic conservation to facilitate cis-regulatory element discovery. Thresholds of phylogenetic conservation and motif stringency can be altered dynamically to facilitate detailed analysis of enhancer architecture. Views can be exported to vector graphics programs to generate high-quality figures for publication. Twine can be extended via Java plugins to manipulate alignments and analyze sequences.

AVAILABILITY

Twine is freely available as a compiled Java .jar package or Java source code at http://labs.bio.unc.edu/crews/twine/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Deep sequencing and expression of microRNAs from early honeybee (Apis mellifera) embryos reveals a role in regulating early embryonic patterning

Background

Recent evidence supports the proposal that the observed diversity of animal body plans has been produced through alterations to the complexity of the regulatory genome rather than increases in the protein-coding content of a genome. One significant form of gene regulation is the contribution made by the non-coding content of the genome. Non-coding RNAs play roles in embryonic development of animals and these functions might be expected to evolve rapidly. Using next-generation sequencing and in situ hybridization, we have examined the miRNA content of early honeybee embryos.

Results

Through small RNA sequencing we found that 28% of known miRNAs are expressed in the early embryo. We also identified developmentally expressed microRNAs that are unique to the Apoidea clade. Examination of expression patterns implied these miRNAs have roles in patterning the anterior-posterior and dorso-ventral axes as well as the extraembryonic membranes. Knockdown of Dicer, a key component of miRNA processing, confirmed that miRNAs are likely to have a role in patterning these tissues.

Conclusions

Examination of the expression patterns of novel miRNAs, some unique to the Apis group, indicated that they are likely to play a role in early honeybee development. Known miRNAs that are deeply conserved in animal phyla display differences in expression pattern between honeybee and Drosophila, particularly at early stages of development. This may indicate miRNAs play a rapidly evolving role in regulating developmental pathways, most likely through changes to the way their expression is regulated.

Response to the BMP gradient requires highly combinatorial inputs from multiple patterning systems in the Drosophila embryo

Pattern formation in the developing embryo relies on key regulatory molecules, many of which are distributed in concentration gradients. For example, a gradient of BMP specifies cell fates along the dorsoventral axis in species ranging from flies to mammals. In Drosophila, a gradient of the BMP molecule Dpp gives rise to nested domains of target gene expression in the dorsal region of the embryo; however, the mechanisms underlying the differential response are not well understood, partly owing to an insufficient number of well-studied targets. Here we analyze how the Dpp gradient regulates expression of pannier (pnr), a candidate low-level Dpp target gene. We predicted that the pnr enhancer would contain high-affinity binding sites for the Dpp effector Smad transcription factors, which would be occupied in the presence of low-level Dpp. Unexpectedly, the affinity of Smad sites in the pnr enhancer was similar to those in the Race enhancer, a high-level Dpp target gene, suggesting that the affinity threshold mechanism plays a minimal role in the regulation of pnr. Our results indicate that a mechanism involving a conserved bipartite motif that is predicted to bind a homeodomain factor in addition to Smads and the Brinker repressor, establishes the pnr expression domain. Furthermore, the pnr enhancer has a highly complex structure that integrates cues not only from the dorsoventral axis, but also from the anteroposterior and terminal patterning systems in the blastoderm embryo.

Huckebein is part of a combinatorial repression code in the anterior blastoderm

The hierarchy of the segmentation cascade responsible for establishing the Drosophila body plan is composed by gap, pair-rule and segment polarity genes. However, no pair-rule stripes are formed in the anterior regions of the embryo. This lack of stripe formation, as well as other evidence from the literature that is further investigated here, led us to the hypothesis that anterior gap genes might be involved in a combinatorial mechanism responsible for repressing the cis-regulatory modules (CRMs) of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) anterior-most stripes. In this study, we investigated huckebein (hkb), which has a gap expression domain at the anterior tip of the embryo. Using genetic methods we were able to detect deviations from the wild-type patterns of the anterior-most pair-rule stripes in different genetic backgrounds, which were consistent with Hkb-mediated repression. Moreover, we developed an image processing tool that, for the most part, confirmed our assumptions. Using an hkb misexpression system, we further detected specific repression on anterior stripes. Furthermore, bioinformatics analysis predicted an increased significance of binding site clusters in the CRMs of h 1, eve 1, run 1 and ftz 1when Hkb was incorporated in the analysis, indicating that Hkb plays a direct role in these CRMs. We further discuss that Hkb and Slp1, which is the other previously identified common repressor of anterior stripes, might participate in a combinatorial repression mechanism controlling stripe CRMs in the anterior parts of the embryo and define the borders of these anterior stripes.

Multiple enhancers ensure precision of gap gene-expression patterns in the Drosophila embryo

Segmentation of the Drosophila embryo begins with the establishment of spatially restricted gap gene-expression patterns in response to broad gradients of maternal transcription factors, such as Bicoid. Numerous studies have documented the fidelity of these expression patterns, even when embryos are subjected to genetic or environmental stress, but the underlying mechanisms for this transcriptional precision are uncertain. Here we present evidence that every gap gene contains multiple enhancers with overlapping activities to produce authentic patterns of gene expression. For example, a recently identified hunchback (hb) enhancer (located 5-kb upstream of the classic enhancer) ensures repression at the anterior pole. The combination of intronic and 5' knirps (kni) enhancers produces a faithful expression pattern, even though the intronic enhancer alone directs an abnormally broad expression pattern. We present different models for "enhancer synergy," whereby two enhancers with overlapping activities produce authentic patterns of gene expression.

Diversity in insect axis formation: two orthodenticle genes and hunchback act in anterior patterning and influence dorsoventral organization in the honeybee (Apis mellifera)

Axis formation is a key step in development, but studies indicate that genes involved in insect axis formation are relatively fast evolving. Orthodenticle genes have conserved roles, often with hunchback, in maternal anterior patterning in several insect species. We show that two orthodenticle genes, otd1 and otd2, and hunchback act as maternal anterior patterning genes in the honeybee (Apis mellifera) but, unlike other insects, act to pattern the majority of the anteroposterior axis. These genes regulate the expression domains of anterior, central and posterior gap genes and may directly regulate the anterior gap gene giant. We show otd1 and hunchback also influence dorsoventral patterning by regulating zerknült (zen) as they do in Tribolium, but that zen does not regulate the expression of honeybee gap genes. This suggests that interactions between anteroposterior and dorsal-ventral patterning are ancestral in holometabolous insects. Honeybee axis formation, and the function of the conserved anterior patterning gene orthodenticle, displays unique characters that indicate that, even when conserved genes pattern the axis, their regulatory interactions differ within orders of insects, consistent with relatively fast evolution in axis formation pathways.

Complex interplay of evolutionary forces in the ladybird homeobox genes of Drosophila melanogaster

Tandemly arranged paralogous genes lbe and lbl are members of the Drosophila NK homeobox family. We analyzed population samples of Drosophila melanogaster from Africa, Europe, North and South America, and single strains of D. sechellia, D. simulans, and D. yakuba within two linked regions encompassing partial sequences of lbe and lbl. The evolution of lbe and lbl is highly constrained due to their important regulatory functions. Despite this, a variety of forces have shaped the patterns of variation in lb genes: recombination, intragenic gene conversion and natural selection strongly influence background variation created by linkage disequilibrium and dimorphic haplotype structure. The two genes exhibited similar levels of nucleotide diversity and positive selection was detected in the noncoding regions of both genes. However, synonymous variability was significantly higher for lbe: no nonsynonymous changes were observed in this gene. We argue that balancing selection impacts some synonymous sites of the lbe gene. Stability of mRNA secondary structure was significantly different between the lbe (but not lbl) haplotype groups and may represent a driving force of balancing selection in epistatically interacting synonymous sites. Balancing selection on synonymous sites may be the first, or one of a few such observations, in Drosophila. In contrast, recurrent positive selection on lbl at the protein level influenced evolution at three codon sites. Transcription factor binding-site profiles were different for lbe and lbl, suggesting that their developmental functions are not redundant. Combined with our previous results on nucleotide variation in esterase and other homeobox genes, these results suggest that interplay of balancing and directional selection may be a general feature of molecular evolution in Drosophila and other eukaryote genomes.

The evolution of dorsal-ventral patterning mechanisms in insects

The gene regulatory network (GRN) underpinning dorsal-ventral (DV) patterning of the Drosophila embryo is among the most thoroughly understood GRNs, making it an ideal system for comparative studies seeking to understand the evolution of development. With the emergence of widely applicable techniques for testing gene function, species with sequenced genomes, and multiple tractable species with diverse developmental modes, a phylogenetically broad and molecularly deep understanding of the evolution of DV axis formation in insects is feasible. Here, we review recent progress made in this field, compare our emerging molecular understanding to classical embryological experiments, and suggest future directions of inquiry.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Organization of developmental enhancers in the Drosophila embryo

Most cell-specific enhancers are thought to lack an inherent organization, with critical binding sites distributed in a more or less random fashion. However, there are examples of fixed arrangements of binding sites, such as helical phasing, that promote the formation of higher-order protein complexes on the enhancer DNA template. Here, we investigate the regulatory ‘grammar’ of nearly 100 characterized enhancers for developmental control genes active in the early Drosophila embryo. The conservation of grammar is examined in seven divergent Drosophila genomes. Linked binding sites are observed for particular combinations of binding motifs, including Bicoid–Bicoid, Hunchback–Hunchback, Bicoid–Dorsal, Bicoid–Caudal and Dorsal–Twist. Direct evidence is presented for the importance of Bicoid–Dorsal linkage in the integration of the anterior–posterior and dorsal–ventral patterning systems. Hunchback–Hunchback interactions help explain unresolved aspects of segmentation, including the differential regulation of the eve stripe 3 + 7 and stripe 4 + 6 enhancers. We also present evidence that there is an under-representation of nucleosome positioning sequences in many enhancers, raising the possibility for a subtle higher-order structure extending across certain enhancers. We conclude that grammar of gene control regions is pervasively used in the patterning of the Drosophila embryo.

Statistical detection of cooperative transcription factors with similarity adjustment

Motivation: Statistical assessment of cis-regulatory modules (CRMs) is a crucial task in computational biology. Usually, one concludes from exceptional co-occurrences of DNA motifs that the corresponding transcription factors (TFs) are cooperative. However, similar DNA motifs tend to co-occur in random sequences due to high probability of overlapping occurrences. Therefore, it is important to consider similarity of DNA motifs in the statistical assessment.

Results: Based on previous work, we propose to adjust the window size for co-occurrence detection. Using the derived approximation, one obtains different window sizes for different sets of DNA motifs depending on their similarities. This ensures that the probability of co-occurrences in random sequences are equal. Applying the approach to selected similar and dissimilar DNA motifs from human TFs shows the necessity of adjustment and confirms the accuracy of the approximation by comparison to simulated data. Furthermore, it becomes clear that approaches ignoring similarities strongly underestimate P-values for cooperativity of TFs with similar DNA motifs. In addition, the approach is extended to deal with overlapping windows. We derive Chen–Stein error bounds for the approximation. Comparing the error bounds for similar and dissimilar DNA motifs shows that the approximation for similar DNA motifs yields large bounds. Hence, one has to be careful using overlapping windows. Based on the error bounds, one can precompute the approximation errors and select an appropriate overlap scheme before running the analysis.

Availability: Software to perform the calculation for pairs of position frequency matrices (PFMs) is available at http://mosta.molgen.mpg.de as well as C++ source code for downloading.

Contact:utz.pape@molgen.mpg.de

Finding evolutionarily conserved cis-regulatory modules with a universal set of motifs

BACKGROUND

Finding functional regulatory elements in DNA sequences is a very important problem in computational biology and providing a reliable algorithm for this task would be a major step towards understanding regulatory mechanisms on genome-wide scale. Major obstacles in this respect are that the fact that the amount of non-coding DNA is vast, and that the methods for predicting functional transcription factor binding sites tend to produce results with a high percentage of false positives. This makes the problem of finding regions significantly enriched in binding sites difficult.

RESULTS

We develop a novel method for predicting regulatory regions in DNA sequences, which is designed to exploit the evolutionary conservation of regulatory elements between species without assuming that the order of motifs is preserved across species. We have implemented our method and tested its predictive abilities on various datasets from different organisms.

CONCLUSION

We show that our approach enables us to find a majority of the known CRMs using only sequence information from different species together with currently publicly available motif data. Also, our method is robust enough to perform well in predicting CRMs, despite differences in tissue specificity and even across species, provided that the evolutionary distances between compared species do not change substantially. The complexity of the proposed algorithm is polynomial, and the observed running times show that it may be readily applied.

MotifCluster: an interactive online tool for clustering and visualizing sequences using shared motifs

MotifCluster finds related motifs in a set of sequences and clusters the sequences into families using the motifs they contain.

MotifCluster finds related motifs in a set of sequences, and clusters the sequences into families using the motifs they contain. MotifCluster, at , lets users test whether proteins are related, cluster sequences by shared conserved motifs, and visualize motifs mapped onto trees, sequences and three-dimensional structures. We demonstrate MotifCluster's accuracy using gold-standard protein superfamilies; using recommended settings, families were assigned to the correct superfamilies with 0.17% false positive and no false negative assignments.

Identification of tissue-specific cis-regulatory modules based on interactions between transcription factors

BACKGROUND

Evolutionary conservation has been used successfully to help identify cis-acting DNA regions that are important in regulating tissue-specific gene expression. Motivated by increasing evidence that some DNA regulatory regions are not evolutionary conserved, we have developed an approach for cis-regulatory region identification that does not rely upon evolutionary sequence conservation.

RESULTS

The conservation-independent approach is based on an empirical potential energy between interacting transcription factors (TFs). In this analysis, the potential energy is defined as a function of the number of TF interactions in a genomic region and the strength of the interactions. By identifying sets of interacting TFs, the analysis locates regions enriched with the binding sites of these interacting TFs. We applied this approach to 30 human tissues and identified 6232 putative cis-regulatory modules (CRMs) regulating 2130 tissue-specific genes. Interestingly, some genes appear to be regulated by different CRMs in different tissues. Known regulatory regions are highly enriched in our predicted CRMs. In addition, DNase I hypersensitive sites, which tend to be associated with active regulatory regions, significantly overlap with the predicted CRMs, but not with more conserved regions. We also find that conserved and non-conserved CRMs regulate distinct gene groups. Conserved CRMs control more essential genes and genes involved in fundamental cellular activities such as transcription. In contrast, non-conserved CRMs, in general, regulate more non-essential genes, such as genes related to neural activity.

CONCLUSION

These results demonstrate that identifying relevant sets of binding motifs can help in the mapping of DNA regulatory regions, and suggest that non-conserved CRMs play an important role in gene regulation.

Exact p-value calculation for heterotypic clusters of regulatory motifs and its application in computational annotation of cis-regulatory modules

Background

cis-Regulatory modules (CRMs) of eukaryotic genes often contain multiple binding sites for transcription factors. The phenomenon that binding sites form clusters in CRMs is exploited in many algorithms to locate CRMs in a genome. This gives rise to the problem of calculating the statistical significance of the event that multiple sites, recognized by different factors, would be found simultaneously in a text of a fixed length. The main difficulty comes from overlapping occurrences of motifs. So far, no tools have been developed allowing the computation of p-values for simultaneous occurrences of different motifs which can overlap.

Results

We developed and implemented an algorithm computing the p-value that s different motifs occur respectively k₁, ..., k_s or more times, possibly overlapping, in a random text. Motifs can be represented with a majority of popular motif models, but in all cases, without indels. Zero or first order Markov chains can be adopted as a model for the random text. The computational tool was tested on the set of cis-regulatory modules involved in D. melanogaster early development, for which there exists an annotation of binding sites for transcription factors. Our test allowed us to correctly identify transcription factors cooperatively/competitively binding to DNA.

Method

The algorithm that precisely computes the probability of simultaneous motif occurrences is inspired by the Aho-Corasick automaton and employs a prefix tree together with a transition function. The algorithm runs with the O(n|Σ|(m|ℋ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFlecsaaa@3762@| + K|σ|^K) ∏_i k_i) time complexity, where n is the length of the text, |Σ| is the alphabet size, m is the maximal motif length, |ℋ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFlecsaaa@3762@| is the total number of words in motifs, K is the order of Markov model, and k_i is the number of occurrences of the ith motif.

Conclusion

The primary objective of the program is to assess the likelihood that a given DNA segment is CRM regulated with a known set of regulatory factors. In addition, the program can also be used to select the appropriate threshold for PWM scanning. Another application is assessing similarity of different motifs.

Availability

Project web page, stand-alone version and documentation can be found at