The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing

Background

We recently showed that enzymes of the TET family convert 5-mC to 5-hydroxymethylcytosine (5-hmC) in DNA. 5-hmC is present at high levels in embryonic stem cells and Purkinje neurons. The methylation status of cytosines is typically assessed by reaction with sodium bisulfite followed by PCR amplification. Reaction with sodium bisulfite promotes cytosine deamination, whereas 5-methylcytosine (5-mC) reacts poorly with bisulfite and is resistant to deamination. Since 5-hmC reacts with bisulfite to yield cytosine 5-methylenesulfonate (CMS), we asked how DNA containing 5-hmC behaves in bisulfite sequencing.

Methodology/Principal Findings

We used synthetic oligonucleotides with different distributions of cytosine as templates for generation of DNAs containing C, 5-mC and 5-hmC. The resulting DNAs were subjected in parallel to bisulfite treatment, followed by exposure to conditions promoting cytosine deamination. The extent of conversion of 5-hmC to CMS was estimated to be 99.7%. Sequencing of PCR products showed that neither 5-mC nor 5-hmC undergo C-to-T transitions after bisulfite treatment, confirming that these two modified cytosine species are indistinguishable by the bisulfite technique. DNA in which CMS constituted a large fraction of all bases (28/201) was much less efficiently amplified than DNA in which those bases were 5-mC or uracil (the latter produced by cytosine deamination). Using a series of primer extension experiments, we traced the inefficient amplification of CMS-containing DNA to stalling of Taq polymerase at sites of CMS modification, especially when two CMS bases were either adjacent to one another or separated by 1–2 nucleotides.

Conclusions

We have confirmed that the widely used bisulfite sequencing technique does not distinguish between 5-mC and 5-hmC. Moreover, we show that CMS, the product of bisulfite conversion of 5-hmC, tends to stall DNA polymerases during PCR, suggesting that densely hydroxymethylated regions of DNA may be underrepresented in quantitative methylation analyses.     

The bisulfite genomic sequencing used in the analysis of epigenetic states, a technique in the emerging environmental genotoxicology research

Methylation at position 5 of cytosine in DNA is an important event in epigenetic changes of cells, the methylation being linked to the control of gene functions. The DNA methylation can be analyzed by bisulfite genomic sequencing, and a large body of data have now been accumulated, based on which causation of diseases, for example cancer, and many other manifestations of cellular activities have been discussed intensively. This article gives an extensive account of the chemical aspects of bisulfite modification of cytosine and 5-methylcytosine in DNA. Various factors that affect the action of bisulfite are discussed, and a recent progress from our laboratory is explained. Conventional procedures for the bisulfite treatment consist incubation of single-stranded DNA with sodium bisulfite under acidic conditions. This treatment converts cytosine into uracil, but 5-methylcytosine remains unchanged. Amplification by polymerase chain reaction (PCR) of the bisulfite-treated DNA followed by sequencing can result in revealing the positions of 5-methylcytosine in the gene. We have discovered that the whole procedure can be significantly speeded up by the use of a highly concentrated bisulfite solution, 10 M ammonium bisulfite. Another recent finding is that urea, which has been often added to the reaction mixture with the purpose of facilitating the bisulfite-mediated deamination of cytosine in DNA, may not work as anticipated: we have observed that urea does not show such promoting actions in our treatments of DNA. A laboratory protocol for quantifying bisulfite, suitably simple for routine practice to ensure valid experiments, is described.     

Identification and resolution of artifacts in bisulfite sequencing

Bisulfite sequencing has become the most widely used application to detect 5-methylcytosine (5-MeC) in DNA, and provides a reliable way of detecting any methylated cytosine at single-molecule resolution in any sequence context. The process of bisulfite treatment exploits the different sensitivity of cytosine and 5-MeC to deamination by bisulfite under acidic conditions, in which cytosine undergoes conversion to uracil while 5-MeC remains unreactive. In this article, we address the more commonly encountered experimental artifacts associated with bisulfite sequencing, and provide methods for the detection and elimination of these artifacts. In particular, we focus on conditions that inhibit complete bisulfite-mediated conversion of cytosines in a target sequence, and demonstrate the necessity of complete protein removal from DNA samples prior to bisulfite treatment. We also include a brief summary of the experimental protocol for bisulfite treatment and tips for designing polymerase chain reaction (PCR) primers to amplify from bisulfite-treated DNA.     

Kinetics of deamination and Cu(II)/H2O2/Ascorbate-induced formation of 5-methylcytosine glycol at CpG sites in duplex DNA

