DNA cytosine methylation and heat-induced deamination

The heat-induced conversion of 5-methylcytosine (m⁵C) residues to thymine residues and of cytosine to uracil residues in single-stranded DNA was studied. The calculated rates for deamination at 37°C and pH 7.4 were 9.5×10^–10 and 2.1×10^–10 sec^–1, respectively. N⁴-Methyldeoxycytidine, which is in the DNA of certain thermophilic bacteria, was more heat-resistant than was deoxycytidine and much more than was 5-methyldeoxycytidine. Thermophilic bacteria which contain N⁴-methylcytosine rather than m⁵C in their genomes may thereby largely avoid heat-induced mutation due to deamination, which is incurred by the many organisms that contain m⁵C in their DNA.     

Potential artifacts in the measurement of DNA deamination

Attack on DNA by some reactive nitrogen species results in deamination of adenine and guanine, leading to the formation of hypoxanthine and xanthine, respectively. Published levels of these products in cellular DNA have varied widely. Although these two deamination products are often measured by GC-MS analysis, the procedure of acid hydrolysis to release DNA bases for derivatization poses a risk of artifactual deamination of the DNA. In this study, we demonstrated the artifactual formation of these two deamination products during acid hydrolysis and hence developed a method for detecting and measuring 2'-deoxyinosine, the nucleoside of hypoxanthine. Our assay for 2'-deoxyinosine employs nuclease P1 and alkaline phosphatase to achieve release of the nucleosides from DNA, followed by HPLC prepurification with subsequent GC-MS analysis of the nucleosides. This assay detected an increase in the levels of 2'-deoxyinosine in DNA when commercial salmon testis DNA was treated with nitrous acid. We also used it to measure levels in various rat tissues of both normal and endotoxin-treated rats, but could not find increased 2'-deoxyinosine formation in tissues even though *NO production was substantially increased.     

High-resolution analysis of cytosine methylation in ancient DNA

Epigenetic changes to gene expression can result in heritable phenotypic characteristics that are not encoded in the DNA itself, but rather by biochemical modifications to the DNA or associated chromatin proteins. Interposed between genes and environment, these epigenetic modifications can be influenced by environmental factors to affect phenotype for multiple generations. This raises the possibility that epigenetic states provide a substrate for natural selection, with the potential to participate in the rapid adaptation of species to changes in environment. Any direct test of this hypothesis would require the ability to measure epigenetic states over evolutionary timescales. Here we describe the first single-base resolution of cytosine methylation patterns in an ancient mammalian genome, by bisulphite allelic sequencing of loci from late Pleistocene Bison priscus remains. Retrotransposons and the differentially methylated regions of imprinted loci displayed methylation patterns identical to those derived from fresh bovine tissue, indicating that methylation patterns are preserved in the ancient DNA. Our findings establish the biochemical stability of methylated cytosines over extensive time frames, and provide the first direct evidence that cytosine methylation patterns are retained in DNA from ancient specimens. The ability to resolve cytosine methylation in ancient DNA provides a powerful means to study the role of epigenetics in evolution.     

Elevated substitution rates estimated from ancient DNA sequences

Ancient DNA sequences are able to offer valuable insights into molecular evolutionary processes, which are not directly accessible via modern DNA. They are particularly suitable for the estimation of substitution rates because their ages provide calibrating information in phylogenetic analyses, circumventing the difficult task of choosing independent calibration points. The substitution rates obtained from such datasets have typically been high, falling between the rates estimated from pedigrees and species phylogenies. Many of these estimates have been made using a Bayesian phylogenetic method that explicitly accommodates heterochronous data. Stimulated by recent criticism of this method, we present a comprehensive simulation study that validates its performance. For datasets of moderate size, it produces accurate estimates of rates, while appearing robust to assumptions about demographic history. We then analyse a large collection of 749 ancient and 727 modern DNA sequences from 19 species of animals, plants and bacteria. Our new estimates confirm that the substitution rates estimated from ancient DNA sequences are elevated above long-term phylogenetic levels.     

Bisulfite-induced cytosine deamination rates in E. coli SSB:DNA complexes

