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

Next-Generation Sequencing and Emerging Technologies

Genetic sequencing technologies are evolving at a rapid pace with major implications for research and clinical practice. In this review, the authors provide an updated overview of next-generation sequencing (NGS) and emerging methodologies. NGS has tremendously improved sequencing output while being more time and cost-efficient in comparison to Sanger sequencing. The authors describe short-read sequencing approaches, such as sequencing by synthesis, ion semiconductor sequencing, and nanoball sequencing. Third-generation long-read sequencing now promises to overcome many of the limitations of short-read sequencing, such as the ability to reliably resolve repeat sequences and large genomic rearrangements. By combining complementary methods with massively parallel DNA sequencing, a greater insight into the biological context of disease mechanisms is now possible. Emerging methodologies, such as advances in nanopore technology, in situ nucleic acid sequencing, and microscopy-based sequencing, will continue the rapid evolution of this area. These new technologies hold many potential applications for hematological disorders, with the promise of precision and personalized medical care in the future.

ROS-induced ribosome impairment underlies ZAKα-mediated metabolic decline in obesity and aging

The ribotoxic stress response (RSR) is a signaling pathway in which the p38- and c-Jun N-terminal kinase (JNK)-activating mitogen-activated protein kinase kinase kinase (MAP3K) ZAKα senses stalling and/or collision of ribosomes. Here, we show that reactive oxygen species (ROS)-generating agents trigger ribosomal impairment and ZAKα activation. Conversely, zebrafish larvae deficient for ZAKα are protected from ROS-induced pathology. Livers of mice fed a ROS-generating diet exhibit ZAKα-activating changes in ribosomal elongation dynamics. Highlighting a role for the RSR in metabolic regulation, ZAK-knockout mice are protected from developing high-fat high-sugar (HFHS) diet-induced blood glucose intolerance and liver steatosis. Finally, ZAK ablation slows animals from developing the hallmarks of metabolic aging. Our work highlights ROS-induced ribosomal impairment as a physiological activation signal for ZAKα that underlies metabolic adaptation in obesity and aging.

Ribosome stalling is a signal for metabolic regulation by the ribotoxic stress response

Impairment of translation can lead to collisions of ribosomes, which constitute an activation platform for several ribosomal stress-surveillance pathways. Among these is the ribotoxic stress response (RSR), where ribosomal sensing by the MAP3K ZAKα leads to activation of p38 and JNK kinases. Despite these insights, the physiological ramifications of ribosomal impairment and downstream RSR signaling remain elusive. Here, we show that stalling of ribosomes is sufficient to activate ZAKα. In response to amino acid deprivation and full nutrient starvation, RSR impacts on the ensuing metabolic responses in cells, nematodes, and mice. The RSR-regulated responses in these model systems include regulation of AMPK and mTOR signaling, survival under starvation conditions, stress hormone production, and regulation of blood sugar control. In addition, ZAK^-/- male mice present a lean phenotype. Our work highlights impaired ribosomes as metabolic signals and demonstrates a role for RSR signaling in metabolic regulation.

Productive visualization of high-throughput sequencing data using the SeqCode open portable platform

Large-scale sequencing techniques to chart genomes are entirely consolidated. Stable computational methods to perform primary tasks such as quality control, read mapping, peak calling, and counting are likewise available. However, there is a lack of uniform standards for graphical data mining, which is also of central importance. To fill this gap, we developed SeqCode, an open suite of applications that analyzes sequencing data in an elegant but efficient manner. Our software is a portable resource written in ANSI C that can be expected to work for almost all genomes in any computational configuration. Furthermore, we offer a user-friendly front-end web server that integrates SeqCode functions with other graphical analysis tools. Our analysis and visualization toolkit represents a significant improvement in terms of performance and usability as compare to other existing programs. Thus, SeqCode has the potential to become a key multipurpose instrument for high-throughput professional analysis; further, it provides an extremely useful open educational platform for the world-wide scientific community. SeqCode website is hosted at http://ldicrocelab.crg.eu, and the source code is freely distributed at https://github.com/eblancoga/seqcode.