Mutation in p53 tumor suppressor gene is a hallmark of human cancers. Six major mutational hotspots in p53 contain methylated CpG (mCpG) sites, and C →T transition is the most common mutation at these sites. It was hypothesized that the formation of 5-methylcytosine glycol induced by reactive oxygen species, its spontaneous deamination to thymine glycol and the miscoding property of the latter may account, in part, for the ubiquitous C →T mutation at CpG site. Here, we assessed the kinetics of deamination for two diastereomers of 5-methylcytosine glycol in duplex DNA. Our results revealed that the half-lives for the deamination of the (5S,6S) and (5R,6R) diastereomers of 5-methylcytosine glycol in duplex DNA at 37°C were 37.4 ± 1.6 and 27.4 ± 1.0 h, respectively. The deamination rates were only slightly lower than those for the two diastereomers in mononucleosides. Next, we assessed the formation of 5-methyl-2′-deoxycytidine glycol in the form of its deaminated product, namely, thymidine glycol (Tg), in methyl-CpG-bearing duplex DNA treated with Cu(II)/H₂O₂/ascorbate. LC-MS/MS quantification results showed that the yield of Tg is similar as that of 5-(hydroxymethyl)-2′-deoxycytidine. Together, our data support that the formation and deamination of 5-methylcytosine glycol may contribute significantly to the C →T transition mutation at mCpG dinucleotide site.     

High-Speed Conversion of Cytosine to Uracil in Bisulfite Genomic Sequencing Analysis of DNA Methylation

Bisulfite genomic sequencing is a widely used technique foranalyzing cytosine-methylation of DNA. By treating DNA withbisulfite, cytosine residues are deaminated to uracil, whileleaving 5-methylcytosine largely intact. Subsequent PCR andnucleotide sequence analysis permit unequivocal determinationof the methylation status at cytosine residues. A major caveatassociated with the currently practiced procedure is that ittakes 16–20 hr for completion of the conversion of cytosineto uracil. Here we report that a complete deamination of cytosineto uracil can be achieved in shorter periods by using a highlyconcentrated bisulfite solution at an elevated temperature.Time course experiments demonstrated that treating DNA with9 M bisulfite for 20 min at 90°C or 40 min at 70°C allcytosine residues in the DNA were converted to uracil. Underthese conditions, the majority of 5-methylcytosines remainedintact. When a high molecular weight DNA derived from a cellline (containing a number of genes whose methylation statuswas known) was treated with bisulfite under the above conditionsand amplified and sequenced, the results obtained were consistentwith those reported in the literature. Although some degradationof DNA occurred during this process, the amount of treated DNArequired for the amplification was nearly equal to that requiredfor the conventional bisulfite genomic sequencing procedure.The increased speed of DNA methylation analysis with this novelprocedure is expected to advance various aspects of DNA sciences.     

Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging

DNA methylation plays a key role in epigenetic regulation of eukaryotic genomes. Hence the genome-wide distribution of 5-methylcytosine, or the methylome, has been attracting intense attention. In recent years, whole-genome bisulfite sequencing (WGBS) has enabled methylome analysis at single-base resolution. However, WGBS typically requires microgram quantities of DNA as well as global PCR amplification, thereby precluding its application to samples of limited amounts. This is presumably because bisulfite treatment of adaptor-tagged templates, which is inherent to current WGBS methods, leads to substantial DNA fragmentation. To circumvent the bisulfite-induced loss of intact sequencing templates, we conceived an alternative method termed Post-Bisulfite Adaptor Tagging (PBAT) wherein bisulfite treatment precedes adaptor tagging by two rounds of random primer extension. The PBAT method can generate a substantial number of unamplified reads from as little as subnanogram quantities of DNA. It requires only 100 ng of DNA for amplification-free WGBS of mammalian genomes. Thus, the PBAT method will enable various novel applications that would not otherwise be possible, thereby contributing to the rapidly growing field of epigenomics.     

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.     

Directional telomeric silencing and lack of canonical B1 elements in two silencer Autonomously Replicating Sequences in S. cerevisiae

Background

Detection of cell-free methylated DNA in plasma is a promising tool for tumour diagnosis and monitoring. Due to the very low amounts of cell-free DNA in plasma, analytical sensitivity is of utmost importance. The vast majority of currently available methods for analysing DNA methylation are based on bisulfite-mediated deamination of cytosine. Cytosine is rapidly converted to uracil during bisulfite treatment, whereas 5-methylcytosine is only slowly converted. Hence, bisulfite treatment converts an epigenetic modification into a difference in sequence, amenable to analysis either by sequencing or PCR based methods. However, the recovery of bisulfite-converted DNA is very poor.

