Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

DNA incisions stimulate genetic instability of triplet repeat sequences related to human hereditary neurological diseases

The molecular mechanisms responsible for the expansion and deletion of trinucleotide repeat sequences (TRS) are the focus of our studies. Several hereditary neurological diseases including Huntington's disease, myotonic dystrophy, and fragile X syndrome are associated with the instability of TRS. Using the well defined and controllable model system of Escherichia coli, the influences of three types of DNA incisions on genetic instability of CTG•CAG repeats were studied: DNA double-strand breaks (DSB), single-strand nicks, and single-strand gaps. The DNA incisions were generated in pUC19 derivatives by in vitro cleavage with restriction endonucleases. The cleaved DNA was then transformed into E. coli parental and mutant strains. Double-strand breaks induced deletions throughout the TRS region in an orientation dependent manner relative to the origin of replication. The extent of instability was enhanced by the repeat length and sequence (CTG•CAG vs. CGG•CCG). Mutations in recA and recBC increased deletions, mutations in recF stabilized the TRS, whereas mutations in ruvA had no effect. DSB were repaired by intramolecular recombination, versus an intermolecular gene conversion or crossover mechanism. 30 nt gaps formed a distinct 30 nt deletion product, whereas single strand nicks and gaps of 15 nts did not induce expansions or deletions. Formation of this deletion product required the CTG•CAG repeats to be present in the single-stranded region and was stimulated by E. coli DNA ligase, but was not dependent upon the RecFOR pathway. Models are presented to explain the DSB induced instabilities and formation of the 30 nucleotide deletion product. In addition to the in vitro creation of DSBs, several attempts to generate this incision in vivo with the use of EcoR I restriction modification systems were conducted. ^

Roles of mismatch repair proteins in the recognition and processing of DNA interstrand crosslinks in human cells

DNA interstrand crosslinks (ICLs) are among the most toxic type of damage to a cell. Many ICL-inducing agents are widely used as therapeutic agents, e.g. cisplatin, psoralen. A bettor understanding of the cellular mechanism that eliminates ICLs is important for the improvement of human health. However, ICL repair is still poorly understood in mammals. Using a triplex-directed site-specific ICL model, we studied the roles of mismatch repair (MMR) proteins in ICL repair in human cells. We are also interested in using psoralen-conjugated triplex-forming oligonucleotides (TFOs) to direct ICLs to a specific site in targeted DNA and in the mammalian genomes. ^ MSH2 protein is the common subunit of two MMR recognition complexes, and MutSα and MutSβ. We showed that MSH2 deficiency renders human cell hypersensitive to psoralen ICLs. MMR recognition complexes bind specifically to triplex-directed psoralen ICLs in vitro. Together with the fact that psoralen ICL-induced repair synthesis is dramatically decreased in MSH2 deficient cell extracts, we demonstrated that MSH2 function is critical for the recognition and processing of psoralen ICLs in human cells. Interestingly, lack of MSH2 does not reduce the level of psoralen ICL-induced mutagenesis in human cells, suggesting that MSH2 does not contribute to error-generating repair of psoralen ICLs, and therefore, may represent a novel error-free mechanism for repairing ICLs. We also studied the role of MLH1, anther key protein in MMR, in the processing of psoralen ICLs. MLH1-deficient human cells are more resistant to psoralen plus UVA treatment. Importantly, MLH1 function is not required for the mutagenic repair of psoralen ICLs, suggesting that it is not involved in the error-generating repair of this type of DNA damage in human cells. ^ These are the first data indicating mismatch repair proteins may participate in a relatively error-free mechanism for processing psoralen ICL in human cells. Enhancement of MMR protein function relative to nucleotide excision repair proteins may reduce the mutagenesis caused by DNA ICLs in humans. ^ In order to specifically target ICLs to mammalian genes, we identified novel TFO target sequences in mouse and human genomes. Using this information, many critical mammalian genes can now be targeted by TFOs.^

Human DNA ligase and DNA polymerase as molecular targets for heavy metals and anticancer drugs

