Improved CUT&RUN chromatin profiling tools

New tools to study renal fibrogenesis

PURPOSE OF REVIEW

Kidney fibrosis is a key pathological aspect and outcome of chronic kidney disease (CKD). The advent of multiomic analyses using human kidney tissue, enabled by technological advances, marks a new chapter of discovery in fibrosis research of the kidney. This review highlights the rapid advancements of single-cell and spatial multiomic techniques that offer new avenues for exploring research questions related to human kidney fibrosis development.

RECENT FINDINGS

We recently focused on understanding the origin and transition of myofibroblasts in kidney fibrosis using single-cell RNA sequencing (scRNA-seq) [1] . We analysed cells from healthy human kidneys and compared them to patient samples with CKD. We identified PDGFRα+/PDGFRβ+ mesenchymal cells as the primary cellular source of extracellular matrix (ECM) in human kidney fibrosis. We found several commonly shared cell states of fibroblasts and myofibroblasts and provided insights into molecular regulators. Novel single-cell and spatial multiomics tools are now available to shed light on cell lineages, the plasticity of kidney cells and cell-cell communication in fibrosis.

SUMMARY

As further single-cell and spatial multiomic approaches are being developed, opportunities to apply these methods to human kidney tissues expand similarly. Careful design and optimisation of the multiomic experiments are needed to answer questions related to cell lineages, plasticity and cell-cell communication in kidney fibrosis.

A conserved Pol II elongator SPT6L mediates Pol V transcription to regulate RNA-directed DNA methylation in Arabidopsis

In plants, the plant-specific RNA polymerase V (Pol V) transcripts non-coding RNAs and provides a docking platform for the association of accessory proteins in the RNA-directed DNA methylation (RdDM) pathway. Various components have been uncovered that are involved in the process of DNA methylation, but it is still not clear how the transcription of Pol V is regulated. Here, we report that the conserved RNA polymerase II (Pol II) elongator, SPT6L, binds to thousands of intergenic regions in a Pol II-independent manner. The intergenic enrichment of SPT6L, interestingly, co-occupies with the largest subunit of Pol V (NRPE1) and mutation of SPT6L leads to the reduction of DNA methylation but not Pol V enrichment. Furthermore, the association of SPT6L at Pol V loci is dependent on the Pol V associated factor, SPT5L, rather than the presence of Pol V, and the interaction between SPT6L and NRPE1 is compromised in spt5l. Finally, Pol V RIP-seq reveals that SPT6L is required to maintain the amount and length of Pol V transcripts. Our findings thus uncover the critical role of a Pol II conserved elongator in Pol V mediated DNA methylation and transcription, and shed light on the mutual regulation between Pol V and II in plants.

Chromosome fusion and programmed DNA elimination shape karyotypes of nematodes

An increasing number of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and increased numbers of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the nematode Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among other nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these nematodes, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes. These karyotype changes may lead to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms are dynamic and may play a role during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

Heterochromatic 3D genome organization is directed by HP1a- and H3K9-dependent and independent mechanisms

Whether and how histone post-translational modifications and the proteins that bind them drive 3D genome organization remains unanswered. Here, we evaluate the contribution of H3K9-methylated constitutive heterochromatin to 3D genome organization in Drosophila tissues. We find that the predominant organizational feature of wild-type tissues is the segregation of euchromatic chromosome arms from heterochromatic pericentromeres. Reciprocal perturbation of HP1a⋅H3K9me binding, using a point mutation in the HP1a chromodomain or replacement of the replication-dependent histone H3 with H3K9R mutant histones, revealed that HP1a binding to methylated H3K9 in constitutive heterochromatin is required to limit contact frequency between pericentromeres and chromosome arms and regulate the distance between arm and pericentromeric regions. Surprisingly, the self-association of pericentromeric regions is largely preserved despite the loss of H3K9 methylation and HP1a occupancy. Thus, the HP1a⋅H3K9 interaction contributes to but does not solely drive the segregation of euchromatin and heterochromatin inside the nucleus.

Rapid metabolic regulation of a novel arginine methylation of KCa3.1 attenuates T cell exhaustion

T cell exhaustion is linked to persistent antigen exposure and perturbed activation events, correlating with poor disease prognosis. Tumor-mediated T cell exhaustion is well documented; however, how the nutrient-deprived tumor niche affects T cell receptor (TCR) activation is largely unclear. We show that methionine metabolism licenses optimal TCR signaling by regulating the protein arginine methylome, and limiting methionine availability during early TCR signaling promotes subsequent T cell exhaustion. We discovered a novel arginine methylation of a Ca 2+ -activated potassium transporter, KCa3.1, prevention of which results in increased Ca 2+ -mediated NFAT1 activation, NFAT1 promoter occupancy, and T cell exhaustion. Furthermore, methionine supplementation reduces nuclear NFAT1 in tumor-infiltrating T cells and augments their anti-tumor activity. These findings demonstrate metabolic regulation of T cell exhaustion determined during TCR engagement.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Small-molecule-induced epigenetic rejuvenation promotes SREBP condensation and overcomes barriers to CNS myelin regeneration