Results

Here we introduce an alternative method for the crucial steps of bisulfite treatment with high recovery. The method is based on an accelerated deamination step and alkaline desulfonation in combination with magnetic silica purification of DNA, allowing preparation of deaminated DNA from patient samples in less than 2 hours.

Conclusions

The method presented here allows low levels of DNA to be easily and reliably analysed, a prerequisite for the clinical usefulness of cell-free methylated DNA detection in plasma.     

A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements

We report a method for studying global DNA methylation based on using bisulfite treatment of DNA and simultaneous PCR of multiple DNA repetitive elements, such as Alu elements and long interspersed nucleotide elements (LINE). The PCR product, which represents a pool of approximately 15000 genomic loci, could be used for direct sequencing, selective restriction digestion or pyrosequencing, in order to quantitate DNA methylation. By restriction digestion or pyrosequencing, the assay was reproducible with a standard deviation of only 2% between assays. Using this method we found that almost two-thirds of the CpG methylation sites in Alu elements are mutated, but of the remaining methylation target sites, 87% were methylated. Due to the heavy methylation of repetitive elements, this assay was especially useful in detecting decreases in DNA methylation, and this assay was validated by examining cell lines treated with the methylation inhibitor 5-aza-2′deoxycytidine (DAC), where we found a 1–16% decrease in Alu element and 18–60% LINE methylation within 3 days of treatment. This method can be used as a surrogate marker of genome-wide methylation changes. In addition, it is less labor intensive and requires less DNA than previous methods of assessing global DNA methylation.     

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k_2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.     

Thymine DNA glycosylase exhibits negligible affinity for nucleobases that it removes from DNA

Thymine DNA Glycosylase (TDG) performs essential functions in maintaining genetic integrity and epigenetic regulation. Initiating base excision repair, TDG removes thymine from mutagenic G·T mispairs caused by 5-methylcytosine (mC) deamination and other lesions including uracil (U) and 5-hydroxymethyluracil (hmU). In DNA demethylation, TDG excises 5-formylcytosine (fC) and 5-carboxylcytosine (caC), which are generated from mC by Tet (ten–eleven translocation) enzymes. Using improved crystallization conditions, we solved high-resolution (up to 1.45 Å) structures of TDG enzyme–product complexes generated from substrates including G·U, G·T, G·hmU, G·fC and G·caC. The structures reveal many new features, including key water-mediated enzyme–substrate interactions. Together with nuclear magnetic resonance experiments, the structures demonstrate that TDG releases the excised base from its tight product complex with abasic DNA, contrary to previous reports. Moreover, DNA-free TDG exhibits no significant binding to free nucleobases (U, T, hmU), indicating a K_d >> 10 mM. The structures reveal a solvent-filled channel to the active site, which might facilitate dissociation of the excised base and enable caC excision, which involves solvent-mediated acid catalysis. Dissociation of the excised base allows TDG to bind the beta rather than the alpha anomer of the abasic sugar, which might stabilize the enzyme–product complex.     

Cytosine-to-Uracil Deamination by SssI DNA Methyltransferase

The prokaryotic DNA(cytosine-5)methyltransferase M.SssI shares the specificity of eukaryotic DNA methyltransferases (CG) and is an important model and experimental tool in the study of eukaryotic DNA methylation. Previously, M.SssI was shown to be able to catalyze deamination of the target cytosine to uracil if the methyl donor S-adenosyl-methionine (SAM) was missing from the reaction. To test whether this side-activity of the enzyme can be used to distinguish between unmethylated and C5-methylated cytosines in CG dinucleotides, we re-investigated, using a sensitive genetic reversion assay, the cytosine deaminase activity of M.SssI. Confirming previous results we showed that M.SssI can deaminate cytosine to uracil in a slow reaction in the absence of SAM and that the rate of this reaction can be increased by the SAM analogue 5’-amino-5’-deoxyadenosine. We could not detect M.SssI-catalyzed deamination of C5-methylcytosine (^m5C). We found conditions where the rate of M.SssI mediated C-to-U deamination was at least 100-fold higher than the rate of ^m5C-to-T conversion. Although this difference in reactivities suggests that the enzyme could be used to identify C5-methylated cytosines in the epigenetically important CG dinucleotides, the rate of M.SssI mediated cytosine deamination is too low to become an enzymatic alternative to the bisulfite reaction. Amino acid replacements in the presumed SAM binding pocket of M.SssI (F17S and G19D) resulted in greatly reduced methyltransferase activity. The G19D variant showed cytosine deaminase activity in E. coli, at physiological SAM concentrations. Interestingly, the C-to-U deaminase activity was also detectable in an E. coli ung ⁺ host proficient in uracil excision repair.     

