Enhancer Regulation of Transcriptional Bursting Parameters Revealed by Forced Chromatin Looping

Regulation of RNA polymerase II transcription through re-initiation and bursting

The regulation of RNA polymerase II (RNAPII) activity requires orchestrated responses among genomic regulatory sequences and an expansive set of proteins and protein complexes. Despite intense study over five decades, mechanistic insights continue to emerge. Within the past 10 years, live-cell imaging and single-cell transcriptomics experiments have yielded new information about enhancer-promoter communication, transcription factor dynamics, and the kinetics of RNAPII transcription activation. These insights have established RNAPII re-initiation and bursting as a common regulatory phenomenon with widespread implications for gene regulation in health and disease. Here, we summarize regulatory strategies that help control RNAPII bursting in eukaryotic cells, which is defined as short periods of active transcription followed by longer periods of inactivity. We focus on RNAPII re-initiation (i.e., a "burst" of two or more polymerases that initiate from the same promoter), with an emphasis on molecular mechanisms, open questions, and controversies surrounding this distinct regulatory stage.

Exploring the central dogma through the lens of gene expression noise

Over the past two decades, cell-to-cell heterogeneity has garnered increasing attention due to its critical role in both developmental and pathological processes. This growing interest has been driven, in part, by the advancements in live-cell and single-molecule imaging techniques. These techniques have provided mechanistic insights into how processes, transcription in particular, contribute to gene expression noise and, ultimately, cell-to-cell heterogeneity. More recently, however, research has expanded to explore how downstream steps in the central dogma influence gene expression noise. In this review, we mostly examine the impact of transcriptional processes on the generation of gene expression noise but also discuss how post-transcriptional mechanisms modulate noise and its propagation to the protein level. This evaluation emphasizes the need for further investigation into how processes beyond transcription shape gene expression noise, highlighting unanswered questions that remain in the field.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single-nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, and with Polycomb occupancy. We show that histone genes are expressed with allelic imbalance in mESCs, and are involved in haplotype-specific chromatin contacts marked by H3K27me3. Conditional knockouts of Polycomb enzymatic subunits, Ezh2 or Ring1, show that one-third of ASE genes, including histone genes, is regulated through Polycomb repression. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Predicting gene expression changes from chromatin structure modification

Spatial organization of chromatin plays a critical role in gene transcription, but connecting population-averaged HiC data to functional outcomes remains a challenge. We present a computational framework linking HiC contact map to gene transcription. Utilizing a bead-spring polymer model informed by HiC contact maps, we generate an ensemble of 3D conformations for a given genomic locus. These conformations are then coupled to gene transcription levels through a Markov chain model, with transition rates derived from molecular dynamics simulations. The efficacy of this framework is demonstrated by simulating the perturbation of a CTCF-mediated TAD boundary, impacting the expression of sox9 and kcnj2. Our model quantitatively reproduces experimentally observed changes in gene expression, revealing that the increased kcnj2 transcription is a consequence of enhancers within the sox9 TAD becoming accessible upon boundary disruption. Quantifying enhancer impact, our model can also identify functional enhancers. This framework enhances our understanding of the relationship between chromosome spatial architecture and gene regulation.

The maternal-to-zygotic transition: reprogramming of the cytoplasm and nucleus

A fertilized egg is initially transcriptionally silent and relies on maternally provided factors to initiate development. For embryonic development to proceed, the oocyte-inherited cytoplasm and the nuclear chromatin need to be reprogrammed to create a permissive environment for zygotic genome activation (ZGA). During this maternal-to-zygotic transition (MZT), which is conserved in metazoans, transient totipotency is induced and zygotic transcription is initiated to form the blueprint for future development. Recent technological advances have enhanced our understanding of MZT regulation, revealing common themes across species and leading to new fundamental insights about transcription, mRNA decay and translation.

Identification of molecular determinants of gene-specific bursting patterns by high-throughput imaging screens

Stochastic transcriptional bursting is a universal property of active genes. While different genes exhibit distinct bursting patterns, the molecular mechanisms that govern gene-specific stochastic bursting are largely unknown. We have developed a high-throughput-imaging-based screening strategy to identify cellular factors that determine the bursting patterns of native genes in human cells. We identify protein acetylation as a prominent effector of burst frequency and burst size acting via decreasing off-times and gene-specific changes in the on-time. These effects are not correlated with promoter acetylation. Instead, we demonstrate acetylation of the Integrator complex as a key determinant of gene bursting that alters Integrator interactions with transcription elongation and RNA processing factors but without affecting pausing. Our results suggest a prominent role for non-histone acetylation of a transcription cofactors as a mechanism for modulation of bursting via a far-downstream checkpoint.

The spatial choreography of mRNA biosynthesis

