ColoWeb: a resource for analysis of colocalization of genomic features

Combined epigenetic and immunotherapy for blastic and classical mantle cell lymphoma

Classical MCL (cMCL) constitutes 6-8% of all B cell NHL. Despite recent advances, MCL is incurable except with allogeneic stem cell transplant. Blastic mantle cell lymphoma (bMCL) is a rarer subtype of cMCL associated with an aggressive clinical course and poor treatment response, frequent relapse and poor outcomes. We treated 13 bMCL patients with combined epigenetic and immunotherapy treatment consisting of vorinostat, cladribine and rituximab (SCR). We report an increased OS greater than 40 months with several patients maintaining durable remissions without relapse for longer than 5 years. This is remarkably better then current treatment regimens which in bMCL range from 14.5-24 months with conventional chemotherapy regimens. We demonstrate that the G/A870 CCND1 polymorphism is predictive of blastic disease, nuclear localization of cyclinD1 and response to SCR therapy. The major resistance mechanisms to SCR therapy are loss of CD20 expression and evasion of treatment by sanctuary in the CNS. These data indicate that administration of epigenetic agents improves efficacy of anti-CD20 immunotherapies. This approach is promising in the treatment of MCL and potentially other previously treatment refractory cancers.

GREAP: a comprehensive enrichment analysis software for human genomic regions

The rapid development of genomic high-throughput sequencing has identified a large number of DNA regulatory elements with abundant epigenetics markers, which promotes the rapid accumulation of functional genomic region data. The comprehensively understanding and research of human functional genomic regions is still a relatively urgent work at present. However, the existing analysis tools lack extensive annotation and enrichment analytical abilities for these regions. Here, we designed a novel software, Genomic Region sets Enrichment Analysis Platform (GREAP), which provides comprehensive region annotation and enrichment analysis capabilities. Currently, GREAP supports 85 370 genomic region reference sets, which cover 634 681 107 regions across 11 different data types, including super enhancers, transcription factors, accessible chromatins, etc. GREAP provides widespread annotation and enrichment analysis of genomic regions. To reflect the significance of enrichment analysis, we used the hypergeometric test and also provided a Locus Overlap Analysis. In summary, GREAP is a powerful platform that provides many types of genomic region sets for users and supports genomic region annotations and enrichment analyses. In addition, we developed a customizable genome browser containing >400 000 000 customizable tracks for visualization. The platform is freely available at http://www.liclab.net/Greap/view/index.

Convergence of SIRT1 and ATR signaling to modulate replication origin dormancy

During routine genome duplication, many potential replication origins remain inactive or 'dormant'. Such origin dormancy is achieved, in part, by an interaction with the metabolic sensor SIRT1 deacetylase. We report here that dormant origins are a group of consistent, pre-determined genomic sequences that are distinguished from baseline (i.e. ordinarily active) origins by their preferential association with two phospho-isoforms of the helicase component MCM2. During normal unperturbed cell growth, baseline origins, but not dormant origins, associate with a form of MCM2 that is phosphorylated by DBF4-dependent kinase (DDK) on serine 139 (pS139-MCM2). This association facilitates the initiation of DNA replication from baseline origins. Concomitantly, SIRT1 inhibits Ataxia Telangiectasia and Rad3-related (ATR)-kinase-mediated phosphorylation of MCM2 on serine 108 (pS108-MCM2) by deacetylating the ATR-interacting protein DNA topoisomerase II binding protein 1 (TOPBP1), thereby preventing ATR recruitment to chromatin. In cells devoid of SIRT1 activity, or challenged by replication stress, this inhibition is circumvented, enabling ATR-mediated S108-MCM2 phosphorylation. In turn, pS108-MCM2 enables DDK-mediated phosphorylation on S139-MCM2 and facilitates replication initiation at dormant origins. These observations suggest that replication origin dormancy and activation are regulated by distinct post-translational MCM modifications that reflect a balance between SIRT1 activity and ATR signaling.

Regulation of RNA Polymerase II Transcription Initiation and Elongation by Transcription Factor TFII-I

Transcription by RNA polymerase II (Pol II) is regulated by different processes, including alterations in chromatin structure, interactions between distal regulatory elements and promoters, formation of transcription domains enriched for Pol II and co-regulators, and mechanisms involved in the initiation, elongation, and termination steps of transcription. Transcription factor TFII-I, originally identified as an initiator (INR)-binding protein, contains multiple protein–protein interaction domains and plays diverse roles in the regulation of transcription. Genome-wide analysis revealed that TFII-I associates with expressed as well as repressed genes. Consistently, TFII-I interacts with co-regulators that either positively or negatively regulate the transcription. Furthermore, TFII-I has been shown to regulate transcription pausing by interacting with proteins that promote or inhibit the elongation step of transcription. Changes in TFII-I expression in humans are associated with neurological and immunological diseases as well as cancer. Furthermore, TFII-I is essential for the development of mice and represents a barrier for the induction of pluripotency. Here, we review the known functions of TFII-I related to the regulation of Pol II transcription at the stages of initiation and elongation.