Computational identification and experimental characterization of preferred downstream positions in human core promoters

Metazoan core promoters, which direct the initiation of transcription by RNA polymerase II (Pol II), may contain short sequence motifs termed core promoter elements/motifs (e.g. the TATA box, initiator (Inr) and downstream core promoter element (DPE)), which recruit Pol II via the general transcription machinery. The DPE was discovered and extensively characterized in Drosophila, where it is strictly dependent on both the presence of an Inr and the precise spacing from it. Since the Drosophila DPE is recognized by the human transcription machinery, it is most likely that some human promoters contain a downstream element that is similar, though not necessarily identical, to the Drosophila DPE. However, only a couple of human promoters were shown to contain a functional DPE, and attempts to computationally detect human DPE-containing promoters have mostly been unsuccessful. Using a newly-designed motif discovery strategy based on Expectation-Maximization probabilistic partitioning algorithms, we discovered preferred downstream positions (PDP) in human promoters that resemble the Drosophila DPE. Available chromatin accessibility footprints revealed that Drosophila and human Inr+DPE promoter classes are not only highly structured, but also similar to each other, particularly in the proximal downstream region. Clustering of the corresponding sequence motifs using a neighbor-joining algorithm strongly suggests that canonical Inr+DPE promoters could be common to metazoan species. Using reporter assays we demonstrate the contribution of the identified downstream positions to the function of multiple human promoters. Furthermore, we show that alteration of the spacing between the Inr and PDP by two nucleotides results in reduced promoter activity, suggesting a spacing dependency of the newly discovered human PDP on the Inr. Taken together, our strategy identified novel functional downstream positions within human core promoters, supporting the existence of DPE-like motifs in human promoters.

Transcription of genes by the RNA polymerase II enzyme initiates at a genomic region termed the core promoter. The core promoter is a regulatory region that may contain diverse short DNA sequence motifs/elements that confer specific properties to it. Interestingly, core promoter motifs can be located both upstream and downstream of the transcription start site. Variable compositions of core promoter elements were identified. The initiator (Inr) motif and the downstream core promoter element (DPE) is a combination of elements that has been identified and extensively characterized in fruit flies. Although a few Inr+DPE -containing human promoters were identified, the presence of transcriptionally important downstream core promoter positions within human promoters has been a matter of controversy in the literature. Here, using a newly-designed motif discovery strategy, we discovered preferred downstream positions in human promoters that resemble fruit fly DPE. Clustering of the corresponding sequence motifs in eight additional species indicated that such promoters could be common to multicellular non-plant organisms. Importantly, functional characterization of the newly discovered preferred downstream positions supports the existence of Inr+DPE-containing promoters in human genes.

AnnoMiner is a new web-tool to integrate epigenetics, transcription factor occupancy and transcriptomics data to predict transcriptional regulators

Gene expression regulation requires precise transcriptional programs, led by transcription factors in combination with epigenetic events. Recent advances in epigenomic and transcriptomic techniques provided insight into different gene regulation mechanisms. However, to date it remains challenging to understand how combinations of transcription factors together with epigenetic events control cell-type specific gene expression. We have developed the AnnoMiner web-server, an innovative and flexible tool to annotate and integrate epigenetic, and transcription factor occupancy data. First, AnnoMiner annotates user-provided peaks with gene features. Second, AnnoMiner can integrate genome binding data from two different transcriptional regulators together with gene features. Third, AnnoMiner offers to explore the transcriptional deregulation of genes nearby, or within a specified genomic region surrounding a user-provided peak. AnnoMiner’s fourth function performs transcription factor or histone modification enrichment analysis for user-provided gene lists by utilizing hundreds of public, high-quality datasets from ENCODE for the model organisms human, mouse, Drosophila and C. elegans. Thus, AnnoMiner can predict transcriptional regulators for a studied process without the strict need for chromatin data from the same process. We compared AnnoMiner to existing tools and experimentally validated several transcriptional regulators predicted by AnnoMiner to indeed contribute to muscle morphogenesis in Drosophila. AnnoMiner is freely available at http://chimborazo.ibdm.univ-mrs.fr/AnnoMiner/.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