Properly timed gene expression is essential for all aspects of organismal physiology. Despite significant progress, our understanding of the complex mechanisms governing the dynamics of gene regulation in response to internal and external signals remains incomplete. Over the past decade, advances in technologies like light and cryo-electron microscopy (Cryo-EM), cryo-electron tomography (Cryo-ET) and high-throughput sequencing have spurred new insights into traditional paradigms of gene expression. In this Review, we delve into recent concepts addressing 'where' and 'when' gene transcription and RNA splicing occur within cells, emphasizing the dynamic spatial and temporal organization of the cell nucleus.

Dual α-globin-truncated erythropoietin receptor knockin restores hemoglobin production in α-thalassemia-derived erythroid cells

The most severe form of α-thalassemia results from loss of all four copies of α-globin. Postnatally, patients face challenges similar to β-thalassemia, including severe anemia and erythrotoxicity due to the imbalance of β-globin and α-globin chains. Despite progress in genome editing treatments for β-thalassemia, there is no analogous curative option for α-thalassemia. To address this, we designed a Cas9/AAV6-mediated genome editing strategy that integrates a functional α-globin gene into the β-globin locus in α-thalassemia patient-derived hematopoietic stem and progenitor cells (HSPCs). Incorporation of a truncated erythropoietin receptor transgene into the α-globin integration cassette significantly increased erythropoietic output from edited HSPCs and led to the most robust production of α-globin, and consequently hemoglobin tetramers. By directing edited HSPCs toward increased production of clinically relevant erythroid cells, this approach has the potential to mitigate the limitations of current treatments for the hemoglobinopathies, including low genome editing and low engraftment rates.

3D chromatin hubs as regulatory units of identity and survival in human acute leukemia

Cancer progression involves genetic and epigenetic changes that disrupt chromatin 3D organization, affecting enhancer-promoter interactions and promoting growth. Here, we provide an integrative approach, combining chromatin conformation, accessibility, and transcription analysis, validated by in silico and CRISPR-interference screens, to identify relevant 3D topologies in pediatric T cell leukemia (T-ALL and ETP-ALL). We characterize 3D hubs as regulatory centers for oncogenes and disease markers, linking them to biological processes like cell division, inflammation, and stress response. Single-cell mapping reveals heterogeneous gene activation in discrete epigenetic clones, aiding in patient stratification for relapse risk after chemotherapy. Finally, we identify MYB as a 3D hub regulator in leukemia cells and show that the targeting of key regulators leads to hub dissolution, thereby providing a novel and effective anti-leukemic strategy. Overall, our work demonstrates the relevance of studying oncogenic 3D hubs to better understand cancer biology and tumor heterogeneity and to propose novel therapeutic strategies.

ParSite is a multicolor DNA labeling system that allows for simultaneous imaging of triple genomic loci in living cells

The organization of the human genome in space and time is critical for transcriptional regulation and cell fate determination. However, robust methods for tracking genome organization or genomic interactions over time in living cells are lacking. Here, we developed a multicolor DNA labeling system, ParSite, to simultaneously track triple genomic loci in the U2OS cells. The tricolor ParSite system is derived from the T. thermophilus ParB/ParSc (TtParB/ParSc) system by rational design. We mutated the interface between TtParB and ParSc and generated a new pair of TtParBm and ParSm for genomic DNA labeling. The insertions of 16 base-pair palindromic ParSc and ParSm into genomic loci allow dual-color DNA imaging in living cells. A pair of genomic loci labeled by ParSite could be colocalized with p53-binding protein 1 (53BP1) in response to CRISPR/Cas9-mediated double-strand breaks (DSBs). The ParSite permits tracking promoter and terminator dynamics of the APP gene, which spans 290 kilobases in length. Intriguingly, the hybrid ParS (ParSh) of half-ParSc and half-ParSm enables for the visualization of a third locus independent of ParSc or ParSm. We simultaneously labeled 3 loci with a genomic distance of 36, 89, and 352 kilobases downstream the C3 repeat locus, respectively. In sum, the ParSite is a robust DNA labeling system for tracking multiple genomic loci in space and time in living cells.

The cells are all-right: Regulation of the Lefty genes by separate enhancers in mouse embryonic stem cells

Enhancers play a critical role in regulating precise gene expression patterns essential for development and cellular identity; however, how gene-enhancer specificity is encoded within the genome is not clearly defined. To investigate how this specificity arises within topologically associated domains (TAD), we performed allele-specific genome editing of sequences surrounding the Lefty1 and Lefty2 paralogs in mouse embryonic stem cells. The Lefty genes arose from a tandem duplication event and these genes interact with each other in chromosome conformation capture assays which place these genes within the same TAD. Despite their physical proximity, we demonstrate that these genes are primarily regulated by separate enhancer elements. Through CRISPR-Cas9 mediated deletions to remove the intervening chromatin between the Lefty genes, we reveal a distance-dependent dosage effect of the Lefty2 enhancer on Lefty1 expression. These findings indicate a role for chromatin distance in insulating gene expression domains in the Lefty locus in the absence of architectural insulation.

Sequestration of gene products by decoys enhances precision in the timing of intracellular events

