Identifying differential transcription factor binding in ChIP-seq

Nuclear lamin A/C phosphorylation by loss of androgen receptor leads to cancer-associated fibroblast activation

Alterations in nuclear structure and function are hallmarks of cancer cells. Little is known about these changes in Cancer-Associated Fibroblasts (CAFs), crucial components of the tumor microenvironment. Loss of the androgen receptor (AR) in human dermal fibroblasts (HDFs), which triggers early steps of CAF activation, leads to nuclear membrane changes and micronuclei formation, independent of cellular senescence. Similar changes occur in established CAFs and are reversed by restoring AR activity. AR associates with nuclear lamin A/C, and its loss causes lamin A/C nucleoplasmic redistribution. AR serves as a bridge between lamin A/C and the protein phosphatase PPP1. Loss of AR decreases lamin-PPP1 association and increases lamin A/C phosphorylation at Ser 301, a characteristic of CAFs. Phosphorylated lamin A/C at Ser 301 binds to the regulatory region of CAF effector genes of the myofibroblast subtype. Expression of a lamin A/C Ser301 phosphomimetic mutant alone can transform normal fibroblasts into tumor-promoting CAFs.

Characterization of the distribution and dynamics of chromatin states in the C. elegans germline reveals substantial H3K4me3 remodeling during oogenesis

Chromatin organization in the C. elegans germline is tightly regulated and critical for germ cell differentiation. Although certain germline epigenetic regulatory mechanisms have been identified, how they influence chromatin structure and ultimately gene expression remains unclear, in part because most genomic studies have focused on data collected from intact worms comprising both somatic and germline tissues. We therefore analyzed histone modification and chromatin accessibility data from isolated germ nuclei representing undifferentiated proliferating and meiosis I populations to define chromatin states. We correlated these states with overall transcript abundance, spatiotemporal expression patterns, and the function of small RNA pathways. Because the essential role of the germline is to transmit genetic information and establish gene expression in the early embryo, we compared epigenetic and transcriptomic profiles from undifferentiated germ cells to those of embryos to define the epigenetic changes during this developmental transition. The active histone modification H3K4me3 shows particularly dynamic remodeling as germ cells differentiate into oocytes, which suggests a mechanism for establishing early transcription of essential genes during zygotic genome activation. This analysis highlights the dynamism of the chromatin landscape across developmental transitions and provides a resource for future investigation into epigenetic regulatory mechanisms in germ cells.

The single-cell transcriptomic atlas and RORA-mediated 3D epigenomic remodeling in driving corneal epithelial differentiation

Proper differentiation of corneal epithelial cells (CECs) from limbal stem/progenitor cells (LSCs) is required for maintenance of ocular homeostasis and clear vision. Here, using a single-cell transcriptomic atlas, we delineate the comprehensive and refined molecular regulatory dynamics during human CEC development and differentiation. We find that RORA is a CEC-specific molecular switch that initiates and drives LSCs to differentiate into mature CECs by activating PITX1. RORA dictates CEC differentiation by establishing CEC-specific enhancers and chromatin interactions between CEC gene promoters and distal regulatory elements. Conversely, RORA silences LSC-specific promoters and disrupts promoter-anchored chromatin loops to turn off LSC genes. Collectively, our work provides detailed and comprehensive insights into the transcriptional dynamics and RORA-mediated epigenetic remodeling underlying human corneal epithelial differentiation.

Expression Analysis of Genes Corresponding to Mucins and Their Glycans from Cervical Tissue Using RNA Sequencing

In the female reproductive tract, mucin proteins are the main component of mucus secreted by cervical goblet cells and play an essential role in many biological functions. They act as a medium for lubrication and mainting a cervical mucosal barrier against ascending pathogens from the vagina while also allowing sperm migration. The expression of mucin genes as well as the levels of O-glycosylation changes across the oestrous cycle. Detection and characterization of mucins and their glycans is important to understand the interface between the external and the internal environment, as the cervical epithelium represents the first line of defense against infections of the upper reproductive tract. Advances in the field of molecular biology have made possible to study differences in mucin and glycan gene expression which can help to understand impeded sperm transport as well as variation in the susceptibility to infection. This chapter discusses procedures relevant for both animals and humans on how to recover cervical tissue and perform a gene expression analysis of genes corresponding to mucins and their glycans using RNA sequencing.

