Comparison of the consensus sequence flanking translational start sites in Drosophila and vertebrates.

The previously presented consensus sequence for eukaryotic translation initiation sites by Kozak was derived substantially from vertebrate mRNA sequences. Drosophila nuclear genes exhibit a significantly different translation start consensus sequence. These differences probably do not represent mechanistic differences in translation initiation inasmuch as both taxa exhibit identical preferences and restrictions at the crucial -3 position. Using more conservative criteria for the assignment of consensus the following consensus sequences were derived: vertebrate--CANCAUG and Drosophila--CAAAACAUG.     

Accumulating variation at conserved sites in potyvirus genomes is driven by species discovery and affects degenerate primer design

Unknown and foreign viruses can be detected using degenerate primers targeted at conserved sites in the known viral gene sequences. Conserved sites are found by comparing sequences and so the usefulness of a set of primers depends crucially on how well the known sequences represent the target group including unknown sequences. METHODOLOGY/PRINCIPAL FINDINGS: We developed a method for assessing the apparent stability of consensus sequences at sites over time using deposition dates from Genbank. We tested the method using 17 conserved sites in potyvirus genomes. The accumulation of knowledge of sequence variants over 20 years caused 'consensus decay' of the sites. Rates of decay were rapid at all sites but varied widely and as a result, the ranking of the most conserved sites changed. The discovery and reporting of sequences from previously unknown and distinct species, rather than from strains of known species, dominated the decay, indicating it was largely a sampling effect related to the progressive discovery of species, and recent virus mutation was probably only a minor contributing factor. CONCLUSION/SIGNIFICANCE: We showed that in the past, the sampling bias has misled the choice of the most conserved target sites for genus specific degenerate primers. The history of sequence discoveries indicates primer designs should be updated regularly and provides an additional dimension for improving the design of degenerate primers.     

A symmetrical six-base-pair target site sequence determines Tn10 insertion specificity

Transposon Tn10 inserts at many sites in the bacterial chromosome, but preferentially inserts at particular hotspots. We believe we have identified the target DNA signal responsible for this specificity. We have determined the DNA sequences of 11 Tn10 insertion sites and identified a particular 6 base pair (bp) symmetrical consensus sequence (GCTNAGC) common to those sites. The sequences at some sites differ from the consensus sequence but only in limited and well defined ways. The sequences at some sites differ from the consensus sequence than do sequences at other sites, and the consensus sequence and closely related sequences are generally absent from potential target regions where Tn10 is known not to insert. Other aspects of the target DNA can significantly influence the efficiency with which a particular target site sequence is used. The 6 bp consensus sequence is symmetrically located within the 9 bp target DNA sequence that is cleaved and duplicated during Tn10 insertion. This juxtaposition of recognition and cleavage sites plus the symmetry of the perfect consensus sequence suggest that the target DNA may be both recognized and cleaved by the symmetrically disposed subunits of a single protein, as suggested for type II restriction endonucleases. There is plausible homology between the consensus sequence and the very ends of Tn10, compatible with recognition of transposon ends and target DNA by the same protein. The sequences of actual insertion sites deviate from the perfect consensus sequence in a way which suggests that the 6 bp specificity determinant may be recognized through protein-DNA contacts along the major groove of the DNA double helix.     

Target specificity of insertion element IS30

The Escherichia coli resident mobile element IS 30 has pronounced target specificity. Upon transposition, the element frequently inserts exactly into the same position of a preferred target sequence. Insertion sites in phages, plasmids and in the genome of E. coli are characterized by an exceptionally long palindromic consensus sequence that provides strong specificity for IS 30 insertions, despite a relatively high level of degeneracy. This 24-bp-long region alone determines the attractiveness of the target DNA and the exact position of IS 30 insertion. The divergence of a target site from the consensus and the occurrence of 'non-permitted' bases in certain positions influence the target activity. Differences in attractiveness are emphasized if two targets are present in the same replicon, as was demonstrated by quantitative analysis. In a system of competitive targets, the oligonucleotide sequence representing the consensus of genomic IS 30 insertion sites proved to be the most efficient target. Having compared the known insertion sites, we suppose that IS 30 -like target specificity, which may represent an alternative strategy in target selection among mobile elements, is characteristic of the insertion sequences IS 3 , IS 6 and IS 21 , too.     

Analysis of the ends of bacteriophage Mu using site-directed mutagenesis