Expressed gene products often interact ubiquitously with binding sites at nucleic acids and macromolecular complexes, known as decoys. The binding of transcription factors (TFs) to decoys can be crucial in controlling the stochastic dynamics of gene expression. Here, we explore the impact of decoys on the timing of intracellular events, as captured by the time taken for the levels of a given TF to reach a critical threshold level, known as the first passage time (FPT). Although nonlinearity introduced by binding makes exact mathematical analysis challenging, employing suitable approximations and reformulating FPT in terms of an alternative variable, we analytically assess the impact of decoys. The stability of the decoy-bound TFs against degradation impacts FPT statistics crucially. Decoys reduce noise in FPT, and stable decoy-bound TFs offer greater timing precision with less expression cost than their unstable counterparts. Interestingly, when both bound and free TFs decay at the same rate, decoy binding does not directly alter FPT noise. We verify these results by performing exact stochastic simulations. These results have important implications for the precise temporal scheduling of events involved in the functioning of biomolecular clocks, development processes, cell-cycle control, and cell-size homeostasis.

Enhancers: A Focus on Synthetic Biology and Correlated Gene Expression

Enhancers are central for the regulation of metazoan transcription but have proven difficult to study, primarily due to a myriad of interdependent variables shaping their activity. Consequently, synthetic biology has emerged as the main approach for dissecting mechanisms of enhancer function. We start by reviewing simple but highly parallel reporter assays, which have been successful in quantifying the complexity of the activator/coactivator mechanisms at enhancers. We then describe studies that examine how enhancers function in the genomic context and in combination with other enhancers, revealing that they activate genes through a variety of different mechanisms, working together as a system. Here, we primarily focus on synthetic reporter genes that can quantify the dynamics of enhancer biology through time. We end by considering the consequences of having many genes and enhancers within a 'local environment', which we believe leads to correlated gene expression and likely reports on the general principles of enhancer biology.

Single-cell new RNA sequencing reveals principles of transcription at the resolution of individual bursts

Analyses of transcriptional bursting from single-cell RNA-sequencing data have revealed patterns of variation and regulation in the kinetic parameters that could be inferred. Here we profiled newly transcribed (4-thiouridine-labelled) RNA across 10,000 individual primary mouse fibroblasts to more broadly infer bursting kinetics and coordination. We demonstrate that inference from new RNA profiles could separate the kinetic parameters that together specify the burst size, and that the synthesis rate (and not the transcriptional off rate) controls the burst size. Importantly, transcriptome-wide inference of transcriptional on and off rates provided conclusive evidence that RNA polymerase II transcribes genes in bursts. Recent reports identified examples of transcriptional co-bursting, yet no global analyses have been performed. The deep new RNA profiles we generated with allelic resolution demonstrated that co-bursting rarely appears more frequently than expected by chance, except for certain gene pairs, notably paralogues located in close genomic proximity. Altogether, new RNA single-cell profiling critically improves the inference of transcriptional bursting and provides strong evidence for independent transcriptional bursting of mammalian genes.

A Compendium of G-Flipon Biological Functions That Have Experimental Validation

As with all new fields of discovery, work on the biological role of G-quadruplexes (GQs) has produced a number of results that at first glance are quite baffling, sometimes because they do not fit well together, but mostly because they are different from commonly held expectations. Like other classes of flipons, those that form G-quadruplexes have a repeat sequence motif that enables the fold. The canonical DNA motif (G3N1-7)3G3, where N is any nucleotide and G is guanine, is a feature that is under active selection in avian and mammalian genomes. The involvement of G-flipons in genome maintenance traces back to the invertebrate Caenorhabditis elegans and to ancient DNA repair pathways. The role of GQs in transcription is supported by the observation that yeast Rap1 protein binds both B-DNA, in a sequence-specific manner, and GQs, in a structure-specific manner, through the same helix. Other sequence-specific transcription factors (TFs) also engage both conformations to actuate cellular transactions. Noncoding RNAs can also modulate GQ formation in a sequence-specific manner and engage the same cellular machinery as localized by TFs, linking the ancient RNA world with the modern protein world. The coevolution of noncoding RNAs and sequence-specific proteins is supported by studies of early embryonic development, where the transient formation of G-quadruplexes coordinates the epigenetic specification of cell fate.

Transcriptional bursting dynamics in gene expression

