Population genetics of dinucleotide (dC-dA)n.(dG-dT)n polymorphisms in world populations.

We have characterized eight dinucleotide (dC-dA)n.(dG-dT)n repeat loci located on human chromosome 13q in eight human populations and in a sample of chimpanzees. Even though there is substantial variation in allele frequencies at each locus, at a given locus the most frequent alleles are shared by all human populations. The level of heterozygosity is reduced in isolated or small populations, such as the Pehuenche Indians of Chile, the Dogrib of Canada, and the New Guinea highlanders. On the other hand, larger average heterozygosities are observed in large and cosmopolitan populations, such as the Sokoto population from Nigeria and German Caucasians. Conformity with Hardy-Weinberg equilibrium is generally observed at these loci, unless (a) a population is isolated or small or (b) the repeat motif of the locus is not perfect (e.g., D13S197). Multilocus genotype probabilities at these microsatellite loci do not show departure from the independence rule, unless the loci are closely linked. The allele size distributions at these (CA)n loci do not follow a strict single-step stepwise-mutation model. However, this features does not compromise the ability to detect population affinities, when these loci are used simultaneously. The microsatellite loci examined here are present and, with the exception of the locus D13S197, are polymorphic in the chimpanzees, showing an overlapping distribution of allele sizes with those observed in human populations.     

Variable (dG-dT)n.(dC-dA)n sequences in the porcine genome.

One of the more widely studied simple repeat sequences in the mammalian genome is the (dG-dT)n.(dC-dA)n dinucleotide repeat sequence. As these repeats are highly polymorphic and fairly evenly distributed in diverse mammalian genomes, they constitute a very powerful tool for genetic mapping in a wide variety of species. So far, the knowledge about repeat sequences in the porcine genome is sparse and only a few areas of this genome have been sequenced. We have isolated and characterized 108 porcine (dG-dT)n.(dC-dA)n sequences and studied the distribution of these, both by investigating random clones and by performing in situ hybridization. A remarkable correlation between humans and pigs was found with respect to the structure, to the number of repeat blocks, and to the chromosomal distribution.     

Informativeness of human (dC-dA)n · (dG-dT)n polymorphisms

Abundant human interspersed repetitive DNA sequences of the form (dC-dA)_n · (dG-dT)_n have been shown to exhibit length polymorphisms. Examination of over 100 human (dC-dA)_n · (dG-dT)_n sequences revealed that the sequences differed from each other both in numbers of repeats and in repeat sequence type. Using a set of precise classification rules, the sequences were divided into three categories: perfect repeat sequences without interruptions in the runs of CA or GT dinucleotides (64% of total), imperfect repeat sequences with one or more interruptions in the run of repeats (25%), and compound repeat sequences with adjacent tandem simple repeats of a different sequence (11%). Informativeness of (dC-dA)_n · (dG-dT)_n markers in the perfect sequence category was found to increase with increasing average numbers of repeats. PIC values ranged from 0 at about 10 or fewer repeats to above 0.8 for sequences with about 24 or more repeats. (dC-dA)_n · (dG-dT)_n polymorphisms in the imperfect sequence category showed lower informativeness than expected on the basis of the total numbers of repeats. The longest run of uninterrupted CA or GT repeats was found to be the best predictor of informativeness of (dC-dA)_n · (dG-dT)_n polymorphisms regardless of the repeat sequence category.     

(dC-dA)n.(dG-dT)n sequences have evolutionarily conserved chromosomal locations in Drosophila with implications for roles in chromosome structure and function.

