Peak calling by Sparse Enrichment Analysis for CUT&RUN chromatin profiling

Suppression of GATA3 promotes epithelial-mesenchymal transition and simultaneous cellular senescence in human extravillous trophoblasts

The regulatory mechanism of the transcription factor GATA3 in the differentiation and maturation process of extravillous trophoblasts (EVT) in early pregnancy placenta, as well as its relevance to the occurrence of pregnancy disorders, remains poorly understood. This study leveraged single-cell RNA sequencing data from placental organoid models and placental tissue to explore the dynamic changes in GATA3 expression during EVT maturation. The expression pattern exhibited an initial upregulation followed by subsequent downregulation, with aberrant GATA3 localization observed in cases of recurrent miscarriage (RM). By identifying global targets regulated by GATA3 in primary placental EVT cells, JEG3, and HTR8/SVneo cell lines, this study offered insights into its regulatory mechanisms across different EVT cell models. Shared regulatory targets among these cell types and activation of trophoblast cell marker genes emphasized the importance of GATA3 in EVT differentiation and maturation. Knockdown of GATA3 in JEG3 cells led to repression of GATA3-induced epithelial-mesenchymal transition (EMT), as evidenced by changes in marker gene expression levels and enhanced migration ability. Additionally, interference with GATA3 accelerated cellular senescence, as indicated by reduced proliferation rates and increased activity levels for senescence-associated β-galactosidase enzyme, along with elevated expression levels for senescence-associated genes. This study provides comprehensive insights into the dual role of GATA3 in regulating EMT and cellular senescence during EVT differentiation, shedding light on the dynamic changes in GATA3 expression in normal and pathological placental conditions.

IκBα controls dormancy in hematopoietic stem cells via retinoic acid during embryonic development

Recent findings suggest that Hematopoietic Stem Cells (HSC) and progenitors arise simultaneously and independently of each other already in the embryonic aorta-gonad mesonephros region, but it is still unknown how their different features are established. Here, we uncover IκBα (Nfkbia, the inhibitor of NF-κB) as a critical regulator of HSC proliferation throughout development. IκBα balances retinoic acid signaling levels together with the epigenetic silencer, PRC2, specifically in HSCs. Loss of IκBα decreases proliferation of HSC and induces a dormancy related gene expression signature instead. Also, IκBα deficient HSCs respond with superior activation to in vitro culture and in serial transplantation. At the molecular level, chromatin regions harboring binding motifs for retinoic acid signaling are hypo-methylated for the PRC2 dependent H3K27me3 mark in IκBα deficient HSCs. Overall, we show that the proliferation index in the developing HSCs is regulated by a IκBα-PRC2 axis, which controls retinoic acid signaling.

Bhlhe40 Promotes CD4+ T Helper 1 Cell and Suppresses T Follicular Helper Cell Differentiation during Viral Infection

In response to acute infection, naive CD4+ T cells primarily differentiate into T helper 1 (Th1) or T follicular helper (Tfh) cells that play critical roles in orchestrating cellular or humoral arms of immunity, respectively. However, despite the well established role of T-bet and BCL-6 in driving Th1 and Tfh cell lineage commitment, respectively, whether additional transcriptional circuits also underlie the fate bifurcation of Th1 and Tfh cell subsets is not fully understood. In this article, we study how the transcriptional regulator Bhlhe40 dictates the Th1/Tfh differentiation axis in mice. CD4+ T cell-specific deletion of Bhlhe40 abrogates Th1 but augments Tfh differentiation. We also assessed an increase in germinal center B cells and Ab production, suggesting that deletion of Bhlhe40 in CD4+ T cells not only alters Tfh differentiation but also their capacity to provide help to B cells. To identify molecular mechanisms by which Bhlhe40 regulates Th1 versus Tfh lineage choice, we first performed epigenetic profiling in the virus specific Th1 and Tfh cells following LCMV infection, which revealed distinct promoter and enhancer activities between the two helper cell lineages. Furthermore, we identified that Bhlhe40 directly binds to cis-regulatory elements of Th1-related genes such as Tbx21 and Cxcr6 to activate their expression while simultaneously binding to regions of Tfh-related genes such as Bcl6 and Cxcr5 to repress their expression. Collectively, our data suggest that Bhlhe40 functions as a transcription activator to promote Th1 cell differentiation and a transcription repressor to suppress Tfh cell differentiation.

Methylation of the chromatin modifier KMT2D by SMYD2 contributes to therapeutic response in hormone-dependent breast cancer

Activating mutations in PIK3CA are frequently found in estrogen-receptor-positive (ER+) breast cancer, and the combination of the phosphatidylinositol 3-kinase (PI3K) inhibitor alpelisib with anti-ER inhibitors is approved for therapy. We have previously demonstrated that the PI3K pathway regulates ER activity through phosphorylation of the chromatin modifier KMT2D. Here, we discovered a methylation site on KMT2D, at K1330 directly adjacent to S1331, catalyzed by the lysine methyltransferase SMYD2. SMYD2 loss attenuates alpelisib-induced KMT2D chromatin binding and alpelisib-mediated changes in gene expression, including ER-dependent transcription. Knockdown or pharmacological inhibition of SMYD2 sensitizes breast cancer cells, patient-derived organoids, and tumors to PI3K/AKT inhibition and endocrine therapy in part through KMT2D K1330 methylation. Together, our findings uncover a regulatory crosstalk between post-translational modifications that fine-tunes KMT2D function at the chromatin. This provides a rationale for the use of SMYD2 inhibitors in combination with PI3Kα/AKT inhibitors in the treatment of ER+/PIK3CA mutant breast cancer.

Single-cell multi-cohort dissection of the schizophrenia transcriptome

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks. We additionally observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Finally, we demonstrated that loss of ATRX in a mesenchymal malignancy, undifferentiated pleomorphic sarcoma, results in similar epigenetic disruption and de-repression of transposable elements. Together, our results reveal a role for ATRX in maintaining epigenetic states and transcriptional repression in mesenchymal progenitors and tumor cells and in preventing aberrant differentiation in the progenitor context.

RNA Polymerase II hypertranscription in cancer FFPE samples

