Predicting regulons and their cis-regulatory motifs by comparative genomics

We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caernorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.     

Combining phylogenetic data with co-regulated genes to identify regulatory motifs

MOTIVATION: Discovery of regulatory motifs in unaligned DNA sequences remains a fundamental problem in computational biology. Two categories of algorithms have been developed to identify common motifs from a set of DNA sequences. The first can be called a 'multiple genes, single species' approach. It proposes that a degenerate motif is embedded in some or all of the otherwise unrelated input sequences and tries to describe a consensus motif and identify its occurrences. It is often used for co-regulated genes identified through experimental approaches. The second approach can be called 'single gene, multiple species'. It requires orthologous input sequences and tries to identify unusually well conserved regions by phylogenetic footprinting. Both approaches perform well, but each has some limitations. It is tempting to combine the knowledge of co-regulation among different genes and conservation among orthologous genes to improve our ability to identify motifs. RESULTS: Based on the Consensus algorithm previously established by our group, we introduce a new algorithm called PhyloCon (Phylogenetic Consensus) that takes into account both conservation among orthologous genes and co-regulation of genes within a species. This algorithm first aligns conserved regions of orthologous sequences into multiple sequence alignments, or profiles, then compares profiles representing non-orthologous sequences. Motifs emerge as common regions in these profiles. Here we present a novel statistic to compare profiles of DNA sequences and a greedy approach to search for common subprofiles. We demonstrate that PhyloCon performs well on both synthetic and biological data. AVAILABILITY: Software available upon request from the authors. http://ural.wustl.edu/softwares.html     

酵母基因上游与内含子可能存在的转录协同作用

酵母基因上游与内含子可能存在的转录协同作用     

Annotating regulatory DNA based on man-mouse genomic comparison

Non-coding DNA segments that are conserved between the human and mouse genomic sequence are good indicators of possible regulatory sequences. Here we report on a systematic approach to delineate such conserved elements from upstream regions of orthologous gene pairs from man and mouse. We focus on orthologous genes in order to maximize our chances to find functionally similar regulatory elements. The identification of conserved elements is effected using the Waterman-Eggert local suboptimal alignment algorithm. We have modified an implementation of this algorithm such that it integrates the determination of statistical significance for the local suboptimal alignments. This has the effect of outputting a dynamically determined number of suboptimal alignments that are deemed statistically significant. Comparison with experimentally determined annotation shows a striking enrichement of regulatory sites among the conserved regions. Furthermore, the conserved regions tend to cover the promotor region described in the EPD database.     

Characterization of pLP18, a novel cryptic plasmid of Lactobacillus plantarum PC518 isolated from Chinese pickle

A cryptic plasmid of Lactobacillus plantarum PC518 isolated from Chinese pickle, designated pLP18, was sequenced and characterized. It is a 1806-bp circular molecule with a G+C content of 37.5%. Sequence analysis of pLP18 revealed three putative open reading frames (ORFs), in which ORF1 contained conserved motifs of pMV158-family Rep proteins and showed 60% similarity with the Rep protein of pPSC22, a member of rolling-circle replication (RCR) pMV158 family. The double strand origin (dso) of pMV158 family and the single strand origin A (ssoA) located upstream of the rep gene. The putative cop and rnaII genes were predicted to be regulatory genes controlling copy number of pLP18. The results of Southern hybridization suggested that pLP18 replicate via the RCR mechanism. Furthermore, the relative copy number of pLP18 was estimated to be about 24 copies per chromosome equivalent by quantitative PCR.     

Algorithms for phylogenetic footprinting.

Phylogenetic footprinting is a technique that identifies regulatory elements by finding unusually well conserved regions in a set of orthologous noncoding DNA sequences from multiple species. We introduce a new motif-finding problem, the Substring Parsimony Problem, which is a formalization of the ideas behind phylogenetic footprinting, and we present an exact dynamic programming algorithm to solve it. We then present a number of algorithmic optimizations that allow our program to run quickly on most biologically interesting datasets. We show how to handle data sets in which only an unknown subset of the sequences contains the regulatory element. Finally, we describe how to empirically assess the statistical significance of the motifs found. Each technique is implemented and successfully identifies a number of known binding sites, as well as several highly conserved but uncharacterized regions. The program is available at http://bio.cs.washington.edu/software.html.     

Pigment cell-specific expression of the tyrosinase gene in ascidians has a different regulatory mechanism from vertebrates

Tyrosinase is the key enzyme required for the synthesis of melanin pigments. Sequence comparison and functional analysis of the 5' upstream regions of vertebrate tyrosinase genes have revealed the importance of conserved E-box motifs in regulating their specific expression in pigment cells, optic cup-derived retinal pigment epithelium (RPE) and neural crest-derived melanocytes. In ascidians (more basal protochordates), two pigment cells that resemble vertebrate RPE cells are formed and specifically express the orthologous tyrosinase gene (HrTyr) in the cerebral vesicle located at the anterior end of the neural tube. To define regulatory sequences required for pigment cell-lineage-specific expression of HrTyr during embryogenesis, a series of mutations of the 5' upstream region of HrTyr were fused to the lacZ reporter gene and were microinjected into fertilized eggs. We found that the -152bp upstream of the translational start site is essential for expression in pigment cell precursors of tailbud-stage embryos. Further, additional positive and unique restriction elements were identified in the region up to -1.8kb. Surprisingly, in the -152bp minimal promoter or in other regions with regulatory activities, there are no E-box motifs or sequences correlating with other conserved elements regulating vertebrate tyrosinase promoters. The possibility that Pax proteins regulate HrTyr expression is also discussed.     

High sequence turnover in the regulatory regions of the developmental gene hunchback in insects

Extensive sequence analysis of the developmental gene hunchback and its 5' and 3' regulatory regions in Drosophila melanogaster, Drosophila virilis, Musca domestica, and Tribolium castaneum, using a variety of computer algorithms, reveals regions of high sequence simplicity probably generated by slippage-like mechanisms of turnover. No regions are entirely refractory to the action of slippage, although the density and composition of simple sequence motifs varies from region to region. Interestingly, the 5' and 3' flanking regions share short repetitive motifs despite their separation by the gene itself, and the motifs are different in composition from those in the exons and introns. Furthermore, there are high levels of conservation of motifs in equivalent orthologous regions. Detailed sequence analysis of the P2 promoter and DNA footprinting assays reveal that the number, orientation, sequence, spacing, and protein-binding affinities of the BICOID-binding sites varies between species and that the 'P2' promoter, the nanos response element in the 3' untranslated region, and several conserved boxes of sequence in the gene (e.g., the two zinc-finger regions) are surrounded by cryptically-simple-sequence DNA. We argue that high sequence turnover and genetic redundancy permit both the general maintenance of promoter functions through the establishment of coevolutionary (compensatory) changes in cis- and trans-acting genetic elements and, at the same time, the possibility of subtle changes in the regulation of hunchback in the different species.