Recommendations for the FAIRification of genomic track metadata

Background: Many types of data from genomic analyses can be represented as genomic tracks, i.e. features linked to the genomic coordinates of a reference genome. Examples of such data are epigenetic DNA methylation data, ChIP-seq peaks, germline or somatic DNA variants, as well as RNA-seq expression levels. Researchers often face difficulties in locating, accessing and combining relevant tracks from external sources, as well as locating the raw data, reducing the value of the generated information.

Description of work: We propose to advance the application of FAIR data principles (Findable, Accessible, Interoperable, and Reusable) to produce searchable metadata for genomic tracks. Findability and Accessibility of metadata can then be ensured by a track search service that integrates globally identifiable metadata from various track hubs in the Track Hub Registry and other relevant repositories. Interoperability and Reusability need to be ensured by the specification and implementation of a basic set of recommendations for metadata. We have tested this concept by developing such a specification in a JSON Schema, called FAIRtracks, and have integrated it into a novel track search service, called TrackFind. We demonstrate practical usage by importing datasets through TrackFind into existing examples of relevant analytical tools for genomic tracks: EPICO and the GSuite HyperBrowser.

Conclusion: We here provide a first iteration of a draft standard for genomic track metadata, as well as the accompanying software ecosystem. It can easily be adapted or extended to future needs of the research community regarding data, methods and tools, balancing the requirements of both data submitters and analytical end-users.

Analytical Approaches for ATAC-seq Data Analysis

ATAC-seq, the assay for transposase-accessible chromatin using sequencing, is a quick and efficient approach to investigating the chromatin accessibility landscape. Investigating chromatin accessibility has broad utility for answering many biological questions, such as mapping nucleosomes, identifying transcription factor binding sites, and measuring differential activity of DNA regulatory elements. Because the ATAC-seq protocol is both simple and relatively inexpensive, there has been a rapid increase in the availability of chromatin accessibility data. Furthermore, advances in ATAC-seq protocols are rapidly extending its breadth to additional experimental conditions, cell types, and species. Accompanying the increase in data, there has also been an explosion of new tools and analytical approaches for analyzing it. Here, we explain the fundamentals of ATAC-seq data processing, summarize common analysis approaches, and review computational tools to provide recommendations for different research questions. This primer provides a starting point and a reference for analysis of ATAC-seq data. © 2020 Wiley Periodicals LLC.

BAMscale: quantification of next-generation sequencing peaks and generation of scaled coverage tracks

Background

Next-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches.

Results

We have developed BAMscale, a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale, we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~ 100 Gb size) in minutes, outperforming currently available tools.

Conclusions

BAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github (https://github.com/ncbi/BAMscale).

LOLAweb: a containerized web server for interactive genomic locus overlap enrichment analysis

The past few years have seen an explosion of interest in understanding the role of regulatory DNA. This interest has driven large-scale production of functional genomics data and analytical methods. One popular analysis is to test for enrichment of overlaps between a query set of genomic regions and a database of region sets. In this way, new genomic data can be easily connected to annotations from external data sources. Here, we present an interactive interface for enrichment analysis of genomic locus overlaps using a web server called LOLAweb. LOLAweb accepts a set of genomic ranges from the user and tests it for enrichment against a database of region sets. LOLAweb renders results in an R Shiny application to provide interactive visualization features, enabling users to filter, sort, and explore enrichment results dynamically. LOLAweb is built and deployed in a Linux container, making it scalable to many concurrent users on our servers and also enabling users to download and run LOLAweb locally.

Metazoan DNA replication origins

DNA replication starts with the opening of DNA at sites called DNA replication origins. From the single sequence-specific DNA replication origin of the small Escherichia coli genome, up to thousands of origins that are necessary to replicate the large human genome, strict sequence specificity has been lost. Nevertheless, genome-wide analyses performed in the recent years, using different mapping methods, demonstrated that there are precise locations along the metazoan genome from which replication initiates. These sites contain relaxed sequence consensus and epigenetic features. There is flexibility in the choice of origins to be used during a given cell cycle, probably imposed by evolution and developmental constraints. Here, we will briefly describe their main features.

The replication initiation determinant protein (RepID) modulates replication by recruiting CUL4 to chromatin

Cell cycle progression in mammals is modulated by two ubiquitin ligase complexes, CRL4 and SCF, which facilitate degradation of chromatin substrates involved in the regulation of DNA replication. One member of the CRL4 complex, the WD-40 containing protein RepID (DCAF14/PHIP), selectively binds and activates a group of replication origins. Here we show that RepID recruits the CRL4 complex to chromatin prior to DNA synthesis, thus playing a crucial architectural role in the proper licensing of chromosomes for replication. In the absence of RepID, cells rely on the alternative ubiquitin ligase, SKP2-containing SCF, to progress through the cell cycle. RepID depletion markedly increases cellular sensitivity to SKP2 inhibitors, which triggered massive genome re-replication. Both RepID and SKP2 interact with distinct, non-overlapping groups of replication origins, suggesting that selective interactions of replication origins with specific CRL components execute the DNA replication program and maintain genomic stability by preventing re-initiation of DNA replication.

RepID has previously been shown to promote origin firing. Here the authors reveal that RepID regulates replication origins via the recruitment of the CRL4 complex, and prevents re-initiation and unscheduled DNA replication.

Epigenomic annotation-based interpretation of genomic data: from enrichment analysis to machine learning

Motivation

One of the goals of functional genomics is to understand the regulatory implications of experimentally obtained genomic regions of interest (ROIs). Most sequencing technologies now generate ROIs distributed across the whole genome. The interpretation of these genome-wide ROIs represents a challenge as the majority of them lie outside of functionally well-defined protein coding regions. Recent efforts by the members of the International Human Epigenome Consortium have generated volumes of functional/regulatory data (reference epigenomic datasets), effectively annotating the genome with epigenomic properties. Consequently, a wide variety of computational tools has been developed utilizing these epigenomic datasets for the interpretation of genomic data.

Results

The purpose of this review is to provide a structured overview of practical solutions for the interpretation of ROIs with the help of epigenomic data. Starting with epigenomic enrichment analysis, we discuss leading tools and machine learning methods utilizing epigenomic and 3D genome structure data. The hierarchy of tools and methods reviewed here presents a practical guide for the interpretation of genome-wide ROIs within an epigenomic context.

Contact

mikhail.dozmorov@vcuhealth.org.

Supplementary information

Supplementary data are available at Bioinformatics online.

FACT is a sensor of DNA torsional stress in eukaryotic cells

Transitions of B-DNA to alternative DNA structures (ADS) can be triggered by negative torsional strain, which occurs during replication and transcription, and may lead to genomic instability. However, how ADS are recognized in cells is unclear. We found that the binding of candidate anticancer drug, curaxin, to cellular DNA results in uncoiling of nucleosomal DNA, accumulation of negative supercoiling and conversion of multiple regions of genomic DNA into left-handed Z-form. Histone chaperone FACT binds rapidly to the same regions via the SSRP1 subunit in curaxin-treated cells. In vitro binding of purified SSRP1 or its isolated CID domain to a methylated DNA fragment containing alternating purine/pyrimidines, which is prone to Z-DNA transition, is much stronger than to other types of DNA. We propose that FACT can recognize and bind Z-DNA or DNA in transition from a B to Z form. Binding of FACT to these genomic regions triggers a p53 response. Furthermore, FACT has been shown to bind to other types of ADS through a different structural domain, which also leads to p53 activation. Thus, we propose that FACT acts as a sensor of ADS formation in cells. Recognition of ADS by FACT followed by a p53 response may explain the role of FACT in DNA damage prevention.

Temporal association of ORCA/LRWD1 to late-firing origins during G1 dictates heterochromatin replication and organization

DNA replication requires the recruitment of a pre-replication complex facilitated by Origin Recognition Complex (ORC) onto the chromatin during G1 phase of the cell cycle. The ORC-associated protein (ORCA/LRWD1) stabilizes ORC on chromatin. Here, we evaluated the genome-wide distribution of ORCA using ChIP-seq during specific time points of G1. ORCA binding sites on the G1 chromatin are dynamic and temporally regulated. ORCA association to specific genomic sites decreases as the cells progressed towards S-phase. The majority of the ORCA-bound sites represent replication origins that also associate with the repressive chromatin marks H3K9me3 and methylated-CpGs, consistent with ORCA-bound origins initiating DNA replication late in S-phase. Further, ORCA directly associates with the repressive marks and interacts with the enzymes that catalyze these marks. Regions that associate with both ORCA and H3K9me3, exhibit diminished H3K9 methylation in ORCA-depleted cells, suggesting a role for ORCA in recruiting the H3K9me3 mark at certain genomic loci. Similarly, DNA methylation is altered at ORCA-occupied sites in cells lacking ORCA. Furthermore, repressive chromatin marks influence ORCA's binding on chromatin. We propose that ORCA coordinates with the histone and DNA methylation machinery to establish a repressive chromatin environment at a subset of origins, which primes them for late replication.

Phosphorylated SIRT1 associates with replication origins to prevent excess replication initiation and preserve genomic stability

Chromatin structure affects DNA replication patterns, but the role of specific chromatin modifiers in regulating the replication process is yet unclear. We report that phosphorylation of the human SIRT1 deacetylase on Threonine 530 (T530-pSIRT1) modulates DNA synthesis. T530-pSIRT1 associates with replication origins and inhibits replication from a group of 'dormant' potential replication origins, which initiate replication only when cells are subject to replication stress. Although both active and dormant origins bind T530-pSIRT1, active origins are distinguished from dormant origins by their unique association with an open chromatin mark, histone H3 methylated on lysine 4. SIRT1 phosphorylation also facilitates replication fork elongation. SIRT1 T530 phosphorylation is essential to prevent DNA breakage upon replication stress and cells harboring SIRT1 that cannot be phosphorylated exhibit a high prevalence of extrachromosomal elements, hallmarks of perturbed replication. These observations suggest that SIRT1 phosphorylation modulates the distribution of replication initiation events to insure genomic stability.

The ChIP-Seq tools and web server: a resource for analyzing ChIP-seq and other types of genomic data

BACKGROUND

ChIP-seq and related high-throughput chromatin profilig assays generate ever increasing volumes of highly valuable biological data. To make sense out of it, biologists need versatile, efficient and user-friendly tools for access, visualization and itegrative analysis of such data.

RESULTS

Here we present the ChIP-Seq command line tools and web server, implementing basic algorithms for ChIP-seq data analysis starting with a read alignment file. The tools are optimized for memory-efficiency and speed thus allowing for processing of large data volumes on inexpensive hardware. The web interface provides access to a large database of public data. The ChIP-Seq tools have a modular and interoperable design in that the output from one application can serve as input to another one. Complex and innovative tasks can thus be achieved by running several tools in a cascade.

CONCLUSIONS

The various ChIP-Seq command line tools and web services either complement or compare favorably to related bioinformatics resources in terms of computational efficiency, ease of access to public data and interoperability with other web-based tools. The ChIP-Seq server is accessible at http://ccg.vital-it.ch/chipseq/ .

Distinct epigenetic features of differentiation-regulated replication origins

BACKGROUND

Eukaryotic genome duplication starts at discrete sequences (replication origins) that coordinate cell cycle progression, ensure genomic stability and modulate gene expression. Origins share some sequence features, but their activity also responds to changes in transcription and cellular differentiation status.

RESULTS

To identify chromatin states and histone modifications that locally mark replication origins, we profiled origin distributions in eight human cell lines representing embryonic and differentiated cell types. Consistent with a role of chromatin structure in determining origin activity, we found that cancer and non-cancer cells of similar lineages exhibited highly similar replication origin distributions. Surprisingly, our study revealed that DNase hypersensitivity, which often correlates with early replication at large-scale chromatin domains, did not emerge as a strong local determinant of origin activity. Instead, we found that two distinct sets of chromatin modifications exhibited strong local associations with two discrete groups of replication origins. The first origin group consisted of about 40,000 regions that actively initiated replication in all cell types and preferentially colocalized with unmethylated CpGs and with the euchromatin markers, H3K4me3 and H3K9Ac. The second group included origins that were consistently active in cells of a single type or lineage and preferentially colocalized with the heterochromatin marker, H3K9me3. Shared origins replicated throughout the S-phase of the cell cycle, whereas cell-type-specific origins preferentially replicated during late S-phase.

CONCLUSIONS

These observations are in line with the hypothesis that differentiation-associated changes in chromatin and gene expression affect the activation of specific replication origins.

LOLA: enrichment analysis for genomic region sets and regulatory elements in R and Bioconductor

Genomic datasets are often interpreted in the context of large-scale reference databases. One approach is to identify significantly overlapping gene sets, which works well for gene-centric data. However, many types of high-throughput data are based on genomic regions. Locus Overlap Analysis (LOLA) provides easy and automatable enrichment analysis for genomic region sets, thus facilitating the interpretation of functional genomics and epigenomics data.

AVAILABILITY AND IMPLEMENTATION

R package available in Bioconductor and on the following website: http://lola.computational-epigenetics.org.