Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci U S A. 2010;107:21931-21936. [PMID: 21106759 DOI: 10.1073/pnas.1016071107]

CRISPR Activation Reverses Haploinsufficiency and Functional Deficits Caused by TTN Truncation Variants

BACKGROUND

TTN truncation variants (TTNtvs) are the most common genetic lesion identified in individuals with dilated cardiomyopathy, a disease with high morbidity and mortality rates. TTNtvs reduce normal TTN (titin) protein levels, produce truncated proteins, and impair sarcomere content and function. Therapeutics targeting TTNtvs have been elusive because of the immense size of TTN, the rarity of specific TTNtvs, and incomplete knowledge of TTNtv pathogenicity.

METHODS

We adapted CRISPR activation using dCas9-VPR to functionally interrogate TTNtv pathogenicity and develop a therapeutic in human cardiomyocytes and 3-dimensional cardiac microtissues engineered from induced pluripotent stem cell models harboring a dilated cardiomyopathy-associated TTNtv. We performed guide RNA screening with custom TTN reporter assays, agarose gel electrophoresis to quantify TTN protein levels and isoforms, and RNA sequencing to identify molecular consequences of TTN activation. Cardiomyocyte epigenetic assays were also used to nominate DNA regulatory elements to enable cardiomyocyte-specific TTN activation.

RESULTS

CRISPR activation of TTN using single guide RNAs targeting either the TTN promoter or regulatory elements in spatial proximity to the TTN promoter through 3-dimensional chromatin interactions rescued TTN protein deficits disturbed by TTNtvs. Increasing TTN protein levels normalized sarcomere content and contractile function despite increasing truncated TTN protein. In addition to TTN transcripts, CRISPR activation also increased levels of myofibril assembly-related and sarcomere-related transcripts.

CONCLUSIONS

TTN CRISPR activation rescued TTNtv-related functional deficits despite increasing truncated TTN levels, which provides evidence to support haploinsufficiency as a relevant genetic mechanism underlying heterozygous TTNtvs. CRISPR activation could be developed as a therapeutic to treat a large proportion of TTNtvs.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

KDM5C-Mediated Recruitment of BRD4 to Chromatin Regulates Enhancer Activation and BET Inhibitor Sensitivity

The BET family member BRD4 is a bromodomain-containing protein that plays a vital role in driving oncogene expression. Given their pivotal role in regulating oncogenic networks in various cancer types, BET inhibitors (BETi) have been developed, but the clinical application has been impeded by dose-limiting toxicity and resistance. Understanding the mechanisms of BRD4 activity and identifying predictive biomarkers could facilitate the successful clinical use of BETis. Herein, we show that KDM5C and BRD4 cooperate to sustain tumor cell growth. Mechanistically, KDM5C interacted with BRD4 and stimulated BRD4 enhancer recruitment. Moreover, binding of the BRD4 C-terminus to KDM5C stimulated the H3K4 demethylase activity of KDM5C. The abundance of both KDM5C-associated BRD4 and H3K4me1/3 determined the transcriptional activation of many oncogenes. Notably, depletion or pharmacologic degradation of KDM5C dramatically reduced BRD4 chromatin enrichment and significantly increased BETi efficacy across multiple cancer types in both tumor cell lines and patient-derived organoid models. Furthermore, targeting KDM5C in combination with BETi suppressed tumor growth in vivo in a xenograft mouse model. Collectively, this work reveals a KDM5C-mediated mechanism by which BRD4 regulates transcription, providing a rationale for incorporating BETi into combination therapies with KDM5C inhibitors to enhance treatment efficacy.

SIGNIFICANCE

BRD4 is recruited to enhancers in a bromodomain-independent manner by binding KDM5C and stimulates KDM5C H3K4 demethylase activity, leading to synergistic effects of BET and KDM5C inhibitor combinations in cancer.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N = 29635, FDR < 0.01 and fold change > 50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Chromatin modifiers in human disease: from functional roles to regulatory mechanisms

The field of transcriptional regulation has revealed the vital role of chromatin modifiers in human diseases from the beginning of functional exploration to the process of participating in many types of disease regulatory mechanisms. Chromatin modifiers are a class of enzymes that can catalyze the chemical conversion of pyrimidine residues or amino acid residues, including histone modifiers, DNA methyltransferases, and chromatin remodeling complexes. Chromatin modifiers assist in the formation of transcriptional regulatory circuits between transcription factors, enhancers, and promoters by regulating chromatin accessibility and the ability of transcription factors to acquire DNA. This is achieved by recruiting associated proteins and RNA polymerases. They modify the physical contact between cis-regulatory factor elements, transcription factors, and chromatin DNA to influence transcriptional regulatory processes. Then, abnormal chromatin perturbations can impair the homeostasis of organs, tissues, and cells, leading to diseases. The review offers a comprehensive elucidation on the function and regulatory mechanism of chromatin modifiers, thereby highlighting their indispensability in the development of diseases. Furthermore, this underscores the potential of chromatin modifiers as biomarkers, which may enable early disease diagnosis. With the aid of this paper, a deeper understanding of the role of chromatin modifiers in the pathogenesis of diseases can be gained, which could help in devising effective diagnostic and therapeutic interventions.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

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.

Distinct Histone H3 Lysine 27 Modifications Dictate Different Outcomes of Gene Transcription

Site-specific histone modifications have long been recognized to play an important role in directing gene transcription in chromatin in biology of health and disease. However, concrete illustration of how different histone modifications in a site-specific manner dictate gene transcription outcomes, as postulated in the influential "Histone code hypothesis", introduced by Allis and colleagues in 2000, has been lacking. In this review, we summarize our latest understanding of the dynamic regulation of gene transcriptional activation, silence, and repression in chromatin that is directed distinctively by histone H3 lysine 27 acetylation, methylation, and crotonylation, respectively. This represents a special example of a long-anticipated verification of the "Histone code hypothesis."

dCas9 Tells Tales: Probing Gene Function and Transcription Regulation in Cancer

Clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing is evolving into an essential tool in the field of biological and medical research. Notably, the development of catalytically deactivated Cas9 (dCas9) enzyme has substantially broadened its traditional boundaries in gene editing or perturbation. The conjugation of dCas9 with various molecular effectors allows precise control over transcriptional processes, epigenetic modifications, visualization of chromosomal dynamics, and several other applications. This expanded repertoire of CRISPR-Cas9 applications has emerged as an invaluable molecular tool kit that empowers researchers to comprehensively interrogate and gain insights into health and diseases. This review delves into the advancements in Cas9 protein engineering, specifically on the generation of various dCas9 tools that have significantly enhanced the CRISPR-based technology capability and versatility. We subsequently discuss the multifaceted applications of dCas9, especially in interrogating the regulation and function of genes that involve in supporting cancer pathogenesis. In addition, we also delineate the designing and utilization of dCas9-based tools as well as highlighting its current constraints and transformative potentials in cancer research.

