Evolution of the primate beta-globin gene region. High rate of variation in CpG dinucleotides and in short repeated sequences between man and chimpanzee

A 5500 base-pair fragment including the beta-globin gene downstream from codon 122 and about 4000 base-pairs of its 5' flanking sequence was cloned from chimpanzee DNA and thoroughly sequenced before being compared with the corresponding human sequence: 88 point differences (83 substitutions and 5 deletions or insertions of 1 base-pair) were detected as well as seven more important deletion/insertion events. These changes occur preferentially in two kinds of structure. First, 40% of the CpG dinucleotides present in either human or chimpanzee sequences are affected by nucleotide variations. This corresponds to a divergence level considerably higher than that expected. Second, most short repeated sequences found in the 5' extragenic sequence are involved in mutational events (amplification or contraction of the number of basic motifs as well as point substitutions or deletions/insertions of 1 base-pair). Considering the very low level of nucleotide sequence divergence between these two closely related species, our data provide direct evidence for CpG and tandem array instability.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

Bent DNA is needed for recombinational enhancer activity in the site-specific recombination system Cin of bacteriophage P1. The role of FIS protein

A series of recombinational enhancer mutants was constructed by manipulating the ClaI site between the two FIS binding sites of the Hin enhancer. These mutants include insertions from two to 12 base-pairs and two deletions of one or two base-pairs. Recombinational enhancer activity was found only with four mutants carrying either a four base-pair substitution, ten base-pair insertions or a one base-pair deletion, respectively; two other ten base-pair insertion mutants, however, were inactive, although FIS protein binding was unaffected. So, besides binding of FIS protein to its specific sites within the enhancer sequence and the correct helical positioning of these sites on the DNA, another criterion for enhancer activity must be fulfilled. DNA bending assays identify this requirement as a change of the enhancer DNA conformation, which FIS protein is able to induce and to stabilize. This conformational change of the DNA can be blocked by mutations in the central segment between the two FIS binding sites of the Hin enhancer. This sequence has special functions for the recombinational enhancer activity.     

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.     

Contribution of opening and bending dynamics to specific recognition of DNA damage

Guanine-uracil (G.U) wobble base-pairs are a detrimental lesion in DNA. Previous investigations have shown that such wobble base-pairs are more prone to base-opening than the normal G.C base-pairs. To investigate the sequence-dependence of base-pair opening we have performed 5ns molecular dynamics simulations on G.U wobble base-pairs in two different sequence contexts, TGT/AUA and CGC/GUG. Furthermore, we have investigated the effect of replacing the guanine base in each sequence with a fluorescent guanine analogue, 6-methylisoxanthopterin (6MI). Our results indicate that each sequence opens spontaneously towards the major groove in the course of the simulations. The TGT/AUA sequence has a greater proportion of structures in the open state than the CGC/GUG sequence. Incorporation of 6MI yields wobble base-pairs that open more readily than their guanine counterparts. In order of increasing open population, the sequences are ordered as CGC相似文献   

Effect of neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.     

Solution structure of a DNA double helix incorporating four consecutive non-Watson-Crick base-pairs

A series of DNA 21-mers containing a variety of the 4 x 4 internal loop sequence 5'-CAAG-3'/3'-ACGT-5' were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers exhibit excellent resolution in the NMR spectra and reveal many unusual NOEs (nuclear Overhauser effect) that allow for the detailed characterization of a DNA hairpin incorporating a track of four different non-Watson-Crick base-pairs in the stem. These include a wobble C.A base-pair, a sheared A.C base-pair, a sheared A.G base-pair, and a wobble G.T base-pair. Significantly different twisting angles were observed between the base-pairs in internal loop that results with excellent intra-strand and inter-strand base stacking within the four consecutive mismatches and the surrounding canonical base-pairs. This explains why it melts at 52 degrees C even though five out of ten base-pairs in the stem adopt non-Watson-Crick pairs. However, the 4 x 4 internal loop still fits into a B-DNA double helix very well without significant change in the backbone torsion angles; only zeta torsion angles between the tandem sheared base-pairs are changed to a great extent from the gauche(-) domain to the trans domain to accommodate the cross-strand base stacking in the internal loop. The observation that several consecutive non-canonical base-pairs can stably co-exist with Watson-Crick base-pairs greatly increases the limited repertoire of irregular DNA folds and reveals the possibility for unusual structural formation in the functionally important genomic regions that have potential to become single-stranded.     

Nucleotide sequence of terminal repeats of 412 transposable elements of Drosophila melanogaster: A similarity to proviral long terminal repeats and its implications for the mechanism of transposition

