Machine and Deep Learning Methods for Predicting 3D Genome Organization.

Three-dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, topologically associating domains (TADs), and A/B compartments, play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers and transcription factor binding site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, and TAD boundaries) and analyze their pros and cons. We also point out obstacles to the computational prediction of 3D interactions and suggest future research directions.

Prediction of Enhancer-Gene Interactions Using Chromatin-Conformation Capture and Epigenome Data Using STARE.

Disentangling the relationship of enhancers and genes is an ongoing challenge in epigenomics. We present STARE, our software to quantify the strength of enhancer-gene interactions based on enhancer activity and chromatin contact data. It implements the generalized Activity-by-Contact (gABC) score, which allows predicting putative target genes of candidate enhancers over any desired genomic distance. The only requirement for its application is a measurement of enhancer activity. In addition to regulatory interactions, STARE calculates transcription factor (TF) affinities on gene level. We illustrate its usage on a public single-cell data set of the human heart by predicting regulatory interactions on cell type level, by giving examples on how to integrate them with other data modalities, and by constructing TF affinity matrices.

Learning Enhancer-Gene associations from Bulk Transcriptomic and Epigenetic Sequencing Data with STITCHIT.

To reveal gene regulation mechanisms, it is essential to understand the role of regulatory elements, which are possibly distant from gene promoters. Integrative analysis of epigenetic and transcriptomic data can be used to gain insights into gene-expression regulation in specific phenotypes. Here, we discuss STITCHIT, an approach to dissect epigenetic variation in a gene-specific manner across many samples for the identification of regulatory elements without relying on peak calling algorithms. The obtained genomic regions are then further refined using a regularized linear model approach, which can also be used to predict gene expression. We illustrate the use of STITCHIT using H3k27ac ChIP-seq and RNA-seq data from the International Human Epigenome Consortium (IHEC).

Enhancer landscape of lung neuroendocrine tumors reveals regulatory and developmental signatures with potential theranostic implications.

Well-differentiated low-grade lung neuroendocrine tumors (lung carcinoids or LNETs) are histopathologically classified as typical and atypical LNETs, but each subtype is still heterogeneous at both the molecular level and its clinical manifestation. Here, we report genome-wide profiles of primary LNETs' cis-regulatory elements by H3K27ac ChIP-seq with matching RNA-seq profiles. Analysis of these regulatory landscapes revealed three regulatory subtypes, independent of the typical/atypical classification. We identified unique differentiation signals that delineate each subtype. The "proneuronal" subtype emerges under the influence of ASCL1, SOX4, and TCF4 transcription factors, embodying a pronounced proneuronal signature. The "luminal-like" subtype is characterized by gain of acetylation at markers of luminal cells and GATA2 activation and loss of LRP5 and OTP. The "HNF+" subtype is characterized by a robust enhancer landscape driven by HNF1A, HNF4A, and FOXA3, with notable acetylation and expression of FGF signaling genes, especially FGFR3 and FGFR4, pivotal components of the FGF pathway. Our findings not only deepen the understanding of LNETs' regulatory and developmental diversity but also spotlight the HNF+ subtype's reliance on FGFR signaling. We demonstrate that targeting this pathway with FGF inhibitors curtails tumor growth both in vitro and in xenograft models, unveiling a potential vulnerability and paving the way for targeted therapies. Overall, our work provides an important resource for studying LNETs to reveal regulatory networks, differentiation signals, and therapeutically relevant dependencies.

Developmental and housekeeping transcriptional programs display distinct modes of enhancer-enhancer cooperativity in Drosophila.

Genomic enhancers are key transcriptional regulators which, upon the binding of sequence-specific transcription factors, activate their cognate target promoters. Although enhancers have been extensively studied in isolation, a substantial number of genes have more than one simultaneously active enhancer, and it remains unclear how these cooperate to regulate transcription. Using Drosophila melanogaster S2 cells as a model, we assay the activities of more than a thousand individual enhancers and about a million enhancer pairs toward housekeeping and developmental core promoters with STARR-seq. We report that housekeeping and developmental enhancers show distinct modes of enhancer-enhancer cooperativity: while housekeeping enhancers are additive such that their combined activity mirrors the sum of their individual activities, developmental enhancers are super-additive and combine multiplicatively. Super-additivity between developmental enhancers is promiscuous and neither depends on the enhancers' endogenous genomic contexts nor on specific transcription factor motif signatures. However, it can be further boosted by Twist and Trl motifs and saturates for the highest levels of enhancer activity. These results have important implications for our understanding of gene regulation in complex multi-enhancer developmental loci and genomically clustered housekeeping genes, providing a rationale to interpret the transcriptional impact of non-coding mutations at different loci.