Nescient helix-loop-helix 1 (Nhlh1) is a novel activating transcription factor 5 (ATF5) target gene in olfactory and vomeronasal sensory neurons in mice

Activating transcription factor 5 (ATF5) is a transcription factor that belongs to the cAMP-response element-binding protein/ATF family and is essential for the differentiation and survival of sensory neurons in mouse olfactory organs. However, transcriptional target genes for ATF5 have yet to be identified. In the present study, chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) experiments were performed to verify ATF5 target genes in the main olfactory epithelium and vomeronasal organ in the postnatal pups. ChIP-qPCR was conducted using hemagglutinin (HA)-tagged ATF5 knock-in olfactory organs. The results obtained demonstrated that ATF5-HA fusion proteins bound to the CCAAT/enhancer-binding protein-ATF response element (CARE) site in the enhancer region of nescient helix-loop-helix 1 (Nhlh1), a transcription factor expressed in differentiating olfactory and vomeronasal sensory neurons. Nhlh1 mRNA expression was downregulated in ATF5-deficient (ATF5-/-) olfactory organs. The LIM/homeobox protein transcription factor Lhx2 co-localized with ATF5 in the nuclei of olfactory and vomeronasal sensory neurons and bound to the homeodomain site proximal to the CARE site in the Nhlh1 gene. The CARE region of the Nhlh1 gene was enriched by the active enhancer marker, acetyl-histone H3 (Lys27). The present study identified Nhlh1 as a novel target gene for ATF5 in murine olfactory organs. ATF5 may upregulate Nhlh1 expression in concert with Lhx2, thereby promoting the differentiation of olfactory and vomeronasal sensory neurons.

Electroacupuncture regulates histone acetylation of Bcl-2 and Caspase-3 genes to improve ischemic stroke injury

Background

Imbalances between Bcl-2 and caspase-3 are significant evidence of apoptosis, which is considered an influential factor in rapidly occurring neuronal cell death and the decline of neurological function after stroke. Studies have shown that acupuncture can reduce poststroke brain cell damage via either an increase in Bcl-2 or a reduction in caspase-3 exposure. The current study aimed to investigate whether acupuncture could modulate Bcl-2 and caspase-3 expression through histone acetylation modifications, which could potentially serve as a neuroprotective mechanism.

Methods

This study used TTC staining, Nissl staining, Clark neurological system score, and Evans Blue (EB) extravasation to evaluate neurological damage following stroke. The expression of Bcl-2/caspase-3 mRNA was detected by real-time fluorescence quantification of PCR (real-time PCR), whereas the protein expression levels of Bcl-2, Bax, caspase-3, and cleaved caspase-3 were assessed using western blotting. TUNEL staining of the ischemic cortical neurons determined apoptosis in the ischemic cortex. Histone acetyltransferase (HAT) and histone deacetylase (HDAC) activities, along with the protein performance of AceH3, H3K9ace, and H3K27ace, were detected to evaluate the degree of histone acetylation. The acetylation enrichment levels of H3K9 and K3K27 in the Bcl-2/caspase-3 gene were assessed using Chromatin Immunoprecipitation (ChIP) assay.

Results

Our data demonstrated that electroacupuncture (EA) exerts a significant neuroprotective effect in middle cerebral artery occlusion (MCAO) rats, as evidenced by a reduction in infarct volume, neuronal damage, Blood-Brain Barrier (BBB) disruption, and decreased apoptosis of ischemic cortical neurons. EA treatment can promote the mRNA and protein expression of the Bcl-2 gene in the ischemic brain while reducing the mRNA and protein expression levels of caspase-3 and effectively decreasing the protein expression levels of Bax and cleaved caspase-3. More importantly, EA treatment enhanced the level of histone acetylation, including Ace-H3, H3K9ace, and H3K27ace, significantly enhanced the occupancy of H3K9ace/H3K27ace at the Bcl-2 promoter, and reduced the enrichment of H3K9ace and H3K27ace at the caspase-3 promoter. However, the Histone Acetyltransferase inhibitor (HATi) treatment reversed these effects.

Conclusions

Our data demonstrated that EA mediated the expression levels of Bcl-2 and caspase-3 in MCAO rats by regulating the occupancy of acetylated H3K9/H3K27 at the promoters of these two genes, thus exerting a cerebral protective effect in ischemic reperfusion (I/R) injury.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Alcohol-Associated Liver Disease Outcomes: Critical Mechanisms of Liver Injury Progression

Alcohol-associated liver disease (ALD) is a substantial cause of morbidity and mortality worldwide and represents a spectrum of liver injury beginning with hepatic steatosis (fatty liver) progressing to inflammation and culminating in cirrhosis. Multiple factors contribute to ALD progression and disease severity. Here, we overview several crucial mechanisms related to ALD end-stage outcome development, such as epigenetic changes, cell death, hemolysis, hepatic stellate cells activation, and hepatic fatty acid binding protein 4. Additionally, in this review, we also present two clinically relevant models using human precision-cut liver slices and hepatic organoids to examine ALD pathogenesis and progression.

Enhancer-promoter activation by the Kaposi sarcoma-associated herpesvirus episome maintenance protein LANA

Higher-order genome structure influences the transcriptional regulation of cellular genes through the juxtaposition of regulatory elements, such as enhancers, close to promoters of target genes. While enhancer activation has emerged as an important facet of Kaposi sarcoma-associated herpesvirus (KSHV) biology, the mechanisms controlling enhancer-target gene expression remain obscure. Here, we discover that the KSHV genome tethering protein latency-associated nuclear antigen (LANA) potentiates enhancer-target gene expression in primary effusion lymphoma (PEL), a highly aggressive B cell lymphoma causally associated with KSHV. Genome-wide analyses demonstrate increased levels of enhancer RNA transcription as well as activating chromatin marks at LANA-bound enhancers. 3D genome conformation analyses identified genes critical for latency and tumorigenesis as targets of LANA-occupied enhancers, and LANA depletion results in their downregulation. These findings reveal a mechanism in enhancer-gene coordination and describe a role through which the main KSHV tethering protein regulates essential gene expression in PEL.