We have sequenced the long terminal direct repeats (and adjacent DNA) of two members of the 412 family of transposable elements of Drosophila melanogaster cloned on fragments of DNA from strain Oregon R. The repeats of the first element are identical and 481 base-pairs long; the repeats of the second are also identical but are 571 base-pairs long. The first 482 base-pairs of the 571 base-pair sequence correspond to the 481 base-pair repeat differing by five base substitutions and one addition/deletion. The 571 base-pair repeats are rare. Each of these 412 elements is flanked by a four base-pair direct repeat, suggesting that insertion of a 412 element is associated with duplication of four base-pairs. Analysis of the “empty site” from strain Canton S corresponding to one of these elements supports this conclusion. The sequence of 481 base-pair repeats and of 412 DNA immediately adjacent to them show striking similarities to corresponding regions of vertebrate proviruses and we discuss the implications this may have for the mechanism of transposition.     

Spectrum of spontaneous frameshift mutations. Sequences of bacteriophage T4 rII gene frameshifts

The DNA sequences of 185 independent spontaneous frameshift mutations in the rIIB gene of bacteriophage T4 are described. Approximately half of the frameshifts, including those at hot spot sites, are fully consistent with classical proposals that frameshift mutations are produced by a mechanism involving the misaligned pairing of repeated DNA sequences. However, the remaining frameshifts are inconsistent with this model. Correlations between the positions of two base-pair frameshifts and the bases of DNA hairpins suggest that local DNA topology might influence frameshift mutation. Warm spots for larger deletions share the property of having endpoints adjacent to DNA sequences whose complementarity to sequences a few base-pairs away suggest that non-classical DNA misalignments may participate in deletion mutation. A model for duplication mutation as a consequence of strand displacement synthesis is discussed. In all, 15 frameshifts were complex combinations of frameshifts and base substitutions. Three of these were identical, and have extended homology to a sequence 256 base-pairs away that is likely to participate in the mutational event; the remainder are unique combinations of frameshifts and transversions. The frequency and diversity of complex mutants suggest a challenge to the assumption that the molecular evolution of DNA must depend primarily upon the accumulation of single nucleotide changes.     

A-form conformational motifs in ligand-bound DNA structures

