Genome-wide mapping of DNA strand breaks

Determination of cellular DNA damage has so far been limited to global assessment of genome integrity whereas nucleotide-level mapping has been restricted to specific loci by the use of specific primers. Therefore, only limited DNA sequences can be studied and novel regions of genomic instability can hardly be discovered. Using a well-characterized yeast model, we describe a straightforward strategy to map genome-wide DNA strand breaks without compromising nucleotide-level resolution. This technique, termed "damaged DNA immunoprecipitation" (dDIP), uses immunoprecipitation and the terminal deoxynucleotidyl transferase-mediated dUTP-biotin end-labeling (TUNEL) to capture DNA at break sites. When used in combination with microarray or next-generation sequencing technologies, dDIP will allow researchers to map genome-wide DNA strand breaks as well as other types of DNA damage and to establish a clear profiling of altered genes and/or intergenic sequences in various experimental conditions. This mapping technique could find several applications for instance in the study of aging, genotoxic drug screening, cancer, meiosis, radiation and oxidative DNA damage.     

Spatial organization of the mouse genome and its role in recurrent chromosomal translocations

The extent to which the three-dimensional organization of the genome contributes to chromosomal translocations is an important question in cancer genomics. We generated a high-resolution Hi-C spatial organization map of the G1-arrested mouse pro-B cell genome and used high-throughput genome-wide translocation sequencing to map translocations from target DNA double-strand breaks (DSBs) within it. RAG endonuclease-cleaved antigen-receptor loci are dominant translocation partners for target DSBs regardless of genomic position, reflecting high-frequency DSBs at these loci and their colocalization in a fraction of cells. To directly assess spatial proximity contributions, we normalized genomic DSBs via ionizing radiation. Under these conditions, translocations were highly enriched in cis along single chromosomes containing target DSBs and within other chromosomes and subchromosomal domains in a manner directly related to pre-existing spatial proximity. By combining two high-throughput genomic methods in a genetically tractable system, we provide a new lens for viewing cancer genomes.     

Asp-Ala-His-Lys (DAHK) inhibits copper-induced oxidative DNA double strand breaks and telomere shortening

Both DNA and the telomeric sequence are susceptible to copper-mediated reactive oxygen species (ROS) damage, particularly damage attributed to hydroxyl radicals. In this study, ROS-induced DNA double strand breaks and telomere shortening were produced by exposure to copper and ascorbic acid. Asp-Ala-His-Lys (DAHK), a specific copper chelating tetrapeptide d-analog of the N-terminus of human albumin, attenuated DNA strand breaks in a dose dependent manner. d-DAHK, at a ratio of 4:1 (d-DAHKCu), provided complete protection of isolated DNA from double strand breaks and, at a ratio of 2:1 (d-DAHKCu), completely protected DNA in Raji cells exposed to copper/ascorbate. Southern blots of DNA treated with copper/ascorbate showed severe depletion and shortening of telomeres and Raji cell treated samples showed some conservation of telomere sequences. d-DAHK provided complete telomere length protection at a ratio of 2:1 (d-DAHKCu). The human albumin N-terminus analog, d-DAHK, protects DNA and telomeres against copper-mediated ROS damage and may be a useful therapeutic adjunct in ROS disease processes.     

RADAR-seq: A RAre DAmage and Repair sequencing method for detecting DNA damage on a genome-wide scale

RAre DAmage and Repair sequencing (RADAR-seq) is a highly adaptable sequencing method that enables the identification and detection of rare DNA damage events for a wide variety of DNA lesions at single-molecule resolution on a genome-wide scale. In RADAR-seq, DNA lesions are replaced with a patch of modified bases that can be directly detected by Pacific Biosciences Single Molecule Real-Time (SMRT) sequencing. RADAR-seq enables dynamic detection over a wide range of DNA damage frequencies, including low physiological levels. Furthermore, without the need for DNA amplification and enrichment steps, RADAR-seq provides sequencing coverage of damaged and undamaged DNA across an entire genome. Here, we use RADAR-seq to measure the frequency and map the location of ribonucleotides in wild-type and RNaseH2-deficient E. coli and Thermococcus kodakarensis strains. Additionally, by tracking ribonucleotides incorporated during in vivo lagging strand DNA synthesis, we determined the replication initiation point in E. coli, and its relation to the origin of replication (oriC). RADAR-seq was also used to map cyclobutane pyrimidine dimers (CPDs) in Escherichia coli (E. coli) genomic DNA exposed to UV-radiation. On a broader scale, RADAR-seq can be applied to understand formation and repair of DNA damage, the correlation between DNA damage and disease initiation and progression, and complex biological pathways, including DNA replication.     