DNA ligase and DNA polymerase play important roles in DNA replication, repair, and recombination. Frequencies of spontaneous and chemical- and physical-induced mutations are correlated to the fidelity of DNA replication. This dissertation elucidates the mechanisms of the DNA ligation reaction by DNA ligases and demonstrates that human DNA ligase I and DNA polymerase $\alpha$ are the molecular targets for two metal ions, Zn$\sp{2+}$ and Cd$\sp{2+},$ and an anticancer drug, F-ara-ATP.^ Human DNA ligases were purified to homogeneity and their AMP binding domains were mapped. Although their AMP-binding domains are similar, there could be difference between the two ligases in their DNA binding domains.^ The formation of the AMP-DNA intermediate and the successive ligation reaction by human DNA ligases were analyzed. Both reactions showed their substrate specificity for ligases I and II, required Mg2+, and were inhibited by ATP.^ A protein inhibitor from HeLa cells and specific for human DNA ligase I but not ligase II and T4 ligase was discovered. It reversibly inhibited DNA ligation activity but not the AMP-binding activity due to the formation of a reversible ligase I-inhibitor complex.^ F-ara-ATP inhibited human DNA ligase I activity by competing with ATP for the AMP-binding site of DNA ligase I, forming a ligase I-F-ara-AMP complex, as well as when it was incorporated at 3$\sp\prime$-terminus of DNA nick by DNA polymerase $\alpha.$^ All steps of the DNA ligation reaction were inhibited by Zn$\sp{2+}$ and Cd$\sp{2+}$ in a concentration-dependent manner. Both ions did not show the ability to change the fidelity of DNA ligation reaction catalyzed by human DNA ligase I. However, Zn$\sp{2+}$ and Cd$\sp{2+}$ showed their contradictory effects on the fidelity of the reaction by human DNA polymerase $\alpha.$ Zn$\sp{2+}$ decreased the frequency of misinsertion but less affected that of mispair extension. On the contrary, Cd$\sp{2+}$ increased the frequencies of both misinsertion and mispair extension at very low concentration. Our data provided strong evidence in the molecular mechanisms for the mutagenicity of zinc and cadmium, and were comparable with the results previously reported. ^

A MOLECULAR APPROACH FOR DETECTING GENETIC DAMAGE USING THE ECO R1 FAMILY OF HUMAN DNA (SEQUENCING)

A clone of the primary Eco R1 family of human DNA sequences has been used as an indicator sequence for detecting alterations induced by a toxic agent. Specific clones of this family have been examined and compared to the consensus sequence to determine the normal variability of this family. Though variations were observed, data indicated that such clones can be used to study induced DNA modifications. This DNA was exposed to the toxic agent dimethyl sulfate under various conditions and a distinct pattern of aberrations was shown to occur. It is suggested that this approach be used to characterize patterns of damage induced by various agents in the ultimate development of a system capable of monitoring human genotoxic exposure. ^

The C57Bl/6 mouse serves as a suitable model of human skeletal muscle mitochondrial function

Función mitocondrial en ratas

Crystal structure of human p32, a doughnut-shaped acidic mitochondrial matrix protein

Human p32 (also known as SF2-associated p32, p32/TAP, and gC1qR) is a conserved eukaryotic protein that localizes predominantly in the mitochondrial matrix. It is thought to be involved in mitochondrial oxidative phosphorylation and in nucleus–mitochondrion interactions. We report the crystal structure of p32 determined at 2.25 Å resolution. The structure reveals that p32 adopts a novel fold with seven consecutive antiparallel β-strands flanked by one N-terminal and two C-terminal α-helices. Three monomers form a doughnut-shaped quaternary structure with an unusually asymmetric charge distribution on the surface. The implications of the structure on previously proposed functions of p32 are discussed and new specific functional properties are suggested.

A human Cds1-related kinase that functions downstream of ATM protein in the cellular response to DNA damage