Hypertranscription is common in human cancers and predicts poor prognosis. However detection of hypertranscription is indirect, relying on accurately quantifying mRNA levels and estimating cell numbers. Previously, we introduced FFPE-CUTAC, a genome-wide method for mapping RNA Polymerase II (RNAPII) in formalin-fixed paraffin-embedded (FFPE) sections. Here we use FFPE-CUTAC to demonstrate genome-wide hypertranscription both in transgene-driven mouse gliomas and in assorted human tumors at active regulatory elements and replication-coupled histone genes with reduced mitochondrial DNA abundance. FFPE-CUTAC identified RNAPII-bound regulatory elements shared among diverse cancers and readily categorized human tumors despite using very small samples and low sequencing depths. Remarkably, RNAPII FFPE-CUTAC identified de novo and precisely mapped HER2 amplifications punctuated by likely selective sweeps including genes encoding direct positive regulators of RNAPII itself. Our results demonstrate that FFPE-CUTAC measurements of hypertranscription and classifications of tumors using small sections provides an affordable and sensitive genome-wide strategy for personalized medicine.

Enhancer switching in cell lineage priming is linked to eRNA, Brg1's AT-hook, and SWI/SNF recruitment

RNA transcribed from enhancers, i.e., eRNA, has been suggested to directly activate transcription by recruiting transcription factors and co-activators. Although there have been specific examples of eRNA functioning in this way, it is not clear how general this may be. We find that the AT-hook of SWI/SNF preferentially binds RNA and, as part of the esBAF complex, associates with eRNA transcribed from intronic and intergenic regions. Our data suggest that SWI/SNF is globally recruited in cis by eRNA to cell-type-specific enhancers, representative of two distinct stages that mimic early mammalian development, and not at enhancers that are shared between the two stages. In this manner, SWI/SNF facilitates recruitment and/or activation of MLL3/4, p300/CBP, and Mediator to stage-specific enhancers and super-enhancers that regulate the transcription of metabolic and cell lineage priming-related genes. These findings highlight a connection between ATP-dependent chromatin remodeling and eRNA in cell identity and typical- and super-enhancer activation.

Targeting BRD4 with PROTAC degrader ameliorates LPS-induced acute lung injury by inhibiting M1 alveolar macrophage polarization

OBJECTIVES

Acute lung injury (ALI) is a highly inflammatory condition with the involvement of M1 alveolar macrophages (AMs) polarization, eventually leading to the development of non-cardiogenic edema in alveolar and interstitial regions, accompanied by persistent hypoxemia. Given the significant mortality rate associated with ALI, it is imperative to investigate the underlying mechanisms of this condition so as to identify potential therapeutic targets. The therapeutic effects of the inhibition of bromodomain containing protein 4 (BRD4), an epigenetic reader, has been proven with high efficacy in ameliorating various inflammatory diseases through mediating immune cell activation. However, little is known about the therapeutic potential of BRD4 degradation in acute lung injury.

METHODS

This study aimed to assess the protective efficacy of ARV-825, a novel BRD4-targeted proteolysis targeting chimera (PROTAC), against ALI through histopathological examination in lung tissues and biochemical analysis in bronchoalveolar lavage fluid (BALF). Additionally, the underlying mechanism by which BRD4 regulated M1 AMs was elucidated by using CUT & Tag assay.

RESULTS

In this study, we found the upregulation of BRD4 in a lipopolysaccharide (LPS)-induced ALI model. Furthermore, we observed that intraperitoneal administration of ARV-825, significantly alleviated LPS-induced pulmonary pathological changes and inflammatory responses. These effects were accompanied by the suppression of M1 AMs. In addition, our findings revealed that the administration of ARV-825 effectively suppressed M1 AMs by inhibiting the expression of IRF7, a crucial transcriptional factor involved in M1 macrophages.

CONCLUSION

Our study suggested that targeting BRD4 using ARV-825 is a potential therapeutic approach for ALI.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

Targeting BIRC5 as a therapeutic approach to overcome ASXL1-associated decitabine resistance

Hypomethylating agents (HMAs) are widely employed in the treatment of myeloid malignancies. However, unresponsive or resistant to HMAs occurs in approximately 50 % of patients. ASXL1, one of the most commonly mutated genes across the full spectrum of myeloid malignancies, has been reported to predict a lower overall response rate to HMAs, suggesting an essential need to develop effective therapeutic strategies for the patients with HMA failure. Here, we investigated the impact of ASXL1 on cellular responsiveness to decitabine treatment. ASXL1 deficiency increased resistance to decitabine treatment in AML cell lines and mouse bone marrow cells. Transcriptome sequencing revealed significant alterations in genes regulating cell cycle, apoptosis, and histone modification in ASXL1 deficient cells that resistant to decitabine. BIRC5 was identified as a potential target for overcoming decitabine resistance in ASXL1 deficient cells. Furthermore, our experimental evidence demonstrated that the small-molecule inhibitor of BIRC5 (YM-155) synergistically sensitized ASXL1 deficient cells to decitabine treatment. This study sheds light on the molecular mechanisms underlying the ASXL1-associated HMA resistance and proposes a promising therapeutic strategy for improving treatment outcomes in affected individuals.

Histone demethylase KDM5 regulates cardiomyocyte maturation by promoting fatty acid oxidation, oxidative phosphorylation, and myofibrillar organization

AIMS

Human pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) provide a platform to identify and characterize factors that regulate the maturation of CMs. The transition from an immature foetal to an adult CM state entails coordinated regulation of the expression of genes involved in myofibril formation and oxidative phosphorylation (OXPHOS) among others. Lysine demethylase 5 (KDM5) specifically demethylates H3K4me1/2/3 and has emerged as potential regulators of expression of genes involved in cardiac development and mitochondrial function. The purpose of this study is to determine the role of KDM5 in iPSC-CM maturation.

METHODS AND RESULTS

KDM5A, B, and C proteins were mainly expressed in the early post-natal stages, and their expressions were progressively downregulated in the post-natal CMs and were absent in adult hearts and CMs. In contrast, KDM5 proteins were persistently expressed in the iPSC-CMs up to 60 days after the induction of myogenic differentiation, consistent with the immaturity of these cells. Inhibition of KDM5 by KDM5-C70 -a pan-KDM5 inhibitor, induced differential expression of 2372 genes, including upregulation of genes involved in fatty acid oxidation (FAO), OXPHOS, and myogenesis in the iPSC-CMs. Likewise, genome-wide profiling of H3K4me3 binding sites by the cleavage under targets and release using nuclease assay showed enriched of the H3K4me3 peaks at the promoter regions of genes encoding FAO, OXPHOS, and sarcomere proteins. Consistent with the chromatin and gene expression data, KDM5 inhibition increased the expression of multiple sarcomere proteins and enhanced myofibrillar organization. Furthermore, inhibition of KDM5 increased H3K4me3 deposits at the promoter region of the ESRRA gene and increased its RNA and protein levels. Knockdown of ESRRA in KDM5-C70-treated iPSC-CM suppressed expression of a subset of the KDM5 targets. In conjunction with changes in gene expression, KDM5 inhibition increased oxygen consumption rate and contractility in iPSC-CMs.

