Identification of consensus patterns in unaligned DNA sequences known to be functionally related

We have developed a method for identifying consensus patternsin a set of unaligned DNA sequences known to bind a common proteinor to have some other common biochemical function. The methodis based on a tnatrix representation of binding site patterns.Each row of the matrix represents one of the four possible bases,each column represents one of the positions of the binding siteand each element is determined by the frequency the indicatedbase occurs at the indicated position. The goal of the methodis to find the most significant matrix-i.e. the one with thelowest probability of occurring by chance-out of all the matricesthat can be formed from the set of related sequences. The reliabilityof the method improves with the number of sequences, while thetime required increases only linearly with the number of sequences.To test this method, we analysed 11 DNA sequences containingpromoters regulated by the Escherichia coli LexA protein. Thematrices we' found were consistent with the known consensussequence, and could distinguish the generally accepted LexAbinding sites from other DNA sequences. Received on November 6, 1989; accepted on December 20, 1989     

Translocation in yeast and mammalian cells: not all signal sequences are functionally equivalent

In Saccharomyces cerevisiae, nascent carboxypeptidase Y (CPY) is directed into the endoplasmic reticulum by an NH2-terminal signal peptide that is removed before the glycosylated protein is transported to the vacuole. In this paper, we show that this signal peptide does not function in mammalian cells: CPY expressed in COS-1 cells is not glycosylated, does not associate with membranes, and retains its signal peptide. In a mammalian cell-free protein-synthesizing system, CPY is not translocated into microsomes. However, if the CPY signal is either mutated to increase its hydrophobicity or replaced with that of influenza virus hemagglutinin, the resulting precursors are efficiently translocated both in vivo and in vitro. The implications of these results for models of signal sequence function are discussed.     

Finding composite regulatory patterns in DNA sequences

Pattern discovery in unaligned DNA sequences is a fundamental problem in computational biology with important applications in finding regulatory signals. Current approaches to pattern discovery focus on monad patterns that correspond to relatively short contiguous strings. However, many of the actual regulatory signals are composite patterns that are groups of monad patterns that occur near each other. A difficulty in discovering composite patterns is that one or both of the component monad patterns in the group may be 'too weak'. Since the traditional monad-based motif finding algorithms usually output one (or a few) high scoring patterns, they often fail to find composite regulatory signals consisting of weak monad parts. In this paper, we present a MITRA (MIsmatch TRee Algorithm) approach for discovering composite signals. We demonstrate that MITRA performs well for both monad and composite patterns by presenting experiments over biological and synthetic data.     

Recognition of specific DNA sequences in eukaryotic chromosomes.

The packaging of DNA into chromatin probably places certain restrictions on how specific DNA sequences can be recognized by DNA sequence specific recognition proteins (SRP). Several unique features of this type of interaction are discussed. Specifically, as a consequence of the coiling of the DNA about a histone core, it is proposed that DNA recognition sites will be compound and that each element of the compound recognition site will be about 10 - 20 b.p. in length and distributed at approximately 80 b.p. intervals--the periodicity of the DNA wrapping around the nucleosome.     

A phyletic perspective on the allometry of plant biomass-partitioning patterns and functionally equivalent organ-categories

Biomass-partitioning patterns influence the functioning of aquatic and terrestrial vegetation at all levels, ranging from individual growth and reproduction to the flow of mass and energy through entire communities. For this reason, leaf, stem and root dry biomass-partitioning patterns across taxonomically and ecologically diverse seed plants (spermatophytes) have been intensively investigated, both empirically and theoretically. By contrast, phyletically disparate plants (e.g. green and brown algal macrophytes, mosses and pteridophytes) have not been examined to determine whether the partitioning of their body parts into 'leaf', 'stem' and 'root' analogs accords with that of spermatophytes. In this review, the biomass-partitioning patterns of siphonous and brown algal macrophytes, mosses and pteridophytes were compared allometrically with those of spermatophytes and were shown to be largely in statistical accordance (thus lending support to the hypothesis that a single scaling relationship exists across eukaryotic photoautotrophs). This concordance is argued to support the hypothesis of functional equivalence across analogous, but developmentally different, body parts, a feature that permits the use of simpler biological model systems with which to derive analytical explanations for the biomass-partitioning patterns reported for more complex seed plants.     