Checkpoints maintain the order and fidelity of the eukaryotic cell cycle, and defects in checkpoints contribute to genetic instability and cancer. Much of our current understanding of checkpoints comes from genetic studies conducted in yeast. In the fission yeast Schizosaccharomyces pombe (Sp), SpRad3 is an essential component of both the DNA damage and DNA replication checkpoints. The SpChk1 and SpCds1 protein kinases function downstream of SpRad3. SpChk1 is an effector of the DNA damage checkpoint and, in the absence of SpCds1, serves an essential function in the DNA replication checkpoint. SpCds1 functions in the DNA replication checkpoint and in the S phase DNA damage checkpoint. Human homologs of both SpRad3 and SpChk1 but not SpCds1 have been identified. Here we report the identification of a human cDNA encoding a protein (designated HuCds1) that shares sequence, structural, and functional similarity to SpCds1. HuCds1 was modified by phosphorylation and activated in response to ionizing radiation. It was also modified in response to hydroxyurea treatment. Functional ATM protein was required for HuCds1 modification after ionizing radiation but not after hydroxyurea treatment. Like its fission yeast counterpart, human Cds1 phosphorylated Cdc25C to promote the binding of 14-3-3 proteins. These findings suggest that the checkpoint function of HuCds1 is conserved in yeast and mammals.

MED1, a novel human methyl-CpG-binding endonuclease, interacts with DNA mismatch repair protein MLH1

The DNA mismatch repair (MMR) is a specialized system, highly conserved throughout evolution, involved in the maintenance of genomic integrity. To identify novel human genes that may function in MMR, we employed the yeast interaction trap. Using the MMR protein MLH1 as bait, we cloned MED1. The MED1 protein forms a complex with MLH1, binds to methyl-CpG-containing DNA, has homology to bacterial DNA repair glycosylases/lyases, and displays endonuclease activity. Transfection of a MED1 mutant lacking the methyl-CpG-binding domain (MBD) is associated with microsatellite instability (MSI). These findings suggest that MED1 is a novel human DNA repair protein that may be involved in MMR and, as such, may be a candidate eukaryotic homologue of the bacterial MMR endonuclease, MutH. In addition, these results suggest that cytosine methylation may play a role in human DNA repair.

Assembly, subunit composition, and footprint of human DNA repair excision nuclease

The assembly and composition of human excision nuclease were investigated by electrophoretic mobility shift assay and DNase I footprinting. Individual repair factors or any combination of up to four repair factors failed to form DNA–protein complexes of high specificity and stability. A stable complex of high specificity can be detected only when XPA/RPA, transcription factor IIH, XPC⋅HHR23B, and XPG and ATP are present in the reaction mixture. The XPF⋅ERCC1 heterodimer changes the electrophoretic mobility of the DNA–protein complex formed with the other five repair factors, but it does not confer additional specificity. By using proteins with peptide tags or antibodies to the repair factors in electrophoretic mobility shift assays, it was found that XPA, replication protein A, transcription factor IIH, XPG, and XPF⋅excision repair cross-complementing 1 but not XPC⋅HHR23B were present in the penultimate and ultimate dual incision complexes. Thus, it appears that XPC⋅HHR23B is a molecular matchmaker that participates in the assembly of the excision nuclease but is not present in the ultimate dual incision complex. The excision nuclease makes an assymmetric DNase I footprint of ≈30 bp around the damage and increases the DNase I sensitivity of the DNA on both sides of the footprint.

A human homolog of the Saccharomyces cerevisiae REV3 gene, which encodes the catalytic subunit of DNA polymerase ζ

To get a better understanding of mutagenic mechanisms in humans, we have cloned and sequenced the human homolog of the Saccharomyces cerevisiae REV3 gene. The yeast gene encodes the catalytic subunit of DNA polymerase ζ, a nonessential enzyme that is thought to carry out translesion replication and is responsible for virtually all DNA damage-induced mutagenesis and the majority of spontaneous mutagenesis. The human gene encodes an expected protein of 3,130 residues, about twice the size of the yeast protein (1,504 aa). The two proteins are 29% identical in an amino-terminal region of ≈340 residues, 39% identical in a carboxyl-terminal region of ≈850 residues, and 29% identical in a 55-residue region in the middle of the two genes. The sequence of the expected protein strongly predicts that it is the catalytic subunit of a DNA polymerase of the pol ζ type; the carboxyl-terminal domain possesses, in the right order, the six motifs characteristic of eukaryotic DNA polymerases, most closely resembles yeast pol ζ among all polymerases in the GenBank database, and is different from the human α, δ, and ɛ enzymes. Human cells expressing high levels of an hsREV3 antisense RNA fragment grow normally, but show little or no UV-induced mutagenesis and are slightly more sensitive to killing by UV. The human gene therefore appears to carry out a function similar to that of its yeast counterpart.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