Gene transcription is a stochastic process that occurs in all organisms. Transcriptional bursting, a critical molecular dynamics mechanism, creates significant heterogeneity in mRNA and protein levels. This heterogeneity drives cellular phenotypic diversity. Currently, the lack of a comprehensive quantitative model limits the research on transcriptional bursting. This review examines various gene expression models and compares their strengths and weaknesses to guide researchers in selecting the most suitable model for their research context. We also provide a detailed summary of the key metrics related to transcriptional bursting. We compared the temporal dynamics of transcriptional bursting across species and the molecular mechanisms influencing these bursts, and highlighted the spatiotemporal patterns of gene expression differences by utilizing metrics such as burst size and burst frequency. We summarized the strategies for modeling gene expression from both biostatistical and biochemical reaction network perspectives. Single-cell sequencing data and integrated multiomics approaches drive our exploration of cutting-edge trends in transcriptional bursting mechanisms. Moreover, we examined classical methods for parameter estimation that help capture dynamic parameters in gene expression data, assessing their merits and limitations to facilitate optimal parameter estimation. Our comprehensive summary and review of the current transcriptional burst dynamics theories provide deeper insights for promoting research on the nature of cell processes, cell fate determination, and cancer diagnosis.

Transcriptional bursting: from fundamentals to novel insights

Transcription occurs as irregular bursts in a very wide range of systems, including numerous different species and many genes within these. In this review, we examine the underlying theories, discuss how these relate to experimental measurements, and explore some of the discrepancies that have emerged among various studies. Finally, we consider more recent works that integrate novel concepts, such as the involvement of biomolecular condensates in enhancer-promoter interactions and their effects on the dynamics of transcriptional bursting.

Transcriptional noise, gene activation, and roles of SAGA and Mediator Tail measured using nucleotide recoding single-cell RNA-seq

We describe a time-resolved nascent single-cell RNA sequencing (RNA-seq) approach that measures gene-specific transcriptional noise and the fraction of active genes in S. cerevisiae. Most genes are expressed with near-constitutive behavior, while a subset of genes show high mRNA variance suggestive of transcription bursting. Transcriptional noise is highest in the cofactor/coactivator-redundant (CR) gene class (dependent on both SAGA and TFIID) and strongest in TATA-containing CR genes. Using this approach, we also find that histone gene transcription switches from a low-level, low-noise constitutive mode during M and M/G1 to an activated state in S phase that shows both an increase in the fraction of active promoters and a switch to a noisy and bursty transcription mode. Rapid depletion of cofactors SAGA and MED Tail indicates that both factors play an important role in stimulating the fraction of active promoters at CR genes, with a more modest role in transcriptional noise.

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

IVEA: an integrative variational Bayesian inference method for predicting enhancer-gene regulatory interactions

Motivation

Enhancers play critical roles in cell-type-specific transcriptional control. Despite the identification of thousands of candidate enhancers, unravelling their regulatory relationships with their target genes remains challenging. Therefore, computational approaches are needed to accurately infer enhancer-gene regulatory relationships.

Results

In this study, we propose a new method, IVEA, that predicts enhancer-gene regulatory interactions by estimating promoter and enhancer activities. Its statistical model is based on the gene regulatory mechanism of transcriptional bursting, which is characterized by burst size and frequency controlled by promoters and enhancers, respectively. Using transcriptional readouts, chromatin accessibility, and chromatin contact data as inputs, promoter and enhancer activities were estimated using variational Bayesian inference, and the contribution of each enhancer-promoter pair to target gene transcription was calculated. Our analysis demonstrates that the proposed method can achieve high prediction accuracy and provide biologically relevant enhancer-gene regulatory interactions.

Availability and implementation

The IVEA code is available on GitHub at https://github.com/yasumasak/ivea. The publicly available datasets used in this study are described in Supplementary Table S4.

3D chromatin architecture, BRD4, and Mediator have distinct roles in regulating genome-wide transcriptional bursting and gene network

Discontinuous transcription is evolutionarily conserved and a fundamental feature of gene regulation; yet, the exact mechanisms underlying transcriptional bursting are unresolved. Analyses of bursting transcriptome-wide have focused on the role of cis-regulatory elements, but other factors that regulate this process remain elusive. We applied mathematical modeling to single-cell RNA sequencing data to infer bursting dynamics transcriptome-wide under multiple conditions to identify possible molecular mechanisms. We found that Mediator complex subunit 26 (MED26) primarily regulates frequency, MYC regulates burst size, while cohesin and Bromodomain-containing protein 4 (BRD4) can modulate both. Despite comparable effects on RNA levels among these perturbations, acute depletion of MED26 had the most profound impact on the entire gene regulatory network, acting downstream of chromatin spatial architecture and without affecting TATA box-binding protein (TBP) recruitment. These results indicate that later steps in the initiation of transcriptional bursts are primary nodes for integrating gene networks in single cells.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.

Programmable promoter editing for precise control of transgene expression

Subtle changes in gene expression direct cells to distinct cellular states. Identifying and controlling dose-dependent transgenes require tools for precisely titrating expression. To this end, we developed a highly modular, extensible framework called DIAL for building editable promoters that allow for fine-scale, heritable changes in transgene expression. Using DIAL, we increase expression by recombinase-mediated excision of spacers between the binding sites of a synthetic zinc finger transcription factor and the core promoter. By nesting varying numbers and lengths of spacers, DIAL generates a tunable range of unimodal setpoints from a single promoter. Through small-molecule control of transcription factors and recombinases, DIAL supports temporally defined, user-guided control of transgene expression that is extensible to additional transcription factors. Lentiviral delivery of DIAL generates multiple setpoints in primary cells and iPSCs. As promoter editing generates stable states, DIAL setpoints are heritable, facilitating mapping of transgene levels to phenotypes. The DIAL framework opens new opportunities for tailoring transgene expression and improving the predictability and performance of gene circuits across diverse applications.