CONCLUSION

KDM5 inhibition enhances maturation of iPSC-CMs by epigenetically upregulating the expressions of OXPHOS, FAO, and sarcomere genes and enhancing myofibril organization and mitochondrial function.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation

Enteroendocrine cells (EECs) secrete serotonin (enterochromaffin [EC] cells) or specific peptide hormones (non-EC cells) that serve vital metabolic functions. The basis for terminal EEC diversity remains obscure. By forcing activity of the transcription factor (TF) NEUROG3 in 2D cultures of human intestinal stem cells, we replicated physiologic EEC differentiation and examined transcriptional and cis-regulatory dynamics that culminate in discrete cell types. Abundant EEC precursors expressed stage-specific genes and TFs. Before expressing pre-terminal NEUROD1, post-mitotic precursors oscillated between transcriptionally distinct ASCL1+ and HES6hi cell states. Loss of either factor accelerated EEC differentiation substantially and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and non-EC cell features. These TFs mainly bind cis-elements that are accessible in undifferentiated stem cells, and they tailor subsequent expression of TF combinations that underlie discrete EEC identities. Thus, early TF oscillations retard EEC maturation to enable accurate diversity within a medically important cell lineage.

Cysteine Rich Intestinal Protein 2 is a copper-responsive regulator of skeletal muscle differentiation

Copper (Cu) is an essential trace element required for respiration, neurotransmitter synthesis, oxidative stress response, and transcriptional regulation. Imbalance in Cu homeostasis can lead to several pathological conditions, affecting neuronal, cognitive, and muscular development. Mechanistically, Cu and Cu-binding proteins (Cu-BPs) have an important but underappreciated role in transcription regulation in mammalian cells. In this context, our lab investigates the contributions of novel Cu-BPs in skeletal muscle differentiation using murine primary myoblasts. Through an unbiased synchrotron X-ray fluorescence-mass spectrometry (XRF/MS) metalloproteomic approach, we identified the murine cysteine rich intestinal protein 2 (mCrip2) in a sample that showed enriched Cu signal, which was isolated from differentiating primary myoblasts derived from mouse satellite cells. Immunolocalization analyses showed that mCrip2 is abundant in both nuclear and cytosolic fractions. Thus, we hypothesized that mCrip2 might have differential roles depending on its cellular localization in the skeletal muscle lineage. mCrip2 is a LIM-family protein with 4 conserved Zn2+-binding sites. Homology and phylogenetic analyses showed that mammalian Crip2 possesses histidine residues near two of the Zn2+-binding sites (CX2C-HX2C) which are potentially implicated in Cu+-binding and competition with Zn2+. Biochemical characterization of recombinant human hsCRIP2 revealed a high Cu+-binding affinity for two and four Cu+ ions and limited redox potential. Functional characterization using CRISPR/Cas9-mediated deletion of mCrip2 in primary myoblasts did not impact proliferation, but impaired myogenesis by decreasing the expression of differentiation markers, possibly attributed to Cu accumulation. Transcriptome analyses of proliferating and differentiating mCrip2 KO myoblasts showed alterations in mRNA processing, protein translation, ribosome synthesis, and chromatin organization. CUT&RUN analyses showed that mCrip2 associates with a select set of gene promoters, including MyoD1 and metallothioneins, acting as a novel Cu-responsive or Cu-regulating protein. Our work demonstrates novel regulatory functions of mCrip2 that mediate skeletal muscle differentiation, presenting new features of the Cu-network in myoblasts.

Human ASXL1-Mutant Hematopoiesis Is Driven by a Truncated Protein Associated with Aberrant Deubiquitination of H2AK119

Mutations in additional sex combs like 1 (ASXL1) confer poor prognosis both in myeloid malignancies and in premalignant clonal hematopoiesis (CH). However, the mechanisms by which these mutations contribute to disease initiation remain unresolved, and mutation-specific targeting has remained elusive. To address this, we developed a human disease model that recapitulates the disease trajectory from ASXL1-mutant CH to lethal myeloid malignancy. We demonstrate that mutations in ASXL1 lead to the expression of a functional, truncated protein and determine that truncated ASXL1 leads to global redistribution of the repressive chromatin mark H2AK119Ub, increased transposase-accessible chromatin, and activation of both myeloid and stem cell gene-expression programs. Finally, we demonstrate that H2AK119Ub levels are tied to truncated ASXL1 expression levels and leverage this observation to demonstrate that inhibition of the PRC1 complex might be an ASXL1-mutant-specific therapeutic vulnerability in both premalignant CH and myeloid malignancy.

SIGNIFICANCE

Mutant ASXL1 is a common driver of CH and myeloid malignancy. Using primary human HSPCs, we determine that truncated ASXL1 leads to redistribution of H2AK119Ub and may affect therapeutic vulnerability to PRC1 inhibition.

Histone lactylation in macrophages is predictive for gene expression changes during ischemia induced-muscle regeneration

OBJECTIVES

We have previously shown that lactate is an essential metabolite for macrophage polarisation during ischemia-induced muscle regeneration. Recent in vitro work has implicated histone lactylation, a direct derivative of lactate, in macrophage polarisation. Here, we explore the in vivo relevance of histone lactylation for macrophage polarisation after muscle injury.

METHODS

To evaluate macrophage dynamics during muscle regeneration, we subjected mice to ischemia-induced muscle damage by ligating the femoral artery. Muscle samples were harvested at 1, 2, 4, and 7 days post injury (dpi). CD45+CD11b+F4/80+CD64+ macrophages were isolated and processed for RNA sequencing, Western Blotting, and CUT&Tag-sequencing to investigate gene expression, histone lactylation levels, and histone lactylation genomic localisation and enrichment, respectively.

RESULTS