Measurement of cellular repair activities for oxidative DNA damage

Two related assays capable of determining cell extract repair activities for different oxidative lesions in DNA are described. Both assays measure the incorporation of radiolabeled nucleotides during repair of an oxidatively damaged template in a cell-free system. The assays differ in the type of oxidative damage present in the DNA. In one, singlet oxygen is used to generate predominantly 8-oxo-2'-deoxyguanosine lesions. In the other, hydroxyl radicals are used to generate a broad spectrum of damage including oxidized bases and strand breaks. Assay conditions were adjusted to ensure that radiolabel incorporation was directly proportional to cell extract repair activity. These assays represent sensitive tools for investigating the regulation of repair systems for oxidative DNA damage.     

Genome-wide analysis of DNA replication and DNA double-strand breaks using TrAEL-seq

Faithful replication of the entire genome requires replication forks to copy large contiguous tracts of DNA, and sites of persistent replication fork stalling present a major threat to genome stability. Understanding the distribution of sites at which replication forks stall, and the ensuing fork processing events, requires genome-wide methods that profile replication fork position and the formation of recombinogenic DNA ends. Here, we describe Transferase-Activated End Ligation sequencing (TrAEL-seq), a method that captures single-stranded DNA 3′ ends genome-wide and with base pair resolution. TrAEL-seq labels both DNA breaks and replication forks, providing genome-wide maps of replication fork progression and fork stalling sites in yeast and mammalian cells. Replication maps are similar to those obtained by Okazaki fragment sequencing; however, TrAEL-seq is performed on asynchronous populations of wild-type cells without incorporation of labels, cell sorting, or biochemical purification of replication intermediates, rendering TrAEL-seq far simpler and more widely applicable than existing replication fork direction profiling methods. The specificity of TrAEL-seq for DNA 3′ ends also allows accurate detection of double-strand break sites after the initiation of DNA end resection, which we demonstrate by genome-wide mapping of meiotic double-strand break hotspots in a dmc1Δ mutant that is competent for end resection but not strand invasion. Overall, TrAEL-seq provides a flexible and robust methodology with high sensitivity and resolution for studying DNA replication and repair, which will be of significant use in determining mechanisms of genome instability.

TrAEL-seq provides genome-wide base pair resolution maps of exposed DNA 3’ ends; this reveals replication fork stalling and normal replication profiles in asynchronous, unlabelled wildtype cell populations, along with the sites of resected DNA breaks.     

HiCORE: Hi-C Analysis for Identification of Core Chromatin Looping Regions with Higher Resolution

Genome-wide chromosome conformation capture (3C)-based high-throughput sequencing (Hi-C) has enabled identification of genome-wide chromatin loops. Because the Hi-C map with restriction fragment resolution is intrinsically associated with sparsity and stochastic noise, Hi-C data are usually binned at particular intervals; however, the binning method has limited reliability, especially at high resolution. Here, we describe a new method called HiCORE, which provides simple pipelines and algorithms to overcome the limitations of single-layered binning and predict core chromatin regions with three-dimensional physical interactions. In this approach, multiple layers of binning with slightly shifted genome coverage are generated, and interacting bins at each layer are integrated to infer narrower regions of chromatin interactions. HiCORE predicts chromatin looping regions with higher resolution, both in human and Arabidopsis genomes, and contributes to the identification of the precise positions of potential genomic elements in an unbiased manner.     

Cigarette smoke-induced DNA damage in cultured human lung cells: role of hydroxyl radicals and endonuclease activation.