In situ hybridization of (dC-dA)n.(dG-dT)n to the polytene chromosomes of Drosophila melanogaster reveals a clearly non-random distribution of chromosomal sites for this sequence. Sites are distributed over most euchromatic regions but the density of sites along the X chromosome is significantly higher than the density over the autosomes. All autosomes show approximately equal levels of hybridization except chromosome 4 which has no detectable stretches of (dC-dA)n.(dG-dT)n. Another striking feature is the lack of hybridization of the beta-heterochromatin of the chromocenter. The specific sites are conserved between different strains of D. melanogaster. The same overall chromosomal pattern of hybridization is seen for the other Drosophila species studied, including D. simulans, a sibling species with a much lower content of middle repetitive DNA, and D. virilis, a distantly related species. The evolutionary conservation of the distribution of (dC-dA)n.(dG-dT)n suggests that these sequences are of functional importance. The distribution patterns seen for D. pseudoobscura and D. miranda raise interesting speculations about function. In these species a chromosome equivalent to an autosomal arm of D. melanogaster has been translocated onto the X chromosome and acquired dosage compensation. In each species the new arm of the X also has a higher density of (dC-dA)n.(dG-dT)n similar to that seen on other X chromosomes. In addition to correlations with dosage compensation, the depletion of (dC-dA)n.(dG-dT)n in beta-heterochromatin and chromosome 4 may also be related to the fact that these regions do not normally undergo meiotic recombination.     

A rapidly evolving region in the immunoglobulin heavy chain loci of rat and mouse: postulated role of (dC-dA)n.(dG-dT)n sequences

The nucleotide sequences of the introns that are located between the C4 exon and the first membrane exon of mouse and rat immunoglobulin epsilon-chain genes have been determined. The rat intron sequence was found to contain four separate clusters of repetitive sequences all of which consisted of (dC-dA)n.(dG-dT)n dinucleotide repeats. A comparison between this chromosomal region in mouse and rat revealed four deletions or duplications, three of which have occurred inside or at the borders of the CA clusters. Rearrangements have occurred inside or at the borders of all four repeats after the evolutionary separation of mouse and rat. The sequence comparison reveals in addition a duplication, connected to the CA repeats, which has occurred early in evolution, before the evolutionary divergence of mouse and rat. These findings suggest that (dC-dA)n.(dG-dT)n sequences are potential targets for recombination events.     

Spectroscopic studies on acetylaminofluorene-modified (dT-dG)n . (dC-dA)n suggest a left-handed conformation

CD spectroscopy on the double-stranded strictly alternating dinucleotide polymer (dT-dG)n . (dC-dA)n partially modified by N-acetoxy-N-acetyl-2-aminofluorene suggests a left-handed conformation in concentrated NaCl solutions. Modification of the (dT-dG)n . (dC-dA)n polymer with acetylaminofluorene is required to promote formation of the left-handed helix since high salt concentrations and several other ionic conditions, which cause a similar transition for (dG-dC)n . (dG-dC)n, are ineffective. Furthermore, substitution of dC with 5-methyl dC in (dT-dG)n . (dC-dA)n does not facilitate formation of a left-handed helix, also in contrast to results found for (dG-dC)n . (dG-dC)n. A 62-base pair tract of almost perfectly alternating (dT-dG)n . (dC-dA)n from the 3'-side of the mouse kappa immunoglobulin gene modified with acetylaminofluorene undergoes the salt-induced transition to a left-handed helix when studied within a 140-base pair restriction fragment. High NaCl concentrations alone will not cause the transition for this 62-base pair tract in this fragment nor in the recombinant plasmid pRW777, which contains this fragment.     

Locations and contexts of sequences that hybridize to poly(dG-dT).(dC-dA) in mammalian ribosomal DNAs and two X-linked genes.

Sequences located several kilobases both 5' and 3' of the stably transcribed portion of several genes hybridize to radio-labeled pure fragments of the alternating sequence poly (dG-dT) (dC-dA) ["poly(GT)"]. The genes include the ribosomal DNA of mouse, rat, and human, and also human glucose-6-phosphate dehydrogenase (G6PD) and mouse hypoxanthine-guanine phosphoribosyl transferase (HPRT). HPRT has additional hybridizing sequences in introns. Fragments that include the hybridizing sequences and up to 300 bp of adjoining DNA show perfect runs of poly(GT) (greater than 30bp) in all but the human 5' region of rDNA, which shows a somewhat different alternating purine:pyrimidine sequence, poly(GTAT) (36bp). Within 150 bp of these sequences in various instances are found a number of other sequences reported to affect DNA conformation in model systems. Most marked is an enhancement of sequences matching at least 67% to the consensus binding sequence for topoisomerase II. Two to ten-fold less of such sequences were found in other sequenced portions of the nontranscribed spacer or in the transcribed portion of rDNA. The conservation of the locations of tracts of alternating purine:pyrimidine between evolutionarily diverse species is consistent with a possible functional role for these sequences.     