We show that, over time, macrophages in the injured muscle undergo extensive gene expression changes, which are similar in nature and in timing to those seen after other types of muscle-injuries. We find that the macrophage histone lactylome is modified between 2 and 4 dpi, which is a crucial window for macrophage polarisation. Absolute histone lactylation levels increase, and, although subtly, the genomic enrichment of H3K18la changes. Overall, we find that histone lactylation is important at both promoter and enhancer elements. Lastly, H3K18la genomic profile changes from 2 to 4 dpi were predictive for gene expression changes later in time, rather than being a reflection of prior gene expression changes.

CONCLUSIONS

Our results suggest that histone lactylation dynamics are functionally important for the function of macrophages during muscle regeneration.

An imbalance between proliferation and differentiation underlies the development of microRNA-defective pineoblastoma

Mutations in the microRNA processing genes DICER1 and DROSHA drive several cancers that resemble embryonic progenitors. To understand how microRNAs regulate tumorigenesis, we ablated Drosha or Dicer1 in the developing pineal gland to emulate the pathogenesis of pineoblastoma, a brain tumor that resembles undifferentiated precursors of the pineal gland. Accordingly, these mice develop pineal tumors marked by loss of microRNAs, including the let-7/miR-98-5p family, and de-repression of microRNA target genes. Pineal tumors driven by loss of Drosha or Dicer1 mimic tumors driven by Rb1 loss, as they exhibit upregulation of S-phase genes and homeobox transcription factors that regulate pineal development. Blocking proliferation of these tumors facilitates expression of pinealocyte maturation markers, with a concomitant reduction in embryonic markers. Select embryonic markers remain elevated, however, as the microRNAs that normally repress these target genes remain absent. One such microRNA target gene is the oncofetal transcription factor Plagl2, which regulates expression of pro-growth genes, and inhibiting their signaling impairs tumor growth. Thus, we demonstrate that tumors driven by loss of microRNA processing may be therapeutically targeted by inhibiting downstream drivers of proliferation.

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.

Class I HDAC inhibitors enhance antitumor efficacy and persistence of CAR-T cells by activation of the Wnt pathway

Epigenetic modification shapes differentiation trajectory and regulates the exhaustion state of chimeric antigen receptor T (CAR-T) cells. Limited efficacy induced by terminal exhaustion closely ties with intrinsic transcriptional regulation. However, the comprehensive regulatory mechanisms remain largely elusive. Here, we identify class I histone deacetylase inhibitors (HDACi) as boosters of CAR-T cell function by high-throughput screening of chromatin-modifying drugs, in which M344 and chidamide enhance memory maintenance and resistance to exhaustion of CAR-T cells that induce sustained antitumor efficacy both in vitro and in vivo. Mechanistically, HDACi decrease HDAC1 expression and enhance H3K27ac activity. Multi-omics analyses from RNA-seq, ATAC-seq, and H3K27ac CUT&Tag-seq show that HDACi upregulate expression of TCF4, LEF1, and CTNNB1, which subsequently activate the canonical Wnt/β-catenin pathway. Collectively, our findings elucidate the functional roles of class I HDACi in enhancing CAR-T cell function, which provides the basis and therapeutic targets for synergic combination of CAR-T cell therapy and HDACi treatment.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Transcription elongation defects link oncogenic SF3B1 mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human diseases remains unexplored. Using isogenic cell lines, patient samples, and a mutant mouse model, we investigated how cancer-associated mutations in SF3B1 alter transcription. We found that these mutations reduce the elongation rate of RNA polymerase II (RNAPII) along gene bodies and its density at promoters. The elongation defect results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1. The decreased promoter-proximal RNAPII density reduces both chromatin accessibility and H3K4me3 marks at promoters. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC/H3K4me pathway, which, when modulated, reverse both transcription and chromatin changes. Our findings reveal how splicing factor mutant states behave functionally as epigenetic disorders through impaired transcription-related changes to the chromatin landscape. We also present a rationale for targeting the Sin3/HDAC complex as a therapeutic strategy.

Widespread impact of nucleosome remodelers on transcription at cis-regulatory elements

Nucleosome remodeling complexes and other regulatory factors work in concert to build a chromatin environment that directs the expression of a distinct set of genes in each cell using cis-regulatory elements (CREs), such as promoters and enhancers, that drive transcription of both mRNAs and CRE-associated non-coding RNAs (ncRNAs). Two classes of CRE-associated ncRNAs include upstream antisense RNAs (uaRNAs), which are transcribed divergently from a shared mRNA promoter, and enhancer RNAs (eRNAs), which are transcribed bidirectionally from active enhancers. The complicated network of CRE regulation by nucleosome remodelers remains only partially explored, with a focus on a select, limited number of remodelers. We endeavored to elucidate a remodeler-based regulatory network governing CRE-associated transcription (mRNA, eRNA, and uaRNA) in murine embryonic stem (ES) cells to test the hypothesis that many SNF2-family nucleosome remodelers collaborate to regulate the coding and non-coding transcriptome via alteration of underlying nucleosome architecture. Using depletion followed by transient transcriptome sequencing (TT-seq), we identified thousands of misregulated mRNAs and CRE-associated ncRNAs across the remodelers examined, identifying novel contributions by understudied remodelers in the regulation of coding and noncoding transcription. Our findings suggest that mRNA and eRNA transcription are coordinately co-regulated, while mRNA and uaRNAs sharing a common promoter are independently regulated. Subsequent mechanistic studies suggest that while remodelers SRCAP and CHD8 modulate transcription through classical mechanisms such as transcription factors and histone variants, a broad set of remodelers including SMARCAL1 indirectly contribute to transcriptional regulation through maintenance of genomic stability and proper Integrator complex localization. This study systematically examines the contribution of SNF2-remodelers to the CRE-associated transcriptome, identifying at least two classes for remodeler action.

ARID1A orchestrates SWI/SNF-mediated sequential binding of transcription factors with ARID1A loss driving pre-memory B cell fate and lymphomagenesis