Assessing NGS-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators with query gene sets. Identification of transcriptional regulators is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Key points

An introduction to available computational methods for predicting functional TRs from a query gene set.A detailed walk-through along with practical concerns and limitations.A systematic benchmark of NGS-based methods in terms of accuracy, sensitivity, coverage, and usability, using 570 TR perturbation-derived gene sets.NGS-based methods outperform motif-based methods. Among NGS methods, those utilizing larger databases and adopting region-centric approaches demonstrate favorable performance. BART, ChIP-Atlas, and Lisa are recommended as these methods have overall better performance in evaluated scenarios.

Pathogenic variants in CRX have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

Dozens of variants in the gene for the homeodomain transcription factor (TF) cone-rod homeobox (CRX) are linked with human blinding diseases that vary in their severity and age of onset. How different variants in this single TF alter its function in ways that lead to a range of phenotypes is unclear. We characterized the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in mouse retina explants carrying knock-ins of two variants, one in the DNA-binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation in these mutant Crx retinas corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, and p.E168d2 has distinct effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are derepressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci are partially predictive of episomal MPRA activity, and distal elements whose accessibility increases later in retinal development are enriched for CREs with silencer activity. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers while having a qualitatively different impact on silencers.

The transcription factor NF-κB orchestrates nucleosome remodeling during the primary response to Toll-like receptor 4 signaling

Inducible nucleosome remodeling at hundreds of latent enhancers and several promoters shapes the transcriptional response to Toll-like receptor 4 (TLR4) signaling in macrophages. We aimed to define the identities of the transcription factors that promote TLR-induced remodeling. An analysis strategy based on ATAC-seq and single-cell ATAC-seq that enriched for genomic regions most likely to undergo remodeling revealed that the transcription factor nuclear factor κB (NF-κB) bound to all high-confidence peaks marking remodeling during the primary response to the TLR4 ligand, lipid A. Deletion of NF-κB subunits RelA and c-Rel resulted in the loss of remodeling at high-confidence ATAC-seq peaks, and CRISPR-Cas9 mutagenesis of NF-κB-binding motifs impaired remodeling. Remodeling selectivity at defined regions was conferred by collaboration with other inducible factors, including IRF3- and MAP-kinase-induced factors. Thus, NF-κB is unique among TLR4-activated transcription factors in its broad contribution to inducible nucleosome remodeling, alongside its ability to activate poised enhancers and promoters assembled into open chromatin.

Decoding gene regulatory circuitry underlying TNBC chemoresistance reveals biomarkers for therapy response and therapeutic targets

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype characterised by extensive intratumoral heterogeneity, high rates of metastasis and chemoresistance, leading to poor clinical outcomes. Despite progress, the mechanistic basis of chemotherapy resistance in TNBC patients remains poorly understood. Here, leveraging single-cell transcriptome datasets of matched longitudinal TNBC chemoresponsive and chemoresistant patient cohorts, we unravel distinct cell subpopulations intricately associated with chemoresistance and the signature genes defining these populations. Notably, using genome-wide mapping of the H3K27ac mark, we show that the expression of these chemoresistance genes is driven via a set of TNBC super-enhancers and associated transcription factor networks across TNBC subtypes. Furthermore, genetic screens reveal that a subset of these transcription factors is essential for the survival of TNBC cells, and their loss increases sensitivity to chemotherapeutic agents. Overall, our study has revealed epigenetic and transcription factor networks underlying chemoresistance and suggests novel avenues to stratify and improve the treatment of patients with a high risk of developing resistance.

Chromatin profiling reveals TFAP4 as a critical transcriptional regulator of bovine satellite cell differentiation

BACKGROUND

Satellite cells are myogenic precursor cells in adult skeletal muscle and play a crucial role in skeletal muscle regeneration, maintenance, and growth. Like embryonic myoblasts, satellite cells have the ability to proliferate, differentiate, and fuse to form multinucleated myofibers. In this study, we aimed to identify additional transcription factors that control gene expression during bovine satellite cell proliferation and differentiation.

RESULTS

Using chromatin immunoprecipitation followed by sequencing, we identified 56,973 and 54,470 genomic regions marked with both the histone modifications H3K4me1 and H3K27ac, which were considered active enhancers, and 50,956 and 59,174 genomic regions marked with H3K27me3, which were considered repressed enhancers, in proliferating and differentiating bovine satellite cells, respectively. In addition, we identified 1,216 and 1,171 super-enhancers in proliferating and differentiating bovine satellite cells, respectively. Analyzing these enhancers showed that in proliferating bovine satellite cells, active enhancers were associated with genes stimulating cell proliferation or inhibiting myoblast differentiation whereas repressed enhancers were associated with genes essential for myoblast differentiation, and that in differentiating satellite cells, active enhancers were associated with genes essential for myoblast differentiation or muscle contraction whereas repressed enhancers were associated with genes stimulating cell proliferation or inhibiting myoblast differentiation. Active enhancers in proliferating bovine satellite cells were enriched with binding sites for many transcription factors such as MYF5 and the AP-1 family transcription factors; active enhancers in differentiating bovine satellite cells were enriched with binding sites for many transcription factors such as MYOG and TFAP4; and repressed enhancers in both proliferating and differentiating bovine satellite cells were enriched with binding sites for NF-kB, ZEB-1, and several other transcription factors. The role of TFAP4 in satellite cell or myoblast differentiation was previously unknown, and through gene knockdown and overexpression, we experimentally validated a critical role for TFAP4 in the differentiation and fusion of bovine satellite cells into myofibers.

CONCLUSIONS

Satellite cell proliferation and differentiation are controlled by many transcription factors such as AP-1, TFAP4, NF-kB, and ZEB-1 whose roles in these processes were previously unknown in addition to those transcription factors such as MYF5 and MYOG whose roles in these processes are widely known.

A SWI/SNF-dependent transcriptional regulation mediated by POU2AF2/C11orf53 at enhancer