SGF29 nuclear condensates reinforce cellular aging

Phase separation, a biophysical segregation of subcellular milieus referred as condensates, is known to regulate transcription, but its impacts on physiological processes are less clear. Here, we demonstrate the formation of liquid-like nuclear condensates by SGF29, a component of the SAGA transcriptional coactivator complex, during cellular senescence in human mesenchymal progenitor cells (hMPCs) and fibroblasts. The Arg 207 within the intrinsically disordered region is identified as the key amino acid residue for SGF29 to form phase separation. Through epigenomic and transcriptomic analysis, our data indicated that both condensate formation and H3K4me3 binding of SGF29 are essential for establishing its precise chromatin location, recruiting transcriptional factors and co-activators to target specific genomic loci, and initiating the expression of genes associated with senescence, such as CDKN1A. The formation of SGF29 condensates alone, however, may not be sufficient to drive H3K4me3 binding or achieve transactivation functions. Our study establishes a link between phase separation and aging regulation, highlighting nuclear condensates as a functional unit that facilitate shaping transcriptional landscapes in aging.

Comprehensive single cell aging atlas of mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. We investigate how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory, or cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveal co-localization of aged immune and epithelial cells in situ. Lastly, transcriptional signatures of aging mammary cells are found in human breast tumors, suggesting mechanistic links between aging and cancer. Together, these data uncover that epithelial, immune, and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment, and neoplasia risk.

Nuclear speckles regulate HIF-2α programs and correlate with patient survival in kidney cancer

Nuclear speckles are membrane-less bodies within the cell nucleus enriched in RNA biogenesis, processing, and export factors. In this study we investigated speckle phenotype variation in human cancer, finding a reproducible speckle signature, based on RNA expression of speckle-resident proteins, across >20 cancer types. Of these, clear cell renal cell carcinoma (ccRCC) exhibited a clear correlation between the presence of this speckle expression signature, imaging-based speckle phenotype, and clinical outcomes. ccRCC is typified by hyperactivation of the HIF-2α transcription factor, and we demonstrate here that HIF-2α drives physical association of a select subset of its target genes with nuclear speckles. Disruption of HIF-2α-driven speckle association via deletion of its speckle targeting motifs (STMs)-defined in this study-led to defective induction of speckle-associating HIF-2α target genes without impacting non-speckle-associating HIF-2α target genes. We further identify the RNA export complex, TREX, as being specifically altered in speckle signature, and knockdown of key TREX component, ALYREF, also compromises speckle-associated gene expression. By integrating tissue culture functional studies with tumor genomic and imaging analysis, we show that HIF-2α gene regulatory programs are impacted by specific manipulation of speckle phenotype and by abrogation of speckle targeting abilities of HIF-2α. These findings suggest that, in ccRCC, a key biological function of nuclear speckles is to modulate expression of a specific subset of HIF-2α-regulated target genes that, in turn, influence patient outcomes. We also identify STMs in other transcription factors, suggesting that DNA-speckle targeting may be a general mechanism of gene regulation.

Tissue memory relies on stem cell priming in distal undamaged areas

Epithelial cells that participated in wound repair elicit a more efficient response to future injuries, which is believed to be locally restricted. Here we show that cell adaptation resulting from a localized tissue damage has a wide spatial impact at a scale not previously appreciated. We demonstrate that a specific stem cell population, distant from the original injury, originates long-lasting wound memory progenitors residing in their own niche. Notably, these distal memory cells have not taken part in the first healing but become intrinsically pre-activated through priming. This cell state, maintained at the chromatin and transcriptional level, leads to an enhanced wound repair that is partially recapitulated through epigenetic perturbation. Importantly wound memory has long-term harmful consequences, exacerbating tumourigenesis. Overall, we show that sub-organ-scale adaptation to injury relies on spatially organized memory-dedicated progenitors, characterized by an actionable cell state that establishes an epigenetic field cancerization and predisposes to tumour onset.

Statistical Analysis in ChIP-seq-Related Applications