Cigarette smoke can cause DNA single strand breaks in cultured human lung cells (T. Nakayama et al., Nature, 314 (1985) 462-464) but the mechanisms behind this DNA damage have not been clearly elucidated. In the present study we have investigated the possibility that one of the major constituents in cigarette smoke, hydroquinone, may be important for mediating smoke-induced DNA damage in the human epithelial lung cell line, A 549, and the mechanisms behind this damage. Cells were exposed to cigarette smoke, hydrogen peroxide, or hydroquinone, in the absence and presence of different inhibitors, and the resulting DNA damage was assessed either as DNA single strand break formation or formation of the oxidative DNA adduct, 8-hydroxydeoxyguanosine. It was found that (i) exposure to cigarette smoke, hydrogen peroxide or hydroquinone causes a rapid decrease in the intracellular thiol level and a considerable DNA single strand break formation, (ii) the formation of DNA single strand breaks in cells exposed to cigarette smoke is inhibited by catalase, dimethylthiourea, and o-phenantroline, suggesting that hydroxyl radicals generated from iron-catalyzed hydrogen peroxide dissociation are involved in the DNA damage, (iii) hydroquinone causes considerable DNA strand break formation that is blocked by aurintricarboxylic acid, an inhibitor of endonuclease activation, and by BAPTA, an intracellular calcium chelator, (iv) addition of hydroquinone to a smoke condensate greatly enhances its ability to cause DNA single strand breaks, and (v) smoke, but not hydroquinone, causes formation of 8-hydroxydeoxyguanosine, a DNA damage product induced by the action of hydroxyl radicals on the DNA base, deoxyguanosine. These findings suggest that the ability of cigarette smoke to cause DNA single strand breaks in cultured lung cells is due to mechanisms involving hydroxyl radical attack on DNA and endonuclease activation. They also suggest that hydroquinone is an important contributor to the DNA damaging effect of cigarette smoke on human lung cells.     

Arsenite induces oxidative DNA adducts and DNA-protein cross-links in mammalian cells.

Arsenic is generally recognized as a nonmutagenic carcinogen because sodium arsenite induces DNA damage only at very high concentrations. In this study we demonstrate that arsenite concentrations above 0.25 microM induce DNA strand breaks in both human leukemia cells and Chinese hamster ovary cells. Therefore, DNA damage may be involved in arsenic-induced carcinogenesis. Formamidopyrimidine-DNA glycosylase and proteinase K greatly increased DNA strand breaks in arsenite-treated cells, providing evidence that a large portion of arsenite-induced DNA strand breaks come from excision of oxidative DNA adducts and DNA-protein cross-links. Because DNA strand breaks appear only temporarily during excision repair, the level of detectable DNA strand breaks will be low at any given time point. For this reason many previous studies have only detected low levels of DNA strand breaks. We also show that catalase, and inhibitors of calcium, nitric oxide synthase, superoxide dismutase, and myeloperoxidase, could modulate arsenite-induced DNA damage. We conclude that arsenite induces DNA adducts through calcium-mediated production of peroxynitrite, hypochlorous acid, and hydroxyl radicals.     

Quantitative measurement of hydroxyl radical induced DNA double-strand breaks and the effect of N-acetyl-L-cysteine

Reactive oxygen species, such as hydroxyl or superoxide radicals, can be generated by exogenous agents as well as from normal cellular metabolism. Those radicals are known to induce various lesions in DNA, including strand breaks and base modifications. These lesions have been implicated in a variety of diseases such as cancer, arteriosclerosis, arthritis, neurodegenerative disorders and others. To assess these oxidative DNA damages and to evaluate the effects of the antioxidant N-acetyl-L-cysteine (NAC), atomic force microscopy (AFM) was used to image DNA molecules exposed to hydroxyl radicals generated via Fenton chemistry. AFM images showed that the circular DNA molecules became linear after incubation with hydroxyl radicals, indicating the development of double-strand breaks. The occurrence of the double-strand breaks was found to depend on the concentration of the hydroxyl radicals and the duration of the reaction. Under the conditions of the experiments, NAC was found to exacerbate the free radical-induced DNA damage.     

Low- and high-resolution mapping of DNA damage at specific sites

