PATTY corrects open-chromatin bias for improved bulk and single-cell CUT&Tag profiling

Precise profiling of epigenomes is essential for better understanding chromatin biology and gene regulation. Cleavage Under Targets & Tagmentation (CUT&Tag) is an efficient epigenomic profiling technique that can be performed on a low number of cells and at the single-cell level. With its growing adoption, CUT&Tag datasets spanning diverse biological systems are rapidly accumulating in the field. CUT&Tag assays use the hyperactive transposase Tn5 for DNA tagmentation. Tn5’s preference toward accessible chromatin alters CUT&Tag sequence read distributions in the genome and introduces open chromatin bias that can confound downstream analysis, an issue more substantial in sparse single-cell data. We show that open chromatin bias extensively exists in published CUT&Tag datasets, including those generated with recently optimized high-salt protocols. To address this challenge, we present PATTY (Propensity Analyzer for Tn5 Transposase Yielded bias), a comprehensive computational method that corrects open chromatin bias in CUT&Tag data by leveraging accompanying ATAC-seq. By integrating transcriptomic and epigenomic data using machine learning and integrative modeling, we demonstrate that PATTY enables accurate and robust detection of occupancy sites for both active and repressive histone modifications, including H3K27ac, H3K27me3, and H3K9me3, with experimental validation. We further develop a single-cell CUT&Tag analysis framework built on PATTY and show improved cell clustering when using bias-corrected single-cell CUT&Tag data compared to using uncorrected data. Beyond CUT&Tag, PATTY sets a foundation for further development of bias correction methods for improving data analysis for all Tn5-based high-throughput assays.

Small cell prostate cancer (SCPC) is an androgen receptor (AR) negative variant that often emerges in the castration-resistant progression of typical AR-positive prostate adenocarcinomas. Its presence predicts for an aggressive clinical course with poor response to hormonal therapies and a short median survival. Evidence suggests that the emergence of (SCPC) occurs through cellular transdifferentiation. Therefore, we hypothesized that loss of AR gene transcription may depend on epigenetic events. Given the central role of AR in prostate cancer, we reasoned that identifying the mechanisms involved in its downregulation could improve our understanding of how the lethal SCPC subtype arises. Here we investigated the DNA methylation and histone modifications of AR in a group of six prostate tumor xenografts developed from men with CRPC (two AR-positive [MDA PCa 170.4 and MDA PCa 180.30] and four AR-negative SCPC [MDA PCa 144.13, MDA PCa 144.4, MDA PCa 155 and MDA PCa 146.10]) and three established prostate cancer cell lines (one AR-positive [LNCaP] and two AR-negative [PC3 and DU145]). We first evaluated the methylation status of selected CpG sites in the AR promoter-associated CpG island using bisulfite-sequencing followed by pyrosequencing analysis. Except for the cancer cell lines PC-3 and DU145, all other cell lines and xenografts were unmethylated regardless of the AR expression status. Next, we used chromatin immunoprecipitation (ChIP) to evaluate marking of the AR promoter and a putative proximal enhancer by active (H3K4me3 and H3K9ac) and repressive (H3K9me2 and H3K27me3) histone modifications. Marking of the constitutively expressed genes (ACTB and GAPDH) and a repressed gene (HBB) served as control for these experiments. As expected, the only samples with AR marking by H3K4me3 and H3K9ac were the two AR-positive xenografts and LNCaP. Marking by repressive histone modifications was more variable: H3K27me3 was a universal finding in AR-negative samples except for Du-145, and was accompanied by H3K9me2 in two out of four xenografts. H3K27me3 and H3K9me2 also marked the putative enhancer region when AR was repressed. Altogether, these data indicate that promoter DNA methylation of AR is an infrequent finding in prostate cancer, while polycomb protein-based silencing is nearly universal in AR-negative variants. Also, it shows that H3K9me2 promotes reinforcement of silencing. Thus, we propose that reactivation of AR may by achieved in selected cases by inhibiting H3K27me3 alone, while concomitant inhibition of both polycomb and SET domain proteins may be necessary for double-positive H3K27/K9 methylation cases. Given the purported differentiating and tumor suppressor effects of AR, its reactivation in SCPC may be of therapeutic benefit. Citation Format: Brittany N. Kleb, Marcos RH Estecio, Guanglin Wu, Jing-Fang Lu, Christopher J. Logothetis, Sankar Maity, Ana M. Aparicio. Mechanism of androgen receptor (AR) silencing in small cell prostate carcinoma. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2986. doi:10.1158/1538-7445.AM2013-2986