E. coli single-stranded binding protein (SSB) has been examined for its ability to modulate bisulfite-induced cytosine deamination rates in single-stranded DNA (ssDNA). We used a lacZ alpha-complementation reversion assay to detect C-->U rates at a single codon in M13mp2 DNA, whether in free ssDNA or in an SSB:ssDNA complex. When incubated at 37 degrees C, the average bisulfite-induced reversion rate constant was four-fold less in SSB:ssDNA complexes than in ssDNA, at a single codon. Across a 250 base pair target and over 23 scorable C-->U sites, the forward rate constant was 4.9-fold less in SSB:ssDNA complexes than in ssDNA alone. After treatment with N-uracil glycosylase, ssDNA incubated with bisulfite had reversion frequencies at the background rate of ssDNA incubated without bisulfite, indicating that virtually all mutations scored were due to C-->U events. The decrease in cytosine deamination rates occurred both in a single codon and over a 250 bp target, indicating that interactions between SSB and ssDNA reduce bisulfite-catalyzed mutations. The structural role of SSB is well recognized in multiple cellular processes; SSB can also function to minimize bisulfite-induced ssDNA mutations.     

Heat- and alkali-induced deamination of 5-methylcytosine and cytosine residues in DNA

5-Methylcytosine residues in DNA underwent deamination at high temperatures. Furthemore, their rate of deamination at neutral or alkaline pH was greater than that of cytosine residues in DNA. As sources of [¹⁴C]5-methylcytosine-containing DNA, we used bacteriophage XP-12 DNA, in which 5-methylcytosine residues completely replace C residues, and calf thymus DNA experimentally substituted with [¹⁴C]5-methylcytosine residues. Upon incubation at 95°C in a physiological buffer or at 60°C in 1 M NaOH, the respective rates of deamination of 5-methylcytosine residues were about 3- and 1.5-times those of cytosine residues. Under the same conditions, the free 5-methyldeoxycytidine was converted to thymidine more rapidly than deoxycytidine was converted to deoxyuridine. The reactions at physiological pH and elevated temperature suggest that deamination of 5-methylcytosine residues may yield a significant portion of spontaneous mutations in vivo, especially in view of the lack of thymine-specific mismatch repair systems with specificity and efficiency comparable to that of uracil excision repair systems.     

Protection of DNA by alpha/beta-type small,acid-soluble proteins from Bacillus subtilis spores against cytosine deamination

Spores of Bacillus subtilis contain high levels of proteins, termed alpha/beta-type small, acid-soluble proteins (SASP), that protect the spore's DNA against different types of DNA damage. We tested one such protein, SspC, and two of its variants for their ability to protect plasmid DNA against hydrolytic deamination of cytosine to uracil. If unrepaired, such damage to DNA causes C to T mutations. We found that one SspC variant, SspC(Delta 11-D13K), protected DNA against cytosine deamination at two different temperatures (45 and 70 degrees C) and pH values (5.2 and 7.9), reducing the rate of deamination by as much as 10-fold. At 70 degrees C, pH 7.9, the wild-type SspC and its variant, SspC(Delta 11), provided little protection against deamination but were effective in protecting DNA at 45 degrees C, pH 7.9. Parallel studies of the abilities of these proteins to protect DNA against restriction digestion revealed that there was a good correlation between the abilities of the proteins to protect against restriction endonucleases and reductions in cytosine deaminations. These results show that the binding of SspC variants to DNA can prevent attack on DNA bases by water and suggest a new general mechanism by which DNA-binding proteins in cells may be able to protect chromosomes from endogenous and exogenous reactive chemicals by excluding them from the vicinity of DNA.     

Integration of multiple PCR amplifications and DNA mutation analyses by using oligonucleotide microchip

We have developed a method for parallel independent on-chip amplification and the following sequence variation analysis of multiple DNA regions directly using microchip with an array of nanoliter gel pads containing specific sets of tethered primers. The method has three key features. First, DNA to be amplified is enriched at gel pads by its hybridization with immobilized primers. Second, different sets of specific primers are immobilized within various gel pads, and primers are detached within gel pads just before polymerase chain reaction to enhance the amplification. A gel pad may contain an additional permanently immobilized dormant primer that is activated to carry out the allele-specific primer extension reaction to detect mutations. Third, multiple polymerase chain reactions are confined within nanoliter gel pads covered and separated from each other with mineral oil. The method was applied to simultaneously identify several abundant drug-resistant mutations in three genes of Mycobacterium tuberculosis.     