Pausing of transcribing RNA polymerase is regulated and creates opportunities to control gene expression. Research in metazoans has so far mainly focused on RNA polymerase II (Pol II) promoter-proximal pausing leaving the pervasive nature of pausing and its regulatory potential in mammalian cells unclear. Here, we developed a pause detecting algorithm (PDA) for nucleotide-resolution occupancy data and a new native elongating transcript sequencing approach, termed nested NET-seq, that strongly reduces artifactual peaks commonly misinterpreted as pausing sites. Leveraging PDA and nested NET-seq reveal widespread genome-wide Pol II pausing at single-nucleotide resolution in human cells. Notably, the majority of Pol II pauses occur outside of promoter-proximal gene regions primarily along the gene-body of transcribed genes. Sequence analysis combined with machine learning modeling reveals DNA sequence properties underlying widespread transcriptional pausing including a new pause motif. Interestingly, key sequence determinants of RNA polymerase pausing are conserved between human cells and bacteria. These studies indicate pervasive sequence-induced transcriptional pausing in human cells and the knowledge of exact pause locations implies potential functional roles in gene expression.

Plasma Membrane Calcium ATPase Regulates Stoichiometry of CD4+ T-Cell Compartments

Immune responses involve mobilization of T cells within naïve and memory compartments. Tightly regulated Ca²⁺ levels are essential for balanced immune outcomes. How Ca²⁺ contributes to regulating compartment stoichiometry is unknown. Here, we show that plasma membrane Ca²⁺ ATPase 4 (PMCA4) is differentially expressed in human CD4⁺ T compartments yielding distinct store operated Ca²⁺ entry (SOCE) profiles. Modulation of PMCA4 yielded a more prominent increase of SOCE in memory than in naïve CD4⁺ T cell. Interestingly, downregulation of PMCA4 reduced the effector compartment fraction and led to accumulation of cells in the naïve compartment. In silico analysis and chromatin immunoprecipitation point towards Ying Yang 1 (YY1) as a transcription factor regulating PMCA4 expression. Analyses of PMCA and YY1 expression patterns following activation and of PMCA promoter activity following downregulation of YY1 highlight repressive role of YY1 on PMCA expression. Our findings show that PMCA4 adapts Ca²⁺ levels to cellular requirements during effector and quiescent phases and thereby represent a potential target to intervene with the outcome of the immune response.

Transcriptional Profiling of STAT1 Gain-of-Function Reveals Common and Mutation-Specific Fingerprints