We showed previously that two regions at the left end (L1 and L3) and one at the right end (R2) of bacteriophage Mu are essential for transposition. These regions all contain a 22 base-pair sequence with the consensus YGtTTCAYtNNAARYRCGAAAR, where Y and R represent any pyrimidine and purine, respectively. The Mu A protein binds to these regions in vitro, and weakly to sequences between nucleotides 1 and 30 of the right end (R1) and between nucleotides 110 and 135 of the left end (L2). These weak A binding sites contain the sequence AARYRCGAAAR. Here we show, using site-directed mutagenesis, that the weak A binding sites are essential for transposition. Mutations in these weak A binding sites have a greater effect on transposition than mutations of corresponding base-pairs in the stronger A binding sites, located adjacent to these weak A binding sites. We confirm the importance of several of the conserved base-pairs in the consensus sequence YGtTTCAYtNNAARYRCGAAAR. The base-pairs in the A binding sites that are shown to be essential for transposition are all conserved in the ends of the related bacteriophage D108. Furthermore, it is shown that the distance of 90 base-pairs between the two regions at the left end (L1 and L2L3) is essential.     

Combinatorial determination of sequence specificity for nanomolar DNA-binding hairpin polyamides

Development of sequence-specific DNA-binding drugs is an important pharmacological goal, given the fact that numerous existing DNA-directed chemotherapeutic drugs rely on the strength and selectivity of their DNA interactions for therapeutic activity. Among the DNA-binding antibiotics, hairpin polyamides represent the only class of small molecules that can practically bind any predetermined DNA sequence. DNA recognition by these ligands depends on their side-by-side amino acid pairings in the DNA minor groove. Extensive studies have revealed that these molecules show extremely high affinity for sequence-directed, minor groove interaction. However, the specificity of such interactions in the presence of a large selection of sequences such as the human genome is not known. We used the combinatorial selection method restriction endonuclease protection, selection, and amplification (REPSA) to determine the DNA binding specificity of two hairpin polyamides, ImPyPyPy-gamma-PyPyPyPy-beta-Dp and ImPyPyPy-gamma-ImPyPyPy-beta-Dp, in the presence of more than 134 million different sequences. These were verified by restriction endonuclease protection assays and DNase I footprinting analysis. Our data showed that both hairpin polyamides preferentially selected DNA sequences having consensus recognition sites as defined by the Dervan pairing rules. These consensus sequences were rather degenerate, as expected, given that the stacked pyrrole-pyrrole amino acid pairs present in both polyamides are unable to discriminate between A.T and T.A base pairs. However, no individual sequence within these degenerate consensus sequences was preferentially selected by REPSA, indicating that these hairpin polyamides are truly consensus-specific DNA-binding ligands. We also discovered a preference for overlapping consensus binding sites among the sequences selected by the hairpin polyamide ImPyPyPy-gamma-PyPyPyPy-beta-Dp, and confirmed by DNase I footprinting that these complex sites provide higher binding affinity. These data suggest that multiple hairpin polyamides can cooperatively bind to their highest-affinity sites.     

A simulated annealing algorithm for finding consensus sequences

MOTIVATION: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. RESULTS: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.     

Insertion site specificity of the transposon Tn3.

The Tn3-deletion method [Davies and Hutchison, Nucleic Acids Res. 19, 5731-5738, (1991)] was used to sequence a 9.4 kb DNA fragment. Transpositional 'warm' spots were not a limiting factor but a 935 bp 'cold' spot was completed using a synthetic oligonucleotide primer. Two hundred and twenty three miniTn3 insertion sites from three sequencing projects were aligned and a 19 bp asymmetric consensus site was identified. There is no absolute sequence requirement at any position in this consensus, so insertion occurs promiscuously (approximately 37% of sites are potential targets). In our sequencing projects, multiply targeted sites always closely matched the consensus, although not all close matches were targeted frequently. The 935 bp cold spot showed no unusual features when analysed with the consensus sequence. The consensus can be used to accurately predict likely insertion sites in a new sequence. Synthetic oligonucleotides based on the consensus and a known hot spot for Tn3 were mutagenised. These sequences were not hot spots in our vectors, suggesting that the primary sequence alone is not sufficient to create an insertional hot spot. We conclude that some other factor, such as DNA secondary structure, also plays an important role in target site selection for the transposon Tn3.     

Base sequence studies of 300 nucleotide renatured repeated human DNA clones