Regulatory genome annotation of 33 insect species.

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival.

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

Methylation patterns at the adjacent CpG sites within enhancers are a part of cell identity.

BACKGROUND: It is generally accepted that methylation status of CpG sites spaced up to 50 bp apart is correlated, and accumulation of locally disordered methylation at adjacent CpG sites is involved in neoplastic transformation, acting in similar way as stochastic accumulation of mutations. RESULTS: We used EPIC microarray data from 596 samples, representing 12 healthy tissue and cell types, as well as 572 blood cancer specimens to analyze methylation status of adjacent CpG sites across human genome, and subsequently validated our findings with NGS and Sanger sequencing. Our analysis showed that there is a subset of the adjacent CpG sites in human genome, with cytosine at one CpG site methylated and the other devoid of methyl group. These loci map to enhancers that are targeted by families of transcription factors involved in cell differentiation. Moreover, our results suggest that the methylation at these loci differ between alleles within a cell, what allows for remarkable level of heterogeneity of methylation patterns. However, different types of specialized cells acquire only one specific and stable pattern of methylation at each of these loci and that pattern is to a large extent lost during neoplastic transformation. CONCLUSIONS: We identified a substantial number of adjacent CpG loci in human genome that display remarkably stable and cell type specific methylation pattern. The methylation pattern at these loci appears to reflect different methylation of alleles in cells. Furthermore, we showed that changes of methylation status at those loci are likely to be involved in regulation of the activity of enhancers and contribute to neoplastic transformation.

Super-enhancer-driven ZFP36L1 promotes PD-L1 expression in infiltrative gastric cancer.

Gastric cancer (GC) is a major cause of cancer-related mortality worldwide. Despite the widespread recognition of tumor immunotherapy in treating unresectable GC, challenges, including ineffective immunotherapy and drug resistance, persist. Therefore, understanding the regulatory mechanisms of PD-L1, particularly in the context of super-enhancers (SEs) and zinc finger protein 36 ring finger protein-like 1 (ZFP36L1) RNA-binding protein, is crucial. In this study, we performed H3K27ac Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing, investigated the heterogeneity of SEs between two GC subtypes with differential growth patterns, and revealed the immune escape signatures driven by ZFP36L1-SE in infiltrative GC through SEs inhibitors treatment. The regulation of ZFP36L1 to PD-L1 was evaluated by quantitative PCR, western blot, flow cytometry, and immunohistochemistry. Furthermore, we explored its regulatory mechanisms using a combination of molecular biology techniques, including luciferase reporter assay, GST/RNA pull-down, chromatin immunoprecipitation (ChIP)/RIP experiments, and in vivo functional assays. We demonstrated that ZFP36L1, driven by an SE, enhances IFN-γ-induced PD-L1 expression, with SPI1 identified as the specific transcription factor binding to ZFP36L1-SE. Mechanistically, ZFP36L1 binds to the adenylate uridylate-rich element in the 3' untranslated region (3'UTR) of HDAC3 mRNA, exacerbating its mRNA decay, and thereby facilitating PD-L1 abnormal transcriptional activation. Collectively, our findings provide mechanistic insights into the role of the SPI1-ZFP36L1-HDAC3-PD-L1 signaling axis in orchestrating immune escape mechanisms in GC, thereby offering valuable insights into the potential targets for immune checkpoint therapy in GC management.

ZIC2 and ZIC3 promote SWI/SNF recruitment to safeguard progression towards human primed pluripotency.

The primed epiblast acts as a transitional stage between the relatively homogeneous naïve epiblast and the gastrulating embryo. Its formation entails coordinated changes in regulatory circuits driven by transcription factors and epigenetic modifications. Using a multi-omic approach in human embryonic stem cell models across the spectrum of peri-implantation development, we demonstrate that the transcription factors ZIC2 and ZIC3 have overlapping but essential roles in opening primed-specific enhancers. Together, they are essential to facilitate progression to and maintain primed pluripotency. ZIC2/3 accomplish this by recruiting SWI/SNF to chromatin and loss of ZIC2/3 or degradation of SWI/SNF both prevent enhancer activation. Loss of ZIC2/3 also results in transcriptome changes consistent with perturbed Polycomb activity and a shift towards the expression of genes linked to differentiation towards the mesendoderm. Additionally, we find an intriguing dependency on the transcriptional machinery for sustained recruitment of ZIC2/3 over a subset of primed-hESC specific enhancers. Taken together, ZIC2 and ZIC3 regulate highly dynamic lineage-specific enhancers and collectively act as key regulators of human primed pluripotency.