Remyelination failure in diseases like multiple sclerosis (MS) was thought to involve suppressed maturation of oligodendrocyte precursors; however, oligodendrocytes are present in MS lesions yet lack myelin production. We found that oligodendrocytes in the lesions are epigenetically silenced. Developing a transgenic reporter labeling differentiated oligodendrocytes for phenotypic screening, we identified a small-molecule epigenetic-silencing-inhibitor (ESI1) that enhances myelin production and ensheathment. ESI1 promotes remyelination in animal models of demyelination and enables de novo myelinogenesis on regenerated CNS axons. ESI1 treatment lengthened myelin sheaths in human iPSC-derived organoids and augmented (re)myelination in aged mice while reversing age-related cognitive decline. Multi-omics revealed that ESI1 induces an active chromatin landscape that activates myelinogenic pathways and reprograms metabolism. Notably, ESI1 triggered nuclear condensate formation of master lipid-metabolic regulators SREBP1/2, concentrating transcriptional co-activators to drive lipid/cholesterol biosynthesis. Our study highlights the potential of targeting epigenetic silencing to enable CNS myelin regeneration in demyelinating diseases and aging.

The emerging H3K9me3 chromatin landscape during zebrafish embryogenesis

The structural organization of eukaryotic genomes is contingent upon the fractionation of DNA into transcriptionally permissive euchromatin and repressive heterochromatin. However, we have a limited understanding of how these distinct states are first established during animal embryogenesis. Histone 3 lysine 9 trimethylation (H3K9me3) is critical to heterochromatin formation and bulk establishment of this mark is thought to help drive large-scale remodeling of an initially naive chromatin state during animal embryogenesis. However, a detailed understanding of this process is lacking. Here, we leverage CUT&RUN to define the emerging H3K9me3 landscape of the zebrafish embryo with high sensitivity and temporal resolution. Despite the prevalence of DNA transposons in the zebrafish genome, we found that LTR transposons are preferentially targeted for embryonic H3K9me3 deposition, with different families exhibiting distinct establishment timelines. High signal-to-noise ratios afforded by CUT&RUN revealed new, emerging sites of low-amplitude H3K9me3 that initiated before the major wave of zygotic genome activation (ZGA). Early sites of establishment predominated at specific subsets of transposons and were particularly enriched for transposon sequences with maternal piRNAs and pericentromeric localization. Notably, the number of H3K9me3 enriched sites increased linearly across blastula development, while quantitative comparison revealed a >10-fold genome-wide increase in H3K9me3 signal at established sites over just 30 minutes at the onset of ZGA. Continued maturation of the H3K9me3 landscape was observed beyond the initial wave of bulk establishment.

Article Summary

Eukaryotic genomes are organized into transcriptionally permissive euchromatin and repressive heterochromatin. However, we have a limited understanding of how these distinct states are first established during animal embryogenesis. Here, we define the emerging landscape of the heterochromatic histone tail modification H3K9me3 across zebrafish embryogenesis. Our results uncover distinct kinetics of H3K9me3 establishment and growth, define unique characteristics associated with initial sites of H3K9me3, and demonstrate additional maturation of the H3K9me3 landscape beyond the window of bulk H3K9me3 deposition.

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.

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.

greenPipes: an integrated data analysis pipeline for greenCUT&RUN and CUT&RUN genome-localization datasets

MOTIVATION

To study gene regulation through transcription factors and chromatin modifiers, a variety of genome-wide techniques are used. Recently, CUT&RUN-based technologies have become popular, but a pipeline for the comprehensive analysis of CUT&RUN datasets is currently lacking. Here, we present the "greenPipes" package, which includes fine-tuned parameters specifically for bioinformatic analyses of greenCUT&RUN and CUT&RUN datasets. greenPipes provides additional functionalities for data analysis and data integration with other -omics technologies, which are either not available in other pipelines developed for CUT&RUN datasets or scattered in the literature as individual packages.

AVAILABILITY AND IMPLEMENTATION

