We present a genome-wide method to map DNA double-strand breaks (DSBs) at nucleotide resolution by direct in situ breaks labeling, enrichment on streptavidin, and next-generation sequencing (BLESS). We comprehensively validated and tested BLESS using different human and mouse cells, DSBs-inducing agents, and sequencing platforms. BLESS was able to detect telomere ends, Sce endonuclease-induced DSBs, and complex genome-wide DSBs landscapes. As a proof of principle, we characterized the genomic landscape of sensitivity to replication stress in human cells, and identified over two thousand non-uniformly distributed aphidicolin-sensitive regions (ASRs) overrepresented in genes and enriched in satellite repeats. ASRs were also enriched in regions rearranged in human cancers, with many cancer-associated genes exhibiting high sensitivity to replication stress. Our method is suitable for genome-wide mapping of DSBs in various cells and experimental conditions with a specificity and resolution unachievable by current techniques.

SummaryDouble-strand breaks (DSBs) are extremely detrimental DNA lesions that can lead to cancer-driving mutations and translocations. Non-homologous end joining (NHEJ) and homologous recombination (HR) represent the two main repair pathways operating in the context of chromatin to ensure genome stability. Despite extensive efforts, our knowledge of DSB-induced chromatin still remains fragmented. Here, we describe the distribution of 20 chromatin features at multiple DSBs spread throughout the human genome using ChIP-seq. We provide the most comprehensive picture of the chromatin landscape set up at DSBs and identify NHEJ- and HR-specific chromatin events. This study revealed the existence of a DSB-induced monoubiquitination-to-acetylation switch on histone H2B lysine 120, likely mediated by the SAGA complex, as well as higher-order signaling at HR-repaired DSBs whereby histone H1 is evicted while ubiquitin and 53BP1 accumulate over the entire γH2AX domains.

The ability of DNA Double Strand Breaks (DSBs) to cluster in mammalian cells has been subjected to intense debate over the past few years. Here we used a high throughput chromosome conformation capture assay (Capture Hi-C) to investigate clustering of DSBs induced at defined loci in the human genome. We unambiguously found that DSBs do cluster but only when induced in transcriptionally active genes. Clustering of damaged genes mainly occurs during the G1 cell cycle phase and coincides with delayed repair. Moreover DSB clustering depends on the MRN complex, as well as the Formin 2 (FMN2) nuclear actin organizer and the LINC (LInker of Nuclear and Cytoplasmic skeleton) complex, suggesting that active mechanisms promote DSB clustering. This work reveals that when damaged, active genes exhibit a very peculiar behavior compared to the rest of the genome, being mostly left unrepaired and clustered in G1 while being repaired by homologous recombination in post-replicative cells.

BackgroundClinical progression of colorectal cancers (CRC) may occur in parallel with distinctive signaling alterations. We designed multidirectional analyses integrating microarray-based data with biostatistics and bioinformatics to elucidate the signaling and metabolic alterations underlying CRC development in the adenoma-carcinoma sequence.Methodology/Principal FindingsStudies were performed on normal mucosa, adenoma, and carcinoma samples obtained during surgery or colonoscopy. Collections of cryostat sections prepared from the tissue samples were evaluated by a pathologist to control the relative cell type content. The measurements were done using Affymetrix GeneChip HG-U133plus2, and probe set data was generated using two normalization algorithms: MAS5.0 and GCRMA with least-variant set (LVS). The data was evaluated using pair-wise comparisons and data decomposition into singular value decomposition (SVD) modes. The method selected for the functional analysis used the Kolmogorov-Smirnov test. Expressional profiles obtained in 105 samples of whole tissue sections were used to establish oncogenic signaling alterations in progression of CRC, while those representing 40 microdissected specimens were used to select differences in KEGG pathways between epithelium and mucosa. Based on a consensus of the results obtained by two normalization algorithms, and two probe set sorting criteria, we identified 14 and 17 KEGG signaling and metabolic pathways that are significantly altered between normal and tumor samples and between benign and malignant tumors, respectively. Several of them were also selected from the raw microarray data of 2 recently published studies (GSE4183 and GSE8671).Conclusion/SignificanceAlthough the proposed strategy is computationally complex and labor–intensive, it may reduce the number of false results.