Unpaired structures in SCA10 (ATTCT)n.(AGAAT)n repeats

A number of human hereditary diseases have been associated with the instability of DNA repeats in the genome. Recently, spinocerebellar ataxia type 10 has been associated with expansion of the pentanucleotide repeat (ATTCT)(n).(AGAAT)(n) from a normal range of ten to 22 to as many as 4500 copies. The structural properties of this repeat cloned in circular plasmids were studied by a variety of methods. Two-dimensional gel electrophoresis and atomic force microscopy detected local DNA unpairing in supercoiled plasmids. Chemical probing analysis indicated that, at moderate superhelical densities, the (ATTCT)(n).(AGAAT)(n) repeat forms an unpaired region, which further extends into adjacent A+T-rich flanking sequences at higher superhelical densities. The superhelical energy required to initiate duplex unpairing is essentially length-independent from eight to 46 repeats. In plasmids containing five repeats, minimal unpairing of (ATTCT)(5).(AGAAT)(5) occurred while 2D gel analysis and chemical probing indicate greater unpairing in A+T-rich sequences in other regions of the plasmid. The observed experimental results are consistent with a statistical mechanical, computational analysis of these supercoiled plasmids. For plasmids containing 29 repeats, which is just above the normal human size range, flanked by an A+T-rich sequence, atomic force microscopy detected the formation of a locally condensed structure at high superhelical densities. However, even at high superhelical densities, DNA strands within the presumably compact A+T-rich region were accessible to small chemicals and oligonucleotide hybridization. Thus, DNA strands in this "collapsed structure" remain unpaired and accessible for interaction with other molecules. The unpaired DNA structure functioned as an aberrant replication origin, in that it supported complete plasmid replication in a HeLa cell extract. A model is proposed in which unscheduled or aberrant DNA replication is a critical step in the expansion mutation.     

Genetic recombination destabilizes (CTG)n.(CAG)n repeats in E. coli

The expansion of trinucleotide repeats has been implicated in 17 neurological diseases to date. Factors leading to the instability of trinucleotide repeat sequences have thus been an area of intense interest. Certain genes involved in mismatch repair, recombination, nucleotide excision repair, and replication influence the instability of trinucleotide repeats in both Escherichia coli and yeast. Using a genetic assay for repeat deletion in E. coli, the effect of mutations in the recA, recB, and lexA genes on the rate of deletion of (CTG)n.(CAG)n repeats of varying lengths were examined. The results indicate that mutations in recA and recB, which decrease the rate of recombination, had a stabilizing effect on (CAG)n.(CTG)n repeats decreasing the high rates of deletion seen in recombination proficient cells. Thus, recombination proficiency correlates with high rates of genetic instability in triplet repeats. Induction of the SOS system, however, did not appear to play a significant role in repeat instability, nor did the presence of triplet repeats in cells turn on the SOS response. A model is suggested where deletion during exponential growth may result from attempts to restart replication when paused at triplet repeats.     

Chromatin structure of the potential Z-forming sequence (dT-dG)n X (dC-dA)n. Evidence for an "alternating-B" conformation

The sequence (dT-dG)n X (dC-dA)n is the most abundant purine-pyrimidine dinucleotide repeat in eukaryotic genomes. This sequence and certain others that contain alternating purine-pyrimidine residues have been shown to adopt the left-handed, Z-DNA conformation in vitro when subjected to negative torsional stress or elevated ionic strengths. We have asked whether (dT-dG)n X (dC-dA)n tracts exist in topologically constrained Z-form structures in vivo by examining the chromatin organization of these sequences in cultured mouse cell nuclei. We find that these elements are quantitatively packaged into typical core particles which are embedded in canonical polynucleosomal arrays. In addition, these sequences neither flank nor reside within regions of chromatin that are preferentially sensitive to S1 nuclease. These characteristics suggest that these tracts do not exist predominantly in the Z-form in vivo. Furthermore, employing techniques that permit prominent hybridization to DNA fragments as short as 18 bases, we provide evidence that in vivo, most (dT-dG)n X (dC-dA)n elements instead adopt an "alternating-B" conformation on the nucleosomal surface.     