STAT1 gain-of-function (GOF) is a primary immunodeficiency typically characterized by chronic mucocutaneous candidiasis (CMC), recurrent respiratory infections, and autoimmunity. Less commonly, also immunodysregulation polyendocrinopathy enteropathy X-linked (IPEX)-like syndromes with CMC, and combined immunodeficiency without CMC have been described. Recently, our group and others have shown that different mutation-specific mechanisms underlie STAT1 GOF in vitro, including faster nuclear accumulation (R274W), and reduced mobility (R321, N574I) to near immobility in the nucleus (T419R) upon IFNγ stimulation. In this work, we evaluated the transcriptomic fingerprint of the aforementioned STAT1 GOF mutants (R274W, R321S, T419R, and N574I) relative to STAT1 wild-type upon IFNγ stimulation in an otherwise isogenic cell model. The majority of genes up-regulated in wild-type STAT1 cells were significantly more up-regulated in cells expressing GOF mutants, except for T419R. In addition to the common interferon regulated genes (IRG), STAT1 GOF mutants up-regulated an additional set of genes, that were in part shared with other GOF mutants or mutation-specific. Overall, R274W and R321S transcriptomes clustered with STAT1 WT, while T419R and N574I had a more distinct fingerprint. We observed reduced frequency of canonical IFNγ activation site (GAS) sequences in promoters of genes up-regulated by all the STAT1 GOF mutants, suggesting loss of DNA binding specificity for the canonical GAS consensus. Interestingly, the T419R mutation, expected to directly increase the affinity for DNA, showed the most pronounced effects on the transcriptome. T419R STAT1 dysregulated more non-IRG than the other GOF mutants and fewer GAS or degenerate GAS promotor sequences could be found in the promoter regions of these genes. In conclusion, our work confirms hyperactivation of common sets of IFNγ-induced genes in STAT1 GOF with additional dysregulation of mutation-specific genes, in line with the earlier observed mutation-specific mechanisms. Binding to more degenerate GAS sequences is proposed as a mechanism toward transcriptional dysregulation in R274W, R321S, and N574I. For T419R, an increased interaction with the DNA is suggested to result in a broader and less GAS-specific response. Our work indicates that multiple routes leading to STAT1 GOF are associated with common and private transcriptomic fingerprints, which may contribute to the phenotypic variation observed in vivo.

Readthrough of stop codons under limiting ABCE1 concentration involves frameshifting and inhibits nonsense-mediated mRNA decay

To gain insight into the mechanistic link between translation termination and nonsense-mediated mRNA decay (NMD), we depleted the ribosome recycling factor ABCE1 in human cells, resulting in an upregulation of NMD-sensitive mRNAs. Suppression of NMD on these mRNAs occurs prior to their SMG6-mediated endonucleolytic cleavage. ABCE1 depletion caused ribosome stalling at termination codons (TCs) and increased ribosome occupancy in 3′ UTRs, implying enhanced TC readthrough. ABCE1 knockdown indeed increased the rate of readthrough and continuation of translation in different reading frames, providing a possible explanation for the observed NMD inhibition, since enhanced readthrough displaces NMD activating proteins from the 3′ UTR. Our results indicate that stalling at TCs triggers ribosome collisions and activates ribosome quality control. Collectively, we show that improper translation termination can lead to readthrough of the TC, presumably due to ribosome collisions pushing the stalled ribosomes into the 3′ UTR, where it can resume translation in-frame as well as out-of-frame.

A Functional Link between Nuclear RNA Decay and Transcriptional Control Mediated by the Polycomb Repressive Complex 2

Pluripotent embryonic stem cells (ESCs) constitute an essential cellular niche sustained by epigenomic and transcriptional regulation. Any role of post-transcriptional processes remains less explored. Here, we identify a link between nuclear RNA levels, regulated by the poly(A) RNA exosome targeting (PAXT) connection, and transcriptional control by the polycomb repressive complex 2 (PRC2). Knockout of the PAXT component ZFC3H1 impairs mouse ESC differentiation. In addition to the upregulation of bona fide PAXT substrates, Zfc3h1^−/− cells abnormally express developmental genes usually repressed by PRC2. Such de-repression is paralleled by decreased PRC2 binding to chromatin and low PRC2-directed H3K27 methylation. PRC2 complex stability is compromised in Zfc3h1^−/− cells with elevated levels of unspecific RNA bound to PRC2 components. We propose that excess RNA hampers PRC2 function through its sequestration from DNA. Our results highlight the importance of balancing nuclear RNA levels and demonstrate the capacity of bulk RNA to regulate chromatin-associated proteins.

•

Depletion of ZFC3H1 in mouse ESCs results in differentiation defects

•

PRC2 target genes are deregulated in Zfc3h1^−/− cells

•

