Strong Patterns in Homooligomer Tracts Occurrences in Non-Coding and in Potential Regulatory Sites in Eukaryotic Genomes

Abstract

Previous studies of the dinucleotides flanking both the 5′ and 3′ ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g., ATA_n, TTA_n, CCG_n, where n=2–5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group A/T_n are flanked by C and/or G; G_n/C_n are flanked by A/T, e.g., CGA_n, T_nGG, G., AT). The frequencies of runs flanked by AorT, and G or C (“mixed” group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g., TTA_n, C_nGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.     

Distinct patterns in homooligomer tract sequence context in prokaryotic and eukaryotic DNA

The distributions of the junction sequences of homooligomer tracts of various lengths have been examined in prokaryotic DNA sequences and compared with those of eukaryotes. The general trends in the nearest and next to nearest neighbors to the tracts are similar for both groups. In both prokaryotes and eukaryotes A/T runs are preferentially flanked on either the 5' or the 3' ends by A and/or T. G/C runs are preferentially flanked by G and/or C. There is discrimination against A/T runs flanked by G or C and G/C runs flanked by A or T. However, whereas the distribution of prokaryotic homooligomer tract junction sequences was quite homogeneous, large variations were observed in the 5-fold larger eukaryotic database, increasing in magnitude from tracts of length 2 to 3 to 4 base pairs long. Possible DNA conformational implications and in particular DNA curvature and packaging aspects of prokaryotes and eukaryotes are discussed.     

The ordering of nucleotides in the DNA: strong pyrimidine-purine patterns near homooligomer tracts

Here, we study the frequencies of occurrence of homooligomers flanked by one base, XnU or UXn, where X = A, C, G, T and U not equal to X. Specifically, we search for preferences (or discriminations) in their nearest neighbor doublet, VV. Extensive analysis of the data base reveals striking patterns in such VVUXn or UXn VV oligomers (V = A, C, G, T). With very few exceptions, if the VV and Xn are composed of complementary nucleotides, those oligomers having a pyrimidine (Y)-purine (R) junction are preferred over those with an RY one. If the VV and Xn nucleotides are not complementary, the RY junction oligomers are preferred over their YR counterparts. These trends are observed consistently in eukaryotic and prokaryotic sequences. They are particularly striking in the YR greater than RY oligomers containing complementary nucleotides. The general preferences and discriminations described here are in the same direction as our previous results for homooligomer tracts. These recurrences, along with some additional universal "rules", aid in our understanding of the ordering of nucleotides in the DNA.     

Distinct patterns in the dinucleotide nearest neighbors to G/C and A/T oligomers in eukaryotic sequences

Summary The eukaryotic and prokaryotic databases are scanned for potential nearest-neighbor doublet preferences at the 5 and 3 flanks of some oligomers. Here we focus on oligomers containing alternating nucleotides, i.e., UV, UVUV, and UUVV where UV. Strong, consistent trends are observed in eukaryotic sequences. A/T alternation oligomers are preferentially flanked by A/T. G/C flanks are disfavored. G/C alternation oligomers are preferentially flanked by G/C. A/T flanks are disfavored. These trends are consistent with those observed previously for homooligomer tracts (Nussinov et al. 1989a,b). G/C tracts are preferentially flanked by G/C. A/T nearest neighbors are disfavored. The reverse holds for A/T tracts. Additional patterns are described here as well. The possible origin of these DNA composition and sequence trends is discussed. These trends are suggested to stem from protein-DNA interaction constraints.     

Sequence context of oligomer tracts in eukaryotic DNA: biological and conformational implications

Recent studies of homooligomer tracts suggest different characteristics from random sequence DNA (dA).(dT) and (dG).(dC) tracts are frequent in upstream regions and in some cases have been shown to be essential for regulation. Here we examine homooligomer occurrences in non-coding and coding eukaryotic sequences, focusing on the context in which the homooligomers occur. This analysis of sequences in the junction areas yields distinct and consistent characteristics. In particular, the nucleotide interrupting a run is most frequently complementary to the run. The base next to it is most frequently identical to the one constituting the run. For A or T runs the least frequent nearest and next to nearest neighbors are G or C. For G or C tracts the least frequent are A or T. Complementary oligomers behave similarly. These and additional trends are strongest for run lengths greater than or equal to 3. The computations are carried out on the whole eukaryotic database of greater than 4 x 10(6) nucleotides, separately for coding and non-coding regions. These same trends are evident for both groups, but are somewhat stronger for the non-coding regions. The context in which the homooligomers occur may yield some clues to DNA conformation and its biological implications.     

Conserved signals around the 5' splice sites in eukaryotic nuclear precursor mRNAs: G-runs are frequent in the introns and C in the exons near both 5' and 3' splice sites

It is known that the GT doublet is well conserved at the 5' exon/intron splice junction and is frequently embedded in the AGGT quartet. Although only the underlined G is invariable, splicing and ligation are accurately executed. In this work we search for additional conserved potential signals which may aid in 5' splice site recognition. Extensive searches which are not limited to a preconceived consensus sequence are carried out. We investigate the distributions of the 256 quartets in a 1000 nucleotide span around the 5' splice sites in approximately 1700 eukaryotic nuclear precursor mRNAs. Several potential signals are noted. Of particular interest are quartets containing runs of G, e.g., G4, G3T, G3C, G3A and AG3 in the intron immediately downstream and some C-containing quartets in the exon upstream of the 5' splice site. In an analogous calculation, (A)GGG(A) has also been found to be frequent in the intron, 60 nucleotides upstream and (A)CCC(A) in the exon downstream of the 3' splice site. These results are consistent with the recent indications that exon sequences may play a role in efficient splicing. Some models are proposed.     

Consecutive GA pairs stabilize medium-size RNA internal loops

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.     

A-tract polarity dominate the curvature in flanking sequences

It is well-known, but little understood, that the nucleotide sequences between phased A(4-6)-tracts (at 10-11 bp intervals) have only a slight effect on overall curvature. To explore this phenomenon, we have examined the gel-migration properties of sequences containing both A-tracts as well as G-tracts (i.e., sequences of the form G(n)C(m) or C(n)G(m), n + m > 4) in various relative positioning. We show that the composite bend of these sequences depends on their relative arrangement. When G-tracts are placed between two A-tracts, such that both tracts are repeated in phase to themselves (e.g., G(5)A(6)G(5)A(5)), or adjacent to the 3'-side of A-tracts (e.g., A(6)G(5)N(10)), they have minimal influence on the extent of bending of the composite sequence. When G-tracts are placed one helical repeat away from A-tracts (e.g., G(5)N(5)A(6)N(6)), or are adjacent only to the 5'-side of A-tracts (e.g., G(5)A(6)N(10)) their influence on the composite bend is larger. The differential behavior of AG- versus GA-tracts means that A-tracts influence their flanking sequences in a polar manner. Whereas they suppress, or make constant, the intrinsic bending characteristics of any sequence placed immediately 3' to them (and hence by definition any sequence placed between two phased A-tracts), sequences adjoining them on their 5'-side are free to modulate the overall curvature. We interpret these results as evidence for the dominant nature of the unique and nonuniform structure adopted by tracts of four adenines or more. The effects of A-tracts extend at least five base pairs into the adjoining 3' region. This is further evidence for the complexity of DNA structure and the inadequacy of simple nearest-neighbor models to explain all its manifestations.     

The cytoskeletal adapter protein 4.1G organizes the internodes in peripheral myelinated nerves

Myelinating Schwann cells regulate the localization of ion channels on the surface of the axons they ensheath. This function depends on adhesion complexes that are positioned at specific membrane domains along the myelin unit. Here we show that the precise localization of internodal proteins depends on the expression of the cytoskeletal adapter protein 4.1G in Schwann cells. Deletion of 4.1G in mice resulted in aberrant distribution of both glial adhesion molecules and axonal proteins that were present along the internodes. In wild-type nerves, juxtaparanodal proteins (i.e., Kv1 channels, Caspr2, and TAG-1) were concentrated throughout the internodes in a double strand that flanked paranodal junction components (i.e., Caspr, contactin, and NF155), and apposes the inner mesaxon of the myelin sheath. In contrast, in 4.1G(-/-) mice, these proteins "piled up" at the juxtaparanodal region or aggregated along the internodes. These findings suggest that protein 4.1G contributes to the organization of the internodal axolemma by targeting and/or maintaining glial transmembrane proteins along the axoglial interface.     

Photoaffinity approaches to determining the sequence selectivities of DNA-small molecule interactions: actinomycin D and ethidium.

The DNA photoaffinity ligands, 7-azidoactinomycin D and 8-azidoethidium, form DNA adducts that cause chain cleavage upon treatment with piperidine. Chemical DNA sequencing techniques were used to detect covalent binding. The relative preferences for modifications of all possible sites defined by a base pair step (e.g. GC) were determined within all quartet contexts such as (IGCJ). These preferences are described in terms of 'effective site occupations', which express the ability of a ligand to covalently modify some base in the binding site. Ideally, the effective site occupations measured for photoaffinity agents can also be related to site-specific, non-covalent association constants of the ligand. The sites most reactive with 7-azidoactinomycin D were those preferred for non-covalent binding of unsubstituted actinomycin D. GC sites were most reactive, but next-nearest neighbors exerted significant influences on reactivity. GC sites in 5'-(pyrimidine)GC(purine)-3' contexts, particularly TGCA, were most reactive, while reactivity was strongly suppressed for GC sites with a 5'-flanking G, or a 3'-flanking C. High reactivities were also observed for bases in the first (5') GG steps in TGGT, TGGG and TGGGT sequences recently shown to bind actinomycin D with high affinity. Pyrimidine-3',5'-purine steps and GG steps flanked by a T were most preferred by 8-azidoethidium, in agreement with the behavior of unsubstituted ethidium. The good correspondence between expected and observed covalent binding preferences of these two azide analogs demonstrates that photoaffinity labeling can identify highly preferred sites of non-covalent DNA binding by small molecules.     

The 5'' splice site: phylogenetic evolution and variable geometry of association with U1RNA.

The 5' splice site sequences of 3294 introns from various organisms (1-672) were analyzed in order to determine the rules governing evolution of this sequence, which may shed light on the mechanism of cleavage at the exon-intron junction. The data indicate that, currently, in all organisms, a common sequence 1GUAAG6U and its derivatives are used as well as an additional sequence and its derivatives, which differ in metazoa (G/1GUgAG6U), lower eucaryotes (1GUAxG6U) and higher plants (AG/1GU3A). They all partly resemble the prototype sequence AG/1GUAAG6U whose 8 contigous nucleotides are complementary to the nucleotides 4-11 of U1RNA, which are perfectly conserved in the course of phylogenetic evolution. Detailed examination of the data shows that U1RNA can recognize different parts of 5' splice sites. As a rule, either prototype nucleotides at position -2 and -1 or at positions 4, 5 or 6 or at positions 3-4 are dispensable provided that the stability of the U1RNA-5' splice site hybrid is conserved. On the basis of frequency of sequences, the optimal size of the hybridizable region is 5-7 nucleotides. Thus, the cleavage at the exon-intron junction seems to imply, first, that the 5' splice site is recognized by U1RNA according to a "variable geometry" program; second, that the precise cleavage site is determined by the conserved sequence of U1RNA since it occurs exactly opposite to the junction between nucleotides C9 and C10 of U1RNA. The variable geometry of the U1RNA-5' splice site association provides flexibility to the system and allows diversification in the course of phylogenetic evolution.     

Initiation of RNA-primed DNA synthesis in vitro by DNA polymerase alpha-primase.

The initiation of new DNA strands at origins of replication in animal cells requires de novo synthesis of RNA primers by primase and subsequent elongation from RNA primers by DNA polymerase alpha. To study the specificity of primer site selection by the DNA polymerase alpha-primase complex (pol alpha-primase), a natural DNA template containing a site for replication initiation was constructed. Two single-stranded DNA (ssDNA) molecules were hybridized to each other generating a duplex DNA molecule with an open helix replication 'bubble' to serve as an initiation zone. Pol alpha-primase recognizes the open helix region and initiates RNA-primed DNA synthesis at four specific sites that are rich in pyrimidine nucleotides. The priming site positioned nearest the ssDNA-dsDNA junction in the replication 'bubble' template is the preferred site for initiation. Using a 40 base oligonucleotide template containing the sequence of the preferred priming site, primase synthesizes RNA primers of 9 and 10 nt in length with the sequence 5'-(G)GAAGAAAGC-3'. These studies demonstrate that pol alpha-primase selects specific nucleotide sequences for RNA primer formation and suggest that the open helix structure of the replication 'bubble' directs pol alpha-primase to initiate RNA primer synthesis near the ssDNA-dsDNA junction.     

Regularities of the nucleotide sequence at the 5'-end of the codon in Escherichia coli genes

The frequencies of occurrence of nucleotides at the 5' side of codons have been determined in highly and weakly expressed genes from E. coli. Significant constraints on the nucleotide 5' to some codons were found in highly expressed genes. Certain rules of synonymous codon usage depending on the amino acid 3' of the codon were established. E. g., codon possessing quanosine in the third position (NNG) are preferred over NNA if the next amino acid is lysine (P less than 10(-5)). On the other hand, rules of synonymous codon usage in relation to 5' flanking nucleotide were found. For example, when coding for aspartic acid, GAC codon is preferred over GAU (P less than 0.001) if uridine is 5' to codon and on the contrary GAU is favoured (P less than 0.0001) if quanosine is at the 5' side of aspartic acid codon. These rules can be used in the chemical synthesis of genes designed for expression in E. coli.     

Kinetics of extension of O6-methylguanine paired with cytosine or thymine in defined oligonucleotide sequences.

The frequency of extending m6G.C or m6G.T pairs, when the 3' and 5' flanking neighbors of m6G are either cytosines or thymines, was investigated using primed 25-base-long oligonucleotides and the Klenow fragment of Escherichia coli DNA polymerase I (Kf). The efficiency, Vmax/Km, of extension to the following normal base pair was up to 40-fold greater than for the formation of the m6G.T or m6G.C pair. The frequencies of inserting either dCMP or dTMP opposite these m6G bases did not appear to be different in the two sequences, C-m6G-C and T-m6G-T, but extension was favored in the C-m6G-C sequence. The m6G.T pair extended to a C.G pair most efficiently, indicating that it was not a strong block to continued replication past the template lesion. Thus, m6G.T flanked by cytosines replicates more readily than when flanked by thymines, increasing G----A transitions. These data lend further support to the importance of sequence context in mutagenesis.     

Purification and RNA binding properties of a C-type hnRNP protein from HeLa cells

A protein of the C group, most likely C3 (Mr approximately 42,000, pI approximately 6, corresponding to IEF 48m,n of the HeLa protein catalogue (Celis, J. E., Bravo, R., Arenstorf, H. P., and LeStourgeon, W. M. (1986) FEBS Lett. 194, 101-109)), a minor hnRNP protein was purified to near homogeneity under nondenaturing conditions from 40 S heterogeneous nuclear ribonucleoprotein particles. Type C protein stoichiometrically disrupts the residual secondary structure of natural and synthetic RNAs, e.g. HeLa hnRNA, coliphage MS2 RNA, and poly(rU)-spermine, and decreases the Tm of duplex structures, e.g. poly[r(A + U)], by about 30 degrees C. Binding of the protein to polynucleotides is not highly cooperative and has a stoichiometry of one protein per about 10 nucleotides. Binding experiments with a variety of synthetic and natural poly- and oligonucleotides, including those containing consensus splice site sequences, indicate that the protein has a high affinity for G-rich and U-rich regions, G-rich regions being preferred. Base analogs I and T have affinities for the protein that are similar to G and U. There is little or no affinity for A- and C-rich regions. The presence of A residues in a G- or U-rich sequence does not interfere with binding while C-rich regions decrease or prevent the binding of the protein. The nucleotide specificity of type C protein, e.g. selective binding to an oligonucleotide from the 3' end of an intron, is discussed in relationship to the abundance of G and U and the relative scarcity of C residues in the processing signals in pre-mRNA.     

Duplicated rDNA sequences of variable lengths flanking the short type I insertions in the rDNA of Drosophila melanogaster.

We describe cloned segments of rDNA that contain short type I insertions of differing lengths. These insertions represent a coterminal subset of sequences from the right hand side of the major 5kb type I insertion. Three of these shorter insertions are flanked on both sides by a short sequence present as a single copy in uninterrupted rDNA units. The duplicated segment is 7, 14 and 15 nucleotides in the different clones. In this respect, the insertions differ from the 5kb type I insertion, where the corresponding sequence is found only at the right hand junction and where at the left hand side there is a deletion of 9 nucleotides of rDNA (Roiha et al.,1981). One clone is unusual in that it contains two type I insertions, one of which is flanked by a 14 nucleotide repeat. The left hand junction of the second insertion occurs 380 nucleotides downstream in the rDNA unit from the first. It has an identical right hand junction to the other elements and the 380 nucleotide rDNA sequence is repeated on both sides of the insertion. We discuss the variety of sequence rearrangements of the rDNA which flank type I insertions.