Massively parallel characterization of insulator activity across the genome.

A key question in regulatory genomics is whether cis-regulatory elements (CREs) are modular elements that can function anywhere in the genome, or whether they are adapted to certain genomic locations. To distinguish between these possibilities we develop MPIRE (Massively Parallel Integrated Regulatory Elements), a technology for recurrently assaying CREs at thousands of defined locations across the genome in parallel. MPIRE allows us to separate the intrinsic activity of CREs from the effects of their genomic environments. We apply MPIRE to assay three insulator sequences at thousands of genomic locations and find that each insulator functions in locations with distinguishable properties. All three insulators can block enhancers, but each insulator blocks specific enhancers at specific locations. However, only ALOXE3 appears to block heterochromatin silencing. We conclude that insulator function is highly context dependent and that MPIRE is a robust method for revealing the context dependencies of CREs.

Dynamic chromatin architecture identifies new autoimmune-associated enhancers for IL2 and novel genes regulating CD4+ T cell activation.

Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter-cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.

Network Analysis of Enhancer-Promoter Interactions Highlights Cell-Type-Specific Mechanisms of Transcriptional Regulation Variation.

Gene expression is orchestrated by a complex array of gene regulatory elements that govern transcription in a cell-type-specific manner. Though previously studied, the ability to utilize regulatory elements to identify disrupting variants remains largely elusive. To identify important factors within these regions, we generated enhancer-promoter interaction (EPI) networks and investigated the presence of disease-associated variants that fall within these regions. Our study analyzed six neuronal cell types across neural differentiation, allowing us to examine closely related cell types and across differentiation stages. Our results expand upon previous findings of cell-type specificity of enhancer, promoter, and transcription factor binding sites. Notably, we find that regulatory regions within EPI networks can identify the enrichment of variants associated with neuropsychiatric disorders within specific cell types and network sub-structures. This enrichment within sub-structures can allow for a better understanding of potential mechanisms by which variants may disrupt transcription. Together, our findings suggest that EPIs can be leveraged to better understand cell-type-specific regulatory architecture and used as a selection method for disease-associated variants to be tested in future functional assays. Combined with these future functional characterization assays, EPIs can be used to better identify and characterize regulatory variants' effects on such networks and model their mechanisms of gene regulation disruption across different disorders. Such findings can be applied in practical settings, such as diagnostic tools and drug development.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance.

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remain unclear. Here, we systematically identify enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We find that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further show that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

Sexually dimorphic renal expression of mouse Klotho is directed by a kidney-specific distal enhancer responsive to HNF1b.

Transcription enhancers are genomic sequences regulating common and tissue-specific genes and their disruption can contribute to human disease development and progression. Klotho, a sexually dimorphic gene specifically expressed in kidney, is well-linked to kidney dysfunction and its deletion from the mouse genome leads to premature aging and death. However, the sexually dimorphic regulation of Klotho is not understood. Here, we characterize two candidate Klotho enhancers using H3K27ac epigenetic marks and transcription factor binding and investigate their functions, individually and combined, through CRISPR-Cas9 genome engineering. We discovered that only the distal (E1), but not the proximal (E2) candidate region constitutes a functional enhancer, with the double deletion not causing Klotho expression to further decrease. E1 activity is dependent on HNF1b transcription factor binding site within the enhancer. Further, E1 controls the sexual dimorphism of Klotho as evidenced by qPCR and RNA-seq. Despite the sharp reduction of Klotho mRNA, unlike germline Klotho knockouts, mutant mice present normal phenotype, including weight, lifespan, and serum biochemistry. Lastly, only males lacking E1 display more prominent acute, but not chronic kidney injury responses, indicating a remarkable range of potential adaptation to isolated Klotho loss, especially in female E1 knockouts, retaining renoprotection despite over 80% Klotho reduction.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements.

CRISPR-based gene activation (CRISPRa) is a strategy for upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis-regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis-regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans-acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate genes in a cell type-specific manner.

Notch1 Phase Separation Coupled Percolation facilitates target gene expression and enhancer looping.