APOBEC3G cytosine deamination hotspots are defined by both sequence context and single-stranded DNA secondary structure

Apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3G (i.e., APOBEC3G or A3G) is an evolutionarily conserved cytosine deaminase that potently restricts human immunodeficiency virus type 1 (HIV-1), retrotransposons and other viruses. A3G has a nucleotide target site specificity for cytosine dinucleotides, though only certain cytosine dinucleotides are ‘hotspots’ for cytosine deamination, and others experience little or no editing by A3G. The factors that define these critical A3G hotspots are not fully understood. To investigate how A3G hotspots are defined, we used an in vitro fluorescence resonance energy transfer-based oligonucleotide assay to probe the site specificity of A3G. Our findings strongly suggest that the target single-stranded DNA (ssDNA) secondary structure as well as the bases directly 3′ and 5′ of the cytosine dinucleotide are critically important A3G recognition. For instance, A3G cannot readily deaminate a cytosine dinucleotide in ssDNA stem structures or in nucleotide base loops composed of three bases. Single-stranded nucleotide loops up to seven bases in length were poor targets for A3G activity unless cytosine residues flanked the cytosine dinucleotide. Furthermore, we observed that A3G favors adenines, cytosines and thymines flanking the cytosine dinucleotide target in unstructured regions of ssDNA. Low cytosine deaminase activity was detected when guanines flanked the cytosine dinucleotide. Taken together, our findings provide the first demonstration that A3G cytosine deamination hotspots are defined by both the sequence context of the cytosine dinucleotide target as well as the ssDNA secondary structure. This knowledge can be used to better trace the origins of mutations to A3G activity, and illuminate its impact on processes such as HIV-1 genetic variation.     

Specific targeting of cytosine methylation to DNA sequences in vivo

Development of methods that will allow exogenous imposition of inheritable gene-specific methylation patterns has potential application in both therapeutics and in basic research. An ongoing approach is the use of targeted DNA methyltransferases, which consist of a fusion between gene-targeted zinc-finger proteins and prokaryotic DNA cytosine methyltransferases. These enzymes however have so far demonstrated significant and unacceptable levels of non-targeted methylation. We now report the development of second-generation targeted methyltransferase enzymes comprising enhanced zinc-finger arrays coupled to methyltransferase mutants that are functionally dominated by their zinc-finger component. Both in vitro plasmid methylation studies and a novel bacterial assay reveal a high degree of target-specific methylation by these enzymes. Furthermore, we demonstrate for the first time transient expression of targeted cytosine methyltransferase in mammalian cells resulting in the specific methylation of a chromosomal locus. Importantly, the resultant methylation pattern is inherited through successive cell divisions.     

Assessing the fidelity of ancient DNA sequences amplified from nuclear genes

To date, the field of ancient DNA has relied almost exclusively on mitochondrial DNA (mtDNA) sequences. However, a number of recent studies have reported the successful recovery of ancient nuclear DNA (nuDNA) sequences, thereby allowing the characterization of genetic loci directly involved in phenotypic traits of extinct taxa. It is well documented that postmortem damage in ancient mtDNA can lead to the generation of artifactual sequences. However, as yet no one has thoroughly investigated the damage spectrum in ancient nuDNA. By comparing clone sequences from 23 fossil specimens, recovered from environments ranging from permafrost to desert, we demonstrate the presence of miscoding lesion damage in both the mtDNA and nuDNA, resulting in insertion of erroneous bases during amplification. Interestingly, no significant differences in the frequency of miscoding lesion damage are recorded between mtDNA and nuDNA despite great differences in cellular copy numbers. For both mtDNA and nuDNA, we find significant positive correlations between total sequence heterogeneity and the rates of type 1 transitions (adenine --> guanine and thymine --> cytosine) and type 2 transitions (cytosine --> thymine and guanine --> adenine), respectively. Type 2 transitions are by far the most dominant and increase relative to those of type 1 with damage load. The results suggest that the deamination of cytosine (and 5-methyl cytosine) to uracil (and thymine) is the main cause of miscoding lesions in both ancient mtDNA and nuDNA sequences. We argue that the problems presented by postmortem damage, as well as problems with contamination from exogenous sources of conserved nuclear genes, allelic variation, and the reliance on single nucleotide polymorphisms, call for great caution in studies relying on ancient nuDNA sequences.     