Measurement of DNA damage and repair at the nucleotide level in intact cells has provided compelling evidence for the molecular details of these events as they occur in intact organisms. Furthermore, these measurements give the most accurate picture of the rates of repair in different structural domains of DNA in chromatin. In this report, we describe two methods currently used in our laboratories to map DNA lesions at (or near) nucleotide resolution in yeast cells. The low-resolution method couples damage-specific strand breaks in DNA with indirect end-labeling to measure DNA lesions over a span of 1.5 to 2 kb of DNA sequence. The resolution of this method is limited by the resolution of DNA length measurements on alkaline agarose gels (about +/-20 bp on average). The high-resolution method uses streptavidin magnetic beads and special biotinylated oligonucleotides to facilitate end-labeling of DNA fragments specifically cleaved at damage sites. The latter method maps DNA damage sites at nucleotide resolution over a shorter distance (<500 bp), and is constrained to the length of DNA resolvable on DNA sequencing gels. These methods are used in tandem for answering questions regarding DNA damage and repair in different chromatin domains and states of gene expression.     

Human Xip1 (C2orf13) is a novel regulator of cellular responses to DNA strand breaks

DNA strand breaks arise continuously as the result of intracellular metabolism and in response to a multitude of genotoxic agents. To overcome such challenges to genomic stability, cells have evolved genome surveillance pathways that detect and repair damaged DNA in a coordinated fashion. Here we identify the previously uncharacterized human protein Xip1 (C2orf13) as a novel component of the checkpoint response to DNA strand breaks. Green fluorescent protein-tagged Xip1 was rapidly recruited to sites of DNA breaks, and this accumulation was dependent on a novel type of zinc finger motif located in the C terminus of Xip1. The initial recruitment kinetics of Xip1 closely paralleled that of XRCC1, a central organizer of single strand break (SSB) repair, and its accumulation was both delayed and sustained when the detection of SSBs was abrogated by inhibition of PARP-1. Xip1 and XRCC1 stably interacted through recognition of CK2 phosphorylation sites in XRCC1 by the Forkhead-associated (FHA) domain of Xip1, and XRCC1 was required to maintain steady-state levels of Xip1. Moreover, Xip1 was phosphorylated on Ser-116 by ataxia telangiectasia-mutated in response to ionizing radiation, further underscoring the potential importance of Xip1 in the DNA damage response. Finally, depletion of Xip1 significantly decreased the clonogenic survival of cells exposed to DNA SSB- or double strand break-inducing agents. Collectively, these findings implicate Xip1 as a new regulator of genome maintenance pathways, which may function to organize DNA strand break repair complexes at sites of DNA damage.     

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

Clusters of S1 nuclease-hypersensitive sites induced in vivo by DNA damage.

DNA end-labeling procedures were used to analyze both the frequency and distribution of DNA strand breaks in mammalian cells exposed or not to different types of DNA-damaging agents. The 3' ends were labeled by T4 DNA polymerase-catalyzed nucleotide exchange carried out in the absence or presence of Escherichia coli endonuclease IV to cleave abasic sites and remove 3' blocking groups. Using this sensitive assay, we show that DNA isolated from human cells or mouse tissues contains variable basal levels of DNA strand interruptions which are associated with normal bioprocesses, including DNA replication and repair. On the other hand, distinct dose-dependent patterns of DNA damage were assessed quantitatively in cultured human cells exposed briefly to menadione, methylmethane sulfonate, topoisomerase II inhibitors, or gamma rays. In vivo induction of single-strand breaks and abasic sites by methylmethane sulfonate was also measured in several mouse tissues. The genomic distribution of these lesions was investigated by DNA cleavage with the single-strand-specific S1 nuclease. Strikingly similar cleavage patterns were obtained with all DNA-damaging agents tested, indicating that the majority of S1-hypersensitive sites detected were not randomly distributed over the genome but apparently were clustered in damage-sensitive regions. The parallel disappearance of 3' ends and loss of S1-hypersensitive sites during post-gamma-irradiation repair periods indicates that these sites were rapidly repaired single-strand breaks or gaps (2- to 3-min half-life). Comparison of S1 cleavage patterns obtained with gamma-irradiated DNA and gamma-irradiated cells shows that chromatin structure was the primary determinant of the distribution of the DNA damage detected.     