The accurate detection of transcription factor (TF) binding sites and histone modifications (HM) on the genome-wide scale is essential for studying functional epigenetics and gene regulation. Cleavage Under Targets & Tagmentation (CUT&Tag) is a low-cost and easy-to-implement epigenomic profiling method that can be performed on a low number of cells and on the single-cell level. CUT&Tag experiments use the hyperactive transposase Tn5 for tagmentation. We find that Tn5 is subject to intrinsic sequence insertion bias (intrinsic bias). Additionally, the preference of Tn5 insertion toward accessible chromatin also affects the distribution of CUT&Tag reads (open chromatin bias). Both types of biases can significantly confound bulk and single-cell CUT&Tag data analysis, which requires careful assessment and new analytical methods. To address this challenge, we present PATTY (Propensity Analyzer for Tn5 Transposase Yielded bias), a computational method for systematic characterization and correction of biases in CUT&Tag data. We demonstrate that histone modification signals (H3K27me3 and H3K27ac) detected from CUT&Tag data after bias correction using PATTY are better associated with orthogonal biological features such as gene expression. PATTY-corrected single-cell CUT&Tag signals for histone modification can better cluster cells into their true cell types. This new computational method can improve both bulk and single-cell CUT&Tag data analysis. Citation Format: Sheng'en Shawn Hu, Qingying Chen, Megan Grieco, Chongzhi Zang. PATTY: A bias estimation and correction model for bulk and single-cell CUT&Tag data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7435.

Poised with purpose: Cell plasticity enhances tumorigenicity

Although a growing body of evidence indicates that phenotypic plasticity exhibited by glioblastoma cells plays a central role in tumor development and post-therapy recurrence, the master drivers of their aggressiveness remain elusive. Here we mapped the changes in active (H3K4me3) and repressive (H3K27me3) histone modifications accompanying the repression of glioblastoma stem-like cells tumorigenicity. Genes with changing histone marks delineated a network of transcription factors related to cancerous behavior, stem state, and neural development, highlighting a previously unsuspected association between repression of ARNT2 and loss of cell tumorigenicity. Immunohistochemistry confirmed ARNT2 expression in cell sub-populations within proliferative zones of patients’ glioblastoma. Decreased ARNT2 expression was consistently observed in non-tumorigenic glioblastoma cells, compared to tumorigenic cells. Moreover, ARNT2 expression correlated with a tumorigenic molecular signature at both the tissue level within the tumor core and at the single cell level in the patients’ tumors. We found that ARNT2 knockdown decreased the expression of SOX9, POU3F2 and OLIG2, transcription factors implicated in glioblastoma cell tumorigenicity, and repressed glioblastoma stem-like cell tumorigenic properties in vivo. Our results reveal ARNT2 as a pivotal component of the glioblastoma cell tumorigenic signature, located at a node of a transcription factor network controlling glioblastoma cell aggressiveness.