The Notch receptor is a pleiotropic signaling protein that translates intercellular ligand interactions into changes in gene expression via the nuclear localization of the Notch intracellular Domain (NICD). Using a combination of immunohistochemistry, RNA in situ, Optogenetics and super-resolution live imaging of transcription in human cells, we show that the N1ICD can form condensates that positively facilitate Notch target gene expression. We determined that N1ICD undergoes Phase Separation Coupled Percolation (PSCP) into transcriptional condensates, which recruit, enrich, and encapsulate a broad set of core transcriptional proteins. We show that the capacity for condensation is due to the intrinsically disordered transcriptional activation domain of the N1ICD. In addition, the formation of such transcriptional condensates acts to promote Notch-mediated super enhancer-looping and concomitant activation of the MYC protooncogene expression. Overall, we introduce a novel mechanism of Notch1 activity in which discrete changes in nuclear N1ICD abundance are translated into the assembly of transcriptional condensates that facilitate gene expression by enriching essential transcriptional machineries at target genomic loci.

Anatomy of a superenhancer.

Interferon regulatory factor-8 (IRF8) is the lineage determining transcription factor for the type one classical dendritic cell (cDC1) subset, a terminal selector for plasmacytoid dendritic cells and important for the function of monocytes. Studies of Irf8 gene regulation have identified several enhancers controlling its activity during development of progenitors in the bone marrow that precisely regulate expression at distinct developmental stages. Each enhancer responds to distinct transcription factors that are expressed at each stage. IRF8 is first expressed in early progenitors that form the monocyte dendritic cell progenitor (MDP) in response to induction of the transcription factor CCAAT/enhancer-binding protein alpha (C/EBPα) acting at the Irf8 +56 kb enhancer. IRF8 levels increase further as the MDP transits into the common dendritic cell progenitor (CDP) in response to E protein activity at the Irf8 +41 kb enhancer. Upon Nfil3-induction in CDPs leading to specification of the cDC1 progenitor, abrupt induction of BATF3 forms the JUN/BATF3/IRF8 heterotrimer that activates the Irf8 +32 kb enhancer that sustains Irf8 autoactivation throughout the cDC1 lifetime. Deletions of each of these enhancers has revealed their stage dependent activation. Surprisingly, studies of compound heterozygotes for each combination of enhancer deletions revealed that activation of each subsequent enhancer requires the successful activation of the previous enhancer in strictly cis-dependent mechanism. Successful progression of enhancer activation is finely tuned to alter the functional accessibility of subsequent enhancers to factors active in the next stage of development. The molecular basis for these phenomenon is still obscure but could have implications for genomic regulation in a broader developmental context.

LXR-dependent enhancer activation regulates the temporal organization of the liver's response to refeeding leading to lipogenic gene overshoot.

Transitions between the fed and fasted state are common in mammals. The liver orchestrates adaptive responses to feeding/fasting by transcriptionally regulating metabolic pathways of energy usage and storage. Transcriptional and enhancer dynamics following cessation of fasting (refeeding) have not been explored. We examined the transcriptional and chromatin events occurring upon refeeding in mice, including kinetic behavior and molecular drivers. We found that the refeeding response is temporally organized with the early response focused on ramping up protein translation while the later stages of refeeding drive a bifurcated lipid synthesis program. While both the cholesterol biosynthesis and lipogenesis pathways were inhibited during fasting, most cholesterol biosynthesis genes returned to their basal levels upon refeeding while most lipogenesis genes markedly overshoot above pre-fasting levels. Gene knockout, enhancer dynamics, and ChIP-seq analyses revealed that lipogenic gene overshoot is dictated by LXRα. These findings from unbiased analyses unravel the mechanism behind the long-known phenomenon of refeeding fat overshoot.

Progress on SOX9 and its enhancers in mammalian sex determination.

The sex determination in mammals refers to the development of an initial bipotential organ, termed the bipotential gonad/genital ridge, into either a testis or an ovary at the early stages of embryonic development, under the precise regulation of transcription factors. SOX9 (SRY-box transcription factor 9) is a multifunctional transcription factor in mammalian development and plays a critical role in sex determination and subsequent male reproductive organs development. Recent studies have shown that several enhancers upstream of SOX9 also play an important role in the process of sex determination. In this review, we summarize the progress on the role of SOX9 and its gonadal enhancers in sex determination. This review will facilitate to understand the regulatory mechanism of sex determination of SOX9 and provides a theoretical basis for the further development of animal sex manipulation technologies.