Condensin controls recruitment of RNA polymerase II to achieve nematode X-chromosome dosage compensation

A nematode-specific ribonucleoprotein complex mediates interactions between the major nematode spliced leader snRNP and its target pre-mRNAs

Spliced leader trans-splicing of pre-mRNAs is a critical step in the gene expression of many eukaryotes. How the spliced leader RNA and its target transcripts are brought together to form the trans-spliceosome remains an important unanswered question. Using immunoprecipitation followed by protein analysis via mass spectrometry and RIP-Seq, we show that the nematode-specific proteins, SNA-3 and SUT-1, form a complex with a set of enigmatic non-coding RNAs, the SmY RNAs. Our work redefines the SmY snRNP and shows for the first time that it is essential for nematode viability and is involved in spliced leader trans-splicing. SNA-3 and SUT-1 are associated with the 5' ends of most, if not all, nascent capped RNA polymerase II transcripts, and they also interact with components of the major nematode spliced leader (SL1) snRNP. We show that depletion of SNA-3 impairs the co-immunoprecipitation between one of the SL1 snRNP components, SNA-2, and several core spliceosomal proteins. We thus propose that the SmY snRNP recruits the SL1 snRNP to the 5' ends of nascent pre-mRNAs, an instrumental step in the assembly of the trans-spliceosome.

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.

Cohesin mediated loop extrusion from active enhancers form chromatin jets in C. elegans

In mammals, cohesin and CTCF organize the 3D genome into topologically associated domains (TADs) to regulate communication between cis-regulatory elements. However, many organisms, including S. cerevisiae , C. elegans , and A. thaliana lack CTCF. Here, we use C. elegans as a model to investigate the function of cohesin in 3D genome organization in an animal without CTCF. We use auxin-inducible degradation to acutely deplete SMC-3 or its negative regulator WAPL-1 from somatic cells. Using Hi-C data, we identify a cohesin-dependent 3D genome organization feature called chromatin jets (aka fountains). These are population average reflections of DNA loops that are ∼20-40 kb in scale and often cover a few transcribed genes. The jets emerge from NIPBL occupied segments, and the trajectory of the jets coincides with cohesin binding. Cohesin translocation from jet origins depends on a fully intact complex and is extended upon WAPL-1 depletion. Hi-C results support the idea that cohesin is preferentially loaded at NIPBL occupied sites and loop extrudes in an effectively two-sided manner. The location of putative loading sites coincide with active enhancers and the strength of chromatin jet pattern correlates with transcription. Hi-C analyses upon WAPL-1 depletion reveal unequal loop extrusion processivity on each side and stalling due to cohesin molecules colliding. Compared to mammalian systems, average processivity of C. elegans cohesin is ∼10-fold shorter and NIPBL binding does not depend on cohesin. We conclude that the processivity of cohesin scales with genome size; and regardless of CTCF presence, preferential loading of cohesin at enhancers is a conserved mechanism of genome organization that regulates the interaction of gene regulatory elements in 3D.

An Epigenomic fingerprint of human cancers by landscape interrogation of super enhancers at the constituent level

Super enhancers (SE), large genomic elements that activate transcription and drive cell identity, have been found with cancer-specific gene regulation in human cancers. Recent studies reported the importance of understanding the cooperation and function of SE internal components, i.e., the constituent enhancers (CE). However, there are no pan-cancer studies to identify cancer-specific SE signatures at the constituent level. Here, by revisiting pan-cancer SE activities with H3K27Ac ChIP-seq datasets, we report fingerprint SE signatures for 28 cancer types in the NCI-60 cell panel. We implement a mixture model to discriminate active CEs from inactive CEs by taking into consideration ChIP-seq variabilities between cancer samples and across CEs. We demonstrate that the model-based estimation of CE states provides improved functional interpretation of SE-associated regulation. We identify cancer-specific CEs by balancing their active prevalence with their capability of encoding cancer type identities. We further demonstrate that cancer-specific CEs have the strongest per-base enhancer activities in independent enhancer sequencing assays, suggesting their importance in understanding critical SE signatures. We summarize fingerprint SEs based on the cancer-specific statuses of their component CEs and build an easy-to-use R package to facilitate the query, exploration, and visualization of fingerprint SEs across cancers.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

The NELF pausing checkpoint mediates the functional divergence of Cdk9

Promoter-proximal pausing by RNA Pol II is a rate-determining step in gene transcription that is hypothesized to be a prominent point at which regulatory factors act. The pausing factor NELF is known to induce and stabilize pausing, but not all kinds of pausing are NELF-mediated. Here, we find that NELF-depleted Drosophila melanogaster cells functionally recapitulate the NELF-independent pausing we previously observed in fission yeast (which lack NELF). Critically, only NELF-mediated pausing establishes a strict requirement for Cdk9 kinase activity for the release of paused Pol II into productive elongation. Upon inhibition of Cdk9, cells with NELF efficiently shutdown gene transcription, while in NELF-depleted cells, defective, non-productive transcription continues unabated. By introducing a strict checkpoint for Cdk9, the evolution of NELF was likely critical to enable increased regulation of Cdk9 in higher eukaryotes, as Cdk9 availability can be restricted to limit gene transcription without inducing wasteful, non-productive transcription.

Evolution of promoter-proximal pausing enabled a new layer of transcription control