Pre-clinical models have identified diverse epigenetic alterations in prostate cancer that associate with treatment resistance and may also be therapeutic targets. However, these assays typically require thousands of cells and are difficult to scale for patient tumor biopsies. This report seeks to adapt the Cleavage Under Targets and Tagmentation (CUT&Tag) assay for use on a microfluidic technology to eliminate cell loss associated with complex, multi-step epigenetic assays to facilitate analysis of low input primary tumor samples. Exclusive Liquid Repellant (ELR) platforms pass magnetic bead bound analytes through immiscible phases to separate and extract analytes of interest while eliminating traditional centrifugation and wash steps that lead to sample loss. This platform enables adaptation of complex molecular assays such as CUT&Tag for lossless sample manipulation. We focus on the analysis of histone modifications in rare cells utilizing antibodies for H3K27 acetylation and methylation. Comparison of prostate cancer cell lines with LNCaP, V16D, 49F, 42D, and H660s with serial dilutions was performed to evaluate the sensitivity and specificity of this integrated assay. The CUT&Tag assay was performed in standard fashion on a serial dilution of 50-100,000 cells from LNCaP and V16D calls and showed dropout of recoverable DNA at 50 cells. Sequencing of DNA at these dilutions confirmed identification of differential histone modifications consistent with previously published data. CUT&Tag was adapted to ELR through the serial patterning of each buffer in distinct wells connected by an oil interface. In order to optimize the assay for reproducibility and increased sensitivity, we optimized antibody binding/washing to minimize the protein A-Transposase inhibition. This optimization showed increased antibody binding by 36% compared to protocol defined concentrations. We have now successfully performed this assay with a sensitivity of ~100 cells with recoverable DNA for sequencing to identify heterogenous histone patterns in each cell line. CUT&Tag has improved sensitivity to detect epigenetic alterations that associate with aggressive prostate cancer. The integration of CUT&Tag with ELR improves both sensitivity and reproducibility for epigenetic analysis at low input cell numbers. This assay is now being utilized to analyze patient tumor samples and prostate cancer heterogeneity. Citation Format: Zachary J. Kauffman, Kevin Koesser, Kyle T. Helzer, Jamie M. Sperger, Marina N. Sharifi, Chao Li, Duane S. Juang, Cole S. Gilsdorf, Shuang G. Zhao, David J. Beebe, Joshua M. Lang. Microscale analysis of histone modifications in rare cell populations in prostate cancer [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr A057.

In mammals, genome-wide chromatin maps and immunofluorescence studies show that broad domains of repressive histone modifications are present on pericentromeric and telomeric repeats and on the inactive X chromosome. However, only a few autosomal loci such as silent Hox gene clusters have been shown to lie in broad domains of repressive histone modifications. Here we present a ChIP-chip analysis of the repressive H3K27me3 histone modification along chr 17 in mouse embryonic fibroblast cells using an algorithm named broad local enrichments (BLOCs), which allows the identification of broad regions of histone modifications. Our results, confirmed by BLOC analysis of a whole genome ChIP-seq data set, show that the majority of H3K27me3 modifications form BLOCs rather than focal peaks. H3K27me3 BLOCs modify silent genes of all types, plus flanking intergenic regions and their distribution indicates a negative correlation between H3K27me3 and transcription. However, we also found that some nontranscribed gene-poor regions lack H3K27me3. We therefore performed a low-resolution analysis of whole mouse chr 17, which revealed that H3K27me3 is enriched in mega-base-pair-sized domains that are also enriched for genes, short interspersed elements (SINEs) and active histone modifications. These genic H3K27me3 domains alternate with similar-sized gene-poor domains. These are deficient in active histone modifications, as well as H3K27me3, but are enriched for long interspersed elements (LINEs) and long-terminal repeat (LTR) transposons and H3K9me3 and H4K20me3. Thus, an autosome can be seen to contain alternating chromatin bands that predominantly separate genes from one retrotransposon class, which could offer unique domains for the specific regulation of genes or the silencing of autonomous retrotransposons.

Mycophenolic acid (MPA) is an important immunosuppressive drug (ISD) prescribed to prevent graft rejection in the organ transplanted patients. However, its use is sporadically linked to the leak flux diarrhea and other gastrointestinal (GI) disturbances in around 75% patients through an unknown mechanism. Recently, we reported increased tight junctions (TJs) permeability through p38MAPK mediated MLCK pathway in MPA treated Caco-2 cells monolayer. On the other hand, certain drugs are known to alter gene expression via epigenetic changes causing un-expected side effects. Whether MPA-induced GI complication, especially diarrhea, has any link with alteration at epigenetic level is unknown. In the present study, we investigated epigenetic changes and possible involvement of midkine dependent PI3K pathway in alteration of TJs under MPA treatment. Caco-2 cells were grown in monolayers to develop the TJs and were treated for 72 hours with MPA or MPA+PI3K inhibitors (/LY/AMG) or MPA+midkine inhibitors (iMDK) and or DMSO. ChIP-O-Proteomics was performed to identify proteins associated with active (H3K4me3) and repressive (H3K27me3) chromatin histone modifications after MPA treatment. The role of appropriate candidate(s) in the context of epithelial monolayer integrity maintenance was analyzed. Caco-2 monolayer integrity was assessed by TEER and FITC-dextran value. We identified a total of 333 and 306 proteins associated with active and repressive histone modification marks, respectively. However, a number of novel proteins have also been identified including midkine. Midkine is a growth factor which is involved in the regulation of various diseases of inflammatory background. Our further analysis identified that midkine functions as the upstream regulator and alters the expression of various genes involved in the TJs maintenance. Epigenetic analysis shows that MPA treatment increases the global histone acetylation levels as well as the enrichment for transcriptional active histone modification mark (H3K4me3) not only at the promoter region of midkine, but also at the promoter regions of midkine effector genes such as PI3K, cdx-2, claudin-2, p38MAPK, ATF-2, MLCK and MLC-2. In contrast, the promoter region of occludin and claudin-1 was enriched for transcriptional repressive histone modification mark (H3K27me3) after MPA treatment. In line with the chromatin status, MPA treatment increased the expression of midkine, PI3K, Cdx-2, claudin-2, p38MAPK, ATF-2, MLCK, and MLC-2 both at transcriptional and translational level, while occludin and claudin-1 expression were negatively influenced. Interestingly, the MPA induced gene expression changes and functional properties of Caco-2 cells could be preserved by the inhibition of midkine mediated PI3K pathway using a chemical inhibitor (PI3K inhibitors or iMDK inhibitors). Collectively, our results highlight that MPA disrupts the structure of TJs via midkine dependent activation of PI3K pathway that results in decreased integrity of Caco-2 monolayer. These results led us to suggest that midkine-mediated lose integrity of epithelial monolayer could be the possible cause of GI disturbance (barrier dysfunction) in the intestine, leading to leaky style diarrhea observed in the organ-transplanted patients treated with MPA.

Mycophenolic acid (MPA) is prescribed to maintain allografts in organ-transplanted patients. However, gastrointestinal (GI) complications, particularly diarrhea, are frequently observed as a side effect following MPA therapy. We recently reported that MPA altered the tight junction (TJ)-mediated barrier function in a Caco-2 cell monolayer model system. This study investigates whether MPA induces epigenetic changes which lead to GI complications, especially diarrhea. Methods: We employed a Chromatin Immunoprecipitation-O-Proteomics (ChIP-O-Proteomics) approach to identify proteins associated with active (H3K4me3) as well as repressive (H3K27me3) chromatin histone modifications in MPA-treated cells, and further characterized the role of midkine, a H3K4me3-associated protein, in the context of epithelial monolayer permeability. Results: We identified a total of 333 and 306 proteins associated with active and repressive histone modification marks, respectively. Among them, 241 proteins were common both in active and repressive chromatin, 92 proteins were associated exclusively with the active histone modification mark, while 65 proteins remained specific to repressive chromatin. Our results show that 45 proteins which bind to the active and seven proteins which bind to the repressive chromatin region exhibited significantly altered abundance in MPA-treated cells as compared to DMSO control cells. A number of novel proteins whose function is not known in bowel barrier regulation were among the identified proteins, including midkine. Our functional integrity assays on the Caco-2 cell monolayer showed that the inhibition of midkine expression prior to MPA treatment could completely block the MPA-mediated increase in barrier permeability. Conclusions: The ChIP-O-Proteomics approach delivered a number of novel proteins with potential implications in MPA toxicity. Consequently, it can be proposed that midkine inhibition could be a potent therapeutic approach to prevent the MPA-mediated increase in TJ permeability and leak flux diarrhea in organ transplant patients.

Epigenetic modifications and small RNAs play an important role in gene regulation. Here, we discuss results of our Solexa/Illumina 1G sequencing-based survey of DNA methylation, activating and repressive histone modifications, small RNAs and mRNA in the maize genome. We analyze tissue-specific epigenetic patterns, discuss antagonistic relationships between repressive epigenetic marks and highlight synergistic relationships between activating histone modifications. We discuss our observation that small RNAs show a tissue-specific distribution in maize. Whereas 24-nucleotide long small interfering RNAs (siRNAs) accumulated preferentially in shoots, 21-nucleotide long micro RNAs (miRNAs) were the most abundant group in roots, which follows the transcript level of mop1. Furthermore, we discuss the possibility that a novel class of 22-nucleotide siRNAs might originate from long double-stranded RNAs in an RNA-dependent RNA polymerase (RdRP)-independent manner. This supports the intriguing possibility that maize possesses at least two distinct pathways to generate siRNAs, one of which relies on RdRP and a second one that might be RdRP-independent.

Histone modifications (HM) and transcription factors (TF) are aberrantly regulated in the pathogenesis of humans and are a major class of cancer cell dependencies. Precise detection of HM and TF binding genome-wide is essential for a better understanding of cancer epigenetics. Cleavage Under Targets & Tagmentation (CUT&Tag) is an easy and low-cost epigenomic profiling method that can be performed on a low number of cells and on the single-cell level. A large number of CUT&Tag datasets have been generated in cancer samples, providing a valuable resource. Unbiased analysis of CUT&Tag data is important for cancer epigenetics research. CUT&Tag experiments rely on the hyperactive transposase Tn5 for tagmentation. We performed integrative analysis of CUT&Tag and other high-throughput data, and found that Tn5 is subject to intrinsic sequence insertion biases (intrinsic biases). Preference of Tn5 captured reads toward chromatin accessibility regions (open chromatin biases) also confound the distribution of CUT&Tag reads. Both sources of biases confound the analysis of CUT&Tag data, especially for factors that are not expected to associate with accessible chromatin. High sparsity of single-cell data makes the effect of both biases more substantial compared to bulk data. We present a computational model for accurate characterization and quantitative estimation of biases in CUT&Tag data. This model paved the way for further development of computational tools for improving both bulk and single-cell CUT&Tag data analysis. Citation Format: Sheng'en Hu, Qingying Chen, Chongzhi Zang. Accurate estimation of open chromatin and intrinsic biases in bulk and single-cell CUT&Tag data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 871.

Posttranslational histone modifications serve critical roles in gene regulation by determining the functional status of chromatin. Histone‐modifying enzymes often work in large multiprotein complexes. A paper in this issue of The EMBO Journal describes a new chromatin‐modifying complex called PEAT that acts via histone deacetylation. The PEAT complex is involved in heterochromatin formation and gene repression but also appears to have a locus‐specific activating role, possibly through promoting histone acetylation.

During central nervous system development, spatiotemporal gene expression programs mediate specific lineage decisions to generate neuronal and glial cell types from neural stem cells (NSCs). However, little is known about the epigenetic landscape underlying these highly complex developmental events. Here, we perform ChIP-seq on distinct subtypes of Drosophila FACS-purified NSCs and their differentiated progeny to dissect the epigenetic changes accompanying the major lineage decisions in vivo By analyzing active and repressive histone modifications, we show that stem cell identity genes are silenced during differentiation by loss of their activating marks and not via repressive histone modifications. Our analysis also uncovers a new set of genes specifically required for altering lineage patterns in type II neuroblasts (NBs), one of the two main Drosophila NSC identities. Finally, we demonstrate that this subtype specification in NBs, unlike NSC differentiation, requires Polycomb-group-mediated repression.

A unique property of the mammalian embryo is that stem cells can be derived from its early tissue lineages. These lineages will give rise to the fetus as well as essential extraembryonic tissues. Understanding how chromatin regulation participates in establishment of these lineages in the embryo and their derived stem cells provides insight that will critically inform our understanding of embryogenesis and stem cell biology. Here, we compare the genomewide location of active and repressive histone modifications in embryonic stem cells, trophoblast stem cells, and extraembryonic endoderm stem cells from the mouse. Our results show that the active modification H3K4me3 has a similar role in the three stem cell types, but the repressive modification H3K27me3 varies in abundance and genomewide distribution. Thus, alternative mechanisms mediate transcriptional repression in stem cells from the embryo. In addition, using carrier chromatin immunoprecipitation we show that bivalent histone domains seen in embryonic stem cells exist in pluripotent cells of the early embryo. However, the epigenetic status of extraembryonic progenitor cells in the embryo did not entirely reflect the extraembryonic stem cell lines. These studies indicate that histone modification mechanisms may differ between early embryo lineages and emphasize the importance of examining in vivo and in vitro progenitor cells.

Term human myometrial expression of progesterone receptor (PR)-A is increased relative to PR-B, and as PR-A is a repressor of progesterone action mediated through PR-B, this increase may mediate the withdrawal of progesterone action and precipitate the onset of labour. PR-A and PR-B expression is regulated by two separate promoters of the PR gene. We hypothesized that epigenetic histone modifications at the two promoters contribute to the labour-associated regulation of PR-A and PR-B expression in term myometrium. PR total, PR-B and PR-A mRNA levels were determined using quantitative real-time PCR, and chromatin immunoprecipitation was used to determine the levels of activating and repressive histone modifications at the PR-A and PR-B promoters in human myometrial samples not in labour (n = 4) and in labour (n = 4). Chromatin extracts were immunoprecipitated with antibodies against activating (histone H3 and H4 acetylation and histone H3 lysine 4 trimethylation), and repressive (histone H3 lysine 9 trimethylation, histone H3 lysine 27 trimethylation and asymmetrical histone H3 arginine 2 dimethylation) histone modifications. PR-A mRNA levels increased during labour, while PR-B mRNA levels remained constant resulting in an increase of PR-A/PR-B mRNA ratio, as expected. Regardless of labour status, significantly higher levels of the activating histone modifications were found at the PR-A promoter compared with the PR-B promoter (P <0.001). H3K4me3 increased significantly at both promoters with labour onset (P =0.001). Low levels of the repressive histone modifications were also present at both promoters, with no labour-associated changes observed. Our data indicate that the PR-A promoter is epigenetically marked for activation in term myometrium more extensively than the PR-B promoter, and that labour is associated with an increase in H3K4me3 activating modification, consistent with the previously described increase in PR protein at this time.