Magdalena Skrzypczak

Nucleotide-resolution DNA double-strand break mapping by next-generation sequencing

We present a genome-wide approach to map DNA double-strand breaks (DSBs) at nucleotide resolution by a method we termed BLESS (direct in situ breaks labeling, enrichment on streptavidin and next-generation sequencing). We validated and tested BLESS using human and mouse cells and different DSBs-inducing agents and sequencing platforms. BLESS was able to detect telomere ends, Sce endonuclease–induced DSBs and complex genome-wide DSB landscapes. As a proof of principle, we characterized the genomic landscape of sensitivity to replication stress in human cells, and we identified >2,000 nonuniformly 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.

Modeling oncogenic signaling in colon tumors by multidirectional analyses of microarray data directed for maximization of analytical reliability.

Background Clinical 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 Findings Studies 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/Significance Although the proposed strategy is computationally complex and labor–intensive, it may reduce the number of false results.

Genome-wide mapping of long-range contacts unveils clustering of DNA double-strand breaks at damaged active genes

The ability of DNA double-strand breaks (DSBs) to cluster in mammalian cells has been a subject of intense debate in recent 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. The results unambiguously demonstrated that DSBs cluster, but only when they are induced within transcriptionally active genes. Clustering of damaged genes occurs primarily 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 linker of nuclear and cytoplasmic skeleton (LINC) complex, thus suggesting that active mechanisms promote clustering. This work reveals that, when damaged, active genes, compared with the rest of the genome, exhibit a distinctive behavior, remaining largely unrepaired and clustered in G1, and being repaired via homologous recombination in postreplicative cells.

Strategies for Achieving High Sequencing Accuracy for Low Diversity Samples and Avoiding Sample Bleeding Using Illumina Platform

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.

Functional features of gene expression profiles differentiating gastrointestinal stromal tumours according to KIT mutations and expression

BackgroundGastrointestinal stromal tumours (GISTs) represent a heterogeneous group of tumours of mesenchymal origin characterized by gain-of-function mutations in KIT or PDGFRA of the type III receptor tyrosine kinase family. Although mutations in either receptor are thought to drive an early oncogenic event through similar pathways, two previous studies reported the mutation-specific gene expression profiles. However, their further conclusions were rather discordant. To clarify the molecular characteristics of differentially expressed genes according to GIST receptor mutations, we combined microarray-based analysis with detailed functional annotations.MethodsTotal RNA was isolated from 29 frozen gastric GISTs and processed for hybridization on GENECHIP® HG-U133 Plus 2.0 microarrays (Affymetrix). KIT and PDGFRA were analyzed by sequencing, while related mRNA levels were analyzed by quantitative RT-PCR.ResultsFifteen and eleven tumours possessed mutations in KIT and PDGFRA, respectively; no mutation was found in three tumours. Gene expression analysis identified no discriminative profiles associated with clinical or pathological parameters, even though expression of hundreds of genes differentiated tumour receptor mutation and expression status. Functional features of genes differentially expressed between the two groups of GISTs suggested alterations in angiogenesis and G-protein-related and calcium signalling.ConclusionOur study has identified novel molecular elements likely to be involved in receptor-dependent GIST development and allowed confirmation of previously published results. These elements may be potential therapeutic targets and novel markers of KIT mutation status.

Ssb1 and Ssb2 cooperate to regulate mouse hematopoietic stem and progenitor cells by resolving replicative stress

Hematopoietic stem and progenitor cells (HSPCs) are vulnerable to endogenous damage and defects in DNA repair can limit their function. The 2 single-stranded DNA (ssDNA) binding proteins SSB1 and SSB2 are crucial regulators of the DNA damage response; however, their overlapping roles during normal physiology are incompletely understood. We generated mice in which both Ssb1 and Ssb2 were constitutively or conditionally deleted. Constitutive Ssb1/Ssb2 double knockout (DKO) caused early embryonic lethality, whereas conditional Ssb1/Ssb2 double knockout (cDKO) in adult mice resulted in acute lethality due to bone marrow failure and intestinal atrophy featuring stem and progenitor cell depletion, a phenotype unexpected from the previously reported single knockout models of Ssb1 or Ssb2 Mechanistically, cDKO HSPCs showed altered replication fork dynamics, massive accumulation of DNA damage, genome-wide double-strand breaks enriched at Ssb-binding regions and CpG islands, together with the accumulation of R-loops and cytosolic ssDNA. Transcriptional profiling of cDKO HSPCs revealed the activation of p53 and interferon (IFN) pathways, which enforced cell cycling in quiescent HSPCs, resulting in their apoptotic death. The rapid cell death phenotype was reproducible in in vitro cultured cDKO-hematopoietic stem cells, which were significantly rescued by nucleotide supplementation or after depletion of p53. Collectively, Ssb1 and Ssb2 control crucial aspects of HSPC function, including proliferation and survival in vivo by resolving replicative stress to maintain genomic stability.