Promoter-proximal pausing of RNA polymerase II (Pol II) is a key regulatory step during transcription. Despite the central role of pausing in gene regulation, we do not understand the evolutionary processes that led to the emergence of Pol II pausing or its transition to a rate-limiting step actively controlled by transcription factors. Here we analyzed transcription in species across the tree of life. We found that unicellular eukaryotes display a slow acceleration of Pol II near transcription start sites. This proto-paused-like state transitioned to a longer, focused pause in derived metazoans which coincided with the evolution of new subunits in the NELF and 7SK complexes. Depletion of NELF reverts the mammalian focal pause to a proto-pause-like state and compromises transcriptional activation for a set of heat shock genes. Collectively, this work details the evolutionary history of Pol II pausing and sheds light on how new transcriptional regulatory mechanisms evolve.

X-chromosome target specificity diverged between dosage compensation mechanisms of two closely related Caenorhabditis species

An evolutionary perspective enhances our understanding of biological mechanisms. Comparison of sex determination and X-chromosome dosage compensation mechanisms between the closely related nematode species Caenorhabditis briggsae (Cbr) and Caenorhabditis elegans (Cel) revealed that the genetic regulatory hierarchy controlling both processes is conserved, but the X-chromosome target specificity and mode of binding for the specialized condensin dosage compensation complex (DCC) controlling X expression have diverged. We identified two motifs within Cbr DCC recruitment sites that are highly enriched on X: 13 bp MEX and 30 bp MEX II. Mutating either MEX or MEX II in an endogenous recruitment site with multiple copies of one or both motifs reduced binding, but only removing all motifs eliminated binding in vivo. Hence, DCC binding to Cbr recruitment sites appears additive. In contrast, DCC binding to Cel recruitment sites is synergistic: mutating even one motif in vivo eliminated binding. Although all X-chromosome motifs share the sequence CAGGG, they have otherwise diverged so that a motif from one species cannot function in the other. Functional divergence was demonstrated in vivo and in vitro. A single nucleotide position in Cbr MEX can determine whether Cel DCC binds. This rapid divergence of DCC target specificity could have been an important factor in establishing reproductive isolation between nematode species and contrasts dramatically with the conservation of target specificity for X-chromosome dosage compensation across Drosophila species and for transcription factors controlling developmental processes such as body-plan specification from fruit flies to mice.

Topoisomerases I and II facilitate condensin DC translocation to organize and repress X chromosomes in C. elegans

Condensins are evolutionarily conserved molecular motors that translocate along DNA and form loops. To address how DNA topology affects condensin translocation, we applied auxin-inducible degradation of topoisomerases I and II and analyzed the binding and function of an interphase condensin that mediates X chromosome dosage compensation in C. elegans. TOP-2 depletion reduced long-range spreading of condensin-DC (dosage compensation) from its recruitment sites and shortened 3D DNA contacts measured by Hi-C. TOP-1 depletion did not affect long-range spreading but resulted in condensin-DC accumulation within expressed gene bodies. Both TOP-1 and TOP-2 depletion resulted in X chromosome derepression, indicating that condensin-DC translocation at both scales is required for its function. Together, the distinct effects of TOP-1 and TOP-2 suggest two distinct modes of condensin-DC association with chromatin: long-range DNA loop extrusion that requires decatenation/unknotting of DNA and short-range translocation across genes that requires resolution of transcription-induced supercoiling.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

Functional genomic assays to annotate enhancer-promoter interactions genome wide

Enhancers are pivotal for regulating gene transcription that occurs at promoters. Identification of the interacting enhancer-promoter pairs and understanding the mechanisms behind how they interact and how enhancers modulate transcription can provide fundamental insight into gene regulatory networks. Recently, advances in high-throughput methods in three major areas-chromosome conformation capture assay, such as Hi-C to study basic chromatin architecture, ectopic reporter experiments such as self-transcribing active regulatory region sequencing (STARR-seq) to quantify promoter and enhancer activity, and endogenous perturbations such as clustered regularly interspaced short palindromic repeat interference (CRISPRi) to identify enhancer-promoter compatibility-have further our knowledge about transcription. In this review, we will discuss the major method developments and key findings from these assays.

Combinatorial clustering of distinct DNA motifs directs synergistic binding of Caenorhabditis elegans dosage compensation complex to X chromosomes

Diverse regulatory mechanisms balance X-chromosome gene expression between sexes in mammals, fruit flies, and nematodes (XY/XO males and XX females/hermaphrodites). We identify DNA motifs on X that recruit dosage compensation complexes (DCCs) in nematode hermaphrodites to reduce X-chromosome expression. Recruitment sites on X, but not regions on autosomes, contain diverse combinations of different motifs or multiple copies of one motif. DCC binding studies in vivo and in vitro of wild-type and mutant X-recruitment sites validate motif usage. We find that clustering of motifs in different combinations with appropriate orientation and spacing promotes synergy in DCC binding, thereby triggering DCC assembly specifically along X. We demonstrate how regulatory complexes can be recruited across an entire chromosome to control its gene expression.

