The CUT&RUN suspect list of problematic regions of the genome

The scaffolding function of LSD1 controls DNA methylation in mouse ESCs

Lysine-specific histone demethylase 1 (LSD1), which demethylates mono- or di- methylated histone H3 on lysine 4 (H3K4me1/2), is essential for early embryogenesis and development. Here we show that LSD1 is dispensable for mouse embryonic stem cell (ESC) self-renewal but is required for mouse ESC growth and differentiation. Reintroduction of a catalytically-impaired LSD1 (LSD1MUT) recovers the proliferation capability of mouse ESCs, yet the enzymatic activity of LSD1 is essential to ensure proper differentiation. Indeed, increased H3K4me1 in Lsd1 knockout (KO) mouse ESCs does not lead to major changes in global gene expression programs related to stemness. However, ablation of LSD1 but not LSD1MUT results in decreased DNMT1 and UHRF1 proteins coupled to global hypomethylation. We show that both LSD1 and LSD1MUT control protein stability of UHRF1 and DNMT1 through interaction with HDAC1 and the ubiquitin-specific peptidase 7 (USP7), consequently, facilitating the deacetylation and deubiquitination of DNMT1 and UHRF1. Our studies elucidate a mechanism by which LSD1 controls DNA methylation in mouse ESCs, independently of its lysine demethylase activity.

Initial Characterization of WDR5B Reveals a Role in the Proliferation of Retinal Pigment Epithelial Cells

The chromatin-associated protein WDR5 has been widely studied due to its role in histone modification and its potential as a pharmacological target for the treatment of cancer. In humans, the protein with highest sequence homology to WDR5 is encoded by the retrogene WDR5B, which remains unexplored. Here, we used CRISPR-Cas9 genome editing to generate WDR5B knockout and WDR5B-FLAG knock-in cell lines for further characterization. In contrast to WDR5, WDR5B exhibits low expression in pluripotent cells and is upregulated upon neural differentiation. Loss or shRNA depletion of WDR5B impairs cell growth and increases the fraction of non-viable cells in proliferating retinal pigment epithelial (RPE) cultures. CUT&RUN chromatin profiling in RPE and neural progenitors indicates minimal WDR5B enrichment at established WDR5 binding sites. These results suggest that WDR5 and WDR5B exhibit several divergent biological properties despite sharing a high degree of sequence homology.

Identification and Characterization of ATOH7-Regulated Target Genes and Pathways in Human Neuroretinal Development

The proneural transcription factor atonal basic helix-loop-helix transcription factor 7 (ATOH7) is expressed in early progenitors in the developing neuroretina. In vertebrates, this is crucial for the development of retinal ganglion cells (RGCs), as mutant animals show an almost complete absence of RGCs, underdeveloped optic nerves, and aberrations in retinal vessel development. Human mutations are rare and result in autosomal recessive optic nerve hypoplasia (ONH) or severe vascular changes, diagnosed as autosomal recessive persistent hyperplasia of the primary vitreous (PHPVAR). To better understand the role of ATOH7 in neuroretinal development, we created ATOH7 knockout and eGFP-expressing ATOH7 reporter human induced pluripotent stem cells (hiPSCs), which were differentiated into early-stage retinal organoids. Target loci regulated by ATOH7 were identified by Cleavage Under Targets and Release Using Nuclease with sequencing (CUT&RUN-seq) and differential expression by RNA sequencing (RNA-seq) of wildtype and mutant organoid-derived reporter cells. Additionally, single-cell RNA sequencing (scRNA-seq) was performed on whole organoids to identify cell type-specific genes. Mutant organoids displayed substantial deficiency in axon sprouting, reduction in RGCs, and an increase in other cell types. We identified 469 differentially expressed target genes, with an overrepresentation of genes belonging to axon development/guidance and Notch signaling. Taken together, we consolidate the function of human ATOH7 in guiding progenitor competence by inducing RGC-specific genes while inhibiting other cell fates. Furthermore, we highlight candidate genes responsible for ATOH7-associated optic nerve and retinovascular anomalies, which sheds light to potential future therapy targets for related disorders.

The histone deacetylase HDAC1 controls dendritic cell development and anti-tumor immunity

Dendritic cell (DC) progenitors adapt their transcriptional program during development, generating different subsets. How chromatin modifications modulate these processes is unclear. Here, we investigate the impact of histone deacetylation on DCs by genetically deleting histone deacetylase 1 (HDAC1) or HDAC2 in hematopoietic progenitors and CD11c-expressing cells. While HDAC2 is not critical for DC development, HDAC1 deletion impairs pro-pDC and mature pDC generation and affects ESAM+cDC2 differentiation from tDCs and pre-cDC2s, whereas cDC1s are unchanged. HDAC1 knockdown in human hematopoietic cells also impairs cDC2 development, highlighting its crucial role across species. Multi-omics analyses reveal that HDAC1 controls expression, chromatin accessibility, and histone acetylation of the transcription factors IRF4, IRF8, and SPIB required for efficient development of cDC2 subsets. Without HDAC1, DCs switch immunologically, enhancing tumor surveillance through increased cDC1 maturation and interleukin-12 production, driving T helper 1-mediated immunity and CD8+ T cell recruitment. Our study reveals the importance of histone acetylation in DC development and anti-tumor immunity, suggesting DC-targeted therapeutic strategies for immuno-oncology.

Aberrant DNA methylation distorts developmental trajectories in atypical teratoid/rhabdoid tumors