Detection of 400-year-old Yersinia pestis DNA in human dental pulp: An approach to the diagnosis of ancient septicemia

Ancient septicemic plague epidemics were reported to have killed millions of people for 2 millenniums. However, confident diagnosis of ancient septicemia solely on the basis of historical clinical observations is not possible. The lack of suitable infected material has prevented direct demonstration of ancient septicemia; thus, the history of most infections such as plague remains hypothetical. The durability of dental pulp, together with its natural sterility, makes it a suitable material on which to base such research. We hypothesized that it would be a lasting refuge for Yersinia pestis, the plague agent. DNA extracts were made from the dental pulp of 12 unerupted teeth extracted from skeletons excavated from 16th and 18th century French graves of persons thought to have died of plague (“plague teeth”) and from 7 ancient negative control teeth. PCRs incorporating ancient DNA extracts and primers specific for the human β-globin gene demonstrated the absence of inhibitors in these preparations. The incorporation of primers specific for Y. pestis rpoB (the RNA polymerase β-subunit-encoding gene) and the recognized virulence-associated pla (the plasminogen activator-encoding gene) repeatedly yielded products that had a nucleotide sequence indistinguishable from that of modern day isolates of the bacterium. The specific pla sequence was obtained from 6 of 12 plague skeleton teeth but 0 of 7 negative controls (P < 0.034, Fisher exact test). A nucleic acid-based confirmation of ancient plague was achieved for historically identified victims, and we have confirmed the presence of the disease at the end of 16th century in France. Dental pulp is an attractive target in the quest to determine the etiology of septicemic illnesses detected in ancient corpses. Molecular techniques could be applied to this material to resolve historical outbreaks.

A Z-DNA binding domain present in the human editing enzyme, double-stranded RNA adenosine deaminase

Editing of RNA changes the read-out of information from DNA by altering the nucleotide sequence of a transcript. One type of RNA editing found in all metazoans uses double-stranded RNA (dsRNA) as a substrate and results in the deamination of adenosine to give inosine, which is translated as guanosine. Editing thus allows variant proteins to be produced from a single pre-mRNA. A mechanism by which dsRNA substrates form is through pairing of intronic and exonic sequences before the removal of noncoding sequences by splicing. Here we report that the RNA editing enzyme, human dsRNA adenosine deaminase (DRADA1, or ADAR1) contains a domain (Zα) that binds specifically to the left-handed Z-DNA conformation with high affinity (KD = 4 nM). As formation of Z-DNA in vivo occurs 5′ to, or behind, a moving RNA polymerase during transcription, recognition of Z-DNA by DRADA1 provides a plausible mechanism by which DRADA1 can be targeted to a nascent RNA so that editing occurs before splicing. Analysis of sequences related to Zα has allowed identification of motifs common to this class of nucleic acid binding domain.

The solution structure of the Zα domain of the human RNA editing enzyme ADAR1 reveals a prepositioned binding surface for Z-DNA

Double-stranded RNA deaminase I (ADAR1) contains the Z-DNA binding domain Zα. Here we report the solution structure of free Zα and map the interaction surface with Z-DNA, confirming roles previously assigned to residues by mutagenesis. Comparison with the crystal structure of the (Zα)2/Z-DNA complex shows that most Z-DNA contacting residues in free Zα are prepositioned to bind Z-DNA, thus minimizing the entropic cost of binding. Comparison with homologous (α+β)helix–turn–helix/B-DNA complexes suggests that binding of Zα to B-DNA is disfavored by steric hindrance, but does not eliminate the possibility that related domains may bind to both B- and Z-DNA.