Organisms that count X-chromosome number to determine sex utilize dosage compensation mechanisms to balance X-gene expression between sexes. Typically, a regulatory complex is recruited to X chromosomes of one sex to modulate gene expression. A major challenge is to determine the mechanisms that target regulatory complexes specifically to X. Here, we identify critical X-sequence motifs in Caenorhabditis elegans that act synergistically in hermaphrodites to direct X-specific recruitment of the dosage compensation complex (DCC), a condensin complex. We find two DNA motifs that collaborate with a previously defined 12-bp motif called MEX (motif enriched on X) to mediate binding: MEX II, a 26-bp X-enriched motif and Motif C, a 9-bp motif that lacks X enrichment. Inserting both MEX and MEX II into a new location on X creates a DCC binding site equivalent to an endogenous recruitment site, but inserting only MEX or MEX II alone does not. Moreover, mutating MEX, MEX II, or Motif C in endogenous recruitment sites with multiple different motifs dramatically reduces DCC binding in vivo to nearly the same extent as mutating all motifs. Changing the orientation or spacing of motifs also reduces DCC binding. Hence, synergy in DCC binding via combinatorial clustering of motifs triggers DCC assembly specifically on X chromosomes. Using an in vitro DNA binding assay, we refine the features of motifs and flanking sequences that are critical for DCC binding. Our work reveals general principles by which regulatory complexes can be recruited across an entire chromosome to control its gene expression.

Increased gene dosage and mRNA expression from chromosomal duplications in Caenorhabditis elegans

Isolation of copy number variations and chromosomal duplications at high frequency in the laboratory suggested that Caenorhabditis elegans tolerates increased gene dosage. Here, we addressed if a general dosage compensation mechanism acts at the level of mRNA expression in C. elegans. We characterized gene dosage and mRNA expression in 3 chromosomal duplications and a fosmid integration strain using DNA-seq and mRNA-seq. Our results show that on average, increased gene dosage leads to increased mRNA expression, pointing to a lack of genome-wide dosage compensation. Different genes within the same chromosomal duplication show variable levels of mRNA increase, suggesting feedback regulation of individual genes. Somatic dosage compensation and germline repression reduce the level of mRNA increase from X chromosomal duplications. Together, our results show a lack of genome-wide dosage compensation mechanism acting at the mRNA level in C. elegans and highlight the role of epigenetic and individual gene regulation contributing to the varied consequences of increased gene dosage.

Auxin-Inducible Degron System Reveals Temporal-Spatial Roles of HSF-1 and Its Transcriptional Program in Lifespan Assurance

HSF-1 is a key regulator of cellular proteotoxic stress response and is required for animal lifespan. In C. elegans, HSF-1 mediated heat shock response (HSR) declines sharply on the first day of adulthood, and HSF-1 was proposed to function primarily during larval stages for lifespan assurance based on studies using RNAi. The tissue requirement for HSF-1 in lifespan, however, is not well understood. Using the auxin-inducible degron (AID) system, we manage to uncouple the roles of HSF-1 in development and longevity. In wild-type animals, we find HSF-1 is required during the whole self-reproductive period for lifespan. This period is extended in long-lived animals that have arrested germline stem cells (GSC) or reduced insulin/IGF-1 signaling (IIS). While depletion of HSF-1 from any major somatic tissues during development results in severe defects, HSF-1 primarily functions in the intestine and likely neural system of adults to support lifespan. Finally, by combining AID and genome-wide transcriptional analyses, we find HSF-1 directly activates the transcription of constitutively-expressed chaperone and co-chaperone genes among others in early adulthood, which underlies its roles in longevity assurance.

The X chromosome in C. elegans sex determination and dosage compensation

Abnormalities in chromosome dose can reduce organismal fitness and viability by disrupting the balance of gene expression. Unlike imbalances in chromosome dose that cause pathologies, differences in X-chromosome dose that determine sex are well tolerated. Dosage compensation mechanisms have evolved in diverse species to balance X-chromosome gene expression between sexes. Mechanisms underlying nematode X-chromosome counting to determine sex revealed how small quantitative differences in molecular signals are translated into dramatically different developmental fates. Mechanisms underlying X-chromosome dosage compensation revealed the interplay between chromatin modification and three-dimensional chromosome structure imposed by an X-specific condensin complex to regulate gene expression over vast chromosomal territories. In a surprising twist of evolution, this dosage-compensation condensin complex also regulates lifespan and tolerance to proteotoxic stress.

Dual roles for nuclear RNAi Argonautes in Caenorhabditis elegans dosage compensation

Dosage compensation involves chromosome-wide gene regulatory mechanisms which impact higher order chromatin structure and are crucial for organismal health. Using a genetic approach, we identified Argonaute genes which promote dosage compensation in Caenorhabditis elegans. Dosage compensation in C. elegans hermaphrodites is initiated by the silencing of xol-1 and subsequent activation of the dosage compensation complex which binds to both hermaphrodite X chromosomes and reduces transcriptional output by half. A hallmark phenotype of dosage compensation mutants is decondensation of the X chromosomes. We characterized this phenotype in Argonaute mutants using X chromosome paint probes and fluorescence microscopy. We found that while nuclear Argonaute mutants hrde-1 and nrde-3, as well as mutants for the piRNA Argonaute prg-1, exhibit derepression of xol-1 transcripts, they also affect X chromosome condensation in a xol-1-independent manner. We also characterized the physiological contribution of Argonaute genes to dosage compensation using genetic assays and found that hrde-1 and nrde-3 contribute to healthy dosage compensation both upstream and downstream of xol-1.

A comparison of experimental assays and analytical methods for genome-wide identification of active enhancers

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. Here, we compared 13 genome-wide RNA sequencing assays in K562 cells and showed that the nuclear run-on followed by cap-selection assay (GRO/PRO-cap) has advantages in eRNA detection and active enhancer identification. We also introduced a tool, Peak Identifier for Nascent Transcript Starts (PINTS), to identify active promoters and enhancers genome-wide and pinpoint the precise location of the 5′ transcription start sites. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected eRNA TSSs available in 120 cell and tissue types that can be accessed at https://pints.yulab.org. With the knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the way for selection and characterization of their functions in a time- and labor-efficient manner in the future.