Stabilization of (dG-dC)n.(dG-dC)n in the Z conformation by a crosslinking reaction

(dG-dC)n.(dG-dC) was converted to the Z conformer by heating in the presence of Mn++n. Reaction of this preparation with the crosslinking reagent, DL-diepoxybutane (DEB), stabilized this conformer so that it retained its structure even when returned to conditions that favored reversion to the B conformation. Treatment of the crosslinked Z conformer with periodate caused scission of the crosslink, allowing reversion to the B conformer. Reaction of (dG-dC)n.(dG-dC)n in the B conformation with DEB did not prevent conversion to the Z conformer in 4M NaC1; dialysis of the high salt solution against low ionic strength buffer allowed return to the B conformer. The Z in equilibrium B transitions were followed by circular dichroism studies and immunochemical procedures. The results suggest the feasibility of stabilizing Z sequences of DNA in chromatin by crosslinking, so that they could then be identified after DNA isolation.     

Long (dA)n.(dT)n tracts can form intramolecular triplexes under superhelical stress.

Plasmids containing long tracts of (dA)n.(dT)n have been prepared and their conformations examined in linear and supercoiled DNA using a series of chemical and enzymic probes which are known to be sensitive to unusual DNA structures. Under superhelical stress and in the presence of magnesium the sequence T69.A69 adopts a conformation at pH 8.0 consistent with the formation of an intramolecular DNA triplex. Site specific cleavage of the supercoiled plasmid by single-strand specific nucleases occurs within the A.T insert; the 5'-end of the purine strand is sensitive to reaction with diethylpyrocarbonate while the central 5-6 bases of the pyrimidine strand are reactive to osmium tetroxide. By contrast shorter inserts of A33.T33 and A23.T23 do not appear to form unusual structures.     

Magnesium-dependent supercoiling-induced transition in (dG)n.(dC)n stretches and formation of a new G-structure by (dG)n strand.

Plasmids containing (dG)27.(dC)27 inserts (pPG27), (dG)37.(dC)37 inserts (pPG37), and (dG)24C(dG)21.(dC)24G(dC)21 inserts (pPG46C) were constructed for the study of structural transitions within (dG)n.(dC)n stretches. Two-dimensional gel electrophoresis has shown that a Mg2+-dependent supercoiling-induced structural transition takes place at pH 8 in plasmid pPG46C. The transition occurs at -0=0.06 and involves a supercoiling release corresponding to 5 superhelical turns. After denaturation of the restriction fragments containing (dG)n.(dC)n inserts, the strands do not renature completely and (dG)n-containing strand migrates in PAGE much faster than the (dC)n-containing one. Chemical modification experiments with the (dG)n-strand have revealed the periodic nature of the protection of guanines against dimethyl sulfate methylation. The (dG)n strand in the presence of Mg2+ forms complexes with the complementary (dC)n strand, which differ from the native duplex in mobility. We believe these effects to be due to the formation of an intrastrand structure within the (dG)n strand stabilized by G.G interactions (we called it G-structure), which in the presence of Mg2+ forms an interstrand complex. with the (dC)n strand.     

Recognition of (dG)n.(dC)n sequences by endonuclease G. Characterization of the calf thymus nuclease

