CHIPIN: ChIP-seq inter-sample normalization based on signal invariance across transcriptionally constant genes

Multiplexed chromatin immunoprecipitation sequencing for quantitative study of histone modifications and chromatin factors

ChIP-seq is a widely used technique for studying histone post-translational modifications and DNA-binding proteins. DNA fragments associated with a specific protein or histone modification epitope are captured by using antibodies, sequenced and mapped to a reference genome. Albeit versatile and popular, performing many parallel ChIP-seq experiments to compare different conditions, replicates and epitopes is laborious, is prone to experimental variation and does not allow quantitative comparisons unless adequate spike-in chromatin is included. We present a detailed protocol for performing and analyzing a multiplexed quantitative chromatin immunoprecipitation-sequencing experiment (MINUTE-ChIP), in which multiple samples are profiled against multiple epitopes in a single workflow. Multiplexing not only dramatically increases the throughput of ChIP-seq experiments (e.g., profiling 12 samples against multiple histone modifications or DNA-binding proteins in a single experiment), but also enables accurate quantitative comparisons. The protocol consists of four parts: sample preparation (i.e., lysis, chromatin fragmentation and barcoding of native or formaldehyde-fixed material), pooling and splitting of the barcoded chromatin into parallel immunoprecipitation reactions, preparation of next-generation sequencing libraries from input and immunoprecipitated DNA and data analysis using our dedicated analysis pipeline. This pipeline autonomously generates quantitatively scaled ChIP-seq tracks for downstream analysis and visualization, alongside necessary quality control indicators. The entire workflow requires basic knowledge in molecular biology and bioinformatics and can be completed in 1 week. MINUTE-ChIP empowers biologists to perform every ChIP-seq experiment with an appropriate number of replicates and control conditions, delivering more statistically robust, exquisitely quantitative and biologically meaningful results.

Epigenetic characterization of adult rhesus monkey spermatogonial stem cells identifies key regulators of stem cell homeostasis

Spermatogonial stem cells (SSCs) play crucial roles in the preservation of male fertility. However, successful ex vivo expansion of authentic human SSCs remains elusive due to the inadequate understanding of SSC homeostasis regulation. Using rhesus monkeys (Macaca mulatta) as a representative model, we characterized SSCs and progenitor subsets through single-cell RNA sequencing using a cell-specific network approach. We also profiled chromatin status and major histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K27me3 and H3K9me3), and subsequently mapped promoters and active enhancers in TSPAN33+ putative SSCs. Comparing the epigenetic changes between fresh TSPAN33+ cells and cultured TSPAN33+ cells (resembling progenitors), we identified the regulatory elements with higher activity in SSCs, and the potential transcription factors and signaling pathways implicated in SSC regulation. Specifically, TGF-β signaling is activated in monkey putative SSCs. We provided evidence supporting its role in promoting self-renewal of monkey SSCs in culture. Overall, this study outlines the epigenetic landscapes of monkey SSCs and provides clues for optimization of the culture condition for primate SSCs expansion.

Dissecting transposable elements and endogenous retroviruses upregulation by HDAC inhibitors in leiomyosarcoma cells: Implications for the interferon response

Transposable elements (TEs) are of interest as immunomodulators for cancer therapies. TEs can fold into dsRNAs that trigger the interferon response. Here, we investigated the effect of different HDAC inhibitors (HDACIs) on the expression of TEs in leiomyosarcoma cells. Our data show that endogenous retroviruses (ERVs), especially ERV1 elements, are upregulated after treatment with HDAC1/2/3-specific inhibitors. Surprisingly, the interferon response was not activated. We observed an increase in A-to-I editing of upregulated ERV1. This could have an impact on the stability of dsRNAs and the activation of the interferon response. We also found that H3K27ac levels are increased in the LTR12 subfamilies, which could be regulatory elements controlling the expression of proapoptotic genes such as TNFRSF10B. In summary, we provide a detailed characterization of TEs modulation in response to HDACIs and suggest the use of HDACIs in combination with ADAR inhibitors to induce cell death and support immunotherapy in cancer.

TCR-independent CD137 (4-1BB) signaling promotes CD8+-exhausted T cell proliferation and terminal differentiation