Recent studies have identified a previously uncharacterized protein C11orf53 (now named POU2AF2/OCA-T1), which functions as a robust co-activator of POU2F3, the master transcription factor which is critical for both normal and neoplastic tuft cell identity and viability. Here, we demonstrate that POU2AF2 dictates opposing transcriptional regulation at distal enhance elements. Loss of POU2AF2 leads to an inhibition of active enhancer nearby genes, such as tuft cell identity genes, and a derepression of Polycomb-dependent poised enhancer nearby genes, which are critical for cell viability and differentiation. Mechanistically, depletion of POU2AF2 results in a global redistribution of the chromatin occupancy of the SWI/SNF complex, leading to a significant 3D genome structure change and a subsequent transcriptional reprogramming. Our genome-wide CRISPR screen further demonstrates that POU2AF2 depletion or SWI/SNF inhibition leads to a PTEN-dependent cell growth defect, highlighting a potential role of POU2AF2-SWI/SNF axis in small cell lung cancer (SCLC) pathogenesis. Additionally, pharmacological inhibition of SWI/SNF phenocopies POU2AF2 depletion in terms of gene expression alteration and cell viability decrease in SCLC-P subtype cells. Therefore, impeding POU2AF2-mediated transcriptional regulation represents a potential therapeutic approach for human SCLC therapy.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Targeting RORγ inhibits the growth and metastasis of hepatocellular carcinoma

Approximately 80%-90% of hepatocellular carcinomas (HCC) occur in a premalignant environment of fibrosis and abnormal extracellular matrix (ECM), highlighting an essential role of ECM in the tumorigenesis and progress of HCC. However, the determinants of ECM in HCC are poorly defined. Here, we show that nuclear receptor RORγ is highly expressed and amplified in HCC tumors. RORγ functions as an essential activator of the matrisome program via directly driving the expression of major ECM genes in HCC cells. Elevated RORγ increases fibronectin-1 deposition, cell-matrix adhesion, and collagen production, creating a favorable microenvironment to boost liver cancer metastasis. Moreover, RORγ antagonists effectively inhibit tumor growth and metastasis in multiple HCC xenografts and immune-intact models, and they effectively sensitize HCC tumors to sorafenib therapy in mice. Notably, elevated RORγ expression is associated with ECM remodeling and metastasis in patients with HCC. Taken together, we identify RORγ as a key player of ECM remodeling in HCC and as an attractive therapeutic target for advanced HCC.

Stromal Pbrm1 mediates chromatin remodeling necessary for embryo implantation in the mouse uterus

Early gestational loss occurs in approximately 20% of all clinically recognized human pregnancies and is an important cause of morbidity. Either embryonic or maternal defects can cause loss, but a functioning and receptive uterine endometrium is crucial for embryo implantation. We report that the switch/sucrose nonfermentable (SWI/SNF) remodeling complex containing polybromo-1 (PBRM1) and Brahma-related gene 1 (BRG1) is essential for implantation of the embryonic blastocyst on the wall of the uterus in mice. Although preimplantation development is unaffected, conditional ablation of Pbrm1 in uterine stromal cells disrupts progesterone pathways and uterine receptivity. Heart and neural crest derivatives expressed 2 (Hand2) encodes a basic helix-loop-helix (bHLH) transcription factor required for embryo implantation. We identify an enhancer of the Hand2 gene in stromal cells that requires PBRM1 for epigenetic histone modifications/coactivator recruitment and looping with the promoter. In Pbrm1cKO mice, perturbation of chromatin assembly at the promoter and enhancer sites compromises Hand2 transcription, adversely affects fibroblast growth factor signaling pathways, prevents normal stromal-epithelial crosstalk, and disrupts embryo implantation. The mutant female mice are infertile and provide insight into potential causes of early pregnancy loss in humans.

The pleiotropic enhancer enh9 promotes cell proliferation and migration in non-small cell lung cancer via ERMP1 and PD-L1

Enhancers, cis-acting DNA elements for transcriptional regulation, are important regulators of cell identity and disease. However, of the hundreds of thousands of enhancers annotated in the human genome, only a few have been studied for their regulatory mechanisms and functions in cancer progression and therapeutic resistance. Here, we report the pleiotropy of one enhancer (named enh9) in both cell proliferation and migration in non-small cell lung cancer (NSCLC) cells. By integrating multi-genomic data, ERMP1 and PD-L1 were screened out as potential targets of enh9. CUT&Tag sequencing demonstrated that enh9 was involved in the genomic interactions between the transcription factor RELA and the promoters of ERMP1 and PD-L1. In addition, ERMP1 and PD-L1 were validated to be involved in cell proliferation and migration, respectively. Our study fully elucidated the function and transcriptional regulation mechanisms of enh9 in NSCLC. The exploration on enhancers is promising to provide new insights for cancer diagnosis and therapy.

Promoter-centred chromatin interactions associated with EVI1 expression in EVI1+3q- myeloid leukaemia cells

EVI1 expression is associated with poor prognosis in myeloid leukaemia, which can result from Chr.3q alterations that juxtapose enhancers to induce EVI1 expression via long-range chromatin interactions. More often, however, EVI1 expression occurs unrelated to 3q alterations, and it remained unclear if, in these cases, EVI1 expression is similarly caused by aberrant enhancer activation. Here, we report that, in EVI1+3q- myeloid leukaemia cells, the EVI1 promoter interacts via long-range chromatin interactions with promoters of distally located, active genes, rather than with enhancer elements. Unlike in 3q+ cells, EVI1 expression and long-range interactions appear to not depend on CTCF/cohesin, though EVI1+3q- cells utilise an EVI1 promoter-proximal site to enhance its expression that is also involved in CTCF-mediated looping in 3q+ cells. Long-range interactions in 3q- cells connect EVI1 to promoters of multiple genes, whose transcription correlates with EVI1 in EVI1+3q- cell lines, suggesting a shared mechanism of transcriptional regulation. In line with this, CRISPR interference-induced silencing of two of these sites minimally, but consistently reduced EVI1 expression. Together, we provide novel evidence of features associated with EVI1 expression in 3q- leukaemia and consolidate the view that EVI1 in 3q- leukaemia is largely promoter-driven, potentially involving long-distance promoter clustering.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Epigenetic Approaches to Identifying Asthma Endotypes

The prevalence of asthma escalated rapidly in the late 20th century. In 2019, the World Health Organization estimated the global number of people affected by the condition to be approximately 260 million, causing 450,000 deaths during that year. While there have been advances in therapeutics with the emergence of biologics targeting T2-high asthma, there is still little clarity on the mechanisms underlying the origins of both the condition and all of its endotypes. Several biomarkers for particular asthma phenotypes have been documented. These are generally identified from transcriptomics and proteomics protocols and tend to be biased to T2-high phenotypes. In this review, we summarize some suggestions that analysis of epigenomes may provide alternative datasets that inform of broader asthma endotypes and might highlight pathways amenable for therapeutic intervention.

