Transcription pausing regulates mouse embryonic stem cell differentiation

Kinetic organization of the genome revealed by ultra-resolution, multiscale live imaging

In the last decade, sequencing methods like Hi-C have made it clear the genome is intricately folded, and that this organization contributes significantly to the control of gene expression and thence cell fate and behavior. Single-cell DNA tracing microscopy and polymer physics-based simulations of genome folding have proposed these population-scale patterns arise from motor- driven, heterogeneous movement, rather than stable 3D genomic architecture, implying that motion, rather than structure, is key to understanding genome function. However, tools to directly observe this motion in vivo have been limited in coverage and resolution. Here we describe TRansposon Assisted Chromatin Kinetic Imaging Technology (TRACK-IT), which combines a suite of imaging and labeling improvements to achieve ultra-resolution in space and time, with self-mapping transposons to distribute labels across the chromosome, uncovering dynamic behaviors across four orders of magnitude of genomic separation. We find that sequences separated by sub-megabase distances, typically 200-500 nm of nanometers apart, can transition to close proximity in tens of seconds - faster than previously hypothesized. This rapid motion is dependent upon cohesin and is exhibited only within certain genomic domains. Domain borders act as kinetic impediments to this search process, substantially slowing the rate and frequency of the transition to proximity. The genomic separation-dependent scaling of the search time for cis-interactions within a domain violates predictions of diffusion, suggesting motor driven folding. This distinctive scaling is lost following cohesin depletion, replaced with a behavior consistent with diffusion. Finally, we found cohesin containing cells exhibited rare, processive movements, not seen in cohesin depleted cells. These processive trajectories exhibit extrusion rates of ∼2.7 kb/s across three distinct genomic intervals, faster than recent in vitro measurements and prior estimates from in vivo data. Taken together, these results reveal a genome in motion across multiple genomic and temporal scales, where motor-dependent extrusion divides the sequence, not into spatially separate domains, but into kinetically separated domains that experience accelerated local search.

Loss-of-Function Variants in SUPT5H as Modifying Factors in Beta-Thalassemia

It is well known that modifiers play a role in ameliorating or exacerbating disease phenotypes in patients and carriers of recessively inherited disorders such as sickle cell disease and thalassemia. Here, we give an overview of the literature concerning a recently described association in carriers of SUPT5H Loss-of-Function variants with a beta-thalassemia-like phenotype including the characteristic elevated levels of HbA2. That SUPT5H acts as modifier in beta-thalassemia carriers became evident from three reported cases in whom combined heterozygosity of SUPT5H and HBB gene variants was observed to resemble a mild beta-thalassemia intermedia phenotype. The different SUPT5H variants and hematologic parameters reported are collected and reviewed to provide insight into the possible effects on hematologic expression, as well as potential disease mechanisms in carriers and patients.

RNA polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

Controlled release of promoter-proximal paused RNA polymerase II (RNA Pol II) is crucial for gene regulation. However, studying RNA Pol II pausing is challenging, as pause-release factors are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H, which encodes SPT5, in individuals with β-thalassemia. During erythropoiesis in healthy human cells, cell cycle genes were highly paused as cells transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, RNA Pol II pause release was globally disrupted, and as cells began transitioning from progenitors to precursors, differentiation was delayed, accompanied by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, identifying a role for RNA Pol II pausing in temporally coordinating the cell cycle and erythroid differentiation.

Orphan quality control shapes network dynamics and gene expression

All eukaryotes require intricate protein networks to translate developmental signals into accurate cell fate decisions. Mutations that disturb interactions between network components often result in disease, but how the composition and dynamics of complex networks are established remains poorly understood. Here, we identify the E3 ligase UBR5 as a signaling hub that helps degrade unpaired subunits of multiple transcriptional regulators that act within a network centered on the c-Myc oncoprotein. Biochemical and structural analyses show that UBR5 binds motifs that only become available upon complex dissociation. By rapidly turning over unpaired transcription factor subunits, UBR5 establishes dynamic interactions between transcriptional regulators that allow cells to effectively execute gene expression while remaining receptive to environmental signals. We conclude that orphan quality control plays an essential role in establishing dynamic protein networks, which may explain the conserved need for protein degradation during transcription and offers opportunities to modulate gene expression in disease.

RNA Polymerase II pausing temporally coordinates cell cycle progression and erythroid differentiation