Chromatin immunoprecipitation sequencing (ChIP-seq) has been widely performed to identify protein binding information along the genome. The sequencing protocol is quite flexible and mature to measure different types of protein binding as long as sequencing parameters are properly tailored to accommodate protein features. Two distinct types of protein binding are point-source-like binding by transcription factors and diffused-distribution binding by histone modifications. Consequently, statistical approaches have been proposed to address ChIP-seq-related questions according to different protein features. In this chapter, we briefly summarize statistical principles, approaches, and tools that are widely implemented in modeling ChIP-seq data, from raw data quality control to final result reporting. We discuss the key solutions in addressing eight routine questions in ChIP-seq applications. We also include discussion on approaches fitting unique data features in different ChIP-seq types. We hope this chapter will serve as a brief guide, especially for ChIP-seq beginners, to provide them with a high-level overview to understand and design processing plans for their ChIP-seq experiments.

Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

Cis-regulatory chromatin loops analysis identifies GRHL3 as a master regulator of surface epithelium commitment

Understanding the regulatory network of cell fate acquisition remains a major challenge. Using the induction of surface epithelium (SE) from human embryonic stem cells as a paradigm, we show that the dynamic changes in morphology-related genes (MRGs) closely correspond to SE fate transitions. The marked remodeling of cytoskeleton indicates the initiation of SE differentiation. By integrating promoter interactions, epigenomic features, and transcriptome, we delineate an SE-specific cis-regulatory network and identify grainyhead-like 3 (GRHL3) as an initiation factor sufficient to drive SE commitment. Mechanically, GRHL3 primes the SE chromatin accessibility landscape and activates SE-initiating gene expression. In addition, the evaluation of GRHL3-mediated promoter interactions unveils a positive feedback loop of GRHL3 and bone morphogenetic protein 4 on SE fate decisions. Our work proposes a concept that MRGs could be used to identify cell fate transitions and provides insights into regulatory principles of SE lineage development and stem cell-based regenerative medicine.

Bridging Glycomics and Genomics: New Uses of Functional Genetics in the Study of Cellular Glycosylation

All living cells are coated with a diverse collection of carbohydrate molecules called glycans. Glycans are key regulators of cell behavior and important therapeutic targets for human disease. Unlike proteins, glycans are not directly templated by discrete genes. Instead, they are produced through multi-gene pathways that generate a heterogenous array of glycoprotein and glycolipid antigens on the cell surface. This genetic complexity has sometimes made it challenging to understand how glycosylation is regulated and how it becomes altered in disease. Recent years, however, have seen the emergence of powerful new functional genomics technologies that allow high-throughput characterization of genetically complex cellular phenotypes. In this review, we discuss how these techniques are now being applied to achieve a deeper understanding of glyco-genomic regulation. We highlight specifically how methods like ChIP-seq, RNA-seq, CRISPR genomic screening and scRNA-seq are being used to map the genomic basis for various cell-surface glycosylation states in normal and diseased cell types. We also offer a perspective on how emerging functional genomics technologies are likely to create further opportunities for studying cellular glycobiology in the future. Taken together, we hope this review serves as a primer to recent developments at the glycomics-genomics interface.

Comprehensive assessment of differential ChIP-seq tools guides optimal algorithm selection

Background

The analysis of chromatin binding patterns of proteins in different biological states is a main application of chromatin immunoprecipitation followed by sequencing (ChIP-seq). A large number of algorithms and computational tools for quantitative comparison of ChIP-seq datasets exist, but their performance is strongly dependent on the parameters of the biological system under investigation. Thus, a systematic assessment of available computational tools for differential ChIP-seq analysis is required to guide the optimal selection of analysis tools based on the present biological scenario.

Results

We created standardized reference datasets by in silico simulation and sub-sampling of genuine ChIP-seq data to represent different biological scenarios and binding profiles. Using these data, we evaluated the performance of 33 computational tools and approaches for differential ChIP-seq analysis. Tool performance was strongly dependent on peak size and shape as well as on the scenario of biological regulation.

Conclusions

Our analysis provides unbiased guidelines for the optimized choice of software tools in differential ChIP-seq analysis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02686-y.

The acetyltransferase p300 is recruited in trans to multiple enhancer sites by lncSmad7

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

Differential ATAC-seq and ChIP-seq peak detection using ROTS

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

Reduction of lamin B receptor levels by miR-340-5p disrupts chromatin, promotes cell senescence and enhances senolysis