An AP Endonuclease Functions in Active DNA Demethylation and Gene Imprinting in Arabidopsis

Active DNA demethylation in plants occurs through base excision repair, beginning with removal of methylated cytosine by the ROS1/DME subfamily of 5-methylcytosine DNA glycosylases. Active DNA demethylation in animals requires the DNA glycosylase TDG or MBD4, which functions after oxidation or deamination of 5-methylcytosine, respectively. However, little is known about the steps following DNA glycosylase action in the active DNA demethylation pathways in plants and animals. We show here that the Arabidopsis APE1L protein has apurinic/apyrimidinic endonuclease activities and functions downstream of ROS1 and DME. APE1L and ROS1 interact in vitro and co-localize in vivo. Whole genome bisulfite sequencing of ape1l mutant plants revealed widespread alterations in DNA methylation. We show that the ape1l/zdp double mutant displays embryonic lethality. Notably, the ape1l^+/−zdp^−/− mutant shows a maternal-effect lethality phenotype. APE1L and the DNA phosphatase ZDP are required for FWA and MEA gene imprinting in the endosperm and are important for seed development. Thus, APE1L is a new component of the active DNA demethylation pathway and, together with ZDP, regulates gene imprinting in Arabidopsis.     

Preferential CEBP binding to T:G mismatches and increased C-to-T human somatic mutations

DNA cytosine methylation in mammals modulates gene expression and chromatin accessibility. It also impacts mutation rates, via spontaneous oxidative deamination of 5-methylcytosine (5mC) to thymine. In most cases the resulting T:G mismatches are repaired, following T excision by one of the thymine DNA glycosylases, TDG or MBD4. We found that C-to-T mutations are enriched in the binding sites of CCAAT/enhancer binding proteins (CEBP). Within a CEBP site, the presence of a T:G mismatch increased CEBPβ binding affinity by a factor of >60 relative to the normal C:G base pair. This enhanced binding to a mismatch inhibits its repair by both TDG and MBD4 in vitro. Furthermore, repair of the deamination product of unmethylated cytosine, which yields a U:G DNA mismatch that is normally repaired via uracil DNA glycosylase, is also inhibited by CEBPβ binding. Passage of a replication fork over either a T:G or U:G mismatch, before repair can occur, results in a C-to-T mutation in one of the daughter duplexes. Our study thus provides a plausible mechanism for accumulation of C-to-T human somatic mutations.     

The rate of hydrolytic deamination of 5-methylcytosine in double-stranded DNA.

The modified base, 5-methylcytosine, constitutes approximately 1% of human DNA, but sites containing 5-methylcytosine account for at least 30% of all germline and somatic point mutations. A genetic assay with a sensitivity of 1 in 10(7), based on reversion to neomycin resistance of a mutant pSV2-neo plasmid, was utilized to determine and compare the deamination rates of 5-methylcytosine and cytosine in double-stranded DNA for the first time. The rate constants for spontaneous hydrolytic deamination of 5-methylcytosine and cytosine in double-stranded DNA at 37 degrees C were 5.8 x 10(-13) s-1 and 2.6 x 10(-13) s-1, respectively. These rates are more than sufficient to explain the observed frequency of mutation at sites containing 5-methylcytosine and emphasize the importance of hydrolytic deamination as a major source of human mutations.     

Effect of Temperature on Structure and Function of the Methanogenic Archaeal Community in an Anoxic Rice Field Soil