Hierarchical classification of functionally equivalent genes in prokaryotes

Functional classification of genes represents a fundamental problem to many biological studies. Most of the existing classification schemes are based on the concepts of homology and orthology, which were originally introduced to study gene evolution but might not be the most appropriate for gene function prediction, particularly at high resolution level. We have recently developed a scheme for hierarchical classification of genes (HCGs) in prokaryotes. In the HCG scheme, the functional equivalence relationships among genes are first assessed through a careful application of both sequence similarity and genomic neighborhood information; and genes are then classified into a hierarchical structure of clusters, where genes in each cluster are functionally equivalent at some resolution level, and the level of resolution goes higher as the clusters become increasingly smaller traveling down the hierarchy. The HCG scheme is validated through comparisons with the taxonomy of the prokaryotic genomes, Clusters of Orthologous Groups (COGs) of genes and the Pfam system. We have applied the HCG scheme to 224 complete prokaryotic genomes, and constructed a HCG database consisting of a forest of 5339 multi-level and 15 770 single-level trees of gene clusters covering approximately 93% of the genes of these 224 genomes. The validation results indicate that the HCG scheme not only captures the key features of the existing classification schemes but also provides a much richer organization of genes which can be used for functional prediction of genes at higher resolution and to help reveal evolutionary trace of the genes.     

Repetitive plasmid sequences generate DNA fingerprinting patterns in mammals

Bacterial plasmids with stringently regulated copy numbers have directly repeated DNA sequences, termed iterons, in the vicinity of their replication origins. These sequences bind a specific protein exerting a key role in the initiation of plasmid replication. Plasmids P1, pSC101 and RFS1010 have different iteron sequences and belong to three different incompatibility groups. Used as DNA probes each of these plasmids generates specific patterns in mammals similar to those obtained by the DNA fingerprinting technique. The iteron-containing regions were identified as the part of the plasmids responsible for those patterns by using polymerase chain reaction (PCR) amplified DNA segments that contained the iteron regions as probes.     

Recognition of DNA sequences by the repressor of bacteriophage 434

The structure of a complex between the DNA-binding domain of phage 434 repressor and a 14 base-pair synthetic DNA operator reveals the molecular interactions important for sequence-specific recognition. A set of contacts with DNA backbone, notably involving hydrogen bonds between peptide-NH groups and DNA phosphates, position the repressor and fix the DNA configuration. Direct interactions between amino acid side chains and DNA bases involve nonpolar van der Waals contacts as well as hydrogen bonds. The structures of the repressor domain and of the 434 cro protein are extremely similar. There appear to be no major conformational changes in the proteins when they bind to DNA.     

Recognition of specific DNA sequences by mitomycin C for alkylation.

Synthetic oligodeoxyribonucleotides were reacted with mitomycin C (MC) under conditions which restricted MC to monofunctional alkylating activity. The yields of monofunctional alkylation of oligonucleotides with variable sequence were determined by enzymatic digestion of the reaction mixture to unreacted nucleosides and the product of alkylation, a MC-deoxyguanosine adduct (2), followed by quantitative analysis by HPLC. The relative yields of 2 reflected relative monoalkylation reactivities. They were compared in a series of oligonucleotides having the sequence 5'-NGN' in which the 5'-base was varied while the 3'-base was kept constant as T. Under Na2S2O4 activation conditions a striking enhancement of the yield was observed at the 5'-CG sequence: 36%, compared to 2% at 5'-AG and 4.1% at 5'-TG. The 5'-GG sequence also showed enhanced reactivity although to a lesser extent (14.7%). The enhancements were specific to the duplex state of the oligonucleotides. Using NADPH:cytochrome c reductase as the reducing agent gave similar results. MC activated by acidic pH also displayed 5'-CG alkylation specificity. 10-Decarbamoyl-MC activated by Na2S2O4 showed the same 5'-CG specificity as MC. Replacement of deoxyguanosine by deoxyinosine in the opposite strand at a 5'-CG site abolished the enhancement of alkylation. Such replacement at a 5'-GG site had a similar effect. It was found that the base 3' to the guanine had only a relatively modest modulating effect on the enhanced reactivity of the G at the 5'-CG sequence. This 3'-base effect appeared to be independent of the 5'-base of the 5'-NGN' triplet. The order of reactivity is 3'-(C greater than T greater than A).(ABSTRACT TRUNCATED AT 250 WORDS)     