A band of 300 nucleotide long duplex DNA is released by treating renatured repeated human DNA with the single strand-specific endonuclease S₁. Since many of the interspersed repeated sequences in human DNA are 300 nucleotides long, this band should be enriched in such repeats. We have determined the nucleotide sequences of 15 clones constructed from these 300 nucleotide S₁-resistant repeats. Ten of these cloned sequences are members of the Alu family of interspersed repeats. These ten sequences share a recognizable consensus sequence from which individual clones have an average divergence of 12.8%. The 300 nucleotide Alu family consensus sequence has a dimeric structure and was evidently formed from a head to tail duplication of an ancestral monomeric sequence. Three of the remaining clones are variations on a simple pentanucleotide sequence previously reported for human satellite III DNA. Two of the 15 clones have distinct and complex sequences and may represent other families of interspersed repeated sequences.     

Structural requirements for additional N-linked carbohydrate on recombinant human erythropoietin

N-Linked glycosylation is a post-translational event whereby carbohydrates are added to secreted proteins at the consensus sequence Asn-Xaa-Ser/Thr, where Xaa is any amino acid except proline. Some consensus sequences in secreted proteins are not glycosylated, indicating that consensus sequences are necessary but not sufficient for glycosylation. In order to understand the structural rules for N-linked glycosylation, we introduced N-linked consensus sequences by site-directed mutagenesis into the polypeptide chain of the recombinant human erythropoietin molecule. Some regions of the polypeptide chain supported N-linked glycosylation more effectively than others. N-Linked glycosylation was inhibited by an adjacent proline suggesting that sequence context of a consensus sequence could affect glycosylation. One N-linked consensus sequence (Asn123-Thr125) introduced into a position close to the existing O-glycosylation site (Ser126) had an additional O-linked carbohydrate chain and not an additional N-linked carbohydrate chain suggesting that structural requirements in this region favored O-glycosylation over N-glycosylation. The presence of a consensus sequence on the protein surface of the folded molecule did not appear to be a prerequisite for oligosaccharide addition. However, it was noted that recombinant human erythropoietin analogs that were hyperglycosylated at sites that were normally buried had altered protein structures. This suggests that carbohydrate addition precedes polypeptide folding.     

Isolation and characterization of a vinculin cDNA from chick-embryo fibroblasts.

A chick-embryo fibroblast lambda gt11 cDNA library was screened with affinity-purified antibodies to chick gizzard vinculin. One recombinant was purified to homogeneity and the fusion protein expressed in Escherichia coli strain C600. The fusion protein was unstable, but polypeptides that reacted with vinculin antibodies, but not non-immune immunoglobulin, were detected by Western blotting. The recombinant contained a single EcoRI fragment of 2891 bp with a single open reading frame. The deduced protein sequence could be aligned with that of six CNBr-cleavage peptides and two tryptic peptides derived from chicken gizzard vinculin. AUG-247 has tentatively been identified as the initiation codon, as it is contained within the consensus sequence for initiation sites of higher eukaryotes. The cDNA lacks 3' sequence and encodes 74% of the vinculin sequence, presuming the molecular mass of vinculin to be 130,000 Da. Analysis of the deduced sequence showed no homologies with other protein sequences, but it does display a triple internal repeat of 112 amino acid residues covering residues 259-589. The sequences surrounding the seven tyrosine residues in the available sequence were aligned with the tyrosine autophosphorylation consensus sequence found in protein tyrosine kinases. Tyr-822 showed a good match to this consensus, and may represent one of the two major sites of tyrosine phosphorylation by pp60v-sre. Northern blots showed that the 2.89 kb vinculin cDNA hybridized to one size of mRNA (approx. 7 kb) in chick-embryo fibroblasts, chick smooth muscle and chick skeletal muscle. Southern blots revealed multiple hybridizing bands in genomic DNA.     

Mitochondrial diversity within modern human populations

With the recent increase in the available number of high-quality, full-length mitochondrial sequences, it is now possible to construct and analyze a comprehensive human mitochondrial consensus sequence. Using a data set of 827 carefully selected sequences, it is shown that modern humans contain extremely low levels of divergence from the mitochondrial consensus sequence, differing by a mere 21.6nt sites on average. Fully 84.1% of the mitochondrial genome was found to be invariant and ‘private’ mutations accounted for 43.8% of the variable sites. Ninety eight percent of the variant sites had a primary nucleotide with an allele frequency of 0.90 or greater. Interestingly, the few truly ambiguous nucleotide sites could all be reliably assigned to either a purine or pyrimidine ancestral state. A comparison of this consensus sequence to several ancestral sequences derived from phylogenetic studies reveals a great deal of similarity, where, as expected, the most phylogenetically informative nucleotides in the ancestral studies tended to be the most variable nucleotides in the consensus. Allowing for this fact, the consensus approach provides variation data on the positions that do not contribute to phylogenetic reconstructions, and these data provide a baseline for measuring human mitochondrial variation in populations worldwide.     