Soil temperatures in Italian rice fields typically range between about 15 and 30°C. A change in the incubation temperature of anoxic methanogenic soil slurry from 30°C to 15°C typically resulted in a decrease in the CH₄ production rate, a decrease in the steady-state H₂ partial pressure, and a transient accumulation of acetate. Previous experiments have shown that these changes were due to an alteration of the carbon and electron flow in the methanogenic degradation pathway of organic matter caused by the temperature shift (K. J. Chin and R. Conrad, FEMS Microbiol. Ecol. 18:85–102, 1995). To investigate how temperature affects the structure of the methanogenic archaeal community, total DNA was extracted from soil slurries incubated at 30 and 15°C. The archaeal small-subunit (SSU) rRNA-encoding genes (rDNA) of these environmental DNA samples were amplified by PCR with an archaeal-specific primer system and used for the generation of clone libraries. Representative rDNA clones (n = 90) were characterized by terminal restriction fragment length polymorphism (T-RFLP) and sequence analysis. T-RFLP analysis produced for the clones terminally labeled fragments with a characteristic length of mostly 185, 284, or 392 bp. Sequence analysis allowed determination of the phylogenetic affiliation of the individual clones with their characteristic T-RFLP fragment lengths and showed that the archaeal community of the anoxic rice soil slurry was dominated by members of the families Methanosarcinaceae (185 bp) and Methanosaetaceae (284 bp), the kingdom Crenarchaeota (185 or 284 bp), and a novel, deeply branching lineage of the (probably methanogenic) kingdom Euryarchaeota (392 bp) that has recently been detected on rice roots (R. Großkopf, S. Stubner, and W. Liesack, Appl. Environ. Microbiol. 64:4983–4989, 1998). The structure of the archaeal community changed when the temperature was shifted from 30°C to 15°C. Before the temperature shift, the clones (n = 30) retrieved from the community were dominated by Crenarchaeota (70%), “novel Euryarchaeota” (23%), and Methanosarcinacaeae (7%). Further incubation at 30°C (n = 30 clones) resulted in a relative increase in members of the Methanosarcinaceae (77%), whereas further incubation at 15°C (n = 30 clones) resulted in a much more diverse community consisting of 33% Methanosarcinaceae, 23% Crenarchaeota, 20% Methanosaetaceae, and 17% novel Euryarchaeota. The appearance of Methanosaetaceae at 15°C was conspicuous. These results demonstrate that the structure of the archaeal community in anoxic rice field soil changed with time and incubation temperature.     

Using structural motif templates to identify proteins with DNA binding function

This work describes a method for predicting DNA binding function from structure using 3-dimensional templates. Proteins that bind DNA using small contiguous helix–turn–helix (HTH) motifs comprise a significant number of all DNA-binding proteins. A structural template library of seven HTH motifs has been created from non-homologous DNA-binding proteins in the Protein Data Bank. The templates were used to scan complete protein structures using an algorithm that calculated the root mean squared deviation (rmsd) for the optimal superposition of each template on each structure, based on C_α backbone coordinates. Distributions of rmsd values for known HTH-containing proteins (true hits) and non-HTH proteins (false hits) were calculated. A threshold value of 1.6 Å rmsd was selected that gave a true hit rate of 88.4% and a false positive rate of 0.7%. The false positive rate was further reduced to 0.5% by introducing an accessible surface area threshold value of 990 Ų per HTH motif. The template library and the validated thresholds were used to make predictions for target proteins from a structural genomics project.     

Chemical synthesis and translesion replication of a cis-syn cyclobutane thymine-uracil dimer

The cytosine base in DNA undergoes hydrolytic deamination at a considerable rate when UV radiation induces formation of a cyclobutane pyrimidine dimer (CPD) with an adjacent pyrimidine base. We have synthesized a phosphoramidite building block of a cis–syn cyclobutane thymine–uracil dimer (T[]U), which is the deaminated form of the CPD at a TC site, and incorporated it into oligodeoxyribonucleotides. The previously reported method for synthesis of the thymine dimer (T[]T) was applied, using partially protected thymidylyl-(3′–5′)-2′-deoxyuridine as the starting material, and after triplet- sensitized irradiation, the configuration of the base moiety in the major product was determined by NMR spectroscopy. Presence of the cis–syn cyclobutane dimer in the obtained oligonucleotides was confirmed by UV photoreversal and reaction with T4 endonuclease V. Using a 30mer containing T[]U, translesion synthesis by human DNA polymerase η was analyzed. There was no difference in the results between the templates containing T[]T and T[]U and pol η bypassed both lesions with the same efficiency, incorporating two adenylates. This enzyme showed fidelity to base pair formation, but this replication causes a C→T transition because the original sequence is TC.     

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.     

Determination of 5-methylcytosine by acid hydrolysis of DNA with hydrofluoric acid

Quantitation of 5-methylcytosine in DNA after acid hydrolysis has been inaccurate because deamination of cytosine and 5-methylcytosine occurs during the hydrolysis procedure. There is little information in the literature regarding the use of hydrofluoric acid (HF) for DNA hydrolysis and we have therefore undertaken a systematic study of this process. The deoxyribonucleotides of cytosine and 5-methylcytosine were shown not to undergo detectable levels of deamination during prolonged periods (up to 24 h) at 80 degrees C in 48% HF. Kinetic studies show that the release of purine and pyrimidine bases was complete by 4 h under these conditions. Analysis of the 5-methylcytosine content of DNA from various tissues gave levels that were very close to the values reported in the literature. This method is ideally suited for the determination of the overall cytosine methylation levels in DNA.