Gene activity fully predicts transcriptional bursting dynamics

Transcription commonly occurs in bursts, with alternating productive (ON) and quiescent (OFF) periods, governing mRNA production rates. Yet, how transcription is regulated through bursting dynamics remains unresolved. Here, we conduct real-time measurements of endogenous transcriptional bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional activities in early fly embryos, we uncover stringent relationships between bursting parameters. Specifically, we find that the durations of ON and OFF periods are linked. Regardless of the developmental stage or body-axis position, gene activity levels predict individual alleles' average ON and OFF periods. Lowly transcribing alleles predominantly modulate OFF periods (burst frequency), while highly transcribing alleles primarily tune ON periods (burst size). These relationships persist even under perturbations of cis-regulatory elements or trans-factors and account for bursting dynamics measured in other species. Our results suggest a novel mechanistic constraint governing bursting dynamics rather than a modular control of distinct parameters by distinct regulatory processes.

N-terminal tagging of RNA Polymerase II shapes transcriptomes more than C-terminal alterations

RNA polymerase II (Pol II) has a C-terminal domain (CTD) that is unstructured, consisting of a large number of heptad repeats, and whose precise function remains unclear. Here, we investigate how altering the CTD's length and fusing it with protein tags affects transcriptional output on a genome-wide scale in mammalian cells at single-cell resolution. While transcription generally appears to occur in burst-like fashion, where RNA is predominantly made during short bursts of activity that are interspersed with periods of transcriptional silence, the CTD's role in shaping these dynamics seems gene-dependent; global patterns of bursting appear mostly robust to CTD alterations. Introducing protein tags with defined structures to the N terminus cause transcriptome-wide effects, however. We find the type of tag to dominate characteristics of the resulting transcriptomes. This is possibly due to Pol II-interacting factors, including non-coding RNAs, whose expression correlates with the tags. Proteins involved in liquid-liquid phase separation appear prominently.

Identification of molecular determinants of gene-specific bursting patterns by high-throughput imaging screens

Stochastic transcriptional bursting is a universal property of active genes. While different genes exhibit distinct bursting patterns, the molecular mechanisms for gene-specific stochastic bursting are largely unknown. We have developed and applied a high-throughput-imaging based screening strategy to identify cellular factors and molecular mechanisms that determine the bursting behavior of human genes. Focusing on epigenetic regulators, we find that protein acetylation is a strong acute modulator of burst frequency, burst size and heterogeneity of bursting. Acetylation globally affects the Off-time of genes but has gene-specific effects on the On-time. Yet, these effects are not strongly linked to promoter acetylation, which do not correlate with bursting properties, and forced promoter acetylation has variable effects on bursting. Instead, we demonstrate acetylation of the Integrator complex as a key determinant of gene bursting. Specifically, we find that elevated Integrator acetylation decreases bursting frequency. Taken together our results suggest a prominent role of non-histone proteins in determining gene bursting properties, and they identify histone-independent acetylation of a transcription cofactor as an allosteric modulator of bursting via a far-downstream bursting checkpoint.

Analysis of developmental gene expression using smFISH and in silico staging of C. elegans embryos

Regulation of transcription during embryogenesis is key to development and differentiation. To study transcript expression throughout Caenorhabditis elegans embryogenesis at single-molecule resolution, we developed a high-throughput single-molecule fluorescence in situ hybridization (smFISH) method that relies on computational methods to developmentally stage embryos and quantify individual mRNA molecules in single embryos. We applied our system to sdc-2, a zygotically transcribed gene essential for hermaphrodite development and dosage compensation. We found that sdc-2 is rapidly activated during early embryogenesis by increasing both the number of mRNAs produced per transcription site and the frequency of sites engaged in transcription. Knockdown of sdc-2 and dpy-27, a subunit of the dosage compensation complex (DCC), increased the number of active transcription sites for the X chromosomal gene dpy-23 but not the autosomal gene mdh-1, suggesting that the DCC reduces the frequency of dpy-23 transcription. The temporal resolution from in silico staging of embryos showed that the deletion of a single DCC recruitment element near the dpy-23 gene causes higher dpy-23 mRNA expression after the start of dosage compensation, which could not be resolved using mRNAseq from mixed-stage embryos. In summary, we have established a computational approach to quantify temporal regulation of transcription throughout C. elegans embryogenesis and demonstrated its potential to provide new insights into developmental gene regulation.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Biallelic non-productive enhancer-promoter interactions precede imprinted expression of Kcnk9 during mouse neural commitment