Overactive BRCA1 Affects Presenilin 1 in Induced Pluripotent Stem Cell-Derived Neurons in Alzheimer's Disease

The BRCA1 protein, one of the major players responsible for DNA damage response has recently been linked to Alzheimers disease (AD). Using primary fibroblasts and neurons reprogrammed from induced pluripotent stem cells (iPSC) derived from familial AD (FAD) patients, we studied the role of the BRCA1 protein underlying molecular neurodegeneration. By whole-transcriptome approach, we have found wide range of disturbances in cell cycle and DNA damage response in FAD fibroblasts. This was manifested by significantly increased content of BRCA1 phosphorylated on Ser1524 and abnormal ubiquitination and subcellular distribution of presenilin 1 (PS1). Accordingly, the iPSC-derived FAD neurons showed increased content of BRCA1(Ser1524) colocalized with degraded PS1, accompanied by an enhanced immunostaining pattern of amyloid-β. Finally, overactivation of BRCA1 was followed by an increased content of Cdc25C phosphorylated on Ser216, likely triggering cell cycle re-entry in FAD neurons. This study suggests that overactivated BRCA1 could both influence PS1 turnover leading to amyloid-β pathology and promote cell cycle re-entry-driven cell death of postmitotic neurons in AD.

Genome-wide co-localization of active EGFR and downstream ERK pathway kinases mirrors mitogen-inducible RNA polymerase 2 genomic occupancy

Genome-wide mechanisms that coordinate expression of subsets of functionally related genes are largely unknown. Recent studies show that receptor tyrosine kinases and components of signal transduction cascades including the extracellular signal-regulated protein kinase (ERK), once thought to act predominantly in the vicinity of plasma membrane and in the cytoplasm, can be recruited to chromatin encompassing transcribed genes. Genome-wide distribution of these transducers and their relationship to transcribing RNA polymerase II (Pol2) could provide new insights about co-regulation of functionally related gene subsets. Chromatin immunoprecipitations (ChIP) followed by deep sequencing, ChIP-Seq, revealed that genome-wide binding of epidermal growth factor receptor, EGFR and ERK pathway components at EGF-responsive genes was highly correlated with characteristic mitogen-induced Pol2-profile. Endosomes play a role in intracellular trafficking of proteins including their nuclear import. Immunofluorescence revealed that EGF-activated EGFR, MEK1/2 and ERK1/2 co-localize on endosomes. Perturbation of endosome internalization process, through the depletion of AP2M1 protein, resulted in decreased number of the EGFR containing endosomes and inhibition of Pol2, EGFR/ERK recruitment to EGR1 gene. Thus, mitogen-induced co-recruitment of EGFR/ERK components to subsets of genes, a kinase module possibly pre-assembled on endosome to synchronize their nuclear import, could coordinate genome-wide transcriptional events to ensure effective cell proliferation.

Hypermethylation of TRIM59 and KLF14 Influences Cell Death Signaling in Familial Alzheimer’s Disease

Epigenetic mechanisms play an important role in the development and progression of various neurodegenerative diseases. Abnormal methylation of numerous genes responsible for regulation of transcription, DNA replication, and apoptosis has been linked to Alzheimers disease (AD) pathology. We have recently performed whole transcriptome profiling of familial early-onset Alzheimers disease (fEOAD) patient-derived fibroblasts. On this basis, we demonstrated a strong dysregulation of cell cycle checkpoints and DNA damage response (DDR) in both fibroblasts and reprogrammed neurons. Here, we show that the aging-correlated hypermethylation of KLF14 and TRIM59 genes associates with abnormalities in DNA repair and cell cycle control in fEOAD. Based on the resulting transcriptome networks, we found that the hypermethylation of KLF14 might be associated with epigenetic regulation of the chromatin organization and mRNA processing followed by hypermethylation of TRIM59 likely associated with the G2/M cell cycle phase and p53 role in DNA repair with BRCA1 protein as the key player. We propose that the hypermethylation of KLF14 could constitute a superior epigenetic mechanism for TRIM59 hypermethylation. The methylation status of both genes affects genome stability and might contribute to proapoptotic signaling in AD. Since this study combines data obtained from various tissues from AD patients, it reinforces the view that the genetic methylation status in the blood may be a valuable predictor of molecular processes occurring in affected tissues. Further research is necessary to define a detailed role of TRIM59 and KLF4 in neurodegeneration of neurons.

Comprehensive Mapping of Histone Modifications at DNA Double-Strand Breaks Deciphers Repair Pathway Chromatin Signatures

Summary Double-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.