Chromatin binding of PRC2 and H3K27me3 is reduced in Zfc3h1^−/− cells

•

Increased binding of RNA impairs PRC2 complex stability

EPD in 2020: enhanced data visualization and extension to ncRNA promoters

The Eukaryotic Promoter Database (EPD), available online at https://epd.epfl.ch, provides accurate transcription start site (TSS) information for promoters of 15 model organisms plus corresponding functional genomics data that can be viewed in a genome browser, queried or analyzed via web interfaces, or exported in standard formats (FASTA, BED, CSV) for subsequent analysis with other tools. Recent work has focused on the improvement of the EPD promoter viewers, which use the UCSC Genome Browser as visualization platform. Thousands of high-resolution tracks for CAGE, ChIP-seq and similar data have been generated and organized into public track hubs. Customized, reproducible promoter views, combining EPD-supplied tracks with native UCSC Genome Browser tracks, can be accessed from the organism summary pages or from individual promoter entries. Moreover, thanks to recent improvements and stabilization of ncRNA gene catalogs, we were able to release promoter collections for certain classes of ncRNAs from human and mouse. Furthermore, we developed automatic computational protocols to assign orphan TSS peaks to downstream genes based on paired-end (RAMPAGE) TSS mapping data, which enabled us to add nearly 9000 new entries to the human promoter collection. Since our last article in this journal, EPD was extended to five more model organisms: rhesus monkey, rat, dog, chicken and Plasmodium falciparum.

Anti-adipogenic signals at the onset of obesity-related inflammation in white adipose tissue

Chronic inflammation that affects primarily metabolic organs, such as white adipose tissue (WAT), is considered as a major cause of human obesity-associated co-morbidities. However, the molecular mechanisms initiating this inflammation in WAT are poorly understood. By combining transcriptomics, ChIP-seq and modeling approaches, we studied the global early and late responses to a high-fat diet (HFD) in visceral (vWAT) and subcutaneous (scWAT) AT, the first being more prone to obesity-induced inflammation. HFD rapidly triggers proliferation of adipocyte precursors within vWAT. However, concomitant antiadipogenic signals limit vWAT hyperplastic expansion by interfering with the differentiation of proliferating adipocyte precursors. Conversely, in scWAT, residing beige adipocytes lose their oxidizing properties and allow storage of excessive fatty acids. This phase is followed by tissue hyperplastic growth and increased angiogenic signals, which further enable scWAT expansion without generating inflammation. Our data indicate that scWAT and vWAT differential ability to modulate adipocyte number and differentiation in response to obesogenic stimuli has a crucial impact on the different susceptibility to obesity-related inflammation of these adipose tissue depots.

CSA: a web service for the complete process of ChIP-Seq analysis

BACKGROUND

Chromatin immunoprecipitation sequencing (ChIP-seq) is a technology that combines chromatin immunoprecipitation (ChIP) with next generation of sequencing technology (NGS) to analyze protein interactions with DNA. At present, most ChIP-seq analysis tools adopt the command line, which lacks user-friendly interfaces. Although some web services with graphical interfaces have been developed for ChIP-seq analysis, these sites cannot provide a comprehensive analysis of ChIP-seq from raw data to downstream analysis.

RESULTS

In this study, we develop a web service for the whole process of ChIP-Seq Analysis (CSA), which covers mapping, quality control, peak calling, and downstream analysis. In addition, CSA provides a customization function for users to define their own workflows. And the visualization of mapping, peak calling, motif finding, and pathway analysis results are also provided in CSA. For the different types of ChIP-seq datasets, CSA can provide the corresponding tool to perform the analysis. Moreover, CSA can detect differences in ChIP signals between ChIP samples and controls to identify absolute binding sites.

CONCLUSIONS

The two case studies demonstrate the effectiveness of CSA, which can complete the whole procedure of ChIP-seq analysis. CSA provides a web interface for users, and implements the visualization of every analysis step. The website of CSA is available at http://CompuBio.csu.edu.cn.