Mechanisms of sex determination and X-chromosome dosage compensation

Abnormalities in chromosome number have the potential to disrupt the balance of gene expression and thereby decrease organismal fitness and viability. Such abnormalities occur in most solid tumors and also cause severe developmental defects and spontaneous abortions. In contrast to the imbalances in chromosome dose that cause pathologies, the difference in X-chromosome dose used to determine sexual fate across diverse species is well tolerated. Dosage compensation mechanisms have evolved in such species to balance X-chromosome gene expression between the sexes, allowing them to tolerate the difference in X-chromosome dose. This review analyzes the chromosome counting mechanism that tallies X-chromosome number to determine sex (XO male and XX hermaphrodite) in the nematode Caenorhabditis elegans and the associated dosage compensation mechanism that balances X-chromosome gene expression between the sexes. Dissecting the molecular mechanisms underlying X-chromosome counting has revealed how small quantitative differences in intracellular signals can be translated into dramatically different fates. Dissecting the process of X-chromosome dosage compensation has revealed the interplay between chromatin modification and chromosome structure in regulating gene expression over vast chromosomal territories.

Condensin DC loads and spreads from recruitment sites to create loop-anchored TADs in C. elegans

Condensins are molecular motors that compact DNA via linear translocation. In Caenorhabditis elegans, the X-chromosome harbors a specialized condensin that participates in dosage compensation (DC). Condensin DC is recruited to and spreads from a small number of recruitment elements on the X-chromosome (rex) and is required for the formation of topologically associating domains (TADs). We take advantage of autosomes that are largely devoid of condensin DC and TADs to address how rex sites and condensin DC give rise to the formation of TADs. When an autosome and X-chromosome are physically fused, despite the spreading of condensin DC into the autosome, no TAD was created. Insertion of a strong rex on the X-chromosome results in the TAD boundary formation regardless of sequence orientation. When the same rex is inserted on an autosome, despite condensin DC recruitment, there was no spreading or features of a TAD. On the other hand, when a 'super rex' composed of six rex sites or three separate rex sites are inserted on an autosome, recruitment and spreading of condensin DC led to the formation of TADs. Therefore, recruitment to and spreading from rex sites are necessary and sufficient for recapitulating loop-anchored TADs observed on the X-chromosome. Together our data suggest a model in which rex sites are both loading sites and bidirectional barriers for condensin DC, a one-sided loop-extruder with movable inactive anchor.

Insulin/IGF-1 signaling and heat stress differentially regulate HSF1 activities in germline development

Germline development is sensitive to nutrient availability and environmental perturbation. Heat shock transcription factor 1 (HSF1), a key transcription factor driving the cellular heat shock response (HSR), is also involved in gametogenesis. The precise function of HSF1 (HSF-1 in C. elegans) and its regulation in germline development are poorly understood. Using the auxin-inducible degron system in C. elegans, we uncovered a role of HSF-1 in progenitor cell proliferation and early meiosis and identified a compact but important transcriptional program of HSF-1 in germline development. Interestingly, heat stress only induces the canonical HSR in a subset of germ cells but impairs HSF-1 binding at its developmental targets. Conversely, insulin/insulin growth factor 1 (IGF-1) signaling dictates the requirement for HSF-1 in germline development and functions through repressing FOXO/DAF-16 in the soma to activate HSF-1 in germ cells. We propose that this non-cell-autonomous mechanism couples nutrient-sensing insulin/IGF-1 signaling to HSF-1 activation to support homeostasis in rapid germline growth.

RNA polymerase II speed: a key player in controlling and adapting transcriptome composition

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co‐transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read‐through transcripts. RNA Pol II speed itself is regulated in response to intra‐ and extra‐cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co‐transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

A role for condensin in mediating transcriptional adaptation to environmental stimuli

Nuclear organisation shapes gene regulation; however, the principles by which three-dimensional genome architecture influences gene transcription are incompletely understood. Condensin is a key architectural chromatin constituent, best known for its role in mitotic chromosome condensation. Yet at least a subset of condensin is bound to DNA throughout the cell cycle. Studies in various organisms have reported roles for condensin in transcriptional regulation, but no unifying mechanism has emerged. Here, we use rapid conditional condensin depletion in the budding yeast Saccharomyces cerevisiae to study its role in transcriptional regulation. We observe a large number of small gene expression changes, enriched at genes located close to condensin-binding sites, consistent with a possible local effect of condensin on gene expression. Furthermore, nascent RNA sequencing reveals that transcriptional down-regulation in response to environmental stimuli, in particular to heat shock, is subdued without condensin. Our results underscore the multitude by which an architectural chromosome constituent can affect gene regulation and suggest that condensin facilitates transcriptional reprogramming as part of adaptation to environmental changes.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Evaluation of various RNA-seq approaches for identification of gene outrons in the flatworm Opisthorchis felineus