It is only partially understood how constitutive allelic methylation at imprinting control regions (ICRs) interacts with other regulation levels to drive timely parental allele-specific expression along large imprinted domains. The Peg13-Kcnk9 domain is an imprinted domain with important brain functions. To gain insights into its regulation during neural commitment, we performed an integrative analysis of its allele-specific epigenetic, transcriptomic, and cis-spatial organization using a mouse stem cell-based corticogenesis model that recapitulates the control of imprinted gene expression during neurodevelopment. We found that, despite an allelic higher-order chromatin structure associated with the paternally CTCF-bound Peg13 ICR, enhancer-Kcnk9 promoter contacts occurred on both alleles, although they were productive only on the maternal allele. This observation challenges the canonical model in which CTCF binding isolates the enhancer and its target gene on either side and suggests a more nuanced role for allelic CTCF binding at some ICRs.

Epigenomic states contribute to coordinated allelic transcriptional bursting in iPSC reprogramming

Two alleles of a gene can be transcribed independently or coordinatedly, which can lead to temporal expression heterogeneity with potentially distinct impacts on cell fate. Here, we profiled genome-wide allelic transcriptional burst kinetics during the reprogramming of MEF to induced pluripotent stem cells. We show that the degree of coordination of allelic bursting differs among genes, and alleles of many reprogramming-related genes burst in a highly coordinated fashion. Notably, we show that the chromatin accessibility of the two alleles of highly coordinated genes is similar, unlike the semi-coordinated or independent genes, suggesting the degree of coordination of allelic bursting is linked to allelic chromatin accessibility. Consistently, we show that many transcription factors have differential binding affinity between alleles of semi-coordinated or independent genes. We show that highly coordinated genes are enriched with chromatin accessibility regulators such as H3K4me3, H3K4me1, H3K36me3, H3K27ac, histone variant H3.3, and BRD4. Finally, we demonstrate that enhancer elements are highly enriched in highly coordinated genes. Our study demonstrates that epigenomic states contribute to coordinated allelic bursting to fine-tune gene expression during induced pluripotent stem cell reprogramming.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

Time will tell: comparing timescales to gain insight into transcriptional bursting

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

Super-enhancers and the super-enhancer reader BRD4: tumorigenic factors and therapeutic targets

Transcriptional super-enhancers and the BET bromodomain protein BRD4 are emerging as critical drivers of tumorigenesis and therapeutic targets. Characterized by substantial accumulation of histone H3 lysine 27 acetylation (H3K27ac) signals at the loci of cell identity genes and critical oncogenes, super-enhancers are recognized, bound and activated by BRD4, resulting in considerable oncogene over-expression, malignant transformation, cancer cell proliferation, survival, tumor initiation and progression. Small molecule compound BRD4 BD1 and BD2 bromodomain inhibitors block BRD4 binding to super-enhancers, suppress oncogene transcription and expression, reduce cancer cell proliferation and survival, and repress tumor progression in a variety of cancer types. Like other targeted therapy agents, BRD4 inhibitors show moderate anticancer effects on their own, and exert synergistic anticancer effects in vitro and in preclinical models, when combined with other anticancer agents including CDK7 inhibitors, CBP/p300 inhibitors and histone deacetylase inhibitors. More recently, BRD4 BD2 bromodomain selective inhibitors, proteolysis-targeting chimera (PROTAC) BRD4 protein degraders, and dual BRD4 and CBP/p300 bromodomain co-inhibitors have been developed and shown better anticancer efficacy and/or safety profile. Importantly, more than a dozen BRD4 inhibitors have entered clinical trials in patients with cancer of various organ origins. In summary, super-enhancers and their reader BRD4 are critical tumorigenic drivers, and BRD4 BD1 and BD2 bromodomain inhibitors, BRD4 BD2 bromodomain selective inhibitors, PROTAC BRD4 protein degraders, and dual BRD4 and CBP/p300 bromodomain co-inhibitors are promising novel anticancer agents for clinical translation.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

Effects of microRNA-mediated negative feedback on gene expression noise

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally in eukaryotes by binding with target mRNAs and preventing translation. miRNA-mediated feedback motifs are ubiquitous in various genetic networks that control cellular decision making. A key question is how such a feedback mechanism may affect gene expression noise. To answer this, we have developed a mathematical model to study the effects of a miRNA-dependent negative-feedback loop on mean expression and noise in target mRNAs. Combining analytics and simulations, we show the existence of an expression threshold demarcating repressed and expressed regimes in agreement with earlier studies. The steady-state mRNA distributions are bimodal near the threshold, where copy numbers of mRNAs and miRNAs exhibit enhanced anticorrelated fluctuations. Moreover, variation of negative-feedback strength shifts the threshold locations and modulates the noise profiles. Notably, the miRNA-mRNA binding affinity and feedback strength collectively shape the bimodality. We also compare our model with a direct auto-repression motif, where a gene produces its own repressor. Auto-repression fails to produce bimodal mRNA distributions as found in miRNA-based indirect repression, suggesting the crucial role of miRNAs in creating phenotypic diversity. Together, we demonstrate how miRNA-dependent negative feedback modifies the expression threshold and leads to a broader parameter regime of bimodality compared to the no-feedback case.