Source code and a manual of the greenPipes are freely available on GitHub website (https://github.com/snizam001/greenPipes). The test datasets, comprehensive annotation files, and other datasets are available at https://osf.io/ruhj9/.

CONTACT

n.sheikh@dkfz-heidelberg.de or m.timmers@dkfz-heidelberg.de.

SUPPLEMENTARY INFORMATION

The handbook of greenPipes is available online at Bioinformatics as Supplementary text.

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.

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.

The MYCN oncoprotein is an RNA-binding accessory factor of the nuclear exosome targeting complex

The MYCN oncoprotein binds active promoters in a heterodimer with its partner protein MAX. MYCN also interacts with the nuclear exosome, a 3'-5' exoribonuclease complex, suggesting a function in RNA metabolism. Here, we show that MYCN forms stable high-molecular-weight complexes with the exosome and multiple RNA-binding proteins. MYCN binds RNA in vitro and in cells via a conserved sequence termed MYCBoxI. In cells, MYCN associates with thousands of intronic transcripts together with the ZCCHC8 subunit of the nuclear exosome targeting complex and enhances their processing. Perturbing exosome function results in global re-localization of MYCN from promoters to intronic RNAs. On chromatin, MYCN is then replaced by the MNT(MXD6) repressor protein, inhibiting MYCN-dependent transcription. RNA-binding-deficient alleles show that RNA-binding limits MYCN's ability to activate cell growth-related genes but is required for MYCN's ability to promote progression through S phase and enhance the stress resilience of neuroblastoma cells.

The transcription factor ZNF469 regulates collagen production in liver fibrosis

Non-alcoholic fatty liver disease (NAFLD) - characterized by excess accumulation of fat in the liver - now affects one third of the world's population. As NAFLD progresses, extracellular matrix components including collagen accumulate in the liver causing tissue fibrosis, a major determinant of disease severity and mortality. To identify transcriptional regulators of fibrosis, we computationally inferred the activity of transcription factors (TFs) relevant to fibrosis by profiling the matched transcriptomes and epigenomes of 108 human liver biopsies from a deeply-characterized cohort of patients spanning the full histopathologic spectrum of NAFLD. CRISPR-based genetic knockout of the top 100 TFs identified ZNF469 as a regulator of collagen expression in primary human hepatic stellate cells (HSCs). Gain- and loss-of-function studies established that ZNF469 regulates collagen genes and genes involved in matrix homeostasis through direct binding to gene bodies and regulatory elements. By integrating multiomic large-scale profiling of human biopsies with extensive experimental validation we demonstrate that ZNF469 is a transcriptional regulator of collagen in HSCs. Overall, these data nominate ZNF469 as a previously unrecognized determinant of NAFLD-associated liver fibrosis.

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.

Iterative oxidation by TET1 is required for reprogramming of imprinting control regions and patterning of mouse sperm hypomethylated regions

Ten-eleven translocation (TET) enzymes iteratively oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxylcytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during mammalian germline reprogramming remains unresolved due to the inability to decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 (Tet1-HxD) and TET1 that stalls oxidation at 5hmC (Tet1-V). Tet1 knockout and catalytic mutant primordial germ cells (PGCs) fail to erase methylation at select imprinting control regions and promoters of meiosis-associated genes, validating the requirement for the iterative oxidation of 5mC for complete germline reprogramming. TET1V and TET1HxD rescue most hypermethylation of Tet1-/- sperm, suggesting the role of TET1 beyond its oxidative capability. We additionally identify a broader class of hypermethylated regions in Tet1 mutant mouse sperm that depend on TET oxidation for reprogramming. Our study demonstrates the link between TET1-mediated germline reprogramming and sperm methylome patterning.

Single-genome analysis reveals a heterogeneous association of the herpes simplex virus genome with H3K27me2 and the reader PHF20L1 following infection of human fibroblasts

The fate of herpesvirus genomes following entry into different cell types is thought to regulate the outcome of infection. For the Herpes simplex virus 1 (HSV-1), latent infection of neurons is characterized by association with repressive heterochromatin marked with Polycomb silencing-associated lysine 27 methylation on histone H3 (H3K27me). However, whether H3K27 methylation plays a role in repressing lytic gene expression in non-neuronal cells is unclear. To address this gap in knowledge, and with consideration that the fate of the viral genome and outcome of HSV-1 infection could be heterogeneous, we developed an assay to quantify the abundance of histone modifications within single viral genome foci of infected fibroblasts. Using this approach, combined with bulk epigenetic techniques, we were unable to detect any role for H3K27me3 during HSV-1 lytic infection of fibroblasts. By contrast, we could detect the lesser studied H3K27me2 on a subpopulation of viral genomes, which was consistent with a role for H3K27 demethylases in promoting lytic gene expression. In addition, viral genomes co-localized with the H3K27me2 reader protein PHF20L1, and this association was enhanced by inhibition of the H3K27 demethylases UTX and JMJD3. Notably, targeting of H3K27me2 to viral genomes was enhanced following infection with a transcriptionally defective virus in the absence of Promyelocytic leukemia nuclear bodies. Collectively, these studies implicate a role for H3K27me2 in fibroblast-associated HSV genome silencing in a manner dependent on genome sub-nuclear localization and transcriptional activity.

IMPORTANCE

Investigating the potential mechanisms of gene silencing for DNA viruses in different cell types is important to understand the differential outcomes of infection, particularly for viruses like herpesviruses that can undergo distinct types of infection in different cell types. In addition, investigating chromatin association with viral genomes informs on the mechanisms of epigenetic regulation of DNA processes. However, there is a growing appreciation for heterogeneity in the outcome of infection at the single cell, and even single viral genome, level. Here we describe a novel assay for quantifying viral genome foci with chromatin proteins and show that a portion of genomes are targeted for silencing by H3K27me2 and associate with the reader protein PHF20L1. This study raises important questions regarding the mechanism of H3K27me2-specific targeting to viral genomes, the contribution of epigenetic heterogeneity to herpesvirus infection, and the role of PHF20L1 in regulating the outcome of DNA virus infection.

The SAGA acetyltransferase module is required for the maintenance of MAF and MYC oncogenic gene expression programs in multiple myeloma

Despite recent advances in therapeutic treatments, multiple myeloma (MM) remains an incurable malignancy. Epigenetic factors contribute to the initiation, progression, relapse, and clonal heterogeneity in MM, but our knowledge on epigenetic mechanisms underlying MM development is far from complete. The SAGA complex serves as a coactivator in transcription and catalyzes acetylation and deubiquitylation. Analyses of datasets in the Cancer Dependency Map Project revealed many SAGA components are selective dependencies in MM. To define SAGA-specific functions, we focused on ADA2B, the only subunit in the lysine acetyltransferase (KAT) module that specifically functions in SAGA. Integration of RNA-seq, ATAC-seq, and CUT&RUN results identified pathways directly regulated by ADA2B include MTORC1 signaling, MYC, E2F, and MM-specific MAF oncogenic programs. We discovered that ADA2B is recruited to MAF and MYC gene targets, and that MAF shares a majority of its targets with MYC in MM cells. Furthermore, we found the SANT domain of ADA2B is required for interaction with both GCN5 and PCAF acetyltransferases, incorporation into SAGA, and ADA2B protein stability. Our findings uncover previously unknown SAGA KAT module-dependent mechanisms controlling MM cell growth, revealing a vulnerability that might be exploited for future development of MM therapy.

Chaetocin-mediated SUV39H1 inhibition targets stemness and oncogenic networks of diffuse midline gliomas and synergizes with ONC201

BACKGROUND

Diffuse intrinsic pontine gliomas (DIPG/DMG) are devastating pediatric brain tumors with extraordinarily limited treatment options and uniformly fatal prognosis. Histone H3K27M mutation is a common recurrent alteration in DIPG and disrupts epigenetic regulation. We hypothesize that genome-wide H3K27M-induced epigenetic dysregulation makes tumors vulnerable to epigenetic targeting.

METHODS

We performed a screen of compounds targeting epigenetic enzymes to identify potential inhibitors for the growth of patient-derived DIPG cells. We further carried out transcriptomic and genomic landscape profiling including RNA-seq and CUT&RUN-seq as well as shRNA-mediated knockdown to assess the effects of chaetocin and SUV39H1, a target of chaetocin, on DIPG growth.

RESULTS

High-throughput small-molecule screening identified an epigenetic compound chaetocin as a potent blocker of DIPG cell growth. Chaetocin treatment selectively decreased proliferation and increased apoptosis of DIPG cells and significantly extended survival in DIPG xenograft models, while restoring H3K27me3 levels. Moreover, the loss of H3K9 methyltransferase SUV39H1 inhibited DIPG cell growth. Transcriptomic and epigenomic profiling indicated that SUV39H1 loss or inhibition led to the downregulation of stemness and oncogenic networks including growth factor receptor signaling and stemness-related programs; however, D2 dopamine receptor (DRD2) signaling adaptively underwent compensatory upregulation conferring resistance. Consistently, a combination of chaetocin treatment with a DRD2 antagonist ONC201 synergistically increased the antitumor efficacy.

CONCLUSIONS

Our studies reveal a therapeutic vulnerability of DIPG cells through targeting the SUV39H1-H3K9me3 pathway and compensatory signaling loops for treating this devastating disease. Combining SUV39H1-targeting chaetocin with other agents such as ONC201 may offer a new strategy for effective DIPG treatment.

Loss of CREBBP and KMT2D cooperate to accelerate lymphomagenesis and shape the lymphoma immune microenvironment

Despite regulating overlapping gene enhancers and pathways, CREBBP and KMT2D mutations recurrently co-occur in germinal center (GC) B cell-derived lymphomas, suggesting potential oncogenic cooperation. Herein, we report that combined haploinsufficiency of Crebbp and Kmt2d induces a more severe mouse lymphoma phenotype (vs either allele alone) and unexpectedly confers an immune evasive microenvironment manifesting as CD8+ T-cell exhaustion and reduced infiltration. This is linked to profound repression of immune synapse genes that mediate crosstalk with T-cells, resulting in aberrant GC B cell fate decisions. From the epigenetic perspective, we observe interaction and mutually dependent binding and function of CREBBP and KMT2D on chromatin. Their combined deficiency preferentially impairs activation of immune synapse-responsive super-enhancers, pointing to a particular dependency for both co-activators at these specialized regulatory elements. Together, our data provide an example where chromatin modifier mutations cooperatively shape and induce an immune-evasive microenvironment to facilitate lymphomagenesis.

A GREB1-steroid receptor feedforward mechanism governs differential GREB1 action in endometrial function and endometriosis

Cellular responses to the steroid hormones, estrogen (E2), and progesterone (P4) are governed by their cognate receptor's transcriptional output. However, the feed-forward mechanisms that shape cell-type-specific transcriptional fulcrums for steroid receptors are unidentified. Herein, we found that a common feed-forward mechanism between GREB1 and steroid receptors regulates the differential effect of GREB1 on steroid hormones in a physiological or pathological context. In physiological (receptive) endometrium, GREB1 controls P4-responses in uterine stroma, affecting endometrial receptivity and decidualization, while not affecting E2-mediated epithelial proliferation. Of mechanism, progesterone-induced GREB1 physically interacts with the progesterone receptor, acting as a cofactor in a positive feedback mechanism to regulate P4-responsive genes. Conversely, in endometrial pathology (endometriosis), E2-induced GREB1 modulates E2-dependent gene expression to promote the growth of endometriotic lesions in mice. This differential action of GREB1 exerted by a common feed-forward mechanism with steroid receptors advances our understanding of mechanisms that underlie cell- and tissue-specific steroid hormone actions.

ETS transcription factors regulate precise matrix metalloproteinase expression and follicle rupture in Drosophila

Drosophila matrix metalloproteinase 2 (MMP2) is specifically expressed in posterior follicle cells of stage-14 egg chambers (mature follicles) and is crucial for the breakdown of the follicular wall during ovulation, a process that is highly conserved from flies to mammals. The factors that regulate spatiotemporal expression of MMP2 in follicle cells remain unknown. Here, we demonstrate crucial roles for the ETS-family transcriptional activator Pointed (Pnt) and its endogenous repressor Yan in the regulation of MMP2 expression. We found that Pnt is expressed in posterior follicle cells and overlaps with MMP2 expression in mature follicles. Genetic analysis demonstrated that pnt is both required and sufficient for MMP2 expression in follicle cells. In addition, Yan was temporally upregulated in stage-13 follicle cells to fine-tune Pnt activity and MMP2 expression. Furthermore, we identified a 1.1 kb core enhancer that is responsible for the spatiotemporal expression of MMP2 and contains multiple pnt/yan binding motifs. Mutation of pnt/yan binding sites significantly impaired the Mmp2 enhancer activity. Our data reveal a mechanism of transcriptional regulation of Mmp2 expression in Drosophila ovulation, which could be conserved in other biological systems.

INO80 regulates chromatin accessibility to facilitate suppression of sex-linked gene expression during mouse spermatogenesis

The INO80 protein is the main catalytic subunit of the INO80-chromatin remodeling complex, which is critical for DNA repair and transcription regulation in murine spermatocytes. In this study, we explored the role of INO80 in silencing genes on meiotic sex chromosomes in male mice. INO80 immunolocalization at the XY body in pachytene spermatocytes suggested a role for INO80 in the meiotic sex body. Subsequent deletion of Ino80 resulted in high expression of sex-linked genes. Furthermore, the active form of RNA polymerase II at the sex body of Ino80 -null pachytene spermatocytes indicates incomplete inactivation of sex-linked genes. A reduction in the recruitment of initiators of meiotic sex chromosome inhibition (MSCI) argues for INO80-facilitated recruitment of DNA repair factors required for silencing sex-linked genes. This role of INO80 is independent of a common INO80 target H2A.Z. Instead, in the absence of INO80, a reduction in chromatin accessibility at DNA repair sites occurs on the sex chromosomes. These data suggest a role for INO80 in DNA repair factor localization, thereby facilitating the silencing of sex-linked genes during the onset of pachynema.

Summary Statement

Chromatin accessibility and DNA repair factor localization at the sex chromosomes are facilitated by INO80, which regulates sex-linked gene silencing during meiotic progression in spermatocytes.

Selection of epigenetically privileged HIV-1 proviruses during treatment with panobinostat and interferon-α2a

CD4+ T cells with latent HIV-1 infection persist despite treatment with antiretroviral agents and represent the main barrier to a cure of HIV-1 infection. Pharmacological disruption of viral latency may expose HIV-1-infected cells to host immune activity, but the clinical efficacy of latency-reversing agents for reducing HIV-1 persistence remains to be proven. Here, we show in a randomized-controlled human clinical trial that the histone deacetylase inhibitor panobinostat, when administered in combination with pegylated interferon-α2a, induces a structural transformation of the HIV-1 reservoir cell pool, characterized by a disproportionate overrepresentation of HIV-1 proviruses integrated in ZNF genes and in chromatin regions with reduced H3K27ac marks, the molecular target sites for panobinostat. By contrast, proviruses near H3K27ac marks were actively selected against, likely due to increased susceptibility to panobinostat. These data suggest that latency-reversing treatment can increase the immunological vulnerability of HIV-1 reservoir cells and accelerate the selection of epigenetically privileged HIV-1 proviruses.

Cellular reprogramming in vivo initiated by SOX4 pioneer factor activity

Tissue damage elicits cell fate switching through a process called metaplasia, but how the starting cell fate is silenced and the new cell fate is activated has not been investigated in animals. In cell culture, pioneer transcription factors mediate "reprogramming" by opening new chromatin sites for expression that can attract transcription factors from the starting cell's enhancers. Here we report that SOX4 is sufficient to initiate hepatobiliary metaplasia in the adult mouse liver, closely mimicking metaplasia initiated by toxic damage to the liver. In lineage-traced cells, we assessed the timing of SOX4-mediated opening of enhancer chromatin versus enhancer decommissioning. Initially, SOX4 directly binds to and closes hepatocyte regulatory sequences via an overlapping motif with HNF4A, a hepatocyte master regulatory transcription factor. Subsequently, SOX4 exerts pioneer factor activity to open biliary regulatory sequences. The results delineate a hierarchy by which gene networks become reprogrammed under physiological conditions, providing deeper insight into the basis for cell fate transitions in animals.

PAF1c links S-phase progression to immune evasion and MYC function in pancreatic carcinoma

In pancreatic ductal adenocarcinoma (PDAC), endogenous MYC is required for S-phase progression and escape from immune surveillance. Here we show that MYC in PDAC cells is needed for the recruitment of the PAF1c transcription elongation complex to RNA polymerase and that depletion of CTR9, a PAF1c subunit, enables long-term survival of PDAC-bearing mice. PAF1c is largely dispensable for normal proliferation and regulation of MYC target genes. Instead, PAF1c limits DNA damage associated with S-phase progression by being essential for the expression of long genes involved in replication and DNA repair. Surprisingly, the survival benefit conferred by CTR9 depletion is not due to DNA damage, but to T-cell activation and restoration of immune surveillance. This is because CTR9 depletion releases RNA polymerase and elongation factors from the body of long genes and promotes the transcription of short genes, including MHC class I genes. The data argue that functionally distinct gene sets compete for elongation factors and directly link MYC-driven S-phase progression to tumor immune evasion.

Adipocyte p53 coordinates the response to intermittent fasting by regulating adipose tissue immune cell landscape

In obesity, sustained adipose tissue (AT) inflammation constitutes a cellular memory that limits the effectiveness of weight loss interventions. Yet, the impact of fasting regimens on the regulation of AT immune infiltration is still elusive. Here we show that intermittent fasting (IF) exacerbates the lipid-associated macrophage (LAM) inflammatory phenotype of visceral AT in obese mice. Importantly, this increase in LAM abundance is strongly p53 dependent and partly mediated by p53-driven adipocyte apoptosis. Adipocyte-specific deletion of p53 prevents LAM accumulation during IF, increases the catabolic state of adipocytes, and enhances systemic metabolic flexibility and insulin sensitivity. Finally, in cohorts of obese/diabetic patients, we describe a p53 polymorphism that links to efficacy of a fasting-mimicking diet and that the expression of p53 and TREM2 in AT negatively correlates with maintaining weight loss after bariatric surgery. Overall, our results demonstrate that p53 signalling in adipocytes dictates LAM accumulation in AT under IF and modulates fasting effectiveness in mice and humans.

Matrin3 mediates differentiation through stabilizing chromatin loop-domain interactions and YY1 mediated enhancer-promoter interactions

Although emerging evidence indicates that alterations in proteins within nuclear compartments elicit changes in chromosomal architecture and differentiation, the underlying mechanisms are not well understood. Here we investigate the direct role of the abundant nuclear complex protein Matrin3 (Matr3) in chromatin architecture and development in the context of myogenesis. Using an acute targeted protein degradation platform (dTAG-Matr3), we reveal the dynamics of development-related chromatin reorganization. High-throughput chromosome conformation capture (Hi-C) experiments revealed substantial chromatin loop rearrangements soon after Matr3 depletion. Notably, YY1 binding was detected, accompanied by the emergence of novel YY1-mediated enhancer-promoter loops, which occurred concurrently with changes in histone modifications and chromatin-level binding patterns. Changes in chromatin occupancy by Matr3 also correlated with these alterations. Overall, our results suggest that Matr3 mediates differentiation through stabilizing chromatin accessibility and chromatin loop-domain interactions, and highlight a conserved and direct role for Matr3 in maintenance of chromosomal architecture.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

The transcriptional cofactor Tle3 reciprocally controls effector and central memory CD8+ T cell fates

Antigen-experienced CD8+ T cells form effector and central memory T cells (TEM and TCM cells, respectively); however, the mechanism(s) controlling their lineage plasticity remains incompletely understood. Here we show that the transcription cofactor Tle3 critically regulates TEM and TCM cell fates and lineage stability through dynamic redistribution in antigen-responding CD8+ T cell genome. Genetic ablation of Tle3 promoted CD8+ TCM cell formation at the expense of CD8+ TEM cells. Lineage tracing showed that Tle3-deficient CD8+ TEM cells underwent accelerated conversion into CD8+ TCM cells while retaining robust recall capacity. Tle3 acted as a coactivator for Tbet to increase chromatin opening at CD8+ TEM cell-characteristic sites and to activate CD8+ TEM cell signature gene transcription, while engaging Runx3 and Tcf1 to limit CD8+ TCM cell-characteristic molecular features. Thus, Tle3 integrated functions of multiple transcription factors to guard lineage fidelity of CD8+ TEM cells, and manipulation of Tle3 activity could favor CD8+ TCM cell production.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and 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.

Inhibition of the YAP-MMB interaction and targeting NEK2 as potential therapeutic strategies for YAP-driven cancers

YAP activation in cancer is linked to poor outcomes, making it an attractive therapeutic target. Previous research focused on blocking the interaction of YAP with TEAD transcription factors. Here, we took a different approach by disrupting YAP's binding to the transcription factor B-MYB using MY-COMP, a fragment of B-MYB containing the YAP binding domain fused to a nuclear localization signal. MY-COMP induced cell cycle defects, nuclear abnormalities, and polyploidization. In an AKT and YAP-driven liver cancer model, MY-COMP significantly reduced liver tumorigenesis, highlighting the importance of the YAP-B-MYB interaction in tumor development. MY-COMP also perturbed the cell cycle progression of YAP-dependent uveal melanoma cells but not of YAP-independent cutaneous melanoma cell lines. It counteracted YAP-dependent expression of MMB-regulated cell cycle genes, explaining the observed effects. We also identified NIMA-related kinase (NEK2) as a downstream target of YAP and B-MYB, promoting YAP-driven transformation by facilitating centrosome clustering and inhibiting multipolar mitosis.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

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.

Supergene evolution via gain of auto-regulation

The development of complex phenotypes requires the coordinated action of many genes across space and time, yet many species have evolved the ability to develop multiple discrete, alternate phenotypes1-3. Such polymorphisms are often controlled by supergenes, sets of tightly-linked mutations in one or more loci that function together to produce a complex phenotype4. Although theories of supergene evolution are well-established, the mutations that cause functional differences between supergene alleles remain essentially unknown. doublesex is the master regulator of insect sexual differentiation but functions as a supergene in multiple Papilio swallowtail butterflies, where divergent dsx alleles control development of discrete non-mimetic or mimetic female wing color patterns5-7. Here we demonstrate that the functional elements of the mimetic allele in Papilio alphenor are six new cis-regulatory elements (CREs) spread across 150 kb that are bound by DSX itself. Our findings provide experimental support to classic supergene theory and suggest that the evolution of auto-regulation may provide a simple route to supergene origination and to the co-option of pleiotropic genes into new developmental roles.

Histone lactylation couples cellular metabolism with developmental gene regulatory networks

Embryonic cells exhibit diverse metabolic states. Recent studies have demonstrated that metabolic reprogramming drives changes in cell identity by affecting gene expression. However, the connection between cellular metabolism and gene expression remains poorly understood. Here we report that glycolysis-regulated histone lactylation couples the metabolic state of embryonic cells with chromatin organization and gene regulatory network (GRN) activation. We found that lactylation marks genomic regions of glycolytic embryonic tissues, like the neural crest (NC) and pre-somitic mesoderm. Histone lactylation occurs in the loci of NC genes as these cells upregulate glycolysis. This process promotes the accessibility of active enhancers and the deployment of the NC GRN. Reducing the deposition of the mark by targeting LDHA/B leads to the downregulation of NC genes and the impairment of cell migration. The deposition of lactyl-CoA on histones at NC enhancers is supported by a mechanism that involves transcription factors SOX9 and YAP/TEAD. These findings define an epigenetic mechanism that integrates cellular metabolism with the GRNs that orchestrate embryonic development.

Combined small-molecule treatment accelerates maturation of human pluripotent stem cell-derived neurons

The maturation of human pluripotent stem cell (hPSC)-derived neurons mimics the protracted timing of human brain development, extending over months to years for reaching adult-like function. Prolonged in vitro maturation presents a major challenge to stem cell-based applications in modeling and treating neurological disease. Therefore, we designed a high-content imaging assay based on morphological and functional readouts in hPSC-derived cortical neurons which identified multiple compounds that drive neuronal maturation including inhibitors of lysine-specific demethylase 1 and disruptor of telomerase-like 1 and activators of calcium-dependent transcription. A cocktail of four factors, GSK2879552, EPZ-5676, N-methyl-D-aspartate and Bay K 8644, collectively termed GENtoniK, triggered maturation across all parameters tested, including synaptic density, electrophysiology and transcriptomics. Maturation effects were further validated in cortical organoids, spinal motoneurons and non-neural lineages including melanocytes and pancreatic β-cells. The effects on maturation observed across a broad range of hPSC-derived cell types indicate that some of the mechanisms controlling the timing of human maturation might be shared across lineages.

MICA: a multi-omics method to predict gene regulatory networks in early human embryos

Recent advances in single-cell omics have transformed characterisation of cell types in challenging-to-study biological contexts. In contexts with limited single-cell samples, such as the early human embryo inference of transcription factor-gene regulatory network (GRN) interactions is especially difficult. Here, we assessed application of different linear or non-linear GRN predictions to single-cell simulated and human embryo transcriptome datasets. We also compared how expression normalisation impacts on GRN predictions, finding that transcripts per million reads outperformed alternative methods. GRN inferences were more reproducible using a non-linear method based on mutual information (MI) applied to single-cell transcriptome datasets refined with chromatin accessibility (CA) (called MICA), compared with alternative network prediction methods tested. MICA captures complex non-monotonic dependencies and feedback loops. Using MICA, we generated the first GRN inferences in early human development. MICA predicted co-localisation of the AP-1 transcription factor subunit proto-oncogene JUND and the TFAP2C transcription factor AP-2γ in early human embryos. Overall, our comparative analysis of GRN prediction methods defines a pipeline that can be applied to single-cell multi-omics datasets in especially challenging contexts to infer interactions between transcription factor expression and target gene regulation.

Nuclear growth and import can be uncoupled

What drives nuclear growth? Studying nuclei assembled in Xenopus egg extract and focusing on importin α/β-mediated nuclear import, we show that, while import is required for nuclear growth, nuclear growth and import can be uncoupled when chromatin structure is manipulated. Nuclei treated with micrococcal nuclease to fragment DNA grew slowly despite exhibiting little to no change in import rates. Nuclei assembled around axolotl chromatin with 20-fold more DNA than Xenopus grew larger but imported more slowly. Treating nuclei with reagents known to alter histone methylation or acetylation caused nuclei to grow less while still importing to a similar extent or to grow larger without significantly increasing import. Nuclear growth but not import was increased in live sea urchin embryos treated with the DNA methylator N-nitrosodimethylamine. These data suggest that nuclear import is not the primary driving force for nuclear growth. Instead, we observed that nuclear blebs expanded preferentially at sites of high chromatin density and lamin addition, whereas small Benzonase-treated nuclei lacking DNA exhibited reduced lamin incorporation into the nuclear envelope. In summary, we report experimental conditions where nuclear import is not sufficient to drive nuclear growth, hypothesizing that this uncoupling is a result of altered chromatin structure.

In Vitro Determination of Temperature-Dependent DNA Binding of the Evening Complex Using Electrophoretic Mobility Shift Assays

Electrophoretic mobility shift assays (EMSAs) of DNA-binding proteins and labeled DNA allow the qualitative and quantitative characterization of protein-DNA complex formation using native (nondenaturing) polyacrylamide or agarose gel electrophoresis. By varying the incubation temperature of the protein-DNA binding reaction and maintaining this temperature during electrophoresis, temperature-dependent protein-DNA interactions can be investigated. Here, we provide examples of the binding of a transcriptional repressor complex called the Evening Complex, comprising the DNA-binding protein LUX ARRYTHMO (LUX), the scaffold protein EARLY FLOWERING 3 (ELF3), and the adapter protein ELF4, to its cognate DNA and demonstrate direct detection and visualization of thermoresponsive binding in vitro. As negative controls we use the LUX DNA-binding domain and LUX full length protein, which do not exhibit temperature-dependent DNA binding.

Cleavage Under Targets & Release Using Nuclease (CUT&RUN) of Histone Modifications in Epidermal Stem Cells of Adult Murine Skin

Cleavage Under Targets & Release Using Nuclease (CUT&RUN) has swiftly become the preferred procedure over the past few years for genomic mapping and detecting interactions between chromatin and its bound proteins. CUT&RUN is now being widely used for characterizing the epigenetic landscape in many cell types as it utilizes far less cell numbers when compared to Chromatin Immunoprecipitation-sequencing (ChIP-seq), thereby making it a powerful tool for researchers working with limited material. This protocol has been specifically optimized for detecting histone modifications in fluorescence-activated cell sorting (FACS)-isolated epidermal stem cells from adult mice.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataracts. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq and CUT&RUN-seq to discover previously unreported mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Furthermore, we identify an epigenetic paradigm of cellular differentiation, defined by progressive loss of the H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

Chromosome fusion and programmed DNA elimination shape karyotypes of parasitic nematodes

A growing list of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and an increased number of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the horse parasite Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among ascarid nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these parasites, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes, leading to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms play a dynamic role in the Parascaris germline chromosome during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

A histone deacetylase network regulates epigenetic reprogramming and viral silencing in HIV-infected cells

A long-lived latent reservoir of HIV-1-infected CD4 T cells persists with antiretroviral therapy and prevents cure. We report that the emergence of latently infected primary CD4 T cells requires the activity of histone deacetylase enzymes HDAC1/2 and HDAC3. Data from targeted HDAC molecules, an HDAC3-directed PROTAC, and CRISPR-Cas9 knockout experiments converge on a model where either HDAC1/2 or HDAC3 targeting can prevent latency, whereas all three enzymes must be targeted to achieve latency reversal. Furthermore, HDACi treatment targets features of memory T cells that are linked to proviral latency and persistence. Latency prevention is associated with increased H3K9ac at the proviral LTR promoter region and decreased H3K9me3, suggesting that this epigenetic switch is a key proviral silencing mechanism that depends on HDAC activity. These findings support further mechanistic work on latency initiation and eventual clinical studies of HDAC inhibitors to interfere with latency initiation.