We report the purification of endonuclease G (Ruiz-Carrillo, A., and Renaud, J. (1987) EMBO J. 6, 401-407) from calf thymus nuclei and whole tissue. The enzyme has been enriched 29,000-fold, and the activity was unambiguously identified with a 26-kDa protein after renaturation following sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The native nuclease behaves as a 50-kDa species by gel filtration, suggesting that it is composed of two subunits, presumably identical. In terms of absolute amounts, endonuclease G (endo G) is a nuclear enzyme although it was also detected in purified mitochondria. Endo G is highly specific for (dG)n.(dC)n tracts in DNA, nicking either strand of relaxed substrates with similar kinetics. The sensitivity of the homopolymer tracts is proportional to their length (from n = 8 to 29), insofar as the flanking sequences are constant. However, the overall rate of cleavage is influenced by the composition of the flanking DNA. Minor cleavage sites contain shorter (dG)n.(dC)n clusters (n = 3-7). Endo G efficiently cleaves (dC)n but not (dG)n runs in single-stranded DNA, suggesting that it may recognize an asymmetric strand conformation of the homopolymer tracts. Endo G does not recognize other homo(co)-polymer sequences or cruciform structures in DNA.     

Solution structures of the individual single strands of the fragile X DNA triplets (GCC)n.(GGC)n.

Three-dimensional structures of the fragile X triplet repeats (GCC)n and (GGC)n are derived by using one- dimensional/two-dimensional NMR. Under a wide range of solution conditions (10-150 mM NaCl,pH6-7)(GCC)5-7 strands form exclusively slipped hairpins with a 3' overhanging C. The slipped hairpins of (GCC)n strands show the following structural characteristics: (i) maximization of Watson-Crick G.C pairs; (ii) formation of C.C mispairs at the CpG steps in the stem; (iii) C2'-endo, anti conformations for all the nucleotides. The ability of (GCC)n strands to form hairpin structures more readily than complementary (GGC)n strands suggests preferential slippage during replication and subsequent expansion of the (GCC)n strands. In addition, the C.C. mispairs at the CpG site of (GCC)n hairpins account for their exceptional substrate efficiencies for human methyltransferase. Gel electrophoresis data show that (GGC)n strands form both hairpin and mismatched duplex structures in 10-150 mM NaCl (ph 6-7) for n < 10, but for n > or + 11 hairpin structures are exclusively present. However, (GGC)n strands remain predominantly in the duplex state for n=4-11 under NMR solution conditions, which require DNA concentrations 100- to 1000-fold higher than in gel electrophoresis. NMR analyses of [(GGC)n]2 duplexes for n=4-6 show the presence of Watson-Crick G.C and mismatched G anti G syn pairs. The mismatches adjacent to the CpG step introduce local structural flexibility in these duplexes. Similar structural properties are also expected in the stem of the hairpins formed by (GGC)n strands.     

Right-handed and left-handed helices of poly(dA-dC) X (dG-dT).

The secondary structures of poly(dA-dC) X (dG-dT) were studied using CD and IR spectroscopies. We give spectroscopic evidence of secondary structure transitions of poly(dA-dC) X (dG-dT) from a B to a Z-like helix, induced by transition metal ions (Ni2+) in presence of high concentrations of Cs+ and Na+. In the presence of Na+, the B in equilibrium Z transition occurs at any temperature, whereas premelting conditions are required in presence of Cs+. For these two alkali ions the Z-like form is only induced by Ni2+ ions through their specific interactions at N7 of purines, under conditions of low water activity due to the high alkali salt concentration. We also show that the CD spectrum obtained in presence of Cs+ ions and characterized by a negative band at 275 nm, cannot be interpreted in terms of Z-like left-handed helix but reflects a modified B right-handed helix.     

The ubiquitous potential Z-forming sequence of eucaryotes, (dT-dG)n . (dC-dA)n, is not detectable in the genomes of eubacteria, archaebacteria, or mitochondria.

The potential Z-forming sequence (dT-dG)n . (dC-dA)n is an abundant, interspersed repeat element that is ubiquitous in eucaryotic nuclear genomes. We report that in contrast to eucaryotic nuclear DNA, the genomes of eubacteria, archaebacteria, and mitochondria lack this sequence, since even a single tract of greater than or equal to 14 base pairs in length is not detectable through either hybridization or sequence analysis. Interestingly, the phylogenetic distribution of the (dT-dG)n . (dC-dA)n repeat exhibits a striking parallel to that of (dT-dC)n . (dG-dA)n, but not to other homocopolymeric sequences such as (dC-dG)n . (dC-dG)n or (dT-dA)n . (dT-dA)n.