Ultrafine mapping of chromosome conformation at hundred basepair resolution reveals regulatory genome architecture

The resolution limit of chromatin conformation capture methodologies (3Cs) has restrained their application in detection of fine-level chromatin structure mediated by cis-regulatory elements (CREs). Here, we report two 3C-derived methods, Tri-4C and Tri-HiC, which utilize multirestriction enzyme digestions for ultrafine mapping of targeted and genome-wide chromatin interaction, respectively, at up to one hundred basepair resolution. Tri-4C identified CRE loop interaction networks and quantitatively revealed their alterations underlying dynamic gene control. Tri-HiC uncovered global fine-gauge regulatory interaction networks, identifying >20-fold more enhancer:promoter (E:P) loops than in situ Hi-C. In addition to vastly improved identification of subkilobase-sized E:P loops, Tri-HiC also uncovered interaction stripes and contact domain insulation from promoters and enhancers, revealing their loop extrusion behaviors resembling the topologically associating domain boundaries. Tri-4C and Tri-HiC provide robust approaches to achieve the high-resolution interactome maps required for characterizing fine-gauge regulatory chromatin interactions in analysis of development, homeostasis, and disease.

Predicting active enhancers with DNA methylation and histone modification

BACKGROUND

Enhancers play a crucial role in gene regulation, and some active enhancers produce noncoding RNAs known as enhancer RNAs (eRNAs) bi-directionally. The most commonly used method for detecting eRNAs is CAGE-seq, but the instability of eRNAs in vivo leads to data noise in sequencing results. Unfortunately, there is currently a lack of research focused on the noise inherent in CAGE-seq data, and few approaches have been developed for predicting eRNAs. Bridging this gap and developing widely applicable eRNA prediction models is of utmost importance.

RESULTS

In this study, we proposed a method to reduce false positives in the identification of eRNAs by adjusting the statistical distribution of expression levels. We also developed eRNA prediction models using joint gene expressions, DNA methylation, and histone modification. These models achieved impressive performance with an AUC value of approximately 0.95 for intra-cell prediction and 0.9 for cross-cell prediction.

CONCLUSIONS

Our method effectively attenuates the noise generated by stochastic RNA production, resulting in more accurate detection of eRNAs. Furthermore, our eRNA prediction model exhibited significant accuracy in both intra-cell and cross-cell validation, highlighting its robustness and potential application in various cellular contexts.

Regulatory Mechanisms for Transcriptional Bursting Revealed by an Event-Based Model

Gene transcription often occurs in discrete bursts, and it can be difficult to deduce the underlying regulatory mechanisms for transcriptional bursting with limited experimental data. Here, we categorize numerous states of single eukaryotic genes and identify 6 essential transcriptional events, each comprising a series of state transitions; transcriptional bursting is characterized as a sequence of 4 events, capable of being organized in various configurations, in addition to the beginning and ending events. By associating transcriptional kinetics with mean durations and recurrence probabilities of the events, we unravel how transcriptional bursting is modulated by various regulators including transcription factors. Through analytical derivation and numerical simulation, this study reveals key state transitions contributing to transcriptional sensitivity and specificity, typical characteristics of burst profiles, global constraints on intrinsic transcriptional noise, major regulatory modes in individual genes and across the genome, and requirements for fast gene induction upon stimulation. It is illustrated how biochemical reactions on different time scales are modulated to separately shape the durations and ordering of the events. Our results suggest that transcriptional patterns are essentially controlled by a shared set of transcriptional events occurring under specific promoter architectures and regulatory modes, the number of which is actually limited.

Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

BACKGROUND

Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states.

METHODS

We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing.

RESULTS

We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease.

CONCLUSION

This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

Nuclear speckles regulate HIF-2α programs and correlate with patient survival in kidney cancer

Nuclear speckles are membrane-less bodies within the cell nucleus enriched in RNA biogenesis, processing, and export factors. In this study we investigated speckle phenotype variation in human cancer, finding a reproducible speckle signature, based on RNA expression of speckle-resident proteins, across >20 cancer types. Of these, clear cell renal cell carcinoma (ccRCC) exhibited a clear correlation between the presence of this speckle expression signature, imaging-based speckle phenotype, and clinical outcomes. ccRCC is typified by hyperactivation of the HIF-2α transcription factor, and we demonstrate here that HIF-2α drives physical association of a select subset of its target genes with nuclear speckles. Disruption of HIF-2α-driven speckle association via deletion of its speckle targeting motifs (STMs)-defined in this study-led to defective induction of speckle-associating HIF-2α target genes without impacting non-speckle-associating HIF-2α target genes. We further identify the RNA export complex, TREX, as being specifically altered in speckle signature, and knockdown of key TREX component, ALYREF, also compromises speckle-associated gene expression. By integrating tissue culture functional studies with tumor genomic and imaging analysis, we show that HIF-2α gene regulatory programs are impacted by specific manipulation of speckle phenotype and by abrogation of speckle targeting abilities of HIF-2α. These findings suggest that, in ccRCC, a key biological function of nuclear speckles is to modulate expression of a specific subset of HIF-2α-regulated target genes that, in turn, influence patient outcomes. We also identify STMs in other transcription factors, suggesting that DNA-speckle targeting may be a general mechanism of gene regulation.