Atypical teratoid/rhabdoid tumors (AT/RTs) are pediatric brain tumors known for their aggressiveness and aberrant but still unresolved epigenetic regulation. To better understand their malignancy, we investigated how AT/RT-specific DNA hypermethylation was associated with gene expression and altered transcription factor binding and how it is linked to upstream regulation. Medulloblastomas, choroid plexus tumors, pluripotent stem cells, and fetal brain were used as references. A part of the genomic regions, which were hypermethylated in AT/RTs similarly as in pluripotent stem cells and demethylated in the fetal brain, were targeted by neural transcriptional regulators. AT/RT-unique DNA hypermethylation was associated with polycomb repressive complex 2 and linked to suppressed genes with a role in neural development and tumorigenesis. Activity of the several NEUROG/NEUROD pioneer factors, which are unable to bind to methylated DNA, was compromised via the suppressed expression or DNA hypermethylation of their target sites, which was also experimentally validated for NEUROD1 in medulloblastomas and AT/RT samples. These results highlight and characterize the role of DNA hypermethylation in AT/RT malignancy and halted neural cell differentiation.

Exhaustive identification of genome-wide binding events of transcriptional regulators

Genome-wide binding assays aspire to map the complete binding pattern of gene regulators. Common practice relies on replication-duplicates or triplicates-and high stringency statistics to favor false negatives over false positives. Here we show that duplicates and triplicates of CUT&RUN are not sufficient to discover the entire activity of transcriptional regulators. We introduce ICEBERG (Increased Capture of Enrichment By Exhaustive Replicate aGgregation), a pipeline that harnesses large numbers of CUT&RUN replicates to discover the full set of binding events and chart the line between false positives and false negatives. We employed ICEBERG to map the full set of H3K4me3-marked regions, the targets of the co-factor β-catenin, and those of the transcription factor TBX3, in human colorectal cancer cells. The ICEBERG datasets allow benchmarking of individual replicates, comparing the performance of peak calling and replication approaches, and expose the arbitrary nature of strategies to identify reproducible peaks. Instead of a static view of genomic targets, ICEBERG establishes a spectrum of detection probabilities across the genome for a given factor, underlying the intrinsic dynamicity of its mechanism of action, and permitting to distinguish frequent from rare regulation events. Finally, ICEBERG discovered instances, undetectable with other approaches, that underlie novel mechanisms of colorectal cancer progression.

Mechanisms of Epigenomic and Functional Convergence Between Glucocorticoid- and IL4-Driven Macrophage Programming

Macrophages adopt distinct phenotypes in response to environmental cues, with type-2 cytokine interleukin-4 promoting a tissue-repair homeostatic state (M2IL4). Glucocorticoids, widely used anti-inflammatory therapeutics, reportedly impart a similar phenotype (M2GC), but how such disparate pathways may functionally converge is unknown. We show using integrative functional genomics that M2IL4 and M2GC transcriptomes share a striking overlap mirrored by a shift in chromatin landscape in both common and signal-specific gene subsets. This core homeostatic program is enacted by transcriptional effectors KLF4 and the GC receptor, whose genome-wide occupancy and actions are integrated in a stimulus-specific manner by the nuclear receptor cofactor GRIP1. Indeed, many of the M2IL4:M2GC-shared transcriptomic changes were GRIP1-dependent. Consistently, GRIP1 loss attenuated phagocytic activity of both populations in vitro and macrophage tissue-repair properties in the murine colitis model in vivo. These findings provide a mechanistic framework for homeostatic macrophage programming by distinct signals, to better inform anti-inflammatory drug design.

The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

The histone chaperone ANP32B regulates chromatin incorporation of the atypical human histone variant macroH2A

All vertebrate genomes encode for three large histone H2A variants that have an additional metabolite-binding globular macrodomain module, macroH2A. MacroH2A variants impact heterochromatin organization and transcription regulation and establish a barrier for cellular reprogramming. However, the mechanisms of how macroH2A is incorporated into chromatin and the identity of any chaperones required for histone deposition remain elusive. Here, we develop a split-GFP-based assay for chromatin incorporation and use it to conduct a genome-wide mutagenesis screen in haploid human cells to identify proteins that regulate macroH2A dynamics. We show that the histone chaperone ANP32B is a regulator of macroH2A deposition. ANP32B associates with macroH2A in cells and in vitro binds to histones with low nanomolar affinity. In vitro nucleosome assembly assays show that ANP32B stimulates deposition of macroH2A-H2B and not of H2A-H2B onto tetrasomes. In cells, depletion of ANP32B strongly affects global macroH2A chromatin incorporation, revealing ANP32B as a macroH2A histone chaperone.

excluderanges: exclusion sets for T2T-CHM13, GRCm39, and other genome assemblies

SUMMARY

Exclusion regions are sections of reference genomes with abnormal pileups of short sequencing reads. Removing reads overlapping them improves biological signal, and these benefits are most pronounced in differential analysis settings. Several labs created exclusion region sets, available primarily through ENCODE and Github. However, the variety of exclusion sets creates uncertainty which sets to use. Furthermore, gap regions (e.g. centromeres, telomeres, short arms) create additional considerations in generating exclusion sets. We generated exclusion sets for the latest human T2T-CHM13 and mouse GRCm39 genomes and systematically assembled and annotated these and other sets in the excluderanges R/Bioconductor data package, also accessible via the BEDbase.org API. The package provides unified access to 82 GenomicRanges objects covering six organisms, multiple genome assemblies, and types of exclusion regions. For human hg38 genome assembly, we recommend hg38.Kundaje.GRCh38_unified_blacklist as the most well-curated and annotated, and sets generated by the Blacklist tool for other organisms.

AVAILABILITY AND IMPLEMENTATION

https://bioconductor.org/packages/excluderanges/. Package website: https://dozmorovlab.github.io/excluderanges/.