Chemical carcinogens in non-enzymatic cytosine deamination: 3-isocyanatoacrylonitrile

Uracil has long been known as the main product of nitrosative cytosine deamination in aqueous solution. Recent mechanistic studies of cytosinediazonium ion suggest that the cation formed by its dediazoniation can ring-open to N-protonated (Z,s-cis)-3-isocyanatoacrylonitrile 7. Stereochemical preferences are discussed of the 3-isocyanatoacrylonitriles (Z,s-cis)-10, (E,s-cis)-11, (Z,s-trans)-12, and (E,s-trans)-13. The electronic structures of 7 and 10–13 have been analyzed and a rationale is provided for the thermodynamic preference for (Z,s-cis)-10. It is shown that s-cis/s-trans-interconversion occurs via C−N rotation–inversion paths with barriers below 3 kcal mol⁻¹. The proton affinities of 3-isocyanatoacrylonitrile 10 and water are nearly identical and, thus, 3-isocyanatoacrylonitriles can and should be formed in aqueous media from 7 along with 3-aminoacrylonitriles 9. The results highlight the relevance of the chemistry of 3-isocyanatoacrylonitriles for the understanding of the chemical toxicology of nitrosation of the nucleobase cytosine. Electronic Supplementary Material Supplementary Material is available for this article at Dedicated to professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Unravelling migrations in the steppe: mitochondrial DNA sequences from ancient central Asians

This study helps to clarify the debate on the Western and Eastern genetic influences in Central Asia. Thirty-six skeletal remains from Kazakhstan (Central Asia), excavated from different sites dating between the fifteenth century BC to the fifth century AD, have been analysed for the hypervariable control region (HVR-I) and haplogroup diagnostic single nucleotide polymorphisms (SNPs) of the mitochondrial DNA genome. Standard authentication criteria for ancient DNA studies, including multiple extractions, cloning of PCR products and independent replication, have been followed. The distribution of east and west Eurasian lineages through time in the region is concordant with the available archaeological information: prior to the thirteenth-seventh century BC, all Kazakh samples belong to European lineages; while later an arrival of east Eurasian sequences that coexisted with the previous west Eurasian genetic substratum can be detected. The presence of an ancient genetic substratum of European origin in West Asia may be related to the discovery of ancient mummies with European features in Xinjiang and to the existence of an extinct Indo-European language, Tocharian. This study demonstrates the usefulness of the ancient DNA in unravelling complex patterns of past human migrations so as to help decipher the origin of present-day admixed populations.     

Mitochondrial DNA from various organisms does not contain internally methylated cytosine in -CCGG- sequences.

Mitochondrial DNAs from yeast, Neurospora, rat and calf do not contain internally methylated cytosine in -CCGG- sequences.     

A reanalysis of the ancient mitochondrial DNA sequences recovered from Neandertal bones

Recent reports analyzing mitochondrial DNA sequences from Neandertal bones have claimed that Neandertals and modern humans are different species. The phylogenetic analyses carried out in these articles did not take into account the high substitution rate variation among sites observed in the human mitochondrial D-loop region and also lack an estimation of the parameters of the nucleotide substitution model. The separate phylogenetic position of Neandertals is not supported when these factors are considered. Our analysis shows that Neandertal-Human and Human-Human pairwise distance distributions overlap more than what previous studies suggested. We also show that the most ancient Neandertal HVI region is the most divergent when compared with modern human sequences. However, the opposite would be expected if the sequence had not been modified since the death of the specimen. Such incongruence is discussed in the light of diagenetic modifications in ancient Neandertal DNA sequences.