Recognition and biochemical processing of DNA requires that proteins and other ligands are able to distinguish their DNA binding sites from other parts of the molecule. In addition to the direct recognition elements embedded in the linear sequence of bases (i.e. hydrogen bonding sites), these molecular agents seemingly sense and/or induce an "indirect" conformational response in the DNA base-pairs that facilitates close intermolecular fitting. As part of an effort to decipher this sequence-dependent structural code, we have analyzed the extent of B-->A conformational conversion at individual base-pair steps in protein and drug-bound DNA crystal complexes. We take advantage of a novel structural parameter, the position of the phosphorus atom in the dimer reference frame, as well as other documented measures of local helical structure, e.g. torsion angles, base-pair step parameters. Our analysis pinpoints ligand-induced conformational changes that are difficult to detect from the global perspective used in other studies of DNA structure. The collective data provide new structural details on the conformational pathway connecting A and B-form DNA and illustrate how both proteins and drugs take advantage of the intrinsic conformational mechanics of the double helix. Significantly, the base-pair steps which exhibit pure A-DNA conformations in the crystal complexes follow the scale of A-forming tendencies exhibited by synthetic oligonucleotides in solution and the known polymorphism of synthetic DNA fibers. Moreover, most crystallographic examples of complete B-to-A deformations occur in complexes of DNA with enzymes that perform cutting or sealing operations at the (O3'-P) phosphodiester linkage. The B-->A transformation selectively exposes sugar-phosphate atoms, such as the 3'-oxygen atom, ordinarily buried within the chain backbone for enzymatic attack. The forced remodeling of DNA to the A-form also provides a mechanism for smoothly bending the double helix, for controlling the widths of the major and minor grooves, and for accessing the minor groove edges of individual base-pairs.     

DNA base sequence changes induced by ultraviolet light mutagenesis of a gene on a chromosome in Chinese hamster ovary cells

The DNA base sequence changes induced by mutagenesis with ultraviolet light have been determined in a gene on a chromosome of cultured Chinese hamster ovary (CHO) cells. The gene was the Escherichia coli gpt gene, of which a single copy was stably incorporated and expressed in the CHO cell genome. The cells were irradiated with ultraviolet light and gpt- colonies were selected by resistance to 6-thioguanine. The gpt gene was amplified from chromosomal DNA by use of the polymerase chain reaction (PCR), and the amplified DNA sequenced directly by the dideoxy method. Of the 58 sequenced mutants of independent origin 53 were base change mutations. Forty-one base substitutions were single base changes, ten had two adjacent (or tandem) base changes, and one had two base changes separated by a single base-pair. Only one mutant had a multiple base change mutation with two or more well separated base changes. In contrast much higher levels of such mutations were reported in ultraviolet mutagenesis of genes on a shuttle vector in primate cells. Two deletions of a single base-pair were observed and three deletions ranging from 6 to 37 base-pairs. The mutation spectrum in the gpt gene had similarities to the ultraviolet mutation spectra for several genes in prokaryotes, which suggests similarities in mutational mechanisms in prokaryotes and eukaryotes.     

Base sequence and helix structure variation in B and A DNA

The observed propeller twist in base-pairs of crystalline double-helical DNA oligomers improves the stacking overlap along each individual helix strand. But, as proposed by Calladine, it also leads to clash or steric hindrance between purines at adjacent base-pairs on opposite strands of the helix. This clash can be relieved by: (1) decreasing the local helix twist angle between base-pairs; (2) opening up the roll angle between base-pairs on the side on which the clash occurs; (3) separating purines by sliding base-pairs along their long axes so that the purines are partially pulled out of the stack (leading to equal but opposite alterations in main-chain torsion angle delta at the two ends of the base-pair); and (4) flattening the propeller twist of the offending base-pairs. Simple sum functions, sigma 1 through sigma 4, are defined, by which the expected local variation in helix twist, base roll angle, torsion angle delta and propeller twist may be calculated from base sequence. All four functions are quite successful in predicting the behavior of B DNA. Only the helix twist and base roll functions are applicable to A DNA, and the helix twist function begins to fail for an A helical RNA/DNA hybrid. Within these limits, the sequence-derived sum functions match the observed helix parameter variation quite closely, with correlation coefficients greater than 0.900 in nearly all cases. Implications of this sequence-derived helix parameter variation for repressor-operator interactions are considered.     

Analysis of local helix geometry in three B-DNA decamers and eight dodecamers

Local variations in B-DNA helix structure are compared among three decamers and eight dodecamers, which contain examples of all ten base-pair step types. All pairwise combinations of helix parameters are compared by linear regression analysis, in a search for internal relationships as well as correlations with base sequence. The primary conclusions are: (1) Three-center hydrogen bonds between base-pairs occur frequently in the major groove at C-C, C-A, A-A and A-C steps, but are less convincing at C-C and C-T steps in the minor groove. The requirements for large base-pair propeller are (1) that the base-pair should be A.T rather than G.C, and (2) that it be involved in a major groove three-center hydrogen bond with the following base-pair. Either condition alone is insufficient. Hence, a large propeller is expected at the leading base-pair of A-A and A-C steps, but not at A-T, T-A, C-A or C-C steps. (2) A systematic and quantitative linkage exists between helix variables twist, rise, cup and roll, of such strength that the rise between base-pairs can hardly be described as an independent variable at all. Two typical patterns of behavior are observed at steps from one base-pair to the next: high twist profile (HTP), characterized by high twist, low rise, positive cup and negative roll, and low twist profile (LTP), marked by low twist, high rise; negative cup and positive roll. Examples of HTP are steps G-C, G-A and Y-C-A-R, where Y is pyrimidine and R is purine. Examples of LTP steps are C-G, G-G, A-G and C-A steps other than Y-C-A-R. (3) The minor groove is especially narrow across the two base-pairs of the following steps: A-T, T-A, A-A and G-A. (4) In general, base step geometry cannot be correlated solely with the bases that define the step in question; the two flanking steps also must be taken into account. Hence, local helix structure must be studied in the context, not of two base-pairs: A-B, but of four: x-A-B-y. Calladine's rules, although too simple in detail, were correct in defining the length of sequence over which a given perturbation is expressed. Whereas ten different two-base steps are possible, allowing for the identity of complementary sequences, there are 136 different four-base steps. Only 33 of these 136 four-base steps are represented in the decamer and dodecamer structures solved to date, and hence it is premature to try to set up detailed structural algorithms. (5) The sugar-phosphate backbone chains of B-DNA place strong limits on sequence-induced structural variation, damping down most variables within four or five base-pairs, and preventing purine-purine anti-anti mismatches from causing bulges in the double helix. Hence, although short-range sequence-induced deformations (or deformability) are observed, long-range deformations propagated down the helix are not to be expected.     

Targeting DNA polymerase to DNA double-strand breaks reduces DNA deletion size and increases templated insertions generated by CRISPR/Cas9

Most insertions or deletions generated by CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9) endonucleases are short (<25 bp), but unpredictable on-target long DNA deletions (>500 bp) can be observed. The possibility of generating long on-target DNA deletions poses safety risks to somatic genome editing and makes the outcomes of genome editing less predictable. Methods for generating refined mutations are desirable but currently unavailable. Here, we show that fusing Escherichia coli DNA polymerase I or the Klenow fragment to Cas9 greatly increases the frequencies of 1-bp deletions and decreases >1-bp deletions or insertions. Importantly, doing so also greatly decreases the generation of long deletions, including those >2 kb. In addition, templated insertions (the insertion of the nucleotide 4 nt upstream of the protospacer adjacent motif) were increased relative to other insertions. Counteracting DNA resection was one of the mechanisms perturbing deletion sizes. Targeting DNA polymerase to double-strand breaks did not increase off-targets or base substitution rates around the cleavage sites, yet increased editing efficiency in primary cells. Our strategy makes it possible to generate refined DNA mutations for improved safety without sacrificing efficiency of genome editing.     