The parasitic flatworm Opisthorchis felineus is one of the causative agents of opisthorchiasis in humans. Recently, we assembled the O. felineus genome, but the correct genome annotation by means of standard methods was hampered by the presence of spliced leader trans-splicing (SLTS). As a result of SLTS, the original 5’-end (outron) of the transcripts is replaced by a short spliced leader sequence donated from a specialized SL RNA. SLTS is involved in the RNA processing of more than half of O. felineus genes, making it hard to determine the structure of outrons and bona fide transcription start sites of the corresponding genes and operons, being based solely on mRNA-seq data. In the current study, we tested various experimental approaches for identifying the sequences of outrons in O. felineus using massive parallel sequencing. Two of them were developed by us for targeted sequencing of already processed branched outrons. One was based on sequence-specific reverse transcription from the SL intron toward the 5’-end of the Y-branched outron. The other used outron hybridization with an immobilized single-stranded DNA probe complementary to the SL intron. Additionally, two approaches to the sequencing of rRNA-depleted total RNA were used, allowing the identification of a wider range of transcripts compared to mRNAseq. One is based on the enzymatic elimination of overrepresented cDNAs, the other utilizes exonucleolytic degradation of uncapped RNA by Terminator enzyme. By using the outron-targeting methods, we were not able to obtain the enrichment of RNA preparations by processed outrons, which is most likely indicative of a rapid turnover of these trans-splicing intermediate products. Of the two rRNA depletion methods, a method based on the enzymatic normalization of cDNA (Zymo-Seq RiboFree) showed high efficiency. Compared to mRNA-seq, it provides an approximately twofold increase in the fraction of reads originating from outrons and introns. The results suggest that unprocessed nascent transcripts are the main source of outron sequences in the RNA pool of O. felineus.

Starvation Responses Throughout the Caenorhabditiselegans Life Cycle

Caenorhabditis elegans survives on ephemeral food sources in the wild, and the species has a variety of adaptive responses to starvation. These features of its life history make the worm a powerful model for studying developmental, behavioral, and metabolic starvation responses. Starvation resistance is fundamental to life in the wild, and it is relevant to aging and common diseases such as cancer and diabetes. Worms respond to acute starvation at different times in the life cycle by arresting development and altering gene expression and metabolism. They also anticipate starvation during early larval development, engaging an alternative developmental program resulting in dauer diapause. By arresting development, these responses postpone growth and reproduction until feeding resumes. A common set of signaling pathways mediates systemic regulation of development in each context but with important distinctions. Several aspects of behavior, including feeding, foraging, taxis, egg laying, sleep, and associative learning, are also affected by starvation. A variety of conserved signaling, gene regulatory, and metabolic mechanisms support adaptation to starvation. Early life starvation can have persistent effects on adults and their descendants. With its short generation time, C. elegans is an ideal model for studying maternal provisioning, transgenerational epigenetic inheritance, and developmental origins of adult health and disease in humans. This review provides a comprehensive overview of starvation responses throughout the C. elegans life cycle.

Population-scale study of eRNA transcription reveals bipartite functional enhancer architecture

Enhancer RNAs (eRNA) are unstable non-coding RNAs, transcribed bidirectionally from active regulatory sequences, whose expression levels correlate with enhancer activity. We use capped-nascent-RNA sequencing to efficiently capture bidirectional transcription initiation across several human lymphoblastoid cell lines (Yoruba population) and detect ~75,000 eRNA transcription sites with high sensitivity and specificity. The use of nascent-RNA sequencing sidesteps the confounding effect of eRNA instability. We identify quantitative trait loci (QTLs) associated with the level and directionality of eRNA expression. High-resolution analyses of these two types of QTLs reveal distinct positions of enrichment at the central transcription factor (TF) binding regions and at the flanking eRNA initiation regions, both of which are associated with mRNA expression QTLs. These two regions-the central TF-binding footprint and the eRNA initiation cores-define a bipartite architecture of enhancers, inform enhancer function, and can be used as an indicator of the significance of non-coding regulatory variants.

How affinity of the ELT-2 GATA factor binding to cis-acting regulatory sites controls Caenorhabditis elegans intestinal gene transcription

We define a quantitative relationship between the affinity with which the intestine-specific GATA factor ELT-2 binds to cis-acting regulatory motifs and the resulting transcription of asp-1, a target gene representative of genes involved in Caenorhabditis elegans intestine differentiation. By establishing an experimental system that allows unknown parameters (e.g. the influence of chromatin) to effectively cancel out, we show that levels of asp-1 transcripts increase monotonically with increasing binding affinity of ELT-2 to variant promoter TGATAA sites. The shape of the response curve reveals that the product of the unbound ELT-2 concentration in vivo [i.e. (ELT-2free) or ELT-2 'activity'] and the largest ELT-XXTGATAAXX association constant (Kmax) lies between five and ten. We suggest that this (unitless) product [Kmax×(ELT-2free) or the equivalent product for any other transcription factor] provides an important quantitative descriptor of transcription-factor/regulatory-motif interaction in development, evolution and genetic disease. A more complicated model than simple binding affinity is necessary to explain the fact that ELT-2 appears to discriminate in vivo against equal-affinity binding sites that contain AGATAA instead of TGATAA.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

EnhancerAtlas 2.0: an updated resource with enhancer annotation in 586 tissue/cell types across nine species