R-loops have both positive and negative impacts on chromosome functions. To identify toxic R-loops in the human genome, here, we map RNA:DNA hybrids, replication stress markers and DNA double-strand breaks (DSBs) in cells depleted for Topoisomerase I (Top1), an enzyme that relaxes DNA supercoiling and prevents R-loop formation. RNA:DNA hybrids are found at both promoters (TSS) and terminators (TTS) of highly expressed genes. In contrast, the phosphorylation of RPA by ATR is only detected at TTS, which are preferentially replicated in a head-on orientation relative to the direction of transcription. In Top1-depleted cells, DSBs also accumulate at TTS, leading to persistent checkpoint activation, spreading of γ-H2AX on chromatin and global replication fork slowdown. These data indicate that fork pausing at the TTS of highly expressed genes containing R-loops prevents head-on conflicts between replication and transcription and maintains genome integrity in a Top1-dependent manner.

In S. cerevisiae, replication timing is controlled by epigenetic mechanisms restricting the accessibility of origins to limiting initiation factors. About 30% of these origins are located within repetitive DNA sequences such as the ribosomal DNA (rDNA) array, but their regulation is poorly understood. Here, we have investigated how histone deacetylases (HDACs) control the replication program in budding yeast. This analysis revealed that two HDACs, Rpd3 and Sir2, control replication timing in an opposite manner. Whereas Rpd3 delays initiation at late origins, Sir2 is required for the timely activation of early origins. Moreover, Sir2 represses initiation at rDNA origins, whereas Rpd3 counteracts this effect. Remarkably, deletion of SIR2 restored normal replication in rpd3Δ cells by reactivating rDNA origins. Together, these data indicate that HDACs control the replication timing program in budding yeast by modulating the ability of repeated origins to compete with single-copy origins for limiting initiation factors.

Sequencing microRNA, reduced representation sequencing, Hi-C technology and any method requiring the use of in-house barcodes result in sequencing libraries with low initial sequence diversity. Sequencing such data on the Illumina platform typically produces low quality data due to the limitations of the Illumina cluster calling algorithm. Moreover, even in the case of diverse samples, these limitations are causing substantial inaccuracies in multiplexed sample assignment (sample bleeding). Such inaccuracies are unacceptable in clinical applications, and in some other fields (e.g. detection of rare variants). Here, we discuss how both problems with quality of low-diversity samples and sample bleeding are caused by incorrect detection of clusters on the flowcell during initial sequencing cycles. We propose simple software modifications (Long Template Protocol) that overcome this problem. We present experimental results showing that our Long Template Protocol remarkably increases data quality for low diversity samples, as compared with the standard analysis protocol; it also substantially reduces sample bleeding for all samples. For comprehensiveness, we also discuss and compare experimental results from alternative approaches to sequencing low diversity samples. First, we discuss how the low diversity problem, if caused by barcodes, can be avoided altogether at the barcode design stage. Second and third, we present modified guidelines, which are more stringent than the manufacturer’s, for mixing low diversity samples with diverse samples and lowering cluster density, which in our experience consistently produces high quality data from low diversity samples. Fourth and fifth, we present rescue strategies that can be applied when sequencing results in low quality data and when there is no more biological material available. In such cases, we propose that the flowcell be re-hybridized and sequenced again using our Long Template Protocol. Alternatively, we discuss how analysis can be repeated from saved sequencing images using the Long Template Protocol to increase accuracy.

Initiation of eukaryotic chromosome replication follows a spatiotemporal program. The current model suggests that replication origins compete for a limited pool of initiation factors. However, it remains to be answered how these limiting factors are preferentially recruited to early origins. Here, we report that Dbf4 is enriched at early origins through its interaction with forkhead transcription factors Fkh1 and Fkh2. This interaction is mediated by the Dbf4 C terminus and was successfully reconstituted in vitro. An interaction-defective mutant, , phenocopies alleles in terms of origin firing. Remarkably, genome-wide replication profiles reveal that the direct fusion of the DNA-binding domain (DBD) of Fkh1 to Dbf4 restores the Fkh-dependent origin firing but interferes specifically with the pericentromeric origin activation. Furthermore, Dbf4 interacts directly with Sld3 and promotes the recruitment of downstream limiting factors. These data suggest that Fkh1 targets Dbf4 to a subset of noncentromeric origins to promote early replication in a manner that is reminiscent of the recruitment of Dbf4 to pericentromeric origins by Ctf19.