PWMScan: a fast tool for scanning entire genomes with a position-specific weight matrix

Summary

Transcription factors regulate gene expression by binding to specific short DNA sequences of 5-20 bp to regulate the rate of transcription of genetic information from DNA to messenger RNA. We present PWMScan, a fast web-based tool to scan server-resident genomes for matches to a user-supplied PWM or transcription factor binding site model from a public database.

Availability and implementation

The web server and source code are available at http://ccg.vital-it.ch/pwmscan and https://sourceforge.net/projects/pwmscan, respectively.

Supplementary information

Supplementary data are available at Bioinformatics online.

MGA repository: a curated data resource for ChIP-seq and other genome annotated data

The Mass Genome Annotation (MGA) repository is a resource designed to store published next generation sequencing data and other genome annotation data (such as gene start sites, SNPs, etc.) in a completely standardised format. Each sample has undergone local processing in order the meet the strict MGA format requirements. The original data source, the reformatting procedure and the biological characteristics of the samples are described in an accompanying documentation file manually edited by data curators. 10 model organisms are currently represented: Homo sapiens, Mus musculus, Danio rerio, Drosophila melanogaster, Apis mellifera, Caenorhabditis elegans, Arabidopsis thaliana, Zea mays, Saccharomyces cerevisiae and Schizosaccharomyces pombe. As of today, the resource contains over 24 000 samples. In conjunction with other tools developed by our group (the ChIP-Seq and SSA servers), it allows users to carry out a great variety of analysis task with MGA samples, such as making aggregation plots and heat maps for selected genomic regions, finding peak regions, generating custom tracks for visualizing genomic features in a UCSC genome browser window, or downloading chromatin data in a table format suitable for local processing with more advanced statistical analysis software such as R. Home page: http://ccg.vital-it.ch/mga/.

TrawlerWeb: an online de novo motif discovery tool for next-generation sequencing datasets

Background

A strong focus of the post-genomic era is mining of the non-coding regulatory genome in order to unravel the function of regulatory elements that coordinate gene expression (Nat 489:57–74, 2012; Nat 507:462–70, 2014; Nat 507:455–61, 2014; Nat 518:317–30, 2015). Whole-genome approaches based on next-generation sequencing (NGS) have provided insight into the genomic location of regulatory elements throughout different cell types, organs and organisms. These technologies are now widespread and commonly used in laboratories from various fields of research. This highlights the need for fast and user-friendly software tools dedicated to extracting cis-regulatory information contained in these regulatory regions; for instance transcription factor binding site (TFBS) composition. Ideally, such tools should not require prior programming knowledge to ensure they are accessible for all users.

Results

We present TrawlerWeb, a web-based version of the Trawler_standalone tool (Nat Methods 4:563–5, 2007; Nat Protoc 5:323–34, 2010), to allow for the identification of enriched motifs in DNA sequences obtained from next-generation sequencing experiments in order to predict their TFBS composition. TrawlerWeb is designed for online queries with standard options common to web-based motif discovery tools. In addition, TrawlerWeb provides three unique new features: 1) TrawlerWeb allows the input of BED files directly generated from NGS experiments, 2) it automatically generates an input-matched biologically relevant background, and 3) it displays resulting conservation scores for each instance of the motif found in the input sequences, which assists the researcher in prioritising the motifs to validate experimentally. Finally, to date, this web-based version of Trawler_standalone remains the fastest online de novo motif discovery tool compared to other popular web-based software, while generating predictions with high accuracy.

Conclusions

TrawlerWeb provides users with a fast, simple and easy-to-use web interface for de novo motif discovery. This will assist in rapidly analysing NGS datasets that are now being routinely generated. TrawlerWeb is freely available and accessible at: http://trawler.erc.monash.edu.au.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4630-0) contains supplementary material, which is available to authorized users.