The controlled release of promoter-proximal paused RNA polymerase II (Pol II) into productive elongation is a major step in gene regulation. However, functional analysis of Pol II pausing is difficult because factors that regulate pause release are almost all essential. In this study, we identified heterozygous loss-of-function mutations in SUPT5H , which encodes SPT5, in individuals with β-thalassemia unlinked to HBB mutations. During erythropoiesis in healthy human cells, cell cycle genes were highly paused at the transition from progenitors to precursors. When the pathogenic mutations were recapitulated by SUPT5H editing, Pol II pause release was globally disrupted, and the transition from progenitors to precursors was delayed, marked by a transient lag in erythroid-specific gene expression and cell cycle kinetics. Despite this delay, cells terminally differentiate, and cell cycle phase distributions normalize. Therefore, hindering pause release perturbs proliferation and differentiation dynamics at a key transition during erythropoiesis, revealing a role for Pol II pausing in the temporal coordination between the cell cycle and differentiation.

Distinct patterns of RNA polymerase II and transcriptional elongation characterize mammalian genome activation

How transcription is regulated as development commences is fundamental to understand how the transcriptionally silent mature gametes are reprogrammed. The embryonic genome is activated for the first time during zygotic genome activation (ZGA). How RNA polymerase II (Pol II) and productive elongation are regulated during this process remains elusive. Here, we generate genome-wide maps of Serine 5 and Serine 2-phosphorylated Pol II during and after ZGA in mouse embryos. We find that both phosphorylated Pol II forms display similar distributions across genes during ZGA, with typical elongation enrichment of Pol II emerging after ZGA. Serine 2-phosphorylated Pol II occurs at genes prior to their activation, suggesting that Serine 2 phosphorylation may prime gene expression. Functional perturbations demonstrate that CDK9 and SPT5 are major ZGA regulators and that SPT5 prevents precocious activation of some genes. Overall, our work sheds molecular insights into transcriptional regulation at the beginning of mammalian development.

CTCF and cohesin promote focal detachment of DNA from the nuclear lamina

Background

Lamina-associated domains (LADs) are large genomic regions that are positioned at the nuclear lamina. It has remained largely unclear what drives the positioning and demarcation of LADs. Because the insulator protein CTCF is enriched at LAD borders, it was postulated that CTCF binding could position some LAD boundaries, possibly through its function in stalling cohesin and hence preventing cohesin invading into the LAD. To test this, we mapped genome–nuclear lamina interactions in mouse embryonic stem cells after rapid depletion of CTCF and other perturbations of cohesin dynamics.

Results

CTCF and cohesin contribute to a sharp transition in lamina interactions at LAD borders, while LADs are maintained after depletion of these proteins, also at borders marked by CTCF. CTCF and cohesin may thus reinforce LAD borders, but do not position these. CTCF binding sites within LADs are locally detached from the lamina and enriched for accessible DNA and active histone modifications. Remarkably, despite lamina positioning being strongly correlated with genome inactivity, this DNA remains accessible after the local detachment is lost following CTCF depletion. At a chromosomal scale, cohesin depletion and cohesin stabilization by depletion of the unloading factor WAPL quantitatively affect lamina interactions, indicative of perturbed chromosomal positioning in the nucleus. Finally, while H3K27me3 is locally enriched at CTCF-marked LAD borders, we find no evidence for an interplay between CTCF and H3K27me3 on lamina interactions.

Conclusions

These findings illustrate that CTCF and cohesin are not primary determinants of LAD patterns. Rather, these proteins locally modulate NL interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02754-3.

RNA polymerase II pausing in development: orchestrating transcription

The coordinated regulation of transcriptional networks underpins cellular identity and developmental progression. RNA polymerase II promoter-proximal pausing (Pol II pausing) is a prevalent mechanism by which cells can control and synchronize transcription. Pol II pausing regulates the productive elongation step of transcription at key genes downstream of a variety of signalling pathways, such as FGF and Nodal. Recent advances in our understanding of the Pol II pausing machinery and its role in transcription call for an assessment of these findings within the context of development. In this review, we discuss our current understanding of the molecular basis of Pol II pausing and its function during organismal development. By critically assessing the tools used to study this process we conclude that combining recently developed genomics approaches with refined perturbation systems has the potential to expand our understanding of Pol II pausing mechanistically and functionally in the context of development and beyond.

Transcription factor OTX2 silences the expression of cleavage embryo genes and transposable elements

Upon mammalian fertilization, zygotic genome activation (ZGA) and activation of transposable elements (TEs) occur in early embryos to establish totipotency and support embryogenesis. However, the molecular mechanisms controlling the expression of these genes in mammals remain poorly understood. The 2-cell-like population of mouse embryonic stem cells (mESCs) mimics cleavage-stage embryos with transient Dux activation. In this study, we demonstrated that deficiency of the transcription factor OTX2 stimulates the expression of ZGA genes in mESCs. Further analysis revealed that OTX2 is incorporated at the Dux locus with corepressors for transcriptional inhibition. We also found that OTX2 associates with TEs and silences the subtypes of TEs. Therefore, OTX2 protein plays an important role in ZGA and TE expression in mESCs to orchestrate the transcriptional network.

