DNA sequence insertion and evolutionary variation in gene regulation.

Current evidence on the long-term evolutionary effect of insertion of sequence elements into gene regions is reviewed, restricted to cases where a sequence derived from a past insertion participates in the regulation of expression of a useful gene. Ten such examples in eukaryotes demonstrate that segments of repetitive DNA or mobile elements have been inserted in the past in gene regions, have been preserved, sometimes modified by selection, and now affect control of transcription of the adjacent gene. Included are only examples in which transcription control was modified by the insert. Several cases in which merely transcription initiation occurred in the insert were set aside. Two of the examples involved the long terminal repeats of mammalian endogenous retroviruses. Another two examples were control of transcription by repeated sequence inserts in sea urchin genomes. There are now six published examples in which Alu sequences were inserted long ago into human gene regions, were modified, and now are central in control/enhancement of transcription. The number of published examples of Alu sequences affecting gene control has grown threefold in the last year and is likely to continue growing. Taken together, all of these examples show that the insertion of sequence elements in the genome has been a significant source of regulatory variation in evolution.

[Leu5]enkephalin-encoding sequences are targets for a specific DNA-binding factor.

A DNA-binding factor with high affinity and specificity for the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes has been characterized. The factor has the highest affinity for the [Leu5]-enkephalin-encoding sequence in the dynorphin B-encoding region of the prodynorphin gene, has relatively high affinity for other [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes, but has no apparent affinity for similar DNA sequences coding for [Met5]-enkephalin in the prodynorphin or proopiomelanocortin genes. The factor has been named [Leu5]enkephalin-encoding sequence DNA-binding factor (LEF). LEF has a nuclear localization and is composed of three subunits of about 60, 70, and 95 kDa, respectively. The highest levels were observed in rat testis, cerebellum, and spleen and were generally higher in late embryonal compared to newborn or adult animals. LEF activity was also recorded in human clonal tumor cell lines. LEF inhibited the transcription of reporter genes in artificial gene constructs where a [Leu5]enkephalin-encoding DNA fragment had been inserted between the transcription initiation site and the coding region of the reporter genes. These observations suggest that the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes also have regulatory functions realized through interaction with a specific DNA-binding factor.

Efficient high-resolution genetic mapping of mouse interspersed repetitive sequence PCR products, toward integrated genetic and physical mapping of the mouse genome.

The ability to carry out high-resolution genetic mapping at high throughput in the mouse is a critical rate-limiting step in the generation of genetically anchored contigs in physical mapping projects and the mapping of genetic loci for complex traits. To address this need, we have developed an efficient, high-resolution, large-scale genome mapping system. This system is based on the identification of polymorphic DNA sites between mouse strains by using interspersed repetitive sequence (IRS) PCR. Individual cloned IRS PCR products are hybridized to a DNA array of IRS PCR products derived from the DNA of individual mice segregating DNA sequences from the two parent strains. Since gel electrophoresis is not required, large numbers of samples can be genotyped in parallel. By using this approach, we have mapped > 450 polymorphic probes with filters containing the DNA of up to 517 backcross mice, potentially allowing resolution of 0.14 centimorgan. This approach also carries the potential for a high degree of efficiency in the integration of physical and genetic maps, since pooled DNAs representing libraries of yeast artificial chromosomes or other physical representations of the mouse genome can be addressed by hybridization of filter representations of the IRS PCR products of such libraries.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.

Inhibition of RNA polymerase II transcription in human cells by synthetic DNA-binding ligands

Sequence-specific DNA-binding small molecules that can permeate human cells potentially could regulate transcription of specific genes. Multiple cellular DNA-binding transcription factors are required by HIV type 1 for RNA synthesis. Two pyrrole–imidazole polyamides were designed to bind DNA sequences immediately adjacent to binding sites for the transcription factors Ets-1, lymphoid-enhancer binding factor 1, and TATA-box binding protein. These synthetic ligands specifically inhibit DNA-binding of each transcription factor and HIV type 1 transcription in cell-free assays. When used in combination, the polyamides inhibit virus replication by >99% in isolated human peripheral blood lymphocytes, with no detectable cell toxicity. The ability of small molecules to target predetermined DNA sequences located within RNA polymerase II promoters suggests a general approach for regulation of gene expression, as well as a mechanism for the inhibition of viral replication.

A zinc finger directory for high-affinity DNA recognition

We have used two monovalent phage display libraries containing variants of the Zif268 DNA-binding domain to obtain families of zinc fingers that bind to alterations in the last 4 bp of the DNA sequence of the Zif268 consensus operator, GCG TGGGCG. Affinity selection was performed by altering the Zif268 operator three base pairs at a time, and simultaneously selecting for sets of 16 related DNA sequences. In this way, only four experiments were required to select for all possible 64 combinations of DNA triplet sequences. The results show that (i) for high-affinity DNA binding in the range observed for the Zif268 wild-type complex (Kd = 0.5–5 nM), finger 1 specifically requires the arginine at the carboxy terminus of its recognition helix that forms a bidentate hydrogen-bond with the guanine base (G) in the crystal structure of Zif268 complexed to its DNA operator sequence GCG TGG GCG; (ii) when the guanine base (G) is replaced by A, C, or T, a lower-affinity family (Kd ⩾ 50 nM) can be detected that shows an overall tendency to bind G-rich DNA; (iii) the residues at position 2 on the finger 2 recognition helix do not appear to interact strongly with the complementary 5′ base in the finger 1 binding site; and (iv) unexpected substitutions at the amino terminus of finger 1 can occasionally result in specificity for the 3′ base in the finger 1 binding site. A DNA recognition directory was constructed for high-affinity zinc fingers that recognize all three bases in a DNA triplet for seven sequences of the type GNN. Similar approaches may be applied to other zinc fingers to broaden the scope of the directory.

Conserved characteristics of heterochromatin-forming DNA at the 15q11-q13 imprinting center

Nuclear matrix binding assays (NMBAs) define certain DNA sequences as matrix attachment regions (MARs), which often have cis-acting epigenetic regulatory functions. We used NMBAs to analyze the functionally important 15q11-q13 imprinting center (IC). We find that the IC is composed of an unusually high density of MARs, located in close proximity to the germ line elements that are proposed to direct imprint switching in this region. Moreover, we find that the organization of MARs is the same at the homologous mouse locus, despite extensive divergence of DNA sequence. MARs of this size are not usually associated with genes but rather with heterochromatin-forming areas of the genome. In contrast, the 15q11-q13 region contains multiple transcribed genes and is unusual for being subject to genomic imprinting, causing the maternal chromosome to be more transcriptionally silent, methylated, and late replicating than the paternal chromosome. We suggest that the extensive MAR sequences at the IC are organized as heterochromatin during oogenesis, an organization disrupted during spermatogenesis. Consistent with this model, multicolor fluorescence in situ hybridization to halo nuclei demonstrates a strong matrix association of the maternal IC, whereas the paternal IC is more decondensed, extending into the nuclear halo. This model also provides a mechanism for spreading of the imprinting signal, because heterochromatin at the IC on the maternal chromosome may exert a suppressive position effect in cis. We propose that the germ line elements at the 15q11-q13 IC mediate their effects through the candidate heterochromatin-forming DNA identified in this study.

A knob-associated tandem repeat in maize capable of forming fold-back DNA segments: Are chromosome knobs megatransposons?

A class of tandemly repeated DNA sequences (TR-1) of 350-bp unit length was isolated from the knob DNA of chromosome 9 of Zea mays L. Comparative fluorescence in situ hybridization revealed that TR-1 elements are also present in cytologically detectable knobs on other maize chromosomes in different proportions relative to the previously described 180-bp repeats. At least one knob on chromosome 4 is composed predominantly of the TR-1 repeat. In addition, several small clusters of the TR-1 and 180-bp repeats have been found in different chromosomes, some not located in obvious knob heterochromatin. Variation in restriction fragment fingerprints and copy number of the TR-1 elements was found among maize lines and among maize chromosomes. TR-1 tandem arrays up to 70 kilobases in length can be interspersed with stretches of 180-bp tandem repeat arrays. DNA sequence analysis and restriction mapping of one particular stretch of tandemly arranged TR-1 units indicate that these elements may be organized in the form of fold-back DNA segments. The TR-1 repeat shares two short segments of homology with the 180-bp repeat. The longest of these segments (31 bp; 64% identity) corresponds to the conserved region among 180-bp repeats. The polymorphism and complex structure of knob DNA suggest that, similar to the fold-back DNA-containing giant transposons in Drosophila, maize knob DNA may have some properties of transposable elements.

Ash1p is a site-specific DNA-binding protein that actively represses transcription

ASH1 encodes a protein that is localized specifically to the daughter cell nucleus, where it has been proposed to repress transcription of the HO gene. Using Ash1p purified from baculovirus-infected insect cells, we have shown that Ash1p binds specific DNA sequences in the HO promoter. DNase I protection analyses showed that Ash1p recognizes a consensus sequence, YTGAT. Mutation of this consensus abolishes Ash1p DNA binding in vitro. We have shown that Ash1p requires an intact zinc-binding domain in its C terminus for repression of HO in vivo and that this domain may be involved in DNA binding. A heterologous DNA-binding domain fused to an N-terminal segment of Ash1p functions as an active repressor of transcription. Our studies indicate that Ash1p is a DNA-binding protein of the GATA family with a separable transcriptional repression domain.

Target DNA chromatinization modulates nicking by L1 endonuclease

L1 elements are human transposons which replicate via an RNA intermediate. At least 15% of the human genome is composed of L1 sequence. An important initial step in the transposition reaction is nicking of the genomic DNA by L1 endonuclease (L1 EN). In vivo much of the genome exists in the form of chromatin or is undergoing biochemical transactions such as transcription, replication or repair, which may alter the accessibility of the L1 transposition machinery to DNA. To investigate this possibility we have examined the effect of substrate chromatinization on the ability of L1 EN to nick DNA. We find that DNA incorporated into nucleosomes is generally refractory to nicking by L1 EN. Interestingly, nicking of a minority of DNA sequences is enhanced when included in chromatin. Thus, dynamic epigenetic factors such as chromatinization are likely to influence the relatively permanent placement of L1 and other retroelements in the human genome.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Visualization of oligonucleotide probes and point mutations in interphase nuclei and DNA fibers using rolling circle DNA amplification

Rolling circle amplification (RCA) is a surface-anchored DNA replication reaction that can be exploited to visualize single molecular recognition events. Here we report the use of RCA to visualize target DNA sequences as small as 50 nts in peripheral blood lymphocytes or in stretched DNA fibers. Three unique target sequences within the cystic fibrosis transmembrane conductance regulator gene could be detected simultaneously in interphase nuclei, and could be ordered in a linear map in stretched DNA. Allele-discriminating oligonucleotide probes in conjunction with RCA also were used to discriminate wild-type and mutant alleles in the cystic fibrosis transmembrane conductance regulator, p53, BRCA-1, and Gorlin syndrome genes in the nuclei of cultured cells or in DNA fibers. These observations demonstrate that signal amplification by RCA can be coupled to nucleic acid hybridization and multicolor fluorescence imaging to detect single nucleotide changes in DNA within a cytological context or in single DNA molecules. This provides a means for direct physical haplotyping and the analysis of somatic mutations on a cell-by-cell basis.

Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5,000–50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.

Exploring the DNA-binding specificities of zinc fingers with DNA microarrays

A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA–protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.