2-D graphical representation for characteristic sequences of DNA and its application

DNA sequencing has resulted in an abundance of data on DNA sequences for various species. Hence, the characterization and comparison of sequences become more important but still difficult tasks. In this paper, we first give a 2-D ladderlike graphical representation for the characteristic sequences of a DNA sequence, and then construct a 3-component vector, in which the normalized ALE-indices extracted from such three 2-D graphs via D/D matrices are individual components, to characterize the DNA sequence. The examination of similarities/dissimilarities among sequences of the beta-globin genes of different species illustrates the utility of the approach.     

Thermodynamic properties of DNA sequences: characteristic values for the human genome

MOTIVATION: Central to many molecular biology techniques as ubiquitous as PCR and Southern blotting is the design of oligonucleotide (oligo) probes and/or primers possessing specific thermodynamic properties. Here, we use validated theoretical methods to generate distributions of predicted thermodynamic properties for DNA oligos of various lengths. These distributions facilitate immediate appreciation of typical thermodynamic values for oligos of various lengths. RESULTS: Distributions of melting temperature (Tm), free energy (DeltaG(T)o), and fraction hybridized or fraction bound (Fb), are presented for oligos of length 10-50 bases sampled from the human genome. The effects of changing temperature, oligo and salt concentrations, constraining G+C content, and introducing mismatches are exemplified. Our results provide the first survey of typical and limiting thermodynamic values evaluated on a genomic scale. Described numbers comprise useful 'rules of thumb' that are applicable to most technologies dependent upon DNA oligo design.     

Spatial,temporal and interindividual epigenetic variation of functionally important DNA methylation patterns

DNA methylation is an epigenetic modification that plays an important role in gene regulation. It can be influenced by stochastic events, environmental factors and developmental programs. However, little is known about the natural variation of gene-specific methylation patterns. In this study, we performed quantitative methylation analyses of six differentially methylated imprinted genes (H19, MEG3, LIT1, NESP55, PEG3 and SNRPN), one hypermethylated pluripotency gene (OCT4) and one hypomethylated tumor suppressor gene (APC) in chorionic villus, fetal and adult cortex, and adult blood samples. Both average methylation level and range of methylation variation depended on the gene locus, tissue type and/or developmental stage. We found considerable variability of functionally important methylation patterns among unrelated healthy individuals and a trend toward more similar methylation levels in monozygotic twins than in dizygotic twins. Imprinted genes showed relatively little methylation changes associated with aging in individuals who are >25 years. The relative differences in methylation among neighboring CpGs in the generally hypomethylated APC promoter may not only reflect stochastic fluctuations but also depend on the tissue type. Our results are consistent with the view that most methylation variation may arise after fertilization, leading to epigenetic mosaicism.     

Virulence-related DNA sequences and adherence patterns in strains of enteropathogenic Escherichia coli

The presence of the genes for Escherichia coli adherence factor (EAF), attaching and effacing lesion (eae) and bundle-forming pili (bfp) in 72 strains identified as enteropathogenic E. coli (EPEC) by slide agglutination was evaluated using hybridization and PCR. The adherence property of these strains was assayed using 3h HeLa cells adherence assay. The results obtained indicated that virulence-associated genes were present in 65% of the strains but only ten (13.9%) isolates were positive for all the three markers (typical EPEC), 37 (51.4%) isolates carried either one or two of these determinants (atypical EPEC) and the remaining 25 (34.7%) were negative for all these genes. In vitro adherence assay showed that 44 (61.1%) strains adhered to HeLa cells with a defined pattern, 13 (18.1%) isolates adhered loosely with no definite pattern and the remaining 15 (20.8%) were non-adherent. Analysis of the results showed a statistically significant association between the presence of the virulence-related genes with adherence of the strains with a defined pattern (P相似文献