LAST-seq: single-cell RNA sequencing by direct amplification of single-stranded RNA without prior reverse transcription and second-strand synthesis

Existing single-cell RNA sequencing (scRNA-seq) methods rely on reverse transcription (RT) and second-strand synthesis (SSS) to convert single-stranded RNA into double-stranded DNA prior to amplification, with the limited RT/SSS efficiency compromising RNA detectability. Here, we develop a new scRNA-seq method, Linearly Amplified Single-stranded-RNA-derived Transcriptome sequencing (LAST-seq), which directly amplifies the original single-stranded RNA molecules without prior RT/SSS. LAST-seq offers a high single-molecule capture efficiency and a low level of technical noise for single-cell transcriptome analyses. Using LAST-seq, we characterize transcriptional bursting kinetics in human cells, revealing a role of topologically associating domains in transcription regulation.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

Enhancer dynamics: Unraveling the mechanism of transcriptional bursting

Transcriptional bursting is a prevalent feature of gene expression. The transient assembly of transcription factor clusters at regulatory DNAs is critical to control bursting dynamics.

Defining global strategies to improve outcomes in sickle cell disease: a Lancet Haematology Commission

All over the world, people with sickle cell disease (an inherited condition) have premature deaths and preventable severe chronic complications, which considerably affect their quality of life, career progression, and financial status. In addition, these people are often affected by stigmatisation or structural racism, which can contribute to stress and poor mental health. Inequalities affecting people with sickle cell disease are also reflected in the distribution of the disease—mainly in sub-Saharan Africa, India, and the Caribbean—whereas interventions, clinical trials, and funding are mostly available in North America, Europe, and the Middle East. Although some of these characteristics also affect people with other genetic diseases, the fate of people with sickle cell disease seems to be particularly unfair. Simple, effective interventions to reduce the mortality and morbidity associated with sickle cell disease are available. The main obstacle preventing better outcomes in this condition, which is a neglected disease, is associated with inequalities impacting the patient populations. The aim of this Commission is to highlight the problems associated with sickle cell disease and to identify achievable goals to improve outcomes both in the short and long term.

The ambition for the management of people with sickle cell disease is that curative treatments become available to every person with the condition. Although this would have seemed unrealistic a decade ago, developments in gene therapy make this potentially achievable, albeit in the distant future. Until these curative technologies are fully developed and become widely available, health-care professionals (with the support of policy makers, funders, etc) should make sure that a minimum standard of care (including screening, prophylaxis against infection, acute medical care, safe blood transfusion, and hydroxyurea) is available to all patients.

In considering what needs to be achieved to reduce the global burden of sickle cell disease and improve the quality of life of patients, this Commission focuses on five key areas: the epidemiology of sickle cell disease (Section 1 ); screening and prevention (Section 2 ); established and emerging treatments for the management of the disease (Section 3 ); cellular therapies with curative potential (Section 4 ); and training and education needs (Section 5 ). As clinicians, researchers, and patients, our objective to reduce the global burden of sickle cell disease aligns with wider public health aims to reduce inequalities, improve health for all, and develop personalised treatment options. We have observed in the past few years some long-awaited momentum following the development of innovative point-of-care testing devices, new approved drugs, and emerging curative options. Reducing the burden of sickle cell disease will require substantial financial and political commitment, but it will impact the lives of millions of patients and families worldwide and the lessons learned in achieving this goal would unarguably benefit society as a whole.

Live-cell imaging of chromatin contacts opens a new window into chromatin dynamics

Our understanding of the organization of the chromatin fiber within the cell nucleus has made great progress in the last few years. High-resolution techniques based on next-generation sequencing as well as optical imaging that can investigate chromatin conformations down to the single cell level have revealed that chromatin structure is highly heterogeneous at the level of the individual allele. While TAD boundaries and enhancer-promoter pairs emerge as hotspots of 3D proximity, the spatiotemporal dynamics of these different types of chromatin contacts remain largely unexplored. Investigation of chromatin contacts in live single cells is necessary to close this knowledge gap and further enhance the current models of 3D genome organization and enhancer-promoter communication. In this review, we first discuss the potential of single locus labeling to study architectural and enhancer-promoter contacts and provide an overview of the available single locus labeling techniques such as FROS, TALE, CRISPR-dCas9 and ANCHOR, and discuss the latest developments and applications of these systems.