A major stress response influenced by microRNAs (miRNAs) is senescence, a state of indefinite growth arrest triggered by sublethal cell damage. Here, through bioinformatic analysis and experimental validation, we identified miR-340-5p as a novel miRNA that foments cellular senescence. miR-340-5p was highly abundant in diverse senescence models, and miR-340-5p overexpression in proliferating cells rendered them senescent. Among the target mRNAs, miR-340-5p prominently reduced the levels of LBR mRNA, encoding lamin B receptor (LBR). Loss of LBR by ectopic overexpression of miR-340-5p derepressed heterochromatin in lamina-associated domains, promoting the expression of DNA repetitive elements characteristic of senescence. Importantly, overexpressing miR-340-5p enhanced cellular sensitivity to senolytic compounds, while antagonization of miR-340-5p reduced senescent cell markers and engendered resistance to senolytic-induced cell death. We propose that miR-340-5p can be exploited for removing senescent cells to restore tissue homeostasis and mitigate damage by senescent cells in pathologies of human aging.

CoBRA: Containerized Bioinformatics Workflow for Reproducible ChIP/ATAC-seq Analysis

Chromatin immunoprecipitation sequencing (ChIP-seq) and the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) have become essential technologies to effectively measure protein–DNA interactions and chromatin accessibility. However, there is a need for a scalable and reproducible pipeline that incorporates proper normalization between samples, correction of copy number variations, and integration of new downstream analysis tools. Here we present Containerized Bioinformatics workflow for Reproducible ChIP/ATAC-seq Analysis (CoBRA), a modularized computational workflow which quantifies ChIP-seq and ATAC-seq peak regions and performs unsupervised and supervised analyses. CoBRA provides a comprehensive state-of-the-art ChIP-seq and ATAC-seq analysis pipeline that can be used by scientists with limited computational experience. This enables researchers to gain rapid insight into protein–DNA interactions and chromatin accessibility through sample clustering, differential peak calling, motif enrichment, comparison of sites to a reference database, and pathway analysis. CoBRA is publicly available online at https://bitbucket.org/cfce/cobra

Unsupervised tensor decomposition-based method to extract candidate transcription factors as histone modification bookmarks in post-mitotic transcriptional reactivation

The histone group added to a gene sequence must be removed during mitosis to halt transcription during the DNA replication stage of the cell cycle. However, the detailed mechanism of this transcription regulation remains unclear. In particular, it is not realistic to reconstruct all appropriate histone modifications throughout the genome from scratch after mitosis. Thus, it is reasonable to assume that there might be a type of “bookmark” that retains the positions of histone modifications, which can be readily restored after mitosis. We developed a novel computational approach comprising tensor decomposition (TD)-based unsupervised feature extraction (FE) to identify transcription factors (TFs) that bind to genes associated with reactivated histone modifications as candidate histone bookmarks. To the best of our knowledge, this is the first application of TD-based unsupervised FE to the cell division context and phases pertaining to the cell cycle in general. The candidate TFs identified with this approach were functionally related to cell division, suggesting the suitability of this method and the potential of the identified TFs as bookmarks for histone modification during mitosis.

p53 mediates target gene association with nuclear speckles for amplified RNA expression

Nuclear speckles are prominent nuclear bodies that contain proteins and RNA involved in gene expression. Although links between nuclear speckles and gene activation are emerging, the mechanisms regulating association of genes with speckles are unclear. We find that speckle association of p53 target genes is driven by the p53 transcription factor. Focusing on p21, a key p53 target, we demonstrate that speckle association boosts expression by elevating nascent RNA amounts. p53-regulated speckle association did not depend on p53 transactivation functions but required an intact proline-rich domain and direct DNA binding, providing mechanisms within p53 for regulating gene-speckle association. Beyond p21, a substantial subset of p53 targets have p53-regulated speckle association. Strikingly, speckle-associating p53 targets are more robustly activated and occupy a distinct niche of p53 biology compared with non-speckle-associating p53 targets. Together, our findings illuminate regulated speckle association as a mechanism used by a transcription factor to boost gene expression.

Core transcription regulatory circuitry orchestrates corneal epithelial homeostasis