ARID1A, a subunit of the canonical BAF nucleosome remodeling complex, is commonly mutated in lymphomas. We show that ARID1A orchestrates B cell fate during the germinal center (GC) response, facilitating cooperative and sequential binding of PU.1 and NF-kB at crucial genes for cytokine and CD40 signaling. The absence of ARID1A tilts GC cell fate toward immature IgM+CD80-PD-L2- memory B cells, known for their potential to re-enter new GCs. When combined with BCL2 oncogene, ARID1A haploinsufficiency hastens the progression of aggressive follicular lymphomas (FLs) in mice. Patients with FL with ARID1A-inactivating mutations preferentially display an immature memory B cell-like state with increased transformation risk to aggressive disease. These observations offer mechanistic understanding into the emergence of both indolent and aggressive ARID1A-mutant lymphomas through the formation of immature memory-like clonal precursors. Lastly, we demonstrate that ARID1A mutation induces synthetic lethality to SMARCA2/4 inhibition, paving the way for potential precision therapy for high-risk patients.

The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Affinity-tagged SMAD1 and SMAD5 mouse lines reveal transcriptional reprogramming mechanisms during early pregnancy

Endometrial decidualization, a prerequisite for successful pregnancies, relies on transcriptional reprogramming driven by progesterone receptor (PR) and bone morphogenetic protein (BMP)-SMAD1/SMAD5 signaling pathways. Despite their critical roles in early pregnancy, how these pathways intersect in reprogramming the endometrium into a receptive state remains unclear. To define how SMAD1 and/or SMAD5 integrate BMP signaling in the uterus during early pregnancy, we generated two novel transgenic mouse lines with affinity tags inserted into the endogenous SMAD1 and SMAD5 loci (Smad1HA/HA and Smad5PA/PA). By profiling the genome-wide distribution of SMAD1, SMAD5, and PR in the mouse uterus, we demonstrated the unique and shared roles of SMAD1 and SMAD5 during the window of implantation. We also showed the presence of a conserved SMAD1, SMAD5, and PR genomic binding signature in the uterus during early pregnancy. To functionally characterize the translational aspects of our findings, we demonstrated that SMAD1/5 knockdown in human endometrial stromal cells suppressed expressions of canonical decidual markers (IGFBP1, PRL, FOXO1) and PR-responsive genes (RORB, KLF15). Here, our studies provide novel tools to study BMP signaling pathways and highlight the fundamental roles of SMAD1/5 in mediating both BMP signaling pathways and the transcriptional response to progesterone (P4) during early pregnancy.

Altered DNA methylation underlies monocyte dysregulation and immune exhaustion memory in sepsis

Monocytes can develop an exhausted memory state characterized by reduced differentiation, pathogenic inflammation, and immune suppression that drives immune dysregulation during sepsis. Chromatin alterations, notably via histone modifications, underlie innate immune memory, but the contribution of DNA methylation remains poorly understood. Using an ex vivo sepsis model, we show altered DNA methylation throughout the genome of exhausted monocytes, including genes implicated in immune dysregulation during sepsis and COVID-19 infection (e.g., Plac8). These changes are recapitulated in septic mice induced by cecal slurry injection. Methylation profiles developed in septic mice are maintained during ex vivo culture, supporting the involvement of DNA methylation in stable monocyte exhaustion memory. Methylome reprogramming is driven in part by Wnt signaling inhibition in exhausted monocytes and can be reversed with DNA methyltransferase inhibitors, Wnt agonists, or immune training molecules. Our study demonstrates the significance of altered DNA methylation in the maintenance of stable monocyte exhaustion memory.

Targeting the chromatin binding of exportin-1 disrupts NFAT and T cell activation

Exportin-1 (XPO1/CRM1) plays a central role in the nuclear-to-cytoplasmic transport of hundreds of proteins and contributes to other cellular processes, such as centrosome duplication. Small molecules targeting XPO1 induce cytotoxicity, and selinexor was approved by the Food and Drug Administration in 2019 as a cancer chemotherapy for relapsed multiple myeloma. Here, we describe a cell-type-dependent chromatin-binding function for XPO1 that is essential for the chromatin occupancy of NFAT transcription factors and thus the appropriate activation of T cells. Additionally, we establish a class of XPO1-targeting small molecules capable of disrupting the chromatin binding of XPO1 without perturbing nuclear export or inducing cytotoxicity. This work defines a broad transcription regulatory role for XPO1 that is essential for T cell activation as well as a new class of XPO1 modulators to enable therapeutic targeting of XPO1 beyond oncology including in T cell-driven autoimmune disorders.

Prediction of DNA i-motifs via machine learning

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine.

PRC1 directs PRC2-H3K27me3 deposition to shield adult spermatogonial stem cells from differentiation

Spermatogonial stem cells functionality reside in the slow-cycling and heterogeneous undifferentiated spermatogonia cell population. This pool of cells supports lifelong fertility in adult males by balancing self-renewal and differentiation to produce haploid gametes. However, the molecular mechanisms underpinning long-term stemness of undifferentiated spermatogonia during adulthood remain unclear. Here, we discover that an epigenetic regulator, Polycomb repressive complex 1 (PRC1), shields adult undifferentiated spermatogonia from differentiation, maintains slow cycling, and directs commitment to differentiation during steady-state spermatogenesis in adults. We show that PRC2-mediated H3K27me3 is an epigenetic hallmark of adult undifferentiated spermatogonia. Indeed, spermatogonial differentiation is accompanied by a global loss of H3K27me3. Disruption of PRC1 impairs global H3K27me3 deposition, leading to precocious spermatogonial differentiation. Therefore, PRC1 directs PRC2-H3K27me3 deposition to maintain the self-renewing state of undifferentiated spermatogonia. Importantly, in contrast to its role in other tissue stem cells, PRC1 negatively regulates the cell cycle to maintain slow cycling of undifferentiated spermatogonia. Our findings have implications for how epigenetic regulators can be tuned to regulate the stem cell potential, cell cycle and differentiation to ensure lifelong fertility in adult males.

Cancer-associated DNA Hypermethylation of Polycomb Targets Requires DNMT3A Dual Recognition of Histone H2AK119 Ubiquitination and the Nucleosome Acidic Patch

During tumor development, promoter CpG islands (CGIs) that are normally silenced by Polycomb repressive complexes (PRCs) become DNA hypermethylated. The molecular mechanism by which de novo DNA methyltransferase(s) catalyze CpG methylation at PRC-regulated regions remains unclear. Here we report a cryo-EM structure of the DNMT3A long isoform (DNMT3A1) N-terminal region in complex with a nucleosome carrying PRC1-mediated histone H2A lysine 119 monoubiquitination (H2AK119Ub). We identify regions within the DNMT3A1 N-terminus that bind H2AK119Ub and the nucleosome acidic patch. This bidentate interaction is required for effective DNMT3A1 engagement with H2AK119Ub-modified chromatin in cells. Furthermore, aberrant redistribution of DNMT3A1 to Polycomb target genes inhibits their transcriptional activation during cell differentiation and recapitulates the cancer-associated DNA hypermethylation signature. This effect is rescued by disruption of the DNMT3A1-acidic patch interaction. Together, our analyses reveal a binding interface critical for countering promoter CGI DNA hypermethylation, a major molecular hallmark of cancer.