Splice site consensus sequences are preferentially accessible to nucleases in isolated adenovirus RNA.

The conformation of RNA sequences spanning five 3' splice sites and two 5' splice sites in adenovirus mRNA was probed by partial digestion with single-strand specific nucleases. Although cleavage of nucleotides near both 3' and 5' splice sites was observed, most striking was the preferential digestion of sequences near the 3' splice site. At each 3' splice site a region of very strong cleavage is observed at low concentrations of enzyme near the splice site consensus sequence or the upstream branch point consensus sequence. Additional sites of moderately strong cutting near the branch point consensus sequence were observed in those sequences where the splice site was the preferred target. Since recognition of the 3' splice site and branch site appear to be early events in mRNA splicing these observations may indicate that the local conformation of the splice site sequences may play a direct or indirect role in enhancing the accessibility of sequences important for splicing.     

The cleavage specificity of RNase III.

We determined sites in lambda cII mRNA that are cleaved by RNase III in the presence of lambda OOP antisense RNA, using a series of OOP RNAs with different internal deletions. In OOP RNA-cII mRNA structures containing a potential region of continuous double-stranded RNA bounded by a non-complementary unpaired region, RNase III cleaved the cII mRNA at one or more preferred sites located 10 to 14 bases from the 3'-end of the region of continuous complementarity. Cleavage patterns were almost identical when the presumptive structure was the same continuously double-stranded region followed by a single-stranded bulge and a second short region of base pairing. The sequences of the new cleavage sites show generally good agreement with a consensus sequence derived from thirty-five previously determined cleavage sequences. In contrast, four 'non-sites' at which cleavage is never observed show poor agreement with this consensus sequence. We conclude that RNase III specificity is determined both by the distance from the end of continuous pairing and by nucleotide sequence features within the region of pairing.     

The presence of interstitial telomeric sequences in constitutional chromosome abnormalities.

We describe a novel chromosome structure in which telomeric sequences are present interstitially, at the apparent breakpoint junctions of structurally abnormal chromosomes. In the linear chromosomes with interstitial telomeric sequences, there were three sites of hybridization of the telomere consensus sequence within each derived chromosome: one at each terminus and one at the breakpoint junction. Telomeric sequences also were observed within a ring chromosome. The rearrangements examined were constitutional chromosome abnormalities with a breakpoint assigned to a terminal band. In each case (with the exception of the ring chromosome), an acentric segment of one chromosome was joined to the terminus of an apparently intact recipient chromosome. One case exhibited apparent instability of the chromosome rearrangement, resulting in somatic mosaicism. The rearrangements described here differ from the telomeric associations observed in certain tumors, which appear to represent end-to-end fusion of two or more intact chromosomes. The observed interstitial telomeric sequences appear to represent nonfunctional chromosomal elements, analogous to the inactivated centromeres observed in dicentric chromosomes.     

High affinity YY1 binding motifs: identification of two core types (ACAT and CCAT) and distribution of potential binding sites within the human beta globin cluster.

PCR-assisted binding site selection was used to define the sequence characteristics of high affinity YY1 binding sites. Compilation of the sequences of 189 selected oligonucleotides containing high affinity YY1 binding sites revealed two types of core sequence: ACAT and CCAT. ACAT cores were surrounded by other invariant nucleotides, forming the consensus GACATNTT. A search of the 73 kb human beta-like globin cluster with this consensus revealed eight matching motifs, six of which were located within 1-3 kb upstream of the gamma and beta genes. CCAT-type cores were more variable in surrounding sequence context; the consensus VDCCATNWY was found to fit 89% of the selected CCAT-containing oligonucleotides. A search of the human beta globin cluster with CCAT consensus sequences revealed 171 potential YY1 binding sites. Several of these were tested directly in gel shift assays and confirmed as high affinity YY1 binding sites. Finally, a strategy called motif-based phylogenetic analysis was employed to determine which of the 179 total sites are evolutionarily conserved. This analysis permits the detection of functionally conserved binding sites despite sequence differences present between the two species. The 21 conserved sites identified will serve as important starting points in further dissection of the possible role of YY1 in globin gene regulation.