Three-Dimensional Chromosomal Landscape Revealing miR-146a Dysfunctional Enhancer in Lupus and Establishing a CRISPR-Mediated Approach to Inhibit the Interferon Pathway

OBJECTIVE

The diminished expression of microRNA-146a (miR-146a) in systemic lupus erythematosus (SLE) contributes to the aberrant activation of the interferon pathway. Despite its significance, the underlying mechanism driving this reduced expression remains elusive. Considering the integral role of enhancers in steering gene expression, our study seeks to pinpoint the SLE-affected enhancers responsible for modulating miR-146a expression. Additionally, we aim to elucidate the mechanisms by which these enhancers influence the contribution of miR-146a to the activation of the interferon pathway.

METHODS

Circular chromosome conformation capture sequencing and epigenomic profiles were applied to identify candidate enhancers of miR-146a. CRISPR activation was performed to screen functional enhancers. Differential analysis of chromatin accessibility was used to identify SLE-dysregulated enhancers, and the mechanism underlying enhancer dysfunction was investigated by analyzing transcription factor binding. The therapeutic value of a lupus-related enhancer was further evaluated by targeting it in the peripheral blood mononuclear cells (PBMCs) of patients with SLE through a CRISPR activation approach.

RESULTS

We identified shared and cell-specific enhancers of miR-146a in distinct immune cells. An enhancer 32.5 kb downstream of miR-146a possesses less accessibility in SLE, and its chromatin openness was negatively correlated with SLE disease activity. Moreover, CCAAT/enhancer binding protein α, a down-regulated transcription factor in patients with SLE, binds to the 32.5-kb enhancer and induces the epigenomic change of this locus. Furthermore, CRISPR-based activation of this enhancer in SLE PBMCs could inhibit the activity of interferon pathway.

CONCLUSION

Our work defines a promising target for SLE intervention. We adopted integrative approaches to define cell-specific and functional enhancers of the SLE critical gene and investigated the mechanism underlying its dysregulation mediated by a lupus-related enhancer.

LINC00355 promotes gastric carcinogenesis by scaffolding p300 to activate CDC42 transcription and enhancing HNRNPA2B1 to stabilize CDC42 mRNA dependent on m6A

LINC00355 is involved in the tumorigenesis of several types of cancer. We verified that LINC00355 is upregulated in gastric cancer (GC) and contributes to GC cells' proliferation and metastasis. RNA sequencing (RNA-seq) and rescue assays suggested that LINC00355 controls gastric carcinogenesis by regulating the expression of cell division cycle 42 (CDC42) guanosine triphosphatase (GTPases), thereby activating their downstream pathways. Most previous studies have shown that LINC00355 acts as a ceRNA by sponging miRNAs to modulate downstream gene expression. Our group focus on epigenetic regulatory potential of LINC00355 in gene expression. Mechanistically, LINC00355 binds to p300 histone acetyltransferase, specifying the histone modification pattern on the CDC42 promoter to activate CDC42 transcription, thereby altering GC cell biology. In addition, HNRNPA2B1, which is upregulated by LINC00355, recognizes the N6-methyladenosine (m6A) sites of CDC42 and enhances the stability of CDC42 mRNA transcripts. Therefore, LINC00355 is mechanistically, functionally, and clinically oncogenic in GC cells.

An ancient polymorphic regulatory region within the BDNF gene associated with obesity modulates anxiety-like behaviour in mice and humans

Obesity and anxiety are morbidities notable for their increased impact on society during the recent COVID-19 pandemic. Understanding the mechanisms governing susceptibility to these conditions will increase our quality of life and resilience to future pandemics. In the current study, we explored the function of a highly conserved regulatory region (BE5.1) within the BDNF gene that harbours a polymorphism strongly associated with obesity (rs10767664; p = 4.69 × 10-26). Analysis in primary cells suggested that the major T-allele of BE5.1 was an enhancer, whereas the obesity-associated A-allele was not. However, CRISPR/CAS9 deletion of BE5.1 from the mouse genome (BE5.1KO) produced no significant effect on the expression of BDNF transcripts in the hypothalamus, no change in weight gain after 28 days and only a marginally significant increase in food intake. Nevertheless, transcripts were significantly increased in the amygdala of female mice and elevated zero maze and marble-burying tests demonstrated a significant increase in anxiety-like behaviour that could be reversed by diazepam. Consistent with these observations, human GWAS cohort analysis demonstrated a significant association between rs10767664 and anxiousness in human populations. Intriguingly, interrogation of the human GTEx eQTL database demonstrated no effect on BDNF mRNA levels associated with rs10767664 but a highly significant effect on BDNF-antisense (BDNF-AS) gene expression and splicing. The subsequent observation that deletion of BE5.1 also significantly reduced BDNF-AS expression in mice suggests a novel mechanism in the regulation of BDNF expression common to mice and humans, which contributes to the modulation of mood and anxiety in both species.

Epigenetic Effects of Resveratrol on Oncogenic Signaling in Breast Cancer

The crosstalk between oncogenic signaling pathways plays a crucial role in driving cancer development. We previously demonstrated that dietary polyphenols, specifically resveratrol (RSV) and other stilbenoids, epigenetically target oncogenes for silencing via DNA hypermethylation in breast cancer. In the present study, we identify signal transduction regulators among RSV-hypermethylated targets and investigate the functional role of RSV-mediated DNA hypermethylation in the regulation of Hedgehog and Wnt signaling. Non-invasive ER-positive MCF-7 and highly invasive triple-negative MCF10CA1a human breast cancer cell lines were used as experimental models. Upon 9-day exposure to 15 µM RSV, pyrosequencing and qRT-PCR were performed to assess DNA methylation and expression of GLI2 and WNT4, which are upstream regulators of the Hedgehog and Wnt pathways, respectively. Our results showed that RSV led to a DNA methylation increase within GLI2 and WNT4 enhancers, which was accompanied by decreases in gene expression. Consistently, we observed the downregulation of genes downstream of the Hedgehog and Wnt signaling, including common targets shared by both pathways, CCND1 and CYR61. Further analysis using chromatin immunoprecipitation identified increased H3K27 trimethylation and decreased H3K9 and H3K27 acetylation, along with abolishing OCT1 transcription factor binding. Those changes indicate a transcriptionally silent chromatin state at GLI2 and WNT4 enhancers. The inhibition of the Wnt signal transduction was confirmed using a phospho-antibody array that demonstrated suppression of positive and stimulation of negative Wnt regulators. In conclusion, our results provide scientific evidence for dietary polyphenols as epigenetics-modulating agents that act to re-methylate and silence oncogenes, reducing the oncogenic signal transduction. Targeting such an action could be an effective strategy in breast cancer prevention and/or adjuvant therapy.