A fluorescence-based protocol to quantitatively titrate CUT&RUN buffer components

Cleavage under targets & release using nuclease (CUT&RUN) is a technique for identifying genomic sites where proteins or histone modifications are present in chromatin in permeabilized cells. Here, we present a fluorescence-based protocol to quantitatively titrate CUT&RUN buffer components, for efficient cell permeabilization and retention of target epitopes on chromatin. We describe steps for capturing cells on concanavalin A beads and using a fluorescently labeled secondary antibody to titrate concentrations of digitonin and NaCl in CUT&RUN buffers. We then detail procedures for fluorescence imaging to identify optimal conditions. For complete details on the use and execution of this protocol, please refer to Lerner et al.1.

Krüppel-like Factor (KLF) family members control expression of genes required for serous cavity and alveolar macrophage identities

Tissue-resident macrophages adopt distinct gene expression profiles and exhibit functional specialization based on their tissue of residence. Recent studies have begun to define the signals and transcription factors that induce these identities. Here we describe an unexpected and specific role for the broadly expressed transcription factor Kruppel-like Factor 2 (KLF2) in the development of embryonically derived Large Cavity Macrophages (LCM) in the serous cavities. KLF2 not only directly regulates the transcription of genes previously shown to specify LCM identity, such as retinoic acid receptors and GATA6, but also is required for induction of many other transcripts that define the identity of these cells. We identify a similar role for KLF4 in regulating the identity of alveolar macrophages in the lung. These data demonstrate that broadly expressed transcription factors, such as Group 2 KLFs, can play important roles in the specification of distinct identities of tissue-resident macrophages.

BET protein-dependent E2F pathway activity confers bell-shaped type resistance to tankyrase inhibitors in APC-mutated colorectal cancer

WNT/β-catenin signaling is aberrantly activated in colorectal cancer (CRC) mainly by loss-of-function mutations in adenomatous polyposis coli (APC) and is involved in tumor progression. Tankyrase inhibitors, which suppress WNT/β-catenin signaling, are currently in pre-clinical and clinical trials. However, the mechanisms of resistance to tankyrase inhibitors remain unclear. In this study, we established tankyrase inhibitor-resistant CRC cells, JC73-RK100, from APC-mutated patient-derived CRC cells. JC73-RK100 cells and several CRC cell lines were sensitive to tankyrase inhibitors at low concentrations but were resistant at high concentrations, showing an intrinsic/acquired bell-shaped dose response. Mechanistically, tankyrase inhibitors at high concentrations promoted BRD3/4-dependent E2F target gene transcription and over-activated cell cycle progression in these cells. BET inhibitors canceled the bell-shaped dose response to tankyrase inhibitors. Combination of tankyrase and BET inhibitors significantly suppressed tumor growth in a mouse xenograft model. These observations suggest that the combination of tankyrase and BET inhibitors may be a useful therapeutic approach to overcome the resistance of a subset of CRCs to tankyrase inhibitors.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

SOX17 enables immune evasion of early colorectal adenomas and cancers

A hallmark of cancer is the avoidance of immune destruction. This process has been primarily investigated in locally advanced or metastatic cancer1-3; however, much less is known about how pre-malignant or early invasive tumours evade immune detection. Here, to understand this process in early colorectal cancers (CRCs), we investigated how naive colon cancer organoids that were engineered in vitro to harbour Apc-null, KrasG12D and Trp53-null (AKP) mutations adapted to the in vivo native colonic environment. Comprehensive transcriptomic and chromatin analyses revealed that the endoderm-specifying transcription factor SOX17 became strongly upregulated in vivo. Notably, whereas SOX17 loss did not affect AKP organoid propagation in vitro, its loss markedly reduced the ability of AKP tumours to persist in vivo. The small fraction of SOX17-null tumours that grew displayed notable interferon-γ (IFNγ)-producing effector-like CD8+ T cell infiltrates in contrast to the immune-suppressive microenvironment in wild-type counterparts. Mechanistically, in both endogenous Apc-null pre-malignant adenomas and transplanted organoid-derived AKP CRCs, SOX17 suppresses the ability of tumour cells to sense and respond to IFNγ, preventing anti-tumour T cell responses. Finally, SOX17 engages a fetal intestinal programme that drives differentiation away from LGR5+ tumour cells to produce immune-evasive LGR5- tumour cells with lower expression of major histocompatibility complex class I (MHC-I). We propose that SOX17 is a transcription factor that is engaged during the early steps of colon cancer to orchestrate an immune-evasive programme that permits CRC initiation and progression.

N-MYC regulates cell survival via eIF4G1 in inv(16) acute myeloid leukemia

N-MYC (encoded by MYCN) is a critical regulator of hematopoietic stem cell function. While the role of N-MYC deregulation is well established in neuroblastoma, the importance of N-MYC deregulation in leukemogenesis remains elusive. Here, we demonstrate that N-MYC is overexpressed in acute myeloid leukemia (AML) cells with chromosome inversion inv(16) and contributes to the survival and maintenance of inv(16) leukemia. We identified a previously unknown MYCN enhancer, active in multiple AML subtypes, essential for MYCN mRNA levels and survival in inv(16) AML cells. We also identified eukaryotic translation initiation factor 4 gamma 1 (eIF4G1) as a key N-MYC target that sustains leukemic survival in inv(16) AML cells. The oncogenic role of eIF4G1 in AML has not been reported before. Our results reveal a mechanism whereby N-MYC drives a leukemic transcriptional program and provides a rationale for the therapeutic targeting of the N-MYC/eIF4G1 axis in myeloid leukemia.

Single-cell chromatin profiling reveals genetic programs activating proregenerative states in nonmyocyte cells