Adult stem cell identity, plasticity, and homeostasis are precisely orchestrated by lineage-restricted epigenetic and transcriptional regulatory networks. Here, by integrating super-enhancer and chromatin accessibility landscapes, we delineate core transcription regulatory circuitries (CRCs) of limbal stem/progenitor cells (LSCs) and find that RUNX1 and SMAD3 are required for maintenance of corneal epithelial identity and homeostasis. RUNX1 or SMAD3 depletion inhibits PAX6 and induces LSCs to differentiate into epidermal-like epithelial cells. RUNX1, PAX6, and SMAD3 (RPS) interact with each other and synergistically establish a CRC to govern the lineage-specific cis-regulatory atlas. Moreover, RUNX1 shapes LSC chromatin architecture via modulating H3K27ac deposition. Disturbance of RPS cooperation results in cell identity switching and dysfunction of the corneal epithelium, which is strongly linked to various human corneal diseases. Our work highlights CRC TF cooperativity for establishment of stem cell identity and lineage commitment, and provides comprehensive regulatory principles for human stratified epithelial homeostasis and pathogenesis.

Corneal epithelium shares similar molecular signatures to other stratified epithelia. Here, the authors map super-enhancers and accessible chromatin in corneal epithelium, identifying a transcription regulatory circuit, including RUNX1, PAX6, and SMAD3, required for corneal epithelial identity and homeostasis.

Novel ChIP-seq simulating program with superior versatility: isChIP

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) is recognized as an extremely powerful tool to study the interaction of numerous transcription factors and other chromatin-associated proteins with DNA. The core problem in the optimization of ChIP-seq protocol and the following computational data analysis is that a 'true' pattern of binding events for a given protein factor is unknown. Computer simulation of the ChIP-seq process based on 'a-priory known binding template' can contribute to a drastically reduce the number of wet lab experiments and finally help achieve radical optimization of the entire processing pipeline. We present a newly developed ChIP-sequencing simulation algorithm implemented in the novel software, in silico ChIP-seq (isChIP). We demonstrate that isChIP closely approximates real ChIP-seq protocols and is able to model data similar to those obtained from experimental sequencing. We validated isChIP using publicly available datasets generated for well-characterized transcription factors Oct4 and Sox2. Although the novel software is compatible with the Illumina protocols by default, it can also successfully perform simulations with a number of alternative sequencing platforms such as Roche454, Ion Torrent and SOLiD as well as model ChIP -Exo. The versatility of isChIP was demonstrated through modelling a wide range of binding events, including those of transcription factors and chromatin modifiers. We also performed a comparative analysis against a few existing ChIP-seq simulators and showed the fundamental superiority of our model. Due to its ability to utilize known binding templates, isChIP can potentially be employed to help investigators choose the most appropriate analytical software through benchmarking of available ChIP-seq programs and optimize the experimental parameters of ChIP-seq protocol. isChIP software is freely available at https://github.com/fnaumenko/isChIP.

Quantitative analysis of ChIP-seq data uncovers dynamic and sustained H3K4me3 and H3K27me3 modulation in cancer cells under hypoxia

Background

A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understanding the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treatment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets.

Results

We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around centromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia-induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells.

Conclusions

In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0090-4) contains supplementary material, which is available to authorized users.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a “gradient effect” where the regulatory factor occupancy levels follow transcription and ii) a “threshold effect” where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.

A comprehensive comparison of tools for differential ChIP-seq analysis

ChIP-seq has become a widely adopted genomic assay in recent years to determine binding sites for transcription factors or enrichments for specific histone modifications. Beside detection of enriched or bound regions, an important question is to determine differences between conditions. While this is a common analysis for gene expression, for which a large number of computational approaches have been validated, the same question for ChIP-seq is particularly challenging owing to the complexity of ChIP-seq data in terms of noisiness and variability. Many different tools have been developed and published in recent years. However, a comprehensive comparison and review of these tools is still missing. Here, we have reviewed 14 tools, which have been developed to determine differential enrichment between two conditions. They differ in their algorithmic setups, and also in the range of applicability. Hence, we have benchmarked these tools on real data sets for transcription factors and histone modifications, as well as on simulated data sets to quantitatively evaluate their performance. Overall, there is a great variety in the type of signal detected by these tools with a surprisingly low level of agreement. Depending on the type of analysis performed, the choice of method will crucially impact the outcome.