Impaired bone morphogenetic protein (BMP) signaling pathways disrupt decidualization in endometriosis

Endometriosis is linked to increased infertility and pregnancy complications due to defective endometrial decidualization. We hypothesized that identification of altered signaling pathways during decidualization could identify the underlying cause of infertility and pregnancy complications. Our study reveals that transforming growth factor β (TGFβ) pathways are impaired in the endometrium of individuals with endometriosis, leading to defective decidualization. Through detailed transcriptomic analyses, we discovered abnormalities in TGFβ signaling pathways and key regulators, such as SMAD4, in the endometrium of affected individuals. We also observed compromised activity of bone morphogenetic proteins (BMP), a subset of the TGFβ family, that control endometrial receptivity. Using 3-dimensional models of endometrial stromal and epithelial assembloids, we showed that exogenous BMP2 improved decidual marker expression in individuals with endometriosis. Our findings reveal dysfunction of BMP/SMAD signaling in the endometrium of individuals with endometriosis, explaining decidualization defects and subsequent pregnancy complications in these individuals.

Dynamic and distinct histone modifications facilitate human trophoblast lineage differentiation

The placenta serves as an essential organ for fetal growth throughout pregnancy. Histone modification is a crucial regulatory mechanism involved in numerous biological processes and development. Nevertheless, there remains a significant gap in our understanding regarding the epigenetic regulations that influence trophoblast lineage differentiation, a fundamental aspect of placental development. Here, through comprehensive mapping of H3K4me3, H3K27me3, H3K9me3, and H3K27ac loci during the differentiation of trophoblast stem cells (TSCs) into syncytiotrophoblasts (STs) and extravillous trophoblasts (EVTs), we reveal dynamic reconfiguration in H3K4me3 and H3K27ac patterns that establish an epigenetic landscape conducive to proper trophoblast lineage differentiation. We observe that broad H3K4me3 domains are associated with trophoblast lineage-specific gene expression. Unlike embryonic stem cells, TSCs lack robust bivalent domains. Notably, the repression of ST- and EVT-active genes in TSCs is primarily attributed to the weak H3K4me3 signal rather than bivalent domains. We also unveil the inactivation of TSC enhancers precedes the activation of ST enhancers during ST formation. Our results provide a comprehensive global map of diverse histone modifications, elucidating the dynamic histone modifications during trophoblast lineage differentiation.

Effects of BRD4 inhibitor JQ1 on the expression profile of super-enhancer related lncRNAs and mRNAs in cervical cancer HeLa cells

Objective

To investigate the effects of bromine domain protein 4 (BRD4) inhibitor JQ1 on the expression profile of super-enhancer-related lncRNAs (SE-lncRNAs) and mRNAs in cervical cancer (CC) HeLa-cells.

Methods

The CCK8 method was implemented to detect the inhibitory effect of JQ1 on HeLa cells and explore the best inhibitory concentration. Whole transcriptome sequencing was performed to detect the changes of lncRNAs and mRNAs expression profiles in cells of the JQ1 treatment group and control group, respectively. The differentially expressed SE-lncRNAs were obtained by matching, while the co-expressed mRNAs were obtained by Pearson correlation analysis.

Results

The inhibitory effect of JQ1 on HeLa cell proliferation increased significantly with increasing concentration and treatment time (P < 0.05). Under the experimental conditions of three concentrations of 0.01, 0.1 and 1 μmol/L of JQ1 on HeLa cells at 24, 48, 72 and 120 h, 1 μmol/L of JQ1 at 72 and 120 h had the same cell viability and the strongest cell proliferation inhibition. In order to understand the inhibitory mechanism of JQ1 on HeLa cells, this study analyzed the expression profile differences from the perspective of SE-lncRNAs and mRNAs. A total of 162 SE-lncRNAs were identified, of which 8 SE-lncRNAs were down-regulated and seven SE-lncRNAs were up-regulated. A total of 418 differentially expressed mRNAs related to SE-lncRNAs were identified, of which 395 mRNAs had positive correlation with 12 SE-lncRNAs and 408 mRNAs had negative correlation with 15 SE-lncRNAs.

Conclusion

JQ1 can significantly inhibit the proliferation of HeLa cells and affect the expression profile of SE-lncRNAs and mRNAs.

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.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

Epigenetic regulation and factors that influence the effect of iPSCs-derived neural stem/progenitor cells (NS/PCs) in the treatment of spinal cord injury

Spinal cord injury (SCI) is a severe neurological disorder that causes neurological impairment and disability. Neural stem/progenitor cells (NS/PCs) derived from induced pluripotent stem cells (iPSCs) represent a promising cell therapy strategy for spinal cord regeneration and repair. However, iPSC-derived NS/PCs face many challenges and issues in SCI therapy; one of the most significant challenges is epigenetic regulation and that factors that influence this mechanism. Epigenetics refers to the regulation of gene expression and function by DNA methylation, histone modification, and chromatin structure without changing the DNA sequence. Previous research has shown that epigenetics plays a crucial role in the generation, differentiation, and transplantation of iPSCs, and can influence the quality, safety, and outcome of transplanted cells. In this study, we review the effects of epigenetic regulation and various influencing factors on the role of iPSC-derived NS/PCs in SCI therapy at multiple levels, including epigenetic reprogramming, regulation, and the adaptation of iPSCs during generation, differentiation, and transplantation, as well as the impact of other therapeutic tools (e.g., drugs, electrical stimulation, and scaffolds) on the epigenetic status of transplanted cells. We summarize our main findings and insights in this field and identify future challenges and directions that need to be addressed and explored.

Kaposi's sarcoma-associated herpesvirus terminal repeat regulates inducible lytic gene promoters