A nuclease- and bisulfite-based strategy captures strand-specific R-loops genome-wide

R-loops are three-stranded nucleic acid structures with essential roles in many nuclear processes. However, their unchecked accumulation is associated with genome instability and is observed in neurodevelopmental diseases and cancers. Genome-wide profiling of R-loops in normal and diseased cells can help identify locations of pathogenic R-loops and advance efforts to attenuate them. We present an antibody-independent R-loop detection strategy, BisMapR, that combines nuclease-based R-loop isolation with non-denaturing bisulfite chemistry to produce genome-wide profiles that retain strand information. BisMapR achieves greater resolution and is faster than existing strand-specific R-loop profiling strategies. In mouse embryonic stem cells, we apply BisMapR to find that gene promoters form R-loops in both directions and uncover a subset of active enhancers that, despite being bidirectionally transcribed, form R-loops exclusively on one strand. BisMapR reveals a previously unnoticed feature of active enhancers and provides a tool to systematically examine their mechanisms in gene expression.

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE⁺⁴⁷) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE⁺⁴⁷ significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE⁺⁴⁷ deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE⁺⁴⁷ organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

The Paf1 complex and P-TEFb have reciprocal and antagonist roles in maintaining multipotent neural crest progenitors

Multipotent progenitor populations are necessary for generating diverse tissue types during embryogenesis. We show the RNA polymerase-associated factor 1 complex (Paf1C) is required to maintain multipotent progenitors of the neural crest (NC) lineage in zebrafish. Mutations affecting each Paf1C component result in near-identical NC phenotypes; alyron mutant embryos carrying a null mutation in paf1 were analyzed in detail. In the absence of zygotic paf1 function, definitive premigratory NC progenitors arise but fail to maintain expression of the sox10 specification gene. The mutant NC progenitors migrate aberrantly and fail to differentiate appropriately. Blood and germ cell progenitor development is affected similarly. Development of mutant NC could be rescued by additional loss of positive transcription elongation factor b (P-TEFb) activity, a key factor in promoting transcription elongation. Consistent with the interpretation that inhibiting/delaying expression of some genes is essential for maintaining progenitors, mutant embryos lacking the CDK9 kinase component of P-TEFb exhibit a surfeit of NC progenitors and their derivatives. We propose Paf1C and P-TEFb act antagonistically to regulate the timing of the expression of genes needed for NC development.

RNA structure maps across mammalian cellular compartments

RNA structure is intimately connected to each step of gene expression. Recent advances have enabled transcriptome-wide maps of RNA secondary structure, termed RNA structuromes. However, previous whole-cell analyses lacked the resolution to unravel the landscape and also the regulatory mechanisms of RNA structural changes across subcellular compartments. Here we reveal the RNA structuromes in three compartments — chromatin, nucleoplasm and cytoplasm in human and mouse cells. The cytotopic structuromes substantially expand RNA structural information, and enable detailed investigation of the central role of RNA structure in linking transcription, translation, and RNA decay. We develop a resource to visualize the interplay of RNA-protein interactions, RNA modifications, and RNA structure, and predict both direct and indirect reader proteins of RNA modifications. We also validate a novel role of the RNA binding protein LIN28A as an N6-methyladenosine modification “anti-reader”. Our results highlight the dynamic nature of RNA structures and its functional significance in gene regulation.

NETSeq reveals heterogeneous nucleotide incorporation by RNA polymerase I

DNA sequence motifs that affect RNA polymerase transcription elongation are well studied in prokaryotic organisms and contribute directly to regulation of gene expression. Despite significant work on the regulation of eukaryotic transcription, the effect of DNA template sequence on RNA polymerase I (Pol I) transcription elongation remains unknown. In this study, we examined the effects of DNA sequence motifs on Pol I transcription elongation kinetics in vitro and in vivo. Specifically, we characterized how the spy rho-independent terminator motif from Escherichia coli directly affects Saccharomyces cerevisiae Pol I activity, demonstrating evolutionary conservation of sequence-specific effects on transcription. The insight gained from this analysis led to the identification of a homologous sequence in the ribosomal DNA of S. cerevisiae We then used native elongating transcript sequencing (NETSeq) to determine whether Pol I encounters pause-inducing sequences in vivo. We found hundreds of positions within the ribosomal DNA (rDNA) that reproducibly induce pausing in vivo. We also observed significantly lower Pol I occupancy at G residues in the rDNA, independent of other sequence context, indicating differential nucleotide incorporation rates for Pol I in vivo. These data demonstrate that DNA template sequence elements directly influence Pol I transcription elongation. Furthermore, we have developed the necessary experimental and analytical methods to investigate these perturbations in living cells going forward.