rII cistrons of bacteriophage T4. DNA sequence around the intercistronic divide and positions of genetic landmarks

An 873 base-pair DNA sequence from the rII region of bacteriophage T4 is presented. The sequence encodes 139 carboxyl-terminal amino acids of rIIA and the amino-terminal 146 amino acids of rIIB. Eleven base-pairs separate the rIIA stop codon (UAA) and the rIIB AUG.An extensive genetic map is superimposed on the DNA sequence, showing the deduced locations of many of the mutations (base-pair substitutions, frameshifts, deletions) found in previous rII genetic studies.     

Sequence specificity of 125I-labelled Hoechst 33258 damage in six closely related DNA sequences

The sequence selectivity of [125I]Hoechst 33258 in six 340 base-pair DNA sequences has been investigated. [125I]Hoechst 33258, which is a bis-benzimidazole and binds to the minor groove of B-DNA, preferentially binds to A + T-rich regions of DNA. Six out of nine strong binding sites contained four or more consecutive A.T base-pairs, while the other three strong binding sites were AAGGATT, TATAGAAA (the peak of damage was in the run of 3 A residues) and AAA. One of the six weak binding sites had five consecutive A.T base-pairs, two of the weak binding sites had three, and three did not have any. In addition to genomic 340 base-pair alpha RI-DNA (which is a tandem repeat in human cells), five 340 base-pair alpha RI-DNA clones were generated that differed from the genomic "consensus" sequence by a number of random base alterations. The effect of these base changes on the sequence specificity of [125I]Hoechst 33258 damage indicated that of the base changes that interrupted 14 binding sites, six decreased and eight did not change the extent of damage, while two sites changed position. Of the base alterations that augmented 17 binding sites, five increased, two decreased and ten did not alter the degree of cleavage, while ten sites changed position. It was concluded from the data that, while runs of consecutive A.T base-pairs was the most important parameter that determines [125I]Hoechst 33258 binding, other factors including position in the DNA sequence, nearest neighbour and long-range interactions were also important.     

Flanking A·T Basepairs Destabilize the B? Conformation of DNA A-Tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA⁺) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA⁺ and slower than the reference in Li⁺. Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA⁺ suggest that the large hydrophobic TBA⁺ ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li⁺ ions decrease the mobility of A-tract oligomers because of the preferential localization of Li⁺ ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B^∗ conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A·T basepairs in the flanking sequences.     

Mutagenic specificity of ultraviolet light

Genetic and sequencing studies of ultraviolet light (u.v.)-induced mutations in the lacI gene of Escherichia coli show the following: u.v. stimulates many types of mutations. In lacI, base substitutions account for 60 to 65% of the observed mutations, small frameshifts 30 to 35%, and deletions of more than several base-pairs approximately 5%. A comparison of the mutational spectrum of u.v.-induced mutations with those of other SOS-dependent mutagens and with the mutations produced by inducing the SOS system in the absence of mutagenic treatment indicates that most u.v.-induced base substitutions are "targeted", resulting from premutational lesions across from the site of the mutations. Among base substitutions, both transitions and transversions occur, although the most favored mutational sites involve G X C----A X T transitions. G X C----A X T transitions are induced preferentially at sites of adjacent pyrimidines. In one case the conversion of a site from -A-C-A- to -T-C-A- results in a 15-fold increase in u.v.-induced C----T transitions. Frameshifts at certain sites are well-induced by u.v., and the largest hotspot in the I gene involves the loss of an (sequence in text) base pair from a (sequence in text) sequence. Of 25 frameshifts detected by DNA sequencing, 23 mutations at seven different sites result from the elimination of a single base-pair, and two mutations result from the elimination of two base-pairs. No additions were detected. The use of a lacI-Z fusion system, which allows direct selection of frameshifts of either sign, reveals that throughout the entire gene frameshifts that eliminate a single base-pair (-1) predominate by a factor of 20 or more over frameshifts that add a single base-pair (+1). In one case a two-base-pair elimination occurs frequently, resulting in the loss of a -C-T- sequence (on one strand), or a -T-C- sequence, from a -C-T-C-T-C-T-C- sequence. For both frameshifts and base substitutions, some aspect of the larger surrounding sequence beyond the nearest neighbors can influence mutation rates by as much as 50-fold, thus determining which sites are seen as hotspots. The bearing of these and other data on the detailed mechanism of mutagenesis is considered in the Discussion.