The Kaposi's sarcoma-associated herpesvirus (KSHV) genome consists of an approximately 140-kb unique coding region flanked by 30-40 copies of a 0.8-kb terminal repeat (TR) sequence. A gene enhancer recruits transcription-related enzymes by having arrays of transcription factor binding sites. Here, we show that KSHV TR possesses transcription regulatory function with latency-associated nuclear antigen (LANA). Cleavage under targets and release using nuclease demonstrated that TR fragments were occupied by LANA-interacting histone-modifying enzymes in naturally infected cells. The TR was enriched with histone H3K27 acetylation (H3K27Ac) and H3K4 tri-methylation (H3K4me3) modifications and also expressed nascent RNAs. The sites of H3K27Ac and H3K4me3 modifications were also conserved in the KSHV unique region among naturally infected primary effusion lymphoma cells. KSHV origin of lytic replication (Ori-Lyt) showed similar protein and histone modification occupancies with that of TR. In the Ori-Lyt region, the LANA and LANA-interacting proteins colocalized with an H3K27Ac-modified nucleosome along with paused RNA polymerase II. The KSHV transactivator KSHV replication and transcription activator (K-Rta) recruitment sites franked the LANA-bound nucleosome, and reactivation evicted the LANA-bound nucleosome. Including TR fragments in reporter plasmid enhanced inducible viral gene promoter activities independent of the orientations. In the presence of TR in reporter plasmids, K-Rta transactivation was drastically increased, while LANA acquired the promoter repression function. KSHV TR, therefore, functions as an enhancer for KSHV inducible genes. However, in contrast to cellular enhancers bound by multiple transcription factors, perhaps the KSHV enhancer is predominantly regulated by the LANA nuclear body.IMPORTANCEEnhancers are a crucial regulator of differential gene expression programs. Enhancers are the cis-regulatory sequences determining target genes' spatiotemporal and quantitative expression. Here, we show that Kaposi's sarcoma-associated herpesvirus (KSHV) terminal repeats fulfill the enhancer definition for KSHV inducible gene promoters. The KSHV enhancer is occupied by latency-associated nuclear antigen (LANA) and its interacting proteins, such as CHD4. Neighboring terminal repeat (TR) fragments to lytic gene promoters drastically enhanced KSHV replication and transcription activator and LANA transcription regulatory functions. This study, thus, proposes a new latency-lytic switch model in which TR accessibility to the KSHV gene promoters regulates viral inducible gene expression.

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

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

Chromatin organization of muscle stem cell

The proper functioning of skeletal muscles is essential throughout life. A crucial crosstalk between the environment and several cellular mechanisms allows striated muscles to perform successfully. Notably, the skeletal muscle tissue reacts to an injury producing a completely functioning tissue. The muscle's robust regenerative capacity relies on the fine coordination between muscle stem cells (MuSCs or "satellite cells") and their specific microenvironment that dictates stem cells' activation, differentiation, and self-renewal. Critical for the muscle stem cell pool is a fine regulation of chromatin organization and gene expression. Acquiring a lineage-specific 3D genome architecture constitutes a crucial modulator of muscle stem cell function during development, in the adult stage, in physiological and pathological conditions. The context-dependent relationship between genome structure, such as accessibility and chromatin compartmentalization, and their functional effects will be analysed considering the improved 3D epigenome knowledge, underlining the intimate liaison between environmental encounters and epigenetics.

SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion

Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 limits endogenous DNA damage, thereby suppressing cGAS-STING-dependent signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a PD-L1 transcriptional regulatory element, thereby promoting PD-L1 expression in cancer cells. SMARCAL1 loss hinders the ability of tumor cells to induce PD-L1 in response to genomic instability, enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a promising target for cancer immunotherapy.

The transcriptional regulatory network modulating human trophoblast stem cells to extravillous trophoblast differentiation

During human pregnancy, extravillous trophoblasts play crucial roles in placental invasion into the maternal decidua and spiral artery remodeling. However, regulatory factors and their action mechanisms modulating human extravillous trophoblast specification have been unknown. By analyzing dynamic changes in transcriptome and enhancer profile during human trophoblast stem cell to extravillous trophoblast differentiation, we define stage-specific regulators, including an early-stage transcription factor, TFAP2C, and multiple late-stage transcription factors. Loss-of-function studies confirm the requirement of all transcription factors identified for adequate differentiation, and we reveal that the dynamic changes in the levels of TFAP2C are essential. Notably, TFAP2C pre-occupies the regulatory elements of the inactive extravillous trophoblast-active genes during the early stage of differentiation, and the late-stage transcription factors directly activate extravillous trophoblast-active genes, including themselves as differentiation further progresses, suggesting sequential actions of transcription factors assuring differentiation. Our results reveal stage-specific transcription factors and their inter-connected regulatory mechanisms modulating extravillous trophoblast differentiation, providing a framework for understanding early human placentation and placenta-related complications.

Integrating leiomyoma genetics, epigenomics, and single-cell transcriptomics reveals causal genetic variants, genes, and cell types

Uterine fibroids (UF), that can disrupt normal uterine function and cause significant physical and psychological health problems, are observed in nearly 70% of women of reproductive age. Although heritable genetics is a significant risk factor, specific genetic variations and gene targets causally associated with UF are poorly understood. Here, we performed a meta-analysis on existing fibroid genome-wide association studies (GWAS) and integrated the identified risk loci and potentially causal single nucleotide polymorphisms (SNPs) with epigenomics, transcriptomics, 3D chromatin organization from diverse cell types as well as primary UF patient's samples. This integrative analysis identifies 24 UF-associated risk loci that potentially target 394 genes, of which 168 are differentially expressed in UF tumors. Critically, integrating this data with single-cell gene expression data from UF patients reveales the causal cell types with aberrant expression of these target genes. Lastly, CRISPR-based epigenetic repression (dCas9-KRAB) or activation (dCas9-p300) in a UF disease-relevant cell type further refines and narrows down the potential gene targets. Our findings and the methodological approach indicate the effectiveness of integrating multi-omics data with locus-specific epigenetic editing approaches for identifying gene- and celt type-targets of disease-relevant risk loci.

Decoding the Epigenetics and Chromatin Loop Dynamics of Androgen Receptor-Mediated Transcription

Androgen receptor (AR)-mediated transcription plays a critical role in normal prostate development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigated the AR-driven epigenetic and chromosomal chromatin looping changes. We collected a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we found that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we showed that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then proposed and experimentally validated an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide new insight into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Enhancer-mediated FOXO3 expression promotes MSC adipogenic differentiation by activating autophagy

BACKGROUND

Mesenchymal stem cells (MSCs) are pluripotent stem cells capable of differentiating into osteocytes, adipocytes and chondrocytes. However, in osteoporosis, the balance of differentiation is tipped toward adipogenesis and the key mechanism is controversial. Researches have shown that, as upstream regulatory elements of gene expression, enhancers ar involved in the expression of identity genes. In this study, we identified enhancers-mediated gene FOXO3 promoting MSC adipogenic differentiation by activating autophagy.