Enhancers are distal cis-regulatory elements that activate the transcription of their target genes. They regulate a wide range of important biological functions and processes, including embryogenesis, development, and homeostasis. As more and more large-scale technologies were developed for enhancer identification, a comprehensive database is highly desirable for enhancer annotation based on various genome-wide profiling datasets across different species. Here, we present an updated database EnhancerAtlas 2.0 (http://www.enhanceratlas.org/indexv2.php), covering 586 tissue/cell types that include a large number of normal tissues, cancer cell lines, and cells at different development stages across nine species. Overall, the database contains 13 494 603 enhancers, which were obtained from 16 055 datasets using 12 high-throughput experiment methods (e.g. H3K4me1/H3K27ac, DNase-seq/ATAC-seq, P300, POLR2A, CAGE, ChIA-PET, GRO-seq, STARR-seq and MPRA). The updated version is a huge expansion of the first version, which only contains the enhancers in human cells. In addition, we predicted enhancer–target gene relationships in human, mouse and fly. Finally, the users can search enhancers and enhancer–target gene relationships through five user-friendly, interactive modules. We believe the new annotation of enhancers in EnhancerAtlas 2.0 will facilitate users to perform useful functional analysis of enhancers in various genomes.

Serotonin signaling by maternal neurons upon stress ensures progeny survival

Germ cells are vulnerable to stress. Therefore, how organisms protect their future progeny from damage in a fluctuating environment is a fundamental question in biology. We show that in Caenorhabditis elegans, serotonin released by maternal neurons during stress ensures the viability and stress resilience of future offspring. Serotonin acts through a signal transduction pathway conserved between C. elegans and mammalian cells to enable the transcription factor HSF1 to alter chromatin in soon-to-be fertilized germ cells by recruiting the histone chaperone FACT, displacing histones, and initiating protective gene expression. Without serotonin release by maternal neurons, FACT is not recruited by HSF1 in germ cells, transcription occurs but is delayed, and progeny of stressed C. elegans mothers fail to complete development. These studies uncover a novel mechanism by which stress sensing by neurons is coupled to transcription response times of germ cells to protect future offspring.

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.

Dosage Compensation: How to Be Compensated…Or Not?

Diverse dosage compensation mechanisms have evolved across species to equalize gene expression between sexes and between the sex chromosomes and autosomes. New results show that two opposite modes of dosage compensation can occur within one species, the monarch butterfly.

Identification and dynamic quantification of regulatory elements using total RNA

The spatial and temporal regulation of transcription initiation is pivotal for controlling gene expression. Here, we introduce capped-small RNA-seq (csRNA-seq), which uses total RNA as starting material to detect transcription start sites (TSSs) of both stable and unstable RNAs at single-nucleotide resolution. csRNA-seq is highly sensitive to acute changes in transcription and identifies an order of magnitude more regulated transcripts than does RNA-seq. Interrogating tissues from species across the eukaryotic kingdoms identified unstable transcripts resembling enhancer RNAs, pri-miRNAs, antisense transcripts, and promoter upstream transcripts in multicellular animals, plants, and fungi spanning 1.6 billion years of evolution. Integration of epigenomic data from these organisms revealed that histone H3 trimethylation (H3K4me3) was largely confined to TSSs of stable transcripts, whereas H3K27ac marked nucleosomes downstream from all active TSSs, suggesting an ancient role for posttranslational histone modifications in transcription. Our findings show that total RNA is sufficient to identify transcribed regulatory elements and capture the dynamics of initiated stable and unstable transcripts at single-nucleotide resolution in eukaryotes.

CAGEfightR: analysis of 5'-end data using R/Bioconductor

BACKGROUND

5'-end sequencing assays, and Cap Analysis of Gene Expression (CAGE) in particular, have been instrumental in studying transcriptional regulation. 5'-end methods provide genome-wide maps of transcription start sites (TSSs) with base pair resolution. Because active enhancers often feature bidirectional TSSs, such data can also be used to predict enhancer candidates. The current availability of mature and comprehensive computational tools for the analysis of 5'-end data is limited, preventing efficient analysis of new and existing 5'-end data.

RESULTS

We present CAGEfightR, a framework for analysis of CAGE and other 5'-end data implemented as an R/Bioconductor-package. CAGEfightR can import data from BigWig files and allows for fast and memory efficient prediction and analysis of TSSs and enhancers. Downstream analyses include quantification, normalization, annotation with transcript and gene models, TSS shape statistics, linking TSSs to enhancers via co-expression, identification of enhancer clusters, and genome-browser style visualization. While built to analyze CAGE data, we demonstrate the utility of CAGEfightR in analyzing nascent RNA 5'-data (PRO-Cap). CAGEfightR is implemented using standard Bioconductor classes, making it easy to learn, use and combine with other Bioconductor packages, for example popular differential expression tools such as limma, DESeq2 and edgeR.

CONCLUSIONS

CAGEfightR provides a single, scalable and easy-to-use framework for comprehensive downstream analysis of 5'-end data. CAGEfightR is designed to be interoperable with other Bioconductor packages, thereby unlocking hundreds of mature transcriptomic analysis tools for 5'-end data. CAGEfightR is freely available via Bioconductor: bioconductor.org/packages/CAGEfightR .

X Chromosome Domain Architecture Regulates Caenorhabditis elegans Lifespan but Not Dosage Compensation

Mechanisms establishing higher-order chromosome structures and their roles in gene regulation are elusive. We analyzed chromosome architecture during nematode X chromosome dosage compensation, which represses transcription via a dosage-compensation condensin complex (DCC) that binds hermaphrodite Xs and establishes megabase-sized topologically associating domains (TADs). We show that DCC binding at high-occupancy sites (rex sites) defines eight TAD boundaries. Single rex deletions disrupted boundaries, and single insertions created new boundaries, demonstrating that a rex site is necessary and sufficient to define DCC-dependent boundary locations. Deleting eight rex sites (8rexΔ) recapitulated TAD structure of DCC mutants, permitting analysis when chromosome-wide domain architecture was disrupted but most DCC binding remained. 8rexΔ animals exhibited no changes in X expression and lacked dosage-compensation mutant phenotypes. Hence, TAD boundaries are neither the cause nor the consequence of DCC-mediated gene repression. Abrogating TAD structure did, however, reduce thermotolerance, accelerate aging, and shorten lifespan, implicating chromosome architecture in stress responses and aging.

The RNA Polymerase II Core Promoter in Drosophila

Transcription by RNA polymerase II initiates at the core promoter, which is sometimes referred to as the "gateway to transcription." Here, we describe the properties of the RNA polymerase II core promoter in Drosophila The core promoter is at a strategic position in the expression of genes, as it is the site of convergence of the signals that lead to transcriptional activation. Importantly, core promoters are diverse in terms of their structure and function. They are composed of various combinations of sequence motifs such as the TATA box, initiator (Inr), and downstream core promoter element (DPE). Different types of core promoters are transcribed via distinct mechanisms. Moreover, some transcriptional enhancers exhibit specificity for particular types of core promoters. These findings indicate that the core promoter is a central component of the transcriptional apparatus that regulates gene expression.

A 50 year history of technologies that drove discovery in eukaryotic transcription regulation

Transcription regulation is critical to organism development and homeostasis. Control of expression of the 20,000 genes in human cells requires many hundreds of proteins acting through sophisticated multistep mechanisms. In this Historical Perspective, I highlight the progress that has been made in elucidating eukaryotic transcriptional mechanisms through an array of disciplines and approaches, and how this concerted effort has been driven by the development of new technologies.

The Caenorhabditis elegans Transgenic Toolbox

The power of any genetic model organism is derived, in part, from the ease with which gene expression can be manipulated. The short generation time and invariant developmental lineage have made Caenorhabditis elegans very useful for understanding, e.g., developmental programs, basic cell biology, neurobiology, and aging. Over the last decade, the C. elegans transgenic toolbox has expanded considerably, with the addition of a variety of methods to control expression and modify genes with unprecedented resolution. Here, we provide a comprehensive overview of transgenic methods in C. elegans, with an emphasis on recent advances in transposon-mediated transgenesis, CRISPR/Cas9 gene editing, conditional gene and protein inactivation, and bipartite systems for temporal and spatial control of expression.

Binding of an X-Specific Condensin Correlates with a Reduction in Active Histone Modifications at Gene Regulatory Elements

Condensins are evolutionarily conserved protein complexes that are required for chromosome segregation during cell division and genome organization during interphase. In Caenorhabditis elegans, a specialized condensin, which forms the core of the dosage compensation complex (DCC), binds to and represses X chromosome transcription. Here, we analyzed DCC localization and the effect of DCC depletion on histone modifications, transcription factor binding, and gene expression using chromatin immunoprecipitation sequencing and mRNA sequencing. Across the X, the DCC accumulates at accessible gene regulatory sites in active chromatin and not heterochromatin. The DCC is required for reducing the levels of activating histone modifications, including H3K4me3 and H3K27ac, but not repressive modification H3K9me3. In X-to-autosome fusion chromosomes, DCC spreading into the autosomal sequences locally reduces gene expression, thus establishing a direct link between DCC binding and repression. Together, our results indicate that DCC-mediated transcription repression is associated with a reduction in the activity of X chromosomal gene regulatory elements.

Measuring RNA polymerase activity genome-wide with high-resolution run-on-based methods

The biogenesis of RNAs is a multi-layered and highly regulated process that involves a diverse set of players acting in an orchestrated manner throughout the transcription cycle. Transcription initiation, elongation and termination factors act on RNA polymerases to modulate their movement along the DNA template in a very precise manner, more complex than previously anticipated. Genome-scale run-on-based methodologies have been developed to study in detail the position of transcriptionally-engaged RNA polymerases. Genomic run-on (GRO), and its many variants and refinements made over the years, are helping the community to address an increasing amount of scientific questions, spanning an increasing range of organisms and systems. In this review, we aim to summarize the most relevant high throughput methodologies developed to study nascent RNA by run-on methods, compare their main features, advantages and limitations, while putting them in context with alternative ways of studying the transcriptional process.

Networks of mRNA Processing and Alternative Splicing Regulation in Health and Disease

mRNA processing events introduce an intricate layer of complexity into gene expression processes, supporting a tremendous level of diversification of the genome's coding and regulatory potential, particularly in vertebrate species. The recent development of massive parallel sequencing methods and their adaptation to the identification and quantification of different RNA species and the dynamics of mRNA metabolism and processing has generated an unprecedented view over the regulatory networks that are established at this level, which contribute to sustain developmental, tissue specific or disease specific gene expression programs. In this chapter, we provide an overview of the recent evolution of transcriptome profiling methods and the surprising insights that have emerged in recent years regarding distinct mRNA processing events - from the 5' end to the 3' end of the molecule.

Sex and death: from cell fate specification to dynamic control of X-chromosome structure and gene expression

Determining sex is a binary developmental decision that most metazoans must make. Like many organisms, Caenorhabditis elegans specifies sex (XO male or XX hermaphrodite) by tallying X-chromosome number. We dissected this precise counting mechanism to determine how tiny differences in concentrations of signals are translated into dramatically different developmental fates. Determining sex by counting chromosomes solved one problem but created another-an imbalance in X gene products. We found that nematodes compensate for the difference in X-chromosome dose between sexes by reducing transcription from both hermaphrodite X chromosomes. In a surprising feat of evolution, X-chromosome regulation is functionally related to a structural problem of all mitotic and meiotic chromosomes: achieving ordered compaction of chromosomes before segregation. We showed the dosage compensation complex is a condensin complex that imposes a specific three--dimensional architecture onto hermaphrodite X chromosomes. It also triggers enrichment of histone modification H4K20me1. We discovered the machinery and mechanism underlying H4K20me1 enrichment and demonstrated its pivotal role in regulating higher-order X-chromosome structure and gene expression.

Chromatin accessibility dynamics across C. elegans development and ageing

An essential step for understanding the transcriptional circuits that control development and physiology is the global identification and characterization of regulatory elements. Here, we present the first map of regulatory elements across the development and ageing of an animal, identifying 42,245 elements accessible in at least one Caenorhabditis elegans stage. Based on nuclear transcription profiles, we define 15,714 protein-coding promoters and 19,231 putative enhancers, and find that both types of element can drive orientation-independent transcription. Additionally, more than 1000 promoters produce transcripts antisense to protein coding genes, suggesting involvement in a widespread regulatory mechanism. We find that the accessibility of most elements changes during development and/or ageing and that patterns of accessibility change are linked to specific developmental or physiological processes. The map and characterization of regulatory elements across C. elegans life provides a platform for understanding how transcription controls development and ageing.

Determinants of promoter and enhancer transcription directionality in metazoans

Divergent transcription from promoters and enhancers is pervasive in many species, but it remains unclear if it is a general feature of all eukaryotic cis regulatory elements. To address this, here we define cis regulatory elements in C. elegans, D. melanogaster and H. sapiens and investigate the determinants of their transcription directionality. In all three species, we find that divergent transcription is initiated from two separate core promoter sequences and promoter regions display competition between histone modifications on the + 1 and −1 nucleosomes. In contrast, promoter directionality, sequence composition surrounding promoters, and positional enrichment of chromatin states, are different across species. Integrative models of H3K4me3 levels and core promoter sequence are highly predictive of promoter and enhancer directionality and support two directional classes, skewed and balanced. The relative importance of features to these models are clearly distinct for promoters and enhancers. Differences in regulatory architecture within and between metazoans are therefore abundant, arguing against a unified eukaryotic model.

Divergent transcription from promoters and enhancers occurs in many species, but it is unclear if it is a general feature of all eukaryotic cis regulatory elements. Here the authors define cis regulatory elements in worms, flies, and human; and identify several differences in regulatory architecture among metazoans.

Non-canonical Drosophila X chromosome dosage compensation and repressive topologically associated domains

Background

In animals with XY sex chromosomes, X-linked genes from a single X chromosome in males are imbalanced relative to autosomal genes. To minimize the impact of genic imbalance in male Drosophila, there is a dosage compensation complex (MSL) that equilibrates X-linked gene expression with the autosomes. There are other potential contributions to dosage compensation. Hemizygous autosomal genes located in repressive chromatin domains are often derepressed. If this homolog-dependent repression occurs on the X, which has no pairing partner, then derepression could contribute to male dosage compensation.

Results

We asked whether different chromatin states or topological associations correlate with X chromosome dosage compensation, especially in regions with little MSL occupancy. Our analyses demonstrated that male X chromosome genes that are located in repressive chromatin states are depleted of MSL occupancy; however, they show dosage compensation. The genes in these repressive regions were also less sensitive to knockdown of MSL components.

Conclusions

Our results suggest that this non-canonical dosage compensation is due to the same transacting derepression that occurs on autosomes. This mechanism would facilitate immediate compensation during the evolution of sex chromosomes from autosomes. This mechanism is similar to that of C. elegans, where enhanced recruitment of X chromosomes to the nuclear lamina dampens X chromosome expression as part of the dosage compensation response in XX individuals.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0232-y) contains supplementary material, which is available to authorized users.