Cardiac regeneration requires coordinated participation of multiple cell types whereby their communications result in transient activation of proregenerative cell states. Although the molecular characteristics and lineage origins of these activated cell states and their contribution to cardiac regeneration have been studied, the extracellular signaling and the intrinsic genetic program underlying the activation of the transient functional cell states remain largely unexplored. In this study, we delineated the chromatin landscapes of the noncardiomyocytes (nonCMs) of the regenerating heart at the single-cell level and inferred the cis-regulatory architectures and trans-acting factors that control cell type-specific gene expression programs. Moreover, further motif analysis and cell-specific genetic manipulations suggest that the macrophage-derived inflammatory signal tumor necrosis factor-α, acting via its downstream transcription factor complex activator protein-1, functions cooperatively with discrete transcription regulators to activate respective nonCM cell types critical for cardiac regeneration. Thus, our study defines the regulatory architectures and intercellular communication principles in zebrafish heart regeneration.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Transcription factor co-expression mediates lineage priming for embryonic and extra-embryonic differentiation

In early mammalian development, cleavage stage blastomeres and inner cell mass (ICM) cells co-express embryonic and extra-embryonic transcriptional determinants. Using a protein-based double reporter we identify an embryonic stem cell (ESC) population that co-expresses the extra-embryonic factor GATA6 alongside the embryonic factor SOX2. Based on single cell transcriptomics, we find this population resembles the unsegregated ICM, exhibiting enhanced differentiation potential for endoderm while maintaining epiblast competence. To relate transcription factor binding in these cells to future fate, we describe a complete enhancer set in both ESCs and naive extra-embryonic endoderm stem cells and assess SOX2 and GATA6 binding at these elements in the ICM-like ESC sub-population. Both factors support cooperative recognition in these lineages, with GATA6 bound alongside SOX2 on a fraction of pluripotency enhancers and SOX2 alongside GATA6 more extensively on endoderm enhancers, suggesting that cooperative binding between these antagonistic factors both supports self-renewal and prepares progenitor cells for later differentiation.

KDM6A/UTX promotes spermatogenic gene expression across generations and is not required for male fertility†

Paternal chromatin undergoes extensive structural and epigenetic changes during mammalian spermatogenesis, producing sperm with an epigenome optimized for the transition to embryogenesis. Lysine demethylase 6a (KDM6A, also called UTX) promotes gene activation in part via demethylation of H3K27me3, a developmentally important repressive modification abundant throughout the epigenome of spermatogenic cells and sperm. We previously demonstrated increased cancer risk in genetically wild-type mice derived from a paternal germ line lacking Kdm6a (Kdm6a cKO), indicating a role for KDM6A in regulating heritable epigenetic states. However, the regulatory function of KDM6A during spermatogenesis is not known. Here, we show that Kdm6a is transiently expressed in spermatogenesis, with RNA and protein expression largely limited to late spermatogonia and early meiotic prophase. Kdm6a cKO males do not have defects in fertility or the overall progression of spermatogenesis. However, hundreds of genes are deregulated upon loss of Kdm6a in spermatogenic cells, with a strong bias toward downregulation coinciding with the time when Kdm6a is expressed. Misregulated genes encode factors involved in chromatin organization and regulation of repetitive elements, and a subset of these genes was persistently deregulated in the male germ line across two generations of offspring of Kdm6a cKO males. Genome-wide epigenetic profiling revealed broadening of H3K27me3 peaks in differentiating spermatogonia of Kdm6a cKO mice, suggesting that KDM6A demarcates H3K27me3 domains in the male germ line. Our findings highlight KDM6A as a transcriptional activator in the mammalian male germ line that is dispensable for spermatogenesis but important for safeguarding gene regulatory state intergenerationally.

Affinity-tagged SMAD1 and SMAD5 mouse lines reveal transcriptional reprogramming mechanisms during early pregnancy

Endometrial decidualization, a prerequisite for successful pregnancies, relies on transcriptional reprogramming driven by progesterone receptor (PR) and bone morphogenetic protein (BMP)-SMAD1/SMAD5 signaling pathways. Despite their critical roles in early pregnancy, how these pathways intersect in reprogramming the endometrium into a receptive state remains unclear. To define how SMAD1 and/or SMAD5 integrate BMP signaling in the uterus during early pregnancy, we generated two novel transgenic mouse lines with affinity tags inserted into the endogenous SMAD1 and SMAD5 loci (Smad1HA/HA and Smad5PA/PA). By profiling the genome-wide distribution of SMAD1, SMAD5, and PR in the mouse uterus, we demonstrated the unique and shared roles of SMAD1 and SMAD5 during the window of implantation. We also showed the presence of a conserved SMAD1, SMAD5, and PR genomic binding signature in the uterus during early pregnancy. To functionally characterize the translational aspects of our findings, we demonstrated that SMAD1/5 knockdown in human endometrial stromal cells suppressed expressions of canonical decidual markers (IGFBP1, PRL, FOXO1) and PR-responsive genes (RORB, KLF15). Here, our studies provide novel tools to study BMP signaling pathways and highlight the fundamental roles of SMAD1/5 in mediating both BMP signaling pathways and the transcriptional response to progesterone (P4) during early pregnancy.

Nucleosome conformation dictates the histone code

Histone post-translational modifications (PTMs) play a critical role in chromatin regulation. It has been proposed that these PTMs form localized 'codes' that are read by specialized regions (reader domains) in chromatin-associated proteins (CAPs) to regulate downstream function. Substantial effort has been made to define [CAP: histone PTM] specificities, and thus decipher the histone code and guide epigenetic therapies. However, this has largely been done using the reductive approach of isolated reader domains and histone peptides, which cannot account for any higher-order factors. Here, we show that the [BPTF PHD finger and bromodomain: histone PTM] interaction is dependent on nucleosome context. The tandem reader selectively associates with nucleosomal H3K4me3 and H3K14ac or H3K18ac, a combinatorial engagement that despite being in cis is not predicted by peptides. This in vitro specificity of the BPTF tandem reader for PTM-defined nucleosomes is recapitulated in a cellular context. We propose that regulatable histone tail accessibility and its impact on the binding potential of reader domains necessitates we refine the 'histone code' concept and interrogate it at the nucleosome level.

HRDE-2 drives small RNA specificity for the nuclear Argonaute protein HRDE-1