METHODS

We integrated data of RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq) and ATAC-sequencing (ATAC-seq) to find the identity gene FOXO3. The expression of FOXO3 protein, adipogenic transcription factors and the substrate of autophagy were measured by western blotting. The Oil Red O (ORO) staining was used to visualize the adipogenesis of MSCs. Immunohistochemistry was used to visualize the FOXO3 expression in adipocytes in bone marrow. Immunofluorescence was used to detect the expression of PPARγ and LC3B.

RESULTS

During adipogenesis, enhancers redistribute to genes associated with adipogenic differentiation, among which we identified the pivotal identity gene FOXO3. FOXO3 could promote the expression of the adipogenic transcription factors PPARγ, CEBPα, and CEBPβ during adipogenic differentiation, while PPARγ, CEBPα, and CEBPβ could in turn bind to FOXO3 and continue to promote FOXO3 expression to form a positive feedback loop. Consistently elevated FOXO3 expression promotes autophagy by activating the PI3K-AKT pathway which mediates adipogenic differentiation.

CONCLUSIONS

Pivotal identity gene FOXO3 promotes autophagy by activating PI3K-AKT pathway, which provokes adipogenic differentiation of MSCs. Enhancer-regulated adipogenic identity gene FOXO3 could be an attractive treatment for osteoporosis.

PARP-1 is a transcriptional rheostat of metabolic and bivalent genes during development

PARP-1 participates in various cellular processes, including gene regulation. In Drosophila, PARP-1 mutants undergo developmental arrest during larval-to-pupal transition. In this study, we investigated PARP-1 binding and its transcriptional regulatory role at this stage. Our findings revealed that PARP-1 binds and represses active metabolic genes, including glycolytic genes, whereas activating low-expression developmental genes, including a subset of "bivalent" genes in third-instar larvae. These bivalent promoters, characterized by dual enrichment of low H3K4me3 and high H3K27me3, a unimodal H3K4me1 enrichment at the transcription start site (conserved in C. elegans and zebrafish), H2Av depletion, and high accessibility, may persist throughout development. In PARP-1 mutant third-instar larvae, metabolic genes typically down-regulated during the larval-to-pupal transition in response to reduced energy needs were repressed by PARP-1. Simultaneously, developmental and bivalent genes typically active at this stage were activated by PARP-1. In addition, glucose and ATP levels were significantly reduced in PARP-1 mutants, suggesting an imbalance in metabolic regulation. We propose that PARP-1 is essential for maintaining the delicate balance between metabolic and developmental gene expression programs to ensure proper developmental progression.

HNF4A guides the MLL4 complex to establish and maintain H3K4me1 at gene regulatory elements

Hepatocyte nuclear factor 4A (HNF4A/NR2a1), a transcriptional regulator of hepatocyte identity, controls genes that are crucial for liver functions, primarily through binding to enhancers. In mammalian cells, active and primed enhancers are marked by monomethylation of histone 3 (H3) at lysine 4 (K4) (H3K4me1) in a cell type-specific manner. How this modification is established and maintained at enhancers in connection with transcription factors (TFs) remains unknown. Using analysis of genome-wide histone modifications, TF binding, chromatin accessibility and gene expression, we show that HNF4A is essential for an active chromatin state. Using HNF4A loss and gain of function experiments in vivo and in cell lines in vitro, we show that HNF4A affects H3K4me1, H3K27ac and chromatin accessibility, highlighting its contribution to the establishment and maintenance of a transcriptionally permissive epigenetic state. Mechanistically, HNF4A interacts with the mixed-lineage leukaemia 4 (MLL4) complex facilitating recruitment to HNF4A-bound regions. Our findings indicate that HNF4A enriches H3K4me1, H3K27ac and establishes chromatin opening at transcriptional regulatory regions.

The ESCRT protein CHMP5 promotes T cell leukemia by controlling BRD4-p300-dependent transcription

Oncogene activity rewires cellular transcription, creating new transcription networks to which cancer cells become addicted, by mechanisms that are still poorly understood. Using human and mouse models of T cell acute lymphoblastic leukemia (T-ALL), we identify an essential nuclear role for CHMP5, a cytoplasmic endosomal sorting complex required for transport (ESCRT) protein, in establishing and maintaining the T-ALL transcriptional program. Nuclear CHMP5 promoted the T-ALL gene program by augmenting recruitment of the co-activator BRD4 by the histone acetyl transferase p300 selectively at enhancers and super-enhancers, an interaction that potentiated H3K27 acetylation at these regulatory enhancers. Consequently, loss of CHMP5 diminished BRD4 occupancy at enhancers and super-enhancers and impaired RNA polymerase II pause release, which resulted in downregulation of key T-ALL genes, notably MYC. Reinforcing its importance in T-ALL pathogenesis, CHMP5 deficiency mitigated chemoresistance in human T-ALL cells and abrogated T-ALL induction by oncogenic NOTCH1 in vivo. Thus, the ESCRT protein CHMP5 is an essential positive regulator of the transcriptional machinery promoting T-ALL disease.

Non-CG DNA methylation and MeCP2 stabilize repeated tuning of long genes that distinguish closely related neuron types

The extraordinary diversity of neuron types in the mammalian brain is delineated at the highest resolution by subtle gene expression differences that may require specialized molecular mechanisms to be maintained. Neurons uniquely express the longest genes in the genome and utilize neuron-enriched non-CG DNA methylation (mCA) together with the Rett syndrome protein, MeCP2, to control gene expression, but the function of these unique gene structures and machinery in regulating finely resolved neuron type-specific gene programs has not been explored. Here, we employ epigenomic and spatial transcriptomic analyses to discover a major role for mCA and MeCP2 in maintaining neuron type-specific gene programs at the finest scale of cellular resolution. We uncover differential susceptibility to MeCP2 loss in neuronal populations depending on global mCA levels and dissect methylation patterns and intragenic enhancer repression that drive overlapping and distinct gene regulation between neuron types. Strikingly, we show that mCA and MeCP2 regulate genes that are repeatedly tuned to differentiate neuron types at the highest cellular resolution, including spatially resolved, vision-dependent gene programs in the visual cortex. These repeatedly tuned genes display genomic characteristics, including long length, numerous intragenic enhancers, and enrichment for mCA, that predispose them to regulation by MeCP2. Thus, long gene regulation by the MeCP2 pathway maintains differential gene expression between closely-related neurons to facilitate the exceptional cellular diversity in the complex mammalian brain.