Super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers

Super-enhancers and stretch enhancers represent classes of transcriptional enhancers that have been shown to control the expression of cell identity genes and carry disease- and trait-associated variants. Specifically, super-enhancers are clusters of enhancers defined based on the binding occupancy of master transcription factors, chromatin regulators, or chromatin marks, while stretch enhancers are large chromatin-defined regulatory regions of at least 3,000 base pairs. Several studies have characterized these regulatory regions in numerous cell types and tissues to decipher their functional importance. However, the differences and similarities between these regulatory regions have not been fully assessed. We integrated genomic, epigenomic, and transcriptomic data from ten human cell types to perform a comparative analysis of super and stretch enhancers with respect to their chromatin profiles, cell type-specificity, and ability to control gene expression. We found that stretch enhancers are more abundant, more distal to transcription start sites, cover twice as much the genome, and are significantly less conserved than super-enhancers. In contrast, super-enhancers are significantly more enriched for active chromatin marks and cohesin complex, and more transcriptionally active than stretch enhancers. Importantly, a vast majority of super-enhancers (85%) overlap with only a small subset of stretch enhancers (13%), which are enriched for cell type-specific biological functions, and control cell identity genes. These results suggest that super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers, and importantly, most of the stretch enhancers that are distinct from super-enhancers do not show an association with cell identity genes, are less active, and more likely to be poised enhancers.