CD137 (4-1BB)-activating receptor represents a promising cancer immunotherapeutic target. Yet, the cellular program driven by CD137 and its role in cancer immune surveillance remain unresolved. Using T cell-specific deletion and agonist antibodies, we found that CD137 modulates tumor infiltration of CD8+-exhausted T (Tex) cells expressing PD1, Lag-3, and Tim-3 inhibitory receptors. T cell-intrinsic, TCR-independent CD137 signaling stimulated the proliferation and the terminal differentiation of Tex precursor cells through a mechanism involving the RelA and cRel canonical NF-κB subunits and Tox-dependent chromatin remodeling. While Tex cell accumulation induced by prophylactic CD137 agonists favored tumor growth, anti-PD1 efficacy was improved with subsequent CD137 stimulation in pre-clinical mouse models. Better understanding of T cell exhaustion has crucial implications for the treatment of cancer and infectious diseases. Our results identify CD137 as a critical regulator of Tex cell expansion and differentiation that holds potential for broad therapeutic applications.

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Cell-type specific profiling of histone post-translational modifications in the adult mouse striatum

Epigenetic gene regulation in the heterogeneous brain remains challenging to decipher with current strategies. Bulk tissue analysis from pooled subjects reflects the average of cell-type specific changes across cell-types and individuals, which obscures causal relationships between epigenetic modifications, regulation of gene expression, and complex pathology. To address these limitations, we optimized a hybrid protocol, ICuRuS, for the isolation of nuclei tagged in specific cell-types and histone post translational modification profiling from the striatum of a single mouse. We combined affinity-based isolation of the medium spiny neuron subtypes, Adenosine 2a Receptor or Dopamine Receptor D1, with cleavage of histone-DNA complexes using an antibody-targeted micrococcal nuclease to release DNA complexes for paired end sequencing. Unlike fluorescence activated cell sorting paired with chromatin immunoprecipitation, ICuRuS allowed for robust epigenetic profiling at cell-type specific resolution. Our analysis provides a framework to understand combinatorial relationships between neuronal-subtype-specific epigenetic modifications and gene expression.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

HDAC1 and PRC2 mediate combinatorial control in SPI1/PU.1-dependent gene repression in murine erythroleukaemia

Although originally described as transcriptional activator, SPI1/PU.1, a major player in haematopoiesis whose alterations are associated with haematological malignancies, has the ability to repress transcription. Here, we investigated the mechanisms underlying gene repression in the erythroid lineage, in which SPI1 exerts an oncogenic function by blocking differentiation. We show that SPI1 represses genes by binding active enhancers that are located in intergenic or gene body regions. HDAC1 acts as a cooperative mediator of SPI1-induced transcriptional repression by deacetylating SPI1-bound enhancers in a subset of genes, including those involved in erythroid differentiation. Enhancer deacetylation impacts on promoter acetylation, chromatin accessibility and RNA pol II occupancy. In addition to the activities of HDAC1, polycomb repressive complex 2 (PRC2) reinforces gene repression by depositing H3K27me3 at promoter sequences when SPI1 is located at enhancer sequences. Moreover, our study identified a synergistic relationship between PRC2 and HDAC1 complexes in mediating the transcriptional repression activity of SPI1, ultimately inducing synergistic adverse effects on leukaemic cell survival. Our results highlight the importance of the mechanism underlying transcriptional repression in leukemic cells, involving complex functional connections between SPI1 and the epigenetic regulators PRC2 and HDAC1.

Suppression of p53 response by targeting p53-Mediator binding with a stapled peptide

DNA-binding transcription factors (TFs) remain challenging to target with molecular probes. Many TFs function in part through interaction with Mediator, a 26-subunit complex that controls RNA polymerase II activity genome-wide. We sought to block p53 function by disrupting the p53-Mediator interaction. Through rational design and activity-based screening, we characterize a stapled peptide, with functional mimics of both p53 activation domains, that blocks p53-Mediator binding and selectively inhibits p53-dependent transcription in human cells; importantly, this "bivalent" peptide has negligible impact, genome-wide, on non-p53 target genes. Our proof-of-concept strategy circumvents the TF entirely and targets the TF-Mediator interface instead, with desired functional outcomes (i.e., selective inhibition of p53 activation). Furthermore, these results demonstrate that TF activation domains represent viable starting points for Mediator-targeting molecular probes, as an alternative to large compound libraries. Different TFs bind Mediator through different subunits, suggesting this strategy could be broadly applied to selectively alter gene expression programs.