RNA interference (RNAi) is a conserved gene silencing process that exists in diverse organisms to protect genome integrity and regulate gene expression. In C. elegans, the majority of RNAi pathway proteins localize to perinuclear, phase-separated germ granules, which are comprised of sub-domains referred to as P granules, Mutator foci, Z granules, and SIMR foci. However, the protein components and function of the newly discovered SIMR foci are unknown. Here we demonstrate that HRDE-2 localizes to SIMR foci and interacts with the germline nuclear Argonaute HRDE-1 in its small RNA unbound state. In the absence of HRDE-2, HRDE-1 exclusively loads CSR-class 22G-RNAs rather than WAGO-class 22G-RNAs, resulting in inappropriate H3K9me3 deposition on CSR-target genes. Thus, our study demonstrates that the recruitment of unloaded HRDE-1 to germ granules, mediated by HRDE-2, is critical to ensure that the correct small RNAs are used to guide nuclear RNA silencing in the C. elegans germline.

Pycallingcards: an integrated environment for visualizing, analyzing, and interpreting Calling Cards data

MOTIVATION

Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data.

RESULTS

We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data.

AVAILABILITY AND IMPLEMENTATION

Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

Chromatin remodeling protein CHD4 regulates axon guidance of spiral ganglion neurons in developing cochlea

The chromodomain helicase binding protein 4 (CHD4) is an ATP-dependent chromatin remodeler. De-novo pathogenic variants of CHD4 cause Sifrim-Hitz-Weiss syndrome (SIHIWES). Patients with SIHIWES show delayed development, intellectual disability, facial dysmorphism, and hearing loss. Many cochlear cell types, including spiral ganglion neurons (SGNs), express CHD4. SGNs are the primary afferent neurons that convey sound information from the cochlea, but the function of CHD4 in SGNs is unknown. We employed the Neurog1(Ngn1) CreERT2 Chd4 conditional knockout animals to delete Chd4 in SGNs. SGNs are classified as type I and type II neurons. SGNs lacking CHD4 showed abnormal fasciculation of type I neurons along with improper pathfinding of type II fibers. CHD4 binding to chromatin from immortalized multipotent otic progenitor-derived neurons was used to identify candidate target genes in SGNs. Gene ontology analysis of CHD4 target genes revealed cellular processes involved in axon guidance, axonal fasciculation, and ephrin receptor signaling pathway. We validated increased Epha4 transcripts in SGNs from Chd4 conditional knockout cochleae. The results suggest that CHD4 attenuates the transcription of axon guidance genes to form the stereotypic pattern of SGN peripheral projections. The results implicate epigenetic changes in circuit wiring by modulating axon guidance molecule expression and provide insights into neurodevelopmental diseases.

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.

Disruption of H3K36 methylation provokes cellular plasticity to drive aberrant glandular formation and squamous carcinogenesis

Chromatin organization is essential for maintaining cell-fate trajectories and developmental programs. Here, we find that disruption of H3K36 methylation dramatically impairs normal epithelial differentiation and development, which promotes increased cellular plasticity and enrichment of alternative cell fates. Specifically, we observe a striking increase in the aberrant generation of excessive epithelial glandular tissues, including hypertrophic salivary, sebaceous, and meibomian glands, as well as enhanced squamous tumorigenesis. These phenotypic and gene expression manifestations are associated with loss of H3K36me2 and rewiring of repressive H3K27me3, changes we also observe in human patients with glandular hyperplasia. Collectively, these results have identified a critical role for H3K36 methylation in both in vivo epithelial cell-fate decisions and the prevention of squamous carcinogenesis and suggest that H3K36 methylation modulation may offer new avenues for the treatment of numerous common disorders driven by altered glandular function, which collectively affect large segments of the human population.

Epstein-Barr virus induces germinal center light zone chromatin architecture and promotes survival through enhancer looping at the BCL2A1 locus

IMPORTANCE

Epstein-Barr virus has evolved with its human host leading to an intimate relationship where infection of antibody-producing B cells mimics the process by which these cells normally recognize foreign antigens and become activated. Virtually everyone in the world is infected by adulthood and controls this virus pushing it into life-long latency. However, immune-suppressed individuals are at high risk for EBV+ cancers. Here, we isolated B cells from tonsils and compare the underlying molecular genetic differences between these cells and those infected with EBV. We find similar regulatory mechanism for expression of an important cellular protein that enables B cells to survive in lymphoid tissue. These findings link an underlying relationship at the molecular level between EBV-infected B cells in vitro with normally activated B cells in vivo. Our studies also characterize the role of a key viral control mechanism for B cell survival involved in long-term infection.

Role of histone variants H2BC1 and H2AZ.2 in H2AK119ub nucleosome organization and Polycomb gene silencing

Ubiquitination of histone H2A at lysine 119 residue (H2AK119ub) plays critical roles in a wide range of physiological processes, including Polycomb gene silencing 1,2 , replication 3-5 , DNA damage repair 6-10 , X inactivation 11,12 , and heterochromatin organization 13,14 . However, the underlying mechanism and structural basis of H2AK119ub remains largely elusive. In this study, we report that H2AK119ub nucleosomes have a unique composition, containing histone variants H2BC1 and H2AZ.2, and importantly, this composition is required for H2AK119ub and Polycomb gene silencing. Using the UAB domain of RSF1, we purified H2AK119ub nucleosomes to a sufficient amount and purity. Mass spectrometry analyses revealed that H2AK119ub nucleosomes contain the histone variants H2BC1 and H2AZ.2. A cryo-EM study resolved the structure of native H2AK119ub nucleosomes to a 2.6A resolution, confirming H2BC1 in one subgroup of H2AK119ub nucleosomes. Tandem GST-UAB pulldown, Flag-H2AZ.2, and HA-H2BC1 immunoprecipitation revealed that H2AK119ub nucleosomes could be separated into distinct subgroups, suggesting their composition heterogeneity and potential dynamic organization. Knockout or knockdown of H2BC1 or H2AZ.2 reduced cellular H2AK119ub levels, establishing H2BC1 and H2AZ.2 as critical determinants of H2AK119ub. Furthermore, genomic binding profiles of H2BC1 and H2AZ.2 overlapped significantly with H2AK119ub binding, with the most significant overlapping in the gene body and intergenic regions. Finally, assays in developing embryos reveal an interaction of H2AZ.2, H2BC1, and RING1A in vivo . Thus, this study revealed, for the first time, that the H2AK119ub nucleosome has a unique composition, and this composition is required for H2AK119ub and Polycomb gene silencing.