Immuno-capture of UVDE generated 3’-OH ends at UV photoproducts

A strategy amenable to the genome-wide study of DNA damage and repair kinetics is described. The ultraviolet damage endonuclease (UVDE) generates 3’-OH ends at the two major UV induced DNA lesions, cyclobutane pyrimidine dimers (CPDs) and 6,4 pyrimidine-pyrimidone dimers (6,4 PPs), allowing for their capture after biotin end-labeling. qPCR amplification of biotinylated DNA enables parallel measuring of DNA damage in several loci, which can then be combined with high-throughput screening of cell survival to test genotoxic reagents. Alternatively, a library of captured sequences could be generated for a genome wide study of damage sites and large-scale assessment of repair kinetics in different regions of the genome, using next-generation sequencing. The assay is suitable to study any DNA lesion that can be converted into 3’-OH by UVDE, or other enzymes. Toward these goals, we compared UVDE with the classical T4 endonuclease V (T4V) assay. We showed that there is a linear correlation between UV dose, 3’-OH formation and capture by immunoprecipitation, together with its potential application for in vivo studies.     

Genotoxicity of singlet oxygen

Singlet oxygen, ¹O₂(¹Δ_g), fulfills essential prerequisites for a genotoxic substance, like hydroxyl radicals and other oxygen radicals: it can react efficiently with DNA and it can be generated inside cells, e.g. by photosensitization and enzymatic oxidation. As might be anticipated from the non-radical character of singlet oxygen, the pattern of DNA modifications it produces is very different from that caused by hydroxyl radicals. While hydroxyl radicals produce DNA strand breaks and sites of base loss (AP sites) in high yield and react with all four bases of DNA, singlet oxygen generates predominantly modified guanine residues and few strand breaks and AP sites. There is now convincing evidence that a major product of base modification caused by singlet oxygen is 8-hydroxyguanine (7,8-dihydro-8-oxoguanine). Indeed, the recently reported miscoding properties of 8-hydroxyguanine can explain the predominant type of mutations observed when DNA modified by singlet oxygen is replicated in cells. There are also strong indications that singlet oxygen generated by photosensitization can act as an ultimate DNA modifying species inside cells. However, indirect genotoxic mechanisms involving other reactive oxygen species produced from singlet oxygen are also possible and appear to predominate in some cases. The cellular defense system against oxidants consists of effective singlet oxygen scavengers such as carotenoids. The observation that carotenoids can inhibit neoplastic cell transformation when administered not only together with but also after the application of chemical or physical carcinogens might indicate a role of singlet oxygen in tumor promotion that could be independent of the direct or indirect DNA damaging properties.     

Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling

Genome-wide mapping of 5-methylcytosine is of broad interest to many fields of biology and medicine. A variety of methods have been developed, and several have recently been advanced to genome-wide scale using arrays and next-generation sequencing approaches. We have previously reported reduced representation bisulfite sequencing (RRBS), a bisulfite-based protocol that enriches CG-rich parts of the genome, thereby reducing the amount of sequencing required while capturing the majority of promoters and other relevant genomic regions. The approach provides single-nucleotide resolution, is highly sensitive and provides quantitative DNA methylation measurements. This protocol should enable any standard molecular biology laboratory to generate RRBS libraries of high quality. Briefly, purified genomic DNA is digested by the methylation-insensitive restriction enzyme MspI to generate short fragments that contain CpG dinucleotides at the ends. After end-repair, A-tailing and ligation to methylated Illumina adapters, the CpG-rich DNA fragments (40-220 bp) are size selected, subjected to bisulfite conversion, PCR amplified and end sequenced on an Illumina Genome Analyzer. Note that alignment and analysis of RRBS sequencing reads are not covered in this protocol. The extremely low input requirements (10-300 ng), the applicability of the protocol to formalin-fixed and paraffin-embedded samples, and the technique's single-nucleotide resolution extends RRBS to a wide range of biological and clinical samples and research applications. The entire process of RRBS library construction takes ～9 d.